it is well known that relativistic quantum field theory models with four - fermion ( 4f ) interactions serve as effective theories for low energy considerations of different real phenomena in a variety of physical branches .
for example , the meson spectroscopy , neutron star and heavy - ion collision physics are often investigated in the framework of ( 3 + 1)-dimensional 4f theories @xcite , known as nambu  jona - lasinio ( njl ) models @xcite .
in particular , low dimensional 4f field theories provide a powerful tool for investigations in condensed - matter physics .
indeed , physics of ( quasi)one - dimensional organic peierls insulators ( the best known material of this kind is polyacetylene ) is well described in terms of the ( 1 + 1)-dimensional 4f gross - neveu ( gn ) model @xcite .
the quasirelativistic treatment of electrons in planar systems like high - temperature superconductors or in graphene ( a planar monoatomic layer of carbon atoms ) is also possible in terms of ( 2 + 1)-dimensional gn models @xcite .
notice that there are physical effects , which were observed for the first time just in the framework of njl- and gn - type models .
in particular , the phenomenon of a dynamical generation of a fermion mass is well - known for strong interaction physics since the time , when nambu and jona - lasinio @xcite proved a dynamical breaking of the continuous @xmath2-symmetry on the basis of a generic 4f interaction theory .
later on , this effect served as the basis of a qualitatively successful description of the low - energy meson spectrum of quantum chromodynamics ( qcd ) @xcite .
the effect of dynamical symmetry breaking and generation of a fermion mass is also known in low dimensional , d=1 + 1 and d=2 + 1 , gn theories @xcite , where the four - fermion theory is renormalizable and asymptotic free @xcite in the case of d=1 + 1 , whereas for d=2 + 1 the 4f gn models are perturbatively nonrenormalizable but become renormalizable in the framework of the @xmath3 expansion technique @xcite ( @xmath4 being the number of fermion fields ) .
another example of this kind is the effect of spontaneous chiral symmetry breaking induced by external magnetic or chromomagnetic fields .
this effect was for the first time studied also in terms of the ( 2 + 1)-dimensional gn model @xcite .
note in addition that due to a rather simple structure , different low dimensional gn models provide a good laboratory for a deeper study of dense baryon matter and , in particular , for the consideration of such phenomena as ( color)superconductivity @xcite , charged pion condensation @xcite , etc .
moreover , these theories are very useful in developing new quantum field theoretical techniques like the @xmath3 expansion @xcite , the optimized perturbation theory @xcite , and so on .
it is necessary to point out that there is one more , not yet mentioned above , phenomenon that is still a well - known feature of only some ( 1 + 1)-dimensional 4f theories .
this is the duality correspondence between chiral symmetry breaking and superconductivity @xcite .
( in order to avoid the prohibition of cooper pairing as well as spontaneous breaking of a continuous symmetry in ( 1 + 1)-dimensional models @xcite , the consideration should be performed in the leading order of the @xmath3-technique , i.e. in the large-@xmath4 limit assumption . in this case
quantum fluctuations , which would otherwise destroy a long - range order corresponding to spontaneous symmetry breaking , are suppressed by @xmath3 factors . ) to formulate this phenomenon , let us imagine that there is a microscopic 4f theory ( for d=1 + 1 see , e.g. , @xcite ) which describes a competition between fermion - antifermion ( or chiral ) and difermion ( or superconducting ) channels of interaction .
moreover , we suppose that the ground state of the model is characterized by nonzero fermion number- and chiral charge densities ( the last is an imbalance between densities of left- and right - handed fermions ) , i.e. there are two external parameters , fermion number- @xmath0 and chiral charge- @xmath1 chemical potentials , respectively .
then the duality phenomenon means that the phase structure of this 4f model is symmetric with respect to the following simultaneous transformations @xmath5 , @xmath6 and chiral symmetry breaking @xmath7 superconductivity ( here @xmath8 and @xmath9 are the coupling constants in the chiral and superconducting channels , respectively ) .
thus , there is a correspondence between properties ( phase structure ) of the model , e.g. , at @xmath10 and @xmath11 , etc .
moreover , knowing condensates and other dynamical and thermodynamical quantities of the system , e.g. in the chirally broken phase , one can obtain the corresponding quantities in the superconducting phase of the model , by simply performing there the duality transformation .
it is worth to note that in recent years properties of media with nonzero chiral chemical potential @xmath1 , i.e. chiral media , attracted considerable interest ( see , e.g. , @xcite and references therein ) . in nature , chiral media might be realized in heavy - ion collisions , compact stars , condensed matter systems , etc .
@xcite ( see also the review @xcite ) . in particular , one can expect that in the quark - gluon - plasma phase of qcd a chirality imbalance is produced , i.e. there appears a nonzero chiral charge density @xmath12 .
the combined effect of the chiral imbalance and of the very strong magnetic field , which can be produced at noncentral heavy ion collision , results in the so - called chiral magnetic effect ( see , e.g. , the paper @xcite ) .
it means that there might be induced an electric current along the direction of the magnetic field .
it is important to note that this phenomenon is described effectively in the framework of 4f models with a chiral chemical potential @xmath1 , i.e. in terms of a quantity conjugated to the chiral charge density @xmath12 .
recently , it was established that in planar condensed matter systems there may exists a phenomenon , which can be regarded as an analogue of the qcd chiral magnetic effect @xcite .
this is the pseudo chiral magnetic effect , which can be observed , e.g. , in a distorted graphene sheet under the influence of external in - plane magnetic field .
indeed , due to a mechanical distortion of the lattice structure of graphene , there might appears a nonzero chiral charge density @xmath12 ( chiral imbalance ) as well as a conserved electric current along the in - plane external magnetic field ( for details , see the paper @xcite ) .
notice that , similar to the qcd chiral magnetic effect , one might expect that the pseudo chiral magnetic effect can be effectively described in terms of ( 2 + 1)-dimensional 4f models with nonzero chiral chemical potential @xmath1 ( in addiion to the usual chemical potential @xmath0 ) . finally , it is important to remark that a mechanical distortion of the lattice structure in graphene - like materials ( as well as other external perturbations ) can also lead to opening of different superconducting channels in the system @xcite .
thus , it makes sense to study a competition between chiral symmetry breaking and superconductivity in ( 2 + 1)-dimensional quantum field theories with chiral chemical potential @xmath1 .
in particular , it is important to predict / observe the possible duality relations between these qualitatively different phenomena accompanied by @xmath1 in field theories in spacetime dimensions d@xmath131 + 1 . in this paper
we demonstrate that there exists a dual correspondence ( or dual symmetry ) between the phenomena of the chiral symmetry breaking and superconductivity in the framework of some ( 2 + 1)-dimensional 4f models .
the consideration is performed at zero temperature @xmath14 , however it can be easily generalized to the case @xmath15 .
we hope that our investigations shed some new light on the physical effects in planar condensed matter systems .
the paper is organized as follows . in sec .
ii a ( 2 + 1)-dimensional 4f gn - type model , containing both fermion - antifermion ( or chiral ) and difermion ( or superconducting ) interaction channels and including two kinds of chemical potentials , @xmath16 , is presented . here
we show that under the pauli - grsey transformations of fermi fields there is a dual relationship between the relevant 4f structures .
next , the unrenormalized thermodynamic potential ( tdp ) of the gn - type model is given in the leading order of the large - n expansion . in sec
iii the dual symmetry of the model tdp is established .
it means that it is invariant under the interchange of coupling constants , @xmath16 chemical potentials and chiral and superconducting condensates .
moreover , the renormalization of the tdp is performed .
iv contains a detailed numerical investigation of various phase portraits with particular emphasis on the role of the duality symmetry of the tdp .
some technical details are relegated to appendices . moreover , in appendix [ apd ] we present an alternative ( 2 + 1)-dimensional gn model with dual relationship between other chiral and superconducting channels of 4f interaction .
our investigation is based on a ( 2 + 1)-dimensional gn  type model with massless fermions belonging to a fundamental multiplet of the auxiliary @xmath17 flavor group .
its lagrangian describes the interaction both in the fermion  antifermion and difermion ( or superconducting ) channels : @xmath18\psi_k&+ & \frac { g_1}n\left ( 4f\right ) _
{ ch}+\frac { g_2}n\left ( 4f\right ) _ { sc},\label{1}\end{aligned}\ ] ] where the four - fermion structures @xmath19 and @xmath20 are used , @xmath21 in addition , @xmath0 and @xmath1 in ( 1 ) denote a fermion number chemical potential and a chiral ( axial ) chemical potential , respectively .
@xmath0 is conjugated to a fermion number density @xmath22 , whereas @xmath1 is conjugated to a nonzero density of chiral charge @xmath23 , which represents an imbalance in densities of the right- and left - handed fermions .
as it is noted above , all fermion fields @xmath24 ( @xmath25 ) form a fundamental multiplet of the @xmath17 group . moreover ,
each field @xmath24 is a four - component ( reducible ) dirac spinor ( the symbol @xmath14 denotes the transposition operation ) .
therefore , in eqs .
( 1)-(2 ) the quantities @xmath26 ( @xmath27 ) and @xmath28 are matrices in the four - dimensional spinor space .
moreover , @xmath29 is the charge conjugation matrix .
the algebra of these matrices as well as their particular representations are given in appendix [ apa ] .
clearly , the lagrangian @xmath30 is invariant under transformations from the internal auxiliary @xmath17 group , which is introduced here in order to make it possible to perform all the calculations in the framework of the nonperturbative large-@xmath4 expansion method .
physically more interesting is that the model ( 1 ) is invariant under transformations from the @xmath31 group , where @xmath32 is the fermion number conservation group , @xmath33 ( @xmath34 ) , and @xmath35 is the group of continuous chiral transformations , @xmath36 ( @xmath25 ) .
it means that in the framework of the model ( 1 ) both the particle number density @xmath37 and the density of chiral charge @xmath38 are conserved quantities .
, which anticommutes with @xmath26 ( @xmath39 ) and @xmath28 , one can consider another continuous chiral @xmath40 transformation group of the spinor fields , @xmath41 ( @xmath25 ) . alternatively , there exists a 4f model with another fermion  antifermion and difermion channels of interaction , which , in addition to @xmath32 , is invariant under the continuous chiral @xmath40 group ( see appendix [ apd ] ) .
it is worth mentioning that in the case of graphene reducible four - component spinors just describe the two sublattice ( or pseudospin ) and two valley ( or dirac point ) degrees of freedom of the hexagonal honeycomb lattice of carbon atoms ( see e.g. @xcite and references therein ) . ]    before studying the thermodynamics of the model , we would like to point out that there is a so - called duality correspondence between chiral and superconducting channels of the model ( 1)-(2 ) . to see this
, it is very useful to form an infinite set @xmath42 composed of all lagrangians @xmath43 when the free model parameters @xmath44 and @xmath1 take arbitrary admissible values , i.e. @xmath45 at arbitrary fixed values of coupling constants @xmath46 , @xmath47 and chemical potentials @xmath48 . then , let us perform in ( 1)-(2 ) the so - called pauli - grsey ( pg ) transformation of spinor fields @xcite , @xmath49 taking into account that all spinor fields anticommute with each other , it is easy to see that under the action of the transformations ( [ 0004 ] ) the 4f structures of the lagrangian ( 1 ) are converted into themselves , i.e. @xmath50 and , moreover , each element ( lagrangian ) @xmath43 of the set @xmath51 is transformed into another element of the set @xmath51 according to the following rule @xmath52 i.e. the set @xmath51 is invariant under the field transformations ( [ 0004 ] ) .
owing to the relations ( [ 040 ] ) and ( [ 004 ] ) there is a connection between properties of the model at some fixed free model parameters @xmath53 and properties of the model in the case , when @xmath5 and @xmath6 . in particular ,
if at some fixed @xmath44 and @xmath1 we have the chiral symmetry breaking ( csb ) phase , then at @xmath5 and @xmath6 one can definitely predict the superconducting ( sc ) phase , and vice versa . due to this reason
, we will call the relations ( [ 040 ] ) and ( [ 004 ] ) the duality property of the model ( or duality correspondence between csb and sc ) .
note also that in appendix [ apd ] we present another ( 2 + 1)-dimensional 4f model , in which the duality correspondence between alternative chiral and superconducting channels is realized .
further on , we will study the role and the influence of the duality property of the model ( 1 ) on its phase structure . to this end , we introduce the semi - bosonized version of lagrangian ( [ 1 ] ) that contains only quadratic powers of fermionic fields as well as auxiliary bosonic fields @xmath54 , @xmath55 , @xmath56 and @xmath57 , @xmath58\psi_k   -\frac{n(\sigma^2+\pi^2)}{4g_1 } -\frac n{4g_2}\delta^{*}\delta-   \frac{\delta^{*}}{2}[\psi_k^tc\psi_k ] -\frac{\delta}{2}[\bar\psi_k c\bar\psi_k^t ] .
\label{2}\end{aligned}\ ] ] ( here and in what follows the summations over repeated indices @xmath25 and @xmath59 are implied . ) clearly , the lagrangians ( 1 ) and ( [ 2 ] ) are equivalent , as can be seen by using the euler - lagrange equations of motion for bosonic fields which take the form @xmath60 one can easily see from ( [ 3 ] ) that the neutral fields @xmath61 and @xmath55 are real quantities , i.e. @xmath62 and @xmath63 ( the superscript symbol @xmath64 denotes the hermitian conjugation ) , but the ( charged ) difermion fields @xmath56 and @xmath65 are hermitian conjugated complex quantities , so that @xmath66 and vice versa . moreover , under the chiral @xmath35 group the fields @xmath67 are singlets , but the fields @xmath68 are transformed in the following way @xmath69 clearly , all the fields ( [ 3 ] ) are also singlets with respect to the auxiliary @xmath17 group , since the representations of this group are real . moreover , with respect to the parity transformation @xmath70 , @xmath71 the fields @xmath61 , @xmath72 and @xmath57 are even quantities , i.e. scalars , but @xmath55 is a pseudoscalar .
if the difermion field @xmath56 has a nonzero ground state expectation value , i.e.   @xmath73 , then the abelian fermion number conservation @xmath32 symmetry of the model is spontaneously broken down and the superconducting phase is realized in the model . however , if @xmath74 then the continuous @xmath35 chiral symmetry of the model is spontaneously broken .
let us now study the phase structure of the four - fermion model ( 1 ) starting from the equivalent semi - bosonized lagrangian ( [ 2 ] ) . in the leading order of the large-@xmath4 approximation ,
the effective action @xmath75 of the considered model is expressed by means of the path integral over fermion fields : @xmath76[d\psi_l]\exp\bigl(i\int { \cal    l}\,d^3 x\bigr),\ ] ] leading to @xmath77 + \widetilde { \cal s}_{\rm { eff}}.   \label{5}\end{aligned}\ ] ] the term @xmath78 in ( [ 5 ] ) is the fermion contribution to the effective action and is given by : @xmath79[d\psi_l]\exp\bigl\{i\int\big [ \bar\psi_k(\gamma^\nu i\partial_\nu + \mu\gamma^0+\mu_5\gamma^0\gamma^5 -\sigma - i\gamma^5\pi ) \psi_k -   \frac{\delta^{*}}{2}(\psi_k^t c\psi_k ) -\frac{\delta}{2}(\bar\psi_kc \bar\psi_k^t)\big ] d^3 x\bigr\}. \label{6}\ ] ] the ground state expectation values @xmath80 , @xmath81 , etc .
of the composite bosonic fields are determined by the saddle point equations , @xmath82 for simplicity , throughout the paper we suppose that the above mentioned ground state expectation values do not depend on spacetime coordinates , i.e. @xmath83 where @xmath84 are constant quantities .
in fact , they are coordinates of the global minimum point of the thermodynamic potential ( tdp ) @xmath85 . in the leading order of the large-@xmath4 expansion and using ( 13 ) it is defined by the following expression : @xmath86 the tdp ( [ 08 ] ) is invariant with respect to chiral @xmath35 symmetry group .
so , as it is clear from ( [ 003 ] ) , it depends on the quantities @xmath87 and @xmath88 through the combination @xmath89 .
moreover , without loss of generality , one can suppose that @xmath90 .
thus , to find the other ground state expectation values @xmath80 etc .
, it is enough to study the global minimum point of the tdp @xmath91 , @xmath92 taking into account the relations ( [ 5 ] ) , ( [ 6 ] ) and ( [ 08 ] ) , we have from ( [ 09 ] ) @xmath93[d\psi_l]\exp\big ( i\int d^3 x\big [ \bar\psi_k d\psi_k\right.\nonumber\\ & & \left.- \frac{\delta^*}{2}(\psi_k^tc\psi_k ) -\frac{\delta}{2}(\bar\psi_k c\bar\psi_k^t)\big ] \big ) \right ) , \label{9}\end{aligned}\ ] ] where @xmath94 .
to proceed further , let us point out again that without loss of generality the quantities @xmath95 might be considered as real ones .
might be eliminated by an appropriate transformation of fermion fields in the path integral ( [ 9 ] ) .
] so , in what follows we will suppose that @xmath96 , where @xmath97 now is already a real quantity .
the path integration in the expression ( [ 9 ] ) is evaluated in appendix [ apb ] , we consider for simplicity the case @xmath98 , however the procedure is easily generalized to the case with @xmath99 . ]
so we have for the tdp ( [ 9 ] ) the following expression @xmath100 , \label{11}\end{aligned}\ ] ] where @xmath101 are presented in ( [ b12 ] ) and superscription  un  denotes the unrenormalized quantity .
note that the tdp ( [ 11 ] ) describes thermodynamics of the model at zero temperature .
taking into account the expressions ( [ b12 ] ) , the tdp ( [ 11 ] ) can be presented in the form @xmath102 where @xmath103 and @xmath104 it is clear from ( [ 013 ] ) that the tdp ( [ 012 ] ) is an even function of each of the quantities @xmath0 , @xmath1 , @xmath87 , and @xmath97 , i.e. without loss of generality we can consider in the following only @xmath105 , @xmath106 , @xmath107 , and @xmath108 values of these quantities . moreover ,
the tdp ( [ 012 ] ) is invariant with respect to the so - called duality transformation @xcite , @xmath109 notice that this invariance of the tdp is a consequence of the rule ( [ 004 ] ) , according to which the 4f lagrangian ( 1 ) is modified under the pauli - gursey transformation ( [ 0004 ] ) of the spinor field .    in powers of @xmath97 the fourth - degree polynomial @xmath110 has the following form @xmath111 expanding the right - hand side of ( [ 17 ] ) in powers of @xmath87
, one can obtain an equivalent alternative expression for this polynomial .
namely , @xmath112 note also that according to the general theorem of algebra , the polynomial @xmath110 can be presented in the form @xmath113 where @xmath114 , @xmath115 , @xmath116 and @xmath117 are the roots of this polynomial .
the fourth - order polynomial with similar coefficients @xmath118 ( [ 013 ] ) was studied in our paper @xcite , where it was shown that all its roots @xmath119 ( @xmath120 ) are real quantities ( see appendix b in @xcite ) .
the roots @xmath119 are the energies of quasiparticle or quasiantiparticle excitations of the system . in particular
, it follows from ( [ 17 ] ) that at @xmath121 @xmath122 whereas it is clear from ( [ 18 ] ) that at @xmath123 @xmath124 taking into account in ( [ 012 ] ) the relation ( [ 170 ] ) as well as the formula @xmath125 ( obtained rigorously , e.g. , in appendix b of @xcite and being true up to an infinite term independent of the real quantity @xmath126 ) , it is possible to integrate there over @xmath127 .
so the unrenormalized tdp ( [ 012 ] ) can be presented in the following form , @xmath128      first of all , let us obtain a finite , i.e. renormalized , expression for the tdp ( [ 28 ] ) at @xmath130 and @xmath131 , i.e. in vacuum .
since in this case a thermodynamic potential is usually called effective potential , we use for it the notation @xmath132 .
it follows from ( [ 012 ] ) and ( [ 013 ] ) that at @xmath130 and @xmath131 @xmath132 looks like @xmath133\nonumber\\ & \equiv&\frac{m^2}{4g_1}+\frac{\delta^2}{4g_2}-\int\frac{d^2p}{(2\pi)^2}\big ( \sqrt{|\vec p|^2+(m+\delta)^2}+\sqrt{|\vec p|^2+(m-\delta)^2}\big ) , \label{25}\end{aligned}\ ] ] where @xmath134 .
( to obtain the second line of this formula , we have integrated there over @xmath127 according to the relation ( [ int ] ) . )
it is evident that the effective potential ( [ 25 ] ) is an ultraviolet divergent quantity .
so , we need to renormalize it .
this procedure consists of two steps : ( i ) first of all we need to regularize the divergent integral in ( [ 25 ] ) , i.e. we suppose there that @xmath135 , @xmath136 .
( ii ) second , we must suppose also that the bare coupling constants @xmath8 and @xmath9 depend on the cutoff parameter @xmath137 in such a way that in the limit @xmath138 one obtains a finite expression for the effective potential .    before performing the steps ( i ) and ( ii ) of the renormalization procedure
, it is useful to take into account the following asymptotic expansion at @xmath139 @xmath140 then , after construction of the regularized effective potential @xmath141 ( see below ) , we use there the asymptotic expansion ( [ 19 ] ) and integrate it over @xmath142 and @xmath143 term - by - term .
the result reads @xmath144+\delta^2\left
[ \frac 1{4g_2}-\frac{2\lambda\ln(1+\sqrt{2})}{\pi^2}\right ] -\frac{2\lambda^3(\sqrt{2}+\ln(1+\sqrt{2}))}{3\pi^2}+{\cal o}(\lambda^0 ) , \label{12}\end{aligned}\ ] ] where @xmath145 denotes an expression which is finite in the limit @xmath146 .
now , we should suppose that the bare coupling constants @xmath8 and @xmath9 depend on the cutoff parameter @xmath137 in such a way that in the limit @xmath138 one obtains a finite expressions in the square brackets of ( [ 12 ] ) . clearly , to fulfil this requirement it is sufficient to require that @xmath147 where @xmath148 are finite and @xmath137-independent model parameters with dimensionality of inverse mass . since bare couplings @xmath8 and @xmath9 do not depend on a normalization point ,
the same property is also valid for @xmath148 .
hence , taking into account in ( [ 12 ] ) the relations ( [ 13 ] ) and ignoring there an infinite @xmath87- and @xmath97-independent constant , one obtains the following _ renormalized _ , i.e. finite , expression @xmath149 for the effective potential , @xmath150 it should also be mentioned that the effective potential ( [ 14 ] ) is a renormalization group invariant quantity .
the fact that it is possible to renormalize the effective potential of the initial model ( 1 ) in the leading order of the large @xmath4-expansion is the reflection of a more general property of ( 2 + 1)-dimensional theories with four - fermion interactions .
indeed , it is well known that in the framework of the  naive  perturbation theory ( over coupling constants ) these models are not renormalizable . however , as it was proved in @xcite , in the framework of nonperturbative large @xmath4-technique these models are renormalizable in each order of the @xmath3-expansion .    the @xmath145 term in ( [ 12 ] ) can be calculated explicitly , so we have for the renormalized effective potential @xmath149 ( [ 14 ] ) the following expression @xmath151 the global minimum point @xmath152 , where @xmath153 and @xmath154 , of this function was already investigated in @xcite , although in the framework of another ( 2 + 1)-dimensional gn model .
so , we present at once the phase structure of the initial model ( 1 ) at @xmath130 and @xmath131 ( see fig .
1 ) .    in fig .
1 the phase portrait of the model is depicted depending on the values of the free model parameters @xmath155 and @xmath156 . there
the plane @xmath157 is divided into several areas . in each area one of the phases
i , ii or iii is implemented . in the phase
i , i.e. at @xmath158 and @xmath159 , the global minimum of the effective potential @xmath160 is arranged at the origin .
so in this case we have @xmath161 and @xmath162 . as a result ,
in the phase i both continuous symmetries , chiral @xmath35 and electromagnetic @xmath32 , remain intact and fermions are massless . due to this reason
the phase i is called symmetric . in the phase ii , which is allowed only for @xmath163 , at the global minimum point
the relations @xmath164 and @xmath162 are valid .
so in this phase chiral @xmath35 symmetry is spontaneously broken down and fermions acquire dynamically the mass @xmath165 . finally , in the superconducting phase iii , where @xmath166 , we have spontaneous breaking of the @xmath32 symmetry .
this phase corresponds to the following values of the gaps : @xmath161 , @xmath167 .
note also that if @xmath168 and , in addition , @xmath169 ( it is just the line _
l _ in fig .
1 ) , then the effective potential ( [ 15 ] ) has two equivalent global minima . the first one , the point @xmath170 , corresponds to a phase with chiral symmetry breaking . the second one , i.e. the point @xmath171 , corresponds to superconductivity .
clearly , if the cutoff parameter @xmath137 is fixed , then the phase structure of the model can be described in terms of bare coupling constants @xmath172 instead of finite quantities @xmath173 .
indeed , let us first introduce a critical value of the couplings , @xmath174 . then , as it follows from fig . 1 and ( [ 13 ] ) , at @xmath175 and @xmath176 the symmetric phase i of the model is located . if @xmath177 , @xmath176 ( @xmath175 , @xmath178 ) , then the chiral symmetry broken phase ii ( the superconducting phase iii ) is realized .
finally , let us suppose that both @xmath177 and @xmath178 . in this case at @xmath11 ( @xmath10 )
we have again the chiral symmetry broken phase ii ( the superconducting phase iii ) .
now , few comments about the nature of the superconductivity ( sc ) and chiral symmetry breaking ( csb ) are in order . in the framework of the model ( 1 )
there are two well - known mechanisms for appearing of these phenomena , ( i ) dynamical symmetry breaking , which occurs at strong couplings , and ( ii ) cooper instability of fermi surface , which takes place , in contrast to the case ( i ) , at weak couplings .
if @xmath130 and @xmath131 , then the initial model ( 1 ) describes effectively the undoped regime of high temperature cuprate superconductors , etc ( see , e.g. , the papers @xcite , where the phase structure , similar to the phase diagram of our fig .
1 , was obtained ) . in this case fermi surface
is absent , so both csb and sc phases appear dynamically in the system at a _ rather strong attraction _ in the fermion - antifermion or difermion channels , i.e. at @xmath177 and @xmath11 or at @xmath178 and @xmath179 , respectively .
hence , at zero values of chemical potentials only the mechanism ( i ) for generating both sc and csb is realized in our model .
is well known since the papers by nambu and jona - lasinio [ 3 ] .
it is also valid in the framework of other ( 3 + 1)-dimensional 4f models with chiral and color superconducting channels of interaction at @xmath130 @xcite . in @xcite
the bare coupling constants @xmath8 and @xmath9 are free model parameters .
moreover , strong constraints were obtained there on @xmath8 and @xmath9 , at which color sc or csb is generated dynamically in the models . ]    the second way , i.e. the mechanism ( ii ) , to break spontaneously a symmetry can be realized only at nonzero chemical potentials .
for example , if the doping or impurities are present in the system , then one must introduce a finite chemical potential @xmath0 . in this case
the fermi surface can be created . as a result ,
if there is an _ arbitrary small attraction _ in the sc channel , i.e. even at @xmath176 , there appears the so - called cooper instability of the fermi surface , which destroys the normal ground state in favor of superconductivity .
( below , when discussing the phase portrait of fig . 5 ( see sect [ ivb ] ) , we will see that at @xmath131 and @xmath180 both ( i ) and ( ii ) mechanisms take part in forming the phase structure . ) if @xmath130 but @xmath181 , then csb can also be generated in the model due to the mechanism ( ii ) at arbitrary small values of @xmath8 ( see in @xcite ) .
this fact supports the duality correspondence between sc and csb in the model ( 1 )
.      +      to renormalize the tdp ( [ 28 ] ) in the most general case , it is necessary to imagine how the integrand of ( [ 28 ] ) , i.e. the quantity @xmath182 , behaves at asymptotically high values of @xmath183 .
the properties of the roots @xmath119 ( @xmath120 ) of the polynomial @xmath110 appearing in ( [ 012 ] ) can be obtained both numerically and analytically using the methods of appendix b of @xcite .
so it is easy to show that at @xmath139 @xmath184 since the asymptotic expansions ( [ 016 ] ) and ( [ 19 ] ) coincide , one can subtract from the integrand of ( [ 28 ] ) the expression @xmath185 , thereby obtaining a convergent integral .
taking into account this circumstance , we have identically @xmath186 where @xmath187 is the effective potential ( [ 25 ] ) of the model in vacuum , i.e. at @xmath130 and @xmath131 , and @xmath188.\label{018}\end{aligned}\ ] ] it is clear that the term @xmath189 in ( [ 017 ] ) is a convergent integral and all ultraviolet divergences of the tdp ( [ 28])-([017 ] ) are located in the term @xmath187 .
hence , we can renormalize only the first term of ( [ 017 ] ) in the way of the previous subsection and obtain finally : @xmath190 where @xmath149 and @xmath189 are presented in ( [ 15 ] ) and ( [ 018 ] ) , correspondingly .
it was already mentioned above that using the method presented in appendix b of the paper @xcite , the roots @xmath119 can be calculated numerically ( at fixed values of @xmath191 , etc ) . as a result , it is also possible to study numerically the whole tdp ( [ 20 ] ) . carrying out this procedure ,
we have seen that the global minimum of the tdp @xmath192 is always at the point of the form @xmath193 or @xmath194 .
so , to get a more detailed information about the phase structure of the model , it is sufficient to investigate the reductions of the tdp ( [ 20 ] ) on the @xmath87- and @xmath97 axes on the @xmath87- and @xmath97 axes is known ( see ( [ 26 ] ) and ( [ 27 ] ) ) .
] , i.e. to consider the functions @xmath195\nonumber\\ & = & v_1(m)-\int\frac{d^2p}{(2\pi)^2 } \sum_{\eta = \pm}\big [ { \cal e}_{\eta}-\sqrt{|\vec p|^2+m^2}+(\mu-{\cal e}_{\eta})\theta(\mu-{\cal e}_{\eta})\big ] , \label{31}\\ f_2(\delta)&\equiv&\omega^{ren}(m=0,\delta)= v_2(\delta)-\int\frac{d^2p}{(2\pi)^2 } \sum_{\eta = \pm}\big [ \sqrt{\delta^2+(|\vec p|+\eta\mu)^2}-\sqrt{|\vec p|^2+m^2}\nonumber\\&+&\left ( \mu_5-\sqrt{\delta^2+(|\vec p|+\eta\mu)^2}\right ) \theta\left ( \mu_5-\sqrt{\delta^2+(|\vec p|+\eta\mu)^2}\right ) \big ] , \label{32}\end{aligned}\ ] ] respectively , where @xmath196 and @xmath197 is the reduction of the vacuum effective potential ( [ 15 ] ) on the @xmath87 ( in this case @xmath198 , @xmath199 ) or @xmath97 ( in this case @xmath200 , @xmath201 ) axes .
moreover , to obtain the second line in ( [ 31 ] ) we use the evident relations , @xmath202 and @xmath203 . after tedious but straightforward calculations , the tdp ( [ 31 ] ) can be presented in the following form ( see appendix [ apc ] ) : @xmath204\nonumber\\ & & -\frac{\theta\left ( \sqrt{m^2+\mu_5 ^ 2}-\mu\right ) } { 2\pi}\left [ \mu_5 ^ 2\sqrt{\mu_5 ^ 2+m^2}+\mu_5 m^2\ln\left ( \frac{\mu_5+\sqrt{\mu_5 ^ 2+m^2}}{m}\right ) \right ] \nonumber\\ & & -\frac{\theta\left ( \mu - m\right ) \theta\left ( \sqrt{m^2+\mu_5 ^ 2}-\mu\right ) } { 2\pi}\left [ \mu_5\mu\sqrt{\mu^2-m^2}-\mu_5 m^2\ln\left ( \frac{\mu+\sqrt{\mu^2-m^2}}{m}\right ) \right ] \label{c18}\end{aligned}\ ] ] ( at @xmath131 this expression coincides with the corresponding tdp ( 22 ) from the paper @xcite ) . moreover , comparing the expressions ( [ 31 ] ) and ( [ 32 ] ) , we see that the quantity @xmath205 can be easily obtained from the tdp @xmath206 by the substitutions @xmath207 , @xmath208 and @xmath6 , i.e. @xmath209 clearly , the connection ( [ 35 ] ) between the reductions of the tdp @xmath192 on the @xmath87- and @xmath97 axes is the consequence of the duality property ( [ 16 ] ) of the model .
next , in order to find the phase structure of the model , we will determine the global minimum points of the tdps @xmath210 ( [ c18 ] ) and @xmath211 ( [ 35 ] ) and then compare the minimum values of these functions vs external parameters @xmath212 .
suppose now that at some fixed particular values of the model parameters , i.e. at @xmath213 , @xmath214 and @xmath215 , @xmath216 , the global minimum of the tdp ( [ 20 ] ) lies at the point , e.g. , @xmath217 .
it means that for such fixed values of the parameters the chiral symmetry breaking ( csb ) phase is realized .
then it follows from the duality invariance of the unrenormalized tdp with respect to the transformation @xmath218 ( [ 16 ] ) that the permutation of the coupling constant and chemical potential values for the renormalized tdp ( [ 20 ] ) means @xmath219 , @xmath220 , @xmath221 , under which this tdp is invariant . ]
( i.e. at @xmath222 , @xmath223 and @xmath224 , @xmath225 ) moves the global minimum of the tdp @xmath226 to the point @xmath227 , and the superconducting ( sc ) phase is originated ( and vice versa ) .
this is the so - called duality correspondence between csb and sc phases in the framework of the model under consideration .
hence , the knowledge of a phase of the model ( 1 ) at some fixed values of external free model parameters @xmath228 is sufficient to understand what phase is realized at rearranged values of external parameters , @xmath229 .
moreover , we would like to emphasize once again that there exists an equality of the order parameters ( condensates ) , which characterizes both the initial phase and the phase corresponding to rearranged external parameters .
in other words the chiral condensate of the csb phase at fixed @xmath228 is equal to the sc condensate of the phase at @xmath229 ( and vice versa ) .
the duality transformation ( [ 16 ] ) of the tdp can also be applied to the arbitrary phase portrait of the model ( see below ) . in particular , it is clear that if we have a most general phase portrait , i.e. the correspondence between any point @xmath230 of the four - dimensional space of external parameters and possible model phases ( csb , sc and symmetric phase ) , then under the duality transformation ( @xmath231 , @xmath6 , csb@xmath7sc ) this phase portrait is mapped to itself , i.e. the most general phase portrait is selfdual . in practice , usually , there are constraints on the model parameters . as a result ,
if the constraint is dually ( non)invariant , then the phase portrait is also a dually ( non)invariant .
below , we will use the dual symmetry of the tdp in order to explain and construct the phase structure of the model in different particular cases .
+      first of all , taking into account the duality correspondence between csb and sc , let us determine its characteristic features and then find ( analytically and numerically ) the phase structure of the model ( 1 ) in two particular cases , ( i ) @xmath232 and ( ii ) @xmath233 .
it is evident that both of these constraints , ( i ) and ( ii ) , are invariant with respect to the duality transformation @xmath231 , @xmath6 .
hence , the corresponding phase portraits should be selfdual ( a more detailed explanation of this fact is presented in the next subsection ) . in the case ( i )
we want to study the @xmath157 phase portrait , whereas in the case ( ii )  the @xmath234 phase portrait .
it follows from the duality correspondence that in the case ( i ) the region of the csb phase is a mirror image of the sc phase with respect to the line @xmath233 of the @xmath157 plane .
however , in the case ( ii ) the regions of these phases are arranged mirror symmetrically with respect to the line @xmath232 of the plane @xmath234 .    of course , the specific form of the csb and sc areas depends on the values of the external parameters .
for example , at @xmath235 the phase structure is presented in fig
( note that in order to obtain this phase portrait it is sufficient to investigate the behavior of the global minimum point of the tdp @xmath149 ( [ 15 ] ) vs the coupling constants @xmath155 and @xmath156 . ) however , at fixed @xmath236 the phase structure is presented in fig .
2 , where one can see the mirror symmetrical arrangement of the phases ii and iii with respect to the line @xmath233 .    in figs 3 , 4
the @xmath234-phase portraits of the model ( 1 ) are presented at fixed @xmath237 and fixed @xmath238 , respectively .
the csb and sc phases in these figures are arranged , as defined above , symmetrically with respect to the line @xmath232 of the @xmath234 plane . if @xmath239 , then at the points of the boundary of the phase i in fig .
3 , the second order phase transition takes place . in this case it is possible to give exact analytical expressions for the boundaries of the region i in fig .
indeed , the global minimum point @xmath165 of the tdp ( [ c18 ] ) is defined by the equation @xmath240 . inside the region ii of fig .
3 we have @xmath241 and @xmath242 .
however , on the boundary between phases i and ii one can see that @xmath161 and , moreover , that @xmath243 . it is the last relation that defines the boundary equation between i and ii phases .
this equation looks like @xmath244 in a similar way it is possible to find the equation for the boundary between i and iii phases of fig . 3 : @xmath245 in contrast , on the boundary of the phase i of fig .
4 , a phase transition of the first order takes place . since on the boundary between ii and iii phases of this figure
we also have a first order phase transitions , it is clear that the point t of the phase portrait fig .
4 is a so - called triple point , i.e. the point where three different phases , i , ii and iii , coexist .
the triple point t of fig . 4 corresponds to the chemical potential values @xmath246 .        in these two particular cases the restrictions on the external model parameters are already not dually invariant , so the corresponding @xmath157-phase portraits are not selfdual .
however , as we shall see , each of these two phase portraits is a dual mapping of another one .    the @xmath157-phase structure of the model at arbitrary fixed nonzero value of @xmath0 ( of @xmath1 ) and at @xmath131 ( at @xmath130 )
is presented in fig . 5 ( in fig . 6 )
let us first compare fig . 1 and fig
it is easy to see that if @xmath158 , @xmath159 and @xmath235 , then the system is in the symmetric phase i ( see fig .
1 ) . however , an arbitrary small nonzero value of the chemical potential @xmath0 induces ( at @xmath131 ) in this case the superconducting phase iii ( see fig .
5 as well as fig . 3 for the particular case @xmath248 ) . note that there is a very simple explanation of this fact , which is based on the symmetry breaking mechanism ( ii ) ( see the end of the section iii a ) , i.e. on the cooper instability of a fermi surface .
indeed , at @xmath158 , @xmath159 and @xmath180 we have a nonzero particle density ( see , e.g. , in @xcite ) , so there is a fermi sea of particles with energies less or equal to @xmath0 ( fermi surface ) .
evidently , in this case there is no energy cost for creating a pair of particles with opposite momenta just over the fermi surface .
then , due to an arbitrary weak attraction between these particles ( @xmath249 ) , the cooper pair is formed and @xmath32 symmetry is spontaneously broken , as a result of bose
einstein condensation of cooper pairs .
since in the energy spectrum of fermions the gap @xmath250 appears , rather small external forces are not able to destroy the superconducting condensate and it is a stable one .    if @xmath158 , @xmath166 , then we have a rather strong interaction in the sc channel , i.e. @xmath178 and @xmath175 . as a result , in this case both at @xmath130 and @xmath180 the dynamical symmetry breaking mechnism ( i ) ( see the end of the section iii a ) plays a decisive role in the forming of sc phenomenon .
now suppose that @xmath251 and @xmath252 , i.e. @xmath177 .
( in particular , the point @xmath157 might belong to the phase ii of fig .
1 , where chiral symmetry is spontaneously broken at @xmath130 , @xmath131 . )
then it follows from fig .
5 that at @xmath131 there is a sufficiently high critical value @xmath253 of the chemical potential @xmath0 , such that at @xmath254 the point @xmath255 will certainly be in one of the regions iii of fig .
hence , the growth ( at @xmath131 ) of the chemical potential @xmath0 leads to the appearance of the sc phase in the model ( 1 ) at arbitrary @xmath256 and @xmath252 .
, @xmath131 without duality correspondence . ] analysing in the same manner fig . 6 , we see that in the opposite case when @xmath130 but @xmath257 , the growth ( at @xmath130 ) of the chiral chemical potential @xmath1 leads to the appearance of the csb phase in the model ( 1 ) at arbitrary fixed @xmath256 and @xmath252 .    indeed , these two properties of the model ( 1 ) , i.e. ( a ) the appearance of the sc at @xmath131 and at growing values of @xmath0 as well as ( b ) the appearance of the csb at @xmath130 and at growing values of @xmath1 ( in these two cases it is supposed that @xmath256 and @xmath252 are fixed ) , are dually connected ( or dually conjugated ) . to understand this fact , it is necessary to emphasize once again that ( a ) follows from the phase portrait of fig . 5 , whereas ( b ) is the consequence of the phase portrait presented in fig .
6 . now we will demonstrate that under the dual transformation @xmath218 ( @xmath231 , @xmath6 , csb@xmath7sc ) the phase portrait of fig .
5 is transformed into the phase portrait of fig .
6 ( and vice versa ) .
note that the application of @xmath218 to fig .
5 can be divided into three more simple steps .
( i ) first , under permutation @xmath258 in ( [ 16 ] ) we have renaming of the phases , i.e. ii@xmath7iii .
( for example , in this case the global minimum point of the csb phase , i.e. the point @xmath259 , is transformed into the point @xmath260 and , as a result , the csb phase is transformed into the sc phase . )
( ii ) second , performing the @xmath6 and @xmath231 transformations in fig . 5
, we rename the coordinate axes of the figure ( after the duality transformation we have along the vertical axis the quantity @xmath261 and along the horizontal one the quantity @xmath262 ) and change its caption or external conditions , at which the phase portrait is obtained .
( in our case the caption * @xmath131 with @xmath0 being an arbitrary fixed quantity * for fig .
5 is replaced by * @xmath130 with @xmath1 being an arbitrary fixed quantity * , i.e. we obtain the caption for fig .
( iii ) finally , it is necessary to direct the axis corresponding to the quantity @xmath261 ( to the quantity @xmath262 ) horizontally ( vertically ) . as a result of this dual transformation of the fig .
5 we obtain just fig . 6 .      +    in a similar way it is possible to apply the duality transformation @xmath218 to each of the phase portraits presented in figs .
14 . as a result
, we see that these phase structures are selfdual . in the next subsubsections
we obtain the phase structure of the model in other more general nonselfdual cases , i.e. the @xmath157-phase portraits at arbitrary fixed values of chemical potentials and the @xmath234-phase portraits at arbitrary fixed values of the coupling constants .
-phase portrait of the model at arbitrary fixed values of @xmath0 and for different values of the chiral chemical potential @xmath1 .
( a ) the case @xmath263 .
( b ) the case @xmath264 .
( c ) the case @xmath265 .
( d ) the case @xmath266 .
we use the same designations of the phases as in fig . 1 . ,
title="fig:",scaledwidth=45.0% ] -phase portrait of the model at arbitrary fixed values of @xmath0 and for different values of the chiral chemical potential @xmath1 .
( a ) the case @xmath263 .
( b ) the case @xmath264 .
( c ) the case @xmath265 .
( d ) the case @xmath266 .
we use the same designations of the phases as in fig . 1 . , title="fig:",scaledwidth=45.0% ] + -phase portrait of the model at arbitrary fixed values of @xmath0 and for different values of the chiral chemical potential @xmath1 .
( a ) the case @xmath263 .
( b ) the case @xmath264 .
( c ) the case @xmath265 .
( d ) the case @xmath266 .
we use the same designations of the phases as in fig . 1 . , title="fig:",scaledwidth=45.0% ] -phase portrait of the model at arbitrary fixed values of @xmath0 and for different values of the chiral chemical potential @xmath1 .
( a ) the case @xmath263 .
( b ) the case @xmath264 .
( c ) the case @xmath265 .
( d ) the case @xmath266 .
we use the same designations of the phases as in fig . 1 . , title="fig:",scaledwidth=45.0% ]
+      let us now fix the fermion number chemical potential @xmath247 and consider how the @xmath157-phase portrait of the model is evolved vs the chiral chemical potential @xmath1 .
this phase portrait at @xmath131 is presented in fig .
5 . in figs . 7 we have drawn several @xmath157-phase portraits at different relations between @xmath1 and @xmath0 , ( a ) @xmath263 , ( b ) @xmath264 , ( c ) @xmath265 , ( d ) @xmath266 . analyzing figs . 5 and 7
, we see the following phase evolution of the system vs @xmath1 at fixed @xmath247 . at @xmath131
the sc phase iii fills almost the whole @xmath255 plane ( see fig .
it means that at arbitrary fixed @xmath256 and @xmath252 we have superconductivity in the system at sufficiently high values of @xmath0 ( see section [ ivb ] ) . if @xmath1 begins to grow , then for a lot of @xmath157 points there might occur a restoration of the initial symmetry ( the symmetrical phase i arises ) in the model
however , the further growth of @xmath1 leads to the final appearance of the csb phase .
it is the general property of the chiral chemical potential @xmath1 , which promotes the chiral symmetry breaking in the system .    since the above mentioned constraints ( a), .. ,(d ) at which figs .
7 were drawn are not invariant under the duality transformation @xmath267 ( [ 16 ] ) , we would like to note that the phase portraits of fig . 7
are not selfdual . in this case
the possible duality transformation ( see the previous subsubsection [ ivb ] ) of figs .
7 might result in the set of four phase portraits , which are a good illustration of the evolution of the model phase structure vs the particle number chemical potential @xmath0 at arbitrary fixed @xmath181 .
it is straightforwardly clear from the dual mapping of figs .
7 that at arbitrary fixed values of @xmath257 and @xmath256 , @xmath252 we have in the system the appearance of superconductivity at sufficiently high @xmath0 .
this property is dually conjugated to the appearing of csb phase at sufficiently high values of @xmath1 , which is a consequence of the phase portraits fig . 7
( see the above consideration ) .
we have already considered the @xmath234-phase structure of the model in the case of dually symmetrical constraints on the coupling constants , i.e. when @xmath233 ( see figs . 3 , 4 ) . in the present section , special attention
is paid to a more general case when @xmath268 .
here we construct @xmath234-phase portraits of the model for qualitatively different representative relations between @xmath155 and @xmath156 .
namely , we fix in each of the regions of fig . 1 some representative @xmath157 point and draw a corresponding @xmath234-phase diagram .
( note that due to a duality symmetry of the model , it is sufficient to reduce the number of the representative @xmath157 points to three . )    in fig .
8 the @xmath234-phase portrait corresponds to arbitrary fixed @xmath158 and @xmath269 ( it is clear that at @xmath235 the point @xmath270 lies in the symmetrical region i of fig .
1 ) . on the phase boundaries of this figure
there are phase transitions of the second order . in fig .
9 the @xmath234-phase portrait corresponds to arbitrary fixed @xmath163 and @xmath271 ( in this case at @xmath235 the point @xmath272 lies in the csb region ii of fig .
1 ) . finally , in fig . 10 the @xmath234-phase portrait corresponds to arbitrary fixed @xmath163 and @xmath273 ( in this case at @xmath235 the point @xmath274 lies in the sc region iii of fig .
1 ) . to obtain the @xmath234-phase portraits of the model when the point @xmath157 is fixed in the two remaining areas of fig .
1 , it is enough to perform the dual mapping of the figs 9 and 10 .    in figs .
9 , 10 all phase boundaries are the curves of the first order phase transitions , so the point t on these figures is a triple point , where three phases coexist .
the main conclusion from figs 8 - 10 is that the growth of @xmath1 ( the growth of @xmath0 ) induces the chiral symmetry breaking ( induces the superconductivity ) .
+    -phase portrait of the model at arbitrary fixed @xmath163 and @xmath273 .
the notations i , ii and iii are the same as in fig .
the letter t denotes a triple point , @xmath275 , @xmath276.,title="fig:",scaledwidth=45.0% ] +
in this paper , the duality correspondence phenomenon between csb and sc , which was found earlier in ( 1 + 1)-dimensional 4f theories @xcite , is demonstrated to take place also in the framework of the ( 2 + 1)-dimensional 4f model ( 1 ) .
( an alternative 4f model with other dually conjugated chiral and superconducting channels of interaction is considered in appendix [ apd ] ) .
it contains both fermion - antifermion ( with @xmath8 a bare coupling constant ) and difermion ( with @xmath9 a bare coupling constant ) interaction channels and two kinds of chemical potentials , @xmath16 .
the duality property of the model means that its thermodynamic potential @xmath277 ( [ 11 ] ) is invariant under the duality transformation ( [ 16 ] ) . as a result
, we have established the duality correspondence between csb and sc phases of the model ( 1 ) ( as well as of the alternative model of appendix [ apd ] ) .
let us suppose , e.g. , that for some fixed set @xmath230 of external parameters ( the coupling constants @xmath155 and @xmath156 are connected with bare couplings @xmath278 by the relation ( [ 13 ] ) ) the chiral symmetry breaking phase is realized in the model .
then for rearranged values of external parameters , @xmath231 , @xmath6 , ( or for dually conjugated set @xmath279 ) we have the so - called dually conjugated superconducting phase ( and vice versa ) .
moreover , it must be emphasized that the chiral condensate of the csb phase , realized for the set @xmath230 , is equal to the superconducting condensate of the dually conjugated sc phase , corresponding to the set of external parameters with @xmath231 , @xmath6 ( and vice versa ) . in this way , it is sufficient to have the information about the ground state of the initial phase , which is realized for the set @xmath230 , in order to determine the properties of the ground state of the dually conjugated phase , corresponding to the rearranged external parameter set @xmath279 .    in order to study the role and influence of the duality property on the phase structure of the model , for comparison and illustrations , we have demonstrated a variety of phase portraits in the ( @xmath16)- and @xmath157 planes . in particular , if the constraint , under which we obtain a phase portrait , is dually invariant , i.e. at @xmath232 or at @xmath233 , then we have a selfdual phase diagram ( see , e.g. , figs 14 ) .
the selfdual phase portrait is mapped into itself under the duality transformation ( it is introduced in sec .
[ ivb ] ) , and its most characteristic feature is the mirror symmetrical arrangement of the csb and sc phases with respect to the line @xmath232 ( or @xmath233 ) of the phase diagram .
we have also presented a series of nonselfdual phase portraits , which do not transform into themselves under the duality mapping ( see figs .
the results obtained lead to the conclusion that the growth of the chiral chemical potential @xmath1 promotes the chiral symmetry breaking , were studied , correspondingly , in the recent papers @xcite and @xcite .
it was shown there that the chiral chemical potential @xmath1 plays a role of the catalyst of dynamical chiral symmetry breaking . ]
whereas particle number chemical potential @xmath0 induces superconductivity in the system ( see also in @xcite ) .
we hope that our investigations may shed some new light on physical effects in planar systems like high - temperature superconductors or graphene .
the two - dimensional irreducible representation of the ( 2 + 1)-dimensional lorentz group so(2,1 ) is realized by the following @xmath281 @xmath282-matrices : @xmath283 acting on two - component dirac spinors .
they have the properties : @xmath284=-2i\varepsilon^{\mu\nu\alpha } \tilde\gamma_{\alpha};~ ~\tilde\gamma^{\mu}\tilde\gamma^{\nu}=-i\varepsilon^{\mu\nu\alpha } \tilde\gamma_{\alpha}+g^{\mu\nu } , \label{a2}\end{aligned}\ ] ] where @xmath285 .
there is also the relation : @xmath286 note that the definition of chiral symmetry is slightly unusual in ( 2 + 1)-dimensions ( spin is here a pseudoscalar rather than a ( axial ) vector ) .
the formal reason is simply that there exists no other @xmath281 matrix anticommuting with the dirac matrices @xmath287 which would allow the introduction of a @xmath28-matrix in the irreducible representation .
the important concept of chiral symmetries and their breakdown by mass terms can nevertheless be realized also in the framework of ( 2 + 1)-dimensional quantum field theories by considering a four - component reducible representation for dirac fields . in this case
the dirac spinors @xmath288 have the following form : @xmath289 with @xmath290 being two - component spinors . in the reducible four - dimensional spinor representation one
deals with 4@xmath2914 @xmath280-matrices : @xmath292 , where @xmath293 are given in ( [ a1 ] ) .
one can easily show , that ( @xmath294 ) : @xmath295 = diag(-i\varepsilon^{\mu\nu\alpha}\tilde\gamma_\alpha , -i\varepsilon^{\mu\nu\alpha}\tilde\gamma_\alpha ) .
\label{a5}\end{aligned}\ ] ] in addition to the dirac matrices @xmath296 there exist two other matrices , @xmath297 and @xmath28 , which anticommute with all @xmath296 and with themselves @xmath298 with @xmath299 being the unit @xmath281 matrix .
it is obvious that @xmath300 and @xmath301 .
each of the matrices @xmath297 and @xmath28 ( see appendix [ apa ] ) can be selected as a generator for the corresponding @xmath40 and @xmath35 chiral group of spinor field transformations .
for example , the lagrangian ( 1 ) is invariant with respect to @xmath35 such that @xmath302 .
alternatively , it is possible to construct a 4f model with fermion - antifermion and superconducting channels , symmetric under @xmath40 continuous chiral transformations , @xmath303 ( @xmath34 ) .
its lagrangian , e.g. , reads @xmath304\psi_k&+ & \frac { g_1}n\left ( \widetilde { 4f}\right ) _ { ch}+\frac { g_2}n\left ( \widetilde { 4f}\right ) _ { sc},\label{a7}\end{aligned}\ ] ] where the four - fermion structures @xmath305 and @xmath306 are used , @xmath307 here @xmath308 and @xmath0 is the usual particle number chemical potential ( as in ( 1 ) ) .
since this lagrangian is invariant under @xmath40 , there exist a corresponding conserved density of chiral charge @xmath309 as well as its thermodynamically conjugate quantity , the chiral ( or axial ) chemical potential @xmath310 .
note that in the framework of the model ( [ a7 ] ) there is also the duality between chiral and superconducting channels of interaction , i.e. between @xmath311 and @xmath312 four - fermion structures ( [ a71 ] ) .
indeed , performing in ( [ a7])-([a71 ] ) the modified pauli - gursey transformation of spinor fields , @xmath313 we see that @xmath314 the relations ( [ d2 ] ) are the basis for the duality invariance of the thermodynamic potential of the model ( [ a7 ] ) .
then , the thermodynamic potential @xmath322 , which describes the ground state of the model ( [ a7])-([a8 ] ) , is analogously defined by the following relation ( for convenience , we omit here and below the tilde in fields ) @xmath323[d\psi_l]\exp\big ( i\int d^3 x\big [ \bar\psi_k d_{\gamma^3}\psi_k\right.\nonumber\\ & & \left.- \frac{\delta}{2}(\psi_k^t\widetilde c\psi_k ) -\frac{\delta}{2}(\bar\psi_k\widetilde c\bar\psi_k^t)\big ] \big ) \right ) , \label{90}\end{aligned}\ ] ] where @xmath324 . in writing down this expression
, we used the same arguments that were used in the derivation of the thermodynamic potential in the framework of the model ( 1 ) ( see the derivation of ( [ 9 ] ) in section ii ) .
path integration in the expression ( [ 90 ] ) is evaluated in appendix [ apb ] so we have for the tdp ( [ 90 ] ) the following expression @xmath325 , \label{91}\end{aligned}\ ] ] where @xmath326 are presented in ( [ b14 ] ) . comparing this thermodynamic potential with the tdp @xmath327 ( [ 11 ] ) of the model ( 1 ) , we see that @xmath328 as a result , it is clear that the tdp @xmath329 of the 4f model ( [ a7 ] ) is invariant under the following dual transformation @xmath330 furthermore , to find phase portraits of the model ( [ a7 ] ) , it is sufficient to perform in all figs 1, .. ,10 the replacement @xmath331 .
let us calculate the following path integral over anticommuting four - component dirac spinor fields @xmath332 , @xmath333 : @xmath334[d\psi ] \exp\big ( i\int d^3 x\big [ \bar\psi   d_\gamma\psi-\frac{\delta}{2}(\psi^tc\gamma\psi ) -\frac{\delta}{2}(\bar\psi c\gamma\bar\psi^t)\big ] \big ) , \label{b1}\end{aligned}\ ] ] where @xmath335 is the charge conjugation matrix , @xmath336 .
since in the 4@xmath2914 spinor space the product matrix @xmath337 should be an antisymmetric one , it is clear that the @xmath338 matrix in ( [ b1 ] ) might be , e.g. , the unit matrix @xmath339 in the spinor space or @xmath340 ( for notations see appendix [ apa ] ) , etc . note in addition that at @xmath341 and @xmath342 ( see in ( [ 9 ] ) ) the integral @xmath343 is equal to the argument of the @xmath344-function in the formula ( [ 9 ] ) in the particular case @xmath98 , whereas at @xmath340 and @xmath345 ( see ( [ 90 ] ) ) the integral @xmath343 is equal to the argument of the @xmath344-function in the formula ( [ 90 ] ) also in the particular case @xmath98 .
recall , there are general gaussian path integrals @xcite : @xmath346\exp\big ( i\int d^3 x\big [ -\frac{1}{2}\psi^t a \psi + \eta^t\psi \big ] \big ) & = & \left(\det(a)\right ) ^{1/2}\exp\big ( -\frac{i}{2}\int d^3 x\big [ \eta^t a^{-1}\eta\big ] \big ) ,   \label{b2}\\ \int[d\bar\psi]\exp\big ( i\int d^3 x\big [ -\frac{1}{2}\bar\psi a \bar\psi^t + \bar\eta\bar\psi^t\big ] \big ) & = & \left ( \det(a)\right ) ^{1/2}\exp\big ( -\frac{i}{2}\int d^3 x\big [ \bar\eta a^{-1}\bar\eta^t\big ] \big ) , \label{b3}\end{aligned}\ ] ] where @xmath347 is an antisymmetric operator in coordinate and spinor spaces , and @xmath348 , @xmath349 are anticommuting spinor sources which also anticommute with @xmath288 and @xmath350 .
first , let us integrate in ( [ b1 ] ) over @xmath288-fields with the help of the relation ( [ b2 ] ) supposing there that @xmath351 , @xmath352 , i.e. @xmath353
. then @xmath354\exp\big ( -\frac{i}{2}\int d^3 x \bar\psi \big [ \delta c\gamma+d_\gamma(\delta c\gamma)^{-1}d_\gamma^t\big ] \bar\psi^t\big ) .
\label{b4}\end{aligned}\ ] ] second , the integration over @xmath350-fields in ( [ b4 ] ) can be easily performed with the help of the formula ( [ b3 ] ) , where one should put @xmath355 and @xmath356 . as a result
, we have @xmath357\right ) ^{1/2}=\left ( \det [ \delta^2(c\gamma)^2+d_\gamma ( c\gamma)^{-1 } d_\gamma^t c\gamma]\right ) ^{1/2}. \label{b5}\end{aligned}\ ] ] we evaluate the expression ( [ b5 ] ) in two cases , ( i ) when @xmath341 and @xmath342 , where @xmath358 is presented in ( [ 9 ] ) , and ( ii ) when @xmath340 and @xmath345 , where @xmath359 is given in ( [ 90 ] ) .    *
( i ) the case @xmath360 and @xmath361*. taking into account the relations @xmath362 and @xmath363 ( @xmath59 ) , we obtain from ( [ b5 ] ) @xmath364\right ) ^{1/2}\equiv\left ( \det b\right ) ^{1/2 } , \label{b7}\end{aligned}\ ] ] where @xmath365 . using the general relation @xmath366
, we get from ( [ b7 ] ) : @xmath367 ( a more detailed consideration of operator traces is presented in appendix a of the paper @xcite . ) in this formula the symbol tr means the trace of an operator both in the coordinate and internal spaces .
moreover , @xmath368 ( @xmath120 ) in ( [ b8 ] ) are four eigenvalues of the 4@xmath2914 fourier transformation matrix @xmath369 of the operator @xmath370 from ( [ b7 ] ) .
namely , @xmath371 where @xmath372 .
it is clear from these relations that @xmath101 are even functions vs @xmath0 and/or @xmath1 .
so , it is enough to take into account only nonnegative values of the chemical potentials @xmath0 and @xmath1 .    *
( ii ) the case @xmath373 and @xmath374*. in this case @xmath375 .
then it follows from ( [ b5 ] ) that @xmath376\right ) ^{1/2}\equiv\left ( \det { \cal b}\right ) ^{1/2 } , \label{b10}\end{aligned}\ ] ] where @xmath377 . then , similar to ( [ b8 ] ) , it is easy to find that in the case under consideration @xmath378 where @xmath379 ( @xmath120 ) are four eigenvalues of the 4@xmath2914 fourier transformation matrix @xmath380 of the operator @xmath381 from ( [ b10 ] ) .
it turns out that the eigenvalues @xmath379 are connected with the eigenvalues @xmath368 ( [ b12 ] ) by the simple relations @xmath382
to calculate the improper convergent integral in ( [ 31 ] ) we first use there a polar coordinate system , i.e. @xmath383 , and then restrict the @xmath384-integration region by @xmath137 ( suppose in addition that @xmath385 ) . as a result , we come to the regularized expression @xmath386 of the tdp ( [ 31 ] ) , @xmath387 , where @xmath388,\label{c9}\\ i_\pm&=&-\int_0^\lambda\frac{pdp}{2\pi}\left ( \mu-\sqrt{m^2+(p\pm\mu_5)^2}\right ) \theta\left ( \mu-\sqrt{m^2+(p\pm\mu_5)^2}\right ) . \label{c10}\end{aligned}\ ] ] due to the presence of the @xmath389-functions in ( [ c10 ] ) and sufficiently high values of the cutoff parameter @xmath390 , the quantities @xmath391
indeed do not depend on @xmath137 .
moreover , it is evident that @xmath392 is a nonzero quantity only in the case @xmath393 .
hence , substituting @xmath394 for the @xmath392-integration in ( [ c10 ] ) , we have @xmath395 ( the upper limit in the integral ( [ c12 ] ) corresponds to a value of the momentum @xmath384 , where the @xmath396-function of the integrand for @xmath392 from ( [ c10 ] ) is equal to zero , i.e. to a value of @xmath384 determined by the condition @xmath397 . )    it is convenient to present the quantity @xmath398 from ( [ c10 ] ) as a sum of two integrals , @xmath399 , where @xmath400 it is clear that @xmath401 is a nonzero quantity only at @xmath402 .
then , performing in the integral ( [ c14 ] ) the substitution @xmath403 , we obtain @xmath404 obviously , we have @xmath405 at @xmath406 .
to evaluate @xmath407 at other values of @xmath0 , we should , first , substitute @xmath408 in the integral ( [ c13 ] ) and then consider two different regions of the parameter @xmath0 , ( i ) @xmath409 , ( ii ) @xmath410 . as a result
, we have @xmath411 summing the expressions ( [ c12 ] ) , ( [ c15 ] ) and ( [ c16 ] ) , we have @xmath412    to calculate the @xmath413-term of the regularized tdp @xmath386 , it is useful to present @xmath413 as a sum of three integrals , each one corresponds to some square root expression of the integrand in ( [ c9 ] ) .
then , substituting @xmath394 and @xmath403 in two of these integrals , it is possible to present @xmath413 in the form @xmath414 since @xmath415 , we an use @xmath416 in the last two terms of ( [ c4 ] ) .
then , after integrations one can see that the sum of these two last integrals in ( [ c4 ] ) is equal to zero in the limit @xmath146 .
hence , @xmath417.\label{c11}\end{aligned}\ ] ] finally , performing in ( [ c17 ] ) trivial table integrations and using the relation @xmath418 we obtain from ( [ c0 ] ) , ( [ c17 ] ) and ( [ c11 ] ) the expression ( [ c18 ] ) ( the function @xmath419 is given in the text below ( [ 32 ] ) ) .      d.  ebert and m.k .
volkov , z.  phys .
c * 16 * , 205 ( 1983 ) ; d.  ebert , h.  reinhardt and m.k .
volkov , prog .
phys .   * 33 * , 1 ( 1994 ) .
m.  buballa , phys .
* 407 * , 205 ( 2005 ) .
y. nambu and g. jona - lasinio , phys .
rev.*122 * , 345 ( 1961 ) ; phys . rev .
* 124 * , 246 ( 1961 ) .
a. chodos and h. minakata , phys .
a * 191 * , 39 ( 1994 ) ; lect .  notes phys .
* 508 * , 231 ( 1998 ) [ hep - th/9709197 ] ; h.  caldas , j.l
.  kneur , m.b .
pinto and r.o .
ramos , phys .
b * 77 * , 205109 ( 2008 ) . v. schn and m. thies , at the frontier of particle physics : handbook of qcd :  boris ioffe festschrift  , vol . * 3 * , 1945 , world scientific ( 2001 ) [ hep - th/0008175 ] ; m.  thies , j.  phys .
a * 39 * , 12707 ( 2006 ) .
semenoff , i.a .
shovkovy and l.c.r .
wijewardhana , mod .
a * 13 * , 1143 ( 1998 ) .
zhukovsky , k.g .
klimenko , v.v .
khudyakov and d.  ebert , jetp lett .
* 73 * , 121 ( 2001 ) ; v.c .
zhukovsky and k.g .
klimenko , theor .  math
.  phys .   * 134 * , 254 ( 2003 ) ; e.j .
ferrer , v.p .
gusynin and v.  de la incera , mod .
b * 16 * , 107 ( 2002 ) ; eur .
j.  b * 33 * , 397 ( 2003 ) .
marino and l.h.c.m .
nunes , nucl .
b * 741 * , 404 ( 2006 ) ; l.h.c.m .
nunes , r.l.s .
farias and e.c .
marino , phys .
a * 376 * , 779 ( 2012 ) .
klimenko , r.n .
zhokhov and v.c .
zhukovsky , phys .
d * 86 * , 105010 ( 2012 ) .
klimenko , r.n .
zhokhov and v.c .
zhukovsky , mod .
lett .  a * 28
* , 1350096 ( 2013 ) .
d. mesterhazy , j. berges , and l. von smekal , phys .
b * 86 * , 245431 ( 2012 )
.    d.  ebert , k.g .
klimenko , p.b .
kolmakov and v.c .
zhukovsky , arxiv:1509.08093 [ cond-mat.mes-hall ] .
klimenko , z. phys .
c * 37 * , 457 ( 1988 ) ; b. rosenstein , b.j .
warr , and s.h .
park , phys .
d * 39 * , 3088 ( 1989 ) ; phys .
lett . * 62 * , 1433 ( 1989 ) .
klimenko , z.  phys .
c * 54 * , 323 ( 1992 ) ; theor .
* 89 * , 1161 ( 1992 ) ; theor .  math .
* 90 * , 1 ( 1992 ) ; v.p .
gusynin , v.a .
miransky and i.a .
shovkovy , phys .
lett .   *
73 * , 3499 ( 1994 ) ; a.s .
vshivtsev , b.v .
magnitsky and k.g .
klimenko , phys .  atom .
nucl .   * 57 * , 2171 ( 1994 ) ; theor .
phys .   * 106 * , 319 ( 1996 ) ; v.p .
gusynin , d.k .
hong and i.a .
shovkovy , phys .
d * 57 * , 5230 ( 1998 ) .
a.  chodos , h.  minakata , f.  cooper , a.  singh , and w.  mao , phys .
d * 61 * , 045011 ( 2000 ) .
d.  ebert , k.g .
klimenko and h.  toki , phys .
d * 64 * , 014038 ( 2001 ) ; h.  kohyama , phys .
d * 77 * , 045016 ( 2008 ) ; phys .
d * 78 * , 014021 ( 2008 ) .
d.  ebert , k.g .
klimenko , a.v .
tyukov and v.c .
zhukovsky , phys .
d * 78 * , 045008 ( 2008 ) ; d.  ebert , k.g .
klimenko , phys .
* d80 * , 125013 ( 2009 ) ; d.  ebert , n.v .
gubina , k.g .
klimenko , s.g .
kurbanov and v.c .
zhukovsky , phys .
d * 84 * , 025004 ( 2011 ) .
gubina , k.g .
klimenko , s.g .
kurbanov and v.c .
zhukovsky , phys .
d * 86 * , 085011 ( 2012 ) ; moscow univ .
phys .  bull .
* 67 * , 131 ( 2012 ) .
kneur , m.b .
pinto , r.o .  ramos and e.  staudt ,
d * 76 * , 045020 ( 2007 ) ; phys .  lett .
b * 657 * , 136 ( 2007 ) .
klimenko , z.  phys .
c * 50 * , 477 ( 1991 ) ; mod .
lett .  a * 9 * , 1767 ( 1994 ) .
i.  ojima and r.  fukuda , prog .
phys .   * 57 * , 1720 ( 1977 ) .
a.n .  vasiliev and g.y .
panasyuk , theor .
* 103 * , 570 ( 1995 ) [ teor
.  mat .
* 103 * , 295 ( 1995 ) ] .
m.  thies , phys .
d * 68 * , 047703 ( 2003 ) ; phys .
d * 90 * , no .
10 , 105017 ( 2014 ) .
d.  ebert , t.g .
khunjua , k.g .
klimenko and v.c .
zhukovsky , phys .
d * 90 * , 045021 ( 2014 ) .
n.d . mermin and h. wagner , phys .
* 17 * , 1133 ( 1966 ) ; s. coleman , commun .
31 * , 259 ( 1973 ) .
andrianov , d.  espriu and x.  planells , eur .
j.  c * 73 * , 2294 ( 2013 ) ; eur .
j.  c * 74 * , 2776 ( 2014 ) ; r.  gatto and m.  ruggieri , phys .
d * 85 * , 054013 ( 2012 ) ; m.  ruggieri , arxiv:1110.4907 ; l.  yu , h.  liu and m.  huang , phys .  rev .
d * 90 * , 074009 ( 2014 ) ; l.  yu , h.  liu and m.  huang , arxiv:1511.03073 [ hep - ph ] ; g.  cao and p.  zhuang , phys .
d * 92 * , 105030 ( 2015 ) .
braguta and a.y .
kotov , arxiv:1601.04957 [ hep - th ] .
zhukovsky , k.g .
klimenko and v.v .
khudyakov , theor .
phys .   * 124 * , 1132 ( 2000 ) [ teor .  mat .
* 124 * , 323 ( 2000 ) ] .
b.  vanderheyden and a.d .
jackson , phys .
d * 61 * , 076004 ( 2000 ) ; v.c .
zhukovsky , v.v .
khudyakov , k.g .
klimenko and d.  ebert , jetp lett .   * 74 * , 523 ( 2001 ) [ pisma zh .
teor .  fiz .
* 74 * , 595 ( 2001 ) ] .
g.  cao , l.  he and p.  zhuang , phys .
d * 90 * , 056005 ( 2014 ) .
a.n . vasiliev,``functional methods in quantum field theory and statistical physics '' , leningrad univ . press , leningrad , 1976 .