The SOFIA Massive (SOMA) Star Formation Survey. V. Clustered Protostars

We present $\sim8-40\,\mu$m SOFIA-FORCAST images of seven regions of ``clustered" star formation as part of the SOFIA Massive (SOMA) Star Formation Survey. We identify a total of 34 protostar candidates and build their spectral energy distributions (SEDs). We fit these SEDs with a grid of radiative transfer models based on the Turbulent Core Accretion (TCA) theory to derive key protostellar properties, including initial core mass, $M_c$, clump environment mass surface density, $\Sigma_{\rm cl}$, and current protostellar mass, $m_*$. We also carry out empirical graybody (GB) estimation of $\Sigma_{\rm cl}$, which allows a case of restricted SED fitting within the TCA model grid. We also release version 2.0 of the open-source Python package \emph{sedcreator}, designed to automate the aperture photometry and SED building and fitting process for sources in clustered environments, where flux contamination from close neighbors typically complicates the process. Using these updated methods, SED fitting yields values of $M_c\sim30-200\:M_{\odot}$, $\Sigma_{\text{cl,SED}}\sim0.1-3\:{\rm{g\:cm}}^{-2}$, and $m_*\sim4-50\:M_{\odot}$. The graybody fitting yields smaller values of $\Sigma_{\text{cl,GB}}\lesssim1\:{\rm{g\:cm}}^{-2}$. From these results, we do not find evidence for a critical $\Sigma_{\rm{cl}}$ needed to form massive ($\gtrsim 8\:M_\odot$) stars. However, we do find tentative evidence for a dearth of the most massive ($m_*\gtrsim30\:M_\odot$) protostars in the clustered regions suggesting a potential impact of environment on the stellar initial mass function.