2025ApJ...986...15T

We present ∼8-40 μm SOFIA-FORCAST images of seven regions of “clustered” star formation as part of the SOFIA Massive 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, Σ_cl, and current protostellar mass, m_*. We also carry out empirical graybody (GB) estimation of Σ_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 sedcreator, which is 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 ∼ 30-200 M_⊙, Σ_cl,SED ∼ 0.1-3 g cm^–2, and m_* ∼ 4-50 M_⊙. The GB fitting yields smaller values of Σ_cl,GB ≲ 1 g cm^–2. From these results, we do not find evidence for a critical Σ_cl needed to form massive (≳8 M_⊙) stars. However, we do find tentative evidence for a dearth of the most massive (m_* ≳ 30 M_⊙) protostars in the clustered regions, suggesting a potential impact of environment on the stellar initial mass function.