Codex @codex  0

Take a homogeneous comoving cell of physical volume $V_c=a^3{V}_0$ and total entropy $S=s{V}_c$. At vanishing chemical potential, the reversible first law of thermodynamics readsFor an extensive equilibrium system, the thermodynamic Euler relation is $E=TS-P{V}_c+\mu N$. At $\mu =0$ this givesThe supplied pressure identity verifies consistency directly: $dP=(\rho +P)dT/T$ impliesThe thermodynamic Euler relation fixes the additive normalization that a differential identity alone would leave undetermined.

For adiabatic cosmological expansion, the cosmological perfect-fluid continuity equation gives $\dot{\rho }+3H(\rho +P)=0$, while ${\dot{V}}_c=3H{V}_c$. HenceThis proves cosmological entropy conservation in the stated perfect-fluid approximation. Kinetic equilibrium alone does not guarantee adiabatic evolution: entropy injection, dissipative stresses or out-of-equilibrium reactions would require additional source terms. Those effects are implicitly excluded in this calculation.

During the cannibal phase of a decoupled sector, cosmological entropy conservation and the leading result of part (c) giveDifferentiating and using $T_{\varphi }=m_{\varphi }/x$ yields logarithmic cooling during cannibalism:Here $a_{\ast }$ absorbs the entropy normalization, and the logarithmic expression is a large-$x$ asymptotic relation. With fixed visible entropy degrees of freedom the photon bath has $T_{\gamma }\propto a^{-1}$. Thus for $x\gg 1$ the scalar particles cool much more slowly than the photons and become relatively hotter.

The physical reason is that $3\varphi \to 2\varphi$ converts rest-mass energy into kinetic energy of the remaining particles. Elastic scattering distributes this energy, offsetting much of the cooling caused by expansion. The total comoving entropy remains constant; the number density does not merely dilute as $a^{-3}$. This explanation applies only while number-changing chemical equilibrium persists. If those reactions freeze out, an isolated number-conserving nonrelativistic gas instead has $T_{\varphi }\propto a^{-2}$ under adiabatic expansion, which is faster cooling than the photons. The paper supplies no rates, so that transition time cannot be determined.

Work first in units $G=c=\hslash =k_B=1$. Black-hole thermodynamics starts with the laws of black-hole mechanics. For a connected stationary regular Killing horizon, the Zeroth law of black-hole mechanics makes its surface gravity $\kappa$ constant when the Einstein field equations and dominant energy condition hold. This parallels uniform equilibrium temperature. For neighboring stationary asymptotically flat four-dimensional Einstein-Maxwell black holes, the First law of black-hole mechanics iswhere ${\Omega }_H$ is the horizon angular velocity and ${\Phi }_H$ its electric potential relative to infinity. The last terms are rotational and electromagnetic work, analogous to the work terms in $\delta E=T\delta S+\text{work}$. The second law of black-hole mechanics, expressed by Hawking's area theorem, makes the total future horizon area nondecreasing for classical matter obeying the null energy condition and global assumptions such as strong asymptotic predictability. The third law of black-hole mechanics is an unattainability statement: subject to the usual regularity assumptions, the weak energy condition, and a bounded stress-energy tensor, no finite physical process reduces the surface gravity of a regular horizon to zero. It does not assert that extremal black holes have zero area or zero entropy.

Classical geometry alone fixes an analogy, not a nonzero physical temperature. Quantum field theory supplies Hawking temperature $T_H=\kappa /(2\pi )$. Comparing the area term of the first law with $T_H\delta S$ gives $\delta S=\delta A/4$ for nonextremal holes, and hence the standard Bekenstein-Hawking entropy. Restoring physical constants,Here restored $\kappa$ is the physical acceleration surface gravity and ${\ell }_P^2=G\hslash /c^3$ is the squared Planck length. The first-law comparison determines the area coefficient, leaving an additive entropy constant unspecified; the displayed formula is the usual convention. Its area scaling, rather than ordinary volume scaling, signals that the available thermodynamic degrees of freedom of a gravitating system are constrained by the horizon geometry.

Particle production explains the quantum input. For a real scalar field obeying the massless covariant wave equation, the Klein-Gordon inner product on solutions isThe wave equation makes its current divergence-free, so the product is independent of a Cauchy hypersurface when the boundary flux vanishes. Use suitably normalized wave packets or the appropriate continuum distributions. Asymptotic past and future Minkowski spacetimes give preferred positive-frequency solution mode bases $u_j^{\mathrm{in}}$ and $u_i^{\mathrm{out}}$, normalized by $(u_i,u_j)={\delta }_{ij}$, $(u_i^{\ast },u_j^{\ast })=-{\delta }_{ij}$ and $(u_i,u_j^{\ast })=0$. Between those regions, a nonstationary geometry generally has no preferred positive-frequency splitting.

The mode bases are related by a Bogoliubov transformation,The last identities follow from conservation of the Klein-Gordon inner product and are the canonical identities for a bosonic Bogoliubov transformation. Expanding the field in either complete basis and extracting its positive-frequency coefficient givesThus the state annihilated by every $a_j^{\mathrm{in}}$ hasThis particle number from Bogoliubov coefficients is nonzero when time evolution mixes positive and negative frequencies. It describes one state as vacuum in the early particle basis and populated in the late basis; it does not require an arbitrary choice of a vacuum at each intermediate time. For a finite implementable transformation the state remains a squeezed pure quantum state. In infinitely many modes, a common unitary bosonic Fock space implementation requires additional conditions, such as the beta map being a Hilbert-Schmidt operator; finite wave packet observables need not share every global divergence of an ideal continuum calculation.

For gravitational collapse, the future is not globally Minkowski: it contains a black hole. The relevant late out-basis includes modes reaching future null infinity together with modes entering the future event horizon. Modes on infinity alone are not a complete Cauchy basis. Nevertheless the same mode-mixing calculation determines the outgoing particle flux. Near a nonextremal horizon, the logarithm in the tortoise coordinate converts regular early null coordinates into the Hawking exponential ray mapwhere $U$ is the early affine null coordinate and $u$ is late retarded time at infinity. Backward propagation of a late outgoing mode $e^{-i\omega u}$ gives a factor $(U_H-U{)}^{i\omega /\kappa }$ for $U<U_H$. Its positive- and negative-frequency Fourier integrals differ by analytic continuation around the logarithmic branch. With a convergence regulator they reduce toTaking $ϵ\downarrow 0$ gives the thermal ratio of Hawking Bogoliubov coefficients, $|\beta {|}^2/|\alpha {|}^2=e^{-2\pi \omega /\kappa }$. Combining this with the canonical normalization yields the bosonic occupation $[e^{2\pi \omega /\kappa }-1{]}^{-1}$. Late-time wave packets turn formal continuum coefficients into a finite number flux.

A Schwarzschild black hole has $\kappa =1/(4M)$, so this spectrum has $T_H=1/(8\pi M)$. The exterior curvature potential partly reflects the outgoing modes. Its transmission probabilities are the greybody factors, giving, for one massless scalar species in the stationary late-time approximation,Thus the horizon temperature is universal, while the spectrum received at infinity is not a perfect featureless blackbody radiation spectrum. The derivation assumes the near-horizon quantum state inherited from a regular collapse vacuum and ignores rapid backreaction over the timescale of the packets. Exponentially blueshifted precursor frequencies indicate the usual short-distance assumption in the semiclassical calculation, not an independently demonstrated quantum-gravity description of the endpoint.

The outgoing positive energy flux drives black-hole evaporation. At leading order for a large isolated Schwarzschild black hole, area scales as $M^2$ and temperature as $M^{-1}$, giving $dM/du\sim -\gamma /M^2$ and a lifetime of order $M^3$, with the coefficient depending on species and transmission factors. Its negative heat capacity of a Schwarzschild black hole, $dM/d{T}_H=-8\pi M^2$, means that it heats up as it loses mass and cannot be in stable canonical equilibrium with an unlimited thermal reservoir. The approximation fails when quantum-gravitational scales are reached and does not fix whether the endpoint is complete evaporation, a remnant, or something else.

Area loss does not contradict the classical area theorem, because the quantum stress-energy tensor need not obey the classical null energy condition; negative horizon energy accompanies the positive outgoing flux. The appropriate thermodynamic statement is the generalized second law, involving $S_{\mathrm{BH}}+S_{\mathrm{outside}}$. Finally, Hawking radiation raises the black hole information paradox: if a pure quantum state completely evaporates into an exactly thermal mixed quantum state and its partners disappear, the result conflicts with ordinary unitary time evolution. Early outgoing radiation can be mixed simply because of its entanglement with interior modes while the complete state remains a pure quantum state; that fact alone is not information loss. Resolving the fate of all correlations through the evaporation endpoint requires more than the leading semiclassical flux calculation.

Assume a single-phase ideal gas mixture in a closed parcel, with common temperature $T$, total pressure $P$, and molar species amounts $N_i$. Define $N={\sum }_i{N}_i$ and volume mixing ratios $x_i=N_i/N$, so ${\sum }_i{x}_i=1$. Ideal mixing makes partial pressure $P_i=x_{i}P$.

For a pure ideal species, $d{\mu }_i=-{\bar{S}}_{i}dT+{\bar{V}}_{i}d{P}_i$, with molar volume ${\bar{V}}_i=\mathcal{R}T/P_i$. At fixed temperature, integration from standard pressure $P^{\circ }$ givesThe same logarithmic composition term follows from the ideal mixing entropy $\Delta S_{\mathrm{mix}}=-\mathcal{R}{\sum }_i{N}_i\log x_i$. Extensivity gives $G={\sum }_i{N}_i{\mu }_i$. Hence the Gibbs free energy of an ideal-gas mixture isEquivalently, write each standard chemical term as ${\mu }_i^{\circ }=h_i^{\circ }-T{s}_i^{\circ }$, using molar enthalpy and entropy. For zero abundance the limit $x\log x\to 0$ defines the extensive expression continuously.

The fractions alone do not specify the extensive Gibbs free energy: the total amount $N$ or an equivalent elemental inventory must also be fixed. Moreover, in a chemically closed reacting system $N$ need not remain constant; for example, hydrogenation of carbon monoxide changes the total molecular count while conserving atoms. Minimization is therefore over species amounts subject to elemental constraints, rather than over fractions with an artificially fixed molecular count. Nonideal gases require activities or fugacities; condensates contribute separate phase terms.

Assume the planet and present-day Jupiter have the same mass, comparable composition and initial cooling normalization, and both can be described by the supplied power-law planetary cooling over the relevant ages. Take the age of present-day Jupiter to be $t_J\simeq 4.5\mathrm{Gyr}$ and the young planet's age to be $0.050\mathrm{Gyr}$. The normalization cancels:The 5-AU orbit is used to regard stellar irradiation as modest compared with a hot-Jupiter orbit; the ratio is specifically intrinsic cooling, not the sum of intrinsic and reradiated luminosity. Initial entropy and irradiation can change the assumed normalization, so this estimate is not independent of formation conditions.

Such a young self-luminous giant is a favourable target for near-infrared exoplanet direct imaging. The increased intrinsic luminosity improves its contrast with the host star, and a wide physical orbit is more readily separated on the sky than a hot-Jupiter orbit. Two essential observing considerations are:

Thus direct imaging of a nearby young system, with sufficient angular resolution and contrast, is the natural discovery method, subject to the actual stellar brightness and instrumental limits.

Hot-Jupiter radius inflation is the excess radius of many irradiated gas giants relative to standard cooling models of the same mass, composition and age. A young gas giant starts with high entropy and a large radius, but ordinarily contracts through Kelvin-Helmholtz contraction. An old inflated object must either retain that thermal reservoir unusually well or receive power that affects the deep interior.

The first class is delayed cooling of an inflated giant planet. Two specific mechanisms are increased atmospheric opacity, which makes internal radiation escape less readily, and inhibited interior convection caused by composition gradients, potentially with layered transport. Both slow entropy loss. Stellar irradiation can also maintain an outer radiative blanket, but an appropriate irradiated boundary is already included in many baseline models and does not alone explain every extreme radius.

The second class is heating of an inflated giant planet. Two mechanisms are tidal heating from eccentricity, obliquity or other time-dependent tidal forcing, and Joule heating from wind-driven electric currents in an ionized atmosphere coupled to the magnetic field. Mechanical energy carried inward from atmospheric circulation is another candidate. Heating must be deposited sufficiently deeply and with sufficient power; superficial absorption that is promptly reradiated is not equivalent to deep interior heating. A circular synchronized isolated orbit does not automatically provide a persistent tidal source.

The classes can coexist and are tested through radius-age-irradiation trends, orbital properties and the required energy budget. Merely changing an observed transit altitude by a few atmospheric scale heights is generally distinct from inflating the bulk giant-planet interior.

At fixed mass, three broad controls of a super-Earth-sized body's observed radius are:

Bulk composition, envelope fraction/composition, and thermal/irradiation history are three independent controls. If mass is not fixed, the mass itself is an additional major variable. The label super-Earth does not ensure a rocky composition or an Earth-like atmosphere, and a radius alone does not determine which of these effects dominates.

For a quasistatic spherical star with fixed total mass, the stellar structure equations in radius areHere $s$ is the specific entropy of a moving mass shell, and the entropy term accounts for release or storage of thermal and gravitational energy. In thermal equilibrium it vanishes. An equation of state, opacity, energy-generation law and composition evolution close the equations. The absence of mass loss fixes the outer mass coordinate; it does not require the star to be in thermal equilibrium.

The radiative temperature gradient needed to carry the full luminosity isHere the actual temperature gradient is $\nabla =d\log T/d\log P$. For uniform composition the Schwarzschild criterion says that a radiative layer is stable if ${\nabla }_{\mathrm{rad}}<{\nabla }_{\mathrm{ad}}$, with marginal stability at equality, where the adiabatic temperature gradient is $(\partial \log T/\partial \log P{)}_s$. Use $\nabla ={\nabla }_{\mathrm{rad}}$ there. If ${\nabla }_{\mathrm{rad}}>{\nabla }_{\mathrm{ad}}$, convection carries some of the flux. Efficient convection gives $\nabla \simeq {\nabla }_{\mathrm{ad}}$; inefficient surface convection needs a transport prescription and can be superadiabatic. For composition gradients, use the Ledoux criterion, with threshold ${\nabla }_{\mathrm{ad}}+\varphi {\nabla }_{{\mu }_{\mathrm{mol}}}/{\delta }_{\mathrm{th}}$ instead. Here ${\nabla }_{{\mu }_{\mathrm{mol}}}=d\log {\mu }_{\mathrm{mol}}/d\log P$, $\varphi =(\partial \log \rho /\partial \log {\mu }_{\mathrm{mol}}{)}_{P,T}$ and ${\delta }_{\mathrm{th}}=-(\partial \log \rho /\partial \log T{)}_{P,{\mu }_{\mathrm{mol}}}$, using the mean molecular weight ${\mu }_{\mathrm{mol}}$.

Above the thin burning shell, take constant luminosity $L$, constant mean molecular weight ${\mu }_{\mathrm{mol}}$, ideal-gas pressure $P=\mathcal{R}T\rho$ with $\mathcal{R}=k_B/({\mu }_{\mathrm{mol}}{m}_u)$, and the Kramers opacity law $\kappa ={\kappa }_0\rho T^{-7/2}$. Then the temperature equation isTo implement the printed assumption that $P,\rho$ and $M_r$ are all power laws while retaining mass conservation, write $M_r\propto r^m$, $\rho \propto r^b$, $P\propto r^p$ and $T\propto r^d$. The mass equation gives $m=b+3$; the ideal gas law gives $d=p-b$; hydrostatic equilibrium gives $d=m-1$. The diffusion equation gives $d-1=2b-2-13d/2$. Solving this linear system gives the self-gravitating Kramers power-law envelope:The coefficients can also be matched consistently. If $r_c$ is the core boundary and $M_c$ its enclosed mass,Radiative transfer fixes the constant luminosity throughwhose right-hand side is independent of radius for these exponents. Thus this is a solution of all four envelope equations, not just a dimensional estimate.

Matching the boundary temperature to the isothermal core and extrapolating the idealized power law to the specified photosphere givesFor reference $M_R/M_c=(R/r_c{)}^{1/11}=2500^{1/10}\simeq 2.19$, so the envelope mass is not negligible in this particular solution. The given $M_c\simeq M_{\odot }$ fixes the core-radius normalization $r_c=(11/42)G{M}_c/(\mathcal{R}10^7\mathrm{K})$; a numerical absolute radius needs the unspecified molecular weight. The exponent ratio also gives $\nabla =d/p=5/21<2/5$, so the fully ionized monatomic version of this idealization is convectively stable.

A different common approximation neglects the envelope's self-gravity and sets $M_r\simeq M_c$ in the force equation. It cannot obey the exact mass equation with a nonzero density and exactly constant $M_r$. Under that additional approximation, the same calculation instead givesThis is the constant-core-mass power-law opacity radiative envelope limit, not the full power-law solution with changing enclosed mass. Stating which approximation is used resolves the otherwise different numerical radius ratios.

A compact nearly isothermal core, a thin luminosity-producing shell and a much larger cool envelope describe a shell-burning red giant. The full power-law profile is a deliberately simple model. The Kramers opacity law and an ionized ideal gas are not reliable at a $4000\mathrm{K}$ photosphere; partial ionization, other opacity sources and a convective envelope usually modify real cool giants. The quoted radius ratios are extrapolations within the stipulated model.

The microcanonical ensemble assigns equal probability to accessible quantum states in a narrow energy shell of an isolated system with fixed energy, particle number and external parameters. To derive the canonical ensemble, put a small system in weak thermal contact with a much larger isolated bath. If the small-system state has energy $E_i$, its probability is proportional to the number ${\Omega }_b(E_{\mathrm{total}}-E_i)$ of bath states. Expanding the bath entropy gives $S_b(E_{\mathrm{total}}-E_i)=S_b(E_{\mathrm{total}})-E_i/T+\cdots$, since $\partial S_b/\partial E=1/T$. Negligible interaction energy and an effectively constant bath temperature therefore give $p_i=e^{-E_i/(k_{B}T)}/Z$. The microcanonical ensemble suits an isolated fixed-energy system; the canonical ensemble suits a system exchanging energy with a thermostat at fixed temperature, with particle number still fixed.

For one atom, summing the finite geometric series over $m=-J,\dots ,J$ gives the partition functionAt $x=0$ the ratio is understood by continuity as $2J+1$. Logarithmic differentiation gives the mean magnetic dipole momentThis is odd in $x$, has small-$x$ expansion ${\mu }_{B}J(J+1)x/3+O(x^3)$, and saturates at $\pm J{\mu }_B$. The requested common graph is plotted in units of ${\mu }_B$: its ordinate is $⟨{\mu }_z⟩/(J{\mu }_B)$, so multiplication by ${\mu }_B$ gives $⟨{\mu }_z⟩/J$. The normalized slopes at zero are $(J+1)/3$ and the limits are $\pm 1$.

The magnetic susceptibility isEquivalently it is $N{\mu }_B^2\mathrm{Var}(m)/(k_{B}T)\ge 0$. For $k_{B}T\gg {\mu }_B|B|$, the apparent singular terms cancel, yielding Curie lawFor fixed nonzero $B$, $T\to 0$ and integer $J\ge 1$, the first excited magnetic level is separated from the aligned ground state by ${\mu }_B|B|$. Its exponentially small occupation givesThe moments are saturated, so changing the magnetic field hardly changes their alignment. This is a fixed-nonzero-field low-temperature limit. At exactly $B=0$, differentiating first gives the Curie expression for every $T>0$, which diverges as $T\to 0$; the two limits do not commute in this noninteracting model. For $J=0$ both moment and susceptibility vanish identically.

For noninteracting fermions, each single-particle level has occupation $n_i\in \{0,1\}$ by the Pauli exclusion principle. A microstate has $E={\sum }_i{ϵ}_i{n}_i$, $N={\sum }_i{n}_i$, so its grand-canonical weight is the product of its independent level weights. The grand potential is $\Omega =-kT\log \mathcal{Z}$; the particle number, entropy and pressure follow from $N=-{\partial }_{\mu }\Omega$, $S=-{\partial }_T\Omega$ and $P=-{\partial }_V\Omega$. This occupation description is the basis for the product, continuum and adiabatic calculations.

Summing separately over the two allowed occupations of every level givesThusThe factorization uses the absence of interactions; each level, including any degeneracy labels, is counted separately.

At thermal equilibrium, microscopic reversibility equates the probabilities of a path and its reversed path when both include their Boltzmann distribution initial weights. Denote their endpoint states by $z_1,z_2$, including velocity if needed. Since the Hamiltonian is even under time reversal in classical mechanics,This detailed balance condition and the energy balance for an autonomous Lagrangian yieldThis is the fluctuation-dissipation relation for a Langevin particle: the strength of Gaussian white noise is fixed by the damping and temperature, with $k_B$ the Boltzmann constant.

For an equilibrium coarse-grained variable, the unresolved microscopic states contribute entropy; their statistical weight is encoded in the Helmholtz free energy, rather than in a single microscopic energy. Relative to the same reference measure, ${\rho }_{\mathrm{eq}}(x)\propto e^{-\beta F(x)}$, so microscopic reversibility becomesThis extension assumes an equilibrium coarse-grained description with reversible path statistics; externally driven dynamics need not obey this relation.

Write $m_i={\mu }_i^2$ and $r^2={\varphi }_1^2+{\varphi }_2^2$. The uniform Landau free energy isso stationarity requiresFor $m_1,m_2>0$, the unique ground state is $(0,0)$ and the ${\mathbb{Z}}_2\times {\mathbb{Z}}_2$ discrete symmetry is unbroken. If $m_1<0$ and $|m_1|>|m_2|$, thenwhich breaks the first ${\mathbb{Z}}_2$ and leaves the second intact. This includes the case in which both masses are negative, because $m_1$ is then the more negative one. The remaining possible ordered case under the stated inequality is $m_1>0>m_2$, for whichand only the second ${\mathbb{Z}}_2$ is broken. The Hessian matrix in each ordered state is positive because the uncondensed direction has squared mass $m_k-m_j>0$, where $m_j$ is the condensed, more negative mass.

Near either continuous transition the nonzero order parameter is proportional to $(-m_i{)}^{1/2}$. Hence the mean-field critical exponents are

The lower critical dimension is the dimension at or below which fluctuations destroy the proposed finite-temperature ordered phase. Here the broken symmetry is discrete, so $d_{\mathrm{l}}=1$. In one dimension a domain wall interpolating between the two signs has finite energy, whereas its possible position gives an entropy growing as $\log L$. Domain walls therefore occur with nonzero density at every positive temperature and split the system into domains of finite typical length. Thus there is no finite-temperature ordered phase in one dimension.

Let $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ be the material derivative, and write $h$ for specific enthalpy. The ideal magnetohydrodynamic induction equation and Gauss's law for magnetism giveDotting the momentum equation with $\mathbf{B}$ eliminates the Lorentz force density. Therefore the cross-helicity density satisfiesThe first law of thermodynamics gives $dh=Tds+dp/\rho$ for specific entropy $s$, without requiring uniform entropy. ConsequentlySince $\nabla \cdot \mathbf{B}=0$, the first term on the right is a divergence. This proves the cross-helicity conservation law, with fluxFor the perfect gas used here, $h=\gamma p/[(\gamma -1)\rho ]$.

The calculation in part (a) leaves the cross-helicity conservation law sourceAt positive temperature, this vanishes exactly when the magnetic field is tangent to constant-specific entropy surfaces: $\mathbf{B}\cdot \nabla s=0$. A homentropic flow is a sufficient special case; uniform entropy throughout space is not necessary. Mere advection of specific entropy, $Ds/Dt=0$, constrains its variation along the velocity rather than along the magnetic field, so it does not by itself eliminate this source. To conserve total cross-helicity in a volume, the net boundary flux must also vanish.

The chemical potential is the change in internal energy when one particle is added at fixed entropy and volume:Equivalently, $\mu =(\partial F/\partial N{)}_{T,V}$ for the Helmholtz free energy.

Let a microstate $s$ have energy $E_s$ and particle number $N_s$. Maximizing the Gibbs entropy subject to normalization and fixed mean values of $E$ and $N$ is the maximum-entropy derivation of equilibrium ensembles. With Lagrange multipliers $\alpha$, $\beta$, and $-\beta \mu$, varyThe stationarity equation is $-\log p_s-1-\alpha -\beta (E_s-\mu N_s)=0$. Normalization therefore gives the grand canonical ensembleHere $\mathcal{Z}$ is the grand canonical partition function and $k_B$ is the Boltzmann constant.

For one fermionic quantum state of energy $\epsilon$, the Pauli exclusion principle permits occupation numbers only $n=0,1$. Its two grand-canonical weights are $1$ and $e^{-\beta (\epsilon -\mu )}$, so its mean occupation is the Fermi-Dirac distribution

For a free nonrelativistic particle, $\epsilon ={\hslash }^2{k}^2/(2m)$. In a region of area $A$, the number of wave-vector states in the annulus from $k$ to $k+dk$, including the two spin angular momentum states, isSince $kdk=(m/{\hslash }^2)d\epsilon$, the two-dimensional free-electron density of states is constant:

We now use units in which $k_B=1$, as in the question. At zero temperature, the Fermi-Dirac distribution is a step function, and henceAt positive temperature the same fixed particle number satisfies the exact relationThusThis is the low-temperature particle-number cancellation for constant density of states. Therefore, for $T\ll {\epsilon }_F$, the mean number of particles in the energy interval $[\epsilon ,\epsilon +d\epsilon ]$ is

Finally, compare the internal energy with its zero-temperature value. The thermally excited particles above ${\epsilon }_F$ and holes below ${\epsilon }_F$ have equal leading particle numbers, so their terms proportional to ${\epsilon }_F$ cancel. With $x=|\epsilon -{\epsilon }_F|/T$,where the integral follows from the fermion-to-boson thermal integral ratio. Differentiation gives the linear low-temperature heat capacity of a two-dimensional Fermi gas:The requested power law is therefore linear in $T$, with exponent one.

Let two systems exchange energy inside an otherwise isolated composite system, so $E_1+E_2=E$ is fixed. By the Second law of thermodynamics, equilibrium maximizes the total entropyAt an interior maximum,Thus $T_1=T_2$, which is entropy maximization under thermal contact and the condition for thermal equilibrium.

Stability requires this entropy maximum to be locally strict, so each ordinary subsystem has a concave entropy-energy relation. Sinceone has $\partial T/\partial E>0$, or equivalentlyThis is the entropy concavity and positive heat capacity criterion.

Now let $N_{\uparrow }=\alpha N$ and $N_{\downarrow }=(1-\alpha )N$. Because the spins are distinguishable, the number of microstates is the binomial coefficientThe Boltzmann entropy $S=k_B\log \Omega$ and the Stirling formula give, to leading order for large $N$,The total energy isThereforeSolving for the up-spin fraction gives the independent spin-one-half two-level system resultFor ordinary positive absolute temperature,with $\alpha \to 0$ as $T\to 0^{+}$ and $\alpha \to 1/2$ as $T\to +\infty$. Since the spectrum is bounded above, the population-inverted range $1/2<\alpha \le 1$ is also mathematically possible and corresponds to negative temperature; $\alpha =1/2$ is the infinite-temperature state.

Substitution givesand direct differentiation yieldsThus the energy increases with temperature throughout either finite-temperature branch, in particular throughout the required positive-temperature range.

An adiabatic process exchanges no heat with its surroundings, so $\delta Q=0$. For a reversible process it is also isentropic. The first law then givesUsing $p=N{k}_{B}T/V$ and the Mayer relation,Integration yieldsSubstituting $T=pV/(N{k}_B)$ gives the reversible ideal-gas adiabat

At fixed volume, a reversible infinitesimal heating supplies $\delta Q_{\mathrm{rev}}=C_{V}dT$. Since entropy satisfies $dS=\delta Q_{\mathrm{rev}}/T$,