# CO+ as a probe of the origin of CO in diffuse interstellar clouds

M. Gerin, H. Liszt

► To cite this version:

M. Gerin, H. Liszt. CO+ as a probe of the origin of CO in diffuse interstellar clouds. Astronomy and Astrophysics - A&A, 2021, 648, pp.A38. 10.1051/0004-6361/202039915. hal-03194930

HAL Id: hal-03194930

https://hal.science/hal-03194930

Submitted on 9 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
# CO⁺ as a probe of the origin of CO in diffuse interstellar clouds*

M. Gerin<fnref n="1" /> and H. Liszt<fnref n="2" />

Received 15 November 2020 / Accepted 28 January 2021

## ABSTRACT

*Context.* The chemistry of the diffuse interstellar medium is driven by the combined influences of cosmic rays, ultraviolet (UV) radiation, and turbulence. Previously detected at the outer edges of photodissociation regions and formed from the reaction of C⁺ and OH, CO⁺ is the main chemical precursor of HCO⁺ and CO in a thermal, cosmic-ray, and UV-driven chemistry.

*Aims.* Our aim was to test whether the thermal cosmic-ray and UV-driven chemistry is producing CO in diffuse interstellar molecular gas through the intermediate formation of CO⁺.

*Methods.* We searched for CO⁺ absorption with the Atacama Large Millimeter Array toward two quasars with known Galactic foreground absorption from diffuse interstellar gas, J1717-3342 and J1744-3116, targeting the two strongest hyperfine components of the J = 2–1 transition near 236 GHz.

*Results.* We could not detect CO⁺ but obtained sensitive upper limits toward both targets. The derived upper limits on the CO⁺ column densities represent about 4% of the HCO⁺ column densities. The corresponding upper limit on the CO⁺ abundance relative to H₂ is <1.2 × 10⁻¹⁰.

*Conclusions.* The non-detection of CO⁺ confirms that HCO⁺ is mainly produced in the reaction between oxygen and carbon hydrides, CH₂ or CH₃, induced by suprathermal processes, while CO⁺ and HOC⁺ result from reactions of C⁺ with OH and H₂O. The densities required to form CO molecules at low extinction are consistent with this scheme.

**Key words.** ISM: clouds – ISM: molecules – radio lines: ISM

## 1. Introduction

The diffuse interstellar medium (ISM) hosts a rich chemistry where many species reach molecular abundances similar to those in denser and darker regions despite the lower densities (a few tens to a few hundred particles per cm⁻³) and the relatively unattenuated far-ultraviolet (FUV) illumination by the interstellar radiation field (Snow & McCall 2006; Gerin et al. 2016, 2019; Liszt et al. 2014a,b, 2018). Diffuse and translucent interstellar clouds (Snow & McCall 2006) have long been used as testbeds of interstellar chemistry (Glassgold & Langer 1975, 1976; Black & Dalgarno 1977; Van Dishoeck & Black 1986, 1988). These efforts were successful in explaining the measured column densities for H₂ and small trace molecules such as CH, OH, and CN using chemical reactions at temperatures 30–80 K corresponding to kinetic temperatures measured in H₂ (Savage et al. 1977), but in some cases assuming regions of density ≥1000 cm⁻³ along the line of sight. However, the observed column densities of CO and the common presence of high column densities of CH⁺ presented challenges to such models of thermal, cosmic-ray, and UV-driven ion-chemistry in diffuse molecular gas.

The observed CO column densities are empirically explained by the serendipitous detection of unexpectedly high column densities of HCO⁺ seen in absorption at 89.2 GHz (Lucas & Liszt 1996). A nearly constant abundance of HCO⁺ relative to H₂, X(HCO⁺) = N(HCO⁺)/N(H₂) = 3 × 10⁻⁹ ± 0.2 dex is found in

comparison with OH and CH (Gerin et al. 2019), using their fixed abundances relative to H₂ established through direct absorption observations at optical and UV wavelengths (Sheffer et al. 2008; Weselak et al. 2009, 2010). With such significant column densities, the observed HCO⁺ ions will recombine with ambient thermal electrons at densities n(H) ≥ 100 cm⁻³ to produce the observed CO column densities (Liszt & Lucas 2000; Visser et al. 2009).

The gas phase chemistry developed during the 1970s–1980s era was qualitatively correct in predicting a fixed abundance ratio N(HCO⁺)/N(OH) via a dominant chemical chain that proceeds from C⁺ + OH → CO⁺ + H (Glassgold & Langer 1975, 1976; Dagdigian 2019) to CO⁺ + H₂ → HCO⁺ + H and HCO⁺ + e⁻ → CO + H. However, the observed N(HCO⁺)/N(OH) ≈ 1/30 is roughly thirty times higher than the predicted value through this simple chemical scheme. Some HCO⁺ can also be produced in this chemistry by the reaction of C⁺ and H₂O, of which the abundance is 28% of that of OH (Gerin et al. 2016; Wiesemeyer et al. 2016). Both reactions, CO⁺ + H₂ and C⁺ + H₂O, produce equal amounts of HCO⁺ and the slightly less tightly bonded isomer HOC⁺ that is observed with N(HOC⁺)/N(HCO⁺) = 0.015 ± 0.003 (Gerin et al. 2019). This disparity in the abundances of HCO⁺ and HOC⁺ is much too large to be explained by the isomerization reaction H₂ + HOC⁺ → H₂ + HCO⁺ as discussed below.

To reproduce the observed HCO⁺ abundance at the low mean densities of diffuse clouds, the most successful models include transient production mechanisms that involve dynamical processes in magnetized gas such as low-velocity shocks or

<footnote n="1">¹ LERMA, Observatoire de Paris, PSL Research University, CNRS, Ecole Normale Supérieure, Sorbonne Université, 75005 Paris, France.</footnote>

<footnote>e-mail: maryvonne.gerin@lra.ens.fr</footnote>

<footnote n="2">² National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903, USA</footnote>

<footnote>* Spectra are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/648/A38.</footnote>
turbulent vortices, in which ions and neutral species are partially decoupled (Godard et al. 2014; Lesaffre et al. 2020). The resulting slow velocity drift can therefore be used as an additional energy source to drive endothermic chemical reactions that could not operate otherwise (Sheffer et al. 2008; Visser et al. 2009), especially the reaction C⁺ + H₂ → CH⁺ + H, which is endothermic by 0.36 eV or about 4300 K. These models therefore include a specific chemical pathway to CO, which starts from CH⁺ and is followed by a series of rapid hydrogen abstraction reactions producing CH₂⁺ and CH₃⁺. These two ions react with atomic oxygen to produce the required HCO⁺, and not HOC⁺. In this scheme, HCO⁺ is once again the precursor of CO, but the predicted CO and HCO⁺ abundances fit the observations well, and the relative rarity of HOC⁺ arises very naturally.

Therefore, the relative abundances of the three molecular ions CO⁺, HCO⁺, and HOC⁺ provide key information concerning the mechanisms at the origin of CO in diffuse gas. Despite the pivotal role played by CO⁺ in some versions of the chemistry, its abundance relative to H₂ is much less well constrained than those of HCO⁺ and HOC⁺. CO⁺ was first identified in the interstellar medium toward the bright photodissociation region (PDR) M17SW and the young, high-excitation planetary nebula NGC 7027 (Latter et al. 1993). CO⁺ is now routinely observed in dense PDRs such as the Orion Bar (Goicoechea et al. 2017) or MonR2 (Treviño-Morales et al. 2016). In such regions, CO⁺ is located at the very edge close to the HI/H₂ transition, and its abundance relative to HCO⁺ ranges from one to several tens of percent (Fuente et al. 2003). In dense PDRs, the intense FUV radiation can pump H₂ into vibrationally excited levels that can reach a significant population (Agúndez et al. 2010). The internal energy of this vibrationally excited H₂ can trigger the endothermic reaction of C⁺ + H₂ forming CH⁺, as shown by the widespread distribution of CH⁺ emission in such regions (Goicoechea et al. 2019). Paradoxically, the presence of CO⁺ in dense PDRs is a sign that nonthermal chemical processes are working, while nonthermal processes are invoked in models of diffuse molecular gas to remove the need to produce CO⁺ on the path to CO.

There are no reported detections of CO⁺ in the diffuse and translucent interstellar medium: in this paper, we present observations with the Atacama Large Millimeter Array (ALMA) that constrain the abundance of CO⁺ there for the first time. The observation strategy is presented in Sect. 2, and the results are discussed in Sect. 3 together with the implications for the CO chemistry that are emphasized in Sects. 4 and 5. The conclusions are summarized in Sect. 6.

## 2. Spectroscopy and observations

### 2.1. Spectroscopy

Absorption spectroscopy represents the best method to probe the molecular content of diffuse clouds, because the molecular excitation is weak at the low densities (a few tens to a few hundred cm⁻³) and pressure (p/k ≈ 3 × 10³–10⁴ K-cm⁻³) of these regions (Jenkins & Tripp 2011; Gerin et al. 2015; Goldsmith et al. 2018), and the level populations are concentrated in the lowest rotational levels. As a ²Σ molecular ion, the energy levels of CO⁺ are described by three quantum numbers; N, the rigid body angular momentum quantum number; S = 1/2, the electron spin angular momentum quantum number; and J = N + S = N ± 1/2, the total angular momentum quantum number (Sastry et al. 1981).

The CO⁺ ground state rotational transitions at 117.7 GHz and 118.1 GHz are close to a strong atmospheric line from molecular

<table><caption>**Table 1.** Spectroscopic parameters of observed lines.</caption><thead><tr><th>Transition <sup>(a)</sup><br>N, J</th><th>Frequency<br>(MHz)</th><th>A<br>(s<sup>-1</sup>)</th><th>N(CO<sup>+</sup>)/∫<sub>τ</sub> dv <sup>(b)</sup><br>(cm<sup>-2</sup> km<sup>-1</sup>s)</th></tr></thead><tbody><tr><td>2,3/2–1 3/2</td><td>235380.046</td><td>0.689 × 10<sup>-4</sup></td><td>1.01 × 10<sup>14</sup></td></tr><tr><td>2,3/2–1 1/2</td><td>235789.641</td><td>3.465 × 10<sup>-4</sup></td><td>2.00 × 10<sup>13</sup></td></tr><tr><td>2,5/2–1 3/2</td><td>236062.553</td><td>4.172 × 10<sup>-4</sup></td><td>1.12 × 10<sup>13</sup></td></tr></tbody></table><table><thead><tr><th>Transition <sup>(a)</sup><br>J</th><th>Frequency<br>(MHz)</th><th>A<br>(s<sup>-1</sup>)</th><th>N(C<sup>13</sup>CO)/∫<sub>τ</sub> dv <sup>(c)</sup><br>(cm<sup>-2</sup> km<sup>-1</sup>s)</th></tr></thead><tbody><tr><td>2–1</td><td>220398.684</td><td>6.076 × 10<sup>-7</sup></td><td>2.42 × 10<sup>15</sup></td></tr></tbody></table><table><thead><tr><th>Transition <sup>(a)</sup><br>J</th><th>Frequency<br>(MHz)</th><th>A<br>(s<sup>-1</sup>)</th><th>N(C<sup>18</sup>O)/∫<sub>τ</sub> dv <sup>(c)</sup><br>(cm<sup>-2</sup> km<sup>-1</sup>s)</th></tr></thead><tbody><tr><td>2–1</td><td>219560.354</td><td>6.012 × 10<sup>-7</sup></td><td>2.42 × 10<sup>15</sup></td></tr></tbody></table>

Notes. <sup>(a)</sup>The spectroscopic data were extracted from the Cologne Database for Molecular Spectroscopy (CDMS; Müller et al. 2001, 2005; Endres et al. 2016). <sup>(b)</sup>For an excitation temperature of 2.73 K. <sup>(c)</sup>For an excitation temperature of 5 K.

oxygen at 118.75 GHz, rendering their observation from the ground very difficult. Hence, observations of CO⁺ have generally targeted excited transitions where the sky transmission is much better. We chose to search for the N = 2 → 1 transitions because the frequencies of this first excited transition near 236 GHz are easily accessible from ground-based radio observatories, and the level population remains significant even in diffuse interstellar gas. This N = 2 → 1 transition is split into three hyperfine components, with the strongest component at 236.06 GHz. Table 1 lists the three targeted transitions, their frequencies, Einstein A-coefficient, and the conversion factor between the integrated line opacity and the molecule column density N(CO⁺)/∫<sub>τ</sub> dv assuming that the energy-level population is determined by the cosmic microwave background at 2.73 K.

Apart from CO⁺, the setting of the ALMA correlator included spectral windows targeting the ¹³CO(2–1), C¹⁸O(2–1) and H₂CO(3₀,₃–2₀,₂) transitions that were accessible with the same receiver tuning while including a dedicated continuum spectral window for phase calibration. Properties of the CO lines are given in Table 1, where the column density/optical depth conversion is calculated for an excitation temperature of 5 K that is typical for this more easily excited species in diffuse molecular gas (Goldsmith et al. 2018).

### 2.2. Target directions

The targeted sources are J1717-3342 and J1744-3116, two bright quasars situated at small galactic latitude behind the Galactic bulge that were known to have high intervening column densities of Galactic neutral gas detected in molecular absorption (Gerin & Liszt 2017; Liszt & Gerin 2018; Riquelme et al. 2018). Their positions, continuum flux densities S<sub>v</sub>, channel-channel baseline rms line/continuum noise and line profile integrals, and velocity intervals used here as integration intervals for the targeted CO⁺ transitions are given in Table 2. Comparable quantities for the ¹³CO(2-1) and C¹⁸O(2-1) lines are given in Table 3.

### 2.3. ALMA observations

The observations were performed with ALMA Band-6 receivers tuned near 236 GHz during ALMA Cycle 7 under the project
<table><caption>**Table 2.** Summary of CO⁺ observations.</caption><thead><tr><th>Source</th><th><i>l</i><br>(°)</th><th><i>b</i><br>(°)</th><th>LSR velocity <sup>(a)</sup><br>(km s<sup>-1</sup>)</th><th>Line</th><th>Flux density <i>S</i><sub>v</sub><br>(Jy)</th><th>σ<sub>l/c</sub></th><th>∫ <i>τ</i> dv <sup>(b)</sup><br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td rowspan="2">J1717-3342</td><td rowspan="2">352.7333</td><td rowspan="2">2.3911</td><td rowspan="2">-14..10</td><td>CO<sup>+</sup></td><td><i>J</i>=5/2-3/2</td><td>0.780</td><td>0.0027</td><td>&lt;0.027</td></tr><tr><td>CO<sup>+</sup></td><td><i>J</i>=3/2-1/2</td><td>0.780</td><td>0.0034</td><td>&lt;0.033</td></tr><tr><td rowspan="2">J1744-3116</td><td rowspan="2">357.8635</td><td rowspan="2">-0.9968</td><td rowspan="2">-20..20</td><td>CO<sup>+</sup></td><td><i>J</i>=5/2-3/2</td><td>0.234</td><td>0.0084</td><td>&lt;0.087</td></tr><tr><td>CO<sup>+</sup></td><td><i>J</i>=3/2-1/2</td><td>0.251</td><td>0.0107</td><td>&lt;0.111</td></tr></tbody></table>

Notes. <sup>(a)</sup>Velocity interval used to integrate the line profile. <sup>(b)</sup>Upper limits are 3σ.

<figure><img src="image_1.png" /><figcaption>**Fig. 1.** Absorption spectra observed toward J1717-3342 (left) and J1744-3116 (right). The spectra were normalized by the continuum flux density. The CO⁺ spectra were shifted vertically and multiplied by 10. HCO⁺ is shown in black, CO⁺ in red, and <sup>13</sup>CO(2–1) in blue for comparison. The HCO⁺ data are taken from Liszt & Gerin (2018) and Gerin & Liszt (2017).</figcaption></figure>

<table><caption>**Table 3.** CO isotopologue *J* = 2–1 observations and column densities.</caption><thead><tr><th>Source</th><th>Line</th><th>Flux<br>(Jy)</th><th>σ<sub>l/c</sub></th><th>∫ <i>τ</i> dv<br>(km s<sup>-1</sup>)</th><th><i>N</i> <sup>(a)</sup><br>(10<sup>14</sup> cm<sup>-2</sup>)</th></tr></thead><tbody><tr><td rowspan="2">J1717-3342</td><td><sup>13</sup>CO</td><td>0.818</td><td>0.0024</td><td>0.260(0.004)</td><td>6.3(0.4)</td></tr><tr><td>C<sup>18</sup>O</td><td>0.822</td><td>0.0123</td><td>&lt;0.055</td><td>&lt;1.3</td></tr><tr><td rowspan="2">J1744-3116</td><td><sup>13</sup>CO</td><td>0.244</td><td>0.0073</td><td>2.975(0.023)</td><td>72.0(0.6)</td></tr><tr><td>C<sup>18</sup>O</td><td>0.244</td><td>0.0084</td><td>0.102(0.018)</td><td>2.5(0.4)</td></tr></tbody></table>

Notes. <sup>(a)</sup>Upper limit is 3σ.

code 2019.1.00120.S. The CO⁺ lines were observed with a channel spacing of 244 kHz, corresponding to a velocity resolution of 0.31 km s⁻¹. The channel spectral resolution is twice the channel spacing. As in our earlier observations toward these sources, the bandpass calibrator was fixed to J1924-2914. Absorption spectra were extracted from the standard pipeline-processed data products at the peak of the continuum map in each spectral window. The ALMA spectra of <sup>13</sup>CO, CO⁺, and HCO⁺ are displayed in Fig. 1. The H₂CO (3₀,₃ – 2₀,₂) line was detected in both directions and will be discussed in a future paper.

## 2.4. Other data discussed here

Also shown in Figs. 1 and 2 are 89.2 GHz ALMA HCO⁺ *J*=1–0 absorption line profiles with 0.205 km s⁻¹ channel spacing and 0.41 km s⁻¹ spectral resolution from our earlier works (Gerin & Liszt 2017; Liszt & Gerin 2018). Profile integrals are quoted in Table 4, assuming excitation in equilibrium with the CMB and using N(HCO⁺) = 1.12 × 10¹² cm⁻² ∫ τ dv as before. Also shown in Fig. 2 are the nearest λ2.6mm CO *J*=1–0 emission profiles

<figure><img src="image_2.png" /><figcaption>**Fig. 2.** Comparison of the HCO⁺ (1–0) (black) and <sup>13</sup>CO (2–1) (green) absorption with the <sup>12</sup>CO (blue) and scaled HI (gray) emission. J1717-3342 is shown at the top and J1744-3116 at the bottom. C<sup>18</sup>O(2–1) is also shown in orange for J1744-3116.</figcaption></figure>
<table><caption>**Table 4.** HCO⁺ and CO⁺ column densities.</caption><thead><tr><th>Source</th><th>N(HCO<sup>+</sup>)<br>(10<sup>12</sup> cm<sup>-2</sup>)</th><th>N(CO<sup>+</sup>)<br>(10<sup>12</sup> cm<sup>-2</sup>)</th><th>N(CO<sup>+</sup>)/N(HCO<sup>+</sup>)</th><th>N(CO<sup>+</sup>)/N(H<sub>2</sub>)<sup>(a)</sup></th></tr></thead><tbody><tr><td>J1717-3342</td><td>7.49 ± 0.03</td><td>&lt;0.30</td><td>&lt;0.040</td><td>&lt;1.2 × 10<sup>-10</sup></td></tr><tr><td>J1744-3416</td><td>23.35 ± 0.06</td><td>&lt;0.98</td><td>&lt;0.042</td><td>&lt;1.3 × 10<sup>-10</sup></td></tr></tbody></table>

Notes. (a) $N(H_2) = N(HCO^+)/3 \times 10^{-9}$.

with an 8.8 beamsize from the work of Bitran et al. (1997) and a λ21cm H I emission profile from Kalberla et al. (2005) taken with a FWHM beam of 0.6°. The CO emission profile for J1744-3116 at l, b = 357.875°, -1° is very nearly along the same sightline; for J1717-3342, the CO profile at 352.75°, 2.5° is centered 6.6 away.

# 3. Results and discussion

## 3.1. The sight lines toward J1717-3342 and J1744-3116

Absorption spectra of HCO⁺, ¹³CO, and CO⁺ toward both sources are shown in Fig. 1, and Fig. 2 presents a comparison of the absorption spectra of CO and HCO⁺ with emission spectra of HI and ¹²CO(1–0) described in Sect. 2.4. Column densities of HCO⁺, ¹³CO, C¹⁸O, and H₂ ($N(H_2) = N(HCO^+)/3 \times 10^{-9}$) are given in Tables 4 and 3. To compute $N(¹³CO)$ and $N(C¹⁸O)$, we assumed an excitation temperature of 5 K (Table 1), which is consistent with the observed brightness of the ¹²CO emission that is seen toward or near these sources, 1–2.5 K.

The diffuse nature of the gas along these lines of sight has previously been discussed by Gerin & Liszt (2017), Liszt & Gerin (2018), and Riquelme et al. (2018). High gas column densities accumulate over long paths through the Galactic disk at low galactic latitude, even without encountering dense clouds. This is indicated in our new data by the large ratio of ¹³CO and C¹⁸O integrated optical depths or column densities, $N(¹³CO)/N(C¹⁸O) = 29 \pm 5$ toward J1744-3116 (Table 3). If the oxygen and carbon ratios have the usual interstellar values, ¹⁶O/¹⁸O ~520 and ¹²C/¹³C ~65 (Wilson & Rood 1994; Milam et al. 2005; Keene et al. 1998; Giannetti et al. 2014), a fully molecular gas would have $N(¹³CO)/N(C¹⁸O) = 8$ toward J1744-3116. The much larger value we observe arises through a combination of selective photodissociation of C¹⁸O, and, more likely, fractionation of ¹³CO through endothermic carbon isotope exchange when the majority of the free gas phase carbon is in C⁺ (Liszt 2017). Our value for $N(¹³CO)$ toward J1744-3116 in Table 3 agrees with that given by Riquelme et al. (2018) based on $J=1-0$ absorption, $N(¹³CO) = 6.3 \times 10^{15}$ cm⁻².

In our earlier work we also derived carbon isotope ratios $H^{12}CO^+/H^{13}CO^+ = 58 \pm 9$ and $64 \pm 4$ in the low-velocity gas toward J1717-3342 and J1744-3116, respectively. These values are typical of Galactic disk gas near the Sun and quite different from what is observed inside the Galactic bulge. Therefore, although the sight lines employed here cross the Galactic bulge, the low-velocity gas that is analyzed in CO and CO⁺ resides relatively nearby in the disk and is discussed under the assumption that local conditions are applicable.

Overall, the CO(1–0) emission is spread over the same velocity range as the HCO⁺ absorption. Toward J1744-3116, there is a fair correspondence of the main CO emission peaks with the HCO⁺ absorption features, but the HCO⁺ absorption profiles are more complex. The positive velocity component toward

J1744-3116 is detected in ¹³CO(2–1) and C¹⁸O(2–1) absorption and is associated with a self-absorption feature in the H I spectrum that was discussed in Liszt & Gerin (2018). The presence of the three CO isotopologues indicates that this velocity component probably has the largest extinction and mean density and could be associated with translucent rather than diffuse gas. There is no self-absorption in the HI profile coincident with the negative velocity component, but the main ¹³CO(2–1) absorption is associated with a faint C¹⁸O(2–1) feature.

The CO(1–0) emission differs more from the HCO⁺ absorption toward J1717-3342, where the direction of the emission spectrum is 6.6' away. Of the eight main HCO⁺ absorption features, only one is detected in ¹³CO(2–1), and none in C¹⁸O(2–1), showing that the gas is more diffuse along this line of sight. The CO emission is weaker than toward J1744-3116, barely reaching 1 K. As discussed by Liszt et al. (2010), the CO(1–0) emission is highly variable in diffuse gas as it traces the regions where CO reaches column densities above $N(CO) = 10^{15}$ cm⁻² that are high enough to produce detectable emission. This transition occurs over a small range of H₂ column densities and depends on the local physical conditions, but CO emission at the level of 1 K–km s⁻¹ integrated intensity is not observed in regions where the molecular fraction of H-nuclei in H₂ is much below 0.6.

The comparison of the emission and absorption profiles therefore demonstrates that the two observed lines of sight encounter matter with variable physical conditions, ranging from diffuse to translucent gas. The upper limits for the CO⁺ integrated opacities shown in Table 2 provide constraints on the CO⁺ column densities and abundances for the whole range of physical conditions. For the diffuse molecular gas sampled in absorption the molecule excitation is dominated by the cosmic microwave background because of the low gas densities. Therefore, the CO⁺ column densities have been derived assuming an excitation temperature of 2.73 K. They are reported in Table 4 together with the HCO⁺ column densities in the same velocity intervals. The achieved upper limit on CO⁺ column densities leads to a very low value of about 0.04 for the $N(CO^+)/N(HCO^+)$ ratio. The corresponding upper limit on the CO⁺ abundance relative to H₂ is $X(CO^+) < 1.2 \times 10^{-10}$, using the standard value for the HCO⁺ abundance relative to H₂ in the diffuse molecular gas, $3 \times 10^{-9}$ (Lucas & Liszt 1996; Gerin et al. 2019). These values are of a factor a few times higher than those established for HOC⁺, $N(HOC^+)/N(HCO^+) = 0.015$, and $X(HOC^+) = 0.45 \times 10^{-11}$ (Gerin et al. 2019).

## 3.2. CO⁺ as a source of HCO⁺

Despite the excellent sensitivity with achieved fractional rms flux density levels $\Delta S_ν/S_ν$ of 0.22% and 0.83% for J1717-3342 and J1744-3116, respectively, CO⁺ was not detected toward these sources. The derived $3\sigma$ upper limits $N(CO^+)/N(HCO^+) \le 0.04$ and $N(CO^+)/N(H_2) < 1.2–1.3 \times 10^{-10}$ are similar for both sources in Table 2, because J1717-3342, with a smaller
<table><caption>**Table 5.** Rates and rate coefficients.</caption><thead><tr><th>Reaction</th><th>Symbol</th><th>Rate coefficient<br>(cm<sup>-3</sup> s<sup>-1</sup>)</th><th>Reference</th></tr></thead><tbody><tr><td>C<sup>+</sup> + OH → CO<sup>+</sup> + H</td><td>k<sub>1</sub></td><td>2.2 × 10<sup>-9</sup>(300/T)<sup>0.107</sup></td><td>Dagdigian (2019)</td></tr><tr><td>CO<sup>+</sup> + H<sub>2</sub> → HCO<sup>+</sup> + H</td><td>k<sub>2</sub></td><td>7.5 × 10<sup>-10</sup></td><td>KIDA</td></tr><tr><td>C<sup>+</sup> + H<sub>2</sub>O → HCO<sup>+</sup> + H</td><td>k<sub>3</sub></td><td>2.1 × 10<sup>-9</sup></td><td>KIDA</td></tr><tr><td>C<sup>+</sup> + H<sub>2</sub>O → HOC<sup>+</sup> + H</td><td>k<sub>4</sub></td><td>2.1 × 10<sup>-9</sup></td><td>KIDA, Martinez et al. (2008)</td></tr><tr><td>HOC<sup>+</sup> + H<sub>2</sub> → HCO<sup>+</sup> + H<sub>2</sub></td><td>k<sub>5</sub></td><td>4 × 10<sup>-10</sup></td><td>Smith et al. (2002)</td></tr><tr><td>HOC<sup>+</sup> + H<sub>2</sub> → HCO<sup>+</sup> + H<sub>2</sub></td><td>k<sub>6</sub></td><td>1 × 10<sup>-11</sup></td><td>KIDA</td></tr><tr><td>HCO<sup>+</sup> + e<sup>-</sup> → CO + H</td><td>α(T)</td><td>1.43 × 10<sup>-5</sup>/T<sup>0.69</sup></td><td>Hamberg et al. (2014)</td></tr><tr><th>Reaction</th><th>Symbol</th><th>Rate<br>(s<sup>-1</sup>)</th><th>Reference</th></tr><tr><td>CO + hν → C + O</td><td>k<sub>d</sub></td><td>2.43 × 10<sup>-10</sup>exp(-3.88A<sub>V</sub>)</td><td>KIDA, Heays et al. (2017)</td></tr><tr><td>CO<sup>+</sup> + hν → C<sup>+</sup> + O</td><td></td><td>1.04 × 10<sup>-10</sup></td><td>KIDA</td></tr></tbody></table>

intervening column density, is much brighter in the continuum. The achieved upper limit on the relative abundance is somewhat above the values seen in dense PDRs where emission lines from CO⁺ are detected and X(CO⁺) reaches a few times 10⁻¹¹. They are below the relative abundances seen in the exceptional case of the young planetary nebula NGC 7027, where N(CO⁺)/N(H₂) = 5 × 10⁻⁹ (Fuente et al. 2003).

It is straightforward to show that the relative abundance of CO⁺ is too small to provide a major source of the observed HCO⁺ in diffuse molecular gas, even independent of the HCO⁺ abundance. In chemical terms, the reaction CO⁺ + H₂ → HCO⁺ + H forming HCO⁺ can be compared with the rate at which HCO⁺ recombines with electrons with rates as given in Table 5.

Equating formation and recombination of HCO⁺, neglecting photodissociation that is four orders of magnitude slower according to the Kinetic Database for Astrochemistry (KIDA; Wakelam et al. 2015):

$$k_2 n(\text{CO}^+) n(\text{H}_2) = \alpha(T) n(\text{HCO}^+) n(e),$$

and when rearranging, we have

$$n(\text{CO}^+) / n(\text{HCO}^+) = (\alpha(T) / k_2) n(e) / n(\text{H}_2),$$

or equivalently

$$n(\text{CO}^+) / n(\text{HCO}^+) = (2 / f_{\text{H}_2}) (\alpha(T) / k_2) x(e),$$

where $f_{\text{H}_2}$ is the fraction of H-nuclei in H₂, $x(e)$ is the electron fraction, $n(\text{H}_2)$ is the density of H₂ molecules, and $n(\text{H})$ is the total density of H-nuclei. Taking $x(e) = 2 \times 10^{-4}$ corresponding to the observed free gas-phase carbon abundance $n(\text{C}^+)/n(\text{H}) = 1.6 \times 10^{-4}$ (Sofia et al. 2004) with a small added contribution from cosmic ray ionization of atomic hydrogen in largely molecular gas, the required amount of CO⁺ needed to restore the HCO⁺ lost to recombination is

$$n(\text{CO}^+) / n(\text{HCO}^+) = (0.6 / f_{\text{H}_2}) \times (40 / T)^{0.69}.$$

This is some 25 times higher than our upper limit if $f_{\text{H}_2} = 0.6$ and $T = 40$ K. We note above that CO emission, even at the low levels shown in Fig. 2, does not arise in regions of very small molecular fraction and at temperatures high enough to trigger the formation of CH⁺: these are only characteristic of regions of strong turbulent energy dissipation.

Alternatively, equating the HCO⁺ formation and recombination rates once again but setting $n(\text{HCO}^+)/n(\text{H}_2) = 3 \times 10^{-9}$ and $x(e) = 2 \times 10^{-4} = 2 \times 10^{-4}(2n(\text{H}_2)/f_{\text{H}_2})$, for the required relative abundance of CO⁺ one derives

$$X(\text{CO}^+) = (1.8 \times 10^{-9} / f_{\text{H}_2}) \times (40 / T)^{0.69}.$$

This is comparable to the CO⁺ relative abundance seen in NGC 7027 (~5 × 10⁻⁹, Fuente et al. 2003) and is at least 15 times above the upper limit achieved here, 1.2 × 10⁻¹⁰, if $T = 40$ K (Table 2). In this simplified analysis, reducing $X(\text{CO}^+)$ below our upper limit would require temperatures of about 2000 K, which correspond to the order of magnitude achieved in the turbulence-driven chemistry, which can only be present in small and intermittent regions of space and time, where other thermal and nonthermal processes are operating. CO⁺ is responsible for producing at most a few percent of the HCO⁺ that recombines with ambient thermal electrons to form CO in diffuse molecular gas.

## 4. The thermal chemistry of CO⁺, HOC⁺, and HCO⁺

Our upper limits show that CO⁺ cannot be responsible for forming the observed CO in diffuse molecular gas, because CO⁺ does not exist in a sufficient quantity to replenish HCO⁺ as it recombines with ambient electrons. However, beyond that, our limits combined with the observed relative abundance $N(\text{HOC}^+)/N(\text{HCO}^+) = 1.4 \pm 0.2 \times 10^{-4}$ (Gerin et al. 2019) can be used to gain further insights into the thermal chemistry that is working in diffuse molecular gas. Most of the reactions that form HCO⁺ also form HOC⁺, yet their observed abundances differ by a factor of 70. To explore this, we extended the thermal CO⁺ and HOC⁺ chemistry discussed by Gerin et al. (2019) to include some hypotheticals corresponding to uncertainties in reaction rates in the KIDA reaction rate database (Wakelam et al. 2015) and molecular abundances (Gerin et al. 2016). As discussed in Gerin et al. (2019), some of the rates of the reactions controlling the formation of CO⁺ and HCO⁺ were not well known, especially the reaction C⁺ + OH → CO⁺ + H/CO + H⁺ and the isomerization reaction destroying HOC⁺ and forming HCO⁺, HOC⁺ + H₂ → HCO⁺ + H₂. The rate of this reaction is reported as $4 \times 10^{-10} \text{ cm}^3 \text{ s}^{-1}$ at 25 and 300 K (Smith et al. 2002), which is about 1/5 the Langevin capture value, and, in this temperature range, it is probably independent of the temperature.
<figure><img src="image_3.png" /><figcaption>**Fig. 3.** Models of thermal chemistry for the relative abundances $N(\text{CO}^+)/N(\text{HCO}^+)$ and $N(\text{HOC}^+)/N(\text{HCO}^+)$ with default rates to the left, and deviations in the middle and right panels. In all panels, results for HOC$^+$ are shown in green and for CO$^+$ in orange. In all cases, $N(\text{HCO}^+)$ is taken to be the observed value $N(\text{HCO}^+) = 3 \times 10^{-9}N(\text{H}_2)$. Left panel: baseline models with the observed upper limit for $N(\text{CO}^+)/N(\text{HCO}^+)$ from this work and the mean $N(\text{HOC}^+)/N(\text{HCO}^+)$ ratio from Gerin et al. (2019) are illustrated schematically as horizontal dashed lines. Middle panel: the solid lines represent the case that the rate constant of $\text{H}_2 + \text{HOC}^+ \rightarrow \text{HOC}^+ + \text{H}_2$ interconversion is taken as $10^{-11}\text{cm}^3\text{s}^{-1}$. The dashed-dotted lines represent model results when the rate constant for the reaction $\text{H}_2 + \text{CO}^+ \rightarrow \text{HOC}^+ + \text{H}$ is set to 0. Right panel: the solid lines show models in which the reaction of $\text{C}^+$ and $\text{H}_2\text{O}$ does not form HOC$^+$ and the dashed-dotted lines show model results when $\text{CH}^+$ is as abundant as OH, $X(\text{CH}^+) = 10^{-7}$.</figcaption></figure>

As originally suggested by Glassgold & Langer (1975, 1976), $\text{CO}^+$ can be formed by the reaction of $\text{C}^+$ with OH, and also by the reaction of O with $\text{CH}^+$. $\text{CO}^+$ is destroyed by recombination with electrons and in chemical reactions with atomic and molecular hydrogen, and photo-dissociated with a free space rate $1 \times 10^{-10} \text{ s}^{-1}$ that is always at least 300 times slower than the other destruction rates in the calculations discussed here. The OH relative abundance to $\text{H}_2$, $X(\text{OH}) = 1.0 \pm 0.14 \times 10^{-7}$ varies only narrowly and is well-determined in local diffuse clouds (Weselak et al. 2009, 2010) such as those probed in this work. The excellent correspondence of the OH and $\text{H}_2\text{O}$ absorption line profiles and of the $\text{H}_2\text{O}$ and HOC$^+$ profiles along long sight lines across the Galactic plane (Wiesemeyer et al. 2016; Gerin et al. 2019) further indicates that this value of the OH abundance holds for the Galactic disk as well. The rate of the reaction of OH and $\text{C}^+$ was recently calculated by Dagdigian (2019), providing accurate information on the rate coefficient and branching ratio over a temperature range 10–1000 K. We approximated the rate constant for $\text{C}^+ + \text{OH}$ forming $\text{CO}^+$ as $k = 2.2 \times 10^{-9} (300/\text{T})^{0.107} \text{cm}^3 \text{s}^{-1}$ (see Table 5). By contrast, the abundance of $\text{CH}^+$ is locally variable, and $\text{CH}^+$ may not co-exist with the other species if it is predominantly formed in the non-thermal chemistry (Godard et al. 2014; Valdivia et al. 2017). We take $X(\text{CH}^+)$ as a free parameter with $X(\text{CH}^+) = 10^{-9}$ by default in the following discussion.

HOC$^+$ can be formed in the reaction of $\text{H}_2$ with $\text{CO}^+$ and in the reaction of $\text{C}^+$ with $\text{H}_2\text{O}$, of which we take the fixed relative abundance $X(\text{H}_2\text{O}) = 2.4 \times 10^{-8}$ from Gerin et al. (2016). We ignored HOC$^+$ formation by the reaction of O with $\text{CH}_3^+$ because the formation of $\text{CH}_3^+$ is an aspect of the nonthermal chemistry. HOC$^+$ recombines with electrons and is susceptible to photodissociation, but most importantly is isomerized to HOC$^+$ in reaction with $\text{H}_2$. The rate constant for the isomerization reaction of $\text{H}_2$ and HOC$^+$ has alternatives in the KIDA database as discussed previously. Here, we considered the two extreme values for the rate constant, $10^{-11} \text{cm}^3 \text{s}^{-1}$ given for 10–280 K and $4.7 \times 10^{-10} \text{cm}^3 \text{s}^{-1}$ given at 305 K.

The $\text{CO}^+$ and HOC$^+$ chemistry is simple when other molecular abundances are held fixed, but it involves the thermal and ionization balance for the temperature-sensitive recombination

rates and the atomic/molecular hydrogen balance for various reactions with $\text{CO}^+$ and HOC$^+$. To address this, we performed calculations similarly to those described earlier to account for the observations of HF and $\text{CF}^+$ (Liszt et al. 2015). In brief, we calculated the self-consistent global atomic/molecular, ionization and thermal equilibrium in gas spheres of constant total hydrogen number and column density, subject to the usual interstellar radiation fields. The hydrogen column density was varied, and molecular abundances were integrated along the central line of sight to form the results that are displayed here. Plotting the results against $N(\text{H}_2)$ removes some model sensitivities; for instance, to the impact parameter about the cloud center used for the line of sight integration.

The results are shown in Fig. 3. To the left is a baseline model using default values (Table 5). The model produces $\text{CO}^+$ relative abundances about 2.5 times smaller than our observed upper limits, and about twice as much HOC$^+$ as observed. Declines in the abundances of $\text{CO}^+$ and HOC$^+$ for very high $N(\text{H}_2)$ arise from the recombination of the free gas-phase carbon from $\text{C}^+$ to neutral carbon and CO. In the middle panel, the solid lines show results when $\text{H}_2$-catalyzed HOC$^+ \rightarrow \text{HOC}^+$ isomerization has the rate constant $10^{-11} \text{cm}^3 \text{s}^{-1}$, a factor of 2.5 above the default value. The green curve for HOC$^+$ shows that the predicted abundance of HOC$^+$ grows even further above what is observed when the HOC$^+ \rightarrow \text{HOC}^+$ interconversion is slow. The orange curves show what happens when the rate constant for HOC$^+$ formation via $\text{H}_2 + \text{CO}^+ \rightarrow \text{HOC}^+ + \text{H}$ is set to 0. The predicted amount of $\text{CO}^+$ grows well above our upper limits, and the abundance of HOC$^+$ falls far below the observed mean. We conclude that HOC$^+$ is forming mostly via $\text{C}^+ + \text{OH}$ followed by $\text{CO}^+ + \text{H}_2$, but much or most of the HOC$^+$ produced by these chemical reactions is lost to isomerization into HOC$^+$.

In the right panel of Fig. 3, the solid curves show model results when the reaction of $\text{C}^+$ and $\text{H}_2\text{O}$ does not form HOC$^+$: all of the observed HOC$^+$ can indeed be made by the reaction of $\text{C}^+$ and OH, and agreement with the observed HOC$^+/\text{HCO}^+$ ratio is somewhat better when the contribution from $\text{H}_2\text{O}$ is neglected. However, $\text{H}_2\text{O}$ is observed to be present with $X(\text{H}_2\text{O}) = 2.4 \times 10^{-8}$ (Gerin et al. 2016) and would be expected to form some amount of HOC$^+$. The default model may predict
an overabundance of CO⁺ or perhaps some aspect of the reaction of C⁺ and H₂O is not fully understood: a slower rate or a branching ratio that favors HCO⁺ would lead to a lower contribution of this reaction to the HOC⁺ formation.

Moreover, CH⁺ may not co-exist with the other species if it is predominantly formed in the nonthermal chemistry (Godard et al. 2014; Valdivia et al. 2017) and X(CH⁺) is negligible (10⁻⁹) in the default model. The final deviation from the default model is represented by the dash-dot lines in the right-hand panel showing results when X(CH⁺) = X(OH) = 10⁻⁷ is large. The reaction of O + CH⁺ becomes an important source of CO⁺ in the thermal chemistry when X(CH⁺) ≥ 5 × 10⁻⁸ and our upper limits on CO⁺ by themselves imply X(CH⁺)/X(OH) ≤ 3. At such high CH⁺ abundances the discrepancy in HOC⁺ grows even larger.

## 5. Conditions for the growth of CO

Balancing the formation of CO through the dissociative recombination of HCO⁺ and the destruction of CO by photodissociation can constrain the physical conditions where CO might form in diffuse molecular gas. With a free-space photodissociation rate $k_d = 2.43 \times 10^{-10}$ s⁻¹ for the Draine radiation field (defined as $G_0 = 1$) and an attenuation scaling as $e^{-\gamma A_V}$, this balance locally implies

$$ \frac{n(\text{CO})}{n(\text{HCO}^+)} = \frac{\alpha(T)x(e)n(\text{H})}{G_0 k_d e^{-\gamma A_V}} $$

For a uniform slab, integrating the attenuation term from 0 up to half the slab extinction,

$$ \frac{N(\text{CO})}{N(\text{HCO}^+)} = \frac{\alpha(T)x(e)n(\text{H})}{G_0 k_d} \frac{0.5\gamma A_V}{1 - e^{-0.5\gamma A_V}} $$

$n(H)/G_0$ can therefore be expressed as

$$ n(H)/G_0 = \frac{N(\text{CO})}{N(\text{HCO}^+)} \frac{k_d(1 - e^{-0.5\gamma A_V})}{\alpha(T)x(e)0.5\gamma A_V} $$

where we approximate $A_V = 3.1E(B-V)$ and take $E(B-V)$ from Green et al. (2019), replacing local quantities with line-of-sight averages. The error bar in $E(B-V)$ is a nominal value 0.02 mag. We take $\gamma = 3.88$ and use the reaction rates and rate coefficients given in Table 5. Figure 4 presents the variation of $n(H)/G_0$ as a function of the reddening $E(B-V)$ for the existing sample of quasars sight lines with HCO⁺ and CO column densities determined in absorption (Liszt & Lucas 1998; Liszt et al. 2019). We used conditions as described in Sect. 2, a uniform temperature of 40 K, and a uniform electron fraction $x(e) = 2 \times 10^{-4}$.

The mean value for the sample is $(n(H)/G_0) \sim 300 \text{ cm}^{-3}$. For a standard interstellar radiation field of $G_0 = 1$, the calculated densities $n(H) \sim 100-300 \text{ cm}^{-3}$ are modest for $E(B-V) \ge 0.5$ mag, even in the absence of any consideration of CO shielding by itself or by H₂. The higher densities $n(H) \sim 1000-2000 \text{ cm}^{-3}$ calculated at $E(B-V) \le 0.2$ mag are higher by a factor of ~5 than those inferred from the brightness of HCO⁺ emission in these directions (Liszt 2020). Thermal pressures $p/k \sim 40 \text{ K} \times 1000 \text{ cm}^{-3}$ are an order of magnitude higher than typical values in the diffuse molecular gas (Jenkins & Tripp 2011) and would more likely represent the conditions under which a turbulence-driven chemistry forms HCO⁺. A lower value of the radiation field toward these sight lines would bring the derived densities closer to the mean values. The data shown in this plot include the sight lines with CO and HCO⁺ detections, which may bias

<figure><img src="image_4.png" /><figcaption>Fig. 4. Number density $n(H)/G$ vs. reddening $E(B-V)$ (Green et al. 2019) derived by equating CO formation and photo-destruction for sight lines with HCO⁺ and CO column densities measured in absorption (Liszt & Lucas 1998; Liszt et al. 2019), as discussed in Sect. 5.</figcaption></figure>

the derived $n(H)/G_0$ to somewhat high values as a higher density implies a faster CO production rate.

The inverse variation of $n(H)/G_0$ with $E(B-V)$ in Fig. 4, mirroring the functional dependence on $A_V$, reflects the fact that the observed values of $N(CO)/N(HCO^+)$ are not correlated with $E(B-V)$ and show a moderate scatter of at most a factor of 1. A majority of the selected sight lines have $N(CO)/N(HCO^+)$ between 1000 and 2000. We also plotted $n(H)/G_0$ against the molecular fraction $f_{H_2} = 2N(H_2)/N(H)$ taking $N(H_2) = N(HCO^+)/3 \times 10^{-9}$ and $N(H) = a \times E(B-V)$. Requiring $f_{H_2} \le 1$ in all cases implies $a \ge 8 \times 10^{21} \text{ H} - \text{nuclei cm}^{-2} \text{ mag}^{-1}$, but there is no trend for $n(H)/G_0$ to vary with $f_{H_2}$ calculated in this way.

## 6. Summary and conclusions

We sought to test the origin of CO in diffuse molecular gas within the generally accepted framework that CO is predominantly formed through the dissociative electron recombination of HCO⁺. HCO⁺ is widely observed in diffuse molecular gas, with a relative abundance $N(HCO^+)/N(H_2) = 3 \times 10^{-9}$.

In a quiescent, thermal, cosmic-ray, and UV-driven chemistry, HCO⁺ is predominantly formed in the reactions C⁺ + OH → CO⁺ + H and CO⁺ + H₂ → HCO⁺ + H. Using Cycle 7 ALMA Band-6 observations at 236 GHz (Table 1), we searched for CO⁺ absorption toward two bright compact extragalactic millimeter-wave continuum sources seen at low latitude in the inner Galaxy and known to be occulted by high column densities of HCO⁺-bearing diffuse molecular gas (see Sect. 2 and Table 2). We failed to detect CO⁺ at levels sufficient to demonstrate in Sect. 3 that the reaction CO⁺ + H₂ → HCO⁺ + H cannot replenish more than a few percent of the observed HCO⁺ that recombines. The non-detection of CO⁺ confirms that HCO⁺ is mainly produced in the reaction between oxygen and carbon hydrides CH₂ or CH₃ induced by suprathermal processes, while CO⁺ and HOC⁺ result from reactions between C⁺ and OH and H₂O in the thermal chemistry that occurs in quiescent diffuse molecular gas. The observed CO⁺ and HOC⁺ abundances relative to H₂ set an upper limit on the rate coefficient of these reactions.

In Sect. 4, we have explored the coupled thermal chemistries of CO⁺ and HOC⁺ given our observational upper limits on CO⁺ and existing observations of HOC⁺ showing
$N(HOC^+)/N(HCO^+)$ = 0.015 ± 0.003. A baseline cosmic-ray and UV-driven thermal chemical model overproduced $HOC^+$ by a factor of two, but by as much as a factor of six if the isomerization reaction $HOC^+ + H_2 \rightarrow HCO^+ + H_2$ does proceed rapidly. The option of a low rate for this reaction in the KIDA database should probably be ignored. The predicted $HOC^+$ abundance is in better agreement with the observations if the reaction of $C^+$ and $H_2O$ preferentially forms $HCO^+$, rather than $HOC^+$. We placed a weak limit on the relative abundance of $CH^+$ in quiescent gas, $N(CH^+)/N(OH) \le 3$ or $X(CH^+)/X(H_2) \le 3 \times 10^{-7}$, such that the reaction $O + CH^+$ does not become a significant source of $CO^+$, substantially overproducing it.

In Sect. 5, we have derived characteristic number densities $n(H)$ at which CO would form in diffuse molecular gas with observed abundances. For $E(B-V) \ge 0.5$ mag, the densities are modest, $n(H) \sim 100-300 \text{ cm}^{-3}$, but for CO-bright regions of low extinction $E(B-V) \sim 0.1-0.2$ mag they are an order of magnitude higher $n(H) \sim 1-2 \times 10^3 \text{ cm}^{-3}$, and such heavily over-pressured gas is more properly considered in the context of a turbulent dissipation chemistry. The presence of CO in sight lines of low reddening could also indicate a lower radiation field intensity.

*Acknowledgements.* The authors thank the referee Helmut Wiesemeyer for his thoughtful comments and suggestions. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2019.1.00120.S.. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This work was supported by the Programme National "Physique et Chimie du Milieu Interstellaire" (PCMI) of CNRS/INSU with INC/INP co-funded by CEA and CNES.

## References

Agúndez, M., Goicoechea, J. R., Cernicharo, J., Faure, A., & Roueff, E. 2010, ApJ, 713, 662

Bitran, M., Alvarez, H., Bronfman, L., May, J., & Thaddeus, P. 1997, ApJS, 125, 99

Black, J. H., & Dalgarno, A. 1977, ApJS, 34, 405

Dagdigian, P. J. 2019, J. Chem. Phys., 151, 054306

Endres, C. P., Schlemmer, S., Schilke, P., Stutzki, J., & Müller, H. S. P. 2016, J. Mol. Spectr., 327, 95

Fuente, A., Rodriguez-Franco, A., Garcia-Burillo, S., Martín-Pintado, J., & Black, J. H. 2003, A&A, 406, 899

Gerin, M., & Liszt, H. 2017, A&A, 600, A48

Gerin, M., Ruaud, M., Goicoechea, J. R., et al. 2015, A&A, 573, A30

Gerin, M., Neufeld, D. A., & Goicoechea, J. R. 2016, ARA&A, 54, 181

Gerin, M., Liszt, H., Neufeld, D., et al. 2019, A&A, 622, A26

Giannetti, A., Wyrowski, F., Brand, J., et al. 2014, A&A, 570, A65

Glassgold, A. E., & Langer, W. D. 1975, ApJ, 197, 347

Glassgold, A. E., & Langer, W. D. 1976, ApJ, 206, 85

Godard, B., Falgarone, E., & Pineau des Forêts, G. 2014, A&A, 570, A27

Goicoechea, J. R., Cuadrado, S., Pety, J., et al. 2017, A&A, 601, L9

Goicoechea, J. R., Santa-Maria, M. G., Bron, E., et al. 2019, A&A, 622, A91

Goldsmith, P. F., Pineda, J. L., Neufeld, D. A., et al. 2018, ApJ, 856, 96

Green, G. M., Schlafly, E., Zucker, C., Speagle, J. S., & Finkbeiner, D. 2019, ApJ, 887, 93

Hamberg, M., Kashperkat, I., Thomas, R. D., et al. 2014, J. Phys. Chem. A., 118, 6034

Heays, A. N., Bosman, A. D., & van Dishoeck, E. F. 2017, A&A, 602, A105

Jenkins, E. B., & Tripp, T. M. 2011, ApJ, 734, 65

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775

Keene, J., Schilke, P., Kooi, J., et al. 1998, ApJ, 494, L107

Latter, W. B., Walker, C. K., & Maloney, P. R. 1993, ApJ, 419, L97

Lesaffre, P., Todorov, P., Levrier, F., et al. 2020, MNRAS, 495, 816

Liszt, H. S. 2017, ApJ, 835, 138

Liszt, H. S. 2020, ApJ, 897, 104

Liszt, H., & Gerin, M. 2018, A&A, 610, A49

Liszt, H. S., & Lucas, R. 1998, A&A, 339, 561

Liszt, H. S., & Lucas, R. 2000, A&A, 355, 333

Liszt, H. S., Pety, J., & Lucas, R. 2010, A&A, 518, A45

Liszt, H., Lucas, R., Pety, J., & Gerin, M. 2014a, IAU Symp., 297, 163

Liszt, H. S., Pety, J., Gerin, M., & Lucas, R. 2014b, A&A, 564, A64

Liszt, H. S., Guzmán, V. V., Pety, J., et al. 2015, A&A, 579, A12

Liszt, H., Gerin, M., Beasley, A., & Pety, J. 2018, ApJ, 856, 151

Liszt, H., Gerin, M., & Grenier, I. 2019, A&A, 627, A95

Lucas, R., & Liszt, H. 1996, A&A, 307, 237

Martinez, Jr. O., Betts, N. B., Villano, S. M., et al. 2008, ApJ, 686, 1486

Milam, S. N., Savage, C., Brewster, M. A., Ziurys, L. M., & Wyckoff, S. 2005, ApJ, 634, 1126

Müller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 2001, A&A, 370, L49

Müller, H. S. P., Schlöder, F., Stutzki, J., & Winnewisser, G. 2005, J. Mol. Struct., 742, 215

Riquelme, D., Bronfman, L., Mauersberger, R., et al. 2018, A&A, 610, A43

Sastry, K. V. L. N., Helminger, P., Herbst, E., & De Lucia, F. C. 1981, ApJ, 250, L91

Savage, B. D., Drake, J. F., Budich, W., & Bohlin, R. C. 1977, ApJ, 216, 291

Sheffer, Y., Rogers, M., Federman, S. R., et al. 2008, ApJ, 687, 1075

Smith, M. A., Schlemmer, S., von Richthofen, J., & Gerlich, D. 2002, ApJ, 578, L87

Snow, T. P., & McCall, B. J. 2006, ARA&A, 44, 367

Sofia, U. J., Lauroesch, J. T., Meyer, D. M., & Cartledge, S. I. B. 2004, ApJ, 605, 272

Treviño-Morales, S. P., Fuente, A., Sánchez-Monge, Á., et al. 2016, A&A, 593, L12

Valdivia, V., Godard, B., Hennebelle, P., et al. 2017, A&A, 600, A114

Van Dishoeck, E. F., & Black, J. H. 1986, ApJS, 62, 109

van Dishoeck, E. F., & Black, J. H. 1988, ApJ, 334, 771

Visser, R., van Dishoeck, E. F., & Black, J. H. 2009, A&A, 503, 323

Wakelam, V., Loison, J. C., Herbst, E., et al. 2015, ApJS, 217, 20

Weselak, T., Galazutdinov, G., Beletsky, Y., & Krelowski, J. 2009, A&A, 499, 783

Weselak, T., Galazutdinov, G. A., Beletsky, Y., & Krelowski, J. 2010, MNRAS, 402, 1991

Wiesemeyer, H., Güsten, R., Heyminck, S., et al. 2016, A&A, 585, A76

Wilson, T. L., & Rood, R. 1994, ARA&A, 32, 191