{"id":"08fb5b97-294e-4255-a6d2-ca3089cd14e5","arxiv_id":"2607.24460","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":5,"one_line_summary":"After bias correction, fast-rotating YSOs have exponentially distributed cold-spot coverages while slow rotators lack small long-lived spots, implying rotation-dependent spot lifetimes.","lead":"Multi-region HOYS photometry of 144 young stars shows cold-spot coverage on fast rotators follows an exponential distribution, while slow rotators lack small long-lived spots. This implies spot lifetimes and magnetic surface structure depend on rotation rate, beyond simple disc braking.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The claimed ~150–200 d spot-lifetime threshold for slow rotators is asserted, not measured: the injection–recovery simulations used strictly coherent sinusoids, so the pipeline's actual sensitivity as a function of spot lifetime × period was never tested.","rationale":"The reader correctly identified the lifetime interpretation as the weakest assumption; I am sharpening where the soft spot actually is. It is not merely that alternatives (contrast, geometry) were considered qualitatively — the paper does address those in Sect. 4.2.1 — it is that the quantitative centerpiece of the interpretation (lifetimes ≲150–200 d) is an unmeasured assertion, and the only tool that could measure it (finite-lifetime injection–recovery) was not used, while the coherent-sinusoid injections that were used are structurally incapable of probing it. This matters because the abstract elevates the lifetime result to a conclusion (\"we conclude that small spots on slowly rotating YSOs have significantly shorter lifetimes\") and the supergranule-size discussion in Sect. 4.2.2 builds further on it. That said, the paper is a careful observational study: selection is documented, bias simulations cover the amplitude–period plane, Teff systematics are tested, alternatives are named rather than hidden, and the empirical distributions stand independent of the interpretation. I therefore do not recommend REJECT. I move from ACCEPT to CONDITIONAL: accept the empirical coverage distributions and the existence of the small-spot deficit, but condition the lifetime-threshold claim (and its supergranule corollary) on the finite-lifetime recovery test described above. If that test confirms a period-independent detection threshold near 150–200 d and still requires shorter intrinsic lifetimes for slow rotators, the full claim stands as published. A secondary, non-load-bearing note for any revision: state clearly whether the Fig. 9 fit treats slices or stars as the unit of independence, and provide the details of the Monte-Carlo fluke test (p<6×10⁻⁵) that are currently absent.","tokens_in":26790,"tokens_out":4417,"duration_ms":166039,"concrete_test":"Re-run the injection–recovery experiment of Sect. 3.1 with evolving spots: inject sinusoids whose amplitude decays (or whose phase drifts) with a finite lifetime τ drawn from 30–600 d, over the same period grid (0.75–19.6 d) and amplitudes corresponding to coverages 0.05–0.25, into the real HOYS light curves, and pass them through the identical pipeline (slice cuts, four periodogram methods, clustering, amplitude S/N>3, phase-alignment cuts). Report the recovery surface R(P, coverage, τ). Then forward-model Fig. 8 assuming slow rotators share the fast rotators' exponential coverage distribution with a period-independent τ distribution. If this reproduces the observed small-spot deficit, the \"shorter lifetimes on slow rotators\" conclusion weakens; if the deficit requires τ to be genuinely shorter at P>5.5 d, the claim is confirmed and the measured detection threshold replaces the asserted","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The strongest claim has two parts. The empirical part (exponential cold-spot coverage distribution for fast rotators; deficit of small spots at P>5.5 d in Fig. 8) is reasonably supported by the bias simulations of Sect. 3.1/Fig. 3 and survives the Teff and contrast checks in Sect. 4.2–4.2.1. The interpretive part — that the deficit implies small-spot lifetimes ≲150–200 d on slow rotators — rests on Sect. 4.2.2, where the threshold appears as an assertion: \"In practice, this implies that only spots with lifetimes ≳150–200 d are included in our analysis.\" No simulation in the paper establishes this number. The injection–recovery tests injected *strictly periodic, infinitely coherent* sinusoids into real light curves. Such a test measures detectability as a function of amplitude and period, but is blind by construction to spot evolution: a spot living 100 d on a P=15 d star produces a phase-unstable, non-sinusoidal signal across the ~9–12 cycles of a 6-month slice, which will fail the multi-method, multi-filter coherence and clustering cuts (Sect. 2.3) far more readily than an equally long signal on a P=2 d star (~90 cycles to average over). The recovery probability is therefore a function of (P, coverage, lifetime), and the 2-D (P, amplitude) grid of Fig. 3 cannot constrain the third axis. Two consequences: (1) the quoted 150–200 d threshold is uncalibrated — the real threshold is likely period-dependent, so the clean dichotomy \"fast rotators keep small spots, slow rotators lose them\" could be quantitatively distorted even if qualitatively real; (2) the exclusion of the \"residual selection effects\" alternative is weaker than stated, because the one simulation that would capture lifetime-dependent selection was not run. Note also that the Fig. 9 exponential fit is built from slices, not objects (620 amplitude sets over 144 YSOs; IC 5070 alone contributes ~1/3 of slices), so repeated measurements of the same star are treated as independent, overstating the fit's discrim","agreement_with_reader":"agree"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The paper presents a homogeneous analysis of HOYS multi-band (V,R,I) photometry across 25 star-forming regions, identifying 144 YSOs with robust periodic signals from 2047 candidate members. The period distribution is bimodal (55% fast, P<5.5 d; 45% slow), with fast rotators predominantly disc-less and slow rotators split between disc-bearing and disc-free, indicating disc braking alone is insufficient. Spot temperatures and coverages are derived from multi-band amplitudes and photometrically calibrated effective temperatures (Eq. 1). The two headline results are: (i) after bias correction via injection-recovery simulations, the intrinsic cold-spot coverage distribution of fast rotators is exponential, consistent with stochastic flux emergence; (ii) slow rotators show a deficit of small cold spots, which the authors attribute to small-spot lifetimes shorter than ~150-200 d, and interpret via turbulent magnetic diffusivity as evidence for rotation-dependent supergranule sizes.","tokens_in":27254,"tokens_out":4607,"duration_ms":163984,"significance":"If the results hold, the paper delivers the largest homogeneous, multi-region census of YSO spot properties to date (144 rotators across 12 fields, a ~2.5x increase over the IC 5070-only analysis), the first empirically bias-corrected intrinsic cold-spot coverage distribution for young stars, and a testable, falsifiable prediction — that small-spot lifetimes depend on rotation — directly checkable with continuous-cadence light curves from TESS or future PLATO data. The injection-recovery quantification of detection bias (Fig. 3) and the explicit robustness tests against Teff systematics and temperature-contrast selection (Sect. 4.2.1) are real methodological strengths that should be credited. The rotation–disc correlation analysis (Table 2, Fig. 4) is a useful, if incremental, confirmation that disc braking alone cannot explain the period bimodality.","major_comments":[{"comment":"The 150-200 d lifetime threshold that anchors the paper's headline interpretation is asserted, not measured. The injection-recovery simulations (Sect. 3.1, Fig. 3) injected strictly coherent sinusoids, so they characterize detectability only on the (period, amplitude) plane and are blind by construction to spot evolution. A spot living ~100 d on a P=15 d star produces a phase-unstable signal over ~9-12 cycles that will fail the multi-method, multi-filter clustering cuts of Sect. 2.3 far more readily than the same lifetime on a P=2 d star (~90 cycles). Detection probability is therefore a function of (P, amplitude, lifetime), and the effective lifetime threshold is period-dependent. This is not a cosmetic issue: if small spots have shorter lifetimes everywhere (consistent with the size-lifetime relation the authors themselves cite, Giles et al. 2017), then part of the small-spot deficit a","section":"Sect. 4.2.2 (also Sect. 3.1, Fig. 3)"},{"comment":"The exponential characterization of the intrinsic cold-spot coverage distribution for fast rotators is a central claim, but no fit parameters, parameter uncertainties, or goodness-of-fit statistic are reported, and the rejection of log-normal and power-law forms is stated without any quantitative model comparison (chi-square, KS, AIC, or similar). Please report the exponential scale and its uncertainty, the number of objects (or slices) entering the fit, the binning sensitivity, and a statistical comparison of the three tested functional forms.","section":"Sect. 4.2.3, Fig. 9"},{"comment":"The sample comprises 620 amplitude sets from 144 objects, with 6-month slices generated every 3 months, so adjacent slices share half their data and multiple slices per object trace the same long-lived spot. The coverage distributions (Figs. 7-9) and the p<6e-5 Monte-Carlo significance (Sect. 4.2) appear to treat slices as independent. Please state explicitly how slices vs. objects were counted, and demonstrate (e.g., by re-computing the significance at the object level or with one randomly selected slice per object) that the period-coverage trend and the exponential fit are not inflated by the correlated repeated measurements.","section":"Sect. 2.2, Sect. 4.2"},{"comment":"The Teff calibration of Eq. (1) never states the actual RMS of the fit in Kelvin; only that the worst colour combination is 10% higher than the best. Since the entire spot-fitting step hinges on Teff, the numerical RMS (and its dependence on spectral type within the ~200-star NAP calibration sample) should be given. Relatedly, the calibration absorbs the NAP region's mean reddening into the fit, but the target regions span a range of extinctions; while the Herbert et al. (2023) +-400 K robustness argument mitigates this, a sentence quantifying the differential-extinction risk across regions (e.g., comparing fitted Teff for the IC 5070 subsample against the Fang et al. 2020 spectroscopic values, which Fig. 6 implies exists) would strengthen Sect. 2.4.","section":"Sect. 2.4, Eq. (1)"}],"minor_comments":[{"comment":"The y-axis label reads 'Fraction of Objcts' (typo). The same figure would benefit from stating N and the fit function explicitly in the caption.","section":"Fig. 9"},{"comment":"The bias-correction grid spans V-band amplitudes 0.02-0.16 mag. Please state the mapping from (coverage, temperature contrast) to V amplitude assumed when applying the Fig. 3 correction to the coverage distribution, and confirm the grid covers the full observed coverage range (up to ~0.45 in Fig. 8). Also clarify what contrast is assumed in quoting the detection limit as a coverage of 0.05.","section":"Sect. 3.1 / Sect. 4.2.3"},{"comment":"Several entries are surprising and warrant a brief explanation: NGC 2264 contributes 285 light curves but only 1 usable slice, and Berkeley 86 contributes 229 light curves and 25 slices with zero periodic detections. A sentence on why these fields yield so few qualifying slices/detections (cadence, crowding, extinction) would help the reader assess cross-region selection effects.","section":"Table 1"},{"comment":"The claim of the 'first empirical constraint on the intrinsic spot coverage distribution of YSOs' should be tempered or more carefully situated relative to cited prior work on spot coverage evolution (e.g., Morris et al. 2020) and the authors' own IC 5070 analysis (Herbert et al. 2024), of which this is an extension by a factor ~2.5.","section":"Sect. 5 (Conclusion)"},{"comment":"The supergranule-size estimates are read off 'their Fig. 1' (Bradshaw & Hartigan 2014). Please describe this procedure explicitly (which curves, what assumed stellar parameters), since the factor-of-five supergranule scaling is a notable secondary claim.","section":"Sect. 4.2.2"},{"comment":"Notation G_G for the absolute G magnitude in Eq. (1) is non-standard; consider M_G. The HS:CS quality-indicator thresholds used to flag degenerate solutions are referenced to Herbert et al. (2024) but never restated; a one-line reminder of the cut values would make Fig. 7 self-contained.","section":"Sect. 2.4, Eq. (1); Fig. 7"}],"recommendation":"major_revision","confidential_remarks":"The methodology leans heavily on the authors' previous HOYS papers (Froebrich et al. 2021; Herbert et al. 2023, 2024), which is efficient and properly disclosed, but referees cannot fully audit the periodogram combination or spot-fitting degeneracies from this manuscript alone. The main revision needed is confined to Sect. 4.2.2 and the framing of the lifetime result; since \"Spot Lifetimes\" is in the title, the editor may wish to ensure this point is genuinely resolved rather than merely reworded."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new empirical pieces are real: after bias correction the cold-spot coverage distribution for the fast rotators (P < 5.5 d) is exponential, and Fig. 8 shows a clear deficit of small cold spots among the slow rotators. Both survive the Teff and contrast checks they ran. That is the first homogeneous multi-region constraint of its kind and is worth having.\n\nThey did the obvious things carefully. Injection-recovery (Fig. 3) shows detection rate is almost period-independent for coherent signals; the Monte-Carlo fluke probability is tiny; systematic Teff shifts move spot temperature by roughly half the Teff error while leaving coverage almost untouched. The disc-braking discussion is measured: they note the usual exceptions (slow disc-less, fast disc-bearing) without over-claiming. Methods are the same pipelines they already published for IC 5070, now applied to a ~2.5\times larger sample, so the work is reproducible in principle.\n\nThe soft spot is exactly the one the stress-test flags. The 150–200 d lifetime threshold is asserted from the 6-month slice length, not measured. All injections were infinitely coherent sinusoids, so the pipeline’s actual sensitivity as a joint function of period, coverage and finite lifetime was never mapped. A 100 d spot on a 15 d rotator produces far fewer stable cycles than the same spot on a 2 d rotator and will fail the multi-method clustering cuts more readily. The clean fast/slow dichotomy could therefore be quantitatively distorted even if the qualitative trend is physical. They do consider reduced contrast and geometry and rule them out reasonably, but residual selection remains under-tested. Minor additional note: the exponential fit treats the 620 slices as independent when many come from the same stars (IC 5070 alone is ~1/3), so the formal goodness-of-fit is a bit optimistic.\n\nNone of that sinks the paper. The empirical distributions are still the best we have for YSOs in this mass and age range, and the lifetime interpretation is offered as the most plausible reading rather than a hard measurement. Anyone working on pre-MS angular-momentum or magnetic evolution will want the numbers. I would send it to referees; they will ask for a short lifetime-injection test or a clearer caveat, which is fair and fixable.","headline":"Solid multi-region HOYS expansion that delivers a usable exponential coverage distribution for fast rotators; the lifetime claim for the slow-rotator deficit is plausible but uncalibrated because the injection tests never varied coherence time.","tokens_in":27947,"tokens_out":574,"would_cite":true,"duration_ms":11879,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Small cold spots on slowly rotating young stars live less than six months; on fast rotators they follow an exponential coverage distribution set by stochastic magnetic flux.","keywords":["young stellar objects","starspots","stellar rotation","spot lifetimes","spot coverage","disc braking","magnetic activity","T Tauri stars"],"falsifier":"High-cadence photometry that recovers the same slow-rotator sample on timescales of weeks instead of months and either detects the missing small-coverage spots or still fails to find them.","tokens_in":27405,"feed_emoji":"⭐","tokens_out":717,"duration_ms":14760,"temperature":0.7,"pith_summary":"Young stars spin at two distinct rates even at roughly one million years of age. This paper measures multi-band light curves of 144 such stars across many star-forming regions and converts the amplitudes into cold-spot coverages and temperatures. Once detection biases are removed, the coverage distribution for fast rotators (periods under 5.5 days) is exponential: tiny spots are far more common than large ones, exactly as expected if magnetic flux emerges stochastically. Slow rotators, by contrast, almost never show those small spots. The authors argue the missing spots are not absent or too cool; they simply decay faster than the six-month windows used to find periods. The result is direct evidence that the surface magnetic pattern of a young star depends on how fast it is spinning, giving models of early angular-momentum loss a new observational target.","feed_headline":"Small spots die fast on slow-spinning young stars","feed_subtitle":"Fast rotators show an exponential coverage law; slow ones lose spots in under six months","key_machinery":"Six-month light-curve slices analysed for multi-band peak-to-peak amplitudes, converted to spot coverage and temperature contrast; injection-recovery simulations then correct the observed coverage histogram, revealing the exponential form for fast rotators and the lifetime cut-off for slow rotators.","core_discovery":"After bias correction the intrinsic cold-spot coverage distribution of fast-rotating young stars is exponential, while slow rotators show a clear deficit of small-coverage cold spots that is best explained by lifetimes shorter than roughly 150–200 days.","pith_inferences":["If supergranule size really tracks rotation, the same lifetime–coverage relation should appear in older, fully convective M dwarfs once comparable multi-band baselines exist.","The exponential coverage law supplies a ready prior for population synthesis of photometric jitter in young exoplanet-host candidates.","Short-lived small spots on slow rotators may still modulate X-ray or UV emission on weekly timescales even when optical periods vanish."],"forward_implications":["Models of YSO magnetic activity must produce an exponential coverage distribution once spots live longer than a few months.","Supergranule (or equivalent diffusion) scales are predicted to shrink by a factor of several when a star spins up after disc dispersal.","Disc braking alone cannot set the observed period bimodality; an additional rotation-dependent magnetic process is required.","Warm/hot-spot solutions separate into two populations (small accretion footprints and larger plage-like features) that can now be tested with phase-curve modelling."],"fun_headline_variants":["Small cold spots fade faster on slow-spinning YSOs","Fast rotators show exponential spot coverage; slow lack small spots","Slow young stars lose small spots in under 200 days","Spot lifetimes drop sharply for slow-rotating young stars","Bias-corrected coverages: exponential in fast YSOs, truncated in slow"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The absence of small cold spots on slow rotators is caused by lifetimes shorter than the six-month detection window rather than by cooler contrasts, different geometry, or residual selection effects the recovery tests missed.","fun_headline_variants_meta":{"raw":{"variants":["Small cold spots fade faster on slow-spinning YSOs","Fast rotators show exponential spot coverage; slow lack small spots","Slow young stars lose small spots in under 200 days","Spot lifetimes drop sharply for slow-rotating young stars","Bias-corrected coverages: exponential in fast YSOs, truncated in slow"]},"model":"grok-4.5","effort":"low","cost_usd":0.004092,"raw_usage":{"total_tokens":1273,"prompt_tokens":832,"num_sources_used":0,"completion_tokens":70,"cost_in_usd_ticks":40924000,"prompt_tokens_details":{"text_tokens":832,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":371,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":832,"tokens_out":70,"duration_ms":6755,"temperature":1.0,"reasoning_tokens":371,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T14:17:24.533130+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"High-cadence photometry that recovers the same slow-rotator sample on timescales of weeks instead of months and either detects the missing small-coverage spots or still fails to find them.","supporting_citations":[],"review_version":1}