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 ) .