- Fluctuation dissipation theorem
In

statistical physics , the**fluctuation dissipation theorem**is a powerful tool for predicting the non-equilibrium behavior of a system — such as the irreversibledissipation of energy intoheat — from its reversible fluctuations in thermal equilibrium. The fluctuation dissipation theorem applies both to classical and quantum mechanical systems. Although formulated originally byNyquist in 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 motion noted 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 mechanics to derive a previously unexpected connection, the Einstein-Smoluchowski relation:: $D\; =\; \{mu\; ,\; k\_BT\}$

linking "D", the diffusion constant, and "μ", the

mobility of the particles. ("μ" is the ratio of the particle's terminal drift velocity to an applied force, "μ = v_{d}/ F"). "k"_{"B"}≈ 1.38065 × 10^{−23}m² kg s^{−2}"K"^{−1}isBoltzmann's constant , and "T" is theabsolute temperature .**Example: Thermal noise in a resistor**In 1928,

John B. Johnson discovered andHarry Nyquist explainedJohnson–Nyquist noise . With no applied current, the mean-square voltage depends on the resistance "R", $k\_BT$, and the bandwidth $Delta\; u$ over which the voltage is measured:: $langle\; V^2\; angle\; =\; 4Rk\_BT,Delta\; u.$

**General applicability**The examples above are consequences of the "fluctuation dissipation theorem", a very general result of

statistical thermodynamics which quantifies the relation between the fluctuations in a system atthermal equilibrium and 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 $langle\; x\; angle\_0$with fluctuations characterized by apower spectrum $S\_x(omega)$.Suppose that we can switch on a (scalar) field "f" which alters the Hamiltionianto $H(x)=H\_0(x)+x\; f$.The response of the observable "x" to a field "f"("t") changing with time is characterized (to first order) by thesusceptibility orlinear response function $chi(t)$ of the system: $langle\; x(t)\; angle\; =\; langle\; x\; angle\_0\; +\; int\_\{-infty\}^\{t\}\; chi(\; au)\; f(t-\; au)\; d\; au,,$

where the perturbation is adiabatically switched on at $au\; =-infty$.

Now the fluctuation dissipation theorem relates the power spectrum to the imaginary part of the

Fourier transform $ilde\{chi\}(omega)=\; m\{Re\}\; [\{\; ilde\{chi(omega\; )]\; +\; i~\; m\{Im\}\; [\{\; ilde\{chi(omega\; )]$ of the susceptibility $chi(t)$: $S\_x(omega)\; =\; frac\{2\; k\_B\; T\}\{omega\}\; m\{Im\}\; [\; ilde\{chi\}(omega)]$.

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 $f(t)\; =\; f\_0\; sin(omega\; t).$.

(This is already the classical form of the theorem; quantum fluctuations are taken into account byreplacing the prefactor by $\{hbar\}cdot\; \{\; m\; cotanh\}frac\{etahbar\; omega\}\{2\},.$. 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.

**Derivation I**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

: $f(t)=f\_0\; heta(-t)\; .$

We can express the expectation value of "x" by the probability distribution"W"("x",0) and the transition probability $P(x\text{'},t\; |\; x,0)$

: $langle\; x(t)\; angle\; =\; int\; dx\text{'}\; int\; dx\; ,\; x\text{'}\; P(x\text{'},t|x,0)\; W(x,0)\; .$

The probability distribution function "W"("x",0) is an equilibrium distribution and hencegiven by the

Boltzmann distribution for the Hamiltionian $H(x)=H\_0(x)\; +\; x\; f\_0$: $W(x,0)=\; frac\{exp(-eta\; H(x))\}\{int\; dx\; ,\; exp(-eta\; H(x))\}.$

For a weak field $eta\; f\_0\; ll\; 1$, we can expand the right-hand side

: $W(x,0)\; approx\; W\_0(x)\; (1+eta\; f\_0\; x),$

here $W\_0(x)$ is the equilibrium distribution in the absence of a field.Plugging this approximation in the formula for $langle\; x(t)\; angle$ yields

: $(*)\; langle\; x(t)\; angle\; =\; langle\; x\; angle\_0\; +\; eta\; f\_0\; A(t),$

where "A"("t") is the auto-correlation function of "x" in the absence of a field.

: $A(t)=langle\; x(t)\; x(0)\; angle.$

Note that in the absence of a field the system is invariant under time-shifts.We can rewrite $langle\; x(t)\; angle\; -\; langle\; x\; angle\_0$ using the susceptibilityof the system and hence find with the above equation (*)

: $f\_0\; int\_0^\{infty\}\; d\; au\; ,\; chi(\; au)\; heta(\; au-t)\; =\; eta\; f\_0\; A(t)$

Consequently,

: $(**)\; -chi(t)\; =\; eta\; frac\{d\}\{dt\}\; (\; A(t)\; heta(t)\; )\; .$

For

stationary process es, theWiener-Khinchin theorem states, thatthe power spectrum equals twice theFourier transform of the auto-correlationfunction: $S\_x(omega)\; =\; 2\; ilde\{A\}(omega).$

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 $frac\{d\}\{dt\}\; (\; A(t)\; heta(t)\; )$ into a symmetric and ananti-symmetric part

: $2\; frac\{d\}\{dt\}\; (\; A(t)\; heta(t)\; )\; =\; frac\{d\}\{dt\}\; A(t)\; +\; frac\{d\}\{dt\}(\; A(t)\; \{\; m\; sign\}(t)\; )\; .$

Now the fluctuation dissipation theorem follows.

**Derivation II**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:$mathrm\{Tr\}\; X\; =\; int\; dGamma\; X\; =\; sum\_\{mathrm\{quantum\; states\}\; n\}\; langle\; n|X|n\; angle$

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:$langle\; A\; angle\; =\; mathrm\{Tr\}\; left\{\; A(q,\; p)\; f(q,\; p)\; ight\}$

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:$langle\; A\; angle\; =\; mathrm\{Tr\}\; left\{\; A(q,\; p)\; f(q,\; p)\; ight\}\; =\; frac\{mathrm\{Tr\}\; left\{\; A\; e^\{-eta\; H\}\; ight\{mathrm\{Tr\}\; e^\{-eta\; H$

where β = 1/"k

_{B}T", "k"_{"B"}is theBoltzmann constant and "T" is thetemperature inKelvin .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 "k

_{B}T". 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; givenergodicity , 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

:$c(t)\; =\; langle\; delta\; A\; (t)\; delta\; A\; (0)\; angle\; =\; lim\_\{T\; ightarrowinfty\}\; frac\{1\}\{2T\}\; int\_\{-T\}^\{T\}\; dt^\{prime\}\; delta\; A\; (t^\{prime\}\; +\; t)\; delta\; A\; (t^\{prime\})$

which may be written as

:$c(t)\; =\; int\; dGamma\; delta\; A(t;\; q,\; p)\; delta\; A(0;\; q,\; p)\; f(q,\; p)$

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

:$ar\{A\}(t)\; =\; frac\{mathrm\{Tr\}\; left\{\; A(t;\; p,\; q)\; e^\{-eta\; left(H\; +\; h\; ight)\}\; ight\{mathrm\{Tr\}\; e^\{-eta\; left(H\; +\; h\; ight)$

Since "h" is much smaller than the thermal energy "k

_{B}T", we may expand the numerator:$mathrm\{Tr\}\; left\{\; A(t;\; p,\; q)\; e^\{-eta\; left(H\; +\; h\; ight)\}\; ight\}\; =mathrm\{Tr\}\; left\{\; A(t;\; p,\; q)\; e^\{-eta\; H\}\; left(1\; -\; eta\; h\; +\; cdots\; ight)\; ight\}\; approx\; mathrm\{Tr\}\; left\{\; A(t;\; p,\; q)\; e^\{-eta\; H\}\; ight\}\; -\; mathrm\{Tr\}\; left\{\; A(t;\; p,\; q)\; e^\{-eta\; H\}\; eta\; h\; ight\}$

We may likewise expand the denominator

:$frac\{1\}\{mathrm\{Tr\}\; e^\{-eta\; left(H\; +\; h\; ight)\; approx\; frac\{1\}\{mathrm\{Tr\}\; left\{\; e^\{-eta\; H\; \}\; ight\}\; -\; mathrm\{Tr\}\; left\{\; e^\{-eta\; H\; \}\; eta\; h\; ight\; =\; frac\{1\}\{mathrm\{Tr\}\; e^\{-eta\; H\; cdot\; frac\{1\}\{1\; -\; langle\; eta\; h\; angle\}\; approx\; frac\{1\; +\; langle\; eta\; h\; angle\}\{mathrm\{Tr\}\; e^\{-eta\; H$

where we have used

:$langle\; eta\; h\; angle\; =\; frac\{mathrm\{Tr\}\; left\{\; eta\; h\; e^\{-eta\; H\}\; ight\{mathrm\{Tr\}\; e^\{-eta\; H$

which is much less than one, by our assumption that "h" is much smaller than the thermal energy 1/β = "k

_{B}T".Combining the numerator and denominator, dropping quadratic and high-order terms in <β"h">, and using the indifference of equilibrium to time, we obtain

:$ar\{A\}(t)\; -\; langle\; A\; angle\; =\; -\; eta\; langle\; A\; h\; angle\; +\; eta\; langle\; A\; angle\; langle\; h\; angle$

Let the perturbation "h = −gA" be proportional to the variable "A" with a constant "−g". Then this formula becomes

:$ar\{A\}(t)\; -\; langle\; A\; angle\; =\; eta\; g\; left\{\; langle\; A(t)\; A(0)\; angle\; -\; langle\; A\; angle^\{2\}\; ight\}\; =\; eta\; g\; langle\; delta\; A(t)\; delta\; A(0)\; angle\; =\; eta\; g\; c(t).$

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 asFourier transform s andLaplace transform s.**ee also***

Non-equilibrium thermodynamics

*Green-Kubo relations

*Onsager reciprocal relations

*Equipartition theorem

*Boltzmann factor

*Dissipative system **References*** 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*]**Further resources*** [

*http://physics416.blogspot.com/2005/12/lecture-24-fluctuation-dissipation.html Audio recording*] of a lecture by Prof. E. W. Carlson ofPurdue University *

*

*

*

*

*

*

*

* [*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.*