heat conductivity of 1d lattices is well known classical problem related to microscopic foundation of fourier law .
the problem started from the famous work of fermi , pasta and ulam @xcite , where abnormal process of heat transfer has been detected in the first time .
non - integrability of the system is necessary condition for normal heat conductivity . as it was demonstrated recently for the fpu lattice @xcite , disordered harmonic chain @xcite , diatomic 1d gas of colliding particles @xcite and the diatomic toda lattice @xcite , the non - integrability is not sufficient in order to get normal heat conductivity .
it leads to linear distribution of the temperature along the chain , but the value of heat flux is proportional to @xmath1 , where @xmath2 and @xmath3 is the number of particles in the chain .
thus , the coefficient of the heat conductivity diverges in thermodynamical limit @xmath4 .
analytical estimations @xcite have demonstrated that any chain possessing acoustic phonon branch should have infinite heat conductivity in the limit of low temperatures . from the other side , there are some artificial systems with on - site potential having normal heat conductivity @xcite
. heat conductivity of frenkel - kontorova chain in the first time was considered in paper @xcite .
finite heat conductivity for certain set of parameters was obtained for frenkel - kontorova chain @xcite , for the chain with sinh - gordon on - site potential @xcite and for the chain with @xmath5 on - site potential @xcite .
these models are not invariant with respect to translation and the momentum is not conserved .
it was supposed that the on - site potential is extremely significant for normal heat conduction @xcite and that the anharmonicity of the on - site potential is sufficient to ensure the validity of fourier law @xcite .
detailed review of the problem is presented in recent paper @xcite .
peculiarities of the heat conduction of the frenkel - kontorova model for complete set of parameters and temperature of the system are not known .
the chains with zero average pressure were demonstrated to have normal heat conductivity @xcite . in papers @xcite
the transition from abnormal to normal heat conductivity has been detected at certain temperature .
paper @xcite contains detailed review of the heat conductivity peculiarities in 1d molecular systems .
there are no detailed investigations similar to mentioned above and concerning the properties of the chains with on - site potential in the whole temperature range .
as it was mentioned above , there exists certain incompleteness of the knowledge concerning even the most popular and paradigmatic discrete frenkel - kontorova chain .
this lattice is of special interest as its continuous counterpart is famous sin - gordon system ( having , of course , divergent heat conductivity ) .
the transition between two regimes with the growth of temperature is expected for discrete system ; however it may be dependent also on other parameters of the lattice    the question of special interest is also the mechanism of heat flow scattering which gives rise to finite heat conductivity . for the chain with periodic nearest - neighbor interaction
it was demonstrated @xcite that the transition to normal heat conductivity corresponds to abrupt growth of concentration of rotation solitons ( rotobreathers ) , demonstrating certain similarity with phase transition .
namely , the region of the transition temperature corresponds to maximum region of the heat capacity of the lattice .
similarly , it is reasonable to suppose that every lattice with finite heat conductivity has its peculiar mechanism of scattering the heat flow .
the paper is devoted to the detailed simulation of the discrete lattices with on - site nonlinearity and quadratic potential of nearest - neighbor interaction and investigation of their heat conductivity .
the lattices are frenkel - kontorova , sinh - gordon , and discrete @xmath0 - 4 . for every case the dependence of the heat conductivity on the temperature and parameters of the lattice will be explored and concrete elementary excitation responsible for the change of regimes will be revealed .
let us consider one - dimensional atomic chain arranged along @xmath6 axis .
all particles are of equal mass @xmath7 , and the nearest - neighbor interaction is described by harmonic potential having rigidity @xmath8 .
then the hamiltonian of the lattice will take a form @xmath9 where the dot denotes the differentiation with respect to time @xmath10 , @xmath11 is the displacement of the @xmath12-th particle from its equilibrium position and @xmath13 is the on - site potential .
the dimensionless variables are introduced as @xmath14 ( @xmath15 is the equilibrium distance between the particles ) for the displacement , @xmath16 for the time and @xmath17 for the energy .
hamiltonian ( [ f1 ] ) takes the form @xmath18 where the stroke denotes the differentiation with respect to the dimensionless time @xmath19 and the dimensionless on - site potential is introduced as @xmath20 .
natural definition for the dimensionless temperature is @xmath21 , where @xmath22 is boltzmann constant and @xmath23 is the temperature in common units .
we are going to consider four widely used models for on - site potential : harmonic potential @xmath24 sin - gorgon potential @xmath25 ;         \label{f4}\ ] ] @xmath0 - 4 potential @xmath26 ^ 2         \label{f5}\ ] ] and sinh - gordon potential @xmath27 .
\label{f6}\ ] ] parameter @xmath28 determines the value of potential barrier between neighboring wells and its inverse @xmath29 characterizes the cooperativeness of the system .
potentials ( [ f4 ] ) and ( [ f5 ] ) have the same distance between neighboring wells equal to 2@xmath30 and equal value of the potential barrier 2@xmath31 .
parameter @xmath32 in ( [ f3 ] ) and ( [ f6 ] ) corresponds to minimal frequency of harmonic vibrations of the lattice .
the goal is to simulate the process of heat conduction in finite chain containing n particles . for this sake
the left side of the chain ( @xmath33 ) has to be connected to a thermostat with temperature @xmath34 , and the right side ( @xmath35 )  to a thermostat with temperature @xmath36 ( @xmath37 ) . for the sake of the simulation we consider the chain of @xmath38 particles , where the first @xmath39 particles are attached to the thermostat @xmath34 , an the last @xmath40 particles  to the thermostat @xmath36 ( fig .
[ fig1 ] ) .
the potential of the nearest - neighbor interaction is harmonic , therefore the equilibrium length of the chain does not depend on the temperature .
it implies that the boundary conditions at the ends of the lattice has no noticeable effect on the process of the heat conduction and both the conditions of free ( fig.[fig1 ] ( a ) ) and fixed ( fig .
[ fig1 ] ( b ) ) may be used .
numerical simulations with @xmath41 has demonstrated the absence of the effect depending on concrete choice of the boundary conditions .
we ll use the condition of free ends with @xmath41 for all simulations .
model of the chain of @xmath42 particles with left @xmath39 particles attached to @xmath43 thermostat and right @xmath40 particles attached to @xmath44 thermostat .
boundary conditions correspond to free ( a ) and fixed ( b ) end particles .
the potential @xmath45 corresponds to discrete frenkel - kontorova model . , scaledwidth=85.0% ]    majority of papers devoted to the topic of the heat conduction @xcite uses deterministic nose - hoover thermostat @xcite with @xmath46 .
however this thermostat has been designed for the description of the equilibrium thermalized system and is not universally suitable for the description of non - equilibrium processes .
therefore our choice is well - known stochastic langevin thermostat .
detailed comparison of these two thermostats is presented in appendix [ pr3 ] .
let us consider the chain with free ends ( @xmath47 ) with @xmath48 particles at both ends attached to langevin thermostats .
the dynamics of the system is described by equations @xmath49 where @xmath50 , the damping coefficient @xmath51 , @xmath52 is the characteristic relaxation time of the particles attached to the thermostat , @xmath53 is the random external force corresponding to gaussian white noise normalized as @xmath54 details of numeric realization of the langevin thermostat and random forces are presented in appendix [ pr1 ] .    at every moment
the dimensionless temperature of the @xmath12-th particle @xmath55 . in order to determine the value of the local heat flux
@xmath56 the energy distribution among the particles of the chain is considered : @xmath57 .
\label{f8}\ ] ] by differentiating equation ( [ f8 ] ) with respect to time @xmath19 we get @xmath58             + u'_{n+1}[u''_{n+1}+f(u_{n+1})]\right\ }             + ( u_{n+1}-u_n)(u'_{n+1}-u'_n).\ ] ] with account of ( [ f7 ] ) we obtain @xmath59 .
\label{f9}\ ] ] taking into account the continuity condition @xmath60 we get the expression for the energy flux : @xmath61 .
system of equations ( [ f7 ] ) has been integrated numerically .
we used the values of @xmath62 , @xmath41 , @xmath63 , 20 , 40 , 80 , 160 , 320 , 640 and initial conditions corresponding to the ground state of the chain .
after the time @xmath64 has elapsed , the end particles achieved thermal equilibrium with the thermostat and stationary heat flux has been formed .
afterwards the dynamics of system ( [ f7 ] ) has been simulated at the time scale of order @xmath65 .
the average temperatures of the particles @xmath66 and average values of the heat flux @xmath67 were computed for the fragment of the chain between the thermostats .
if the temperature gradient @xmath68 is small , this method allows avoiding the temperature jumps at the ends of the free fragment of the chain @xcite .
characteristic distributions of the heat flux @xmath69 and local temperature @xmath70 are demonstrated at fig .
[ fig2 ] . at the inner fragment of the chain @xmath71 the heat flux is constant @xmath72 and the temperature profile is linear .
the coefficient of the heat conductivity is determined using the information concerning the inner fragment of the chain : @xmath73 if the function @xmath74 is the best linear approximation for @xmath75 data at the inner fragment of the chain @xmath71 , then the value @xmath76 may be calculated with better accuracy by taking @xmath77 .
distribution of the local heat flux @xmath69 ( a ) and local temperature @xmath70 ( b ) in the chain with periodic on - site potential ( [ f4 ] ) , @xmath78 , @xmath79 , @xmath41 , @xmath80 , @xmath81 .
time of averaging @xmath65 .
fragments of the chain interacting with the thermostats are embedded in gray .
, scaledwidth=85.0% ]    limit value @xmath82 will correspond to the coefficient of the heat conductivity at temperature @xmath83 .
the question regarding the finiteness of the heat conductivity is reduced to the existence of finite limit ( [ f13 ] ) .
i the sequence @xmath84 diverges as @xmath4 then the chain has infinite heat conductivity at this value of the temperature .
alternative way to compute the heat conductivity @xmath85 is related to well - known green - kubo formula @xcite : @xmath86 where @xmath3 is already the number of particles in the chain with periodic boundary conditions ,    @xmath87 is the general heat flux and the averaging @xmath88 is performed over all thermalized states of the chain .
consequently , the finiteness of the heat conductivity is related to the convergence of the integral @xmath89 e.g. to the evaluation of the descending rate of the function @xmath90    numerically the above autocorrelation function may be found only for finite chain @xmath91 for large enough values of @xmath3 the correlation function @xmath92 is believed to approximate the function @xmath93 with acceptable accuracy . in order to get stable results the value @xmath94 is usually sufficient .
more details concerning the computation of the autocorrelation function are presented in appendix [ pr2 ] .
the methods for computing the heat conductivity coefficient ( [ f13 ] ) , ( [ f14 ] ) are complementary and allow mutual verification of the results .
the chain with harmonic on - site potential ( [ f3 ] ) is described by linear equations and therefore is completely integrable .
the energy transport is performed by non - interacting phonon modes .
the heat flux @xmath95 does not depend on the chain length @xmath3 , but only on the temperature difference @xmath96 .
linear thermal profile is not formed . at the inner part of the chain
the temperature is nearly constant @xmath97 ( fig .
[ fig3 ] ) . therefore according to ( [ f12 ] ) , the heat conductivity coefficient diverges . correspondingly , the average correlation function @xmath98 is constant and integral ( [ f15 ] ) diverges .
distribution of local heat flux @xmath69 ( a ) and local temperature @xmath70 ( b ) in the chain with harmonic on - site potential ( [ f3 ] ) , @xmath99 , @xmath100 , @xmath41 , @xmath80 , @xmath81 , averaging time @xmath65 .
the fragments of the chain interaction with the thermostats are embedded in grey .
thin lines ( 1 , 3 ) are obtained by using the nose - hoover thermostat with @xmath101 , and thick ( 2 , 4 )  by using the langevine thermostat with @xmath102 .
, scaledwidth=85.0% ]
characteristic features of dynamics of the chain with periodic on - site potential ( [ f4 ] ) depend on the values of the temperature .
as the temperature is small @xmath103 , the on - site potential may be approximated by harmonic single - well potential ( [ f3 ] ) with @xmath104 .
the heat transport is governed by weakly interacting phonons . at the temperature
@xmath105 the chaotic superlattice of topological solitons is formed and the transport properties change drastically . at very high temperatures @xmath106 the chain
is effectively detached from the site and again weakly interacting phonons govern the heat transport . therefore it is reasonable to investigate the dependence of the heat conductivity on the reduced temperature @xmath107 .
the behavior of the chain also depends on the cooperativeness parameter @xmath108
. the more the cooperativeness , the less is the density of the soliton superlattice and the phonon scattering effects are less significant .
the limit @xmath109 ( @xmath110 ) corresponds to completely integrable continuum sin - gordon equation .
generally , three limits of discrete frenkel - kontorova system correspond to completely integrable systems : at @xmath111 the system reduces to the harmonic chain with harmonic on - site potential ; at @xmath112  to isolated harmonic chain ; at @xmath109  to continuum sin - gordon equation .
all these limit systems have diverging heat conductivity .
the behavior of the system in the vicinity of these limits is a natural question to be addressed .
dependence of the logarithm of the heat conductivity coefficient @xmath113 ( a ) and @xmath76 ( b ) on the logarithm of the inner fragment length @xmath114 ( @xmath115 ) for the chain with periodic on - site potential ( [ f4 ] ) , @xmath78 , @xmath116 ( markers 1 ) , @xmath117 ( markers 2 ) , @xmath118 ( markers 3 ) and @xmath119 ( markers 4 ) .
the markers denote the computed values and the lines correspond to the best linear approximations .
, scaledwidth=85.0% ]    let us start from @xmath120 @xmath121 , @xmath122 and investigate the sequence @xmath76 as @xmath3 grows @xmath123 20 , 40 , 80 , 160 , 320 , 640 ) and different values of @xmath124 .
as it is may be suggested from fig .
[ fig4 ] at small ( @xmath116 ) and large ( @xmath117 ) temperatures the heat conductivity coefficient @xmath76 grows as @xmath125 , at @xmath119  as @xmath114 , and at @xmath118 converges to finite value @xmath126 .
therefore it may be concluded that at @xmath118 the chain has finite heat conductivity .
the data related to the other values of the temperature does not allow to draw any conclusions about the behavior of the heat conductivity at larger values of @xmath3 .
generally speaking , it may happen that for longer chains @xmath84 will attain certain finite value .
computational tools we use do not allow to investigate higher values of @xmath3 .
stil , it is possible to get additional information from the behavior of the autocorrelation function @xmath98 at @xmath127 .     exponential decrease of the autocorrelation function @xmath98 in the chain with periodic on - site potential ( [ f4 ] ) , @xmath78 , @xmath118 ( curve 1 ) and @xmath119 ( curve 2 ) ( semilogarythmic coordinates ) . ,
scaledwidth=85.0% ]     autocorrelation functions @xmath98 in the chain with periodic on - site potential ( [ f4 ] ) , @xmath78 , @xmath116 ( curve 1 ) and @xmath117 ( curve 2 ) . ,
scaledwidth=85.0% ]     the zone in the space of parameters @xmath128 , where for finite chains of length @xmath129 with on - site potential ( [ f4 ] ) he heat conductivity converges ( a , grey zone ) and diverges ( b , white zone ) .
curve 1 divides these two zones .
interval 2 corresponds to the parameters used in @xcite . for finite chains
@xmath130 with on - site potential @xmath04 ( [ f5 ] ) finite heat conductivity is detected only above line 3 . , scaledwidth=85.0% ]    numerical simulation demonstrates that for @xmath118 the autocorrelation function decreases exponentially ( fig .
[ fig5 ] , curve 1 ) .
integral ( [ f15 ] ) converges and the green - kubo formula ( [ f14 ] ) gives @xmath131 , in good correspondence with @xmath126 obtained from direct simulation of the heat flux . at @xmath119
the autocorrelation function at time scale @xmath132 also decreases exponentially ( fig .
[ fig5 ] , curve 2 ) .
if this trend will persist also for @xmath133 , the green - kubo formula will give @xmath134 .
it is reasonable to compare this value with the result for @xmath84 presented at fig .
[ fig4 ] ( curve 4 ) . maximum value of @xmath135 and no trend towards any finite limit of @xmath84 may be detected .
therefore the likely result is divergence . in order to verify this result the simulation for larger values of @xmath3 ( 1280 , 2560 , 5120 , 10240 )
is required , which is beyond our computational possibilities .
the problem for @xmath116 and @xmath117 is even more difficult .
the autocorrelation function is presented at fig .
the decrease of the function is very slow and no unambiguous conclusion concerning its character may be drawn out . while extrapolating @xmath136 for @xmath137 by exponent ,
the green - kubo formula yields @xmath138 for @xmath116 and @xmath139 for @xmath117 . in order to get additional information
another consuming simulation is required . still , from the other side , for @xmath117 at @xmath140 the logarithm of the heat conductivity @xmath141 , and the dependence @xmath142 ( fig .
[ fig4 ] , curve 2 ) does not demonstrate any trend towards convergence .
therefore the most likely result in this case is also the divergence of the heat conductivity .
let us consider the sequence @xmath84 ( @xmath63 , 20 , 40 , 80 , 160 , 320 , 640 ) at other values of the cooperativeness .
the results are summarized at fig .
the space of parameters @xmath128 is divided to two zones denoted by different colors . in the first ( gray ) zone
the sequence @xmath84 converges ( @xmath143 ) , and in the second ( white ) zone the sequence grows monotonously .
then , in the first zone frenkel - kontorova model has finite heat conductivity , and in the second zone the heat conductivity is either divergent or finite but very high ( for @xmath129 the fourier law is not valid ) .
the first zone is limited by certain finite value of @xmath144 : for some @xmath145 and for all @xmath146 no convergence of @xmath84 was detected .
the explanation is that for growing @xmath144 the system becomes closer to continuum integrable sin - gordon equation . at any fixed @xmath147 for @xmath129
the heat conductivity converges only for some finite temperature interval @xmath148 .
as the cooperativeness decreases ( @xmath149 ) , the upper boundary of this interval tends to infinity ( @xmath150 ) , and the lower boundary tends to zero ( @xmath151 ) proportionally to @xmath144 .
dependence of the heat conductivity coefficient @xmath85 from the reduced temperature @xmath107 for the chain with periodic on - site potential ( [ f4 ] ) for @xmath152 ( curve 1 ) , @xmath153 ( curve 2 ) and @xmath154 ( curve 3 ) .
, scaledwidth=85.0% ]    the dependence of @xmath85 on the reduced temperature @xmath155 is presented at fig .
[ fig8 ] . within the interval
@xmath156 $ ] there exists a critical value @xmath157 corresponding to the minimum of heat conductivity .
the dependence of the dimensionless heat capacity @xmath158 on the reduced temperature @xmath107 ( a ) for the chain with periodic on - site potential ( [ f4 ] ) and ( b ) for the chain with @xmath04 potential ( [ f5 ] ) for @xmath154 ( curves 1,6 ) , @xmath153 ( curves 2 , 7 ) , @xmath152 ( curves 3 , 8) , @xmath78 ( curves 4 , 9 ) and @xmath159 ( curves 5 , 10 ) . dashed curve 11 gives similar dependence for the chain with on - site sinh - gordon potential ( [ f6 ] ) , for @xmath99 .
, scaledwidth=85.0% ]    in order to reveal the mechanism of the heat conduction it is reasonable to explore the behavior of heat capacity @xmath160 ( @xmath161 is the average energy of cyclic @xmath3-atomic chain at the temperature @xmath124 ) on the reduced temperature @xmath155 ( fig .
[ fig9 ] ) .
the heat capacity of classic harmonic chain is always unity , therefore the discrepancy of this value from unity characterizes the significance of nonlinear effects at given temperature .
the lattice considered has negative anharmonism and therefore its heat capacity must be more than unity for all temperature diapason .
the heat capacity tends to unity as @xmath111 and @xmath112 and has single maximum at certain temperature @xmath162 .
this value fairly coincides with the temperature @xmath157 , which corresponds to the minimum of the heat conductivity .
moreover , the increase and decrease of the heat capacity is clearly correlated with the decrease and increase of the heat conductivity respectively .
this fact allows suggesting the same physical effects as responsible for both processes . for zero temperature
the heat capacity is equal to unity .
the increase of the heat capacity at higher temperatures is related to thermal activation of topological solitons ( kinks and antikinks ) which represent additional degrees of freedom for this system . as a result the dynamical superlattice of solitons appears .
the density of this superlattice approaches its maximum at the temperature @xmath157 .
further growth of the temperature results in the decrease of the number of degrees of freedom , which is manifested as effective detaching of the chain from the on - site potential .
therefore the heat capacity decreases and tends to unity as the temperature grows .
zones in the space of parameters @xmath128 , where for finite chains @xmath129 with periodic on - site potential ( [ f4 ] ) and period @xmath163 the heat conductivity is normal ( a , gray ) and abnormal ( b , white ) . dashed line
denotes the same boundary for commensurate frenkel - kontorova model ( @xmath164 ) .
, scaledwidth=85.0% ]    correlations between the behavior of the heat capacity and the heat conductivity and especially fair coincidence of @xmath157 and @xmath162 allow us to suppose that the heat transfer is limited by phonon scattering on the soliton superlattice .
the effectiveness of such scattering depends on the density of the superlattice as well as on the ability of single kink to scatter phonons . in the strongly cooperative regime
@xmath146 the interaction between solitons and phonons is nearly elastic ( close to the case of complete integrability ) and therefore the heat conductivity has the trend to diverge . for lower cooperativeness
the soliton - phonon interaction is less elastic and finite temperature diapason @xmath156 $ ] of converging heat conductivity appears . for the cases of low @xmath165 and high @xmath166 temperatures
the convergence of the heat conductivity can not be detected in the framework of current experiment .
the suggested reason of this effect is that the soliton superlattice effectively disappears .
let us consider now incommensurate frenkel - kontorova chain where the period of the chain is different from the period of on - site potential .
the dimensionless on - site potential is periodic function([f4 ] ) with period @xmath167 , and the chain has period @xmath168 .
then in system of equations ( [ f7 ] ) function @xmath169 will take the form @xmath170 for the sake of simulation we choose @xmath171 , corresponding in certain sense to extremely incommensurate case .
it is well - known @xcite that such a lattice in its ground state already has soliton superlattice of nonzero density .
therefore the convergence of the heat conductivity is expected to be facilitated as compared to the commensurate case .
[ f10 ] demonstrates the zone in the space of parameters @xmath128 where the sequence @xmath84 , @xmath63 , 20 , 40 , 80 , 160 , 320 , 640 .
converges . for the sake of comparison the boundary for the commensurate case
is also presented ( @xmath164 ) .
the result is than no qualitative change of the behavior occurs .
the only difference is that the zone with normal heat conductivity moves downwise .
this effect is related to presence of superlattice of solitons at any temperature .
the transition to normal heat conduction occurs at lower temperature since less solitons should be thermally activated in order to achieve convergence . from the other side ,
the soliton superlattice facilitates effective detachment of the lattice from the on - site potential ( the average coupling energy in the ground state is less ) and therefore the upper boundary for the normal heat conduction is also achieved at lower temperatures .
let us consider the heat conductivity of the chain with on - site potential @xmath0 - 4 ( [ f5 ] ) .
for this case the analysis of the sequence @xmath84 , @xmath63 , 20 , 40 , 80 , 160 , 320 , 640 , demonstrates that the heat conductivity converges as @xmath172 ( @xmath173 ) - see fig .
[ fig7 ] .     exponential decrease of the autocorrelation function @xmath98 in the chain with on - site potential @xmath04 ( [ f5 ] ) , @xmath78 , @xmath119 ( curve 1 ) , @xmath174 ( curve 2 ) , @xmath175 ( curve 3 ) and @xmath118 ( curve 4 ) .
( semilogarythmic coordinates @xmath176 versus @xmath19 ) .
, scaledwidth=85.0% ]    in order to investigate the character of the heat conduction in the temperature range @xmath177 let us consider the temperature behavior of the autocorrelation function @xmath98 . for @xmath120 ( @xmath78 )
this behavior is demonstrated at fig .
[ fig11 ] . as @xmath127 the autocorrelation function decreases exponentially .
the decrease rate grows as the temperature increases and therefore the conclusion concerning finite heat conductivity at @xmath178 is confirmed . at lower temperatures
the decrease rate satisfies a power law @xmath179  see fig .
[ fig12 ] .
the degree @xmath180 decreases with the decrease of the temperature . at @xmath181 @xmath182 ,
therefore integral ( [ f15 ] ) converges and the heat conductivity is finite , at @xmath183 @xmath184 . within the accuracy available for current numerical possibilities this value corresponds to transition to abnormal heat conduction .
it is extremely difficult to obtain reliable data for lower temperatures in order to substantiable this conclusion because of huge computation time required .
the reason is that the system is rather close to completely integrable case .
new numerical methods based on the latter fact are desirable for investigation of this kind of systems .
power - law decrease of the autocorrelation function @xmath98 in the chain with @xmath04 on - site potential ( [ f5 ] ) , @xmath78 , @xmath181 ( curve 1 ) and @xmath183 ( curve 2 ) . ( double logarithmic coordinates , @xmath176 versus @xmath185 ) .
the angle coefficient @xmath180 determines the decrease rate . for @xmath181 @xmath186 , for @xmath183 @xmath184 .
, scaledwidth=85.0% ]    the dependence of the heat conductivity @xmath85 on reduced temperature @xmath155 is presented at fig .
[ fig13 ] . for low cooperativeness ( @xmath187 )
the heat conductivity approaches local minimum and afterwards local maximum at @xmath188 and monotonously decreases to zero as @xmath112 .
the relative value of the maximum decreases as the cooperativeness grows and disappears for certain critical value of @xmath31 .
heat conductivity @xmath85 versus reduced temperature @xmath107 for the chain with double - well on - site potential ( [ f5 ] ) , @xmath189 ( curve 1 ) , @xmath190 ( curve 2 ) and @xmath78 ( curve 3 ) and for the chain with sinh - gordon on - site potential ( [ f6 ] ) , @xmath99 ( dashed curve 4 , @xmath191 ) . ,
scaledwidth=85.0% ]    in order to reveal the physical reasons of such behavior of the heat conductivity it is also reasonable to investigate the behavior of the heat capacity @xmath158 ( fig .
[ fig9]b ) . as @xmath111 the heat capacity @xmath192 .
as the temperature grows , the heat capacity grows , achieves its maximum at the temperature @xmath162 and then decreases monotonously to the value less than unity . the value @xmath162 is situated near the maximum point of the heat conductivity @xmath193 .
such behavior is related with the peculiarities of @xmath04 potential . at low temperatures
the main effect is due to negative anharmonism near the ground state ( therefore the heat capacity exceeds unity ) and for high temperatures ( @xmath194 ) the process is governed by positive anharmonism bringing the heat capacity to the value below unity .
let us now consider the frequency spectrum of vibrations of the chain .
the spectrum is computed for @xmath189 @xmath195 and three characteristic temperatures @xmath196 , 10 , 100 .
the spectrum of the chain with harmonic on - site potential ( [ f3 ] ) does not depend on the temperature and has the form @xmath197 where maximum frequency @xmath198 . for @xmath189 , @xmath199 , @xmath200 . as it is demonstrated at fig .
[ fig14]a , for temperature @xmath196 the spectrum of the chain with on - site @xmath04 potential nearly coincides with the spectrum of purely harmonic chain ( [ f17 ] ) .
such spectrum means that at low temperatures only thermalized phonons contribute to the frequency spectrum and other excitations do not play any significant role .
for @xmath201 the distribution crosses the lower boundary of the propagation zone @xmath32 ( fig . [
fig14]b ) .
such low - frequency component may be associated with intrinsic vibrations of the solitons superlattice . for even higher temperatures
@xmath202 the spectrum crosses also the upper boundary of the propagation zone @xmath203 ( fig . [
. such effect may be attributed only to thermalization of high - frequency discrete breathers .
frequency spectrum of energy of vibrations in the chain with double - well on - site potential ( [ f5 ] ) at temperatures @xmath196 ( a ) , @xmath174 ( b ) and @xmath204 ( c ) .
dashed line denotes the spectrum of harmonic chain ( [ f3 ] ) with @xmath205 .
, scaledwidth=85.0% ]    therefore , for low temperatures @xmath206 the dynamics of the system is close to that of harmonic chain .
the heat transport is governed by weakly interacting phonons and heat conductivity diverges . for higher temperatures the heat conductivity converges . in the intermediate diapason
@xmath207 the effective phonon scattering mechanism exists due to the superlattice of topological solitons , and for high temperatures @xmath208  due to high - frequency discrete breathers .
interplay of two different mechanisms of the phonon scattering explains also the dependence of the heat conductivity on the cooperativeness of the system ( fig .
[ fig13 ] ) .
the minimum and maximum of heat conductivity disappear with growth of the cooperativeness since the soliton mechanism of scattering becomes less effective ( the soliton - phonon interaction is closer to elastic ) and simultaneously the excitation of the discrete breathers becomes easier .
heat conductivity of this system has been investigated in paper @xcite and here it is reasonable to elucidate the details related to physical mechanisms of the process .
the on - site potential ( [ f6 ] ) is single - well function with positive anharmonism .
the analysis of the finite sequence @xmath84 , @xmath63 , 20 , 40 , 80 160 , 320 , 640 , demonstrates that the heat conductivity converges for high temperatures ( @xmath209 ) .
this observation is supported by the fact that the autocorrelation function @xmath98 at high temperatures for @xmath127 decreases exponentially ( fig .
[ fig15 ] ) , and for low temperatures - by power law ( fig .
[ fig16 ] ) .     exponential decrease of the autocorrelation function @xmath98 in the chain with sinh - gordon on - site potential ( [ f6 ] ) , @xmath99 , @xmath174 ( curve 1 ) , @xmath210 ( curve 2 ) and @xmath175 ( curve 3 ) .
( semilogarythmic coordinates @xmath176 versus @xmath19 ) . , scaledwidth=75.0% ]     power - law decrease of the autocorrelation function @xmath98 in the chain with sinh - gordon on - site potential ( [ f6 ] ) , @xmath99 , @xmath211 .
solid line corresponds to the number of particles @xmath212 , dotted  to @xmath213 , dashed - dotted  to @xmath214 .
, scaledwidth=75.0% ]    the heat conductivity decreases monotonously and for @xmath215 exponentially tends to zero ( fig .
[ f13 ] , curve 4 ) .
positive anharmonism of the potential leads to monotonous decrease of the heat capacity ( fig .
[ f9 ] , curve 11 ) .
the frequency spectrum of vibrations moves towards the upper boundary of the propagation zone with growth of the temperature .
these facts allow concluding that the high - frequency discrete breathers provide effective phonon scattering in this model and facilitate the convergence of the heat conductivity .
growing concentration of these breathers with the growth of the temperature leads to monotonous decrease of the heat conductivity coefficient
.    chain with on - site potential @xmath216 ( positive @xmath5 model ) also has finite heat conductivity @xcite .
potential ( [ f18 ] ) as well as sinh - gordon on - site potential ( [ f6 ] ) is single - well symmetric function with positive anharmonism .
therefore the mechanism of the phonon scattering is also related to the discrete breathers and @xmath217 for @xmath215 .
for @xmath218 the heat conductivity @xmath219 @xcite .
the investigation presented above demonstrates that the anharmonicity of the on - site potential does not constitute sufficient condition for the convergence of the heat conductivity coefficient .
the behavior of any concrete model in the above respect depends on its peculiar nonlinear excitations which determine the process of the heat transfer and phonon scattering .
two typical mechanisms of the phonon scattering were revealed in the paper  thermalized soliton superlattice ( discrete sin - gordon and @xmath04 models ) and discrete high - frequency breathers ( @xmath04 and sinh - gordon models ) .
phonon scattering mechanism may switch with the change of the temperature ( @xmath04  model ) .
for the discrete frenkel - kontorova model the zone of the converging heat conductivity for given chain length is limited by low and high temperatures and by high cooperativeness .
the numerical possibilities available up to date does not allow establishing unambiguously the character of the heat conductivity outside the zone designated at fig .
still there is a reason to suggest that infinite chain for certain parameters has diverging heat conductivity , although the zone corresponding to finite heat conductivity will be larger that computed above .    unlike frenkel - kontorova model , for @xmath04 model
it is possible to demonstrate that for low temperatures the boundary of the transition to abnormal heat conduction may be achieved .
it is possible to suppose that there exist a transition from infinite to finite heat conductivity with growth of temperature for any cooperativeness .
the probable reason for the divergence to be detectable is the presence of odd - power terms in the on - site potential in the vicinity of the extrema of the potential wells .
the sinh - gordon model does not allow to detect the divergence of the heat conductivity in current experiment ; still , the transition also may be suggested for any cooperativeness .
it is possible to suggest that for any analytic on - site potential for low temperatures the heat conductivity will diverge .
the authors are grateful to russian foundation of basic research ( grant 01 - 03 - 33122 ) , to ras commission for support of young scientists ( 6th competition , grant no .
123 ) and to fund for support of young scientists for financial support .
savin is grateful to international association of assistance for the promotion of co - operation with scientists from the new independent states of the former soviet union ( project intas no .
96 - 158 ) for financial support .
system of equations describing the dynamics of the chain attached to thermostats ( [ f7 ] ) has been integrated numerically by standard fouth - order runge - kutta method with constant step of integration @xmath220 .
numeric realization of delta - function is performed as @xmath221 for @xmath222 and @xmath223 for @xmath224 , i.e. the step of integration corresponds to the correlation time of the random forces .
that is why in order to get correct description of the langevin thermostat we must guarantee that the relaxation time @xmath225 . in order to fulfill this condition the relaxation time
was chosen as @xmath102 , and the step of integration for different values of @xmath3 , was chosen as @xmath226 , 0.025 , 0.0125 .
for every step of integration the random forces @xmath53 were taken to be constant .
they were computed as independent realizations of the random value @xmath227 , normally distributed with zero average @xmath228 and dispersion @xmath229 . for generating the random value @xmath227 program package zufall @xcite
was used .
the initial state for the integration of equations ( [ f7 ] ) was chosen to be equal to ground state of the chain : @xmath230 where @xmath231 for on - site potentials ( [ f3 ] ) , ( [ f6 ] ) and @xmath232 for potentials ( [ f4 ] ) , ( [ f5 ] ) .
it is convenient to control the accuracy of the simulation through behavior of sequence of average local heat fluxes @xmath233 .
if the choice of the integration step @xmath220 is correct then this sequence should be constant . if the local average heat flux changes from particle to particle then the integration step should be reduced . for growing chain length @xmath3 the step of integration
should be also reduced in order to provide sufficient accuracy ; the averaging time also grows ( see @xcite ) and therefore the time of simulation necessary for obtaining reliable results for large @xmath3 turns out to be extremely large .      in order to compute the autocorrelation function of the heat flux @xmath92 dynamics of cyclic @xmath3-particle chain
was simulated .
the thermalized chain with temperature @xmath124 was obtained by integrating langevin system of equations @xmath234 where @xmath235 for @xmath236 and @xmath237 for @xmath238 , @xmath62 ( relaxation time @xmath102 ) , @xmath239  white gaussian noise normalized as @xmath240 system ( [ f20 ] ) has been integrated numerically with initial conditions corresponding to the ground state of the chain .
after time @xmath241 the chain approached equilibrium with the thermostat and the coordinates @xmath242 corresponding to the thermalized state at temperature @xmath124 .
afterwards the dynamics of isolated thermalized chain was simulated . for this sake system ( [ f20 ] ) was integrated with zero damping @xmath243 and zero external force @xmath244 ) .
thermalized state ( [ f21 ] ) was used as initial condition .
the result was the dependence of the general heat flux * j * on time @xmath19 . afterwards with the help of ( [ f16 ] ) the autocorrelation function @xmath92 was computed for given thermalized state of the chain .
the autocorrelation function depends significantly on concrete realization of the thermalized chain .
that is why in order to improve the accuracy this procedure was performed @xmath245 with independent initial realizations of the thermalized state .
finally the shape of the correlation function was computed as average over all these realizations .
it is worth while mentioning that the alternative way of computation ( performing of one very large simulation ) would not bring about any sufficient gain in the accuracy because of growing integration errors .    in order to verify the independence of the correlation function on the chain length
the appropriate calculations were performed for different values of @xmath3 .
[ fig16 ] demonstrates the function @xmath92 for the chain with sinh - gordon on - site potential for @xmath99 , @xmath211 and @xmath212 , 1000 , 2000 .
it is clear that the autocorrelation function is nearly independent on @xmath3 ( the differences are noticeable only for large times and are reduced as the number of realizations used for averaging grows ) . for given set of parameters @xmath213 provides sufficient accuracy .
unlikely the langevin thermostat , the nose - hoover thermostat ( nht)@xcite is not stochastic .
its dynamics is completely determined by the initial conditions .
it turned out to be good choice for simulations of fpu system @xcite but its deterministic nature can bring about artifacts in the behavior of the system .
we compare this thermostat with the langevin thermostat ( lt ) we use for the case of frenkel - kontorova model .
let us consider the chain with fixed ends ( @xmath47 ) with @xmath48 particles attached to nht having the temperature @xmath246 .
dynamics of the system is described by equations @xmath247 where @xmath50 , and @xmath52 is the relaxation time of the thermostat .
distribution of local heat flux @xmath69 ( a ) and and local temperature @xmath70 ( b ) in the chain with periodic on - site potential ( [ f4 ] ) , @xmath78 , @xmath100 , @xmath41 , @xmath248 , @xmath249 , averaging time @xmath65 .
grey zones denote the chain fragments embedded in the thermostats .
thick lines ( 1 , 3 ) correspond to nht ( @xmath101 ) , and thin ( 2 , 4 )  to lt ( @xmath102 ) . , scaledwidth=85.0% ]     frequency distribution of the energy of particle having number @xmath250 in the chain with periodic on - site potential ( [ f4 ] ) , @xmath78 , @xmath100 , @xmath41 , @xmath248 , @xmath249 .
thick line corresponds to use of nht ( @xmath101 ) , and thin  to use of lt ( @xmath102 ) .
, scaledwidth=85.0% ]    usually the simulations of the heat conductivity @xcite take @xmath101 , and @xmath251 ( only end particles are attached to the thermostat @xmath252 and @xmath253 ) .
but , as stated in @xcite , such thermostats are not enough random
they cover only a part of the phase space and correspond to strange attractors . in order to reduce this effect
we attach to the thermostat @xmath115 particles from every side of the chain .
the dynamics of system ( [ f22 ] ) is also completely deterministic .
it should be mentioned that it is impossible to use the initial condition ( [ f19 ] ) corresponding to ground state of the system ( it is stationary point of system ( [ f22 ] ) ) .
we take the initial condition @xmath254 where @xmath239
 independent realizations of the random variable over the interval [ 0,1 ] .
we choose @xmath78 ( @xmath120 ) , @xmath248 , @xmath249 , @xmath255 and integrate system ( [ f22 ] ) numerically with initial condition ( [ f23 ] ) .
the distribution of heat fluxes @xmath69 and local temperatures @xmath70 is presented at fig . [ fig17 ] ( for the sake of comparison we present also the results received bu using lt - thin lines ) . within the left thermostat the heat flux grows linearly and within the other thermostat it decreases linearly with @xmath12 . at central part of the chain the value of the heat flux
does not depend on @xmath12 .
linear temperature profile is formed and the heat conductivity coefficient may be computed according to ( [ f12 ] )  @xmath256 .
use of lt gives @xmath257 ( see above ) , i.e. the value of @xmath85 does not depend on the type of the thermostat .
in addition , it is possible to conclude from fig .
[ fig18 ] that the frequency distribution of the energy of vibrations also does not depend on the type of thermostat used .
it means that for the case of the temperatures close to the value of the potential barrier the choices of nht or lt bring about equivalent results    the situation is strikingly different if the temperature is lower and the chain is closer to the linear case .
the nose - hoover thermostat is not effective in this case . in order to illustrate this fact we use the model of harmonic chain . as it is clear from fig .
[ fig3 ] nht gives values of the heat flow substantially different from the correct values ; at the same times the use of lt secures much better results .
that is why in the present paper we used more complicated and consuming lt .
99 e.  fermi , j.  pasta , and s.  ulam , los alamos rpt la-1940 , 1955 .
s.  lepri , l.  roberto , and a.  politi , phys .
lett . * 78 * , 1896 ( 1997 ) .
s.  lepri , r.  livi , and a.  politi , physica d * 119 * , 140 ( 1998 ) .
s.  lepri , r.  livi , and a.  politi , evrophys . lett . * 43 * , 271 ( 1998 ) .
r.  rubin and w.  greer , j. math .
* 12 * , 1686 ( 1971 ) .
a.  casher and j.l .
lebowitz , j. math .
( n.y . ) * 12 * , 1701 ( 1971 )
. a.  dhar , phys .
lett . * 86 * , 5882 ( 2001 ) .
a.  dhar , phys .
lett . * 86 * , 3554 ( 2001 ) .
savin , g.p .
tsironis , and a.v .
zolotaryuk , phys .
* 88 * , ( 15)/154301 ( 2002 ) .
g.  casati and t.  prosen , arxiv : cond - mat/0203331 . t.  hatano , phys . rev .
e * 59 * , r1 ( 1999 ) .
g.  casati , j.  ford , f.  vivaldi , and v.m .
visscher , phys .
* 52 * , 1861 ( 1984 ) .
t.  prosen , m.  robnik , j. phys .
a * 25 * , 3449 ( 1992 ) .
m.  j.  gillan , r.  w.  holloway , j. phys .
c * 18 * , 5705 ( 1985 ) .
b.  hu , b.  li , and h.  zhao , phys .
e * 57 * , 2992 ( 1998 ) .
tsironis , a.r .
bishop , a.v .
savin , and a.v .
zolotaryuk , phys .
e * 60 * , 6610 ( 1999 ) .
b.  hu , b.  li , and h.  zhao , phys .
e * 61 * , 3828 ( 2000 ) .
k.  aoki and d.  kuznezov , phys .
a * 265 * , 250 ( 2000 ) ; arxiv : chao - dyn/9910015 . f.  bonetto , j.l .
lebowitz and l.  ray - bellet in _
mathematical physics 2000 _ , ed .
a. fokas , a. grigoryan , t. kibble and b. zegarlins s.  lepri , r.  livi , and a.  politi , arxiv : cond - mat/0112193 . c.  giardina , r.  livi , a.  politi , and m.  vassalli , phys .
lett . * 84 * 2144 ( 2000 ) .
o.v .  gendelman and a.v .
savin , phys .
lett . * 84 * , 2381 ( 2000 ) .
a.v .  savin and o.v .
gendelman , fizika tverdogo tela * 43 * , 341 ( 2001 ) [ physics of the solid state * 43 * , 355 ( 2001 ) ] . s. nose , j. chem .
phys . * 81 * , 511 ( 1984 ) ; w. g. hoover , phys .
a * 31 * , 1695 ( 1985 ) .
a.  fillipov , b.  hu , b.  li , and a.  zeltser , j. phys .
a * 31 * , 7719 , ( 1998 ) .
k.  aoki and d.  kusnezov , phys .
* 86 * , 4029 ( 2001 ) .
r.  kubo , m.  toda , and n.  hashitsume , statistical physics ii , springer ser .
solid state sci .
, * 31 * ( 1991 ) . v.l .
pokrovsky , a.l .
talapov , zh .
75 * , 1151 ( 1978 ) .
petersen , `` lagged fibonacci random number generators for the nec sx-3 ''  is to be published in the international journal of high speed computing ( 1994 ) ; ftp://ftp.netlib.org/random/zufall.f a.  dhar , phys . rev
. lett . * 87 * , 069401 ( 2001 ) . m.  terraneo , m.  peyrard , and g.  casati , phys . rev . lett . * 88 * , 094302 ( 2002 ) .