Barger, C. J., Lam, K. H., Li, Z.-Y., Davis, S. W., Herbst, E. & Garrod, R. T.
Context. Gas-grain models have long been employed to simulate hot-core chemistry; however, these simulations have traditionally
neglected to couple chemical evolution in tandem with a rigorous physical evolution of a source. This over-simplification particularly
lacks an accurate treatment of temperature and spatial distribution, which are needed for realistic simulations of hot cores.
Aims. We aim to combine radiation hydrodynamics (RHD) with hot-core chemical kinetics in one dimension to produce a set of
astrochemical models that evolve according to explicitly calculated temperature, density, and spatial profiles.
Methods. We solve radiation hydrodynamics for three mass-accretion-rate models using Athena++. We then simulate the chemistry
using the hot-core chemical kinetic code MAGICKAL according to the physics derived from the RHD treatment.
Results. We find that as the mass-accretion rate decreases, the overall gas density of the source decreases. In particular, the gas density
for the lowest mass-accretion rate is low enough to restrict the proper formation of many complex organic molecules.We also compare
our chemical results in the form of calculated column densities to those of observations toward Sgr B2(N2). We find a generally good
agreement for oxygen-bearing species, particularly for the two highest mass-accretion rates.
Conclusions. Although we introduce hot-core chemical modeling using a self-consistent physical treatment, the adoption of a twodimensional
model may better reproduce chemistry and physics toward real sources and thus achieve better chemical comparisons with
Astronomy & Astrophysics, 2021, 651, A43
Full-text URL: https://ui.adsabs.harvard.edu/abs/2021A%26A...651A..43B/abstract