according to the instanton liquid model ( ilm ) instantons are largely responsible for the spontaneous breaking of chiral symmetry and the low energy properties of light hadrons @xcite .
this has been demonstrated using phenomenological instanton liquid models .
the two main tasks involved in the construction are    * producing an equilibrium instanton ensemble assuming some interaction between the pseudoparticles .
this is typically done with a variational calculation @xcite . * computing hadronic observables using the above instanton ensemble as the gauge background .
this is a formidable task and the computation necessarily involves several approximations .
it would thus be very desirable to check the physical picture emerging from the ilm , starting from first principles .    in the last few years
there has been a considerable effort in this direction , using lattice simulations of qcd as a tool .
the properties of the instanton ensemble have been extracted from monte carlo generated lattice configurations using several different `` smoothing '' techniques .
the most important parameters , the topological susceptibility , the instanton density , and the average instanton size were found to be in qualitative agreement with the ilm . on the other hand
 although there is a considerable amount of indirect evidence @xcite  so far the lattice has not given any direct results concerning the second part of the above programme .    in ref .
@xcite cooled lattice configurations were used to show that low eigenmodes of the dirac operator saturate the point - to - point hadronic correlators .
more recently a strong correlation was found between the chiral density of low modes of the staggered dirac operator and the topological charge density @xcite .
these results suggest that instantons probably really play an important role in qcd .
all this work however has involved monte carlo generated or cooled lattices that contained other gauge field fluctuations in addition to the instantons .
therefore , based solely on these studies it is impossible to separate the part of the effects due to the instantons and the part ( if any ) that is due to the rest of the gauge field fluctuations .    to my knowledge , the only lattice study so far that makes such a separation possible has appeared in ref .
@xcite . there
, the instanton content of smoothed gauge configurations was reconstructed .
hadron spectroscopy done on these artificial instanton backgrounds was compared to the spectroscopy on the corresponding smoothed configurations .
while the smoothed ensemble yielded physically reasonable results for the pion and the rho mass , the instantons failed to reproduce this .
in fact wilson spectroscopy on the instanton ensemble was hardly different from free field theory ( trivial gauge field background ) . in view of the positive results cited above
, this was highly unexpected , especially that at first sight the only difference between the compared two ensembles appeared to be the presence versus the lack of confinement on the smoothed and the artificial instanton ensembles respectively . according to the ilm , confinement is not essential for the low energy properties of light hadrons @xcite
. the findings of ref .
@xcite indeed suggest that the physics in the chiral limit is rather insensitive to the string tension . with less than half of the original string tension ,
the low fermion modes on the cooled configurations still produced physically sensible hadron correlators .    in the present paper
i address the question of how to reconcile the insensitivity of light hadrons to the string tension ( also expected from the ilm ) with the failure of the instantons alone to reproduce the correct physics ( found in @xcite ) .
i will exclusively consider the su(2 ) gauge group . along the way i will also obtain some useful insight into the structure of the low end of the spectrum of lattice dirac operators and in particular i will show the importance of using a ( close to ) chiral lattice dirac operator .
this is essential for a faithful reproduction of continuum - like quark ( quasi)zero modes .
the plan of the paper is as follows . in section [ se : zmz ]
i briefly discuss the instanton liquid picture of physics in the chiral limit , and in particular the role of the zero modes and the so called `` zero mode zone '' . in section
[ se : cs ] i point out the role of the explicit chiral symmetry breaking of the lattice dirac operator and the importance of minimising it .
i show that the optimised clover dirac operator introduced in @xcite performs much better in this respect than the wilson operator .
i explicitly study the zero mode zone in the simplest nontrivial case , that of an instanton antiinstanton pair and find that the mixing of zero modes closely follows the pattern expected in the continuum .    in section [ se : spectrum ]
i find that the optimised clover operator produces spectroscopy results markedly different from the wilson operator but it still falls short of reproducing the correct physics on the instanton backgrounds .
instead , it yields a peculiar mixture of qcd and free field theory .
i demonstrate that this is due to the presence of a large number of low lying free field modes in the spectrum of the dirac operator with the instanton backgrounds that produce artificially light hadrons .
the free modes are completely absent from the smoothed configurations .
they can be also eliminated from the instanton configurations by superimposing an ensemble of zero momentum gauge field fluctuations ( torons ) over the instantons .
these gauge field fluctuations survive even deep in the perturbative regime .
i argue that they are probably needed to compensate for the anomalously large small - volume effects present in the instanton ensemble due to the lack of confinement . by comparing the density of eigenvalues around zero on the physical ( smoothed ) , the toron , and toron plus instanton ensembles ,
i obtain direct evidence that the instantons are responsible for chiral symmetry breaking . after the elimination of the free quark modes ,
the instanton backgrounds are shown to reproduce the correct physical pion and rho correlators .
section [ se : con ] contains my conclusions along with some speculations on still unanswered questions .
finally , in the appendix , for reference i collected some general properties of wilson type lattice dirac operators that i used in the main text .
some of these can be also found in the literature .
in this section i describe how instantons can produce low lying modes of the dirac operator and why this is important for the physics in the chiral limit . for simplicity
i start with the continuum theory and in the second part of this section i discuss the complications that appear on the lattice .
let @xmath0 be the continuum massless dirac operator which implicitly depends on a gauge field background .
the basic building block of the physics in the chiral limit is the @xmath1 limit of the massive quark propagator , @xmath2 . from the spectral decomposition in terms of eigenmodes of @xmath0 @xmath3
it is easily seen that in the chiral limit the most important modes are the ones with eigenvalues @xmath4 close to 0 .
the crucial assumption of the ilm is that it is the instantons that are responsible for generating the bulk of these lowest eigenmodes .
this can be qualitatively understood as follows .
according to the atiyah - singer index theorem , in the background of an instanton ( antiinstanton ) @xmath0 has at least one negative ( positive ) chirality zero mode . in the presence of an infinitely separated instanton antiinstanton pair
we still expect to find opposite chirality zero modes localised on the two objects .
if the members of the pair are brought closer to one another , the two degenerate zero eigenvalues will in general split into two complex ones still close to the origin .
a very simplistic but qualitatively correct description of this can be given as follows .
let @xmath5 and @xmath6 be the zero modes of the dirac operators @xmath7 and @xmath8 respectively . since we consider @xmath0 in different gauge backgrounds , we explicitly indicate this ; @xmath9 and @xmath10 refer to an instanton and an antiinstanton .
we now want to describe the spectrum of @xmath11 where @xmath12 is some superposition of the gauge fields @xmath9 and @xmath10 . since @xmath5 and @xmath6 are of opposite chirality , @xmath13 .
let us arbitrarily complete the set of these two vectors into an orthonormal basis and construct the matrix of @xmath11 in this basis .
the elements of this matrix in the subspace spanned by @xmath5 and @xmath6 are @xmath14 where @xmath15 the diagonal matrix elements vanish because @xmath0 maps left handed vectors into right handed ones and vice versa . by a suitable choice of the phases ,
the off - diagonal elements can always be made real and due to the anti - hermiticity of @xmath0 they are of equal magnitude and have opposite signs .
if we now assume that @xmath16 and @xmath17 lie approximately in this two - dimensional subspace , i.e. the matrix of @xmath11 contains the above @xmath18 block - diagonal part , then @xmath11 can be easily seen to have two complex conjugate eigenvalues @xmath19 with the corresponding ( unnormalised ) eigenvectors being @xmath20 .
this mechanism can be generalised to gauge fields which are superpositions of several instantons and antiinstantons .
the basis in this case consists of the zero modes of the individual ( anti)instantons . the only additional approximation involved is the assumption that the zero modes corresponding to different ( anti)instantons are orthogonal .
if the pseudoparticles are well separated , this is approximately true .
one can thus consider the matrix of @xmath0 restricted to the subspace of zero modes  called the `` zero mode zone ''  in a general instanton background .
if the off diagonal elements of this matrix are small then it will have complex eigenvalues close to the origin , in addition to a number of exact zero eigenvalues corresponding to the total topological charge .
moreover , the complex eigenmodes will be approximately linear combinations of several instanton and antiinstanton zero modes and thus become highly delocalised , making it possible to propagate quarks to large distances .    the most important quantities in this construction
are the off diagonal instanton antiinstanton matrix elements , the @xmath21 s . in the simplest case of only one pair
, @xmath21 is known to depend on the distance and the relative orientation of the members of the pair as @xmath22 where @xmath23 describes the relative orientation in group space , @xmath24 and @xmath25 is the relative position of the instanton and the antiinstanton @xcite .
@xmath26 is a function of the instanton and antiinstanton scale parameters and the distance between them . at small distances
it depends on the ansatz used to combine the field of the instanton and the antiinstanton , at large separation it is expected to fall off as @xmath27 .
we have seen how instanton zero modes can give rise to the low - lying eigenmodes that are the most important ones in the chiral limit .
close to the chiral limit the number of modes giving a substantial contribution to the quark propagator becomes smaller and smaller .
it is thus crucial to reproduce these low eigenmodes and eigenvalues as faithfully as possible on the lattice if we want to test the ilm .
the most important obstruction to this is that due to the nielsen - ninomiya theorem it is impossible to construct a local lattice dirac operator describing one fermion species with exact chiral symmetry @xcite . in the present paper
i shall limit the discussion to wilson type lattice dirac operators that describe one fermion species but explicitly break chiral symmetry . in the absence of chiral symmetry
the fermion zero modes corresponding to the topology of the gauge field are not protected ; any small fluctuation of the gauge field can shift the zero eigenvalues .    in order to understand how this happens , a useful starting point is the following property shared by the lattice and the continuum dirac operator .
let @xmath28 be an eigenvector of @xmath0 corresponding to the eigenvalue @xmath29 . if @xmath0 has no degeneracies then
@xmath30 if and only if @xmath29 is real .
this is a simple consequence of the so - called `` @xmath31 hermiticity '' , @xmath32 , a common property of the continuum and the lattice dirac operator .
a proof of the above statement is given in the appendix .    in the continuum , @xmath0 is anti - hermitian and has eigenvalues on the imaginary axis .
thus the only real eigenvalue it can have is 0 .
chirality makes a sharp distinction between zero and non - zero modes in the continuum .
zero modes have a chirality of @xmath33 while all the non - zero eigenmodes have zero chirality .
since the eigenmodes and their chiralities depend continuously on the gauge field , zero modes are protected , they can not be shifted away from zero by continuous deformations of the gauge field .
on the other hand , the lattice dirac operator is not anti - hermitian , and its eigenvalues are not constrained to be on the imaginary axis . while lattice zero modes are still protected from being shifted off the real axis by smooth deformations of the gauge field ( this would cause their chirality to jump to zero ) ,
there is nothing preventing them from moving continuously along the real axis .
indeed , a lattice dirac operator will in general not have any exact zero modes unless the gauge field is fine tuned .
it is now clear that the lattice analogues of the continuum zero modes are real modes of the dirac operator . these have non - vanishing chiralities but in general their magnitude is smaller than 1 .
the crucial role chiral symmetry plays in this discussion can be understood by noting that chiral symmetry of @xmath0 ( i.e. @xmath34 ) together with @xmath31 hermiticity would imply that @xmath0 is anti - hermitian and has protected zero modes .
the smaller the explicit chiral symmetry breaking of @xmath0 is , the more its low lying spectrum will resemble that of an anti - hermitian operator . recently it has been shown that the closest a lattice dirac operator can come to being chirally symmetric is by obeying the ginsparg - wilson relation and thus having an ultralocal chirality breaking in the quark propagator @xcite . in this case
the spectrum of @xmath0 lies on a circle in the complex plane passing through the origin .
indeed the low end of this spectrum is almost on the imaginary axis and there are protected zero modes .
i shall now discuss how the continuum description of the zero mode zone presented in the previous subsection has to be modified on the lattice .
assume that there is an instanton on the lattice .
in the limit when the instanton size goes to infinity ( in units of the lattice spacing ) , the instanton `` does not see '' the lattice at all , and the continuum description applies .
the lattice dirac operator , @xmath7 has an exact zero eigenvalue with a corresponding @xmath35 chirality eigenmode @xmath5 . if the instanton is not infinitely large , the real eigenvalue gets shifted away from zero , its magnitude depends on the instanton size and the chirality breaking of @xmath7 @xcite .
the corresponding eigenmode also becomes only an approximate eigenvector of @xmath31 and thus @xmath36 .
let us now consider the dirac operator in the presence of an instanton antiinstanton pair .
following the continuum discussion in the previous subsection , its matrix elements between the ( lattice approximate ) `` zero modes '' @xmath5 and @xmath6 can be parametrised as @xmath37 where @xmath38 is real , @xmath39 and @xmath40 are non - negative real numbers . in the su(2 ) case the diagonal elements are automatically real and the off - diagonal ones can be also made real by a suitable choice of the phases of the eigenvectors .
this will be proved in the appendix . due to the @xmath31 hermiticity of @xmath0 , since @xmath5 and @xmath6 are approximate eigenvectors of @xmath31 , the off - diagonal terms have opposite signs and about the same magnitude .
the eigenvalues of this matrix are @xmath41 with the corresponding ( unnormalised ) eigenvectors being @xmath42 the most important feature of these eigenvectors is that  as in the continuum  they are mixtures of the instanton and antiinstanton zero modes with roughly the same magnitude ( the ratio being @xmath43 ) and with a relative phase @xmath44 , up to a correction proportional to @xmath38 .
this mechanism produces highly delocalised eigenmodes in the zero mode zone and thus facilitates the propagation of quarks by jumping from instanton to antiinstanton .
maintaining this feature on the lattice is thus crucial for the description of the zero mode zone .
the requirement for the proper mixing of zero modes is the inequality @xmath45 for producing qualitatively continuum - like eigenmodes , the difference of the diagonal elements has to be much smaller than the magnitude of the off - diagonal mixing matrix elements of @xmath11 .
otherwise the two eigenvalues will be real and the corresponding eigenmodes rapidly become localised on the instanton and the antiinstanton respectively . recall that in the continuum ,
the mixing matrix elements @xmath46 are proportional to @xmath47 , where @xmath48 is the invariant angle of the instanton antiinstanton relative orientation ( see eq .
( [ eq : relor ] ) ) .
@xmath38 ( and @xmath49 ) is expected to be of the order of the real eigenvalues of @xmath7 and @xmath8 , both small but nonzero numbers .
this means that inequality ( [ eq : ineq ] ) will be violated for @xmath48 sufficiently close to @xmath50 .
we can minimise the range in relative orientation , where this happens by making @xmath38 small , i.e. having the would - be zero modes close to zero .
in this section i study how the general features of the zero mode zone described in the previous section are realised by two commonly used lattice dirac operators , the wilson and the clover operator .
the gauge background i use for this test contains an instanton of size @xmath51 and an antiinstanton of size @xmath52 , separated by a distance @xmath53 in the time direction .
the lattice size is @xmath54 and i always use antiperiodic boundary conditions in the time direction , periodic in all other directions .
these parameters are typical in the spectroscopy calculations that i shall present in the next section .
let me start with the wilson operator . in the presence of only the ( anti)instanton with the above parameters
, @xmath55 has a real eigenvalue at ( 0.17 ) 0.51 with chirality ( 0.39 ) -0.53 .
it is clear that for these small instantons the wilson operator does not even come close to satisfying the conditions of the previous section : the `` zero modes '' are far away from zero with chiralities of magnitude much smaller than 1 .
the simplified description of the instanton - antiinstanton configuration in terms of the @xmath18 matrix is then inadequate .
if the instanton and antiinstanton are close to being parallel @xmath56 then the real part of the ( anti)instanton related eigenvalues is around 0.3 .
moreover already at @xmath57 the two complex eigenvalues become real and the corresponding modes are rather localised at the instanton and the antiinstanton respectively .
a much better description of the zero mode zone is provided by the optimised clover action proposed in ref .
( see also @xcite for more discussion and a recent scaling test of this action . )
the clover coefficient was chosen to minimise the range in which the real eigenvalues occur in the physical branch of the spectrum on locally smooth gauge backgrounds . as explained in the previous section ,
this is exactly what is needed for a good description of the mixing of zero modes .
the optimal value for @xmath58 was found to be around 1.2 on a set of ape smeared @xmath59 quenched su(3 ) configurations .
physically these are quite similar to the su(2 ) configurations appearing in the present study , therefore in what follows i shall always choose @xmath60 .
the real modes occurring in the single ( anti)instanton background already reveal that the clover action is much closer to the continuum than the wilson action .
the eigenvalues are located at -0.014 and -0.002 and the chiralities of the corresponding eigenmodes are 0.993 and -0.995 respectively . in the instanton plus antiinstanton background
the lowest two modes are always complex for almost any relative orientation of the pair , even at @xmath61 .
we can hope that the simple description in terms of the mixing matrix works well in this case .
indeed , the magnitude of the diagonal elements is small , neither @xmath38 nor @xmath29 goes above 0.02 in magnitude for any relative orientation . in fig .
[ fig : t12 ] i plotted the two off - diagonal matrix elements , @xmath62 versus @xmath48 along with the one - parameter fits of the form const.@xmath63 . as expected ,
the mixing matrix elements are very well described by the continuum ansatz of eq.([eq : relor ] ) .
a comparison of the eigenvalues of the @xmath18 mixing matrix and the explicitly computed eigenvalues of the lattice dirac operator in the corresponding gauge backgrounds shows agreement to within @xmath64 in the whole range of @xmath65 .
we can conclude that the @xmath60 clover operator should give a satisfactory treatment of the zero mode zone as long as the instanton size does not drop below around @xmath51 .
this is of course true only if the gauge configuration is locally smooth .
otherwise the fermion mass acquires an additive renormalisation which completely destroys the above features important for the description of the zero mode zone .
having a fermion action that is expected to describe the zero modes and their mixing reasonably well , we can discuss the main topic of the paper , namely , what part of qcd is reproduced by instantons alone . in this section
i present quenched hadron spectroscopy results obtained with the @xmath60 clover fermion action .
the gauge field backgrounds on which the spectroscopy is done are the ones used in ref .
the starting point is a set of 28 @xmath54 su(2 ) configurations generated with a fixed point action .
the lattice spacing is @xmath66fm , as fixed by the sommer parameter of the heavy quark potential . these monte carlo generated configurations are not yet suitable for our purposes since they are not locally `` smooth '' .
both the instantons and the details of the zero mode zone are obscured by short distance fluctuations .
cycling , a smoothing technique based on the renormalisation group , has been shown to preserve the main physical features of the gauge configurations .
after 9 cycling steps the string tension @xcite and @xmath67 versus @xmath68 are essentially unchanged .
moreover , these lattices are smooth enough that their instanton content can be unambiguously identified .
in fact , about @xmath69 of the action is accounted for by instantons ( assuming there is no interaction between them ) .
we can create artificial lattice configurations that have the same instanton content as the smoothed ones .
the instanton sizes and locations are reproduced but the relative orientation in group space is distributed according to the su(2 ) haar measure .
recently a study of the relative orientation appeared in @xcite , however reproducing the relative orientations would make our procedure much more cumbersome .
i compare the spectroscopy done on the smoothed ( cycled ) lattices and the corresponding artificial instanton configurations with the same instanton content .
these ensembles are both locally smooth enough that the optimised clover action should give a satisfactory description of their zero mode zone . in fig .
[ fig : pirho ] the vector meson mass is plotted as a function of the pseudoscalar mass for both ensembles .
the correlators were always fitted with the assumption that there is only one lightest particle dominating both the pseudoscalar and the vector channels .
i also show the @xmath70 free field line .
in fact , the wilson fermion action produces exactly degenerate free field like pions and rhos on the instanton backgrounds .
compared to the wilson action , the clover action shows a marked improvement , but it still fails to reproduce the @xmath67 vs.  @xmath68 observed on the smoothed lattices .
[ fig : pi ] shows the pion mass sqared vs.  the bare quark mass for both ensembles . on the smoothed ensemble @xmath71 , as expected from pcac .
this is clearly not the case on the instanton ensemble , where a best fit to the form @xmath72 gives @xmath73 which is between @xmath74 and the free field value , @xmath75 .
we also note that since both ensembles are locally very smooth , the additive mass renormalisations are very close to zero .
therefore , it is also meaningful to compare hadron masses obtained at the same bare quark mass on the two gauge ensembles .
this comparison reveals that the instanton ensemble typically produces much lighter hadrons than the smoothed ensemble .      in spite of its substantially improved chiral properties ,
even the clover dirac operator fails to reproduce the correct physics of qcd on the instanton backgrounds .
instead , it yields a peculiar mixture of qcd and free field theory .
is there a problem with the instanton liquid model or the physics of these instanton lattices is still not close enough to the continuum ? to answer this question we look at the low lying eigenmodes of the dirac operator in more detail .    in fig .
[ fig : dspec ] i plot the lowest 30 eigenvalues in the complex plane for the instanton and the smoothed ensemble .
the eigenvalues of all 28 configurations are superimposed .
the boundary condition is antiperiodic in the time direction , periodic in all other directions . in both cases
most of the eigenvalues lie very close to a circle passing through the origin .
this is a sign of the approximate chiral symmetry of the action .
the most striking difference between the two ensembles is the appearance of a gap at around im@xmath29=0.5 on the instanton ensemble .
it is also instructive to compare the density of eigenvalues . since the low eigenvalues
are almost imaginary i project them on the imaginary axis and plot the densities as a function of the imaginary part of the eigenvalues .
this corresponds to the density of eigenvalues around zero in the continuum .    in fig .
[ fig : ev_hist1 ] we can compare the densities corresponding to the two ensembles .
the real modes have been removed , these would show up as spikes at the origin .
the densities at zero seem to agree quite well . according to the banks - casher relation
the chiral condensate is proportional to the density of modes around zero @xcite ( excluding exact zero modes , which do not contribute in the infinite volume limit ) .
this shows that on both the instanton and the smoothed ensemble chiral symmetry is broken and the value of the chiral condensate is approximately the same .
apart from the vicinity of the origin , however , the two distributions are quite different . on the instanton ensemble
there is a substantial `` piling up '' of modes above im(@xmath76 and then a `` thinning '' of modes at 0.5 .
what is the reason of this huge difference between the two distributions ?
the answer can be easily given by looking at the eigenvectors of the dirac operator .
it turns out that the quark density of the eigenvectors on the instanton ensemble corresponding to the peak of the distribution are much more delocalised than all the other eigenvectors occurring in either ensemble .
in fact , these modes are essentially free quark modes .
this can be seen as follows . on a trivial gauge field configuration ( all links are 1 ) with periodic bondary conditions in all directions , there are @xmath77 trivial ( constant ) zero modes , this is the number of c - number degrees of freedom corresponding to one fermion flavour . in the presence of antiperiodic boundary condition in the time direction , the lowest eigenmodes are shifted away from zero , to @xmath78 since in our case @xmath79 .
both eigenvalues are @xmath77-fold degenerate .
the mixing of the free field modes into a particular eigenmode @xmath28 can be characterised by @xmath80 , where @xmath81 is the projection of the normalised eigenmode @xmath28 onto the eigenspace corresponding to the eigenvalue @xmath82 .
this quantity is not gauge invariant therefore i work in landau gauge . any generic fermion mode on a non - trivial but locally smooth configuration
will have a nonzero projection on the @xmath82 eigenspaces .
the question is whether this is just an accidental mixing or the given ensemble really contains close to free field modes . to decide this ,
a good quantity to look at is @xmath83 .
if the mixing is accidental , we expect @xmath84 to fluctuate around 1 , independently of the corresponding eigenvalue @xmath29 . on the other hand , if there are free field like modes on a given configuration then @xmath84 will increase substantially when @xmath29 approaches the free field eigenvalue @xmath85 .    in fig .
[ fig : proj ] i plotted @xmath84 versus the distance of the given eigenvalue from the free field mode @xmath85 .
the figure shows a random selection of eigenvalues with imaginary parts between 0.0 and 0.2 . in the smoothed ensemble @xmath84
fluctuates around 1 everywhere , there is no trace of the free field modes . on the other hand , in the instanton ensemble
, the modes close to @xmath85 have a substantially larger projection on the @xmath85 eigenspace than on the @xmath86 subspace .
in fact , these modes , close to @xmath85 have @xmath87 , so they are essentially free field modes .
a detailed study of the chiral density @xmath88 reveals that all the modes with @xmath89 look like mixtures of instanton zero modes with the density concentrated in several lumps .
this is to be contrasted with the chiral density of the modes with @xmath90 that spreads almost homogenously over the whole lattice , as expected of the lowest free field eigenmodes .
now we can understand why the instanton ensemble gave substantially smaller hadron masses than the smoothed ensemble .
the reason is that the former had a large number of free quark modes .
even if the zero mode zones of the two ensembles are similar , the free field modes provide a very efficient way to propagate quarks to large distances .
this is why the instanton ensemble shows a peculiar mixture of qcd and free field characteristics . in the smoothed ensemble confinement
completely eliminates the freely propagating quarks .      according to the instanton liquid model , instantons alone , without confinement ,
can reproduce most of the properties of the light hadrons .
our result shows that in the presence of instantons only , without confinement , free quark modes provide a very effective way to propagate quarks and they contaminate the hadron correlators producing unusually light hadrons .    one possible way to suppress the free quark modes
is to add to the instantons some locally very smooth background that contains only long wavelength fluctuations .
there has to be enough long distance structure to suppress the free propagation of quarks but on the other hand they have to be locally smooth enough not to distort the instantons .
cycling is a very efficient way to create such backgrounds . to this end
i started with an ensemble of wilson @xmath91 ( @xmath92fm ) @xmath93 su(2 ) gauge configurations , performed 4 cycling steps and finally inverse blocked them to size @xmath54 .
this produced a set of lattices extremely smooth on the few lattice spacing scale but containing the longest wavelength fluctuations characteristic to a @xmath94@xmath95 lattice .
their spatial size was a bit smaller , their temporal size a bit larger than the confinement scale .
moreover these configurations had a small enough physical size that they almost never contained instantons .
after superimposing these lattices on the instanton configurations by simply multiplying the corresponding links ( in landau gauge ) , i checked that the original topological charge of the instanton configurations changed very little .    to show how the addition of the smooth background affects the hadron correlators , in fig .
[ fig : pirho_b2.4 ] i plotted typical pseudoscalar and vector correlators ( at @xmath96 ) .
the change is dramatic .
while on the instanton configurations the pion and the rho are very light ( crosses ) , the addition of the smooth background brings the correlators ( bursts ) very close to the physical ones obtained on the smoothed lattices ( boxes ) .
the good agreement of the physical correlators and the instanton plus smooth background correlators persists in the whole range of masses where i tested it ( @xmath97 )
.    it would be tempting to interpret this as a consequence of confinement ; the superimposed smooth backgrounds still contained the longest scale fluctuations of the confining lattices and these removed the free field modes that contaminated the spectroscopy before .
but is confinement really needed for that ?
this can be easily tested by replacing the superimposed backgrounds with ones of much smaller physical size .
the lattices that i used for this purpose were generated in exactly the same way as before , except that the starting @xmath98 lattices were produced at wilson @xmath99 .
these lattices have a tiny physical size , they are expected to be completely perturbative . in fig .
[ fig : pirho_b3.0 ] i plotted again the pion and rho correlators at @xmath96 .
quite surprisingly , the correlators have absolutely no dependence on the @xmath100 at which the background was created , they are the same on the instantons plus confining and instantons plus perturbative configurations .    apparently , there is still something nontrivial on these physically very tiny perturbative lattices that has the same effect on the correlators as a confining background . what could this be ?
the answer can be given by noting that the @xmath101 limit of yang - mills theory on the torus is not completely trivial because there is not only one gauge field configuration with zero action .
the set of flat ( minimal action ) gauge field configurations can be characterised by the su(2 ) conjugacy classes represented by four mutually commuting gauge group elements corresponding to the ( constant ) polyakov loops in the four perpendicular directions .
these are essentially zero momentum gauge field configurations , with a constant gauge potential ( link ) , sometimes referred to as torons .
their contribution to physical quantities is a finite size effect , and has been estimated in ref .
@xcite .    to explicitly check that it is really the torons that affect the correlators so strongly , i extracted the `` toron content '' of the perturbative background fields .
this was done in the following way .
after landau gauge fixing i averaged the polyakov loops in all four directions . since the polyakov loops along any given direction were almost constant ( up to small perturbative fluctuations ) , the averages were close to being su(2 ) elements .
i projected the averages back onto su(2 ) and then distributed the average polyakov loop evenly among the links in the given direction .
this resulted in configurations with constant gauge fields ( links in any direction ) that carried the average polyakov loops of the original perturbative configurations .
i shall refer to these lattices as `` toron '' configurations .
we can now compare the correlators in the instantons plus toron background and the instantons plus the full perturbative background . the result
 again in the typical case of @xmath96  is shown in fig .
[ fig : toron ] .
the instantons and the torons together fully reproduce the correlators , demonstrating that the most important fluctuations on the perturbative lattices are indeed the torons .
it is also quite instructive to compare the distribution of the lowest eigenvalues of the dirac operator on the different types of gauge configurations .
[ fig : bc_hist ] shows the densities on the smoothed lattices ( solid lines ) , the perturbative ones ( dotted line ) and the lattices containing both the instantons and the perturbative backgrounds ( dashed line ) .
the three distributions agree quite well away from zero .
however around zero , the perturbative ensemble produces only a very low density ( compatible with zero ) , while the two other ensembles yield qualitatively similar densities , both nonzero .
this demonstrates very clearly that it is the instantons that are responsible for creating a nonzero eigenvalue density around zero which through the banks - casher relation @xcite implies the spontaneous breaking of chiral symmetry .
we also note that the peak corresponding to the free field modes ( see fig .  [ fig : ev_hist1 ] ) has been completely eliminated by the perturbative backgrounds .
this again can be attributed to the torons , as can be checked by looking at the eigenvalue distribution on the toron lattices ( not shown ) which is very similar to the dotted line in fig .
[ fig : bc_hist ] .
in the present paper i addressed the question of how much of qcd in the chiral limit is reproduced by restricting the fluctuations of the gauge field to instantons only .
i reconstructed the instanton content of smoothed monte carlo lattice configurations and compared hadron spectroscopy on this instanton ensemble to the spectroscopy on the original `` physical '' smoothed configurations .
the fermion action used for the spectroscopy was a `` chirally optimised '' clover action .
i explicitly studied the fermion zero modes in the presence of an instanton and also the mixing of the zero modes corresponding to an instanton and an antiinstanton .
i found that in these simple cases the optimised action described the continuum features of the zero mode zone rather well . in spite of this , the optimised action still failed to reproduce the physical hadron spectrum .
it yielded anomalously light hadrons with a mixture of qcd and free - field like features .    a closer look at the spectrum of the dirac operator revealed that on the instanton ensemble it had a large number of free quark modes .
these were absent on the physical smoothed configurations .
they provided a very efficient way of propagating quarks to large distances thereby substantially reduced the hadron masses and contaminated the spectroscopy .
superimposing an ensemble of very smooth , essentially perturbative gauge field fluctuations on the instantons , was enough to eliminate the free quark modes and to restore the physical hadron correlators . by comparing the density of low eigenvalues of the dirac operator
i also obtained direct evidence that it is the instantons that are responsible for creating a nonzero density of modes around zero and by the banks - casher relation also for chiral symmetry breaking .
the important fluctuations on the superimposed perturbative lattices , needed to recover the physical hadron correlators , turned out to be zero momentum ( constant link or gauge potential ) configurations , torons .
i would like to emphasise here that this toron ensemble was a perturbative one , produced at large @xmath100 i.e. small physical lattice size .
the meson correlators appeared to be largely independent of the @xmath100 at which the superimposed background was created , as long as it was above the finite temperature phase transition on the given lattice size .
it seems rather surprising that such `` mild '' gauge field fluctuations as the torons , can have so profound an impact on physical quantities .
after all , torons exist only on the torus , they are absent in an infinite space - time and thus constitute only finite size effects . on the other hand we know that the instantons by themselves do not confine @xcite .
this means that the instanton ensemble does not produce a mass gap and it is expected to exhibit much stronger finite size effects than qcd , from which the instantons have been extracted . in fact , the free field modes that are eliminated by the torons are also absent in the infinite volume limit , since they are not normalisable .
in free field theory the density of modes around zero is proportional to the linear size of the box @xmath102 , whereas the density of instanton related modes grows proportionally to the volume .
it is thus not inconceivable that in increasingly bigger volumes , although the torons have less and less influence , at the same time they become less and less needed since the relative importance of the free field modes also dies out .
the torons might be needed only for compensating the unnaturally large finite size effects due to the absence of confinement in the instanton ensemble .    what is the picture emerging from this lattice study concerning the instanton liquid model ?
the qcd vacuum contains an ensemble of instantons that is generated by the non - perturbative gauge field dynamics .
the instantons by themselves break chiral symmetry but hadron correlators in the instanton backgrounds are strongly contaminated by freely propagating quarks .
this yields anomalously light mesons and a small splitting between the pseudoscalar and the vector channel .
the free quarks can be eliminated and the physical hadron correlators can be restored by superimposing a perturbative ensemble of zero momentum gauge configurations ( torons ) on the instantons .
this suggests that the contamination by free quarks is most likely a finite volume effect which , in qcd , is completely suppressed by confinement .    more evidence in favour of this scenario could be gathered by studying the volume dependence of these effects on larger volumes .
another important question , not answered by the present study is why it is exactly the perturbative toron ensemble that is needed to restore the physical meson correlators .
it would be also very desirable to extend this investigation to a larger class of observables , in particular baryon correlators , for which the su(3 ) case need be considered .
finally , for more precise tests a better chiral action , such as the overlap @xcite could be used .
in the appendix , for reference , i collected some properties of wilson type lattice dirac operators that were used in the main body of the paper .
some of these results can be found in the literature , in particular in refs .
@xcite .
i assume that the dirac operator satisfies @xmath31 hermiticity , i.e. @xmath103    let @xmath0 be a @xmath31-hermitian dirac operator with no degeneracy in its spectrum , i.e.  having as many different eigenvalues as the dimension of the space it acts on .
let @xmath29 , and @xmath49 be two eigenvalues of @xmath0 with the corresponding eigenvectors being @xmath28 and @xmath104 .
then @xmath105 if and only if @xmath106 .
_ proof : _ a simple consequence of @xmath31 hermiticity is that @xmath107 which implies
@xmath108 it immediately follows that if @xmath109 then @xmath106 . to prove the converse we first note that since @xmath0 has as many different eigenvalues as the dimension of the space it acts on , its eigenvectors form a basis .
we have already seen that the only eigenvector on which @xmath110 can have a nonzero projection , is @xmath111 .
if this projection were also zero , then @xmath110 would be orthogonal to all the vectors in a complete set and consequently it would be zero .
this is impossible since @xmath31 has a trivial kernel .    for proving that @xmath112 it was essential that the eigenvalues of @xmath0 be non - degenerate and its eigenvectors form a complete set .
generally this is expected to be the case unless the gauge field is fine tuned .
there is however an important exception . when topology
 as seen by the fermions  changes through a smooth deformation of the gauge fields , two opposite chirality real eigenvalues collide and leave the real axis as a complex conjugate pair .
it can be shown that when the two eigenvalues coincide , the corresponding eigenspace is one - dimensional , the eigenvector has zero chirality , @xmath0 does not possess a complete set of eigenvectors and therefore has no spectral decomposition .
a simple consequence of statement 1 is that if @xmath0 has a complete non - degenerate spectrum , then the chirality of an eigenvector , @xmath88 is nonzero if and only if the corresponding eigenvalue , @xmath29 , is real .    in the remainder of the appendix
i shall prove the following property of wilson type su(2 ) lattice dirac operators .
let @xmath113 and @xmath114 be real eigenvectors of two ( not necessarily the same ) su(2 ) lattice dirac operators such that the corresponding eigenvalues and their complex conjugates are non - degenerate .
it is then possible to choose the phases of @xmath113 and @xmath114 such @xmath115 , @xmath116 , @xmath117 , and @xmath118 are all real for any ( third ) lattice dirac operator @xmath0 .
_ proof : _ besides @xmath31 hermiticity , this property depends on an additional symmetry of the dirac operator relating @xmath0 to its complex conjugate @xcite .
this is specific to the su(2 ) case .
let @xmath119 be the charge conjugation operator and @xmath120 , where @xmath121 is a pauli matrix acting in colour space . using that @xmath122 ( where the last equality follows because we work in a representation with all @xmath123 s hermitian ) , and that for any @xmath124 @xmath125 , it is not hard to prove that @xmath126 now if @xmath28 is an eigenvector of @xmath0 with eigenvalue @xmath29 then @xmath127 which means that @xmath128 is an eigenvector of @xmath0 with eigenvalue @xmath129 .
due to the non - degeneracy of the eigenvalues and that @xmath130 , it follows that @xmath131 in the special case when @xmath29 is real , the phase choice @xmath132 eliminates the extra phase factor on the r.h.s .
( [ eq : phase ] ) .    now writing @xmath133 , and making use of eqs .
( [ eq : dstar ] ) and ( [ eq : phase ] ) with the above phase choice , the non - diagonal matrix elements can be easily seen to be real .
the same argument also shows that the diagonal matrix elements are automatically real , regardless of the phase choice .
this completes the proof .
we would like to note that while any wilson type dirac operator has a spectrum symmetric with respect to the real axis ( this follows from @xmath31 hermiticity ) , in general there is no simple relation between the eigenvectors corresponding to complex conjugate eigenvalues .
the only exception is the su(2 ) case , when they are related by eq .
( [ eq : phase ] ) .
i thank pierre van baal for discussions , tom degrand and anna hasenfratz for correspondence .
the computations were performed using a computer code based on the milc @xcite code and the arpack @xcite package .
i thank the milc collaboration and the authors of the arpack package for making their code publicly available .
this work has been supported by fom .
de  forcrand , m. garca prez , and i .- o .  stamatescu , nucl .
phys .  * b499 * ( 1997 ) 409 ; c.  michael , and p.s .
spencer , phys .
* d52 * ( 1995 ) 4691 ; d.a .
smith , and m.j .
teper , phys .
* d58 * ( 1998 ) 014505 ; b.  alles , m.  campostrini , a.  di  giacomo , y.  gunduc , e.  vicari , phys .
* d48 * ( 1993 ) 2284 ; a.  hasenfratz , and c.  nieter , phys .
* b439 * ( 1998 ) 366 .
edwards , u.m .
heller , j.  kiskis , and r.  narayanan , phys .
* 82 * ( 1999 ) 4188 ; r.g .
edwards , u.m .
heller , j.  kiskis , and r.  narayanan , chiral condensate in the deconfined phase of quenched gauge theories , hep - lat/9919941 .


a double exchange model was introduced in @xcite to describe the motion of conduction electrons in the background of magnetic ions . in its simplest form
the hamiltonian for the model can be written as : @xmath0 where @xmath1 is the electron annihilation operator on the site @xmath2 of the lattice , @xmath3 is a spin @xmath4 of an electron with implicit summation over repeated spin indices @xmath5 .
] at site @xmath2 , @xmath6 is a core spin .
the first term of ( [ eq : de ] ) describes the hopping of conduction electrons between nearest neighbor sites @xmath2 and @xmath7 of the lattice .
the hopping amplitude is @xmath8 .
the second term is a ferromagnetic hund s interaction between the core spin and the spin of the conduction electron at the same lattice site . to specify the model completely we have to fix the total number of conduction electrons in the sample .
we denote filling factor @xmath9 , where @xmath10 is a total number of sites of the lattice . in this paper
we are interested @xcite in the limit of a very strong hund s interaction @xmath11 and of very large core spins @xmath12 . in ( [ eq : de ] )
we did not include any ( bare ) core spin interactions , e.g. , due to superexchange .
the only interaction between core spins in the model ( [ eq : de ] ) occurs through the exchange by conduction electrons .    because of an infinite hund s coupling @xmath13 an electron spin is always aligned with a core spin on any given lattice site .
therefore , the effective hopping amplitude of an electron between two nearest neighbors is the largest when corresponding core spins are parallel . in this case
an electron can hop between lattice sites without changing its spin orientation .
one concludes that conduction electrons induce a ferromagnetic interaction between core spins @xcite .    in the limit of large core spins
@xmath12 we consider one can use an adiabatic approximation and think of core spins as of slow variables so that the conduction electrons ( fast variables ) are always in the ground state in the background of slowly moving core spins .
it is tempting to average over the motion of conduction electrons ( to integrate conduction electrons out ) and obtain an effective spin model describing the dynamics of core spins ( and electron spins ) induced by charge degrees of freedom . in an almost ferromagnetic state
the average spin per lattice site is @xmath14 corresponding to the sum of a core spin and an average electron spin .
therefore , heuristically , one arrives to the model of ferromagnetically interacting _ fractional _ spins .
however , this immediately raises the question of what one means by a fractional spin .
although the average spin per lattice site is fractional , the total spin operator on a given site may take only values @xmath15 ( no electron on the site ) or @xmath16 ( there is an electron on the site ) .
moreover , a fractional value of @xmath17 is incompatible with quantum mechanics of spins ( see section [ sec : singlespin ] ) . resolving
this apparent controversy is a main goal of this paper .
the paper is organized in the following way .
the aim of the section [ sec : singlespin ] is mainly pedagogical .
we establish notations and remind the reader how to write the path integral for a single spin .
we show how the ambiguity of the berry phase for a single spin requires @xmath17 be an integer . for completeness
we derive the berry phase for a single spin in appendix [ app : pathint ] using a fermionic determinant representation for the action .
the reader familiar with the subject can start from the section [ sec : deinfj ] where we consider the double exchange model ( [ eq : de ] ) in the limit @xmath11 and @xmath12 . in section [ sec : de2 ] we formulate a double exchange model on two core spins and obtain our main result  an effective action for core spins induced by a single conduction electron .
this action ( namely , its kinetic part ) answers the question of how to write berry phase for a system of ( two ) fractional spins .
appendix [ app : diagrams ] contains more accurate calculation of the effective action which uses the method of gradient expansion .
we discuss the topology behind berry phase for fractional spins in section [ sec : tabp ] and conclude in section [ sec : concl ] .
in this section we remind the reader how to write down the path integral for a single quantum spin and discuss topological reasons for spin quantization .
quantum spin @xmath18 is defined as a representation of dimension @xmath19 of an @xmath20 lie algebra with commutation relations @xmath21 & = & 2s^{z } ,       \nonumber \\
\left[s^{z},s^{\pm}\right ] & = & \pm s^{\pm } ,       \label{eq : su2}\end{aligned}\ ] ] where as usual @xmath22 .
an @xmath20 lie algebra has irreducible finite dimensional representations of dimension @xmath19 labelled by spin @xmath15 taking integer or half - integer values so that @xmath23 .
we would like to construct the classical action for the time dependent unit vector @xmath24 which upon quantization by path integral reproduces ( [ eq : su2 ] ) .
one could start directly from ( [ eq : su2 ] ) and construct the path integral using coherent states corresponding to an @xmath20 group ( see e.g. , ref.@xcite ) . in appendix [ app : pathint ]
we give an alternative derivation of this path integral using the fermionic determinant representation .
the result is given by partition function @xmath25 } ,       \label{eq : partition}\end{aligned}\ ] ] where @xmath26 & = &         i ( 2s ) 2\pi \int_{\cal b } d^{2}x\ ,         \frac{1}{8\pi}\epsilon^{abc}\epsilon_{\mu\nu }         n^{a}\partial_{\mu}n^{b}\partial_{\nu}n^{c}.      \label{eq : wznw}\end{aligned}\ ] ] let us explain what ( [ eq : partition],[eq : wznw ] ) mean .
first of all we used an euclidean time representation replacing real physical time @xmath8 by an imaginary time @xmath27 . the change to an imaginary time changes quantum amplitude @xmath28 used in path integral into a `` boltzmann '' factor @xmath29 .
we impose periodic boundary conditions ( in imaginary time ) on @xmath30 so that @xmath31 and @xmath32 is a maximal span of an imaginary time . with these boundary conditions
one can think of time as of a one - dimensional circle with identification of point @xmath33 with @xmath34 .
physical results will be obtained in the limit @xmath35 after an analytical continuation back to the real time @xmath36 .
imaginary time representation is slightly more convenient for calculations and for a discussion of topological aspects of spin quantization .
secondly , the two - dimensional integration in ( [ eq : wznw ] ) is performed over some closed domain @xmath37 of two dimensional plane such that the boundary of this domain is identified with an imaginary time .
the values of @xmath38 field in the area @xmath37 are defined as an extension from the physical trajectory @xmath30 so that @xmath39 is a smooth function of @xmath40 $ ] and auxiliary variable @xmath41 $ ] such that @xmath42 and @xmath43 .
we shall refer to the action ( [ eq : wznw ] ) as to a berry phase @xcite for a spin ( see appendix [ app : pathint ] ) .
although the action ( [ eq : wznw ] ) is given in terms of two - dimensional integral over auxiliary space , the variation of the integrand is a full derivative .
therefore , the variation of berry phase ( [ eq : wznw ] ) is given by a one - dimensional integral over the `` physical '' ( imaginary ) time as @xmath44 we see that the variation of ( [ eq : wznw ] ) does not depend on the choice of an extension of @xmath30 in the unphysical domain @xmath37 . using ( [ eq : variation ] ) one can show that the equation @xmath45 , where @xmath46 is a coupling to an external magnetic field gives a classical equation of motion @xmath47 $ ] .
the latter is the equation of precession of spin in magnetic field @xmath48 .
partition function ( [ eq : partition ] ) with the action ( [ eq : wznw ] ) solve the problem of path integral representation for the algebra ( [ eq : su2 ] ) .
namely , the commutation relations ( [ eq : su2 ] ) can be obtained as a result of quantization @xcite of `` classical action '' ( [ eq : wznw ] ) .      we have already checked that an infinitesimal variation of berry phase @xmath49 $ ] does not depend on the extension of @xmath30 into the auxiliary domain @xmath37 .
how about @xmath50 $ ] itself ?
one can easily show ( see appendix [ app : pathint ] ) that the action ( [ eq : wznw ] ) is equal to @xmath51 where @xmath52 is a solid angle swept by a unit vector @xmath38 during its motion ( see fig.[fig : wznw ] ) .
this angle is not uniquely defined : one can always add or subtract any angle which is a multiple of a full solid angle @xmath53 .
the unit vector @xmath30 draws a closed line on the surface of a sphere with unit radius during its motion in imaginary time .
berry phase is proportional to the solid angle ( shaded region ) swept by the vector @xmath30 .
one can calculate this solid angle by extending @xmath38 into the two - dimensional domain @xmath37 as @xmath54 and calculating ( [ eq : wznw ] ) . ]    therefore , @xmath50 $ ] is not unambiguously defined and is a multi - valued functional of @xmath30 .
multi - valued functionals of this type were studied in both mathematics and quantum field theory ( see e.g. , @xcite ) and are often referred to as wznw term  wess - zumino - novikov - witten term . although @xmath55 is a multiple - valued functional , in partition function ( [ eq : partition ] ) one needs @xmath56}$ ] which is a single - valued functional of @xmath30 under the condition that the coefficient @xmath17 is an integer .
then @xmath57 and the boltzmann weights @xmath58 is defined unambiguously by physical configurations @xmath30 .
one can think of this constraint on the value of @xmath17 as of the _ reason _ for spin to be an integer or a half - integer .
having non - integer @xmath17 in ( [ eq : wznw ] ) amounts to having _ inconsistent _ theory .
there is a one more representation of berry phase which is particularly convenient for the purposes of this paper .
it is a so - called @xmath59 representation .
we write @xmath60 , where @xmath61 is a complex two - component vector satisfying the constraint @xmath62 .
the latter reproduces @xmath63 .    in terms of @xmath64-variable we have for berry phase @xmath26 & = & -i ( 2s ) \int_{0}^{\beta }
d\tau \ ,         z^{\dagger}i\partial_{\tau}z      \label{eq : wznwcp1}\end{aligned}\ ] ] which is apparently independent from the auxiliary coordinate @xmath65 .
this `` simplification '' is , however , deceptive because the equation @xmath60 does not define @xmath64 unambiguously .
one can always change @xmath66 with any function @xmath67 satisfying @xmath68 .
the latter is necessary to conserve periodic boundary conditions @xmath69 .
this `` gauge '' transformation obviously does not change the vector @xmath38 .
the change of wznw action ( [ eq : wznwcp1 ] ) under this transformation is @xmath70= i ( 2s)2\pi k$ ] where @xmath71 is any integer .
we arrive to the same result : the action @xmath72 itself is multivalued while @xmath29 is single - valued under _ the condition _ that @xmath17 is an _ integer number _ !
let us now briefly discuss topological aspects of berry phase ( wznw term ) for a single spin ( [ eq : wznw ] ) .
the classical unit vector @xmath38 from ( [ eq : wznw ] ) takes values on a unit two - dimensional sphere @xmath73  the target space of the problem .
there exists nontrivial expression ( the integrand of ( [ eq : wznw ] ) ) depending on vector @xmath38 and its derivatives which is ( i ) local ( ii ) being integrated over two - dimensional sphere gives integer values for all possible smooth configurations of @xmath38 defined on @xmath73 .
mathematically we say that there exists a closed but non - exact 2-form on @xmath73 ( target space ) or even more formally that the second cohomology group of @xmath73 is @xmath74 .
this allows one to write down the topological term ( [ eq : wznw ] ) where an integration of the mentioned 2-form is performed over a domain @xmath37 with physical time as the boundary of this disk .
the special properties mentioned above guarantee that ( [ eq : wznw ] ) unambiguously defines the path integral ( [ eq : partition ] ) .
let us now consider the double exchange model ( [ eq : de ] ) and take a limit of infinite hund s constant @xmath75
. we will do it in the lagrangian formulation .
first we use the wznw action ( [ eq : wznwcp1 ] ) to write the lagrangian ( in this section we use real time representation ) of the double exchange model corresponding to the hamiltonian ( [ eq : de ] ) @xmath76 - h.         \label{eq : delagr}\ ] ] here @xmath77 is a chemical potential which is chosen to guarantee the correct filling factor of electrons and @xmath78 is a complex vector corresponding to the core spin @xmath79 at the site @xmath2 of the lattice .
we will make the following change of variables in the lagrangian @xmath80 . here
@xmath81 are spinless fermions , @xmath82 is a core spin variable , and @xmath83 is a complex vector orthogonal to @xmath82 , namely , @xmath84 .
physically , @xmath85 is a component of an electron on the site @xmath2 with a spin parallel to the core spin at the same site and @xmath86 is the one with a spin antiparallel to the core spin . taking limit
@xmath87 corresponds to neglecting @xmath88 component of an electron .
we substitute @xmath89 into ( [ eq : delagr ] ) and omitting constant terms obtain @xmath90    \nonumber \\         & + & \sum_{i}(2s+\psi_{i}^{\dagger}\psi_{i } )         z_{i}^{\dagger}i\partial_{t}z_{i}.    \label{eq : delagrproj}\end{aligned}\ ] ]    the first term of ( [ eq : delagrproj ] ) is a kinetic term for projected electrons . let us take a close look at the second and the third terms .
we assume that the motion of spins is much slower then the motion of electrons .
this is justified in the limit of large core spins @xmath12
. then one can think of electrons as moving in the background of almost static @xmath64 s
. the second term in ( [ eq : delagrproj ] ) describes the hopping of these electrons with the effective hopping amplitude @xmath91 .
this amplitude is complex number and has both the magnitude and the phase .
the magnitude of the effective hopping amplitude is given by @xmath92 or @xmath93 where @xmath94 is the angle between two core spins .
in this derivation we used a relation @xmath95 and the known identity on pauli matrices @xmath96 .
the absolute value of the effective hopping is proportional to the cosine of the half of the angle between core spins  @xcite .
it is maximal and equal @xmath8 for parallel core spins and vanishes for antiparallel spins .
the phase of the effective amplitude @xmath97 is not defined unambiguously because of the gauge freedom ( see the previous section ) .
consider , however , an electron moving along some closed trajectory in a background of a smooth configuration of core spins
. then along the path we can write @xmath98 , where @xmath99 is small , and @xmath100 .
the product of all the effective amplitudes along the path is @xmath101 , where we use parameterization @xmath102 of @xmath82 in terms of euler angles @xmath103 , and @xmath85 . along the closed path @xmath104
, so only the second term contributes to the phase picked up by an electron .
it is easy to see that this second term is equal to the half of the solid angle enclosed by vectors @xmath38 along the trajectory of an electron .
we see that ( projected ) electrons moving along some closed trajectory acquire the phase equal to the half of solid angle formed by spins along this trajectory .
one can say , therefore , that they `` feel '' the solid angle formed by core spins as a flux of magnetic field .
this flux is ( modulo @xmath105 ) a gauge invariant and well defined quantity .
finally , the last term in ( [ eq : delagrproj ] ) is a `` modified '' kinetic term for core spins .
one might expect that after the averaging over electrons this term will describe the dynamics of effective spins of magnitude @xmath106 .
however , having non - integer coefficient in front of berry phase results in an inconsistent theory .    to make this seeming paradox more precise we average over electrons accurately by integrating fermions @xcite out of ( [ eq : delagrproj ] ) .
we obtain an effective action for core spins in terms of @xmath64-variables @xmath107    \nonumber \\         & + & \int dt\ ,         \sum_{i } 2s         z_{i}^{\dagger}i\partial_{t}z_{i } .
\label{eq : deeff}\end{aligned}\ ] ] here the expression in the determinant is both an operator in time and @xmath108 matrix in lattice sites .
the matrix element @xmath109 is present only if @xmath110 are nearest neighbors . in the adiabatic limit
there are two main effects produced by an integration over fermionic degrees of freedom .
first , the energy of the filled ( to the filling factor @xmath111 ) fermi sea is added to the action .
this leads to an effective ferromagnetic interaction between core spins .
second , the correction to the berry phase of core spins must reflect the fact that spins of conduction electrons are effectively added to core spins .
the latter correction is the matter of consideration of this paper .
the variation of ( [ eq : deeff ] ) over @xmath112 immediately gives @xmath113 , where @xmath114 are some fractional numbers which in general depend on the configuration of core spins ( @xmath82 s ) .
this implies the presence of the term @xmath115 in the action ( [ eq : deeff ] ) and confirms the result we expected from the heuristic `` averaging '' over electrons .
however , the presence of this term alone renders theory inconsistent .    in the next section
we calculate the fermionic determinant of the type ( [ eq : deeff ] ) albeit for a simplified double exchange model on two lattice sites .
we show that there is an additional term present in the effective action for core spins .
this additional term will make the effective theory consistent as it is shown in section [ sec : tabp ] .
the paradox of having fractional spin in the effective action for a double exchange model is present already for a much simpler double exchange model on two lattice sites with a single conduction electron .
in addition the electron spectrum in such model is discrete which makes integration over fermions in adiabatic approximation to be a well - defined procedure .
therefore , in this section we consider the double exchange model on two lattice sites with one conduction electron and show how to resolve the issue of fractional spin on this simple example .    the projected ( @xmath87 ) double exchange model ( [ eq : delagrproj ] ) on two lattice sites ( in imaginary time ) @xmath116    \nonumber \\         & - & t\left [         \psi_{1}^{\dagger}(z_{1}^{\dagger}z_{2})\psi_{2 }         + h.c.\right ] ,    \label{eq : delagrproj2}\end{aligned}\ ] ] where @xmath117 is a euclidean action .
chemical potential @xmath77 must be chosen to guarantee that there is exactly one conduction electron in the system , i.e. , @xmath118 .
one can repeat the heuristic argument of the previous section replacing @xmath119 in the second term of ( [ eq : delagrproj2 ] ) by its average value @xmath4 .
one obtains two effective spins of magnitude @xmath120 instead of bare core spins .
again , we have the problem of having the fractional coefficient in the berry phase term .    we rewrite the lagrangian ( [ eq : delagrproj2 ] ) as @xmath121         + \psi^{\dagger}(\partial_{\tau}+i\mu+\hat{a})\psi .
\label{eq : proj2}\ ] ] here we defined @xmath122 and @xmath123 matrix @xmath124 @xmath125 where we introduced explicitly the phase and the magnitude of an effective hopping amplitude @xmath126 the absolute value of an effective hopping amplitude @xmath127 , where @xmath128 is an angle between core spins @xmath129 and @xmath130 .
it is equal to zero only for strictly antiparallel core spins @xmath131 . for fixed directions of core
spins the matrix @xmath124 becomes constant matrix with eigenvalues @xmath132 corresponding to the two - level spectrum of an electron hopping between two core sites in double exchange model with infinite @xmath13 .
two - level spectrum is symmetric and one can choose the chemical potential @xmath133 so that there is only one filled energy level at any given moment of time .
then the spinless fermion @xmath134 ( electron with projected spin ) occupies the lowest energy level @xmath135 which is minimal ( @xmath136 ) when core spins are parallel .
changing the angle between core spins requires the energy of the order of @xmath8 while we are interested in low energy dynamics with typical values of energy @xmath137 .
therefore , we assume that @xmath138 does not vanish and is of the order of @xmath8 .
the grassman variable @xmath139 corresponding to the conduction electron must satisfy antiperiodic boundary conditions in imaginary time @xmath140 : @xmath141 .
we represent @xmath142 as a product of spinless fermion @xmath143 and spin @xmath4 boson @xmath144 , which after projection @xmath87 becomes @xmath64-boson corresponding to core spins .
the boundary conditions for @xmath145 and @xmath146 are antiperiodic and periodic respectively .
then we have for @xmath147 and @xmath148 ( see ( [ eq : proj2],[eq : notation ] ) ) @xmath149 where @xmath71 is an integer  the winding number of the phase @xmath150 in imaginary time @xmath151 .
the effective action of core spins is given by @xmath152 \,d\tau         + w_{e}^{ind } ,    \label{eq : ea}\ ] ] where the effective action of core spins induced by fermions @xmath153 here the operator @xmath154 is linear differential operator with explicit dependence on time @xmath140 through the @xmath123 matrix @xmath155 and we put @xmath133 corresponding to having only one conduction electron in the system . to calculate the functional trace ( [ eq : ft2 ] ) we find the eigenvalues of this operator by solving differential equation @xmath156 with boundary conditions ( [ eq : boundarypsi],[eq : boundarya ] ) and
calculate the sum of logarithms of all eigenvalues @xmath157 .
matrix ( [ eq : a2 ] ) can be rewritten as @xmath158 where @xmath159 we would like to solve ( [ eq : equation1 ] ) with the matrix @xmath155 from ( [ eq : a20 ] ) in an adiabatic approximation , e.g. , assuming that all time derivatives are much smaller than the energy scale @xmath160 .
first we `` unwind '' wavefunction @xmath161 by a unitary transformation @xmath162 where @xmath163 means the derivative of @xmath150 over @xmath164 .
then , instead of ( [ eq : equation1 ] ) the function @xmath165 satisfies @xmath166 where we defined `` zero order '' operator @xmath167 and perturbation @xmath168 which is proportional to time derivatives of core spins as @xmath169 now we can solve ( [ eq : ptequation ] ) using the perturbation theory in @xmath168 .
the boundary conditions of @xmath165 are given by ( [ eq : boundarypsi ] ) and ( [ eq : transformation ] ) as @xmath170 where berry phase @xmath171 where we used that @xmath172 ( @xmath71 is an integer , see ( [ eq : boundarya ] ) ) and that @xmath173 .
we are going to find the eigenvalues of ( [ eq : ptequation ] ) using the perturbation theory in @xmath168 .
this is justified in an adiabatic limit when @xmath174 proportional to the rate of change of core spins is much smaller than the typical fermionic energy @xmath138 or in other words when the conduction electron is the fast degree of freedom compared to core spins .
the solutions of the zero order equation @xmath175 are labelled by an integer number @xmath176 and by plus or minus sign @xmath177 with eigenvalues @xmath178 where @xmath179 there are no corrections to these eigenvalues in the first order of perturbation theory in @xmath168 . therefore ,
up to the second order of perturbation theory we obtain @xmath180 where we have restored the chemical potential @xmath77 which should be @xmath133 for the problem with a single conduction electron .    the effective action ( [ eq : ft2 ] ) is given by @xmath181 .
\label{eq : ft3}\ ] ] to calculate the obviously divergent sum in ( [ eq : ft3 ] ) we regularize it by subtracting the effective action for the system without conduction electrons .
that is we calculate the difference between ( [ eq : ft3 ] ) calculated at @xmath133 and the one calculated at @xmath182 . @xmath183 .
\label{eq : ft4}\end{aligned}\ ] ] this sum is still logarithmically divergent but one can show that the divergent part is a constant not depending on physical parameters @xmath184 and @xmath185 . calculating the sum in ( [ eq : ft4 ] )
we obtain    @xmath186    in the limit @xmath182 we obtain @xmath187 this is the final result for the effective action calculated at finite `` temperature '' @xmath32 .
notice that the only imaginary part of a euclidean effective action is the berry phase @xmath185 .
we have obtained this phase in the first order of perturbation theory . however , the same expression is also valid in all orders of perturbation theory .
it is easy to show for our simple model because the higher orders of perturbation theory can only give corrections to the real part of ( [ eq : firstorder ] ) which does not contribute to the imaginary part of an euclidean action .
we consider the higher orders of perturbation theory using the regular gradient expansion in appendix [ app : diagrams ] .
in fact , the imaginary part of a euclidean action has a topological nature and , therefore , can not be changed by small perturbations ( see the discussion in section [ sec : tabp ] ) .    in the `` zero temperature '' limit @xmath188
we have @xmath189 or explicitly @xmath190 the latter result can of course be obtained just by replacing the discrete sum over @xmath176 in ( [ eq : ft4 ] ) by an integral . for the full effective action ( [ eq : ea ] ) we obtain @xmath191 ,    \label{eq : kin } \\
w_{e}^{dyn } & = & -\int_{0}^{\beta}\delta\,d\tau       + \ldots .
\label{eq : dyn}\end{aligned}\ ] ]    the first term of the dynamic part of an effective action @xmath192 is just the time integral of the energy of the ground state of the conduction electron in a fixed instantaneous configuration of core spins .
the dots denote the higher order corrections to ( [ eq : dyn ] ) .
these corrections are small if the background core spins are moving adiabatically and they are diverging if derivatives @xmath193 are bigger than @xmath138 . for such processes an adiabatic approximation fails .
in particular , if core spins have opposite directions @xmath194 we have @xmath195 and an adiabatic approximation is not applicable at any rate of the change of core spins .
the berry phase term ( [ eq : kin ] ) is the main result of this paper .
we see that the the correction to the berry phase of bare core spins induced by the conduction electron is twofold .
first , as we expected averaging over fermions in ( [ eq : delagrproj2 ] ) there is a change of the coefficient in berry phases of core spins from @xmath17 to @xmath196 .
in addition , there is a correction @xmath197 , where the integer number @xmath71 is a winding number ( [ eq : boundarya ] ) of the phase @xmath150 in imaginary time .
the latter term is a full time derivative and , therefore , does not change the classical dynamics of spins .
classically , the effective `` core '' spins behave as if they had the value of @xmath120 !
let us show now that the full berry phase ( [ eq : kin ] ) is consistent with quantum mechanics of spins and , therefore , answers the question of how to write down the berry phase for two fractional spins @xmath120 .
as we have seen at the end of the section [ sec : singlespin ] the inconsistency occurs if the partition function is not invariant under gauge transformations of @xmath64 variables .
we perform the gauge transformation @xmath198 independently for first and second core spins . in order not to change the periodic boundary conditions of @xmath199 we require @xmath200 , where @xmath201 are arbitrary integer numbers . under this transformation
the three terms of ( [ eq : kin ] ) depending on @xmath202 , @xmath203 , and @xmath150 change by @xmath204 , @xmath205 , and @xmath206 respectively .
it is easy to see that neither of these three terms gives an invariant contribution to a partition function .
e.g. , for @xmath207 the change of the first term is @xmath208 which gives @xmath209 being exponentiated .
however , the total kinetic term ( [ eq : kin ] ) is equal to the sum of all these terms and changes under the gauge transformation by @xmath210 $ ] which is a multiple of @xmath105 and , therefore , does not change the weight @xmath211 .
this means that neither of three terms considered makes any sense separately from the others .
only their gauge invariant ( modulo @xmath105 ) sum has a physical meaning of the combined berry phase of the system of two core spins plus one conduction electron .
we have seen in section [ sec : singlespin ] that the existence of wznw term for a single unit vector ( single spin ) is due to the fact that the target space @xmath73 has a nontrivial second cohomology group @xmath74 .
the addition of wznw term to the action with the coefficient @xmath17 is consistent only if this coefficient is quantized ( is an integer number ) . for two independent spins the target space becomes @xmath212 and
one obviously can write two independent wznw terms corresponding to @xmath213 .
this is apparently true for two independent spins .
however , the result is different for the conduction electron induced dynamics of core spins in the double exchange model .
let us discuss why this is so . to obtain the induced dynamics of core spins
we have excluded the charge degree of freedom of conduction electron integrating it out in an adiabatic approximation .
this approximation , however , breaks when the the rate of the change of core spins becomes comparable to the energy scale @xmath214 given by the levels of conduction electron in the background of time independent core spins .
there is a special point though when core spins are opposite to each other . at this point @xmath215 and an adiabatic
approximation breaks no matter how slow the motion of core spins is .     two degenerate states with antiparallel core spins .
an electron with spin aligned along the core spin can not hop from one lattice site to another and can be either on the right or on the left lattice site . ]
this is quite obvious from the physical point of view .
consider two time independent core spins opposite to each other .
the effective hopping amplitude of a conduction electron is @xmath216 and one has two different configurations .
the electron can be either on the left lattice site or on the right one ( see figure [ fig : oppspins ] ) and the matrix element of the hopping term is zero between those states .
this means that the ground state of the system is doubly degenerate .
the charge degree of freedom ( left or right position of the electron ) is essential for this particular configuration of core spins and can not be integrated out . in our derivation
we used an adiabatic approximation and neglected such configurations .
this is justified in the limit @xmath12 for ferromagnetically interacting core spins .
however , this approximation drastically changes the topology of the system . essentially , we prohibited the configurations @xmath194 . the relevant target space changes from @xmath212 for two independent spins to @xmath217 .
the latter expression denotes the direct product @xmath212 with points @xmath194 ( topologically equivalent to one more @xmath73 ) removed .
the resulting target space @xmath218 is topologically equivalent to @xmath73 . to prove it one can show
is , in fact , homotopic to a tangent bundle to @xmath73 which is in turn can be contracted onto @xmath73 . ] that it can be contracted along itself onto @xmath73 .
we illustrate the similar phenomenon on the corresponding one - dimensional example in figure [ fig : cutting ] .
the contraction @xmath219 is shown as a one - dimensional illustration of @xmath220 .
a ) the torus space @xmath221 is represented as a square with the opposite sides identified according to the direction of arrows .
the diagonal shows the points given by @xmath194 .
b ) the space @xmath219 is obtained by the removal of the diagonal .
c ) one pair of opposite sides is glued .
d ) the second pair is glued making a cylinder .
e ) the cylinder is contracted onto the unit circle . ]    therefore , the relevant target space @xmath222 . to check for allowed wznw terms
we calculate @xmath223 .
there is only one wznw term allowed for the chosen target space !
the existence of only one allowed topological term explains why only the full berry phase ( [ eq : kin ] ) is well defined . the most general form ( in @xmath59 representation ) of the topological term for two spins
is given by @xmath224 .
\label{eq : tterm}\end{aligned}\ ] ] coefficients @xmath225 satisfy quantization conditions ( or gauge invariance conditions ) .
namely , @xmath226 and @xmath227 must be integers .
we see that for the double exchange model of two core spins and a single conduction electron we have ( [ eq : kin ] ) which corresponds to particular integer values @xmath228 and @xmath229 .
if by some reason the conduction electron spends unequal time at two sites ( e.g. , the local magnetic field acting on electron spin on one of sites is added ) we have an asymmetric situation @xmath230
. however , the combination @xmath226 and @xmath227 must stay integer to guarantee the consistency of the theory .
in this paper we considered the infinite hund s coupling limit of the double exchange model with a given concentration @xmath111 of conduction electrons . in the limit of large core spins
@xmath12 one can think of electrons and core spins as of fast and slow degrees of freedom respectively .
the superficial averaging over the fast motion of electrons leads to an inconsistent theory with effective fractional spins . using the example of the double exchange model on two lattice sites we showed that in addition to the change of the value of the effective core spin on site from @xmath15 to @xmath231 a new topological term arises in the effective action . with this term
the exponent of the effective action becomes single - valued and well defined .
we showed that the full berry phase ( [ eq : kin ] ) reflects the change of topology of the system which occurs due to the adiabaticity condition  the condition justifying the averaging over the charge degrees of freedom .
the generalization of our results to the double exchange model for many spins is relatively straightforward .
the only major difference with the case of two spins is that the average occupation number of a given lattice site by conduction electrons depends on the configuration of core spins .
this makes the magnitude of a `` fractional '' spin on the given site a fluctuating quantity with correspondent modifications of berry phase terms . the latter effect is not expected to be strong in the ferromagnetically ordered state .
we are grateful to boris spivak who stimulated us to think about fractional spins in ferromagnetic metals .
the discussion of topological aspects of berry phase with dmitri ivanov is deeply appreciated .
a.  g.  a. is alfred p.  sloan research fellow .
the work of ar .
a. is supported by dr project 200153 at los alamos national laboratory .
in this section for the sake of completeness we rederive the well - known path integral representation for the quantum mechanics of a single spin .
we reduce the problem to a calculation of a simple fermionic determinant .
the method we use differs from the commonly used coherent state representation ( see e.g. , @xcite ) for the path integral .
it is conceptually simpler but involves the gradient expansion calculation of a fermionic determinant .
it turns out that the similar determinant is needed for our analysis of a double exchange model for two core spins .
it is easy to check that the representation @xmath232 with constraint @xmath233 is a faithful representation of ( [ eq : su2 ] ) with @xmath234 . here @xmath235 is a set of pauli matrices and @xmath236 with @xmath5 are annihilation ( creation ) fermi operators respectively which satisfy the following ( anti)commutation relations @xmath237 and @xmath238 .
in fact , this representation reflects how spin operators appear in nature with @xmath134 being an annihilation operator of the physical electron .
we consider the following hamiltonian @xmath239 , where @xmath240 is a spin of a fermion and @xmath138 is a strength of the coupling of this spin with the unit vector @xmath38 .
the correspondent action is @xmath241\psi ,   \label{eq : actionspin}\ ] ] where we introduced a notation @xmath242 .
one can think of @xmath243 as of the magnetic field acting on the fermionic spin @xmath18 . in the limit of an infinite coupling @xmath244 ( infinite magnetic field )
the fermionic spin @xmath18 will be `` frozen '' along the classical vector @xmath38 so that @xmath38 in the limit of infinite @xmath138 becomes a direction of quantization axis of a quantum spin of a fermion .
if now we change the direction of the classical vector @xmath38 with time , the fermionic spin will follow the direction of @xmath38 .
therefore , we expect that in the limit of an infinite @xmath138 an effective dynamics of classical vector @xmath38 induced by a fermion will be that of the quantum spin @xmath4 .
the path integral representation for spin @xmath4 is then given by @xmath245 } ,       \label{eq : pathint}\end{aligned}\ ] ] where @xmath246=\lim_{\delta\to\infty}w^{eff}[\vec{n}]$ ] is a classical action corresponding to a quantum spin @xmath4 .
this classical action is a functional of the trajectory of the time dependent unit vector @xmath24 .
the effective action @xmath247 $ ] is defined as @xmath248 } & = & \int d\psi d\psi^{\dagger } \ ;            e^{i\int dt\ , \psi^{\dagger }            \left[i\partial_{t}+\delta\hat{n}\right]\psi}.       \label{eq : fermint}\end{aligned}\ ] ] the functional integration in ( [ eq : fermint ] ) is performed over grassman variables@xcite @xmath249  classical analogs of fermionic annihilation ( creation ) operators .
integrating over fermions we obtain @xmath250 & = & -i \ln \,{\rm det}\ ,            \left[i\partial_{t}+\delta\hat{n}\right ] .
\label{eq : effdet}\end{aligned}\ ] ]    to generalize the result to an arbitrary spin @xmath15 we can consider the same fermionic theory ( [ eq : fermint ] ) but with @xmath17 independent species of fermions coupled to the same vector @xmath38 . in this case
all spins @xmath4 of fermions will be aligned along @xmath38 so that ( [ eq : pathint ] ) will correspond to quantum spin @xmath15 .
integration over all fermionic species is independent and we have for an effective action @xmath250 & = & -i 2s \ln\,{\rm det}\ ,            \left[i\partial_{t}+\delta\hat{n}\right ] .
\label{eq : effaction}\end{aligned}\ ] ] the only problem which is left is technical .
it is to calculate the fermionic determinant ( [ eq : effaction ] ) in the limit of an infinite @xmath138 .
one can obtain the result without ( almost ) any calculations .
let us introduce the complex , two - component , normalized vector @xmath251 , @xmath252 so that it is an eigenvector of @xmath253 with the eigenvalue @xmath254 : @xmath255 .
for a time independent @xmath253 this vector is a ground state of our hamiltonian @xmath256 with the ground state energy @xmath257 .
the other eigenstate has an energy @xmath258 and in the limit of an infinite @xmath138 is never excited by a motion of the unit vector @xmath38 .
therefore , in the limit of an infinite @xmath138 we can use @xmath259 with @xmath165  spinless ( projected ) fermion . substituting this into the action ( [ eq : actionspin ] )
we obtain @xmath260.\ ] ] the variable @xmath165 does not have any dynamics in the limit of an infinite @xmath138 .
we average over @xmath165 using @xmath261 and obtain @xmath262 for spin @xmath4 . in case of spin @xmath15
one has @xmath17 species of fermions , @xmath263 and @xmath264 .
this is the desired result for the action for a single spin @xmath15 .
however , the main point of this paper is to show that the `` naive '' averaging over fermions instead of the careful determinant calculation can result in error . in this case
one can show that due to topological reasons the calculation we performed is a legitimate one . to add even more weight to this statement
we proceed with a formal calculation of the fermionic determinant ( [ eq : effaction ] ) and show that it gives the same result for a single spin action .      before going to the calculation of fermionic determinant ( [ eq : effaction ] )
we will perform wick rotation of time and go to an imaginary time representation .
imaginary time representation is slightly more convenient for calculations for reasons that will be clear later .
the change to an imaginary time changes quantum amplitude @xmath28 into a `` boltzmann '' factor @xmath29 , where subscript `` e '' stands for `` euclidean '' . in imaginary time @xmath265
we have instead of ( [ eq : effaction ] ) @xmath266 & = & - 2s \ln\,{\rm det}\ ,            \left[\partial_{\tau}-\delta\hat{n}\right ]       \label{eq : effactione}\end{aligned}\ ] ] with the partition function for a single spin @xmath25 } ,       \label{eq : pathinteapp}\end{aligned}\ ] ] where @xmath50=\lim_{\delta\to\infty}w_{e}^{eff}[\vec{n}]$ ] .
we impose periodic boundary conditions ( in imaginary time ) on @xmath30 so that @xmath31 and @xmath32 is a maximal span of an imaginary time . with these boundary conditions
we can think of time as of the one - dimensional circle with identification of points @xmath33 with @xmath34 .
physical results will be obtained in the limit @xmath35 after analytical continuation back to real time @xmath36 .    using the method presented in appendix [ app : diagrams ]
one can calculate the effective action ( [ eq : effactione ] ) in adiabatic approximation with an arbitrary accuracy .
however , we are interested in the limit @xmath244 . in this limit all terms of gradient expansion in @xmath267 except for the first ( topological ) one will vanish . here
we show how to obtain the most important ( topological ) @xmath184-independent part of the effective action ( [ eq : effactione ] ) which is the exact result for the action for a single spin .
determinant of an operator does not change if one applies a unitary transformation to it .
this is , however , a potentially dangerous procedure in the context of gradient expansion as one has to make sure the boundary conditions for electrons do not change in the process of such transformation .
it turns out that this is not an issue for the single spin problem .
however , e.g. , for the double exchange problem the boundary condition of fermions do change in a non - trivial way under the unitary rotation ( see section [ sec : de2 ] ) .
let us use the transformation @xmath268 such that @xmath269 .
this transformation @xmath268 is easily found using the so called @xmath59 representation . in this representation
we introduce a two - component complex vector @xmath64 with the constraint @xmath270 such that @xmath271 .
then the matrix @xmath268 can be written as @xmath272 we write the equation ( [ eq : effactione ] ) in the form @xmath273 & = & - 2s \,{\rm{tr}\,}\ln            \left[\partial_{\tau}+\hat{a}-\delta\sigma ^{3}\right ] ,       \label{eq : effactione2}\end{aligned}\ ] ] where @xmath274 and @xmath275 means both functional trace and the trace over @xmath276-matrices .
the matrix elements of @xmath155 are proportional to the first time derivative of spin variables and are small in comparison to @xmath138 .
we can expand the logarithm in ( [ eq : effactione2 ] ) in powers of @xmath155 .
we keep only the first @xmath155 dependent term as the other terms will have additional powers of @xmath138 in denominators .
( one can obtain those terms in a regular way explained in appendix [ app : diagrams ] ) @xmath277 & = & 2s \,{\rm{tr}\,}\frac{\partial _ { \tau }       + \delta \sigma ^{3 } }       { -\partial ^{2}_{\tau } + \delta ^{2}}\hat{a }    \nonumber \\       & = & i\int \frac{d\omega } { 2\pi }       \,\mbox{tr}\,\frac{-\omega + i\delta \sigma ^{3 } }       { \omega ^{2}+\delta ^{2}}\hat{a}(\omega = 0 ) \nonumber \\       & = & 2s\frac{1}{2}\int_{0}^{\beta } d\tau \,\mbox{tr}\ ,       \sigma ^{3}\hat{a }    \nonumber \\       & = & - i 2s\int_{0}^{\beta } d\tau\ , z^{\dagger } i\partial _ { \tau } z.       \label{eq : effactione3}\end{aligned}\ ] ]    it is instructive to rewrite ( [ eq : effactione3 ] ) in terms of the original variables @xmath38 .
we parameterize complex vector @xmath278 .
it is easy to see that the angles @xmath128 and @xmath279 are the polar and azimuthal angles corresponding to the direction of the unit vector @xmath38 .
the phase @xmath134 can be chosen arbitrarily as it does not affect @xmath38 .
up to the full time derivative of @xmath134 we obtain @xmath280 , where @xmath52 is a solid angle enclosed by the unit vector @xmath38 during its time evolution .
one can check that this solid angle can be rewritten in the form ( [ eq : wznw ] ) which is explicitly rotationally invariant .
in this appendix we calculate the effective action ( [ eq : ft2 ] ) which can be written as @xmath281 where the differential operator @xmath282 is defined as @xmath283 with matrix @xmath155 given by ( [ eq : a2],[eq : a20 ] ) .
we rewrite @xmath284 where @xmath285 , @xmath286 , and @xmath287 .
we use an adiabatic approximation using the ( gradient ) expansion in @xmath267 .
the operator @xmath282 contains the magnitude of the effective hopping @xmath138 which is slowly varying but it is not small itself . to avoid this difficulty we calculate instead of effective action its variation and
then restore the result from this variation @xmath288 ,    \nonumber\ ] ] where @xmath289 let us now calculate the functional trace in the basis of plane waves @xmath290 here the integral over time @xmath140 is the calculation of the diagonal matrix element @xmath291 and integration over @xmath292 is the calculation of the functional trace . the residual trace @xmath293 is the trace in the space of @xmath123 matrices .
`` pulling '' @xmath294 from right to the left we have @xmath295 in the latter expression the partial derivative @xmath296 is assumed to act `` to the right '' .
this means that when we expand @xmath297    \nonumber \\          & = &        \int dt\left[\delta
l_{0}+\delta l_{1}+\delta l_{2}+\dots        \right ] ,    \label{eq : expansion}\end{aligned}\ ] ] all derivatives are assumed to act on all terms to the right .
e.g. , @xmath298 .
we introduced @xmath299 and @xmath300 we notice that @xmath301 because @xmath302 by adiabatic approximation ( it contains time derivative of core spins ) . therefore , the denominator of ( [ green ] ) is never zero and integrals ( [ deltal ] ) are not singular .
the expansion ( [ eq : expansion ] ) is essentially the gradient expansion , i.e. , the expansion in the number of time derivatives .
one can see that @xmath303 contains at least one time derivative , @xmath304  at least two etc .
the calculation of @xmath305 is straightforward and gives @xmath306     \label{deltal0}\ ] ] or removing variation @xmath307 to obtain the latter expression we expanded up to the second order in time derivatives using @xmath308 .
let us now prove that the whole dependence of the full result ( [ eq : expansion ] ) on @xmath309 comes from @xmath310 .
varying the action ( [ eq : proj2 ] ) over @xmath309 we obtain @xmath311 .
this statement is exact because the total number of electrons in the system is fixed @xmath312 . for the effective action we have , therefore , @xmath313 .
however , from ( [ deltal0 ] ) we have @xmath314 which gives @xmath315 and all higher order terms @xmath303 , @xmath304 etc .
do not depend on @xmath309 .
this simplifies further calculations allowing to use @xmath316 for all higher order terms .    as a result of calculations one obtains from ( [ deltal ] ) @xmath317 this is identical in form to a berry phase of a single spin @xmath4 given by ( [ eq : variation ] ) with identification @xmath318 .
therefore , @xmath303 is given by formula analogous to the berry phase for a single spin ( [ eq : wznw ] ) and is to be interpreted as a half of solid angle swept by vector @xmath319 during its time evolution .
assume for a moment now that @xmath320 .
then the vector @xmath321 moves along the equator of unit sphere and covers the solid angle @xmath322 during its motion . in this approximation
we obtain @xmath323 which gives half of the solid angle after time integration .
is known as the method of embedding in quantum field theory @xcite ] calculating corrections due to nonzero @xmath324 from ( [ eq : l1 ] ) we obtain @xmath325 where we expanded the latter expression up to the terms of the second order in time derivatives .
finally , combining terms ( [ fulll0],[fulll1],[fulll2 ] ) we obtain an effective action up to the second order terms @xmath327 it is interesting to notice that although neither of contributions ( [ fulll0],[fulll1],[fulll2 ] ) is gauge invariant , the full action ( [ fullle ] ) is gauge invariant .
this means that the final result ( [ fullle ] ) can be rewritten in terms of two unit vectors @xmath328 and @xmath329 which represent physical degrees of freedom
core spins of double exchange model ( [ eq : de ] ) .