skip to main content


Title: The SOFIA Massive (SOMA) Star Formation Survey. IV. Isolated Protostars
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
NSF-PAR ID:
10388093
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
942
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 7
Size(s):
["Article No. 7"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We present a catalog of 315 protostellar outflow candidates detected in SiOJ= 5 − 4 in the ALMA-IMF Large Program, observed with ∼2000 au spatial resolution, 0.339 km s−1velocity resolution, and 2–12 mJy beam−1(0.18–0.8 K) sensitivity. We find median outflow masses, momenta, and kinetic energies of ∼0.3M, 4Mkm s−1, and 1045erg, respectively. Median outflow lifetimes are 6000 yr, yielding median mass, momentum, and energy rates ofṀ= 10−4.4Myr−1,Ṗ= 10−3.2Mkm s−1yr−1, andĖ= 1L. We analyze these outflow properties in the aggregate in each field. We find correlations between field-aggregated SiO outflow properties and total mass in cores (∼3σ–5σ), and no correlations above 3σwith clump mass, clump luminosity, or clump luminosity-to-mass ratio. We perform a linear regression analysis and find that the correlation between field-aggregated outflow mass and total clump mass—which has been previously described in the literature—may actually be mediated by the relationship between outflow mass and total mass in cores. We also find that the most massive SiO outflow in each field is typically responsible for only 15%–30% of the total outflow mass (60% upper limit). Our data agree well with the established mechanical force−bolometric luminosity relationship in the literature, and our data extend this relationship up toL≥ 106LandṖ≥ 1Mkm s−1yr−1. Our lack of correlation with clumpL/Mis inconsistent with models of protocluster formation in which all protostars start forming at the same time.

     
    more » « less
  2. Abstract We present a Spitzer/Herschel focused survey of the Aquila molecular clouds ( d ∼ 436 pc) as part of the eHOPS (extension of the Herschel orion protostar survey, or HOPS, Out to 500 ParSecs) census of nearby protostars. For every source detected in the Herschel/PACS bands, the eHOPS-Aquila catalog contains 1–850 μ m SEDs assembled from the Two Micron All Sky Survey, Spitzer, Herschel, the Wide-field Infrared Survey Explorer, and James Clerk Maxwell Telescope/SCUBA-2 data. Using a newly developed set of criteria, we classify objects by their SEDs as protostars, pre-main-sequence stars with disks, and galaxies. A total of 172 protostars are found in Aquila, tightly concentrated in the molecular filaments that thread the clouds. Of these, 71 (42%) are Class 0 protostars, 54 (31%) are Class I protostars, 43 (25%) are flat-spectrum protostars, and four (2%) are Class II sources. Ten of the Class 0 protostars are young PACS bright red sources similar to those discovered in Orion. We compare the SEDs to a grid of radiative transfer models to constrain the luminosities, envelope densities, and envelope masses of the protostars. A comparison of the eHOPS-Aquila to the HOPS protostars in Orion finds that the protostellar luminosity functions in the two star-forming regions are statistically indistinguishable, the bolometric temperatures/envelope masses of eHOPS-Aquila protostars are shifted to cooler temperatures/higher masses, and the eHOPS-Aquila protostars do not show the decline in luminosity with evolution found in Orion. We briefly discuss whether these differences are due to biases between the samples, diverging star formation histories, or the influence of environment on protostellar evolution. 
    more » « less
  3. Abstract

    The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work, we focus on the most massive galaxies in the nearby universelog(M/M)>11.2. We obtain high-quality Magellan/LDSS-3 long-slit spectroscopy with a wide wavelength coverage of 0.4–1.01μm for 41 early-type galaxies (ETGs) in the MASSIVE survey and derive high signal-to-noise spectra within an aperture ofRe/8. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter,αIMF= (M/L)/(M/L)MW, ranges from 1.1 to 3.1, with an average of 〈αIMF〉 = 1.84, suggesting that on average, the IMF is more bottom heavy than Salpeter. Comparing the estimated stellar masses with the dynamical masses, we find that most galaxies have stellar masses that are smaller than their dynamical masses within the 1σuncertainty. We complement our sample with lower-mass galaxies from the literature and confirm thatlog(αIMF)is positively correlated withlog(σ),log(M), andlog(Mdyn). From the combined sample, we show that the IMF in the centers of more massive ETGs is more bottom heavy. In addition, we find thatlog(αIMF)is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density ΣKroupamight be responsible for the variation ofαIMF. We conclude thatσ, [Mg/Fe], and [Z/H] are the primary drivers of the global stellar IMF variation.

     
    more » « less
  4. Abstract

    The unprecedented angular resolution and sensitivity of the Atacama Large Millimeter/submillimeter Array make it possible to unveil disk populations in distant (>2 kpc), embedded young cluster environments. We have conducted an observation toward the central region of the massive protocluster G286.21+0.16 at 1.3 mm. With a spatial resolution of 23 mas and a sensitivity of 15μJy beam−1, we detect a total of 38 protostellar disks. These disks have dust masses ranging from about 53 to 1825M, assuming a dust temperature of 20 K. This sample is not closely associated with previously identified dense cores, as would be expected for disks around Class 0 protostars. Thus, we expect our sample, being flux-limited, to be mainly composed of Class I/flat-spectrum source disks, since these are typically more massive than Class II disks. Furthermore, we find that the distributions of disk masses and radii are statistically indistinguishable from those of the Class I/flat-spectrum objects in the Orion molecular cloud, indicating that similar processes are operating in G286.21+0.16 to regulate disk formation and evolution. The cluster center appears to host a massive protostellar system composed of three sources within 1200 au, including a potential binary with 600 au projected separation. Relative to this center, there is no evidence for widespread mass segregation in the disk population. We do find a tentative trend of increasing disk radius versus distance from the cluster center, which may point to the influence of dynamical interactions being stronger in the central regions.

     
    more » « less
  5. Abstract

    We present Stratospheric Observatory For Infrared Astronomy (SOFIA) + Atacama Large Millimeter/submillimeter Array (ALMA) continuum and spectral-line polarization data on the massive molecular cloud BYF 73, revealing important details about the magnetic field morphology, gas structures, and energetics in this unusual massive star formation laboratory. The 154μm HAWC+ polarization map finds a highly organized magnetic field in the densest, inner 0.55 × 0.40 pc portion of the cloud, compared to an unremarkable morphology in the cloud’s outer layers. The 3 mm continuum ALMA polarization data reveal several more structures in the inner domain, including a parsec-long, ∼500M“Streamer” around the central massive protostellar object MIR 2, with magnetic fields mostly parallel to the east–west Streamer but oriented north–south across MIR 2. The magnetic field orientation changes from mostly parallel to the column density structures to mostly perpendicular, at thresholdsNcrit= 6.6 × 1026m−2,ncrit= 2.5 × 1011m−3, andBcrit= 42 ± 7 nT. ALMA also mapped Goldreich–Kylafis polarization in12CO across the cloud, which traces, in both total intensity and polarized flux, a powerful bipolar outflow from MIR 2 that interacts strongly with the Streamer. The magnetic field is also strongly aligned along the outflow direction; energetically, it may dominate the outflow near MIR 2, comprising rare evidence for a magnetocentrifugal origin to such outflows. A portion of the Streamer may be in Keplerian rotation around MIR 2, implying a gravitating mass 1350 ± 50Mfor the protostar+disk+envelope; alternatively, these kinematics can be explained by gas in free-fall toward a 950 ± 35Mobject. The high accretion rate onto MIR 2 apparently occurs through the Streamer/disk, and could account for ∼33% of MIR 2's total luminosity via gravitational energy release.

     
    more » « less