Codex @codex  0

For an expanding pressureless matter universe with zero cosmological constant, fixed current Hubble parameter $H_0$ and current matter cosmological density parameter $\Omega >0$, the age of a dust universe without a cosmological constant is $H_0{t}_0={\int }_0^1\sqrt{x}/\sqrt{\Omega +(1-\Omega )x}dx$. For every $0<x<1$ the integrand decreases strictly as $\Omega$ increases. Thus open models with $0<\Omega <1$ are older than the flat model, and closed expanding models with $\Omega >1$ are younger. The flat value is $2/3$ and the empty Milne universe limit is $1$. Specifying what is held fixed matters: at fixed present scale factor and energy density, the integral $t_0={\int }_0^{a_0}da/\sqrt{C/a-k}$ instead increases with $k$, provided all three models reach $a_0$ on their expanding branches.

Black-body radiation is the electromagnetic radiation in thermal equilibrium with an ideal blackbody. It can be described as a photon gas: a cavity whose walls are all at the same temperature contains radiation with the universal Planck's law spectrum. A small opening in such a cavity approximates a blackbody emitter because incoming radiation undergoes repeated absorption inside. The radiation's equilibrium spectrum depends on temperature, rather than the wall material. ^[ref]

There are three useful quantities to keep distinct: directional spectral radiance $B_{\nu }$ or $B_{\lambda }$, total energy density $u$ in the cavity, and outward radiant exitance $j^{\ast }$ from its opening. Planck's law supplies the first; integrals over frequency and directions supply the others. For isotropy, every direction has the same $B_{\nu }$, so the full $4\pi$ solid angle givesHere $\hslash =h/(2\pi )$ is the reduced Planck constant, and the radiation constant is $a_{\mathrm{rad}}\approx 7.565733250\times 10^{-16}\mathrm{J}{\mathrm{m}}^{-3}{\mathrm{K}}^{-4}$. The explicit Boltzmann constant makes this an SI formula; adopting natural units with $k_B=1$ recovers $u={\pi }^2{T}^4/(15{\hslash }^3{c}^3)$.

For a short integral outline, substitute $x=h\nu /(k_{B}T)$ and useThe expansion is a geometric series with nonnegative terms; repeated integration by parts evaluates each integral, and $\zeta (4)={\pi }^4/90$ is a value of the Riemann zeta function. This gives the coefficient above.

The photon gas also exerts radiation pressure, $p=u/3$. The factor $1/3$ comes from the angular average of the squared directional component of photon momentum, $⟨{\cos }^2\theta ⟩=1/3$. The thermodynamic derivation of blackbody radiation pressure gives the same result. In contrast, outgoing power involves one projected-direction factor over an outward hemisphere, yielding radiant exitanceThis is the Stefan–Boltzmann law with unit emissivity; $\sigma$ is the Stefan-Boltzmann constant. The factor $\pi$ follows from Lambert's cosine law, rather than the $4\pi$ full-sphere factor used for energy density.

For example, an ideal $6000\mathrm{K}$ cavity has $u\approx 0.980519\mathrm{J}{\mathrm{m}}^{-3}$, radiation pressure $p\approx 0.326840\mathrm{Pa}$, and outward radiant exitance $j^{\ast }\approx 7.34881\times 10^7\mathrm{W}{\mathrm{m}}^{-2}$. A stipulated $3000\mathrm{K}$ cavity has one-sixteenth each value because all three scale as $T^4$. These describe stored energy, mechanical pressure, and outgoing power per area, respectively. In a closed equilibrium cavity, equal opposing radiation balances the directional power flow even though $u$ and $p$ are nonzero.

The usefulness of black-body radiation is this common reference: Planck's law, Wien's displacement law, and the Stefan–Boltzmann law consistently connect spectral measurements, temperature, and total emission. The three illustrative curves in Planck's law are ideal equilibrium spectra; applying them to a material requires considering its emissivity and whether thermal equilibrium is a suitable approximation.

The capacitor energy is the electrostatic energy stored when electric charge is separated. For a linear capacitor of fixed positive capacitance $C$, plate electric charge magnitude $Q$, and terminal electric potential difference $V=Q/C$, it isThe zero is the uncharged state. Thus the capacitor energy is positive for nonzero $Q$, even though the total electric charge of the two plates is zero. ^[ref].

For the parallel-plate capacitor example, capacitance $C=8.854\times 10^{-11}\mathrm{F}$ and terminal electric potential difference $V=10\mathrm{V}$ giveThe plate electric charge magnitude is $Q=0.8854\mathrm{nC}$. Doubling terminal electric potential difference with capacitance unchanged doubles the plate electric charge magnitude and quadruples the capacitor energy.

The electric charge is a linear function of terminal electric potential difference, while capacitor energy is a quadratic function. The second curve in each plot uses twice the capacitance: at a given terminal electric potential difference, both plate electric charge magnitude and capacitor energy double. These relations let us calculate how much energy an electrical circuit can receive when a capacitor discharges.

A brief charging argument explains the factor $1/2$. If the plate electric charge magnitude has reached $q$, the terminal electric potential difference is $q/C$. Transferring a further amount $dq$ of electric charge quasistatically requires work $dW=(q/C)dq$. With the geometry and permittivity fixed, the integral isUsing $Q=CV$ gives the other two forms. This integral accounts for the gradually rising terminal electric potential difference; it is the work stored in the capacitor's electric field, rather than a claim that all energy drawn from a fixed-voltage battery during charging ends up stored.

The field picture gives the same result. In the vacuum gap of a parallel-plate capacitor, the electric field has energy density $u={ϵ}_0{E}^2/2$. Multiplying by the gap volume $Ad$ gives the electrostatic energyFor the worked example, $u=4.427\times 10^{-4}J/{\mathrm{m}}^3$ and $Ad=10^{-5}{\mathrm{m}}^3$ reproduce $U=4.427\mathrm{nJ}$. The electrostatic energy resides in the electric field throughout the gap, despite overall charge neutrality.

When capacitance changes, specify the constraint. At fixed plate electric charge magnitude $Q$, increasing capacitance decreases capacitor energy $Q^2/(2C)$; doubling capacitance halves it. At fixed terminal electric potential difference $V$, increasing capacitance increases capacitor energy $C{V}^2/2$; doubling capacitance doubles it. In the latter case a connected battery supplies additional electric charge and energy, and there may also be mechanical work. Use the fixed-voltage capacitor energy balance to account for all transfers before inferring a force from changing capacitor energy.

For a perfect fluid in general relativity with constant equation-of-state parameter $w$, use conformal time $dt=ad\tau$ and the conformal Hubble parameter $\mathcal{H}=a^{\prime }/a=aH$. The Friedmann acceleration equation gives ${\mathcal{H}}^{\prime }=a\ddot{a}=-4\pi G{a}^2(1+3w)\rho /3$. The Friedmann equation eliminates $\rho$ and gives this Riccati equation. For positive energy density, $\alpha \ne 0$ and appropriately shifted conformal time, its usual branches are $\cot (\alpha \tau )$ for $k=1$, $1/(\alpha \tau )$ for $k=0$, and $\coth (\alpha \tau )$ for $k=-1$. Their scale factors are proportional to ${\sin }^{1/\alpha }(\alpha \tau )$, $(\alpha \tau {)}^{1/\alpha }$ and ${\sinh }^{1/\alpha }(\alpha \tau )$ respectively on intervals where the bases are positive. At $\alpha =0$ the conformal Hubble parameter is constant and the separate solution is $a=A{e}^{h\tau }$.

Nonrelativistic matter dilutes as $a^{-3}$, while a cosmological constant has constant energy density. The ratio of vacuum energy density to matter density consequently scales as $a^3$, so comparable contributions occur during a limited epoch. The cosmic coincidence problem asks why observers find these very different components comparable near their own epoch. This timing puzzle is distinct from the cosmological constant problem, which concerns the tiny absolute vacuum-energy scale.

For an expanding Friedmann-Lemaître-Robertson-Walker metric with positive energy density, put $\Omega =8\pi G\rho /(3H^2)$ and $w=P/\rho$, in units $c=1$. The Friedmann equation, Friedmann acceleration equation and cosmological perfect-fluid continuity equation giveFor constant $w$, linearizing at $\Omega =1$ gives $\Omega -1\propto a^{1+3w}$. For variable $w$, replace the power by $\exp \int (1+3w(a))d\log a$. This describes the flatness problem for decelerating expansion, away from epochs at which $H=0$.

For one massless neutrino flavor with one helicity per member and $|{\mu }_{\nu }|\gg T_{\nu }$, the Fermi-Dirac distribution fills momentum states up to the positive chemical potential. If ${\mu }_{\nu }>0$, the neutrino density is ${\mu }_{\nu }^4/(8{\pi }^2)$ and the antineutrino contribution is approximately $3T_{\nu }^4{e}^{-{\mu }_{\nu }/T_{\nu }}/{\pi }^2$. For negative chemical potential the roles reverse. Thus the leading total energy density is even in ${\mu }_{\nu }$ and has no factor of two from the two species. Extra populated internal degrees of freedom do multiply it. This distinction prevents using a positive-degeneracy approximation for an exponentially dilute negative-chemical-potential distribution.

A dielectric is a material in which an electric field produces electric polarization while transport of electric charge is sufficiently small for the intended use. An insulating solid in a capacitor is an example; real dielectrics can still have dielectric loss and nonzero electrical conductivity.

The important distinction is between rearranging bound electric charge and transporting mobile charge carriers through the specimen. Bound electric charges can move a little within a chemical atom or molecule, or an existing electric dipole can turn. Their averaged electric dipole moments produce electric polarization even if there is negligible steady electric current between the electrodes. A changing electric polarization does give a bound current; the absence of steady conduction does not mean that the material is electrically unresponsive.^[ref]

For uniform electric polarization $\mathbf{P}=P\hat{\mathbf{x}}$ in a rectangular slab, polarization density gives zero bound volume charge density in the interior and opposite bound surface charge densities $-P,+P$ on the left and right faces. The electric charge has been redistributed, not created. More generally, the divergence theorem makes the cancellation explicit:Here $\mathbf{n}$ is the outward normal vector on the material boundary. This accounts for both bound volume charge density and bound surface charge density; omitting either can falsely suggest that an initially neutral dielectric acquired net electric charge.

An electrical insulator between capacitor plates supports a useful combination: appreciable electric polarization, little leakage electric current, and tolerance of the intended electric field. This lets a small parallel-plate capacitor store separated electric charge and recoverable capacitor energy. The comparison in dielectric capacitor response explains why specifying fixed electric charge or fixed electric potential difference changes the answer. High relative permittivity alone does not establish good insulation or a large usable energy density; dielectric loss and dielectric breakdown matter as well.

Dielectric strength characterizes the electric field at dielectric breakdown under specified test conditions. In an ideal uniform gap, it sets a voltage scale $V_{\mathrm{bd}}\simeq E_{\mathrm{bd}}d$.^[ref] The parallel-plate capacitor model consequently gives the stored energy density scaleThis is a model scale, not an operating rating: fringing field of a capacitor make local electric fields differ from $V/d$, while leakage and dielectric loss spend some supplied energy as heat. Comparing usable capacitor energy therefore requires insulation and breakdown performance as well as permittivity.

With signature $(-,+,+,+)$ and vanishing cosmological constant, the Einstein field equations areIn four spacetime dimensions, contraction with $g^{\mu \nu }$ gives the useful trace relation $R=-\kappa T$.

In the usual commutator convention for the Riemann curvature tensor, $\kappa =8\pi G/c^4$: $G$ is the gravitational constant, whereas $G_{\mu \nu }$ is the Einstein tensor. The Ricci tensor $R_{\mu \nu }$ and scalar curvature $R$ encode contractions of curvature. The stress-energy tensor $T_{\mu \nu }$ includes energy density, momentum flow and stresses, so pressure also gravitates. A cosmological constant adds $\Lambda g_{\mu \nu }$ to the left-hand side. See Tong's Einstein equations and matter coupling.

For a concrete consequence, set $T_{\mu \nu }=0$ and $\Lambda =0$. The trace gives $R=0$, and substitution gives $R_{\mu \nu }=0$. Yet the exterior Schwarzschild spacetime still has the tidal forces computed under geodesic deviation: vanishing Ricci tensor does not imply vanishing Riemann curvature tensor. Geometry carries gravitational effects through empty space. Einstein presented the final matter equations on 25 November 1915, including the essential trace correction, in Die Feldgleichungen der Gravitation, pages 844-847, especially equation (2a). His symbols and sign conventions differ from modern ones; the displayed formula uses modern notation.

For a homogeneous equilibrium system with energy density $\rho$, pressure $P$, chemical potential $\mu$, and number density $n$, the entropy density is

The entropy density of a photon gas is its entropy per physical volume, $s=S/V$. For homogeneous equilibrium vacuum radiation with two states of electromagnetic polarization and zero photon chemical potential,The first equality is the entropy density at zero chemical potential, obtained from the thermodynamic Euler relation $U=TS-pV+{\mu }_{\gamma }{N}_{\gamma }$ with internal energy $U=uV$. The second uses the photon gas equation of state, $u=a_R{T}^4$ and $p=u/3$. If an independently constrained photon number gives a nonzero chemical potential, the general entropy density is $(u+p-{\mu }_{\gamma }{n}_{\gamma })/T$ and $u=a_R{T}^4$ need not hold. This distinction explains why the absence of a conserved number matters for entropy as well as for the spectral shape. ^[ref]

A short first law of thermodynamics check, for neighboring equilibrium states with ${\mu }_{\gamma }=0$, isThe thermodynamic Euler relation fixes the additive normalization. Equivalently, this radiation entropy tends to zero with temperature at fixed volume; no contribution from the containing walls is included. Interpreting $TdS$ as transferred heat and $pdV$ as boundary work requires a reversible thermodynamic process. For a reversible thermodynamic process that is also an adiabatic process, there is no entropy production or entropy exchange, soAn irreversible adiabatic process can produce entropy even without heat exchange, so the word adiabatic alone is insufficient to infer $T^{3}V=\text{constant}$.

For an analytic example, expand an initially equilibrated photon gas reversibly from $V_0$ to $8V_0$, keeping only the same two photon states of electromagnetic polarization. Its temperature halves, its energy density becomes $u_0/16$, its entropy density of a photon gas becomes $s_0/8$, and its total entropy stays $S_0$. The internal energy halves rather than remaining fixed: the missing half leaves as work against radiation pressure,This is a stipulated reversible radiation model, not a report of a piston experiment. Maintaining a real cavity along this path requires controlling wall heat exchange and including the wall heat capacity in any actual apparatus balance.

For a computed scale at stipulated $T=2.725\mathrm{K}$, $s\simeq 2.04122\times 10^{-14}\mathrm{J}{\mathrm{K}}^{-1}{\mathrm{m}}^{-3}$. Combining this with the photon number density givesThis is the total equilibrium radiation entropy divided by its mean photon number, not a separately assigned entropy carried by each labeled particle. During the ideal expansion, $n_{\gamma }\propto T^3$ and $T^{3}V$ is fixed, so the mean number $N_{\gamma }=n_{\gamma }V$ also stays constant. This follows from the path; it was not imposed as the independent number constraint used to establish chemical potential.

The beautiful connection is that expansion lowers each typical photon energy and dilutes the radiation in just the proportions needed to retain its equilibrium spectral shape and total entropy. For the cosmic microwave background, homogeneous expansion of the universe has physical volume proportional to the cube of the scale factor; fixed-species cosmological entropy conservation therefore gives $T\propto a^{-1}$. This describes preservation of an already established Planck photon distribution and connects the FIRAS spectral observation in photon gas to an entropy prediction; FIRAS did not perform a direct entropy measurement. When other equilibrated species disappear and transfer energy into radiation, use the effective number of relativistic entropy degrees of freedom, $g_{\ast s}$: the conserved quantity for the combined sector is $g_{\ast s}{T}^{3}V$, under the corresponding equilibrium and no-entropy production assumptions. A changing species population, separate decoupled sectors or entropy production requires an appropriate entropy balance rather than the pure photon gas formula. ^[ref]

A finite-energy field configuration has integrable energy density after the vacuum energy has been subtracted. The boundary conditions used for a localized classical field-theory soliton require approaching a vacuum at spatial infinity. In one dimension the two ends may approach distinct scalar-field vacua; in a gauge theory the asymptotic field can have nontrivial winding even when its local energy density tends to zero.

For an internal gravity wave ray at angle $\theta$ to the horizontal reflecting from a slope at angle $\alpha$, conservation of tangential wavenumber givesAt a subcritical internal-wave reflection, $\gamma >1$: the reflected wavelength shortens and its wavelength-averaged energy density increases by ${\gamma }^2$.

The front velocity is the propagation speed of the earliest boundary of a newly introduced disturbance's support. It asks when an initially unaffected observer can first respond, whereas group velocity asks how a narrow wave envelope translates. A detector threshold defines a practical arrival time, but can miss arbitrarily small earlier responses; a packet peak is another distinct timing convention. In a causal dielectric model, a high-angular frequency contour argument gives exactly zero field for $t<x/c$, even when phase velocity exceeds $c$. This is a statement about that model's relativistic causality, not a universal identification of group velocity with front velocity.^[ref]

For the stipulated Klein-Gordon equation with compactly supported real Cauchy data, front velocity is bounded by $c$, and $c$ is the model's maximal speed of influence. A short energy method outline usesIntegration outside an expanding interval shows that zero initial energy there remains zero. Thus initial support in $[-a,a]$ cannot influence $x>a$ before $t=(x-a)/c$. The positive energy density includes the restoring term ${\omega }_0^2{u}^2/2$: this adds oscillation behind the front without extending its domain of dependence. A particular disturbance may fail to excite the leading boundary, so the maximal speed is not a promise of a detectable response for every source.

A Gaussian function is nonzero arbitrarily far away, so the Gaussian envelope in a Klein-Gordon medium has no sharp initial support boundary. Its small distant tails are already in the initial data. Likewise, keeping only a narrow interval of the frequency spectrum cannot reproduce an exact abrupt onset. Cutting off a Gaussian tail changes its Fourier transform and the accompanying initial velocity; it is a different initial-value problem. The exponential initial-velocity tail for the Klein-Gordon equation is another concrete example of a pre-existing tail outside the main disturbance.

The water-tank measurement in phase velocity records free surface elevation with calibrated resistance probes and uses a ruler and stopwatch for an independent check. The inspected protocol calls the front or back of a group a timing feature. That visible edge establishes an envelope transit estimate, not the first nonzero disturbance in the sense of front velocity.^[ref] A proposed additional control would repeat the transit fit at several detection thresholds and train durations, using matched probe calibration and the same sampling rate. Dependence on the threshold diagnoses wave envelope reshaping or limited sensitivity; records should be compared before reflections return. Such a test could expose limitations of a group velocity estimate, but would not establish an exact causal front or validate the stipulated Klein-Gordon equation for water waves.

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.

For a spatially homogeneous minimally coupled scalar field with zero scalar potential, the stress-energy tensor has $\rho =P={\dot{\varphi }}^2/2$. Its field equation is $\ddot{\varphi }+3H\dot{\varphi }=0$, so $a^3\dot{\varphi }$ is constant and its energy density decreases as $a^{-6}$. Thus its homogeneous stress is precisely that of a stiff fluid. Spatial gradients would introduce directional stresses, so this equivalence requires homogeneity.

For a fixed Newtonian gravitational potential, total energy density is $\mathcal{E}=\rho (u^2/2+\Phi )+p/(\gamma -1)+B^2/(2{\mu }_0)$. Its flux combines advected kinetic/potential energy, specific enthalpy and the Poynting vector. Dotting the ideal magnetohydrodynamic momentum equation with velocity, adding the adiabatic internal energy equation and the magnetic energy equation cancels the magnetic work. The pressure terms combine into $-\nabla \cdot (p\mathbf{u})$, proving the conservation law. A time-dependent imposed potential instead contributes $\rho {\partial }_t\Phi$.

A slowly established upward plane internal gravity wave has $\overline{u^{\prime }{w}^{\prime }}=-m|a{|}^2/(2k)$. The mean horizontal conservation of momentum equation gives ${\overline{u}}_t=-{\partial }_z\overline{u^{\prime }{w}^{\prime }}$. Using envelope transport and an initially vanishing mean gives $\overline{u}=-m|a{|}^2/(2k{c}_{gz})$. Its energy density per unit mass is $\mathcal{E}=(k^2+m^2)|a{|}^2/(2k^2)$, so $\overline{u}=(k/\omega )\mathcal{E}$. This is the mean flow left behind the advancing wave envelope, rather than a claim that an established uniform wavetrain has nonzero local stress divergence.