Stellar evolution

Stellar evolution is the change of a star's composition, structure and observable properties over time. Nuclear fusion changes the core's composition and available fuel; gravity, energy transport and mass loss then alter the stellar radius, luminosity and effective temperature. The initial mass strongly influences the sequence, but composition, rotation and binary mass transfer also matter.

The evolution of a solar-mass star connects the main sequence to giant phases and eventually a white dwarf. More massive stars can reach stellar core collapse. An evolutionary track in the Hertzsprung-Russell diagram makes these changes visible, while comparisons with coeval stellar populations test the predicted lifetimes. Two stars with similar present appearance can have different prior histories, especially after binary interaction.

Rapid contraction of a stellar core when its pressure support becomes inadequate. An evolved iron-group core can lose support through electron capture and photodisintegration and collapse towards a neutron star or black hole. A successful outgoing explosion produces a core-collapse supernova; core collapse need not always give a bright successful supernova.

Core hydrogen fusion supports the main sequence. Hydrogen-shell fusion around an inert helium core then produces the subgiant and red giant phases. Core helium fusion gives the red clump, followed by hydrogen- and helium-shell burning on the Asymptotic giant branch. Envelope loss can expose a hot planetary nebula nucleus before a carbon-oxygen white dwarf cools. This sequence concerns a roughly solar-mass star; initial mass, composition and mass loss affect detailed tracks and timescales.

A compact hot evolved star with little remaining envelope. A hot subdwarf B is commonly interpreted as a core-helium-burning object with a very thin hydrogen envelope; spectral subtype and evolutionary state must be specified.

A blue straggler lies above and blueward of the normal main-sequence turnoff of its associated old stellar population. Binary mass transfer and stellar mergers can give an old system a rejuvenated, more massive hydrogen-burning star. Its position therefore need not imply recent ordinary star formation.

A post-helium-ignition stage of low-mass stellar evolution with core helium nuclear fusion and usually a hydrogen-burning shell. The branch spans a broad range of effective temperatures at comparatively similar luminosities, largely because envelope masses vary.

The concentration of core-helium-burning stars on the cool part of the horizontal branch. It must be distinguished from the ascending hydrogen-shell-burning red giant branch.

The phase of stellar evolution in which hydrogen-shell nuclear fusion surrounds an inert helium core and the envelope expands and cools. Low-mass stars climb toward helium ignition with growing core mass and luminosity; their cores become supported by electron degeneracy pressure. A red giant is a broader observational category and need not be on this particular evolutionary branch.

After substantial envelope loss on the Asymptotic giant branch, a low- or intermediate-mass star's exposed core heats while its radius shrinks, often at approximately constant luminosity during continued shell burning. It may ionize the expelled envelope as a planetary nebula. After the remaining nuclear fuel is depleted, its luminosity decreases and it joins a white dwarf cooling track. Evolutionary tracks on the Hertzsprung-Russell diagram connect the central-star temperature and luminosity with core mass and elapsed time.

The time for stellar nuclear fusion to consume the accessible fuel at characteristic luminosity $L$. If a fraction $f$ of a star's mass $M$ supplies fuel abundance $X$ and released energy per unit fuel mass $E_H$, hydrogen burning gives $t_{\mathrm{nuc}}\sim fXM{E}_H/L$. Accessible core fuel, rather than the entire star, sets the lifetime unless mixing is complete. This is longer than a structural Kelvin-Helmholtz cooling time when the nuclear reservoir greatly exceeds the gravitational contraction reservoir.

A stellar remnant is the object left after a star finishes its principal nuclear-burning evolution, such as a white dwarf, neutron star or black hole. Stellar remnants retain mass while often contributing little optical light, affecting an old stellar population's mass-to-light ratio.

A dense stellar-scale object such as a white dwarf, neutron star or stellar-mass black hole. accretion releases a large gravitational binding energy per unit mass because the accreting surface or inner orbit is small. A black hole has no material surface, so its radiative yield depends on the accretion flow rather than surface impact.

The stellar age is the elapsed time since a star formed. A constant number formation rate on a finite interval gives a uniform distribution of present ages if every formed object remains counted, including stellar remnants. An age distribution inferred only from surviving main sequence stars generally differs because their lifetimes depend on mass.

A cosmological age must exceed the ages of its oldest stellar populations plus their formation interval. globular cluster ages use the main-sequence turnoff with modeled distances, composition, nuclear rates and stellar transport. Their errors are correlated with these inputs, so a lower age bound must propagate the stellar and distance uncertainties rather than treating an isochrone as an exact clock.^[ref]

Estimating astronomical ages from radioactive abundance ratios. A ratio evolves as $Q=Q_0{e}^{-({\lambda }_1-{\lambda }_2)t}$ for decay constants ${\lambda }_1,{\lambda }_2$. The initial production ratio and subsequent enrichment or mixing must be constrained; an observed ratio alone does not specify a unique age.

A planetary nebula is expelled stellar gas illuminated and ionized by a hot central remnant after an evolved star loses its outer envelope. It belongs to late stellar evolution, despite its historical name suggesting a planet. Ultraviolet photons from the exposed hot core cause photoionization; subsequent atomic processes produce prominent emission lines.^[ref]

The planetary-nebula spectral energy distribution can combine those lines with continua and thermal emission from dust. Spectroscopy therefore measures more than a colour: line ratios diagnose gas conditions, while Doppler shifts trace expansion. The Zanstra method connects nebular light to the central star's ionizing output. The central remnant subsequently follows white-dwarf cooling as the surrounding gas disperses.

The spectral energy distribution of a planetary nebula combines hot central-star radiation, nebular free-free/free-bound/two-photon continua, recombination and collisionally excited lines, infrared dust emission and sometimes molecular emission. Ionizing-photon supply, gas thermal/ionization balance, density, composition, geometry and extinction determine their relative strengths. A line-dominated broad-band flux is not by itself a continuum temperature measurement.

Dust contributes a broad infrared continuum together with composition-dependent vibrational structure. Carbon-rich material may show aromatic bands, while oxygen-rich silicates can show features near 9.7 and 18 micrometres. Grain size, charging, transient heating, optical depth and chemical processing affect these features, so one dust temperature cannot characterize every band carrier.

A nebular free-bound continuum is emitted when a free electron is captured into a bound atomic level and its excess energy leaves as a photon. The threshold energies produce recombination edges such as the Balmer discontinuity. Its spectrum depends on electron temperature, ionic abundances, capture cross-sections and emission measure; it is distinct from bound-bound recombination lines and free-free bremsstrahlung.

The Zanstra method estimates a central star's effective temperature by comparing its observed nonionizing continuum with the absorbed ionizing-photon rate inferred from nebular recombination lines. If ${\Phi }_{E_0}(T)={\int }_{E_0/h}^{\infty }{B}_{\nu }(T)/(h\nu )d\nu$, then $[F_{\ell }/(h{\nu }_{\ell }{F}_{{\nu }_o,\ast })][\beta /{\alpha }_{\ell }^{\mathrm{eff}}]={\Phi }_{E_0}(T_s)/B_{{\nu }_o}(T_s)$ when the line traces the full appropriate ionizing budget. Both the distance and stellar radius cancel. Density inhomogeneity also cancels from the integrated recombination estimate when the coefficient ratio is known; photon leakage, dust and mixed-species transfer can invalidate the simple budget.

A carbon-oxygen core is the product of core helium burning. After central helium exhaustion it has little hydrogen or helium, while helium-rich and hydrogen-rich layers surround it.

A blue loop is an excursion to higher effective temperature, followed by a return redward, during core helium burning in some evolved stars. Its extent depends on core and envelope structure, composition gradients and mixing. Only loops reaching the Cepheid instability strip produce the associated additional Cepheid variable crossings.

The Hertzsprung gap is a sparsely populated region of the Hertzsprung-Russell diagram traversed rapidly by many stars evolving away from the main sequence toward the red giant region.

Kelvin-Helmholtz contraction releases gravitational binding energy as a self-gravitating body cools and shrinks. It powers young stars and giant planets before nuclear burning or degeneracy dominates their thermal evolution.

For a monatomic gas obeying the ideal gas law, the specific internal energy is $e=3P/(2\rho )$. With no stellar nuclear fusion heating, the first law of thermodynamics makes the effective stellar energy-generation rate $ϵ=-{\partial }_{t}e-P{\partial }_t(1/\rho )$, with the derivatives taken at fixed enclosed mass. Differentiating gives $ϵ=-3{\partial }_{t}P/(2\rho )+5P{\partial }_t\rho /(2{\rho }^2)$. In a fixed-mass stellar homology Kelvin-Helmholtz contraction, $\rho \propto R^{-3}$ and $P\propto R^{-4}$ at each mass label. Hence $ϵ=-3P\dot{R}/(2\rho R)>0$ when $\dot{R}<0$. This includes internal energy storage as well as pressure work; using the gravitational potential-energy decrease alone would overestimate the outgoing luminosity.

The intrinsic planetary luminosity $L_{\mathrm{int}}$ is the power supplied by interior cooling, contraction, nuclear reactions if present, or internal heating, excluding reflected or passively reradiated stellar light. For fixed mass and no accretion or nuclear source, global conservation of energy gives $L_{\mathrm{int}}=Q_{\mathrm{deep}}-d(U+\Omega )/dt$, where $Q_{\mathrm{deep}}$ is any additional power deposited in the interior, such as tidal heating or Joule heating. For $Q_{\mathrm{deep}}=0$ and a gas-supported quasi-static sphere of fixed specific-heat ratio $\gamma$ and negligible surface pressure, the stellar virial theorem gives $\Omega =-3(\gamma -1)U$, henceFor a dynamically stable gas-supported model take $\gamma >4/3$. For $\gamma =5/3$ this is half the gravitational release. Degenerate giants require a different internal energy reservoir; the half-factor is not universal.

With negligible nuclear or external heating, the stellar structure equations give $L_{\mathrm{int}}=-{\int }_0^{M}T\partial s/\partial tdm$. A nearly isentropic convective interior therefore has $\dot{s}=-L_{\mathrm{int}}/\int Tdm<0$ whenever it radiates. Decreasing specific entropy does not require temperature at every depth to decrease immediately: contraction of an ideal gas can raise its central temperature while total energy is lost.

If otherwise comparable isolated planets of one mass follow $L_{\mathrm{int}}=C{t}^{-\alpha }$, their intrinsic luminosity ratio is $L_{\mathrm{int}}(t_1)/L_{\mathrm{int}}(t_2)=(t_2/t_1{)}^{\alpha }$. The normalization depends on composition and initial entropy. Strong stellar irradiation, deep heating, or accretion can invalidate extrapolating the same cooling law.

The Kelvin-Helmholtz cooling time estimates how long a self-gravitating body can radiate a reservoir of energy of order its gravitational binding energy. It uses intrinsic luminosity $L_{\mathrm{int}}$, excluding reradiated stellar light. At fixed radius and structure, this reservoir scales as $M^2$; the approximation changes when degeneracy or other internal physics changes the accessible energy.

A pre-main-sequence star contracts toward the main sequence before sustained core hydrogen burning supplies its luminosity. Convection, accretion, and deuterium burning shape its track across the Hertzsprung-Russell diagram.

A pre-main-sequence star with a developing radiative core can contract along a relatively horizontal path toward higher effective temperature in the Hertzsprung-Russell diagram. This differs from the nearly vertical fully convective Hayashi track.

A Hayashi track is the nearly vertical path of a fully convective pre-main-sequence star on the Hertzsprung-Russell diagram. Photospheric opacity fixes a weak relation between luminosity and effective temperature.

Match a fully convective stellar polytrope of index $3/2$ to the grey-atmosphere convection onset with eighth-power temperature opacity. At the fixed matching optical depth, $T_b\propto T_e$ and $P_b\propto g^{2/3}{T}_e^{-16/3}$. The ideal-gas adiabat $P=K{\rho }^{5/3}$ gives $K\propto T_b^{5/3}{P}_b^{-2/3}\propto T_e^{47/9}{g}^{-4/9}$. The polytropic mass-radius relation gives $K\propto R{M}^{1/3}$; with $g\propto M/R^2$, this yields $T_e^{47}\propto R{M}^7$. The Stefan–Boltzmann law then gives the displayed luminosity. Composition and opacity normalization are held fixed between models.

Match an efficient, monatomic gas, fully convective star to the grey-atmosphere convection onset with thirteenth-power opacity. Its convective relation is $P=K_T{T}^{5/2}$. At the matching temperature $T_b\propto T_e$, the stellar atmosphere gives $K_T^2\propto M{L}^{-1}{T}_e^{-13}\propto M{R}^{13/2}{L}^{-17/4}$ using the Stefan–Boltzmann law. The ideal-gas relation converts this to $P=K_{\rho }{\rho }^{5/3}$ with $K_{\rho }=(\mathcal{R}/\mu {)}^{5/3}{K}_T^{-2/3}$. The polytropic mass-radius relation at index $3/2$ gives $K_{\rho }\propto G{M}^{1/3}R$, hence $K_T^2\propto M^{-1}{R}^{-3}$. Equating the two expressions proves $L^{17/4}\propto M^2{R}^{19/2}$, or $L\propto M^{8/17}{R}^{38/17}$. The entropy coefficient is constant throughout the convective envelope, but can change between models.

For a monatomic fully convective ideal gas star, uniform specific entropy gives $P=K{T}^{5/2}$ with $K\propto {\mu }^{-5/2}{M}^{-1/2}{R}^{-3/2}$. Combining this with the stellar surface boundary condition $\kappa P=2g/3$ and $\kappa ={\kappa }_0{\rho }^a{T}^b$ givesAll composition comparisons here hold ${\kappa }_0$ fixed. For $a=1,b=4$, $T_e^8\propto {\mu }^4{M}^{2}R$, and the Stefan–Boltzmann law gives $L\propto {\mu }^{2}M{R}^{5/2}$. The stellar virial theorem and Kelvin-Helmholtz contraction then give $R^{-7/2}-R_0^{-7/2}\propto {\mu }^2{M}^{-1}(t-t_0)$, with a positive coefficient. Once the initial-radius term is negligible, $R\propto {\mu }^{-4/7}{M}^{2/7}(t-t_0{)}^{-2/7}$.

For a fixed-mass, fully convective monatomic gas protostar, write $P=K{T}^{5/2}$ and $K=K_0{\mu }^{-5/2}{M}^{-1/2}{R}^{-3/2}$. Matching the photosphere to $\kappa P=2g/3$ with $\kappa ={\kappa }_0\rho T^4$ gives $T_e^8=2G\mathcal{R}{\mu }^4{M}^{2}R/(3{\kappa }_0{K}_0^2)$. The Stefan–Boltzmann law gives $L=\pi ac[2G\mathcal{R}/(3{\kappa }_0){]}^{1/2}{K}_0^{-1}{\mu }^{2}M{R}^{5/2}$. The stellar virial theorem gives total energy $-3G{M}^2/(7R)$, so $L=-3G{M}^2\dot{R}/(7R^2)$. Integrating proves $C=(49\pi ac/6)[2\mathcal{R}/(3{\kappa }_{0}G){]}^{1/2}{K}_0^{-1}$. This holds before stellar nuclear fusion, changing opacity or electron degeneracy pressure invalidates the assumed contraction law.

Match an $n=3/2$ stellar polytrope to the grey-atmosphere convection onset with seventeenth-power opacity. At each instant the convective entropy coefficient obeys $P=K{T}^{5/2}$, $K\propto M^{-1/2}{R}^{-3/2}$ at fixed composition, while the boundary gives $K\propto g^{1/2}{T}_e^{-21/2}$. Thus $T_e^{21/2}\propto M{R}^{1/2}$ and the first displayed luminosity law follows. The stellar virial theorem for a monatomic gas gives total energy $-3G{M}^2/(7R)$, so Kelvin-Helmholtz contraction has $\dot{R}=-7L{R}^2/(3G{M}^2)$. Integrating at fixed mass gives $R^{-67/21}-R_0^{-67/21}\propto t$. The second displayed relation holds after the initial-radius term is negligible. The entropy coefficient is spatially constant but changes during contraction.

For convective contraction with power-law opacity $\kappa \propto \rho T^4$ at fixed stellar mass and opacity normalization, $T_e\propto {\mu }^{1/2}{R}^{1/8}$ and $L\propto {\mu }^2{R}^{5/2}$. Eliminating $R$ gives $L\propto {\mu }^{-8}{T}_e^{20}$. Thus a larger mean molecular weight moves the track to larger effective temperature at a given luminosity, with $T_e{|}_L\propto {\mu }^{2/5}$. Contraction decreases both temperature and luminosity for this particular opacity exponent.

Matter accreting at rate $\dot{M}$ onto a star of mass $M$ and radius $R$ releases gravitational power of orderup to an efficiency and an inner-boundary correction.

Integrating $L_{\mathrm{acc}}=GM\dot{M}/r_0$ at fixed radius gives this result. With negligible seed mass it supplies the radiative requirement of monatomic cloud collapse only if ${\Omega }_f=-G{M}^2/r_0$ under the adopted zero-initial-energy approximation. A different density-profile coefficient, retained accretion heat, a wind or a changing radius modifies the budget. Positive binding-energy magnitude and negative potential energy must not be interchanged silently.

For constant supply, negligible seed mass and fixed emitting radius, $M=\dot{M}t$ and $L_{\mathrm{acc}}\propto t$. Thermal emission with $L=4\pi r_0^2{\sigma }_{\mathrm{SB}}{T}_{\mathrm{eff}}^4$ therefore gives the quarter-power law. It describes an effective surface temperature. A homologous ideal-gas interior at the same fixed radius instead has virial temperature proportional to mass and hence to time; its absolute temperature requires a structural model.

The radiative efficiency is the escaping radiative power divided by an adopted accreted rest-mass energy rate. Its meaning depends on where the mass rate is measured and which emission is counted. For efficient Newtonian accretion onto a stellar surface it scales as $GM/(R_{\ast }{c}^2)$, whereas advection or outflows can reduce luminosity for the same supplied mass.

At the temperatures relevant to pre-main-sequence stars, a useful steep power-law approximation to the specific deuterium-burning rate is ${ϵ}_D\propto \rho T^{15}$. Its sensitivity to central temperature makes deuterium burning act as a temporary thermostat.

A subgiant has begun evolving away from the main sequence after central hydrogen exhaustion, increasing its radius before the extended red giant stage.

A red giant is an evolved star with a large radius and low effective temperature. On the first giant branch, a low- or intermediate-mass star is powered primarily by shell hydrogen burning around a helium-rich core.

Use a Kramers radiative-zero envelope around a stellar core and specific heating $ϵ={ϵ}_0\rho T^{67/4}$. Since $\rho \propto C{T}^{13/4}$ and $T\propto M_c/r$, $ϵ\propto C{r}^{-20}$ and $ϵ(1.05R_c)/ϵ(R_c)=1.05^{-20}\simeq 0.377$. The volume heating is $\rho ϵ\propto C^2{T}^{93/4}$. Integrating $4\pi r^2\rho ϵ$ from $R_c$ outwards gives $L\propto C^2{M}_c^{93/4}{R}_c^{-81/4}$, while radiative transport gives $C^2\propto M_c/L$. Therefore $L^2\propto M_c^{97/4}{R}_c^{-81/4}$, or $L\propto M_c^{97/8}{R}_c^{-81/8}$. For a core sequence with fixed $R_c$, the mass exponent is $97/8$. The steep integrand confines most heating to a thin shell; using the exterior constant-luminosity profile throughout it estimates the coefficient rather than resolving the shell structure.

In a thin-shell giant model with a nonrelativistic degenerate core, $R_c\propto M_c^{-1/3}$. For a gas-pressure radiative layer with $\kappa \propto \rho T^{-3}$ and specific energy generation $ϵ\propto \rho T^{10}$, the deep profile has $P\propto T^4$, $\rho \propto T^3$ and $T\propto M_c/r$. Integrating the strongly concentrated shell power gives $L\propto C^2{M}_c^{16}{R}_c^{-13}$, while radiative transport gives $C^2\propto M_c/L$. Therefore $L^2\propto M_c^{17}{R}_c^{-13}$, yielding the displayed exponent. This is a model-specific scaling, not a universal core-mass law; it fails when envelope stripping removes the shell fuel or the confining structure.

A negligible-mass radiative envelope surrounding a core of mass $M_c$, carrying constant luminosity and having constant opacity, obeys $dP/dT=16\pi acG{M}_c{T}^3/(3{\kappa }_{0}L)$. Neglecting the outer pressure integration constant gives $\Lambda =4\pi acG{M}_c/(3{\kappa }_{0}L)$. The gas pressure fraction $\beta$ is constant because $P_{\mathrm{rad}}=a{T}^4/3$, and $\rho =\mu \beta \Lambda T^3/\mathcal{R}$. Hydrostatic balance gives $T=T_s+\mu \beta G{M}_c(1/r-1/R_s)/(4\mathcal{R})$; dropping this outer-boundary constant gives the displayed deep-envelope approximation. A positive gas-pressure fraction requires $L<4\pi cG{M}_c/{\kappa }_0$.

With specific heating $ϵ={ϵ}_0\rho T^{16}$, the constant-opacity radiative-zero giant envelope has $T\propto \mu \beta M_c/r$ and $\rho \propto \mu \beta M_c{T}^3/L$. Therefore the shell power $4\pi r^2\rho ϵ$ is proportional to $(\mu \beta M_c{)}^{24}{L}^{-2}{r}^{-20}$. Integrating from the core radius gives $L^3\propto (\mu \beta M_c{)}^{24}{R}_c^{-19}$. The steep radial dependence justifies a thin-shell estimate dominated by the base. Using the outer radiative profile throughout the emitting layer is an approximation; a resolved shell can change the numerical coefficient. Also $\beta =1-{\kappa }_{0}L/(4\pi cG{M}_c)$, so retaining radiation pressure makes the relation implicit in luminosity.

Stellar rejuvenation occurs when mixing brings fresh hydrogen into the core of a mass-gaining or merged main sequence star, making its remaining nuclear lifetime resemble that of a younger star. The amount of rejuvenation depends on the accretor's internal mixing and prior composition.

A helium star has lost most of its hydrogen-rich envelope and is made predominantly of helium. If massive enough, its core can ignite helium; a lower-mass stripped core instead becomes a helium white dwarf.

An asymptotic-giant-branch star has an inert carbon-oxygen core surrounded by helium- and hydrogen-burning shells and a deep convective envelope. Alternating shell burning produces thermal pulses and dredge-up, while strong mass loss removes the envelope.

The fate of an Asymptotic giant branch star depends on whether its envelope is removed before its core reaches the next ignition or instability threshold. The approximate envelope-loss time is $t_{\mathrm{env}}=M_{\mathrm{env}}/|\dot{M}|$, while the time to add a required core mass $\Delta M_c$ is $t_{\mathrm{grow}}=\Delta M_c/{\dot{M}}_c$. Shell stellar nuclear fusion grows the core; third dredge-up can reduce the net growth. If $t_{\mathrm{env}}\ll t_{\mathrm{grow}}$, shell fuel disappears first and a white dwarf remnant is favored. The opposite ordering permits further core growth but does not by itself specify whether carbon burning is stable, explosive, or followed by further evolution. Both rates change with time, so the estimates guide rather than replace a stellar model.

When the bottom of an Asymptotic giant branch star's convective envelope becomes hot enough for hydrogen burning, the CNO cycle can convert newly dredged-up carbon into nitrogen. Its onset and strength depend on mass, composition and envelope structure; it is not a universal property of all stars at one initial mass.

Stellar dredge-up is mixing that brings material altered by nuclear reactions from deeper layers into a star's outer envelope and hence its observable surface. As a convective envelope changes its depth during stellar evolution, it can reach layers with a different composition. Surface spectroscopy can then reveal processing that occurred far below the photosphere.

The numbered events distinguish evolutionary circumstances rather than a single repeated clock. First dredge-up exposes products of prior hydrogen burning during giant evolution; third dredge-up can carry carbon and neutron-capture products outward after a thermal pulse of an asymptotic-giant-branch star. The measured abundance therefore records both prior reactions and subsequent transport. It need not reproduce the composition with which the star was born.

First dredge-up occurs as a star ascends the red-giant branch. Its deepening convective envelope exposes material processed by hydrogen burning, especially enhanced helium and nitrogen and depleted carbon.

Second dredge-up occurs after core helium exhaustion in sufficiently massive intermediate-mass stars. The convective envelope moves inward, mixes helium- and hydrogen-burning products to the surface, and reduces the hydrogen-exhausted core mass.

Third dredge-up follows a thermal pulse on the asymptotic giant branch, when the convective envelope reaches into the intershell and can carry freshly made carbon and slow-neutron-capture elements to the surface.

A thermal pulse is a helium-shell flash caused by thermally unstable burning in a thin shell. It briefly drives intershell convection, expands and cools the hydrogen-burning shell, and can be followed by third dredge-up.

The Härm–Schwarzschild instability is the thin-shell thermal instability of helium burning in an asymptotic-giant-branch star. Heating greatly raises the nuclear rate while the geometrically thin shell initially expands too little to provide stabilizing cooling.

A super-AGB star is massive enough to ignite carbon and form a degenerate oxygen-neon core, but initially too light to proceed directly through all advanced burning stages. Its fate depends on competition between core growth and envelope mass loss.

An oxygen-neon-magnesium core is the degenerate product of carbon burning, composed mainly of oxygen and neon with a smaller magnesium abundance. It can become an oxygen-neon white dwarf or collapse after electron captures if it approaches the Chandrasekhar mass.

An electron-capture supernova occurs when electron captures, especially on neon and magnesium nuclei in a degenerate oxygen-neon core, remove pressure support and trigger collapse before ordinary iron-core formation.

A stellar wind is a continuous outflow of gas from a star. Radiation, gas pressure, waves, or magnetic fields can supply its acceleration, and the resulting mass loss can determine a star's final remnant.

A steady outward radial flow with negligible pressure force conserves $v^2/2-GM/r$. If its speed at a reference radius equals the mechanical escape velocity, this constant is zero and $v\propto r^{-1/2}$. Conservation of mass gives $\rho \propto r^{-3/2}$. An adiabatically evolving ideal gas carried along that leading flow has $p\propto r^{-3\gamma /2}$ and $T\propto r^{-3(\gamma -1)/2}$.

Finite specific enthalpy changes the dynamical problem: the full conserved Bernoulli equation is $v^2/2+h-GM/r=B$. At the same boundary mechanical escape speed, $B=h_e$ rather than zero. A warm wind therefore needs its enthalpy and equation of state specified before the explicit ballistic velocity can be replaced by a consistent pressure-driven solution.

An episode of particularly strong mass loss near the end of Asymptotic giant branch evolution, capable of removing the remaining envelope rapidly compared with its prior evolution. Radiation acting on dust can transfer momentum to the gas; pulsations help lift material into cooler dust-forming layers. Envelope removal terminates sustained hydrogen-burning shell and helium-burning shell activity.

A magnetized stellar wind extracts angular momentum through matter transport and Maxwell stress tensor. The Alfvén radius supplies an effective lever arm larger than the stellar radius, with geometry factors absorbed into an averaged lever arm. Dynamo action couples the torque to rotation by supplying the magnetic field. For a simplified radial wind of fixed speed, $R_A^2$ scales as $B_{\ast }^2{R}_{\ast }^4/({\mu }_0{\dot{M}}_w{v}_w)$, so unsaturated $B_{\ast }\propto \Omega$ gives a torque proportional to $-{\Omega }^3$. Magnetic topology and wind thermodynamics can alter this scaling.

If the moment of inertia I changes slowly and wind-driven magnetic braking of a solar-type star has torque $-K{\Omega }^3$, then ${\Omega }^{-2}(t)={\Omega }^{-2}(t_0)+2K(t-t_0)/I$. At sufficiently late times within this unsaturated approximation, rotation scales as $t^{-1/2}$. This empirical main-sequence trend is not a universal law for saturated rapid rotators, changing stellar structure or every magnetic geometry.

Using an axisymmetric hydrodynamic mass flux function, set $s=s(\psi )$, $\ell =R{u}_{\varphi }=\ell (\psi )$ and $\epsilon =\epsilon (\psi )$. The azimuthal vorticity is ${\omega }_{\varphi }=-R\nabla \cdot (\nabla \psi /(\rho R^2))$. The meridional projection of the Crocco form of the unsteady Euler equation gives the displayed equation. The mass density also satisfies the Bernoulli relation $\epsilon =|\nabla \psi {|}^2/(2{\rho }^2{R}^2)+{\ell }^2/(2R^2)+w+\Phi$, so an equation of state is part of closing the system.

A steady spherical polytropic flow escaping a point-mass potential conserves the outward mass loss rate and the Bernoulli function $C=u^2/2+c_s^2/(\gamma -1)-GM/r$. For positive terminal speed, a finite-radius smooth sonic point requires $C>0$ and $1<\gamma <5/3$. Its sound speed obeys $c_{s,c}^2=2C(\gamma -1)/(5-3\gamma )$, and evaluating the conserved flux there selects the transonic mass loss rate. Other parameter limits need separate treatment.

Put $y=\rho r^{4/(\gamma +1)}$ and $\alpha =4(\gamma -1)/(\gamma +1)$. The polytropic stellar wind equations become $F(r)=C{r}^{\alpha }+GM{r}^{\alpha -1}=G(y)=\gamma K{y}^{\gamma -1}/(\gamma -1)+{\dot{M}}^2/(32{\pi }^2{y}^2)$. At a regular sonic point both derivatives vanish. Matching their minima allows the two positive-$y$ branches to join with finite slope; unequal minima either leave a forbidden radial interval or separate the subsonic and supersonic branches. This reduces a singular differential-equation crossing to a geometric minimum comparison.

In an optically thick steady spherical stellar wind of constant mean molecular weight ${\mu }_{\mathrm{mol}}$, let $v$ be the outward speed and $a_g^2=k_{B}T/({\mu }_{\mathrm{mol}}{m}_u)$ the isothermal sound speed squared. Mass conservation gives ${\rho }^{\prime }/\rho =-2/r-v^{\prime }/v$, where a prime denotes differentiation with respect to radius $r$. Use $P_{\mathrm{gas}}=\rho a_g^2$ and radiative diffusion to write $P_{\mathrm{rad}}^{\prime }/\rho =-g\Gamma$, where $g=G{M}_r/r^2$ and $\Gamma =L_r/L_{\mathrm{Edd}}$ uses the local Eddington luminosity $L_{\mathrm{Edd}}=4\pi cG{M}_r/\kappa$ and opacity $\kappa$. Substituting into the momentum equation yields the displayed wind equation. Smooth sonic point crossing requires both its denominator and numerator to vanish. The radiation-pressure gradient already supplies radiative acceleration; adding that force separately would count it twice.

For an optically thick diffusion-controlled wind, put $\beta =P_{\mathrm{gas}}/(P_{\mathrm{gas}}+P_{\mathrm{rad}})$. Diffusion gives $(\log T{)}^{\prime }=-g\Gamma \beta /[4a_g^2(1-\beta )]$. At $\Gamma =1$, the numerator of the gas and radiation pressure stellar wind equation is $2a_g^2/r+g\beta /[4(1-\beta )]>0$ for finite-temperature gas and $0<\beta <1$. A subsonic branch then decelerates and cannot pass a regular sonic point. This statement is specific to subsonically launched winds in this transport model, not to every supersonic solution or every form of stellar-wind driving.

A neutron star is a compact stellar remnant whose gravity compresses roughly a stellar mass into a radius of order ten kilometres. Degenerate matter supplies most of its pressure support, while rapid rotation and intense magnetic fields can power observable emission.

The mass-shedding limit is the spin at which equatorial surface matter can no longer remain attached in equilibrium. Newtonian gravity gives a Keplerian frequency scale at the surface. Rotation-induced deformation, general relativity and the equation of state change the numerical bound. A magnetic-equilibrium formula is not valid indefinitely as the exterior magnetosphere shrinks inside the star.

A pulsar is a rotating magnetized neutron star whose directed emission is observed as periodic pulses. Rotation supplies the energy of an ordinary radio pulsar.

A pulsar turns stellar rotation into an observable clock: the neutron star's emission can continue while its beam sweeps past the observer only briefly. The pulse spacing measures a repetition period; identifying that period with stellar rotation requires a physical model. The discovery of the first pulsar established a remarkably regular astronomical signal before its engine was understood. Today, pulsar dispersion delay also makes these clocks probes of the intervening interstellar medium.

Loss of a rotating neutron star's angular momentum and rotational energy through electromagnetic torque and a magnetized outflow. For constant moment of inertia $I$, differentiating $E_{\mathrm{rot}}=I(2\pi /P{)}^2/2$ gives the displayed positive loss rate when $\dot{P}>0$. It is not the same as the measured radiation luminosity in any one observing band.

Extrapolate an anchored magnetic dipole field to the light cylinder $R_c=c/\Omega$ and estimate outward power as magnetic energy density times $c{R}_c^2$. The result fixes the displayed scaling but not the geometric coefficient. An aligned vacuum dipole does not radiate, while a current-carrying plasma magnetosphere can have nonzero aligned spin-down. Corotation at the light cylinder is a notional limiting construction, not a material surface moving rigidly at light speed.

For a power-law torque $\dot{\Omega }=-K{\Omega }^n$ with constant $K$, the displayed combination equals the exponent. A measured value need not be three: magnetospheric torques, evolving fields and timing variability complicate inference.

For $\dot{\Omega }=-K(t){\Omega }^q$, differentiation gives the displayed braking index. Inferring the torque power $q$ from a measured index therefore requires an effectively constant coefficient. Field evolution, changing inclination and moment-of-inertia changes can all contribute to the difference.

For a constant-coefficient braking law $\dot{\Omega }=-K{\Omega }^3$, integration gives $t={\tau }_c[1-(P_{\mathrm{birth}}/P{)}^2]$. The characteristic age approximates the actual elapsed time only when the birth period was short and torque evolution and kinematic timing contributions are negligible.

The mature experiment uses channel-template comparison, pulsar dedispersion, GPS clock-offset transfer, Solar system barycenter corrections and connected pulse counts.^[ref]

For a geometric illustration, let a perfectly uniform source emit pulse $n$ at $t_n=t_0+n{P}_s$, with spin period $P_s$, and let $z_1$ increase away from the observer. At leading order in orbital velocity divided by the speed of light, an extra path length $z_1$ adds ${\Delta }_R=z_1/c$. Consecutive arrivals then obeyA receding source has $v_{\parallel }>0$ and a longer apparent period. This derivation isolates the geometric Doppler effect; a precision relativistic clock model also needs gravitational redshift and time dilation.

The difference between an observed arrival time and its model prediction. Positive pulsar timing residuals mean late arrivals in this convention. A periodic residual can expose an omitted orbit; a smooth long-term residual can also come from a wrong spin model or clock correction. Fitting and checking alternative causes are part of the inference, rather than interpreting every residual as gravitational radiation.

The binary-pulsar Einstein delay is the periodic clock-rate contribution to pulse arrival times from a pulsar moving in an eccentric binary. Near closest approach, the pulsar's orbital speed is larger and the companion's gravitational potential is deeper. Time dilation and gravitational redshift therefore make its clock run at a different rate relative to a distant time coordinate.

At leading order the variable timing term is ${\Delta }_E=\gamma \sin E$, with eccentric anomaly $E$ and time amplitude $\gamma$. Constant clock-rate offsets are absorbed into the fitted pulsar spin. This term is distinct from the geometric binary-pulsar Roemer delay and the propagation Shapiro time delay. Fitting it with other orbital observables constrains the component masses and tests general relativity.^[ref]

A binary pulsar is a pulsar in a binary star: the companion's gravity changes the emitting neutron star's motion and hence its pulse arrival times. The companion need not itself be a detectable pulsar. Pulsar timing turns that motion into a measurable orbit, making the system useful for stellar mass determination and tests of general relativity.

Russell Alan Hulse and Joseph Hooton Taylor Jr. discovered this binary pulsar at Arecibo in 1974; their January 1975 paper reports a roughly $59\mathrm{ms}$ pulse period and $7.75\mathrm{h}$ orbital period. Its unseen companion was not detected as a second radio pulsar. Discovery of the orbit preceded the later binary-pulsar orbital decay test.^[ref]

The 1993 Nobel Prize in Physics recognized the discovery.^[ref]

By 1993, the system was interpreted as a pair of neutron stars.^[ref]

The surprise was a pulsar whose apparent period would not stay fixed. Russell Alan Hulse's Arecibo search used a multichannel receiver, pulsar dedispersion and computer searches over period and pulse width. Follow-up at higher time resolution retained the changing period, excluding inadequate sampling as its explanation. The observed variation exceeded the expected terrestrial Doppler shift; a predicted reversal from decreasing to increasing period appeared on 16 September 1974.^[ref]

The 1975 analysis fitted an eccentric Kepler orbit to Doppler shifts inferred from about 200 five-minute period measurements over ten days. That discovery fit was not yet the later precision arrival-time experiment.^[ref] An orbit explains both signs of the apparent period change without requiring the stellar rotation to alternate between speeding up and slowing down.

The semimajor axis $a_1$ of the pulsar's orbit about the binary barycenter, projected along the line of sight and divided by the speed of light. The orbital inclination $i$ is zero for a face-on orbit. Thus $x$ is a time, often quoted as the numerical number of light-seconds in $a_1\sin i$, not the true binary separation. The relative-orbit semimajor axis is instead $a=a_1+a_2$; conservation of momentum gives $a_1=a{m}_2/(m_1+m_2)$.

Combining Kepler's third law with that barycenter relation gives the timing form of the binary mass function, also valid for eccentric Kepler orbits:Choose $x=2\mathrm{s}$ and $P_b=8\mathrm{h}$ for an original example: $a_1\sin i=599,585\mathrm{km}$ and $f\simeq 0.0773M_{\odot }$, where $M_{\odot }$ is a solar mass. For a specified pulsar mass $m_1=1.4M_{\odot }$, the smallest allowed companion mass occurs at $\sin i=1$ and is approximately $0.698M_{\odot }$. A different assumed $m_1$ changes this bound. One measured $f$ does not determine both masses and orbital inclination.

The geometric light-travel contribution from the pulsar's position within the binary orbit. To construct an eccentric Kepler orbit, use eccentric anomaly $E$, orbital eccentricity $e$, mean anomaly $M=2\pi (t-T_0)/P_b$, and argument of periapsis $\omega$. Here $T_0$ is passage through periapsis. Coordinates in the orbital plane are $X=a_1(\cos E-e)$ and $Y=a_1\sqrt{1-e^2}\sin E$. Rotate through $\omega$ and project by $\sin i$ to obtainThe first equation is Kepler's equation: equal time intervals give equal increments of $M$, not of $E$. Since $1-e\cos E>0$ for $0\le e<1$, each $M$ has a unique $E$ on the corresponding unwrapped orbit. A constant delay is absorbed into the arbitrary arrival-time origin. The circular limit is ${\Delta }_R=x\sin (M+\omega )$, so its amplitude measures $x$. For $e>0$, the delay's offset and amplitude also depend on $\omega$; its maximum is not generally $x$.