Autocorrelation studies in two-flavour Wilson Lattice QCD using DD-HMC algorithm

Transcription

1 Autocorrelation studies in two-flavour Wilson Lattice QCD using DD-HMC algorithm Abhishek Chowdhury, Asit K. De, Sangita De Sarkar, A. Harindranath, Jyotirmoy Maiti, Santanu Mondal, Anwesa Sarkar June 25, Lattice 212, Cairns, Australia

2 Dynamical Wilson fermion simulations at smaller quark masses, smaller lattice spacings and and larger lattice volumes on currently available computers have become feasible with recent developments such as DD-HMC algorithm (Luscher). Measurements of autocorrelation times help us to evaluate the performance of an algorithm in overcoming critical slowing down. In addition, an accurate determination of the uncertainty associated with the measurement of an observable requires a realistic estimation of the autocorrelation of the observable. In this work we study autocorrelations of several observables measured with DD-HMC algorithm using naive Wilson fermion and gauge action.

3 Definitions The normalized autocorrelation function for an observable O is defined as, Γ O (t) = C O (t)/c O () where C O (t) = 1 N t N t r=1 The integrated autocorrelation time is, (O r O )(O r+t O ), O = 1 N τint O = 1 window 2 + Γ O (t) t=1 N O t. t=1

4 The error, δ O = 2τint (O)var[O] For any stationary Markov chain satisfiying ergodicity and detailed balance, N Γ O (t) = e t/τ n η n (O) 2 n 1 where τ n = 1 lnλ n. λ n s are real eigenvalues of symmetrized probability transition matrix with λ = 1 and λ n < 1 for n 1. τ 1 = 1 lnλ 1 exponential autocorrelation time (τ exp ). Different obsevables couples differently to the eigenmodes. Γ O (t) cannot be negative.

5 Simulation Details β = 5.6 tag lattice κ block N 2 N trj τ B 1b B 3a,, B 3b,, B 4a,, B 4b,, B 5a,, B 5b,, C C 3,, C 4,, β = 5.8 tag lattice κ block N 2 N trj τ D D 5,, Table: Here block, N 2, N trj, τ refers to HMC block, step number for the force F 2, number of HMC trajectories and the Molecular Dynamics trajectory length respectively.

6 Negativity of autocorrelation function for plaquette B 5b : L24T48, β = 5.6, κ = N trj = 5 B 5b : L24T48, β = 5.6, κ = N trj = Γ (t) Γ (t) t t B 5b : L24T48, β = 5.6, κ = N trj = Γ (t) t

7 S. Schaefer et al. [ALPHA Collaboration], Nucl. Phys. B 845, 93 (211); S. Schaefer and F. Virotta, PoS LATTICE 21, 42 (21). Improved Estimation of τ int : Let τ be the best estimate of dominant time constant. If for an observable O all relevant time scales are smaller or of the same order of τ then the upper bound of τ int, τint u = 1 W u 2 + Γ O (W u ) + A O (W u )τ where A O = max(γ O (W u ),2δΓ O (W u )), t=1 W u chosen where autocorrelation function is still significant. Estimate of τ : τ exp eff τ eff = = Max O t 2ln ΓO (t/2) Γ O (t) t 2ln ΓO (t/2) Γ O (t)

8 Γ(t) τ eff t Γ(t) τ eff t Figure: Normalized autocorrelation function and effective autocorrelation time for P (left) Q 2 2 (right) for the ensemble B 3b.

9 Autocorrelations for topological susceptibility (Q 2 2 ) 3 L24T48, β = 5.6, κ =.158 L24T48, β = 5.6, κ = B 3b 15 D 1 : L32T64, β = 5.8, κ =.1543 D 5 : L32T64, β = 5.8, κ = τ int τ int D 1 D B 4b W/ W/32 Figure: Integrated autocorrelation times for topological suceptibilities (Q2 2 ) at β = 5.6 (a =.69 fm) (left) and at β = 5.8 (right) (a =.53 fm). τ int (Q 2 2 ) as κ

10 4 3 L24T48, β = 5.6, κ =.158 L24T48, β = 5.6, κ = B 3b 2 15 D 1 : L32T64, β = 5.8, κ =.1543 D 5 : L32T64, β = 5.8, κ = D 1 τ int 2 τ int B 4b W/6 D W/ Figure: Integrated autocorrelation times and their upper bounds (τint u ) for topological suceptibilities (Q2 2 ) at β = 5.6 (a =.7 fm) (left) and at β = 5.8 (right) (a =.55 fm). τ int (Q 2 2 ) & τu int (Q2 2 ) as κ

11 Topological Charge Density Correlator (C(r)) 1.8 Tag : B : β = 5.6, Tag : D : β = 5.8 D 1 1 Tag B : β = 5.6, Tag D : β = 5.8 Γ 2 Q2 (t).6.4 B 3b t/ Γ(t).5 D 1, Q 2 2 D 1, C(r=12.) B 3b, Q 2 2 B 3b, C(r=12.) t /96 Figure: (left) Comaparison of normalized autocorrelation functions for Q2 2 at β = 5.6 and 5.8. (right) Comparison of β dependence of normalized autocorrelation functions for Q2 2 and C(r = 12.). C(r) is affected only slightly by critical slowing down. More about the properties of C(r) talk by SM, Session: Chiral Symmetry

12 Effect of Size and Smearing τ int D 1 : L32T64, β = 5.8, κ =.1543 R=1, T=1 R=3, T=3 R=5, T=4 τ int C 2 : L32T64, β = 5.6, κ =.158 HYP 5 HYP 15 HYP 25 HYP 35 HYP W W Figure: Integrated autocorrelation times of Wilson loops for different sizes (left) and for different levels of HYP smearing (right). Autocorrelation increases with increasing size and smearing level.

13 Pion and Nucleon Propagators β = 5.6 tag κ τint Pion τ Nucleon int B 3a (19) 75(18) B 4a (9) 34(9) B 5a (1) 25(9) C (13) 33(17) C (15) 26(7) C (11) 18(6) Table: Integrated autocorrelation times for pion (PP) and nucleon propagators with wall sources. Autocorrelation decreases with increasing κ. About low lying hadron hadrons and chiral condensate talk by Asit De, session: Chiral Symmetry

14 6 4 PP 6 4 AP 2 2 τ int PA AA Window Figure: Integrated autocorrelation times for PP, AP, PA and AA correlators with wall source for the ensemble B 3a. Measurements are done with a gap of 24 trajectories. P = qγ 5 q and A = qγ 4 γ 5 q q = u/d quark.

15 Conclusions Autocorrelations of topological susceptibility, pion and nucleon propagators decrease with decreasing quark mass. The topological charge density correlator is affected only slightly by critical slowing down compared to topological susceptibility. Increasing the size and the smearing level increases the autocorrelation of Wilson loop.

16 β = 5.6 MeV κ m q m pi β = 5.8, Volume= MeV κ m q m pi

17 The statistical variance of the measured value O is, σ 2 = 2τ int σ 2 For any stationary Markov chain Γ O (t) = (λ n ) t η n (O) 2. n 1 λ n, are the eigenvalues of the matrix, 2 T x,y = π 1 x P xy π 2 1 y P xy probability transition matrix of the Markov chain, π x is the stationary distribution. η n (O) = x O(x)χ n (x)π 1 2 x corresponding to λ n. where χ n (x) is the eigenfunction

18 T x,y is positive definite. Now if the algorithm satisfies detailed balance T x,y is symmetric. Then by Perron-Frobenius theorem, T x,y has real positive eigenvalues λ n, n with λ = 1 and λ n < 1 for n 1. Hence, Γ O (t) = e t/τ n η n (O) 2 n 1 where τ n = 1 lnλ n. τ 1 = 1 lnλ 1 exponential autocorrelation time (τ exp ).

19 Q 2 2 β κ τ int τ u int (27) 276(3) (93) 156(19) C(r = 12.) β κ τ int τ u int (53) 314(76) (83) 85(168) Table: Integrated autocorrelation times (τ int ) and their upper bounds (τint u ) for Q2 2 and C(r) in two β s.