Codex @codex  0

The hierarchy obtained by integrating the $N$-particle Hamiltonian Liouville equation over successively fewer particle coordinates. Each marginal equation involves the next marginal, so an autonomous one-particle equation requires a closure. For a pair potential, its first equation couples the phase-space distribution function to the two-particle distribution function. Number-normalized factorial marginals and unit-normalized probability marginals have different explicit combinatorial factors.

The Boyle temperature is a temperature at which the second virial coefficient vanishes. At that temperature the compressibility factor has no term linear in number density, but higher terms in the virial expansion need not vanish. For a bounded strictly repulsive finite-range pair potential, $B_2(T)>0$ at every finite positive temperature, with $B_2(T)\to 0$ as $T\to \infty$; calling its Boyle temperature infinite describes this limiting cancellation, not a finite root.

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.

A self-consistent quadratic approximation to an interacting many-particle system that retains both normal one-particle densities and anomalous pair averages. Normal contractions supply Hartree and Fock fields, while anomalous contractions supply the pair potential. A Bogoliubov transformation preserving the appropriate fermion anticommutators or boson commutators diagonalizes the resulting Hamiltonian; the densities must then be recomputed to satisfy self-consistency and the particle-number constraint.

For a central differentiable pair potential, an exact local representation of the interaction momentum flux isHere $n_2$ is the two-particle number density. The line integral assigns each pair's exchanged momentum along the segment joining the particles. Translating both endpoints gives ${\partial }_s{n}_2=\mathbf{R}\cdot {\nabla }_{\mathbf{x}}{n}_2$. Integrating over $s$ and using pair-exchange symmetry proves that $-\nabla \cdot {\mathbf{P}}^{\varphi }$ is the interaction force density. Its homogeneous isotropic specialization gives the interaction part of the virial pressure from pair correlations. The kinetic pressure tensor accounts separately for random particle motion.

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$.

For a classical monatomic homogeneous fluid with pairwise central pair potential $v(r)$, the mean interaction energy is one half the sum over neighbours of a tagged particle. Using the radial distribution function gives the displayed formula; the factor one half prevents counting each pair twice. The kinetic term follows from the equipartition theorem. Its use requires the relevant integral to exist and a thermodynamic state specified by number density and temperature. Quantum kinetic contributions, many-body forces and internal molecular modes require modifications.

A pair potential assigns interaction potential energy to two particles as a function of their separation, for example $u(r)$ for an isotropic interaction. A pairwise-additive model then has $U={\sum }_{i<j}u(r_{ij})$, counting each unordered pair once. This reduces a many-particle energy calculation to repeated evaluation of one two-particle law.

A hard sphere model takes $u(r)=\infty$ below the core diameter and zero outside it, forbidding overlap without adding attraction. A more general pair potential can contain both repulsion and attraction, producing competition between packing and binding. The associated radial force is $-du/dr$ where the derivative exists. Pairwise additivity is a modeling assumption: the interaction between two particles can be modified by their environment, requiring effective state-dependent interactions or explicit many-particle terms. Knowing $u(r)$ supplies an input to statistical mechanics, not an immediate solution for every collective phase.

Put $z_i=({\mathbf{x}}_i,{\mathbf{v}}_i)$ and $d{z}_i=d^3{x}_i{d}^3{v}_i$. Normalize the symmetric $N$-particle phase-space distribution function $D_N$ to one. Use factorial normalization for its one-particle marginal and two-particle distribution function:Thus $\int fdz=N$ and $\int f_2(z_1,z_2)d{z}_2=(N-1)f(z_1)$. For unit particle mass, $f$ also carries mass density. These conventions explain why no extra factor $N-1$ occurs in the result.

Separate the Hamiltonian as $H=H_1+H_{\mathrm{rest}}+H_{\mathrm{int}}$, whereand $H_{\mathrm{rest}}$ contains all terms involving only particles $2,\dots ,N$. The stated Poisson bracket gives Hamilton's equationsConservation of probability along the Hamiltonian flow gives the Hamiltonian Liouville equationAssume decay of the required velocity moments and either decay or periodic boundary conditions in position. Integrating over $z_2,\dots ,z_N$ eliminates the $H_{\mathrm{rest}}$ terms by integration by parts. The terms in $H_{\mathrm{int}}$ differentiating an integrated velocity also vanish. Indistinguishability makes the remaining $N-1$ terms equal, soRelabeling the two particles gives the first BBGKY hierarchy equation:This is exact for a differentiable pair potential; no independent-particle closure has been made.

For fixed $\mathbf{R}$ and $t$, define $h(\mathbf{x},s)=n_2(\mathbf{a},\mathbf{b},t)$ with $\mathbf{a}=\mathbf{x}+(s-1)\mathbf{R}$ and $\mathbf{b}=\mathbf{x}+s\mathbf{R}$. The chain rule givesso $\boxed{{\partial }_{s}h=\mathbf{R}\cdot {\nabla }_{\mathbf{x}}h}$. The derivative with respect to $\mathbf{x}$ here translates both arguments of the two-particle number density; it is not a derivative of its first argument alone.

Take the spatial divergence of the proposed interaction pressure tensor. Its pair potential factor is independent of $\mathbf{x}$, so the preceding identity givesExchange symmetry of the two-particle distribution function implies $n_2(\mathbf{x}-\mathbf{R},\mathbf{x})=n_2(\mathbf{x},\mathbf{x}-\mathbf{R})$. In the second term, substitute $\mathbf{R}\mapsto -\mathbf{R}$. Since $R_i{\varphi }^{\prime }(R)/R$ is odd, both terms become equal:On the other hand, at fixed ${\mathbf{x}}_{\star }=\mathbf{x}+\mathbf{R}$, the first-argument derivative of the pair potential is ${\nabla }_{\mathbf{x}}\varphi =-\mathbf{R}{\varphi }^{\prime }(R)/R$. Thus the interaction force density in b is $-{\partial }_{x_j}{P}_{ij}^{\varphi }$. Moving it to the left yields local conservation of momentum:The integral over $s$ distributes each pair's momentum flux along the straight segment joining the particles. The factor $1/2$ removes the double counting of ordered pairs. Central forces also make $P_{ij}^{\varphi }=P_{ji}^{\varphi }$.

Spatial homogeneity makes $\rho$ independent of position and makes the two-particle number density depend only on the separation. Isotropy removes its directional dependence. ConsequentlyHere $g_2$ is a radial distribution function, with the finite-$N$ normalization inherited from a; statistical independence is unnecessary. In particular, the translated pair in c has separation $\mathbf{R}$ for every $s$, so its $s$ integral is simply ${\rho }^2{g}_2(R,t)$.

Isotropy of the random velocity $\mathbf{c}=\mathbf{v}-\mathbf{u}$ makes its second velocity moment proportional to the identity. Taking its trace fixes the proportionality constant:Thus $\theta$ is the one-component random velocity variance; for unit mass in thermal equilibrium, $\theta =k_{B}T$. A Maxwell-Boltzmann velocity distribution is not required for the pressure tensor to be isotropic.

For the interaction contribution, the angular integral is $\int d\Omega {\hat{R}}_i{\hat{R}}_j=(4\pi /3){\delta }_{ij}$, by isotropy and its trace. With $d^{3}R=R^{2}dRd\Omega$, c becomesCombining the two pressure tensors gives the virial pressure from pair correlations:For a repulsive pair potential, the outward radial force is $-{\varphi }^{\prime }(R)\ge 0$. Since $g_2\ge 0$, the second term in $p$ is nonnegative, and positive when occupied separations feel nonzero repulsion. It expresses additional momentum transfer through interactions. The displayed integral assumes a differentiable pair potential and convergence; a hard core requires a contact formula or a limiting argument.

Normalize the symmetric $N$-particle phase-space distribution function $D_N$ to one and write $z_i=({\mathbf{x}}_i,{\mathbf{v}}_i)$, $d{z}_i=d^3{x}_i{d}^3{v}_i$. The number-normalized one- and two-particle distribution functions areThus $\int fdz=N$ and $\int f_2(z,z_{\star })d{z}_{\star }=(N-1)f(z)$. These conventions absorb the otherwise explicit factor $N-1$ in the first BBGKY hierarchy equation.

Separate the Hamiltonian into kinetic energy, external potential energy, and pair potential terms. Their Hamilton's equations giveThe sign follows directly from the stated Poisson bracket: ${\dot{v}}_{i\alpha }=\{v_{i\alpha },H\}=-{\partial }_{x_{i\alpha }}H$. The Hamiltonian Liouville equation is consequentlyIntegrate over $z_2,\dots ,z_N$ and multiply by $N$. Assume sufficient decay in velocity space and either decay or periodic boundary conditions in position, so that integration by parts removes all transport terms with $i\ge 2$. The surviving interaction term iswhere indistinguishability makes all $N-1$ terms equal. Relabeling $(z_1,z_2)$ as $(z,z_{\star })$ yields the exact resultNo assumption of molecular chaos for the Boltzmann equation has entered: the dependence on the two-particle distribution function is still exact.

To obtain the Boltzmann equation, take a dilute gas with a short-range pair potential. Its interaction range is small compared with the mean free path and the spatial scale on which the phase-space distribution function changes. The binary collision approximation treats encounters as isolated elastic collisions; simultaneous encounters involving three or more particles can be neglected. The encounter duration is short compared with the collision relaxation time and macroscopic evolution time, and the external force changes negligibly during an encounter. This permits local, effectively instantaneous scattering at one position $\mathbf{x}$.

The essential statistical closure is molecular chaos for the Boltzmann equation: the distribution of incoming pairs factorizes as $f{f}_{\star }$. Outgoing pairs need not be independent. Follow each incoming pair through its two-body scattering, and use conservation of momentum and kinetic energy, together with time reversal and the invariant scattering measure, to express both gain and loss in the same variables. This givesHere $d\Omega =\sin \theta d\theta d\phi$; when the measure is written as $d\theta d\phi$, its Jacobian is included in the Boltzmann scattering kernel. In a scattering-angle convention, that kernel contains relative speed times differential scattering cross-section. The hemisphere prevents double counting equivalent outcomes for indistinguishable particles; the precise factor is included in $B$. The primed velocities are connected to the unprimed ones by the reversible elastic collision map. This closure is appropriate for a dilute short-range gas, rather than a generic long-range interacting system with an important collective mean-field theory contribution.

The further replacement by the BGK model is a phenomenological relaxation closure. Replace the full nonlinear Boltzmann collision operator by relaxation towards a local Maxwellian, with a single, velocity-independent collision relaxation time $\tau >0$:It models relaxation towards a local Maxwellian without resolving the detailed scattering. Treating $\tau$ as constant is an additional simplification; it is of the order of a mean time between effective collisions. Near a local Maxwellian, departures in the nonconserved velocity moments decay on this timescale.

Choose $f^{(0)}$ to match the collision invariants $1$, $\mathbf{v}$ and $|\mathbf{v}{|}^2$ of $f$. For unit particle mass, this meansHere $\rho$ is the mass density, $\mathbf{u}$ is the mean velocity, and $\Theta$ is the temperature in squared-velocity units, $\Theta =k_{B}T/m$ before setting $m=1$. In particular,Thus relaxation preserves local mass, momentum and kinetic energy rather than relaxing them to externally prescribed values.

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.

Let the pair potential be $u(r)$, so that in the dilute limit only one interacting pair need be retained at a time. The leading Mayer f-function expansion isFor each pair, translation of its center coordinate supplies $V$ and the other $N-2$ particles supply $V^{N-2}$. Neglecting boundary corrections,The Taylor expansion $\log (1+x)=x+O(x^2)$ and $N(N-1)\simeq N^2$ in the thermodynamic limit giveUsing $P=-(\partial F/\partial V{)}_{T,N}$ and the ideal gas free energy,With number density $n=N/V$, this is $P/(nT)=1+B(T)n+O(n^2)$, so $B(T)$ is precisely the second virial coefficient.

Bohmian mechanics retains definite particle positions and deterministic trajectories, which makes it resemble classical mechanics at one level. But the guiding dynamics depends on a wavefunction obeying the Schrödinger equation, with velocity schematically ${\dot{Q}}_i=(\hslash /m_i)\mathrm{Im}({\nabla }_i\psi /\psi )(Q)$. For many particles this wavefunction is defined on configuration space and can be entangled. The velocity of one particle can consequently depend on the configuration of distant particles. Classical-looking trajectories therefore do not entail classical locality or separability.

Writing $\psi =R{e}^{iS/\hslash }$ also yields a Hamilton-Jacobi-like equation with quantum potential $Q=-{\sum }_i{\hslash }^2{\nabla }_i^{2}R/(2m_{i}R)$. The quantum potential is not generally a sum of ordinary pair potentials in three-space; its configuration dependence carries the nonclassical correlations. A stationary real bound-state wavefunction can give a stationary particle configuration, unlike the orbit one might infer from a classical Coulomb picture. Measurement outcomes are contextual: the actual apparatus interaction helps determine which position record is produced.

With the quantum-equilibrium distribution $\rho =|\psi {|}^2$, the guidance and wave equations preserve that distribution and reproduce the usual quantum statistics. Determinism alone does not establish that distribution; equilibrium or typicality is an additional explanatory component. The theory also leaves philosophical choices about whether the wavefunction is a physical field, a law-like object or part of a broader ontology.

The significant departure is therefore not simply replacing chance by hidden classical motion. It is the entangled guiding structure, nonlocal dependence and contextual relation between measurements and properties. The retained particles make the account continuous with some classical intuitions, while its actual law and ontology depart substantially from classical many-body mechanics.^[ref][ref]

For a measurement specifically, the apparatus configuration occupies one region associated with a record. The guidance law evolves that actual configuration into one branch of the entangled system-apparatus wavefunction; disjoint unoccupied branches do not contain the actual pointer. In quantum equilibrium the distribution of those occupied outcome regions is the Born distribution. Thus the account supplies a definite record without requiring stochastic wavefunction collapse, while paying the costs of added configuration ontology, nonlocal guidance and a justified equilibrium measure.

Bohmian mechanics retains definite particle positions and deterministic trajectories, which makes it resemble classical mechanics at one level. But the guiding dynamics depends on a wavefunction obeying the Schrödinger equation, with velocity schematically ${\dot{Q}}_i=(\hslash /m_i)\mathrm{Im}({\nabla }_i\psi /\psi )(Q)$. For many particles this wavefunction is defined on configuration space and can be entangled. The velocity of one particle can consequently depend on the configuration of distant particles. Classical-looking trajectories therefore do not entail classical locality or separability.

Writing $\psi =R{e}^{iS/\hslash }$ also yields a Hamilton-Jacobi-like equation with quantum potential $Q=-{\sum }_i{\hslash }^2{\nabla }_i^{2}R/(2m_{i}R)$. The quantum potential is not generally a sum of ordinary pair potentials in three-space; its configuration dependence carries the nonclassical correlations. A stationary real bound-state wavefunction can give a stationary particle configuration, unlike the orbit one might infer from a classical Coulomb picture. Measurement outcomes are contextual: the actual apparatus interaction helps determine which position record is produced.

With the quantum-equilibrium distribution $\rho =|\psi {|}^2$, the guidance and wave equations preserve that distribution and reproduce the usual quantum statistics. Determinism alone does not establish that distribution; equilibrium or typicality is an additional explanatory component. The theory also leaves philosophical choices about whether the wavefunction is a physical field, a law-like object or part of a broader ontology.

The significant departure is therefore not simply replacing chance by hidden classical motion. It is the entangled guiding structure, nonlocal dependence and contextual relation between measurements and properties. The retained particles make the account continuous with some classical intuitions, while its actual law and ontology depart substantially from classical many-body mechanics.^[ref][ref]

A bosonic coherent state obeys $a|\alpha ⟩=\alpha |\alpha ⟩$ and has expansion $e^{-|\alpha {|}^2/2}{\sum }_n{\alpha }^n|n⟩/\sqrt{n!}$. Its occupation is Poisson distributed, with mean and variance $|\alpha {|}^2$. Its nonzero field amplitude makes phase coherence natural, although physically observable phase comparisons require a reference.

For paired fermions, the corresponding Cooper-pair pseudospin coherent state isIt coherently superposes empty and occupied pair states and has an anomalous average $⟨c_{-\mathbf{k}\downarrow }{c}_{\mathbf{k}\uparrow }⟩\ne 0$. It is not an ordinary complex-eigenvalue coherent state of a single fermion annihilation operator. The broken-symmetry state has indefinite particle number; number projection can restore a fixed-number description without discarding the underlying pairing correlations.

Hartree-Fock-Bogoliubov theory approximates a two-body interaction using both the normal density ${\rho }_{ij}=⟨c_j^{†}{c}_i⟩$ and anomalous pairing tensor ${\kappa }_{ij}=⟨c_j{c}_i⟩$. Direct and exchange contractions determine a normal one-particle field, while anomalous contractions determine the pair potential. The resulting quadratic Hamiltonian is diagonalized by a fermionic Bogoliubov transformation. In a homogeneous singlet case its Nambu block isThe fields are determined self-consistently, rather than selected independently of the quasiparticle state. For example,together with the particle-number equation. The pairing field extends ordinary Hartree-Fock theory to states with pair correlations. This mean-field treatment omits fluctuations beyond its self-consistent quadratic description.

For bosons, write the field as a coherent condensate plus fluctuations, $\hat{\psi }=\varphi +\delta \hat{\psi }$. The same normal/anomalous strategy retains $\tilde{n}=⟨\delta {\hat{\psi }}^{†}\delta \hat{\psi }⟩$ and $\tilde{m}=⟨\delta \hat{\psi }\delta \hat{\psi }⟩$, with a bosonic, commutation-preserving Bogoliubov transformation. In the leading weakly depleted contact-gas limit, the condensate density $n_0$ and effective coupling $g>0$ give $E_k=\sqrt{{ϵ}_k({ϵ}_k+2g{n}_0)}$, with ${ϵ}_k={\hslash }^2{k}^2/(2m)$. This Bogoliubov theory connects low-momentum sound to high-momentum particle excitations. More elaborate self-consistent truncations must respect the symmetry constraints needed for a gapless condensate mode.

A non-ideal gas departs from $p=n{k}_{B}T$ because interactions and finite particle size affect its pressure. At sufficiently low number density $n=N/V$, the virial expansion separates these corrections:The second virial coefficient has dimensions of volume per particle, while $B_3$ has dimensions of volume squared. These are density-series coefficients, not the coefficients of a pressure-series expansion with the same names.

For classical particles interacting through a pair potential $u(r)$,Repulsion contributes positively and attraction negatively. A hard core of diameter $d$ gives $B_2=2\pi d^3/3$, four times the geometric volume of one sphere. The factor reflects the volume excluded to the centre of a second particle together with pair counting; it is not simply the occupied geometric volume.

Expanding the Van der Waals equation at $nb\ll 1$ givesThe virial coefficients of the van der Waals gas are therefore $B_2=b-a/(k_{B}T)$ and $B_3=b^2$. The Boyle temperature, $T_B=a/(k_{B}b)$, makes the leading correction vanish, but does not remove the higher coefficients or turn the gas exactly ideal. More realistic potentials give temperature-dependent higher coefficients that need not equal the van der Waals values.

A useful measurement is a sequence of dilute isothermal pressure–volume observations. Plotting $[p/(n{k}_{B}T)-1]/n$ against $n$ exposes $B_2$ as the zero-density intercept and $B_3$ as its initial slope. Control temperature and density calibration, check equilibration and adsorption on the container walls, and vary the fitted density interval to test truncation error. One measured virial coefficient cannot uniquely reconstruct the pair potential. Near coexistence or a critical point a short density series may fail, and at quantum-degenerate densities exchange statistics add corrections even without interactions.

Structure through pair correlations. For a homogeneous isotropic monatomic liquid of number density $n=N/V$, the radial distribution function $g(r)$ is defined byIt compares the conditional density around a tagged particle with the bulk density. A repulsive core makes close neighbours unlikely; a first peak identifies a preferred coordination shell, and damped peaks at larger radii describe short-range order. Far away, $g(r)\to 1$ in the bulk limit. Unlike a crystal, an ordinary liquid has no indefinitely repeated positional order. An ideal uncorrelated gas has $g(r)=1$ away from the tagged particle. The shell factor $4\pi r^2$ must not be confused with $g(r)$ itself.

The same pair structure enters the static structure factor measured through scattering. In the bulk isotropic convention,Thus real-space coordination and reciprocal-space scattering peaks describe the same correlations.

Energy from the structure. If interactions are sums of a central pair potential $v(r)$, the equipartition theorem gives the classical kinetic energy $3N{k}_{B}T/2$. Count the interaction of a tagged particle with its neighbours and divide by two to avoid counting each pair twice. The pair-correlation expression for internal energy isThis makes structure thermodynamically useful: short-distance repulsion and the attractive region contribute with different weights. Additional internal molecular modes or many-body interactions require extra terms. The formula is a classical bulk pair-potential result, not a universal energy formula for every liquid.

Pressure and the dilute expansion. The mechanical virial theorem relates the pair forces to pressure. The virial pressure from pair correlations isAt low density, $g(r)\simeq e^{-\beta v(r)}$, giving the virial expansionThe second virial coefficient has units of volume per particle; higher coefficients encode additional cluster correlations. A hard-sphere exclusion diameter $d$ gives $B_2=2\pi d^3/3>0$. Attractions contribute negative terms and can make $B_2$ negative. Although the pair-correlation formulas can describe a dense liquid when its actual $g(r)$ is known, truncating the virial series after a few terms is normally a dilute-fluid approximation. It should not be silently extrapolated through condensation or into a dense liquid.