chiral symmetry plays an essential role in particle physics . in qcd
it is a global symmetry which forbids an additive quark mass renormalization and whose spontaneous breakdown provides us with ( nearly ) goldstone bosons with their specific interactions .    on the other hand , in electroweak interactions
chiral symmetry plays an even more fundamental role : it is a local gauged symmetry .
its presence is necessary to provide a renormalizable theory .
accordingly , when trying to regularize these theories on the lattice ( which is the only non - perturbative regularization known for gauge theories ) one meets two basic problems . as it is ( or rather _ was _ ) well known one can not realize chiral symmetry exactly on the lattice without sacrificing some even more important principles , like locality or absence of extra physical particles . this is the content of the famous no - go theorem by nielsen and ninomiya @xcite . in wilson s formulation of lattice qcd chiral symmetry
is explicitly broken by an irrelevant term which does not affect the continuum limit .
it causes , however , o@xmath0 lattice artifacts and other unwanted effects at finite lattice spacing .    for the electroweak theory , on the other side , the problem of chiral symmetry on the lattice is a principal one : is it possible to define a chiral gauge theory non - perturbatively ?
this second problem is much more difficult , and assumes that the first one  a solution for global chirality  has been successfully solved .
several ways to solve these problems have been proposed .
these were mostly covered on previous lattice conferences . here
i mention only the domain wall fermions @xcite and the related overlap formalism @xcite . on
can state that these approaches provide a solution to the problem of global chirality .
the numerical results using domain wall fermions @xcite are encouraging . in the overlap formalism also
important steps were made towards the solution of the gauged chiral symmetry ( for references see ref .
@xcite ) .
it was unclear , however , to what extent are the theoretical constructions used in these approaches necessary , or what is the basic ingredient behind the chiral properties observed . as it became clear recently ,
the basic relation which has to be satisfied to realize chiral symmetry at finite lattice spacing , is the ginsparg - wilson ( gw ) relation @xcite suggested a long time ago !    in this talk
i shall concentrate on the role of the gw relation and its consequences .
some of the results were obtained earlier in the approaches mentioned above , but even these are more transparent within the new framework .
consider the free fermionic action @xmath1 describing a massless fermion on the lattice .
( colour , dirac and flavour indices are suppressed . )    the desirable properties of the dirac operator @xmath2 should be    * @xmath3 is local ( bounded by @xmath4 ) * @xmath5 for @xmath6 * @xmath7 is invertible for @xmath8 ( no massless doublers ) * @xmath9    here @xmath7 is the fourier transform of @xmath3 .    the theorem by nielsen and ninomiya @xcite states that properties ( a)-(d ) can not hold simultaneously .
following wilson , to avoid the massless doublers present in the naive discretization of the dirac operator , one introduces a chiral symmetry breaking , irrelevant operator .
the corresponding wilson - dirac operator is given by @xmath10 where @xmath11 ( @xmath12 ) is the forward ( backward ) lattice derivative .
due to the explicit breaking of chiral invariance , however , some extra problems appear :    \(1 ) chiral symmetry is broken at finite lattice spacing by o(@xmath13 ) lattice artifacts and it is recovered only in the continuum limit .
\(2 ) one has to fine tune the bare quark mass since there is an additive quark mass renormalization .
\(3 ) operators of different chiral representations get mixed , etc .
one of the conditions of the no - go theorem was locality .
the natural definition ( used here ) is that the dirac operator in an arbitrary gauge field background falls off exponentially i.e. satisfies the bound @xmath14 where @xmath15 and @xmath16 are independent of the gauge field @xmath17 .
in addition , it is also assumed that @xmath18 depends negligibly on @xmath17 s far apart from @xmath19 and @xmath20 , i.e. that @xmath21 also falls exponentially in @xmath22 , @xmath23 and @xmath24 .
it is generally assumed that such interactions satisfy universality .
note that it is physically too restrictive to consider only nearest neighbour interactions or even those with strictly finite support . in a real system
the interaction coefficients die away exponentially rather then having a finite number of strictly nonzero terms .    on the other hand ,
interactions decreasing more slowly then exponentially should be considered as non - local and non - acceptable , since universality is not expected to hold for these
. some non - local interaction can give the same correlations as a local one ( e.g. when a bad rg transformation is used , or when some physical fields are integrated out ) .
this is , however , no sufficient excuse to use them  a slight change in the interaction could modify the picture drastically .
theoretically any decay rate @xmath25 in eq .
( [ bd ] ) is allowed . in practice , however , one works at finite lattice spacing and @xmath26 must hold .
for that reason it is important to have a sufficiently small interaction range .    in principle
it could happen that the decay rate of @xmath18 is not uniformly bounded by an exponential , but there exists a bound @xmath27 ( with @xmath28 ) .
in general , such interactions should be considered non - local .
however , it may happen that those configurations when @xmath29 are so strongly suppressed that they practically ` never ' occur in a mc simulation .
this is a potentially dangerous possibility , but could perhaps be accepted with extra caution .
back in 1982 ginsparg and wilson @xcite suggested a way to avoid the no - go theorem and preserve consequences of the chiral symmetry .
they suggested to require instead of relation @xmath9 the following milder condition @xmath30 here @xmath13 is the lattice spacing and @xmath31 is a _ local _ operator . for simplicity it will be assumed also that @xmath31 is trivial in dirac indices .
that @xmath31 commutes with @xmath32 follows from eq .
( [ 1 ] ) .
the locality of @xmath31 expresses the requirement that the propagating states are effectively chiral , so at distances larger than the range of @xmath31 its presence is not felt . obviously , this is a highly nontrivial condition since the propagator @xmath33 on the lhs .
is non - local .
( the coefficient 2 on the rhs . is for historical reasons only . )    accordingly , the dirac operator @xmath2 should satisfy the ginsparg - wilson relation @xmath34    in ref .
@xcite this relation originates from rg considerations applied to the free fermionic action and has been generalized to fermions in a gauge field .
the local operator @xmath31 comes from a chiral symmetry breaking term in the _ block transformation _ only , not from the original action of the rg procedure .
for that reason the physics described by @xmath2 is expected to be chirally invariant .
( for a chirally symmetric block transformation @xmath2 becomes chirally invariant but non - local , in accordance with the no - go theorem . )
two comments are in order here :    \(1 ) the naive commutation property for @xmath2 is recovered in the continuum limit @xmath35 .
\(2 ) the rhs . of eq .  (
[ 2 ] ) is zero on solutions , i.e. for @xmath36 .
ginsparg and wilson suggested that a solution to eq .
( [ 2 ] ) is the _
mildest _ way to break chiral invariance : it should maintain the physical consequences of chiral invariance ( soft pion theorem , etc . ) .
however , no solution was found for the interacting case ( qcd ) , and the paper has been practically forgotten for 15 years !    last year peter hasenfratz realized @xcite that the fixed point ( fp ) action ( or classically perfect action ) for qcd @xcite satisfies the gw relation .
this observation revived the interest in this relation .
let me recall the definition of the fp action .
for the pure gauge part of the action one has the saddle point ( i.e. classical ) equation @xmath37 here the gauge field @xmath38 lives on the coarse lattice , @xmath17 on the fine one , and @xmath39 is the blocking kernel .
this is a recursive definition , but its recursive evaluation converges rapidly since the minimizing field @xmath40 is much smoother than the original @xmath38 .    in the classical approximation
the fermionic part remains quadratic , and the corresponding propagator satisfies the recursion relation @xmath41 here @xmath42 label the sites on the coarse lattice , @xmath43 on the fine lattice , while on the rhs .
@xmath40 , the minimizing configuration of the pure gauge problem is taken .
the term proportional to @xmath44 comes from the chiral breaking term in the blocking kernel for fermions ; @xmath45 characterizes the averaging procedure for the fermions in rg transformation
it is assumed to be trivial in dirac indices and usually it only extends over the hypercube .
applying eq .
( [ 4 ] ) recursively and using the fact that the original @xmath2 ( in the continuum ) anti - commutes with @xmath32 , it is easy to see that @xmath46 satisfies eq .  ( [ 1 ] ) with @xmath47 ( in fact , for the rg transformation used , @xmath48 or it lives on the hypercube , i.e. is even ultra - local . )    in ref .
@xcite hasenfratz has demonstrated that , indeed , the gw relation implies the chiral ward identities .
the argument is the following : when the symmetry breaking part appears in some expressions ( integrated out in fermionic fields ) @xmath49 then the @xmath33 factors from contractions of the fermionic fields are cancelled by the @xmath2 s on both sides in the breaking term , leaving only the local operator @xmath31 . as a consequence , the symmetry relation @xmath50 is satisfied up to a contact term  i.e. one obtains the ward identities . in this paper
it has also been shown that the proper order parameter for spontaneous chiral symmetry breaking is given by @xmath51 it has been proven as well that the current is not renormalized ( @xmath52 ) and operators in different chiral representations do not mix .
i shall discuss these issues later , in the light of newer developments .    in ref .
@xcite it has been shown that there is an exact atiyah - singer index theorem on the lattice , and the fermionic spectrum ( in arbitrary gauge background ) has nice chiral properties .
the discussion was based on the fp action , but it used mostly only the gw relation to establish these facts .    as mentioned in the introduction ,
the domain wall and the related overlap formalism also produced dirac operators with nice chiral properties .
neuberger @xcite has pointed out that the dirac operator obtained earlier @xcite also satisfies the gw relation . to be on more general ground ,
i discuss first the overlap dirac operator .
this is an explicit construction .
introduce @xcite @xmath53 where @xmath54 is the standard massless wilson - dirac operator given by eq .
( [ dw ] ) .
note that @xmath55 then define @xmath56 it satisfies the gw relation @xmath57 ( i.e. eq .
( [ 2 ] ) with @xmath58 ; sometimes this special case will be referred to as the gw relation ) . introducing @xmath38 through @xmath59 together with the relation @xmath60 one
concludes that eq .
( [ 12 ] ) is equivalent to @xmath61 i.e. @xmath38 is unitary . in other words ,
@xmath62 with a unitary @xmath38 defines a general solution to eq .
( [ 12 ] ) . of course , eq .  (
[ 14 ] ) in itself is not enough to define an acceptable dirac operator . for the free case the overlap dirac operator has the properties ( a)-(c ) listed in the no - go theorem .
it will be argued later that @xmath2 remains local in a gauge field background as well . as a consequence
, @xmath63 defines an _ acceptable _ lattice dirac operator ( properties ( a)-(c ) of no - go theorem ) .
let me make a few simple comments .
( some of these straightforward observations in this or modified form have been made independently by several people and i shall not always give a reference .
see , however , refs .
@xcite . )
\(1 ) instead of @xmath54 in eq .
( [ 8 ] ) one could start with any acceptable dirac operator @xmath64 ( massless , no doublers , local ) .
\(2 ) if the original @xmath64 satisfies the gw relation eq .
( [ 12 ] ) then @xmath65 hence @xmath66 , i.e. the operator reproduces itself .
\(3 ) @xmath67 is hermitian and @xmath68 . as a consequence @xmath69 where @xmath70 is the sign function
\(4 ) for @xmath71 ( where @xmath72 ) one has @xmath73 which satisfies a slightly modified relation @xmath74 the choice of @xmath75 will be useful as discussed later .
note that for @xmath76 there are no massless fermions at all , while for @xmath77 there are too many , since the modes around @xmath78 ( the massive doublers ) also become massless modes of @xmath2 .
therefore only @xmath72 produces an acceptable dirac operator .
\(5 ) it is also easy to generalize to the case of arbitrary local @xmath31 : @xmath79 and @xmath80 this satisfies eq .
( [ 2 ] ) . for simplicity
, however , we shall consider mostly the case @xmath58 .
for a massless dirac operator in the continuum theory the atiyah - singer index theorem @xcite holds : @xmath81 here @xmath82 is the topological charge of the gauge field , and @xmath83 where @xmath84 , @xmath85 is the number of solutions to @xmath86 with @xmath87 chiralities , i.e. @xmath88 .    the spectrum of @xmath2 in the continuum consists of    * @xmath89 ; the corresponding modes have definite chirality * @xmath90 ( @xmath91 real , non - zero ) ; for the corresponding modes @xmath92 .
obviously , @xmath93 is topological invariant : an unpaired eigenvalue at @xmath89 can not be continuously produced or destroyed .    in the lattice formulation one
has problems with both sides of eq .
( [ 17 ] ) . on the lhs .
there are several ways to assign @xmath94 to a given gauge configuration ( field theoretical definition , geometrical definition ) . to avoid ambiguities for strongly fluctuating gauge fields ,
cooling is often applied .
an additional problem is that , in general , dislocations could be present , i.e. configurations which violate the inequality @xmath95    on the rhs . in an arbitrary gauge field @xmath54
has no exactly zero eigenvalues
they are distorted by lattice artifacts .
the spectrum of @xmath54 consists of    * real eigenvalues , with @xmath96 but not exactly @xmath97 * complex eigenvalues , in cc .
pairs , @xmath98 and @xmath99 with @xmath92 .
for smooth instanton background ( @xmath100 ) @xmath54 has an isolated real , nearly zero eigenvalue with @xmath101 or @xmath102 .
however , for a small instanton , or a typical mc configuration ( say at @xmath103 ) there are a lot of small real eigenvalues . therefore the definition of topological charge based on investigating the spectrum of @xmath54 is quite ambiguous @xcite .
let me remind you that the rg approach has a theoretically appealing ( although difficult in practice ) way to define the topological charge of a given configuration , in spin models @xcite and in gauge theories @xcite .
the solution is the following : take the minimizing configuration on the fine lattice @xmath104 and define the fp topological charge on @xmath105 through @xmath106 it is easy to show that with this definition there are no dislocations , i.e. @xmath107 holds for any configuration .    as has been shown in ref .
@xcite the fp fermion operator has nice chiral properties which allow one to establish the exact index theorem on the lattice .
let me discuss the main points .
the gw relation eq .
( [ 12 ] ) ( with @xmath108 ) and @xmath109 implies @xmath110 hence @xmath2 and @xmath111 commute , i.e. @xmath2 is normal .
( as a consequence , the eigenvectors corresponding to different eigenvalues are orthogonal to each other . ) the spectrum of @xmath2 ( cf .
( [ 14 ] ) ) is on the circle @xmath112 .
there are three types of eigenvalues    * @xmath89 , with @xmath113 ( denote their numbers by @xmath84 , @xmath85 ) , * @xmath114 , with @xmath113 ( with multiplicities @xmath115 , @xmath116 ) , * @xmath117 complex , with @xmath118 and @xmath119    note that @xcite since @xmath120 the index of @xmath2 and @xmath121 is opposite : @xmath122    define the topological charge density by @xmath123 for the topological charge @xmath124 we have @xmath125 here we have used the properties of eigenvectors of @xmath2 listed above .    obviously , @xmath126 is topological invariant since the integer @xmath93 can not change smoothly .
this fact can also be obtained directly from the gw relation .
equation ( [ 28 ] ) is the index theorem on the lattice : it states that one can define a ( gauge invariant , real ) topological charge density @xmath127 on a given gauge configuration which leads to an integer topological charge .
the topological charge defined this way equals to the index of @xmath2 .
both sides are defined through @xmath2 .
is this a tautology ?
no : the important point here is that one has a local topological charge density not just an integer number defining the charge .    it is easy to see that for smooth configurations @xmath128 since @xmath129 is a gauge invariant , pseudoscalar , dimension 4 operator and the prefactor is fixed by the fact that in the continuum limit @xmath35 one should recover the continuum topological charge .
i would like to stress that _ any acceptable _ solution to the gw relation defines a topological charge for which the exact index theorem holds .
note , however , that different solutions will give , in general , different values for the topological charge  they need to agree only for sufficiently smooth configurations .
from a practical point of view the operator @xmath2 has to provide a topological charge which is robust enough , giving a reasonable answer also for relatively rough configurations . to illustrate this point , consider neuberger s operator in eq .
( [ 15 ] ) with a small @xmath130 , say @xmath131 . for this choice
only very large instantons are counted as such when defined through eq .
( [ 28 ] )  perhaps only those with @xmath132 .
smaller instantons `` fall through the lattice '' since the smallest eigenvalue of @xmath54 would shift beyond @xmath130 . obviously , this would be a bad dirac operator and a bad topological charge  in spite of the fact that they satisfy the exact index theorem , eq .
( [ 28 ] ) .
any _ acceptable _ discretization of the dirac operator ( without requiring to be a solution to the gw relation ) could be used to determine the correct topological charge for _ sufficiently small _
lattice spacing @xmath13
the exactly zero eigenvalue of the continuum dirac operator will move only slightly away from zero , remaining real , and no other nearly zero eigenvalues will appear .
the extra requirement of gw relation forbids the eigenvalue to move away from zero for such configurations .
hence no extra conditions ( beyond acceptability and gw relation ) are required for the index theorem to hold . however , as pointed out above , these conditions do not guarantee that the topological charge defined this way is the ` correct ' one  except for sufficiently large instantons .
of course , there is no unique way to define the topological charge on the lattice , but for not - too - small instantons ( with a radius of a few lattice spacing , say ) the geometrical definition is natural . for practical reasons , one would like to have a gw dirac operator @xmath2 which reproduces the natural topological charge of not too large instantons as well .
i discuss these points perhaps more than necessary , because there is some controversy in the literature @xcite .    in ref .
@xcite it has also been shown that for the fp dirac operator one has in addition @xmath133 i.e. the topological charge defined by @xmath134 _ coincides _ with the old definition in pure gauge theory . in this case , for example , there are no dislocations , and one expects small lattice artifacts .
an important point is that @xmath18 is discontinuous in the gauge field .
the reason is that on the lattice one can connect continuously any two configurations , but @xmath135 has a jump on the path connecting , say , @xmath136 and 1 .    with the fp action
the origin of this discontinuity is that the minimizing configuration @xmath137 becomes ambiguous when the rhs . of eq .
( [ 3 ] ) has a degenerate minimum ( when the instanton `` falls through the lattice '' ) . for neuberger s dirac operator
the discontinuity appears when a real eigenvalue of @xmath54 passes 1 . in this case
an eigenvalue of @xmath138 changes sign , and @xmath139 in eq .
( [ 14b ] ) is discontinuous .
note that a similar discontinuity appears when the fermions are defined in the continuum , in a continuous gauge background obtained by interpolating the lattice gauge field @xcite . in this case
the process of interpolation becomes discontinuous , similarly to the situation with the fp action .
an important question : does @xmath2 become non - local when this discontinuity occurs ?
this would be a very bad news for the present approach .
fortunately , this is not true in general . in the fp action
the two solutions @xmath140 and @xmath141 corresponding to the degenerate minima are expected to differ only locally .
neuberger s @xmath2 , as shown by lscher @xcite , remains also local when an isolated eigenvalue of @xmath142 passes zero , so a discontinuity of d does not imply a loss of locality .
restoring the lattice spacing @xmath13 , the spectrum of @xmath2 is given by a circle @xmath143 with @xmath144 , crossing the real axis at @xmath89 and @xmath145 .
correspondingly , the spectrum of @xmath33 is a straight line parallel to the imaginary axis , of the form @xmath146 , @xmath147 for arbitrary gauge field .
denote the unnormalized expectation value of @xmath148 in a fixed gauge background by @xmath149 then we have @xmath150 obviously , the @xmath151 part of @xmath152 is independent of the dynamics and has to be subtracted to define a proper chiral order parameter .
this yields to @xmath153 the same as suggested in ref .
@xcite .    adding a mass term @xmath154 to @xmath2 would shift the whole circle of eigenvalues . since the real eigenvalues at @xmath145 also contribute to chiral relations it is better to leave them unchanged .
this is simply achieved by the replacement @xmath155 which not only shifts the circle but also rescales its radius appropriately . with this definition of @xmath154 one obtains @xmath156 this expression is identical to eq .
( [ 33 ] ) , i.e. the proper order parameter .
this equivalent form on the rhs . has been suggested first by chandrasekharan @xcite , who obtained it from lscher s exact symmetry .
at this moment we are in a peculiar situation : all signs of a symmetry are present ( ward identities , index theorem , chiral spectrum ,  ) but the action is not chirally invariant .
lscher made an important observation that in fact the action possesses an exact symmetry .
perform an infinitesimal change of variables @xmath157 and @xmath158 with a global flavour singlet transformation @xmath159 from the gw relation eq .
( [ 12 ] ) it follows that the action is invariant @xmath160 the corresponding flavour non - singlet transformation is given by @xmath161    now it seems that we have more symmetry than possible  the flavour singlet transformation should be anomalous . where is the anomaly hidden ? as pointed out by lscher ,
the fermionic integration measure is _ not _ invariant under the singlet transformation in eq .
( [ 35 ] ) @xmath162 & = & { \rm tr}(a\gamma_5 d ) [ d\overline{\psi}d\psi ] \label{37 } \\   & = & 2n_{\rm f}(n_- - n_+)[d\overline{\psi}d\psi ] \nonumber\,.\end{aligned}\ ] ] that is , the measure breaks the flavour singlet chiral symmetry in a topologically non - trivial gauge field where @xmath163 .
the corresponding ( global ) ward identity for the flavour singlet transformation is @xmath164 where @xmath165 .
note that the non - singlet symmetry is not anomalous , @xmath166 since eq .
( [ 37 ] ) for this case contains @xmath167 .    using these exact symmetries of the action ,
the ward identities follow , of course , directly . by taking @xmath168 and a non - singlet transformation one
has @xmath166 i.e. @xcite @xmath169 this is the order parameter which is expected to be broken spontaneously in the thermodynamic limit @xcite : @xmath170 in ref .
@xcite there is also a relation connected with the u(1 ) problem ( massive @xmath171 meson )  this will be discussed later .
whether the chiral symmetry is indeed broken spontaneously in qcd is a delicate dynamical question .
the chiral invariant dirac operator provides a firm basis to investigate this question by avoiding an explicit breaking .
note that the nielsen - ninomiya no - go theorem is avoided here by the non - standard form of the lattice chiral symmetry : @xmath2 is not invariant under the naive @xmath32 symmetry , nevertheless it is invariant under the generalized transformation .
in the continuum , chiral symmetry allows one to decompose the fermion field into left and right chiral components , transforming independently . besides yielding a convenient formalism for vector theories like qcd
, this decomposition is a necessary first step towards the chiral gauge theories .
the matrix @xmath32 anti - commutes with the continuum dirac operator and @xmath172 .
this allows to introduce the chiral projectors @xmath173 and the action decomposes as @xmath174 .    on the lattice
it is useful to introduce @xcite the operator @xmath175 where @xmath2 satisfies the gw relation eq .
( [ 12 ] ) .
remember that @xmath38 is unitary and @xmath176 . then @xmath177 one can introduce the projectors @xmath178 they satisfy the relations @xmath179 and allow the decomposition @xmath180 to conclude the definitions let us introduce the chiral components of @xmath181 : @xmath182 in terms of these components we have @xmath183 analogously to the continuum case .
note that eqs .
( [ 50],[51 ] ) are not symmetric in @xmath181 and @xmath184 . as discussed later , this asymmetric definition seems to play an essential role .
analogously to the continuum case , the mass term should couple the l and r components ( but not the l , l or r , r components ) : @xmath185 the corresponding pseudoscalar is given by @xmath186 the scalar ( pseudoscalar ) densities @xmath187 and @xmath188 will often appear in the formulae below .    performing a left chiral transformation
means @xmath189 using this transformation one has @xmath190    the combinations in eqs .
( [ 53],[54 ] ) are the natural generalizations of the corresponding continuum expressions .
note that the @xmath191 factor cancels the contribution of the unphysical real eigenmode at @xmath145 . using these combinations of fields ,
one obtains very simple , compact relations , for example : @xmath192 note that only the @xmath89 eigenmodes contribute when evaluating the trace .
it is instructive to consider the flavour singlet ward identity which leads to a relation for the @xmath171 propagator : @xmath193 if the @xmath171 stays massive in the chiral limit @xmath194 then the lhs .
( after averaging over the gauge fields ) should go to zero .
this leads to the interesting relation @xcite between the topological susceptibility and chiral condensate in full qcd @xmath195 this relation has been obtained previously in the continuum @xcite but it is important to have its lattice version as well .
it shows , for example that the topological susceptibility must be a well behaved physical quantity on the lattice .    for the flavour non - singlet variant of eq .
( [ 58 ] ) the term proportional to @xmath196 is absent since this symmetry is not anomalous .
therefore a non - vanishing chiral condensate , eq .
( [ 42 ] ) , yields to ( quasi ) goldstone bosons with mass @xmath197 .
we conclude this section with a few remarks .
note that @xmath198 coincides with @xmath139 introduced in the overlap formalism .
another point worth mentioning is that for the modes with @xmath89 one has @xmath199 while for the modes at @xmath145 they differ in sign , @xmath200 , as seen from the definition eq .
. this observation will be important when we consider chiral fermions .
we can now introduce the @xmath201 parameter , as in the continuum . here
i shall restrict the discussion to some basic points only .    define the lattice lagrangean with a complex mass @xmath154 and the @xmath201-parameter : @xmath202 where @xmath203 is the gauge action and @xmath129 is defined through @xmath2 by eq .
( [ 27 ] ) . by changing the variables
@xmath204 one can restore the invariance of the action by @xmath205 and the non - invariance of the measure can be compensated by @xmath206 .
therefore the free energy should depend only on the combination @xmath207 , as in the continuum .
one can proceed now e.g. by repeating ( now on the lattice ) the steps of ref .
@xcite where e.g. useful sum rules have been obtained for the spectrum of the dirac operator . to illustrate this ,
consider the partition function @xmath208 under some natural assumptions @xcite one has for the free energy @xmath209 taking then the second derivative with respect to @xmath201 ( for real @xmath154 ) one obtains @xmath210 which is equivalent to eq .
( [ 59 ] ) .    as another example , it is easy to show that the fermion determinant is given by @xmath211 for @xmath212 , and the product is over the complex eigenvalues only . for @xmath213 the prefactor
is replaced by @xmath214 .
( [ 67 ] ) differs from its continuum counterpart by the presence of the factor @xmath215 .
an important consequence of the exact chiral symmetry on the lattice is that there are no o(@xmath13 ) lattice artifacts , to any order of @xmath216 .
the reason is the following .
the o(@xmath13 ) lattice artifacts of any action could be compensated by adding an extra term which is a lattice version of @xmath217 with an appropriately chosen coefficient @xcite . the original action in our case has the exact ( generalized ) chiral symmetry of eq .
( [ 35 ] ) and its artifacts could be compensated only by a similarly invariant extra term .
any lattice regularization of @xmath217 will violate , however , this symmetry , hence the o(@xmath13 ) artifacts should be absent for the original lattice action as well .
note also that the clover term needed for o(@xmath13 ) improvement has been determined by requiring that some ward identity holds independently of the gauge field .
this is , however , automatically satisfied here for @xmath218 hence there is no need for extra clover term .
in addition , in the presence of a mass @xmath154 there are also no o(@xmath219 ) artifacts , if the mass is introduced properly , as in eq .
. then the symmetry @xmath220 forbids such term .
indeed , under a finite chiral transformation @xmath221 @xmath222 is invariant while the mass term @xmath223 changes sign .
although due to the anomaly the partition function is not invariant for @xmath224 , this fact is not relevant for our considerations .
indeed , for @xmath225 and sufficiently large volume one can restrict the system to the topologically trivial sector @xmath226 @xcite .
note , however , that the o(@xmath227 ) lattice artifacts could still be large at the relevant lattice spacings , and could be quite different for different solutions of the gw relation .
finally , it is interesting to see ( although it is not necessary after the previous discussion ) that the absence of o(@xmath13 ) artifacts on the tree level follows directly from the gw relation : expanding @xmath2 in powers of @xmath13 it implies that @xmath228 ( with some constant @xmath229 ) i.e. a term with @xmath230 is automatically generated .
here we shall discuss some basic properties of neuberger s dirac operator and its straightforward generalizations .
consider first the free case .
the spectrum @xmath231 of @xmath2 defined in eq .
( [ 11 ] ) is shown in fig .
[ fig : spect ] . together with the spectrum for the standard wilson - dirac operator for @xmath232 .
the latter is denoted by w and given by @xmath233 while n denotes the two branches for neuberger s dirac operator where @xmath234 is given by @xmath235 ( note that the lower curve extends only up to the point where the two branches join each other  this peculiar behaviour is due to the non - analyticity in the square root . )
the spectrum for @xmath236 is worse than that for @xmath54 .
this illustrates that a solution of gw relation can have large o(@xmath227 ) artifacts  the fact that @xmath2 is spread over some distance does not necessarily make the artifacts small !
the locality of @xmath2 for the free case is guaranteed by the inequality @xmath237 where @xmath238 and @xmath239 .
this asymptotic behaviour is shown on fig .
[ fig : loc ] by dashed line .
the decay of @xmath240 for neuberger s @xmath2 is shown by circles .
two more solutions to gw relation are added .
crosses show the case with @xmath241 where @xmath242 is a fp action truncated to the hypercube . since there is a non - trivial @xmath31 in this case the more natural choice , @xmath243 , ( cf .
( [ 16],[16a ] ) ) is also shown ( + signs ) .
finally , the true fp action is also indicated ( triangles ) .
it is obvious that there are rather large differences in the range of interactions for different gw dirac operators , and this is an important factor in making a choice for numerical simulations ( see also ref .
@xcite ) .
a less trivial point is the question of locality in a gauge field background . in ref .
@xcite it is proven that for small gauge fields @xmath244 , the corresponding @xmath18 is local with exponentially decaying tail , and the localization length is uniformly bounded .
although eq .
( [ 74 ] ) should hold sufficiently close to the continuum limit , it is violated at practical values of @xmath13 .
theoretically one can restrict the gauge fields to satisfy eq .
( [ 74 ] ) .
this is sufficient to show that there exists an exactly chiral invariant , acceptable lattice regularization of qcd .
however , the restriction leads to an astronomical correlation length , hence it is only of theoretical interest . as mentioned before , in ref .  @xcite it is also shown that an isolated zero of @xmath245 does not spoil locality .
an important question is the locality in gauge fields generated by mc simulations .
first results on this are presented in refs .
@xcite . in fig .
[ fig : f2r ] results from ref .
@xcite are shown , for su(3 ) gauge fields generated at @xmath246 , @xmath247 and @xmath248 ( triangles , squares , circles ) .
the full symbols are for @xmath249 while the open ones for ( roughly ) optimized values , @xmath250 for @xmath246 and @xmath251 for @xmath252 .
it is also important to have a sufficiently effective way of evaluating the inverse square root .
this issue is discussed in refs .
@xcite .
i only mention here that it seems to be a feasible task but still much has to be done in this direction .
a few comments are in order :    \(1 ) unlike evaluating @xmath33 for massless quarks , taking the inverse square root is not a critical problem  appearance of zero modes of @xmath245 is accidental .
\(2 ) starting with a good approximation of the fp action instead of @xmath54 ( as in eq .
( [ 16 ] ) ) the probability to have a zero mode of @xmath245 is strongly suppressed , since with the exact fp action @xmath245 would be unitary .
note , however , that for any continuous parametrization of the fp action , @xmath245 must have a zero mode for specific configurations , when an instanton `` falls through the lattice '' .
this happens on the boundary separating regions of gauge field configurations with different topological charges .
\(3 ) using a gw dirac operator there will be no `` exceptional configurations '' in the usual sense @xcite  these would correspond to eigenvalues @xmath253 .
( one should distinguish these from configurations where @xmath254 appearing in this approach . )
\(4 ) one can construct , in principle , configurations when many `` instantons fall through '' simultaneously at different sites of a large lattice .
this would produce a situation when @xmath245 has many nearly zero eigenvalues .
will @xmath2 be local on such configurations ?
one expects that such configurations have negligible weight in the partition function .
will their effect be also negligible ?
i think , these questions have to be investigated .
there are several equivalent ways to measure the topological charge defined by @xmath2 through eqs .
( [ 27],[28 ] ) :    \(1 ) find the zero eigenmodes of @xmath2 ( or of hermitian @xmath255 ) ,    \(2 ) count the real eigenvalues of @xmath54 in the region @xmath256 , with chirality signs ,    \(3 ) trace the `` level crossings '' of @xmath257 while @xmath258 runs from @xmath259 to @xmath130 .
note that although we rely here on the exact index theorem on the lattice , the value of the topological charge determined through @xmath2 with @xmath260 is the same as the method using @xmath54 suggested a long time ago @xcite .
the last method has also been used in ref .
@xcite in the framework of the overlap formalism .
this definition of the topological charge works well for small instantons , as expected . but how does it work for typical mc configurations ?
does @xmath2 ( generated using @xmath54 ) give a ` better ' definition than the geometrical one ?
hernandez @xcite argues that the answer is no .
gattringer and hipp @xcite simulated su(2 ) on a @xmath261 lattice at @xmath262 and found that the value of the topological charge , @xmath263 on individual configurations depends strongly on the parameter @xmath130 since there are usually many real eigenvalues in the region of interest , and only very few configurations have a gap in the spectrum where @xmath263 does not depend on @xmath130 chosen in a reasonable range .
they also note that a naive inclusion of the clover term makes the situation worse . the same conclusion is found by degrand , a.  hasenfratz and kovcs @xcite .
this last fact is not too surprising since the corresponding artifact is not an o(@xmath13 ) effect .
note , however , that in the average this definition of charge works better : edwards , heller and narayanan @xcite found that the topological susceptibility stays roughly constant ( independent of @xmath130 ) for @xmath264 . here
i would like to point out that the @xmath130-independence of @xmath265 does not necessarily mean the ` true ' value : if @xmath266 then @xmath267 . a good definition of the topological charge ( at a given lattice spacing ) should give relatively unambiguous definition of @xmath263 on individual configurations .
in other words , the value of the lattice spacing @xmath13 when a gw dirac operator gives a ` good ' definition of the topological charge could strongly depend on the choice of @xmath2 . with the presently used lattice spacings the original choice of neuberger ( using @xmath54 ) seems to produce a topological charge which is too ambiguous .    on the other hand
, measuring the topological charge with the fp action in the schwinger model @xcite gives quite convincing results . for 4d gauge theories ,
however , the situation is less satisfactory @xcite  mostly because of technical problems to perform the required multigrid minimization on the gauge configurations .    for qcd a good strategy would be a compromise : take a dirac operator @xmath64 with good properties ( small cut - off effects , nearly satisfying the gw relation ,  ) and build @xmath2 from it , which in turn could be used in the definition of the topological charge .
it would have smaller o(@xmath227 ) artifacts , better locality , fewer ` dangerous ' configurations , larger gap in the real eigenvalues hence less ambiguity of @xmath126 , but will also need more programming effort and extra computational overhead .    measuring
the topological charge density using @xmath268 is perhaps too time consuming to be feasible .
the question whether it is possible to define chiral gauge theories ( like the standard model ) on the lattice is a purely theoretical problem , but of great importance .
much has been done in this direction and has been reported earlier @xcite .
nevertheless , there remain several open questions .
a gw dirac operator seems to be a good starting point : it has exact gauge invariance and one can define weyl fermions ( cf .
( [ 50],[51 ] ) ) . here
i shall report briefly about new developments ( partially yet unpublished ) by lscher @xcite based on the properties of gw dirac operators .
the discussion will touch only a few interesting points which could be easily understood on the basis of previous sections .
consider fermion fields which are purely left handed ( see eq .
( [ 50 ] ) ) @xmath269 these are local , gauge invariant conditions .
introduce a basis for @xmath270 and @xmath271 through relations @xmath272 where @xmath273 and @xmath274 are grassmann variables . in this basis
the fermion action reads @xmath275 with @xmath276 .
the fermion propagator ( when there are no zero modes of @xmath2 ) is given by @xmath277 it is easy to show that this is chiral in the generalized sense .
note that @xmath278 the first term on the rhs .
is chiral in the usual sense ( with @xmath32 ) while @xmath279 is a contact term .
consider a chiral gauge theory with fermion multiplet @xmath280 , with integer charges @xmath281 .
for the number of zero modes @xmath282 , @xmath283 the index theorem implies @xmath284 since @xmath285 is replaced by @xmath286 for smooth gauge fields when @xmath287 .
a simple but important observation is that the number of @xmath181 and @xmath184 fields are not necessarily equal .
they are given by @xmath288 and @xmath289 , respectively .
for their difference we have @xmath290 to illustrate this , consider a background field with @xmath291 , where @xmath292 , @xmath293 and @xmath294 , @xmath295 .
then for the @xmath89 mode @xmath296 while for the eigenvector @xmath297 at @xmath114 @xmath298 according to eqs .
( [ 76]-[78 ] ) , at @xmath89 we have only @xmath181 ( but not @xmath184 ) while at @xmath114 both @xmath181 and @xmath184 are present since @xmath32 and @xmath198 differ by sign at this point . obviously for complex values of @xmath98 both @xmath181 and @xmath184
are present again .
altogether there is one more @xmath181 than @xmath184 .
of course , this statement depends on the topological charge of the gauge field .    as a consequence the fermionic average @xmath299 of any operator in a given background gauge field vanishes unless the fermionic number of type @xmath91 ( the number of @xmath300 s minus the number of @xmath301 s ) in @xmath302 equals to @xmath303 .
perform now ( for the u(1 ) case ) a global gauge transformation @xmath304 .
this modifies @xmath302 by a phase factor @xmath305 where @xmath306 . from the other side a global u(1 )
gauge transformation does not change the gauge fields at all , hence invariance under this transformation requires @xmath307 i.e. the condition of anomaly cancellation of the continuum theory .
these observations , however , constitute only the starting point of a difficult mathematical problem  to find a satisfactory choice of the fermionic integration measure ( `` chiral determinant '' ) , i.e. of the basis @xmath308 .
lscher has proven that it is possible to construct a basis for the u(1 ) case that satisfies the following requirements :    \(1 ) the fermionic expectation values @xmath309 are smooth functions of @xmath17    \(2 ) the field equations ( obtained by some change of variables ) are _ local _ : @xmath310 if @xmath20 is far from @xmath311 .
\(3 ) @xmath309 is gauge invariant if @xmath302 is .
note that the last property may not be necessary to require . in principle , it is possible that after the leading violation of gauge invariance is cancelled by the condition @xmath312 , the rest is small and goes away in the continuum limit , without the need of extra counter terms .
this possibility is based on a simple but surprising result of foerster , nielsen and ninomiya @xcite , who show that gauge invariance on the lattice is robust in the sense that the effect of a sufficiently small violating term vanishes in the continuum limit .
the constructions of the overlap and gauge fixing @xcite formalisms assume this scenario .
it is clear , however , that it is much safer  if possible  to preserve chiral gauge invariance at finite lattice spacing .
it has been illustrated that dirac operators satisfying the gw relation provide an elegant solution to the chirality problem in lattice qcd .
in particular , there is no additive quark mass renormalization , and ward identities , soft pion theorem , etc .
are valid , there are no o(@xmath13 ) artifacts .      for practical implementations of these ideas
new algorithms are needed . that the task is feasible is illustrated by mc calculations in the domain wall approach @xcite which could be considered as a particular approximation to a gw dirac operator .
not restricting itself to the given ideology could open ways to more efficient algorithms .
t.  degrand , a.  hasenfratz and d.  zhu , nucl .
b475 ( 1996 ) 321 ; nucl .
b478 ( 1996 ) 349 ; nucl .
phys . b(proc.suppl . ) 53 ( 1997 ) 945 .
+ t.  degrand , a.  hasenfratz and t.  kovcs , nucl .
phys . b(proc.suppl . ) 63 ( 1998 ) 528 .


geomagnetism is still one of the unsolved fundamental problems of the earth since william gilbert discovered the magnetized earth ( gilbert 1600 ) .
not only the very existence of the dipole magnetism , motonori matsuyama discovered the polarity reversal events ( matsuyama 1927 ) and vivid dynamics of geomagnetism .
the qualitative nature of this polarity flip was explained by a simple disk model by tsuneji rikitake ( rikitake 1957 ) .
recent numerical simulations of magnetohydrodynamics ( mhd ) for geomagnetism ( roberts & glatzmaier 2000 ) , ( kono & roberts 2002 ) seem to describe such dynamics .
inspired by these simulations , we recently studied to clarify the physics behind this remarkable polarity reversals proposing a simple model ( domino model ) ( mori , schmitt , ferriz - mas et al .
2011 ) composed from many macro - spins which are interacting with each other .
this _ coupled spin model _ is based on the idea that the whole dynamo mechanism is described by global interactions of many small dynamo elements ( macro - spins ) .
this model naturally yields many of the observed features of geomagnetism : its time evolution , the power spectrum , the frequency distribution of stable polarity periods , etc .
( mori et al . 2011 ) . in the case of the earth , the dynamo element , that a macro - spin describes ,
is considered to be the taylor cell in the iron fluid core produced and supported by the coriolis force .    in this paper
, we would like to develop the above idea by elaborating the domino model in a form that is _ minimal and general_.    in order to construct the _ minimal _ model , we introduce a _ long - range coupling _ of macro - spins , i.e. all the spins interact with each other with the same amplitude of interactions .
this is essentially different from the previous domino model in which only neighboring spins interact with each other .
it turns out that this long - range coupling model is also effective to describe the observed features of geomagnetism as the domino model .    in order to show the _ generality _ of the model , we consider it in wider contexts .
the spin model , so far , is at most a phenomenological model to describe geomagnetism .
therefore it needs some _ supplement from the basic mhd equations _ for better descriptions . on the other hand ,
magnetic fields associated with celestial objects are quite common in the universe ; many planets , satellites , and the stars including our sun .
although the magnetic fields of these objects have variety that is quite different from the geomagnetism , all the magnetic fields are thought to be generated and supported by the robust dynamo mechanism .
therefore we try to fit our model , supplemented with mhd equations , in _ other celestial objects _ and see how they are consistent with each other .
the long - range coupled spin model can also describe the _ synchronization _ physics as well if we simply change the parameter value in the model .
this property is particularly useful to describe the quasi - periodic variation of solar magnetic field .
furthermore the long - range coupling in our model itself can produce sufficient chaoticity and randomness without introducing explicit random force as in the domino model .
thus the long - range coupling spin model has rich physics and generality .
we start our study by the basic description for geomagnetism .
the _ fluid iron _ in the center of the earth is now considered to be the place where geomagnetism is created ( inglis 1981 ) . the conductive fluid motion and its interaction with electromagnetic fields should control the dynamics of geomagnetism through dynamo mechanism .
they are described by the _ magneto - hydrodynamics _ ( mhd )  for incompressible ( @xmath0 ) fluid , @xmath1@xmath2@xmath3@xmath4with appropriate boundary conditions ( roberts & glatzmaier 2000 ) .
the dynamics of fluid velocity @xmath5  in eq.([eqforv ] ) is governed by the nonlinearity @xmath6 the pressure @xmath7 the centrifugal force potential @xmath8 , gravity @xmath9 , and the lorentz force @xmath10 where @xmath11 is the electromagnetic current and @xmath12 is the mass density .
the coriolis force @xmath13 yields the vorticity and convective pattern of the fluid .
the magnetic field @xmath14 is amplified or reduced by the flow @xmath15 as described by eq.([eqforb ] ) .
the dissipations ( @xmath16 ) and the energy injections ( from the heat generation @xmath17 and possibly the rotational driving force )  should balance with each other in the stationary state .
main non - dimensional quantities which characterize the above set of equations , and their typical values in the iron fluid core in the earth , are as follows:@xmath18@xmath19@xmath20@xmath21these values are important to estimate the scale of dynamo elements in later sections .
straightforward approach based on these highly complicated equations requires sophisticated numerical calculations .
such analysis have recently been developed and the generation , maintenance and even reversals of geomagnetism have been obtained ( kida & kitauchi 1998),(roberts & glatzmaier 2000),(kono & roberts 2002 ) . these works have provided us with valuable knowledge to develop fundamental understanding of geomagnetism .
however the parameter range in numerical calculations is still far from the above realistic values
. therefore it would be very favorable if our coupled spin model is complementary to the mhd calculations and if it describes some basic features of dynamo mechanism .
these are the aim of this paper .
the plan of this paper is as follows . in _ section 2 _
, we describe our idea of the coupled dynamo elements which plays the central role in the present paper .
we emphasize the necessary condition for the elements in order to clarify what the macro - spin represents in our model . readers who immediately want to see
the model can skip this section . in _ section 3 _
, we briefly review our previous study on the domino model ( short - range coupled spin model ( scs ) ) , which is a first realization of the idea of the coupled dynamo elements to describe geomagnetism . in _ section 4 _ , we introduce the generalized long - range coupled - spin model ( lcs ) in which macro - spins have long range interactions
. then we apply this model to geomagnetism and successfully describe many of the characteristic features of geomagnetism . in _ section 5 _
, we study the physical relevance of lcs model showing the similarity of it to the mean field model ( hmf ) and kuramoto model .
readers who immediately want to see further applications of the model can skip this section . in _
section 6 _ , we apply our model , supplemented with mhd equations , to other planets and satellites in order to estimate the scale of macro - spins . in _ section 7 _
, we apply the same model , with different parameter values , to the solar magnetism and successfully describe the quasi - periodic dynamics of the solar magnetism .
the last _ section 8 _ is devoted to summarize our work , and speculate possible generalizations of our idea of coupled dynamo elements .
we now describe the idea of coupled dynamo elements ( mori et al .
2011 ) which plays the fundamental role in this paper .
the central equation to describe the dynamo mechanism in mhd is eq.([eqforb ] ) , in which the flow @xmath22  is determined by eq.([eqforv ] ) .
this is practically a very complicated set of equations to describe geomagnetism and the dynamo mechanism though we can extract some basic insight from it .
geomagnetism is characterized by the coexistence of many time scales , from million years to thousand years , with clear power - law power spectrum .
the dynamo mechanism is ubiquitous in the sense that most of the celestial objects have magnetism generated from this mechanism .
the ubiquitousness implies the basic dynamo mechanism is simple and the many time scales with power - law power spectrum imply the dynamo system is composed from many components .
therefore we expect that _ the whole dynamo system may be composed of many elements , each of which has minimal dynamo function , that interact with other elements in a simple form_. it is often happen that a simple interaction yields very complicated structures .
then what is the element which has the minimal dynamo function in the case for geomagnetism ?
we notice that the above eq.([eqforb ] ) is similar to the general equation for vorticity @xmath23 @xmath24which is derived by taking the rotation of eq.([eqforv ] ) , neglecting the lorentz force term .
this fact reveals the apparent duality of the magnetic field @xmath14 and the vorticity @xmath25 .
the iron fluid core is estimated to be highly turbulent from eq.([re ] ) and is full of vorticities .
the coherent vorticity @xmath25  is generated from the convection and coriolis force .
these effects seem to be significant from the values in eqs.([ta]-[ra ] ) .
therefore the iron fluid is full of vorticity , and possibly of magnetic fields .
it is clear that the coriolis force dominates in the iron fluid as in eq.([ta ] ) . in general ,
if the coriolis force dominates in almost the stationary flow , then eq.([eqforv ] ) becomes @xmath26 by taking the rotation of this form , we have @xmath27 .
therefore the flow has a tendency to be two - dimensional perpendicular to the rotational axes @xmath28 , _
i.e. _ the taylor - proudman theorem
. therefore this theorem suggests that a convection of iron fluid forms convective columns ( usually called taylor cell ) parallel to the rotation axes .
the scale of the taylor cell will be determined by geometry and various parameters in eqs.([re]-[ra ] ) .
thus coherent vorticity @xmath29  is naturally expected in the earth core .
however the vorticity @xmath25  is not sufficient to produce magnetism @xmath14 ( tsinober 2007 ) .
this is because the circular vorticity itself is not endowed with energy production .
any deviation from the complete circular motion would be necessary for the production of magnetic energy .
we would like to emphasize that the inward - winding vorticity is essential for dynamo mechanism .
suppose that we prepare a typical inward - winding solution for the navier - stokes equation , the burgers vortex , whose velocity field is expressed in the cylindrical coordinate @xmath30 as    @xmath31    with some constants @xmath32 , and @xmath33 for inward - winding and @xmath34 for outward - winding . then putting the above eq.([bergers ] ) into eq.([eqforb ] ) yields , in the cartesian coordinate ,    @xmath35    for the magnetic field @xmath36 parallel to the symmetry axes of the vorticity . only the first term in rhs can definitely enhance @xmath37 if @xmath38 the inward - winding flow squeezes the magnetic force lines against the positive pressure from @xmath37 .
summarizing all the above , the iron fluid in the earth core may allow _ multiple vortex columns parallel to the earth rotation axes and their interaction may cooperatively yield the whole geomagnetism , if the columns are inward - winding_.    on the other hand in the super computer simulations , such column structure is commonly formed ( kida & kitauchi 1998 ) and robustly persists .
these columns are called taylor cell(tc ) . a close look into tc in the numerical data , reveals that tc are either inward - winding ( _ i.e. _ high - pressure due to the compression ) or outward - winding ( _ i.e. _ low - pressure due to the spread ) .
they are called , respectively , anti - cyclone and cyclone .
that is the anti - cyclone counter rotates to the earth and the cyclone co - rotates .
the coriolis force associated with the earth rotation makes these difference and the smooth flow may favor the alternating configuration of cyclone and anti - cyclone .
according to the above argument , only the anti - cyclone , which has inward - winding flow , can produce magnetic fields ( kageyama & sato 1997 ) . in the real situations in the earth ,
the flow is highly turbulent and non - linear on top of such tc structures reflecting the large reynolds number eq.([re ] ) .
complicated magnetic fields in the outer core are captured by inward - winding tc and aligned with the rotation axes , and then compressed to enhance the poloidal magnetic fields .
these non - linear flow , as well as huge electric current , yields the interaction among tc , either short - range or long - range .
each dynamo element of inward - winding vorticity , anti - cyclone , accompanies generated magnetic fields which can be characterized by a direction and a strength in the first approximation .
therefore we characterize each dynamo element as a _ vector _ quantity @xmath39 .
this vector represents the direction of the magnetic fields associated with tc and not the direction of tc itself .
since each element is supposed to have the minimum dynamo function , this vector can be identified as a macroscopic spin or magnetic moment .
then the whole system should be a set of such macro - spins @xmath40 , which is located on each site @xmath41 in the iron fluid core .
although the interactions between the macro - spins should be very complicated , the potential energy of it must be a _
scaler _ or the vector inner product of the spins @xmath42 .
this simplicity may guarantee the robustness and universality of the geomagnetism .
however we can not determine the amplitude of the coupling parameter @xmath43 by a simple argument .
it depends on at least the iron flow pattern , electric currents , and configuration of cyclone and anti - cyclone . on the other hand
we can guess the signature of @xmath44 .
electric currents are smoothly winding each tcs , and cyclone and anti - cyclone align alternately .
therefore the smooth electric current yields the winding direction is the same for all anti - cyclones ( and opposite direction for all cyclones ) .
this suggests the magnetic fields of anti - cyclones themselves ( and also cyclones themselves ) have a tendency to align to the same direction .
this suggests the _ _  negative value for the coupling strength : _ _ @xmath45 . ] , like in the ferromagnetism in ferromagnetism has the quantum origin , which has nothing to do with our macro - spin model . ] . then the effective  curie point
 of these macro - spin system can be high enough and therefore , the whole system may actually show  ferromagnetism  .
this makes a quite contrast to the case of the ( micro ) spins of the iron in the earth ; no ferromagnetism is produced in the core of temperature @xmath46k which well exceeds the curie point about @xmath47k .
we have emphasized the necessary condition for the macro - spin @xmath48to be relevant and their interaction although the construction of the macro - spin @xmath48 is not our aim in this paper .
we would like to study the inter - relation between many of such elements @xmath48  which yields various non - trivial dynamics .    before our study
, there have been some works ( mazaud & laj 1989 ) , ( seki & ito 1993 ) , ( dias , franco , & papa 2008 ) in which the spin models with the interaction @xmath42 were studied to reveal the polarity flip dynamics of geomagnetism .
they are based on the two - dimensional ising model and have found the polarity flip near the critical temperature . our aim to study geomagnetism in this paper
is not simply to extend the similarity with phase transitions but to study wider point of view such as _ phase coexistence _ and _
synchronization_.    there are two types of coupling between the spins according to the range of the interaction
. if the fluid motion is responsible to the coupling , then only the neighboring spins interact with each other and yield short - range coupled - spin ( _ scs _ ) system . on the other hand ,
if the electric current or magnetic fields are responsible to the coupling , then all the spins interact with the same strength and yields long - range coupled - spin system ( _ lcs _ ) .
these are the two extreme ends of the series of models with finite interaction rage .
we have already studied scs model in our paper ( mori 2011 ) , which successfully describes many characteristic features of geomagnetism as we now review in the next section before we study lcs model is the subsequent sections .
we now briefly review our previous model for geomagnetic dynamics ( mori et al . 2011 ) , in which the concept of coupled dynamo elements is realized as short - range interacting macro - spins ( scs model ) .
our intention was to establish a simple model which realizes various basic properties of geomagnetic observations based on the idea that the coupled dynamics of dynamo elements .
we assume that each element has dynamo function to produce element magnetism putting the very origin of the dynamo mechanism aside .
we first introduce the scs model , which is described by the lagrangian @xmath49 where@xmath50@xmath51and by appropriate dissipation and fluctuation terms .
the spins are assumed to be located on a circle with equal separation ( therefore @xmath52 ) .
this circle is set on the plane which includes the equator and fully immersed inside the fluid core region of the earth .
the spin @xmath53 at the cite @xmath41 may be represented by a single angle parameter @xmath54 as @xmath55 .
we measure the angle @xmath56  from the angular velocity vector @xmath57 . ) .
the angular velocity @xmath58 points to the top .
] this plane the spin @xmath48 moves is set perpendicular to the radius toward the cite @xmath41 .
the full evolution equation then becomes the stochastic differential equation ( langevin equation),@xmath59for @xmath60 .
the parameter @xmath61 is the dissipation coefficient .
the fluctuation force @xmath62 is assumed to be gaussian white noise as usual , characterized by the correlation:@xmath63then we have , so far , four parameters in our model @xmath64 .
the parameter @xmath65 represents the tendency for each spin to be parallel to the rotational axes .
larger value of @xmath65 means the dominance of the coriolis force over the other effects and the taylor - proudman theorem better holds .
the parameter @xmath44 measures the strength of interaction between the element spins .
it reflects the nature of the fluid flow and electric current across the neighboring spins .
the parameter @xmath61 represents the whole energy dissipation @xmath66 and @xmath67 in eqs.([eqforv],[eqforb ] , [ eqfort ] ) .
the parameter @xmath68 represents the inhomogeneous heat generation from the inner iron core and the random perturbation from the other element spins .
we also have to specify the number @xmath69 of spins .
it may depend on the geometry of the fluid core region of the earth in which cyclones and anti - cyclones are formed .
we can save the number of parameters if we focus on the dynamics in long time scale and neglect the inertial term @xmath70 .
then the reduced equation of motion becomes ( @xmath60),@xmath71which now has three parameters @xmath72 .
the full equation eq.([eomlcs ] ) is solved in ( mori et al .
the neglect of the inertial term is justified by the fact that both the methods yield almost the same results .
each dynamo element represented by the spin @xmath48 naturally corresponds to a magnetic dipole located at the site @xmath41  of the dynamo element .
we are interested in the order parameter @xmath73 of the set of equations eq.([redlcs ] ) defined by @xmath74this quantity @xmath73(@xmath75 $ ] ) is a good indicator of the whole magnetic field when we analyze the history of the geomagnetism though it is not the magnetic field itself .
it is also possible to describe the realistic configuration of the magnetic field on the earth by using the set of @xmath76 , ( @xmath60 ) .
as we have explained before , the magnetic moments @xmath48 ( @xmath60 ) are located on the ring with equal separation , on the surface which includes the equator .
then the whole magnetic field @xmath14 becomes @xmath77where each spin @xmath48 is located at the position @xmath78 , and @xmath79 is the unit vector @xmath80 .
this simple model naturally describes many of the observed properties of geomagnetism including the intermittent polarity reversal , power - law power spectrum , chron interval distribution , etc .
the results are summarized in the paper ( mori et al .
the short - range coupled spin ( scs ) model reviewed in the last section was simple and effective to describe the statistical properties of geomagnetism .
we now study a slightly modified model introducing the long - range interaction coupled spin ( lcs ) model .
the simplest one would be described by the same lagrangian @xmath81 with eq.([k ] , [ v ] ) but the interaction of macro spins is long ranged , in the denominator in the coupling term is simply from a technical reason so that all the terms in @xmath82 formally becomes additive . ]
@xmath83    according to the naive extension of scs , the full evolution equation seems to become ( @xmath60),@xmath84as in the previous section , introducing the fluctuation force @xmath62 and the dissipation term @xmath85 .
however , thanks to the long - range coupling , these _ fluctuation - dissipation terms are not necessary in lcs _ model
. the long - range coupling by itself can yield sufficient chaoticity which is necessary to describe the intermittent polarity flip .
it is still possible to reserve the fluctuation - dissipation terms , which does not alter the characteristic feature of the system .
therefore we would like to choose the simpler description which has no fluctuation - dissipation terms .
then the system becomes conservative system in which the total energy is conserved .
the evolution equation becomes simple ( @xmath60),@xmath86which has only two parameters @xmath65 and @xmath44 through the potential @xmath82 as eq.([vingcs ] ) .
the order parameter and the magnetic fields are described by the same expressions eqs.([m ] , [ b ] ) as in the scs model .
we next report the numerical results of lcs model with comparison with observations .
we have solved the set of equations for @xmath87 spins which are located on a circle .
we have set parameters of the model as @xmath88 , all of which is chosen to be of order one and we did not perform systematic fine tuning of them .
further we set the initial condition so that all the spins are rest and their direction angles are equally separated within the width of angle @xmath89 .
the reason that we did not distribute the initial spin angles in the full range @xmath90 $ ] is because we had to choose low energy to make the coexistence possible of the cluster and expanded spins , as we will explain in the next section .
time evolution of the total magnetic moment projected on the rotation axes @xmath91 is shown in fig.(fig4 - 1 ) top .
the end time @xmath92 is chosen so that the total number of polarity - reversal becomes several hundreds as is observed number @xmath93 in the history of geomagnetism within @xmath94 year in the past .
the total number of the polarity reversal is @xmath95 in this calculation .
therefore the unit calculation time ( @xmath96 ) corresponds to about @xmath97 year .
we have obtained the following results .    1
.   unit calculation time : @xmath98 @xmath97 year   2 .
average flipping time : @xmath99 yeas .
( observation @xmath100year ) 3 .
average chron time : @xmath101 @xmath102 year .
( observation @xmath103 year . ) 4 .
superchron : @xmath104 year ( observation @xmath105 year ) + the power spectrum of the time series @xmath106 is shown in fig.([fig4 - 1 ] ) bottom .
5 .   the power index : @xmath107 and @xmath108 , respectively , for right and left of the knee .
( observation : index @xmath109 in high frequency regime @xmath110 , and is @xmath111 in the low frequency regime @xmath112 . ) 6 .
the location of the knee : @xmath113 year ( observation : @xmath114my ) + these values are not much changed for the projected time series @xmath115 because the original time series @xmath116  is already almost the same as the projected shape , relatively long steady chrons and sharp reversal .
the frequency of chron intervals is calculated . the power index : @xmath117 ( observation : @xmath118 )
next we pick up a typical chron and calculated the power spectrum within this period .
+ the power index : @xmath119 for @xmath120year  ( odp : @xmath121 for @xmath122 )     in the numerical calculation of lcs model .
the positive and negative values in the vertical axis correspond , respectively , to the normal and the reversed polarity .
( botom ) the power spectrum of this time series . the data is sucessfully fitted by the the two power laws with the indeces -0.5 for low frequency and -1.7 for high frequency regions . , width=453 ]
. the left - right symmetry , or the symmetry in normal - reversed polarites , is not respected .
this indicates the model has strong fluctuation in the long time variation .
this accord with the ever increasing power spectrum toward low frequency region fig.([fig4 - 1 ] ) .
, width=453 ]    the apparent asymmetry in fig.([fig4 - 3 ] ) is not caused by the initial condition but suggests long range temporal fluctuations in our model .
this seems to be consistent with the ever increasing power spectrum in the low - frequency direction in fig.([fig4 - 1 ] ) .     in some chron .
the valuea are mostly negative but perpetually fluctuating strongly .
( botom ) the power spectrum of this time series .
the data is sucessfully fitted by the the two power laws with the indeces -0.46 for low grequency and -2.0 for high frequency regions . , width=453 ]     for each spin with shifted zero points for visibility .
the first spin from top was rapidly rotating before the flip but the nineth spin took over this role after the flip .
these perpetually rotating spins form the gas phase .
see the next section for detail .
( bottom ) the countour / density graph for horizontal component of the magnetic field eq.([b ] ) on the earth surface around this polarity flipping era .
each ellipse represents the entire earth surface drown in mollweide mapping , before , right in the middle , and after the flip from left to right .
the radius of the spin ring is assumed to be one - third of the earth radius .
, width=453 ]    the results of lcs model for geomagnetism is almost the same as our previous study on scs model and are not so far from the observational values .
thus the range of interaction is shown to be irrelevant for geomagnetic characteristics . in the scs model ,
even the inertial term was not necessary to yield the similar results .
it is interesting to notice that this lcs system , we explored in the last section , is a slight extension of the _ hamiltonian mean field model _ ( hmf ( antoni & ruffo 1995 ) , ( campa , dauxois , & ruffo 2009 ) by adding the potential term @xmath123 . hmf model a good tool to describes phase transitions and statistical fluctuations of a deterministic system composed from many elements .
slightly modified version of hmf model was applied to describe the self gravitating systems(sota , iguchi , morikawa , tatekawa , & maeda 2001 ) .
however this hmf model itself can not describe the polarity flip of geomagnetism because of the lack of the preferred angle .
the polarity flip is made possible simply by introducing the above potential term . in this sense
our lcs model is the _ minimal model _ to describe geomagnetism .
acute readers might have some concern that the geodynamics is a dissipative system rather than the conservative system like hmf model .
however the difference is not relevant in our context .
hmf and lcs systems have strong chaoticity and the mean field @xmath124 has apparent statistical fluctuations . in short , these systems have the relevant central degrees of freedom ( _ i.e. _ mean field ) and the rest ( _ i.e. _ environmental degrees of freedom ) together .
we focus on the mean field . the mean field immersed in
the rest behaves as the dissipative system .
this point will be clarified in the rest of this section .
dynamics of both hmf and lcs models can be visualized by the interacting @xmath69-particles which are moving on a circle of unit radius as shown in fig.(fig5 - 1 ) .
this is a typical snapshot of the spin angles @xmath125 represented by the points on a circle .
the whole ring represents all the directions of spins @xmath126 and the top position of the ring is the direction of the earth rotation axes @xmath127 @xmath128 .    .
the whole ring represents @xmath129 and the top position of the ring is @xmath130 .
the earth rotation axes @xmath28 has the fixed angle @xmath130 . as is shown the most spins are condenced to the almost fixed angle ( i.e. fixed polarity ) , and few rapidly rotating spins exit .
the latter trigges the flip of the condensed spins intermittently .
, width=453 ]    furthermore it is remarkable that both hmf and lcs models have two classes ( phases ) of elements : the spins in _ expanded phase _ which are almost freely moving all the angle and the spins in _ condensed phase _ which are bounded in a localized angle .
there exists a critical energy in hmf model in the limit @xmath131 ; the condensed phase appears only when the system energy is below this critical energy and the portion of the condensed phase increases for lower energy .
similarly in our lcs model , expanded and condensed phases coexist in wide rage of parameters is not so large . ] .
we have chosen our parameter so that the system energy is low enough to guarantee that the both phases coexist in one system .
the condensed phase yields finite order parameter defined by eq.([m ] ) .
free spins in the expanded phase perpetually disturb the condensed cluster spins and is considered to trigger the intermittent signature change of the order parameter eq.([m ] ) .
it is interesting that the expanded and condensed phases often exchange their constituent spins when the angles of the spins coincide with each other , _ i.e. _  collide  in the above visualization .    in terms of geomagnetism ,
coexistence of the expanded and condensed phases are essential ; the latter yields the well defined polarity in magnetism eq.([m ] ) and the former triggers the intermittent flip of the polarity .
this phase coexistence is naturally realized in lcs model when we choose mildly low energy values .
thus , this model naturally deduces that the system energy required for intermittent polarity flip necessitates _
the coexistence of the relatively steady component and the rapidly changing component in geomagnetism_. the former is of course from the condensed phase and the latter from the expanded phase .
since the macro - spin has magnetic dipole , these freely moving spins in the expanded phase yields _ rapidly moving pairs of local magnetic spots of reverse polarity _
superposed on the average steady dipole magnetism caused by the spins in condensed phase .
these pairs also have a tendency to align on north - south direction according to the potential form in eq.([vingcs ] ) .
these pairs have some similarity with the south atlantic magnetic anomaly ( sama ) observed on the core - mantle boundary ( korte 2010 ) .
the above interpretation that the spins in the expanded phase trigger the flip of the global magnetic polarity suggests that a possible _ correlation between the activity of the spins in expanded phase and the duration of the fixed polarity_. actually in our numerical calculations , the expanded phase temporally disappears in the superchron , i.e. the longest period of fixed polarity .
this correlation should be further studied .
anomalous acceleration in the movement of the local magnetic spots of reverse polarity may predict a polarity flip in the near future .
it is also interesting to notice that this set of equations eq.([eomgcs ] ) is a slight extension of the _ kuramoto model for synchronization_(kuramoto 2003 , acebron , bonilla , perez , ritort , & spigler 2005 ) by the change of potential : @xmath132 , and the reduction of the order of time derivative : @xmath133 .
the potential @xmath134 forces each angle @xmath56  to rotate with the fixed frequency @xmath135 despite the lack of the inertial term .
therefore the latter simplification @xmath136 is possible in kuramoto model .
kuramoto model is a set of long - range coupled oscillators with different frequencies .
if there were no interaction between the oscillators , then each phase of the oscillator behaves independently .
however if the long - range interaction is sufficiently strong , then the phases of oscillators _ synchronize _ with each other despite the difference in each frequencies @xmath137 .
kuramoto model very generally describes the synchronization processes of many type of oscillators and limit cycles including chemical and biological systems .
it is remarkable that lcs model can also describe the synchronization of element spins and yield _ quasi - periodic motion of the mean field _ , by simply changing the parameters .
the first term of rhs in eq.([vingcs ] ) gives independent non - linear oscillations of each spin angle . since the oscillation is nonlinear and the amplitudes are different from spin to spin , the phases of the spins behave very chaotically .
the second term of rhs in eq.([vingcs ] ) gives the cosine potential between all the spin pairs ( @xmath138 ) and thus naturally yields the tendency for all the pairs to synchronize with each other .
whether the whole synchronization actually takes place or not depends on the balance between the two terms . if the pairwise interaction @xmath139  is sufficiently large in lcs model , then the synchronization takes place as in the case of kuramoto model .
however , reflecting the strong non - linearity of the cosine potential , the synchronized oscillation becomes quasi - periodic in general .
this potential ability of lcs model to describe synchronization physics is especially important when we apply lcs to the solar magnetism in later sections .
we would like to argue , in this section , to what extent the lcs  model can be _ general _ in wider contexts .
the concept of the coupled dynamo elements itself , as we have discussed in section 2 , is general and is not in principle restricted to geomagnetism .
we expect the same for lcs model , which is the minimal realization of this concept .
the magnetic fields associated with celestial objects are quite common in the universe ; many planets , satellites , and the stars including our sun .
although the magnetic fields of these objects have variety that is quite different from the geomagnetism , all the magnetic fields are thought to be generated and supported by the robust dynamo mechanism .
however the lcs model is at most a phenomenological model to describe dynamo systems .
the lcs model is simply a macroscopic description and can not be directly deduced from the microscopic mhd equations eqs.([eqforv]-jb ) , which are especially important for quantitative argument in wider contexts .
therefore we _ need to supplement lcs model with the mhd equations_.    first we examine the _ energy flow _ for geomagnetism .
geomagnetic fields would be amplified by the inward - winding motion of the convection vortex , anti - cyclone .
this inward motion is caused by the coriolis force induced by the earth rotation although the rotation itself can not directly transfer its energy to geomagnetism .
on the other hand , the convective motion of the iron fluid is supported by the thermal flow from the central region of the earth .
this flow is mainly supported by the release of the gravitational potential by the core sinking .
thus the thermal flow accelerates the convective motion and supports the inward - winding force to amplify the magnetic fields through the coriolis force .
one essential fact for this thermal energy flow is that the main dissipation is due to the joule heat .
this is because the magnetic prandtl number @xmath140 , i.e. the ratio of the time scale of magnetic diffusion ( @xmath141 )  and that of convective diffusion ( @xmath142 ) , is actually very small : @xmath143 .
we now closely look at the basic equations eqs.([eqforv ] , [ eqforb ] ) in the above picture of energy flow . from eq.([eqforb ] ) , we first notice that the steady state for the magnetic field requires the _ energy balance on average _ : @xmath144 , which reduces to the relation which is independent of @xmath14 , @xmath145where @xmath146 is the fe - core radius and @xmath147 is the radius of the taylor cell .
the existence of this special scale is essential to our model to characterize the dynamo elements ; most of the relevant process of dynamo takes place on this scale .
this small scale compared with the core size is foreseen from the order in eq.([rm ] ) .
the magnetic field is amplified by the inward winding motion of the convection vortex through the term @xmath148 until this steady state eq.([v - steady ] ) is realized . on the other hand ,
the inward - winding motion of flow is generated by the term @xmath149 in eq.(eqforv ) , and this term should eventually _ balances _ with the backreaction from the generated magnetic pressure through the term @xmath150  in the same equation .
therefore in the steady state , we have the estimate of the typical strength of the magnetic field @xmath151 at the scale @xmath147 , the scale of the taylor cell , as    @xmath152    the observed magnetic field @xmath153  at the earth surface is related with this estimated @xmath151  through the magnetic flux conservation as , @xmath154where @xmath69 is the number of taylor cells in the core and @xmath155  is the earth radius .
we simply assumed that all the taylor cell is the same scale and the magnetic fields associated with tc are parallel to the rotation axes of the earth for simplicity . combining the above three equations
, we have @xmath156and the corresponding ( virtual ) magnetic moment @xmath157  becomes@xmath158    in the geomagnetic case , we have chosen @xmath87 and putting reasonable values , @xmath159 km , @xmath160 km , @xmath161emens / m , @xmath162kg / m@xmath163 , and assuming @xmath164 , we have @xmath165tesla from eq.([bout ] ) , and @xmath166m/@xmath167 from eq.([v - steady ] ) .
although the value of @xmath153is consistent with observations the value of @xmath168  is only @xmath169 of the estimated speed of the convection flow simply from the west - ward drift motion .
if this discrepancy is real , the time scale of the magnetic field variation or of the electric current change is about ten times larger than the convective flow speed .
this point should be further studied .
the estimate of the effective magnetic moment in eq.([d ] ) is general since it is deduced simply from the basic equations ( [ eqforv ] , eqforb ) based on the lcs model .
the expression eq.([d ] ) is consistent with the general claim in the literature ( christensen2010 ) , ( stevenson 2010 ) that the planetary magnetic fields are essentially determined from the factor @xmath170 .
we can step further and estimate the unknown factor @xmath171 in eq.([d ] ) . by using various data for celestial objects , but assuming @xmath172 is a constant
, we can draw a graph for @xmath157 against @xmath173 . ] .
this turns out to show a power law relation with the index about 3/4 and the angular moment @xmath174  for celestial objects @xmath175 where @xmath176meter@xmath177ampere)@xmath178 @xmath179 . ] : @xmath180 , though this is not what we needed .  simply motivated by this scaling relation
, we can draw a similar graph for @xmath181 against @xmath124  in fig.([fig6 - 1 ] ) .
this relation turns out to be almost a constant within a factor @xmath182 for the whole mass range of 8-digits as shown in fig.([fig6 - 1 ] ) .
is shown in the vertical axes against the mass @xmath124 in the horizontal axes for various celestial objects .
these quantities are normalized by the values of earth .
the data are from stevenson ( 2010 ) .
eulopa and callisto are excluded from the figure because their magnetic fields are probably induced from the external fields ( stevenson 2010 ) .
the value @xmath183 t is used for the solar dipole magnetic field .
the quantity @xmath184 turns out to be an almost constant within a factor @xmath182 for the mass range over 8 digits .
this fact suggests the scaling @xmath185 , which determins the number and the radius of the dynamo element in general .
, width=453 ]    the factor @xmath186 is important for this constancy yields the smallest variation of the data as far as we have checked .
when we choose @xmath124 for the extensive variable , the best index was @xmath187 which is @xmath188 within an error . ] .
we can deduce the scaling formula , from this phenomenological relation , as @xmath189 in other words , the scaling eq.([d ] ) holds if we choose eq.(taylorcellscaling ) .
this phenomenological scaling is helpful to reveal the number of dynamo elements .
for example , if we suppose the parameter @xmath190 is almost a constant for various objects , and @xmath87  for earth , we have @xmath191 for the sun .
however , we must keep in mind that the planets , satellites and stars are very different with each other and the structure of dynamo mechanism is different ( stanley 2010 ) despite the above generality .    keeping the above caution in mind , it is interesting to try the unified view of planetary magnetism within our model . in
the context of the lcs model , observed large dipole tilt more than @xmath192 and the quadrupole moment of magnetic fields in _ uranus and neptune _ are not exceptional .
the magnetic field described by lcs model naturally becomes irregular and asymmetric when occasional large excursion takes place .
we may observing these irregular and asymmetric fields in the present uranus and neptune .
moreover in the context of lcs model , almost axis symmetric magnetic field of _ saturn _ does not conflict with the anti - dynamo theorem , which inhibits dynamo function in  the axisymmetic steady configuration .
this is because individual element dynamo yields magnetic field independently from others and easily violates axis symmetry of the whole system .
the magnetic field in lcs model is rapidly changing and steady condition is also violated .
as we have seen in the previous section , a simple argument on the mhd equations based on the lcs model can yield a general scaling law for planets and satellites .
this was useful to estimate the scale of dynamo elements . in this section
, we would like to explore this generality of the lcs model by applying it to the solar magnetic dynamics .
lcs model does not specify the macro - spin or dynamo element .
the macro - spin can express any subsystem which has minimal dynamo property . typically the element is the inward - winding vorticity . in case of our sun
, the macro - spin may be the convecting vortex deep inside the solar surface .
it may happen that this region is turbulent and the hierarchy of vortices yield a network of many dynamo elements .
the very end of this hierarchy may appear on the entire surface of the sun , possibly as the supergranulations .
recently the pattern of the supergranulation is observed by using the local helioseismology ( nagashima 2010 ) . on top of these supergranulation , strong horizontal magnetic field
is observed everywhere on the surface ( ishikawa 2008 ) .
we speculate that these vortices associated with dynamo function may fill the entire convecting region inside the sun .
the solar magnetism changes its polarity quasi - regularly with the period about 22 years , which synchronizes with the solar activity of period 11 years .
the solar activity can be roughly measured simply by counting the sunspots ( sidc 2011 ) .
we show this daily count profile in fig .
( [ fig7 - 1 ] ) top , during the past two hundred years .
the power spectrum , shown in fig .
( [ fig7 - 1 ] ) bottom , reveals the clear periodicity of about 11 years , as well as the period of about a month .
the latter period is an artifact reflecting the systematic sunspot change on the solar surface which rotates with this period .
what is interesting for us is the scaling property with index @xmath193 , near to @xmath194 , on top of the two peaks associated with the above periods .
one - over - f noise or the pink noise , _
i.e. _ its power spectrum scales with index @xmath194 , has been observed everywhere in nature .
( for a comprehensive list of references on 1/f noise , see ( li 2011 ) ) .
similar 1/f noise in the monthly averages of the sunspot numbers has been reported in the high frequency ( @xmath195 ) regions ( franchiotti , sciutto , garcia , & hojvat 2004 ) .     and @xmath196 .
the former peak corresponds to the 11 year solar activity cycle and the latter peak fictitiously appears corresponds to a month the solar rotation period . , width=453 ]
it is remarkable that the lcs model shows quasi - periodic polarity flip , which resembles the solar magnetism , as a result of synchronization of many spins .
for the synchronization process to set in , we need large number of spins ( @xmath197 ) and low potential barrier ( @xmath198 ) .
then the spins in lcs model can naturally synchronize to yield quasi - periodic motion like in the kuramoto model .
for example in a numerical demonstrations , we set @xmath199 @xmath200 @xmath201 , and obtained the result as shown in fig.([fig7 - 2 ] ) .
the top figure shows the time evolution of the order parameter @xmath202 defined in eq.([m ] ) .
quasi - periodicity in its evolution is apparent .
the amplitude irregularly varies and even becomes almost flat in some periods .
the bottom figure shows its power spectrum . beside the apparent peak corresponding to the quasi - periodicity , one - over - f noise , actually the power @xmath203 ,
is observed in low frequency region .
however this result needs to be refined before serious discussions .
for example , the number of spins @xmath204 is apparently too small compared with the previous estimate @xmath205 in the last section .
if we suppose that the solar magnetic activity is linked with the number of sunspots , then the above results suggest that the lcs model can reproduce the observed solar activity especially _ the quasi - periodicity as well as 1/f like power law in the low frequency region in its power spectrum_.     in the numerical calculation with the parameter @xmath199
@xmath206 @xmath207 , and the initial condition is chosed so that all the spin angles are equally separated within the 70% of the whole angle @xmath129 .
each spin has the interaction with nearest @xmath182 spins .
the quai - periodic oscillation is apparent .
this oscillation appears as a result of synchronization of many spins .
( botom ) the power spectrum of this time series .
there is a peak at @xmath208 reflecting the above quasi - periodicity on top of the rough power law behavior with index @xmath209 in low frequency region .
these behaviors accord with the solar observation in the previous figure .
, width=453 ]
we have developed a long - range coupled macro - spin lcs model for geomagnetism , solar magnetism , and others .
we have seen this model is _ minimal and general_.    in _ section 2 _ , we developed the idea of the coupled dynamo elements for general dynamo mechanism in magnetohydrodynamical systems such as geomagnetism .
the element should have inward - winding vorticity to amplify the magnetic fields .
the we introduced a macro - spin representing such element . in _ section 3 _ ,
we have reviewed the short - range coupled spin ( scs )  model which successfully describes geomagnetism . in _ section 4 _
, we introduced the long - range coupled - spin ( lcs ) model which also successfully describes geomagnetism .
thus the spin models are shown to be useful to describe geomagnetism in the context of statistical fluctuations .    in _ section 5 _
, we studied the physical relevance of lcs model .
first of all in this model , _ two phases of spins can coexist _ : the expanded phase in which the spins can rapidly move all the directions almost freely and the condensed phase in which the spin directions are bounded into narrow region to form the dominant dipole moment which is almost steady .
the spins in expanded phase perpetually `` collide '' ( see the figure [ fig5 - 1 ] ) with the condensed phases and trigger its intermittent polarity reversal . on
the other hand lcs model shows quasi - periodicity as well as 1/f noise on top of it .
this is the _ synchronization _ of the constituent spins . in both cases ,
the system shows power law behavior reflecting the long range interactions .
the phase coexistence is essential for geomagnetism and the synchronization for solar magnetism .
these study taught us how the general models such as hamiltonian mean field ( hmf ) model and kuramoto model , if slightly modified , can describe the universality and variety of general dynamo mechanism .    in _ section 6 _ , we applied lcs model , supplemented with mhd equations , to other planets and satellites . finding a useful scaling low
, we could estimate the scale and the total number of the dynamo element , or macro - spin , in terms of the mass @xmath124 of the body . in _
section 7 _ , we applied lcs model , with different parameter values , to the solar magnetism and successfully describe the quasi - periodic oscillation of the solar magnetism on top of the 1/f noise property in low frequency regions .
thus we have actually demonstrated the universality and variety of lcs model .
although the lcs model can describe most of the geomagnetic observations , it can not be deduced from the basic mhd equations .
it is also possible to generalize lcs model to endow each spin with dynamo function .
this is realized for example by setting the homopolar generator ( faraday disk ) for each dynamo element , @xmath210 , and introduce appropriate interactions between them .
this is an extension of the original rikitake model ( rikitake 1957 ) .
if we introduce nearest neighbor interaction , we obtain the short - range coupled rikitake model ( scr ) , for the currents @xmath211  and rotation angles @xmath212,@xmath213@xmath214 we can obtain the occasional polarity flip and power law power spectrum also in this model .
however even @xmath215 case , the time evolution of scr model inherits chaotic spiky profile from the original rikitake model .
the long - range coupled version of rikitake model ( lcr )  is also possible,@xmath216@xmath217 where @xmath218 is the mean field . in this case
we obtain a regular and spiky oscillation due to the strong synchronization of the elements .
we have not yet performed a full parameter study .
further extension of such models would be possible and they will be worth extensive research .