the study of quantum chromodynamics at high density is almost as old as the theory itself @xcite . in recent years
the field has attracted wide interest in the wake of a revival of the idea of color superconductivity ( csc ) @xcite . the stimulus for this revival was the observation @xcite that the instanton - induced quark
quark interaction can be much stronger than that induced by simple one - gluon exchange , and can thus give a transition temperature on the order of 100  mev .
subsequent work @xcite showed that the perturbative color - magnetic interaction also gives rise to a strong pairing interaction .    these and other dynamical considerations @xcite underlie a picture of the ground state of high - density qcd in which the @xmath3 gauge symmetry is spontaneously broken by a bcs - like condensate .
the details of the breaking , which include both the higgs ( or meissner ) effect and the rearrangement of global symmetries and goldstone bosons , depend on quark masses , chemical potentials , and temperature .
prominent in the list of possibilities are those of color - flavor locking in three - flavor qcd @xcite and crystalline superconductivity  with broken translation invariance  when there are two flavors with different densities @xcite .
for a review see @xcite .    as noted , csc at high density is so far a prediction of weak - coupling analysis .
one expects the coupling to become weak only at high densities , and in fact it turns out that reliable calculations demand extremely high densities @xcite .
the use of weak - coupling methods to make predictions for moderate densities is thus not an application of qcd , but of a model based on it .
it is imperative to confirm these predictions by non - perturbative methods .
standard lattice monte carlo methods , unfortunately , fall afoul of well - known technical problems when the chemical potential is made nonzero , although we do note remarkable progress made recently in the small-@xmath4 regime @xcite .    in this paper
we initiate a study of high - density quark matter based on lattice qcd in the strong - coupling limit .
we work in the hamiltonian formalism , which is more amenable than the euclidean formalism to strong - coupling perturbation theory and to qualitative study of the ensuing effective theory @xcite .
the fermion kinetic hamiltonian is a perturbation that mixes the zero - flux states that are the ground - state sector of the electric term in the gauge hamiltonian . in second order
, it moves color - singlet fermion pairs around the lattice ; the effective hamiltonian for these pairs is a generalized antiferromagnet , with spin operators constructed of fermion bilinears .
we depart from studies of the vacuum by allowing a background baryon density , which is perforce static in second order in perturbation theory .
our aim at this stage is to discover the ground state of the theory with this background . in third order ( when @xmath5 ) the baryons become dynamical ; we display the effective hamiltonian but make no attempt to treat it .
the symmetry group of the effective antiferromagnet is the same as the global symmetry group of the original gauge theory .
this depends on the formulation chosen for the lattice fermions .
following @xcite , we begin with naive , nearest - neighbor fermions , which suffer from species doubling @xcite and possess a global @xmath0 symmetry group that contains the ordinary chiral symmetries [ as well as the axial @xmath6 as subgroups .
we subsequently break the too - large symmetry group with next - nearest - neighbor ( ) couplings along the axes in the fermion hopping hamiltonian .
a glance at the menu of fermion formulations reveals the reasons for our choice .
wilson fermions @xcite have no chiral symmetry and make comparison of results to continuum csc difficult if not impossible . staggered fermions @xcite likewise possess only a reduced axial symmetry while suffering a reduced doubling problem .
the overlap action @xcite is non - local in time and hence possesses no hamiltonian ; attempts @xcite to construct an overlap hamiltonian directly have not borne fruit . finally , domain - wall fermions @xcite have been shown @xcite to lose chiral symmetry and regain doubling when the coupling is strong .
as we discuss below , while the  theory still exhibits doubling in the free fermion spectrum , we are not interested in the perturbative fermion propagator but in the spectrum of the confining theory .
we take it as a positive sign that the unbroken symmetry is now @xmath7 .
this symmetry is what we want for the continuum theory , except for the axial @xmath8 .
the latter can still be broken by hand .
our emphasis on the global symmetries is a consequence of the fact that the gauge field is not present in the ground - state sector and does not reappear in strong - coupling perturbation theory . in other words , confinement is a _
kinematic _ feature of the theory , leaving no possibility of seeing the higgs - meissner effect directly .
this is but an instance of confinement - higgs duality , typical of gauge theories with matter fields in the fundamental representation @xcite .
our aim is thus to identify the pattern of spontaneous breaking of global symmetries .
for various values of @xmath9 and @xmath10 , this can be compared to weak - coupling results @xcite .
this paper is largely an exposition of formalism , along with partial results .
we study the nearest - neighbor antiferromagnetic hamiltonian , both with and without a uniform baryonic background density .
we transform the quantum hamiltonian into a path integral for a non - linear @xmath11 model , where the manifold of the @xmath11 field depends on the baryon background .
we then investigate the limit of large @xmath9 and show that the global @xmath0 symmetry is indeed spontaneously broken .
adding in the  couplings is a problem of vacuum alignment @xcite .
we do this in the vacuum sector and recover the result @xcite that the @xmath1 chiral symmetry is broken to the vector @xmath12 .
the analysis for the finite - density theory is more involved and we defer it to a future publication .
other groups have recently studied the strong - coupling effective hamiltonian for naive and wilson fermions at non - zero chemical potential @xcite .
we differ from their approaches in eschewing mean field theory in favor of the exact transformation to the @xmath11 model , which is amenable to semiclassical treatment .
as noted above , we base our program on  fermions ; we also work at fixed baryon density .
let us walk through the paper .
we review in sec .
[ sec : heff ] the derivation of the effective hamiltonian of lattice gauge theory in strong - coupling perturbation theory @xcite .
the second - order hamiltonian [ @xmath13 is an antiferromagnet with @xmath0 spins ; the global symmetry group is @xmath0 for the nearest - neighbor theory , broken to @xmath1 by  terms .
the baryon number at each site determines the representation of @xmath0 carried by the spin at that site . in second order ,
baryon number is static ; it becomes mobile in the next order , where ( for @xmath5 ) the new term in the effective hamiltonian is a baryon hopping term .
the baryon operators responsible for the hopping are composite operators that do not obey canonical anti - commutation relations .
if this were not the case , then the effective hamiltonian in third order would strongly resemble that of the @xmath14@xmath15 model @xcite , @xmath16 here @xmath17 is an annihilation operator for an electron at site @xmath18 with spin @xmath19 , and the number operators @xmath20 and spin operators @xmath21 are constructed from it .
the added term @xmath22 is a more complicated hopping and interaction term .
the @xmath14@xmath15 model describes a doped antiferromagnet ; it arises as the strong - binding limit of hubbard model , a popular model for itinerant magnetism and possibly for high-@xmath23 superconductivity .
the model is not particularly tractable and , absent new theoretical developments , does not offer much hope for progress in our finite - density problem .
it is nonetheless worth pondering the fact that a model connected , however tentatively , with superconductivity appears in a study of high - density nuclear matter .    in the remainder of this paper
, we work only to @xmath24 , where the baryons are fixed in position .
motivated by the similarity of our hamiltonian to the heisenberg antiferromagnet , we apply condensed matter methods developed for that problem .
indeed , condensed matter physicists have generalized the @xmath25 , spin-@xmath26 heisenberg model to @xmath27 in many representations @xcite , which corresponds to adding flavor and color degrees of freedom to the electrons . plane . ]
these are exactly the generalizations needed for our effective hamiltonian . with @xmath9 colors and @xmath28 ( single - component ) flavors , a site of the lattice can be constrained to contain a color - singlet combination of @xmath29 particles .
the flavor indices of the spin then make up a representation of @xmath27 whose young diagram has @xmath9 columns and @xmath30 rows ( see fig .
[ fig : young1 ] ) .
we set @xmath31 and the correspondence is complete ( until we include  terms in the hamiltonian ) .    in sec .
[ sec : nls ] we derive a @xmath11 model representation for the partition function of the antiferromagnet . following read and sachdev @xcite ,
we employ spin coherent states @xcite to define the @xmath11 field . @xmath28 and
@xmath30 determine the target space of the @xmath11 model to be the symmetric space @xmath32 $ ] ; the number of colors @xmath9 becomes an overall coefficient of the action . multiplies the quantum hamiltonian , and hence serves only to set the energy scale .
] as for the quantum hamiltonian , the nearest - neighbor theory is symmetric under @xmath33 while the  terms break the symmetry to @xmath1 ( while leaving the manifold unchanged ) .
the @xmath9 multiplying the action invites a large-@xmath9 analysis , and in sec .
[ sec : zerob ] we study the vacuum sector , meaning zero baryon number , thereby .
we return to an exercise proposed and solved by smit @xcite , in generalizing the vacuum sector to allow baryon number @xmath34 on alternating sites ; this means specifying conjugate representations of @xmath33 on alternating sites , with respectively @xmath30 and @xmath35 rows .
as shown by read and sachdev @xcite , in this situation one can carry out an alternating @xmath33 rotation to convert the antiferromagnet into a ferromagnet with _
identical _ spins on alternating sites , and the classical ( @xmath36 ) analysis gives a homogeneous ground state .
the result is , as one might expect , that @xmath33 is broken to @xmath37 in the classical vacuum ; the ground state energy is independent of @xmath30 .
the @xmath2 corrections to the energy do depend on @xmath30 , however , and they select the self - conjugate @xmath38 configuration ( _ i.e. _ , @xmath39 everywhere ) as the lowest - energy vacuum .
thus the true ground state breaks @xmath40 .
calculation includes the effect of the time - derivative terms in the action that were dropped in the leading order , and these terms _ do _ remember the difference between the ferromagnet and the antiferromagnet .
] when we add  terms to the action as a perturbation , we find that the ground state breaks @xmath41 to the vector @xmath12 , as expected .
we turn to non - zero baryon density in sec .
[ sec : nonzerob ] .
we study homogeneous states , in which all sites carry the same representation of @xmath33 , with @xmath42 .
the classical vacuum of the @xmath11 model is more elusive than for the vacuum sector , since now there are identical manifolds on adjacent sites but the coupling is antiferromagnetic .
we begin by studying the two - site problem , and we learn that when one of the classical spins is fixed then when the energy is minimized the other spin is still free to wander a submanifold of the original symmetric space . if we replicate this to the infinite lattice then we have a situation where the even spins , say , are fixed in direction while each odd spin wanders the submanifold , independent of the other odd spins .
this means a ground state whose degeneracy is exponential in the volume , similar to some frustrated models or the antiferromagnetic potts model @xcite .
the cure to this disease comes from the @xmath43 fluctuations , which couple the odd spins to each other and make them align . in the end we find that the @xmath33 symmetry is broken by the vacuum to @xmath44 .
perturbing this ground state with the  terms is technically difficult , and we do not attempt it here despite its obvious physical interest .
we close with a brief summary and discussion .
the @xmath43 calculation in the @xmath45 case is relegated to an appendix , as are other ( but not all ) technical details .
for an @xmath46 gauge theory with @xmath10 flavors of fermions , we write the lattice hamiltonian @xmath47 here @xmath48 is the electric term , a sum over links @xmath49 of the @xmath46 casimir operator on each link , @xmath50 next is the magnetic term , a sum over plaquettes , @xmath51 finally we have the fermion hamiltonian , @xmath52 the fermion field @xmath53 carries color @xmath54 and flavor @xmath55 at site @xmath56 .
the function @xmath57 is a kernel that defines the lattice fermion derivative .
it can yield a naive , nearest - neighbor action if @xmath58 ; a long - range slac derivative @xcite if @xmath59 ; or anything in between , such as a  action obtained by truncating the slac kernel to its @xmath60 terms .
for @xmath61 the ground state of @xmath62 is determined by @xmath48 alone to be any state with zero electric field , whatever its fermion content , @xmath63    and are degenerate with respect to all the fermionic degrees of freedom .
we consider perturbation theory in @xmath64 .
both @xmath65 and @xmath66 are sums of operators that are strictly bounded , independent of @xmath67 except for the explicit coefficient in @xmath65 .
we can dismiss first - order perturbations by noting that @xmath65 and @xmath66 are multilinear in link operators @xmath68 and @xmath69 , which are raising / lowering operators for the electric field ; thus there are no non - zero matrix elements within the zero - field sector .
we proceed to higher orders , and seek an effective hamiltonian that acts in the zero - field sector @xcite . define @xmath70 to be the projector onto the subspace of all the @xmath71 states .
then perturbation theory in @xmath72 gives an effective hamiltonian , @xmath73 here @xmath74 projects onto the subspace orthogonal to the @xmath71 states ; the operator @xmath75 supplies energy denominators , so that @xmath76    the intermediate states @xmath77 contain flux excitations . in second and third order the patterns of flux can only be strings of length @xmath18 in the fundamental representation of the color group . thus the energy denominators are @xmath78 where @xmath79 is the quadratic casimir of the fundamental representation of @xmath46 .
the perturbations @xmath65 and @xmath66 are explicitly of @xmath24 and @xmath80 , respectively ; each energy denominator gives a factor of @xmath81 .
thus to @xmath82 we can forget about @xmath65 .
our result to this order is @xmath83 since @xmath66 has no non - zero matrix elements within the @xmath71 sector , we have dispensed with @xmath84 in .
the first term in arises for any value of @xmath9 and is @xmath24 ; the case @xmath85 must be treated carefully , but all cases @xmath86 are generic .
the second term is special to @xmath5 and is @xmath82 .
we calculate explicitly the first term in @xmath87 .
each term in @xmath66 creates a string of flux of length @xmath18 , which must be destroyed by the conjugate term .
thus @xmath88\sum_{\bn\mu } \left(\psi_\bn^{\dag\alpha f}\alpha_\mu\psi_{\bn+j\muhat}^{\beta f}\right ) \langle0|\left(\prod u\right)_{\alpha\beta } \left(\prod u^\dag\right)_{\gamma\delta}|0\rangle \left(\psi_{\bn+j\muhat}^{\dag\gamma g}\alpha_\mu\psi_\bn^{\delta g}\right ) , \label{h2eff1}\ ] ] where we define @xmath89 ^ 2 } { \frac12g^2c_f|j|}>0.\ ] ] the matrix element of the gauge fields yields @xmath90 , independent of @xmath18 .    as they appear in ,
each @xmath91 is next to a @xmath92 on a different site .
this invites a fierz transformation on the product of fermion fields , which we write generally as @xmath93 here @xmath94 are combined site , flavor , and color indices , and we have assumed that @xmath95 and @xmath96 are always different while @xmath18 and @xmath95 might be equal [ as in ] .
the matrices @xmath97 are the 16 dirac matrices , normalized to @xmath98 , and we have defined @xmath99 this sign factor is @xmath100 according to whether @xmath97 commutes or anticommutes with @xmath101 ; it will be a constant companion in our calculations . as they appear in @xmath102 ,
the indices @xmath103 are the same site and color but different flavors , and likewise @xmath104 .
leaving the flavor indices explicit , we obtain @xmath105 each fermion bilinear in parentheses is a color singlet located at a given site .
the second term contains the baryon density is zero in the half - filled state , and thus @xmath106 . ] @xmath107 , and the sum @xmath108 is the total baryon number @xmath109 .
we now combine the dirac indices with the flavor indices and write @xmath110 we have defined new matrices @xmath111 as direct products of the @xmath112 dirac matrices and the @xmath12 flavor generators , @xmath113 and we have normalized them conventionally according to @xmath114 the @xmath111 generate a @xmath33 algebra , with @xmath115 .    an alternating flip @xmath116\psi_\bn\ ] ] ( spin diagonalization @xcite ) removes the @xmath101 matrices from the odd-@xmath18 terms in @xmath66 , and hence removes the sign factors @xmath117 from the odd-@xmath18 terms in @xmath102 .
we have finally @xmath118 the odd-@xmath18 terms are of the form @xmath119 which can be written in any basis for the @xmath0 algebra .
the even-@xmath18 terms , however , contain @xmath120 which is defined only in the original basis ( [ dirprod ] ) .
in the zero - field sector in which we work , gauss law constrains the fermion state at each site to be a color singlet .
the drained state @xmath121 , with @xmath122 for all @xmath123 , is the unique state with @xmath124 .
the other color singlet states may be generated by repeated application of the baryon creation operator , @xmath125 with @xmath9 operators @xmath92 .
( here and henceforth , the indices @xmath126 combine the flavor and dirac indices . )
as noted above , at each site @xmath56 the operators @xmath127 generate a @xmath33 algebra , with @xmath128 .
the drained state is obviously a singlet under this algebra .
the creation operator @xmath129 is in the symmetric representation of @xmath33 with one row and @xmath9 columns ( see fig .  [
fig : baryon ] ) .    repeated application of @xmath129 to the drained state gives the state @xmath130 if there are @xmath30 operators @xmath131 , then the state @xmath132 lies in the representation with @xmath9 columns and @xmath30 rows ( see fig .  [
fig : young1 ] ) .
its baryon number is @xmath133 .
the second - order effective hamiltonian @xmath102 preserves @xmath134 , the baryon number on each site .
thus any distribution of @xmath134 defines a sector within which @xmath102 is to be diagonalized . in other words , baryons constitute a fixed background in which to study `` mesonic '' dynamics .
the baryon number at each site fixes the representation of @xmath33 at that site , which is the space of states in which the charges @xmath135 act .
the @xmath136 terms in are of the form @xmath137 , and they commute with the generators @xmath138 of a global @xmath33 symmetry group .
this symmetry is in fact familiar from the lattice hamiltonian of naive , nearest - neighbor fermions : spin diagonalization of @xmath10 naive dirac fermions transforms the hamiltonian into that of @xmath139 staggered fermions . in the weak coupling limit , there are in fact @xmath140 fermion flavors  the doubling problem .
this doubling is partially reflected in the accidental @xmath141 symmetry , which is intact in the @xmath142 limit and is respected by the effective hamiltonian . retaining terms in the fermion hamiltonian ( and thus in @xmath143 ) that involve odd separations @xmath18 does not break this symmetry .    the nielsen - ninomiya theorem @xcite guarantees that any fermion hamiltonian of finite range will possess the full doubling problem .
this is a statement , however , about weak coupling only , since the dispersion relation of free fermions is irrelevant if the coupling is strong and the fermions are confined .
it is interesting that the accidental @xmath0 symmetry nonetheless survives into strong coupling as a vestige of doubling .    the terms in with _ even _ @xmath18 ,
on the other hand , break the @xmath33 symmetry , as do even-@xmath18 terms in the original fermion hamiltonian .
it is easy to see via spin diagonalization , which leaves the even-@xmath18 terms unchanged , that the only generators left unbroken are the @xmath144 corresponding to @xmath145 which form the @xmath146 chiral algebra .
this of course makes no difference to the nielsen - ninomiya theorem , which will enforce 8-fold doubling in the perturbative propagator even without the @xmath0 symmetry .
if we are interested in the realization of the global symmetries of the continuum theory , though , we can study this lattice theory which has the same symmetry .
the simplest theory one may study is thus one containing nearest - neighbor and  terms .
we shall proceed to discard terms with longer range ; we shall begin with the nearest - neighbor theory , with its accidental doubling symmetry , and later break this symmetry to @xmath146 with the  terms .
two essential differences will always remain between this lattice theory and the continuum theory .
one is the presence of the axial @xmath8 symmetry on the lattice .
this symmetry is exact , broken by no anomaly , and may make the drawing of conclusions for the continuum theory less than straightforward unless it is broken by hand .
the other difference is the fact that the effective hamiltonian for baryons ( see below ) is also a short - ranged hopping hamiltonian .
if the baryons were almost free , we would say that they are surely doubled like the original quarks .
the fact that the simplicity of the hopping terms is only apparent , and that the baryons are still coupled strongly to mesonic excitations , offers the possibility that doubling may not return .
the third - order term in @xmath87 , which only exists in the case of @xmath5 , is calculated via @xmath147 for a single link , we have @xmath148 and the same can be proven for a chain of links , @xmath149 thus [ @xmath126 are here ( temporarily ) flavor indices ] , @xmath150 the kernel is @xmath151 again , spin diagonalization simplifies the odd-@xmath18 terms , but not the even-@xmath18 terms .
the result is @xmath152 with @xmath153 where @xmath154 is the usual staggered - fermion sign factor , and @xmath155_{ii ' } b_{\bn+j\muhat}^{i'}\zeta_\mu(\bn)+\text{\em h.c.},\ ] ] where @xmath156 .
the baryon operators are @xmath157 where we have written @xmath158 to represent the compound index @xmath159 , taking values in the symmetric three - index representation of @xmath33 ( @xmath126 once more combine flavor and dirac indices ) . the odd-@xmath18 part of @xmath160 , like that of @xmath143 , is symmetric under the @xmath33 doubling symmetry .
the even-@xmath18 part breaks @xmath33 to @xmath146 .
@xmath160 is a baryon hopping term . as mentioned in the introduction , however , its simplicity is deceptive .
the baryon operators @xmath161 are composite and hence do not obey canonical anticommutation relations , _
i.e. , _ @xmath162 the separation of @xmath160 into a canonical kinetic energy and an interaction term is a challenge for the future .
because of the complexity of the third - order effective hamiltonian , we restrict ourselves henceforth to the second - order theory , in which baryons are a fixed background .
the theory defined by @xmath143 is a generalized spin model , with spins chosen to be in representations of @xmath33 according to the baryon distribution .
we review @xcite in this section how to convert the spin model into a @xmath11 model .
the @xmath11 field at each site will move in a manifold determined by the baryon number at that site .
we employ a generalization of spin coherent states @xcite to derive a path integral for the spin model of @xmath143 .
we recall that a given site carries generators @xmath135 of @xmath33 in a representation with @xmath9 columns and @xmath30 rows , with @xmath133 .
first we choose a basis for the lie algebra of @xmath33 .
this consists of the generators @xmath163 , with @xmath164 , whose matrix elements in the fundamental rep are @xmath165 the corresponding charges are @xmath166 where we have subtracted a constant for convenience .
the cartan subalgebra consists of the operators @xmath167    we build the coherent states from the state of highest weight .
the highest - weight state @xmath168 in the representation is an eigenstate of the cartan generators , @xmath169 -(n_c/2)|\psi_0\rangle&{\rm for}\ i = m+1,\ldots , n .
\label{cartanh}\ ] ] in this state , the generators take the simple form @xmath170 with @xmath171 the state @xmath168 is invariant ( up to a phase ) under the subgroup of @xmath33 that commutes with @xmath172 ; this is @xmath37 .
the most general rotation of @xmath168 is carried out with the generators @xmath173 that are _ not _ in the corresponding subalgebra .
we choose parameters @xmath174 , with @xmath175 $ ] and @xmath176 $ ] , and write @xmath177 the only generators @xmath178 that appear in are those that lower an @xmath179 that starts from @xmath180 in while raising another @xmath179 that starts from @xmath181 .
any other generator would annihilate @xmath168 and thus give no effect in the exponential .
the coherent states are normalized , @xmath182 and over - complete , @xmath183 in the integral is over the coset space @xmath32 $ ] ( see below ) .
matrix elements of the generators are given by @xmath184 where the matrix @xmath185 is given by a unitary rotation from @xmath172 , @xmath186 the matrix @xmath187 is built out of the @xmath188 matrix @xmath174 , @xmath189 .
\label{ua}\ ] ] @xmath11 is both hermitian and unitary .
the manifold of matrices @xmath11 is the coset space @xmath32 $ ] , a sub - manifold of @xmath33 .
this is because for any matrix @xmath187 , one can generate an orbit @xmath190 by multiplying with a matrix @xmath191 where @xmath192 and @xmath193 . all matrices in the orbit
will give the same matrix @xmath11 when inserted into , and thus in integrating over the configuration space of @xmath11 one must choose only a single representative of each orbit .
this set of representatives , the coset space @xmath194 $ ] , is the quotient space of the non - invariant subgroup @xmath37 .
the measure over the coset space must be invariant under unitary rotations , @xmath195 where @xmath196 represents the rotation @xmath72 in hilbert space .
equation  ( [ sq ] ) shows that this is a rotation @xmath197 and by , this means that a measure in @xmath68 must be invariant under @xmath198 .
this fixes the measure uniquely to be the haar measure in @xmath33 , and thus one can integrate over the coset space by integrating with respect to @xmath68 over @xmath33 and using .    a representation whose young diagram has @xmath35 rows is the conjugate to the representation with @xmath30 rows .
its coherent state space can be constructed to look the same , with only a sign difference .
to do this we start with the _ lowest_-weight state , which satisfies [ cf .
] @xmath199 ( n_c/2)|\psi_0\rangle&{\rm for}\ i = m+1,\ldots , n .
\end{array}\right.\ ] ] this introduces a minus sign into .
the subsequent steps are identical , with only the replacement of by @xmath200 here , too , @xmath11 is given in terms of @xmath172 and @xmath68 by .
the partition function @xmath201 can be written as a path integral by inserting the completeness relation ( [ complete ] ) at every slice of imaginary time .
this gives @xmath202 where the action is @xmath203.\ ] ] the hamiltonian @xmath204 is a transcription of the quantum hamiltonian to the classical @xmath11 matrices .
starting with the quantum operator @xmath205 , we have @xmath206 expressed in these variables , the quantum hamiltonian is term from indeed disappears .
we have dropped an additive constant that is independent of @xmath109 . ]
@xmath207 where @xmath208 .
we transcribe this according to to obtain the classical hamiltonian , @xmath209 where @xmath210 recall that each @xmath211 is an @xmath212 matrix ranging over the coset space appropriate to site @xmath56 .
the time - derivative term in @xmath213 is a berry phase @xcite .
it can be expressed in terms of the matrix @xmath68 that determines @xmath11 via .
the result is is of the form given in , and in that case @xmath68 can not be integrated over all of @xmath33 . ]
@xmath214 . \label{action_sigma}\ ] ] @xmath215 will vary from site to site if @xmath30 does . if one takes the route of for a representation with @xmath35 rows , then the kinetic term for that site acquires a minus sign ( see below ) .
the number of colors has largely dropped out of the problem , since @xmath11 is just an @xmath212 unitary matrix field .
the explicit factors of @xmath9 in eqs .
( [ hsigma ] ) and  ( [ action_sigma ] ) invite a semiclassical approximation in the large-@xmath9 limit .
this of course neglects the @xmath9-dependence of the couplings @xmath216 , but the common scale of the couplings only serves to set an energy scale .
the ground state will be independent of this scale , although correlations and the temperature scale will reflect it .
we take the point of view that after all @xmath5 , and we are interested in properties of the effective theory for this value only .
the @xmath217 limit , for fixed couplings @xmath218 , will be a device for investigating the properties of a generalized effective theory .
the simplest way to set the baryon density to zero is just to choose @xmath219 on each site , meaning @xmath220 .
it turns out to be just as easy to consider a slightly generalized case @xcite , in which @xmath221 is chosen to alternate , @xmath222 , on even and odd sublattices .
this means to choose a representation with @xmath223 rows on even sites and @xmath224 rows on odd sites , which gives a pair of conjugate representations of @xmath33 . in view of
, we can substitute @xmath225 on the odd sites and thus have identical manifolds on all sites .
the hamiltonian is then @xmath226 @xmath30 is a new parameter in the theory , and we can ask what value of @xmath30 gives the lowest energy for the ground state .
we will see that @xmath38 is indeed preferred , but only in the next - to - leading order in @xmath2 .
in the large-@xmath9 limit , we seek the classical saddle point of @xmath213 .
we assume the saddle is at a configuration @xmath227 that is independent of time , and so we drop the time derivative .
we begin with the nearest - neighbor hamiltonian , ) is preserved by the alternating signs in the time - derivative term . ] @xmath228 again , the matrices @xmath211 are hermitian and unitary , and the expansion coefficients @xmath229 are real and satisfy @xmath230 .
the minimum of @xmath62 is clearly at a constant field , @xmath231 , which can be diagonalized to @xmath232 by a @xmath33 rotation .
this is a nel state in the original variables .
the @xmath33 symmetry is broken to @xmath37 and there are @xmath233 goldstone bosons .
the classical energy density ( per link ) is @xmath234 , independent of @xmath30 .
thus at leading order in @xmath2 , the optimal value of @xmath30 is undetermined , and any alternating background of baryon number is equally good .    since it turns out that the @xmath2 corrections select @xmath38 , let us consider the effect of the  term in @xmath62 for this case only .
the perturbation is @xmath235 it breaks the @xmath33 symmetry to @xmath236 . assuming that @xmath237 , we again seek the minimum energy configuration in the form of a constant field ; we minimize @xmath238 among the @xmath33-equivalent @xmath239 states that minimize the nearest - neighbor action .
it is not hard to show ( see appendix [ app : nnn ] ) that @xmath240 is minimized for @xmath241 .
this is a condensate that is symmetric under the vector generators @xmath242 but not under the axial generators @xmath243 , and thus it breaks the chiral symmetry to the vector subgroup , @xmath244 .      returning to the nearest - neighbor theory , we consider fluctuations around the @xmath232 minimum of @xmath213 .
first we rescale @xmath245 in order to put the kinetic and potential terms on an equal footing , giving @xmath246 , \label{srescaled}\ ] ] with @xmath247 .
recalling , we can write @xmath248 where @xmath249 is anti - hermitian and anticommutes with @xmath172 .
it is more convenient to work with the hermitian matrix @xmath250 in terms of which we expand @xmath251 if we further expand @xmath252 in the basis of generators of @xmath33 , @xmath253 we find that the @xmath254 corresponding to generators of @xmath37 vanish ; this is the subgroup under which the vacuum is symmetric .
the field @xmath252 thus contains @xmath233 real degrees of freedom , corresponding to the goldstone bosons .
we expand @xmath255 and @xmath211 in powers of @xmath252 ; using , we obtain to second order @xmath256.\ ] ] the coefficient matrix is @xmath257 ) , \label{cetaeta}\ ] ] and the classical energy is @xmath258 .
@xmath259 is antisymmetric and purely imaginary ; we show in appendix [ app : cd ] that @xmath259 has eigenvalues @xmath100 , each with degeneracy @xmath260 .
we change basis so as to diagonalize @xmath259 , and write the index @xmath261 as the compound @xmath262 with @xmath263 and the @xmath264 corresponding to the eigenvalue of @xmath259 .
since the original @xmath254 are real , we have @xmath265 thus we eliminate @xmath266 and write @xmath267 . \label{sln}\ ] ]    the alternating sign in is what makes the theory antiferromagnetic .
it forces us to differentiate between even and odd sites , and we transform to momentum space as follows ( dropping the @xmath268 superscript ) : @xmath269 l^\alpha_{2,\bk}e^{i\bk\cdot(\bn-\hat z)}&\quad \text{$\bn$ odd . } \end{array}\right.\ ] ] the even sites comprise an _ fcc _ lattice with lattice constant 2 , and the momenta @xmath270 take values in its brillouin zone .
we obtain @xmath271 here @xmath272 -jd\gamma(\bk)&jd+\partial_\tau \end{array}}\right)\ ] ] and @xmath273 .
the gaussian path integral over the @xmath96 field now gives the free energy , @xmath274,\ ] ] where @xmath275 and @xmath276 .
for the ground state energy , we take @xmath277 and obtain ( restoring all constants ) @xmath278 .
\label{eq : e0}\ ] ] this is exactly the result of smit @xcite .
the @xmath43 corrections lift the degeneracy of the ground states with different values of @xmath30 .
the integral in is positive and its coefficient contains the number of goldstone bosons .
thus the state of lowest energy is that with @xmath38 , and the symmetry breaking scheme is @xmath40 .
further breaking by the  terms was discussed above .
the zero - density theories considered in the preceding section were defined by selecting representations with @xmath30 and @xmath35 rows on alternating sites . for any @xmath30 , this led to a @xmath11 model with identical degrees of freedom on all sites
after redefinition of the spins on the odd sublattice  and ferromagnetic couplings .
factor in the kinetic energy retained information about the antiferromagnetic nature of the quantum problem ; it did not affect the classical analysis .
] we eventually settled on @xmath38 as the background that gives the ground state of lowest energy .    introducing non - zero baryon density
means changing @xmath30 on some sites of the lattice .
since in general there will be representations on different sites that are not mutually conjugate , different sites will carry @xmath11 variables that do not live in the same submanifold of @xmath33 .
we here limit ourselves to the simpler case of uniform @xmath30 , where adjacent sites carry identical spins  but the coupling is _
antiferromagnetic_.    in order to learn how to work with such a theory ,
we begin by studying the two - site problem . the results of this study will lead directly to an _ ansatz _ for the ground state of a lattice with a fixed density of baryons .
consider , therefore , two sites with quantum spins @xmath279 and  @xmath280 that carry representations of @xmath33 with @xmath281 and  @xmath282 rows , and @xmath9 columns ( see fig .
[ fig : twospins ] ) .
the quantum hamiltonian is @xmath283 an antiferromagnetic coupling .
the corresponding classical @xmath11 model has the interaction hamiltonian @xmath284 where @xmath285 the two @xmath172 matrices reflect the different values of @xmath286 according to @xmath287    the @xmath217 limit is the classical limit , in which we seek values of @xmath288 that minimize @xmath204 . a global @xmath33 rotation , _ viz .
_ , @xmath289 can be used to diagonalize @xmath290 so that @xmath291
. now we have to minimize @xmath292 .
the case of conjugate representations , @xmath293 , is easy : @xmath294 is a unitary rotation of @xmath295 , which ( in this case ) can be rotated into @xmath296 .
this is the unique antiferromagnetic ground state .
@xmath290 and  @xmath294 can be copied to the odd and even sublattices of an infinite lattice to give the classical nel state considered in the preceding section .
the case @xmath297 is more complex .
we consider @xmath42 for definiteness ; the other case is similarly handled .
we write explicitly [ from eqs .
( [ qulu])([ua ] ) ] @xmath298 -\frac{\displaystyle\sin\left(2\sqrt{a^\dag a}\right ) }              { \displaystyle\sqrt{a a^\dag}}\,a^\dag & -\cos\left(2\sqrt{a^\dag a}\right ) \end{array}\right ) .
\label{eq : general_q}\ ] ] @xmath299 is a square matrix of dimension @xmath35 and @xmath300 is a square matrix of dimension @xmath30 . since @xmath294 is a rotation of @xmath172 , @xmath301 and hence the energy is @xmath302.\end{aligned}\ ] ] @xmath303 is minimized when all the eigenvalues of @xmath304 are equal to @xmath305 .
this means that the @xmath35 column vectors @xmath306 form an orthogonal set in @xmath30 dimensions , with @xmath307 since @xmath308 by assumption , such a set of vectors can always be found .    since @xmath309
, we have @xmath310 and so the off - diagonal blocks of vanish .
the lower - right block of @xmath294 is the unit matrix @xmath311 .
we know that @xmath312 since @xmath313 and thus the upper - left block must have eigenvalues @xmath100 . equating traces of @xmath294 and @xmath172
, we find that the upper - left @xmath314 block must take the form @xmath315 where @xmath316 and @xmath317 .
@xmath318 represents the coset @xmath319 $ ] .
we conclude that the classical ground state of this @xmath45 two - site problem is degenerate , even beyond breaking the overall @xmath33 symmetry .
the solutions can be written as @xmath320 \sigma_2&=&\left ( \begin{array}{cc } \sigma^{(m ) }      &        0        \\ 0        &        \bm1_{n - m } \end{array } \right ) \label{eq : q2}\end{aligned}\ ] ] ( to which a global @xmath33 rotation can be applied ) .
a particular instance of @xmath318 is @xmath321 , given by .
the symmetry of the vacuum is the set of rotations that leaves both @xmath290 and @xmath294 invariant , namely , @xmath322 .
the classical solution of the two - site problem will guide us in approaching the problem of an infinite lattice below .
we expect that spontaneous symmetry breaking will give a vacuum of the same character , with continuous degeneracy .
there are , however , two kinds of degeneracy in the two - site problem : that which results from breaking the global @xmath33 to @xmath37 , and that which comes of breaking the @xmath323 subgroup to @xmath324 the latter degeneracy is connected with freedom in choosing the orientation of @xmath294 _ relative to _ @xmath290 .
it is instructive to see how quantum mechanical fluctuations lift the classical degeneracies .
the quantum two - site problem is easy to solve .
we rewrite the hamiltonian ( [ 2site_qh ] ) as @xmath325 . \label{2site_qh2}\ ] ] @xmath326 and @xmath327 are constants , the quadratic casimir operator in the @xmath30-row , @xmath9-column representation of @xmath33 .
the first term in is minimized by coupling @xmath279 and @xmath280 to the representation that minimizes the casimir , which is the representation with @xmath328 rows and @xmath9 columns ( see fig .  [
fig : twosite ] ) .
the ground state has discrete degeneracy equal to the dimension of this representation .
the exact quantum solution naturally shows no sign of spontaneous symmetry breaking and hence it is not of much relevance to the infinite volume problem .
more interesting is the problem where the state of @xmath279 is _ fixed _ and @xmath280 is allowed to vary .
this breaks by hand the global @xmath33 while allowing quantum fluctuations to lift any remaining degeneracy in the relative orientation of the two spins , so it can be regarded as quantization in the presence of spontaneous symmetry breaking . in effect , this is mean field theory .
we replace the hamiltonian ( [ 2site_qh ] ) by @xmath329 to minimize the energy we maximize the mean field by choosing the state of @xmath279 to be the highest - weight state .
this state diagonalizes the generators @xmath179 of the cartan subalgebra while other generators of @xmath33 have expectation value zero .
thus @xmath330 the operators @xmath179 all commute , and their eigenvalues make up the weight diagram of the representation .
traceless diagonal generators of @xmath27 . these can be obtained by isolating the @xmath8 member of the set @xmath179 and taking linear combinations of the rest .
] @xmath331 is a dot product of the weight vectors of the two spins .
the energy is minimized by choosing for @xmath280 a state that lies opposite the highest - weight state in the weight diagram . as shown in the example of fig .
[ fig : twositeweights ] , this still leaves a degeneracy , albeit a discrete one .
we stress that this degeneracy comes from freedom in the relative orientation of @xmath279 and  @xmath280 ; it remains after quantum fluctuations lift the continuous degeneracy of the classical system .    in the @xmath217 limit
, the discrete degeneracy becomes infinite and presumably it is well described by the continuous degeneracy of the classical problem .      at @xmath36
we seek the saddle point of the action , which we assume to be a time - independent configuration .
the classical hamiltonian of the @xmath11 model is @xmath332 seeking an antiferromagnetic ground state , we set @xmath232 on the sublattice of even sites .
the odd sites are then governed by @xmath333 this is just the two - site problem studied above , replicated over the lattice . as we saw above ,
the ground state configuration is degenerate with respect to the configuration at each odd site , @xmath334 a uniform choice for the odd sites , @xmath335 for instance , breaks the @xmath33 symmetry to @xmath322 ; a non - uniform choice can break the symmetry all the way to @xmath336 .
the entropy of this classical ground state is evidently proportional to the volume .
as noted for the two - site problem , the continuous degeneracy of the ground state is an artifact of the classical , @xmath217 limit .
quantum fluctuations will spread each odd spin s wave function over the @xmath319 $ ] manifold . a mean - field _ ansatz _ for the even spins , as noted in the discussion of the two - spin problem , will still leave a discrete degeneracy for the odd spins ; the symmetry breaking scheme will depend on how the odd spins are allocated to the available states .
furthermore , one can contemplate making a non - uniform _ ansatz _ for the _ even _ spins as well , reducing the phase space available for odd spins with unequal neighbors but adding entropy on the even sublattice .
the situation is reminiscent of that in the antiferromagnetic potts model @xcite , the phase structure of which is not yet understood .
this potts - like discrete degeneracy , however , is an artifact of the mean - field approach that , like the classical _ ansatz _ , assumes a fixed state for the spins on the even sites .
an essential difference between our @xmath11 model and the potts model is that our degrees of freedom are continuous and will fluctuate as soon as they are allowed to do so .
a given odd spin will _ not _ be surrounded by a uniform fixed background ; the neighboring even spins will be influenced by _ all _ their odd neighbors , and will induce an interaction among the odd spins that makes them rotate together .
this should reduce the entropy of the ground state to zero . the systematic way to see
this effect is to carry out a @xmath2 expansion around the classical _ ansatz _ , which we do in appendix [ app : corrections ] .
the result is a ferromagnetic interaction among the @xmath211 on the odd sites .
thus the ground state turns out to be the two - site solution , replicated _ uniformly _ over the lattice : @xmath337 \sigma_\bn&=&\left(\begin{array}{cc } \sigma^{(m)}&0\\ 0&\bm1_{n - m } \end{array}\right),&\ \text{$\bn$ odd } , \end{array}\ ] ] where @xmath338 is a _
global _ degree of freedom .
the symmetry group of the vacuum is @xmath322 .
let us summarize the results presented in this paper . in the vacuum sector
, we have rederived smit s result for the lowest - energy configuration of alternating @xmath339 sites . the result is indeed @xmath219 ; the @xmath0 symmetry of the nearest - neighbor theory is spontaneously broken to @xmath340 .
we extended this result to the  theory and found that its @xmath1 chiral symmetry is broken to the vector @xmath12 flavor subgroup . adding net baryon number to the system
, we examined the case of uniform baryon density , @xmath341 . here our study was limited to the nearest - neighbor theory , and we found a nel - like ground state that breaks @xmath33 to @xmath342 .
the number of goldstone bosons @xmath343 thus depends on the baryon density @xmath221 as @xmath344    directions for future work begin with adding  interactions to the @xmath345 theory and extracting from it a prediction for the breaking of the continuum - like chiral symmetry .
( the axial @xmath8 symmetry can be broken by hand . )
another direction is to gain greater freedom in fixing the baryon density .
a constant value of @xmath346 means a baryon density that is close to the maximum allowed on the lattice ; the density can be lowered by setting @xmath347 only on a sparse sublattice , along the lines shown in @xcite . an ultimate goal , as for the hubbard model , is the incorporation of the third - order term in the effective hamiltonian in order to have a theory with dynamical baryons .
perhaps an instructive half - measure would be to study the second - order theory in the presence of a disordered baryon background .
the strong - coupling effective theory can be regarded as a qcd - like model , possessing gauge invariance and the correct degrees of freedom . in that case
the lattice spacing is merely a parameter , an overall scale .
more insight can be gained by considering the strong - coupling theory to represent qcd at large distances , derived by some renormalization - group transformation from a weak - coupling short - distance hamiltonian . on the one hand
, one would expect any such effective hamiltonian to contain many terms of great complexity ; on the other hand , a simple lattice theory such as ours might offer a qualitative approximation to the real theory ( as long as one accepts the loss of lorentz invariance ) .
we can estimate the lattice spacing to be some scale at which the running qcd coupling is large , certainly greater than the radius of a proton .
the limitation that the lattice puts on the density then becomes a physical issue .
taking the lattice spacing @xmath348 to be on the order of 1  fm , the highest baryon density allowed by the lattice is @xmath349  @xmath350 .
for low values of @xmath10 this may not be enough to see finite - density phase transitions , in particular a transition to color superconductivity .
perhaps a way out is to consider an unphysically large number of flavors .
this work was begun while b.s .
was on sabbatical leave at the center for theoretical physics at mit .
he thanks the members of the center for their hospitality .
thanks the organizers of the 19th jerusalem winter school in theoretical physics for their hospitality .
we thank assa auerbach , richard brower , krishna rajagopal , subir sachdev , and uwe - jens wiese for helpful discussions .
this work was supported by the u.s .
department of energy under grant no .
df - fc02 - 94er40818 and by the israel science foundation under grant no .  222/02 - 1 .
the @xmath351 signs are defined only when the @xmath352 are written in the basis @xmath353 .
we choose a chiral basis for the gamma matrices , so that @xmath354 , @xmath355 , and @xmath356 .
the energy ( [ epsprime ] ) is a sum of squares , @xmath357 with the constraint @xmath358 .
the coefficients @xmath359 take on the values @xmath360 .
the minimum of @xmath240 occurs when all @xmath352 are zero except those corresponding to @xmath361 , namely , those for @xmath362 and @xmath363 ; the energy is independent of these @xmath352 .
thus the set of solutions can be written in the form @xmath364 recalling that @xmath365 , we have @xmath366 , so @xmath367 .
a chiral rotation @xmath368 , with @xmath369 turns @xmath370 into @xmath371 , which is invariant only under vector transformations generated by @xmath372 .
to analyze the matrices @xmath259 and @xmath373 , given by eqs .
( [ cetaeta ] ) and  ( [ dmatrix ] ) , we begin with the @xmath33 generators @xmath111 that lie outside the subalgebra @xmath374 that commutes with @xmath375 we choose for them a basis @xmath376 and @xmath377 with @xmath378 and @xmath379 , given by @xmath380 since the coset space @xmath381 is a symmetric space , the commutator @xmath382 $ ] lies in @xmath62 ; in order for @xmath383)$ ] to be nonzero , the commutator must have a nonzero component in the cartan subalgebra of @xmath62 .
this is only possible if @xmath384 and @xmath385 .
thus in this basis @xmath259 takes the form @xmath386 diagonalizing @xmath259 gives eigenvalues @xmath100 .
the generators corresponding to the basis that diagonalizes @xmath259 are @xmath387    as for @xmath388 : the only anticommutators among the @xmath389 that do not vanish [ note the bounds on @xmath390 are between @xmath391 and @xmath392 , _ viz . , _
@xmath393 noting that @xmath394 , we find that @xmath373 takes the block - diagonal form @xmath395      \\ & = & \delta_{qq ' } \sum_{\bmm(\bn ) } ( \bm1_{m}-\sigma^{(m)}_\bmm)_{pp'}.\end{aligned}\ ] ] we summarize this by writing @xmath396 where @xmath397 is an @xmath314 matrix given by @xmath398 it is easy to prove that the eigenvalues of @xmath303 range from @xmath399 to @xmath400 .
in particular , @xmath388 is positive .
we build on the @xmath36 vacua described in sec .  [
sec : nonzerob ] by allowing fluctuations around them .
we let the @xmath211 on the even sites fluctuate around @xmath172 ; we let the @xmath211 on odd sites roll freely around the @xmath319 $ ] manifold covered by , and also execute small oscillations off the manifold into the @xmath32 $ ] coset space .
our goal is an effective action for the classical part ( [ submanifold ] ) of the odd spins . to reach this ,
we integrate out the even spins ; the off - manifold fluctuations of the odd spins must be included for consistency in the @xmath2 expansion .
the counterpart of the action ( [ srescaled ] ) for our problem has an antiferromagnetic spin - spin interaction , with no @xmath401 factors .
we separate it into odd , even , and coupled terms , @xmath402 here @xmath403 , where @xmath404 are the nearest neighbors of the even site @xmath56 .
we expand the field on the even sites around @xmath232 in the manner of , @xmath405 while for the odd sites we write ( see appendix [ app : oddspins ] ) @xmath406 with @xmath407 and @xmath408 @xmath252 describes the fluctuations of the even spins around their classical value @xmath172 .
@xmath255 rotates the odd spins within the manifold of their classical values , while @xmath409 describes their fluctuations outside that manifold .
we further define @xmath410 the classical field on the odd sites .
we leave @xmath411 alone and expand @xmath412 and  @xmath413 around the classical values of the fields , @xmath414 here @xmath415 is the rotated fluctuation field on the odd sites , and the hermitian matrices @xmath416 and @xmath417 are sums over the odd neighbors of the even site @xmath56 , @xmath418 since both @xmath416 and @xmath172 are block diagonal , the first trace in each integral in is zero .
now we organize the partition function as follows : @xmath419 z_{\text{even } } , \label{eq : z}\ ] ] with @xmath420.\end{aligned}\ ] ] @xmath421 is a product of decoupled single - site integrals .
again we expand in the group algebra , @xmath422 where the sum is over the @xmath233 generators of @xmath33 that are not in @xmath37 .
@xmath417 can be expanded similarly and we obtain the following form for the integral over the even fields : @xmath423,\ ] ] where @xmath424 the matrix @xmath259 is the same as in , @xmath425),\ ] ] while the new matrix @xmath388 varies with the site @xmath56 according to the average @xmath416 of its neighboring spins , @xmath426 we study the two matrices in appendix [ app : cd ] .
diagonalizing @xmath259 as before , we arrive at @xmath427 \right\},\ ] ] where @xmath428 is the @xmath314 matrix @xmath429 the quantities @xmath430 ( and @xmath431 ) for each @xmath432 are complex @xmath30-component vectors ; they are rotations of the @xmath233 real components @xmath254 ( and @xmath433 ) into the basis that diagonalizes @xmath259 .
the matrix @xmath397 is given by ; it carries the dependence on @xmath416 . performing the gaussian integration we get @xmath434 .
\label{eq : zeven}\ ] ]    finally we separate the integral ( [ eq : z ] ) over the odd spins into an integral over the classical field @xmath435 and an integral over the fluctuations around it .
we obtain @xmath436 \right\}. \label{eq : z1}\end{aligned}\ ] ] here @xmath437 stands for an odd site , @xmath56 for an even one .
equation  ( [ eq : z1 ] ) gives an effective action for the classical odd spins @xmath438 .
these enter the exponents through @xmath439 [ via eqs .
( [ eq : e ] ) and  ( [ eq : m ] ) ] and through @xmath440 [ via eqs .
( [ tildel ] ) and  ( [ eq : lo ] ) ] .
the action in the first exponent is minimized when each matrix @xmath441 has the largest number of zero eigenvalues , each of which makes @xmath442 approach @xmath443 .
it is easy to check that @xmath397 has @xmath328 zero eigenvalues ( the maximal number ) when the @xmath438 on _ all _ the odd sites @xmath404 align , _
\label{eq : minimum}\ ] ] moreover , when holds , all the @xmath445 s align parallel to each other and @xmath446 is maximized ; also the eigenvalues of @xmath447 are maximized ( to @xmath448 ) .
thus the action in the second exponent also has a minimum at this point in configuration space .
these effects add up to an effective _ ferromagnetic _ interaction among the @xmath449 nearest neighbors @xmath437 of each even site @xmath56 .
this effective interaction will align the classical spins on the odd sublattice to the same direction in their submanifold , @xmath319 $ ] .
the divergences in the effective action have their origin in the fact that the semiclassical corrections are calculated as gaussian integrals in the even fluctuation fields @xmath252 , and the coefficient matrix @xmath450 acquires zero eigenvalues .
the correct range of integration over @xmath252 is of course not infinite , but rather the volume of the @xmath32 $ ] manifold .
this will regulate the divergences , but leave the effective action for the odd spins attractive .
in the classical analysis , the fields on the odd sites take values in the sub - manifold @xmath319 $ ] of the manifold @xmath32 $ ] .
we denote these values @xmath451 , @xmath452 here @xmath453 with @xmath454 , and @xmath455    @xmath318 contains @xmath456 independent degrees of freedom .
any @xmath457 $ ] can be written as @xmath458 -\frac{\displaystyle\sin\left(2\sqrt{a^{\dag}a}\right ) }      { \displaystyle\sqrt{a^{\dag}a}}\,a^{\dag }    & -\cos\left(2\sqrt{a^{\dag}a}\right ) \end{array } \right ) \label{eq : sig_gen}\ ] ] [ cf .  ] , which coincides with if @xmath459 recall that @xmath348 is an @xmath188 matrix , so the zero block has dimensions @xmath460 .
we allow motion out of the sub - manifold by allowing @xmath348 to vary further , @xmath461 the @xmath462 degrees of freedom in @xmath463 complement the @xmath456 degrees of freedom inherent in @xmath464 to give @xmath233 , the dimensionality of the entire @xmath194 $ ] coset space .
writing @xmath11 with the generalized @xmath348 we have @xmath465 0    & \cos\left(2\sqrt{\bar{a}\bar{a}^{\dag}}\right ) & -\bar{a}\,\frac{\displaystyle\sin\left(2\sqrt{\bar{a}^{\dag}\bar{a}}\right ) }      { \displaystyle\sqrt{\bar{a}^{\dag}\bar{a } } } \\[2pt ] 0 & -\frac{\displaystyle\sin\left(2\sqrt{\bar{a}^{\dag}\bar{a}}\right ) } { \displaystyle\sqrt{\bar{a}^{\dag}\bar{a}}}\,\bar a^\dag     & -\cos\left(2\sqrt{\bar{a}^{\dag}\bar{a}}\right ) \end{array } \right )      u^{\dag},\ ] ] with @xmath466 we can also write this as @xmath467}(\bar{a},\bar{a}^{\dag } ) \end{array } \right ) u^{\dag}.\ ] ] @xmath468}$ ] is a matrix in the manifold @xmath469 $ ] . indeed for @xmath470
, we have @xmath471 } = \left ( \begin{array}{cc }      -\bm1_{n - m } &    0    \\      0    &    \bm1_{n - m }   \end{array } \right ) \equiv -\bar{\lambda}.\ ] ] since @xmath469 $ ] is a self - conjugate manifold , its structure near @xmath468}=-\bar\lambda$ ] is the same as its structure near @xmath468}=\bar\lambda$ ] , which corresponds to @xmath472 .
expanding @xmath468}$ ] around @xmath473 gives .
j.  c.  collins and m.  j.  perry , phys .
lett .  * 34 * , 1353 ( 1975 ) .
b.  c.  barrois , nucl .
b * 129 * , 390 ( 1977 ) .
d.  bailin and a.  love , nucl .
b * 205 * , 119 ( 1982 ) ; phys .
b * 107 * , 377 ( 1981 ) ; phys .
rept .  * 107 * , 325 ( 1984 ) .
m.  alford , k.  rajagopal , and f.  wilczek , phys .
b * 422 * , 247 ( 1998 ) .
r.  rapp , t.  schfer , e.  v.  shuryak , and m.  velkovsky , phys .
* 81 * , 53 ( 1998 ) .
d.  t.  son , phys .  rev .
d * 59 * , 094019 ( 1999 ) .
t.  schfer and f.  wilczek , phys .
b * 450 * , 325 ( 1999 ) [ arxiv : hep - ph/9810509 ] .
m.  alford , k.  rajagopal and f.  wilczek , nucl .
b * 537 * , 443 ( 1999 ) .
m.  g.  alford , j.  a.  bowers , and k.  rajagopal , phys .
d * 63 * , 074016 ( 2001 ) [ arxiv : hep - ph/0008208 ] ; j.  a.  bowers and k.  rajagopal , phys .
d * 66 * , 065002 ( 2002 ) [ arxiv : hep - ph/0204079 ] .
k.  rajagopal and f.  wilczek , arxiv : hep - ph/0011333 . k.  rajagopal and e.  shuster , phys .
d * 62 * , 085007 ( 2000 ) [ arxiv : hep - ph/0004074 ] .
z.  fodor and s.  d.  katz , jhep * 0203 * , 014 ( 2002 ) [ arxiv : hep - lat/0106002 ] ; z.  fodor , s.  d.  katz , and k.  k.  szabo , arxiv : hep - lat/0208078 .
de forcrand and o.  philipsen , nucl .
b * 642 * , 290 ( 2002 ) [ arxiv : hep - lat/0205016 ] .
b.  bringoltz and b.  svetitsky , arxiv : hep - lat/0209005 .
t.  banks , l.  susskind , and j.  kogut , phys .
d * 13 * , 1043 ( 1976 ) .
j.  smit , nucl .
b * 175 * , 307 ( 1980 ) .
r.  c.  brower and b.  svetitsky , phys .
d * 61 * , 114511 ( 2000 ) [ arxiv : hep - lat/9912019 ] .
h.  b.  nielsen and m.  ninomiya , nucl .
b * 185 * , 20 ( 1981 ) [ erratum-_ibid .
_  * 195 * , 541 ( 1982 ) ] ; _ ibid . _
* 193 * , 173 ( 1981 ) .
k.  g.  wilson , in _ new phenomena in subnuclear physics _ ,
erice 1975 , edited by a.  zichichi ( plenum , new york , 1977 ) , part a p.  69 .
l.  susskind , phys .
d * 16 * , 3031 ( 1977 ) .
r. narayanan and h. neuberger , phys .
b * 302 * , 62 ( 1993 ) [ arxiv : hep - lat/9212019 ] ; nucl . phys .
* b412 * , 574 ( 1994 ) [ arxiv : hep - lat/9307006 ] ; * b443 * , 305 ( 1995 ) [ arxiv : hep - th/9411108 ] ; h.  neuberger , phys
.  rev .
d * 57 * , 5417 ( 1998 ) .
m.  creutz , i.  horvath , and h.  neuberger , nucl .
suppl .   * 106 * , 760 ( 2002 ) [ arxiv : hep - lat/0110009 ] .
d.  b.  kaplan , phys .
b * 288 * , 342 ( 1992 ) [ arxiv : hep - lat/9206013 ] .
y.  shamir , nucl .
b406 * , 90 ( 1993 ) .
s.  d.  drell , m.  weinstein , and s.  yankielowicz , phys .
d * 14 * , 1627 ( 1976 ) .
s.  myint and c.  rebbi , nucl .
b * 421 * , 241 ( 1994 ) [ arxiv : hep - lat/9401009 ] .
g.  t  hooft , in _ recent developments in gauge theories _ ,
cargse 1979 , ed .
g.  t  hooft _ et al .
_ ( plenum , new york , 1980 ) ; e.  h.  fradkin and s.  h.  shenker , phys .
d * 19 * , 3682 ( 1979 ) ; l.  f.  abbott and e.  farhi , phys .
b * 101 * , 69 ( 1981 ) ; nucl .
b * 189 * , 547 ( 1981 ) .
t.  schfer , nucl .
b * 575 * , 269 ( 2000 ) .
a.  le yaouanc , l.  oliver , o.  pene , j.  c.  raynal , m.  jarfi , and o.  lazrak , phys .
d * 37 * , 3691 ( 1988 ) ; * 37 * , 3702 ( 1988 ) ; * 38 * , 3256 ( 1988 ) .
e.  b.  gregory , s.  h.  guo , h.  kroger , and x.  q.  luo , phys .
d * 62 * , 054508 ( 2000 ) [ arxiv : hep - lat/9912054 ]
. y.  umino , phys .
d * 66 * , 074501 ( 2002 ) [ arxiv : hep - ph/0101144 ] .
m.  e.  peskin , nucl .
b * 175 * , 197 ( 1980 ) .
a.  auerbach , _ interacting electrons and quantum magnetism _
( springer - verlag , 1994 ) .
s.  randjbar - daemi , a.  salam , and j.  strathdee , phys .
b * 48 * , 3190 ( 1993 ) [ arxiv : hep - th/9210145 ] .
t.  kato , prog .
phys .  * 4 * , 154 ( 1949 ) ; a.  messiah , _ quantum mechanics _
( north - holland , amsterdam , 1961 ) , ch .  16 .
b.  schlittgen and u .- j .
wiese , nucl .
suppl .   * 83 * , 718 ( 2000 ) [ arxiv : hep - lat/0005020 ] ; phys .
d * 63 * , 085007 ( 2001 ) [ arxiv : hep - lat/0012014 ] .
a.  n.  berker and l.  p.  kadanoff , j.  phys .
a * 13 * , l259 ( 1980 ) ; j.  r.  banavar , g.  s.  grest , and d.  jasnow , phys .  rev .
b * 25 * , 4639 ( 1982 ) .