{"id":"4c5b2aa7-8477-41a2-ab0e-8e7b19c7125e","arxiv_id":"2508.20891","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A new public code (SIRIO) predicts that reconnection-driven star-planet interaction can explain tentative radio detections from Proxima, YZ Ceti, and GJ 1151, while Alfvén-wing emission is too faint and free-free absorption is significant at high mass-loss rates.","lead":"This paper models radio emission from magnetic interactions between three close-in exoplanets and their M dwarf stars using a new public code. It finds the planets likely sit in the sub-Alfvénic regime and that reconnection-based mechanisms are more detectable than Alfvén wings, with stellar wind absorption mattering in YZ Ceti.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"PFSS source-surface radius is adopted, not derived; the sub-Alfvénic claims hinge on R_SS=4.5 R* and the paper's own sensitivity note shows nearby values change the conclusions qualitatively.","rationale":"The paper is best read as a conditional modeling claim: under a hybrid PFSS wind with R_SS = 4.5 R*, these systems are sub-Alfvénic and some SPI models predict detectable radio flux. The public SIRIO code, the successful benchmark against Turnpenney et al. (2018), and the addition of free-free absorption are genuine strengths. The central weakness is not the code internals but the adopted source-surface radius. The paper's own Sect. 3.1 sensitivity sentence and the Sect. 4.1 statement that Proxima b is super-Alfvénic for R_SS = 5.5 R* at B* = 200 G show that the qualitative conclusion is not robust to plausible changes in this parameter. A system-specific R_SS derived from pressure balance or MHD would settle whether 4.5 R* is physically appropriate for these stars. This is the same concern identified by the reader's weakest_assumption, so I agree with the CONDITIONAL verdict. I would not move to REJECT because the authors explicitly present the sub-Alfvénic conclusion as contingent on the hybrid PFSS geometry and the code is available for re-running; the appropriate scientific status is conditional on validating R_SS. The beta = 1e-3 efficiency uncertainty is also real but secondary: it scales flux predictions uniformly and does not affect the sub-Alfvénic vs super-Alfvénic regime.","tokens_in":27532,"tokens_out":7750,"duration_ms":80997,"concrete_test":"Compute R_SS for each star from the pressure-balance estimate of Réville et al. (2015): B(R_SS)^2/8π = ρ v_sw^2 + p_th, using the paper's own Parker-wind density/temperature profiles and the Table 1 values of B* and Mdot, including B* = 200 G for Proxima. Then re-run SIRIO's M_A and flux calculations with these self-consistent R_SS values. If the pressure-balanced R_SS is above the threshold where M_A = 1 at the planet's orbit, or if Proxima with B* = 200 G gives M_A > 1, the central sub-Alfvénic claim fails for that system; if the derived R_SS values are consistent with 4.5 R*, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is the choice R_SS = 4.5 R* in the hybrid PFSS geometry (Sect. 3.1). The central conclusion that Proxima b, YZ Cet b, and short-period GJ 1151 b are sub-Alfvénic and should produce SPI radio emission is obtained in the PFSS geometry; for a pure dipole, which is the limit of very large R_SS, the planets are super-Alfvénic (Sect. 4.1, Fig. 3). The paper adopts R_SS to match Pineda & Villadsen (2023), 'thus easily allowing comparisons,' rather than deriving it for each star, and immediately notes that 'lower (higher) values of R_SS would result in significantly higher (lower) values of the radio emission.' This is not a minor flux calibration: for Proxima with B* = 200 G the planet is always super-Alfvénic if R_SS = 5.5 R* (Sect. 4.1), and for GJ 1151 the PFSS sub-Alfvénic window closes for orbital periods above about 2 days (Sect. 4.3). Thus the qualitative claim that SPI is at work in these systems is contingent on R_SS lying below a threshold that is never computed from system-specific parameters. This concern is sharpest for Proxima and GJ 1151; YZ Cet is less sensitive because its Parker-spiral model is also sub-Alfvénic. The same fixed R_SS is nevertheless applied to all three stars despite different B*, Mdot, and temperatures, so the headline conclusion is carried by an unvalidated parameter choice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces SIRIO, a public Python code that models sub-Alfvénic magnetic star-planet interaction (SPI) radio emission using a 1D isothermal Parker wind with three magnetic geometries (pure dipole, Parker spiral, hybrid PFSS), three Poynting-flux prescriptions (Alfvén wing, reconnection, stretch-and-break), and free-free absorption. The authors benchmark SIRIO against the M-dwarf systems Proxima Centauri, YZ Ceti, and GJ 1151, comparing with previous modeling by Turnpenney et al. (2018), Kavanagh et al. (2021), Réville et al. (2024), and Pineda & Villadsen (2023). They find that, assuming a hybrid PFSS geometry with R_SS = 4.5 R*, Proxima b, YZ Cet b, and a short-period GJ 1151 b lie in the sub-Alfvénic regime, that the Alvén wing model predicts negligible fluxes, and that reconnection/stretch-and-break models predict detectable fluxes, with the stretch-and-break model matching the Proxima tentative detection and the GJ 1151 Stokes I flux to within a factor of two. Free-free absorption suppresses YZ Cet model A emission by ~75%.","tokens_in":27902,"tokens_out":4306,"duration_ms":44706,"significance":"If the central conclusions hold, the paper provides a computationally efficient, publicly available tool for target selection and for constraining stellar mass-loss rates and planetary magnetic fields from SPI radio searches. The paper's benchmarks reproduce several previously published results, and the inclusion of free-free absorption is a genuine addition that can affect detectability, as shown for YZ Cet. These are strengths that make the code and parameter exploration useful to the community. However, the headline conclusions are quantitatively contingent on at least two externally calibrated parameters (β and R_SS), and the qualitative sub-Alfvénic/super-Alfvénic boundary is not robust to the adopted R_SS, so the paper's central claim is currently conditional rather than established.","major_comments":[{"comment":"The choice R_SS = 4.5 R* is adopted to match Pineda & Villadsen (2023), not derived for each star. The paper's own sensitivity note states that lower/higher R_SS would significantly change the emission, and the text quantifies the fragility: for Proxima with B* = 200 G, Proxima b is always super-Alfvénic if R_SS = 5.5 R* (Sect. 4.1); for GJ 1151 the PFSS sub-Alfvénic window closes for orbital periods above ~2 days (Sect. 4.3); and a pure dipole (R_SS → ∞) gives super-Alfvénic conditions for all three systems. Since the central claim that SPI is at work in all three systems rests entirely on this single adopted value, the paper should either derive R_SS from system-specific pressure balance (e.g., the Réville et al. 2015 estimate cited in Sect. 3.1) or present the sub/super-Alfvénic regimes and predicted fluxes as a function of R_SS over its plausible range (≈2.5–5.5 R*). At present the q","section":"Sect. 3.1, 4.1, 4.3"},{"comment":"The conversion efficiency β is fixed at 10^-3, with the acknowledged range 10^-4 to 10^-2. Because F_R is linear in β, all quantitative comparisons to observations are conditional on this nominal value. Specifically, the 'factor of two' agreement with Vedantham et al. (2020) for GJ 1151 (Sect. 4.3) and the statement that the stretch-and-break model is 'the only one predicting values in agreement' with the Proxima detection (Sect. 4.1) would not hold at β = 10^-4 (fluxes an order of magnitude or more below the detections), and at β = 10^-2 the Alfvén-wing model would also become detectable. The paper should propagate this factor-100 uncertainty through the flux predictions and adjust the model-comparison claims accordingly, or state them explicitly as nominal-β results.","section":"Sect. 3.2, Eq. (9)"},{"comment":"The text notes that for Proxima in the PFSS geometry 'M_A is so close to unity that star-planet interaction could be very inefficient,' yet the flux predictions in Fig. 4 are computed with the full Poynting-flux formulas without any reduction for this marginal sub-Alfvénicity. This internal tension directly affects the Proxima detection claim: the predicted detectability assumes an efficiency that the same section casts into doubt. Please either quantify how the near-unity M_A modifies the Poynting flux/radio conversion, or explicitly caveat that the Proxima flux should be treated as an upper limit.","section":"Sect. 4.1"}],"minor_comments":[{"comment":"The summary mislabels the YZ Cet models: the high mass-loss, open-field model is called 'model B' and the low mass-loss PFSS model is also called 'model B'; from Sect. 4.2 the former should be model A and the latter model B.","section":"Summary, Sect. 5"},{"comment":"Typo: 'GJ 1551 b' and 'GJ 1551 c' should read 'GJ 1151 b' and 'GJ 1151 c'.","section":"Sect. 4.3"},{"comment":"Several minor typos: 'Fig, 5' (Sect. 4.1), 'with with' (Sect. 4.1), 'fells short' (Sect. 4.1), 'or and' (Sect. 4.1), and 'Potentially' capitalization in the abstract. These do not affect the science but should be cleaned up.","section":"Throughout"},{"comment":"The sentence 'We note that lower (higher) values of R_SS would result in significantly higher (lower) values of the radio emission' is not always true for all three models; for the reconnection and stretch-and-break models the dependence is more complex because R_mp and the magnetosphere condition also enter. Please qualify this statement or refer the reader to the relevant parameter study.","section":"Sect. 3.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful code paper with credible benchmarks, and the qualitative conclusions are interesting but over-stated relative to the adopted R_SS and β. The referee report asks for a sensitivity analysis on R_SS and β; if the authors can show that the sub-Alfvénic conclusions hold over a plausible system-specific range, the paper would be suitable for publication. I do not see grounds for rejection: the issues are fixable within the scope of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a genuinely useful paper, and the value is the tool more than the conclusions. SIRIO is public, easy to run, and the benchmarks hold up: they reproduce Turnpenney et al. profiles and Pineda & Villadsen's YZ Ceti architectures with the published parameters. The new physics is the free-free absorption treatment, which suppresses the YZ Ceti high-mass-loss model by ~75% and can kill the signal entirely at mass-loss rates above ~10 solar; that is a concrete, testable result. The GJ 1151 prediction — under a PFSS geometry the sub-Alfvénic window closes for periods above ~2 days — is also genuinely new and useful, and the stretch-and-break model's agreement with the tentative Proxima detection is shown as a comparison, not a fit.\n\nThe soft spots are real but are mostly on the table. The headline sub-Alfvénic claim rests on R_SS = 4.5 R*, adopted from Pineda & Villadsen rather than derived for each star. Under a pure dipole all three systems are super-Alfvénic, and for Proxima with the 200 G field, R_SS = 5.5 R* flips the planet into the super-Alfvénic regime. The paper states this sensitivity in a sentence but does not show it in the main flux figures. The abstract says \"assuming a hybrid PFSS geometry,\" so the conditional framing is there, but the title's \"iconic M dwarfs\" may overstate the security of the conclusions. Second, beta = 10^-3 is a nominal value spanning two orders of magnitude (10^-4 to 10^-2), and it scales the predicted fluxes linearly. The paper does not propagate that uncertainty into any of the quoted flux densities or the comparisons with detections. That is a real gap for a paper whose purpose is to guide observers.\n\nNone of this is disqualifying. The 1D isothermal wind is a simplification the authors acknowledge, and the code lets anyone re-run with different assumptions. I disagree a little with the stress-test's framing of R_SS as a load-bearing flaw: it is a sensitivity, clearly flagged, that the code itself lets you explore. For observational planning and for comparison of model families, this paper is a step up from models that fix a single geometry or ignore free-free absorption.\n\nWho benefits: anyone writing a proposal to look for SPI radio emission, and modelers who want a quick estimate before going to 3D MHD. I would bring it to the next reading group, and I would cite the code.\n\nRecommendation: send it to peer review. Ask for a sensitivity table showing how the fluxes and the sub-Alfvénic boundaries shift across the allowed beta range and across plausible R_SS values; that is a light revision and would make the conditional claims precise.","headline":"A useful public-code paper with honest conditional claims; the fluxes are order-of-magnitude only because beta is unconstrained and the R_SS choice is adopted, but the tool deserves peer review and a light revision.","tokens_in":28480,"tokens_out":3461,"would_cite":true,"duration_ms":36307,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Under a hybrid PFSS wind geometry, Proxima b, YZ Cet b, and the putative GJ 1151 b are sub-Alfvénic and should emit radio via star-planet interaction, with the stretch-and-break mechanism matching tentative detections.","keywords":["star-planet interaction","radio emission from exoplanets","M dwarf stars","Alfvén surface","stellar wind","cyclotron maser emission","free-free absorption","PFSS magnetic field"],"falsifier":"Measure each star's surface magnetic field topology and wind mass-loss rate (e.g., with spectropolarimetry and Lyman-alpha absorption) and locate the Alfvén surface. If Proxima b orbits beyond that surface for the measured ~600 G field, the sub-Alfvénic claim collapses. On the flux side, a sensitive low-frequency radio observation at the predicted cyclotron frequency and orbital phase that reaches an upper limit below the stretch-and-break prediction—roughly 0.1 mJy for Proxima—and sees nothing would falsify the model's match.","tokens_in":27378,"feed_emoji":"📡","tokens_out":9763,"duration_ms":87459,"temperature":0.7,"pith_summary":"This paper asks whether the tentative radio detections from three nearby M-dwarf systems—Proxima Centauri, YZ Ceti, and GJ 1151—are genuine magnetic star-planet interaction. It argues yes, provided the stellar wind field is a hybrid PFSS geometry with closed field lines out to 4.5 stellar radii and open beyond. Under that geometry all three planets sit in the sub-Alfvénic regime, so energy from the planet's motion can travel back to the star and generate cyclotron radio emission. The paper's code, SIRIO, predicts that the Alfvén wing mechanism is too faint, while reconnection and stretch-and-break mechanisms produce detectable fluxes, with stretch-and-break alone matching the Proxima detection and coming within a factor of two of the GJ 1151 detection. It also shows that free-free absorption in the stellar wind can suppress emission, explaining why some systems are seen and others not.","feed_headline":"Three M-dwarf planets may glow in radio from star-planet coupling","feed_subtitle":"A fast model crosses the Alfvén boundary for Proxima b, YZ Cet b and GJ 1151 b, matching tentative detections.","key_machinery":"SIRIO, a public Python code, couples a 1D isothermal Parker stellar wind (solved with the Lambert W function) to three magnetic field geometries: closed dipole, open Parker spiral, and hybrid PFSS, where field lines are closed dipolar out to R_SS = 4.5 R* and open radial beyond. It computes Poynting flux with three star-planet interaction models—Alfvén wing, magnetic reconnection, and stretch-and-break—whose relative sizes are set by the ratios S_rec/S_Alf = 2.5/M_A and S_rec/S_sb = 2.5 ξ^2/f_ap, with ξ = B_sw/B_pl and f_ap the connected-area fraction. The planet's magnetopause radius comes from pressure balance, and free-free absorption is integrated along the line of sight to the observer.","core_discovery":"The paper's central claim is that all three iconic M-dwarf systems are most likely sub-Alfvénic, meaning their close-in planets move through the stellar wind slower than the local Alfvén speed, so magnetic energy can flow back to the star and power electron-cyclotron radio emission. This holds only under a hybrid PFSS wind geometry with a source surface at 4.5 stellar radii; a pure dipole geometry puts all the planets in the super-Alfvénic regime and kills the interaction. Among the three emission mechanisms tested, the Alfvén wing model yields negligible flux, while reconnection and stretch-and-break predict detectable emission. The stretch-and-break model is the only one matching Proxima's","pith_inferences":["Measuring each star's actual magnetic source-surface radius would convert the model into a sharp test: radio detections or non-detections would then pin down the stellar wind mass-loss rate at the planet's orbit.","The planetary magnetic fields are estimated from a tidal-lock scaling law that gives lower limits, so stronger planetary fields would boost reconnection and stretch-and-break fluxes further, widening the detectability window.","If confirmed, SPI radio emission would provide a detection channel for close-in Earth-mass planets around M dwarfs that is independent of transits and radial velocities.","The same machinery could be applied to other M-dwarf systems with tentative detections to decide whether free-free absorption or magnetic geometry explains null results."],"forward_implications":["Under the hybrid PFSS geometry, all three systems are sub-Alfvénic, so SPI radio emission is physically expected and can be searched at the cyclotron frequency.","The stretch-and-break model sets the most promising flux scale: it alone matches Proxima's tentative detection and sits within a factor of two of GJ 1151's, so future observations should target its predicted frequency and orbital phase.","Free-free absorption can erase the signal at high wind mass-loss rates; for YZ Ceti model A about 75% is absorbed, so non-detections in high-wind stars do not rule out SPI.","GJ 1151 b, if real, must be low-mass (about 0.73 Earth masses) with an orbital period under two days for the PFSS geometry; the five-day high-mass scenario is ruled out.","SIRIO's fast 1D approach lets observers compute detectability curves as functions of orbital separation, mass-loss rate, and planetary magnetic field, guiding target selection and constraining stellar winds or planetary fields."],"supporting_citations":[{"why":"Supplies the Proxima Parker-spiral benchmark whose diagnostic plots SIRIO reproduces exactly.","marker":"Turnpenney et al. (2018)"},{"why":"Defines YZ Ceti Models A and B and the adopted R_SS = 4.5 R*, the baseline to which SIRIO is compared.","marker":"Pineda & Villadsen (2023)"},{"why":"Reports the GJ 1151 120-167 MHz detection whose flux SIRIO's stretch-and-break model matches within a factor of two.","marker":"Vedantham et al. (2020)"},{"why":"Reports Proxima's tentative SPI detection that only the stretch-and-break model reproduces.","marker":"Pérez-Torres et al. (2021)"},{"why":"Provides the 10^-3 Poynting-to-radio conversion efficiency and the Io-Jupiter interaction context.","marker":"Zarka (2007)"},{"why":"Supplies the Alfvén wing Poynting flux formula used in SIRIO.","marker":"Saur et al. (2013)"},{"why":"Supplies the reconnection model Poynting flux expression.","marker":"Lanza (2009)"},{"why":"Supplies the stretch-and-break Poynting flux expression SIRIO implements.","marker":"Strugarek et al. (2022)"},{"why":"Supplies the free-free absorption coefficient and Gaunt factor that SIRIO integrates along the line of sight.","marker":"Cox (2000)"},{"why":"Provides an alternative low-mass-loss Proxima scenario and R_SS estimate whose Alfvén surface SIRIO reproduces.","marker":"Réville et al. (2024)"}],"fun_headline_variants":["Sub-Alfvénic trio: Proxima, YZ Cet, GJ 1151 may shine in radio","M dwarf planets' magnetic dance may spark detectable radio","Radio whispers from Proxima b? New model says possible","Three M dwarfs may go radio: Proxima, YZ Cet, GJ 1151","Sub-Alfvénic planets: Proxima, YZ Cet, GJ 1151 could emit radio"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The load-bearing premise is that the magnetic field of each of the three stars opens up at 4.5 stellar radii (the hybrid PFSS geometry); if any field stays closer to a closed dipole, its planet becomes super-Alfvénic, no Poynting flux reaches the star, and the paper's radio predictions vanish.","fun_headline_variants_meta":{"raw":{"variants":["Sub-Alfvénic trio: Proxima, YZ Cet, GJ 1151 may shine in radio","M dwarf planets' magnetic dance may spark detectable radio","Radio whispers from Proxima b? New model says possible","Three M dwarfs may go radio: Proxima, YZ Cet, GJ 1151","Sub-Alfvénic planets: Proxima, YZ Cet, GJ 1151 could emit radio"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001,"raw_usage":{"total_tokens":4126,"prompt_tokens":859,"completion_tokens":3267,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":603,"completion_tokens_details":{"reasoning_tokens":3158}},"tokens_in":603,"tokens_out":3267,"duration_ms":23623,"temperature":1.0,"reasoning_tokens":3158,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T14:44:09.297385+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure each star's surface magnetic field topology and wind mass-loss rate (e.g., with spectropolarimetry and Lyman-alpha absorption) and locate the Alfvén surface. If Proxima b orbits beyond that surface for the measured ~600 G field, the sub-Alfvénic claim collapses. On the flux side, a sensitive low-frequency radio observation at the predicted cyclotron frequency and orbital phase that reaches an upper limit below the stretch-and-break prediction—roughly 0.1 mJy for Proxima—and sees nothing would falsify the model's match.","supporting_citations":[{"cited_title":"N., 2000, Allen's astrophysical quantities","cited_arxiv_id":null,"evidence_quote":"Supplies the free-free absorption coefficient and Gaunt factor that SIRIO integrates along the line of sight."}],"review_version":1}