Jacob B. Simon

ATHENA: A NEW CODE FOR ASTROPHYSICAL MHD

A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher order Godunov methods with the constrained transport (CT) technique to enforce the divergence-free constraint on the magnetic field. Discretization is based on cell-centered volume averages for mass, momentum, and energy, and face-centered area averages for the magnetic field. Novel features of the algorithm include (1) a consistent framework for computing the time- and edge-averaged electric fields used by CT to evolve the magnetic field from the time- and area-averaged Godunov fluxes, (2) the extension to MHD of spatial reconstruction schemes that involve a dimensionally split time advance, and (3) the extension to MHD of two different dimensionally unsplit integration methods. Implementation of the algorithm in both C and FORTRAN95 is detailed, including strategies for parallelization using domain decomposition. Results from a test suite which includes problems in one-, two-, and three-dimensions for both hydrodynamics and MHD are given, not only to demonstrate the fidelity of the algorithms, but also to enable comparisons to other methods. The source code is freely available for download on the web.

LOCALITY OF MHD TURBULENCE IN ISOTHERMAL DISKS

We numerically evolve turbulence driven by the magnetorotational instability (MRI) in a three-dimensional, unstratified shearing box and study its structure using two-point correlation functions. We confirm Fromang and Papaloizous result that shearing box models with zero net magnetic flux are not converged; the dimensionless shear stress α is proportional to the grid scale. We find that the two-point correlation of shows that it is composed of narrow filaments that are swept back by differential rotation into a trailing spiral. The correlation lengths along each of the correlation function principal axes decrease monotonically with the grid scale. For mean azimuthal field models, which we argue are more relevant to astrophysical disks than the zero net field models, we find that: α increases weakly with increasing resolution at fixed box size; α increases slightly as the box size is increased; α increases linearly with net field strength, confirming earlier results; the two-point correlation function of the magnetic field is resolved and converged, and is composed of narrow filaments swept back by the shear; the major axis of the two-point increases slightly as the box size is increased; these results are code independent, based on a comparison of ATHENA and ZEUS runs. The velocity, density, and magnetic fields decorrelate over scales larger than ~H, as do the dynamical terms in the magnetic energy evolution equations. We conclude that MHD turbulence in disks is localized, subject to the limitations imposed by the absence of vertical stratification, the use of an isothermal equation of state, finite box size, finite run time, and finite resolution.

LONG-TERM EVOLUTION OF PROTOSTELLAR AND PROTOPLANETARY DISKS. I. OUTBURSTS

As an initial investigation into the long-term evolution of protostellar disks, we explore the conditions required to explain the large outbursts of disk accretion seen in some young stellar objects. We use one-dimensional time-dependent disk models with a phenomenological treatment of the magnetorotational instability (MRI) and gravitational torques to follow disk evolution over long timescales. Comparison with our previous two-dimensional disk model calculations indicates that the neglect of radial effects and two-dimensional disk structure in the one-dimensional case makes only modest differences in the results; this allows us to use the simpler models to explore parameter space efficiently. We find that the mass infall rates typically estimated for low-mass protostars generally result in AU-scale disk accretion outbursts, as predicted by our previous analysis. We also confirm quasi-steady accretion behavior for high mass infall rates if the values of α-parameter for the MRI are small, while at this high accretion rate convection from the thermal instability may lead to some variations. We further constrain the combinations of the α-parameter and the MRI critical temperature, which can reproduce observed outburst behavior. Our results suggest that dust sublimation may be connected with full activation of the MRI. This is consistent with the idea that small dust captures ions and electrons to suppress the MRI. In a companion paper, we will explore both long-term outburst and disk evolution with this model, allowing for infall from protostellar envelopes with differing angular momenta

WEAK TURBULENCE IN THE HD 163296 PROTOPLANETARY DISK REVEALED BY ALMA CO OBSERVATIONS

Turbulence can transport angular momentum in protoplanetary disks and influence the growth and evolution of planets. With spatially and spectrally resolved molecular emission line measurements provided by (sub)millimeter interferometric observations, it is possible to directly measure non-thermal motions in the disk gas that can be attributed to this turbulence. We report a new constraint on the turbulence in the disk around HD 163296, a nearby young A star, determined from ALMA Science Verification observations of four CO emission lines (the CO(3-2), CO(2-1), 13CO(2-1), and C18O(2-1) transitions). The different optical depths for these lines permit probes of non-thermal line-widths at a range of physical conditions (temperature and density) and depths into the disk interior. We derive stringent limits on the non-thermal motions in the upper layers of the outer disk such that any contribution to the line-widths from turbulence is 30AU) disk than has been previously considered.

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 ? 103. In the intermediate regime (10 ? Am ? 103) 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.

SIMULATIONS OF MAGNETOROTATIONAL TURBULENCE WITH A HIGHER-ORDER GODUNOV SCHEME

We apply a new, second-order Godunov code, Athena, to studies of the magnetorotational instability (MRI) using unstratified shearing box simulations with a uniform net vertical field and a sinusoidally varying zero net vertical field. The Athena results agree well with similar studies that used different numerical algorithms, including the observation that the turbulent energy decreases with increasing resolution in the zero net field model. We conduct analyses to study the flow of energy from differential rotation to turbulent fluctuations to thermalization. A study of the time correlation between the rates of change of different volume-averaged energy components shows that energy injected into turbulent fluctuations dissipates on a timescale of Ω–1, where Ω is the orbital frequency of the local domain. Magnetic dissipation dominates over kinetic dissipation, although not by as great a factor as the ratio of magnetic to kinetic energy. We Fourier-transform the magnetic and kinetic energy evolution equations and, using the assumption that the time-averaged energies are constant, determine the level of numerical dissipation as a function of length scale and resolution. By modeling numerical dissipation as if it were physical in origin, we characterize numerical resistivity and viscosity in terms of effective Reynolds and Prandtl numbers. The resulting effective magnetic Prandtl number is ~2, independent of resolution or initial field geometry. MRI simulations with effective Reynolds and Prandtl numbers determined by numerical dissipation are not equivalent to those where these numbers are set by physical resistivity and viscosity. These results serve, then, as a baseline for future shearing box studies where dissipation is controlled by the inclusion of explicit viscosity and resistivity.

Magnetically driven accretion in protoplanetary discs

We characterize magnetically driven accretion at radii between 1 au and 100 au in protoplanetary discs, using a series of local non-ideal magnetohydrodynamic (MHD) simulations. The simulations assume a Minimum Mass Solar Nebula (MMSN) disc that is threaded by a net vertical magnetic field of specified strength. Confirming previous results, we find that the Hall effect has only a modest impact on accretion at 30 au, and essentially none at 100 au. At 1-10 au the Hall effect introduces a pronounced bi-modality in the accretion process, with vertical magnetic fields aligned to the disc rotation supporting a strong laminar Maxwell stress that is absent if the field is anti-aligned. In the anti-aligned case, we instead find evidence for bursts of turbulent stress at 5-10 au, which we tentatively identify with the non-axisymmetric Hall-shear instability. The presence or absence of these bursts depends upon the details of the adopted chemical model, which suggests that appreciable regions of actual protoplanetary discs might lie close to the borderline between laminar and turbulent behaviour. Given the number of important control parameters that have already been identified in MHD models, quantitative predictions for disc structure in terms of only radius and accretion rate appear to be difficult. Instead, we identify robust qualitative tests of magnetically driven accretion. These include the presence of turbulence in the outer disc, independent of the orientation of the vertical magnetic fields, and a Hall-mediated bi-modality in turbulent properties extending from the region of thermal ionization to 10 au.

Turbulence in global simulations of magnetized thin accretion discs

We use a global magnetohydrodynamic simulation of a geometrically thin accretion disc to investigate the locality and detailed structure of turbulence driven by the magnetorotational instability (MRI). The model disc has an aspect ratio H/R≃ 0.07, and is computed using a higher order Godunov magnetohydrodynamics (MHD) scheme with accurate fluxes. We focus the analysis on late times after the system has lost direct memory of its initial magnetic flux state. The disc enters a saturated turbulent state in which the fastest growing modes of the MRI are well resolved, with a relatively high efficiency of angular momentum transport 〈〈α〉〉≈ 2.5 × 10−2. The accretion stress peaks at the disc mid-plane, above and below which exists a moderately magnetized corona with patches of superthermal field. By analysing the spatial and temporal correlations of the turbulent fields, we find that the spatial structure of the magnetic and kinetic energy is moderately well localized (with correlation lengths along the major axis of 2.5H and 1.5H, respectively), and generally consistent with that expected of homogeneous incompressible turbulence. The density field, conversely, exhibits both a longer correlation length and a long correlation time, results that we ascribe to the importance of spiral density waves within the flow. Consistent with prior results, we show that the mean local stress displays a well-defined correlation with the local vertical flux, and that this relation is apparently causal (in the sense of the flux stimulating the stress) during portions of a global dynamo cycle. We argue that the observed flux–stress relation supports dynamo models in which the structure of coronal magnetic fields plays a central role in determining the dynamics of thin-disc accretion.

The Mass and Size Distribution of Planetesimals Formed by the Streaming Instability. I. The Role of Self-Gravity

We study the formation of planetesimals in protoplanetary disks from the gravitational collapse of solid over-densities generated via the streaming instability. To carry out these studies, we implement and test a particle-mesh self-gravity module for the Athena code that enables the simulation of aerodynamically coupled systems of gas and collisionless self-gravitating solid particles. Upon employment of our algorithm to planetesimal formation simulations, we find that (when a direct comparison is possible) the Athena simulations yield predicted planetesimal properties that agree well with those found in prior work using different numerical techniques. In particular, the gravitational collapse of streaming-initiated clumps leads to an initial planetesimal mass function that is well-represented by a power-law, dN/dM ~ M^{-p}, with p = 1.6 +/- 0.1, which equates to a differential size distribution dN/dR ~ R^{-q}, with q = 2.8 +/- 0.1. We find no significant trends with resolution from a convergence study of up to 512^3 grid zones and Npar = 1.5x10^8 particles. Likewise, the power-law slope appears indifferent to changes in the relative strength of self-gravity and tidal shear, and to the time when (for reasons of numerical economy) self-gravity is turned on, though the strength of these claims is limited by small number statistics. For a typically assumed radial distribution of minimum mass solar nebula solids (assumed here to have dimensionless stopping time \tau = 0.3), our results support the hypothesis that bodies on the scale of large asteroids or Kuiper Belt Objects could have formed as the high-mass tail of a primordial planetesimal population.

Accretion disc dynamo activity in local simulations spanning weak-to-strong net vertical magnetic flux regimes

Strongly magnetized accretion discs around black holes have attractive features that may explain enigmatic aspects of X-ray binary behaviour. The structure and evolution of these discs are governed by a dynamo-like mechanism, which channels part of the accretion power liberated by the magnetorotational instability (MRI) into an ordered toroidal magnetic field. To study dynamo activity, we performed three-dimensional, stratified, isothermal, ideal magnetohydrodynamic shearing box simulations. The strength of the self-sustained toroidal magnetic field depends on the net vertical magnetic flux, which we vary across almost the entire range over which the MRI is linearly unstable. We quantify disc structure and dynamo properties as a function of the initial ratio of mid-plane gas pressure to vertical magnetic field pressure,