BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Stephen R. Sharpe appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
I review recent progress in calculating scattering amplitudes and resonance properties involving three particles using results from lattice QCD. The necessary input is the finite-volume spectrum, and the outputs -- via solutions of integral equations -- are scattering amplitudes that can be continued into the complex plane to search for resonance poles. I describe the outlook for future extensions and applications of this work.
We study pion-pion scattering in partially-quenched twisted-mass lattice QCD using chiral perturbation theory. The specific partially-quenched setup corresponds to that used in numerical lattice QCD calculations of the $I=0$ scattering length. We study the discretization errors proportional to $a^2$, with $a$ the lattice spacing, and the errors that arise due to the use of Lüscher's two-particle quantization condition in a theory that is not unitary. We argue that the former can be as large as $\sim 100\%$, but explain how they can be systematically subtracted using a calculation of the $I=2$ scattering amplitude in the same partially-quenched framework. We estimate the error from the violation of unitarity to be $\sim 25\%$, and argue that this error will be difficult to reduce in practice.
The status and prospects for investigations of exotic and conventional hadrons with lattice QCD are discussed. The majority of hadrons decay strongly via one or multiple decay-channels, including most of the experimentally discovered exotic hadrons. Despite this difficult challenge, the properties of several hadronic resonances have been determined within lattice QCD. To further discern the spectroscopic properties of various hadrons and to help resolve their nature we present our suggestions for future analytic and lattice studies.
Recently, the formalism needed to relate the finite-volume spectrum of systems of nondegenerate spinless particles has been derived. In this work we discuss a range of issues that arise when implementing this formalism in practice, provide further theoretical results that can be used to check the implementation, and make available codes for implementing the three-particle quantization condition. Specifically, we discuss the need to modify the upper limit of the cutoff function due to the fact that the left-hand cut in the scattering amplitudes for two nondegenerate particles moves closer to threshold; we describe the decomposition of the three-particle amplitude $\mathcal{K}_\text{df,3}$ into the matrix basis used in the quantization condition, including both $s$ and $p$ waves, with the latter arising in the amplitude for two nondegenerate particles; we derive the threshold expansion for the lightest three-particle state in the rest frame up to $\mathcal{O}(1/L^5)$; and we calculate the leading-order predictions in chiral perturbation theory for $\mathcal{K}_\text{df,3}$ in the $π^+π^+K^+$ and $π^+K^+K^+$ systems. We focus mainly on systems with two identical particles plus a third that is different ("2+1" systems). We describe the formalism in full detail, and present numerical explorations in toy models, in particular checking that the results agree with the threshold expansion, and making a prediction for the spectrum of $π^+π^+K^+$ levels using the two- and three-particle interactions predicted by chiral perturbation theory.
We present a generalization of the relativistic, finite-volume, three-particle quantization condition for non-identical pions in isosymmetric QCD. The resulting formalism allows one to use discrete finite-volume energies, determined using lattice QCD, to constrain scattering amplitudes for all possible values of two- and three-pion isospin. As for the case of identical pions considered previously, the result splits into two steps: The first defines a non-perturbative function with roots equal to the allowed energies, $E_n(L)$, in a given cubic volume with side-length $L$. This function depends on an intermediate three-body quantity, denoted $\mathcal{K}_{\mathrm{df},3}$, which can thus be constrained from lattice QCD input. The second step is a set of integral equations relating $\mathcal{K}_{\mathrm{df},3}$ to the physical scattering amplitude, $\mathcal M_3$. Both of the key relations, $E_n(L) \leftrightarrow \mathcal{K}_{\mathrm{df},3}$ and $\mathcal{K}_{\mathrm{df},3}\leftrightarrow \mathcal M_3$, are shown to be block-diagonal in the basis of definite three-pion isospin, $I_{πππ}$, so that one in fact recovers four independent relations, corresponding to $I_{πππ}=0,1,2,3$. We also provide the generalized threshold expansion of $\mathcal{K}_{\mathrm{df},3}$ for all channels, as well as parameterizations for all three-pion resonances present for $I_{πππ}=0$ and $I_{πππ}=1$. As an example of the utility of the generalized formalism, we present a toy implementation of the quantization condition for $I_{πππ}=0$, focusing on the quantum numbers of the $ω$ and $h_1$ resonances.
In this article we complete our formalism relating the finite-volume energy spectrum of a scalar quantum field theory to the three-to-three scattering amplitude, ${\cal M}_3$. In previous work we found a quantization condition relating the spectrum to a non-standard infinite-volume quantity, denoted ${\cal K}_{{\rm df},3}$. Here we present the relation between ${\cal K}_{{\rm df},3}$ and ${\cal M}_3$. We then discuss briefly how our now completed formalism can be practically implemented to extract ${\cal M}_3$ from the finite-volume energy spectrum.
We report progress on extending the relativistic model-independent quantization condition for three particles, derived previously by two of us, to a broader class of theories, as well as progress on checking the formalism. In particular, we discuss the extension to include the possibility of 2->3 and 3->2 transitions and the calculation of the finite-volume energy shift of an Efimov-like three-particle bound state. The latter agrees with the results obtained previously using non-relativistic quantum mechanics.
In a recent paper we used chiral perturbation theory to determine the phase diagram and pion spectrum for Wilson and twisted-mass fermions at non-zero lattice spacing with non-degenerate up and down quarks. Here we extend this work to include the effects of electromagnetism, so that it is applicable to recent simulations incorporating all sources of isospin breaking. For Wilson fermions, we find that the phase diagram is unaffected by the inclusion of electromagnetism---the only effect is to raise the charged pion masses. For maximally twisted fermions, we previously took the twist and isospin-breaking directions to be different, in order that the fermion determinant is real and positive. However, this is incompatible with electromagnetic gauge invariance, and so here we take the twist to be in the isospin-breaking direction, following the RM123 collaboration. We map out the phase diagram in this case, which has not previously been studied. The results differ from those obtained with different twist and isospin directions. One practical issue when including electromagnetism is that the critical masses for up and down quarks differ. We show that one of the criteria suggested to determine these critical masses does not work, and propose an alternative.
We present results for the matrix elements of the additional $ΔS=2$ operators that appear in models of physics beyond the Standard Model (BSM), expressed in terms of four BSM $B$-parameters. Combined with experimental results for $ΔM_K$ and $ε_K$, these constrain the parameters of BSM models. We use improved staggered fermions, with valence HYP-smeared quarks and $N_f=2+1$ flavors of "asqtad" sea quarks. The configurations have been generated by the MILC collaboration. The matching between lattice and continuum four-fermion operators and bilinears is done perturbatively at one-loop order. We use three lattice spacings for the continuum extrapolation: $a\approx 0.09$, $0.06$ and $0.045\;$fm. Valence light-quark masses range down to $\approx m_s^{\rm phys}/13$ while the light sea-quark masses range down to $\approx m_s^{\rm phys}/20$. Compared to our previous published work, we have added four additional lattice ensembles, leading to better controlled extrapolations in the lattice spacing and sea-quark masses. We report final results for two renormalization scales, $μ=2\;\text{GeV}$ and $3\;\text{GeV}$, and compare them to those obtained by other collaborations. Agreement is found for two of the four BSM $B$-parameters ($B_2$ and $B_3^\text{SUSY}$). The other two ($B_4$ and $B_5$) differ significantly from those obtained using RI-MOM renormalization as an intermediate scheme, but are in agreement with recent preliminary results obtained by the RBC-UKQCD collaboration using RI-SMOM intermediate schemes.
We calculate the energy of the state closest to threshold for two and three identical, spinless particles confined to a cubic spatial volume with periodic boundary conditions and with zero total momentum in the finite-volume frame. The calculation is performed in relativistic quantum field theory with particles coupled via a $λϕ^4$ interaction, and we work through order $λ^3$. The energy shifts begin at ${\cal O}(1/L^3)$, and we keep subleading terms proportional to $1/L^4$, $1/L^5$ and $1/L^6$. These terms allow a non-trivial check of the results obtained from quantization conditions that hold for arbitrary interactions, namely that of Lüscher for two particles and our recently developed formalism for three particles. We also compare to previously obtained results based on non-relativistic quantum mechanics.
We determine the phase diagram and chiral condensate for lattice QCD with two flavors of twisted-mass fermions in the presence of nondegenerate up and down quarks, discretization errors and a nonzero value of $Θ_{\rm QCD}$. Although such a theory has a complex action and cannot, at present, be simulated, the results are needed to understand how to tune to maximal twist in the presence of electromagnetism, a topic discussed in a companion paper. We find that, in general, the only phase structure is a first-order transition of finite length. Pion masses are nonvanishing throughout the phase plane except at the endpoints of the first-order line. Only for extremal values of the twist angle and $Θ_{\rm QCD}$ ($ω=0$ or $π/2$ and $Θ_{\rm QCD}=0$ or $π$) are there second-order transitions.
We present a generalization of Luescher's relation between the finite-volume spectrum and scattering amplitudes to the case of three particles. We consider a relativistic scalar field theory in which the couplings are arbitrary aside from a Z2 symmetry that removes vertices with an odd number of particles. The theory is assumed to have two-particle phase shifts that are bounded by π/2 in the regime of elastic scattering. We determine the spectrum of the finite-volume theory from the poles in the odd-particle-number finite-volume correlator, which we analyze to all orders in perturbation theory. We show that it depends on the infinite-volume two-to-two K-matrix as well as a nonstandard infinite-volume three-to-three K-matrix. A key feature of our result is the need to subtract physical singularities in the three-to-three amplitude and thus deal with a divergence-free quantity. This allows our initial, formal result to be truncated to a finite dimensional determinant equation. At present, the relation of the three-to-three K-matrix to the corresponding scattering amplitude is not known, although previous results in the non-relativistic limit suggest that such a relation exists.
We present updated results for kaon B-parameters for operators arising in models of new physics. We use HYP-smeared staggered quarks on the $N_f = 2+1$ MILC asqtad lattices. During the last year we have added new ensembles, which has necessitated chiral-continuum fitting with more elaborate fitting functions. We have also corrected an error in a two-loop anomalous dimension used to evolve results between different scales. Our results for the beyond-the-Standard-Model B-parameters have total errors of $5-10$\%. We find that the discrepancy observed last year between our results and those of the RBC/UKQCD and ETM collaborations for some of the B-parameters has been reduced from $4\!-\!5\,σ$ to $2\!-\!3\,σ$.
We present results for the kaon mixing parameter $B_K$ obtained using improved staggered fermions on a much enlarged set of MILC asqtad lattices. Compared to our previous publication, which was based largely on a single ensemble at each of the three lattice spacings $a\approx 0.09\;$fm, $0.06\;$fm and $0.045\;$fm, we have added seven new fine and four new superfine ensembles, with a range of values of the light and strange sea-quark masses. We have also increased the number of measurements on one of the original ensembles. This allows us to do controlled extrapolations in the light and strange sea-quark masses, which we do simultaneously with the continuum extrapolation. This reduces the extrapolation error and improves the reliability of our error estimates. Our final result is $\hat{B}_K = 0.7379 \pm 0.0047 (\text{stat}) \pm 0.0365 (\text{sys})$.
We determine the phase diagram and pion spectrum for Wilson and twisted-mass fermions in the presence of non-degeneracy between the up and down quark and discretization errors, using Wilson and twisted-mass chiral perturbation theory. We find that the CP-violating phase of the continuum theory (which occurs for sufficiently large non-degeneracy) is continuously connected to the Aoki-phase of the lattice theory with degenerate quarks. We show that discretization effects can, in some cases, push simulations with physical masses closer to either the CP-violating phase or another phase not present in the continuum, so that at sufficiently large lattice spacings physical-point simulations could lie in one of these phases.
In this talk we determine the phase diagram and pion spectrum for Wilson and twisted-mass fermions in the presence of non-degeneracy between the up and down quark and discretization errors. We find that the CP-violating phase of the continuum theory, which occurs for sufficiently large non-degeneracy, is continuously connected to the Aoki phase found in the lattice theory with degenerate quarks. Both for the Aoki and first-order scenarios, this results in a critical surface along which at least one of the pions is massless. In the pion spectrum, we focus mainly on the untwisted case, where there is competition between the effects of non-degeneracy and discretization errors. A more extensive analysis can be found in our recent paper [11].
We give an update on our derivation of a quantization condition relating the finite-volume spectrum of three particles in a cubic box to infinite-volume scattering quantities. We have discovered and fixed technical problems in the derivation sketched in the proceedings of last year's lattice conference [arXiv:1311.4848], and have presented a detailed description of the corrected derivation in Ref. [arXiv:1408.5933]. Here we give an overview of the problems and their solutions, and describe open questions.
A recent numerical lattice calculation of the kaon mixing matrix elements of general $ΔS=2$ four-fermion operators using staggered fermions relied on two auxiliary theoretical calculations. Here we describe the methodology and present the results of these two calculations. The first concerns one-loop matching coefficients between staggered lattice operators and the corresponding continuum operators. Previous calculations with staggered fermions have used a non-standard regularization scheme for the continuum operators, and here we provide the additional matching factors needed to connect to the standard regularization scheme. This is the scheme in which two-loop anomalous dimensions are known. We also observe that all previous calculations of this operator matching using staggered fermions have overlooked one matching step in the continuum. This extra step turns out to have no impact on three of the five operators (including that relevant for $B_K$), but does affect the other two operators. The second auxiliary calculation concerns the two-loop renormalization group (RG) evolution equations for the $B$-parameters of the $ΔS=2$ operators. For one pair of operators, the standard analytic solution to the two-loop RG equations fails due to a spurious singularity introduced by the approximations made in the calculation. We give a non-singular expression derived using analytic continuation, and check the result using a numerical solution to the RG equations. We also describe the RG evolution for "golden" combinations of $B$-parameters, and give numerical results for RG evolution matrices needed in the companion lattice calculation.
We calculate the kaon mixing B-parameters for operators arising generically in theories of physics beyond the standard model. We use HYP-smeared improved staggered fermions on the $N_f = 2+1$ MILC asqtad lattices. Operator matching is done perturbatively at one-loop order. Chiral extrapolations are done using "golden combinations" in which one-loop chiral logarithms are absent. For the combined sea-quark mass and continuum extrapolation, we use three lattice spacings: $a \approx 0.045, 0.06$ and $0.09 \text{fm}$. Our results have a total error of 5-6%, which is dominated by the systematic error from matching and continuum extrapolation. For two of the BSM $B$-parameters, we agree with results obtained using domain-wall and twisted-mass dynamical fermions, but we disagree by $(4-5)σ$ for the other two.
We present results for matrix elements of $ΔS=2$ four-fermion operators arising generically in models of new physics. These are needed to constrain such models using the measured values of $\varepsilon_K$ and $ΔM_K$. We use lattice QCD with 2+1 flavors of improved staggered fermions on lattices generated by the MILC collaboration. We extrapolate to the continuum from three lattice spacings ranging down to $a\approx 0.045\;$fm. Total errors are $\sim 5-6%$, arising primarily from our use of one-loop matching between lattice and continuum operators. For two of the matrix elements, our results disagree significantly from those obtained using different fermion discretizations.
We apply non-perturbative renormalization to bilinears composed of improved staggered fermions. We explain how to generalize the method to staggered fermions in a way which is consistent with the lattice symmetries, and introduce a new type of lattice bilinear which transforms covariantly and avoids mixing. We derive the consequences of lattice symmetries for the propagator and vertices. We implement the method numerically for hypercubic-smeared (HYP) and asqtad valence fermion actions, using lattices with asqtad sea quarks generated by the MILC collaboration. We compare the non-perturbative results so obtained to those from perturbation theory, using both scale-independent ratios of bilinears (of which we calculate 26), and the scale-dependent bilinears themselves. Overall, we find that one-loop perturbation theory provides a successful description of the results for HYP-fermions if we allow for a truncation error of roughly the size of the square of the one-loop term (for ratios) or of size O(1) \times α^2 (for the bilinears themselves). Perturbation theory is, however, less successful at describing the non-perturbative asqtad results.
This is a combined write-up for two talks which were given consecutively and which described different aspects of the same topic. We present a generalization of Lüscher's relation between the finite-volume spectrum and S-matrix to three particles. Specifically, we consider a scalar field theory, which has a $\mathbb{Z}_2$ symmetry that prevents even/odd coupling. The theory is assumed to have no two-particle bound states and to have two-particle phase shifts that are bounded by $π/2$ in the regime of elastic scattering. Considering center of mass energies between one and five particle masses, we evaluate a three-to-three finite-volume correlator to all orders in perturbation theory. Only terms which are exponentially suppressed in volume are neglected. From poles in the correlator we then determine the relation between finite-volume spectrum and scattering quantities. In this analysis one must carefully treat the unitary cusp at two-particle threshold. This point, which was neglected in the conference talks, is described in some detail here. We also describe an important check on our main result by reproducing the large volume expansion of the energy shift from the three-particle threshold. This is found to be consistent with previous work through four non-trivial orders.
We calculate the B-parameters for operators arising in theories of new physics beyond the standard model (BSM) using HYP-smeared improved staggered fermions on the MILC asqtad lattices with N_f = 2+1 flavors. We use three different lattice spacings (a ~ 0.045, 0.06 and 0.09 fm) at a fixed ratio of light to strange quarks, m_l/m_s = 1/5, to obtain the continuum results. Operator matching is done using perturbative matching at one-loop order, and results are run to 2 or 3 GeV using two-loop running in the MS-bar scheme. For the chiral and continuum extrapolations, we use SU(2) staggered chiral perturbation theory. We present preliminary results with only statistical errors.
Models of new physics induce K_0-K_0bar mixing through operators having Dirac structures other than the "left-left" form of the Standard Model. To carry out the chiral-continuum extrapolation of results from numerical simulations, one needs to know the quark mass and lattice spacing dependence of the corresponding B-parameters in the partially quenched theory at least at next-to-leading order. For simulations using staggered fermions (such as that we are doing with HYP-smeared valence fermions on the MILC asqtad lattices) one must determine this dependence using staggered chiral perturbation theory (SChPT). We have calculated the required dependence in both SU(3) and SU(2) SChPT, working at next-to-leading order, and we give here an overview of the methodology and results. The SU(3) SChPT result turns out to be much simpler than that for the Standard Model B_K operator, due to the absence of chiral suppression for the new operators. The SU(2) SChPT result turns out to be closely related to that for B_K: the chiral logarithms are identical, up to an operator-dependent sign. Our results are also useful for fermions with chiral symmetry as they provide, in the continuum limit, the partially quenched generalization of existing continuum results.
We present a calculation of the kaon $B$-parameter, $B_K$, using lattice QCD. We use improved staggered valence and sea fermions, the latter generated by the MILC collaboration with $N_f=2+1$ light flavors. To control discretization errors, we use four different lattice spacings ranging down to $a\approx 0.045\;$fm. The chiral and continuum extrapolations are done using SU(2) staggered chiral perturbation theory. Our final result is $\hat{B}_K = 0.727 \pm 0.004 (\text{stat}) \pm 0.038 (\text{sys})$, where the dominant systematic error is from our use of truncated (one-loop) matching factors.
Models of new physics induce K-Kbar mixing operators having Dirac structures other than the "left-left" form of the Standard Model. We calculate the functional form of the corresponding B-parameters at next-to-leading order in both SU(3) and SU(2) staggered chiral perturbation theory (SChPT). Numerical results for these matrix elements are being generated using improved staggered fermions; our results can be used to extrapolate these matrix elements to the physical light and strange quark masses. The SU(3) SChPT results turn out to be much simpler than that for the Standard Model B_K operator, due to the absence of chiral suppression in the new operators. The SU(2) SChPT result is of similar simplicity to that for B_K. In fact, in the latter case, the chiral logarithms for two of the new B-parameters are identical to those for B_K, while those for the other two new B-parameters are of opposite sign. In addition to providing results for the 2+1 flavor theory in SU(3) SChPT and the 1+1+1 flavor theory in SU(2) SChPT, we present the corresponding continuum partially quenched results, as these are not available in the literature.
We use lattice simulations to study the single-site version of SU(N) lattice gauge theory with two flavors of Wilson-Dirac fermions in the adjoint representation, a theory whose large volume correspondent is expected to be conformal or nearly conformal. Working with N as large as 53, we map out the phase diagram in the plane of bare `t Hooft coupling, g^2 N, and of the lattice quark mass, a*m, and look for the region where the Z_N^4 center symmetry of the theory is intact. In this region one expects the large-N equivalence of the single site and infinite volume theories to be valid. As for the N_f=1 case (see Phys. Rev. D80: 065031), we find that the center-symmetric region is large and includes both light fermion masses and masses at the cutoff scale. We study the N-dependence of the width of this region and find strong evidence that it remains of finite width as N goes to infinity. Simulating with couplings as small as g^2 N = 0.005, we find that the width shrinks slowly with decreasing g^2 N, at a rate consistent with analytic arguments. Within the center-symmetric region our results for the phase structure, when extrapolated to infinite N, apply also for the large volume theory, which is minimal walking technicolor with N=infinity. We find a first-order transition as a function of a*m for all values of b, which we argue favors that the theory is confining in the infrared. Finally, we measure the eigenvalue densities of the Wilson-Dirac operator and its hermitian version, and use large Wilson loops to study the utility of reduction for extracting physical observables.
We analyze the large-N behavior of SU(N) lattice gauge theories with adjoint fermions by studying volume-reduced models, as pioneered by Eguchi and Kawai. We perform simulations on a single-site lattice for Nf = 1 and Nf = 2 Wilson Dirac fermions with values of N up to 53. We show for both values of Nf that in the large-N limit there is a finite region, containing both light and heavy fermions, of unbroken center symmetry where the theory exhibits volume independence. Using large-N reduction we attempt to calculate physical quantities such as the string tension and meson masses.
We generalize the Lellouch-Luscher formula, relating weak matrix elements in finite and infinite volumes, to the case of multiple strongly-coupled decay channels into two scalar particles. This is a necessary first step on the way to a lattice QCD calculation of weak decay rates for processes such as D -> pi pi and D -> KK. We also present a field theoretic derivation of the generalization of Luscher's finite volume quantization condition to multiple two-particle channels. We give fully explicit results for the case of two channels, including a form of the generalized Lellouch-Luscher formula expressed in terms of derivatives of the energies of finite volume states with respect to the box size. Our results hold for arbitrary total momentum and for degenerate or non-degenerate particles.
We describe a generalization of the Lellouch-Lüscher formula to the case of multiple strongly-coupled decay channels. As in the original formula, our final result is a relation between weak matrix elements in finite and infinite volumes. Our extension is limited to final states with two scalar particles, with center of mass energies below the lowest three- or four-particle threshold. Otherwise the extension is general, accommodating any number of channels, arbitrary strong coupling between channels, as well as any form of weak decay operators in the matrix elements. Among many possible applications, we emphasize that this is a necessary first step on the way to a lattice-QCD calculation of weak decay rates for D -> pi pi and D -> K K-bar. Our results allow for arbitrary total momentum and hold for degenerate or non-degenerate particles.
We compute Z-factors for general staggered bilinears on fine (a \approx 0.09 fm) MILC ensembles using both asqtad and HYP-smeared valence actions, comparing the results to the predictions of one-loop perturbation theory. This is an extension of previous work on the coarse (a \approx 0.12 fm) MILC ensembles. It provides a laboratory for studying NPR methodology in the staggered context, and is an important stepping stone for fully non-perturbative matching factors in ongoing computations of B_K and other weak matrix elements. We also implement non-exceptional RI/SMOM renormalization conditions using the asqtad action and present first results.
We report on recent progress in the calculation of $B_K$ using HYP-smeared improved staggered fermions on the MILC asqtad lattices. We have added measurements on fine ($a\sim 0.09 $fm) and superfine ($a\sim 0.06 $fm) ensembles at different values of the light sea quark mass ($a m_\ell$), as well as increased the statistics on some other ensembles. We find that the results on the fine lattices show a significantly stronger $am_\ell$ dependence than those on the superfine and coarse ($a\sim 0.12 $fm) lattices. We discuss different methods for accounting for these new results when doing the $a m_\ell$ and continuum extrapolations.
We present recent progress in our calculation of $B_K$ with improved staggered fermions using chiral extrapolations based on SU(3) staggered chiral perturbation theory. We have accumulated significantly higher statistics on the coarse, fine, and ultrafine MILC asqtad lattices. This leads to a reduction in statistical error and an improved continuum extrapolation. Our updated result is $\hat{B}_K = B_K(\text{RGI}) = 0.737 \pm 0.003(\text{stat}) \pm 0.046 (\text{sys})$. This is consistent with the result obtained using chiral extrapolations based on SU(2) staggered chiral perturbation theory, although the total error is somewhat larger with the SU(3) analysis.
Wilson chiral perturbation theory (WChPT) is the effective field theory describing the long- distance properties of lattice QCD with Wilson or twisted-mass fermions. We consider here WChPT for the theory with two light flavors of Wilson fermions or a single light twisted-mass fermion. Discretization errors introduce three low energy constants (LECs) into partially quenched WChPT at O(a^2), conventionally called W'_6, W'_7 and W'_8 . The phase structure of the theory at non-zero a depends on the sign of the combination 2W'_6 + W'_8, while the spectrum of the lattice Hermitian Wilson-Dirac operator depends on all three constants. It has been argued, based on the positivity of partition functions of fixed topological charge, and on the convergence of graded group integrals that arise in the epsilon-regime of ChPT, that there is a constraint on the LECs arising from the underlying lattice theory. In particular, for W'_6 = W'_7 = 0, the constraint found is W'_8 \le 0. Here we provide an alternative line of argument, based on mass inequalities for the underlying partially quenched theory. We find that W'_8 \le 0, irrespective of the values of W'_6 and W'_7. Our constraint implies that 2W'_6 > |W'_8| if the phase diagram is to be described by the first-order scenario, as recent simulations suggest is the case for some choices of action.
We report on recent progress in the calculation of $B_K$ using HYP-smeared staggered fermions on the MILC asqtad lattices. Our main focus is on the continuum extrapolation, which is done using (up to) four different lattice spacings---$a\approx$ 0.12, 0.09, 0.06 and 0.045 fm. Since Lattice 2010, we have reduced the statistical errors on the $a\approx 0.09\;$fm lattices by a factor of $\sim 3$, and roughly doubled the size of the $a\approx0.045\;$fm ensemble. We find that these improvements have a very significant impact on the continuum extrapolation, with the $a\approx 0.12\;$fm data lying outside the range of applicability of simple functional forms. Hence we use only the three smallest lattice spacings to perform the extrapolation, finding $\hat{B}_K = B_K(\text{RGI}) = 0.725 \pm 0.004(\text{stat}) \pm 0.038(\text{sys}) $. This value is consistent with our published value from 2010 (based the three coarsest lattice spacings), but has smaller errors.
In the low energy effective theory describing the partially quenched extension of two light Wilson fermions, three low energy constants (LECs) appear in terms proportional to a^2 (a being the lattice spacing). We propose methods to separately calculate these LECs, typically called W'_6, W'_7 and W'_8. While only one linear combination of these constants enters into physical quantities, different combinations enter into the description of the spectral density and eigenvalue distributions of the lattice Dirac operator and its Hermitian counterpart. Thus it is useful to be able to determine the LECs separately. Our methods require studying certain correlation functions for either two or three pion scattering, which are accessible only in the partially quenched extension of the theory.
We extend our recent unquenched ($N_f=2+1$ flavor) calculation of $B_K$ using improved staggered fermions by including in the fits the finite volume shift predicted by one-loop staggered chiral perturbation theory. The net result is to lower the result in the continuum limit by 0.6%. This shift is slightly smaller than our previous estimate of finite volume effects based on a direct comparison between different volumes. To include the finite volume effects in a reasonable time, we found it necessary to calculate them using Graphics Processing Units.
We discuss finite volume errors in our calculations of $B_K$ using improved staggered fermions on the MILC asqtad lattices. Using GPUs, we are now able to extrapolate using next-to-leading order (NLO) staggered SU(2) chiral perturbation theory including the finite volume corrections arising from pion loops. We find that the impact of FV fitting is very small, giving a 0.5% shift in the continuum limit.
We study the single site SU(N) lattice gauge theory with N_f=2 adjoint Wilson fermions for values of N up to 53. We determine the phase diagram of the theory as a function of the hopping parameter kappa and the inverse 't Hooft coupling b, searching for the region in which the Z_N^4 center symmetry is unbroken. In this region the theory is equivalent to the infinite volume theory when N goes to infinity. We find a region of values of kappa on both sides of kappa_c for which the symmetry is unbroken, including both light physical quarks and masses ~O(1/a). This is surrounded by a region with a complicated sequence of partially broken phases. We calculate Wilson loop expectation values and find that using N <= 53 it is possible to extract the heavy-quark potential at small distances (1-3 links) but not at longer distances. For this, larger values of N, or lattices with more sites, are needed.
We present results for non-perturbative renormalization (NPR) factors for staggered fermion bilinears of arbitrary spin and taste. We use "covariant" bilinears which transform irreducibly under the lattice translation and rotation group, and thus do not mix. We form ~ 30 ratios which have no anomalous dimensions, and compare the NPR results to those from 1-loop perturbation theory. We also compare the absolute renormalization factors (which, in general, do have anomalous dimensions) to 1-loop perturbation theory. We use asqtad and HYP-smeared staggered valence fermions on the coarse MILC asqtad lattices.
We present results for one-loop matching factors of four-fermion operators composed of HYP-smeared staggered fermions. We generalize previous calculations by using the tree-level improved Symanzik gauge action. These results are needed for our companion numerical calculation of $B_K$ and related matrix elements. We find that the impact on one-loop matching factors of using the improved gluon action is much smaller than that from the use of either HYP smearing or mean-field improvement. The one-loop coefficients for mean-field improved, HYP-smeared operators with the Symanzik gauge action have a maximum magnitude of $O(1)\times α_s$, indicating that perturbation theory is reasonably convergent.
We present results for kaon mixing parameter $B_K$ calculated using HYP-smeared improved staggered fermions on the MILC asqtad lattices. We use three lattice spacings ($a\approx 0.12$, $0.09$ and $0.06\;$fm), ten different valence quark masses ($m\approx m_s/10-m_s$), and several light sea-quark masses in order to control the continuum and chiral extrapolations. We derive the next-to-leading order staggered chiral perturbation theory (SChPT) results necessary to fit our data, and use these results to do extrapolations based both on SU(2) and SU(3) SChPT. The SU(2) fitting is particularly straightforward because parameters related to taste-breaking and matching errors appear only at next-to-next-to-leading order. We match to the continuum renormalization scheme (NDR) using one-loop perturbation theory. Our final result is from the SU(2) analysis, with the SU(3) result providing a (less accurate) cross check. We find $B_K(\text{NDR}, μ= 2 \text{GeV}) = 0.529 \pm 0.009 \pm 0.032$ and $\hat{B}_K =B_K(\text{RGI})= 0.724 \pm 0.012 \pm 0.043$, where the first error is statistical and the second systematic. The error is dominated by the truncation error in the matching factor. Our results are consistent with those obtained using valence domain-wall fermions on lattices generated with asqtad or domain-wall sea quarks.
We report updated results for $B_K$ calculated using HYP-smeared staggered fermions on the MILC asqtad 2+1 flavor lattices. We use four different lattice spacings ($a \approx$ 0.12, 0.09, 0.06 and 0.045 fm) to control the continuum extrapolation. We use SU(2) staggered chiral perturbation theory to do the data analysis. We find that $B_K(\text{NDR}, μ=2 \text{GeV}) = 0.526 \pm 0.007 \pm 0.024$ and $\hat{B}_K = B_K(\text{RGI}) = 0.720 \pm 0.010 \pm 0.033$. Here the first error is statistical and the second systematic. The dominant source of error is that due to our use of a truncated (one-loop) matching factor.
We report updated results for $B_K$ using HYP-smeared staggered valence quarks on MILC asqtad lattices based on an analysis using SU(3) staggered chiral perturbation theory. The most important new feature of our data sample is the inclusion of a fourth ("ultrafine") lattice spacing. This improves the control over the continuum extrapolation and errors due our use of one-loop perturbative matching. We present a complete updated error budget, which leads to $B_K(\text{NDR}, μ= 2 \text{GeV}) = 0.5309 \pm 0.0051 \pm 0.0424$ and $\hat{B}_K = B_K(\text{RGI}) = 0.727 \pm 0.07 \pm 0.058$. The results of the SU(3) analysis are inferior to those based on SU(2) staggered chiral perturbation theory, primarily because of the dependence on the Bayesian priors we use in the SU(3) fits.
We present preliminary results for $B_K$ from the MILC ultrafine lattices, based on a partial ensemble of 305 configurations. We use HYP-smeared improved staggered valence quarks. The analysis is done using fitting forms based on both SU(2) and SU(3) staggered chiral perturbation thery. For the SU(2) analysis, we find that the result using the NLO fit function is consistent with that from a partial NNLO fit. For the SU(3) analysis, where we have to use partially constrained fits due to the number of fit parameters, we find that our two preferred fits ("N-BB1" and "N-BB2") are also consistent, both with each other and with the results of the SU(2) fits. These results are used in companion proceedings to improve the control over the continuum extrapolation.
We present results for matching factors for staggered four-fermion operators constructed using HYP-smeared fat links both in the action and the operators. We use perturbation theory to calculate the matching factors and work to one-loop order. The new feaure of this work is the use of the Symanzik-improved gauge action, as opposed to the Wilson gauge action. Our results are needed for our ongoing calculation of weak matrix elements using HYP-smeared staggered valence quarks and operators on MILC lattices. We give explicit results for matching factors of the operator needed to calculate $B_K$. We compare the impact of the improvement of the gauge action on one-loop coefficients with that of mean-field improvement of the operators.
We present results for one-loop perturbative matching factors using bilinear operators composed of improved staggered fermions, using unimproved (Wilson) and improved (Symanzik, Iwasaki, and DBW2) gluon actions. We consider two fermions actions---HYP/$\bar{\text{Fat7}}$-smeared and "asqtad". The former is being used in calculations of electroweak matrix elements, while the latter have been used extensively by the MILC collaboration. We observe that using the improved gluon action leads to small reductions in the perturbative corrections, but that these reductions are smaller than those obtained when moving from the tadpole-improved naive staggered action to either HYP-smeared or asqtad action.
We study three sources of error in our calculation of $B_K$ using HYP-smeared staggered fermions on the MILC asqtad lattices. These are (1) dependence on the light sea quark mass; (2) finite volume effects; and (3) the impact of an order of magnitude increase in the number of measurements. Our main results are (1) the dependence on the light sea-quark mass is weaker than expected by naive dimensional analysis, (2) including finite volume effects in SU(2) staggered chiral perturbation theory fits leads to a very small change in $B_K$, of size $\approx 0.1%$, and (3) increasing the statistics on one of the coarse MILC lattices resolves a potential discrepancy with other coarse results.
We present the results for $B_K$ calculated using HYP-smeared staggered fermions using lattices generated by the MILC collaboration using the asqtad staggered action. We have done the calculation on 8 ensembles of these lattices, including three different lattice spacings ($a=0.12, 0.09, 0.06$ fm). We fit the data to forms based on those predicted by SU(3) mixed-action partially-quenched staggered chiral perturbation theory. Our preliminary result is $ B_K(\text{NDR},μ=2\text{GeV}) = 0.528 \pm 0.011 \pm 0.048 $, where the first error is statistical and the second systematic. The error turns out to be larger than that from an analysis using SU(2) chiral perturbation theory.
We present results for $B_K$ calculated using HYP-smeared improved staggered fermions on the MILC asqtad lattices. In this report, the data is analyzed using the results of SU(2) staggered chiral perturbation theory (SChPT). We outline the derivation of the NLO SU(2) SChPT result, explain our fitting procedure, and outline how we estimate systematic errors. We also show the light sea-quark mass and lattice spacing dependence for both SU(2) and SU(3)-based analyses. Our preliminary result from the SU(2) analysis is $B_K(\text{NDR}, μ= 2 \text{GeV}) = 0.512 \pm 0.014 \pm 0.034$ and $\hat{B}_K = B_K(\text{RGI})= 0.701 \pm 0.019 \pm 0.047$. This is somewhat more accurate than our result from the SU(3) analysis. It is consistent with results obtained using valence domain-wall fermions
We study the finite-volume effects in our calculation of $B_K$ using HYP-smeared improved staggered valence fermions. We calculate the predictions of both SU(3) and SU(2) staggered chiral perturbation theory at one-loop order. We compare these to the results of a direct calculation, using MILC coarse lattices with two different volumes: $20^3$ and $28^3$. From the direct calculation, we find that the finite volume effect is $\approx 2%$ for the SU(3) analysis and $\approx 0.9%$ for the SU(2) analysis. We also show how the statistical error depends on the number of measurements made per configuration, and make a first study of autocorrelations.
We discuss the impact of using one-loop matching on the calculation of $B_K$ using HYP-smeared improved staggered fermions. We give estimates of size of the truncation errors from the missing two-loop corrections.
We present results for matching factors for bilinear operators composed of HYP-smeared staggered fermions and constructed using HYP-smeared fat links. The matching factors are calculated perturbatively at one-loop order. The new feature of our calculation compared to previous work on HYP-smeared staggered fermions is the use of the Symanzik-improved gluon propagator, which allows our results to be applied to our ongoing simulations based on configurations generated by the MILC collaboration. We address the issue of the relative efficiency of various improvement schemes in reducing one-loop corrections to the matching factors.
We reconsider the construction and matching of the vector and axial currents in Wilson Chiral Perturbation Theory (WChPT), the low-energy effective theory for lattice QCD with Wilson fermions. In particular, we discuss in detail the impact of the finite renormalization of the currents on their matching from the lattice theory to WChPT.We explicitly show that imposing chiral Ward-Takahashi identities on the currents leads, in general, to additional terms of O(a) in the axial current. We illustrate the impact on physical quantities by computing the pion decay constant to one-loop order in the two flavor theory. Our result differs from previously published ones.
We reconsider the construction of the vector and axial-vector currents in Wilson Chiral Perturbation Theory (WChPT), the low-energy effective theory for lattice QCD with Wilson fermions. We discuss in detail the finite renormalization of the currents that has to be taken into account in order to properly match the currents. We explicitly show that imposing the chiral Ward identities on the currents does, in general, affect the axial-vector current at O(a). As an application of our results we compute the pion decay constant to one loop in the two flavor theory. Our result differs from previously published ones.
It has been proposed that four-dimensional QCD with fermions in the adjoint representation exhibits volume-independence in the large-N limit. If correct, this would mean that results for physical quantities could be obtained from the single-site version of the model. A necessary condition for volume-independence is that the (Z_N)^4 center-symmetry of the single-site theory is unbroken. We explore the phase diagram of the theory with a single Dirac fermion using Wilson fermions for a number of colors in the range N=8-15, and identify the region in the parameter space of quark mass and gauge coupling where the symmetry appears to be unbroken. Our evidence suggests that this region includes both light and heavy quarks, and our results are consistent with this region extending to the continuum limit.
We study the validity of the large-N equivalence between four-dimensional SU(N) lattice gauge theory and its momentum quenched version -- the Quenched Eguchi-Kawai (QEK) model. We have found strong evidence that this equivalence does not hold in the weak-coupling regime (and thus in the continuum limit). This is based on weak-coupling analytic arguments and Monto-Carlo simulations at intermediate couplings with 20 <= N <= 200. Since detailed expositions of our arguments, methods and results have already appeared in Phys. Rev. D78:034507 (2008) and Phys. Rev. D78:074503 (2008), we attempt here to give a more intuitive explanation of our results. The breakdown of reduction that we find is due to a dynamically generated correlation between different Euclidean components of the gauge fields.
We give an update on our calculation of $B_K$ using HYP-smeared valence staggered quarks. We have results for $B_K$ at tree-level on several coarse MILC lattices ($a\approx 0.12 $fm) and one of the fine lattices ($a\approx 0.09 $fm), using 10 light valence quarks ranging down to $m_s^{\rm phys}/10$. We have generalized staggered chiral perturbation theory to our mixed action setup, and outline the results. We explain our present fitting strategy, and give some preliminary results.
We study the pion spectrum (and in particular taste-symmetry breaking within it) using HYP-smeared valence staggered fermions on the coarse and fine MILC lattices (which have asqtad staggered sea quarks). We focus on the dependence on lattice spacing and sea-quark mass. We also update our results on source dependence. Our main conclusion is that on the MILC fine lattices the appropriate power-counting for SU(3) staggered chiral perturbation theory may have discretization errors entering at next-to-leading order rather than at leading-order.
I give a status report on the validity of the so-called ``fourth-root trick'', i.e. the procedure of representing the determinant for a single fermion by the fourth root of the staggered fermion determinant. This has been used by the MILC collaboration to create a large ensemble of lattices using which many quantities of physical interest have been and are being calculated. It is also used extensively in studies of QCD thermodynamics. The main question is whether the theory so defined has the correct continuum limit. There has been significant recent progress towards answering this question. After recalling the issue, and putting it into a broader context of results from statistical mechanics, I critically review the new work. I also address the related issue of the impact of treating valence and sea quarks differently in rooted simulations, discuss whether rooted simulations at finite temperature and density are subject to additional concerns, and briefly update results for quark masses using the MILC configurations. An answer to the question in the title is proposed in the summary.
I make a number of observations concerning discretization errors in twisted-mass lattice QCD that can be deduced by applying chiral perturbation theory including lattice artifacts. (1) The line along which the PCAC quark mass vanishes in the twisted mass-twisted mass plane makes an angle to the untwisted mass axis which is a direct measure of O(a) terms in the chiral Lagrangian, and is found numerically to be large; (2) Numerical results for pionic quantities in the mass plane show the qualitative properties predicted by chiral perturbation theory, in particular an asymmetry in slopes between positive and negative untwisted quark masses; (3) By extending the description of the ``Aoki regime'' (where m_q is of size a^2 Lambda_QCD^3) to next-to-leading order in chiral perturbation theory I show how the phase transition lines and lines of maximal twist (using different definitions) extend into this region, and give predictions for the functional form of pionic quantities; (4) I argue that the recent claim that lattice artifacts at maximal twist have apparent infrared singularities in the chiral limit results from expanding about the incorrect vacuum state. Shifting to the correct vacuum (as can be done using chiral perturbation theory) the apparent singularities are summed into non-singular, and furthermore predicted, forms. I further argue that there is no breakdown in the Symanzik expansion in powers of lattice spacing, and no barrier to simulating at maximal twist in the Aoki regime.
We study taste and Euclidean rotational symmetry violation for staggered fermions at nonzero lattice spacing using staggered chiral perturbation theory. We extend the staggered chiral Lagrangian to O(a^2 p^2), O(a^4) and O(a^2 m), the orders necessary for a full next-to-leading order calculation of pseudo-Goldstone boson masses and decay constants including analytic terms. We then calculate a number of SO(4) taste-breaking quantities, which involve only a small subset of these NLO operators. We predict relationships between SO(4) taste-breaking splittings in masses, pseudoscalar decay constants, and dispersion relations. We also find predictions for a few quantities that are not SO(4) breaking. All these results hold also for theories in which the fourth-root of the fermionic determinant is taken to reduce the number of quark tastes; testing them will therefore provide evidence for or against the validity of this trick.
We study the improvement of staggered fermions using hypercubically smeared (HYP) links. We calculate the strange quark mass and the kaon B-parameter, $B_K$, in quenched QCD on a $16^3 \times 64$ lattice at $β=6.0$. We find $m_s(\bar{\rm MS},2 {\rm GeV})=101.2\pm1.3\pm4 $MeV and $B_K(\bar{\rm MS},2 {\rm GeV}) = 0.578 \pm 0.018\pm 0.042$, where the first error is from statistics and fitting, and the second from using one-loop matching factors. The scale ($1/a=1.95$GeV) is set by $M_ρ$, and $m_s$ is determined using the kaon mass. Comparing to quenched results obtained using unimproved staggered fermions and other discretizations, we argue that the size of discretization errors in $B_K$ is substantially reduced by improvement.
We have extended staggered chiral perturbation theory to O(a^2 p^2), O(a^4), and O(a^2 m), the orders necessary for a full next-to-leading order calculation of pseudo-Goldstone boson masses and decay constants including taste-symmetry violations. We present predictions relating SO(4) taste-breaking splittings in masses, pseudoscalar decay constants, and dispersion relations. These can be used to test the fourth-root trick.
This is an approximate reconstruction of the panel discussion on chiral extrapolation of physical observables. The session consisted of brief presentations from panelists, followed by responses from the panel, and concluded with questions and comments from the floor with answers from panelists. In the following, the panelists have summarized their statements, and the ensuing discussion has been approximately reconstructed from notes.
I discuss the properties of pions in ``partially quenched'' theories, i.e. those in which the valence and sea quark masses, $m_V$ and $m_S$, are different. I point out that for lattice fermions which retain some chiral symmetry on the lattice, e.g. staggered fermions, the leading order prediction of the chiral expansion is that the mass of the pion depends only on $m_V$, and is independent of $m_S$. This surprising result is shown to receive corrections from loop effects which are of relative size $m_S \ln m_V$, and which thus diverge when the valence quark mass vanishes. Using partially quenched chiral perturbation theory, I calculate the full one-loop correction to the mass and decay constant of pions composed of two non-degenerate quarks, and suggest various combinations for which the prediction is independent of the unknown coefficients of the analytic terms in the chiral Lagrangian. These results can also be tested with Wilson fermions if one uses a non-perturbative definition of the quark mass.
Some phase transitions of cosmological interest may be weakly first-order and cannot be analyzed by a simple perturbative expansion around mean field theory. We propose a simple two-scalar model--the cubic anisotropy model--as a foil for theoretical techniques to study such transitions, and we review its similarities and dissimilarities to the electroweak phase transition in the early universe. We present numerical Monte Carlo results for various discontinuities across very weakly first-order transitions in this model and, as an example, compare them to epsilon-expansion results. For this purpose, we have computed through next-to-next-to-leading order in epsilon.