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

Xue-Ning Bai

Xue-Ning Bai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.HEastro-ph.GAastro-ph.IMphysics.plasm-phastro-ph.COphysics.comp-ph

Catalog footprint

What is connected

37works
8topics
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

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

0 citations0 reviewsNo code linkedOpen access
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 "cold" ($<12k_{\rm{B}}/m_{\rm H}$) and "hot" ($>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.

0 citations0 reviewsNo code linkedOpen access
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 "streaming box" 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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"$ (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.

0 citations0 reviewsNo code linkedOpen access
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 "double gap" features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

ALMA Observations of a Gap and a Ring in the Protoplanetary Disk around TW Hya

We report the first detection of a gap and a ring in 336 GHz dust continuum emission from the protoplanetary disk around TW Hya, using the Atacama Large Millimeter/Submillimeter Array (ALMA). The gap and ring are located at around 25 and 41 au from the central star, respectively, and are associated with the CO snow line at ~30 au. The gap has a radial width of less than 15 au and a mass deficit of more than 23%, taking into account that the observations are limited to an angular resolution of ~15 au. In addition, the 13CO and C18O J = 3 - 2 lines show a decrement in CO line emission throughout the disk, down to ~10 au, indicating a freeze-out of gas-phase CO onto grain surfaces and possible subsequent surface reactions to form larger molecules. The observed gap could be caused by gravitational interaction between the disk gas and a planet with a mass less than super-Neptune (2M_{Neptune}), or could be the result of the destruction of large dust aggregates due to the sintering of CO ice.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Ringed Substructure and a Gap at 1 AU in the Nearest Protoplanetary Disk

We present long-baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870 micron continuum emission from the nearest gas-rich protoplanetary disk, around TW Hya, that trace millimeter-sized particles down to spatial scales as small as 1 AU (20 mas). These data reveal a series of concentric ring-shaped substructures in the form of bright zones and narrow dark annuli (1-6 AU) with modest contrasts (5-30%). We associate these features with concentrations of solids that have had their inward radial drift slowed or stopped, presumably at local gas pressure maxima. No significant non-axisymmetric structures are detected. Some of the observed features occur near temperatures that may be associated with the condensation fronts of major volatile species, but the relatively small brightness contrasts may also be a consequence of magnetized disk evolution (the so-called zonal flows). Other features, particularly a narrow dark annulus located only 1 AU from the star, could indicate interactions between the disk and young planets. These data signal that ordered substructures on ~AU scales can be common, fundamental factors in disk evolution, and that high resolution microwave imaging can help characterize them during the epoch of planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Towards a Global Evolutionary Model of Protoplanetary Disks

A global evolution picture of protoplanetary disks (PPDs) is key to understanding almost every aspect of planet formation, where standard alpha-disk models have been constantly employed for its simplicity. In the mean time, disk mass loss has been conventionally attributed to photoevaporation, which controls disk dispersal. However, a paradigm shift towards accretion driven by magnetized disk winds has been realized in the recent years, thanks to studies of non-ideal magneto-hydrodynamic effects in PPDs. I present a framework of global PPD evolution aiming to incorporate these advances, highlighting the role of wind-driven accretion and wind mass loss. Disk evolution is found to be largely dominated by wind-driven processes, and viscous spreading is suppressed. The timescale of disk evolution is controlled primarily by the amount of external magnetic flux threading the disks, and how rapidly the disk loses the flux. Rapid disk dispersal can be achieved if the disk is able to hold most of its magnetic flux during the evolution. In addition, because wind launching requires sufficient level of ionization at disk surface (mainly via external far-UV radiation), wind kinematics is also affected by far-UV penetration depth and disk geometry. For typical disk lifetime of a few Myrs, the disk loses approximately the same amount of mass through the wind as through accretion onto the protostar, and most of the wind mass loss proceeds from the outer disk via a slow wind. Fractional wind mass loss increases with increasing disk lifetime. Significant wind mass loss likely substantially enhances the dust to gas mass ratio, and promotes planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magneto-thermal Disk Wind from Protoplanetary Disks

Global evolution and dispersal of protoplanetary disks (PPDs) is governed by disk angular momentum transport and mass-loss processes. Recent numerical studies suggest that angular momentum transport in the inner region of PPDs is largely driven by magnetized disk wind, yet the wind mass-loss rate remains unconstrained. On the other hand, disk mass loss has conventionally been attributed to photoevaporation, where external heating on the disk surface drives a thermal wind. We unify the two scenarios by developing a 1D model of magnetized disk winds with a simple treatment of thermodynamics as a proxy for external heating. The wind properties largely depend on 1) the magnetic field strength at the wind base, characterized by the poloidal Alfvén speed $v_{Ap}$, 2) the sound speed $c_s$ near the wind base, and 3) how rapidly poloidal field lines diverge (achieve $R^{-2}$ scaling). When $v_{Ap}\gg c_s$, corotation is enforced near the wind base, resulting in centrifugal acceleration. Otherwise, the wind is accelerated mainly by the pressure of the toroidal magnetic field. In both cases, the dominant role played by magnetic forces likely yields wind outflow rates that well exceed purely hydrodynamical mechanisms. For typical PPD accretion-rate and wind-launching conditions, we expect $v_{Ap}$ to be comparable to $c_s$ at the wind base. The resulting wind is heavily loaded, with total wind mass loss rate likely reaching a considerable fraction of wind-driven accretion rate. Implications for modeling global disk evolution and planet formation are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetohydrodynamic-Particle-in-Cell Method for Coupling Cosmic Rays with a Thermal Plasma: Application to Non-relativistic Shocks

We formulate a magnetohydrodynamic-particle-in-cell (MHD-PIC) method for describing the interaction between collisionless cosmic ray (CR) particles and a thermal plasma. The thermal plasma is treated as a fluid, obeying equations of ideal MHD, while CRs are treated as relativistic Lagrangian particles subject to the Lorentz force. Backreaction from CRs to the gas is included in the form of momentum and energy feedback. In addition, we include the electromagnetic feedback due to CR-induced Hall effect that becomes important when the electron-ion drift velocity of the background plasma induced by CRs approaches the Alfvén velocity. Our method is applicable on scales much larger than the ion inertial length, bypassing the microscopic scales that must be resolved in conventional PIC methods, while retaining the full kinetic nature of the CRs. We have implemented and tested this method in the Athena MHD code, where the overall scheme is second-order accurate and fully conservative. As a first application, we describe a numerical experiment to study particle acceleration in non-relativistic shocks. Using a simplified prescription for ion injection, we reproduce the shock structure and the CR energy spectra obtained with more self-consistent hybrid-PIC simulations, but at substantially reduced computational cost. We also show that the CR-induced Hall effect reduces the growth rate of the Bell instability and affects the gas dynamics in the vicinity of the shock front. As a step forward, we are able to capture the transition of particle acceleration from non relativistic to relativistic regimes, with momentum spectrum $f(p)\sim p^{-4}$ connecting smoothly through the transition, as expected from the theory of Fermi acceleration.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical Simulation of Hot Accretion Flows (III): Revisiting wind properties using trajectory approach

Previous MHD simulations have shown that wind must exist in black hole hot accretion flows. In this paper, we continue our study by investigating the detailed properties of wind, such as mass flux and poloidal speed, and the mechanism of wind production. For this aim, we make use of a three dimensional GRMHD simulation of hot accretion flows around a Schwarzschild black hole. The simulation is designed so that the magnetic flux is not accumulated significantly around the black hole. To distinguish real wind from turbulent outflows, we track the trajectories of the virtual Largrangian particles from simulation data. We find two types of real outflows, i.e., a quasi-relativistic jet close to the axis and a sub-relativistic wind subtending a much larger solid angle. Most of the wind originates from the surface layer of the accretion flow. The poloidal wind speed almost remains constant once they are produced, but the flux-weighted wind speed roughly follows $v_{\rm p, wind}(r)\approx 0.25 v_k(r)$. The mass flux of jet is much lower but the speed is much higher, $v_{\rm p,jet}\sim (0.3-0.4) c$. Consequently, both the energy and momentum fluxes of the wind are much larger than those of the jet. We find that the wind is produced and accelerated primarily by the combination of centrifugal force and magnetic pressure gradient, while the jet is mainly accelerated by magnetic pressure gradient. Finally, we find that the wind production efficiency $ε_{\rm wind}\equiv\dot{E}_{\rm wind}/\dot{M}_{\rm BH}c^2\sim 1/1000$, in good agreement with the value required from large-scale galaxy simulations with AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Grain-Modified Magnetic Diffusivities in Protoplanetary Disks

Weakly ionized protoplanetary disks (PPDs) are subject to non-ideal-magnetohydrodynamic (MHD) effects including Ohmic resistivity, the Hall effect and ambipolar diffusion (AD), and the resulting magnetic diffusivities ($η_O, η_H$ and $η_A$) largely control the disk gas dynamics. The presence of grains not only strongly reduces disk ionization fraction, but also modify the scalings of $η_H$ and $η_A$ with magnetic field strength. We derive analytically asymptotic expressions of $η_H$ and $η_A$ in both strong and weak field limits and show that towards strong field, $η_H$ can change sign (at a threshold field strength $B_{\rm th}$), mimicking a flip of field polarity, and AD is substantially reduced. Applying to PPDs, we find that when small $\sim0.1$ ($0.01$)$μ$m grains are sufficiently abundant [mass ratio $\sim0.01$ ($10^{-4}$)], $η_H$ can change sign up to $\sim2-3$ scale heights above midplane at modest field strength (plasma $β\sim100$) over a wide range of disk radii. Reduction of AD is also substantial towards the AD dominated outer disk and may activate the magneto-rotational instability. We further perform local non-ideal MHD simulations of the inner disk (within 10 AU) and show that with sufficiently abundant small grains, magnetic field amplification due to the Hall-shear instability saturates at very low level near the threshold field strength $B_{\rm th}$. Together with previous studies, we conclude by discussing the grain-abundance-dependent phenomenology of PPD gas dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Signatures of MRI-Driven Turbulence in Protoplanetary Disks: Predictions for ALMA Observations

Spatially resolved observations of molecular line emission have the potential to yield unique constraints on the nature of turbulence within protoplanetary disks. Using a combination of local non-ideal magnetohydrodynamic simulations and radiative transfer calculations, tailored to properties of the disk around HD 163296, we assess the ability of ALMA to detect turbulence driven by the magnetorotational instability (MRI). Our local simulations show that the MRI produces small-scale turbulent velocity fluctuations that increase in strength with height above the mid-plane. For a set of simulations at different disk radii, we fit a Maxell-Boltzmann distribution to the turbulent velocity and construct a turbulent broadening parameter as a function of radius and height. We input this broadening into radiative transfer calculations to quantify observational signatures of MRI-driven disk turbulence. We find that the ratio of the peak line flux to the flux at line center is a robust diagnostic of turbulence that is only mildly degenerate with systematic uncertainties in disk temperature. For the CO(3-2) line, which we expect to probe the most magnetically active slice of the disk column, variations in the predicted peak-to-trough ratio between our most and least turbulent models span a range of approximately 15%. Additional independent constraints can be derived from the morphology of spatially resolved line profiles, and we estimate the resolution required to detect turbulence on different spatial scales. We discuss the role of lower optical depth molecular tracers, which trace regions closer to the disk mid-plane where velocities in MRI-driven models are systematically lower.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dust Transport in MRI Turbulent Disks: Ideal and Non-ideal MHD with Ambipolar Diffusion

We study dust transport in turbulent protoplanetary disks using three-dimensional global unstratified magnetohydrodynamic (MHD) simulations including Lagrangian dust particles. The turbulence is driven by the magnetorotational instability (MRI) with either ideal or non-ideal MHD that includes ambipolar diffusion (AD). In ideal MHD simulations, the surface density evolution (except for dust that drifts fastest), turbulent diffusion, and vertical scale height of dust can all be reproduced by simple one-dimensoinal and/or analytical models. However, in AD dominated simulations which simulate protoplanetary disks beyond 10s of AU, the vertical scale height of dust is larger than previously predicted. To understand this anomaly in more detail, we carry out both unstratified and stratified local shearing box simulations with Lagrangian particles, and find that turbulence in AD dominated disks has very different properties (e.g., temporal autocorrelation functions and power spectra) than turbulence in ideal MHD disks, which leads to quite different particle diffusion efficiency. For example, MRI turbulence with AD has a longer correlation time for the vertical velocity, which causes significant vertical particle diffusion and large dust scale height. In ideal MHD the Schmidt numbers ($Sc$) for radial and vertical turbulent diffusion are $Sc_{r}\sim 1$ and $Sc_{z}\gtrsim 3$, but in the AD dominated regime both $Sc_{r}$ and $Sc_{z}$ are $\lesssim 1$. Particle concentration in pressure bumps induced by MRI turbulence has also been studied. Since non-ideal MHD effects dominate most regions in protoplanetary disks, our study suggests that modeling dust transport in turbulence driven by MRI with non-ideal MHD effects is important for understanding dust transport in realistic protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hall-effect Controlled Gas Dynamics in Protoplanetary Disks: I. Wind Solutions at the Inner Disk

The gas dynamics of protoplanetary disks (PPDs) is largely controlled by non-ideal magnetohydrodynamic (MHD) effects including Ohmic resistivity, the Hall effect and ambipolar diffusion. Among these the role of the Hall effect is the least explored and most poorly understood. We have included all three non-ideal MHD effects in a self-consistent manner to investigate the role of the Hall effect on PPD gas dynamics using local shearing-box simulations. In this first paper, we focus on the inner region of PPDs, where previous studies excluding the Hall effect have revealed that the inner disk up to ~10 AU is largely laminar, with accretion driven by a magnetocentrifugal wind. We confirm this basic picture and show that the Hall effect introduces modest modifications to the wind solutions, depending on the polarity of the large-scale poloidal magnetic field B_0 threading the disk. When B_0.Omega>0, the horizontal magnetic field is strongly amplified toward the disk interior, leading to a stronger disk wind (by ~50% or less in terms of the wind-driven accretion rate). The enhanced horizontal field also leads to much stronger large-scale Maxwell stress (magnetic braking) that contributes to a considerable fraction of the wind-driven accretion rate. When B_0.Omega<0, the horizontal magnetic field is reduced, leading to a weaker disk wind (by ~20%) and negligible magnetic braking. Moreover, we find that when B_0.Omega>0, the laminar region extends farther to ~15 AU before the magneto-rotational instability sets in, while for B_0.Omega<0, the laminar region extends only to ~3-5 AU for a typical PPD accretion rates. Scaling relations for the wind properties, especially the wind-driven accretion rate, are provided for aligned and anti-aligned field geometries. Issues with the symmetry of the wind solutions and grain abundance are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hall-effect Controlled Gas Dynamics in Protoplanetary Disks: II. Full 3D Simulations toward the Outer Disk

We perform 3D stratified shearing-box MHD simulations on the gas dynamics of protoplanetary disks threaded by net vertical magnetic field Bz. All three non-ideal MHD effects, Ohmic resistivity, the Hall effect and ambipolar diffusion are included in a self-consistent manner based on equilibrium chemistry. We focus on regions toward outer disk radii, from 5-60AU, where Ohmic resistivity tends to become negligible, ambipolar diffusion dominates over an extended region across disk height, and the Hall effect largely controls the dynamics near the disk midplane. We find that around R=5AU, the system launches a laminar/weakly turbulent magnetocentrifugal wind when the net vertical field Bz is not too weak, as expected. Moreover, the wind is able to achieve and maintain a configuration with reflection symmetry at disk midplane. The case with anti-aligned field polarity (Omega. Bz<0) is more susceptible to the MRI when Bz drops, leading to an outflow oscillating in radial directions and very inefficient angular momentum transport. At the outer disk around and beyond R=30AU, the system shows vigorous MRI turbulence in the surface layer due to far-UV ionization, which efficiently drives disk accretion. The Hall effect affects the stability of the midplane region to the MRI, leading to strong/weak Maxwell stress for aligned/anti-aligned field polarities. Nevertheless, the midplane region is only very weakly turbulent. Overall, the basic picture is analogous to the conventional layered accretion scenario applied to the outer disk. In addition, we find that the vertical magnetic flux is strongly concentrated into thin, azimuthally extended shells in most of our simulations beyond 15AU. This is a generic phenomenon unrelated to the Hall effect, and leads to enhanced zonal flow. Observational and theoretical implications, as well as future prospects are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetic Flux Concentration and Zonal Flows in Magnetorotational Instability Turbulence

Accretion disks are likely threaded by external vertical magnetic flux, which enhances the level of turbulence via the magnetorotational instability (MRI). Using shearing-box simulations, we find that such external magnetic flux also strongly enhances the amplitude of banded radial density variations known as zonal flows. Moreover, we report that vertical magnetic flux is strongly concentrated toward low-density regions of the zonal flow. Mean vertical magnetic field can be more than doubled in low-density regions, and reduced to nearly zero in high density regions in some cases. In ideal MHD, the scale on which magnetic flux concentrates can reach a few disk scale heights. In the non-ideal MHD regime with strong ambipolar diffusion, magnetic flux is concentrated into thin axisymmetric shells at some enhanced level, whose size is typically less than half a scale height. We show that magnetic flux concentration is closely related to the fact that the magnetic diffusivity of the MRI turbulence is anisotropic. In addition to a conventional Ohmic-like turbulent resistivity, we find that there is a correlation between the vertical velocity and horizontal magnetic field fluctuations that produces a mean electric field that acts to anti-diffuse the vertical magnetic flux. The anisotropic turbulent diffusivity has analogies to the Hall effect, and may have important implications for magnetic flux transport in accretion disks. The physical origin of magnetic flux concentration may be related to the development of channel flows followed by magnetic reconnection, which acts to decrease the mass-to-flux ratio in localized regions. The association of enhanced zonal flows with magnetic flux concentration may lead to global pressure bumps in protoplanetary disks that helps trap dust particles and facilitates planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Local Study of Accretion Disks with a Strong Vertical Magnetic Field: Magnetorotational Instability and Disk Outflow

We perform 3D vertically-stratified local shearing-box ideal MHD simulations of the magnetorotational instability (MRI) that include a net vertical magnetic flux, which is characterized by beta_0 (ratio of gas pressure to magnetic pressure of the net vertical field at midplane). We have considered beta_0=10^2, 10^3 and 10^4 and in the first two cases the most unstable linear MRI modes are well resolved in the simulations. We find that the behavior of the MRI turbulence strongly depends on beta_0: The radial transport of angular momentum increases with net vertical flux, achieving alpha=0.08 for beta_0=10^4 and alpha>1.0 for beta_0=100, where alpha is the Shakura-Sunyaev parameter. A critical value lies at beta_0=10^3: For beta_0>10^3, the disk consists of a gas pressure dominated midplane and a magnetically dominated corona. The turbulent strength increases with net flux, and angular momentum transport is dominated by turbulent fluctuations. The magnetic dynamo that leads to cyclic flips of large-scale fields still exists, but becomes more sporadic as net flux increases. For beta_0<10^3, the entire disk becomes magnetic dominated. The turbulent strength saturates, and the magnetic dynamo is quenched. Stronger large-scale fields are generated with increasing net flux, which dominates angular momentum transport. A strong outflow is launched from the disk by the magnetocentrifugal mechanism, and the mass flux increases linearly with net vertical flux and shows sign of saturation at beta_0=10^2. However, the outflow is unlikely to be directly connected to a global wind: for beta_0>10^3, the large-scale field has no permanent bending direction due to dynamo activities, while for beta_0<10^3, the outflows from the top and bottom sides of the disk bend towards opposite directions, inconsistent with a physical disk wind geometry. Global simulations are needed to address the fate of the outflow.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics

We describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for the study of astrophysical plasma dynamics. The code incorporates an energy-conserving particle integrator into a stable, second-order--accurate, three-stage predictor-predictor-corrector integration algorithm. The constrained transport method is used to enforce the divergence-free constraint on the magnetic field. A delta-f scheme is included to facilitate a reduced-noise study of systems in which only small departures from an initial distribution function are anticipated. The effects of rotation and shear are implemented through the shearing-sheet formalism with orbital advection. These algorithms are embedded within an architecture similar to that used in the popular astrophysical magnetohydrodynamics code Athena, one that is modular, well-documented, easy to use, and efficiently parallelized for use on thousands of processors. We present a series of tests in one, two, and three spatial dimensions that demonstrate the fidelity and versatility of the code.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Turbulence In the Outer Regions of Protoplanetary Disks. I. Weak Accretion with No Vertical Magnetic Flux

We use local numerical simulations to investigate the strength and nature of magnetohydrodynamic (MHD) turbulence in the outer regions of protoplanetary disks, where ambipolar diffusion is the dominant non-ideal MHD effect. The simulations include vertical stratification and assume zero net vertical magnetic flux. We employ a super time-stepping technique to ameliorate the Courant restriction on the diffusive time step. We find that in idealized stratified simulations, with a spatially constant ambipolar Elsasser number Am, turbulence driven by the magnetorotational instability (MRI) behaves in a similar manner as in prior unstratified calculations. Turbulence dies away for Am < 1, and becomes progressively more vigorous as ambipolar diffusion is decreased. Near-ideal MHD behavior is recovered for Am > 1000. In the intermediate regime (10 < Am < 1000) ambipolar diffusion leads to substantial increases in both the period of the MRI dynamo cycle and the characteristic scales of magnetic field structures. To quantify the impact of ambipolar physics on disk accretion, we run simulations at 30 AU and 100 AU that include a vertical Am profile based upon far ultraviolet (FUV) ionized disk models. These models develop a vertically layered structure analogous to the Ohmic dead zone that is present at smaller radii. We find that, although the levels of surface turbulence can be strong (and consistent with constraints on turbulent line widths at these radii), the inferred accretion rates are at least an order of magnitude smaller than those observed in T Tauri stars. This discrepancy is very likely due to the assumption of zero vertical magnetic field in our simulations and suggests that vertical magnetic fields are essential for MRI-driven accretion in the outer regions of protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Turbulence In the Outer Regions of Protoplanetary Disks. II. Strong Accretion Driven by a Vertical Magnetic Field

We carry out a series of local, vertically stratified shearing box simulations of protoplanetary disks that include ambipolar diffusion and a net vertical magnetic field. The ambipolar diffusion profiles we employ correspond to 30AU and 100AU in a minimum mass solar nebula (MMSN) disk model, which consists of a far-UV-ionized surface layer and low-ionization disk interior. These simulations serve as a follow up to Simon et al. (2013), in which we found that without a net vertical field, the turbulent stresses that result from the magnetorotational instability (MRI) are too weak to account for observed accretion rates. The simulations in this work show a very strong dependence of the accretion stresses on the strength of the background vertical field; as the field strength increases, the stress amplitude increases. For gas to magnetic pressure ratios of 1e4 and 1e5, we find accretion rates between 1e-8 and 1e-7 solar masses per year. These accretion rates agree with observational constraints, suggesting a vertical magnetic field strength between 60 and 200 microgauss at 30AU and 10 and 30 microgauss at 100AU in a MMSN disk. Furthermore, the stress has a non-negligible component due to a magnetic wind. For sufficiently strong vertical field strengths, MRI turbulence is quenched, and the flow becomes largely laminar, with accretion proceeding through large scale correlations in the radial and toroidal field components as well as through the magnetic wind. In all simulations, the presence of a low ionization region near the disk mid-plane, which we call the ambipolar damping zone, results in reduced stresses there.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Wind-driven Accretion in Protoplanetary Disks --- II: Radial Dependence and Global Picture

Non-ideal magnetohydrodynamical effects play a crucial role in determining the mechanism and efficiency of angular momentum transport as well as the level of turbulence in protoplanetary disks (PPDs), which are key to understanding PPD evolution and planet formation. It was shown in our previous work that at 1 AU, the magnetorotational instability (MRI) is completely suppressed when both Ohmic resistivity and ambipolar diffusion (AD) are taken into account, resulting in a laminar flow with accretion driven by magnetocentrifugal wind. In this work, we study the radial dependence of the laminar wind solution using local shearing-box simulations. Scaling relation on the angular momentum transport for the laminar wind is obtained, and we find that the wind-driven accretion rate can be approximated as M_dot~0.91x10^(-8)R_AU^(1.21)(B_z/10mG)^(0.93)M_Sun/yr, where B_z is the strength of the large-scale vertical magnetic field threading the disk. The result is independent of disk surface density. Four criteria are outlined for the existence of the laminar wind solution: 1). Ohmic resistivity dominated midplane region; 2). AD dominated disk upper layer; 3). Presence of (not too weak) net vertical magnetic flux. 4). Sufficiently well ionized gas beyond disk surface. All these criteria are likely to be met in the inner region of the disk from ~0.3 AU to about 5-10 AU for typical PPD accretion rates. Beyond this radius, angular momentum transport is likely to proceed due to a combination of the MRI and disk wind, and eventually dominated by the MRI (in the presence of strong AD) in the outer disk. Our simulation results provide key ingredients for a new paradigm on the accretion processes in PPDs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Wind-driven Accretion in Protoplanetary Disks. I: Suppression of the Magnetorotational Instability and Launching of the Magnetocentrifugal Wind

We perform local, vertically stratified shearing-box MHD simulations of protoplanetary disks (PPDs) at a fiducial radius of 1 AU that take into account the effects of both Ohmic resistivity and ambipolar diffusion (AD). The magnetic diffusion coefficients are evaluated self-consistently from a look-up table based on equilibrium chemistry. We first show that the inclusion of AD dramatically changes the conventional picture of layered accretion. Without net vertical magnetic field, the system evolves into a toroidal field dominated configuration with extremely weak turbulence in the far-UV ionization layer that is far too inefficient to drive rapid accretion. In the presence of a weak net vertical field (plasma beta~10^5 at midplane), we find that the MRI is completely suppressed, resulting in a fully laminar flow throughout the vertical extent of the disk. A strong magnetocentrifugal wind is launched that efficiently carries away disk angular momentum and easily accounts for the observed accretion rate in PPDs. Moreover, under a physical disk wind geometry, all the accretion flow proceeds through a strong current layer with thickness of ~0.3H that is offset from disk midplane with radial velocity of up to 0.4 times the sound speed. Both Ohmic resistivity and AD are essential for the suppression of the MRI and wind launching. The efficiency of wind transport increases with increasing net vertical magnetic flux and the penetration depth of the FUV ionization. Our laminar wind solution has important implications on planet formation and global evolution of PPDs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hot and Cold Galactic Gas in the NGC 2563 Galaxy Group

The role of environmentally induced gas stripping in driving galaxy evolution in groups remains poorly understood. Here we present extensive Chandra and Very Large Array mosaic observations of the hot and cold interstellar medium within the members of the nearby, X-ray bright NGC 2563 group, a prime target for studies of the role of gas stripping and interactions in relatively small host halos. Our observations cover nearly all group members within a projected radius of 1.15 Mpc (~1.4 R_vir) of the group center, down to a limiting X-ray luminosity and HI mass of 3e39 erg/s and 2e8 M_sun, respectively. The X-ray data are consistent with efficient ram pressure stripping of the hot gas halos of early-type galaxies near the group core, but no X-ray tails are seen and the limited statistics preclude strong conclusions. The HI results suggest moderate HI mass loss from the group members when compared to similar field galaxies. Six of the 20 HI-detected group members show HI evidence of ongoing interactions with other galaxies or with the intragroup medium. Suggestive evidence is further seen for galaxies with close neighbors in position-velocity space to show relatively low HI content, consistent with tidal removal of HI. The results thus indicate removal of both hot and cold gas from the group members via a combination of ram pressure stripping and tidal interactions. We also find that 16 of the 20 HI detections occur on one side of the group, reflecting an unusual morphological segregation whose origin remains unclear.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Effect of Ambipolar Diffusion on the Non-linear Evolution of Magnetorotational Instability in Weakly Ionized Disks

We study the role of ambipolar diffusion (AD) on the non-linear evolution of the MRI in protoplanetary disks using the strong coupling limit, which applies when the electron recombination time is much shorter than the orbital time. The effect of AD in this limit is characterized by the dimensionless number Am, the frequency of which neutral particles collide with ions normalized to the orbital frequency. We perform three-dimensional unstratified shearing-box simulations of the MRI over a wide range of Am as well as different magnetic field strengths and geometries. The saturation level of the MRI turbulence depends on the magnetic geometry and increases with the net magnetic flux. There is an upper limit to the net flux for sustained turbulence, corresponding to the requirement that the most unstable vertical wavelength be less than the disk scale height. Correspondingly, at a given Am, there exists a maximum value of the turbulent stress alpha_max. For Am<1, the largest stress is associated with a field geometry that has both net vertical and toroidal flux. In this case, we confirm the results of linear analyses that show the fastest growing mode has a non-zero radial wave number with growth rate exceeding the pure vertical field case. We find there is a very tight correlation between the turbulent stress (alpha) and the plasma beta=P_gas/P_mag~1/(2alpha) at the saturated state of the MRI turbulence regardless of field geometry, and alpha_max rapidly decreases with decreasing Am. In particular, we quote alpha_max~0.007 for Am=1 and alpha_max~0.0006 for Am=0.1.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetorotational Instability Driven Accretion in Protoplanetary Disks

The gas dynamics in protoplanetary disks (PPDs) is strongly affected by non-ideal MHD effects. Using a complex chemical reaction network with standard prescriptions for X-ray and cosmic-ray ionizations, as well as the most up-to-date results from numerical simulations, we study the non-ideal MHD effects on the magnetorotational instability (MRI) and angular momentum transport in PPDs. We first show that no matter grains are included or not, the recombination time is always shorter than the orbital time in the bulk of PPDs, justifying the validity of local ionization equilibrium. The full conductivity tensor at different disk radii and heights is evaluated, with the MRI active region determined by requiring that (1) the Ohmic Elsasser number be greater than 1; (2) the ratio of gas to magnetic pressure beta be greater than beta_min(Am) as identified by Bai & Stone (2011), where Am is the Elsasser number for ambipolar diffusion. With full flexibility as to the magnetic field strength, we provide a general framework for estimating the MRI-driven accretion rate M_dot and the magnetic field strength in the MRI-active layer. We find that the MRI-active layer always exists at any disk radii as long as the magnetic field in PPDs is sufficiently weak. However, the optimistically predicted M_dot in the inner disk (r=1-10 AU) appears insufficient to account for the observed range of accretion rate in PPDs (around 1e-8M_Sun/yr) even in the grain-free calculation, and the presence of solar abundance sub-micron grains further reduces M_dot by one to two orders of magnitude. Our results suggest that stronger sources of ionization, and/or additional mechanisms such as magnetized wind are needed to explain the observed accretion rates in PPDs. In contrast, our predicted M_dot is on the order of 1e-9M_Sun/yr in the outer disk, consistent with the observed accretion rates in transitional disks. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Role of Tiny Grains on the Accretion Process in Protoplanetary Disks

Tiny grains such as PAHs have been thought to dramatically reduce the coupling between gas and magnetic fields in weakly ionized gas such as in protoplanetary disks (PPDs) because they provide tremendous surface area to recombine free electrons. The presence of tiny grains in PPDs thus raises the question of whether the magnetorotational instability (MRI) is able to drive rapid accretion to be consistent with observations. Charged tiny grains have similar conduction properties as ions, whose presence leads to qualitatively new behaviors in the conductivity tensor, characterized by n_bar/n_e>1, where n_e and n_bar denote the number densities of free electrons and all other charged species respectively. In particular, Ohmic conductivity becomes dominated by charged grains rather than electrons when n_bar/n_e exceeds about 10^3, and Hall and ambipolar diffusion (AD) coefficients are reduced by a factor of (n_bar/n_e)^2 in the AD dominated regime relative to that in the Ohmic regime. Applying the methodology of Bai (2011), we find that in PPDs, when PAHs are sufficiently abundant (>1e-9 per H_2), there exists a transition radius r_trans of about 10-20 AU, beyond which the MRI active layer extends to the disk midplane. At r<r_trans, the optimistically predicted MRI-driven accretion rate M_dot is one to two orders of magnitude smaller than that in the grain-free case, which is too small compared with the observed rates, but is in general no smaller than the predicted M_dot with solar-abundance 0.1 micron grains. At r>r_trans, we find that remarkably, the predicted M_dot exceeds the grain-free case due to a net reduction of AD by charged tiny grains, and reaches a few times 1e-8M_Sun/yr. This is sufficient to account for the observed M_dot in transitional disks. Larger grains (>0.1 micron) are too massive to reach such high abundance as tiny grains and to facilitate the accretion process.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamics of Solids in the Midplane of Protoplanetary Disks: Implications for Planetesimal Formation

(Abridged) We present local 2D and 3D hybrid numerical simulations of particles and gas in the midplane of protoplanetary disks (PPDs) using the Athena code. The particles are coupled to gas aerodynamically, with particle-to-gas feedback included. Magnetorotational turbulence is ignored as an approximation for the dead zone of PPDs, and we ignore particle self-gravity to study the precursor of planetesimal formation. Our simulations include a wide size distribution of particles, ranging from strongly coupled particles with dimensionless stopping time tau_s=Omega t_stop=1e-4 to marginally coupled ones with tau_s=1 (where Omega is the orbital frequency, t_stop is the particle friction time), and a wide range of solid abundances. Our main results are: 1. Particles with tau_s>=0.01 actively participate in the streaming instability, generate turbulence and maintain the height of the particle layer before Kelvin-Helmholtz instability is triggered. 2. Strong particle clumping as a consequence of the streaming instability occurs when a substantial fraction of the solids are large (tau_s>=0.01) and when height-integrated solid to gas mass ratio Z is super-solar. 3. The radial drift velocity is reduced relative to the conventional Nakagawa-Sekiya-Hayashi (NSH) model, especially at high Z. We derive a generalized NSH equilibrium solution for multiple particle species which fits our results very well. 4. Collision velocity between particles with tau_s>=0.01 is dominated by differential radial drift, and is strongly reduced at larger Z. 5. There exist two positive feedback loops with respect to the enrichment of local disk solid abundance and grain growth. All these effects promote planetesimal formation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling of Gamma-Ray Pulsar Light Curves with Force-Free Magnetic Field

(Abridged) Gamma-ray emission from pulsars has long been modeled using a vacuum dipole field. This approximation ignores changes in the field structure caused by the magnetospheric plasma and strong plasma currents. We present the first results of gamma-ray pulsar light curve modeling using the more realistic field taken from 3D force-free magnetospheric simulations. Having the geometry of the field, we apply several prescriptions for the location of the emission zone, comparing the light curves to observations. We find that the conventional two-pole caustic model fails to produce double-peak pulse profiles, mainly because the size of the polar cap in force-free magnetosphere is larger than the vacuum field polar cap. The conventional outer-gap model is capable of producing only one peak under general conditions, because a large fraction of open field lines does not cross the null charge surface. We propose a novel "separatrix layer" model, where the high-energy emission originates from a thin layer on the open field lines just inside of the separatrix that bounds the open flux tube. The emission from this layer generates two strong caustics on the sky map due to the effect we term "Sky Map Stagnation" (SMS). It is related to the fact that force-free field asymptotically approaches the field of a rotating split monopole, and the photons emitted on such field lines in the outer magnetosphere arrive to the observer in phase. The double-peak light curve is a natural consequence of SMS. We show that most features of the currently available gamma-ray pulsar light curves can be reasonably well reproduced and explained with the sepatratrix model using the force-free field. Association of the emission region with the current sheet will guide more detailed future studies of the magnetospheric acceleration physics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Particle-Gas Dynamics with Athena: Method and Convergence

The Athena MHD code has been extended to integrates the motion of particles coupled with the gas via aerodynamic drag, in order to study the dynamics of gas and solids in protoplanetary disks and the formation of planetesimals. Our particle-gas hybrid scheme is based on a second order predictor-corrector method. Careful treatment of the momentum feedback on the gas guarantees exact conservation. The hybrid scheme is stable and convergent in most regimes relevant to protoplanetary disks. We describe a semi-implicit integrator generalized from the leap-frog approach. In the absence of drag force, it preserves the geometric properties of a particle orbit. We also present a fully-implicit integrator that is unconditionally stable for all regimes of particle-gas coupling. Using our hybrid code, we study the numerical convergence of the non-linear saturated state of the streaming instability. We find that gas flow properties are well converged with modest grid resolution (128 cells per pressure length ηr for dimensionless stopping time tau_s=0.1), and equal number of particles and grid cells. On the other hand, particle clumping properties converge only at higher resolutions, and finer resolution leads to stronger clumping before convergence is reached. Finally, we find that measurement of particle transport properties resulted from the streaming instability may be subject to error of about 20%.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Effect of the Radial Pressure Gradient in Protoplanetary Disks on Planetesimal Formation

The streaming instability (SI) provides a promising mechanism for planetesimal formation because of its ability to concentrate solids into dense clumps. The degree of clumping strongly depends on the height-integrated solid to gas mass ratio Z in protoplanetary disks (PPDs). In this letter, we show that the magnitude of the radial pressure gradient (RPG) which drives the SI (characterized by $q\equivηv_K/c_s$, where $ηv_K$ is the reduction of Keplerian velocity due to the RPG and $c_s$ is the sound speed) also strongly affects clumping. We present local two-dimensional hybrid numerical simulations of aerodynamically coupled particles and gas in the midplane of PPDs. Magnetic fields and particle self-gravity are ignored. We explore three different RPG values appropriate for typical PPDs: $q=0.025, 0.05$ and 0.1. For each $q$ value, we consider four different particle size distributions ranging from sub millimeter to meter sizes and run simulations with solid abundance from Z=0.01 up to Z=0.07. We find that a small RPG strongly promotes particle clumping in that: 1) At fixed particle size distribution, the critical solid abundance $Z_{\rm crit}$ above which particle clumping occurs monotonically increases with $q$; 2) At fixed Z, strong clumping can occur for smaller particles when $q$ is smaller. Therefore, we expect planetesimals to form preferentially in regions of PPDs with a small RPG.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Uncertainties of Modeling Gamma-Ray Pulsar Light Curves with Vacuum Dipole Magnetic Field

Current models of pulsar gamma-ray emission use the magnetic field of a rotating dipole in vacuum as a first approximation to the shape of plasma-filled pulsar magnetosphere. In this paper we revisit the question of gamma-ray light-curve formation in pulsars in order to ascertain the robustness of the "two-pole caustic" (TPC) and "outer gap" (OG) models based on the vacuum magnetic field. We point out an inconsistency in the literature on the use of the relativistic aberration formula, where in several works the vacuum field was treated as known in the instantaneously corotating frame, rather than in the laboratory frame. With the corrected formula, we find that the peaks in the light curves predicted from the TPC model using the vacuum field are less sharp. The sharpness of the peaks in the OG model is less affected by this change, but the range of magnetic inclination angles and viewing geometries resulting in double-peaked light curves is reduced. In a realistic magnetosphere, the modification of field structure near the light cylinder due to plasma effects may change the shape of the polar cap and the location of the emission zones. We study the sensitivity of the light curves to different shapes of the polar cap for static and retarded vacuum dipole fields. In particular, we consider polar caps traced by the last open field lines and compare them to circular polar caps. We find that the TPC model is very sensitive to the shape of the polar cap, and a circular polar cap can lead to four peaks of emission. The OG model is less affected by polar cap shapes, but is subject to big uncertainties of applying the vacuum field near the light cylinder. We conclude that deviations from vacuum field can lead to large uncertainties in the pulse shapes, and a more realistic force-free field should be applied to the study of pulsar high energy emission.

0 citations0 reviewsNo code linkedOpen access