Magnetic Star-Planet Interactions

• Magnetic star-planet interactions (MSPI) are electrodynamic and magnetohydrodynamic couplings where a planet in a sub‐Alfvénic regime induces Alfvén wings that transfer energy and momentum.
• MSPI drives observable phenomena such as aurorae, chromospheric modulation, radio bursts, and atmospheric escape, with interaction strength governed by magnetic topology and relative motion.
• Researchers use multi-wavelength observations and MHD simulations to infer planetary magnetic fields and assess implications for orbital migration in close-in exoplanetary systems.

Searching arXiv for recent and foundational papers on magnetic star–planet interactions to support the article. arxiv_search(query="magnetic star-planet interactions Alfvén wings exoplanets", max_results=10, sort_by="relevance") Magnetic star–planet interactions (MSPI) are electrodynamic and magnetohydrodynamic couplings between a planet and the magnetized plasma environment of its host star, especially strong for close-in systems in which the planet orbits inside the stellar sub-Alfvénic region. In that regime, the planet acts as a conducting obstacle embedded in the stellar wind or corona, launching Alfvénic perturbations, driving Poynting flux, and enabling angular-momentum exchange. The resulting phenomena include aurorae and shocks by analogy with the Solar System, and in exosystems have been associated with planetary heating or inflation, atmospheric escape, migration, and changes in the host star’s apparent activity (Strugarek, 2021). Across current theory and observation, MSPI is treated as a topology-dependent, time-variable process whose observables span chromospheric lines, photometric modulation, X-rays, ultraviolet tracers, and coherent radio bursts (Strugarek et al., 18 Feb 2025).

1. Physical regime and basic MHD description

The fundamental condition for strong magnetic coupling is sub-Alfvénic relative motion. For a planet moving through a magnetized stellar wind of density $\rho_{0}$, field strength $B_{0}$, and local bulk velocity $\mathbf v_{0}$ in the planet’s rest frame, the Alfvén speed is

$v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$

and the Alfvén Mach number is

$M_A=\frac{|\mathbf v_0|}{v_A}.$

When $M_A<1$, the disturbance remains causally connected along the magnetic field to the star, and a stationary interaction with two Alfvén wings can form (Strugarek, 2021).

The Alfvén wings lie in the plane spanned by the ambient field direction $\hat{\mathbf b}$ and the flow direction $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$. If $\Theta$ is the angle between $\mathbf v_0$ and $B_{0}$0, then each wing makes an angle $B_{0}$1 with the ambient field given by

$B_{0}$2

which reduces to $B_{0}$3 for the simple case $B_{0}$4 (Strugarek, 2021). This geometry underlies the Jupiter–Io analogy that recurs throughout the MSPI literature, including rocky-planet cases such as YZ Ceti b, where the oblique-wing interaction is explicitly treated as an Alfvén-wing interaction (Pineda et al., 26 Nov 2025).

Several dimensionless parameters organize the interaction. Beyond $B_{0}$5, the magnetic Reynolds number

$B_{0}$6

compares advection and magnetic diffusion, while the conductance ratio

$B_{0}$7

controls how effectively current closes through the obstacle (Strugarek, 2021). In the synthesized overview of the field, plasma $B_{0}$8,

$B_{0}$9

is also used to quantify the relative importance of thermal and magnetic pressure in the wind (Strugarek et al., 18 Feb 2025).

A central distinction is therefore between sub-Alfvénic and super-Alfvénic coupling. In the former, Alfvén wings form and energy and angular-momentum transport are focused; in the latter, an MHD bow shock appears and the coupling is weaker and confined to the wake (Strugarek, 2021). This regime structure is the starting point for interpreting both steady-state models and intermittent detections.

2. Energy transfer, Poynting flux, and magnetic torques

In ideal MHD, the electric field in the stellar-wind frame is

$\mathbf v_{0}$0

so the Poynting flux is

$\mathbf v_{0}$1

Projected along an Alfvén-wing characteristic $\mathbf v_{0}$2, the Alfvén-wing energy flux can be approximated as

$\mathbf v_{0}$3

Integrating over the obstacle cross-section and including the coupling efficiency $\mathbf v_{0}$4 yields the total Alfvén-wing power

$\mathbf v_{0}$5

with $\mathbf v_{0}$6 the effective radius of the planet-plus-magnetosphere obstacle (Strugarek, 2021).

A widely used simplified scaling takes $\mathbf v_{0}$7, $\mathbf v_{0}$8, $\mathbf v_{0}$9, and $v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$0, giving

$v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$1

Equivalent condensed summaries write the available power as

$v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$2

or, in unipolar-inductor form,

$v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$3

with the precise prefactors and angular factors depending on the adopted model (Strugarek et al., 18 Feb 2025, Figueira et al., 29 Jun 2026).

Part of this electromagnetic power is accompanied by angular-momentum transport. A simple estimate relates the magnetic torque to the power through

$v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$4

which leads to migration timescales

$v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$5

For sufficiently strong stellar fields and small orbital separations, $v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$6 can be as short as $v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$7–$v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$8 yr, whereas in the Solar System $v_A=\frac{B_0}{\sqrt{\mu_0\,\rho_0}},$9 yr (Strugarek, 2021). Self-consistent MHD summaries similarly state that magnetic torques can lead to orbital migration on $M_A=\frac{|\mathbf v_0|}{v_A}.$0 yr for hot Jupiters (Strugarek et al., 18 Feb 2025).

Alternative energetic channels exist. In the “stretch-and-break” picture, the planet’s dipole reconnects with the stellar field, is twisted by orbital motion, and intermittently releases energy; for HD 189733 this channel was found able to reach $M_A=\frac{|\mathbf v_0|}{v_A}.$1–$M_A=\frac{|\mathbf v_0|}{v_A}.$2 W for $M_A=\frac{|\mathbf v_0|}{v_A}.$3–$M_A=\frac{|\mathbf v_0|}{v_A}.$4 G, whereas the Alfvén-wing scenario was limited to $M_A=\frac{|\mathbf v_0|}{v_A}.$5–$M_A=\frac{|\mathbf v_0|}{v_A}.$6 W in the modeled states (Strugarek et al., 2022). This distinction is important because not all observed candidate powers are reproducible with the same interaction mechanism.

3. Magnetic topology as the controlling variable

A recurring result of global MHD calculations is that MSPI strength is controlled not only by field amplitude but by topology. In three-dimensional simulations of a magnetized planet in a self-consistent stellar wind, three representative configurations were studied: aligned, anti-aligned, and perpendicular (Strugarek et al., 2015). In the aligned case, the planetary dipole is locally parallel to the wind field; in the anti-aligned case it is anti-parallel; in the perpendicular case the fields are orthogonal.

These configurations produce markedly different magnetic obstacles. For a planet at $M_A=\frac{|\mathbf v_0|}{v_A}.$7, the aligned topology opens the magnetosphere on the day side and yields an obstacle up to $M_A=\frac{|\mathbf v_0|}{v_A}.$8 in cross-section, whereas the anti-aligned topology keeps the magnetosphere closed and limits the effective radius to $M_A=\frac{|\mathbf v_0|}{v_A}.$9 (Strugarek et al., 2015). The Poynting powers differ correspondingly: $M_A<1$0 W for the aligned case, $M_A<1$1 W for the anti-aligned case, and $M_A<1$2 W for the perpendicular case. Reversing the planetary field therefore reduces $M_A<1$3 by a factor $M_A<1$4 (Strugarek et al., 2015).

The same topological sensitivity appears in angular-momentum exchange. Normalized to the stellar-wind torque $M_A<1$5, the reported torques are $M_A<1$6 for aligned, $M_A<1$7 for anti-aligned, and $M_A<1$8 for perpendicular geometries, with migration timescales of $M_A<1$9 Myr, $\hat{\mathbf b}$0 Myr, and $\hat{\mathbf b}$1 Myr, respectively, in the adopted T Tauri-like setup (Strugarek et al., 2015). A plausible implication is that observational intermittency need not imply absence of coupling; it can arise from the planet sampling different local magnetic connectivities as the stellar field evolves.

This view is consistent with earlier axisymmetric and two-dimensional PLUTO calculations that contrasted unipolar and dipolar interactions. In that proof-of-concept study, the unipolar configuration produced a torque about $\hat{\mathbf b}$2 the steady-wind torque, while the dipolar interaction produced about $\hat{\mathbf b}$3 the wind-only torque (Strugarek et al., 2013). Later 2.5D simulation grids generalized this dependence, finding that for slowly rotating young stars with strong fields, magnetic torques can become comparable to or exceed tides, with migration timescales as short as $\hat{\mathbf b}$4 Myr at $\hat{\mathbf b}$5 for $\hat{\mathbf b}$6 G (Strugarek et al., 2014).

Topology also determines whether phase-locked emission is expected to be simple or confused. In force-free coronal models, chromospheric hot spots or flaring activity phased to the orbital motion are found only when the stellar field is axisymmetric; in non-axisymmetric fields, the modulation becomes multiperiodic and can be easily confused with intrinsic stellar variability (Lanza, 2012). This provides a theoretical basis for the long-standing “on/off” phenomenology of MSPI candidates.

4. Interaction morphologies and atmospheric coupling

MSPI is not limited to Alfvén wings and stellar hot spots. In three-dimensional MHD simulations including stellar and planetary outflows, the interaction morphology can be classified by the ordering of three characteristic length scales: the magnetospheric standoff radius $\hat{\mathbf b}$7, the wind interaction radius $\hat{\mathbf b}$8, and the tidal radius $\hat{\mathbf b}$9 (Matsakos et al., 2015). Four generic types are identified: a magnetically dominated bow-shock regime, a wind-dominated bow-shock and tail regime, an inspiraling-flow overflow regime, and a magnetically choked Roche-lobe overflow regime.

In Type I, $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$0, the stellar wind meets the closed planetary magnetosphere, producing a bow shock and a comet-like tail. In Type II, $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$1, a strong planetary wind opens the magnetosphere and forms a broader shock and tail. In Type III, $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$2, part of the shocked gas passes through $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$3, spirals inward, fragments through Kelvin–Helmholtz and Rayleigh–Taylor instabilities, and accretes onto localized hot spots well ahead of the sub-planetary point. In Type IV, $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$4, gas leaks from beyond the Hill radius but is guided by magnetic tension into a nearly radial accretion stream (Matsakos et al., 2015). These regimes connect the magnetic problem to transit asymmetries, ultraviolet absorption, and accretion-like signatures.

Planetary atmospheric escape is another major consequence. For close-in systems irradiated strongly enough to drive ionized outflows, the escaping plasma can carry additional current or mass-loading in the Alfvén wings, modify the local density $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$5, and thereby alter both $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$6 and the wing geometry (Strugarek, 2021). Lanza’s evaporation models quantify the atmospheric impact of magnetic energy deposition. Energetic electrons accelerated by reconnection can reach column densities of $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$7–$\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$8, comparable with or deeper than EUV photons, and can increase the mass-loss rate up to a factor of $\hat{\mathbf h}=\mathbf v_0/|\mathbf v_0|$9–$\Theta$0 in close-in $\Theta$1, massive $\Theta$2 Jupiter masses$\Theta$3 planets (Lanza, 2013).

The same work reports mass-loss rates up to $\Theta$4–$\Theta$5 for atmospheres heated by electrons accelerated at the planetary magnetosphere boundary, and average mass-loss rates up to $\Theta$6–$\Theta$7 in the case of magnetic loops interconnecting the planet with the star (Lanza, 2013). It further states that the star–planet magnetic interaction provides a source of energy generally comparable with or exceeding stellar EUV radiation for close-in planets. This makes MSPI relevant to atmospheric chemistry, escape, and ultraviolet transit signatures, not merely to stellar-side activity diagnostics.

5. Observational diagnostics across the spectrum

The observational literature treats MSPI as a multi-wavelength phenomenon. A methodological synthesis distinguishes radio diagnostics, chromospheric line monitoring, precision photometry, and transmission spectroscopy, and emphasizes that no single diagnostic suffices for unambiguous confirmation (Figueira et al., 29 Jun 2026). Instead, consistent phase-locked variability, recurrence when the field topology is favorable, and self-consistent energetics are required.

In radio, coherent Electron-Cyclotron Maser emission is the canonical signature. Its characteristic frequency is

$\Theta$8

for $\Theta$9 in gauss (Trigilio et al., 2023, Figueira et al., 29 Jun 2026). Nearly 100% circular polarization and a spectral cutoff tied to the local magnetic field are therefore key discriminants. Radio fluxes are commonly related to magnetic power by

$\mathbf v_0$0

with $\mathbf v_0$1 a conversion efficiency (Figueira et al., 29 Jun 2026).

Chromospheric searches most often use Ca II H&K and H$\mathbf v_0$2. Standard practice involves computing residual line-core fluxes relative to an epoch mean or baseline profile and then searching for orbital, synodic, or anti-synodic periodicities rather than merely stellar rotation (Figueira et al., 29 Jun 2026). Reported amplitudes in Ca II H&K are of order a few $\mathbf v_0$3 of the continuum, while H$\mathbf v_0$4 variations can reach $\mathbf v_0$5–$\mathbf v_0$6 in index during flares or plages (Figueira et al., 29 Jun 2026). The crucial claim in MSPI searches is not variability per se but orbit locking.

Precision photometry searches for a phase-locked orbital component in decompositions of the form

$\mathbf v_0$7

where $\mathbf v_0$8 captures variability at $\mathbf v_0$9 or the synodic period (Figueira et al., 29 Jun 2026). Expected amplitudes from magnetic hot spots or planet-driven faculae are $B_{0}$00–$B_{0}$01 ppm. Rolling periodograms, harmonic fits, pre-whitening of stellar rotation, and bootstrap false-alarm estimation are now standard tools in the field (Figueira et al., 29 Jun 2026).

Spectropolarimetry is additionally important because Zeeman–Doppler Imaging provides the large-scale stellar topology required to infer magnetic connectivity and favorable footpoint phasing (Strugarek et al., 18 Feb 2025). This has become particularly relevant in systems such as YZ Ceti, where the plausibility of SPI depends jointly on orbital phase and the stellar magnetic geometry (Pineda et al., 26 Nov 2025).

6. Representative systems, detections, and controversies

Selected case studies

System	Reported signature	Main constraint or ambiguity
HD 189733	Ca II K modulation at $B_{0}$02 d in 2013 August	Strong epoch dependence; other campaigns showed no clear SPI [(Cauley et al., 2018); (Fares et al., 2010)]
YZ Ceti	Polarized radio bursts recurring in orbital phase windows	Frequency, polarization sense, and field evolution remain debated (Trigilio et al., 2023, Pineda et al., 26 Nov 2025)
HD 118203	TESS variability at the planet’s orbital period in an eccentric system	Stellar rotation not fully excluded observationally (Castro-González et al., 2024)
Proxima Centauri	Phase-locked flare clustering for Proxima d; chromospheric periodicities near Proxima b and d	Magnetic-field estimate depends on geometry and flare interpretation (Osorio et al., 21 May 2026)

HD 189733 is a benchmark case because it combines chromospheric monitoring, spectropolarimetry, coronal extrapolation, and wind modeling. A uniform reanalysis of six Ca II K epochs found significant modulation only in August 2013, with a best-fit period $B_{0}$03 d, consistent with the orbital period, and a peak near $B_{0}$04, corresponding to a phase lead of about $B_{0}$05 ahead of the sub-planetary point (Cauley et al., 2018). The coronal field at the planet’s orbit in that epoch was reported as $B_{0}$06 mG, the largest among the studied epochs, strengthening the SPI interpretation because the released power scales with the stellar field (Cauley et al., 2018). Yet an earlier spectropolarimetric study reported no clear evidence of magnetospheric interactions in activity indicators, with rotation dominating the variability (Fares et al., 2010). Later 3D wind modeling reconciled some of this intermittency by showing that only the stretch-and-break mechanism could explain the observed $B_{0}$07 W Ca II K residual, and that the observational cadence implied a detection probability of only $B_{0}$08 to $B_{0}$09 (Strugarek et al., 2022).

YZ Ceti has become the principal rocky-planet radio case. uGMRT observations detected radio emission four times in nine epochs, with two detections showing $B_{0}$10–$B_{0}$11 circular polarization; when combined with earlier VLA detections, the phase clustering yields a $B_{0}$12 confidence against a random-flare origin (Trigilio et al., 2023). Modeling of the auroral radio emission inferred a stellar field of about $B_{0}$13 kG and a planetary polar field lower limit of $B_{0}$14 G (Trigilio et al., 2023). A later Zeeman–Doppler study measured the large-scale stellar field directly, obtaining $B_{0}$15 G, with $B_{0}$16 of magnetic energy in poloidal modes, $B_{0}$17 in the dipole, and $B_{0}$18 in axisymmetric modes; it concluded that the measured topology does not rule out SPI scenarios, but also identified tensions involving ECM frequencies, polarization sense, and the required $B_{0}$19 G (Pineda et al., 26 Nov 2025). This is an example of a genuine controversy rather than a settled detection.

HD 118203 provides a distinct eccentric hot-Jupiter case. TESS periodograms in four sectors showed a dominant peak at about $B_{0}$20 d with false-alarm probability $B_{0}$21, matching the orbital period of HD 118203 b, and phase-folded amplitudes of about $B_{0}$22–$B_{0}$23 ppt (Castro-González et al., 2024). No consistent stellar-rotation signal was found in ELODIE FWHM or ASAS-SN, and the evolved star’s projected rotation was argued to be incompatible with a $B_{0}$24 d spin period (Castro-González et al., 2024). The interpretation advanced is that eccentricity, through pseudo-synchronization, can maintain a larger planetary magnetic moment and thereby enhance detectability.

The Proxima system extends MSPI claims to terrestrial planets through optical high-resolution spectroscopy. In 117 ESPRESSO spectra, flare epochs identified via Fe I lines were found to cluster with significant statistical evidence at the orbital phase of Proxima d, while prewhitened chromospheric time series showed peaks near the orbital periods of Proxima b and d (Osorio et al., 21 May 2026). Modeling through helicity-driven reconnection and Poynting-flux formalism gave a likely polar magnetic field of $B_{0}$25 G for Proxima d, with a plausible range of $B_{0}$26–$B_{0}$27 G depending on geometry, radius, and flare intensity (Osorio et al., 21 May 2026). This suggests that MSPI may provide an indirect route to terrestrial exoplanet magnetometry, though the estimate remains model-dependent.

7. Inference of planetary magnetic fields and current limits

One of the principal motivations for MSPI research is that it offers one of the few indirect routes to planetary magnetic-field estimation. For hot Jupiters, orbitally modulated Ca II K emission has been used to infer surface magnetic fields by linking observed chromospheric power to a total dissipated power through an assumed fractional Ca II K radiative yield. In the flux-tube model,

$B_{0}$28

and, using an adopted Ca II K energy fraction of $B_{0}$29, inferred planetary surface fields for HD 189733 b, HD 179949 b, $B_{0}$30 Boo b, and $B_{0}$31 And b were reported in the range $B_{0}$32–$B_{0}$33 G (Cauley et al., 2019). The same work notes that these values are about $B_{0}$34–$B_{0}$35 times larger than classical rotation-based dynamo predictions for tidally locked hot Jupiters, but are consistent with scaling laws tied to internal heat flux (Cauley et al., 2019).

Radio-based inferences proceed differently. For YZ Ceti b, the observed cutoff near $B_{0}$36 MHz and the ECM relation implied a source-region field of about $B_{0}$37 G and, under a dipolar extrapolation, a stellar polar field of about $B_{0}$38 G; matching the radiated power to the incident power further required a planetary magnetosphere with $B_{0}$39–$B_{0}$40 G (Trigilio et al., 2023). For Proxima d, inversion of flare energetics under a helicity-driven reconnection model yielded the previously noted likely field of $B_{0}$41 G (Osorio et al., 21 May 2026). These estimates are explicitly model-contingent, since they depend on field geometry, emission efficiency, and obstacle size.

Current modeling emphasizes those uncertainties. For M-dwarf systems, Alfvén-wave-driven wind models indicate that several TRAPPIST-1 planets likely orbit within the Alfvén surface and that Proxima Cen b may lie at the edge of or just inside it, but the mass-loss rate, turbulent correlation length, and large-scale field strength remain uncertain enough to propagate one to two orders of magnitude uncertainty into predicted MSPI power and radio flux (Réville et al., 2024). The public SIRIO framework reaches similar conclusions for Proxima Centauri, YZ Ceti, and GJ 1151: the systems are likely sub-Alfvénic under a hybrid PFSS geometry, the pure Alfvén-wing model predicts very low radio emission, magnetic reconnection or stretch-and-break scenarios are more favorable for detection, and free-free absorption may be especially relevant in YZ Ceti (Peña-Moñino et al., 28 Aug 2025).

A common misconception is that any orbit-phased variability constitutes proof of MSPI. The observational syntheses argue the opposite: intrinsic stellar variability, changing large-scale topology, sparse cadence, and ambiguous periodicities can all mimic or erase the signal (Figueira et al., 29 Jun 2026). Another misconception is that a single mechanism explains all candidate systems. The HD 189733 wind-model study explicitly found that Alfvén wings could not deliver enough power for the reported 2013 chromospheric event, whereas stretch-and-break could (Strugarek et al., 2022). Conversely, other systems are framed primarily in Alfvén-wing terms. This suggests that “MSPI” is best understood as a family of related magnetically mediated couplings rather than a single canonical process.

The field’s near-term direction is correspondingly methodological: simultaneous spectropolarimetry and radio campaigns, denser phase-resolved chromospheric and photometric coverage, Zeeman-broadening constraints on total surface fields, and 3D MHD wind-plus-magnetosphere modeling are all identified as necessary for turning candidate signatures into firm characterizations (Pineda et al., 26 Nov 2025, Strugarek et al., 2022). As an overview of current results, MSPI is now established as a plausible and in some systems quantitatively constrained channel for energy, mass, and angular-momentum exchange in compact exoplanetary systems, but its empirical confirmation remains strongly dependent on topology, cadence, and multi-wavelength consistency.