Codex @codex  0

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.

For a spherical, nonstreaming stellar population with isotropic velocity dispersion, each component of the random pressure tensor equals $p=n⟨v_r^2⟩$. The observed second stellar velocity moment satisfies $N(R){\sigma }_{\mathrm{los}}^2(R)=2{\int }_R^{\infty }p(r)r/\sqrt{r^2-R^2}dr$. Composing this Abel transform with itself gives the displayed inverse under sufficient decay and regularity. The Spherical Jeans equation then gives ${\psi }^{\prime }=p^{\prime }/n$, with gravitational acceleration $\nabla \psi$. The Poisson equation for Newtonian gravity recovers total gravitating mass density, not just the tracer density. Without velocity isotropy the same projected moment also contains anisotropic stresses.

In Gaussian units, the equation in the title uses the gyrotropic pressure tensor $\mathbf{P}=p_{\perp }\mathbf{I}+(p_{\parallel }-p_{\perp })\mathbf{bb}$. This allows pressure anisotropy. Its magnetic force contains both magnetic pressure and magnetic tension.

A magnetized long-wavelength plasma description with a common leading perpendicular flow, magnetic flux freezing and a gyrotropic pressure tensor. The kinetic magnetohydrodynamic momentum equation couples to a drift-kinetic equation for each species; parallel plasma pressure and perpendicular plasma pressure are calculated as velocity moments, rather than prescribed by an isotropic adiabatic closure. Relative species drifts enter at higher order in the ratio of gyroradius to macroscopic length.

Use the source's Gaussian units and neglect the displacement current, as in nonrelativistic magnetohydrodynamics. Write $\Delta p=p_{\perp }-p_{\parallel }$ for the pressure anisotropy. The pressure tensor is $\mathsf{P}=p_{\perp }\mathsf{I}-\Delta p\mathbf{bb}$, so its force density isAmpère's law and the solenoidal magnetic-field constraint give the Lorentz force decompositionAdding these two force densities and using $\mathbf{BB}=B^2\mathbf{bb}$ yieldsThe first term is the gradient of the modified magnetohydrodynamic total pressure; the second is the anisotropic Maxwell stress tensor contribution together with the directional pressure force. Their effective coefficient is the pressure-anisotropy-modified magnetic tension, $\mathcal{T}=B^2/(4\pi )+\Delta p$. At fixed magnetic field, positive pressure anisotropy increases resistance to bending, while negative pressure anisotropy reduces it. When $\mathcal{T}$ becomes negative, bending experiences an amplifying force rather than a restoring force. Gradients of $\Delta p$ can also drive motion, so the curvature interpretation refers specifically to the bending contribution.

Define mean streaming velocity $u_i=\overline{v_i}$ and the velocity dispersion tensor ${\sigma }_{ij}^2=\overline{(v_i-u_i)(v_j-u_j)}=\overline{v_i{v}_j}-u_i{u}_j$. The zeroth velocity moment of the steady Collisionless Boltzmann equation gives the continuity equation ${\sum }_i{\partial }_i(n{u}_i)=0$. ThereforeSubstituting in a yields the standard steady Jeans equations:The unsteady version adds ${\partial }_t{u}_j$ on the left. Defining the stellar pressure tensor $P_{ij}=\rho {\sigma }_{ij}^2$ makes the analogy with fluid momentum balance explicit. The off-diagonal terms are correlations between random velocity components; they need not vanish or equal the diagonal terms.

The Jeans equations express local momentum conservation for a collisionless stellar system, just as the Euler equations do for a gas. The mean streaming velocity plays the role of fluid flow, the gravitational force is a body force, and random stellar motions supply a pressure tensor through $\rho {\sigma }_{ij}^2$. Taking moments smooths over individual stellar orbits and produces fields of density, mean velocity and random-motion stress; this is the sense of “stellar hydrodynamics”.

Ordinary gas pressure is usually nearly isotropic because collisions relax its local velocity distribution. A collisionless stellar system can retain an anisotropic velocity ellipsoid, so its pressure tensor is not generally a scalar times the identity. There are more unknown second moments than momentum equations, and their evolution introduces higher moments. The Jeans equations therefore require further dynamical assumptions or a specified one-particle distribution function for closure; they do not by themselves establish local thermal equilibrium or a gas equation of state.

Define the local mean velocity and momentum flux byHere $\rho =\int f{d}^{3}v$ and $\mathbf{u}$ is defined where $\rho >0$. The pressure tensor separates random motion from the bulk velocity:Integrating the mean-field alignment equation over velocity space removes the force divergence and gives the continuity equationIts first velocity moment follows by multiplying by $v_i$ and using $\int v_i{\partial }_{v_j}(f{L}_j)d^{3}v=-\int f{L}_i{d}^{3}v$:Inserting the mean-field alignment equation's force and integrating each particle velocity separately givesThe right-hand side exchanges momentum between positions and aligns their bulk velocities. Its spatial integral vanishes by exchange symmetry, recovering total momentum conservation. No closure for the pressure tensor is needed for these exact velocity moment identities.

Take the zeroth velocity moment of a. The right side is a velocity space divergence with vanishing boundary contribution, giving the continuity equationFor the first velocity moment, use integration by parts componentwise:Write $\mathbf{c}=\mathbf{v}-\mathbf{u}$; then $\int \mathbf{c}f{d}^{3}v=0$, and the momentum flux splits into bulk transport and the pressure tensor:Thereforewhere repeated Cartesian indices are summed and $n_2=\int f_2{d}^{3}v{d}^3{v}_{\star }$. The negative sign is the attractive or repulsive force $-\nabla \varphi$, rather than an additional sign convention for the pressure tensor.

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.

Use the BGK model from b. Multiply it by $v_i$ and integrate over velocity space. Matching the collision invariants makes the collision contribution vanish. For the external force term, integration by parts givesDefine the momentum flux and force density byThe resulting momentum conservation equation isThe pressure tensor is the second central velocity moment,The zeroth velocity moment also supplies the continuity equation ${\partial }_t\rho +{\partial }_j(\rho u_j)=0$, which will be needed in d.

Write $c_i=v_i-u_i$ and introduce the raw third velocity moment $M_{ijk}$ and the third central velocity moment $Q_{ijk}$:The expansion of their product givesThe heat flux is $q_i=Q_{ijj}/2$ for unit mass.

Multiply the BGK model by $v_i{v}_j$. The external force term isTherefore the raw second velocity moment satisfiesThe matching local Maxwellian has ${\Pi }_{ij}^{(0)}=\rho u_i{u}_j+\rho \Theta {\delta }_{ij}$, soTo remove the bulk motion, let $D_t={\partial }_t+u_k{\partial }_k$ be the material derivative. The continuity equation and the momentum conservation equation implyand henceSubtract this identity from the raw second velocity moment equation and insert the expansion of $M_{ijk}$. The terms $u_i{\partial }_k{P}_{jk}$ and $u_j{\partial }_k{P}_{ik}$ cancel by the product rule. The external force terms also cancel, since a common acceleration changes the local mean velocity without directly changing the relative velocities. What remains isThe trace of the relaxation term is zero because $P_{ii}=3\rho \Theta$. Thus the BGK model relaxes anisotropy of the pressure tensor while preserving the local random kinetic energy. The evolution is still a velocity-moment hierarchy, since $Q_{ijk}$ requires further information about $f$.

The pressure tensor has Cartesian components $p_{ab}=\int f(v_a-u_a)(v_b-u_b)d^{3}v$. With zero mean velocity and an isotropic galactic distribution function, off-diagonal components vanish by reflection symmetry and diagonal components are equal by rotational symmetry. ThereforeHere $p$ is a tracer number-weighted random-motion stress, rather than thermodynamic pressure per volume. Its tensor divergence is $(\mathrm{div}\mathbf{p}{)}_a={\partial }_b(p{\delta }_{ab})={\partial }_{a}p$. The time-independent, nonstreaming Jeans equation consequently reduces toThe gravitational acceleration convention is $\nabla \psi$: this $\psi$ differs by a sign from the usual Newtonian gravitational potential $\Phi$.

Let $z$ measure distance along the line of sight. Spherical density projection givesBecause isotropy makes the local second stellar velocity moment along every line of sight equal to $p/n$, and because there is no streaming motion, the projected moment isThis is the forward relation for isotropic stellar pressure deprojection.

Both inversions follow from one Abel transform calculation. Write $F(u)=N(\sqrt{u})$ and $a(v)=n(\sqrt{v})$, so $F(u)={\int }_u^{\infty }a(v)(v-u{)}^{-1/2}dv$. Compose the transforms, using Tonelli's theorem for nonnegative densities, or absolute convergence for signed profiles:Assume sufficient decay to make these integrals finite and sufficient regularity for the last derivative. Then $J^{\prime }(x)=-\pi a(x)$. Setting $x=r^2$ gives the spherical Abel deprojectionReplacing $F(u)$ by $N(\sqrt{u}){\sigma }^2(\sqrt{u})$ and $a(v)$ by $p(\sqrt{v})$ provesThese are the requested expressions, with $u=R^2$. The composed-transform argument avoids an invalid separate differentiation of a divergent lower-end kernel.

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.

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.

The boundaries on a stellar equation-of-state regime diagram concern the electron component. Define the electron relativistic density threshold and thermal electron relativistic threshold byA degenerate electron gas is nonrelativistic well below ${\rho }_{\ast }$ and ultrarelativistic well above it: this is an approximately vertical division on a log-density plot. A nondegenerate electron gas instead becomes thermally relativistic near $T_{\ast }$: this is an approximately horizontal division. The regimes have broad crossovers rather than a discontinuity at either line.

The exact electron Fermi temperature, subtracting electron rest energy, isFor electrons with negligible thermal pairs, strong degeneracy requires $T\ll T_F$, whereas $T\gg T_F$ gives a nondegenerate gas. In the pair-rich regime the actual electron and Positron distributions must instead be treated with their chemical potentials, as discussed below. The two limiting degeneracy boundaries areThese slopes $2/3$ and $1/3$ explain the bent degeneracy boundary on the logarithmic graph. A thermal-wavelength test $n_e[2\pi {\hslash }^2/(m_e{k}_{B}T){]}^{3/2}\sim 2$ gives the same nonrelativistic density/temperature scaling, with an order-one definition of the crossover.

The radiation-to-gas pressure boundary in the nondegenerate fully ionized regime isRadiation dominates above this line. Once electrons are strongly degenerate, compare radiation pressure with $P_e(\rho )$ instead of continuing the ideal-gas comparison into that regime. The radiation-to-degeneracy pressure boundary is $T=[3P_e(\rho )/a_{\mathrm{rad}}{]}^{1/4}$, with logarithmic slopes $5/12$ and $1/3$ in the nonrelativistic and ultrarelativistic limits. Degeneracy of the electrons and dominance of their pressure are distinct criteria.

For the pair curve, distinguish baryonic net electrons from thermally created electrons and Positrons. Chemical equilibrium with photons requires opposite electron/Positron chemical potentials when the one-particle energy includes rest energy, as in the distributions below. At low temperature and low degeneracy the zero-chemical-potential density per charge species isCharge neutrality gives $n_{-}-n_{+}=n_b=\rho /({\mu }_e{m}_u)$, while the nondegenerate equilibrium product is $n_{-}{n}_{+}=n_0^2$. HencePairs become important when $n_0$ is comparable with or exceeds $n_b$, approximately the electron-positron thermal pair abundance curve ${\rho }_{\mathrm{pair}}(T)={\mu }_e{m}_u{n}_0(T)$. Below its density at a fixed temperature, pairs dominate over the net charge electrons. The exponential makes the low-temperature portion steep in a log-log plot. For ${\mu }_e=2$, the full zero-potential integral gives approximate pair-marker densities $1.80\times 10^3\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$ at $10^9\mathrm{K}$ and $7.52\times 10^5\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$ at $3\times 10^9\mathrm{K}$. At relativistic temperature use the full zero-chemical-potential Fermi-Dirac distribution instead; it gives $n_0\to [3\zeta (3)/(2{\pi }^2)](k_{B}T/\hslash c{)}^3$, so the high-temperature pair curve approaches slope three in log-density versus log-temperature. Strong net-electron degeneracy suppresses Positrons and requires the full chemical-potential-dependent distribution.