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.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: 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:
10404115
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
521
Issue:
3
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 4579-4597
Size(s):
p. 4579-4597
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    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
  2. ; ; ; ; ; (, The Astrophysical Journal)
    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
  3. ; ; ; ; ; ; ; ; ; ; et al (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT Young massive clusters (YMCs) are compact (≲1 pc), high-mass (>104 M⊙) stellar systems of significant scientific interest. Due to their rarity and rapid formation, we have very few examples of YMC progenitor gas clouds before star formation has begun. As a result, the initial conditions required for YMC formation are uncertain. We present high resolution (0.13 arcsec, ∼1000 au) ALMA observations and Mopra single-dish data, showing that Galactic Centre dust ridge ‘Cloud d’ (G0.412 + 0.052, mass = 7.6 × 104 M⊙, radius = 3.2 pc) has the potential to become an Arches-like YMC (104 M⊙, r ∼ 1 pc), but is not yet forming stars. This would mean it is the youngest known pre-star-forming massive cluster and therefore could be an ideal laboratory for studying the initial conditions of YMC formation. We find 96 sources in the dust continuum, with masses ≲3 M⊙ and radii of ∼103 au. The source masses and separations are more consistent with thermal rather than turbulent fragmentation. It is not possible to unambiguously determine the dynamical state of most of the sources, as the uncertainty on virial parameter estimates is large. We find evidence for large-scale (∼1 pc) converging gas flows, which could cause the cloud to grow rapidly, gaining 104 M⊙ within 105 yr. The highest density gas is found at the convergent point of the large-scale flows. We expect this cloud to form many high-mass stars, but find no high-mass starless cores. If the sources represent the initial conditions for star formation, the resulting initial mass function will be bottom heavy. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al (, The Astrophysical Journal)
    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
  5. ; ; ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT We study the formation, evolution, and collapse of dense cores by tracking structures in a magnetohydrodynamic simulation of a star-forming cloud. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including Neural Gas prototype learning and Fuzzy c-means clustering to analyse the density and velocity dispersion profiles of cores together with six bulk properties. We produce a 2-d visualization using a Uniform Manifold Approximation and Projection (UMAP), which facilitates the connection between physical properties and three partially-overlapping phases: i) unbound turbulent structures (Phase I), ii) coherent cores that have low turbulence (Phase II), and iii) bound cores, many of which become protostellar (Phase III). Within Phase II, we identify a population of long-lived coherent cores that reach a quasi-equilibrium state. Most prestellar cores form in Phase II and become protostellar after evolving into Phase III. Due to the turbulent cloud environment, the initial core properties do not uniquely predict the eventual evolution, i.e. core evolution is stochastic, and cores follow no one evolutionary path. The phase lifetimes are 1.0 ± 0.1 × 105 yr, 1.3 ± 0.2 × 105 yr, and 1.8 ± 0.3 × 105 yr for Phase I, II, and III, respectively. We compare our results to NH3 observations of dense cores. Known coherent cores predominantly map into Phase II, while most turbulent pressure-confined cores map to Phase I or III. We predict that a significant fraction of observed starless cores have unresolved coherent regions and that ≳20 per cent of observed starless cores will not form stars. Measurements of core radial profiles in addition to the usual bulk properties will enable more accurate predictions of core evolution. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email