Optically thick medium 2026-10-05
An optically thick medium has optical depth along the relevant path. Photons undergo many interactions before escaping. Radiative diffusion applies under suitable near-isotropy conditions, but a large scattering depth alone does not guarantee thermalization when true absorption is weak.
Past exam of the mathematics course of the University of Cambridge 2001 iii Paper 70 2 g Solution Created 2026-10-03 Updated 2026-10-14
For the preceding result reduces to . Using and givesAt fixed and , a lower matter density lowers relative to the flat dust value at positive redshift. The same Electron density then persists for a longer scattering time, increasing the optical depth. Hence the reionization redshift giving unit depth is lower in an open universe.
In the high-redshift matter regime the scaling is at fixed optical depth; for this suggests a redshift near fifty rather than seventy. This comparison holds the baryon density fixed, not the baryon fraction , and is not an estimate of the actual observed reionization epoch.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 63 1 a Solution Created 2026-10-03 Updated 2026-10-14
Efficient convection keeps the interior at nearly uniform specific entropy. For the fully ionized monatomic gas, the heat capacity ratio is , so the adiabatic process gives , with spatially constant. The ideal gas relation then impliesThe constant is spatial, not temporal: loss of specific entropy during contraction changes it. This is a stellar polytrope with polytropic index . We neglect radiation pressure and take the thin photosphere as an atmospheric matching layer rather than a significant fraction of the stellar mass.
For a fixed dimensionless profile of the Lane-Emden equation, the stellar homology relations have the formwhere are positive dimensionless structure constants. More explicitly, write and with and fixed dimensionless profiles. Mass conservation gives . Integrating the hydrostatic pressure support equation, with negligible surface pressure on the interior scale, gives . Both integrals are fixed positive numbers, proving the two central scalings. Substituting these relations into givesThus the same polytropic index fixes throughout this homologous sequence, while the dimensional stellar polytrope constant changes.
At the photosphere, whose optical depth is , identify the temperature with the effective temperature . The ideal gas relation and the interior adiabat give . With the prescribed opacity and stellar surface boundary condition,Eliminating therefore provesUsing the Stefan–Boltzmann law, with , gives the luminosity
The gravitational energy of a stellar polytrope is for . The stellar virial theorem for a monatomic gas gives , and hence the total energy is . With no appreciable stellar nuclear fusion or accretion, energy conservation givesCombining the two expressions for yieldsMultiplication by and integration proves the contraction law with density-linear fourth-power stellar opacity:When the accumulated contraction term dominates the initial-radius term, . Substitution into the central ideal gas relation giveswhich has the stated behavior when the time origin is negligible. This late-time approximation still requires the assumed fully convective star, fixed composition and opacity law to remain applicable.
Past exam of the mathematics course of the University of Cambridge 2004 iii Paper 63 2 Solution Created 2026-10-03 Updated 2026-10-14
Write for the direction cosine and for the radiative transfer source function. Integrating the radiative transfer equation over solid angle givesRadiative equilibrium makes the radiative flux constant, so .
For , the radiation-field moments areConsequently 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 isPutting this equal to zero fixes , hence . Using and the effective temperature definition gives the grey atmosphere lawThe 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 givesTaking and , the grey-atmosphere pressure with power-law opacity is thereforeThe 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 givesThis 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 .
Past exam of the mathematics course of the University of Cambridge 2004 iii Paper 67 3 a Solution Created 2026-10-03 Updated 2026-10-14
In the Maxwell-Boltzmann velocity distribution, each Cartesian component has mean-square velocity . Thus . The original PDF calls this quantity an rms velocity, but the rms speed is its square root. Also, its must mean the dimensionless mean occupation number , since its equilibrium value is a Planck photon distribution; it is not the frequency-integrated photon number density.
Take at a fixed local radiation temperature. Because is proportional to frequency,The hot-electron Kompaneets equation therefore becomesFor traversal length , the Compton y parameter isIt measures electron pressure integrated along the line of sight, or the Thomson scattering optical depth weighted by . In this hot-electron limit the mean fractional energy gain is approximately . The derivation assumes nonrelativistic electrons, Thomson-limit photon energies and . The original PDF uses , as required for ordinary frequency; the converted TeX incorrectly substitutes .
Past exam of the mathematics course of the University of Cambridge 2004 iii Paper 67 3 b Solution Created 2026-10-03 Updated 2026-10-14
Write , with and . Differentiating the Planck photon distribution givesDividing by proves the first-order thermal Sunyaev-Zeldovich effect:In the Rayleigh-Jeans law limit, and , so . In the Wien approximation, the prefactor approaches and the bracket approaches , so . The last limit is an asymptotic description of the first-order coefficient: a perturbative prediction at large still requires .
Hot electrons transfer energy to the cosmic microwave background, removing low-frequency photons and enhancing the high-frequency tail. The first-order distortion crosses zero where , at ; for today's CMB this is near . Photon conservation follows directly from the boundary termIn contrast, integration by parts gives , so radiation gains a fraction of its original energy. The thermal Sunyaev-Zeldovich effect is thus a spectral distortion, not simply a new blackbody temperature.
A galaxy cluster with bulk peculiar velocity along the line of sight additionally produces the kinetic Sunyaev-Zeldovich effect. Define for recession and . To first order in velocity and optical depth, : receding gas gives a decrement and approaching gas an increment. Its occupation number perturbation is , the spectrum of a small blackbody temperature change. It remains nonzero at the thermal null and therefore complicates measurements near that frequency.
Past exam of the mathematics course of the University of Cambridge 2005 iii Paper 74 2 c Solution Created 2026-10-03 Updated 2026-10-14
Typical weak Lyman-alpha forest absorbers have neutral hydrogen column densities roughly –, with stronger forest absorption extending towards . Common Doppler broadening parameters are a few tens of , consistent with gas at approximately K together with bulk velocity gradients. The thermal contribution isA line width is consequently a temperature upper bound if all other broadening is neglected, not a unique thermometer. The gas commonly lies in diffuse sheets and filaments at modest density contrast; only the strongest systems need substantial self-shielding or a galactic environment.
For optically thin gas illuminated by the cosmic ionizing background, photoionization equilibrium requireswhere the last approximation uses almost complete hydrogen ionization. For example, , and give . A kpc path then has despite its tiny neutral fraction. This illustrates the Lyman-alpha forest photoionization scaling: low density makes the two-body radiative recombination rate weak relative to the one-body photoionization rate, but a long path and a strong resonance still produce measurable absorption.
Independent information is important because a column density alone is not a volume density. Similar absorption in close pairs of quasars indicates substantial transverse coherence, while an ionizing optical depth shows that these weak systems do not shield themselves from the background. The continuum-shielding threshold is . Together with the observed line widths and ionization balance, these facts favor extended, low-density, highly ionized gas over small cold neutral clouds. The physical interpretation and characteristic scales are consistent with Rauch's forest analysis and the measured column and line-width distributions.
Past exam of the mathematics course of the University of Cambridge 2006 iii Paper 71 2 Solution Created 2026-10-03 Updated 2026-10-14
Consider a small outward displacement of a fluid element. It adjusts to the ambient pressure while retaining its specific entropy and stellar composition, so an adiabatic process gives . Its excess mass density over the new surroundings is, to first order,Restoring buoyancy requires this bracket to be positive. With , the Schwarzschild criterion is thereforeEquality is neutral stability. For uniform mean molecular weight, the ideal gas law gives . Since hydrostatic equilibrium pressure decreases outward, the inequality is equivalent to . Dividing the radiative diffusion in a star and hydrostatic equilibrium stellar structure equations givesThis is the stellar radiative temperature gradient test for stellar convective stability.
In the thin upper stellar atmosphere take , and use gas pressure dominance. The prescribed opacity becomes . Dividing the two structure equations then yieldsAt zero optical depth, the atmospheric boundary has and negligible gas pressure. Integrating from that boundary, with , givesTo locate the onset, put . Logarithmic derivative gives , henceThis starts at zero at the outer boundary and increases inward. It first reaches atThe corresponding optical depth is , so the onset lies inside the stellar atmosphere. This is grey-atmosphere convection onset with thirteenth-power opacity.
Let the interior of the fully convective star have , with spatially constant. At its boundary is fixed, so the atmospheric pressure formula givesThe Stefan–Boltzmann law, with , consequently impliesThe ideal gas equation transforms the interior relation intoFor completeness, the mass-radius scaling follows by taking and with . Hydrostatic equilibrium and mass conservation reduce to the Lane-Emden equation , with , . Its fixed dimensionless profile gives and . Eliminating gives , or . Thus at fixed stellar composition. Equating this with the atmospheric expression gives , and henceThis Hayashi relation with thirteenth-power opacity concerns the interior of a fully convective star beneath a thin stellar atmosphere. The coefficient may vary between stars; spatially constant specific entropy does not mean equal specific entropy across the whole sequence.
Past exam of the mathematics course of the University of Cambridge 2008 iii Paper 70 2 Solution Created 2026-10-03 Updated 2026-10-14
Let be the gravitational energy. Multiply the stellar hydrostatic equation by and integrate through the star. Integration by parts and zero surface pressure givewhile its hydrostatic right-hand side is . Hence the stellar virial theorem is
There is a genuine reference-sheet error to resolve. A monatomic ideal gas with adiabatic exponent has specific internal energy of an ideal gas , not the printed on PDF page 2. The latter would give total internal energy and zero total energy, inconsistent with radiating contraction. The following derivation uses the consistent monatomic value.
A fully convective star is approximately isentropic and has an adiabatic stellar polytrope with index . The given polytropic gravitational-energy formula yields . The virial relation gives total internal energy , so total energy is . At constant mass with no nuclear heating, Kelvin-Helmholtz contraction obeys , whenceHere denotes physical time.
In the plane-parallel atmosphere, the inward optical depth satisfies , so hydrostatic equilibrium in optical depth gives . Neglecting radiation pressure, and the specified opacity becomes . ThereforeThe grey atmosphere relation and the boundary integrate toTo apply the Schwarzschild criterion, put . Logarithmic differentiation givesThis increases inward, starting from zero. Marginal convection has , so . Since , the grey-atmosphere convection onset with seventeenth-power opacity isThe purely radiative atmosphere is used only up to this boundary.
The convective relation is spatially uniform in its entropy coefficient at each instant; can change as the star radiates and contracts. At the convection boundary is a fixed number, so at fixed composition and opacity normalization. ThusOn the other hand, the central scalings for the stellar polytrope give and . Hence, at fixed composition,Equating the interior and atmosphere expressions yields . Applying the Stefan–Boltzmann law then proves the Hayashi contraction with seventeenth-power opacity relation
At fixed mass write , with constant and . The contraction equation is , where . Direct integration givesMultiplying through by and substituting the current luminosity givesFor the stated very large initial radius, or after that initial-radius term becomes negligible,This is the luminosity-age relation within the contraction model, rather than an assertion that a physical star begins at literally infinite radius.
Past exam of the mathematics course of the University of Cambridge 2009 iii Paper 66 3 iv Solution Created 2026-10-03 Updated 2026-10-14
The quasar proximity effect is weaker Lyman-alpha forest absorption near a quasar's emission redshift. Extra photons from the quasar increase the photoionization rate and reduce the neutral fraction without necessarily reducing the total gas density.
For isotropic specific intensity and a steady, isotropically radiating quasar with spectral luminosity , the background and local rates areHere is a proper separation small enough that attenuation and expansion across it can initially be neglected, and is the hydrogen photoionization cross-section. At fixed density and temperature, photoionization equilibrium makes the neutral column density, and hence a given line's optical depth, scale asThe observed quasar luminosity supplies ; the distance at which the absorption changes significantly estimates where . Adopting a background spectral shape then converts the inferred rate to its specific intensity, often quoted near the hydrogen ionization energy threshold.
Important uncertainties include an enhanced gas density contrast around the host, the unobserved ionizing spectrum, anisotropic emission, source lifetime and variability, the systemic redshift and continuum placement, attenuation between source and absorber, and fluctuations between lines of sight. Higher gas density enhances recombination and can partially hide the effect, biasing a simple background estimate upward. Dall'Aglio, Wisotzki and Worseck's 2008 measurements illustrate the observed sightline scatter and the roles of variability and clustering.
Past exam of the mathematics course of the University of Cambridge 2015 iii Paper 59 1 a Solution Created 2026-10-03 Updated 2026-10-06
Take the optical depth to increase inward and let denote an outward ray. For thermal absorption and emission in local thermodynamic equilibrium, the radiative transfer source function is the Planck function, so the radiative transfer equation isThis assumes negligible scattering, or a source function genuinely thermalized to ; LTE by itself does not turn an arbitrary scattering source into a Planck function. Multiply by , integrate between the two boundaries, and solve for the outward specific intensity:The bottom boundary value is the quantity denoted in the supplied notation; it is not an optical-depth-zero boundary. The first term is attenuated incident radiation, and the second is emission from the intervening layers. This is the formal solution of the radiative transfer equation.
For an isothermal atmosphere at temperature , the Planck function is constant and the integral can be evaluated:For a bounded bottom intensity and , in every outgoing direction, and . Thus the emergent spectrum is a blackbody spectrum. The blackbody limit of an isothermal atmosphere does not require a temperature gradient: it follows from complete thermalization in an optically thick medium.
Past exam of the mathematics course of the University of Cambridge 2015 iii Paper 59 2 a Solution Created 2026-10-03 Updated 2026-10-14
Let increase outward and define the inward optical depth by . With constant outward internal radiative flux , radiative diffusion givesThe constant is a boundary condition; diffusion alone does not determine it. For an unirradiated grey atmosphere with the Eddington closure approximation and Eddington surface boundary condition, , soThe internal effective temperature of a planet is defined by its intrinsic cooling flux, not by its incident stellar heating. Since and , temperature decreases outward. Towards the thin upper layers, ; if the density and optical-depth gradient vanish there, . This is the upper nearly isothermal atmosphere. The diffusion approximation itself fails at small optical depth: the grey boundary closure supplies the approximate continuation. Small irradiation perturbs this intrinsic-flux solution.
Young, self-luminous gas giants at wide orbital separations can have intrinsic cooling dominate their photospheric budget. They are favorable for exoplanet direct imaging, especially in the infrared, where their own thermal radiation is easier to separate from the host light. A strongly irradiated hot Jupiter instead has a large stable outer radiative region set mainly by stellar heating; it can be nearly isothermal over a broad pressure range or develop an atmospheric thermal inversion if stellar light is absorbed sufficiently high. Inversion is not inevitable for every strongly irradiated planet. Deep convection begins below its radiative-convective boundary.
Intrinsic grey-atmosphere cooling profile compared with illustrative irradiated hot-Jupiter profiles
. Details. Past exam of the mathematics course of the University of Cambridge 2016 iii Paper 315 1 a Solution Created 2026-10-03 Updated 2026-10-14
Assume , , uniform disk-averaged thermal emission, a nearly blackbody stellar continuum, negligible reflected light at , and a flat response over the narrow band. If the eclipse depth is normalized to the out-of-eclipse star-plus-planet flux, the planet/star flux ratio is ; using changes the estimate by less than one percent. The exoplanet secondary eclipse then givesUsing the Planck law at gives , or about at the precision justified by the radius assumption. The Rayleigh-Jeans law instead gives , a useful rough check. This is a brightness temperature, not the bolometric planetary equilibrium temperature. With zero Bond albedo and global redistribution, the latter would be only about at ; a thermal window can sample hotter, deeper gas.
A continuum photon normally emerges from an optical depth of order unity. In hydrostatic equilibrium,The emission-pressure degeneracy of a brightness temperature means that pressure cannot be uniquely recovered from the measured brightness temperature without and an opacity model. For an explicit nominal estimate, take and an effective continuum mass absorption coefficient . Then , or approximately one bar. In a hydrogen-rich atmosphere, collision-induced absorption by transient pairs supplies continuum opacity even without ordinary molecular lines; its density dependence makes the constant-opacity calculation an illustrative estimate. With a normal deep atmospheric pressure-temperature profile, temperature increases at larger pressures, eventually following an interior adiabatic temperature gradient.
Two mechanisms for are wavelength-dependent sampling of a vertically varying temperature, and wavelength-dependent weighting of a horizontally nonuniform dayside. In the first, carbon monoxide-rich opacity near can sample cooler gas around bar while the continuum samples hotter gas around one bar. In the second, hot and cool areas contribute different proportions of the total flux in different bands because the Planck function depends nonlinearly on temperature. The vertical mechanism alone produces the illustrated pair of temperatures; a single isothermal, horizontally uniform blackbody would have the same brightness temperature in both bands.
Illustrative pressure-temperature profile with different infrared brightness-temperature contribution depths
. Details. If the J-band observation is thermal emission, its low-opacity window normally samples at least as deeply as the continuum at . On the illustrated non-inverted profile, . The ordering is conditional: molecular-line absence alone does not specify the relative continuum opacities, and an inversion or reflected-light contamination can change the inference. The exoplanet emission spectrum measures weighted contribution regions rather than three exact thermometers at known pressures.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 315 1 b Solution Created 2026-10-03 Updated 2026-10-06
Take to be the radius of an opaque planetary disc and the geometric thickness of the model atmosphere, not necessarily one atmospheric scale height. Assume a uniform stellar specific intensity, a fully projected non-grazing transit, negligible planetary emission in the measured band, and no scattering or refraction returning light to the beam. The cylindrical approximation assigns the same slant optical depth to every ray through the annulus.
The opaque disc blocks area . The annulus has projected area , and the radiative transfer equation transmits fraction through it. Its blocked fraction is therefore . Dividing the missing light by the unobscured stellar-disc light gives the exoplanet transmission spectrumFor , the annulus model for transmission spectroscopy becomesThe optically thin excess is ; the optically thick limit is the area ratio before the thin-annulus approximation. For a realistic atmosphere, the slant optical depth varies with the ray's impact parameter , giving insteadLimb darkening replaces the simple area weighting by the local stellar specific intensity. The constant-depth model is therefore an explicit geometric approximation, not the slant-depth law of a spherical hydrostatic atmosphere.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 315 2 c ii Solution Created 2026-10-03 Updated 2026-10-06
The two measurements sample different atmospheric regions and different radiative transfer geometries. An exoplanet emission spectrum probes the dayside; an exoplanet transmission spectrum probes the morning and evening limbs of the day-night terminator. Several alternatives to an exoplanet cloud deck can therefore suppress the transmission feature.
- A cooler terminator has a smaller atmospheric scale height , so the same change in slant optical depth produces a smaller transit-area change. A larger mean molecular weight has the same geometric effect. The assumption that the nightside temperature sets both limbs may itself be wrong.
- Genuine chemical heterogeneity can reduce terminator carbon monoxide. At lower temperature, thermochemical equilibrium can favour methane through , if local reaction times are sufficiently short. Atmospheric photochemistry, elemental abundance differences, and horizontal chemical quenching can modify this simple picture.
- The feature depends on a band-to-continuum opacity ratio, not abundance alone. A continuum enhanced by other gas absorption or a partly covered heterogeneous limb can mute the feature; the two-limb transmission spectrum illustrates this averaging. Stellar surface heterogeneity and measurement systematics can also bias an inferred feature amplitude.
- A strong CO emission feature can be amplified by an atmospheric thermal inversion, because the line forms in hotter gas than the continuum. It is consequently not an unambiguous abundance measurement without a simultaneous atmospheric pressure-temperature profile fit.
Temperature, mean molecular weight, limb chemistry, continuum opacity and the dayside temperature gradient are alternatives or degeneracies to test. In particular, the three-scale-height amplitude alone cannot establish different CO abundances between the hemispheres.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 315 4 d Solution Created 2026-10-03 Updated 2026-10-06
Let an isothermal non-scattering atmospheric layer of temperature lie above an optically thick continuum source with brightness temperature . Assume local thermodynamic equilibrium, so the layer emits with Planck function . Its vertical optical depth is ; a ray with direction cosine has slant depth .
The formal solution of the radiative transfer equation givesThe two terms are attenuated background and layer emission. Compare a line and nearby continuum at the same wavelength to adequate approximation, with . The one-layer spectral contrast isThe second bracket is positive, and increases with temperature. Therefore a cooler upper layer produces absorption, a hotter upper layer produces emission, and an isothermal source-plus-layer produces no line contrast:This is why higher opacity sampling a cooler altitude creates absorption in an outward-cooling atmosphere, whereas an atmospheric thermal inversion can create emission. If both line and continuum are very optically thick in the same layer, their contrast tends to vanish even with a different deeper temperature. Scattering, non-LTE excitation, and nonuniform layers require a more general radiative transfer treatment.
Past exam of the mathematics course of the University of Cambridge 2017 iii Paper 317 1 Solution Created 2026-10-03 Updated 2026-10-06
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 , the enclosed mass is , the luminosity is , and the outward radiative flux defines the effective temperature through . With no external illumination the incoming specific intensity satisfies for . A deeper boundary is matched to the nearly isotropic radiative diffusion field. A grey atmosphere commonly matches near ; under the Eddington surface boundary condition this gives . For roughly constant gravity and opacity, the total pressure rises by relative to its outer boundary value; if radiation supplies appreciable support, the gas-pressure gradient uses the effective gravity instead.
Take increasing outward and the outward direction cosine. If is emission per unit path length and solid angle, the radiative transfer equation isDefine inward-increasing optical depth by and the source function by . Then . This states the emission-coefficient convention; if is defined per unit mass instead, the emission term is .
For frequency-integrated intensity, introduce the radiation-field momentsThe first angular approximation writes . Its odd part does not contribute to , 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 and the incoming hemisphere to be dark. Its moments give and . Integrating the moment equations with the closure then gives ; in thermal equilibrium, gives . This derives the usual approximate photospheric matching condition rather than imposing a zero surface temperature. For conservative coherent isotropic scattering, . Integrating the transfer equation gives and , 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 , take and defineThe continuous zero limits are , , . Direct integration yieldsFor a decaying nontrivial correction, and . Exact constant flux then requires . But for every nonzero real , so forces . 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 ,Comparing the constant and linear coefficients gives and . Hence regular real coefficients require , 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, and supplies an explicit counterexample to the claimed exact flux constancy, because .
Several requested conclusions remain useful as formal approximations, provided they are identified as such. If , the large-depth exponential vanishes, giving , hence for the asymptotic physical flux. It also gives . Since , neglecting second-order angular errors yields the formal relation . However , 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 used above, the Hopf function of a grey atmosphere is defined byThe PDF's coefficient would require redefining its as physical flux divided by ; it cannot simultaneously use the earlier 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 , 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, for . Nor does the simple linear-plus-exponential ansatz satisfy the exact no-incoming-radiation condition for every negative . Thus is a clearly labelled approximate disk-centre-to-limb result, not an exact consequence of inconsistent premises.

