Abstract
We perform ideal magnetohydrodynamics high-resolution adaptive mesh refinement simulations with driven turbulence and self-gravity and find that long filamentary molecular clouds are formed ...at the converging locations of large-scale turbulence flows and the filaments are bounded by gravity. The magnetic field helps shape and reinforce the long filamentary structures. The main filamentary cloud has a length of ∼4.4 pc. Instead of a monolithic cylindrical structure, the main cloud is shown to be a collection of fibre/web-like substructures similar to filamentary clouds such as L1495. Unless the line-of-sight is close to the mean field direction, the large-scale magnetic field and striations in the simulation are found roughly perpendicular to the long axis of the main cloud, similar to L1495. This provides strong support for a large-scale moderately strong magnetic field surrounding L1495. We find that the projection effect from observations can lead to incorrect interpretations of the true three-dimensional physical shape, size, and velocity structure of the clouds. Helical magnetic field structures found around filamentary clouds that are interpreted from Zeeman observations can be explained by a simple bending of the magnetic field that pierces through the cloud. We demonstrate that two dark clouds form a T-shaped configuration that is strikingly similar to the infrared dark cloud SDC13, leading to the interpretation that SDC13 results from a collision of two long filamentary clouds. We show that a moderately strong magnetic field (${{\cal M}_{\rm A}}\sim 1$) is crucial for maintaining a long and slender filamentary cloud for a long period of time ∼0.5 Myr.
ABSTRACT
While magnetic fields are important in contemporary star formation, their role in primordial star formation is unknown. Magnetic fields of the order of 10−16 G are produced by the Biermann ...battery due to the curved shocks and turbulence associated with the infall of gas into the dark matter minihaloes that are the sites of formation of the first stars. These fields are rapidly amplified by a small-scale dynamo until they saturate at or near equipartition with the turbulence in the central region of the gas. Analytical results are given for the outcome of the dynamo, including the effect of compression in the collapsing gas. The mass-to-flux ratio in this gas is two to three times the critical value, comparable to that in contemporary star formation. Predictions of the outcomes of simulations using smooth particle hydrodynamics (SPH) and grid-based adaptive mesh refinement are given. Because the numerical viscosity and resistivity for the standard resolution of 64 cells per Jeans length are several orders of magnitude greater than the physical values, dynamically significant magnetic fields affect a much smaller fraction of the mass in simulations than in reality. An appendix gives an analytical treatment of free-fall collapse, including that in a constant-density background. Another appendix presents a new method of estimating the numerical viscosity; results are given for both SPH and grid-based codes.
Abstract
Star formation in a filamentary infrared dark cloud (IRDC) is simulated over the dynamic range of 4.2 pc to 28 au for a period of 3.5 × 105 yr, including magnetic fields and both radiative ...and outflow feedback from the protostars. At the end of the simulation, the star formation efficiency is 4.3 per cent and the star formation rate per free-fall time is εff ≃ 0.04, within the range of observed values. The total stellar mass increases as ∼t2, whereas the number of protostars increases as ∼t1.5. We find that the density profile around most of the simulated protostars is ∼ρ ∝ r−1.5. At the end of the simulation, the protostellar mass function approaches the Chabrier stellar initial mass function. We infer that the time to form a star of median mass 0.2 M⊙ is about 1.4 × 105 yr from the median mass accretion rate. We find good agreement among the protostellar luminosities observed in the large sample of Dunham et al., our simulation and a theoretical estimate, and we conclude that the classical protostellar luminosity problem is resolved. The multiplicity of the stellar systems in the simulation agrees, to within a factor of 2, with observations of Class I young stellar objects; most of the simulated multiple systems are unbound. Bipolar protostellar outflows are launched using a subgrid model, and extend up to 1 pc from their host star. The mass–velocity relation of the simulated outflows is consistent with both observation and theory.
Similar to their low-mass counterparts, massive stars likely form via the collapse of prestellar molecular cores. Recent observations suggest that most massive cores are subvirial (i.e., not ...supported by turbulence) and therefore are likely unstable to gravitational collapse. Here we perform radiation-hydrodynamic simulations to follow the collapse of turbulent massive prestellar cores with subvirial and virialized initial conditions to explore how their dynamic state affects the formation of massive stars and core fragmentation into companion stars. We find that subvirial cores undergo rapid monolithic collapse, resulting in higher accretion rates at early times as compared to the collapse of virialized cores that have the same physical properties. In contrast, we find that virialized cores undergo a slower, gradual collapse and significant turbulent fragmentation at early times, resulting in numerous companion stars. In the absence of strong magnetic fields and protostellar outflows, we find that the faster growth rate of massive stars that are born out of subvirial cores leads to an increase in the radiative heating of the core, thereby further suppressing fragmentation at early times when turbulent fragmentation occurs for virialized cores. Regardless of initial condition, we find that the massive accretion disks that form around massive stars dominant the accretion flow onto the star at late times and eventually become gravitationally unstable and fragment to form companion stars at late times.
The most accurate measurements of magnetic fields in star-forming gas are based on the Zeeman observations analysed by Crutcher et al. We show that their finding that the 3D magnetic field scales ...approximately as density0.65 can also be obtained from analysis of the observed line-of-sight fields. We present two large-scale adaptive-mesh-refinement magnetohydrodynamic simulations of several thousand M⊙ of turbulent, isothermal, self-gravitating gas, one with a strong initial magnetic field (Alfv
$\acute{\rm e}$
n Mach number
${{\cal M}_{\rm A,0}}=1$
) and one with a weak initial field (
${{\cal M}_{\rm A,0}}=10$
). We construct samples of the 100 most massive clumps in each simulation and show that they exhibit a power-law relation between field strength and density (
$\bar{n}_{\rm H}$
) in excellent agreement with the observed one. Our results imply that the average field in molecular clumps in the interstellar medium (ISM) is
${\langle B_{\rm tot}(\bar{n}_{\rm H}) \rangle }\approx 42 \bar{n}_{\rm H,\, 4}^{0.65} \,\mu {\rm G}$
. Furthermore, the median value of the ratio of the line-of-sight field to density0.65 in the simulations is within a factor of about (1.3, 1.7) of the observed value for the strong- and weak-field cases, respectively. The median value of the mass-to-flux ratio, normalized to the critical value, is 70 per cent of the line-of-sight value. This is larger than the 50 per cent usually cited for spherical clouds because the actual mass-to-flux ratio depends on the volume-weighted field, whereas the observed one depends on the mass-weighted field. Our results indicate that the typical molecular clump in the ISM is significantly supercritical (∼ factor of 3). The results of our strong-field model are in very good quantitative agreement with the observations of Li et al., which show a strong correlation in field orientation between small and large scales. Because there is a negligible correlation in the weak-field model, we conclude that molecular clouds form from strongly magnetized (although magnetically supercritical) gas, in agreement with the conclusion of Li et al.
ABSTRACT
Beginning with cosmological initial conditions at z = 100, we simulate the effects of magnetic fields on the formation of Population III stars and compare our results with the predictions of ...Paper I. We use gadget-2 to follow the evolution of the system while the field is weak. We introduce a new method for treating kinematic fields by tracking the evolution of the deformation tensor. The growth rate in this stage of the simulation is lower than expected for diffuse astrophysical plasmas, which have a very low resistivity (high magnetic Prandtl number); we attribute this to the large numerical resistivity in simulations, corresponding to a magnetic Prandtl number of order unity. When the magnetic field begins to be dynamically significant in the core of the minihalo at z = 27, we map it on to a uniform grid and follow the evolution in an adaptive mesh refinement, MHD simulation in orion2. The non-linear evolution of the field in the orion2 simulation violates flux-freezing and is consistent with the theory proposed by Xu & Lazarian. The fields approach equipartition with kinetic energy at densities ∼1010–1012 cm−3. When the same calculation is carried out in orion2 with no magnetic fields, several protostars form, ranging in mass from ∼1 to 30 M⊙; with magnetic fields, only a single ∼30 M⊙ protostar forms by the end of the simulation. Magnetic fields thus suppress the formation of low-mass Pop III stars, yielding a top-heavy Pop III initial mass function and contributing to the absence of observed Pop III stars.
Stars rarely form in isolation. Nearly half of the stars in the Milky Way have a companion, and this fraction increases in star-forming regions. However, why some dense cores and filaments form bound ...pairs while others form single stars remains unclear. We present a set of three-dimensional, gravo-magnetohydrodynamic simulations of turbulent star-forming clouds, aimed at understanding the formation and evolution of multiple-star systems formed through large-scale ( 103 au) turbulent fragmentation. We investigate three global magnetic field strengths, with global mass-to-flux ratios of φ = 2, 8, and 32. The initial separations of protostars in multiples depend on the global magnetic field strength, with stronger magnetic fields (e.g., φ = 2) suppressing fragmentation on smaller scales. The overall multiplicity fraction (MF) is between 0.4 and 0.6 for our strong and intermediate magnetic field strengths, which is in agreement with observations. The weak field case has a lower fraction. The MF is relatively constant throughout the simulations, even though stellar densities increase as collapse continues. While the MF rarely exceeds 60% in all three simulations, over 80% of all protostars are part of a binary system at some point. We additionally find that the distribution of binary spin misalignment angles is consistent with a randomized distribution. In all three simulations, several binaries originate with wide separations and dynamically evolve to 102 au separations. We show that a simple model of mass accretion and dynamical friction with the gas can explain this orbital evolution.
ABSTRACT
The ATOMS, standing for ALMA Three-millimeter Observations of Massive Star-forming regions, survey has observed 146 active star-forming regions with ALMA band 3, aiming to systematically ...investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the ATOMS survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS J = 2−1, HCO+J = 1−0, and HCN J = 1−0, are found to show extended gas emission in low-density regions within the clump; less than 25 per cent of their emission is from dense cores. SO, CH3OH, H13CN, and HC3N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over ∼1 pc), which may be caused by slow shocks from large–scale colliding flows or H ii regions. Stellar feedback from an expanding H ii region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analysed with other survey data, e.g. MALT90, Orion B, EMPIRE, ALMA_IMF, and ALMAGAL, to deepen our understandings of ‘dense gas’ star formation scaling relations and massive protocluster formation.
ABSTRACT
We present new 3-mm continuum and molecular lines observations from the ATOMS survey towards the massive protostellar clump, MM1, located in the filamentary infrared dark cloud (IRDC), ...G034.43+00.24 (G34). The lines observed are the tracers of either dense gas (e.g. HCO+/H13CO+ J= 1–0) or outflows (e.g. CS J= 2–1). The most complete picture to date of seven cores in MM1 is revealed by dust continuum emission. These cores are found to be gravitationally bound, with virial parameter, αvir < 2. At least four outflows are identified in MM1 with a total outflowing mass of ∼45 M⊙, and a total energy of 1 × 1047 erg, typical of outflows from a B0-type star. Evidence of hierarchical fragmentation, where turbulence dominates over thermal pressure, is observed at both the cloud and the clump scales. This could be linked to the scale-dependent, dynamical mass inflow/accretion on clump and core scales. We therefore suggest that the G34 cloud could be undergoing a dynamical mass inflow/accretion process linked to the multiscale fragmentation, which leads to the sequential formation of fragments of the initial cloud, clumps, and ultimately dense cores, the sites of star formation.
The Davis–Chandrasekhar–Fermi method revisited Chen, Che-Yu; Li, Zhi-Yun; Mazzei, Renato R ...
Monthly notices of the Royal Astronomical Society,
06/2022, Letnik:
514, Številka:
2
Journal Article
Recenzirano
Odprti dostop
ABSTRACT
Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical ...scales, because direct measurement of the field strength is observationally difficult. The Davis–Chandrasekhar–Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength from polarization data. It is based on the assumption that gas turbulent motion is the driving source of field distortion via linear Alfvén waves. In this work, using MHD simulations of star-forming clouds, we test the validity of the assumption underlying the DCF method by examining its accuracy in the real 3D space. Our results suggest that the DCF relation between turbulent kinetic energy and magnetic energy fluctuation should be treated as a statistical result instead of a local property. We then develop and investigate several modifications to the original DCF method using synthetic observations, and propose new recipes to improve the accuracy of DCF-derived magnetic field strength. We further note that the biggest uncertainty in the DCF analysis may come from the linewidth measurement instead of the polarization observation, especially since the line-of-sight gas velocity can be used to estimate the gas volume density, another critical parameter in the DCF method.