{"id":"25e574bf-6aa4-4d3a-81a6-afdcee266147","arxiv_id":"2411.13247","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Using far-UV H2 fluorescent line ratios from HST ULLYSES spectra, the authors derive interstellar AV and RV for 34 classical T Tauri stars and find broad agreement with optical extinction measurements, with caveats from line self-absorption.","lead":"Astronomers measured how much dust blocks the ultraviolet light of 34 young stars by comparing far-ultraviolet hydrogen molecule emission lines with theoretical line ratios. The method, which assumes the lines are not absorbed by the emitting gas itself, mostly agrees with earlier optical measurements but can be biased by self-absorption.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim rests on an optically thin assumption that is already contradicted by the TW Hya control, and the Table 2 tau lower limits do not rule out optical depth in the 'safe' progressions.","rationale":"I read the paper in good faith: it transparently presents a new method to derive AV/RV from H2 fluorescent line ratios, including a discussion of self-absorption and a control object. The simulations (Sect. 3.1, Fig. 1) show that the fitting pipeline can recover injected extinction values on optically thin synthetic data, which supports internal consistency but cannot validate the central physical assumption. The load-bearing condition for the paper's conclusion is that the observed H2 line ratios are set only by branching ratios and dust, i.e., that the lines are optically thin. This condition is not met for the excluded [0,1] and [0,2] progressions, and the paper's Table 2 tau values are lower limits that may be orders of magnitude too small. The strongest evidence that even the 'safe' progressions are affected is TW Hya: with a well-established AV of 0.0 mag, the method returns AV=1.2±0.3 at RV=3.1. The paper's own language ('our procedure failed' for TW Hya) acknowledges this failure. Since TW Hya is the cleanest control, the anomaly cannot be dismissed, and it raises the possibility that AV values for the other stars are systematically biased high by self-absorption that mimics dust. The post hoc exclusion of the [0,1] progression and the lack of a published line-flux table or code make the level of agreement with literature (Fig. 4) harder to assess independently. A concrete radiative-transfer-based correction for the optically thin assumption would settle whether the central claim holds. I therefore agree with the reader's CONDITIONAL verdict and recommend no change.","tokens_in":19605,"tokens_out":7398,"duration_ms":78149,"concrete_test":"Re-derive Method I for TW Hya after correcting the [1,7] and [1,4] line fluxes for self-absorption using the McJunkin et al. (2016) radiative-transfer column densities (log N_tot ~ 19, T_exc = 2500 K) instead of the lower-limit log N_tot ~ 15.5 used in Table 2. If the corrected AV returns to about 0 (or to within its uncertainty), the quoted AV=1.2 is a self-absorption artifact and the optically thin assumption fails for the 'safe' progressions; if the corrected AV remains above 0.5 mag, then either the column is lower than assumed or there is genuine extra UV extinction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The method's core assumption is that the used H2 fluorescent lines are optically thin, so the observed flux ratios depend only on Einstein A coefficients, wavelengths, and dust extinction (Eq. 1, Sect. 3.1). This is load-bearing because self-absorption mimics dust, biasing AV high and RV low. The paper's own tau estimates (Table 2) are computed in Sect. 3.2 using a column density log N_tot ~ 15.5 that the authors explicitly describe as a lower limit. Adopting McJunkin et al. (2016) values (log N_tot ~ 19; TW Hya at 2500 K) scales tau for the [1,7] and [1,4] progressions (listed as < 10^-3 in Table 2) by roughly 10^3.5, making them of order unity or larger. Thus the claim that progressions [1,7] and [1,4] are usable for all stars is not secured. The single decisive control object, TW Hya, has a literature AV of 0.0 mag but yields AV=1.2±0.3 at fixed RV=3.1 (Table 3) and AV=0.3±0.1, RV=1.65 in the box method; in Sect. 5 the authors state that 'our procedure failed' for TW Hya. This is direct evidence that the optically thin assumption is violated even in the progressions used for the other 33 stars, so the central claim of reliable AV values and 'no additional UV extinction' is not independently established. Additionally, the post hoc exclusion of the [0,1] progression when it disagrees (Sect. 4) weakens the comparison with literature values in Fig. 4.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a method to measure interstellar extinction, A_V and R_V, for classical T Tauri stars (CTTSs) using far-UV H2 fluorescent line fluxes. The method assumes the selected H2 lines are optically thin, so that observed flux ratios within a progression depend only on known Einstein A coefficients, wavelengths, and the extinction curve (Eq. 1 in Sect. 3.1). The authors apply the method to 34 stars from the ULLYSES sample, combine results from up to three H2 progressions, and compare the resulting A_V values with literature values. They report a correlation with literature A_V (r = 0.70, or 0.80 without TW Hya) at fixed R_V = 3.1, and they discuss the degeneracy between A_V and R_V and the possible impact of self-absorption. The main conclusion is that H2-based A_V values are good estimates of the extinction toward the inner warm disk and that no significant additional UV extinction exists between the star and the disk.","tokens_in":19957,"tokens_out":4438,"duration_ms":47172,"significance":"If the central assumption is valid, the method provides an independent extinction probe of the inner disk environment, which is valuable for SED construction and accretion-rate estimates in CTTSs. The paper's strengths include the use of external atomic data, test simulations that recover input A_V within errors, and a transparent discussion of limitations. The method is falsifiable and the comparison to literature is a useful external check. However, the load-bearing assumption of optically thin lines is not secured by the paper's own estimates: the tau values in Table 2 are lower limits, and the single decisive control object (TW Hya) fails in the sense that the method recovers a significant nonzero A_V for a star with well-established near-zero extinction. These issues must be resolved before the central claims about reliable A_V and the absence of additional UV extinction can be accepted.","major_comments":[{"comment":"The tau estimates used to justify treating progressions [1,7] and [1,4] as optically thin are lower limits. In Sect. 3.2 the authors derive log N_tot ~ 15.5 and explicitly state that this is a lower limit because only some pumping lines were considered and a Voigt profile was set to unity. Using the average column density from McJunkin et al. (2016), log N_tot ~ 19.0, scales the tau values for the [1,7] and [1,4] lines (listed as < 10^-3 in Table 2) by roughly 10^3.5, which would make several of them order unity or larger. Consequently, the statement in Sect. 6 that progressions [1,4] and [1,7] 'seem to be usable for all stars' is not secured by the present analysis. Since self-absorption removes blue photons and mimics dust, the paper should either propagate the full plausible range of column densities into the tau estimates and show which individual lines remain optically thin, or benchmark the method against radiative-transfer calculations such as those in McJunkin et al. (2016).","section":"Sect. 3.2 and Table 2"},{"comment":"The TW Hya control is not merely a small discrepancy; it is a failure of the method on the only star with a securely known near-zero extinction. TW Hya yields A_V = 0.3 +/- 0.1 with R_V = 1.65 in the box search and A_V = 1.2 +/- 0.3 at fixed R_V = 3.1 (Table 3), while the literature value is A_V = 0.0. In Sect. 5 the authors state that 'our procedure failed' for TW Hya and that self-absorption may affect all progressions for this star. Because self-absorption mimics dust, this failure is exactly the expected signature of the mechanism the paper seeks to exclude, and it prevents the conclusions in Sect. 6 that the H2-derived A_V values are reliable for the other stars and that no significant additional extinction exists between the star and the disk. The paper should treat TW Hya as a quantitative bound on the systematic error of the method and should investigate whether the same bias, even if smaller, afflicts the other 33 stars.","section":"Sect. 4, Sect. 5, and Table 3"},{"comment":"The decision to exclude progression [0,1] from the A_V determination is made post hoc: in Sect. 4 the authors state that the progression 'led to deviating results and was omitted because the deviation of this progression might be self-absorption.' This creates a selection effect in the comparison with literature values shown in Fig. 4: whenever the [0,1] progression disagrees with the other two progressions, its data are removed, which removes exactly the data that would test the optically thin assumption. The correlation with literature (r = 0.70, or 0.80 without TW Hya) is therefore computed on a selected subset. The paper should define an a priori criterion for excluding a progression (for example, based on a tau threshold computed with the full column-density range) and should report the correlation both with and without the [0,1] progression, as well as the number of stars for which the exclusion changes the result.","section":"Sect. 4, Sect. 5, and Fig. 4"}],"minor_comments":[{"comment":"The sample count is presented inconsistently: the text says that the cut by H2 luminosity 'left 51 stars,' but then states that 'all 52 preselected stars (including TW Hya)' are listed in Table 1. Please clarify whether the initial cut left 51 or 52 stars before the S/N cut.","section":"Sect. 2.1"},{"comment":"The sentence 'For TW Hya and fixing R_V=3.1, we found a significant A_V=1.2 +/- 0.3 mag, without any indications of contamination by self-absorption' is difficult to reconcile with the later statement in Sect. 5 that 'our procedure failed' for TW Hya. Please rephrase to make clear that the lack of indication is not evidence against self-absorption, especially given the lower-limit nature of the tau estimates.","section":"Sect. 6"},{"comment":"The text notes that for stars where all three progressions were used, the H2-based A_V values are equal to or higher than the literature values, and that face-on disks tend to show higher H2-based extinction. This trend is in tension with the conclusion that no additional extinction exists between the star and the disk; the discussion would benefit from a quantitative assessment of how large a systematic offset could be hidden by the current error bars.","section":"Sect. 5 and Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about its limitations, and the literature correlation suggests that the method is not purely noise. However, the TW Hya failure and the lower-limit nature of the tau estimates strike at the central claim. In my view the manuscript needs a substantive revision that either (a) restricts the claims to lines whose optical depth is securely small using literature column densities, or (b) includes a quantitative treatment of the systematic bias from self-absorption, rather than a cosmetic rewrite. The fit to A&A is otherwise appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read: the paper is a solid, transparent application of the classic Miller line-ratio idea to far-UV H2 fluorescent lines in 34 CTTSs from ULLYSES, and it produces a new AV/RV catalog that broadly tracks optical literature values. The genuinely new piece is the sample size and the systematic multi-progression comparison; the best parts are the simulation test and the honest handling of bad lines. This deserves a serious referee.\n\nWhat is new and good: nobody has done a no-radiative-transfer H2 line-ratio extinction measurement on a sample this size. The AV values at fixed RV=3.1 are directly useful for SED and accretion-rate work in the ULLYSES sample, and the paper is careful about blends, weak lines, and the AV-RV degeneracy. The box search combining progressions is a reasonable way to break that degeneracy, and the simulation recovers input values.\n\nThe soft spot is exactly what the stress-test flags: the optically thin assumption is load-bearing, and the paper's own tau estimates are lower limits. Table 2 lists tau < 10^-3 for [1,7] and [1,4] using log N_tot ~15.5, which the authors explicitly call a lower limit. If the true column is closer to McJunkin's log N_tot ~19, those tau values could be order unity. TW Hya makes this concrete: a zero-AV star gives AV=1.2 at fixed RV=3.1, and the box method only recovers zero by forcing RV=1.65. The paper says 'our procedure failed' for TW Hya, yet then uses [1,7] and [1,4] for the rest and concludes there is no significant extra UV extinction. That is a real logical gap. It does not kill the empirical catalog, but it means the physical conclusion is not independently secured.\n\nMinor complaints: no machine-readable line flux table and no code, so verifying per-object results is annoying; and the post hoc dropping of [0,1] when it disagrees weakens the Fig. 4 comparison, even if it is defensible.\n\nBottom line: this is a useful paper for people working on CTTS SEDs and accretion diagnostics, and it is also a good teaching example of how a nice method can strain against its own assumption. It deserves peer review. I would encourage acceptance with major comments: state clearly that the tau lower limits do not rule out optical depth, and soften the conclusion about additional UV extinction.","headline":"New H2 line-ratio extinction catalog for 34 CTTSs with an honest but load-bearing optically thin assumption; worth a serious referee, though the TW Hya failure should force softer conclusions.","tokens_in":20549,"tokens_out":4541,"would_cite":true,"duration_ms":45245,"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":"Far-UV fluorescent H2 line ratios measure interstellar extinction toward classical T Tauri stars, largely matching optical values.","keywords":["interstellar extinction","classical T Tauri stars","H2 fluorescent lines","far-ultraviolet spectroscopy","reddening law","pre-main-sequence stars","protoplanetary disks","Hubble Space Telescope spectra"],"falsifier":"Measure the H2 line ratios of a second star with independently known zero extinction, at higher signal-to-noise than TW Hya, and apply the paper's method; if the inferred $A_V$ at the canonical $R_V=3.1$ again comes out near 1 mag while optical data say zero, the optically thin assumption fails. A more direct check is to compute the expected flux deficit from the paper's own optical-depth estimates, which reach $\\tau\\sim1.6$ for the bluest $[0,2]$ lines, and show that correcting for it removes the discrepancy.","tokens_in":19395,"feed_emoji":"🔭","tokens_out":10531,"duration_ms":101630,"temperature":0.7,"pith_summary":"This paper argues that the far-ultraviolet fluorescent H2 lines emitted by the warm inner disks of classical T Tauri stars can be used as an independent measure of interstellar and circumstellar extinction. Because lines sharing the same upper level must appear in fixed ratios set only by transition probabilities and wavelengths, any observed deviation from those ratios can be attributed to dust reddening. The authors apply this idea to Hubble Space Telescope spectra of 72 classical T Tauri stars and obtain $A_V$ and $R_V$ values for the 34 objects with sufficient line quality, finding broad agreement with extinction values derived from optical data. This matters because modeled ultraviolet spectra and accretion rates depend on how much ultraviolet light is assumed to be absorbed, and the agreement suggests no hidden extra ultraviolet extinction lies between the star and the inner disk.","feed_headline":"Far-UV H2 lines reveal dust toward 34 young stars","feed_subtitle":"Fluorescent hydrogen lines from the inner disk give extinction values that match optical measurements—and show no extra UV dust.","key_machinery":"The central object is a set of H2 Lyman-band fluorescence progressions: groups of emission lines that share the same upper rovibrational level in the $B^1\\Sigma_u^+$ state, so their intrinsic flux ratios are fixed by branching ratios. The load-bearing identity is $r_{\\mathrm{theo}} = F_{ik}/F_{im} = (\\lambda_{im} A_{ik})/(\\lambda_{ik} A_{im})$, which converts observed line ratios into an extinction measurement. The authors deredden measured fluxes on a grid of $A_V$ and $R_V$ using a standard parameterized extinction law, compute a chi-square-like statistic over all flux-ratio permutations, and then use a box search over the overlapping best-fit regions from the three usable progressions to locate the best $A_V$–$R_V$ pair. The $[0,2]$ progression is discarded because of blended lines and strong self-absorption, and $[0,1]$ is dropped for individual stars when its solution deviates systematically from the others.","core_discovery":"The paper's central claim is that optically thin far-UV H2 line ratios provide reliable extinction estimates for classical T Tauri stars, needing no radiative-transfer modeling. For each fluorescence progression, the theoretical ratio of any two lines is fixed by their Einstein $A$ coefficients and wavelengths; the observed ratios are dereddened across a grid of $A_V$ and $R_V$, and the pair that brings all ratios into best agreement is the extinction solution. Combining the overlap of the best-fit regions from the $[1,7]$, $[1,4]$, and $[0,1]$ progressions partly breaks the known $A_V$–$R_V$ degeneracy. The resulting values largely agree with optical determinations, and low-extinction calibration stars recover near-zero $A_V$; the notable exception is TW Hya, whose H2-inferred value at fixed $R_V=3.1$ is $A_V=1.2\\pm0.3$ mag despite a well-established zero extinction, which the authors attribute to self-absorption that mimics dust. On this basis the authors conclude that standard $A_V$ values are good estimates of extinction toward the inner warm disk, without a significant extra ultraviolet extinction component.","pith_inferences":["The same ratio technique could be extended to FUV-pumped lines redward of 1600 Å, or to other molecules such as CO, where the $A_V$–$R_V$ degeneracy is weaker, making it possible to pin down both parameters for individual disks.","The paper's tentative finding that face-on disks tend to show higher H2-inferred than optical extinction, though based on only three stars, predicts that a larger sample of face-on systems will show a systematic offset that distinguishes disk-rim extinction from self-absorption.","Comparing H2-based $A_V$ with extinction derived from Ly-alpha reconstructions would test where the dust sits: the two probes sample different disk heights, so agreement would place the dust in the intervening interstellar medium, while disagreement would reveal a circumstellar component."],"forward_implications":["For the 34 stars with sufficient data, far-UV H2 line ratios yield $A_V$ and $R_V$ values that largely agree with optical extinction measurements.","Stars with literature $A_V=0$ mag recover near-zero extinction, confirming that the redder H2 progressions are not seriously contaminated by self-absorption for those objects.","Assuming the canonical $R_V=3.1$ gives a mean $A_V$ for every star, providing a consistent set of extinction values for reddening modeled ultraviolet emission from accretion columns.","The general agreement between H2-based and optical $A_V$ implies that no significant additional ultraviolet extinction exists between the star and the inner disk.","Combining multiple H2 progressions breaks much of the $A_V$–$R_V$ degeneracy, because random noise affects the progressions differently and their best-fit regions overlap near the true solution."],"supporting_citations":[{"why":"Supplies the Einstein A coefficients and wavelengths for the H2 lines that enter the theoretical flux ratios.","marker":"Abgrall et al. (1993)"},{"why":"Provides the extinction curve parameterization used to deredden the measured line fluxes on the AV and RV grid.","marker":"Fitzpatrick (1999)"},{"why":"Earlier H2-based extinction study with radiative transfer; source of the self-absorption estimates and comparison values for six stars.","marker":"McJunkin et al. (2016)"},{"why":"Establishes the use of all flux-ratio permutations and the statistical properties of the chi-square-like statistic.","marker":"Hartigan & Morse (2007)"},{"why":"Introduced the shared-upper-level line-ratio technique that the method applies to H2 progressions.","marker":"Miller (1968)"},{"why":"Provides the far-UV H2 line library and defines the sample of classical T Tauri stars used as the starting point.","marker":"France et al. (2023)"},{"why":"Identifies the Ly-alpha-pumped H2 progressions in classical T Tauri stars and the optical depth behavior of their lines.","marker":"Herczeg et al. (2006)"},{"why":"Supplies the unred function used to apply the extinction law in the dereddening step.","marker":"Czesla et al. (2019)"}],"fun_headline_variants":["Far-UV H2 lines size up dust toward young stars","No extra UV dust: H2 lines match optical extinction","34 T Tauri stars get extinction from H2 line ratios","H2 fluorescence measures extinction in T Tauri stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole measurement rests on the assumption that the H2 fluorescent lines are optically thin, so dust is the only thing that changes their flux ratios; if the gas absorbs some of its own bluer fluorescent photons, that self-absorption mimics dust, inflating the inferred $A_V$ and lowering the inferred $R_V$.","fun_headline_variants_meta":{"raw":{"variants":["Far-UV H2 lines size up dust toward young stars","No extra UV dust: H2 lines match optical extinction","34 T Tauri stars get extinction from H2 line ratios","H2 fluorescence measures extinction in T Tauri stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00087,"raw_usage":{"total_tokens":3791,"prompt_tokens":991,"completion_tokens":2800,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":607,"completion_tokens_details":{"reasoning_tokens":2732}},"tokens_in":607,"tokens_out":2800,"duration_ms":21900,"temperature":1.0,"reasoning_tokens":2732,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:39:54.902919+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the H2 line ratios of a second star with independently known zero extinction, at higher signal-to-noise than TW Hya, and apply the paper's method; if the inferred $A_V$ at the canonical $R_V=3.1$ again comes out near 1 mag while optical data say zero, the optically thin assumption fails. A more direct check is to compute the expected flux deficit from the paper's own optical-depth estimates, which reach $\\tau\\sim1.6$ for the bluest $[0,2]$ lines, and show that correcting for it removes the discrepancy.","supporting_citations":[{"cited_title":"Y ., & Subtil, J","cited_arxiv_id":null,"evidence_quote":"Supplies the Einstein A coefficients and wavelengths for the H2 lines that enter the theoretical flux ratios."},{"cited_title":"2016, ApJ, 828, 69","cited_arxiv_id":null,"evidence_quote":"Earlier H2-based extinction study with radiative transfer; source of the self-absorption estimates and comparison values for six stars."},{"cited_title":"& Morse, J","cited_arxiv_id":null,"evidence_quote":"Establishes the use of all flux-ratio permutations and the statistical properties of the chi-square-like statistic."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduced the shared-upper-level line-ratio technique that the method applies to H2 progressions."},{"cited_title":"2023, AJ, 166, 67","cited_arxiv_id":null,"evidence_quote":"Provides the far-UV H2 line library and defines the sample of classical T Tauri stars used as the starting point."},{"cited_title":"J., Linsky, J","cited_arxiv_id":null,"evidence_quote":"Identifies the Ly-alpha-pumped H2 progressions in classical T Tauri stars and the optical depth behavior of their lines."}],"review_version":1}