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
What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates
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- NSF-PAR ID:
- 10446938
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 954
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 93
- Size(s):
- ["Article No. 93"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
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
Abstract Magnetic fields of molecular clouds in the central molecular zone (CMZ) have been relatively under-observed at sub-parsec resolution. Here, we report JCMT/POL2 observations of polarized dust emission in the CMZ, which reveal magnetic field structures in dense gas at ∼0.5 pc resolution. The 11 molecular clouds in our sample include two in the western part of the CMZ (Sgr C and a farside cloud candidate), four around the Galactic longitude 0 (the 50 km s−1cloud, CO 0.02−0.02, the
Stone , and theSticks and Straw among the Three Little Pigs), and five along the Dust Ridge (G0.253+0.016, clouds b, c, d, and e/f), for each of which we estimate the magnetic field strength using the angular dispersion function method. The morphologies of magnetic fields in the clouds suggest potential imprints of feedback from expanding Hii regions and young massive star clusters. A moderate correlation between the total viral parameter versus the star formation rate (SFR) and the dense gas fraction of the clouds is found. A weak correlation between the mass-to-flux ratio and the SFR, and a weak anticorrelation between the magnetic field and the dense gas fraction are also found. Comparisons between magnetic fields and other dynamic components in clouds suggest a more dominant role of self-gravity and turbulence in determining the dynamical states of the clouds and affecting star formation at the studied scales. -
null (Ed.)Context. Recent surveys of the Galactic plane in the dust continuum and CO emission lines reveal that large (≳50 pc) and massive (≳10 5 M ⊙ ) filaments, know as giant molecular filaments (GMFs), may be linked to Galactic dynamics and trace the mid-plane of the gravitational potential in the Milky Way. Yet our physical understanding of GMFs is still poor. Aims. We investigate the dense gas properties of one GMF, with the ultimate goal of connecting these dense gas tracers with star formation processes in the GMF. Methods. We imaged one entire GMF located at l ~ 52–54° longitude, GMF54 (~68 pc long), in the empirical dense gas tracers using the HCN(1–0), HNC(1–0), and HCO + (1–0) lines, and their 13 C isotopologue transitions, as well as the N 2 H + (1–0) line. We studied the dense gas distribution, the column density probability density functions (N-PDFs), and the line ratios within the GMF. Results. The dense gas molecular transitions follow the extended structure of the filament with area filling factors between 0.06 and 0.28 with respect to 13 CO(1–0). We constructed the N-PDFs of H 2 for each of the dense gas tracers based on their column densities and assumed uniform abundance. The N-PDFs of the dense gas tracers appear curved in log–log representation, and the HCO + N-PDF has the flattest power-law slope index. Studying the N-PDFs for sub-regions of GMF54, we found an evolutionary trend in the N-PDFs that high-mass star-forming and photon-dominated regions have flatter power-law indices. The integrated intensity ratios of the molecular lines in GMF54 are comparable to those in nearby galaxies. In particular, the N 2 H + / 13 CO ratio, which traces the dense gas fraction, has similar values in GMF54 and all nearby galaxies except Ultraluminous Infrared Galaxies. Conclusions. As the largest coherent cold gaseous structure in our Milky Way, GMFs, are outstanding candidates for connecting studies of star formation on Galactic and extragalactic scales. By analyzing a complete map of the dense gas in a GMF we have found that: (1) the dense gas N-PDFs appear flatter in more evolved regions and steeper in younger regions, and (2) its integrated dense gas intensity ratios are similar to those of nearby galaxies.more » « less
-
ABSTRACT In typical environments of star-forming clouds, converging supersonic turbulence generates shock-compressed regions, and can create strongly magnetized sheet-like layers. Numerical magnetohydrodynamic simulations show that within these post-shock layers, dense filaments and embedded self-gravitating cores form via gathering material along the magnetic field lines. As a result of the preferred-direction mass collection, a velocity gradient perpendicular to the filament major axis is a common feature seen in simulations. We show that this prediction is in good agreement with recent observations from the CARMA Large Area Star Formation Survey (CLASSy), from which we identified several filaments with prominent velocity gradients perpendicular to their major axes. Highlighting a filament from the north-west part of Serpens South, we provide both qualitative and quantitative comparisons between simulation results and observational data. In particular, we show that the dimensionless ratio Cv ≡ Δvh2/(GM/L), where Δvh is half of the observed perpendicular velocity difference across a filament, and M/L is the filament’s mass per unit length, can distinguish between filaments formed purely due to turbulent compression and those formed due to gravity-induced accretion. We conclude that the perpendicular velocity gradient observed in the Serpens South north-west filament can be caused by gravity-induced anisotropic accretion of material from a flattened layer. Using synthetic observations of our simulated filaments, we also propose that a density-selection effect may explain observed subfilaments (one filament breaking into two components in velocity space) as reported in recent observations.more » « less
