Subtractive renormalization of the NN interaction in chiral effective theory up to next-to-next-to-leading order: S waves

We extend our subtractive-renormalization method in order to evaluate the 1S0 and 3S1-3D1 NN scattering phase shifts up to next-to-next-to-leading order (NNLO) in chiral effective theory. We show that, if energy-dependent contact terms are employed in the NN potential, the 1S0 phase shift can be obtained by carrying out two subtractions on the Lippmann-Schwinger equation. These subtractions use knowledge of the the scattering length and the 1S0 phase shift at a specific energy to eliminate the low-energy constants in the contact interaction from the scattering equation. For the J=1 coupled channel, a similar renormalization can be achieved by three subtractions that employ knowledge of the 3S1 scattering length, the 3S1 phase shift at a specific energy and the 3S1-3D1 generalized scattering length. In both channels a similar method can be applied to a potential with momentum-dependent contact terms, except that in that case one of the subtractions must be replaced by a fit to one piece of experimental data. This method allows the use of arbitrarily high cutoffs in the Lippmann-Schwinger equation. We examine the NNLO S-wave phase shifts for cutoffs as large as 5 GeV and show that the presence of linear energy dependence in the NN potential creates spurious poles in the scattering amplitude. In consequence the results are in conflict with empirical data over appreciable portions of the considered cutoff range. We also identify problems with the use of cutoffs greater than 1 GeV when momentum-dependent contact interactions are employed. These problems are ameliorated, but not eliminated, by the use of spectral-function regularization for the two-pion exchange part of the NN potentialComment: 40 pages, 21 figure

Charmonium properties from lattice QCD + QED: hyperfine splitting, J/ψJ/\psi leptonic width, charm quark mass and aμca_{\mu}^c

We have performed the first $n_f = 2+1+1$ lattice QCD computations of the properties (masses and decay constants) of ground-state charmonium mesons. Our calculation uses the HISQ action to generate quark-line connected two-point correlation functions on MILC gluon field configurations that include $u/d$ quark masses going down to the physical point, tuning the $c$ quark mass from $M_{J/\psi}$ and including the effect of the $c$ quark's electric charge through quenched QED. We obtain $M_{J/\psi}-M_{\eta_c}$ (connected) = 120.3(1.1) MeV and interpret the difference with experiment as the impact on $M_{\eta_c}$ of its decay to gluons, missing from the lattice calculation. This allows us to determine $\Delta M_{\eta_c}^{\mathrm{annihiln}}$ =+7.3(1.2) MeV, giving its value for the first time. Our result of $f_{J/\psi}=$ 0.4104(17) GeV, gives $\Gamma(J/\psi \rightarrow e^+e^-)$=5.637(49) keV, in agreement with, but now more accurate than experiment. At the same time we have improved the determination of the $c$ quark mass, including the impact of quenched QED to give $\overline{m}_c(3\,\mathrm{GeV})$ = 0.9841(51) GeV. We have also used the time-moments of the vector charmonium current-current correlators to improve the lattice QCD result for the $c$ quark HVP contribution to the anomalous magnetic moment of the muon. We obtain $a_{\mu}^c = 14.638(47) \times 10^{-10}$, which is 2.5$\sigma$ higher than the value derived using moments extracted from some sets of experimental data on $R(e^+e^- \rightarrow \mathrm{hadrons})$. This value for $a_{\mu}^c$ includes our determination of the effect of QED on this quantity, $\delta a_{\mu}^c = 0.0313(28) \times 10^{-10}$.Comment: Added extra discussion on QED setup, some new results to study the effects of strong isospin breaking in the sea (including new Fig. 1) and a fit stability plot for the hyperfine splitting (new Fig. 7). Version accepted for publication in PR

QCD Factorization for B→ππB\to\pi\pi Decays: Strong Phases and CP Violation in the Heavy Quark Limit

We show that, in the heavy quark limit, the hadronic matrix elements that enter $B$ meson decays into two light mesons can be computed from first principles, including `non-factorizable' strong interaction corrections, and expressed in terms of form factors and meson light-cone distribution amplitudes. The conventional factorization result follows in the limit when both power corrections in $1/m_b$ and radiative corrections in $\alpha_s$ are neglected. We compute the order-$\alpha_s$ corrections to the decays $B_d\to\pi^+\pi^-$, $B_d\to\pi^0\pi^0$ and $B^+\to\pi^+\pi^0$ in the heavy quark limit and briefly discuss the phenomenological implications for the branching ratios, strong phases and CP violation.Comment: 6 pages, 1 figur

Mesonic decay constants in lattice NRQCD

Lattice NRQCD with leading finite lattice spacing errors removed is used to calculate decay constants of mesons made up of heavy quarks. Quenched simulations are done with a tadpole improved gauge action containing plaquette and six-link rectangular terms. The tadpole factor is estimated using the Landau link. For each of the three values of the coupling constant considered, quarkonia are calculated for five masses spanning the range from charmonium through bottomonium, and one set of quark masses is tuned to the B(c). "Perturbative" and nonperturbative meson masses are compared. One-loop perturbative matching of lattice NRQCD with continuum QCD for the heavy-heavy vector and axial vector currents is performed. The data are consistent with the vector meson decay constants of quarkonia being proportional to the square root of their mass and the B(c) decay constant being equal to 420(13) MeV.Comment: 25 pages in REVTe

Update: Accurate Determinations of alpha_s from Realistic Lattice QCD

We use lattice QCD simulations, with MILC configurations (including vacuum polarization from u, d, and s quarks), to update our previous determinations of the QCD coupling constant. Our new analysis uses results from 6 different lattice spacings and 12 different combinations of sea-quark masses to significantly reduce our previous errors. We also correct for finite-lattice-spacing errors in the scale setting, and for nonperturbative chiral corrections to the 22 short-distance quantities from which we extract the coupling. Our final result is alpha_V(7.5GeV,nf=3) = 0.2120(28), which is equivalent to alpha_msbar(M_Z,n_f=5)= 0.1183(8). We compare this with our previous result, which differs by one standard deviation.Comment: 12 pages, 2 figures, 4 table

A Study of Gluon Propagator on Coarse Lattice

We study gluon propagator in Landau gauge with lattice QCD, where we use an improved lattice action. The calculation of gluon propagator is performed on lattices with the lattice spacing from 0.40 fm to 0.24 fm and with the lattice volume from $(2.40 fm)^4$ to $(4.0 fm)^4$. We try to fit our results by two different ways, in the first one we interpret the calculated gluon propagators as a function of the continuum momentum, while in the second we interpret the propagators as a function of the lattice momentum. In the both we use models which are the same in continuum limit. A qualitative agreement between two fittings is found.Comment: Revtex 14pages, 11 figure

Subtractive renormalization of the chiral potentials up to next-to-next-to-leading order in higher NN partial waves

We develop a subtractive renormalization scheme to evaluate the P-wave NN scattering phase shifts using chiral effective theory potentials. This allows us to consider arbitrarily high cutoffs in the Lippmann-Schwinger equation (LSE). We employ NN potentials computed up to next-to-next-to-leading order (NNLO) in chiral effective theory, using both dimensional regularization and spectral-function regularization. Our results obtained from the subtracted P-wave LSE show that renormalization of the NNLO potential can be achieved by using the generalized NN scattering lengths as input--an alternative to fitting the constant that multiplies the P-wave contact interaction in the chiral effective theory NN force. However, in order to obtain a reasonable fit to the NN data at NNLO the generalized scattering lengths must be varied away from the values extracted from the so-called high-precision potentials. We investigate how the generalized scattering lengths extracted from NN data using various chiral potentials vary with the cutoff in the LSE. The cutoff-dependence of these observables, as well as of the phase shifts at $T_{lab} \approx 100$ MeV, suggests that for a chiral potential computed with dimensional regularization the highest LSE cutoff it is sensible to adopt is approximately 1 GeV. Using spectral-function regularization to compute the two-pion-exchange potentials postpones the onset of cutoff dependence in these quantities, but does not remove it.Comment: 27 pages, 14 figure

Direct determination of the strange and light quark condensates from full lattice QCD

We determine the strange quark condensate from lattice QCD for the first time and compare its value to that of the light quark and chiral condensates. The results come from a direct calculation of the expectation value of the trace of the quark propagator followed by subtraction of the appropriate perturbative contribution, derived here, to convert the non-normal-ordered mψ̅ ψ to the MS̅ scheme at a fixed scale. This is then a well-defined physical “nonperturbative” condensate that can be used in the operator product expansion of current-current correlators. The perturbative subtraction is calculated through O(αs) and estimates of higher order terms are included through fitting results at multiple lattice spacing values. The gluon field configurations used are “second generation” ensembles from the MILC collaboration that include 2+1+1 flavors of sea quarks implemented with the highly improved staggered quark action and including u/d sea quarks down to physical masses. Our results are ⟨s̅ s⟩MS̅ (2  GeV)=-(290(15)  MeV)3, ⟨l̅ l⟩MS̅ (2  GeV)=-(283(2)  MeV)3, where l is a light quark with mass equal to the average of the u and d quarks. The strange to light quark condensate ratio is 1.08(16). The light quark condensate is significantly larger than the chiral condensate in line with expectations from chiral analyses. We discuss the implications of these results for other calculations

Light pseudoscalar decay constants, quark masses, and low energy constants from three-flavor lattice QCD

As part of our program of lattice simulations of three flavor QCD with improved staggered quarks, we have calculated pseudoscalar meson masses and decay constants for a range of valence quark masses and sea quark masses on lattices with lattice spacings of about 0.125 fm and 0.09 fm. We fit the lattice data to forms computed with staggered chiral perturbation theory. Our results provide a sensitive test of the lattice simulations, and especially of the chiral behavior, including the effects of chiral logarithms. We find: f_\pi=129.5(0.9)(3.5)MeV, f_K=156.6(1.0)(3.6)MeV, and f_K/f_\pi=1.210(4)(13), where the errors are statistical and systematic. Following a recent paper by Marciano, our value of f_K/f_\pi implies |V_{us}|=0.2219(26). Further, we obtain m_u/m_d= 0.43(0)(1)(8), where the errors are from statistics, simulation systematics, and electromagnetic effects, respectively. The data can also be used to determine several of the constants of the low energy effective Lagrangian: in particular we find 2L_8-L_5=-0.2(1)(2) 10^{-3} at chiral scale m_\eta. This provides an alternative (though not independent) way of estimating m_u; 2L_8-L_5 is far outside the range that would allow m_u=0. Results for m_s^\msbar, \hat m^\msbar, and m_s/\hat m can be obtained from the same lattice data and chiral fits, and have been presented previously in joint work with the HPQCD and UKQCD collaborations. Using the perturbative mass renormalization reported in that work, we obtain m_u^\msbar=1.7(0)(1)(2)(2)MeV and m_d^\msbar=3.9(0)(1)(4)(2)MeV at scale 2 GeV, with errors from statistics, simulation, perturbation theory, and electromagnetic effects, respectively.Comment: 86 pages, 22 figures. v3: Remarks about m_u=0 and the strong CP problem modified; reference added. Figs 5--8 modified for clarity. Version to be published in Phys. Rev. D. v2: Expanded discussion of finite volume effects, normalization in Table I fixed, typos and minor errors correcte

Tadpole-improved SU(2) lattice gauge theory

A comprehensive analysis of tadpole-improved SU(2) lattice gauge theory is made. Simulations are done on isotropic and anisotropic lattices, with and without improvement. Two tadpole renormalization schemes are employed, one using average plaquettes, the other using mean links in Landau gauge. Simulations are done with spatial lattice spacings $a_s$ in the range of about 0.1--0.4 fm. Results are presented for the static quark potential, the renormalized lattice anisotropy $a_t/a_s$ (where $a_t$ is the ``temporal'' lattice spacing), and for the scalar and tensor glueball masses. Tadpole improvement significantly reduces discretization errors in the static quark potential and in the scalar glueball mass, and results in very little renormalization of the bare anisotropy that is input to the action. We also find that tadpole improvement using mean links in Landau gauge results in smaller discretization errors in the scalar glueball mass (as well as in the static quark potential), compared to when average plaquettes are used. The possibility is also raised that further improvement in the scalar glueball mass may result when the coefficients of the operators which correct for discretization errors in the action are computed beyond tree level.Comment: 14 pages, 7 figures (minor changes to overall scales in Fig.1; typos removed from Eqs. (3),(4),(15); some rewording of Introduction