{"id":"d90df62d-948a-4ac1-99a3-440a8fdc0bdb","arxiv_id":"2505.04699","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"An all-sky Gaia XP-spectra catalogue of 145,975 Hα-emitting young stellar object candidates within 500 pc yields homogeneous accretion luminosities and rates and reveals a spatially dispersed, low-accreting population.","lead":"Using Gaia's low-resolution spectra, the authors identify 145,975 nearby stars with Hα emission, treat them as candidate young stars with accreting disks, and estimate their accretion luminosities and rates. It is the first all-sky, homogeneously processed accretion survey of the local 500 pc volume, and it points to a scattered population of old, low-accreting young stars.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed dispersed low-accreting YSO population rests on the unfiltered Table 2, yet Sect. 2.7 reports 55% of that table lies below the chromospheric emission locus; until that subset is removed or validated, the low-accreting population claim is not established.","rationale":"The reader's weakest-assumption diagnosis is accurate and is the load-bearing issue: the pipeline's conversion of Halpha into accretion luminosity is applied uniformly to the full 145,975-source table, while 55% of those sources fall below the chromospheric emission level by the paper's own flag. This does not invalidate the catalogue construction or the purified sample results, which are independently supported by the comparison with 341 X-Shooter YSOs and by the very low contamination fractions in samples A, B, and C. But it does mean the abstract's headline claim of a large, previously untraced population of dispersed low-accreting YSO candidates is not supported by the evidence as presented for the unfiltered sample. The appropriate remedy is not rejection but a conditional acceptance that requires either a chromospheric cut in the headline population or an explicit demonstration that the low-accreting component survives that cut. Since the reader already arrived at CONDITIONAL, my stress-test pass does not change the verdict.","tokens_in":33253,"tokens_out":3374,"duration_ms":37457,"concrete_test":"Recompute all low-accreting population statistics after applying flag_above_chromospheric_level=True as a hard cut to every source in Table 2, including the 20,274 without XP spectra, then re-derive the counts, spatial distribution, and 'Peter Pan'-type classification of sources with log Lacc < -3.5 or Mdot < 1e-11. If most of the dispersed low-accreting population lies below the chromospheric locus or disappears once M-dwarf contaminants are excluded, the abstract claim should be restricted to the purified subsets A, B, and C.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central physical inference — a large previously untraced population of low-accreting, often dispersed YSO candidates, including 'Peter Pan' analogs — is extracted from the full Table 2 sample selected at pEW < -0.5 nm (Sect. 2.1), a threshold deliberately chosen with a completeness-weighted F-beta score (beta=2) rather than for purity (Appendix B). For every source in this table, Lacc is computed by converting pEW to EWHalpha (Eq. 1), estimating Fcont (Eq. 3), and applying the Alcala et al. (2017) Lacc-LHalpha relation (Eq. 5). That relation is calibrated on classical T Tauri magnetospheric accretors and does not distinguish accretion-powered Halpha from chromospheric Halpha. The paper's own flag_above_chromospheric_level (Sect. 2.7) places 55% of the full table below the Manara et al. (2017a) chromospheric emission locus. For those sources the reported accretion luminosities and rates could be dominated by chromospheric activity, so the abstract's 'large population of low-accreting YSO candidates untraced by previous surveys' is not actually validated for the unfiltered catalogue. The purified samples A/B/C are credible (8.44%, 0.15%, and ~0 below the locus, respectively), but those are not the samples used for the headline low-accretor discovery claim, and the 341-source X-Shooter comparison in Sect. 2.8 tests known, relatively strong accretors rather than the low-accretion regime. The additional 20,274 sources without public XP spectra, filtered only through pEW and a manually tuned Lacc threshold in scanning-law regions (Appendix E), add further contamination risk.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first all-sky, homogeneous derivation of H-alpha-based accretion properties for young stellar objects within 500 pc using Gaia DR3 XP spectra. The authors construct a catalogue of 145,975 H-alpha emitter candidates (Table 2), derive equivalent widths from ESP-ELS pseudo-equivalent widths, continuum fluxes from an RP-magnitude relation, extinction corrections via two methods (GSP-Phot and median GSP-Phot), and convert H-alpha luminosities to accretion luminosities and mass accretion rates using the Alcala et al. (2017) calibration. They define three purer sub-samples: flag_IR (4,208 sources), flag_pEW (6,170 sources), and flag_combined (1,945 sources), and use these to derive Lacc-Lstar and Macc-Mstar relations, a Sco-Cen accretion timescale of 2.7±0.4 Myr, and the claim of a large, spatially dispersed population of low-accreting YSO candidates. The catalogue and methodology are the main deliverables.","tokens_in":33686,"tokens_out":3804,"duration_ms":39773,"significance":"If the catalogue is accepted as reliable, it represents a substantial advance: a homogeneous all-sky local census of YSO accretion, with transparent, reproducible methodology and explicit quality flags. The comparison against 341 X-Shooter YSOs is a genuine strength, as are the machine-readable catalogue and the explicit documentation of selection thresholds. However, the central physical claim—the existence of a large, previously unseen low-accreting population—rests on the unfiltered catalogue, for which the paper itself shows that 55% of sources fall below the chromospheric emission locus. The purified samples are more credible but are not the basis for the headline discovery claim. The small formal errors on the derived power-law slopes also need to be reconciled with the known systematic scatter in the calibration chain.","major_comments":[{"comment":"The abstract's claim of a 'large population of low-accreting YSO candidates untraced by previous surveys' is drawn from the full Table 2, yet Sect. 2.7 reports that 55% of sources in that table lie below the chromospheric emission locus of Manara et al. (2017a). For those sources, the derived Lacc and Macc values could be dominated by chromospheric emission, making the population claim unsubstantiated as presented. The authors should re-derive the low-accretor statistics after applying flag_above_chromospheric_level, or explicitly quantify how many low-accreting sources in Table 2 remain above the chromospheric locus and show that the spatial-dispersion result survives that cut.","section":"Sect. 2.7 and Abstract"},{"comment":"The validation against 341 X-Shooter YSOs tests relatively strong, known accretors and shows agreement only to within an order of magnitude. The low-accretion regime (log Lacc below about -4) is not probed by this comparison, yet it is exactly the regime where the new dispersed population is claimed to reside. The paper should either validate the pipeline on known chromospherically active non-accreting stars (to demonstrate that the method does not overproduce low-Lacc sources) or state explicitly that the low-accretion rates in the catalogue are unvalidated and should be treated with caution until spectroscopic follow-up is available.","section":"Sect. 2.8, Fig. 9"},{"comment":"The quoted uncertainties on the power-law slopes (0.02 for the Lacc-Lstar slope and 0.1 for the Macc-Mstar slope) reflect only formal propagation of random errors. Table 3 shows that switching from GSP-Phot to med-GSP-Phot extinction changes the combined-sample slope from 1.41±0.02 to 1.50±0.02, and switching samples changes it further; the intrinsic scatter in Eq. (1) and Eq. (5) and the choice of extinction treatment are not included in the reported errors. The paper should report a systematic error budget or a combined uncertainty that accounts for these choices, since the claim '1.41±0.02' is likely an underestimate of the true uncertainty.","section":"Sect. 3.2, Table 3"},{"comment":"The full Table 2 includes 20,274 sources without public XP spectra, which receive no linefinder-based M-dwarf filtering (width_lf cut) and are retained on the basis of pEW alone. The pEW<−0.5 nm threshold was chosen with an F_beta score (beta=2) that prioritizes completeness over purity. The paper should state explicitly how many sources in Table 2 lack XP spectra and quantify the impact of their inclusion on the low-accretor statistics, especially for the spatially dispersed population claimed in the abstract.","section":"Sect. 2.1, Appendix E"}],"minor_comments":[{"comment":"The text contains several spacing errors ('di fferent', 'e ffects') and a grammatical error in Sect. 2.8 ('the lower limit to the mass accretion rate that we are sensitive too' should be 'sensitive to'). A careful proofread would improve readability.","section":"Throughout"},{"comment":"The fit of Eq. (1) discards 18 points with pEW > -0.1 nm or EWHalpha > -0.1 nm; the paper should justify more explicitly that these are non-emitting or peculiar sources and confirm that the retained 120 sources span the parameter space of the final catalogue, including the low-pEW regime.","section":"Sect. 2.2"},{"comment":"Several Sco-Cen clusters have only upper limits or no median accretion values (e.g., clusters 11, 12, 24, 25, 27, 30), and the table is not self-explanatory about why. A short note in the caption clarifying that medians require at least four accretors would help.","section":"Sect. 4.1, Table 4"},{"comment":"The F_beta-score definition and the choice of beta=2 are described, but the paper does not state the numerical values of the accuracy, completeness, and F_beta for the chosen thresholds (pEW<−0.5 and <−1.0 nm). Reporting those numbers would allow readers to judge the purity/completeness tradeoff quantitatively.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"The manuscript presents a potentially valuable catalogue and a generally transparent pipeline. The main reservation is the mismatch between the headline low-accreting population claim and the paper's own 55% flag on the unfiltered table. This is fixable within the scope of a revision, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The catalogue is the real product here, and it is a good one. First all-sky, distance-limited (500 pc), homogeneous derivation of Lacc and Macc from Gaia XP spectra, covering 145,975 candidate Hα emitters. The pipeline is transparent, the flags are sensible, and the validation against 341 X-Shooter YSOs is honest: order-of-magnitude agreement for GSP-Phot extinction, with the med-GSP-Phot values deliberately flagged as lower limits. That alone makes this a subfield-level resource that roughly triples the number of YSOs with accretion estimates.\n\nWhat is genuinely new: the XP-spectra-based method itself, the catalogue, and the claim of a dispersed low-accreting population that previous surveys missed. The scaling relations (Lacc–L* slope 1.41, Macc–M* slope 2.4) agree with the range of published values, and the Sco-Cen accretion decay analysis is a sensible application.\n\nThe soft spots are real but mostly in the packaging. The abstract's \"large population of low-accreting YSO candidates untraced by previous surveys\" is drawn from the full Table 2, yet Sect. 2.7 reports 55% of that table lies below the chromospheric emission locus. The purified samples A/B/C are much cleaner (8.44%, 0.15%, and ~0 below the locus), and the paper does show a low-accreting population persists in those samples, so the scientific claim is not dead. But the abstract should not rest on the unfiltered table. Second, the headline τ=2.7 Myr comes from an exponential fit that the paper itself says fits worse than a power law; reporting both is fine, but the abstract leads with the worse-fitting model. Third, the 70%/2.8% power-law extrapolations in the abstract carry no error bars. These are wording and emphasis problems, not fatal flaws. A more substantive concern is that the 341-object validation mostly probes relatively strong accretors, so the low-accretion regime is less well calibrated.\n\nCitation pattern is solid — Alcalá 2017, Manara 2023, and the Sco-Cen cluster papers are all handled properly. The catalogue will be used widely, especially after DR4.\n\nWho is this for? Anyone working on YSO accretion, disc evolution, or the local star-forming population. It deserves a serious referee. Send it out, but ask the authors to either move the low-accreting population claim to the filtered samples or gate it on the chromospheric flag, and to present the power-law fit as the primary timescale or at least give both equal weight in the abstract.","headline":"A genuinely useful all-sky YSO accretion catalogue from Gaia XP spectra, with an abstract that oversells the low-accreting population by leaning on the unfiltered table.","tokens_in":34284,"tokens_out":2525,"would_cite":true,"duration_ms":27005,"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 argues that Gaia DR3 XP spectra can yield reliable accretion rates for 145,975 young stellar objects within 500 pc, revealing a dispersed low-accreting population.","keywords":["accretion rates","young stellar objects","Gaia DR3 XP spectra","H-alpha emission","protoplanetary disks","Sco-Cen complex","accretion timescale","low-accreting YSO population"],"falsifier":"Take a random subset of the 145,975 candidates that fall below the chromospheric emission locus (55% of the table) and observe them with medium-resolution spectroscopy (resolving power ~10,000) measuring both UV excess and $\\mathrm{H}\\alpha$. If nearly all show chromospheric $\\mathrm{H}\\alpha$ with no UV continuum excess and no veiling, the claimed large dispersed low-accreting population is mostly a selection artifact; if a substantial fraction show genuine UV excess, the population is real.","tokens_in":33019,"feed_emoji":"🌟","tokens_out":10404,"duration_ms":83906,"temperature":0.7,"pith_summary":"This paper tries to establish that the low-resolution Gaia DR3 XP spectra can be used as a homogeneous, all-sky probe of accretion onto young stellar objects within 500 pc. Using the $\\mathrm{H}\\alpha$ pseudo-equivalent width and the line information extracted from the spectra, it derives accretion luminosities and mass accretion rates for 145,975 candidate $\\mathrm{H}\\alpha$ emitters, roughly an order of magnitude more sources than previous targeted surveys. A sympathetic reader would care because the catalogue turns a set of individually studied star-forming regions into a single statistically tractable sample, and it exposes a population of faint, dispersed accretors that earlier surveys missed. The paper also argues that this population's existence changes how accretion timescales, cluster environments, and the so-called 'Peter Pan' discs are understood.","feed_headline":"Gaia spectra map accretion for 145,975 nearby young stars","feed_subtitle":"First all-sky homogeneous census finds a dispersed low-accreting population and a 2.7 Myr accretion timescale.","key_machinery":"The load-bearing object is the $\\mathrm{H}\\alpha$ pseudo-equivalent width (pEW) from the Gaia DR3 ESP-ELS module, together with line measurements from the linefinder tool on the XP spectra. The pEW is calibrated against medium-resolution equivalent widths, converted to a line flux using a continuum flux derived from the Gaia $RP$ magnitude, corrected for extinction in three ways (none, GSP-Phot, and med-GSP-Phot), and finally turned into $L_\\mathrm{acc}$ and $\\dot{M}_\\mathrm{acc}$ through the empirical relation $\\log(L_\\mathrm{acc}/L_\\odot) = 1.13 + 1.74\\,\\log(L_{\\mathrm{H}\\alpha}/L_\\odot)$ and the standard magnetospheric-accretion formula $\\dot{M}_\\mathrm{acc} \\approx 1.25\\,L_\\mathrm{acc} R_\\star/(G M_\\star)$. Linefinder line-width and line-depth cuts ($\\mathrm{width} < 25$ nm, $\\mathrm{depth} > 10^{-17}$ W/nm/m$^2$) remove M-dwarf TiO contamination, while three quality flags (IR excess, pEW strength, and their combination) define purified subsamples. This machinery is what lets a low-resolution, all-sky survey stand in for targeted high-resolution spectroscopy.","core_discovery":"The central claim is that $\\mathrm{H}\\alpha$ emission measured in Gaia XP spectra can be converted into reliable accretion luminosities and mass accretion rates for essentially all YSOs within 500 pc, not just those with infrared-bright discs. After calibrating the Gaia pseudo-equivalent width to literature equivalent widths, deriving $\\mathrm{H}\\alpha$ line fluxes from a continuum flux–$RP$ relation and three extinction treatments, and converting through the empirical $L_\\mathrm{acc}$–$L_{\\mathrm{H}\\alpha}$ relation, the authors produce a table of 145,975 objects with accretion properties and stellar parameters. From the cleanest subsample (1,945 objects passing both IR-excess and strong-$\\mathrm{H}\\alpha$ cuts) they obtain $L_\\mathrm{acc} \\propto L_\\star^{1.41\\pm0.02}$ and $\\dot{M}_\\mathrm{acc} \\propto M_\\star^{2.4\\pm0.1}$, and from the Sco-Cen clusters an exponential accretion timescale of $2.7\\pm0.4$ Myr with a power-law decay giving 70% accretors at 2 Myr and 2.8% at 10 Myr. They further claim that the newly found population of low-accreting candidates is mostly spatially dispersed, often unconnected to known star-forming regions, and may contain numerous 'Peter Pan' discs.","pith_inferences":["The 55% of sources below the chromospheric emission level are probably not all real accretors; medium-resolution spectroscopy of a random subset of those sources would likely shrink the dispersed low-accreting population, though the purified flag_combined sample suggests a real core remains.","If the 'Peter Pan' interpretation is right, the dispersed low accretors should show old ages, thin-disc signatures, and kinematics unconnected to any association; checking their lithium abundances or space motions would test this.","The same XP-based pipeline could be applied to the epoch-level spectra expected in future data releases, turning accretion variability and burst statistics into a general diagnostic rather than a single-epoch census.","The method's extinction limitations (med-GSP-Phot underestimates accretion luminosities) imply that individual-source accretion rates in the table should be treated as lower limits when only the filled-in extinction is available."],"forward_implications":["Accretion rates can now be assigned to any $\\mathrm{H}\\alpha$-emitting source within 500 pc with public XP spectra, so future studies can build statistically complete YSO samples without new spectroscopic campaigns.","The $L_\\mathrm{acc}$–$L_\\star$ and $\\dot{M}_\\mathrm{acc}$–$M_\\star$ relations are recovered on a homogeneous all-sky sample, with slopes consistent with most literature values, which supports using these relations to interpret larger and fainter samples.","The fraction of accretors in Sco-Cen declines on a timescale of $2.7\\pm0.4$ Myr (exponential) or as a power law with 70% at 2 Myr and 2.8% at 10 Myr, giving a comparative benchmark for disc dispersal models.","A large population of low-accreting, spatially dispersed YSO candidates, many with no clear young association, is now available for follow-up, including candidate 'Peter Pan' discs.","Quality-flag subsamples (4,208 by IR excess, 6,170 by strong $\\mathrm{H}\\alpha$, 1,945 by both) let users balance completeness against contamination for different science cases."],"supporting_citations":[{"why":"Calibrates the $L_\\mathrm{acc}$–$L_{\\mathrm{H}\\alpha}$ relation (A=1.13, B=1.74) used to convert every $\\mathrm{H}\\alpha$ line luminosity into an accretion luminosity.","marker":"Alcalá et al. 2017"},{"why":"Provides the ESP-ELS $\\mathrm{H}\\alpha$ pseudo-equivalent widths in Gaia DR3 that define the input sample.","marker":"Fouesneau et al. 2023"},{"why":"Supplies the geometric distances used to select the <500 pc sample and to compute luminosities.","marker":"Bailer-Jones et al. 2021"},{"why":"The reference set of 341 well-characterised T Tauri stars used to validate the derived accretion luminosities and rates.","marker":"Manara et al. 2023"},{"why":"Supplies the homogeneous Sco-Cen cluster ages used to derive the accretion timescale.","marker":"Ratzenböck et al. 2023a"},{"why":"Supplies the Sco-Cen cluster memberships used to measure the fraction of accretors per cluster.","marker":"Ratzenböck et al. 2023b"},{"why":"Pre-main-sequence evolutionary tracks used to derive stellar masses and luminosities from the dereddened colour-magnitude diagram.","marker":"Baraffe et al. 2015"},{"why":"Provides the previous accretion timescale and fraction-of-accretors method the Sco-Cen analysis is compared with.","marker":"Fedele et al. 2010"},{"why":"Bayesian regression method that accounts for measurement uncertainties in all power-law and decay fits.","marker":"Kelly 2007"},{"why":"Linefinder tool on XP spectra used to measure $\\mathrm{H}\\alpha$ line widths, depths, and continuum fluxes and to remove M-dwarf contaminants.","marker":"Weiler et al. 2023"}],"fun_headline_variants":["Gaia XP spectra unveil accretion for 145,975 young stars","All-sky Gaia census finds dispersed low-accreting YSOs","Accretion timescale 2.7 Myr from Gaia's all-sky YSO map","First homogeneous accretion rates for all nearby young stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that $\\mathrm{H}\\alpha$ pseudo-equivalent width traces magnetospheric accretion through the literature $L_\\mathrm{acc}$–$L_{\\mathrm{H}\\alpha}$ relation for every source in the table, even though the paper finds 55% of sources fall below the chromospheric emission level.","fun_headline_variants_meta":{"raw":{"variants":["Gaia XP spectra unveil accretion for 145,975 young stars","All-sky Gaia census finds dispersed low-accreting YSOs","Accretion timescale 2.7 Myr from Gaia's all-sky YSO map","First homogeneous accretion rates for all nearby young stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000246,"raw_usage":{"total_tokens":1690,"prompt_tokens":1247,"completion_tokens":443,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":863,"completion_tokens_details":{"reasoning_tokens":363}},"tokens_in":863,"tokens_out":443,"duration_ms":4228,"temperature":1.0,"reasoning_tokens":363,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:24:08.592983+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a random subset of the 145,975 candidates that fall below the chromospheric emission locus (55% of the table) and observe them with medium-resolution spectroscopy (resolving power ~10,000) measuring both UV excess and $\\mathrm{H}\\alpha$. If nearly all show chromospheric $\\mathrm{H}\\alpha$ with no UV continuum excess and no veiling, the claimed large dispersed low-accreting population is mostly a selection artifact; if a substantial fraction show genuine UV excess, the population is real.","supporting_citations":[{"cited_title":"2023, , 674, A28","cited_arxiv_id":null,"evidence_quote":"Provides the ESP-ELS $\\mathrm{H}\\alpha$ pseudo-equivalent widths in Gaia DR3 that define the input sample."}],"review_version":1}