{"id":"dd19ab83-7848-4685-ae57-83a3d317bb3d","arxiv_id":"2412.14072","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Ohmic heating from a time-varying stellar magnetic field can rival or exceed XUV heating in the upper atmospheres of close-in exoplanets and can screen the field above the 1 microbar level.","lead":"The paper calculates how much heat a time-changing magnetic field from a host star deposits in the upper atmospheres of hot exoplanets. It finds that for planets like Trappist-1 b and pi Men c, this Ohmic heating can rival or exceed X-ray heating and can block the field from reaching deeper layers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The heating rates assume a monochromatic orbital-period oscillation and omit the motional v×B induction of the stellar wind; the central numbers are not robust until advection is tested.","rationale":"The reader's weakest_assumption identifies the same load-bearing premise: the sinusoidal, orbital-period, non-advective 1D induction model. My stress-test agrees and sharpens it with a quantitative point: the magnetic Reynolds number in the peak-conductivity region is of order 10–500, so neglecting ∇×(u×B) is not a benign approximation there, and the frequency choice directly sets the headline heating rate through Eq. (14). This does not refute the paper's existence claim that Ohmic heating can reach 1e-3 erg/s/cm3 in favorable cases; if anything, including motional induction could add more heating. But it does mean the specific profiles, the Q/Qph ratios at the sonic point, and the inferred BP ranges in Figs. 5–6 are not self-consistent predictions of the full electromagnetic interaction. Since the paper explicitly frames the model as a first simplified step and discusses the advective term as future work, the conditional verdict remains appropriate. I see no need to move to reject or accept: the mechanism is plausible and the math is internally consistent, but the central quantitative claims are contingent on an assumption that a concrete numerical experiment can settle.","tokens_in":26387,"tokens_out":18616,"duration_ms":185379,"concrete_test":"Implement a 1D time-dependent induction code solving ∂B/∂t = ∇×(u×B) + ∇×(η∇×B) for the García Muñoz (2023) atmosphere profiles of Trappist-1 b and π Men c, with u taken from the outflow solution and with the stellar-wind field advected at the relative wind speed (orbital plus wind velocity) imposed at the top. Compare the time-averaged Ohmic heating profile Q(z) and the Q/Qph ratio at the sonic point with the paper's Figs. 5–6. If the advective contribution changes Q(rc) by more than a factor of about 3 or shifts the intermediate-BP range, then the central quantitative claim is conditional on a term the model omits; if the change is small, the simplified formalism is adequate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on Eq. (14), Qmax = Ω Bsw^2/(8π), which presumes that the only source of currents is a monochromatic ∂B/∂t at Ω = 2π/Porb and that the atmosphere is a stationary conductor. In the planet's rest frame, the stellar wind is a flowing magnetized plasma: the electromagnetic driver is the motional electric field E = -(v_orb + v_wind)×B/c, not a pre-specified sinusoidal B(t). The full induction equation is ∂B/∂t = ∇×(u×B) + η∇²B, and the paper's Eq. (10) drops the ∇×(u×B) term from the outset. For a rotating tilted stellar dipole, the field at the planet varies at P_rot or a beat period rather than Porb; for advected wind structures, the characteristic frequency is v/L. Because Qmax ∝ Ω, the assumed frequency directly sets the claimed 1e-3 erg/s/cm3. Moreover, the motional E persists even when ∂B/∂t = 0, so the formalism omits a potentially dominant Ohmic channel. The Sect. 6 justification that a radial outflow makes advection negligible is not quantitatively supported: with σP ≈ 0.4 S/m (π Men c), u ≈ 10 km/s near the sonic point, and L ≈ 100 km, the magnetic Reynolds number is Rm = uL/η ≈ 500, so advection and diffusion are comparable in exactly the high-conductivity region where the heating peaks.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a 1D model for the penetration of an external time-varying magnetic field into the upper atmospheres of hot exoplanets and for the associated Ohmic heating. It solves a diffusion equation for the vector potential with a Pedersen conductivity (Eq. 12), validates the solver against an analytic constant-conductivity solution and against Earth ionospheric conductivities, and applies the framework to Trappist-1 b and π Men c using photochemical upper-atmosphere models. It introduces a saturation heating rate Qmax = ΩBsw^2/(8π) (Eq. 14), compares Ohmic heating with XUV heating at the sonic point, and produces a population-level map of Qmax. The central claim is that Ohmic heating can reach ~10^-3 erg s^-1 cm^-3 and can surpass XUV heating for intermediate planetary field strengths of roughly 0.01-1 G, making it a potentially important term in the thermal budget of escaping atmospheres of close-in planets.","tokens_in":26682,"tokens_out":18225,"duration_ms":160188,"significance":"If the central claim holds, the paper establishes a new heating mechanism that must be included in upper-atmosphere energy balance models of close-in exoplanets. The manuscript has several clear strengths: the diffusion problem is solved with a numerically validated scheme; the conductivity calculations are benchmarked against Earth ionospheric data (Appendix B) and against the analytic constant-conductivity solution (Appendix C); the paper explicitly reconciles its heating-rate saturation with the earlier Cohen et al. (2024) result (Appendix D); and the population-level Qmax map is a simple, falsifiable prediction. These strengths make the paper a useful contribution to the star-planet interaction literature. However, the physical driver of the time variation (the assumed monochromatic orbital-period oscillation) and the neglect of advective induction are load-bearing assumptions that need to be justified or relaxed before the quantitative claims can be considered robust.","major_comments":[{"comment":"The induction equation used in the paper, ∂xx A0 + (4πiΩσP/c^2) A0 = 0, omits the advective term ∇×(u×B). In the planet's rest frame, the stellar wind generates a motional electric field E = -u×B/c that is present even when ∂B/∂t = 0. The statement in Sect. 6 that a radial outflow makes this term negligible is not quantitatively supported. Using the paper's own numbers for π Men c (σP≈0.4 S/m, u≈10 km/s near the sonic point, L≈100 km), the magnetic Reynolds number is Rm≈500, so advection dominates over diffusion precisely in the region where the heating peaks. Consequently, the computed heating rates and the 'upper limit' Qmax are not robust for a planet embedded in a flowing stellar wind. This limitation is acknowledged in Sect. 7 ('we have neglected so far inductive effects'), but it affects the central quantitative claim of the paper and should be addressed with a quantitative ordering argument or by including the advective term.","section":"Eq. (10), Sect. 6"},{"comment":"The population estimate and the claimed maximum heating rate assume a monochromatic oscillation at Ω = 2π/Porb. For a planet orbiting a rotating tilted stellar dipole, the magnetic-field component in the planetary frame varies at the stellar rotation period or at a beat period, and for advected wind structures the characteristic timescale is L/v_wind rather than the orbital period. Since Qmax ∝ Ω, the frequency choice directly sets the peak values in Fig. 7. The paper uses rotation periods of 3 and 30 days only in the scaling law for B⋆ (Eq. 20) and not for Ω, so the population-level heating rates can change by factors of several depending on which physical timescale actually modulates the field. The authors should specify the driver of the time variation and justify the use of the orbital frequency in the population analysis.","section":"Eq. (14), Sect. 5.2"},{"comment":"The statement in Sect. 5.2 that 'this Qmax can never be surpassed' is contradicted by the paper's own Fig. 3c, where a localized conductivity enhancement (green curve) produces a peak Ohmic heating rate eight times larger than the corresponding constant-conductivity case, which represents Qmax for those parameters. Appendix C provides the δP >> d limit Q = (σP/σP0)^3 Qmax, and while that particular limit lies below Qmax, the non-uniform-conductivity case in Fig. 3 shows that local volumetric heating can exceed Qmax. Therefore the population estimate in Fig. 7, which uses Eq. (14) as a nominal upper bound, is not a rigorous upper limit for a stratified atmosphere; it is the value for a uniform conductor with skin depth smaller than the layer depth. The text should be corrected to state the conditions under which Qmax applies.","section":"Sect. 5.2 and Fig. 3c"}],"minor_comments":[{"comment":"The formula for the electron-ion collision frequency is dimensionally unclear as printed: the text states n_i is in g cm^-3, but a collision frequency should scale with number density; please clarify the units or the intended formula.","section":"Eq. (19)"},{"comment":"The label 'Maedea 1977' is a typo for 'Maeda 1977'.","section":"Fig. B.1"},{"comment":"The choice of Prot = 3 and 30 days for all stars in the population estimate is not clearly justified for the actual stellar sample; the magnetic field scaling of Eq. (20) is sensitive to this parameter, and the paper does not discuss the resulting uncertainty in the Qmax map.","section":"Sect. 5.2"},{"comment":"The statement that a CME 'sweep[s] an Earth-like planet in a few seconds to a few minutes' is imprecise: for v=3000 km/s and a planet radius of ~6400 km, the crossing time is about 4 seconds; a few minutes would require a much larger structure or a slower speed.","section":"Sect. 6"}],"recommendation":"major_revision","confidential_remarks":"The advection caveat is likely to be the main point of contention for space-physics readers; even though the authors acknowledge it, the quantitative magnitude of the effect (Rm≈500 in the peak-heating region) means that the central numbers may be materially different when the motional term is included. The internally contradictory statement about Qmax being an absolute upper limit, contradicted by the paper's own Fig. 3, is an important issue to fix in revision. The paper is within the scope of A&A and the model development is a useful contribution, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look. The paper gives a clean, self-consistent 1D model for how a time-varying external magnetic field penetrates a hot exoplanet's upper atmosphere and deposits heat, using photochemical conductivity profiles from García Muñoz's models for Trappist-1 b and π Men c. The saturation bound Qmax = ΩBsw²/8π is derived clearly (and matches Chyba & Hand 2021 in a different geometry), and the comparison to Cohen et al. 2024 explains why earlier estimates were too high by 2-3 orders of magnitude: they ignored skin-depth screening. The validation against the analytic constant-conductivity solution and against Earth ionospheric conductivities is honest and useful. If the mechanism operates as described, Ohmic heating belongs in the thermal budget of close-in rocky planets around M dwarfs, and screening of external fields by the upper atmosphere is a real effect.\n\nThe soft spots are the usual ones for a first application. The external field is treated as a monochromatic sinusoid at the orbital period. That is a convenient idealization, but the real driver is the planet moving through a stellar wind and a non-axisymmetric stellar field; the characteristic frequency could be the stellar rotation period, a beat period, or v/L for advected wind structures. Since Qmax scales linearly with Ω, the frequency choice matters. The paper acknowledges the neglect of the ∇×(u×B) advection term in Sect. 6, but does not quantify it; the stress-test's estimate (Rm on the order of 10²-10³ near the sonic point for π Men c) suggests it is not obviously negligible and should have been addressed explicitly. The conclusions for the two planets also depend on BP and Bsw, both poorly constrained; the parameter maps are honest about this, but the stated ranges (e.g., BP∈[0.03,0.1] G for π Men c) are not robust yet. The atmosphere models themselves carry no error bars.\n\nNone of this breaks the central claim. The paper is a legitimate first step, and the formalism is reusable. The right next step is a version that includes advection and a more realistic driver spectrum. I'd send it to a knowledgeable referee; the issues are addressable, and the mechanism deserves to be tested rather than dismissed.","headline":"A careful 1D treatment that makes a plausible case for Ohmic heating as a thermal term in hot exoplanet upper atmospheres, but the central numbers rest on an idealized sinusoidal driver and unconstrained field strengths.","tokens_in":27243,"tokens_out":4725,"would_cite":true,"duration_ms":43580,"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":"Ohmic heating from an orbital-timescale, time-varying stellar-wind magnetic field can reach $10^{-3}$ erg s$^{-1}$ cm$^{-3}$ in hot-exoplanet upper atmospheres and can rival or exceed XUV heating for intermediate planetary fields.","keywords":["Ohmic heating","hot exoplanets","upper atmosphere","time-varying magnetic field","magnetic screening","Pedersen conductivity","Trappist-1 b","π Men c"],"falsifier":"Run a time-dependent 3D stellar-wind model for Trappist-1 and decompose the magnetic-field fluctuations at the planet's orbit in Fourier space: if the variance is dominated by timescales much shorter than the 1.51-day orbital period (CME crossings or advection past the planet), then the sinusoidal, orbital-frequency assumption at the heart of the skin-depth calculation fails, and the predicted penetration and heating profiles do not apply.","tokens_in":26192,"feed_emoji":"🪐","tokens_out":9410,"duration_ms":79693,"temperature":0.7,"pith_summary":"Close-in exoplanets are heated by more than the XUV radiation of their host star. This paper argues that the magnetic field of the stellar wind, varying as the planet orbits, induces currents in the upper atmosphere whose Ohmic dissipation can deposit up to $10^{-3}$ erg s$^{-1}$ cm$^{-3}$ locally, comparable to or larger than photochemical XUV heating in favorable cases. That matters because models of thermally escaping atmospheres that include only XUV heating would omit this term. The paper shows the effect is strongest for close-in planets around low-mass, fast-rotating stars, and that strong Ohmic heating also implies efficient screening of the external field by the upper atmosphere.","feed_headline":"Stellar magnetic swings can heat exoplanet air as much as XUV","feed_subtitle":"Induced currents can deposit up to 10^-3 erg/s/cm3 in hot-planet skies, rivaling starlight heating.","key_machinery":"The machinery is the 1D induction equation for the vector potential amplitude $A_0(x)$, $\\partial_{xx} A_0 + (4\\pi i \\Omega \\sigma_P/c^2) A_0 = 0$, with the time-varying field imposed at the top and zero current at the bottom. The penetration is set by the skin depth $\\delta_P = c/\\sqrt{2\\pi \\Omega \\sigma_P}$. The volumetric heating is $Q = [\\sigma_P/(\\sigma_H^2+\\sigma_P^2)] (c/(4\\pi) \\Re(\\partial_{xx} A_0))^2$, and in the thin-skin limit it saturates at $Q_{\\max} = \\Omega B_{\\mathrm{sw}}^2/(8\\pi)$, independent of the Pedersen conductivity. A key result is that the heating peak need not sit at the conductivity maximum, and that localized conductivity enhancements pull both the screening and the heat deposition toward them.","core_discovery":"The paper's central claim is that the upper atmospheres of hot exoplanets ohmically dissipate a time-varying external magnetic field, and that this dissipation can be a leading term in the local energy budget. Solving the 1D induction equation with ab-initio conductivity profiles, the authors find maximum volumetric heating rates up to $10^{-3}$ erg s$^{-1}$ cm$^{-3}$, with a theoretical cap $Q_{\\max} = \\Omega B_{\\mathrm{sw}}^2/(8\\pi)$ that is independent of the Pedersen conductivity once the skin depth is small. For intermediate steady planetary fields, roughly 0.01 to 1 G, the heating can peak near or below the sonic point of the escaping atmosphere; for very small fields the atmosphere screens the external field high up, and for very large fields the field passes through without being dissipated. Applied to Trappist-1 b and $\\pi$ Men c, both planets are predicted to experience substantial Ohmic heating.","pith_inferences":["Editorial inference: if the $Q_{\\max}$ saturation is general, then searches for induced radio emission or magnetic star-planet interaction signatures around close-in planets can use $Q_{\\max} = \\Omega B_{\\mathrm{sw}}^2/(8\\pi)$ as an upper bound on the energy available to any Ohmic dissipation channel.","Editorial inference: the screening result implies that interior-induction heating scenarios for close-in exoplanets must first pass through the upper atmosphere; planets with strong XUV-driven electron densities may block the varying field before it reaches the interior, weakening interior Ohmic heating in exactly the systems where it was previously predicted to be strongest.","Editorial inference: atmospheric composition becomes a first-order control on this heating channel—an H2O-dominated upper atmosphere can suppress electron density and therefore conductivity relative to an H-dominated one at comparable XUV flux, so escape and habitability models for M-dwarf planets may need composition-dependent Ohmic heating terms."],"forward_implications":["Ohmic heating from an orbital-timescale, time-varying external field can reach $10^{-3}$ erg s$^{-1}$ cm$^{-3}$ in hot-exoplanet upper atmospheres and can locally exceed XUV photochemical heating.","For planetary magnetic fields in the roughly 0.01–1 G range, the heating peak can sit close to or below the sonic point of the escaping atmosphere, so it should enter the thermal budget that sets atmospheric mass loss.","When Ohmic heating is strong, the upper atmosphere screens the external time-varying field almost completely, so little of it reaches deeper atmospheric layers or the planetary interior.","Trappist-1 b and $\\pi$ Men c are both predicted to experience substantial Ohmic heating, with the effect becoming important for ambient fields of about 1 G and 0.1 G, respectively.","For the known exoplanet population, the maximal heating rate is largest for close-in planets around fast-rotating, low-mass stars."],"supporting_citations":[{"why":"Supplies the induction-equation framework and the analytical skin-depth solution this paper generalizes to arbitrary conductivity profiles.","marker":"Parkinson (1983)"},{"why":"Provides the ab-initio photochemical models of the Trappist-1 b and π Men c upper atmospheres from which electron densities and conductivities are computed.","marker":"García Muñoz (2023)"},{"why":"Supplies the parametrized collision frequencies used to build the parallel, Pedersen, and Hall conductivities.","marker":"Schunk & Nagy (1980)"},{"why":"Gives the conductivity and ionization framework applied to exoplanet atmospheres.","marker":"Johnstone et al. (2018)"},{"why":"Established induction heating of exoplanet interiors from time-varying fields, the baseline the upper-atmosphere screening result qualifies.","marker":"Kislyakova et al. (2017)"},{"why":"Prior upper-atmosphere Ohmic-heating model whose heating rates are compared and bounded by the Qmax saturation.","marker":"Cohen et al. (2024)"},{"why":"Provides the stellar magnetic-field scaling used to estimate Qmax across the exoplanet population.","marker":"Ahuir et al. (2020)"},{"why":"Dynamo scaling-law estimates of planetary magnetic moments that define the intermediate BP range of interest.","marker":"McIntyre et al. (2019)"},{"why":"Observational estimate of Trappist-1's surface magnetic field used to infer plausible Bsw at the planet.","marker":"Reiners & Basri (2010)"}],"fun_headline_variants":["Magnetic swings rival XUV heating in hot exoplanet skies","Time-varying fields heat exoplanet air as much as XUV","Stellar magnetism can match XUV in heating hot exoplanets","Induced currents heat exoplanet atmospheres near starlight level"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything hangs on treating the magnetic field felt by the planet as a single smooth oscillation repeated once per orbit, with no motion of the wind itself; if the real field varies in sharper bursts or gets swept past the planet, the calculated heating and screening profiles change.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic swings rival XUV heating in hot exoplanet skies","Time-varying fields heat exoplanet air as much as XUV","Stellar magnetism can match XUV in heating hot exoplanets","Induced currents heat exoplanet atmospheres near starlight level"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000184,"raw_usage":{"total_tokens":1405,"prompt_tokens":1117,"completion_tokens":288,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":733,"completion_tokens_details":{"reasoning_tokens":212}},"tokens_in":733,"tokens_out":288,"duration_ms":3437,"temperature":1.0,"reasoning_tokens":212,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:31:55.498435+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a time-dependent 3D stellar-wind model for Trappist-1 and decompose the magnetic-field fluctuations at the planet's orbit in Fourier space: if the variance is dominated by timescales much shorter than the 1.51-day orbital period (CME crossings or advection past the planet), then the sinusoidal, orbital-frequency assumption at the heart of the skin-depth calculation fails, and the predicted penetration and heating profiles do not apply.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the induction-equation framework and the analytical skin-depth solution this paper generalizes to arbitrary conductivity profiles."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the parametrized collision frequencies used to build the parallel, Pedersen, and Hall conductivities."},{"cited_title":"G., Noack, L., Johnstone, C","cited_arxiv_id":null,"evidence_quote":"Established induction heating of exoplanet interiors from time-varying fields, the baseline the upper-atmosphere screening result qualifies."},{"cited_title":"2024, ApJ, 962, 157, publisher: IOP ADS Bibcode: 2024ApJ...962..157C","cited_arxiv_id":null,"evidence_quote":"Prior upper-atmosphere Ohmic-heating model whose heating rates are compared and bounded by the Qmax saturation."},{"cited_title":"From Stellar Coron{\\ae} to Gyrochronology: a theoretical and observational exploration","cited_arxiv_id":"2002.00696","evidence_quote":"Provides the stellar magnetic-field scaling used to estimate Qmax across the exoplanet population."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Dynamo scaling-law estimates of planetary magnetic moments that define the intermediate BP range of interest."},{"cited_title":"& Basri, G","cited_arxiv_id":null,"evidence_quote":"Observational estimate of Trappist-1's surface magnetic field used to infer plausible Bsw at the planet."}],"review_version":1}