Let increase outward and define the inward optical depth by . With constant outward internal radiative flux , radiative diffusion gives
The 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, , so
The 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.
Figure 1.
Intrinsic grey-atmosphere cooling profile compared with illustrative irradiated hot-Jupiter profiles
. Details.
The intrinsic curve follows the grey formula. The irradiated curves illustrate possible shapes only; they are not solutions for a specified opacity model. Smaller optical depth corresponds to greater altitude.
Assume a circular orbit of radius , a spherical planet, bolometric Bond albedo , negligible internal heat, unit thermal emissivity and complete redistribution of absorbed heat over the sphere. The stellar luminosity is , and the planet absorbs the incident radiative flux through its projected cross-section:
Thermal reradiation is . Equating these gives the planetary equilibrium temperature
The planet's radius cancels. A non-black thermal emissivity divides the absorbed flux by in the fourth-power balance. For uniform dayside-only reradiation, the emitting area is and . With no local redistribution, the substellar point has , while other points depend on incidence angle. These are different temperature conventions, not contradictory formulas.
The planetary equilibrium temperature is not the surface greenhouse temperature or the internal effective temperature of a planet. If intrinsic cooling matters and both powers escape through the same emitting area, the total effective temperature obeys .
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 , and the upper boundary supplies the usual term. Treating the incoming radiation with representative direction cosine gives its attenuation . Intrinsic flux enters from below, while the imposed external stellar flux is absorbed from above.
With constant gravity and hydrostatic equilibrium, , 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.
For a plane-parallel atmosphere with optical depth increasing downward, the radiative transfer equation is
Integrating over solid angle gives the first radiation-field moment
In radiative equilibrium, matter has no net local radiative heating, so the opacity-weighted frequency integral of vanishes. After converting each optical-depth derivative to physical depth and integrating over frequency,
Thus the bolometric internal flux is constant with depth, as stated by constant flux in a plane-parallel radiative-equilibrium atmosphere.
In local thermal equilibrium . Define the planet's internal effective temperature of a planet by the constant outward flux. The standard Planck law integral gives
The corresponding intrinsic luminosity is .
The internal effective temperature of a planet parametrizes intrinsic cooling, while the irradiation temperature parametrizes incident stellar flux before the redistribution factor . The quantity
is the visible-to-thermal mean-opacity ratio, and is downward thermal optical depth. Under the Eddington two-stream boundary condition,
while the common semi-grey choice is .
For a young giant at , stellar heating is weak and a still-large dominates. Its pressure-temperature profile rises steadily inward, approximately as , and joins a deep convective adiabat.
For a hot Jupiter close to a Sun-like star, . It has a broad, nearly isothermal irradiated radiative layer, followed by a deep rise where intrinsic flux and increasing opacity matter. If , absorption of starlight above the thermal photosphere can create an atmospheric thermal inversion.
For a temperate sub-Neptune around an M dwarf, irradiation and internal cooling can be more comparable. Its profile generally has a moderate radiative layer above a convective interior; near-infrared stellar radiation, molecular opacity, clouds, and hazes determine whether the upper profile is weakly inverted or decreases outward.