Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2024 Aug;632(8023):50-54.
doi: 10.1038/s41586-024-07688-3. Epub 2024 Jul 17.

A hot-Jupiter progenitor on a super-eccentric retrograde orbit

Affiliations

A hot-Jupiter progenitor on a super-eccentric retrograde orbit

Arvind F Gupta et al. Nature. 2024 Aug.

Abstract

Giant exoplanets orbiting close to their host stars are unlikely to have formed in their present configurations1. These 'hot Jupiter' planets are instead thought to have migrated inward from beyond the ice line and several viable migration channels have been proposed, including eccentricity excitation through angular-momentum exchange with a third body followed by tidally driven orbital circularization2,3. The discovery of the extremely eccentric (e = 0.93) giant exoplanet HD 80606 b (ref. 4) provided observational evidence that hot Jupiters may have formed through this high-eccentricity tidal-migration pathway5. However, no similar hot-Jupiter progenitors have been found and simulations predict that one factor affecting the efficacy of this mechanism is exoplanet mass, as low-mass planets are more likely to be tidally disrupted during periastron passage6-8. Here we present spectroscopic and photometric observations of TIC 241249530 b, a high-mass, transiting warm Jupiter with an extreme orbital eccentricity of e = 0.94. The orbit of TIC 241249530 b is consistent with a history of eccentricity oscillations and a future tidal circularization trajectory. Our analysis of the mass and eccentricity distributions of the transiting-warm-Jupiter population further reveals a correlation between high mass and high eccentricity.

PubMed Disclaimer

Conflict of interest statement

The authors declare no competing interests.

Figures

Fig. 1
Fig. 1. TIC 241249530 b transit measurements.
a, TESS photometric measurements (blue), shown at their native 30-min cadence, and the best-fit transit model (black curve) and 3σ confidence region (grey). b, Residuals to the best-fit model for the TESS transit detection. c, Diffuser-assisted ARCTIC photometric measurements (red) and the best-fit transit model (black curve) and 3σ confidence region (grey). We show both the raw 30-s cadence and binned 15-min cadence measurements, along with the residuals to the fitted transit signal. d, Residuals to the best-fit model for the ARCTIC transit detection. All brightnesses are given in parts per thousand (ppt). Error bars on individual data points indicate the 1σ measurement uncertainties. Source Data
Fig. 2
Fig. 2. Phase-folded RV measurements for TIC 241249530.
a, RV measurements from NEID (blue), HPF (black), HARPS-N (yellow) and best-fit orbit model (black curve). b, Residuals to the RV orbit fit. c, Best-fit Rossiter–McLaughlin model (solid curve) and 1σ confidence region (grey), with the signal for aligned orbit shown for comparison (dashed curve). d, Residuals to the Rossiter–McLaughlin model fit. The dashed red box in a highlights the in-transit RV measurements, which are shown in b. Error bars on individual data points indicate the 1σ measurement uncertainties. Source Data
Fig. 3
Fig. 3. Mass–eccentricity distribution for transiting warm Jupiters.
a, Best-fit beta distributions for transiting warm Jupiters more massive (red) and less massive (blue) than 1.935 MJ. The shaded regions represent the 1σ (darkest), 2σ and 3σ (lightest) posteriors on each fit. b, Masses and eccentricities of transiting warm Jupiters. The super-eccentric hot-Jupiter progenitors HD 80606 b and TIC 241249530 b are labelled with upright and inverted stars, respectively. The horizontal dashed line indicates the median mass of the population, which is the threshold used for the fits shown in a. Our results are insensitive to the exact value of the threshold (see Methods); the median is chosen simply for visualization purposes. Error bars on individual data points indicate the 1σ uncertainties from the literature. Source Data
Extended Data Fig. 1
Extended Data Fig. 1. TESS light curves for TIC 241249530.
a, Sector 20 TESS-SPOC data (30-min cadence). The vertical dashed line marks the best-fit transit midpoint. b, Sector 60 SPOC data (2-min cadence). Error bars on individual points represent 1σ measurement uncertainties.
Extended Data Fig. 2
Extended Data Fig. 2. Reconstructed NESSI speckle images and 5σ contrast curves for TIC 241249530.
Observations were taken simultaneously at 562 nm with the blue camera (upper-left inset image) and at 832 nm with the red camera (upper-right inset image). The contrast curves indicate the limiting magnitude difference at which bound or background sources could be detected for separations between 0.2″ and 1.2″.
Extended Data Fig. 3
Extended Data Fig. 3. RV time series for TIC 241249530.
a, RV measurements from NEID (blue), HPF (black) and HARPS-N (gold) and best-fit orbit model (black curve). b, Residuals to the RV orbit fit. Vertical red lines in both panels mark the predicted transit times. The grey-shaded region bounds the 3σ confidence intervals for the fit. The horizontal axis is in units of days relative to Barycentric Julian Date (BJD) 2457000. Error bars on individual points represent 1σ measurement uncertainties.
Extended Data Fig. 4
Extended Data Fig. 4. SED for TIC 241249530.
Red symbols represent the observed photometric measurements, for which the vertical error bars represent the 1σ measurement uncertainties and the horizontal bars represent the effective width of the passband. Blue symbols are the model fluxes from the best-fit PHOENIX atmosphere model (black). The Gaia spectrophotometry is represented as a grey swathe (see also inset plot).
Extended Data Fig. 5
Extended Data Fig. 5. Constraints on the parameters of the first vZLK oscillation for TIC 241249530 b.
a, Constraints on initial inclination. b, Constrains on initial maximum eccentricity. The red region (0 au < ai < 4.2 au) indicates the absence of vZLK oscillations. The orange region (4.2 au < ai < 7 au) indicates the presence of vZLK oscillations but with insufficient ei,max to reach the present-day orbit owing to strong short-range forces (primarily GR, with further contributions from tides and rotational distortions). The blue region (ai > 7 au) indicates the presence of vZLK oscillations that can drive the planet to reach present-day parameters.
Extended Data Fig. 6
Extended Data Fig. 6. Simulated evolution of the orbit of TIC 241249530 b resulting from high-eccentricity migration driven by vZLK oscillations.
a, Evolution of the eccentricity of the planetary orbit over time. b, Evolution of the mutual inclination between the planet and binary orbits over time. c, Evolution of the semimajor axis, a, and periastron separation, a(1 − e), of the planetary orbit over time. The vertical red lines at 3.2 Gyr mark the age at which the orbit reaches the present-day conditions (e = 0.94, a = 0.64 au). We adopt ai = 10 au, ei = 0.1 for the initial orbital parameters for illustrative purposes.

References

    1. Dawson, R. I. & Johnson, J. A. Origins of hot Jupiters. Annu. Rev. Astron. Astrophys.56, 175–221 (2018).
    1. Holman, M., Touma, J. & Tremaine, S. Chaotic variations in the eccentricity of the planet orbiting 16 Cygni B. Nature386, 254–256 (1997).
    1. Fabrycky, D. & Tremaine, S. Shrinking binary and planetary orbits by Kozai cycles with tidal friction. Astrophys. J.669, 1298–1315 (2007).
    1. Naef, D. et al. HD 80606 b, a planet on an extremely elongated orbit. Astron. Astrophys.375, L27–L30 (2001).
    1. Wu, Y. & Murray, N. Planet migration and binary companions: the case of HD 80606b. Astrophys. J.589, 605–614 (2003).

LinkOut - more resources