Madelung equations explained

In theoretical physics, the Madelung equations, or the equations of quantum hydrodynamics, are Erwin Madelung's equivalent alternative formulation of the Schrödinger equation for a spinless non relativistic particle, written in terms of hydrodynamical variables, similar to the Navier–Stokes equations of fluid dynamics. The derivation of the Madelung equations is similar to the de Broglie–Bohm formulation, which represents the Schrödinger equation as a quantum Hamilton–Jacobi equation.

History

In the fall of 1926, Erwin Madelung reformulated Schrödinger's quantum equation in a more classical and visualizable form resembling hydrodynamics. His paper was one of numerous early attempts at different approaches to quantum mechanics, including those of Louis de Broglie and Earle Hesse Kennard. The most influential of these theories was ultimately de Broglie's through the 1952 work of David Bohm now called Bohmian mechanics

Equations

The Madelung equations are quantum Euler equations:\begin& \partial_t \rho_m + \nabla\cdot(\rho_m \mathbf v) = 0, \\[4pt]& \frac = \partial_t\mathbf v + \mathbf v \cdot \nabla\mathbf v = -\frac \mathbf(Q + V),\endwhere

v

is the flow velocity,

\rhom=m\rho=m|\psi|2

is the mass density,

Q=-

\hbar2
2m
\nabla2\sqrt{\rho
} = -\frac \frac is the Bohm quantum potential,

The Madelung equations answer the question whether

v(x,t)

obeys the continuity equations of hydrodynamics and, subsequently, what plays the role of the stress tensor.

The circulation of the flow velocity field along any closed path obeys the auxiliary quantization condition \Gamma \doteq \oint = 2\pi n\hbar for all integers .

Derivation

The Madelung equations are derived by first writing the wavefunction in polar form\psi(\mathbf, t) = R(\mathbf, t) e^, with

R\geq0

and

S

both real and \rho(\mathbf,t)=\psi(\mathbf,t)^*\psi(\mathbf,t)=R^2(\mathbf,t),the associated probability density. Substituting this form into the probability current gives:\mathbf = \frac(\psi^* \nabla \psi - \psi \nabla \psi^*) = \frac\rho(\mathbf,t)\nabla S(\mathbf,t) = \rho(\mathbf,t)\mathbf(\mathbf,t),where the flow velocity is expressed as\mathbf(\mathbf,t)=\frac\nabla S(\mathbf,t).However, the interpretation of

v

as a "velocity" should not be taken too literal, because a simultaneous exact measurement of position and velocity would necessarily violatethe uncertainty principle.

Next, substituting the polar form into the Schrödinger equationi\hbar\frac \psi(\mathbf,t) = \left[\frac{-\hbar^2}{2m} \nabla^2 + V(\mathbf{x}) \right] \psi(\mathbf, t),and performing the appropriate differentiations, dividing the equation by

e

and separating the real and imaginary parts, one obtains a system of two coupled partial differential equations:\begin&\partial_R(\mathbf,t) + \frac\nabla R(\mathbf,t)\cdot\nabla S(\mathbf,t) + \frac R(\mathbf,t)\Delta S(\mathbf,t) = 0,\\&\partial_S(\mathbf,t) + \frac\left[\nabla S(\mathbf{x},t)\right]^2 + V(\mathbf) = \frac\frac.\endThe first equation corresponds to the imaginary part of Schrödinger equation and can be interpreted as the continuity equation. The second equation corresponds to the real part and is also referred to as the quantum Hamilton-Jacobi equation.Multiplying the first equation by

2R

and calculating the gradient of the second equation results in the Madelung equations:\begin&\partial_\rho(\mathbf,t) + \nabla\cdot\left[\rho(\mathbf{x},t)v(\mathbf{x},t) \right]= 0,\\&\frac\mathbf(\mathbf,t)=\partial_v(\mathbf,t) + \left[v(\mathbf{x},t)\cdot \nabla\right]v(\mathbf,t) = -\frac\nabla \left[V(\mathbf{x}) - \frac{\hbar^2}{2m}\frac{\Delta \sqrt{\rho(\mathbf{x},t)}}{\sqrt{\rho(\mathbf{x},t)}}\right] =-\frac\nabla \left[V(\mathbf{x}) + Q(\mathbf{x},t)\right].\endwith quantum potentialQ(\mathbf,t) = - \frac\frac.

Alternatively, the quantum Hamilton-Jacobi equation can be written in a form similar to the Cauchy momentum equation:\frac\mathbf= \mathbf - \frac \nabla \cdot \mathbf_Q,with an external force defined as\mathbf(\mathbf)=-\frac\nabla V(\mathbf),and a quantum pressure tensor \mathbf_Q = - (\hbar/2m)^2 \rho_m \nabla \otimes \nabla \ln \rho_m.

The integral energy stored in the quantum pressure tensor is proportional to the Fisher information, which accounts for the quality of measurements. Thus, according to the Cramér–Rao bound, the Heisenberg uncertainty principle is equivalent to a standard inequality for the efficiency of measurements.

Quantum energies

The thermodynamic definition of the quantum chemical potential\mu = Q + V = \frac \widehat H \sqrtfollows from the hydrostatic force balance above:\nabla \mu = \frac \nabla \cdot \mathbf p_Q + \nabla V.According to thermodynamics, at equilibrium the chemical potential is constant everywhere, which corresponds straightforwardly to the stationary Schrödinger equation. Therefore, the eigenvalues of the Schrödinger equation are free energies, which differ from the internal energies of the system. The particle internal energy is calculated as\varepsilon = \mu - \operatorname(\mathbf p_Q) \frac = -\frac (\nabla \ln \rho_m)^2 + Uand is related to the local Carl Friedrich von Weizsäcker correction.

See also

References