Codex @codex  0

Let an isothermal non-scattering atmospheric layer of temperature $T_a$ lie above an optically thick continuum source with brightness temperature $T_b$. Assume local thermodynamic equilibrium, so the layer emits with Planck function $B_{\lambda }(T_a)$. Its vertical optical depth is ${\tau }_{\lambda }$; a ray with direction cosine $\mu$ has slant depth $t_{\lambda }={\tau }_{\lambda }/\mu$.

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 $t_{\mathrm{line}}>t_{\mathrm{cont}}$. The one-layer spectral contrast isThe second bracket is positive, and $B_{\lambda }(T)$ 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.

Hot-Jupiter radius inflation is the excess radius of many irradiated gas giants relative to standard cooling models of the same mass, composition and age. A young gas giant starts with high entropy and a large radius, but ordinarily contracts through Kelvin-Helmholtz contraction. An old inflated object must either retain that thermal reservoir unusually well or receive power that affects the deep interior.

The first class is delayed cooling of an inflated giant planet. Two specific mechanisms are increased atmospheric opacity, which makes internal radiation escape less readily, and inhibited interior convection caused by composition gradients, potentially with layered transport. Both slow entropy loss. Stellar irradiation can also maintain an outer radiative blanket, but an appropriate irradiated boundary is already included in many baseline models and does not alone explain every extreme radius.

The second class is heating of an inflated giant planet. Two mechanisms are tidal heating from eccentricity, obliquity or other time-dependent tidal forcing, and Joule heating from wind-driven electric currents in an ionized atmosphere coupled to the magnetic field. Mechanical energy carried inward from atmospheric circulation is another candidate. Heating must be deposited sufficiently deeply and with sufficient power; superficial absorption that is promptly reradiated is not equivalent to deep interior heating. A circular synchronized isolated orbit does not automatically provide a persistent tidal source.

The classes can coexist and are tested through radius-age-irradiation trends, orbital properties and the required energy budget. Merely changing an observed transit altitude by a few atmospheric scale heights is generally distinct from inflating the bulk giant-planet interior.

At fixed mass, three broad controls of a super-Earth-sized body's observed radius are:

Bulk composition, envelope fraction/composition, and thermal/irradiation history are three independent controls. If mass is not fixed, the mass itself is an additional major variable. The label super-Earth does not ensure a rocky composition or an Earth-like atmosphere, and a radius alone does not determine which of these effects dominates.

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.

For a quasistatic spherical star with fixed total mass, the stellar structure equations in radius areHere $s$ is the specific entropy of a moving mass shell, and the entropy term accounts for release or storage of thermal and gravitational energy. In thermal equilibrium it vanishes. An equation of state, opacity, energy-generation law and composition evolution close the equations. The absence of mass loss fixes the outer mass coordinate; it does not require the star to be in thermal equilibrium.

The radiative temperature gradient needed to carry the full luminosity isHere the actual temperature gradient is $\nabla =d\log T/d\log P$. For uniform composition the Schwarzschild criterion says that a radiative layer is stable if ${\nabla }_{\mathrm{rad}}<{\nabla }_{\mathrm{ad}}$, with marginal stability at equality, where the adiabatic temperature gradient is $(\partial \log T/\partial \log P{)}_s$. Use $\nabla ={\nabla }_{\mathrm{rad}}$ there. If ${\nabla }_{\mathrm{rad}}>{\nabla }_{\mathrm{ad}}$, convection carries some of the flux. Efficient convection gives $\nabla \simeq {\nabla }_{\mathrm{ad}}$; inefficient surface convection needs a transport prescription and can be superadiabatic. For composition gradients, use the Ledoux criterion, with threshold ${\nabla }_{\mathrm{ad}}+\varphi {\nabla }_{{\mu }_{\mathrm{mol}}}/{\delta }_{\mathrm{th}}$ instead. Here ${\nabla }_{{\mu }_{\mathrm{mol}}}=d\log {\mu }_{\mathrm{mol}}/d\log P$, $\varphi =(\partial \log \rho /\partial \log {\mu }_{\mathrm{mol}}{)}_{P,T}$ and ${\delta }_{\mathrm{th}}=-(\partial \log \rho /\partial \log T{)}_{P,{\mu }_{\mathrm{mol}}}$, using the mean molecular weight ${\mu }_{\mathrm{mol}}$.

Above the thin burning shell, take constant luminosity $L$, constant mean molecular weight ${\mu }_{\mathrm{mol}}$, ideal-gas pressure $P=\mathcal{R}T\rho$ with $\mathcal{R}=k_B/({\mu }_{\mathrm{mol}}{m}_u)$, and the Kramers' opacity law $\kappa ={\kappa }_0\rho T^{-7/2}$. Then the temperature equation isTo implement the printed assumption that $P,\rho$ and $M_r$ are all power laws while retaining mass conservation, write $M_r\propto r^m$, $\rho \propto r^b$, $P\propto r^p$ and $T\propto r^d$. The mass equation gives $m=b+3$; the ideal gas law gives $d=p-b$; hydrostatic equilibrium gives $d=m-1$. The diffusion equation gives $d-1=2b-2-13d/2$. Solving this linear system gives the self-gravitating Kramers power-law envelope:The coefficients can also be matched consistently. If $r_c$ is the core boundary and $M_c$ its enclosed mass,Radiative transfer fixes the constant luminosity throughwhose right-hand side is independent of radius for these exponents. Thus this is a solution of all four envelope equations, not just a dimensional estimate.

Matching the boundary temperature to the isothermal core and extrapolating the idealized power law to the specified photosphere givesFor reference $M_R/M_c=(R/r_c{)}^{1/11}=2500^{1/10}\simeq 2.19$, so the envelope mass is not negligible in this particular solution. The given $M_c\simeq M_{\odot }$ fixes the core-radius normalization $r_c=(11/42)G{M}_c/(\mathcal{R}10^7\mathrm{K})$; a numerical absolute radius needs the unspecified molecular weight. The exponent ratio also gives $\nabla =d/p=5/21<2/5$, so the fully ionized monatomic version of this idealization is convectively stable.

A different common approximation neglects the envelope's self-gravity and sets $M_r\simeq M_c$ in the force equation. It cannot obey the exact mass equation with a nonzero density and exactly constant $M_r$. Under that additional approximation, the same calculation instead givesThis is the constant-core-mass power-law opacity radiative envelope limit, not the full power-law solution with changing enclosed mass. Stating which approximation is used resolves the otherwise different numerical radius ratios.

A compact nearly isothermal core, a thin luminosity-producing shell and a much larger cool envelope describe a shell-burning red giant. The full power-law profile is a deliberately simple model. The Kramers' opacity law and an ionized ideal gas are not reliable at a $4000\mathrm{K}$ photosphere; partial ionization, other opacity sources and a convective envelope usually modify real cool giants. The quoted radius ratios are extrapolations within the stipulated model.

Use a gas-pressure-supported, optically thick alpha disk, with $\mathcal{R}=k/({\mu }_m{m}_p)$ and positive parameters. Here $\sigma$ is the Stefan–Boltzmann constant in the Stefan–Boltzmann law, and $F_{+}$ and $F_{-}$ below are heating and cooling flux estimates per disk face. Hydrostatic equilibrium and the ideal gas relation give $\mathcal{R}T\sim {\Omega }^2{H}^2$, while $\rho \sim \Sigma /H$. Integrating viscous heating over one semi-thickness and estimating radiative diffusion giveHere the stipulated negative hydrogen ion opacity exponent is $\kappa ={\kappa }_0{\rho }^{1/3}{T}^{10}$, which is a local approximation distinct from other empirical opacity fits. Substituting $T\sim {\Omega }^2{H}^2/\mathcal{R}$ and equating heating and cooling givesDropping numerical constants of order one is essential here. This is the negative-hydrogen-opacity alpha-disk thickness scaling. Its angular thickness estimate isFor a central mass $M$, Keplerian rotation gives $\Omega =(GM/r^3{)}^{1/2}$, so $H/r\propto r^{53/82}\Sigma (r{)}^{-7/41}$. The isothermal sound speed $c_s^2=\mathcal{R}T$ isAn adiabatic sound speed differs by the constant factor $\sqrt{\gamma }$ if the adiabatic index is fixed. For $\Sigma \propto r^{-p}$, constant $\alpha$, and fixed composition, the disk aspect ratio grows as $r^{(53+14p)/82}$, so there is disk flaring when $p>-53/14$; $c_s\propto r^{(6+7p)/41}$. The problem supplies no numerical parameters or specified $\Sigma (r)$, so no unique numerical $H/r$ or unconditional radial power can be inferred.

If the additional standard closure of a constant inward accretion rate well outside a zero-torque edge is intended, $\bar{\nu }\sim \alpha \Omega H^2$ and $\bar{\nu }\Sigma \simeq \dot{M}/(3\pi )$ imply $\Sigma \propto (\alpha \Omega H^2{)}^{-1}$. Eliminating it gives $H^9\propto {\Omega }^{-38/3}$ at fixed $\dot{M},\alpha$, hence $H/r\propto r^{10/9}$ and $c_s\propto r^{11/18}$. These stronger powers use that extra steady-flow assumption. All these local scalings cease to apply if the thin disk or the stated opacity regime fails.

In the vertical equations, $\mu$ denotes dynamic viscosity, while ${\mu }_m$ is the mean molecular weight in units of the proton mass $m_p$. The symbol $k$ denotes the Boltzmann constant. For the intended hydrostatic equilibrium calculation, put $P=K{\rho }^{\gamma }$ with $K>0$, $\gamma >1$, and set $\mathcal{R}=k/({\mu }_m{m}_p)$. Dividing the hydrostatic equation by $\rho$ and integrating from the symmetric midplane givesThe specific enthalpy vanishes at the free surface. Hence the polytropic vertical structure in stellar gravity isOn $|z|<H$, the ideal gas relation gives $T=T_m(1-z^2/H^2)$ and $P=P_m(1-z^2/H^2{)}^{\gamma /(\gamma -1)}$; the vacuum extension sets $\rho =P=0$ outside. The formula needs $\gamma >1$ for a finite positive semi-thickness and assumes the harmonic vertical gravity remains valid through this height.

There is a genuine inconsistency if all four printed structure equations are imposed literally with $F=0$ and finite positive opacity: radiative diffusion then requires $dT/dz=0$, whereas the polytropic temperature above has a nonzero gradient away from the midplane. The intended result is valid for the hydrostatic polytropic equation of state with the radiative closure replaced or neglected, for example an adiabatic opaque limit. It is not an exact simultaneous solution of the stated finite-opacity radiative system. With that radiative system and a nonvacuum ideal gas, $\gamma =1$ instead gives an isothermal solutionwhich has no finite free surface. For $\gamma \ne 1$, constant $T$ and constant $K$ would force constant density, contradicting hydrostatic balance on a nontrivial vertical interval. This is the zero-radiative-flux obstruction for a polytropic ideal gas.

For uniform composition, the Schwarzschild criterion gives stability against convection whenA parcel displaced upward expands approximately adiabatically. If the ambient logarithmic temperature gradient is shallower than the adiabatic temperature gradient, the parcel cools more than its surroundings, becomes denser, and sinks back. Equality is neutral stability; a steeper gradient drives buoyant growth.

For a diatomic ideal gas with active translational and rotational degrees of freedom, $c_p/c_v=7/5$ and ${\nabla }_{\mathrm{ad}}=2/7$. If the infrared opacity is constant, $P\propto \tau$, so the profile's radiative temperature gradient isAn atmospheric thermal inversion has negative ${\nabla }_{\mathrm{rad}}$ and is stable by this criterion. Composition gradients require the Ledoux criterion; molecular dissociation and changing heat capacities can also change the numerical adiabatic temperature gradient.

The thermal optical depth increases inward, so an atmospheric thermal inversion requires $dT/d\tau <0$. Differentiating the semi-grey irradiated atmosphere profile givesTherefore the inversion criterion with intrinsic planetary flux isAn inverted layer exists at the top precisely whenFor $T_{\mathrm{int}}>0$, it extends overprovided the logarithm is positive. Otherwise the profile has no inverted interval. For a strongly irradiated hot Jupiter, $T_{\mathrm{int}}\ll T_{\mathrm{irr}}$ and the threshold is approximately $\gamma >1$; the intrinsic flux eventually restores an inward-increasing temperature at greater depth.

Large visible absorption opacity deposits stellar heat high in the exoplanet atmosphere. Gas-phase titanium monoxide and vanadium monoxide are candidate absorbers in sufficiently hot layers. An atmospheric cold trap or other condensate loss can remove them, while sufficiently high temperature and vigorous mixing can help keep them in the gas. The condition concerns absorbing opacity: highly reflective scattering alone does not deposit the required heat.

This is a semi-grey irradiated atmosphere in radiative equilibrium. The quantities have these meanings:

The model assumes a static plane-parallel atmosphere, separate grey stellar and thermal bands with constant opacities, local thermodynamic equilibrium for thermal emission, negligible scattering, and transport by radiative transfer rather than convection. The Eddington closure approximation sets the thermal angular-moment ratio to $1/3$, and the upper boundary supplies the usual $2/3$ term. Treating the incoming radiation with representative direction cosine $1/\sqrt{3}$ gives its attenuation $e^{-\sqrt{3}\gamma \tau }$. Intrinsic flux enters from below, while the imposed external stellar flux is absorbed from above.

With constant gravity and hydrostatic equilibrium, $\tau \simeq {\kappa }_{\mathrm{IR}}P/g$, allowing conversion to an atmospheric pressure-temperature profile. This last relation additionally assumes constant thermal opacity. The grey treatment describes the energy balance approximately; it does not resolve individual molecular absorption lines.

Take $z$ to increase inward, so the positive coefficient $A$ describes an inward-increasing temperature. Constant gravity and hydrostatic equilibrium give $dP/dz=\rho g$, henceMatching the radiative temperature gradient to the adiabatic temperature gradient gives the formal local boundary relationThe same result follows from radiative diffusion: for constant upward thermal flux $F$ and Rosseland mean opacity $\kappa$, $A=3\kappa F/(16{\sigma }_{\mathrm{SB}})$.

There is an important limitation to treating $A$ as constant over the entire radiative layer. Integration from an irradiated outer boundary givesFor a diatomic ideal gas, ${\nabla }_{\mathrm{ad}}=2/7>1/4$, so this profile cannot actually reach a radiative-convective boundary. Formally, imposing a constant ${\nabla }_{\mathrm{ad}}$ would givewhich is positive only for $0<{\nabla }_{\mathrm{ad}}<1/4$. A finite boundary for a normal molecular atmosphere requires additional opacity, flux, or thermodynamic variation. The local matching formula is usable near a real boundary, but constant $A$ is not a complete global model of it. This is the convective stability of a constant-opacity irradiated atmosphere.

For the intended order-of-magnitude scaling, suppose the local values of $g/A$ and ${\nabla }_{\mathrm{ad}}$ are comparable for Jupiter and a hot Jupiter, and assume $T_b$ scales with planetary equilibrium temperature. Equal absorbed-flux factors around the same stellar luminosity give $T_b\propto a^{-1/2}$. Taking $a_J=5.2\mathrm{AU}$ and a representative close-in orbit $a_h=0.05\mathrm{AU}$ yieldsThis illustrates how irradiation can push a boundary much deeper. It is a conditional estimate calibrated from the supplied $0.1\mathrm{bar}$ reference, not a self-consistent prediction of the globally constant-$A$ model. Different intrinsic cooling flux, opacity, gravity, or atmospheric metallicity of a giant planet can substantially alter it.