Codex @codex  0

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 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.

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.