The cosmic-ray ionisation rate in the pre-stellar core L1544
Journal
ASTRONOMY & ASTROPHYSICS
Date Issued
2021
Author(s)
DOI
10.1051/0004-6361/202141776
Abstract
Context. Cosmic rays (CRs), which are energetic particles mainly composed of protons and electrons, play an important role in the chemistry and dynamics of the interstellar medium. In dense environments, they represent the main ionising agent, hence driving the rich chemistry of molecular ions. Furthermore, they determine the ionisation fraction, which regulates the degree of coupling between the gas and the interstellar magnetic fields, and the heating of the gas. Estimates of the CR ionisation rate of molecular hydrogen (ζ2) span several orders of magnitude, depending on the targeted sources and on the method used.
Aims: Recent theoretical models have characterised the CR attenuation with increasing density. We aim to test these models for the attenuation of CRs in the low-mass pre-stellar core L1544.
Methods: We used a state-of-the-art gas-grain chemical model, which accepts the CR ionisation rate profile as input, to predict the abundance profiles of four ions: N2H+, N2D+, HC18O+, and DCO+. Non-local thermodynamic equilibrium radiative transfer simulations were run to produce synthetic spectra based on the derived abundances. These were compared with observations obtained with the Institut de Radioastronomie Millimétrique 30 m telescope.
Results: Our results indicate that a model with high ζ2 (>10−16 s−1) is excluded by the observations. Also the model with the standard ζ2 = 1.3 × 10−17 s−1 produces a worse agreement with respect to the attenuation model based on Voyager observations, which is characterised by an average ⟨ ζ2 ⟩ = 3 × 10−17 s−1 at the column densities typical of L1544. The single-dish data, however, are not sensitive to the attenuation of the CR profile, which changes only by a factor of two in the range of column densities spanned by the core model (N = 2−50 × 1021 cm−2). Interferometric observations at higher spatial resolution, combined with observations of transitions with lower critical density - hence tracing the low-density envelope - are needed to observe a decrease in the CR ionisation rate with density.
Aims: Recent theoretical models have characterised the CR attenuation with increasing density. We aim to test these models for the attenuation of CRs in the low-mass pre-stellar core L1544.
Methods: We used a state-of-the-art gas-grain chemical model, which accepts the CR ionisation rate profile as input, to predict the abundance profiles of four ions: N2H+, N2D+, HC18O+, and DCO+. Non-local thermodynamic equilibrium radiative transfer simulations were run to produce synthetic spectra based on the derived abundances. These were compared with observations obtained with the Institut de Radioastronomie Millimétrique 30 m telescope.
Results: Our results indicate that a model with high ζ2 (>10−16 s−1) is excluded by the observations. Also the model with the standard ζ2 = 1.3 × 10−17 s−1 produces a worse agreement with respect to the attenuation model based on Voyager observations, which is characterised by an average ⟨ ζ2 ⟩ = 3 × 10−17 s−1 at the column densities typical of L1544. The single-dish data, however, are not sensitive to the attenuation of the CR profile, which changes only by a factor of two in the range of column densities spanned by the core model (N = 2−50 × 1021 cm−2). Interferometric observations at higher spatial resolution, combined with observations of transitions with lower critical density - hence tracing the low-density envelope - are needed to observe a decrease in the CR ionisation rate with density.
This work is based on observations carried out with the IRAM 30 m telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain).
Volume
656
Start page
A109
File(s)![Thumbnail Image]()
![Thumbnail Image]()
Loading...
Name
aa41776-21.pdf
Description
PDF editoriale
Size
1.41 MB
Format
Adobe PDF
Checksum (MD5)
5f8f82de986ef2fa81d45148e58553a6
Loading...
Name
2109.08169.pdf
Description
postprint
Size
1.01 MB
Format
Adobe PDF
Checksum (MD5)
e4240137d086908e2e4c1d9b00336439