The SOMA MM Survey. I. An Astrochemical Census of Massive Protostars

Abstract

During massive star formation, dense gas undergoes chemical evolution, producing both simple and complex organic molecules (COMs) characteristic of hot molecular cores. How this evolution depends on protostellar physical properties remains unclear. We investigate the chemical content of 22 well-studied massive protostars from the SOFIA Massive (SOMA) Star Formation survey, aiming to identify correlations between chemical and physical parameters. We analyzed Atacama Compact Array and Total Power 1.3 mm (Band 6) data, deriving column densities, line widths, and excitation temperatures of multiple molecular species by modeling detected lines under local thermodynamic equilibrium (LTE) using MADCUBA. Spectra show 35 species, from simple molecules (e.g., CO, SO, SiO) to complex organic molecules (COMs), with seven sources exhibiting high chemical complexity (> 100 transitions). Average excitation temperatures vary across the sample: $T_\text{ex}>100~\text{K}$ for eight sources, $50-100~\text{K}$ for four, and $T_\text{ex} < 50~\text{K}$ for the remainder. Sources with $T_\text{ex} < 50~\text{K}$ trace lukewarm, chemically simple gas, while those with $T_\text{ex}>100~\text{K}$ indicate the presence of typical hot cores where thermal desorption is efficient, resulting in line-rich spectra. Comparing these chemical properties with the bolometric luminosity to envelope mass ratio ($L_\text{bol}/M_\text{env}$), an evolutionary tracer, we find tentative correlations with line widths, excitation temperature, and column densities. These data provide important constraints for chemodynamical models of massive protostellar cores.