For faster navigation, this Iframe is preloading the Wikiwand page for Green–Kubo relations.

Green–Kubo relations

This article needs attention from an expert in Physics. Please add a reason or a talk parameter to this template to explain the issue with the article. WikiProject Physics may be able to help recruit an expert. (April 2023)

The Green–Kubo relations (Melville S. Green 1954, Ryogo Kubo 1957) give the exact mathematical expression for a transport coefficient in terms of the integral of the equilibrium time correlation function of the time derivative of a corresponding microscopic variable (sometimes termed a "gross variable", as in [1]):

One intuitive way to understand this relation is that relaxations resulting from random fluctuations in equilibirum are indistinguishable from those due to an external perturbation in linear response.[2]

Green-Kubo relations are important because they relate a macroscopic transport coefficient to the correlation function of a microscopic variable. In addition, they allow one to measure the transport coefficient without perturbing the system out of equilibrium, which has found much use in molecular dynamics simulations.[3]

Thermal and mechanical transport processes

Thermodynamic systems may be prevented from relaxing to equilibrium because of the application of a field (e.g. electric or magnetic field), or because the boundaries of the system are in relative motion (shear) or maintained at different temperatures, etc. This generates two classes of nonequilibrium system: mechanical nonequilibrium systems and thermal nonequilibrium systems.

The standard example of an electrical transport process is Ohm's law, which states that, at least for sufficiently small applied voltages, the current I is linearly proportional to the applied voltage V,

As the applied voltage increases one expects to see deviations from linear behavior. The coefficient of proportionality is the electrical conductance which is the reciprocal of the electrical resistance.

The standard example of a mechanical transport process is Newton's law of viscosity, which states that the shear stress is linearly proportional to the strain rate. The strain rate is the rate of change streaming velocity in the x-direction, with respect to the y-coordinate, . Newton's law of viscosity states

As the strain rate increases we expect to see deviations from linear behavior

Another well known thermal transport process is Fourier's law of heat conduction, stating that the heat flux between two bodies maintained at different temperatures is proportional to the temperature gradient (the temperature difference divided by the spatial separation).

Linear constitutive relation

Regardless of whether transport processes are stimulated thermally or mechanically, in the small field limit it is expected that a flux will be linearly proportional to an applied field. In the linear case the flux and the force are said to be conjugate to each other. The relation between a thermodynamic force F and its conjugate thermodynamic flux J is called a linear constitutive relation,

L(0) is called a linear transport coefficient. In the case of multiple forces and fluxes acting simultaneously, the fluxes and forces will be related by a linear transport coefficient matrix. Except in special cases, this matrix is symmetric as expressed in the Onsager reciprocal relations.

In the 1950s Green and Kubo proved an exact expression for linear transport coefficients which is valid for systems of arbitrary temperature T, and density. They proved that linear transport coefficients are exactly related to the time dependence of equilibrium fluctuations in the conjugate flux,

where (with k the Boltzmann constant), and V is the system volume. The integral is over the equilibrium flux autocovariance function. At zero time the autocovariance is positive since it is the mean square value of the flux at equilibrium. Note that at equilibrium the mean value of the flux is zero by definition. At long times the flux at time t, J(t), is uncorrelated with its value a long time earlier J(0) and the autocorrelation function decays to zero. This remarkable relation is frequently used in molecular dynamics computer simulation to compute linear transport coefficients; see Evans and Morriss, "Statistical Mechanics of Nonequilibrium Liquids", Academic Press 1990.

Nonlinear response and transient time correlation functions

In 1985 Denis Evans and Morriss derived two exact fluctuation expressions for nonlinear transport coefficients—see Evans and Morriss in Mol. Phys, 54, 629(1985). Evans later argued that these are consequences of the extremization of free energy in Response theory as a free energy minimum.[4]

Evans and Morriss proved that in a thermostatted system that is at equilibrium at t = 0, the nonlinear transport coefficient can be calculated from the so-called transient time correlation function expression:

where the equilibrium () flux autocorrelation function is replaced by a thermostatted field dependent transient autocorrelation function. At time zero but at later times since the field is applied .

Another exact fluctuation expression derived by Evans and Morriss is the so-called Kawasaki expression for the nonlinear response:

The ensemble average of the right hand side of the Kawasaki expression is to be evaluated under the application of both the thermostat and the external field. At first sight the transient time correlation function (TTCF) and Kawasaki expression might appear to be of limited use—because of their innate complexity. However, the TTCF is quite useful in computer simulations for calculating transport coefficients. Both expressions can be used to derive new and useful fluctuation expressions quantities like specific heats, in nonequilibrium steady states. Thus they can be used as a kind of partition function for nonequilibrium steady states.

Derivation from the fluctuation theorem and the central limit theorem [clarification needed]

For a thermostatted steady state, time integrals of the dissipation function are related to the dissipative flux, J, by the equation

We note in passing that the long time average of the dissipation function is a product of the thermodynamic force and the average conjugate thermodynamic flux. It is therefore equal to the spontaneous entropy production in the system. The spontaneous entropy production plays a key role in linear irreversible thermodynamics – see de Groot and Mazur "Non-equilibrium thermodynamics" Dover.

The fluctuation theorem (FT) is valid for arbitrary averaging times, t. Let's apply the FT in the long time limit while simultaneously reducing the field so that the product is held constant,

Because of the particular way we take the double limit, the negative of the mean value of the flux remains a fixed number of standard deviations away from the mean as the averaging time increases (narrowing the distribution) and the field decreases. This means that as the averaging time gets longer the distribution near the mean flux and its negative, is accurately described by the central limit theorem. This means that the distribution is Gaussian near the mean and its negative so that

Combining these two relations yields (after some tedious algebra!) the exact Green–Kubo relation for the linear zero field transport coefficient, namely,

Here are the details of the proof of Green–Kubo relations from the FT.[5] A proof using only elementary quantum mechanics was given by Robert Zwanzig.[6]


This shows the fundamental importance of the fluctuation theorem (FT) in nonequilibrium statistical mechanics. The FT gives a generalisation of the second law of thermodynamics. It is then easy to prove the second law inequality and the Kawasaki identity. When combined with the central limit theorem, the FT also implies the Green–Kubo relations for linear transport coefficients close to equilibrium. The FT is, however, more general than the Green–Kubo Relations because, unlike them, the FT applies to fluctuations far from equilibrium. In spite of this fact, no one has yet been able to derive the equations for nonlinear response theory from the FT.

The FT does not imply or require that the distribution of time-averaged dissipation is Gaussian. There are many examples known when the distribution is non-Gaussian and yet the FT still correctly describes the probability ratios.

See also


This article includes a list of general references, but it lacks sufficient corresponding inline citations. Please help to improve this article by introducing more precise citations. (December 2010) (Learn how and when to remove this message)
  1. ^ Green, Melville S. (1954). "Markoff Random Processes and the Statistical Mechanics of Time-Dependent Phenomena. II. Irreversible Processes in Fluids". The Journal of Chemical Physics. 22 (3): 398–413. Bibcode:1954JChPh..22..398G. doi:10.1063/1.1740082. ISSN 0021-9606.
  2. ^ Evans DJ, Morriss G (2008). Statistical Mechanics of Nonequilibrium Liquids (Second ed.). Cambridge University Press. ISBN 978-0-521-85791-8.
  3. ^ Nevins, D.; Spera, F. J. (December 2007). "Accurate computation of shear viscosity from equilibrium molecular dynamics simulations". Molecular Simulation. 33 (15): 1261–1266. doi:10.1080/08927020701675622. ISSN 0892-7022. Retrieved 8 November 2023.
  4. ^ Evans, Denis J. (1985-11-01). "Response theory as a free-energy extremum". Physical Review A. 32 (5): 2923–2925. Bibcode:1985PhRvA..32.2923E. doi:10.1103/physreva.32.2923. ISSN 0556-2791. PMID 9896433.
  5. ^ Evans, Denis J.; Searles, Debra J.; Rondoni, Lamberto (2005). "Application of the Gallavotti-Cohen fluctuation relation to thermostated steady states near equilibrium". Physical Review E. 71 (5): 056120. arXiv:cond-mat/0312353. Bibcode:2005PhRvE..71e6120E. doi:10.1103/PhysRevE.71.056120. PMID 16089615. S2CID 4617097.
  6. ^ Zwanzig, R. (1965). "Time-Correlation Functions and Transport Coefficients in Statistical Mechanics". Annual Review of Physical Chemistry. 16: 67–102. Bibcode:1965ARPC...16...67Z. doi:10.1146/annurev.pc.16.100165.000435.
{{bottomLinkPreText}} {{bottomLinkText}}
Green–Kubo relations
Listen to this article

This browser is not supported by Wikiwand :(
Wikiwand requires a browser with modern capabilities in order to provide you with the best reading experience.
Please download and use one of the following browsers:

This article was just edited, click to reload
This article has been deleted on Wikipedia (Why?)

Back to homepage

Please click Add in the dialog above
Please click Allow in the top-left corner,
then click Install Now in the dialog
Please click Open in the download dialog,
then click Install
Please click the "Downloads" icon in the Safari toolbar, open the first download in the list,
then click Install

Install Wikiwand

Install on Chrome Install on Firefox
Don't forget to rate us

Tell your friends about Wikiwand!

Gmail Facebook Twitter Link

Enjoying Wikiwand?

Tell your friends and spread the love:
Share on Gmail Share on Facebook Share on Twitter Share on Buffer

Our magic isn't perfect

You can help our automatic cover photo selection by reporting an unsuitable photo.

This photo is visually disturbing This photo is not a good choice

Thank you for helping!

Your input will affect cover photo selection, along with input from other users.


Get ready for Wikiwand 2.0 🎉! the new version arrives on September 1st! Don't want to wait?