We have observed the massive protostar AFGL 2136 IRS 1 in multiple wavelength windows in the near-to-mid-infrared at high ($\sim3$ km s$^{-1}$) spectral resolution using VLT+CRIRES, SOFIA+EXES, and Gemini North+TEXES. There is an abundance of H$_2$O absorption lines from the $\nu_1$ and $\nu_3$ vibrational bands at 2.7 $\mu$m, from the $\nu_2$ vibrational band at 6.1 $\mu$m, and from pure rotational transitions near 10-13 $\mu$m. Analysis of state-specific column densities derived from the resolved absorption features reveals that an isothermal absorbing slab model is incapable of explaining the relative depths of different absorption features. In particular, the strongest absorption features are much weaker than expected, indicating optical depth effects resulting from the absorbing gas being well-mixed with the warm dust that serves as the "background" continuum source at all observed wavelengths. The velocity at which the strongest H$_2$O absorption occurs coincides with the velocity centroid along the minor axis of the compact disk in Keplerian rotation recently observed in H$_2$O emission with ALMA. We postulate that the warm regions of this dust disk dominate the continuum emission at near-to-mid infrared wavelengths, and that H$_2$O and several other molecules observed in absorption are probing this disk. Absorption line profiles are not symmetric, possibly indicating that the warm dust in the disk that produces the infrared continuum has a non-uniform distribution similar to the substructure observed in 1.3 mm continuum emission.

In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. In this review we summarise in a coherent picture the main results obtained by observations and by theoretical models of propagation and generation of cosmic rays, from the smallest scales of protostars and circumstellar discs, to young stellar clusters, up to Galactic and extragalactic scales. We also discuss the new fields that will be explored in the near future thanks to new generation instruments, such as: CTA, for the $\gamma$-ray emission from high-mass protostars; SKA and precursors, for the synchrotron emission at different scales; and ELT/HIRES, JWST, and ARIEL, for the impact of cosmic rays on exoplanetary atmospheres and habitability.

The molecular gas in the Central Molecular Zone (CMZ) of the Galaxy has been studied using infrared absorption spectra of H$_3^+$ lines at 3.5-4.0 $\mu$m and CO lines near 2.34 $\mu$m. In addition to the previously reported spectra of these lines toward 8 stars located within 30 pc of Sgr A$^\ast$, there are now spectra toward $\sim$30 bright stars located from 140 pc west to 120 pc east of Sgr A$^\ast$. The spectra show the presence of warm ($T\sim 200$ K) and diffuse ($n < 100 $cm$^{-3}$) gas with $N$(H$_3^+$) $\sim 3 \times 10^{15} $cm$^{-2}$ on majority of sightlines. Instead of our previous analysis in which only electrons from photoionization of carbon atoms were considered, we have developed a simple model calculation in which the cosmic ray ionization of H$_2$ and H are also taken into account. We conclude: (1) Warm and diffuse gas dominates the volume of the CMZ. The volume filling factor of dense gas must be much less than 0.1 and the CMZ is not as opaque as previously considered. The X-ray emitting ultra-hot 10$^8 $K plasma, which some thought to dominate the CMZ, does not exist over extended regions. (2) The cosmic ray ionization rate is $\zeta \sim 2 \times 10^{-14} $s$^{-1}$, higher than in Galactic dense clouds and diffuse clouds by factors of $\sim$1000 and $\sim$100, respectively. If the equipartition law stands, this suggests a pervading magnetic field on the order of $\sim$100 $\mu$G.

Sulfur has been observed to be severely depleted in dense clouds leading to uncertainty in the molecules that contain it and the chemistry behind their evolution. Here, we aim to shed light on the sulfur chemistry in young stellar objects (YSOs) by using high-resolution infrared spectroscopy of absorption by the $\nu_3$ rovibrational band of SO$_2$ obtained with the Echelon-Cross-Echelle Spectrograph on the Stratospheric Observatory for Infrared Astronomy. Using local thermodynamic equilibrium models we derive physical parameters for the SO$_2$ gas in the massive YSO MonR2 IRS3. This yields a SO$_2$/$\mathrm{H}$ abundance lower limit of $5.6\pm0.5\times10^{-7}$, or $>\!4\%$ of the cosmic sulfur budget, and an intrinsic line width (Doppler parameter) of $b<3.20\;\mathrm{km\;s}^{-1}$. The small line widths and high temperature ($T_\mathrm{ex}=234\pm15\;\mathrm{K}$) locate the gas in a relatively quiescent region near the YSO, presumably in the hot core where ices have evaporated. This sublimation unlocks a volatile sulfur reservoir (e.g., sulfur allotropes as detected abundantly in comet 67P/Churyumov--Gerasimenko), which is followed by SO$_2$ formation by warm, dense gas-phase chemistry. The narrowness of the lines makes formation of SO$_2$ from sulfur sputtered off grains in shocks less likely toward MonR2 IRS3.

We have performed a 5-8 $\mu$m spectral line survey of the hot molecular core associated with the massive protostar AFGL 2591, using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We have supplemented these data with a ground based study in the atmospheric M band around 4.5 $\mu$m using the iSHELL instrument on the Infrared Telescope Facility (IRTF), and the full N band window from 8-13 $\mu$m using the Texas Echelon Cross Echelle Spectrograph (TEXES) on the IRTF. Here we present the first detection of ro-vibrational transitions of CS in this source. The absorption lines are centred on average around -10 kms$^{-1}$ and the line widths of CS compare well with the hot component of $^{13}$CO (around 10 kms$^{-1}$). Temperatures for CS, hot $^{13}$CO and $^{12}$CO v=1-2 agree well and are around 700 K. We derive a CS abundance of 8$\times$10$^{-3}$ and 2$\times$10$^{-6}$ with respect to CO and H$_2$ respectively. This enhanced CS abundance with respect to the surrounding cloud (1$\times$10$^{-8}$) may reflect sublimation of H$_2$S ice followed by gas-phase reactions to form CS. Transitions are in LTE and we derive a density of $>$10$^7$ cm$^{-3}$, which corresponds to an absorbing region of $<$0.04$''$. EXES observations of CS are likely to probe deeply into the hot core, to the base of the outflow. Submillimeter and infrared observations trace different components of the hot core as revealed by the difference in systemic velocities, line widths and temperatures, as well as the CS abundance.

Orion BN has a large proper motion and radial velocity with respect to the gas and other stars in the region where it is presumed to have formed. Multiple dynamical interaction scenarios have been proposed to explain this motion. In one case BN is thought to have interacted with stars in the Trapezium cluster, while in another it is thought to have interacted with source I while deeply embedded in molecular gas. If there is dense gas that has been retained in close proximity to BN, it may be evidence that the latter scenario is favored. We observed BN at high spectral resolution in three windows near 6 $\mu$m using SOFIA/EXES targeting the $\nu_2$ vibrational band of H$_2$O. Absorption from only three transitions of H$_2$O is detected, and through kinematic analysis is associated with cool, dense foreground gas, not BN itself. We find no evidence for H$_2$O absorption or emission at the systemic velocity of BN.

Cosmic rays are predominantly accelerated in shocks associated with star formation such as supernova remnants and stellar wind bubbles, so the cosmic-ray flux and thus cosmic-ray ionization rate, $\zeta_{\rm H}$, should correlate with the star-formation rate in a galaxy. Submillimeter bright galaxies (SMGs) are some of the most prolific star forming galaxies in the Universe, and gravitationally lensed SMGs provide bright continuum sources suitable for absorption line studies. Abundances of OH$^+$ and H$_2$O$^+$ are useful for inferring $\zeta_{\rm H}$ when combined with chemical models, and have been used for this purpose within the Milky Way. At redshifts $z\gtrsim2$ transitions out of the ground rotational states of OH$^+$ and H$_2$O$^+$ are observable with ALMA, and we present observations of both molecules in absorption toward the lensed SMGs SMM J2135$-$0102 and SDP 17b. These detections enable an exploration of $\zeta_{\rm H}$ in galaxies with extreme star formation and high supernova rates, both of which should significantly enhance cosmic-ray production. The observed OH$^+$ and H$_2$O$^+$ absorption is thought to arise in massive, extended halos of cool, diffuse gas that surround these galaxies. Using a chemical model designed to focus on the reaction network important to both species, we infer cosmic-ray ionization rates of $\zeta_{\rm H}\sim10^{-16}$-$10^{-14}$ s$^{-1}$ in these extended gaseous halos. As our estimates come from gas that is far away from the sites of cosmic-ray acceleration, they imply that cosmic-ray ionization rates in the compact regions where star formation occurs in these galaxies are orders of magnitude higher.

Recent submillimeter and far-infrared wavelength observations of absorption in the rotational ground-state lines of various simple molecules against distant Galactic continuum sources have opened the possibility of studying the chemistry of diffuse molecular clouds throughout the Milky Way. In order to calculate abundances, the column densities of molecular and atomic hydrogen, HI, must be known. We aim at determining the atomic hydrogen column densities for diffuse clouds located on the sight lines toward a sample of prominent high-mass star-forming regions that were intensely studied with the HIFI instrument onboard Herschel. Based on Jansky Very Large Array data, we employ the 21 cm HI absorption-line technique to construct profiles of the HI opacity versus radial velocity toward our target sources. These profiles are combined with lower resolution archival data of extended HI emission to calculate the HI column densities of the individual clouds along the sight lines. We employ Bayesian inference to estimate the uncertainties of the derived quantities. Our study delivers reliable estimates of the atomic hydrogen column density for a large number of diffuse molecular clouds at various Galactocentric distances. Together with column densities of molecular hydrogen derived from its surrogates observed with HIFI, the measurements can be used to characterize the clouds and investigate the dependence of their chemistry on the molecular fraction, for example.

The relative populations in rotational transitions of CO can be useful for inferring gas conditions and excitation mechanisms at work in the interstellar medium. We present CO emission lines from rotational transitions observed with Herschel/HIFI in the star-forming cores Orion S, Orion KL, Sgr B2(M), and W49N. Integrated line fluxes from these observations are combined with those from Herschel/PACS observations of the same sources to construct CO spectral line energy distributions (SLEDs) from $5\leq J_u\leq48$. These CO SLEDs are compared to those reported in other galaxies, with the intention of empirically determining which mechanisms dominate excitation in such systems. We find that CO SLEDs in Galactic star-forming cores cannot be used to reproduce those observed in other galaxies, although the discrepancies arise primarily as a result of beam filling factors. The much larger regions sampled by the Herschel beams at distances of several Mpc contain significant amounts of cooler gas which dominate the extragalactic CO SLEDs, in contrast to observations of Galactic star-forming regions which are focused specifically on cores containing primarily hot molecular gas.

Hydrides are simple compounds containing one or a few hydrogen atoms bonded to a heavier atom. They are fundamental precursor molecules in cosmic chemistry and many hydride ions have become observable in high quality for the first time thanks to the Herschel Space Observatory. Ionized hydrides, such as CH+ and OH+, and also HCO+ that affect the chemistry of molecules such as water, provide complementary information on irradiation by far UV (FUV) or X-rays and gas temperature. The targeted lines of CH+, OH+, H2O+, C+ and CH are detected mostly in blue-shifted absorption. H3O+ and SH+ are detected in emission and only toward some high-mass objects. The observed line parameters and correlations suggest two different origins, related to gas entrained by the outflows and to the circumstellar envelope. The column density ratios of CH+/OH+ are estimated from chemical slab models, assuming that the H2 density is given by the specific density model of each object at the beam radius. For the low-mass YSOs the observed ratio can be reproduced for an FUV flux of 2-400 times the ISRF at the location of the molecules. In two high-mass objects, the UV flux is 20-200 times the ISRF derived from absorption lines, and 300-600 ISRF using emission lines. If the FUV flux required for low-mass objects originates at the central protostar, a substantial FUV luminosity, up to 1.5 L_sun, is required. There is no molecular evidence for X-ray induced chemistry in the low-mass objects on the observed scales of a few 1000 AU. For high-mass regions, the FUV flux required to produce the observed molecular ratios is smaller than the unattenuated flux expected from the central object(s) at the Herschel beam radius. This is consistent with an FUV flux reduced by circumstellar extinction or by bloating of the protostar.

We report additional detections of the chloronium molecular ion, H$_2$Cl$^+$, toward four bright submillimeter continuum sources: G29.96, W49N, W51, and W3(OH). With the use of the HIFI instrument on the Herschel Space Observatory, we observed the $2_{12}-1_{01}$ transition of ortho-H$_2^{35}$Cl$^+$ at 781.627 GHz in absorption toward all four sources. Much of the detected absorption arises in diffuse foreground clouds that are unassociated with the background continuum sources and in which our best estimates of the $N({\rm H_2Cl^+})/N({\rm H})$ ratio lie in the range $(0.9 - 4.8) \times 10^{-9}$. These chloronium abundances relative to atomic hydrogen can exceed the predictions of current astrochemical models by up to a factor of 5. Toward W49N, we have also detected the $2_{12}-1_{01}$ transition of ortho-H$_2^{37}$Cl$^+$ at 780.053 GHz and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ at 485.418 GHz. These observations imply $\rm H_2^{35}Cl^+/H_2^{37}Cl^+$ column density ratios that are consistent with the solar system $^{35}$Cl/$^{37}$Cl isotopic ratio of 3.1, and chloronium ortho-to-para ratios consistent with 3, the ratio of spin statistical weights.

We present high spectral resolution (~3 km/s) observations of the nu_2 ro-vibrational band of H2O in the 6.086--6.135 micron range toward the massive protostar AFGL 2591 using the Echelon-Cross-Echelle Spectrograph (EXES) on the Stratospheric Observatory for Infrared Astronomy (SOFIA). Ten absorption features are detected in total, with seven caused by transitions in the nu_2 band of H2O, two by transitions in the first vibrationally excited nu_2 band of H2O, and one by a transition in the nu_2 band of H2{18}O. Among the detected transitions is the nu_2 1(1,1)--0(0,0) line which probes the lowest lying rotational level of para-H2O. The stronger transitions appear to be optically thick, but reach maximum absorption at a depth of about 25%, suggesting that the background source is only partially covered by the absorbing gas, or that the absorption arises within the 6 micron emitting photosphere. Assuming a covering fraction of 25%, the H2O column density and rotational temperature that best fit the observed absorption lines are N(H2O)=(1.3+-0.3)*10^{19} cm^{-2} and T=640+-80 K.

In diffuse interstellar clouds the chemistry that leads to the formation of the oxygen bearing ions OH+, H2O+, and H3O+ begins with the ionization of atomic hydrogen by cosmic rays, and continues through subsequent hydrogen abstraction reactions involving H2. Given these reaction pathways, the observed abundances of these molecules are useful in constraining both the total cosmic-ray ionization rate of atomic hydrogen (zeta_H) and molecular hydrogen fraction, f(H2). We present observations targeting transitions of OH+, H2O+, and H3O+ made with the Herschel Space Observatory along 20 Galactic sight lines toward bright submillimeter continuum sources. Both OH+ and H2O+ are detected in absorption in multiple velocity components along every sight line, but H3O+ is only detected along 7 sight lines. From the molecular abundances we compute f(H2) in multiple distinct components along each line of sight, and find a Gaussian distribution with mean and standard deviation 0.042+-0.018. This confirms previous findings that OH+ and H2O+ primarily reside in gas with low H2 fractions. We also infer zeta_H throughout our sample, and find a log-normal distribution with mean log(zeta_H)=-15.75, (zeta_H=1.78x10^-16 s^-1), and standard deviation 0.29 for gas within the Galactic disk, but outside of the Galactic center. This is in good agreement with the mean and distribution of cosmic-ray ionization rates previously inferred from H3+ observations. Ionization rates in the Galactic center tend to be 10--100 times larger than found in the Galactic disk, also in accord with prior studies.

The chemistry in the diffuse interstellar medium initiates the gradual increase of molecular complexity during the life cycle of matter. A key molecule that enables build-up of new molecular bonds and new molecules via proton-donation is H3+. Its evolution is tightly related to molecular hydrogen and thought to be well understood. However, recent observations of ortho and para lines of H2 and H3+ in the diffuse ISM showed a puzzling discrepancy in nuclear spin excitation temperatures and populations between these two key species. H3+, unlike H2, seems to be out of thermal equilibrium, contrary to the predictions of modern astrochemical models. We conduct the first time-dependent modeling of the para-fractions of H2 and H3+ in the diffuse ISM and compare our results to a set of line-of-sight observations, including new measurements presented in this study. We isolate a set of key reactions for H3+ and find that the destruction of the lowest rotational states of H3+ by dissociative recombination largely control its ortho/para ratio. A plausible agreement with observations cannot be achieved unless a ratio larger than 1:5 for the destruction of (1,1)- and (1,0)-states of H3+ is assumed. Additionally, an increased CR ionization rate to 10(-15) 1/s further improves the fit whereas variations of other individual physical parameters, such as density and chemical age, have only a minor effect on the predicted ortho/para ratios. Thus our study calls for new laboratory measurements of the dissociative recombination rate and branching ratio of the key ion H3+ under interstellar conditions.

We describe the assignment of a previously unidentified interstellar absorption line to ArH$^+$ and discuss its relevance in the context of hydride absorption in diffuse gas with a low H$_2$ fraction. The column densities along several lines of sight are determined and discussd in the framework of chemical models. The column densities of ArH$^+$ are compared to those of other species, tracing interstellar medium (ISM) components with different H$_2$ abundances. Chemical models are constructed, taking UV radiation and cosmic ray ionization into account. Due to the detection of two isotopologues, $^{36}$ArH$^+$ and $^{38}$ArH$^+$, we are confident about the carrier assignment to ArH$^+$. NeH$^+$ is not detected with a limit of [NeH$^+$]/[ArH$^+$] $\le$ 0.1. The derived column densities agree well with the predictions of chemical models. ArH$^+$ is a unique tracer of gas with a fractional H$_2$ abundance of $10^{-4}- 10^{-3}$ and shows little correlation with H$_2$O$^+$, which traces gas with a fractional H$_2$ abundance of $\approx $0.1. A careful analysis of variations in the ArH$^+$, OH$^+$, H$_2$O$^+$ and HF column densities promises to be a faithful tracer of the distribution of the H$_2$ fractional abundance, providing unique information on a poorly known phase in the cycle of interstellar matter, its transition from atomic diffuse gas to dense molecular gas traced by CO emission. Abundances of these species put strong observational constraints upon magnetohydrodynamical (MHD) simulations of the interstellar medium, and potentially could evolve into a tool to characterize the ISM. Paradoxically, the ArH$^+$ molecule is a better tracer of \new{almost} purely atomic hydrogen gas than H{\sc i} itself, since H{\sc i} can also be present in gas with a significant molecular content, but ArH$^+$ singles out gas that is $>99.9$\% atomic.

We have obtained improved spectra of key fundamental band lines of H3+, R(1,1)l, R(3,3)l, and R(2,2)l, and ro-vibrational transitions of CO on sightlines toward the luminous infrared sources GCIRS 3 and GCIRS 1W, each located in the Central Cluster of the Galactic center within several arcseconds of Sgr A*. The spectra reveal absorption occurring in three kinds of gaseous environments: (1) cold dense and diffuse gas associated with foreground spiral/lateral arms; (2) warm and diffuse gas absorbing over a wide and mostly negative velocity range, which appears to fill a significant fraction of the Galaxy's Central Molecular Zone (CMZ); and (3) warm, dense and compact clouds with velocities near +50 km s^-1 probably within 1-2 pc of the center. The absorptions by the first two cloud types are nearly identical for all the sources in the Central Cluster, and are similar to those previously observed on sightlines from Sgr A* to 30 pc east of it. Cloud type (3), which has only been observed toward the Central Cluster, shows distinct differences between the sightlines to GCIRS 3 and GCIRS 1W, which are separated on the sky by only 0.33 pc in projection. We identify this material as part of an inward extension of the Circumnuclear Disk previously known from HCN mapping. Lower limits on the products of the hydrogen ionization rate zeta and the path length L are 2.3 x 10^5 cm s^-1 and 1.5 x 10^3 cm s^-1 for the warm and diffuse CMZ gas and for the warm and dense clouds in the core, respectively. The limits indicate that the ionization rates in these regions are well above 10^-15 s^-1.

We would like to understand the chemistry of dense clouds and their hot cores more quantitatively by obtaining more complete knowledge of the chemical species present in them. We have obtained high-resolution infrared absorption spectroscopy at 3-4 um toward the bright infrared source CRL 2136. The fundamental vibration-rotation band of HCl has been detected within a dense cloud for the first time. The HCl is probably located in the warm compact circumstellar envelope or disk of CRL 2136. The fractional abundance of HCl is (4.9-8.7)e-8, indicating that approximately 20 % of the elemental chlorine is in gaseous HCl. The kinetic temperature of the absorbing gas is 250 K, half the value determined from infrared spectroscopy of 13CO and water. The percentage of chlorine in HCl is approximately that expected for gas at this temperature. The reason for the difference in temperatures between the various molecular species is unknown.

We present ground-based observations of the nu_1 and nu_3 fundamental bands of H2O toward the massive protostar AFGL 2136 IRS 1, identifying absorption features due to 47 different ro-vibrational transitions between 2.468 micron and 2.561 micron. Analysis of these features indicates the absorption arises in warm (T=506+-25 K), very dense (n(H2)>5*10^9 cm^-3) gas, suggesting an origin close to the central protostar. The total column density of warm water is estimated to be N(H2O)=(1.02+-0.02)*10^19 cm^-2, giving a relative abundance of N(H2O)/N(H2)~10^-4. Our study represents the first extensive use of water vapor absorption lines in the near infrared, and demonstrates the utility of such observations in deriving physical parameters.

We report observations and analysis of infrared spectra of H3+ and CO lines in the Galactic center, within a few parsecs of the central black hole, Sgr A*. We find a cosmic ray ionization rate typically an order of magnitude higher than outside the Galactic center. Notwithstanding, the elevated cosmic ray ionization rate is 4 orders of magnitude too short to match the proton energy spectrum as inferred from the recent discovery of the TeV gamma-ray source in the vicinity of Sgr A*.

We present a 52-671um spectral scan toward SgrA* taken with the PACS and SPIRE spectrometers onboard Herschel. The achieved angular resolution allows us to separate, for the first time at far-IR wavelengths, the emission toward the central cavity (gas in the inner central parsec of the galaxy) from that of the surrounding circum-nuclear disk. The spectrum toward SgrA* is dominated by strong [OIII], [OI], [CII], [NIII], [NII], and [CI] fine structure lines (in decreasing order of luminosity) arising in gas irradiated by UV-photons from the central stellar cluster. In addition, rotationally excited lines of 12CO (from J=4-3 to 24-23), 13CO, H2O, OH, H3O+, HCO+ and HCN, as well as ground-state absorption lines of OH+, H2O+, H3O+, CH+, H2O, OH, HF, CH and NH are detected. The excitation of the 12CO ladder is consistent with a hot isothermal component at Tk ~ 10^{3.1} K and n(H2)< 10^4 cm^{-3}. It is also consistent with a distribution of temperature components at higher density with most CO at Tk<300 K. The detected molecular features suggest that, at present, neither very enhanced X-ray, nor cosmic-ray fluxes play a dominant role in the heating of the hot molecular gas. The hot CO component (either the bulk of the CO column or just a small fraction depending on the above scenario) results from a combination of UV- and shock-driven heating. If irradiated dense clumps/clouds do not exist, shocks likely dominate the heating of the hot molecular gas. This is consistent with the high-velocity gas detected toward SgrA*.

We report the first detection of the v=1--0, R(0) ro-vibrational transition of HF at 2.499385 microns arising from interstellar gas. The line is seen in absorption toward 3 background sources---HD 154368, Elias 29, and AFGL 2136 IRS 1---all of which have reported H2 column densities determined from observations of H2. This allows for the first direct determination of the HF/H2 abundance ratio. We find values of N(HF)/N(H2)=1.15*10^-8 and 0.69*10^-8 for HD 154368 and Elias 29, respectively. The sight line toward AFGL 2136 IRS 1 also shows absorption from the v=1--0, R(1) transition of HF, indicating warm, dense (n_H > 10^9 cm^-3) gas, likely very close to the central protostar. Ascribing portions of the HF absorption to warm and cold gas, we find N(HF)/N(H2)=(1.7--2.9)*10^-8 and (0.33--0.58)*10^-8 for the two components, respectively. Except for the warm component toward AFGL 2136 IRS 1, all observed HF/H2 ratios are well below N(HF)/N(H2)=3.6*10^-8, the value predicted if nearly all gas-phase fluorine is in the form of HF. Models of fluorine chemistry that account for depletion onto grains are able to reproduce the results toward HD 154368, but not in the cold, dense gas toward AFGL 2136 IRS 1 and Elias 29. Most likely, some combination of simplifying assumptions made in the chemical models are responsible for these discrepancies.

Absorption lines from the molecules OH+, H2O+, and H3+ have been observed in a diffuse molecular cloud along a line of sight near W51 IRS2. We present the first chemical analysis that combines the information provided by all three of these species. Together, OH+ and H2O+ are used to determine the molecular hydrogen fraction in the outskirts of the observed cloud, as well as the cosmic-ray ionization rate of atomic hydrogen. H3+ is used to infer the cosmic-ray ionization rate of H2 in the molecular interior of the cloud, which we find to be zeta_2=(4.8+-3.4)x10^-16 per second. Combining the results from all three species we find an efficiency factor---defined as the ratio of the formation rate of OH+ to the cosmic-ray ionization rate of H---of epsilon=0.07+-0.04, much lower than predicted by chemical models. This is an important step in the future use of OH+ and H2O+ on their own as tracers of the cosmic-ray ionization rate.

Using the Herschel Space Observatory's Heterodyne Instrument for the Far-Infrared (HIFI), we have observed para-chloronium (H2Cl+) toward six sources in the Galaxy. We detected interstellar chloronium absorption in foreground molecular clouds along the sight-lines to the bright submillimeter continuum sources Sgr A (+50 km/s cloud) and W31C. Both the para-H2-35Cl+ and para-H2-37Cl+ isotopologues were detected, through observations of their 1(11)-0(00) transitions at rest frequencies of 485.42 and 484.23 GHz, respectively. For an assumed ortho-to-para ratio of 3, the observed optical depths imply that chloronium accounts for ~ 4 - 12% of chlorine nuclei in the gas phase. We detected interstellar chloronium emission from two sources in the Orion Molecular Cloud 1: the Orion Bar photodissociation region and the Orion South condensation. For an assumed ortho-to-para ratio of 3 for chloronium, the observed emission line fluxes imply total beam-averaged column densities of ~ 2.0E+13 cm-2 and ~ 1.2E+13 cm-2, respectively, for chloronium in these two sources. We obtained upper limits on the para-H2-35Cl+ line strengths toward H2 Peak 1 in the Orion Molecular cloud and toward the massive young star AFGL 2591. The chloronium abundances inferred in this study are typically at least a factor ~10 larger than the predictions of steady-state theoretical models for the chemistry of interstellar molecules containing chlorine. Several explanations for this discrepancy were investigated, but none has proven satisfactory, and thus the large observed abundances of chloronium remain puzzling.

Observations of H3+ in the Galactic diffuse interstellar medium (ISM) have led to various surprising results, including the conclusion that the cosmic-ray ionization rate (zeta_2) is about 1 order of magnitude larger than previously thought. The present survey expands the sample of diffuse cloud sight lines with H3+ observations to 50, with detections in 21 of those. Ionization rates inferred from these observations are in the range (1.7+-1.3)x10^-16 s^-1<zeta_2<(10.6+-8.2)x10^-16 s^-1 with a mean value of zeta_2=(3.5^+5.3_-3.0)x10^-16 s^-1. Upper limits (3 sigma) derived from non-detections of H3+ are as low as zeta_2<0.4x10^-16 s^-1. These low upper-limits, in combination with the wide range of inferred cosmic-ray ionization rates, indicate variations in zeta_2 between different diffuse cloud sight lines. A study of zeta_2 versus N_H (total hydrogen column density) shows that the two parameters are not correlated for diffuse molecular cloud sight lines, but that the ionization rate decreases when N_H increases to values typical of dense molecular clouds. Both the difference in ionization rates between diffuse and dense clouds and the variation of zeta_2 among diffuse cloud sight lines are likely the result of particle propagation effects. The lower ionization rate in dense clouds is due to the inability of low-energy (few MeV) protons to penetrate such regions, while the ionization rate in diffuse clouds is controlled by the proximity of the observed cloud to a site of particle acceleration.

We present follow-up observations to those of Geballe & Oka (2010), who found high column densities of H3+ ~100 pc off of the Galactic center (GC) on the lines of sight to 2MASS J17432173-2951430 (J1743) and 2MASS J17470898-2829561 (J1747). The wavelength coverages on these sightlines have been extended in order to observe two key transitions of H3+, R(3,3)l and R(2,2)l, that constrain the temperatures and densities of the environments. The profiles of the H3+ R(3,3)l line, which is due only to gas in the GC, closely matches the differences between the H3+ R(1,1)l and CO line profiles, just as it does for previously studied sightlines in the GC. Absorption in the R(2,2)l line of H3+ is present in J1747 at velocities between -60 and +100 km/s. This is the second clear detection of this line in the interstellar medium after GCIRS 3 in the Central Cluster. The temperature of the absorbing gas in this velocity range is 350 K, significantly warmer than in the diffuse clouds in other parts of the Central Molecular Zone. This indicates that the absorbing gas is local to Sgr B molecular cloud complex. The warm and diffuse gas revealed by Oka et al. (2005) apparently extends to ~100 pc, but there is a hint that its temperature is somewhat lower in the line of sight to J1743 than elsewhere in the GC. The observation of H3+ toward J1747 is compared with the recent Herschel observation of H2O+ toward Sgr B2 and their chemical relationship and remarkably similar velocity profiles are discussed.

The excitation temperature T_01 derived from the relative intensities of the J = 0 (para) and J = 1 (ortho) rotational levels of H2 has been assumed to be an accurate measure of the kinetic temperature in interstellar environments. In diffuse molecular clouds, the average value of T_01 is ~70 K. However, the excitation temperature T(H3+) derived from the (J,K) = (1,1) (para) and (1,0) (ortho) rotational levels of H3+ has been observed to be ~30 K in the same types of environments. In this work, we present observations of H3+ in three additional diffuse cloud sight lines for which H2 measurements are available, showing that in 4 of 5 cases T_01 and T(H3+) are discrepant. We then examine the thermalization mechanisms for the ortho:para ratios of H3+ and H2, concluding that indeed T_01 is an accurate measure of the cloud kinetic temperature, while the ortho:para ratio of H3+ need not be thermal. By constructing a steady-state chemical model taking into account the nuclear-spindependence of reactions involving H3+, we show that the ortho:para ratio of H3+ in diffuse molecular clouds is likely governed by a competition between dissociative recombination with electrons and thermalization via reactive collisions with H2.

Observational and theoretical evidence suggests that high-energy Galactic cosmic rays are primarily accelerated by supernova remnants. If also true for low-energy cosmic rays, the ionization rate near a supernova remnant should be higher than in the general Galactic interstellar medium (ISM). We have searched for H3+ absorption features in 6 sight lines which pass through molecular material near IC 443---a well-studied case of a supernova remnant interacting with its surrounding molecular material---for the purpose of inferring the cosmic-ray ionization rate in the region. In 2 of the sight lines (toward ALS 8828 and HD 254577) we find large H3+ column densities, N(H3+)~3*10^14 cm^-2, and deduce ionization rates of zeta_2~2*10^-15 s^-1, about 5 times larger than inferred toward average diffuse molecular cloud sight lines. However, the 3 sigma upper limits found for the other 4 sight lines are consistent with typical Galactic values. This wide range of ionization rates is likely the result of particle acceleration and propagation effects, which predict that the cosmic-ray spectrum and thus ionization rate should vary in and around the remnant. While we cannot determine if the H3+ absorption arises in post-shock (interior) or pre-shock (exterior) gas, the large inferred ionization rates suggest that IC 443 is in fact accelerating a large population of low-energy cosmic rays. Still, it is unclear whether this population can propagate far enough into the ISM to account for the ionization rate inferred in diffuse Galactic sight lines.

The formation of CH+ in the interstellar medium has long been an outstanding problem in chemical models. In order to probe the physical conditions of the ISM in which CH+ forms, we propose the use of CH3+ observations. The pathway to forming CH3+ begins with CH+, and a steady state analysis of CH3+ and the reaction intermediary CH2+ results in a relationship between the CH+ and CH3+ abundances. This relationship depends on the molecular hydrogen fraction, f_H2, and gas temperature, T, so observations of CH+ and CH3+ can be used to infer the properties of the gas in which both species reside. We present observations of both molecules along the diffuse cloud sight line toward Cyg OB2 No. 12. Using our computed column densities and upper limits, we put constraints on the f_H2 vs. T parameter space in which CH+ and CH3+ form. We find that average, static, diffuse molecular cloud conditions (i.e. f_H2>0.2, T~60 K) are excluded by our analysis. However, current theory suggests that non-equilibrium effects drive the reaction C+ + H_2 --> CH+ + H, endothermic by 4640 K. If we consider a higher effective temperature due to collisions between neutrals and accelerated ions, the CH3+ partition function predicts that the overall population will be spread out into several excited rotational levels. As a result, observations of more CH3+ transitions with higher signal-to-noise ratios are necessary to place any constraints on models where magnetic acceleration of ions drives the formation of CH+.

The ionization rate of interstellar material by cosmic rays has been a major source of controversy, with different estimates varying by three orders of magnitude. Observational constraints of this rate have all depended on analyzing the chemistry of various molecules that are produced following cosmic-ray ionization, and in many cases these analyses contain significant uncertainties. Even in the simplest case (H3+) the derived ionization rate depends on an (uncertain) estimate of the absorption path length. In this paper we examine the feasibility of inferring the cosmic-ray ionization rate using the 10830 A absorption line of metastable helium. Observations through the diffuse clouds toward HD 183143 are presented, but yield only an upper limit on the metastable helium column density. A thorough investigation of He+ chemistry reveals that only a small fraction of He+ will recombine into the triplet state and populate the metastable level. In addition, excitation to the triplet manifold of helium by secondary electrons must be accounted for as it is the dominant mechanism which produces He* in some environments. Incorporating these various formation and destruction pathways, we derive new equations for the steady state abundance of metastable helium. Using these equations in concert with our observations, we find zeta_He < 1.2*10^-15 s^-1, an upper limit about 5 times larger than the ionization rate previously inferred for this sight line using H3+. While observations of interstellar He* are extremely difficult at present...

Diffuse interstellar clouds show large abundances of H_3^+ which can be maintained only by a high ionization rate of H_2. Cosmic rays are the dominant ionization mechanism in this environment, so the large ionization rate implies a high cosmic-ray flux, and a large amount of energy residing in cosmic rays. In this paper we find that the standard propagated cosmic-ray spectrum predicts an ionization rate much lower than that inferred from H_3^+. Low-energy (~10 MeV) cosmic rays are the most efficient at ionizing hydrogen, but cannot be directly detected; consequently, an otherwise unobservable enhancement of the low-energy cosmic-ray flux offers a plausible explanation for the H_3^+ results. Beyond ionization, cosmic rays also interact with the interstellar medium by spalling atomic nuclei and exciting atomic nuclear states. These processes produce the light elements Li, Be, and B, as well as gamma-ray lines. To test the consequences of an enhanced low-energy cosmic-ray flux, we adopt two physically-motivated cosmic-ray spectra which by construction reproduce the ionization rate inferred in diffuse clouds, and investigate the implications of these spectra on dense cloud ionization rates, light element abundances, gamma-ray fluxes, and energetics. One spectrum proposed here provides an explanation for the high ionization rate seen in diffuse clouds while still appearing to be broadly consistent with other observables, but the shape of this spectrum suggests that supernovae remnants may not be the predominant accelerators of low-energy cosmic rays.

Infrared absorption lines of H3+, including the metastable R(3,3)l line, have been observed toward eight bright infrared sources associated with hot and massive stars located in and between the Galactic Center Cluster and the Quintuplet Cluster 30 pc to the east. The absorption lines with high velocity dispersion arise in the Galaxy's Central Molecular Zone (CMZ) as well as in foreground spiral arms. The temperature and density of the gas in the CMZ, as determined from the relative strengths of the H3+ lines, are T=200-300K and n=50-200cm^-3. The detection of high column densities of H3+ toward all eight stars implies that this warm and diffuse gaseous environment is widespread in the CMZ. The products of the ionization rate and path length for these sight lines are 1000 and 10 times higher than in dense and diffuse clouds in the Galactic disk, respectively, indicating that the ionization rate, zeta, is not less than 10^-15 s^-1 and that L is at least on the order of 50 pc. The warm and diffuse gas is an important component of the CMZ, in addition to the three previously known gaseous environments: (1) cold molecular clouds observed by radio emission of CO and other molecules, (2) hot (T=10^4-10^6K) and highly ionized diffuse gas (n_e=10-100cm^-3) seen in radio recombination lines, far infrared atomic lines, and radio-wave scattering, and (3) ultra-hot (T=10^7-10^8K) X-ray emitting plasma. Its prevalence significantly changes the understanding of the environment of the CMZ. The sight line toward GC IRS 3 is unique in showing an additional H3+ absorption component, which is interpreted as due to either a cloud associated with circumnuclear disk or the "50 km s^-1 cloud" known from radio observations. An infrared pumping scheme is examined as a mechanism to populate the (3,3) metastable level in this cloud.

Using high resolution infrared spectroscopy we have surveyed twenty sightlines for H3+ absorption. H3+ is detected in eight diffuse cloud sightlines with column densities varying from 0.6x10^14 cm^-2 to 3.9x10^14 cm^-2. This brings to fourteen the total number of diffuse cloud sightlines where H3+ has been detected. These detections are mostly along sightlines concentrated in the Galactic plane, but well dispersed in Galactic longitude. The results imply that abundant H3+ is common in the diffuse interstellar medium. Because of the simple chemistry associated with H3+ production and destruction, these column density measurements can be used in concert with various other data to infer the primary cosmic-ray ionization rate, zeta_p. Values range from 0.5x10^-16 s^-1 to 3x10^-16 s^-1 with an average of 2x10^-16 s^-1. Where H3+ is not detected the upper limits on the ionization rate are consistent with this range. The average value of zeta_p is about an order of magnitude larger than both the canonical rate and rates previously reported by other groups using measurements of OH and HD. The discrepancy is most likely due to inaccurate measurements of rate constants and the omission of effects which were unknown when those studies were performed. We believe that the observed column density of H3+ is the most direct tracer for the cosmic-ray ionization rate due to its simple chemistry. Recent models of diffuse cloud chemistry require cosmic-ray ionization rates on the order of 10^-16 s^-1 to reproduce observed abundances of various atomic and molecular species, in rough accord with our observational findings.