High-eccentricity Migration Research Articles

High-eccentricity Migration with Disk-induced Spin–Orbit Misalignment: A Preference for Perpendicular Hot Jupiters

High-eccentricity migration is a likely formation mechanism for many observed hot Jupiters, particularly those with a large misalignment between the stellar spin axis and orbital angular momentum axis of the planet. In one version of high-eccentricity migration, an inclined stellar companion excites von Zeipel–Lidov–Kozai (ZLK) eccentricity oscillations of a cold Jupiter, and tidal dissipation causes the planet’s orbit to shrink and circularize. Throughout this process, the stellar spin can evolve chaotically, resulting in highly misaligned hot Jupiters (HJs). Previous population studies of this migration mechanism have assumed that the stellar spin is aligned with the planetary orbital angular momentum when the companion begins to induce ZLK oscillations. However, in the presence of a binary companion, the star’s obliquity may be significantly excited during the dissipation of its protoplanetary disk. We calculate the stellar obliquities produced in the protoplanetary disk phase and use these to perform an updated population synthesis of ZLK-driven high-eccentricity migration with an F-type host star. We find that the resulting obliquity distribution of HJ systems is predominantly retrograde with a broad peak near 90°. The distribution we obtain has intriguing similarities to the recently observed preponderance of perpendicular planets close to their host stars.

DREAM

The distribution of close-in exoplanets is shaped by a complex interplay between atmospheric and dynamical processes. The Desert-Rim Exoplanets Atmosphere and Migration (DREAM) program aims at disentangling those processes through the study of the hot Neptune desert, whose rim hosts planets that are undergoing, or survived, atmospheric evaporation and orbital migration. In this first paper, we use the Rossiter-McLaughlin revolutions (RMR) technique to investigate the orbital architecture of 14 close-in planets ranging from mini-Neptune to Jupiter-size and covering a broad range of orbital distances. While no signal is detected for the two smallest planets, we were able to constrain the sky-projected spin-orbit angle of six planets for the first time, to revise its value for six others, and, thanks to constraints on the stellar inclination, to derive the 3D orbital architecture in seven systems. These results reveal a striking three-quarters of polar orbits in our sample, all being systems with a single close-in planet but of various stellar and planetary types. High-eccentricity migration is favored to explain such orbits for several evaporating warm Neptunes, supporting the role of late migration in shaping the desert and populating its rim. Putting our measurements in the wider context of the close-in planet population will be useful to investigate the various processes shaping their architectures.

Long-term evolution of exoplanet systems

AbstractIn this Review, I discuss recent developments on the long-term dynamical evolution of exoplanet systems, focusing on how distinctive dynamical processes may have shaped the orbital architectures of observed populations. I include three applications that highlight part of my own work. First, I examine the high-eccentricity tidal migration of hot Jupiters from a phase of dynamical instability and subsequent secular interactions in two-planet systems. Second, secular chaos as the origin of ultra-short-period planets with extreme period ratios. Third, secular resonance sweeping driven by a dispersing protoplanetary disk as the origin hot Neptunes residing in polar orbits. Finally, I discuss how upcoming observations will allow further constraining the prevalence of these dynamical processes.

An old warm Jupiter orbiting the metal-poor G-dwarf TOI-5542

We report the discovery of a 1.32−0.10+0.10 MJup planet orbiting on a 75.12 day period around the G3V 10.8−3.6+2.1 Gyr old star TOI-5542 (TIC 466206508; TYC 9086-1210-1). The planet was first detected by the Transiting Exoplanet Survey Satellite (TESS) as a single transit event in TESS Sector 13. A second transit was observed 376 days later in TESS Sector 27. The planetary nature of the object has been confirmed by ground-based spectroscopic and radial velocity observations from the CORALIE and HARPS spectrographs. A third transit event was detected by the ground-based facilities NGTS, EulerCam, and SAAO. We find the planet has a radius of 1.009−0.035+0.036 RJup and an insolation of 9.6−0.8+0.9 S⊕, along with a circular orbit that most likely formed via disk migration or in situ formation, rather than high-eccentricity migration mechanisms. Our analysis of the HARPS spectra yields a host star metallicity of [Fe/H] = −0.21 ± 0.08, which does not follow the traditional trend of high host star metallicity for giant planets and does not bolster studies suggesting a difference among low- and high-mass giant planet host star metallicities. Additionally, when analyzing a sample of 216 well-characterized giant planets, we find that both high masses (4 MJup < Mp < 13 MJup) and low masses (0.5 MJup < Mp < 4 MJup), as well as both both warm (P > 10 days) and hot (P < 10 days) giant planets are preferentially located around metal-rich stars (mean [Fe/H] > 0.1). TOI-5542b is one of the oldest known warm Jupiters and it is cool enough to be unaffected by inflation due to stellar incident flux, making it a valuable contribution in the context of planetary composition and formation studies.

Making hot Jupiters in stellar clusters: The importance of binary exchange

ABSTRACT It has been suggested that the occurrence rate of hot Jupiters (HJs) in open clusters might reach several per cent, significantly higher than that of the field (∼a per cent). In a stellar cluster, when a planetary system scatters with a stellar binary, it may acquire a companion star, which may excite large-amplitude von Zeipel–Lidov–Kozai oscillations in the planet’s orbital eccentricity, triggering high-eccentricity migration, and the formation of an HJ. We quantify the efficiency of this mechanism by modelling the evolution of a gas giant around a solar mass star under the influence of successive scatterings with binary and single stars. We show that the chance that a planet ∈ (1, 10) au becomes an HJ in a Gyr in a cluster of stellar density n* = 50 pc−3, and binary fraction fbin = 0.5 is about 2 per cent and an additional 4 per cent are forced by the companion star into collision with or tidal disruption by the central host. An empirical fit shows that the total percentage of those outcomes asymptotically reaches an upper limit determined solely by fbin (e.g. 10 per cent at fbin = 0.3 and 18 per cent at fbin = 1) on a time-scale inversely proportional to n* (∼Gyr for n* ∼ 100 pc−3). The ratio of collisions to tidal disruptions is roughly a few, and depends on the tidal model. Therefore, if the giant planet occurrence rate is 10 per cent, our mechanism implies an HJ occurrence rate of a few times 0.1 per cent in a Gyr and can thus explain a substantial fraction of the observed rate.

Production of hot Jupiter candidates from high-eccentricity mechanisms for different initial planetary mass configurations

ABSTRACT Hot Jupiters (HJs) are giant planets with orbital periods of the order of a few days with semimajor axis within ∼0.1 au. Several theories have been invoked in order to explain the origin of this type of planets, one of them being the high-eccentricity migration. This migration can occur through different high-eccentricity mechanisms. Our investigation focused on six different kinds of high-eccentricity mechanisms, namely, direct dispersion, coplanar, Kozai–Lidov, secular chaos, E1 and E2 mechanisms. We investigated the efficiency of these mechanisms for the production of HJ candidates in multiplanet systems initially tightly-packed in the semimajor axis, considering a large set of numerical simulations of the exact equations of motion in the context of the N-body problem. In particular, we analyzed the sensitivity of our results to the initial number of planets, the initial semimajor axis of the innermost planetary orbit, the initial configuration of planetary masses, and to the inclusion of general relativity (GR) effects. We found that the E1 mechanism is the most efficient in producing HJ candidates both in simulations with and without the contribution of GR, followed by the Kozai–Lidov and E2 mechanisms. Our results also revealed that, except for the initial equal planetary mass configuration, the E1 mechanism was notably efficient in the other initial planetary mass configurations considered in this work. Finally, we investigated the production of HJ candidates with prograde, retrograde, and alternating orbits. According to our statistical analysis, the Kozai–Lidov mechanism has the highest probability of significantly exciting the orbital inclinations of the HJ candidates.

The Discovery of a Planetary Companion Interior to Hot Jupiter WASP-132 b

Hot Jupiters are generally observed to lack close planetary companions, a trend that has been interpreted as evidence for high-eccentricity migration. We present the discovery and validation of WASP-132 c (TOI-822.02), a 1.85 ± 0.10 R ⊕ planet on a 1.01 day orbit interior to the hot Jupiter WASP-132 b. Transiting Exoplanet Survey Satellite and ground-based follow-up observations, in conjunction with vetting and validation analysis, enable us to rule out common astrophysical false positives and validate the observed transit signal produced by WASP-132 c as a planet. Running the validation tools vespa and TRICERATOPS on this signal yields false-positive probabilities of 9.02 × 10−5 and 0.0107, respectively. Analysis of archival CORALIE radial velocity data leads to a 3σ upper limit of 28.23 ms−1 on the amplitude of any 1.01 day signal, corresponding to a 3σ upper mass limit of 37.35 M ⊕. Dynamical simulations reveal that the system is stable within the 3σ uncertainties on the planetary and orbital parameters for timescales of ∼100 Myr. The existence of a planetary companion near the hot Jupiter WASP-132 b makes the giant planet’s formation and evolution via high-eccentricity migration highly unlikely. Being one of just a handful of nearby planetary companions to hot Jupiters, WASP-132 c carries with it significant implications for the formation of the system and hot Jupiters as a population.

Nodal Precession and Tidal Evolution of Two Hot Jupiters: WASP-33 b and KELT-9 b

Hot Jupiters orbiting rapidly rotating stars on inclined orbits undergo tidally induced nodal precession measurable over several years of observations. The Hot Jupiters WASP-33 b and KELT-9 b are particularly interesting targets because they are among the hottest planets found to date, orbiting relatively massive stars. Here, we analyze archival and new data that span 11 and 5 yr for WASP-33 b and KELT-9 b, respectively, in order to model and improve upon their tidal precession parameters. Our work confirms the nodal precession for WASP-33 b and presents the first clear detection of the precession of KELT-9 b. We determine that WASP-33 and KELT-9 have gravitational quadrupole moments and , respectively. We estimate the planets’ precession periods to be yr and yr, respectively, and that they will cease to transit their host stars around the years CE and CE, respectively. Additionally, we investigate both planets’ tidal and orbital evolution, suggesting that a high-eccentricity tidal migration scenario is possible to produce both system architectures and that they will most likely not be engulfed by their hosts before the end of their main-sequence lifetimes.

The Discovery of a Planetary Companion Interior to Hot Jupiter WASP-132 b

Hot Jupiters are generally observed to lack close planetary companions, a trend that has been interpreted as evidence for high-eccentricity migration. We present the discovery and validation of WASP-132 c (TOI-822.02), a 1.85 $\pm$ 0.10 $R_{\oplus}$ planet on a 1.01 day orbit interior to the hot Jupiter WASP-132 b. Transiting Exoplanet Survey Satellite (TESS) and ground-based follow-up observations, in conjunction with vetting and validation analysis, enable us to rule out common astrophysical false positives and validate the observed transit signal produced by WASP-132 c as a planet. Running the validation tools \texttt{vespa} and \texttt{triceratops} on this signal yield false positive probabilities of $9.02 \times 10^{-5}$ and 0.0107, respectively. Analysis of archival CORALIE radial velocity data leads to a 3$\sigma$ upper limit of 28.23 ms$^{-1}$ on the amplitude of any 1.01-day signal, corresponding to a 3$\sigma$ upper mass limit of 37.35 $M_{\oplus}$. Dynamical simulations reveal that the system is stable within the 3$\sigma$ uncertainties on planetary and orbital parameters for timescales of $\sim$100 Myr. The existence of a planetary companion near the hot Jupiter WASP-132 b makes the giant planet's formation and evolution via high-eccentricity migration highly unlikely. Being one of just a handful of nearby planetary companions to hot Jupiters, WASP-132 c carries with it significant implications for the formation of the system and hot Jupiters as a population.

Secular chaos in white dwarf planetary systems: origins of metal pollution and short-period planetary companions

ABSTRACT Secular oscillations in multiplanet systems can drive chaotic evolution of a small inner body through non-linear resonant perturbations. This ‘secular chaos’ readily pushes the inner body to an extreme eccentricity, triggering tidal interactions or collision with the central star. We present a numerical study of secular chaos in systems with two planets and test particles using the ring-averaging method, with emphasis on the relationship between the planets’ properties and the time-scale and efficiency of chaotic diffusion. We find that secular chaos can excite extreme eccentricities on time-scales spanning several orders of magnitude in a given system. We apply our results to the evolution of planetary systems around white dwarfs (WDs), specifically the tidal disruption and high-eccentricity migration of planetesimals and planets. We find that secular chaos in a planetesimal belt driven by large (≳10 M⊕), distant ($\gtrsim 10 \, \mathrm{au}$) planets can sustain metal accretion on to a WD over Gyr time-scales. We constrain the total mass of planetesimals initially present within the chaotic zone by requiring that the predicted mass delivery rate to the Roche limit be consistent with the observed metal accretion rates of WDs with atmospheric pollution throughout the cooling sequence. Based on the occurrence of long-period exoplanets and exo-asteroid belts, we conclude that secular chaos can be a significant (perhaps dominant) channel for polluting solitary WDs. Secular chaos can also produce short-period planets and planetesimals around WDs in concert with various circularization mechanisms. We discuss prospects for detecting exoplanets driving secular chaos around WDs using direct imaging and microlensing.

A low-eccentricity migration pathway for a 13-h-period Earth analogue in a four-planet system

It is commonly accepted that exoplanets with orbital periods shorter than 1 day, also known as ultra-short period (USP) planets, formed further out within their natal protoplanetary disk, before migrating to their current-day orbits via dynamical interactions. One of the most accepted theories suggests a violent scenario involving high-eccentricity migration followed by tidal circularization. Here, we present the discovery of a four planet system orbiting the bright (V=10.5) K6 dwarf star TOI-500. The innermost planet is a transiting, Earth-sized USP planet with an orbital period of $\sim$ 13 hours, a mass of 1.42 $\pm$ 0.18 M$_{\oplus}$, a radius of $1.166^{0.061}_{-0.058}$ R$_{\oplus}$, and a mean density of 4.89$^{+1.03}_{-0.88}$ gcm$^{-3}$. Via Doppler spectroscopy, we discovered that the system hosts three outer planets on nearly circular orbits with periods of 6.6, 26.2, and 61.3d and minimum masses of 5.03 $\pm$ 0.41 M$_{\oplus}$, 33.12 $\pm$ 0.88 M$_{\oplus}$ and 15.05$^{+1.12}_{-1.11}$ M$_{\oplus}$, respectively. The presence of both a USP planet and a low-mass object on a 6.6-day orbit indicates that the architecture of this system can be explained via a scenario in which the planets started on low-eccentricity orbits, then moved inwards through a quasi-static secular migration. Our numerical simulations show that this migration channel can bring TOI-500 b to its current location in 2 Gyrs, starting from an initial orbit of 0.02au. TOI-500 is the first four planet system known to host a USP Earth analog whose current architecture can be explained via a non-violent migration scenario.

Kepler-1656b’s Extreme Eccentricity: Signature of a Gentle Giant

Highly eccentric orbits are one of the major surprises of exoplanets relative to the solar system and indicate rich and tumultuous dynamical histories. One system of particular interest is Kepler-1656, which hosts a sub-Jovian planet with an eccentricity of 0.8. Sufficiently eccentric orbits will shrink in the semimajor axis due to tidal dissipation of orbital energy during periastron passage. Here our goal was to assess whether Kepler-1656b is currently undergoing such high-eccentricity migration, and to further understand the system’s origins and architecture. We confirm a second planet in the system with M c = 0.40 ± 0.09 M jup and Pc = 1919 ± 27 days. We simulated the dynamical evolution of planet b in the presence of planet c and find a variety of possible outcomes for the system, such as tidal migration and engulfment. The system is consistent with an in situ dynamical origin of planet b followed by subsequent eccentric Kozai–Lidov perturbations that excite Kepler-1656b’s eccentricity gently, i.e., without initiating tidal migration. Thus, despite its high eccentricity, we find no evidence that planet b is or has migrated through the high-eccentricity channel. Finally, we predict the outer orbit to be mutually inclined in a nearly perpendicular configuration with respect to the inner planet orbit based on the outcomes of our simulations and make observable predictions for the inner planet’s spin–orbit angle. Our methodology can be applied to other eccentric or tidally locked planets to constrain their origins, orbital configurations, and properties of a potential companion.

TOI-1670 b and c: An Inner Sub-Neptune with an Outer Warm Jupiter Unlikely to Have Originated from High-eccentricity Migration

We report the discovery of two transiting planets around the bright (V = 9.9 mag) main-sequence F7 star TOI-1670 by the Transiting Exoplanet Survey Satellite. TOI-1670 b is a sub-Neptune ( R ⊕) on a 10.9 day orbit, and TOI-1670 c is a warm Jupiter ( R Jup) on a 40.7 day orbit. Using radial velocity observations gathered with the Tull Coudé Spectrograph on the Harlan J. Smith telescope and HARPS-N on the Telescopio Nazionale Galileo, we find a planet mass of M Jup for the outer warm Jupiter, implying a mean density of g cm−3. The inner sub-Neptune is undetected in our radial velocity data (M b < 0.13 M Jup at the 99% confidence level). Multiplanet systems like TOI-1670 hosting an outer warm Jupiter on a nearly circular orbit () and one or more inner coplanar planets are more consistent with “gentle” formation mechanisms such as disk migration or in situ formation rather than high-eccentricity migration. Of the 11 known systems with a warm Jupiter and a smaller inner companion, eight (73%) are near a low-order mean-motion resonance, which can be a signature of migration. TOI-1670 joins two other systems (27% of this subsample) with period commensurabilities greater than 3, a common feature of in situ formation or halted inward migration. TOI-1670 and the handful of similar systems support a diversity of formation pathways for warm Jupiters.

Orbital alignment of HD 332231 b

Context. Contrary to the orthodox picture of planet formation resulting in a neatly ordered Solar System, exoplanet systems exhibit highly diverse orbits: short and long periods, circular and eccentric, well- and misaligned, and even retrograde orbits. In order to understand this diversity it is essential to probe key orbital parameters. Spin-orbit alignment is such a parameter and can provide information about the formation and migration history of the system. However, tidal circularisation and alignment might hamper interpretations of orbital eccentricity and obliquities in the context of planet formation and evolution for planets on orbits shorter than about 10 days. Aims. Here we aim to measure the projected stellar obliquity in the HD 332231system in which a warm (period ≈ 18.7 days) giant planet orbits a bright F star on a circular orbit. Methods. We observed the system during a transit with the HARPS-N spectrograph and obtained data on the Rossiter-McLaughlin effect. We analysed the spectroscopic transit data together with new TESS photometry employing three different analysis methods. Results. The results from the different approaches are fully consistent. We find a projected obliquity of −2 ± 6°, indicating the stellar spin axis to be well-aligned with the orbit of the planet. We furthermore find evidence for transit timing variations suggesting the presence of an additional third body in the system. Conclusions. Together with the low orbital eccentricity, the good alignment suggests that this warm giant planet has not undergone high-eccentricity migration.

Were the Obliquities in DI Herculis Excited by an Unseen Tertiary Companion?

Abstract The eclipsing binary DI Herculis has garnered interest for several decades because of an apparent disagreement between the observed and calculated values of the apsidal precession rate. The problem was resolved when both stars were found to have high obliquities, but the reason for the high obliquities is unknown. Here, we investigate the possibility that the obliquities are (or were) excited by an unseen tertiary star. Obliquity excitation in the current orbital configuration can be ruled out with existing data; any tertiary star that is sufficiently close or massive to overcome the strong spin–orbit coupling of the binary would have been detected through various dynamical effects. It remains possible that the orbit of DI Herculis was initially wider and the obliquity was excited during high-eccentricity tidal migration driven by a tertiary companion, but in this scenario it would be difficult to explain why the observed spin rates are much faster than the pseudo-synchronous rate. In addition, inward migration is most likely to arise when the mass of the perturbing star is comparable to the binary mass, and such a bright tertiary would have been detected in imaging or spectroscopic data. Alternative explanations that do not invoke a tertiary star should be sought for the large obliquities in DI Herculis.

Origins of Hot Jupiters from the Stellar Obliquity Distribution

Abstract The obliquity of a star, or the angle between its spin axis and the average orbit normal of its companion planets, provides a unique constraint on that system’s evolutionary history. Unlike the solar system, where the Sun’s equator is nearly aligned with its companion planets, many hot-Jupiter systems have been discovered with large spin–orbit misalignments, hosting planets on polar or retrograde orbits. We demonstrate that, in contrast to stars harboring hot Jupiters on circular orbits, those with eccentric companions follow no population-wide obliquity trend with stellar temperature. This finding can be naturally explained through a combination of high-eccentricity migration and tidal damping. Furthermore, we show that the joint obliquity and eccentricity distributions observed today are consistent with the outcomes of high-eccentricity migration, with no strict requirement to invoke the other hot-Jupiter formation mechanisms of disk migration or in situ formation. At a population-wide level, high-eccentricity migration can consistently shape the dynamical evolution of hot-Jupiter systems.

NEID Rossiter–McLaughlin Measurement of TOI-1268b: A Young Warm Saturn Aligned with Its Cool Host Star

Abstract Close-in gas giants present a surprising range of stellar obliquity, the angle between a planet’s orbital axis and its host star’s spin axis. It is unclear whether the obliquities reflect the planets’ dynamical history (e.g., aligned for in situ formation or disk migration versus misaligned for high-eccentricity tidal migration) or whether other mechanisms (e.g., primordial misalignment or planet–star interactions) are more important in sculpting the obliquity distribution. Here we present the stellar obliquity measurement of TOI-1268 (TIC-142394656, V mag ∼ 10.9), a young K-type dwarf hosting an 8.2 day period, Saturn-sized planet. TOI-1268’s lithium abundance and rotation period suggest the system age between the ages of the Pleiades cluster (∼120 Myr) and the Prasepe cluster (∼670 Myr). Using the newly commissioned NEID spectrograph, we constrain the stellar obliquity of TOI-1268 via the Rossiter–McLaughlin effect from both radial velocity and Doppler tomography signals. The 3σ upper bounds of the projected stellar obliquity ∣λ∣ from both models are below 60°. The large host star separation (a/R ⋆ ∼ 17), combined with the system’s young age, makes it unlikely that the planet has realigned its host star. The stellar obliquity measurement of TOI-1268 probes the architecture of a young gas giant beyond the reach of tidal realignment (a/R ⋆ ≳ 10) and reveals an aligned or slightly misaligned system.

A 20 Second Cadence View of Solar-type Stars and Their Planets with TESS: Asteroseismology of Solar Analogs and a Recharacterization of π Men c

Abstract We present an analysis of the first 20 second cadence light curves obtained by the TESS space telescope during its extended mission. We find improved precision of 20 second data compared to 2 minute data for bright stars when binned to the same cadence (≈10%–25% better for T ≲ 8 mag, reaching equal precision at T ≈ 13 mag), consistent with pre-flight expectations based on differences in cosmic-ray mitigation algorithms. We present two results enabled by this improvement. First, we use 20 second data to detect oscillations in three solar analogs (γ Pav, ζ Tuc, and π Men) and use asteroseismology to measure their radii, masses, densities, and ages to ≈1%, ≈3%, ≈1%, and ≈20% respectively, including systematic errors. Combining our asteroseismic ages with chromospheric activity measurements, we find evidence that the spread in the activity–age relation is linked to stellar mass and thus the depth of the convection zone. Second, we combine 20 second data and published radial velocities to recharacterize π Men c, which is now the closest transiting exoplanet for which detailed asteroseismology of the host star is possible. We show that π Men c is located at the upper edge of the planet radius valley for its orbital period, confirming that it has likely retained a volatile atmosphere and that the “asteroseismic radius valley” remains devoid of planets. Our analysis favors a low eccentricity for π Men c (&lt;0.1 at 68% confidence), suggesting efficient tidal dissipation (Q/k 2,1 ≲ 2400) if it formed via high-eccentricity migration. Combined, these early results demonstrate the strong potential of TESS 20 second cadence data for stellar astrophysics and exoplanet science.

A Uniform Search for Nearby Planetary Companions to Hot Jupiters in TESS Data Reveals Hot Jupiters Are Still Lonely

Abstract We present the results of a uniform search for additional planets around all stars with confirmed hot Jupiters observed by the Transiting Exoplanet Survey Satellite (TESS) in its Cycle 1 survey of the southern ecliptic hemisphere. Our search comprises 184 total planetary systems with confirmed hot Jupiters with R p &gt; 8 R ⊕ and orbital period &lt;10 days. The Transit Least Squares algorithm was utilized to search for periodic signals that may have been missed by other planet search pipelines. While we recovered 169 of these confirmed hot Jupiters, our search yielded no new statistically validated planetary candidates in the parameter space searched (P &lt; 14 days). A lack of planet candidates nearby hot Jupiters in the TESS data supports results from previous transit searches of each individual system, now down to the photometric precision of TESS. This is consistent with expectations from a high-eccentricity migration formation scenario, but additional formation indicators are needed for definitive confirmation. We injected transit signals into the light curves of the hot Jupiter sample to probe the pipeline’s sensitivity to the target parameter space, finding a dependence proportional to R p 2.32 P − 0.88 for planets within 0.3 ≤ R p ≤ 4 R ⊕ and 1 ≤ P ≤ 14 days. A statistical analysis accounting for this sensitivity provides a median and 90% confidence interval of 7.3 − 7.3 + 15.2 % for the rate of hot Jupiters with nearby companions in this target parameter space. This study demonstrates how TESS uniquely enables comprehensive searches for nearby planetary companions to nearly all the known hot Jupiters.

Hot Jupiter formation in dense clusters: secular chaos in multiplanetary systems

ABSTRACT Exoplanetary observations reveal that the occurrence rate of hot Jupiters is correlated with star clustering. In star clusters, interactions between planetary systems and close flyby stars can significantly change the architecture of primordially coplanar, circular planetary systems. Flybys can impact hot Jupiter formation via activation of high-eccentricity excitation mechanisms such as the Zeipel–Lidov–Kozai (ZLK) effect and planet–planet scattering. Previous studies have shown that, for a two-planet system, close flybys, especially at high incidence angles, can efficiently activate the ZLK mechanism, thus triggering high-eccentricity tidal migration and ultimately form hot Jupiters. Here, we extend our previous study with a multiplanet (triple) system. We perform high-precision, high-accuracy few-body simulations of stellar flybys and subsequent planetary migration within the perturbed planetary systems using the code spacehub. Our simulations demonstrate that a single close flyby on a multiplanet system can activate secular chaos and ultimately lead to hot Jupiter formation via high-eccentricity migration. We find that the hot Jupiter formation rate per system increases with both the size of the planetary system and the mass of the outer planet, and we quantify the relative formation fractions for a range of parameters. Hot Jupiters formed via secular chaos are expected to be accompanied by massive companions with very long periods. Our study further shows that flyby-induced secular chaos is preferred in low-density clusters where multiplanet systems are more likely to survive, and that it contributes a significant fraction of hot Jupiter formation in star clusters compared to the flyby-induced ZLK mechanism.