Mesonic and nucleon fluctuation effects at finite baryon density Research Center for Nuclear Physics Osaka University Workshop on Strangeness and charm in hadrons and dense matter Yukawa Institute for Theoretical Physics, Kyoto University 22 May, 2017 GF & A. Hosaka, Phys. Rev. D 94, 036005 (2016) aaa GF & A. Hosaka, arxiv:1701.03717
Outline aaa Motivation Functional Renormalization Group Chiral effective nucleon-meson theory at finite µ B Numerical results Summary
Motivation FLUCTUATION EFFECTS IN FIELD THEORY: 2nd order transitions in statistical field theory diverging correlation length invalidates PT solution: Wilson s momentum space RG explanation of universality, critical exponents, etc.
Motivation FLUCTUATION EFFECTS IN FIELD THEORY: 2nd order transitions in statistical field theory diverging correlation length invalidates PT solution: Wilson s momentum space RG explanation of universality, critical exponents, etc. Large logarithms in continuum quantum field theory PT is organized in terms of λ log(e/m 2 ) PT can fail even for small couplings resummation! solution: Gell-Mann - Low RG
Motivation FLUCTUATION EFFECTS IN FIELD THEORY: 2nd order transitions in statistical field theory diverging correlation length invalidates PT solution: Wilson s momentum space RG explanation of universality, critical exponents, etc. Large logarithms in continuum quantum field theory PT is organized in terms of λ log(e/m 2 ) PT can fail even for small couplings resummation! solution: Gell-Mann - Low RG Problem in QCD: running coupling grows as E decreases RG cannot provide a generic solution effective models can help, but are also strongly coupled
Motivation FLUCTUATION EFFECTS IN FIELD THEORY: 2nd order transitions in statistical field theory diverging correlation length invalidates PT solution: Wilson s momentum space RG explanation of universality, critical exponents, etc. Large logarithms in continuum quantum field theory PT is organized in terms of λ log(e/m 2 ) PT can fail even for small couplings resummation! solution: Gell-Mann - Low RG Problem in QCD: running coupling grows as E decreases RG cannot provide a generic solution effective models can help, but are also strongly coupled If one is interested in the finite T and density behavior of the aa strongly interacting matter, fluctuation effects are important.
Motivation AXIAL ANOMALY OF QCD: U A (1) anomaly: anomalous breaking of the U A (1) subgroup of chiral symmetry vacuum-to-vacuum topological fluctuations (instantons) µ j µa A = g 2 16π 2 ɛµνρσ Tr [T a F µν F ρσ ] Induced U A (1) breaking interactions depend on instanton density suppressed at high T calculations are trustworthy only beyond T C is the anomaly present around and below T C? Recent lattice QCD simulations do not seem to have agreement on the issue 1,2 1 S. Sharma et al., Nucl. Phys. A 956, 793 (2016) 2 A. Tomiya et al., arxiv:1612.01908
Motivation η - NUCLEON BOUND STATE: Effective models at finite T and/or density: effective models (NJL 3, linear sigma models 4 ) predict a aaaadrop in m η at finite T and µ B Effective description of the mass drop: attractive potential in medium η N bound state Analogous to Λ(1405) KN bound state 3 P. Costa, M. C. Ruivo & Yu. L. Kalinovsky, Phys. Lett. B 560, 171 (2003). 4 S. Sakai & D. Jido, Phys. Rev. C88, 064906 (2013).
Motivation η - NUCLEON BOUND STATE: Effective models at finite T and/or density: effective models (NJL 3, linear sigma models 4 ) predict a aaaadrop in m η at finite T and µ B Effective description of the mass drop: attractive potential in medium η N bound state Analogous to Λ(1405) KN bound state Problem with effective model calculations: they treat model parameters as environment independent constants a v type of terms decrease (a-constant, v-decreases) evolution of a at finite T and µ B? 3 P. Costa, M. C. Ruivo & Yu. L. Kalinovsky, Phys. Lett. B 560, 171 (2003). 4 S. Sakai & D. Jido, Phys. Rev. C88, 064906 (2013).
Motivation CHIRAL EFFECTIVE NUCLEON-MESON MODEL: Chiral symmetry of QCD: U(N f ) U(N f ) (Ψ L/R U L/R Ψ L/R ) Effective model of mesons (M) and the nucleon (N): [M: π, K, η, η and a 0, κ, f 0, σ, N: n, p] L = Tr [ µ M µ M ] µ 2 Tr (MM ) g 1 9 [ Tr (MM )] 2 g 2 3 Tr (MM ) 2 Tr [H(M + M )] a(det M + det M ) + N( i γ i + µ B γ 0 m N )N g N M 5 N
Motivation CHIRAL EFFECTIVE NUCLEON-MESON MODEL: Chiral symmetry of QCD: U(N f ) U(N f ) (Ψ L/R U L/R Ψ L/R ) Effective model of mesons (M) and the nucleon (N): [M: π, K, η, η and a 0, κ, f 0, σ, N: n, p] L = Tr [ µ M µ M ] µ 2 Tr (MM ) g 1 9 [ Tr (MM )] 2 g 2 3 Tr (MM ) 2 Tr [H(M + M )] a(det M + det M ) + N( i γ i + µ B γ 0 m N )N g N M 5 N Model parameters: meson mass parameter (µ 2 ), quartic couplings (g 1, g 2 ), explicit breaking (H = h 0 T 0 + h 8 T 8 ), U A (1) anomaly (a), nucleon mass parameter (m N ), Yukawa coupling (g) Short-range N N interactions: ω and ρ mesons
Functional Renormalization Group Fluctuation effects are included in the partition function and/or in the quantum effective action Z[J] = Dφe (S[φ]+ Jφ) Γ[ φ] = log Z[J] J φ
Functional Renormalization Group Fluctuation effects are included in the partition function and/or in the quantum effective action Z[J] = Dφe (S[φ]+ Jφ) Γ[ φ] = log Z[J] J φ Functional RG approach: generalized Wilsonian RG Similarities: provides a way to handle IR singularities gradually integrates out high momentum modes Differences: flow of the complete effective action is considered IR regulator is not fixed optimization!
Functional Renormalization Group Fluctuation effects are included in the partition function and/or in the quantum effective action Z[J] = Dφe (S[φ]+ Jφ) Γ[ φ] = log Z[J] J φ Functional RG approach: generalized Wilsonian RG Similarities: provides a way to handle IR singularities gradually integrates out high momentum modes Differences: flow of the complete effective action is considered IR regulator is not fixed optimization! Functional Renormalization Group provides a non-perturbative approach to resummation of fluctuation effects.
Functional Renormalization Group Mathematical implementation: Scale dependent partition function: aa Z k [J] = Dφe (S[φ]+ Jφ) aaaaaaaa e 1 2 φrk φ R k (q) k 2 Scale dependent effective action: 0 0 k q Γ k [ φ] = log Z k [J] J φ 1 2 φr k φ k Λ: no fluctuations included aaaa Γ k [ φ] k=λ = S[ φ] k = 0: all fluctuations included aaaa Γ k [ φ] k=0 = Γ[ φ]
Chiral effective nucleon-meson model Flow equation of the effective action: k Γ k = 1 2 (T ) q,p k R k (q, p)(γ (2) k + R k ) 1 (p, q) = 1 2 One-loop structure with dressed and regularized propagators RG change in the n-point vertices are aaaadescribed solely by one-loop diagrams exact relation, approximations are necessary
Chiral effective nucleon-meson model Flow equation of the effective action: k Γ k = 1 2 (T ) q,p k R k (q, p)(γ (2) k + R k ) 1 (p, q) = 1 2 One-loop structure with dressed and regularized propagators RG change in the n-point vertices are aaaadescribed solely by one-loop diagrams exact relation, approximations are necessary Derivative expansion (local potential approximation): [ Γ k = Tr [ µ M µ M ] + N( ] i γ i + µ B γ 0 m N )N V k x V k = µ 2 k [M] Tr (MM ) + g 1,k[M] [ Tr (MM )] 2 + g 2,k[M] Tr (MM ) 2 9 3 + Tr [H k (M + M )] + A k [M](det M + det M ) + g k [M] N M 5 N
Chiral effective nucleon-meson model Flow equation: k Γ k = 1 2 (T ) q,p k R k (q, p)(γ (2) k + R k ) 1 (p, q) projection to various operators: individual flow equations aaaafor field dependent couplings
Chiral effective nucleon-meson model Flow equation: k Γ k = 1 2 (T ) q,p k R k (q, p)(γ (2) k + R k ) 1 (p, q) projection to various operators: individual flow equations aaaafor field dependent couplings Step I.: solve equations at T = 0 determine parameters µ, g 1, g 2, a: determined by fitting masses of π, K, η, η g: determined by fitting nucleon mass with m N minimal PCAC relations: (α = π, K) m 2 αf αˆπ α = µ J 5µ α = θa α Tr (H(M + M )) h 0 (286MeV ) 3 h 8 (311MeV ) 3
Chiral effective nucleon-meson model Flow equation: k Γ k = 1 2 (T ) q,p k R k (q, p)(γ (2) k + R k ) 1 (p, q) projection to various operators: individual flow equations aaaafor field dependent couplings Step I.: solve equations at T = 0 determine parameters µ, g 1, g 2, a: determined by fitting masses of π, K, η, η g: determined by fitting nucleon mass with m N minimal PCAC relations: (α = π, K) m 2 αf αˆπ α = µ J 5µ α = θa α Tr (H(M + M )) h 0 (286MeV ) 3 h 8 (311MeV ) 3 Step II.: solve the same equations at finite T and µ B spectrum, symmetry restoration, U A (1) anomaly
Numerical results: mass spectrum at finite T 1100 1000 900 800 masses [MeV] 700 600 500 400 f 0κ 300 a 0 η 200 η K 100 σ π N 0 0 100 200 300 400 500 600 700 T [MeV]
Numerical results: η η system at finite T 1000 900 η' mass η, η' masses [MeV] 800 700 600 η mass 500 400 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 T/T C solid: full solution, dashed: field- and T independent anomaly
Numerical results: anomaly at finite T 9 8 A k 0 [GeV] 7 6 5 T = 1.2 T c T = 0.6 T c T = 0 4 3 0 100 200 300 I 1 [MeV] I 1 = (v 2 ns + v 2 s )/2
Numerical results: anomaly at finite T
Numerical results: anomaly at finite µ B A (µ B ;T) [MeV] 1400 1200 1000 800 600 400 200 0 T=0 MeV T=30 MeV T=50 MeV T=100 MeV T=200 MeV 600 800 1000 1200 µ B [MeV]
Numerical results: heavy spectrum at finite µ B T = 0 MeV 1100 f 0 κ masses [MeV] 1000 900 N η a 0 800 400 500 600 700 800 900 1000 1100 1200 µ B [MeV]
Numerical results: heavy spectrum at finite µ B T = 100 MeV 1100 κ f 0 masses [MeV] 1000 900 N η a 0 800 400 500 600 700 800 900 1000 1100 1200 µ B [MeV]
Numerical results Mesonic and nucleon fluctuations seems to have a strong effect on the anomaly evolution they cause a relative change of about 20% aaaaat T T C and at µ B Λ ( 1 GeV )
Numerical results Mesonic and nucleon fluctuations seems to have a strong effect on the anomaly evolution they cause a relative change of about 20% aaaaat T T C and at µ B Λ ( 1 GeV ) Caution 1.) due to large anomaly the nucleon mass piece arises from aaaachiral symmetry breaking has an upper limit aaaa nucleon mass parameter m N becomes large and aaaaaa violates chiral symmetry
Numerical results Mesonic and nucleon fluctuations seems to have a strong effect on the anomaly evolution they cause a relative change of about 20% aaaaat T T C and at µ B Λ ( 1 GeV ) Caution 1.) due to large anomaly the nucleon mass piece arises from aaaachiral symmetry breaking has an upper limit aaaa nucleon mass parameter m N becomes large and aaaaaa violates chiral symmetry Caution 2.) large anomaly ruins the value of T C at µ B = 0 it is off by almost a factor of 2!
Numerical results 1 0.8 strange v s/ns (T)/v s/ns (0) 0.6 0.4 non strange 0.2 0 0 100 200 300 400 500 T [MeV] dashed: without anomaly, solid: with anomaly
Numerical results Mesonic and nucleon fluctuations seems to have a strong effect on the anomaly evolution they cause a relative change of about 20% aaaaat T T C and at µ B Λ ( 1 GeV ) Caution 1.) due to large anomaly the nucleon mass piece arises from aaaachiral symmetry breaking has an upper limit aaaa nucleon mass parameter m N becomes large and aaaaaa violates chiral symmetry Caution 2.) large anomaly ruins the value of T C at µ B = 0 it is off by almost a factor of 2! Way out: the bare anomaly parameter a has to depend explicitly on T and µ B. It should represent the underlying instanton dynamics of QCD.
Summary Finite T and density properties of the axial anomaly and mesonic spectra in a chiral effective nucleon-meson model Fluctuations (quantum and thermal) via the Functional Renormalization Group (FRG) approach thermal evolution of the mass spectrum and condensates temperature dependence of the U A (1) anomaly factor Findings: meson fluctuations strengthen the anomaly with respect aaaato the temperature no recovery at T C putting the system into nuclear medium also aaaaaffects the anomaly it increases as µ B grows η mass is increasing with µ B and T Important: T C of chiral restoration and m N comes out high! T -dependence of the bare anomaly coeff. can be relevant!