The dual properties of chiral and isospin asymmetric dense quark matter formed of two-color quarks

In this paper the phase structure of dense baryon matter composed of u and d quarks with two colors has been investigated in the presence of baryon μB , isospin μI and chiral isospin μI5 chemical potentials in the framework of Nambu-Jona-Lasinio model with quark-antiquark and quark-quark interaction channels. In the chiral limit, it has been shown in the mean-field approximation that the duality between phases with spontaneous chiral symmetry breaking and condensation of charged pions, found in the three color case, remains valid in the two color case. In addition, it has been shown that there are two more dualities in the phase diagram in two color case, namely (as in the case with μI5 = 0), at μI5 ≠ 0 the general (μ, μI, μI5)-phase portrait of the model has dual symmetry between the phase with condensation of charged pions and the phase with diquark condensation. This duality stays exact even in the physical point, m0 ≠ 0. And at m0 = 0 the (μ, μI, μI5)- phase portrait becomes even more symmetrical, since dual symmetry between phases with spontaneous chiral symmetry breaking and diquark condensation appears. It is shown that due to the dualities the phase diagram is extremely symmetric and has interlacing structure. One can show that the phase portrait of two-color NJL model can be obtained just by duality properties from the results of investigations of three-color NJL model (it was noticed only after the numerical calculations have been performed). Three-color case shares only one duality of the two color one, and one can only see a facet of this enormously symmetric picture in the case of three colors. Using dualities only, it is possible to show that there are no mixed phases (phases with two non-zero condensates). This prediction of dualities is of great use, because for sure it can be shown by the direct calculations but it would be enormously more complicated and time-consuming numerically.


Introduction
It is well known that quantum chromodynamics (QCD) is the theoretical basis for the investigation of strongly interacting matter. However, to study the properties of dense quark (baryon) matter, which can exist in the cores of compact stars or arise in heavy ion collisions, the perturbative QCD methods are not applicable because the coupling constant of strong interactions is quite large in these cases. Another well-known (and nonperturbative) method of QCD research, the lattice approach, is also not applicable in this case due to a notorious sign problem. Thus, interest in the phenomena occurring in dense baryonic matter serves as a trigger for studying QCD-like models and/or theories that are free from the above disadvantages and contribute to a deeper understanding of the properties of the phase diagram of matter at nonzero temperature and baryon density. Among such theories, for example, there are Nambu-Jona-Lasinio (NJL)-like models [1] (see also the reviews [2][3][4][5][6][7][8] and references therein), the QCD and QCD-like models with two-color quarks [9][10][11][12][13][14][15][16][17][18][19][20][21], etc.
In addition to the non-zero density of the baryon charge, quark matter can have a nonzero density of isospin (isospin asymmetry), which is an evident characteristic of neutron stars. Such a phenomenon as the charged pion condensation (PC) is closely related to the isotopic asymmetry of matter. And the question of whether the charged PC phase can JHEP06(2020)148 exist in dense quark matter has been discussed . Note that usually dense baryonic matter with isospin imbalance is described in terms of baryon µ B -and isospin µ I chemical potentials.
Quite recently, it was understood that chiral asymmetry (or chiral imbalance), i.e. unequal densities n L and n R of all left-and all right-handed quarks, is also one of the properties of dense quark matter. Usually, chiral asymmetry is characterized by a quantity n 5 called chiral density, n 5 ≡ n R − n L . It can be generated dynamically at high temperatures, for example, in the fireball after heavy ion collision, by virtue of the Adler-Bell-Jackiw anomaly and quarks interacting with gauge (gluon) field configurations with nontrivial topology, named sphalerons. In the presence of an external strong magnetic field, which can be produced in heavy ion collisions as well, this can lead to the so-called chiral magnetic effect [47]. Moreover, in the presence of external magnetic field chiral density n 5 can be produced (even at rather low temperature) in dense quark matter due to the so-called chiral separation effect [48] (it can be also produced under fast rotations of the system due to the so-called chiral vortical effect). Now, let us notice that usually when one talk about chiral density n 5 one implies that chiral density n u5 of u quarks and chiral density n d5 of d quarks are equal to each other (it is evident that n 5 = n u5 + n d5 ). Indeed, that is the case when one has in mind the mechanism of generation of chiral imbalance at high temperatures due to nontrivial topology of gauge field configuration. In this case it is quite plausible that n u5 = n d5 due to the fact that gluon field interacts with different quark flavors in exactly the same way and does not feel the difference between flavors. But another mechanism, the chiral separation effect, is sensitive to the flavor of quarks (as it was shown in appendices A to the papers [49][50][51]). So in dense quark matter a strong magnetic field separates u and d quarks in different ways. As a result, we see that, e.g., in such astrophysical objects as magnetars there might exist areas, in which the quantity n I5 = n u5 − n d5 , called the chiral isospin density, is not zero. Moreover, it has been argued that chiral imbalance is generated by parallel magnetic and electric fields [52,53], one can generalize these arguments to chiral isospin imbalance as well.
So in the most general case, chiral imbalance of quark matter is described by two chemical potentials, chiral µ 5 and chiral isospin (or isotopic) µ I5 chemical potentials, which are thermodynamically conjugated to n 5 and n I5 , respectively. The first, µ 5 , is usually used when isospin asymmetry of quark matter is absent, i.e. in the case µ I = 0 [54][55][56][57][58][59][60][61][62][63][64][65]. The second, µ I5 , might be taken into account when, in addition to chiral, there is also isotopic asymmetry of matter, in which charged PC phenomenon may occur, etc. [66]. In particular, it was established in the framework of NJL models that µ 5 catalyzes the chiral symmetry breaking (CSB) [61][62][63][64][65][66], whereas µ I5 promotes the charged PC in dense quark matter [66]. It should be noted once more that chiral asymmetry of baryonic matter can occur under the influence of a strong magnetic field due to chiral separation and chiral vortical effects, i.e. chiral imbalance is an inevitable characteristic of a dense baryonic (quark) medium in compact stars, as well as in the collision of heavy ions.
In our recent studies [49][50][51][66][67][68][69][70], we showed that chiral asymmetry contributes to the generation of charged PC phenomenon in dense quark matter. Moreover, as it was demonstrated in these papers, in the presence of isospin (µ I = 0) and chiral isospin JHEP06(2020)148 asymmetry (µ I5 = 0) there is a duality between CSB and charged PC phenomena. It means that there is a symmetrical arrangement of these phases on the full phase diagram of quark matter. Our recent studies were carried out in the framework of the NJL model in which the diquark interaction channel, and hence the color superconductivity (CSC) phenomenon, was not taken into account. It is well known that CSC can appear in quark medium at densities several times higher than the density of ordinary nuclear matter, i.e. for µ B 1 GeV (see, e.g., in reviews [71][72][73][74][75] devoted to the CSC phenomenon). Therefore, despite the fact that our results are in well agreement with the lattice QCD analysis [49][50][51], we can trust the duality properties of dense quark matter only in the range of rather low values of the baryonic chemical potential µ B < 1 GeV. On the other hand, inside of neutron stars baryon densities can reach rather large values (corresponding to higher values of µ B ), which means that CSC can be present there. Hence, it is interesting to study how the chiral asymmetry of baryonic matter, which, as mentioned above, is one of the essential features of compact stars with a strong magnetic field, correlates with the phenomenon of color superconductivity. Just the goal of the present paper is to shed new light on the role of both isospin and chiral isospin asymmetries in the formation of diquark condensation, i.e. color superconductivity, at zero temperature. At the same time, to simplify the analysis and avoid unnecessary computational difficulties, we use an approach based on the NJL model, in which quarks are two color.
As an additional motivation for the above proposed two-color approach to quarks, it can be noted that in the QCD with two-color gauge symmetry there is no sign problem of the Dirac operator determinant, so in the framework of this theory first principle lattice calculations can be used to study dense quark-hadron matter (see, e.g., in refs. [9,10]). Moreover, an analysis of dense hadron matter carried out in the two-color NJL model by the mean-field method, showed a good agreement with corresponding lattice results [12-17, 20, 21]. Since in the two-color NJL model not only mesons, but also the simplest colorless baryons are composed of quarks, we have a good opportunity to simulate the meson-baryon interaction both in vacuum and in dense medium, relying on a uniform quark model. Furthermore, in this case the properties of mesons, baryons as well as different phenomena of dense hadron matter, such as diquark condensation, etc, can be investigated with the help of well-known analytical methods. In principle, these results can then be generalized (or extrapolated) to real physical situation with three color quarks, etc.
The paper is organized as follows. In section 2 a (3+1)-dimensional NJL model with two quark flavors (u and d quarks), each is a color doublet, is introduced. It contains both quark-antiquark and diquark channels of interaction and, as it was shown in ref. [12], is a low energy effective model of the two-color QCD. Then, in order to simulate the properties of real dense quark-hadron matter in the framework of this model, we introduce three chemical potentials, baryon µ B , isospin µ I and chiral isospin µ I5 and find in the meanfield approximation the thermodynamic potential (TDP) dependent on chiral, diquark, and charged pion condensates. In section 3 three duality properties (dual symmetries) of the model TDP are established. Each duality property of the model means that its TDP is invariant under some interchange of chemical potentials as well as simultaneous interchange of corresponding condensates. Using dual symmetries of the TDP, in this section a phase JHEP06(2020)148 structure for some simple chemical potential contents of the model is established. For example, if the phase structure of the model is known at µ B = 0 and µ I = µ I5 = 0 [12], then using a dual mapping it is possible to find its phase structure without any calculations at µ I5 = 0 and µ B = µ I = 0, etc. In section 4 the general (µ B , µ I , µ I5 )-phase diagram of the model is analyzed on the basis of dual symmetries of the TDP. In particular, it is shown here that the inclusion of a diquark interaction channel in the analysis does not destroy the dual symmetry of dense quark matter between CSB and charged PC phases, previously found in refs. [49][50][51][66][67][68][69][70]. Moreover, the new symmetries between CSB, charged PC and diquark condensation phases of quark matter appear in the general phase portrait of the model. In section 5 summary and conclusions are given. Some technical details are relegated to three appendices.
2 Two-color (3+1)-dimensional NJL model and its thermodynamic potential It is well known that in the chiral limit the two-color QCD with u and d quarks is symmetrical not only with respect to the usual flavor SU(2) L × SU(2) R group, but also with an enlarged flavor SU(4) group of transformations of quark fields [9][10][11]. Since the SU(4) symmetry connects quarks and antiquarks, it is usually called Pauli-Gursey symmetry. In two-color QCD the colorless baryons are formed by two quarks, i.e. baryons in this theory are bosons.
To obtain an effective 4-fermion NJL model that would reproduce the basic low-energy properties of the initial two-color QCD theory and had the same symmetry group, the authors of ref. [12] integrated out the gluon fields in the generated functional of the QCD. Then, after "approximating" the nonperturbative gluon propagator by a δ−function, one arrives at an effective local chiral four-quark interaction of the form (color current)×(color current) of the NJL type describing low-energy hadron physics. Finally, by performing a Fierz transformation of this interaction term and taking into account only scalar and pseudo-scalar (qq)-as well as scalar (qq)-interaction channels, one obtains a four-fermionic model given by the following Lagrangian (in Minkowski space-time notation) 1 L =q i∂ − m 0 q + H (qq) 2 + (qiγ 5 τ q) 2 + qiγ 5 σ 2 τ 2 q c q c iγ 5 σ 2 τ 2 q . (2.1) In (2.1),∂ ≡ γ ρ ∂ ρ ; q c = Cq T , q c = q T C are charge-conjugated spinors, and C = iγ 2 γ 0 is the charge conjugation matrix (the symbol T denotes the transposition operation). The quark field q ≡ q iα is a flavor and color doublet as well as a four-component Dirac spinor, where i = 1, 2 or u, d; α = 1, 2. (Latin and Greek indices refer to flavor and color indices, respectively; spinor indices are omitted.) Furthermore, we use the notations τ ≡ (τ 1 , τ 2 , τ 3 ) and σ 2 for usual Pauli matrices acting in the two-dimensional flavor and color spaces, respectively. The Lagrangian (2.1) is invariant under color SU(2) c and baryon U(1) B

JHEP06(2020)148
symmetries. Moreover, at m 0 = 0, it has the same Pauli-Gursey flavor SU(4) symmetry as two-color QCD. Generally speaking, Lagrangian (2.1) describes physical processes in a vacuum. In order to study the physics of dense quark matter, it is necessary to add in eq. (2.1) several terms with chemical potentials, The Lagrangian L dense contains baryon µ B -, isospin µ I -and chiral isospin µ I5 chemical potentials. In other words, this model is able to describe the properties of quark matter with nonzero baryon n B = (n u + n d )/2 ≡ n/2, isospin n I = (n u − n d )/2 and chiral isospin n I5 = (n u5 − n d5 )/2 densities, which are the quantities thermodynamically conjugated to chemical potentials µ B , µ I and µ I5 , respectively. (Above we use the notations n f and n f L(R) for density of quarks as well as density of left(right)-handed quarks with individual flavor f = u, d, respectively. Moreover, n f 5 = n f R − n f L .) Note that at m 0 = 0 the Lagrangian (2.2), due to the terms with chemical potentials, is no longer invariant with respect to Pauli-Gursey SU(4) symmetry. So, in the chiral limit the Lagrangian (2.2), apart from the color SU(2) c , is invariant with respect to the U(1) 3) Moreover, the quantities n B , n I and n I5 are the ground state expectation values of the densities of conserved charges corresponding to U(1) B , U(1) I 3 and U(1) AI 3 symmetry groups. So we have from (2.3) that n B = qγ 0 q /2, n I = qγ 0 τ 3 q /2 and n I5 = qγ 0 γ 5 τ 3 q /2. We would like also to remark that, in addition to (2.3), the Lagrangian (2.2) is invariant with respect to the electromagnetic U(1) Q group, U(1) Q : q → exp(iQα)q, (2.4) where Q = diag(q u , q d ) (here q u , q d are the electric charges of u and d quarks, respectively). It is natural to assume that in this theory u and d quarks have a baryon charge equal to 1/2, and their electric charges are 1/2 and (−1/2), correspondingly [13]. In this case scalar colorless diquarks with unit baryon charge have zero electric charge. Their condensation leads to a phase, which is naturally called the (baryon) superfluid phase (by analogy with the phenomenon of superfluidity of helium). At the same time, in some papers on two-color QCD, it is assumed that the electric charge of u quark is 2/3, and that of the d quark is (−1/3) [14]. In this case colorless diquarks of the model have electric charge 1/3. However, since in our consideration both the electromagnetic interaction between quarks and strong external magnetic (electric) fields are ignored, the results of the study of the phase structure of the model do not depend on how to interpret the electric charges of quarks. The ground state expectation values of n B , n I and n I5 can be found by differentiating the thermodynamic potential (TDP) of the model (2.2) with respect to the corresponding chemical potentials. The goal of the present paper is the investigation of the ground state JHEP06(2020)148 properties (or phase structure) of the system (2.2) and its dependence on the chemical potentials µ B , µ I and µ I5 .
Introducing the Nambu-Gorkov bispinor field Ψ, where one can present the auxiliary Lagrangian (2.5) in the following form where matrix Y is given in eq. (2.7) and (2.9)

JHEP06(2020)148
Notice that matrix elements of the 2×2 matrix Z, i.e. the quantities D ± , K and K * , are the nontrivial operators in the (3+1)-dimensional coordinate, four-dimensional spinor, 2dimensional flavor and (N c = 2)-dimensional color spaces. Then, in the one fermion-loop (or mean-field) approximation, the effective action S eff (σ, π, ∆, ∆ * ) of the model (2.2) (this quantity is the generating functional of a one-particle irreducible Green functions of bosonic fields (2.6)) is expressed by means of the path integral over quark fields: where N is a normalization constant and The quark contribution to the effective action, i. e. the term S eff in (2.11), is given by: where the last equality is valid due to the evident relation det Y = 1. Then, using the eqs. (2.11) and (2.13) one can obtain the following expression for the effective action (2.11): (2.14) Starting from (2.14), one can define in the mean-field approximation the thermodynamic potential (TDP) Ω(σ, π, ∆, ∆ * ) of the model (2.2)-(2.5) at zero temperature, The ground state expectation values (mean values) of the fields, σ(x) ≡ σ, π(x) ≡ π, ∆(x) ≡ ∆, ∆ * (x) ≡ ∆, are solutions of the gap equations for the TDP Ω (below, in our approach all ground state expectation values σ, π, ∆, ∆ * do not depend on coordinates x): Since the matrix Z in eq. (2.14) has an evident 2×2 block structure (see in eq. (2.9)), one can use there a general formula

JHEP06(2020)148
and find that (taking into account the relation τ 2 τ τ 2 = − τ T and assuming that all bosonic fields do not depend on x) where (in the following we use also the notations µ ≡ µ B /2, ν ≡ µ I /2, Obviously, the quantity which is in the square brackets of eq. (2.18) is proportional to the unit operator in the N c -color space. (Below, in all numerical calculations we put N c = 2.) Hence, where D is the 2×2 matrix in the 2-dimensional flavor space (its matrix elements D kl are the nontrivial operators in the 4-dimensional spinor and in the (3+1)-dimensional coordinate spaces). Using this expression for det(Z) in eq. (2.14) when σ, π, ∆, ∆ * do not depend on coordinates x, and taking into account the general relation (A.6) of appendix A, we find where the 2×2 matrix D(p) is the momentum space representation of the matrix D of eq. (2.20). Its matrix elements D kl (p) have the following form where M ≡ m 0 + σ. Using in eq. (2.21) again the general relation (2.17), we have

JHEP06(2020)148
Notice that the matrix L(p), i.e. the expression in square brackets of eq. (2.23), is indeed a 4×4 matrix in 4-dimensional spinor space only, which is composed of 4×4 matrices D ij (p) (see in eq. (2.22)). Now, taking into account the definition (2.15) and using the eqs. (2.21)-(2.23), it is possible to obtain in the mean-field approximation and at N c = 2 the following expression for the TDP of the model, . . , 4) are four eigenvalues of the 4×4 matrix L(p), in the following, in order to find the TDP of the model in various cases, first of all, we will find the eigenvalues of the matrix L(p). Then, after integration in eq. (2.24) over p 0 , this TDP is used in some numerical calculations with sharp threemomentum cutoff Λ = 657 MeV (i.e. it is assumed below that the integration over threemomentum p occurs over the region | p| < Λ) at H = 7.23 GeV −2 and m 0 = 5.4 MeV [15,16]. In order to simplify numerical calculations, we study the phase structure of the model in the chiral limit, m 0 = 0, at the same values of Λ and H.

Calculation of the TDP (2.24) and its duality properties
3.1 The case µ = 0, but ν = 0 and ν 5 = 0 First of all let us consider the TDP (2.24) in the most simple case when only µ = 0. In this case, using any package of analytical calculations, it is easy to establish that the 4×4 matrix L(p) of eqs. (2.23) and (2.24) has a unique fourfold degenerated eigenvalue λ(p), (3.1) (It means that L(p) is proportional to the identity 4×4 matrix.) As a result, in this case the TDP (2.24) looks like Taking into account the general formula (being true up to an infinite term independent of the real quantities a and b), it is possible to integrate in (3.2) over p 0 . Then, the TDP takes the form
The case m 0 = 0. Below, just in the chiral limit, we will mainly study the properties of the NJL model (2.2). Therefore, notice that at m 0 = 0 and for all values of µ such that 0 < µ <μ c ≈ 0.88 GeV the BSF phase is realized. Based on this phase portrait and using the dual properties of the model, we are able to predict the phase structure of the model for other special cases of chemical potentials (see the next sections).

3.2
The case with µ = 0, ν = 0 and ν 5 = 0 In the previous sections we have considered the TDP (2.24) of the initial model as a function of six variables (condensates), M , π, ∆ and ∆ * . But due to a symmetry of the model, the number of condensates may be reduced. Indeed, at m 0 = 0 and zero chemical potentials the If chemical potentials are taken into consideration, then in the chiral limit the symmetry of the Lagrangian (2.2)

JHEP06(2020)148
reduces to U(1) B ×U(1) I 3 ×U(1) AI 3 symmetry (plus color SU(2) c , which in our consideration is not violated at all). As a result, we see that at m 0 = 0 and nonzero chemical potentials the TDP of the model depends only on the |∆| 2 , π 2 1 + π 2 2 and M 2 + π 2 0 field combinations, correspondingly. It is important to note that the matrix L(p) as well as the one-loop contribution to the TDP, i.e. the last term in eq. (2.24), do not depend on the bare mass m 0 . Hence in any case, both at m 0 = 0 and m 0 = 0, these quantities depend only on the above mentioned combinations of bosonic fields. As a result, we see that in the chiral limit one can put π 0 = 0 and π 2 = 0 without loss of generality of consideration. However, at m 0 = 0 and at nonzero chemical potentials the symmetry of the Lagrangian (2.2) reduces to U B (1) × U I 3 (1), i.e. in this case the TDP depends on |∆| 2 , M, π 0 and π 2 1 + π 2 2 . So in this case without loss of generality we can also use the symmetry properties of the model and put π 2 = 0. Moreover, since at m 0 = 0 the discussion presented just after the eq. (3.8) is valid, we can attract these dynamical reasons and put π 0 = 0 in addition. Hence, below throughout the paper we suppose that the TDP (2.24) is a function of only M, π 1 and |∆|.
In the case when µ = 0, ν = 0 and ν 5 = 0 the matrix L(p) of eqs. (2.23) and (2.24) has four different eigenvalues λ i (p) (they can be found with the help of any package of analytical calculations), 12) and the TDP (2.24) has the form It is easy to see from eqs. (3.10)-(3.12) that each of the eigenvalues λ i (3.9) is invariant with respect to the so-called dual transformation D 1 , (3.14)

JHEP06(2020)148
As a result, the whole TDP (3.13) is also invariant under the transformation (3.14). (For the first time, the discrete symmetry (3.14) of the thermodynamic potential of the NJL model with two-color quarks was noted, but only at ν 5 = 0, in refs. [11,16]. It was also shown there that due to this symmetry, the charged PC and BSF phases are located symmetrically on the (µ, ν)-phase diagram of the model.) Moreover, using any program of analytical calculations, it is possible to establish that the quantities λ 1 (p)λ 2 (p) and λ 3 (p)λ 4 (p) are invariant in addition with respect to the dual discrete transformations D 2 and D 3 , where Indeed, it is possible to show that Note that in the general case when m 0 = 0 the investigation of the TDP (3.13) is a rather hard task. So in the following we study the NJL model only in the chiral limit, m 0 = 0, as a rule. It is possible because in general the chiral limit often is a good approximation. In particular, this conclusion is supported by the results of refs. [49][50][51], where it was shown that for the three-color NJL model with µ, ν, ν 5 chemical potentials and m 0 = 0, the phase structure obtained in the chiral limit gives a rather good approximation for the phase portrait, except the region of low chemical potentials, smaller than the value of the half of the pion mass. Moreover, the exact dual symmetry D 3 (3.15) of this massless JHEP06(2020)148 three-color NJL model is a very good approximation at m 0 = 0. Therefore, we think that the two-color NJL model (2.2) can be investigated with a sufficient degree of reliability in the chiral limit, as well.
Let us argue that in this case it is possible to assume that for sufficiently low values of the chemical potentials (say at µ, ν, ν 5 < 1 GeV) at the global minimum point (GMP) (M, π 1 , |∆|) of the TDP (3.13), there can be no more than one nonzero coordinates, i.e. condensates or order parameters. The argument could be that in previous investigations there have not been found (mixed) phases with several nonzero condensates, e.g., in the three-color NJL model [49][50][51] and in the two-color one [11,16]. But it is not overwhelmingly convincing argument anyway. There is a more persuasive argument indicating that it is the case. Dualities stand for that. It will be discussed below that the phase structure in terms of chiral symmetry breaking and charged pion condensation is the same for two-color and three-color NJL models. The phase structure of the three-color NJL model has been studied in refs. [49][50][51][66][67][68][69][70]. And it has been shown that there is no mixed phase with simultaneously non-zero chiral and charged pion condensates. So, one can conclude that there is no such a mixed phase in two-color case as well. Then one can apply the duality D 1 (3.14) and argue that there is no mixed phase, where there are simultaneously non-zero chiral and diquark condensates. Or use the duality D 2 (3.15) and show that there is no mixed phase, where there are simultaneously non-zero pion and diquark condensates. 2 Let us note that this does not exclude the case where all three (chiral, pion and diquark) condensates are simultaneously non-zero, but let us be honest it is quite unlikely, especially considering the fact that mixed phases are not common. So it was rigorously proved that there is no mixed phases with only two non-zero condensates and one can undisturbed assume that there are no mixed phases (mixed phase with three condensates are extremely unlikely and even if it is present it would probably be small region in the phase diagram). To show the absence of mixed phases with two non-zero condensates without using the dualities, would require a way more time-consuming numerical calculations.
Therefore, with such a restriction on chemical potentials, in the chiral limit only four different phases can be realized in the system. (i) If GMP has the form (M = 0, π 1 = 0, |∆| = 0), then the chiral symmetry breaking (CSB) phase appears in the model. (ii) If it has the form (M = 0, π 1 = 0, |∆| = 0), the charged pion condensation (PC) phase is realized. (iii) When the GMP looks like (M = 0, π 1 = 0, |∆| = 0), it corresponds to 2 Strictly speaking, the statement that it is impossible to have two or three condensates that are not equal to zero in the ground state of the system refers to the points of a phase diagram that do not lie on the boundary between phases. In the last case there might occurs another situation. For example, in refs. [18,19] it was shown in the framework of the chiral perturbation theory for two-color QCD that at µB = µI (in figure 1(a) it is the boundary between charged PC and BSF phases) there is a multicomponent superfluid phase in which two condensates are nonzero, π1 = 0 and |∆| = 0. The same is valid in the framework of the NJL model (2) as well. Indeed, it is easy to check that at µB = µI the eigenvalues λi (3.9) and the TDP (3.13) depend on the combination π 2 1 + |∆| 2 but not on π1 and |∆| separately. Hence, in this case in the GMP of the TDP (3.13) the relations π1 = 0 and |∆| = 0 can be realized. But in the present paper we do not study the properties of the system on the boundaries between phases, so the corresponding exceptional sets of chemical potentials, e.g., such as µB = µI , are outside the scope of our consideration.

Consequences of dual symmetries of the TDP in some simple cases
By the duality property (or symmetry, or relation, etc) of any model, we understand any discrete symmetry of its TDP with respect to transformations as order parameters (in our case, condensates M , π 1 and |∆|) and free external parameters of the system (these may be chemical potentials, coupling constants, etc). The presence of the dual symmetry of the model TDP means that in its phase portrait there is some symmetry between phases with respect to the transformation of external parameters, which can greatly simplify the construction of the full phase diagram of the system. (The invariance of the TDP (3.5) under the changing of a sign of the quark chemical potential µ is the simplest example of the dual symmetry of the model (2.2). Due to this kind of duality, it is enough to JHEP06(2020)148 study the phase structure of the model only, e.g., at µ > 0 when other chemical potentials are zero, etc.) Below, we investigate the phase portrait of the model (2.2) in the meanfield approximation in the presence of three nonzero chemical potentials, µ, ν and ν 5 , in the chiral limit. In this case, the problem is greatly simplified due to the fact that the TDP (3.13) of the model has three dual symmetries, D 1 (3.14) and D 2 , D 3 (3.15). Indeed, let us suppose that m 0 = 0 and that at the point (µ = a, ν = b, ν 5 = c) of the phase portrait the GMP of the TDP (3.13) lies, e.g., in the point of condensate space of the form (M = A, π 1 = 0, |∆| = 0), i.e. in this case the CSB phase is realized in the system. Then, according to the symmetries D 2 and D 3 (3.15), the TDP will have the same meaning if we interchange the values of chemical potentials and simultaneously appropriately transpose the values of the condensates. As a result we see that, e.g., at µ = c, ν = b, ν 5 = a and in the point (M = 0, π 1 = 0, |∆| = A) (it is the result of the action of the D 2 dual transformation on the TDP) as well as that at µ = a, ν = c, ν 5 = b and in the point (M = 0, π 1 = A, |∆| = 0) (it is the D 3 dual transformation of the TDP) it has the initial meaning. Moreover, it is evident that these new points of the condensate space are nothing but the GMPs of the TDP after its D 2 and D 3 transformations. 3 Consequently, at the points (µ = c, ν = b, ν 5 = a) and (µ = a, ν = c, ν 5 = b) of the phase diagram of the model, which we call dually D 2 and dually D 3 conjugated to the starting point (µ = a, ν = b, ν 5 = c) of the phase portrait, there are BSF and charged PC phases that are respectively dually D 2 and dually D 3 conjugated to the initial CSB phase of the model. Thus, knowing the phase of the model, which is realized at some point of its phase portrait, we can predict which phases are arranged at the dually conjugated points of a phase diagram.
At m 0 = 0 each duality transformation D i (i = 1, 2, 3) (3.14) and (3.15) of the TDP can also be applied to an arbitrary phase portrait of the model as a whole. In particular, it is clear that if we have a most general (µ, ν, ν 5 )-phase portrait, then the action, e.g., of the D 3 on the TDP can be understood as the following dual D 3 transformation of the model (µ, ν, ν 5 )-phase portrait. It is necessary to rename both the diagram axes and phases in such a way, that ν ↔ ν 5 and CSB↔charged PC. At the same time the µaxis and BSF and symmetrical phases should not change their names. It is evident that after such D 3 transformation the (µ, ν, ν 5 )-phase diagram is mapped to itself, i.e. the most general (µ, ν, ν 5 )-phase portrait is self-D 3 -dual. Furthermore, the self-D 3 -duality of the (µ, ν, ν 5 )-phase portrait means that in the three-dimensional (µ, ν, ν 5 ) space the regions of the CSB and charged PC phases are arranged mirror-symmetrically with respect to the plane ν = ν 5 of this space. In a similar way it is possible to describe the action of other, D 1 and D 2 , duality transformations on the (µ, ν, ν 5 )-phase portrait of the model, which is, of course, invariant, or self-dual, under these mappings. But different cross-sections of the full (µ, ν, ν 5 )-phase diagram, e.g., the (µ, ν)-phase portrait at some fixed value of   figure 1(a). In all panels BSF means the baryon superfluid (or diquark condensation) phase in which M = 0, π 1 = 0, |∆| = 0. PC is the shorthand for the charged pion condensation phase with M = 0, π 1 = 0, |∆| = 0. In the chiral symmetry breaking (CSB) phase M = 0, π 1 = 0, |∆| = 0. SYM is the notation for the symmetrical phase, where M = 0, π 1 = 0, |∆| = 0. phase remains intact, i.e. it does not change its position on the phase diagram.
As a result, based on this mechanism of dual transformation of phase diagrams and using previously obtained particular (simple) phase portraits of the NJL model (2.2) with two quark colors, it is possible to describe its phase structure at qualitatively different chemical potential contents at m 0 = 0. Namely, • For example, as it is noted at the end of subsection 3.1, the phase structure of the NJL model (2.2), when only quark chemical potential µ is nonzero, was studied in ref. [12] at m 0 = 0 and nonzero temperature T . Reducing this analysis to the case m 0 = 0, one can see that at 0 < µ < 0.88 GeV and T = 0 the BSF phase is realized in the model. It means that at µ > 0, ν = 0 and ν 5 = 0 the condensates are the following, (M = 0, π 1 = 0, |∆| = 0). Applying to this phase diagram the D 1 duality tranformation (3.14), we replace µ by ν and BSF phase by charged PC phase and find the phase portrait of the model at ν > 0, µ = 0 and ν 5 = 0. Hence, we see that if only isospin asymmetry is presented in the model (i.e. only ν = 0), then in the chiral limit and for sufficiently low values of ν the charged PC phase is observed. However, acting on the original phase portrait by duality transformation D 2 (3.15), one can find that if ν 5 > 0 (and small enough), then usual CSB phase is realized in the model in the chiral limit.
• There is another simple phase diagram of the model (2.2) from which it is possible to find its phase structure for qualitatively different chemical potential sets, using only dual transformations of this phase diagram (and without any numerical calculations). Earlier, the (µ, ν)-phase portrait of the model was obtained in refs. [11,16] at ν 5 = 0 and m 0 = 0. In figure 1(a) this phase portrait is presented in the chiral limit, and it is clear that there, due to the D 1 -invariance (3.14) of the model TDP (3.13), the BSF and charged PC phases of the model are arranged mirror-symmetrically with respect JHEP06(2020)148 to the line µ = ν (it is one of the conclusions of the papers [11,16]), i.e. this diagram is a self-dual with respect to D 1 mapping. Then, performing in each point of this diagram the D 2 dual transformation of the TDP, one can get the so-called D 2 dual conjugation (or mapping) of this (µ, ν)-phase portrait at ν 5 = 0, which is no more than the (ν 5 , ν)-phase portrait of the model at µ = 0 (see in figure 1(b)). In a similar way, applying to the diagram of figure 1(a) the D 3 dual transformation, one can obtain its D 3 dual conjugation, i.e. the (µ, ν 5 )-phase portrait of the model at ν = 0, where instead of the charged PC the CSB phase appears. And mirror-symmetrically to it (with respect to the line ν 5 = µ) the BSF phase is arranged. Note that under any dual transformation D i the position of the symmetrical phase is not changed.
Taking into account these simple phase diagrams, e.g., presented in figure 1, it is possible to say that different phenomena are generated by different chemical potentials. So, charged PC is caused by isospin chemical potential µ I ≡ 2ν, baryon chemical potential µ B ≡ 2µ leads to the diquark condensation in the system, whereas at non-zero chiral isospin chemical potential µ I5 ≡ 2ν 5 the chiral symmetry gets broken.
4 Phase structure in the general case µ = 0, ν = 0 and ν 5 = 0 In the present section we study the general (µ, ν, ν 5 )-phase portrait of the model in the chiral limit. In this case the investigation is simplified due to the dual symmetries D i (3.14) and of it, it is very convenient to consider different cross-sections of this general (µ, ν, ν 5 )-phase portrait by the planes of the form µ = const, ν = const and ν 5 = const. Below, these different cross-sections of the most general (µ, ν, ν 5 )-phase portrait will be presented in the chiral limit. Moreover, the solution of the problem is greatly facilitated because of the dual symmetries of the thermodynamic potential. Indeed, it is sufficient to construct, for example, (µ, ν)-phase diagrams for some fixed values of ν 5 , and then apply duality transformations D 2,3 (3.15) to them according to the scheme described in the subsection 3.3. As a result, one can obtain the missing (ν, ν 5 )-and (µ, ν 5 )-phase diagrams at the corresponding fixed values of µ and ν, respectively.  3(a,b)). Their characteristic feature is the fact that the diquark condensation BSF phase is located in the region, where µ > ν, ν 5 , and, in contrast, the charged PC phase is located at ν > µ, ν 5 . For simplicity, let us assume that ν 5 ≈ 0. Then this property of the model can be explained qualitatively using the same arguments as in ref. [16], where the phase structure of this model has been investigated at ν 5 = 0. Indeed, if µ > ν, then both u and d quarks form their own Fermi seas. As a result, the system gives the possibilities of the formation of ud particle-particle Cooper pairs with quantum numbers of auxiliary field ∆(x) (see in eq. (2.6)), their condensation, and therefore of the formation of the ground state of the BSF phase. At ν > µ, u andd quarks form Fermi seas. This leads both to the possibility of the formation and condensation of ud Cooper pairs, and to the formation of the charged PC phase in the system.
• The situation is drastically changed in the case when ν 5 is rather large, e.g., at ν 5 greater than 0.25-0.3 GeV. In this case the typical (µ, ν)-phase portrait is presented in Fig 4(a), where we fix ν 5 to be 0.6 GeV. But to get the idea of the phase structure at other rather high fixed values of ν 5 , one can imagine that the charged PC phase (the BSF phase) of this diagram has the form of a boot sole pointing by its tip to the value µ = ν 5 = 0.6 GeV of the µ-axis (to the value ν = ν 5 = 0.6 GeV of the ν-axis). (And in the figure, these "boot soles" intersect with each other, as well as with the strip of the CSB phase along the line µ = ν.) Then, to find (approximately) the (µ, ν)-phase diagram at another rather high value of ν 5 = 0.6 GeV, one should, starting from the diagram of figure 4(a), simply shift the "boot sole" of the charged PC phase parallel to the µ-axis in a position, in which its tip is directed to the point of this axis, where µ = ν 5 . In a similar way the BSF phase should be shifted along the ν-axis. Note, that already in figures 3(a,b) one can see the beginning of the process of forming the charged PC and BSF phases in the form of the "boot soles". At rather large value of ν 5 , say at ν 5 = 0.4 GeV, a small part of the BSF phase, a tip of the "boot sole", penetrates into the region where µ < ν (see in figure 3(c)). And at even higher values of ν 5 in all the (µ, ν)-phase diagrams obtained in this way, the BSF phase (the charged PC phase) will be located mainly in the region µ < ν ≈ ν 5 (in the JHEP06(2020)148 To explain this qualitative difference in the phase diagrams, it must be borne in mind that for large values of ν 5 , the Fermi seas and surfaces for left-and right-handed quarks of the same flavor differ significantly from each other. Indeed, it was shown in our previous paper [66] that in the presence of µ, ν, ν 5 the chemical potentials of right-and left-handed u and d quarks look like It is clear from (4.1) that at µ < ν ≈ ν 5 the chemical potential µ dR of the d R quarks is negative. So the Fermi sea with positive energies can be formed by charge conjugated d c R quarks, and just the d c R quark with some momentum p might be born above the corresponding Fermi surface. The chemical potential µ uR of right-handed u R quarks is positive, so their Fermi sea also exists, over the surface of which the u R quark with the same momentum p (as for the d c R quark) can also be born. But in the last case inside the u R -Fermi sea the holeū R with the momentum (− p) appears, and theū R d c R particle-hole Cooper pair with quantum numbers of the auxiliary field ∆ * (x) (see in eq. (2.6)) can be created. Condensation of these Cooper pairs leads to a rearrangement of the ground state of the system and its transition to the BSF phase. Note once more that in the region µ > ν, ν 5 and rather small values of ν 5 (as, e.g., in figures 2(a) and 3(a,b)) the diquark condensation BSF phase appears due to a particle-particle Cooper pairing, whereas at µ < ν ≈ ν 5 and rather high values of ν 5 (see in figure 4(a)) the BSF phase exists due to a particle-hole Cooper pairing mechanism.
In each of figures 1(a), 2(a), 3 and 4(a) it is easy to see the duality D 1 (3.14) between charged PC and BSF phases, i.e. their symmetrical arrangement with respect to the line µ = ν. Moreover, it is clear, for example, from figure 4(a) that in quark matter, which is characterized by chiral imbalance with ν 5 = const, the diquark condensation BSF phase can appear even at rather low values of quark number chemical potential µ 0.2 GeV, but only under condition that ν ≈ ν 5 . In contrast, in this case the charged PC phenomenon can be realized at rather small values of ν if µ ≈ ν 5 . This means that in order for a BSF phase with diquark condensate to appear in dense quark matter with chiral imbalance (ν 5 = 0), an isotopic imbalance with ν ≈ ν 5 of this dense medium is also necessary.
• Applying to these (µ, ν)-phase portraits (each at some fixed value of ν 5 ) the duality transformations D 2 and D 3 (3.15), one can find the corresponding dually conjugated (ν, ν 5 )-and (µ, ν 5 )-phase diagrams at corresponding fixed values of µ and ν, respectively (see, e.g., in figures 2 and 4). So in figure 4 we present the (µ, ν)-phase diagram at fixed ν 5 = 0.6 GeV as well as its dual D 2 mapping, i.e. the (ν, ν 5 )-diagram at fixed µ = 0.6 GeV, and its dual D 3 mapping, which is the (µ, ν 5 )-diagram at fixed ν = 0.6 GeV. Precisely the diagram of figure 4(b) that most clearly supports the assertion of the previous paragraph that just at ν ≈ ν 5 the diquark condensation BSF phase can be generated in quark matter for a rather wide interval of µ. Moreover, it is clear from figure 4(c) that for a wide interval of ν the charged PC phase can be realized in the system if µ ≈ ν 5 . And in the region, where µ ≈ ν ≈ ν 5 , a rather complicated interplay between BSF, CSB and charged PC phases takes place demonstrating a very rich phase structure of the two-color NJL model. In this region several first order phase transitions can occur (see figure 5, where the behavior of the condensates vs chemical potentials are presented). It follows, e.g., from figure 5 (left panel) that at fixed ν = 0.57 GeV, ν 5 = 0.6 GeV and µ less than 0.49 GeV, the BSF phase is realized in the system. Then, with the growth of µ (in this case one goes deeper in the neutron star, or baryon density increases in heavy ion collision, etc), the first order phase transition takes place to the CSB phase, and then, at even larger µ, there is a second first order phase transition to the charged PC phase, after which the system returns to the BSF state again. In figure 5 (right panel), which is dually D 1 conjugated to the left panel of the same figure, the behavior of the gaps vs ν are presented. It can be used to simulate other, dually D 1 conjugate modes of physical processes in dense baryonic medium. Moreover, the right and left panels of figure 5 demonstrate the fact that under the dual D 1 transformation the diquark |∆| and charged pion π 1 condensates pass one into another.

JHEP06(2020)148
• Finally, we would like to note another interesting feature of all phase portraits in figures 1-4, and the two-color NJL model (2.2) itself. In ref. [66] the phase structure of the NJL model with three-color quarks was studied in the mean-field approximation and the same set of chemical potentials, µ, ν, ν 5 . There, in contrast with the twocolor NJL model (2.2), only theqq interaction channel was taken into consideration with coupling constant G, where G = 15.03 GeV −2 (i.e. the diquark qq channel was ignored). Moreover, it was shown there that in this model the duality D 3 (3.15) between CSB and charged PC phases exists (in the chiral limit, m 0 = 0). It turns out that at G = 2H the TDP of the three-color NJL model of ref. [66] is the half of the TDP (2.24)-(3.13) of the present two-color NJL model, if the diquark condensate ∆ is zero (see the remark below eq. (C.5) of appendix C), and it can be viewed as sometimes one divides coupling constant G by number of colors. (Indeed, in ref. [66] in the three-color considerations the coupling constant G was divided by number of colors. However, in this paper the coupling constant H was not divided by N c in order to be consistent and have the same fit of parameters as in ref. [16].) In this case, i.e. at G = 2H, it is clear that if for some set µ, ν, ν 5 of chemical potential values the phase with ∆ = 0 is observed in the framework of the two-color NJL model (2.2), then the same phase is realized in the three-color NJL model of ref. [66] (but generally speaking, not vice versa).
Let us discuss the fit of the model parameters in more detail (see also appendix C). Note that the fit in refs. [15,16], which is also used in this paper, was constructed from the fit of three color case by the following N c scaling of physical quantities f π ∼ √ N c and qq ∼ N c . It can be easily shown that it is equivalent to the scaling of coupling constant G ∼ 1/N c . So from the fit used in three-color case of ref. [66], one could directly obtain the fit for two-color NJL model (2.2), H = G/2 = 7.515 GeV −2 and Λ = 650 MeV. However, in order to be consistent with earlier investigations, we use in all numerical calculations throughout the paper the fit employed in the previous considerations (see in refs. [15,16]) of the twocolor case, H = 7.23 GeV −2 and Λ = 657 MeV (compare with the parameter fixing in the two-color NJL model under consideration given just after eq. (2.24)) that corresponds to f π = 75.4 MeV and − qq 1 3 = 218 MeV. Nevertheless, a simple comparison shows that in this case the relation 2H = G between the coupling constants of these two NJL models is performed with a rather high accuracy. Therefore, it can be expected that in each of the diagrams of figures 1-4, the same phase structure is realized outside the region occupied by the BSF phase as predicted by the three-color NJL model of ref. [66]. Comparing the corresponding phase portraits obtained in these two models, we see that this is the case.

Summary and conclusions
It is well known that at low temperatures and rather high baryon densities the appearance of the CSC (or diquark condensation) phenomenon in quark matter is expected [71][72][73][74][75]. The influence of chiral imbalance (when chiral chemical potential µ 5 is nonzero) of dense cold quark matter on CSC was investigated, e.g., in ref. [60], where it was shown that an JHEP06(2020)148 increase in µ 5 , the CSC phase in quark matter is suppressed. This effect can be explained by the fact that the chiral chemical potential µ 5 is a catalyst for spontaneous breaking of chiral symmetry [61][62][63][64][65]. However, there is another mechanism of chiral asymmetry, which can be realized in dense quark matter in strong magnetic fields (see the discussion in Introduction) and is characterized by the chiral isospin chemical potential µ I5 . And it was shown that chiral isospin imbalance (when µ I5 = 0) promotes the generation of the charged PC phenomenon in dense quark matter [49][50][51][66][67][68][69][70] (see also the recent reviews [7,8] on the topic). Of course, this result can be trusted in the region of relatively low baryon densities of quark matter. What happens at high densities when it is necessary to take CSC into account remains unclear.
In the present paper an interplay between various physical phenomena, CSB, charged PC and diquark condensation, in dense (µ B ≡ 2µ = 0) quark matter with isospin (µ I ≡ 2ν = 0) and chiral isospin (µ I5 ≡ 2ν 5 = 0) imbalances is investigated at zero temperature. To simplify the consideration, we study the problem within the framework of the NJL model (2.1)-(2.2), in which quarks are two color. In addition, a diquark interaction channel is also taken into account in this model and, therefore, its thermodynamic potential (TDP) Ω(M, π 1 , ∆) (3.13) depends on three order parameters (or condensates), chiral M , charged pion π 1 , as well as diquark ∆. The last, diquark condensate, both in two-color QCD and in two-color NJL model is a colorless and electrically neutral quantity (if electric charges of quarks q u,d = ±1/2, as in ref. [13]) with nonzero baryon charge. So if in the ground state of the system the diquark condensate ∆ is nonzero, then the baryon superfluid (BSF) phase is realized (it is an analog of helium superfluidity), and the baryon U(1) B symmetry of the model is spontaneously broken. However, if q u = 2/3 and q d = −1/3, as in ref. [14], then diquarks are colorless baryons with electric charge 1/3. Their condensation leads to a phase with spontaneous breaking of the baryon U(1) B and electromagnetic U(1) Q symmetries. This phase resembles the color superconductive phase of real QCD (in three-color QCD the diquarks have in addition color charge [71][72][73][74][75].) Since in the present consideration the JHEP06(2020)148 electromagnetic interaction is not taken into account, the results, i.e. the phase structure of the model (2.1)-(2.2), do not depend on which of these two ways of fixing the electric charges of quarks is chosen. So in our paper the diquark condensation phase is denoted by BSF phase, in each case.
We have calculated the TDP of this two-color NJL model in the mean-field approximation and found that in the chiral limit (bare mass of quarks is zero) it is symmetric under three discrete D 1,2,3 (3.14) and (3.15) simultaneous transformations of both condensates and chemical potentials. Finally note that the consideration has been performed within the assumption, which is confirmed by a numerical analysis in the region of sufficiently small values of µ, ν, ν 5 < 1 GeV, that in the chiral limit there is no mixed phase. It means that in the global minimum point of the TDP two or more condensates cannot have nonzero values (see the discussion in section 3.2).
Let us summarize the main results on our investigation of dense quark matter properties in the presence of isospin and chiral isospin imbalances within the framework of the two-color NJL model (2.2) withqq and qq interaction channels (quarks are massless).
• Considering the (µ, ν, ν 5 )-phase portrait of two-color quark matter as a whole, we see that for m 0 = 0 it has a surprising (dual) symmetry between CSB, charged PC and BSF phenomena. It is a consequence of the invariance of the TDP (3.13) with respect to three discrete D 1,2,3 (3.14) and (3.15) transformations. (One of them, D 3 is similar to the one observed in the three-color case (N c = 3), others exist only when one has two colors.) In the simplest phase portraits, this dual symmetry manifests itself in the symmetric arrangement of different phases in phase diagrams. For example, in figure 4(a) the charged PC and BSF phases are mirror symmetric to each other with respect to the line µ = ν, which is a consequence of the dual invariance D 1 (3.14) of the thermodynamic potential of the system, etc.
• Relying on the dual symmetries D 1,2,3 of the thermodynamic potential, one can obtain qualitatively different phase portraits of the system without any additional numerical calculations. For example, from different (µ, ν)-phase diagrams, each of which with a fixed value of ν 5 , it is possible to obtain the corresponding (ν, ν 5 )-and (µ, ν 5 )phase portraits of the model which are the dually D 2 and dually D 3 conjugated to the initial (µ, ν)-phase diagrams (see in figures 1, 2 and 4), etc. The corresponding technical procedure for constructing a dual mapping of a phase diagram is described in section 3.3.
• Based on the duality properties of the phase portrait of the model (2.2), we have shown that in quark matter composed of two-color quarks, a mixed phase, i.e. phase with two non-zero simultaneously condensates, cannot be realized (see the discussion in the section 3.2). Knowing this greatly simplifies the numeric calculations making them several orders of magnitude less machine-time consuming.
The found dualities are very interesting properties of the phase structure and greatly simplify the investigation of the phase diagram. The dualities interconnect the whole phase JHEP06(2020)148 diagram (different parts of the phase diagram) and the knowledge of only a part of the phase diagram can give us the knowledge about the whole phase structure.
• It has been shown that in a sense each phenomena is inextricably linked with the corresponding chemical potential: µ I is linked with charged PC phenomenon and it leads to its appearance in the system, BSF (or diquark condensation) phenomenon is caused by µ B (in three-color case the CSC appears at large baryon chemical potential), whereas chiral isospin chemical potential is connected with CSB phenomenon (catalysis of chiral symmetry breaking by chiral imbalance shown as in three-color as well as in two-color cases). This picture can be clearly seen in the range of chemical potentials smaller than 200-300 MeV (see, e.g., the phase diagrams of figures 1, 2, 3(a)). At even larger chemical potentials picture changes and becomes more rich and new interesting features appears (see below).
• We show that in quark matter with rather high baryon density the generation of the charged PC phase by ν 5 is not prohibited at all by the processes of diquark condensation (see in figures 1-4). Moreover, it can be seen from these phase diagrams that at low values of ν 5 the charged PC phase is located in the region ν > µ (see, e.g., in figures 1-3), whereas at large ν 5 it can appear even at ν < µ (see in figures 3-4).
• And vice versa, as it is clear from the same diagrams, at rather high values of ν 5 the BSF phase can appear even at sufficiently low values of µ compared with ν, i.e., in a sense, chiral isospin asymmetry of quark matter with isospin density promotes the formation of diquark condensate and shifts the BSF phase to the region of low baryon densities ν > µ, whereas at low ν 5 the BSF phase is realized in the region where ν < µ. Chiral imbalance leads to the generation of BSF phase in quark matter with non-zero isospin density (if isospin density is larger than the baryon density). If there is no chiral imbalance one expects PC phase to appear in this regime (if ν > µ).
• It is interesting also to note that outside of the BSF regions of all phase diagrams in figures 1-4, the predictions of the phase structure, made in the present two-color NJL model (2.2), correspond to a phase structure which is observed in the framework of the three-color NJL model [49][50][51][66][67][68][69][70] without a diquark channel (see the remark at the end of section 4). The converse is generally not true, i.e. if in the NJL model with three-color quarks and without a diquark channel, some phase is realized (it can be CSB, charged PC or symmetric phase) at some set µ, ν, ν 5 = 0, then this does not mean that in the 2-color NJL model (2.2) approach this phase will also be realized at the same values of chemical potentials. The appearance of the phase with diquark condensation ∆ = 0 is also possible, which requires additional consideration.
As it was noted above, the phase diagram in the two-color NJL model at nonzero µ, ν and ν 5 is highly symmetric due to the existent dualities. One can also notice that the phase diagram investigated in the three color case [49][50][51]66] was in a sense of dualities not complete and only one facet of this brilliant diamond. This facet is shared by two and three color cases. There was only one duality and in the case of three colors symmetric JHEP06(2020)148 was only one section of the phase diagram, namely, (ν, ν 5 ) and it was hard to imagine different sections and that sections were not such symmetric. But in the two color case the phase diagram is extremely symmetric and one can easily imagine all the sections and they are very similar (have the same structure and can be mapped between each other by three dual symmetries). It is symmetric to such a degree (to tell the truth it has been only realized after the phase portrait was explored numerically) that the whole phase structure of dense quark matter with isospin and chiral imbalances in the two color case can be easily obtained and comprehended from the part of it (bearing in mind the connection with three color case, from the results of three color NJL model (without considering diquark channel)). So the results concerning diquark condensation in two color case can be obtained just by considerations of three color case without any diquark condensation, which is quite interesting in itself. So one has enormously symmetric phase structure of dense quark matter with isospin and chiral imbalances in the case of two color QCD one facet of which was explored and noticed in the three color case investigations.
Our analysis performed in the framework of two-color NJL model shows that there is a very rich phase structure of dense quark matter with isospin and chiral isospin imbalances. And there could occur a lot of phase transitions that can be interesting in various applications. For example, at high baryon density, e.g., in neutron star or intermediate energy heavy ion collision experiment, etc, several first order phase transitions can occur with rather moderate change of parameters. It is clear from figure 5 where behaviors of the order parameters (i.e. gaps or condensates) are shown.
We performed our study mainly in the chiral limit, m 0 = 0, when the general phase portrait of the model has three dual symmetries D 1,2,3 (3.14)-(3.15). It is worth noting that when m 0 = 0, only the D 1 (3.14) is preserved. In this case, in our opinion, the other two dual symmetries D 2,3 of the phase portrait of two-color quark matter are fulfilled only approximately, and then in the region where µ B , µ I , µ I5 > m π . We draw such a conclusion based on an analysis of the duality at m 0 = 0 between the CSB and charged PC phenomena carried out in refs. [49][50][51].
Note that our study confirms the previously known fact that chiral asymmetry induces the phenomenon of condensation of charged pions in a dense baryonic medium [67][68][69][70]. In addition, recall that it was previously shown in refs. [30][31][32] that in an electrically neutral medium (without chiral asymmetry) a charged PC phase is impossible. Therefore, in the future it would be of course interesting to consider the phase structure of the two-color and chirally asymmetric NJL model (2.1)-(2.2) taking into account other real parameters characterizing a dense baryon medium in neutron stars. Namely, to take into account a strong external magnetic field, as well as the fact that in the cores of compact stars, dense matter is electrically neutral and in beta equilibrium (for this it is necessary to introduce electrons into consideration) and see their effect on the charged PC phase. Of course, in this case, in the framework of the two-color QCD the results will depend on how the charges of the quarks are selected (see above comments in the text). At N c = 3, some information on the investigations in this direction is already available [76]. So we can conclude that in an electrically neutral medium in beta equilibrium, the chemical potential µ I5 is able to induce a charged PC phase. However, the range of parameters at which this phenomenon can be JHEP06(2020)148 realized is much wider if both µ I5 and µ 5 are nonzero (the details of the investigation have been presented in a recent paper [76]). Although, probably consideration of these conditions is slightly more pertinent in the case N c = 3, it could be useful to generalize this study to chirally asymmetric two-color quark matter as well.
In the same way, it would be interesting to note that we have also some preliminary information that the inclusion of one more chiral chemical potential µ 5 into account does not spoil the dual properties of the model (2.1)-(2.2), but these dualities can be nicely demonstrated in the case of zero µ 5 without considering the most general case. Otherwise in the most general case there are a lot of parameters and there would be a lot of figures that are rather hard to track. Of course, the phase structure becomes more complicated, but µ 5 does not break the dualities between CSB, charged PC and BSF phases of two-colored quark matter. The influence of µ 5 on different phases will be considered in a forthcoming studies. Moreover, we are going to consider the properties of dense baryonic matter with isospin and chiral asymmetries in the framework of a more realistic three-color NJL model withqq and diquark qq channels. where z = x − y, it is possible to obtain from the above formulae

JHEP06(2020)148
If there is an operator function F (Â), whereÂ is a translationally invariant operator, then in the coordinate representation its matrix elements depend on the difference (x − y). Obviously, it is possible to define the Fourier transformations F (A)(p) of its matrix elements, and the following relations are valid (A(p) is the Fourier transformation (A.2) for the matrix element A(x − y)): Finally, suppose thatÂ is an operator in some internal n-dimensional vector space, in addition. Evidently, the same is valid for the Fourier transformation A(p) which is now some n×n matrix. Let λ i (p) are eigenvalues of the n×n matrix A(p), where i = 1, 2, . . . , n. Then In this formula we use the notation tr for the trace of any operator in the internal n-
(C. 4) In this case, using the definition (2.15) of the TDP, it is possible to get from eq. (C.2) the following expression for the TDP (3.13) at ∆ = 0, Ω(M, π 1 , |∆| = 0) = (M − m 0 ) 2 + π 2 1 4H + i d 4 p (2π) 4 ln P − 4 (η)P + 4 (η)P + 4 (η)P − 4 (η) . (C.5) In ref. [66] the phase structure of the NJL model without diquark channel was investigated in the mean-field approximation (or in the leading order of the large-N c expansion) and with the same chemical potential content, i.e. at µ B = 0, µ I = 0, µ I5 = 0. There only the single (qq) channel with coupling constant G was taken into account and quarks are considered as three color. It is possible to show that at 2H = G the TDP (C.5) is actually the doubled TDP obtained in ref. [66]. It can be explained by the following. If one divides the coupling constant G by the number of colors N c , it can be shown that the TDP divided by N c does not depend on the number of colors if one keep G the same for all number of colors. In the present paper neither H nor the model TDP (2.15) is divided by N c . That is why the TDP (C.5) in our case is equal to doubled TDP in three color case obtained in ref. [66] if one substitute G = 2H. Now let us discuss the fit of the model parameters.
The fit in the two color case can be constructed from the fit in three color case by the following N c scaling of physical quantities f π ∼ √ N c and ψ ψ ∼ N c [15]. It can be easily shown that it is equivalent to the scaling of coupling constant G ∼ 1/N c . Together with the above statement that TDP divided by N c does not depend on the number of colors, if G is divided by N c (scaling G ∼ 1/N c ), one can see that the phase structure of NJL model with quark-antiquark channel should not depend on the number of colors. If one divide the coupling constant G by the number of colors N c as it was done, for example, in [66] then the TDP of the system divided by N c does not depend on the number of colors.
Hence, if diquark condensate ∆ is equal to zero in the framework of two-color NJL model, then the phase structures of this model and the ordinary three-color NJL model coincide at 2H = G and the same values of the cutoff parameter Λ. Note that at ∆ = 0 the quantities (products) λ 1 (p)λ 2 (p) and λ 3 (p)λ 4 (p) from the expression (3.13) are just the products P − 4 (η)P + 4 (η) and P + 4 (η)P − 4 (η), respectively. Moreover, the results presented after eq. (C.2) can be obtained with the help of any program of analytical calculations. Let us denote by x 1 , . . . , x 4 the roots of the polynomial P − 4 (x). Then it is clear that the polynomial P + 4 (x) has the same roots but with opposite sign, −x i (i = 1, . . . , 4). So,
The case M = 0, π 1 = 0, ∆ = 0. To obtain the TDP in this case, one can take into account that at m 0 = 0 it is invariant under the dual transformation D 3 (3.15). So, in this case the TDP is equal to the projection F 1 (M ) (3.18), in which we should perform the changes M → π 1 and ν ↔ ν 5 . The result is the expression F 2 (π 1 ) (3.19).
The case M = 0, π 1 = 0, ∆ = 0. Finally, to get the TDP (3.13) in this particular case, one can simply take into account that at m 0 = 0 it is invariant under the dual transformation D 2 (3.15). So, in this case it is no more than the expression (3.18), in which M → |∆| and ν 5 ↔ µ changes should be performed. The result is the expression F 3 (|∆|) (3.20). Note that in all projections F i of the TDP on the condensate axes, the first two terms in square brackets of the expressions (3.18)- (3.20) correspond to the contribution of the eigenvalues λ 1 (p) and λ 2 (p) to the expression (3.13) of the TDP. However, the last two terms in square brackets there correspond to the contribution of λ 3 (p) and λ 4 (p) to it.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.