A Thousand Alien Atmospheres: How Many Worlds Does It Take to See a Trend?

The next big question in exoplanet science may not be what one alien atmosphere contains, but what a thousand of them reveal together.

ESA’s Ariel Space Mission is planned to characterise the atmospheres of roughly 1,000 exoplanets. Its goal is not just to catalog strange worlds one by one, but to quantify population-level trends that encode how planets form and evolve. A new preprint by Wasi M. F. Naqvi and Nicolas B. Cowan at McGill University describes a tool built for that coming flood of data. It is called HERMES—HiERarchical Modelling for Exoplanet Science—and it is designed to probe whether astronomers can detect a three-way relationship between a planet’s mass, its atmospheric metallicity, and the metallicity of its host star.

This is a simulation and methods paper, not an observation of actual exoplanet atmospheres. But it asks a question that will matter when Ariel flies: how large and how diverse must a survey be before a subtle cosmic trend emerges from the noise?

Why Host-Star Metallicity Leaves a Trail

Metallicity is astronomy’s shorthand for the abundance of elements heavier than hydrogen and helium. In this study, the abundance of water in a planetary atmosphere is used as a proxy for atmospheric metallicity. In a star, the usual shorthand is [Fe/H], the iron abundance relative to hydrogen compared with the Sun.

The solar system offers a tantalising hint. Jupiter is enriched in heavy elements relative to hydrogen by about 3 to 6 times solar, while Uranus and Neptune are enriched by roughly 70 to 100 times solar. More massive giant planets tend to be less enriched, because they pull in proportionally more hydrogen and helium during runaway gas accretion.

Does that inverse mass–metallicity trend hold across exoplanets? Earlier transmission spectroscopy studies have found hints, but the picture is debated. Clouds, carbon-to-oxygen ratios, formation history, and updated stellar abundance measurements can all weaken or complicate the signal. Host-star metallicity adds another possible thread: planets forming in a metal-rich disk may have a larger reservoir of heavy elements, possibly leaving a detectable imprint in their atmospheres.

The problem is that planets are not identical. Even at the same mass and around the same kind of star, atmospheric metallicities can vary because of different accretion histories, migration paths, atmospheric mixing, and cloud properties. That planet-to-planet spread is called intrinsic astrophysical scatter. Disentangling it from measurement noise is hard. HERMES is built to do that job across multiple dimensions at once.

Building Fake Surveys to Test a Real Idea

Naqvi and Cowan began with the Ariel Mission Candidate Sample, a curated list of planets suitable for Ariel observations. From 977 confirmed planets, they selected 858 that have measured host-star metallicities and physically valid planetary-mass uncertainty bounds.

They then injected a plausible relationship into those planets: atmospheric water abundance—used as a proxy for atmospheric metallicity—decreases with planetary mass, increases with host-star metallicity, and includes random intrinsic scatter. They generated 180 mock surveys with different sample sizes, from 50 to 600 planets, and different mass ranges. Finally, they fitted a hierarchical Bayesian model to each fake survey.

Hierarchical Bayesian models are well suited to this problem because they infer both the individual planets and the population trend simultaneously, while accounting for uncertainty. The team compared a 3D model, which includes host-star metallicity, with a 2D model, which omits it. If the 3D model predicts the mock data better, that suggests the stellar signal can be separated from intrinsic scatter.

Leverage Is a Compass, Not a Luxury

The first result is reassuring: HERMES recovered the injected trends across the full range of survey designs tested. Its posterior distributions were also well calibrated, meaning the model was neither overconfident nor overly cautious about its own uncertainties.

The second result concerns survey “leverage.” In this context, leverage is the spread of planets along an axis of diversity; the paper also uses a normalized version that accounts for measurement uncertainties. A survey with planets spanning a wide range of masses has more mass leverage. A survey with a wide range of host-star metallicities has more stellar-metallicity leverage.

The team found that leverage along the relevant axis remains a reliable predictor of how precisely a slope can be measured. For the planetary mass–metallicity slope, mass leverage is the best predictor of precision. The stellar-metallicity slope improves with stellar-metallicity information, but more weakly, because the range of [Fe/H] values is narrower and its measurement uncertainties are often large.

Sample size plays a different role. Increasing the number of planets lowers the overall uncertainty for every parameter. But for the planetary mass–metallicity slope, adding mass leverage at fixed sample size is more efficient than simply adding more planets of the same mass range. For the intercept and intrinsic scatter, sample size is the primary control.

When Planet-to-Planet Scatter Takes Over

The most consequential finding is a threshold. If intrinsic scatter is small, even moderate surveys can recover the correlation between stellar and planetary metallicity. But as scatter grows, small surveys lose the signal first.

For surveys with 150 planets or fewer, sensitivity drops once intrinsic scatter reaches roughly 0.8 to 1.0 dex. For an Ariel-scale Tier 2 transit survey of at least 400 planets, HERMES robustly recovers the stellar–planetary metallicity correlation even when intrinsic scatter is as large as 1.2 dex. At 600 planets, the 3D model was favoured in every mock survey tested across the full scatter range, up to 2.0 dex.

That has a counterintuitive implication for survey design. The paper argues that the breadth of Ariel’s Tier 1 survey may be more valuable for population-level trend detection than the depth of a smaller Tier 2 survey, because additional planets contribute leverage on multiple axes at once. In multiple dimensions, optimising for one kind of diversity does not guarantee another.

What HERMES Has Not Yet Seen

These results are forecasts, not detections. The trends were injected by the researchers, not measured from real planets. The atmospheric water abundance was assigned an idealised precision of 0.2 dex, and planetary mass uncertainties were not explicitly propagated in the baseline model. Host-star metallicities were treated with symmetric Gaussian uncertainties, and the z-score reference values came from a population-level fit to the full candidate sample.

Intrinsic scatter is also a catch-all term. It can absorb genuine planet-to-planet diversity as well as any unmodelled physical effect, such as equilibrium temperature, age, irradiation history, cloud properties, or formation location. The paper argues that if the true astrophysical scatter exceeds current estimates, Ariel’s planned sample size should still provide enough statistical power, though smaller precursor surveys may struggle unless they target low-scatter populations or achieve unusually high stellar-metallicity leverage.

Still, the paper establishes HERMES as a practical tool for survey design and science-yield forecasting. It suggests that Ariel’s planned breadth should give astronomers a fighting chance to distinguish a stellar-metallicity signal from intrinsic planet-to-planet scatter—provided stellar metallicities are measured precisely and homogeneously.

The Long View: A Murmur of Worlds

If HERMES is right, the signature of planet formation may not be written in any single atmosphere. It may live in the statistical murmur of hundreds: a slight tilt in a graph, a correlation that only becomes visible when enough worlds are compared. Ariel’s survey is being designed to listen for that murmur. The new work shows how hard that will be—and why sample size becomes decisive once planet-to-planet scatter is large.