Charged pion condensation in dense quark matter: Nambu--Jona-Lasinio model study

In this short review we tried to give an outline of investigations of charged pion condensation (PC) in dense baryonic (quark) matter in the framework of effective Nambu--Jona-Lasinio (NJL) type models. The possibility of charged PC phase in dense quark matter with isospin asymmetry is investigated. First, it is demonstrated that this phase can be realized in the framework of massless NJL model. But the existence of this phase is enormously fragile to the values of current quark mass %.Then, and we show that charged PC phase is forbidden in electrically neutral dense quark matter with $\beta$-equilibrium when current quark masses are close to their physical value of 5.5 MeV. Nevertheless, then it is shown that in real physical systems there could be conditions promoting the appearance of charged PC phenomenon in dense quark matter, namely, it was shown that if one includes into consideration the fact that system can have finite size, then a dense charged PC phase can be realized there. It was also demonstrated that the possibility of inhomogeneous pion condensate might allow this phase to appear. And more recently it was revealed that there is another interesting factor that can induce a charged PC phase in dense quark matter even without isospin imbalance. It is a chiral imbalance of the system (non-zero difference between densities of left- and right-handed quarks). This results can be interesting in heavy ion collision experiments, where it is expected to get high baryon densities. It is of interest also in the context of the neutron stars, where quark matter might be realized in the core and very high baryon and isospin densities are attained.


I. INTRODUCTION
Recently, great efforts have been made in trying to understand the properties and the phase diagram of strongly interacting (quark) matter under high temperatures and/or baryon densities (or high values of the baryon chemical potential µ B ).It is believed that quark matter under such extreme conditions can exist in cores of neutron stars, where baryon density is much higher than the density ρ 0 of ordinary nuclear matter, where ρ 0 is 0.15 baryon per fm −3 , or can be formed as a result of a collision of heavy ions [1].(In this case there are no reasons to speak of protons and neutrons as particles that make up baryonic matter, and it would be more correct to say that we are dealing with dense quark matter.)However, in reality, within the framework of quantum chromodynamics (QCD), the fundamental theory designed to describe strongly interacting systems, it is rather problematic to get at least some information about the phase structure of dense quark matter.The fact is that despite asymptotic freedom and feasibility of perturbative calculations at extremely high temperature and the baryon density in real physical conditions (in compact stars or heavy ion collision experiments) the temperature and the baryon density of quark matter are usually not that great, which leads to the fact that the effective interaction constant of the particles is quite large.Therefore, one of the main research methods in the framework of QCD, perturbative analysis over a small coupling constant, is not applicable.Another well-known first-principle calculation approach to QCD is the method of lattice QCD simulations.But, unfortunately, due to the sign problem (complex fermion determinant) it encounters up to now insurmountable difficulties at µ B = 0.
Due to these reasons, in order to study the QCD phase diagram at nonzero baryon density, effective field theories are employed and the most frequently used one is the NambuJona-Lasinio (NJL) model [2] (for a review one can see Refs.[3,4]).This model uses quarks as the degrees of freedom instead of hadrons, and it can be used as an effective field theory for QCD in the low energy regime in order to describe dense quark matter.The most appealing feature of NJL model is the description of the phenomenon of the dynamical chiral symmetry breaking when quarks get a rather large mass compared to the value of the current quark mass.The NJL model can be used to get insights not only in meson physics [5], but also in color superconductivity phenomenon [6,7], etc.
As it was mentioned above, dense quark matter can be formed inside neutron stars, i.e. it consists of different numbers of u and d quarks (it is the so-called isotopic asymmetry of quark matter).Hence, the ground state of such matter is characterized by baryon and isotopic (or isospin) densities and can be investigated in the framework of QCD or other alternative effective approaches using quark degrees of freedom and extended by baryon µ B and isotopic (or isospin) µ I chemical potentials.Then, as it was shown in Refs.[8][9][10][11][12][13][14][15][16][17][18][19], in the corresponding systems there might appear a Bose-Einstein condensate (even at µ B = 0) of quark-antiquark pairs of the form ūγ 5 d, which has quantum numbers of charged pions.In this case the charged pion condensation (PC) phase is realized in matter.However, the presence of this phase of dense quark matter is not a well-established fact up to now.Indeed, it was shown in Refs.[8,9] that charged PC phase with nonzero quark density occupies a very narrow part of the phase diagram of the massless NJL model when µ B = 0 and µ I = 0.But this phase is extremely fragile to the effect of explicit chiral symmetry breaking via a non-zero current quark mass and is ruled out from the phase diagram of electrically neutral dense quark matter with β-equilibrium, when quark masses are larger than or of the order of 10 keV (see, e.g., in Refs.[17,[20][21][22][23]).The same results were found in the framework of the Polyakov NJL model [22], where it was shown that the value of isospin chemical potential µ I turns out to be not large enough for charged PC phase to takes place in electrically neutral dense quark matter.Similar qualitative picture of the charged PC has been also found in an analysis in the framework of the toy (1+1)-dimensional massive Gross-Neveu model with isospin and baryon chemical potentials [24].In particular, for example, it was found that baryon density is equal to zero in the charged PC phase of this model, i.e. this model does not predict dense quark matter with charged PC phenomenon.
Before continuing our discussion on the fate of the charged PC phase in dense quark matter, let us say a few words about the PC phenomenon in dense nuclear matter.It is clear that in this case the baryon density is around or greater than the normal baryon density ρ 0 , i.e. much less than in dense quark matter.It is clear that the main elementary objects forming nuclear matter are nucleons and pions.Therefore, a long time ago, in the early seventies, it was understood that pions and pion condensates can play a significant role in the physical processes of dense nuclear matter [25] (see also, for example, the reviews [26,27]).As for physics of neutron stars, where in nuclear matter there is a large isotopic asymmetry (the predominance of neutrons over protons), here, of course, processes with the formation of charged pion condensate can play a prominent role.Indeed, in a series of papers [28][29][30] it was shown that charged PC phase has a significant effect on the cooling scenario of a neutron star, on its equation of state, on the phenomenon of superfluidity of nucleons in neutron stars, etc.The main thing for us is that different effective approaches (such as the σ-model, etc.) used in these studies to describe neutron matter, in principle, do not reject the possibility of the existence of a phase with charged PC. (Though let us note that the s-wave charged PC is considered highly unlikely to be realized in dense matter [31] but it was argued that p-wave charged PC is possible [32,33].)In contrast, as it follows from the above discussion, the possibility of the charged PC phenomenon in dense quark matter, in which baryon density is much higher than ρ 0 , described in the framework of the effective NJL model remains in question.
In this mini-review we discuss the present status of the charged PC phase in dense quark matter and, especially, pay attention to the factors promoting the emergence of charged PC phenomenon in dense quark matter.Namely, in the paper [34] it was found that this phase might be realized in dense baryonic system that has finite size.It was also demonstrated in [37] that in dense baryonic matter spatially inhomogeneous charged PC might appear, and more recently it was found in [38][39][40][41][42] that chiral imbalance is another interesting factor that can induce a charged PC phase.Let us recall that chiral imbalance is a nonzero difference between densities of left-handed and right-handed fermions and it may appear in the fireball after heavy ion collisions and possibly lead to the so-called chiral magnetic effect [43,44].It is a remarkable phenomenon that stems from highly nontrivial interplay of chiral symmetry, chiral (axial) anomaly of QCD and the topological structure of gluon configurations.It might be also realized in compact stars [42] or even in condensed matter systems [45,46].
The promotion of the charged PC phase in dense quark matter by finite size and spatial inhomogeneity of the charged pion condensate was demonstrated using a (1+1)-dimensional toy model with four-quark interactions, whereas the generation of the charged PC phase by chiral imbalance was shown as in (1+1)-dimensional toy model [38,39] as well as in more realistic effective QCD motivated NJL model [40,41].
The paper is organized as follows.In Sec.II NJL model with two quark flavors (u and d quarks) and with two kinds of chemical potentials, µ B and µ I , is presented and the charged PC phenomenon in dense isospin asymmetric quark matter is discussed.Here we show that it is unlikely to appear in the simplest case.Then in Sec.III it is revealed that there exist the conditions and factors promoting the appearance of charged PC in dense quark matter and different possibilities are discussed.In Sec.III A the influence of finite size effect and nontrivial space topology is considered, then in Sec.III B the possibility of inhomogeneous charged PC is discussed, finally in Sec.III C it is shown that charged PC can appear in dense chiral imbalanced quark matter.Sec.IV contains the summary and conclusions.Throughout the paper, for simplicity, all the considerations are made in the case of zero temperature T = 0.

A. The model and its thermodynamic potential
Since the considered studies are mostly performed in the framework of the (3+1)-dimensional NJL model, let us recall the two flavored NJL model Lagrangian (there are no gluons and it is symmetric under global color SU (N c ) group).It has the following form where q = (q u , q d ) T is a flavor doublet and q u and q d (u and d quarks) are four-component Dirac spinors as well as color N c -plets.The summation over flavor, color, and spinor indices is implied; And m is current quark mass.Moreover, since in (3+1)-spacetime dimensions the four-fermion interaction operator has a mass dimension 6, the NJL model is not renormalizable and may be considered only as an effective QCD-like model.Hence, some regularization is needed in order to avoid usual ultraviolet divergences and throughout the paper we will use the simplest momentum cutoff regularization, i.e. the integration region in all momentum integrals is restricted by a cutoff Λ, | p| < Λ.And the model parameters are fixed as follows: G = 15.03GeV −2 , Λ = 0.65 GeV.
Since our paper is devoted to the dense baryonic (quark) matter at various conditions, now, as a starting point, let us first consider the dense quark matter with isospin asymmetry.The Lagrangian of two-flavored NJL model with baryon and isopin chemical potentials has the following form Since the Lagrangian (2) contains baryon µ B and isospin µ I chemical potentials, it can describe quark matter with nonzero baryon n B and isospin n I densities.Baryon n B and isospin n I densities are the quantities, thermodynamically conjugated to chemical potentials µ B and µ I , respectively, and can be obtained if one takes a derivative of the thermodynamic potential (TDP) of the system (2) with respect to µ B and µ I chemical potentials.The quantities n B and n I are densities of conserved charges corresponding to the invariance of Lagrangian (2) with respect to the abelian groups U B (1), U I3 (1), where1 U B (1) : q → exp(iα/3)q; U I3 (1) : q → exp(iατ 3 /2)q. (3) So one can show that n B = qγ 0 q/3 and n I = qγ 0 τ 3 q/2.One can also see that, in addition to (3), Lagrangian (2) is invariant under the transformation of the electromagnetic U Q (1) group, where In order to find the thermodynamic potential of the system it is convenient to use a semibosonized version of the Lagrangian (2), with auxiliary composite bosonic fields σ(x) and π a (x) (a = 1, 2, 3), Here the summation over repeated indices is implied.Hereafter, we will use the notations µ ≡ µ B /3 for quark chemical potential and ν ≡ µ I /2.One can get the equations of motion for the bosonic fields from the auxiliary Lagrangian (5), Utilizing the equations (6), it is possible to show that the semibosonized Lagrangian L (5) is equivalent to the initial Lagrangian (2).Note that the auxiliary composite bosonic field π 3 (x) can be identified with the physical π 0 (x)meson field, whereas the physical π ± (x)-meson fields are the following combinations of the composite bosonic fields, In addition, let us note that the composite bosonic fields (6) are transformed under the isospin U I3 (1) group in the following manner (it can be shown from (6) and footnote 1): From the auxiliary Lagrangian (5), it is possible to obtain in the leading order of the large-N c expansion (i.e. in the one loop approximation) the effective action S eff (σ, π a ) of the model.This quantity is a functional of the auxiliary bosonic σ(x) and π a (x) fields and has the following form: It is easy to obtain for S eff (σ, π a ) the following expression in which the last term (i.e. the term Seff in (8)) is a quark contribution to the effective action, where N ′ is a normalization constant and D ≡ γ ν i∂ ν − m + µγ 0 + ντ 3 γ 0 − σ(x) − iγ 5 π a (x)τ a is the Dirac operator acting in the flavor, spinor and coordinate spaces.Now let us use the general formula Det D = exp Tr ln D and obtain the following expression for the effective action (8) where the Tr operation is performed in spinor (s), flavor (f ) and four-dimensional space-time coordinate (x) spaces, respectively.Now the ground state expectation values σ(x) and π a (x) of composite bosonic fields can be found from the saddle point equations, where a = 1, 2, 3.The knowledge of values σ(x) and π a (x) , namely their behaviour with respect to chemical potentials, provides us with phase structure of the model.It is clear that if there is a non-zero ground state expectation values of composite bosonic σ(x) and/or π 3 (x) fields, i. e. σ(x) = 0 and/or π 3 (x) = 0, then in the chiral limit we have a spontaneous breaking of the axial isospin U AI3 (1) symmetry in the model (if we consider the chiral limit m = 0, when the Lagrangian (2) is symmetric in addition with respect to U AI3 (1) group, which is defined as U AI3 (1) : q → exp(iαγ 5 τ 3 /2)q), whereas if other component of pion field acquire non-zero ground state expectation value, i. e. π 1 (x) = 0 and/or π 2 (x) = 0, then, as it is obvious from (7), the isospin symmetry (U I3 (1) group) is spontaneously broken down.Since in this latter case the condensates, or the ground state expectation values, of the charged pion fields π + (x) and π − (x) are non-zero, this phase is usually called the charged PC phase.Furthermore, one can see from (6) that the nonzero condensates π 1,2 (x) (or π ± (x) ) are not invariant with respect to the transformations of the electromagnetic group U Q (1) (4) and this leads to spontaneous breaking of electromagnetic U Q (1) symmetry.Hence one can say that the electromagnetic invariance in the charged PC phase is broken down and superconductivity is an unavoidable property of the charged PC phase.
As a rule one supposes that in the ground state of the system expectation values σ(x) and π a (x) do not depend on space-time coordinates x, where M and π a (a = 1, 2, 3) are already constant quantities and they can be found as coordinates of the global minimum point of the thermodynamic potential Ω(M, π a ).The TDP is defined by the following expression in the leading order of the large-N c expansion: To find the phase portrait of the model one needs to investigate the dependence on µ, ν, of the global minimum point (GMP) of the TDP Ω(σ, π a ) vs σ, π a .To this end let us simplify the task and note that due to a U I3 (1) invariance of the model, the TDP (13) depends really only on the three combinations of bosonic fields, σ, π 3 and π 2 1 + π 2 2 (this can be seen from ( 7)).So without loss of generality, it can be assumed that π 2 = 0.Moreover, it is also possible to prove that π 3 = 0 at m = 0 (see, e.g., in [47]).Hence, one can study the TDP as a function of only two variables, M ≡ σ + m and ∆ ≡ π 1 .It means that we use the following ansatz The expression for the TDP (13) of the dense (baryon density) and isospin asymmetric quark matter can be easily obtained in the mean field approximation and it has the form where there have been introduced the following notations , and Let us recall also that the integration in the momentum space is restricted by a cutoff | p| < Λ.From (15) it is possible to obtain the gap equations Note that the quark gap M is just the dynamical quark mass.First, let us discuss the phase diagram and the charged PC in dense quark matter with isospin asymmetry in the chiral limit.This case has been considered and the phase portrait of the model at both µ I = 0 and µ = 0 has been obtained in [8,40].Since, in general, it is a rather hard task to consider the phase diagram at the physical point, in this paper, for simplicity, the current quark mass m was assumed to be zero, i. e. the chiral limit was used.The corresponding (ν, µ)-phase portrait is depicted in Figure 1.Let us discuss its content.One can see there is a huge region of PC phase at this phase diagram.Note that it has been shown already in [10][11][12]50] that at m = 0 the charged PC phase appears at µ I > m π .But in the chiral limit (see in Figure 1) the same relation µ I > m π is also valid since m π = 0 in this case.Now let us discuss all possible phases in this case.The system of the gap equations ( 16) at m = 0 has three different types of solutions: (i) M = 0, ∆ = 0, (ii) M = 0, ∆ = 0, (iii) M = 0, ∆ = 0.The (i) solution corresponds to the symmetric phase, where all condensates are zero.If the GMP of the TDP is located at the point of the form (ii), then the charged pion condensate is equal to zero in the ground state and the quark condensate is nonzero.This phase is called a chiral symmetry breaking (CSB) phase.Finally, the (iii) solution corresponds to charged PC phase, where the quark condensate is equal to zero, but the charged pion one has a nonzero value.
If the GMP of the TDP is achieved on the solution of type (i), then one can see that the ground state of the system is invariant with respect to the initial symmetry group of the model, which at m = 0 (in the chiral limit) and (1).The (ii)-type solution (CSB phase) corresponds to the phase symmetric with respect to U B (1) × U I3 (1) group.In the charged PC phase the ground state is invariant under transformations from U B (1) × U AI3 (1) group.Moreover, it is clear that in the ground state of quark matter with non-zero charged pion condensate parity is also broken down.
Let us return to the discussion of the (ν, µ)-phase portrait of the model (see in Figure 1).One can see that apart from the PC phase there is a region of the PC d phase at the phase diagram of Figure 1.This phase denotes the ground state of quark matter, in which besides non-zero charged pion condensate the baryon density is also non-zero.This consideration in terms of NJL model shows that the charged PC can be realized in dense quark matter with isospin asymmetry.But let us recall once again that considerations have been performed in the chiral limit when pions are a genuine Goldstone bosons and have zero masses.
C. Account of electric neutrality and β-equilibrium conditions in the chiral limit After the prediction of the charged PC phenomenon in the dense and isospin asymmetric quark matter (that can be realized, for example, in relativistic heavy ion collisions or inside the cores of the compact stars) it was interesting to take into account the charge neutrality and β-equilibrium constraints that are important inside compact stars.The possibility of the charged PC phenomenon in cold and electrically neutral dense baryonic matter was investigated in β-equilibrium in [9], where the consideration was also performed in the chiral limit.
Let us first elaborate on how to implement all these constraints.As it has been said above, the electric charge is conserved in the model and it can be written as Q = diag(2/3, −1/3) = I 3 + B/2, where I 3 = τ 3 /2 is the generator of the third isospin component in the flavor space and B = diag(1/3, 1/3) is the baryon charge generator.Due to the β-equilibrium requirement, one should incorporate electrons in our consideration.They are necessary to neutralize the electric charge of quarks.Moreover, β-equilibrium of quark matter means that all β-processes that include quarks and leptons should go in both directions with equal rates.As it was argued above, full Lagrangian of the system should include the electron field and electron chemical potential, but since our system is in β-equilibrium, one can show that the electron chemical potential should be equal to charge chemical potential of quarks with opposite sign.
It was also supposed that one can neglect the mass of the electron and put it to zero for simplicity for their mass is quite small in the considered scale.Since in the system composed of quarks and electrons both baryon charge and electric charge are conserved quantities, its properties in equilibrium are described by the following Lagrangian where L N JL is the Lagrangian of NJL model ( 1), e is the electron spinor field; N B and N Q are baryon and electric charge density operators, respectively, and µ B , µ Q are corresponding chemical potentials.It is clear that Hence the Lagrangian ( 17) can be presented in the form In particular, it is clear from this relation that µ Q is just the isospin chemical potential µ I .The TDP in this case was obtained in [9] in the mean field approximation.In contrary to the previous sections where the large N c approximation was used, here (since in neutral quark matter the electrons are taken into account and the count of degrees of freedom matters) it is important that there are three colors of quarks (N c = 3) and one should work in the mean field approximation.One should be careful as in this paper slightly different definitions are used compared to [9].In the chiral limit the TDP has the following form (the definition of the TDP is the same as in (13)) where the first term in the right hand side of ( 20) is the TDP of free massless electrons.The summation in the third term of (20) runs over all quasiparticles (a = u, d, ū, d), where and As always, in order to avoid usual ultraviolet divergences, the integration region in (20) is restricted by a cutoff Λ, i.e. | p| < Λ.The gap equations that can be obtained from (20) In order to impose the neutrality constraint locally, we should look for the ground state of the system with the electric charge density n Q ≡ −N c ∂Ω/∂µ Q = −3∂Ω/∂µ Q equals identically to zero.One can rephrase it in other words as the study of the GMP of the TDP Ω(M, ∆) under the constraint n Q = 0 that has the following form Then, it was shown in [9] by numerical investigations of the gap and neutrality equations ( 22) and ( 23) that at sufficiently small quark chemical potential, namely, at µ < µ 1c ≈ 301 MeV, the ground state of the system corresponds to CSB phase.As it has been discussed above, in this case the GMP has the form (ii) and the dynamical quark mass M ≈ 301 MeV.There is no need to neutralize the quark electric charge because in this phase both the baryon and the isospin densities equals to zero.It means that µ I = µ Q = 0 in this phase.At µ 1c < µ < µ 2c ≈ 323. 5 MeV the GMP has already the form (iii), meaning that in this case the charged PC phase is realized.In this phase the baryon density is non-zero so it is PC d phase.At Figure 2 the pion condensate ∆ 0 is depicted as a function of µ.Then, at larger values of the chemical potential µ (i.e. at µ 2c < µ at Figure 2) symmetric phase is realized, where both condensates are zero.So one can see that in electrically neutral dense baryonic matter in β-equilibrium there might appear charged PC phenomenon.Let us stress once more that in [9] this result was obtained in the chiral limit.Though current quark mass is small compared to the scales studied in the paper, it was shown in [10][11][12] that charged PC phenomenon at not so large values of isospin imbalance is strongly influenced by current quark mass values.As we have already said above, charged PC phenomenon takes place only at values of µ I > m π .And at physical point m π = 0 it happens only at values of isospin chemical potential larger than approximately 139 MeV.So non-zero current quark mass can be of importance for the investigations of the phase structure.
After the prediction of the charged PC phase of dense matter, it was necessary to consider this phenomenon at the physical point.It means that one needs to use the Lagrangian (1) and ( 2) with non-zero physical value of current (bare) quark mass m ≈ 5.5 MeV.The consideration at the physical point but without electric neutrality condition showed that there remains marginal region of PC d phase at the phase portrait.The corresponding (ν, µ) phase diagram is shown at Figure 3 and one can see that in this case the region of PC d phase is much smaller than at Figure 1 corresponding to the chiral limit.Although, still charged PC phenomenon can take place in this case even at the physical point (even though one should bear in mind that very small size of this phase and the fact that it is obtained in the framework of the effective model makes it not a very robust prediction).Now let us turn to the discussion on what happens if one includes into play the electric neutrality and betaequilibrium conditions at the physical point.In [17,20,21] the consideration of the electrically neutral dense baryonic matter in β-equilibrium has been performed in the framework of NJL model at the physical point.
In this case the expression for the TDP (20) should be slightly changed, namely, the non-zero current quark mass should be added and one obtains Numerical studies of this case have been performed in [17,20,21] and the phase structure has been obtained.In these studies the current quark mass m has been treated as a free parameter, therefore, the corresponding pion mass m π at µ = T = 0 was a free parameter as well.
Let us say a couple of words on the results obtained.At Figure 4 the (µ, m π ) phase diagram is presented.The solid line denotes the border between the two phases, CSB phase (which is arranged at rather small values of µ) and ApprSYM one (this phase lies in the region of large µ).ApprSYM means an approximate symmetrical phase and corresponds to a GMP of the TDP with M ≈ m = 0 and ∆ = 0, but in contrast to the CSB phase, dynamical quark mass M in ApprSYM drops rapidly and incessantly to the current quark mass value m at increasing temperature or baryon chemical potential.The bold dot at the end of the line is the critical endpoint of the first order phase transition.The shaded region denotes the charged PC phase.Now if we go to the chiral limit (m π = 0) (the line µ at m π = 0) these results are in good agreement with the ones discussed above.There are two critical values of the chemical potential, µ c1 and µ c2 that corresponding to the onset and disappearance of charged PC, respectively.But increasing the value of the current quark mass the shaded region shrinks until these two points meet, µ c1 ≡ µ c2 .They meet at the value of pion mass m c π ∼ 9 MeV that corresponds to a current quark mass of m ∼ 10 keV, much less than the physical values, m ∼ 5 MeV.The (µ,µ Q ) phase portrait at the physical point is depicted at Figure 5 [20], where one can also see that the electric neutrality line does not cross the PC phase.Hence, one can infer that charged PC is an extremely fragile phenomenon with respect to the explicit chiral symmetry breaking effect (non-zero current quark mass).
In the papers [20,21] it was also shown that the equality between the electric chemical potential (in essence the isospin one) and the pion mass in medium, |µ Q | = m π , is a threshold to the charged PC phase provided the transition to the charged PC phase is of second order.And the absence of the charged PC can be understood in the following way.Note that increase of current quark mass leads to a drastic amplification of the vacuum pion mass because at small masses m pion mass depends on quark mass as m π ∝ √ m, whereas the change in the value of the quark mass in keV scale brings about no considerable modification to the value of electric chemical potential µ Q .
So let us note that not only PC d phase is absent at the charge neutrality line, but there is no PC at all and electrically neutral quark matter never meets the regions occupied by the charged PC phase.Thus, one can summarize that at the physical quark mass the NJL model consideration of neutral quark matter forbids the charged PC to take place in the ground state.

III. CONDITIONS PROMOTING THE APPEARANCE OF CHARGED PION CONDENSATION IN DENSE QUARK MATTER
In the previous Section 2 (see 2.4) it was shown that charged PC does not take place in neutral quark matter with isospin asymmetry.But recently there have been found conditions that can be realized in real physical situations and can lead to the appearance of charged PC in dense quark matter [34,[37][38][39][40][41].Some of these conditions are elaborated in the framework of the (3+1)-dimensional NJL model discussed in the previous section [40,41].Other physical constraints that may lead to a generation of the charged PC phase of dense quark matter are investigated, for simplicity, within the framework of a toy (1+1)-dimensional NJL model (NJL 2 ) [34,[37][38][39].So first of all, let us remind the basic properties of the NJL 2 model in the presence of baryon µ B ≡ 3µ and isospin µ I ≡ 2ν chemical potentials (see, e.g., in Ref. [24]).
A. The NJL2 model and its phase structure As it was noted above, the the NJL 2 model can be used in order to get quite common qualitative information about the properties of real dense quark matter with isospin asymmetry.Formally, its Lagrangian has the same form as the Lagrangian (1)-( 2) of the NJL model in (3+1)-dimensional spacetime.But u and d quark fields are now (in two dimensions) the two-component Dirac spinors and, moreover, the γ ν (ν = 0, 1) and γ 5 matrices are the 2×2 matrices in the two-dimensional spinor space, It is clear that the NJL 2 model is a direct generalization of the (1+1)-dimensional Gross-Neveu model [48] with a single flavor (i.e. single massless color N c -plet) to the case of two flavors of massive quarks.
To find the TDP of the NJL 2 model in the leading order of the large-N c expansion, one can use the technique of the previous section, where we have obtained the corresponding TDP of the (3+1)-dimensional NJL model.Namely, it is possible to introduce the auxiliary composite bosonic fields σ(x) and π a (x) as well as their ground state expectation values M − m and π a , respectively (see in Eqs. ( 6) and ( 12), for more details).Moreover, without loss of generality, it can be assumed that π 1 ≡ ∆ ≥ 0 and π 2,3 = 0. Then the expression Ω(M, ∆) for the unrenormalized TDP of the NJL 2 model can be easily shown to have the following form where It is possible to perform in Eq. ( 26) the integration with respect to p 0 , as in the (3+1)-dimensional case.Then, one can see that the TDP ( 26) is ultraviolet divergent, hence to get any physical information, one should first renormalize it (in contrast to the (3+1)-dimensional NJL model, the NJL 2 model is renormalizable).For this purpose one can use a special assumption that the bare coupling constant G depends on the cutoff parameter Λ (where Λ restricts the region of integration in the divergent integral of (26), |p 1 | < Λ).At zero bare quark mass, m = 0, the renormalization procedure was discussed in the literature, see, e.g., in Refs.[34,38,49].In short, it is evident in this case that cutting off the divergent integral in (26) and then using the following dependence of the bare coupling constant on the cutoff parameter, G ≡ G(Λ), one has a new free finite renormalization group invariant (does not depend on the cutoff Λ) massive parameter M 0 (and it can be shown that it is just a dynamically generated quark mass in vacuum, i.e. at µ = 0 and ν = 0).Then in the limit Λ → ∞ one can obtain a finite renormalization group invariant expression for the TDP.In particular, in vacuum, where we denote the TDP as V 0 (M, ∆), it looks like It can be shown that the global minimum of the TDP ( 28) is reached at M = M 0 and ∆ = 0, and chiral symmetry of the model is broken down.In the chiral limit (m = 0) the phase structure of NJL 2 model was studied in detals in Refs.[34,49], and it follows from this considerations that at arbitrary values of µ and ν the phase portrait of the NJL 2 model does not contain the charged PC phase with nonzero density of baryons.The same property is valid for the phase structure of the massive NJL 2 model, which does not predict the charged PC phenomenon in dense quark matter (see, e.g., in Refs.[24,37]).
Let us note that (1+1)-dimensional NJL model has the same chiral structure as the corresponding (3+1)-dimensional NJL model and it possesses a lot of common properties and features with real QCD.For example, QCD and NJL 2 models are both renormalizable, both have asymptotic freedom, dimensional transmutation and the spontaneous chiral symmetry breaking (in vacuum) [50][51][52][53].Moreover, they have similar µ B − T phase diagrams.Hence, NJL 2 model can be used as a laboratory for the qualitative simulation of specific features of QCD and, since it is renormalizable, it can be done at arbitrary energies in contrast to effective (3+1)-dimensional NJL model.
Let us now discuss the breaking of continuous symmetry and the famous no-go theorem [54] that, in general, forbids in two-dimensional space-time the spontaneous breaking of any continuous symmetry.Currently it is well understood (see, e.g., the discussion in [50,53,55]) that in the limit N c → ∞ (N c is the number of quark colors) the usual no-go theorem [54] does not work.This fact follows directly from the observation that at finite N c the quantum fluctuations, which at finite N c usually washes away a long range order corresponding to a spontaneous breakdown of symmetry, are suppressed by the powers of 1/N c factor in the large-N c limit.Thus, it is possible to simulate the phenomena taking place in real hadronic (quark) matter, such as spontaneous breaking of continuous chiral symmetry or charged pion condensation (spontaneous breaking of the continuous isospin symmetry), in the framework of a comparatively simple (1+1)-dimensional NJL-type models, although only in the limit of the large N c .

B. Finite size effect and nontrivial topology
Now let us return to the factors that can generate the appearance of charged PC in dense quark matter.Physical phenomena, in general, usually occur in a restricted space.For instance, the quark gluon matter droplets produced in heavy ion collisions (for example, at RHIC) always has a finite size.A volume of the droplets can range and can be as large as of order of several dozens or even hundred of fm, while size of the smallest system produced is estimated to be as low as 2 fm [56].So it is interesting to account for the influence of the finite size effects on (µ,µ I ) phase diagram and in particular on charged PC and it has been done in [34].The analysis has been performed in the framework of NJL 2 model.In order to study the influence of finite size effects on the phase structure of quark matter, one needs to assume from the beginning some boundary conditions.In finite temperature field theory one can introduce temperature by replacing integral over the temporal direction to a summation of the Matsubara frequencies with the requirement that fermion field satisfies the antiperiodic boundary conditions.But in contrast to the temporal direction (thermal case), where boundary conditions are fixed just by the statistics, there is no such restriction in the case of spatial directions.So for the system with the finite size there is always important issue concerning the choice of the boundary conditions for the fermion fields.
In the papers on this subject, both antiperiodic and periodic boundary conditions have been used and neither of them has been decisively excluded.Let us emphasize what the most substantial difference between these two boundary conditions is.It is actually whether one includes zero-momentum mode contribution or not.The zero-mode contribution is taken into account in the periodic boundary condition, whereas it is excluded in the case of antiperiodic one.For bosons the most natural choice of boundary condition is periodic.However, it has been always not clear what boundary conditions one should use for fermions, for example, quarks.
In some cases, the antiperiodic spatial boundary condition is used to keep symmetry between the time and space directions or to get a consistent picture with the dependence of pion mass on the finite size obtained from chiral perturbation theory (ChPT).
Various boundary conditions induce opposite results on vacuum properties of QCD: the antiperiodic spatial boundary condition for quarks induces the chiral symmetry restoration in the small system (analogously to the thermal effects), while the periodic boundary condition enhances the chiral symmetry breaking in vacuum (similar to the catalysis of chiral symmetry breaking under the influence of magnetic field) [57][58][59].
The above discussion was devoted to the quark models such as NJL model but finite size effects have been also studied in the framework of different hadronic models such as ChPT, Walecka model, hadron resonance gas model etc.. Let us elaborate on that a little bit.In order to get the extrapolation to infinite volume limit for lattice QCD finite size effects have been investigated in the framework of ChPT ( [68] and references therein), where the periodic boundary conditions were imposed.While twisted boundary conditions (similar to the general case of the ( 29)) have been used in ChPT in [69].In terms of hadron resonance gas model, the finite size effects were investigated in [70] employing a lower momentum cutoff, and the multiple reflection expansion formalism has been used to study influence of finite size on the chemical freeze-out [71].Walecka model in a finite volume with antiperiodic spatial boundary conditions have been considered in [72] and both types of boundary conditions (periodic and antiperiodic) have been investigated in the framework of parity-doublet and Walecka model in [73].Other methods not directly involving boundary conditions were used in studies of finite size matter in nuclei [74,75].
To simulate the influence of finite size effect, in Ref. [34] the (1+1)-dimensional system with Lagrangian of the form (1)-( 2) was placed into a confined space region of the form 0 ≤ x ≤ L (where x is the space coordinate).It is well-known that such a constraint is tantamount to the investigation in a space-time with nontrivial topology, for example in our case, one of the space coordinates are compactified.It means that one can consider the model in space-time with R 1 × S 1 topology and quark quantum fields should satisfy the boundary conditions of the following form where 0 ≤ α ≤ 2 is the parameter fixing the boundary conditions, for example, α = 0 corresponds to periodic boundary condition and α = 1 to the antiperiodic one, the parameter L is the length of the circumference S 1 .As a result, in order to obtain the TDP Ω L (M, ∆) of the NJL 2 system with Lagrangian (1)-( 2) in the confined domain, one should make a standard replacement of the integration over momentum in the expression (26) to the summation according to the rule: Let us briefly discuss the analogy of parameter α with imaginary baryon chemical potential.For temporal dimension the Matsubara frequencies are p 0l = π β (2l+1)−iµ.Here the antiperiodic boundary conditions have been used.One can note that in the case of imaginary chemical potential the Matsubara frequencies are p 0l = π β (2l + 1 + Im(µ)β/π) and the temporal boundary condition can be transformed into the imaginary chemical potential divided by temperature µ/T .In recent years a great effort has been made in studying the phase diagram at imaginary chemical potential as in lattice QCD [35] as well as in effective models including NJL model ( [36] and references therein).
Let us briefly discuss the analogy of parameter α with imaginary baryon chemical potential.For temporal dimension the Matsubara frequencies are p 0l = π β (2l+1)−iµ.Here the antiperiodic boundary conditions have been used.One can note that in the case of imaginary chemical potential the Matsubara frequencies are p 0l = π β (2l + 1 + Im(µ)β/π) and the temporal boundary condition can be transformed into the imaginary chemical potential divided by temperature µ/T .In recent years a great effort has been made in studying the phase diagram at imaginary chemical potential as in lattice QCD [35] as well as in effective models including NJL model ( [36] and references therein).Now let us introduce the notations used in [34].Instead of dimensionfull quantities such as M , ∆, µ, ν (see section 3.1) and L it is convenient to use the dimensionless quantities: So λ is just the inverse length of the circumference S 1 and shows how large the volume of space, in which our system is confined, is.The case λ = 0 corresponds to the topology R 1 × R 1 (infinite space) and the case λ = 0 corresponds to R 1 × S 1 , where the length of compactified dimension is π αM0 .One can think in the following way: the larger λ corresponds to smaller volume and stronger finite size effects.
The influence of nontrivial topology on the phase structure of NJL 2 model has been studied in [34], where both periodic and antiperiodic boundary conditions have been studied.Let us present a couple of phase portraits.For example, in Figures 6 and 7 the (ν,µ) phase diagrams are presented for various values of λ in the case of periodic boundary conditions.And one can see that there is a region of the PC d phase at both plots, it spans around the small values of isospin chemical potential.One can also notice that PC d phase is presented in a wide interval of λ and it is enlarged if the system becomes smaller.In the case of antiperiodic boundary conditions this phase has not been found and it is probably natural, as in this case the finite size effect is similar to the effect of the finite temperature on the system (it usually restores all the symmetries and leads to the vanishing condensates).But PC d phase has been found in quark matter system of a finite size at least for periodic boundary conditions, so probably one can conclude that finite size and non-trivial topology effects might lead to the appearance of the charged PC phase with nonzero baryon density at least in the framework of the NJL 2 model.It would be interesting to consider this situation in the framework of more realistic (3+1)-dimensional model approach, but unfortunately, to the best of our knowledge, this consideration has not been performed up to now.Now let us mention another interesting feature of the phase structure characteristic to the finite size systems at zero temperature.One can see at Figures 6 and 7 that there are small spikes at the border line between PC d and symmetrical phases.At the first glance one can think that they are numerical calculations artifacts.But they are actually physical and comes from the fact that in the finite systems the momentum spectrum became discrete and when one changes external parameters new momentum eigenvalues go into play and this can lead to the non-analyticity of critical lines, condensates and some thermodynamical quantities.Similar oscillations were observed in [76].

C. Inhomogeneous pion condensation in dense baryonic matter
Now let us turn our attention and discuss another factor that can lead to the generation of charged PC in dense quark matter, namely, the possibility of inhomogeneity of condensates i. e. dependence of condensates on the spacial space-time dimensions.The phase structure of dense quark matter with isospin asymmetry has been studied in the framework of NJL 2 model with possibility of inhomogeneous charged pion condensate [37].The consideration was performed at the physical point, i. e. for physical values of current quark mass.
Let us say a couple of words about why one should consider inhomogeneous condensates.In vacuum, i. e. in empty space without particle density and, equivalently, at zero values of the chemical potentials µ and µ I , the vacuum expectation values σ(x) and π a (x) do not depend on space-time coordinates.Although, in dense quark medium or, in other words, at non-zero values of the chemical potentials µ = 0, µ I = 0, the ground state expectation values of bosonic fields might assume a nontrivial dependence on spatial coordinate x.For example, in the considerations of Ref. [37] there has been used the following ansatz where M, b, and ∆ are quantities that do not depend on space-time coordinates.In [37] the phase structure in this approach was explored and it was found that in the inhomogeneous case there is charged PC phase with non-zero baryon density at the phase portrait and it is quite extensive.At Figure 8 the (ν, µ) phase portrait of the model NJL 2 model is shown.There the phase IPC d are inhomogeneous phases with charged PC and with non-zero baryon density, and one can see that it occupies the large part of the phase diagram.
So one can conclude that spatial inhomogeneity of charged pion condensate is also a factor promoting the appearance of PC d phases at the phase diagram of dense quark matter.

D. Chiral imbalanced dense quark matter
Now let us return to the (3+1)-dimensional NJL (1) model and talk about another factor that can generate charged pion condensation in dense quark matter.Besides non-zero temperature and baryon density (nonzero baryon chemical potential), there exist additional parameters that may be relevant to the properties of the QCD phase diagram.We have already talked about isotopic chemical potential.There is another captivating phenomenon that fell into spotlight quite recently, it is the chiral imbalance (different average densities of right-(n R ) and left-handed (n L ) quarks).This phenomenon is a remarkable effect originating from highly non-trivial interplay of different features of QCD such as chiral symmetry, chiral (axial) anomaly, and possible gluon configuration topology.It is expected to appear in the processes in heavy ion collisions violating P and CP on the event-by-event basis [44] and possibly may lead to interesting phenomena such as chiral magnetic effect [43].Chiral imbalance might be also realized in interiors of compact stars [42].In addition, chiral media, i. e. media with non-zero chiral imbalance, can be attained in condensed matter systems such as Dirac and Weyl semimetals [45,46].Also, note that phenomena regarding a chiral imbalance in the context of QCD are usually studied in the framework of NJL models with a chiral chemical potential [77].
The chiral imbalance of the system is described by the chiral charge density n 5 = n R − n L but instead of it one can use the corresponding chiral chemical potential µ 5 in the grand canonical ensemble approach.In general, the chiral charge n 5 is a more relevant quantity, but due to technical reasons when studying the phase diagram, it is easier to work with chemical potential µ 5 , and we will follow this way in this paper.Let us also note that the chiral charge n 5 is not a strictly conserved quantity, because of non-zero quark condensate as well as quantum chiral anomaly.Therefore, µ 5 chemical potential is not a quantity conjugated to a strictly conserved charge.But it is possible to treat µ 5 as the chemical potential describing a system in thermodynamical equilibrium in the state with nonzero value of n 5 , although, on a larger time scale than the one needed for chirality changing processes to occur.Moreover, let us talk about other possibility and note that n 5 can be introduced for individual quark flavor, n u5 = n uR − n uL and n d5 = n dR − n dL , and it is evident that n 5 = n u5 + n d5 .Then it is possible to consider chiral isospin charge n I5 = (n u5 − n d5 ) as well as the corresponding chiral isospin chemical potential µ I5 .Unlike the chiral charge n 5 , the chiral isospin charge n I5 turns out to be a conserved quantity at least in the massless NJL model.Since gluons interact with different quark flavors in exactly the same way, it is usually assumed that in QCD matter the chiral charges n u5 and n d5 are equal and therefore n I5 = 0 and µ I5 = 0.But the chiral isospin chemical potential µ I5 can be generated due to chiral separation effect in dense quark matter.
In [40][41][42] the phase diagram of dense quark matter both with chiral and isospin imbalances has been considered and the influence of chiral µ 5 and chiral isospin µ I5 chemical potentials on the phase structure and charged PC, in particular, has been studied.The investigations have been performed in the framework of two flavored (3+1)dimensional massless NJL model with baryon µ B , isospin µ I , chiral µ 5 and chiral isospin µ I5 chemical potentials.Its Lagrangian has the form where L N JL and other notations are presented in Eq. ( 1).In the following we use the notations µ ≡ µ B /3, ν = µ I /2 and ν 5 = µ I5 /2.The phase diagram of NJL model with baryon µ B , isospin µ I , and chiral isospin µ I5 chemical potentials (at µ 5 = 0) has been considered in [40] and the (ν, ν 5 ) phase portrait of the model is shown at Figure 9.It is clear from this figure that chiral imbalance (chiral isospin chemical potential µ I5 ) leads to the generation of charged PC in dense quark matter with isotopic asymmetry.But let us note that for the PC d phase to be generated the matter should have non-zero isospin imbalance (isospin chemical potential µ I ).One can also observe that in contrast to the first prediction of PC d phase in dense quark matter with isotopic asymmetry (see Sections 2.2, 2.4 and Figures 1, 3), the generation is quite extensive and the region of PC d phase is quite large (not marginal in any way).
Then in [41] the influence of chiral chemical potential µ 5 was studied in addition, and it was shown that chiral chemical potential µ 5 generates charged PC phase in dense baryonic (quark) matter as well.It is easily seen from Figure 10, where the (ν, µ 5 ) phase portrait of the NJL model is presented at ν 5 = 0.But this generation is not that vast and occurs at not so large baryon densities.The case of the joint effect of two chiral chemical potentials, µ 5 = 0 and ν 5 = 0, on the phase structure of the NJL model was also considered in this work.The (ν, ν 5 )-phase diagram (at µ 5 = 0) is depicted at Figure 11  phase.Moreover, this phase can be generated even at zero values of isospin chemical potential µ I .The fact that there is no need for isospin imbalance in order to generate the PC d phase is clearly presented at Figure 12, where the (ν 5 , µ 5 )-phase diagram is depicted at zero value of isospin chemical potential, ν = 0, and it exhibits a large region containing PC d phase.So one can conclude that chiral and chiral isospin chemical potentials together generate PC d phase very efficiently.Summarizing, one can say that chiral imbalance is a factor that leads to the generation of charged PC in dense quark matter.

IV. SUMMARY AND CONCLUSIONS
In this short review we tried to give an outline of investigations of charged PC in dense baryonic (quark) matter in the framework of effective (3+1)-dimensional NJL model and toy QCD motivated (1+1)-dimensional NJL model (NJL 2 model).
The possibility of charged PC phase with non-zero baryon density was shown for the first time in the framework of the massless NJL model in Ref. [8] and then in Ref. [9] it was found that this phase persists even in electrically neutral dense baryonic (quark) matter.But later it was shown in Refs.[17,[20][21][22][23] that charged PC condensation in neutral dense quark matter is enormously frail with respect to explicit chiral symmetry breaking effect (non-zero current quark mass), and it was shown, in particular, that it is forbidden for the physical values of current quark masses.
Nevertheless, recently there have been found factors and conditions that can be realized in real physical systems and can promote the appearance of charged PC phenomenon in dense baryonic (quark) matter.Namely, in the paper [34] it was shown that if one includes into consideration the fact that system can have finite size, then charged PC phase with non-zero baryon density can be realized in the system.It was also revealed in [37] that the possibility of spatially inhomogeneous charged PC condensate allows this phase to appear as well.And more recently [38][39][40][41] it was found that there is another interesting factor that can induce a charged PC phase in dense baryonic matter, it is a chiral imbalance of the system (nonzero difference between densities of left-handed and right-handed quarks).
The finite size effect and spatially inhomogeneous charged PC have been considered in Refs.[34,37] in the framework of a toy QCD motivated NJL 2 model, and the influence of chiral imbalance has been investigated both in the framework of NJL 2 model [38,39] and in more realistic (3+1)-dimensional NJL model [40,41].These models allow one to consider the region of non-zero baryon densities of QCD phase diagram, which up to now cannot be studied in lattice QCD (simulations are almost impossible due to the infamous sign problem).For simplicity, the whole review is devoted to the consideration of the case of zero temperature only.However, our recent investigation [42], made in the framework of the NJL 4 model with non-zero bare quark masses, shows that in chirally asymmetric dense quark matter charged PC phenomenon can be realized even at sufficiently high temperatures.
Let us once more recapitulate the conditions that can lead to the generation of charged PC phenomenon in dense baryonic matter: 1) finite size effects and non-trivial topology, 2) possibility of inhomogeneous charged pion condensates, 3) chiral imbalance.
These results can be interesting in the context of heavy ion collision experiments such as NICA and FAIR, where it is expected to get high baryon densities.It is of interest also in the neutron star physics, since quark matter might be produced in their cores, where very high baryon and isospin densities are attained.

FIG. 1 .
FIG.1.The (ν, µ) phase portrait in the chiral limit in the framework of NJL model.PC phase is charged pion condensation phase.Subscript d means that in this phase baryon density is nonzero.The sym denotes the symmetric phase.

FIG. 2 .
FIG.2.The value of charged pion condensate ∆0 as a function of quark number chemical potential µ in electric charge neutral quark matter.

FIG. 3 .FIG. 4 .FIG. 5 .
FIG.3.The (ν, µ) phase portrait of the model at the physical point (m = 0).CSB phase denotes chiral symmetry breaking phase and the notation ApprSYM means the phase where dynamical quark mass is close to the current quark mass value M ≈ m, i.e. the chiral symmetry is approximately restored in this phase.There is no a real phase transition between CSB and ApprSYM phases but a continuous crossover.Other notations are the same as in Figure1.

FIG. 8 .
FIG. 8.The (ν, μ) phase portrait of the model with the possibility of inhomogeneous pion condensation.IPC stand for inhomogeneous charged PC phases.All other notations are the same.