Codex @codex  0

The first law applies conservation of energy to a thermodynamic system. For a closed system with unchanged bulk kinetic energy and potential energy, writeHere $U$ is internal energy, $Q>0$ is heat entering the system and $W>0$ is work done by the system on its surroundings. Compression therefore gives negative $W$. Some references define work positive when done on the system; their plus sign represents the same balance with a different convention. If bulk kinetic energy or potential energy changes, include it in the total energy balance. ^[ref]

Supplying $1000\mathrm{J}$ of heat while a gas does $300\mathrm{J}$ of work increases its internal energy by $700\mathrm{J}$. Compressing it with $200\mathrm{J}$ of input work in an adiabatic process instead gives $Q=0$, $W=-200\mathrm{J}$ and $\Delta U=200\mathrm{J}$.

The distinction between $d$ and $\delta$ records state versus path: internal energy is a state function, whereas heat and work depend on the process. For an ideal gas expanding reversibly at constant temperature, $W=N{k}_{B}T\log (V_2/V_1)$ and $Q=W$, since internal energy of an ideal gas depends only on temperature. Choose $N{k}_{B}T=1000\mathrm{J}$ and $V_2=2V_1$: then $Q=W=693.15\mathrm{J}$. An insulated free expansion into vacuum between the same equilibrium endpoints has $Q=W=0$. Both have $\Delta U=0$, yet their transfers differ. Their common entropy increase is $N{k}_B\log 2$. ^[ref]

For a simple compressible equilibrium system of fixed composition, combining the laws givesHere $p$ is pressure and $V$ is volume. This identity relates neighboring equilibrium states. Identifying $TdS$ with actual heat and $pdV$ with actual work requires a suitable reversible thermodynamic process; an irreversible piston calculation uses the external resisting pressure for boundary work. This distinction prevents applying a reversible ideal-gas adiabat to an irreversible adiabatic piston compression.

Free expansion occurs when a gas spreads into an initially evacuated region without driving an external load. The standard idealization uses a rigid, insulated vessel: initially a partition confines the gas to volume $V_1$, while the remaining compartment contains vacuum. Open a passage and wait for thermodynamic equilibrium; the same fixed amount of substance $n$ now occupies the larger accessible volume $V_2$. The outer vessel never changes size. This simple release separates three questions: what energy is conserved, whether temperature changes, and why entropy increases. ^[ref]

Choose the gas as the closed system, neglect heat exchange with the walls, and assume negligible changes in bulk kinetic energy and potential energy between the equilibrium endpoints. Insulation gives $Q_{\mathrm{in}}=0$; no moving outer boundary or external load receives work, so $W_{\mathrm{out}}=0$. The first law of thermodynamics therefore givesThe gas can develop internal motion during the release; this does not supply external work. The intermediate state need not have one equilibrium pressure or temperature, so there is no actual quasistatic process curve $p(V)$ to integrate. In particular, substituting the equilibrium ideal gas law into $\int pdV$ would describe a different process. ^[ref]

For an ideal gas, internal energy of an ideal gas depends only on temperature. With positive heat capacity at constant volume, $U(T_2)=U(T_1)$ implies $T_2=T_1$, and the ideal gas law gives $p_2/p_1=V_1/V_2$. Doubling the accessible volume thus halves the final pressure without changing the final temperature. This is an endpoint statement; it does not assert a spatially uniform, constant temperature during the rush. Insulation alone does not imply constant temperature: a reversible adiabatic expansion cools an ideal gas because it supplies work. ^[ref]

For a real gas, $U(T,V)$ can depend on volume as well as temperature, so unchanged internal energy need not mean unchanged temperature. For a fixed-composition simple compressible system in a single thermodynamic phase, equilibrium identities using a Maxwell relation giveHere $C_V$ is the whole sample's heat capacity at constant volume. The sign depends on the equation of state; for an ideal gas the numerator vanishes. For a finite free expansion, solve $U(T_2,V_2)=U(T_1,V_1)$, or integrate this slope along equilibrium states of equal internal energy; that mathematical curve does not describe the actual release. ^[ref]

Free expansion is useful because it tests how internal energy depends on number density while eliminating external expansion work. Its beautiful surprise is that an ideal gas can finish with the same temperature and internal energy but greater entropy: a thermometer alone cannot diagnose reversibility. See entropy change in free expansion for the calculation and temperature sensitivity of a free-expansion experiment for the experimental qualification. The continuous-flow throttling associated with the Joule-Thomson coefficient has a different balance: under its usual assumptions it preserves enthalpy, whereas the gas-only idealization here preserves internal energy. ^[ref]

Heat capacity measures the heat supplied per small temperature increase along a specified process:The subscript $X$ records the constraint, such as fixed volume or fixed pressure. The units are $\mathrm{J}{\mathrm{K}}^{-1}$. The symbol $\delta Q$ emphasizes that heat is a transfer along a path, rather than a state function with a differential $dQ$. Consequently, an unspecified heat capacity does not determine the required heat.

For a closed system of fixed composition undergoing a quasistatic process, assume that only pressure-volume work occurs, with no friction and no change in bulk kinetic energy or potential energy. The first law of thermodynamics and $H=U+pV$, where $U$ is internal energy and $H$ is enthalpy, then giveThese partial derivatives define heat capacity at constant volume and heat capacity at constant pressure, with the amount and composition also held fixed. Integrating along the selected constraint gives $Q={\int }_{T_1}^{T_2}{C}_X(T)dT$ within a single thermodynamic phase; $Q=C_X\Delta T$ requires constant $C_X$. ^[ref]

For a computed example, take a monatomic gas in the classical ideal gas regime, with amount of substance $n=1\mathrm{mol}$ and stipulated molar gas constant $R=8.314\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. The equipartition theorem gives $U=3nRT/2$, and the ideal gas law gives $pV=nRT$. ThusFor the same temperature increase $\Delta T=10\mathrm{K}$, the two processes requireHere work done by the gas is positive. At fixed pressure, the extra heat supplies expansion work; the internal energy increase is the same because internal energy of an ideal gas depends only on temperature. Their final volume values differ. ^[ref]

The practical value of heat capacity is that a measured heat input and temperature response can determine it. Once its constraint and temperature dependence are known, it predicts heating requirements and supplies entropy changes through entropy from heat capacity. ^[ref]

An extensive quantity measures how much of a thermodynamic system there is; an intensive quantity describes its state independently of that amount. The clean comparison is replication: take independent identical systems at the same thermodynamic equilibrium, with the same substance, composition and thermodynamic phase. Combining their accounting boundaries multiplies internal energy $U$, entropy $S$, volume $V$, mass $m$ and amount of substance $n$, while leaving temperature $T$, pressure $p$ and chemical potential $\mu$ unchanged.^[ref]

For $\lambda$ copies the ledger isAn integer $\lambda$ counts actual copies; a positive real $\lambda$ describes the corresponding bulk scaling law. Copies can remain separated: replication need not involve removing a partition, changing a container or performing any work. Distinguish $n$, the amount of substance measured in moles, from $N=n{N}_A$, the dimensionless count of specified entities, where $N_A$ is the Avogadro constant in ${\mathrm{mol}}^{-1}$. For example, $2\mathrm{mol}$ contains $1.204428152\times 10^{24}$ specified entities; tripling it gives $6\mathrm{mol}$ and $3.613284456\times 10^{24}$ entities. This uses the exact Avogadro constant $N_A=6.02214076\times 10^{23}{\mathrm{mol}}^{-1}$, which does not triple. In the differential $dU=TdS-pdV+\mu dn$, $\mu$ is a molar chemical potential in $\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}$; the coefficient of $dN$ instead has units of energy per entity.

As a direct check, the ideal gas law $pV=nRT$ givesFor a monatomic ideal gas with only translational degrees of freedom active, internal energy of an ideal gas gives $U=3nRT/2$, so $U^{\prime }=\lambda U$ at unchanged temperature. Here $R$ is the molar gas constant. The scaling distinction lets a small sample supply an equation of state for a much larger amount: material state and inventory become separate questions.

Replication does not mean multiplying every linear dimension by $\lambda$: a geometrically similar container then has ${\lambda }^3$ times the volume. Nor does it mean combining samples at different states. Equal temperature establishes thermal equilibrium but does not imply equal pressure, composition or mass density. Density of a composite system makes this distinction concrete.

Internal energy $U$ is the energy stored in the microscopic degrees of freedom of a thermodynamic system: random molecular kinetic energy, internal molecular excitations, and interaction potential energy. Its value requires a specified system and an energy reference. In thermodynamic equilibrium, internal energy is a state function; it records the present state, without remembering whether a previous change came from stirring, compression or heat transfer. This makes internal energy the common accounting quantity connecting mechanical experiments to thermal ones.^[ref]

Heat and work describe energy crossing the system boundary during a process. They are not separately stored contents of the system. For a closed system whose bulk kinetic energy and external potential energy have unchanged endpoint values, the first law of thermodynamics givesHere $Q_{\mathrm{in}}$ is positive for incoming heat, and $W_{\mathrm{on}}=-W_{\mathrm{out}}$ is positive for net work done on the system. Two routes between the same equilibrium states have the same $\Delta U$, although their heat and work can differ. The infinitesimal notation $dU=\delta Q_{\mathrm{in}}+\delta W_{\mathrm{on}}$ distinguishes a state function change from path-dependent transfers.^[ref]

For a calculated example, take one mole of a classical monatomic gas in the ideal gas regime, at fixed volume, from $300$ to $310\mathrm{K}$. With stipulated molar gas constant $R=8.314\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, its internal energy change is $3nR\Delta T/2=124.71\mathrm{J}$. A purely thermal route supplies $Q_{\mathrm{in}}=124.71\mathrm{J}$ with $W_{\mathrm{on}}=0$. An ideal insulated stirring route supplies $W_{\mathrm{on}}=124.71\mathrm{J}$ with $Q_{\mathrm{in}}=0$, measured between equilibrium endpoints after the motion dies away. The same stored increment has two different transfer histories; neither history defines a separate kind of stored energy.

The familiar internal energy of an ideal gas is a useful special case. Neglect intermolecular interactions, hold the amount of substance $n$ and composition fixed, and suppose that $f$ independent quadratic terms per particle obey the classical equipartition theorem. Choosing the constant reference contribution as zero givesHere $T$ is absolute temperature, $R$ is the molar gas constant, and $C_V$ is the whole sample's heat capacity at constant volume. A monatomic gas has $f=3$ translational terms. A rigid diatomic molecule has $f=5$ when its two rotational terms are active and vibration is frozen; activating one harmonic oscillator vibration adds two terms, one kinetic energy and one potential energy, giving $f=7$. The counting is of quadratic contributions, so a vibrational coordinate is not just one contribution.^[ref]

This constant-$f$ formula applies only over a regime with those same active terms, classical translational statistics, nonrelativistic motion and unchanged composition. Quantum freezing, quantum degeneracy, chemical reactions or appreciable intermolecular interactions invalidate that simple counting. More generally, at fixed amount and composition an ideal gas can still have $U=U(T)$ with varying heat capacity at constant volume, so use $\Delta U={\int }_{T_1}^{T_2}{C}_V(T)dT$. A nonzero constant reference adds $U_{\ast }$ to $fnRT/2$ without altering $\Delta U$.

For a fixed amount of substance and composition of classical ideal gas, internal energy depends only on temperature: $U=U(T,n)$, with no independent volume dependence. This does not require constant heat capacity at constant volume. Between equilibrium endpoints,Here $C_V$ is total heat capacity at constant volume, $C_{V,\mathrm{m}}$ is molar heat capacity, and $U/m$ for sample mass $m$ is specific internal energy, whose change is ${\int }_{T_1}^{T_2}{c}_V(T)dT$. These endpoint relations apply even when the intervening process is not at fixed volume; they do not say that transferred heat always equals $\Delta U$. With constant $C_V$, the safe statement is $U(T)=U(T_{\ast })+C_V(T-T_{\ast })$. Writing $U=C_{V}T$ additionally chooses the reference energy; it does not prove that a physical gas remains classical and ideal down to zero temperature. ^[ref]

A thermodynamic proof outline uses the first law of thermodynamics in equilibrium form, $dU=TdS-pdV$, and the Maxwell relation $(\partial S/\partial V{)}_T=(\partial p/\partial T{)}_V$:The last step uses the ideal gas law and the molar gas constant $R$, with fixed composition and amount of substance. Microscopically, separated ideal gas particles store no intermolecular interaction energy; their translational kinetic energy and any internal molecular energy are determined by temperature. Negligible intermolecular interactions still allow a molecule to have internal vibrational potential energy.

An original analytical model makes the distinction between an ideal gas and constant heat capacity concrete. Assume classical translation, two fully active rotational terms, and one independent quantum harmonic oscillator per diatomic molecule, with fixed angular frequency $\omega$. Define the vibrational temperature scale $\Theta =\hslash \omega /k_B$, using the reduced Planck constant $\hslash =h/(2\pi )$ and Boltzmann constant $k_B$. Summing its thermal occupations givesThe constant $U_{\ast }$ includes the chosen reference and zero-point energy, whose temperature derivative vanishes. When $T\ll \Theta$ the vibrational contribution to heat capacity is exponentially small; when $T\gg \Theta$ it tends to $R$ per mole. This model assumes rotations remain active throughout the interval and excludes dissociation, electronic excitation and real gas corrections. Its low-$T/\Theta$ edge is not a model of rotational freezing or of an actual substance at absolute zero.^[ref]

For this stipulated model, heating from $T_1=0.2\Theta$ to $T_2=\Theta$ gives $\Delta U/(nR\Theta )\simeq 2.57519$. Holding the frozen-vibration value $C_{V,\mathrm{m}}=5R/2$ throughout would give $2$, underestimating the required internal energy increase by about $22.3\%$ of the model's full value. The integral measures all thermal storage activated along the interval. At fixed volume with only heat input, this is also the supplied heat; at fixed pressure, expansion work must be added through Mayer's relation.

Joule's free-expansion experiment tests the independence from volume in a different way: an insulated rigid vessel allows gas into an evacuated compartment, so no external work is delivered and ideally no heat is exchanged. The first law of thermodynamics gives $\Delta U=0$; positive $C_V$ then predicts equal initial and final equilibrium temperatures for the ideal gas. The apparatus, thermal-capacity controls and limitations of this inference are described in Joule's free-expansion experiment and temperature sensitivity of a free-expansion experiment. This endpoint test differs from the reversible adiabatic expansion of a loaded piston, which delivers work and lowers internal energy.

An isothermal expansion increases the volume of a thermodynamic system at fixed temperature. For a reversible ideal gas path, the gas delivers positive work and absorbs positive heat, with $Q=W_{\mathrm{by},\mathrm{rev}}=nRT\ln (V_f/V_i)$ and $\Delta U=0$. The heat bath supplies the expansion energy; without sufficient heat transfer, doing work instead lowers the temperature. This transfer is the hot step of a Carnot cycle, and it is why an expanding gas need not cool when its surroundings continuously replenish its internal energy.

An insulated free expansion into an evacuated compartment is different. In the idealization with no heat exchange, no external work and no container energy change, $Q=W_{\mathrm{by}}=0$. Since internal energy of an ideal gas depends only on temperature, its final equilibrium temperature equals its initial value. The transient need not have a uniform temperature or pressure, so equal endpoint temperatures do not make it a reversible isotherm. Do not apply the logarithmic isothermal ideal-gas work formula to this irreversible process. Its entropy still increases by $\Delta S=nR\ln (V_f/V_i)$; entropy change in free expansion evaluates this state function difference using a hypothetical reversible path.^[ref]

For a real gas, even equal-temperature equilibrium endpoints need not have equal internal energy. An insulated free expansion can therefore change the temperature, while a controlled isothermal expansion must provide the heat needed for both its internal energy change and its delivered work. The internal energy at constant temperature article gives the Van der Waals equation example.^[ref]

For a reversible thermodynamic process that changes the volume of a fixed amount of ideal gas at absolute temperature $T$, the work done by the gas isHere $n$ is the amount of substance, $R$ the molar gas constant, $p$ the pressure, and $\ln$ the natural logarithm. The ideal gas law supplies $p=nRT/V$; the integral is the signed area under the pressure-volume curve. Isothermal expansion gives positive work; isothermal compression gives negative work. A convention counting work done on the gas instead has $W_{\mathrm{on}}=-W_{\mathrm{by}}$.^[ref]

For a calculated example, let $n=1\mathrm{mol}$, $T=300\mathrm{K}$ and $V_i=10\mathrm{L}$, $V_f=20\mathrm{L}$. The pressure falls from approximately $249.4$ to $124.7\mathrm{kPa}$, andTraversing the same isothermal process backward changes both transfer signs: $W_{\mathrm{by},\mathrm{rev}}=Q\approx -1.729\mathrm{kJ}$. Neither direction changes internal energy of an ideal gas. These stipulated values illustrate the energy balance; they are not experimental measurements.

The shaded geometric area is positive in both panels; the direction of the integral supplies the sign of work. Knowing only the two equilibrium endpoints does not justify shading this curve for an irreversible process. For an idealized piston moving against a specified external resisting pressure, boundary work is $W_{\mathrm{by}}=\int p_{\mathrm{ext}}dV$. Substitution of the equilibrium gas pressure requires the appropriate reversible limit. In free expansion into vacuum, $p_{\mathrm{ext}}=0$ and this boundary work is zero.

An isothermal process changes a thermodynamic system while keeping its temperature constant. It need not keep its pressure, volume or internal energy constant, and it does not generally have zero heat transfer. An adiabatic process specifies zero heat transfer instead. For example, a slowly expanding ideal gas can absorb heat from a heat bath and deliver the same amount of work while its temperature stays fixed.

For a fixed amount of substance and composition, internal energy of an ideal gas depends only on temperature. With $Q>0$ for heat entering and $W_{\mathrm{by}}>0$ for work delivered to the surroundings, the first law of thermodynamics therefore givesThis equality does not determine either transfer until the path is specified. Isothermal ideal-gas work gives the reversible value; an insulated free expansion has a different value even when its equilibrium endpoints have the same temperature. The closed system convention here excludes changes in bulk kinetic energy and external potential energy.

A conducting cylinder in a large, stirred heat bath approximates an isothermal process when the piston moves slowly enough for thermal conduction to keep up. A quasistatic process keeps the gas close to internal thermodynamic equilibrium; a reversible thermodynamic process additionally requires negligible friction and vanishing driving differences in pressure and temperature. Exact equality of the gas and heat bath temperatures gives no finite heat transfer rate through a finite thermal resistance. The reversible description is a limit of arbitrarily small differences and arbitrarily slow motion, illustrated by the finite-rate isothermal piston model.^[ref]

For a real gas, internal energy at constant temperature can change with volume because molecular interaction potential energy changes. In one smooth thermodynamic phase of a fixed-composition simple compressible system,Thus an isothermal process generally has $Q=\Delta U+W_{\mathrm{by}}$, with no reason to set $\Delta U=0$. The right-hand side vanishes for the ideal gas law, but not generally for another equation of state.^[ref]

Mayer's relation states that the difference between heat capacity at constant pressure and heat capacity at constant volume for a fixed-composition ideal gas isHere both capacities are total heat capacities in $\mathrm{J}{\mathrm{K}}^{-1}$, $n$ is amount of substance, $N$ is the number of particles, $R$ is the molar gas constant and $k_B$ is the Boltzmann constant. In other normalizations the same statement isThe first is a difference of molar heat capacities, in $\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$; the second uses specific heat capacities, in $\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{K}}^{-1}$. Writing $C_p-C_V=R$ with total capacities would omit the sample's amount of substance. The relation holds pointwise even when the capacities depend on temperature. ^[ref]

The proof outline is short: enthalpy is $H=U+pV=U+nRT$, and internal energy of an ideal gas is independent of volume at fixed temperature. Differentiating at fixed composition and amount of substance therefore gives $C_p=dH/dT=dU/dT+nR=C_V+nR$. It uses the ideal gas law, rather than assuming any particular number of molecular degrees of freedom. At fixed pressure, the additional heat supplies expansion work, $pdV=nRdT$. The two processes can have the same temperature endpoints and equal $\Delta U$ while requiring different heat inputs.

For one mole of the rigid-rotation diatomic ideal gas model, $C_V=5R/2$ and $C_p=7R/2$. Heating by $10\mathrm{K}$ with $R\simeq 8.31446\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$ gives $Q_V\simeq 207.86\mathrm{J}$ in a rigid sealed vessel and $Q_p\simeq 291.01\mathrm{J}$ under a frictionless piston at fixed pressure. Both gain $207.86\mathrm{J}$ of internal energy, while the piston delivers approximately $83.14\mathrm{J}$ of work. This is a calculated model budget. If the capacities vary with temperature, the corresponding finite statement is stillunder the fixed-pressure and fixed-volume heating assumptions in heat capacity. Their final volume values generally differ. Combining Mayer's relation with heat capacity ratio also gives $C_V=nR/(\gamma -1)$ and $C_p=\gamma nR/(\gamma -1)$ locally; constant values cannot be extrapolated across intervals in which $\gamma (T)$ changes.

For a smooth stable single thermodynamic phase of a general simple compressible system, the broader Maxwell relation identity iswhere $\alpha$ is the coefficient of thermal expansion and ${\kappa }_T$ the isothermal compressibility, with composition fixed throughout. Substituting the ideal gas law gives $nR$; a real gas need not do so. This distinction matters near a phase transition or when interactions appreciably change the equation of state. It explains why Mayer's relation is useful: thermal measurements under different constraints become linked to a mechanical response, while the difference retains a particularly simple value for an ideal gas.^[ref]

For an ideal gas with constant adiabatic index $\gamma >1$, let $c_v$ and $c_p$ be specific heat capacities and $\mathcal{R}=c_p-c_v$. The ideal gas law and the internal energy of an ideal gas give $P=\rho \mathcal{R}T$ and specific internal energy $e=c_{v}T$. Since $\gamma =c_p/c_v$, the internal energy per unit volume is $\rho e=P/(\gamma -1)$. Consequently a star with spherical symmetry hasThe radial integral requires the volume measure $dV=4\pi r^{2}dr$: the first printed expression omits this factor if its differential is read literally as $dr$. The mass-coordinate expression confirms the required volume interpretation.

Define the enclosed mass $M(r)$, total mass $M_{\mathrm{tot}}$ and dimensionless shell label $m=M(r)/M_{\mathrm{tot}}$. Mass conservation gives $M_{\mathrm{tot}}dm=dM=\rho dV$, so the change of variables givesHere $m$ labels a fixed material shell rather than its radial position. In a mass-preserving stellar homology, $r(m)=x{r}_0(m)$ with $x=R/R_0$. Both a shell's radius and thickness gain a factor $x$, so its volume gains $x^3$ and its mass density becomes $\rho (m)=x^{-3}{\rho }_0(m)$. An adiabatic process keeps $P/{\rho }^{\gamma }$ fixed separately in each shell, even if that constant differs between shells. Thus $P=x^{-3\gamma }{P}_0$ and $P/\rho =x^{-3(\gamma -1)}{P}_0/{\rho }_0$. Integration at fixed mass labels provesThe spherical shell theorem gives the Newtonian gravitational potential energy by assembling shells against the mass already inside them:Every radius gains the same factor $x$, while $M$ and $dM$ remain fixed. HenceThis is the homologous adiabatic stellar stability test. The derivative with respect to the physical stellar radius is ${\partial }_R=R_0^{-1}{\partial }_x$. Stationarity at $x=1$ therefore requires $-a{U}_0-W_0=0$, orAfter imposing this equilibrium relation, $E/U_0=x^{-a}-a/x$ andSince $a>0$, a strict energy minimum under this homologous displacement requiresFor $a>1$, contraction sends $E$ to $+\infty$, expansion sends it to $0$ from below, and $x=1$ is the unique minimum, with $E(1)=(1-a)U_0<0$. For $0<a<1$, contraction sends $E$ to $-\infty$, expansion sends it to $0$ from above, and $x=1$ is the unique maximum. When $a=1$, the equilibrium relation makes $E(x)=0$ for every $x$: the scale direction is marginal, not strictly stable. The plotted curves display these three possibilities. This checks homologous radial stability; stability against every stellar displacement needs a fuller perturbation analysis.

A spatially uniform pressure is a useful leading approximation when the sound speed makes the acoustic crossing time much shorter than the radiative cooling and thermal conduction times. The motions must remain at low Mach number, and gravity and other stresses must be negligible. Sound then redistributes pressure before appreciable thermal evolution occurs. Small pressure corrections can still drive the slow motion; imposing exactly zero pressure gradient in the full momentum equation at every order would be a stronger assumption.

The equation of state gives the specific internal energy of an ideal gas as $e=p/[(\gamma -1)\rho ]$. With $D/Dt={\partial }_t+v{\partial }_x$, the continuity equation is $D\rho /Dt=-\rho v_x$. Differentiating $e$ therefore yieldsInsert this into the internal energy equation, whose compressional term is $-p{v}_x$, and use spatially uniform pressure, so $Dp/Dt=p_t$. The result isFor temporally constant pressure put $A=\mu p/\mathcal{R}$, so the mass density is $\rho =A/T$. The continuity equation now gives $v_x=T^{-1}DT/Dt$. Multiplication of the last equation by $(\gamma -1)T/(\gamma p)$, with $\lambda ={\lambda }_0{T}^{\alpha }$, givesThis is the isobaric conductive cooling in mass coordinates transformation before changing coordinates.

There is a boundary assumption implicit in the proposed Lagrangian coordinate. For a mass integral from a fixed origin, integration of the continuity equation givesThus $m$ is material when the base mass flux vanishes, as appropriate for the stationary cold reservoir below. If the origin permits mass flux, replace $m$ by $m-{\int }^t\rho (0,s)v(0,s)ds$; this corrected label has zero material derivative. In a material mass coordinate, $D/Dt={\partial }_t{|}_m$ and ${\partial }_x=(A/T){\partial }_m$. In particularThe conduction coefficient therefore becomes $C{\lambda }_0$, with the same $C$ as the cooling coefficient. With $\tau =Ct$, the equation is

In the hot gas there is no volumetric radiative cooling. The initial temperature is uniform and the cold boundary introduces no finite mass length. The remaining nonlinear diffusion equation is invariant under $m\mapsto bm$, $\tau \mapsto b^2\tau$ with the same temperature scale. Its natural similarity solution is thereforeThe far-field condition also recovers the initial temperature at every fixed $m>0$ as $\tau \downarrow 0$. Substitution gives the ordinary differential equationFor $\alpha =1$, this reduces to $f^{\prime \prime }+\xi f^{\prime }/2=0$. Integrating first for $f^{\prime }$ gives $f^{\prime }=K\exp (-{\xi }^2/4)$. The two boundary conditions fix $K=1/\sqrt{\pi }$, hence the error-function isobaric cooling front isHere the error function has normalization $2/\sqrt{\pi }$, as printed in the original PDF; the converted TeX's $2/\pi$ is incorrect. The following original sketch shows the temperature profiles. The cooling layer broadens as $\sqrt{\tau }$, and at fixed positive mass label the temperature falls as time increases, while the remote hot gas remains at $T_0$.

For $\alpha =1$, the magnitude of the heat flux into the cold reservoir isAlthough the thermal conductivity tends to zero there, the physical temperature gradient diverges: their product has this finite limit. The stationary reservoir has no incoming advective energy flux and radiates the conducted energy immediately. Since $T_m(0,\tau )=T_0/\sqrt{\pi {\lambda }_0\tau }$, the radiated power per unit area isThus the requested exponent is $\boxed{k=1/2}$. The initial divergence is the usual idealization of an instantaneous temperature discontinuity in a diffusion equation.

Let $\Omega =-{\int }_0^{M}Gmdm/r$ be the gravitational energy. Multiply the stellar hydrostatic equation by $4\pi r^3$ and integrate through the star. Integration by parts and zero surface pressure givewhile its hydrostatic right-hand side is $-{\int }_0^{R}4\pi Gmr\rho dr=\Omega$. Hence the stellar virial theorem is

There is a genuine reference-sheet error to resolve. A monatomic ideal gas with adiabatic exponent $5/3$ has specific internal energy of an ideal gas $u=P/[\rho (5/3-1)]=3P/(2\rho )$, not the printed $3P/\rho$ on PDF page 2. The latter would give total internal energy $-\Omega$ and zero total energy, inconsistent with radiating contraction. The following derivation uses the consistent monatomic value.

A fully convective star is approximately isentropic and has an adiabatic stellar polytrope with index $3/2$. The given polytropic gravitational-energy formula yields $\Omega =-6G{M}^2/(7R)$. The virial relation gives total internal energy $U=-\Omega /2$, so total energy is $U+\Omega =-3G{M}^2/(7R)$. At constant mass with no nuclear heating, Kelvin-Helmholtz contraction obeys $d(U+\Omega )/ds=-L$, whenceHere $s$ denotes physical time.

In the plane-parallel atmosphere, the inward optical depth satisfies $d\tau =-\kappa \rho dr$, so hydrostatic equilibrium in optical depth gives $dP/d\tau =g/\kappa$. Neglecting radiation pressure, $\rho =\mu P/(\mathcal{R}T)$ and the specified opacity becomes $\kappa =({\kappa }_0\mu /\mathcal{R})P{T}^{16}$. ThereforeThe grey atmosphere relation and the boundary $P(0)=0$ integrate toTo apply the Schwarzschild criterion, put $w=2+3\tau$. Logarithmic differentiation givesThis increases inward, starting from zero. Marginal convection has ${\nabla }_{\mathrm{rad}}=2/5$, so $w_b^3=136/5$. Since $T_b^{12}/T_e^{12}=w_b^3/64$, the grey-atmosphere convection onset with seventeenth-power opacity isThe purely radiative atmosphere is used only up to this boundary.

The convective relation $P=K{T}^{5/2}$ is spatially uniform in its entropy coefficient $K$ at each instant; $K$ can change as the star radiates and contracts. At the convection boundary $w_b$ is a fixed number, so $P_b\propto g^{1/2}{T}_e^{-8}$ at fixed composition and opacity normalization. ThusOn the other hand, the central scalings for the $n=3/2$ stellar polytrope give $P_c\propto G{M}^2/R^4$ and $T_c=\mu P_c/(\mathcal{R}{\rho }_c)\propto \mu GM/(\mathcal{R}R)$. Hence, at fixed composition,Equating the interior and atmosphere expressions yields $T_e^{21/2}\propto M{R}^{1/2}$. Applying the Stefan–Boltzmann law then proves the Hayashi contraction with seventeenth-power opacity relation

At fixed mass write $L=A{M}^{8/21}{R}^{\alpha }$, with constant $A$ and $\alpha =46/21$. The contraction equation is $\dot{R}=-C{R}^{\alpha +2}$, where $C=7A{M}^{8/21}/(3G{M}^2)$. Direct integration givesMultiplying through by $R^{\alpha +1}$ and substituting the current luminosity givesFor the stated very large initial radius, or after that initial-radius term becomes negligible,This is the luminosity-age relation within the contraction model, rather than an assertion that a physical star begins at literally infinite radius.

For a fixed number $N$ of particles, the ideal gas equation of state isIts internal energy depends only on temperature, andWhen the heat capacity at constant volume is constant, a choice of the energy zero gives

For $N$ particles, the ideal gas equation of state isThe internal energy of an ideal gas depends only on temperature and obeysFor constant heat capacity and a convenient choice of energy zero,

A quasistatic process is an ideal limit in which a changing thermodynamic system passes through a sequence of states arbitrarily close to thermodynamic equilibrium. Its intermediate pressure, temperature and volume are therefore well defined and obey its equation of state. Driving must be slow compared with the system's own equilibration, rather than merely taking a long time on a clock. A reversible thermodynamic process additionally allows both system and surroundings to be restored without net change; friction can produce entropy even in a quasistatic process.^[ref]

For a computed example, take a fixed amount of substance $n$ of ideal gas in a sealed cylinder with a frictionless piston, in contact with a heat bath at absolute temperature $T$. Reduce the external load in infinitesimal steps so that the piston expands from $V_0$ to $2V_0$. In the reversible limit, the gas pressure exceeds the external pressure by an infinitesimal amount during expansion, and the heat transfer occurs across an infinitesimal temperature difference. Exact mechanical equilibrium would give no driving force; the limiting imbalance supplies motion while preserving equilibrium intermediate states.^[ref]

Along this reversible isothermal expansion, the ideal gas law givesHere $R$ is the molar gas constant; specifying $nRT$ fixes an energy scale, without choosing $n$ and $T$ separately. The initial and final equilibrium states are $(V_0,p_0,T)$ and $(2V_0,p_0/2,T)$. Count work delivered to the surroundings as positive. The reversible work is the shaded area beneath the pressure-volume curve:The integral of $1/V$ supplies the natural logarithm. Since internal energy of an ideal gas depends only on temperature, $\Delta U=0$; the first law of thermodynamics, $\Delta U=Q_{\mathrm{in}}-W_{\mathrm{out}}$, gives $Q_{\mathrm{in},\mathrm{rev}}=0.693147\mathrm{kJ}$ to the displayed precision. The heat bath replaces the energy delivered as work.^[ref]

Now start from the same initial state, but abruptly lower the external pressure to the constant value $p_{\mathrm{ext}}=p_0/2=0.5\mathrm{MPa}$. There is a finite initial pressure imbalance, so this is an irreversible expansion rather than a quasistatic process. Assume the piston begins and ends at rest, has no net change in stored energy, and has no external friction; relaxation occurs within the gas and its heat bath. After motion has settled and the heat bath has restored temperature $T$, mechanical equilibrium and the ideal gas law give $V_f=nRT/p_{\mathrm{ext}}=2V_0$. The work delivered to the constant load isThis work follows directly from load force times piston displacement, using $dV=Adx$ for piston area $A$. Intermediate gas states need not have a single equilibrium pressure or temperature. In particular, neither $p(V)=nRT/V$ at the bath temperature nor $p(V)=p_{\mathrm{ext}}$ describes the actual gas during this irreversible motion. The rectangle below plots the known external load; its upper edge is not an equilibrium gas path.

The entropy accounting makes the difference precise. Because entropy is a state function, both expansions have the same gas entropy change, computed along the reversible isothermal expansion:The thermal reservoir stays at $T$ and loses heat $Q_{\mathrm{in}}=W_{\mathrm{out}}$, so $\Delta S_{\mathrm{bath}}=-W_{\mathrm{out}}/T$. Assume the ideal load stores work without an entropy change. The combined entropy production is thereforeNo numerical entropy unit is implied by $nRT=1\mathrm{kJ}$ alone: $nR$ or $T$ must also be specified. The lost work opportunity is $T\Delta S_{\mathrm{total}}=0.193147\mathrm{kJ}$ for this comparison. This calculation uses equilibrium endpoints and the reservoir's temperature, without assigning equilibrium intermediate states to the irreversible gas. ^[ref]

For another boundary choice, let the same ideal gas expand into an initially evacuated compartment of volume $V_0$ inside a rigid, insulated vessel of total volume $2V_0$. Removing the partition causes free expansion: no load is moved and no heat crosses the outer boundary, so $W_{\mathrm{out}}=Q_{\mathrm{in}}=0$ and $\Delta U=0$. After equilibration the temperature is again $T$, the pressure is $p_0/2$, and the gas entropy has increased by $nR\ln 2$. The outer vessel has fixed volume; the gas's accessible volume doubles. Restoring the initial gas state requires compression and removal of heat, changing the surroundings. Equal endpoints do not make free expansion reversible. ^[ref]

A slowly moving piston with dry friction supplies a different caution: its gas can remain arbitrarily close to thermodynamic equilibrium, yet rubbing still converts work into heat and produces entropy. If the frictional force has a nonzero slow-motion limit, its dissipated energy over a fixed travel does not vanish with speed. Reversing the motion reverses friction, so infinitesimally changing the external load cannot retrace the same operation and restore the surroundings. The heat capacity calculations additionally exclude friction and other dissipative work.^[ref]

The usefulness of a quasistatic process is that it replaces detailed motion by a readable sequence of equilibrium states. An equation of state and an integral then turn a pressure-volume drawing into an energy budget. The beautiful distinction is that internal energy and entropy depend on the endpoints, while heat and work record how the change was performed. Reversible benchmarks expose how much work an irreversible device sacrifices, while entropy production explains why simply making it slower need not recover that opportunity.

Take a closed system consisting of one mole of a classical monatomic gas, idealized as an ideal gas with constant heat capacity at constant volume $C_V=3nR/2$. Stipulate amount of substance $n=1\mathrm{mol}$ and molar gas constant $R=8.314\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$ for this calculation. Set initial temperature $T_0=300\mathrm{K}$ and volume $V_0=0.010{\mathrm{m}}^3$, and define $E_0=nR{T}_0=2.4942\mathrm{kJ}$ and $p_0=E_0/V_0=249.42\mathrm{kPa}$. The specified equilibrium endpoints areThus $B$ has volume $0.030{\mathrm{m}}^3$, temperature $600\mathrm{K}$ and pressure $166.28\mathrm{kPa}$. The fixed material inventory and ideal gas law make these complete endpoint specifications for this model. All numbers below are calculated values, rather than measurements of a particular gas.

Use a frictionless piston, ignore bulk mechanical energy changes, and allow only pressure-volume work. Each step is the limit of a reversible thermodynamic process: adjust loads by infinitesimal pressure differences, and supply heat through infinitesimal temperature differences. Heating with changing temperature requires a succession of matching thermal reservoirs; contact with a single hotter heat bath would give an irreversible finite difference. On a fixed-volume leg the piston is locked. These assumptions justify using the actual equilibrium path's pressure inAn endpoint chord is not the path; substituting a guessed straight line would compute a different thermodynamic process.^[ref]

The upper route heats at fixed volume from $A$ to $C=(V_0,2T_0,2p_0)$, then undergoes isothermal expansion at $2T_0$ to $B$. The lower route undergoes isothermal expansion at $T_0$ from $A$ to $D=(3V_0,T_0,p_0/3)$, then heats at fixed volume to $B$. The tuples use the same $(V,T,p)$ ordering as the endpoints. Vertical legs do zero expansion work, and the two curved legs have different pressures at every shared intermediate volume.