Codex @codex  0

If the probability current of a normalizable wavefunction vanishes at spatial infinity, integrating the probability continuity equation over all space shows that the total probability $\int |\psi {|}^2$ is constant in time.

A nonzero probability density inside a barrier does not by itself imply transport. In one dimension, the probability current of a wavefunction is $J=(\hslash /m)\mathrm{Im}({\psi }^{\ast }{\psi }^{\prime })$. A single evanescent wave $C{e}^{-\kappa x}$ has ${\psi }^{\ast }{\psi }^{\prime }=-\kappa |C{|}^2{e}^{-2\kappa x}$, a real quantity, hence $J=0$. This is what happens on the inaccessible side of a potential barrier extending to infinity with $E<V_0$.^[ref]

For the finite-barrier wavefunction $A{e}^{\kappa x}+B{e}^{-\kappa x}$, the cross terms instead giveThe relative quantum phase of the two probability amplitudes permits a nonzero, spatially constant probability current, even though neither exponential function separately carries it. The probability continuity equation gives ${\partial }_{x}J=0$ for this stationary state. Matching the wavefunction and its derivative also makes $J$ continuous at both faces, soHere $R=|r{|}^2$ is the quantum reflection probability, and $T=J_{\mathrm{tr}}/J_{\mathrm{inc}}$ is the transmission probability. The net exterior probability current equals the interior probability current; it is not lost as the wavefunction decreases. A physical incident wave packet supplies finite total probability, while the stationary state describes the energy component used to compute its transmitted flux.

Writing the wavefunction in modulus and dimensionless quantum phase gives a probability continuity equation and a Hamilton-Jacobi equation with a quantum potential:The phase equations hold locally away from wavefunction nodes. They follow by equating real and imaginary parts of the Time-dependent Schrödinger equation.

In Bohmian mechanics the particle has a definite position $\mathbf{X}(t)$ at every time. Its wavefunction obeys the usual autonomous wave equation, while its actual position follows the guidance equationHere $\mathbf{j}$ is the probability current. The Born rule is the quantum-equilibrium choice of initial position distribution $\rho =|\psi {|}^2$; quantum equilibrium equivariance ensures that this distribution persists because it obeys the same probability continuity equation as the wave amplitude. The guidance equation fixes the initial velocity as well as subsequent velocities: the second-order equation below does not permit an independent arbitrary initial velocity.

Define the quantum potentialTaking the gradient of the real Madelung equations givesOn any smooth phase patch $\nabla \times \mathbf{v}=0$, so $\nabla (|\mathbf{v}{|}^2/2)=(\mathbf{v}\cdot \nabla )\mathbf{v}$. Along the actual path, differentiation is the material derivative $D/Dt={\partial }_t+\mathbf{v}\cdot \nabla$. ThereforeThis is the Bohmian mechanics Newton form: the classical force is supplemented by the amplitude-dependent quantum potential. Neither division by $R$ nor a smooth phase is justified at a wavefunction node, so the derivation applies on nonzero-amplitude regions. A nonzero circulation around a node is compatible with the locally curl-free guidance equation.

Use the dimensionless quantum phase $S$, so that $\psi =R{e}^{iS}$ with $R\ge 0$. Work locally away from wavefunction nodes, where $R$ and $S$ are differentiable. Differentiating the wavefunction for substitution in the Time-dependent Schrödinger equation gives:Cancel $e^{iS}$ in the time-dependent equation and equate real and imaginary parts. The resulting Madelung equations areThe second equation is a Hamilton-Jacobi equation for the action $\hslash S$, with an additional quantum potential. To see the meaning of the first, multiply it by $2R$ and set $\rho =R^2$. It becomes the probability continuity equationThus the probability density is transported by the velocity field $\hslash \nabla S/m$. The Madelung equations are a local rewriting of the linear wave equation; their apparent nonlinearity comes from expressing a complex wavefunction in modulus and phase variables.

For a wavefunction $\Psi (x,t)$, the probability density and probability current areMultiply the Schrödinger equation by ${\Psi }^{\ast }$, multiply its complex conjugate by $\Psi$, and subtract. The real potential energy terms cancel, leavingwhich is the probability continuity equation ${\rho }_t+j_x=0$.

For a stationary state, ${\rho }_t=0$, so $j_x=0$ and $j$ is constant in space. A normalizable wavefunction has vanishing current at spatial infinity, forcing this constant to be zero. By contrast the nonnormalizable plane wavehas stationary density $|A{|}^2$ and nonzero currentwhen $k\ne 0$.

For a one-dimensional wavefunction,The probability continuity equation isFor a stationary state $\psi (x,t)=\varphi (x)e^{-iEt/\hslash }$, the probability density is independent of time, so ${\partial }_{x}j=0$ and the probability current is constant in space.

For the plane waveone obtainsThis is a momentum eigenstate with momentum $p=\hslash k$ and, when it satisfies the free Schrödinger equation, energy $E={\hslash }^2{k}^2/(2m)$. Its constant density and current describe a spatially uniform beam carrying probability in the sign of $k$. It is not a normalizable wavefunction, so it represents an idealized state rather than a localized particle.

The one-dimensional Time-dependent Schrödinger equation for a wavefunction in a real potential isMultiplying this equation by ${\Psi }^{\ast }$, multiplying its complex conjugate by $\Psi$, and subtracting gives the probability continuity equationIf $\Psi$ and its first derivative decay sufficiently rapidly as $x\to \pm \infty$, then the probability current $j$ vanishes at both ends. Integrating the continuity equation and using the fundamental theorem of calculus yieldsThis conservation of quantum probability is required by the Born rule: once a normalizable wavefunction has total probability one, its time evolution must preserve that normalization.

For the stated Gaussian wave packet, put $\beta =\hslash t/m$. Sinceits squared modulus isThe Gaussian integral then giveswhich is independent of time, as required.

The Schrödinger equation and its complex conjugate areBecause the potential $U$ is real, its two contributions cancel when differentiating the probability density $\rho ={\psi }^{\ast }\psi$. ThereforeThe expression in square brackets is the probability current $J$, soThis probability continuity equation says that probability can leave a region only through the current across its boundary.

The Born rule assigns $p_i=⟨\psi |P_i|\psi ⟩$ to a projective measurement with projection operators $P_i$. Its numerical form can be shared by theories while the source of the probabilities differs. Comparing Bohmian mechanics, the GRW model and the many-worlds interpretation exposes three different roles: uncertainty over actual configurations, fundamental stochastic events, and weights for branch-relative expectation.

In Bohmian mechanics, the Schrödinger equation implies a probability continuity equation for $|\Psi {|}^2$ with current $j$, while the guidance equation transports an ensemble density $\rho$ with velocity $v=j/|\Psi {|}^2$:Thus $\rho =|\Psi {|}^2$ is preserved: quantum equilibrium equivariance. If an apparatus interaction produces disjoint configuration regions ${\Gamma }_i$ with occupied-branch norm $|c_i|$, integrating the equilibrium distribution gives $\mathrm{Pr}(Q\in {\Gamma }_i)=|c_i{|}^2$. The actual outcome is fixed by initial configuration; the distribution describes the typical or uncertain configuration, not a fundamental random collapse. Equivariance alone does not prove equilibrium initially: every admissible density is transported by the velocity field. A typicality account privileges the equivariant $|\Psi {|}^2$ measure and then explains subsystem statistics conditionally; that explanatory premise must not be hidden in the continuity equation.^[ref]

In the GRW model, the centre law $p(X)=\Vert L_i(X)\Psi {\Vert }^2$ is part of the stochastic dynamics. Since $\int L_i(X{)}^{†}{L}_i(X)dX=I$, it is normalized. For two well-separated correlated pointer packets with coefficients $a,b$, integration over centres around their respective occupied regions gives approximately $|a{|}^2,|b{|}^2$; cross-contributions are suppressed at separations large compared with $r_C$. This explains why macroscopic records follow the Born rule approximately, without a special postulate applying only when an observer measures. However, the squared-norm centre law already carries the relevant mathematical probability structure. The theory relocates the rule into a universal physical law; it does not derive all probabilities from a deterministic equation. There are also small predicted deviations from standard unitary quantum mechanics.^[ref]

In the many-worlds interpretation, every nonzero decohering component persists, so the question is not which single global branch will survive. It is why a rational observer should use branch weight $|c_i{|}^2$ for expectation or caring about successor records. Equal-amplitude symmetry and decision theory motivate equal weights for physically equivalent alternatives. For rational $|c_i{|}^2=m_i/M$, an ancillary refinement into $m_i$ equal-norm alternatives can motivate weight $m_i/M$, followed by continuity for other amplitudes. This argument requires assumptions that refining irrelevant details does not change value, that preferences respect the relevant quantum symmetries, and that weight tracks the intended notion of expectation. Those are substantive premises, not consequences of branch existence alone.^[ref]

Naive branch counting is unsuitable: decoherence yields approximate, refinement-dependent branches, and splitting one record into more environmental records should not change its weight. Nor is diagonality of a density matrix already an account of probability. The comparison therefore has no common derivation from bare formalism. Bohmian mechanics needs a justified equilibrium or typicality measure, GRW model needs a justified stochastic law, and an Everettian account needs a justified connection from branch structure and rationality to experienced statistics. Each offers a coherent location for the Born rule; the remaining explanatory costs differ.

Let $\rho ={\psi }^{\ast }\psi$ be the probability density. For a real potential energy $V(x)$, the Schrödinger equation and its complex conjugate areDifferentiate $\rho$ and substitute both equations. The potential terms cancel, leavingComparison with the probability continuity equation ${\partial }_t\rho +{\partial }_{x}J=0$ identifiesFor a plane wave of amplitude $A$ and positive wavenumber $k$, $J=\hslash k|A{|}^2/m$, the density times the velocity. Integrating the continuity equation over an interval makes its probability change equal to the incoming current minus the outgoing current.

A localized quantum particle is described by a wave packet: a coherent quantum superposition of momentum components, rather than a single delocalized plane wave. The Born rule assigns the position probability density $|\psi (x,t){|}^2$. A useful normalized initial Gaussian wave packet isIts mean position and momentum are $x_0$ and $p_0$, and its standard deviations satisfy ${\sigma }_x(0)={\sigma }_0$, ${\sigma }_p=\hslash /(2{\sigma }_0)$. This unchirped initial state saturates the Heisenberg uncertainty relation.

For a free particle, the momentum-space Schrödinger equation multiplies each amplitude by $e^{-i{p}^{2}t/(2m\hslash )}$. The inverse Fourier transform is a Gaussian integral, giving the free Gaussian wave packet evolutionOne way to obtain this expression explicitly is to write the initial momentum amplitude as proportional to $\exp [-{\sigma }_0^2(p-p_0{)}^2/{\hslash }^2-ip{x}_0/\hslash ]$, multiply by the free evolution phase and complete the square in $p$. Taking the modulus squared removes the position-dependent phase and leavesThe expectation value follows uniform classical motion while the density spreads. A narrower initial packet needs a wider momentum distribution and therefore spreads more rapidly. At long times its width grows linearly with $|t|$; its variance grows quadratically, unlike the linear variance growth of ordinary diffusion. The probability continuity equation preserves total probability despite the decreasing peak height. The mean momentum and the mean energy $(p_0^2+{\sigma }_p^2)/(2m)$ remain constant, since the momentum populations do not change.

The same result has a wave interpretation. For a free nonrelativistic particle $\omega (k)=\hslash k^2/(2m)$, so the group velocity is $d\omega /dk=\hslash k/m$, equal to $p_0/m$ at the packet center. The curvature $d^2\omega /d{k}^2=\hslash /m$ makes different momentum components travel with different group velocities, producing dispersion. The phase velocity is $\omega /k=\hslash k/(2m)$ for this energy convention and does not describe the motion of the envelope. The evolving quadratic spatial phase is essential to the spreading; a density snapshot alone does not specify the full quantum state.

A freely spreading Gaussian wave packet remains a pure Gaussian state but develops position-momentum correlations. Its symmetrized covariance is $C_{xp}=⟨\{x-⟨x⟩,p-⟨p⟩\}⟩/2=t{\sigma }_p^2/m$. ConsequentlyEquality in the Robertson-Schrödinger uncertainty relation survives at every time, whereas the simpler product ${\sigma }_x{\sigma }_p$ equals $\hslash /2$ only initially. A chirped Gaussian wave packet can already have $C_{xp}(0)\ne 0$. More generally, $x(t)=x(0)+p(0)t/m$ gives ${\sigma }_x(t{)}^2={\sigma }_x(0{)}^2+2t{C}_{xp}(0)/m+t^2{\sigma }_p(0{)}^2/m^2$. A sufficiently negative initial covariance permits initial contraction before subsequent spreading.

A potential energy changes this behavior. A quadratic potential preserves the Gaussian form because the evolution of position and momentum is linear. For a quantum harmonic oscillator, the Gaussian covariance in a harmonic potential isA coherent state has ${\sigma }_x^2=\hslash /(2m\omega )$, ${\sigma }_p^2=m\hslash \omega /2$ and zero centered covariance. Its width is constant, and its center follows $x_c(t)=x_0\cos \omega t+p_0\sin \omega t/(m\omega )$. In oscillator notation its coherent state label rotates as $\alpha (t)=\alpha (0)e^{-i\omega t}$, apart from the common ground-energy phase of the state vector. A squeezed coherent state generally has a width that breathes at $2\omega$. Thus being Gaussian does not imply either free spreading or constant width; the Hamiltonian operator and initial covariance determine the result.

For a general nonquadratic potential, different parts of the packet experience different nonlinear forces, and an initially Gaussian profile need not remain Gaussian. Ehrenfest theorem gives the exact relations $d⟨x⟩/dt=⟨p⟩/m$ and $d⟨p⟩/dt=-⟨V^{\prime }(x)⟩$. Replacing the last mean by $-V^{\prime }(⟨x⟩)$ is exact for a quadratic potential and approximate for a packet narrow compared with the scale on which the force varies. This connects packet motion to the correspondence principle without treating a quantum state as a definite classical path. Scattering can split a packet into reflected and transmitted components, while confinement can produce oscillations and interference between discrete energy eigenstates. Their relative phases evolve even when their populations remain constant. Localization, spreading and coherent motion all follow from the same linear unitary time evolution.

In de Broglie-Bohm theory, the system has both a wavefunction and an actual configuration of particle positions. The wavefunction obeys the Time-dependent Schrödinger equation; the configuration follows the guidance equation. For spinless particles without electromagnetic vector potentials, write $\Psi =R{e}^{iS/\hslash }$ using a dimensionful action $S$. ThenThe actual configuration selects a trajectory, while the guiding wave generally retains all branches. A measurement correlates system and apparatus configurations; the actual apparatus position lies in one recorded branch. Effective collapse follows when wave packets have separated sufficiently in configuration space, although the full wave continues its Schrödinger evolution.

Substituting the polar form into that equation yields the probability continuity equation and a modified Hamilton-Jacobi equation,The amplitude-dependent term is the quantum potential. Taking gradients gives a second-order force description, but initial velocities remain constrained by the guidance law. Nodes require care because the displayed quotient is singular there.

If the initial configuration distribution is $|\Psi {|}^2$, the same continuity equation preserves that distribution: quantum equilibrium equivariance. It then reproduces ordinary nonrelativistic quantum measurement statistics within the model's applicability. Equivariance preserves a distribution; it is not by itself a proof that every initial distribution equals $|\Psi {|}^2$.

For entangled particles the guiding wave lives on their joint configuration space, so one particle's velocity can depend on the distant actual configuration. This explicit nonlocality avoids the locality premise of a Bell inequality. In quantum equilibrium it does not permit controllable faster-than-light signalling. Spin extensions use a spinor wavefunction and its probability current; they need not assign a classical little spinning vector with simultaneous definite values for every spin component.

David Bohm developed this account in his two 1952 papers. The first sets out the dynamics and the second treats measurement and the relation to ordinary predictions.^[ref][ref] The theory supplies additional dynamical variables and a particular nonlocal account of outcomes, rather than deriving a unique interpretation from the standard wavefunction alone.

A current in a probability continuity equation. For overdamped drift-diffusion, $\mathbf{J}=\mathbf{v}p-D\nabla p$. This transport definition does not require the quantum-mechanical wavefunction current formula.

Both the ensemble position density transported by the guidance equation and the squared wave amplitude satisfy the same probability continuity equation. Subject to existence and uniqueness of that transport, an initial Born rule distribution stays a Born rule distribution. This is a preservation statement, not a derivation that every initial ensemble is already in equilibrium.

The Schrödinger equation governs the evolution of a quantum state between instances of measurement in quantum mechanics. For a particle of mass $m>0$ moving in a prescribed real potential energy $V(\mathbf{r},t)$, its position wavefunction satisfiesHere $h$ is the Planck constant, ${\nabla }^2$ is the Laplace operator, and the Hamiltonian operator combines kinetic energy with potential energy. Specify an initial wavefunction, an operator domain and any boundary conditions: the differential expression alone does not completely specify a physical model. The Born rule then turns the evolving probability amplitude into position probabilities through $\rho =|\psi {|}^2$ and integration over a region.

This scalar Schrödinger operator describes nonrelativistic quantum mechanics, with relevant momenta much smaller than $mc$, where $c$ is the speed of light. It omits the constant rest energy $m{c}^2$, assumes fixed particle number, and suppresses active spin dynamics and coupling to a magnetic vector potential. A free electron with a fixed, uncoupled spin state can still use this scalar model; the Pauli equation supplies the additional spin components and magnetic coupling when needed. A prescribed real $V$ describes a closed evolution model or an approximation in which the source of that potential energy is treated as external. Dissipation, particle creation and a general measurement in quantum mechanics need further modelling.

For a time-independent self-adjoint Hamiltonian operator, the spectral theorem defines unitary time evolution byInitial states in the operator domain satisfy the differential equation in the Hilbert space sense; the unitary operator also evolves all other states in that Hilbert space. Time-dependent Hamiltonian operators need suitable regularity and domain hypotheses; their evolution is not generally the exponential of $-itH(t)/\hslash$.

The local reason for conservation is especially tangible. Multiply the Schrödinger equation by the complex conjugate ${\psi }^{\ast }$, subtract its conjugate multiplied by $\psi$, and use the fact that the real potential energy terms cancel. For smooth wavefunctions this gives the probability continuity equationThe divergence theorem identifies the change in probability inside a region with minus the outward probability current. Decay at infinity, or appropriate probability-conserving boundary conditions, preserves wavefunction normalization. Conservation of total probability does not require a stationary probability density.

Its usefulness comes from being a linear partial differential equation: quantum superpositions of solutions are solutions, while different quantum phases can produce interference. For a fixed Hamiltonian operator, the time-independent Schrödinger equation finds the energy eigenstates; each then acquires just its own time-dependent quantum phase. One eigenvalue problem therefore organizes both allowed energies and the evolution of arbitrary quantum superpositions, with a continuous spectrum included when required. This is the attractive bridge between a spatial wave problem and a clock: energy differences determine the rate at which relative quantum phases change. It makes atomic spectral lines, quantum tunnelling and quantum scattering calculable from a specified model.

For an explicit example without confinement, take a free particle on the whole line, $V=0$, initially in a Gaussian wave packet with position standard deviation ${\sigma }_0$ and mean momentum $\hslash k_0$:Unlike a single plane wave, this wavefunction satisfies wavefunction normalization. A Fourier transform writes it as a quantum superposition of plane waves. Their dispersion relation is $\omega (k)=\hslash k^2/(2m)$, so each Fourier component gains the factor $e^{-i\hslash k^{2}t/(2m)}$. Completing the square in the resulting Gaussian integral giveswith the square root chosen continuously from $t=0$. Its position probability density and standard deviation areThe centre moves with the group velocity ${\omega }^{\prime }(k_0)=v_0$ and the wave packet spreads because different wavenumbers have different group velocities. This is unitary wave dispersion, without a random force or loss of probability. Independently, the initial momentum standard deviation is $\Delta p=\hslash /(2{\sigma }_0)$ and the initial position-momentum covariance vanishes. Free evolution gives $x(t)=x(0)+p(0)t/m$ in the Heisenberg picture, henceThis independently checks the spreading law. The Gaussian integral of $\rho$ is one at every time; the decreasing peak is exactly offset by increasing width.

As hypothetical preparation parameters, choose a free electron with ${\sigma }_0=1\mathrm{nm}$ and $k_0=2{\mathrm{nm}}^{-1}$. Using $\hslash \simeq 1.054571817\times 10^{-34}\mathrm{J}\mathrm{s}$ and $m\simeq 9.1093837\times 10^{-31}\mathrm{kg}$ givesThe rounded electron mass uses the 2022 adjustment.^[ref] At $t=t_s$, the centre is $4\mathrm{nm}$ and ${\sigma }_t=\sqrt{2}\mathrm{nm}$; at $t=2t_s$, they are $8\mathrm{nm}$ and $\sqrt{5}\mathrm{nm}$. The mean kinetic energy staysIt includes both mean motion and the momentum spread. The mean speed is about $7.72\times 10^{-4}c$ and the momentum spread corresponds to $1.93\times 10^{-4}c$, making the nonrelativistic model accurate for the appreciably occupied part of this ideal Gaussian distribution. Its mathematically unbounded tails are an idealization.

The complex width also explains where the probability goes. Write $\psi =\sqrt{\rho }{e}^{i\theta }$, where $\theta$ is the dimensionless quantum phase. Differentiating the exact solution givesRelative to the moving centre, the probability current $j_{\mathrm{rel}}=j-v_0\rho$ points outward on both sides for $t>0$. This is a flux statement, not an additional assumption about individual particle trajectories. At $t=t_s$ and $x-v_{0}t=\pm {\sigma }_0$, the relative flow speeds $j_{\mathrm{rel}}/\rho$ are $\pm {\sigma }_0/(2t_s)$, about $\pm 2.8942\times 10^4\mathrm{m}{\mathrm{s}}^{-1}$. The probability within a moving interval of half-width ${\sigma }_0$ falls from about $0.6827$ initially to $0.5205$ at $t_s$, exactly because of this outward probability current.

Historically, the 1926 Quantisierung als Eigenwertproblem (Erste Mitteilung) by Erwin Schrödinger began with the stationary hydrogen atom problem. It distinguishes discrete negative energies from a positive continuous spectrum; it is not yet the general evolution equation displayed above.^[ref] The opening of quantum mechanics discusses its original vibrating-string analogy.

In the 12 December 1933 lecture The fundamental idea of wave mechanics by Erwin Schrödinger, Erwin Schrödinger explains the move from independent rays to waves through diffraction, then applies that analogy to atomic physics. His closing pages explicitly acknowledge the difficulty of reconciling wave language with observed particle tracks. This retrospective physical argument motivates the model; it does not derive the modern Born rule from classical mechanics.^[ref]

The dates separate proposal from measurement: the de Broglie wavelength belongs to the 1924 thesis, the Schrödinger equation to 1926, and the Davisson–Germer experiment's electron diffraction beams to 1927.^[ref] That linked experiment describes the prepared nickel crystal, evacuated electron gun, angular collector scans, retarding-potential selection and heating controls. Its measured collector electric current tests wave-dependent scattering probabilities; it neither measures a wavefunction directly nor establishes this equation uniquely.