The local conservation laws. The first equation is the radial Euler equation: the material acceleration is caused by the pressure gradient force per unit mass, . It neglects gravity and viscosity. The second equation is mass conservation in spherical coordinates,
The term appears because the area of a spherical surface grows as .
To obtain the third equation, let be the specific internal energy of an ideal gas with constant heat capacity ratio . The first law of thermodynamics for a smooth adiabatically moving element is
Differentiating and collecting terms gives
and hence
Thus the entropy parameter is conserved along each smooth pathline; for fixed composition, specific entropy differs from only by a constant. This is local isentropic flow behind the shock. A physical shock wave itself generates entropy and must be treated with jump conditions, rather than by applying this smooth equation across the discontinuity.
The front-radius scaling. In the self-similar blast wave approximation there is no scale set by ambient pressure, cooling, gravity, or the initial source radius. Write . Using and , dimensional analysis gives
Thus , , and
The dimensionless coefficient depends on and the resolved similarity profiles; dimensions alone do not fix it. Define the length scale , so and the front is at the fixed coordinate .
Kinetic energy. Substitute and into the volume integral. Since and ,
But , so the first of the spherical blast-wave similarity energy integrals is
The conclusion requires a fixed, integrable similarity profile, as in the finite-energy blast solution. Expansion increases the swept-up mass while reducing the characteristic squared speed, leaving this product constant.
Pressure and internal energy. The derivatives of the similarity variable are
Writing a prime for differentiation with respect to , the material acceleration becomes
The radial momentum equation therefore requires a pressure gradient on the scale . Integrating over radius gives pressure similarity in a spherical blast wave:
where
The prefactor has units , as required for pressure. Strictly, the momentum equation alone permits an additional pressure uniform in space and arbitrary in time. In the scale-free strong-blast solution, shock matching and the remaining thermodynamic equation fix that freedom with the same similarity scaling; any allowed constant of integration is included in . An independently imposed external pressure would introduce another scale and invalidate the simple dimensional ansatz.
For constant , the internal energy is consequently
If is the total excess explosion energy, its normalization is
The two constant fractions need not be equal. No explicit solution for is needed to establish their time independence.
The entropy form of dry potential temperature follows from its parcel definition. The potential temperature is the temperature a parcel would reach if brought reversibly and adiabatically to a reference pressure . For dry air of fixed composition, treated as an ideal gas with constant specific heat capacity,
Per unit mass, the first law of thermodynamics, together with for specific volume, gives . The reversible form of the Second law of thermodynamics then gives
Hence dry reversible adiabatic motion has . Pressure changes cause expansion cooling or compression warming, while the pressure factor in theta compensates exactly. Adiabaticity alone is insufficient if irreversible mixing or viscous dissipation produces entropy. Condensation or freezing supplies latent heating, so ordinary dry potential temperature need not be conserved in a moist adiabatically moving parcel.
Take a fixed amount of dry air, with no phase change, moving quasistatically so that its pressure matches the surrounding atmosphere. From the first law of thermodynamics and ,
For adiabatic motion . Treating the parcel as an ideal gas of fixed composition makes its enthalpy depend only on temperature, so for moles even though the pressure changes. Hence . The environment is in hydrostatic equilibrium, , and the parcel is taken to have the ambient density at the displacement point. Using gives
Here is molar and is mass-specific heat capacity. Positive points upward; the commonly quoted positive dry adiabatic lapse rate is . Approximate constancy of and heat capacity gives a constant gradient. Condensation, external heating, entrainment or composition changes invalidate this dry-parcel relation.
For nonrelativistic particles the energy density is dominated by rest energy, while their isotropic kinetic pressure is approximately . Therefore
The matter pressure is negligible in this approximation, including in the gravitational combination . For isotropic radiation in cosmology, implies one third of the energy density in each diagonal spatial stress:
Because the two fluids do not exchange energy or momentum, they are separately conserved cosmological fluids: for each, rather than merely conservation of their sum. In a comoving volume proportional to , the first law of thermodynamics gives
Equivalently this follows by projecting the individual stress-energy conservation equation along the common comoving four-velocity. Differentiating and dividing by yields
For constant equation-of-state parameter ,
The extra power for radiation in cosmology represents the cosmological redshift of photon energies in addition to number dilution. Differentiating the total Friedmann equation gives only the total continuity equation; the separate equations also use the stated absence of interactions between components.
For a perfect fluid in general relativity, let be rest-frame total energy density divided by , its isotropic pressure, and its four-velocity, normalized by . Then
The contracted Bianchi identity and the Einstein field equations give stress-energy conservation, , for a constant cosmological constant. Contract with and use the normalization to obtain
In comoving coordinates and the covariant divergence is , where a dot is a cosmic time derivative and the Hubble parameter. Therefore
For a fixed region in comoving coordinates, let be its physical volume, where is a fixed coordinate-volume integral, and let . Since ,
Thus the cosmological perfect-fluid continuity equation is exactly equivalent to the adiabatic first law of thermodynamics, , for every nonzero such region. Conversely, expanding that first law of thermodynamics and dividing by recovers the same equation. No heat or particle flux crosses the comoving boundary; interacting components separately need source terms, although the total still obeys this conservation equation.
There is a sign inconsistency in the supplied formula sheet worth isolating. Its Riemann curvature tensor and Ricci tensor conventions, together with its stated metric signature, give for with constant . Its vacuum Einstein field equations with then require , whereas its printed vacuum Friedmann equation gives . For example, , and satisfy the latter but leave the former's component equal to . Changing the geometric cosmological constant term to repairs this mismatch. The conservation derivation is unaffected; below take the geometric cosmological constant to be zero and include any vacuum energy as a fluid, so no ambiguous sign enters.
There are two first-order effects: an intrinsic temperature change from adiabatic compression, and the gravitational redshift of the emerging photons. Evaluate the excess Newtonian gravitational potential at the boundary of the spherical fluctuation, with the exterior reference potential subtracted:
In the weak-field approximation, gravitational time dilation gives . A local emitter's successive wave crests therefore arrive farther apart at a reference observer outside the potential well. Their frequency changes by . Since a redshifted black-body radiation spectrum has its temperature scaled by the same factor as its frequency, the gravitational contribution is
For the intrinsic effect, compare the local scale factor with the unperturbed value at the reference epoch. An equilibrium photon gas has energy density and pressure . For adiabatic compression, the first law of thermodynamics applied to a comoving volume gives
Thus is constant along this adiabatic process, and the initial intrinsic contribution is . The smaller local expansion associated with the fluctuation leaves its radiation hotter; this is the compression or local expansion-history contribution, distinct from the clock-rate effect on the escaping radiation.
Multiplying the two first-order factors, and cancelling the common subsequent background redshift, gives
For the overdense region produced by adiabatic compression, , so the final term is positive. It competes with the negative gravitational redshift term. Their sum is the intrinsic-plus-endpoint part of the Sachs-Wolfe combination; computing a full observed Cosmic microwave background anisotropy can also require Doppler CMB anisotropy and the Integrated Sachs-Wolfe effect.[ref]
A stationary black hole behaves mechanically like an equilibrium thermodynamic system. Its horizon area A, asymptotic energy E, angular momentum J and electric charge Q are linked by the laws of black-hole mechanics. Let generate its Killing horizon, with the time translation normalized at infinity; its surface gravity is defined by on that horizon. We use and keep G explicit. The horizon electrostatic potential is measured relative to infinity.
The Zeroth law of black-hole mechanics states that kappa is constant on each connected stationary horizon, under the Einstein equations and suitable energy/regularity assumptions. It parallels the zeroth law of thermodynamics, which makes temperature uniform in equilibrium. Kappa need not be constant in a dynamical geometry. For an asymptotically flat static hole, local observers see gravitationally redshifted temperatures; the equilibrium quantity comparable to kappa is the temperature referenced to infinity.
The First law of black-hole mechanics relates nearby stationary solutions:
The same leading relation describes a small physical process which perturbs a horizon and subsequently settles to equilibrium, with the appropriate hypotheses on fluxes and final state. The work terms correspond to rotational and electromagnetic work. The structural analogy with the first law of thermodynamics, , identifies area as an entropy-like variable and surface gravity as a temperature-like variable. Classically this fixes their product in the first law, rather than their separate normalizations.
The second law of black-hole mechanics, or black-hole area theorem, states that future event-horizon area does not decrease classically. Its essential geometric ingredient is the Raychaudhuri equation for affinely parametrized null generators. In four dimensions a hypersurface-orthogonal congruence obeys
The null energy condition and Einstein equations make the last term nonpositive. If the expansion became negative, the inequality would focus the generators to a conjugate point in finite affine parameter. Under future completeness and predictability assumptions, that would conflict with their remaining on an achronal event horizon. Thus , and the area element, satisfying , cannot decrease. For a merger the final area is at least the sum of the initial component areas. These assumptions matter: the theorem is a classical statement, not a prohibition of quantum evaporation.
The third law of black-hole mechanics is an unattainability statement: a nonextremal hole cannot be driven to zero surface gravity in a finite physical process satisfying the usual bounded-stress and regularity conditions. This parallels the unattainability form of the third law of thermodynamics. Extremal stationary solutions exist, so the law concerns how an initially nonextremal hole approaches them. The stronger entropy form of Nernst's theorem does not transfer unchanged: extremal holes can have finite area and entropy, depending on their conserved charges.
Quantum field theory turns the mechanical analogy into a physical thermodynamic identification. Hawking radiation is thermal at infinity to leading semiclassical order, with
These give exactly, so the mechanical first law becomes the thermodynamic first law with Bekenstein-Hawking entropy. Scattering through the exterior adds greybody factors to the emitted spectrum but leaves this temperature identification intact. Equilibrium with an external bath and isolated evaporation are different physical states of the field.
During black-hole evaporation, the area can decrease because the semiclassical stress tensor need not satisfy the classical null energy condition. The proposed generalized second law applies instead to , asserting nondecrease of the combined entropy. A Schwarzschild hole illustrates an unusual thermodynamic property: , so its heat capacity is negative. It becomes hotter as it loses energy and is unstable in a simple canonical heat bath, while entropy and energy laws still describe an isolated system. The quantum temperature, area entropy and generalized second law supply the correspondence; the mechanical laws alone do not establish a thermal emission spectrum or a microscopic entropy count.
Use for the stellar gas constant, reserving for the stellar radius. The specific internal energy of a monatomic gas obeying the ideal gas law is . At fixed enclosed mass, the first law of thermodynamics gives the heat available to the outgoing luminosity as when there is no stellar nuclear fusion heating. Therefore
The second coefficient includes both the mass density dependence of and pressure work. This is gravitational energy generation in a homologously contracting ideal-gas star.
Write the characteristic scales as
Then , , , and . Since , direct substitution in the stellar hydrostatic equation and mass conservation gives
For radiative diffusion in a star, dividing the temperature equation by gives
In particular the fourth power here is of , not of the radius.
Because is fixed and is constant, a fixed mass label has fixed . Thus stellar homology gives and . The heating rate reduces to , positive during Kelvin-Helmholtz contraction. Substitution in now yields
The dimensionless profiles fix and . In a common stellar homology family with fixed stellar composition and opacity normalization, the expression for implies
It is independent of radius and constant in time for a particular fixed-mass star. The formula for gives , which integrates to
When the initial-radius term is negligible,
This is the Kelvin-Helmholtz contraction time law in the prescribed stellar homology model. Infinite initial radius is a limiting idealization; a finite initial radius retains the first term.
At main sequence arrival, CNO cycle heating with fixed stellar composition has luminosity
The integral is dimensionless and fixed within the stellar homology family. The constant-opacity radiative scaling still gives . Equating the stellar nuclear fusion and transported luminosities gives
Finally the Kelvin-Helmholtz contraction law evaluated at gives . In a coeval star cluster, the mass just arriving at the main sequence consequently satisfies
These results comprise constant-opacity homologous contraction to CNO ignition; their mass exponents compare fixed-composition models with common dimensionless profiles.
Take any fixed material subvolume of the reference configuration, with outward normal vector , current particle velocity , and positive temperature . The first law of thermodynamics, including kinetic energy, is
The heat-flux sign is negative because is outward. The Second law of thermodynamics supplies the corresponding integral Lagrangian entropy inequality:
The temperature in each supply or flux term is the local temperature, not a common temperature pulled outside the integrals.
Apply the divergence theorem and use the reference momentum conservation equation . In the total-energy equation, the kinetic derivative cancels the body-force work and the term. Since , the remaining local Lagrangian internal energy balance and Lagrangian entropy inequality are
They hold pointwise because was arbitrary. In particular,
The specific Helmholtz free energy is . Substitution yields its balance
and elimination of the heat supply yields the Clausius-Duhem inequality
For , expand the derivative to obtain
The usual Coleman-Noll procedure admits independent reversible choices of and at a fixed state. Their coefficients must vanish, since either sign and arbitrary magnitude would otherwise violate the inequality. Thus
Internal-variable evolution and heat conduction must satisfy this remaining inequality. Substituting the constitutive relations back into the energy balance gives
The last printed equation is missing a factor on its internal-variable term. The PDF defines with per unit mass, so the expression above is the dimensionally consistent consequence. Equivalently define the volumetric thermodynamic force conjugate to an internal variable ; the term is then with no extra density. The two conventions must not be mixed.
For thermoelasticity with a temperature-dependent constraint, admissible rates obey
The reversible coefficient functional need now vanish only on this tangent hyperplane. Hence it can be a scalar multiple of the constraint normal:
Therefore
The reaction terms cancel on admissible rates, leaving the same internal-variable/conduction inequality and entropy heating equation. The multiplier is determined by equilibrium and boundary conditions rather than by the free energy alone.
For linear thermoelasticity with prescribed thermal volume, write , , and . Hold the internal variables fixed in this reversible linearization, and evaluate all derivatives below at and zero reaction. The stress-free reference state has . Since there, first-order nominal and Cauchy stress tensors agree. Their constitutive laws become
where equality of mixed free-energy derivatives gives
The linearized constraint is . Thus volume is thermally prescribed; strict temperature-independent incompressibility corresponds to . The coefficient is an entropy slope; the corresponding volumetric heat capacity at fixed deformation is .
Together with , the first-order thermal balance is . A heat-flux constitutive law completes the evolution problem. The paired signs of the coupling tensor follow from the same free energy, rather than from separately postulated mechanical and thermal equations.
A Lagrangian submanifold of a -dimensional symplectic manifold is an -dimensional submanifold such that the pullback of to is zero. Equivalently, each tangent space is its own symplectic orthogonal. The dimension condition makes it maximally isotropic: in a nondegenerate alternating vector space no isotropic subspace has dimension greater than .
In geometric quantization, a line bundle with a connection on a vector bundle whose curvature is proportional to provides the initial quantum-state space. A real polarization foliates phase space locally by Lagrangian leaves; requiring sections to be covariantly constant along these leaves reduces their independent variables by half. The curvature restricts to zero on a Lagrangian leaf, allowing local flat sections. Globally its holonomy can impose the familiar Bohr-Sommerfeld quantization condition. Thus Lagrangian submanifolds describe choices of quantum configuration data and the support of semiclassical states, rather than requiring simultaneous sharp values of all position and momentum coordinates.
For the thermodynamic application, regard as base coordinates and as fiber coordinates on . The canonical one-form on a cotangent bundle and the position-first symplectic form are
The negative sign in the pressure coordinate is essential. A smooth potential determines the embedding
where the last three coordinates are ordered as : explicitly , , and , so the second fiber coordinate is . Accordingly, in the actual coordinates the graph is
Projection to the base is the identity, so this is an embedded three-dimensional submanifold, even if the Hessian of is degenerate. The first law of thermodynamics gives . Hence
It has half the dimension of the ambient cotangent bundle and is therefore Lagrangian. This proves the thermodynamic potential as a Lagrangian graph construction for every smooth choice of , locally on any domain on which the potential is defined. No invertibility, extensivity or stability assumption is needed for this geometric conclusion.
In components the same vanishing follows from equality of mixed derivatives. The thermodynamic conjugate variables satisfy
These cancel the three coefficients of the pulled-back two-form. They are the integrability relations behind the local thermodynamic potential, and demonstrate directly why the first-law graph is Lagrangian.
Let be the specific internal energy, the specific entropy, the absolute temperature, and the specific enthalpy. The first law of thermodynamics for a simple fluid of fixed composition is , so
The unmagnetized Euler equations for an inviscid fluid are . Define the vorticity and use the vorticity cross-product identity . Substitution gives the Crocco form of the unsteady Euler equation:
Adiabaticity means for smooth flow without dissipative heating. It does not require to be uniform between different fluid elements, so the entropy-gradient term generally remains. For a perfect gas with constant adiabatic exponent , one may write , but the thermodynamic derivation does not require that particular equation of state.
Write the photon energy density as , where by the Stefan–Boltzmann law. The photon gas has pressure and zero photon chemical potential. The first law of thermodynamics therefore gives
Dividing by and integrating gives the entropy density of a photon gas:
The integration constant is fixed to zero at zero temperature, consistently with the third law of thermodynamics.
The inequality says interactions act much faster than the expansion time, maintaining local thermal equilibrium. For adiabatic reversible expansion with no additional source of entropy, cosmological entropy conservation fixes the entropy in a comoving volume. Since that volume scales as , and the relativistic degrees of freedom are fixed here, is constant. Thus
During radiation domination, the Friedmann equation gives , and hence . Fast interactions alone would not exclude an additional entropy-producing process; changing relativistic species likewise changes the simple temperature scaling.
Take a fixed comoving volume, whose physical volume is . For adiabatic cosmic expansion, the first law of thermodynamics gives . Therefore
Spatial curvature does not change the scaling of a fixed comoving volume. With and , the Friedmann equation becomes
Differentiating this expression, rather than dividing by at a fixed point, yields the cosmological density parameter flow
This holds on an expanding branch with . The spatially flat solution is a fixed point. For , . Consequently curvature deviations grow as during radiation domination and as during matter domination. A small present curvature therefore requires a much smaller initial deviation: this is the flatness problem of a decelerating Big Bang. The issue applies to either sign of curvature, not just an open universe. Conversely, accelerated expansion with suppresses small deviations, which is the inflationary solution of the flatness problem.
The stellar adiabatic exponents are fixed-composition, constant-specific entropy derivatives:
The chain rule gives . Define the pressure derivatives and . The first law of thermodynamics, with , gives
Combining this with on an adiabat gives
Thus the specific-heat ratio is not generally equal to the three stellar adiabatic exponents.
For the mixture, and . The specific heat capacity at constant volume, obtained by differentiating at fixed mass density, is
Therefore the adiabatic exponents of a monatomic gas-radiation mixture and its specific-heat ratio are
If a relation involving only the exponents is wanted, eliminate from :
For a pure monatomic perfect gas, and . For any calorically perfect gas with constant heat capacities, the same equality holds with its own . In a genuine gas-radiation mixture, , while the exponents are given separately above.
In the radiation limit, and the adiabatic temperature gradient is . This follows independently from photon entropy: a comoving volume has , so an adiabat obeys and . However, is singular in the pure-radiation limit, not . This is the pure-radiation constant-pressure heat-capacity singularity. The equation fixes temperature whenever pressure is fixed, so an ordinary constant-pressure temperature derivative is not available; along the mixture limit and . For photons alone, mass-specific quantities additionally require a material mass label. The often quoted radiation index is its pressure-density adiabatic exponent, not a finite constant-pressure/constant-volume heat-capacity ratio.
For a fixed amount of gas, the first law of thermodynamics gives . Define the Gibbs free energy , so . Equality of its mixed partial derivatives gives the Maxwell relation
The pipe is adiabatic, so the gas receives no heat. The work done on it at the inlet is , and the work it does at the outlet is . Neglecting kinetic and gravitational energy changes, the first law of thermodynamics gives . Hence the throttling process is isenthalpic:
The pressure drop is irreversible; constant enthalpy does not imply constant entropy.
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 . The first law of thermodynamics gives
Since ,
The differential is therefore exact:
After fixing the irrelevant additive entropy zero, the entropy density at zero chemical potential is
More precisely the first law leaves an additive constant in ; extensivity fixes that constant for the usual entropy-density normalization. A nonzero chemical potential would instead require , 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, and the cosmological perfect-fluid continuity equation gives . Differentiating the extensive expression explicitly, or using the exact differential above, gives
This is cosmological entropy conservation for the closed equilibrium gas with no entropy-producing energy injection.
For radiation in cosmology, , so . Comparing and yields the coefficient ratio
When all relativistic species share the same temperature and count, the familiar normalizations are and . More generally the energy and entropy effective counts can differ; the common 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 is
For an adiabatic process, . Using hydrostatic equilibrium, , the parcel follows
After an upward displacement , its temperature relative to its new surroundings is
Here is the specific heat capacity at constant pressure, and is the buoyancy frequency. At the same pressure, a warmer parcel has lower density and accelerates upward. Thus strict instability requires . Equivalently, the parcel equation is , with
Negative 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.
Let denote specific entropy, reserving from the first question for envelope thickness. The stellar adiabatic exponents are defined at fixed composition by
It is useful to introduce , so .
Take the stellar gas to be a fully ionized, nonrelativistic monatomic ideal gas, with fixed specific gas constant . Its mixture with thermal radiation has
Here is specific internal energy, and the stellar gas-pressure fraction is . The assumption of monatomic gas fixes its heat capacity; the perfect-gas pressure law alone would not determine it.
At fixed and fixed , respectively, the pressure derivatives are
For an isentropic process, the first law of thermodynamics gives . Differentiate the internal energy rather than artificially holding fixed during the perturbation:
The parenthesis is , giving
Since ,
Using the relation between the exponents gives the adiabatic exponents of a monatomic gas-radiation mixture
The requested values, including the adiabatic temperature gradient, are
The pure-radiation and pure-gas rows are understood as limits of the mixture.
For uniform composition, the Schwarzschild criterion for stellar convective instability is
To test whether a radiative configuration becomes unstable, use its required stellar radiative temperature gradient for . Since , the equivalent radial condition is
An outward-displaced parcel then cools less than its new surroundings, remains less dense at the same pressure, and is further accelerated outward.
The requested Eddington-model convective-core mass fraction follows from the usual global constant- closure. At the outer radiative surface, put , , and . Dividing the radiation-pressure gradient by the stellar hydrostatic equation gives
There is no energy generation outside the convective core, so use there. The local radiative-gradient formula is
The same representative has been used in the global luminosity closure and at the boundary. Setting this gradient equal to at produces
This is , , and for , respectively, with a marginal radiation-dominated limit.
There is a consistency qualification to this last model estimate. It cannot be an exact stellar solution if all the printed assumptions are enforced pointwise. An exactly uniform nonzero radiation fraction would give at every radiative point. The exact equations would then require
If is constant throughout an envelope with positive mass density, increases with radius, so this equality cannot hold there. Equivalently uniform fixes the actual gradient to , whereas for . The boxed core fraction is the intended Eddington closure and boundary estimate, which relaxes exact constancy of in the detailed envelope. The constant-beta radiative-envelope obstruction identifies the missing approximation; it is not legitimate to silently assert exact compatibility.
For a slowly evolving spherical star, use a Lagrangian enclosed mass coordinate . Let be specific thermal internal energy and follow a mass element. The stellar energy balance equation is
The outward luminosity is energy transported through the mass shell, is nuclear rest-mass energy released per unit mass and time, and is escaping neutrino power per unit mass. Define the gravothermal stellar energy generation by . At fixed composition it is by the first law of thermodynamics. During composition changes the full internal-energy derivative must use the appropriate equation of state, with nuclear rest-mass changes counted only once. If a tabulated nuclear rate already subtracts reaction-neutrino energy, that loss must not be subtracted a second time in .
In a nearly steady main sequence star, stellar nuclear fusion supplies most of the luminosity. Before stable hydrogen ignition, Kelvin-Helmholtz contraction releases gravitational energy. In a gas-supported hydrostatic star the stellar virial theorem gives , so about half the change in gravitational binding energy heats the star and half is available for radiation. The associated timescale is . After nuclear burning ends, a white dwarf can shine by losing stored thermal energy; crystallization releases latent heat, and composition separation can add gravitational energy. These are heat and gravitational reservoirs rather than sustained hydrogen fusion.
Stellar nuclear fusion is possible because lighter nuclei can combine into products with larger nuclear binding energy. The energy per reaction is , including the relevant particles consistently. Although thermal energies are below the Coulomb barrier, quantum tunnelling permits fusion. The thermonuclear reaction rate averages a cross-section over the thermal relative-speed distribution:
For nonresonant charged-particle reactions the energy weighting contains ; the compromise between the thermal tail and penetration produces the Gamow peak. A resonance can greatly enhance the rate. Density, abundances and temperature therefore all affect the stellar energy-generation rate, and power laws such as or are local approximations.
For the proton–proton chain, first produce deuterium and helium-3:
The initial reaction is slow because it converts a proton to a neutron through the weak interaction. It is the bottleneck that allows long hydrogen-burning lifetimes. The alternative pep reaction also feeds the same chain. There are three principal proton-proton chain branches. In pp I, two helium-3 nuclei terminate the chain:
In pp II, helium-3 first captures an existing helium-4 nucleus, then electron capture and proton capture finish the branch:
In pp III, the beryllium-7 instead captures a proton:
Relative branch weights change with temperature and composition. All branches convert four protons into a net helium-4 nucleus, with two weak conversions and two neutrinos, allowing for the electron consumed in pp II. After including positron annihilation, the atomic-mass energy budget is approximately per net helium-4 nucleus, but the deposited heat is smaller by the branch-dependent escaping-neutrino energy. Near ordinary low-mass main-sequence conditions, is a useful local approximation.
The CNO cycle provides an alternative catalyzed hydrogen burning route in hotter cores. The CNO-I cycle, also called the CN cycle, consists of
The carbon seed is regenerated: the net reaction again consumes four protons and produces one helium-4 nucleus, two positrons and two neutrinos. The slow reaction governs the ordinary cycle and causes nitrogen-14 to accumulate. The CNO-II cycle branches through , , , and , after which the main cycle continues. Thus CNO nuclei are catalytic in the closed cycles, although their relative abundances change. A local rate is with a much steeper exponent than the PP chain; the preceding homology problem specifies . The concentration of this heating near the centre favors a convective core in hotter, more massive hydrogen-burning stars.
Once hydrogen is exhausted in a core, contraction can raise its temperature enough for core helium burning. The Triple-alpha process overcomes the absence of stable mass-5 and mass-8 nuclei by maintaining a tiny transient beryllium-8 population:
The Hoyle state, an excited carbon-12 resonance near the three-alpha threshold, makes the second capture efficient enough despite beryllium-8's very short lifetime. The net conversion is , releasing about . Because three alpha particles are required, the specific rate scales as ; a schematic resonant rate is
Its local temperature exponent is , of order forty near . A degenerate core cannot expand promptly to regulate this increase, giving a helium flash; in a nondegenerate core the stellar thermostat permits stable helium burning. The competing reaction determines much of the eventual carbon-oxygen mixture.
If the core becomes hot enough, carbon burning follows, for example through and . Neon burning begins with photodisintegration, , followed by alpha capture such as . The first step consumes heat, but the combined rearrangement can release net energy. Oxygen burning includes and channels producing phosphorus and other nearby nuclei. At still higher temperatures, silicon burning proceeds through photodisintegration and particle captures in a reaction network approaching quasi-statistical equilibrium, producing iron-group nuclei. It is not simply direct silicon-plus-silicon fusion; the products depend on the electron fraction and weak-interaction timescale.
The increasing nuclear binding energy per nucleon supplies progressively less energy as products approach the iron group. Further fusion past that region cannot provide sustained net heat to support a core. Massive stars can therefore end with an unstable iron-group core; lower-mass stars do not attain all these burning stages and instead leave remnants such as carbon-oxygen white dwarfs. Heavy-element neutron captures can synthesize nuclei beyond the iron group without being the principal hydrostatic power source.
Finally, stellar neutrino energy loss competes with every heat source. Reaction neutrinos accompany the proton–proton and CNO chains. Hot or dense material also emits neutrino pairs through , plasmon decay, , and electron-ion bremsstrahlung. These neutrinos usually escape, unlike the photons whose transport is diffusive; neutrino trapping requires much more extreme collapse conditions. Late burning has smaller fuel-energy reservoirs and strong neutrino cooling, so its duration is much shorter than core hydrogen burning. The surface luminosity is the escaping photon power, , not the sum of photons and unobserved neutrino power. Radiative diffusion in a star, convection and electron conduction redistribute energy, while the energy-balance equation distinguishes true sources, losses and storage. Accretion or mass loss adds boundary energy and mechanical work when present, rather than changing the nuclear reaction budget.