BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Pak Shing Li

Pak Shing Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.COastro-phastro-ph.IMphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

15works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): How do dense core properties affect the multiplicity of protostars?

During the transition phase from a prestellar to a protostellar cloud core, one or several protostars can form within a single gas core. The detailed physical processes of this transition, however, still remain unclear. We present 1.3 mm dust continuum and molecular line observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward 43 protostellar cores in the Orion Molecular Cloud Complex ($λ$ Orionis, Orion B, and Orion A) with an angular resolution of $\sim$ 0.35" ($\sim$ 140 au). In total, we detect 13 binary/multiple systems. We derive an overall multiplicity frequency (MF) of 28$\%$ $\pm$ 4$\%$ and a companion star fraction (CSF) of 51$\%$ $\pm$ 6$\%$, over a separation range of 300-8900 au. The median separation of companions is about 2100 au. The occurrence of stellar multiplicity may depend on the physical characteristics of the dense cores. Notably, those containing binary/multiple systems tend to show higher gas density and Mach number than cores forming single stars. The integral-shaped filament (ISF) of Orion A giant molecular cloud (GMC), which has the highest gas density and hosts high-mass star formation in its central region (the Orion Nebula cluster), shows the highest MF and CSF among the Orion GMCs. In contrast, the $λ$ Orionis Giant Molecular Cloud (GMC) has a lower MF and CSF than the Orion B and Orion A GMCs, indicating that feedback from HII regions may suppress the formation of multiple systems. We also find that the protostars comprising a binary/multiple system are usually at different evolutionary stages.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Magnetic fields in the formation of the first stars.--II Results

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 onto a uniform grid and follow the evolution in an adaptive mesh refinement, MHD simulation in Orion2. The nonlinear 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 ~ 10^10 - 10^12 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_sol with magnetic fields, only a single ~ 30 M_sol 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 IMF and contributing to the absence of observed Pop III stars.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Davis-Chandrasekhar-Fermi Method Revisited

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Role of Magnetic Fields in the Stability and Fragmentation of Filamentary Molecular Clouds: Two Case Studies at OMC-3 and OMC-4

We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained $214~μ$m HAWC+/SOFIA polarization data, and we make use of archival data for the column density and C$^{18}$O (1-0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anti-correlated with the column density and the polarization-angle dispersion function. We find that the filamentary cloud and dense clumps in OMC-3 are magnetically supercritical and strongly subvirial. This region should be in the gravitational collapse phase and is consistent with many young stellar objects (YSOs) forming in the region. Our histogram of relative orientations (HROs) analysis shows that the magnetic field is dynamically sub-dominant in the dense gas structures of OMC-3. We present the first polarization map of OMC-4. We find that the observed region is generally magnetically subcritical except for an elongated dense core, which could be a result of projection effect of a filamentary structure aligned close to the line-of-sight. The relative large velocity dispersion and the unusual positive shape parameters at high column densities in the HROs analysis suggest that our viewing angle may be close to axes of filamentary substructures in OMC-4. The dominating strong magnetic field in OMC-4 is unfavorable for star formation and is consistent with much fewer YSOs than in OMC-3.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic Fields in the Formation of the First Stars. I. Theory vs. Simulation

While magnetic fields are important in contemporary star formation, their role in primordial star formation is unknown. Magnetic fields of order 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 minihalos 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. Analytic 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 2-3 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 (AMR) 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 analytic 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Massive Star Formation via the Collapse of Subvirial and Virialized Turbulent Massive Cores

Similar to their low-mass counterparts, massive stars likely form via the collapse of pre-stellar 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 pre-stellar 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Formation and Evolution of Wide-Orbit Stellar Multiples In Magnetized Clouds

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 (>~$10^3$ 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 depends 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-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 mis-alignment angles is consistent with a randomized distribution. In all three simulations, several binaries originate with wide separations and dynamically evolve to <~ $10^2$ AU separations. We show that a simple model of mass accretion and dynamical friction with the gas can explain this orbital evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetized Interstellar Molecular Clouds. I. Comparison Between Simulations and Zeeman Observations

The most accurate measurements of magnetic fields in star-forming gas are based on the Zeeman observations analyzed by Crutcher et al. (2010). We show that their finding that the 3D magnetic field scales approximately as density$^{0.65}$ can also be obtained from analysis of the observed line-of-sight fields. We present two large-scale AMR MHD simulations of several thousand $M_\odot$ of turbulent, isothermal, self-gravitating gas, one with a strong initial magnetic field (Alfven Mach number $M_{A,0}= 1$) and one with a weak initial field ($M_{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 in excellent agreement with the observed one. Our results imply that the average field in molecular clumps in the interstellar medium is $<B_{tot}> \sim 42 n_{H,4}^{0.65} μ$G. Furthermore, the median value of the ratio of the line-of-sight field to density$^{0.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% of the line-of-sight value. This is larger than the 50% 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. (2009), 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The CH+ Abundance in Turbulent, Diffuse Molecular Clouds

The intermittent dissipation of interstellar turbulence is an important energy source in the diffuse ISM. Though on average smaller than the heating rates due to cosmic rays and the photoelectric effect on dust grains, the turbulent cascade can channel large amounts of energy into a relatively small fraction of the gas that consequently undergoes significant heating and chemical enrichment. In particular, this mechanism has been proposed as a solution to the long-standing problem of the high abundance of CH+ along diffuse molecular sight lines, which steady-state, low temperature models under-produce by over an order of magnitude. While much work has been done on the structure and chemistry of these small-scale dissipation zones, comparatively little attention has been paid to relating these zones to the properties of the large-scale turbulence. In this paper, we attempt to bridge this gap by estimating the temperature and CH+ column density along diffuse molecular sight-lines by post-processing 3-dimensional MHD turbulence simulations. Assuming reasonable values for the cloud density (30 / cm^3), size (20 pc), and velocity dispersion (2.3 km / s), we find that our computed abundances compare well with CH+ column density observations, as well as with observations of emission lines from rotationally excited H2 molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ambipolar Diffusion Heating in Turbulent Systems

The temperature of the gas in molecular clouds is a key determinant of the characteristic mass of star formation. Ambipolar diffusion (AD) is considered one of the most important heating mechanisms in weakly ionized molecular clouds. In this work, we study the AD heating rate using 2-fluid turbulence simulations and compare it with the overall heating rate due to turbulent dissipation. We find that for observed molecular clouds, which typically have Alfven Mach numbers of ~1 (Crutcher 1999) and AD Reynolds numbers of ~20 (McKee et al. 2010), about 70% of the total turbulent dissipation is in the form of AD heating. AD has an important effect on the length scale where energy is dissipated: when AD heating is strong, most of the energy in the cascade is removed by ion-neutral drift, with a comparatively small amount of energy making it down to small scales. We derive a relation for the AD heating rate that describes the results of our simulations to within a factor of two. Turbulent dissipation, including AD heating, is generally less important that cosmic-ray heating in molecular clouds, although there is substantial scatter in both.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Stable, Accurate Methodology for High Mach Number, Strong Magnetic Field MHD Turbulence with Adaptive Mesh Refinement: Resolution and Refinement Studies

Performing a stable, long duration simulation of driven MHD turbulence with a high thermal Mach number and a strong initial magnetic field is a challenge to high-order Godunov ideal MHD schemes because of the difficulty in guaranteeing positivity of the density and pressure. We have implemented a robust combination of reconstruction schemes, Riemann solvers, limiters, and Constrained Transport EMF averaging schemes that can meet this challenge, and using this strategy, we have developed a new Adaptive Mesh Refinement (AMR) MHD module of the ORION2 code. We investigate the effects of AMR on several statistical properties of a turbulent ideal MHD system with a thermal Mach number of 10 and a plasma $β_0$ of 0.1 as initial conditions; our code is shown to be stable for simulations with higher Mach numbers ($M_rms = 17.3$) and smaller plasma beta ($β_0 = 0.0067$) as well. Our results show that the quality of the turbulence simulation is generally related to the volume-averaged refinement. Our AMR simulations show that the turbulent dissipation coefficient for supersonic MHD turbulence is about 0.5, in agreement with unigrid simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sub-Alfvenic Non-Ideal MHD Turbulence Simulations with Ambipolar Diffusion: III. Implications for Observations and Turbulent Enhancement

Ambipolar diffusion (AD) is believed to be a crucial process for redistributing magnetic flux in the dense molecular gas that occurs in regions of star formation. We carry out numerical simulations of this process in regions of low ionization using the heavy ion approximation. The simulations are for regions of strong field (plasma β=0.1) and mildly supersonic turbulence (M=3, corresponding to an Alfven mach number of 0.67). The velocity power spectrum of the neutral gas changes from an Iroshnikov-Kraichnan spectrum in the case of ideal MHD to a Burgers spectrum in the case of a shock-dominated hydrodynamic system. The magnetic power spectrum shows a similar behavior. We use a 1D radiative transfer code to post-process our simulation results; the simulated emission from the CS J=2-1 and H13CO+ J=1-0 lines shows that the effects of AD are observable in principle. Linewidths of ions are observed to be less than those of neutrals, and we confirm previous suggestions that this is due to AD. We show that AD is unlikely to affect the Chandrasekhar-Fermi method for inferring field strengths unless the AD is stronger than generally observed. Finally, we present the first fully 3D study of the enhancement of AD by turbulence, finding that AD is accelerated by factor 2-4.5 for non self-gravitating systems with the level of turbulence we consider.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Two regimes of Turbulent Fragmentation and the stellar IMF from Primordial to Present Day Star Formation

The Padoan and Nordlund model of the stellar initial mass function (IMF) is derived from low order statistics of supersonic turbulence, neglecting gravity (e.g. gravitational fragmentation, accretion and merging). In this work the predictions of that model are tested using the largest numerical experiments of supersonic hydrodynamic (HD) and magneto-hydrodynamic (MHD) turbulence to date (~1000^3 computational zones) and three different codes (Enzo, Zeus and the Stagger Code). The model predicts a power law distribution for large masses, related to the turbulence energy power spectrum slope, and the shock jump conditions. This power law mass distribution is confirmed by the numerical experiments. The model also predicts a sharp difference between the HD and MHD regimes, which is recovered in the experiments as well, implying that the magnetic field, even below energy equipartition on the large scale, is a crucial component of the process of turbulent fragmentation. These results suggest that the stellar IMF of primordial stars may differ from that in later epochs of star formation, due to differences in both gas temperature and magnetic field strength. In particular, we find that the IMF of primordial stars born in turbulent clouds may be narrowly peaked around a mass of order 10 solar masses, as long as the column density of such clouds is not much in excess of 10^22 cm^-2.

0 citations0 reviewsNo code linkedOpen access