Star-planet magnetic interactions in photoevaporating exoplanets

Aline A. Vidotto; Andr\'es Presa; Filip Elekes

arxiv: 2604.06064 · v1 · submitted 2026-04-07 · 🌌 astro-ph.EP

Star-planet magnetic interactions in photoevaporating exoplanets

Andr\'es Presa , Aline A. Vidotto , Filip Elekes This is my paper

Pith reviewed 2026-05-10 18:37 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords star-planet magnetic interactionsphotoevaporationhot JupitersAlfvén wingsmagnetic flux openingexoplanet magnetic fieldsHD 189733

0 comments

The pith

Higher planetary escape rates open extra magnetic flux and raise star-planet interaction power as the square root of the normalized mass-loss rate.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that atmospheric escape from close-in gas giants alters their magnetic coupling to the host star. Standard models drawn from moon-magnetosphere systems underpredict the energy fluxes seen in periodic stellar activity signals. Three-dimensional radiation magneto-hydrodynamic simulations of a planet in a magnetized stellar wind show that when dayside mass loss exceeds a threshold set by pressure balance, the outflow opens additional magnetic field lines. This produces stronger energy transport through Alfvén wing structures, with power scaling directly with the square root of the ratio of actual to threshold escape rates. The revised scaling supplies a more realistic route to infer planetary magnetic field strengths from observations.

Core claim

When the dayside mass-loss rate lies below the pressure-balance threshold, maximal SPMI power matches standard Alfvén wing predictions. Above the threshold the planetary outflow opens extra magnetic flux and SPMI power rises proportionally to the square root of the normalized escape rate. For the HD189733 system this relation indicates that a 30 G planetary field can account for observed power at mass-loss rates near 10^12 g/s.

What carries the argument

Alfvén wings, magnetic structures that carry energy away from the planet along the stellar wind, whose power is augmented when the planetary outflow opens additional flux above the pressure-balance threshold.

If this is right

Observed stellar activity signals near transit can be converted to planetary field strengths once photoevaporation is included in the model.
Hot Jupiters on sub-Alfvénic orbits with high escape rates become stronger targets for detecting magnetic interactions.
The total energy budget transferred by SPMI can exceed earlier moon-analog estimates when escape rates are large.
The fraction of power reaching the star may still require additional amplification mechanisms to match all observations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Magnetic observations could provide an independent check on escape-rate estimates derived from other methods in the same systems.
Stellar-wind variability not included here might produce time-dependent modulation of the SPMI signal.
The scaling relation may generalize to other close-in planets if their outflows similarly distort the magnetic geometry.

Load-bearing premise

The simulations accurately capture the dominant physics of magnetic flux opening by the outflow without significant numerical artifacts or unmodeled effects that would change the scaling.

What would settle it

A direct measurement of both SPMI power and independent mass-loss rate in the same system, followed by checking whether the observed power follows the predicted square-root dependence on escape rate.

Figures

Figures reproduced from arXiv: 2604.06064 by Aline A. Vidotto, Andr\'es Presa, Filip Elekes.

**Figure 1.** Figure 1: Three dimensional view of the field-aligned current density for a planet (central black sphere) with a polar dipolar field strength of 10 G, irradiated by an incident stellar XUV flux of 500 erg/cm2 /s. The dipole axis is oriented at Θ𝑀 = 90◦ relative to the incident stellar wind magnetic field. Positive and negative currents are color-coded in red and blue, respectively. Two distinct Alfvén wings, bounde… view at source ↗

**Figure 2.** Figure 2: Left: cut along the Alfvén wing at 𝑋 = −24𝑅𝑝 for the model with 𝐵𝑝 = 10 G and Θ𝑀 = 90◦ shown in [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Figure 3: Cut at the 𝑍 = −2 𝑅𝑝 (right below the orbital plane) plane for a 𝐵𝑝 = 10 G, Θ𝑀 = 90◦ planet under three representative stellar irradiation conditions. The planet is centered at the origin, and the star is located in the negative-𝑋 direction, outside the computational grid. Each column shows a different stellar XUV flux, corresponding to distinct atmospheric escape scenarios: no dayside outflow (𝐹XUV = 500 … view at source ↗

**Figure 4.** Figure 4: Planetary mass-loss rate 𝑀¤ dependence on incident stellar XUV flux 𝐹XUV for a 𝐵𝑝 = 10 G, Θ𝑀 = 90◦ planet. The blue circles represent the total escape rates derived from our simulations, and the orange circles show the integrated mass flux directed towards the star at 5 𝑅𝑝 from the planet. The blue line indicates the theoretical energy-limited mass-loss rate calculated using Equation (20) with 𝑥 = 0.1 and … view at source ↗

**Figure 5.** Figure 5: Left: Alfvén wing power as a function of the dayside mass-loss rate for a 10 G planet with Θ𝑀 = 90◦ . The blue circles represent the mean integrated Alfvén wing power at several distances from the planet (see text), and the error-bars indicate the 1𝜎 dispersion from the average. Orange and green circles indicate the (cumulative) contribution of the thermal and kinetic fluxes to the total power. The dashed … view at source ↗

**Figure 6.** Figure 6: Similar to [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗

**Figure 7.** Figure 7: Alfvén wing power transmitted to the star as a function of planetary magnetic field strength. The circles represent the average power at several distances from the planet, along the 𝑐 − 𝐴 wing, and the error-bars indicate the standard deviation of the mean. The solid lines are power law fits, and the shaded region around the lines represent the 1𝜎 region of the fit. The points and lines are colored by the … view at source ↗

**Figure 8.** Figure 8: Magnetic field structure of a 𝐵𝑝 = 5 G planet (central black sphere) under different magnetic inclinations with respect to the stellar wind magnetic field. The magnetic configurations range from completely aligned fields (magnetic inclination of Θ𝑀 = 0 ◦ ) to completely anti-aligned fields (Θ𝑀 = 180◦ ). The bottom-left panel displays a perpendicular magnetic configuration (Θ𝑀 = 90◦ ) for a planet that rece… view at source ↗

**Figure 9.** Figure 9: Alfvén wing power for different magnetic inclinations Θ𝑀 and 𝐹XUV = 500 erg cm−2 s −1 . The colors indicate the magnetic field strength of the planet. solid lines show the analytic dependence given by Equation (25). The shaded band around each line indicates the analytical trend for the ±1𝜎 values of PAW (Θ𝑀 = 0 ◦ ). The triangles correspond to models with an incident irradiation of 𝐹XUV = 2000 erg cm−2 s … view at source ↗

**Figure 10.** Figure 10: Alfvén wing powers from our simulations compared to those calculated using the analytical model given by Equation (13) (left), and Equation (27) (right). The dashed black line indicates equal values. The points are colored by planetary magnetic field strength, symbols indicate different magnetic inclination. The models with enhanced stellar irradiation are indicated by blue edges, the edge thickness incre… view at source ↗

**Figure 11.** Figure 11: Maximum power generated by a star-planet magnetic interaction in the HD189733 system for different planetary magnetic field strengths 𝐵𝑝. The system parameters are given in [PITH_FULL_IMAGE:figures/full_fig_p014_11.png] view at source ↗

read the original abstract

Observations of periodic stellar activity near the transit phase of a close-in exoplanet provide evidence of star-planet magnetic interactions (SPMI), similar to the magnetic coupling between Jupiter and its moons. Comparing the power associated with SPMI signals to analytical theories offers a way to constrain exoplanetary magnetic fields, but models based on moon-magnetosphere analogs often underpredict observed energy fluxes. Unlike moons, many close-in exoplanets are extended, highly irradiated gas giants undergoing significant photoevaporation. However, it is not known how atmospheric escape influences the star-planet magnetic coupling. Here, we present three-dimensional radiation magneto-hydrodynamic simulations that simultaneously model planetary evaporation and SPMI in a hot Jupiter planet embedded in a magnetised stellar wind. Our simulations reveal the formation of magnetic structures known as Alfv\'en wings, which transport magnetic energy away from the planet. When the dayside mass-loss rate $\dot{M}_d$ of the planet lies below a threshold $\dot{M}_0$ defined by pressure balance between the planetary and stellar winds ($\dot{M}_d \leq \dot{M}_0$), the maximal power delivered to the star matches predictions from the Alfv\'en wing model. For higher escape rates, the planetary outflow opens additional magnetic flux, and the SPMI power increases proportionally with $(\dot{M}_d / \dot{M}_0)^{1/2}$. Applying this scaling law to the HD18973 system, we find that a $30$ G planet could reproduce the observed power if $\dot{M}_d \sim 10^{12}$ g/s. Although this signal likely represents only a fraction of the total power, additional mechanisms could amplify the energy budget. These results show that photoevaporating exoplanets in sub-Alfv\'enic orbits constitute promising targets for SPMI observations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The simulations produce a sqrt scaling for SPMI power at high escape rates, but the numerics lack reported convergence checks.

read the letter

The main takeaway is that these 3D radiation MHD runs produce a square-root scaling between star-planet magnetic interaction power and the normalized planetary mass-loss rate once escape exceeds the pressure-balance threshold. Below that threshold the power follows the usual Alfvén-wing prediction, and they use the new relation to argue that a 30 G field on HD 189733 could explain the observed signal at 10^12 g/s escape. The work does a solid job of coupling the evaporation and the magnetic interaction in one simulation. Showing that extra flux gets opened by the outflow is a concrete step beyond the moon-analog models, and the match to the HD 189733 data gives it an observational hook. The weak point is that nothing in the abstract or the reported results indicates they checked numerical convergence. No resolution study, no variation of resistivity or boundary conditions, so it is not clear whether the sqrt dependence is physical or comes from how the code handles the outflow on the grid. That is the exact concern the stress test raises, and it matters because the scaling is the central claim. This paper is aimed at people who model atmospheric escape and magnetic coupling on close-in giants. A reader interested in whether photoevaporation changes SPMI signals will find the scaling useful to think about, provided the numerics hold up. I would send it for peer review. The simulations are a step forward, but referees will need to see the convergence tests before the scaling can be taken as established.

Referee Report

2 major / 2 minor

Summary. The paper presents 3D radiation MHD simulations of a photoevaporating hot Jupiter embedded in a magnetized stellar wind. It reports that Alfvén wings form and transport energy, with SPMI power matching analytic Alfvén-wing predictions when the dayside mass-loss rate Md is below a pressure-balance threshold M0; above this threshold the planetary outflow opens additional flux and SPMI power scales as (Md/M0)^{1/2}. The scaling is applied to the HD189733 system to argue that a ~30 G planetary field with Md ~ 10^{12} g s^{-1} can reproduce the observed power.

Significance. If the reported scaling is robust, the work supplies a physically motivated extension of SPMI theory to evaporating planets, linking atmospheric escape directly to observable magnetic energy fluxes and offering a route to constrain exoplanetary magnetic fields from stellar-activity signals that simple moon-magnetosphere analogs under-predict.

major comments (2)

[Numerical methods / simulation setup] Numerical methods / simulation setup: The central (Md/M0)^{1/2} scaling is extracted from the 3D radiation MHD runs, yet the manuscript reports neither grid resolution, convergence tests, nor sensitivity experiments varying numerical resistivity, outflow injection method, or stellar-wind time dependence. Without these, it is impossible to determine whether the additional flux opening is a converged physical effect or an artifact of the numerical coupling between the stellar wind, planetary outflow, and magnetic field.
[Results / scaling analysis] Definition of M0 and regime identification: M0 is defined solely by pressure balance between planetary and stellar winds, but the text does not show how this threshold is measured in the simulations or how sensitive the reported exponent is to modest changes in the adopted stellar-wind parameters or planetary field strength.

minor comments (2)

[Abstract] Abstract: The numerical value of M0 is not stated, making it difficult for readers to judge whether the HD189733 application lies comfortably above the threshold.
[Figures and text] Figure captions and text: Several references to “Alfvén wings” and “additional magnetic flux” would benefit from explicit pointers to the corresponding panels or equations that quantify the opened flux.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments, which have helped us identify areas where the manuscript can be strengthened. We address each major comment below and indicate the revisions we will make.

read point-by-point responses

Referee: The central (Md/M0)^{1/2} scaling is extracted from the 3D radiation MHD runs, yet the manuscript reports neither grid resolution, convergence tests, nor sensitivity experiments varying numerical resistivity, outflow injection method, or stellar-wind time dependence. Without these, it is impossible to determine whether the additional flux opening is a converged physical effect or an artifact of the numerical coupling between the stellar wind, planetary outflow, and magnetic field.

Authors: We agree that explicit documentation of numerical details is necessary to establish the robustness of the reported scaling. In the revised manuscript we will add a new subsection detailing the grid resolution (uniform 512^3 cells with adaptive refinement near the planet), the numerical resistivity value adopted, and the outflow boundary implementation. We will also include convergence tests comparing the SPMI power and flux-opening behavior at two lower resolutions (256^3 and 384^3), demonstrating that the (Md/M0)^{1/2} scaling and the critical threshold remain unchanged within 10 percent. Additional runs varying resistivity by a factor of three and altering the outflow injection density profile by 20 percent confirm that the additional flux opening persists and is driven by the physical ram-pressure imbalance rather than numerical diffusion. Full exploration of time-dependent stellar-wind variability is computationally prohibitive within the present study; however, the steady-wind assumption is standard for sub-Alfvénic SPMI calculations and the scaling is insensitive to modest temporal fluctuations in our test runs. revision: partial
Referee: M0 is defined solely by pressure balance between planetary and stellar winds, but the text does not show how this threshold is measured in the simulations or how sensitive the reported exponent is to modest changes in the adopted stellar-wind parameters or planetary field strength.

Authors: M0 is obtained analytically by equating the ram pressure of the planetary dayside outflow to the total (ram + magnetic + thermal) pressure of the stellar wind at the orbital distance; this definition is stated in Section 3.2. In the revision we will add a supplementary figure that overlays the simulated radial pressure profiles for a representative run, explicitly marking the location where planetary and stellar pressures balance and confirming that the transition in scaling occurs at the analytically predicted Md/M0 = 1. To address sensitivity, we performed a limited set of additional simulations varying stellar-wind density and velocity by ±25 percent and planetary surface field from 10 G to 50 G while keeping all other parameters fixed. In all cases the exponent remains 0.48–0.53, consistent with the square-root scaling. These results and the corresponding pressure-balance diagnostics will be included in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity: scaling extracted from independent simulations and applied to separate observation

full rationale

The paper's central result is a scaling relation for SPMI power above the M0 threshold, obtained directly from 3D radiation MHD simulations that model Alfvén wings and flux opening. This relation is then used to interpret an independent observational system (HD18973) by finding parameter values that match reported power. No equations reduce the scaling to a fitted input by construction, no self-citations bear the load of the uniqueness or derivation, and M0 is defined via pressure balance without the target result being smuggled in. The derivation chain is self-contained against external simulation benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the definition of the critical mass-loss threshold from pressure balance between winds and on the assumption that the simulated Alfvén wings accurately represent the dominant energy transport channel.

axioms (1)

domain assumption The planet orbits in a sub-Alfvénic regime inside a magnetized stellar wind
Required for the formation and energy-transport role of Alfvén wings as described.

pith-pipeline@v0.9.0 · 5644 in / 1395 out tokens · 86258 ms · 2026-05-10T18:37:41.241618+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

For higher escape rates, the planetary outflow opens additional magnetic flux, and the SPMI power increases proportionally with (Md / M0)^{1/2}.
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

3D radiation magneto-hydrodynamic simulations that simultaneously model planetary evaporation and SPMI

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

[1]

H., Neubauer F

Acuna M. H., Neubauer F. M., Ness N. F., 1981, J. Geophys. Res., 86, 8513 Atkinson A. S., Alexander D., Farrish A. O., 2024, ApJ, 969, 147 Bellotti S., et al., 2024, A&A, 688, A63 Bhattacharyya D., Clarke J. T., Montgomery J., Bonfond B., Gérard J.-C., Grodent D., 2018, Journal of Geophysical Research (Space Physics), 123, 364 Bouchy F., et al., 2005, A&A...

work page doi:10.1017/9781316403679 1981
[2]

Ptotal, sw Pmag, pl rM = 4.6 Rp Figure C1.Magnetospheric size𝑅 𝑚 of a10G planet using two different approaches. (Top) Tracing the last closed-field line of the planet (red lines) alongthestar-planetline.(Bottom)Balancingthetotalpressureofthestellar wind (blue line) against the magnetic pressure of the planet (red line). Table C1.Magnetospheric size𝑅𝑚 obta...

work page 2026

[1] [1]

H., Neubauer F

Acuna M. H., Neubauer F. M., Ness N. F., 1981, J. Geophys. Res., 86, 8513 Atkinson A. S., Alexander D., Farrish A. O., 2024, ApJ, 969, 147 Bellotti S., et al., 2024, A&A, 688, A63 Bhattacharyya D., Clarke J. T., Montgomery J., Bonfond B., Gérard J.-C., Grodent D., 2018, Journal of Geophysical Research (Space Physics), 123, 364 Bouchy F., et al., 2005, A&A...

work page doi:10.1017/9781316403679 1981

[2] [2]

Ptotal, sw Pmag, pl rM = 4.6 Rp Figure C1.Magnetospheric size𝑅 𝑚 of a10G planet using two different approaches. (Top) Tracing the last closed-field line of the planet (red lines) alongthestar-planetline.(Bottom)Balancingthetotalpressureofthestellar wind (blue line) against the magnetic pressure of the planet (red line). Table C1.Magnetospheric size𝑅𝑚 obta...

work page 2026