Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: nanotechnology , as well as biology , biophysics and chemistry are using or studying setups and objects which are smaller and smaller . in these systems , one is usually interested in mean values , but thermal fluctuations play an important role because their amplitude are often comparable to the mean values . this is for example the case for quantities such as the energy injected in the system or the energy dissipated by the system . these fluctuations can lead to unexpected and undesired effects : for instance , the instantenous energy transfer can be reversed by a large fluctuation , leading energy to flow from a cold source to a hot one . these events , although rare , are quantitatively studied by the recent fluctuations theorems ( fts ) . these theorems give fluctuation relations ( frs ) that quantify the probability of these rare events in systems which can be arbitrarily far from equilibrium . fts have been demonstrated in both deterministic systems  @xcite and stochastic dynamics  @xcite . experiments searching for frs have been performed in dynamical systems  @xcite , but interpretations are very difficult because a quantitative comparison with theoretical prediction is impossible . other experiments have been performed in stochastic systems described by a first order langevin equation : a brownian particle in a moving optical trap  @xcite and an out - of - equilibrium electrical circuit  @xcite in which existing theoretical predictions  @xcite were verified . interesting comments on the langevin equation can be found in @xcite .    in the present article , we study a thermostated harmonic oscillator described by a second order langevin equation . we experimentally search frs for the work done by an external operator and for the heat dissipated by the system , and present analytical derivations of fts based on experimental observations . this paper is organized as follows . in section [ sec : system : description ] , we present the experimental system , write its energy balance to define the work given to the system together with the heat dissipated . we then introduce the fluctuation relations ( frs ) and the fluctuation theorems ( fts ) . in sections [ sec : tftexp ] , [ sec : ssftrampexp ] and [ sec : sinusexp ] , we present experimental results on the fluctuations of first the work and then the heat . a short discussion on experimental results in given in [ sec : expconc ] . then , in sections [ sec : worktheo ] and [ sec : heattheo ] , we present some analytical derivations of fts based on hypothesis inspired by experimental observations . we compare these analytical predictions to the experimental observations and finally conclude in section [ sec : conc ] . our system is a harmonic oscillator and we measure the non - equilibrium fluctuations of its position degree of freedom . the oscillator is damped due to the viscosity of a surrounding fluid that acts as a thermal bath at temperature @xmath0 . our oscillator , depicted in fig . [ fig : pendulum]a , is a torsion pendulum composed of a brass wire ( length @xmath1  @xmath2 , width @xmath3  @xmath2 , thickness @xmath4  @xmath5 ) and a glass mirror glued in the middle of this wire ( length @xmath6  @xmath2 , width @xmath7  @xmath2 , thickness @xmath8  @xmath2 ) . the elastic torsional stiffness of the wire is @xmath9 @xmath10 . it is enclosed in a cell filled by a water - glycerol mixture at @xmath11 concentration . the system is a harmonic oscillator with resonant frequency @xmath12  @xmath13 and a relaxation time @xmath14  @xmath15 . @xmath16 is the total moment of inertia of the displaced masses ( _ i.e. _ the mirror and the mass of displaced fluid ) @xcite . the damping has two contributions : the viscous damping @xmath17 of the surrounding fluid and the viscoelasticity of the brass wire which can be neglected here . the angular displacement of the pendulum @xmath18 is measured by a differential interferometer  @xcite . the measurement noise is two orders of magnitude smaller than thermal fluctuations of the pendulum . @xmath19 is acquired with a resolution of @xmath20 bits at a sampling rate of @xmath21  @xmath13 , which is about 40 times @xmath22 . we drive the system out - of - equilibrium by forcing it with an external torque @xmath23 by means of a small electric current @xmath24 flowing in a coil glued behind the mirror ( fig . [ fig : pendulum]b ) . the coil is inside a static magnetic field . the displacements of the coil and therefore the angular displacements of the mirror are much smaller than the spatial scale of inhomogeneity of the magnetic field . so the torque is proportional to the injected current : @xmath25 ; the slope @xmath26 depends on the geometry of the system . the angular displacement @xmath18 of this harmonic oscillator is very well described by a second order langevin equation : @xmath27 where @xmath28 is the thermal noise , delta - correlated in time of variance @xmath8 and @xmath29 the boltzmann constant and @xmath0 the temperature of the system which is the one of the surrounding fluid . the fluctuation dissipation theorem ( fdt ) gives a relation between the amplitude of the thermal angular fluctuations of the oscillator at equilibrium and its response function . for a harmonic oscillator , the equilibrium thermal fluctuation power spectral density ( psd ) is : @xmath30 where @xmath31 . using fdt ( eq .  [ eq : fdt ] ) , we measure the coefficient @xmath26 and test the calibration accuracy of the apparatus which is better than @xmath32 . more details on the set - up can be found in  @xcite . when the system is driven out of equilibrium using a deterministic torque , it receives some work and a fraction of this energy is dissipated into the heat bath . multiplying eq . ( [ eq : langevin_oscillator ] ) by @xmath33 and integrating between @xmath34 and @xmath35 , we obtain a formulation of the first law of thermodynamics between the two states at time @xmath34 and @xmath35 ( eq .  ( [ eq : energy_conservation ] ) ) . the change in internal energy @xmath36 of the oscillator over a time @xmath37 , starting at a time @xmath34 , is written as : @xmath38 where @xmath39 is the work done on the system over a time @xmath37 : @xmath40 and @xmath41 is the heat given to the system . equivalently , @xmath42 is the heat dissipated by the system . @xmath36 , @xmath39 and @xmath41 are defined as energy in @xmath43 units . the internal energy is the sum of the potential energy and the kinetic energy : @xmath44^{2 } + \frac{1}{2 } c \theta(t)^2 \right\}. \label{eq : udef}\ ] ] the heat transfer @xmath41 is deduced from equation ( [ eq : energy_conservation ] ) ; it has two contributions : @xmath45^{2}{\mathrm{d}}t ' + \frac{1}{k_b \ t } \int_{t_i}^{t_i+\tau } \eta(t ' ) \frac{{\mathrm{d}}\theta}{{\mathrm{d}}t}(t ' ) { \mathrm{d}}t ' . \label{eq : qdef}\end{aligned}\ ] ] the first term corresponds to the opposite of viscous dissipation and is always negative , whereas the second term can be interpreted as the work of the thermal noise which have a fluctuating sign . the second law of thermodynamics imposes @xmath46 to be positive . we rescale the work @xmath39 ( the heat @xmath41 ) by the average work @xmath47 ( the average heat @xmath48 ) and define : @xmath49 ( @xmath50 ) . the brackets are ensemble averages . in the present article , @xmath51 , respectively @xmath52 , stands for either @xmath53 or @xmath54 , respectively @xmath39 or @xmath41 . there are two classes of fts . the _ stationary state fluctuation theorem _ ( ssft ) considers a non - equilibrium steady state . the _ transient fluctuation theorem _ ( tft ) describes transient non - equilibrium states where @xmath37 measures the time since the system left the equilibrium state . a fluctuation relation ( fr ) examines the symmetry of the probability density function ( pdf ) @xmath55 of a quantity @xmath51 around @xmath56 ; @xmath51 is an average value over a time @xmath37 . it compares the probability to have a positive event ( @xmath57 ) versus the probability to have a negative event ( @xmath58 ) . we quantify the fr using a function @xmath59 ( symmetry function ) : @xmath60    the _ transient fluctuation theorem _ ( tft ) states that the symmetry function is linear with @xmath51 for any values of the time integration @xmath37 and the proportionality coefficient is equal to @xmath8 for any value of @xmath37 . @xmath61 contrary to tft , the _ stationary state fluctuation theorem _ ( ssft ) holds only in the limit of infinite time ( @xmath37 ) . @xmath62    the questions we ask are whether fluctuation relations for finite time satisfy the two theorems and what are the finite time corrections . in a first time , we test the correction to the proportionality between the symmetry function @xmath63 and @xmath51 . in the region where the symmetry function is linear with @xmath51 , we define the slope @xmath64 : @xmath65 . in a second time we measure finite time corrections to the value @xmath66 which is the asymptotic value expected by the two theorems . to the applied torque @xmath67 . ]    for the transient fluctuation theorem , we choose the torque @xmath67 depicted in fig .  [ fig : lineartorque]a ) . it is a linear function of time : @xmath68 with @xmath69  @xmath70 and @xmath71  @xmath72 . the value of @xmath73 is chosen such that the mean response of the oscillator is of order of the thermal noise , as can be seen in fig . [ fig : lineartorque]b ) where @xmath19 is plotted during the same time interval of fig . [ fig : lineartorque]a ) . the system is at equilibrium at @xmath74 ( @xmath75 pn.m and @xmath76 pn.m @xmath77 ) . in this section the starting time @xmath34 of integration of all quantities defined before ( @xmath39 , @xmath36 and @xmath41 ) is @xmath74 . so the work is : @xmath78       @xmath79 , @xmath80 @xmath81 and @xmath82 @xmath83 plotted as a function of @xmath37 . b ) pdfs of @xmath53 for various @xmath84 : 0.31 @xmath79 , 1.015 @xmath81 , 2.09 @xmath83 and 4.97 @xmath85 . continuous lines are theoretical predictions with no adjustable parameters . c ) corresponding functions @xmath86 . the straight continuous line is a line with slope @xmath8 . d ) pdfs of @xmath80 for two values of @xmath84 : 4.97 @xmath79 and 8,96 @xmath81 . e ) corresponding pdfs of @xmath54 . continuous lines are gaussian fits . f ) corresponding functions @xmath87 . the straight continuous line is a line with slope @xmath8 . ]    in fig . [ fig : tft]a ) , we represent the time average ( @xmath88 ) of the power injected into the system , the internal energy difference @xmath89 and the time average ( @xmath88 ) of the power dissipated by the system . @xmath88 and @xmath89 are linear in @xmath37 after some short relaxation time @xmath90 defined in the langevin equation : for @xmath84 smaller than @xmath8 , some oscillations around the linear behavior can be seen . the average value of work @xmath47 is therefore quadratic in @xmath37 and is equal to @xmath91  @xmath43 for @xmath92 . the difference between @xmath47 and @xmath93 corresponds to the mean value of dissipated heat @xmath46 ( eq . ( [ eq : qdef ] ) ) . as can be seen in fig.[fig : tft]a ) , @xmath94 is larger than @xmath93 for all times @xmath37 . the average of the dissipated power ( @xmath95 ) is therefore positive for all times @xmath37 as expected from the second principle . for @xmath37 larger than several @xmath90 , the dissipated power is constant and equal to a few @xmath43 per second because @xmath96 and @xmath89 have the same slope after some @xmath90 . so we have the following behavior : the work done by the external work is used by the system to increase its internal energy but a small amount of energy is lost at a constant rate by viscous dissipation and exchange with thermostat . the probability density functions ( pdfs ) @xmath97 of @xmath53 is plotted in fig . [ fig : tft]b ) for different values of @xmath84 . four typical value of @xmath37 are presented : the first is smaller than the relaxation time and the last equals five relaxation times ; the results are the same for any value of @xmath37 . the pdfs of @xmath53 are gaussian for any @xmath37 . we observe that @xmath53 takes negative values as long as @xmath37 is not too large . the probability of having negative values of @xmath53 decreases when @xmath37 is increased . from the pdfs , we compute the symmetry functions . they are plotted in fig .  [ fig : tft]c ) as a function of @xmath53 . they all collapse on the same linear function of @xmath53 for any @xmath37 , which implies they all have the same slope @xmath98 . the straight line in fig  [ fig : tft]c ) has slope @xmath8 . within experimental error bars , @xmath98 is equal to @xmath8 for all time @xmath37 . therefore work fluctuations for a harmonic oscillator under a linear forcing satisfy the tft . we checked that this property is true for other values of @xmath73 and @xmath99 . the pdfs of @xmath80 are plotted in fig . [ fig : tft]d ) for two values of @xmath84 : they are not symmetric and have exponential tails . the pdfs of @xmath54 can be seen in fig . [ fig : tft]e ) for the same values of @xmath84 . they are qualitatively different from the ones of the work . we have plotted in the same figure the gaussian fit of the two pdfs of the dissipated heat . it is clear that the pdfs of @xmath54 are not gaussian . extreme events of @xmath54 are distributed on exponential tails . these tails can be interpreted noticing that @xmath100 and @xmath36 have exponential tails . the variance of the pdfs of @xmath54 is also much larger than the variance of the pdfs of @xmath53 . we plot on fig . [ fig : tft]f ) symmetry functions @xmath101 for the same times @xmath84 . only the behavior of large events can be analyzed here because the variance is much larger than the mean @xmath102 . as it can be seen in fig . [ fig : tft]f , @xmath101 is not proportional to @xmath54 , therefore tfr is not satisfied for finite time . within experimental resolution , @xmath101 is constant for extreme events and equal to @xmath6 . this behavior can be interpreted by writing for large @xmath54 , @xmath103 where @xmath104 and @xmath105 are the decrease rate on the exponential tails . each coefficient depends on @xmath37 . there is a simple expression of @xmath101 for large fluctuations : @xmath106 in fig . [ fig : tft]c ) , the pdfs of @xmath54 are symmetric around the mean value for the two values of @xmath37 . it is not the case for small @xmath84 . thus we can conclude that @xmath107 and that the symmetry function is so equal to the constant : @xmath108 .    as it can be seen in fig . [ fig : tft]e ) , the pdfs become more and more gaussian when @xmath37 tends to infinity . it is expected that for infinite time , the pdf of @xmath54 is a gaussian . thus , tft appears to be satisfied experimentally in the limit of infinite @xmath37 . our interesting finding is that , for @xmath41 tft if not valid for any times . we call a steady state a state in which both forcing and response to the forcing do not depend on the initial time @xmath34 , but only on @xmath37 . this implies that @xmath109 i ndependant of @xmath34 ; and so is @xmath110 . if the torque drifts along time , the mean of @xmath111 is linear with @xmath35 . thus we have to change the definition of the work done on the system to be in a steady state . this is equivalent to a galilean change of reference frame . the work is now defined as : @xmath112 \frac{\textrm d \theta}{\textrm dt}(t ' ) dt ' . \label{eq : work_ssft_def}\ ] ] with this definition , the forcing is @xmath113 and the response to the forcing @xmath114 . when we impose a forcing linear in time ( @xmath115 ) , the first condition ( @xmath116 independent of @xmath34 ) is satisfied . the second ( @xmath110 independent of @xmath34 ) is also satisfied if @xmath117 , _ i.e. _ after a transient state . thus the system is on a steady state . we remark that , in the transient state , this definition of the work reduces to the usual one , because @xmath118  @xmath70 . the average of @xmath39 is quadratic in @xmath37 for any value of @xmath84 . there are no oscillations in time for small @xmath84 . the pdfs of @xmath53 are gaussian for any value of @xmath84 ( fig . [ fig : ssftramp]a ) . probability of negative values is high and decreases with @xmath37 , like in the transient case . the symmetry functions @xmath119 are again proportional to @xmath53 ( fig . [ fig : ssftramp]b ) but the slope @xmath120 is not equal to @xmath8 for smaller @xmath37 and tends to @xmath8 for @xmath121 only , as can be seen in fig . [ fig : ssftramp]c . thus we obtain a fluctuation relation for the work done on the system in this steady state and this relation satisfy the ssft . the slope at finite time is slightly oscillating at a frequency , close to @xmath122 .     for various @xmath84 : 0.019 @xmath79 , 0.31 @xmath81 , 2.09 @xmath83 and 4.97 @xmath85 . b ) corresponding functions @xmath123 . c ) the slope @xmath98 of @xmath123 is plotted versus @xmath37 ( @xmath124 : experimental values ; continuous line : theoretical prediction eq.([eq : epsilon_linear ] ) with no adjustable parameters ) . ]     for various @xmath84 : 0.019 @xmath79 , 0.31 @xmath81 , 2.09 @xmath83 and 4.97 @xmath85 . b ) corresponding functions @xmath123 . c ) the slope @xmath98 of @xmath123 is plotted versus @xmath37 ( @xmath124 : experimental values ; continuous line : theoretical prediction eq.([eq : epsilon_linear ] ) with no adjustable parameters ) . ] the heat dissipated during this linear forcing has a behavior very similar to the one observed in the transient case ( section [ sec : heattft ] ) . we can so transpose here all what we said in section [ sec : heattft ] . we now consider a periodic forcing @xmath125 . this is a very common kind of forcing but it has never been studied in this context . using fourier transform , any periodical forcing can be decomposed in a sum of sinusoidal forcing . we explain here the behavior of a single mode . we choose @xmath126  @xmath70 and @xmath127  @xmath13 . this torque is plotted in fig . [ fig : sinustorque]a . the mean of the response to this torque is sinusoidal , with the same frequency , as can be seen in fig . [ fig : sinustorque]b . we studied other frequencies @xmath128 . the system is clearly in a steady state . we choose the integration time @xmath37 to be a multiple of the period of the driving ( @xmath129 with @xmath130 integer ) . the starting phase @xmath131 is averaged over all possible @xmath34 in order to increase statistics ; in the remaining of this section , we drop the brackets @xmath132 .    . ]       integrated over @xmath130 periods of forcing , with @xmath133 ( @xmath134 ) , @xmath135 ( @xmath124 ) , @xmath136 ( @xmath137 ) and @xmath138 ( @xmath139 ) . b ) the function @xmath140 measured at @xmath141  @xmath13 is plotted as a function of @xmath142 for several @xmath130 : @xmath143 ; @xmath144 @xmath145 ; @xmath146 . for these two plots , continuous lines are theoretical predictions with no adjustable paramaters ( eq .  ( [ eq : ssft_sinus_mean_pdf_work ] ) and eq . ( [ eq : ssft_sinus_var_pdf_work ] ) ) . c ) the slopes @xmath147 , plotted as a function of @xmath130 for two different driving frequencies @xmath128 = 64 hz ( @xmath124 ) and 256 hz ( @xmath134 ) ; continuous lines are theoretical predictions from eq . ( [ eq : epsilon_sinus ] ) with no adjustable parameters . ] integrated over @xmath130 periods of forcing , with @xmath133 ( @xmath134 ) , @xmath135 ( @xmath124 ) , @xmath136 ( @xmath137 ) and @xmath138 ( @xmath139 ) . b ) the function @xmath140 measured at @xmath141  @xmath13 is plotted as a function of @xmath142 for several @xmath130 : @xmath143 ; @xmath144 @xmath145 ; @xmath146 . for these two plots , continuous lines are theoretical predictions with no adjustable paramaters ( eq .  ( [ eq : ssft_sinus_mean_pdf_work ] ) and eq .  ( [ eq : ssft_sinus_var_pdf_work ] ) ) . c ) the slopes @xmath147 , plotted as a function of @xmath130 for two different driving frequencies @xmath128 = 64 hz ( @xmath124 ) and 256 hz ( @xmath134 ) ; continuous lines are theoretical predictions from eq . ( [ eq : epsilon_sinus ] ) with no adjustable parameters . ] the work is written as a function of @xmath130 , the number of periods of the forcing : @xmath148 the pdfs of @xmath142 are plotted in fig . [ fig : sinuswork]a . work fluctuations are gaussian for all values of @xmath130 as in previous cases . thus symmetry functions are again linear in @xmath142 ( fig . [ fig : sinuswork]b ) . the slope @xmath147 is not equal to @xmath8 for all @xmath130 but there is a correction at finite time ( fig . [ fig : sinuswork]c ) . nevertheless , @xmath147 tends to @xmath8 for large @xmath130 , so ssft is satisfied . the convergence is very slow and we have to wait a large number of periods of forcing for the slope to be @xmath8 ( after @xmath149 periods , the slope is still @xmath150 ) . this behavior is independent of the amplitude of the forcing @xmath73 and consequently of the mean value of the work @xmath151 . the system satisfies the ssft for all forcing frequencies @xmath128 but finite time corrections depends on @xmath128 , as can be seen in fig . [ fig : sinuswork]c . ( @xmath134 ) and @xmath152 ( @xmath124 ) . in the next plots , the integration time @xmath37 is a multiple of the period of forcing , @xmath153 , with @xmath133 ( @xmath134 ) , @xmath135 ( @xmath124 ) , @xmath136 ( @xmath137 ) and @xmath138 ( @xmath139 ) . continuous lines are theoretical predictions with no adjustable parameters . b ) pdfs of @xmath36 . c ) pdfs of @xmath54 . d ) symmetry functions @xmath101 . e ) the slope @xmath154 of @xmath101 for @xmath155 , plotted as a function of @xmath130 ( @xmath134 ) . the slope @xmath147 of @xmath119 plotted as a function of @xmath130(@xmath124 ) . continuous line is theoretical prediction . ] ( @xmath134 ) and @xmath152 ( @xmath124 ) . in the next plots , the integration time @xmath37 is a multiple of the period of forcing , @xmath153 , with @xmath133 ( @xmath134 ) , @xmath135 ( @xmath124 ) , @xmath136 ( @xmath137 ) and @xmath138 ( @xmath139 ) . continuous lines are theoretical predictions with no adjustable parameters . b ) pdfs of @xmath36 . c ) pdfs of @xmath54 . d ) symmetry functions @xmath101 . e ) the slope @xmath154 of @xmath101 for @xmath155 , plotted as a function of @xmath130 ( @xmath134 ) . the slope @xmath147 of @xmath119 plotted as a function of @xmath130(@xmath124 ) . continuous line is theoretical prediction . ] we first do some comments on the average values . the average of @xmath36 is obviously vanishing because the time @xmath37 is a multiple of the period of the forcing . @xmath156 and @xmath157 have consequently the same behavior and they are linear in @xmath37 , as can be seen in fig . [ fig : sinusheat]a ) but the pdfs of heat fluctuations @xmath158 have exponential tails ( fig . [ fig : sinusheat]c ) . this can be understood noticing that , from eq . ( [ eq : qdef ] ) , @xmath159 and that @xmath36 has an exponential pdf independent of @xmath130 ( fig . [ fig : sinusheat]b ) . therefore , in a first approximation , the pdf of @xmath54 is a convolution between an exponential distribution ( pdf of @xmath36 ) and a gaussian distribution ( pdf of @xmath53 ) . symmetry functions @xmath160 are plotted in fig . [ fig : sinusheat]d ) for different values of @xmath130 ; three different regions appear :    \(i ) for large fluctuations @xmath158 , @xmath160 equals @xmath6 . when @xmath37 tends to infinity , this region spans from @xmath161 to infinity . \(ii ) for small fluctuations @xmath158 , @xmath160 is a linear function of @xmath158 . we then define @xmath154 as the slope of the function @xmath160 , _ i.e. _ @xmath162 . this slope is plotted in fig . [ fig : sinusheat]e ) where we see that it tends to @xmath8 when @xmath37 is increased . so , ssft holds in this region ii which spans from @xmath163 up to @xmath164 for large @xmath37 . \(iii ) a smooth connection between the two behaviors . we observe that @xmath147 matches experimentally @xmath154 , for all values of @xmath130 ( fig . [ fig : sinusheat]e ) . so the finite time corrections to the ft for the heat are the same than the ones of ft for work : @xmath165 . these regions define the fluctuation relation from the heat dissipated by the oscillator . the limit for large @xmath37 of the symmetry function @xmath101 is rather delicate and we will discuss it in section [ sec : heat_sinus_theo ] . in the previous sections , we have presented experimental results on a harmonic oscillator driven out of equilibrium by an external deterministic forcing @xmath23 . we operated with two different time - prescriptions : one in which @xmath23 is a linear function of time , and one in which @xmath23 is a sinusoidal function of time . the energy injected into the system is the work @xmath166 of the torque @xmath23 . the pdfs of the work @xmath166 are gaussian whatever the time prescription of @xmath23 is , and work fluctuations satisfy a tft ( @xmath23 linear in time ) and a ssft ( @xmath23 linear or sinusoidal in time ) . the energy dissipated by the system is represented by the heat @xmath167 , and we measured it using the first principle of thermodynamics ( eq .  [ eq : qdef ] ) . heat probability distributions are not gaussian and are very different from the ones of the work . they nevertheless satisfy a ssft in both the case of a sinusoidal forcing and a linear forcing . but they do not satisfy a tft in the case of a linear forcing , because the symmetry functions are not linear for all values of dissipated heat @xmath54 .    in the next two sections , we use some experimental evidences to derive analytical expressions of the pdfs of work and the heat exchanged on an arbitrary time interval @xmath37 . we then derive fts together with their finite time corrections . in this section , we derive the analytical expression of the pdf of the work given to the system , and defined as the work of the torque applied to the pendulum , which is either linear or sinusoidal in time . experimentally , we observed that the pdfs are always gaussian , so we restrict our task to deriving expressions for the first two moments of the work distribution . to do so , we use experimental observations on the fluctuations of the angle @xmath18 , as described in section  [ sec : exp_observations ] below . we then compute in section [ sec : pdf_w ] the mean and the variance of the work @xmath168 in the different experimental situations , and then write formally the corresponding fluctuations relations , from which we obtain analytical expressions of the finite time corrections to the fluctuation theorems . we discuss here the angular fluctuations . we decompose the angle @xmath18 into a mean value @xmath169 and a fluctuating part @xmath170 , writing @xmath171 . the mean value corresponds to an ensemble average . it is obtained experimentally by averaging over realisations of the forcing , and it is presented in fig . [ fig : lineartorque ] and [ fig : sinustorque ] . a first experimental observation is as follows . the measured mean response @xmath169 is exactly equal to the solution of the deterministic second order equation obtained when removing the noise term ( @xmath172 ) in the langevin equation  ( [ eq : langevin_oscillator ] ) . we checked this from our data , and found this way a value of the calibration @xmath26 ( see section  [ sec : pendulum ] ) in perfect agreement with the one obtained from the application of the fluctuation dissipation theorem . a second experimental observation concerns the probability distribution of @xmath170 in out - of - equilibrium conditions . we know and observed that at equilibrium , @xmath170 has a gaussian distribution with variance @xmath173 , and the associated momentum @xmath33 has fluctuations @xmath174 which also have a gaussian distribution , with a variance @xmath175 . we observe that the statistical properties of angular fluctuations @xmath170 when a torque @xmath67 linear in time is applied are the same as the statistical properties at equilibrium , when no torque is applied . in figure [ fig : fluct_ang_comparison]a , we plot the pdf of @xmath170 measured at @xmath176 together with the gaussian fit of the pdf at equilibrium ( continuous line ) . the two curves matches perfectly within experimental accuracy . thus we conclude that the external driving does not perturb the equilibrium distribution of angular fluctuations , so we use : @xmath177     when the torque is applied ( @xmath134 ) , compared with a gaussian fit of the pdf at equilibrium ( continuous line ) . b ) the measured spectrum of @xmath170 ( @xmath134 ) is compared with the prediction of fluctuation dissipation theorem in equilibrium ( continuous line ) . ]    the third experimental observation concerns time correlations . in figure [ fig : fluct_ang_comparison]b , we plot the power spectral density function of @xmath170 when applying an external forcing ( @xmath134 ) . we compare it to the prediction of fluctuation dissipation theorem at equilibrium ( eq . [ eq : fdt ] ) computed using the oscillator parameters . the two spectra are identical , so we can confidently use for our system a description in terms of a second order langevin dynamic where the noise term is not perturbed by the presence of the driving . from the power spectral density function of @xmath18 ( eq .  [ eq : fdt ] ) , we derive the autocorrelation function @xmath178 of @xmath170 during a time interval @xmath37 . it is the same at equilibrium and out of equilibrium , and decreases exponentially : @xmath179 where @xmath180 and @xmath181 is defined by @xmath182 and @xmath183 . thus we observe experimentally that when we drive the system out of equilibrum , the angular fluctuations @xmath170 are identical ( with respect to the expressions above ) to those at equilibrium . we verify the same properties for the sinusoidal time prescription of the torque , and use equilibrium expression for the correlation function in the next sections .      in figs . [ fig : tft ] , [ fig : ssftramp ] and [ fig : sinuswork ] , we see that the pdfs of the work are gaussian for any integration time @xmath37 and whatever the forcing is . so these distributions are fully characterized by their mean value @xmath47 and their variance @xmath184 . the external torque @xmath23 is determistic , so the mean value of the work done on the system can be written as : @xmath185 we have defined @xmath186 . the value @xmath187 depends on the time - prescription of the torque we apply to the oscillator . choosing @xmath188 describes the linear ramp and @xmath189 corresponds to the sinusoidal forcing . the variance of the pdfs is : @xmath190 this expression involves the autocorrelation function of the angular speed @xmath174 , ( @xmath191 ) . using the expression of the autocorrelation function of angular fluctuations ( @xmath192 , eq . ( [ eq : autocorrelation ] ) ) , we can calculate exactly the expression of @xmath193 : @xmath194 we have calculated the mean value and the variance of the pdfs in the three situations of interest : stationary and transient cases with a forcing linear in time , and stationary case with a forcing sinusoidal in time . details and results can be found in the appendix . in all of the cases , we compare the theoretical pdfs and the symmetry functions with the experimental results . we have plotted in fig . [ fig : tft ] , [ fig : ssftramp ] and [ fig : sinuswork ] our theoretical pdfs and the corresponding symmetry functions with no adjustable parameters . within experimental error bars , our analytical and experimental results are in excellent agreement . @xmath119 is linear in @xmath53 because the pdfs of the work are gaussian . we now want to calculate analytically the corrections to the slope @xmath98 for finite time @xmath37 . for a gaussian distribution , the symmetry function is : @xmath195 the expression of the slope @xmath98 uses only the mean value and the variance of the gaussian distribution . we define @xmath196 , where the correction @xmath197 is a decreasing function of @xmath37 . we obtain @xmath198 for the transient case , which is in agreement with a tft . for the two steady states , there are corrections to the value @xmath8 ; we find :    a.   linear forcing @xmath199 . \label{eq : epsilon_linear}\ ] ] b.   sinusoidal forcing @xmath200    exact values of the coefficients @xmath201 are given in the appendix . these two expressions are in perfect agreement with experimental results as can be seen in fig . [ fig : ssftramp ] and fig . [ fig : sinuswork ] . these corrections depend on the kind of forcing and it is difficult to predict their form for an arbitrary forcing . nevertheless , the two situations we consider are useful as building blocks of such an arbitrary forcing , and they provide a very nice test of our method . we now determine an analytical expression of the pdf of the dissipated heat . to do so , we make the same hypothesis as in the case of the work ( see section  [ sec : exp_observations ] above ) , and we complete them by additional assumptions to simplify our derivations . we are interested by pdfs of the heat for integration time @xmath37 large compared to @xmath90 , so that exponential corrections which are scaling like @xmath202 can be neglected . in the case of sinusoidal forcing , this is correct after @xmath203 or @xmath204 periods of forcing ( @xmath205 ) . within this assumption , @xmath206 and @xmath207 are independent , and so are @xmath208 and @xmath209 . additionally , as the equation of motion of the oscillator is second order in time , @xmath18 and @xmath210 are independent at any given time @xmath211 . we use the technique proposed in  @xcite . to obtain the pdf @xmath212 of the heat , we define its fourier transform , the characteristic function , as @xmath213 we then write @xmath214 using eq . ( [ eq : qdef ] ) as : @xmath215 where @xmath216 is the joint distribution of the work @xmath168 , @xmath18 and @xmath217 at the beginning and at the end of the time interval @xmath37 . this distribution is expected to be gaussian because @xmath168 is linear in @xmath33 and additionally @xmath18 , @xmath33 and @xmath168 are gaussian . the details of the calculation are given in the appendix . the fourier transform of the pdf of dissipated heat can be exactly calculated : @xmath218 \right.\right . \nonumber \\ & + \left.\left . \frac{-16\cos(\varphi)^2 + 4 + 4 i s(4\cos(\varphi)^2 + 1)}{1+s^2}\right)\right\}. \nonumber \label{eq : fourierpdfheat_linear}\end{aligned}\ ] ] as far as we know , there is no analytic expression for inverse fourier transform of this function , so for the pdf of dissipated heat . however we can do some comments . this expression is very similar to the one found in the case of a brownian particle  @xcite . the factor @xmath219 is the fourier transform of an exponential pdf and this is directly connected to the exponential tails of the pdf . moreover the pdf is not symmetric around its mean , because there is a non vanishing third moment . in this expression , only two terms depend on @xmath37 . for large @xmath37 , this expression reduces to : @xmath220 this expression will turn out to be similar to the one obtain with a sinusoidal forcing , as we will comment in the next section . both expressions depend on the non - dimensional factor @xmath221 defined as : @xmath222 all moments of the distribution of @xmath223 are linear with @xmath224 and @xmath225 . so @xmath224 compares the mean value of the angular speed to the root mean square of the angular speed fluctuations . this coefficient @xmath224 increases when the system is driven further from equilibrium . we consider it as a measure of the distance to equilibrium . in our system @xmath221 is positive , but smaller than @xmath8 , so we are out - of - equilibrium but not very far from it ( @xmath226 ) .      just like in the experiments , we choose the integration time @xmath37 to be a multiple of the period of the forcing , so @xmath227 and therefore @xmath228 . within this framework , we find that the pdf of @xmath229 is exponential : @xmath230 it is independant of @xmath37 because @xmath229 depends only on @xmath18 and @xmath217 at times @xmath34 and @xmath35 which are uncorrelated . this expression is in perfect agreement with the experimental pdfs for all times ( see fig .  [ fig : sinusheat]b ) . some algebra then yields for the characteristic function of @xmath167 : @xmath231 the characteristic function of heat fluctuations is therefore the product of the characteristic function of an exponential distribution ( @xmath232 ) with the one of a gaussian distribution ( @xmath233 ) . thus the pdf of heat fluctuations is nothing but the convolution of a gaussian and an exponential pdf , just as if @xmath168 and @xmath229 were independent . the inverse fourier transform can be computed exactly : @xmath234 \ , , \label{eq : pdfq}\end{aligned}\ ] ] where @xmath235 stands for the complementary erf function . in fig . [ fig : sinusheat]c , we have plotted the analytical pdf from eq . ( [ eq : pdfq ] ) together with the experimental ones , using values of @xmath236 and @xmath237 from the experiment and no adjustable parameters . the agreement is perfect for all values of @xmath130 , _ i.e. _ for any time @xmath37 . from eq . ( [ eq : pdfq ] ) , we isolate three different regions for @xmath87 :    \(i ) if @xmath238 , then @xmath239 . this domain of @xmath240 corresponds to fluctuations larger then three times the average value . the pdf has exponential tails , corresponding to an exponential distribution with a non - vanishing mean . \(ii ) if @xmath241 , then @xmath242 with @xmath243 . in this domain , values of the heat are small and heat fluctuations behave like work fluctuations . the slope @xmath244 is the same as the one found for work fluctuations . the exact correction to the asymptotic value 1 is plotted in fig . [ fig : sinusheat]e and again it describes perfectly the experimental behavior . \(iii ) for @xmath245 , there is an intermediate region connecting domains ( i ) and ( ii ) by a second order polynomial : @xmath246 . these three domains offer a perfect description of the three regions observed experimentally ( fig . [ fig : sinusheat]d ) . now , we examine the limit of infinite @xmath37 in which ssft is supposed to hold . to do so , we distinguish two variables : the heat @xmath223 or the normalized heat @xmath247 . their asymptotic behavior are different because the average heat @xmath237 depends on @xmath37 , more precisely it is linear in @xmath37 .    we discuss first @xmath223 . the asymptotic shape of the pdf of @xmath223 ( eq .  ( [ eq : pdfq ] ) ) for large @xmath37 is a gaussian whose variance is @xmath236 , the variance of the pdf of @xmath168 . thus , the pdf of @xmath223 coincides with the pdf of @xmath168 for @xmath37 strictly infinite . as we have already shown , work fluctuations satisfy the conventional ssft ; therefore heat fluctuations also satisfy the conventional ssft ( eq . ( [ eq : ssft ] ) ) . we have found three different regions separated by two limit values : the mean and three times the mean . but in the limit of large times @xmath37 , the pdf shrinks and only region ( ii ) is relevant . region ( ii ) corresponds to small fluctuations and it is bounded from above by @xmath248 with the average @xmath237 being linear in @xmath37 . so all the behavior of the fluctuations of @xmath223 for large @xmath37 lays in region ( ii ) where the symetry function is linear and ssft holds . we turn now to the normalized heat @xmath247 . as the average value of @xmath223 is linear in @xmath37 , rescaling by @xmath237 is equivalent to a division by @xmath37 ; the mean of @xmath247 is then @xmath8 . this normalization makes the two limit values constant . the boundary between regions ( ii ) and ( iii ) is @xmath249 and the boundary between ( iii ) and ( i ) is @xmath250 . the function @xmath87 is not linear for large values of @xmath251 but it is linear only in region ( ii ) , for @xmath252 , _ i.e. _ for small fluctuations . so ssft is satisfied only for small fluctuations but not for all values of @xmath247 , and we obtain for @xmath247 a fluctuation relation which prescribes a symetry function that is non - linear in @xmath247 . these two different pictures , in terms of @xmath223 or @xmath247 , results from taking two non - commutative limits differently . the first description using @xmath223 implies that the limit @xmath37 infinite is taken before the limit of large @xmath223 . the second description does the opposite . however , the probability to have large fluctuations decreases with @xmath37 and experimentally , for large @xmath37 , only the region ( ii ) can be seen , and it is the region in which where ssft holds .    as we have done in the case of the linear forcing , we introduce a non - dimensional factor @xmath221 such as : @xmath253 the moments of the distribution of @xmath41 are linear with @xmath224 and , like the linear torque , @xmath221 is equal to the amplitude of @xmath33 divided by @xmath254 . we consider it also as a measure of the distance to equilibrium . in our system @xmath221 is positive , but smaller than @xmath8 , so we are out - of - equilibrium but not very far it : here @xmath255 . in this paper , we have studied the fluctuations of energy input and energy dissipation in a harmonic oscillator driven out of equilibrium . this oscillator is very well described by a second order langevin equation . we have performed experiments using a torsion pendulum driven out of equilibrium following a stationary protocol in which either the torque increases linearly in time , or oscillates at a given frequency . we have also studied transient evolutions from the equilibrium state . we have defined the work given to the system as the work of the torque applied during a time @xmath37 . accordingly we have defined the heat dissipated by the pendulum during this time @xmath37 , by writing the first principle of thermodynamics between the two states separated by time @xmath37 . fluctuations relations are obtained experimentally for both the work and the heat , for the stationary and transient evolutions . we have experimentally observed that angle fluctuations of the brownian pendulum have the same statistical and dynamical properties at equilibrium and for any non - equilibrium driving . from this observation , we have derived expressions for the probability density functions of the work and the heat . in our system , fluctuations of the angle are gaussian , and so are fluctuations of the work @xmath53 . so the symmetry functions @xmath119 of the work are linear , and we have calculated exactly the time - correction to the proportionality coefficient between @xmath119 and @xmath53 . these corrections match perfectly the experimental results , both in the case of a forcing linear in time , and sinusoidal in time . we have also computed the analytic expression of the fourier transform of the pdfs of the dissipated heat . for the sinusoidal forcing , we have obtained for the first time an analytic expression of the pdf of the heat . this expression is in excellent agreement with the experimental measurements . for a torque linear in time , the pdf of the heat has no simple expression but its fourier transform gives insight on the behavior of the symmetry function of the heat . it is very similar to the one obtained in the case of a first order langevin dynamics  @xcite . we emphasize here that our analytical derivations are strongly connected to experimental observations on the properties of the noise ; and are therefore different from any previous theoretical approach . we have introduced a dimensionless variable @xmath221 which we think is a measure of the distance from equilibrium : the average dissipation rate is proportional to @xmath221 , and it increases when the system is further from equilibrium . @xmath221 is also proportional to the strength of the driving and in the fluctuation relations , it gives a proper unit to measure the amplitude of fluctuations . so @xmath221 plays the same role as the dissipation coefficient ( the viscosity in our case ) in the fluctuation dissipation theorem at equilibrium . we have an expression of @xmath221 for the two different time - prescriptions we have used . these expressions can be generalized : @xmath256 the numerator corresponds to the solution of the langevin equation when removing the thermal noise term ( @xmath257 ) . the denominator corresponds to the variance of thermal fluctuations of angular speed @xmath258 . we thank g. gallavotti for useful discussions . this work has been partially supported by anr-05-blan-0105 - 01 . in this section , we will calculate the mean and the variance of the work given to the system in the cases :    1 .   transient state , linear forcing 2 .   steady state , linear forcing 3 . steady state , sinusoidal forcing      [ tftramp]the torque is : @xmath259 . the mean value of the angular displacement is the solution of eq . [ eq : langevin_oscillator ] : @xmath260 for the work done on the system , the pdfs are gaussian for all integration time @xmath37 . the mean of the pdf of @xmath168 for a given @xmath37 is : @xmath261 \label{eq : meanworkramp}\end{aligned}\ ] ] and its variance is : @xmath262 , \label{eq : varianceworkramp}\end{aligned}\ ] ]    @xmath263      [ ssftramp]the torque is : @xmath259 . the mean value of the angular displacement is the solution of eq . [ eq : langevin_oscillator ] after some @xmath90 . thus the exponential term is vanished : @xmath264 for the work done on the system , the pdfs are gaussian for all integration time @xmath37 . the mean of the pdf is : @xmath265 and the variance is : @xmath266 . \label{eq : varianceworkrampstat}\end{aligned}\ ] ] from this , we deduce : @xmath267 where @xmath268      [ ssftsinus]the torque is @xmath269 . the mean value of the angular displacement is : @xmath270 where @xmath271 for the work done on the system , the pdfs are gaussian for all integration time @xmath37 . the mean of the pdf is : @xmath272 and the variance is : @xmath273 where @xmath274 @xmath275\,.\label{eq : ssft_sinus_var_pdf_work}\end{aligned}\ ] ] [ appendixheat ] in this section , we will calculate the fourier transform of the pdf of the dissipated heat in two cases :    1 . linear forcing 2 . sinusoidal forcing      we introduce non - dimensional parameters in order to simplify calculations : @xmath276 the mean value and the variance of @xmath277 and @xmath278 can be simply expressed : @xmath279 where @xmath221 is a non - dimensional value . integrating by part , the work @xmath168 can be rewritten : @xmath280-\frac{(d\omega_0)^2}{2}\left[(t_i+\tau)^2-t_i^2\right]+w^{*},\nonumber\\ w^ { * } & = & -(d\omega_0)\int_{t_i}^{t_i+\tau}\tilde{x}(t'){\mathrm{d}}t'.\end{aligned}\ ] ] with these definitions , we obtain : @xmath281 and @xmath282 . like the distribution of @xmath168 , the distribution of @xmath283 is gaussian for all values of @xmath37 and we find : @xmath284 for notational convenience , we introduce a five dimensional vector : @xmath285 . @xmath216 is gaussian and is so fully characterized by the covariance matrix @xmath286 defined as : @xmath287 where @xmath288 denotes the complex conjugate of @xmath289 . so the distribution @xmath216 is written : @xmath290 where @xmath291 denotes the transpose of @xmath289 . we suppose that the integration time is larger than the relaxation time . within this assumption , @xmath206 and @xmath207 are independent , and so are @xmath292 and @xmath293 . as the equation of motion of the oscillator is second order in time , @xmath18 and @xmath33 are independent at any given times @xmath211 . with these hypotheses , we get : @xmath294 the other coefficients of the covariance matrix are : @xmath295    the fourier transform of the pdf of the heat will be calculated here . we define two quantities : @xmath296 one can write @xmath223 as : @xmath297 . the fourier transform can be so written : @xmath298 we use a new variable defined as : @xmath299 with this definition , the argument in the exponential @xmath23 can be rewritten : @xmath300.\end{aligned}\ ] ] changing the integration variable to @xmath301 yields : @xmath302 to make eq . [ eq : gen_heat_lin_exp ] into an explicit expression for @xmath303 , the inverse of matrix @xmath304 is required in the expression for @xmath305 and its determinant . these are obtained as follows . we find : @xmath306 the determinant of the matrix is @xmath307 . for the inverse of this matrix , we get : @xmath308 we now have the material needed to calculate @xmath305 . we find : @xmath309.\end{aligned}\ ] ] the analytic expression of the fourier transform of the pdf of the heat dissipated during a linear forcing is : @xmath310 \right.\right . \nonumber \\ & + \left.\left . \frac{-16\cos(\varphi)^2 + 4 + 4 i \,s(4\cos(\varphi)^2 + 1)}{1+s^2}\right)\right\}.\end{aligned}\ ] ]      [ ap : sinus_heat ] we will determine in a first time the gaussian joint distribution @xmath216 of @xmath168 , @xmath207 , @xmath206 , @xmath293 and @xmath292 . for notational convenience , we introduce a five dimensional vector : @xmath311 . the pdf @xmath216 is fully characterized by the covariance matrix @xmath286 . @xmath312 where @xmath313 denotes the complex conjugate of @xmath289 . we suppose that the integration time is larger than the relaxation time . within this assumption , @xmath206 and @xmath207 are independent , and so are @xmath292 and @xmath293 . as the equation of motion of the oscillator is second order in time , @xmath18 and @xmath33 are independent at any given times @xmath211 . with these hypotheses , we get : @xmath314 the covariance matrix is so a diagonal matrix : @xmath315 @xmath36 is a function of the positions and velocities at the beginning ( @xmath34 ) and at the end ( @xmath35 ) . thus , @xmath36 and @xmath39 can be considered as independent . the pdf of @xmath41 is so the convolution between the pdf of @xmath39 which is gaussian and the pdf of @xmath36 : @xmath316    we first calculate exactly the pdf of the variation of internal energy . we have shown that @xmath207 , @xmath206 , @xmath293 and @xmath292 are independent . the fourier transform of the pdf is : @xmath317 where @xmath318 and @xmath319 . the distribution of @xmath18 is gaussian with variance @xmath320 . the distribution of @xmath321 and the distribution of @xmath322 are the same : @xmath323 the fourier transform of this distribution is : @xmath324 . this distribution is the same for @xmath321 and @xmath322 at @xmath34 and @xmath35 . thus the fourier transform of the variation of internal energy is @xmath325 and the probability is : @xmath326    as @xmath36 and @xmath39 are independent , the fourier transform of the dissipated heat can be calculated : @xmath327 this expression can be inversed because it is simply the convolution between a gaussian distribution and an exponential distribution . so we find eq . ( [ eq : pdfq ] ) . 0        d.  j.  evans and d.  j.  searles , _ phys . e _ * 50 * , 1645 ( 1994 ) . d.  j.  evans , d.  j.  searles , _ advances in physics _ * 51 * ( 7 ) 1529 ( 2002 ) . d.  j.  evans , d.  j.  searles and l.  rondoni , _ phys . e _ * 71 * ( 5 ) 056120 ( 2005 ) . r.  van zon , s.  ciliberto , e.  g.  d.  cohen , _ phys . lett . _ * 92 * ( 13 ) 130601 ( 2004 ) . r.  van zon and e.  g.  d.  cohen , _ phys . * 91 * ( 11 ) 110601 ( 2003 ) ; _ phys . e _ * 67 * 046102 ( 2003 ) ; _ phys . * 69 * 056121 ( 2004 ) . g.  m.  wang , e.  m.  sevick , e.  mittag , d.  j.  searles , d.  j.  evans , _ phys .  rev . lett . _ * 89 * 050601 ( 2002 ) ; g.  m.  wang , j.  c.  reid , d.  m.  carberry , d.  r.  m.  williams , e.  m.  sevick , d.  j.  evans , _ phys . e _ * 71 * 046142 ( 2005 ) .
Abstract:

Answer: we study experimentally the thermal fluctuations of energy input and dissipation in a harmonic oscillator driven out of equilibrium , and search for fluctuation relations . we study transient evolution from the equilibrium state , together with non equilibrium steady states . fluctuations relations are obtained experimentally for both the work and the heat , for the stationary and transient evolutions . a stationary state fluctuation theorem is verified for the two time prescriptions of the torque . but a transient fluctuation theorem is satisfied for the work given to the system but not for the heat dissipated by the system in the case of linear forcing . experimental observations on the statistical and dynamical properties of the fluctuation of the angle , we derive analytical expressions for the probability density function of the work and the heat . we obtain for the first time an analytic expression of the probability density function of the heat . agreement between experiments and our modeling is excellent .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the essential point of stochastic resonance @xcite is the following : if a stochastic resonator is subjected to an external influence , commonly referred to as the signal , its response exhibits signal features , which are most pronounced at a certain level of noise present in the system . the fascinating aspect of stochastic resonance is therefore that , counter - intuitively , an increasing noise level does not steadily deteriorate the transmission of the signal through the resonator . on the contrary , noise may rather be used to optimise the transmission process @xcite . the implications of these effects are intensively investigated . in biological systems , for example , stochastic resonance apparently plays a role in the neural transmission of information @xcite . from a technical point of view one might possibly see the emergence of a novel type of detector which incorporates an optimal amount of noise to perform best . in this regard superconducting quantum interference devices ( squids ) have been studied @xcite . they may be made to detect weak magnetic fields reasonably well without the usual costly shielding from environmental noise .    in recent years it was shown that the performance of a single stochastic resonator can be enhanced , if it is embedded into an ensemble of other stochastic resonators which are properly coupled @xcite@xcite@xcite@xcite@xcite : compared to being operated isolated the response of the resonator to the signal increases within the coupled ensemble . however , if the coupling becomes too strong this response was found to deteriorate again . thus in an ensemble of stochastic resonators the coupling strength turns out to be a second design parameter : apart from the noise level it can be tuned to achieve an optimal performance of the embedded resonator . this effect was called _ array - enhanced stochastic resonance _ @xcite . it will possibly find technical exploitation and might also be relevant to biological systems , for example , coupled neurons . the aim of the present paper is a further investigation of the phenomena related to array - enhanced stochastic resonance . to this end we study a simple prototype system of coupled two - state stochastic resonators under periodic modulation . the model we propose allows to analytically calculate the weak - signal limits of two prominent stochastic - resonance quantities , spectral power amplification ( spa ) and signal - to - noise ratio ( snr ) , respectively . following the concept of array - enhanced stochastic resonance both quantities refer to the response of a resonator as part of the coupled ensemble . the impact of the coupling on both response measures will be studied . there is a close link of the present analysis to the work presented in @xcite and @xcite , which to our knowledge are the most comprehensive in the field . however , compared to @xcite , where the spa of a system of globally coupled bistable elements was studied analytically , too , the present approach offers a more refined description , since it allows to avoid a mean - field approximation . compared to @xcite , on the other hand , where the snr of a chain of locally coupled elements was investigated in a ( rather extensive ) simulation , the snr can now be obtained analytically . for a quick review of the results the reader is referred to the 3rd and 4th passage of the summary section .    from a conceptual point of view the model proposed here can be seen as an extension of the two - state model developed by mcnamara and wiesenfeld @xcite to study stochastic resonance in noisy bistable systems . instead of considering individual two - state elements as in @xcite these elements are now arranged in a chain . a simple next - neighbour interaction is introduced which brings the model close to the system studied in @xcite . the interaction is chosen in a way that the resulting evolution of the elements is given by glauber s stochastic ising model @xcite . a detailed description of this approach is given in the next section which also provides the necessary background . we note that the connection to the glauber model is made for mathematical convenience . we are not concerned with the observation of stochastic resonance in ising systems which was the central theme in @xcite and @xcite . our purpose will also lead to the assumption of arrhenius - type transition rates of uncoupled unmodulated elements which is unusual in the context of ising systems . first results of the present analysis had been published in @xcite .    a further interesting approach to arrays of stochastic resonators is c entered around a response of a more collective nature . here stochastic resonance is studied in the summed output of @xmath0 resonators which need not necessarily be coupled @xcite@xcite@xcite . recently , it was shown analytically that the snr of this output approaches the input snr , if a sufficiently large number of uncoupled resonators is used @xcite . whether in some way coupling may still be beneficial for the collective response of the present model will be the subject of a future study . only the summed output of a very large number of resonators shall be briefly discussed here . it will be shown that in this case the snr always deteriorates under coupling . this indicates that the coupling - induced improvement of the performance of stochastic resonators associated with array - enhanced stochastic resonance is a local rather than a global effect . an early theory of stochastic resonance in noisy bistable systems with one ( generally continuous ) variable was worked out by mcnamara and wiesenfeld @xcite . they studied the effect of a periodic modulation of these systems in a way which is independent of the precise dynamics involved . to this end the behaviour of the bistable system was approximated by a random telegraph process . for simplicity this process was taken to be symmetric , randomly switching between two states @xmath1 . for convenience we shall consider here @xmath2 . the probabilities @xmath3 to find the process in state @xmath4 satisfy @xmath5 and their time - evolution is governed by @xmath6    where @xmath7 denotes the rate of the transition @xmath8 . these rates must be extracted from the precise dynamics at hand . to build a theory of stochastic resonance they were assumed to be approximately given by @xcite    @xmath9\ ] ]    with @xmath10 being the small parameter of the theory . obviously , at this level the periodic modulation of the bistable system with ( angular ) frequency @xmath11 and phase @xmath12 is taken into account by the cosine term , only . this simple model allows to calculate the long - time limit of the power spectrum of the random telegraph process averaged over a uniformly distributed phase @xmath12 . the spectrum consists of a continuous part and a delta function at modulation frequency @xmath11 . the continuous part of the spectrum is the lorentzian of the unperturbed process times a frequency - dependent prefactor . the latter is close to 1 and governs the modulation - induced transfer of broad - band power to the delta peak in the spectrum . from this spectrum mcnamara and wiesenfeld obtained their central theoretical result with respect to stochastic resonance : an exact analytical expression for the signal - to - noise ratio ( snr ) of the two - state process . this snr is defined as the ratio of the weight of the delta function to the continuous part of the spectrum at modulation frequency . in view of the present work we shall neglect the signal - induced suppression of the continuous part of the spectrum and consider the linear response approximation of the mcnamara - wiesenfeld result with respect to the small parameter @xmath10 . doing so one finds for the snr    @xmath13    if the approximation ( [ mwrate ] ) is performed on a particular noisy bistable system subject to a periodic modulation , the parameters @xmath14 and @xmath10 become functions of noise strength and modulation amplitude , respectively . the dependence of the snr on the noise intensity can then be studied and the occurrence of stochastic resonance may be established for that particular system .    as an example the overdamped double - well system was given in @xcite . the model equation reads    @xmath15    where @xmath16 is gaussian white noise with @xmath17 and @xmath18 . using a modified kramers formula for the transition rates valid for sufficiently low modulation frequencies mcnamara and wiesenfeld found    @xmath19 { } \delta = \frac{a}{d}.\ ] ]    the central idea of the present work is to pass from the single bistable element to a set of coupled bistable elements by extending the two - state model of mcnamara and wiesenfeld . to this end we consider a chain of two - state elements which for convenience is taken to be of infinite length . if the elements interact the simple gain - loss balance ( [ mwmaster ] ) has to be modified to describe the evolution of the set of probabilities @xmath20 to find the chain in a particular configuration @xmath21 at time @xmath22 . introducing a formal operator @xmath23 defined for any function @xmath24 by @xmath25 the new gain - loss balance reads    @xmath26    with the initial condition @xmath27 the @xmath20 are interpreted as the transition probabilities @xmath28 of a markovian process with infinitely many discrete components . the statistics of the process is fully determined by ( [ master ] ) .    to introduce interactions we assume that the transition rates ( [ mwrate ] ) of an element depend on the states of its next neighbours . a simple choice for this coupling controlled by a parameter @xmath29 is    @xmath30,\ ] ]    where @xmath31 is the mcnamara - wiesenfeld rate ( [ mwrate ] ) . with positive @xmath29 neighbouring elements prefer to be in the same state , whereas they tend to be in opposite states , if @xmath29 is negative . both tendencies grow with growing coupling strength @xmath32 . to avoid negative transition rates @xmath33 has to be required . we note that the assumed type of coupling is thus able to model ferromagnetic - type interactions for which a coupling - induced improvement of the performance of stochastic resonators was found in @xcite and @xcite . for positive @xmath29 the present model is particularly close to the model simulated in @xcite , since both are linear arrays of next - neighbour - coupled elements . for convenience we shall often use the terms ferromagnetic or antiferromagnetic coupling instead of coupling with positive or negative @xmath29 .    of course , other types of transition rates ( [ transrate ] ) could be considered as well . the advantage of the present choice is that most of the relevant stochastic properties of the resulting model are already known . they were studied by glauber who introduced this model as a stochastic form of the ising model @xcite . ( we note that in the glauber model the term @xmath34 in ( [ transrate ] ) is replaced by a general time - dependent parameter @xmath35 . ) as in the mcnamara - wiesenfeld model it is again possible to find an analytical expression for the snr in leading order of the modulation parameter @xmath10 which is the subject of the next section . the respective analysis is essentially an exploitation of glauber s work .    it may be possible to reduce the dynamics of some coupled systems to the two - state model given by ( [ master ] ) and ( [ transrate ] ) and hence to directly employ the snr formula derived in the next section however , this approach is certainly more involved than the rate expansion ( [ mwrate ] ) needed to make use of the mcnamara - wiesenfeld snr . we shall not pursue this more general aspect of the model in further detail . instead we are interested in devising a simple prototype system which allows to study the impact of the coupling on the stochastic - resonance effect . to this end we are looking for simple assumptions on the dependence of the model parameter @xmath14 , @xmath29 , and @xmath10 on some noise intensity and signal amplitude , respectively . we proceed in two steps :    first , we retain the glauber s original relation of the present model to the ising model given by the hamiltonian    @xmath36    this relation is defined by an assumption on the probability @xmath37 to find the @xmath38-th element in state @xmath39 , if all other elements are fixed : in an adiabatic limit the @xmath37 of both models are to be identical . this implies that    @xmath40    holds which allows to establish a relation between the parameters of both models . one finds @xmath41 and @xmath42 , respectively @xcite . the role of the noise intensity is thus played by the temperature @xmath43 . to skip unnecessary parameters we set @xmath44 and @xmath45 . in the present paper we assume that @xmath46 , which shall be referred to as the signal , is given by    @xmath47    for weak amplitudes @xmath48 this results in @xmath49 . ( in section [ srsnr ] we shall also consider periodic signals ( [ forcing ] ) with an additional noisy background . ) second , a plausible assumption has to be made on the temperature dependence of the time - scale parameter @xmath50 which is not effected by requiring ( [ detailedbal ] ) . it is assumed that it retains the qualitative dependence ( [ mwparam ] ) found by mcnamara and wiesenfeld for the double - well system , i.e. , @xmath51 with properly scaled temperature @xmath43 . this dependence is typical for many rate processes . an appropriate temperature dependence of @xmath14 is also necessary to turn the elements into stochastic resonators for vanishing coupling . for simplicity @xmath52 is considered . the prototype system we are looking for is thus specified by ( [ master ] ) and ( [ transrate ] ) together with    @xmath53 { }    \gamma = \tanh \left(\frac{2j}{t } \right ) , \makebox[0.3 cm ] { }    \delta = \frac{h_0}{t}.\ ] ]    we emphasise that this system is not an approximation of the ising system ( [ hamiltonian ] ) in the same way as the double - well system ( [ doublewell ] ) is approximated by ( [ mwmaster]),([mwrate ] ) , and ( [ mwparam ] ) , respectively : the ising system and the stochastic glauber model are based on totally different dynamical concepts . nevertheless the devised prototype system appears to capture in a simple way the essence of the effect of coupling and forcing on coupled continuous - variable systems like the double well .    to motivate this , let @xmath54 denote the height of the barrier to be surmounted in order to escape from state @xmath4 . in a very simple approximation the effect of next - neighbour coupling and forcing on the barrier height could be modelled as    @xmath55    @xmath56    comparing this to ( [ hamiltonian ] ) one finds that within this simple picture the effect of coupling and forcing on the barrier heights is the same as their effect on the energy levels of the ising model . ( we note that for coupled double - well systems the chosen approximation would only hold for weak coupling . for strong coupling the systems may no longer be bistable and the entire concept of a simple - minded two - state approximation breaks down . )    within an adiabatic limit and with some noise strength @xmath43 the transition rate @xmath57 of the barrier system would roughly read    @xmath58    from here rates of the form ( [ transrate ] ) together with parameters ( [ params ] ) can be obtained in two steps . first , the factor related to @xmath59 is kept as the parameter @xmath14 with a respective choice of @xmath59 . second , the factor related to @xmath60 , which describes the effect of forcing and coupling , is replaced by the respective term of ( [ transrate ] ) , which is easier to handle . a connection between both types of rates is again made by requiring that they lead to identical stationary distributions in the sense of ( [ detailedbal ] ) . a particular close relation between the glauber - dynamics model and the barrier system should thus arise within the regime of quasistatic response . nevertheless we shall not restrict the following investigation of the simple prototype system to the quasistatic regime . whether or not the suggested model provides insights which are also relevant to more realistic settings and possibly helpful in the design of coupled stochastic resonators remains to be seen . the aim of this section is to calculate the snr of the response of a two - state element embedded into the chain of coupled elements given by ( [ master ] ) and ( [ transrate ] ) . this calculation will be done in terms of the parameters @xmath14 , @xmath29 , and @xmath10 and is essentially based on glauber s work . we shall only take account of terms in leading order of the modulation amplitude @xmath10 . the resulting snr will then be an extension of the simplified mcnamara - wiesenfeld result ( [ mwsnr ] ) . at the end of the section the general result will be applied to the prototype system specified by ( [ params ] ) . in addition to the snr we shall also give an analytical expression for the so - called spectral power amplification ( spa ) . following glauber we start by calculating the averages @xmath61 . from ( [ master ] ) one derives @xmath62 . it results in    @xmath63 \nonumber \\ & & +   \left[1 - \frac{\gamma}{2 } ( r_{k-1,k } + r_{k , k+1})\right ] \delta \cos ( \omega t + \phi)\end{aligned}\ ] ]    with @xmath64 . a closed set of equations for the @xmath65 is found by linearising ( [ mean ] ) in @xmath10 . in this case the @xmath66 can be taken from the unperturbed model ( @xmath67 ) . after long times these @xmath66 read @xmath68 with @xmath69 @xcite . the resulting long - time - limit set of equations can be further simplified . since all elements are forced uniformly , a particular element can not be distinguished from any other , once the initial distribution has been forgotten . hence , all elements will have identical statistics and the indices in ( [ mean ] ) can be skipped . together with the previous assumption one obtains    @xmath70    obviously , the long - time dynamics of the average @xmath71 is identical to the long - time dynamics of the average of an uncoupled element with rescaled relaxation rate @xmath72 and rescaled modulation term @xmath73 . the resulting long - time limit of the averaged state simply reads @xmath74 @xmath75 { } \tan \psi = - \frac{\omega}{\alpha(1-\gamma)},\ ] ]    where @xmath76 is the response to static signals @xmath77 :    @xmath78    formally this static response might even exceed @xmath79 . we conclude that in this case @xmath10 is not sufficiently small and the linearisation of ( [ mean ] ) is no longer justified . later on it will be shown that this does not amount to a considerable restriction . the power spectrum of an element is determined from the average @xmath80 . like in @xcite its @xmath22-dependence is removed by averaging @xmath81 over an uniformly distributed initial phase @xmath12 of the modulation term in ( [ transrate ] ) . introducing the correlation function @xmath82 [ \sigma_k(t+\tau)-\langle \sigma_k(t+\tau ) \rangle ] \rangle $ ] one finds    @xmath83    the second term on the r.h.s . contributes a delta function with weight @xmath84 at signal frequency @xmath11 to the ( one - sided ) power spectrum as defined below . the first term , on the other hand , forms the continuous part of the spectrum .    for the purpose of calculating the snr in leading order of the modulation amplitude @xmath10 it is sufficient to approximate this continuous part by the power spectrum of the unperturbed model . this implies that we neglect any possible signal - induced transformation of the continuous part of the spectrum as we did to obtain the simplified mcnamara - wiesenfeld result ( [ mwsnr ] ) . for the unperturbed model the long - time correlation function @xmath85 is not effected by the average @xmath81 . it also no longer depends on the index of the element . according to glauber @xmath85 reads    @xmath86    with modified bessel functions @xmath87 . using the relation @xcite    @xmath88    which holds for @xmath89 and @xmath90 the power spectrum @xmath91 is calculated . for the one - sided spectrum defined by @xmath92 at @xmath93 we eventually obtain    @xmath94    with @xmath95 and @xmath96 .    now the snr of the response of a two - state element embedded into the chain can be calculated , which is the central result of this section . it reads    @xmath97    for vanishing coupling it reduces to the simplified mcnamara - wiesenfeld result ( [ mwsnr ] ) . in addition to the effect of the coupling on the snr we also wish to study its impact on the delta peak in the power spectrum . in this regard the spectral power amplification ( spa ) is a convenient measure ( cf . e.g. jung in @xcite ) . it is defined as the ratio of power contained in the signal peaks of output to input spectrum , i.e. , the ratio of the weights of the respective delta functions . because of its dependence on the input spectrum the spa depends on the precise dynamics to be modelled by the two - state chain .    for the prototype system specified by ( [ params ] ) the spa reads @xmath98 , or explicitly    @xmath99    @xmath100    where @xmath101 is the spa of static signals . obviously , the spa does no longer depend on the signal amplitude @xmath102 . for convenience we also remove the @xmath102-dependence of the snr by considering    @xmath103    for the prototype system the parameter @xmath104 involved in the unperturbed spectrum ( [ spectrum ] ) simplifies to @xmath105 @xcite .    both quantities , spa ( [ spa ] ) and rescaled snr ( [ snr ] ) , only depend on three parameters : temperature @xmath43 , coupling strength @xmath106 , and signal frequency @xmath11 . restrictions arise from @xmath48 , which lead to a simple expression for @xmath10 , as well as from the linearisation of ( [ mean ] ) , which results in an upper bound on @xmath10 as discussed in connection with ( [ meanstatic ] ) . together one finds    @xmath107    it implies that finite @xmath102 place a lower bound on @xmath43 and an upper bound on @xmath106 . however , these restrictions are rather weak : @xmath102 can be made arbitrarily small because its size is immaterial within the present weak - signal approximation . we are now in a position to study the impact of the coupling on the response measures spa and snr of a single resonator embedded into the chain . this will be the subject of the following two sections . the spa ( [ spa ] ) has an unique maximum over temperature @xmath43 and coupling parameter @xmath106 for any time - dependent signal ( [ forcing ] ) . the maximum spa is obtained for a frequency - dependent value @xmath108 which is always positive . in other words : a properly tuned ferromagnetic - type coupling yields the best spa performance of the two - state resonator element which is embedded into the chain .    for our simple model this maximum can be studied analytically . the partial derivatives of the spa with respect to @xmath43 and @xmath106 are found to vanish at pairs @xmath109 given by    @xmath110    and    @xmath111    respectively . for @xmath108 to be real and finite , @xmath112 has to hold . with this restriction ( [ simple2 ] ) has a unique solution @xmath113 for any given @xmath114 . the corresponding @xmath108 is found via ( [ simple ] ) . within the given range of @xmath115 one can easily check on ( [ simple ] ) that @xmath108 is indeed always positive . an inspection of the second derivatives of the spa finally reveals that there is a maximum at @xmath109 . its peak value is implicitly given by    @xmath116    from ( [ simple ] ) , ( [ simple2 ] ) , and ( [ spamax ] ) the following features of the maximum can be derived which are shown in fig.1 : tuning the signal frequency @xmath11 from very large to vanishing small values the temperature @xmath115 falls from @xmath117 to @xmath118 . at the same time , the coupling strength @xmath108 as well as the peak height @xmath119 increase from vanishing small to very large values . in addition , an analysis of the curvatures @xmath120 and @xmath121 , both expressed in terms of @xmath115 , shows that the sharpness of the maximum grows as its height increases . this indicates that to achieve optimal performance the system parameters @xmath43 and @xmath106 have to be tuned with increasing accuracy as the signal frequency decreases . fig.1 also includes a comparison to the spa of an uncoupled element . at given @xmath114 this spa has a maximum located at a temperature @xmath122 , implicitly given by    @xmath123    on the interval @xmath124 , whereby its peak value is found to be @xmath125 . it can be shown from ( [ simple2 ] ) and ( [ simple3 ] ) that @xmath115 always exceeds @xmath122 . the fact that @xmath108 was found to be never zero also implies that @xmath119 always exceeds @xmath126 . however , looking at fig.1 it is obvious that at frequencies @xmath127 both compared quantities of the coupled element approach those of the uncoupled one and @xmath108 approaches zero . the coupling - induced increase in the spa is thus vanishingly small at sufficiently high frequencies of the signal . this comparison allows to distinguish between two different types of spa behaviour : at low frequencies the spa is enhanced under ferromagnetic coupling whereas at high frequencies it is basicly not , although a tiny increase still occurs . both situations are illustrated at selected frequencies in fig.2 and fig.3 , respectively . the low frequency spa exactly reproduces the qualitative effect found analytically by jung et al . @xcite in a system of globally interacting elements . the high frequency behaviour , on the other hand , where the spa does practically not increase under coupling , has to our knowledge not been reported before .    in both cases fig.2 and fig.3 clearly demonstrate that at fixed coupling parameter @xmath106 stochastic resonance occurs : the spa has a maximum over temperature @xmath43 which , unfortunately , can not be established analytically . fig.2(bottom ) shows that while ferromagnetic coupling improves the spa at any fixed temperature @xmath43 this improvement is lost , if the coupling becomes too strong . there is thus an optimal coupling strength @xmath128 for every given temperature , which can even be exactly calculated at any @xmath114 . one finds    @xmath129    which is always positive . since the spa does not have a further extremum over @xmath106 , it implies that antiferromagnetic coupling always decreases the spa . moreover , @xmath128 increases with growing temperature t and decreases as the signal frequency grows . in fig.2 the spa at the optimal coupling strength is included ( dash - dotted curves ) . in fig.3(bottom ) , where the high - frequency spa is illustrated , the coupling - induced improvement of the spa is hardly detectable and @xmath128 is almost zero .    considering the spa of time - independent signals ( @xmath130 ) given by ( [ spa_static ] ) one finds that it has a maximum over temperature for antiferromagnetic coupling @xmath131 , only , whereby it decreases as the coupling strength @xmath132 grows . for @xmath133 the spa increases with increasing @xmath106 as well as with decreasing @xmath43 . it formally diverges for @xmath134 and @xmath135 , respectively . in both cases the weak - signal limit ( [ range ] ) breaks down .    in general , the spa ( [ spa ] ) is given by its static value ( [ spa_static ] ) times a dynamical factor @xmath136 ( cf.([meanresult ] ) ) . here @xmath137 is the ratio @xmath138 of signal frequency to long - time relaxation rate , which also governs the phase shift @xmath139 in ( [ meanresult ] ) . with growing @xmath137 the elements gradually loose their ability to follow the signal : the spa weakens and the phase shift grows . this effect occurs , for example , if the signal frequency @xmath11 increases . subsequently , the spa decreases with growing @xmath11 . the impact of @xmath137 is also responsible for the occurrence of stochastic resonance and optimal coupling in the spa . here a decrease of the long - time relaxation rate @xmath72 plays the crucial role : the elements dynamics slows down as @xmath106 increases or as @xmath43 decreases . thus this slow - down occurs whenever the static spa grows . hence , the dynamical factor always counteracts the static spa as @xmath43 or @xmath106 is changed . eventually , the increase of the monotonous static spa is outperformed by the decrease of the dynamical factor which results into a maximum of the spa over @xmath43 and @xmath106 , respectively . in other words , stochastic resonance and optimal coupling occur .    as illustrated in fig.2 the spa can be seen as a transition between two limits , the static spa given by @xmath140 ( @xmath141 , dashed curves ) and @xmath142 ( @xmath143 , dotted curves ) , respectively . since both limits intersect at @xmath144 , i.e. at @xmath145 , the plots nicely show the well - known approximate matching of time scales at the spa peak . this matching does not only occur over @xmath43 ( fig.2 , top ) , but over @xmath106 ( fig.2 , bottom ) , too . if the signal frequency is changed , the peaks of the spa shift whereby all curves share the static spa as a limit . this results into the qualitative frequency dependence found by mcnamara and wiesenfeld for the delta function part of the spectrum of the double - well system ( [ doublewell ] ) .    roughly speaking , an improvement of the spa under ferromagnetic coupling only occurs at @xmath146 , where the impact of the static spa is not yet outperformed by the dynamical factor . this explains why there is no improvement at signal frequencies @xmath147 : since under ferromagnetic coupling @xmath148 holds , @xmath146 can not be fulfilled in this case . at @xmath149 ( fig.3 ) the spa is already within line width of the limit @xmath150 , where the impact of the dynamical factor prevails . this limit is basicly given by the modulation term in ( [ simplemean ] ) . its non - monotonous temperature dependence is known to be no longer associated with a matching of time scales @xcite . the static response , which is thus the origin of the coupling - induced enhancement of the spa , can be interpreted as follows . without coupling and for a given fixed signal the two - state elements prefer to be in the state of lowest energy ( [ hamiltonian ] ) or highest barrier ( [ deltau ] ) , respectively . since the signal is homogeneous , this state is the same for all elements . without a signal but under ferromagnetic coupling neighbouring elements tend to be in the same state , too , although they do not favour a particular state @xmath151 or @xmath152 . taking now signal and coupling together , both effects add up . the tendency to find the elements in the state favoured by the signal grows and hence @xmath76 and the static spa grow compared to the uncoupled case .    with antiferromagnetic coupling , on the other hand , neighbouring elements prefer to be in opposite states . this counteracts the effect of the signal and leads to a decrease in the spa . finally , the decrease of the static spa with growing temperature directly follows from ( [ detailedbal ] ) . there the imbalance in the distribution of probability between the two states , which dependent on the view taken either differ in energy or barrier height , decreases as the system heats up : the mean response @xmath76 weakens due to increasing fluctuations . the investigation of the snr presented here basicly relies on numerical evaluations of ( [ snr ] ) . an explicit calculation of the power spectrum ( [ spectrum ] ) already yields a rather complicated expression which does not lend itself to a detailed analytical study .    at finite signal frequencies @xmath114 the snr displays the same qualitative dependence on temperature @xmath43 and coupling parameter @xmath106 as shown for the spa in fig.2 and fig.3 , respectively . at low frequencies the snr is enhanced under ferromagnetic coupling whereas at high frequencies it is not . due to this close similarity we omitted the respective plots for the snr . we note , however , that this similarity is by no means a trivial result , since the spectrum ( [ spectrum ] ) is itself a non - monotonous function of @xmath43 and @xmath106 . the low - frequency snr qualitatively reproduces the snr behaviour found in a chain of next - neighbour - coupled overdamped double - well systems simulated in @xcite . it is the first analytical confirmation of this behaviour . we note that this correspondence occurs although the simulation in @xcite was performed with strong forcing while the present model is studied within a weak - signal limit . as for the spa , the high - frequency behaviour of the snr has to our knowledge not been reported before . if the elements are not coupled , their snr ( [ mwsnr ] ) does not depend on the signal frequency . this was already found in @xcite , if the signal - induced reduction of the continuous part of the spectrum is neglected , as ( [ mwsnr ] ) clearly shows . under coupling this frequency independence of the snr for ferromagnetic coupling numerical evaluations of ( [ snr ] ) predict a decrease of the snr with growing signal frequency . only in the limit of low and high frequencies the snr approaches constant values . in both cases the snr formula ( [ snr ] ) simplifies significantly . one finds    @xmath153    @xmath154    where @xmath155 is the simplified mcnamara - wiesenfeld snr ( [ mwsnr ] ) of uncoupled elements . @xmath156 is the static snr while @xmath157 represents the leading order term of its high - frequency expansion . fig.4 illustrates the typical qualitative behaviour of the snr as a function of the signal frequency . at various frequencies it shows a set of snr curves over temperature @xmath43 and coupling parameter @xmath106 , respectively . the dashed curves represent @xmath156 and @xmath157 , respectively . the apparent fact that the snr decreases with growing signal frequency implies that its static value can not be exceeded at any other signal frequency . then it follows immediately from ( [ snrlim1 ] ) that the coupling - induced enhancement of the snr possesses an upper limit : it can not be better than a factor four . this is a rigorous result for the general chain model and not limited to the prototype system , where a special choice of the parameters was made . due to the simplicity of the chain model this result may well reflect a limit of more general nature for the improvement of a weak - signal snr . fig.5 shows the static snr over temperature @xmath43 . it differs qualitatively from the static spa which according to ( [ spa_static ] ) has a monotonous behaviour .    turning briefly to antiferromagnetic coupling one finds that the continuous part of the spectrum is insensitive to the sign of the coupling parameter @xmath106 . now even the static snr decreases under coupling . as shown in fig.4(bottom ) it is found below the high - frequency snr curve , i.e. , here the snr increases with growing signal frequency . since at medium frequencies again a transition occurs between both curves similar to the transition in fig.4(top ) , a second maximum in the snr over @xmath43 may emerge in this situation . the investigation of the impact of the coupling on the snr can also be extended to the more natural situation , where the signal itself is embedded into noise . the question to address is , whether the improvement of the snr , which was established for independent internal noise sources , can still be found with additional coherent external noise . ( this problem does not occur for the spa , where only the height of the signal peak is of interest . )    to study this case it is assumed that the input spectrum consists of the previous signal peak described by @xmath158 and a noise part @xmath159 . within the weak - signal limit it was shown that the signal peak of the input spectrum leads to a respective peak in the output spectrum at signal frequency , only : no additional peaks at multiples of that frequency occur . therefore , any additional contribution to the continuous part of the output spectrum at signal frequency can only arise from @xmath160 . in analogy to ( [ meanresult ] ) one finds that this additional contribution is given by @xmath161 .    together with ( [ snr ] ) the resulting snr @xmath162 reads , again in units of @xmath163 , @xmath164 it can be expressed in terms of the snr @xmath165 ( cf.([snr ] ) ) and the input snr @xmath166 written in units of @xmath163 , respectively . one finds @xmath167    @xmath162 is thus a steadily growing function of the previously studied snr @xmath165 , whereby it can not exceed the input snr @xmath168 . ( the latter is a meanwhile well - known result of linear response theory @xcite which is in fact the limit we are taking here . ) since @xmath168 is constant at fixed signal frequency @xmath11 , a coupling - induced increase of @xmath165 will lead to an increase in @xmath162 , too . the improvement of the snr under coupling is thus preserved with external coherent noise . from ( [ noisysnr2 ] ) one can easily show that the maximum enhancement of @xmath162 compared to the uncoupled chain reaches the previously found factor four for vanishing input noise , only . in general , this factor is smaller approaching 1 , if the input noise is so strong that the input snr and hence @xmath162 go to zero .    at the end we shall make a brief excursion to the collective response of the chain , which was mentioned in the introduction . we consider the following sum which involves the states of @xmath0 elements :    @xmath169    it is not difficult to show that the spa of this new quantity is simply @xmath170 times the previously studied spa ( [ spa ] ) . for the snr , however , this new situation is completely different . previously the continuous part of the spectrum was determined from the auto - correlation function @xmath171 of an element alone ( cf . ( [ correl1 ] ) ) . now this spectrum will involve cross - correlation contributions @xmath172 of different elements @xmath173 , too . the new contributions to the continuous part of the spectrum change the qualitative behaviour of the snr under coupling . we shall demonstrate this in the limit @xmath174 . for this case glauber calculated the unperturbed spectrum @xcite which reads in its one - sided version    @xmath175    inserting this spectrum into ( [ snr ] ) , multiplied by @xmath170 due to the mentioned increase of the spa , one finds for the new snr per element @xmath176    which , of course , again reduces to the simplified mcnamara - wiesenfeld snr ( [ mwsnr ] ) for vanishing coupling . clearly , the snr is now a decreasing function of the coupling parameter and the sign of the latter is no longer important . hence , there is a drastic difference in the snr of single and collective response with respect to coupling . of course , we have so far merely investigated the limits of ( [ global0 ] ) , @xmath177 and @xmath174 , respectively . we expect that the snr of the collective response of only a few elements still increases under coupling @xcite . a detailed analysis of the present model with respect to these collective effects will be the subject of a further investigation . in this paper we extended the two - state model of stochastic resonance introduced by mcnamara and wiesenfeld to a chain of infinitely many coupled two - state elements which are periodically modulated . the interaction of the elements was chosen in a way that the chain evolves according to glauber s stochastic ising model . in analogy to the work of mcnamara and wiesenfeld and based glauber s results analytical expressions for the signal - to - noise ratio ( snr ) and the spectral power amplification ( spa ) have been obtained in the limit of weak modulations . here both quantities refer to the response of a single element as part of the chain .    instead of approximating the dynamics of a particular coupled system by the chain model , we used the latter to build a prototype system which hopefully captures the essential features of an entire class of coupled stochastic resonators . to this end additional assumptions were made on the dependence of the model parameters on some noise intensity and signal amplitude , respectively . the prototype system was used to study the effect of the coupling on the response of a single resonator embedded into the chain . the results show that array - enhanced stochastic resonance occurs for ferromagnetic - type coupling in spa and snr . the qualitative features of the effects reproduce those previously found in coupled stochastic resonators . the simplicity of the chain model allowed for a detailed analytical investigation of the spa . for the snr the observed effects have been confirmed analytically for the first time . in addition , it was found that the improvement of the single - resonator response compared to the response of the uncoupled resonator still occurs , if the signal is embedded into noise . for the snr this improvement was shown to be limited by a factor four which is reached for vanishing input noise , only . a brief excursion into the collective response of @xmath0 resonators , on the other hand , disclosed that coupling can not improve the snr , if @xmath0 is very large .    a closer look at the mechanisms behind the effects revealed that in the present model an improvement of the stochastic resonators under coupling is associated with the regime of quasistatic response . since the model studied here possesses an upper bound to the time scale of its dynamics , the desired improvement is essentially restricted to sufficiently slow signals . for the spa this improvement is based on a stronger tendency of the two - state elements to align in parallel , if signal and ferromagnetic coupling act together , compared to this tendency caused by the signal alone . the reason why the improvement of the spa is lost , if the coupling is too strong , was found to be the slow - down of the system dynamics under coupling : it simply prevents the resonators from responding quasistatically .          for review - type articles on stochastic resonance see , e.g. : + p. jung , phys . * 234 * , 175 ( 1993 ) . + f. moss , d. pierson , d. ogorman , int . 4*(6 ) , 1383 ( 1994 ) . + k. wiesenfeld , f. moss , nature * 373 * , 33 ( 1995 ) . + a. bulsara and l. gammaitoni , physics today , march 1996 , 39 ( 1996 ) . + l. gammaitoni , p. hnggi , p. jung , and f. marchesoni , rev . * 70 * , 223 ( 1998 ) + a bibliography on stochastic resonance can be found at http://www.pg.infn.it / sr/.    j. k. douglas , l. wilkens , e. pantazelou , and f. moss , nature * 365 * , 337 ( 1993 ) . + s. m. besrukov and i. vodyanoy , nature * 378 * , 362 ( 1995 ) . + j. e. levin and j. p. miller , nature * 380 * , 165 ( 1996 ) . + p. cordo , j. t. ingils , s. verschueren , j. j. collins , d. m. merfeld , s. rosenblum , s. buckley , and f. moss , nature * 383 * , 769 ( 1996 ) . + j. j. collins , t. t. imhoff , and p. grigg nature * 383 * , 770 ( 1996 ) . r. rouse , siyuan han , and j. e. lukens , appl . phys . lett . * 66 * , 108 ( 1995 ) . + a. d. hibbs , a. l. singsaas , e. w. jacobs , a. r. bulsara , j. j. bekkedahl , and f. moss , j. appl . * 77 * , 2582 ( 1995 ) . + a. d. hibbs and b. r. whitecotton , in _ applied nonlinear dynamics and stochastic systems near the millennium _ , ed . by j. b. kadtke and a. bulsara , ajp conference proceedings , vol . 411 , p.227 ( aip , woodbury , new york , 1997 ) . a. bulsara , g. schmera , phys . rev . * e47 * , 3734 ( 1993 ) . + a. bulsara , a. maren , and g. schmera , biol . cybern . * 70 * , 145 ( 1993 ) . + m. e. inchiosa , a. r. bulsara : phys . lett . * a200 * , 283 ( 1995 ) and phys . rev . * e52 * , 327 ( 1995 ) .
Abstract:

Answer: the two - state model of stochastic resonance is extended to a chain of coupled two - state elements governed by the dynamics of glauber s stochastic ising model . appropriate assumptions on the model parameters turn the chain into a prototype system of coupled stochastic resonators . in a weak - signal limit analytical expressions are derived for the spectral power amplification and the signal - to - noise ratio of a two - state element embedded into the chain . the effect of the coupling between the elements on both quantities is analysed and array - enhanced stochastic resonance is established for pure as well as noisy periodic signals . the coupling - induced improvement of the snr compared to an uncoupled element is shown to be limited by a factor four which is only reached for vanishing input noise .