A Search for Transit Duration Variations in M dwarf Multi-Planet Systems

The nominal habitable zone for exoplanets orbiting M dwarfs lies close to the host star, making dynamical considerations especially important. One consequence of this proximity is the expectation of spin synchronization, with implications for atmospheric circulation. Several mechanisms can maintain non-zero obliquities over long timescales in compact multi-planet systems, including capture into Cassini State 2 (CS2) and other forms of secular spin-orbit coupling; such pathways are plausible in the orbital architectures of close-in M-dwarf planets. In this study, we search for transit duration variations (TDVs) consistent with the nodal precession rates predicted by Laplace-Lagrange secular theory in compact M-dwarf multi-planet systems. Our sample includes 23 exoplanets orbiting 12 stars. We compare recent, high-precision transit durations obtained from JWST white-light curves with measurements published at the discovery epoch and afterward. The resulting transit duration variation ranges from seconds to minutes, and we fit a linear trend to duration versus time for each planet. All systems are consistent with flat (no TDV) at the 3σ level. The strongest candidate is TRAPPIST-1d, whose fitted slope differs from zero with 2.2σ confidence. We calculate the expected TDV signals predicted by secular precession and compare them to the observed limits. Our null detection is consistent with the low-impact-parameter regime, where theoretical TDVs are only a few seconds per decade and below our sensitivity. Higher-impact-parameter configurations predict substantially larger TDVs and are disfavored: under uniformly distributed geometries, at least half of the allowed configurations would be excluded.

💡 Research Summary

This paper presents the first systematic search for transit‑duration variations (TDVs) in compact multi‑planet systems orbiting M‑dwarf stars, using high‑precision white‑light curves from the James Webb Space Telescope (JWST). The authors motivate the study by noting that the habitable zones of M dwarfs lie very close to their host stars, making tidal forces strong enough to synchronize planetary spins. However, several mechanisms—most notably capture into Cassini State 2 (CS2) and other secular spin‑orbit resonances—can maintain non‑zero obliquities over gigayear timescales. If a planet’s orbital node precesses at the rate predicted by Laplace‑Lagrange secular theory, the line‑of‑sight impact parameter changes slowly, leading to a measurable drift in the transit duration.

The sample consists of 23 transiting planets in 12 M‑dwarf multi‑planet systems for which JWST/NIRSpec BOTS observations are publicly available. The time baseline between the discovery epoch (2014–2019) and the JWST observations spans 3–10 years, with an average of five years. The authors process the raw NIRSpec data with the JWST Science Calibration Pipeline (v1.17.1), correct for 1/f noise and bad pixels, and extract a single white‑light curve from the NRS1 detector to avoid inter‑detector offsets. After detrending with a second‑order polynomial baseline, they fit each light curve using the batman transit model within an MCMC framework (emcee), allowing the planet‑to‑star radius ratio, mid‑transit time, orbital period, scaled semi‑major axis, impact parameter, and two quadratic limb‑darkening coefficients to float. Eccentricities are fixed to zero, consistent with the low‑eccentricity nature of M‑dwarf multis.

Transit durations are computed from the fitted parameters via an analytic expression (Eq. 1) and compared to literature values from the discovery papers and subsequent follow‑up studies. For each planet the authors fit a linear trend of duration versus time, interpreting the slope as a potential TDV signal. All 23 planets are consistent with a flat (zero‑slope) model at the 3σ level. The most notable outlier is TRAPPIST‑1d, whose slope deviates from zero at 2.2σ—insufficient for a claim of detection but suggestive given its relatively high impact parameter.

The authors then calculate the expected TDV amplitudes from Laplace‑Lagrange secular theory. Nodal precession frequencies (g) for compact M‑dwarf multis are typically 10⁻⁴–10⁻² rad yr⁻¹ (≈0.005–0.5° yr⁻¹). In the low‑impact‑parameter regime (b ≲ 0.2) these frequencies translate into TDVs of only a few seconds per decade, well below the JWST precision of ~5–10 s. In contrast, higher‑impact‑parameter configurations (b ≈ 0.5–0.8) could produce TDVs of minutes per decade, which would be readily detectable. By assuming a uniform distribution of impact parameters, the authors find that at least half of the allowed high‑b configurations are ruled out by their non‑detections.

The discussion emphasizes that the null result is compatible with the majority of M‑dwarf planets having low impact parameters, as suggested by transit geometry studies. It also highlights the need for longer baselines (≥10 yr) and improved timing precision (≤1 s) to probe the modest TDVs expected for low‑b systems. Future JWST observations, as well as upcoming missions such as ARIEL or Twinkle, could target the subset of systems with higher impact parameters or provide the temporal coverage needed to detect the subtle signatures of secular nodal precession. Detecting such TDVs would offer a rare indirect probe of planetary obliquity states, potentially confirming the existence of long‑lived Cassini State 2 configurations and informing models of atmospheric dynamics on tidally influenced worlds.

In summary, the paper demonstrates that (1) current JWST white‑light data are insufficient to detect the small TDVs expected for most M‑dwarf multi‑planet systems, (2) the observed lack of TDVs is consistent with low‑impact‑parameter geometries, and (3) a substantial fraction of high‑impact‑parameter configurations can already be excluded. The work sets quantitative limits on secular nodal precession in these systems and outlines the observational requirements for future detection of TDVs, thereby advancing our ability to infer spin‑orbit dynamics and atmospheric circulation regimes on potentially habitable exoplanets.