In a groundbreaking study poised to reshape our understanding of terrestrial exoplanet evolution, researchers have unveiled a comprehensive simulation framework that illuminates the volatile-rich history of the molten super-Earth L 98-59 d. This innovative effort combines interior dynamics with atmospheric processes to reveal how such planets evolve over billions of years while retaining vast inventories of critical volatile elements. The study, conducted by Nicholls, Lichtenberg, Chatterjee, and colleagues, leverages a modular simulation architecture named PROTEUS, which synergistically integrates complex geophysical and atmospheric components, thus offering an unprecedented window into planetary evolution beyond our solar system.
At the heart of the PROTEUS framework lies a seamless coupling between SPIDER, a sophisticated model of planetary interior evolution, and AGNI, a cutting-edge atmosphere model grounded in radiative–convective equilibrium physics. The synthesis of these models permits a fully self-consistent portrayal of how mantle dynamics, phase transitions, and volatile partitioning interact with atmospheric composition and escape mechanisms across geologic timescales. Crucially, the approach incorporates state-of-the-art thermochemical equilibrium calculations accounting for volatile solubility in magma oceans, which dominate the interiors of molten super-Earths like L 98-59 d.
The research team assumes that the primordial protoplanetary nebula has already dissipated, focusing instead on the secondary atmosphere’s evolution driven by outgassing, volatile partitioning, and ongoing atmospheric escape. The volatile budget, including elements essential for habitability such as carbon, hydrogen, nitrogen, and sulfur (CHNS), is initialized based on plausible planetary parameters consistent with observational constraints. By simulating these complex interactions, the study circumvents the need to model volcanic outgassing from fully solidified mantles, a process orders of magnitude less efficient than degassing from magma oceans, which remain prevalent in the scenarios examined.
SPIDER is tasked with tracking the planet’s molten mantle, composed primarily of MgSiO3 melt, along an initial adiabatic temperature profile before evolving it in response to various heating mechanisms. These include radiogenic decay, core cooling, and particularly tidal heating, which they model with the Maxwell viscoelastic rheology using LOVEPY software. This rheological treatment accounts for mantle viscosity variations dependent on melt fraction and temperature, yielding realistic heating profiles that can maintain large-scale mantle melt fractions over billions of years. Interestingly, the mantle is considered chemically inert in the current simulations, thus simplifying the thermodynamic calculations while focusing on physical processes.
The core-mantle structure draws analogies with Earth but incorporates adjustments for exoplanet-specific conditions such as mantle pressure and potential incorporation of lighter elements into the core, which modulate core density relative to pure iron. The researchers highlight that, for L 98-59 d, the core mass is approximately 0.66 Earth masses, with the mantle mass closer to 1.48 Earth masses. By holding the planetary interior radius constant over time, they effectively disentangle changes in mantle phase from radius evolution, streamlining the model without sacrificing realism in capturing thermal evolution dynamics.
Atmosphere modeling is managed by AGNI which simulates radiative and convective heat transport through a thick, outgassed atmosphere. The team deploys SOCRATES, a state-of-the-art correlated-k radiative transfer code with 48 spectral bands and angle-dependent treatment, to accurately calculate the net atmospheric energy flux even under dense, volatile-rich atmospheric conditions. Rayleigh scattering and basaltic surface emissivity are integrated, adding correction factors that realistically simulate planetary albedo and infrared emission. The code employs mixing-length theory for convective transport, ensuring that atmospheric heat fluxes balance with the mantle’s thermal output, closing the energy budget at each time step.
Stellar evolution is intimately tied to planetary atmospheric dynamics, and the researchers accordingly incorporate L 98-59’s properties through the MORS stellar evolution model. MORS realistically tracks the star’s luminosity, radius, and extreme ultraviolet (XUV) flux evolution, critical parameters influencing atmospheric escape and surface temperature. With a stellar mass of 0.273 solar masses and a rotation period of approximately 81 days, L 98-59 evolves in a way that gradually reduces its high-energy emissions. This temporal modulation directly affects the rate of hydrodynamic escape from the planet’s atmosphere, providing a dynamic context for volatile retention or loss over billions of years.
The simulations commence at an assumed planetary age of 50 million years, a juncture where the initial rapid atmospheric boil-off has subsided and volatile inventories stabilize. This choice reflects current understanding that early escape processes strongly deplete primordial hydrogen envelopes, after which secondary atmospheres evolve via outgassing and escape. The simulations continue until either mantle solidification occurs or the system reaches the present estimated age of about 4.94 billion years. This timescale ensures that the model encompasses the crucial phases during which volatile partitioning and escape shape the planet’s final atmospheric composition and bulk density.
A pivotal component of the volatile evolution modeled involves atmospheric escape driven by stellar XUV flux. Employing the classic energy-limited photoevaporation formula, the framework calculates the mass-loss rate based on the planet’s gravitational potential and the effective radius at which XUV photons are absorbed. The escape efficiency is fixed near 10%, consistent with high-fidelity hydrodynamic simulations that account for cooling effects in the upper atmosphere. Notably, the atmospheric composition is assumed to escape en masse without fractionation under intense stellar irradiation, reflecting rapid hydrodynamic flow regimes anticipated for planets receiving high XUV fluxes.
Importantly, although the escape process here is non-fractionating, the elemental composition of the escaping gas reflects that of the atmosphere rather than the mantle or core. This distinction leads to ongoing volatile fractionation within the planet over time, as partitioning between interior and atmosphere varies by element. Such a framework naturally captures how high-molecular-weight elements like sulfur and nitrogen may become enriched or depleted relative to hydrogen, informing interpretations of observational data on atmospheric compositions derived from transit and spectroscopic measurements.
The model’s fidelity extends to accurately estimating the observable planetary radius, often measured near the 20 mbar pressure level, which is where transit measurements effectively probe the atmosphere. By coupling radiative-convective atmospheric profiles with hydrodynamic escape simulations, the approach aligns theoretical predictions closely with observable parameters. The planetary bulk density, an essential metric for assessing volatile content and internal structure, is computed using precise mass-radius relationships driven by the model outputs, enabling direct comparisons with measured densities from transit timing variations and radial velocity data.
One of the most compelling aspects of the study is the comparative analysis between the simulated data and recently revised bulk density estimates of L 98-59 d. Earlier estimates placed the density near 3.45 g/cm³, while newer analyses suggest even lower values around 2.2 g/cm³. The authors caution that their volatile inventory estimates are conservative under the older density assumption, and that the lower density findings only reinforce the conclusion that L 98-59 d harbors a substantial volatile component. This finding challenges traditional models of rocky super-Earths as dry, refractory bodies, instead painting a picture of planets that retain deep magma oceans enriched with volatiles over geological timescales.
Hearteningly, the sensitivity analyses embedded within the study confirm robustness against variations in core size assumptions, escape efficiency, and stellar irradiation conditions. This rigorous testing ensures that the volatile-rich evolutionary pathway is not an artifact of narrow parameter choices but rather a natural outcome given plausible physical inputs. The results broadly support emerging paradigms wherein small, close-in exoplanets sustain magma oceans and thick atmospheres, shaped by their complex interplay with host star environments.
The implications of these findings for exoplanet habitability and characterization are profound. Volatile retention, particularly of water and carbon-bearing species, directly impacts surface conditions and the potential for life-supporting environments. By illustrating how intense early irradiation does not inevitably strip volatile reservoirs completely, the study opens new avenues for probing whether super-Earths orbiting M-dwarf stars could host detectable atmospheres amenable to biosignature searches with forthcoming telescopes.
By meticulously combining interior thermodynamics, atmospheric physics, and stellar evolution into a unified model, this research sets a new benchmark for interrogating the histories of terrestrial exoplanets. It showcases how interdisciplinary approaches, grounded in planetary science and astrophysics, can decrypt the subtle processes shaping alien worlds light-years away. As exoplanet discovery marches forward, tools like PROTEUS and AGNI promise to guide interpretations of ever more detailed observational data, ushering in an era of nuanced understanding of planetary architecture and habitability in the cosmos.
The study’s integration of high-fidelity physical modeling with observational constraints marks a decisive advance in exoplanet science—a realm where bold theories and detailed data increasingly intersect. By demonstrating that L 98-59 d, a world initially thought potentially barren, likely evolved with a thick, volatile-rich atmosphere sustained by a persistent magma ocean, it compels reconsideration of how common such evolutionary trajectories may be throughout the galaxy. This insight challenges prevailing assumptions and invites further exploration into the diversity of rocky exoplanet pathways.
Ultimately, this work exemplifies how modern planetary science deploys computational innovation to unravel the inner lives of distant worlds. With telescopes poised to examine exoplanet atmospheres in exquisite detail, understanding the coupled evolution of interiors and atmospheres will prove indispensable. The volatile-rich evolution scenario posited here therefore not only reshapes our knowledge of a single super-Earth but also propels the broader quest to comprehend planetary habitability, informing strategies to identify promising targets for future observational campaigns.
Subject of Research: Evolution and volatile retention of molten super-Earth exoplanet L 98-59 d through coupled interior-atmosphere modeling.
Article Title: Volatile-rich evolution of molten super-Earth L 98-59 d.
Article References:
Nicholls, H., Lichtenberg, T., Chatterjee, R.D. et al. Volatile-rich evolution of molten super-Earth L 98-59 d. Nat Astron (2026). https://doi.org/10.1038/s41550-026-02815-8
Image Credits: AI Generated
