Cosmic-ray Ionization Rate Research Articles

Magnetic fields do not suppress global star formation in low metallicity dwarf galaxies

ABSTRACT Many studies concluded that magnetic fields suppress star formation in molecular clouds and Milky Way like galaxies. However, most of these studies are based on fully developed fields that have reached the saturation level, with little work on investigating how an initial weak primordial field affects star formation in low metallicity environments. In this paper, we investigate the impact of a weak initial field on low metallicity dwarf galaxies. We perform high-resolution arepo simulations of five isolated dwarf galaxies. Two models are hydrodynamical, two start with a primordial magnetic field of 10$^{-6} \, \mu$G and different sub-solar metallicities, and one starts with a saturated field of 10$^{-2} \, \mu$G. All models include a non-equilibrium, time-dependent chemical network that includes the effects of gas shielding from the ambient ultraviolet field. Sink particles form directly from the gravitational collapse of gas and are treated as star-forming clumps that can accrete gas. We vary the ambient uniform far ultraviolet field, and cosmic ray ionization rate between 1 per cent and 10 per cent of solar values. We find that the magnetic field has little impact on the global star formation rate (SFR), which is in tension with some previously published results. We further find that the initial field strength has little impact on the global SFR. We show that an increase in the mass fractions of both molecular hydrogen and cold gas, along with changes in the perpendicular gas velocity dispersion and the magnetic field acting in the weak-field model, overcome the expected suppression in star formation.

Dependence of Chemical Abundance on the Cosmic-Ray Ionization Rate in IC 348

Ions (e.g., , H2O+) have been used extensively to quantify the cosmic-ray ionization rate (CRIR) in diffuse sightlines. However, measurements of the CRIR in low-to-intermediate density gas environments are rare, especially when background stars are absent. In this work, we combine molecular line observations of CO, OH, CH, and HCO+ in the star-forming cloud IC 348, and chemical models to constrain the value of the CRIR and study the response of the distribution of chemical abundances. The cloud boundary is found to have an A V of approximately 4 mag. From the interior to the exterior of the cloud, the observed 13CO line intensities drop by an order of magnitude. The calculated average abundance of 12CO (assuming 12C/13C=65) is (1.2 ± 0.9) × 10−4, which decreases by a factor of 6 from the interior to the outside regions. The average abundance of CH (3.3 ± 0.7 × 10−8) is in good agreement with previous findings in diffuse and translucent clouds (A V < 5 mag). However, we did not find a decline in CH abundance in regions of high extinction (A V ≃ 8 mag) as previously reported in Taurus. By comparing the observed molecular abundances and chemical models, we find a decreasing trend of the CRIR as A V increases. The inferred CRIR of ζ cr = (4.7 ± 1.5) × 10−16 s−1 at low A V is consistent with measurements toward two nearby massive stars.

Magnetically regulated collapse in the B335 protostar?

Context. Whether or not magnetic fields play a key role in dynamically shaping the products of the star formation process is still largely debated. For example, in magnetized protostellar formation models, magnetic braking plays a major role in the regulation of the angular momentum transported from large envelope scales to the inner envelope, and is expected to be responsible for the resulting protostellar disk sizes. However, non-ideal magnetohydrodynamic effects that rule the coupling of the magnetic field to the gas also depend heavily on the local physical conditions, such as the ionization fraction of the gas. Aims. The purpose of this work is to observationally characterize the level of ionization of the gas at small envelope radii and to investigate its relation to the efficiency of the coupling between the star-forming gas and the magnetic field in the Class 0 protostar B335. Methods. We obtained molecular line emission maps of B335 with ALMA, which we use to measure the deuteration fraction of the gas, RD, its ionization fraction, χe, and the cosmic-ray ionization rate, ζCR, at envelope radii ≲1000 au. Results. We find large fractions of ionized gas, χe ≃ 1–8 × 10−6. Our observations also reveal an enhanced ionization that increases at small envelope radii, reaching values up to ζCR ≃ 10−14 s−1 at a few hundred astronomical units (au) from the central protostellar object. We show that this extreme ζCR can be attributed to the presence of cosmic rays accelerated close to the protostar. Conclusions. We report the first resolved map of ζCR at scales ≲1000 au in a solar-type Class 0 protostar, finding remarkably high values. Our observations suggest that local acceleration of cosmic rays, and not the penetration of interstellar Galactic cosmic rays, may be responsible for the gas ionization in the inner envelope, potentially down to disk-forming scales. If confirmed, our findings imply that protostellar disk properties may also be determined by local processes that set the coupling between the gas and the magnetic field, and not only by the amount of angular momentum available at large envelope scales and the magnetic field strength in protostellar cores. We stress that the gas ionization we find in B335 significantly stands out from the typical values routinely used in state-of-the-art models of protostellar formation and evolution. If the local processes of ionization uncovered in B335 are prototypical to low-mass protostars, our results call for a revision of the treatment of ionizing processes in magnetized models for star and disk formation.

PDFchem: A new fast method to determine ISM properties and infer environmental parameters using probability distributions

ABSTRACT Determining the atomic and molecular content of the interstellar medium (ISM) is of fundamental importance to understand the star-formation process. Although there exist various three-dimensional hydro-chemical codes modelling the ISM, they are computationally expensive and inefficient for studies over a large parameter space. Building on our earlier approach, we present PDFchem, a novel algorithm that models the cold ISM at moderate and large scales using functions connecting the quantities of the local (AV, eff) and the observed (AV, obs) visual extinctions, and the local number density, nH, with probability density functions (PDF) of AV, obs on cloud scales typically tens-to-hundreds of pc as an input. For any given AV, obs-PDF, the algorithm instantly computes the average abundances of the most important species (H i, H2, C ii, C i, CO, OH, OH+, H2O+, CH, HCO+) and performs radiative transfer calculations to estimate the average emission of the most commonly observed lines ([C ii] 158$\mu$m, both [C i] fine-structure lines and the first five rotational transitions of 12CO). We examine two AV, obs-PDFs corresponding to a non-star-forming and a star-forming ISM region, under a variety of environmental parameters combinations. These cover far-ultraviolet intensities in the range of χ/χ0 = 10−1 − 103, cosmic ray ionization rates in the range of $\zeta _{\rm CR}=10^{-17}-10^{-13}\, {\rm s}^{-1}$ and metallicities in the range of $Z=0.1-2\, {\rm Z}_{\odot }$. PDFchem is fast, easy to use, reproduces the PDR quantities of the hydrodynamical models, and can be used directly with observed data of the cold ISM.

Tracing Interstellar Heating: An ALCHEMI Measurement of the HCN Isomers in NGC 253

We analyze HCN and HNC emission in the nearby starburst galaxy NGC 253 to investigate its effectiveness in tracing heating processes associated with star formation. This study uses multiple HCN and HNC rotational transitions observed using the Atacama Large Millimeter/submillimeter Array via the ALCHEMI Large Program. To understand the conditions and associated heating mechanisms within NGC 253's dense gas, we employ Bayesian nested sampling techniques applied to chemical and radiative transfer models, which are constrained using our HCN and HNC measurements. We find that the volume density and cosmic-ray ionization rate (CRIR) ζ are enhanced by about an order of magnitude in the galaxy’s central regions as compared to those further from the nucleus. In NGC 253's central giant molecular clouds (GMCs), where observed HCN/HNC abundance ratios are the lowest, n ∼ 105.5 cm−3 and ζ ∼ 10−12 s−1 (greater than 104 times the average Galactic rate). We find a positive correlation in the association of both density and CRIR with the number of star formation-related heating sources (supernova remnants, H ii regions, and super hot cores) located in each GMC, as well as a correlation between CRIRs and supernova rates. Additionally, we see an anticorrelation between the HCN/HNC ratio and CRIR, indicating that this ratio will be lower in regions where ζ is higher. Though previous studies suggested HCN and HNC may reveal strong mechanical heating processes in NGC 253's CMZ, we find cosmic-ray heating dominates the heating budget, and mechanical heating does not play a significant role in the HCN and HNC chemistry.

Cosmic-Ray Ionization Rate in Protoplanetary Disks with Sheared Magnetic Fields

We investigate the effects of magnetic-field configurations on the ionization rate by cosmic rays in protoplanetary disks. First, we consider cosmic-ray propagation from the interstellar medium (ISM) to the protoplanetary disks and showed that the cosmic-ray density around the disk should be 2 times lower than the ISM value. Then, we compute the attenuation of cosmic rays in protoplanetary disks. The magnetic fields in the disk are stretched to the azimuthal directions, and cosmic rays need to detour while propagating to the midplane. Our results show that the detouring effectively enhances the column density by about two orders of magnitude. We employ a typical ionization rate by cosmic rays in diffuse ISM, which is considered too high to be consistent with observations of protoplanetary disks, and find that the cosmic rays are significantly shielded at the midplane. In the case of the disk around IM Lup, the midplane ionization rate is very low for r ≲ 100 au, while the value is as large as a diffuse ISM in the outer radii. Our results are consistent with the recent Atacama Large Millimeter/submillimeter Array observation that indicates the radial gradient in the cosmic-ray ionization rate. The high ionization rate in the outer radii of disks may activate the magnetorotational instability that was thought to be suppressed due to ambipolar diffusion. These results will have a strong influence on the dynamical and chemical evolutions of protoplanetary disks.

Cosmic-ray-induced H2 line emission

Context. It has been proposed that H2 near-infrared lines may be excited by cosmic rays and thus allow for a determination of the cosmic-ray ionization rate in dense gas. One-dimensional models show that measuring both the H2 gas column density and H2 line intensity enables a constraint on the cosmic-ray ionization rate as well as on the spectral slope of low-energy cosmic-ray protons in the interstellar medium. Aims. We aim to investigate the impact of certain assumptions regarding the H2 chemical models and interstellar medium density distributions on the emission of cosmic-ray-induced H2 emission lines. This is of particular importance for utilizing observations of these lines with the James Webb Space Telescope to constrain the cosmic-ray ionization rate. Methods. We compare the predicted emission from cosmic-ray-induced, rovibrationally excited H2 emission lines for different one-and three-dimensional models with varying assumptions regarding the gas chemistry and density distribution. Results. We find that the model predictions of the H2 line intensities for the (1-0)S(0), (1-0)Q(2), (1-0)O(2), and (1-0)O(4) transitions at 2.22, 2.41, 2.63, and 3.00 μm, respectively, are relatively independent of the astro-chemical model and the gas density distribution when compared against the H2 column density, making them robust tracers of the cosmic-ray ionization rate. Conclusions. We recommend the use of rovibrational H2 line emission in combination with estimations of the cloud’s H2 column density to constrain the ionization rate and the spectrum of low-energy cosmic rays.

First predicted cosmic ray spectra, primary-to-secondary ratios, and ionization rates from MHD galaxy formation simulations

ABSTRACT We present the first simulations evolving resolved spectra of cosmic rays (CRs) from MeV–TeV energies (including electrons, positrons, (anti)protons, and heavier nuclei), in live kinetic-magnetohydrodynamics galaxy simulations with star formation and feedback. We utilize new numerical methods including terms often neglected in historical models, comparing Milky Way analogues with phenomenological scattering coefficients ν to Solar-neighbourhood [Local interstellar medium (LISM)] observations (spectra, B/C, e+/e−, $\mathrm{\bar{p}}/\mathrm{p}$, 10Be/9Be, ionization, and γ-rays). We show it is possible to reproduce observations with simple single-power-law injection and scattering coefficients (scaling with rigidity R), similar to previous (non-dynamical) calculations. We also find: (1) The circumgalactic medium in realistic galaxies necessarily imposes an $\sim 10\,$ kpc CR scattering halo, influencing the required ν(R). (2) Increasing the normalization of ν(R) re-normalizes CR secondary spectra but also changes primary spectral slopes, owing to source distribution and loss effects. (3) Diffusive/turbulent reacceleration is unimportant and generally sub-dominant to gyroresonant/streaming losses, which are sub-dominant to adiabatic/convective terms dominated by $\sim 0.1-1\,$ kpc turbulent/fountain motions. (4) CR spectra vary considerably across galaxies; certain features can arise from local structure rather than transport physics. (5) Systematic variation in CR ionization rates between LISM and molecular clouds (or Galactic position) arises naturally without invoking alternative sources. (6) Abundances of CNO nuclei require most CR acceleration occurs around when reverse shocks form in SNe, not in OB wind bubbles or later Sedov–Taylor stages of SNe remnants.

Astrochemical modelling of infrared dark clouds

Context. Infrared dark clouds (IRDCs) are cold, dense regions of the interstellar medium (ISM) that are likely to represent the initial conditions for massive star and star cluster formation. It is thus important to study the physical and chemical conditions of IRDCs to provide constraints and inputs for theoretical models of these processes. Aims. We aim to determine the astrochemical conditions, especially the cosmic ray ionisation rate (CRIR) and chemical age, in different regions of the massive IRDC G28.37+00.07 by comparing observed abundances of multiple molecules and molecular ions with the predictions of astrochemical models. Methods. We have computed a series of single-zone, time-dependent, astrochemical models with a gas-grain network that systematically explores the parameter space of the density, temperature, CRIR, and visual extinction. We have also investigated the effects of choices of CO ice binding energy and temperatures achieved in the transient heating of grains when struck by cosmic rays. We selected ten positions across the IRDC that are known to have a variety of star formation activity. We utilised mid-infrared extinction maps and sub-millimetre (sub-mm) emission maps to measure the mass surface densities of these regions needed for abundance and volume density estimates. The sub-mm emission maps were also used to measure temperatures. We then used Instituto de Radioas-tromía Milimétrica (IRAM) 30 m observations of various tracers, especially C18O(1-0), H13CO+(1-0), HC18O+(1-0), and N2H+(1-0), to estimate column densities and thus abundances. Finally, we investigated the range of astrochemical conditions that are consistent with the observed abundances. Results. The typical physical conditions of the IRDC regions are nH ~ 3 × 104 to 105 cm−3 and T ≃ 10 to 15 K. Strong emission of H13CO+(1-0) and N2H+(1-0) is detected towards all the positions and these species are used to define relatively narrow velocity ranges of the IRDC regions, which are used for estimates of CO abundances, via C18O(1-0). We would like to note that CO depletion factors are estimated to be in the range fD ~ 3 to 10. Using estimates of the abundances of CO, HCO+, and N2H+, we find consistency with astrochemical models that have relatively low CRIRs of ζ ~ 10−18 to ~10−17 s−1, with no evidence for systematic variation with the level of star formation activity. Astrochemical ages, which are defined with a reference to an initial condition of all H in H2, all C in CO, and all other species in atomic form, are found to be &lt;1 Myr. We also explore the effects of using other detected species, that is HCN, HNC, HNCO, CH3OH, and H2CO, to constrain the models. These generally lead to implied conditions with higher levels of CRIRs and older chemical ages. Considering the observed fD versus nH relation of the ten positions, which we find to have relatively little scatter, we discuss potential ways in which the astrochemical models can match such a relation as a quasi-equilibrium limit valid at ages of at least a few free-fall times, that is ≳0.3 Myr, including the effect of CO envelope contamination, small variations in temperature history near 15 K, CO-ice binding energy uncertainties, and CR-induced desorption. We find general consistency with the data of ~0.5 Myr-old models that have ζ ~ 2-5 × 10−18 s−1 and CO abundances set by a balance of freeze-out with CR-induced desorption. Conclusions. We have constrained the astrochemical conditions in ten regions in a massive IRDC, finding evidence for relatively low values of CRIR compared to diffuse ISM levels. We have not seen clear evidence for variation in the CRIR with the level of star formation activity. We favour models that involve relatively low CRIRs (≲10−17 s−1) and relatively old chemical ages (≳0.3 Myr, i.e. ≳3tff). We discuss potential sources of systematic uncertainties in these results and the overall implications for IRDC evolutionary history and astrochemical models.

The Young Embedded Disk L1527 IRS: Constraints on the Water Snowline and Cosmic-Ray Ionization Rate from HCO+ Observations

The water snowline in circumstellar disks is a crucial component in planet formation, but direct observational constraints on its location remain sparse owing to the difficulty of observing water in both young embedded and mature protoplanetary disks. Chemical imaging provides an alternative route to locate the snowline, and HCO+ isotopologues have been shown to be good tracers in protostellar envelopes and Herbig disks. Here we present ∼0.″5 resolution (∼35 au radius) Atacama Large Millimeter/submillimeter Array (ALMA) observations of HCO+ J = 4 − 3 and H13CO+ J = 3 − 2 toward the young (Class 0/I) disk L1527 IRS. Using a source-specific physical model with the midplane snowline at 3.4 au and a small chemical network, we are able to reproduce the HCO+ and H13CO+ emission, but for HCO+ only when the cosmic-ray ionization rate is lowered to 10−18 s−1. Even though the observations are not sensitive to the expected HCO+ abundance drop across the snowline, the reduction in HCO+ above the snow surface and the global temperature structure allow us to constrain a snowline location between 1.8 and 4.1 au. Deep observations are required to eliminate the envelope contribution to the emission and to derive more stringent constraints on the snowline location. Locating the snowline in young disks directly with observations of H2O isotopologues may therefore still be an alternative option. With a direct snowline measurement, HCO+ will be able to provide constraints on the ionization rate.

Chemical Evolution During the Formation of Molecular Clouds

To study the chemical evolution during the formation of molecular clouds, we model three types of clouds with different density structures: collapsing spherical, collapsing ellipsoidal, and static spherical profiles. The collapsing models are better than the static models in matching the observational characteristics in typical molecular clouds. This is mainly because the gravity can speed up the formation of some important molecules (e.g., H2, CO, OH) by increasing the number density during collapse. The different morphologies of prolate, oblate, and spherical clouds lead to differences in chemical evolution, which are mainly due to their different evolution of number density. We also study the effect of initial chemical compositions on chemical evolution, and find that H atoms can accelerate OH formation by two major reactions: O + H → OH in gas phase and on dust grain surfaces, leading to the models in which hydrogen is mainly atomic initially better match observations than the models in which hydrogen is mainly molecular initially. Namely, to match observations, initially hydrogen must be mostly atomic. The CO molecules are able to form even without the pre-existence of H2. We also study the influence of gas temperature, dust temperature, intensity of interstellar radiation field and cosmic-ray ionization rate on chemical evolution in static clouds. The static CO clouds with high dust temperature, strong radiation field, and intensive cosmic rays are transient due to rapid CO destruction.

Gas phase Elemental abundances in Molecular cloudS (GEMS)

Context. In the context of the IRAM 30 m Large Program Gas phase Elemental abundances in Molecular CloudS (GEMS), we present a study of thioformaldehyde (H2CS) and its deuterated versions (HDCS and D2CS) in several starless cores located in a selected set of star-forming filaments of Taurus, Perseus, and Orion. These regions have different star formation activities and, therefore, distinct physical and chemical conditions. Aims. Our goal is to investigate the influence of the environmental conditions on the abundances of these molecules in the cores, as well as the effect of time evolution. Methods. We have modelled the observed lines of H2CS, HDCS, and D2CS using the radiative transfer code RADEX. We have also used the chemical code Nautilus to model the evolution of these species depending on the characteristics of the starless cores. Results. We derive column densities and abundances for all the cores. We also derive deuterium fractionation ratios, Dfrac, which allow us to determine and compare the evolutionary stage between different parts of each star-forming region. Our results indicate that the north region of the B 213 filament in Taurus is more evolved than the south, while the north-eastern part of Perseus presents an earlier evolutionary stage than the south-western zone. Model results also show that Dfrac decreases with the cosmic-ray ionisation rate, while it increases with density and with the degree of sulphur depletion. In particular, we can only reproduce the observations when the initial sulphur depletion in the starless cores is at least one order of magnitude lower than the solar elemental sulphur abundance. Conclusions. The progressive increase in HDCS/H2CS and D2CS/H2CS with time makes these ratios powerful tools for deriving the chemical evolutionary stage of starless cores. However, they cannot be used to derive the temperature of these regions, since both ratios present a similar evolution at two different temperature ranges (~7–11 K and ~ 15–19 K). Regarding chemistry, (deuterated) thioformaldehyde is mainly formed through gas-phase reactions (double-replacement and neutral-neutral displacement reactions), while surface chemistry plays an important role as a destruction mechanism.

Unusually High HCO+/CO Ratios in and outside Supernova Remnant W49B

Galactic supernova remnants (SNRs) and their environments provide the nearest laboratories to study SN feedback. We performed molecular observations toward SNR W49B, the most luminous Galactic SNR in the X-ray band, aiming to explore signs of multiple feedback channels of SNRs on nearby molecular clouds (MCs). We found very broad HCO+ lines with widths of dv ∼ 48–75 km s−1 in the SNR southwest, providing strong evidence that W49B is perturbing MCs at a systemic velocity of V LSR = 61–65 km s−1, and placing the W49B at a distance of 7.9 ± 0.6 kpc. We observed unusually high-intensity ratios of HCO+ J=1–0/CO J=1–0 not only at shocked regions (1.1 ± 0.4 and 0.70 ± 0.16) but also in quiescent clouds over 1 pc away from the SNR’s eastern boundary (≥0.2). By comparing with the magnetohydrodynamics shock models, we interpret that the high ratio in the broad-line regions can result from a cosmic-ray (CR) induced chemistry in shocked MCs, where the CR ionization rate is enhanced to around 10–102 times of the Galactic level. The high HCO+/CO ratio outside the SNR is probably caused by the radiation precursor, while the luminous X-ray emission of W49B can explain a few properties in this region. The above results provide observational evidence that SNRs can strongly influence the molecular chemistry in and outside the shock boundary via their shocks, CRs, and radiation. We propose that the HCO+/CO ratio is a potentially useful tool to probe an SNR’s multichannel influence on MCs.

Energizing Star Formation: The Cosmic-Ray Ionization Rate in NGC 253 Derived from ALCHEMI Measurements of H3O+ and SO

The cosmic-ray ionization rate (CRIR) is a key parameter in understanding the physical and chemical processes in the interstellar medium. Cosmic rays are a significant source of energy in star formation regions, impacting the physical and chemical processes that drive the formation of stars. Previous studies of the circum-molecular zone of the starburst galaxy NGC 253 have found evidence for a high CRIR value: 103–106 times the average CRIR within the Milky Way. This is a broad constraint, and one goal of this study is to determine this value with much higher precision. We exploit ALMA observations toward the central molecular zone of NGC 253 to measure the CRIR. We first demonstrate that the abundance ratio of H3O+ and SO is strongly sensitive to the CRIR. We then combine chemical and radiative transfer models with nested sampling to infer the gas properties and CRIR of several star-forming regions in NGC 253 from emission from their transitions. We find that each of the four regions modeled has a CRIR in the range (1–80) × 10−14 s−1 and that this result adequately fits the abundances of other species that are believed to be sensitive to cosmic rays, including C2H, HCO+, HOC+, and CO. From shock and photon-dominated/X-ray dominated region models, we further find that neither UV-/X-ray-driven nor shock-dominated chemistry is a viable single alternative as none of these processes can adequately fit the abundances of all of these species.

Ionize Hard: Interstellar PO+ Detection

We report the first detection of the phosphorus monoxide ion (PO+) in the interstellar medium. Our unbiased and very sensitive spectral survey toward the G+0.693–0.027 molecular cloud covers four different rotational transitions of this molecule, two of which (J = 1–0 and J = 2–1) appear free of contamination from other species. The fit performed, assuming local thermodynamic equilibrium conditions, yields a column density of N=(6.0 ± 0.7) × 1011 cm−2. The resulting molecular abundance with respect to molecular hydrogen is 4.5 × 10–12. The column density of PO+ normalized by the cosmic abundance of P is larger than those of NO+ and SO+, normalized by N and S, by factors of 3.6 and 2.3, respectively. The N(PO+)/N(PO) ratio is 0.12 ± 0.03, more than one order of magnitude higher than that of N(SO+)/N(SO) and N(NO+)/N(NO). These results indicate that P is more efficiently ionized than N and S in the ISM. We have performed new chemical models that confirm that the PO+ abundance is strongly enhanced in shocked regions with high values of cosmic-ray ionization rates (10–15 − 10–14 s−1), as occurring in the G+0.693–0.027 molecular cloud. The shocks sputter the interstellar icy grain mantles, releasing into the gas phase most of their P content, mainly in the form of PH3, which is converted into atomic P, and then ionized efficiently by cosmic rays, forming P+. Further reactions with O2 and OH produces PO+. The cosmic-ray ionization of PO might also contribute significantly, which would explain the high N(PO+)/N(PO) ratio observed. The relatively high gas-phase abundance of PO+ with respect to other P-bearing species stresses the relevance of this species in the interstellar chemistry of P.

ArH+ and H2O+ absorption towards luminous galaxies

Context. Along several sight lines within the Milky Way ArH+ has been ubiquitously detected with only one detection in extragalactic environments, namely along two sight lines in the redshift z = 0.89 absorber towards the lensed blazar PKS 1830-211. Being formed in predominantly atomic gas by reactions between Ar+, which were initially ionised by cosmic rays and molecular hydrogen, ArH+ has been shown to be an excellent tracer of atomic gas as well as the impinging cosmic-ray ionisation rates. Aims. In this work, we attempt to extend the observations of ArH+ in extragalactic sources to examine its use as a tracer of the atomic interstellar medium (ISM) in these galaxies. Methods. We report the detection of ArH+ towards two luminous nearby galaxies, NGC 253 and NGC 4945, and the non-detection towards Arp 220 observed using the SEPIA660 receiver on the APEX 12 m telescope. In addition, the two sidebands of this receiver allowed us to observe the NKaKc = 11,0 − 10,1 transitions of another atomic gas tracer p-H2O+ at 607.227 GHz with the ArH+ line, simultaneously. We modelled the optically thin spectra of both species and compared their observed line profiles with that of other well-known atomic gas tracers such as OH+ and o-H2O+ and diffuse and dense molecular gas tracers HF and CO, respectively. Results. Assuming that the observed absorption from the ArH+, OH+, and H2O+ molecules are affected by the same flux of cosmic rays, we investigate the properties of the different cloud layers. Based on a steady-state analysis of the chemistry of these three species and using statistical equilibrium calculations, we estimate the molecular fraction traced by ArH+ to be ∼​10−3 and find that ArH+ resides in gas volumes with low electron densities. We further study the ortho-to-para ratio of H2O+ and find that the derived ratios do not significantly deviate from the equilibrium value of three with spin temperatures greater than 15 and 24 K.

The impact of cosmic-ray attenuation on the carbon cycle emission in molecular clouds

Context. Observations of the emission of the carbon cycle species (C, C+, CO) are commonly used to diagnose gas properties in the interstellar medium, but they are significantly sensitive to the cosmic-ray ionization rate. The carbon-cycle chemistry is known to be quite sensitive to the cosmic-ray ionization rate, ζ, controlled by the flux of low-energy cosmic rays which get attenuated through molecular clouds. However, astrochemical models commonly assume a constant cosmic-ray ionization rate in the clouds. Aims. We investigate the effect of cosmic-ray attenuation on the emission of carbon cycle species from molecular clouds, in particular the [CII] 158 μm, [CI] 609 μm, and CO (J = 1–0) 115.27 GHz lines. Methods. We used a post-processed chemical model of diffuse and dense simulated molecular clouds and quantified the variation in both column densities and velocity-integrated line emission of the carbon cycle with different cosmic-ray ionization rate models. Results. We find that the abundances and column densities of carbon cycle species are significantly impacted by the chosen cosmic-ray ionization rate model: no single constant ionization rate can reproduce the abundances modeled with an attenuated cosmic-ray model. Further, we show that constant ionization rate models fail to simultaneously reproduce the integrated emission of the lines we consider, and their deviations from a physically derived cosmic-ray attenuation model is too complex to be simply corrected. We demonstrate that the two clouds we modeled exhibit a similar average AV,eff – nH relationship, resulting in an average relation between the cosmic-ray ionization rate and density ζ(nH). Conclusions. We conclude by providing a number of implementation recommendations for cosmic rays in astrochemical models, but emphasize the necessity for column-dependent cosmic-ray ionization rate prescriptions.

Cosmic rays in molecular clouds probed by H2 rovibrational lines

Context. Low-energy cosmic rays (&lt;1 TeV) play a fundamental role in the chemical and dynamical evolution of molecular clouds, as they control the ionisation, dissociation, and excitation of H2. Their characterisation is therefore important both for the interpretation of observations and for the development of theoretical models. However, the methods used so far for estimating the cosmic-ray ionisation rate in molecular clouds have several limitations due to uncertainties in the adopted chemical networks. Aims. We refine and extend a previously proposed method to estimate the cosmic-ray ionisation rate in molecular clouds by observing rovibrational transitions of H2 at near-infrared wavelengths, which are mainly excited by secondary cosmic-ray electrons. Methods. Combining models of interstellar cosmic-ray propagation and attenuation in molecular clouds with the rigorous calculation of the expected secondary electron spectrum and updated electron-H2 excitation cross sections, we derive the intensity of the four H2 rovibrational transitions observable in cold dense gas: (1−0)O(2), (1−0)Q(2), (1−0)S(0), and (1−0)O(4). Results. The proposed method allows the estimation of the cosmic-ray ionisation rate for a given observed line intensity and H2 column density. We are also able to deduce the shape of the low-energy cosmic-ray proton spectrum impinging upon the molecular cloud. In addition, we present a look-up plot and a web-based application that can be used to constrain the low-energy spectral slope of the interstellar cosmic-ray proton spectrum. We finally comment on the capability of the James Webb Space Telescope to detect these near-infrared H2 lines, which will make it possible to derive, for the first time, spatial variation in the cosmic-ray ionisation rate in dense gas. Besides the implications for the interpretation of the chemical-dynamic evolution of a molecular cloud, it will finally be possible to test competing models of cosmic-ray propagation and attenuation in the interstellar medium, as well as compare cosmic-ray spectra in different Galactic regions.

The Role of Neutral Hydrogen in Setting the Abundances of Molecular Species in the Milky Way’s Diffuse Interstellar Medium. II. Comparison between Observations and Theoretical Models

Abstract We compare observations of H i from the Very Large Array (VLA) and the Arecibo Observatory and observations of HCO+ from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA) in the diffuse (A V ≲ 1) interstellar medium (ISM) to predictions from a photodissociation region (PDR) chemical model and multiphase ISM simulations. Using a coarse grid of PDR models, we estimate the density, FUV radiation field, and cosmic-ray ionization rate (CRIR) for each structure identified in HCO+ and H i absorption. These structures fall into two categories. Structures with T s &lt; 40 K, mostly with N(HCO+) ≲ 1012 cm−2, are consistent with modest density, FUV radiation field, and CRIR models, typical of the diffuse molecular ISM. Structures with spin temperature T s &gt; 40 K, mostly with N(HCO+) ≳ 1012 cm−2, are consistent with high density, FUV radiation field, and CRIR models, characteristic of environments close to massive star formation. The latter are also found in directions with a significant fraction of thermally unstable H i. In at least one case, we rule out the PDR model parameters, suggesting that alternative mechanisms (e.g., nonequilibrium processes like turbulent dissipation and/or shocks) are required to explain the observed HCO+ in this direction. Similarly, while our observations and simulations of the turbulent, multiphase ISM agree that HCO+ formation occurs along sight lines with N(H I) ≳ 1021 cm−2, the simulated data fail to explain HCO+ column densities ≳ few × 1012 cm−2. Because a majority of our sight lines with HCO+ had such high column densities, this likely indicates that nonequilibrium chemistry is important for these lines of sight.

Constraining the cosmic-ray ionization rate and spectrum with NIR spectroscopy of dense clouds

Context.Low-energy cosmic rays (CRs) control the thermo-chemical state and the coupling between gas and magnetic fields in dense molecular clouds. However, current estimates of the low-energy CR spectrum (E ≲ 1 GeV) and the associated CR ionization rate are highly uncertain.Aims.We apply, for the first time, a new method for constraining the CR ionization rate and the CR spectral shape using H2rovibrational lines from cold molecular clouds.Methods.Using the MMIRS instrument on the MMT, we obtained deep near-infrared (NIR) spectra in six positions within four dense cores, namely, G150, G157, G163, G198, with column densities ofNH2 ≈ 1022cm−2.Results.We derived 3σupper limits on the H2(1 − 0)S(0) line (2.22 μm) brightness in the rangeI = 5.9 × 10−8to 1.2 × 10−7erg cm−2s−1sr−1for the different targets. Using both an analytic model and a numerical model of CR propagation, we convert these into upper limits on the CR ionization rate in the clouds’ interior,ζ = 1.5 to 3.6 × 10−16s−1, and lower limits on the low-energy spectral slope of interstellar CR protons,α = −0.97 to −0.79. We show that while MMT was unable to detect the H2lines due to high atmospheric noise, JWST/NIRSpec will be able to efficiently detect the CR-excited H2lines, making it the ideal method for constraining the otherwise elusive low-energy CRs and shedding light on the sources and propagation modes of CRs.