arxiv: v1 [hep-lat] 24 Dec 2008

Transcription

1 of hadrons from improved staggered quarks in full QCD arxiv: v1 [hep-lat] 4 Dec 008, a A. Bazavov, b C. Bernard, c C. DeTar, d W. Freeman, b Steven Gottlieb, a U.M. Heller, e J.E. Hetrick, f J. Laiho, c L. Levkova, d J. Osborn, g R. Sugar, h and D. Tosaint b (MILC Collaboration) a Department of Physics, Indiana University Bloomington, IN 47405, USA b Physics Department, University of Arizona Tucson, AZ 8571, USA c Physics Department, Washington University St. Louis, MO 610, USA d Physics Department, University of Utah Salt Lake City, UT 8411, USA e American Physical Society One Research Road, Box 9000, Ridge, NY , USA f Physics Department, University of the Pacific Stockton, CA 9511, USA g Argonne Leadership Computing Facility, Argonne National Laboratory Argonne, IL 6049, USA h Physics Department, University of California Santa Barbara, CA 9106, USA sbasak@indiana.edu We present our initial sty of the electromagnetic splittings of charged and neutral mesons, and the violation of Dashen s theorem. Hadron maes are calculated on MILC N f = +1 QCD ensembles at lattice spacing 0.15 fm, together with quenched non-compactu(1) configurations. The O(a ) tadpole improved staggered quark (asqtad) action is ed both for the sea quarks and for six different valence quark maes. Chiral extrapolations are performed ing partially quenched chiral perturbation theory incorporating electromagnetic corrections. The XXVI International Symposium on Lattice Field Theory July 14-19, 008 Williamsburg, Virginia, USA Speaker. Presently at NISER, Bhubaneswar, Oria , India c Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence.

2 1. Introduction Lattice QCD has long been ed to determine nonperturbative hadronic properties; the spectrum, in particular, has been extensively investigated. However, electromagnetic and isospinbreaking effects are ually not incled in those investigations. Although the strong interaction is blind to the electromagnetic charges, we do in fact have ma splittings of mesons and baryons in the real world which depend on both isospin-breaking and electromagnetic interactions. Further, it has been pointed out that the current evaluation of the light quark maes, particularly the ratio m u /m d, suffers significant uncertainties coming from the electromagnetic contributions to the maes of π and especially K mesons [1]. It is therefore important to incle electromagnetic effects in the lattice QCD simulations for more realistic spectrum calculations. The pioneering attempt to calculate charged and neutral pion splittings and the light quark maes from lattice QCD with electromagnetic interactions was by Duncan, Eichten and Thacker []. The photon fields were introduced in non-compact form and were treated in the quenched approximation. The pseoscalar meson maes were calculated with the Wilson action and at different values of electric charge. Recently, Blum et al. [] have calculated electromagnetic contributions to π and K mesons ma splittings and determined the light quark maes with domain wall quarks and N f = QCD configurations at the physical value of electric charge. They have also estimated the correction to the Dashen s theorem [4] at O(αm q ), and found em = 0.7(40) or 0.64(4) depending on their fitting range. Electromagnetic splittings of π and ρ mesons have also been calculated ing a RG-improved gauge action and a mean-field improved clover quark action at two different volumes and three lattice spacings in Ref. [5]. Electromagnetic effects on baryons have been discsed in Refs. [6, 7]; and calculation with unquenched photon fields in Ref. [8]. In this work, we sty the electromagnetic ma splittings of pseoscalar mesons in the presence of a quenched electromagnetic field and calculate the correction to the Dashen s theorem, parametrized following the chiral perturbation theory calculation of Urech [9]. The lattice data is fitted ing partially quenched chiral perturbation theory formulas at O(p 4, e m q ) given by Bijnens et al. [10]. We have ed a non-compact action for the photon fields following []. The valence quarks are asqtad [11] staggered quarks, and the configurations are the MILC N f = + 1 QCD configurations at lattice spacing 0.15 fm, with varying sea quark maes [1].. Lattice QCD with electromagnetic interaction The ma differences among the members of hadron isomultiplets arise from two sources: (a) strong isospin-breaking contributions due to the difference in maes of the light quarks and (b) different electromagnetic charges of the quarks. Apart from a small isospin-breaking contribution of order (m d m u ), the π ± π 0 ma splitting is almost purely electromagnetic in origin; whereas it is the isospin-breaking contributions that dominate for K 0 K ±. Dashen s theorem [4] summarizes the electromagnetic effects on meson maes: in the chiral limit to O(e ), electromagnetic interactions modify the leading order M π ± and M K ± maes while the π 0 and K 0 maes remain

3 unaffected, M π ± = ˆmB 0 + A (1) e (q u q d ), M K ± = ( ˆm+m s)b 0 + A (1) e (q u q s ) ( M π ± M π 0 )em = ( M K ± M K 0 )em, (.1) where q u, q d, and q s are the u, d, and s quark charges in units of e. At O(e m), however, M π ±,0 and M K ±,0 can receive large and different contributions, and the correction to Dashen s theorem can parametrized as, M D = M K M π = ( M K ± M K 0 ) em ( M π ± M π 0 )em, (.) or as E, defined by MK = (1+ E) Mπ. The partially quenched chiral perturbation theory relevant for QCD + QED with +1 dynamical flavors at NLO has recently been worked out by Bijnens et al. [10]; it can be ed to perform fits to the lattice Mπ,K data in order to determine the electromagnetic low energy constants. The pure electromagnetic correction, relevant for calculating the correction to Dashen s theorem, can be computed from the expreion, M em = M ps(χ 1, χ,q 1,q ) M ps(χ 1, χ,q,q ) M ps(χ 1, χ 1,q 1,q )+M ps(χ 1, χ 1,q,q ). (.) In this notation [10], the normalized quark ma is χ i = B 0 m i, where B 0 is related to the quarkantiquark vacm expectation value in the chiral limit, and Mps (χ 1, χ,q 1,q ) denotes the squared pseoscalar meson ma having valence quark maes χ 1, χ and valence charges q 1,q in units of electron charge e. In the isospin limit m u = m d, when χ 1 = m u, χ = m s and the quark charges q 1 = q u, q = q s, we have MD = M em. We fit the lattice data for MD according to MD = A 1(χ 1 χ 11 )+A [χ 1 log ( ) χ1 f + A [χ 1 f log µ χ f log + A 4 [ χ 1 { 1 1 log ( χ1 µ ( χ1 µ ( χ f µ ) χ 11 log )] ( )] χ11 µ )} χ 11 { 1 1 log ( χ11 µ )}], (.4) where, χ i j (χ i + χ j )/ and A i are the constants to be determined from the fit. For the scale µ we have ed 1000MeV. The sea-quark dependence of the log-terms appears through χ i f, where the sea-quark index f is summed over all the sea quarks. Terms with the sea-quark charges contributing to MD do not involve unknown low energy constants (LECs) at this order and hence are computable without lattice simulation [10]. In our expreion we have incled the parameters A,,4, which are not present in Ref. [10], in front of the logarithm terms. This is becae we expect the continm expreion to be modified by discretization effects on the quite coarse lattices we are ing. Finite volume effects, due to the malene of the photon, may also be rather significant. At this stage, therefore, we can obtain at best a rough estimate of the quantity MD ; a precise determination will require lattice data at finer lattice spacings and larger volumes.

4 . Numerical results Phenomenologically relevant quantities, such as the ma-squared differences of charged and neutral mesons and the correction to Dashen s theorem, can be calculated from partially quenched lattice QCD where photons are treated in the quenched approximation. Therefore, we can e the existing SU() configurations generated without dynamical photons. For unquenched QCD gluon configurations, we have ed four ensembles of existing a 0.15fm MILC lattices with +1 dynamical flavors, the details of which are provided in Table 1. The lattice spacing is set by the static quark potential [1]. Size β = 10/g a ˆm /am s a (fm) # Cfgs / / / / Table 1: N f = +1 MILC lattices ed in this project. The quenched photon configurations {A µ (n)} are generated from the non-compact U(1) action, S em = 1 4 nµν ( r µ A ν (n) r νa µ (n) ), (.1) subjected to Coulomb gauge fixing i la i(n) = 0. Aitional global gauge fixing is done to ensure that Gas s law is satisfied. In momentum space, the action is Gasian distributed, and the scalar and vector potentials [A 0 (p), A(p)] can be generated independently of each other from Gasian distributed random numbers. The coordinate-space Coulomb-gauge photon configurations are recovered by FFT. The neceary valence quark propagators are calculated with wall sources in these SU() U(1) background fields. We calculate the pseoscalar meson propagator with nine different valence quark maes, 0.1m s m q 1.0m s. All the fits we report here are obtained by correlated χ -fit; the errors are obtained from a jackknife analysis. The meson maes are extracted from the exponential fall-off of the meson propagators in the time range 9 4, taking into account the correlations among the time slices. We have ignored the contribution from the disconnected diagram that effects the neutral pion ma. We first try to estimate the O(e ) and O(e m) contribution to pseoscalar maes at three different (physical and non-physical) values of electric charge, M π (e 0) M π (e = 0) = A 0e (q u q d ) +O(e m). (.) In Fig. 1, we plot the dependence of π ma-squared splittings corresponding to Eq. (.). A straight line fit describes the O(e ) behavior fairly well over the full range of electric charges that are examined. In a further test to ascertain the O(e m) or higher order effects, in Fig., we plot the variation of the same pion ma-squared splittings against quark maes. The slope of each of line, corresponding to different electric charges, nicely matches to its e value. 4

5 m π (e 0) - mπ (e=0) (lattice unit) m u = m d = m u = m d = e /4π Figure 1: The e dependence of ma-squared electromagnetic splittings of the pion. m π (e 0)- mπ (e=0) (lattice unit) e = 0.0 e = e = m u + m d Figure : The O(e m) behavior of ma-squared electromagnetic splittings of the pion. In a series of four plots in Fig., we show the meson ma splittings between mesons computed with physical electric charge e = 4πα em and with e = 0. The labels q 1 q in the plots denote quark charges, i.e., indicates quarks with charge q u = e/ and q d = 1e/. The points ing the same symbol and color are obtained by varying the q quark ma. The lines are merely to guide the eye. We have not yet fit the squared meson maes to the Bijnens form based on partially quenched chiral perturbation theory (and hence not yet extracted the electromagnetic LECs) [10]. These LECs will be important in determining the physical isomultiplet splittings in mesons and therefore in calculating quark maes. The electromagnetic LECs can be ed to obtain the correction to Dashen s theorem, but we have estimated this correction directly by chiral extrapolation of Eq. (.) ing Eq. (.4). We put q 1 = q u, q = q s in Eq. (.) and extrapolate/interpolate to χ 1 = B 0 m u, χ = B 0 m s. Figure 4 is our preliminary result for the violation of Dashen s theorem, where we have plotted the M D for two ensembles, β = and 6.57, as a function of M π, obtained with e = 0, in physical units. The data from the remaining two QCD ensembles at heavier quark ma are still being analyzed. We have determined the coefficients A i in Eq. (.4) ing our lattice data for pseoscalar maes 5

6 q - mq1 (e=0) (lattice unit) q β = 6.566, m l /m s = 84/ q - mq1 (e=0) (lattice unit) q β = 6.57, m l /m s = / m q m q q - mq1 (e=0) (lattice unit) q β = 6.586, m l /m s = / q - mq1 (e=0) (lattice unit) q β = 6.600, m l /m s = 0.090/ m q m q Figure : The ma-squared splittings between pseoscalar mesons with physical value of electric charge e = 4πα em and e = 0. Each panel corresponds to different β and light sea-quark ma am l for fixed sea-strange quark ma am s = The labels q 1 q denote the quark charges. We write s instead of d whenever m q is equal to m s, the strange sea-quark ma in the simulation. χ i j, after which chiral extrapolation is performed at the experimental maes for π and K. Our preliminary result for the deviation from Dashen s theorem is < M D /(GeV ) < This may be compared with the estimate in Ref [10]: Although we find that our value is consistent with the phenomenological estimate, the error in our result is large and does not yet incle the systematic errors due to discretization and finite volume effects. We hope to reduce the error by increasing the statistics and extending the calculation to finer lattices and larger volumes. Fitting the meson splittings in Fig. to determine the electromagnetic LECs directly may also provide better control over our errors. 4. Conclions and outlook In this work we have calculated the electromagnetic ma splittings of pseoscalar maes with O(a ) improved staggered fermions on a 0.015fm N f = +1 MILC SU() ensembles and Coulomb-gauge U(1) configurations. We have determined the correction to Dashen s theorem at the physical point (in the isospin limit), but with relatively large errors, i.e., a wide range for M D. Our determination of meson and baryon isomultiplet splittings and quark maes is in progre. Our immediate goal is ing finer lattices to reduce the error in the correction to Dashen s theorem enough to bring down the present error in m u /m d. 6

7 M D (GeV ) β = (m l /m s = 84/0.0484) β = 6.57 (m l /m s = /0.0484) m π (GeV ) Figure 4: Correction to Dashen s theorem, the difference of electromagnetic contributions between K and π maes, as a function of the LO π ma squared (equivalent to the π ma squared with e = 0). References [1] MILC Collaboration: C. Aubin et al., Phys. Rev. D70, (004) [hep-lat/040708]. [] A. Duncan, E. Eichten, H. Thacker, Phys. Rev. Lett. 76, 894 (1996) [hep-lat/960005]. [] T. Blum et al., Phys. Rev. D76, (007) [ (hep-lat)]. [4] R.F. Dashen, Phys. Rev. 18, 145 (1969). [5] ] Y. Namekawa, Y. Kikukawa, PoS (LAT005) 090, 005 [hep-lat/050910]. [6] A. Duncan, E. Eichten, H. Thacker, Phys. Lett. B409, 87 (1997) [hep-lat/96070]. [7] T. Doi et al., PoS (LAT006) 174, 006 [hep-lat/ ]. [8] A. Duncan, E. Eichten, R. Sedgewick, Phys. Rev. D71, (005) [hep-lat/ ]. [9] R. Urech, Nucl. Phys. B4, 4 (1995) [hep-ph/940541]. [10] J. Bijnens, N. Danielon, Phys. Rev. D75, (007) [hep-lat/061017]. [11] C. Bernard et al., Phys. Rev. D58, (1998) [hep-lat/971010]. [1] C. Bernard et al., PoS (LAT007) 090, 007 [ (hep-lat)]. 7