REVIEW 3 major objections 3 minor 41 references
Measurement of interstellar extinction for classical T Tauri stars using far-UV H2 line fluxes
T0 review · 3 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Far-UV fluorescent H2 line ratios measure interstellar extinction toward classical T Tauri stars, largely matching optical values.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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$.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Sect. 3.2 and Table 2] 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).
- [Sect. 4, Sect. 5, and Table 3] 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.
- [Sect. 4, Sect. 5, and Fig. 4] 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.
minor comments (3)
- [Sect. 2.1] 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.
- [Sect. 6] 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.
- [Sect. 5 and Fig. 4] 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.
Circularity Check
No significant circularity: the H2-based AV/RV inversion is driven by external Einstein A values and a fixed extinction law, and the main caveats are physical validity of the optically thin assumption rather than circular reductions.
full rationale
The central derivation is self-contained. Eq. 1 computes theoretical flux ratios from externally tabulated Einstein A coefficients and wavelengths (Abgrall et al. 1993), and the observed ratios are dereddened with the Fitzpatrick (1999) extinction law; AV and RV are then obtained by minimizing a chi-square-like statistic on a grid (Eq. 2, Sect. 3.1). No parameter is fitted to literature AV values and then renamed as a prediction; the comparison to literature is a genuine external benchmark. The paper does rely on some co-authored prior work (France et al. 2023 for line lists, McJunkin et al. 2016 for tau estimates), but these are used as data sources and testable published results, not as unverified uniqueness claims that force the conclusion. The post hoc omission of the [0,1] progression when it deviates (Sect. 4) is a model-selection step that weakens the validation, but it does not make the reported AV values equal to the inputs by construction. The paper itself flags the strongest caveat: for TW Hya, a star with literature AV=0, the method gives AV=1.2 at fixed RV=3.1 and the authors state in Sect. 5 that 'our procedure failed'. Similarly, the tau values in Table 2 are based on a total column density that the authors explicitly call a lower limit, so the optically thin assumption is not fully secured. These are correctness and robustness risks, not circularity: the derivation can disagree with external data, and it visibly does for the TW Hya control. The central claim therefore retains independent empirical content, and no equation reduces to its own input.
Assumptions & free parameters
free parameters (4)
- AV (best-fit total extinction per star) =
0.0 to 3.0 mag depending on star (Table 3)
- RV (best-fit extinction curve shape per star) =
1.5 to 5.0 depending on star (Table 3)
- Search box half-widths for error estimation =
0.25 mag in AV, 0.3 mag in RV
- Number of best-fitting C models per progression =
160
assumptions (5)
- domain assumption All lines within a progression share the same upper-level population, so intrinsic flux ratios are set by Einstein A coefficients and wavelengths (Eq. 1).
- domain assumption The chosen H2 lines are optically thin.
- domain assumption The Fitzpatrick (1999) extinction law, via PyAstronomy's unred function, describes the wavelength dependence of extinction in these sight lines.
- domain assumption Measured Gaussian line flux represents the full H2 line flux from the same emitting region.
- standard math Atomic data (wavelengths and Einstein A values) from Abgrall et al. (1993) are accurate.
Cite this review
Pith. "Pith review of Measurement of interstellar extinction for classical T Tauri stars using far-UV H2 line fluxes." pith.science (2026). https://pith.science/paper/HKZX43HV
@misc{pith2026241113247,
author = {Pith},
title = {Pith review of: Measurement of interstellar extinction for classical T Tauri stars using far-UV H2 line fluxes},
year = {2026},
howpublished = {\url{https://pith.science/paper/HKZX43HV}},
note = {Machine review of arXiv:2411.13247}
}
abstract
Understanding the interstellar and potentially circumstellar extinction in the sight lines of classical T Tauri stars is an important ingredient for constructing reliable spectral energy distributions, which catalyze protoplanetary disk chemistry, for example. Therefore, some attempts of measuring $A_{V}$ toward individual stars have been made using partly different wavelength regimes and different underlying assumptions. We used strong lines of Ly{\alpha} fluorescent H2 and derived the extinction based on the assumption of optically thin transitions. We investigated a sample of 72 classical T Tauri stars observed with the Hubble Space Telescope in the framework of the ULLYSES program. We computed $A_{V}$ and $R_{V}$ values for the 34 objects with sufficient data quality and an additionally $A_{V}$ value for the canonical $R_{V}$ = 3.1 value. Our results agree largely with values obtained from optical data. Moreover, we confirm the degeneracy between $A_{V}$ and $R_{V}$ and present possibilities to break this. Finally, we discuss whether the assumption of optical thin lines is valid.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
Abgrall, H., Roueff, E., Launay, F., Roncin, J. Y ., & Subtil, J. L. 1993, A&AS, 101, 273 Alcalá, J. M., Manara, C. F., Natta, A., et al. 2017, A&A, 600, A20
work page 1993
-
[2]
Ardila, D. R., Basri, G., Walter, F. M., Valenti, J. A., & Johns-Krull, C. M. 2002, ApJ, 566, 1100
work page 2002
- [3]
-
[4]
A., Clayton, G
Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245
1989
-
[5]
Carvalho, A. S. & Hillenbrand, L. A. 2022, ApJ, 940, 156
work page 2022
-
[6]
P., et al
Czesla, S., Schröter, S., Schneider, C. P., et al. 2019, PyA: Python astronomy- related packages
2019
-
[7]
Draine, B. T. 2003, ARA&A, 41, 241
2003
-
[8]
Espaillat, C. C., Herczeg, G. J., Thanathibodee, T., et al. 2022, AJ, 163, 114
work page 2022
Show all 41 references
-
[9]
Fitzpatrick, E. L. 1999, PASP, 111, 63
1999
-
[10]
2023, AJ, 166, 67
France, K., Arulanantham, N., Maloney, E., et al. 2023, AJ, 166, 67
2023
-
[11]
F., et al
Gangi, M., Nisini, B., Manara, C. F., et al. 2023, A&A, 675, A153
2023
-
[12]
Grankin, K. N. 2017, in Astronomical Society of the Pacific Conference Series, V ol. 511, Non-Stable Universe: Energetic Resources, Activity Phenomena, and Evolutionary Processes, ed. A. M. Mickaelian, H. A. Harutyunian, & E. H. Nikoghosyan, 37
2017
-
[13]
C., Froning, C
Green, J. C., Froning, C. S., Osterman, S., et al. 2012, ApJ, 744, 60
2012
-
[14]
& Morse, J
Hartigan, P. & Morse, J. 2007, ApJ, 660, 426
2007
-
[15]
Herczeg, G. J. & Hillenbrand, L. A. 2014, ApJ, 786, 97
2014
-
[16]
J., Linsky, J
Herczeg, G. J., Linsky, J. L., Valenti, J. A., Johns-Krull, C. M., & Wood, B. E. 2002, ApJ, 572, 310
2002
-
[17]
J., Linsky, J
Herczeg, G. J., Linsky, J. L., Walter, F. M., Gahm, G. F., & Johns-Krull, C. M. 2006, ApJS, 165, 256
2006
-
[18]
J., Wood, B
Herczeg, G. J., Wood, B. E., Linsky, J. L., Valenti, J. A., & Johns-Krull, C. M. 2004, ApJ, 607, 369
2004
-
[19]
D., McJunkin, M., & Schneider, P
Hoadley, K., France, K., Alexander, R. D., McJunkin, M., & Schneider, P. C. 2015, ApJ, 812, 41
2015
-
[20]
2013, ApJ, 767, 112
Ingleby, L., Calvet, N., Herczeg, G., et al. 2013, ApJ, 767, 112
2013
-
[21]
E., Bartoe, J
Jordan, C., Brueckner, G. E., Bartoe, J. D. F., Sandlin, G. D., & van Hoosier, M. E. 1977, Nature, 270, 326
1977
-
[22]
Kenyon, S. J. & Hartmann, L. 1995, ApJS, 101, 117
1995
-
[23]
L., Mamajek, E
Luhman, K. L., Mamajek, E. E., Shukla, S. J., & Loutrel, N. P. 2017, AJ, 153, 46
2017
-
[24]
F., Testi, L., Herczeg, G
Manara, C. F., Testi, L., Herczeg, G. J., et al. 2017, A&A, 604, A127
2017
-
[25]
Mathis, J. S. 1990, ARA&A, 28, 37
1990
-
[26]
2016, ApJ, 828, 69
McJunkin, M., France, K., Schindhelm, R., et al. 2016, ApJ, 828, 69
2016
-
[27]
C., et al
McJunkin, M., France, K., Schneider, P. C., et al. 2014, ApJ, 780, 150
2014
-
[28]
Miller, J. S. 1968, ApJ, 154, L57
1968
-
[29]
& Schmitt, J
Predehl, P. & Schmitt, J. H. M. M. 1995, A&A, 293, 889
1995
-
[30]
Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618
1985
-
[31]
R., Taylor, J
Roman-Duval, J., Proffitt, C. R., Taylor, J. M., et al. 2020, Research Notes of the American Astronomical Society, 4, 205
2020
-
[32]
Rucinski, S. M. & Krautter, J. 1983, A&A, 121, 217
1983
-
[33]
2018, A&A, 609, A70
Rugel, M., Fedele, D., & Herczeg, G. 2018, A&A, 609, A70
2018
-
[34]
J., et al
Schindhelm, R., France, K., Herczeg, G. J., et al. 2012, ApJ, 756, L23
2012
-
[35]
C., France, K., Günther, H
Schneider, P. C., France, K., Günther, H. M., et al. 2015, A&A, 584, A51
2015
-
[36]
C., Manara, C
Schneider, P. C., Manara, C. F., Facchini, S., et al. 2018, A&A, 614, A108
2018
-
[37]
Schultz, G. V . & Wiemer, W. 1975, A&A, 43, 133
1975
-
[38]
Shull, J. M. 1978, ApJ, 224, 841
1978
-
[39]
1988, ApJ, 332, 400
Sternberg, A. 1988, ApJ, 332, 400
1988
-
[40]
Welty, D. E. & Fowler, J. R. 1992, ApJ, 393, 193
1992
-
[41]
E., Karovska, M., & Raymond, J
Wood, B. E., Karovska, M., & Raymond, J. C. 2002, ApJ, 575, 1057 Article number, page 10 of 13 B. Fuhrmeister et al.: Reddening determination by H2 measurements Appendix A: Total flux in a single progression (Method II) Here we discuss an alternative method that is conceptuall...
2002
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.