{"id":"d6a273d2-0256-49bf-89e3-35527bfa4415","arxiv_id":"1908.06695","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A framework that combines stellar wind pressure, planetary magnetic field, and orbital distance to estimate the likelihood that an exoplanet retains its atmosphere.","lead":"This paper estimates how much of an exoplanet's atmosphere is exposed to erosion by its star's wind, based on the planet's magnetic field and orbit. It reports that tidally locked planets around small, active stars have very low chances of keeping their atmospheres, which matters for deciding which worlds could be habitable.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The tidal-locking headline rests on an unstated dynamo scaling: Table 3's Bpp collapse from 3 G to 0.07–0.52 G is never derived from the cited Olson & Christensen / Zuluaga scalings, so the central result cannot be checked.","rationale":"The reader's weakest-assumption analysis identified the same gate: the tidal-locking result depends on planetary magnetic moments computed from a dynamo scaling that is cited but never shown. My reading confirms this is the single most load-bearing concern. The paper has genuine independent support—Table 2 reproduces observed solar-system magnetopause radii, and the Monte Carlo uncertainty propagation for real cases is a reasonable treatment of observational errors—but none of that validates the 0.07–0.52 G Bpp values assigned to synchronous rotators. I considered the hand-calibrated likelihood function in Eq. (28) as an alternative concern; it is arbitrary but monotonic, so it would shift thresholds rather than overturn the qualitative α0 behavior. The unstated dynamo scaling, by contrast, determines whether the tidal-locking increase occurs at all. The appropriate remedy is the one the reader already states: write the magnetic-moment equation explicitly, show that it reproduces Table 3, and test sensitivity to reasonable alternative scalings. Since the reader's verdict is already CONDITIONAL and my concern reinforces that condition, I recommend no change to the verdict.","tokens_in":20092,"tokens_out":7535,"duration_ms":78888,"concrete_test":"Recompute Table 3's Bpp column by implementing the magnetic-moment scaling from Olson & Christensen (2006) and Zuluaga et al. (2013) referenced in §2.3 and applying it to the Table 3 grid (masses 0.08–1 M⊙, 1:1 spin-orbit, orbital periods listed). Then rerun Eq. (21) with a floor of Bpp=1 G for the tidally locked cases. If any α0 drops below 40° (or P rises above 0.12), the missing dynamo formula is load-bearing; if not, the claim is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that tidal locking raises the unprotected angle above 40° and drives retention probability below 0.12—is produced in Table 3 by a 40-fold drop in polar field Bpp, from 3 G for a freely rotating super-Earth to 0.07–0.52 G for synchronous cases. Section 2.3 says only that the magnetic moment is estimated 'following Olson & Christensen (2006)' and Zuluaga et al. (2013), depending on rotation state, but the actual dynamo scaling is never written down. No equation connects Bpp to the tidal-locking rotation period. Table A.3 similarly lists dimensionless 'Magnetic Moment' values (0.70 for GJ 2232 b vs 92.0 for Kepler-62 e) without units or formula. Because Eq. (21) enters as rM/RP ∝ Bpp^{1/3}, the entire tidal-locking result hinges on whether slow synchronous rotation really suppresses the dipole moment this strongly. The solar-system comparison in Table 2 supports the pressure-balance part of the pipeline but does not test the dynamo scaling, since all those planets rotate in known, non-synchronous states. The manuscript itself (§4.4) admits model uncertainties 'can be hardly evaluated.'","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents an analytic formalism for estimating the magnetopause radius rM/RP of Earth-like exoplanets orbiting low-mass stars (0.08–1.3 M⊙), based on pressure balance between the planetary dipole field and the stellar wind's dynamic and magnetic pressures. From rM/RP the paper computes an auroral aperture angle α0 and defines a Gaussian likelihood P(α0) for atmospheric retention, calibrated to a young-Earth scenario and to a Lammer et al. critical magnetopause radius. The authors validate the magnetopause calculation against solar system planets (Table 2), apply the formalism to Earth-like planets at the habitable-zone boundaries of 33 stars, simulate tidally locked super-Earths around synthetic stars, and analyze 26 known telluric exoplanets. The headline results are a saturated/unsaturated stellar-regime split at α0 ≈ 36°, a dramatic increase in unprotected angle for tidally locked planets, and near-zero retention likelihood for planets in sub-Alfvénic winds.","tokens_in":20366,"tokens_out":4861,"duration_ms":51020,"significance":"If the magnetopause calculation is accepted, the paper offers an inexpensive first-order screening tool for exoplanet atmospheric retention and a concrete component for the SEPHI habitability index. The strengths include the explicit pressure-balance equations, the favorable comparison with observed solar system magnetopause radii in Table 2, the comparison of predicted Alfvén radii with 3D MHD results for Proxima Cen and TRAPPIST-1, and the Monte Carlo propagation of observational uncertainties in Section 4.4. The central tidal-locking conclusion, however, rests on an unstated planetary dynamo scaling that connects synchronous rotation to a 40-fold reduction in polar field strength; until that step is written down and validated, the dramatic tidal-locking effect cannot be considered established.","major_comments":[{"comment":"The planetary magnetic induction for tidally locked planets is the load-bearing input of the central claim, but the paper never states the dynamo scaling that maps rotation to magnetic moment. Section 2.3 says the magnetic moment is estimated 'following Olson & Christensen (2006)' and 'depending on whether the planet is tidally locked ... or rotates freely,' yet no equation connects Bpp to the synchronous rotation period. In Table 3, Bpp drops from 3 G for free rotation to 0.07–0.52 G for tidally locked cases; in Table A.3, the dimensionless 'Magnetic Moment' spans 0.03 to 92.0 with no formula or units. Because Eq. (21) enters as rM/RP ∝ Bpp^{1/3}, the reported jump of the unprotected angle from ≲30° to >42° for tidal locking is entirely controlled by this omitted scaling. Please supply the explicit Olson–Christensen/Zuluaga formulation, the assumed core radius and mass scalings, and a validation or sensitivity test of the factor-of-ten reduction.","section":"Section 2.3, Tables 3 and A.3"},{"comment":"The likelihood values presented as headline results (probabilities below 0.12 in Table 3 and 'almost zero' in the abstract) are largely a restatement of the two hand-set thresholds α0,young earth = 25° and α0,max = 50°. With σ = (50−25)/3 = 8.33°, any planet with α0 above about 40° automatically receives P ≲ 0.12, and once α0 exceeds about 45° the probability is close to zero. Thus the likelihood is a calibrated output rather than an independent physical prediction. The paper should either present P as explicitly conditional on these choices with a sensitivity analysis, or replace the Gaussian with a more physical escape criterion; without this, the conclusion that tidally locked planets have a 'low likelihood' of retaining their atmospheres is not a standalone result.","section":"Section 3, Eq. (28)"},{"comment":"The Monte Carlo uncertainty analysis propagates only observational errors. The text in Section 4.4 concedes that uncertainties from the assumptions and simplifications 'can be hardly evaluated,' but this concession is not reflected in the precision with which likelihoods are quoted in Table A.4. The largest systematic uncertainties are in the stellar mass-loss and magnetic-field scalings (Eqs. 5, 6, 10) and in the planetary dynamo input, and none of these enter the 25,000-realization Gaussian propagation. In addition, Eq. (10) is stated to be valid only for stars older than 600 Myr with FX < 1e6 erg cm^-2 s^-1, yet it is applied to active M dwarfs such as TRAPPIST-1 and Ross 128 without an age check. Please quantify or at least bound these systematics before presenting the likelihood values to two significant figures.","section":"Section 4.4"}],"minor_comments":[{"comment":"Typo: 'superssonic' should be 'supersonic'.","section":"Section 4.2"},{"comment":"Typo: 'magentized' should be 'magnetized'.","section":"Section 5"},{"comment":"Typo: 'GRENN' should be 'GREEN'.","section":"Table A.3"},{"comment":"The entry for DS Leo ('15.7 7 3770') appears malformed; the rotation-period column contains an extra '7'.","section":"Table A.1"},{"comment":"The header dates ('Received September 15, 1996; accepted March 16, 1997') are inconsistent with the arXiv submission and should be corrected.","section":"Manuscript header"},{"comment":"The footnote 'e' used for estimated masses and radii should be defined in the table caption rather than left as an implicit marker.","section":"Table A.3"}],"recommendation":"major_revision","confidential_remarks":"The decisive issue for me is the missing dynamo scaling in Section 2.3. If the authors can provide the explicit Olson–Christensen/Zuluaga equation and show that the tidal-locking reduction in Bpp is robust, the paper is publishable as a first-order estimator. If that scaling cannot be supplied or validated, the central tidal-locking claim would collapse to an assumption. I saw no misconduct or citation problem, but the omission of the exact scaling is unusual for a result that states 'low likelihood of keeping its atmosphere' as a quantitative finding."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You'll want to know this about Rodriguez-Mozos & Moya: the core pipeline is actually useful, and the solar-system check in Table 2 is a real asset. The authors assemble standard pieces—Sadeghi Ardestani's wind scalings, Vidotto's unprotected angle, Olson-Christensen for magnetic moments—into a transparent first-order estimator for magnetopause radius and the fraction of atmosphere exposed to stellar wind. They validate it against Mercury through Neptune and get rM/RP within the observed ranges. That's not nothing. The new bit is the retention likelihood (Eq. 28) and the reported split between saturated and unsaturated stars, where auroral angles sit above or below 36 degrees with apparently different slopes in Fig. 1.\n\nThe soft spots are real, though, and one is load-bearing. The headline claim—tidally locked planets have unprotected angles above 40 degrees with retention likelihood below 0.12—depends entirely on the planetary magnetic moment collapsing from 3 G for a free rotator to 0.07–0.52 G for synchronous rotators in Table 3. Section 2.3 cites Olson & Christensen and Zuluaga et al. but never writes down the dynamo scaling that connects Bpp to rotation period. You literally cannot check the central result from the paper as it stands. If that scaling is wrong, the tidal-locking conclusion collapses. The solar-system validation doesn't help here, because all those planets rotate non-synchronously, so it never exercises this part of the model.\n\nSecond, the likelihood function is calibrated on Earth's history with two hand-picked angles—25 and 50 degrees—and a sigma set to one-third of their difference. That makes Eq. 28 a reasonable ranking tool, but the output probabilities are calibration-dependent, not absolute. The authors don't oversell this explicitly, but the abstract's \"low likelihood\" phrasing invites readers to treat it as a physical probability.\n\nThird, the linear relations in Fig. 1 have no error bars, no regression fit, no significance test. The 36-degree split is a consequence of the chosen Gaussian threshold, not an independent discovery. The manuscript even concedes (§4.4) that uncertainties from the model assumptions \"can be hardly evaluated.\"\n\nSo this is a promising estimator with a missing equation in the middle. The fix is straightforward: write down the Bpp scaling, justify the parameters, add scatter and a significance test to Fig. 1, and present the probabilities as relative rankings.\n\nWho is this for? People building habitability indices like SEPHI and anyone doing target selection for atmospheric observations who wants a first cut at magnetic shielding. It's a tool paper, not a mechanism paper.\n\nMy take: send it to a serious referee. The flaws are addressable, the solar-system check earns it a real look, and the field could use a transparent estimator like this—once the dynamo scaling is on the table.","headline":"A genuinely useful first-order estimator for magnetopause radius and retention likelihood, but the headline tidal-locking result rests on a dynamo scaling the paper never writes down.","tokens_in":20901,"tokens_out":2369,"would_cite":true,"duration_ms":22892,"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":"The paper claims that a first-order magnetosphere-radius estimate from stellar-wind pressure can rank exoplanets by atmospheric retention, and predicts that tidally locked planets in the habitable zones of low-mass stars have retention…","keywords":["exoplanet atmospheres","stellar wind erosion","magnetosphere radius","tidal locking","habitable zone","M dwarfs","planetary magnetic fields","atmospheric retention"],"falsifier":"Take a tidally locked rocky exoplanet around a $0.1$--$0.3\\,M_\\odot$ star and measure its polar magnetic field through radio emission or related diagnostics. If the field comes out near $1$ G or higher rather than the $0.07$--$0.5$ G the model assigns, the sharp drop in retention probability for synchronized planets is falsified; likewise, a confirmed thick atmosphere on a close-in sub-Alfvénic-impact planet such as TRAPPIST-1e would contradict the near-zero retention claim.","tokens_in":19870,"feed_emoji":"🪐","tokens_out":10562,"duration_ms":99966,"temperature":0.7,"pith_summary":"This paper is trying to establish that the fate of an exoplanet's atmosphere under stellar-wind erosion can be estimated, at first order, from a pressure-balance calculation: the radius where the planetary magnetic field balances the stellar wind sets the fraction of the surface left unprotected through open polar field lines. The authors carry this through for Earth-like planets at both edges of the habitable zone around G, K, M, and very low-mass stars, and translate the unprotected fraction into a likelihood of retaining an atmosphere. Their central result is a division by rotation state: freely rotating planets with Earth-strength or stronger fields mostly keep their atmospheres, while tidally locked planets show unprotected angles above 40 degrees and retention likelihoods below about 0.12; planets whose magnetospheres are struck in the sub-Alfvénic regime have near-zero retention likelihoods. If these predictions hold, the formalism gives a cheap screening tool for habitability that can be folded into the SEPHI index and applied to known exoplanet catalogs.","feed_headline":"Tidally locked exoplanets probably cannot keep atmospheres","feed_subtitle":"Pressure-balance model gives near-zero retention odds for synchronous planets and sub-Alfvénic impacts.","key_machinery":"The load-bearing object is the magnetopause stand-off radius ratio $r_M/R_P$, obtained from pressure balance: $r_M/R_P = [B_{p,\\mathrm{eq}}^2/(8\\pi(P_{\\mathrm{dyn}}+B_{sa}^2))]^{1/6}$, with $B_{p,\\mathrm{eq}}\\approx B_{pp}/2$. This radius sets the unprotected polar cap through $\\alpha_0 = \\arcsin\\sqrt{R_P/r_M}$, and the paper converts $\\alpha_0$ into a retention likelihood anchored to the young Earth. The stellar inputs (mass loss, magnetic induction, wind speed, Alfvén radius) come from rotation-based scalings that split stars into saturated and unsaturated regimes; the planetary field for tidally locked planets comes from a dynamo scaling tied to rotation. Every prediction in the paper flows through this single radius, which is why the authors first validate it against solar-system magnetopause measurements.","core_discovery":"On the paper's own terms, the discovery is that the magnetopause radius of an exoplanet, computed from a balance between the planet's equatorial magnetic pressure and the sum of the stellar wind's dynamic pressure and the stellar magnetic field at the orbit, controls atmospheric survival through the unprotected polar-cap angle $\\alpha_0 = \\arcsin\\sqrt{R_P/r_M}$. Applying this to an Earth-like planet at the habitable-zone boundaries, the authors find two clean behavioral regimes: fast-rotating (saturated) host stars push Earth-like planets above $36^\\circ$ at the inner edge with retention likelihoods below 0.4, while slower (unsaturated) hosts keep them below $36^\\circ$; and tidal locking collapses the magnetopause to roughly $1.1$--$2.2\\,R_P$, driving $\\alpha_0$ above $40^\\circ$ and the retention likelihood below 0.12. In sub-Alfvénic impacts, which occur for very close-in planets, $r_M/R_P$ falls to about $1.2$--$1.9$ and the retention likelihood essentially vanishes regardless of host saturation.","pith_inferences":["Editorial inference: if slow rotation suppresses planetary dynamos as strongly as assumed, then target lists for atmospheric characterization around M dwarfs should exclude or down-weight synchronized planets and instead prioritize planets on wider, non-synchronous orbits or spin-orbit resonances like 3:2.","Editorial inference: the Gaussian likelihood anchored to the young Earth is a placeholder; the same magnetopause machinery could be coupled to hydrodynamic escape codes to replace the $25^\\circ$ and $50^\\circ$ thresholds with physically computed erosion timescales.","Editorial inference: transients such as flares and coronal mass ejections, deliberately excluded here, would only shrink the magnetopause further and lengthen the time spent at large unprotected angles, so the retention likelihoods are likely upper bounds for active M dwarfs."],"forward_implications":["For stars in the saturated regime, Earth-like planets at the inner edge of the habitable zone have unprotected angles above $36^\\circ$ and retention likelihoods below 0.4, making host-star rotation state a first-order habitability filter.","Tidally locked super-Earths in the simulated sample have $r_M/R_P$ between 1.1 and 2.2 and retention likelihoods at most 0.12, implying that synchronous rotation is a strong anti-habitability factor for close-in planets.","When the stellar wind hits the magnetopause in the sub-Alfvénic regime, retention likelihood is essentially zero for every case studied, singling out very close-in planets such as the inner TRAPPIST-1 planets as atmosphere-poor.","The likelihood function, combined with the SEPHI index, ranks Kepler-186 f, Kepler-1229 b, and Kepler-442 b as the known telluric planets most likely to retain their atmospheres.","The same equations applied to 26 real telluric exoplanets reproduce the division: free-rotating planets have likelihood near unity, while tidally locked planets fall below 0.17."],"supporting_citations":[{"why":"Supplies the magnetopause radius and unprotected-angle formulas ($\\alpha_0 = \\arcsin\\sqrt{R_P/r_M}$) and the sample of M-type stars used in Section 4.2.","marker":"Vidotto et al. (2013)"},{"why":"Source of the pressure-balance equation for $r_M/R_P$ used throughout the paper.","marker":"Vidotto (2018)"},{"why":"Provides the dipole scaling for planetary magnetic induction from core radius and magnetic moment.","marker":"Olson & Christensen (2006)"},{"why":"Used for the planetary magnetic moment estimates, including cases of tidally locked planets.","marker":"Zuluaga et al. (2013)"},{"why":"Refines the planetary magnetic moment estimate for slow rotators.","marker":"Zuluaga & Bustamante (2018)"},{"why":"Gives the Alfvén radius, angular momentum loss, and Rossby-number-based mass-loss and magnetic-field scalings.","marker":"Sadeghi Ardestani et al. (2017)"},{"why":"Provides the X-ray flux--mass-loss relation used for low-mass stars below 0.5 solar masses.","marker":"Wood et al. (2002, 2005)"},{"why":"Supplies the saturation threshold and coronal temperature scaling used for wind speed and sound radius.","marker":"Johnstone & Güdel (2015)"},{"why":"Defines the habitable-zone boundaries D2 and D3 used in the planet placement calculations.","marker":"Kopparapu et al. (2013, 2014)"},{"why":"Defines the SEPHI habitability index into which the retention likelihood is incorporated.","marker":"Rodríguez-Mozos & Moya (2017)"}],"fun_headline_variants":["Sub-Alfvenic winds make atmosphere retention nearly impossible","Tidal locking collapses magnetopause, slashes retention odds","Magnetopause radius dictates exoplanet atmospheric survival","Habitable-zone edges show extreme atmospheric erosion regimes","Pressure balance model forecasts which planets keep atmospheres"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The tidal-locking result rests on the assumption that synchronous rotation weakens the planet's dynamo-generated magnetic field to roughly $B_{pp} \\approx 0.07$--$0.5$ G; the paper adopts this scaling without displaying the equation, and if a slow rotator kept a field near $1$ G the predicted loss of atmosphere would not occur.","fun_headline_variants_meta":{"raw":{"variants":["Sub-Alfvenic winds make atmosphere retention nearly impossible","Tidal locking collapses magnetopause, slashes retention odds","Magnetopause radius dictates exoplanet atmospheric survival","Habitable-zone edges show extreme atmospheric erosion regimes","Pressure balance model forecasts which planets keep atmospheres"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000773,"raw_usage":{"total_tokens":3507,"prompt_tokens":1119,"completion_tokens":2388,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":735,"completion_tokens_details":{"reasoning_tokens":2311}},"tokens_in":735,"tokens_out":2388,"duration_ms":19832,"temperature":1.0,"reasoning_tokens":2311,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:37:11.551149+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a tidally locked rocky exoplanet around a $0.1$--$0.3\\,M_\\odot$ star and measure its polar magnetic field through radio emission or related diagnostics. If the field comes out near $1$ G or higher rather than the $0.07$--$0.5$ G the model assigns, the sharp drop in retention probability for synchronized planets is falsified; likewise, a confirmed thick atmosphere on a close-in sub-Alfvénic-impact planet such as TRAPPIST-1e would contradict the near-zero retention claim.","supporting_citations":[{"cited_title":"A., Jardine, M., Morin, J., et al.\\ 2013, , 557, A67","cited_arxiv_id":null,"evidence_quote":"Supplies the magnetopause radius and unprotected-angle formulas ($\\alpha_0 = \\arcsin\\sqrt{R_P/r_M}$) and the sample of M-type stars used in Section 4.2."},{"cited_title":"A.\\ 2018, Handbook of Exoplanets, 26","cited_arxiv_id":null,"evidence_quote":"Source of the pressure-balance equation for $r_M/R_P$ used throughout the paper."},{"cited_title":"R.\\ 2006, Earth and Planetary Science Letters, 250, 561","cited_arxiv_id":null,"evidence_quote":"Provides the dipole scaling for planetary magnetic induction from core radius and magnetic moment."},{"cited_title":"I., Bustamante, S., Cuartas, P","cited_arxiv_id":null,"evidence_quote":"Used for the planetary magnetic moment estimates, including cases of tidally locked planets."},{"cited_title":"I., & Bustamante, S.\\ 2018, , 152, 55","cited_arxiv_id":null,"evidence_quote":"Refines the planetary magnetic moment estimate for slow rotators."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Alfvén radius, angular momentum loss, and Rossby-number-based mass-loss and magnetic-field scalings."},{"cited_title":"E., M \\\"u ller, H.-R., Zank, G","cited_arxiv_id":null,"evidence_quote":"Provides the X-ray flux--mass-loss relation used for low-mass stars below 0.5 solar masses."},{"cited_title":"P., & G \\\"u del, M.\\ 2015, , 578, A129","cited_arxiv_id":null,"evidence_quote":"Supplies the saturation threshold and coronal temperature scaling used for wind speed and sound radius."},{"cited_title":"M., & Moya, A.\\ 2017, , 471, 4628","cited_arxiv_id":null,"evidence_quote":"Defines the SEPHI habitability index into which the retention likelihood is incorporated."}],"review_version":1}