An official website of the United States government
ABSTRACT Stars form from the gravitational collapse of turbulent, magnetized molecular cloud cores. Our non-ideal MHD simulations reveal that the intrinsically anisotropic magnetic resistance to gravity during the core collapse naturally generates dense gravomagneto sheetlets within inner protostellar envelopes – disrupted versions of classical sheet-like pseudo-discs. They are embedded in a magnetically dominant background, where less dense materials flow along the local magnetic field lines and accumulate in the dense sheetlets. The sheetlets, which feed the disc predominantly through its upper and lower surfaces, are the primary channels for mass and angular momentum transfer from the envelope to the disc. The protostellar disc inherits a small fraction (up to 10 per cent) of the magnetic flux from the envelope, resulting in a disc-averaged net vertical field strength of 1–10 mG and a somewhat stronger toroidal field, potentially detectable through ALMA Zeeman observations. The inherited magnetic field from the envelope plays a dominant role in disc angular momentum evolution, enabling the formation of gravitationally stable discs in cases where the disc field is relatively well-coupled to the gas. Its influence remains significant even in marginally gravitationally unstable discs formed in the more magnetically diffusive cases, removing angular momentum at a rate comparable to or greater than that caused by spiral arms. The magnetically driven disc evolution is consistent with the apparent scarcity of prominent spirals capable of driving rapid accretion in deeply embedded protostellar discs. The dense gravomagneto sheetlets observed in our simulations may correspond to the ‘accretion streamers’ increasingly detected around protostars.
more »
« less
Fragmentation of dense rotation-dominated structures fed by collapsing gravo-magneto-sheetlets and origin of misaligned 100 au-scale binaries and multiple systems
ABSTRACT The majority of stars are in binary/multiple systems. How such systems form in turbulent, magnetized cores of molecular clouds in the presence of non-ideal magnetohydrodynamic (MHD) effects remains relatively underexplored. Through athena++-based non-ideal MHD adaptive mesh refinement simulations with ambipolar diffusion, we show that the collapsing protostellar envelope is dominated by dense gravo-magneto-sheetlets, a turbulence-warped version of the classic pseudodisc produced by anisotropic magnetic resistance to the gravitational collapse, in agreement with previous simulations of turbulent, magnetized single-star formation. The sheetlets feed mass, magnetic fields, and angular momentum to a Dense ROtation-Dominated (DROD) structure, which fragments into binary/multiple systems. This DROD fragmentation scenario is a more dynamic variant of the traditional disc fragmentation scenario for binary/multiple formation, with dense spiral filaments created by inhomogeneous feeding from the highly structured larger-scale sheetlets rather than the need for angular momentum transport, which is dominated by magnetic braking. Provided that the local material is sufficiently demagnetized, with a plasma-$$\beta$$ of 10 or more, collisions between the dense spiralling filaments play a key role in facilitating gravitational collapse and stellar companion formation by pushing the local magnetic Toomre parameter $$Q_\mathrm{m}$$ below unity. This mechanism can naturally produce in situ misaligned systems on the 100-au scale, often detected with high-resolution Atacama Large Millimeter Array (ALMA) observations. Our simulations also highlight the importance of non-ideal MHD effects, which affect whether fragmentation occurs and, if so, the masses and orbital parameters of the stellar companions formed.
more »
« less
- Award ID(s):
- 1910106
- PAR ID:
- 10525945
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 532
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 3135-3150
- Size(s):
- p. 3135-3150
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)ABSTRACT Non-ideal magnetohydrodynamic (MHD) effects have been shown recently as a robust mechanism of averting the magnetic braking ‘catastrophe’ and promoting protostellar disc formation. However, the magnetic diffusivities that determine the efficiency of non-ideal MHD effects are highly sensitive to microphysics. We carry out non-ideal MHD simulations to explore the role of microphysics on disc formation and the interplay between ambipolar diffusion (AD) and Hall effect during the protostellar collapse. We find that removing the smallest grain population (≲10 nm) from the standard MRN size distribution is sufficient for enabling disc formation. Further varying the grain sizes can result in either a Hall-dominated or an AD-dominated collapse; both form discs of tens of au in size regardless of the magnetic field polarity. The direction of disc rotation is bimodal in the Hall-dominated collapse but unimodal in the AD-dominated collapse. We also find that AD and Hall effect can operate either with or against each other in both radial and azimuthal directions, yet the combined effect of AD and Hall is to move the magnetic field radially outward relative to the infalling envelope matter. In addition, microphysics and magnetic field polarity can leave profound imprints both on observables (e.g. outflow morphology, disc to stellar mass ratio) and on the magnetic field characteristics of protoplanetary discs. Including Hall effect relaxes the requirements on microphysics for disc formation, so that prestellar cores with cosmic ray ionization rate of ≲2–3 × 10−16 s−1 can still form small discs of ≲10 au radius. We conclude that disc formation should be relatively common for typical prestellar core conditions, and that microphysics in the protostellar envelope is essential to not only disc formation, but also protoplanetary disc evolution.more » « less
-
ABSTRACT We study the flow structure in 3D magnetohydrodynamic (MHD) simulations of accretion on to Sagittarius A* via the magnetized winds of the orbiting Wolf–Rayet stars. These simulations cover over 3 orders of magnitude in radius to reach ≈300 gravitational radii, with only one poorly constrained parameter (the magnetic field in the stellar winds). Even for winds with relatively weak magnetic fields (e.g. plasma β ∼ 106), flux freezing/compression in the inflowing gas amplifies the field to β ∼ few well before it reaches the event horizon. Overall, the dynamics, accretion rate, and spherically averaged flow profiles (e.g. density, velocity) in our MHD simulations are remarkably similar to analogous hydrodynamic simulations. We attribute this to the broad distribution of angular momentum provided by the stellar winds, which sources accretion even absent much angular momentum transport. We find that the magneto-rotational instability is not important because of (i) strong magnetic fields that are amplified by flux freezing/compression, and (ii) the rapid inflow/outflow times of the gas and inefficient radiative cooling preclude circularization. The primary effect of magnetic fields is that they drive a polar outflow that is absent in hydrodynamics. The dynamical state of the accretion flow found in our simulations is unlike the rotationally supported tori used as initial conditions in horizon scale simulations, which could have implications for models being used to interpret Event Horizon Telescope and GRAVITY observations of Sgr A*.more » « less
-
null (Ed.)ABSTRACT Understanding the evolution of self-gravitating, isothermal, magnetized gas is crucial for star formation, as these physical processes have been postulated to set the initial mass function (IMF). We present a suite of isothermal magnetohydrodynamic (MHD) simulations using the gizmo code that follow the formation of individual stars in giant molecular clouds (GMCs), spanning a range of Mach numbers found in observed GMCs ($$\mathcal {M} \sim 10\!-\!50$$). As in past works, the mean and median stellar masses are sensitive to numerical resolution, because they are sensitive to low-mass stars that contribute a vanishing fraction of the overall stellar mass. The mass-weighted median stellar mass M50 becomes insensitive to resolution once turbulent fragmentation is well resolved. Without imposing Larson-like scaling laws, our simulations find $$M_\mathrm{50} \,\, \buildrel\propto \over \sim \,\,M_\mathrm{0} \mathcal {M}^{-3} \alpha _\mathrm{turb}\, \mathrm{SFE}^{1/3}$$ for GMC mass M0, sonic Mach number $$\mathcal {M}$$, virial parameter αturb, and star formation efficiency SFE = M⋆/M0. This fit agrees well with previous IMF results from the ramses, orion2, and sphng codes. Although M50 has no significant dependence on the magnetic field strength at the cloud scale, MHD is necessary to prevent a fragmentation cascade that results in non-convergent stellar masses. For initial conditions and SFE similar to star-forming GMCs in our Galaxy, we predict M50 to be $$\gt 20 \, \mathrm{M}_{\odot }$$, an order of magnitude larger than observed ($$\sim 2 \, \mathrm{M}_\odot$$), together with an excess of brown dwarfs. Moreover, M50 is sensitive to initial cloud properties and evolves strongly in time within a given cloud, predicting much larger IMF variations than are observationally allowed. We conclude that physics beyond MHD turbulence and gravity are necessary ingredients for the IMF.more » « less
-
Abstract Gravitational-wave observations provide the unique opportunity of studying black hole formation channels and histories—but only if we can identify their origin. One such formation mechanism is the dynamical synthesis of black hole binaries in dense stellar systems. Given the expected isotropic distribution of component spins of binary black holes in gas-free dynamical environments, the presence of antialigned or in-plane spins with respect to the orbital angular momentum is considered a tell-tale sign of a merger’s dynamical origin. Even in the scenario where birth spins of black holes are low, hierarchical mergers attain large component spins due to the orbital angular momentum of the prior merger. However, measuring such spin configurations is difficult. Here, we quantify the efficacy of the spin parameters encoding aligned-spin (χeff) and in-plane spin (χp) at classifying such hierarchical systems. Using Monte Carlo cluster simulations to generate a realistic distribution of hierarchical merger parameters from globular clusters, we can infer mergers’χeffandχp. The cluster populations are simulated using Advanced LIGO-Virgo sensitivity during the detector network’s third observing period and projections for design sensitivity. Using a “likelihood-ratio”-based statistic, we find that ∼2% of the recovered population by the current gravitational-wave detector network has a statistically significantχpmeasurement, whereas noχeffmeasurement was capable of confidently determining a system to be antialigned with the orbital angular momentum at current detector sensitivities. These results indicate that measuring spin-precession throughχpis a more detectable signature of hierarchical mergers and dynamical formation than antialigned spins.more » « less
