pith. sign in

arxiv: 2604.25590 · v1 · submitted 2026-04-28 · 🌌 astro-ph.GA · astro-ph.SR

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules. III. The absorption lines of [O I], CH, and OH

W.-J. Kim , A. M. Jacob , D. A. Neufeld , P. Schilke , H. Wiesemeyer , M. Gerin , M. G. Wolfire , V. Ossenkopf-Okada
show 12 more authors
V. Valdivia E. Falgarone D. C. Lis S. Bialy M. R. Rugel \'A. S\'anchez-Monge M. Busch T. M\"oller F. Wyrowski D. Seifried K. M. Menten A. Saintonge
This is my paper

Pith reviewed 2026-05-07 15:44 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.SR
keywords interstellar mediumabsorption lineshydridesmolecular tracersoxygen abundanceGalactic ISMdiffuse gasmolecular fraction
0
0 comments X

The pith

CH, OH, HCO+, and CCH absorption lines show strong mutual correlations with uniform column density ratios across the Galactic ISM, and the gas-phase oxygen abundance is (3.09 ± 0.64) × 10^{-4}.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This work measures absorption lines of neutral atomic oxygen, CH, and OH toward star-forming regions to probe the makeup of gas in the Milky Way. It reports tight correlations among the column densities of CH, OH, HCO+, and CCH, with the ratios between them staying the same in different locations and at different velocities. These ratios match earlier results from nearby diffuse clouds. The amount of oxygen in the gas relative to all hydrogen is found to be slightly below the solar value yet consistent with prior gas-phase measurements. Neutral hydrogen column density falls once the molecular fraction passes about 0.5, while the molecular species grow more abundant and HCO+ samples gas with higher molecular fractions than HI or the hydrides.

Core claim

The paper establishes that CH, OH, HCO+, and CCH exhibit strong mutual correlations, with column density ratios that remain uniform across Galactic environments and velocity intervals. It reports the average gas-phase oxygen abundance relative to total hydrogen as ⟨X(O)⟩ = N(O)/N(H_total) = (3.09 ± 0.64) × 10^{-4}. Atomic oxygen abundance stays roughly constant, while the abundances of OH, HCO+, and CCH increase with molecular fraction; HCO+ traces gas at higher molecular fractions than HI or the hydride ions.

What carries the argument

Mutual correlations among column densities of CH, OH, HCO+, and CCH from absorption lines, together with [O I] measurements to determine oxygen abundance and its variation with molecular fraction.

If this is right

  • The consistent ratios allow these species to serve as reliable H2 tracers throughout the diffuse Galactic ISM.
  • Neutral hydrogen declines and molecular abundances rise once the molecular fraction exceeds ~0.5, marking the atomic-to-molecular transition.
  • Atomic oxygen abundance remains steady while OH, HCO+, and CCH abundances grow with increasing molecular content.
  • HCO+ absorption samples gas at higher molecular fractions than HI or the hydrides, indicating density variations along lines of sight.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The uniformity suggests the underlying chemistry operates similarly across wide ranges of Galactic conditions.
  • The results could guide interpretation of hydride and molecular lines observed in external galaxies.
  • Multi-tracer observations at higher resolution could test whether small ratio deviations appear near energetic regions such as supernova remnants.

Load-bearing premise

Column densities from the absorption lines accurately reflect true abundances without major effects from optical-depth saturation, non-LTE excitation, or blended velocity components, and total hydrogen is correctly divided between atomic and molecular phases.

What would settle it

New spectra along additional sightlines that show CH/OH or OH/HCO+ column density ratios varying strongly with Galactic position, velocity, or environment, or that yield an oxygen abundance outside roughly 2.45–3.73 × 10^{-4}.

read the original abstract

The HyGAL Stratospheric Observatory for Infrared Astronomy (SOFIA) legacy program aims at characterizing the interstellar medium in the Milky Way using hydrides, [C II], and [O I] absorption lines with the 2.7 m SOFIA telescope toward twenty-five submillimeter-bright Galactic star-forming regions. As part of HyGAL, we investigated correlations among the known H$_2$ tracers -- CH and OH from SOFIA observations, and HCO$^+$ and CCH from ancillary absorption line data from ground-based telescopes. We also examined the abundance variation of neutral atomic oxygen, [O I], observed in absorption. CH, OH, HCO$^+$, and CCH all exhibit strong mutual correlations. OH in particular shows tight correlations with HCO$^+$ and CCH, reflecting their linked chemical and physical pathways. Column density ratios among these H$_2$ tracers are consistent with previous measurements in local diffuse clouds and remain uniform across Galactic environments and velocity intervals. The gas phase oxygen abundance relative to total hydrogen, $\langle X$(O)$\rangle=N$(O)/$N$(H$_{\rm total}$), is $(3.09\pm0.64)\times10^{-4}$, slightly below the elemental solar value but consistent with the previous observations measuring gas-phase abundances. We also find that $N$(HI) decreases toward the regions where the molecular fraction exceeds $f_{H_2}^N \sim 0.5$, marking the onset of the molecular phase. While the atomic oxygen abundance remains roughly constant, the abundances of OH, HCO$^+$, and CCH increase with the molecular fraction. Gas traced by the HCO$^+$ absorption corresponds to higher molecular fractions than that traced by HI and hydride ions, highlighting density variations in the diffuse-to-translucent ISM along different lines of sight.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper reports SOFIA observations of [O I], CH, and OH absorption lines toward 25 Galactic star-forming regions as part of the HyGAL program. It finds strong mutual correlations among column densities of CH, OH, HCO+, and CCH, with ratios that are uniform across Galactic environments and velocity intervals and consistent with prior local diffuse-cloud measurements. The gas-phase oxygen abundance is given as ⟨X(O)⟩ = N(O)/N(H_total) = (3.09 ± 0.64) × 10^{-4}, slightly below the solar value, and trends are noted with molecular fraction f_H2^N, including a decrease in N(H I) above f_H2^N ≈ 0.5.

Significance. If the column-density measurements hold, the work supplies a valuable set of empirical correlations and an abundance benchmark for the diffuse-to-translucent ISM over a wider range of lines of sight than previous studies. The reported uniformity of tracer ratios and the quantitative oxygen abundance with uncertainty provide concrete constraints for chemical models of the Galactic ISM.

major comments (3)
  1. [§3] §3 (Data reduction and column-density derivation): No curve-of-growth analysis, optical-depth checks, or multi-component profile modeling is described for the [O I], CH, or OH lines. Because the central claims—mutual correlations, uniform ratios, and ⟨X(O)⟩ = (3.09 ± 0.64) × 10^{-4}—rest directly on the accuracy of these N values, the absence of saturation or blending verification is load-bearing.
  2. [§4.3] §4.3 (Oxygen abundance): N(H_total) is constructed as N(H I) + 2N(H2), yet the partitioning and the ancillary N(H2) values used to compute ⟨X(O)⟩ are not cross-checked against the same absorption data; any systematic offset in the molecular fraction would propagate directly into the reported abundance and its uncertainty.
  3. [§4.1] §4.1 (Correlations): The statement that ratios remain uniform “across Galactic environments and velocity intervals” is presented without tabulated scatter, reduced-χ² values, or explicit treatment of upper limits and measurement errors on the individual column densities, making it impossible to judge whether the uniformity is statistically supported.
minor comments (2)
  1. [Figures] Figure captions and axis labels should explicitly state whether plotted quantities are total column densities or per-velocity-component values.
  2. [Abstract] The abstract states the sample comprises twenty-five sources but does not list the specific transitions observed for each species; a concise table of observed lines would improve clarity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thorough review and valuable comments, which have helped us identify areas where the manuscript can be strengthened. We address each major comment point by point below, indicating where revisions will be made to improve clarity and rigor without altering the core scientific conclusions.

read point-by-point responses

  1. Referee: §3 (Data reduction and column-density derivation): No curve-of-growth analysis, optical-depth checks, or multi-component profile modeling is described for the [O I], CH, or OH lines. Because the central claims—mutual correlations, uniform ratios, and ⟨X(O)⟩ = (3.09 ± 0.64) × 10^{-4}—rest directly on the accuracy of these N values, the absence of saturation or blending verification is load-bearing.

    Authors: We agree that explicit documentation of these procedures would strengthen the paper. While the column densities were derived using standard methods that included initial optical-depth assessments and checks for saturation (particularly for the stronger [O I] and OH lines), these steps were not described in sufficient detail in the submitted manuscript. In the revised version, we will expand §3 to include a dedicated subsection on the curve-of-growth analysis, optical-depth verification, and profile fitting. Where multi-component structure was evident in the spectra, we will note how it was handled; for lines showing potential saturation, we will clarify the criteria used to adopt the optically thin limit or apply corrections. revision: yes

  2. Referee: §4.3 (Oxygen abundance): N(H_total) is constructed as N(H I) + 2N(H2), yet the partitioning and the ancillary N(H2) values used to compute ⟨X(O)⟩ are not cross-checked against the same absorption data; any systematic offset in the molecular fraction would propagate directly into the reported abundance and its uncertainty.

    Authors: The N(H2) values are drawn from ancillary literature data (primarily UV absorption and other tracers, as referenced in the text), since the SOFIA observations target [O I], CH, and OH rather than direct H2 lines. We acknowledge that a direct cross-check using only the present absorption dataset is not feasible. However, we will revise §4.3 to include an explicit comparison of molecular fractions inferred from our CH and OH column densities against the adopted N(H2) values, and we will discuss potential systematic offsets. The quoted uncertainty on ⟨X(O)⟩ already incorporates scatter from the sample; we will make the propagation of molecular-fraction uncertainties more transparent and add a sensitivity test. revision: partial

  3. Referee: §4.1 (Correlations): The statement that ratios remain uniform “across Galactic environments and velocity intervals” is presented without tabulated scatter, reduced-χ² values, or explicit treatment of upper limits and measurement errors on the individual column densities, making it impossible to judge whether the uniformity is statistically supported.

    Authors: We concur that quantitative statistical support is needed. In the revised manuscript, we will add a table in §4.1 (or an appendix) that reports the mean column-density ratios, their standard deviations, reduced-χ² values for the linear fits, and the number of detections versus upper limits. Upper limits will be handled using appropriate censored-data statistics (e.g., Kaplan–Meier or survival-analysis methods). This will allow readers to assess the degree of uniformity directly. revision: yes

Circularity Check

0 steps flagged

No circularity: direct observational measurements and empirical ratios

full rationale

The paper reports column densities derived from observed absorption line equivalent widths and profile fits for [O I], CH, and OH, together with ancillary data for HCO+ and CCH. It then computes empirical correlations and the simple ratio ⟨X(O)⟩ = N(O)/N(H_total) directly from those measured columns. No model is fitted to a subset of the data and then used to 'predict' a related quantity; no ansatz or uniqueness theorem is imported via self-citation to force the result; and the uniformity statements are presented as observed facts rather than derived predictions. The derivation chain is therefore self-contained against external benchmarks and contains no load-bearing step that reduces to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard astrophysical assumptions for converting absorption-line equivalent widths to column densities and for estimating total hydrogen columns; no new free parameters, ad-hoc entities, or invented quantities are introduced.

axioms (2)
  • domain assumption Absorption-line equivalent widths can be converted to column densities using standard curve-of-growth or excitation assumptions without significant saturation or blending effects
    Required to report N(O), N(CH), N(OH) and their correlations from the observed spectra
  • domain assumption Total hydrogen column density can be reliably partitioned into atomic and molecular components using HI and H2 tracers
    Underlies the molecular-fraction trends and the statement that N(HI) decreases above f_H2^N ~ 0.5

pith-pipeline@v0.9.0 · 5779 in / 1576 out tokens · 45518 ms · 2026-05-07T15:44:28.137017+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

5 extracted references

  1. [1]

    Akritas, M. G. & Bershady, M. A. 1996, ApJ, 470, 706 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 Balashev, S. A., Gupta, N., & Kosenko, D. N. 2021, MNRAS, 504, 3797 Balser, D. S., Rood, R. T., Bania, T. M., & Anderson, L. D. 2011, ApJ, 738, 27 Baumgartner, W. H. & Mushotzky, R. F. 2006, ApJ, 639, 929 Bellomi, E., Godard, B.,...

    1996

  2. [2]

    Another column density table for the velocity interval for VI(M) is available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5)

    Appendix C: Column density spectra Figure C shows the column density profile per velocity interval for each detected species in each observed source, and TableC.1 lists the fraction of the full derived column densities over the specific velocity intervals for VI(A). Another column density table for the velocity interval for VI(M) is available at the CDS v...

    2014

  3. [3]

    ) are determined by adopting the abundance of CH relative to H2,N(CH)/N(H 2)=3.5×10 −8 (Sheffer et al. 2008). For the case (a) in the absence of CH data, we used the measured column density ratio of N(OH)/N(CH)=3.29 to convert toN(H 2). (†) The atomic hydrogen column densities are obtained from Rugel et al. (in prep.). Article number, page 16 of 20 W.-J. ...

    2008

  4. [4]

    C.1: Channel-wise column density spectra (Nd/d3) as a function of velocity

    Fig. C.1: Channel-wise column density spectra (Nd/d3) as a function of velocity. Same as Fig. 2 Article number, page 17 of 20 A&A proofs:manuscript no. aa58863-26 Appendix D: Chemical networks of the observed species C2H+ 2 CH2 C2H+ C+ 2 CH+ 2 CH+ 3 CCH CH C+ 3 C+ C+ H2 H2 H2 e− e− (A) C+ HCO+ C3H+ C2 e− H2 H2 HorOH C C photons orCR O H2 e− e− e− e− C+ C ...

    1986

  5. [5]

    The large dispersions in our data points may result from line-of-sight integration through multiple ISM phases with distinct metallicities, particularly between 4 and 8 kpc

    (solid red line). The large dispersions in our data points may result from line-of-sight integration through multiple ISM phases with distinct metallicities, particularly between 4 and 8 kpc. Because absorption-line measurements are integrated over velocity intervals defined by line profiles, they inevitably sample gas components with different physical c...

    2011