Bulk composition of the transiting hot Neptune around GJ 436

The hot Neptune orbiting around GJ 436 is a unique example of an intermediate mass planet. Its close-in orbit suggests that the planet has undergone migration and its study is fundamental to understanding planet formation and evolution. As it transits its parent star, it is the only Neptune-mass extrasolar planet of known mass and radius, being slightly larger and more massive than Neptune (M=22.6 M_Earth, R=4.19R_Earth). In this regime, several bulk compositions are possible: from an Earth-like core with a thick hydrogen envelope to a water-rich planet with a thin hydrogen envelope comprising a Neptune-like structure. We combine planet-structure modeling with an advanced planet-formation model to assess the likelihood of the different possible bulk compositions of GJ 436 b. We find that both an envelope-free water planet (“Ocean planet”) as well as a diminute version of a gaseous giant planet are excluded. Consisting of a rocky core with a thick hydrogen/helium envelope, a “dry” composition produces not only too small a radius but is also a very unlikely outcome of planet formation around such a low-mass star. We conclude that GJ 436 b is probably of much higher rock content than Neptune (more than 45% in mass), with a small H-He envelope (10 - 20% in mass). This is the expected outcome of the gathering of materials during the migration process in the inner disk, creating a population of which the hot Neptune is representative.

💡 Research Summary

The paper presents a comprehensive investigation of the bulk composition and formation history of GJ 436b, the only known transiting Neptune‑mass exoplanet with both measured mass (22.6 M⊕) and radius (4.19 R⊕). Because its mass and radius place it in a regime where many internal configurations are mathematically possible, the authors combine detailed interior structure modeling with a sophisticated planet‑formation simulation to evaluate which configurations are physically plausible.

The interior model treats the planet as a three‑layer object: (1) a rocky/iron core, (2) a high‑pressure water/ice mantle, and (3) a hydrogen‑helium envelope. For each layer the authors adopt state‑of‑the‑art equations of state (EOS) that incorporate recent high‑pressure laboratory data and theoretical predictions, especially for water at megabar pressures where super‑ionic and plasma phases appear. By varying the mass fractions of the three layers they generate a family of models that reproduce the observed mass–radius point. This exercise reveals a strong degeneracy: a planet with a massive rocky core (≈70 % of the total mass) and a thin H‑He envelope (≈10 % mass) can have the same radius as a planet with a smaller core (≈45 % mass), a substantial water layer (≈35 % mass), and a thicker envelope (≈20 % mass). Consequently, a purely structural approach cannot uniquely determine the composition.

To break this degeneracy the authors embed the interior models within a formation framework. They simulate the growth of a planetary embryo in a protoplanetary disk around a low‑mass M dwarf (the host star GJ 436). The simulation follows three main phases: (i) solid accretion of planetesimals (rock, iron, and water ice) that builds the core, (ii) runaway gas accretion once the core reaches a critical mass, and (iii) Type I/II migration driven by disk torques. The disk model reflects the reduced solid surface density expected around an M dwarf, and the temperature profile determines where water ice is present. As the embryo migrates inward, it traverses regions with decreasing solid‑to‑gas ratios, which limits further water delivery and gas capture.

The formation results show that, under realistic disk lifetimes (a few Myr) and migration rates, embryos that end up at the observed 0.03 AU orbit typically acquire a core that is at least 45 % of the final planetary mass and are dominated by refractory material (rock and iron). Water delivery during migration is modest, contributing at most ~35 % of the total mass, while the hydrogen‑helium envelope remains limited to 10‑20 % of the mass because the gas accretion phase is truncated either by rapid inward migration or by early disk dispersal.

By comparing the outcomes of the formation simulations with the families of interior models, the authors can discard two extreme scenarios: (a) an “ocean planet” composed almost entirely of water with a negligible gas envelope, which is inconsistent with the limited water delivery predicted by the migration pathway; and (b) a “dry” Neptune‑like planet with a massive rocky core and an almost absent envelope, which would be too small in radius and is statistically unlikely given the disk’s solid inventory. The most probable configuration is therefore a high‑density core (>45 % of the mass) of rock and iron, a substantial water/ice mantle (≈30‑35 % mass), and a modest H‑He envelope accounting for 10‑20 % of the total mass. This composition is richer in refractory material than the Solar System’s Neptune, reflecting the inner‑disk environment in which the planet assembled.

The study concludes that GJ 436b represents a distinct class of “hot Neptunes” that formed by gathering a large fraction of refractory solids in the inner disk, migrated inward while accreting a limited amount of volatile ices and gas, and finally settled into a close‑in orbit with a comparatively thin hydrogen‑helium envelope. The combined structural‑formation approach not only narrows down the plausible bulk composition of GJ 436b but also provides a framework for interpreting the growing population of transiting Neptune‑mass planets around low‑mass stars.