the problem of heat conduction is a well - studied field .
more than two centuries ago , joseph fourier summarized the behavior of heat conduction by the law that bears his name .
this law describes phenomenologically that the heat current is proportional to the temperature gradient .
the detailed atomistic theories of heat conduction appeared only much later . for heat conduction in gas ,
the simple kinetic theory gives the result @xmath3 , where @xmath4 is specific heat , @xmath5 is average velocity and @xmath6 is mean - free path .
peierls theory of heat conduction in insulating solids @xcite is a classic on this subject .
these early theories deal with mostly the relevant three dimensions .
it turns out that low dimensional systems are more interesting and in some sense strange .
an analysis of a simple one - dimensional ( 1d ) harmonic oscillator model shows @xcite that there is no well - defined temperature gradient , the thermal conductivity diverges with system sizes as @xmath7 or @xmath1 , depending on the boundary conditions .
there is a general argument , for momentum conserving systems , that the thermal conduction in 1d is necessarily divergent @xcite .
there have been many analytical and numerical studies of 1d heat conduction ( see ref .
@xcite and @xcite for review ) .
we ll mention some of the most relevant papers to current work .
the work of lepri _ et al _ @xcite by mode - coupling theory and molecular dynamics suggests a divergent thermal conductivity exponent of @xmath8 , i.e. , @xmath9 for a 1d chain model with fermi - pasta - ulam ( fpu ) interactions .
mode - coupling theory is usually applied in the dynamics of liquids @xcite .
the first use of this theory in the context of heat conduction appears only recently , mostly due to lepri and his collaborators @xcite .
perverzev analyzed the same problem with the peierls theory for phonon gas and gave the same conclusion of @xmath8 exponent @xcite .
the result of @xmath8 is also supported by numerical simulation from several groups @xcite .
these results are supposed to be universal to some extent .
however , it is challenged by a different result of @xmath10 by narayan and ramaswamy @xcite , based on fluctuating hydrodynamics and renormalization group analysis . the numerical result for this 1/3 law is not convincing , as for the same model
 the hard - particle gas model  some obtained 1/3 @xcite , while others obtained different value 1/4 @xcite .
but for the fpu model , there is no good evidence for an exponent of 1/3 @xcite .
when momentum conservation is broken such as the one with on - site potential , the heat conduction can become normal again like the frenkel - kontorova model @xcite and the @xmath11 model @xcite .    in order to understand the underlying microscopic dynamical mechanism of the fourier law , a different class of models  billiard channels
have been introduced and studied in recent years @xcite .
various exponent values are found in such systems .
thus , it is believed that a universal constant does not exist at all . instead
, the divergent ( convergent ) exponent of the thermal conductivity is found related to the power of super ( sub ) diffusion @xcite .
besides the theoretical significance of heat conduction research in low dimensional systems , it is also of practical importance .
recent development of nanotechnology will enable us to manufacture devices with feature sizes at molecular level .
the understanding of heat conduction mechanism will allow us to control and manipulate heat current , and eventually to design novel thermal devices with certain function @xcite . to this end
, more realistic physical models are necessary . among many others , nanotubes and polymer chains are most promising .
there have been a number of numerical works to compute the thermal conductivity of the carbon nanotubes @xcite .
recent molecular dynamics ( md ) study of carbon nanotubes with realistic interaction potential suggested a divergent thermal conductivity for narrow diameter tubes @xcite .
the quantum effect of such systems is also very interesting @xcite .
we study the heat conduction of a 1d solid , as a classical system .
a brief version of this paper is reported in ref .
@xcite . in the rest of the paper ,
we introduce the quasi - one - dimensional chain model in sec .  [
sec : chain ] .
we discuss the basis of the mode - coupling theory , the projection method in sec .
[ sec : projection ] .
the mode - coupling approximations and their numerical and analytical solutions are discussed in sec .
[ sec : mode - coupling ] and [ sec : solution ] .
the basic output of the mode - coupling analysis is the dependence of damping of the modes with the wave - vector of the modes .
we find that the transverse modes are diffusive , with @xmath12 , while the longitudinal modes are super - diffusive , @xmath13 , where @xmath14 is the decay rate for mode with momentum or lattice wave number @xmath15 .
we discuss the relationship between the damping of the modes with the heat conductivity through green - kubo formula in sec .
[ sec : green - kubo ] .
our mode - coupling theory predicts that the heat conductance diverges with the @xmath10 exponent when the transverse motions are important , while @xmath8 is recovered if the transverse motions can be neglected . in sec .
[ sec : md ] , we present nonequilibrium molecular dynamics results ( with heat baths ) of heat conductance and compare with mode - coupling theory .
we conclude in the last section .
most of the previous studies considered only strictly 1d models , with the fermi - pasta - ulam ( fpu ) model as the most representative .
the strictly 1d models may not be applicable to real systems such as the nanotubes .
real systems of nanotubes or wires live in three - dimensional space . the added transverse motion and the flexibility of the tube at long length scales will certainly scatter phonons , and thus should have a profound effect on thermal transport .    while a direct simulation of a realistic system such as a polyethelene chain with empirical force fields such as that in refs .
@xcite or nanotube with tersoff potential @xcite is possible , we think it is useful to consider a simplified model which captures one of the important features of the real systems  transverse degrees of freedom . therefore we propose to study the following chain model in two dimensions @xcite : @xmath16 where the position vector @xmath17 and momentum vector @xmath18 are two - dimensional ; @xmath19 is lattice constant .
the minimum energy state is at @xmath20 for @xmath21 to @xmath22 .
if the system is restricted to @xmath23 ( corresponding to @xmath24 ) , it is essentially a 1d gas with harmonic interaction .
the coupling @xmath25 is the spring constant ; @xmath26 signifies bending or flexibility of the chain , while @xmath27 is the bond angle formed with two neighboring sites , @xmath28 , and unit vector @xmath29 , @xmath30 .
unlike the fpu model which does not have an energy scale , the second bond - angle bending term introduces an energy scale . in this work , we take mass @xmath31 , spring constant @xmath32 , and the boltzmann constant @xmath33 , thus the most important parameters are @xmath26 and temperature @xmath34 .
we follow the formulation of the projection method in ref .  @xcite .
let @xmath35 be a column vector of @xmath36 components of some arbitrary functions of dynamical variables @xmath37 .
each of the function @xmath38 can be complex .
later , we shall choose @xmath39 to be the canonical coordinates of the system .
we use @xmath40 to denote the hermitian conjugate of @xmath41 .
the equation of motion for @xmath41 is @xmath42 where @xmath43 is the liouville operator @xmath44 equation  ( [ eqmotion ] ) can be viewed as a partial differential equation with variables @xmath37 and time @xmath45 .
@xmath46 , i.e. , the quantity @xmath41 at time @xmath45 when the initial condition at @xmath47 is @xmath37 . quantity without a subscript @xmath45 will be understood to be evaluated at time @xmath47 , e.g. , @xmath48 .
a formal solution to eq .
( [ eqmotion ] ) is simply @xmath49 .
the projection operator on a column vector @xmath50 is defined by @xmath51 where @xmath52 and @xmath53 are @xmath54 matrices .
the angular brackets denote thermodynamical average in a canonical ensemble at temperature @xmath34 .
the comma separating the two terms is immaterial , but we use a notation of inner product .
one can verify that @xmath55 is indeed a projection operator , i.e. , @xmath56 .
if we apply the projection operator @xmath57 and @xmath58 to the equation of motion , we get two coupled equations . solving formally the second equation associated with @xmath59 , and substituting it back into the first equation
, we obtain an equation for the projected variable that formally resembles a langevin equation : @xmath60 where @xmath61 is the complex unit , and @xmath62 this set of equations is deterministic and formally exact .
the only assumption made is that equilibrium distribution can be realized .
what is more important is the correlation functions , which are physical observables .
we define the normalized correlation function ( correlation matrix ) as @xmath63 it is an identity matrix at @xmath47 and has the property that @xmath64 . from eq .
( [ langevin ] ) we have , @xmath65 this equation can be solved formally using proper initial condition with a fourier - laplace transform , @xmath66 = \int_0^\infty e^{-i z t } g(t)\ , dt,\ ] ] and similarly for @xmath67 .
the solution is @xmath66 = \bigl ( i(z - \omega ) + \gamma[z ] \bigr)^{-1}. \label{eqgw}\ ] ] to simplify notation , we have used parentheses ( e.g. , @xmath68 ) for functions in time domain , and square brackets with the same symbol ( e.g. , @xmath69 $ ] ) for their corresponding fourier - laplace transform in frequency domain .    the information about the system is in the memory kernel @xmath67 .
however , such a correlation function is difficult to calculate , since the evolution of the `` random force '' @xmath70 does not follow the dynamics of the original hamiltonian system .
for example , it is impossible to compute directly @xmath70 from molecular dynamics .
for this reason , a `` true force '' can be introduced which obeys the normal evolution , i.e. , @xmath71 the correlation function of the true force , @xmath72 and that of the random force are related in fourier space as @xcite @xmath73^{-1 } = \gamma_f[z]^{-1 } -                         \bigl ( i ( z - \omega ) \bigr)^{-1}. \label{eqf2r}\ ] ] this completes the formal theory of projection due originally to zwanzig @xcite and mori @xcite .
these results are formal and exact .
they give us relation between correlation of the `` force '' and correlation of dynamical variables .
they are the starting point for mode - coupling theory . in the next subsections
, we apply it to our chain model and introduce a series of approximations to solve it .
we now apply the projection method to our chain model .
we choose the normal modes as the basic quantities @xmath39 that we are going to project out .
there are several reasons for choosing the normal modes . to zeroth order approximation ,
each mode is nearly independent .
the slowest process corresponds to long wave - length modes .
the effect of short wave - length modes can be treated as stochastic noise ( the random force @xmath70 ) .
we choose @xmath41 to be the complete set of canonical momenta and coordinates : @xmath74 where @xmath75 we have defined the position vector @xmath76 , so that @xmath77 and @xmath78 are deviations from zero - temperature equilibrium position . because the fourier transform is a periodic function , the index @xmath79 is unique only modulo @xmath1 . as a result
, we can also let @xmath79 vary in the range @xmath80 .
we also note that @xmath81 .    with these definitions
, we can compute the matrix @xmath82 and expression for the true force @xmath83 in the general theory .
we find that for @xmath84 , the components are @xmath85 we have used equal - partition theorem for the average kinetic energy expression , and the last equation merely defines the effective frequencies @xmath86 for each mode .
note that due to translational invariance , correlation between different @xmath79 modes vanishes .
correlation between transverse and longitudinal modes also vanishes due to the reflection symmetry of @xmath87 for the chain .
thus equal - time correlation for @xmath41 is diagonal , @xmath88 we have defined @xmath89 as a @xmath90 diagonal matrix with elements @xmath91 ; @xmath92 is a @xmath93 identity matrix .
similarly , the correlation @xmath94 is found from @xmath95 the second equation is from a general virial theorem @xcite .
we have @xmath96 the expression for the true force @xmath97 is then @xmath98 note that only the momentum sector has a nonzero value , and @xmath99 is given approximately by eq .
( [ eqmotion1 ] ) and ( [ eqmotion2 ] ) below . with this special form of @xmath83 ,
the damping matrix @xmath100 is also diagonal and is nonzero only in the @xmath101 components . with these results , eq .
( [ eqgw ] ) becomes @xmath66   = \begin{pmatrix }                 iz d   &   - \tilde \omega^2 d \cr                 d           &   ( iz + \tilde \gamma[z ] ) d \cr                \end{pmatrix } , \label{general - gz}\ ] ] where @xmath102)^{-1}$ ] is a @xmath90 diagonal matrix , and @xmath103 $ ] is the fourier - laplace transform of the correlation @xmath104 , which is also diagonal ( we ll denote as @xmath105 ) .
in particular , we have the usual expression for the normalized coordinate correlation , @xmath106 & = & { i z + \gamma_k^\mu[z ] \over              \tilde \omega_k^\mu{}^2 - z^2 + iz \gamma_k^\mu[z]}. \label{eqvtog}\end{aligned}\ ] ] for simplicity , we drop the @xmath107 subscript for the coordinate correlation for the rest of the paper .    finally , the relation between the random force correlation and true force correlation , eq .  ( [ eqf2r ] ) , becomes , @xmath108 } = \frac{1}{\gamma_{f , k}^\mu[z ] }        - { i z \over \tilde \omega_k^\mu{}^2 - z^2 } .
\label{eqvftov}\ ] ] the force correlation @xmath109 can be computed from the correlation function @xmath110 .
it is convenient to separate the linear term in the force from the nonlinear contribution .
so we write @xmath111 where @xmath112 is effective angular frequency and @xmath113 is bare angular frequency of the mode .
we have dropped the indices @xmath79 and @xmath114 since these equations apply for any of the modes .
@xmath115 is at least quadratic in @xmath116 .
we note that the correlation of @xmath117 is linearly related to the coordinate correlation function @xmath118 $ ] , @xmath119 = { - iz(\tilde\omega^2-z^2 ) \over \tilde\omega^2 }   + { ( \tilde\omega^2 - z^2)^2 \over \tilde\omega^2 } g[z ] . \label{eqgtogamf}\ ] ] this is simply a consequence of eqs .
( [ eqvtog ] ) and ( [ eqvftov ] ) , but can also be derived directly from the definition .
the second derivative of @xmath120 in frequency domain is @xmath121 $ ] multiplied by @xmath122 .
using the fact that @xmath123 - \dot q(0 ) - iz q(0),\ ] ] we can also derive eq .
( [ eqgtogamf ] ) with the understanding that @xmath124 is an average over the initial conditions .
since @xmath118 $ ] is finite or at least should not diverge precisely at @xmath125 , this implies @xmath126 = 0 $ ] .    with the above results
, we can derive an expression of the true force correlation in terms of nonlinear force correlation , @xmath127 where @xmath128 .
the mixed term can be expressed in terms of @xmath118 $ ] by noting that @xmath129 .
the two mixed terms @xmath130 and @xmath131 are equal due to time - reversal symmetry . we find @xmath132 - iz.\ ] ] finally , we have @xmath119 =   { \delta \tilde\omega^2 \over \tilde\omega^2 }   \bigl((2z^2-\omega_0 ^ 2-\tilde\omega^2)g[z ]   + 2iz \bigr ) + \gamma_n[z].\ ] ] we can also express @xmath118 $ ] in terms of the nonlinear part of the force correlation , @xmath133 ( @xmath134 ) : @xmath135 = { iz(2\omega_0 ^ 2 - \tilde\omega^2-z^2 ) + \tilde\omega^2 \gamma_n[z ]     \over ( z^2-\omega_0 ^ 2)^2 } . \label{eqvntog}\ ] ] again , since @xmath118 $ ] can not diverge precisely at @xmath136 , this implies that @xmath137 $ ] must take a special form to cancel the apparent divergence .
thus , if we do not take care of these superficial divergences , we will not be able to make correct prediction for the correlation function .
we can also relate the original damping function @xmath100 to the nonlinear one , @xmath73 = { -i(\tilde \omega^2 - \omega_0 ^ 2)^2 z +                 \tilde\omega^2 ( \tilde\omega^2 - z^2 ) \gamma_n[z ]        \over \omega_0 ^ 4 - \tilde \omega^2 z^2 - iz\tilde \omega^2 \gamma_n[z]}. \label{eqvntov}\ ] ] this last equation is useful for approximating the damping function .
all of these relations are exact .
this ends our formal application of the projection method to the chain model .
to make some progress for analytic and numerical treatment , we have to make some approximations .
first , we ll consider small oscillations valid at relatively low temperatures .
an approximate hamiltonian for small oscillations near zero - temperature equilibrium position , keeping only leading cubic non - linearity in the hamiltonian , is then given by @xmath138 where @xmath139 are the ` bare ' dispersion relations , and @xmath140 the absence of @xmath141 term in hamiltonian ( [ h - perturb ] ) is due to the quadratic nature of the potential , while the absence of the terms of the form @xmath142 and @xmath143 is due to the reflective symmetry about @xmath144-axis of the hamiltonian .
we view @xmath79 and @xmath145 as independent component when taking the derivatives . a slightly modified hamilton s equation ( because of the use of complex numbers ) describes the dynamics : @xmath146 this gives the following equations of motion : @xmath147 where @xmath148 with these expressions , we are ready to compute the force correlation function in terms of dynamic variables @xmath149 .
we write @xmath150 where @xmath151 is the difference between bare dispersion relation and effective dispersion relation . the second term @xmath152 denotes the rest of the nonlinear force ( we take only the quadratic terms in @xmath116 ) .
due to translational and reflective symmetries , the correlation matrix formed by @xmath153 is diagonal without any approximation .
the time - displaced correlation for the diagonal terms defines true force correlation .
the nonlinear part of the contributions is @xmath154    in order to have a closed system of equations for the normalized correlation functions , we use the standard mean - field type approximation , @xmath155 .
owing to the @xmath156 correlation between different @xmath79 , the double summation can be reduced to a single one . we introduce @xmath157 and similarly @xmath158 associated with @xmath159 . in terms of @xmath160 , we obtain @xmath161 where @xmath162 equations ( [ eqmodekp - full - para ] ) and ( [ eqmodekp - full - perp ] ) together with the relations among @xmath163 , @xmath159 , and @xmath164 , eqs .
( [ eqvtog ] ) , ( [ eqvntov ] ) , and ( [ eqgammatonu ] ) , form a system of close equations , which can be solved in principle . however , because of the singular nature at @xmath136 in eq .
( [ eqvntog ] ) , any approximation to @xmath133 will destroy a subtle cancellation of the singularity , rendering the problem impossible to solve .
damping functions of the slowest modes , i.e. , real part of @xmath165 $ ] ( solid line ) and @xmath166 $ ] ( dashed line ) vs @xmath167 , computed from molecular dynamics for parameters of set e : @xmath32 , @xmath168 , @xmath169 , @xmath170 , @xmath31 .
the top figure is for size @xmath171 and bottom figure for @xmath172 . ]
the damping function @xmath173 $ ] is the central function that a successful theory needs to be able to calculate . in fig .
[ fig - mdgm ] we present examples of such functions determined from equilibrium molecular dynamics ( md ) simulation in a microcanonical ensemble with periodic boundary condition
. a more correct comparison of md with mode - coupling theory should use an ensemble of initial conditions distributed according to canonical weight .
this may be unimportant when @xmath1 is large .
we compute @xmath174 $ ] from the ratio of two correlation functions , @xmath175/g_{qp}[z ] = iz + \gamma[z]$ ] . for small systems
, there are a lot of peak structures associated with the low frequency modes ; the feature appears washed out for large systems .
to make the problem tractable , we make a bold approximation . instead of eq .
( [ eqvntov ] ) , we take @xmath176 \approx \nu_{n , k}^\mu[z ] .
\label{our - approximation}\ ] ] this is equivalent to say @xmath177 .
this appears justified for the longitudinal modes at sufficiently low temperatures but problematic for the transverse mode , as @xmath178 is linear in wave number @xmath179 but @xmath113 is quadratic in @xmath179 .
we can consider the limit of small @xmath179 .
in such limit , the difference between @xmath180 and @xmath100 is dropped .
lepri s treatment @xcite is similar to the above approximation .
we also note that in the work of scheipers and schirmacher for damping of anharmonic crystals @xcite , effective cubic coupling is used instead of the ` bare ' coupling , @xmath181 .
the standard operation of projecting the random force onto bilinear form @xmath182 to get the mode - coupling equations also agrees with the approximation , eq .
( [ our - approximation ] ) , but with somewhat different `` vertex functions '' replacing eq .
( [ eqkpara ] ) and ( [ eqkperp ] ) .    in order to obtain a numerical solution , besides the model parameters ( mass @xmath183 , lattice spacing @xmath19 , couplings @xmath25 and @xmath26 ) , we also need the effective dispersion relation .
we used md data for this purpose .
it turns out that a two - parameter fit of the form @xmath184 characterizes the effective dispersion relation very well , where @xmath185 is sound velocity at small @xmath79 and @xmath186 is maximum frequency at @xmath187 .
we used a fast fourier transform to solve the equations iteratively .
given an initial @xmath188 $ ] , correlation function in frequency domain is calculated by eq .
( [ eqvtog ] ) .
inverse transform gives us @xmath189 , which is used to calculate the discrete sum in eqs .
( [ eqmodekp - full - para ] ) and ( [ eqmodekp - full - perp ] ) to obtain @xmath190 .
we transform back to frequency domain for the next iteration
.     real part of @xmath165 $ ] ( solid line ) and @xmath166 $ ] ( dashed line ) vs @xmath167 from a full mode - coupling theory for parameters of set e ( same as that in fig.[fig - mdgm ] ) at @xmath171 and 64 .
the input effective dispersion relation parameters , @xmath191 , and @xmath192 , are given in table  [ tb-1 ] . ]
real part of @xmath165 $ ] ( solid line ) and @xmath166 $ ] ( dashed line ) vs @xmath167 from a full mode - coupling theory for parameters of set e , but at different temperature @xmath34 .
other parameters are the same as in fig .
[ fig - mkgm ] . ]
figure  [ fig - mkgm ] is our theoretical calculation with parameters exactly at set e. this figure should be compared with the md data in fig .
[ fig - mdgm ] .
we observe that the predicted results are too smooth without much structure .
the mode - coupling theory predicts a damping that is a factor of two larger for the longitudinal modes and much too large for the transverse modes . by adjusting the temperature to a lower value ( one half for @xmath171 and a quarter for @xmath172 )
, we obtain curves which closely resemble the md results . however , the mode - coupling theory still overestimates the value of transverse damping .
it is clear that @xmath193 is considerably smaller in actually system , but our mode - coupling theory can not produce a small enough @xmath193 .
the results presented in fig .  [ fig - mkgmh ] give us some hints that the mode - coupling equations are essentially correct , but the parameters of the model need to be adjusted .
in fact , simple perturbative expansion of the hamiltonian can not correctly predict @xmath194 . according to md result , @xmath195 , but the leading contribution from perturbation calculation should give the ` bare ' frequency , @xmath196 .
the adjusting of the parameters can be made more rigorous @xcite with an introduction of an effective hamiltonian , @xmath197 the form of the hamiltonian is dictated by symmetry .
translational invariance requires that @xmath198 .
the system is symmetric under reflection about @xmath144-axis or @xmath199-axis .
thus the hamiltonian should be invariant under the transformation @xmath200 or @xmath201 .
in addition , @xmath202 is symmetric under permutation of @xmath15 and @xmath179 , and @xmath203 is symmetric under permutation of all three indices .
we also have @xmath204 and similarly for @xmath203 .
although the original hamiltonian does not have @xmath205 term , such a term can present .
one of the reasons that such term can present is due to the non - analytic behavior of the potential ( because of the absolute value @xmath206 ) .
ldddd set e & bare & n=8 & n=64 & n=1024 + @xmath207 & 2 & 1.435 & 1.341 & 1.329 + @xmath208 & 0 & 0.621 & 0.669 & 0.674 + @xmath209 & 2 & 1.495 & 1.543 & 1.553 + @xmath210 & 1.095 & 1.087 & 1.182 & 1.194 + @xmath211 & 2 & 0.7113 & 0.6333 & 0.6337 + @xmath212 & 0 & -0.1060 & -0.0221 & -0.0076 + @xmath213 & -0.6 & -0.1099 & -0.1242 & -0.1293 + @xmath214 & 0 & 0.6362 & 0.8965 & 0.9312 + @xmath215 & 0 & -0.0424 & -0.1346 & -0.1469 +    we determine the parameters of the effective hamiltonian , @xmath216 , @xmath217 , @xmath218 , and @xmath205 , by fitting them to the observed time - independent correlation functions from md . in order to be able to carry out the calculation analytically
, we treat the interactions as perturbations , i.e. , @xmath218 and @xmath205 are assumed small .
all the equilibrium averages are approximated by the leading contribution from the perturbation , @xmath219 where @xmath220 is thermodynamical average with respect to the non - interacting harmonic oscillators ( product of gaussian integrals ) . under the above approximation , we find @xmath221 and the interaction parameters @xmath222 the actual form of the @xmath89 is fitted according to eq .
( [ omega_tilde_fit ] ) . for @xmath218 and @xmath205
we fit into a functional form @xmath223 where @xmath224 , @xmath225 , and @xmath226 .
the values of @xmath211 to @xmath215 , together with @xmath185 and @xmath227 are listed in table  [ tb-1 ] .
it is surprising that not only is the @xmath205 term present , but also is its magnitude comparable to @xmath218 .
the column marked ` bare ' are the parameters corresponding to the ` bare ' hamiltonian , eq .
( [ h - perturb ] ) .
the bare parameters are reached only at very low temperatures , such as @xmath228 .
since we have factored out the leading size dependence , we expect the parameters weakly depending on size @xmath1 .    before presenting our results with the effective parameters ,
it is worth pointing out that the standard procedure of deriving mode - coupling equations gives identical result as our approximation , eq .
( [ our - approximation ] ) .
the standard procedure is to project the random force , @xmath229 , onto a bilinear form of the basic variable @xmath230 . applying the general projection operator , eq .
( [ projector ] ) , to our case , this bilinear projector is @xmath231 where @xmath232 or @xmath233 is a pair of indices , e.g. , @xmath234 , and @xmath235 is @xmath236 . in evaluating the four - point correlation functions ,
we have neglected the perturbation term .
it turns out that the last term projects out only the @xmath237 component .
thus , the specific values of @xmath238 are not needed .    after applying projector @xmath239 onto @xmath240
, we neglect the difference between normal dynamics @xmath241 and anomalous dynamics for @xmath242 , @xmath243 .
we have @xmath244 the mode - coupling equations are then , @xmath245 where @xmath246 note that the vertex coupling @xmath247 is proportional to the three - point correlation independent of the specific form of random force .
this version of mode - coupling equation need not have the limitation of small oscillations , and can be applied to high temperatures as well .
real part of @xmath165 $ ] ( solid line ) and @xmath166 $ ] ( dashed line ) vs @xmath167 from a full mode - coupling theory for set e , using effective parameters given in table  [ tb-1 ] . ]
figure  [ fig - mkvgm ] is a calculation of the real part of @xmath248 $ ] for @xmath249 . comparing with fig .
[ fig - mkgm ] and [ fig - mkgmh ] , we see that using effective parameters brings into much better agreement with the md data .
in fact , comparing to fig .
[ fig - mdgm ] , most of the features are reproduced , such as the locations and heights of the peaks . the most important improvement is the ratio of parallel to perpendicular damping .    the normalized correlation functions @xmath250 ( upper part ) and @xmath251 ( lower part ) on a @xmath252 system for data set e. the solid curves are from full mode - coupling theory , while the dotted dashes are from equilibrium molecular dynamics . ]    in figure  [ fig - mkgt ] , we compare the normalized dynamic correlation functions for size @xmath252 and the slowest mode @xmath249 .
the frequencies of oscillations are slightly different in md and mode - coupling calculation , thus there are phase shifts at long times .
the longitudinal damping agrees with each other extremely well .
a fit of the logarithm of amplitudes ( maxima and minima ) versus time for the md data gives the decay rate @xmath253 , @xmath254 , while mode - coupling predicts @xmath255 and @xmath256 , respectively .
mode - coupling theory with effective parameters still overestimates the transverse damping by a factor of 2 to 3 .
the decay rate @xmath257 for each mode @xmath79 for system size @xmath258 .
circles ( @xmath259 ) and triangles ( @xmath260 ) are from molecular dynamics by fitting the correlation function , while the continuous curves are from a full mode - coupling theory for set e , using effective parameters given in table  [ tb-1 ] . ]    in fig .
[ fig - gamma - compare ] , we show the damping of the mode in terms of the decay constant @xmath257 , which is defined by the exponential decay in @xmath261 .
the symbols are obtained from md simulation , while the curves are from the full mode - coupling theory with effective parameters . for md data ,
@xmath257 is obtained by a least - square fit to the amplitudes . for mode - coupling theory
, we use an approximation , @xmath262 $ ] .
the bar means we average over a window centered around @xmath263 with width about @xmath257 .
we observe that excellent agreement between md data and mode - coupling theory is obtained for the parallel component .
the mode - coupling theory overestimates the transverse mode damping by some constant factor . in any case
, the slopes of the curves agree very well with the expected results ( to be discussed in the next subsection ) , @xmath264 and @xmath265 .
to simplify the equations further , we consider the large size limit .
since the large size asymptotic behavior is the most interesting , such limiting results are justified .
we thus consider the limit of @xmath266 and keep only leading contributions in lattice momentum @xmath267 . in this limit
the three kernel functions @xmath268 become constants . in addition , we define a @xmath79-independent @xmath269 function by taking the limit @xmath270 .
the leading @xmath79-dependence in @xmath271 is @xmath272 .
when this term is factored out , we expect that @xmath273 becomes independent of @xmath79 in the limit of large @xmath1 and small @xmath15 . with these simplification and approximation , and converting the discrete summation to integral , we have in the limit of large sizes : @xmath274 where @xmath275 if we use the perturbation expansion hamiltonian , or @xmath276 if we use the effective hamiltonian .
we also linearize the dispersion relation so that @xmath277 , @xmath278 or @xmath208 are effective sound velocities for the longitudinal and transverse modes .    we shall refer to eqs .
( [ eqsimple - mode - para ] ) and ( [ eqsimple - mode - perp ] ) as simple mode - coupling theory and the finite @xmath1 version as the full mode - coupling theory .
assuming that the simplified mode - coupling equations represent the essence of the physics regarding the damping and time - dependent correlation functions , we now consider analytic solution to eqs .
( [ eqsimple - mode - para ] ) and ( [ eqsimple - mode - perp ] ) .
first we notice some constraints on the functions .
because @xmath279 is a real function , we must have @xmath188^ * = \nu[-z]$ ] . since @xmath280 , we must have @xmath281 , for @xmath282 or @xmath283 ( assuming @xmath284 ) .
this implies that the fourier - laplace transform must be integrable : @xmath285 \,dz = \frac{k}{a}.\ ] ] next , the large @xmath167 behavior can be obtained by integrating by part few times , using the boundary conditions @xmath286 , @xmath287 , @xmath288 , and @xmath289 . for simplicity , we have omitted the lattice momentum index @xmath179 and mode index @xmath290 or @xmath283 , since it is true for any @xmath291 . with these results ,
we find @xmath292 & = & \int_0^\infty \!\!\nu(t)\ , e^{-izt } dt   \nonumber\\   & = & \frac{k}{2\pi } \int_{-\pi / a}^{\pi / a}\!\!\!\!\!dq \int_0^\infty\!\ !        g_q^1(t)\ , g_q^2(t)\ , e^{-izt } dt \nonumber \\ & = & \frac{k}{iaz } + \frac{k}{2\pi iz^3 } \int_{-\pi / a}^{\pi / a}\!\!\!\!\!\!dq(\tilde \omega_1 ^ 2 + \tilde \omega_2 ^
2 ) + o(\frac{1}{z^{5}}),\end{aligned}\ ] ] where for @xmath293 , @xmath294 and @xmath295 are all @xmath296 , while for @xmath297 , it is @xmath298 and @xmath296 .
we can establish that @xmath299 \propto ( iz)^{-1/2 } , \quad \nu^\perp[z ] \propto { \rm const},\ ] ] in the limit of small @xmath193 .
we proceed as follows : first we assume that @xmath300 = \nu_0^\perp$ ] or @xmath301 .
we then derive an expression for @xmath302 and @xmath303 , and show that indeed @xmath303 approaches a @xmath156-function in a proper limit .
when @xmath300 $ ] is a constant , the inverse fourier transform for the correlation function @xmath118 $ ] can be performed exactly , to give the transverse correlation function as @xmath304 where @xmath305 , and @xmath306 .
we have linearized the dispersion relation , @xmath307 .
we only look at the dominant term in time dependence of @xmath302 . for @xmath308 , since large @xmath179-mode decays rapidly , we ignore the term with amplitude proportional to @xmath309 , and also drop the oscillatory term and approximate , @xmath310 . after integrating over @xmath179 ( we also extend the limit as from @xmath311 to @xmath312 ) , we obtain the leading @xmath45 dependence , as @xmath313 if the oscillatory terms are kept , they only contribute an exponentially small term , @xmath314 ; and the @xmath309 terms give a contribution that decays much faster , as @xmath315 . the contribution from @xmath316 term can be neglected , because it decays as @xmath317 ( due to eq .
( [ gpara - of - t ] ) ) .
but this term does provide a crossover from @xmath318 for intermediate @xmath167 to the asymptotic @xmath319 for very small @xmath167 .
the fourier - laplace transform gives @xmath299
\approx \frac{k^\parallel}{4 }    \sqrt { \frac{1}{i\nu_0^\perp z } } = b ( iz)^{-1/2}. \label{nu - para - z}\ ] ]    we now get an expression for the longitudinal correlation function . formally , it is given by the inverse transform @xmath320 \over   z^2 - c^\parallel{}^2 q^2 - i z \gamma_q^\parallel[z ] } e^{izt } dz,\ ] ] where @xmath321 = q^2 \nu^\parallel[z ] = bq^2/\sqrt{iz}$ ] , @xmath322 .
the dominant contribution is from @xmath167 when the denominator is close to zero .
the integral can be approximately estimated by the residue theorem . by location
the poles @xmath323 we obtain the dispersion relation for @xmath324 . in the limit of small @xmath179
, we find @xmath325 where @xmath326 .
therefore , @xmath327 similarly with the dispersion relation for the transverse mode , @xmath328 , we can compute @xmath329 we have neglected the @xmath330 term . since @xmath331 is smaller than @xmath332 , we drop the second term .
we symmetrically extend the function . note that @xmath303 can be casted into the functional form @xmath333 which has the property that @xmath334 thus @xmath335 behaves as a @xmath156-function as @xmath336 :
@xmath337 this allows us to identify the constant @xmath338 .
the condition for @xmath336 is the same as @xmath339 or @xmath340 . since @xmath341 ( c.f .
eq.([kratio ] ) ) , small @xmath193 also implies small @xmath208 .
although the above derivation suggests that the result is valid only for the special limit .
but in fact , this is also the asymptotic result in the limit @xmath342 . to show this
, we note that the asymptotic behavior is picked up by the scaling @xmath343 \approx \nu[z]$ ] .
if this is true , then @xmath188 \propto z^{-\delta}$ ] .
this scaling limit coincides with the limit of @xmath340 .
for the simple mode - coupling theory that we have already taken the large @xmath1 limit , we can also apply the fast fourier transform method .
however , we find that a direct numerical integration in frequency space is much more accurate .
for example , @xmath300 $ ] is expressed as @xmath344 & = & { k^\perp \over ( 2\pi)^3 }   \int_{-\infty}^{+\infty}\!\!\!\!\!\!\ ! d\omega   \int_{-\infty}^{+\infty}\!\!\!\!\!\!\ ! d\omega ' \,2\,f_\perp(\omega , \omega ' ) \nonumber \\ & & \left ( { \rm p } { 1 \over i(z-\omega-\omega ' ) }   + \pi \delta(z-\omega-\omega ' ) \right ) ,   \label{principle - value - int}\end{aligned}\ ] ] where p stands for cauchy principal value and @xmath345 q^2 \over             \omega^2 - ( c^{\parallel,\perp})^2 q^2 - i\omega \nu^{\parallel,\perp}[\omega ] q^2 } .\end{aligned}\ ] ] when the linear approximation is made to the dispersion relations , the integral over @xmath179 can be performed analytically .
we only need to do a two - dimensional integral over @xmath346 and @xmath347 .
the principal value integral is taken care by locating exactly the singularity .
since the integration routine needs a smooth function @xmath188 $ ] , we fit the results by a pad approximation ( in variable @xmath348 or @xmath349 ) .
the procedure is iterated several times for convergence .
this is programmed in _ mathematica_.
it turns out essential to have more than double precision accuracy in the integration routine in order for the singular integrals properly converged .
some results of this calculation are already presented in ref .
.     real part of @xmath350 $ ] ( top ) and @xmath300 $ ] ( bottom ) for data set e. the smooth curves are from simple mode - coupling theory , the ( red ) curves with spikes are computed from full theory at @xmath258 , while the straight , dash lines are analytic asymptotic results . ]
real part of @xmath350 $ ] ( solid line ) and @xmath300 $ ] ( dashed line ) for data set j ( @xmath351 , @xmath352 , and @xmath353 ) .
the insert shows the effective exponent @xmath354 . ]    in fig .
[ fig - nu ] , we compare three levels of approximations of @xmath355 $ ] : the simple mode - coupling theory , the full theory computed on @xmath258 by @xmath355 = \gamma_1^\mu[z]/(2\pi/(na))^2 $ ] , and the asymptotic result of @xmath300 = \frac{1}{2 } k^\perp \bigl(1/|c^\parallel - c^\perp| + 1/|c^\parallel + c^\perp|\bigr)$ ] , and eq .  ( [ nu - para - z ] ) .
noticeable differences are seen between full theory and simplified .
this is partly due to the fact that we are comparing a finite @xmath1 with an infinite @xmath1 result .
we also note that the asymptotic slope of @xmath356 for @xmath357 $ ] is approached rather slowly .
figure  [ fig - nu - crossover ] shows a stronger crossover effect from slope of @xmath10 ( corresponding to @xmath9 , to be discussed later ) to that of the true asymptotic law of 1/2 .
this set of curves corresponds to data set j.       real part ( top ) and imaginary part ( bottom ) of @xmath358 vs @xmath167 for @xmath359 and @xmath360 ( solid line ) , 1/8 ( open circles ) , and 8 ( crosses ) . ]
we now look at some of the scaling properties that are implied from eqs .
( [ eqsimple - mode - para ] ) and ( [ eqsimple - mode - perp ] ) .
first , we consider the simple case of @xmath361 , @xmath284 , and @xmath362 . in this case
the two equations degenerate into one equation given by lepri @xcite . by measuring frequency in terms of @xmath185 (
i.e. , consider variable @xmath363 ) we can scale away @xmath185 , thus the following equation is an exact scaling @xmath364 where @xmath269 is a function of @xmath167 with parameters @xmath365 , @xmath185 , and @xmath19 that we have written out explicitly .
this equation tells us that out of the three parameters that characterize the solution , we can pick two as independent .
the `` shapes '' of the solutions are all the same for @xmath366 . without loss of generality , we can fix @xmath367 .
next , we consider a general scaling solution of the form @xmath368 substituting this scaling form into eq .
( [ principle - value - int ] ) , by requiring that the result must be consistent , we obtain the following conditions @xmath369 where @xmath370 is scaling exponent associated with integration variable @xmath179 , @xmath371 . a unique solution to the set of linear equations
is obtained , @xmath372 , @xmath373 , @xmath374 , and @xmath375 . since eq .
( [ lepri - scaling ] ) is ( approximately ) valid for any @xmath376 , we can choose @xmath377 .
this scaling solution implies that @xmath378 power - law behavior @xmath318 is obtained in the limit of small @xmath167 , relatively large @xmath185 , and small @xmath19 , if @xmath379 is finite .
the crossover to other behavior occurs at large @xmath380 and @xmath381 .
a simple dimension analysis also leads to the @xmath318 factor .
let the dimension of length and time be @xmath382 and @xmath34 , respectively .
then the dimensions of relevant quantities are @xmath383 = t^{-1}$ ] , @xmath384 = lt^{-1}$ ] , @xmath385 = l$ ] , @xmath386 = l^2 t^{-1}$ ] , and @xmath387 = l^3 t^{-2}$ ] . from the five quantities ,
we can construct three dimensionless variables @xmath388 if there is any relation between the @xmath389 s , it must be in the form @xmath390 ( buckingham pi theorem @xcite ) , or @xmath391 = z^{-1/3 } k^{2/3 } f\left ( { z a\over c } , { k \over a c^2 } \right),\ ] ] where @xmath392 is an arbitrary , dimensionless function .
this result , of course , is consistent with eq .
( [ lepri - scaling ] ) .
this also suggests that the scaling is not approximate , but an exact result .
[ fig - lepri - scaling ] is a test of the above scaling by numerical solutions . for small @xmath167 ,
the power - law @xmath318 is verified to high accuracy . since this scaling is exact , the deviations are due purely to numerical errors in solving the equations .
the case of the coupled longitudinal and transverse equations is somewhat difficult to analyze .
we can require a very general scaling of the form @xmath393 and similarly for the perpendicular component .
if we require for consistency of scaling for both longitudinal and transverse components , we find that the only scaling solution is the symmetric solution , i.e. , the scaling discussed above with identical longitudinal and transverse scaling exponents .
if we abandon the exact scaling for the transverse component and look only at longitudinal component , we can require that @xmath394 as before we can require that there is no scaling for the coupling constant @xmath395 without loss of generality ( i.e. , @xmath396 ) , the above equations imply a relation @xmath397 in particular , if @xmath398 , we must have @xmath399 and @xmath400 .
this is consistent with the fact that the scaling region of @xmath401 is obtained for small @xmath208 .
a complete and clear scaling picture is still lacking . from the numerical solution
, we observe that besides well characterized scaling regimes , there are also plenty of intermediate regions and crossovers . for very large @xmath208 and small @xmath193 , or small @xmath208 and large @xmath193 ,
the behavior of the solutions are difficult to characterize .
in this section , we make the connection of the damping of modes with thermal conduction .
the starting point is the green - kubo formula for transport coefficient : @xmath402 for the special case of heat conduction in one - dimensional chain , it is better to call @xmath403 heat conductance rather than heat conductivity . in analogous to electric circuit , the heat conductance relates the temperature gradient to energy current ( fourier law ) , @xmath404 the quantity @xmath405 is related to energy current by @xmath406 .
the central quantity that we compute in equilibrium and nonequilibrium heat conduction is the ( total ) heat current @xmath407 .
since the heat current is a macroscopic concept determined by the conservation of ( internal ) energy , a microscopic version of it is not unique .
the current expression must satisfy the energy continuity equation in the long - wave limit .
we derive an expression starting from @xmath408 , where the local energy per particle is @xmath409 . by regrouping some of the terms using translational invariance
, we arrive at the heat current per particle : @xmath410 where @xmath411 , @xmath412 , @xmath413 is a unit vector in direction of @xmath414 . the total heat current in @xmath199 direction , @xmath415 is the quantity appearing in the green - kubo formula .
it is equal to the macroscopic heat current density ( energy per unit area per unit time ) integrated over a volume .    for theoretical analysis
, we need to express @xmath407 in terms of the dynamical variables @xmath416 and @xmath230 . a general expression will be rather complicated .
again as in mode - coupling approach , we ll consider small oscillation expansions , and consider only the leading contributions .
neglecting the nonlinear contribution of @xmath417 , we obtain @xmath418 where @xmath419 is the bare dispersion relation given by eqs .
( [ bare-1 ] ) and ( [ bare-2 ] ) .
more general expression in a quantum - mechanical framework including the cubic terms is given in ref .
@xcite .
an approximation for the correlation function of the current can be obtained , again using a dynamic mean - field approximation , as @xmath420    the above expression can be further simplified by the approximation @xmath421 .
we find @xmath422 we have dropped the contribution from the perpendicular component , because it decays much faster ( due to an extra @xmath423 factor in the integrand ) .     the green - kubo integrand from md ( circles ) and mode - coupling theory ( curve ) for set e with size @xmath258 .
the straight has a slope @xmath424 .
the insert shows the same data on a linear scale . ]
we compare the mode - coupling result of green - kubo integrand with that of md in fig .
[ fig - jt - compare ] .
we have used eq .
( [ gk - mkp ] ) for the mode - coupling calculation .
note that the correlation functions @xmath425 ,
@xmath426 , and @xmath427 are simply related in frequency domain through eq .
( [ general - gz ] ) .
the agreement is reasonable with the largest deviation about a factor of 2 .
the asymptotic behavior of @xmath317 is consistent with both sets of results .
the finite lattice conductance defined by green - kubo formula at parameter set e. solid circles are from full mode - coupling theory ; the triangles are from molecular dynamics with periodic boundary condition ; the straight lines have slope of 1/2 . ]
the heat conductance on a finite lattice is sensitive to boundary conditions .
clearly , in the mode - coupling formulation , we have used periodic boundary condition .
if we integrate over time from 0 to @xmath312 first on the second line of eq .
( [ mk - t - dependence ] ) , we obtain @xmath428 . on a finite lattice
we have a lower momentum cut - off , @xmath429 .
the size dependence of the conductance on finite lattice is then @xmath430 since @xmath401 , @xmath431 .
this result is in agreement with that of deutsch and narayan for a model ( the random collision model ) where transverse motion is taken into account only stochastically @xcite .    in fig .
[ fig - kappa ] , the mode - coupling result of @xmath432 is compared with equilibrium molecular dynamics result with periodic boundary conditions .
excellent 1/2 power is observed for the mode - coupling result .
we found it is rather difficult to get converged value from molecular dynamics .
but the results are consistent with a @xmath433 law .
the relation of the mode - coupling theory with the results of nonequilibrium situation of low and high temperature heat - baths at the ends is not clear cut .
the standard assumption is that the correlation should be cut off by a time scale of order @xmath1 due to the interaction with the heat baths .
finite size result is obtained by integrating the power - law decay to a time of @xmath434 .
this gives us the asymptotic behavior for the conductance at large @xmath1 as @xmath435 since we find @xmath401 for small @xmath193 , the thermal conduction diverges as @xmath0 . when @xmath316 is sufficiently large , we observe @xmath2
we note that the interaction potential is not smooth at the point when two particles overlap .
this can cause numerical instability .
thus , we replaced the original potential with a modified one , @xmath436 which smooths out the discontinuous derivatives at @xmath437 with a small correction of order @xmath438 .
in addition to replacing the spring potential by @xmath439 , we also replace the cosine term by @xmath440 so that the value is well - defined for all positions @xmath441 .    in actual simulation
, we have used very small @xmath442 to @xmath443 so that its effect should be comparable to error caused by finite time - step @xmath444 .
we also used large @xmath445 as a way of simulating a slightly different model to study the robustness of the results .
we solve the system of equations @xcite @xmath446 where @xmath447 is the total force acting on the @xmath232-th particle , @xmath448 and @xmath449 obey the equation @xmath450 where @xmath451 and @xmath452 are the temperatures of the two heat baths at the ends .
the summation is over the particles belonging to the heat bath .
we have used 4 particles for the nos - hoover heat baths , with extra first two and last two particles at fixed positions .
the coupling parameter @xmath453 is taking to be 1 .
for the central part , we used a second - order symplectic algorithm ( or equivalently , the velocity verlet algorithm ) , while for the heat - bath , we used a simple difference scheme accurate to second order in @xmath444 for @xmath454 .
although the total energy ( hamiltonian ) is no longer a conserved quantity when the heat - bath is introduced , the above equation still has a conserved quantity of similar character : @xmath455 this quantity can be used to monitor the stability of the algorithm .    since we run for very long time steps of @xmath456 to @xmath457 , it is essential that the algorithm is stable over extended period of time .
even with a symplectic algorithm , stability is not guaranteed by merely taking small @xmath444 @xmath458 .
thus , we adjusted the conserved quantity , eq .
( [ conserved ] ) , to its starting value after certain number of steps .
temperature of @xmath232-th particle , @xmath459 vs scaled position , @xmath460 , for data set e with @xmath461 ( plus ) , 256 ( dash line ) , and 1024 ( solid line ) . ]    in fig .
[ fig - temperature - scaling ] , we observe that good temperature profiles are established .
due to relatively large temperature difference between low and high temperature heat baths , and perhaps also due to the nature of heat baths , the profile is not linear .
however , the scaling with @xmath1 is approximately obeyed .
@xmath462 vs @xmath1 for different parameters @xmath463 . set c : @xmath464 ,
the chain is compressed to have an average distance between particle @xmath465 ; set g : ( 10 , 0.2 , 0.4 , 0 ) ; set i : ( 0.1 , 0.3 , 0.5 , 0.2 ) ; set k : ( 0.5 , 1.2 , 2 , 0.4 ) ; set l : ( 25 , 1 , 1.5 , 0.2 ) .
all of the sets have @xmath32 , mass @xmath31 , and lattice constant @xmath169 .
the number indicates the slope of least - squares fit . ]
in addition to the thermal conductance results reported in ref .
@xcite , fig .
[ fig - jl - extra ] gives additional data for variety of parameters .
one of the aim of these additional runs is to check the robustness of the @xmath10 power law for thermal conductance .
it is found in ref  @xcite for set e exponent @xmath466 .
this is an excellent confirmation of the mode - coupling theory .
we note that at parameter set @xmath365 , the @xmath10 law is not destroyed by introducing large @xmath445 .
thus we believe that the @xmath10 law is not specific to the potential used in ref .  @xcite .
in fact , the fermi - pasta - ulam model with transverse motion was studied @xcite and result is also consistent with the @xmath10 power law .    at low temperatures for set i and c
, we appear to observe exponent of 0.4 .
what is particularly interesting for set c is that the chain is compressed to an average distance of 1.5 from the equilibrium distance of 2 .
the behavior is not changed much comparing to uncompressed chain . set g and
l represent very large angular coupling @xmath26 .
we note that @xmath467 corresponds to a situation that equilibrium can not be established .
as @xmath26 becomes larger , the equilibration becomes increasingly difficult . in ref .
@xcite we reported a logarithmic behavior for @xmath468 ( set b ) .
this logarithmic behavior does not maintain when @xmath26 is increased further .
the exponent @xmath469 for set l and g is closed to 1/2 .
this may be related the results of disordered harmonic chain .
lepri _ et al _ first introduced mode - coupling theory into the problem of heat conduction in low dimensional systems to interpret the md results qualitatively .
we have developed the theory further by introducing a ` full theory ' and effective parameters for theory .
the full theory is solved numerically . at this level of approximation
, we find that mode - coupling theory gives excellent prediction on a quantitative level .
the longitudinal damping ( decay rate ) @xmath259 agrees with md data with few per cents of deviation .
the mode - coupling theory somewhat overestimated the transverse damping @xmath260 by a factor of two to three .
the full theory assumes that the damping function ( memory kernel ) @xmath470 $ ] is an explicit function of two variables , @xmath79 and @xmath167 . in a simple mode - coupling theory ,
we assume @xmath470 \approx ( 2\pi k/(na))^2 \nu[z]$ ] , valid for small @xmath471 .
the simple theory is computationally efficient , and contains essential asymptotic features , i.e. , @xmath264 and @xmath472 . through detailed comparison between md and mode - coupling theory on memory kernel , decay of the modes , the time - dependent correlation functions , green - kubo integrand , and finally the finite - size conductance
, we have strong evidence that the mode - coupling theory is correct , and it captures the essential features of the system .
we have used parameter set e as the main data set for comparison .
the reason for choosing this particular set of parameters is that it gives the cleanest power - law behavior .
other parameters will be either qualitatively similar , or crossovers will be observed .
in fact mode - coupling theory naturally predicts crossover , due to the presence of @xmath316 term . when this term is dominant , or the transverse effect
can be neglected , we should see behavior of @xmath473 .
thus we feel that the controversy of the result of lepri _ et al _ @xcite and that of narayan and ramaswamy ( @xmath474 ) @xcite can be reconciled within the current mode - coupling theory .
w. is supported by an academic research grant of national university of singapore and by singapore - mit alliance .
b. l. is supported by academic research grant of national university of singapore .
b. li , l. wang , and b. hu , phys .
lett . * 88 * , 223901 ( 2002 ) ; b. li , g. casati , and j. wang , phys .
e * 67 * , 021204 ( 2003 ) ; b. li , g. casati , j. wang , and t. prosen , cond - mat/0307692 ; d. alonso , r. artuso , g. casati , and i. guarneri phys .
lett . * 82 * , 1859 ( 1999 ) ; d. alonso , a. ruiz , and i. de .
vega , phys .
e * 66 * , 066131 ( 2002 ) .
l. i. manevitch and a. v. savin , phys .
e * 55 * , 4713 ( 1997 ) ; a. v. savin and l. i. manevitch , _ ibid .
_ * 61 * , 7065 ( 2000 ) ; a. v. savin , l. i. manevitch , p. l. christiansen , and a. v. zolotaryuk , physics - uspekhi * 42 * , 245 ( 1999 ) .