846 J. Bulava et al. / Nuclear Physics B 910 (2016) 842–867
evident in the measured correlation functions.
2
In order to mitigate effects of this autocorrelation
on estimates of statistical uncertainties, we average measurements on pairs of subsequent con-
figurations. Statistical errors are estimated using the bootstrap technique [26] on this rebinned
ensemble with N
B
=800 bootstrap samples.
In this anisotropic setup the lattice re
gulator is fully specified by the temporal lattice spacing
a
t
and renormalized anisotropy ξ = a
s
/a
t
. The anisotropy is determined from the (continuum)
pion dispersion relation
a
t
E
π
(d
2
)
2
=(a
t
m
π
)
2
+
2πa
s
ξL
2
d
2
, (2.1)
where d ∈Z
3
is the quantized finite-volume three-momentum of the pion.
Determination of ξ requires the single-pion ener
gies a
t
E
π
(d
2
) in Eq. (2.1). Periodic tempo-
ral boundary conditions are used for this ensemble, potentially complicating the extraction of
finite-volume energies from the fall-off of temporal correlation functions. In particular, a sin-
gle zero-momentum pion correlation function has the ‘cosh’ form in the limit of ground-state
saturation
lim
t1/E
1
,
T −t 1/E
1
C
π
(t) = Ae
−m
π
t
1 + e
−m
π
(T −2t)
, (2.2)
where E
1
is the relevant first excited-state energy. Two-pion correlation functions with zero total
momentum have the more complicated form (ignoring small energy shifts due to pion interac-
tions)
lim
t1/E
1
,
T −t 1/E
1
C
2π
(t) = Ae
−2m
π
t
1 + e
−2m
π
(T −2t)
+Be
−m
π
(T −2t)
, (2.3)
while two-pion correlation functions with non-zero total momenta have a similar but more com-
plicated additional exponential term.
Since our finite-v
olume energies are extracted from fits of temporal correlation functions to
an exponential form, these additional terms add potentially significant complication as has been
discussed in (e.g.) Ref. [9]. Fortunately, the large temporal extent of our lattice (m
π
T ≈ 10) sup-
presses such terms below the statistical accuracy of the energy levels. This can be demonstrated
by performing two-parameter correlated-χ
2
fits of the single zero-momentum pion correlation
function to both a single exponential and the cosh of Eq. (2.2). The second exponential in
Eq. (2.2) is larger than or equal to the additional problematic exponential terms which appear
in two-pion correlation functions, apart from small hadronic interaction effects.
As the single pion at rest is our most precisely determined correlation function, it is most
sensiti
ve to these thermal effects. The absence of these effects, such as the second exponential
in Eq. (2.2), indicates that additional exponentials in two-pion correlation functions may be ne-
glected. The comparison of single-exponential and cosh fits for v
arious fitting ranges is shown
in Fig. 1. Clearly, no effect from the finite temporal extent is evident for these temporal separa-
tions. All subsequent correlated-χ
2
fits to temporal correlation functions will thus ignore finite-T
effects. The extraction of these energies will be discussed in more detail in Sec. 2.3.
The fits to the single-pion zero-momentum correlation function sho
wn in Fig. 1, as well as all
other fits to correlated data in this work, minimize a correlated-χ
2
to properly treat the covariance
2
In lattice QCD the largest autocorrelations are typically observed for ‘smoothed’ observables such as the topological
charge and smoothed action [25], which are not examined here.