Einstein relation (kinetic theory) explained

In physics (specifically, the kinetic theory of gases), the Einstein relation is a previously unexpected connection revealed independently by William Sutherland in 1904,^{[1] [2] [3]} Albert Einstein in 1905,^[4] and by Marian Smoluchowski in 1906^[5] in their works on Brownian motion. The more general form of the equation in the classical case is^[6]

$$D\; =\; \backslash mu\; \backslash ,\; k\_\backslash text\; T,$$where

• is the diffusion coefficient;
• is the "mobility", or the ratio of the particle's terminal drift velocity to an applied force, ;
• is the Boltzmann constant;
• is the absolute temperature.

This equation is an early example of a fluctuation-dissipation relation.^[7] Note that the equation above describes the classical case and should be modified when quantum effects are relevant.

Two frequently used important special forms of the relation are:

• Einstein–Smoluchowski equation, for diffusion of charged particles:^[8] $$D\; =\; \backslash frac$$
• Stokes–Einstein–Sutherland equation, for diffusion of spherical particles through a liquid with low Reynolds number: $$D\; =\; \backslash frac$$

Here

• is the electrical charge of a particle;
• is the electrical mobility of the charged particle;
• is the dynamic viscosity;
• is the radius of the spherical particle.

Special cases

Electrical mobility equation (classical case)

For a particle with electrical charge, its electrical mobility is related to its generalized mobility by the equation . The parameter is the ratio of the particle's terminal drift velocity to an applied electric field. Hence, the equation in the case of a charged particle is given as$$D\; =\; \backslash frac,$$

where

D is the diffusion coefficient (m ^{2} s ^{-1} ).\mu _{q} is the electrical mobility (m ^{2} V ^{-1} s ^{-1} ).q is the electric charge of particle (C, coulombs)T is the electron temperature or ion temperature in plasma (K).^[9]

If the temperature is given in volts, which is more common for plasma:$$D\; =\; \backslash frac,$$where

Z is the charge number of particle (unitless)T is electron temperature or ion temperature in plasma (V).

Electrical mobility equation (quantum case)

For the case of Fermi gas or a Fermi liquid, relevant for the electron mobility in normal metals like in the free electron model, Einstein relation should be modified:$$D\; =\; \backslash frac,$$where

E _{F} is Fermi energy.

Stokes–Einstein–Sutherland equation

In the limit of low Reynolds number, the mobility μ is the inverse of the drag coefficient

\zeta . A damping constant \gamma = \zeta / m is frequently used for the inverse momentum relaxation time (time needed for the inertia momentum to become negligible compared to the random momenta) of the diffusive object. For spherical particles of radius r, Stokes' law gives$$\backslash zeta\; =\; 6\; \backslash pi\; \backslash ,\; \backslash eta\; \backslash ,\; r,$$where η is the viscosity of the medium. Thus the Einstein–Smoluchowski relation results into the Stokes–Einstein–Sutherland relation$$D\; =\; \backslash frac.$$This has been applied for many years to estimating the self-diffusion coefficient in liquids, and a version consistent with isomorph theory has been confirmed by computer simulations of the Lennard-Jones system.^[10]

In the case of rotational diffusion, the friction is

\zeta _{r} = 8 \pi η r ^{3} , and the rotational diffusion constant D _{r} is$$D\_\backslash text\; =\; \backslash frac.$$This is sometimes referred to as the Stokes–Einstein–Debye relation.

Semiconductor

In a semiconductor with an arbitrary density of states, i.e. a relation of the form

p = p(\varphi) between the density of holes or electrons p and the corresponding quasi Fermi level (or electrochemical potential) \varphi , the Einstein relation is^{[11] [12]} $$D\; =\; \backslash frac,$$where \mu _{q} is the electrical mobility (see for a proof of this relation). An example assuming a parabolic dispersion relation for the density of states and the Maxwell–Boltzmann statistics, which is often used to describe inorganic semiconductor materials, one can compute (see density of states):$$p(\backslash varphi)\; =\; N\_0\; e^,$$where N _{0} is the total density of available energy states, which gives the simplified relation:$$D\; =\; \backslash mu\_q\; \backslash frac.$$

Nernst–Einstein equation

By replacing the diffusivities in the expressions of electric ionic mobilities of the cations and anions from the expressions of the equivalent conductivity of an electrolyte the Nernst–Einstein equation is derived:$$\backslash Lambda\_e\; =\; \backslash frac(D\_+\; +\; D\_-).$$were R is the gas constant.

Proof of the general case

The proof of the Einstein relation can be found in many references, for example see the work of Ryogo Kubo.^[13]

U generates a conservative force F( x ) =- \nabla U( x ) (for example, an electric force) on a particle located at a given position x . We assume that the particle would respond by moving with velocity v( x ) = \mu( x ) F( x ) (see Drag (physics)). Now assume that there are a large number of such particles, with local concentration \rho( x ) as a function of the position. After some time, equilibrium will be established: particles will pile up around the areas with lowest potential energy U , but still will be spread out to some extent because of diffusion. At equilibrium, there is no net flow of particles: the tendency of particles to get pulled towards lower U , called the drift current, perfectly balances the tendency of particles to spread out due to diffusion, called the diffusion current (see drift-diffusion equation).

The net flux of particles due to the drift current is$$\backslash mathbf\_\backslash mathrm(\backslash mathbf)\; =\; \backslash mu(\backslash mathbf)\; F(\backslash mathbf)\; \backslash rho(\backslash mathbf)\; =\; -\backslash rho(\backslash mathbf)\; \backslash mu(\backslash mathbf)\; \backslash nabla\; U(\backslash mathbf),$$i.e., the number of particles flowing past a given position equals the particle concentration times the average velocity.

The flow of particles due to the diffusion current is, by Fick's law,$$\backslash mathbf\_\backslash mathrm(\backslash mathbf)=-D(\backslash mathbf)\; \backslash nabla\; \backslash rho(\backslash mathbf),$$where the minus sign means that particles flow from higher to lower concentration.

Now consider the equilibrium condition. First, there is no net flow, i.e.

J _{drift} + J _{diffusion} = 0 . Second, for non-interacting point particles, the equilibrium density \rho is solely a function of the local potential energy U , i.e. if two locations have the same U then they will also have the same \rho (e.g. see Maxwell-Boltzmann statistics as discussed below.) That means, applying the chain rule,$$\backslash nabla\backslash rho\; =\; \backslash frac\; \backslash nabla\; U.$$

Therefore, at equilibrium:$$0\; =\; \backslash mathbf\_\backslash mathrm\; +\; \backslash mathbf\_\backslash mathrm\; =\; -\backslash mu\; \backslash rho\; \backslash nabla\; U\; -\; D\; \backslash nabla\; \backslash rho\; =\; \backslash left(-\backslash mu\; \backslash rho\; -\; D\; \backslash frac\backslash right)\backslash nabla\; U.$$

As this expression holds at every position

x , it implies the general form of the Einstein relation:$$D\; =\; -\backslash mu\; \backslash frac.$$

The relation between

\rho and U for classical particles can be modeled through Maxwell-Boltzmann statistics$$\backslash rho(\backslash mathbf)\; =\; A\; e^,$$where A is a constant related to the total number of particles. Therefore$$\backslash frac\; =\; -\backslash frac\backslash rho.$$

Under this assumption, plugging this equation into the general Einstein relation gives:$$D\; =\; -\backslash mu\; \backslash frac\; =\; \backslash mu\; k\_\backslash text\; T,$$which corresponds to the classical Einstein relation.

See also

• Smoluchowski factor
• Conductivity (electrolytic)
• Stokes radius
• Ion transport number

External links

• Einstein relation calculators
• ion diffusivity

Notes and References

1. http://www.ph.unimelb.edu.au/~dnj/wyop/wyop2005-sutherland-essay.html World Year of Physics – William Sutherland at the University of Melbourne
2. 10.1080/14786440509463331 . 9 . 54 . LXXV. A dynamical theory of diffusion for non-electrolytes and the molecular mass of albumin . 1905 . Philosophical Magazine . Series 6 . 781–785 . Sutherland William.
3. P. Hänggi, "Stokes–Einstein–Sutherland equation".
4. Einstein, A.. Über die von der molekularkinetischen Theorie der Wärme geforderte Bewegung von in ruhenden Flüssigkeiten suspendierten Teilchen. Annalen der Physik. 322. 8. 549–560. 1905. de. 10.1002/andp.19053220806 . 1905AnP...322..549E . free.
5. von Smoluchowski, M.. Zur kinetischen Theorie der Brownschen Molekularbewegung und der Suspensionen. Annalen der Physik. 326. 14. 756–780. 1906. de. 10.1002/andp.19063261405 . 1906AnP...326..756V .
6. Book: Molecular Driving Forces: Statistical Thermodynamics in Chemistry and Biology . Dill . Ken A. . Bromberg . Sarina . 2003 . Garland Science . 9780815320517 . 327 . en.
7. Umberto Marini Bettolo Marconi, Andrea Puglisi, Lamberto Rondoni, Angelo Vulpiani, "Fluctuation-Dissipation: Response Theory in Statistical Physics".
8. Van Zeghbroeck, "Principles of Semiconductor Devices", Chapter 2.7 .
9. Book: Raizer, Yuri . Gas Discharge Physics . Springer . 2001 . 978-3540194620 . 20–28.
10. Costigliola. Lorenzo. Heyes. David M.. Schrøder. Thomas B.. Dyre. Jeppe C.. 2019-01-14. Revisiting the Stokes-Einstein relation without a hydrodynamic diameter. The Journal of Chemical Physics . en. 150. 2. 021101. 10.1063/1.5080662. 30646717. 2019JChPh.150b1101C . 0021-9606. free.
11. Book: Ashcroft, N. W. . Mermin, N. D.. Solid State Physics. 826. 1988. Holt, Rineheart and Winston. New York (USA).
12. Book: Bonnaud, Olivier. Composants à semiconducteurs. 78. 2006. Ellipses. Paris (France). fr.
13. Kubo, R.. The fluctuation-dissipation theorem. Rep. Prog. Phys.. 29. 1. 255–284. 1966. 1966RPPh...29..255K . 10.1088/0034-4885/29/1/306 . 250892844. 0710.4394.