Consider a planar circular aperture of diameter in constant refractive index . Suppose every position in its area is filled with directions inside the same rotationally symmetric cone about the unit normal, with half-angle . Using spherical coordinates, the exact etendue is
where is the appropriate side's numerical aperture. The ordinary cone solid angle is ; its cosine-weighted value is . A full forward hemisphere therefore gives , rather than . [ref]
For an original model example, choose , , and . Then
If the radiance is uniform across this filled bundle and equals , its radiant flux is . Neither the geometry nor the chosen radiance is an experimental measurement.
A detector collects a range of solid angles, not a single direction. If its nonnegative angular response is , define
Integrating the polarized beta-decay angular distribution gives a count proportional to . For a uniform forward cone with half-angle , put . Direct integration gives
The opposite cone has mean cosine . A hemisphere gives ; a narrow cone tends to . Hence a spin-reversal cross-ratio for beta decay measures , and an acceptance correction is needed to obtain .
With unequal angular responses, let the two laboratory mean cosines be and , including their signs. Then the actual expectation is
which reduces to the simple square only when . The cancellation of stable multiplicative detection efficiencies survives, but their angular shapes still matter. An energy-integrated response likewise averages with the accepted spectrum; a mean speed cannot in general be substituted independently of acceptance.
Additive backgrounds, unequal polarization magnitudes or detection efficiencies that change on reversal do not have the assumed factorization. For example, a reversal-dependent efficiency cross-ratio multiplies by . Background subtraction also creates statistical complications: positive expected signal counts are assumed in this derivation, whereas finite samples with Poisson distributions can contain zero counts. Taking square roots of negative background-subtracted counts is not a valid estimator. These limits explain why reversal and independent preparation checks belong together.
Focal ratio Created 2026-10-06 Updated 2026-10-18
The focal ratio, or f-number, of a converging optical instrument focused at infinity is the dimensionless quantity
where is effective focal length and is the on-axis entrance pupil diameter. The notation f/4 means , not in a particular unit. A focal length with gives f/4. Closing the aperture stop until gives f/8 while retaining the same focal length. The mechanical opening need not equal . [ref].
For a thin lens in air with a circular aperture stop at the thin lens, the converging bundle from a distant on-axis point has half-angle satisfying
The last step uses the small-angle approximation. For a compound optical lens, the cone is set by the image-side optical pupil and its distance from the receiving plane; the simple triangle is not an exact high-angle rule. At finite object distance a thin lens instead gives , by its imaging equation and transverse magnification, so its working cone corresponds to . Nominal focal ratio therefore does not alone specify close-range image irradiance. [ref].
Two useful scalings follow under stated conditions. For a fixed uniform incident beam and unchanged transmission, collected radiant flux scales as . For an extended source with fixed radiance, the central image irradiance scales approximately as in the paraxial approximation: the solid angle of the cone scales as . For an ideal uniformly illuminated unobscured circular entrance pupil, the first dark ring of the point spread function has angular radius and focal plane radius
Here is wavelength. These are first-zero radii, not diameters or a universal rule for resolving arbitrary objects. A proof outline evaluates the circular aperture's Fraunhofer diffraction using its Fourier transform, then uses the first zero of the resulting Bessel function. [ref].
At , the example's f/4 and f/8 first-ring radii are and . Closing the aperture stop halves , quarters collected radiant flux and doubles this diffraction radius. It can also reduce optical aberrations by excluding outer light rays; actual sharpness reflects both effects. This explains why a smaller aperture need not produce a sharper image. [ref].
Figure 1.
Aperture and focal ratio: a common optical focus, reduced collection and increased diffraction radius
. Original analytical illustration. Left: equal-scale thin lens geometry with and or . Right: collected radiant flux relative to f/4 is , and first-ring radius relative to f/4 is . The wavelength is ; the entrance pupil is circular, uniformly illuminated and unobscured. These are calculated quantities, not measurements. Details.
A proposed optical bench check puts a fixed collimator before a converging thin lens and a photodetector large enough to collect the entire image behind it. Compare two centred aperture stops of known diameters while fixing source output, detector gain and integration time. The integrated signal tests the collection law; a spatially resolved optical camera at the focal plane can separately measure the point spread function. Repeat a dark reading and the first aperture stop, and check that neither the incoming beam nor detector clips either bundle. A resolved pinhole source, imperfect collimation, optical aberrations and detector pixel size can dominate the apparent spot width; a single peak reading does not measure total collected radiant flux.
For a narrow angular power pattern about an equatorial beam centre, the solid angle measure is locally . Extending the small lobe to the plane gives its beam solid angle . Both sides of the azimuthal lobe must be included. A lossless antenna then has peak antenna gain and peak antenna effective area .
For tracers sampled uniformly on a Galactocentric shell of radius , with the Sun at radius , the correction in the heliocentric projection of a spherical velocity dispersion is
The limiting value at is ; it tends to as and to at large radius. Uniform heliocentric solid angle at fixed instead gives exactly. Specify the selection measure when applying the correction.
For box-normalized plane waves in volume , a potential has matrix element . The final-state density per energy and solid angle is . Dividing the golden-rule rate by incident flux gives . Elastic scattering has ; the arbitrary box volume cancels.
Lumen (unit) 2026-10-14
The lumen is the International System of Units unit of luminous flux: . A source with constant luminous intensity within a cone of solid angle supplies into that cone. A source with the same luminous intensity in every direction supplies over a full sphere. [ref]
Luminous intensity 2026-10-14
Luminous intensity is luminous flux per outgoing solid angle, measured in candelas. It is the visually weighted counterpart of radiant intensity. The same total luminous flux concentrated into a narrower cone can give a greater luminous intensity within that cone.[ref]
How much does a small directional uncertainty help? Consider an original idealized model: a source direction fixed in advance, uniform background directions over the full sphere, and circular acceptance radius . The accepted solid angle and single-event background probability are
The derivation integrates the area element on a unit sphere. For small in radians, : halving the radius quarters the accidental acceptance, to leading order. At acceptance vanishes; at it covers the full sphere. The solid angle has units of steradians; is dimensionless.
For , take a deliberately crude point spread function: stipulate that every signal direction lies uniformly inside that cap. The signal and background probability density functions per steradian are inside the cap and zero outside, and everywhere. Both integrate to one. For an event inside the cap their likelihood ratio is
This ratio compares two specified directional models; it is not the probability that the event came from the source. Such a conditional probability would additionally require prior source odds and the other selection information. For fixed source fraction among selected events, a normalized mixture distribution has density . For independent events the likelihood is the product of those densities; its ratio to background alone is
Outside-cap events contribute factors rather than being silently discarded. At , ; for a full-sky cap, and every gives . A directional model that cannot distinguish source from background carries no directional evidence.
Figure 1.
Chosen angular acceptance, directional likelihood and repeated background opportunities
. Original exact spherical-cap calculation and uniform-cap signal model. Radii and event counts are chosen illustrations, not IceCube angular resolutions or discovery data. The fixed source direction and independent uniform background assumptions are essential. Details.
As an original numerical example, radii and give and , respectively: nearly a hundredfold difference. For a fixed number of independent background events, the binomial distribution gives
With chosen , those radii give approximately and . A small single-event acceptance can become a common coincidence after many opportunities. This probability concerns a fixed cap, rather than a scan for the most favorable cap anywhere. For independent searches with common local tail probability , the corresponding global tail probability is ; overlapping sky searches generally violate that independence. Actual source tests must include directional response, nonuniform exposure, energy information and the look-elsewhere effect in their null calibration.
Nuclear elastic electron scattering cross-section Created 2026-10-18 Updated 2026-10-20
For an ultrarelativistic electron scattered elastically by a static, spherical, spin-zero atomic nucleus, a plane-wave Born approximation retains a single electromagnetic interaction. The Fourier transform of its Coulomb potential is proportional to . Relative to the same point charge, the scattering amplitude is therefore multiplied by and the differential scattering cross-section by :
Here is the incident total electron energy, the scattering angle, a solid angle element, and the fine-structure constant. The cosine factor already accounts for the projectile's spin; a spinless projectile would give a different point-charge factor. The formula neglects nuclear recoil, electron rest-mass corrections, radiative corrections and Coulomb distortion in electron scattering. The latter can be substantial at large : electron waves then require treatment in the target's Coulomb potential, rather than undistorted plane waves.
Nonzero target spin can introduce magnetic and additional electric contributions, so division by this point-charge expression does not generally isolate one charge form factor for a proton or deuteron. Even for spin zero, finite energy resolution, recoil and Coulomb distortion in electron scattering can change a measured minimum. A zero of the uniform-sphere model is a clear illustration of amplitude cancellation, not a universal prediction of a zero experimental counting rate.
Observed flux of an unresolved spherical emitter Created 2026-10-06 Updated 2026-10-14
For a uniformly emitting spherical body of radius at Euclidean distance , integration of its specific intensity over the apparent solid angle gives . For a blackbody, and the bolometric result is . A grey emissivity factor multiplies this expression; isothermality alone does not imply blackbody emission.
Fermi's golden rule converts a weak transition matrix element and the available density of final states into a transition rate. For a static perturbation and a continuum,
When the matrix element is nearly constant on the energy shell, this becomes . The energy delta function is to be integrated, rather than interpreted as a finite rate to one exactly resonant discrete level. The observation time must resolve a smooth final-state continuum while remaining short enough that first-order depletion is small. A periodically driven perturbation instead gives the corresponding absorption and emission energy shells.
For elastic scattering from a weak localized potential , use box-normalized free states . Then
The number of final states is . With , the density per energy and solid angle is . Dividing the resulting golden-rule rate by the incident single-particle flux gives
This is the golden-rule derivation of the Born cross-section. For the elastic problem , so the first Born scattering amplitude and cross-section are
The box volume cancels, as a physical cross-section must. The minus sign is fixed by the outgoing-wave scattering convention; the golden-rule rate alone determines only the modulus.
For example, has , hence , with . Replacing the exact scattering state by the incident plane wave is the essential Born approximation. It can fail for strong scattering, low-energy resonances and situations requiring repeated scattering; a small perturbation and a suitable continuum are separate conditions.[ref]
For elastic scattering by the Yukawa potential, take and use the first Born approximation in the plane-wave normalization
Choose the polar axis along the momentum transfer . The angular integral is
The radial volume element cancels the remaining powers of , leaving
For real , the differential scattering cross-section is therefore
This also has a finite forward limit when . The Born amplitude for a Yukawa potential provides a screened route to the Coulomb potential by taking at fixed .
For an alpha particle of charge scattering from a heavy nucleus of charge , set . Neglect recoil, so that the reduced mass is approximately the projectile mass, , and elastic geometry gives
The Coulomb-limit scattering amplitude becomes
which gives the Rutherford differential cross-section
The unscreened Coulomb potential has long range, so taking this limit does not justify the Born approximation for its full scattering state, particularly for a strongly charged target. The exact Rutherford scattering amplitude differs by a phase and has the same modulus. Thus the resulting differential scattering cross-section remains applicable to the ideal nonrelativistic Coulomb problem beyond the weak-potential regime.
To relate the detector counts to the atomic number, first compute the number of target atoms per illuminated area. With molar mass , the number density integrated through the foil is
The detector's transverse area is already its area normal to the scattered ray. Its small solid angle is
There is no extra factor: the foil is traversed normally and the detector area is transverse to its line of sight.
For a thin foil with negligible multiple scattering and energy loss, the detected fraction is . Thus the measurements imply
Using in the Rutherford differential cross-section gives
The stated measurements do not reproduce gold's atomic number 79 in this model. The values and units above agree with the original PDF. Substituting instead would predict with the other inputs unchanged, rather than . No stated correction accounts for that discrepancy; forcing the answer to 79 would require changing the data or adding an experimental correction not specified in the problem.
Choose the baseline vector from the first antenna to the second. Apart from a sign convention for receiver delay, the geometric delay is , and the phase difference is . For equal received source powers , adding the voltages and averaging over many wave periods gives
Uncorrelated receiver noise adds a mean level but not the source's coherent cross term. For an east-west baseline and one choice of hour-angle sign, . Earth rotation therefore sweeps the phase; locally the fringe frequency is . Tracking the antenna pointing keeps the source in the primary beam. Electronic fringe stopping is a separate operation: it can remove this phase evolution, making a source at the phase centre constant in the corrected complex visibility. The oscillating power above is the geometrical fringe before that correction.
Figure 1.
Projected interferometer baseline, geometric delay and illustrative fringes before electronic fringe stopping
. Details.
For a small angular displacement along the projected baseline, . One full fringe therefore spans . Here . This is the angular period of the measured Fourier component and the characteristic resolution scale when such long baselines are synthesized; precise beam width and image ambiguity also depend on baseline coverage and weighting. A single pair by itself does not give an unambiguous complete image.
Writing , the expansion immediately gives the Rayleigh-Jeans law,
The factor two counts the two independent modes of polarization of waves. A small brightness-temperature perturbation thus has in this limit.
Hot intracluster electrons inverse-Compton scatter the microwave-background photons. The scattering optical depth is and the characteristic thermal fractional energy transfer is of order . The integrated distortion is therefore controlled by
Scattering conserves photon number while shifting photons toward higher frequencies, giving a low-frequency decrement and a high-frequency increment. Thus the printed estimate correctly gives a negative change of order . Its numerical coefficient, however, is not the standard Rayleigh-Jeans coefficient. For nonrelativistic electrons and small , the Kompaneets equation gives . At small , , so and
Equivalently the thermodynamic-temperature spectral factor is . The original PDF really prints the estimate with , not ; both numerical conventions are retained below. Relativistic corrections at the supplied are small compared with this factor-of-two difference.[ref]
The uniform cube has a constant line-of-sight depth . Using the supplied cross-section rather than replacing it with a more precise value gives
Adopt , explicitly an external measured constant, not a temperature supplied in the question. The values , and set the physical units. Hence the printed relation gives ; the standard Rayleigh-Jeans result gives . Keeping the answer as , or twice that, avoids needing a numerical background temperature.[ref]
The face is a square, so its solid angle is , not the area of a circular six-arcminute aperture. At ,
Thus using the printed estimate, or using the standard factor two. The sign denotes missing brightness relative to the surrounding microwave background. At this frequency ; the exact nonrelativistic spectral intensity differs from the strict Rayleigh-Jeans approximation by about one percent. This precision does not affect the estimates.
For the array noise, the system temperature sets the receiver noise scale; let be the effective area of each dish, and let describe correlation losses. No efficiencies or polarization of waves count are supplied, so first take and a single independent polarization of waves. An unpolarized point source produces antenna temperature , and the system-equivalent flux density is . The interferometric radiometer equation gives one baseline an rms flux error . There are cross-correlation baselines; averaging their independent thermal noise for a point source gives
The ideal single-dish effective area is , giving and
Combining two independent parallel-hand polarization of waves measurements gives . Nonunit efficiencies increase these estimates by . Here signal equal to noise means a weak-source signal-to-noise ratio of one after integration; it does not mean equating the source to the entire unaveraged receiver-noise power.[ref]
At the band centre, . Taking a representative maximum projected baseline of gives arcminutes. A corner-to-corner baseline of would give about arcminutes before projection and weighting. The exact antenna locations and source declination are not supplied, so neither a unique synthesized beam shape nor an exact FWHM is determined. The individual-dish primary-beam scale, roughly arcminutes, is the field of view, not the synthesized resolution.
To extract the cluster signal, retain frequency-resolved complex visibilities and fit an extended SZ template together with contaminating compact sources. For the ideal uniform square centred at the phase centre and with a nearly constant primary beam across it, the Fourier transform is
This is why the integrated decrement must not be treated as a point-source flux on every baseline. The first axial visibility null occurs at . Short baselines carry most of the large-scale decrement; the longest baselines are useful for separating unresolved positive radio sources. natural visibility weighting or a suitable taper gives a clearer extended-source image than emphasizing the longest baselines, while a visibility-domain matched filter uses all samples with weights proportional to template amplitude divided by noise variance. Its amplitude estimate is for a real template amplitude. A joint fit propagates radio-source uncertainty instead of subtracting exact, noise-free source models.
Figure 2.
Normalized visibilities of the six-arcminute square cluster and their suppression on long baselines
. Details.
Earth-rotation and multifrequency synthesis improve Fourier coverage, but each channel must use its own baseline in wavelengths, primary beam and SZ spectrum. Use the longer baselines and spectra to identify variable point sources, flag bad data, and test the decrement with time, antenna and frequency splits. Primary-beam and absolute-flux calibration must be included. Large-scale microwave-background fluctuations and residual ground or atmospheric signals can remain even below the thermal-noise estimate. The narrow low-frequency band alone gives limited separation of SZ and primary microwave-background temperature structure. There is no zero-spacing interferometric datum: integrated flux is model dependent unless sufficiently short baselines or additional total-power information are available. Exact cluster detection significance cannot be obtained from the point-source noise alone without the actual projected baseline distribution and weighting.
There are four relevant tree-level Feynman diagrams: -channel photon, Z boson and Higgs boson exchange, and -channel electron-neutrino exchange. The electron and Positron are incoming, and and outgoing. The Higgs diagram is normally tiny because its electron coupling is proportional to ; it should nevertheless be included in the requested list.
Figure 1.
The four tree-level diagrams for electron–positron annihilation into a charged W pair
. Original schematics. Fermion arrows show particle-number flow, which opposes the momentum of an incoming positron. W lines are labelled by charge. Details.
Let . The Higgs-electron coupling is , while the Higgs boson coupling to W bosons is . Multiplying the vertices and scalar propagator gives, up to an irrelevant overall phase,
A scalar propagator has no vector index, the Yukawa interaction supplies the scalar spinor bilinear, and the second vertex contracts the two polarization vectors. We retain the tree propagator and neglect widths, consistently with the stated leading-order calculation.
For unpolarized incoming beams the spin average is , not a sum without normalization. The fermion spin sums and gamma-matrix trace identities give
The terms with one gamma matrix have zero trace; the two-gamma trace gives . Both final polarization states are specified as longitudinal, so no final polarization average is appropriate.
In the centre-of-momentum frame, take and , where and . A pair of longitudinal polarization vectors is
Each obeys and . Their contraction is angle independent:
Consequently
For the two-body scattering cross section, the outgoing-to-incoming momentum ratio is
The source's final sentence interchanges the particle names attached to these two momenta; physical final-state phase space and initial-state flux fix the ratio as written here. Since the Higgs-only spin-averaged expression is isotropic, integrating the solid angle gives :
There is no identical-particle factor: and are distinguishable. The threshold factor makes this contribution vanish when the outgoing momentum vanishes. This is the isolated squared Higgs amplitude for a specified longitudinal pair, not the complete physical cross section: the full result squares the coherent sum of all four amplitudes and includes their interference.
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 satisfy
together with conservation of the element's total abundance. Thus the low-density ion fraction depends primarily on . 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 is
The rate per ion is . Maxwellian detailed balance gives, for ,
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 . A stationary collisional-radiative model solves
where, in the simplest electron-plus-radiative model, with 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 distributes isotropic emission over solid angle. photon intensities omit . A practical source of the necessary atomic rates and population calculations is the CHIANTI reference guide.
Label the two ground-configuration levels and the two indicated excited levels :
The common core is omitted. Si X is silicon with nine electrons removed, hence five bound electrons; the left superscript on or is the spin multiplicity, not another electron-occupancy exponent.
The observed lines correspond to decays and . For approximately uniform plasma at the indicated , form the theoretical energy-intensity ratio
The common abundance, Si X fraction and path length cancel, leaving a ratio governed by level populations. The frequency factor is 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 has only a weak radiative decay to . Ignoring higher-level feeding for this explanation, their balance gives
At low density ; above the critical density of an atomic transition it approaches , approximately two when the gap due to fine structure is small compared with . This saturation can occur while the strong upper-level transitions still remain radiatively dominated.
Let be the total radiative loss rates of levels , and , their line branching fractions. In that upper-level low-density regime,
Thus
The 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 , 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 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 instead
A 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.
Use the luminosity convention in the supplied flux relation: , where is the total spectral luminosity. For a general spectrum , the cosmological spectral flux-density relation gives
In a spatially flat universe, a shell subtending the solid angle contains sources, because is their comoving number density. The cosmological spectral background intensity per unit solid angle is consequently
This is the cosmological background intensity from comoving emissivity specialized to identical, nonevolving sources in an Einstein-de Sitter universe. For the integrand reduces to , or directly . The finite comoving particle horizon therefore gives
For , instead,
The emitted frequency grows with redshift, and the rising spectral luminosity now defeats the redshift dimming. The Einstein-de Sitter spectral-background convergence criterion is ; equality produces a logarithmic divergence. This divergence describes the unbounded power-law, eternal-population model. A finite formation epoch, a high-frequency spectral break or absorption changes that model and can make its background finite.
An isotropic source's bolometric luminosity is distributed over a wavefront of physical area . Each photon loses an energy factor through cosmological redshift, and successive photons arrive with a further rate factor through cosmological time dilation. Therefore the radiative flux and luminosity distance are
These are bolometric formulas; a fixed observed frequency would require the shifted emission spectrum as well.
Normalize the scale factor by . Summing the radiative flux over the source counts on an FLRW past light cone makes the transverse shell-area factors cancel:
Thus is power per effective detector area, summed over the whole sky; an all-sky detector of effective area receives . The corresponding intensity per solid angle is . Restoring the speed of light multiplies these flux integrals by .
Spatially uniform by itself does not specify its time evolution. For a persistent population comoving with the expansion, with neither creation nor destruction, the cosmological continuity equation for source number is . Its conserved comoving number density is , so the cosmological bolometric background from conserved sources becomes . For the Einstein-de Sitter universe,
The last equality uses . The early-time improper integral converges despite the diverging physical source density. If instead the proper number density were maintained at a constant value by continuous source creation, the same geometry would give . For a completely unspecified source history the general integral above is the determined answer; constancy of individual luminosity alone does not imply number conservation.
Use units , and let be the areal radial coordinate of the FRW metric, rather than the radial comoving distance when . On a constant-cosmic time slice the spatial metric tensor is
The determinant of the spatial metric tensor supplies the factor . Integrating over the solid angle gives . An incoming radial null geodesic obeys
Consequently the proper shell volume sampled by an emission-time interval of positive width is . Multiplying by the proper number density proves the source counts on an FLRW past light cone:
Here is evaluated on the observer's past light cone, not along a fixed source worldline. The curvature factor cancels between the spatial volume and the radial light speed; it has not been set to one.
Write for the direction cosine and for the radiative transfer source function. Integrating the radiative transfer equation over solid angle gives
Radiative equilibrium makes the radiative flux constant, so .
For , the radiation-field moments are
Consequently the Eddington closure approximation holds. Substituting into the radiative transfer equation, with , gives . Equality for every requires and , with . The stipulated inward hemispheric radiative flux at the surface is
Putting this equal to zero fixes , hence . Using and the effective temperature definition gives the grey atmosphere law
The positivity limitation of a linear Eddington intensity matters here: the formal surface intensity is negative for . Thus this angular ansatz gives the requested Eddington surface boundary condition for approximate moments; vanishing integrated inward flux is not an exact, nonnegative, pointwise no-incoming radiation field.
Neglecting radiation pressure, the stellar hydrostatic equation and the definition of optical depth give hydrostatic equilibrium in optical depth, . Since , putting gives
Taking and , the grey-atmosphere pressure with power-law opacity is therefore
The pressure formula printed in the PDF is missing the factor . Differentiating its right-hand side would otherwise give . For the correct limit is .
At the specified matching location, the luminosity relation makes , hence , and . Multiplying the corrected pressure relation by gives
This agrees with the PDF's final surface boundary condition for a power-law-opacity grey atmosphere. It also extends to negative wherever the model applies; at , its continuous limit is .