- Fluctuation dissipation theorem
statistical physics, the fluctuation dissipation theorem is a powerful tool for predicting the non-equilibrium behavior of a system — such as the irreversible dissipationof energy into heat— from its reversible fluctuations in thermal equilibrium. The fluctuation dissipation theorem applies both to classical and quantum mechanical systems. Although formulated originally by Nyquistin 1928 [cite journal | author = Nyquist H | date = 1928 | title = Thermal Agitation of Electric Charge in Conductors | journal = Physical Review | volume = 32 | pages = 110–113 | doi = 10.1103/PhysRev.32.110] , the fluctuation-dissipation theorem was first proved by Herbert B. Callen and Theodore A. Welton in 1951. [cite journal | author = Callen HB, Welton TA | date = 1951 | title = Irreversibility and Generalized Noise | journal = Physical Review | volume = 83 | pages = 34–40 | doi = 10.1103/PhysRev.83.34]
The fluctuation dissipation theorem relies on the assumption that the response of a system in thermodynamic equilibrium to a small applied force is the same as its response to a spontaneous fluctuation. Therefore, there is a direct relation between the fluctuation properties of the thermodynamic system and its linear response properties. Often the linear response takes the form of one or more exponential decays.
Example: Brownian motion
For example, Einstein in his 1905 paper on
Brownian motionnoted that the same random forces which cause the erratic motion of a particle in Brownian motion would also cause drag if the particle were pulled through the fluid. In other words, the "fluctuation" of the particle at rest has the same origin as the "dissipative" frictional force one must do work against, if one tries to perturb the system in a particular direction.
From this observation he was able to use
statistical mechanicsto derive a previously unexpected connection, the Einstein-Smoluchowski relation:
linking "D", the diffusion constant, and "μ", the
mobilityof the particles. ("μ" is the ratio of the particle's terminal drift velocity to an applied force, "μ = vd / F"). "k""B" ≈ 1.38065 × 10−23 m² kg s−2 "K"−1 is Boltzmann's constant, and "T" is the absolute temperature.
Example: Thermal noise in a resistor
John B. Johnsondiscovered and Harry Nyquistexplained Johnson–Nyquist noise. With no applied current, the mean-square voltage depends on the resistance "R", , and the bandwidth over which the voltage is measured:
The examples above are consequences of the "fluctuation dissipation theorem", a very general result of
statistical thermodynamicswhich quantifies the relation between the fluctuations in a system at thermal equilibriumand the response of the system to applied perturbations. It thus allows, for example, the use of molecular models to predict material properties in the context of linear response theory. The theorem assumes that applied perturbations (mechanical forces, electric fields, "etc.") are weak enough that rates of relaxation remain unchanged.
General form of the fluctuation dissipation theorem
The fluctuation-dissipation theorem can be formulated in many ways;one particularly useful form is the following:
Let "x" be an observable of a dynamical system with
Hamiltonian"H""0"("x") subject to thermal fluctuations.The observable "x" will fluctuate around its mean value with fluctuations characterized by a power spectrum.Suppose that we can switch on a (scalar) field "f" which alters the Hamiltionianto .The response of the observable "x" to a field "f"("t") changing with time is characterized (to first order) by the susceptibilityor linear response functionof the system
where the perturbation is adiabatically switched on at .
Now the fluctuation dissipation theorem relates the power spectrum to the imaginary part of the
Fourier transformof the susceptibility
The left-hand side describes fluctuations in x, the right-hand side is closely related to the energy dissipated by the system when pumped by an oscillatory field .
(This is already the classical form of the theorem; quantum fluctuations are taken into account byreplacing the prefactor by . A proof can be found by means of the
LSZ reduction, an identity from quantum field theory.)
The fluctuation-dissipation theorem can be generalized in a straight-forward way to the case of space-dependent fields,to the case of several variables or to a quantum-mechanics setting.
Violations of FDT in glassy systems
While the FDT provides a general relation between the response of equilibrium systems to small external perturbations and their spontaneous fluctuations, no general relation is known for systems out of equilibrium.
In the mid 1990s, in the study of non-equilibrium dynamics of spin glass models it was discovered a generalization of FDT valid for asymptotic non-stationary states, where the temperature appearing in the equilibrium relation is substituted by an effective temperature with a non-trivial dependence on the time scales. This relation is proposed to hold in glassy systems beyond the models for which it was initially found.
We derive the fluctuation dissipation theorem in the form given above, using the same notation.Consider the following test case: The field "f" has been on for infinite time and is switched offat "t"=0
We can express the expectation value of "x" by the probability distribution"W"("x",0) and the transition probability
The probability distribution function "W"("x",0) is an equilibrium distribution and hencegiven by the
Boltzmann distributionfor the Hamiltionian
For a weak field , we can expand the right-hand side
here is the equilibrium distribution in the absence of a field.Plugging this approximation in the formula for yields
where "A"("t") is the auto-correlation function of "x" in the absence of a field.
Note that in the absence of a field the system is invariant under time-shifts.We can rewrite using the susceptibilityof the system and hence find with the above equation (*)
stationary processes, the Wiener-Khinchin theoremstates, thatthe power spectrum equals twice the Fourier transformof the auto-correlationfunction
The last step is to Fourier transform equation (**) and to take the imaginary part. For this it is useful to recall that the Fourier transformof a real symmetric function is real, while the Fourier transform of a realantisymmetric function is purely imaginary.We can split into a symmetric and ananti-symmetric part
Now the fluctuation dissipation theorem follows.
The following general derivation of the fluctuation-dissipation theorem uses averaging in
phase space. The derivation applies equally well to classical as well as quantum mechanical systems, although the former uses a continuous integral over phase space, whereas the latter uses a sum over quantum states. To represent both, we introduce the "trace notation", which applies both to classical and quantum systems
where "dΓ" represents an infinitesimal volume in
phase space. Thus, if a system is described by a probability distribution "f"("q", "p") in phase space, the average value of an arbitrary function "A" of the system's state is given by
where angular brackets are used to denote the averaging over the ensemble. In particular, if the probability distribution is given by the equilibrium
Boltzmann distribution, the ensemble average equals
where β = 1/"kBT", "k""B" is the
Boltzmann constantand "T" is the temperaturein Kelvin.
Having defined our notation and basic variables, we now derive the fluctuation-dissipation theorem. Consider a system that has reached equilibrium under the Hamiltonian "H + h", where "h" is much smaller than the thermal energy "kBT". Being in thermal equilibrium, the probability of any state is proportional to its Boltzmann factor "e−β(H + h)". At time "t" = 0, let the perturbation "h" be turned off; given
ergodicity, the system will gradually relax to a new equilibrium, which has Boltzmann factors "e−βH". The fluctuation-dissipation theorem addresses the question of how quickly the system reaches its new equilibrium.
Let the correlation function "c"("t") be defined
which may be written as
where "δA"("t"; "q", "p") is the deviation from its mean at a time "t", given that the system began at time "t" = 0 at position ("q", "p") in phase space. In other words, the integration is over all "initial" positions of the system in phase space.
The mean value of "A" as it evolves towards its new equilibrium is given by
Since "h" is much smaller than the thermal energy "kBT", we may expand the numerator
We may likewise expand the denominator
where we have used
which is much less than one, by our assumption that "h" is much smaller than the thermal energy 1/β = "kBT".
Combining the numerator and denominator, dropping quadratic and high-order terms in <β"h">, and using the indifference of equilibrium to time, we obtain
Let the perturbation "h = −gA" be proportional to the variable "A" with a constant "−g". Then this formula becomes
Note that the system's relaxation is independent of "A" and linear in "g". These results imply that perturbations will relax independently of one another; if two perturbations, "g"1 and "g"2 are applied, the net relaxation will be the sum of the individual relaxations to "g"1 and "g"2 taken separately. Such continuous linear systems have been well-studied, and many methods developed for their solution, such as
Fourier transforms and Laplace transforms.
Onsager reciprocal relations
* H. B. Callen and T. A. Welton, Phys. Rev. 83, 34 (1951)
* L. D. Landau et E. M. Lifshitz, Cours de physique théorique t.5 Physique Statistique (Mir)
* "Fluctuation-Dissipation: Response Theory in Statistical Physics" by Umberto Marini Bettolo Marconi, Andrea Puglisi, Lamberto Rondoni, Angelo Vulpiani, [http://arxiv.org/abs/0803.0719]
* [http://physics416.blogspot.com/2005/12/lecture-24-fluctuation-dissipation.html Audio recording] of a lecture by Prof. E. W. Carlson of
* [http://query.nytimes.com/gst/fullpage.html?res=9F0CE7DB163CF934A15756C0A965958260 "Herbert B. Callen, 73, Theoretical Physicist", "The New York Times" obituary published on]
May 27 1993
Wikimedia Foundation. 2010.