Next Article in Journal
Editorial for the Special Issue “The Modern Physics of Compact Stars and Relativistic Gravity 2017”
Next Article in Special Issue
Strange Stars in the Vector Interaction Enhanced Bag Model
Previous Article in Journal
Nonlocal Scalar Quantum Field Theory—Functional Integration, Basis Functions Representation and Strong Coupling Expansion
Previous Article in Special Issue
Simulating Binary Neutron Stars with Hybrid Equation of States: Gravitational Waves, Electromagnetic Signatures and Challenges for Numerical Relativity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Meson Condensation

by
Massimo Mannarelli
INFN-Laboratori Nazionali del Gran Sasso, Via G. Acitelli, 22, I-67100 Assergi (AQ), Italy
Particles 2019, 2(3), 411-443; https://doi.org/10.3390/particles2030025
Submission received: 5 August 2019 / Revised: 4 September 2019 / Accepted: 9 September 2019 / Published: 13 September 2019
(This article belongs to the Special Issue Dense QCD and neutron stars)

Abstract

:
We give a pedagogical review of the properties of the various meson condensation phases triggered by a large isospin or strangeness imbalance. We argue that these phases are extremely interesting and powerful playground for exploring the properties of hadronic matter. The reason is that they are realized in a regime in which various theoretical methods overlap with increasingly precise numerical lattice QCD simulations, providing insight on the properties of color confinement and of chiral symmetry breaking.

1. Introduction

The great success of the Standard Model of particle physics relies on the possibility of making accurate and testable predictions that are in agreement with increasingly precise experimental data. Despite this success, many aspects of the Standard Model are still not completely clear. Among these, there are the mechanisms of color confinement and of chiral symmetry breaking ( χ SB) of the strong interaction. The typical energy scales of confinement and χ SB pertain to the nonperturbative region of quantum chromodynamics (QCD), which makes their study extremely challenging.
For some theoreticians, the confinement and the χ SB mechanisms are uninteresting because they are details of a robust theoretical construction, thus sooner or later they will be fully understood; for others (including the author), unraveling the origin of these mechanisms is of the utmost importance for a comprehensive understanding of QCD; for all, it is still unclear which is the path that can bring us to a full understanding of these mechanisms. My view is that any path, as far it is physically sound, should be explored and tested. This brief review is about one those paths, exploring the behavior of matter when there is an asymmetry in the number of particles with different isospin and/or strangeness. This seems a promising direction because a number of theoretical methods can be used for studying these phases. Comparing the results obtained by different methods helps us to check their consistency and the degree of the reached accuracy.
To clarify the setting we report in Figure 1 a sketch of the so-called QCD phase diagram: a grand canonical description of the phases of hadronic matter as a function of the hadronic temperature, T, of the isospin chemical potential, μ I , and of the baryonic chemical potentials, μ B . The total baryonic density is determined by μ B , while μ I describes the isospin asymmetry, say due to a different number of up and down quarks. If it were possible we would have added a further axis, μ S , indicating the strangeness content. The blue region corresponds to a gas of confined hadrons with a chiral broken symmetry. At large energy scales quarks and gluons should be liberated [1] realizing different phases. We have indicated three of them. The quark-gluon plasma (QGP), realized at large temperature, is asymptotically a gas of quarks and gluons that becomes strongly interacting for the temperature reachable in heavy-ion collisions, see for instance [2,3,4]. At large μ B we expect that deconfined quarks fill their Fermi spheres and that the color interaction drives the formation of Cooper pairs in a BCS-like color superconducting phase (CSC), see [5,6,7] for reviews. At large and positive μ I we expect to populate the u and d ¯ states with the color interaction inducing the formation of π + states that will eventually condense. A negative isospin chemical potential does instead favor the formation of π states. Whether these phases persist to the hadron gas surface of Figure 1 depends on the non-perturbative properties of QCD. This diagram is somehow our starting point to set the stage, we will review it in Section 5 feeding in the up to date results. For the time being we note that there are two uncontroversial results: The critical temperature T c 160 MeV has been experimentally investigated at RHIC and at LHC and precisely determined by LQCD simulations [8,9] to correspond to a analytic crossover. The transition to the pion condensed phase at T = 0 is a second order phase transition happening exactly at μ I = m π [10,11].

Outside the Beta-Equilibrated Sheet

The three chemical potentials, μ B , μ I and μ S are not three independent quantities because the weak interactions regulate the isospin and strangeness content of matter at a given baryonic density. Moreover, matter must be electrically neutral and the strong interactions can modify the dispersion laws of quasiparticles. Therefore, we should indicate in Figure 1 a beta-equilibrated sheet μ I = f ( μ B , T ) , where f is some equation of state giving the isospin asymmetry at any temperature and baryonic chemical potential. The beta-equilibrated sheet corresponds to the configuration realized in long lived systems, as in compact stars. It is however instructive to consider configurations outside this surface, for three main reasons. The first is that we do not actually know f, except in restricted energy regions: at small energy scales by nuclear experiments, and at asymptotically large energy scales by perturbative QCD; the second is that in the theoretical investigation we can turn off the weak interactions, thus we can compare the outcomes of different theoretical methods outside the beta-equilibrated sheet to test their robustness and consistency. Finally, working on the beta-equilibrated sheet is interesting for studying the properties of dense nuclear matter in compact stars, but makes the problem so complicated that it is presently hard to make predictions.
In the first papers discussing the pion condensation [12,13], Migdal tried to work on the beta-equilibrated sheet considering how nuclear matter mechanisms could make the π 0 stable states. Then, different authors [14,15,16] considered the possible mechanisms for in-medium stabilization of charged pions by a softening of the pion spectrum by the p-wave pion-nucleus interaction. The in-medium pion dispersion law was assumed to be
ω 2 ( k ) = m π 2 + 0.7 n m π ω + k 2 ( 1 6 n ) ,
where n is the baryonic density in units of fm 3 . The exciting result is that a gapless mode appears for n > 1 / 6 fm 3 , thus very close to the nuclear saturation density n 0 0.16 fm 3 , at a momentum k = m π / 6 n 1 . Therefore, a transition to a superfluid phase was expected just above the nuclear saturation density. However, it was pointed out by Migdal [17] that the π π repulsive interaction may qualitatively change this result and he argued in favor of π 0 condensation as well as stable π + π molecules. Still today we have not solved this problem, but most of the theoretical works are now directed to understanding what happens outside the beta-equilibrated sheet. This new approach has allowed Son and Stephanov [10] to qualitatively and quantitatively assess the main properties of the pion condensed phase by the use of a simple approach based on the modelization of QCD by chiral perturbation theory ( χ PT) [10,11].
Actually, there is a number of theoretical approaches that can be used. In principle, any information on the phase diagram of Figure 1 could be obtained introducing in the QCD action a chemical potential for the charge of interest. The QCD Lagrangian turns to be
L QCD = ψ ¯ γ μ i D μ M ψ 1 4 F μ ν a F a , μ ν ,
where F μ ν a with a = 1 , , 8 are the color gauge field strengths, ψ T = ( u , d , s ) is the spinor describing up, down and strange quarks (with suppressed color and spinorial indices), the bare quark masses are collected in the mass matrix
M = diag ( m , m , m s ) ,
where we have assumed degenerate light quark masses, and the covariant derivative
D μ = μ + i g A μ i 2 v μ ,
includes both the minimal interaction with the gauge fields, A μ , and with the external source
v μ = 2 μ δ μ 0 ,
with the chemical potentials collected in the matrix
μ = diag μ u , μ d , μ s = μ B μ S 3 I + μ I T 3 + 2 μ S 3 T 8 ,
where T 3 and T 8 are the two diagonal S U ( 3 ) generators, I is the 3 × 3 identity matrix and we have parameterized the quark chemical potentials as μ u , d = μ B / 3 ± μ I / 2 and μ s = μ B / 3 μ S . To determine the phase diagram in Figure 1 one should obtain the behavior of the chiral, pion and diquark condensates, respectively given by
σ ψ ¯ ψ ,
π a ψ ¯ σ a γ 5 ψ ,
Δ ψ C γ 5 ψ ,
where σ a with a = 1 , 2 , 3 are the Pauli matrices, as a function of T, μ B and μ I .
Given the nonperturbative character of QCD at the energy scales of the various phase transitions in Figure 1, the Lagrangian in Equation (2) is of little direct use. Various theoretical approaches have been developed, including linear sigma models and chiral perturbation theory ( χ PT) [10,11,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33], the Nambu-Jona Lasinio (NJL) models [34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56], the quark-meson models [57,58,59,60,61,62], the random matrix model [63,64], the AdS/QCD model [65] and perturbative QCD (pQCD) (with diagrams resummation) [66,67]. A guiding role in this forest of theoretical approaches is played by the lattice QCD (LQCD) simulations [68,69,70,71,72,73,74,75,76,77,78,79,80,81,82], which provide a powerful tool for a numerical check and for exploring non-perturbative QCD. As we shall discuss in some detail in Section 4.3, the grand-canonical LQCD simulations at finite baryonic density and/or strangeness density are hampered by the so-called sign problem, but are feasible at μ B = μ S = 0 and μ I 0 [68]; moreover it is possible to simulate an ensemble of kaons by the canonical LQCD approach [74,75], corresponding to a system at nonvanishing strangeness density, as discussed in more detail in Section 4.3. This places the study of the meson condensation on a firmer ground with respect to the study of the phases at large baryonic density: the various theoretical models (we will mostly focus on χ PT and the NJL model) give a qualitative and semiquantitative description of the meson condensed phase; the LQCD simulations provide numerical evidence for the proposed phase transitions and for the meson properties. The theoretical understanding at large μ B cannot count on experimental data nor on numerical simulations. The phases realized at high baryonic densities could be relevant for dense stellar objects, but it is hard to obtain constraints on the microscopic properties of matter from the macroscopic properties of compact stars [83,84].
This review is organized as follows. In Section 2 we investigate the stability of pions in nuclear matter by a simplified non-interacting model description of hadronic matter. This model, mainly used in the 1960s and in the 1970s, serves as a guide for understanding by simple qualitative reasoning how meson condensation can occur and why it is unclear whether it be realized in compact stars or any other stellar object. In Section 3 we use an argument based on group theory to derive the phase diagram of the meson condensed phases. This is a useful result because any other theoretical modeling is expected to reproduce this phase diagram. In Section 4 we report and compare the results obtained by three different approaches: χ PT, NJL and LQCD, showing that the obtained results are in qualitative and quantitative agreement for T = 0 and μ I 2 m π . At nonvanishing temperature the three approaches give similar qualitative results, but more work is needed to reconcile the χ PT and NJL methods with the precise numerical results of the LQCD simulations. We conclude, with a new discussion of the QCD phase diagram, in Section 5.

2. The Early Works and Models

To understand why the pion condensation in dense hadronic matter is controversial [85] we scrutinize the condition for the meson stabilization against weak decays by a simple non-interacting gas approximation (NGA). This is a mean field model based on the assumption that neutrons, protons and electrons behave as independent Fermi gases of quasiparticles. As we shall see below, the result of these kind of models is that fermions are favored with respect to mesons [86,87,88,89,90]. The pion states can only be populated at non extreme densities if the strong interaction modifies the pion spectrum.

2.1. The Equilibrium Configuration

The condensation mechanism of any type of particle relies on three basic requirements:
  • The particles must be bosons, as 4 He atoms, or boson-like, as Cooper pairs in the BCS theory
  • The system has to be sufficiently cold: the particle condensation can be disrupted by the thermal disorder
  • The particles must be stable.
Since mesons are bosons, they satisfy the first requirement. We can also imagine that they can be produced in a relatively cold environment, as in the core of neutron stars [83]. The third point is typically neglected in ultracold atom physics, because experiments are done with stable atoms [91]. In contrast, all mesons in vacuum are unstable. The point is whether medium effects can stabilize mesons or not.
The first papers discussing stable pions in dense nuclear matter appeared in the 1960s [86,87,88,89] considering a medium of catalyzed nuclear matter: neutral matter consisting of nucleons and leptons in electroweak equilibrium. In the Fermi gas approximation, the equilibrium distribution is determined by the Urca process
n p + e + ν ¯ e p + e n + ν e ,
that with the assumption of neutrino transparency implies that
μ n = μ p + μ e ,
where μ i , with i = n , p , e , are the appropriate chemical potentials, see [92] for the temperature corrections to the Fermi gas approximation. The electrical neutrality requires that
n p = n e ,
where n p and n e are the number densities of protons and electrons, respectively.
In the NGA, the number density of each component is
n i = 2 0 p f , i d 3 p ( 2 π ) 3 = p f , i 3 3 π 2 ,
with p f , i the Fermi momentum. The neutrality condition Equation (12) then implies that the Fermi momenta of electrons and protons are the same, that is
p f , p = p f , e μ e ,
where in the last expression we have neglected the electron mass. For free Fermi gases, without in-medium effects, we obtain using Equation (11) that
p f , n 2 μ e ( μ e + μ e 2 + m n 2 ) ,
where m n is the nucleon mass (we have neglected the small mass difference between protons and neutrons). Upon substituting these expressions in Equation (13) we can express any number density as a function of the electron chemical potential; in particular, the total number density as an increasing function of the electron chemical potential. We can now include different hadronic states. We determine and compare, for reasons that will soon become clear, the threshold electron chemical potential for the appearance of π and Σ .
In vacuum π decays to leptons, mostly in muon and in its antineutrino. Since muons decay in electrons and neutrinos, we simplify the discussion considering the process (solid red line in Figure 2)
π e + ν ¯ e ,
where we assume that neutrinos are not trapped. This process can be Pauli blocked as depicted in the right side of Figure 2, meaning that it is in equilibrium with the electron decay process (dashed blue line in Figure 2)
e π + ν e ,
in the configuration with
μ e π = m π 135 MeV .
We now proceed with a similar analysis for the stability of the Σ baryonic resonance. In the quark model it is a ( d d s ) state with mass, m Σ 1.2 GeV, weakly decaying by
Σ n π n + e + ν ¯ e ,
in about 10 10 s. The Σ states can be populated by the electron capture process
e + n Σ + ν e ,
with eventually an additional spectator neutron to ensure the energy-momentum conservation. The equilibrium is reached at
μ e Σ = 2 m Σ m Σ 2 + 3 m p 2 3 126 MeV ,
and since μ e π > μ e Σ the Σ states appear at a lower density than the π states. This result is somehow surprising: at high density the heavier Σ states are favored over the lighter mesonic π states. This happens because the production channels of these two hadrons are different and because the Σ is not much heavier than nucleons. The density for the appearance of the Σ states in the NGA is pretty large, about 4 n 0 . Once the Σ states are populated they must be included in the neutrality condition in Equation (12), forbidding the appearance of pions up to n 300 n 0 , see for instance [90] and references therein.

2.2. Including in-Medium Effects

Although in the free gas approximation the π is not energetically favored with respect to Σ , the difference between the two critical chemical potentials is small, μ e π μ e Σ 10 MeV m π . Thus, a slight change of their effective masses may invert this result: charged pion states could be populated first.
To gain insight, we assume that the the strong interactions cause a constant Fermi energy shift [90], independent of the particle momentum. In the NGA framework, for any baryon i, we define the effective chemical potential
μ eff , i = μ i B ¯ i ,
where B ¯ i is the constant chemical energy shift due to the strong interactions, while for leptons B ¯ lepton = 0 . The chemical equilibrium, Equation (11), now reads
μ eff , n = μ eff , p + μ e ,
where we have assumed as before that neutrinos escape: the medium does not trap neutrinos. Turning to the problem of populating π and Σ states, the charged pion states are now populated (do not decay in leptons) at
μ e π = m π B ¯ π ,
thus a positive B π favors the appearance of π . The Σ appears at the critical electron chemical potential in Equation (21) with
m Σ m Σ B ¯ Σ + B ¯ p ,
therefore assuming small binding energies we find that pions appear first ( μ e π μ e Σ ) for
B ¯ π 0.47 ( B ¯ Σ B ¯ p ) 10 MeV ,
which is a condition depending on three in-medium quantities.
The point of this simple model is that the appearance of pions may or may not happen depending on quite small parameters that are not under quantitive control. In particular, it is not clear in which direction the medium effects go [90]. Any result obtained in this way can hardly stand firm against scrutiny and it is doomed to be troublesome, because based on a too naive modeling and/or on extrapolating nuclear matter properties at least to 2 3 n 0 [19].

3. Group Theory Analysis

In recent years a simpler approach to the mechanism of meson condensation has been developed, disentangling it from weak equilibrium and strong interaction effects. We illustrate it for pions. Pions are an isospin triplet: they are the three eigenstates of an I = 1 multiplet with different I 3 projection. In vacuum the masses of the charged pions are degenerate, however this degeneracy is removed by the isospin chemical potential as in the Stark-Lo Surdo effect. In particular, we expect that the energy levels split as follows
E π 0 = m π 2 + p 2 ,
E π = + μ I + m π 2 + p 2 ,
E π + = μ I + m π 2 + p 2 ,
which are valid for | μ I | < m π , because for | μ I | > m π one of the two charged mesons becomes massless. This is an indication of a spontaneous symmetry breaking (SSB), with the resulting massless mode to be identified with the Nambu-Goldstone boson (NGB) associated to a broken generator. The best tool for exploring the (global) symmetry breaking is group theory. The good thing about this method is that it gives robust results. Its main limitation is that it does not provide a microscopic mechanism for the occurrence of the symmetry breaking.

3.1. Global Symmetries of QCD

In Figure 3 we report the global symmetry breaking pattern which is relevant for meson condensation. We start assuming three flavor massless quarks with μ I = μ S = 0 described by an Hamiltonian, H , with global symmetries
G = S U ( 3 ) L × S U ( 3 ) R U ( 1 ) Q × U ( 1 ) B ,
where we have not included the U ( 1 ) A anomalous group and we have specified that the electromagnetic gauge group, U ( 1 ) Q , is a subgroup of the considered symmetries. The baryon symmetry group, U ( 1 ) B , can be broken by the formation of quark Cooper pairs in the CSC phase, see Figure 1. In the following we will not consider this possibility, focusing on the chiral symmetry, S U ( 3 ) L × S U ( 3 ) R , corresponding to the invariance of massless QCD with respect to left- and right-handed rotations, indicated respectively with U L and U R in Figure 3, of the quark fields.
In vacuum the chiral condensate locks the chiral rotation to the vectorial flavor group by the SSB
G S U ( 3 ) V U ( 1 ) Q × U ( 1 ) B ,
corresponding to the simultaneous rotations of left- and right-handed quark fields. The ground state is still invariant under rotation of quark flavors but the left and right handed rotations are now locked: the quark transformation leaving the vacuum invariant is U L = U R . According to the Goldstone’s theorem this symmetry breaking pattern results in 8 NGBs associated to the broken generators, see the left side of Figure 4. Actually, the bare quark masses explicitly break the chiral symmetry, meaning that these modes are massive pseudo NGBs, identified with the pseudoscalar meson octet. Assuming that the light quark masses are degenerate the isospin symmetry is preserved, meaning that the resulting symmetry is now (second row of Figure 3)
S U ( 2 ) I × U ( 1 ) Y U ( 1 ) Q × U ( 1 ) B ,
where S U ( 2 ) I is the isospin symmetry group and U ( 1 ) Y is the symmetry group generated by T 8 . Thanks to the isospin invariance the pions are degenerate, with mass m π 135 MeV. The kaons are grouped in two isodoublets with degenerate masses m K 500 MeV. The kaon masses differ from the pion masses because kaons involve a strange quark. These meson masses are reported on the central part of Figure 4. The last member of the octet, which is not reported in Figure 4, is the η field, with mass
m η 4 m K 2 m π 2 3 ,
by the Gell Mann-Okubo relation.
The chemical potentials induce two different symmetry breakings, one is explicit, at the Lagrangian level, while the second one is a SSB, causing the meson condensation. At the Lagrangian level, the chemical potentials explicitly break the Lorentz boost invariance, due to the presence of a privileged reference frame, the one in which particles are at rest with the medium. Space rotations and translations are instead unaffected: the considered medium is isotropic and homogeneous. The chemical potentials also explicitly break the charge symmetry, meaning that charged conjugated states will have different masses. This effect can also be seen in a more detailed way scrutinizing how the chemical potentials in Equation (6) break the flavor symmetry. Since
[ H , T 3 ] = 0 and [ H , T 8 ] = 0 ,
while the commutator of the Hamiltonian with any other S U ( 3 ) generator is nonzero, the explicit symmetry breaking is (third row in Figure 3)
S U ( 2 ) I × U ( 1 ) Y U ( 1 ) Q × U ( 1 ) B N = U ( 1 ) I × U ( 1 ) Y U ( 1 ) Q × U ( 1 ) B
where U ( 1 ) I and U ( 1 ) Y are the symmetries generated by T 3 and T 8 , respectively. The remaining symmetry group can be viewed as generated by the independent phase rotations of the three flavor fields. This residual symmetry is extremely important, because it proves possible a SSB and a phase transition between the normal phase, characterized by the symmetry N in Equation (35), to a superfluid phase, with a reduced global symmetry. The meson condensation, last part of the diagram in Figure 3, is indeed related to the SSB that locks the phases of quarks with different flavors. It is induced by a large μ I or μ S .

3.2. Phases of Condensed Mesons

We can figure out the conditions for the final SSB of Figure 3 by inspecting the meson mass spectrum. The isospin and strangeness chemical potentials remove the degeneracy within isomultiplets and between states with different hypercharge. Since
[ H , T 3 ] = 0 and [ H , T 2 ] = 0 ,
in the normal phase we can still group eigenstates as isospin multiplets. By the lemmas of Schur, these states will not be degenerate in mass: there is a mass splitting within members of isomultiplets, proportional to μ I , and there is a mass splitting between states with different hypercharge, proportional to μ S . The only unaffected states are the π 0 and η states, because they have both vanishing isospin and hypercharge. Since the mass splitting is given by the corresponding charges we have that
m π 0 = m π ,
m π ± = m π μ I ,
m η = 4 m K 2 m π 2 3 ,
m K ± = m K 1 2 μ I μ S ,
m K 0 / K ¯ 0 = m K ± 1 2 μ I μ S ,
where m π and m K indicate the meson masses for μ I = μ S = 0 .
To elucidate the mass splittings we show on the right of Figure 4 the results obtained for μ I = 80 MeV and μ S = 90 MeV. As discussed above, the first mass splitting is induced by the nonvanishing quark masses. The second splitting is due to the considered values of the chemical potentials, which completely remove the level degeneracies. The shown hierarchy is only one of the possible ones: different values of μ I and/or μ S may imply different level splittings, possibly with kaons lighter than pions.
The pion masses are independent of μ S and we see from Equation (38) that for
μ I = | m π | ,
one of the charged mesons becomes massless. Analogously, for | μ S | = m K | μ I | / 2 one of the kaons becomes massless. The normal phase will persists until no mode is massless, that is for
| μ I | < m π and | μ S | < m k | μ I | 2 ,
which corresponds to the irregular hexagon (solid red line) in Figure 5.
This phase diagram was first derived in [11] in the χ PT framework. At the boundary of the hexagon we expect a massless mode, signaling a SSB. This massless mode eventually condenses if the temperature is below the relevant critical temperature. When meson condensation occurs, the vacuum will have a nonvanishing charge associated to one of the non-diagonal generators of S U ( 3 ) characterizing the condensed meson.
For the description of the broken phases, we first note that the μ I and μ S Lagrangian terms are proportional to T 3 and T 8 , respectively, see Equation (6). Since these two generators form the Cartan subalgebra of S U ( 3 ) , it turns out useful to consider the three SU ( 2 ) Lie subalgebras of SU ( 3 ) generated by
T ± , T 3 , V ± , V 3 , U ± , U 3 ,
where the step operators are respectively
T ± = T 1 ± i T 2 , V ± = T 4 ± i T 5 , U ± = T 6 ± i T 7 ,
and
Y = 2 3 T 8 , K = T 3 1 3 T 8 , Q = T 3 + 1 3 T 8 ,
are the corresponding weights.
Let us now consider the case in which one of the charged pions, say the π + , becomes massless and condenses. The spontaneously broken generator is T + , which locks the phases of up and down quarks. Since [ T + , Y ] = 0 , it follows the SSB
U ( 1 ) I × U ( 1 ) Y U ( 1 ) Q × U ( 1 ) B U ( 1 ) Y × U ( 1 ) B U ( 1 ) Q ,
meaning that the massless mode is the NGB associated to U ( 1 ) I breaking. We remark that this superfluid mode actually signals the transition to a superconducting phase, because the U ( 1 ) Q symmetry is broken. Since [ T + , Q ] 0 the charged pion states can mix and indeed the massless mode coincides with the π + only at the phase transition point. In general it will be a superposition of the two charged pion states. Analogous results hold if the π becomes massless and condenses, with T the broken generator.
The second possibility is that the K + becomes massless and condenses; the spontaneously broken generator is V + , which locks the phases of up and strange quarks. In this case [ V + , K ] = 0 , with the resulting SSB
U ( 1 ) I × U ( 1 ) Y U ( 1 ) Q × U ( 1 ) B U ( 1 ) K × U ( 1 ) B U ( 1 ) Q ,
meaning that there is a residual global symmetry associated to the weight K. In this case the isospin and hypercharge phases are locked. The system is a kaon superconductor, because the condensed meson is charged. Analogous results hold if the K becomes massless and condenses.
The third case is the condensation of one of the neutral kaons. If the K 0 becomes massless and condenses, then the spontaneously broken generator is U + , with down-strange quark phase locking. Since [ U + , Q ] = 0 , the resulting SSB is
U ( 1 ) I × U ( 1 ) Y U ( 1 ) Q × U ( 1 ) B U ( 1 ) Q × U ( 1 ) B ,
meaning that the electromagnetic gauge field is unbroken. Therefore, the system is a kaon superfluid where the superfluid mode is given by the mixing of the two neutral kaons.
Let us analyze the order of the transition lines in Figure 5. The χ SB and the meson condensation mechanisms are independent: the χ SB is related to the locking of the phases of left- and right-handed quarks, while the meson condensation locks the phases of quarks with different flavors. This means that the chiral condensate and the meson condensate can coexist, therefore the irregular hexagon should be a second order phase transition line. On the other hand, it is not possible to have the simultaneous condensation of say π + and K + because both condensates involve an up quark. In group theory this is related to the fact that the two generators T + and V + do not commute. We expect that the presence of a meson condensate excludes the other, implying the first order phase transition shown in Equation (5) as dashed lines. The simultaneous condensation can only happen at the vertices of the hexagon or close to the first order phase transition lines if inhomogeneous phases are realized. This has been preliminarily explored by canonical LQCD simulation in [75] and by χ PT in [30].
To summarize, any meson condensate tilts the vacuum in a certain direction having a residual U ( 1 ) B × U ( 1 ) symmetry which is generated by the baryonic charge and by one of the weight operators in Equation (46). The low-energy spectrum consists of one NGB. This NGB has the same quantum numbers of one of the standard mesons at the phase transition point, but then in the superfluid phase it mixes with its charged conjugated state as shown in Figure 6, see also [26].

4. Modern Approaches

We now present two different approaches to the meson condensation, one is based on the effective field theory χ PT description of mesons, the second on the NJL modelization of the strong interaction by a contact term. Then we compare the results of these methods with those of the pertinent LQCD numerical simulations. We restrict to consider an homogeneous and static medium; for inhomogeneous phases see [7,61,93,94,95].

4.1. Chiral Perturbation Theory

The χ PT Lagrangian is an extremely powerful tool for systematically describe the strong interactions between hadrons [96,97,98,99,100,101,102,103,104]. Remarkably, χ PT can also be used to study a variety of gauge theories with isospin asymmetry, including 2 color QCD with different flavors [105,106,107,108,109,110,111]. Quite generally, any χ PT realization is based on two key ingredients: the global symmetries of the studied theory and an appropriate low momentum expansion.
The relevant global symmetry of QCD for constructing the χ PT Lagrangian is the chiral symmetry
S U ( N f ) L × S U ( N f ) R ,
with N f the number of relevant quark flavors. The meson fields are collected in the unimodular Σ field, transforming as
Σ L Σ R ,
where L S U ( N f ) L and R S U ( N f ) R . Based on the chiral symmetry and this transformation property one builds the most general Lagrangian at the given order in the momentum expansion, assuming that the meson momenta satisfy
p Λ χ ,
where Λ χ 1 GeV is the χ PT breaking scale. At each order in the momentum expansion the chiral symmetry fixes the form of the various Lagrangian pieces but the pre-factors, the so-called low energy constants (LECs), must be determined by different means. The leading O ( p 2 ) χ PT Lagrangian [11,97,103] describing the in-medium pseudoscalar mesons is given by
L = f π 2 4 Tr ( D ν Σ D ν Σ ) + f π 2 4 Tr ( X Σ + Σ X ) ,
where the trace is in flavor space and the auxiliary field X transforms as Σ . The locking of the chiral rotations to the vector S U ( N f ) V group is induced by the vev of X, see for example the discussion in [98,103], which is usually written as X = 2 B M with M the quark mass matrix defined in Equation (3). The two LECS, B and f π , can be fixed by the vacuum properties [97,99,100,103,104]; for instance B from the mass relations m π 2 = 2 B m and m K 2 = B ( m + m s ) , while f π 93 MeV from the weak pion decays.
The adjoint covariant derivative in Equation (53) takes into account the minimal coupling of the meson fields with gauge fields, external currents and the effect of different chemical potentials [97,99,101]. Pretty much as the covariant derivative in the quark sector, see Equation (4), we define it as
D ν Σ = ν Σ + i 2 [ v ν , Σ ] ,
where the external current, v ν , is given in Equation (5); here we can take μ = diag ( μ I / 2 , μ I / 2 , μ S ) because mesons have no baryonic charge. This leading order (LO) χ PT Lagrangian is sufficient to accurately describe the phase structure of QCD at μ I 2 m π [10,11,29], including finite temperature effects [22,23,24] and pion in-medium stability [26,28].

4.1.1. Ground State

The ground state can be variationally determined by replacing Σ Σ ¯ in Equation (53), where Σ ¯ is the time independent and homogeneous vacuum expectation value of the meson fields. The resulting static Lagrangian
L static = f π 2 4 Tr μ , Σ ¯ μ , Σ ¯ + f π 2 B 2 Tr M ( Σ ¯ + Σ ¯ ) ,
must be maximized to determine the vacuum configuration. Specifically, one can parameterize the most general S U ( 3 ) vev as
Σ ¯ = e i α n · λ ,
where λ i , for i = 1 , , 8 are the Gell-Mann matrices, and the tilting angle α and the eight-dimensional unit vector n are the variational parameters that should be obtained maximizing Equation (55). The normal phase is easily described by α = 0 , thus Σ ¯ = Σ ¯ N = I , but in general one should maximize Equation (55) with respect to eight independent parameters, which is a rather formidable task. As far as I know, nobody has ever considered this procedure, and for a good reason. From the insight gained in Section 3 by the group theory analysis, we expect that in each different meson condensed phase the vacuum is rotated in a specific way. In the pion condensed phase the only nonvanishing components should be ( n 1 , n 2 ) , in the charged kaon condensed phase ( n 4 , n 5 ) , while in the neutral kaon condensed phase ( n 6 , n 7 ) . One can further simplify the ansatz observing that the potential must have a flat direction, the one spanned by the NGB.
To elucidate these aspects let us first focus on the N f = 2 case, with the most general ansatz
Σ ¯ = e i α n · σ = cos α + i n · σ sin α ,
where n = ( n 1 , n 2 , n 3 ) . The static Lagrangian now reads
L static = f π 2 m π 2 cos α + f π 2 2 ( sin α ) 2 μ I 2 ( 1 n 3 2 ) ,
showing that it is independent of n 1 and n 2 , corresponding to the flat direction spanned by the NGB that interpolates between the π + and π fields, see Section 3. For maximizing the Lagrangian one has to take n 3 = 0 , then we have the freedom to take say n 1 = 0 and n 2 = 1 . In summary, the ground state has vanishing projection along the isospin direction while the rotations around the direction of the chemical potential leave the vacuum invariant.
We now return to the three-flavor case. In the three different meson condensed phases we expect that it suffices to take only one entry of the unit vector n nonzero. In particular, in the pion condensed phase
Σ ¯ π = e i α λ 2 = cos α sin α 0 sin α cos α 0 0 0 1 = 1 + 2 cos α 3 I + i λ 2 sin α + cos α 1 3 λ 8 ,
while in the charged kaon condensed phase
Σ ¯ K = e i α λ 5 = cos α 0 sin α 0 1 0 sin α 0 cos α = 1 + 2 cos α 3 I + cos α 1 2 3 3 λ 3 λ 8 + i λ 5 sin α ,
finally in the neutral kaon condensed phase
Σ ¯ K 0 = e i α λ 7 = 1 0 0 0 cos α sin α 0 sin α cos α = 1 + 2 cos α 3 I + 1 cos α 2 3 3 λ 3 + λ 8 + i λ 7 sin α ,
where α assumes different values in the three different phases.
To allow the transition between the three phases one can consider the more general ansatz
Σ ¯ = e i γ λ 2 e i β λ 7 e i α λ 2 e i β λ 7 e i γ λ 2 ,
with α , β , γ [ 0 , π / 2 ] are three different angles. Once again, the normal phase corresponds to α = 0 and it is insensitive to the values of β and γ . The pion condensed phase corresponds to β = 0 , γ = 0 , the K + condensed phase to β = π / 2 , γ = 0 , the K 0 condensed phase to β = π / 2 , γ = π / 2 . Any value of β and γ different from the above ones indicates a phase with simultaneous condensation of two or more meson fields. We have argued in Section 3 that this is not the case. This has been explicitly shown in χ PT for the π + and K + condensation phases [11]. In this case we can simplify the ground state ansatz by taking γ = 0 in Equation (62), obtaining the same ansatz of [11]. Any value β ( 0 , π / 2 ) corresponds to the simultaneous π + - K + condensation. The analysis of [11] shows that at the phase transition the angle β discontinuously jumps from 0 to π / 2 resulting in a first order phase transition (as depicted in Figure 5): no simultaneous condensation is possible.
Below we summarize the main properties of the possible phases obtained by the χ PT analysis. The phase diagram is in agreement with Figure 5, indeed it was first derived by χ PT in [11]. Given the symmetry of the phase diagram we focus on the μ I > 0 and μ S > 0 part of Figure 5 characterized by the π + and K + condensates. The red solid line is exactly determined, while the χ PT prediction for the dashed blue line is
μ S = m π 2 + ( m π 2 μ I 2 ) 2 + 4 m K 2 μ I 2 2 μ I .
We also report below the relevant thermodynamic quantities. Once Σ ¯ is determined, the pressure is obtained from
P = L static ( Σ ¯ ) ,
and then the number densities and the energy density follow from the thermodynamic relations
n I , S = P μ I , S , ϵ = μ I n I + μ S n S P .
Finally one can obtain the equation of state (EoS) by appropriately expressing the chemical potentials as a function of the pressure.
  • The normal phase is favored for
    μ I < m π , μ S < m K 1 2 μ I ,
    with the trivial vev Σ ¯ N = diag ( 1 , 1 , 1 ) . The nonvanishing condensates are the three chiral condensates, see Equation (7),
    σ u = σ d = σ s σ 0 ,
    where the subscript indicates the quark flavor and σ 0 is the value of the chiral condensate in vacuum. In the normal phase the values of the chiral condensates are not affected by the chemical potentials. The pressure is given by
    P = f π 2 m π 2 for N f = 2 , f π 2 m π 2 1 2 + m K 2 m π 2 for N f = 3 ,
    and thus the isospin and strangeness number densities vanish, that is
    n I = n S = 0 .
  • The π + condensed phase is favored for
    μ I > m π , μ S < m π 2 + ( m π 2 μ I 2 ) 2 + 4 m K 2 μ I 2 2 μ I ,
    resulting in the vacuum in Equation (59) with
    cos α π = m π μ I 2 ,
    determined by maximizing the static Lagrangian in Equation (58). The condensates are given by
    σ u = σ d = σ 0 cos α π σ s = σ 0 π + = σ 0 sin α π K + = 0 ,
    and the pressure produced by the condensation of pions is given by [10,11]
    P = f π 2 μ I 2 2 1 m π 2 μ I 2 2 ,
    where the normal phase pressure has been subtracted. This expression is valid for both N f = 2 and N f = 3 ; it is of course insensitive to the kaon mass and the strange quark chemical potential. The number densities are
    n I = f π 2 μ I 1 m π 4 μ I 4 n s = 0 ,
    and the O ( p 2 ) equation of state [27] is
    ϵ ( P ) = P + 2 P ( 2 f π 2 m π 2 + P ) .
  • The K + condensed phase is favored for
    μ S > m K 1 2 μ I , μ S > m π 2 + ( m π 2 μ I 2 ) 2 + 4 m K 2 μ I 2 2 μ I ,
    resulting in the vacuum in Equation (60) with
    cos α K = m K μ K 2 ,
    where μ K = μ I / 2 + μ S is the relevant combination of chemical potentials, because K + has isospin 1 / 2 and strangeness 1. The condensates are given by
    σ u = σ s = σ 0 cos α K σ d = σ 0 π + = 0 K + = σ 0 sin α K ,
    and the normalized pressure by
    P = f π 2 μ K 2 2 1 m K 2 μ K 2 2 ,
    which consents to obtain the number densities
    n I = f π 2 μ K 2 1 m K 4 μ K 4 n s = 2 n I .
    The EoS is
    ϵ ( P ) = P + 2 P ( 2 f π 2 m K 2 + P ) .
As we shall see in Section 4.3, these χ PT results are in good agreement with the NJL and LQCD results close to the second order phase transitions.
Let us briefly comment on the difference between the charged meson condensed phases and the neutral kaon condensed phase. The expression of the various thermodynamic quantities in the neutral kaon condensed phase can be obtained from those of the charged kaon condensed phase by replacing μ I μ I , due to the fact that K 0 has isospin 1 / 2 . The relevant difference, as we have already noted is Section 3, is that the K 0 and the K ¯ 0 condensed phases are superfluid, while the charged meson condensed phases are superconductors. In the latter case one can determine the screening masses of the electromagnetic field by gauging the U ( 1 ) Q subgroup of the chiral group, see Equation (30), resulting in the Debye and Meissner screening masses [26]
M D 2 = M M 2 = f π 2 e 2 ( sin α ) 2 ,
where the Debye mass is related to the electric charge susceptibility, see for instance [112], while the nonvanishing value of the Meissner mass implies that the system is a superconductor.

4.1.2. Low-Energy Excitations

Once the ground state has been identified, one can determine the low-energy fluctuations by an appropriate expansion. This is quite useful also because it allows us to identify the NGB. We briefly illustrate the procedure for the two-flavor case. A useful parameterization is
Σ = cos ρ + i φ ^ · σ sin ρ ,
where the radial field, ρ , and the unit vector field, φ ^ , encode in a nontrivial way the three pion fields. By this parameterization the LO χ PT Lagrangian takes the form obtained in [29]
L = f π 2 2 μ ρ μ ρ + sin 2 ρ μ φ ^ i μ φ ^ i 2 m π γ sin 2 ρ ϵ 3 i k φ ^ i 0 φ ^ k V ( ρ ) ,
where
V ( ρ ) = f π 2 m π 2 cos ρ + γ 2 2 sin 2 ρ ,
is the potential and γ = μ I / m π is the control parameter. From the ground state analysis we know that the pion condensed phase is favored for | γ | > 1 , and the minimum of the potential is attained for the radial field vev, ρ ¯ = α π , see Equation (71).
The low-energy radial and angular excitations can now be introduced as follows [26,29,30]
ρ = α π + χ , φ ^ = ( φ 1 , φ 2 ) ,
where we have neglected the fluctuation of the φ 3 π 0 field, because it decouples. We can parameterize the angular field by
φ ^ = ( cos θ , sin θ ) ,
where θ is the Bogolyubov mode, and rescaling the radial field as χ χ / f π and the Bogolyubov mode as θ θ / f π sin ρ ¯ , one obtains the quadratic Lagrangian
L = 1 2 μ χ μ χ 1 2 m χ 2 χ 2 + 1 2 μ α μ α g χ 0 α ,
where m χ = m π γ sin α π is the mass of the radial field fluctuations and g = 2 μ I is the coupling between the oscillations of the radial and the angular fields. The mass of the radial mode vanishes at the phase transition to the normal phase because it is a second order phase transition. The Bogolyubov field seems to propagate at the speed of light, but integrating out the radial fluctuations one obtains the actual NGB with a phonon-like dispersion law
E ph = c s p ,
thus propagating at the sound speed [10],
c s = P n I = γ 4 1 γ 4 + 3 .
Alternatively, by diagonalizing the quadratic Lagrangian one obtains the dispersion laws
E ± = p 2 + m eff 2 2 ± m eff 2 2 2 + g 2 p 2 ,
where the low momentum expansion of the E field coincides with the NGB dispersion law and the other mode with mass
m eff 2 = m χ 2 + g 2 = m π 2 γ 4 + 3 γ 2 ,
is the rotated radial mode. In conclusion, the low-energy modes correspond to a NGB with dispersion law in Equation (89) and to a radial mode with mass m eff .

4.2. The Nambu-Jona Lasinio Model

The meson condensed phases can also be studied by a modeling of the strong interaction by contact interaction terms [34,35,36,37,38,39,40,41,45,46,47,50,52,55], see [54] for a brief recent review. These models stem from the original work by Nambu and Jona Lasinio [113,114,115] of a pre-QCD Lagrangian for the description of the strong interaction by contact interaction terms:
L = ψ ¯ i γ μ μ M ψ + G [ ( ψ ¯ ψ ) 2 + ( ψ ¯ i γ 5 σ ψ ) 2 ] ,
where ψ is the two-nucleon isodoublet, M is the pertinent mass matrix and G is a dimensional coupling. The interaction preserves the global symmetry group G, see Equation (30), for a proper description of hadronic matter. The nucleons emerge as quasiparticle states and the spontaneously breaking of the chiral symmetry leads to the appearance of mesons. The model is based on an analogy with the BCS theory of superconductivity describing the electron interaction by means of a local interaction term with no gauge fields. It is not completely specified until a regularization scheme is provided and the value of the coupling constant is fixed.
In the modern view [35,57,116,117,118], the model describes quark matter with an effective contact interaction term that preserve the chiral symmetries of QCD. The spinor ψ in Equation (93) now represents the quark fields and M the corresponding mass matrix, see Equation (3). The NJL model (eventually supplemented by Polyakov loop terms) has been applied to study the entire QCD phase diagram in Figure 1. The major phenomenological shortcoming of the quark NJL model is that it does not provide a confinement mechanism, indeed it has no gauge dynamics. Moreover, the presence of a dimension 6 operator requires an ultra-violet regularization scheme [57], which in the most used approximations is a hard cutoff at the Λ 1 GeV scale or a form factor of the form [119]
F ( p 2 ) = Λ 2 p 2 + Λ 2 ,
to mimic the asymptotic freedom property of QCD. The coupling constant and the bare quark masses are then fixed to reproduce the low energy physics. Typical values of these quantities are
G Λ 2 6 m 1.5 MeV m s 50 MeV ,
see however [57,118]. The NJL model is a useful tool for a qualitative and semiquantitative exploration of the properties of hadronic matter, however the obtained results depend on the choice of these parameters and on the regularization scheme employed. Unfortunately, one cannot systematically improve the model because no expansion parameter can be identified. Despite these limitations, the NJL Lagrangian is Lorentz invariant, with the chiral symmetry realized and spontaneously broken exactly as it is expected to happen in QCD: by a chiral condensate. Moreover, the chiral symmetry can be explicitly broken by the inclusion of small current quark masses.
There is a certain degree of uncertainty in the form of the NJL Lagrangian. In the two-flavor case most of the authors retain the form in Equation (93), although different chirally symmetric interactions can be written. This increases the number of phenomenological parameters that have to be fixed. Following [57,118,120], the NJL Lagrangian can be generalized to
L = ψ ¯ i γ μ μ M ψ + L 1 + L 2 ,
where the two interaction terms are
L 1 = G 1 [ ( ψ ¯ ψ ) 2 + ( ψ ¯ i γ 5 σ ψ ) 2 + ( ψ ¯ σ ψ ) 2 + ( ψ ¯ i γ 5 ψ ) 2 ] L 2 = G 2 [ ( ψ ¯ ψ ) 2 + ( ψ ¯ i γ 5 σ ψ ) 2 ( ψ ¯ σ ψ ) 2 ( ψ ¯ i γ 5 ψ ) 2 ] .
The first term preserves the U ( 1 ) A symmetry while the second term explicitly breaks it. Whether or not the latter is comparable with the first depends on nonperturbative effects, indeed the U ( 1 ) A breaking term is supposed to describe the interaction mediated by some instanton configurations of the gauge fields. Considering different values of the coupling constants G 1 and G 2 implies distinct values of the chiral condensates of different flavors and a different phase diagram. This issue emerges also in the three-flavor case, where the NJL Lagrangian takes a slightly different form [121],
L = ψ ¯ i γ μ μ M ψ + G a = 0 8 [ ( ψ ¯ λ a ψ ) 2 + ( ψ ¯ i γ 5 λ a ψ ) 2 ] K [ det ψ ¯ ( 1 + γ 5 ) ψ + det ψ ¯ ( 1 γ 5 ) ψ ] ,
where now ψ T = ( u , d , s ) and G and K are the two coupling constants analogous of G 1 and G 2 ; the determinant term removes the U ( 1 ) A symmetry. Considering the U ( 1 ) A symmetric Lagrangian with K = 0 , the authors of [36] find the same transition hexagon (solid red line) reported in Figure 5, however the up and down chiral condensates split and at large values of the chemical potentials they find phase transition lines not present in Figure 5.
In the following we will focus on the traditional two-flavor NJL model, with G 1 = G 2 = G / 2 corresponding to the Lagrangian in Equation (93), and thus maximally violated U ( 1 ) A symmetry. The presence of a medium can be described by a covariant derivative analogous to Equation (4), which takes into account the baryon and isospin chemical potentials. With the NJL model one can explore the entire QCD phase diagram (with the limitations discussed above), a clear advantage with respect to the χ PT approach which can hardly investigate the effect of the baryon chemical potential. To obtain the properties of the vacuum and of the low-energy excitations one can perform a Hubbard-Stratanovich transformation introducing the collective boson variables
σ f ( x ) = 4 G Λ ( ψ ¯ f ψ f ) π a = 2 G Λ ( ψ ¯ γ 5 i λ a ψ + h . c . ) ,
corresponding to scalar and pseudoscalar fields. Their expectation values are determined by minimizing the one loop effective potential, or equivalently, by solving the coupled gap equations. By this analysis it has been confirmed that the pion condensed phase sets in at μ I = m π , see [38]. Moreover, it has been determined the dependence of the condensates on the chemical potentials and on the temperature.
At vanishing temperature, in the two flavor case the grand potential has the particularly simple expression
Ω = G ( σ 2 + π 2 ) 3 2 π 2 0 Λ d k k 2 ( E + + E ) ,
where we have used the hard cutoff procedure and the quark quasiparticle dispersion laws are given by
E q , ± = E k ± μ I 2 2 + 4 G 2 π 2 ,
where E k = k 2 + m ¯ 2 with m ¯ = m 2 G σ the effective quark mass. Therefore, the effect of the chiral condensate is a shift of the quark masses [57], while the pion condensate opens a gap between the quasiparticle dispersion laws. This is the typical effect of condensation on the quasiparticle spectrum, as it indicates the formation of correlated pairs of fermions. In the pion condensed phase it costs additional energy to produce quasiparticle fermionic excitations because of quark-antiquark pairs. Clearly, this picture of meson condensation cannot be directly compared with the χ PT results of the previous section, because the considered degrees of freedom are different. However, as we shall discuss below, the values of the pion and chiral condensates can be compared, as well as various thermodynamic quantities. Moreover, the three flavor NJL model phase diagram obtained in [38] is in agreement with the theory group expectation in Figure 5 and quantitatively very similar to that obtained in χ PT.
In the NJL model it is possible to include an electron (or positron) background to neutralize the pion electric charge. When requiring the electrical neutrality [39,40], the NJL models tend to disfavor the appearance of the pion condensed phase [45,47]. At the physical point, corresponding to a neutral configuration in hydrostatic equilibrium, the pions do not condense [45].

4.3. Comparison with Lattice QCD

The LQCD simulations are numerical implementation of the QCD action on a discretized grid; the relevant physical results are then obtained performing the limit to the continuum. These simulations can lead to the precise determination of many hadronic quantities and can provide numerical evidence for conjectured properties of strongly interacting matter. The LQCD simulations have been very successfully used for simulating hadronic matter in vacuum, but dealing with in medium effects poses a series of problems. The most important one is that the LQCD simulations at finite baryonic density are hampered by the so-called sign problem. Very briefly, in the LQCD simulations with dynamical quarks the Dirac degrees of freedom are typically integrated out, see for instance [122,123], resulting in a partition function that can be written as the euclidean path integral
Z d A μ e S ( A μ ) det Δ D ,
where S is the euclidean action, A μ are the gauge fields and det Δ D is the determinant of the Dirac operators. The standard Monte Carlo simulations are based on importance sampling of the possible gauge configurations. This procedure works if det Δ D > 0 , that is with a real and positive Euclidean path integral measure. At nonvanishing baryonic density the LQCD numerical technique becomes problematic because the Dirac determinant is complex. Although continuous progress for facing this problem has been reported over the years, see for example [124,125,126,127,128], as of yet it is not a feasible tool for exploring the QCD phase diagram at large μ B and, in particular, the transition from the confined phase to the CSC phase in Figure 1. See however [68,129,130,131] for different LQCD approaches to the region with nonvanishing baryonic and isospin chemical potentials.
Since it is hard to manage baryons in LQCD simulations, people decided to ignore baryons. This poses the LQCD simulations outside the beta-equilibrated sheet, as discussed in Section 1, to explore a part of the QCD phase diagram where the outcomes of the numerical simulations can be compared with different methods, in particular with the χ PT and the NJL results. The key point is indeed that the LQCD simulations at nonvanishing isospin chemical potential and zero baryonic density are not affected by the sign problem [68]. This does not mean that this direction is without obstacles: the realization of multi-hadron systems in LQCD is an extremely challenging problem, see [132] for a review. There is a wealth of LQCD results for the pion condensation [70,71,72,73,76,77,79,82] while there has been little progress on kaons [74,75]. As we will see, the LQCD results on pion condensation are reliable only for μ I 2 m π , but these simulations are steadily improving and becoming more accurate, even with physical quark masses and external magnetic fields [77,79,82]
The LQCD simulations can be performed in the canonical or in the grand canonical ensembles. In the grand canonical simulation one discretizes on a lattice the actual QCD Lagrangian in Equation (2) with the isospin chemical potential as external source. In this approach the strangeness density has to be zero, as the strangeness chemical potential makes the measure complex. Moreover, the QCD Lagrangian has to be supplemented with a pionic source [69,70]
L λ = i λ ψ ¯ γ 5 σ 2 ψ ,
to trigger the breaking of the U ( 1 ) I symmetry and to stabilize the numerical simulations. Since the λ term explicitly breaks the U ( 1 ) I symmetry, in the pion condensed phase there is a pseudo-NGB with vanishing mass in the λ 0 limit. Therefore, the physical interesting results are obtained doing both the continuum and the λ 0 limits. The first quenched LQCD simulations [69] already reported the expected behavior of the condensates; these results were soon improved considering N f = 2 dynamical quarks in [70]. To obtain more precise results one has to consider that both the pion and the chiral condensates depend in a rather non-trivial way on λ . Moreover, the first simulations employed a large pion mass. Progress with respect to these aspects has been reported in [133] where the λ 0 limit has been tackled by a reweighing technique in simulations with 2 light flavors and a heavy strange quark at the physical pion mass.
In Figure 7 we compare the pion and chiral condensates obtained by the LO χ PT, see Section 4.1, by the two-flavor NJL model, see Section 4.2, and by the LQCD simulations of [133].
The chiral condensate obtained with the three approaches agree for any considered (or available) value of μ I . On the other hand, the pion condensates deviate at μ I 2 m π . In χ PT the the two condensates obey the relation
σ 2 + π + 2 = σ 0 2 ,
and therefore the pion condensate quickly saturates at μ I 2 m π . Both the NJL and the LQCD results for the pion condensate indicate that it exceeds σ 0 and it does not saturate at μ I 2 m π . This seems a robust result, although for larger values of μ I both the NJL and the LQCD approaches become problematic. The NJL results show a non-monotonic behavior due to the hard cutoff Λ , which serves to mimic the asymptotic behavior of QCD, but that also signals the scale at which the NJL results are not under control. Similarly, the LQCD simulations feel the finite size lattice effects, indeed for large value of μ I saturation effects become important, see for example [70]. Anyway, one can certainly regard this comparison as successful, in the sense that in the range μ I 2 m π , where all the three approaches are supposed to work, they give very similar results.
We now turn to the canonical approach [73,74,75,76]. The canonical LQCD simulations can explore quark matter at nonvanishing isospin and strangeness density, while grand-canonical LQCD simulations can only deal with finite isospin density. In the canonical LQCD simulations the isospin and strangeness density are fixed and the corresponding chemical potentials are determined by thermodynamic relations. The description of mesons in the canonical LQCD simulation is attained by the introduction of external sources with a fixed isospin or strangeness charge. In these simulations the calculation of the meson field correlator requires the computation of a large number of Wick contraction of the quark fields on the lattice, leading to time consuming and expensive calculations. This is the main limitation of the canonical lattice simulations. Various different algorithms for reducing theses costs have been developed in [76], resulting in the simulation of up to 72 π + in configurations with spatial extents L 2 , 2.5 and 3 fm, resulting in isospin chemical potentials up to 4.5 m π [76].
In Figure 8 we compare the energy density obtained with the canonical LQCD simulations with that of the χ PT and the NJL approaches.
More precisely, in this figure it is shown the normalized energy density ϵ / ϵ SB , where ϵ S B = 9 μ I 4 / ( 4 π 2 ) is the Stefan-Boltzmann limit, as a function of the normalized isospin chemical potential, μ I / m π . The dots with error bars are the results of [76] obtained with three different spatial volumes L 3 . With increasing μ I the error bars increase, signaling that lattice simulations cease to be reliable at μ I 2 m π , as in the grand-canonical LQCD simulations discussed above. In this regime, the normalized energy well agrees with the χ PT results (dashed black line), and with the NJL results (solid blue line). Both the χ PT and the NJL curves perfectly capture the peak structure at low μ I , while they begin to depart from the LQCD results at μ I 2 m π . The χ PT and the NJL peak positions are respectively at
μ I , χ PT peak = 13 2 1 / 2 m π 1.27 m π ,
μ I , NJL peak 1.27 m π ,
where the χ PT results are independent of f π , while the NJL results are not very sensitive to the parameter set used. The LQCD results of [76] are peaked at
μ I , LQCD peak = { 1.20 , 1.25 , 1.275 } m π ,
where the different values are obtained for lattice sides L = { 16 , 20 , 24 } , respectively; the continuum-linearly-extrapolated peak is at μ I , LQCD peak 1.30 m π . Therefore, also the canonical LQCD simulations are in agreement with the χ PT and the NJL results for μ I 2 m π .

Nonvanishing Temperature

Given the successful comparison of the χ PT, NJL and LQCD approaches at T = 0 , one may expect a similar agreement at small temperatures. As we will see, the agreement between the three methods at T 0 is much worse. Herein we report on the investigation of the phase diagram at T 0 and μ I 0 comparing the results for the transition lines at μ B = 0 separating the pion condensed phase from the normal phase at low T, and between the pion condensed phase and the quark-gluon plasma, at high T. We show in Figure 9 the results obtained with the different approaches. The LQCD simulations of [133] indicate that at μ B = 0 there is in the μ I , T plane a chiral crossover line (shaded blue region) joining the points ( 0 , 160 ) to the (pseudo) tricritical point ( 140 , 151 ) (orange dot with error bars). At the tricritical point the chiral crossover line joins the second order phase transition line (shaded green region). The LQCD results for the second order phase transition are almost insensitive to the temperature for T 150 MeV, then the phase transition line becomes strongly temperature dependent, with a sort of “T-like” phase diagram shape. The mean-field NJL second order phase transition (solid blue line) [38] shows a behavior similar to that of the LQCD simulations for T 100 MeV, then for higher temperatures the NJL results show a more pronounced temperature dependence. Eventually, the NJL critical curve saturates with a critical temperature that is not sensitive to the isospin chemical potential for 500 MeV < μ I < 1 GeV (not shown in the figure). The analytic χ PT temperature dependence of the second order phase transition has been obtained in [22]
μ I ( T ) = m π + 1 4 f π 2 m π 3 T 3 2 π 3 ζ 3 2 ,
and is reported in Figure 9 with a dashed black line. The T 3 / 2 behavior does not agree with the LQCD nor with the NJL results. The χ PT results of [24] indicate an even stronger temperature dependence. These results are somehow surprising, as one would expect χ PT to work up to T 100 MeV while Figure 9 shows that it is inconsistent with the LQCD low temperature behavior.
Given the rather precise LQCD data, one should understand what are the origins of the discrepancies. The NJL results have been obtained by a hard cutoff scheme, maybe one can relax this requirement by a Pauli-Villars regularization scheme or by a form factor, as in (94), that does not completely eliminate the hard scale contribution. The improved χ PT results of [24] do not match the LQCD behavior at low μ I but indicate a critical temperature that is independent of the chemical potential for μ I > T , which is in agreement with both the LQCD and NJL simulations. One should certainly try to understand what is the χ PT missing ingredient at lower temperatures.
Quite remarkably, the LQCD simulations are now tackling the color deconfinement transition as a function of the isospin chemical potential. Color deconfinement can be characterized by the behavior of the so-called Polyakov loop, see for instance [4,134],
P ( r ) = 1 3 Tr Π τ = 1 N τ U τ ( r ) ,
at large lattice spacing. The expectation value of the Polyakov loop is related to the correlation function between two static heavy quarks, therefore it is a measure of the strength of the color interaction. It has the important property to vanish in the color confined phase of pure gauge theories [4]. A first inspection of the Polyakov loop dependence at nonvanishing μ I and T has been done in [133], by a 24 3 × 6 lattice. Various lines of constant values of the Polyakov loop have been obtained to infer the position of the deconfinement critical temperature. In Figure 9 we report the line of [133] (shaded red line) corresponding to P = 1 , which in their notation can be taken as indicating the color deconfinement transition. This line partially overlaps with the chiral crossover line, but then it starts to bend inside the pion condensed phase. These preliminary results should be tested with different lattice spacings. Quite interestingly, the results of [133] indicate that the deconfinement line is quite insensitive to the presence of the pion condensate, or equivalently, to the melting of the chiral condensate at large μ I . We expect that for sufficiently large μ I the BEC pion condensate turns in a BCS condensate. In this case the deconfined quarks should form quark-antiquark Cooper pairs, pretty much as in the CSC phase, but with an important difference: the BCS pairs in this case can be color singlets. Therefore, in this case there is no need to have color deconfinement nor any phase transition at all. Quite generally, indeed, there is no phase transition between the BEC and the BCS phases [91]. Moreover, it would be interesting to see the behavior of the energy density of the system as a function of T for a fixed μ I , as we expect that the deconfinement phase transition should induce a rapid increase of the energy density due to the liberation of the quark and gluon degrees of freedom.

5. Conclusions

We have briefly reviewed the meson condensation phenomenon happening when the isospin or the strange chemical potentials exceed a critical value. We have clarified that it is unclear whether or not these phases can be realized in Nature. In vacuum all mesons are unstable, therefore a stable meson can only exist in a dense medium, as in compact stars. In the core of these stellar objects the large number of electrons may stabilize the π , however the problem is that with increasing density other particles compete with π to share the excess electron charge. By a simple noninteracting model we have seen how the electron negative charge is drained off into Σ states favoring the strangeness production. The strong interactions can modify this picture, but it is unclear whether they favor or disfavor the appearance of stable pions. Quite recently, it has been proposed that pion stars consisting of pions and charged leptons may exist [29,135,136]. The astrophysical observation of this exotic star would certainly be a smoking gun of a macroscopic coherent state of pions.
Although it is unclear whether the meson condensed phases are realized in compact stars or in any other physical setting, they are interesting by themselves. The reason is that they allow us to explore the properties of QCD in a regime in which various methods overlap. In particular, the χ PT, the NJL and the LQCD approaches give similar results at vanishing temperature for μ I 2 m π . The μ I μ S phase diagram in Figure 5 finds χ PT and NJL in excellent agreement, and the LQCD numerical results have confirmed that the phase transition between the normal phase and the pion condensation phase is of the second order. Unfortunately, the entire μ I μ S phase diagram has not been completely explored by canonical LQCD simulations; such simulations could allow to figure out whether mixed phases are realized.
These findings allow us to improve the first version of the QCD phase diagram shown in Figure 1. That diagram was based on naive arguments on the strong interaction. We now draw in Figure 10 the QCD phase diagram in which we have fed the acquired knowledge.
The solid thick lines correspond to the second order phase transitions that are most tenable. The LQCD simulations give numerical support to the second order phase transition line at μ I = m π , which should be temperature independent up to T 150 MeV, see also Figure 9. The dotted lines indicate a chiral and (quite probably) deconfinement crossover. Since there is a chiral crossover at small μ B and μ I = 0 [125,126,127,128] and since we have seen in the previous section that there is a chiral crossover at μ B = 0 and μ I m π , the most simple possibility is that there is an almost temperature independent chiral crossover region at T T c (dashed area). There is indeed growing evidence that the chiral crossover extends in the μ B , μ I plane at an almost constant temperature T T c [35,38,63,137,138]. The general result of these works is that the transition temperature smoothly decreases with μ I (or μ S ), but the order of the phase transition is hard to establish. From Figure 10 it seems like the hadron gas occupies a first octant sphere with the edges cut by a vertical and an horizontal plane. However, there are still uncertain transitions, marked with a thin blue line. We have added a phase transition line between the hadron gas phase and the color superconducting phase, although there is no experimental data nor any LQCD simulation that supports it and may as well be a smooth crossover [139], see also [5,6]. We have also assumed that the hadron gas is separated from the pion condensed phase and/or the color superconducting phase by a transition (thin blue) line extending in the T = 0 plane, which is just a guess. We have not shown in Figure 9 any transition line between the color superconducting phase and the pion condensed phase. This happens in a region where three different quark condensates compete and it is not at all obvious that the phase diagram has a simple form.

Funding

This research received no external funding.

Acknowledgments

I would like to thank Gergely Endrodi and Bastian Brandt for sharing their LQCD data, Mark Alford and Gergely Endrodi for very useful suggestions and Jens Oluf Andersen for discussion. I thank the University of Bari and INFN for the support during the completion of this work.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Cabibbo, N.; Parisi, G. Exponential Hadronic Spectrum and Quark Liberation. Phys. Lett. B 1975, 59, 67–69. [Google Scholar] [CrossRef]
  2. Gyulassy, M. The QGP Discovered at RHIC. In Structure and Dynamics of Elementary Matter, Proceedings of the NATO Advanced Study Institute, Camyuva-Kemer, Turkey, 22 September–2 October 2003; Springer: Berlin/Heidelberg, Germany, 2004; pp. 159–182. [Google Scholar]
  3. Shuryak, E. Physics of Strongly coupled Quark-Gluon Plasma. Prog. Part. Nucl. Phys. 2009, 62, 48–101. [Google Scholar] [CrossRef]
  4. Satz, H. Extreme states of matter in strong interaction physics: An introduction. Lect. Notes Phys. 2012, 841, 1–239. [Google Scholar] [CrossRef]
  5. Rajagopal, K.; Wilczek, F. The Condensed matter physics of QCD. Front. Part. Phys. 2001, 3, 2061–2151. [Google Scholar]
  6. Alford, M.G.; Schmitt, A.; Rajagopal, K.; Schafer, T. Color superconductivity in dense quark matter. Rev. Mod. Phys. 2008, 80, 1455–1515. [Google Scholar] [CrossRef] [Green Version]
  7. Anglani, R.; Casalbuoni, R.; Ciminale, M.; Ippolito, N.; Gatto, R.; Mannarelli, M.; Ruggieri, M. Crystalline color superconductors. Rev. Mod. Phys. 2014, 86, 509–561. [Google Scholar] [CrossRef] [Green Version]
  8. Borsanyi, S.; Fodor, Z.; Hoelbling, C.; Katz, S.D.; Krieg, S.; Ratti, C.; Szabo, K.K. Is there still any Tc mystery in lattice QCD? Results with physical masses in the continuum limit III. JHEP 2010, 09, 073. [Google Scholar] [CrossRef]
  9. Bazavov, A.; Bhattacharya, T.; Cheng, M.; DeTar, C.; Ding, H.-T.; Gottlieb, S.; Gupta, R.; Hegde, P.; Heller, U.M.; Laermann, E.; et al. The chiral and deconfinement aspects of the QCD transition. Phys. Rev. D 2012, 85, 054503. [Google Scholar] [CrossRef]
  10. Son, D.; Stephanov, M.A. QCD at finite isospin density. Phys. Rev. Lett. 2001, 86, 592–595. [Google Scholar] [CrossRef]
  11. Kogut, J.; Toublan, D. QCD at small nonzero quark chemical potentials. Phys. Rev. D 2001, 64, 034007. [Google Scholar] [CrossRef] [Green Version]
  12. Migdal, A.B. Stability of vacuum and limiting fields. Zh. Eksp. Teor. Fiz. 1971, 61, 2209–2224. [Google Scholar] [CrossRef]
  13. Migdal, A.B. Vacuum Stability and Limiting Fields. Soviet Phys. Uspekhi 1972, 14, 813. [Google Scholar] [CrossRef]
  14. Sawyer, R.F. Condensed pi- phase in neutron star matter. Phys. Rev. Lett. 1972, 29, 382–385. [Google Scholar] [CrossRef]
  15. Scalapino, D.J. Pi-condensate in dense nuclear matter. Phys. Rev. Lett. 1972, 29, 386–388. [Google Scholar] [CrossRef]
  16. Kogut, J.; Manassah, J.T. πcondensation and neutron star cooling. Phys. Lett. A 1972, 41, 129–131. [Google Scholar] [CrossRef]
  17. Migdal, A.B. Pi condensation in nuclear matter. Phys. Rev. Lett. 1973, 31, 257–260. [Google Scholar] [CrossRef]
  18. Baym, G.; Campbell, D.K. Chiral Symmetry and Pion Condensation. In Mesons in Nuclei; Rho, M., Wilkinson, D., Eds.; North Holland Pub. Co.: Amsterdam, The Netherlands, 1978; p. 1031. [Google Scholar]
  19. Kaplan, D.B.; Nelson, A.E. Strange Goings on in Dense Nucleonic Matter. Phys. Lett. B 1986, 175, 57–63. [Google Scholar] [CrossRef]
  20. Dominguez, C.; Loewe, M.; Rojas, J. Pion and nucleon thermal widths in the linear sigma model. Phys. Lett. B 1994, 320, 377–380. [Google Scholar] [CrossRef]
  21. Birse, M.C.; Cohen, T.D.; McGovern, J.A. Phases of QCD with nonvanishing isospin density. Phys. Lett. B 2001, 516, 27–32. [Google Scholar] [CrossRef]
  22. Splittorff, K.; Toublan, D.; Verbaarschot, J.J.M. Thermodynamics of chiral symmetry at low densities. Nucl. Phys. B 2002, 639, 524–548. [Google Scholar] [CrossRef] [Green Version]
  23. Loewe, M.; Villavicencio, C. Thermal pions at finite isospin chemical potential. Phys. Rev. D 2003, 67, 074034. [Google Scholar] [CrossRef]
  24. Loewe, M.; Villavicencio, C. Thermal pion masses in the second phase: |mu(I)|>m(pi). Phys. Rev. D 2004, 70, 074005. [Google Scholar] [CrossRef]
  25. Loewe, M.; Villavicencio, C. Pion stability in a hot dense media. arXiv 2011, arXiv:1107.3859. [Google Scholar]
  26. Mammarella, A.; Mannarelli, M. Intriguing aspects of meson condensation. Phys. Rev. D 2015, 92, 085025. [Google Scholar] [CrossRef]
  27. Carignano, S.; Mammarella, A.; Mannarelli, M. Equation of state of imbalanced cold matter from chiral perturbation theory. Phys. Rev. D 2016, 93, 051503. [Google Scholar] [CrossRef] [Green Version]
  28. Loewe, M.; Raya, A.; Villavicencio, C. Metastable Pions in Dense Media. Phys. Rev. D 2016, 95, 096013. [Google Scholar] [CrossRef]
  29. Carignano, S.; Lepori, L.; Mammarella, A.; Mannarelli, M.; Pagliaroli, G. Scrutinizing the pion condensed phase. Eur. Phys. J. A 2017, 53, 35. [Google Scholar] [CrossRef] [Green Version]
  30. Lepori, L.; Mannarelli, M. Multicomponent meson superfluids in chiral perturbation theory. Phys. Rev. D 2019, 99, 096011. [Google Scholar] [CrossRef] [Green Version]
  31. Adhikari, P.; Andersen, J.O.; Kneschke, P. QCD at finite isospin density: chiral perturbation theory confronts lattice data. arXiv 2019, arXiv:1909.01131. [Google Scholar]
  32. Tawfik, A.N.; Diab, A.M.; Ghoneim, M.T.; Anwer, H. SU(3) Polyakov Linear-Sigma Model With Finite Isospin Asymmetry: QCD Phase Diagram. arXiv 2019, arXiv:1904.09890. [Google Scholar]
  33. Mishustin, I.N.; Anchishkin, D.V.; Satarov, L.M.; Stashko, O.S.; Stoecker, H. Condensation of interacting scalar bosons at finite temperatures. Phys. Rev. C 2019, 100, 022201. [Google Scholar] [CrossRef] [Green Version]
  34. Barducci, A.; Casalbuoni, R.; De Curtis, S.; Gatto, R.; Pettini, G. Pion Decay Constant at Finite Temperature and Density. Phys. Rev. D 1990, 42, 1757–1763. [Google Scholar] [CrossRef] [PubMed]
  35. Toublan, D.; Kogut, J.B. Isospin chemical potential and the QCD phase diagram at nonzero temperature and baryon chemical potential. Phys. Lett. B 2003, 564, 212–216. [Google Scholar] [CrossRef] [Green Version]
  36. Barducci, A.; Casalbuoni, R.; Pettini, G.; Ravagli, L. A Calculation of the QCD phase diagram at finite temperature, and baryon and isospin chemical potentials. Phys. Rev. D 2004, 69, 096004. [Google Scholar] [CrossRef]
  37. Barducci, A.; Casalbuoni, R.; Pettini, G.; Ravagli, L. Pion and kaon condensation in a 3-flavor NJL model. Phys. Rev. D 2005, 71, 016011. [Google Scholar] [CrossRef]
  38. He, L.Y.; Jin, M.; Zhuang, P.F. Pion superfluidity and meson properties at finite isospin density. Phys. Rev. D 2005, 71, 116001. [Google Scholar] [CrossRef] [Green Version]
  39. Ebert, D.; Klimenko, K.G. Gapless pion condensation in quark matter with finite baryon density. J. Phys. G 2006, 32, 599–608. [Google Scholar] [CrossRef]
  40. Ebert, D.; Klimenko, K.G. Pion condensation in electrically neutral cold matter with finite baryon density. Eur. Phys. J. C 2006, 46, 771–776. [Google Scholar] [CrossRef] [Green Version]
  41. Mukherjee, S.; Mustafa, M.G.; Ray, R. Thermodynamics of the PNJL model with nonzero baryon and isospin chemical potentials. Phys. Rev. D 2007, 75, 094015. [Google Scholar] [CrossRef]
  42. He, L.; Zhuang, P. Phase structure of Nambu-Jona-Lasinio model at finite isospin density. Phys. Lett. B 2005, 615, 93–101. [Google Scholar] [CrossRef]
  43. He, L.; Jin, M.; Zhuang, P. Pion Condensation in Baryonic Matter: from Sarma Phase to Larkin-Ovchinnikov- Fudde-Ferrell Phase. Phys. Rev. D 2006, 74, 036005. [Google Scholar] [CrossRef]
  44. Sun, G.F.; He, L.; Zhuang, P. BEC-BCS crossover in the Nambu-Jona-Lasinio model of QCD. Phys. Rev. D 2007, 75, 096004. [Google Scholar] [CrossRef] [Green Version]
  45. Andersen, J.O.; Kyllingstad, L. Pion Condensation in a two-flavor NJL model: the role of charge neutrality. J. Phys. G 2009, 37, 015003. [Google Scholar] [CrossRef]
  46. Abuki, H.; Ciminale, M.; Gatto, R.; Ippolito, N.D.; Nardulli, G.; Ruggieri, M. Electrical neutrality and pion modes in the two flavor PNJL model. Phys. Rev. D 2008, 78, 014002. [Google Scholar] [CrossRef]
  47. Abuki, H.; Anglani, R.; Gatto, R.; Pellicoro, M.; Ruggieri, M. The Fate of pion condensation in quark matter: From the chiral to the real world. Phys. Rev. D 2009, 79, 034032. [Google Scholar] [CrossRef]
  48. Mu, C.f.; He, L.y.; Liu, Y.x. Evaluating the phase diagram at finite isospin and baryon chemical potentials in the Nambu-Jona-Lasinio model. Phys. Rev. D 2010, 82, 056006. [Google Scholar] [CrossRef]
  49. Xia, T.; He, L.; Zhuang, P. Three-flavor Nambu–Jona-Lasinio model at finite isospin chemical potential. Phys. Rev. D 2013, 88, 056013. [Google Scholar] [CrossRef]
  50. Xia, T.; Zhuang, P. Quark-antiquark Scattering Phase Shift and Meson Spectral Function in Pion Superfluid. Chin. Phys. D 2014, 43, 054103. [Google Scholar] [CrossRef]
  51. Chao, J.; Huang, M.; Radzhabov, A. Charged pion condensation under parallel electromagnetic fields. arXiv 2018, arXiv:1805.00614. [Google Scholar]
  52. Khunjua, T.G.; Klimenko, K.G.; Zhokhov, R.N. Chiral imbalanced hot and dense quark matter: NJL analysis at the physical point and comparison with lattice QCD. Eur. Phys. J. C 2019, 79, 151. [Google Scholar] [CrossRef]
  53. Khunjua, T.G.; Klimenko, K.G.; Zhokhov, R.N. Dualities and inhomogeneous phases in dense quark matter with chiral and isospin imbalances in the framework of effective model. JHEP 2019, 06, 006. [Google Scholar] [CrossRef]
  54. Khunjua, T.; Klimenko, K.; Zhokhov, R. Charged Pion Condensation in Dense Quark Matter: Nambu–Jona-Lasinio Model Study. Symmetry 2019, 11, 778. [Google Scholar] [CrossRef]
  55. Avancini, S.S.; Bandyopadhyay, A.; Duarte, D.C.; Farias, R.L.S. Cold QCD at finite isospin density: confronting effective models with recent lattice data. arXiv 2019, arXiv:1907.09880. [Google Scholar]
  56. Lu, Z.Y.; Xia, C.J.; Ruggieri, M. Thermodynamics and susceptibilities of isospin imbalanced QCD matter. arXiv 2019, arXiv:1907.11497. [Google Scholar]
  57. Klevansky, S.P. The Nambu-Jona-Lasinio model of quantum chromodynamics. Rev. Mod. Phys. 1992, 64, 649–708. [Google Scholar] [CrossRef]
  58. Andersen, J.O.; Naylor, W.R.; Tranberg, A. Phase diagram of QCD in a magnetic field: A review. Rev. Mod. Phys. 2016, 88, 025001. [Google Scholar] [CrossRef]
  59. Adhikari, P.; Andersen, J.O.; Kneschke, P. On-shell parameter fixing in the quark-meson model. Phys. Rev. D 2017, 95, 036017. [Google Scholar] [CrossRef] [Green Version]
  60. Adhikari, P.; Andersen, J.O.; Kneschke, P. Pion condensation and phase diagram in the Polyakov-loop quark-meson model. Phys. Rev. D 2018, 98, 074016. [Google Scholar] [CrossRef] [Green Version]
  61. Andersen, J.O.; Kneschke, P. Chiral density wave versus pion condensation at finite density and zero temperature. Phys. Rev. D 2018, 97, 076005. [Google Scholar] [CrossRef] [Green Version]
  62. Andersen, J.O.; Adhikari, P.; Kneschke, P. Pion Condensation and QCD Phase Diagram at Finite Isospin Density. In Proceedings of the 13th Conference on Quark Confinement and the Hadron Spectrum (Confinement XIII), Maynooth, Ireland, 31 July–6 August 2018. [Google Scholar]
  63. Klein, B.; Toublan, D.; Verbaarschot, J.J.M. The QCD phase diagram at nonzero temperature, baryon and isospin chemical potentials in random matrix theory. Phys. Rev. D 2003, 68, 014009. [Google Scholar] [CrossRef]
  64. Klein, B.; Toublan, D.; Verbaarschot, J. Diquark and pion condensation in random matrix models for two color QCD. Phys. Rev. D 2005, 72, 015007. [Google Scholar] [CrossRef]
  65. Lv, M.; Li, D.; He, S. Pion condensation in a soft-wall AdS/QCD model. arXiv 2018, arXiv:1811.03828. [Google Scholar]
  66. Graf, T.; Schaffner-Bielich, J.; Fraga, E.S. Perturbative thermodynamics at nonzero isospin density for cold QCD. arXiv 2015, arXiv:1511.09457. [Google Scholar]
  67. Andersen, J.O.; Haque, N.; Mustafa, M.G.; Strickland, M. Three-loop HTLpt thermodynamics at finite temperature and isospin chemical potential. arXiv 2015, arXiv:1511.04660. [Google Scholar]
  68. Alford, M.G.; Kapustin, A.; Wilczek, F. Imaginary chemical potential and finite fermion density on the lattice. Phys. Rev. D 1999, 59, 054502. [Google Scholar] [CrossRef] [Green Version]
  69. Kogut, J.B.; Sinclair, D.K. Quenched lattice QCD at finite isospin density and related theories. Phys. Rev. D 2002, 66, 014508. [Google Scholar] [CrossRef] [Green Version]
  70. Kogut, J.B.; Sinclair, D.K. Lattice QCD at finite isospin density at zero and finite temperature. Phys. Rev. D 2002, 66, 034505. [Google Scholar] [CrossRef] [Green Version]
  71. Kogut, J.B.; Sinclair, D.K. The Finite temperature transition for 2-flavor lattice QCD at finite isospin density. Phys. Rev. D 2004, 70, 094501. [Google Scholar] [CrossRef]
  72. Beane, S.R.; Detmold, W.; Luu, T.C.; Orginos, K.; Savage, M.J.; Torok, A. Multi-Pion Systems in Lattice QCD and the Three-Pion Interaction. Phys. Rev. Lett. 2008, 100, 082004. [Google Scholar] [CrossRef]
  73. Detmold, W.; Savage, M.J.; Torok, A.; Beane, S.R.; Luu, T.C.; Orginos, K.; Parreno, A. Multi-Pion States in Lattice QCD and the Charged-Pion Condensate. Phys. Rev. D 2008, 78, 014507. [Google Scholar] [CrossRef]
  74. Detmold, W.; Orginos, K.; Savage, M.J.; Walker-Loud, A. Kaon Condensation with Lattice QCD. Phys. Rev. D 2008, 78, 054514. [Google Scholar] [CrossRef]
  75. Detmold, W.; Smigielski, B. Lattice QCD study of mixed systems of pions and kaons. Phys. Rev. D 2011, 84, 014508. [Google Scholar] [CrossRef]
  76. Detmold, W.; Orginos, K.; Shi, Z. Lattice QCD at non-zero isospin chemical potential. Phys. Rev. D 2012, 86, 054507. [Google Scholar] [CrossRef]
  77. Endrödi, G. Magnetic structure of isospin-asymmetric QCD matter in neutron stars. Phys. Rev. D 2014, 90, 094501. [Google Scholar] [CrossRef] [Green Version]
  78. Janssen, O.; Kieburg, M.; Splittorff, K.; Verbaarschot, J.J.M.; Zafeiropoulos, S. Phase Diagram of Dynamical Twisted Mass Wilson Fermions at Finite Isospin Chemical Potential. Phys. Rev. D 2016, 93, 094502. [Google Scholar] [CrossRef]
  79. Brandt, B.B.; Endrodi, G. QCD phase diagram with isospin chemical potential. arXiv 2016, arXiv:1611.06758. [Google Scholar] [CrossRef]
  80. Brandt, B.B.; Endrodi, G. Reliability of Taylor expansions in QCD. Phys. Rev. D 2019, 99, 014518. [Google Scholar] [CrossRef] [Green Version]
  81. Brandt, B.B.; Endrodi, G.; Schmalzbauer, S. QCD at finite isospin chemical potential. EPJ Web Conf. 2018, 175, 07020. [Google Scholar] [CrossRef] [Green Version]
  82. Brandt, B.B.; Endrodi, G.; Schmalzbauer, S. QCD at nonzero isospin asymmetry. arXiv 2018, arXiv:1811.06004. [Google Scholar]
  83. Shapiro, S.L.; Teukolsky, S.A. Black Holes, White Dwarfs, and Neutron Stars: The Physics of Compact Objects; WILEY-VCH Verlag GmbH & Co. KgaA: Weinheim, Germany, 1983. [Google Scholar]
  84. Glendenning, N.K. Compact Stars: Nuclear Physics, Particle Physics, and General Relativity; Springer: New York, NY, USA, 1997. [Google Scholar]
  85. Migdal, A.B.; Saperstein, E.; Troitsky, M.; Voskresensky, D. Pion degrees of freedom in nuclear matter. Phys. Rept. 1990, 192, 179–437. [Google Scholar] [CrossRef]
  86. Cameron, A.G. Neutron Star Models. Astrophys. J. 1959, 130, 884. [Google Scholar] [CrossRef]
  87. Ambartsumyan, V.A.; Saakyan, G.S. The Degenerate Superdense Gas of Elementary Particles. Sov. Astron. 1960, 4, 187. [Google Scholar]
  88. Salpeter, E.E. Matter at high densities. Ann. Phys. 1960, 11, 393–413. [Google Scholar] [CrossRef]
  89. Ambartsumyan, V.A.; Saakyan, G.S. Internal Structure of Hyperon Configurations of Stellar Masses. Sov. Astron. 1962, 5, 779. [Google Scholar]
  90. Bahcall, J.N.; Wolf, R.A. Neutron Stars. 1. Properties at Absolute Zero Temperature. Phys. Rev. 1965, 140, B1445–B1451. [Google Scholar] [CrossRef]
  91. Giorgini, S.; Pitaevskii, L.P.; Stringari, S. Theory of ultracold atomic Fermi gases. Rev. Mod. Phys. 2008, 80, 1215–1274. [Google Scholar] [CrossRef] [Green Version]
  92. Alford, M.G.; Harris, S.P. Beta equilibrium in neutron star mergers. Phys. Rev. C 2018, 98, 065806. [Google Scholar] [CrossRef]
  93. Sadzikowski, M. Coexistence of pion condensation and color superconductivity in two flavor quark matter. Phys. Lett. B 2003, 553, 45–50. [Google Scholar] [CrossRef] [Green Version]
  94. Buballa, M.; Carignano, S. Inhomogeneous chiral condensates. Prog. Part. Nucl. Phys. 2015, 81, 39–96. [Google Scholar] [CrossRef] [Green Version]
  95. Carignano, S.; Mannarelli, M.; Anzuini, F.; Benhar, O. Crystalline phases by an improved gradient expansion technique. Phys. Rev. D 2018, 97, 036009. [Google Scholar] [CrossRef] [Green Version]
  96. Weinberg, S. Phenomenological Lagrangians. Physica 1979, 96, 327. [Google Scholar] [CrossRef]
  97. Gasser, J.; Leutwyler, H. Chiral Perturbation Theory to One Loop. Ann. Phys. 1984, 158, 142. [Google Scholar] [CrossRef]
  98. Georgi, H. Weak Interactions and Modern Particle Theory; Dover Publications: New York, NY, USA, 1984. [Google Scholar]
  99. Leutwyler, H. On the foundations of chiral perturbation theory. Ann. Phys. 1994, 235, 165–203. [Google Scholar] [CrossRef]
  100. Ecker, G. Chiral perturbation theory. Prog. Part. Nucl. Phys. 1995, 35, 1–80. [Google Scholar] [CrossRef] [Green Version]
  101. Leutwyler, H. Phonons as goldstone bosons. Helv. Phys. Acta 1997, 70, 275–286. [Google Scholar]
  102. Pich, A. Effective Field Theory: Course. In Probing the Standard Model of Particle Interactions, Proceedings of the Summer School in Theoretical Physics, NATO Advanced Study Institute, 68th session, Les Houches, France, 28 July–5 September 1997; North Holland: Amsterdam, The Netherlands, 1998; pp. 949–1049. [Google Scholar]
  103. Scherer, S. Introduction to chiral perturbation theory. Adv. Nucl. Phys. 2003, 27, 277. [Google Scholar]
  104. Scherer, S.; Schindler, M.R. A Chiral perturbation theory primer. arXiv 2005, arXiv:hep-ph/0505265. [Google Scholar]
  105. Kogut, J.B.; Stephanov, M.A.; Toublan, D. On two color QCD with baryon chemical potential. Phys. Lett. B 1999, 464, 183–191. [Google Scholar] [CrossRef]
  106. Kogut, J.B.; Stephanov, M.A.; Toublan, D.; Verbaarschot, J.J.M.; Zhitnitsky, A. QCD—Like theories at finite baryon density. Nucl. Phys. B 2000, 582, 477–513. [Google Scholar] [CrossRef]
  107. Hands, S.; Montvay, I.; Morrison, S.; Oevers, M.; Scorzato, L.; Skullerud, J. Numerical study of dense adjoint matter in two color QCD. Eur. Phys. J. C 2000, 17, 285–302. [Google Scholar] [CrossRef] [Green Version]
  108. Kogut, J.B.; Sinclair, D.K.; Hands, S.J.; Morrison, S.E. Two color QCD at nonzero quark number density. Phys. Rev. D 2001, 64, 094505. [Google Scholar] [CrossRef]
  109. Brauner, T. On the chiral perturbation theory for two-flavor two-color QCD at finite chemical potential. Mod. Phys. Lett. A 2006, 21, 559–570. [Google Scholar] [CrossRef]
  110. Braguta, V.V.; Ilgenfritz, E.M.; Kotov, A.Yu.; Molochkov, A.V.; Nikolaev, A.A. Study of the phase diagram of dense two-color QCD within lattice simulation. Phys. Rev. D 2016, 94, 114510. [Google Scholar] [CrossRef] [Green Version]
  111. Adhikari, P.; Beleznay, S.B.; Mannarelli, M. Finite Density Two Color Chiral Perturbation Theory Revisited. Eur. Phys. J. C 2018, 78, 441. [Google Scholar] [CrossRef]
  112. Le Bellac, M. Thermal Field Theory. In Cambridge Monographs on Mathematical Physics; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  113. Nambu, Y. Axial vector current conservation in weak interactions. Phys. Rev. Lett. 1960, 4, 380–382. [Google Scholar] [CrossRef]
  114. Nambu, Y.; Jona-Lasinio, G. Dynamical Model of Elementary Particles Based on an Analogy with Superconductivity. 1. Phys. Rev. 1961, 122, 345–358. [Google Scholar] [CrossRef]
  115. Nambu, Y.; Jona-Lasinio, G. Dynamical Model of Elementary Particles Based on An Analogy with Superconductivity. II. Phys. Rev. 1961, 124, 246–254. [Google Scholar] [CrossRef]
  116. Ebert, D.; Volkov, M.K. Composite Meson Model with Vector Dominance Based on U(2) Invariant Four Quark Interactions. Z. Phys. C 1983, 16, 205. [Google Scholar] [CrossRef]
  117. Ebert, D.; Reinhardt, H. Effective Chiral Hadron Lagrangian with Anomalies and Skyrme Terms from Quark Flavor Dynamics. Nucl. Phys. B 1986, 271, 188–226. [Google Scholar] [CrossRef]
  118. Buballa, M. NJL model analysis of quark matter at large density. Phys. Rept. 2005, 407, 205–376. [Google Scholar] [CrossRef]
  119. Alford, M.G.; Rajagopal, K.; Wilczek, F. Color flavor locking and chiral symmetry breaking in high density QCD. Nucl. Phys. B 1999, 537, 443–458. [Google Scholar] [CrossRef]
  120. Asakawa, M.; Yazaki, K. Chiral Restoration at Finite Density and Temperature. Nucl. Phys. A 1989, 504, 668–684. [Google Scholar] [CrossRef]
  121. Bernard, V.; Jaffe, R.L.; Meissner, U.G. Strangeness Mixing and Quenching in the Nambu-Jona-Lasinio Model. Nucl. Phys. B 1988, 308, 753–790. [Google Scholar] [CrossRef]
  122. Smit, J. Introduction to quantum fields on a lattice: A robust mate. Camb. Lect. Notes Phys. 2002, 15, 1–271. [Google Scholar]
  123. Gattringer, C.; Lang, C.B. Quantum chromodynamics on the lattice. Lect. Notes Phys. 2010, 788, 1–343. [Google Scholar] [CrossRef]
  124. Muroya, S.; Nakamura, A.; Nonaka, C.; Takaishi, T. Lattice QCD at finite density: An Introductory review. Prog. Theor. Phys. 2003, 110, 615–668. [Google Scholar] [CrossRef]
  125. Schmidt, C. Lattice QCD at finite density. PoS 2006, LAT2006, 021. [Google Scholar] [CrossRef]
  126. De Forcrand, P. Simulating QCD at finite density. PoS 2009, LAT2009, 010. [Google Scholar] [CrossRef]
  127. Philipsen, O. The QCD equation of state from the lattice. Prog. Part. Nucl. Phys. 2013, 70, 55–107. [Google Scholar] [CrossRef] [Green Version]
  128. Aarts, G. Introductory lectures on lattice QCD at nonzero baryon number. J. Phys. Conf. Ser. 2016, 706, 022004. [Google Scholar] [CrossRef]
  129. Lombardo, M.P. Finite density (might well be easier) at finite temperature. Nucl. Phys. Proc. Suppl. 2000, 83, 375–377. [Google Scholar] [CrossRef]
  130. Cea, P.; Cosmai, L.; D’Elia, M.; Papa, A.; Sanfilippo, F. The critical line of two-flavor QCD at finite isospin or baryon densities from imaginary chemical potentials. Phys. Rev. D 2012, 85, 094512. [Google Scholar] [CrossRef]
  131. Nishida, Y. Phase structures of strong coupling lattice QCD with finite baryon and isospin density. Phys. Rev. D 2004, 69, 094501. [Google Scholar] [CrossRef] [Green Version]
  132. Detmold, W. Nuclear Physics from Lattice QCD. Lect. Notes Phys. 2015, 889, 153–194. [Google Scholar] [CrossRef]
  133. Brandt, B.B.; Endrodi, G.; Schmalzbauer, S. QCD phase diagram for nonzero isospin-asymmetry. Phys. Rev. D 2018, 97, 054514. [Google Scholar] [CrossRef] [Green Version]
  134. Cheng, T.; Li, L. Gauge Theory of Elementary Particle Physics; Oxford science publications, Clarendon Press: Oxford, UK, 1984. [Google Scholar]
  135. Brandt, B.B.; Endrődi, G.; Fraga, E.S.; Hippert, M.; Schaffner-Bielich, J.; Schmalzbauer, S. New class of compact stars: Pion stars. Phys. Rev. D 2018, 98, 094510. [Google Scholar] [CrossRef] [Green Version]
  136. Andersen, J.O.; Kneschke, P. Bose-Einstein condensation and pion stars. arXiv 2018, arXiv:hep-ph/1807.08951. [Google Scholar]
  137. Barducci, A.; Pettini, G.; Ravagli, L.; Casalbuoni, R. Ladder QCD at finite isospin chemical potential. Phys. Lett. B 2003, 564, 217–224. [Google Scholar] [CrossRef]
  138. Toublan, D.; Kogut, J.B. The QCD phase diagram at nonzero baryon, isospin and strangeness chemical potentials: Results from a hadron resonance gas model. Phys. Lett. B 2005, 605, 129–136. [Google Scholar] [CrossRef]
  139. Schafer, T.; Wilczek, F. Continuity of quark and hadron matter. Phys. Rev. Lett. 1999, 82, 3956–3959. [Google Scholar] [CrossRef]
Figure 1. Cartoon of the grand-canonical phase diagram of hadronic matter as a function of temperature and of the isospin and baryonic chemical potentials. The shaded green regions are explored by heavy ion collisions (HIC) or are possibly realized in compact stars. For each phase we have indicated the relevant quark condensates, see Equations (7)–(9) below.
Figure 1. Cartoon of the grand-canonical phase diagram of hadronic matter as a function of temperature and of the isospin and baryonic chemical potentials. The shaded green regions are explored by heavy ion collisions (HIC) or are possibly realized in compact stars. For each phase we have indicated the relevant quark condensates, see Equations (7)–(9) below.
Particles 02 00025 g001
Figure 2. Sketch of the Pauli blocking mechanism for the charged pion decay (solid red line) induced by an increasing electron chemical potential. With increasing electron density, the electron chemical potential grows. The equilibrium configuration is reached for sufficiently high electron density, last figure on the right, with the pion decay process equilibrated by the electron decay process (blue dashed line), corresponding to μ e = m π .
Figure 2. Sketch of the Pauli blocking mechanism for the charged pion decay (solid red line) induced by an increasing electron chemical potential. With increasing electron density, the electron chemical potential grows. The equilibrium configuration is reached for sufficiently high electron density, last figure on the right, with the pion decay process equilibrated by the electron decay process (blue dashed line), corresponding to μ e = m π .
Particles 02 00025 g002
Figure 3. Symmetry breaking path of three flavor quark matter; the arrows indicate the various symmetry breakings. On the right we report the emergent low-energy degrees of freedom. The normal phase is defined by the corresponding symmetry group, see Equation (35). See the text below for a through description of the symmetry breakings.
Figure 3. Symmetry breaking path of three flavor quark matter; the arrows indicate the various symmetry breakings. On the right we report the emergent low-energy degrees of freedom. The normal phase is defined by the corresponding symmetry group, see Equation (35). See the text below for a through description of the symmetry breakings.
Particles 02 00025 g003
Figure 4. Sketch of the energy levels of pions and kaons. On the left we assume massless quarks and vanishing chemical potentials. The effect of the explicit symmetry breaking induced by the quark masses is shown in the central region. On the right we show the mass splitting within isomultiplets, obtained for μ I = 80 MeV and μ S = 90 MeV. The isospin and the strangeness chemical potentials completely remove the level degeneracy.
Figure 4. Sketch of the energy levels of pions and kaons. On the left we assume massless quarks and vanishing chemical potentials. The effect of the explicit symmetry breaking induced by the quark masses is shown in the central region. On the right we show the mass splitting within isomultiplets, obtained for μ I = 80 MeV and μ S = 90 MeV. The isospin and the strangeness chemical potentials completely remove the level degeneracy.
Particles 02 00025 g004
Figure 5. Phase diagram of the meson condensed phase, see [11]. The irregular hexagon (solid red line), obtained from the meson mass spectrum, corresponds to the second order phase transitions between the normal phase and the meson condensed phases. The dashed lines have been drawn by hand and represent the expected first order phase transitions between two different meson condensed phases.
Figure 5. Phase diagram of the meson condensed phase, see [11]. The irregular hexagon (solid red line), obtained from the meson mass spectrum, corresponds to the second order phase transitions between the normal phase and the meson condensed phases. The dashed lines have been drawn by hand and represent the expected first order phase transitions between two different meson condensed phases.
Particles 02 00025 g005
Figure 6. Weight diagram of the mesonic octet. We have enclosed with the same symbol the states that can mix in the meson condensed phase. They have the same T-spin, U-spin and V-spin quantum numbers. The π 0 and the η are not simultaneous T-spin, U-spin and V-spin eigenstates; their mixing depends on the spontaneously induced charge of the vacuum, see [26] for more details.
Figure 6. Weight diagram of the mesonic octet. We have enclosed with the same symbol the states that can mix in the meson condensed phase. They have the same T-spin, U-spin and V-spin quantum numbers. The π 0 and the η are not simultaneous T-spin, U-spin and V-spin eigenstates; their mixing depends on the spontaneously induced charge of the vacuum, see [26] for more details.
Particles 02 00025 g006
Figure 7. Comparison of the chiral and pion condensates obtained by χ PT (dashed black line), by the two-flavor NJL model (solid blue line) and by LQCD simulations on a 6 × 24 3 lattice by [133] (orange squares). The condensates have been normalized at the value of the scalar condensate in the normal phase (note that the LQCD data of [133] have been obtained at T = 113 MeV), while the chemical potential is in units of the pion mass in the normal phase. Quite remarkably, the condensates obtained with the three methods overlap for μ I 2 m π . For larger values of the isospin chemical potentials the evaluation of the pion condensate becomes problematic: the NJL model feels the hard cutoff Λ ; the χ PT approaches the critical scale Λ χ ; the LQCD simulations start to feel the lattice saturation effects [70].
Figure 7. Comparison of the chiral and pion condensates obtained by χ PT (dashed black line), by the two-flavor NJL model (solid blue line) and by LQCD simulations on a 6 × 24 3 lattice by [133] (orange squares). The condensates have been normalized at the value of the scalar condensate in the normal phase (note that the LQCD data of [133] have been obtained at T = 113 MeV), while the chemical potential is in units of the pion mass in the normal phase. Quite remarkably, the condensates obtained with the three methods overlap for μ I 2 m π . For larger values of the isospin chemical potentials the evaluation of the pion condensate becomes problematic: the NJL model feels the hard cutoff Λ ; the χ PT approaches the critical scale Λ χ ; the LQCD simulations start to feel the lattice saturation effects [70].
Particles 02 00025 g007
Figure 8. Comparison of the energy density over the Stefan-Boltzmann limit obtained by different methods. We report the lattice data points of the canonical simulations of [76] obtained at T = 20 MeV with three different lattice volumes. The pQCD results of [66] (orange dotted line) indicate a constant asymptotic energy density. The χ PT leading order results of [27] (dashed black line) and the mean field NJL results (solid blue line), see [38,49], perfectly reproduce the peak structure at μ I 1.27 m π .
Figure 8. Comparison of the energy density over the Stefan-Boltzmann limit obtained by different methods. We report the lattice data points of the canonical simulations of [76] obtained at T = 20 MeV with three different lattice volumes. The pQCD results of [66] (orange dotted line) indicate a constant asymptotic energy density. The χ PT leading order results of [27] (dashed black line) and the mean field NJL results (solid blue line), see [38,49], perfectly reproduce the peak structure at μ I 1.27 m π .
Particles 02 00025 g008
Figure 9. Phase diagram of hadronic matter at μ B = 0 . The green shaded area marked with LQCD corresponds to the second order phase transition separating the pion condensed phase from the normal phase, at low temperature, and the chirally restored quark gluon plasma, at high temperature. These three phases meet at a (pseudo) tricritical point (orange dot with error bars), at T tricritical 151 MeV and μ I , tricritical 140 MeV. The shaded blue area marks the chiral crossover line which in part overlaps with a (probable) deconfinement phase transition (shaded red area), see the text for more details. These results have been obtained in [133] by LQCD simulations. The NJL second order transition line (solid blue), see [38], overlaps with the LQCD data for T 100 MeV. Then it shows a stronger temperature dependence. The χ PT second order transition line (dashed black), corresponding to Equation (108), shows a temperature dependence that disagrees with both the LQCD and NJL results.
Figure 9. Phase diagram of hadronic matter at μ B = 0 . The green shaded area marked with LQCD corresponds to the second order phase transition separating the pion condensed phase from the normal phase, at low temperature, and the chirally restored quark gluon plasma, at high temperature. These three phases meet at a (pseudo) tricritical point (orange dot with error bars), at T tricritical 151 MeV and μ I , tricritical 140 MeV. The shaded blue area marks the chiral crossover line which in part overlaps with a (probable) deconfinement phase transition (shaded red area), see the text for more details. These results have been obtained in [133] by LQCD simulations. The NJL second order transition line (solid blue), see [38], overlaps with the LQCD data for T 100 MeV. Then it shows a stronger temperature dependence. The χ PT second order transition line (dashed black), corresponding to Equation (108), shows a temperature dependence that disagrees with both the LQCD and NJL results.
Particles 02 00025 g009
Figure 10. Sketch of the phase diagram of hadronic matter. The solid thick lines have been determined in the LQCD simulations of [133]. For μ B = 0 and T 150 MeV, the transition to the pion condensed phase happens at μ I = m π . The shaded area on the top corresponds to a chiral and deconfinement crossover. We have indicated with thin blue lines the possible transitions from the hadron gas phase to the color superconducting and to the pion condensed phases, although there are neither experimental data nor LQCD simulations that support their existence. How the color superconducting phase turns in a quark-gluon plasma or in a pion condensed phase is not known.
Figure 10. Sketch of the phase diagram of hadronic matter. The solid thick lines have been determined in the LQCD simulations of [133]. For μ B = 0 and T 150 MeV, the transition to the pion condensed phase happens at μ I = m π . The shaded area on the top corresponds to a chiral and deconfinement crossover. We have indicated with thin blue lines the possible transitions from the hadron gas phase to the color superconducting and to the pion condensed phases, although there are neither experimental data nor LQCD simulations that support their existence. How the color superconducting phase turns in a quark-gluon plasma or in a pion condensed phase is not known.
Particles 02 00025 g010

Share and Cite

MDPI and ACS Style

Mannarelli, M. Meson Condensation. Particles 2019, 2, 411-443. https://doi.org/10.3390/particles2030025

AMA Style

Mannarelli M. Meson Condensation. Particles. 2019; 2(3):411-443. https://doi.org/10.3390/particles2030025

Chicago/Turabian Style

Mannarelli, Massimo. 2019. "Meson Condensation" Particles 2, no. 3: 411-443. https://doi.org/10.3390/particles2030025

APA Style

Mannarelli, M. (2019). Meson Condensation. Particles, 2(3), 411-443. https://doi.org/10.3390/particles2030025

Article Metrics

Back to TopTop