Interior dynamics of envelopes around disk-embedded planets

Using three-dimensional hydrodynamical simulations, this study reveals that the thermal balance between heating and cooling in disk-embedded planets creates three distinct envelope regimes, where intermediate cooling rates form a radiative layer that traps volatiles while slow cooling in inner disks leads to fully convective envelopes that stall super-Earth growth and deplete volatiles.

The Gist

Imagine a young planet, like a toddler in a cosmic nursery, trying to grow up. It starts as a rocky core, but to become a giant like Jupiter or even a large "Super-Earth," it needs to grab a massive blanket of gas from the swirling disk of dust and gas surrounding its star. This gas blanket is called an atmospheric envelope.

For decades, scientists thought this process was like a quiet, still pond: the planet sits there, the gas slowly settles down, and the planet gets heavier. But this new paper by Kuwahara and Lambrechts shows that the reality is much more like a bustling, windy city square. The gas isn't just sitting there; it's swirling, flowing in, flowing out, and constantly exchanging with the surroundings.

Here is the story of their discovery, broken down into simple concepts:

1. The Battle: Heating vs. Cooling

Think of the planet's atmosphere as a room with a heater and an air conditioner.

• The Heater (Accretion): As the planet eats up rocks and ice (pebbles), the impact releases energy, heating up the gas. This makes the gas want to expand and puff up.
• The Air Conditioner (Cooling): The gas tries to radiate that heat away into space to cool down and settle.

The paper asks: What happens when we turn the heater up or down, and when the air conditioner works fast or slow?

2. The Three "Moods" of the Atmosphere

The authors ran super-computer simulations to see how the atmosphere behaves under different conditions. They found the atmosphere has three distinct "personalities" or regimes, depending on how fast it can cool down:

A. The "Fast Cools" (The Calm, Thin Blanket)

• When: The air conditioner is super efficient (cooling time is very short).
• What happens: The gas stays cool and thin. It's like a calm, still lake.
• The Result: The gas near the planet is "shielded." It doesn't mix much with the outside world. The planet keeps its own little pocket of gas, and the outer gas just flows around it like a river around a rock.

B. The "Slow Cools" (The Chaotic, Boiling Pot)

• When: The air conditioner is broken or very weak (cooling time is very long). This usually happens close to the star where it's hot.
• What happens: The heat builds up, making the gas boil and churn. It becomes fully convective.
• The Result: Imagine a pot of boiling water. Everything mixes! Gas from the deep interior gets thrown all the way to the surface, and outside gas gets sucked all the way down. The planet's atmosphere is a giant, churning mixer.

C. The "Just Right" (The Three-Layer Cake)

• When: The cooling is in the middle range.
• What happens: This is the most complex and interesting state. The atmosphere forms a three-layer cake:
  1. Inner Core: A boiling, convective layer (like the bottom of the pot).
  2. Middle Layer: A calm, radiative layer (like a lid). This layer is stable and doesn't mix much.
  3. Outer Layer: A recycling layer where gas flows in and out from the surrounding disk.
• The Result: The middle "lid" acts as a shield. It traps things inside.

3. The Great Trapper: Why This Matters for Planets

This is the most exciting part of the paper. The "lid" (the middle radiative layer) changes what the planet is made of.

• The Volatile Trap: When the planet eats icy pebbles, the ice turns into gas (vapor) inside the hot atmosphere.
  • In the "Boiling Pot" (Slow Cooling): The churning gas throws this new vapor right back out into space. The planet loses its water and other volatiles. It ends up dry and rocky.
  • In the "Three-Layer Cake" (Intermediate Cooling): The middle "lid" traps the vapor. The gas can't escape easily, so the water and other ingredients get stuck inside the planet's deep interior. The planet becomes rich in water and volatiles.

4. The Big Picture: Where You Live Matters

The authors realized that where a planet forms in the solar system determines its personality:

• Inner Solar System (Close to the Star): It's hot, so cooling is slow. Planets here tend to become "Boiling Pots." They mix everything up, lose their water, and might stop growing early because the churning gas pushes away new pebbles. This might explain why our inner planets (like Earth and Mars) are rocky and relatively dry compared to giants.
• Outer Solar System (Far from the Star): It's cold, so cooling is fast. Planets here form "Three-Layer Cakes." They trap their water and volatiles deep inside. This helps them grow massive and become gas giants or water-rich worlds.

Summary Analogy

Imagine two people trying to build a snowman in a blizzard:

• Person A (Inner Disk): The wind is so hot and chaotic that every time they try to pack snow, it melts and blows away. They can't build a big, wet snowman; they end up with a small, dry rock.
• Person B (Outer Disk): The wind is cold and steady. They can build a protective wall (the radiative layer) that keeps the snow from melting and blowing away. They can pack the snow deep inside, creating a massive, water-rich snowman.

The Takeaway:
This paper teaches us that a planet isn't just a static ball of rock and gas. It's a dynamic, breathing system. How fast it cools down decides whether it stays dry and small or becomes a wet, massive giant. It's a crucial piece of the puzzle in understanding why our solar system looks the way it does.

Technical Summary

Here is a detailed technical summary of the paper "Interior dynamics of envelopes around disk-embedded planets" by Kuwahara and Lambrechts (2026).

1. Problem Statement

In the core accretion scenario of planet formation, protoplanets accrete solid cores and subsequently capture gaseous envelopes. Conventional one-dimensional (1D) models treat these envelopes as static, isolated, and spherically symmetric structures in hydrostatic equilibrium. However, recent three-dimensional (3D) hydrodynamical simulations have revealed that envelopes are dynamic systems interacting with the surrounding protoplanetary disk through recycling flows (gas exchange between the envelope and the disk).

The efficiency of this recycling process, and consequently the envelope's ability to cool and contract (leading to runaway gas accretion), depends on the balance between:

• Heating: Primarily from the gravitational potential energy released by accreting solids (pebbles).
• Cooling: Regulated by radiative cooling, which is highly uncertain due to poorly constrained opacities (dust content) within the envelope.

Previous studies have not systematically explored the wide range of cooling and heating rates simultaneously. This paper aims to fill that gap by investigating how the interplay between accretion heating and radiative cooling (parameterized by the cooling timescale) dictates the internal structure, dynamics, and material transport of planetary envelopes.

2. Methodology

The authors performed a comprehensive suite of 3D hydrodynamical simulations using the Athena++ code.

• Physical Setup:
  • Geometry: Spherical polar coordinates centered on a planet embedded in a non-self-gravitating disk.
  • Fluid Dynamics: Solved the continuity, Euler, and energy equations for a compressible, inviscid fluid.
  • Forces: Included Coriolis and tidal forces, as well as the planet's gravitational potential (smoothed to avoid numerical heating at the inner boundary).
  • Heating ($Q_{acc}$): Modeled as a volumetric heating source derived from pebble accretion luminosity ($L_{acc}$), distributed radially to avoid artificial temperature jumps at the surface.
  • Cooling ($Q_{cool}$): Implemented as a thermal relaxation process (Newtonian cooling) with a characteristic timescale $t_{cool}$.
• Key Parameters:
  • Cooling Regime: Controlled by the dimensionless cooling time $\beta \equiv t_{cool}\Omega$ (where $\Omega$ is the orbital frequency). The study explored $\beta$ ranging from $0.01$to$10,000$.
  • Planetary Mass: Varied around the thermal mass ($m \approx 0.1 - 0.3$).
  • Accretion Rate: Varied by two orders of magnitude ($1$to$100 M_\oplus$/Myr) to test luminosity dependence.
  • Tracer Particles: Passive scalars were used to track the transport of gas, small dust, and volatiles within the envelope.
• Numerical Resolution: High radial resolution ($N_r = 256$) to resolve the Bondi radius ($R_B$) with $\sim118$ grids, ensuring accurate capture of hydrostatic balance and flow structures.

3. Key Contributions and Results

A. Identification of Three Distinct Cooling Regimes

The simulations revealed that the internal structure of the envelope is determined by the cooling timescale $\beta$, leading to three distinct regimes:

1. Fast Cooling Regime ($\beta \lesssim 1$):

   • Structure: The envelope is nearly isothermal.
   • Dynamics: An inner radiative layer forms, shielded from the outer disk by a buoyancy barrier (positive entropy gradient).
   • Recycling: The outer envelope undergoes recycling flows (horseshoe orbits entering from poles, exiting at the midplane), but the inner radiative layer is effectively isolated from these flows.
   • Flow Speed: Low, dominated by shear at the outer edge.

2. Intermediate Cooling Regime ($1 \lesssim \beta \lesssim 300$):

   • Structure: A three-layer structure emerges:
     1. Inner convective layer.
     2. Middle radiative layer.
     3. Outer recycling layer.
   • Dynamics: The radiative layer acts as a barrier. While the inner convective zone mixes material, the radiative layer traps small dust and vapors released from sublimating solids, preventing them from being advected out to the disk.
   • Flow Speed: Convective velocities are a few percent of the sound speed; flow speeds drop significantly in the radiative layer.

3. Slow Cooling Regime ($\beta \gtrsim 1000$):

   • Structure: The envelope becomes fully convective.
   • Dynamics: No stable radiative barrier exists. Material enters and leaves the envelope isotropically.
   • Recycling: Highly efficient mixing throughout the entire envelope.
   • Flow Speed: Convective velocities reach $\sim10\%$ of the sound speed.

B. Analytic Models for Flow and Transport

The authors derived analytic formulae for:

• Convective Flow Speed: Based on mixing-length theory, $v_{con} \propto (L_{acc}/\rho)^{1/3}$. This successfully predicted the flow speeds observed in simulations.
• Material Transport Timescales:
  • Fully Convective: Tracers are exchanged with the disk on a timescale of $\tau \sim 3$ orbits.
  • Radiative/Recycling: Tracers in the radiative layer are trapped for extremely long timescales ($\tau > 10^4$ orbits), effectively decoupling the deep interior from the disk.

C. Implications for Volatile Delivery and Planetary Growth

By comparing the volatile enrichment timescale ($\tau_{vol}$) with the material exchange timescale ($\tau_{rec}$), the authors identified a spatial dichotomy in planet formation:

• Inner Disk ($\lesssim 1$ au): Cooling times are long (due to high temperatures/density), leading to fully convective envelopes. Here, $\tau_{rec} < \tau_{vol}$. Sublimated volatiles (e.g., water) are rapidly recycled back into the disk before they can be retained. This implies inner-disk super-Earths may be volatile-depleted and have stalled growth.
• Outer Disk ($\gtrsim 10$ au): Cooling is efficient, allowing radiative layers to form. Here, $\tau_{rec} \gg \tau_{vol}$. Sublimated volatiles are trapped in the deep interior, leading to volatile-rich planets.

4. Significance

This work fundamentally challenges the static 1D view of planetary envelopes by demonstrating that 3D dynamics and thermodynamics are inextricably linked.

• Resolution of the Opacity Problem: By parameterizing cooling via $\beta$, the study provides a framework to understand how uncertain dust opacities translate into distinct structural regimes (convective vs. radiative).
• Planet Composition: It offers a physical explanation for the potential compositional differences between inner (volatile-poor) and outer (volatile-rich) planets, suggesting that envelope dynamics, not just disk chemistry, dictate final planetary composition.
• Growth Stalling: The findings suggest that in the inner disk, the formation of fully convective envelopes may inhibit the accumulation of heavy elements and volatiles, potentially stalling the growth of super-Earths before they can undergo runaway gas accretion.
• Future Modeling: The paper highlights the necessity of coupling hydrodynamics with dust evolution (growth, fragmentation, settling) and radiative transfer to accurately model the co-evolution of gas, solids, and volatiles during planet formation.

In conclusion, Kuwahara and Lambrechts establish that the thermal state of a planetary envelope is not merely a passive background condition but an active regulator of mass and volatile accretion, with profound implications for the diversity of exoplanets observed in our galaxy.