What could be learned about phase transitions, meson and nucleon structure in hot medium from a chiral quark-meson theory?
Christo V. Christov

TL;DR
This paper investigates phase transitions and the structure of mesons and nucleons in hot matter using a chiral quark-meson theory, revealing different transition types and thermodynamic behaviors at varying densities and temperatures.
Contribution
It introduces a hybrid NJL model combining Dirac and Fermi seas to analyze phase transitions and particle structures in hot medium, highlighting the nature of these transitions.
Findings
Quark matter exhibits a smooth chiral phase transition at high temperatures.
Nucleon matter shows a sharp, first-order phase transition with large discontinuities.
Nucleon solitons become less bound and eventually disappear at critical temperatures.
Abstract
In this contribution I summarize and discuss the results of the bulk thermodynamic characteristics, meson and nucleon structure in hot matter obtained in the framework of a chiral quark-meson theory. A hybrid NJL model is used in which a Dirac sea of quarks is combined with a Fermi sea of quarks or of nucleons. In the model mesons are described as collective excitations and the nucleon appears as a baryon-number-one soliton of valence quarks coupled to both Dirac and Fermi sea. According to the model at some critical density and/or temperature phase transitions from nucleons to quarks as well as from Goldstone to Wigner phase are expected. At finite density the chiral order parameter and the constituent quark mass have a non-monotonic temperature dependence-at temperatures not close to the critical one they are less affected than in cold matter. The quark matter is…
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4
Figure 5
Figure 6
Figure 7
Figure 8
Figure 9
Figure 10
Figure 11Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsHigh-pressure geophysics and materials · Quantum, superfluid, helium dynamics · High-Energy Particle Collisions Research
What could be learned about phase transitions, meson and nucleon structure in hot medium from a chiral quark-meson theory?
Christo V. Christov
1Institüt für Threoretische Physik II, Ruhr-Universität Bochum, Germany
2Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, Sofia 1784, Bulgaria
\runningheads
Christo V. ChristovWhat could be learned about phase transitions, meson and nucleon structure in medium…
{start}
1,2
1
2
{Abstract}
In this contribution I summarize and discuss the results of the bulk thermodynamic characteristics, meson and nucleon structure in hot matter obtained in the framework of a chiral quark-meson theory. A hybrid NJL model is used in which a Dirac sea of quarks is combined with a Fermi sea of quarks or of nucleons. In the model mesons are described as collective excitations and the nucleon appears as a baryon-number-one soliton of valence quarks coupled to both Dirac and Fermi sea. According to the model at some critical density and/or temperature phase transitions from nucleons to quarks as well as from Goldstone to Wigner phase are expected. At finite density the chiral order parameter and the constituent quark mass have a non-monotonic temperature dependence - at temperatures not close to the critical one they are less affected than in cold matter. The quark matter is rather soft against thermal fluctuations and the corresponding chiral phase transition is smooth. The nucleon matter is much stiffer and the phase transition is very sharp. In the case of quark matter a first-order transition is suggested at low temperatures ( MeV) which changes to a second-order one at higher temperatures. In contrast to the quark matter in the case of nucleon matter the thermodynamic variables show large discontinuities which is a clear indication for a first-order phase transition. In hot medium at intermediate temperature the nucleon soliton is more bound and less swelled than in the case of cold matter. At some critical temperature, which for nucleon matter coincides with the critical temperature for the phase transition, no more a localized solution is found. According to this model scenario one should expect a first-order phase transition from nucleon to quark matter. The results show that the hybrid model provides a consistent picture where the Fermi sea of quarks is reasonable only for hot matter with temperature and/or density around/above the critical values.
1 Introduction
At some finite density and/or temperature one generally expects a restoration of the chiral symmetry and a deconfinement, and hence a change of the structure of the hadrons immersed in a hot and dense medium. The bulk properties of hot nuclear and quark matter and especially the phase diagram of QCD is a topic of increasing interest since it is related to the evolution of the early universe after the Big Bang as well as to the processes in the interior of neutron stars. Rather encouragingly, direct experimental studies of such phenomena are now possible in the ultra-relativistic heavy-ion reactions accomplished at high-energy accelerators like RHIC, NA61 at the SPS as well as FAIR and NICA (see e.g. [1] and for a review [2]). It is illustrated on Figure 1 where the phase diagram of QCD is schematically shown.
Obviously a detailed understanding of this phase diagram is one of the contemporary challenges of both particle and nuclear physics. In particular, the order of the phase transition is still an open question. It seems to be generally accepted that below some critical point there is a first order-phase transition which above the critical point undergoes a crossover transition to a second-order one. Around the critical point one expect a mix of the two phases with different order parameters (see e.g. for an overview [3] and also [6]).
From the side of theory these expected low critical values mean that we have to deal with a non-perturbative phenomena. So far the analytical as well as the numerical (lattice QCD) methods (see review [2] and references therein) are still not developed enough to allow the solution of low-energy non-perturbative phenomena, especially if nucleons are involved. It is the reason for applying in the last two decades different effective models like e.g. NJL [7] and linear sigma model [5] to study these phenomena. These models are rather attractive since although lacking confinement they incorporate, similar to QCD, the chiral symmetry and also allows for its spontaneous breaking. It was the motivation to apply a chiral quark-meson model (see for review [6] and references therein) for these analyses. The model is able to reproduce quite reasonably the nucleon and delta properties as well as the corresponding formfactors in vacuum. Within the model the nucleon in vacuum is considered as a bound state of N, quarks coupled to the polarized Dirac sea. The mean-field approximation is used, which means that the meson quantum (loop) effects are not considered. Since the meson loop effects are dominant at low temperatures and vanishing density (pions are the lightest mode) the results are restricted mainly to the case of finite density and relatively large temperatures where the nucleon and quark degrees of freedom are most relevant. In order to consider the medium effects in the model a Fermi sea is added.
2 Quark-meson chiral model at finite temperature and density
The SU(2)-version of the NJL Lagrangian [7] contains chirally invariant local scalar and pseudoscalar four-quark interaction:
[TABLE]
where is the quark field, is the coupling constant, are the Pauli matrices in the isospin space and is the current quark mass taken equal for both up and down quarks. Applying the well-known bosonization procedure [8] the NJL model is expressed in terms of the auxiliary meson fields :
[TABLE]
Using the functional integral technique the quantized theory at finite temperature and density can be written in terms of the corresponding euclidean grand canonical partition function [9]:
[TABLE]
where is the volume of the system, is the inverse temperature and is the chemical potential. The integration over the quarks can be done exactly, whereas for the integration over the mesons we use a large saddle-point (mean-field) approximation. This means that the meson fields are treated classically - no meson loops are taken into account. Following [9] the integration over the imaginary time is replaced by a sum over fermionic Matsubara frequencies. Finally the effective action is expressed as
[TABLE]
The energies are eigenvalues of the one-particle hamiltonian .
[TABLE]
and are eigenfunctions. The saddle-point solution makes the effective action stationary
[TABLE]
The number of particles in the volume is kept fixed
[TABLE]
Here and are the “classical” values of the meson fields and is the chemical potential related to the number of particles .
The thermodynamic characteristics of a many-body system are specified by the thermodynamic potential
[TABLE]
It should be noticed that the saddle-point solution ( , ) minimizes not but the Helmholtz free energy
[TABLE]
with a constraint (7) and playing a role of a Lagrange multiplier.
In the mean-field approximation (leading order in the inverse meson propagator is given by the second variation of the effective action at the stationary point ( , ):
[TABLE]
The on-shell meson masses correspond to the poles of the meson propagator at . The physical quark meson coupling constants are given by the residue of the propagator at the pole
[TABLE]
Due to the local four-fermion interaction the lagrangian (1) is not renormalizable and a regularization procedure with an appropriate cut-off is needed to make the effective action finite:
[TABLE]
Actually only the part of the effective action coming from the Dirac sea (negative-energy part of the spectrum), is divergent and one needs to regularize it. Here the proper-time regularization scheme is used. The difference is finite (Fermi sea contribution) and does not need any regularization. Moreover, any regularization of the medium part would suppress the medium effects ([10]), since the positive part of the spectrum would be affected by the cutoff as well.
Two different scenario are considered: Fermi sea of quarks as well as Fermi sea of nucleons. In both cases the Dirac sea consists of quarks and it determines the vacuum sector. In both pictures the mesons appear as excitations but they are also directly coupled to the Fermi sea as well. In the case of Fermi sea of quarks at finite temperature/density single quarks are allowed to be excited leaving holes (antiquarks) in the Dirac sea. Apparently because of confinement this picture is applicable only at temperature/density above the critical ones - after the phase transition.
Since the Fermi sea contribution is finite it is straightforward to write the Fermi sea part in terms of nucleon.
[TABLE]
The energies are the solutions of the corresponding Dirac equation.
[TABLE]
The meson fields are coupled to the nucleons with a coupling constant which relates the nucleon mass to the non-zero expectation value of the scalar meson field (constituent quarks mass ) in vacuum.
[TABLE]
As in the quark spectrum, there is a gap of in the nucleon spectrum which separates the negative part of the spectrum from the positive one. The chemical potential is fixed by the baryon number (7).
3 Fixing the model parameters
In vacuum () the stationary conditions lead to a transitionally invariant solution and .
The parameters of the model, namely the current mass , the cutoff and the coupling constant In the vacuum sector, are fixed reproducing the physical pion mass MeV and the pion decay constant MeV. It leads to the well-known Goldberger-Treiman (GT) relation on the quark level
[TABLE]
and one also recovers the Gell-Mann-Oakes-Renner (GMOR) relation:
[TABLE]
The only free parameter remained is the constituent quark mass MeV is taken to reproduce properly the properties of a free nucleon (see [6] and references therein). The corresponding value of the proper-time cutoff used is MeV. The value is taken large enough not to influence to results.
4 Phase transition and meson properties at finite temperature and density
The results [11] for the constituent mass M as a function of temperature and density and are presented in Figure 2 and Figure 3. At both vanishing temperature and density it is fixed to MeV
At finite density the constituent mass is a non-monotonic function of temperature. It means that in hot matter the mass is less affected by the medium compared to the case of cold matter. At low values the curves are close to those of quark matter. However already at intermediate temperatures MeV they start to deviate significantly. At some critical values of temperature and/or density the constituent mass is reduced to the current mass which is an indication for a transition from Goldstone to the Wigner phase where the chiral symmetry is not spontaneously broken. Whereas the quark matter is quite soft against thermal fluctuations, the nucleon matter is much stiffer and the corresponding chiral phase transition is rather sharp. In order to determine the order of the phase transition we follow the Ehrenfest classification. Both the chiral condensate and the constituent quark mass show a discontinuity at the critical temperature of about 200 MeV. In fact, the system jumps between two minima. The latter suggests that in the case of nucleon matter we have to deal with a first order phase transition even in case of vanishing density and high temperature in contrast to the case of quark matter. The results of at and finite suggest a critical density ( with .) which in fact, as we will see later, is also the critical density for the delocalization of the soliton in cold nucleon matter.The same is valid for the temperature which makes this model picture consistent.
According to the present model picture, at some critical temperature and/or density one expects a chiral phase transition together with a delocalization transition from nucleon matter to quark matter. It can be clearly seen from the corresponding EOS (pressure versus baryon density) for different temperatures plotted in Figure 4. In this Figure we combine the results from nucleon matter (below the critical temperatures) in the hybrid model with those of quark matter after the transition. Since we do not include vector mesons in the hybrid model, we are not able to reproduce the nuclear matter saturation at zero temperature and finite density which in the Walecka approach [12] is due to the interplay between the -meson attraction and the -meson repulsion. All curves in Figure 4 show rapid change discontinuity at the delocalization transitions from nucleon to quark matter. Similar behavior can be seen in Figure 5 [11] for the energy density as a function of for different densities.
The phase diagrams for quark matter as well as for nucleon matter are shown in Figure 6 [11]. As can be seen in quark matter at temperatures MeV a second-order phase transition is predicted. At lower temperatures in quark matter a first-order transition is expected. For the nuclen matter the model predicts always a first-order transition. It should be noted also that below the critical temperature the critical density for the phase transition in nucleon matter is always smaller than one in quark matter. It means that at least according to the present model picture at temperatures below the critical two phase transitions are expected - a phase transition from nucleon to quark matter, where the chiral symmetry is still partially broken, followed by a phase transition from Goldstone to Wigner phase at a higher density in quark matter.
Concerning the order of the phase transition the model picture suggests at low temperatures ( MeV) a first-order transition from both the nucleon and from quark matter. At higher temperatures the picture is more complicated. For the phase transition from nucleon matter to quark one expects first-order whereas from quark matter the second-order phase transition is predicted. It means that at higher temperatures one could expect a mix of different phases with different order parameters.
Apart from the medium part written in terms of nucleons the inverse meson propagators in the nucleon medium have the same structure (29) as in the quark medium. The meson masses in the medium defined as the lowest zero solution of the inverse meson propagators at q = 0 are presented in Figure 4.
In chiral limit the pions remain Goldstone bosons and their mass remains almost unchanged also in the medium, whereas the sigma mass follows the constituent quark mass. Both GMOR and GT relations are also valid in the medium. Since in the Wigner phase the constituent quark mass vanishes up to the current mass , the mesons become degenerated in mass and appear as parity-doubled mesons.In the Wigner phase, after the deconfinement transition (suggested also by the lattice results) the meson states appear as resonances in the continuum and the poles of the meson propagators are complex. In the nucleon medium similar to the constituent mass M the mass of the sigma meson has a discontinuity at the critical temperature.
5 Nucleon as a non-topological soliton in the medium
The nucleon in a hot medium appears as B = 1 localized bound solution (soliton) of valence quarks interacting with Fermi and Dirac sea both getting polarized due to the interaction. The thermodynamic potential (effective action) includes an explicit valence quark contribution as well as Dirac and Fermi sea contributions:
[TABLE]
where is the proper-time regularization function of the form (12). The meson fields are assumed to be in a hedgehog form
[TABLE]
restricted on the chiral circle.
[TABLE]
We use a numerical self-consistent iterative procedure (see [6] and references therein) solving in an iterative way the Dirac equation together with the equations of motion of the meson fields with the constraint = const. The latter is a condition fixing the chemical potential . In the case of fixed and the proper way to describe the equilibrium state of a thermodynamic system is to use the free energy (9). Hence, the energy of the B = 1 soliton (effective soliton mass) is given by the change of the free energy when the valence quarks are added to the medium. Subtracting the free energy of the unperturbed Fermi and Dirac sea (translationally invariant medium solution), the soliton energy is given by the sum of the energy of the valence quarks and the contributions due to the polarization of both continua:
[TABLE]
where is the baryon density. Accordingly the soliton baryon density distribution is split in valence, sea and medium parts.
[TABLE]
We do not find a localized solution (soliton) at temperatures larger than a critical value MeV. It means that at sufficiently large values the temperature effects simply disorder the system and destroy the soliton. Within the present model picture this may be interpreted as an indication for a delocalization of the nucleon in a hot medium. However, one should keep in mind that the model lacks confinement.
All detailed results can be found in [11]. In Figure 8 the temperature dependence of the calculated soliton energy are presented for the different quark matter density. In case of finite and zero the soliton energy shows almost no reduction for temperatures not close to the critical one, whereas in the opposite case of finite and zero the soliton energy is linearly decreasing with the density. In a hot medium (both and finite) the temperature clearly suppresses the finite density effects and stabilizes the soliton. Further, at densities larger than two times the soliton exists only at intermediate temperatures 100 MeV MeV. It means that according to the present model calculations, the soliton is more stable in hot matter than in cold matter. In fact, this can be easily understood. Because of a gap in the quark spectrum, , the Dirac sea is much less affected by the temperature than the Fermi sea (positive-energy part of the spectrum). At temperatures close to the critical one, the thermal fluctuations become comparable with the chiral order parameter (chiral condensate ). They completely disorder the system and destroy the soliton.
In Figure 9 the temperature dependence of the calculated soliton energy for the nuclear matter density is compared with those of quark matter. The two curves have similar trends at intermediate values of and start to deviate at close to the critical temperature . The soliton in nucleon matter is more bound and less affected by the temperature.
In order to illustrate the change of the soliton structure in a hot medium we also plot the soliton square radius as a function of density for different temperatures in Figure 10. All curves show a clear trend to grow rapidly at temperatures close to the critical values which is an indication for a delocalization of the soliton. At both finite density and temperature, however, the radius is smaller than in the case of cold matter which is a sign for stabilization of the soliton in a hot medium compared to the case of a cold medium.
The calculated m.s. soliton radius for both quark and nucleon matter is presented in Figure 11. Close to the critical temperature it starts to grow but it is much less pronounced than in the case of quark matter.
It should be noted also that the coupling constant defined as
[TABLE]
stays almost constant (see [11]) which means that the relation (23) is a good approximation also in a hot medium.
6 Conclusion
The bulk thermodynamic characteristics, meson properties and the properties of the nucleon as a soliton of valence quarks with Dirac and Fermi seas are studied in a hot medium in the framework of an effective quark-meson theory. The Fermi sea of nucleons as well as of quarks is considered. At some critical values of temperature and/or density chiral phase transitions from the Goldstone to the Wigner phase are suggested. The quark matter is quite soft against thermal fluctuations whereas the nucleon matter is much stiffer and the corresponding phase transition is rather sharp. In the quark matter at low temperatures (below MeV) a first-order phase transition is expected. At higher temperatures it changes to a second-order phase transition. In contrast to the quark matter in the case of nucleon matter the model suggests a first-order phase transition even in case of vanishing density and high temperature. According to this model picture at higher temperatures a mix of different phases with different order parameters is expected. The bulk thermodynamic characteristics in nucleon and quark matter are different and at least according to the present results one has to consider also the baryon degrees of freedom in order to get a complete picture.
The nucleon as a soliton in medium is getting swelled and its mass is reduced as well. In the baryon medium the soliton is less affected by the medium. At finite density the temperature stabilizes the soliton. At some critical values of temperature and/or density the nucleon as a soliton disappears. In the model this delocalization means a transition from nucleon to quark matter. The critical values for the delocalization of the soliton are same as for the phase transition from nucleons to quarks in the nucleon matter which makes this model picture rather consistent.
Acknowledgements
The author is greatful for the financial support of the Bulgarian Science Fund under Contract No. DFNI–T02/19.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] The Frontiers of Nuclear Science, A Long Range Plan (2008)
- 2[2] C. Ratti, Rept.Prog.Phys. 81 (2018) 84301
- 3[3] M. Stephanov, XXIV International Symposium on Lattice Field Theory, July 23-28, Tucson Arizona, USA, Po S LAT 2006 (2006) 024
- 4[4] A. Ohnishi, J.Phys.Conf.Ser. 503 (2014) 12009
- 5[5] M. Gell-Mann and M. Levy, Nuovo Cim. 16 (1960) 705
- 6[6] Chr. V. Christov, A. Blotz, H.-C. Kim, PV. Pobylitsa, T. Watabe, Th. Meissner, E. Ruiz Arriola and K. Goeke, Prog. Part. Nucl. Phys. 37 (1996) 1
- 7[7] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 354
- 8[8] T. Eguchi and H. Sugawara, Phys. Rev. D 10 (1974) 4257; T. Eguchi, Phys. Rev. D 14 (1976) 2755
