the phenomena of diffusion and transport over a potential energy landscape play a key role in a number of processes in physics , chemistry and biology @xcite .
josephson tunneling junctions , superionic conductors , phase - locked - loop frequency control systems , charge density waves are a few examples of systems in which these processes in periodic potential are important @xcite .    in recent years
the interest has been growing to experimental studies of manipulating the particle diffusion by external fields .
one can effectively control the diffusion processes by varying the field parameters .
for instance , huge enhancement in diffusivity ( exceeding its value in the same system without extra forcing by a factor of order hundreds ) was observed in studying particle diffusion in colloids with optical traps ( optical vortices ) @xcite . large increase in diffusivity
was also observed in studies of paramagnetic particle diffusion on surfaces of garnet ferrites influenced by external time - periodic magnetic field @xcite . in the same manner diffusivity growth with shaking strength was noticed in experiments with granular gas @xcite .
it turned possible to increase the diffusivity of ions in membrane channels by varying the external electromagnetic field @xcite .
one can increase the rates of diffusion - limited physical processes , and effectively separate micro- and nano - particles of different nature by varying diffusion coefficients in different directions @xcite . despite these remarkable achievements , quantitative description of such enhanced diffusion under
the influence of external forces remains fragmentary to date .
special cases of underdamped and overdamped physical systems are often considered . in overdamped systems
inertial effects may be neglected , significantly simplifying the mathematical treatment of the problem . in underdamped systems on the other hand
the viscous decay of oscillations occurs on times long compared with the oscillation period in the potential ; such systems are harder to investigate , and the progress in studying these is more limited compared with overdamped ones .
first detailed studies of a brownian particle motion in a washboard potential ( _ i.e. _ a periodic potential tilted with an additional constant force field ) were carried out by h.  risken and colleagues @xcite , for all values of the friction .
it was shown that at low friction the appearance of locked  ( in which the particle oscillates around one minimum of the potential ) and running  ( where the particle travels through many potential periods , not getting trapped at separate minima ) solutions was important for the particle ensemble behavior . at small temperatures
running solutions appear abruptly at the external force value exceeding certain critical @xmath3 ; at smaller forces a single class of solutions ( locked ) is realized .
systematization of the results of these studies may be found in authoritative monograph @xcite .
approximate solution of the fokker - planck equation ( fpe ) for the coordinate and velocity distribution function was the central moment of these studies .
two different approaches were used , that were checked to produce compatible results .    in the first of these @xcite by expanding the distribution function in hermite functions in velocity and fourier - series in space coordinate , the fokker - planck equation was transformed into an infinite set of coupled equations for the expansion coefficients .
these equations were solved by a matrix continued fraction method .
particle flux was then expressed through the expansion coefficients , after appropriate series truncation . at smaller friction progressively more terms
must be kept , making this method less efficient .
another approach @xcite is best suited for low friction , when the energy becomes a slow variable ( changing on friction dissipation timescales ) .
the distribution function was rewritten in the energy and configuration space coordinate variables , fpe was solved in these new variables with the help of physically - motivated self - consistent ansatz .
many illuminating results were obtained for the distribution functions at various parameter values , stability domains of locked and running solutions , plots of the mean particle velocity , etc .
diffusion properties however remained not studied .
this matter was partly addressed in other works , based on simulations of the stochastic langevin equation
.    early systematic study of diffusion by means of numerical simulation of the langevin equation in periodic potential with extra constant forcing was undertaken in @xcite ( similar simulations were reported in @xcite for forcing with a time - periodic component ) .
costantini & marchesoni analyzed both overdamped and underdamped settings . at low dissipation
they observed significant enhancement of the spatial diffusion near the critical force @xmath3 ( @xcite , fig .  1 : diffusivity ten times larger than its post - maximum value at @xmath4 , at the temperature equal to 0.6 of the potential well depth @xmath5 ) , and related that to features of the particle jump statistics and locked - to - running transition .
further progress in studies of the diffusion under the action of constant force was made in @xcite .
time - dependence of the particle ensemble dispersion was studied in @xcite .
it was shown for the first time that in underdamped systems a special regime of dispersionless transport was realized , in which dispersion virtually does not change with time , on certain limited interval of time .
the authors explained this phenomenon , they showed that strongly non - equilibrium distribution of particles in space ( with steep front and exponential tail ; formed as the particles first exit their original potential well ) persists for long times , before eventually broadening and assuming normal gaussian shape .
the temperature dependence of the diffusivity was studied as well @xcite .
the authors observed that the maximal diffusivity @xmath6 was achieved near the @xmath3 . at low friction @xmath6
was shown to depend abnormally on the temperature  increase with decreasing temperature @xmath7 .
fitting ( for the three temperatures studied ) yielded @xmath8 relation .
the presence of such abnormal temperature dependence is one important aspect in which underdamped systems differ from overdamped ones .
in fact , certain peculiarities in diffusivity at low temperatures are observed ( for driven brownian particles in space - periodic potentials ) in overdamped situation as well .
p. reimann and colleagues @xcite showed analytically that for a sinusoidal in space potential at external constant force value near @xmath9 , the value required for direct particle pull over the potential barrier ( eq .  [ eq : flat ] below , p.  ) , the ratio of @xmath10 to the diffusivity value @xmath11 in the viscous medium without lattice and bias forces grows at temperature decreasing @xmath12 . yet
the diffusivity itself still vanishes at @xmath13 , as @xmath14 @xcite , and grows monotonically for all @xmath7 .
regions with @xmath15 were obtained in overdamped problem in @xcite ( cf .
fig .  6 ) for specially crafted exotic potential , flat apart from narrow peaks .
the authors linked this abnormal diffusivity behavior to the large ratio of relaxation to escape time in such a system .
similar abnormal @xmath16 behavior was found in @xcite ( cf .
24 ) for piecewise linear potential at extreme asymmetry , essentially sawtooth .
so , while possible at large friction , special atypical conditions must be met for existence of parameters @xmath17 at which @xmath15 . on the contrary ,
such parameters exist universally in the problem with small friction ( @xcite , more below ) .
another distinction of the diffusivity in overdamped problem is that the maximum in @xmath18 is achieved close to the @xmath9 . whereas at low friction @xmath19 is achieved at ( typically much smaller ) @xmath20 .
experimental data for diffusion of colloidal particles in a tilted potential created with laser traps @xcite , fig .  1 , agreed reasonably well with this analytic @xcite @xmath18 dependence . the potential barrier height and the free diffusion coefficient were treated as fit parameters .
similarly , diffusivity in @xcite ( fig .  7 ) showed peaked at @xmath9 behavior , and agreed fairly well with analytic predictions from @xcite with the potential form inferred from independent study of the colloidal particle spatial distribution function . in this work colloidal particles diffused over the potential created by the bottom grid of colloidal spheres , and extra constant forcing was due to gravity ( controlled by the sample inclination angle ) .
different limiting cases of analytic predictions @xcite for the mean particle velocity and diffusivity were compared with experimental results at different settings , and showed good agreement .
the physical reason behind the abnormal temperature dependence of diffusivity in underdamped systems is traced to increasing jump length of running particles ( before getting retrapped at another potential minimum ) with the temperature decrease @xcite .
this phenomenon being absent in overdamped systems ( as these have no proper running states ) .
the maximal diffusivity is achieved at forcing near @xmath20 , as at this value populations of locked and running particles are about equal , and it is the mutual motion between these two populations ( with temperature - independent running particle speed ) that leads to fast spreading of the particle packet , manifested as giant diffusion @xcite .
more thorough analysis of ( the temperature- and force - dependence of ) the diffusion in underdamped systems was carried out in @xcite . it was shown by numeric simulations of the langevin equation that the phenomenon of diffusivity growth with the temperature decreasing ( in the current work called tad , `` temperature - abnormal diffusivity '' ) was only manifested in a narrow interval of applied external forcing .
study in a wide range of temperatures demonstrated that the diffusivity grew with decreasing temperature as @xmath21}$ ] for certain @xmath22 , checked valid for @xmath23 $ ] ; whereas the power - law fit @xmath8 proposed in @xcite could approximate the results only in the narrow temperature range studied by the authors .
we showed that the main exponential growth resulted from such of the correlation time at the temperature decreasing ( found from the simulation results ) .
the diffusion properties were further elucidated and systematized for a broad range of the forces @xmath24 , friction coefficients @xmath25 and temperatures @xmath7 in @xcite , in sinusoidal 1d potential in underdamped setting .
we introduced a velocity potential @xmath26 , which in part governed time evolution of certain slow ( not changing much on times of order one period @xmath27 of small oscillations in the potential , rather evolving on viscous dissipation timescales ) velocity variable , defined through a discrete quantity , maximal velocity in such consecutive oscillations .
with certain reservations that @xmath26 also defined the velocity distribution function @xmath28 agreeing well with the one found numerically .
the latter showed clear bimodal structure with maxima at @xmath29 ( locked states ) and @xmath30 ( running ones ) at intermediate forcing values . in this region of forcing values tad
was observed .
analytic expressions for the mean particle velocity @xmath31 and diffusivity were obtained with the model @xmath28 , good agreement was shown with direct langevin equation simulation results .
simple fits were proposed determining the model parameters [ entering @xmath26 and @xmath28 ] through @xmath24 and @xmath25 , which  apart from usefulness for experimentalists / applications  allowed one to comprehend the features of transport and diffusivity in all the physical parameter space .
physical reason for the @xmath33}$ ] scaling was explained . more detail in the introductory part of sec .
[ s : results ] , p.   below .      in practical applications periodic in time external forces ( electromagnetic , acoustic , etc .
) are more straightforwardly realized than constant forcing .
it is natural to study how the findings for the constant forcing problem are modified for the time - periodic forcing setup , both from theoretical and practical point of view . by changing the diffusivity with the aid of periodic external fields it is possible to create surface structures with prescribed properties ; influence the dynamics of point defects , under irradiation in part ; alter the dislocation motion dynamics ; increase chemical reaction rates in selected directions ; manipulate the efficiency of processes in biological systems , to mention just a few applications @xcite .
most of the works on periodic forcing deal with ratchet - type systems @xcite or stochastic resonance .
studies focusing on diffusion properties are scarce .
we are not aware of any works that systematically study the combined dependence of the diffusivity on the temperature , the forcing amplitude and frequency in a representatively broad range of these parameters .
dependence of the diffusivity on the interplay between the external periodic ( piecewise constant , eq .
[ f : pcp ] on p.   below ) forcing and the lattice force was studied in @xcite in 2d overdamped setup .
the forcing frequency was fixed at one value .
it was shown that in a system with a sum of yukawa potentials at square grid nodes the particle diffusivity was significantly enhanced for certain intricate set of bands ( `` a - windows '' ) of the forcing amplitudes and orientations with respect to the grid . at special forcing arrangement tad
was observed for some diffusivity tensor components , and one of its principal values increased exponentially with inverse temperature .
the other diffusivity tensor principal value decreased exponentially at the same time .
it is likely that multi - dimensionality effects were crucial for tad . instead of having two , locked and running , populations in 1d accounting for ( already sophisticated in this case ) diffusion features , flows in higher dimensional problems are far more complex , known to result in further peculiarities in transport and diffusion @xcite .
underdamped 1d system was studied in @xcite , with a sinusoidal potential ( same as eq .  [ ux ] below ) and
sinusoidal in time forcing .
diffusivity enhancement by orders of magnitude was seen .
tad was observed for some ( low ) forcing frequencies , at the forcing amplitude @xmath34 studied .
dependence of the diffusivity on the forcing frequency was studied in more detail in @xcite , and was shown to be nonmonotonic at large forcing amplitudes ( @xmath34 and @xmath35 , at @xmath36  k , @xmath37 . ) at @xmath34 the maximum in the frequency dependence of the diffusivity @xmath38 persisted for all the temperatures studied , corresponding frequency @xmath39 increasing from @xmath40 to @xmath41 ( in units of the potential eigenfrequency ) at @xmath7 varying from 120 to 360  k , @xmath42 $ ] .
the goal of the present work is to investigate in detail the dependence of the diffusivity on the amplitude and the frequency of uniform time - periodic external forcing in a 1d underdamped space - periodic system , and conditions under which tad is realized .
we mainly study piecewise constant periodic ( pcp : eq .
[ f : pcp ] ) forcing as this case is more straightforward to relate to the results for constant forcing .    the equations solved and the methods used are described in sec .  [ s : setup ] .
[ s : results ] we show how the diffusivity changes with the pcp forcing frequency @xmath0 , starting from the understood results at constant forcing , @xmath43 .
we observe that tad is realized for some ( intermediate ) forcing amplitudes .
[ s : dvst ] is devoted to studying the @xmath16 behavior at fixed @xmath0 .
the results for the diffusivity in the case of sinusoidal in time forcing are shown in sec .
[ s : fsin ] .
we conclude in sec .
[ s : concl ] .
in this paper the 1d dynamics of a particle subject to external forcing @xmath44 is described by the langevin equation : @xmath45 where @xmath46 is the time , @xmath47 is the particle coordinate , @xmath48 is its mass , @xmath25 is the friction coefficient .
overdot stands for time differentiation .
@xmath49 represents thermal fluctuations , described by gaussian white noise with correlation @xmath50 here @xmath51 is the boltzmann constant , @xmath7 is the temperature .
the potential energy @xmath52 is taken @xmath53 where @xmath54 is the lattice constant and @xmath5 is the potential barrier height .
the lattice force @xmath55 acting upon the particle is @xmath56 quantity @xmath57 corresponds to the minimal external force needed to drag the particle over the potential barrier separating potential minima on the 1d lattice at large friction .
we use the same physical parameters as in @xcite , typical for adatom diffusion on close - packed metal surfaces .
namely , values @xmath58 for the activation barrier and the lattice constant are adopted .
two possibilities for the external periodic forcing are studied : piecewise constant periodic ( * pcp * ) driving @xmath59\ ] ] and sinusoidal one , @xmath60 @xmath0 is the angular frequency of the external force and @xmath61 is its amplitude .
it is convenient to use dimensionless time @xmath62 and space coordinate @xmath63 , normalizing the physical values to the period of small oscillations @xmath64 at the potential minimum in the absence of friction , and to the lattice spatial period @xmath54 respectively : @xmath65    dimensionless temperature @xmath66 and friction coefficient @xmath67 are further introduced : @xmath68 definitions ( [ dimtx],[dimtgamma ] ) differ from such in @xcite by extra factors of 2 appearing in the definition of @xmath62 , and consequently in @xmath66 and @xmath67 as well .
current definition makes interpreting the results easier ( @xmath69 corresponds to one oscillation period at small @xmath67 ) . for more straightforward comparison with our previous works we also use another dimensionless temperature in the figures , @xmath70 the parameters of the dimensionless problem are @xmath71 , @xmath72 , @xmath67 , @xmath66 . in the new variables
the equation of motion has form eq .  1 with unit new mass , potential @xmath73 with no prefactors and
the external forcing amplitude made dimensionless with a @xmath74 multiplier .
below lowercase @xmath75 s ( frequencies ) , and corresponding periods @xmath76 denote dimensionless quantities . physical frequencies and periods are obtained by respectively dividing and multiplying these by the the period @xmath77 of small oscillations in the lattice potential at @xmath78 .
overdamped problem is characterized by @xmath79 , whereas @xmath80 ( equivalently , @xmath81 ) corresponds to underdamped dynamics .
we study the underdamped case here .
friction coefficient @xmath82 we use corresponds to the same physical value as in earlier studies @xcite .
we solve stochastic equation ( [ langevin],[xi_cor ] ) numerically using a verlet - type algorithm @xcite with a time step @xmath83 somewhat shorter than @xmath84 of the period of oscillation @xmath77 .
the statistic averaging is performed over the ensemble consisting of at least @xmath85 particles . to verify the modeling consistency
some computations were performed with @xmath86 .
the initial conditions are set as follows .
each particle is placed at @xmath87 , its velocity chosen at random , with maxwellian distribution for a given @xmath7 .
thermalization over 100 oscillation periods , in the potential @xmath52 but without the @xmath88 , is then used to get equilibrium particle distribution over both coordinate and velocity .
tests show that the distribution function does not change after that time . in the process of such thermalization
the particles can jump to nearby elementary cells of the lattice .
such relocated particles are returned back to the initial cell by translation over an integer multiple of the lattice constant , to get diffusion of all particles starting from the initial cell of the lattice .
the diffusivity is computed based on the dispersion @xmath89 of the particle distribution in the limit of large times : @xmath90 angle brackets @xmath91 with no subscripts mean averaging over the ensemble . @xmath92 and @xmath93 are dimensionless in our notation .
the physical diffusivity is @xmath94 .    for each diffusivity calculation
we find time @xmath95 , after which the dispersion grows linearly with time ( if averaged over the forcing period ) .
the @xmath92 is calculated as @xmath96 at @xmath97 .
it is known that in systems with low dissipation transient regimes of anomalous diffusion can be realized @xcite , characterized by @xmath98 the process with @xmath99 is referred to as superdiffusion , whereas subdiffusion is characterized by @xmath100 .
stages of anomalous diffusion were observed in our simulations . the choice of initial conditions above aided in having these stages in clear form .
we determine the exponents @xmath101 by the least - square linear fitting between @xmath102 and @xmath103 on the time - interval of interest .
in @xcite it was shown that depending on the value of the external force ( constant in that study ) different functional dependences of the diffusivity on the temperature were realized .
plots of the diffusivity as a function of the external force at several temperatures ( based on @xcite , recomputed for the current units ) are shown in fig .
[ f : difconstf ] .
there are three regions of @xmath61 values ( i  iii ) , with qualitatively different dependence of the diffusivity @xmath104 on @xmath7 . in region
i diffusion is enhanced as external force @xmath61 value increases , whereas in region iii diffusion is inhibited as @xmath61 increases . in regions
i and iii @xmath104 grows with @xmath7 increasing .
contrary to this ( normal temperature behavior ) in region ii @xmath104 grows as @xmath7 decreases .
this is the region of `` temperature - abnormal diffusivity '' ( tad ) .
the _ diffusion process _ at sufficiently late times is normal , @xmath105 , however the _ diffusion coefficient
_ @xmath106 ( abnormally ) increases as the temperature decreases .    .
1 : @xmath107 [ @xmath108k for the potential characteristic of the h - adatom diffusion on metal surface ( [ u0a ] ) ] , 2 : @xmath109 , 3 : @xmath110 , 4 : @xmath111 , 5 : @xmath112 , 6 : @xmath113 , 7 : @xmath114 .
+ bottom : the effective velocity potential in the three force regions .
left , middle and right diagrams correspond to @xmath115 , 0.095 and 0.25 respectively.,title="fig:",scaledwidth=47.0% ] + .
1 : @xmath107 [ @xmath108k for the potential characteristic of the h - adatom diffusion on metal surface ( [ u0a ] ) ] , 2 : @xmath109 , 3 : @xmath110 , 4 : @xmath111 , 5 : @xmath112 , 6 : @xmath113 , 7 : @xmath114 .
+ bottom : the effective velocity potential in the three force regions
. left , middle and right diagrams correspond to @xmath115 , 0.095 and 0.25 respectively.,title="fig:",scaledwidth=12.3% ] @xmath116 . 1 : @xmath107 [ @xmath108k for the potential characteristic of the h - adatom diffusion on metal surface ( [ u0a ] ) ] , 2 : @xmath109 , 3 : @xmath110 , 4 : @xmath111 , 5 : @xmath112 , 6 : @xmath113 , 7 : @xmath114 .
+ bottom : the effective velocity potential in the three force regions .
left , middle and right diagrams correspond to @xmath115 , 0.095 and 0.25 respectively.,title="fig:",scaledwidth=10.0% ] @xmath116 . 1 : @xmath107 [ @xmath108k for the potential characteristic of the h - adatom diffusion on metal surface ( [ u0a ] ) ] , 2 : @xmath109 , 3 : @xmath110 , 4 : @xmath111 , 5 : @xmath112 , 6 : @xmath113 , 7 : @xmath114 .
+ bottom : the effective velocity potential in the three force regions .
left , middle and right diagrams correspond to @xmath115 , 0.095 and 0.25 respectively.,title="fig:",scaledwidth=10.0% ] @xmath117    in @xcite it was seen that tad is also possible at time - periodic external forcing . at such forcing tad
was only observed in a limited interval of the forcing frequencies .
studies were carried out at @xmath118 , that would correspond to region iii for the same amplitude of constant forcing ( see fig .
[ f : difconstf ] ) , so tad would not be observed at such @xmath61 value were the forcing @xmath88 constant .    to understand physical reasons of the abnormal temperature dependence of the diffusivity in these two different ( time - periodic , and constant ) forcing regimes we consider the periodic forcing regime in more detail below . for piecewise constant periodic ( pcp )
forcing one can use intuition gained in the constant forcing problem for the time intervals in which @xmath44 stays constant , the findings are thus easier to interpret .
most of the results in this work are for the pcp forcing setup .
in @xcite we showed that the transport properties of an ensemble of underdamped particles in a washboard potential can be found by considering a simpler problem of the overdamped motion in effective potential @xmath26 in velocity space .
the problem showed similarities with the problem of active brownian particle motion @xcite . whether the particle belongs to locked [ tending with time to @xmath29 , one minimum of @xmath26 ] or running [ tending to @xmath30 , another minimum of the asymmetric double - well potential @xmath26 ] population in the absence of noise is determined by the particle s initial state ( its velocity , if placed at the potential minimum at @xmath120 ) .
it was shown that region ii of the force values , the region in which tad is observed , is precisely the region where both locked and running solutions are simultaneously possible .
the @xmath26 was also instrumental for proposing a model distribution function @xmath28 , and for defining diffusivity through kubo relation .
the analytic predictions for @xmath104 and @xmath31 based on the model @xmath28 agreed well with direct numerical simulation results for all values of @xmath121 , temperatures in the range @xmath122 $ ] , friction @xmath123 $ ] ( in the definitions of dimensionless quantities adopted here , eqs .
( [ dimtx],[dimtgamma ] ) ; agreement getting progressively better at smaller @xmath67 ) . in the expression for the diffusivity ,
( the analytic prediction for ) the correlation time @xmath124 factor was shown to be responsible for the main features in @xmath125 dependence , numerically found in earlier study @xcite .
namely , @xmath124 showed proper nonmonotonic dependence on @xmath24 , with a maximum near @xmath20 , progressively more pronounced at smaller @xmath7 . @xmath124
increased at temperature decreasing @xmath126 $ ] ( for certain @xmath22 ) ; this increase at low @xmath7 reflected the kramers rate for transitions in the double - well velocity potential @xmath26 , between the wells corresponding to locked and running states .
schematics of the effective velocity potential in the three regions of the applied force are shown in the bottom insets in fig .
[ f : difconstf ] .
the first inset corresponds to @xmath127 , at which only locked solutions exist .
in the second inset @xmath128 , corresponding to the maximal diffusivity ; both types of solutions coexist at this force value .
finally @xmath129 corresponding to region iii is shown in the third inset ; at this @xmath61 only running solutions are realized .
the boundaries between regions i  iii depend on @xmath25 and @xmath48 ( on @xmath67 in the dimensionless problem ) . at smaller @xmath67
these boundaries move to the left , to the lower force values , and the width of region ii decreases . in the current study
we use the same fixed values of @xmath25 and @xmath48 as in @xcite , and investigate the features of @xmath130 at the three values of the external forcing amplitude listed above .
it may be expected that the same picture for different temperature behaviors of the diffusivity , with analogous three regions i  iii of the external force amplitudes , holds for pcp driving , at least at small frequencies @xmath0 . the description with the effective potential is applicable , the potential @xmath26 keeping constant form for each of the two half - periods of the external forcing [ on which @xmath44 stays constant , this @xmath24 value defining corresponding @xmath131 . ] below we show that in fact the situation is more subtle .
similar three regions of the external force amplitudes @xmath61 exist , differing in qualitative @xmath16 dependence . tad is realized in certain range of intermediate @xmath61 , smoothly dependent on @xmath0 . at decreasing temperatures the dependence
@xmath16 starts progressively deviating from such at constant forcing .
for any fixed @xmath0 at temperatures below certain critical @xmath132 , normal temperature dependence of the diffusivity is restored .
we demonstrate these features by analyzing @xmath130 in sequence for the above values of @xmath61 ( @xmath115 , 0.095 , 0.25 ) , that correspond to regions i , ii , iii in the constant external force problem .
we observe transient anomalous diffusion phase , and see how its characteristics relate to the diffusivity in the asymptotic late time regime , when normal diffusion regime sets in .
the properties of this anomalous diffusion stage are different in the three regions of the applied force amplitudes ; this translates into different functional dependence of the diffusivity on the temperature .
here we investigate how the diffusivity changes with ( dimensionless ) frequency @xmath72 of the external forcing at forcing amplitude @xmath127 .
we show that the diffusivity stays nearly constant at small frequencies , up to reciprocal of the superdiffusion stage duration . at frequencies above this
the diffusivity drops according to a power - law with exponent related to the exponent of the superdiffusion .
this power - law dependence holds up to frequencies a few times less than the lattice potential eigenfrequency .
the diffusivity grows with the temperature ( normal temperature behavior ) , related to the superdiffusion stage duration increasing with @xmath7 .
figure  [ f : sigma006 ] shows the growth of the particle dispersion with time . at constant external forcing this growth demonstrates a clear superdiffusion stage , @xmath133 , @xmath134 , that switches to normal diffusion ( same power - law growth but with @xmath135 ) after certain time @xmath136 .
the @xmath136 as well as the superdiffusion exponent @xmath101 depend on @xmath71 and other problem parameters , @xmath67 and @xmath66 . for the values used @xmath137 , @xmath138 .
lowercase @xmath139 s everywhere are dimensionless , normalized to the period of small oscillations at the lattice potential minimum at @xmath140 , @xmath141 @xmath89 and @xmath92 are dimensionless dispersion and diffusivity .
( region i ) for different force frequencies .
@xmath136 denotes the end of the superdiffusion stage at constant forcing .
@xmath76 is the period of the external periodic forcing .
1 : @xmath143 , 2 : @xmath144 , 3 : @xmath145 , 4 : @xmath146 ( constant external force . )
straight dashed and dotted lines show the power - law fitting behavior of the dispersion.,scaledwidth=47.0% ]    curves 13 correspond to three different frequencies .
oscillations in @xmath147 are clearly seen with double forcing frequency , @xmath148 .
these are particularly large at @xmath75 around @xmath149 ( curve 2 ) . at these values ,
interestingly , intervals of @xmath150 decreasing with time are observed .
this effect can be seen in earlier works @xcite . at frequencies below @xmath149
the curves @xmath150 closely follow that at @xmath146 ( constant external force , curve 4 ) . at higher frequencies the superdiffusion phase effectively ends after one period @xmath139 of the external force , normal diffusion sets in ( in a sense that averaged over the forcing period @xmath151 ) .
higher frequencies of the forcing correspondingly result in lower diffusivity , as the stage of @xmath89 growing with time faster than linearly ends earlier than at lower @xmath75 s .
these results suggest the following approximations for the smoothed over the forcing period dispersion @xmath152 , valid not too close to the transition moments between different diffusion regimes : @xmath153 here @xmath154 ( for the approximately boltzmannian initial distribution of the particles in @xmath155 well of the potential ) , @xmath156 , @xmath157 is close to the first potential well exit time .
constants are chosen in such a way that the spline approximation ( [ sigma2fiti ] ) is continuous .
quantities @xmath101 , @xmath158 depend on the problem parameters , @xmath66 , @xmath71 and @xmath67 . at fixed values of these parameters ,
the last approximation yields the dependence of the diffusivity @xmath159 on the forcing period @xmath139 @xmath160 or on its frequency @xmath161 @xmath162 the diffusivity thus should stay about constant at @xmath163 , and at larger forcing frequencies drop with @xmath75 according to the power law , with exponent equal to @xmath164 , determined by the superdiffusion exponent .
figure  [ f : domega006 ] shows simulation results for @xmath165 .
dotted lines show approximations eq .  [ domegafiti ] .
dashed lines show the diffusivity at constant forcing with the same @xmath61 value .
these results confirm the features suggested by eq .
[ domegafiti ] , away from transition frequency @xmath149 and up to @xmath166 .
@xmath167 at any @xmath75 , and decreases with @xmath75 , up to the frequency values a few times smaller than the potential eigenfrequency .     on the forcing frequency at @xmath168 ( triangles ) and @xmath169 ( circles ) .
the dashed line asymptotes show the diffusivity values at constant forcing with the same @xmath115 .
the slanting dotted line shows the power law drop of @xmath92 at intermediate @xmath75 s .
, scaledwidth=47.0% ]    further ramification might have been added in eqs .
( [ sigma2fiti][domegafiti ] ) for the case @xmath170 , changing @xmath157 to @xmath171 ; that would indeed lead to @xmath172 increasing at @xmath173 . with our parameters , however , resulting approximations are an overstretch ; the diffusion gets markedly nonequilibrium at such high @xmath75 s , resulting in the complex @xmath172 behavior at @xmath75 approaching @xmath174 .
more on this in sec .
[ s : highomega ] .    finally , fig .
[ f : sigmat006 ] shows @xmath150 for a range of temperatures at constant forcing . despite the superdiffusion exponent @xmath101
is seen to decrease with the temperature , the duration of the superdiffusion stage @xmath175 grows fast with the temperature . as a result the dispersion is larger for larger temperatures at any given time . at late times
this translates into diffusivity being larger at larger temperatures ; so normal ( increasing ) @xmath16 dependence is realized at force amplitudes in region i.    ) . 1 : @xmath176 , 2 : @xmath168 , 3 : @xmath142 , 4 : @xmath177 .
the short - dashed straight lines show power - law fits for the transient superdiffusion stage , and the late - time normal diffusion regime.,scaledwidth=47.0% ]      here we repeat analysis of the previous section for periodic forcing with amplitude @xmath128 ( the value at which the maximal diffusivity is achieved at constant forcing ) . in the same manner as for the forcing amplitude in region i , transition from initial superdiffusion stage to normal diffusion at late times is observed , at time @xmath136 . at forcing frequency @xmath75 growing
, the diffusivity again decreases only slightly from @xmath179 till @xmath75 reaches @xmath180 .
it decreases fast , according to the same power law @xmath181 , at larger @xmath75 s , up to @xmath182 .
contrary to the case of weaker forcing considered in the previous section , superdiffusion duration @xmath175 now _ decreases _ with the temperature .
this leads to diffusivity growing at temperature decrease ( tad ) , at moderately low frequencies .
this is the main effect we intended to find in the periodic driving setup , extending our previous findings at constant forcing @xcite .
figure  [ f : sigma0095 ] shows growth of the particle dispersion with time .
again the superdiffusion stage is observed , and it ends after one period of pcp forcing or after @xmath136 ( the end time of superdiffusion at constant forcing ) , whichever is smaller .
thus the same approximations eq .  ( [ sigma2fiti][domegafiti ] )
are applicable for the ( averaged over @xmath139 ) dispersion growth @xmath152 and dependence of the diffusivity on the forcing frequency @xmath165 .
simulation results for @xmath165 ( fig .
[ f : domega0095 ] ) again agree with prediction ( eq .  [ domegafiti ] ) .
( region ii ) for different force frequencies .
curves 15 , bottom to top : ( dimensionless ) @xmath183 @xmath184 , @xmath185 , @xmath186 , @xmath187 , 0 ( constant forcing ) .
, scaledwidth=47.0% ]     on the external force frequency @xmath75 , at @xmath168 ( triangles ) and @xmath169 ( circles ) .
the dashed line asymptotes show the diffusivity values for constant forcing with the same @xmath178 .
dotted line shows the power law drop of @xmath92 at higher @xmath75 s .
the dot - dashed horizontal line asymptote ( for high @xmath75 ) is at the diffusivity value in the lattice periodic potential in the absence of extra driving , @xmath188 , for @xmath142.,scaledwidth=47.0% ]    as @xmath172 does not vary much at @xmath163 , tad must be observed at small frequencies , at it is observed at constant forcing ( @xcite : @xmath189 decreases with the temperature when @xmath61 is in region ii ) .
such an abnormal dependence on the temperature is due to longer jumps the particles perform on average at lower temperatures .
this in turn follows from lower probability at lower temperatures to turn the running particle into a locked state : such a transition requires certain threshold energy transfer to such a particle , and at lower temperatures the probability of such a transfer decreases exponentially with inverse temperature . hence the running particle travels at about the stationary speed @xmath190 ( the speed fluctuating around this average due to the lattice potential ) for increasingly longer distances before getting trapped again in some next well of the potential , resulting in larger diffusivity at smaller temperatures .
this mechanism is not applicable in region i , where only locked solutions are realized .    in terms of the dispersion evolution features ,
this phenomenon of the diffusivity growing at the temperature decreasing results from the superdiffusion duration @xmath191 growing at temperature decreasing , as fig .
[ f : sigmat0095 ] demonstrates .
this qualitatively differs from the situation at forcing amplitude in region i ( end of the previous section ) .
we note that @xmath101 increases with the temperature decreasing , same as in region i.    ) . 1 :
@xmath176 , 2 : @xmath168 , 3 : @xmath142 , 4 : @xmath177 .
@xmath178.,scaledwidth=47.0% ]    at higher frequencies , on the other hand , the diffusivity starts increasing with the temperature ( normal behavior ) .
this results from @xmath136 decreasing with the temperature , as can be seen in fig .
[ f : sigmat0095 ] . despite @xmath179 is larger at smaller temperature @xmath192 , as @xmath75 grows the @xmath193 starts decreasing according to the power law earlier [ at @xmath75 s near the smaller @xmath194 than @xmath195 for @xmath196 does .
and the exponent of this power - law decreases with @xmath7 decreasing , @xmath197  therefore at larger frequencies , commensurate with @xmath198 , the diffusivity starts increasing with the temperature , _
tad disappears_.    as it is the case for the forcing amplitude in region i , @xmath172 decreases monotonically with @xmath75 from @xmath199 , up to the frequency values @xmath200 ( at frequencies so high the particle transport becomes substantially non - equilibrium , that affecting the functional dependence @xmath172 , as discussed in sec .
[ s : highomega ] below ) .      for strong forcing ( region iii , in which the diffusivity decreases with @xmath61 in constant forcing setup , and increases with the temperature )
we demonstrate in this subsection that the dependence of dispersion on time shows new features .
after initial superdiffusion stage at times @xmath202 , prolonged dispersionless stage sets in , till @xmath203 . in this stage
the dispersion stays nearly constant .
as at weaker external forcing ( regions i , ii ) , both these stages may be interrupted earlier in the case of high forcing frequency ( after about a period of the forcing ) . at late times
normal diffusion takes place .
dependence @xmath172 is nonmonotonic here , agreeing with findings in @xcite .
@xmath104 reaches a maximum near @xmath204 .
tad is observed at intermediate frequencies .
figure  [ f : sigma025 ] shows the particle dispersion as a function of time . at low frequencies of the forcing ,
superdiffusion is observed till @xmath205 .
dispersionless stage is realized later , at @xmath206 .
physically @xmath136 is the time at which the particle distribution function in velocities assumes its stationary form .
the distribution in space is still strongly non - equilibrium at @xmath136 , it has exponential tail and sharp front @xcite ( in the direction of @xmath24 ) .
it takes till @xmath207 for the distribution function to assume approximately gaussian shape in space ( if averaged over the potential spatial period ) . on interval @xmath208 dispersion does not increase siubstantially .
after @xmath207 normal diffusion sets in , @xmath209 at @xmath210 .
( region iii ) for different force frequencies .
curves 15 : ( dimensionless ) @xmath183 @xmath211 , @xmath184 , @xmath185 , @xmath186 , @xmath187 .
dotted lines show power - law fitting behavior of the dispersion.,scaledwidth=47.0% ]    the whole curve @xmath147 differs insignificantly from @xmath147 in constant forcing setup if the forcing frequency @xmath212 .
this @xmath147 dependence is modified at higher frequencies of the external forcing ( in a similar way to the behavior discussed in sec .
[ s : regioni][s : regionii ] at weaker forcing ) . for @xmath213
the dispersionless stage is interrupted after about one period of external force , after which time approximately linear dependence ( if averaged over times @xmath214 ) sets in , @xmath215 , @xmath92 equal to its asymptotic ( at @xmath216 ) value . at
yet higher forcing frequencies , @xmath217 , dispersionless stage is not realized at all , normal diffusion sets in right after the prematurely interrupted superdiffusion stage .    similarly to the previous sections , described above dependence of @xmath218 on @xmath219
can be schematized by the following continuous piecewise linear approximation ( meaningful not too close to switching time - points @xmath220 ) : @xmath221 here @xmath222 , @xmath156 ; @xmath223 assumed .
as before , @xmath172 can be found from this , showing the same power law decay as in eq .
[ domegafiti ] at @xmath224 , and additional , linearly growing with @xmath75 feature at intermediate frequencies : @xmath225    simulation results for @xmath172 are shown in fig .
[ f : domega025 ] for two temperatures .
they indeed show the features suggested in eq .
[ domegafitiii ] . at high frequencies ,
@xmath224 , the superdiffusion stage switches to normal diffusion earlier than at @xmath136 , at time @xmath226 .
this behavior is the same as for the forcing amplitude in regions i  ii .
thus in the same manner the diffusivity grows as @xmath75 decreases , as the longer interval ( till @xmath227 ) of superdiffusive dispersion growth is used at the forcing period @xmath139 getting longer .
however at smaller frequencies , @xmath228 , it is the dispersionless phase that gets interrupted prematurely ( _ i.e. _ it is not realized in full , till @xmath207 , as it would have at @xmath146 ) , so the earlier such a switch occurs ( _ i.e. _ the shorter the @xmath139 ) , the higher diffusivity is obtained in the late - time normal diffusion stage .
this explains the growth of @xmath172 with @xmath75 at @xmath229 .
the maximum in @xmath172 is thus observed , at @xmath230 .     on the forcing frequency , at @xmath168 ( circles ) and 0.194 ( triangles ) .
the dashed lines show approximately linear growth at @xmath231.,scaledwidth=47.0% ]    this differs from monotonically decreasing @xmath172 in regions i , ii .
such a nonmonotonic dependence @xmath172 was observed in @xcite ; the forcing studied was sinusoidal in time , its amplitude ( indeed ) corresponded to region iii in terminology of this work .    @xmath232 together with @xmath233 increase with the temperature .
the maximal value of the diffusivity ( achieved at @xmath232 ) increases as the temperature decreases , due to the dispersionless stage duration @xmath234 growing fast at temperature decreasing . and
so does @xmath172 at a given @xmath75 near @xmath198 in a limited range of temperatures .
so abnormal temperature dependence of the diffusivity ( tad ) is realized in intermediate region of the forcing frequencies around @xmath198 .
the temperature dependence of the diffusivity is studied in more detail in sec .
[ s : dvst ] below .      as the simulation results suggest , the power law decay of @xmath172 at large @xmath75 s ( figs .
[ f : domega006 ] , [ f : domega0095 ] , [ f : domega025 ] ) does not hold at frequencies that are too high , above @xmath235 of the eigenfrequency at the lattice potential minima .
this is due to the diffusion process becoming markedly non - equilibrium .
the distribution function in velocities differs significantly from such at constant forcing at all times ; description in terms of the velocity potential @xmath26 is no longer useful .    to analyze this
we introduce kinetic temperature @xmath236 ( _ i.e. _ , made dimensionless with the same multiplier @xmath237 as @xmath238 ) , and study its evolution at different @xmath75 s : fig .
[ f : tkin ] , bottom . at low frequency , @xmath239 does not change ( only shows random noise ) after one period @xmath76 of the forcing . at larger @xmath240 regular oscillations in @xmath239
are observed at all times , with amplitude of order 5% of the average @xmath239 value . that average is about 10% higher than the average @xmath239 at @xmath145 .
@xmath239 averaged over the forcing oscillation period grows for about @xmath241 until it reaches its stationary value .
@xmath242 and @xmath185 .
@xmath115 . note that time normalized to one _ forcing _
period @xmath139 ( rather than to @xmath77 as in the previous @xmath46-dependence graphs ) is shown along the abscissa axis.,title="fig:",scaledwidth=24.1354% ] @xmath242 and @xmath185 .
note that time normalized to one _ forcing _
period @xmath139 ( rather than to @xmath77 as in the previous @xmath46-dependence graphs ) is shown along the abscissa axis.,title="fig:",scaledwidth=24.0% ]    together with @xmath239 , @xmath147 also shows transient behavior at @xmath243 .
we illustrate this in the top graphs in fig .
[ f : tkin ] , by plotting @xmath244 in settled `` normal diffusion '' regime this quantity must coincide ( in average over the forcing period of duration @xmath139 ) with the diffusivity , as in such a regime theoretically @xmath245 for @xmath246 .
these results in part mean that the rule of superdiffusion or dispersionless stage ( however these are modified in the high - frequency regime ) terminating and turning into normal diffusion after one period of forcing no longer applies at high forcing frequencies .
this leads to the @xmath172 derivation based on approximations eqs .
( [ sigma2fiti],[sigma2fitiii ] ) progressively less accurate at higher @xmath75 s , the power - law decay first slowing down , and then going into a completely different regime at @xmath75 approaching @xmath174 .    at @xmath247
the diffusivity tends to the diffusivity of particles in the lattice periodic potential in the absence of extra driving @xmath88 .
such an asymptote is shown by the horizontal dot - dashed line in fig .
[ f : domega0095 ] .
having studied how the temperature behavior of diffusion changes at gradual forcing frequency increase , here we investigate in more detail how the diffusivity changes with the temperature at fixed frequency . besides theoretical importance , the results are of interest to experimentalists .
oftentimes changing the frequency in the system is problematic , while changing the temperatures is simpler ; so the findings in this section are more straightforward to test .
@xmath248 at + ( left ) : @xmath145 , force amplitudes @xmath249 0 ( squares , and linear fit ) , 0.06 ( circles ) , 0.095 ( diamonds ) , 0.25 ( triangles ) ; + ( right ) : @xmath178 , top to bottom : @xmath183 0 ( squares ) , @xmath186 ( triangles ) , @xmath185 ( diamonds ) , @xmath250 ( circles ) , @xmath184 ( pluses).,title="fig:",scaledwidth=23.8% ]   at + ( left ) : @xmath145 , force amplitudes @xmath249 0 ( squares , and linear fit ) , 0.06 ( circles ) , 0.095 ( diamonds ) , 0.25 ( triangles ) ; + ( right ) : @xmath178 , top to bottom : @xmath183 0 ( squares ) , @xmath186 ( triangles ) , @xmath185 ( diamonds ) , @xmath250 ( circles ) , @xmath184 ( pluses).,title="fig:",scaledwidth=24.0% ] @xmath251    @xmath252 @xmath117
@xmath117 @xmath253    left part of fig .
[ f : dtomega0001f ] shows how the diffusivity changes with @xmath238 at fixed @xmath254 . for the forcing amplitude in region
i the diffusivity monotonically grows with the temperature , this ( normal ) behavior agreeing with such at constant forcing @xcite . for @xmath61 in regions ii and iii tad is observed , as already noted in sec .
[ s : regionii][s : regioniii ] .
tad is confined to limited ( from below ) temperature intervals , contrary to the tad behavior observed in @xcite for constant forcing . in @xcite
the diffusivity was numerically checked growing ( only in region ii at @xmath146 ) at temperature decreasing down to @xmath176 .
such a behavior was understood theoretically , dependence @xmath255 $ ] was predicted at low temperatures going all the way down to 0 , and agreed well with simulation results .    in the current setup with periodic forcing ,
similar temperature dependence is observed too , @xmath256 but only above certain temperature .
the growth slows down as the temperature keeps decreasing , @xmath104 reaches maximum @xmath257 at certain @xmath258 , and starts decreasing with further decrease in @xmath7 . the behavior at lower temperatures is well approximated by @xmath259 in the absence of forcing we observe @xmath260 ( corresponding dimensionful quantity @xmath261 )  classical arrhenius temperature behavior of the diffusivity .
@xmath262 decreases with @xmath61 , whereas @xmath263 increases . @xmath258 and @xmath10 increase with @xmath61 as well .
temperature behavior of the diffusivity is related to such of the superdiffusion stage duration @xmath264 and @xmath157 .
these times may be inferred from simulations of the particle dispersion evolution with time , fig .
[ f : sigma3f ] . for the forcing amplitude in region i , @xmath157 ( the time at which superdiffusion starts , about the potential well escape time ) increases with the temperature decreasing .
@xmath175 at the same time decreases .
this is seen to result in the diffusivity growing with the temperature , agreeing with the results shown in the @xmath145 part of fig .
[ f : dtomega0001f ] .    in region
ii @xmath157 increases with the temperature decreasing as well .
@xmath264 however also increases with the temperature decreasing ( middle graph in fig .
[ f : sigma3f ] , curves 35 )  contrary to the behavior for @xmath61 in region i. tad is realized at these temperatures . at temperature @xmath265 ( @xmath266 in this graph ) @xmath136 becomes equal to @xmath76 , the period of forcing . at this temperature
the maximal diffusivity is achieved , corresponding to the longest actual superdiffusion stage at the given forcing frequency .
tad is only realized at temperatures above @xmath267 : even though @xmath264 ( superdiffusion duration at constant forcing ) still grows at further temperature decrease together with @xmath136 , in periodic forcing setup superdiffusion effectively ends earlier , at @xmath227 .
it is this superdiffusion interruption at @xmath227 that makes @xmath16 behavior at lower temperatures to diverge in the periodic forcing setup from that at constant forcing . in the constant forcing problem , based on the same logic , tad
must be realized all the way down to @xmath268 ; in agreement with conclusions in @xcite .
it is clear from this argument that there is no uniform convergence with @xmath75 of @xmath269 curves to @xmath270 . at @xmath75
getting smaller @xmath269 gets close to such a @xmath16 curve for constant forcing on longer @xmath7-interval , down to progressively lower temperatures .
however at any @xmath75 there exists @xmath271 such that normal temperature dependence of the diffusivity is restored at @xmath272 .
this non - uniform convergence is shown in the right plot in fig .
[ f : dtomega0001f ] .
the same behavior is observed in region iii .
the maximal diffusivity is achieved at @xmath267 corresponding to the longest superdiffusion stage at the given @xmath75 .
the lowest temperature at which tad is realized is the one at which @xmath136 reaches the period of external force @xmath139 ( more precisely , the time of superdiffusion interruption at a given @xmath75 ) .
at this temperature @xmath273 crosses zero value .
the interval of temperatures , in which tad is observed , is wider than at forcing in region ii , due to the presence of dispersionless phase ( curve 4 in the right plot in fig .
[ f : sigma3f ] ) .
experimentally it is simpler to observe tad at larger @xmath61 , at which within the period of @xmath88 the superdiffusion stage is used to larger extent ; and the temperature interval of tad is thus wider .
@xmath117 @xmath274 @xmath274 @xmath242 + @xmath242 @xmath275 @xmath275 @xmath242    presented study of the diffusion under external forcing can be summarized by the following schematics , fig .
[ f : dtomega ] .
force amplitudes fall into three regions , i  iii . in region
i the diffusivity grows with the temperature increasing . in regions ii
iii the dependence @xmath16 is nonmonotonic , it shows a maximum at certain @xmath258 . below @xmath258
the temperature dependence of the diffusivity is normal , @xmath276 .
temperature interval of tad exists above @xmath258 , @xmath15 .    in the forcing frequency dependence of the diffusivity , on the other hand , it is regions i and ii that look alike : at the frequency decreasing from a fraction of the potential eigenfrequency the diffusivity grows monotonically , getting to a flat region at @xmath163 , in which it virtually reaches its limiting @xmath179 value . in region
iii the @xmath172 curve passes through a maximum , near @xmath204 , corresponding to the end time of superdiffusion stage @xmath136 at constant forcing . at frequencies getting smaller , the forcing period contains progressively longer interval of `` would - be '' ( in the constant forcing problem ) dispersionless and ( at yet smaller frequencies ) normal - diffusion time - interval , resulting in @xmath104 decreasing at such further @xmath75 decrease .
@xmath277 @xmath242 @xmath242 @xmath277
here we show the frequency dependence of the diffusivity for sinusoidal in time external forcing .
we observe the same features in @xmath172 as we saw for pcp forcing .
main conclusions of the sections above thus appear not to depend critically on specific form of time - dependence of the ( smooth enough , symmetric ) forcing .
figure  [ f : dfsin ] presents a comparison between diffusivities at sinusoidal in time and pcp forcing . at small frequencies the diffusivity stays close to its value at @xmath146 . in regions
i and ii of the external forcing amplitude @xmath172 starts decaying according to a power law at @xmath278 , @xmath136 being the superdiffusion stage ending time at constant external forcing with the same amplitude @xmath61 . in region
iii @xmath172 grows linearly with @xmath75 at intermediate @xmath75 s , reaches a maximum , then shows a power - law decay at larger frequencies . for all force
amplitudes the power law decay switches to more complex behavior at @xmath75 reaching about @xmath279 of the potential eigenfrequency .
presence of periodic forcing of any form of time - dependence leads to significant enhancement in diffusivity .
this enhancement is somewhat smaller for sinusoidal forcing in regions i , ii ( can be understood , as at the same @xmath61 the average over period forcing @xmath280 is smaller for sinusoidal than for pcp forcing ) . in region iii the maximum in diffusivity is achieved at smaller @xmath232 , and the value of @xmath10 is larger than such a maximal diffusivity at corresponding @xmath281 under pcp driving with the same amplitude @xmath61 .
this can be understood from the same argument , as in region iii the diffusivity _ decreases _ with the forcing increasing , all other parameters kept constant .
we investigated enhancement of diffusion in 1d space - periodic underdamped systems by external time - periodic fields .
we showed that the diffusion can be enhanced by orders of magnitude at certain choice of temperatures @xmath7 , forcing amplitudes @xmath61 and frequencies @xmath75 .
three regions of @xmath61 exist , i
iii , in which dependences of diffusivity @xmath104 on @xmath7 and @xmath75 differ qualitatively .
at all @xmath61
there is an interval of small frequencies , in which @xmath104 depends only weakly on the frequency , approaching the limiting value @xmath179 .
this value coincides with @xmath104 at constant bias force with the same value @xmath61 in the case of piecewise constant in time periodic ( pcp ) forcing , the form @xmath88 of forcing we focused on in this study . at larger frequencies , from @xmath149 corresponding to the superdiffusion stage ending time @xmath136
, to about @xmath282 ( dimensionful frequency @xmath283 of the potential eigenfrequency @xmath284 ) @xmath285 decays with @xmath0 according to a power - law .
the exponent of this decay is related to the superdiffusion stage exponent @xmath101 . at yet larger frequencies non - equilibrium effects slow down the power - law decay of @xmath285 , and , after a local maximum in @xmath285 near @xmath286 , @xmath104 tends to its asymptotic value coinciding with the diffusivity in the absence of external forcing .    at small forcing amplitudes , in regions i and ii , @xmath285 decreases monotonically with frequency @xmath0 from its value at constant forcing , all the way till the non - equilibrium effects become critical , at @xmath287 . at stronger
forcing , region iii , @xmath172 is nonmonotonic , a maximum is reached at @xmath288 .
@xmath172 grows monotonically on @xmath289 , this growth is approximately linear on certain interval to the left of @xmath198 .
we studied the temperature dependence of the diffusivity at fixed @xmath75 .
contrary to the constant forcing problem @xcite the diffusivity increases with the temperature for all @xmath61 when the temperature is below certain frequency - dependent threshold value .
limited temperature intervals exist for the force amplitudes in regions ii and iii , in which the diffusivity decreases with the temperature ( `` temperature - abnormal diffusivity '' , tad ) .
the physical reason behind the strong diffusivity enhancement is emergence of two populations of particles under the action of external force , locked and running ones . at optimal @xmath61 and
@xmath7 the number of particles in the two populations is comparable , long flights of the running particles relative to the locked population take place with significant probability , resulting in giant enhancement of diffusion .    for comparison we studied diffusivity under sinusoidal in time driving .
we saw that the same features in @xmath172 are observed in the three regions of force amplitudes as in pcp driving case .
we thus conjecture that the qualitative features of @xmath290 behavior are insensitive to specific functional dependence of the ( symmetric , smooth enough in @xmath46 ) external force on time .
the effects investigated allow for simple experimental verification .
diffusion of particles on solid body surface is one natural arena for this .
such systems are characterized by low dissipation @xcite .
another field to study tad in is in propagation of magnetic particles on non - magnetic substrate , acted upon with electromagnetic fields .
abnormal diffusion enhancement could be manifested in , _
e.g. _ , enhanced growth of islands of the new phase at decreasing temperatures .
once verified experimentally , the effect of abnormal diffusion enhancement can find applications in a number of new technologies : in sorting of particles , manufacturing surface structures with required properties , controlling penetration of particles through biological and artificial membranes , in memristors , devices with charge density waves , etc .
50ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]
 + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop _ _ ( ,  ,  ) @noop _ _ ( ,  ,  ) , ,  and  ,  eds . ,  @noop _ _ ,  , vol .
( ,  ,  ) ,  ed .
, link:\doibase 10.1007/978 - 1 - 4419 - 1543 - 6 [ _ _ ]  ( , ,  ) link:\doibase 10.1016/s0370 - 1573(98)00029 - 5 [ * * ,   ( ) ] link:\doibase 10.1039/c0sm01051k [ * * ,   ( ) ] link:\doibase 10.3762/bjnano.7.31 [ * * ,   ( ) ] , link:\doibase 10.1088/0034 - 4885/76/4/046602 [ * * ,   ( ) ] ,   link:\doibase 10.1103/revmodphys.81.387 [ * * ,   ( ) ] , link:\doibase
10.1140/epjst / e2013 - 02071 - 2 [ * * ,   ( ) ] ,   @noop _ _  ( ,  ,  ) link:\doibase 10.1103/physrevlett.96.190601 [ * * ,   ( ) ] ,   link:\doibase 10.1103/physrevlett.105.230602 [ * * , ( ) ] ,   link:\doibase 10.1103/physrevlett.104.248001 [ * * ,   ( ) ] link:\doibase 10.1002/adma.201203500 [ * * ,   ( ) ] link:\doibase 10.1007/s10404 - 013 - 1291 - 9 [ * * ,   ( ) ] link:\doibase 10.1007/bf01325627 [ * * ,   ( ) ] link:\doibase 10.1007/bf01323506 [ * * ,   ( ) ] link:\doibase 10.1016/0375 - 9601(79)90384 - 0 [ * * ,   ( ) ] link:\doibase 10.1007/bf01321244 [ * * ,   ( ) ] link:\doibase 10.1007/bf01352745 [ * * ,   ( ) ] link:\doibase 10.1007/bf01485833 [ * * ,   ( ) ] link:\doibase 10.1209/epl / i1999 - 00510 - 7 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.76.4874 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.82.2820 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.98.020602 [ * * ,   ( ) ] ,   link:\doibase 10.1103/physrevlett.92.250601 [ * * ,   ( ) ] ,   link:\doibase 10.1103/physreve.80.021123 [ * * , ( ) ] http://stacks.iop.org/1367-2630/7/i=1/a=029 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.87.010602 [ * * , ( ) ] link:\doibase 10.1103/physreve.65.031104 [ * * ,   ( ) ] @noop * * , ( ) link:\doibase 10.1016/j.physa.2004.04.019 [ * * ,   ( ) ] link:\doibase 10.1140/epjb / e2013 - 40866 - 7 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.77.041107 [ * * ,   ( ) ] link:\doibase 10.1039/c4sm02387k [ * * ,   ( ) ] link:\doibase 10.1209/0295 - 5075/100/50005 [ * * ,   ( ) ] , ,  and  ,  eds . ,
@noop _ _ , vol .
( ,  ,  ) ,  ed .
, @noop _ _  ( ,  ,  ) @noop _ _  ( ,  ,  )   link:\doibase 10.1209/0295 - 5075/97/60004 [ * * , ( ) ] ,   link:\doibase 10.1140/epjst / e2010 - 01270 - 7 [ * * ,   ( ) ] link:\doibase 10.1134/s0021364012030083 [ * * ,   ( ) ] link:\doibase 10.1088/1742 - 6596/514/1/012045 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.70.051104 [ * * , ( ) ] ,   @noop _ _ ( ,  ,  )   link:\doibase 10.1103/physrevlett.101.190603 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.80.062102 [ * * , ( ) ] ,   link:\doibase 10.1103/revmodphys.79.1381 [ * * ,   ( ) ] link:\doibase 10.1080/00018732.2012.706401 [ * * ,   ( ) ]