Inhibiting Conduction by He Mixing in Interiors of Jupiter and Saturn

Accurate knowledge of the electrical and thermal conductivities and structural properties of hydrogen-helium mixtures under thermodynamic conditions within and beyond the immiscibility range is very important to predict the thermal evolution and internal structure of gas giant planets like Jupiter and Saturn. Here, we propose a novel method to determine the immiscibility boundary accurately without the need for free energy calculations, while providing consistent insights into structural and transport properties of mixtures. We show with direct large-scale ab initio simulations that the insulator-metal transition (IMT) of the hydrogen subsystem is strongly affected by an admixture with a small fraction of helium and occurs at temperatures significantly higher than those of pure hydrogen. At pressures below 150 GPa, the IMT boundary is not related anymore to the H2 subsystem dissociation, the system remains insulating even after the full dissociation of H2 molecules and its transition to an H-He mixture. The offset of the IMT in the H-He mixture relative to the dissociation region in the hydrogen subsystem and the significant reduction of static electrical and thermal conductivity by a factor between two and a few thousand relative to pure hydrogen found in mixtures have consequences for Jupiter and Saturn’s thermal evolution, internal structure, and dynamo action, affecting a large fraction of the interior of both planets.

💡 Research Summary

This paper addresses a critical gap in our understanding of the electrical, thermal, and structural properties of hydrogen‑helium (H‑He) mixtures under the extreme pressure–temperature conditions that prevail inside gas‑giant planets such as Jupiter and Saturn. While previous work has largely relied on free‑energy differences (ΔG) to map the H‑He immiscibility region, those methods provide little insight into the microscopic structure or transport behavior of the mixtures and often suffer from finite‑size limitations.

The authors perform large‑scale ab‑initio molecular dynamics (AIMD) simulations in the NPT ensemble using 1 024 electrons for two helium concentrations (x = 0.11304 and x = 0.27522). They employ the temperature‑dependent KDT16 generalized‑gradient approximation (GGA) functional to capture thermal exchange‑correlation effects, and for selected state points they cross‑check band gaps and conductivities with the hybrid HSE06 functional. Electrical conductivity is obtained via the Kubo‑Greenwood formalism, and the insulator‑to‑metal transition (IMT) is identified using Mott’s criterion (σ ≈ 2000 S cm⁻¹).

A key methodological innovation is the use of the first peak of the H‑He radial distribution function (RDF) as a direct, “mechanical” indicator of phase separation. In a demixed system, helium‑rich and helium‑poor regions form a surface‑like interface, dramatically reducing the probability of finding a helium atom near a hydrogen atom; this suppression appears as a lowered first RDF peak. As temperature rises, the peak first increases (partial mixing) then decreases due to thermal expansion. The temperature at which the peak reaches its maximum corresponds to the transition from a demixed to a perfectly mixed state. This approach eliminates the need for costly free‑energy integrations while delivering immiscibility boundaries with an estimated uncertainty of ±500 K.

The simulations reveal that even a modest helium admixture (≈10 % by number) shifts the IMT of the hydrogen subsystem to significantly higher temperatures—by several hundred to a few thousand kelvin—compared with pure hydrogen. At 150 GPa, the H₂‑He mixture metallizes concurrently with H₂ dissociation (≈1500 K), but at 75 GPa the metallization temperature (≈4000 K) far exceeds the H₂ dissociation temperature (≈2750 K). Consequently, a wide region of atomic hydrogen‑helium remains insulating. The static electrical conductivity of the mixtures is reduced by factors ranging from two to several thousand relative to pure hydrogen at the same pressure and temperature. Helium also delays H₂ dissociation by ≈250 K, reflecting a strengthening of the H₂ bond in the presence of helium.

These findings have profound implications for planetary modeling. The reduced conductivity implies a lower thermal conductivity in the deep interiors, slowing the outward transport of heat and extending the cooling timescales of Jupiter and Saturn. The insulating character of H‑He mixtures in a substantial fraction of the planetary interior also modifies the conditions required for dynamo action, potentially explaining observed differences in magnetic field morphology between the two planets. Moreover, the immiscibility boundary obtained with KDT16 aligns closely with recent evolutionary models that incorporate helium rain, supporting the hypothesis that helium phase separation contributes significantly to Saturn’s excess luminosity.

In summary, the paper delivers three major contributions: (1) a novel, free‑energy‑independent method to locate the H‑He immiscibility line using RDF peak analysis; (2) a detailed quantification of how helium admixture raises the IMT temperature and suppresses both electrical and thermal conductivities; and (3) a comprehensive discussion of how these microscopic material changes affect macroscopic planetary properties such as thermal evolution, interior structure, and magnetic field generation. The work sets a new benchmark for high‑pressure planetary material simulations and provides a robust framework for future experimental and theoretical investigations of H‑He mixtures.