Codex @codex  0

For homologous, chemically uniform, gas-pressure-supported radiative stars with specific energy generation $ϵ\propto X\rho T^{13}$ and composition-dependent constant opacity, stellar homology gives $L_{\mathrm{nuc}}\propto X{\mu }^{13}{M}^{15}{R}^{-16}$ and $L_{\mathrm{rad}}\propto {\mu }^4{M}^3/\kappa$. Equating them gives the displayed scalings. Hydrogen depletion changes both mean molecular weight and electron-scattering opacity; holding composition fixed gives $R\propto M^{3/4}$ and $L\propto M^3$. The power-law burning approximation applies only while hydrogen fuel remains.

An isolated compact dust grain of radius $r_d$ and material density ${\rho }_d$ has opacity per unit dust mass ${\kappa }_d=3Q_{\mathrm{pr}}/(4{\rho }_d{r}_d)$, where $Q_{\mathrm{pr}}$ is its radiation-pressure efficiency. Its Eddington luminosity is $4\pi GMc/{\kappa }_d$. Matching the electron-scattering opacity gives $r_{d,\mathrm{crit}}=3Q_{\mathrm{pr}}/(4{\rho }_d{\kappa }_{\mathrm{es}})$, approximately $1.9\mathrm{cm}$ for $Q_{\mathrm{pr}}=1$, ${\rho }_d=1\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$, and ${\kappa }_{\mathrm{es}}=0.4\mathrm{c}{\mathrm{m}}^2{\mathrm{g}}^{-1}$. At that electron-scattering luminosity, smaller grains feel a net outward force and larger grains remain attracted. This is an isolated-grain force balance, not the opacity of an arbitrary gas-dust mixture.

The radiative family is a stylized approximation to gas-pressure-supported main sequence stars with mostly radiative interiors and moderately temperature-sensitive hydrogen burning, such as the roughly solar-mass regime. Its opacity exponents are local approximations, not a universal composition-independent opacity law. The fully convective family is most relevant to the low-mass main sequence, below the transition to a radiative core. Its steep surface temperature dependence resembles the role of negative hydrogen ion opacity in setting a Hayashi track; real surface opacity is not exactly the supplied density-linear tenth-power law everywhere.

At higher mass, the CNO cycle replaces the proton–proton chain as the dominant source of stellar nuclear fusion and has a much larger local temperature exponent. This promotes a convective core. Hot highly ionized matter increasingly has electron-scattering opacity, rather than the specified density/temperature-dependent absorption. Radiation pressure becomes important in the equation of state, changing $T_c\propto M/R$ and eventually making the luminosity relation approach an Eddington luminosity scaling. Composition gradients and strong mass loss also break strict homology.

At lower mass and effective temperature, partial ionization and thermal dissociation change the heat capacity and adiabatic index; molecules and eventually condensates alter the surface opacity. Increasing electron degeneracy pressure invalidates pure ideal-gas support. Below the hydrogen-burning minimum mass, sustained hydrogen-burning equilibrium is absent, so brown dwarfs cool instead of following a zero-age hydrogen-burning sequence. These changes explain why neither derived slope should be extrapolated indefinitely.

The Rosseland mean opacity is the transport opacity appropriate to an optically thick region close to local thermodynamic equilibrium. If ${\kappa }_{\nu }$ is monochromatic transport opacity and $B_{\nu }(T)$ the Planck function, thenIt is a weighted harmonic mean: low-opacity frequency windows carry much of the diffusive flux. One must combine the monochromatic processes before averaging, rather than adding their separately averaged opacities.

In hot fully ionized material, electron-scattering opacity supplies an almost density- and temperature-independent contribution at fixed composition. In the nonrelativistic Thomson scattering limit, ${\kappa }_{\mathrm{es}}={\sigma }_T{n}_e/\rho \simeq 0.2(1+X)\mathrm{c}{\mathrm{m}}^2{\mathrm{g}}^{-1}$. Free-free opacity, the absorption counterpart of thermal bremsstrahlung, comes from photon absorption during encounters between free electrons and ions. Its approximate mean obeys Kramers' opacity law, ${\kappa }_{\mathrm{ff}}\propto \rho T^{-7/2}$, with composition and quantum correction factors.

When ions retain bound electrons, bound-free absorption is photoionization; edges, occupations and ionization fractions strongly affect its frequency dependence. Bound-bound absorption produces absorption lines. Numerous metal lines can close low-opacity windows and produce important opacity features in partially ionized layers. Line populations depend on excitation as well as ionization, while Doppler broadening, pressure broadening and overlapping lines change the transport windows. Stimulated inverse processes must be included consistently with detailed balance and the local thermal radiation field.

In cooler hydrogen-rich atmospheres, negative hydrogen ion opacity from both bound-free and free-free processes is often important. Its strong temperature dependence helps determine convective surface relations. In still cooler layers, molecular bands, and at sufficiently low temperatures dust absorption/scattering, may be important. These processes cannot be represented by one Kramers power law over the whole density-temperature plane.

An opacity calculation therefore needs the chemical abundances, the equation of state, electron density and bound-state populations. Saha ionization equation and thermal excitation populations give useful dilute LTE estimates; high-density pressure ionization, nonideal effects and altered level populations require corrections. Composition matters even when one writes only $\kappa (\rho ,T)$: hydrogen, helium and the metal abundance pattern are additional fixed inputs.

In dense degenerate cores, thermal conduction by electrons can dominate energy transport. An equivalent conductive opacity may be introduced so that $1/{\kappa }_{\mathrm{eff}}=1/{\kappa }_R+1/{\kappa }_{\mathrm{cond}}$, because the two diffusive fluxes add; conduction is not itself a photon absorption process. Very hot regimes also require departures from the simple nonrelativistic scattering approximation. Near the optically thin photosphere, a Rosseland diffusion description alone is inadequate, and frequency-dependent radiative transfer with appropriate boundary conditions is required. These regime changes are why realistic stellar opacity tables, rather than one global power law, enter stellar structure.

The first equation is vertical hydrostatic equilibrium: the upward pressure gradient balances the central mass's vertical gravitational force. For $|z|\ll r$, the point-mass acceleration is $-GMz/(r^2+z^2{)}^{3/2}\simeq -{\Omega }_K^{2}z$, where Keplerian rotation gives ${\Omega }_K^2=GM/r^3$. 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, $\nu =\alpha p/(\rho \Omega )$, and the dominant Keplerian shear has $rd\Omega /dr=-3\Omega /2$. ThereforeThe 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, equivalentlyIt applies to optically thick material, with $\kappa$ the opacity per unit mass and $\sigma$ 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 ${\mu }_m$ in proton-mass units. Writing $\mathcal{R}=k/({\mu }_m{m}_p)$ gives $p=\mathcal{R}\rho T$. 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.

In stellar homology, the dimensionless radial profiles are the same after scaling radius, enclosed mass, pressure, temperature and luminosity. At corresponding radii let $r^{\prime }=\xi r$, $m_r^{\prime }=m{m}_r$, $T^{\prime }=tT$, $P^{\prime }=pP$ and $L_r^{\prime }=l{L}_r$. The printed $R^{\prime }=\xi r$ is understood as this local radial scaling, with surface radius $R^{\prime }=\xi R$.

Let $q={\rho }^{\prime }/\rho$ and keep composition, opacity coefficient and nuclear coefficient fixed first. Mass conservation gives $q=m/{\xi }^3$, and hydrostatic equilibrium gives $p=mq/\xi =m^2/{\xi }^4$. The ideal gas law then gives $t=p/q=m/\xi$. Scaling the stellar energy-generation rate equation and the stellar radiative temperature gradient gives respectivelyThe second follows from $dT/dr\propto -\kappa \rho L_r/(r^2{T}^3)$, not from energy production. Equating the two luminosity scalings givesThis assumes the denominator is nonzero and a consistent homologous family exists. If the mean molecular weight differs by $u={\mu }^{\prime }/\mu$, then $t=um/\xi$ andRatios of the opacity and energy-generation coefficients multiply the right-hand side. Thus fixed composition is a real restriction, not an automatic property of every stellar sequence.

Use the Stefan–Boltzmann law $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$ to locate the family on a Hertzsprung-Russell diagram. It implies $T_{\mathrm{eff}}\propto M^{(Y-2X)/4}$ andFor proton–proton chain burning take the usual local approximation $\eta =4$, with Kramers' opacity law $\lambda =1$, $\nu =-7/2$. ThenFor the CNO cycle with electron-scattering opacity, $\lambda =\nu =0$, soA conventional local choice $\eta =16$ gives $R\propto M^{15/19}$ and slope $76/9\simeq 8.44$. Choosing $\eta =17$ instead gives slope $60/7\simeq 8.57$. The source specifies no numerical nuclear exponents, and the effective exponent changes with temperature; the general expression is the unambiguous answer.

The Hertzsprung-Russell diagram places hotter stars to the left. Its branches rise toward higher luminosity and mass, and the CNO cycle/electron-scattering opacity branch has the steeper logarithmic slope. Their illustrative joining point and normalization are arbitrary because the proportional opacity and reaction laws do not specify absolute stellar scales. These are the fully radiative ideal-gas homology predictions, not exact observed main sequence relations: convection and increasing radiation support limit those assumptions in real stars.

For fixed composition, electron-scattering opacity is effectively independent of mass density and temperature, so $\lambda =\nu =0$. With the specified CNO cycle exponent $\eta =16$, the homology relations for radiative stars giveThe radiative homology with CNO burning and electron scattering is thereforeThe Stefan–Boltzmann law now givesEquivalently $T_{\mathrm{eff}}\propto L^{9/76}$. These slopes are illustrated with arbitrary separate reference-model normalizations:

The proton–proton/Kramers model applies to the lower-mass regime; the CNO/electron-scattering model applies to the higher-mass regime. Higher-mass cores are hotter, making the more temperature-sensitive CNO cycle efficient, and highly ionized material favors electron-scattering opacity. These are the stated idealized radiative, gas-pressure-supported sequences: a fully convective low-mass star or a radiation-pressure-dominated massive star needs a different structural model. The intersection of independently normalized lines in the figure does not specify a physical crossover mass.

Take the gas pressure $P$ to vanish at the exterior boundary and the atmosphere to be thin enough that surface gravity of a star is constant. With constant electron-scattering opacity, negligible radiative acceleration, and hydrostatic equilibrium in optical depth,If the radiative force is retained, $dP/dr=-\rho (g-\kappa F/c)$ for gas pressure, and the same result holds with the constant positive effective gravity $g_{\mathrm{eff}}=g-\kappa F/c$. A nonzero exterior gas pressure instead gives $d\log P/d\tau =1/(\tau +\kappa P(0)/g_{\mathrm{eff}})$. Total gas-plus-radiation pressure has the nonzero radiative surface term from part (ii), so the advertised exact $1/\tau$ formula presumes gas pressure with zero exterior gas pressure.

Retain the same stellar homology and gas-pressure assumptions, but use the specified CNO cycle law $ϵ\propto \rho T^{16}$ and constant electron-scattering opacity. The nuclear luminosity now scales asThe radiative diffusion in a star scaling becomesEquating the two yields $R^{19}\propto M^{15}$. Therefore radiative homology with CNO burning and electron scattering givesSince $R\propto L^{5/19}$, the Stefan–Boltzmann law givesThis is a model comparison using the exponent 16 prescribed here. Actual massive stars may have appreciable radiation pressure and convective cores, which violate the assumptions of this purely radiative gas-pressure family.

Treat each grain as an isolated compact sphere and take its wavelength-independent radiation-pressure efficiency to be $Q_{\mathrm{pr}}=1$. At distance $R$, its mass and radiation-force cross-section areThe inward gravitational force and outward radiation pressure force areThus the opacity per unit dust mass and the dust-grain Eddington limit areTo equal the Eddington luminosity based on electron-scattering opacity, ${\kappa }_d$ must equal ${\kappa }_{\mathrm{es}}$. ThereforeFor fully ionized hydrogen, ${\kappa }_{\mathrm{es}}={\sigma }_T/m_p\simeq 0.40\mathrm{c}{\mathrm{m}}^2{\mathrm{g}}^{-1}$, and the stated ${\rho }_d=1\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$ givesUsing $0.34\mathrm{c}{\mathrm{m}}^2{\mathrm{g}}^{-1}$ for a roughly solar-composition gas instead gives $2.2\mathrm{cm}$. This composition-dependent normalization should be stated because the problem does not fix it. For a different constant radiation-pressure efficiency, the critical radius is multiplied by $Q_{\mathrm{pr}}$.

At $L=L_{\mathrm{Edd},\mathrm{es}}$, the grain's radiation-pressure coefficient isConsequently, smaller grains experience a net outward force, larger grains remain gravitationally attracted, and critical-size grains have zero net radial force. Grains released at rest are therefore size-sorted into expelled and infalling populations. With an initial inward velocity, an outward force does not guarantee reversal before the grain reaches the hole; trajectories also depend on initial conditions. Gas drag and dust destruction have been omitted from this isolated-grain calculation. Ordinary interstellar grains are much smaller than a centimetre, so their direct radiative acceleration is very strong in this idealized model.

The usual Pressure-opacity zones of a Shakura--Sunyaev thin disk, in order of increasing cylindrical radius, are an inner region dominated by radiation pressure and electron-scattering opacity, a middle region dominated by gas pressure and electron-scattering opacity, and an outer region dominated by gas pressure and the Kramers' opacity law, often approximated by free-free opacity. The boundaries depend on mass, mass accretion rate, and the alpha disk viscosity parameter. At sufficiently large radius, accretion-disk self-gravity or changes in ionization can invalidate this three-zone model.

For the inner region writeHere $F$ is the standard thin-disk dissipation flux through one face, and $R_{\ast }$ is the innermost stable circular orbit treated as a zero-torque inner boundary. The given factor $1/2$ in the dissipation rate is appropriate to one face. For Keplerian rotation, it is $(9/8)\nu \Sigma {\Omega }_K^2$.

Let $H$ denote the disk scale height or vertical half-thickness, not the full thickness. In vertical hydrostatic equilibrium, the leading vertical acceleration is ${\Omega }_K^{2}z$, so a one-zone estimate givesIn the inner zone, $P_c\simeq P_{\mathrm{rad},\mathrm{c}}=a_{\mathrm{rad}}{T}_c^4/3$. For an optically thick medium with constant electron-scattering opacity, radiative diffusion gives the emergent one-face fluxEquivalently, combine $F=-c({\kappa }_{\mathrm{es}}\rho {)}^{-1}d{P}_{\mathrm{rad}}/dz$ with vertical pressure balance. Equating the radiative flux to the dissipated flux yields the radiation-supported inner-disk height:Order-unity vertical-profile factors depend on the precise height convention. The central mass cancels at fixed physical mass accretion rate, because both vertical gravity and local dissipation are proportional to $G{M}_{\mathrm{BH}}/R^3$.

If ${\dot{m}}_{\mathrm{Edd}}=L_{\mathrm{Edd},\mathrm{es}}/({\eta }_0{c}^2)$ and $r_g=G{M}_{\mathrm{BH}}/c^2$, this can also be writtenThe radiation-supported height rises from the formal zero at $R_{\ast }$ toward the approximately constant value $H_0=3{\kappa }_{\mathrm{es}}\dot{m}/(8\pi c)$. That zero is not a reliable description of the physical plunging region: the radiation-dominated approximation and the Newtonian zero-torque formula fail very close to the inner edge.

Beyond the inner zone, the height grows slowly with radius rather than staying on that plateau forever. Far from the inner boundary, the standard gas-pressure scalings give $H\propto R^{21/20}$ in the electron-scattering zone and $H\propto R^{9/8}$ in the Kramers zone. These follow from $H\propto T_c^{1/2}{R}^{3/2}$ with $T_c\propto R^{-9/10}$ and $R^{-3/4}$, respectively. A schematic profile is

vertical half-thickness H
  ^
  |                                               / outer gas/Kramers
  |                                          ____/
  |                                     ____/ middle gas/scattering
  |              ______________________/
  |          ___/ inner radiation/scattering plateau
  |       __/
  |     _/
  +----|----------------------------------------------> radius R
      R_*              zone boundaries are schematic

The height formula also tests the thin disk assumption. For a nonrotating hole with $R_{\ast }=6r_g$, the largest $H/R$ in this inner approximation occurs at $R=(9/4)R_{\ast }$. At $\dot{m}={\dot{m}}_{\mathrm{Edd}}$ and ${\eta }_0=0.1$, it is about $0.37$, so a disc this close to the conventional Eddington accretion rate is only marginally thin. Slim accretion disk effects and advected heat can then matter.

Analyze a local temperature perturbation at fixed radius and fixed surface density of a disk $\Sigma$. The thermal timescale of an accretion disk is shorter than its viscous timescale, so mass redistribution is negligible during the perturbation. Let $Q^{+}$ be viscous heating and $Q^{-}$ radiative cooling, measured consistently through one face or both faces, and define net cooling by $\dot{Q}=Q^{-}-Q^{+}$. A temperature increase is unstable if it decreases net cooling.

With gas pressure dominant, vertical hydrostatic equilibrium and the ideal gas law giveFor the alpha disk prescription with temperature-independent $\alpha$, $\nu =\alpha c_{s}H\propto T_c$, and the Keplerian rotation heating law isOne must not hold the steady-state $\dot{m}$ fixed while perturbing the temperature: the stress and heating respond through $\nu$, even though $\Sigma$ is temporarily fixed.

At an equilibrium temperature $T_0$ where $Q^{+}=Q^{-}=Q_0$, if $Q^{-}\propto T_c^q$ and $Q^{+}\propto T_c$, thenThe cooling exponent therefore decides the sign.

For hot, ionized, optically thin gas, the usual intended cooling model is thermal bremsstrahlung, whose volume emissivity scales as ${\rho }^2{T}_c^{1/2}$. Integrated vertically,Thus $q=0$, andA hotter annulus expands vertically, lowering its density enough to offset the direct increase in bremsstrahlung emission; viscous heating still rises.

For an optically thick medium with the Kramers' opacity law, $\kappa \propto {\rho }_c{T}_c^{-7/2}\propto T_c^{-4}$ at fixed $\Sigma$. Thermalized radiative diffusion givesIts falling opacity lets a hotter annulus radiate much more efficiently.

For an optically thick medium with dominant electron-scattering opacity, $\kappa$ is nearly independent of temperature. Provided enough true absorption is present to thermalize the radiation,This conclusion differs from the thermal instability of a radiation-pressure-dominated alpha disk, where the heating has a different temperature dependence. The question explicitly excludes that pressure regime.

The first conclusion needs its cooling assumption stated. $T>10^4\mathrm{K}$ alone does not guarantee that thermal bremsstrahlung dominates; atomic-line cooling, ionization changes, external illumination, or energy advection can change the answer. For a general optically thin emissivity ${\rho }^2\Lambda (T_c)$, the corresponding exponent isAccordingly, under the same fixed-$\Sigma$ alpha prescription,The unstable/stable/stable classification is the standard free-free/diffusion answer, rather than a universal result specified by temperature and optical depth alone.

At $x=1/2$, defineDifferentiating the temperature profile givesSince the luminosity is already constant there, the radiative diffusion in a star equation yields

For the Kramers' opacity law $\kappa ={\kappa }_0\rho T^{-7/2}$, substitution giveswhereThus $\alpha =11/2$ and $\beta =-1/2$ in the notation of the question.

For constant electron-scattering opacity ${\kappa }_{\mathrm{es}}$,whereHence $\delta =3$ and $\gamma =0$: the radius cancels, and the mass dependence is shallower than under Kramers opacity.

The standard Shakura--Sunyaev thin disk has three qualitative radial zones. Its hot inner part is dominated by radiation pressure and electron-scattering opacity; farther out, gas pressure overtakes radiation pressure while electron scattering can remain the main opacity; in the cool outer zone, gas pressure remains dominant and free-free opacity becomes important. The transition radii vary with $M$, $\dot{M}$, and $\alpha$.

For the inner zone, hold the surface density $\Sigma$ fixed during a local thermal perturbation. Vertical hydrostatic equilibrium givesso $H\propto T_c^4/\Sigma$. The alpha disk prescription then yieldswhereas optically thick radiative diffusion with nearly constant electron-scattering opacity givesAt equilibrium $Q^{+}=Q^{-}=Q_0$. For net cooling $\dot{Q}=Q^{-}-Q^{+}$,Thus the total-pressure alpha prescription predicts thermal instability of a radiation-pressure-dominated alpha disk: a temperature increase makes heating outrun cooling. The absence of ubiquitous, large-amplitude thermal limit cycles in luminous AGN light curves indicates that this local model omits stabilizing effects, plausibly magnetic pressure and stress, vertical advection, winds, or a stress law that does not simply track total pressure.

The inner Shakura--Sunyaev thin disk is supported mainly by radiation pressure and has nearly constant electron-scattering opacity. Equating radiative diffusion with local dissipation and imposing vertical hydrostatic equilibrium gives the characteristic half-thicknessThe radiation-supported height approaches a plateau away from the inner boundary. The formal zero at that boundary is outside the reliable radiation-dominated approximation. At an Eddington accretion rate defined using efficiency $0.1$, this height can make the inner flow only marginally thin.

Using the idealized CNO cycle law $ϵ\propto \rho T^{16}$ and constant electron-scattering opacity, homology relations for radiative stars giveThese retain ideal-gas pressure support; significant radiation pressure or convection changes the model.