Source author record

James M. Stone

James M. Stone appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph.EP astro-ph.SR astro-ph astro-ph.IM astro-ph.GA physics.comp-ph astro-ph.CO physics.plasm-ph physics.space-ph gr-qc physics.flu-dyn physics.optics

Catalog footprint

What is connected

59works

13topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2023arXiv

Dissipation of AGN jets in a clumpy interstellar medium

Accreting supermassive black holes (SMBHs) frequently power jets that interact with the interstellar/circumgalactic medium (ISM/CGM), regulating star-formation in the galaxy. Highly supersonic jets launched by active galactic nuclei (AGN) power a cocoon that confines them and shocks the ambient medium. We build upon the models of narrow conical jets interacting with a smooth ambient medium, to include the effect of dense clouds that are an essential ingredient of a multiphase ISM. The key physical ingredient of this model is that the clouds along the supersonic jet-beam strongly decelerate the jet-head, but the subsonic cocoon easily moves around the clouds without much resistance. We propose scalings for important physical quantities -- cocoon pressure, head & cocoon speed, and jet radius. We obtain, for the first time, the analytic condition on clumpiness of the ambient medium for the jet to dissipate within the cocoon and verify it with numerical simulations of conical jets interacting with a uniform ISM with embedded spherical clouds. A jet is defined to be dissipated when the cocoon speed exceeds the speed of the jet-head. We compare our models to more sophisticated numerical simulations, direct observations of jet-ISM interaction (e.g., quasar J1316+1753), and discuss implications for the Fermi/eROSITA bubbles. Our work also motivates effective subgrid models for AGN jet feedback in a clumpy ISM unresolved by the present generation of cosmological galaxy formation simulations.

preprint2020arXiv

Ab Initio Horizon-Scale Simulations of Magnetically Arrested Accretion in Sagittarius A* Fed by Stellar Winds

We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $\sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained properties of the WR winds with limited free parameters. For this initial study we assume a non-spinning black hole. Our simulations naturally produce a $\sim 10^{-8} M_\odot$ yr$^{-1}$ accretion rate, consistent with previous phenomenological estimates. We find that a magnetically arrested flow is formed by the continuous accretion of coherent magnetic field being fed from large radii. Near the event horizon, the magnetic field is so strong that it tilts the gas with respect to the initial angular momentum and concentrates the originally quasi-spherical flow to a narrow disk-like structure. We also present 230 GHz images calculated from our simulations where the inclination angle and physical accretion rate are not free parameters but are determined by the properties of the WR stellar winds. The image morphology is highly time variable. Linear polarization on horizon scales is coherent with weak internal Faraday rotation.

preprint2020arXiv

Global 3-D Radiation Magnetohydrodynamic Simulations for FU Ori's Accretion Disk and Observational Signatures of Magnetic Fields

FU Ori is the prototype of FU Orionis systems which are outbursting protoplanetary disks. Magnetic fields in FU Ori's accretion disks have previously been detected using spectropolarimetry observations for Zeeman effects. We carry out global radiation ideal MHD simulations to study FU Ori's inner accretion disk. We find that (1) when the disk is threaded by vertical magnetic fields, most accretion occurs in the magnetically dominated atmosphere at z$\sim$R, similar to the "surface accretion" mechanism in previous locally-isothermal MHD simulations. (2) A moderate disk wind is launched in the vertical field simulations with a terminal speed of $\sim$300-500 km/s and a mass loss rate of 1-10\% the disk accretion rate, which is consistent with observations. Disk wind fails to be launched in simulations with net toroidal magnetic fields. (3) The disk photosphere at the unit optical depth can be either in the wind launching region or the accreting surface region. Magnetic fields have drastically different directions and magnitudes between these two regions. Our fiducial model agrees with previous optical Zeeman observations regarding both the field directions and magnitudes. On the other hand, simulations indicate that future Zeeman observations at near-IR wavelengths or towards other FU Orionis systems may reveal very different magnetic field structures. (4) Due to energy loss by the disk wind, the disk photosphere temperature is lower than that predicted by the thin disk theory, and the previously inferred disk accretion rate may be lower than the real accretion rate by a factor of $\sim$2-3.

preprint2020arXiv

The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers

The design and implementation of a new framework for adaptive mesh refinement (AMR) calculations is described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic execution model based on a simple design called a "task list" that improves parallel performance by overlapping communication and computation, simplifies the inclusion of a diverse range of physics, and even enables multiphysics models involving different physics in different regions of the calculation. We describe physics modules implemented in this framework for both non-relativistic and relativistic magnetohydrodynamics (MHD). These modules adopt mature and robust algorithms originally developed for the Athena MHD code and incorporate new extensions: support for curvilinear coordinates, higher-order time integrators, more realistic physics such as a general equation of state, and diffusion terms that can be integrated with super-time-stepping algorithms. The modules show excellent performance and scaling, with well over 80% parallel efficiency on over half a million threads. The source code has been made publicly available.

preprint2020arXiv

The Surprisingly Small Impact of Magnetic Fields On The Inner Accretion Flow of Sagittarius A* Fueled By Stellar Winds

We study the flow structure in 3D magnetohydrodynamic (MHD) simulations of accretion onto Sagittarius A* via the magnetized winds of the orbiting Wolf-Rayet stars. These simulations cover over 3 orders of magnitude in radius to reach $\approx$ 300 gravitational radii, with only one poorly constrained parameter (the magnetic field in the stellar winds). Even for winds with relatively weak magnetic fields (e.g., plasma $β$ $\sim$ $10^6$), flux freezing/compression in the inflowing gas amplifies the field to $β$ $\sim$ few well before it reaches the event horizon. Overall, the dynamics, accretion rate, and spherically averaged flow profiles (e.g., density, velocity) in our MHD simulations are remarkably similar to analogous hydrodynamic simulations. We attribute this to the broad distribution of angular momentum provided by the stellar winds, which sources accretion even absent much angular momentum transport. We find that the magneto-rotational instability is not important because of i) strong magnetic fields that are amplified by flux freezing/compression, and ii) the rapid inflow/outflow times of the gas and inefficient radiative cooling preclude circularization. The primary effect of magnetic fields is that they drive a polar outflow that is absent in hydrodynamics. The dynamical state of the accretion flow found in our simulations is unlike the rotationally supported tori used as initial conditions in horizon scale simulations, which could have implications for models being used to interpret Event Horizon Telescope and GRAVITY observations of Sgr A*.

preprint2016arXiv

An Extension of the Athena++ Code Framework for GRMHD Based on Advanced Riemann Solvers and Staggered-Mesh Constrained Transport

We present a new general relativistic magnetohydrodynamics (GRMHD) code integrated into the Athena++ framework. Improving upon the techniques used in most GRMHD codes, ours allows the use of advanced, less diffusive Riemann solvers, in particular HLLC and HLLD. We also employ a staggered-mesh constrained transport algorithm suited for curvilinear coordinate systems in order to maintain the divergence-free constraint of the magnetic field. Our code is designed to work with arbitrary stationary spacetimes in one, two, or three dimensions, and we demonstrate its reliability in a number of tests. We also report on its promising performance and scalability.

preprint2016arXiv

Dust dynamics in 2D gravito-turbulent disks

The dynamics of solid bodies in protoplanetary disks are subject to the properties of any underlying gas turbulence. Turbulence driven by disk self-gravity shows features distinct from those driven by the magnetorotational instability (MRI). We study the dynamics of solids in gravito-turbulent disks with two-dimensional (in the disk plane), hybrid (particle and gas) simulations. Gravito-turbulent disks can exhibit stronger gravitational stirring than MRI-active disks, resulting in greater radial diffusion and larger eccentricities and relative speeds for large particles (those with dimensionless stopping times $t_{stop} Ω> 1$, where $Ω$ is the orbital frequency). The agglomeration of large particles into planetesimals by pairwise collisions is therefore disfavored in gravito-turbulent disks. However, the relative speeds of intermediate-size particles $t_{stop} Ω\sim 1$ are significantly reduced as such particles are collected by gas drag and gas gravity into coherent filament-like structures with densities high enough to trigger gravitational collapse. First-generation planetesimals may form via gravitational instability of dust in marginally gravitationally unstable gas disks.

preprint2016arXiv

Global MHD Simulations of Accretion Disks in Cataclysmic Variables (CVs): I. The Importance of Spiral Shocks

We present results from the first global 3D MHD simulations of accretion disks in Cataclysmic Variable (CV) systems in order to investigate the relative importance of angular momentum transport via turbulence driven by the magnetorotational instability (MRI) compared to that driven by spiral shock waves. Remarkably, we find that even with vigorous MRI turbulence, spiral shocks are an important component to the overall angular momentum budget, at least when temperatures in the disk are high (so that Mach numbers are low). In order to understand the excitation, propagation, and damping of spiral density waves in our simulations more carefully, we perform a series of 2D global hydrodynamical simulations with various equation of states and both with and without mass inflow via the Lagrangian point (L1). Compared with previous similar studies, we find the following new results. 1) Linear wave dispersion relation fits the pitch angles of spiral density waves very well. 2) We demonstrate explicitly that mass accretion is driven by the deposition of negative angular momentum carried by the waves when they dissipate in shocks. 3) Using Reynolds stress scaled by gas pressure to represent the effective angular momentum transport rate alpha_{eff} is not accurate when mass accretion is driven by non-axisymmetric shocks. 4) Using the mass accretion rate measured in our simulations to directly measure alpha defined in standard thin-disk theory, we find 0.02 < alpha_{eff} < 0.05 for CV disks, consistent with observed values in quiescent states of dwarf novae (DNe). In this regime the disk may be too cool and neutral for the MRI to operate and spiral shocks are a possible accretion mechanism. However, we caution that our simulations use unrealistically low Mach numbers in this regime, and therefore future models with more realistic thermodynamics and non-ideal MHD are warranted.

preprint2016arXiv

Magnetorotational Turbulence and Dynamo in a Collisionless Plasma

We present results from the first 3D kinetic numerical simulation of magnetorotational turbulence and dynamo, using the local shearing-box model of a collisionless accretion disc. The kinetic magnetorotational instability grows from a subthermal magnetic field having zero net flux over the computational domain to generate self-sustained turbulence and outward angular-momentum transport. Significant Maxwell and Reynolds stresses are accompanied by comparable viscous stresses produced by field-aligned ion pressure anisotropy, which is regulated primarily by the mirror and ion-cyclotron instabilities through particle trapping and pitch-angle scattering. The latter endow the plasma with an effective viscosity that is biased with respect to the magnetic-field direction and spatio-temporally variable. Energy spectra suggest an Alfvén-wave cascade at large scales and a kinetic-Alfvén-wave cascade at small scales, with strong small-scale density fluctuations and weak non-axisymmetric density waves. Ions undergo non-thermal particle acceleration, their distribution accurately described by a kappa distribution. These results have implications for the properties of low-collisionality accretion flows, such as that near the black hole at the Galactic center.

preprint2016arXiv

Shock-driven Accretion in Circumplanetary Disks: Observables and Satellite Formation

Circumplanetary disks (CPDs) control the growth of planets, supply material for satellites to form, and provide observational signatures of young forming planets. We have carried out two dimensional hydrodynamical simulations with radiative cooling to study CPDs, and suggested a new mechanism to drive the disk accretion. Two spiral shocks are present in CPDs, excited by the central star. We find that spiral shocks can at least contribute to, if not dominate the angular momentum transport and energy dissipation in CPDs. Meanwhile, dissipation and heating by spiral shocks have a positive feedback on shock-driven accretion itself. As the disk is heated up by spiral shocks, the shocks become more open, leading to more efficient angular momentum transport. This shock driven accretion is, on the other hand, unsteady on a timescale of months/years due to production and destruction of vortices in disks. After being averaged over time, a quasi-steady accretion is reached from the planet's Hill radius all the way to the planet surface, and the disk $α$-coefficient characterizing angular momentum transport due to spiral shocks is $\sim$0.001-0.02. The disk surface density ranges from 10 to 1000 g cm$^{-2}$ in our simulations, which is at least 3 orders of magnitude smaller than the "minimum mass sub-nebula" model used to study satellite formation; instead it is more consistent with the "gas-starved" satellite formation model. Finally, we calculate the millimeter flux emitted by CPDs at ALMA and EVLA wavelength bands and predict the flux for several recently discovered CPD candidates, which suggests that ALMA is capable of discovering these accreting CPDs.

preprint2015arXiv

Saturation of the magnetorotational instability in the unstratified shearing box with zero net flux: convergence in taller boxes

Previous studies of the nonlinear regime of the magnetorotational instability in one particular type of shearing box model -- unstratified with no net magnetic flux -- find that without explicit dissipation (viscosity and resistivity) the saturation amplitude decreases with increasing numerical resolution. We show that this result is strongly dependent on the vertical aspect ratio of the computational domain $L_z/L_x$. When $L_z/L_x\lesssim 1$, we recover previous results. However, when the vertical domain is extended $L_z/L_x \gtrsim 2.5$, we find the saturation level of the stress is greatly increased (giving a ratio of stress to pressure $α\gtrsim 0.1$), and moreover the results are independent of numerical resolution. Consistent with previous results, we find that saturation of the MRI in this regime is controlled by a cyclic dynamo which generates patches of strong toroidal field that switches sign on scales of $L_x$ in the vertical direction. We speculate that when $L_z/L_x\lesssim 1$, the dynamo is inhibited by the small size of the vertical domain, leading to the puzzling dependence of saturation amplitude on resolution. We show that previous toy models developed to explain the MRI dynamo are consistent with our results, and that the cyclic pattern of toroidal fields observed in stratified shearing box simulations (leading to the so-called butterfly diagram) may also be related. In tall boxes the saturation amplitude is insensitive to whether or not explicit dissipation is included in the calculations, at least for large magnetic Reynolds and Prandtl number. Finally, we show MRI turbulence in tall domains has a smaller critical $\rm{Pm}_c$, and an extended lifetime compared to $L_z/L_x\lesssim 1$ boxes.

preprint2015arXiv

Spreading layers in accreting objects: role of acoustic waves for angular momentum transport, mixing and thermodynamics

Disk accretion at high rate onto a white dwarf or a neutron star has been suggested to result in the formation of a spreading layer (SL) - a belt-like structure on the object's surface, in which the accreted matter steadily spreads in the poleward (meridional) direction while spinning down. To assess its basic characteristics we perform two-dimensional hydrodynamic simulations of supersonic SLs in the relevant morphology with a simple prescription for cooling. We demonstrate that supersonic shear naturally present at the base of the SL inevitably drives sonic instability that gives rise to large scale acoustic modes governing the evolution of the SL. These modes dominate the transport of momentum and energy, which is intrinsically global and cannot be characterized via some form of local effective viscosity (e.g. $α$-viscosity). The global nature of the wave-driven transport should have important implications for triggering Type I X-ray bursts in low mass X-ray binaries. The nonlinear evolution of waves into a system of shocks drives effective re-arrangement (sensitively depending on thermodynamical properties of the flow) and deceleration of the SL, which ultimately becomes transonic and susceptible to regular Kelvin-Helmholtz instability. We interpret this evolution in terms of the global structure of the SL and suggest that mixing of the SL material with the underlying stellar fluid should become effective only at intermediate latitudes on the accreting object's surface, where the flow has decelerated appreciably. In the near-equatorial regions the transport is dominated by acoustic waves and mixing is less efficient. We speculate that this latitudinal non-uniformity of mixing in accreting white dwarfs may be linked to the observed bipolar morphology of classical novae ejecta.

preprint2015arXiv

Survival and Structure of Dusty Vortices in Protoplanetary Discs

We have studied the impact of dust feedback on the survival and structure of vortices in protoplanetary discs using 2-D shearing box simulations with Lagrangian dust particles. We consider dust with a variety of sizes (stopping time $t_s = 10^{-2}Ω^{-1} - 10^{2}Ω^{-1}$, from fully coupled with the gas to the decoupling limit. We find that a vortex is destroyed by dust feedback when the total dust-to-gas mass ratio within the vortex is larger than 30-50%, independent of the dust size. The dust distribution can still be asymmetric in some cases after the vortex has been destroyed. With smaller amounts of dust, a vortex can survive for at least 100 orbits, and the maximum dust surface density within the vortex can be more than 100 times larger than the gas surface density, potentially facilitating planetesimal formation. On the other hand, in these stable vortices, small ($t_s < Ω^{-1}$) and large ($t_s > Ω^{-1}$) dust grains concentrate differently and affect the gas dynamics in different ways. The distribution of large dust is more elongated than that of small dust. Large dust ($t_s > Ω^{-1}$) concentrates in the centre of the vortex and feedback leads to turn-over in vorticity towards the centre, forming a quiescent region within an anticyclonic vortex. Such a turn-over is absent if the vortex is loaded with small grains. We demonstrate that, in protoplanetary discs where both large and small dust grains are present and under the right condition, the concentration of large dust towards the vortex centre can lead to a quiescent centre, repelling the small dust and forming a small dust ring around the vortex centre. Such anticorrelations between small and large dust within vortices may explain the discrepancy between ALMA and near-IR scattered light observations in the asymmetric region of transitional discs.

preprint2015arXiv

The Structure of Spiral Shocks Excited by Planetary-mass Companions

Direct imaging observations have revealed spiral structures in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain some of these spiral patterns, due to the large pitch angle and high contrast of the spiral arms in observations. We have carried out three dimensional (3-D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the non-linear density wave theory. A secondary (or even a tertiary) spiral arm, especially for inner arms, is also excited by a massive planet. With a more massive planet in the disk, the excited spiral arms have larger pitch angle and the separation between the primary and secondary arms in the azimuthal direction is also larger. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3-D hydrodynamical models than images based on 2-D models assuming vertical hydrostatic equilibrium, indicating the need to model observations with full 3-D hydrodynamics. Overall, companion-induced spiral arms not only pinpoint the companion's position but also provide three independent ways (pitch angle, separation between two arms, and contrast of arms) to constrain the companion's mass.

preprint2014arXiv

A Global Three Dimensional Radiation Magneto-hydrodynamic Simulation of Super-Eddington Accretion Disks

We study super-Eddington accretion flows onto black holes using a global three dimensional radiation magneto-hydrodynamical simulation. We solve the time dependent radiative transfer equation for the specific intensities to accurately calculate the angular distribution of the emitted radiation. Turbulence generated by the magneto-rotational instability provides self-consistent angular momentum transfer. The simulation reaches inflow equilibrium with an accretion rate ~220L_edd/c^2 and forms a radiation driven outflow along the rotation axis. The mechanical energy flux carried by the outflow is ~20% of the radiative energy flux. The total mass flux lost in the outflow is about 29% of the net accretion rate. The radiative luminosity of this flow is ~10L_edd. This yields a radiative efficiency ~4.5%, which is comparable to the value in a standard thin disk model. In our simulation, vertical advection of radiation caused by magnetic buoyancy transports energy faster than photon diffusion, allowing a significant fraction of the photons to escape from the surface of the disk before being advected into the black hole. We contrast our results with the lower radiative efficiencies inferred in most models, such as the slim disk model, which neglect vertical advection. Our inferred radiative efficiencies also exceed published results from previous global numerical simulations, which did not attribute a significant role to vertical advection. We briefly discuss the implications for the growth of supermassive black holes in the early universe and describe how these results provided a basis for explaining the spectrum and population statistics of ultraluminous X-ray sources.

preprint2014arXiv

An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation

(Abridged) We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics when the light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the Athena code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers. The method is ideally suited for problems in which characteristic velocities are non-relativistic, but still within a few percent or more of the speed of light. The method is an intermediate step towards algorithms for fully relativistic flows.

preprint2014arXiv

Characterizing the variation of propagation constants in multicore fibre

We demonstrate a numerical technique that can evaluate the core-to-core variations in propagation constant in multicore fibre. Using a Markov Chain Monte Carlo process, we replicate the interference patterns of light that has coupled between the cores during propagation. We describe the algorithm and verify its operation by successfully reconstructing target propagation constants in a fictional fibre. Then we carry out a reconstruction of the propagation constants in a real fibre containing 37 single-mode cores. We find that the range of fractional propagation constant variation across the cores is approximately $\pm2 \times 10^{-5}$.

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.

preprint2014arXiv

Dust Trapping by Vortices in Transitional Disks: Evidence for Non-ideal MHD Effects in Protoplanetary Disks

We study particle trapping at the edge of a gap opened by a planet in a protoplanetary disk. In particular, we explore the effects of turbulence driven by the magnetorotational instability on particle trapping, using global three-dimensional magnetohydrodynamic (MHD) simulations including Lagrangian dust particles. We study disks either in the ideal MHD limit or dominated by ambipolar diffusion (AD) that plays an essential role at the outer regions of a protoplanetary disk. With ideal MHD, strong turbulence (the equivalent viscosity parameter $α\sim 10^{-2}$) in disks prevents vortex formation at the edge of the gap opened by a 9 $M_{J}$ planet, and most particles (except the particles that drift fastest) pile up at the outer gap edge almost axisymmetrically. When AD is considered, turbulence is significantly suppressed ($α\lesssim 10^{-3}$), and a large vortex forms at the edge of the planet induced gap, which survives $\sim$ 1000 orbits. The vortex can efficiently trap dust particles that span 3 orders of magnitude in size within 100 planetary orbits. We have also carried out two-dimensional hydrodynamical simulations using viscosity as an approximation to MHD turbulence. These hydrodynamical simulations can reproduce vortex generation at the gap edge as seen in MHD simulations. Finally, we use our simulation results to generate synthetic images for ALMA dust continuum observations on Oph IRS 48 and HD 142527, which show good agreement with existing observations. Predictions for future ALMA cycle 2 observations have been made. We conclude that the asymmetry in ALMA observations can be explained by dust trapping vortices and the existence of vortices could be the evidence that the outer protoplanetary disks are dominated by ambipolar diffusion with $α<10^{-3}$ at the disk midplane.

preprint2014arXiv

Firehose and Mirror Instabilities in a Collisionless Shearing Plasma

Hybrid-kinetic numerical simulations of firehose and mirror instabilities in a collisionless plasma are performed in which pressure anisotropy is driven as the magnetic field is changed by a persistent linear shear $S$. For a decreasing field, it is found that mostly oblique firehose fluctuations grow at ion Larmor scales and saturate with energies $\sim$$S^{1/2}$; the pressure anisotropy is pinned at the stability threshold by particle scattering off microscale fluctuations. In contrast, nonlinear mirror fluctuations are large compared to the ion Larmor scale and grow secularly in time; marginality is maintained by an increasing population of resonant particles trapped in magnetic mirrors. After one shear time, saturated order-unity magnetic mirrors are formed and particles scatter off their sharp edges. Both instabilities drive sub-ion-Larmor--scale fluctuations, which appear to be kinetic-Alfvén-wave turbulence. Our results impact theories of momentum and heat transport in astrophysical and space plasmas, in which the stretching of a magnetic field by shear is a generic process.

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.

preprint2014arXiv

Particle Concentration At Planet Induced Gap Edges and Vortices: I. Inviscid 3-D Hydro Disks

We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global 2-D and 3-D inviscid hydrodynamic simulations. Lagrangian particles have been implemented into magnetohydrodynamic code Athena with cylindrical coordinates. We find two distinct outcomes depending on the mass of the embedded planet. In the presence of a low mass planet ($8 M_{\oplus}$), two narrow gaps start to open in the gas on each side of the planet where the density waves shock. These shallow gaps can dramatically affect particle drift speed and cause significant, roughly axisymmetric dust depletion. On the other hand, a more massive planet ($>0.1 M_{J}$) carves out a deeper gap with sharp edges, which are unstable to the vortex formation. Particles with a wide range of sizes ($0.02<Ωt_{s}<20$) are trapped and settle to the midplane in the vortex, with the strongest concentration for particles with $Ωt_{s}\sim 1$. The dust concentration is highly elongated in the $ϕ$ direction, and can be as wide as 4 disk scale heights in the radial direction. Dust surface density inside the vortex can be increased by more than a factor of 10$^2$ in a very non-axisymmetric fashion. For very big particles ($Ωt_{s}\gg 1$) we find strong eccentricity excitation, in particular around the planet and in the vicinity of the mean motion resonances, facilitating gap opening there. Our results imply that in weakly turbulent protoplanetary disk regions (e.g. the "dead zone") dust particles with a very wide range of sizes can be trapped at gap edges and inside vortices induced by planets with $M_{p}<M_{J}$, potentially accelerating planetesimal and planet formation there, and giving rise to distinctive features that can be probed by ALMA and EVLA.

preprint2014arXiv

Radiation Feedback in ULIRGS: Are Photons Movers and Shakers?

We use our variable Eddington tensor (VET) radiation hydrodynamics code to perform two-dimensional simulations to study the impact of radiation forces on atmospheres composed of dust and gas. Our setup closely follows that of Krumholz & Thompson, assuming that dust and gas are well-coupled and that the radiation field is characterized by blackbodies with temperatures >~ 80 K, as might be found in ultraluminous infrared galaxies. In agreement with previous work, we find that Rayleigh-Taylor instabilities develop in radiation supported atmospheres, leading to inhomogeneities that limit momentum exchange between radiation and dusty gas, and eventually providing a near balance of the radiation and gravitational forces. However, the evolution of the velocity and spatial distributions of the gas differs significantly from previous work, which utilized a less accurate flux-limited diffusion (FLD) method. Our VET simulations show continuous net acceleration of the gas, with no steady-state reached by the end of the simulation. In contrast, FLD results show little net acceleration of the gas and settle in to a quasi-steady, turbulent state with low velocity dispersion. The discrepancies result primarily from the inability of FLD to properly model the variation of the radiation field around structures that are less than a few optical depths across. We conclude that radiation feedback remains a viable mechanism for driving high-Mach number turbulence. We discuss implications for observed systems and global numerical simulations of feedback, but more realistic setups are needed to make robust observational predictions and assess the prospect of launching outflows with radiation.

preprint2014arXiv

Radiation Magneto-hydrodynamic Simulations of the Formation of Hot Accretion Disk Coronae

A new mechanism to form a magnetic pressure supported, high temperature corona above the photosphere of an accretion disk is explored using three dimensional radiation magneto-hydrodynamic (MHD) simulations. The thermal properties of the disk are calculated self-consistently by balancing radiative cooling through the surfaces of the disk with heating due to dissipation of turbulence driven by magneto-rotational instability (MRI). As has been noted in previous work, we find the dissipation rate per unit mass increases dramatically with height above the mid-plane, in stark contrast to the alpha-disk model which assumes this quantity is a constant. Thus, we find that in simulations with a low surface density (and therefore a shallow photosphere), the fraction of energy dissipated above the photosphere is significant (about 3.4% in our lowest surface density model), and this fraction increases as surface density decreases. When a significant fraction of the accretion energy is dissipated in the optically thin photosphere, the gas temperature increases substantially and a high temperature, magnetic pressure supported corona is formed. The volume-averaged temperature in the disk corona is more than 10 times larger than at the disk mid-plane. Moreover, gas temperature in the corona is strongly anti-correlated with gas density, which implies the corona formed by MRI turbulence is patchy. This mechanism to form an accretion disk corona may help explain the observed relation between the spectral index and luminosity from AGNs, and the soft X-ray excess from some AGNs. It may also be relevant to spectral state changes in X-ray binaries.

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.

preprint2013arXiv

Low mass planets in protoplanetary disks with net vertical magnetic fields: the Planetary Wake and Gap Opening

We study wakes and gap opening by low mass planets in gaseous protoplanetary disks threaded by net vertical magnetic fields which drive magnetohydrodynamical (MHD) turbulence through the magnetorotational instabilty (MRI), using three dimensional simulations in the unstratified local shearing box approximation. The wakes, which are excited by the planets, are damped by shocks similar to the wake damping in inviscid hydrodynamic (HD) disks. Angular momentum deposition by shock damping opens gaps in both MHD turbulent disks and inviscid HD disks even for low mass planets, in contradiction to the "thermal criterion" for gap opening. To test the "viscous criterion", we compared gap properties in MRI-turbulent disks to those in viscous HD disks having the same stress, and found that the same mass planet opens a significantly deeper and wider gap in net vertical flux MHD disks than in viscous HD disks. This difference arises due to the efficient magnetic field transport into the gap region in MRI disks, leading to a larger effective αwithin the gap. Thus, across the gap, the Maxwell stress profile is smoother than the gap density profile, and a deeper gap is needed for the Maxwell stress gradient to balance the planetary torque density. We also confirmed the large excess torque close to the planet in MHD disks, and found that long-lived density features (termed zonal flows) produced by the MRI can affect planet migration. The comparison with previous results from net toroidal flux/zero flux MHD simulations indicates that the magnetic field geometry plays an important role in the gap opening process. Overall, our results suggest that gaps can be commonly produced by low mass planets in realistic protoplanetary disks, and caution the use of a constant α-viscosity to model gaps in protoplanetary disks.

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.

preprint2013arXiv

Saturation of the MRI in Strongly Radiation Dominated Accretion Disks

The saturation level of the magneto-rotational instability (MRI) in a strongly radiation dominated accretion disk is studied using a new Godunov radiation MHD code in the unstratified shearing box approximation. Since vertical gravity is neglected in this work, our focus is on how the MRI saturates in the optically thick mid-plane of the disk. We confirm that turbulence generated by the MRI is very compressible in the radiation dominated regime, as found by previous calculations using the flux-limited diffusion approximation. We also find little difference in the saturation properties in calculations that use a larger horizontal domain (up to four times the vertical scale height in the radial direction). However, in strongly radiation pressure dominated disks (one in which the radiation energy density reaches 1% of the rest mass energy density of the gas), we find Maxwell stress from the MRI turbulence is larger than the value produced when radiation pressure is replaced with the same amount of gas pressure. At the same time, the ratio between Maxwell stress and Reynolds stress is increased by almost a factor of 8 compared with the gas pressure dominated case. We suggest that this effect is caused by radiation drag, which acts like bulk viscosity and changes the effective magnetic Prandtl number of the fluid. Radiation viscosity significantly exceeds both the microscopic plasma viscosity and resistivity, ensuring that radiation dominated systems occupy the high magnetic Prandtl number regime. Nevertheless, we find radiative shear viscosity is negligible compared to the Maxwell and Reynolds stress in the flow. This may have important implications for the structure of radiation dominated accretion disks.

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.

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.

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.

preprint2012arXiv

A Godunov Method for Multidimensional Radiation Magnetohydrodynamics based on a variable Eddington tensor

We describe a numerical algorithm to integrate the equations of radiation magnetohydrodynamics in multidimensions using Godunov methods. This algorithm solves the radiation moment equations in the mixed frame, without invoking any diffusion-like approximations. The moment equations are closed using a variable Eddington tensor whose components are calculated from a formal solution of the transfer equation at a large number of angles using the method of short characteristics. We use a comprehensive test suite to verify the algorithm, including convergence tests of radiation-modified linear acoustic and magnetosonic waves, the structure of radiation modified shocks, and two-dimensional tests of photon bubble instability and the ablation of dense clouds by an intense radiation field. These tests cover a very wide range of regimes, including both optically thick and thin flows, and ratios of the radiation to gas pressure of at least 10^{-4} to 10^{4}. Across most of the parameter space, we find the method is accurate. However, the tests also reveal there are regimes where the method needs improvement, for example when both the radiation pressure and absorption opacity are very large. We suggest modifications to the algorithm that will improve accuracy in this case. We discuss the advantages of this method over those based on flux-limited diffusion. In particular, we find the method is not only substantially more accurate, but often no more expensive than the diffusion approximation for our intended applications.

preprint2012arXiv

A New Way to Conserve Total Energy for Eulerian Hydrodynamic Simulations with Self-Gravity

We propose a new method to conserve the total energy to round-off error in grid-based codes for hydrodynamic simulations with self-gravity. A formula for the energy flux due to the work done by the the self-gravitational force is given, so the change in total energy can be written in conservative form. Numerical experiments with the code Athena show that the total energy is indeed conserved with our new algorithm and the new algorithm is second order accurate. We have performed a set of tests that show the numerical errors in the traditional, non-conservative algorithm can affect the dynamics of the system. The new algorithm only requires one extra solution of the Poisson equation, as compared to the traditional algorithm which includes self-gravity as a source term. If the Poisson solver takes a negligible fraction of the total simulation time, such as when FFTs are used, the new algorithm is almost as efficient as the original method. This new algorithm is useful in Eulerian hydrodynamic simulations with self-gravity, especially when results are sensitive to small energy errors, as for radiation pressure dominated flow.

preprint2012arXiv

A Radiation Transfer Solver for Athena using Short Characteristics

We describe the implementation of a module for the Athena magnetohydrodynamics (MHD) code which solves the time-independent, multi-frequency radiative transfer (RT) equation on multidimensional Cartesian simulation domains, including scattering and non-LTE effects. The module is based on well-known and well-tested algorithms developed for modeling stellar atmospheres, including the method of short characteristics to solve the RT equation, accelerated Lambda iteration to handle scattering and non-LTE effects, and parallelization via domain decomposition. The module serves several purposes: it can be used to generate spectra and images, to compute a variable Eddington tensor (VET) for full radiation MHD simulations, and to calculate the heating and cooling source terms in the MHD equations in flows where radiation pressure is small compared with gas pressure. For the latter case, the module is combined with the standard MHD integrators using operator-splitting and we describe this approach in detail. Implementation of the VET method for radiation pressure dominated flows is described in a companion paper. We present results from a suite of test problems for both the RT solver itself, and for dynamical problems that include radiative heating and cooling. These tests demonstrate that the radiative transfer solution is accurate, and confirm that the operator split method is stable, convergent, and efficient for problems of interest. We demonstrate there is no need to adopt ad-hoc assumptions of questionable accuracy to solve RT problems in concert with MHD: the computational cost for our general-purpose module for simple (e.g. LTE grey) problems can be comparable to or less than a single timestep of Athena's MHD integrators, and only few times more expensive than that for more general problems. (Abridged)

preprint2012arXiv

Angular Momentum Transport and Variability in Boundary Layers of Accretion Disks Driven by Global Acoustic Modes

Disk accretion onto a weakly magnetized central object, e.g. a star, is inevitably accompanied by the formation of a boundary layer near the surface, in which matter slows down from the highly supersonic orbital velocity of the disk to the rotational velocity of the star. We perform high resolution 2D hydrodynamical simulations in the equatorial plane of an astrophysical boundary layer with the goal of exploring the dynamics of non-axisymmetric structures that form there. We generically find that the supersonic shear in the boundary layer excites non-axisymmetric quasi-stationary acoustic modes that are trapped between the surface of the star and a Lindblad resonance in the disk. These modes rotate in a prograde fashion, are stable for hundreds of orbital periods, and have a pattern speed that is less than and of order the rotational velocity at the inner edge of the disk. The origin of these intrinsically global modes is intimately related to the operation of a corotation amplifier in the system. Dissipation of acoustic modes in weak shocks provides a universal mechanism for angular momentum and mass transport even in purely hydrodynamic (i.e. non-magnetized) boundary layers. We discuss the possible implications of these trapped modes for explaining the variability seen in accreting compact objects.

preprint2012arXiv

Buoyancy Instabilities in a Weakly Collisional Intracluster Medium

The intracluster medium of galaxy clusters is a weakly collisional, high-beta plasma in which the transport of heat and momentum occurs primarily along magnetic-field lines. Anisotropic heat conduction allows convective instabilities to be driven by temperature gradients of either sign, the magnetothermal instability (MTI) in the outskirts of non-isothermal clusters and the heat-flux buoyancy-driven instability (HBI) in their cooling cores. We employ the Athena MHD code to investigate the nonlinear evolution of these instabilities, self-consistently including the effects of anisotropic viscosity (i.e. Braginskii pressure anisotropy), anisotropic conduction, and radiative cooling. We highlight the importance of the microscale instabilities that inevitably accompany and regulate the pressure anisotropies generated by the HBI and MTI. We find that, in all but the innermost regions of cool-core clusters, anisotropic viscosity significantly impairs the ability of the HBI to reorient magnetic-field lines orthogonal to the temperature gradient. Thus, while radio-mode feedback appears necessary in the central few tens of kpc, conduction may be capable of offsetting radiative losses throughout most of a cool core over a significant fraction of the Hubble time. Magnetically-aligned cold filaments are then able to form by local thermal instability. Viscous dissipation during the formation of a cold filament produces accompanying hot filaments, which can be searched for in deep Chandra observations of nearby cool-core clusters. In the case of the MTI, anisotropic viscosity maintains the coherence of magnetic-field lines over larger distances than in the inviscid case, providing a natural lower limit for the scale on which the field can fluctuate freely. In the nonlinear state, the magnetic field exhibits a folded structure in which the field-line curvature and field strength are anti-correlated.

preprint2012arXiv

Planet-Disk interaction in 3D: the importance of buoyancy waves

We carry out local three dimensional (3D) hydrodynamic simulations of planet-disk interaction in stratified disks with varied thermodynamic properties. We find that whenever the Brunt-Vaisala frequency (N) in the disk is nonzero, the planet exerts a strong torque on the disk in the vicinity of the planet, with a reduction in the traditional "torque cutoff". In particular, this is true for adiabatic perturbations in disks with isothermal density structure, as should be typical for centrally irradiated protoplanetary disks. We identify this torque with buoyancy waves, which are excited (when N is non-zero) close to the planet, within one disk scale height from its orbit. These waves give rise to density perturbations with a characteristic 3D spatial pattern which is in close agreement with the linear dispersion relation for buoyancy waves. The torque due to these waves can amount to as much as several tens of per cent of the total planetary torque, which is not expected based on analytical calculations limited to axisymmetric or low-m modes. Buoyancy waves should be ubiquitous around planets in the inner, dense regions of protoplanetary disks, where they might possibly affect planet migration.

preprint2012arXiv

Two-Dimensional Magnetohydrodynamic Simulations of Barred Galaxies

Barred galaxies are known to possess magnetic fields that may affect the properties of bar substructures such as dust lanes and nuclear rings. We use two-dimensional high-resolution magnetohydrodynamic (MHD) simulations to investigate the effects of magnetic fields on the formation and evolution of such substructures as well as on the mass inflow rates to the galaxy center. The gaseous medium is assumed to be infinitesimally-thin, isothermal, non-self-gravitating, and threaded by initially uniform, azimuthal magnetic fields. We find that there exists an outermost x1-orbit relative to which gaseous responses to an imposed stellar bar potential are completely different between inside and outside. Inside this orbit, gas is shocked into dust lanes and infalls to form a nuclear ring. Magnetic fields are compressed in dust lanes, reducing their peak density. Magnetic stress removes further angular momentum of the gas at the shocks, temporarily causing the dust lanes to bend into an 'L' shape and eventually leading to a smaller and more centrally distributed ring than in unmagnetized models. The mass inflow rates in magnetized models correspondingly become larger, by more than two orders of magnitude when the initial fields have an equipartition value with thermal energy, than in the unmagnetized counterparts. Outside the outermost x1-orbit, on the other hand, an MHD dynamo due to the combined action of the bar potential and background shear operates near the corotation and bar-end regions, efficiently amplifying magnetic fields. The amplified fields shape into trailing magnetic arms with strong fields and low density. The base of the magnetic arms has a thin layer in which magnetic fields with opposite polarity reconnect via a tearing-mode instability. This produces numerous magnetic islands with large density which propagate along the arms to turn the outer disk into a highly chaotic state.

preprint2011arXiv

A Second Order Godunov Method for Multidimensional Relativistic Magnetohydrodynamics

We describe a new Godunov algorithm for relativistic magnetohydrodynamics (RMHD) that combines a simple, unsplit second order accurate integrator with the constrained transport (CT) method for enforcing the solenoidal constraint on the magnetic field. A variety of approximate Riemann solvers are implemented to compute the fluxes of the conserved variables. The methods are tested with a comprehensive suite of multidimensional problems. These tests have helped us develop a hierarchy of correction steps that are applied when the integration algorithm predicts unphysical states due to errors in the fluxes, or errors in the inversion between conserved and primitive variables. Although used exceedingly rarely, these corrections dramatically improve the stability of the algorithm. We present preliminary results from the application of these algorithms to two problems in RMHD: the propagation of supersonic magnetized jets, and the amplification of magnetic field by turbulence driven by the relativistic Kelvin-Helmholtz instability (KHI). Both of these applications reveal important differences between the results computed with Riemann solvers that adopt different approximations for the fluxes. For example, we show that use of Riemann solvers which include both contact and rotational discontinuities can increase the strength of the magnetic field within the cocoon by a factor of ten in simulations of RMHD jets, and can increase the spectral resolution of three-dimensional RMHD turbulence driven by the KHI by a factor of 2. This increase in accuracy far outweighs the associated increase in computational cost. Our RMHD scheme is publicly available as part of the Athena code.

preprint2011arXiv

Central Regions of Barred Galaxies: Two-Dimensional Non-self-gravitating Hydrodynamic Simulations

The inner regions of barred galaxies contain substructures such as off-axis shocks, nuclear rings, and nuclear spirals. These substructure may affect star formation, and control the activity of a central black hole (BH) by determining the mass inflow rate. We investigate the formation and properties of such substructures using high-resolution, grid-based hydrodynamic simulations. The gaseous medium is assumed to be infinitesimally-thin, isothermal, and non-self-gravitating. The stars and dark matter are represented by a static gravitational potential with four components: a stellar disk, the bulge, a central BH, and the bar. To investigate various galactic environments, we vary the gas sound speed c_s as well as the mass of the central BH M_BH. Once the flow has reached a quasi-steady state, off-axis shocks tend to move closer to the bar major axis as c_s increases. Nuclear rings shrink in size with increasing c_s, but are independent of M_BH, suggesting that ring position is not determined by the Lindblad resonances. Rings in low-c_s models are narrow since they are occupied largely by gas on x2-orbits and well decoupled from nuclear spirals, while they become broad because of large thermal perturbations in high-c_s models. Nuclear spirals persist only when either c_s is small or M_BH is large; they would otherwise be destroyed completely by the ring material on eccentric orbits. The shape and strength of nuclear spirals depend sensitively on c_s and M_BH such that they are leading if both c_s and M_BH are small, weak trailing if c_s is small and M_BH is large, and strong trailing if both c_s and M_BH are large. While the mass inflow rate toward the nucleus is quite small in low-c_s models because of the presence of a narrow nuclear ring, it becomes larger than 0.01 Msun/yr when c_s is large, providing a potential explanation of nuclear activity in Seyfert galaxies.

preprint2011arXiv

Density Waves Excited by Low-Mass Planets in Protoplanetary Disks I: Linear Regime

Density waves excited by planets embedded in protoplanetary disks play a central role in planetary migration and gap opening processes. We carry out 2D shearing sheet simulations to study the linear regime of wave evolution with the grid-based code Athena, and provide detailed comparisons with the theoretical predictions. Low mass planets (down to ~0.03 Earth mass at 1 AU) and high spatial resolution (256 grid points per scale height) are chosen to mitigate the effects of wave nonlinearity. To complement the existing numerical studies, we focus on the primary physical variables such as the spatial profile of the wave, torque density, and the angular momentum flux carried by the wave, instead of secondary quantities such as the planetary migration rate. Our results show percent level agreement with theory in both physical and Fourier space. New phenomena such as the change of the toque density sign far from the planet are discovered and discussed. Also, we explore the effect of the numerical algorithms, and find that a high order of accuracy, high resolution, and an accurate planetary potential are crucial to achieve good agreement with the theory. We find that the use of a too large time-step without properly resolving the dynamical time scale around the planet produces incorrect results, and may lead to spurious gap opening. Global simulations of planet migration and gap opening violating this requirement may be affected by spurious effects resulting in e.g. the incorrect planetary migration rate and gap opening mass.

preprint2011arXiv

Density Waves Excited by Low-Mass Planets in Protoplanetary Disks II: High-Resolution Simulations of the Nonlinear Regime

We investigate numerically the propagation of density waves excited by a low-mass planet in a protoplanetary disk in the nonlinear regime, using 2D local shearing box simulations with the grid-based code Athena at high spatial resolution (256 grid points per scale height h). The nonlinear evolution results in the wave steepening into a shock, causing damping and angular momentum transfer to the disk. On long timescales this leads to spatial redistribution of the disk density, causing migration feedback and potentially resulting in gap opening. Previous numerical studies concentrated on exploring these secondary phenomena as probes of the nonlinear wave evolution. Here we focus on exploring the evolution of the basic wave properties, such as its density profile evolution, shock formation, post-shock wave behavior, and provide comparison with analytical theory. The generation of potential vorticity at the shock is computed analytically and is subsequently verified by simulations and used to pinpoint the shock location. We confirm the theoretical relation between the shocking length and the planet mass (including the effect of the equation of state), and the post-shock decay of the angular momentum flux carried by the wave. The post-shock evolution of the wave profile is explored, and we quantitatively confirm its convergence to the theoretically expected N-wave shape. The accuracy of various numerical algorithms used to compute the nonlinear wave evolution is also investigated: we find that higher order spatial reconstruction and high resolution are crucial for capturing the shock formation correctly.

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.

preprint2011arXiv

Global Simulations of Accretion Disks I: Convergence and Comparisons with Local Models

Grid-based magnetohydrodynamic (MHD) simulations have proven invaluable for the study of astrophysical accretion disks. However, the fact that angular momentum transport in disks is mediated by MHD turbulence (with structure down to very small scales) raises the concern that the properties of the modeled accretion disks are affected by the finite numerical resolution of the simulation. By implementing an orbital advection algorithm into the Athena code in cylindrical geometry, we have performed a set of global (but unstratified) Newtonian disk simulations extending up to resolutions previously unattained. We study the convergence of these models as a function of spatial resolution and initial magnetic field geometry. The usual viscosity parameter ($α$) or the ratio of thermal-to-magnetic pressure ($β$) are found to be poor diagnostics of convergence, whereas the average tilt angle of the magnetic field in the $(r,ϕ)$-plane is a very good diagnostic of convergence. We suggest that this is related to the saturation of the MHD turbulence via parasitic modes of the magnetorotational instability. Even in the case of zero-net magnetic flux, we conclude that our highest resolution simulations (with 32-zones and 64-zones per vertical scale height) have achieved convergence. Our global simulations reach resolutions comparable to those used in local, shearing box models of MHD disk turbulence. We find that the saturation predictors derived from local simulations correspond well to the instantaneous correlations between local flux and stress found in our global simulations. However, the conservation of magnetic flux implicit in local models is not realized in our global disks. Thus, the magnetic connectivity of an accretion disk represents physics that is truly global and cannot be captured in any ab-initio local model.

preprint2011arXiv

The Effect of Anisotropic Conduction on the Thermal Instability in the Interstellar Medium

Thermal instability (TI) can strongly affect the structure and dynamics of the interstellar medium (ISM) in the Milky Way and other disk galaxies. Thermal conduction plays an important role in the TI by stabilizing small scales and limiting the size of the smallest condensates. In the magnetized ISM, however, heat is conducted anisotropically (primarily along magnetic field lines). We investigate the effects of anisotropic thermal conduction on the nonlinear regime of the TI by performing two-dimensional magnetohydrodynamic simulations. We present models with magnetic fields of different initial geometries and strengths, and compare them to hydrodynamic models with isotropic conduction. We find anisotropic conduction does not significantly alter the overall density and temperature statistics in the saturated state of the TI. However, it can strongly affect the shapes and sizes of cold clouds formed by the TI. For example, for uniform initial fields long filaments of cold gas are produced that are reminiscent of some observed HI clouds. For initially tangled fields, such filaments are not produced. We also show that anisotropic conduction suppresses turbulence generated by evaporative flows from the surfaces of cold blobs, which may have implications for mechanisms for driving turbulence in the ISM.

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.

preprint2010arXiv

Implementation of the Shearing Box Approximation in Athena

We describe the implementation of the shearing box approximation for the study of the dynamics of accretion disks in the Athena magnetohydrodynamics (MHD) code. Second-order Crank-Nicholson time differencing is used for the Coriolis and tidal gravity source terms that appear in the momentum equation for accuracy and stability. We show this approach conserves energy for epicyclic oscillations in hydrodynamic flows to round-off error. In the energy equation, the tidal gravity source terms are differenced as the gradient of an effective potential in a way which guarantees that total energy (including the gravitational potential energy) is also conserved to round-off error. We introduce an orbital advection algorithm for MHD based on constrained transport to preserve the divergence-free constraint on the magnetic field. This algorithm removes the orbital velocity from the time step constraint, and makes the truncation error more uniform in radial position. Modifications to the shearing box boundary conditions applied at the radial boundaries are necessary to conserve the total vertical magnetic flux. In principle similar corrections are also required to conserve mass, momentum and energy, however in practice we find the orbital advection method conserves these quantities to better than 0.03% over hundreds of orbits. The algorithms have been applied to studies of the nonlinear regime of the MRI in very wide (up to 32 scale heights) horizontal domains.

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

preprint2009arXiv

Buoyant Bubbles in Intracluster Gas: Effects of Magnetic Fields and Anisotropic Viscosity

Recent observations by Chandra and XMM-Newton indicate there are complex structures at the cores of galaxy clusters, such as cavities and filaments. One plausible model for the formation of such structures is the interaction of radio jets with the intracluster medium (ICM). To investigate this idea, we use three-dimensional magnetohydrodynamic simulations including anisotropic (Braginskii) viscosity to study the effect of magnetic fields on the evolution and morphology of buoyant bubbles in the ICM. We investigate a range of different initial magnetic field geometries and strengths, and study the resulting x-ray surface brightness distribution for comparison to observed clusters. Magnetic tension forces and viscous transport along field lines tend to suppress instabilities parallel, but not perpendicular, to field lines. Thus, the evolution of the bubble depends strongly on the initial field geometry. We find toroidal field loops initially confined to the interior of the bubble are best able reproduce the observed cavity structures.

preprint2009arXiv

Turbulent resistivity driven by the magnetorotational instability

We measure the turbulent resistivity in the nonlinear regime of the MRI, and evaluate the turbulent magnetic Prandtl number. We perform a set of numerical simulations with the Eulerian finite volume codes Athena and Ramses in the framework of the shearing box model. We consider models including explicit dissipation coefficients and magnetic field topologies such that the net magnetic flux threading the box in both the vertical and azimuthal directions vanishes. We first demonstrate good agreement between the two codes by comparing the properties of the turbulent states in simulations having identical microscopic diffusion coefficients (viscosity and resistivity). We find the properties of the turbulence do not change when the box size is increased in the radial direction, provided it is elongated in the azimuthal direction. To measure the turbulent resistivity in the disk, we impose a fixed electromotive force on the flow and measure the amplitude of the saturated magnetic field that results. We obtain a turbulent resistivity that is in rough agreement with mean field theories like the Second Order Smoothing Approximation. The numerical value translates into a turbulent magnetic Prandtl number Pm_t of order unity. Pm_t appears to be an increasing function of the forcing we impose. It also becomes smaller as the box size is increased in the radial direction, in good agreement with previous results obtained in very large boxes. Our results are in general agreement with other recently published papers studying the same problem but using different methodology. Thus, our conclusion that Pm_t is of order unity appears robust.

preprint2008arXiv

Anisotropic winds from close-in extra-solar planets

We present two-dimensional hydrodynamic models of thermally driven winds from highly irradiated, close-in extra-solar planets. We adopt a very simple treatment of the radiative heating processes at the base of the wind, and instead focus on the differences between the properties of outflows in multidimensions in comparison to spherically symmetric models computed with the same methods. For hot (T > 2 x 10^{4} K) or highly ionized gas, we find strong (supersonic) polar flows are formed above the planet surface which produce weak shocks and outflow on the night-side. In comparison to a spherically symmetric wind with the same parameters, the sonic surface on the day-side is much closer to the planet surface in multidimensions, and the total mass loss rate is reduced by almost a factor of four. We also compute the steady-state structure of interacting planetary and stellar winds. Both winds end in a termination shock, with a parabolic contact discontinuity which is draped over the planet separating the two shocked winds. The planetary wind termination shock and the sonic surface in the wind are well separated, so that the mass loss rate from the planet is essentially unaffected. However, the confinement of the planetary wind to the small volume bounded by the contact discontinuity greatly enhances the column density close to the planet, which might be important for the interpretation of observations of absorption lines formed by gas surrounding transiting planets.

preprint2008arXiv

Dissipation and Heating in Supersonic Hydrodynamic and MHD Turbulence

We study energy dissipation and heating by supersonic MHD turbulence in molecular clouds using Athena, a new higher-order Godunov code. We analyze the dependence of the saturation amplitude, energy dissipation characteristics, power spectra, sonic scaling, and indicators of intermittency in the turbulence on factors such as the magnetic field strength, driving scale, energy injection rate, and numerical resolution. While convergence in the energies is reached at moderate resolutions, we find that the power spectra require much higher resolutions that are difficult to obtain. In a 1024^3 hydro run, we find a power law relationship between the velocity dispersion and the spatial scale on which it is measured, while for an MHD run at the same resolution we find no such power law. The time-variability and temperature intermittency in the turbulence both show a dependence on the driving scale, indicating that numerically driving turbulence by an arbitrary mechanism may not allow a realistic representation of these properties. We also note similar features in the power spectrum of the compressive component of velocity for supersonic MHD turbulence as in the velocity spectrum of an initially-spherical MHD blast wave, implying that the power law form does not rule out shocks, rather than a turbulent cascade, playing a significant role in the regulation of energy transfer between spatial scales.

preprint2008arXiv

Spherical Accretion with Anisotropic Thermal Conduction

We study the effects of anisotropic thermal conduction on magnetized spherical accretion flows using global axisymmetric MHD simulations. In low collisionality plasmas, the Bondi spherical accretion solution is unstable to the magnetothermal instability (MTI). The MTI grows rapidly at large radii where the inflow is subsonic. For a weak initial field, the MTI saturates by creating a primarily radial magnetic field, i.e., by aligning the field lines with the background temperature gradient. The saturation is quasilinear in the sense that the magnetic field is amplified by a factor of $\sim 10-30$ independent of the initial field strength (for weak fields). In the saturated state, the conductive heat flux is much larger than the convective heat flux, and is comparable to the field-free (Spitzer) value (since the field lines are largely radial). The MTI by itself does not appreciably change the accretion rate $\dot M$ relative to the Bondi rate $\dot M_B$. However, the radial field lines created by the MTI are amplified by flux freezing as the plasma flows in to small radii. Oppositely directed field lines are brought together by the converging inflow, leading to significant resistive heating. When the magnetic energy density is comparable to the gravitational potential energy density, the plasma is heated to roughly the virial temperature; the mean inflow is highly subsonic; most of the energy released by accretion is transported to large radii by thermal conduction; and the accretion rate $\dot M \ll \dot M_B$. The predominantly radial magnetic field created by the MTI at large radii in spherical accretion flows may account for the stable Faraday rotation measure towards Sgr A* in the Galactic Center.

preprint2007arXiv

Effect of the Coriolis Force on the Hydrodynamics of Colliding Wind Binaries

Using fully three-dimensional hydrodynamic simulations, we investigate the effect of the Coriolis force on the hydrodynamic and observable properties of colliding wind binary systems. To make the calculations tractable, we assume adiabatic, constant velocity winds. The neglect of radiative driving, gravitational deceleration, and cooling limit the application of our models to real systems. However, these assumptions allow us to isolate the effect of the Coriolis force, and by simplifying the calculations, allow us to use a higher resolution (up to 640^3) and to conduct a larger survey of parameter space. We study the dynamics of collidng winds with equal mass loss rates and velocities emanating from equal mass stars on circular orbits, with a range of values for the ratio of the wind to orbital velocity. We also study the dynamics of winds from stars on elliptical orbits and with unequal strength winds. Orbital motion of the stars sweeps the shocked wind gas into an Archimedean spiral, with asymmetric shock strengths and therefore unequal postshock temperatures and densities in the leading and trailing edges of the spiral. We observe the Kelvin-Helmholtz instability at the contact surface between the shocked winds in systems with orbital motion even when the winds are identical. The change in shock strengths caused by orbital motion increases the volume of X-ray emitting post-shock gas with T > 0.59 keV by 63% for a typical system as the ratio of wind velocity to orbital velocity decreases to V_w/V_o = 2.5. This causes increased free-free emission from systems with shorter orbital periods and an altered time-dependence of the wind attenuation. We comment on the importance of the effects of orbital motion on the observable properties of colliding wind binaries.

preprint2007arXiv

Electron Heating in Hot Accretion Flows

Local (shearing box) simulations of the nonlinear evolution of the magnetorotational instability in a collisionless plasma show that angular momentum transport by pressure anisotropy ($p_\perp \ne p_\parallel$, where the directions are defined with respect to the local magnetic field) is comparable to that due to the Maxwell and Reynolds stresses. Pressure anisotropy, which is effectively a large-scale viscosity, arises because of adiabatic invariants related to $p_\perp$ and $p_\parallel$ in a fluctuating magnetic field. In a collisionless plasma, the magnitude of the pressure anisotropy, and thus the viscosity, is determined by kinetic instabilities at the cyclotron frequency. Our simulations show that $\sim 50$ % of the gravitational potential energy is directly converted into heat at large scales by the viscous stress (the remaining energy is lost to grid-scale numerical dissipation of kinetic and magnetic energy). We show that electrons receive a significant fraction ($\sim [T_e/T_i]^{1/2}$) of this dissipated energy. Employing this heating by an anisotropic viscous stress in one dimensional models of radiatively inefficient accretion flows, we find that the radiative efficiency of the flow is greater than 0.5% for $\dot{M} \gtrsim 10^{-4} \dot{M}_{Edd}$. Thus a low accretion rate, rather than just a low radiative efficiency, is necessary to explain the low luminosity of many accreting black holes. For Sgr A* in the Galactic Center, our predicted radiative efficiencies imply an accretion rate of $\approx 3 \times 10^{-8} M_\odot {\rm yr^{-1}}$ and an electron temperature of $\approx 3 \times 10^{10}$ K at $\approx 10$ Schwarzschild radii; the latter is consistent with the brightness temperature inferred from VLBI observations.

preprint2007arXiv

Faraday Rotation in Global Accretion Disk Simulations: Implications for Sgr A*

These Faraday rotation calculations of hot, thick accretion flows are motivated by the measured steady rotation measure (RM) of $\approx -6 \times 10^5$ rad m$^{-2}$ from Sgr A*. In our numerical simulations, the quasi-steady state structure of the accretion flow, and the RM it produces, depends on the initial magnetic field. In spite of this dependence, we can draw several robust conclusions about Faraday rotation produced by geometrically thick accretion disks: i) the time averaged RM does not depend that sensitively on the viewing angle, but the stability of the RM can. Equatorial viewing angles show significant variability in RM (including sign reversals), while polar viewing angles are relatively stable if there is a large scale magnetic field threading the disk at large radii. ii) Most of the RM is produced at small radii for polar viewing angles while all radii contribute significantly near the midplane of the disk. Our simulations confirm previous analytic arguments that the accretion rate onto Sgr A* must satisfy $\dot M_{\rm in} \ll \dot M_{\rm Bondi} \sim 10^{-5} \mpy$ in order to not over-produce the measured RM. We argue that the steady RM $\approx -6 \times 10^5$ rad m$^{-2}$ from Sgr A* has two plausible explanations: 1) it is produced at $\sim 100$ Schwarzschild radii, requires $\dot{M}_{\rm in} \approx 3 \times 10^{-8} M_\odot$ yr$^{-1}$, and we view the flow at an angle of $\sim 30^\circ$ relative to the rotation axis of the disk; in our simulations, the variation in RM across a finite-sized source is sufficient to depolarize the emission below $\approx$ 100 GHz, consistent with observations. 2) Alternatively, the RM may be produced in the relatively spherical inflowing plasma near the circularization radius at $\sim 10^3-10^4$ Schwarzschild radii.

preprint2006arXiv

Three Dimensional Compressible Hydrodynamic Simulations of Vortices in Disks

We carry out three-dimensional, high resolution (up to $1024^2\times 256$) hydrodynamic simulations of the evolution of vortices in vertically unstratified Keplerian disks using the shearing sheet approximation. The transient amplification of incompressible, linear amplitude leading waves (which has been proposed as a possible route to nonlinear hydrodynamical turbulence in disks) is used as one test of our algorithms; our methods accurately capture the predicted amplification, converges at second-order, and is free from aliasing. Waves expected to reach nonlinear amplitude at peak amplification become unstable to Kelvin-Helmholtz modes when $\mid W_{\rm max}\mid\gtrsim Ω$ (where $W_{\rm max}$ is the local maximum of vorticity and $Ω$ the angular velocity). We study the evolution of a power-law distribution of vorticity consistent with Kolmogorov turbulence; in two-dimensions long-lived vortices emerge and decay slowly, similar to previous studies. In three-dimensions, however, vortices are unstable to bending modes, leading to rapid decay. Only vortices with a length to width ratio smaller than one survive; in three-dimensions the residual kinetic energy and shear stress is at least one order of magnitude smaller than in two-dimensions. No evidence for sustained hydrodynamical turbulence and transport is observed in three-dimensions. Instead, at late times the residual transport is determined by the amplitude of slowly decaying, large-scale vortices (with horizontal extent comparable to the scale height of the disk), with additional contributions from nearly incompressible inertial waves possible. Evaluating the role that large-scale vortices play in astrophysical accretion disks will require understanding the mechanisms that generate and destroy them.

preprint2005arXiv

Diffusion coefficient of a passive contaminant in a local MHD model of a turbulent accretion disc

We calculate the radial diffusion coefficient for a passive contaminant in an accretion disc which is turbulent due to the action of the magnetorotational instability. Numerical MHD simulations are used to follow the evolution of a local patch of the disc using the shearing box formalism. A separate continuity equation for the mass fraction of contaminant is integrated along with the MHD system, and radial profiles of this fraction are obtained as a function of time. Solutions of a linear diffusion equation are fitted to the numerical measured profiles of the contaminant, treating the diffusion coefficient D as the fitting parameter. At early times, the value of D is found to vary, however once the contaminant is spread over scales comparable to the box size, it saturates at a steady value. The ratio of D to the transport coefficient of angular momentum due to shear stress is small. If D can be used as a proxy for the turbulent magnetic diffusivity, the effective magnetic Prandtl number P_eff=ν/D (where νis the coefficient of "efective viscosity" associated with shear stress) would be large.

preprint1997arXiv

Radiation driven winds from luminous accretion disks

We study the 2-D, time-dependent hydrodynamics of radiation-driven winds from accretion disks in which the radiation force is mediated by spectral lines. If the dominant contribution to the total radiation field comes from the disk, then we find the outflow is intrinsically unsteady and characterised by large amplitude velocity and density changes. Both infall and outflow can occur in different regions of the wind at the same time. On the other hand, if the total luminosity of the system is dominated by the central star, then the outflow is steady. In either case, we find the 2-D structure of the wind consists of a dense, slow outflow, typically confined to angles within about 45$^o$ of the equatorial plane, that is bounded on the polar side by a high -velocity, lower density stream. The flow geometry is controlled largely by the geometry of the radiation field. Global properties e.g., the total mass loss rate and terminal velocity depend more on the system luminosity and are insensitive to geometry. Matter is fed into the fast wind from within a few stellar radii of the central star. Our solutions agree qualitatively with the kinematics of outflows in CV systems inferred from spectroscopic observations. We predict that low luminosity systems may display unsteady behavior in wind-formed spectral lines. Our study also has application to winds from active galactic nuclei and from high mass YSOs.

James M. Stone

What is connected

Connect this record

See the researcher in context

Building this map preview

59 published item(s)

Dissipation of AGN jets in a clumpy interstellar medium

Ab Initio Horizon-Scale Simulations of Magnetically Arrested Accretion in Sagittarius A* Fed by Stellar Winds

Global 3-D Radiation Magnetohydrodynamic Simulations for FU Ori's Accretion Disk and Observational Signatures of Magnetic Fields

The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers

The Surprisingly Small Impact of Magnetic Fields On The Inner Accretion Flow of Sagittarius A* Fueled By Stellar Winds

An Extension of the Athena++ Code Framework for GRMHD Based on Advanced Riemann Solvers and Staggered-Mesh Constrained Transport

Dust dynamics in 2D gravito-turbulent disks

Global MHD Simulations of Accretion Disks in Cataclysmic Variables (CVs): I. The Importance of Spiral Shocks

Magnetorotational Turbulence and Dynamo in a Collisionless Plasma

Shock-driven Accretion in Circumplanetary Disks: Observables and Satellite Formation

Saturation of the magnetorotational instability in the unstratified shearing box with zero net flux: convergence in taller boxes

Spreading layers in accreting objects: role of acoustic waves for angular momentum transport, mixing and thermodynamics

Survival and Structure of Dusty Vortices in Protoplanetary Discs

The Structure of Spiral Shocks Excited by Planetary-mass Companions

A Global Three Dimensional Radiation Magneto-hydrodynamic Simulation of Super-Eddington Accretion Disks

An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation

Characterizing the variation of propagation constants in multicore fibre

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

Dust Trapping by Vortices in Transitional Disks: Evidence for Non-ideal MHD Effects in Protoplanetary Disks

Firehose and Mirror Instabilities in a Collisionless Shearing Plasma

Magnetic Flux Concentration and Zonal Flows in Magnetorotational Instability Turbulence

Particle Concentration At Planet Induced Gap Edges and Vortices: I. Inviscid 3-D Hydro Disks

Radiation Feedback in ULIRGS: Are Photons Movers and Shakers?

Radiation Magneto-hydrodynamic Simulations of the Formation of Hot Accretion Disk Coronae

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

Low mass planets in protoplanetary disks with net vertical magnetic fields: the Planetary Wake and Gap Opening

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

Saturation of the MRI in Strongly Radiation Dominated Accretion Disks

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

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

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

A Godunov Method for Multidimensional Radiation Magnetohydrodynamics based on a variable Eddington tensor

A New Way to Conserve Total Energy for Eulerian Hydrodynamic Simulations with Self-Gravity

A Radiation Transfer Solver for Athena using Short Characteristics

Angular Momentum Transport and Variability in Boundary Layers of Accretion Disks Driven by Global Acoustic Modes

Buoyancy Instabilities in a Weakly Collisional Intracluster Medium

Planet-Disk interaction in 3D: the importance of buoyancy waves

Two-Dimensional Magnetohydrodynamic Simulations of Barred Galaxies

A Second Order Godunov Method for Multidimensional Relativistic Magnetohydrodynamics

Central Regions of Barred Galaxies: Two-Dimensional Non-self-gravitating Hydrodynamic Simulations

Density Waves Excited by Low-Mass Planets in Protoplanetary Disks I: Linear Regime

Density Waves Excited by Low-Mass Planets in Protoplanetary Disks II: High-Resolution Simulations of the Nonlinear Regime

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

Global Simulations of Accretion Disks I: Convergence and Comparisons with Local Models

The Effect of Anisotropic Conduction on the Thermal Instability in the Interstellar Medium

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

Implementation of the Shearing Box Approximation in Athena

Particle-Gas Dynamics with Athena: Method and Convergence

Buoyant Bubbles in Intracluster Gas: Effects of Magnetic Fields and Anisotropic Viscosity

Turbulent resistivity driven by the magnetorotational instability

Anisotropic winds from close-in extra-solar planets

Dissipation and Heating in Supersonic Hydrodynamic and MHD Turbulence

Spherical Accretion with Anisotropic Thermal Conduction

Effect of the Coriolis Force on the Hydrodynamics of Colliding Wind Binaries

Electron Heating in Hot Accretion Flows

Faraday Rotation in Global Accretion Disk Simulations: Implications for Sgr A*

Three Dimensional Compressible Hydrodynamic Simulations of Vortices in Disks

Diffusion coefficient of a passive contaminant in a local MHD model of a turbulent accretion disc

Radiation driven winds from luminous accretion disks