-
The change in chiral symmetry is one of the most important properties of quantum chromodynamics (QCD) in a hot and dense medium, which is essential for understanding the light hadrons at finite temperature and density [1-4]. In the chiral limit, the phase transition from chiral symmetry breaking in vacuum and at low temperature to its restoration at high temperature occurs at a critical temperature
Tc , which has been reported in a recent lattice QCD simulation to beTc≃132 MeV [5]. In a real case with non-vanishing current quark mass, the chiral symmetry restoration is no longer a genuine phase transition but rather a smooth crossover. Because the crossover occurs in a region and not at a point, the way to describe it with a fixed temperature is not unique. Considering the maximum fluctuations around a continuous phase transition in the chiral limit, the pseudo-critical temperatureTpc to characterize the chiral crossover in the real case is normally defined by the maximum change in the chiral condensate,∂2⟨ˉψψ⟩/∂T2pc=0 . From the lattice QCD simulation [6], this value is approximatelyTpc≃156 MeV.The mechanism for a continuous phase transition is spontaneous symmetry breaking. It is possible to define an order parameter that changes from a nonzero value to zero or vice versa when the phase transition occurs. Conversely, spontaneous breaking of a global symmetry manifests itself in Goldstone's theorem [7, 8]: whenever a global symmetry is spontaneously broken, massless fields, known as Goldstone bosons, emerge. Corresponding to the spontaneous chiral symmetry breaking, the order parameter is the chiral condensate, and the Goldstone modes are pions. If we take
Tpc defined above as the characteristic temperature of the chiral crossover, the problem is whether the chiral condensate atTpc is already small enough and the pseudo-Goldstone modes atTpc are already heavy enough to guarantee that the system is in chiral restoration phase. According to Goldstone's theorem, in the chiral breaking phase at low temperature, pions as pseudo-Goldstone modes should be in bound states, and in the chiral restoration phase at high temperature, pions should be in resonant states with nonzero width. The connection between the two states is the Mott transition at temperatureTm [9-11], where the decay processπ→qˉq begins. It is clear thatTpc=Tm in the chiral limit. In the real case, however, there is no guarantee of the coincidence of the two temperatures. In this case, Goldstone's theorem breaks down because pions may already be in resonant states with large mass atT<Tpc or still be in bound states with small mass atT>Tpc . To remain consistent with Goldstone's theorem, we propose to pin down the chiral crossover using the Mott transition of the pseudo-Goldstone boson. Taking into account energy conservation for the decay process, the Mott transition temperatureTm is defined through the pion massmπ(T) and quark massmq(T) ,mπ(Tm)=2mq(Tm).
(1) Because of the non-perturbative difficulty in QCD, we calculate the Mott transition temperature in an effective chiral model. One of the models that enables us to see directly how the dynamical mechanism of chiral symmetry breaking and restoration operates is the Nambu-Jona-Lasinio (NJL) model applied to quarks [12-16]. Within this model, the hadronic mass spectrum and the static properties of mesons can be obtained remarkably well. In this paper, we calculate in the model the chiral condensate in mean field approximation and the meson mass beyond mean field through random phase approximation (RPA), which has been proven to guarantee Goldstone's theorem in the chiral breaking phase.
In recent years, the investigation of chiral symmetry has been extended to include external electromagnetic fields. Considering the dimension reduction of fermions in external magnetic fields, chiral symmetry breaking in vacuum is enhanced by the background magnetic field, from both the lattice QCD simulation [17-23] and effective models [24-28]. The surprise is in the behavior of the chiral crossover temperature
Tpc . With increasing magnetic field, it decreases, from lattice QCD simulations [17-23], which is called the inverse magnetic catalysis, but increases from effective models in mean field approximation [24-28]. Many scenarios have been proposed to understand this qualitative difference between lattice QCD and effective models [29-42]. Note that, in previous works, the chiral crossover temperature is usually defined by the variation of the chiral condensate. In this work, we point out that, when defining the crossover temperature by the Mott transition of pseudo-Goldstone bosons, an inverse magnetic catalysis effect appears with decreasing Mott temperature in magnetic fields. Throughout this paper, the increase (decrease) in characteristic temperature with magnetic fields is called the (inverse) magnetic catalysis effect.The magnetized two-flavor NJL model is defined through the Lagrangian density [12-16]
L=ˉψ(iγμDμ−m0)ψ+G2[(ˉψψ)2+(ˉψiγ5τψ)2],
(2) where the covariant derivative
Dμ=∂μ+iQAμ couples quarks with electric chargeQ=diag(Qu,Qd)=diag(2e/3,−e/3) to the external magnetic fieldB=Bez through the potentialAμ=(0,0,Bx,0) , G is the coupling constant in scalar and pseudo-scalar channels, andm0 is the current quark mass characterizing the explicit chiral symmetry breaking.Using the Leung-Ritus-Wang method [43-48], the chiral condensate
⟨ˉψψ⟩ or the dynamical quark massmq=m0−G⟨ˉψψ⟩ at mean field level is controlled by the gap equationmq(1−GJ1)=m0
(3) with
J1=Nc∑f,nαn|QfB|2π∫dpz2πtanhEf2TEf,
(4) where
Nc=3 is the number of colors, which is trivial in the NJL model;αn=2−δn0 is the spin degeneracy; T is the temperature of the quark system; andEf=√p2z+2n|QfB|+m2q is the quark energy with flavorf=u,d , longitudinal momentumpz , and Landau energy level n.Mesons in the model are treated as quantum fluctuations above the mean field and constructed through RPA [12-16]. In the chiral limit with
m0=0 , the isospin tripletπ0 andπ± and isospin singletσ are respectively the Goldstone modes and Higgs mode, corresponding to spontaneous chiral symmetry breaking with a vanishing magnetic field. Turning on the external magnetic field, only the neutral pionπ0 remains as the Goldstone mode.With the RPA method, the meson propagator
Dm can be expressed in terms of the meson polarization function or quark bubbleΠm ,Dm(q)=G1−GΠm(q).
(5) The meson mass
mm is defined as the pole of the propagator at zero momentumq=0 ,1−GΠm(mm,0)=0
(6) with
Πm(q0,0)=J1−(q20−ϵ2m)J2(q0),J2(q0)=−Nc∑f,nαn|QfB|2π∫dpz2πtanhEf2TEf(4E2f−q20),
(7) ϵπ0=0 for the Goldstone mode, andϵσ=2mq for the Higgs mode. In a nonzero magnetic field, the three-dimensional quark momentum integration in the gap equation (3) and pole equation (6) becomes a one-dimensional momentum integration plus a summation over the discrete Landau levels.In the chiral limit with vanishing current quark mass, by comparing the gap equation (3) for quark mass with the pole equation (6) for meson mass, we have the analytic solutions
mπ0=0, mσ=2mq
(8) in the chiral breaking phase with
mq≠0 andmπ0=mσ≠0
(9) in the chiral restoration phase with
mq=0 . A direct consequence of these solutions is that the Mott transition temperatureTm defined bymπ0(Tm)=2mq(Tm) coincides with the critical temperatureTc defined bymq(Tc)=0 . The phase transition from chiral symmetry breaking to its restoration is a second order phase transition.In the physical case with nonzero current quark mass, the chiral restoration becomes a smooth crossover. At low temperature, spontaneous chiral symmetry breaking dominates the system. Considering the fact that the explicit chiral symmetry breaking is slight, we can use
m0 expansion to solve the gap equation for quark mass and pole equation for meson mass. With the notationsmq=mclq+δq andmπ0=mclπ0+δπ0 , wheremclq andmclπ0=0 are the quark and neutral pion masses in the chiral limit, respectively, and keeping only the linear term inδq and quadratic term inδπ0 in the gap and pole equations, we haveδq=−m0mclq1G∂J1∂mq|cl,δ2π0=−m0mclq+δq1GJ2|cl.
(10) It is obvious that, in the chiral limit with
m0=0 , we haveδq=0 andδπ0=0 . The explicit chiral symmetry breaking withm0≠0 modifies the dynamical quark mass, and the Goldstone mode in the chiral limit becomes a pseudo-Goldstone mode with nonzero mass.At high temperature, the quark dimension reduction under an external magnetic field causes an infrared (
pz→0 ) singularity of the quark bubbleΠm(mm,0) [49-51]. For the pseudo-Goldstone modeπ0 , the infrared singularity ofΠπ0(mπ0,0) occurs at the Mott transition temperatureTm , where the massmπ0 jumps up frommπ0<2mq tomπ0>2mq . This indicates a sudden transition from a bound state to a resonant state [49, 51].We next perform numerical calculations on the Mott transition temperature in both the chiral limit and real case. Because of the four-fermion interaction, the NJL model is not a renormalizable theory and requires regularization. To guarantee the law of causality in magnetic fields, we apply the Pauli-Villars regularization scheme, as explained in detail in Ref. [42]. The three parameters in the NJL model, namely the current quark mass
m0 , coupling constant G, and Pauli-Villars mass parameterΛ , are listed in Table 1 by fitting the chiral condensate⟨ˉψψ⟩ , pion massmπ , and pion decay constantfπ in a vacuum atT=B=0 . We take the current quark mass to bem0=0 in the chiral limit and6.4 MeV in the real world.m0 /MeVG/GeV −2 Λ /MeV⟨ˉψψ⟩ /MeV3mπ /MeVfπ /MeV0 5.03 977.3 −2303 0 93 6.4 4.9 977.3 −2303 134 93 Table 1. NJL parameters in Pauli-Villars regularization.
We first discuss the temperature behavior of the phase structure for chiral symmetry at
B=0 . Fig. 1 shows the quark mass and pion mass as functions of temperature in the chiral limit and the real world. For a vanishing magnetic field, the three pions are all Goldstone or pseudo-Goldstone modes. To clearly see the Mott transition temperatureTm and its difference from the pseudo-critical temperatureTpc , we plot2mq instead ofmq itself. In the chiral limit, the quark mass, which is proportional to the order parameter⟨ˉψψ⟩ , continuously decreases at low temperature, reaches zero at the critical temperatureTc=163 MeV, and remains zero at higher temperature. This denotes a second order chiral phase transition. Correspondingly, the Goldstone modesπ remain massless in the chiral breaking phase and begin to have mass atTc . It is clear that the Mott transition temperature defined by the threshold condition (1) is exactly the critical temperatureTm=Tc=163 MeV. Note that the critical temperature obtained here in the NJL model is higher by20−30 MeV than the results of recent lattice QCD [5] and other effective models [52, 53]. In the real world, the quark mass continuously decreases, and the pion mass continuously increases throughout the entire temperature region. In this case, the chiral phase transition becomes a smooth crossover, and there is no strict definition for the crossover temperature. Generalizing the idea of maximum fluctuations around the second order phase transition in the chiral limit, the maximum change in the chiral condensate or dynamical quark mass is commonly used to identify the crossover, referring to the corresponding temperature as the pseudo-critical temperatureTpc of the chiral crossover. From the definition∂2mq/∂T2pc=0 , we numerically haveTpc=162 MeV, which is close to the lattice QCD result of156 MeV [6]. From the definitionmπ(Tm)=2mq(Tm) , denoted by the crossing point of the two dashed lines in Fig. 1, the Mott transition temperature is different from the pseudo-critical temperature:Tm=174 MeV>Tpc=162 MeV. Therefore, pions as pseudo-Goldstone modes can still survive as bound states after the chiral crossover. This means that the definition ofTpc explicitly violates Goldstone's theorem.Figure 1. (color online) Dynamical quark mass
mq and pion massmπ as functions of temperature in the chiral limit (solid lines) with current quark massm0=0 and real world (dashed lines) withm0=6.4 MeV.In the chiral limit, chiral restoration is a genuine phase transition, and the phase transition temperature is unique from either the order parameter or Goldstone's theorem, with
Tc=Tm . In the physical world, which is the focus of this paper, chiral restoration is a smooth crossover. The characteristic temperatureTpc from the maximum change in order parameter andTm from the Mott transition of pseudo-Goldstone bosons will not coincide with each other:Tpc≠Tm . To guarantee Goldstone's theorem, we propose to define the crossover temperature asTm .Because charged pions interact with the magnetic field, they are no longer pseudo-Goldstone modes, and only the neutral pion is the pseudo-Goldstone boson corresponding to spontaneous chiral symmetry breaking. Fig. 2 shows the pseudo-critical temperature
Tpc and Mott transition temperatureTm for the pseudo-Goldstone mode as functions of the magnetic field. WhileTpc is controlled by the magnetic catalysis, i.e., it increases with increasing magnetic field, the Mott transition temperatureTm clearly shows the inverse magnetic catalysis effect, decreasing throughout the entire magnetic field region. The physics to explain this difference is as follows:Tpc is controlled by quarks, which are calculated at mean field level, butTm is governed by mesons, which are treated as quantum fluctuations beyond mean field. It is the quantum fluctuation that changes the magnetic catalysis to inverse magnetic catalysis. This is consistent with the scenario of fluctuation induced inverse magnetic catalysis discussed in Refs. [30, 40-42]. When we consider the feedback effect from mesons to quarks by including the meson contribution to the thermodynamics of the system,Ω=Ωmf+ΩM , the suppressed chiral condensate at high temperature and decreasing pseudo-critical temperatureTpc have been observed.Figure 2. (color online) Pseudo-critical temperature
Tpc (solid line) and Mott transition temperatureTm (dashed line) for the pseudo-Goldstone mode as functions of magnetic field in the physical case withm0=6.4 MeV.Again, the result of
Tpc≠Tm violates Goldstone's theorem in any magnetic field. In a weak magnetic field, witheB/m2π<7 , wheremπ is the pion mass in a vacuum atT=B=0 ,Tm>Tpc , which leads to the survival of the neutral pion as a bound state in the chiral restoration phase. Conversely, in a strong magnetic field witheB/m2π>7 ,Tm<Tpc , which results in the disappearance of the neutral pion in the chiral breaking phase. In either case, Goldstone's theorem is significantly broken.We investigate in this paper the chiral crossover at finite temperature and in external magnetic fields. In the physical world, chiral restoration is a smooth crossover because of the explicit chiral symmetry breaking. Different from the commonly used maximum change in chiral condensate, we propose defining the crossover temperature by the Mott transition of pseudo-Goldstone bosons. This, by definition, guarantees Goldstone's theorem for chiral symmetry. As an analytical example, we calculate the order parameter (dynamical quark mass) in mean field and Goldstone mode (pion mass) beyond mean field in the frame of a Pauli-Villars regularized NJL model. If we take the maximum change in chiral condensate to describe the chiral crossover, the pseudo-Goldstone mode will survive in the chiral restoration phase with weak magnetic fields and disappear in the chiral breaking phase with strong magnetic fields.
While we believe that the idea of choosing the Mott transition temperature
Tm to characterize the chiral crossover is correct, the values of the temperaturesTm,Tpc , andTc we obtained here are model dependent. To precisely fix the Mott transition temperature, we need to use the direct result from lattice QCD simulations.
HTML
-
I am grateful for the hospitality of Professor Dirk H. Rischke of Frankfurt University. Part of this work was done while I was visiting his group as an EMMI visiting professor.