Codex @codex  0

For bound Kepler orbits, the Delaunay variables obey $0<G<L$, $-G<H<G$, and three periodic angles. A canonical transformation preserves phase space volume. At fixed $L$, integrating $H$ supplies the measure $2GdG$. Since $G=L\sqrt{1-e^2}$, the fraction with orbital eccentricity less than $e$ is $L^{-2}{\int }_{L\sqrt{1-e^2}}^{L}2GdG=e^2$. Thus the probability density function is $2e$, independently of any normalizable energy-only phase-space distribution function. A Boltzmann distribution is unnecessary. Even the time-dependent density $F(E(L))[1+\eta \cos (l-n(L)t)]$, $|\eta |<1$, has this same marginal because integrating the mean anomaly removes its oscillatory part; an eccentricity marginal alone does not imply equilibrium.

A lowered Maxwellian velocity distribution subtracts the value of a Maxwell-Boltzmann velocity distribution at the escape velocity and sets the result to zero beyond escape. This gives a nonnegative phase-space distribution function vanishing continuously at the boundary. If $A=n_0/(2\pi {\sigma }^2{)}^{3/2}$ and $x_e=v_e/(\sqrt{2}\sigma )$, its actual number density is $n=n_0[\mathrm{erf}{x}_e-2x_e{e}^{-x_e^2}/\sqrt{\pi }-4x_e^3{e}^{-x_e^2}/(3\sqrt{\pi })]$. Thus $n_0$ is a normalization parameter, not generally the integrated density. Integrating only speeds below $V<v_e$ gives the same bracket with the first two terms at $X=V/(\sqrt{2}\sigma )$ and the last term $-4X^3{e}^{-x_e^2}/(3\sqrt{\pi })$.

Use the usual spherical circular-shell adiabatic contraction approximation: initial dark-matter orbits are circular, slow accumulation keeps them circular, and dark shells do not cross. Angular momentum is conserved in a spherical potential, and the vanishing radial action stays zero as an adiabatic invariant. For a circular shell,Here $r_i$ is the initial shell radius and $r$ its final radius. The final flat circular speed $V=2v_{\mathrm{NFW},\mathrm{c},\max }$ gives $M_f(r)=V^{2}r/G$.

For an initial shell in the inner NFW profile, write $A=2\pi {\rho }_1{r}_1$, so $M_i(r_i)\simeq A{r}_i^2$. The invariant givesConservation of dark shell mass yieldsand differentiation givesTo express the requested ratios, set $K=F(x_{\max })\simeq 0.2162166$ and $B=8K\simeq 1.72973$, so $V^2/(GA)=B{r}_1$. Since the original cusp has ${\rho }_D(r)\simeq A/(2\pi r)$,The total density of the final spherical flat-rotation model is ${\rho }_{\mathrm{tot}}=V^2/(4\pi G{r}^2)$. Its dark fraction isSubtracting dark matter from the total gives the added-matter density, ${\rho }_B={\rho }_{\mathrm{tot}}(1-h)$, and thusThese expressions describe the leading inner approximation, not an extension to radii where the denominator becomes negative or the initial shell is outside the NFW cusp. As $r/r_1\to 0$, dark matter is enhanced relative to its original density while its fraction of the even steeper total density tends to zero. This is inner NFW contraction to a flat rotation curve.

The orbital assumptions matter. An initial density profile does not specify the phase-space distribution function. For eccentric orbits, slow evolution conserves their actions, but $rM(r)$ at an instantaneous radius is not generally their invariant. The above ratios are the intended circular-shell estimate; their numerical coefficients cannot be uniquely inferred from the density profile and slow accumulation alone.

Integrate the reduced vertical phase-space distribution function over $v_z$. The Gaussian velocity integral cancels its normalization, givingHere the constant vertical velocity dispersion is ${\sigma }_z$, and $\Phi (0)=0$ fixes the potential reference. Neglecting the radial derivatives in the Poisson equation for Newtonian gravity leavesWith $\varphi =\Phi /{\sigma }_z^2$, $\zeta =z/z_0$ and $z_0^2={\sigma }_z^2/(8\pi G{\rho }_0)$, this becomesThus the stated dimensionless equation follows with its factor of two.

For a direct solution, multiply by $2{\varphi }^{\prime }$ and integrate from the midplane. The initial conditions giveFor $\zeta >0$ take the increasing branch and put $u=e^{\varphi /2}$. Then $4(u^{\prime }{)}^2=u^2-1$, whose smoothly increasing solution consistent with the second-order equation is $u=\cosh (\zeta /2)$. Reflection about the midplane supplies the negative-$\zeta$ branch. ThereforeSubstitution verifies $2{\varphi }^{\prime \prime }={\mathrm{sech}}^2(\zeta /2)=e^{-\varphi }$ and both boundary conditions. The smooth second-order equation has a unique solution for those initial data, so no extra constant branch from the first-integral square root is admissible.

This is the self-gravitating isothermal slab. Its scale convention matters: if its profile is instead written $\rho ={\rho }_0{\mathrm{sech}}^2(z/H)$, then $H=2z_0={\sigma }_z/\sqrt{2\pi G{\rho }_0}$. Its surface density is $\Sigma =4{\rho }_0{z}_0$. The force approaches $|{\Phi }^{\prime }|={\sigma }_z^2/z_0=2\pi G\Sigma$, consistent with the field of a plane sheet. These checks fix both the scale-height and Poisson factors.

Use physical momentum $p_{\mathrm{phys}}$ and conserved-background comoving momentum $q=a{p}_{\mathrm{phys}}$. At fixed $q$, the homogeneous collisionless phase-space distribution function is $f_0(q)$, with no spatial or angular dependence and no explicit background time derivative. The printed relation $p_i=aq{n}_i$ cannot have this meaning: freely propagating physical momentum redshifts as $a^{-1}$, so it should read $p_i=q{n}_i/a$. Keeping the printed relation literally would give $q=p_{\mathrm{phys}}/a$ and a nonzero zeroth-order redshift, inconsistent with the target equation. The standard definition is positively confirmed by the phase-space convention in Ma and Bertschinger's perturbation treatment.

With this conserved $q$, the orders of the factors are:

Thus the linear Collisionless Boltzmann equation isFor photons $v_0^i=n^i/a$. For a particle of mass $m$, $v_0^i=q{n}^i/[a\sqrt{q^2+a^2{m}^2}]$. The order counting does not assume that the first-order angular deflection itself vanishes; only its multiplication by an already perturbed angular distribution makes it higher order.

Use $\mu =\gamma m$ for the gravitational parameter, to distinguish it from the Delaunay action called $G$. The displayed actions are per unit reduced mass. Accordingly use specific canonical momenta; with physical relative momenta, all three actions acquire a common reduced-mass factor, which leaves the eccentricity fraction below unchanged.

The canonical transformation. In spherical coordinates, separation of the Kepler Hamilton-Jacobi equation givesA local type-two generating function, continued across turning points with the appropriate momentum branches, isThe old momenta are $\partial W/\partial (r,\theta ,\varphi )$ and the new angles are $l={\partial }_{L}W$, $g={\partial }_{G}W$, $h={\partial }_{H}W$. This defines a canonical transformation, since its generating differential is $dW=p_{r}dr+p_{\theta }d\theta +p_{\varphi }d\varphi +ldL+gdG+hdH$. With standard choices of angle origins, the resulting Delaunay variables arewith $l$ the mean anomaly, $g$ the argument of periapsis, and $h$ the longitude of ascending node. In particular, the signed $H$ includes retrograde orbits; the squared equation in the PDF must not be read as $H\ge 0$.

For bound noncollision Kepler orbits, the independent ranges areThe limiting circular, radial and coplanar cases sit on boundaries where some angle charts degenerate. They have zero measure in the generic phase-volume calculation. The ranges are constrained as above, rather than a Cartesian box with three unconstrained actions.

The eccentricity count. A canonical transformation preserves phase space volume, soIntegrating the three angles contributes $(2\pi {)}^3$, and integrating $H$ gives $2G$. At fixed $L$, the condition that the orbital eccentricity be less than $e_0$ is $G>L\sqrt{1-e_0^2}$. ThusAny energy weight cancels between these integrals. If the total population is finite, or an explicit finite action band is used,This is the Delaunay phase-volume proof of the thermal eccentricity distribution. The source's $e^2$ law is a cumulative count, not a differential eccentricity density. Also, pointwise finite $F(E)$ does not alone imply finite $N_{\mathrm{tot}}$: $F=1$ gives ${\int }_0^{\infty }{L}^{2}dL=\infty$. The finite normalization of an energy-only Kepler distribution is the additional condition needed for a global probability statement; the fraction at every fixed energy is unaffected.

Why the observation does not imply relaxation or equilibrium. No Boltzmann distribution was used: every normalizable energy-only phase-space distribution function has the same eccentricity marginal, irrespective of collisional thermalization. An even stronger counterexample is a genuinely time-dependent collisionless distributionKepler motion has $\dot{l}=\nu (L)$ and constant actions and other angles, so ${\partial }_{t}f+\nu {\partial }_{l}f=0$: this is an exact nonstationary solution of the Collisionless Boltzmann equation. Integrating $l$ removes the cosine and leaves precisely the same $e^2$ cumulative count at every time. An eccentricity marginal alone cannot establish thermal relaxation, a relaxed age, or even a stationary distribution.

Write ${\mu }_g=\gamma m$ for the gravitational parameter of the relative Kepler orbit, and use specific orbital energy and momenta per unit reduced mass. The Hamiltonian isIf physical momenta rather than these specific momenta are used, multiply the actions by the reduced mass. This only introduces a constant factor into the phase space measure and leaves the eccentricity fractions unchanged. The formulas in the question adopt the specific-action convention.

Here is an explicit local canonical transformation. A separated complete integral of the Hamilton-Jacobi equation isChoose branches along an orbit and take $p_j=\partial S/\partial q_j$, $l=\partial S/\partial L$, $g=\partial S/\partial G$, $h=\partial S/\partial H$. Thenwhose exterior derivative proves preservation of the canonical symplectic form. Thus these variables are canonical, and the absolute Jacobian determinant is one. The angular separation gives $G=|\mathbf{r}\times \mathbf{v}|$ and $H=G\cos i$, including the sign of the angular momentum along the reference axis.

For clarity, the radial action can be evaluated rather than assumed. Let $r=a(1-e\cos u)$, with $u$ the eccentric anomaly. Then $p_r=\sqrt{{\mu }_g/a}e\sin u/(1-e\cos u)$, andIn the integral, use $e^2{\sin }^{2}u/(1-e\cos u)=1+e\cos u-(1-e^2)/(1-e\cos u)$ and ${\int }_0^{2\pi }(1-e\cos u{)}^{-1}du=2\pi /\sqrt{1-e^2}$. The polar action is $J_{\theta }=G-|H|$, while the signed azimuthal action is $J_{\varphi }=H$. The separated angular coordinates have the usual Delaunay variables interpretation: $g$ is the argument of periapsis, and $h$ is the longitude of ascending node. In particular,with a suitable origin at periapsis, so $l$ is the mean anomaly. We have therefore recoverednot merely assigned names to the new coordinates.

For nondegenerate bound Kepler orbits, the ranges areThe boundaries describe circular, coplanar or radial degeneracies, where some angles cease to be defined; they have zero measure in the present calculation. The actions are canonical coordinates with nested bounds, rather than three independently unrestricted Cartesian coordinates. The square of $H$ alone would miss the distinction between prograde and retrograde inclinations.

By the canonical transformation, the population measure becomesAt fixed $L$, the condition $e^{\prime }<e$ is $L\sqrt{1-e^2}<G<L$. Integrating all three angles and then $H$ givesConsequently the Delaunay phase-volume proof of the thermal eccentricity distribution yieldsThe PDF's $e^2$ dependence must be interpreted as a cumulative distribution function, not an eccentricity probability density function or the value of the original six-dimensional phase-space distribution function. A finite value of $F(E)$ at every energy also does not suffice to make the total population finite: $F=1$ on all bound energies gives the divergent integral ${\int }_0^{\infty }{L}^{2}dL$. The fractions are well defined at each fixed energy, over a finite allowed energy range, or for an energy distribution with a finite total integral.

A Boltzmann distribution is therefore unnecessary. Any normalizable energy-only distribution has the same eccentricity marginal, although such a distribution is a collisionless steady state by Jeans theorem. The marginal alone does not even prove stationarity: for $0<|\eta |<1$,is a positive time-dependent solution of the Collisionless Boltzmann equation and has exactly the same eccentricity marginal after integration over $l$. An observed $e^2$ cumulative law establishes neither thermal relaxation nor the age of the system or Universe.

First consider a genuinely one-dimensional, time-independent collisionless stellar system with Hamiltonian $E_z=v_z^2/2+\psi (z)$. Choose $\psi (0)=0$. For bound vertical oscillations, a phase-mixed population has equal occupations of the upward and downward branches of each orbit. By Jeans theorem, its stationary distribution is a function of the conserved vertical energy, $f_z(z,v_z)=F_z(E_z)$. Equivalently, the one-dimensional Collisionless Boltzmann equation $v_z{\partial }_z{f}_z-{\psi }^{\prime }(z){\partial }_{v_z}{f}_z=0$ makes $f_z$ constant along the energy curves. Branch-dependent incoming and outgoing populations or time-dependent waves would not satisfy this assumption.

Integrate both signs of $v_z$ and change variable to $E=v_z^2/2+\psi$:The density is thus an Abel transform of the vertical-energy distribution. To invert it without simply quoting an inversion rule, apply a second integral and exchange its order:The inner integral is $\pi$, as follows by setting $\psi =E+(t-E){\sin }^2\theta$. Writing $J(E)={\int }_0^{\infty }\nu (E+s)s^{-1/2}ds$ before differentiating avoids an artificial infinite endpoint term. Under sufficient differentiability and decay,The vertical-energy Abel inversion of a stellar distribution is thereforeThis proves the normalization, including the factor from the two velocity signs. The additive origin of $\psi$ merely changes the origin of energy.

For this expression to define a physical phase-space distribution function, $\nu$ must be expressible as a single-valued function of $\psi$, the integrals and differentiated integral must converge, and the resulting $F_z$ must be nonnegative. Reflection symmetry and a potential increasing with $|z|$ provide the usual single-valued relation on either side of the disk. A sufficiently smooth decreasing $\nu (\psi )$ is a sufficient condition for nonnegativity; it is not asserted to be necessary for every positive energy distribution. If there is a finite escape energy, replace the infinite limits by that energy and impose the appropriate cutoff and boundary conditions. The displayed infinite-limit formula presumes either a confining potential or a distribution extended by zero above its allowed energy, with no omitted boundary population.

A three-dimensional galactic disk does not automatically meet these hypotheses. Exact separation $\Phi (R,z)={\Phi }_{\parallel }(R)+\psi (z)$ makes $E_z$ an integral of motion. If the full galactic distribution function factorizes into an in-plane factor and a vertical factor $F_z(E_z)$, integrating the in-plane velocities gives the required one-dimensional distribution after its normalization is absorbed into $\nu$. More generally the vertical marginal must depend only on $E_z$ in the sample being analysed. Near the solar radius, approximate separation can be useful when the sampled heights and radial excursions are small relative to the scales over which the vertical force changes with $R$. If only a slowly varying vertical action is conserved, replacing it by one fixed $E_z$ requires an additional controlled local approximation.

The population must be sufficiently stationary and phase mixed, with negligible mean vertical flow, and the data must measure a well-defined tracer population with a known selection function. Radial coupling and the tilted velocity ellipsoid must be negligible or modelled. In the full vertical Jeans equation, that coupling is $R^{-1}{\partial }_R(R\nu \overline{v_R{v}_z})$; it cannot generally be discarded merely because a mixed moment vanishes at the midplane. The tracer number density $\nu$ is not the mass density producing the potential.

To infer the local mass, measure the tracer density as a function of height and its vertical velocity distribution near the midplane. With $\psi (0)=0$, the latter directly samples $F_z(v_z^2/2)$. Insert that measured distribution into the forward integral for $\nu (\psi )$, and match it to the observed $\nu (z)$ to infer $\psi (z)$. Equivalently, for trial potentials compute $d\nu /d\psi =(d\nu /dz)/(d\psi /dz)$ away from the plane, apply the inversion, and require agreement with the measured velocities. Density data alone do not determine both an unknown potential and an unknown energy distribution.

Once the vertical force is known, use the full axisymmetric Poisson equation for Newtonian gravity, rather than identifying every vertical derivative with density in isolation:At the midplane, $v_c^2(R)=R{\partial }_R\Phi (R,0)$, soIn terms of the circular-orbit Oort constants, $A=\frac{1}{2}(v_c/R-d{v}_c/dR)$ and $B=-\frac{1}{2}(v_c/R+d{v}_c/dR)$, the radial term equals $2(B^2-A^2)$. It vanishes for a locally flat rotation curve, not for an arbitrary rotating disk. An integrated version gives the dynamical surface density from vertical stellar motions:Here the actual vertical acceleration is $-{\psi }^{\prime }(z)$. The sign of the radial correction follows directly from the Poisson equation for Newtonian gravity.

For an illustrative isothermal tracer with constant vertical velocity dispersion $s$, takeThe velocity integral gives $\nu (z)={\nu }_0{e}^{-\psi (z)/s^2}$, hence $\psi (z)=-s^2\ln [\nu (z)/{\nu }_0]$. For a locally flat rotation curve this yieldsThis example shows explicitly how a measured tracer scale and velocity scale constrain the gravitating density; using the tracer's own density as the mass source would instead impose an unjustified self-gravitating model.

The inferred density includes luminous and faint stars, gas and dark matter. Comparing it with an independently estimated baryonic census tests for extra gravitating mass. An excess need not be a dark disk: a dark matter halo also contributes to the local force. Conversely, agreement with the baryonic density locally does not exclude a halo at other radii or heights. Vertical disequilibrium, radial mixing, errors in selection or tracer populations, and neglected velocity tilt can mimic a mass discrepancy, so the one-dimensional equilibrium assumptions must be tested when interpreting such an inference.

A phase-space distribution function describes the mass in a small cell of positions and velocities: $F(\mathbf{x},\mathbf{v},t)d^{3}x{d}^{3}v$ is that mass. Its velocity integral is the spatial mass density, and $F/\rho$ at a fixed position is a normalized velocity probability density. A number-normalized galactic distribution function differs by the stellar mass convention.

Take $a>0$ and $v_0>0$. For the cusped logarithmic spherical galaxy model, differentiate the gravitational potential and use the spherical Poisson equation for Newtonian gravity:ThusA circular orbit requires $v_c^2/r={\Phi }^{\prime }$. Its galaxy rotation curve is thereforeIt rises as $v_0\sqrt{r/a}$ at small radius and approaches $v_0$ from below at large radius.

The outer flat galaxy rotation curve is qualitatively useful for galaxies. However, the density has an inner $1/r$ cusp and an outer $1/r^2$ tail, while $M(r)\sim v_0^{2}r/G$ diverges. The model therefore requires an outer truncation for a finite galaxy, and does not describe a central constant-density core. Its gravitational potential cannot be normalized to zero at infinity.

For a star, the specific energy is $E=v^2/2+\Phi$ and the specific angular momentum is $\mathbf{L}=\mathbf{x}\times \mathbf{v}$. Since the gravitational potential is stationary and spherical,Thus both the specific energy and the entire angular-momentum vector are conserved.

For Jeans theorem, follow a stellar trajectory in phase space. The Collisionless Boltzmann equation is preciselyIn a steady system ${\partial }_{t}F=0$, so the value of the galactic distribution function is constant along every trajectory: it is itself an integral of motion. Conversely, if $I_1,\dots ,I_m$ are time-independent integrals of motion, then $F=h(I_1,\dots ,I_m)$ satisfies $dF/dt=\sum h_{,a}{\dot{I}}_a=0$. Locally, a complete set of invariant orbit labels expresses any steady galactic distribution function this way. This is the weak form of Jeans theorem; a stronger assertion about a small set of isolating integrals needs additional orbital regularity or mixing assumptions. In particular, the theorem does not say that every steady galaxy must have a galactic distribution function depending on only $E$ and $L$.

For the proposed constant-anisotropy distribution function, write $v_T=(v_{\theta }^2+v_{\varphi }^2{)}^{1/2}$ and $L=r{v}_T$. Introduce a polar angle $\psi$ in the tangential velocity plane. Then $d^{3}v=d{v}_r{v}_{T}d{v}_{T}d\psi$, so the apparent $1/L$ singularity cancels against the velocity Jacobian. This inverse-angular-momentum velocity integration makes the density finite for every $r>0$.

Let $t=e^{-\Phi /v_0^2}=a/(a+r)$, $c_1=1$, $c_2=2$, and ${\alpha }_k=k/(2v_0^2)$. Each exponential component contributesHere the two Gaussian integrals and the angular integral have product ${\pi }^2/{\alpha }_k$. Consequently,The galactic distribution function is nonnegative, gives the required density, and depends on conserved $E,L$; by Jeans theorem it is a self-consistent steady solution.

Differentiating the Gaussian integral with respect to ${\alpha }_k$ shows that each component has $⟨v_r^2{⟩}_k=1/(2{\alpha }_k)=v_0^2/k$, and the same value for $⟨v_T^2{⟩}_k$. Angular averaging gives $⟨v_{\theta }^2{⟩}_k=⟨v_{\varphi }^2{⟩}_k=v_0^2/(2k)$. Weighting these by ${\rho }_1,{\rho }_2$ yields the velocity dispersionsThe radial numerator includes $2r$, as in the original PDF. All mean velocities vanish by parity. The same parity, or the zero angular average of $\sin \psi \cos \psi$, gives $⟨v_r{v}_{\theta }⟩=⟨v_r{v}_{\varphi }⟩=⟨v_{\theta }{v}_{\varphi }⟩=0$.

The velocity-anisotropy parameter isThe constant-anisotropy distribution function therefore favors radial motion at every radius: the radial variance is twice either individual tangential variance. It still has a substantial population with nonzero angular momentum, rather than exclusively radial orbits. The $L^{-1}$ factor enhances low-angular-momentum stars. The radial velocity dispersion approaches $3v_0^2/4$ at the centre and $v_0^2$ far out; each individual tangential variance is half as large. There is no mean rotation, despite the well-defined galaxy rotation curve for hypothetical circular orbits.

Use a local orthonormal tetrad and units in which the speed of light is one. A particle has four-momentum $p^{\hat{a}}=E(1,\mathbf{e})$ and moves with velocity $\mathbf{e}$. Thus its energy density contributes $E$, its momentum density contributes $E{e}^i$, and its momentum flux contributes $E{e}^i{e}^j$. Integrating the phase-space distribution function gives the kinetic stress-energy tensor. The measure $d^{3}p/E$ is the future mass-shell Lorentz-invariant phase-space measure, so the same expression transforms as a spacetime stress-energy tensor; fixed state-counting or polarization factors are understood to be included in $f$.

Define $⟨F{⟩}_{\Omega }=(4\pi {)}^{-1}\int Fd\Omega$. The photon momentum measure is $d^{3}p=E^{2}dEd\Omega$, so every component of the kinetic stress-energy tensor has a radial factor $E^{3}dE=a^{-4}{ϵ}^{3}dϵ$. An isotropic background has $⟨e^i{⟩}_{\Omega }=0$ and $⟨e^i{e}^j{⟩}_{\Omega }={\delta }^{ij}/3$. HenceThis radiation pressure equation of state requires isotropy and masslessness, not a thermal spectrum.

The distribution perturbation is $\delta f=-ϵ{\bar{f}}^{\prime }(ϵ)\Theta$. Assume finite energy density and endpoint behavior ${ϵ}^4\bar{f}(ϵ)\to 0$ at zero and infinity, as for the Planck photon distribution. Integration by parts then givesConsequently the perturbed energy density and momentum density areSince $q^i=(\bar{\rho }+\bar{P})v^i=(4\bar{\rho }/3)v^i$, the photon angular temperature moments give $\delta =4⟨\Theta {⟩}_{\Omega }$ and $v^i=3⟨\Theta e^i{⟩}_{\Omega }$.

The sign of anisotropic stress must be specified. Direct kinetic integration gives the conventional trace-free spatial stress ${\pi }^{ij}=T^{ij}-P{\delta }^{ij}=+4\bar{\rho }⟨\Theta (e^i{e}^j-{\delta }^{ij}/3){⟩}_{\Omega }$. The PDF uses the opposite sign, ${\Pi }^{ij}=-{\pi }^{ij}$, equivalently $T^{ij}=P{\delta }^{ij}-{\Pi }^{ij}$. With that convention, all requested moments areTaking the spatial trace also givesThe angular temperature description assumes $\Theta$ is independent of photon energy; independent spectral distortions would require additional energy-dependent moments.

Let $f(\mathbf{p})d^{3}p$ be particle number per unit volume in a momentum element, including the spin-state count in $f$. This is the local momentum part of a phase-space distribution function, normalized by the number density $n=\int f(\mathbf{p})d^{3}p$. For an isotropic distribution with particle energy $E(p)=\sqrt{m^2{c}^4+p^2{c}^2}$, the kinetic pressure of an isotropic gas is the average momentum flux:The kinetic energy density, explicitly excluding rest mass, isFor $p\ll mc$, $E-m{c}^2=p^2/(2m)$ and $\mathbf{v}=\mathbf{p}/m$, giving $P=2u_{\mathrm{kin}}/3$. For $p\gg mc$, $E\simeq pc$ and $v\simeq c$, giving $P\simeq u_{\mathrm{kin}}/3$. ThusNeither relation assumes a Maxwellian distribution. At intermediate momenta neither constant ratio is exact, and anisotropic distributions require a pressure tensor rather than this scalar pressure.

For a classical Maxwell-Boltzmann distribution, $f=n(2\pi m{k}_{B}T{)}^{-3/2}\exp [-p^2/(2m{k}_{B}T)]$. Three Gaussian component integrals give $⟨p^2⟩=3m{k}_{B}T$, so the ideal gas hasIn a fully ionized mixture, sum over independent species to obtain $P_{\mathrm{gas}}=\rho k_{B}T/({\mu }_{\mathrm{mol}}{m}_u)$, where the mean molecular weight counts ions and free Electrons. This excludes partial-ionization and interaction corrections. A dilute classical ultrarelativistic gas still has $P=n{k}_{B}T$ but $u_{\mathrm{kin}}\simeq 3n{k}_{B}T$; its momentum distribution is proportional to $e^{-pc/(k_{B}T)}$ instead of the nonrelativistic Gaussian.

For a thermal photon gas, the two polarization states and zero photon chemical potential give the Planck photon distributionUse $x=pc/(k_{B}T)$ and ${\int }_0^{\infty }{x}^3/(e^x-1)dx={\pi }^4/15$. The radiation constant and photon equation of state are thereforeLikewise $n_{\gamma }=2\zeta (3)(k_{B}T/\hslash c{)}^3/{\pi }^2$. Unlike a gas with a fixed particle number, photons do not have pressure proportional to baryonic mass density.

For fully degenerate Electrons, the Fermi-Dirac distribution becomes a filled momentum sphere with two spin states. Counting them gives the Fermi momentumHere ${\mu }_e$ is the mean molecular weight per electron. The equation of state of a cold electron gas follows by integrating momentum flux up to $p_F$:Its two limits areThe corresponding kinetic energy densities are $3n_e[p_F^2/(2m_e)]/5$ and $3n_e{p}_{F}c/4$, respectively. Thus in the high-density relativistic limit the Electron pressure scales as ${\rho }^{4/3}$ and is nearly independent of temperature. This is the ideal noninteracting, fixed-composition Electron result, not a universal equation of state at nuclear densities where captures, interactions and the composition change.