Constant gas-radiation ratio vertical disc structure Created 2026-10-07 Updated 2026-11-05
For constant positive and molecular gas constant , the gas–radiation equation of state gives . Vertical hydrostatic equilibrium in harmonic central gravity therefore gives with . The mass density and pressure powers are three and four. Integrating both faces yields and . With constant opacity, radiative diffusion produces a linear flux and height-independent dynamic viscosity . This is an interior formal profile; its low-density surface requires a photospheric treatment beyond diffusion.
For homologous gas-pressure-supported stars, hydrostatic equilibrium gives and . Radiative diffusion gives , so the radius cancels and yields the displayed cubic mass scaling. Composition and opacity are held fixed. Strong radiation pressure alters this gas-pressure relation towards the Eddington luminosity scaling.
Constant-pressure-ratio vertical disk model Created 2026-10-05 Updated 2026-11-05
Assume a thin disk has height-independent opacity and radiation-to-gas pressure ratio , with vertical gravity . Hydrostatic equilibrium and radiative diffusion imply
Balancing its derivative against viscous heating fixes the dynamic viscosity
It is independent of height under these assumptions. For a spherically symmetric potential, . Integrating through the full surface-to-surface thickness gives , which determines the thickness from the radial accretion rate solution.
Convective stability of a constant-opacity irradiated atmosphere Created 2026-10-05 Updated 2026-11-05
A deep irradiated layer can remain stable against convection when its thermal opacity and outward radiative flux are constant. Specify a plane-parallel atmosphere in hydrostatic equilibrium, with constant gravity , positive constant Rosseland mean opacity and upward radiative flux . At depths where incident stellar light has already been absorbed, radiative diffusion gives
Consequently
where is the extrapolated boundary temperature of this simplified diffusion solution. The radiative temperature gradient obeys
For a homogeneous diatomic ideal gas with fixed rotational degrees of freedom, gives adiabatic temperature gradient . The Schwarzschild criterion therefore predicts stability throughout this model. Variable opacity, changing thermodynamics or varying radiative flux can change that conclusion. The diffusion approximation also needs large optical depth; it does not establish stability of the entire optically thin irradiated atmosphere.
Energy transport 2026-10-25
Transport of stellar energy from the interior to escaping radiation, principally by radiative diffusion and convection in ordinary stellar models. Carrying a nonzero luminosity generally requires a nonzero temperature gradient.
In an optically thick steady spherical stellar wind of constant mean molecular weight , let be the outward speed and the isothermal sound speed squared. Mass conservation gives , where a prime denotes differentiation with respect to radius . Use and radiative diffusion to write , where and uses the local Eddington luminosity and opacity . 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 alpha disk dominated by gas pressure, vertical hydrostatic equilibrium gives and . The specified opacity therefore scales as . Viscous heating scales as , while radiative diffusion gives . Their balance gives , hence . Fixed and Keplerian rotation yield the displayed law.
The local thermal stability of an accretion disk follows because, at fixed surface density and radius, cooling has logarithmic temperature slope two while heating has slope one. The viscous stability of an accretion disk follows from on the thermal-equilibrium branch. Both conclusions apply within this pressure, opacity and transport closure.
Gas-pressure viscosity closure with mixed pressure support Created 2026-10-14 Updated 2026-11-05
Let , take constant opacity , and prescribe the dynamic viscosity in a Keplerian disk. The pressure supporting the disk may contain both gas pressure and radiation pressure. Combining viscous dissipation and radiative diffusion gives in the vertical mass coordinate of a disk, where . Set , and . With a radiative-zero disk surface, satisfies , . Writing , the density-weighted mean kinematic viscosity is , giving the displayed formula. This calculation does not assume gas-pressure domination: the closure depends on gas pressure while hydrostatic support can be mixed.
Let and choose . Set , , , and . Vertical hydrostatic equilibrium, viscous dissipation, radiative diffusion and the ideal gas law reduce to , , , . Midplane symmetry gives ; a regular radiative-zero surface has at a universal , with . All physical parameters have disappeared. Uniqueness of this normalized boundary-value problem makes every physical profile a rescaling of one dimensionless solution.
// Target: astrophysics.bigb
An imposed spherical profile gives in hydrostatic equilibrium. For an ideal gas with zero surface pressure, and , where , and . The density has a central cusp and temperature initially rises outward, so positive outward radiative diffusion with regular positive nuclear heating cannot satisfy this imposed profile near the centre.
Main-sequence mass-radius relation Created 2026-10-25 Updated 2026-11-05
An approximate relation between the mass and radius of main-sequence stars in a specified composition and mass range. It is not determined by hydrostatic equilibrium alone. For gas-pressure stellar homology with Kramers' opacity law, radiative diffusion gives . Substitution of gives at fixed composition, explaining precisely which extra radius assumption produces any claimed exponent.
Let for constant mean molecular weight in proton-mass units, and let be the Stefan–Boltzmann constant in the Stefan–Boltzmann law. The gas sound speed used here satisfies , and is the local orbital frequency. For the local opacity approximation , a gas-supported optically thick alpha disk has and . Integrated viscous heating is . Radiative diffusion gives . Equating the fluxes proves the thickness scaling. At fixed composition and , it gives and in Keplerian rotation. Radial powers therefore require a surface-density profile; they are not determined by opacity alone. The formula is local to the stated opacity regime and small disk aspect ratio.
Optically thick medium 2026-10-05
An optically thick medium has optical depth along the relevant path. Photons undergo many interactions before escaping. Radiative diffusion applies under suitable near-isotropy conditions, but a large scattering depth alone does not guarantee thermalization when true absorption is weak.
In local thermodynamic equilibrium, write the one-ray radiative transfer equation as
Large total optical depth removes the boundary contribution, but an additional local condition is essential: the temperature and opacity vary little over an absorption length . Iterating once in that small gradient gives
The integrated Planck function and its derivative are
Using the specified Rosseland mean opacity, with measured here in inverse length, the source's definition therefore gives
The printed leading term lacks a factor under that definition. An isothermal optically thick medium is a concrete counterexample: gives exactly, while the printed formula gives .
A physical radiative flux requires an angular integral. For direction cosine , the same gradient expansion gives . Integration over all directions gives the net radiative diffusion flux
The isotropic equilibrium contribution cancels between opposite directions. For comparison the outward-hemisphere flux alone is at this order. None of these consistent definitions yields the printed mixture of coefficients.
Stellar rotation is a problem of coupled force balance, heat transport and angular momentum transport. A star need not rotate as a rigid body: in the Sun, the convective envelope rotates faster near the equator than near the poles, while much of the radiative interior rotates more nearly uniformly. The tachocline is the thin transition near the base of the convection zone. Surface tracking and Doppler effect measurements constrain the exterior flow; solar rotational mode splitting in helioseismology constrains interior stellar rotation through integrals of the form . Multiple kernels are needed, so surface rotation is not a measurement of every depth.
A useful mean velocity decomposition is , where is distance from the axis and is stellar meridional circulation. In a frame of constant angular velocity , use the relative velocity . Its momentum equation includes
The Coriolis force changes the direction of a parcel's motion but does no work. Thermal driving, release of differential rotational kinetic energy, contraction, or external torques supply the energy of the flows. Solid-body stellar rotation in a suitable hydrostatic equilibrium does not by itself require a meridional flow.
For slowly circulating, approximately axisymmetric material, the leading mechanical balance in an inertial frame is . Taking its curl gives the stellar thermal-wind balance
Here is the axial coordinate, is specific entropy and is the specific heat at fixed pressure. The last form uses approximately radial hydrostatic equilibrium and uniform composition. If the material is barotropic, the right-hand side vanishes and angular velocity is constant on cylinders. Small latitudinal entropy differences can instead support the observed noncylindrical differential rotation. Stresses and inertial terms modify this leading balance. The equation relates shear to baroclinicity; it does not give the circulation speed or explain how the shear is maintained.
That second question requires an angular-momentum budget. With and density-weighted mean quantities, conservation of angular momentum can be written schematically as
The three contributions are Reynolds stress, magnetic Maxwell stress tensor and viscosity transport. Appropriate density weighting replaces the simple correlations in a strongly compressible average. Since in a steady star, the stationary equation gives gyroscopic pumping in a star:
A torque that maintains differential rotation therefore drives a meridional response unless another stress cancels it. Conversely the circulation transports angular momentum and feeds back on the rotation. An assumed meridional pattern without its driving torque is not a closed dynamical model.
In the convection zone, outward heat transport drives convection. Rotation deflects the rising and falling parcels and makes their correlations anisotropic. Nondiffusive convective angular momentum transport can drive shear even from an initially uniform rotation rate; treating every Reynolds stress as a positive eddy viscosity would only smooth shear and miss its source. A limiting parcel that conserves when it moves outwards slows relative to uniform rotation, whereas rapid rotational constraint and correlated horizontal motions change the sign and direction of the net transport. Thus the solar equator-fast pattern requires the actual convective correlations, not a universal parcel argument. The Rossby number compares convective inertia with rotational deflection: deep, slower, larger-scale convection is more rotationally constrained than the fastest photospheric motions. Turnover times range from minutes near the visible surface to weeks or longer at greater depth. On these times the Reynolds stresses can redistribute angular momentum and drive circulation, while the star's total spin changes much more slowly.
In a radiative region, rotation distorts equipotential surfaces. Under uniform rotation and uniform composition, leading hydrostatic equilibrium makes pressure, density and temperature constant on each such surface, but effective gravity varies along it. Radiative diffusion then produces a flux . Its divergence contains , which varies over a surface and generally cannot match a nuclear source constant there. This is the radiative-equilibrium obstruction in a rotating barotropic star. Thermal imbalance drives Eddington-Sweet circulation, with heat advection entering
Its global thermal estimate is
A rotational thermal imbalance of order processes the thermal reservoir on this longer time. For the slowly rotating Sun, is of order years and of order , giving years, much longer than the solar age. Stable composition gradients impede the flow further. Thus classical radiative Eddington-Sweet circulation cannot by itself establish the Sun's interior rotation rapidly. In a rapidly rotating star the same estimate is shorter and rotational mixing can compete with its evolutionary lifetime. Torque-driven convection-zone circulation and radiative thermal circulation have different controlling balances.
Magnetic fields introduce a further coupling. A poloidal magnetic field linking different layers is wound by differential rotation into a toroidal magnetic field, roughly . The resulting Maxwell stress tensor transport angular momentum. Communication along a field occurs on an Alfvén speed crossing time , which can be far shorter than molecular diffusion if a coherent connecting field exists. In steady ideal axisymmetric induction, absence of continuing winding requires ; this is constancy along the field, not automatically solid-body rotation throughout the star. The topology and strength of an interior field matter. Solar dynamo action, especially the interaction of shear with poloidal-field regeneration, couples stellar rotation to the magnetic cycle. Magnetic feedback can change shear and circulation on dynamical or cycle times; an approximately eleven-year activity cycle and twenty-two-year magnetic polarity cycle are distinct from secular spin-down.
A magnetized stellar wind also supplies an external torque. It enforces approximate corotation out to an effective Alfvén radius , giving wind-driven magnetic braking of a solar-type star . Its spin-down time is , much longer than a convective turnover time. For slowly changing moment of inertia and an unsaturated field scaling that gives , integration gives , the Skumanich rotation law. A growing or contracting star also changes and can spin up or down even without torque. Differential spin-down can create shear between layers unless internal transport couples them.
Finally, shear can excite Kelvin-Helmholtz instability, and convection can launch internal gravity waves into the stable radiative interior. Wave damping deposits angular momentum; the net transport depends on excitation and selective absorption and is not fixed by a wave-crossing time alone. Such stresses, magnetic coupling and circulation are candidate agents for the nearly uniform radiative rotation and the narrow tachocline. Molecular kinematic viscosity has a diffusion time generally too long to couple the entire solar interior over its age, whereas a turbulent effective kinematic viscosity can act much faster but is not appropriate to every stable layer. A circulation has an advection time , thermal diffusion has time , and magnetic diffusion has time ; these must be compared with the local stellar rotation and evolutionary times rather than conflated with them. The solar dynamical time is about half an hour, its rotation period is of order a month, and secular magnetic braking takes an evolutionary time. Rapid force adjustment, convective maintenance of shear, cyclic magnetic feedback and slow global spin evolution are separate processes. Their coupled conservation equations and timescale hierarchy are the basis of a consistent account of solar rotation.
Set , so and the prescribed dynamic viscosity is . Keplerian shear and radiative diffusion then give
Differentiating the second expression, equating it to the first and dividing by gives
This reduction uses the gas-pressure viscosity closure, not an assumption that gas pressure dominates the hydrostatic support.
For , introduce the vertical mass coordinate of a disk
Let be the full surface density of a reflection-symmetric disk. The upper surface is at . With , the temperature equation becomes
Write and choose
The equation reduces exactly to . Midplane symmetry gives , and the zero surface temperature gives . The stipulated unique positive nontrivial dimensionless solution therefore applies at every radius and surface density, with all physical dependence contained in and the mass-coordinate scaling. A zero temperature at the ideal surface does not mean zero escaping flux: can have a finite nonzero surface value.
The required viscosity is the density-weighted mean kinematic viscosity, not an unweighted vertical mean. Since ,
Substituting the dimensionless solution gives, with ,
For central mass , , so the requested explicit radius and surface-density dependence is
The universal integral need not be evaluated. No separate solution for is needed to find this mean: the mass-coordinate transformation removes it from both the temperature equation and the viscosity integral. The result is the gas-pressure viscosity closure with mixed pressure support.
In the limiting model with negligible gas pressure, . Its gradient is . Combining radiative diffusion and hydrostatic equilibrium gives
Because and are constant with height at the fixed disk radius,
Equating this to the viscous heating rate forces
This proves constant viscosity in radiation-supported vertical balance. It fixes the local effective dynamic viscosity uniquely, independent of height, surface density and radius within this ideal limiting model. It does not fix a height-independent kinematic viscosity: still depends on the local mass density.
The vertical hydrostatic equilibrium equation balances the gradient of the total pressure against gravity. In a thin disk with Keplerian rotation, expanding the central gravitational field about the midplane gives vertical acceleration , so . Disk self-gravity is neglected in this approximation.
The flux-gradient equation is local thermal balance: the upward radiative flux increases with height because differential rotation dissipates energy. For a Newtonian effective dynamic viscosity , the volumetric heating rate is , which becomes for Keplerian rotation. The symbol here is dynamic viscosity, not kinematic viscosity; the latter is .
The radiative diffusion equation applies to an optically thick medium in local thermodynamic equilibrium. Its diffusion coefficient is , and the equilibrium radiation energy density is . Hence
The minus sign means energy flows down the temperature gradient. Here is the Stefan-Boltzmann constant, the opacity per unit mass, and the speed of light. Constant Thomson scattering opacity is the adopted transport approximation; local thermodynamic equilibrium is an additional assumption behind the radiation equation of state.
Finally the equation of state adds ideal gas gas pressure and isotropic radiation pressure . The Boltzmann constant is , and is the mean mass per gas particle. Both contributions supply the total pressure entering vertical hydrostatic equilibrium.
A Type Ia supernova is a thermonuclear disruption associated with a carbon-oxygen white dwarf in a binary system. Electron degeneracy pressure supports the dense star, and initially depends much less on temperature than an ideal gas pressure does. Once carbon burning heats degenerate material, the immediate pressure response is insufficient to produce the usual expansion-and-cooling thermostat. The strongly temperature-sensitive nuclear reactions can therefore undergo a thermonuclear runaway.
One channel involves accretion toward the Chandrasekhar limit, about solar masses for a typical electron composition. Other proposed channels involve merging white dwarfs or an initial helium-layer detonation that triggers carbon burning below that mass. Thus the Chandrasekhar limit motivates an important model, but it is not an assumption that every Type Ia supernova has exactly the same progenitor mass or ignition history. Burning carbon and oxygen into more tightly bound nuclei releases enough nuclear binding energy to overcome the stellar gravitational binding energy. An ordinary successful event ejects the star rather than leaving the compact core characteristic of a core-collapse supernova.
The ejecta contain iron-group material, including radioactive nickel-56, and intermediate-mass elements from incomplete burning. The radioactively powered Type Ia supernova light curve is powered largely by the nickel-56 to cobalt-56 to iron-56 decay chain: its energetic products heat the ejecta, and radiative diffusion delays the emerging radiation. Near maximum light, ordinary Type Ia supernovae show prominent silicon absorption and lack the hydrogen lines defining Type II events. The spectral classification and the thermonuclear mechanism provide complementary descriptions.
The similar progenitors and radioactive power make Type Ia supernovae useful standard candles, but their peak luminosities are not identical. The relation between light curve width and peak brightness, together with color corrections, allows standardization for distance measurements. Extinction, composition and explosion diversity still introduce systematic uncertainties. A Type Ia event is a white-dwarf thermonuclear explosion with a radioactively powered light curve; its corrected luminosity makes it a valuable distance indicator.
The first equation is vertical hydrostatic equilibrium: the upward pressure gradient balances the central mass's vertical gravitational force. For , the point-mass acceleration is , where Keplerian rotation gives . The disc's self-gravity and vertical inertia are neglected.
The second equation is local thermal balance. The divergence of the vertical radiative flux equals viscous dissipation. In the convention of this alpha disk, , and the dominant Keplerian shear has . Therefore
The radiative flux increases from zero at the symmetric midplane to a positive outgoing value at the upper surface. Radial heat transport, thermal advection and rapid storage of thermal energy are excluded from this local balance.
The third equation is radiative diffusion, equivalently
It applies to optically thick material, with the opacity per unit mass and the Stefan-Boltzmann constant. Here constant electron-scattering opacity represents Thomson scattering; the formula uses a diffusive, nearly thermal radiation field. Positive upward flux requires temperature to decrease upwards.
The final equation is the ideal gas pressure law, with mean molecular weight in proton-mass units. Writing gives . The pressure is gas pressure because radiation pressure is assumed negligible. The four equations express vertical support, local heating balance, diffusive cooling and the gas equation of state.