Stellar astrophysics

Stellar astrophysics studies the structure, formation, energy production, observable properties, and evolution of stars.

Orbital dynamics gives the most direct stellar masses: a resolved binary orbit and distance determine the relative semimajor axis, and Kepler's third law gives the total mass. The barycentric size or radial-velocity amplitude ratio gives individual masses. An eclipsing binary inclination removes the sine-cubed projection factor from spectroscopic masses; a single-lined binary mass function alone does not uniquely determine both masses. Surface-gravity/radius and stellar-track estimates are model-dependent and are calibrated against dynamical binary measurements.

Loss of matter from a star or other astronomical system, for example through a stellar wind or explosive ejection. Mass loss rate conventions must distinguish the positive ejected rate $|\dot{M}|$ from the negative change of the retained mass, $\dot{M}$. In a binary star, the fate of an orbit also depends on the specific angular momentum carried away.

The spin of a star about its axis, which can be solid-body rotation or vary through the interior. Angular momentum transport and rotational mixing can change the stellar composition profile and available fuel. An observed surface rotation rate does not determine the entire interior profile.

In a steady axisymmetric star, a torque due to Reynolds stress, Maxwell stress tensor or viscosity must balance advection of specific angular momentum by stellar meridional circulation. Writing the nonadvective angular momentum transport flux as ${\mathbf{F}}_J$ gives the displayed balance. Thus stresses that maintain differential rotation also drive a meridional response. The stellar thermal-wind balance supplies a separate leading force constraint, not a replacement for this torque equation.

Taking the curl of axisymmetric hydrostatic equilibrium with centrifugal acceleration gives $\varpi {\partial }_z{\Omega }^2=(\nabla p\times \nabla \rho {)}_{\varphi }/{\rho }^2$. For an approximately spherical ideal gas with constant composition, use $\nabla \rho =\nabla p/c_s^2-\rho \nabla S/c_p$ to obtain the displayed relation. Latitudinal entropy differences can support noncylindrical differential rotation. A barotropic fluid gives cylindrical rotation in this leading balance; the relation does not determine the transport stresses or circulation rate.

Outside a spherical star, its gravitational acceleration is $-G{M}_{\star }\mathbf{r}/r^3$. It provides the dominant central force for Kepler orbits in a star-dominated planetary system.

A stellar population is an ensemble of stars described by its formation times, initial masses, composition and present stellar evolution states. Number fractions must be distinguished from mass-weighted fractions. Shared formation times can correlate the states of components of a binary star even when their initial masses have independence.

A Population III star is a first-generation star formed from essentially metal-free primordial gas. Limited cooling generally raises its characteristic fragmentation mass and favors a more top-heavy initial mass function than in enriched star formation.

A Population II star is an old metal-poor star formed from gas enriched by earlier stellar generations. Metal and dust cooling let enriched gas fragment to lower masses than primordial gas under otherwise comparable conditions.

Integrated spectra can confuse older stellar populations with more metal-rich ones because both can redden the light and change absorption lines. Line strengths depend on stellar metallicity, temperatures, ages, and abundance ratios. Weak metal lines in a young population dominated by hot bright stars need not indicate low stellar metallicity.

Given a birth-mass distribution and age distribution, an evolutionary state occupies the mass-age set where its starting lifetime has elapsed and its ending lifetime has not. The state number fraction is the probability measure of this set. Initial mass function tail coordinates and the probability integral transform can turn these weighted regions into simple geometric areas.

A phase bounded by mass-dependent lifetimes ${\tau }_a(M)$ and ${\tau }_b(M)$ occupies ${\tau }_a(M)<t<{\tau }_b(M)$, clipped to the population age range. Its boundaries can be mapped into uniform probability coordinates to evaluate evolutionary-state fractions.

If a population forms at a constant number rate for duration $T$ and every birth object is retained in the census, its present ages are uniform on $[0,T]$. A fixed initial mass function with finite mean birth mass also makes constant mass-based star formation rate equivalent to constant number rate.

A stellar atmosphere is the outer region that determines the escaping spectrum and matches the interior stellar structure to its radiation boundary conditions. A photosphere is the layer near unit optical depth, not a zero-temperature boundary. The inward intensity at the outer surface vanishes when there is no external illumination.

A stellar abundance can be expressed as an elemental or isotopic mass fraction, or as a number ratio to hydrogen. Inferring it from a spectral line requires a model for level populations, opacity, radiation and velocity fields. Agreement between suitable ionization stages and weak photospheric lines tests whether the adopted stellar atmosphere is adequate.

The abundance of chemical elements heavier than helium in astronomical material. A metal mass fraction measures their combined mass; spectroscopic abundance ratios are other measures of metallicity and need not be numerically identical to that fraction.

The fraction of total mass in chemical elements heavier than helium, so the hydrogen mass fraction, helium mass fraction and metal mass fraction satisfy $X+Y+Z=1$. Despite the astronomical name, this includes elements such as carbon and oxygen that are not chemical metals.

When radiative rates compete with or dominate collisions, level populations need statistical equilibrium coupled to radiative transfer, rather than thermal population laws at the local temperature. Radiation from other depths can make the local radiation field differ from the Planck function. This changes ionization and spectral-line strengths, so a geometrically thin atmosphere may still require non-LTE abundance modelling.

A heated layer above a photosphere can have an outward temperature increase. Strong lines formed there may show emission where an outward-cooling photosphere would produce absorption. The energy maintaining the temperature inversion is additional atmospheric heating; the exact height and transfer regime vary between stars.

A narrow spectral feature records enhanced emission or absorption near a transition frequency. In a stellar atmosphere, its strength depends on atomic populations, opacity, velocity fields and the source function. Lines therefore diagnose composition as well as the physical atmosphere in which they form.

Collisions and fluctuating neighboring-particle electric fields alter atomic energy levels and coherence times, broadening a spectral line. In a simple impact treatment the coherence-decay rate gives a Lorentzian width of order $\Gamma /(2\pi )$; the precise definition of full or half width must be stated. Unlike Doppler broadening, this depends primarily on collision rates and perturber density. In stellar opacity calculations it alters line overlap and the low-opacity windows entering the Rosseland mean opacity.

A Voigt profile is the convolution of a Gaussian distribution and a Cauchy distribution. For an absorption line, thermal Doppler broadening contributes the Gaussian core while finite excited-state lifetime contributes the Lorentzian wings. Instrumental broadening and unresolved motions must also be modeled when inferring column density and temperature.

An emission line is a localized excess of specific intensity above the continuum. Heated chromospheric material can give a line source above the photospheric continuum source. Optically thin circumstellar gas can also add line photons through recombination or collisional excitation. Such line formation need not obey a thermal photospheric source prescription.

An absorption line is a deficit of observed spectral specific intensity relative to the surrounding continuum. In a photosphere with temperature falling outward, increased line opacity makes photons emerge from cooler layers than the continuum. For an absorption-dominated LTE source, the Planck function then supplies a smaller line specific intensity. A foreground layer with source below the incident continuum gives another direct transfer description.

A low-mass star maintained by sustained hydrogen burning, rather than the continuing cooling of a brown dwarf. Near the hydrogen-burning minimum mass its radius is of order a Jovian radius. Stars below roughly $0.35M_{\odot }$ are approximately fully convective, with a radiative photosphere; mass, composition and age affect the precise transition.

The solar radius is the conventional unit of length based on the radius of the Sun. The nominal value is $R_{\odot }=6.957\times 10^8\mathrm{m}$.

The solar corona is the Sun's hot, dilute outer atmosphere. Its optically thin atomic lines diagnose temperature, density and composition through collisional-radiative models and ionization balance. Its plasma extends into the solar wind; the coronal approximation describes useful low-density equilibrium limits rather than every coronal condition.

A larger-scale horizontal cellular flow pattern on the Sun's surface. Transport by these flows helps arrange small magnetic elements into a solar magnetic network at cell boundaries. The observed pattern does not establish a unique thermal driving depth or an exact hierarchy of isolated steady eddies.

A solar granule is a photospheric convection cell, with relatively hot rising material in its interior and cooler descending material near its edges. The pattern formed by their bright interiors and darker edges is solar granulation. Their evolving horizontal outflows can sweep weak magnetic flux toward surrounding lanes.

The evolving cellular intensity pattern on the Sun's photosphere, produced by solar granules. Its radiative appearance and flow pattern are related but are not identical observables. Near-surface horizontal flow can redistribute weak magnetic field into the surrounding intergranular lanes.

An intergranular lane is a region between solar granules, commonly associated with cooler material and descending flow. Horizontal convergence can concentrate approximately vertical magnetic flux there. The field's geometry, diffusivity and feedback determine the actual concentration rather than the lane label alone.

A solar neutrino is a neutrino produced by a nuclear reaction in the Sun, principally through the proton–proton chain and CNO cycle. Its weak interactions allow it to escape directly from the stellar core. Accounting for neutrino oscillation allows measurements of the different detected flavors to test the nuclear reaction power and rates independently of escaping photons.

The solar tachocline is the rotational transition near the base of the convection zone, approximately at $0.7$ solar radii, between differential rotation above and more nearly uniform rotation in the radiative interior. Its shear can contribute to the Omega effect, and stable stratification can help store strong toroidal magnetic fields. It is a proposed site for important parts of the solar dynamo, rather than a demonstrated unique location for all field generation.

The solar magnetic field is generated and maintained by electrically conducting flows through solar dynamo action. It includes a large-scale approximately dipolar component near activity minimum and strong, intermittent fields in active regions and sunspots. Its organized polarity reversals define the magnetic solar cycle. Magnetic stresses modify the generating flows, so a saturated dynamo requires nonlinear feedback as well as magnetic induction.

A spatially extended pattern of magnetic concentrations near boundaries of supergranulation. Surface transport and convergence redistribute magnetic elements into this network; flux emergence, cancellation and diffusion also continually modify it. It is distinct from the much smaller-scale field pattern in individual intergranular lanes and from the organized strong fields of sunspots.

The solar dynamo is the induction and regeneration of the solar magnetic field by conducting flows powered ultimately by solar convection and rotation. In a mean-field dynamo, differential rotation supplies the Omega effect and helical correlations or the Babcock-Leighton mechanism regenerate a poloidal magnetic field. The mean induction equation iswhere $\mathcal{E}=⟨\mathbf{u}\times \mathbf{b}⟩$ is the mean-field electromotive force. A model must close this correlation and include feedback to predict a saturated magnetic cycle; prescribing a shear alone does not close the regeneration loop.

Rotation influences convection, its helicity and differential rotation, which feed dynamo action. The generated magnetic field produces surface activity and strengthens the angular-momentum lever arm of a stellar wind. Wind braking lowers the rotation rate and ordinarily reduces unsaturated activity. Lorentz force feedback also changes the differential rotation and convective motions, saturating the dynamo on shorter time scales. Secular evolution, cyclic polarity reversal and irregular modulation are distinct time scales of this coupled system.

An interface solar dynamo separates the region of strong differential rotation from the region of poloidal-field regeneration, often placing shear near the tachocline and an alpha effect in the convection zone. Magnetic transport across the interface couples the two regions. The separation can provide toroidal-flux storage in a relatively stable layer, but growth, period and migration depend on the transport, diffusivity contrast, and signs of the regeneration and shear coefficients.

An ideal planar alpha-Omega dynamo places poloidal regeneration above an interface and shear-generated toroidal induction below it. With equal constant diffusivities, decay away from the interface and ordinary electromagnetic matching, its surface-wave branch grows when $D>32$ for positive $k,\alpha ,V$. The model separates the two regeneration regions while allowing diffusive coupling. The threshold refers to this model, not every solar interface dynamo or a full isotropic alpha closure.

Let $p=s-i\omega$ and $q^2=k^2+p/\eta$, with $\mathrm{Re}q>0$. The upper poloidal amplitude is $(A_0+Cz)e^{-qz}$ and the lower toroidal amplitude is $(B_0+cz)e^{qz}$. Their forcing equations give $2\eta qC=\alpha B_0$ and $-2\eta qc=ikV{A}_0$; interface matching gives $C=2q{A}_0$ and $c=-2q{B}_0$. Eliminating amplitudes yields the dispersion relation. Its potentially growing branch has $s=\eta k^2(\sqrt{D/32}-1)$ and $\omega =-\eta k^2\sqrt{D/32}$ with the convention $e^{ikx-i\omega t}$. The opposite branch decays.

For a planar interface, equal permeability and no singular sheet current make all magnetic components continuous. A nonzero horizontal Fourier wavenumber then gives the first three conditions in the title. Continuity of the tangential electric field supplies the last: continuous mean velocity cancels its advective contribution, and the retained alpha mean-field electromotive force is along the toroidal direction, so the other tangential electric component is $-U{\partial }_{x}A-\eta {\partial }_{z}B$. The given alpha-Omega equations retain only toroidal-to-poloidal alpha generation; adding a full isotropic alpha contribution would alter both the bulk equation and the matching problem.

A flux-transport solar dynamo couples toroidal and poloidal field generation through meridional advection, turbulent magnetic diffusion and possibly turbulent pumping. A common model uses the surface Babcock-Leighton mechanism, differential rotation near the tachocline, poleward surface circulation and an equatorward deep return flow. If advection dominates, circulation strongly influences the cycle period and migration; if diffusion dominates, diffusive coupling is important instead. Agreement with a butterfly pattern depends on the assumed transport geometry and coefficients, not merely on including a circulation term.

The Babcock-Leighton mechanism regenerates large-scale poloidal solar flux through the emergence and dispersal of systematically tilted bipolar active regions. With Joy's law tilt, opposite polarities have different latitudes: cancellation near the equator and poleward transport of surviving following-polarity flux alter and reverse the polar field. Thus a toroidal field that produces rising bipolar regions feeds back into a poloidal magnetic field. The source is intrinsically nonlocal between the interior flux reservoir and surface emergence and is stochastic because active-region fluxes and tilts vary.

A sunspot is a relatively cool, dark photospheric region associated with concentrated magnetic flux. The strong magnetic field inhibits convective heat transport, reducing the emergent local intensity. Spots commonly emerge in bipolar active regions, whose polarity and systematic tilt follow Hale's polarity law and Joy's law. Their number and latitude distribution are visible tracers of the solar cycle.

The axis of a bipolar solar active region is statistically tilted relative to the local east-west direction, with its following polarity typically farther from the equator and a mean tilt that grows with latitude. Scatter is substantial. The separation in latitude of the two polarities allows their surface dispersal to generate a net poloidal field through the Babcock-Leighton mechanism.

The solar cycle is the recurrent variation of magnetic activity, with a characteristic sunspot-number period of about eleven years. The large-scale magnetic polarity reverses from one activity cycle to the next, giving an approximately twenty-two-year magnetic cycle. Activity cycles vary in length and amplitude. A satisfactory solar dynamo model must account for polarity reversal, equatorward activity-belt migration in the solar butterfly diagram, the phase relation between polar and active-region fields, and long-term modulation including solar grand minima.

A solar grand minimum is a prolonged interval, lasting several ordinary cycles, of strongly reduced sunspot activity. Historical records identify the Maunder minimum; cosmogenic carbon-14 and beryllium-10 provide longer indirect records because stronger heliospheric magnetic shielding generally reduces cosmic-ray production of these isotopes. Weak activity does not imply a zero solar magnetic field. Proposed causes include stochastic weakening of poloidal-field generation, nonlinear dynamo modulation, and intermittency near an excitation threshold.

The Maunder minimum is the interval approximately 1645–1715 when recorded sunspot activity was exceptionally low. It is the best-known historical solar grand minimum. Historical coverage and cosmogenic-isotope records require care in inferring its precise magnetic amplitude; evidence for weak modulation need not imply that the normal spot-producing dynamo disappeared completely.

Bipolar solar active regions preferentially have one leading magnetic polarity in a given hemisphere and cycle, the opposite leading polarity in the other hemisphere, and reversed assignments in the next cycle. These statistical polarity rules reveal the ordered toroidal component of the solar magnetic field. Individual active regions can violate the preferred pattern, so the law is a statistical regularity rather than an exact rule for every spot pair.

A solar butterfly diagram plots sunspot occurrence or magnetic flux against latitude and time. New-cycle sunspot belts appear at middle latitudes and migrate toward the equator, producing successive wing-shaped patterns. Magnetic butterfly diagrams additionally show poleward flux transport and polar-field reversals. The diagram constrains both Parker dynamo wave propagation and magnetic transport in flux-transport solar dynamos.

The opacity mechanism drives stellar pulsation when compression increases opacity in a layer, trapping heat at a phase that feeds mechanical oscillation. Helium ionization zones provide this driving in Cepheid variables.

A Cepheid variable is a pulsating star whose period is tightly related to its luminosity, allowing its distance to be inferred from its measured period and apparent magnitude.

The pulsation period of a Cepheid variable is correlated with its absolute magnitude in a specified photometric band. Calibrating the relation with stars of known distance allows an observed period and corrected apparent magnitude to give a distance. Extinction and metallicity corrections matter for a precise cosmic distance ladder.

A supernova is a luminous stellar explosion produced by thermonuclear disruption or the collapse of a stellar core.

A hypothetical explosion in which a degenerate carbon-oxygen core within a single evolved star grows toward the Chandrasekhar mass and ignites explosive carbon burning while retaining an envelope. This differs from accretion onto an already exposed white dwarf in standard Type Ia supernova scenarios. Evolutionary calculations by Doherty et al. explain the strong dependence of this proposed channel on envelope mass loss; their adopted winds remove the carbon-oxygen envelope before this outcome. It is a proposed possibility under different assumptions, rather than an established endpoint for normal $5M_{\odot }$ stars.

Collapse of a massive stellar core can launch an explosion that ejects material processed by hydrostatic and explosive burning. Whether an explosion succeeds, and how much material escapes instead of falling back, depend on the core and envelope structure. A failed supernova can retain much of the nucleosynthetic core; initial stellar mass alone is not a complete prediction.

Core collapse can fail to produce an ordinary successful ejecting supernova, leaving a black hole and retaining much of the processed core. Some envelope material can still escape, depending on the event; failure is not equivalent to a universal zero-yield prescription.

In a hot radiation-supported stellar core, electron-positron thermal pair abundance can grow as radiation energy is converted into pair rest energy. The stellar adiabatic exponent can fall below $4/3$, permitting contraction and rapid heating. If the resulting explosive nuclear fusion releases enough energy, it disrupts the star in a pair-instability supernova. Weaker events can produce pulses rather than full disruption, and still heavier cores can collapse without an explosion. Thus the thermal relativity scale alone neither establishes the instability nor predicts its outcome.

A Type Ia supernova is a thermonuclear supernova whose standardized luminosity makes it a useful astronomical distance indicator.

The nickel-56 to cobalt-56 to iron-56 decay chain supplies gamma rays and charged particles that heat supernova ejecta. Early expansion and photon diffusion delay the escaping optical light; after a peak on a scale of weeks, the luminosity falls as decay power declines and gamma rays escape more readily. With comparable diffusion times and deposition efficiencies, peak luminosity scales approximately with the nickel-56 mass. Using one fixed absolute luminosity without accounting for composition-dependent yields can therefore bias standard-candle distances.

Approximate a fixed mass of iron-group freeze-out ashes by nickel-56 and iron-54, with no further weak neutronization and negligible residual light particles. Their electron fractions are $1/2$ and $26/54$. Conservation of the initial electron fraction in stellar matter, $1/2-X_{22}/22$, then gives the displayed nickel mass fraction. Its argument is the initial neon-22 mass fraction, which need not equal the earlier nitrogen metallicity. The two-isotope mass-budget approximation is not an exact full nuclear equilibrium with all free-particle densities set to zero.

A brown dwarf is a self-gravitating object too light to sustain ordinary hydrogen fusion. Electron degeneracy limits its contraction, leaving radii near that of Jupiter across a broad mass range.

The conventional deuterium-burning boundary is near $13$ Jupiter masses, with composition-dependent variation. It is a useful classification scale rather than a sharp formation-based distinction between planets and brown dwarfs.

The hydrogen-burning minimum mass is the lowest mass that can sustain core hydrogen fusion, about $0.075$--$0.08$ solar masses for near-solar composition. Less massive brown dwarfs cool continually after any deuterium burning ends.

A fixed stellar homology profile with $P=\rho \mathcal{R}T/{\mu }_i+K_{\mathrm{nr}}{\rho }^{5/3}$ giveswhere $A,B>0$. As a function of $1/R$ this is a concave quadratic; its maximum is $A^2{G}^2{\mu }_i{M}^{4/3}/(4B{K}_{\mathrm{nr}}\mathcal{R})$, at $R=2B{K}_{\mathrm{nr}}{M}^{-1/3}/(AG)$. At the maximum, the central electron degeneracy pressure and thermal pressure are equal. If this temperature is insufficient for sustained hydrogen burning, the object becomes a brown dwarf. Only the positive-temperature segment represents a physical sequence.

For a homologous stellar polytrope of index $3/2$, take $P=K{\rho }^{5/3}+\mathcal{R}\rho T/\mu$ and zero surface pressure. The stellar virial theorem, with $\Omega =-6G{M}^2/(7R)$, givesThe normalized density profile fixes $A$ independently of mass and radius at fixed composition. Writing ${\overline{T}}_M=\alpha /R-\delta /R^2$ shows that its maximum is ${\alpha }^2/(4\delta )=3\mu G^2{M}^{4/3}/(49\mathcal{R}A)$, at $R=7A{M}^{-1/3}/(3G)$. Below this radius contraction increases the fraction of support provided by electron degeneracy pressure and lowers the thermal temperature. A fixed normalized profile also fixes the ratio of central to mass-averaged temperature, giving a hydrogen-burning minimum mass.

Star formation is the gravitational collapse and fragmentation of cold interstellar gas into protostars and eventually stars. Cooling, turbulence, magnetic fields, rotation, and feedback determine its efficiency and characteristic masses.

A protostar is a forming star whose energy comes chiefly from gravitational contraction and accretion, before sustained core hydrogen burning supplies its luminosity. A fixed-mass contraction model approximates the stage after most accretion has ended; an actively accreting protostar need not have fixed mass or the same Hayashi track.

Gravitational fragmentation divides a self-gravitating gas distribution into condensations that can undergo separate collapse. Pressure support, cooling and rotation influence the scales and outcome; it is not equivalent to an instantaneous equal-mass fission of a rigidly rotating sphere. Isothermal fragmentation is a useful limiting case in which thermal adjustment keeps the gas temperature approximately fixed.

The Kennicutt–Schmidt law relates star-formation-rate surface density to gas surface density, commonly as ${\Sigma }_{\mathrm{SFR}}=A{\Sigma }_g^N$. Disk-averaged total-gas measurements give $N$ near $1.4$, while resolved molecular-gas relations are often closer to linear.

The gas-depletion time is the time required to consume a gas reservoir at its current star-formation rate, $t_{\mathrm{dep}}=M_g/{\dot{M}}_{\ast }$. It ignores replenishment, returned stellar mass, and changes in the rate.

The star-formation efficiency per free-fall time is the fraction of gas converted into stars during one gravitational free-fall time, ${ϵ}_{\mathrm{ff}}=t_{\mathrm{ff}}/t_{\mathrm{dep}}$. Molecular clouds commonly have values of order one per cent.

The Jeans mass is the characteristic gas mass above which self-gravity overcomes pressure support. Up to the convention used to define the collapsing region,for an ideal gas of fixed composition.

During isothermal collapse, $T$ is constant and the Jeans mass decreases as ${\rho }^{-1/2}$. Successively denser regions can therefore become independently unstable and fragment to smaller masses.

For adiabatic monatomic gas, $T\propto {\rho }^{2/3}$ and the Jeans mass increases as ${\rho }^{1/2}$. Once cooling fails and collapse becomes adiabatic, further fragmentation is suppressed.

The initial mass function is the distribution of stellar birth masses. It records the mass scales selected by cloud fragmentation, accretion, feedback, and stellar multiplicity.

For an initial mass function proportional to $M^{-p}$ on $M\ge M_{\min }$ with $p>1$, integration and normalization give the displayed upper-tail fraction. The probability integral transform makes $X$ uniform on $[0,1]$ for a randomly sampled birth mass. Evolutionary lifetime regions can then be measured as ordinary areas when the independent age coordinate is also uniform.

Stellar structure follows from mass conservation, hydrostatic equilibrium, energy generation, and energy transport together with an equation of state and opacity law.

The temperature at a star's centre. In gas-pressure-supported stellar homology, the stellar hydrostatic equation and ideal gas law give $T_c\propto \mu GM/(\mathcal{R}R)$, with a coefficient set by the dimensionless density and pressure profiles. It controls the strongly temperature-dependent stellar energy-generation rate.

A stellar core is the central region distinguished from the surrounding envelope by its composition, energy generation or thermodynamic state. A core boundary at radius $R_c$ encloses mass $M_c=m(R_c)=4\pi {\int }_0^{R_c}\rho (r)r^{2}dr$. In a hydrogen-shell-burning red giant, a helium-rich core can be inert while a narrow surrounding layer supplies the luminosity. Core mass and core radius are separate variables: the relation between them must come from a structural model, such as a polytropic mass-radius relation, rather than from the word core alone.

Uniform solid-body rotation adds the centrifugal potential $-{\Omega }^2{\varpi }^2/2$ to the gravitational potential. In the corotating frame, hydrostatic equilibrium becomes $\nabla P=-\rho \nabla \varphi$. The Poisson equation gives ${\nabla }^2\varphi =4\pi G\rho -2{\Omega }^2$. On regular connected equipotential surfaces, $P$ is constant; taking the curl shows that $\rho$ is constant too. Uniform composition and an invertible equation of state then make $T$ constant on each such surface. These statements do not require $|\nabla \varphi |$ to be constant there.

This approximation represents the gravitational potential by a central point mass while retaining the centrifugal potential. A surface through polar radius $R_p$ satisfies $R_p/r+q(r/R_p{)}^2{\sin }^2\theta /2=1$, where $q={\Omega }^2{R}_p^3/(GM)$. The inner closed surface is regular below $q=8/27$; at the limiting value its equatorial radius is $3R_p/2$ and equatorial effective gravity vanishes. The approximation is useful for illustrating oblateness and gravity darkening, rather than modeling the density distribution of every rotating star.

The divergence theorem and Poisson equation give ${\int }_S\nabla \varphi \cdot d\mathbf{S}=4\pi Gm-2{\Omega }^{2}V$ on a bounding equipotential surface. Equipotential luminosity conservation with stellar circulation then fixes $f=L/(4\pi Gm-2{\Omega }^{2}V)$. Combining $dT/d\varphi =-f/\chi$, $dP/d\varphi =-\rho$ and radiative conductivity gives the displayed logarithmic temperature-pressure gradient. Both appearances of mass mean the enclosed mass $m$; only on the outer boundary can both be replaced by the total mass.

Slow circulation in planes containing the rotation axis can transport entropy through a rotating star while the leading mechanical balance remains nearly hydrostatic equilibrium. Its steady mass flux satisfies the continuity equation. It is distinct from assuming an exactly motionless radiative interior, and from turbulent convection.

The thermal imbalance of a slowly rotating radiative star can drive stellar meridional circulation. Its fractional centrifugal distortion is ${ϵ}_{\Omega }\sim {\Omega }^2{R}^3/(GM)$, so balancing an order-${ϵ}_{\Omega }$ fraction of the luminosity by thermal advection estimates $t_{\mathrm{ES}}\sim t_{\mathrm{KH}}/{ϵ}_{\Omega }$, where $t_{\mathrm{KH}}\sim G{M}^2/(RL)$. Stable composition gradients can further impede the circulation. This radiative thermal estimate is different from rapid torque-driven circulation in a convective envelope.

For steady stellar meridional circulation, $T\nabla s=\nabla (h+\varphi )$ and $\nabla \cdot (\rho \mathbf{v})=0$. Thus $\rho \mathbf{v}\cdot T\nabla s=\nabla \cdot [\rho \mathbf{v}(h+\varphi )]$. In the leading uniformly rotating stellar hydrostatic equilibrium, $h+\varphi$ is constant on an equipotential surface; its total advected flux vanishes because the net mass flux vanishes. Integrating the energy equation then identifies the outward luminosity with enclosed nuclear energy generation, even though local heat advection is nonzero.

On a regular equipotential surface, $f,f^{\prime },\rho$ and ${\nabla }^2\varphi$ are constant, but $|\nabla \varphi {|}^2$ generally varies. Consequently a flux $\mathbf{F}=f(\varphi )\nabla \varphi$ cannot ordinarily satisfy local radiative equilibrium $\nabla \cdot \mathbf{F}=\rho ϵ$ at all latitudes. Special cases such as $f^{\prime }=0$ or spherical symmetry evade this obstruction. Stellar meridional circulation supplies heat advection in the leading hydrostatic treatment.