BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jacob B. Simon

Jacob B. Simon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPastro-ph.HEastro-ph.GAgr-qcMachine Learning

Catalog footprint

What is connected

29works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

29 published item(s)

preprint2022arXiv

Resilience of Planetesimal Formation in Weakly-Reinforced Pressure Bumps

The discovery that axisymmetric dust rings are ubiquitous in protoplanetary disks has provoked a flurry of research on the role of pressure bumps in planet formation. High-resolution simulations by our group have shown that even a modest bump can collect enough dust to trigger planetesimal formation by the streaming instability. In this work, we probe the limits of planetesimal formation when the external source of pressure bump reinforcement is extremely weak. We conduct simulations of radially elongated shearing boxes to capture the entire bump, which is generated and maintained over some timescale $t_{\rm reinf}$ by a Newtonian relaxation scheme. We find that planetesimal formation is extremely resilient for cm-sized grains. We reduced the strength of reinforcement by up to a factor of 100 and the location and initial masses of planetesimals were essentially unaffected. However, we do find that strong reinforcement causes much faster pebble drift compared to the the standard pebble drift rates. The resulting larger pebble flux enhances the planetesimal growth rate by pebble accretion. We hypothesize that to sustain the bump, our code has to extract angular momentum (the strength of this negative torque depends on $t_{\rm reinf}$), and some of this torque is transferred to the particles, causing them to drift faster for a stronger torque (i.e., smaller $t_{\rm reinf}$). Since any physical process that sustains a pressure bump must do so by torquing the gas, we conjecture that the effect on pebble drift is a real phenomenon, motivating further work with physically realistic sources to generate the bump.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Streaming Instability Cannot Form Planetesimals from mm-size Grains in Pressure Bumps

We present evidence that it is unlikely that the streaming instability (SI) can form planetesimals from mm grains inside axisymmetric pressure bumps. We conducted the largest simulation of the SI so far (7 million CPU hours), consisting of a large slice of the disk with mm grains, a solar-like dust-to-gas ratio ($Z = 0.01$), and the largest pressure bump that does not cause gravitational instability (GI) in the particle layer. We used a high resolution of $1000/H$ to resolve as many SI unstable modes as possible. The simulation produced a long-lived particle over-density far exceeding the SI criteria (i.e., a critical solid abundance to headwind parameter ratio $Z/Π$) where strong clumping would occur if these conditions were present over an extended region of the disk; yet we observed none. The likely reason is that the time it takes particles to cross the high-$Z/Π$ region ($t_{\rm cross}$) is shorter than the growth timescale of the SI ($t_{\rm grow}$). We propose an added criterion for planetesimal formation by the SI -- that $t_{\rm cross} > t_{\rm grow}$. We show that any bump larger than the one in this run would form planetesimals by the GI instead of the SI. Our results significantly restrict the pathways to planet formation: Either protoplanetary disks regularly form grains larger than 1~mm, or planetesimals do not form by the SI in axisymmetric pressure bumps. Since bumps large enough to induce the GI are likely Rossby-wave unstable, we propose that mm grains may only form planetesimals in vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Using Bayesian Deep Learning to infer Planet Mass from Gaps in Protoplanetary Disks

Planet induced sub-structures, like annular gaps, observed in dust emission from protoplanetary disks provide a unique probe to characterize unseen young planets. While deep learning based model has an edge in characterizing the planet's properties over traditional methods, like customized simulations and empirical relations, it lacks in its ability to quantify the uncertainty associated with its predictions. In this paper, we introduce a Bayesian deep learning network "DPNNet-Bayesian" that can predict planet mass from disk gaps and provides uncertainties associated with the prediction. A unique feature of our approach is that it can distinguish between the uncertainty associated with the deep learning architecture and uncertainty inherent in the input data due to measurement noise. The model is trained on a data set generated from disk-planet simulations using the \textsc{fargo3d} hydrodynamics code with a newly implemented fixed grain size module and improved initial conditions. The Bayesian framework enables estimating a gauge/confidence interval over the validity of the prediction when applied to unknown observations. As a proof-of-concept, we apply DPNNet-Bayesian to dust gaps observed in HL Tau. The network predicts masses of $ 86.0 \pm 5.5 M_{\Earth} $, $ 43.8 \pm 3.3 M_{\Earth} $, and $ 92.2 \pm 5.1 M_{\Earth} $ respectively, which are comparable to other studies based on specialized simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Protoplanetary Disk Rings as Sites for Planetesimal Formation

Axisymmetric dust rings are a ubiquitous feature of young protoplanetary disks. These rings are likely caused by pressure bumps in the gas profile; a small bump can induce a traffic jam-like pattern in the dust density, while a large bump may halt radial dust drift entirely. The resulting increase in dust concentration may trigger planetesimal formation by the streaming instability (SI), as the SI itself requires some initial concentration. Here we present the first 3D simulations of planetesimal formation in the presence of a pressure bump modeled specifically after those observed by ALMA. In particular, we place a pressure bump at the center of a large 3D shearing box, along with an initial solid-to-gas ratio of $Z = 0.01$, and we include both particle back-reaction and particle self-gravity. We consider both mm-sized and cm-sized particles separately. For simulations with cm-sized particles, we find that even a small pressure bump leads to the formation of planetesimals via the streaming instability; a pressure bump does {\it not} need to fully halt radial particle drift for the SI to become efficient. Furthermore, pure gravitational collapse via concentration in pressure bumps (such as would occur at sufficiently high concentrations and without the streaming instability) is not responsible for planetesimal formation. For mm-sized particles, we find tentative evidence that planetesimal formation does not occur. This result, if it holds up at higher resolution and for a broader range of parameters, could put strong constraints on where in protoplanetary disks planetesimals can form. Ultimately, however, our results suggest that for cm-sized particles, planetesimal formation in pressure bumps is an extremely robust process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Strongly magnetized accretion disks: structure and accretion from global magnetohydrodynamic simulations

We use global magnetohydrodynamic simulations to study the influence of net vertical magnetic fields on the structure of geometrically thin ($H/r \approx 0.05$) accretion disks in the Newtonian limit. We consider initial mid-plane gas to magnetic pressure ratios $β_0 = 1000,\, 300$ and $100$, spanning the transition between weakly and strongly magnetized accretion regimes. We find that magnetic pressure is important for the disks' vertical structure in all three cases, with accretion occurring at $z/R\approx 0.2$ in the two most strongly magnetized models. The disk midplane shows outflow rather than accretion. Accretion through the surface layers is driven mainly by stress due to coherent large scale magnetic field rather than by turbulent stress. Equivalent viscosity parameters measured from our simulations show similar dependencies on initial $β_0$ to those seen in shearing box simulations, though the disk midplane is not magnetic pressure dominated even for the strongest magnetic field case. Winds are present but are not the dominant driver of disk evolution. Over the (limited) duration of our simulations, we find evidence that the net flux attains a quasi-steady state at levels that can stably maintain a strongly magnetized disk. We suggest that geometrically thin accretion disks in observed systems may commonly exist in a magnetically "elevated" state, characterized by non-zero but modest vertical magnetic fluxes, with potentially important implications for disk phenomenology in X-ray binaries (XRBs) and active galactic nuclei (AGN).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Visualizing the Kinematics of Planet Formation

A stunning range of substructures in the dust of protoplanetary disks is routinely observed across a range of wavelengths. These gaps, rings and spirals are highly indicative of a population of unseen planets, hinting at the possibility of current observational facilities being able to capture planet-formation in action. Over the last decade, our understanding of the influence of a young planet on the dynamical structure of its parental disk has progressed significantly, revealing a host of potentially observable features which would betray the presence of a deeply embedded planet. In concert, recent observations have shown that subtle perturbations in the kinematic structure of protoplanetary disks are found in multiple sources, potentially the characteristic disturbances associated with embedded planets. In this work, we review the theoretical background of planet-disk interactions, focusing on the kinematical features, and the current methodologies used to observe these interactions in spatially and spectrally resolved observations. We discuss the potential pit falls of such kinematical detections of planets, providing best-practices for imaging and analysing interferometric data, along with a set of criteria to use as a benchmark for any claimed detection of embedded planets. We finish with a discussion on the current state of simulations in regard to planet-disk interactions, highlighting areas of particular interest and future directions which will provide the most significant impact in our search for embedded planets. This work is the culmination of the 'Visualizing the Kinematics of Planet Formation' workshop, held in October 2019 at the Center for Computational Astrophysics at the Flatiron Institute in New York City.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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, $β_0^{\rm mid} = p_{\rm gas} / p_B$. For $10^5 \geq β_0^{\rm mid} \geq 10$ the effective $α$-viscosity parameter scales as a power-law. Dynamo activity persists up to and including $β_0^{\rm mid} = 10^2$, at which point the entire vertical column of the disc is magnetic pressure-dominated. Still stronger fields result in a highly inhomogeneous disc structure, with large density fluctuations. We show that the turbulent steady state $β^{\rm mid}$ in our simulations is well-matched by the analytic model of Begelman et al. (2015) describing the creation and buoyant escape of toroidal field, while the vertical structure of the disc can be broadly reproduced using this model. Finally, we discuss the implications of our results for observed properties of X-ray binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Prompt planetesimal formation beyond the snow line

We develop a simple model to predict the radial distribution of planetesimal formation. The model is based on the observed growth of dust to mm-sized particles, which drift radially, pile-up, and form planetesimals where the stopping time and dust-to-gas ratio intersect the allowed region for streaming instability-induced gravitational collapse. Using an approximate analytic treatment, we first show that drifting particles define a track in metallicity--stopping time space whose only substantial dependence is on the disk's angular momentum transport efficiency. Prompt planetesimal formation is feasible for high particle accretion rates (relative to the gas, $\dot{M}_p / \dot{M} > 3 \times 10^{-2}$ for $α= 10^{-2}$), that could only be sustained for a limited period of time. If it is possible, it would lead to the deposition of a broad and massive belt of planetesimals with a sharp outer edge. Including turbulent diffusion and vapor condensation processes numerically, we find that a modest enhancement of solids near the snow line occurs for cm-sized particles, but that this is largely immaterial for planetesimal formation. We note that radial drift couples planetesimal formation across radii in the disk, and suggest that considerations of planetesimal formation favor a model in which the initial deposition of material for giant planet cores occurs well beyond the snow line.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Strongly magnetized accretion discs require poloidal flux

Motivated by indirect observational evidence for strongly magnetized accretion discs around black holes, and the novel theoretical properties of such solutions, we investigate how a strong magnetization state can develop and persist. To this end, we perform local simulations of accretion discs with an initially purely toroidal magnetic field of equipartition strength. We demonstrate that discs with zero net vertical magnetic flux and realistic boundary conditions cannot sustain a strong toroidal field. However, a magnetic pressure-dominated disc can form from an initial configuration with a sufficient amount of net vertical flux and realistic boundary conditions. Our results suggest that poloidal flux is a necessary prerequisite for the sustainability of strongly magnetized accretion discs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Influence of Magnetic Field Geometry on the Formation of Close-in Exoplanets

Approximately half of Sun-like stars harbor exoplanets packed within a radius of ~0.3 AU, but the formation of these planets and why they form in only half of known systems are still not well understood. We employ a one-dimensional steady state model to gain physical insight into the origin of these close-in exoplanets. We use Shakura & Sunyaev alpha values extracted from recent numerical simulations of protoplanetary disk accretion processes in which the magnitude of alpha, and thus the steady-state gas surface density, depends on the orientation of large scale magnetic fields with respect to the disk's rotation axis. Solving for the metallicity as a function of radius, we find that for fields anti-aligned with the rotation axis, the inner regions of our model disk often falls within a region of parameter space not suitable for planetesimal formation, whereas in the aligned case, the inner disk regions are likely to produce planetesimals through some combination of streaming instability and gravitational collapse, though the degree to which this is true depends on the assumed parameters of our model. More robustly, the aligned field case always produces higher concentrations of solids at small radii compared to the anti-aligned case. In the in situ formation model, this bimodal distribution of solid enhancement leads directly to the observed dichotomy in exoplanet orbital distances.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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 τ= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Turbulence, Transport and Waves in Ohmic Dead Zones

We use local numerical simulations to study a vertically stratified accretion disk with a resistive mid-plane that damps magnetohydrodynamic (MHD) turbulence. This is an idealized model for the dead zones that may be present at some radii in protoplanetary and dwarf novae disks. We vary the relative thickness of the dead and active zones to quantify how forced fluid motions in the dead zone change. We find that the residual Reynolds stress near the mid-plane decreases with increasing dead zone thickness, becoming negligible in cases where the active to dead mass ratio is less than a few percent. This implies that purely Ohmic dead zones would be vulnerable to episodic accretion outbursts via the mechanism of Martin & Lubow (2011). We show that even thick dead zones support a large amount of kinetic energy, but this energy is largely in fluid motions that are inefficient at angular momentum transport. Confirming results from Oishi & Mac Low (2009), the perturbed velocity field in the dead zone is dominated by an oscillatory, vertically extended circulation pattern with a low frequency compared to the orbital frequency. This disturbance has the properties predicted for the lowest order r mode in a hydrodynamic disk. We suggest that in a global disk similar excitations would lead to propagating waves, whose properties would vary with the thickness of the dead zone and the nature of the perturbations (isothermal or adiabatic). Flows with similar amplitudes would buckle settled particle layers and could reduce the efficiency of pebble accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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 <3% of the local sound speed. These limits are approximately an order of magnitude lower than theoretical predictions for full-blown MHD turbulence driven by the magneto-rotational instability, potentially suggesting that this mechanism is less efficient in the outer (R>30AU) disk than has been previously considered.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficiency of Particle Trapping in the Outer Regions of Protoplanetary Disks

We investigate the strength of axisymmetric local pressure maxima (zonal flows) in the outer regions of protoplanetary disks, where ambipolar diffusion reduces turbulent stresses driven by the magnetorotational instability. Using local numerical simulations we show that in the absence of net vertical magnetic fields, the strength of turbulence in the ambipolar dominated region of the disk is low and any zonal flows that are present are weak. For net fields strong enough to yield observed protostellar accretion rates, however, zonal flows with a density amplitude of 10-20% are formed. These strengths are comparable to those seen in simulations of ideal MHD disk turbulence. We investigate whether these zonal flows are able to reverse the inward radial drift of solids, leading to prolonged and enhanced concentration as a prelude to planetesimal formation. For commonly assumed mean surface density profiles (surface density proportional to radius to the -1/2 power or steeper) we find that the predicted perturbations to the background disk profile do not correspond to local pressure maxima. This is a consequence of radial width of the simulated zonal flows, which is larger than was assumed in prior analytic models of particle trapping. These larger scale flows would only trap particles for higher amplitude fluctuations than observed. We conclude that zonal flows are likely to be present in the outer regions of protoplanetary disks and are potentially large enough to be observable, but are unlikely to lead to strong particle trapping.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantifying energetics and dissipation in magnetohydrodynamic turbulence

We perform a suite of two- and three-dimensional magnetohydrodynamic (MHD) simulations with the Athena code of the non-driven Kelvin-Helmholtz instability in the subsonic, weak magnetic field limit. Focusing the analysis on the non-linear turbulent regime, we quantify energy transfer on a scale-by-scale basis and identify the physical mechanisms responsible for energy exchange by developing the diagnostic known as spectral energy transfer function analysis. At late times when the fluid is in a state of MHD turbulence, magnetic tension mediates the dominant mode of energy injection into the magnetic reservoir, whereby turbulent fluid motions twist and stretch the magnetic field lines. This generated magnetic energy turbulently cascades to smaller scales, while being exchanged backwards and forwards with the kinetic energy reservoir, until finally being dissipated. Incorporating explicit dissipation pushes the dissipation scale to larger scales than if the dissipation were entirely numerical. For scales larger than the dissipation scale, we show that the physics of energy transfer in decaying MHD turbulence is robust to numerical effects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Two timescale dispersal of magnetized protoplanetary disks

Protoplanetary disks are likely to be threaded by a weak net flux of vertical magnetic field that is a remnant of the much larger fluxes present in molecular cloud cores. If this flux is approximately conserved its dynamical importance will increase as mass is accreted, initially by stimulating magnetorotational disk turbulence and subsequently by enabling wind angular momentum loss. We use fits to numerical simulations of ambipolar dominated disk turbulence to construct simplified one dimensional evolution models for weakly magnetized protoplanetary disks. We show that the late onset of significant angular momentum loss in a wind can give rise to "two timescale" disk evolution in which a long phase of viscous evolution precedes rapid dispersal as the wind becomes dominant. The wide dispersion in disk lifetimes could therefore be due to varying initial levels of net flux. Magnetohydrodynamic (MHD) wind triggered dispersal differs from photoevaporative dispersal in predicting mass loss from small (less that 1 AU) scales, where thermal winds are suppressed. Our specific models are based on a limited set of simulations that remain uncertain, but qualitatively similar evolution appears likely if mass is lost from disks more quickly than flux, and if MHD winds become important as the plasma beta decreases.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Emergent Mesoscale Phenomena in Magnetized Accretion Disc Turbulence

We study how the structure and variability of magnetohydrodynamic (MHD) turbulence in accretion discs converge with domain size. Our results are based on a series of vertically stratified local simulations, computed using the Athena code, that have fixed spatial resolution, but varying radial and azimuthal extent (from ΔR = 0.5H to 16H, where H is the vertical scale height). We show that elementary local diagnostics of the turbulence, including the Shakura-Sunyaev α parameter, the ratio of Maxwell stress to magnetic energy, and the ratio of magnetic to fluid stresses, converge to within the precision of our measurements for spatial domains of radial size Lx \geq 2H. We obtain α = 0.02-0.03, consistent with recent results. Very small domains (Lx = 0.5H) return anomalous results, independent of spatial resolution. The convergence with domain size is only valid for a limited set of diagnostics: larger spatial domains admit the emergence of dynamically important mesoscale structures. In our largest simulations, the Maxwell stress shows a significant large scale non-local component, while the density develops long-lived axisymmetric perturbations (zonal flows) at the 20% level. Most strikingly, the variability of the disc in fixed-sized patches decreases strongly as the simulation volume increases. We find generally good agreement between our largest local simulations and global simulations with comparable spatial resolution. There is no direct evidence that the presence of curvature terms or radial gradients in global calculations materially affect the turbulence, except to perhaps introduce an outer radial scale for mesoscale structures. The demonstrated importance of mean magnetic fields, seen in both large local and global simulations implies that the growth and saturation of these fields is likely of critical importance for the evolution of accretion discs. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Turbulent Linewidths as a Diagnostic of Self-Gravity in Protostellar Discs

We use smoothed particle hydrodynamics simulations of massive protostellar discs to investigate the predicted broadening of molecular lines from discs in which self-gravity is the dominant source of angular momentum transport. The simulations include radiative transfer, and span a range of disc-to-star mass ratios between 0.25 and 1.5. Subtracting off the mean azimuthal flow velocity, we compute the distribution of the in-plane and perpendicular peculiar velocity due to large scale structure and turbulence induced by self-gravity. For the lower mass discs, we show that the characteristic peculiar velocities scale with the square root of the effective turbulent viscosity parameter, as expected from local turbulent-disc theory. The derived velocities are anisotropic, with substantially larger in-plane than perpendicular values. As the disc mass is increased, the validity of the locally determined turbulence approximation breaks down, and this is accompanied by anomalously large in-plane broadening. There is also a high variance due to the importance of low-m spiral modes. For low-mass discs, the magnitude of in-plane broadening is, to leading order, equal to the predictions from local disc theory and cannot constrain the source of turbulence. However, combining our results with prior evaluations of turbulent broadening expected in discs where the magnetorotational instability (MRI) is active, we argue that self-gravity may be distinguishable from the MRI in these systems if it is possible to measure the anisotropy of the peculiar velocity field with disc inclination. Furthermore, for large mass discs, the dominant contribution of large-scale modes is a distinguishing characteristic of self-gravitating turbulence versus MRI driven turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resistivity-driven State Changes in Vertically Stratified Accretion Disks

We investigate the effect of shear viscosity and Ohmic resistivity on the magnetorotational instability (MRI) in vertically stratified accretion disks through a series of local simulations with the Athena code. First, we use a series of unstratified simulations to calibrate physical dissipation as a function of resolution and background field strength; the effect of the magnetic Prandtl number, Pm = viscosity/resistivity, on the turbulence is captured by ~32 grid zones per disk scale height, H. In agreement with previous results, our stratified disk calculations are characterized by a subthermal, predominately toroidal magnetic field that produces MRI-driven turbulence for |z| < 2 H. Above |z| = 2 H, magnetic pressure dominates and the field is buoyantly unstable. Large scale radial and toroidal fields are also generated near the mid-plane and subsequently rise through the disk. The polarity of this mean field switches on a roughly 10 orbit period in a process that is well-modeled by an alpha-omega dynamo. Turbulent stress increases with Pm but with a shallower dependence compared to unstratified simulations. For sufficiently large resistivity, on the order of cs H/1000, where cs is the sound speed, MRI turbulence within 2 H of the mid-plane undergoes periods of resistive decay followed by regrowth. This regrowth is caused by amplification of toroidal field via the dynamo. This process results in large amplitude variability in the stress on 10 to 100 orbital timescales, which may have relevance for partially ionized disks that are observed to have high and low accretion states.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Turbulence in Global Simulations of Magnetized Thin Accretion Disks

We use a global magnetohydrodynamic simulation of a geometrically thin accretion disk to investigate the locality and detailed structure of turbulence driven by the magnetorotational instability (MRI). The model disk has an aspect ratio $H / R \simeq 0.07$, and is computed using a higher-order Godunov 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 disk 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 $< < α> > \approx 2.5 \times 10^{-2}$. The accretion stress peaks at the disk midplane, above and below which exists a moderately magnetized corona with patches of superthermal field. By analyzing 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 from homogenous incompressible turbulence. The density field, conversely, exhibits both a longer correlation length and a long correlation time, results which 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-disk accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Turbulent Linewidths in Protoplanetary Disks: Predictions from Numerical Simulations

Sub-mm observations of protoplanetary disks now approach the acuity needed to measure the turbulent broadening of molecular lines. These measurements constrain disk angular momentum transport, and furnish evidence of the turbulent environment within which planetesimal formation takes place. We use local magnetohydrodynamic (MHD) simulations of the magnetorotational instability (MRI) to predict the distribution of turbulent velocities in low mass protoplanetary disks, as a function of radius and height above the mid-plane. We model both ideal MHD disks, and disks in which Ohmic dissipation results in a dead zone of suppressed turbulence near the mid-plane. Under ideal conditions, the disk mid-plane is characterized by a velocity distribution that peaks near v \simeq 0.1cs (where cs is the local sound speed), while supersonic velocities are reached at z > 3H (where H is the pressure scale height). Residual velocities of v \approx 0.01cs persist near the mid-plane in dead zones, while the surface layers remain active. Anisotropic variation of the linewidth with disk inclination is modest. We compare our MHD results to hydrodynamic simulations in which large-scale forcing is used to initiate similar turbulent velocities. We show that the qualitative trend of increasing v with height, seen in the MHD case, persists for forced turbulence and is likely a generic property of disk turbulence. Percent level determinations of v at different heights within the disk, or spatially resolved observations that probe the inner disk containing the dead zone region, are therefore needed to test whether the MRI is responsible for protoplanetary disk turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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 (Zhu et al. 2009b, Z2009b) 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 (Zhu et al. 2009a,Z2009a). We also confirm quasi-steady accretion behavior for high mass infall rates if the values of $α$-parameter for the magnetorotational instability is 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 later paper we will explore both long-term outburst and disk evolution with this model, allowing for infall from protostellar envelopes with differing angular momenta.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Locality of MHD Turbulence in Isothermal Disks

We numerically evolve turbulence driven by the magnetorotational instability (MRI) in a 3D, unstratified shearing box and study its structure using two-point correlation functions. We confirm Fromang and Papaloizou's 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 the magnetic field 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 $\sim 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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Viscous and Resistive Effects on the MRI with a Net Toroidal Field

Resistivity and viscosity have a significant role in establishing the energy levels in turbulence driven by the magnetorotational instability (MRI) in local astrophysical disk models. This study uses the Athena code to characterize the effects of a constant shear viscosity νand Ohmic resistivity ηin unstratified shearing box simulations with a net toroidal magnetic flux. A previous study of shearing boxes with zero net magnetic field performed with the ZEUS code found that turbulence dies out for values of the magnetic Prandtl number, P_m = ν/η, below P_m \sim 1; for P_m \gtrsim 1, time- and volume-averaged stress levels increase with P_m. We repeat these experiments with Athena and obtain consistent results. Next, the influence of viscosity and resistivity on the toroidal field MRI is investigated both for linear growth and for fully-developed turbulence. In the linear regime, a sufficiently large νor ηcan prevent MRI growth; P_m itself has little direct influence on growth from linear perturbations. By applying a range of values for νand ηto an initial state consisting of fully developed turbulence in the presence of a background toroidal field, we investigate their effects in the fully nonlinear system. Here, increased viscosity enhances the turbulence, and the turbulence decays only if the resistivity is above a critical value; turbulence can be sustained even when P_m < 1, in contrast to the zero net field model. While we find preliminary evidence that the stress converges to a small range of values when νand ηbecome small enough, the influence of dissipation terms on MRI-driven turbulence for relatively large ηand νis significant, independent of field geometry.

0 citations0 reviewsNo code linkedOpen access