skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on June 3, 2026

Title: The SOFIA Massive (SOMA) Star Formation Survey. V. Clustered Protostars
Abstract 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,Mc, 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 packagesedcreator, 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 ofMc∼ 30–200M, Σcl,SED∼ 0.1–3 g cm−2, andm*∼ 4–50M. The GB fitting yields smaller values of Σcl,GB≲ 1 g cm−2. From these results, we do not find evidence for a critical Σclneeded to form massive (≳8M) stars. However, we do find tentative evidence for a dearth of the most massive (m*≳ 30M) protostars in the clustered regions, suggesting a potential impact of environment on the stellar initial mass function.  more » « less
Award ID(s):
2206437
PAR ID:
10642481
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Inst. of Physics
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
986
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
15
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present ∼10–40μm SOFIA-FORCAST images of 11isolatedprotostars as part of the SOFIA Massive (SOMA) Star Formation Survey, with this morphological classification based on 37μm imaging. We develop an automated method to define source aperture size using 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, developed within the framework of turbulent core accretion (TCA) theory, to estimate key protostellar properties. Here, we release the sedcreator python package that carries out these methods. The SEDs are generally well fitted by the TCA models, from which we infer initial core massesMcranging from 20–430M, clump mass surface densities Σcl∼ 0.3–1.7 g cm−2, and current protostellar massesm*∼ 3–50M. From a 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 Σclfor massive star formation. However, the upper end of them*−Σcldistribution follows trends predicted by models of internal protostellar feedback that find greater star formation efficiency in higher Σclconditions. We also investigate protostellar far-IR variability by comparison with IRAS data, finding no significant variation over an ∼40 yr baseline. 
    more » « less
  2. Abstract We report high-resolution 1.3 mm continuum and molecular line observations of the massive protostar G28.20-0.05 with Atacama Large Millimeter/submillimeter Array. The continuum image reveals a ring-like structure with 2000 au radius, similar to morphology seen in archival 1.3 cm Very Large Array observations. Based on its spectral index and associated H30αemission, this structure mainly traces ionized gas. However, there is evidence for ∼30Mof dusty gas near the main millimeter continuum peak on one side of the ring, as well as in adjacent regions within 3000 au. A virial analysis on scales of ∼2000 au from hot core line emission yields a dynamical mass of ∼80M. A strong velocity gradient in the H30αemission is evidence for a rotating, ionized disk wind, which drives a larger-scale molecular outflow. An infrared spectral energy distribution (SED) analysis indicates a current protostellar mass ofm*∼ 40Mforming from a core with initial massMc∼ 300Min a clump with mass surface density of Σcl∼ 0.8 g cm−2. Thus the SED and other properties of the system can be understood in the context of core accretion models. A structure-finding analysis on the larger-scale continuum image indicates G28.20-0.05 is forming in a relatively isolated environment, with no other concentrated sources, i.e., protostellar cores, above ∼1Mfound from ∼0.1 to 0.4 pc around the source. This implies that a massive star can form in relative isolation, and the dearth of other protostellar companions within the ∼1 pc environs is a strong constraint on massive star formation theories that predict the presence of a surrounding protocluster. 
    more » « less
  3. Abstract We present James Webb Space Telescope (JWST) Near Infrared Camera observations of the massive star-forming molecular cloud Sagittarius C (Sgr C) in the Central Molecular Zone (CMZ). In conjunction with ancillary mid-IR and far-IR data, we characterize the two most massive protostars in Sgr C via spectral energy distribution (SED) fitting, estimating that they each have current masses ofm*∼ 20Mand surrounding envelope masses of ∼100M. We report a census of lower-mass protostars in Sgr C via a search for infrared counterparts to millimeter continuum dust cores found with the Atacama Large Millimeter/submillimeter Array (ALMA). We identify 88 molecular hydrogen outflow knot candidates originating from outflows from protostars in Sgr C, the first such unambiguous detections in the infrared in the CMZ. About a quarter of these are associated with flows from the two massive protostars in Sgr C; these extend for over 1 pc and are associated with outflows detected in ALMA SiO line data. An additional ∼40 features likely trace shocks in outflows powered by lower-mass protostars throughout the cloud. We report the discovery of a new star-forming region hosting two prominent bow shocks and several other line-emitting features driven by at least two protostars. We infer that one of these is forming a high-mass star given an SED-derived mass ofm*∼ 9Mand associated massive (∼90M) millimeter core and water maser. Finally, we identify a population of miscellaneous molecular hydrogen objects that do not appear to be associated with protostellar outflows. 
    more » « less
  4. Abstract Protostellar disks are an ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks large program, we present high angular resolution dust continuum (∼40 mas) and molecular line (∼150 mas) observations of the Class 0 protostar IRAS 15398–3359. The dust continuum is small, compact, and centrally peaked, while more extended dust structures are found in the outflow directions. We perform a 2D Gaussian fitting and find the deconvolved size and 2σradius of the dust disk to be 4.5 × 2.8 au and 3.8 au, respectively. We estimate the gas+dust disk mass assuming optically thin continuum emission to be 0.6MJ–1.8MJ, indicating a very low mass disk. The CO isotopologues trace components of the outflows and inner envelope, while SO traces a compact, rotating disk-like component. Using several rotation curve fittings on the position–velocity diagram of the SO emission, the lower limits of the protostellar mass and gas disk radius are 0.022Mand 31.2 au, respectively, from our Modified 2 single power-law fitting. A conservative upper limit of the protostellar mass is inferred to be 0.1M. The protostellar mass accretion rate and the specific angular momentum at the protostellar disk edge are found to be in the range of (1.3–6.1) × 10−6Myr−1and (1.2–3.8) × 10−4km s−1pc, respectively, with an age estimated between 0.4 × 104yr and 7.5 × 104yr. At this young age with no clear substructures in the disk, planet formation would likely not yet have started. This study highlights the importance of high-resolution observations and systematic fitting procedures when deriving dynamical properties of deeply embedded Class 0 protostars. 
    more » « less
  5. Abstract We analyze the first cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modeled as a constant fractional energy dissipation (fdiss= 0.75). In this paper, we examine the properties of dwarf galaxies withM*∼ 105–109Min both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm2g−1. On the contrary, models with (σ/m) = 10 cm2g−1produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm2g−1are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the Hirotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with (σ/m) ≳ 0.1 cm2g−1,fdiss= 0.75 are in tension in massive dwarfs (Mhalo∼ 1011M) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures. 
    more » « less