Dense gas in molecular clouds is an important signature of ongoing and future star formation. We identify and track dense cores in the Starforge simulations, following the core evolution from birth through dispersal by stellar feedback for typical Milky Way cloud conditions. Only ∼8% of cores host protostars, and most disperse before forming stars. The median starless and protostellar core lifetimes are ∼0.5–0.6 Myr and ∼0.8–1.1 Myr, respectively, where the protostellar phase lasts 0.1 Myr. While core evolution is stochastic, we find that virial ratios and line widths decline in prestellar cores, coincident with turbulent decay. Collapse occurs over ∼0.1 Myr, once the central density exceeds ≳106cm−3. Starless cores, only, follow line-width–size and mass–size relations,σ∝R0.3andM∝R1. The core median mass, radius, and velocity dispersion scale weakly with the cloud magnetic field strength. We cluster the core properties and find that protostellar cores have >80% likelihood of belonging to three particular groups that are characterized by high central densities, compact radii, and lower virial parameters. Overall, core evolution appears to be universally set by the interplay of gravity and magnetized turbulence, while stellar feedback dictates protostellar core properties and sets the protostellar phase lifetime.
more »
« less
3D radiative transfer modelling and virial analysis of starless cores in the B10 region of the Taurus molecular cloud
ABSTRACT Low-mass stars like our Sun begin their evolution within cold (10 K) and dense (∼105 cm−3) cores of gas and dust. The physical structure of starless cores is best probed by thermal emission of dust grains. We present a high-resolution dust continuum study of the starless cores in the B10 region of the Taurus Molecular Cloud. New observations at 1.2 and 2.0 mm (12 and 18 arcsec resolution) with the NIKA2 instrument on the IRAM 30m have probed the inner regions of 14 low-mass starless cores. We perform sophisticated 3D radiative transfer modelling for each of these cores through the radiative transfer framework pandora, which utilizes RADMC-3D. Model best-fits constrain each cores’ central density, density slope, aspect ratio, opacity, and interstellar radiation field strength. These ‘typical’ cores in B10 span central densities from 5 × 104 to 1 × 106 cm−3, with a mean value of 2.6 × 105 cm−3. We find the dust opacity laws assumed in the 3D modelling, as well as the estimates from Herschel, have dust emissivity indices, β’s, on the lower end of the distribution constrained directly from the NIKA2 maps, which averages to β = 2.01 ± 0.48. From our 3D density structures and archival NH3 data, we perform a self-consistent virial analysis to assess each core’s stability. Ignoring magnetic field contributions, we find nine out of the 14 cores (64 per cent) are either in virial equilibrium or are bound by gravity and external pressure. To push the bounded cores back to equilibrium, an effective magnetic field difference of only ∼15 $$\mu$$G is needed.
more »
« less
- Award ID(s):
- 2205474
- PAR ID:
- 10441799
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 521
- Issue:
- 3
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 4579 to 4597
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)ABSTRACT (Sub)millimetre dust opacities are required for converting the observable dust continuum emission to the mass, but their values have long been uncertain, especially in discs around young stellar objects. We propose a method to constrain the opacity κν in edge-on discs from a characteristic optical depth τ0,ν, the density ρ0, and radius R0 at the disc outer edge through κν = τ0,ν/(ρ0R0), where τ0,ν is inferred from the shape of the observed flux along the major axis, ρ0 from gravitational stability considerations, and R0 from direct imaging. We applied the 1D semi-analytical model to the embedded, Class 0, HH 212 disc, which has high-resolution data in Atacama Large Millimetre/submillimetre Array (ALMA) bands 9, 7, 6, and 3 and Very Large Array Ka band (λ = 0.43, 0.85, 1.3, 2.9, and 9.1 mm). The modelling is extended to 2D through RADMC-3D radiative transfer calculations. We find a dust opacity of κν ≈ 1.9 × 10−2, 1.3 × 10−2, and 4.9 × 10−3 cm2 g−1 of gas and dust for ALMA bands 7, 6, and 3, respectively, with uncertainties dependent on the adopted stellar mass. The inferred opacities lend support to the widely used prescription κλ = 2.3 × 10−2(1.3mm/λ) cm2 g−1 . We inferred a temperature of ∼45 K at the disc outer edge that increases radially inwards. It is well above the sublimation temperatures of ices such as CO and N2, which supports the notion that the disc chemistry cannot be completely inherited from the protostellar envelope.more » « less
-
Abstract We present new Atacama Large Millimeter/submillimeter Array (ALMA) continuum and NH2D, N2D+, and H2D+line emission at matched, ∼100 au resolution toward the dense star-forming cores SM1N and N6 within the Ophiuchus molecular cloud. We determine the density and temperature structure of SM1N based on radiative transfer modeling and simulated observations of the multiwavelength continuum emission at 0.8, 2, and 3 mm. We show that SM1N is best fit by either a broken power-law or Plummer-like density profile with high central densities (n∼ 108cm−3), and an inner transition radius of only ∼80–300 au. The free-fall time of the inner region is only a few ×103yr. The continuum modeling rules out the presence of an embedded first hydrostatic core (FHSC) or protostar. SM1N is therefore a dynamically unstable but still starless core. We find that NH2D is likely depleted at high densities within SM1N. The nonthermal velocity dispersions increase from NH2D to N2H+and H2D+, possibly tracing increasing (but still subsonic) infall speeds at higher densities as predicted by some models of starless core contraction. Toward N6, we confirm the previous ALMA detection of a faint, embedded point source (N6-mm) in 0.8 mm continuum emission. NH2D and N2D+avoid N6-mm within ∼100 au, while H2D+is not strongly detected toward N6. The distribution of these tracers is consistent with heating by a young, warm object. N6-mm thus remains one of the best candidate FHSCs detected so far, although its observed (sub)millimeter luminosity remains below predictions for FHSCs.more » « less
-
We present Atacama Large Millimeter/submillimeter Array Band 6 (1.3 mm) observations of dense cores in three massive molecular clouds within the central molecular zone (CMZ) of the Milky Way, including the Dust Ridge cloud e, Sgr C, and the 20 km s−1cloud, at a spatial resolution of 2000 au. Among the 834 cores identified from the 1.3 mm continuum, we constrain temperatures and linewidths of 253 cores using local thermodynamic equilibrium methods to fit the H2CO and/or CH3CN spectra. We determine their masses using the 1.3 mm dust continuum and derived temperatures, and then evaluate their virial parameters using the H2CO and/or CH3CN linewidths and construct the core mass functions (CMFs). We find that the contribution of external pressure is crucial for the virial equilibrium of the dense cores in the three clouds, which contrasts with the environment in the Galactic disk where dense cores are already bound, even without the contribution of external pressure. With our new temperature estimates we also find that the CMFs show a Salpeter-like slope in the high-mass (≳3–6M⊙) end, a change from previous works. Combined with the possible top-heavy initial mass functions (IMFs) in the CMZ, our result suggests that gas accretion and further fragmentation may play important roles in transforming the CMF to the IMF.more » « less
-
Abstract Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ 214μm observations of polarized thermal dust emission and high-resolution GBT-Argus C18O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion ofB-field orientations, we produce a map of the B-field strength of the IRDC, which exhibits values between ∼0.03 and 1 mG based on a refined Davis–Chandrasekhar–Fermi method proposed by Skalidis & Tassis. Comparing to a map of inferred density, the IRDC exhibits aB–nrelation with a power-law index of 0.51 ± 0.02, which is consistent with a scenario of magnetically regulated anisotropic collapse. Consideration of the mass-to-flux ratio map indicates that magnetic fields are dynamically important in most regions of the IRDC. A virial analysis of a sample of massive, dense cores in the IRDC, including evaluation of magnetic and kinetic internal and surface terms, indicates consistency with virial equilibrium, sub-Alfvénic conditions, and a dominant role forB-fields in regulating collapse. A clear alignment of magnetic field morphology with the direction of the steepest column density gradient is also detected. However, there is no preferred orientation of protostellar outflow directions with theB-field. Overall, these results indicate that magnetic fields play a crucial role in regulating massive star and star cluster formation, and therefore they need to be accounted for in theoretical models of these processes.more » « less
An official website of the United States government

