{"id":"76f08323-7c32-4491-8c7f-6a7d95a74187","arxiv_id":"2506.18972","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Radius inflation in young low-mass stars follows a saturated Rossby-number relation like standard activity indicators, and depends directly on measured starspot covering fraction.","lead":"A study of 261 young, low-mass stars in four clusters finds that stellar radius inflation is tightly linked to rotation and magnetic starspots, following the same pattern as stellar activity. The result supports the idea that magnetism, not rotation alone, is the direct cause of inflated radii in young stars.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Same-fit Teff–fspot covariance may manufacture the inflation–spot correlation; injection-recovery is required to confirm.","rationale":"The paper's central empirical claims are the Rossby-number dependence of radius inflation and a direct dependence on starspot covering fraction. The Rossby relation is supported by measured rotation periods and a four-cluster sample, but its interpretation is complicated by the use of fspot-dependent convective turnover timescales from the same SPOTS models used for the comparison isochrones. The more load-bearing, and less protected, claim is the direct inflation–fspot relation, because the radius and fspot are not measured independently: R comes from the flux-weighted Teff of the same two-temperature fit, with R ∝ Teff^{-2}. The paper does not demonstrate that the fitted Teff and fspot are uncorrelated. Two-temperature fits can have degeneracies among spot contrast, filling fraction, and continuum temperature, and the authors do not report the covariance or an injection-recovery calibration. If such a covariance exists, it would manufacture both the inflation–fspot correlation and the apparent luminosity preservation, which is nearly built into the radius definition. This does not make the paper wrong, but it makes the headline claim conditional on a check that is straightforward to run. The comparison with SPOTS models in Figures 3–4 is an interesting consistency test, but because those models also provide the τ_CZ values and the base isochrones, it cannot independently validate the same-fit relation. An injection-recovery test of the LEOPARD pipeline is the cleanest way to settle the concern. If the covariance is small, the empirical result stands and the paper should be accepted; if not, the central claim needs reanalysis. The reader's CONDITIONAL verdict is therefore appropriate and unchanged.","tokens_in":13271,"tokens_out":7697,"duration_ms":91111,"concrete_test":"Run an injection-recovery test with the LEOPARD two-temperature fitting pipeline. Generate synthetic APOGEE H-band spectra with known T_phot, T_spot, fspot, and vsini spanning the sample's parameter range; add realistic noise; then recover fspot and flux-weighted Teff using the same fitting procedure and compute R with the Section 2.2 SED recipe. Measure the covariance between recovered Teff and fspot, and propagate that covariance into R_infl. If the pipeline covariance alone reproduces a significant fraction of the observed inflation–fspot slope, the direct starspot-inflation claim is unsecured. If the recovered Teff–fspot covariance is negligible, the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 2.2, the stellar radius is computed from the Stefan–Boltzmann relation using the flux-weighted Teff produced by the same two-temperature LEOPARD fit that yields fspot. At fixed Fbol and parallax, R ∝ Teff^{-2}, so any covariance between the fitted Teff and fspot—for example from degeneracies among photospheric temperature, spot temperature contrast, and filling fraction—propagates directly into R and then into R_infl. If high-fspot solutions are preferentially paired with low recovered Teff, the claimed direct dependence of radius inflation on starspot covering fraction (Section 3.3, Figures 3–4) could be generated by the measurement pipeline rather than by stellar physics. Because R is constructed so that 4πR^2σTeff^4 equals the measured Fbol, the 'luminosity preservation' reported in Section 3.2 is also not an independent check. The paper does not report covariance diagnostics or injection-recovery tests for the LEOPARD fits, leaving this confound unquantified. This is load-bearing for the central claim that starspots directly drive the inflation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes 261 low-mass stars in Upper Sco, alpha Per, the Pleiades, and Praesepe, combining APOGEE two-temperature spectroscopic fits (which yield fspot and Teff), SED-based radii, literature rotation periods, and Gaia parallaxes. The authors report a strong relationship between radius inflation (relative to non-magnetic isochrones) and Rossby number, with a saturated regime for log Ro less than about -1.4 and a declining power-law at larger Ro. They also report a direct dependence of inflation on starspot covering fraction, a temperature suppression that preserves bolometric luminosity, and agreement with the SPOTS magnetic stellar model predictions, albeit with a small offset. The sample and methodology are described in Sections 2 and 3, with MCMC fits in Appendices A and B.","tokens_in":13610,"tokens_out":8662,"duration_ms":89616,"significance":"If the central claims hold, this would be a substantial step toward unifying radius inflation with rotation and magnetic activity across the pre-main sequence to ZAMS, and would provide strong empirical support for starspot-based stellar evolution models. The paper's strengths include a large, homogeneous sample spanning four well-studied clusters; use of empirical rotation periods and Gaia parallaxes; a careful MCMC fitting procedure with variance underestimation; and a multi-criterion binary rejection strategy. However, the validity of the fspot-Rossby-inflation relation rests on assumptions about the two-temperature fits and the choice of convective turnover timescale that are not yet demonstrated.","major_comments":[{"comment":"The radius is derived from the flux-weighted Teff produced by the same two-temperature LEOPARD fit that yields fspot. Because R is proportional to Teff^{-2} at fixed bolometric flux and parallax, any covariance between the fitted Teff and fspot--for example, from degeneracies between spot temperature contrast and filling fraction--propagates directly into R and then into R_infl. The claimed direct dependence of radius inflation on starspot covering fraction (Section 3.3, Figures 3-4) could therefore be generated by the measurement pipeline rather than by stellar physics. The paper does not report covariance diagnostics, corner plots for the two-temperature fits, or injection-recovery tests; the error analysis in Section 2.2 assumes Teff and Fbol errors are independent and does not address Teff-fspot covariance. Please add such tests, or recompute radii with an independently determined Teff, to demonstrate that the fspot-inflation relation is not an artifact of the fitting procedure.","section":"Section 2.2; Section 3.3"},{"comment":"The statement that temperature suppression balances radius inflation so as to preserve bolometric luminosity is a tautology of the radius derivation. Since R is computed from the same Fbol and Teff via the Stefan-Boltzmann relation, the luminosity is exactly equal to the measured bolometric flux by construction. Consequently, the observed correlation between R_infl and Teff suppression in Figure 2 (bottom) is not an independent empirical check and cannot be used to support the claim of luminosity preservation. Please reframe this result, and if the goal is to test luminosity conservation, derive radii from independent measurements (e.g., interferometry or eclipsing binaries) or use Teff from a fit that does not adopt the same Fbol.","section":"Section 3.2; Figure 2"},{"comment":"The Rossby number is computed using convective turnover timescales taken from the SPOTS models, which are also the comparison theory in Section 3.3. Any fspot-dependent error in tau_CZ will shift the location of the saturation break and can create or distort a correlation between R_infl and Ro even if rotation is not the underlying driver. To demonstrate that the inflation-Rossby relation is robust, please recompute Ro with an independent tau_CZ calibration (e.g., Noyes et al. 1984 or a non-magnetic standard model) and show that the linear and saturated regimes, and the break near log Ro approximately -1.4, are unchanged. The current approach makes the test of the SPOTS models partly self-referential.","section":"Section 2.2; Section 2.3; Section 3.1"}],"minor_comments":[{"comment":"The equation for the fiducial fit is garbled and appears to have sign errors; the printed expression does not reproduce the stated decline from 28% to 5.5% at Ro = 1. Please correct the equation and the surrounding text.","section":"Figure 2 caption"},{"comment":"The text says the relation declines by a factor of 5 to 5.5% at a Rossby number of unity; with the stated parameters (saturation y = 0.276 and slope m = 1.154) this only holds if the exponent has a negative sign. Please unify the mathematical expression and the textual description.","section":"Appendix B"},{"comment":"The claim that the inflation signal is a universal function of the stellar activity level is stronger than the data support, since Figure 2 shows systematic vertical offsets between clusters that the authors attribute to age differences; please soften or qualify this statement.","section":"Section 4"},{"comment":"The paper states that errors in Teff and Fbol are assumed to be normally distributed and independent; please also provide the covariances among the LEOPARD output parameters or a reference to where they are published, since these covariances are essential for assessing the fspot-Teef degeneracy.","section":"Section 2.2"},{"comment":"The binary rejection completeness of at least 85% is stated as an assumption; please provide the source of this estimate and consider showing the sensitivity of the main relations to a stricter binary cut.","section":"Appendix A"},{"comment":"There are several typographical and formatting issues, including 'inreased' in Section 3.2 and the corrupted phrase 'Rspotted SPOTS fspot = 0/RSPOTS fspot = 0' in the Figure 2 caption; please proofread the manuscript.","section":"Various"}],"recommendation":"major_revision","confidential_remarks":"The paper is potentially interesting, but the central claims require stronger validation of the two-temperature fitting degeneracy and the use of SPOTS-based tau_CZ. The reliance on 'Cao et al. (in prep)' for the LEOPARD catalog is a reproducibility concern; the authors should either include sufficient methodological detail in an appendix or provide a public release of the fitting code and covariance information. I recommend major revision rather than rejection, as the identified issues are addressable with additional analysis and transparency."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I'll be direct. This is the most serious empirical attempt to tie radius inflation to rotation and starspots, and the inflation-Rossby relation is probably real. But the headline claim of a direct starspot-fraction dependence is not yet proven, because the radius and fspot come from the same two-temperature fit.\n\nWhat's new: 261 stars in four clusters spanning 10–670 Myr, with rotation periods, SED-based radii, and APOGEE two-temperature spot filling fractions. The inflation-Rossby relation shows a clean linear-to-saturated break near log Ro ~ -1.4, reminiscent of activity-rotation relations. That's a real advance over the single-cluster studies by Somers & Stassun (2017) and Jaehnig et al. (2019), and extending the signal to 10 Myr in Upper Sco is new. The MCMC fits are careful, the binary rejection is explicit, and the paper doesn't hide the scatter.\n\nSoft spots. The radius uses the flux-weighted Teff from the same LEOPARD fit that yields fspot. With R ∝ Teff^-2 at fixed Fbol and parallax, any covariance between fitted Teff and fspot flows straight into the inflation-spot correlation. The luminosity-preservation check in Section 3.2 is not independent: R is constructed so that 4πR^2σTeff^4 equals the measured Fbol, so L is conserved by definition. And the Rossby numbers use tau_CZ from the SPOTS models, set by the observed fspot, with SPOTS as the comparison theory—so that leg is partly self-referential. These issues don't threaten the empirical inflation-Rossby relation, which uses independent rotation and radius data, but they do hang over the direct starspot-dependence claim. The paper offers no injection-recovery or covariance diagnostics for the fits. A serious referee should ask for that.\n\nProportion: the core result—radius inflation scaling with Rossby number, with saturation—looks robust across clusters. The fspot dependence is plausible but unverified. Age assumptions shift cluster offsets but don't erase the pattern.\n\nThis paper is for people working on stellar activity, pre-main-sequence evolution, and exoplanet radius systematics. It deserves full peer review and likely heavy revision, not a desk reject. I'd want covariances addressed before I'd cite the fspot dependence, but the sample and the Rossby relation are worth having.\n\nRecommendation: send to review; conditional on the same-fit covariance being tested.","headline":"Strong new inflation-Rossby relation; the direct starspot dependence needs covariance controls before it can be trusted.","tokens_in":14000,"tokens_out":3762,"would_cite":false,"duration_ms":41584,"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":"A study of 261 low-mass stars in four open clusters aged 10–670 Myr shows radius inflation is a strong function of the Rossby number—with linear and saturated regimes breaking at log Ro ≈ −1.4—and that inflation depends directly on the…","keywords":["radius inflation","starspots","Rossby number","pre-main-sequence stars","open clusters","stellar rotation","magnetic activity","stellar evolution models"],"falsifier":"Measure radii of a subset of the same stars with a route that does not use the two-temperature effective temperature—eclipsing binaries or long-baseline interferometry—and check whether the inflation–Rossby and inflation–spot correlations persist, or recompute Rossby numbers with convective turnover timescales from non-SPOTS models and see whether the log Ro ≈ −1.4 break survives.","tokens_in":13068,"feed_emoji":"⭐","tokens_out":7732,"duration_ms":79125,"temperature":0.7,"pith_summary":"This paper tries to establish that the mysterious \"radius inflation\" of active low-mass stars—radii 10–15 percent larger than standard models predict—is a magnetic phenomenon driven by starspots. Using 261 stars in clusters aged 10–670 Myr, it reports a tight relation between inflation and the Rossby number, with a linear regime and a saturated regime breaking at log Ro ≈ −1.4, closely resembling the canonical activity–rotation relation. It demonstrates for the first time that the inflation depends directly on the measured starspot covering fraction. It also finds that effective temperature drops as the Rossby number drops, and that this temperature suppression balances the inflated radius so the bolometric luminosity stays fixed. If correct, the result unifies inflation across the whole pre-main-sequence age range and points to magnetic starspot models as the right basis for interpreting young, low-mass stars.","feed_headline":"Magnetic starspots directly inflate young stars' radii","feed_subtitle":"In 261 cluster stars aged 10-670 Myr, inflation follows the Rossby number and saturates like other activity proxies.","key_machinery":"The argument runs on four linked objects: the fractional radius inflation measured against the $f_{\\rm spot}=0$ non-magnetic isochrone; the Rossby number $\\mathrm{Ro} = P_{\\rm rot}/\\tau_{\\rm CZ}$, using rotation periods from the literature and convective overturn timescales $\\tau_{\\rm CZ}$ taken from the spotted SPOTS models; the starspot filling fraction $f_{\\rm spot}$ and two-temperature effective temperature from the spectroscopic fitting of APOGEE H-band spectra; and the SPOTS/ROTEVOL magnetic isochrones that act as the comparison theory. The same two-temperature fit supplies both the temperature that sets the radius (through the SED pseudo-interferometry with Gaia parallaxes) and the spot fraction, so the machinery ties the measurement of inflation to the measurement of its suspected cause.","core_discovery":"The central claim is that the degree of radius inflation across a large sample of young cluster stars is a universal function of magnetic activity: it rises steeply with decreasing Rossby number, saturates at log Ro ≈ −1.4, and depends directly on the starspot filling fraction $f_{\\rm spot}$. The fitted relation is $R_{\\rm obs}/R_{\\rm isoc} = 0.276\\,\\exp(1.154(\\max(\\log \\mathrm{Ro}, -1.396)+1.396))$, so the most active stars are about 28 percent larger than non-magnetic isochrones predict, declining by a factor of five by Rossby number unity. Stars' effective temperatures are suppressed as the Rossby number decreases, and the combination of larger radius and lower temperature preserves luminosity, exactly as expected if the core luminosity is radiated through a spotted, less effective surface. The observed inflation correlates with the radii predicted by the spotted SPOTS stellar evolution models at better than 99.9 percent confidence (Kendall's $\\tau = 0.492$), with a slope consistent with unity ($0.99 \\pm 0.08$) and a small offset of $0.040 \\pm 0.003\\,R_\\odot$; the models capture the inflation but under-predict it by about 0.27 dex for the most active stars. The paper concludes that rotation is the underlying driver, but magnetism—expressed through starspots—is the most likely direct cause.","pith_inferences":["If the inflation–spot relation is genuine, radius inflation should vary over individual activity cycles as spot coverage changes, which long-baseline photometric and spectroscopic monitoring of the same stars could test directly.","The paper's suggestion that different activity proxies saturate at different Rossby numbers could be checked within this same sample using X-ray or H-alpha indices, potentially probing the pressure or depth scale that anchors each magnetic mechanism.","Applying the fitted inflation relation as a correction factor would remove an activity-dependent bias from exoplanet transit radii and from lithium-based age estimates for young clusters, which the paper points to but does not quantify.","Because the Rossby numbers use convective timescales from the very spotted models being tested, an independent determination of $\\tau_{\\rm CZ}$ from non-magnetic models would clarify whether the saturation break at about −1.4 is a physical feature or partly a model artifact."],"forward_implications":["Radii of active young and low-mass stars, and therefore exoplanet radii, carry an activity-dependent bias of up to roughly 25 percent if magnetism is ignored in stellar characterization.","Radius inflation joins the canonical activity-rotation framework, saturating at a lower Rossby number ($\\log \\mathrm{Ro} \\approx -1.4$) than other activity proxies studied so far.","Spotted stellar evolution models can predict the radii of inflated pre-main-sequence stars to within about 0.04 solar radii, supporting their routine use in stellar characterization workflows.","Isochronal age discrepancies reported for pre-main-sequence clusters are largely explained by magnetic starspots rather than by data artifacts, potentially shifting cluster age estimates.","The inflation signal is a universal function of activity level across 10–670 Myr, not a peculiarity of any one cluster."],"supporting_citations":[{"why":"Supplies the SPOTS spotted stellar isochrones used as the magnetic comparison theory and the convective turnover timescales used to define Rossby number.","marker":"Somers et al. 2020"},{"why":"Provides the two-temperature spectroscopic fitting and LEOPARD catalog from which starspot filling fractions and effective temperatures are measured.","marker":"Cao & Pinsonneault 2022"},{"why":"Established the SED-based radius inflation analysis in the Pleiades that this study extends across four clusters and a wider age range.","marker":"Somers & Stassun 2017"},{"why":"Provides the pseudo-interferometric SED fitting machinery used to derive stellar radii from bolometric flux and Gaia parallaxes.","marker":"Stassun & Torres 2016"},{"why":"Supplies the Upper Scorpius rotation period catalog used to compute Rossby numbers.","marker":"Rebull et al. 2018"},{"why":"Supplies the Pleiades rotation period catalog used to compute Rossby numbers.","marker":"Rebull et al. 2016"},{"why":"Supplies the Praesepe rotation period catalog used to compute Rossby numbers.","marker":"Rebull et al. 2017"},{"why":"Supplies the alpha Persei rotation period catalog and the 80 Myr age estimate adopted for the cluster.","marker":"Boyle & Bouma 2023"},{"why":"Established the empirical relation between radius inflation, temperature suppression, and H-alpha activity that this paper re-derives at the sample level.","marker":"Stassun et al. 2012"},{"why":"Provides the canonical rotation–activity saturation relation that the inflation–rotation relation is directly compared to.","marker":"Wright et al. 2011"}],"fun_headline_variants":["Stellar radius inflation tied to starspot coverage","Young stars' radii inflate with magnetic starspots","Rotation and starspots drive stellar radius inflation","Radius inflation in 261 young stars hinges on magnetic activity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The inflation and the spot fraction come out of the same two-temperature fit, so any covariance in that fit could manufacture the inflation–spot correlation, and the Rossby numbers lean on convective timescales from the very spotted models the paper is testing.","fun_headline_variants_meta":{"raw":{"variants":["Stellar radius inflation tied to starspot coverage","Young stars' radii inflate with magnetic starspots","Rotation and starspots drive stellar radius inflation","Radius inflation in 261 young stars hinges on magnetic activity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000243,"raw_usage":{"total_tokens":1620,"prompt_tokens":1127,"completion_tokens":493,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":743,"completion_tokens_details":{"reasoning_tokens":432}},"tokens_in":743,"tokens_out":493,"duration_ms":5546,"temperature":1.0,"reasoning_tokens":432,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T23:14:25.114544+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure radii of a subset of the same stars with a route that does not use the two-temperature effective temperature—eclipsing binaries or long-baseline interferometry—and check whether the inflation–Rossby and inflation–spot correlations persist, or recompute Rossby numbers with convective turnover timescales from non-SPOTS models and see whether the log Ro ≈ −1.4 break survives.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the SPOTS spotted stellar isochrones used as the magnetic comparison theory and the convective turnover timescales used to define Rossby number."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Established the SED-based radius inflation analysis in the Pleiades that this study extends across four clusters and a wider age range."},{"cited_title":"G., & Torres , G","cited_arxiv_id":null,"evidence_quote":"Provides the pseudo-interferometric SED fitting machinery used to derive stellar radii from bolometric flux and Gaia parallaxes."},{"cited_title":"M., Stauffer , J","cited_arxiv_id":null,"evidence_quote":"Supplies the Upper Scorpius rotation period catalog used to compute Rossby numbers."},{"cited_title":"M., Stauffer , J","cited_arxiv_id":null,"evidence_quote":"Supplies the Pleiades rotation period catalog used to compute Rossby numbers."},{"cited_title":"M., Stauffer , J","cited_arxiv_id":null,"evidence_quote":"Supplies the Praesepe rotation period catalog used to compute Rossby numbers."},{"cited_title":"W., & Bouma , L","cited_arxiv_id":null,"evidence_quote":"Supplies the alpha Persei rotation period catalog and the 80 Myr age estimate adopted for the cluster."},{"cited_title":"G., Kratter , K","cited_arxiv_id":null,"evidence_quote":"Established the empirical relation between radius inflation, temperature suppression, and H-alpha activity that this paper re-derives at the sample level."}],"review_version":1}