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.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Friday, May 16 until 2:00 AM ET on Saturday, May 17 due to maintenance. We apologize for the inconvenience.


Title: The Effects of Magnetic Fields and Outflow Feedback on the Shape and Evolution of the Density Probability Distribution Function in Turbulent Star-forming Clouds
Abstract Using a suite of 3D hydrodynamical simulations of star-forming molecular clouds, we investigate how the density probability distribution function (PDF) changes when including gravity, turbulence, magnetic fields, and protostellar outflows and heating. We find that the density PDF is not lognormal when outflows and self-gravity are considered. Self-gravity produces a power-law tail at high densities, and the inclusion of stellar feedback from protostellar outflows and heating produces significant time-varying deviations from a lognormal distribution at low densities. The simulation with outflows has an excess of diffuse gas compared to the simulations without outflows, exhibits an increased average sonic Mach number, and maintains a slower star formation rate (SFR) over the entire duration of the run. We study the mass transfer between the diffuse gas in the lognormal peak of the PDF, the collapsing gas in the power-law tail, and the stars. We find that the mass fraction in the power-law tail is constant, such that the stars form out of the power-law gas at the same rate at which the gas from the lognormal part replenishes the power law. We find that turbulence does not provide significant support in the dense gas associated with the power-law tail. When including outflows and magnetic fields in addition to driven turbulence, the rate of mass transfer from the lognormal to the power law, and then to the stars, becomes significantly slower, resulting in slower SFRs and longer depletion times.  more » « less
Award ID(s):
2009679
PAR ID:
10438218
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
927
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
75
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high- and low-density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations. 
    more » « less
  2. null (Ed.)
    ABSTRACT The velocity gradients technique (VGT) and the probability density functions (PDFs) of mass density are tools to study turbulence, magnetic fields, and self-gravity in molecular clouds. However, self-absorption can significantly make the observed intensity different from the column density structures. In this work, we study the effects of self-absorption on the VGT and the intensity PDFs utilizing three synthetic emission lines of CO isotopologues 12CO (1–0), 13CO (1–0), and C18O (1–0). We confirm that the performance of VGT is insensitive to the radiative transfer effect. We numerically show the possibility of constructing 3D magnetic fields tomography through VGT. We find that the intensity PDFs change their shape from the pure lognormal to a distribution that exhibits a power-law tail depending on the optical depth for supersonic turbulence. We conclude the change of CO isotopologues’ intensity PDFs can be independent of self-gravity, which makes the intensity PDFs less reliable in identifying gravitational collapsing regions. We compute the intensity PDFs for a star-forming region NGC 1333 and find the change of intensity PDFs in observation agrees with our numerical results. The synergy of VGT and the column density PDFs confirms that the self-gravitating gas occupies a large volume in NGC 1333. 
    more » « less
  3. Abstract We have comprehensively studied the multiscale physical properties of the massive infrared dark cloud G28.34 (the Dragon cloud) with dust polarization and molecular line data from Planck, FCRAO-14 m, James Clerk Maxwell Telescope, and Atacama Large Millimeter/submillimeter Array. We find that the averaged magnetic fields of clumps tend to be either parallel with or perpendicular to the cloud-scale magnetic fields, while the cores in clump MM4 tend to have magnetic fields aligned with the clump fields. Implementing the relative orientation analysis (for magnetic fields, column density gradients, and local gravity), velocity gradient technique, and modified Davis–Chandrasekhar–Fermi analysis, we find that G28.34 is located in a trans-to-sub-Alfvénic environment; the magnetic field is effectively resisting gravitational collapse in large-scale diffuse gas, but is distorted by gravity within the cloud and affected by star formation activities in high-density regions, and the normalized mass-to-flux ratio tends to increase with increasing density and decreasing radius. Considering the thermal, magnetic, and turbulent supports, we find that the environmental gas of G28.34 is in a supervirial (supported) state, the infrared dark clumps may be in a near-equilibrium state, and core MM4-core4 is in a subvirial (gravity-dominant) state. In summary, we suggest that magnetic fields dominate gravity and turbulence in the cloud environment at large scales, resulting in relatively slow cloud formation and evolution processes. Within the cloud, gravity could overwhelm both magnetic fields and turbulence, allowing local dynamical star formation to happen. 
    more » « less
  4. ABSTRACT Simulations of isolated giant molecular clouds (GMCs) are an important tool for studying the dynamics of star formation, but their turbulent initial conditions (ICs) are uncertain. Most simulations have either initialized a velocity field with a prescribed power spectrum on a smooth density field (failing to model the full structure of turbulence) or ‘stirred’ turbulence with periodic boundary conditions (which may not model real GMC boundary conditions). We develop and test a new GMC simulation setup (called turbsphere) that combines advantages of both approaches: we continuously stir an isolated cloud to model the energy cascade from larger scales, and use a static potential to confine the gas. The resulting cloud and surrounding envelope achieve a quasi-equilibrium state with the desired hallmarks of supersonic ISM turbulence (e.g. density PDF and a ∼k−2 velocity power spectrum), whose bulk properties can be tuned as desired. We use the final stirred state as initial conditions for star formation simulations with self-gravity, both with and without continued driving and protostellar jet feedback, respectively. We then disentangle the respective effects of the turbulent cascade, simulation geometry, external driving, and gravity/MHD boundary conditions on the resulting star formation. Without external driving, the new setup obtains results similar to previous simple spherical cloud setups, but external driving can suppress star formation considerably in the new setup. Periodic box simulations with the same dimensions and turbulence parameters form stars significantly slower, highlighting the importance of boundary conditions and the presence or absence of a global collapse mode in the results of star formation calculations. 
    more » « less
  5. Abstract We present Atacama Large Millimeter/submillimeter Array observations of the ∼10,000 au environment surrounding 21 protostars in the Orion A molecular cloud tracing outflows. Our sample is composed of Class 0 to flat-spectrum protostars, spanning the full ∼1 Myr lifetime. We derive the angular distribution of outflow momentum and energy profiles and obtain the first two-dimensional instantaneous mass, momentum, and energy ejection rate maps using our new approach: the pixel flux-tracing technique. Our results indicate that by the end of the protostellar phase, outflows will remove ∼2–4 M ⊙ from the surrounding ∼1 M ⊙ low-mass core. These high values indicate that outflows remove a significant amount of gas from their parent cores and continuous core accretion from larger scales is needed to replenish core material for star formation. This poses serious challenges to the concept of cores as well-defined mass reservoirs , and hence to the simplified core-to-star conversion prescriptions. Furthermore, we show that cavity opening angles, and momentum and energy distributions all increase with protostar evolutionary stage. This is clear evidence that even garden-variety protostellar outflows: (a) effectively inject energy and momentum into their environments on 10,000 au scales, and (b) significantly disrupt their natal cores, ejecting a large fraction of the mass that would have otherwise fed the nascent star. Our results support the conclusion that protostellar outflows have a direct impact on how stars get their mass, and that the natal sites of individual low-mass star formation are far more dynamic than commonly accepted theoretical paradigms. 
    more » « less