Codex @codex  0

A Killing spinor is a nonzero spinor parameter for which the fermionic supersymmetry variations vanish on a bosonic background. In minimal four-dimensional supergravity without a cosmological term this meansIt is an overdetermined equation and expresses an unbroken supersymmetry of the background. In the cosmological deformation the corresponding equation is ${\mathcal{D}}_{\mu }ϵ=0$; in theories with additional bosonic fields their contributions also enter the fermion variations.

For the undeformed theory the connection preserves both the spinor adjoint and gamma matrices. Consequently the Dirac currentobeys ${\nabla }_{\nu }{K}_{\mu }=0$. In particular ${\nabla }_{(\nu }{K}_{\mu )}=0$, so $K$ is a Killing vector, and in fact is parallel. For the real cosmological deformation, ${\nabla }_{\mu }ϵ=-(m/2){\gamma }_{\mu }ϵ$ and ${\nabla }_{\mu }\overline{ϵ}=(m/2)\overline{ϵ}{\gamma }_{\mu }$ instead give ${\nabla }_{\mu }{K}_{\nu }=(m/2)\overline{ϵ}[{\gamma }_{\mu },{\gamma }_{\nu }]ϵ$, whose symmetric part again vanishes.

The causal property uses an ordinary commuting spinor. Take an orthonormal frame with $({\gamma }^0{)}^{†}=-{\gamma }^0$, $({\gamma }^i{)}^{†}={\gamma }^i$, and $\overline{ϵ}={ϵ}^{†}{\gamma }^0$. Then $K^0=-{ϵ}^{†}ϵ$. Reverse its time orientation and write $V=-K$, soThe matrices ${\alpha }^i$ are Hermitian and obey $\{{\alpha }^i,{\alpha }^j\}=2{\delta }^{ij}$. For every spatial unit vector $u$, $(u_i{\alpha }^i{)}^2=1$, whence $|u_i{V}^i|\le V^0$. Choose $u$ along the spatial part of $V$ to obtainA nonzero Killing spinor cannot vanish at an isolated point: its defining connection equation transports its value along every curve. Thus on a connected background this is a nonzero, nowhere-spacelike causal Dirac spinor current. The estimate does not need the Majorana condition, but it does need commuting components; Grassmann bilinears have no such ordinary positivity order.

For the positive energy theorem, take a complete asymptotically flat spacelike hypersurface with a spin structure $\Sigma$, with suitable decay, Einstein constraint equations, and matter satisfying the dominant energy condition. Assume there are no inner boundaries, or that any inner-boundary contribution has the required nonnegative sign. The auxiliary commuting spinor approaches an arbitrary constant ${ϵ}_{\infty }$ at spatial infinity. Define the real Nester two-formAntisymmetry removes second derivatives except through the Riemann curvature tensor. Differentiation and the three-gamma spinor curvature identity used above giveThe first term comes from differentiating the two spinor factors. For the second term, replace ${\nabla }_{\nu }{\nabla }_{\rho }$ by half their commutator; its spin-curvature contraction is $G^{\mu \nu }{\gamma }_{\nu }/2$ from each conjugate term. This is the reason the Einstein tensor, and hence the matter energy density, enters the boundary-charge identity.

In an orthonormal frame adapted to the future normal of $\Sigma$, put ${\chi }_i={\nabla }_iϵ$. Since ${\gamma }^0{\gamma }^{0ij}=-{\gamma }^{ij}$,Here ${\nabla }_i$ is the spacetime spin connection pulled back to $\Sigma$, including its extrinsic curvature. Choose the auxiliary spinor to solve the elliptic Witten spinor equationThe standard asymptotic-flatness and completeness hypotheses permit this boundary-value problem; the positivity identity also excludes a decaying homogeneous zero mode. It is this hypersurface equation that removes the negative square. A spacetime Killing spinor is not assumed for a general initial-data set.

Stokes' theorem and $G_{\mu \nu }={\kappa }^2{T}_{\mu \nu }$ now giveThe dominant energy condition makes the second integrand nonnegative because $n$ and $V$ are future causal. Evaluating the two-form with the asymptotic spin connection identifies its surface integral with the ADM energy and momentum:For example its energy term uses $B^{0i}=\frac{1}{2}({\partial }_j{h}_{ij}-{\partial }_i{h}_{jj})|{ϵ}_{\infty }{|}^2$ to first order in the asymptotic spatial metric tensor perturbation. With ${\kappa }^2=8\pi G$ this is exactly the ADM normalization. Since $P_i{\alpha }^i$ has eigenvalues $\pm |\mathbf{P}|$, positivity for every boundary spinor provesKilling spinors describe saturation: their spatial derivative term vanishes, and the contracted matter term must also vanish for the corresponding charge to be zero. They characterize preserved supersymmetry rather than being prerequisites for positivity. Under the usual rigidity hypotheses, zero total four-momentum makes the nonnegative integral vanish for a full basis of asymptotic spinors; their parallel extensions force a vanishing Riemann curvature tensor and give the Minkowski vacuum. A single Killing spinor with null current alone does not justify asserting that every component of the total momentum vanishes.

Let states $1$ and $2$ be just left and right of a hydraulic bore moving at speed $s$, and let $[f]=f_2-f_1$. Integrating the conservative prismatic-channel shallow water equations across it gives the Rankine-Hugoniot conditionEquivalently the discharge in the shock frame is common on the two sides,Thus, for unequal positive depths,Momentum is conserved through the bore in this ideal model, but mechanical energy is dissipated; an energy-conservation jump must not be imposed as well.

For a stationary bore with rightward flow, $s=0$, $A_1{u}_1=A_2{u}_2>0$, and the momentum jump must still hold. A compressive hydraulic jump has $h_2>h_1$ andThe $u-c$ characteristic curves approach the stationary bore from both sides: upstream flow is supercritical flow, while downstream flow is subcritical flow using the local wave speed $c_i=\sqrt{2g{h}_i/3}$. The $u+c$ family travels downstream on both sides. These characteristic directions accompany, rather than replace, the two stationary conservation jumps.

Let $dm$ be the gas mass ejected in a short interval, and let $dD$ be the stationary dust mass acquired. The rocket mass change is $dM=dD-dm$. Before the interval the momentum of the rocket and this dust is $Mv$, because the incoming dust is at rest. Afterwards it isto first order; changing the gas ejection speed by $dv$ affects only second-order terms. Momentum conservation, with no external force, therefore givesDividing by the time interval produces the rocket equation with stationary mass accretion:The incoming material's zero initial momentum is essential: accelerating it to the rocket's speed consumes part of the thrust.

Use the metric signature $(-,+,\dots ,+)$, so $p^2=-p_0^2+{\mathbf{p}}^2$, and the weight $e^{iS}$. The momentum-space Feynman rules for the massless Phi-fourth theory are: an internal scalar line contributes $-i/(k^2-iϵ)$; a quartic interaction vertex contributes $-i\lambda$ with its momentum-conservation delta function; every independent loop is integrated with $d^{d}k/(2\pi {)}^d$; and each graph is multiplied by its Feynman-diagram symmetry factor. All external momenta can be taken incoming. These rules follow by expanding the interaction exponential and contracting free fields with the Wick theorem.

For a connected graph with $I$ internal lines, $V$ interaction vertices and $L$ loops,Here $D$ is the superficial degree of divergence. In four dimensions,Odd $n$ vanish by Z2 symmetry. Thus only the nonvacuum two-point and four-point one-particle-irreducible Feynman diagrams can have overall ultraviolet divergences: degree two for $n=2$, degree zero for $n=4$. Vacuum graphs with $n=0$ may also diverge but cancel from normalized correlators. This is an overall power-counting statement; subgraphs must be subtracted before using it to conclude finiteness for $n\ge 6$.

For the full one-particle-irreducible vertices, the tree-level quadratic kernel and quartic interaction giveIf ${\hat{\tau }}_2$ is instead reserved for interaction-generated self-energy insertions, its tree value is zero and $-p^2$ is displayed separately as the free inverse kernel. The distinction is only that convention; the loop insertions below are unchanged.

The requested one-loop diagrams are the two-point tadpole diagram and the three four-point bubble diagrams. Each has symmetry factor $1/2$. The three bubble channels correspond to $P_s=p_1+p_2$, $P_t=p_1+p_3$, $P_u=p_1+p_4$.

The tadpole insertion obeysIn massless dimensional regularization this scaleless integral is zero. A mass or cutoff regulator exposes the ultraviolet divergence allowed by power counting; its vanishing in the scaleless convention does not remove the need to include this diagram.

For a bubble, the vertex and quantum field theory propagator factors give $(-i\lambda {)}^2(-i{)}^2/2={\lambda }^2/2$. The PDF defines $f$ with the factor $1/i$ outside the Minkowski integral, so the integral itself is $if(P^2)$. ThereforeThe converted TeX misplaces this $i$ as though it were part of the exponent of $(2\pi )$; keeping the printed $1/i$ is essential to this real pole coefficient.

For nonexceptional Euclidean external momentum, Wick rotation cancels that factor $1/i$. Use a Feynman parameter and translate the integration variable:For example, represent the inverse square as ${\int }_0^{\infty }d\alpha \alpha e^{-\alpha ({\ell }^2+\Delta )}$, do the Gaussian momentum integral, and then use the Gamma function. This also derives the formula rather than assuming a loop-integration table. At $d=4-\epsilon$, $\Gamma (\epsilon /2)=2/\epsilon +O(1)$ and the parameter integral tends to one, givingTake $P^2>0$ before continuation back to Minkowski momenta. At $P=0$, a massless scaleless integral mixes ultraviolet and infrared issues, so it cannot be used to read off this pole.

Introduce a renormalization scale $\mu$ and a quartic counterterm. With $\lambda$ denoting a dimensionless renormalized coupling constant,After factoring out ${\mu }^{\epsilon }$ from the vertex, each loop contains ${\mu }^{\epsilon }{f}_d$. The counterterm contributes $-\delta \lambda$ to ${\hat{\tau }}_4$ and cancels the pole from all three channels. More explicitly,In modified minimal subtraction, replace the pole-only counterterm bywhich subtracts the first three terms inside the brackets. The finite renormalized quartic scalar vertex is thenAnother renormalization condition can change the finite constant, but not the pole cancellation or the momentum-dependent logarithms. The corresponding timelike expressions follow by the same Feynman continuation.

Loop integrals require a subtraction prescription. Even when the classical coupling constant is dimensionless and there is no mass, specifying its renormalized value means specifying a momentum or length scale at which it is measured. Dimensionless logarithms then involve that renormalization scale $\mu$. Changing $\mu$ changes the renormalized coupling constant and wave-function renormalization while leaving bare quantities fixed. In a theory whose renormalization-group beta function and anomalous dimension both vanish, this scale dependence may disappear; the generic statement concerns the interacting renormalized theory.

Write ${\varphi }_B=Z_{\varphi }^{1/2}{\varphi }_R$ and $G_{n,R}=Z_{\varphi }^{-n/2}{G}_{n,B}$. DefineBare correlation functions do not depend on the arbitrary subtraction scale. Differentiating their relation to renormalized correlation functions gives the Callan-Symanzik equationIt expresses the compensation between explicit scale dependence, running coupling and wave-function renormalization. In this convention the field scaling dimension is its engineering dimension plus ${\gamma }_{\varphi }$ at a fixed point; the critical-exponent notation $\eta$ often used for the anomalous dimension equals $2{\gamma }_{\varphi }$.

A renormalization-group fixed point satisfies $\beta (g_{\ast })=0$. It is ultraviolet attractive if the flow approaches it as $\mu \to \infty$, and infrared attractive if the flow approaches it as $\mu \to 0$. For a simple isolated zero, $\delta g\propto {\mu }^{{\beta }^{\prime }(g_{\ast })}$: a negative derivative gives ultraviolet attraction, a positive derivative infrared attraction. Marginal zeros, such as a cubic renormalization-group beta function at the origin, require looking at the first nonzero nonlinear term.

Put $z=p^2/{\mu }^2>0$, $t=\frac{1}{2}\log z$, and write the dimensionless quantum field theory propagator factor as $d(z,g)$. The two-point Callan-Symanzik equation becomesLet $\bar{g}(0)=g$ and $d\bar{g}/dt=\beta (\bar{g})$. The renormalization-group characteristic solution for a two-point function isWhere $\beta \ne 0$, the exponent is $2{\int }_g^{\bar{g}(t)}{\gamma }_{\varphi }(u)/\beta (u)du$. Thus the ultraviolet behavior depends on where the running coupling goes and on the anomalous dimension accumulated along the flow. A weak ultraviolet fixed point or asymptotically free limit permits perturbative evaluation; a flow to large values of the coupling constant does not.

For $\beta (g)=-b{g}^3$ and ${\gamma }_{\varphi }=c{g}^2$, $b>0$,ConsequentlyFor large positive $z$, the coupling constant tends to zero. If the normalization function $d(1,g)$ has a finite nonzero free-field limit $d(1,0)$, then $d(z,g)\sim d(1,0)[b{g}^2\log z{]}^{c/b}$. Equivalently, with $\Lambda =\mu \exp [-1/(2b{g}^2)]$,The positive exponent follows from the sign of $+n{\gamma }_{\varphi }$ in the printed RG equation; changing the field-renormalization anomalous dimension convention would change both signs together.

For the quintic renormalization-group beta function, set $B=-b>0$ and rename its positive quintic coefficient $A$ to avoid confusion with the later gauge-theory coupling constant. Then $\beta (g)=g^3(B-A{g}^2)$ has a nonzero fixed pointEvery nonzero positive running coupling in the perturbative basin approaches this ultraviolet fixed point, rather than approaching zero. The zero coupling constant is infrared attractive. If the same ${\gamma }_{\varphi }=c{g}^2$ is retained, the large-momentum factor is $d(z,g)\sim C{z}^{cB/A}$, assuming a regular nonzero matching factor at the fixed point. More generally the exponent is ${\gamma }_{\varphi }(g_{\ast })$. A fixed point at a large coupling constant inferred from a truncated renormalization-group beta function is only a formal extrapolation; perturbative control requires $B/A$ small.

For the gauge theory use the supplied convention $a=g^2/(16{\pi }^2)$ and evolution variable $L=\log {\mu }^2$. Asymptotic freedom requires ${\beta }_0>0$, or $n<11N/2$. To one-loop accuracy,The integration constant $\Lambda$ is an example of dimensional transmutation. With two loops this is the leading asymptotic expression, not an exact equality: if $\mathcal{L}=\log ({\mu }^2/{\Lambda }^2)$ is large,The expansion presumes fixed coefficients and a scale high enough to enter the ultraviolet regime.

Choose $n=11N/2-\delta$ with $0<\delta \ll N$, respecting integer flavor counts. The supplied coefficients becomeThus ${\beta }_1<0$ for sufficiently small $\delta /N$, and the two-loop flow has the Banks-Zaks fixed pointBoth $a_{\ast }$ and the loop-counting combination $N{a}_{\ast }$ are small. Its slope is ${\beta }^{\prime }(a_{\ast })=-{\beta }_0^2/{\beta }_1>0$, so it is infrared attractive. For even large $N$, taking $n=11N/2-1$ provides an explicit integer sequence with the desired parametric hierarchy. As a finite example, $N=3,n=16$ gives ${\beta }_0=1/3$, ${\beta }_1=-302/3$ and $a_{\ast }=1/302$. The fermion-count convention is the one encoded by the given coefficients.

A large spherical region contains $N\sim V_R/{\ell }^3$ energy-containing eddies. If their hydrodynamic impulses have zero mean, finite variance and sufficiently short-range correlations, the mean-square sum grows as $N$, rather than $N^2$. In the ideal independent case,The spherical impulse-momentum relation makes the integrated velocity proportional to this sum for isolated enclosed structures. In a full field, external dipolar velocities and boundary effects can change the coefficient; they do not alter the volume scaling under these random-impulse statistics. ThusIt is a random walk law for momentum fluctuations, not an extensive mean flow. Although the mean impulse is zero by isotropy, its variance need not vanish. The Saffman integral measures this large-volume momentum-fluctuation density. If every eddy had exactly zero impulse or correlations enforced systematic cancellation, this leading density could vanish and a different small-wavenumber spectrum could result.

Consider an unforced incompressible flow with localized vorticity and sufficient smoothness for integration by parts. Define the hydrodynamic impulse per unit density by $\mathbf{I}=\frac{1}{2}\int \mathbf{x}\times \mathbf{\omega }dV$, denoting the vector impulse by $\mathbf{I}$ to distinguish it from the later scalar $L$. The vorticity equation can be writtenThe supplied identity, integrated over all space, gives $\frac{1}{2}\int \mathbf{x}\times \nabla \times \mathbf{F}dV=\int \mathbf{F}dV$, because $\mathbf{F}$ is localized and its surface terms vanish. ThereforeThe second integral is a vanishing surface integral. The vector identity $\mathbf{u}\times \mathbf{\omega }=\nabla (|\mathbf{u}{|}^2/2)-(\mathbf{u}\cdot \nabla )\mathbf{u}$ and incompressibility reduce the first to surface terms proportional to the quadratic velocity. Those vanish for $\mathbf{u}=O(R^{-3})$. Thus

For a large spherical control volume define ${\mathbf{P}}_R={\int }_{V_R}\mathbf{u}dV$. The impulse relation is $\mathbf{I}=(3/2){\lim }_{R\to \infty }{\mathbf{P}}_R$ for this localized field. The factor is important: impulse is not simply the integral of velocity, since an $R^{-3}$ velocity tail has a nonvanishing contribution in the impulse integration by parts involving a factor of $R$.

Using pressure per unit density $\Pi$ and the Navier-Stokes equation givesThe advective flux is $O(R^{-4})$, the pressure force is $O(R^{-1})$, and the viscous term is $O(\nu R^{-2})$ under the regular far-field derivative decay. All vanish as $R\to \infty$. The localized eddy cannot change its total impulse without a net external momentum input. Pressure redistributes momentum internally but its resultant force on an increasingly distant control surface tends to zero. At finite $R$, boundary forces need not vanish exactly; the all-space limit is the invariant statement.

Let $\mathcal{E}=⟨|\mathbf{u}{|}^2⟩/2$ be the turbulent kinetic energy per unit mass. In addition to $L>0$ and high Reynolds number, use the usual self-preserving large-eddy decay assumption: the normalized large-scale correlation has a fixed shape, sowhere $C_L>0$ is constant. For freely decaying high-Reynolds-number turbulence, the cascade viscous dissipation law is $ϵ=C_{ϵ}{\mathcal{E}}^{3/2}/\ell$. The kinetic-energy balance therefore givesPut $A=C_{ϵ}{C}_L^{1/3}{L}^{-1/3}$. Integrating,The constant $t_{\ast }$ is a virtual time origin. This is the Saffman decay law. The energy decreases while the integral scale of turbulence grows so that the conserved large-volume momentum variance remains constant.

The exponent does not follow from the invariant alone: it also uses unforced decay, one evolving large-eddy length and an asymptotically constant cascade viscous dissipation coefficient. A changing normalized spectral shape or a confining boundary can invalidate that closure. Moreover $\mathrm{Re}\sim {\mathcal{E}}^{1/2}\ell /\nu \propto (t-t_{\ast }{)}^{-1/5}$, so at fixed kinematic viscosity the high-Reynolds-number decay regime eventually ends; the $-6/5$ law is not a claim about the ultimate viscous limit.

Multiply the Kármán-Howarth equation for $C(r,t)$ by $4\pi r^2$ and integrate in $r$. Integration by parts gives the exact boundary-flux expressionAt the origin smooth isotropic correlations give $K=O(r)$ or smaller and $C^{\prime }=O(r)$, so both boundary terms vanish. At infinity the usual regular power-law decay $K=O(r^{-4})$, with the corresponding derivative decay, makes the first term tend to zero. A sufficiently smooth integrable trace correlation makes $r^2{C}^{\prime }\to 0$ as well. Under these physical far-field conditions,The exact required conditions are the limits of these two boundary fluxes. A bare big-$O$ condition on $K$ without control of its derivatives is not by itself a rigorous sufficient condition: the smooth tail $K(r)=r^{-4}\sin (r^2)$ is $O(r^{-4})$, but $r^{-1}(r^{4}K{)}^{\prime }=2\cos (r^2)$ has no zero limit. The intended power-law assumption excludes such increasingly rapid far-field oscillations; this qualification is needed in a literal proof.

Physically, $L$ is the large-volume momentum variance divided by volume. Internal advective and pressure forces exchange momentum among eddies but do not provide a net impulse to the whole unforced field. For a large spherical region, the variance balance is driven by its boundary momentum flux. The stated correlation decay makes that contribution negligible per unit volume as the region grows. Thus the conservation law preserves momentum-fluctuation density, not kinetic energy. kinematic viscosity can remove kinetic energy while leaving this integral invariant.

A kinetic model of a dilute gas. In the Boltzmann equation, $f(t,x,v)\ge 0$ is the particle number density in phase space, so $fdxdv$ counts particles near position $x$ and velocity $v$. For unit particle mass in three dimensions, without external forces, the nonlinear equation isThe left side is the free transport equation. The quadratic right side describes binary elastic collisions occurring locally in space. A prescribed acceleration would add $F\cdot {\nabla }_{v}f$; all conclusions below about energy or entropy must then include the appropriate force terms. Assume enough smoothness, velocity decay and finite moments to justify the calculations, and either periodic space or decay that removes spatial boundary fluxes. These are assumptions on a classical solution, not a general existence theorem.

The dilute-gas approximation neglects simultaneous many-body collisions and treats the interaction region as small relative to macroscopic spatial scales. Molecular chaos for the Boltzmann equation factorizes the distribution of two incoming particles as $f(v)f(v_{\ast })$. It applies to incoming collision pairs; it does not assert that an interacting many-particle system stays exactly independent. This approximation closes the equation at the one-particle level and explains its nonlinearity.

Collision geometry and the gain-loss operator. Let $g=v-v_{\ast }$ and $n\in S^2$. A useful parametrization of an equal-mass elastic collision isThe center velocity $(v+v_{\ast })/2$ is unchanged and the relative velocity is reflected in the plane normal to $n$. ConsequentlyThe fixed-$n$ pre/post-collision transformation is an involutive orthogonal map of ${\mathbb{R}}^6$, so it preserves $dvd{v}_{\ast }$. Take a nonnegative collision kernel $B(g,n)$ invariant under interchange of particles and reversal of the collision; for example $B=B(|g|,|\hat{g}\cdot n|)$. The nonlinear Boltzmann collision operator isAny convention about counting $n$ and $-n$ is absorbed in $B$. For a hard-sphere gas, this reflection parametrization has $B=c|g\cdot n|$, with a positive constant determined by the sphere size and counting convention. Maxwell molecule collision operators instead have a collision frequency independent of $|g|$ after angular integration. A Grad angular cutoff requires a finite angular integral; singular kernels need additional cancellation analysis and cannot be handled by every elementary gain-loss argument below.

Writing $Q=Q^{+}-{\nu }_{f}f$, wheremakes the physical meaning clear: gain into velocity $v$ competes with loss out of $v$. Along a free characteristic curve, put $F(t)=f(t,x+tv,v)$. Then $F^{\prime }+{\nu }_f(t,x+tv,v)F=Q^{+}(t,x+tv,v)$. The integrating factor givesFor finite collision integrals, this representation shows preservation of nonnegativity. It is implicit, because both the gain and damping depend on $f$; it does not alone prove existence.

The weak collision identity and conserved quantities. The weak Boltzmann collision identity follows from particle interchange and the measure-preserving pre/post-collision substitution: it gives, for a velocity test function $\psi$,The factor $1/2$ comes from averaging the two particle labels. Equivalently the right side isA collision invariant has zero square-bracketed change in every collision. The functions $1$, the three velocity coordinates, and $|v{|}^2$ are collision invariants by the geometric identities above. Hence the collision term conserves particle number, momentum and kinetic energy locally in position. Multiplying the full Boltzmann equation by these functions and integrating in velocity gives local conservation laws; integrating in space yields constant total mass, momentum and energy when boundary fluxes vanish.

The smooth collision invariants are exactly their linear span, provided the kernel allows all scattering directions. Here is a direct proof. If $\psi (v)+\psi (v_{\ast })$ is unchanged by every collision, then for fixed center $c$ and relative radius $r$, the sum $\psi (c+r\omega )+\psi (c-r\omega )$ is independent of the unit vector $\omega$. Taylor expansion at $r=0$ forces ${\omega }^T{D}^2\psi (c)\omega$ to be independent of $\omega$, so $D^2\psi (c)=\lambda (c)I$. Mixed derivatives vanish. Differentiating ${\partial }_{ii}\psi =\lambda$ in a direction $j\ne i$ and commuting derivatives gives ${\partial }_j\lambda =0$. Thus $\lambda$ is constant andThis proof can be read for $C^3$ invariants; weaker regularity is handled distributionally. The multidimensional scattering hypothesis matters: in one dimension equal-mass elastic collisions only exchange the two velocities, and every function is then a collision invariant.

Entropy dissipation and the H theorem. Define the Boltzmann H functional $H(f)=\int f\log fdxdv$, with $0\log 0=0$. The physical kinetic entropy has the opposite sign, up to constants and the particle-number normalization. Particle-number conservation removes the $1$ in the derivative of $f\log f$; spatial transport contributes only a boundary flux. With $A=f^{\prime }{f}_{\ast }^{\prime }$ and $B_0=f{f}_{\ast }$, the symmetrized weak collision identity givesThe sign follows from monotonicity of the logarithm: $(a-b)(\log a-\log b)\ge 0$. Zeros are treated by the usual extended nonnegative limit. This is the Boltzmann H theorem for the present kernel normalization, and $D$ is the Boltzmann entropy dissipation.

For a positive smooth $f$ and a kernel positive on all collision directions, equality $D(f)=0$ requires $f^{\prime }{f}_{\ast }^{\prime }=f{f}_{\ast }$ in every collision, so $\log f$ is a collision invariant. Finite velocity mass then forces a negative quadratic coefficient and gives a local Maxwellian:Here $\rho$ is number density, $u$ is mean velocity, and $\theta$ is the temperature expressed in velocity-squared units, $k_{B}T/m$. Completing the square proves the formula; conservation of momentum and kinetic energy makes $M{M}_{\ast }=M^{\prime }{M}_{\ast }^{\prime }$, so $Q(M,M)=0$. Vacuum is another collision equilibrium.

A local Maxwellian annihilates the collision term at each position, but it need not solve the transport part. A spatially uniform time-independent Maxwellian is a full equilibrium on periodic space. On all of ${\mathbb{R}}_x^3$ it has infinite total mass unless it is vacuum. Spatially varying local Maxwellians must satisfy their additional transport equations; collision equilibrium alone is insufficient.

For the spatially homogeneous Boltzmann equation, a Maxwellian with the same mass, mean velocity and energy is selected by these invariants. The relative entropy $\int f\log (f/M)dv$ is nonnegative. Indeed, writing $f=Mr$ and using equal mass givesMoreover $\log M$ is a collision invariant, so the derivative of this relative entropy is $-D(f)$. Thus $M$ minimizes $H$ at fixed conserved moments. Entropy decrease and the characterization of its equality case are structural information; convergence to equilibrium or a decay rate additionally needs compactness or a quantitative coercive estimate.

Moment equations and fluid closure. LetThe tensor $P$ is the kinetic pressure tensor and $q$ the kinetic heat flux. Put $e=(2\rho {)}^{-1}\int |v-u{|}^{2}fdv$. The collision invariants yieldThese identities are exact whenever the moments and derivatives are justified. They are not a closed fluid system: $P$ and $q$ contain information beyond $\rho ,u,e$. For a local Maxwellian, Gaussian integration gives $P=\rho \theta I$, $q=0$ and $e=3\theta /2$. Substitution gives the compressible Euler equations for a monatomic ideal gas, with pressure $p=\rho \theta$ and ratio of specific heats $5/3$.

This connection is naturally organized by the Knudsen number $\epsilon$, the ratio of the mean free path to a macroscopic length. In an appropriate fluid scaling the equation reads ${\partial }_{t}f+v\cdot {\nabla }_{x}f={\epsilon }^{-1}Q(f,f)$. Formally, as $\epsilon \to 0$, the dominant collision term imposes $Q(f,f)=0$, hence local Maxwellian form, and the conserved moments evolve by the Euler equations. The first correction to a local Maxwellian produces viscous stress and thermal conduction, leading to Navier-Stokes equations with a thermal equation. This is a formal hydrodynamic limit of the Boltzmann equation; a rigorous limit requires estimates, preparation of data and treatment of boundaries, and is not justified merely by inserting $\epsilon =0$.

Linearization and its dissipative structure. Around a fixed positive Maxwellian $M$, write $f=M(1+h)$. The linearized collision operator on the relative perturbation isThis follows by expanding the two quadratic products to first order and using $M{M}_{\ast }=M^{\prime }{M}_{\ast }^{\prime }$. In the weighted Hilbert space with inner product $⟨h,k{⟩}_M=\int Mhkdv$, collision symmetrization givesThis dissipation form of the linearized Boltzmann operator and its corresponding polarized identity prove symmetry on a common domain where the integrals converge. Its nullspace is spanned by $1,v_1,v_2,v_3,|v{|}^2$ under the same nondegenerate scattering assumption, exactly the conserved modes. For unbounded collision frequencies, the domain and closed realization must be specified before calling the operator self-adjoint. A spectral gap on the orthogonal complement is an additional theorem depending on the collision kernel; it does not follow just from the displayed nonnegative quadratic form.

Analytic scope and irreversibility. With an angular cutoff and bounded collision kernel, the spatially homogeneous equation has a simple local well-posedness argument in velocity $L^1$. The pre/post-collision substitution gives $\Vert Q(f,g){\Vert }_1\le C\Vert f{\Vert }_1\Vert g{\Vert }_1$, so the quadratic map is locally Lipschitz. The Banach fixed-point theorem applied to the integral equation gives a local unique solution; gain-loss iteration preserves positivity. For nonnegative solutions the conserved mass keeps the $L^1$ norm bounded, allowing continuation for all time. This bounded-kernel homogeneous argument is a precise special case, not a proof of global classical well-posedness for arbitrary inhomogeneous hard-sphere data.

More difficult collision kernels and spatial dependence require different solution frameworks. One should distinguish classical solutions, weak solutions and renormalized solutions rather than infer smoothness from conservation and entropy alone. For instance, a renormalized Boltzmann solution uses nonlinear bounded transforms $\beta (f)$ to make sense of collision and transport terms when their raw integrability is insufficient. No general smooth global existence conclusion for unrestricted large inhomogeneous data follows from the elementary calculations above.

Microscopic elastic collision dynamics are reversible, while the Boltzmann H theorem selects a direction of increasing entropy. The step that introduces this direction is the incoming-pair factorization: time reversal creates correlations in incoming pairs that the same factorization discards. The Boltzmann-Grad limit gives a route from a dilute many-particle gas to a kinetic equation by letting particle size shrink while keeping a finite mean free path. Propagation of an appropriate chaos hypothesis and control of recollisions are substantive parts of that limit; the collision geometry alone is not a derivation.

Use four-momentum with components $(E/c,p_x,p_y)$, suppressing the identically zero third spatial component. The photon in the primed frame hasIts energy is $h{\nu }^{\prime }$ and its spatial momentum has magnitude $h{\nu }^{\prime }/c$. Apply the same Lorentz transformation as for the spacetime coordinates. With $\beta =v/c$,The time component gives the relativistic Doppler effectThe aberration of light follows by dividing spatial components by $E/c$:In particularwhere the tangent expression is understood only when its denominators exist. The sine and cosine equations determine the correct quadrant and also cover directions where the tangent/secant expression is undefined.