{"id":"c7b2f06d-fe74-4e72-85a2-44661d6c61c9","arxiv_id":"2506.09826","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Veiling in classical T Tauri stars is wavelength-dependent, showing rises in the blue, mid-optical, and red, and correlates with Hα emission strength.","lead":"This paper measures how much extra light from accretion, called veiling, changes across the optical spectrum for over 1,800 young stars. It finds that veiling is not constant with wavelength and that stronger hydrogen emission corresponds to stronger veiling, which can improve how astronomers measure accretion.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The WTTS zero-point subtraction assumes WTTS veiling is purely systematic and identical for CTTSs and WTTSs; if spot properties or template mismatch differ with accretion activity, the corrected veiling and its wavelength features are biased.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing concern: the WTTS subtraction assumes no class-dependent systematic difference between CTTSs and WTTSs. I agree with the reader that this assumption is untested and central. The manuscript's own Section 3.3 admits that spots contribute to WTTS veiling, and there is no evidence that spot properties are identical in CTTSs. The injection-recovery test in Figure 7 is a useful validation of the fitting machinery, but it only proves that known veiling can be recovered on WTTS spectra; it does not prove that the zero-point is transferable. I also share the reader's concern about the absence of propagated uncertainties, which makes it impossible to assess whether the reported features are statistically significant. However, I do not see an internal inconsistency in the Teff discussion that would rise to rejection: the claim that veiling shows no clear Teff dependence is consistent with the figures, and the Hα-eqw correlation is a separate claim. I would keep the verdict CONDITIONAL because the paper is novel and useful but requires a validation step before its empirical veiling profiles should be used as a reference. The strongest_claim is plausible but not yet secured; the WTTS subtraction is the insecure point.","tokens_in":14819,"tokens_out":1583,"duration_ms":19192,"concrete_test":"Take a sample of CTTSs with independently measured veiling from higher-resolution or multi-epoch data (e.g., the ULLYSES/CVSO sources already cited) and run the same BOSS pipeline on contemporaneous or appropriately matched spectra; compare the WTTS-subtracted veiling to the independent measurements. If the corrected veiling matches independent values in shape and magnitude, the WTTS zero-point assumption is validated. Additionally, recompute the corrected veiling using field-star subtraction instead of WTTS subtraction for stars in the same Teff range; if the three-feature profile changes significantly, the result is not robust to the choice of control sample.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim—empirical veiling in CTTSs has a specific wavelength dependence with features near 3700 Å, 6000–7000 Å, and 8500–10000 Å—rests on subtracting the average WTTS veiling from CTTS veiling (Section 3.3). This assumes that WTTS veiling is entirely systematic (spots, template mismatch, continuum errors) and that this systematic is identical for CTTSs and WTTSs of the same Teff. If CTTSs have different spot filling factors or spot temperature contrasts, or if the PHOENIX template mismatch depends on accretion-related surface inhomogeneities, the subtraction will remove real accretion veiling and bias the wavelength dependence. The manuscript itself attributes WTTS veiling partly to spots (Section 3.3), so this is not a hypothetical: spot properties are physically expected to differ between accreting and non-accreting stars. The injection-recovery test (Figure 7) only shows that known veiling added to WTTS spectra is recovered; it does not test whether the WTTS zero-point applies to CTTSs. Furthermore, uncertainties on the corrected veiling are not propagated, so the significance of the three reported features cannot be assessed. This assumption is load-bearing because every result (Figures 6, 9, 10) uses the corrected veiling.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents empirical veiling measurements for ~1800 classical T Tauri stars (CTTSs) and ~35,000 weak-line T Tauri stars (WTTSs) using SDSS-V/BOSS spectra. The authors fit PHOENIX synthetic spectra to extinction-corrected spectra in 1000 Å-wide bins stepped by 100 Å to derive wavelength-dependent veiling, then subtract the average WTTS veiling per Teff bin as a systematic zero-point correction. They report three characteristic features in the corrected CTTS veiling (a blue rise from 3700 Å, a mid-optical peak near 6000–7000 Å, and a red rise from 8500 Å), a correlation between veiling and Hα equivalent width, no strong Teff dependence, and no clear age trend. The paper also compares two sources with previous studies (Ingleby et al. 2013; ULLYSES) and tests the pipeline by injecting artificial veiling into WTTS spectra.","tokens_in":15073,"tokens_out":3365,"duration_ms":36571,"significance":"If the measured profiles are accurate, this is a valuable large-sample empirical characterization of wavelength-dependent veiling in CTTSs, directly relevant to accretion shock models. The injection-recovery test (Figure 7) provides a useful sanity check on the pipeline, and the homogeneous treatment of a large sample is a strength. However, the central result depends on the assumption that WTTS veiling is purely systematic and transferable to CTTSs, which is not demonstrated; and the significance of the reported spectral features is not quantifiable from the presented figures because uncertainties are not propagated into the mean veiling profiles.","major_comments":[{"comment":"The zero-point subtraction assumes that the WTTS veiling is entirely systematic (spots, template mismatch, continuum errors) and identical for CTTSs and WTTSs at the same Teff. The manuscript itself attributes WTTS veiling partly to spots, and later attributes the red-end rise in CTTSs to photospheric spots (Section 5.1). If spot properties or template mismatch differ between accreting and non-accreting stars, the subtraction will bias the corrected veiling and the inferred wavelength dependence. The injection-recovery test only shows that known continuum added to WTTS spectra is recovered; it does not test whether the WTTS zero-point applies to CTTSs. Please provide additional validation that the systematic component is the same for both populations, or propagate this uncertainty into the corrected veiling values and the final results.","section":"Section 3.3, Figure 7"},{"comment":"The average veiling profiles are plotted without any propagated uncertainties or confidence bands, despite the text stating that weighted means were computed using the veiling uncertainties as weights. Moreover, the wavelength bins are 1000 Å wide with 100 Å steps, so adjacent measurements are strongly correlated. Without error bars and an account of bin-to-bin correlation, the significance of the 'three distinct peaks' (the blue rise, mid-optical peak, and red rise) cannot be assessed. Please include uncertainties on the mean profiles and quantify the correlation, or state the effective number of independent wavelength points.","section":"Figures 9 and 10"},{"comment":"The Results and Discussion (Section 5.3) state that veiling shows no clear dependence on Teff, yet the Conclusions say 'we found veiling to be related to both Hα eqw and Teff' and 'some of the features in veiling were found to be more prominent in some Teff ranges.' These statements are inconsistent. Please reconcile the summary with the analysis and clarify which Teff trends are actually supported by the data shown in Figures 9 and 10.","section":"Section 5.3 and Conclusions"}],"minor_comments":[{"comment":"Equation (1) appears malformed: 'Fsource,λ = Fphot,λ + rλ (1 + rλ)' is missing a fraction or operator, making it impossible to read the intended veiling definition. Please check the typesetting.","section":"Equation (1)"},{"comment":"The description of the wavelength binning says '1000 Å wide with advancing in steps of 100 Å'; please state explicitly how many independent wavelength measurements result from this heavily overlapping binning and why this choice was made rather than non-overlapping bins.","section":"Section 3.2"},{"comment":"Several figure captions contain garbled symbols (e.g., 'H® eqw', 'Te®', 'Wavelength(ºA)') and should be corrected in the final version.","section":"Figure captions"},{"comment":"The statement that 'the peak of the veiling in mid-optical region shifts towards redder wavelengths for lower eqw values' appears only in the Conclusions; if this is a result, it should be shown in a figure or quantified in Section 4.2.","section":"Conclusions, Section 4.2"},{"comment":"The comparison with 'Pittman et al. 2025, in prep' is cited as a private communication; please update this to a public reference if available, or remove it if the comparison cannot be verified by readers.","section":"Section 4.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for a stellar astrophysics journal. The main concern is the untested WTTS subtraction assumption; if the authors can provide further validation (e.g., comparing CTTSs with independent veiling indicators or testing whether spot properties differ systematically between CTTSs and WTTSs), the paper could be publishable. The lack of uncertainty propagation on the mean veiling profiles is also a serious presentation issue that undermines the headline claims about spectral features."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is the first large-sample empirical mapping of veiling vs wavelength in CTTSs, and it's worth engaging with, but the central calibration step—subtracting mean WTTS veiling—carries more weight than the paper acknowledges.\n\nWhat's new: Basri & Batalha and Ingleby et al. did this for a handful of objects; this applies it to ~1800 CTTSs with homogeneous BOSS spectra. The pipeline is clearly described, the injection-recovery test is a real check, and the comparison to CVSO 58/90 shows they engage honestly with prior data. The correlation with Halpha eqw is not circular, since veiling is measured from continuum shape independently of Halpha. If the corrected veiling profiles hold up, the three spectral features are a useful empirical constraint on accretion shock geometries. That is a legitimately new result.\n\nThe soft spots are in the calibration, not in the execution. The WTTS subtraction assumes the systematic veiling is identical for CTTSs and WTTSs at the same Teff, and the paper itself attributes WTTS veiling partly to spots. CTTSs are generally more active and have different spot properties, so this zero-point could absorb real accretion veiling. The injection-recovery test only shows you can recover added artificial veiling in WTTS spectra; it does not test whether the WTTS zero-point transfers. This is load-bearing for all corrected results. Relatedly, figures 9 and 10 do not propagate uncertainties onto the mean veiling profiles, so the significance of the three features is not actually demonstrated. The overlapping wavelength bins (1000 Å wide, 100 Å step) are a presentation concern rather than a fatal one, but they do inflate apparent correlations between adjacent points.\n\nOne smaller internal inconsistency: the conclusions state veiling is related to Teff, while Section 5.3 says no clear dependence across the three Teff ranges. The plots seem to show similar shapes with different normalizations; the paper should settle on what it actually claims.\n\nBottom line: a serious referee should see this. The science is clear and the result would be a useful reference, but the WTTS subtraction assumption and the missing error propagation need to be fixed—or at least argued for much more carefully—before the wavelength features are treated as established. For a reading group it is a good case study in large-sample calibration choices. I'd accept it for review and push for revision.","headline":"First large-sample empirical veiling–wavelength profiles for CTTSs, worth reviewing, but the WTTS zero-point subtraction is load-bearing and needs stronger justification.","tokens_in":15635,"tokens_out":1648,"would_cite":true,"duration_ms":18820,"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":"Empirical veiling in Classical T Tauri stars is wavelength-dependent, with a UV rise, a mid-optical peak, and a red-end rise, and it correlates with H-alpha emission strength.","keywords":["stellar accretion","veiling","Classical T Tauri stars","wavelength dependence","H-alpha emission","pre-main-sequence stars","PHOENIX synthetic spectra","SDSS-V BOSS"],"falsifier":"Measure veiling for the same Classical T Tauri stars on the same nights with high-resolution optical spectra fitted by a full model that includes spots, and compare the inferred $r_{\\lambda}$ profile with the BOSS-derived profile; if the ultraviolet rise, the 6000–7000 Å peak, or the red-end rise does not reproduce, the empirical profile is dominated by template or spot systematics.","tokens_in":14636,"feed_emoji":"🔭","tokens_out":8894,"duration_ms":88821,"temperature":0.7,"pith_summary":"This paper tries to establish that veiling—the excess continuum emission that accretion adds to a young star's photospheric spectrum—is not a single constant but a structured function of wavelength. Analyzing roughly 1,800 Classical T Tauri stars observed by the BOSS spectrograph, the authors measure empirical veiling profiles and find three recurring features: a rise toward the ultraviolet from 3700 Å, a mid-optical peak near 6000–7000 Å, and a red-end rise beyond 8500 Å. They also claim that veiling tracks H-alpha equivalent width, with stronger H-alpha emission corresponding to higher veiling, while age and effective temperature do not strongly organize accretion activity. If the claim holds, veiling becomes a practical diagnostic of where and how accretion energy is released, and the common assumption of wavelength-independent veiling in accretion studies should be abandoned.","feed_headline":"Veiling in young stars is not flat: three features emerge","feed_subtitle":"Analysis of 1,800+ Classical T Tauri stars links a UV rise, mid-optical peak, and red-end excess to H-alpha strength.","key_machinery":"Veiling is measured empirically by comparing each extinction-corrected BOSS spectrum with a PHOENIX synthetic template: the spectrum is divided into 1000 Å wavelength bins stepped by 100 Å, and a least-squares fit solves for the additive excess $r_{\\lambda}$ that best reproduces the observed line depths under the definition $F_{\\mathrm{source},\\lambda} = F_{\\mathrm{phot},\\lambda} + r_{\\lambda}/(1+r_{\\lambda})$. Because weak-line T Tauri stars should have no accretion veiling, their measured veiling is treated as a systematic zero point, binned in effective temperature, and subtracted from the Classical T Tauri veiling. An injection-recovery test, in which continuum of known strength is added to WTTS spectra and recovered through the same pipeline, is used to validate that the measurement procedure returns the veiling that was put in.","core_discovery":"The central claim is that empirical veiling in Classical T Tauri stars varies with wavelength in a reproducible three-feature pattern across 3600–10400 Å, visible in the sample-averaged profiles and stable across temperature bins. The ultraviolet rise is attributed to small, very hot spots where accretion columns shock onto the stellar surface; the stronger mid-optical peak near 6000–7000 Å is attributed to larger, cooler emitting regions; and the red-end rise from 8500 to 10000 Å is tentatively attributed to photospheric spots or noise rather than to accretion. The paper further claims that veiling increases monotonically with H-alpha emission strength, that this trend holds across effective temperature ranges, and that veiling shows no clear dependence on age because disk dissipation is stochastic.","pith_inferences":["Editorial inference: if the measured profile is real, fitting Classical T Tauri spectra with constant veiling will systematically bias effective temperature and surface gravity estimates; a two-component hot-plus-cool spot veiling law would likely improve parameter recovery.","Editorial inference: the red-end veiling rise implies that mass accretion rates derived from red or near-infrared veiling alone could be overestimated, a prediction that independent near-infrared veiling measurements could test.","Editorial inference: a direct test of the WTTS subtraction assumption would be to observe spotted weak-line T Tauri stars over a full rotation; if their systematic veiling changes with spot coverage, the subtraction can distort the accreting-star profiles."],"forward_implications":["Veiling must be treated as wavelength dependent in Classical T Tauri star studies; a single scaling factor hides the three-feature structure and will bias any derived photospheric parameters.","H-alpha equivalent width is a working proxy for accretion strength: more negative equivalent widths correspond to higher veiling across the sampled effective temperature ranges.","The mid-optical veiling peak indicates emission from large, relatively cool regions on the stellar surface, while the ultraviolet rise indicates compact, hot accretion shock spots nested within them.","The red-end veiling rise is most likely not accretion related and should not be used as an accretion tracer in this wavelength range.","Age alone does not predict veiling; disk dissipation is probabilistic, so some older stars keep strong accretion while some young stars are already quiescent."],"supporting_citations":[{"why":"Supplies the theoretical accretion shock models that predict wavelength-dependent veiling and the physical interpretation of the measured profile.","marker":"Calvet & Gullbring 1998"},{"why":"Early empirical veiling measurements across wavelengths; establishes the expectation that veiling varies and provides the comparison baseline.","marker":"Basri & Batalha 1990"},{"why":"Previous source-by-source veiling profiles; CVSO 58 and CVSO 90 from that work are matched and compared with the new measurements.","marker":"Ingleby et al. 2013"},{"why":"PHOENIX synthetic spectral library used as the photospheric template for the veiling fits.","marker":"Husser et al. 2013"},{"why":"The LINEFOREST pipeline that identifies young stars and classifies Classical versus weak-line T Tauri stars, defining the sample.","marker":"Saad et al. 2024"},{"why":"BOSS Net supplies the effective temperatures, surface gravities, and radial velocities used to select the matching synthetic spectra.","marker":"Sizemore et al. 2024"},{"why":"Provides the extinction model used to de-redden each observed spectrum before veiling is measured.","marker":"Gordon 2024"},{"why":"Sagitta neural net produces the individual stellar ages used in the veiling versus age analysis.","marker":"McBride et al. 2021"}],"fun_headline_variants":["Young star veiling shows three wavelength features","Veiling in T Tauri stars: UV rise, optical peak, red excess","Accretion veiling reveals three spectral patterns in young stars","New veiling measurements trace accretion in T Tauri stars","Veiling varies with wavelength in young stars, new study shows"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that veiling measured in weak-line T Tauri stars is purely systematic and identical for Classical and weak-line T Tauri stars of the same effective temperature, so subtracting it removes bias rather than real accretion-related differences.","fun_headline_variants_meta":{"raw":{"variants":["Young star veiling shows three wavelength features","Veiling in T Tauri stars: UV rise, optical peak, red excess","Accretion veiling reveals three spectral patterns in young stars","New veiling measurements trace accretion in T Tauri stars","Veiling varies with wavelength in young stars, new study shows"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000769,"raw_usage":{"total_tokens":3401,"prompt_tokens":932,"completion_tokens":2469,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":548,"completion_tokens_details":{"reasoning_tokens":2384}},"tokens_in":548,"tokens_out":2469,"duration_ms":18101,"temperature":1.0,"reasoning_tokens":2384,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:39:36.771362+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure veiling for the same Classical T Tauri stars on the same nights with high-resolution optical spectra fitted by a full model that includes spots, and compare the inferred $r_{\\lambda}$ profile with the BOSS-derived profile; if the ultraviolet rise, the 6000–7000 Å peak, or the red-end rise does not reproduce, the empirical profile is dominated by template or spot systematics.","supporting_citations":[],"review_version":1}