Codex @codex  0

A smooth phase gives ${\partial }_{X_i}{k}_j={\partial }_{X_j}{k}_i$ and ${\mathbf{k}}_T=-\nabla \omega$. These mixed-derivative identities allow the spatially differentiated dispersion relation to be rewritten as Hamiltonian ray equations for a local dispersion relation. Away from phase singularities they express the curl-free character of a local wavevector.

If a differentiable dispersion relation $\omega (\mathbf{k})$ is homogeneous of degree zero at nonzero $\mathbf{k}$, the Euler homogeneous function theorem gives $\mathbf{k}\cdot {\nabla }_{\mathbf{k}}\omega =0$. Therefore phase velocity $\omega \mathbf{k}/k^2$ is orthogonal to group velocity ${\nabla }_{\mathbf{k}}\omega$. The result applies exactly to limiting inertial waves; finite compressibility can alter the dispersion.

A resonant interaction requires matching both wavenumber and angular frequency: $k_3=k_1+k_2$ and ${\omega }_3={\omega }_1+{\omega }_2$ for a sum interaction, with corresponding sign changes for differences. A common phase velocity alone does not guarantee that a generated harmonic lies on the dispersion relation.

A monochromatic internal gravity wave in uniform stratification has a spatially linear phase. In the nonrotating inviscid Boussinesq approximation, its dispersion relation is ${\omega }^2=N^2|{\mathbf{k}}_h{|}^2/(|{\mathbf{k}}_h{|}^2+m^2)$. The wave vector is normal to its phase planes, while the group velocity lies along them.

For a closed string in a string sigma model in an isotropic plane wave, insert $e^{-i{\omega }_n\tau +2\pi in\sigma /L}$ into the transverse wave equation. Periodicity requires integer $n$, giving the displayed dispersion relation. The zero spatial mode is an oscillator of frequency $|m\kappa |$; when this mass vanishes it becomes a free particle mode and needs the position-plus-linear-time expansion.

A harmonic line force localized at one point of a two-dimensional sheet gives $\hat{\eta }=F/D$ and $\Delta =\gamma D$ in the Fourier transform. The outgoing upper pressure is the displayed amplitude. Dividing by $c_0^2$ converts it to acoustic density perturbation. The method of steepest descent gives cylindrical sound from the saddle point, while crossed zeros of the dispersion relation give guided or leaky modal residues. Retaining $k_0=\omega /c_0$ is essential for dimensionally consistent spreading and prefactors.

A positive jump $\Delta Q$ of shallow-water quasi-geostrophic potential vorticity supports a current $U=U_0{e}^{-|y|/L_R}$ with $U_0=L_R\Delta Q/2$. A disturbance has structure $\varphi =A{e}^{-\sqrt{k^2+L_R^{-2}}|y|}$. Integrating the linear potential-vorticity conservation equation across the interface gives $(U_0-c)[{\varphi }_y]+\Delta Q\varphi (0)=0$, hence the displayed dispersion relation. The wave propagates intrinsically against the mean current because the displaced PV interface induces the velocity that advects it; no planetary beta is required.

The high-pressure limit of the slow magnetosonic wave approaches an incompressible polarization in the plane of the background magnetic field and wavevector, perpendicular to the latter. It has the displayed leading dispersion relation, the same as the out-of-plane Alfvén wave. Magnetic tension is the restoring force, and pressure balances the compressive component of the magnetic force. In exactly incompressible homogeneous ideal magnetohydrodynamics it is one of the two transverse magnetic polarizations.

A radio wave traversing a dilute cold plasma has a group velocity that depends on its frequency. In the ideal unmagnetized limit, the dispersion relation is ${\omega }^2=c^2{k}^2+{\omega }_p^2$, where ${\omega }_p$ is the angular plasma frequency. For $\omega \gg {\omega }_p$, differentiate to obtain $v_g=c\sqrt{1-{\omega }_p^2/{\omega }^2}$, then expand $1/v_g$ to first order. The excess travel time over a vacuum path isHere $f=\omega /(2\pi )$ is ordinary frequency, $e$ is the magnitude of the elementary charge, $m_e$ is the electron mass, ${ϵ}_0$ is the vacuum permittivity, $n_e$ is electron number density, and $c$ is the speed of light. This approximation neglects magnetic corrections, collisions and multipath scattering; its positive sign means lower frequencies arrive later.

A quadratic product of plane waves generates sums and differences of their wavenumbers and angular frequencies. A component drives a secular term when its generated pair also satisfies the linear dispersion relation. A solvability condition in the method of multiple scales then determines the coupled amplitude equations on a slow time.

For the whole-line equation $A_t+a_1{A}_{xx}+a_2{A}_{xxxx}+a_{3}A=0$ with $\mathrm{Re}{a}_2>0$, the dispersion relation is $s=a_1{k}^2-a_2{k}^4-a_3$. The temporal Laplace transform of the impulse response at zero is the real-k integral of $[s+a_3-a_1{k}^2+a_2{k}^4{]}^{-1}/(2\pi )$. Factor the denominator as $a_2(k^2+{\kappa }_1^2)(k^2+{\kappa }_2^2)$ with $\mathrm{Re}{\kappa }_j>0$ for large $\mathrm{Re}s$. Partial fractions give the displayed resolvent. Square-root sheets are chosen by continuation from this causal large-s region; choosing roots independently can introduce false spatial pinch points.

For $V=m^2{\varphi }^2/2+\lambda ({\varphi }^2{)}^2/4$, $m^2<0$, the vacuum manifold is a sphere of radius $v=\sqrt{-m^2/\lambda }$. At a chosen vacuum the Hessian matrix is $2\lambda v^{2}n{n}^T$, where $n$ is its unit direction. The radial scalar has mass squared $2\lambda v^2$; the $n-1$ angular modes are physical Goldstone bosons when the symmetry is global.

For the two-component case, write $\mathbf{\varphi }=(v+\rho )(\cos \theta ,\sin \theta )$. The kinetic term becomes $\frac{1}{2}({\partial }_{\mu }\rho {)}^2+\frac{1}{2}(v+\rho {)}^2({\partial }_{\mu }\theta {)}^2$. At quadratic order the canonically normalized angular field is $\pi =v\theta$, with no mass term. Its Euler-Lagrange equation is the wave equation, so ${\omega }^2=|\mathbf{k}{|}^2$ in units where the propagation speed is one. The radial field instead obeys ${\omega }^2=|\mathbf{k}{|}^2+2\lambda v^2$. Thus a common rotation costs no energy, whereas a slowly varying rotation propagates as a collective wave. With the selecting source of continuous vacuum selection by an external source, use $\pi =w\theta$ and obtain ${\omega }^2=|\mathbf{k}{|}^2+h/w$. This elementary dispersion relation separates a Goldstone boson, a pseudo-Goldstone boson and a massive radial excitation within one example.

For tightly wound perturbations in a locally uniform, razor-thin gas disk, continuity equation and rotating-fluid momentum equations couple the surface-density perturbation to radial and azimuthal velocities. The sheet potential is $\delta \Phi =-2\pi G\delta \Sigma /|k|$. Eliminating these perturbations gives the displayed dispersion relation. radial epicyclic frequency supplies the rotation term, the sound-speed term resists short waves, and self-gravity supplies the negative term. For an axisymmetric mode $n=0$, use ${\omega }^2$ directly; the pattern speed $\omega /n$ is undefined.

A free acoustic membrane wave is a zero of the outgoing dispersion relation $D$. At this zero, the impedance seen from one side equals the negative of that side's outgoing normal acoustic impedance. Thus the continued by analytic continuation reflection coefficient has a pole. A real guided mode is evanescent normal to the sheet; a complex root can describe a leaky mode. The pole describes nontrivial homogeneous motion with no incoming wave, rather than unbounded passive reflection at an ordinary real incidence angle.

For a one-dimensional free particle with $m>0$ and $p>0$, take the positive total energy from the relativistic energy-momentum relation:A fixed plane wave quantum phase obeys $kx-\omega t=\text{constant}$. Its phase velocity is thereforeFor a narrow wave packet centered at $p_0$, expand its dispersion relation as $E(p)\simeq E(p_0)+E^{\prime }(p_0)(p-p_0)$. In the Fourier inversion integral building the wave packet, the envelope then depends on $x-E^{\prime }(p_0)t$. Thus its group velocity isThe superluminal phase velocity follows a repeating phase marker; it does not describe a localized signal or superluminal transport of energy. Higher derivatives of the dispersion relation produce wave packet spreading, so the translated-envelope argument is a narrow-band approximation. For $v=0.6c$, the Lorentz factor is $1.25$, $p=0.75mc$ and $E=1.25m{c}^2$: the group velocity is $0.6c$ while the phase velocity is $5c/3$.

The total-energy convention matters. Define $\psi =e^{im{c}^{2}t/\hslash }\Psi$ to remove the common rest energy from a wavefunction $\Psi$ with total-energy time dependence. ThenThis is exactly the free-particle momentum convention used by the Schrödinger equation here. Removing a common time-dependent global phase changes the reported phase velocity, while leaving the de Broglie wavelength, relative spatial quantum phases and group velocity unchanged. Adding $m{c}^2$ again to a total energy would count the rest energy twice.

A layer with shear modulus ${\mu }_1$ and wave speed $c_1$ lies between a rigid support and a half-space with shear modulus ${\mu }_2$ and wave speed $c_2>c_1$. A trapped antiplane shear normal mode has layer amplitude $A\sin (my)$ and exterior amplitude $A\sin (mh)e^{-\ell (y-h)}$, where $m=k\sqrt{c^2/c_1^2-1}$ and $\ell =k\sqrt{1-c^2/c_2^2}$. Continuity of displacement field and traction gives the displayed dispersion relation. It requires $c_1<c<c_2$ and differs from a Love wave with a free boundary.

A roton is an excitation near a local minimum of a superfluid's dispersion relation at nonzero momentum $p_0$. The displayed local approximation defines a positive energy gap $\Delta$ and positive curvature effective mass $m_r$; it is not a free helium atom or a macroscopic quantum vortex.^[ref] In the stipulated plot, these correspond to $p_0=2p_{\ast }$, $\Delta =0.6E_{\ast }$ and $m_r=p_{\ast }^2/(4E_{\ast })$. The plot illustrates how a low energy at appreciable momentum can lower the Landau criterion below the speed of sound, without assigning any measured parameter to real helium.

Landau's 1947 revision joined the low-momentum linear part to a spectrum with a finite-momentum minimum, retaining the macroscopic superfluid theory while changing its excitation model.^[ref]

On a ray $x/t=V$, a localized wave packet is governed by a physically selected saddle point satisfying $d\omega /dk=V$ in its analytic dispersion relation. The exponent's real growth is $\mathrm{Im}({\omega }_{\star }-V{k}_{\star })$. The fixed-frame case $V=0$ gives the absolute growth rate. This ray formulation makes the frame dependence of convective hydrodynamic instability explicit.

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.

Consider an original ideal one-dimensional monatomic harmonic oscillator chain with neighbor separation $a$, mass $M$ and nearest-neighbor restoring coefficient $K>0$. Its displacement equation isSubstitute the normal mode $u_r=u_0{e}^{i(qra-\omega t)}$. Since $e^{iqa}+e^{-iqa}-2=-4{\sin }^2(qa/2)$,Choose $a=4$ ångströms and ${\epsilon }_{\ast }=10\mathrm{meV}$. This is one acoustic phonon branch, distinct from the two-atom elastic example. It has $\epsilon (\pi /(2a))=5\sqrt{2}\mathrm{meV}$ and $\epsilon (\pi /a)=10\mathrm{meV}$.

For coherent creation of one phonon in a crystal, the excitation's wave vector $q$ and the measured transfer obey $Q=G+q$, with $G$ a reciprocal lattice vector. To connect the dispersion relation to an illustrative intensity, take zero temperature, longitudinal displacements and the first interval $Q=q\in (0,\pi /a]$, choosing $G=0$. Expand the scattering phase to first order, $e^{iQ(ra+u_r)}\simeq e^{iQra}(1+iQ{u}_r)$. The constant term supplies elastic scattering; the linear term creates one phonon. Its displacement matrix element is proportional to $1/\sqrt{\omega }$, so the squared amplitude has weight $Q^2/\epsilon (Q)$, apart from common factors. Normalize this weight to one at $Q=\pi /a$:The limit $w(0)=0$ follows from $\sin (\pi x/2)\sim \pi x/2$. Set $\sigma =0.4\mathrm{meV}$ and replace the ideal sharp line by a normalized Gaussian function:This chosen broadening illustrates readout width; it is not a measured resolution or a full spectrometer resolution function. Integrating over all $\Delta E$ gives $w(x)$. The small negative-energy tail is a broadening artifact, not a zero-temperature phonon-annihilation signal.

These panels separate three kinds of information: allowed reciprocal lattice positions, cancellation within a crystal basis structure factor, and an excitation's energy as a function of wave vector. The chain has no magnetic degrees of freedom, optical branches, damping or multi-phonon processes. It illustrates how one excitation could be read; it proves nothing about the completeness of a real material's excitation spectrum.