The nucleon mass and pion-nucleon sigma term from a chiral analysis of lattice QCD world data

Transcription

1 The nucleon mass and pion-nucleon sigma term from a chiral analysis of lattice QCD world data L. Alvarez-Ruso 1, T. Ledwig 1, J. Martin-Camalich, M. J. Vicente-Vacas 1 1 Departamento de Física Teórica and IFIC, Centro Mixto Universidad de Valencia-CSIC, E Valencia, Spain Department of Physics and Astronomy, University of Sussex, B1 9QH, Brighton, UK MEU 013, Rome, Italy

2 Outline 1 Introduction BχPT, Perturbative ucleon Mass, Chiral Extrapolator 3 Results: Fits to LQCD f = Ensembles 4 Results: Fits to LQCD f = +1 Ensembles

3 Introduction ucleon mass non-perturbative regime of QCD BχPT: Effective theory using the QCD chiral symmetry for m q =0 SU(): u, d quarks pions, nucleons, (13) SU(3): u, d, s quarks p.-scalar mesons, 8/10 baryons LQCD: Lattice simulations using the discretized QCD-action flavors with m u = m d = m +1 flavors with m u = m d = m and m s m s(phys)

4 Introduction > 11 LQCD collab. BχPT LQCD: fixing LECs χlogs + natural LECs Hellmann-Feynman theorem GOR: m M π m m M (m) = σ π = m (p) uu + dd (p) σ π : π scattering, LQCD, BχPT π scattering: σ π = MeV 59(7) MeV Alarcon et al. 01 LQCD: σ π = MeV [aive RQM: σ π = 3m = 1 MeV]

5 Introduction > 11 LQCD collab. BχPT LQCD: fixing LECs χlogs + natural LECs Hellmann-Feynman theorem GOR: m M π m m M (m) = σ π = m (p) uu + dd (p) σ π : π scattering, LQCD, BχPT π scattering: σ π = MeV 59(7) MeV Alarcon et al. 01 LQCD: σ π = MeV [aive RQM: σ π = 3m = 1 MeV]

6 Introduction > 11 LQCD collab. BχPT LQCD: fixing LECs χlogs + natural LECs Hellmann-Feynman theorem GOR: m M π m m M (m) = σ π = m (p) uu + dd (p) σ π : π scattering, LQCD, BχPT π scattering: σ π = MeV 59(7) MeV Alarcon et al. 01 LQCD: σ π = MeV [aive RQM: σ π = 3m = 1 MeV]

7 Introduction > 11 LQCD collab. BχPT LQCD: fixing LECs χlogs + natural LECs Hellmann-Feynman theorem GOR: m M π m m M (m) = σ π = m (p) uu + dd (p) σ π : π scattering, LQCD, BχPT quark mass contribution (strangeness) of the nucleon mass σ-terms dark matter searches [aive RQM: σ π = 3m = 1 MeV]

8 Perturbative ucleon Mass ucleon mass pole position of its full propagator at /p = M 1 /p M 0 Σ(/p) M = M 0 +Σ(M ) Expand Σ around bare mass M 0 Rewrite propagator as Σ(/p) = Σ(M 0 )+( /p M 0 )Σ (M 0 )+R( /p) 1 /p M 0 Σ(/p) = 1 /p M 0 Σ(M 0) 1 Σ (M 0 ) ucleon mass M and its residue Z R( /p) 1 Σ (M 0 ) 1 1 Σ (M 0 ) M = M 0 +Z Σ(M 0 )+Z R(M ) Z = 1 1 Σ (M 0 )

9 BχPT O(p 4 ) Self-Energy: Σ p 4( /p) Covariant baryon chiral perturbation theory up to order O(p 4 ) with explicit (13) degrees of freedom: Σ p 4 ( /p) = Σ () +Σ (3) ( /p)+σ (4) ( /p) ΣC Σ3 ΣC4 ΣT4 Σ Σ 3 Σ 4 Σ 4a Σ 4b Perturbation in the light scales: p/λ, m π /Λ, (M M )/Λ with Λ = 4πf π Fits: including (13) ( ) and excluding it ( / ) (EOMS) EOMS - Gegelia, Fuchs et al. (1999, 003) consistent (13) couplings - Pascalutsa (1998,001,03) L () - Long, Lensky (011) π

10 M (4) and Hellmann-Feynman Theorem BχPT ucleon mass M (m π ) BχPT Scalar Form factor σ π = m (p) uu + dd (p) ΣC Σ σ 1 1 σ Chiral structures from BχPT M (4) ( ) m π = M0 +Σ C +Σ C4 +Σ (3)+(4) loops σ (4) π Hellmann-Feynman Theorem ( ) ] m π = Z [σ C +σ C4 +σ (3)+(4) loops q =0 ( ) ( ) σ π m π = m π M mπ m π

11 Chiral Extrapolator M (4) ( ) M π = M 0 c 1 4Mπ + 1 αm4 π + c 1 8π fπ Mπ 4 ln M π M 0 ( ) M π, M 0, M 0, f π, g A, h A, c 1, c, c 3, c 1 +Σ (3)+(4) loops Parameters: 10 low energy constants Only 3 are fitted L (n) p 1 p p 3 p 4 fixed f π, g A, h A, M 0 c, c 3, c 1 fitted M 0 c 1 α = 4(8e 38 + e e 116 )

12 Fixed Parameters f π, g A experimental values LECs c and c 3 from BχPT fits to π scattering data Alarcon et al. 01 (13) decay width h A =.87 (13) mass up to O(p 3 ) c 1 0.9±0.4 GeV 1 ucleon fits to LQCD Fits are insensitive to variations in c, c 3 Changes in c 1 compensated by changes in α LQCD data from: PRD64(054506), PRD70(094505), PRL107(141601)

13 Discretization and Finite Volume Finite Spacing Corrections O(a) - improved LQCD actions M = M a=0 + c a a +O(a 3 ) c a action specific constant Wilson (1974), Sheikholeslami (1985), Luscher et al.(1996) Tiburzi (005) LQCD data given in the dimensionless form (amπ, am ) Finite Volume Corrections d 4 l (4π) 4 = dl0 π d l (π) 3 dl0 1 π L 3 n with l = π L n, n Z 3 Kahn et al. (004)

14 Ingredients covariant SU() BχPT p 4 EOMS explicit (13) (consistent couplings) 10 parameters, only 3 are fitted finite volume effects finite size effects scaling/correlation of LQCD data points f = +1 LQCD data points for m s m s(phys) strange quark effects integrated out into the LECs

15 Results for f = Ensembles Fitting to LQCD data from ( deg. light quarks): BGR (010), ETMC (011), Mainz (01), QCDSF (013) Data is given in dimensionless form: (am π, am ) Ratio r 0 /a is known for each collaboration r 0 Sommer-scale (1994) (am π, am ) (r 0 M π, r 0 M ) Condition on points: LM π > 3.8 and r 0 M π < 1.11 χ = i ) M (n) ( M π + (n) Σ ( M π, L) + c a ã d i ( M π, L σ i ) with M(n) = r 0M (n), M π = (r 0 M π ), Σ(n) = r 0Σ (n) Fixing r 0 self-consistently and recursively inside the fit r k 0 = M (n) ( r k 1 0 M π(phys) ) M (phys) until r k 0 r k 1 0 < fm

16 Results for f = Ensembles Fitting to LQCD data from ( deg. light quarks): BGR (010), ETMC (011), Mainz (01), QCDSF (013) Data is given in dimensionless form: (am π, am ) Ratio r 0 /a is known for each collaboration r 0 Sommer-scale (1994) (am π, am ) (r 0 M π, r 0 M ) Condition on points: LM π > 3.8 and r 0 M π < 1.11 χ = i ) M (n) ( M π + (n) Σ ( M π, L) + c a ã d i ( M π, L σ i ) with M(n) = r 0M (n), M π = (r 0 M π ), Σ(n) = r 0Σ (n) Fixing r 0 self-consistently and recursively inside the fit r k 0 = M (n) ( r k 1 0 M π(phys) ) M (phys) until r k 0 r k 1 0 < fm

17 Results for f = Ensembles Fitting to LQCD data from ( deg. light quarks): BGR (010), ETMC (011), Mainz (01), QCDSF (013) Data is given in dimensionless form: (am π, am ) Ratio r 0 /a is known for each collaboration r 0 Sommer-scale (1994) (am π, am ) (r 0 M π, r 0 M ) Condition on points: LM π > 3.8 and r 0 M π < 1.11 χ = i ) M (n) ( M π + (n) Σ ( M π, L) + c a ã d i ( M π, L σ i ) with M(n) = r 0M (n), M π = (r 0 M π ), Σ(n) = r 0Σ (n) Fixing r 0 self-consistently and recursively inside the fit r k 0 = M (n) ( r k 1 0 M π(phys) ) M (phys) until r k 0 r k 1 0 < fm

18 no explicit (13) Results for f = Ensembles with explicit (13) one σ π data point at 90 MeV QCDSF p 4 c 1 GeV 1 χ dof r 0 fm σ π MeV / /σ 1.18(14) (13) / σ 0.91(4) (3) /σ 0.77(9) (10) σ 0.80(1) ()

19 no explicit (13) Results for f = Ensembles with explicit (13) one σ π data point at 90 MeV QCDSF p 4 c 1 GeV 1 χ dof r 0 fm σ π MeV / /σ 1.18(14) (13) / σ 0.91(4) (3) /σ 0.77(9) (10) σ 0.80(1) ()

20 Results for f = +1 Ensembles Fitting to LQCD data from ( deg. light quarks + 1 heavy quark): BMW (01), HSC (009), LHPC (009), MILC (01), PLQCD (011), PACS (010), RBCUK-QCD (01) SU() BχPT s-quark effects integrated out into the LECs Data is given in dimensionless form (am π, am ) and a a-uncertainty points from same collaboration are correlated χ A = i [ M (n) ( M π ) +Σ (n) ( M π, L ) ( f i d i M π, L ) ] [ fi 1 + f i σ i σ fi ] χ B = T V 1 [ i = M (n) ( ) Mπ,i ( ( )] + c i a i +Σ (n) Mπ,i, L i ) d i Mπ,i, L i Correlation included: (A) through collaboration specific fit constants f i (B) through correlation matrix V

21 Results for f = +1 Ensembles Fitting to LQCD data from ( deg. light quarks + 1 heavy quark): BMW (01), HSC (009), LHPC (009), MILC (01), PLQCD (011), PACS (010), RBCUK-QCD (01) SU() BχPT s-quark effects integrated out into the LECs Data is given in dimensionless form (am π, am ) and a a-uncertainty points from same collaboration are correlated χ A = i [ M (n) ( M π ) +Σ (n) ( M π, L ) ( f i d i M π, L ) ] [ fi 1 + f i σ i σ fi ] χ B = T V 1 [ i = M (n) ( ) Mπ,i ( ( )] + c i a i +Σ (n) Mπ,i, L i ) d i Mπ,i, L i Correlation included: (A) through collaboration specific fit constants f i (B) through correlation matrix V

22 Results for f = +1 Ensembles Condition on points: LM π > 3.8 and M π < 415 MeV p 4 c 1 GeV 1 χ dof σ π MeV / /a 1.15(3) (3) / a 1.10(5) 1. 55(3) /a 0.89(3) () a 0.84(4) 1. 44(3)

23 Results for f = +1 Ensembles Condition on points: LM π > 3.8 and M π < 415 MeV p 4 c 1 GeV 1 χ dof σ π MeV / /a 1.15(3) (3) / a 1.10(5) 1. 55(3) /a 0.89(3) () a 0.84(4) 1. 44(3)

24 Results f = vs f = +1 f = σ π = 41(5)(4) MeV incl. the σ(90) point f = +1 σ π = 5(3)(8) MeV average of and / results Chiral convergence (comparing only fits with ) σ π MeV f p p 3 p 4 full M () = M 0 4c 1 M π

25 Results f = vs f = +1 f = M f = +1 σ π M π [MeV] M [MeV] M π [GeV ] 41(5)(4) vs 5(3)(8): f = : more M data for < GeV f = +1: already one σ π data point for < 0.1 GeV

26 Summary Fitted SU() BχPT O(p 4 ) mass formula to LQCD nucleon mass data for and +1 flavor ensembles Considered scale-setting, finite volume and finite spacing effects Extracted LECs and calculated the σ π Correlation between c 1 and α f = : Either uncertainty of the σ(90) data-point is too small or inclusion of (13) is necessary f = +1: Consistent descriptions of M data by different strategies, small slope variations σ π changes by 10 MeV More f = M data points in the region m π < 0.1 GeV Direct f = +1 σ π LQCD points for mπ < 0.1 GeV

27 Chiral Extrapolator M (4) ( ) M π = M 0 c 1 4Mπ + 1 αm4 π + c 1 8π fπ Mπ 4 ln M π M 0 ( ) M π, M 0, M 0, f π, g A, h A, c 1, c, c 3, c 1 +Σ (3)+(4) loops Parameters: 10 low energy constants Only 3 are fitted L (n) p 1 p p 3 p 4 fixed f π, g A, h A, M 0 c, c 3, c 1 fitted M 0 c 1 α = 4(8e 38 + e e 116 ) Pion-mass O(p 4 ): α = α + c 16 1 f π l 3 r ( M 0 ) Mπ ( mπ ) (Λ = mπ + lr 3 ) f π0 mπ π f π0 mπ 4 ln m π Λ Shifting: M (mπ ) M (M π ) for σ π l 3 (139 MeV) = 3.(8), Colangelo (011)

28 Results f = vs f = +1 f = σ π = 41(3)(1) MeV σ π = 45(13)(10) MeV f = +1 σ π = 5(3)(8) MeV Chiral convergence σ π MeV full p p 3 p 4 f = σ f =