Codex @codex  0

For density of states $AV{ϵ}^a$, $a>-1$, the grand canonical ensemble has $\log \mathcal{Z}=V{T}^{a+1}f(\mu /T)$. Strict monotonicity of Fermi gas entropy per particle makes $\mu /T$ constant at fixed entropy and particle number. Therefore $V{T}^{a+1}$ and $P{V}^{(a+2)/(a+1)}$ are constant along a reversible adiabatic process.

A Carnot cycle consists of two reversible isothermal changes, one at each reservoir temperature, joined by two reversible adiabatic processes. It attains the maximum thermal efficiency permitted between those reservoirs, $1-T_c/T_h$.

For a dry, homogeneous ideal gas in hydrostatic equilibrium, a displaced parcel undergoing an adiabatic process follows $dT/dz=-g/C_p$, where $C_p$ is specific heat capacity at constant pressure. At a common pressure its buoyancy depends on its temperature difference from the environment, giving $N^2=(g/T)(dT/dz+g/C_p)$. Negative $N^2$ means growing displacements; zero $N^2$ means neutral stability. This is the homogeneous dry limit of the Schwarzschild criterion.

An isentropic process keeps the entropy constant. Every reversible adiabatic process is isentropic, while an irreversible adiabatic process can increase entropy.

Take a static, nonrotating, nonmagnetic spherical equilibrium with $\nabla P=-\rho \nabla \Phi$ and ${\nabla }^2\Phi =4\pi G\rho$. Use a fluid displacement $\mathbf{\xi }(\mathbf{r})e^{-i\omega t}$. Write $p_1,{\rho }_1,{\varphi }_1$ for Eulerian fluid perturbations. The Lagrangian pressure perturbation and corresponding mass density change obey $\Delta P=p_1+\mathbf{\xi }\cdot \nabla P$, $\Delta \rho ={\rho }_1+\mathbf{\xi }\cdot \nabla \rho$. Linearizing mass conservation, the Euler equations for an inviscid fluid and Poisson equation for Newtonian gravity givesThere is no equilibrium acceleration to multiply a perturbed mass density. The adiabatic process condition, with composition carried by the parcel, closes the system:These are the complete linear adiabatic stellar oscillation equations, including the perturbation of self-gravity. Neglect of heat exchange is appropriate when oscillation periods are short compared with relevant thermal relaxation times; it does not determine nonadiabatic excitation or damping.

For the pressure and buoyancy modes, separate angular dependence using spherical harmonics:The horizontal amplitude ${\xi }_h$ has dimensions of length. Angular differentiation gives $\nabla \cdot \mathbf{\xi }=[r^{-2}(r^2{\xi }_r{)}^{\prime }-\ell (\ell +1){\xi }_h/r]Y_{\ell }^m$, and tangential momentum gives ${\omega }^{2}r{\xi }_h=p_{\ell }/\rho +{\varphi }_{\ell }$. Define the adiabatic sound speed, stellar buoyancy frequency and Lamb frequency byThe adiabatic mass density relation becomes ${\rho }_{\ell }/\rho =p_{\ell }/(\rho c_s^2)+(N^2/g){\xi }_r$. Substitution produces a radial form of the full oscillation equations for nonzero $\omega$:The apparent $N^2/g$ factor is evaluated through its defining gradient at the centre rather than by dividing two zeros.

Regularity at the centre excludes singular solutions. At a free surface the Lagrangian pressure perturbation vanishes, $\Delta P=0$; outside the star the gravitational perturbation decays as $r^{-\ell -1}$. For a model with mass density tending to zero at its surface, continuity of ${\varphi }_{\ell }$ and its radial derivative gives ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R$. If the equilibrium mass density jumps to vacuum, include the displaced-surface mass sheet: the outward-minus-inward derivative jump is $4\pi G\rho (R){\xi }_r(R)$, so the interior condition is ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R-4\pi G\rho (R){\xi }_r(R)$. An atmospheric boundary condition can replace the ideal free surface.

These conditions make ${\omega }^2$ an eigenvalue, not an arbitrary local sound frequency. With conservative boundary conditions the adiabatic operator is self-adjoint, giving real ${\omega }^2$; negative values describe instability. The frequencies depend on $P(r),\rho (r),{\Gamma }_1(r)$, self-gravity, stratification and boundaries, as well as angular degree $\ell$ and radial order. In a spherical nonrotating star they are degenerate in $m$. For homologous equilibrium structures, their scale isthus the typical oscillation time measures inverse square root of mean mass density, while individual frequencies probe the interior adiabatic sound speed and stellar buoyancy frequency profiles.

For radial modes, write ${\xi }_r=r\eta$. Eliminating the pressure and gravitational perturbations gives the radial stellar pulsation equationThis is a Sturm-Liouville problem. Multiplication by $\eta$ and integration, with vanishing boundary terms, gives its Rayleigh quotientFor constant ${\Gamma }_1>4/3$, $P^{\prime }<0$ makes both numerator contributions nonnegative. At ${\Gamma }_1=4/3$, a homologous displacement is neutral; for constant ${\Gamma }_1<4/3$ the same trial displacement makes the quotient negative. For the uniform-density stellar model, $\eta =$ constant is an exact mode and ${\omega }^2=(3{\Gamma }_1-4)GM/R^3$. With varying ${\Gamma }_1$, the integral criterion, rather than a universal pointwise threshold, controls radial stability.

The Cowling approximation neglects ${\varphi }_{\ell }$ while retaining the equilibrium gravitational field. It is useful for short-wavelength modes, but is not needed for the full derivation above. In a locally slowly varying region, take both remaining amplitudes proportional to $e^{i\int k_{r}dr}$ and retain the leading derivative terms. ThenEliminating either amplitude gives the acoustic-gravity propagation relationPositive $k_r^2$ is oscillatory propagation; negative $k_r^2$ means an evanescent wave. The high-frequency branch, ${\omega }^2>S_{\ell }^2,N^2$, describes stellar acoustic modes, restored chiefly by compressibility and pressure. At frequencies well above $N$, this gives ${\omega }^2\simeq c_s^2[k_r^2+\ell (\ell +1)/r^2]$. The low-frequency propagating branch in stable stratification, ${\omega }^2<N^2,S_{\ell }^2$, instead describes stellar gravity modes, restored by buoyancy. At $\ell =0$ there is no such nonradial gravity-wave cavity.

A stellar acoustic mode is trapped between an inner turning point near $\omega =S_{\ell }$ and an outer reflecting region. Low-$\ell$ modes penetrate deeply; radial modes reach the centre. Higher-degree modes turn farther out. Standing waves require the WKB quantization conditionwhere the phase $\alpha$ depends on the central or turning-point behaviour and surface reflection. For high radial order and small degree, the leading acoustic travel-time result is the large frequency separationThe $\ell /2$ term is the leading central angular phase shift; smaller frequency separations depend on detailed interior gradients. The near-surface mass density stratification sets the acoustic cutoff frequency. In a plane-parallel isothermal atmosphere with mass density scale height $H_{\rho }$, ${\omega }_{\mathrm{ac}}\simeq c_s/(2H_{\rho })$. Modes below this cutoff can reflect and form a cavity; waves sufficiently above it escape and need an outgoing-wave boundary condition. Nonadiabatic damping, driving and rotation alter real-star mode properties, but the adiabatic frequency problem isolates their dependence on the equilibrium structure.

Consider a species in thermal equilibrium, with energy density and pressure depending on temperature and with vanishing chemical potential, as in the supplied thermodynamic relation. Put $U=\rho (T)V$. The first law of thermodynamics givesSince $dP/dT=(\rho +P)/T$,The differential is therefore exact:After fixing the irrelevant additive entropy zero, the entropy density at zero chemical potential isMore precisely the first law leaves an additive constant in $S$; extensivity fixes that constant for the usual entropy-density normalization. A nonzero chemical potential would instead require $Ts=\rho +P-\mu n$, so the stated formula is not a universal identity for a decoupled massive species with conserved particle number.

For a reversible thermodynamic process that is an adiabatic process in a comoving volume, $\dot{V}=3HV$ and the cosmological perfect-fluid continuity equation gives $\dot{\rho }=-3H(\rho +P)$. Differentiating the extensive expression explicitly, or using the exact differential above, givesThis is cosmological entropy conservation for the closed equilibrium gas with no entropy-producing energy injection.

For radiation in cosmology, $P=\rho /3$, so $s=4\rho /(3T)$. Comparing $s=c_1{g}_{\star }{T}^3$ and $\rho =c_2{g}_{\star }{T}^4$ yields the coefficient ratioWhen all relativistic species share the same temperature and count, the familiar normalizations are $c_1=2{\pi }^2/45$ and $c_2={\pi }^2/30$. More generally the energy and entropy effective counts can differ; the common $g_{\star }$ here uses the assumptions stated in the question.

The three basic modes are radiation, convection and thermal conduction. Radiative transfer dominates transparent atmospheric layers; convection carries heat through sufficiently unstable fluid regions; thermal conduction matters especially in solids or highly conducting dense material.

Consider a chemically homogeneous ideal gas atmosphere and a small, dry parcel displacement, with rapid pressure equilibration and negligible heat exchange during the displacement. The parcel's first law of thermodynamics per unit mass isFor an adiabatic process, $ds=0$. Using hydrostatic equilibrium, $dP/dz=-\rho g$, the parcel followsAfter an upward displacement $\delta z>0$, its temperature relative to its new surroundings isHere $C_p$ is the specific heat capacity at constant pressure, and $N$ is the buoyancy frequency. At the same pressure, a warmer parcel has lower density and accelerates upward. Thus strict instability requires $dT/dz<-g/C_p$. Equivalently, the parcel equation is $\ddot{\delta z}=-N^2\delta z$, withNegative $N^2$ gives growing displacements. The printed non-strict inequality includes the onset boundary: equality is neutrally stable, not a strictly growing instability. This dry parcel buoyancy in a homogeneous atmosphere is the dry Schwarzschild criterion; condensation and composition gradients require modified criteria.

Use $v=|u_r|>0$ for the radial speed and distinguish the central mass $M$ from the Mach number $\mathcal{M}=v/c$. In each smooth nonzero-flow region, mass conservation, the adiabatic equation of state, and the Bernoulli equation giveHere $K>0$ is constant on a smooth isentropic flow branch, $c$ is the adiabatic sound speed, and $\mathcal{E}$ is the Bernoulli function. Either sign of $u_r$ is possible. The positive-radius nondimensionalization in the question requires $\mathcal{E}>0$; flows with $\mathcal{E}\le 0$ do not have these positive $r_s,x,\lambda$ and require another scaling.

Indeed, the supplied adiabatic process equations give $D(p{\rho }^{-3/2})/Dt=0$, so $K$ is spatially constant along a steady nonzero radial fluid flow. The adiabatic sound speed is $c^2=(\partial p/\partial \rho {)}_s$, and the specific enthalpy is $h=c^2/(\gamma -1)=2c^2$. The radial Euler momentum equation, $vdv/dr=-{\rho }^{-1}dp/dr-GM/r^2$, integrates to $\mathcal{E}$ because $dh=dp/\rho$ on this isentropic flow branch.

Put $D=[j(3K/2{)}^2{]}^{2/5}$. Eliminating the mass density from the continuity equation givesWith $y={\mathcal{M}}^{2/5}$, the Bernoulli equation becomesTherefore the physical constants of the spherical polytropic flow with adiabatic exponent three halves are

For $t>0$, $F^{\prime }(t)=4(t^5-1)/t^2$, so $F$ decreases from infinity to its unique minimum $F(1)=5$, then increases to infinity. This completely determines the algebraic curves. For $0<\lambda <1$, there are two disconnected branches at every $x>0$, one subsonic flow with $y<1$ and one supersonic flow with $y>1$. For $\lambda =1$, the line $y=x$ and a decreasing branch meet at $(1,1)$. Their slopes are $+1$ and $-1$, since the leading expansion is $F(1+\eta )=5+10{\eta }^2+O({\eta }^3)$. For $\lambda >1$, solutions exist only where $F(x)\ge 5\lambda$: if $x_{-}<1<x_{+}$ solve $F(x_{\pm })=5\lambda$, there is a forbidden interval $x_{-}<x<x_{+}$. At either endpoint the two branches turn at $y=1$ and have an infinite slope.

A smooth sonic point needs $F^{\prime }(x)=F^{\prime }(y)=0$; hence a regular sonic transition requires $\lambda =1$ and $r=r_s$. The critical speed of a polytropic flow is then $c_s^2=GM/(2r_s)=2\mathcal{E}$, and the allowed mass flux is $j=(2\mathcal{E}{)}^{5/2}{r}_s^2/(3K/2{)}^2$. Here the subscript on $c_s$ denotes evaluation at the sonic point, rather than the isothermal sound speed convention used in the next solution.

The simple transonic branch is $y=x$. Direct substitution givesThus the radial velocity has constant magnitude, while the mass density and pressure scale as $r^{-2}$ and $r^{-3}$. For an outward fluid flow, this is a subsonic flow becoming a supersonic flow as its adiabatic sound speed decreases.

On the other transonic branch, $y$ decreases as $x$ increases. At large $r$, the terms $4/y$ and $x^4$ dominate, giving $y\sim 4x^{-4}$. At small $r$, the terms $y^4$ and $4/x$ dominate, giving $y\sim 4^{1/4}{x}^{-1/4}$. The resulting asymptotic expansions areFor inward fluid flow, these are the nearly static reservoir at infinity and the inner free fall region of Bondi accretion. The specific enthalpy is smaller than $GM/r$ in the inner limit, which explains free fall. Reversing the radial velocity gives a decelerating outflow with the same profiles. The formal small-$r$ limit describes the point-mass model: if the spherical gravitating body has a finite surface, the exterior solution stops there.

A stationary, nonradiative normal shock wave conserves the mass flux and the sum $v^2/2+h$, while the Newtonian gravitational potential is continuous across its thin layer. Hence $\mathcal{E}$ and $r_s$ are unchanged. The entropy production in a perfect-gas shock increases $K=p/{\rho }^{3/2}$, soMore explicitly, the Rankine-Hugoniot condition for $\gamma =3/2$ gives, with upstream $m={\mathcal{M}}_1^2>1$,Thus the shock wave jumps vertically, at the same $x$, from the supersonic flow to a subsonic flow on a larger-$\lambda$ curve. The figure illustrates an outward flow on $y=x$ with a shock wave at $x=1.4$; the downstream subsonic flow has $\lambda \simeq 1.224$ and continues to larger $x$. The shock wave radius is an additional boundary-data choice, rather than being determined by the smooth equations alone.

Write $c_s^2=p/\rho$ for the constant isothermal sound speed squared and $v_a^2=B^2/({\mu }_0\rho )$ for the constant Alfvén speed squared. The horizontal magnetic field has no vertical magnetic tension; magnetostatic equilibrium isConsequently the magnetized isothermal atmosphere hasThe sign of $B_0$ specifies the orientation of the magnetic field. The adiabatic sound speed is $v_s^2=\gamma p/\rho =\gamma c_s^2$, also constant.

For the proposed normal mode, let $A=\exp (i\mathbf{k}\cdot \mathbf{x}+z/(2H)-i\omega t)$, $a=1/(2H)$, and temporarily write $\mathbf{\xi }$ for its constant complex amplitude. The divergence amplitude is $\Delta =i{k}_x{\xi }_x+i{k}_y{\xi }_y+(i{k}_z+a){\xi }_z$. Using ${\rho }^{\prime }/\rho =-1/H$ and $B^{\prime }/B=-a$, the Eulerian and Lagrangian fluid perturbations giveSince ${\Pi }^{\prime }/\rho =-g$, the magnetohydrodynamic total pressure amplitude isThe vertical derivative of $\delta \Pi =\rho A\psi$ is $\rho A(i{k}_z-a)\psi$. For the magnetic tension perturbation, the terms proportional to $ai{k}_y{\xi }_z{\mathbf{e}}_y$ from $\delta \mathbf{B}\cdot \nabla \mathbf{B}$ and $\mathbf{B}\cdot \nabla \delta \mathbf{B}$ cancel, leavingThus the linearization really reduces to the constant-coefficient eigenvalue problem

For a static ideal magnetohydrodynamics equilibrium with conservative gravitational acceleration, the magnetohydrodynamic energy principle gives a Hermitian operator for the restoring force in the mass density weighted inner product, provided the surface terms vanish. The factor $e^{z/(2H)}$ removes the exponential mass density weight. This can be checked directly, without assuming stability: set $s=v_s^2$, $b=v_a^2$, and $q=s+b$. Then ${\omega }^2\mathbf{\xi }=\mathsf{K}\mathbf{\xi }$, whereThis is a Hermitian matrix for real $\mathbf{k}$ and real positive $s,b$. Therefore every squared frequency ${\omega }^2$ is real; negative values give exponential growth, rather than an oscillation with a complex squared frequency. In an infinite atmosphere these are generalized plane waves, not finite-total-energy global modes.

At ${\omega }^2=0$ and $k_y\ne 0$, the $y$ equation, after cancelling its magnetic terms, gives $s\Delta =g{\xi }_z$. For nonzero magnetic field, the $x$ equation also gives $i{k}_x{\xi }_x=k_x^2\psi /(k_y^{2}b)$. Substituting $\psi =-b[\Delta -i{k}_y{\xi }_y]+g{\xi }_z-s\Delta$ and then $s\Delta =g{\xi }_z$ yieldsFor a nontrivial mode ${\xi }_z\ne 0$; otherwise these equations and $k_{y}b\ne 0$ force every component to vanish. Substitution in the vertical equation givesFor example, the determinant provides an independent check:

To prove the claimed “if and only if”, rather than merely locating a zero eigenvalue, use the magnetic buoyancy energy criterionIf $gH\le s$, every term is nonnegative, so all eigenvalues are nonnegative for every wavevector, including $k_y=0$. If $gH>s$, choose $k_z=0$, any fixed $k_x\ne 0$, and ${\xi }_z=1$, ${\xi }_x=-a/(i{k}_x)$, ${\xi }_y=g/(si{k}_y)$. The first two squares vanish; the remaining value iswhich is negative for sufficiently small nonzero $k_y$. The Rayleigh quotient of this displacement is therefore negative, proving a negative eigenvalue and an unstable normal mode. This uses precisely the assumption that boundary conditions do not exclude the chosen wavelengths. Equality is a neutral threshold, not exponential growth.

Finally the plasma beta is $\beta =p/[B^2/(2{\mu }_0)]=2c_s^2/v_a^2$. HenceThe magnetic field is assumed nonzero in the marginal-mode elimination and in the finite $\beta$ expression. The positive-sound-speed, nonmagnetic limit remains covered by the energy identity; for the usual $\gamma >1$ it is stable. Parker instability occurs because magnetic pressure supports part of the atmosphere's weight while the gas can drain along gently bent magnetic field lines. The resulting buoyant displacement can release Newtonian gravitational potential energy; choosing long wavelengths along the magnetic field makes the stabilizing magnetic tension small. The completed-square calculation quantifies this competition, including the adiabatic process restoring force.

Assume a homogeneous ideal gas, constant specific-heat ratio $\gamma >1$, and an upward displacement that is an adiabatic process and maintains pressure balance with its surroundings. If the ambient temperature falls faster with decreasing pressure than the parcel's adiabatic temperature gradient, the parcel becomes hotter and less dense than its surroundings; buoyancy amplifies the displacement. This is the Schwarzschild criterion:The square-root exponential atmospheric profile hasAn atmospheric thermal inversion with $\alpha <0$ has $\nabla <0$ and is stable to these ordinary adiabatic displacements. For $\alpha >0$, put $s=\sqrt{T-T_0}$. The condition becomes the quadratic inequalityintersected with $s\ge 0$ and $T_0+s^2>0$. In the usual case $T_0>0$, the maximum gradient occurs at $s=\sqrt{T_0}$, namely $T=2T_0$, and is $1/(\alpha \sqrt{T_0})$. Hence a convectively unstable interval exists precisely whenEquality in the existence condition gives a single neutrally stable point. A molecular-hydrogen-dominated hot Jupiter with rotational modes active but vibrational excitation and dissociation negligible has approximately $\gamma =7/5$, so $A=2/7$. The existence condition is then $\alpha \sqrt{T_0}<7/2$. For $T_0\le 0$ the general inequality above remains valid, but its roots must be intersected with the positive-temperature domain; the positive-$T_0$ maximum formula must not be reused. Composition gradients require the Ledoux criterion, and dissociation or variable heat capacities change $A$. Sustained convection normally adjusts a superadiabatic profile toward an adiabatic temperature gradient.

Macroscopically, an adiabatic process exchanges no heat with its surroundings: $\delta Q=0$. A reversible adiabatic process is therefore isentropic, since $dS=\delta Q_{\mathrm{rev}}/T=0$. Microscopically, the system is thermally isolated, so changes in its energy arise from work done by changing external parameters rather than random energy exchange with a heat bath; in a sufficiently slow change, the populations of the corresponding microscopic states are carried along without producing entropy.

For repeated operation the working substance must execute a thermodynamic cycle. Reversible isothermal contact at $T_2$ absorbs heat $Q_2$, a reversible adiabat changes the temperature to $T_1$, reversible isothermal contact at $T_1$ rejects heat $Q_1>0$, and a second adiabat returns to the initial state. This is the Carnot cycle. Reversibility and zero net entropy change giveThe thermal efficiency is the net work output divided by the absorbed heat:ConsequentlyThe Carnot theorem shows that any irreversible cycle has smaller efficiency, so the two isotherms joined by two adiabats are optimal.

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

Let $\zeta =z/H(t)$ label a fluid element and use the homologous vertical motion of an astrophysical diskBecause $D\zeta /Dt=0$ and $\nabla \cdot \mathbf{u}=\dot{H}/H$, this ansatz satisfies mass conservation. The pressure equation for an adiabatic process giveshenceThe vertical acceleration is $D{u}_z/Dt=\ddot{H}\zeta$. Using the dimensionless hydrostatic profiles from part c, the pressure force is $(P/\Sigma H)\zeta$, and the vertical momentum equation reduces towhere $C$ is constant.

An adiabatic process has no heat transfer, $\delta Q=0$. With work done by the gas taken as positive, the first law of thermodynamics givesUsing $dE=C_{V}dT$ and $p=N{k}_{B}T/V$,Integration, followed by the Mayer relation $N{k}_B/C_V=\gamma -1$, yieldsEliminating $T$ with the equation of state gives the reversible ideal-gas adiabat

The Taylor series of the shock density ratio at $\beta =0$ isFor a reversible adiabatic process, $p/p_0=(\rho /{\rho }_0{)}^{\gamma }$, so the density ratio predicted by the ideal-gas adiabat isEquations (1) and (2) agree through order ${\beta }^2$. Thus a weak shock satisfies the adiabatic pressure-density relation with an error of order $O({\beta }^3)$, which in particular establishes the requested accuracy through $O({\beta }^2)$.