Soft and Hard Scale QCD Dynamics in Mesons Peter Tandy Center for - - PowerPoint PPT Presentation

soft and hard scale qcd dynamics in mesons
SMART_READER_LITE
LIVE PREVIEW

Soft and Hard Scale QCD Dynamics in Mesons Peter Tandy Center for - - PowerPoint PPT Presentation

Soft and Hard Scale QCD Dynamics in Mesons Peter Tandy Center for Nuclear Research Kent State University Mazatlan Nov09 p. 1/5 Topics Overview of DSE modeling of meson physicsmainly soft scale Masses, decays, form factors Including


slide-1
SLIDE 1

Soft and Hard Scale QCD Dynamics in Mesons

Peter Tandy Center for Nuclear Research Kent State University

Mazatlan Nov09 – p. 1/5

slide-2
SLIDE 2

Topics

Overview of DSE modeling of meson physics—mainly soft scale Masses, decays, form factors Including a hard scale: DIS: quark distributions in π, K mesons Mesons involving a heavy quark Summary

Mazatlan Nov09 – p. 2/5

slide-3
SLIDE 3

Lattice-QCD and DSE-based modeling

Lattice: O =

qqG O(¯ q, q, G) e−S[¯

q,q,G]

Euclidean metric, x-space, Monte-Carlo Issues: lattice spacing and vol, sea and valence mq, fermion Det Large time limit ⇒ nearest hadronic mass pole EOMs (DSEs): 0 =

qqG

δ δq(x) e−S[¯ q,q,G]+(¯ η,q)+(¯ q,η)+(J,G)

Euclidean metric, p-space, continuum integral eqns Issues: truncation and phenomenology—not full QCD Analtyic contin. ⇒ nearest hadronic mass pole Can be quick to identify systematics, mechanisms, · · ·

Mazatlan Nov09 – p. 3/5

slide-4
SLIDE 4

DSE-based modeling of Hadron Physics

Soft physics: truncate DSEs to min: 2-pt, 3-pt fns Should be relativistically covariant—-convenient for decays, Form Factors, etc No boosts needed on wavefns of recoiling bound st. ∞ d.o.f. → few quasi-particle effective d.o.f. Do not make a 3-dimensional reduction Preserve 1-loop QCD renorm group behavior in UV Preserve global symmetries, conserved em currents, etc Preserve PCAC ⇒ Goldstone’s Thm Can’t preserve local color gauge covariance—-just choose Landau gauge [RG fixed pt] Parameterize the deep infrared (large distance) QCD coupling

Mazatlan Nov09 – p. 4/5

slide-5
SLIDE 5

Constraints on Modeling

Preserve vector WTI, and axial vector WTI E.g. −iPµΓ5µ(k; P) = S−1(k+)γ5

τ 2 + γ5 τ 2S−1(k−)

−2 mq(µ) Γ5(k; P) ⇒ kernels of DSEq and KBSE are related Ladder-rainbow is the simplest implementation Goldstone Theorem preserved, ps octet masses good, indep of model details DCSB ⇒ π: Γ0

π(p2) = iγ5 f0

π

[ 1

4 tr S−1 0 (p2)] + · · ·

Here, 1-body and 2-body systems are the same

Mazatlan Nov09 – p. 5/5

slide-6
SLIDE 6

Ladder-Rainbow Model

K

  • KBSE → −γµ

λa 2

4παeff(q2) Dfree

µν (q) γν λa 2

αeff(q2)

→ IR ¯

qqµ=1 GeV = −(240MeV)3 , incl vertex dressing αeff(q2)

→ UV

α1−loop

s

(q2)

= +

p

  • 1
  • 1

p k p-k

P . Maris & P .C. Tandy, PRC60, 055214 (1999) Mρ, Mφ, MK⋆ good to 5%, fρ, fφ, fK⋆ good to 10%

Mazatlan Nov09 – p. 6/5

slide-7
SLIDE 7

Range of light meson observables

Summary of light meson results mu=d = 5.5 MeV, ms = 125 MeV at µ = 1 GeV

Pseudoscalar (PM, Roberts, PRC56, 3369) expt. calc.

  • ¯

qq0

µ

(0.236 GeV)3 (0.241†)3

0.1385 GeV 0.138†

0.0924 GeV 0.093†

mK

0.496 GeV 0.497†

fK

0.113 GeV 0.109 Charge radii (PM, Tandy, PRC62, 055204)

r2

π

0.44 fm2 0.45

r2

K+

0.34 fm2 0.38

r2

K0

  • 0.054 fm2
  • 0.086

γπγ transition (PM, Tandy, PRC65, 045211) gπγγ

0.50 0.50

r2

πγγ

0.42 fm2 0.41 Weak Kl3 decay (PM, Ji, PRD64, 014032)

λ+(e3)

0.028 0.027

Γ(Ke3)

7.6 ·106 s−1 7.38

Γ(Kµ3)

5.2 ·106 s−1 4.90 Vector mesons (PM, Tandy, PRC60, 055214)

mρ/ω

0.770 GeV 0.742

fρ/ω

0.216 GeV 0.207

mK⋆

0.892 GeV 0.936

fK⋆

0.225 GeV 0.241

1.020 GeV 1.072

0.236 GeV 0.259 Strong decay (Jarecke, PM, Tandy, PRC67, 035202)

gρππ

6.02 5.4

gφKK

4.64 4.3

gK⋆Kπ

4.60 4.1 Radiative decay (PM, nucl-th/0112022)

gρπγ/mρ

0.74 0.69

gωπγ/mω

2.31 2.07

(gK⋆Kγ/mK)+

0.83 0.99

(gK⋆Kγ/mK)0

1.28 1.19 Scattering length (PM, Cotanch, PRD66, 116010)

a0

0.220 0.170

a2

0.044 0.045

a1

1

0.038 0.036

bsampl

Mazatlan Nov09 – p. 7/5

slide-8
SLIDE 8

DSE kernel constrained from Lattice QCD

Mazatlan Nov09 – p. 8/5

slide-9
SLIDE 9

Lattice-assisted DSE Results

Evident vertex enhancement Curvature in low mq depn M IR(p2) 40% below linear Chiral Extrapolation

¯ qqqu−lat

µ=1 GeV = −(190 MeV)3

¯ qqqu−lat ≈ ¯ qqexpt/2

fπ 30% low

2 4 6 8 10 p

2 (GeV 2)

1.0 1.2 1.4 1.6 1.8 2.0

v(p

2)

0.000 0.025 0.050 0.075 0.100 0.125 m(ζ=19 GeV) (GeV) 0.100 0.200 0.300 0.400 M(p

2= 0.38 GeV 2) (GeV) Mazatlan Nov09 – p. 9/5

slide-10
SLIDE 10

Qu-lattice S(p), D(q) mapped to a DSE kernel

S(p) = Z(p) [i p + M(p)]−1

1 2 3 4

p [GeV]

0.1 0.2 0.3 0.4 0.5

M (p) [GeV]

Old data New ’improved action’ data mq = 0.168GeV mq = 0.030GeV mq = 0.225GeV mq = 0.055GeV mq = 0.110GeV mq = 0.0GeV

Mazatlan Nov09 – p. 10/5

slide-11
SLIDE 11

Quenched lattice ⇒ mq Depn of DSE Kernel

10

  • 2

10

  • 1

10 10

1

10

2

10

3

q

2 [GeV 2]

10

  • 3

10

  • 2

10

  • 1

10 10

1

10

2

10

3

10

4

4 π αeff(q

2)/q 2

DSE-LR (MT) V(q

2,m=0)*D(q 2)

chiral quark Bhagwat,Pichowsky,Roberts,Tandy, PRC68, 015203 (2003)

Mazatlan Nov09 – p. 11/5

slide-12
SLIDE 12

Quenched lattice ⇒ mq Depn of DSE Kernel

10

  • 2

10

  • 1

10 10

1

10

2

10

3

q

2 [GeV 2]

10

  • 3

10

  • 2

10

  • 1

10 10

1

10

2

10

3

10

4

4 π αeff(q

2)/q 2

DSE-LR (MT) V(q

2,m=0)*D(q 2)

V(q

2, mu)*D(q 2)

u-quark Bhagwat,Pichowsky,Roberts,Tandy, PRC68, 015203 (2003)

Mazatlan Nov09 – p. 11/5

slide-13
SLIDE 13

Quenched lattice ⇒ mq Depn of DSE Kernel

10

  • 2

10

  • 1

10 10

1

10

2

10

3

q

2 [GeV 2]

10

  • 3

10

  • 2

10

  • 1

10 10

1

10

2

10

3

10

4

4 π αeff(q

2)/q 2

DSE-LR (MT) V(q

2,m=0)*D(q 2)

V(q

2, mu)*D(q 2)

V(q

2, ms)*D(q 2)

s-quark Bhagwat,Pichowsky,Roberts,Tandy, PRC68, 015203 (2003)

Mazatlan Nov09 – p. 11/5

slide-14
SLIDE 14

Quenched lattice ⇒ mq Depn of DSE Kernel

10

  • 2

10

  • 1

10 10

1

10

2

10

3

q

2 [GeV 2]

10

  • 3

10

  • 2

10

  • 1

10 10

1

10

2

10

3

10

4

4 π αeff(q

2)/q 2

DSE-LR (MT) V(q

2,m=0)*D(q 2)

V(q

2, mu)*D(q 2)

V(q

2, ms)*D(q 2)

V(q

2, mc)*D(q 2)

c-quark Bhagwat,Pichowsky,Roberts,Tandy, PRC68, 015203 (2003)

Mazatlan Nov09 – p. 11/5

slide-15
SLIDE 15

Quenched lattice ⇒ mq Depn of DSE Kernel

10

  • 2

10

  • 1

10 10

1

10

2

10

3

q

2 [GeV 2]

10

  • 3

10

  • 2

10

  • 1

10 10

1

10

2

10

3

10

4

4 π αeff(q

2)/q 2

DSE-LR (MT) V(q

2,m=0)*D(q 2)

V(q

2, mu)*D(q 2)

V(q

2, ms)*D(q 2)

V(q

2, mc)*D(q 2)

V(q

2, mb)*D(q 2)

b-quark Bhagwat,Pichowsky,Roberts,Tandy, PRC68, 015203 (2003)

Mazatlan Nov09 – p. 11/5

slide-16
SLIDE 16

Quark Confinement—positivity violation

Confinement/positivity analysis (Osterwalder-Schrader axiom No. 3) Fourier transf σS(p4, p = 0) to Eucl time T

5 10 15 20 25 30 T (GeV

  • 1)

10

  • 6

10

  • 5

10

  • 4

10

  • 3

10

  • 2

10

  • 1

10 |∆S(T)|

solid = lattice prop, dashed = MT DSE, dotted = cc pole eg

Mazatlan Nov09 – p. 12/5

slide-17
SLIDE 17

DSE and Lattice results for MV and Mps

Mazatlan Nov09 – p. 13/5

slide-18
SLIDE 18

Pion electromagnetic form factor

Λµ = (P ′ + P)µ Fπ(Q2) = Nc

  • d4q

(2π)4 Tr ¯ Γπ S iΓµ S Γπ S

  • π

π γ

Mazatlan Nov09 – p. 14/5

slide-19
SLIDE 19

Pion F(Q2): Low Q2

Mazatlan Nov09 – p. 15/5

slide-20
SLIDE 20

Kaon F(Q2): Low Q2

Mazatlan Nov09 – p. 16/5

slide-21
SLIDE 21

Pion electromagnetic form factor

1 2 3 4 Q

2 [GeV 2]

0.1 0.2 0.3 0.4 0.5 Q

2 Fπ(Q 2) [GeV 2]

Our prediction VMD ρ pole CERN ’80s Cornell ’70s

PM and Tandy, PRC62,055204 (2000) [nucl-th/0005015]

Mazatlan Nov09 – p. 17/5

slide-22
SLIDE 22

Pion electromagnetic form factor

1 2 3 4 Q

2 [GeV 2]

0.1 0.2 0.3 0.4 0.5 Q

2 Fπ(Q 2) [GeV 2]

Our prediction VMD ρ pole CERN ’80s JLab, 2001

JLab data from Volmer et al, PRL86, 1713 (2001) [nucl-ex/0010009] PM and Tandy, PRC62,055204 (2000) [nucl-th/0005015]

Mazatlan Nov09 – p. 17/5

slide-23
SLIDE 23

Pion electromagnetic form factor

2 4 6 Q

2 [GeV 2]

0.1 0.2 0.3 0.4 0.5 Q

2 Fπ(Q 2) [GeV 2]

Our prediction VMD ρ pole CERN ’80s JLab, 2001 JLab at 12 GeV

  • pert. QCD

JLab, 2006b JLab, 2006a

PM and Tandy, PRC62,055204 (2000) [nucl-th/0005015] 2006a: V. Tadevosyan et al, [nucl-ex/0607007], 2006b: T. Horn et al, [nucl-ex/0607005]

Mazatlan Nov09 – p. 17/5

slide-24
SLIDE 24

1-loop chiral correction to rπ vs mπ

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6

mπ [GeV]

0.1 0.2 0.3 0.4 0.5

2 [fm 2]

Ladder-rainbow DSE Expt C / fπ

2

1--loop Ch PT Ch PT contact/core term 12 L9

r/fπ 2

P . Maris and PCT, in preparation

Mazatlan Nov09 – p. 18/5

slide-25
SLIDE 25

1-loop chiral correction to rπ vs mπ

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6

mπ [GeV]

0.1 0.2 0.3 0.4 0.5

2 [fm 2]

Ladder-rainbow DSE Expt C / fπ

2

1--loop Ch PT Ch PT contact/core term 12 L9

r/fπ 2

P . Maris and PCT, in preparation

Mazatlan Nov09 – p. 18/5

slide-26
SLIDE 26

γ⋆π0 → γ Transition Form Factor

Q P-Q/2 P+Q/2

Abelian axial anomaly + π pole in Γ5µ ⇒ G(0, 0) Chiral limit G(0, 0) = 1 2 ⇒ Γπγγ to 2%

0.0 1.0 2.0 3.0 Q

2 [GeV 2]

0.0 0.2 0.4 0.6 0.8 1.0 f(Q

2)/gexpt

CELLO CLEO all 8 covariants 5 covariants BL monopole

0.0 0.1 0.9 1.0

Mazatlan Nov09 – p. 19/5

slide-27
SLIDE 27

γ⋆πγ⋆ Asymptotic Limit

Lepage and Brodsky, PRD22, 2157 (1980): LC-QCD/OPE ⇒

10

  • 2

10

  • 1

10 10

1

10

2

10

3

Q

2 [GeV 2]

10

  • 4

10

  • 3

10

  • 2

10

  • 1

10 F(Q

2,Q 2)

DSE results VMD dipole bare vertices (4/3) π

2 fπ 2 / Q 2

Mazatlan Nov09 – p. 20/5

slide-28
SLIDE 28

LR: Successes, Problems, Resolutions

Successes: S-wave mesons, PS and V, light quarks and QQ, no spurious thresholds Exact PS mass formula, Goldstone Thm, ∆MHF from DCSB fEW , strong decays, radiative decays, form factors, Q2 < 5GeV 2 Problems: Axial vector (L > 0) mesons (a1, b1, · · · ) too light Physical diquarks, no physical V or PS qQ states Excited states are difficult Probable Resolution: Quark-gluon vertex: Γµ ⇒ Σq ⇒ KBSE Use analysis of spacelike correlators, 3-pt functions

Mazatlan Nov09 – p. 21/5

slide-29
SLIDE 29

From Gluon vertex to BSE Kernel

A symmetry-preserving procedure [Bender, Roberts, von Smekal, PLB380,

(1996), nucl-th/9602012; Munczek 1995] ; Axial vector and vector WTIs, and

Goldstone Thm preserved KBSE(x′, y′; x, y) = −

δ δS(x,y)Σ(x′, y′)

Vertex Γµ(p, q) = diagrams ⇒ KBSE = diagrams If Σ contains: KBSE contains: Independent of model parameters. Model does not fight chiral symmetry, use light vector mesons to fix parameters

Mazatlan Nov09 – p. 22/5

slide-30
SLIDE 30

Deep Inelastic Lepton Scattering

k k’ q X θ P

Bjorken limit: ν = q · P/M → ∞ ; − q2 = Q2 → ∞ 0 < x =

Q2 2P ·q < 1

W αβ =

~ Im

2 q

P

q q

P P

=

1 2π Disc T αβ(ν)

W αβ = −(gαβ − qαqβ

q2 ) F1 + P α

T (q) P β T (q)

P·q

F2 F1(x) = Σq

e2

q

2 fq(x) + · · ·

Mazatlan Nov09 – p. 23/5

slide-31
SLIDE 31

Deep Inelastic Lepton Scattering

Convenient basis in Bj lim: nν = M

2ω(1, −1;

0⊥) ; n2 = 0 = p2 ; p · n = 2 .; ω = M/2 (rest frame) , ω = ∞ (IMF) P µ = M

2 (nµ + pµ) ; qµ→ν nµ + Mx 2 (nµ − pµ) + O( 1 ν)

W αβ → (a ν +b) (F2 −2x F1)+(−gαβ +nα P β

M + P α M nβ) F1 +O( 1 ν)

{W αβ qβ}LO = 0 = W αβ nβ handbag diagram ⇒ W αβ

HB nβ = 0, (LO current consv)

Mazatlan Nov09 – p. 24/5

slide-32
SLIDE 32

Deep Inelastic Lepton Scattering

T µν(LO) = T µν

GHB =

P q q P ζ

µ

q+ = q·n = −Mx, |ξ−| ∼

1 Mx

q− = q · p = 2ν, |ξ+| ∼ 0 DIS is hard and fast—confinement is soft and slow ⇒ S(k + q) →

γ+ 2 (k+−P +x)+iǫ

W µν ∝ {T µν(ǫ) − T µν(−ǫ)} ⇒ Euclidean model elements can be continued EG, π+target : fq(x) = 1 4π Z dξ−eiq+ξ−π(P)|¯ q(ξ−)γ+q(0)|π(P)c = −f¯

q(−x)

fq(x) = 1 2 tr Z d4k (2π)4 δ(k+ − P +x) S(k)γ+S(k) T(k, P) General T(k, P) = ¯ uπ+ scattering amplitude: s-channel structure → ”spectator ¯ d” ⇒ fu(x), 0 < x < 1 u-channel structure → ”spectator uu ¯ d” ⇒ f¯

u(−x), 0 < x < 1

correct x support

Mazatlan Nov09 – p. 25/5

slide-33
SLIDE 33

Deep Inelastic Lepton Scattering

Quark number sum: NV

q =

R 1

0 dx {fq(x) − f¯ q(x)} = 1 2P + π(P)|J+(0)|π(P) = 1

DSE calculation: uπ(x), uK(x), sK(x) [T. Nguyen, PCT, (2009)] BSE q¯ q solutions for π, K DSE solns for dressed quark S(k) Constituent mass approx for spectator propagator Vertex approx via Ward Id Γ

+

k k -P P P k

Mazatlan Nov09 – p. 26/5

slide-34
SLIDE 34

DIS on pion: from DSE-BSE solutions

  • Valence quarks, handbag

diagram, γ+, Γ+

W I(k)

  • Data: J. S. Conway et al,

PRD39, 92 (1989) Ml¯

l = 4.05 GeV

  • Previous: Hecht, Roberts,

Schmidt, PRC63, 025213 (2001) Γπ(k, P) ≈ iγ5 B0(k2)/f 0

π + · · ·

S(p) fit to data

0.0 0.2 0.4 0.6 0.8 1.0

x

0.0 0.1 0.2 0.3 0.4

xuv(x;q=5.2 GeV)

DSE-BSA 5.2 GeV E615 πN Drell-Yan Holt et al. (09) 5.2 GeV Hecht et al. 5.2 GeV

Large x behavior: (1 − x)α , α = ?

Mazatlan Nov09 – p. 27/5

slide-35
SLIDE 35

DIS on pion: large x behavior?

0.70 0.75 0.80 0.85 0.90 0.95 1.00

x

0.00 0.05 0.10 0.15 0.20

xuπ(x)

Full DSE-BSA, q = 4.05 GeV E615 πN Drell-Yan 4GeV Hecht, q = 4.05 GeV fit to DSE (α = 2.41) fit to Hecht (α = 2.32)

Fit: a x (1 − x)α(x) BSE ampls: pQCD behavior sets in at a larger scale

Mazatlan Nov09 – p. 28/5

slide-36
SLIDE 36

DIS on pion: large x behavior?

0.75 0.80 0.85 0.90 0.95 1.00

x

2.00 2.10 2.20 2.30 2.40 2.50

α(x)

Hecht Full DSE-BSE

Global fits to (limited) DIS data produce α ∼ 1.5 Parton model (F-J), pQCD (Brodsky, Ezawa), DSEs, ⇒ α ∼ 2+ Constituent q models, NJL, duality, etc ⇒ α ∼ 1

Mazatlan Nov09 – p. 29/5

slide-37
SLIDE 37

Quark Distributions in π and K

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 x 0.1 0.2 0.3 0.4 xu(x) xuk(x) (kaon) xsk(x) (kaon) xuπ(x) (pion)

Evolved to q = 4.05 GeV

Environmental depn of u(x) in accordance with effective quark mass u(x), s(x) difference in K in accordance with effective quark mass

Mazatlan Nov09 – p. 30/5

slide-38
SLIDE 38

uK(x)/uπ(x) Ratio

0.0 0.2 0.4 0.6 0.8 1.0

x

0.0 0.2 0.4 0.6 0.8 1.0

uK(x)/uπ(x)

Data: (Drell-Yan, CERN-SPS) J. Badier it al., PLB 93, 354 (1980); Ml¯

l = 4 − 8 GeV

u has greater fraction of Pπ than it has of PK, in accord with effective quark mass

Mazatlan Nov09 – p. 31/5

slide-39
SLIDE 39

Axial anomaly and η − η′ states

Ch symm: ∂µ(z)jα

5µ(z) q(x)¯

q(y) involves 2 trf(Fα)Qt(z)q(x)¯ q(y) Matrix elements, amputated ⇒ AV-WTI PµΓα

5µ(k; P) = −2i MαβΓβ 5(k; P) −δα,0 ΓA(k; P)

+S−1(k+) iγ5Fα + iγ5FαS−1(k−) Residues at PS poles ⇒ PS mass formula for arbitrary mq, any flavor: m2

pf α p = 2 Mαβρβ p + δα,0 np

, np = 2 trf(F0) 0|Qt|p ρα

p(µ) = 0|¯

q γ5Fα q|p , p = any PS —–[Bhagwat, Chang, Liu, Roberts, PCT, PRC (76), 2007; arXiv:0708.1118]

Mazatlan Nov09 – p. 32/5

slide-40
SLIDE 40

A Schematic Model: Flavor mixing, η, η′

KN KA

[Bhagwat, Chang, Liu, Roberts, PCT, PRC (76), 2007; arXiv:0708.1118]

Structure: KN = LR vector gluon exch, KA = F(γ5, P /γ5) ⊗ (γ5, P /γ5)F , F = diag(1/Mf) (Munczek-Nemirovsky) t-channel δ4(k) for KN and KA 2 strength parameters: ρ0 ⇒ KN , mη′ ⇒ KA. Fix mu, md, ms · · · via vector mesons

Mazatlan Nov09 – p. 33/5

slide-41
SLIDE 41

π0 − η − η′ mixing: 3 flavors

mu − md causes π0 to be mixed in:

135 MeV : |π0 ∼ 0.72 ¯ uu − 0.69 ¯ dd − 0.013 ¯ ss 455 MeV : |η ∼ 0.53 ¯ uu + 0.57 ¯ dd − 0.63 ¯ ss 922 MeV : |η′ ∼ 0.44 ¯ uu + 0.45 ¯ dd + 0.78 ¯ ss

mu = md ⇒

455 MeV : |η ∼ 0.55 (¯ uu + ¯ dd) − 0.63 ¯ ss, θη = −15.4◦ 924 MeV : |η′ ∼ 0.45 (¯ uu + ¯ dd) + 0.78 ¯ ss, θη′ = −15.7◦ Chiral limit: m2

η′ = (0.852 GeV)2 ≡ 2trf(F0) 0|Qt|η′/f0 η′

cf Witten-Veneziano a-v ghost scenario ⇒ m2

η′ = h2 + m2 GB

It is worth extending to a realistic LR model for KN with separable KA: one obtains access to decay constants, residues, and details of the mass relations

Mazatlan Nov09 – p. 34/5

slide-42
SLIDE 42

Flavor Non-singlet PS Mass Relation

γ5γµ γ5 fH m2

H = 2 mq(µ) ρH(µ)

5 10 15 20

mq / mup/down

0.2 0.4 0.6

pion mass [GeV]

pseudoscalar meson, BSE solution Gell-Man-Oakes-Renner relation

mstrange

PM, Roberts, Tandy, PLB420, 267 (1998) [nucl-th/9707003]

Mazatlan Nov09 – p. 35/5

slide-43
SLIDE 43

Inaccuracy of GMOR

Mazatlan Nov09 – p. 36/5

slide-44
SLIDE 44

Quark mass functions from DSE solutions

10

−2

10

−1

10 10

1

10

2

10

3

10

4

p

2 [GeV 2]

10

−5

10

−4

10

−3

10

−2

10

−1

10 10

1

M(p

2) = B(p 2)/A(p 2) [GeV]

b−quark c−quark s−quark u/d−quark chiral limit

Mazatlan Nov09 – p. 37/5

slide-45
SLIDE 45

Constituent Mass Concept for c- and b-quarks

All GeV D(uc) D∗(uc) Ds(sc) D∗

s(sc)

expt M 1.86 2.01 1.97 2.11 calc M 1.85(FIT) 2.04 1.97 2.17 expt f 0.222 ? 0.294 ? calc f 0.154 0.160 0.197 0.180 All GeV B(ub) B∗(ub) Bs(sb) B∗

s(sb)

Bc(cb) B∗

c(cb)

expt M 5.28 5.33 5.37 5.41 6.29 ? calc M 5.27(FIT) 5.32 5.38 5.42 6.36 6.44 expt f 0.176 ? ? ? ? ? calc f 0.105 0.182 0.144 0.20 0.210 0.18 Fit ⇒ constituent masses: Mcons

c

= 2.0 GeV, Mcons

b

= 5.3 GeV Consistent with MDSE(p2 ∼ −M2) generated by mc = 1.2 ± 0.2, mb = 4.2 ± 0.2, [PDG, µ = 2 GeV] Does heavy quark dressing contribute anything? Too much in this DSE model—no mass shell !

Mazatlan Nov09 – p. 38/5

slide-46
SLIDE 46

Compare Quark Masses with PDG

Mazatlan Nov09 – p. 39/5

slide-47
SLIDE 47

Quarkonia

All GeV Mηc fηc MJ/ψ fJ/ψ expt 2.98 0.340 3.09 0.411 calc with Mcons

c

3.02 0.239 3.19 0.198 calc with ΣDSE

c

(p2) 3.04 0.387 3.24 0.415 All GeV Mηb fηb MΥ fΥ expt 9.4 ? ? 9.46 0.708 calc with Mcons

b

9.6 0.244 9.65 0.210 calc with ΣDSE

b

(p2) 9.59 0.692 9.66 0.682 QQ and qQ decay constants too low by 30-50% in constituent mass approximation Quarkonia decay constants much better for DSE dressed quarks (within 5% of expt.) IR sector (gluon k below ∼ 0.8 GeV) contribute little for bb or cc quarkonia in DSE, BSEs QQ states are more point-like than qq or qQ states

Mazatlan Nov09 – p. 40/5

slide-48
SLIDE 48

Recovery of a qQ Mass Shell

Suppress gluon k below ∼ 0.8 GeV in DSE dressing of b propagator Retain IR sector for dressed "light" quark and BSE kernel Now a mass shell is produced All GeV B(ub) B∗(ub) Bs(sb) B∗

s(sb)

Bc(cb) B∗

c(cb)

expt M 5.28 5.33 5.37 5.41 6.29 ? calc M 4.66 – 4.75 – 5.83 — expt f 0.176 ? ? ? ? ? calc f 0.133 – 0.164 – 0.453 – Masses are ∼ 10 % low It makes sense that Rb < RqQ ⇒ greater limit on low k in Σb May be partial confirmation of Brodsky and Shrock’s suggestion of universal maximum wavelength for quarks/gluons in hadrons [Phys. Lett. B666, (2008)]

Mazatlan Nov09 – p. 41/5

slide-49
SLIDE 49

IR Suppression of Kernel

10

  • 4

10

  • 3

10

  • 2

10

  • 1

10 10

1

k

2(GeV 2)

10

  • 4

10

  • 3

10

  • 2

10

  • 1

10 10

1

10

2

4π α(k

2)/k 2

IR+UV IR UV

1 2 3 4 5 6 7 8 9 10 11 12 13 MH (GeV) 0.2 0.4 0.6 0.8 1 D(MH)(GeV

2), p IR min(GeV)

D(MH) (GeV

2)

p

IR min(MH) (GeV) (Kernel only)

Mazatlan Nov09 – p. 42/5

slide-50
SLIDE 50

The V-A Current Correlator

ΠV

µν(x) = 0| T jµ(x)j† ν(0) |0 ,

isovector currents jµ = ¯ uγµd, j5

µ = ¯

uγ5γµd ΠV

µν(P) = (P 2δµν − PµPν) ΠV (P 2)

ΠA

µν(P) = (P 2δµν − PµPν) ΠA(P 2) + PµPν ΠL(P 2)

Πµν

V (P) = -

q+ q- q γµ

Γν

V(q,P)

P

Z1(µ, Λ) Λ

mq = 0 : ΠV − ΠA = 0 , to all orders in pQCD ΠV − ΠA probes the scale for onset of non-perturbative phenomena in QCD

Mazatlan Nov09 – p. 43/5

slide-51
SLIDE 51

The 4-quark Condensates

Operator product expansion ⇒ leading uv behavior ΠV −A(P 2) = 32παs¯

qq¯ qq 9 P 6

  • 1 + αs

4π[ 247 4π + ln( µ2 P 2)]

  • + O( 1

P 8)

Often vacuum saturation (¯ qq¯ qq ≈ ¯ qq2) is assumed for QCD Sum Rules. Validity not known. Extract ¯ qq¯ qq from lim|P 2→∞ P 6ΠV −A(P 2)

Model − < ¯ qq >µ=19 (GeV )3 < ¯ qq¯ qq >µ=19 (GeV )6 R(µ = 19) Set A (0.5682)3 (0.619)6 1.67 Set B (0.1734)3 (0.1902)6 1.74 Set C (0.2469)3 (0.2695)6 1.69 Set D (0.216)3 (0.235)6 1.65 —–T. Nguyen, PCT, in preparation, 2008

Mazatlan Nov09 – p. 44/5

slide-52
SLIDE 52

DSE Calculation: Weinberg Sum Rules

I:

1 4π2

0 ds[ρv(s) − ρa(s)] = [P 2 ΠV −A(P 2)]P 2→0 = −f 2 π

II: P 2 [P 2 ΠV −A(P 2)]|P 2→∞ = 0 DGMLY: ∞

0 dP 2 [P 2 ΠV −A(P 2)] = − 4πf2

π

3α [m2 π± − m2 π0] Model f2

π (GeV 2)

fπ (MeV ) fexp

π

/fnum

π

∆mπ (MeV ) (∆mπ)exp Set A 0.00456291 67.5 1.37 4.86 Set B 0.00538895 73.4 1.26 5.2 4.43 ± 0.03 Set C 0.00518379 72.0 1.28 4.88

Mazatlan Nov09 – p. 45/5

slide-53
SLIDE 53

Summary

Effective ladder-rainbow model based on QCD -DSEs; ¯ qqµ ⇒ 1 IR parameter Convenient and covariant approach to hadronic form factors: N, π, various transitions Ground state qQ and QQ mesons (V & PS) up to b-quark region Dynamical dressing in S(p) at each stage increases the value of the decay constant [factor of 3 for ¯ bb, factor of 2 for ¯ cc] ! First combination of BSE-DSE solutions for pion and kaon DIS distributions u(x), s(x) Used J J, V-A, to estimate ¯ qq¯ qq as ∼ 70% greater than vac saturation, and npQCD enters at scale 0.5 fm.

Mazatlan Nov09 – p. 46/5

slide-54
SLIDE 54

Collaborators

Craig Roberts, Argonne National Lab Pieter Maris, Iowa State University Yu-xin Liu, Lei Chang, Peking University Nick Souchlas, Trang Nguyen, Kent State University

Thankyou!

Mazatlan Nov09 – p. 47/5

slide-55
SLIDE 55

DIS on pion: from DSE-BSE solutions

  • Valence quarks, handbag

diagram, γ+, Γ+

W I(k)

  • Data: J. S. Conway et al,

PRD39, 92 (1989)

  • Previous: Hecht, Roberts,

Schmidt, PRC63, 025213 (2001) Γπ(k, P) ≈ iγ5 B0(k2)/f 0

π + · · ·

S(p) fit to data

0.0 0.2 0.4 0.6 0.8 1.0

x

0.1 0.2 0.3 0.4

xuπ(x)

reduced DSE-BSE (E

0, F 0, G 1)

Full DSE-BSE, q = 4.05 GeV E615 πN Drell-Yan 4GeV Hecht

Large x behavior: (1 − x)α , α = ?

Mazatlan Nov09 – p. 48/5

slide-56
SLIDE 56

Constituent Quark-like Behavior for c, b-quarks

Mass shell positions marked for ¯ bb and ¯ cc quarkonia qQ mesons sample MQ(p2) ∼ 4 times further into timelike region The same constituent or pole mass is unlikely to suffice for QQ and qQ mesons

Mazatlan Nov09 – p. 49/5

slide-57
SLIDE 57

General Pseudoscalar Mass Formula

Nf = 3, charge neutral states: p = π0, η, η′

m2

p

      f3

p

f8

p

f0

p

      =      np      +

  • 2 M3×3

     ρ3

p

ρ8

p

ρ0

p

     

Isospin breaking: mu = md allows anomaly, F 0, and s¯ s into π0 η′ in SU(Nf) limit: m2

η′f 0 η′ = nη′ + 2m ρ0 η′

Mazatlan Nov09 – p. 50/5

slide-58
SLIDE 58

Model Bethe-Salpeter Kernel for flavor singlet?

Vertex integral eqns do not involve Qt(x) explicitly:

Γα

5µ(k; P) = Z2 γ5γµFα +

Λ K S+Γα

5µS−

DSE models need: KBSE = KN + KA, both are ¯ qq irreducible, KN is also n-gluon irreducible

KA ∼ f1 f2

IS

IS e.g. IS = f1 f2

A scenario that works: Witten-Veneziano massless axial-vector ghost linking pseudoscalar GBs

Mazatlan Nov09 – p. 51/5

slide-59
SLIDE 59

c- and b-Quark Mass Function for BSE

  • 26 -24 -22 -20 -18 -16 -14 -12 -10
  • 8
  • 6
  • 4
  • 2

2 4 6 8 10 q

2(GeV 2)

1 2 3 4 5 6 Re(Mc,b(q

2))(GeV)

c,b quark mass function near the peak of the parabolic region with P

2 near the meson mass shells

mc(19 GeV)=0.88 GeV, mb(19 GeV)=3.8 GeV

Mazatlan Nov09 – p. 52/5

slide-60
SLIDE 60

DSE Calculation: Estimated 4 quark condensate

Model − < ¯ qq >µ=19 (GeV )3 < ¯ qq¯ qq >µ=19 (GeV )6 R(µ = 19) Set A (0.5682)3 (0.619)6 1.67 Set B (0.1734)3 (0.1902)6 1.74 Set C (0.2469)3 (0.2695)6 1.69 Set D (0.216)3 (0.235)6 1.65

Mazatlan Nov09 – p. 53/5