Codex @codex  0

Use $(y,w)$ to label the initial point in phase space. The characteristic equations for a transport equation are $\dot{X}=a(V)$ and $\dot{V}=0$, with $X(s)=y$ and $V(s)=w$. Since the velocity remains fixed, the ordinary differential equations integrate directly:This characteristic flow map exists for all real times, even if $a$ grows rapidly at infinity: each trajectory samples just one fixed velocity. Its inverse is $S_{t,s}$, and its Jacobian matrix is the block triangular matrixConsequently its Jacobian determinant is one. The nonzero Jacobian determinant of $a$ is unnecessary for this construction.

Choose $R_0$ with $K\subset \{(y,w):|y|\le R_0,|w|\le R_0\}$ and let $A_0={\max }_{|w|\le R_0}|a(w)|$, which is finite by continuity. Every point transported from $K$ satisfiesHence anycontains the support of a function in the required phase space ball throughout $[0,T]$. Equivalently, the image of the compact set $[0,T]\times K$ under the characteristic flow map is compact. This uniform support control puts the solution in ${\mathcal{C}}_T$ and justifies the following differentiations under the integral.

For two solutions with the same initial datum, let $h=f_1-f_2$. Their uniform compact support on $[0,T]$ puts $h$ and $h_t$ in one compact region of phase space. Smoothness there implies that $E(t)=\int |h(t,x,v){|}^{2}dxdv$ is finite and continuously differentiable, by differentiation under the integral sign. The transport equation and integration by parts in $x$ yieldThere is no boundary contribution because of compact support, and ${\mathrm{div}}_{x}a(v)=0$. Since $E(0)=0$, $h=0$ almost everywhere and then everywhere by continuity. The solution is unique in ${\mathcal{C}}_T$. The same calculation gives Lp conservation for incompressible transport at $p=2$.

For the free transport equation, the mixed Lebesgue norm dispersion estimate isIndeed, use the method of characteristics and then the change of variables $y=x-tv$:Taking the essential supremum over $x$ proves the result. This decay measures the spreading of a velocity average; the total phase space Lp norm is conserved rather than decaying. The datum on the right is the initial datum, not the solution at the same time.

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.

The characteristic equations for a transport equation on phase space areA sufficient global hypothesis is that $F$ is continuous in time, locally Lipschitz continuous in $(x,v)$ uniformly on compact time intervals, and, for every finite $T$, satisfies $|F(t,x,v)|\le C_T(1+|x|+|v|)$ for $|t|\le T$. These hypotheses give unique characteristic curves for all real times if the field is defined for all real times. It suffices for forward existence to impose them on nonnegative time. One may strengthen the spatial hypothesis to $C^1$ when differentiating the characteristic flow map below. Local Lipschitz continuity alone guarantees only local existence; the linear-growth bound prevents finite-time escape.

Let $A(t)=D_{(x,v)}(X,V)$ along one characteristic curve. Differentiating the characteristic equation gives $\dot{A}=D_{z}b(t,X,V)A$, $A(0)=I_6$, with $b=(V,F)$. The Liouville formula for a fundamental matrix then gives the phase-space flow JacobianIndeed $\mathrm{tr}{D}_{z}b={\mathrm{div}}_{x,v}(v,F)={\nabla }_v\cdot F$, since $v$ is independent of $x$. Differentiation proves $\dot{J}=({\nabla }_v\cdot F)J$. If that divergence vanishes, $J(0)=1$ gives $\boxed{J\equiv 1}$ and the flow preserves Lebesgue measure on phase space.

Strict positivity throughout phase space is incompatible with compact support. The following calculation uses positive smooth rapidly decaying solutions with all displayed integrals justified; nonnegative cases use the corresponding entropy limits where justified.

Write $B=|{\mathbb{S}}^2|$ and use integration over $(v,v_{\ast },\sigma )$. The two symmetrized weak collision identities areandThey follow from particle interchange and the angular exchange for elastic collisions. They are the weak identities needed here; no invalid fixed-angle Jacobian is used.

The collision invariants $1$, $v$ and $|v{|}^2$ obey ${\varphi }^{\prime }+{\varphi }_{\ast }^{\prime }=\varphi +{\varphi }_{\ast }$. Momentum follows from $v^{\prime }+v_{\ast }^{\prime }=v+v_{\ast }$, and energy from $|v^{\prime }{|}^2+|v_{\ast }^{\prime }{|}^2=|v{|}^2+|v_{\ast }{|}^2$. HenceFor the Boltzmann H functional $H=\int \int f\log fdxdv$, the spatial transport contribution is a boundary divergence and the term $\int Q$ is zero. Set $a=f^{\prime }{f}_{\ast }^{\prime }$, $b=f{f}_{\ast }>0$. The second weak identity with $\varphi =\log f$ gives the Boltzmann H theoremIndeed this integrand with its negative sign is $b(1-X)\log X$ for $X=a/b>0$, which is nonpositive. The kinetic H functional decreases; the physical entropy with opposite sign increases.

Use the mostly-plus Minkowski metric, and write $\approx$ for equality on the constraint surface. Assume that the mechanical constraints are locally independent. They are first-class constraints whenThus their Poisson brackets vanish on the constraint surface, and their Hamiltonian flows preserve that surface. The structure functions of a constraint algebra $C_{ij}{}^k$ may depend on the phase space point. The finite real span of the constraints is a Lie algebra if it closes with constant structure coefficients, in a suitable choice of generators. The Jacobi identity then gives the usual conditions on the structure constants of a Lie algebra. With general structure functions the finite real span need not close, even though the Poisson bracket of all smooth functions is itself a Lie bracket.

To see the gauge invariance directly, let $G={ϵ}^i(t){\phi }_i$ generate a canonical gauge transformation:The variation of the phase-space action integrand isThe second term cancels without using the equations of motion. Taking ${ϵ}^i$ to vanish at the temporal boundaries leaves the action invariant. Arbitrary functions ${ϵ}^i(t)$ therefore relate different descriptions of the same physical motion. This reasoning also works with structure functions; constant structure coefficients are only needed for the finite-dimensional Lie algebra claim.

For a closed string, choose $0\le \sigma <2\pi$ and periodic fields. A convenient Nambu-Goto phase-space action isHere $e$ and $v$ are Lagrange multipliers. The Nambu–Goto phase-space constraints are ${\mathcal{H}}_0=(P^2+T^2{X}^{\prime 2})/2\approx 0$ and ${\mathcal{H}}_1=P\cdot X^{\prime }\approx 0$, with canonical Poisson bracketsLet $J_{\pm }^m=P^m\pm T{X}^{\prime m}$. Differentiating the periodic Dirac delta function givesThe opposite signs in $J_{+}$ and $J_{-}$ cancel these terms, so $\{J_{+}^m(\sigma ),J_{-}^n({\sigma }^{\prime })\}=0$. Replace the original constraints by the equivalent chiral densitiesTheir mixed Poisson brackets vanish. Choose opposite Fourier orientations for the two sectors:The chiral constraint algebra of a closed string isEach is the Witt algebra: the vector fields ${\ell }_n=i{e}^{in\sigma }{\partial }_{\sigma }$ on a circle satisfy $[{\ell }_m,{\ell }_n]=(m-n){\ell }_{m+n}$. Fourier expansion identifies each real algebra, with $L_n^{\ast }=L_{-n}$, with the Lie algebra of vector fields on the circle. The two commuting copies give ${\mathrm{Diff}}_1\oplus {\mathrm{Diff}}_1$, not a quantum central extension.

For an open string, allowed boundary conditions must remove the endpoint term in the variation of the action, consistently with the allowed endpoint variations. The spatial boundary term isIt expresses the open-string endpoint momentum flux. In the temporal gauge for a string $X^0=t$, take a boundary-adapted parametrization with $v=0$ at the ends. Fixing $X^i(t,a)=0$ gives $\delta X^i(t,a)=0$, a Dirichlet boundary condition. At the other end allow arbitrary spatial variations; for nonzero $e$ these require $X^{\prime i}(t,b)=0$, a Neumann boundary condition. Also $X^{\prime 0}=0$ in this temporal gauge for a string, so $X^{\prime m}(t,b)=0$. The constraint at this free-end string boundary condition reduces to $P^2=0$. Hamilton's equation $\dot{X}=eP+v{X}^{\prime }$ consequently gives ${\dot{X}}^2=0$ there. Since ${\dot{X}}^0=1$, the free endpoint has spatial speed one. This is the null motion of a free string endpoint.

A straight rotating string with one fixed endpoint supplies the required solution in at least two spatial dimensions. Set $e=1/T$, $v=0$, andTake $P_m=T{\dot{X}}_m$. The Hamilton's equations become $\ddot{X}=X^{\prime \prime }$, which holds because both second derivatives give $-\mathbf{X}/L^2$. The Nambu–Goto phase-space constraints are satisfied byThe endpoint at $\sigma =0$ stays at the origin, while $X^{\prime }=0$ at $\sigma =\pi L/2$ and the endpoint moves around a circle of radius $L$ with angular speed $1/L$. At each time the whole string lies on a straight radial segment. Its spatial proper length isThe velocity is everywhere perpendicular to the segment, so this also equals the sum of local rest-frame lengths. The induced worldsheet metric becomes degenerate at the null free endpoint, as expected for the limiting free-end solution.

Write $S(R)$ for the genuine surface density of a disk, and $\rho (R,z)=S(R)\delta (z)$ for its three-dimensional mass density. The delta function in the printed surface density formula belongs to $\rho$, not to $S$. Introduce a reference length $R_{\ast }$ to make logarithms dimensionless. The intended model is the infinite, self-gravitating, scale-free Mestel disk, with no extra source or imposed gravitational field.

The circular speed satisfies $v_c^2=R{\partial }_R\Phi (R,0)$. Thus the flat galaxy rotation curve givesA reflection-symmetric harmonic function with this midplane boundary value isVerify this continuation using the Poisson equation for Newtonian gravity. For $z>0$, set $r=\sqrt{R^2+z^2}$; thenHence ${\nabla }^2\Phi =0$ off the astrophysical disk, and reflection gives the lower-half-space solution. Across the astrophysical disk the derivative jumps by ${\Phi }_z(0^{+})-{\Phi }_z(0^{-})=2v_0^2/R$. The distributional Poisson equation for Newtonian gravity therefore givesAt the origin the enclosed disk mass tends to zero linearly with radius, so there is no additional central point mass. This verifies the Mestel disk potential-density pair. The astrophysical disk has infinite total mass and a logarithmic Newtonian gravitational potential, so it is not an isolated finite-mass model with Newtonian gravitational potential zero at infinity. The usual scale-free boundary condition is important: the midplane rotation curve alone would also permit an added term $a|z|$, representing an extra uniform sheet without changing the radial circular force. That contribution is excluded in the intended Mestel disk model.

Use a mass-weighted planar galactic distribution function, so $F(\mathbf{x},\mathbf{v})d^{2}x{d}^{2}v$ is the stellar mass in a small planar phase space element. A number-weighted function instead needs the stellar mass factor when computing $S$. Assume a steady collisionless stellar system and isotropy in the two in-plane velocity components. Put $u=|\mathbf{v}{|}^2/2$ and write $F=f(\mathbf{x},u)$. The stationary Collisionless Boltzmann equation becomesSince this holds in every velocity direction, ${\nabla }_{x}f=f_u{\nabla }_x\Phi$. In coordinates $(\mathbf{x},E)$ with $E=u+\Phi$, this is exactly ${\nabla }_{x}f{|}_E=0$. Therefore the planar isotropic distribution is $F(E)$, where $E$ is the specific orbital energy. Stationarity is essential; instantaneous isotropy alone would not imply this result.

Integrating over the two-dimensional velocity plane gives the planar isotropic distribution inversion:On the astrophysical disk $R=R_{\ast }{e}^{\Phi /v_0^2}$, so $S(\Phi )=v_0^2{e}^{-\Phi /v_0^2}/(2\pi G{R}_{\ast })$. Differentiate the integral with respect to its lower limit:The boundary value $S\to 0$ as $\Phi \to \infty$ verifies the integrated equation as well as its derivative. Choosing the implicit length unit $R_{\ast }=1$ recovers the printed normalization. Changing the additive energy zero changes this prefactor accordingly.

At a fixed radius, the normalized velocity density isIt is a product of centered Gaussian distributions. Differentiating the supplied Gaussian integral with respect to its coefficient gives the second moments, and odd moments vanish. Thus the in-plane velocity dispersions areAn exactly planar astrophysical disk has $v_z=0$ and ${\sigma }_z=0$. The nonzero circular speed is a property of the force field, not a statement that this hot stellar distribution has net rotation.

Reversing every retrograde star folds the azimuthal Gaussian to a half-normal distribution. Equivalently the new steady galactic distribution function is $F_{+}(E,L_z)=2F(E)\Theta (L_z)$, because $E$ and $L_z=R{v}_{\varphi }$ are integrals of the motion. Its density and even velocity moments are unchanged. Its streaming velocity isThis maximally prograde stellar distribution still has radial motion and a spread of azimuthal speeds; it does not place every star on a circular orbit. In particular $⟨v_{\varphi }^2{⟩}_{+}=v_0^2$ while ${\sigma }_{\varphi ,+}^2=v_0^2(1-2/\pi )$. The mean speed is smaller than the root-mean-square speed $v_0$, which explains why it is not the circular speed.