the two best studied variants of the massless gross - neveu ( gn ) model @xcite differ by their chiral symmetry groups .
the first model has the lagrangian @xmath3 throughout this paper , we are always in 1 + 1 dimensions and suppress contracted  color " labels on the fermion bilinears ( 1 ...
@xmath0 ) .
the lagrangian ( [ a1 ] ) is invariant under independent sign flips of left- and right - handed dirac fields , corresponding to the discrete chiral group z@xmath4z@xmath5 .
the second model has a continuous u(1)@xmath6u(1)@xmath2 chiral symmetry and is often referred to as two - dimensional nambu  jona - lasinio ( njl@xmath7 ) model @xcite , @xmath8 \qquad ( { \rm u(1)}{\text -}{\rm njl } ) .
\label{a2}\ ] ] here , phases of left- and right - handed fermion fields can be rotated independently . as a phenomenological model , the njl model in 3 + 1 dimensions
is in general considered with @xmath9 , but two additional flavors ( @xmath10 ) corresponding to isospin in strong interaction physics @xcite .
this leads to the lagrangian with su(2)@xmath6su(2)@xmath2 chiral symmetry @xmath11   \qquad ( { \rm su(2)}{\text -}{\rm njl } ) .
\label{a3}\ ] ] the same theory in 1 + 1 dimensions is also exactly solvable in the large @xmath0 limit , see the recent paper @xcite . as shown there and noticed before in the context of the phase diagram @xcite , the physics is closer to that of the gn model ( [ a1 ] ) than to that of the njl@xmath7 model ( [ a2 ] ) .
this reflects the following group - theoretical fact : the su(2 ) chiral symmetry also entails the discrete symmetry of the one - flavor gn model through the center of the group su(2 ) , but obviously not the full u(1 ) symmetry . if lagrangian ( [ a3 ] ) is the generalization of lagrangian ( [ a1 ] ) to two flavors , it is not too hard to identify the generalization of lagrangian ( [ a2 ] ) to two flavors : it is the njl - type model with u(2)@xmath1u(2)@xmath2 chiral symmetry and lagrangian @xmath12        \qquad ( { \rm u(2)}{\text -}{\rm njl } ) .
\label{a4}\ ] ] in which sense model ( [ a4 ] ) is the generalization of model ( [ a2 ] ) and model ( [ a3 ] ) the generalization of model ( [ a1 ] ) will be made more precise below , after having the solutions of all four models at our disposal .
the u(2)-njl model has been considered in 3 + 1 dimensions before @xcite . in 1 + 1 dimensions
, it has been explored using non - abelian bosonization at finite @xmath0 @xcite .
when going through the four models ( [ a1])([a4 ] ) , the complexity increases significantly , as evidenced by the number of ( real ) bilinear condensates which is doubled at each step ( @xmath13 ) . nevertheless , according to the experience with the one - flavor models ,
the solution of the u(2)-njl model is expected to be simpler than that of the su(2)-njl model .
model ( [ a4 ] ) is also interesting for another reason : recently , takahashi has generalized solutions of the bogoliubov - de gennes ( bdg ) equation to multicomponent form in the context of unconventional fermi superfluids @xcite .
this work has been found useful for solving twisted kink dynamics in the su(2 ) symmetric model @xcite .
it was noticed there that the elementary kink of takahashi , the simplest multi - fermion bound state , does not exist in the su(2 ) model , since it is not charge conjugation invariant .
for the same reason the kink of shei in the njl@xmath7 model @xcite does not appear in isolation in the gn model , only as a constituent of kink - antikink bound states ( dashen - hasslacher - neveu baryon @xcite ) .
it is then plausible , and will indeed be confirmed below , that takahashi s kink will appear in the model ( [ a4 ] ) as a physical state .    the present paper deals mostly with lagrangian ( [ a4 ] ) , but cross references to the other variants of the gn model will be frequent .
it is therefore necessary to adopt a simple naming of the four distinct models .
since our emphasis is on the chiral symmetry group , we will refer in the present work to the models ( [ a1 ] )  ( [ a4 ] ) as z@xmath7-gn , u(1)-njl , su(2)-njl and u(2)-njl models , respectively , as already indicated in the equations .
we also omit the subscript 2 for the number of dimensions from now on .
this paper is organized as follows .
[ sect2 ] develops the mean field theory of the u(2)-njl model , i.e. , the hartree - fock ( hf ) and time - dependent hartree - fock ( tdhf ) approaches . in sect .
[ sect3 ] , spontaneous chiral symmetry breaking , the dynamical fermion mass and the vacuum manifold are discussed .
the meson spectrum is the topic of sect .
[ sect4 ] , where the relativistic random phase approximation ( rpa ) is invoked to analyze small fluctuations around the hf vacuum .
massless hadrons , chiral spiral type condensates and the phase diagram as a function of @xmath14 are all related to the chiral anomaly and the fact that baryon number has a topological interpretation , as discussed in sect .
[ sect5 ] . in sect .
[ sect6 ] we generalize our previous solution of the tdhf equations from the u(1)-njl model to the two - flavor case , rederiving takahashi s recent work on the bdg equation in a slightly different language .
this has been done to make sure that the formalism from condensed matter physics really applies to the relativistic quantum field case at hand , and to facilitate comparison with previous works on other variants of gn or njl models @xcite .
a few simple applications to soliton problems with one and two bound states are given in sects .
[ sect7 ] and [ sect8 ] .
we finish with a short concluding section , sect .
[ sect9 ] .
the lagrangian ( [ a4 ] ) has the u(2)@xmath1u(2)@xmath2 chiral symmetry @xmath15 giving rise to 8 conserved noether currents @xmath16 in 1 + 1 dimensions , vector and axial vector currents are not independent , but satisfy @xmath17 adding and subtracting the conservation laws @xmath18 , and introducing light cone coordinates @xmath19 one finds @xmath20 if we take the expectation value of these equations in an arbitrary state , we conclude that the right - handed density @xmath21 depends only on @xmath22 , the left - handed density @xmath23 only on @xmath24 , i.e. , they can only move with the velocity of light to the right or to the left ( or be constant ) . in a localized , massive state like a solitonic bound state or breather
, these densities must therefore vanish identically .
the same argument goes through for the isovector currents and densities , since @xmath25 hence we anticipate that all densities and current densities must be zero inside an arbitrary soliton or multi - soliton state , a strong constraint on multi - fermion states following from chiral symmetry .    the basic equation in the large @xmath0 limit is the relativistic version of the time dependent hartree - fock ( tdhf ) equation . in the present case of lagrangian ( [ a4 ] ) , it reads @xmath26 \psi = 0 \label{b7}\ ] ] together with the self - consistency conditions @xmath27 we choose a chiral representation of the dirac matrices ( diagonal @xmath28 ) @xmath29 the tdhf equation in hamiltonian form then assumes the form @xmath30 here , the hamiltonian has been written in 2@xmath312-block - form , @xmath32 and @xmath33 are 2d isospinors , and @xmath34 is the 2@xmath312 matrix @xmath35 the self - consistency conditions ( [ b8 ] ) are particularly concise in the 2@xmath312 matrix form @xmath36 in light cone coordinates ( [ b4 ] ) the covariant form of the tdhf equation takes on the form @xmath37 a 2-component generalization of the equations familiar from the one - flavor case .
not only the dirac - tdhf equation ( [ b13 ] ) , but also the self - consistency condition ( [ b12 ] ) are manifestly preserved under chiral transformations @xmath38 that is to say that two solutions which differ only by a chiral transformation have to be considered as being physically indistinguishable .
if chiral symmetry is spontaneously broken but lorentz invariance is preserved , the vacuum is characterized by a space - time independent mean field @xmath39 . in order to decide whether this happens , we diagonalize the hf hamiltonian ( a @xmath40 matrix ) with homogeneous @xmath39 , eq .  ( [ b11 ] ) , in momentum space , @xmath41 the 4 eigenvalues are @xmath42 with dynamical fermion masses @xmath43 next we minimize the vacuum energy density , following closely the corresponding steps in the z@xmath7-gn model ( see e.g. @xcite ) , @xmath44 where @xmath45 is the single particle vacuum energy density of the z@xmath7-gn model , @xmath46 , \label{c5}\end{aligned}\ ] ] whereas @xmath47 denotes the double counting correction of the interaction energy density characteristic for the hf approach , @xmath48 minimizing @xmath49 with respect to @xmath50 yields the two conditions @xmath51 not surprisingly , the minimum is at @xmath52 , and we recover the gap equation of the z@xmath7-gn model with @xmath53 flavors .
the renormalized vacuum energy density is accordingly @xmath54 what is the vacuum manifold ?
the condition @xmath55 implies @xmath56 , eq .
( [ c3 ] ) , or @xmath57 the dynamical mass then becomes @xmath58 the vacuum potential matrix @xmath39 [ see ( [ b11 ] ) ] is given by @xmath59 and satisfies @xmath60 choosing units where @xmath61 , we first divide @xmath39 by @xmath62 from eq .
( [ c9 ] ) to get the u(2 ) matrix @xmath63 where the first factor belongs to u(1 ) , the second to su(2 ) .
the condition @xmath61 implies furthermore that @xmath64 since a u(2 ) matrix has 4 real parameters , but @xmath34 in ( [ c12 ] ) is parametrized by 5 real numbers ( @xmath65 ) , one condition is indeed necessary .
alternatively , it is obvious that expression ( [ c12 ] ) depends only on the 4 parameters @xmath66 .
thus the vacuum manifold is u(2 ) or @xmath67 . in cases where only a single vacuum is involved ,
the simplest choice is @xmath68 , which can always be achieved by a chiral transformation . for twisted field configurations
, we shall use @xmath68 for @xmath69 , but we then need the general expression ( [ c12 ] ) for @xmath70 .
in order to derive the spectrum of fermion - antifermion bound states ( mesons ) , one has to go beyond the hf approximation and consider small fluctuations around the hf vacuum .
the appropriate tool is the relativistic rpa .
although the technicalities are fairly involved , we can be brief here since one can follow almost literally the corresponding caculation in the z@xmath7-gn model , described in more detail in refs .  @xcite .
one starts from the equation of motion for the color singlet bilinear operator @xmath71 for earlier applications of this method , mainly to two - dimensional quantum chromodynamics , see also refs .  @xcite . in the present case , @xmath72
are combined dirac- and isospin indices ranging from 1 to 4 .
the @xmath40 matrix operator @xmath73 satisfies the equation of motion @xmath74 o_n q(x , y ) - q(x , y ) o_n { \rm tr } \left [ o_n q(y , y)\right ] \right\ } \label{d2}\end{aligned}\ ] ] where @xmath75 we have dropped terms irrelevant in the large @xmath0 limit . in the case of isovector operators @xmath76 , scalar products are implied in eq .  ( [ d2 ] ) . in the z@xmath7-gn model , only the @xmath77-term was present .
the next steps will not be shown in detail , but only enumerated :    1 .
expand the bilinear operator around the vacuum expectation value @xmath78 2 .
linearize the equation of motion in the fluctuation part @xmath79 .
3 .   choose the vacuum @xmath68 without loss of generality .
4 .   transform @xmath80 and @xmath81 to momentum space .
expand @xmath80 into vacuum spinors @xmath82 . in the large @xmath0 limit ,
only particle - hole components are important , @xmath83 this expansion is related to the dirac indices only , so that @xmath84 are still 2@xmath312 matrices consisting of isoscalar and isovector pieces , @xmath85 6 .
sandwich the equation of motion between the vacuum and one - meson states with momentum @xmath86 , @xmath87    the rpa - amplitudes @xmath88 then satisfy coupled linear integral equations of the standard rpa form .
what one finds is that the scalar and pseudoscalar mesons decouple , as do the 4 different isospin components ( isoscalar and isovector ) .
the kernel of the integral equations is 1-term separable , so that there can be at most one meson bound state per channel .
owing to the separability one can easily obtain an eigenvalue condition from the rpa equations .
all 4 scalar mesons satisfy the same equation leading to the eigenvalue condition @xmath89 identical to the one of the z@xmath7-gn model with @xmath53 flavors .
we use the notation @xmath90 and @xmath91 is the meson energy . for meson mass @xmath92 , the factor in curly brackets in eq .
( [ d8 ] ) becomes 1 and the whole equation reduces to the vacuum gap equation , hence the 4 scalar mesons all have the same mass @xmath92 . the 4 pseudoscalar mesons also satisfy identical equations , leading to the different eigenvalue condition @xmath93 from which we can read off a vanishing meson mass @xmath94 .
these are the would - be goldstone bosons , reflecting the 4 flat directions ( i.e. , the dimension of the vacuum manifold ) .
note that covariance is manifest in the relativistic rpa . since the meson spectrum has not been discussed in ref .
@xcite , let us mention that a corresponding calculation for the su(2)-njl model would have yielded one massive scalar / isoscalar ( @xmath92 ) and three massless pseudoscalar / isovector mesons .
for the present purpose , it is sufficient to consider the case where the internal rotation axis is  frozen " , say in the 3-direction , @xmath95 then the two isospin channels decouple , each one corresponding to a u(1)-njl model with @xmath0 flavors and mean field @xmath96 in the same way as z@xmath7-gn solutions are particular solutions of the su(2)-njl model @xcite , we thus find that u(1)-njl solutions are particular solutions of the u(2)-njl model .
the tdhf equations are evidently satisfied , as are the self - consistency conditions .
the vacuum can be chosen as @xmath68 . to describe dense matter , we combine two u(1)-chiral spirals @xcite in the two isospin channels to the following transformation of the vacuum spinors @xmath97 the two resulting chiral spirals have same radius ( @xmath61 ) but different pitches , @xmath98 ( @xmath99 sign for isospin up , @xmath100 sign for isospin down ) . in this manner
matter with arbitrary @xmath101 can be described most efficiently , @xmath102 in the low density limit we find massless , delocalized baryons with different isospin content .    without detailed calculation
, we can predict the phase diagram in ( @xmath103 space by using the known results for the u(1)-njl and su(2)-njl models @xcite . for temperatures above @xmath104 ,
chiral symmetry is restored ( @xmath105 ) and the fermions are massless for all ( @xmath106 ) . below @xmath107 the order parameter is expected to be @xmath108 ( @xmath109 is the dynamical fermion mass at @xmath110 ) with non - vanishing components @xmath111 and the densities and grand canonical potential density @xmath112 in the su(2)-njl@xmath7 model , we had a spatial modulation of the radius of the chiral spiral . here
the radius is spatially constant , but the pitch gets modulated .
notice the lack of periodicity of @xmath34 if @xmath113 and @xmath114 are incommensurate .
thus the phase diagram is extremely simple and closely reminiscent of the u(1)-njl model , see fig .
[ fig1 ] , in contrast to the more complicated phase diagram of the z@xmath7-gn and su(2)-njl models @xcite .
it is worth noting that the type of order parameter ( [ e6 ] ) has appeared before in the literature as variational ansatz for the su(2)-njl model in 3 + 1 dimensions @xcite . in our case ,
the fact that the chiral group is u(2 ) rather than su(2 ) is crucial for this simple ansatz to work quantitatively , see also ref .
@xcite for a comparison with the phase diagram of the su(2)-njl model in 1 + 1 dimensions .
a framework which enables us to solve soliton dynamics in the u(2)-njl model is available from the theory of fermionic superfluidity and superconductivity in the form of a general solution of the multicomponent bdg equation @xcite . in a previous work ,
we have outlined how to adapt this formalism to the su(2)-njl model @xcite . since it is inconvenient to mix two formalisms with different conventions , and since details of the symmetries of @xmath34 and occupation fractions may be somewhat different in condensed matter and particle physics , we proceed in this section to generalize the one - flavor formalism of refs .
@xcite to two flavors .
we are strongly guided by takahashi s work and will arrive at results equivalent to his results , but in such a way that the equations resemble the previous ones from the one - flavor case and where covariance is more manifest .
we hope that this will make applications to relativistic quantum field theories easier .    after inspecting takahashi s formalism
, one quickly discovers that the generalization of the framework of refs .
@xcite is rather straightforward . due to the additional flavor index ,
the notation becomes more cumbersome , otherwise the whole scheme remains practically unchanged .
this is true notably for the construction of transparent potentials where we follow closely the logic of ref .
@xcite .    in @xcite , the starting point for attacking @xmath115 soliton problems was a @xmath115-dimensional vector @xmath116 with components @xmath117 the @xmath118 are complex numbers ( i m @xmath119 ) characterizing the pole positions of the tdhf continuum wave functions in the complex @xmath120 plane , @xmath120 being the spectral parameter related to light cone momentum and energy ( uniformizing parameter in condensed matter language ) , @xmath121 note that @xmath122 so that @xmath123 is recognized as a plane wave evaluated at a complex spectral parameter corresponding to a bound state pole .
the main generalization when going from one to two flavors consists in introducing two copies of each @xmath123 , @xmath124 here , @xmath125 is a 2-component , constant , complex vector with components @xmath126 .
its precise meaning will be clarified later on when we investigate few - soliton problems in more detail .
the generalization to more than 2 flavors is straightforward , but not needed here . without loss of generality , we can assume that the vectors @xmath125 are normalized ( @xmath127 ) .
we shall use greek indices and the summation convention for flavor ( @xmath128 ) and suppress the indices @xmath129 referring to the bound state poles whenever possible .
then we have to modify the equations of refs .  @xcite as follows : continuum tdhf spinors are now 4-component objects represented as @xmath130 the following ansatz is inspired by the assumed pole structure of the continuum spinors ( @xmath115 poles , corresponding to @xmath115 bound states ) @xmath131 the @xmath132 are the amplitudes of the flavor components of the incoming plane wave @xmath133 when summing over all continuum states , they should be chosen as @xmath134 and @xmath135 to account for incoming waves in the two isospin channels ( at least , if the vacuum at @xmath136 is chosen as @xmath137 ) .
the @xmath138 are closely related to bound state wave functions .
they can be evaluated by linear algebra as follows : define a hermitean @xmath139 matrix @xmath140 , @xmath141 as in refs .
@xcite , we suppress soliton indices ( @xmath142 ) to ease the notation .
the @xmath143 satisfy the following system of linear , algebraic equations @xmath144 where @xmath145 .
a constant , hermitean @xmath139 matrix @xmath146 encoding further information about the soliton configuration ( geometry , initial conditions , details about time dependence for breathers ) has been introduced .
the dimension of the linear systems does not increase with the number of flavors , but depends only on the total number of bound state poles .
what is new is the factor @xmath147 in @xmath148 and the fact that one needs to solve the linear equation for each flavor component @xmath149 in turn , a rather mild complication .
notice that the pairs of exponentials ( @xmath150 ) provide us with @xmath151 non - normalizable solutions of the free , massive dirac equation ( @xmath61 ) , @xmath152 derivatives of @xmath140 are now 2-term separable , as opposed to 1-term separable before , @xmath153 differentiating ( [ f9 ] ) with respect to @xmath154 yields @xmath155 upon applying @xmath156 from the right , we arrive at @xmath157 this is just the 2-component version of the covariant tdhf equation with the following 2@xmath312 matrix potential @xmath158 turning to the continuum states , the dirac equation for @xmath159 introduced in ( [ f5 ] ) assumes the form @xmath160 we write down the three terms of each of these equations in detail . from the first line of ( [ f15 ] ) , @xmath161 as in the one - flavor case , terms containing @xmath162 in the first and 3rd line cancel . terms containing @xmath163 in the first and 2nd lines add up to @xmath164 and cancel the remaining terms in the 2nd and 3rd line .
likewise , from the 2nd line of ( [ f15 ] ) , @xmath165 cancellations work exactly as for one flavor : terms containing @xmath163 in the 1st and 3rd lines cancel .
terms containing @xmath162 in the 1st and 2nd lines add up to @xmath166 and cancel the remaining terms in 2nd and 3rd lines .
next we consider the question of normalization and orthogonality of the bound state spinors .
using @xmath167 we show that @xmath168 this yields the same overlap matrix as in the one flavor case , @xmath169 orthonormal bound states can be constructed as before by linear combinations of the @xmath170 , @xmath171 the resulting condition coincides with the one in the one - flavor case , @xmath172 we now turn to the self - consistency condition .
the mean field @xmath173 receives contributions from the sea and the bound states , @xmath174 with @xmath175 to evaluate the continuum part , we insert the @xmath176 s and integrate over @xmath177 with a cutoff .
the pole at @xmath178 yields the divergent contribution @xmath179 owing to the vacuum gap equation , this part gives self - consistency by itself , as usual in many variants of gn models .
the convergent part of the sea contribution can be written down most easily if one introduces a diagonal matrix @xmath180 , @xmath181 we find @xmath182 the bound state contribution , expressed through the diagonal matrix @xmath115 with matrix elements proportional to bound state occupation fractions , @xmath183 is given by @xmath184 expressions ( [ f27 ] ) and ( [ f29 ] ) cancel provided that @xmath185 remarkably , the self - consistency condition is not affected at all when going from one to two flavors .    as outlined in sect .
[ sect2 ] , due to current conservation , the solitons described by this formalism should have identically vanishing baryon density ( @xmath186 ) and isospin densities ( @xmath187 ) .
this should hold for left- and right - handed fermions separately , or , equivalently , for charge and current densities .
it is a good test of the formalism to prove this in all generality .
define @xmath188 introducing @xmath189 in addition to the ordinary su(2 ) generators @xmath190 , we can treat the baryon density and isovector densities on the same footing .
the densities comprise a part from the dirac sea and a part from the bound states .
the sea part has to be subtracted by the vacuum contribution for the case @xmath191 only ( baryon density ) .
the sea contribution is evaluated as follows .
we take into account the measure @xmath192 and write the total density as the integral @xmath193 here , @xmath194 are the contribution from the continuum state with spectral parameter @xmath120 to the densities , including the measure . inserting the continuum spinors and summing over the two isospin channels , we find @xmath195 we have introduced the vector @xmath196 with @xmath115 components @xmath197 generalizing the corresponding vector @xmath198 in the one - flavor case .
consider the quadratic terms in @xmath199 in eq .
( [ f33 ] ) first . for @xmath200 ,
use the following identity : @xmath201 eq .
( [ f9 ] ) in the form @xmath202 then serves to eliminate @xmath148 from ( [ f35 ] ) , thereby generating two terms which cancel exactly the terms linear in @xmath162 in ( [ f33 ] ) .
the remainder yields @xmath203 or , with the help of the diagonal matrix @xmath204 more compactly @xmath205 along the same lines , one finds for the right - handed density @xmath206 the bound state contribution in a notation similar to ( [ f31 ] ) is @xmath207 since @xmath208 with the matrix @xmath115 defined in ( [ f28 ] ) , the condition which guarantees vanishing of all 4 left- and right - handed densities is @xmath209 upon performing the integration over @xmath177 and noting that @xmath210 we recover the self - consistency condition ( [ f30 ] ) .
this shows that the physical solutions indeed have identically vanishing densities .
the result is non - trivial in the sense that it comes about as a cancellation between non - vanishing contributions from the dirac sea and the valence bound states and only holds for self - consistent solutions .
hence it may be viewed as a complete shielding effect
choose @xmath211 ( one soliton ) in the general formalism and work it out .
after setting @xmath212 we find the mean field @xmath213 this is the most elementary type of kink interpolating between two vacua . in the rest frame , @xmath214 for @xmath136 , @xmath215 .
for @xmath216 , @xmath217 with @xmath218 by a unitary transformation we can map @xmath219 onto the vector @xmath220 , in which case we recover the standard u(1)-njl kink in the isospin - up channel and the vacuum @xmath221 in the diagonal form @xmath222 hence we have a twisted kink in the isospin up channel and the vacuum for isospin down . the general expression ( [ g4 ] )
is nothing but the spectral representation of the unitary matrix @xmath221 . denoting the ( normalized ) eigenvectors corresponding to the eigenvalues @xmath223 and 1 by @xmath219 and @xmath224 , respectively
, we have @xmath225 eliminating the second term with the help of the completeness relation @xmath226 then reproduces ( [ g4 ] ) .
this shows that the vector @xmath219 can be interpreted as follows : @xmath219 is the eigenvector of the unitary matrix @xmath221 , the vacuum at @xmath216 , belonging to the  twisted " eigenvalue @xmath223 . with one bound state
, only one eigenvalue can be twisted .
this is the fundamental twisted kink .
since we can reduce the simplest kink problem to the kink of the u(1)-njl model , the spinors , self - consistency and vanishing baryon and isospin density follow from one - flavor results and need not be repeated here .
this is exactly the kink of takahashi @xcite .
it is special in the following sense : a generic twisted kink in the u(2)-njl model should connect the vacuum @xmath137 with the vacuum @xmath221 , a unitary matrix with eigenvalues @xmath227 . in the present framework
this requires a bound state of two elementary kinks , see next section .
this was already used in a previous work on the su(2)-njl model @xcite , where the eigenvalues have to be complex conjugates @xmath228 .
it is then necessary to bind a kink with its charge conjugate .
no such restriction exists in the u(2)-njl model .
the elementary kink exists as a physical state in the u(2)-njl model , whereas it is confined in the su(2)-njl model .
this is completely analoguous to what happens with shei s twisted kink @xcite in the u(1)-njl and z@xmath7-gn models , respectively @xcite .
we start with the scattering of two elementary kinks of the type discussed in the preceding section .
this requires a diagonal matrix @xmath146 .
we set @xmath229 where @xmath230 will be interpreted below .
then the result of using the general formalism for two poles ( @xmath231 ) follows the same pattern as kink - kink scattering in the u(1)-njl model or su(2)-njl model .
the mean field can be cast into the form @xmath232 the interpretation of the u(2 ) matrices @xmath233 and of the factor @xmath230 follows upon considering the asymptotics of the scattering process , @xmath234 thus @xmath230 accounts for the time delay of the solitons during the collision , and we have made explicit the intrinsic form of the soliton by pulling out the factors @xmath235 ( for @xmath236 ) and @xmath237 ( for @xmath238 ) to the right .
one recognizes @xmath239 and @xmath240 as the intrinsic twist factors of the participating solitons .
the formalism yields the following results @xmath241 in this formulation , the normalized vectors @xmath242 are input and specify the twist of incoming kink 1 and outgoing kink 2 .
the normalized vectors @xmath243 of the twist of outgoing kink 1 and incoming kink 2 can be expressed in terms of @xmath242 and the pole positions @xmath244 as follows , @xmath245 , \nonumber \\ \vec{q}_2 & = & { \cal n}_2 \left [   \sigma_{12 } \zeta_2^ * ( \zeta_1-\zeta_1^ * ) \vec{p}_1 - \zeta_1^ * ( \zeta_1-\zeta_2^ * ) \vec{p}_2 \right ] .
\label{h5}\end{aligned}\ ] ] we have used the notation @xmath246 the normalization factors ensuring that @xmath247 are given by @xmath248 here , @xmath230 is the time delay factor also appearing in eqs .
( [ h2],[h3 ] ) for which we find @xmath249 let us consider some interesting limiting cases .
obviously , the modulus of @xmath250 controls the strength of the kink - kink interaction between the two kinks .
if @xmath242 are parallel , @xmath251 and everything can be rotated into the isospin up channel where we recover scattering of two standard kinks of the u(1)-njl model , with the known time delay factor @xmath252 if @xmath242 are orthogonal , @xmath253 and the two kinks live in the two isospin channels without any interaction ( @xmath254 ) . by varying the angle spanned by @xmath219 and @xmath224 , we can thus vary the strength of the interaction between these two extreme cases .
a special case of the kink - kink scattering problem is the bound state .
this may be viewed as the generic u(2)-twisted kink , since it enables us to connect an arbitrary vacuum @xmath221 at @xmath216 to @xmath137 at @xmath136 . in the rest frame of such a composite kink
, we have to choose @xmath255 some simplifications occur , @xmath256 \nonumber \\ \vec{q}_2 & = & - { \cal n}_2 \left [ 2i \sigma_{12 } e^{i\phi_2 } \sin \phi_1 \vec{p}_1 - \left ( e^{i(\phi_1+\phi_2)}-1 \right ) \vec{p}_2 \right ] \nonumber \\ { \cal n}_1^{-2 } & = & { \cal n}_2^{-2 } = 2 \left ( 1- \cos ( \phi_1+\phi_2 ) - 2 \sigma_{12 } \sigma_{21 } \sin \phi_1 \sin \phi_2 \right ) \nonumber \\ \xi^{-1 } & = & 1 - \sigma_{12 } \sigma_{21 } \frac{2 \sin \phi_1 \sin \phi_2}{1-\cos(\phi_1+\phi_2 ) } \label{h11}\end{aligned}\ ] ] in the one flavor case , it was not possible to construct a bound state of two kinks with the same twist , @xmath257 .
if one specializes the formalism to this case , one recovers a single kink .
in the two flavor case , this restriction does not exist anymore and we can construct a non - trivial bound state out of two constituents with the same twist . in this case , eq .
( [ h11 ] ) reduces to @xmath258 finally , we turn to the twisted breather .
a breather at rest can be generated by choosing @xmath259 and a non - diagonal matrix @xmath146 .
the twisted breather is even more complicated here than in the single flavor case , so we refrain from discussing all possible parameter choices and illustrate the reult for a few simple special cases only .
we choose the same matrix @xmath146 than for one flavor @xcite , @xmath260 and the vectors @xmath242 as @xmath261 thus @xmath262 is the angle between @xmath219 and @xmath224 , @xmath263 if we choose @xmath219 and @xmath224 to be parallel ( @xmath264 ) , we find that @xmath265 where @xmath266 is the twisted breather in the one - flavor case @xcite , @xmath267 choosing @xmath219 and @xmath224 to be orthogonal ( @xmath268 ) one finds @xmath269 with @xmath270 here , the diagonal components @xmath271 are static , whereas the off - diagonal components @xmath272 oscillate with the same frequency as the one - flavor breather , @xmath273 for any other choice of the angle @xmath262 , all components of @xmath34 start to oscillate with the same frequency but different phases .
the corresponding more complicated expressions for @xmath34 can easily be generated using the general framework , but will not be given here .
in this paper , we have studied a variant of the integrable gross - neveu model family which has not yet received much attention so far in 1 + 1 dimensions : the u(2)-njl model with u(2)@xmath1u(2)@xmath2 chiral symmetry .
phenomenologically , in 3 + 1 dimensions , the su(2)@xmath1su(2)@xmath2 model is more relevant .
however , from a theoretical point of view , it is quite instructive to include the u(2 ) model as well . after developing the mean field approach and setting up the tdhf equation for the u(2)-njl model ,
we have studied the vacuum and identified the vacuum manifold as u(2 ) .
small fluctuations in the 4 flat directions give rise to 4 massless pseudoscalar mesons , whereas the other 4 directions yield massive scalars right at threshold @xmath274 .
the chiral spiral construction is especially simple here , as it can be invoked both for baryonic and isospin charge .
the phase diagram is extremely simple , generalizing the known u(1)-njl phase diagram to the ( @xmath275 ) space .
perhaps the most interesting topic is that of soliton dynamics .
here it turns out that a recently developed solution of the matrix bdg system in condensed matter physics fits perfectly the u(2)-njl model .
we have rederived takahashi s results in the language of our previous work on the u(1)-njl model , and confirmed that the fundamental kink is a physical state in the u(2 ) model while appearing only as a confined constituent in the su(2 ) model .
this is another striking example for the close relationship between relativistic quantum field theory toy models and sophisticated , realistic condensed matter problems .    in the present work
together with ref .
@xcite , we have generalized the well - known z@xmath7-gn and u(1)-njl models with abelian chiral groups to two flavors and non - abelian chiral groups su(2 ) and u(2 ) .
the relationship between these various models is summarized in table i.    .relationship among four - fermion models as discussed in main text . [ cols="^,^",options="header " , ]     the first row contains the original variants of the gn model @xcite with one flavor only and discrete or continuous chiral symmetry , respectively .
the most conspicuous differences between these two models are the phase diagram in the ( @xmath276 ) plane and the role played by twisted kinks , which are free in the u(1 ) model but confined into bound states in the z@xmath7 model .
this pattern repeats itself in the 2nd row , the two - flavor generalizations .
in fact , models in the same column share identical phase diagrams in the ( @xmath276 ) plane , as first noticed in ref .
@xcite for the first column .
a further generalization to the groups su(@xmath115 ) and u(@xmath115 ) is straightforward , thereby extending the spectrum of integrable quantum field theories substantially .
it is plausible ( but ought to be checked in future work ) that the marked differences in the phase diagram and in the role of twisted kinks show up for arbitrary numbers of flavors as well .
our final remark concerns the relationship between the leading order large @xmath0 results and the method of non - abelian bosonization @xcite .
it is well known that two - dimensional multi - color and -flavor njl models at any finite @xmath277 can be mapped onto decoupled bosonic field theories of wess - zumino - novikov - witten type @xcite and free bosonic fields .
recently there has been a lot of progress in solving such models using conformal field theory techniques and numerical methods , even for the case when chiral symmetry is explicitly broken by a mass term and more than one coupling constants , see ref .
@xcite and references therein .
it would be interesting to compute space - time dependent condensates in these integrable models at finite @xmath278 and compare them with the results in the limit @xmath279 , as this limit is expected to be rather singular .
the author thanks falk bruckmann for his interest in this work and useful discussions .
99 d. j. gross and a. neveu , phys . rev .
d * 10 * , 3235 ( 1974 ) .
y. nambu and g. jona - lasinio , phys . rev . *
122 * , 345 ( 1961 ) .
y. nambu and g. jona - lasinio , phys . rev . *
124 * , 246 ( 1961 ) .
m. thies , phys .
d * 93 * , 085024 ( 2016 ) .
a. heinz , f. giacosa , m. wagner , d. h. rischke , phys .
d * 93 * , 024512 ( 2016 ) .
d. ebert and m. k. volkov , z. phys .
c * 16 * , 205 ( 1983 ) .
d. ebert and h. reinhardt , nucl .
b * 271 * , 188 ( 1986 ) .
d. ebert , m. nagy , m. k. volkov , phys .
atom . nucl . * 59 * , 140 ( 1996 ) .
a. m. tsvelik , sov .
jetp * 66 * , 754 ( 1987 ) .
d. a. takahashi , phys .
b * 93 * , 024512 ( 2016 )
. d. a. takahashi , prog .
2016 * , 043i01 .
shei , phys .
d * 14 * , 535 ( 1976 ) .
r. f. dashen , b. hasslacher , a. neveu , phys .
d * 12 * , 2443 ( 1975 ) .
g. v. dunne and m. thies , phys .
d * 89 * , 025008 ( 2014 ) .
g. v. dunne and m. thies , phys .
* 111 * , 121602 ( 2013 ) .
g. v. dunne and m. thies , phys .
a * 88 * , 062115 ( 2013 ) . v. schn and m. thies , _ at the frontiers of particle physics : handbook of qcd , boris ioffe festschrift _ , edited by m. shifman ( world scientific , singapore , 2001 ) , vol .
3 , p. 1945 .
r. pausch , m. thies , v. l. dolman , z. phys .
a * 338 * , 441 ( 1991 ) .
w. r. gutierrez , nucl .
b * 176 * , 185 ( 1980 ) .
m. cavicchi , p. di vecchia , i. pesando , mod .
lett . a * 8
* , 2427 ( 1993 ) ; erratum : mod .
a * 8 * , 2909 ( 1993 ) .
r. de mello koch , j. p. rodrigues , phys .
d * 54 * , 7794 ( 1996 ) .
v. schn and m. thies , phys .
d * 62 * , 096002 ( 2000 ) .
g. basar , g. v. dunne , m. thies , phys .
d * 79 * , 105012 ( 2009 ) .
m. thies and k. urlichs , phys .
d * 67 * , 125015 ( 2003 ) .
d. nowakowski , m. buballa , s. carignano , j. wambach , arxiv:1506.04260 [ hep - ph ] .
e. witten , comm .
phys . * 92 * , 455 ( 1984 ) .
j. wess and b. zumino , phys .
b * 37 * , 95 ( 1971 ) .
s. p. novikov , sov . math
* 24 * , 222 ( 1981 ) .
p. azaria , r. m. konik , ph .
lecheminant , t. plmai , g. takcs , a. m. tsvelik , arxiv:1601.02979 [ hep - th ]