The evolution of HCO${}^{+}$ in molecular clouds using a novel chemical post-processing algorithm

Abstract

Modeling the chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational cost, 3D simulations generally restrict their chemistry to species with strong heating and cooling effects. Time-dependent information about the evolution of other species is therefore often neglected. We address this gap by post-processing tracer particles in the SILCC-Zoom molecular cloud simulations. Using a chemical network of 39 species and 301 reactions (including freeze-out of CO and H${}_2$O), and a novel algorithm to reconstruct a density grid from sparse tracer particle data, we produce time-dependent density distributions for various species. We focus upon the evolution of HCO${}^{+}$, which is a critical formation reactant of CO but is not typically modeled on-the-fly. We find that $\sim 90$% of the HCO${}^{+}$ content of the cold molecular gas forms in situ around $n_\textrm{HCO$^+$}\simeq10^3$-$10^4$ cm${}^{-3}$, over a time-scale of approximately 1 Myr. The remaining $\sim 10$% forms at high extinction sites, with minimal turbulent mixing out into the less dense gas. We further show that the dominant HCO${}^{+}$ formation pathway is dependent on the visual extinction, with the reaction H${}_3^{+}$ + CO contributing 90% of the total HCO${}^{+}$ production above $A_{\text{V,3D}}=3$. We produce the very first maps of the HCO${}^{+}$ column density, $N$(HCO${}^{+}$), and show that it reaches values as high as $10^{15}$ cm${}^{-2}$. We find that 50% of the HCO${}^{+}$ mass is located within $A_{\text{V}}\sim 10$-30, in a density range of $10^{3.5}$-$10^{4.5}$ cm${}^{-3}$. Our maps of $N$(HCO${}^{+}$) are shown to be in good agreement with recent observations of the W49A star-forming region.

25 pages including appendix, 22 figures, accepted for publication in MNRAS

References in corpus (26)

• Formation of methyl formate and other organic species in the warm-up phase of hot molecular cores
• Simulating the formation of molecular clouds. I. Slow formation by gravitational collapse from static initial conditions
• Simulating the formation of molecular clouds. II. Rapid formation from turbulent initial conditions
• CO-dark gas and molecular filaments in Milky Way type galaxies
• Molecular Line Emission as a Tool for Galaxy Observations (LEGO). I. HCN as a tracer of moderate gas densities in molecular clouds and galaxies
• ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- I. Survey description and a first look at G9.62+0.19
• HCN versus HCO+ as dense molecular gas mass tracer in Luminous Infrared Galaxies
• The environmental dependence of the X_CO conversion factor
• Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths
• Chemistry in Disks. IX. Observations and modeling of HCO+ and DCO+ in DM Tau
• Metallicity dependence of the H/H${}_2$ and C${}^{+}$/C/CO distributions in a resolved self-regulating interstellar medium
• The natural emergence of the correlation between H2 and star formation rate surface densities in galaxy simulations
• Photodissociation Region Diagnostics Across Galactic Environments
• LEGO II: A 3 mm molecular line study covering 100 pc of one of the most actively star-forming portions within the Milky Way Disc
• The effect of non-equilibrium metal cooling on the interstellar medium
• ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- II. Compact objects in ACA observations and star formation scaling relations
• On the accuracy of HI observations in molecular clouds -- More cold HI than thought?
• Modelling the chemistry of star-forming filaments - II. Testing filament characteristics with synthetic observations
• HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules -- I. Survey description and a first look toward W3(OH), W3 IRS5 and NGC 7538 IRS1
• H${}_2$ ortho-to-para conversion on grains: A route to fast deuterium fractionation in dense cloud cores?
• In search of infall motion in molecular clumps III: HCO+ (1-0) and H13CO+ (1-0) mapping observations toward the confirmed infall sources
• Chemical post-processing of magneto-hydrodynamical simulations of star-forming regions: robustness and pitfalls
• Synthetic CO emission and the $X_{\mathrm{CO}}$ factor of young molecular clouds: a convergence study
• TIMES I: a Systematic Observation in Multiple Molecular Lines Toward the Orion A and Ophiuchus Clouds
• CO+ as a probe of the origin of CO in diffuse interstellar clouds
• ALMA Observations of HCO+ and HCN Emission in a massive star forming region N55 of the Large Magellanic Cloud

Cited by in corpus (6)

• Non-Equilibrium Abundances Treated Holistically (NEATH): the molecular composition of star-forming clouds
• Mapping Synthetic Observations to Prestellar Core Models: An Interpretable Machine Learning Approach
• Observational studies on S-bearing molecules in massive star forming regions
• Shocks in the warm neutral medium II -- Origin of neutral carbon at high pressure
• A Search for 3-mm Molecular Absorption Line Transitions in the Magellanic Stream
• The impact of attenuation on cosmic-ray chemistry: I. Abundances and chemical calibrators in molecular clouds