Researcher profile

Xue-Ning Bai

Xue-Ning Bai contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
6topics
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

10 published item(s)

preprint2022arXiv

A Multi-Fluid Dust Module in Athena++: Algorithms and Numerical Tests

We describe the algorithm, implementation and numerical tests of a multifluid dust module in the Athena++ magnetohydrodynamic (MHD) code. The module can accommodate an arbitrary number of dust species interacting with the gas via aerodynamic drag (characterized by the stopping time), with a number of numerical solvers. In particular, we describe two second-order accurate, two-stage, fully-implicit solvers that are stable in stiff regimes including short stopping time and high dust mass loading, and they are paired with the second-order explicit van-Leer and Runge-Kutta gas dynamics solvers in Athena++, respectively. Moreover, we formulate a consistent treatment of dust concentration diffusion with dust back-reaction, which incorporates momentum diffusion and ensures Galilean invariance. The new formulation and stiff drag solvers are implemented to be compatible with most existing features of Athena++, including different coordinate systems, mesh refinement, shearing-box and orbital advection. We present a large suite of test problems, including the streaming instability in linear and nonlinear regimes, as well as local and global setting, which demonstrate that the code achieves the desired performance. This module will be particularly useful for studies of dust dynamics and planet formation in protoplanetary disks.

preprint2022arXiv

Planetary Accretion Shocks with a Realistic Equation of State

The final stage of gas giant formation involves accreting gas from the parent protoplanetary disk. In general, the infalling gas likely approaches a free-fall velocity, creating an accretion shock, leading to strong shock heating and radiation. We investigate the kinematics and energetics of such accretion shocks using 1D radiation hydrodynamic simulations. Our simulations feature the first self-consistent treatment of hydrogen dissociation and ionization, radiation transport, and realistic grey opacity. By exploring a broad range of giant planet masses (0.1-3 M$_{J}$) and accretion rates ($10^{-3}$-$10^{-2}$M$_{\oplus}\cdot\rm{yr}^{-1}$), we focus on global shock efficiency and the final entropy of the accreted gas. We find that radiation from the accretion shock can fully disassociate the molecular hydrogen of the incoming gas when the shock luminosity is above a critical luminosity. Meanwhile, the post-shock entropy generally fall into &#34;cold&#34; ($<12k_{\rm{B}}/m_{\rm H}$) and &#34;hot&#34; ($>16k_{\rm{B}}/m_{\rm H}$) groups which depends on the extent of the endothermic process of $\rm{H}_2$ dissociation. While 2D or 3D simulations are needed for more realistic understandings of the accretion process, this distinction likely carries over and sheds light on the interpretation of young direct imaging planets.

preprint2022arXiv

Towards First-principle Characterization of Cosmic-ray Transport Coefficients from Multi-scale Kinetic Simulations

A major uncertainty in understanding the transport and feedback of cosmic-rays (CRs) within and beyond our Galaxy lies in the unknown CR scattering rates, which are primarily determined by wave-particle interaction at microscopic gyro-resonant scales. The source of the waves for the bulk CR population is believed to be self-driven by the CR streaming instability (CRSI), resulting from the streaming of CRs downward a CR pressure gradient. While a balance between driving by the CRSI and wave damping is expected to determine wave amplitudes and hence the CR scattering rates, the problem involves significant scale separation with substantial ambiguities based on quasi-linear theory (QLT). Here we propose a novel &#34;streaming box&#34; framework to study the CRSI with an imposed CR pressure gradient, enabling first-principle measurement of the CR scattering rates as a function of environmental parameters. By employing the magnetohydrodynamic-particle-in-cell (MHD-PIC) method with ion-neutral damping, we conduct a series of simulations with different resolutions and CR pressure gradients and precisely measure the resulting CR scattering rates in steady state. The measured rates show scalings consistent with QLT, but with a normalization smaller by a factor of several than typical estimates based on single-fluid treatment of CRs. A momentum-by-momentum treatment provides better estimates when integrated over momentum, but is also subject substantial deviations especially at small momentum. Our framework thus opens up the path towards providing comprehensive subgrid physics for macroscopic studies of CR transport and feedback in broad astrophysical contexts.

preprint2021arXiv

Cosmic Ray Transport in the Ionized and Neutral ISM: MHD-PIC Simulations and Effective Fluid Treatments

Cosmic rays (CRs) have critical impacts in the multiphase interstellar medium (ISM), driving dynamical motions in low-density plasma and modifying the ionization state, temperature, and chemical composition of higher-density atomic and molecular gas. We present a study of CR propagation between the ionized ISM and a neutral cloud. Using one-dimensional magnetohydrodynamic particle-in-cell simulations which include ion-neutral drag to damp Alfv$\acute{\text{e}}$n waves in the cloud, we self-consistently evolve the kinetic physics of CRs and fluid dynamics of the multiphase gas. By introducing the cloud in our periodic domain, our simulations break translational symmetry and allow the emergence of spatial structure in the CR distribution function. A negative spatial gradient forms across the fully-ionized ISM region while a positive gradient forms across the neutral cloud. We connect our results with CR hydrodynamics formulations by computing the wave-particle scattering rates as predicted by quasilinear, fluid, and Fokker-Planck theory. For momenta where the mean free path is short relative to the box size, we find excellent agreement among all scattering rates. By exploring different cloud sizes and ion-neutral collision rates, we show that our results are robust. Our work provides a first-principles verification of CR hydrodynamics when particles stream down their pressure gradient, and opens a pathway toward comprehensive calibrations of transport coefficients from self-generated Alfv$\acute{\text{e}}$n wave scattering with CRs.

preprint2021arXiv

Dust Settling and Clumping in MRI Turbulent Outer Protoplanetary Disks

Planetesimal formation is a crucial yet poorly understood process in planet formation. It is widely believed that planetesimal formation is the outcome of dust clumping by the streaming instability (SI). However, recent analytical and numerical studies have shown that the SI can be damped or suppressed by external turbulence, and at least the outer regions of protoplanetary disks are likely weakly turbulent due to magneto-rotational instability (MRI). We conduct high-resolution local shearing-box simulations of hybrid particle-gas magnetohydrodynamics (MHD), incorporating ambipolar diffusion as the dominant non-ideal MHD effect, applicable to outer disk regions. We first show that dust backreaction enhances dust settling towards the midplane by reducing turbulence correlation time. Under modest level of MRI turbulence, we find that dust clumping is in fact easier than the conventional SI case, in the sense that the threshold of solid abundance for clumping is lower. The key to dust clumping includes dust backreaction and the presence of local pressure maxima, which in our work is formed by the MRI zonal flows overcoming background pressure gradient. Overall, our results support planetesimal formation in the MRI-turbulent outer protoplanetary disks, especially in ring-like substructures.

preprint2021arXiv

Dust Transport in Protoplanetary Disks with Wind-driven Accretion

It has recently been shown that the inner region of protoplanetary disks (PPDs) is governed by wind-driven accretion, and the resulting accretion flow showing complex vertical profiles. Such complex flow structures are further enhanced due to the Hall effect, especially when the background magnetic field is aligned with disk rotation. We investigate how such flow structures impact global dust transport via Monte-Carlo simulations, focusing on two scenarios. In the first scenario, the toroidal magnetic field is maximized in the miplane, leading to accretion and decretion flows above and below. In the second scenario, the toroidal field changes sign across the midplane, leading to an accretion flow at the disk midplane, with decretion flows above and below. We find that in both cases, the contribution from additional gas flows can still be accurately incorporated into the advection-diffusion framework for vertically-integrated dust transport, with enhanced dust radial diffusion up to an effective $α^{\rm eff}\sim10^{-2}$ for strongly coupled dust, even when background turbulence is weak $α<10^{-4}$. Dust radial drift is also modestly enhanced in the second scenario. We provide a general analytical theory that accurately reproduces our simulation results, thus establishing a framework to model global dust transport that realistically incorporates vertical gas flow structures. We also note that the theory is equally applicable to the transport of chemical species.

preprint2021arXiv

History of the Solar Nebula from Meteorite Paleomagnetism

We review recent advances in our understanding of magnetism in the solar nebular and protoplanetary disks (PPDs). We discuss the implications of theory, meteorite measurements, and astronomical observations for planetary formation and nebular evolution. Paleomagnetic measurements indicate the presence of fields of 0.54$\pm$0.21 G at $\sim$1 to 3 astronomical units (AU) from the Sun and $\gtrsim$0.06 G at 3 to 7 AU until >1.22 and >2.51 million years (Ma) after solar system formation, respectively. These intensities are consistent with those predicted to enable typical astronomically-observed protostellar accretion rates of $\sim$10$^{-8}$ M$_\odot$ yr$^{-1}$, suggesting that magnetism played a central role in mass and angular momentum transport in PPDs. Paleomagnetic studies also indicate fields <0.006 G and <0.003 G in the inner and outer solar system by 3.94 and 4.89 Ma, respectively, consistent with the nebular gas having dispersed by this time. This is similar to the observed lifetimes of extrasolar protoplanetary disks.

preprint2020arXiv

Global Simulations of the Vertical Shear Instability with Non-ideal Magnetohydrodynamical Effects

The mechanisms of angular momentum transport and level of turbulence in protoplanetary disks (PPDs) are crucial for understanding many aspects of planet formation. In the recent years, it has been realized that the magneto-rotational instability (MRI) tends to be suppressed in PPDs due to non-ideal MHD effects, and the disk is largely laminar with accretion driven by magnetized disk winds. In parallel, several hydrodynamical mechanisms have been identified that likely also generate vigorous turbulence and drive disk accretion. We study the interplay between MHD winds in PPDs with the vertical shear instability (VSI), one of the most promising hydrodynamical mechanisms, through 2D global non-ideal MHD simulations with ambipolar diffusion and Ohmic resistivity. We find that for typical disk parameters, MHD winds can coexist with the VSI with accretion primarily wind-driven accompanied by vigorous VSI turbulence. The properties of the VSI remain similar to unmagnetized case, and the wind and overall field configuration are not strongly affected by VSI turbulence, showing modest level of variability and corrugation of midplane current sheet. Enhanced coupling between gas and magnetic field weakens the VSI. The VSI is also weakened with increasing magnetization, and we find that corrugation motions characteristic of the VSI transitions to low-amplitude breathing mode oscillations.

preprint2020arXiv

Measuring turbulent motion in planet-forming disks with ALMA: A detection around DM Tau and non-detections around MWC 480 and V4046 Sgr

Turbulence is a crucial factor in many models of planet formation, but it has only been directly constrained among a small number of planet forming disks. Building on the upper limits on turbulence placed in disks around HD 163296 and TW Hya, we present ALMA CO J=2-1 line observations at $\sim0.3&#34;$ (20-50 au) resolution and 80 ms$^{-1}$ channel spacing of the disks around DM Tau, MWC 480, and V4046 Sgr. Using parametric models of disk structure, we robustly detect non-thermal gas motions around DM Tau of between 0.25 c$_s$ and 0.33 c$_s$, with the range dominated by systematic effects, making this one of the only systems with directly measured non-zero turbulence. Using the same methodology, we place stringent upper limits on the non-thermal gas motion around MWC 480 ($<$0.08c$_s$) and V4046 Sgr ($<$0.12 c$_s$). The preponderance of upper limits in this small sample, and the modest turbulence levels consistent with dust studies, suggest that weak turbulence ($α\lesssim10^{-3}$) may be a common, albeit not universal, feature of planet-forming disks. We explore the particular physical conditions around DM Tau that could lead this system to be more turbulent than the others.

preprint2020arXiv

The Disk Substructures at High Angular Resolution Project (DSHARP): II. Characteristics of Annular Substructures

The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few au to tens of au. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures co-exist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and &#34;double gap&#34; features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk.