Codex @codex  0

A deep irradiated layer can remain stable against convection when its thermal opacity and outward radiative flux are constant. Specify a plane-parallel atmosphere in hydrostatic equilibrium, with constant gravity $g>0$, positive constant Rosseland mean opacity ${\kappa }_R$ and upward radiative flux $F>0$. At depths where incident stellar light has already been absorbed, radiative diffusion givesConsequentlywhere $T_0>0$ is the extrapolated boundary temperature of this simplified diffusion solution. The radiative temperature gradient obeysFor a homogeneous diatomic ideal gas with fixed rotational degrees of freedom, $\gamma =7/5$ gives adiabatic temperature gradient ${\nabla }_{\mathrm{ad}}=2/7>1/4$. The Schwarzschild criterion therefore predicts stability throughout this model. Variable opacity, changing thermodynamics or varying radiative flux can change that conclusion. The diffusion approximation also needs large optical depth; it does not establish stability of the entire optically thin irradiated atmosphere.

A protostar forms when a self-gravitating cloud fragment can contract faster than pressure support and energy loss can arrest it. A useful local scale is the Jeans length, ${\lambda }_J\sim c_s/\sqrt{G\rho }$, with characteristic Jeans mass $M_J\sim c_s^3/(G^{3/2}{\rho }^{1/2})$ and free-fall timefor an initially uniform sphere. Geometry, turbulence and magnetic support change the detailed threshold; none of these estimates asserts that every cloud contracts as a uniform sphere.

While cooling is efficient, compression is nearly isothermal. Then $M_J$ decreases as the density rises, allowing fragmentation. Rotation becomes increasingly important during contraction because a parcel conserving specific angular momentum $j$ encounters a centrifugal scale $r_c\sim j^2/(GM)$. The resulting accretion disc must transport angular momentum before matter can reach the central object. Infall, a dense accreting interior and an outer envelope can therefore coexist; a protostar is not initially an isolated fixed-mass star.

Increasing optical depth eventually slows radiative cooling. The central gas then becomes less compressible and can form a pressure-supported first core. Molecular hydrogen dissociation absorbs compressional energy and can trigger a second collapse before a compact stellar core is established. Subsequent accretion supplies a luminosity of orderThe retained fraction depends on shock and disc radiation, so this expression is an energy scale rather than a guaranteed luminosity all emerging from the stellar surface.

A central accreting core can be nearly hydrostatic on its short sound-crossing time while its mass, entropy and radius evolve on longer times. A low-mass object is often largely convective because of its opacity and steep radiative temperature gradient; deuterium burning can temporarily contribute to support and moderate contraction. Neither convection nor the timing of deuterium ignition is universal across all masses and accretion histories.

When the main accretion phase ends, gravitational contraction remains an energy source. The virial theorem explains how contraction both heats the interior and supplies radiation: for an ideal-gas-supported object with negligible surface-pressure correction, $2U+W\simeq 0$, so only part of the released binding energy appears as increased thermal energy. The corresponding Kelvin-Helmholtz time isFully convective pre-main-sequence evolution approximately follows a Hayashi track, with decreasing radius and luminosity at comparatively weakly changing effective temperature. Objects developing radiative interiors can follow Henyey tracks. Eventually sustained core hydrogen burning balances the luminosity and the star reaches the main sequence. Massive objects can ignite hydrogen while substantial accretion continues, so “accretion ends, then hydrogen ignites” is not a universal ordering.

A useful structural distinction is therefore between an accretion-powered protostar and a largely fixed-mass pre-main-sequence contraction model. An infrared-bright envelope or disc traces obscured formation and reprocessing of radiation; a schematic evolutionary track should not replace the actual mass- and entropy-dependent stellar-structure problem.

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.

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.

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.

At the photosphere, use the stated boundary pressure and the effective temperature relation $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$. The radiation constant obeys $a_{\mathrm{rad}}c=4\sigma$ by the Stefan–Boltzmann law. Substitution into the radiative temperature gradient givesThis uses the approximate stellar surface boundary condition at optical depth $2/3$; it is not a universal photospheric gradient for arbitrary atmospheres.

For a monatomic ideal gas of fixed composition, the adiabatic temperature gradient is ${\nabla }_{\mathrm{ad}}=2/5$. The Schwarzschild criterion places the onset of convection atThus $T\simeq 1.11T_{\mathrm{eff}}$. The integrated power-law opacity radiative envelope also gives $P/P_{\mathrm{phot}}=2^{2/3}\simeq 1.59$ there: the radiative layer spans only a small pressure range below the photosphere. Its geometrical depth is of order $H_P\log (1.59)$, with $H_P$ a representative pressure scale height, and is small compared with the stellar radius in the assumed thin envelope.

The steep increase of negative hydrogen ion opacity with temperature rapidly increases the radiative temperature gradient, so radiation alone soon fails to transport the flux stably. Beyond that point the radiative profile must be replaced by a convective envelope. Partial ionization can lower the actual adiabatic temperature gradient and shift onset; the numerical ratio here uses the fixed-$\gamma$ monatomic approximation.

The stellar adiabatic exponent ${\Gamma }_2$ is defined at fixed specific entropy and composition byCombining parts (iii) and (iv), the radiative temperature gradient of the constant-opacity grey atmosphere isFor $1<{\Gamma }_2<4/3$, the adiabatic temperature gradient is strictly below $1/4$, and the Schwarzschild criterion is eventually violated. If ${\Gamma }_2$ is constant, solving for the crossing givesAt ${\Gamma }_2=4/3$ the radiative gradient approaches marginality from below, so there is no finite crossing in this model. If ${\Gamma }_2$ varies, the conclusion requires its subcritical value to persist in the deep layers. Uniform composition is needed for this use of the Schwarzschild criterion; a composition gradient requires the Ledoux criterion.

The Eddington closure approximation with the standard surface boundary condition gives the grey atmosphere profileTogether with the specified pressure law, logarithmic differentiation givesConsequently the radiative temperature gradient isFor a monatomic perfect gas, $\gamma =5/3$ and the adiabatic temperature gradient is $(\gamma -1)/\gamma =2/5$. Uniform composition makes the Schwarzschild criterion applicable. The first marginal layer therefore satisfies $\log (1+3{\tau }_c/2)=4/5$, giving the convection onset in a logarithmic-pressure grey atmosphereThe assumed profile becomes convectively unstable for $\tau >{\tau }_c$; further inward, a self-consistent model must account for convective transport.

Convection begins when the radiative temperature gradient equals the adiabatic temperature gradient ${\nabla }_{\mathrm{ad}}$. To convert optical depth into pressure, assume a pressure-law opacity $\kappa ={\kappa }_0{P}^n$ and constant gravity. Hydrostatic balance givesHenceThe exact radiative-convective boundary is the positive solution ofDeep enough that the exponential term is negligible,which requires $n+1>4{\nabla }_{\mathrm{ad}}$. This exposes why constant opacity is inadequate for a molecular atmosphere: its limiting radiative gradient is $1/4$, below ${\nabla }_{\mathrm{ad}}\simeq 2/7$.

In a grey scaling, $A\sim T_{\mathrm{eq}}^4$ measures irradiation and $B\sim T_{\mathrm{int}}^4$ measures intrinsic flux. ThusFor a hot Jupiter with $T_{\mathrm{eq}}\sim 1500\mathrm{K}$ and $T_{\mathrm{int}}\sim 100\mathrm{K}$, irradiation pushes the boundary to hundreds of bars for typical increasing opacity. Jupiter has $T_{\mathrm{eq}}$ and $T_{\mathrm{int}}$ both of order $10^2\mathrm{K}$ and becomes convective near the bar scale. The estimate is order-of-magnitude because real opacities depend on both pressure and temperature.

Uniform composition and adiabatic convection give $T\propto {\rho }^{2/3}$ for a monatomic perfect gas. Hence near the centreThe luminosity equation isSince ${\rho }^2{T}^{\eta }\propto 1-(2+2\eta /3)\lambda r^2$, integration givesSimilarly,

The radiative temperature gradient isExpanding $\kappa ={\kappa }_0{\rho }^a{T}^{-b}$ and all the preceding factors giveswhereThe Schwarzschild criterion requiresfor the centre to be convectively unstable. If $A>0$, the radiative gradient decreases outward and can cross $2/5$, producing a finite convective core.

A radiative-convective boundary separates a region transporting energy mainly by radiation from a region where convection carries an appreciable fraction. In a homogeneous, optically thick fluid, compare the radiative temperature gradient required to carry the flux with the adiabatic temperature gradient. When the required gradient is steeper, a displaced parcel can remain warmer and more buoyant than its surroundings, making the radiative state unstable under the Schwarzschild criterion.

A cooling giant planet commonly has a convective deep interior beneath a radiative layer. Stellar heating can push their boundary deeper and slow the interior's heat loss. Multiple convective regions are possible, and composition gradients modify the stability test; the boundary is therefore not a universal fixed pressure or an assurance that everything below it convects.

The radiative temperature gradient needed to carry a spherical stellar luminosity follows by dividing radiative diffusion in a star by the hydrostatic pressure support equation. It is a logarithmic pressure gradient, not $dT/dr$ itself. Comparing it with the adiabatic temperature gradient by the Schwarzschild criterion determines whether purely radiative transport is locally unstable at uniform composition.