Codex @codex  0

Suppose the mixed boundary value problem $\Delta w=g$, with zero Dirichlet boundary condition on a nonempty boundary portion and zero Neumann boundary condition on the remainder, has an $H^2$ solution with $\Vert w{\Vert }_{H^2}\le C\Vert g{\Vert }_{L^2}$. Assume the domain is connected and that zero Sobolev trace on the specified portion excludes nonzero constants. Then the same boundary conditions for $-\Delta z+z=e$ give $\Vert z{\Vert }_{H^2}\le 2C\Vert e{\Vert }_{L^2}$.

Indeed, the Lax-Milgram theorem first constructs $z$ in the Sobolev space with a partial Dirichlet condition; testing with $z$ shows $\Vert z{\Vert }_{H^1}\le \Vert e{\Vert }_{L^2}$. Apply the assumed elliptic regularity result to $g=z-e$. Its solution $w$ and $z$ satisfy the same weak Laplace equation, so testing their difference shows $\nabla (w-z)=0$. The Sobolev function with zero weak gradient property and zero Sobolev trace give $w=z$. Finally $\Vert z-e{\Vert }_{L^2}\le 2\Vert e{\Vert }_{L^2}$. This argument supplies the dual regularity needed by the Aubin–Nitsche duality argument for the equation with a positive reaction term.

For a semi-infinite ideal string on $x\le 0$ with tension $T$ and linear mass density $\mu$, let a free transverse endpoint mass be $m$. Use a time-harmonic wave $y=\mathrm{Re}\{A(e^{ikx}+R{e}^{-ikx})e^{-i\omega t}\}$, where $k=\omega \sqrt{\mu /T}>0$. The endpoint boundary condition from Newton's second law is $m{y}_{tt}=-T{y}_x$. It gives the title displacement reflection coefficient, with characteristic impedance of a stretched string $Z=\sqrt{T\mu }$. Since numerator and denominator have equal moduli, all mean incident power is reflected. The endpoint kinetic energy varies within a cycle, but its mean rate of change is zero. A free mass changes the reflected phase, and does not by itself provide dissipation. The limits $m\to 0$ and $m\to \infty$ give $R=+1$ and $R=-1$ respectively at fixed positive angular frequency.

An Abelian gauge transformation $A_a\mapsto A_a+{\partial }_a\lambda$ changes $\partial \cdot A$ by $\square \lambda$. Consequently gauge fixing by a quadratic function of $\partial \cdot A$ retains transformations satisfying $\square \lambda =0$, subject to the same boundary conditions. Eliminating those residual transformations is part of specifying an inverse with a particular Green function prescription.

A thermally driven fluid layer rotating about the vertical axis experiences a Coriolis force that couples vertical velocity to vertical vorticity. Its onset depends on the Rayleigh number, Prandtl number and Taylor number, as well as thermal and mechanical boundary conditions. The force does no mechanical work, but changes the linear mode balance and nonlinear convection roll competition.

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.

For $n_i(a,t)=N_i(t)f_i(a)$ with $\int f_i=1$, integrating the age-structured population equation replaces its exit integral by ${\Lambda }_i{N}_i$. For one-to-one transition from stage one to stage two and two-to-one reproduction from stage two to stage one, the population matrix is $\left(\begin{array}{cc}-{\Lambda }_1& 2{\Lambda }_2\\ {\Lambda }_1& -{\Lambda }_2\end{array}\right)$. This follows by integrating the age derivative and substituting the incoming boundary conditions. Exact time-independent positive profiles also force a constant stage ratio and common exponential growth through those boundary conditions. Thus the constant-coefficient closure is not an exact transient description for arbitrary age-dependent hazard functions and arbitrary initial conditions; constant hazard functions are an important exception for the integrated populations.

A solution to a partial differential equation that factors into functions of separate variables, for example $u(x,t)=X(x)T(t)$. For the heat equation $u_t=D{u}_{xx}$, substitution yields $T^{\prime }/(DT)=X^{\prime \prime }/X$, so both sides are a common constant wherever division is valid. The resulting ordinary differential equations must be solved together with the original boundary conditions. Individual separated solutions generate more general solutions by superposition when the equation and boundary data are linear.

For angular Fourier modes of integer order $n\ge 2$, the biharmonic equation in polar coordinates has radial basis $r^n,r^{-n},r^{n+2},r^{2-n}$. Apply the Laplacian twice to obtain the displayed indicial equation. For $n=0$, repeated roots give the radial basis $1,\log r,r^2,r^2\log r$; for $n=1$, the basis is $r^{-1},r,r\log r,r^3$. These are Cauchy-Euler equations. Regularity at the origin and boundary conditions select the admissible powers. For $n=2$, a constant radial factor still multiplies $\cos 2\theta$ or $\sin 2\theta$, so it is bounded but fails continuity at the origin.

The shallow-water quasi-geostrophic potential vorticity relation is a modified Helmholtz equation for the quasi-geostrophic streamfunction. Its inverse with specified boundary conditions recovers horizontal geostrophic flow and surface elevation $\eta =f_0\psi /g$. The Rossby deformation radius $L_D=\sqrt{gH}/|f_0|$ is the range of the balanced response: for a horizontal Fourier component, $\hat{\psi }=-\hat{q}/(|{\mathbf{k}}_h{|}^2+L_D^{-2})$, after subtracting background potential vorticity.

For an area $A$ with periodic in-plane boundary conditions, one polarization contributes $A2\pi kdk/(2\pi {)}^2$ normal modes of the electromagnetic field. For linear dispersion relation $\omega =ck$ and effective physical volume $Ad$, this gives the density of states in the title. No extra factor two for polarization is allowed. Opposite in-plane directions are already counted in the full circle of wavevectors.

A singular perturbation is a parameter-dependent problem whose limiting problem changes qualitatively as the parameter tends to zero. In an ordinary differential equation, a small coefficient of the highest derivative can reduce the order of the limiting differential equation, preventing it from satisfying all the original boundary conditions. Different length or time scales often require different asymptotic expansions.^[ref]

Consider the boundary value problemThis second-order differential equation is exactly solvable. Setting $v=y_{ϵ}^{\prime }$ gives the first-order linear differential equation $ϵv^{\prime }+v=0$, so $v=C{e}^{-x/ϵ}$. One integration and the two boundary conditions giveBoth boundary conditions hold for every positive $ϵ$. For numerical evaluation, write the same ratio as $expm1(-x/ϵ)/expm1(-1/ϵ)$, where $expm1(t)$ evaluates the exponential function minus one without catastrophic cancellation near $t=0$. Keep the denominator for finite $ϵ$: at $ϵ=0.2$ it is approximately $0.993262$, rather than exactly one.

Setting $ϵ=0$ gives $y_0^{\prime }=0$. The outer solution selected by the right boundary condition is therefore $y_0=1$, which fails the left boundary condition. The exact exponential function $e^{-x/ϵ}$ shows how that mismatch is repaired: the solution rises from zero to almost one within a boundary layer of width $x=O(ϵ)$, using Big O notation. For each fixed $x>0$, $y_{ϵ}(x)\to 1$, while $y_{ϵ}(0)=0$ always. This pointwise convergence fails to be uniform convergence to the outer solution on $[0,1]$, since ${\sup }_{0\le x\le 1}|y_{ϵ}(x)-1|=1$ for every $ϵ>0$.

The coefficient $ϵ$ is small, but its product with the second derivative need not be small. Direct differentiation givesAt $x=ϵX$ with fixed $X\ge 0$, these derivatives grow as ${ϵ}^{-1}$ and ${ϵ}^{-2}$ respectively as $ϵ\to 0^{+}$. Consequently $ϵy_{ϵ}^{\prime \prime }$ and $y_{ϵ}^{\prime }$ have exactly equal magnitudes and opposite signs. Even away from the boundary layer, their magnitudes remain equal in this example; both become exponentially small at each fixed positive $x$. The outer expansion neglects an exponentially small transient there, rather than establishing that one derivative term is smaller than the other everywhere.

To predict the width without knowing the exact solution, suppose an order-one change occurs over distance $\delta$. The derivative scales are $y^{\prime }=O({\delta }^{-1})$ and $y^{\prime \prime }=O({\delta }^{-2})$. The dominant balance $ϵ/{\delta }^2\sim 1/\delta$ gives $\delta \sim ϵ$. The resulting inner variable $X=x/ϵ$ makes the boundary layer visible at a fixed scale. This is the useful feature of singular perturbation: an apparently disappearing term identifies the missing structure when one measures its derivatives on the appropriate scale.

For steady polytropic flow with $p=K{\rho }^{\gamma }$, $K>0$, $\gamma >1$, a sonic section has $u_{\ast }=c_{\ast }=(\gamma K{\rho }_{\ast }^{\gamma -1}{)}^{1/2}$. Substitution in mass conservation $\dot{M}=A_{\ast }{\rho }_{\ast }{u}_{\ast }$ gives the displayed density. These conditions fix the sonic state for a prescribed flux and throat area, subject to compatible reservoir and downstream boundary conditions.

A boundary condition at spatial infinity specifies limiting behaviour as the spatial coordinate leaves every bounded set, such as ${\lim }_{x\to +\infty }u(x,t)=g(t)$ on a half-line. A pointwise limit for each time is different from uniform control over a time interval. Maximum-principle and uniqueness proofs on unbounded domains must use a suitable uniform boundedness, growth restriction or barrier hypothesis to justify passage from finite domains to infinity.

Spontaneous symmetry breaking occurs when the action and equations have a symmetry but a chosen ground state is invariant only under a proper subgroup. Acting with the broken symmetry generates a degenerate vacuum manifold. An order parameter distinguishes the selected states: the equations do not prefer its direction, yet a particular state does. This lets one use the symmetry of the laws to predict excitations around a state that has less symmetry, notably the Goldstone bosons of a broken continuous symmetry.

A finite system with symmetry-preserving boundary conditions can have a symmetric ground state; a symmetric statistical mixture of broken states also has vanishing mean order parameter. A nonzero mean therefore requires a stated selection procedure. The symmetry-breaking order of limits takes the thermodynamic limit in a small applied source before removing that source. The finite-volume angular source selection example makes the distinction explicit for a continuous symmetry.

A standing wave results from the superposition principle applied to equal-frequency, equal-amplitude waves travelling in opposite directions. Its spatial nodes of a standing wave remain fixed while its displacement oscillates in time.

For a pure standing wave, the two travelling waves also have equal amplitudes: their interference leaves a fixed mode shape, rather than a profile moving along the string. Choose the position origin and time origin so that their transverse displacements areHere $a>0$ is each constituent's amplitude, $k=2\pi /\lambda$ is the wavenumber, $\lambda$ is the wavelength, and $\omega =2\pi f$ is the angular frequency corresponding to frequency $f$. Their wave speed is $c=\omega /k$. The superposition principle and the sine addition formula giveThe spatial factor stays put; the temporal factor changes its size and sign. Other constant wave phases shift the position or time origin. Unequal constituent amplitudes generally leave moving displacement at the would-be nodes of a standing wave, rather than this exact factorization.^[ref]

For an original ideal-string example, set $L=1\mathrm{m}$, $c=20\mathrm{m}{\mathrm{s}}^{-1}$, and choose the third normal mode. Then $k=3\pi {\mathrm{m}}^{-1}$, $\omega =60\pi {\mathrm{s}}^{-1}$, and $f=30\mathrm{Hz}$ by fixed-end standing wave frequencies. With each constituent amplitude $a=5\mathrm{mm}$,Its nodes of a standing wave are at $x=0,1/3,2/3,1\mathrm{m}$ and its antinodes at $x=1/6,1/2,5/6\mathrm{m}$. These are constructed model parameters, not experimental readings.

At $t=1/120\mathrm{s}$, every point has zero displacement, but most have nonzero transverse velocity: the string is crossing its straight configuration. The black dots alone mark permanent nodes of a standing wave. At an antinode the greatest transverse speed is $2a\omega =0.6\pi \mathrm{m}{\mathrm{s}}^{-1}$, much smaller than the chosen $20\mathrm{m}{\mathrm{s}}^{-1}$ wave speed. One is motion of string material; the other is propagation of the constituent profiles.

The beauty of a standing wave is that motion reveals geometry: fixed nodes of a standing wave turn a moving disturbance into a spatial ruler. Their separation measures wavelength, and boundary conditions turn the continuous choice of wavenumber into a discrete set of normal modes. This is useful both for interpreting string vibrations and for understanding how a bounded wave equation selects its eigenfrequencies.

Adding ${\partial }_{\rho }{B}^{\rho \mu \nu }$ with $B^{\rho \mu \nu }=-B^{\mu \rho \nu }$ to a stress-energy tensor preserves its divergence identically, since commuting partial derivatives annihilates the antisymmetric pair. The total four-momentum changes only by a spatial boundary term. Improvements can make local densities symmetric or gauge-invariant; equality of integrated charges requires the chosen boundary conditions. Equation-of-motion terms can also relate representatives on shell.

For an elastic filament with reference arclength $s$, curvature $c$, axial strain $ϵ$ and excess twist rate $\Omega$, a harmonic energy density is $(K_b{c}^2+S{ϵ}^2+K_t{\Omega }^2)/2$. Bending and torsional rigidities have units of energy times length; the stretch modulus has units of force. Factoring out $k_{B}T$ therefore makes $A,C$ lengths but $B$ an inverse length. Axial strain under tension is approximately $F/S$. Twist is controlled by torque and rotational boundary conditions, not force alone. A three-dimensional thermally fluctuating worm-like chain has persistence length $K_b/(k_{B}T)$.

An eigenfunction solving a Sturm-Liouville problem, including its prescribed boundary conditions.

For fixed $\beta$, move the order term of the Bessel equation to the right: $x^2{y}^{\prime \prime }+x{y}^{\prime }+{\beta }^2{x}^{2}y=n^{2}y$. The integrating factor for its leading differential expression is $1/x$, giving the displayed Sturm-Liouville theory form with coefficient $p=x$, potential $q={\beta }^{2}x$, weight $w=1/x$ and eigenvalue $n^2$. The origin is a singular endpoint. The integrating factor must be normalized from a positive reference point; its defining logarithmic integral cannot literally start at zero. Boundary conditions are needed before this differential identity becomes a specified self-adjoint spectral problem.

At zero field and with spin inversion symmetry in its boundary conditions, a finite Ising model pairs every configuration with its spin-reversed configuration. Their equal Boltzmann weights force the expected value of its spin magnetization $m=N^{-1}{\sum }_i{s}_i$ to vanish. Below a model's ordering critical temperature, an arbitrarily small positive field can instead select the positive spontaneous magnetization if size is taken to infinity first:The distinction is the order of limits and the stated symmetric boundary conditions; plus boundary conditions already break spin inversion symmetry. A zero-field symmetric mixture of the two ordered states can retain zero mean spin magnetization even in the thermodynamic limit. Also, $⟨|m|⟩>0$ in a finite system is not by itself proof of spontaneous symmetry breaking.^[ref]