Researcher profile

Pak Shing Li

Pak Shing Li contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 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.

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.

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.

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.

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.

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.

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.

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.