Codex @codex  0

A hard-sphere gas is a classical model of identical impenetrable spheres, each of diameter $\sigma$ and radius $\sigma /2$. Its pair potential is $u(r)=+\infty$ for centre separation $r<\sigma$ and $u(r)=0$ for $r\ge \sigma$; the convention at exact contact has no effect on the integral. It isolates the effect of particle size without an attractive potential energy.

The upper-left panel is a cross-section of a three-dimensional forbidden region: the second centre cannot enter a sphere of radius $\sigma$ around the first centre. The particles themselves have radius $\sigma /2$. Consequently, with particle volume $v_{\mathrm{s}}$ and pair excluded volume $v_{\mathrm{ex}}$,The last equality follows directly from the Mayer f-function: $f=-1$ throughout the forbidden volume and zero outside it. Confusing particle radius with contact separation changes the answer by a factor of eight; omitting the unordered-pair factor changes it by a further factor of two.

Define the packing fraction $\eta =n{v}_{\mathrm{s}}$, where $n$ is number density. Then the virial expansion beginsFor example, at the toy value $\eta =0.01$, the leading correction predicts $Z\approx 1.04$; the omitted terms still matter if accuracy beyond this approximation is needed. The classical second virial coefficient is independent of temperature because the geometry of forbidden configurations does not change when the gas is warmed. The pressure itself still scales with temperature at fixed number density.

This model makes the positive correction particularly clear: eliminating overlapping configurations changes the dependence of the available configuration space on volume, even though the allowed configurations all have zero interaction energy. The figure's other panels add a toy attractive shell and show how its Boltzmann factor competes with that geometry. The hard-sphere gas is a useful reference model, while a low-order virial expansion remains a dilute-gas approximation rather than a dense sphere packing formula.

The Knudsen number is the dimensionless ratio of an intermolecular mean free path ${\ell }_{\mathrm{mfp}}$ to a specified characteristic length $L$ of the apparatus or flow. The mean free path measures the distance travelled between gas–gas collisions; $L$ measures the distance over which boundaries or macroscopic conditions matter. Wall encounters must not be counted as gas–gas collisions in this comparison: doing so would hide the very wall-dominated behavior we want to detect. The same gas can therefore have a small Knudsen number in a large vessel and a large Knudsen number in a tiny opening.

For a dilute hard-sphere gas of identical molecular diameter $d$ in thermal equilibrium, use the conventionHere $n$ is the number density, $p$ the pressure, $T$ the absolute temperature, and $k_B$ the Boltzmann constant; the second equality uses the ideal gas law. The collision cross-section is $\pi d^2$, and the $\sqrt{2}$ accounts for both partners moving with independent equilibrium velocities. Mixtures and velocity-dependent collision cross-sections require a stated molecular model and mean free path convention; the displayed expression is not an exact formula for every real gas. ^{[ref] [ref]}

As a numerical toy example, choose $T=300\mathrm{K}$ and $d=0.37\mathrm{nm}$, without claiming that this is a measured molecular diameter for a particular gas. The formula gives ${\ell }_{\mathrm{mfp}}\simeq 68.1\mathrm{nm}$ at $p=10^5\mathrm{Pa}$. It yields $\mathrm{Kn}\simeq 6.81\times 10^{-6}$ for a $1\mathrm{cm}$ vessel but $\mathrm{Kn}\simeq 0.681$ for a $100\mathrm{nm}$ gap. At $p=1\mathrm{Pa}$ the same toy gas has ${\ell }_{\mathrm{mfp}}\simeq 6.81\mathrm{mm}$, comparable to the centimetre vessel. Reducing either pressure or apparatus size can thus expose molecular transport.

Choose $L$ for the question being asked: a tube diameter compares the mean free path with wall separation, while tube length describes a different geometric scale. For a positive macroscopic field $q$, a local gradient scale can instead be $L_q=|q|/|\nabla q|$, where its gradient is nonzero. A large vessel can contain a thin region with a large local Knudsen number. Small geometric Knudsen number alone does not guarantee a slowly varying near-equilibrium flow. The time, velocity and space scalings must also be stated when the Knudsen number appears as a coefficient in a Boltzmann equation. ^[ref]

The Mayer f-function measures how a pair potential changes a pair's Boltzmann factor relative to the noninteracting value one:It vanishes where $u=0$, lies between $-1$ and $0$ for repulsion, and is positive for attraction. Even an infinite hard-core potential energy produces the finite value $f=-1$. This is why expanding in Mayer f-functions handles a hard-sphere gas more naturally than a Taylor series in $u$. ^[ref]

For pairwise-additive potential energy $U={\sum }_{i<j}u(r_{ij})$, the exact identityturns the configurational part of the canonical partition function into a sum of products associated with pairs. Drawing one edge for each chosen factor $f(r_{ij})$ gives the graphs used in the cluster expansion. The second virial coefficient uses just one pair's relative-position integral, $B_2=-\frac{1}{2}\int f{d}^3\mathbf{r}$.

For the square-well toy pair potential in second virial coefficient, $f=-1$ in the core, $f=e^{1/T^{\ast }}-1$ in the attractive shell, and $f=0$ beyond it. Its positive plateau need not be small: low number density, rather than small $f$ at every separation, controls the dilute-gas approximation. The integral weights each plateau by spherical-shell volume; equal-looking unweighted areas on a plot of $f$ versus $r$ need not cancel in $B_2$.

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 Clausius gas equation can be writtenComparison with the virial expansion givesFor a hard-sphere gas of atomic radius $r_0$, the pair potential is infinite when the center separation is below $2r_0$ and zero otherwise. ThusThis is four times the volume of one atom and is the per-particle excluded volume appearing in the dilute equation of state.

Molecular transport across gradients. A gas can transport momentum, thermal energy and species abundance even when its average bulk velocity is zero. Molecules arriving at a point carry properties acquired near their preceding collisions. A small gradient makes the populations arriving from opposite sides slightly different, producing transport from larger to smaller values of the carried property. A continuum transport law requires a mean free path $\lambda$ small compared with the gradient length, and relaxation times short compared with the observation time.

For a dilute hard-sphere gas,where $d$ is molecular collision diameter. Both speed and collision distance matter: high speed alone does not determine a transport coefficient.

A simple kinetic estimate makes the common structure explicit. For a quantity $a(y)$ carried per particle, a molecule moving at angle $\theta$ to $+y$ brings approximately $a(y-\lambda \cos \theta )$. Expanding to first order and averaging an isotropic velocity distribution givessince $⟨\cos \theta ⟩=0$ and $⟨{\cos }^2\theta ⟩=1/3$. This mean-free-path argument estimates the coefficients; a detailed collision calculation changes their numerical prefactors and distinguishes different relaxation rates.

Momentum transport and viscosity. For a shear flow with mean velocity $U_x(y)$, the transported directed momentum per molecule is $m{U}_x$. The net flux of $x$ momentum in the $+y$ direction isHere $\eta$ is dynamic viscosity in $\mathrm{Pa}\mathrm{s}$. The associated shear stress on a material surface has the corresponding traction sign convention; the flux carries momentum from the faster-moving layer towards the slower-moving one. The kinematic viscosity is $\nu =\eta /(nm)$, with units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Viscosity smooths velocity gradients, and the lost ordered mechanical energy becomes heat.

Heat transport. Let $c_v=d\epsilon /dT$ be the thermal heat capacity per molecule, in $\mathrm{J}{\mathrm{K}}^{-1}$. Substituting $a=\epsilon (T)$ gives Fourier's law,Thermal conductivity has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. For a monatomic gas, translational equipartition gives $c_v=3k_B/2$. Rotational or vibrational excitations in other gases can increase the energy carried per molecule, while collisions govern how those modes exchange energy. A temperature gradient can therefore be transported differently from directed momentum even when the dimensional estimates look similar.

Species transport. For a dilute tracer or a mixture with prescribed total density and negligible bulk flow, molecules diffuse down concentration gradients. The first-gradient number flux is Fick's first law,The diffusivity $D$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$; multiplying number flux by the species mass gives mass flux. More generally one uses diffusion flux relative to the mixture's bulk velocity and specifies mass or mole fractions. The approximation is not a claim that an arbitrary pure-gas pressure gradient can remain at rest while only a Fick flux acts. In a mixture, collision cross-sections and relative molecular speeds, including the reduced mass, determine the species diffusivities.

Dependence on density and temperature. The kinetic estimates of gas transport coefficients predict, for fixed collision diameter and molecular heat capacity,Thus dynamic viscosity and thermal conductivity are approximately density independent in the dilute continuum regime. Reducing density reduces the number of carriers but increases their mean free path in the same proportion. At fixed density all three estimates grow as $T^{1/2}$; at fixed pressure, $n=p/(k_{B}T)$ andKinematic viscosity and thermal diffusivity also grow roughly as $T^{3/2}/p$ at fixed pressure in this model. Real interaction potentials make effective collision cross-sections temperature dependent, and internal excitations change $c_v$. Dense-gas correlations and boundaries invalidate the simple cancellation. In particular, when the mean free path approaches the apparatus scale, transport becomes free molecular; it cannot be described by extending the same local, pressure-independent Fourier and viscosity coefficients indefinitely.

Heat diffusion equation. For a stationary gas with constant volumetric heat capacity $c=n{c}_v$, constant thermal conductivity and negligible compression work, the energy balance in any fixed volume isHere $H$ is a possible heat generation rate per volume. The divergence theorem and Fourier's law yieldThermal diffusivity $\alpha$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$ and is of order $\overline{v}\lambda /3$. In a moving or compressing gas, the material derivative, pressure-volume work and viscous heating must also enter the internal-energy equation. The displayed pure diffusion equation describes the approximation in which those effects are absent or negligible.

Representative solutions. With no source, a steady one-dimensional slab between temperatures $T_a,T_b$ at $x=0,L$ has $T=T_a+(T_b-T_a)x/L$ and uniform heat flux $-\kappa (T_b-T_a)/L$. Thus a constant spatial gradient transports heat without changing local temperature in time when the boundaries continually supply and remove it.

For an infinite medium and an initial sinusoidal temperature perturbation,Short spatial wavelengths decay faster, and the characteristic smoothing length after time $t$ is $\sqrt{\alpha t}$. A finite slab with both ends held at $T_b$ has a sine expansion with factors $e^{-\alpha (n\pi /L{)}^{2}t}$; its slowest relaxation time is $L^2/({\pi }^2\alpha )$.

An instantaneous heat input $Q_A$ per unit transverse area at $x=0$ in an otherwise uniform infinite medium gives the Gaussian solutionIts integral $\int c(T-T_b)dx=Q_A$ conserves the deposited heat. The peak falls as $t^{-1/2}$ while its width grows as $\sqrt{t}$; diffusion redistributes energy instead of transporting a sharp front at a single velocity.

Finally, a semi-infinite body initially at $T_i$ and held at surface temperature $T_s$ from $t=0$ hasIt has $T(0,t)=T_s$ and approaches $T_i$ far from the boundary at finite time. These solutions demonstrate steady conduction, mode decay, spreading of localized heat and boundary-driven diffusion. Their idealized behaviour at extremely short times or arbitrarily long distances must be interpreted within the continuum regime rather than as an assertion of instantaneous microscopic communication.

The second virial coefficient is the leading interaction correction to the ideal gas law: at fixed temperature, $B_2(T)={\lim }_{n\to 0}(Z-1)/n$, where $Z=p/(n{k}_{B}T)$ and $n$ is number density. For identical classical particles in a homogeneous three-dimensional gas, with an isotropic pair potential $u(r)$ and $r$ the centre-to-centre separation,The function $f$ is the Mayer f-function. The factor $4\pi r^{2}dr$ counts the volume of a thin spherical shell; the factor $1/2$ counts each unordered pair once. The result has units of volume. Repulsive regions, $u>0$, contribute positively, while attractive regions, $u<0$, contribute negatively. ^[ref]

This is a striking bridge from a microscopic pair potential to a macroscopic pressure reading. At sufficiently small number density, $B_2>0$ raises pressure relative to an ideal gas at the same temperature and number density, while $B_2<0$ lowers it. A single value of $B_2$, however, cannot determine the entire pair potential: many different functions have the same weighted integral. Measurements over several temperatures supply more information.

The integral must exist. For example, a bounded short-range attraction outside a hard core has no difficulty. If instead $u(r)\sim A/r^s$ as $r\to \infty$, then $1-e^{-u/(k_{B}T)}\sim u/(k_{B}T)$ and the radial tail behaves like $\int r^{2-s}dr$: it converges only for $s>3$. A bare unscreened Coulomb potential therefore falls outside this formula's ordinary short-range setting. This condition on the tail does not by itself establish the stability or convergence of the full virial expansion.

For an exactly integrable toy attraction, choose the following square-well pair potential, with diameter $\sigma >0$, attraction depth $ϵ>0$, and range multiplier $a>1$:Set $T^{\ast }=k_{B}T/ϵ$ and $b_{\mathrm{HS}}=2\pi {\sigma }^3/3$, the second virial coefficient of a hard-sphere gas. Integrating the core and attractive shell separately gives the exact toy resultThe core contributes $b_{\mathrm{HS}}$; the shell contributes $-b_{\mathrm{HS}}(a^3-1)(e^{1/T^{\ast }}-1)$. Heating weakens the relative Boltzmann factor enhancement of attractive configurations, so $B_2$ approaches the positive hard-core value as $T^{\ast }\to \infty$. For the illustrative choice $a=3/2$, the zero occurs atAt this temperature, $Z=1+B_3{n}^2+\cdots$: cancellation removes the linear correction, not every correction. At low $T^{\ast }$, the negative shell contribution becomes large, requiring a correspondingly smaller number density for a useful truncated virial expansion. The square well is a deliberately chosen model, not a fit to the carbon dioxide experiment.