The Astrophysical Journal (Jan 2023)
Modeling CO Line Profiles in Shocks of W28 and IC 443
Abstract
Molecular emission arising from the interactions of supernova remnant (SNR) shock waves and molecular clouds provide a tool for studying the dispersion and compression that might kick-start star formation as well as understanding cosmic-ray production. Purely rotational CO emission created by magnetohydrodynamic shock in the SNR–molecular cloud interaction is an effective shock tracer, particularly for slow-moving, continuous shocks into cold inner clumps of the molecular cloud. In this work, we present a new theoretical radiative transfer framework for predicting the line profile of CO with the Paris–Durham 1D shock model. We generated line profile predictions for CO emission produced by slow, magnetized C shocks into gas of density ∼10 ^4 cm ^−3 with shock speeds of 35 and 50 km s ^−1 . The numerical framework to reproduce the CO line profile utilizes the large velocity gradient (LVG) approximation and the omission of optically thick plane-parallel slabs. With this framework, we generated predictions for various CO spectroscopic observations up to J = 16 in SNRs W28 and IC 443, obtained with SOFIA, IRAM-30 m, APEX, and KPNO. We found that CO line profile prediction offers constraints on the shock velocity and pre-shock density independent of the absolute line brightness and requires fewer CO lines than diagnostics using a rotational excitation diagram.
Keywords