Impact of cosmic-ray propagation on the chemistry and ionisation fraction of dark clouds

Abstract

A proper modelling of the cosmic-ray ionisation rate within gas clouds is crucial to describe their chemical evolution accurately. However, this modelling is computationally demanding because it requires the propagation of cosmic rays throughout the cloud over time. We present a more efficient approach that simultaneously guarantees a reliable estimate of the cosmic-ray impact on the chemistry of prestellar cores. We introduce a numerical framework that mimics the cosmic-ray propagation within gas clouds and applies it to magnetohydrodynamic simulations performed with the code GIZMO. It simulates the cosmic-ray attenuation by computing the effective column density of H${}_2$ that is traversed, which is estimated using the same kernel weighting approach as employed in the simulation. The obtained cosmic-ray ionisation rate is then used in post-processing to study the chemical evolution of the clouds. We found that cosmic-ray propagation affects deuterated and non-deuterated species significantly and that it depends on the assumed cosmic-ray spectrum. We explored correlations between the electron abundance, the cosmic-ray ionisation rate, and the abundance of the most relevant ions (HCO${}^{+}$, N${}_2$H${}^{+}$, DCO${}^{+}$, N${}_2$D${}^{+}$, and o-H${}_2$D${}^{+}$), with the purpose of finding simple expressions that link them. We provide an analytical formula to estimate the ionisation fraction, X(e${}^{-}$), from observable tracers and applied it to existing observations of high-mass clumps. We obtained values of about 10${}^{-8}$, which is in line with previous works and with expectations for dense clouds. We also provide a linear fit to calculate the cosmic-ray ionisation rate from the local H${}_2$ density, which is to be employed in three-dimensional simulations that do not include cosmic-ray propagation.

arXiv comment10 pages, 9 figures, 1 table

References in corpus (38)

• The cosmological simulation code GADGET-2
• GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods
• Cold Dark Clouds: The Initial Conditions for Star Formation
• ATLASGAL - The APEX Telescope Large Area Survey of the Galaxy at 870 microns
• Investigating the Cosmic-Ray Ionization Rate in the Galactic Diffuse Interstellar Medium through Observations of H3+
• H3+ in Diffuse Interstellar Clouds: a Tracer for the Cosmic-Ray Ionization Rate
• KROME - a package to embed chemistry in astrophysical simulations
• Cosmic-ray ionisation in circumstellar discs
• Accurate, Meshless Methods for Magneto-Hydrodynamics
• Weeds: a CLASS extension for the analysis of millimeter and sub-millimeter spectral surveys
• Cosmic ray feedback in the FIRE simulations: constraining cosmic ray propagation with GeV gamma ray emission
• The cosmic ray ionization rate in the Galactic disk, as determined from observations of molecular ions
• Impact of low-energy cosmic rays on star formation
• ATLASGAL-selected massive clumps in the inner Galaxy: I. CO depletion and isotopic ratios
• Cosmic-ray ionisation in collapsing clouds
• ATLASGAL-selected massive clumps in the inner Galaxy: V. Temperature structure and evolution
• The Ionization Fraction in Dense Molecular Gas II: Massive Cores
• The Ionization Fraction of Barnard 68: Implications for Star and Planet Formation
• Cosmic rays in molecular clouds probed by H${}_2$ rovibrational lines -- Perspectives for the James Webb Space Telescope
• Cosmic Ray Transport, Energy Loss, and Influence in the Multiphase Interstellar Medium
• Evolution of cosmic ray electron spectra in magnetohydrodynamical simulations
• Chemical Analysis of a Diffuse Cloud along a Line of Sight Toward W51: Molecular Fraction and Cosmic-Ray Ionization Rate
• On the enhanced cosmic-ray ionization rate in the diffuse cloud towards Zeta Persei
• The cosmic-ray ionisation rate in the pre-stellar core L1544
• Survey of ortho-H${}_2$D${}^{+}$ in high-mass star-forming regions
• Magnetic Mirroring and Focusing of Cosmic Rays
• A new proxy to estimate the cosmic-ray ionisation rate in dense cores
• Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. I. Spectrally resolved propagation of CR electrons on Eulerian grids
• First ALMA maps of cosmic ray ionisation rate in high-mass star-forming regions
• A finite volume method for two-moment cosmic-ray hydrodynamics on a moving mesh
• The ALMA Survey of 70 $\mu \mathrm{m}$ Dark High-mass Clumps in Early Stages (ASHES). VII: Chemistry of Embedded Dense Cores
• The 3D structure of CO depletion in high-mass prestellar regions
• The impact of cosmic-ray attenuation on the carbon cycle emission in molecular clouds
• Constraining the cosmic-ray ionization rate and their spectrum with NIR spectroscopy of dense clouds -- A test-bed for JWST
• On the Inference of the Cosmic Ray Ionization Rate $\zeta$ from the HCO${}^{+}$-to-DCO${}^{+}$ Abundance Ratio: The Effect of Nuclear Spin
• Chemical post-processing of magneto-hydrodynamical simulations of star-forming regions: robustness and pitfalls
• Parsec-scale cosmic-ray ionisation rate in Orion
• APEX observations of ortho-H${}_2$D${}^{+}$ towards dense cores in the Orion B9 filament

Cited by in corpus (2)

• Cosmic-ray ionisation rate in low-mass cores: the role of the environment
• The impact of attenuation on cosmic-ray chemistry: I. Abundances and chemical calibrators in molecular clouds