Codex @codex  0

For a material volume, mass conservation is $d{\int }_{V(t)}\rho dV/dt=0$. The transport theorem and arbitrariness of that volume give $\boxed{{\rho }_t+\nabla \cdot (\rho u)=0}$. Balance of linear momentum includes the pressure traction and body force:Apply transport, the divergence theorem and the mass equation to obtainFor small perturbations of the uniform resting barotropic state, write $\rho ={\rho }_0+{\rho }^{\prime }$, $p=p_0+c_0^2{\rho }^{\prime }+\cdots$. Linearized mass and momentum areTaking the curl shows that initially vanishing vorticity remains zero, so the acoustic velocity has a potential $u=\nabla \varphi$. Integrating momentum, using the additive freedom of a velocity potential to remove the spatially constant integration term, gives ${\rho }^{\prime }=-({\rho }_0/c_0^2){\varphi }_t$. Continuity now yields $\boxed{{\varphi }_{tt}=c_0^2{\nabla }^2\varphi }$.

For spherical symmetry, set $\varphi =e^{-i\omega t}F(r)$ and $k=\omega /c_0$. The radial equation is $(r^2{F}^{\prime }{)}^{\prime }+k^2{r}^{2}F=0$, or $(rF{)}^{\prime \prime }+k^2(rF)=0$. HenceRigid boundaries require radial velocity $F^{\prime }=0$ at both radii. A form automatically satisfying the condition at $a$ isAt $b$, the condition $b(rF{)}^{\prime }-(rF)=0$ becomesTherefore the nonzero oscillatory eigenfrequencies obeyThe undivided sine-cosine form avoids artificial issues from dividing by a cosine. The formal $k=0$ root is the spatially constant zero-frequency potential additive freedom of a velocity potential, not a propagating acoustic oscillation. This is the radial acoustic spectrum between rigid concentric spheres.

Use conservation of mass and Newton's momentum balance, in addition to the supplied constitutive law. For a material volume,The Reynolds transport theorem and ${\rho }_t+\nabla \cdot (\rho \mathbf{u})=0$ reduce the left side to $\int \rho D\mathbf{u}/Dt$. The divergence theorem then gives $\rho D{u}_i/Dt={\partial }_j{\sigma }_{ij}+\rho g_i$. For constant $\mu$,Incompressibility removes the last term. ThusThe stress formula and zero divergence alone are not a dynamical law: the force balance supplies the required physical input.

For a material volume $V(t)$ moving with velocity $\mathbf{u}$, the Reynolds transport theorem isIt applies componentwise to vector quantities. Cauchy's stress theorem represents traction on a surface with outward unit normal $\mathbf{n}$ as $t_i={\sigma }_{ij}{n}_j$. Conservation of mass gives $d{\int }_{V(t)}\rho dV/dt=0$ and hence ${\partial }_t\rho +\nabla \cdot (\rho \mathbf{u})=0$. Constant density reduces this to $\boxed{\nabla \cdot \mathbf{u}=0}$.

With no body force, conservation of momentum and the divergence theorem giveApplying transport and then mass conservation cancels the terms proportional to $u_i$, leaving $\rho ({\partial }_t{u}_i+u_j{\partial }_j{u}_i)={\partial }_j{\sigma }_{ij}$. For the stated Newtonian stress, constant viscosity gives ${\partial }_j{\sigma }_{ij}=-{\partial }_{i}p+\mu \Delta u_i+\mu {\partial }_i(\nabla \cdot \mathbf{u})$. The last term vanishes. Since the material volume was arbitrary,A body force, if present, would have to be added to the momentum balance; it is absent from the specified model.

The traction $\mathbf{t}(\mathbf{n})$ is force per unit area across an oriented surface with unit normal $\mathbf{n}$. The stress component ${\sigma }_{ij}$ is its $i$th component on a surface normal to the $j$th coordinate axis. Cauchy's stress theorem states that this dependence is linear:For the boundary of a material volume, take the outward normal and the traction exerted by the surrounding fluid. Newton's law and the divergence theorem giveConstant density and incompressibility make the transport-theorem integrand on the left $\rho D{u}_i/Dt$. The volume is arbitrary, soSubstitute the Newtonian constitutive law. With constant viscosity,Hence

Insert the traction formula into the material-volume balance and apply the divergence theorem to each row of the stress tensor:The volume is arbitrary, so continuity of the fields gives the local Cauchy momentum equationFor an incompressible Newtonian fluid, the constitutive law isHere the pressure enforces incompressibility, and the viscous stress is twice the viscosity times the symmetric rate-of-strain tensor. With constant $\mu$,Using $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ proves the stated Navier-Stokes equations.

If $L^{\prime }=L+dF(q,t)/dt$, the action changes by an endpoint term:In the principle of stationary action with fixed endpoints, $\delta q(t_i)=\delta q(t_f)=0$ and the endpoint times are fixed. Therefore $\delta (S^{\prime }-S)=0$, and the two actions give the same bulk Euler-Lagrange equations.

For a Lagrangian density, adding ${\partial }_{\mu }{K}^{\mu }$ changes the action by ${\int }_{\partial \Omega }{K}^{\mu }{n}_{\mu }d\Sigma$, by the divergence theorem. Its variation vanishes for compactly supported field variations, or when the appropriate boundary data are held fixed. These total divergences in a variational principle can therefore be discarded when deriving bulk equations under those conditions.

The qualification about boundaries matters: if $K^{\mu }$ depends on field derivatives, their variations on the boundary must also be controlled. A total derivative can change the canonical momenta or the chosen representative of a Noether current, even while leaving the bulk dynamics unchanged.

The Lorenz vector field has constant divergence $\nabla \cdot F=-\sigma -1-b$. Liouville formula for a fundamental matrix for the derivative of the flow therefore givesThis concerns transported three-dimensional phase volume, not the Euclidean distance between each pair of trajectories.

A quasiperiodic orbit with the usual smooth two-frequency torus closure would lie on an invariant embedded two-torus. That torus bounds a solid torus of positive finite volume in ${\mathbb{R}}^3$. Since the vector field is tangent to its boundary, its outward flux is zero, contradicting the divergence theorem, whose volume integral is $-(\sigma +1+b)V<0$. This negative-divergence obstruction to invariant tori excludes such quasiperiodic tori. The argument uses global negative divergence and the enclosed volume; saying only that an attracting torus has zero volume would not exclude it.

At an equilibrium the sum of eigenvalues is $-(\sigma +1+b)<0$, so they cannot all have positive real parts: no equilibrium is a fully repelling source. For a periodic orbit of period $T$, one Floquet multiplier equals one in the flow direction. The product of the other two is $e^{-(\sigma +1+b)T}<1$, so they cannot both have modulus greater than one: no periodic orbit is fully repelling either. Saddle equilibria and saddle periodic orbits remain possible; contraction of volume does not imply contraction in every direction.

Let $J$ be population flux. For any fixed control volume $V\subseteq \Omega$, mass conservation with local production givesAssume local, isotropic random dispersal, no imposed advection, and constant diffusion coefficient $D$. Fick's first law is $J=-D\nabla u$. The divergence theorem and arbitrariness of $V$ yieldThe coefficient has units length squared per time; $f$ has units population density per time. Smoothness sufficient to pass from the integral balance to the local equation is understood. This is a reaction-diffusion equation.

The Poincaré-Bendixson theorem states that the omega-limit set of a forward orbit confined to a compact region of the plane, containing no equilibrium points, is a periodic orbit.

Put $x=R\cos \theta$, $y=R\sin \theta$. Direct calculation givesHere $1/2\le H\le 17/16$ and $|\cos \theta {\sin }^3\theta |\le 3\sqrt{3}/16$. Thus on the inner circle $R=\sqrt{3}/2$ the radial velocity is positive, and on the outer circle $R=3/2$ it is negative. In their closed annulus,so there is no equilibrium. The theorem gives a periodic orbit in this trapping annulus.

It is asymptotically stable: the vector-field divergence isthroughout the annulus. Its nontrivial Floquet multiplier is $\exp {\int }_0^T\mathrm{div}Fdt<1$. This proves an attracting periodic orbit exists. There can be only one there: two such nested orbits would enclose an annular region with zero outward flux, contradicting the strictly negative divergence via the divergence theorem.

For an arbitrary material volume $V(t)$, Newton's second law equates the rate of change of linear momentum to the resultant contact force, since there is no body force:Use the Reynolds transport theorem on the left. Expanding the product and using the continuity equation leavesOn the right, the divergence theorem applied to each row of the Cauchy stress tensor givesConsequently the volume integral of $\rho D{u}_i/Dt-{\partial }_j{\sigma }_{ij}$ vanishes on every material volume. Localizing as in part (i) gives the Cauchy momentum equation:The left side is mass density times material acceleration; the right side is the divergence of the Cauchy stress tensor.

Let $V(t)$ have a boundary moving with velocity ${\mathbf{v}}_b$, and let $\mathbf{n}$ be its outward unit normal vector. For a smooth scalar field $q$, the Reynolds transport theorem isIn the last expression, extend ${\mathbf{v}}_b$ smoothly into $V(t)$ and use the divergence theorem. For a material volume, its boundary moves with the fluid velocity $\mathbf{u}$, so ${\mathbf{v}}_b=\mathbf{u}$. Equivalently, the Reynolds transport theorem for a material volume readswhere $D/Dt$ is the material derivative.

For an arbitrary moving region, mass conservation includes the relative mass flow rate across its boundary:Substitution of the Reynolds transport theorem and the divergence theorem givesSince the region is arbitrary, a continuous integrand must vanish: if it had a positive or negative value at a point, a sufficiently small surrounding region would have an integral of the same sign. This gives the continuity equation. With constant positive mass density $\rho$, it reduces toThus the flow is incompressible. For a material volume, the same argument starts directly from conservation of its total mass, with zero relative boundary flux.

For the upper curved half of $x^2/a^2+y^2/b^2+z^2/c^2=1$, take the vector field $(x{y}^2,y^2,C)$ and outward unit normal. Its divergence is $y^2+2y$. The odd term integrates to zero over the half-ellipsoid. Scaling a unit half-ball gives $\int y^{2}dV=2\pi a{b}^{3}c/15$. The equatorial closing disk has downward flux $-C\pi ab$. The divergence theorem therefore gives the displayed curved-surface flux. The constant vertical component is not detected by the divergence, but is recovered by accounting for the closing disk.

The Poynting theorem is the local energy balance of classical electromagnetism. For vacuum field laws, with current density $\mathbf{J}$, the electromagnetic energy per unit volume isandThus a decrease in local electromagnetic energy can supply outward energy flux or work on electric charge. The term $\mathbf{J}\cdot \mathbf{E}$ is power transferred to matter per unit volume. Its sign matters: positive values remove field energy, while negative values describe matter supplying the field. The magnetic part of the Lorentz force supplies no work term.

For a fixed volume $V$, the divergence theorem gives the directly measurable balanceFor example, $2.0\mathrm{W}$ entering a region with $0.50\mathrm{W}$ leaving and $1.20\mathrm{W}$ transferred to matter means stored electromagnetic energy increases at $0.30\mathrm{W}$. In source-free vacuum only the boundary flux changes the stored energy; a nonzero local Poynting vector can still have zero net closed-surface flux.

A brief proof takes the dot product of the Ampère-Maxwell equation with $\mathbf{E}/{\mu }_0$ and uses Faraday's law with $\mathbf{B}/{\mu }_0$. The identity $\nabla \cdot (\mathbf{E}\times \mathbf{B})=\mathbf{B}\cdot (\nabla \times \mathbf{E})-\mathbf{E}\cdot (\nabla \times \mathbf{B})$ turns the remaining terms into the divergence of the Poynting vector. ^[ref].

Let $-1<s<1$ and let $S_s$ be the outward-oriented boundary of $x^2+y^2+z^2\le 1$, $z\le s$. For a continuous function $h:[-1,s]\to \mathbb{R}$, the vector field $\mathbf{F}=h(z)(x,y,z{)}^T$ has outward fluxOn the curved unit sphere, the outward unit normal is $(x,y,z{)}^T$, so the normal component is $h(z)$. A horizontal spherical band has area $2\pi dz$, producing the first term. On the top disk, the outward unit normal is $(0,0,1{)}^T$, the normal component is $sh(s)$, and the area is $\pi (1-s^2)$, producing the second term.

For $h=f^{\prime }$ with a continuous first derivative, the fundamental theorem of calculus givesThis direct surface integral proof needs no second derivative of $f$. For $h=1$, it reduces to $\pi (2+3s-s^3)$; the divergence theorem identifies this with three times the enclosed volume.

The divergence theorem and Poisson equation give ${\int }_S\nabla \varphi \cdot d\mathbf{S}=4\pi Gm-2{\Omega }^{2}V$ on a bounding equipotential surface. Equipotential luminosity conservation with stellar circulation then fixes $f=L/(4\pi Gm-2{\Omega }^{2}V)$. Combining $dT/d\varphi =-f/\chi$, $dP/d\varphi =-\rho$ and radiative conductivity gives the displayed logarithmic temperature-pressure gradient. Both appearances of mass mean the enclosed mass $m$; only on the outer boundary can both be replaced by the total mass.

A vector field singular on an axis need not satisfy the classical divergence theorem on the original closed region. Remove a thin tube $\rho <\epsilon$, apply the theorem to the remaining smooth region, then examine the added cylindrical boundary. If its flux tends to zero and the remaining volume integral converges, the improper boundary flux equals the improper integral of the divergence. For $\mathbf{F}=(\log {\rho }^2,\log {\rho }^2,z/\rho )$ with $\rho =\sqrt{x^2+y^2}$, the cylindrical unit normal has zero $z$ component, so the added flux is $O(\epsilon |\log \epsilon |)$ on a bounded axial interval. The divergence $(2\cos \varphi +2\sin \varphi +1)/\rho$ is locally integrable under the cylindrical coordinates measure $\rho d\rho d\varphi dz$. Hence no finite line-source correction survives, despite the field itself being unbounded.

For a smoothly moving volume whose boundary has normal velocity $v_b\cdot n$, the derivative of an integral equals the local time-change integral plus the boundary-transport flux. The boundary term comes from the thin layer swept by the moving surface during a short time interval. Taking $f=1$ gives $d|V|/dt=\int v_b\cdot ndS$. If the boundary moves with the fluid, this combines with the divergence theorem to give the integral of the material derivative plus $f$ times the velocity divergence.

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.

A solid torus of major radius $A$ and minor radius $R<A$, centred on the horizontal plane, has volume $2{\pi }^{2}A{R}^2$. In circular torus coordinates, $z=r\cos \alpha$. Integrating $z^2$ with the Jacobian determinant $r(A+r\sin \alpha )$ gives ${\int }_V{z}^{2}dV=({\pi }^2/2)A{R}^4$; the part odd in $\sin \alpha$ vanishes. Thus the mean squared height is $R^2/4$. The identity is useful for moments of inertia and flux calculations using the divergence theorem.

Suppose two oriented surfaces have the same oriented boundary and $S_1\cup (-S_2)$ bounds a volume inside a domain where the vector field $\mathbf{F}$ is smooth and solenoidal. The divergence theorem gives zero flux through their closed union, proving the displayed equality. If a global vector potential $\mathbf{F}=\nabla \times \mathbf{G}$ exists on the relevant surfaces, the Stokes theorem proves the equality directly from their common boundary. The domain condition matters: a singularity inside the joining volume can obstruct the divergence argument.