Goldstone bosons in the CFL phase

Transcription

1 Goldstone bosons in the CFL phase Verena Werth 1 Michael Buballa 1 Micaela Oertel 2 1 Institut für Kernphysik, Technische Universität Darmstadt 2 Observatoire de Paris-Meudon Dense Hadronic Matter and QCD Phase Transitions Workshop III Rathen, October 15, 2006

2 Outline 1 Motivation 2 Formalism Model Method 3 Results Equal quark masses Dependence on the strange quark mass Calculation of f π 4 Summary & Outlook

3 Phase diagram of neutral dense quark matter [Rüster, V.W., Buballa, Shovkovy, Rischke Phys. Rev. D 72 (2005)] first order phase transition second order phase transition normal quark matter: ud = us = ds = 0 2SC phase: ud 0, us = ds = 0 usc phase: ud, us 0, ds = 0 CFL phase: ud, us, ds 0

4 Phase diagram of neutral dense quark matter [Rüster, V.W., Buballa, Shovkovy, Rischke Phys. Rev. D 72 (2005)] first order phase transition second order phase transition normal quark matter: ud = us = ds = 0 2SC phase: ud 0, us = ds = 0 usc phase: ud, us 0, ds = 0 CFL phase: ud, us, ds 0

5 Why are we interested in Goldstone bosons? chiral symmetry breaking Goldstone bosons lightest excitations: mesons effective theories predict low meson masses and meson condensation

6 Why are we interested in Goldstone bosons? chiral symmetry breaking Goldstone bosons lightest excitations: mesons effective theories predict low meson masses and meson condensation influence of the pseudoscalar meson condensates studied by Warringa [hep-ph/ ]

7 CFL phase Diquark condensates diquark condensates: q T Ôq Ô antisymmetric scalar diquark condensates ud = 2 q T Cγ 5 τ 2 λ 2 q us = 5 q T Cγ 5 τ 5 λ 5 q ds = 7 q T Cγ 5 τ 7 λ 7 q τ A : flavor space λ A : color space

8 CFL phase Diquark condensates diquark condensates: q T Ôq Ô antisymmetric scalar diquark condensates ud = 2 q T Cγ 5 τ 2 λ 2 q us = 5 q T Cγ 5 τ 5 λ 5 q ds = 7 q T Cγ 5 τ 7 λ 7 q color-flavor locking Condensates invariant under q e iθa τa λt a 2 q τ A : flavor space λ A : color space

9 CFL phase Diquark condensates diquark condensates: q T Ôq Ô antisymmetric scalar diquark condensates ud = 2 q T Cγ 5 τ 2 λ 2 q us = 5 q T Cγ 5 τ 5 λ 5 q ds = 7 q T Cγ 5 τ 7 λ 7 q τ A : flavor space λ A : color space color-flavor locking Condensates invariant under q e iθa τa λt a 2 q chiral symmetry breaking A not invariant under q e iθa τa 2 γ 5 q Goldstone bosons

10 Nambu Jona-Lasinio model Lagrangian L eff = q(i / ˆm)q + H [( qiγ5 τ A λ A C q T ) ( q T Ciγ 5 τ A λ A q ) A=2,5,7 A =2,5,7 + ( qτ A λ A C q T ) ( q T Cτ A λ A q )]

11 Nambu Jona-Lasinio model Lagrangian L eff = q(i / ˆm)q + H [( qiγ5 τ A λ A C q T ) ( q T Ciγ 5 τ A λ A q ) A=2,5,7 A =2,5,7 + ( qτ A λ A C q T ) ( q T Cτ A λ A q )] work in Nambu-Gorkov space ψ = 1 ( ) q 1 2 C q T, ψ = 2 ( q, q T C )

12 Nambu Jona-Lasinio model Lagrangian L eff = q(i / ˆm)q + H [( qiγ5 τ A λ A C q T ) ( q T Ciγ 5 τ A λ A q ) A=2,5,7 A =2,5,7 + ( qτ A λ A C q T ) ( q T Cτ A λ A q )] four vertices in Nambu-Gorkov space ( ) Γ ll 0 0 s = iγ 5 τ A λ A 0 ( ) Γ ll 0 0 ps =, Γ ur ps = τ A λ A 0, Γ ur s = ( ) 0 iγ5 τ A λ A 0 0 ( ) 0 τa λ A 0 0

13 Ingredients inverse propagator in Nambu-Gorkov space p/ + ˆµγ 0 ˆm A γ 5 τ A λ A S 1 = A=2,5,7 Aγ 5 τ A λ A p/ ˆµγ 0 ˆm A=2,5,7 scalar diquark condensates ud = 2 q T Cγ 5 τ 2 λ 2 q us = 5 q T Cγ 5 τ 5 λ 5 q ds = 7 q T Cγ 5 τ 7 λ 7 q

14 Ingredients inverse propagator in Nambu-Gorkov space p/ + ˆµγ 0 ˆm A γ 5 τ A λ A S 1 = A=2,5,7 Aγ 5 τ A λ A p/ ˆµγ 0 ˆm A=2,5,7 scalar diquark condensates ud = 2 q T Cγ 5 τ 2 λ 2 q us = 5 q T Cγ 5 τ 5 λ 5 q ds = 7 q T Cγ 5 τ 7 λ 7 q A and µ 8 derived from self-consistent solution of the CFL gap equation + neutrality conditions

15 Bethe-Salpeter equation Bethe-Salpeter equation = = + i ˆT = i ˆV + i ˆV ( iĵ)i ˆT

16 Bethe-Salpeter equation Bethe-Salpeter equation = = + i ˆT = i ˆV + i ˆV ( iĵ)i ˆT ˆT = Γ i T ij Γ j ˆV = Γ i V ij Γ j Γ i ĴΓ j = J ij

17 Bethe-Salpeter equation Bethe-Salpeter equation = = + i ˆT = i ˆV + i ˆV ( iĵ)i ˆT T -matrix ˆT = Γ i T ij Γ j Γ i ĴΓ j = J ij ˆV = Γ i V ij Γ j T = V + V JT = T = ( 1 V J) 1 V

18 Polarization function Γ j p + q Γ k ij Γ j Γ k (q) = d 4 p (2π) Tr[Γ j is(p + q)γ kis(p)] p

19 Polarization function Γ j p + q Γ k ij Γ j Γ k (q) = d 4 p 1 (2π) 4 2 Tr[Γ j is(p + q)γ kis(p)] p applying Matsubara formalism and restriction to q = 0 q=0 ij Γ j Γ (q) = i k 2 T n d 3 p (2π) Tr[Γ j is(iω n + iω m, p)γ k is(iω n, p)]

20 Polarization function Γ j p + q Γ k ij Γ j Γ k (q) = d 4 p 1 (2π) 4 2 Tr[Γ j is(p + q)γ kis(p)] p applying Matsubara formalism and restriction to q = 0 q=0 ij Γ j Γ (q) = i k 2 T n d 3 p 1 (2π) 3 2 Tr[Γ j is(iω n + iω m, p)γ k is(iω n, p)] only a few vertex-combinations non-vanishing J can be transformed into block-diagonal structure

21 Structure of T choose the basis in which J is block-diagonal V is diagonal T is block-diagonal problem can be decomposed into smaller blocks: six 2x2 blocks (correspond to π +, π, K +, K, K 0, and K 0 ) one 6x6 block (corresponds to π 0, η, and η )

22 Mesons ˆT couples to external meson current Γ M

23 Mesons ˆT couples to external meson current = we have to calculate which diquark vertices contribute to the loop Γ M Γ j Γ M

24 Mesons ˆT couples to external meson current = we have to calculate which diquark vertices contribute to the loop Γ M Γ j Γ M meson vertices Γ π ± = i «τ1 ± iτ γ5 0 τ 1 iτ 2 Γ K ± = i «τ4 ± iτ γ5 0 τ 4 iτ 5 Γ ( ) = i «τ6 ± iτ 7 0 K 0 2 γ5 0 τ 6 iτ 7

25 Mesons ˆT couples to external meson current = we have to calculate which diquark vertices contribute to the loop Γ M Γ j Γ M meson vertices Γ π ± = i «τ1 ± iτ γ5 0 τ 1 iτ 2 Γ K ± = i «τ4 ± iτ γ5 0 τ 4 iτ 5 Γ ( ) = i «τ6 ± iτ 7 0 K 0 2 γ5 0 τ 6 iτ 7 non-vanishing combinations Γ π + Γ ll 75, Γ ur 57 Γ π Γ ll 57, Γ ur 75 Γ K + Γ ll 72, Γ ur 27 Γ K Γ ll 27, Γ ur 72 Γ K 0 Γ ll 52, Γ ur 25 Γ K0 Γ ll 25, Γ ur 52

26 Meson mass for equal quark masses π ±, K ±, K 0, K m meson [MeV] m q [MeV]

27 Meson mass for equal quark masses π ±, K ±, K 0, K 0 m meson [MeV] m q [MeV] prediction from EFT: m M = am q r 8A a = f 2 π A = 3 2 4π 2 fπ ln 2 µ 2 = 18 2π 2 [Son, Stephanov, Phys. Rev. D 61 (2000)]

28 Meson mass for equal quark masses π ±, K ±, K 0, K 0 m meson [MeV] calculated from EFT m q [MeV] MeV : ) a fit 0.38 a EFT 0.58 prediction from EFT: m M = am q r 8A a = f 2 π A = 3 2 4π 2 fπ ln 2 µ 2 = 18 2π 2 [Son, Stephanov, Phys. Rev. D 61 (2000)] fit 35 % smaller

29 Meson mass for equal quark masses π ±, K ±, K 0, K 0 m meson [MeV] calculated from EFT m q [MeV] MeV : ) a fit 0.29 a EFT 0.39 prediction from EFT: m M = am q r 8A a = f 2 π A = 3 2 4π 2 fπ ln 2 µ 2 = 18 2π 2 [Son, Stephanov, Phys. Rev. D 61 (2000)] fit 28 % smaller

30 Meson mass for equal quark masses π ±, K ±, K 0, K 0 m meson [MeV] calculated from EFT m q [MeV] MeV : ) a fit 0.21 a EFT 0.27 prediction from EFT: m M = am q r 8A a = f 2 π A = 3 2 4π 2 fπ ln 2 µ 2 = 18 2π 2 [Son, Stephanov, Phys. Rev. D 61 (2000)] fit 23 % smaller

31 Meson mass for equal quark masses π ±, K ±, K 0, K 0 m meson [MeV] calculated from EFT m q [MeV] MeV : ) a fit 0.14 a EFT 0.16 prediction from EFT: m M = am q r 8A a = f 2 π A = 3 2 4π 2 fπ ln 2 µ 2 = 18 2π 2 [Son, Stephanov, Phys. Rev. D 61 (2000)] fit 14 % smaller

32 Meson mass for equal quark masses π ±, K ±, K 0, K 0 m meson [MeV] calculated from EFT m q [MeV] 9.53 MeV : ) a fit a EFT prediction from EFT: m M = am q r 8A a = f 2 π A = 3 2 4π 2 fπ ln 2 µ 2 = 18 2π 2 [Son, Stephanov, Phys. Rev. D 61 (2000)] deviation < 1%

33 Dependence on the coupling comparison with EFT a calc / a EFT [MeV] for weaker interaction meson masses approach the EFT result

34 Dependence on the strange quark mass m u = m d = 30 MeV q 0 (pole) [MeV] π ± K -, K 0 K +, K m s [MeV]

35 Dependence on the strange quark mass m u = m d = 30 MeV prediction from EFT q 0 (pole) [MeV] π ± K -, K 0 K +, K m s [MeV] [Bedaque, Schaefer, Nucl.Phys. A697 (2002)] q meson(pole) = µ meson + m meson q π ± (pole) = m2 d m2 u 2µ s 4A + fπ 2 m s(m u + m d ) q K ± (pole) = m2 s m 2 u 2µ s 4A + fπ 2 m d (m u + m s) q K 0, K 0(pole) = m2 s m 2 d 2µ s 4A + fπ 2 m u(m d + m s)

36 Dependence on the strange quark mass m u = m d = 30 MeV prediction from EFT q 0 (pole) [MeV] π ± K -, K 0 K +, K m s [MeV] [Bedaque, Schaefer, Nucl.Phys. A697 (2002)] q meson(pole) = µ meson + m meson q π ± (pole) = m2 d m2 u 2µ s 4A + fπ 2 m s(m u + m d ) q K ± (pole) = m2 s m 2 u 2µ s 4A + fπ 2 m d (m u + m s) r 4A f 2 π a fit 2 q K 0, K 0(pole) = m2 s m 2 d 2µ s 4A + fπ 2 m u(m d + m s)

37 Kaon propagator m u = m d = 30 MeV inverse K propagator for m s = 30 MeV 1000 comparison with a free boson propagator with chemical potential µ K and mass m K t K - -1 [MeV 2 ] q 0 [MeV] t 1 K = (q0 + µ K )2 m 2 K g 2 µ K = m2 s m 2 u 2µ r afit m K = m d (m u + m s) 2

38 Kaon propagator m u = m d = 30 MeV inverse K propagator for m s = 30 MeV 1000 comparison with a free boson propagator with chemical potential µ K and mass m K t K - -1 [MeV 2 ] calculated fitted q 0 [MeV] t 1 K = (q0 + µ K )2 m 2 K g 2 µ K = m2 s m 2 u 2µ r afit m K = m d (m u + m s) 2

39 Kaon propagator m u = m d = 30 MeV inverse K propagator for m s = 60 MeV 1000 comparison with a free boson propagator with chemical potential µ K and mass m K t K - -1 [MeV 2 ] calculated fitted q 0 [MeV] t 1 K = (q0 + µ K )2 m 2 K g 2 µ K = m2 s m 2 u 2µ r afit m K = m d (m u + m s) 2

40 Kaon propagator m u = m d = 30 MeV inverse K propagator for m s = 120 MeV 1000 comparison with a free boson propagator with chemical potential µ K and mass m K t K - -1 [MeV 2 ] calculated fitted q 0 [MeV] t 1 K = (q0 + µ K )2 m 2 K g 2 µ K = m2 s m 2 u 2µ r afit m K = m d (m u + m s) 2

41 Meson propagators inverse π propagators inverse K propagators t π - -1 [MeV 2 ] m s =30 MeV m s =60 MeV m s =120 MeV q 0 t K - -1 [MeV 2 ] m s =30 MeV m s =60 MeV m s =120 MeV q 0

42 Pion decay constant A µ π A 0 π = i ( ) + 2 τ1 iτ 2 γ0 γ τ 1 + iτ 2 calculate loop at q 0 = q 0 (pole) d 4 p 1 [ ] f π q 0 = (2π) 4 2 Tr A 0 π S(p + q)(g + 1 Γ ll 75 + g 2 Γ ur 57)S(p)

43 Pion decay constant A µ π A 0 π = i ( ) + 2 τ1 iτ 2 γ0 γ τ 1 + iτ 2 calculate loop at q 0 = q 0 (pole) d 4 p 1 [ ] f π q 0 = (2π) 4 2 Tr A 0 π S(p + q)(g + 1 Γ ll 75 + g 2 Γ ur 57)S(p) for µ = 500 MeV, T = 0, and m u = m d = m s = 30 MeV we find f π 95.1 MeV comparison with effective theory: 21 8 ln 2 µ f π = MeV 36 π

44 Pion decay constant A µ π A 0 π = i ( ) + 2 τ1 iτ 2 γ0 γ τ 1 + iτ 2 calculate loop at q 0 = q 0 (pole) d 4 p 1 [ ] f π q 0 = (2π) 4 2 Tr A 0 π S(p + q)(g + 1 Γ ll 75 + g 2 Γ ur 57)S(p) for µ = 500 MeV, T = 0, and m u = m d = m s = 30 MeV we find f π 95.1 MeV 9 % smaller than in EFT comparison with effective theory: 21 8 ln 2 µ f π = MeV 36 π

45 Dependence on the coupling comparison with EFT f π [MeV] [MeV] for weaker interaction f π approaches the EFT result

46 Summary & Outlook Summary description of mesons with diquark loops good agreement with low-energy effective theory

47 Summary & Outlook Summary description of mesons with diquark loops good agreement with low-energy effective theory Outlook wider range of chemical potentials finite temperature include additional interactions quark-antiquark interactions mesons in the new CFL+meson groundstate