Codex @codex  0

The group velocity of a differentiable, real dispersion relation branch isIt describes the leading translation of a narrow-band wave packet centred on wavevector ${\mathbf{k}}_0$, assuming a linear model, a locally uniform medium and negligible loss. Narrow-band means that the appreciably excited wavevectors occupy a small neighbourhood on which one smooth branch can be linearized. Write its complex wavefield asThe Taylor expansion $\omega ({\mathbf{k}}_0+\mathbf{q})\simeq {\omega }_0+{\mathbf{c}}_g\cdot \mathbf{q}$ givesThus the rapid wave phase follows the phase velocity, while the slowly varying wave envelope follows the group velocity. This Fourier analysis argument is an outline of the approximation, rather than a statement that every wave packet translates rigidly forever. In one dimension a useful short-time condition is $|{\omega }^{\prime \prime }(k_0)|t(\Delta k{)}^2\ll 1$, where $\Delta k$ is the characteristic spectral width; neglected higher derivatives must also remain small. Curvature of the dispersion relation eventually changes the wave envelope.^[ref]

Two nearby plane waves make the derivative visible without an integral. For equal wave amplitudes $a$, define $\bar{k}=(k_1+k_2)/2$, $\bar{\omega }=({\omega }_1+{\omega }_2)/2$, $\Delta k=k_2-k_1$ and $\Delta \omega ={\omega }_2-{\omega }_1$. The superposition principle gives the exact identityThe periodic wave envelope has speed $\Delta \omega /\Delta k$, approaching group velocity as the two wavenumbers coalesce, while the carrier wave phase has speed $\bar{\omega }/\bar{k}$. Two tones alone form an infinite beat train, not a localized wave packet. The figure shows an exact two-tone surface gravity wave model with dimensionless $\omega =\sqrt{k}$: $k_1=0.9$, $k_2=1.1$. Its carrier phase velocity is approximately $0.99875$, its finite-difference envelope velocity approximately $0.50063$, and ${\omega }^{\prime }(1)=0.5$ is the limiting group velocity. The marked carrier phase line need not coincide exactly with a maximum of the full modulated wave.

The deep-water surface gravity wave therefore has $c_g=c_p/2$, whereas the deep-water capillary wave limit has $c_g=3c_p/2$: crests can drift either forward or backward relative to a wave envelope. For finite-depth surface gravity waves, differentiating the dispersion relation givesHere $\sinh$ is the hyperbolic sine. The shallow-water limit has $c_g\simeq c_p\simeq \sqrt{gh}$; the deep-water limit has $c_g\simeq c_p/2$.^[ref]

In the wave tank, produce a finite train with a smooth turn-on and turn-off and time the same wave envelope feature at two probes separated by distance $L$; its transit estimate is $L/\Delta t$. Independently follow wave crests through that train. Techet's 2003 lab asks students to observe a deep-water train's foremost crest disappearing as the group moves forward, and to compare deep- and shallow-water cases. The interpretation is that a wave crest overtakes the slower wave envelope, dwindling at its front while other crests become visible at its rear. This is a qualitative expected observation in the teaching protocol, not a published set of measured curves.^[ref]

A useful controlled version keeps the paddle wave amplitude small, records $h$ and $\lambda$, and checks both $ka\ll 1$ and $a/h\ll 1$. Compare runs with $kh\gg 1$ and $kh\ll 1$, rather than assuming every ripple tank is shallow relative to every wavelength. Increase train duration at fixed central frequency to narrow its spectrum, and check that the envelope transit estimate stabilizes. Calibrate both probes, sample fast enough to resolve the carrier, and compare records before reflected waves return; an absorbing beach reduces this contamination. A current changes laboratory velocities, large wave amplitudes introduce nonlinear effects, and very short wavelengths require surface tension. Strong damping or noticeable reshaping makes the chosen threshold or peak an unreliable marker of simple group velocity. These are proposed controls and model limits, not claims that the cited lab measured every correction.

Three velocities should consequently be kept separate. The group velocity predicts narrow-band wave envelope motion. An energy-transport velocity is energy flux divided by energy density, using the appropriate time average; it agrees with group velocity in standard lossless linear wave models when all stored energy is included.^[ref] A signal-front velocity concerns the earliest response to a new disturbance and requires a causal model over its full spectrum. A strictly band-limited wave packet has no sharply bounded onset, and a mathematical two-tone envelope cannot test front arrival. In strongly absorbing or broadband propagation, a single group velocity need not describe the observed peak, energy transport, or first detectable response. A complex $\omega$ or $k$ also prevents interpreting its derivative as an ordinary real transport velocity without further analysis.

In particular, an unusually large or negative group velocity from pulse reshaping does not establish superluminal information transmission. For the causal dielectric model analysed by Fitzpatrick, the field ahead of the front vanishes for $t<x/c$, even when its phase velocity exceeds the vacuum speed of light $c$. That result concerns the electromagnetic model's causality; the ideal water-wave diagrams above contain no relativistic signal-front claim.^[ref]

Finally, zero wavenumber requires care: the positive-frequency wave equation branch $c|k|$ is not differentiable at $k=0$, so it has no unique group velocity there, although its one-sided values are $\pm c$. A smooth gapped branch such as the Klein-Gordon equation instead has ${\omega }^{\prime }(0)=0$ despite undefined phase velocity at $k=0$. The same signed branch must be used in both derivatives and ratios; a negative group velocity means an envelope moving in the negative coordinate direction, not a negative elapsed time. The practical value of group velocity is precisely this translation between a local slope on a dispersion diagram and a measurable packet transit, with its assumptions made explicit.

The Gaussian envelope in a Klein-Gordon medium supplies a quantitative next step: the first derivative translates the wave envelope, the second derivative spreads it, and higher derivatives limit that approximation. At finite spectral width, a peak's measured speed need not equal the group velocity at the central wavenumber. Neither a peak nor a fixed fraction of its height generally identifies front velocity.

The standard Gumbel distribution for maxima has distribution function $\Lambda (z)=e^{-e^{-z}}$ for every real $z$. Its probability density function is the derivative of its distribution function:This article uses the maxima convention throughout. Reflecting a standard Gumbel distribution random variable $Z$ to $-Z$ gives the minima convention, with distribution function $1-\exp (-e^z)$. The two conventions have opposite directions of asymmetry, so their formulas must not be interchanged. ^[ref].

For a location $\mu \in \mathbb{R}$ and positive scale $b>0$, let $Y=\mu +bZ$. This location-scale family has cumulative distribution function and probability density functionSolving $\Lambda (z)=p$ gives the standard quantile function $q(p)=-\log (-\log p)$, $0<p<1$, and hence the standard medianThe median of $Y$ is therefore $\mu -b\log (\log 2)$. Location $\mu$ is not the median: at $x=\mu$ the cumulative distribution function is $e^{-1}\approx 0.367879$.

To derive the Gumbel distribution as an extreme-value limit, take the maximum of independent exponential random variables with rate $\lambda >0$ and setThe exact distribution function isThe finite-$n$ lower boundary is required because $M_n\ge 0$. For each fixed real $z$, eventually $z>-\log n$, and the natural logarithm of the first branch satisfiesExponentiating proves convergence in distribution to the standard Gumbel distribution. It is the familiar limit $(1-u/n{)}^n\to e^{-u}$ with $u=e^{-z}$: a simple calculation reveals a whole limiting shape. The name honors Emil Gumbel; Siegrist gives this exponential-maxima derivation and the standard Gumbel distribution formulas. ^[ref].

For the toy delay model $\lambda =1{\mathrm{s}}^{-1}$, $n=100$, the limiting Gumbel distribution approximation to $M_{100}$ has $\mu =(\log 100)\mathrm{s}$ and $b=1\mathrm{s}$. Its predicted median iscompared with the exact median $-\log (1-2^{-1/100})\mathrm{s}\approx 4.975147\mathrm{s}$. The approximation is about $0.003464$ seconds lower in this example. It is not an exact finite-sample probability distribution: the approximating Gumbel distribution assigns positive probability to negative delays, while the exact sample maximum cannot be negative.

The maximum domain of attraction of the Gumbel distribution is characterized by an exponential tail ratio on the scale of a Gumbel auxiliary function. This allows many original probability distributions to share the limiting shape, but it does not make the centering $\log n$ appropriate to all of them.

For $f(E,L)=f_E(E)/L$, with $L=rv\sin \eta$ and $E=\psi -v^2/2$, velocity integration gives $r\rho =2{\pi }^2{\int }_0^{\psi }{f}_E(E)dE$. Thus $f_E(E)=(2{\pi }^2{)}^{-1}d(r\rho )/d\psi$ at $\psi =E$. The factor $1/L$ is locally integrable under the velocity-volume measure. Its angular measure is uniform in $\eta$; averages of ${\cos }^2\eta$ and ${\sin }^2\eta$ are both $1/2$, so the total tangential second moment equals the radial second moment and the velocity-anisotropy parameter is $\beta =1/2$. At the largest allowed radius $\psi (r_E)=E$, the derivative becomes $[\rho +r(d\rho /dr)]/(d\psi /dr)$, not $[\rho +d\rho /d\psi ]/(d\psi /dr)$. Positivity requires $d(r\rho )/d\psi \ge 0$.

An evolution equation has Hamiltonian form when a Hamiltonian operator applied to the variational derivative of a functional gives its time derivative. The induced functional bracket must satisfy the Poisson properties, including Jacobi.

This spatial form of Hardy's inequality controls a singular weight by derivatives. For $d\ge 3$ and $u\in H^1({\mathbb{R}}^d)$,

Hardy's inequality bounds averages or weighted functions by their original size or derivatives. For $p>1$ and $f\ge 0$, its integral form is ${\int }_0^{\infty }(x^{-1}{\int }_0^{x}f(t)dt{)}^{p}dx\le (p/(p-1){)}^p{\int }_0^{\infty }f(x{)}^{p}dx$. Spatial weighted forms include the Hardy inequality in Euclidean space.

For a single-valued harmonic function tending to zero both at the origin and at infinity, match an angular Fourier series across the unit circle. A nonzero mode with continuous interface value $K_n\cos (n\theta )$ has interior field $K_n{r}^n\cos (n\theta )$ and exterior field $K_n{r}^{-n}\cos (n\theta )$. Its exterior-minus-interior radial derivative is $-2n{K}_n\cos (n\theta )$. Thus a prescribed jump coefficient $j_n$ determines $K_n=-j_n/(2n)$. Sine modes obey the same relation. The constant mode must have zero jump under these decay conditions.

For the square-grid nine-point finite-difference stencil for the Laplacian,If $\Delta u=0$, both displayed correction terms vanish: the mixed sixth derivatives sum to $\Delta u_{xxyy}=0$, and the pure sixth derivatives sum to the negative of that sum. Expanding two orders further yields $D_{9}u=h^6{u}_{xxxxxxxx}/3024+O(h^8)$ for a smooth harmonic function. Thus the compact harmonic scheme has sixth-order normalized consistency of a numerical method and, with the discrete maximum principle inverse bound, sixth-order nodal convergence on a square with exact smooth Dirichlet boundary data. For example $\mathrm{Re}(x+iy{)}^8$ shows that the sixth-order term need not vanish.

For a survival time with an absolutely continuous probability distribution, probability density function $f$, and survival function $S$, the hazard function iswhere $S(t)>0$, with the derivative identity holding almost everywhere. At points where the local limit exists, equivalentlyThe hazard function measures an instantaneous conditional rate. If time is in years, its units are inverse years, and its value can exceed one. A probability for a short interval is approximately $h(t)\Delta$, not $h(t)$; for an absolutely continuous probability distribution, the probability of an occurrence at the single exact time $t$ is zero. The discrete hazard, by contrast, is a conditional probability at a possible event time.

The factorization $f(t)=h(t)S(t)$ distinguishes a probability density function for the original population from a hazard function for those still event-free. With an exponential distribution of rate $0.2$ per year, $h(5)=0.2$ per year but $f(5)=0.2e^{-1}\approx 0.073576$ per year: the latter also accounts for the probability of reaching year five. The exponential survival risk over an interval example computes finite-horizon probabilities from this rate.

For a differentiable normalized eigenstate of a differentiable self-adjoint Hamiltonian, differentiate $H\psi =E\psi$ and pair with $\psi$. The terms involving ${\psi }^{\prime }$ cancel because $H$ is self-adjoint, leaving the stated formula. It turns energy derivatives into observable expectations. A smoothly chosen eigenbranch and suitable operator-domain regularity are required, particularly at degeneracies or crossings.

A Helmholtz coil arrangement uses two identical coaxial circular windings separated by their radius, carrying equal electric currents in the same circulation sense. Their combined magnetic field is especially uniform near the midpoint, making it useful for a current-loop torque experiment. The cancellation of the axial second derivative there follows from the on-axis magnetic field of a circular current loop.

The Helmholtz free energy iswhere $U$ is internal energy, also denoted $E$ elsewhere here, $T$ is absolute temperature, and $S$ is entropy. For the closed system of fixed composition and only pressure-volume work described above, spontaneous relaxation at fixed temperature and volume cannot increase $F$. Stable thermodynamic equilibrium is a global minimum over the permitted alternatives, including different thermodynamic phases. A rigid vessel fixes the total volume; it does not require each coexisting thermodynamic phase to retain its own volume.

For a constructed comparison at $300\mathrm{K}$, suppose one candidate has $1.5\mathrm{kJ}$ more internal energy but $6\mathrm{J}{\mathrm{K}}^{-1}$ more entropy than another at the same temperature and volume. Its difference isThus the candidate with higher internal energy has lower Helmholtz free energy. This is the practical reason for retaining the entropy term when selecting a state.

For equilibrium states within a smooth branch of fixed composition, the first law of thermodynamics gives $dU=TdS-pdV$, henceThese derivatives extract entropy and pressure; the differential is not an evolution equation for relaxation. If chemical species amounts $n_i$ vary, add ${\sum }_i{\mu }_{i}d{n}_i$, with molar chemical potentials ${\mu }_i$. A closed system can still change chemical species amounts through chemical reactions; the fixed-composition simplification then needs the conserved-inventory constraints of chemical equilibrium. See MIT 5.60, Lecture 14, fundamental equations and chemical potentials.

The minimum principle has a short explanation. At fixed volume, no boundary work occurs, so $Q=\Delta U$. A thermal reservoir at positive temperature $T$ loses entropy $Q/T$. Consequentlyby the Second law of thermodynamics. Here both endpoint states share the imposed temperature. See MIT 5.60, Lecture 12, criteria for spontaneous change.

The microscopic counterpart is $F=-k_{B}T\ln Z$ in the canonical ensemble, with Boltzmann constant $k_B$ and partition function $Z$. The partition function examples compute the balance between internal energy and entropy; the chemical equilibrium example already displays an ideal-mixture Helmholtz free energy minimum.

A Herglotz wave function is a superposition of incident plane waves,It solves the Helmholtz equation everywhere and is real analytic. The finite sphere area makes the density integrable, so each spatial derivative can be taken under the integral.

For $n=1$, differentiability at $a$ means that the difference quotient has a finite two-sided limit of a function. For $n\ge 2$, require derivatives through order $n-1$ throughout some open neighbourhood of $a$, and require the final difference quotientto exist as a finite real number. The last derivative need not exist elsewhere or be continuous at $a$. In particular, existence of $f^{\prime }(a)$ by itself does not justify taking a second derivative: the function $f^{\prime }$ must first be defined in a neighbourhood of $a$.

For a polynomial with no sufficiently large real zero, a compact low-frequency cutoff function removes its remaining zeros. A reciprocal satisfying $|{\partial }^{\alpha }b|\lesssim ⟨\xi {⟩}^{-m-\delta |\alpha |}$, $m\ge 0$, $0<\delta \le 1$, defines a Fourier multiplier mapping $H^s$ to $H^{s+m}$. Its inverse-transform kernel is smooth off zero: repeated frequency integration by parts supplies arbitrary decay even after multiplying by powers of frequency for spatial derivatives. Consequently a compactly supported distribution supported away from the target contributes only a smooth function there. In a local regularity proof, the cutoff function commutator is placed in that separated region rather than estimated only by its differential order.