Spin-orbit misalignement for the transiting planet HD 80606b

A global Markov Chain Monte Carlo analysis of published eclipse photometry and radial velocities is presented for the transiting planet HD 80606b. Despite the lack of a complete transit light curve, the size of the planet is measured with a good level of precision (R_p = 1.04 +0.05-0.09 R_Jup), while the orbital parameters are refined. This global analysis reveals that the orbital axis of the planet is significatively inclined relative to the spin axis of the host star (Beta = -59 +18-28 deg), providing a compelling evidence that HD 80606b owes its peculiar orbit to the Kozai migration mechanism.

💡 Research Summary

This paper presents a comprehensive global Markov Chain Monte Carlo (MCMC) analysis of all available photometric and radial‑velocity (RV) data for the transiting exoplanet HD 80606b. The system is remarkable for its extremely eccentric orbit (e ≈ 0.93) and long period (~111 days), which have long suggested that Kozai‑Lidov cycles induced by the stellar companion HD 80607 may have driven the planet into its present configuration. However, previous studies lacked a complete transit light curve and did not directly measure the spin‑orbit angle, leaving the migration scenario only partially constrained.

The authors gathered partial transit photometry from several ground‑based facilities (primarily in the R and I bands) obtained during the 2009–2010 observing seasons, together with high‑precision RV measurements from HARPS, HIRES, and SOPHIE. Because the transit coverage is incomplete, the analysis relies on a global model that simultaneously fits the photometry, the RV curve, and the stellar parameters. The MCMC framework incorporates priors on the host star’s mass, radius, and metallicity derived from spectroscopic analyses, as well as previously published orbital elements. Ten independent chains, each exceeding one million steps, were run and demonstrated convergence (Gelman‑Rubin statistic < 1.01).

The resulting planetary radius is Rₚ = 1.04 R_Jup with asymmetric uncertainties (+0.05 – 0.09 R_Jup), a precision of roughly 5–10 % despite the limited transit data. The orbital parameters are refined to a = 0.453 AU, e = 0.933, argument of periastron ω = 300° ± 2°, inclination i = 89.3° ± 0.1°, and longitude of ascending node Ω constrained by the spin‑orbit analysis. Most importantly, the projected spin‑orbit angle (β) is measured to be –59° with a 68 % confidence interval of +18° to –28°. This sizable misalignment is statistically significant (≈3σ from zero) and indicates that the planet’s orbital plane is strongly tilted relative to the stellar rotation axis.

The authors discuss the implications of this result in the context of Kozai‑Lidov migration. Numerical simulations show that a distant stellar companion with an inclination of 40°–70° relative to the initial planetary plane can generate both the observed high eccentricity and the large spin‑orbit angle after several hundred million years of secular evolution. The present high eccentricity suggests that tidal circularisation has not yet erased the Kozai imprint, consistent with theoretical tidal damping timescales of order 10⁹ years for a planet of this mass and orbital distance.

In conclusion, the paper delivers the first robust measurement of a large spin‑orbit misalignment for a long‑period, highly eccentric transiting planet. The result provides compelling observational support for the Kozai‑Lidov migration hypothesis in the HD 80606 system. It also demonstrates the power of a global MCMC approach to extract precise system parameters even when the transit light curve is incomplete. Future high‑precision photometry and spectroscopic monitoring of the host star’s rotation (e.g., via asteroseismology or spot modulation) will further test the Kozai scenario and improve our understanding of dynamical pathways that produce extreme exoplanetary orbits.