---
title: Magnetic Star-Planet Interactions
url: https://www.emergentmind.com/topics/magnetic-star-planet-interactions-mspi
type: topic
---

# Magnetic Star-Planet Interactions

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 [2104.05968]. 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 [2502.13262].

## 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 [2104.05968].

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 $\mathbf B_0$, then each wing makes an angle $\Theta_A$ with the ambient field given by
$$
\sin\Theta_A
=
M_A\;\frac{\cos\Theta}{\sqrt{1 + M_A^2 - 2\,M_A\,\sin\Theta}},
$$
which reduces to $\sin\Theta_A=M_A$ for the simple case $\Theta=0$ [2104.05968]. 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 [2511.21853].

Several dimensionless parameters organize the interaction. Beyond $M_A$, the magnetic Reynolds number
$$
{\rm Re}_m=\frac{v_0L}{\eta}
$$
compares advection and magnetic diffusion, while the conductance ratio
$$
\bar\alpha=\frac{\Sigma_P}{\Sigma_P+\Sigma_A},\qquad \Sigma_A=\frac{1}{\mu_0 v_A},
$$
controls how effectively current closes through the obstacle [2104.05968]. In the synthesized overview of the field, plasma $\beta$,
$$
\beta=\frac{2\mu_0\,n\,k_B\,T}{B^2},
$$
is also used to quantify the relative importance of thermal and magnetic pressure in the wind [2502.13262].

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 [2104.05968]. 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 E=-\,\mathbf v\times\mathbf B,
$$
so the Poynting flux is
$$
\mathbf S=\frac{\mathbf E\times \mathbf B}{\mu_0}.
$$
Projected along an Alfvén-wing characteristic $\mathbf c_A^\pm=\mathbf v_0\pm\mathbf v_A$, the Alfvén-wing energy flux can be approximated as
$$
S_{AW}\approx \frac{v_0\,B_0^2}{\mu_0}\,\sin\Theta_A\cos\Theta.
$$
Integrating over the obstacle cross-section and including the coupling efficiency $\bar\alpha$ yields the total Alfvén-wing power
$$
P_{AW}\equiv\mathcal{P}
\simeq
2\pi\,R_\mathrm{eff}^2\,\bar\alpha^2\,
\sqrt{1 + M_A^2 - 2M_A\sin\Theta}\;\;S_{AW},
$$
with $R_\mathrm{eff}$ the effective radius of the planet-plus-magnetosphere obstacle [2104.05968].

A widely used simplified scaling takes $\Theta=0$, $\sin\Theta_A=M_A$, $\bar\alpha=1$, and $R_\mathrm{eff}\simeq R_p$, giving
$$
P_{AW}\sim \frac{\pi\,R_p^2}{\mu_0}\;B_*^2\;v_\mathrm{rel}\;M_A.
$$
Equivalent condensed summaries write the available power as
$$
\mathcal{P}_{\rm in}\sim \frac{B_w^2R_p^2v_{\rm rel}}{\mu_0},
$$
or, in unipolar-inductor form,
$$
P_{\rm mag}\simeq \frac{\pi R_p^2\,B^2\,v_{\rm rel}}{\mu_0},
$$
with the precise prefactors and angular factors depending on the adopted model [2502.13262; 2606.30743].

Part of this electromagnetic power is accompanied by angular-momentum transport. A simple estimate relates the magnetic torque to the power through
$$
\tau_{\rm mag}\approx\frac{P_{AW}}{\Omega_{\rm rel}}
=\frac{d\,P_{AW}}{v_{\rm rel}},
$$
which leads to migration timescales
$$
\tau_{\rm mig}=\frac{L_{\rm orb}}{\tau_{\rm mag}},
\qquad
L_{\rm orb}=M_p\sqrt{G\,M_*\,d}.
$$
For sufficiently strong stellar fields and small orbital separations, $\tau_{\rm mig}$ can be as short as $10^8$–$10^9$ yr, whereas in the Solar System $\tau_{\rm mig}\gg10^{12}$ yr [2104.05968]. Self-consistent MHD summaries similarly state that magnetic torques can lead to orbital migration on $\lesssim10^{8\ldots9}$ yr for hot Jupiters [2502.13262].

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 $10^{21}$–$10^{22}$ W for $B_p\sim1$–$30$ G, whereas the Alfvén-wing scenario was limited to $\sim10^{17}$–$10^{19}$ W in the modeled states [2203.10956]. 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 [1511.02837]. 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 $R_{\rm orb}=5\,R_\star$, the aligned topology opens the magnetosphere on the day side and yields an obstacle up to $\sim3\,R_P$ in cross-section, whereas the anti-aligned topology keeps the magnetosphere closed and limits the effective radius to $\sim1.2\,R_P$ [1511.02837]. The Poynting powers differ correspondingly: $\mathcal{P}\simeq1.4\times10^{19}$ W for the aligned case, $\mathcal{P}\simeq9.7\times10^{17}$ W for the anti-aligned case, and $\mathcal{P}\simeq7.7\times10^{17}$ W for the perpendicular case. Reversing the planetary field therefore reduces $\mathcal{P}$ by a factor $\sim14$ [1511.02837].

The same topological sensitivity appears in angular-momentum exchange. Normalized to the stellar-wind torque $\mathcal{T}_w$, the reported torques are $\mathcal{T}/\mathcal{T}_w\simeq0.11$ for aligned, $0.02$ for anti-aligned, and $0.13$ for perpendicular geometries, with migration timescales of $1.4\times10^2$ Myr, $8.5\times10^2$ Myr, and $2.7\times10^3$ Myr, respectively, in the adopted T Tauri-like setup [1511.02837]. 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 $4\times$ the steady-wind torque, while the dipolar interaction produced about $20\times$ the wind-only torque [1301.5239]. 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 $\sim20$ Myr at $r_{\rm orb}\approx2.5\,r_\star$ for $B_\star=246$ G [1409.5268].

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 [1206.5893]. 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 $R_m$, the wind interaction radius $R_w$, and the tidal radius $R_t$ [1503.03551]. 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, $R_t>R_m>R_w$, the stellar wind meets the closed planetary magnetosphere, producing a bow shock and a comet-like tail. In Type II, $R_t>R_w>R_m$, a strong planetary wind opens the magnetosphere and forms a broader shock and tail. In Type III, $R_w>\max(R_m,R_t)$, part of the shocked gas passes through $L_1$, 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, $R_m>\max(R_w,R_t)$, gas leaks from beyond the Hill radius but is guided by magnetic tension into a nearly radial accretion stream [1503.03551]. 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 $\rho_0$, and thereby alter both $v_A$ and the wing geometry [2104.05968]. Lanza’s evaporation models quantify the atmospheric impact of magnetic energy deposition. Energetic electrons accelerated by reconnection can reach column densities of $10^{23}$–$10^{25}\ {\rm m}^{-2}$, comparable with or deeper than EUV photons, and can increase the mass-loss rate up to a factor of $30$–$50$ in close-in $(<0.10\ {\rm AU})$, massive $(>1.5$ Jupiter masses$)$ planets [1307.2341].

The same work reports mass-loss rates up to $(0.5$–$1.0)\times10^{9}\ {\rm kg\,s^{-1}}$ for atmospheres heated by electrons accelerated at the planetary magnetosphere boundary, and average mass-loss rates up to $(0.3$–$1.0)\times10^{10}\ {\rm kg\,s^{-1}}$ in the case of magnetic loops interconnecting the planet with the star [1307.2341]. 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 [2606.30743]. 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
$$
\nu_c=\frac{eB}{2\pi m_e}\approx2.8\,B\ {\rm MHz}
$$
for $B$ in gauss [2305.00809; 2606.30743]. 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
$$
F_\nu\sim \eta\,\frac{P_{\rm mag}}{4\pi d^2\,\Delta\nu},
$$
with $\eta$ a conversion efficiency [2606.30743].

Chromospheric searches most often use Ca II H&K and H$\alpha$. 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 [2606.30743]. Reported amplitudes in Ca II H&K are of order a few $\times10^{-3}$ of the continuum, while H$\alpha$ variations can reach $\sim10^{-3}$–$10^{-2}$ in index during flares or plages [2606.30743]. 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
$$
y(t)=m_\star(t)+m_{\rm orb}(t)+m_{\rm tr}(t)+s(t)+\epsilon(t),
$$
where $m_{\rm orb}(t)$ captures variability at $P_{\rm orb}$ or the synodic period [2606.30743]. Expected amplitudes from magnetic hot spots or planet-driven faculae are $\sim10$–$10^3$ ppm. Rolling periodograms, harmonic fits, pre-whitening of stellar rotation, and bootstrap false-alarm estimation are now standard tools in the field [2606.30743].

Spectropolarimetry is additionally important because Zeeman–Doppler Imaging provides the large-scale stellar topology required to infer magnetic connectivity and favorable footpoint phasing [2502.13262]. 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 [2511.21853].

## 6. Representative systems, detections, and controversies

### Selected case studies

| System | Reported signature | Main constraint or ambiguity |
|---|---|---|
| HD 189733 | Ca II K modulation at $2.29\pm0.04$ d in 2013 August | Strong epoch dependence; other campaigns showed no clear SPI [1810.05253; 1003.6027] |
| YZ Ceti | Polarized radio bursts recurring in orbital phase windows | Frequency, polarization sense, and field evolution remain debated [2305.00809; 2511.21853] |
| HD 118203 | TESS variability at the planet’s orbital period in an eccentric system | Stellar rotation not fully excluded observationally [2401.17272] |
| 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 [2605.22925] |

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 $2.29\pm0.04$ d, consistent with the orbital period, and a peak near $\phi_{\rm orb}\approx0.9$, corresponding to a phase lead of about $40^\circ$ ahead of the sub-planetary point [1810.05253]. The coronal field at the planet’s orbit in that epoch was reported as $B_*^{\rm orb}\approx39$ mG, the largest among the studied epochs, strengthening the SPI interpretation because the released power scales with the stellar field [1810.05253]. Yet an earlier spectropolarimetric study reported no clear evidence of magnetospheric interactions in activity indicators, with rotation dominating the variability [1003.6027]. Later 3D wind modeling reconciled some of this intermittency by showing that only the stretch-and-break mechanism could explain the observed $\sim5\times10^{19}$ W Ca II K residual, and that the observational cadence implied a detection probability of only $12$ to $23\%$ [2203.10956].

YZ Ceti has become the principal rocky-planet radio case. uGMRT observations detected radio emission four times in nine epochs, with two detections showing $75$–$93\%$ circular polarization; when combined with earlier VLA detections, the phase clustering yields a $4.37\,\sigma$ confidence against a random-flare origin [2305.00809]. Modeling of the auroral radio emission inferred a stellar field of about $2.4$ kG and a planetary polar field lower limit of $0.4$ G [2305.00809]. A later Zeeman–Doppler study measured the large-scale stellar field directly, obtaining $\langle|B|\rangle_{\rm surface}\simeq225$ G, with $99\%$ of magnetic energy in poloidal modes, $71\%$ in the dipole, and $87\%$ 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_p\gtrsim10$ G [2511.21853]. 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 $6.1$ d with false-alarm probability $<0.1\%$, matching the orbital period of HD 118203 b, and phase-folded amplitudes of about $0.2$–$0.4$ ppt [2401.17272]. 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 $6.1$ d spin period [2401.17272]. 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 [2605.22925]. Modeling through helicity-driven reconnection and Poynting-flux formalism gave a likely polar magnetic field of $16$ G for Proxima d, with a plausible range of $3$–$280$ G depending on geometry, radius, and flare intensity [2605.22925]. 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,
$$
P_{\rm SPI}\simeq \frac{2\pi}{\mu}\,f_{AP}\,R_p^2\,B_{p0}^2\,v_{\rm rel},
$$
and, using an adopted Ca II K energy fraction of $0.21\%\pm0.08\%$, inferred planetary surface fields for HD 189733 b, HD 179949 b, $\tau$ Boo b, and $\upsilon$ And b were reported in the range $20$–$120$ G [1907.09068]. The same work notes that these values are about $10$–$100$ times larger than classical rotation-based dynamo predictions for tidally locked hot Jupiters, but are consistent with scaling laws tied to internal heat flux [1907.09068].

Radio-based inferences proceed differently. For YZ Ceti b, the observed cutoff near $\nu_{\min}\approx500$ MHz and the ECM relation implied a source-region field of about $90$ G and, under a dipolar extrapolation, a stellar polar field of about $2.4\times10^3$ G; matching the radiated power to the incident power further required a planetary magnetosphere with $B_{\rm planet}\gtrsim0.4$–$0.9$ G [2305.00809]. For Proxima d, inversion of flare energetics under a helicity-driven reconnection model yielded the previously noted likely field of $16$ G [2605.22925]. 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 [2410.01621]. 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 [2508.20891].

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 [2606.30743]. 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 [2203.10956]. 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 [2511.21853; 2203.10956]. As a synthesis 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.

Source: https://www.emergentmind.com/topics/magnetic-star-planet-interactions-mspi