Scientific Premise

Star-planet magnetic interactions (SPI) have been theorized for two decades: a close-in planet orbiting inside its host star's sub-Alfvénic radius can channel energy along magnetic field lines connecting the planet to the stellar surface, triggering localized chromospheric hot spots and, in some cases, flares that are phase-locked to the planet's orbit. Recent work has confirmed this phenomenon in the young system HIP 67522 (Nature, 2025), where flares are statistically concentrated at specific orbital phases of the inner planet. Earlier, Cauley et al. (2019) inferred magnetic field strengths of 20–120 G for four hot Jupiters via Ca II K emission modulated at the planetary orbital period.

Despite these breakthroughs, no one has performed a systematic, blind survey of the full TESS short-cadence archive to identify phase-correlated flare excesses across the entire hot Jupiter population. The TESS dataset now spans 6+ years and covers most known hot Jupiter hosts at 2-minute or 20-second cadence. A population-level search could (a) dramatically expand the sample of SPI detections, (b) measure or constrain planetary magnetic field strengths as a function of planet mass, age, and orbital separation, and (c) test whether SPI scales with the sub-Alfvénic condition as theory predicts.

Target Datasets

• TESS 2-minute and 20-second cadence light curves (MAST/SPOC) for all confirmed hot Jupiter hosts (~500 systems with a < 0.1 AU).
• TESS Full-Frame Images (30 min / 10 min / 200 s cadence depending on sector) for fainter hosts not observed at short cadence.
• Ephemerides from the NASA Exoplanet Archive (orbital period, epoch, eccentricity).
• Stellar activity indicators: TESS light curve variability (rotation period, spot modulation amplitude), plus archival Ca II H&K and Hα indices from spectroscopic surveys (LAMOST, GALAH DR4, California Legacy Survey) where available.

Novelty

The HIP 67522 detection was a single-system study of an unusually young star. Cauley et al. (2019) used spectroscopic monitoring of four hand-picked systems. This project would be distinct because it: (1) is a blind, population-level flare-phase analysis of hundreds of hot Jupiter systems simultaneously; (2) uses photometric flare detection on TESS data (which is far more abundant than spectroscopic time series); (3) builds a hierarchical model that treats the SPI flare rate enhancement as a function of physical parameters (planet mass, semi-major axis, stellar B-field proxy, sub-Alfvénic radius), enabling the first empirical scaling relation for star-planet magnetic coupling.

Concrete Workflow

1. Compile a target list of all confirmed transiting hot Jupiters with TESS short-cadence data (expected ~300–400 systems with adequate coverage).
2. Apply an automated flare detection pipeline (e.g., AltaiPony or stella, both open-source and trained on TESS data) to every light curve. Record flare start times, durations, and energies.
3. Phase-fold flare occurrence times to the planetary orbital period for each system. For each system, compute the Rayleigh statistic (Z²) to test whether flare times are uniformly distributed in orbital phase or clustered.
4. At the population level, stack the phase-folded flare distributions across all systems (weighted by the number of flares per system). Test for a global excess at the expected SPI phase (near the sub-planetary longitude or shifted by the Parker spiral angle, depending on the Alfvén travel time).
5. Bin the sample by planet mass, orbital period, and stellar activity level. Fit a hierarchical Bayesian model: the SPI flare rate enhancement Δf is parameterized as Δf ∝ (Bp × Rp²) / (a² × v_sw), where Bp is the planetary magnetic field, Rp is planet radius, a is separation, and v_sw is stellar wind speed (proxied by stellar activity). This lets you solve for Bp (or its distribution) as the key latent variable.
6. Validate by checking that the signal disappears for a control sample of planets at a > 0.3 AU (where sub-Alfvénic interaction should be negligible).

Expected Signal / Observable

In HIP 67522, the planet-correlated flare rate was roughly 2–3× the baseline rate at specific orbital phases. For the broader hot Jupiter population, a more modest enhancement (~20–50% in phase-binned flare rate) would be expected, detectable at 3–5σ when stacking across ~100+ active host stars. Systems with young, active hosts and massive, close-in planets should show the strongest signal.

Possible False Positives

• Stellar rotation aliasing: If the stellar rotation period is close to the planetary orbital period (tidal synchronization), activity modulation mimics orbital-phase-locked flaring. Mitigation: measure the stellar rotation period from the light curve; exclude systems where P_rot/P_orb is within 5% of an integer ratio; test whether the flare-phase signal persists after removing the rotational modulation.
• Small-number statistics: Many systems will have few flares. Mitigation: the hierarchical model properly handles low-count systems by shrinking toward the population mean; the detection is made at the population level, not per-system.
• Flare detection false positives: Instrumental artifacts or cosmic rays flagged as flares. Mitigation: use injection-recovery tests to characterize the flare pipeline's false positive rate; require flares to span at least 3 consecutive cadences.

Why It Could Matter

Exoplanetary magnetic fields are among the most important yet poorly constrained properties in exoplanet science. Magnetic fields govern atmospheric escape rates (and thus long-term atmospheric evolution), control the interaction between planetary magnetospheres and stellar winds, and may be a prerequisite for surface habitability. A population-level measurement of hot Jupiter magnetic field strengths as a function of mass and age would be the first empirical constraint on the planetary magnetic dynamo scaling law beyond our solar system — directly relevant to understanding whether terrestrial planets in habitable zones retain protective magnetospheres.