rotational brownian motion of colloidal particles is a classical subject in statistical physics @xcite .
contrary to its translational cousin , _
i.e. _ `` usual '' brownian motion , it does not result in changes of the particle location and is hence not as easily demonstrated in experiments .
if , however , the suspended particles carry an electric or magnetic moment their orientation couples to external fields and the electric or magnetic relaxation properties of the suspension are direct consequences of the rotational diffusion of the constituting particles . in particular , in the well developed theory of ferrofluids the importance of rotational brownian motion for the theoretical description of the often complex and surprising hydrodynamic and magnetic properties was recognized already a long time ago @xcite .
ferrofluids are suspensions of nano - size ferromagnetic grains in carrier liquids like water or oil .
they combine the hydrodynamic properties of newtonian fluids with the magnetic behaviour of superparamagnets @xcite . due to the viscosity of the carrier liquid
there is a coupling between the rotation of the magnetic grains and the local vorticity of the hydrodynamic flow .
this coupling can , _ e.g. _ , be used to spin up ferromagnetic drops floating in a non - magnetic fluid of the same density with the help of a _ rotating magnetic field _ @xcite .
complementary , a rotational ferrofluid flow exposed to a static magnetic field exhibits an enhanced shear viscosity @xcite .
various related effects such as `` negative '' rotational viscosity @xcite , magneto - vortical resonance @xcite , and anomalously enhanced ac - response due to coherent particle rotation @xcite have been investigated .
they rely on the exchange of angular momentum between rotating particles and an oscillating magnetic field . in these cases ,
the imposed non - zero _ flow vorticity _ of the ferrofluid is crucial for breaking the symmetry between clockwise and counter - clockwise particle rotation .    in the present paper
we investigate a much more indirect and subtle aspect of the interplay between rotational brownian motion of ferrofluid particles and their relaxation dynamics in an external magnetic field .
we will show that a suitably designed time dependent external magnetic field _ without net rotating component _ may rectify the fluctuations of the particle orientation and set up a _ noise - induced _ rotation of the ferromagnetic grains .
we will hence investigate how angular momentum can be transferred from an _ oscillating _ magnetic field to a ferrofluid _ at rest_. the effect was predicted theoretically and demonstrated experimentally in a previous short communication @xcite . in the present paper we provide much more details on the theoretical description and add several new results .
the extraction of directed motion from random fluctuations is an old and controversial problem in statistical mechanics with a long and interesting history @xcite .
although excluded by the second law of thermodynamics for equilibrium systems , rectification of fluctuations _ is _ possible in systems driven sufficiently far away from thermal equilibrium @xcite .
the problem has gained renewed attention under the trademarks of `` thermal ratchets '' and `` brownian motors '' due to its possible relevance for biological transport @xcite and the prospects of nano - technology @xcite .
ferrofluids are ideal systems to investigate such fluctuation driven transport phenomena and also to demonstrate them experimentally @xcite : as already discussed above , the rotational dynamics of the ferromagnetic grains is strongly influenced by thermal fluctuations .
appropriate time - dependent potentials can be easily designed with the help of external magnetic fields . finally , directed rotational transport in ferrofluids should manifest itself as systematic rotation of the ferromagnetic nano - particles .
this in turn can be easily detected from the resulting _ macroscopic _ torque on the carrier liquid .
various somewhat related but still quite different phenomena in rotational dynamical systems are treated in @xcite .    throughout the paper
we will use two basic approximations which simplify the analysis considerably and which are rather common for ferrofluids @xcite .
the first is to neglect neel - relaxation of the magnetization , _
i.e. _ the rotation of the magnetization vector with respect to the ferromagnetic particle .
this is justified for particle sizes that are not too small and amounts to assuming that the magnetic moments are firmly attached to the geometry of the particles .
any reorientation of the magnetic moment hence requires a rotation of the particle as a whole .
the second approximation is to neglect dipole - dipole interactions between the particles .
although these interactions may be important in concentrated ferrofluids , they are negligible in sufficiently diluted ones .
the ratchet mechanism of central interest in the present investigation can operate without any particle - particle interaction .
we will hence assume that our ferrofluids are sufficiently diluted and use a _ single particle _ model .
for a more quantitative assessment of the role of the dipole - dipole interactions see _ e.g. _ @xcite and references therein .
the paper is organized as follows . in section [ sec : model ]
we introduce the general framework for describing the rotational brownian motion of a single suspended ferromagnetic particle in an external magnetic field .
section [ sec : symm ] gives a detailed account of the symmetries characteristic of our system and discusses how these symmetries influence the possibility for a ratchet mechanism to operate . in section [ sec : fpe ] the central fokker - planck equation describing our system is analyzed . besides an investigation of the weak noise limit
we describe the numerical solution of the fokker - planck equation and compare the results with two different approximate treatments yielding analytical expressions for the transferred angular momentum .
finally , section [ sec : conc ] contains some conclusions .
we consider the rotational motion of a spherical particle of volume @xmath0 and magnetic moment @xmath1 immersed in a liquid with dynamic viscosity @xmath2 and subject to a horizontal , time dependent , spatially homogeneous magnetic field @xmath3 .
the field is composed of a constant part @xmath4 parallel to the @xmath5-axis and an oscillatory part @xmath6 with period @xmath7 along the @xmath8-direction @xmath9 different choices for the time dependence of @xmath10 are of interest . in the present paper
we will mainly discuss two cases which are each representative for a whole class .
our first standard choice is @xmath11 where the amplitudes @xmath12 , and the phase @xmath13 are control parameters .
the main features of this time dependence are a zero average over one period and the presence of a higher harmonic of the basic frequency . as a second example
, we will also discuss the form @xmath14 for which the average over one period is different from zero .
it arises naturally if the constant field component is not perpendicular to the time dependent one . in any case , the magnetic field is of a pure oscillatory character , _
i.e. _ it _ does not contain a net rotating component_.    the orientation of the particle at time @xmath15 is described by the unit vector @xmath16 where @xmath17 denotes the modulus of the magnetic moment .
the time evolution of @xmath18 is given by @xmath19 where @xmath20 denotes the instantaneous angular velocity of the particle .
changes of @xmath21 are due to torques on the particle .
denoting by @xmath22 the angular part of the three dimensional nabla operator , the magnetic torque @xmath23 derives from the potential energy @xmath24 of a magnetic dipole in an external field @xcite .
further , the viscosity @xmath2 of the carrier liquid gives rise to a viscous torque @xcite @xmath25 additionally , the interaction between the rotating particle and the surrounding liquid also causes thermal fluctuations which generate a stochastic torque @xcite @xmath26 here , @xmath27 is a vector of independent , @xmath13-correlated gaussian noise sources of zero mean , the noise intensity is related to the temperature @xmath28 and the dissipation ( [ defvisctor ] ) of the carrier liquid by the fluctuation - dissipation relation , and @xmath29 stands for boltzmann s constant . denoting the moment of inertia of the particle by @xmath30 , the equation of motion for @xmath21 acquires the form @xmath31 eqs.([evol1 ] ) and ( [ evol2 ] ) form a closed set of equations for the description of the rotational motion of the particle . using experimentally relevant parameter values ( density of the particle @xmath32 kg / m@xmath33 , particle radius @xmath34
nm , viscosity @xmath35 pas , time scale @xmath36 s ) we find that the first term on the l.h.s .
( [ evol2 ] ) is five to seven orders of magnitude smaller than the second one .
we may hence safely neglect inertial effects @xcite and find for @xmath21 in the overdamped limit @xmath37 where we have introduced @xmath38 .
using this result in ( [ evol1 ] ) yields a closed equation for the time evolution of @xmath18 @xmath39 it is convenient to introduce dimensionless units . to this end
we measure time in units of the inverse of the external driving frequency , @xmath40 , use @xmath41 as unit for the magnetic field strength , @xmath42 , and rescale the noise intensity according to @xmath43 . the evolution equation for the orientation
@xmath18 then reads @xmath44 introducing the brownian relaxation time @xmath45 we note that the rescaled , dimensionless noise intensity @xmath46 occurring in ( [ evol4 ] ) just gives the ratio between the relevant deterministic and stochastic time scales in the system : @xmath47    to proceed , we parametrize the orientation @xmath18 of the magnetic particle by two angles @xmath48 and @xmath49 according to @xmath50 from ( [ evol4 ] ) we then find the following langevin equations for the time evolution of these angles @xcite @xmath51 in the dimensionless units adopted , the noise intensity @xmath46 is given by ( [ defd ] ) and the potential ( [ defu ] ) takes the form @xmath52 the thermal fluctuations @xmath53 and @xmath54 are given by two independent , @xmath13-correlated gaussian noise sources of zero mean .
note that in ( [ lan2 ] ) we are dealing with multiplicative white noise but since the multiplicative function is independent of @xmath55 , no ambiguity ( ito- vs. stratonovich - interpretation @xcite ) arises .
it is instructive also to consider the simplified situation of a one - dimensional dynamics in which the vector @xmath18 is assumed to be constrained to the horizontal @xmath5-@xmath8-plane .
accordingly , ( [ lan1 ] ) is replaced by @xmath56 and ( [ lan2 ] ) simplifies to @xmath57 where @xmath58    the observable of foremost interest in the present investigation is the time and ensemble averaged torque ( [ defmagtor ] ) exerted by the magnetic field upon the magnetic particle in the long time limit , _
i.e. _ after initial transient have died out .
since the magnetic field ( [ deffield ] ) is constrained to the @xmath5-@xmath8 plane , only the @xmath59-component of this magnetic torque can be different from zero . suppressing the subscript @xmath60 in ( [ defmagtor ] ) from now on
, we denote the averaged @xmath59-component of the magnetic torque by @xmath61 , where @xmath62 stands for the ensemble average over the different realizations of the noise terms in ( [ lan1 ] ) , ( [ lan2 ] ) and the overbar represents the time average over one period of the magnetic field . using ( [ deffield ] ) and ( [ defmagtor ] ) we get @xmath63 exploiting ( [ defpot ] ) and ( [ lan2 ] )
one readily finds the equivalent expressions @xmath64 for reasons of ergodicity , the ensemble average in ( [ torav ] ) is equivalent to a time average of a single realization over an infinite time interval .
then the extra time average over one period of the external driving drops out and we are left with @xmath65 and similarly for the equivalent expressions in ( [ torav ] ) , ( [ torav ] ) .
space - time plot of the potential @xmath66 from eq.([defpot ] ) for the time dependence ( [ foft1 ] ) with @xmath67 , and @xmath68 .
dark and bright regions correspond to small and large values of @xmath69 , respectively . in the long - time limit ,
the deterministic dynamics ( [ lan1 ] ) , ( [ lan2 ] ) with @xmath70 approaches @xmath71 and a periodical @xmath55 , oscillating back and forth as represented by either of the full black lines . in the presence of a small amount of noise ,
occasional transitions across the unstable deterministic orbits shown as full white lines become possible which are schematically indicated by the dashed lines . the spatial asymmetry and temporal anharmonicity of
the potential conspire to yield slightly different rates for noise induced increments and decrements of @xmath49 , respectively . as a result a noise driven rotation of the particles arises . for a detailed discussion see section [ inst1d ] . ]    in the absence of thermal fluctuations , the particle orientation is governed by the overdamped deterministic relaxation dynamics given by ( [ lan1 ] ) , ( [ lan2 ] ) with @xmath70 .
hence , at any given moment , the orientation @xmath72 tends to align with the instantaneous magnetic field @xmath73 which lies in the @xmath5-@xmath8-plane .
one can then show @xcite that @xmath74 converges to @xmath75 for @xmath76 and that @xmath55 approaches a _ periodic _ long time behaviour .
this implies via ( [ torav  ] ) that @xmath77 _
i.e. _ , in the absence of thermal fluctuations no particle rotation will occur and no average torque can arise .
an explicit example is displayed in fig .  [
fig : fig1 ] .
this scenario changes fundamentally in the presence of fluctuations , _
i.e. _ if @xmath78 . as shown qualitatively in fig .
[ fig : fig1 ] , for small noise intensities , the time dependence of @xmath55 and @xmath74 will still closely follow the deterministic trajectories , except that now also rare , fluctuation induced _ transitions _ between different deterministic solutions may occur . if the rates of forward ( increasing @xmath49 by @xmath79 ) and backward transitions ( decreasing @xmath49 by @xmath79 ) are different from each other , then on the average a net rotation of the particle will occur , implying with ( [ eq24 ] ) that @xmath80 .
this is a manifestation of the so - called ratchet effect @xcite in which an unbiased potential and undirected brownian fluctuations cooperate to produce directed transport .
the detailed operation of the ratchet mechanism for the rotational motion of a colloidal particle is the main focus of the present paper .
symmetry considerations admit conclusions from the equations of motion without actually solving them . for
the investigation of ratchet systems , studies of symmetry turned out to be valuable since directed transport generically occurs if it is not forbidden by symmetries , a statement referred to as curie s principle @xcite . here
we specifically address the question which transformations of the magnetic field components @xmath4 and @xmath6 in ( [ deffield ] ) leave the average @xmath59-component of the torque invariant , @xmath81 , and which lead to its inversion , @xmath82 .    from eq .
( [ torav  ] ) one readily concludes that the following transformations of @xmath49 and @xmath48 imply the indicated transformation behavior of @xmath61 : @xmath83 where @xmath84 and @xmath85 are arbitrary . observing that the values of both the original and the transformed @xmath74 are restricted to the interval @xmath86 $ ] and that the entire interval will indeed be sampled by some realization of the dynamics ( [ lan1 ] ) , ( [ lan2 ] )
, it follows that no transformations of the form @xmath87 with @xmath88 or @xmath89 with @xmath90 are possible .
finally , in order to conclude @xmath82 in ( [ 1e ] ) one has to exploit that in ( [ torav  ] ) the limit @xmath91 is not necessarily tantamount to @xmath92 with @xmath93 kept fixed , but can also be realized by letting @xmath94 with @xmath95 kept fixed .
closer inspection of ( [ torav  ] ) reveals that any further symmetry transformation that leaves @xmath61 invariant or changes its sign can be composed of the `` elementary '' transformations ( [ 1a])-([1e ] ) .
next we ask ourselves what equations of motion govern the transformed @xmath96 , @xmath55 defined by ( [ 1a])-([1e ] ) . in the simplest case ( [ 1a ] )
one readily concludes from ( [ lan1 ] ) and ( [ lan2 ] ) that the transformed @xmath96 , @xmath55 satisfy the same equations of motion ( [ lan1 ] ) and ( [ lan2 ] ) but with a transformed potential @xmath97 with arbitrary @xmath98 .
turning to the second transformation ( [ 1b ] ) one finds similarly that the transformed @xmath96 and @xmath55 satisfy the equations of motion ( [ lan1 ] ) and ( [ lan2 ] ) but now with a transformed potential @xmath99 and in addition with a transformed noise @xmath100 . since the statistical properties of the two noises @xmath101 and @xmath54 are equal , _
i.e. _ they are identical stochastic processes , the same follows for the corresponding processes given by ( [ lan1 ] ) and ( [ lan2 ] ) . in other words , to each realization of the original @xmath96 , @xmath55 corresponds a realization of the transformed @xmath96 , @xmath55 occurring with the same probability , and vice versa .
since the average torque in ( [ torav  ] ) is independent of the specific realization @xmath96 , @xmath55 with probability 1 ( _ i.e. _ up to a set of realizations of measure zero ) we can conclude that in this case @xmath82 .
similar considerations can be done for the cases given by ( [ 1c ] ) and ( [ 1d ] ) . the last transformation ( [ 1e ] ) is insofar special as a mapping between the dynamics of the original and the transformed variables is possible only , if the term @xmath102 in ( [ lan1 ] ) identically vanishes , _
i.e. _ if @xmath103 .
we therefore find the following implications of transformations in the potential @xmath104 : @xmath105 where @xmath98 is an arbitrary constant .    finally , we have to determine all those transformations of the magnetic field components @xmath4 and @xmath6 in ( [ deffield ] ) which correspond via ( [ defpot ] ) to transformations of the potential @xmath106 given in ( [ 2a])-([2e ] ) .
since changes in @xmath4 and @xmath6 always imply changes in the potential @xmath69 that depend on @xmath48 and @xmath49 , we can only induce transformations with @xmath107 .
the only non - trivial way to realize ( [ 2a ] ) is @xmath108 combined with @xmath109 with an arbitrary integer @xmath110 .
similarly , ( [ 2b ] ) and ( [ 2d ] ) are equivalent to @xmath111 and @xmath112 . on the other hand ( [ 2c ] ) only admits the trivial realization @xmath113 .
finally , ( [ 2e ] ) can be implemented by @xmath114 . in summary , we find that @xmath115 where the trivial result ( [ 3c ] ) is listed only for the sake of completeness .    while ( [ 3a])-([3d ] ) are intuitively more or less obvious , ( [ 3e ] ) respectively ( [ 2e ] ) is not .
the additional condition @xmath116 shows that the constrained dynamics ( [ lan3 ] ) , ( [ deff ] ) has an extra symmetry which is lost if the full two - dimensional dynamics ( [ lan1 ] ) , ( [ lan2 ] ) is considered .
we note that ( [ 3a])-([3e ] ) could also have been obtained in a more direct way .
our somewhat more involved line of reasoning has the advantage that we can exclude that there are any transformations other than ( [ 3a])-([3e ] ) and combinations thereof which would leave @xmath61 invariant or change its sign .    in order to see that the above considerations already allow some non - trivial conclusions about the possibility of angular momentum transfer from the oscillating magnetic field to the particle we combine ( [ 3a ] ) and ( [ 3b ] ) , ( [ 3b ] ) and ( [ 3d ] ) , and ( [ 3a ] ) and ( [ 3e ] ) to obtain @xmath117 respectively .
from ( [ 3f ] ) we immediately can infer that the constant field in @xmath5-direction is indispensable for a non - zero average torque , @xmath118 similarly ( [ 3 g ] ) implies that for a time dependence obeying @xmath119 for some @xmath85 the average torque has to vanish as well @xmath120 choosing in particular @xmath121 , _ i.e. _ half the period of the external driving , we realize that an oscillating field @xmath6 with a fourier expansion containing only odd multiples of the basic frequency @xmath122 will result in a zero average torque @xmath61 , irrespective of the particular choices of the amplitudes @xmath123 and phases @xmath124 .
eqs.([4a ] ) and ( [ 4c ] ) motivate our choices ( [ foft1 ] ) and ( [ foft2 ] ) as simple time dependencies of the oscillating field resulting in a non - zero average torque @xmath61 .    for the special time dependence ( [ foft1 ] ) we find from ( [ 3 g ] ) @xmath125 as a consequence , upon continuously varying @xmath13 we can infer that there must exist a @xmath126 such that @xmath127 for @xmath128 or @xmath129 .
it is important to note that this case of zero averaged torque is qualitatively different from the situations described by eqs .
( [ 4a ] ) and ( [ 4b ] ) since it is due to the _ fine tuning of a parameter _ rather than resulting from an underlying symmetry .
consequently , upon variation of @xmath13 around @xmath130 the average torque @xmath61 changes sign , which is a realization of a so - called _ current inversion _ @xcite .
furthermore , by fixing @xmath128 , a sign change of @xmath61 will generically also occur upon variation of any other parameter of the system .
interesting examples for such a parameter include the driving frequency @xmath131 , the noise intensity @xmath46 , and the particle volume @xmath0 . in the latter case
we hence face the intriguing possibility that in a polydisperse ferrofluid under the same experimental conditions the larger and smaller particles will rotate in _ opposite _ directions .
we finally note that in the case where the dynamics is constrained to the @xmath5-@xmath8-axis , the additional symmetry ( [ 3h ] ) implies that @xmath132 .
consequently the @xmath13-values where the torque changes sign are fixed by the additional symmetry and do not depend on the other parameters of the problem .
hence no current reversal upon changes in these parameters is possible .
this is the main difference between the full two - dimensional dynamics ( [ lan1 ] ) , ( [ lan2 ] ) and its simplified one - dimensional counterpart ( [ lan3 ] ) , see also @xcite .
a quantitative analysis of the overdamped brownian rotation of a single ferrofluid particle in a time dependent external field can be built on the numerical simulation of the langevin equations ( [ lan1 ] ) , ( [ lan2 ] ) .
such simulations unambiguously show the occurrence of noise induced rotation of the particle in an oscillating field @xcite . for a detailed study of the influence of the various parameters of the problem on the transferred torque @xmath61 it is , however , more convenient to start with the equivalent fokker - planck equation ( fpe ) for the probability density @xmath133 for the orientation ( [ defe ] ) of the particle @xmath134 here
@xmath135 denotes the angular part of the three dimensional nabla operator .
the equivalence of ( [ lan1 ] ) , ( [ lan2 ] ) and ( [ fpe ] ) is , _
e.g. _ , shown in @xcite . for a time - periodic potential @xmath104 , the solution @xmath136 of ( [ fpe ] ) will also be periodic after initial transients have died out . from this solution
we can determine the average orientation of the particle @xmath137 with the help of which the average torque from ( [ torav])-([torav ] ) can be calculated according to ( [ defmagtor ] ) in yet another way , namely @xmath138_z .\ ] ]    in the present section we discuss the numerical and approximate analytical solution of the fokker - planck equation ( [ fpe ] ) and the corresponding results ( [ eq48 ] ) , ( [ torfpe ] ) for the averaged torque @xmath61 . to first substantiate the intuitive understanding of the ratchet effect in our system advocated with fig .
[ fig : fig1 ] , we start by considering the weak noise limit @xmath139 . to keep the analysis simple
we will restrict ourselves in this part to the simplified one - dimensional model defined by ( [ lan3 ] ) , ( [ deff ] ) . in the second subsection we detail the numerical methods used to solve the fpe for the general case .
finally we present some analytical results from the perturbative solution of the fpe .
our starting point is the simplified one - dimensional model dynamics ( [ lan3 ] ) , ( [ deff ] ) with a @xmath13-correlated , unbiased gaussian noise source @xmath54 .
the deterministic dynamics described by ( [ lan3 ] ) with @xmath70 has a family of stable and unstable periodic orbits @xmath141 and @xmath142 , respectively , where @xmath110 is an arbitrary integer ( cf .
[ fig : fig1 ] ) @xcite .
hence deterministically the particle simply follows the direction of the magnetic field with a certain time lag due to viscous damping and no particle rotation can occur .    for any non - zero noise intensity @xmath143 , there will be few stochastic transitions between the deterministic trajectories and their detailed form in the limit @xmath139 can be determined . comparing the associated rates for forward ( increasing @xmath49 by @xmath79 ) and backward ( decreasing @xmath49 by @xmath79 ) phase slips we may directly obtain direction and magnitude of the noise - driven transport by means of ( [ eq24 ] ) .
the analysis of this section builds on the path integral representation for the transition probability of a markovian stochastic process @xcite .
the main complication is the _ time dependence _ of the potential @xmath104 . in notation and general strategy
we follow the recent detailed analysis of noise - driven escape over oscillating barriers @xcite
. we will only determine the dominating exponential term in the transition probabilities . a more detailed analysis including also the prefactor can be done along the lines of ref .
@xcite .
the transition probability @xmath144 of the markov process described by ( [ lan3 ] ) can be written in the form @xmath145}{d } ) \ , \ ] ] where the action functional @xmath146 $ ] is given by @xmath147=\frac 1 4 \int_{t_i}^{t_f } dt         \big({\partial}_t{\varphi}(t)-f({\varphi}(t),t)\big)^2.\ ] ] with @xmath148 defined in ( [ deff ] ) .
the intuitive understanding behind this representation is that _ all _ possible trajectories @xmath55 starting at initial time @xmath93 in @xmath149 and arriving at final time @xmath95 in @xmath150 contribute to the transition probability @xmath144 , each with a weight related to the value of the action @xmath151 evaluated along this trajectory . in the weak noise limit @xmath139 a laplace argument tells us that only those transitions @xmath152 contribute significantly to @xmath144 for which the action is minimal .
all others have probabilities which are exponentially small for @xmath139 . in the weak noise limit
we are hence able to determine the precise form of the ( dominating ) stochastic transitions by solving a variational problem .    under rather general conditions which are satisfied in our particular setting one can show that the transition rate between two stable orbits of the dynamics is determined by the probability to reach from the starting stable orbit the unstable orbit separating the two stable ones @xcite . the subsequent relaxation to the final stable orbit is then purely deterministic and does not contribute to the transition probability @xmath144 in the weak - noise limit .
we have hence to solve a variational problem for which the initial and final point of the trajectory are not fixed .
instead they are required to lie on two known functions @xmath153 and @xmath154 , respectively . if we perform the optimization with respect to these locations we find @xmath155 hence the optimal transition trajectory @xmath152 must be _ parallel _ to the stable and unstable orbit at the contact points .
this in turn implies @xmath156 and @xmath157 .
the optimal transition trajectory @xmath152 hence starts out at @xmath156 at the stable orbit , moves for a very long time near to it , changes then in a rather short time interval to the immediate neighbourhood of the unstable orbit ( hence the name `` instanton '' ) which it finally reaches at @xmath157 @xcite .    from the minimization of the action
@xmath146 $ ] with respect to the functions @xmath55 we get the usual euler - lagrange equation for the minimizing instanton @xmath152 which for the action given by ( [ action ] ) takes the form @xmath158    choosing for the time dependence of @xmath6 the example specified by ( [ foft1 ] ) we have solved ( [ ele1d ] ) numerically as a boundary value problem for a system of ordinary differential equations .
the initial and final points were chosen slightly off the stable and unstable orbit , respectively .
if the difference between the initial and final time , @xmath159 , is larger than several periods of the external driving , the solution hardly depends on the precise value of these deviations from the deterministic orbits .
having obtained @xmath152 , the corresponding value of the action @xmath151 can be calculated from ( [ action ] ) .    in fig .
[ fig : inst1d ] the resulting instanton trajectories @xmath152 corresponding to backward and forward transitions are displayed .
they show the qualitative behaviour discussed above .
note , however , that only the noise induced transition from the stable to the unstable orbit is shown .
the complete transitions giving rise to @xmath160 comprises also the subsequent relaxations from the unstable orbit to the next stable one .
consequently the stochastic transitions show substantially more structure than the rather poor caricatures used to represent them in fig .
[ fig : fig1 ] .
the corresponding values of the action are @xmath161 for the forward instanton increasing @xmath49 by @xmath79 and @xmath162 for the backward instanton decreasing @xmath49 by @xmath79 .
hence , at least in the weak noise limit , forward transitions are more likely ( their action is smaller ) than backward transitions and correspondingly for the parameter set specified in the caption of fig .  [
fig : inst1d ] we find a noise driven rotation in the direction of increasing @xmath49 .
the mismatch between forward and backward transition clearly demonstrates the operation of the ratchet effect in our system . with the relevant transition probabilities being of order
@xmath163 the effect is however very weak and in particular rapidly disappears with @xmath139 .
numericaly determined instanton trajectories ( full lines ) describing the dominating weak noise stochastic transitions between stable ( dashed lines ) and unstable ( dashed - dotted lines ) deterministic periodic orbits of the one - dimensional rotational dynamics of a ferrofluid particle described by ( [ lan3 ] ) , ( [ deff ] ) , ( [ foft1 ] )
. the parameter values are @xmath164 and @xmath68 . for the forward instanton ,
increasing @xmath49 by @xmath79 ( right ) , the value of the action is @xmath161 , for the backward instanton , decreasing @xmath49 by @xmath79 ( left ) , it is @xmath162 .
hence , a noise induced average net rotation in the `` forward direction '' results in the weak noise limit , implying @xmath165 according to ( [ eq24]).,title="fig : " ] 0.5 cm   numericaly determined instanton trajectories ( full lines ) describing the dominating weak noise stochastic transitions between stable ( dashed lines ) and unstable ( dashed - dotted lines ) deterministic periodic orbits of the one - dimensional rotational dynamics of a ferrofluid particle described by ( [ lan3 ] ) , ( [ deff ] ) , ( [ foft1 ] )
. the parameter values are @xmath164 and @xmath68 . for the forward instanton ,
increasing @xmath49 by @xmath79 ( right ) , the value of the action is @xmath161 , for the backward instanton , decreasing @xmath49 by @xmath79 ( left ) , it is @xmath162 .
hence , a noise induced average net rotation in the `` forward direction '' results in the weak noise limit , implying @xmath165 according to ( [ eq24]).,title="fig : " ]      for the numerical solution of the fpe describing the complete two - dimensional dynamics for general values of the parameters it is convenient to expand the probability distribution @xmath136 in spherical harmonics @xmath166 @xcite @xmath167 with so far unknown time dependent complex expansion coefficients @xmath168 .
since @xmath136 is real we have @xmath169 with the star denoting complex conjugation .
moreover , from the normalization of @xmath136 we find @xmath170 plugging the ansatz ( [ pexp ] ) into the fpe ( [ fpe ] ) , exploiting ( [ defpot ] ) , and using well - known properties of the spherical harmonics we find for @xmath171 the following set of ordinary differential equations @xmath172 where we have introduced the auxiliary factors @xmath173 moreover , we find for the ensemble averages of the components of @xmath72 from ( [ defe ] ) @xmath174 where @xmath175 and @xmath176 denote the real and imaginary part of the complex number @xmath59 , respectively .
for the @xmath59-component of the ensemble averaged but still time - dependent magnetic torque @xmath177 this implies via ( [ defmagtor ] ) @xmath178    we have used two somewhat different procedures to explicitly calculate this torque .
the first is more general , the second is slightly more efficient for the special case of the time dependence ( [ foft1 ] ) .    in the first method we solve the system ( [ odeset ] ) of ordinary differential equations for the expansion coefficients
@xmath168 with @xmath179 by a standard routine .
we have used values of @xmath180 up to 25 but found typically @xmath181 to be sufficient . starting with some initial condition
, we integrate the equation until the periodic time dependence of the solution is reached , which was usually the case after at most two periods of the external field .
we then integrate further for one period of the external driving , use the result for @xmath182 in ( [ torexp ] ) , and calculate the remaining time average numerically .
some care is needed here since this time - averaged torque @xmath61 is typically 2 - 3 orders of magnitude smaller than the typical time - dependent values of @xmath183 . in the time average
we are hence subtracting numbers of equal size from each other giving rise to well - known problems of numerical accuracy .
the advantage of this method is that it works for any numerically sensible time dependence of the magnetic field @xmath6 in ( [ deffield ] ) .
the second method builds on the fact that the long - time solution of the fpe will be periodic in time with the period @xmath79 of the external driving .
it is hence useful to expand the time - dependent expansion coefficients @xmath168 in ( [ pexp ] ) into fourier modes with respect to time @xmath184 instead of a system of ode s for the coefficients @xmath168 we find now from the fpe ( [ fpe ] ) a system of _ algebraic _
equations for the coefficients @xmath185 . moreover , the time averaged torque @xmath61 can be readily expressed as function of some of these coefficients by means of ( [ torexp ] ) .
the advantages of this method are that effectively we are directly dealing with the stationary solution of the fpe , _
i.e. _ no initial conditions are necessary and no equilibration process must be simulated , and that the time average of the torque needs not to be performed numerically .
the disadvantage is that we have to specify @xmath6 in order to find the algebraic system for the @xmath185 and that for general @xmath6 this system will be rather dense with respect to the index @xmath186 .
for the special time dependence ( [ foft1 ] ) the situation is somewhat more gratifying since @xmath6 involves only two fourier components .
hence there are only couplings between coefficients @xmath185 differing in their @xmath186-value by at most 2 and the result for @xmath61 is a linear combination of just 5 different coefficients @xmath185 .
this makes the numerical analysis rather fast .
both methods must of course yield the same results when applied to the same @xmath6 .
we have frequently used both to verify our numerical findings .    using the numerical methods described above we have verified all the predictions derived in section [ sec : symm ] on the basis of symmetry arguments .
in particular we find @xmath127 if @xmath187 in accordance with ( [ 4a ] ) and @xmath127 for @xmath188 , _
i.e. _ for ( [ foft1 ] ) with @xmath189 , or @xmath190 in accordance with ( [ 4b ] ) .    in fig .
[ fig : numres ] we give some examples for current inversions in our system .
the left panel shows the averaged torque @xmath61 as function of the phase angle @xmath13 for a time dependent magnetic field @xmath6 of the form ( [ foft1 ] ) .
the results are consistent with ( [ nofdel ] ) .
moreover , the values of @xmath13 at which the torque changes sign are clearly different from @xmath75 and @xmath191 indicated by the squares on the horizontal axis .
hence the numerical solution verifies the statement of section  [ sec : symm ] that the current reversal in the full two dimensional model ( [ lan1 ] ) , ( [ lan2 ] ) occurs at values of @xmath13 which are not fixed by any symmetry but rather depend in a complicated manner on all the other parameters of the problem .
this is further corroborated by the right panel of fig .
[ fig : numres ] , depicting the torque @xmath61 as function of the relative particle size for @xmath128 .
there are two different ways to scale the particle size that make sense in a ferrofluid . in the first ,
both the ferromagnetic core of the particle as well as the polymer coating are scaled with the same factor .
this gives rise to a proportional change of the magnetic moment @xmath17 and the particle volume @xmath0 . in the dimensionless units adopted
this implies unchanged fields @xmath4 and @xmath6 and a noise intensity @xmath46 rescaled according to ( [ defd ] ) . in the other case ,
only the polymer coating is scaled and therefore only the hydrodynamic radius of the particle changes .
hence the magnetic moment remains the same whereas the volume @xmath0 changes .
correspondingly the dimensionless fields @xmath4 and @xmath6 and the noise intensity @xmath46 change by the same factor . in view of ( [ fpe ] ) this can be absorbed in a rescaled time which in turn is equivalent to a changed frequency @xmath192 of the external field . as anticipated in section  [ sec : symm ] , we find current inversions in both cases .     left : averaged torque @xmath61 for the time modulation ( [ foft1 ] ) with @xmath164 and @xmath193 as function of the phase angle @xmath13 as obtained from the numerical solution of the fpe ( full line ) .
symbols are results from simulations of the langevin equations ( [ lan1 ] ) and ( [ lan2 ] ) as reported in @xcite with statistical errors being smaller than the symbol size .
the squares on the horizontal axis are the points @xmath194 and @xmath195 at which the current reversal would occur in the reduced model ( [ lan3 ] ) .
right : dependence of @xmath61 on the relative particle size for the same values of the parameters and @xmath13 fixed to @xmath196 , the first point of current inversion in the left panel .
the full line corresponds to a proportional change of the magnetic and hydrodynamic radii ( top sketch ) , the dashed line is for a change of the hydrodynamic radius only ( bottom sketch).,title="fig : " ] 0.5 cm   left : averaged torque @xmath61 for the time modulation ( [ foft1 ] ) with @xmath164 and @xmath193 as function of the phase angle @xmath13 as obtained from the numerical solution of the fpe ( full line ) .
symbols are results from simulations of the langevin equations ( [ lan1 ] ) and ( [ lan2 ] ) as reported in @xcite with statistical errors being smaller than the symbol size .
the squares on the horizontal axis are the points @xmath194 and @xmath195 at which the current reversal would occur in the reduced model ( [ lan3 ] ) .
right : dependence of @xmath61 on the relative particle size for the same values of the parameters and @xmath13 fixed to @xmath196 , the first point of current inversion in the left panel .
the full line corresponds to a proportional change of the magnetic and hydrodynamic radii ( top sketch ) , the dashed line is for a change of the hydrodynamic radius only ( bottom sketch).,title="fig : " ]      it is very helpful to have some analytical expression for the average torque @xmath61 since the dependencies on the various parameters of the system can then be identified much more directly . in the present and subsequent subsections we discuss two approximate methods to calculate the time averaged torque in our system .
the first one employs the so - called effective field method to approximately solve the fpe .
it has been a standard tool in the theory of ferrofluids for many years @xcite .
two main ingredients are necessary .
first , let us recall that the stationary ( equilibrium ) solution of the fpe ( [ fpe ] ) for a _ time independent _ homogeneous magnetic field @xmath3 is of the form @xmath197 where @xmath198 denotes the modulus of the magnetic field , @xmath199 .
averages with this distribution are easily calculated , in particular we find the well - known result @xmath200 with the langevin function @xmath201 second , for general , time - dependent fields @xmath73 the fpe ( [ fpe ] ) can be used to derive the _ exact _
equation @xmath202 for the time evolution of the average @xmath203 .
note that this equation for the first moment of @xmath204 has the usual flaw of being not closed but involving higher moments of the probability distribution .
the central idea of the effective field approximation is now to assume that the solution @xmath205 of the fpe for a general time dependent magnetic field @xmath73 can be written as equilibrium distribution for some a - priori unknown _ effective field _
@xmath206 , _ i.e. _ @xmath207 . the average on the r.h.s . of ( [ exequ ] ) can then be performed and this equation gives rise to an evolution equation for the effective field @xmath208 of the form @xcite @xmath209 given @xmath73 , this is a closed equation for the time evolution of @xmath206 .
having obtained @xmath206 we get the desired approximate result for @xmath203 from ( [ enot ] ) @xmath210 the solution of the partial differential equation ( [ fpe ] ) is hence replaced by the solution of the set of three coupled ordinary differential equations ( [ efff ] ) .
there are other , equivalent forms of ( [ efff ] ) , see _
e.g. _ @xcite .
exploiting ( [ efff ] ) , it is also possible to derive the following equivalent closed equation for @xmath203 @xcite @xmath211 where @xmath212 denotes the modulus of @xmath213 and the effective field @xmath208 is to be expressed as function of @xmath213 by the inverse of ( [ eefff ] ) .
note the difference between @xmath214 and @xmath215 .    in general , ( [ efff ] ) can only be solved numerically . in fig .
[ fig : efff ] we give a comparison between results obtained in this way for the time dependence ( [ foft1 ] ) and the corresponding outcome of the numerical solution of the fpe .
it is clearly seen that the effective field approximation yields rather accurate results for the time dependent orientation of the particle and , equivalently , for the time dependent magnetization of the ferrofluid . on the other hand the results for the time averaged torque @xmath61
differ from the numerically exact values by a factor between 2 and 3 .
the reason for this discrepancy lies in the fact that @xmath61 is typically 2 - 3 orders of magnitude smaller than typical values of the time dependent torque @xmath177 . in calculating the time average of the torque we hence subtract quantities of comparable order and
so amplify the inaccuracies of the effective field approximation from a few percent to a few hundred percent .
horizontal components of the average orientation @xmath203 of the particle as function of time for one period of the external magnetic field .
the lines show the results for the @xmath5-and @xmath8-component of @xmath203 ( full and dashed line , respectively ) from the numerical solution of the fokker - planck equation ( [ fpe ] ) with a time dependence of the external field as given in ( [ foft1 ] ) .
the symbols are the corresponding results from the effective field theory .
the parameter values are @xmath216 as well as @xmath217 ( left ) and @xmath193 ( right ) .
the corresponding values for the time averaged torque @xmath61 are also displayed.,title="fig : " ] 0.5 cm   horizontal components of the average orientation @xmath203 of the particle as function of time for one period of the external magnetic field .
the lines show the results for the @xmath5-and @xmath8-component of @xmath203 ( full and dashed line , respectively ) from the numerical solution of the fokker - planck equation ( [ fpe ] ) with a time dependence of the external field as given in ( [ foft1 ] ) .
the symbols are the corresponding results from the effective field theory .
the parameter values are @xmath216 as well as @xmath217 ( left ) and @xmath193 ( right ) .
the corresponding values for the time averaged torque @xmath61 are also displayed.,title="fig : " ]    in order to obtain an analytical expression for @xmath61 within the framework of the effective field approximation , it is useful to consider the case of small values of @xmath218 in which the langevin function @xmath219 or equivalently its inverse may be approximated by the first terms of the respective power series expansion .
using the lowest order approximation @xmath220 we have @xmath221 and from ( [ efff1 ] ) it then follows that @xmath222 for large @xmath15 we therefore find @xmath223 hence @xmath224 tends asymptotically to a constant whereas @xmath225 contains only oscillating parts . from ( [ defmagtor ] )
we may thus infer that the torque @xmath226 harmonically oscillates and therefore its time average vanishes , @xmath127 .    to get a non - trivial result for @xmath61 we have to push the expansion of the langevin function further .
in other words , the non - linearity of the magnetization curve is essential to get a non - zero average torque .
if the next term in the expansion is included , @xmath227 , we already obtain an approximation for @xmath219 that is accurate to within one per cent for the experimentally relevant values of @xmath5 which are generally less than @xmath228 . calculating the corrections to ( [ efffpert1 ] ) and ( [ efffpert2 ] ) induced by these higher order terms we finally find after some algebra the following result for the average torque @xmath229 this result is in accordance with our symmetry considerations of section [ sec : symm ] since it shows that both a non - zero field component in @xmath5-direction and an even higher harmonic in the time dependence of @xmath6 are essential for a non - zero torque to occur , cf .
( [ 4a ] ) , ( [ 4b ] ) .      in the last subsection we saw that the effective field approximation alone is not sufficient to get an explicit analytical expression for the time averaged torque @xmath61 in our system
. additionally we had to expand the langevin function in its argument @xmath218 which amounts to a small field or equivalently large noise expansion . moreover , the approximate result for the torque was typically at variance with the numerically exact result by a factor of about 2 .
it is therefore tempting to avoid the effective field approximation altogether and to use a perturbation expansion in the ratio of magnetic field strength and noise intensity right in the fokker - planck equation ( [ fpe ] ) .
this is conveniently accomplished by making the formal substitution @xmath230 in ( [ fpe ] ) and by expanding the solution of the fpe in a power series in @xmath231 @xmath232 at the end of the calculation @xmath231 is put equal to one .
it will turn out that we have to use the expansion up to third order in @xmath231 to get a non - zero result for the average torque @xmath61 . since the potential @xmath104 is of order @xmath231
we find in zeroth order @xmath233 describing the stationary distribution of pure rotational diffusion .
using the expansion ( [ ppert ] ) in ( [ fpe ] ) we get a system of linear inhomogeneous partial differential equations for the various @xmath234 that can be solved in principle one after the other .
this procedure is , however , rather cumbersome
. it can be somewhat simplified by using the exact equation ( [ exequ ] ) , the r.h.s of which is of order @xmath231 .
hence to find @xmath203 from this equation to order @xmath235 it is sufficient to calculate the average on the r.h.s . with an expression for @xmath136 accurate to order @xmath236 . in our case
the lowest order non - zero result for @xmath61 is @xmath237 , hence we need @xmath203 to order @xmath238 ( cf .
( [ defmagtor ] ) , ( [ defe ] ) ) and therefore the second order result for @xmath136 would be sufficient .
it is , however , more convenient to use again the expansion ( [ pexp ] ) of @xmath136 in terms of spherical harmonics .
the perturbation ansatz ( [ ppert ] ) then translates into an expansion of the coefficients @xmath168 in powers of @xmath231 @xmath239 the ode s ( [ odeset ] ) couple only coefficients @xmath168 which differ in both @xmath240 and @xmath241 by 1 .
moreover , we know from ( [ resa00 ] ) @xmath242 and hence @xmath243 .
as shown by ( [ torexp ] ) the torque is determined by @xmath182 alone .
we hence need an expression for @xmath182 correct up to order @xmath238 .    to this end
we determine from ( [ odeset ] ) and ( [ resa00 ] ) first the @xmath244 terms of @xmath245 . with these
we calculate the necessary @xmath246 to order @xmath247 which in turn are used to determine the @xmath248 terms in @xmath182 .
the procedure is sketched in fig .
[ fig : pertdiag ] . note that it is sufficient to calculate coefficients @xmath168 with @xmath249 because of the symmetry property ( [ asym ] ) .
diagrammatic sketch of the perturbative determination of the expansion coefficients @xmath250 in ( [ pexp ] ) . as shown by the arrows only @xmath250 differing in both @xmath240 and @xmath241 by 1
are coupled by equations ( [ odeset ] ) . in a first step
we start from @xmath251 as given in ( [ resa00 ] ) and obtain the @xmath244 terms of @xmath252 ( dashed arrows ) .
next @xmath253 and @xmath254 are determined to order @xmath247 ( dotted arrows ) .
finally these are used to get the @xmath248 contributions to @xmath255 ( full arrows ) from which the torque can be extracted to order @xmath256 using ( [ torexp ] ) .
all other @xmath250 are irrelevant for @xmath61 at this order in @xmath231 .
note also that the left half of the diagram is redundant due to the symmetry property ( [ asym]).,height=241 ]    in this way we first get from the @xmath244 equations @xmath257 since we are interested in the asymptotic solution valid for large @xmath15 we can put @xmath258 .
next we use the @xmath259 equations @xmath260 to determine the relevant @xmath246 to the desired order .
finally we obtain @xmath261 from @xmath262 the explicit expressions for the various coefficients @xmath263 are rather long , the above system of equations is therefore conveniently solved with the help of a computer algebra . focusing on the specific time dependent field @xmath6 from ( [ foft1 ] ) , exploiting ( [ torexp ] ) and performing the time average over one period of the external driving we finally get the following lowest order perturbative result for the average torque @xmath264 the dependencies on the parameters @xmath265 , and @xmath266 are identical to those of the effective field result ( [ effresn ] ) , whereas the dependence on the noise strength @xmath46 is more complicated .
one readily verifies that the symmetry properties ( [ 3a])-([3d ] ) are satisfied by ( [ pertresn ] ) . in particular , these symmetries imply that the torque @xmath61 must to be an odd function of both the static field and the oscillating field .
moreover , symmetry reasons also imply @xcite that in linear order of the oscillating field one can not expect a finite @xmath61 .
rather , a coupling of several modes of the oscillating field ( harmonic mixing ) is required , see also @xcite .
all together , these arguments explain that the first nontrivial contribution must be ( at least ) of fourth order in the magnetic field strength and that within a linear response theory one will always find @xmath267 @xcite . on the other hand
this does of course not imply that the non - linear magnetization curve alone would `` explain '' the effect . to understand the microscopic origin of the angular momentum transfer from the external field to the particle one really has to go to the description in terms of the fpe as discussed above , see also @xcite .
left : time averaged torque @xmath61 as a function of the noise intensity @xmath46 for the time dependence ( [ foft1 ] ) of the magnetic field .
the parameter values are @xmath268 , and @xmath68 .
the scaling of the magnetic field strength with @xmath46 ensures that the perturbation parameter @xmath231 stays constant .
shown are the perturbative results ( [ pertresn ] ) ( circles ) and ( [ effresn ] ) ( squares ) together with the numerically exact result ( full line ) .
right : same for the time dependence ( [ foft2 ] ) of the oscillating field with @xmath269 .
the approximate results are given by ( [ pertresn2 ] ) and ( [ effresn2 ] ) , respectively.,title="fig : " ]   left : time averaged torque @xmath61 as a function of the noise intensity @xmath46 for the time dependence ( [ foft1 ] ) of the magnetic field .
the parameter values are @xmath268 , and @xmath68 .
the scaling of the magnetic field strength with @xmath46 ensures that the perturbation parameter @xmath231 stays constant .
shown are the perturbative results ( [ pertresn ] ) ( circles ) and ( [ effresn ] ) ( squares ) together with the numerically exact result ( full line ) .
right : same for the time dependence ( [ foft2 ] ) of the oscillating field with @xmath269 .
the approximate results are given by ( [ pertresn2 ] ) and ( [ effresn2 ] ) , respectively.,title="fig : " ]    in the left panel of fig .
[ fig : perttest ] we compare the two approximate expressions ( [ effresn ] ) and ( [ pertresn ] ) for the averaged torque with the results obtained from the numerical solution of the fpe .
displayed is the dependence of the torque on the noise intensity @xmath46 which is different in the two approximations . in order to keep the expansion parameter @xmath270 roughly constant
when changing @xmath46 we have scaled the magnetic fields strengths with @xmath46 . as can be seen , up to values of @xmath271
the accuracy of ( [ pertresn ] ) is rather good whereas ( [ effresn ] ) describes the torque only qualitatively .
proportional change of magnetic field strength and noise intensity is equivalent to a change of the frequency @xmath192 of the external field in ( [ deffield ] ) .
[ fig : perttest ] hence also demonstrates the resonance - like character of the investigated ratchet effect .
@xmath61 is largest when the deterministic time scale of the external driving matches an intrinsic stochastic time scale of the system related to the brownian relaxation time @xmath272 defined in ( [ deftaub ] ) .
the two approximate expressions for the torque ( [ effresn ] ) and ( [ pertresn ] ) differ in their dependence on the phase angle @xmath13 . when using ( [ pertresn ] ) to fit the experimental results reported in @xcite we find for the brownian relaxation time the value @xmath273s instead of @xmath274s as obtained on the basis of ( [ effresn ] ) in @xcite .
the experiments were done with a ferrofluid with @xmath275 pas . using ( [ deftaub ] ) the two results for @xmath272
hence translate in fits for the particle diameter of @xmath276 nm and @xmath277 nm , respectively .
these values exceed the typical diameter of roughly @xmath278 nm by factors between 2 and 3 , respectively .
the main reason for this discrepancy is probably the polydispersedness of real ferrofluids , having a particle size distribution with a long tail ( see _ e.g. _ @xcite , ch.2 ) .
hence , a whole spectrum of relaxation times is necessary to accurately describe the dynamics of the magnetization .
however , this issue is not at the focus of the present investigation .
we also note that the transferred torque typically increases with the particle size and hence the described ratchet effect is likely to be dominated by the larger grains in the ferrofluid .
a similar approximate calculation of the average torque is also possible for the time dependence ( [ foft2 ] ) .
the result of the effective field approximation is @xmath279 whereas we find from the perturbative solution of the fpe @xmath280 the two expression only differ significantly from each other if the noise intensity @xmath46 is very small . in the right panel of fig .
[ fig : perttest ] they are compared with the result from the numerical solution of the fpe ( [ fpe ] ) .
for @xmath271 the accuracy is again seen to be rather satisfactory
in the present paper we have theoretically investigated the rotational brownian motion of colloidal ferromagnetic particles in an oscillating magnetic field .
the central tool was the fokker - planck equation for the probability density @xmath204 of the particle orientation @xmath18 at time @xmath15 . solving this equation either numerically or approximately by using the effective field approximation as well as a perturbative expansion we have determined the time averaged torque @xmath61 exerted by the magnetic field on the particles .
the main and a priori quite unexpected qualitative finding is the fact that a purely oscillating magnetic field without net rotating component can transfer angular momentum to a ferromagnetic grain . as basic mechanism behind this transfer ,
a ratchet effect was identified by which the magnetic field rectifies the thermal fluctuations of the particle orientation that arise due to random collisions with the molecules of the carrier liquid .
the detailed operation of this ratchet effect in the present system was discussed on the basis of a weak - noise analysis of the fokker - planck equation for a closely related one - dimensional system in section [ inst1d ] .    via the viscous coupling to the carrier liquid the torque on the particle is transmitted to the liquid .
the combined _ microscopic _ torques of the huge number of individual nanoparticles then yields a _ macroscopic _ rotation of the ferrofluid as a whole , as observed experimentally in @xcite . in the absence of thermal noise ,
no average torque can arise , _
i.e. _ thermal fluctuations are an indispensable requirement for the rotation of the individual grains and hence of the ferrofluid as a whole to occur .
the results from the fokker - planck equation are completely consistent with rather general symmetry considerations detailed in section [ sec : symm ] . moreover , they agree very well with simulations of the corresponding langevin equations as given in @xcite and quantitatively describe the experimental findings reported in @xcite .
we remark that our present system puts forward a new type of thermal ratchet device which does not fit into any previously known specific class of ratchet systems . adopting the classification scheme from @xcite ,
our present system has some similarity with so - called asymmetrically tilting ratchets as well as with a so - called traveling potential ratchets , however in the generalized sense involving two counterpropagating traveling potentials .
yet there remain significant differences with both these classes . in particular , we note that we are not dealing here with a periodic , asymmetric so - called `` ratchet''-potential .
rather , at any fixed instance of time the relevant potential is perfectly symmetric about the instantaneous direction of the magnetic field and hence no preferential direction of rotation seems to exists .
it is only via the time evolution that a symmetry breaking arise , which is sometimes called a dynamical symmetry breaking .
in summary we hope that the present investigation has demonstrated that ferrofluids are very suitable systems to study various aspects of thermal ratchet behaviour , and that it may stimulate further theoretical , numerical and experimental work in this direction .
special thanks is due to ubbo felderhoff for several interesting suggestions and to martin raible for providing the numerical simulation results in fig .
discussions with hanns walter mller , konstantin morozov , and mark shliomis are gratefully acknowledged .
part of this work was done during a sabbatical stay of a.e . at the laboratoire de physique thorique at the universit louis pasteur in strasbourg .
he acknowledges the kind hospitality as well as financial support from the _
volkswagenstiftung_. this work was further supported by deutsche forschungsgemeinschaft under sfb 613 and re 1344/3 - 1 and by the esf - program stochdyn
.    99 a. einstein , _ ann .
d. physik _ * 19 * , 371 ( 1906 ) m. i. shliomis , _ sov
* 17 * , 153 ( 1974 ) r. e. rosensweig , _ ferrohydrodynamics _ , ( cambridge university press , cambridge , 1985 ) j .- c .
bacri , a. cebers , and r. perzynski , _ phys .
lett . _ * 72 * , 2705 ( 1994 ) a. v. lebedev , a. engel , and k. i. morozov , _ new .
_ * 5 * , 57 ( 2003 ) j. p. mctague , j. chem .
phys . * 51 * , 133 ( 1969 ) m. i. shliomis and k. i. morozov , _ phys .
fluids _ * 6 * , 2855 ( 1994 ) f. gazeau et al .
, _ europhys .
_ * 35 * , 609 ( 1996 ) h. w. mller , _ phys .
lett . _ * 82 * , 3907 ( 1999 ) a. engel , h. w. mller , p. reimann , and a. jung , _ phys .
lett . _ * 91 * , 060602 ( 2003 ) j. c. maxwell , _ theory of heat _ ( longmans , green , and co. , london , 1871 ) m. v. smoluchowski , _ z. physik _ * 13 * , 1069 ( 1912 ) r. p. feynman , r. b. leighton , and m. sands , _ the feynman lectures on physics - vol.1 _ , ch . 46 ( addison - wesley , reading , 1963 ) s. l. harvey and a. f. rex ( eds . ) , _
maxwell s demon : entropy , information , computing _
( adam hilger , bristol , 1990 ) p. reimann , _ phys .
rep . _ * 361 * , 57 ( 2002 ) r. d. astumian and p. hnggi , _ physics today _ * 55 * , 33 ( 2002 ) m. o. magnasco , _ phys .
lett . _ * 71 * , 1477 ( 1993 ) f. jlicher , a. ajdari , and j. prost , _ rev .
phys . _ * 69 * , 1269 ( 1997 ) j. rousselet , l. salome , a. ajdari , and j. prost , _ nature _ * 370 * , 446 ( 1994 ) h. linke ( ed . )
_ special issue on ratchets and brownian motors : basic experiments and applications _
phys . a _ * 75 * , 167 - 352 ( 2002 ) h. linke et al . , _ science _ * 286 * , 2314 ( 1999 ) i. turcu and c. m. lucaciu , j. phys .
a : math . gen .
* 22 * , 995 ( 1989 ) p.
palffy - muhoray and e. weinan , _ mol .
cryst . _ * 320 * , 193 ( 1998 ) s. flach , a. e. miroshnichenko , and a. a. ovchinikov , _ phys .
rev . _ * e65 * , 104438 ( 2002 ) w. t. coffey , y. p. kalmykov , and j. t. waldron , _ the langevin equation with applications in physics , chemistry and electrical engineering _
, ( world scientific , singapore , 1996 ) f. ritort , _ phys . rev .
lett . _ * 80 * , 6 ( 1998 ) z. wang , c. holm , and h. w. mller , _ phys .
* e66 * , 021405 ( 2002 ) l. d. landau and e. m. lifshitz , _ electrodynamics of continuous media _ , 2nd ed . ,   34 ( pergamon , new york , 1984 ) .
l. d. landau and e. m. lifshitz , _ fluid mechanics _
, 2nd ed . ,   20 , ( pergamon , new york , 1987 ) .
m. raible and a. engel , _ langevin equation for the rotation of a magnetic particle _ , to appear in _ applied organometallic chemistry _
h. risken , _ the fokker - planck equation _
, springer , berlin 1984 . v. becker and a. engel , unpublished p.  curie , j. phys . ( paris ) 3 .
srie ( thorique et appliqu ) * t. iii * , 393 ( 1894 ) .
a.  ajdari , d.  mukamel , l.  peliti , and j.  prost , j. phys .
i france * 4 * , 1551 ( 1994 ) r.  kanada and k.  sasaki , j. phys .
* 68 * , 3759 ( 1999 ) s.  weiss , d.  koelle , j.  mller , k.  barthel , and r.  gross , europhys
. lett . * 51 * , 499 ( 2000 ) .
s.  flach , o.  yevtushenko , and y.  zolotaryuk , phys .
* 84 * , 2358 ( 2000 ) s.  cilla and l.  m. floria , phys .
e * 63 * , 031110 ( 2001 ) e.  neumann and a.  pikovsky , eur . phys .
j. b * 26 * , 219 ( 2002 ) o.  yevtushenko , s.  flach , y.  zolotaryuk , and a.  a.  ovchinikov , europhys
. lett . * 54 * , 141 ( 2001 ) b.  yan , r.  m. miura , and y .- d .
chen , j. theor .
210 * , 141 ( 2001 ) p.  reimann , phys .
lett . * 86 * , 4992 ( 2001 ) s. denisov , s. flach , a. a. ovchinikov , o. yevtushenko , and y. zolotaryuk , phys .
e * 66 * , 041104 ( 2002 ) k. seeger and w. maurer , _ solid state commun . _ * 27 * , 603 ( 1978 ) d. r. chialvo and m. m. millonas , _ phys . lett . _ * a209 * , 26 ( 1995 ) i. goychuk and p. hnggi , _ europhys . lett . _ * 43 * , 503 ( 1998 ) r. p. feynman , a. r. hibbs , _ quantum mechanics and path integrals _ ( mcgraw - hill , new york , 1965 ) l. s. schulman , _ techniques and applications of path integration _ ( wiley , new york , 1981 ) m. chaichian and a. demichev , _ path integrals in physics _ , vol.1 ( iop publishing , bristol , 2001 ) j. lehmann , p. reimann , and p. hnggi , _ phys .
_ * 84 * , 1639 ( 2000 ) ; _ phys .
rev . _ * e62 * , 6282 ( 2000 ) ; _ phys .
( b ) * 237 * , 53 ( 2003 ) and further references therein .
m. abramowitz and i. stegun ( eds . ) , _ handbook of mathematical functions _
( national bureau of standards , washington d.c . , 1964 )
m. a. martsenyuk , yu . l. raikher , and m. i. shliomis , _ sov .
. jetp _ * 38 * , 412 ( 1974 ) h. w. mller and m. liu , _ phys .
rev . _ * e64 * , 061405 ( 2001 ) t.  alarcon , a.  perez - madrid , j.  m. rubi , _ phys .
lett . _ * 85 * , 3995 ( 2000 ) m. shliomis , _ phys .
lett . _ * 92 * , 188901 a. engel and p. reimann , _ phys .
_ * 92 * , 188902