Chiral imbalanced hot and dense quark matter: NJL analysis at the physical point and comparison with lattice QCD

TL;DR

This paper uses the NJL model to explore the phase structure of hot, dense quark matter with isospin and chiral imbalances, revealing dualities and predicting phase transitions consistent with lattice QCD results.

Contribution

It demonstrates an approximate duality between chiral symmetry breaking and pion condensation at the physical point and compares NJL predictions with lattice QCD, supporting the duality hypothesis.

Findings

Chiral isospin chemical potential induces pion condensation at high temperatures.

Pseudo-critical temperature increases with chiral chemical potential.

Lattice QCD results support the NJL model's duality predictions.

Abstract

Hot and dense quark matter with isospin and chiral imbalances is investigated in the framework of the (3+1)-dimensional Nambu--Jona-Lasinio model (NJL) in the large-$N_c$ limit ($N_c$ is the number of quark colors). Its phase structure is considered in terms of baryon -- $\mu_B$, isospin -- $\mu_I$ and chiral isospin -- $\mu_{I5}$ chemical potentials. It is shown in the paper that (i) in the chiral limit there is a duality between chiral symmetry breaking (CSB) and charged pion condensation (PC) phenomena. (ii) At the physical point, i.e. at nonzero bare quark mass $m_0$, and temperature this duality relation is only approximate, although rather accurate. (iii) We have shown that the chiral isospin chemical potential $\mu_{I5}$ in dense quark matter generates charged pion condensation both at zero and nonzero $m_0$, and at $\mu_{I5}\ne 0$ this phase might be observed up to temperatures as high as 100 MeV. (iv) Pseudo-critical temperature of the chiral crossover transition rises in the NJL model with increasing $\mu_{I5}$. (v) It has been found an agreement between particular sections of the phase diagram in the framework of NJL model and corresponding ones in lattice QCD simulations. Two different plots from different lattice simulations that are completely independent and are not connected at the first sight are in reality dual to each other, it means that lattice QCD simulations support the hypothesis that in real quark matter there exists the (approximate) duality between CSB and charged PC. Moreover, we can reverse the logic and we can predict the increase of pseudo-critical temperature with chiral chemical potential, the much debated effect recently, just by the duality notion, hence bolster confidence in this result (lattice QCD showed this feature for unphysically large pion mass) and put it on the considerably more solid ground.

Equations

$$\displaystyle L=\bar{q}\Big{[}\gamma^{\nu}\mathrm{i}\partial_{\nu}-m_{0}+\frac{\mu_{B}}{3}\gamma^{0}+\frac{\mu_{I}}{2}\tau_{3}\gamma^{0}\Big{]}q+\frac{G}{N_{c}}\Big{[}(\bar{q}q)^{2}+(\bar{q}\mathrm{i}\gamma^{5}\vec{\tau}q)^{2}\Big{]}.$$

$$U_B(1):q\to \exp (\mathrm{i}\alpha /3)q;U_{I_3}(1):q\to \exp (\mathrm{i}\alpha {\tau }_3/2)q;U_{A{I}_3}(1):q\to \exp (\mathrm{i}\alpha {\gamma }^5{\tau }_3/2)q.$$

$${\hat{n}}_B=\frac{1}{3}\left({\hat{n}}_u+{\hat{n}}_d\right),{\hat{n}}_I=\frac{1}{2}\left({\hat{n}}_u-{\hat{n}}_d\right).$$

$${\hat{n}}_{I5}=\frac{1}{2}\left({\hat{n}}_{uR}-{\hat{n}}_{uL}-{\hat{n}}_{dR}+{\hat{n}}_{dL}\right)=\frac{1}{2}\left({\hat{n}}_{u5}-{\hat{n}}_{d5}\right),$$

$$\bar{L}$$

$$\displaystyle{\cal L}\displaystyle=\bar{q}\Big{[}\gamma^{\rho}\mathrm{i}\partial_{\rho}-m_{0}+\mu\gamma^{0}+\nu\tau_{3}\gamma^{0}+\nu_{5}\tau_{3}\gamma^{0}\gamma^{5}-\sigma-\mathrm{i}\gamma^{5}\pi_{a}\tau_{a}\Big{]}q-\frac{N_{c}}{4G}\Big{[}\sigma\sigma+\pi_{a}\pi_{a}\Big{]}.$$

$$\sigma (x)=-2\frac{G}{N_c}(\bar{q}q);{\pi }_a(x)=-2\frac{G}{N_c}(\bar{q}\mathrm{i}{\gamma }^5{\tau }_{a}q).$$

$$U_{I_3}(1):$$

$$U_{A{I}_3}(1):$$

$$\exp(\mathrm{i}{\cal S}_{\rm{eff}}(\sigma,\pi_{a}))=N^{\prime}\int[d\bar{q}][dq]\exp\Bigl{(}\mathrm{i}\int{\cal L}\,d^{4}x\Bigr{)},$$

$${\mathcal{S}}_{\mathrm{eff}}(\sigma (x),{\pi }_a(x))=-N_c\int d^{4}x\left[\frac{{\sigma }^2+{\pi }_a^2}{4G}\right]+{\tilde{\mathcal{S}}}_{\mathrm{eff}},$$

$$\exp (\mathrm{i}{\tilde{\mathcal{S}}}_{\mathrm{eff}})$$

$$D\equiv {\gamma }^{\nu }\mathrm{i}{\partial }_{\nu }-m_0+\mu {\gamma }^0+\nu {\tau }_3{\gamma }^0+{\nu }_5{\tau }_3{\gamma }^0{\gamma }^5-\sigma (x)-\mathrm{i}{\gamma }^5{\pi }_a(x){\tau }_a,$$

$${\mathcal{S}}_{\mathrm{eff}}(\sigma (x),{\pi }_a(x))=-N_c\int d^{4}x\left[\frac{{\sigma }^2(x)+{\pi }_a^2(x)}{4G}\right]-\mathrm{i}{N}_c{\mathrm{Tr}}_{sfx}\ln D,$$

$$\frac{\delta {\mathcal{S}}_{\mathrm{eff}}}{\delta \sigma (x)}=0,\frac{\delta {\mathcal{S}}_{\mathrm{eff}}}{\delta {\pi }_a(x)}=0,$$

$$⟨\sigma (x)⟩\equiv \sigma ,⟨{\pi }_a(x)⟩\equiv {\pi }_a,$$

$$\int d^{4}x\Omega(\sigma,\pi_{a})=-\frac{1}{N_{c}}{\cal S}_{\rm{eff}}\big{(}\sigma(x),\pi_{a}(x)\big{)}\Big{|}_{\sigma(x)=\sigma,\pi_{a}(x)=\pi_{a}}.$$

$$⟨\sigma (x)⟩=M-m_0,⟨{\pi }_1(x)⟩=\Delta ,⟨{\pi }_2(x)⟩=0,⟨{\pi }_3(x)⟩=0.$$

$$\Omega (M,\Delta )$$

$$\overline{D}(p)=\not\!p+\mu\gamma^{0}+\nu\tau_{3}\gamma^{0}+\nu_{5}\tau_{3}\gamma^{0}\gamma^{5}-M-\mathrm{i}\gamma^{5}\Delta\tau_{1}\equiv\left(\begin{array}[]{cc}A~{},&U\\ V~{},&B\end{array}\right)$$

$$A=/p+\mu {\gamma }^0+\nu {\gamma }^0+{\nu }_5{\gamma }^0{\gamma }^5-M;B=/p+\mu {\gamma }^0-\nu {\gamma }^0-{\nu }_5{\gamma }^0{\gamma }^5-M;U=V=-\mathrm{i}{\gamma }^5\Delta ,$$

$$\displaystyle\mathop{\rm Det}\nolimits\overline{D}(p)\equiv\det\left(\begin{array}[]{cc}A~{},&U\\ V~{},&B\end{array}\right)=\det[-VU+VAV^{-1}B]=\det[BA-BUB^{-1}V]$$

$$\displaystyle\mathop{\rm Det}\nolimits\overline{D}(p)=\big{(}\eta^{4}-2a\eta^{2}-b\eta+c\big{)}\big{(}\eta^{4}-2a\eta^{2}+b\eta+c\big{)}\equiv P_{-}(p_{0})P_{+}(p_{0}),$$

$$a$$

$$c$$

$$\mathcal{D}:M⟷\Delta ,\nu ⟷{\nu }_5.$$

$${\eta }_1=\frac{1}{2}\left(-\sqrt{r^2-4q}-r\right),$$

$${\eta }_3=\frac{1}{2}\left(r-\sqrt{r^2-4s}\right),$$

$${\eta }_5=\frac{1}{2}\left(-\sqrt{r^2-4s}-r\right)=-{\eta }_4,$$

$${\eta }_7=\frac{1}{2}\left(r-\sqrt{r^2-4q}\right)=-{\eta }_2,$$

$$\mathrm{Det}\overline{D}(p)={\Pi }_{i=1}^8(\eta -{\eta }_i),$$

$$\Omega (M,\Delta )=\frac{(M-m_0{)}^2+{\Delta }^2}{4G}+\mathrm{i}\sum _{i=1}^8\int \frac{d^{4}p}{(2\pi {)}^4}\ln (p_0+\mu -{\eta }_i).$$

Full text

I Introduction

At normal (Earth) conditions, protons and neutrons form atomic nuclei, and the latter, together with their orbital electrons, form the ordinary matter of our environment. If matter is subjected to extreme compression, eventually all chemical and nuclear bonds are broken, and the matter is squeezed from the molecular scale to the sub-particle scale with density higher than 0.15 baryon per $\rm{fm^{3}}$. Experimental creation of such dense matter is a very hard problem but such conditions can take place inside compact stars due to compression by gravity into a stable and extremely dense state. As a rule neutron stars have comparatively low temperatures and one can assume that it is zero. What happens at high temperature is hard to probe studying the physics of neutron stars. Nevertheless, due to technology advances, modern accelerators of elementary particles are now able to collide not only single high energy protons, but also heavy ions consisting of many coupled protons and neutrons. It is believed that in the fireball just after heavy-ion collisions there emerges a droplet of quark gluon plasma with very high temperature. Physics of heavy ion collision experiments can shed some light on the conditions that existed a few microseconds after the Big Bang and provide answers to several other questions.

The fundamental theory of matter in such extreme conditions is quantum chromodynamics (QCD) which is a gauge field theory associated with $SU(3)$ group, where gauge bosons (gluons) play the role of interaction carriers of quarks. The main method of QCD analysis is the perturbative technique on the basis of coupling constant. However, it is not always possible to use this technique, as QCD calculations can be too complex or in the low energy region when the coupling constant is too large. In particular, QCD perturbative technique is not applicable in a consideration of physically reachable dense matter, etc. In these cases, non-perturbative methods, such as effective theories or lattice calculations, are usually used.

Thus, the entire QCD phase diagram could not be described currently in the framework of a unified theory. Lattice calculations are very useful for description of the region of zero density and high temperature. However, the so-called sign problem still presents insurmountable difficulties for lattice calculations in the nonzero density region. On the other hand effective theories do not have fundamental background and as a result do not share the main prominent features with QCD such as a gauge invariance, renormalizability, etc. Nevertheless, at this moment, effective models are the best tool for investigating dense quark matter. At this time one of the most widely used effective model is the Nambu–Jona-Lasinio (NJL) model Nambu:1961fr ; Klevansky:1992qe ; Hatsuda:1994pi ; Buballa:2003qv .

It is well known that usually dense baryonic matter in compact stars obeys an isospin asymmetry, i.e. where the densities of up- and down quarks are different (it is characterised by isospin chemical potential $\mu_{I}$). In experiments on heavy-ion collisions, we also have to deal with quark matter which has an evident isospin asymmetry because of different neutron and proton contents of colliding ions. In early 70-th Sawyer Sawyer:1972cq and independently Migdal Migdal:1973 have shown that there might be phase transition from pure neutron matter to mixed hadron matter with protons, neutrons and $\pi^{0}$-pions at superdense matter in the compact stars. Later, using the chiral perturbation theory, it was shown that there is a threshold $\mu_{I}^{c}=m_{\pi}\approx 140MeV$ of a phase transition to the charged pion condensation (PC) phase Son:2000 . This result was ultimately proved in the framework of random matrix model Klein:2003 , Ladder-QCD model Barducci , resonance gas model Toublan , quark-meson model Fraga , NJL model ek1 (including (1+1)-dimensional version of the NJL model ek2 ) and lattice simulations Gupta . Nevertheless, the whole picture is still a matter of debate.

Now the main question is whether the charged pion condensation exists in the real world and how this phenomenon behaves under influence of various external factors. And different factors can have a completely different effect on this phase. For example, in the framework of NJL model the finite-size effects, spatial inhomogeneity of the pion condensate Khunjua or chromomagnetic background field Fedotov could promote the charged PC phase. On the other hand, if the electric charge neutrality and $\beta$-equilibrium constraints are imposed, the charged PC phenomenon in quark matter depends strongly on the bare (current) quark mass values. In particular, it turns out that the charged PC phase with nonzero baryonic density is not realized within NJL models, if the bare quark mass $m_{0}$ reaches the physically acceptable values of $5\div 10$ MeV abuki , i.e. at the physical point. In addition, temperature $T$ and different model parameters such as coupling constants, etc, as well strongly influence on this phase Ebert . It is also worth to note that the phase structure of the isospin imbalanced quark matter below the threshold ($\mu_{I}<m_{\pi}$) is an important question because even small nonzero $\mu_{I}$ could double the critical endpoint of a phase diagram and affects the results of heavy-ion collision experiments Klein:2003 ; Kogut .

Recently, it has been shown in the framework of the massless (3+1)-dimensional NJL model (and in the leading large-$N_{c}$ order, where $N_{c}$ is the number of colors of quarks) that chiral imbalance promotes charged PC phase in dense matter at zero temperature Khunjua:2017mkc ; Khunjua:2018sro and responsible for the existence of the duality between chiral symmetry breaking (CSB) and charged PC phases. The imbalance between densities of left-handed and right-handed quarks (chiral imbalance) is a highly anticipated phenomenon that could occur both in compact stars and heavy ion collisions. This effect could stem from nontrivial interplay of axial anomaly and the topology of gluon configurations.111It is predicted that there is an electrical current in the chiral imbalanced quark matter under strong magnetic field Fukushima:2008xe . This phenomenon is called chiral magnetic effect and it could be an evidence of the chiral imbalance in QCD. Also, there is another mechanism of its origin – chiral separation effect which can be realized in dense matter in the presence of a strong magnetic field. In this case left-handed and right-handed quarks tend to move in opposite directions along the magnetic field, thereby creating regions with chiral imbalance. Moreover, in the case of two-flavored quark matter the chiral separation effect could promote (see below in Appendix A) both nonzero chiral density $n_{5}$ and nonzero isotopic chiral density $n_{I5}$, and quark matter can be described using the corresponding chemical potentials $\mu_{5}$ and $\mu_{I5}$.

It was already mentioned above that nonzero bare quark mass $m_{0}$ and nonzero temperature $T$ could destroy charged PC phase in the physically adequate circumstances. So one of the aims of our present work is to check the robustness of the charged PC phase generated by chiral imbalance under the influence of these destructive factors. Another purpose is to study in the framework of the NJL4 model the fate of the duality observed in the chiral limit Khunjua:2017mkc (where it is an exact symmetry) between CSB and charged PC phenomena in the leading large-$N_{c}$ order: we investigate the influence of the bare quark mass and temperature on this effect, etc. In particular, it is shown in our paper that duality correspondence between CSB and charged PC still is a very good approximate symmetry of a phase portrait of the NJL4 model even at $m_{0}\neq 0$ and $T\neq 0$.

It is interesting to investigate not only the charged PC phase but also hot quark matter itself with chiral asymmetry only. In this case at zero baryon chemical potential, $\mu_{B}=0$, there is no sign-problem and we have solid results from lattice simulations Braguta:2015zta ; Braguta:2015owi . Nevertheless, some key properties of chirally imbalanced quark matter are still under debate. So, in addition to charged PC phase, in the present paper we also investigate in the framework of the NJL4 model at $m_{0}\neq 0$ the dependence of the (pseudo-)critical temperature, which characterizes the chiral cross-over region of the phase diagram, on the chiral isospin chemical potential $\mu_{I5}$ and compare our results with other effective model investigations and lattice simulations on this topic. Note also that at $\mu_{B}=0$ and $m_{0}\neq 0$ the $(\mu_{I},T)$- and $(\mu_{5},T)$-phase diagrams have been obtained both using lattice QCD simulations and in the framework of the NJL model, and the results are in good agreement. Moreover, in the present paper we show that just these phase diagrams are dually conjugated (with a good precision) to each other, so there is a good reason to argue that duality between CSB and charged PC phenomena is confirmed by lattice QCD calculations.

The paper is organized as follows. In Sec. II a (3+1)-dimensional NJL model with two massive quark flavors ($u$ and $d$ quarks) that includes three kinds of chemical potentials, $\mu_{B},\mu_{I},\mu_{I5}$, is introduced. Furthermore, the symmetries of the model are discussed and its thermodynamic potential is presented in the leading order of the large-$N_{c}$ expansion both at zero and nonzero temperature $T$. In particular, it is shown in this section that in the chiral limit ($m_{0}=0$) the phase structure of the model (in the leading order over $1/N_{c}$) has a dual symmetry between CSB and charged PC phenomena. In the next section we formulate the main consequences of the exact dual symmetry (Sec. III A), using which it is possible to decide that dual symmetry is performed approximately in the NJL model at $m_{0}\neq 0$ and $T=0$, but with good accuracy (Sec. III B). It Sec. III C we show that at nonzero values of the chiral isospin chemical potential $\nu_{I5}$ the charged PC phase with nonzero quark density can be realized in the model up to rather high values of temperature, $T\approx 100$ MeV. Moreover, here we show that duality is also fulfilled approximately at $T\neq 0$. In Sec. III D the plot of the pseudo-critical temperature of the chiral crossover transition as a function of $\mu_{I5}$ at $\mu=\mu_{I}=\mu_{5}=0$ is obtained. Here it is compared with results of other effective models and lattice QCD approaches. Sec. IV presents summary and discussion leading to the conclusion that duality between CSB and charged PC observed in the NJL4 model is supported by some phase diagrams obtained by lattice QCD simulations at $\mu_{B}=0$. Some technical details and issues not directly related to this work are relegated to Appendices A and B.

II The model and its thermodynamic potential

II.1 Lagrangian and symmetries

It is well known that in the framework of effective four-fermion field theories dense and isotopically asymmetric quark matter, composed of $u$ and $d$ quarks, can be described by the following (3+1)-dimensional NJL Lagrangian

[TABLE]

Here $q$ is a flavor doublet, $q=(q_{u},q_{d})^{T}$, where $q_{u}$ and $q_{d}$ are four-component Dirac spinors as well as color $N_{c}$-plets of the $u$ and $d$ quark fields, respectively (the summation in Eq. (1) over flavor, color, and spinor indices is implied); $\tau_{k}$ ($k=1,2,3$) are Pauli matrices; $m_{0}$ is the bare quark mass (for simplicity, we assume that $u$ and $d$ quarks have the same mass); $\mu_{B}$ and $\mu_{I}$ are chemical potentials which are introduced in order to study quark matter with nonzero baryon and isospin densities, respectively.

The symmetries of the Lagrangian (1) depends essentially on wether the bare quark mass $m_{0}$ and chemical potentials take zero or nonzero values. For example, in the most particular case, when $m_{0}=\mu_{I}=0$ the Lagrangian (1) is invariant under transformations from chiral $SU(2)_{L}\times SU(2)_{R}$ group, which is also inherent in 2-flavor QCD in the chiral limit. This symmetry is reduced to $U_{B}(1)\times U_{I_{3}}(1)\times U_{AI_{3}}(1)$ group if all chemical potentials are nonzero, and $m_{0}=0$. In this case the abelian baryon $U_{B}(1)$, isospin $U_{I_{3}}(1)$ and chiral isospin $U_{AI_{3}}(1)$ subgroups act on flavor doublet $q$ in the following way

[TABLE]

As a result, we see that in the chiral limit ($m_{0}=0$) the quantities $\hat{n}_{B}\equiv\bar{q}\gamma^{0}q/3$, $\hat{n}_{I}\equiv\bar{q}\gamma^{0}\tau^{3}q/2$ and $\hat{n}_{I5}\equiv\bar{q}\gamma^{0}\gamma^{5}\tau^{3}q/2$ are the density operators of the conserved baryon, isospin and chiral isospin charges of the system (1), respectively. Introducing the particle density operators for $u$ and $d$ quarks, $\hat{n}_{u}\equiv q_{u}\gamma^{0}q_{u}$ and $\hat{n}_{d}\equiv q_{d}\gamma^{0}q_{d}$, we have

[TABLE]

One can also introduce the particle density operators $\hat{n}_{fR}$ and $\hat{n}_{fL}$ for right- and left-handed quarks of each flavor $f=u,d$ (see in Appendix A). In this case the density operator of the chiral isospin charge looks like

[TABLE]

where the quantity $\hat{n}_{f5}\equiv\hat{n}_{fR}-\hat{n}_{fL}$ is usually called the density operator of the chiral charge for the quark flavor $f=u,d$. Below, in Appendix A, we discuss the possibility of the appearance of a nonzero chiral isotopic density in quark matter inside neutron stars. It can be explained on the basis of the chiral separation effect in the presence of a strong magnetic field in a dense baryonic medium.

However, at the physical point ($m_{0}\neq 0$) the symmetry of the Lagrangian (1) under transformations from axial isotopic group $U_{AI_{3}}(1)$ is explicitly broken. So in the most general case with $m_{0}\neq 0$, $\mu_{B}\neq 0$ and $\mu_{I}\neq 0$ the initial model (1) is invariant only under the $U_{B}(1)\times U_{I_{3}}(1)$ group. (We would like also to remark that Lagrangian (1) is invariant with respect to the electromagnetic $U_{Q}(1)$ group, $U_{Q}(1):~{}q\to\exp(\mathrm{i}Q\alpha)q$, at arbitrary values of $m_{0}$, where $Q={\rm diag}(2/3,-1/3)$.)

The ground state (the state of thermodynamic equilibrium) of quark matter with $n_{B}\neq 0$ and $n_{I}\neq 0$, where $n_{B}\equiv\langle\hat{n}_{B}\rangle$, $n_{I}\equiv\langle\hat{n}_{I}\rangle$,222The notation $\langle\hat{O}\rangle$ means the ground state expectation value of the operator $\hat{O}$. both at zero and nonzero values of $m_{0}$ has been investigated in the framework of the NJL model (1), e.g., in Refs. ek1 ; Ebert . However, the fact that quark matter may have a nonzero chiral isotopic charge was ignored in those papers. Recently, this gap in researches was filled in the paper Khunjua:2017mkc , where we have studied the properties of equilibrium quark matter at $n_{B}\neq 0$, $n_{I}\neq 0$ as well as at nonzero chiral isospin charge density $n_{I5}\equiv\langle\hat{n}_{I5}\rangle\neq 0$ in the framework of the massless (3+1)-dimensional two-flavor NJL model (temperature $T$ was taken to be zero in Ref. Khunjua:2017mkc ). In contrast to this, in the present paper we consider the properties of a more realistic quark matter, i.e. at $m_{0}\neq 0$ and $T\neq 0$, for which all densities $n_{B}$, $n_{I}$ and $n_{I5}$ are also nonzero. The solution of this problem can be most conveniently carried out in terms of chemical potentials $\mu_{B}$, $\mu_{I}$ and $\mu_{I5}$, which are the quantities, thermodynamically conjugated to corresponding charge densities $\hat{n}_{B}$, $\hat{n}_{I}$ and $\hat{n}_{I5}$ presented in Eqs. (3) and (4). Therefore, when solving this problem, one can rely on the Lagrangian of the form

[TABLE]

(Generally speaking, in this case the chiral isospin charge is no more a conserved quantity of our system. Therefore, chiral isospin chemical potential $\mu_{I5}$ is not conjugated to a strictly conserved charge. However, denoting by $\tau$ the typical time scale in which all chirality changing processes take place, one can treat $\mu_{I5}$ as the chemical potential that describes a system in thermodynamic equilibrium with a fixed value of $n_{I5}$ on a time scale much larger than $\tau$.)

Our goal is the investigation of the ground state properties (or phase structure) of the system, described by the Lagrangian (5), and its dependence on the chemical potentials $\mu_{B}$, $\mu_{I}$ and $\mu_{I5}$ (both at zero and nonzero temperature). It is well known that all information on the phase structure of the model is contained in its thermodynamic potential (TDP). Namely, in the behavior of its global minimum point vs. chemical potentials. Moreover, the values of charge densities $n_{B}\equiv\langle\hat{n}_{B}\rangle$, $n_{I}\equiv\langle\hat{n}_{I}\rangle$ and $n_{I5}\equiv\langle\hat{n}_{I5}\rangle$ in equilibrium quark matter can be found by differentiating the TDP in the global minimum point with respect to the corresponding chemical potentials $\mu_{B}$, $\mu_{I}$ and $\mu_{I5}$, etc. In order to find the TDP of the model, we start from a semibosonized version of the Lagrangian (5), which contains composite bosonic fields $\sigma(x)$ and $\pi_{a}(x)$:

[TABLE]

Here, $a=1,2,3$ and also we introduced the notations $\mu\equiv\mu_{B}/3$, $\nu\equiv\mu_{I}/2$ and $\nu_{5}\equiv\mu_{I5}/2$. From the auxiliary Lagrangian (6) one gets the equations for the bosonic fields:

[TABLE]

Note that the composite bosonic field $\pi_{3}(x)$ can be identified with the physical $\pi^{0}(x)$-meson field, whereas the physical $\pi^{\pm}(x)$-meson fields are the following combinations of the composite fields, $\pi^{\pm}(x)=(\pi_{1}(x)\mp i\pi_{2}(x))/\sqrt{2}$. Obviously, the semibosonized Lagrangian ${\cal L}$ is equivalent to the initial Lagrangian (5) when using the equations (7). Furthermore, the composite bosonic fields (7) change under the influence of transformations from the isospin $U_{I_{3}}(1)$ and axial isospin $U_{AI_{3}}(1)$ groups in the following manner:

[TABLE]

II.2 Thermodynamical potential. Zero temperature case.

Starting from the auxiliary Lagrangian (6), one obtains in the leading order of the large-$N_{c}$ expansion (i.e. in the one-fermion loop approximation) the following path integral expression for the effective action ${\cal S}_{\rm{eff}}(\sigma,\pi_{a})$ of the bosonic $\sigma(x)$ and $\pi_{a}(x)$ fields:

[TABLE]

where

[TABLE]

The quark contribution to the effective action, i.e. the term $\tilde{\cal S}_{\rm{eff}}$ in (9), is given by:

[TABLE]

where $N^{\prime}$ is a normalization constant. Moreover, in (10) we have introduced the notation $D$,

[TABLE]

for the Dirac operator, which acts in the flavor-, spinor- as well as coordinate spaces only. Using the general formula $\mathop{\rm Det}\nolimits D=\exp{\rm Tr}\ln D$, one obtains for the effective action (9) the following expression

[TABLE]

where the Tr-operation stands for the trace in spinor- ($s$), flavor- ($f$) as well as four-dimensional coordinate- ($x$) spaces, respectively.

The ground state expectation values $\langle\sigma(x)\rangle$ and $\langle\pi_{a}(x)\rangle$ of the composite bosonic fields are determined by the saddle point equations,

[TABLE]

where $a=1,2,3$. Just the knowledge of $\langle\sigma(x)\rangle$ and $\langle\pi_{a}(x)\rangle$ and, especially, of their behaviour vs. chemical potentials supplies us with a phase structure of the model. In the present paper we suppose that in the ground state of the system the quantities $\langle\sigma(x)\rangle$ and $\langle\pi_{a}(x)\rangle$ do not depend on spacetime coordinates $x$,

[TABLE]

where $\sigma$ and $\pi_{a}$ ($a=1,2,3$) are already spatially independent constant quantities. In fact, they are coordinates of the global minimum point of the thermodynamic potential (TDP) $\Omega(\sigma,\pi_{a})$. In the leading order of the large-$N_{c}$ expansion it is defined by the following expression:

[TABLE]

In what follows we are going to investigate the $\mu,\nu,\nu_{5}$-dependence of the global minimum point of the function $\Omega(\sigma,\pi_{a})$ vs $\sigma,\pi_{a}$. Let us note that in the chiral limit (due to a $U_{I_{3}}(1)\times U_{AI_{3}}(1)$ invariance of the model) the TDP (15) depends effectively only on the combinations $\sigma^{2}+\pi_{3}^{2}$ and $\pi_{1}^{2}+\pi_{2}^{2}$. Whereas at the physical point (i.e. at $m_{0}\neq 0$) it depends effectively on the combination $\pi_{1}^{2}+\pi_{2}^{2}$ as well as on $\sigma$ and $\pi_{3}$. Since in this case the relations $\langle\sigma(x)\rangle\neq 0$ and $\langle\pi_{3}(x)\rangle=0$ are always satisfied (see, e.g., in Ref. Fedotov ), at $m_{0}\neq 0$ one can put without loss of generality $\pi_{2}=\pi_{3}=0$ in Eq. (15), and study the TDP as a function of only two variables. For simplicity, we introduce the following $M\equiv\sigma+m_{0}$ and $\Delta\equiv\pi_{1}$ notations, and throughout the paper use the ansatz

[TABLE]

If in the global minimum point of the TDP we have $\Delta\neq 0$, then isospin $U_{I_{3}}(1)$ symmetry of the model is spontaneously broken down. Moreover, since at $m_{0}\neq 0$ chiral symmetry is explicitly broken down in the model, the $M$ coordinate of the global minimum is always a nonzero quantity. Note also that $M$ is a dynamical or constituent quark mass. In terms of $M$ and $\Delta$ the TDP (15) reads

[TABLE]

where

[TABLE]

is the momentum space representation of the Dirac operator $D$ (11) under the constraint (16). The quantities $A,B,U,V$ in Eq. (18) are really the following 4$\times$4 matrices,

[TABLE]

so the quantity $\overline{D}(p)$ from Eq. (18) is indeed a 8$\times$8 matrix whose determinant appears in the expression (17). Based on the following general relations

[TABLE]

and using any program of analytical calculations, one can find from Eqs. (19) and (22)

[TABLE]

where $\eta=p_{0}+\mu$, $|\vec{p}|=\sqrt{p_{1}^{2}+p_{2}^{2}+p_{3}^{2}}$ and

[TABLE]

It is evident from Eq. (24) that the TDP (17) is an even function over the variable $\Delta$, and parameters $\nu$ and $\nu_{5}$. In addition, it is invariant under the transformation $\mu\to-\mu$. 333Indeed, if simultaneously with $\mu\to-\mu$ we perform in the integral (17) the $p_{0}\to-p_{0}$ change of variables, then one can easily see that the expression (17) remains intact. Hence, without loss of generality we can consider in the following only $\mu\geq 0$, $\nu\geq 0$, $\nu_{5}\geq 0$, and $\Delta\geq 0$ values of these quantities. Moreover in the chiral limit, the TDP (17) is invariant with respect to the so-called duality transformation:

[TABLE]

(It is interesting to note that the dual symmetry (25) is also an inherent property of the TDP of the model (5) in the chiral limit and at $N_{c}\to\infty$, but in the (1+1)-dimensional spacetime 2dim .) One can find roots of the polynomials (23) analytically, the procedure is relegated to Appendix B. Four roots of $P_{+}(\eta)$ have the following form

[TABLE]

The roots of $P_{-}(\eta)$ can be obtained by changing $b\to-b$ (changing $b\to-b$ is equivalent to $q\leftrightarrow s$),

[TABLE]

where $q=\frac{1}{2}\left(-2a+r^{2}-\frac{b}{r}\right),\,\,\,s=\frac{1}{2}\left(-2a+r^{2}+\frac{b}{r}\right)$, and $r$ has quite complicated form, but could be always chosen as a real one (all the details can be found in Appendix B). As a result, we have from Eq. (23) that

[TABLE]

where each root $\eta_{i}$ is invariant with respect to the duality transformation (25). So, it is evident from Eqs. (17) and (23) that for the TDP one can obtain the following expression

[TABLE]

Then, taking in account a general formula

[TABLE]

and using the fact that each root $\eta_{i}$ of Eqs. (26) and (27) has a counterpart with opposite sign as well as the relation $|\mu-\eta_{i}|+|\mu+\eta_{i}|=2|\eta_{i}|+2\theta(\mu-|\eta_{i}|)(\mu-|\eta_{i}|)$, one gets

[TABLE]

To obtain the second line of Eq. (31), where $p\equiv|\vec{p}|$ and $\Lambda$ is a three-momentum cutoff parameter, we have integrated in the first line of it over angle variables. If we are interested in knowing the phase structure of the model at zero temperature, we should study just the TDP (31) vs $M$ and $\Delta$ on the global minimum point (GMP). It is clear that at $m_{0}\neq 0$ the GMP of the TDP has the form $(M_{0},\Delta_{0})$, where $M_{0}$ is always a nonzero quantity. If in this case $\Delta_{0}\neq 0$, then we are in the charged PC phase with spontaneous breaking of the isospin $U_{I_{3}}(1)$ symmetry.

II.3 Thermodynamical potential. Non-zero temperature case.

Though, the effect of non-zero temperatures is quite predictable (one can expect that the temperatures just restore all the broken symmetries of the model), here we include nonzero temperatures into consideration because it is important in a number of applications. In heavy ion collisions and early Universe the temperatures are huge and its account looks inevitable, but it even makes sense in other not so apparent situations. We know that compact stars are cold and one can consider their temperatures as zero. But probably there could be scenarios in which the temperatures could be important even in the context of compact stars. For example, their temperatures right after they are born in a supernova explosion can be as high as $T\approx 10$ MeV. So it is instructive to know how robust the charged PC phase under temperature.

To introduce finite temperature into consideration, it is very convenient to use the zero temperature expression (29) for the TDP. Then, to find the temperature dependent TDP $\Omega_{T}(M,\Delta)$ one should replace in Eq. (29) the integration over $p_{0}$ in favor of the summation over Matsubara frequencies $\omega_{n}$ by the rule

[TABLE]

In the expression obtained, it is possible to sum over Matsubara frequencies using the general formula (the corresponding technique is presented, e.g., in jacobs )

[TABLE]

where $\beta=1/T$. As a result, one can obtain the following expression for the TDP $\Omega_{T}(M,\Delta)$

[TABLE]

where $\Omega(M,\Delta)$ is the TDP (31) of the system at zero temperature. Since each root $\eta_{i}$ in Eq. (34) is a dually ${\cal D}$ invariant quantity (see in Eq. (25)), it is clear that in the chiral limit the temperature dependent TDP (34) is also symmetric with respect to the duality transformation ${\cal D}$.

Finally, it is necessary to note that in the framework of the NJL4 model the leading order of the large-$N_{c}$ limit is identical to the mean-field approximation. This suggests that at $T=0$, where fluctuations are expected to be suppressed, the results are likely to be a good approximation to the $N_{c}=3$ case. To be sure that this fact is also valid at nonzero temperature, one can remember, e.g., Refs. Radzhabov , where it was shown that mean-field approximation is a rather good approximation both at zero and finite temperature. So not only in the limit $N_{c}\to\infty$ but also at finite $N_{c}$ thermal fluctuations are not that large and probably cannot destroy the results, obtained in this approximation. Hence, in the following we may compare our NJL4 results with lattice $N_{c}=3$ QCD results at $T>0$. Despite all this arguments one should be very cautious comparing the results obtained in different approaches with different setups and approximations, and it should be mentioned that the agreement can be qualitative and not very precise.

II.4 Technical details

Technically, to define the ground state of the system one should find the coordinates $(M_{0},\Delta_{0})$ of the global minimum point (GMP) of the TDP (31). Since the NJL model is a non-renormalizable theory we have to use fitting parameters for the quantitative investigation of the system. We use the following, widely used parameters:

[TABLE]

In this case at $\mu=\nu=\nu_{5}=0$ one gets for constituent quark mass the value $M=309\,{\rm MeV}$. Moreover, we suppose that quark chemical potentials are varied in the region $\mu<\Lambda$, $\nu<\Lambda$ and $\nu_{5}<\Lambda$. At higher values of $\mu$, $\nu$ and $\nu_{I5}$ the NJL model (5) no longer describes a phase structure of real quark matter. The reason is rather obvious and in addition in this case it is necessary to take into account the condensation of $\rho$ mesons, color superconductivity phenomenon, etc. Actually, even though it is shown in Cao:2015xja that if one includes $\mu_{5}$ into consideration the transition to the color superconducting phase shifts to higher values of $\mu$, it does not in any way forbid this phenomenon and the color superconducting phase can appear in the region under consideration. So we should admit that it is interesting to include the possibility of color superconductivity but for simplicity here we will neglect it.

As our main goal of the present paper is to prove the possibility of the charged PC phenomenon in hot dense quark matter with chiral imbalance, i.e. in the framework of the NJL model (5), the consideration of the physical quantity $n_{q}$, called quark number density, is now in order. This quantity is a very important characteristic of the ground state, especially in dynamical phenomena such as superfluidity. It is related to the baryon number density as $n_{q}=3n_{B}$ because $\mu=\mu_{B}/3$. In the general case this quantity is defined by the relation

[TABLE]

where $M_{0}$ and $\Delta_{0}$ are coordinates of the GMP of a thermodynamic potential. In addition, one can find also the density $n_{I}$ of isospin, $n_{I}=-\partial\Omega(M_{0},\Delta_{0})/\partial\mu_{I}$, as well as the chiral isospin density $n_{I5}$, $n_{I5}=-\partial\Omega(M_{0},\Delta_{0})/\partial\mu_{I5}$.

We distinguish the following phases that could be realized in the chirally asymmetric system under different external circumstances (the quantities $M_{0}$ and $\Delta_{0}$ below are the coordinates of the GMP of the TDP (31) in the corresponding phase):

• •

$M_{0}=0;\Delta_{0}=0$ – symmetrical phase. It could be realized only in the chiral limit, $m_{0}=0$. Usually, in this phase $n_{q}\neq 0$ at $\mu\neq 0$.

• •

$M_{0}\neq 0;\Delta_{0}=0;n_{q}=0$ – chiral symmetry breaking phase (we use for it the notation CSB). Since quark number (baryon) density is zero in this phase, sometimes it is called the ordinary baryonic vacuum.

• •

$M_{0}\neq 0;\Delta_{0}=0;n_{q}\neq 0$ – chiral symmetry breaking phase with nonzero quark density (below it is CSBd phase).

In the CSB phase the order parameter $M_{0}$ is usually greater than quark number chemical potential $\mu$. Moreover, $M_{0}$ is of order of the gap in the energy spectrum of quarks. Due to this reasons quarks cannot be created in this phase and $n_{q}=0$. However, with increasing of chemical potentials, it is advantageous for the system to abruptly decrease the parameter $M_{0}$ (see, e.g., the right panel of Fig. 5) and move into a new CSBd phase. In this case, the gap in the energy spectrum of quarks significantly decreases, which makes it possible to create quarks in the ground state. As a result, the quark number density $n_{q}$ is nonzero in the CSBd phase.

• •

$M_{0}\neq 0;\Delta_{0}\neq 0;n_{q}=0$ – charged pion condensation phase with zero quark density (below in all phase diagrams we use for it the notation PC) ($M_{0}=0$ in the chiral limit). In the charged PC phase $U_{I_{3}}(1)$ symmetry is spontaneously broken down. Since in this phase $n_{q}=0$, sometimes it is called the charged pion gas phase.

• •

$M_{0}\neq 0;\Delta_{0}\neq 0;n_{q}\neq 0$ – charged pion condensation phase with nonzero quark density (PCd). In the PCd phase $U_{I_{3}}(1)$ symmetry is also spontaneously broken down. 444The transition between PC and PCd phases is also a first-order phase transition, as in this case the order parameter $\Delta_{0}$ decreases by a jump (see, e.g., the left panel of Fig. 5), and the possibility for the creation of quarks appears. Therefore, in the PCd phase the quark number density $n_{q}$ is nonzero. Moreover, in both phases the isospin density $n_{I}$ is nonzero.

• •

We use the notation ApprSYM for the approximate symmetrical phase. In the literature this phase is usually called Wigner-Weyl phase Klevansky:1992qe ; Zong . It also corresponds to a GMP of the TDP (31), in which $M_{0}\neq 0$ and $\Delta_{0}=0$. But in contrast to the CSB and CSBd phases, dynamical quark mass $M_{0}$ in the ApprSYM phase drops rapidly and continuously to the current quark mass $m_{0}$ with increasing temperature or chemical potentials. As it follows from Eqs. (7) and (16), under such conditions the chiral condensate $\langle\bar{q}q\rangle$ is almost zero, and the chiral symmetry is approximately restored in the model. Moreover, at $m_{0}\to 0$ this phase turns into an exactly symmetrical phase with $M_{0}=0$. These are the reasons why we use the name ApprSYM in all phase portraits below.

Note that at zero temperature $M_{0}$ changes its value by a jump when there is a phase transition from different CSB or charged PC phases to the ApprSYM phase (see, e.g., in Figs 4, 5). However, at nonzero temperature there is usually a chiral crossover transition between CSB and ApprSYM phases (see in Fig. 8).

Below we present different phase portraits of the model as well as its properties in terms of this notations.

III Phase structure of the model

III.1 Exact duality in the chiral limit $(m_{0}=0)$ at zero temperature $(T=0)$

Let us first consider some equilibrium properties of the model starting from the TDP (17) or (31), i.e. at zero temperature, and in the chiral limit ($m_{0}=0$). Although this case has been investigated in details in the article Khunjua:2017mkc , it is useful to recall the main features of the model phase structure obtained in the leading order of the large-$N_{c}$ expansion.

It was already noted above that in the chiral limit the TDP (17) is invariant under the so-called duality transformation ${\cal D}$, where ${\cal D}:~{}M\leftrightarrow\Delta,~{}~{}\nu\leftrightarrow\nu_{5}$, which could be strictly seen from Eq. (24). It means that if at some fixed values $\mu,\nu=A,\nu_{5}=B$ of the chemical potentials the TDP has a GMP of the form $(M=M_{0},\Delta=\Delta_{0})$, then at the transposed values of the isospin chemical potentials, i.e. at $\nu=B,\nu_{5}=A$, but at the unaltered value of $\mu$, the GMP of the TDP (17) lies already at the point $(M=\Delta_{0},\Delta=M_{0})$. As a result, we see that if at $\mu,\nu=A,\nu_{5}=B$, e.g., the CSB phase is realized with $M=M_{0}\neq 0,\Delta=0$, then at the permuted (we say dually conjugated) values $\mu,\nu=B,\nu_{5}=A$ of chemical potentials the charged PC phase should be realized with $M=0,\Delta=M_{0}$, and vice versa. Hence, in the $(\nu,\nu_{5})$-phase portrait all charged PC phases should be arranged mirror symmetrically to all CSB phases with respect to the line $\nu=\nu_{5}$. However, the symmetrical phase turns into itself under the duality transformation, and on the $(\nu,\nu_{5})$-plane the line $\nu=\nu_{5}$ is the axis of symmetry of this phase. Just these facts are well illustrated by the $(\nu,\nu_{5})$-phase diagrams of Fig. 1. There one can see three $(\nu,\nu_{5})$-phase portraits of the model: the left panel corresponds to $\mu=0\,{\rm MeV}$, at the central panel $\mu=150\,{\rm MeV}$ and at the right one $\mu=200\,{\rm MeV}$. Moreover, It is clear from the phase diagrams of Fig. 1 that in dense quark matter, i.e. at $\mu>0$, $\nu_{5}$-chemical potential does promote the charged PC phase with nonzero quark density (there it is PCd phase).

So, in the presence of duality the knowledge of a phase of the model (5) at some fixed values of external free model parameters $\mu,\nu,\nu_{5}$ (and at $m_{0}=0$) is sufficient to understand what a phase (we call it a dually conjugated) is realized at rearranged values of isospin chemical potentials, $\nu\leftrightarrow\nu_{5}$, at fixed $\mu$. Furthermore, different physical parameters such as condensates, densities, etc, which characterize both the initial phase and the dually conjugated one, are connected by the main duality transformation ${\cal D}$. For example, the chiral condensate of the initial CSB phase at some fixed $\mu,\nu,\nu_{5}$ is equal to the charged-pion condensate of the dually conjugated charged PC phase. The quark number density $n_{q}(\nu,\nu_{5})$ (36) of the initial CSB phase is equal to the quark number density in the dually conjugated charged PC phase, etc.

Perhaps, the duality between CSB and charged PC phases is valid in the framework of the NJL4 model under consideration only in the leading large-$N_{c}$ order (and at $m_{0}=0$). However, we think that some signs of this duality remain at the physical point of the full theory and can be observed, e.g., using lattice calculations. What gives us duality? If exact or approximate dual symmetry between different phenomena exists in the model, then, knowing the phase structure or other thermodynamic characteristics of the model in a certain region of chemical potentials, one can predict its properties in the dual-conjugated domain. For example, due to the duality between CSB and charge PC phenomena, there was no need to investigate numerically the TDP (31) at each point of the $(\nu,\nu_{5})$-plane in order to find the phase diagrams of Fig. 1 (or the similar diagrams at other values of $\mu$). Instead, it would be sufficient to obtain a phase portrait in a more narrow region, e.g., at $\nu\geq\nu_{5}\geq 0$. In this case it is composed of PC, PCd and symmetrical phases (see in Fig. 1). Then one should transform each phase of it, using the mapping $\nu\longleftrightarrow\nu_{5}$, into a dually conjugated phase, which is already located in the region $\nu_{5}\geq\nu\geq 0$. At the same time we should change the name of the phase according to the rule: PC$\to$CSB, PC${}_{d}\to$CSBd and the name of the symmetric phase under the dual transformation does not change. Thus, the duality property of the model can help to save not only the time of numerical calculations but also immediately imagine the properties of the model in previously unexplored regions of the values of chemical potentials.

There is an even more interesting use of duality. So, if we know, for example, the $(\nu,\mu)$-phase portrait of the model at fixed $\nu_{5}=A$, there is no need to perform detailed calculations in order to obtain its $(\nu_{5},\mu)$-phase portrait at fixed $\nu=A$. To do this, it is enough to rename the $\nu$ axis of the initial phase diagram to the $\nu_{5}$ axis and change the name of the phases according to the rule: PC$\to$CSB, PC${}_{d}\to$CSBd (symmetrical phase remains intact). We call this technical procedure as the dual conjugation of a phase diagram. Hence, the $(\nu_{5},\mu)$ and $(\nu,\mu)$-phase portraits are mutually conjugate to each other. However, any $(\nu,\nu_{5})$-phase portrait (such as in Fig. 1) is self-dual, i.e. it is transformed into itself by the dual conjugation.

Finally note that there is another kind of duality, the duality between chiral symmetry breaking and superconductivity phenomena, which is realized in some (1+1)- and (2+1)-dimensional four-fermion theories thies ; ekkz2 . But in these models the duality is a consequence of Pauli–Gürsey symmetry of initial Lagrangians.

III.2 Approximate duality in the case of $m_{0}\neq 0$ and $T=0$

In the present section we study the influence of a nonzero value (35) of the bare quark mass $m_{0}$ on the charged PC phase. Moreover, since at $m_{0}\neq 0$ the TDP (17) is no more invariant with respect to the dual symmetry ${\cal D}$ (25), which is exact only in the chiral limit, we will examine the question whether there are some formal signs indicating that the dual symmetry ${\cal D}$ is at least an approximate symmetry of the NJL model at $m_{0}\neq 0$. Among these signs are the following features of the NJL model at the physical point, when $m_{0}=5.5$ MeV,

1. (i)

At some reliable values of the chemical potentials each $(\nu,\nu_{5})$-phase portrait of the model (at some fixed $\mu$) is approximately self-dual, i.e. approximately all charged PC phases of it are arranged mirror symmetrically to all CSB phases with respect to the line $\nu=\nu_{5}$. 2. (ii)

Each $(\nu,\nu_{5})$-phase diagram has a phase (it is the ApprSYM phase), which is approximately symmetric under the transformation $\nu\leftrightarrow\nu_{5}$, i.e. it is arranged symmetrically with respect to the line $\nu=\nu_{5}$. 3. (iii)

Under the dual transformation, when $\nu\leftrightarrow\nu_{5}$, the order parameter $M_{0}$ of CSB or CSBd phase is approximately equal to the order parameter $\Delta_{0}$ of the dually conjugated charged PC or PCd phase. 4. (iv)

The quark number density $n_{q}$ in any phase, corresponding to the chemical potential point $(\mu,\nu=A,\nu_{5}=B)$, is approximately equal to quark number density $n_{q}$ of its dually conjugated phase that lies at the point $(\mu,\nu=B,\nu_{5}=A)$. 5. (v)

Each $(\nu_{5},\mu)$-phase portrait (at some fixed $\nu=A$) of the model is approximately the dual ${\cal D}$ mapping of a corresponding $(\nu,\mu)$-phase portrait (at some fixed $\nu_{5}=A$) and vice versa.

If these properties are inherent in the model or theory, then we say that in the model (theory) there is an approximate duality between its chiral properties and charged pion condensation phenomena.

Bearing this in mind, let us look at the $(\nu,\nu_{5})$-phase portraits of Fig. 2, which are depicted for the same values of the quark number chemical potential $\mu$ as in Fig. 1. First of all note that at $\nu_{5}=0$ in all diagrams of Fig. 2 there is a threshold $\nu_{c}=m_{\pi}/2\approx 70$ MeV of a second order phase transition to the PC phase, which is also predicted by all known investigations Son:2000 (including lattice calculations Gupta ). Moreover, it is easily seen from these diagrams that $\nu_{5}$ promotes the charged PC phase in dense quark matter (it is the phase PCd in Fig. 2 and 3) even in the case of $m_{0}\neq 0$. The NJL model is believed to work well at rather high baryon densities where quark matter can be realized and for low values of chemical potential $\mu<300$ MeV NJL model is likely to give not very trustworthy results and hadron effective model is needed. So the prediction of the generation of PCd phase at $\mu<300$ MeV is questionable and the use of hadron effective model, taking into account the presence of baryons, is needed in this region but one can consider the prediction of PCd phase generation with rather large baryon density (at $\mu>300$ MeV) to be reliable (see Fig. 3).

Concerning the above-listed duality signs (i)-(v), we see that in the region $\omega=\{(\mu,\nu,\nu_{5}):\mu<\Gamma(m_{0}),\nu<\Gamma(m_{0}),\nu_{5}<\Gamma(m_{0})\}$, where $\Gamma(m_{0})$ is of the order of the pion mass $m_{\pi}$, there is no sense to say about duality (even approximate), because the point (i) of this list is not fulfilled. However, as it follows from Figs. 2 and 3, outside the region $\omega$ and for all values of $\mu$, $\nu$ and $\nu_{5}$ restricted by the conditions $\mu<\Lambda$, $\nu<\Lambda$ and $\nu_{5}<\Lambda$ (the duality is even better symmetry in the region of larger values of chemical potentials but the results of NJL model in this region are not trustworthy) we see that the items (i) and (ii) are satisfied.

To have a more precise picture, let us take a look at the Figs. 4 and 5, where the Gaps $M_{0},\Delta_{0}$ and baryon density $n_{B}$ vs. $\nu_{5}$ and $\nu$ for are depicted. It follows from these pictures that if we go from the phase, corresponding, e.g., to a chemical potential set $(\mu=200,\nu=350,\nu_{5}=A)$ MeV, to the dually ${\cal D}$ conjugated phase with $(\mu=200,\nu=A,\nu_{5}=350)$ MeV (or vice versa), then pion condensate $\Delta_{0}$ in the charged PC phase is approximately the same as dynamical quark mass $M_{0}$ in the dually conjugated CSB phase (compare the left and right panels of Fig. 5) and baryon density $n_{B}$ is not changed (approximately). In the dually conjugated points of the ApprSYM phase both $n_{B}$ and dynamic quark mass $M$ are not changed, approximately. The same conclusions one can obtain from Fig. 4 for $\mu=260$ MeV when two phases, CSB and PC, are present. Hence, the items (iii) and (iv) of the list of duality signs are also satisfied.

Finally, comparing, e.g., the $(\nu,\mu)$-phase diagram at fixed $\nu_{5}=200$ MeV and the $(\nu_{5},\mu)$-phase diagram at fixed $\nu=200$ MeV (see in Fig. 6), we see that qualitatively they are dually ${\cal D}$ conjugated to each other at a rather low values of $\mu\lesssim 200$ MeV, i.e in this region of each diagram of Fig. 6 one can perform the following axis and phase renaming, $\nu\leftrightarrow\nu_{5}$, CSB$\leftrightarrow$PC and CSB${}_{d}\leftrightarrow$PCd (the ApprSYM phase does not change its name by the duality transformation), in order to obtain (approximately) the corresponding region of another diagram of Fig. 6. This conclusion agrees with phase portraits of Fig. 2 for moving along the lines $\nu=200$ MeV (or $\nu_{5}=200$ MeV) of these diagrams we intersect just the phases shown in Fig. 6 at low $\mu$. In addition, it is easy to see that there is a duality ${\cal D}$ between diagrams of Fig. 6 in the regions, where $\nu\gtrsim 200$ MeV (left panel) and $\nu_{5}\gtrsim 200$ MeV (right panel). So the item (v) of the list of duality signs is also satisfied.

In conclusion of this section, we can say that the duality between the phenomena of CSB and a charged PC, inherent for this model in the chiral limit at $N_{c}\to\infty$, is approximately fulfilled even at $m_{0}\neq 0$, but only for the points $(\mu,\nu,\nu_{5})$ of the chemical potential space from the region, in which $\mu,\nu,\nu_{5}<500$ MeV and at the same time $(\mu,\nu,\nu_{5})\not\in\omega$, where $\omega=\{(\mu,\nu,\nu_{5}):\mu,\nu,\nu_{5}<\Gamma(m_{0})\}$ (here $\Gamma(m_{0})\sim m_{\pi}$).

III.3 Phase portrait at the physical point ($m_{0}=5.5$ MeV) and nonzero temperature $(T\neq 0)$

Though, the effect of nonzero temperatures is quite predictable (indeed, one can expect that the temperature just restores all the broken symmetries of the model), we investigate nonzero temperature case because it is important in a number of applications. We know that compact stars are cold and one can consider their temperatures as zero, but probably there could be scenarios in which the temperatures could be important even in the context of compact stars. So it is instructive to know how robust the PCd-phase under the influence of temperature and chiral imbalance.

To clarify this issue, we calculated two $(\nu,T)$-phase diagrams of the model at $\mu=0$ (in order to compare our results with lattice investigations) and at different values of $\nu_{5}$. In Fig. 7 (left panel) one can see this diagram at $\nu_{5}=0$, whereas in the right panel it is at $\nu_{5}=200$ MeV. Note that the phase portrait at $\mu=0$ and $\nu_{5}=0$ is in accordance with the same phase portrait obtained within first principle lattice calculations Gupta . Also, as one could expect, it is clear from Fig. 7 that temperature restores broken $U_{I_{3}}(1)$-symmetry at some rather high critical values $T_{c}^{PC}$, where charged PC phase is disappeared. 555In the points of the boundary between ApprSYM and charged PC phases in Figs 7, 8 there is a second-order phase transition, whereas the dashed line in each of these figures represents the so-called pseudo-critical temperature, which characterizes the so-called cross-over region between CSB and ApprSYM phases. (Of course, at fixed values of $\mu$ and $\nu_{5}$ the critical temperature $T_{c}^{PC}$ depends strongly on the isospin chemical potential $\mu_{I}\equiv 2\nu$.) Moreover, as it follows from Fig. 7, $T^{PC}_{c}$ vs. $\nu_{5}$ (at $\mu=0$) drops from values $T^{PC}_{c}\approx 200$ MeV at $\nu_{5}=0$ to values $T^{PC}_{c}\approx 100$ MeV at $\nu_{5}=200$ MeV (compare left and right panels of Fig. 7), i.e. when $\nu_{5}$ increases the region of the PC phase is shrinked in the phase portrait of the model (this fact is in accordance with the phase diagrams of Fig. 2), but nevertheless the charged PC is a quite robust effect vs temperature at $\mu=0$.

Finally, we would like to note that in addition to the list of signs (i)-(v) (see in the previous section 3) indicating on the presence in the NJL model (5) at $m_{0}\neq 0$ of an approximate dual symmetry between CSB and charged pion condensation at $T=0$, the similar approximate dual correspondence between phase diagrams exists also at nonzero temperature. For example, two diagrams of Fig. 8, the $(\nu,T)$-phase portrait at fixed $\mu=200$ MeV and $\nu_{5}=200$ MeV (left panel) and the $(\nu_{5},T)$-phase portrait at fixed $\mu=200$ MeV and $\nu=200$ MeV (right panel), can be considered as a dually conjugated to each other. Indeed, applying to each of these diagrams the dual mapping, i.e. the following replacements $\nu\leftrightarrow\nu_{5}$, CSB$\leftrightarrow$PC and CSB${}_{d}\leftrightarrow$PCd, it is possible to obtain approximately another diagram. So dual mapping of a well-known phase portraits can be used in order to predict (approximately) a phase structure of the model at $m_{0}\neq 0$ in the dually conjugated region, i.e. at $\nu\leftrightarrow\nu_{5}$.

Last but not least conclusion from Fig. 8 (in addition to the analysis of Fig. 7) is that the charged PC is also a rather temperature stable effect in dense quark matter (at $\mu>0$).

III.4 Pseudo-critical temperature $T_{c}(\nu_{5})$ in the NJL4 model: comparison with lattice QCD and other approaches

A rather significant part of the previous section was devoted to the consideration of the critical temperature $T^{PC}_{c}$ (at $\mu=0$ and different fixed values of $\nu$) of the second-order phase transition from the charged PC to the ApprSYM phase, as well as its dependence on $\nu_{5}$. In addition to this, in the present section we will study in the framework of the NJL4 model the behavior (at $\mu=0$ and $\nu=0$) of the pseudo-critical temperature $T_{c}(\nu_{5})$, which characterizes the so-called chiral cross-over region of the phase diagram (see, e.g., Figs. 7, 8 where this region is arranged around the dashed lines). It separates the low-temperature phase with CSB and a (partially) chirally restored ApprSYM phase, which corresponds at high temperatures to quark-gluon plasma. In the cross-over region different physical parameters, such as the dynamical quark mass $M_{0}$ etc, of the CSB phase smoothly (without jumps), but rather sharply go over to the corresponding parameters of the ApprSYM phase. Therefore, in this region, there occurs not a true phase transition with corresponding critical temperature, etc., but rather a pseudo-phase transition (cross-over) characterized by a pseudo-critical temperature $T_{c}\equiv T_{c}(\mu,\nu,\nu_{5})$, etc. (In Figs. 7, 8 the pseudo-critical temperature $T_{c}$ is represented by dashed lines.) Here we study the behavior of the pseudo-critical temperature only as a function of $\nu_{5}$ and at fixed $\mu=0$, $\nu=0$. That is, we investigate the quantity denoted by $T_{c}(\nu_{5})$,

[TABLE]

In particular, it is clear from Fig. 7 that $T_{c}$ at $\nu_{5}=0$ MeV (left panel of Fig. 7) is slightly smaller than $T_{c}$ at $\nu_{5}=200$ MeV, which can be found at the right panel of Fig. 7. The plot of the function $T_{c}(\nu_{5})$ vs $\nu_{5}$ is presented in Fig. 9. However, before comparing these our results on the pseudo-critical temperature with the predictions, obtained in the framework of other effective models and lattice QCD calculations, it is necessary to make a few remarks.

Strictly speaking, so far nobody has investigated the function $T_{c}(\nu_{5})$ (37) both in the NJL model and other approaches. The matter is that in the most general case, the chiral asymmetry of dense quark matter is described by two chemical potentials, chiral $\mu_{5}$ and chiral isospin $\mu_{I5}\equiv 2\nu_{5}$ chemical potential. 666In general, chiral imbalance of dense quark matter is characterized by two densities, chiral isospin $\hat{n}_{I5}=\frac{1}{2}\left(\hat{n}_{u5}-\hat{n}_{d5}\right)$ and chiral density $\hat{n}_{5}=\hat{n}_{u5}+\hat{n}_{d5}$ (see Introduction for notations). Alternatively, it can be described by corresponding chemical potentials $\mu_{I5}$ and $\mu_{5}$, which are the quantities thermodynamically conjugated to $\hat{n}_{I5}$ and $\hat{n}_{5}$, respectively. The first, $\mu_{5}$, is usually used when isotopic asymmetry of quark matter is absent, i.e. in the case $\mu_{I}=0$ Ruggieri . The second, $\mu_{I5}$, might be taken into account when, in addition to chiral, there is also isotopic asymmetry of matter, in which charged PC phenomenon can be observed, etc. Khunjua:2017mkc ; Khunjua:2018sro . And up to now the behavior of a pseudo-critical temperature of the cross-over region as a function of only the chiral chemical potential $\mu_{5}$ was investigated in different approaches at fixed $\mu=0$, $\mu_{I}=0$ and $\nu_{5}=0$ Fukushima2 ; Ruggieri2 ; Frasca . That is the possibility of the existence of a quark system with nonzero chiral isospin imbalance was ignored in these works. (In this particular case we use for a pseudo-critical temperature of the NJL4 model the notation $T_{c}\equiv\widetilde{T}_{c}(\mu_{5})$. Do not confuse with the expression (37), which in fact corresponds to a pseudo-critical temperature, obtained in another limiting case of an external parameter set of the NJL4 model, $\mu=0$, $\nu=0$, $\mu_{5}=0$ and for arbitrary values of $\nu_{5}$.)

So in order to compare our results on the pseudo-critical temperature with other approaches to this quantity, we need formally to find in the model under consideration the behavior of a pseudo-critical temperature vs $\mu_{5}$, i.e. the quantity $T_{c}\equiv\widetilde{T}_{c}(\mu_{5})$.

Recall that for simplicity, in the present paper we study the phase structure of the NJL4 model (5) only in the case $\mu_{5}=0$ with other nonzero chemical potentials, however in the recent paper Khunjua:2018sro the phase structure of this model was investigated in the chiral limit in a more general case with all four nonzero chemical potentials $\mu,\mu_{5},\mu_{I},\mu_{I5}$. In particular, it was established in Khunjua:2018sro that in addition to the dual symmetry (25) the TDP of the NJL4 model is invariant with respect to a transformation ${\cal D}_{M}$,

[TABLE]

In Ref. Khunjua:2018sro the symmetry (38) is called constrained duality (due to the relation $\Delta=0$). Note that, in contrast to the main duality relation ${\cal D}$ (25), which is only an approximate symmetry between CSB and charged PC phenomena at the physical point (see in the sections above), the constrained duality ${\cal D}_{M}$ is an exact symmetry of phase portraits of the model even at $m_{0}\neq 0$. It means that in all phase diagrams of the present paper, obtained for $\mu_{5}=0$ and $\nu_{5}\geq 0$, one can treat $\nu_{5}$ outside of the charged PC phase as a chiral chemical potential $\mu_{5}$. Moreover, just due to this additional dual symmetry the following relation between pseudo-critical temperatures of the NJL4 model is valid,

[TABLE]

where $T_{c}(\nu_{5})$ and $\widetilde{T}_{c}(\mu_{5})$ are the pseudo-critical temperatures of this model in the cases $\mu=\nu=\mu_{5}=0$ and $\mu=\nu=\nu_{5}=0$, respectively. The duality can be used only in the case of $\Delta=0$ and this assumption was shown to be true in these cases in Khunjua:2018sro . So it can be used in the considered cases, $\mu=\nu=\mu_{5}=0$ and $\mu=\nu=\nu_{5}=0$.

Just the relation (39) gives us the possibility to compare our results with previous predictions for the (pseudo-)critical temperature $\widetilde{T}_{c}(\mu_{5})$ obtained in the framework of different effective models Ruggieri ; Fukushima2 ; Ruggieri2 ; Frasca , using the Dyson-Schwinger equations Xu ; Wang and lattice calculations Braguta:2015zta ; Braguta:2015owi ; Yamamoto . It should be noted that all these studies do not provide a well-defined consistent prediction for the behavior of $\widetilde{T}_{c}(\mu_{5})$ vs $\mu_{5}$, and rather contradict each other. For example, the works Ruggieri ; Fukushima2 predict decreasing of $\widetilde{T}_{c}(\mu_{5})$ with growing $\mu_{5}$, whereas in the works Xu ; Frasca there is an opposite picture. The situation has been partially clarified in the works Yu ; Ruggieri2 ; Farias ; Cui , where it has been shown that in the framework of the effective models, regularisation scheme could play a crucial role in the behaviour of the pseudo-critical temperature $\widetilde{T}_{c}(\mu_{5})$. There are several approaches to regularization of effective models. For example, one can regularize the whole TDP $\Omega_{T}(M,\Delta)$ (34), including both the vacuum $\Omega(M,\Delta)$ and thermal terms of Eq. (34), or regularize only the vacuum term $\Omega(M,\Delta)$, etc. It was shown in Refs. Yu ; Farias ; Cui that the behaviour of the $\widetilde{T}_{c}(\mu_{5})$ in the NJL model depends strongly on the scheme that used. And in the present paper, as it follows from Eq. (34), we regularize the whole TDP $\Omega_{T}(M,\Delta)$ using in Eq. (34) the so-called hard-cutoff regularization scheme when the integration region of the thermal part of the TDP is restricted by the cutoff parameter $\Lambda$.

Taking into account the relation (39), it is easy to see from the plot of Fig. 9, that in this regularization scheme the pseudo-critical temperature $\widetilde{T}_{c}(\mu_{5})$ of the NJL model increases for $\mu_{5}<\mu^{\ast}_{5}\lessapprox 400$ MeV. Above this value it drops down, but at $\mu_{5}>\mu^{\ast}_{5}$ the NJL4 model, in our opinion, does not provide very reliable predictions, because $\mu_{5}$ is near the cutoff $\Lambda$. The similar quantity has been investigated in the NJL4 model with the same regularization scheme but in the chiral limit Yu . In this case $\widetilde{T}_{c}(\mu_{5})$ is no more a pseudo-critical but rather a critical temperature of a 2nd order phase transition from CSB to symmetrical phase. In contrast to Ref. Yu , we study the NJL4 model at the physical point ($m_{0}\neq 0$). However, the behavior of $\widetilde{T}_{c}(\mu_{5})$ vs $\mu_{5}$ in both cases is qualitatively the same.

The selection of such a regularization scheme is supported and justified by several things. Namely, using the first principle lattice calculations, it was shown that the $\widetilde{T}_{c}(\mu_{5})$ increases with $\mu_{5}$. Then, there is a qualitative description of the mechanism leading to an increase in the pseudo-critical temperature $\widetilde{T}_{c}(\mu_{5})$. It is based on the Fermi-sphere treatment Braguta and backed up with the results achieved in the framework of different non-perturbative methods Xu ; Wang . And finally, and maybe most importantly in the context of the present work, we are guided by predictions for $\widetilde{T}_{c}(\mu_{5})$, which follow from the duality symmetries (25) and (38) of the model. Indeed, it is well established that in the $(\nu,T)$-phase portrait the critical temperature $T_{c}^{PC}$ at which isospin symmetry is restored increases with the chemical potential $\nu$ (see, e.g., the left panel in Fig. 7). Applying to this diagram an (approximate) duality transformation ${\cal D}$ (25), we obtain a $(\nu_{5},T)$-phase diagram corresponding to $\nu=0$, $\mu=0$ and $\mu_{5}=0$ with horizontal $\nu_{5}$ axis as well as with the CSB phase at $\nu_{5}\gtrapprox 0.1$ GeV. On the boundary between the ApprSYM and CSB phases there will most likely be a cross-over region with a pseudo-critical temperature $T_{c}(\nu_{5})$ (37) that, due to the approximate dual symmetry ${\cal D}$, should increase vs $\nu_{5}$, as it does $T_{c}^{PC}$ vs $\nu$ in Fig. 7. And finally, applying to this phase diagram the constrained duality transformation ${\cal D}_{M}$ (38), we obtain a $(\mu_{5},T)$-phase diagram corresponding to $\nu=0$, $\mu=0$ and $\nu_{5}=0$ from which it is clear that $\widetilde{T}_{c}(\mu_{5})$ also rises vs $\mu_{5}$. It is this qualitative analysis based on the duality properties of the NJL4 model that is confirmed by Eq. (39) along with the plot of Fig. 9. The duality is only approximate but we also saw in the previous sections that it is a good approximation for values of chemical potential larger than approximately pion mass.

Note that in lattice approach to QCD the simplest $(\nu,T)$- and $(\mu_{5},T)$-phase diagrams at $\mu=0$ are well investigated. Moreover, they are in accordance with the similar phase diagrams, obtained in the framework of the NJL model (although it should be mentioned that the agreement can be not very precise because in NJL model the results were obtained in the large-$N_{c}$ limit and in lattice QCD with $N_{c}=3$). But in the last approach, as it follows from above consideration, these phase portraits are (approximately) dually conjugated to each other. Consequently, the same connection can exist between these phase diagrams in real QCD. So there is a solid foundation, the lattice QCD simulations, which allows us to hope that duality between CSB and charged PC phenomena is one of the properties of real dense quark matter.

IV Summary and Discussion

In this paper the influence of isotopic and chiral imbalance on phase structure of hot/cold dense quark matter has been investigated at the physical point (i.e. at nonzero current quark mass $m_{0}$) in the framework of the (3+1)-dimensional NJL model with two quark flavors in the large-$N_{c}$ limit ($N_{c}$ is the number of colors). Dense matter means that our consideration has been performed at nonzero baryon $\mu_{B}$ chemical potential. Isotopic and chiral imbalance in the system were accounted for by introducing isospin $\mu_{I}$ and chiral isospin $\mu_{I5}$ chemical potentials (see Lagrangian (5)). Of course one knows that current quark masses of $u$ and $d$ quarks (the ones that we considered in the paper) are rather small and, in general, the chiral limit is a very good approximation. But sometimes although small but nonzero masses can change some aspects of the phase diagram. For example, charged pion condensation (PC) phase in the chiral limit and at $T=0$ starts from infinitesimally small values of isospin chemical potential $\mu_{I}$, but when one takes into account quark masses, then it shifts the charged PC to the values of $\mu_{I}$ larger than pion mass $m_{\pi}\approx 140$ MeV (compare diagrams of Figs 1 and 2). The phase structure of cold dense quark matter in the chiral limit has been obtained in Khunjua:2017mkc ; Khunjua:2018sro , where it has been shown that chiral isospin chemical potential $\mu_{I5}$ generates charged PC in dense quark matter and there is a duality correspondence between CSB and charged PC phenomena in the leading order of the large-$N_{c}$ approximation. The goal of the present paper is the extension of this consideration to a more physical case of the NJL4 model with nonzero current quark masses. This allows us to draw more accurate phase diagram and perform comparison with lattice QCD. Moreover, we take into account finite temperatures, which give us a chance to consider the results in the context of heavy ion collision experiments in which temperatures are always rather large. Even in the context of neutron stars it can be interesting to consider the case of finite temperature. It has been found that the duality between CSB and charged PC phenomena observed in Khunjua:2017mkc ; Khunjua:2018sro in the chiral limit (where it was exact) is valid with good accuracy even in the physical point. It has been also shown that temperature does not spoil the duality as well.

We have studied the full ($\mu_{B}$, $\mu_{I}$, $\mu_{I5}$, $T$)-phase diagram of quark matter in terms of the NJL4 model with $m_{0}\neq 0$. This general consideration is not feasible in the lattice QCD simulations, mainly due to the famous sign problem which does not allow for the consideration of finite baryon densities (nonzero baryon chemical potential $\mu_{B}$). But contrary to the case of non-zero baryon chemical potential, simulations with non-vanishing isospin $\mu_{I}$ and chiral $\mu_{5}$ chemical potentials are not hampered by a sign problem and some particular cases have been considered on the lattice. For example, the ($\mu_{I}$, $T$)-phase diagram at zero values of $\mu_{B},\mu_{I5},\mu_{5}$ chemical potentials is well established as in lattice QCD as well as in different effective models and a rather good agreement can be observed between this different approaches. And there are lattice QCD simulations of the quark matter with only nonzero chiral chemical potential $\mu_{5}$ in terms of as $SU_{c}(2)$ QCD (two-colour QCD) Braguta:2015zta as well as real $SU_{c}(3)$ QCD (three-colour QCD) Braguta:2015owi , where the catalysis of chiral symmetry breaking by chiral chemical potential has been established. Namely, it has been shown that chiral condensate and (pseudo)critical temperature (the temperature at which the chiral condensate drops) grows with increase of chiral $\mu_{5}$. In this paper, as well as in Ref. Khunjua:2018sro , we have supported these conclusions by effective NJL model considerations. In paricular, the plot of the pseudo-critical temperature $\widetilde{T}_{c}(\mu_{5})$ vs chiral chemical potential (see Fig. 9 and take into account Eq. (39)) has been drawn and it was shown that this quantity rises with the raise of $\mu_{5}$ and the behaviour is rather similar to the results of the lattice QCD simulations.

So let us gaze at all this from the general picture viewpoint. We have two lattice simulation results, ($\mu_{I}$, $T$)- and ($\mu_{5}$, $T$)-phase diagrams. These phase diagrams have been also obtained in the NJL model and the results are in a good agreement with lattice QCD simulations. But in terms of NJL model we can consider the general case and we know that there is the duality ${\cal D}$ (25) (it is exact in the chiral limit only) between CSB and charged PC phenomena in the leading order of the large-$N_{c}$ approximation. Moreover, there is also the so-called constrained ${\cal D}_{M}$ (38) duality of the NJL4 model phase diagram, which is valid even in the physical point Khunjua:2018sro . So the significant regions (at $\nu_{I},\mu_{5}\gtrapprox m_{\pi}/2$) of the particular ($\mu_{I}$, $T$)- and ($\mu_{5}$, $T$)-phase diagrams should be dually conjugated to each other with respect to a sequential action of two mappings, ${\cal D}$ and ${\cal D}_{M}$ (see the discussion at the end of Sec. 9), but the duality ${\cal D}$ in the case of the physical point is only approximate (see in Secs. 3 and III.3), although it is valid with a good precision. Since the particular ($\mu_{I}$, $T$)- and ($\mu_{5}$, $T$)-phase diagrams in these two approaches agrees, one can conclude that the duality can be observed in the lattice QCD simulations. And this put the notion of the duality on another level of confidence, for it is observed in terms of the toy (1+1)-dimensional NJL model 2dim , effective (3+1)-dimensional NJL model Khunjua:2017mkc ; Khunjua:2018sro , lattice QCD simulations and similar dualities has been observed in the large-$N_{c}$ orbifold equivalences approach. Comparison to the lattice QCD is important not only due to the fact that it is ab initio method for dealing with QCD, but because it does not make use of, for example, large-$N_{c}$ approximation (as in NJL models or in large $N_{c}$ orbifold equivalences approaches).

The question of catalysis of chiral symmetry breaking by chiral chemical potential, i.e. the growth of $T_{c}$ vs $\mu_{5}$, is a rather debated one and there are a number of papers that predicted that (pseudo-)critical temperature decrease with increase of chiral chemical potential $\mu_{5}$ Ruggieri ; Fukushima2 , as well there are a number of papers that support our results Yu ; Xu ; Frasca ; Cui . Different regularization schemes have been applied in Ruggieri2 ; Yu ; Cui ; Farias and it has been stated that if the right one is used there is no catalysis, but lattice QCD results is probably more trustworthy and it disagrees with them. But the catalysis of chiral symmetry breaking by chiral chemical potential $\mu_{5}$ can be established in terms of duality notion, let us elaborate on that. As we have talked about, the ($\mu_{I}$, $T$)-phase diagram is well established one and the duality fails only in the region of small isospin and chiral isospin chemical potentials (smaller than half of the pion mass), but works quite well for the larger values (see, e.g., in Figs 2, 3). But at the ($\mu_{I}$, $T$)-phase diagram in the region of isospin chemical potential larger than half of the pion mass the critical temperature increases when $\mu_{I}$ is raised and the duality here is a good approximation, so the critical temperature at the duality conjugated ($\mu_{5}$, $T$)-phase diagram should increase with rising of $\mu_{5}$ as well (see at the end of Sec. 9).

Let us summarize the core results of our paper.

• •

It has been also demonstrated that the duality correspondence between CSB and charged PC phenomena observed in Refs. Khunjua:2017mkc ; Khunjua:2018sro in the chiral limit (where it was exact) is a very good approximate symmetry of the phase diagram even in the physical point in the framework of the NJL4 model in the leading order of the large-$N_{c}$ approximation (see in Sec. 3). And it stays a very instructive feature of the phase diagram that can be used in different situations.

• •

It has been also shown that temperature does not spoil the duality correspondence between CSB and charged PC phenomena and it stays exact at finite temperature in the chiral limit (see the comment at the end of Sec. II.3) and it is a good approximate symmetry at the physical point (see in Sec. III.3).

• •

We have shown that there is a huge PCd phase region in the phase portrait of the model (1) promoted by $\nu_{5}$ even in the physical point (see, e.g., in Fig. 3). And it has been revealed that PCd phase can exist at rather large temperatures up to even about 100 MeV (see in Fig. 8).

• •

The particular cross sections of the obtained phase portraits are in qualitative accordance with the recent lattice simulations Gupta ; Braguta:2015zta ; Braguta:2015owi . And it has been established that lattice QCD results support the existence of the duality.

• •

The rise of pseudo-critical temperature with increase of chiral chemical potential $\mu_{5}$ has been established in terms of duality notion and the well explored results of lattice QCD and different approaches on phase structure of isotopicaly imbalanced quark matter. It gives additional argument in favour of this behaviour of the pseudo-critical temperature and it is of importance because, although the lattice results confirming this behaviour are conclusive, the pion mass that is used in these simulations is still quite high and well above the physical pion mass and our results are made at the physical point with the right value of the pion mass.

The central result of our paper is the fact that there is an approximate duality that is supported by recent lattice results. Two different lattice simulations that are completely different and are not connected at the first sight are in reality dual to each other. Moreover, the logic can be reversed and we can predict the increase of pseudo-critical temperature with rising of chiral chemical potential, the much debated effect recently, just by the duality notion.

This work is intended to generalize and refine the previously obtained results of Refs. Khunjua:2017mkc ; Khunjua:2018sro to a more physically motivated situation (physical point and finite temperatures). These generalizations require much more computing resources and technically is rather challenging, but it pays off when you can compare the results with lattice QCD and it supports them. Moreover, we hope that our results might shed new light on phase structure of dense quark matter with isotopic and chiral imbalance and hence could be of interest in the context of the heavy ion collision experiments and neutron stars interiors.

V Acknowledgement

The authors are thankful to Igor Shovkovy for the idea of the possibility of generation of chiral isospin imbalance in dense matter due to chiral separation effect.

Appendix A Generation of nonzero chiral isospin charge in dense quark matter

Let us suppose, for simplicity, that dense quark matter consists of two massless $u$ and $d$ quarks, whose chemical potentials, $\mu_{u}=\mu+\nu$ and $\mu_{d}=\mu-\nu$ (see the notations adopted just after the Eq. (6)), are positive. Moreover, we suppose also that quarks do not interact, and there is an external magnetic field $\vec{B}=(0,0,B)$ directed along $z$ axis. In this case in the equilibrium state of quark matter there is a nonzero and nondissipative axial current

[TABLE]

for each quark flavor $f=u,d$. In Eq. (40) $Q_{f}$ is an electric charge of the quark flavor $f$, i.e. $Q_{u}=2/3$, $Q_{d}=-1/3$. In this case it is not difficult to conclude from Eq. (40) that axial currents of $u$ and $d$ quarks are opposite in their directions. Since $\vec{j}_{5f}=\langle\bar{q}_{fR}\vec{\gamma}q_{fR}\rangle-\langle\bar{q}_{fL}\vec{\gamma}q_{fL}\rangle$, where

[TABLE]

we see from Eq. (40) that left- and right-handed quarks of each flavor $f=u,d$ moves in opposite directions of the $z$ axis. As a result, a spatial separation of quark chiralities for each flavor $f$ occurs. It is the so-called chiral separation effect Metlitski . In other words, one can say that in the upper half of the three-dimensional space, i.e. at $z>0$, the density, e.g., $n_{uR}\equiv\langle\bar{q}_{uR}\gamma^{0}q_{uR}\rangle$ of the right-handed $u$ quarks is greater than the density $n_{uL}\equiv\langle\bar{q}_{uL}\gamma^{0}q_{uL}\rangle$ of the left-handed $u$ quarks. Hence, in this case we have at $z>0$ the positive values of the chiral charge density $n_{u5}\equiv n_{uR}-n_{uL}$ for $u$ quarks. (It is evident that at $z<0$ the chiral charge of $u$ quarks is negative.)

On the contrary, since the axial current $\vec{j}_{5d}$ of $d$ quarks differs in its direction from the axial current $\vec{j}_{5u}$ of $u$ quarks, one can see that in this case at $z>0$ (at $z<0$) the density $n_{d5}$ of the chiral charge of $d$ quarks is negative (positive). Consequently, we have at $z>0$ the positive values of the quantity $n_{I5}\equiv n_{u5}-n_{d5}$, which is the ground state expectation value of the density operator for the chiral isospin charge (it is defined in Eq. (4)). Whereas at $z<0$ the chiral isospin charge is negative.

In summary, we can say that in dense quark medium under the influence of a strong magnetic field (as an example we can mention neutral stars), regions with a nonzero chiral isospin charge $n_{I5}$ might appear. Therefore physical processes inside these regions can be described, e.g., in the framework of the Lagrangians of the form (5), containing chiral isospin chemical potential $\mu_{I5}$.

Appendix B Calculation of roots of $P_{\pm}(\eta)$

In this appendix it will be shown how to get roots of the following quartic equation (general quartic equation could be reduced to the one of this form)

[TABLE]

The coefficients $a,b,c$ in Eq. (42) are given by the relations (24). First, we represent the polynomial on the left-hand side of this equation as the product of two quadratic polynomials,

[TABLE]

where

[TABLE]

It follows from these relations that

[TABLE]

Substituting Eq. (44) into Eq. (43), one gets that $r=\sqrt{R}$, where $R$ is one of the solutions of the following cubic equation

[TABLE]

where we used notations $A,B,C$ that are given by

[TABLE]

All three solutions of the cubic equation (45) are

[TABLE]

where

[TABLE]

and $\sqrt[3]{J}$ in Eq. (46) means each of three possible complex valued roots. There is a determinant $D\equiv 4L^{3}-K^{2}>0$ of the equation (45) that can tell us the structure of roots $R_{1,2,3}$. Namely, if $D>0$ then all roots $R_{i}$ are real and different, if $D=0$ all roots are real and at least two are equal. Finally, if $D<0$ then one root is real and two are complex conjugate. So, there is always a real solution of Eq. (45). In numerical simulations it is more handy to work with real solution and it is always possible to choose one. There is a procedure that, depending on values of parameters, chooses a real solution, but it is quite lengthy so we will not present it here.

And when one has found $r$, the roots of Eq. (42) has the following form

[TABLE]

The roots $\eta_{5,6,7,8}$ of the equation $P_{-}\equiv\eta^{4}-2a\eta^{2}-b\eta+c=0$ can be obtained by changing $b\to-b$ in Eq. (42) (or $q\leftrightarrow s$ in Eq. (47) with $r$ unchanged). So, we have

[TABLE]