🔭 astrophysics

Self-Consistent Evolution Models Show Weak Double-Diffusive Mixing in Jupiter and Saturn

This study demonstrates that self-consistent evolutionary models of Jupiter and Saturn reveal double-diffusive convection is too inefficient to significantly redistribute heavy elements or erase primordial compositional gradients, implying that additional mechanisms or formation pathways are required to explain the observed heavy-element distributions.

Original authors: J. R. Fuentes, Ankan Sur, David J. Stevenson, Peter Bodenheimer

Published 2026-07-07

📖 3 min read☕ Coffee break read

CC BY 4.0

Original authors: J. R. Fuentes, Ankan Sur, David J. Stevenson, Peter Bodenheimer

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine Jupiter and Saturn as giant, swirling balls of gas and rock. For a long time, scientists have been puzzled by their insides. We know they have "cores" made of heavy stuff (like rock and ice), but recent data suggests these cores aren't sharp, distinct balls at the center. Instead, they look like "fuzzy" clouds of heavy material that blend gradually into the surrounding gas.

The big question was: How did these fuzzy cores get that way?

One popular theory was that after the planets formed, a process called "double-diffusive convection" acted like a slow-motion blender, mixing the heavy core material up into the gas layers over billions of years. Think of it like stirring honey into coffee; eventually, you'd expect the honey to spread out evenly.

The New Discovery
This paper, written by a team of researchers, ran a massive, 4.5-billion-year simulation of Jupiter and Saturn to test if this "blender" actually works. They used a sophisticated computer code (named APPLE) to track how heat and heavy materials move inside these planets, accounting for the complex physics of deep space.

The Result: The Blender is Broken
The simulation showed that this mixing process is incredibly weak. It's like trying to stir a giant vat of thick molasses with a feather.

The Limit: Over the entire history of the Solar System, the "blender" managed to move less than one Earth-worth of heavy material from the core into the outer layers.
The Reason: The paper explains that in the deep, high-pressure interiors of these planets, the energy required to mix the heavy stuff is huge. The heat energy available to do the "stirring" is mostly used just to keep the planet warm and glowing, leaving very little leftover energy to actually move the heavy rocks and ice around.
The Analogy: Imagine trying to push a heavy boulder up a hill. In a small lab experiment (which previous studies used), the hill is tiny, and you have a strong engine, so the boulder moves easily. But in a real giant planet, the hill is the size of a mountain, and your engine is barely running. You simply don't have enough fuel to push the boulder very far.

What This Means
Because this "blender" is so inefficient, it cannot be the reason why Jupiter and Saturn have fuzzy cores. The heavy material didn't get mixed up after the planets formed.

Instead, the paper suggests that the "fuzziness" must have happened during the planet's birth.

The Birth Theory: When Jupiter and Saturn were forming, they likely swallowed up chunks of rock and ice (planetesimals) that melted and spread out before the planet fully settled down.
Other Possibilities: Maybe giant collisions during their formation smashed the core apart, or perhaps our models of how planets eat up gas and rock are missing some early mixing steps.

In Summary
The paper concludes that double-diffusive convection is too weak to explain the current structure of Jupiter and Saturn. The "fuzzy cores" we see today are likely a fossil of how these planets were born, not the result of a slow, billion-year mixing process. To understand them better, we need to look closer at the violent and chaotic events of their formation, rather than just their quiet cooling afterward.

Technical Summary: Self-Consistent Evolution Models Show Weak Double-Diffusive Mixing in Jupiter and Saturn

Problem Statement
The internal structure of giant planets, particularly the existence of "fuzzy" or extended cores in Jupiter and Saturn, remains a significant challenge for planetary science. While double-diffusive convection has been proposed as a mechanism to redistribute heavy elements from a dense core into the overlying envelope, its efficiency in realistic evolutionary contexts is uncertain. Previous studies have relied on idealized compositional structures and have not treated double-diffusive transport self-consistently within planetary evolution calculations. Furthermore, there is a discrepancy between local hydrodynamical simulations, which often predict rapid layer merging and efficient mixing, and the requirements of global evolutionary models that must reproduce present-day observational constraints.

Methodology
The authors implement a self-consistent treatment of double-diffusive transport within the planetary evolution code APPLE.

Initial Conditions: The study utilizes post-formation interior models for Jupiter (based on Stevenson et al. 2022) and Saturn (based on Bodenheimer et al. 2025). These models feature extended compositional gradients ("fuzzy cores") resulting from the partial disruption and evaporation of accreted planetesimals during formation, rather than sharp core boundaries.
Transport Physics: The authors incorporate the layered-convection prescription derived by Spruit (2013). This approach models transport across convective staircases as an enhanced diffusive process. The effective compositional ( $D_{\text{eff}}$ ) and thermal ( $\kappa_{T,\text{eff}}$ ) diffusivities are derived based on asymptotic regimes where transport is independent of layer thickness, scaling with the stability ratio ( $R_\rho$ ) and microscopic diffusivities.
Evolutionary Simulation: The models are evolved for 4.56 Gyr (the age of the Solar System). The code simultaneously solves for thermal evolution and compositional transport using a Lagrangian mass grid. Regions unstable to overturning convection (Ledoux criterion) are treated with standard mixing-length theory, while Ledoux-stable regions with $1 < R_\rho < \tau^{-1/2}$ utilize the double-diffusive prescription.
Sensitivity Analysis: To test the robustness of the results, the authors run simulations where the effective diffusivities are enhanced by factors of 100 and 1000 relative to the fiducial Spruit (2013) values.

Key Results

Inefficient Mixing: The simulations reveal that double-diffusive convection produces very limited mixing between the deep interior and the envelope. Over the full 4.56 Gyr cooling history, less than $\sim 1 M_\oplus$ $\sim 1 M_{\oplus}$ of heavy material is redistributed in both Jupiter and Saturn.
- In the fiducial model, only $\approx 0.11 M_\oplus$ (Jupiter) and $\approx 0.07 M_\oplus$ (Saturn) are redistributed.
- Even with diffusivities enhanced by a factor of 1000, the redistributed mass remains modest ( $\approx 2 M_\oplus$ for Jupiter and $\approx 1.3 M_\oplus$ for Saturn).
Self-Limiting Mechanism: The inefficiency arises from a negative feedback loop. Enhanced thermal transport modifies the thermal stratification, reducing the stability ratio ( $R_\rho$ ) and thereby suppressing the effective diffusivities. This prevents runaway mixing, ensuring that even unrealistically large transport coefficients fail to homogenize the primordial compositional gradients.
Energy Constraint: The authors identify a fundamental energetic constraint. In the Boussinesq approximation used in many hydrodynamical simulations, mixing is energetically cheap because the layer thickness is small compared to the thermal scale height. However, in compressible giant planet interiors, the stable region spans a substantial fraction of the planetary radius ( $H \sim 0.25\text{--}0.5 R_p$ ), making the energy required to mix the compositional gradient comparable to the available thermal energy budget. Consequently, the buoyancy work available to drive compositional transport is severely constrained.

Significance and Claims
The paper claims that double-diffusive convection, when treated self-consistently in a global evolutionary context, cannot significantly erode the compositional gradients generated during planet formation. The observed heavy-element distributions in Jupiter and Saturn likely require additional transport mechanisms or formation pathways that produce broader initial mixing than standard accretion models predict.

The authors suggest that the "fuzzy" cores may be a direct product of the formation process itself (e.g., more extensive early mixing during planetesimal ablation) or the result of large collisional events, rather than a product of post-formation double-diffusive mixing. The study highlights a critical limitation in applying local hydrodynamical results (often based on the Boussinesq approximation) to the global evolution of giant planets, emphasizing that the energy budget of the deep interior is the primary regulator of mixing efficiency.

Technical Summary: Self-Consistent Evolution Models Show Weak Double-Diffusive Mixing in Jupiter and Saturn

More like this