Codex @codex  0

Consider a narrow bundle of rays between two surface elements $d{A}_1,d{A}_2$ separated by $s$ in a transparent, stationary Euclidean vacuum. Let the ray make angles ${\theta }_1,{\theta }_2$ to the respective surface normals. The receiving element subtends solid angle $d{\Omega }_1=\cos {\theta }_{2}d{A}_2/s^2$ at the emitter. By the definition of specific intensity, the beam power in frequency interval $d\nu$ isAt the receiving end the source subtends $d{\Omega }_2=\cos {\theta }_{1}d{A}_1/s^2$, so the identical power isIn the absence of absorption, emission or frequency shifts, conservation of energy equates these expressions:This vacuum conservation of specific intensity also follows directly from the radiative transfer equation $d{I}_{\nu }/ds=0$. The apparent solid angle of an unresolved source shrinks as $s^{-2}$, so its integrated observed radiative flux still follows the inverse-square law. Specific intensity and unresolved flux are different quantities.

The vacuum and frequency assumptions are essential: absorption and emission change specific intensity, and gravitational or cosmological redshift conserves $I_{\nu }/{\nu }^3$ along a ray instead of $I_{\nu }$ itself. A fixed telescope aperture observing an unresolved object also cannot infer constant total received power from intensity conservation.

At the stellar surface boundary condition, match the interior to a stellar atmosphere, rather than setting the temperature and pressure to zero at an arbitrary radius. At a photosphere of radius $R$, the enclosed mass is $M$, the luminosity is $L$, and the outward radiative flux defines the effective temperature through $F=L/(4\pi R^2)={\sigma }_{\mathrm{SB}}{T}_{\mathrm{eff}}^4$. With no external illumination the incoming specific intensity satisfies $I_{\nu }(0,\mu )=0$ for $\mu <0$. A deeper boundary is matched to the nearly isotropic radiative diffusion field. A grey atmosphere commonly matches near $\tau =2/3$; under the Eddington surface boundary condition this gives $T=T_{\mathrm{eff}}$. For roughly constant gravity and opacity, the total pressure rises by $g\tau /\kappa$ relative to its outer boundary value; if radiation supplies appreciable support, the gas-pressure gradient uses the effective gravity $g-\kappa F/c$ instead.

Take $z$ increasing outward and $\mu$ the outward direction cosine. If $j_{\nu }$ is emission per unit path length and solid angle, the radiative transfer equation isDefine inward-increasing optical depth by $d{\tau }_{\nu }=-\rho {\kappa }_{\nu }dz$ and the source function by $S_{\nu }=j_{\nu }/(\rho {\kappa }_{\nu })$. Then $\mu {\partial }_{{\tau }_{\nu }}{I}_{\nu }=I_{\nu }-S_{\nu }$. This states the emission-coefficient convention; if $j_{\nu }$ is defined per unit mass instead, the emission term is $\rho j_{\nu }$.

For frequency-integrated intensity, introduce the radiation-field momentsThe first angular approximation writes $I\simeq J+3\mu H$. Its odd part does not contribute to $K$, soThis is the Eddington closure approximation, not an exact description of the escaping angular distribution. An approximate surface condition takes the outgoing hemisphere to have direction-independent intensity $I_0$ and the incoming hemisphere to be dark. Its moments give $F=\pi I_0$ and $J(0)=I_0/2=F/(2\pi )$. Integrating the moment equations with the closure then gives $J=3F(\tau +2/3)/(4\pi )$; in thermal equilibrium, $J=a_{\mathrm{rad}}c{T}^4/(4\pi )$ gives $T^4=(3/4)T_{\mathrm{eff}}^4(\tau +2/3)$. This derives the usual approximate photospheric matching condition rather than imposing a zero surface temperature. For conservative coherent isotropic scattering, $S=J$. Integrating the transfer equation gives $dH/d\tau =0$ and $dK/d\tau =H$, the moment conditions for radiative equilibrium.

The exponential ansatz in the PDF has a genuine consistency problem. Compute its moments before making any approximation. For a nonsingular real intensity on all $-1\le \mu \le 1$, take $|ϵ|<1$ and defineThe continuous zero limits are $A_0=1$, $A_1=0$, $A_2=1/3$. Direct integration yieldsFor a decaying nontrivial correction, $\delta >0$ and $\gamma \ne 0$. Exact constant flux then requires $A_1=0$. But $\mathrm{atanh}ϵ/ϵ>1$ for every nonzero real $|ϵ|<1$, so $A_1=0$ forces $ϵ=0$. Thus exact radiative equilibrium cannot imply a small nonzero $ϵ$ for this ansatz.

Substitution into the full transfer equation makes the obstruction sharper. The exponential coefficient must satisfy, for every $\mu$,Comparing the constant and linear coefficients gives $A_0=1$ and $\delta =A_0ϵ=ϵ$. Hence regular real coefficients require $\delta =ϵ=0$, eliminating any decaying correction. This is the obstruction to a single smooth exponential mode in conservative grey transfer. There is no exact nontrivial real solution of the stipulated form with a small nonzero decay parameter. Choosing, for example, $ϵ=\delta =0.1$ and $\gamma \ne 0$ supplies an explicit counterexample to the claimed exact flux constancy, because $A_1\ne 0$.

Several requested conclusions remain useful as formal approximations, provided they are identified as such. If $\delta >0$, the large-depth exponential vanishes, giving $F\to 4\pi \alpha /3$, hence $\alpha =3F/(4\pi )$ for the asymptotic physical flux. It also gives $c{P}_{\mathrm{rad}}\to (4\pi /3)J$. Since $A_0=1+{ϵ}^2/3+O({ϵ}^4)$, neglecting second-order angular errors yields the formal relation $\delta \simeq ϵ$. However $A_1=-ϵ/3+O({ϵ}^3)$, so the finite-depth flux still varies at first order. The formal relation is not an exact transfer or constant-flux solution.

There is also a flux-normalization mismatch in the PDF's last displayed relation. With the same physical $F$ used above, the Hopf function of a grey atmosphere is defined byThe PDF's coefficient $3F/4$ would require redefining its $F$ as physical flux divided by $\pi$; it cannot simultaneously use the earlier $\alpha =3F/(4\pi )$ for the same physical flux. With the consistent physical-flux normalization, formal moment matching givesIn the intended small-$ϵ$ approximation this isThe approximate outgoing surface ansatz then gives the limb darkening ratioFor finite $ϵ$, $\gamma =-0.14/A_0$ gives a parameter-dependent ratio, not a unique exact value. The two Hopf endpoint values alone do not determine the exact emergent intensity: the formal solution of the radiative transfer equation uses the entire source profile, $I(0,\mu )={\int }_0^{\infty }J(t)e^{-t/\mu }dt/\mu$ for $\mu >0$. Nor does the simple linear-plus-exponential ansatz satisfy the exact no-incoming-radiation condition for every negative $\mu$. Thus $2.75$ is a clearly labelled approximate disk-centre-to-limb result, not an exact consequence of inconsistent premises.

The outward radiative flux is $F=L/(4\pi r^2)$. A photon of energy $E$ carries momentum $E/c$, so the momentum transferred per unit volume and time to matter is $\kappa \rho F/c$, when $\kappa$ is the appropriate transport opacity. Momentum conservation makes the radiation stress decrease outward by this amount. In the thin atmosphere, where the radiation stress is treated by its scalar radiation pressure,This is the first moment of the radiative transfer equation. The grey atmosphere approximation uses a frequency-independent transport opacity; in an optically thick diffusion regime it becomes the Rosseland mean opacity. Near the surface the plane-parallel atmosphere approximation removes geometric stress-divergence corrections.

For the plane-parallel radiative transfer equationangular integration gives the zeroth radiation-field moment equationThe net radiative heating per unit volume is thereforeRadiative equilibrium requires it to vanish, equivalently that the frequency-integrated radiative flux be independent of depth:In local thermodynamic equilibrium with coherent isotropic scattering,where ${\omega }_{\nu }$ is the single-scattering albedo. Since ${\alpha }_{\nu }(1-{\omega }_{\nu })={\alpha }_{\nu ,\mathrm{abs}}$, the condition becomesConservative scattering redistributes directions but contributes no net material heating.

The first equation is vertical hydrostatic equilibrium: the pressure gradient balances the vertical gravity of the central mass, ${\Omega }^{2}z$, and the disk's own potential ${\Phi }_d$. The second is the plane-parallel Poisson equation for disk self-gravity. The third balances the vertical increase of radiative flux against local viscous heating in a Keplerian alpha disk, and the fourth is the optically thick radiative diffusion law.

With $P/\rho \sim c_s^2$, $\Sigma \sim \rho H$, and all terms in hydrostatic balance comparable,The first comparison gives $H\sim c_s/\Omega$; inserting it into the second gives $\Omega c_s\sim G\Sigma$. Henceup to the order-one constants deliberately omitted by the scaling argument.

Along a ray, the radiative transfer equation has the formal solution of the radiative transfer equationfor each isothermal region in local thermodynamic equilibrium, neglecting scattering. For a self-luminous atmosphere with no incident intensity from below at the relevant photosphere, and in the optically thick limit, $I_{\nu ,i}=B_{\nu }(T_i)$.

The observed radiative flux is the projected-disc integralThe inner region occupies projected fraction $f={\sin }^2{\theta }_T$, henceFor finite optical depths, each Planck law function in this formula is replaced by its corresponding emergent intensity above.

Put $\mu =\cos \theta$ and $S=j/\kappa$. Integrating the radiative transfer equation over solid angle givesThe first integral is the radiative flux $F$. A stationary atmosphere with no local energy source has $dF/d\tau =0$, soBy definition, the mean intensity is $J=(4\pi {)}^{-1}\int Id\Omega$. Thereforewhich is radiative equilibrium: the angle-integrated emission and absorption rates are equal.

For blackbody star and planet observed at one common wavelength sufficiently long for the Rayleigh-Jeans law to hold for both,This follows by substituting $B_{\lambda }=2c{k}_{B}T/{\lambda }^4$ in the thermal eclipse depth. A radiative flux ratio alone does not determine the radius without a planetary brightness temperature and a stellar temperature. Spectral peaks interpreted with Wien's displacement law must refer to $F_{\lambda }$, not $F_{\nu }$ or $\lambda F_{\lambda }$.

For the affine angular ansatz $I=A+C\cos \theta$, the radiative flux is $F=4\pi C/3$. Setting the inward hemispheric flux to zero gives $A(0)=2C/3$, but then the formal intensity at $\cos \theta =-1$ is negative when $F>0$. Exact no-incoming radiation and nonnegative intensity cannot both be represented by a nonzero affine polynomial on all directions. The Eddington surface boundary condition therefore specifies approximate moments, rather than an exact emergent angular radiation field.

Radiant flux is the power carried by electromagnetic radiation: radiation energy transferred per unit time, measured in watts, $\mathrm{W}=\mathrm{J}{\mathrm{s}}^{-1}$. Specify what is counted: all radiation emitted by a source, radiation within one beam, or radiation intercepted by a photodetector. Writing $Q(t)$ for the cumulative energy crossing the specified boundary givesThe first expression is an instantaneous derivative; the second is an average over the stated time interval. If $0.020\mathrm{J}$ arrives during $0.010\mathrm{s}$, the mean radiant flux is $2\mathrm{W}$. The photodetector's absorbed power equals this incident radiant flux only if it absorbs all of the incident radiation.

Keeping radiant flux fixed while widening a beam reduces its irradiance, the incident power per receiving area. Radiance additionally resolves the propagation direction and has units $\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}$; an optical camera can collect different radiant fluxes from sources with the same radiance but different apparent solid angles. Luminous flux, measured in lumens, weights the optical spectrum by the human eye's standardized response instead of counting energy equally at every wavelength. Equal radiant fluxes need not give equal luminous fluxes.^[ref]

Radiant power also differs from the optical power of an optical lens, which is reciprocal focal length in dioptres. Likewise, astronomical radiative flux is power per area. Stating the units prevents these shared words from hiding different physical quantities.

Radiative flux is the first angular moment $F_{\nu }=\int I_{\nu }\cos \theta d\Omega$ of specific intensity and measures net radiant power crossing unit area. Integrating over frequency gives bolometric radiative flux. The receiving surface's orientation matters through $\cos \theta$; observed radiative flux from a compact distant source conventionally refers to an area perpendicular to its line of sight.

Radiative forcing is a change in the net downward radiative flux caused by an imposed perturbation, evaluated with a specified adjustment convention before the surface climate responds. State the reference period and whether the diagnostic is at the tropopause or top of the atmosphere. A local seasonal ocean heat flux is not directly interchangeable with global annual radiative forcing.

Radiometry measures the energy and power carried by electromagnetic radiation, including light. It asks how much radiation a source supplies, where it travels, and how much reaches a receiving surface. The measurements are physical rather than weighted by the sensitivity of the human eye: equal radiant fluxes at different wavelengths need not look equally bright. ^[ref].

The five quantities below answer different questions about the same radiation. The International System of Units makes the distinctions visible: a watt is a joule per second, while a steradian measures solid angle. ^[ref].

Resolving wavelength or frequency adds a further “per unit interval”, as in spectral radiance. The notation is consistent with NIST's table of spectral quantities. In particular, radiant flux is power, whereas the existing astronomical radiative flux is power per area. ^[ref].

This bookkeeping makes an optical camera's collected signal calculable from a source's radiance and the geometry of the receiving aperture. It also lets radiation measurements from different instruments be compared. For example, NIST describes separate calibrations of spectral radiance and spectral irradiance for instruments used to observe Earth. ^[ref].

For a plane-parallel atmosphere in radiative equilibrium, constant surface gravity of a star $g$ and opacity $\kappa ={\kappa }_0{P}^{\alpha -1}{T}^{4-4\beta }$, match the stellar structure interior where the radiative flux satisfies $F=\sigma T^4$. The grey atmosphere law then puts the surface at optical depth $\tau =2/3$ with $T^4=2T_0^4$. Integrating hydrostatic equilibrium in optical depth from $P(0)=0$ gives $P^{\alpha }=2\alpha g(T^{4\beta }-T_0^{4\beta })/(3\beta {\kappa }_0{T}_0^4)$. Multiplication by ${\kappa }_0{T}^{4-4\beta }/g$ proves the displayed boundary condition for $\alpha >0$, $\beta \ne 0$. At $\beta =0$ the integral is logarithmic and the continuous limit is $P\kappa /g=4\alpha \log 2/3$. The relation assumes negligible radiation pressure and describes the Eddington closure approximation, rather than exact angular transfer.