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

TL;DR

This study uses SOFIA-FORCAST imaging and SED modeling to analyze 34 protostar candidates in clustered star-forming regions, revealing properties and environmental influences on massive star formation.

Contribution

It introduces an improved SED fitting approach with the sedcreator package and provides new insights into the environmental factors affecting massive star formation.

Findings

Protostar properties span initial core masses of 30-200 M_sun and current masses of 4-50 M_sun.

No evidence found for a critical clump surface density for massive star formation.

Potential environmental impact on the upper end of the stellar initial mass function.

Abstract

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.