The SOFIA Massive (SOMA) Star Formation Survey. IV. Isolated Protostars

Ruben Fedriani, Jonathan C. Tan, Zoie Telkamp, Yichen Zhang, Yao-Lun Yang, Mengyao Liu, Chi-Yan Law, Maria T. Beltran, Viviana Rosero, Kei E. I. Tanaka, Giuliana Cosentino, Prasanta Gorai, Juan Farias, Jan E. Staff, James M. De Buizer, Barbara Whitney

Published 2022-05-23Version 1

We present $\sim10-40\,\mu$m \textit{SOFIA}-FORCAST images of 11 "isolated" protostars as part of the \textit{SOFIA} Massive (SOMA) Star Formation Survey, with this morphological classification based on 37\,$\mu$m imaging. We develop an automated method to define source aperture size based on the gradient of its background-subtracted enclosed flux and apply this to build spectral energy distributions (SEDs). We fit the SEDs with radiative transfer models, based on the Turbulent Core Accretion (TCA) theory, to estimate key protostellar properties. Here we release the \textit{sedcreator} python package that carries out these methods. The SEDs are generally well-fit by the TCA models, from which we infer initial core masses $M_c$ ranging from $50-430\:M_\odot$, clump mass surface densities $\Sigma_{\rm cl}\sim0.1-3\:{\rm{g\:cm}}^{-2}$ and current protostellar masses $m_*\sim2-40\:M_\odot$. From an uniform analysis of the 40 sources in the full SOMA survey to date, we find that massive protostars form across a wide range of clump mass surface density environments, placing constraints on theories that predict a minimum threshold $\Sigma_{\rm cl}$ for massive star formation. However, the upper end of the $m_*-\Sigma_{\rm cl}$ distribution follows trends predicted by models of internal protostellar feedback that find greater star formation efficiency in higher $\Sigma_{\rm cl}$ conditions. We also investigate protostellar FIR variability by comparison with IRAS data, finding no significant variation over a $\sim$40-year baseline.