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

Philip J. Armitage

Philip J. Armitage appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.HEastro-phastro-ph.GAgr-qcastro-ph.COMachine Learningphysics.flu-dyn

Catalog footprint

What is connected

54works
9topics
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

54 published item(s)

preprint2026arXiv

On the orbital evolution of binaries with polar circumbinary disks

Binaries occur in many astrophysical systems, from young protostellar binaries in star forming regions to supermassive black hole binaries in galaxy centers. In many cases, a circumbinary disk of gas forms around the binary with an orbit that may be misaligned to the binary plane. Misaligned disks around nearly circular binaries evolve into disks that are either aligned or counteraligned with the binary orbit. However, if the binary is sufficiently eccentric, then it can be more likely that the disk ends up in a polar-aligned configuration in which the disk angular momentum vector aligns with the binary eccentricity vector. We use Smoothed Particle Hydrodynamics simulations, evolved to an approximate steady state under mass injection, to determine the orbital evolution of a binary with a polar-aligned disk for a range of binary-disk parameters. We find that, in all of the cases we have simulated, the binary shrinks with time. The decay rate is larger than for binaries surrounded by aligned or retrograde disks with matched disk parameters. The rate of shrinkage is largely unaltered by the size of the sink radii employed for the binary stars, but for small enough sink radii some of the models exhibit long-lived polar circumprimary disks, which are continually fed mass from the circumbinary disk. We discuss our results in the contexts of planet formation in young polar-aligned disks and merging supermassive black holes in galaxy centers.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning

We investigate the utility of deep learning for modeling the clustering of particles that are aerodynamically coupled to turbulent fluids. Using a Lagrangian particle module within the Athena++ hydrodynamics code, we simulate the dynamics of particles in the Epstein drag regime within a periodic domain of isotropic forced hydrodynamic turbulence. This setup is an idealized model relevant to the collisional growth of micron to mm-sized dust particles in early stage planet formation. The simulation data are used to train a U-Net deep learning model to predict gridded three-dimensional representations of the particle density and velocity fields, given as input the corresponding fluid fields. The trained model qualitatively captures the filamentary structure of clustered particles in a highly non-linear regime. We assess model fidelity by calculating metrics of the density field (the radial distribution function) and of the velocity field (the relative velocity and the relative radial velocity between particles). Although trained only on the spatial fields, the model predicts these statistical quantities with errors that are typically <10%. Our results suggest that, given appropriately expanded training data, deep learning could complement direct numerical simulations in predicting particle clustering within turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Instability from high-order resonant chains in wide-separation massive planet systems

Diversity in the properties of exoplanetary systems arises, in part, from dynamical evolution that occurs after planet formation. We use numerical integrations to explore the relative role of secular and resonant dynamics in the long-term evolution of model planetary systems, made up of three equal mass giant planets on initially eccentric orbits. The range of separations studied is dominated by secular processes, but intersects chains of high-order mean-motion resonances. Over time-scales of $10^8$ orbits, the secular evolution of the simulated systems is predominantly regular. High-order resonant chains, however, can be a significant source of angular momentum deficit (AMD), leading to instability. Using a time-series analysis based on a Hilbert transform, we associate instability with broad islands of chaotic evolution. Previous work has suggested that first-order resonances could modify the AMD of nominally secular systems and facilitate secular chaos. We find that higher-order resonances, when present in chains, can have similar impacts.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lecture notes on accretion disk physics

These notes introduce and review some of the physical principles underlying the theory of astrophysical accretion, emphasizing the central roles of angular momentum transport, angular momentum loss, and radiative cooling in determining the structure and evolution of accretion flows. Additional topics covered include the effective viscous theory of thin disks, classical instabilities of disk structure, the evolution of warped or eccentric disks, and the basic properties of waves within disks.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Kozai-Lidov oscillations triggered by a tilt instability of detached circumplanetary discs

Circumplanetary discs can be linearly unstable to the growth of disc tilt in the tidal potential of the star-planet system. We use three-dimensional hydrodynamical simulations to characterize the disc conditions needed for instability, together with its long term evolution. Tilt growth occurs for disc aspect ratios, evaluated near the disc outer edge, of $H/r\gtrsim 0.05$, with a weak dependence on viscosity in the wave-like regime of warp propagation. Lower mass giant planets are more likely to have circumplanetary discs that satisfy the conditions for instability. We show that the tilt instability can excite the inclination to above the threshold where the circumplanetary disc becomes unstable to Kozai--Lidov (KL) oscillations. Dissipation in the Kozai--Lidov unstable regime caps further tilt growth, but the disc experiences large oscillations in both inclination and eccentricity. Planetary accretion occurs in episodic accretion events. We discuss implications of the joint tilt--KL instability for the detectability of circumplanetary discs, for the obliquity evolution of forming giant planets, and for the formation of satellite systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A fast-growing tilt instability of detached circumplanetary disks

Accretion disks in binary systems can exhibit a tilt instability, arising from the interaction between components of the tidal potential and dissipation. Using a linear analysis, we show that the aspect ratios and outer radii of circumplanetary disks provide favorable conditions for tilt growth. We quantify the growth rate of the instability using particle-based ({\sc phantom}) and grid-based ({\sc athena++}) hydrodynamic simulations. For a disk with outer aspect ratio $H/r \simeq 0.1$, initially moderate tilts double on a time scale of about 15-30 binary orbits. Our results imply that detached circumplanetary disks, whose evolution is not entirely controlled by accretion from the circumstellar disk, may commonly be misaligned to the planetary orbital plane. We discuss implications for planetary spin evolution, and possible interactions between the tilt instability and Kozai-Lidov dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bondi on spherically symmetric accretion

Hermann Bondi's 1952 paper "On spherically symmetrical accretion" is recognized as one of the foundations of accretion theory. Although Bondi later remarked that it was "not much more than an examination exercise", his mathematical analysis of spherical accretion on to a point mass has found broad use across fields of astrophysics that were embryonic or non-existent at the time of the paper's publication. In this non-technical review, I describe the motivations for Bondi's work, and briefly discuss some of the applications of Bondi accretion in high energy astrophysics, galaxy formation, and star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Future Simulations of Tidal Disruption Events

Tidal disruption events involve numerous physical processes (fluid dynamics, magnetohydrodynamics, radiation transport, self-gravity, general relativistic dynamics) in highly nonlinear ways, and, because TDEs are transients by definition, frequently in non-equilibrium states. For these reasons, numerical solution of the relevant equations can be an essential tool for studying these events. In this chapter, we present a summary of the key problems of the field for which simulations offer the greatest promise and identify the capabilities required to make progress on them. We then discuss what has been---and what cannot be---done with existing numerical methods. We close with an overview of what methods now under development may do to expand our ability to understand these events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Predicting the long-term stability of compact multiplanet systems

We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $\approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems' angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Simulation of a Compact Object with Outflows Moving Through a Gaseous Background

A compact object moving relative to surrounding gas accretes material and perturbs the density of gas in its vicinity. In the classical picture of Bondi-Hoyle-Lyttleton accretion, the perturbation takes the form of an overdense wake behind the object, which exerts a dynamical friction drag. We use hydrodynamic simulations to investigate how the accretion rate and strength of dynamical friction are modified by the presence of outflow from the compact object. We show that the destruction of the wake by an outflow reduces dynamical friction, and reverses its sign when the outflow is strong enough, in good quantitative agreement with analytic calculations. For a strong isotropic outflow, the outcome on scales that we have simulated is a negative dynamical friction, i.e., net acceleration. For jet-like outflows driven by reprocessed accretion, both the rate of accretion and the magnitude of dynamical friction drop for more powerful jets. The accretion rate is strongly intermittent when the jet points to the same direction as the motion of the compact object. The dynamical effects of outflows may be important for the evolution of compact objects during the common envelope phase of binary systems, and for accreting compact objects and massive stars encountering AGN discs.

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

EXor outbursts from disk amplification of stellar magnetic cycles

EXor outbursts - moderate-amplitude disk accretion events observed in Class I and Class II protostellar sources - have time scales and amplitudes that are consistent with the viscous accumulation and release of gas in the inner disk near the dead zone boundary. We suggest that outbursts are indirectly triggered by stellar dynamo cycles, via poloidal magnetic flux that diffuses radially outward through the disk. Interior to the dead zone the strength of the net field modulates the efficiency of angular momentum transport by the magnetorotational instability. In the dead zone changes in the polarity of the net field may lead to stronger outbursts because of the dominant role of the Hall effect in this region of the disk. At the level of simple estimates we show that changes to kG-strength stellar fields could stimulate disk outbursts on 0.1 AU scales, though this optimistic conclusion depends upon the uncertain efficiency of net flux transport through the inner disk. The model predicts a close association between observational tracers of stellar magnetic activity and EXor events.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the structure of tidally-disrupted stellar debris streams

A tidal disruption event (TDE) -- when a star is destroyed by the immense gravitational field of a supermassive black hole -- transforms a star into a stream of tidally-shredded debris. The properties of this debris ultimately determine the observable signatures of TDEs. Here we derive a simple, self-similar solution for the velocity profile of the debris streams produced from TDEs, and show that this solution agrees extremely well with numerical results. Using this self-similar solution, we calculate an analytic, approximate expression for the radial density profile of the stream. We show that there is a critical adiabatic index that varies as a function of position along the stream above (below) which the stream is unstable (stable) to gravitational fragmentation. We also calculate the impact of heating and cooling on this stability criterion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Planet-disc evolution and the formation of Kozai-Lidov planets

With hydrodynamical simulations we determine the conditions under which an initially coplanar planet-disc system that orbits a member of a misaligned binary star evolves to form a planet that undergoes Kozai-Lidov (KL) oscillations once the disc disperses. These oscillations may explain the large orbital eccentricities, as well as the large misalignments with respect to the spin of the central star, observed for some exoplanets. The planet is assumed to be massive enough to open a gap in the disc. The planet's tilt relative to the binary orbital plane is subject to two types of oscillations. The first type, present at even small inclination angles relative to the binary orbital plane, is due to the interaction of the planet with the disc and binary companion and is amplified by a secular resonance. The second type of oscillation is the KL oscillation that operates on both the planet and disc at larger binary inclination angles. We find that for a sufficiently massive disc, even a relatively low inclination planet-disc system can force a planet to an inclination above the critical KL angle, as a consequence of the first type of tilt oscillation, allowing it to undergo the second type of oscillation. We conclude that the hydrodynamical evolution of a sufficiently massive and inclined disc in a binary system broadens the range of systems that form eccentric and misaligned giant planets to include a wide range of initial misalignment angles (20 to 160 degrees).

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

Accretion disk dynamo as the trigger for X-ray binary state transitions

Magnetohydrodynamic accretion disk simulations suggest that much of the energy liberated by the magnetorotational instability (MRI) can be channeled into large-scale toroidal magnetic fields through dynamo action. Under certain conditions, this field can dominate over gas and radiation pressure in providing vertical support against gravity, even close to the midplane. Using a simple model for the creation of this field, its buoyant rise, and its coupling to the gas, we show how disks could be driven into this magnetically dominated state and deduce the resulting vertical pressure and density profiles. Applying an established criterion for MRI to operate in the presence of a toroidal field, we show that magnetically supported disks can have two distinct MRI-active regions, separated by a "dead zone" where local MRI is suppressed, but where magnetic energy continues to flow upward from the dynamo region below. We suggest that the relative strengths of the MRI zones, and the local poloidal flux, determine the spectral states of X-ray binaries. Specifically, "intermediate" and "hard" accretion states occur when MRI is triggered in the hot, upper zone of the corona, while disks in "soft" states do not develop the upper MRI zone. We discuss the conditions under which various transitions should take place and speculate on the relationship of dynamo activity to the various types of quasi-periodic oscillations that sometimes appear in the hard spectral components. The model also explains why luminous accretion disks in the "soft" state show no signs of the thermal/viscous instability predicted by standard alpha models.

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

Post-periapsis pancakes: sustenance for self-gravity in tidal disruption events

A tidal disruption event, which occurs when a star is destroyed by the gravitational field of a supermassive black hole, produces a stream of debris, the evolution of which ultimately determines the observational properties of the event. Here we show that a post-periapsis caustic -- a location where the locus of gas parcels comprising the stream would collapse into a two-dimensional surface if they evolved solely in the gravitational field of the hole -- occurs when the pericenter distance of the star is on the order of the tidal radius of the hole. It is demonstrated that this "pancake" induces significant density perturbations in the debris stream, and, for stiffer equations of state (adiabatic index $γ\gtrsim 5/3$), these fluctuations are sufficient to gravitationally destabilize the stream, resulting in its fragmentation into bound clumps. The results of our findings are discussed in the context of the observational properties of tidal disruption events.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A mechanism for hysteresis in black hole binary state transitions

We suggest that the hysteretic cycle of black hole state transitions arises from two established properties of accretion disks: the increase in turbulent stress in disks threaded by a net magnetic field and the ability of thick (but not thin) disks to advect such a field radially. During quiescence, magnetic field loops are generated by the magnetorotational instability at the interface between the inner hot flow and outer thin disk. Vertical flux is advected into and accumulates stochastically within the inner flow, where it stimulates the turbulence so that $α\sim 1$. The transition to a geometrically thin inner disk occurs when $L \sim α^2 L_{\rm Edd} \sim L_{\rm Edd}$, and the first "thin" disk to form is itself moderately thick, strongly magnetized, and able to advect field inward. These properties favor episodic jet production. As the accretion rate declines magnetic flux escapes, $α$ decreases to $α\sim 0.01 - 0.1$, and a hot inner flow is not re-established until $L \ll L_{\rm Edd}$. We discuss possible observational consequences of our scenario.

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
preprint2014arXiv

Giant Outbursts in Be/X-ray Binaries

Be/X-ray binary systems exhibit both periodic (Type I) X-ray outbursts and giant (Type II) outbursts, whose origin has remained elusive. We suggest that Type II X-ray outbursts occur when a highly misaligned decretion disk around the Be star becomes eccentric, allowing the compact object companion to capture a large amount of material at periastron. Using 3D smoothed particle hydrodynamics simulations we model the long term evolution of a representative Be/X-ray binary system. We find that periodic (Type I) X-ray outbursts occur when the neutron star is close to periastron for all disk inclinations. Type II outbursts occur for large misalignment angles and are associated with eccentricity growth that occurs on a timescale of about 10 orbital periods. Mass capture from the eccentric decretion disk results in an accretion disk around the neutron star whose estimated viscous time is long enough to explain the extended duration of Type II outbursts. Previous studies suggested that the outbursts are caused by a warped disk but our results suggest that this is not sufficient, the disk must be both highly misaligned and eccentric to initiate a Type II accretion event.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Importance of thermal diffusion in the gravo-magnetic limit cycle

We consider the role of thermal diffusion due to turbulence and radiation on accretion bursts that occur in protoplanetary discs which contain dead zones. Using 1D viscous disc models we show that diffusive radial transport of heat is important during the gravo-magnetic limit cycle, and can strongly modify the duration and frequency of accretion outbursts. When the Prandtl number is large - such that turbulent diffusion of heat is unimportant - radial radiative diffusion reduces the burst duration compared to models with no diffusive transport of heat. When the Prandtl number is small ($\lesssim 25$) we find that turbulent diffusion dominates the radial transport of heat, reducing the burst duration to $\lesssim 10^3$ years as well as increasing the burst frequency. Furthermore, inclusion of radial transport of heat extends the range of infall rates under which the disc undergoes accretion bursts from $10^{-8}$ to $10^{-6}$ M$_\odot$ yr$^{-1}$ with no diffusion, to $10^{-8}$ to $\gtrsim10^{-4}$ M$_\odot$ yr$^{-1}$ with radiative and strong turbulent diffusion. The relative roles of radiative and turbulent thermal diffusion are likely to vary during an accretion burst, but simple estimates suggest that the expected Prandtl numbers are of the order of 10 in protoplanetary discs, and hence that turbulent diffusion is likely to be an important process during accretion outbursts due to the gravo-magnetic limit cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Kozai-Lidov Mechanism in Hydrodynamical Disks

We use three dimensional hydrodynamical simulations to show that a highly misaligned accretion disk around one component of a binary system can exhibit global Kozai-Lidov cycles, where the inclination and eccentricity of the disk are interchanged periodically. This has important implications for accreting systems on all scales, for example, the formation of planets and satellites in circumstellar and circumplanetary disks, outbursts in X-ray binary systems and accretion on to supermassive black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Formation of circumbinary planets in a dead zone

Circumbinary planets have been observed at orbital radii where binary perturbations may have significant effects on the gas disk structure, on planetesimal velocity dispersion, and on the coupling between turbulence and planetesimals. Here, we note that the impact of all of these effects on planet formation is qualitatively altered if the circumbinary disk structure is layered, with a non-turbulent midplane layer (dead zone) and strongly turbulent surface layers. For close binaries, we find that the dead zone typically extends from a radius close to the inner disk edge up to a radius of around 10-20 au from the centre of mass of the binary. The peak in the surface density occurs within the dead zone, far from the inner disk edge, close to the snow line, and may act as a trap for aerodynamically coupled solids. We suggest that circumbinary planet formation may be easier near this preferential location than for disks around single stars. However, dead zones around wide binaries are less likely and hence planet formation may be more difficult there.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The interplay between X-ray photoevaporation and planet formation

We assess the potential of planet formation instigating the early formation of a photoevaporation driven gap, up to radii larger than typical for photoevaporation alone. For our investigation we make use of hydrodynamics models of photoevaporating discs with a giant planet embedded. We find that, by reducing the mass accretion flow onto the star, discs that form giant planets will be dispersed at earlier times than discs without planets by X-ray photoevaporation. By clearing the portion of the disc inner of the planet orbital radius, planet formation induced photoevaporation (PIPE) is able to produce transition disc that for a given mass accretion rate have larger holes when compared to standard X-ray photoevaporation. This constitutes a possible route for the formation of the observed class of accreting transition discs with large holes, which are otherwise difficult to explain by planet formation or photoevaporation alone. Moreover, assuming that a planet is able to filter dust completely, PIPE produces a transition disc with a large hole and may provide a mechanism to quickly shut down accretion. This process appears to be too slow however to explain the observed desert in the population of transition disc with large holes and low mass accretion rates.

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

Debris disks as signposts of terrestrial planet formation. II Dependence of exoplanet architectures on giant planet and disk properties

We present models for the formation of terrestrial planets, and the collisional evolution of debris disks, in planetary systems that contain multiple unstable gas giants. We previously showed that the dynamics of the giant planets introduces a correlation between the presence of terrestrial planets and debris disks. Here we present new simulations that show that this connection is qualitatively robust to changes in: the mass distribution of the giant planets, the width and mass distribution of the outer planetesimal disk, and the presence of gas in the disk. We discuss how variations in these parameters affect the evolution. Systems with equal-mass giant planets undergo the most violent instabilities, and these destroy both terrestrial planets and the outer planetesimal disks that produce debris disks. In contrast, systems with low-mass giant planets efficiently produce both terrestrial planets and debris disks. A large fraction of systems with low-mass outermost giant planets have stable gaps between these planets that are frequently populated by planetesimals. Planetesimal belts between outer giant planets may affect debris disk SEDs. If Earth-mass seeds are present in outer planetesimal disks, the disks radially spread to colder temperatures. We argue that this may explain the very low frequency of > 1 Gyr-old solar-type stars with observed 24 micron excesses. Among the (limited) set of configurations explored, the best candidates for hosting terrestrial planets at ~1 AU are stars older than 0.1-1 Gyr with bright debris disks at 70 micron but with no currently-known giant planets. These systems combine evidence for rocky building blocks, with giant planet properties least likely to undergo destructive dynamical evolution. We predict an anti-correlation between debris disks and eccentric giant planets, and a positive correlation between debris disks and terrestrial planets.

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

Global variation of the dust-to-gas ratio in evolving protoplanetary discs

Recent theories suggest planetesimal formation via streaming and/or gravitational instabilities may be triggered by localized enhancements in the dust-to-gas ratio, and one hypothesis is that sufficient enhancements may be produced in the pile-up of small solid particles inspiralling under aerodynamic drag from the large mass reservoir in the outer disc. Studies of particle pile-up in static gas discs have provided partial support for this hypothesis. Here, we study the radial and temporal evolution of the dust-to-gas ratio in turbulent discs, that evolve under the action of viscosity and photoevaporation. We find that particle pile-ups do not generically occur within evolving discs, particularly if the introduction of large grains is restricted to the inner, dense regions of a disc. Instead, radial drift results in depletion of solids from the outer disc, while the inner disc maintains a dust-to-gas ratio that is within a factor of ~2 of the initial value. We attribute this result to the short time-scales for turbulent diffusion and radial advection (with the mean gas flow) in the inner disc. We show that the qualitative evolution of the dust-to-gas ratio depends only weakly upon the parameters of the disc model (the disc mass, size, viscosity, and value of the Schmidt number), and discuss the implications for planetesimal formation via collective instabilities. Our results suggest that in discs where there is a significant level of midplane turbulence and accretion, planetesimal formation would need to be possible in the absence of large-scale enhancements. Instead, trapping and concentration of particles within local turbulent structures may be required as a first stage of planetesimal formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mini-Oort clouds: Compact isotropic planetesimal clouds from planet-planet scattering

Starting from planetary systems with three giant planets and an outer disk of planetesimals, we use dynamical simulations to show how dynamical instabilities can transform planetesimal disks into 100-1000 AU-scale isotropic clouds. The instabilities involve a phase of planet-planet scattering that concludes with the ejection of one or more planets and the inward-scattering of the surviving gas giant(s) to remove them from direct dynamical contact with the planetesimals. "Mini-Oort clouds" are thus formed from scattered planetesimals whose orbits are frozen by the abrupt disappearance of the perturbing giant planet. Although the planetesimal orbits are virtually isotropic, the surviving giant planets tend to have modest inclinations (typically ~10 degrees) with respect to the initial orbital plane. The collisional lifetimes of mini-Oort clouds are long (10 Myr to >10 Gyr) and there is a window of ~100 Myr or longer during which they produce spherical clouds of potentially observable dust at 70 microns. If the formation channel for hot Jupiters commonly involves planetary close encounters, we predict a correlation between this subset of extrasolar planetary systems and mini-Oort clouds.

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

Debris disks as signposts of terrestrial planet formation

Circumstantial evidence suggests that most known extra-solar planetary systems are survivors of violent dynamical instabilities. Here we explore how giant planet instabilities affect the formation and survival of terrestrial planets. We simulate planetary system evolution around Sun-like stars from initial conditions that comprise: an inner disk of planetesimals and planetary embryos, three giant planets at Jupiter-Saturn distances, and a massive outer planetesimal disk. We then calculate dust production rates and debris disk SEDs assuming that each planetesimal particle represents an ensemble of smaller bodies in collisional equilibrium. We predict a strong correlation between the presence of terrestrial planets and debris disks, mediated by the giant planets. Strong giant planet instabilities destroy all rocky material - including fully-formed terrestrial planets if the instabilities occur late - along with the icy planetesimals. Stable or weakly unstable systems allow terrestrial planets to accrete and significant dust to be produced in their outer regions. Stars older than ~100 Myr with bright cold dust emission (at ~70 microns) signpost the dynamically calm environments conducive to efficient terrestrial accretion. We predict that while the typical eccentricities of terrestrial planets are small, there should exist a novel class of terrestrial planet system whose single planet undergoes large amplitude oscillations in eccentricity and inclination. By scaling to the observed semimajor axis distribution of giant exoplanets, we estimate that terrestrial exoplanets in the same systems should be a few times more abundant at 0.5 AU than giant or terrestrial exoplanets at 1 AU. Finally, we discuss the Solar System, which appears to be unusual in combining a rich terrestrial planet system with a low dust content.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics of Protoplanetary Disks

Protoplanetary disks are quasi-steady structures whose evolution and dispersal determine the environment for planet formation. I review the theory of protoplanetary disk evolution and its connection to observations. Substantial progress has been made in elucidating the physics of potential angular momentum transport processes - including self-gravity, magnetorotational instability, baroclinic instabilities, and magnetic braking - and in developing testable models for disk dispersal via photoevaporation. The relative importance of these processes depends upon the initial mass, size and magnetization of the disk, and subsequently on its opacity, ionization state, and external irradiation. Disk dynamics is therefore coupled to star formation, pre-main-sequence stellar evolution, and dust coagulation during the early stages of planet formation, and may vary dramatically from star to star. The importance of validating theoretical models is emphasized, with the key observations being those that probe disk structure on the scales, between 1 AU and 10 AU, where theory is most uncertain.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrodynamic outcomes of planet scattering in transitional discs

A significant fraction of unstable multiple planet systems likely scatter during the transitional disc phase as gas damping becomes ineffectual. Using an ensemble of FARGO hydrodynamic simulations and MERCURY n-body integrations, we directly follow planet-disc and planet-planet interactions through the clearing phase and on through 50 Myr of dynamical evolution. Disc clearing occurs via X-ray driven photoevaporation. The hydrodynamic evolution of individual scattering systems is complex, and involves phases in which massive planets orbit within eccentric gaps, or accrete directly from the disc without a gap. Comparing the results to a gas-free model, we find that the n-body dynamics and hydrodynamics of scattering into one- and two-planet final states are almost identical. The eccentricity distributions in these channels are almost unaltered by the presence of gas. The hydrodynamic simulations, however, also form low eccentricity three-planet systems in long-term stable configurations, and the admixture of these systems results in modestly lower eccentricities in hydrodynamic as opposed to gas-free simulations. The incidence of these three-planet systems is likely a function of the initial conditions; different planet setups (number or spacing) may change the character of this result. We analyze the properties of surviving multiple planet systems, and show that only a small fraction (a few percent) enter mean-motion resonances after scattering, while a larger fraction form stable resonant chains and avoid scattering entirely. Our results remain consistent with the hypothesis that exoplanet eccentricity results from scattering, though the detailed agreement between observations and gas-free simulation results is likely coincidental. We discuss the prospects for testing scattering models by observing planets or non-axisymmetric gas structure in transitional discs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The butterfly effect in the extreme-mass ratio inspiral problem

Measurements of gravitational waves from the inspiral of a stellar-mass compact object into a massive black hole are unique probes to test General Relativity (GR) and MBH properties, as well as the stellar distribution about these holes in galactic nuclei. Current data analysis techniques can provide us with parameter estimation with very narrow errors. However, an EMRI is not a two-body problem, since other stellar bodies orbiting nearby will influence the capture orbit. Any deviation from the isolated inspiral will induce a small, though observable deviation from the idealised waveform which could be misinterpreted as a failure of GR. Based on conservative analysis of mass segregation in a Milky Way like nucleus, we estimate that the possibility that another star has a semi-major axis comparable to that of the EMRI is non-negligible, although probably very small. This star introduces an observable perturbation in the orbit in the case in which we consider only loss of energy via gravitational radiation. When considering the two first-order non-dissipative post-Newtonian contributions (the periapsis shift of the orbit), the evolution of the orbital elements of the EMRI turns out to be chaotic in nature. The implications of this study are twofold. From the one side, the application to testing GR and measuring MBHs parameters with the detection of EMRIs in galactic nuclei with a millihertz mission will be even more challenging than believed. From the other side, this behaviour could in principle be used as a signature of mass segregation in galactic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The debris disk - terrestrial planet connection

The eccentric orbits of the known extrasolar giant planets provide evidence that most planet-forming environments undergo violent dynamical instabilities. Here, we numerically simulate the impact of giant planet instabilities on planetary systems as a whole. We find that populations of inner rocky and outer icy bodies are both shaped by the giant planet dynamics and are naturally correlated. Strong instabilities -- those with very eccentric surviving giant planets -- completely clear out their inner and outer regions. In contrast, systems with stable or low-mass giant planets form terrestrial planets in their inner regions and outer icy bodies produce dust that is observable as debris disks at mid-infrared wavelengths. Fifteen to twenty percent of old stars are observed to have bright debris disks (at wavelengths of ~70 microns) and we predict that these signpost dynamically calm environments that should contain terrestrial planets.

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

Connections Between Local and Global Turbulence in Accretion Disks

We analyze a suite of global magnetohydrodynamic (MHD) accretion disk simulations in order to determine whether scaling laws for turbulence driven by the magnetorotational instability, discovered via local shearing box studies, are globally robust. The simulations model geometrically-thin disks with zero net magnetic flux and no explicit resistivity or viscosity. We show that the local Maxwell stress is correlated with the self-generated local vertical magnetic field in a manner that is similar to that found in local simulations. Moreover, local patches of vertical field are strong enough to stimulate and control the strength of angular momentum transport across much of the disk. We demonstrate the importance of magnetic linkages (through the low-density corona) between different regions of the disk in determining the local field, and suggest a new convergence requirement for global simulations -- the vertical extent of the corona must be fully captured and resolved. Finally, we examine the temporal convergence of the average stress, and show that an initial long-term secular drift in the local flux-stress relation dies away on a time scale that is consistent with turbulent mixing of the initial magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Planet-planet scattering in planetesimal disks II: Predictions for outer extrasolar planetary systems

We develop an idealized dynamical model to predict the typical properties of outer extrasolar planetary systems, at radii beyond 5 AU. Our hypothesis is that dynamical evolution in outer planetary systems is controlled by a combination of planet-planet scattering and planetary interactions with an exterior disk of small bodies ("planetesimals"). Using 5,000 long duration N-body simulations, we follow the evolution of three planets surrounded by a 50 Earth mass primordial planetesimal disk. For large planet masses (above that of Saturn) the influence of the disk is modest, and we recover the observed eccentricity distribution of extrasolar planets (observed primarily at smaller radii). We explain the observed mass dependence of the eccentricity by invoking strong correlations between planet masses in the same system. For lower mass planets we observe diverse dynamical behavior: strong scattering events, sudden jumps in eccentricity due to resonance crossings, and re-circularization of scattered low-mass planets in the disk. We present distributions of the final eccentricity and inclination, and discuss how they vary with planet mass and initial system architecture. We predict a transition to lower eccentricities for low mass planets at radii where disks influence the dynamics. Radial velocity measurements capable of detecting planets with K~5 m/s and periods in excess of 10 years will constrain this regime. We also study the population of resonant and non-resonant multiple planet systems. We show that, among systems with Jupiter-mass planets that avoid close encounters, the planetesimal disk acts as a damping mechanism that frequently populates mean motion resonances. Resonant chains ought to be common among massive planets in outer planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Planet-planet scattering leads to tightly packed planetary systems

The known extrasolar multiple-planet systems share a surprising dynamical attribute: they cluster just beyond the Hill stability boundary. Here we show that the planet-planet scattering model, which naturally explains the observed exoplanet eccentricity distribution, can reproduce the observed distribution of dynamical configurations. We calculated how each of our scattered systems would appear over an appropriate range of viewing geometries; as Hill stability is weakly dependent on the masses, the mass-inclination degeneracy does not significantly affect our results. We consider a wide range of initial planetary mass distributions and find that some are poor fits to the observed systems. In fact, many of our scattering experiments overproduce systems very close to the stability boundary. The distribution of dynamical configurations of two-planet systems actually may provide better discrimination between scattering models than the distribution of eccentricity. Our results imply that, at least in their inner regions which are weakly affected by gas or planetesimal disks, planetary systems should be "packed", with no large gaps between planets.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Extrasolar planetary dynamics with a generalized planar Laplace-Lagrange secular theory

The dynamical evolution of nearly half of the known extrasolar planets in multiple-planet systems may be dominated by secular perturbations. The commonly high eccentricities of the planetary orbits calls into question the utility of the traditional Laplace-Lagrange (LL) secular theory in analyses of the motion. We analytically generalize this theory to fourth-order in the eccentricities, compare the result with the second-order theory and octupole-level theory, and apply these theories to the likely secularly-dominated HD 12661, HD 168443, HD 38529 and Ups And multi-planet systems. The fourth-order scheme yields a multiply-branched criterion for maintaining apsidal libration, and implies that the apsidal rate of a small body is a function of its initial eccentricity, dependencies which are absent from the traditional theory. Numerical results indicate that the primary difference the second and fourth-order theories reveal is an alteration in secular periodicities, and to a smaller extent amplitudes of the planetary eccentricity variation. Comparison with numerical integrations indicates that the improvement afforded by the fourth-order theory over the second-order theory sometimes dwarfs the improvement needed to reproduce the actual dynamical evolution. We conclude that LL secular theory, to any order, generally represents a poor barometer for predicting secular dynamics in extrasolar planetary systems, but does embody a useful tool for extracting an accurate long-term dynamical description of systems with small bodies and/or near-circular orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Theory of disk accretion onto supermassive black holes

Accretion onto supermassive black holes produces both the dramatic phenomena associated with active galactic nuclei and the underwhelming displays seen in the Galactic Center and most other nearby galaxies. I review selected aspects of the current theoretical understanding of black hole accretion, emphasizing the role of magnetohydrodynamic turbulence and gravitational instabilities in driving the actual accretion and the importance of the efficacy of cooling in determining the structure and observational appearance of the accretion flow. Ongoing investigations into the dynamics of the plunging region, the origin of variability in the accretion process, and the evolution of warped, twisted, or eccentric disks are summarized.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Accretion during the merger of supermassive black holes

We study the evolution of disk accretion during the merger of supermassive black hole binaries in galactic nuclei. In hierarchical galaxy formation models, the most common binaries are likely to arise from minor galactic mergers, and have unequal mass black holes. Once such a binary becomes embedded in an accretion disk at a separation of the order of 0.1 pc, the merger proceeds in two distinct phases. During the first phase, the loss of orbital angular momentum to the gaseous disk shrinks the binary on a timescale of ~10 Myr. The accretion rate onto the primary black hole is not increased, and can be substantially reduced, during this disk-driven migration. At smaller separations, gravitational radiation becomes the dominant angular momentum loss process, and any gas trapped inside the orbit of the secondary is driven inwards by the inspiralling black hole. The implied accretion rate just prior to coalescence exceeds the Eddington limit, so the final merger is likely to occur within a common envelope formed from the disrupted inner disk, and be accompanied by high velocity (~10,000 km/s) outflows.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Simulations of accretion flows crossing the last stable orbit

We use three dimensional magnetohydrodynamic simulations, in a pseudo-Newtonian potential, to study geometrically thin accretion disc flows crossing the marginally stable circular orbit around black holes. We concentrate on vertically unstratified and isothermal disk models, but also consider a model that includes stratification. In all cases, we find that the sonic point lies just inside the last stable orbit, with a modest increase in the importance of magnetic field energy, relative to the thermal energy, observed interior to this radius. The time-averaged gradient of the specific angular momentum of the flow, (dl/dr), is close to zero within the marginally stable orbit, despite the presence of large fluctuations and continuing magnetic stress in the plunging region. The result that the specific angular momentum is constant within the last stable orbit is in general agreement with traditional disk models computed using a zero-torque boundary condition at the last stable orbit.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Hydrodynamics of the stream-disk impact in interacting binaries

We use hydrodynamic simulations to provide quantitative estimates of the effects of the impact of the accretion stream on disks in interacting binaries. For low accretion rates, efficient radiative cooling of the hotspot region can occur, and the primary consequence of the stream impact is stream overflow toward smaller disk radii. The stream is well described by a ballistic trajectory, but larger masses of gas are swept up and overflow at smaller, but still highly supersonic, velocities. If cooling is inefficient, overflow still occurs, but there is no coherent stream inward of the disk rim. Qualitatively, the resulting structure appears as a bulge extending downstream along the disk rim. We calculate the mass fraction and velocity of the overflowing component as a function of the important system parameters, and discuss the implications of the results for X-ray observations and doppler tomography of cataclysmic variables, low-mass X-ray binaries and supersoft X-ray sources.

0 citations0 reviewsNo code linkedOpen access