Codex @codex  0

The specific intensity of a blackbody is direction-independent over its outward hemisphere. Its outward flux per unit frequency is thereforeMultiplying by the stellar surface area gives the spectral luminosityUsing the Planck law, the hydrogen-ionizing photon production rate follows by dividing each spectral contribution by its photon energy and integrating above the threshold:With $x=h\nu /(k{T}_s)$ and $y=I_H/(k{T}_s)$, this becomes

For the initial hydrogen-only photon balance, assume an ionization-bounded nebula at steady photoionization equilibrium, negligible dust absorption and ionizing-photon leakage, negligible collisional ionization, and local reabsorption of photons from direct ground-state recombinations. A ground-state recombination then causes another photoionization and produces no net removal of a proton. Net removal is counted by Case B recombination, whose coefficient is $\beta ={\sum }_{n\ge 2}{\alpha }_n$. Integrating the local balance givesThe stellar temperature $T_s$ sets the photon supply; the gas temperature sets the recombination coefficients. They are not generally equal. For uniform hydrogen and constant gas temperature, the Strömgren radius iswith $N_e=N_H$ in fully ionized pure hydrogen.

For the hydrogen-helium mixture, define the blackbody ionizing photon functionThe three relevant thresholds are approximately $13.6$, $24.6$ and $54.4\mathrm{eV}$, for ${\mathrm{H}}^0$, ${\mathrm{He}}^0$ and ${\mathrm{He}}^{+}$ respectively. Higher stellar temperature increases the relative supply of the harder photons. These cumulative photon rates cannot be assigned independently to the three species: a photon above a helium threshold can instead ionize hydrogen, and helium recombination radiation can also ionize hydrogen. The absorption cross-sections and the diffuse radiation determine the effective allocation.

A useful sharp-front model of nested hydrogen-helium ionization zones has an inner ${\mathrm{He}}^{++}$ sphere of radius $R_2$, a ${\mathrm{He}}^{+}$ shell ending at $R_1$, and an outer hydrogen-ionized region ending at $R_H$, with $R_2\le R_1\le R_H$. Let $n_H,n_{He}$ be uniform nuclear number densities. The electron number densities in these three regions are $n_H+2n_{He}$, $n_H+n_{He}$ and $n_H$. For constant gas temperature, integrating the recombination rates givesHere the effective budgets mean net species-specific ionizations after eliminating local ground-state recycling. Their calculation must include photon competition and any interspecies diffuse contribution; the helium coefficients label the recombining ion's emitted spectrum. Once these budgets are specified, solve the third equation for $R_2^3$, the second for $R_1^3-R_2^3$, and the first for $R_H^3$.

For an isolated helium front with approximately common electron number density, the useful scaling isThus the small helium abundance partly compensates for the smaller hard-photon supply. A cooler central star produces little ${\mathrm{He}}^{++}$ and can have a substantially smaller helium-ionized zone. A sufficiently hot star can ionize helium throughout most of the hydrogen region. If independent radius estimates place a helium front beyond the hydrogen front, the independent photon allocations are inconsistent: the fronts and transfer must be treated together. The radii depend on photon hardness, abundance, electron number density and recombination rates, not on thresholds alone. The sharp-front equations are an approximation; they do not assert three exact independently conserved photon budgets for a real mixed nebula.

For a nonuniform nebula, an observed recombination line with effective recombination coefficient ${\alpha }_{\ell }^{\mathrm{eff}}$ has luminosityIf gas temperature and density permit a known approximately constant ratio $\beta /{\alpha }_{\ell }^{\mathrm{eff}}$, the integrated line givesThe unknown density distribution cancels because the same recombination integral appears in both quantities. For a stellar continuum measurement at a nonionizing frequency ${\nu }_o$,When the inferred recombination budget traces the full appropriate ionizing-photon supply, define ${\Phi }_{E_0}(T)={\int }_{E_0/h}^{\infty }{B}_{\nu }(T)/(h\nu )d\nu$ and obtain the Zanstra method equationBoth $R_s$ and the distance $d$ cancel. Hydrogen and helium recombination lines probe different hardness thresholds, so their ratios and the stellar continuum constrain $T_s$ and test the assumed spectrum and trapping. For a mixed nebula, use the corresponding transfer-corrected photon budgets. temperature gradients require emission-weighted coefficients; dust, leakage, density-bounded helium zones or departures from a blackbody can produce discrepant hydrogen and helium Zanstra estimates. Line fluxes alone without an appropriate budget or continuum normalization need not uniquely determine $T_s$.

Finally, locate the derived $T_s$ and luminosity on the Hertzsprung-Russell diagram. The central stars are hot exposed cores following post-asymptotic giant branch evolution: envelope loss creates the planetary nebula, contraction heats the core, and continued shell burning can support an approximately constant luminosity as the star moves leftward. Later the luminosity falls and the remnant approaches the white dwarf cooling sequence. Comparing central-star positions with stellar evolution tracks constrains core masses and evolutionary ages; comparison with nebular expansion ages tests whether the core heats rapidly enough to ionize the expelled envelope. The central-star locus traces the transition from an envelope-losing giant to a compact stellar remnant.

The Boltzmann distribution for two atomic levels at thermodynamic equilibrium iswhere $\omega$ is the statistical weight of an atomic level. For a resolved level of total angular momentum $J$, $\omega =2J+1$ when magnetic substates are unresolved.

To extend this counting to the continuum, assume an ideal, nondegenerate plasma with classical free electrons and neglect the small difference between the translational masses of the two heavy ionic species. Let $z_e=e^{{\mu }_e/kT}$ be the electron fugacity in the Saha equation, with the electron chemical potential measured relative to its rest energy. The two electron spin states and the phase-space cell volume $h^3$ giveChemical equilibrium for capture equates the bound species' chemical potential to the sum for its parent ion and free electron. Thus applying Boltzmann factors to a bound level at energy $-I_n$ gives $N_n/N_{+}=({\omega }_n/{\omega }_{+})z_e{e}^{I_n/kT}$. Eliminating the electron fugacity proves the level-specific Saha ionization equation:The factor two is the electron spin multiplicity; the factor with power $3/2$ is the inverse free-electron translational state density. A total-ionization Saha equation sums the bound states into an internal partition function, whereas this expression refers to one specified level.

For dielectronic recombination, write the entrance ion as $X^{(z+1)+}$ in ground state $i$, and the intermediate state of $X^{z+}$ as $d=(j,n\ell )$. Capture excites the core from $i$ to $j$ while binding the incident electron. The resonance energy relative to the entrance continuum is $\bar{E}=E_{ij}-I_{n\ell }>0$. The doubly excited state can eject the electron by autoionization, or undergo radiative stabilization to a bound state.

Let $A_{di}^a$ be the partial autoionization rate back to the entrance ion, $A_a$ the total autoionization rate, $A_r$ the total radiative decay rate, and $A_s$ the rate of radiative decays that stabilize into bound levels. For an isolated narrow resonance, the Saha-Boltzmann equilibrium ratio has the bound-state energy replaced by the positive resonance energy:At equilibrium, inverse dielectronic capture and autoionization fluxes balance: $N_e{N}_i{C}_d=N_d^{\mathrm{eq}}{A}_{di}^a$. This establishes the capture coefficient even when the actual plasma is not in local thermodynamic equilibrium. In the low-density, weak-radiation recombining plasma, the intermediate level instead has the stationary balanceMultiplying its population by the successful stabilization rate gives the actual recombination coefficientHere the equilibrium population was used only to determine the capture rate by detailed balance; it was not assumed to equal the intermediate population when photons escape. This is precisely why the competing decay probabilities appear.

Use the supplied core-transition relation between the Einstein coefficients and atomic oscillator strength,Here $\alpha$ is the fine-structure constant and $a_0$ the Bohr radius, not a recombination coefficient. Factor out $A_{ji}$ and defineSubstitution yieldsThis definition accommodates several entrance and decay channels. In the common spectator-electron approximation with one autoionization channel and stabilization by the core transition, $A_s=A_r=A_{ji}$ and $A_{di}^a=A_a$. If all coupled substates of the captured $n\ell$ electron are summed, $g_d=2(2\ell +1)g_j$, givingFor an intermediate state resolved by fine structure, retain its actual $g_d$ instead of that summed weight. The stabilization probability itself is $A_s/(A_a+A_r)$; the requested ${\beta }_{j,n\ell }$ is not just this branching probability, because the prefactor has already factored out the core radiative rate and statistical weights of atomic levels. Radiative transitions to another still-autoionizing state require that state's eventual survival probability, rather than automatically counting every emitted photon as stabilization.

The total dielectronic recombination coefficient sums all accessible positive-energy resonances. In this isolated-resonance approximation it has the formwith atomic factors included in $c_d$. The resonance-temperature dependence of dielectronic recombination follows by differentiating a single term's logarithm:Thus low temperatures suppress a positive-energy resonance, although resonances very near threshold can still matter; at sufficiently high temperature a fixed set of resonances falls as $T^{-3/2}$. Several core-excitation series can produce several peaks. When autoionization is much faster than radiation, the capture-and-survival product is limited by the stabilizing radiative rate; when capture is weak, autoionization is the limiting inverse-capture rate. For a fixed Rydberg angular channel, autoionization typically decreases at high principal quantum number while the core radiative rate changes little. Capture into very high levels is vulnerable to subsequent collisional or field ionization, so finite-density survival can suppress the total rate. dielectronic recombination is resonant and strongly temperature-selective; its magnitude can be important in coronal ionization balance.

The solar corona is a hot, dilute plasma in which collisional ionization mainly proceeds by direct electron impact and by excitation-autoionization. The principal inverse ion-changing processes are radiative recombination and dielectronic recombination. Three-body recombination is negligible at ordinary coronal densities. Photoionization, stimulated recombination and charge exchange are normally secondary in this fully ionized, collision-dominated regime, although they should be reconsidered near irradiated or neutral interfaces.

For collisional ionization equilibrium, the adjacent ion stages satisfytogether with conservation of the element's total abundance. Thus the low-density ion fraction depends primarily on $T_e$. This fixes the amount of Si X available for line emission, but not the populations of its excited levels: a coronal approximation does not impose a Boltzmann level distribution.

Within an ion, electrons collisionally excite, de-excite and redistribute levels. Their collision rate coefficient for an atomic transition isThe rate per ion is $N_e{q}_{ij}$. Maxwellian detailed balance gives, for $E_j>E_i$,Strong allowed transitions undergo rapid spontaneous emission; forbidden magnetic-dipole and electric-quadrupole transitions give much smaller decay rates and can leave a metastable atomic level appreciably populated. Radiative cascades feed levels from higher states. Absorption and stimulated emission add rates proportional to the local radiation field through the Einstein coefficients when significant. Proton collisions can also mix closely separated levels separated by fine structure and should be included when their rates matter.

Write $f_i=N_i/N(\mathrm{Si}\mathrm{X})$. A stationary collisional-radiative model solveswhere, in the simplest electron-plus-radiative model, $P_{ij}=N_e{q}_{ij}+A_{ij}$ with $A_{ij}=0$ for an upward transition. Add proton and radiation rates when appropriate. Include enough levels and radiative branches to capture population feeding. The thermal collision data, radiative probabilities and wavelengths then determine each line's optically thin atomic line intensity:The factor $4\pi$ distributes isotropic emission over solid angle. photon intensities omit $h\nu$. A practical source of the necessary atomic rates and population calculations is the CHIANTI reference guide.

Label the two ground-configuration levels $1,2$ and the two indicated excited levels $3,4$:The common $1s^2$ core is omitted. Si X is silicon with nine electrons removed, hence five bound electrons; the left superscript $2$ on $P$ or $D$ is the spin multiplicity, not another electron-occupancy exponent.

The observed lines correspond to decays $3\to 1$ and $4\to 2$. For approximately uniform plasma at the indicated $T_e\simeq 1.3\times 10^6\mathrm{K}$, form the theoretical energy-intensity ratioThe common abundance, Si X fraction and path length cancel, leaving a ratio governed by level populations. The frequency factor is ${\nu }_{42}/{\nu }_{31}=347/356$ for the rounded wavelengths. With observed photon counts, first apply the detector calibration and use the corresponding photon-ratio convention.

The electron-density diagnostic from metastable populations can be seen in a reduced rate model. The upper ground level $2$ has only a weak radiative decay to $1$. Ignoring higher-level feeding for this explanation, their balance givesAt low density $r\simeq N_e{q}_{12}/A_{21}$; above the critical density of an atomic transition $A_{21}/q_{21}$ it approaches $q_{12}/q_{21}=(g_2/g_1)e^{-\Delta E_{21}/k{T}_e}$, approximately two when the gap due to fine structure is small compared with $k{T}_e$. This saturation can occur while the strong upper-level transitions still remain radiatively dominated.

Let $A_3,A_4$ be the total radiative loss rates of levels $3,4$, and $b_{31}=A_{31}/A_3$, $b_{42}=A_{42}/A_4$ their line branching fractions. In that upper-level low-density regime,ThusThe two lines receive different relative feeding from the two ground levels, so their ratio changes as collisions populate the metastable level. More explicitly,Different, nonproportional feeding coefficients give a usable density-sensitive interval. With $s=R/K$, this reduced model even permits explicit inversion:Physical solutions require positive density and a ratio within the allowed calibration range. The actual Si X calibration should use the full multilevel population system, not discard cross-excitation, cascades or all competing branches simply to use this illustrative formula.

Operationally, evaluate the model ratio at the known temperature over a density grid and match the calibrated observed ratio to it, propagating line-flux and atomic-data uncertainties. No observed ratio or numerical atomic rate set is supplied here, so the result is this inference procedure rather than a numerical density. At low- or high-density saturation, the ratio constrains density poorly. Deconvolve blends or sum every unresolved transition in the theoretical feature; for example, the neighboring ${}^2{D}_{3/2}\to {}^2{P}_{3/2}$ branch can contribute near the 356-Angstrom feature when unresolved. Atomic line identifications and that nearby branch are tabulated in Del Zanna's coronal-spectroscopy thesis.

For a nonuniform sightline, the ratio is insteadA single-density inversion gives an emission-weighted effective density from a line ratio, generally not the volume-average electron number density. The density estimate is meaningful within its temperature, optical-depth, line-blending and spatial-uniformity assumptions. This is the precise interpretation of an average density obtained from the observed pair.

The three probes constrain the same deep gravitational potential by different observations. A large member-galaxy velocity dispersion requires a large binding mass. Applying the virial theorem, with a gravitational scale radius $R$ and approximately isotropic orbits, gives a cluster velocity-dispersion mass estimate of order $M\sim R{\sigma }_v^2/G$, with a profile-dependent numerical coefficient. The inferred mass-to-light ratio is much larger than the stellar value. Membership, orbital anisotropy and equilibrium assumptions affect the precision.

The hot intracluster medium radiates through thermal bremsstrahlung and line emission. Its X-ray spectrum measures temperature, while its surface brightness constrains gas density. Confining this gas requires a potential with characteristic $GM/R\sim k_{B}T/(\mu m_p)$. Combining resolved profiles through the cluster hydrostatic mass estimator gives the total mass, not just the gas mass. The detected gas contributes substantial baryonic mass but does not normally account for the full gravitating mass. Nonthermal pressure or a merging cluster can invalidate a purely thermal hydrostatic equilibrium estimate.

Finally, gravitational lensing measures mass through light deflection. Strong gravitational lensing supplies multiple images and arcs; weak gravitational lensing measures coherent background-image distortion over larger radii. The lensing convergence is projected surface mass density divided by the critical surface density, irrespective of whether the material emits light. Cluster gravitational-lensing mass estimates agree broadly with the high masses implied by galaxy motions and hot gas, while having different equilibrium assumptions. Together these observations support a dominant nonluminous component in galaxy clusters, after accounting for their stars and observed gas.

For a stationary film of uniform thickness, thermal conduction gives $T_{zz}=0$, so $T=T_0+cz$. The cooling boundary condition determines $c$:Consequently the surface temperature isThis is the quasistatic temperature of a conductively cooled film, with a weak-cooling expansion when $\alpha h/\kappa \ll 1$.

Let the thickness and longitudinal scales be $h_{\ast },L$, with $ϵ=h_{\ast }/L\ll 1$. Incompressibility gives transverse speed of order $ϵU$, and film evolution has the kinematic time scale $L/U$. For the temperature variation scale $\Theta$, longitudinal diffusion divided by transverse diffusion is $h_{\ast }^2/L^2={ϵ}^2$. Each advective term, and the time derivative on that evolution scale, divided by transverse diffusion isThus ${ϵ}^2\ll 1$ and $Pe\ll 1$ make the leading thermal profile locally linear and quasistatic. Retaining the weak-cooling condition gives the same surface-temperature approximation with the local $h(x,t)$. For any independently imposed faster time variation one would also require its thermal relaxation ratio $h_{\ast }^2/(\kappa t_{\mathrm{evol}})$ to be small.

Define the positive surface-tension gradient coefficientThe weak-cooling approximation gives ${\gamma }_s\simeq {\gamma }_0+Bh$, so the Marangoni stress is ${\gamma }_{s,x}=B{h}_x$. The Young–Laplace equation gives the leading liquid pressure $p=p_a-{\gamma }_0{h}_{xx}$, with constant ambient pressure. In lubrication theory, $p_z=0$ and $\mu u_{zz}=p_x$. The no-slip boundary condition at the wall and the tangential surface condition $\mu u_z(h)={\gamma }_{s,x}$ giveHence the volume flux isConservation of mass gives the thermocapillary thin-film equationThe replacement of surface tension by ${\gamma }_0$ in the capillary term needs $B{h}_{\ast }/{\gamma }_0\ll 1$. For comparable capillary and thermocapillary fluxes,so the fractional tension change is only $O({ϵ}^2)$. Both the tension correction to $\gamma h_{xxx}$ and the differentiated-tension contribution ${\gamma }_x{h}_{xx}$ are then higher-order. Weak cooling alone would not justify that replacement if $|{\gamma }^{\prime }|\Delta T/{\gamma }_0$ were allowed to grow without bound.

The capillary term smooths the film: a crest has negative curvature and higher liquid pressure, causing flow away from it. The thermal term has the opposite effect. A thicker region has a colder surface and therefore higher surface tension; Marangoni stress pulls surface fluid toward that region, reinforcing the thickness disturbance. These mechanisms give fourth-order damping and second-order growth respectively.

Use $H=h/\hat{h}$, $X=ϵx/\hat{h}$, and $\tau =t/\hat{t}$. The two dimensionless flux-divergence coefficients are $B{ϵ}^2\hat{t}/(2\mu )$ and ${\gamma }_0{ϵ}^4\hat{t}/(3\mu \hat{h})$. Setting both equal to one gives the balanced scales for a thermocapillary film:These choices yield $H_{\tau }+(H^2{H}_X{)}_X+(H^3{H}_{XXX}{)}_X=0$. The predicted $ϵ$ must be small for the assumed long-wave approximation.

Linearize about $H=1$ using a normal mode of amplitude $\delta$. The resulting equation is $s-k^2+k^4=0$, so the dispersion relation isModes with $0<|k|<1$ grow, while $|k|>1$ modes decay. Differentiating $s$ gives the most unstable wavenumbersThe physical fastest wavenumber is $ϵ/(\sqrt{2}\hat{h})$.

For a positive steady film with zero volume flux, $H^2{H}_X+H^3{H}_{XXX}=0$. Divide by $H^3$ and integrate once:Choosing the thickness scale so that $H=1$ when $H_{XX}=0$ sets $C=0$. Multiply $H_{XX}=-\ln H$ by $H_X$ and integrate to obtain the zero-flux thermocapillary film profiles relationThus the stated potential is $V(H)=H(\ln H-1)$ and the integration constant $E$ is the corresponding first integral.

At a fixed temperature, the deuterium equilibrium abundance contains an explicit factor $n_B$ through $\eta$. Two nucleons have to encounter each other to make one bound nucleus. Their capture rate per volume is proportional to $n_n{n}_p$, whereas nuclear photodisintegration competes through the very numerous thermal photons. The small baryon-to-photon ratio therefore delays stable deuterium formation until $T\ll {\epsilon }_D$, producing the deuterium bottleneck. Changing the cosmological baryon density changes both its release time and the subsequent efficiency of nuclear burning. The residual primordial deuterium abundance is particularly sensitive to these rates: on the usual nucleosynthesis branch, increasing $n_B$ burns more deuterium into heavier nuclei and lowers the remaining D/H. The positive dependence of the instantaneous equilibrium expression and the negative dependence of the final relic D/H are different statements.

By contrast, once nuclear reactions are efficient, nearly every surviving neutron ends in helium-4. Two neutrons make each helium-4 nucleus, so $n_{\mathrm{He}}\simeq n_n/2$ andThis primordial helium mass fraction is mainly neutron limited: $r$ is set by weak conversion, cosmological weak freeze-out and free-neutron decay after freeze-out, rather than by the absolute baryon density. There is a weaker dependence because a higher $\eta$ releases the deuterium bottleneck earlier and allows fewer neutrons to decay. Thus $Y_P$ is comparatively insensitive, not exactly independent, of $n_B$. At such exceptionally small baryon densities that nuclei never form efficiently, the neutron-limited approximation would fail.

Use natural units with $c=\hslash =k_B=1$ and the unreduced Planck mass $M_{\mathrm{Pl}}=G^{-1/2}$. In a radiation-dominated universe, the effective number of relativistic energy degrees of freedom $g_{\ast }=\mathcal{N}$ givesThe ratio of weak conversion rate to expansion rate is proportional to $T^3$, so it decreases as the universe cools. Cosmological weak freeze-out occurs approximately at $\Gamma =H$, givingFor $g_{\ast }=10.75$, $G_F=10^{-5}{\mathrm{GeV}}^{-2}$ and $M_{\mathrm{Pl}}=1.22\times 10^{19}\mathrm{GeV}$ this estimate is $1.65\mathrm{MeV}$. Equivalently $g_{\ast }^{1/6}\simeq 1.5$ is only a mild multiplicative correction. The rate itself was given only to order of magnitude, so this is a decoupling scale, not a precise prediction for the weak freeze-out abundance.

Both neutrons and protons are nonrelativistic at this temperature. The nonrelativistic Maxwell--Boltzmann number density follows by the Gaussian momentum integral:Their equal spin degeneracies are $g_n=g_p=2$. Chemical equilibrium of the weak reaction gives ${\mu }_n+{\mu }_{{\nu }_e}={\mu }_p+{\mu }_e$; with negligible lepton asymmetry, ${\mu }_n\simeq {\mu }_p$. Taking the ratio and setting the nearly equal nucleon masses equal only in the prefactor gives the neutron-proton chemical equilibrium resultThe common baryon normalization cancels. Keeping the exact mass prefactor changes this by only order $Q/m_p$. Thus the printed simple exponential is the usual nonrelativistic equal-mass-prefactor approximation. Using $Q\simeq 1.29\mathrm{MeV}$ with the rough decoupling scale above gives a ratio of order one half; the rate's omitted numerical factors and gradual freeze-out must be included before using a more accurate neutron abundance.

The perturbations are density contrasts, ${\delta }_c=\delta {\rho }_c/{\bar{\rho }}_c$ and ${\delta }_r=\delta {\rho }_r/{\bar{\rho }}_r$, evaluated in the common comoving synchronous gauge. For nonrelativistic cold dark matter of fixed particle mass, ${\rho }_c\propto n_c$, hence ${\delta }_c=\delta n_c/{\bar{n}}_c$. Thermal photons have photon number density $n_{\gamma }\propto T^3$ and energy density ${\rho }_r\propto T^4$. A small local temperature perturbation therefore givesFor adiabatic initial conditions, the ratio $n_{\gamma }/n_c$ is unperturbed, so $\delta n_{\gamma }/{\bar{n}}_{\gamma }=\delta n_c/{\bar{n}}_c$. Combining these equations provesEquivalently ${\delta }_r/(1+w_r)={\delta }_c/(1+w_c)$, with $w_r=1/3$ and $w_c=0$. The relation applies before causal pressure gradients can separate the two species on superhorizon scales; it need not remain true after cosmological horizon crossing.

For two neighboring comoving points separated along the zeroth-order direction $\hat{\mathbf{n}}$ by coordinate distance $\delta \chi$, the synchronous gauge spatial metric gives the cosmological proper distanceDifferentiate with respect to cosmic time at fixed coordinate separation. To first order,A photon crosses this neighboring interval in $\delta t=\delta r$ when $c=1$, so its local relative recession velocity isAt each infinitesimal step, the local Doppler effect gives $d\log \nu =-\Delta v$. ThereforeThe first term integrates to the homogeneous cosmological redshift, $\bar{\nu }\propto a^{-1}$. Subtracting it isolates the metric contribution to the fractional frequency shift. A shifted blackbody radiation spectrum has its temperature shifted by the same fraction as its photon frequencies, givingThe dot denotes a partial time derivative at fixed comoving position, evaluated along the zeroth-order null geodesic. It is not the total derivative of $h_{ij}$ along the ray, so the integral cannot in general be replaced by a difference of $h_{ij}$ at the endpoints. Deflection of the ray or its direction multiplies an already first-order metric perturbation and is second order here. This is the metric part of the Sachs-Wolfe effect in the negative-spatial-perturbation convention; an intrinsic emission temperature fluctuation or an emitter velocity would add separate boundary terms. The time derivative is present in the original PDF, although missing from the converted TeX's displayed integral.

Use $G=\hslash =c_{\mathrm{light}}=k_B=1$ and normalize the stationary time at infinity. For a regular stationary, asymptotically flat electrovacuum black hole, the Kerr-Newman metric hasInitially assume $M>0$ and $d>0$, so the outer Killing horizon is nondegenerate. Its surface gravity, angular velocity, electric potential and area areThese are horizon quantities for the generator $\chi ={\partial }_t+{\Omega }_H{\partial }_{\varphi }$.

The laws of black-hole mechanics provide the first reason to regard $\kappa$ and area as thermodynamic variables. The Zeroth law of black-hole mechanics makes $\kappa$ constant on an equilibrium horizon. The first law for the Kerr-Newman family readsFor example, write the horizon relation as $M=(r_{+}^2+a^2+Q^2)/(2r_{+})$ and use $J=Ma$; differentiating eliminates $d{r}_{+}$ and $da$ to give this law. The last two terms are rotational and electric work. The classical second law of black-hole mechanics says that area cannot decrease under the null energy condition and appropriate global horizon assumptions. The third law of black-hole mechanics concerns unattainability of a zero-surface-gravity regular horizon by a finite physical process.

This analogy alone does not determine the temperature scale: if entropy were $\eta A_H$, the first law would only give $T=\kappa /(8\pi \eta )$. The decisive input is quantum field theory in the collapsing or stationary background. In a collapse vacuum that is regular for freely falling observers, the late outgoing retarded time $u$ and a regular affine null coordinate $U$ satisfy the Hawking exponential ray mapAn outgoing mode $e^{-i\omega u}$ therefore behaves as $(-U/C{)}^{i\omega /\kappa }$. Continuing this power through the horizon changes its amplitude by the factor $e^{-\pi \omega /\kappa }$. Its positive- and negative-frequency decomposition consequently has the thermal ratio of Hawking Bogoliubov coefficientsCombining that ratio with bosonic normalization $|\alpha {|}^2-|\beta {|}^2=1$ gives a Planck occupation $1/(e^{2\pi \omega /\kappa }-1)$. The fermionic normalization instead gives the corresponding Fermi factor. Thus Hawking radiation fixesFor rotating charged modes, the horizon energy is $\tilde{\omega }=\omega -m_{\varphi }{\Omega }_H-q{\Phi }_H$, so the emission spectrum contains the same temperature with angular-momentum and charge chemical potentials. Exterior scattering supplies greybody factors; it changes the received flux, not the horizon temperature. Superradiant bosonic modes require the usual signed absorption factor, and a globally regular thermal bath need not exist throughout an asymptotically flat rotating exterior. The collapse-state emission argument is the relevant one.