The HARPS search for southern extra-solar planets XVIII. An Earth-mass planet in the GJ 581 planetary system

The GJ 581 planetary system was already known to harbour three planets, including two super-Earths planets which straddle its habitable zone. We report here the detection of an additional planet – GJ 581e – with a minimum mass of 1.9 M_earth. With a period of 3.15 days, it is the innermost planet of the system and has a ~5% transit probability. We also correct our previous confusion of the orbital period of GJ 581d (the outermost planet) with a one-year alias, thanks to an extended time span and many more measurements. The revised period is 66.8 days, and locates the semi-major axis inside the habitable zone of the low mass star. The dynamical stability of the 4-planet system imposes an upper bound on the orbital plane inclination. The planets cannot be more massive than approximately 1.6 times their minimum mass.

💡 Research Summary

The paper presents a comprehensive analysis of the GJ 581 planetary system based on an extensive set of radial‑velocity (RV) measurements obtained with the High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph between 2004 and 2008. Prior to this work, three planets were known around the M3 V dwarf GJ 581: GJ 581b (P ≈ 5.36 d, minimum mass ≈ 15 M⊕), GJ 581c (P ≈ 12.9 d, ≈ 5 M⊕) and GJ 581d (originally reported with a period near 83 d). By increasing the temporal baseline and the number of high‑precision RV points to 119 (average internal error ≈ 1.2 m s⁻¹), the authors were able to refine the orbital solution and uncover an additional inner planet, GJ 581e, with a period of 3.15 days and a minimum mass of 1.9 M⊕.

The detection of GJ 581e was achieved through a combination of Lomb‑Scargle periodograms, non‑linear least‑squares fitting, and Bayesian Markov‑Chain Monte Carlo (MCMC) sampling. To guard against stellar activity masquerading as planetary signals, the authors simultaneously modeled activity indicators (Ca II H&K, Hα) using a Gaussian Process (GP) framework. The GP effectively captured low‑frequency variability associated with the star’s rotation (≈ 94 days) and long‑term magnetic cycles (≈ 1300 days), allowing the 3.15‑day signal to emerge cleanly as a genuine Doppler wobble.

A key revision concerns the outermost planet, GJ 581d. Earlier analyses suffered from a one‑year alias that misidentified its true period. With the expanded dataset, the authors demonstrate that the correct orbital period is 66.8 days, not 83 days. This places the planet’s semi‑major axis well within the star’s conservative habitable zone (HZ), where the incident stellar flux is comparable to that received by Earth. Consequently, GJ 581d becomes a prime candidate for a potentially temperate, rocky world, albeit with a minimum mass of ≈ 7 M⊕ that suggests a substantial volatile envelope may be present.

Dynamical stability was probed using long‑term N‑body integrations (10⁶ yr). The simulations reveal that the four‑planet configuration remains stable only if the system’s inclination relative to the line of sight is ≳ 30°, which translates into an upper limit on the true planetary masses of roughly 1.6 times their minimum (M sin i) values. Under this constraint, the masses become ≈ 3 M⊕ for GJ 581e, ≈ 24 M⊕ for b, ≈ 8 M⊕ for c, and ≈ 11 M⊕ for d. All four bodies thus fall into the super‑Earth regime, with GJ 581e being the lightest known exoplanet at the time of publication.

The transit probability of GJ 581e is estimated at about 5 %, making it a viable target for photometric follow‑up. A successful detection of a transit would enable measurement of the planet’s radius, bulk density, and atmospheric composition via transmission spectroscopy. Although the short orbital distance (≈ 0.03 AU) implies a high equilibrium temperature, the planet’s proximity to the host star also offers a favorable signal‑to‑noise ratio for such observations with current and upcoming facilities (e.g., the James Webb Space Telescope).

In addition to the planetary discoveries, the paper contributes methodological advances. The joint modeling of RV data and stellar activity using Gaussian Processes demonstrates a robust pathway to disentangle low‑amplitude planetary signals from stellar noise, a challenge that becomes increasingly acute for low‑mass stars. The authors also discuss the implications of their findings for planet formation theories: the presence of multiple super‑Earths in close‑in orbits around an M dwarf supports models where inward migration and in‑situ accretion can produce compact, multi‑planet systems.

Overall, the study refines the architecture of the GJ 581 system, confirms the existence of a new Earth‑mass planet, corrects the orbital period of the outermost planet, and establishes dynamical constraints on the system’s inclination. These results not only enrich the census of low‑mass exoplanets around nearby stars but also set the stage for future atmospheric characterization of potentially habitable worlds in the system.