sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Armitage, Phil

Normalized to: Armitage, P.

145 article(s) in total. 267 co-authors, from 1 to 20 common article(s). Median position in authors list is 2,0.

[1]  [pdf] - 2038484
Turbulence Regulates the Rate of Planetesimal Formation via Gravitational Collapse
Comments: Submitted to ApJ, Comments welcome
Submitted: 2020-01-27
We study how the interaction between the streaming instability and intrinsic gas-phase turbulence affects planetesimal formation via gravitational collapse. Turbulence impedes the formation of dense particle clumps by acting as an effective turbulent diffusivity, but it can also promote planetesimal formation by concentrating solids, for example in zonal flows. We quantify the effect of turbulent diffusivity using numerical simulations of the streaming instability in small local domains, forced with velocity perturbations that establish approximately Kolmogorov-like fluid turbulence. We find that planetesimal formation is strongly suppressed by turbulence once velocity fluctuations exceed a threshold value of $\delta v^2 \simeq 10^{-3.5} - 10^{-3} c_s^2$. Turbulence whose strength is just below the threshold reduces the rate of solid material being converted into bound clumps. The main effect of turbulence is to thicken the mid-plane solid layer. Our results are thus consistent with planetesimal formation requiring a mid-plane solid-to-gas ratio $\epsilon \gtrsim 0.5$. We describe a method for tracking bound clumps in our simulations, and use this method to construct a mass function of planetesimals that is measured for each clump shortly after its collapse. Adopting this definition of the initial mass, instead of measuring masses at a fixed time, reduces planetesimal masses by a factor of three. For models in which planetesimals form, we show that the initial mass function is well-described by a broken power law, whose parameters are robust to the inclusion and strength of imposed turbulence. Turbulence in protoplanetary disks is likely to substantially exceed the threshold for planetesimal formation at radii where temperatures $T \gtrsim 10^3 \ {\rm K}$ lead to thermal ionization. Fully in situ planetesimal and planet formation may therefore not be viable for the closest-in exoplanets.
[2]  [pdf] - 2035774
Strongly magnetized accretion disks: structure and accretion from global magnetohydrodynamic simulations
Submitted: 2019-07-21, last modified: 2020-01-21
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 $\beta_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 $\beta_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).
[3]  [pdf] - 2015169
Simulation of a Compact Object with Outflows Moving Through a Gaseous Background
Comments: Submitted to MNRAS
Submitted: 2019-12-14
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 in some cases reverses its sign, 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. For some outflow configurations, the accretion rate is strongly intermittent. 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.
[4]  [pdf] - 1979784
The Influence of Black Hole Binarity on Tidal Disruption Events
Comments: 22 pages, 5 figures. Accepted chapter for Springer Space Science Reviews book on tidal disruption events
Submitted: 2019-09-27
Mergers are fundamental to the standard paradigm of galaxy evolution, and provide a natural formation mechanism for supermassive black hole binaries. The formation process of such a binary can have a direct impact on the rate at which stars are tidally disrupted by one or the other black hole, and the luminous signature of the tidal disruption itself can have distinct imprints of a binary companion. In this chapter we review our current understanding of the influence of black hole binarity on the properties of tidal disruption events. We discuss the rates of tidal disruption by supermassive black hole binaries, the impact of a second black hole on the fallback of debris and the formation of an accretion flow, and the prospects for detection of tidal disruption events by supermassive black hole binaries.
[5]  [pdf] - 1975191
The Mass and Size Distribution of Planetesimals Formed by the Streaming Instability. II. The Effect of the Radial Gas Pressure Gradient
Comments: 18 pages, 10 figures, 3 tables, accepted to ApJ
Submitted: 2018-10-23, last modified: 2019-08-29
The streaming instability concentrates solid particles in protoplanetary disks, leading to gravitational collapse into planetesimals. Despite its key role in producing particle clumping and determining critical length scales in the instability's linear regime, the influence of the disk's radial pressure gradient on planetesimal properties has not been examined in detail. Here, we use streaming instability simulations that include particle self-gravity to study how the planetesimal initial mass function depends on the radial pressure gradient. Fitting our results to a power-law, ${\rm d}N / {\rm d}M_p \propto M_p^{-p}$, we find a single value $p \approx 1.6$ describes simulations in which the pressure gradient varies by $\gtrsim 2$. An exponentially truncated power-law provides a significantly better fit, with a low mass slope of $p^\prime \approx 1.3$ that weakly depends on the pressure gradient. The characteristic truncation mass is found to be $\sim M_G = 4 \pi^5 G^2 \Sigma_p^3 / \Omega^4$. We exclude the cubic dependence of the characteristic mass with pressure gradient suggested by linear considerations, finding instead a linear scaling. These results strengthen the case for a streaming-derived initial mass function that depends at most weakly on the aerodynamic properties of the disk and participating solids. A simulation initialized with zero pressure gradient---which is {\em not} subject to the streaming instability---also yields a top-heavy mass function but with modest evidence for a different shape. We discuss the consistency of the theoretically predicted mass function with observations of Kuiper Belt planetesimals, and describe implications for models of early stage planet formation..
[6]  [pdf] - 1848654
Envisioning the next decade of Galactic Center science: a laboratory for the study of the physics and astrophysics of supermassive black holes
Comments: 9 pages, 4 figures, submitted for Astro2020 White Paper
Submitted: 2019-03-12
As the closest example of a galactic nucleus, the Galactic center (GC) presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environment. We describe several exciting new research directions that, over the next 10 years, hold the potential to answer some of the biggest scientific questions raised in recent decades: Is General Relativity (GR) the correct description for supermassive black holes? What is the nature of star formation in extreme environments? How do stars and compact objects dynamically interact with the supermassive black hole? What physical processes drive gas accretion in low-luminosity black holes? We describe how the high sensitivity, angular resolution, and astrometric precision offered by the next generation of large ground-based telescopes with adaptive optics will help us answer these questions. First, it will be possible to obtain precision measurements of stellar orbits in the Galaxy's central potential, providing both tests of GR in the unexplored regime near a SMBH and measurements of the extended dark matter distribution that is predicted to exist at the GC. Second, we will probe stellar populations at the GC to significantly lower masses than are possible today, down to brown dwarfs. Their structure and dynamics will provide an unprecedented view of the stellar cusp around the SMBH and will distinguish between models of star formation in this extreme environment. This increase in depth will also allow us to measure the currently unknown population of compact remnants at the GC by observing their effects on luminous sources. Third, uncertainties on the mass of and distance to the SMBH can be improved by a factor of $\sim$10. Finally, we can also study the near-infrared accretion onto the black hole at unprecedented sensitivity and time resolution, which can reveal the underlying physics of black hole accretion.
[7]  [pdf] - 1790659
A brief overview of planet formation
Comments: Short introduction, aimed at new students, for a forthcoming survey volume
Submitted: 2018-03-28
The initial conditions, physics, and outcome of planet formation are now constrained by detailed observations of protoplanetary disks, laboratory experiments, and the discovery of thousands of extrasolar planetary systems. These developments have broadened the range of processes that are considered important in planet formation, to include disk turbulence, radial drift, planet migration, and pervasive post-formation dynamical evolution. The N-body collisional growth of planetesimals and protoplanets, and the physics of planetary envelopes - key ingredients of the classical model - remain central. I provide an overview of the current status of planet formation theory, and discuss how it connects to observations.
[8]  [pdf] - 1659751
Interstellar object 'Oumuamua as an extinct fragment of an ejected cometary planetesimal
Comments: ApJ Letters, in press
Submitted: 2018-03-07
'Oumuamua was discovered passing through our Solar System on a hyperbolic orbit. It presents an apparent contradiction, with colors similar to those of volatile-rich Solar System bodies but with no visible outgassing or activity during its close approach to the Sun. Here we show that this contradiction can be explained by the dynamics of planetesimal ejection by giant planets. We propose that 'Oumuamua is an extinct fragment of a comet-like planetesimal born in a planet-forming disk that also formed Neptune- to Jupiter-mass giant planets. On its pathway to ejection 'Oumuamua's parent body underwent a close encounter with a giant planet and was tidally disrupted into small pieces, similar to comet Shoemaker-Levy 9's disruption after passing close to Jupiter. We use dynamical simulations to show that 0.1-1% of cometary planetesimals undergo disruptive encounters prior to ejection. Rocky asteroidal planetesimals are unlikely to disrupt due to their higher densities. After disruption, the bulk of fragments undergo enough close passages to their host stars to lose their surface volatiles and become extinct. Planetesimal fragments such as 'Oumuamua contain little of the mass in the population of interstellar objects but dominate by number. Our model makes predictions that will be tested in the coming decade by LSST.
[9]  [pdf] - 1656365
Circumbinary, not transitional: On the spiral arms, cavity, shadows, fast radial flows, streamers and horseshoe in the HD142527 disc
Comments: 16 pages, 12 figures, accepted to MNRAS. Movies at
Submitted: 2018-03-06
We present 3D hydrodynamical models of the HD142527 protoplanetary disc, a bright and well studied disc that shows spirals and shadows in scattered light around a 100 au gas cavity, a large horseshoe dust structure in mm continuum emission, together with mysterious fast radial flows and streamers seen in gas kinematics. By considering several possible orbits consistent with the observed arc, we show that all of the main observational features can be explained by one mechanism - the interaction between the disc and the observed binary companion. We find that the spirals, shadows and horseshoe are only produced in the correct position angles by a companion on an inclined and eccentric orbit approaching periastron - the 'red' family from Lacour et al. (2016). Dust-gas simulations show radial and azimuthal concentration of dust around the cavity, consistent with the observed horseshoe. The success of this model in the HD142527 disc suggests other mm-bright transition discs showing cavities, spirals and dust asymmetries may also be explained by the interaction with central companions.
[10]  [pdf] - 1648708
Implications of the interstellar object 1I/'Oumuamua for planetary dynamics and planetesimal formation
Comments: MNRAS, in press
Submitted: 2017-11-27, last modified: 2018-02-26
'Oumuamua, the first bona-fide interstellar planetesimal, was discovered passing through our Solar System on a hyperbolic orbit. This object was likely dynamically ejected from an extrasolar planetary system after a series of close encounters with gas giant planets. To account for 'Oumuamua's detection, simple arguments suggest that ~1 Earth mass of planetesimals are ejected per Solar mass of Galactic stars. However, that value assumes mono-sized planetesimals. If the planetesimal mass distribution is instead top-heavy the inferred mass in interstellar planetesimals increases to an implausibly high value. The tension between theoretical expectations for the planetesimal mass function and the observation of 'Oumuamua can be relieved if a small fraction (~0.1-1%) of planetesimals are tidally disrupted on the pathway to ejection into 'Oumuamua-sized fragments. Using a large suite of simulations of giant planet dynamics including planetesimals, we confirm that 0.1-1% of planetesimals pass within the tidal disruption radius of a gas giant on their pathway to ejection. 'Oumuamua may thus represent a surviving fragment of a disrupted planetesimal. Finally, we argue that an asteroidal composition is dynamically disfavoured for 'Oumuamua, as asteroidal planetesimals are both less abundant and ejected at a lower efficiency than cometary planetesimals.
[11]  [pdf] - 1608362
Tidal disruption by extreme mass ratio binaries and application to ASASSN-15lh
Comments: 10 pages, 3 figures; updated to include revisions from referee; fixed figures
Submitted: 2017-05-12, last modified: 2017-12-01
Tidal disruption events (TDEs) observed in massive galaxies with inferred central black hole masses $M_h > 10^8 \ M_\odot$ are presumptive candidates for TDEs by lower mass secondaries in binary systems. We use hydrodynamic simulations to quantify the characteristics of such TDEs, focusing on extreme mass ratio binaries and mpc separations where the debris stream samples the binary potential. The simulations are initialised with disruption trajectories from 3-body integrations of stars with parabolic orbits with respect to the binary center of mass. The most common outcome is found to be the formation of an unbound debris stream, with either weak late-time accretion or no accretion at all. A substantial fraction of streams remain bound, however, and these commonly yield structured fallback rate curves that exhibit multiple peaks or sharp drops. We apply our results to the superluminous supernova candidate ASASSN-15lh and show that its features, including its anomalous rebrightening at $\sim 100$ days after detection, are consistent with the tidal disruption of a star by a supermassive black hole in a binary system.
[12]  [pdf] - 1596320
Physical processes in protoplanetary disks
Comments: 150 pages, write up of lectures given at the 45th Saas-Fee Advanced Course "From Protoplanetary Disks to Planet Formation". Illustrations and movies at
Submitted: 2015-09-21, last modified: 2017-11-30
This review introduces physical processes in protoplanetary disks relevant to accretion and the initial stages of planet formation. After a brief overview of the observational context, I introduce the elementary theory of disk structure and evolution, review the gas-phase physics of angular momentum transport through turbulence and disk winds, and discuss possible origins for the episodic accretion observed in Young Stellar Objects. Turning to solids, I review the evolution of single particles under aerodynamic forces, and describe the conditions necessary for the development of collective gas-particle instabilities. Observations show that disks can exhibit pronounced large-scale structure, and I discuss the types of structures that may form from gas and particle interactions at ice lines, vortices and zonal flows, prior to the formation of large planetary bodies. I conclude with disk dispersal.
[13]  [pdf] - 1583307
Evidence for universality in the initial planetesimal mass function
Comments: 7 pages, 4 figures, accepted to ApJ Letters after minor modifications, including two new figures and some new text that better clarify our results
Submitted: 2017-05-10, last modified: 2017-09-13
Planetesimals may form from the gravitational collapse of dense particle clumps initiated by the streaming instability. We use simulations of aerodynamically coupled gas-particle mixtures to investigate whether the properties of planetesimals formed in this way depend upon the sizes of the particles that participate in the instability. Based on three high resolution simulations that span a range of dimensionless stopping time $6 \times 10^{-3} \leq \tau \leq 2$ no statistically significant differences in the initial planetesimal mass function are found. The mass functions are fit by a power-law, ${\rm d}N / {\rm d}M_p \propto M_p^{-p}$, with $p=1.5-1.7$ and errors of $\Delta p \approx 0.1$. Comparing the particle density fields prior to collapse, we find that the high wavenumber power spectra are similarly indistinguishable, though the large-scale geometry of structures induced via the streaming instability is significantly different between all three cases. We interpret the results as evidence for a near-universal slope to the mass function, arising from the small-scale structure of streaming-induced turbulence.
[14]  [pdf] - 1586079
Circumbinary discs from tidal disruption events
Comments: 5 pages, 2 figures, MNRAS Letters accepted
Submitted: 2017-07-18
Tidal disruption events, which occur when a star is shredded by the tidal field of a supermassive black hole, provide a means of fueling black hole accretion. Here we show, using a combination of three body orbit integrations and hydrodynamic simulations, that these events are also capable of generating circumbinary rings of gas around tight supermassive black hole binaries with small mass ratios. Depending on the thermodynamics, these rings can either fragment into clumps that orbit the binary, or evolve into a gaseous circumbinary disc. We argue that tidal disruptions provide a direct means of generating circumbinary discs around supermassive black hole binaries and, more generally, can replenish the reservoir of gas on very small scales in galactic nuclei.
[15]  [pdf] - 1584678
A Three-Dimensional View of Turbulence: Constraints on Turbulent Motions in the HD 163296 Protoplanetary Disk using DCO$^+$
Comments: Accepted to ApJ
Submitted: 2017-06-14
Gas kinematics are an important part of the planet formation process. Turbulence influences planetesimal growth and migration from the scale of sub-micron dust grains through gas-giant planets. Radio observations of resolved molecular line emission can directly measure this non-thermal motion and, taking advantage of the layered chemical structure of disks, different molecular lines can be combined to map the turbulence throughout the vertical extent of a protoplanetary disk. Here we present ALMA observations of three molecules (DCO$^+$(3-2), C$^{18}$O(2-1) and CO(2-1)) from the disk around HD 163296. We are able to place stringent upper limits ($v_{\rm turb}<$0.06c$_s$, $<$0.05c$_s$ and $<$0.04c$_s$ for CO(2-1), C$^{18}$O(2-1) and DCO$^+$(3-2) respectively), corresponding to $\alpha\lesssim$3$\times$10$^{-3}$, similar to our prior limit derived from CO(3-2). This indicates that there is little turbulence throughout the vertical extent of the disk, contrary to theoretical predictions based on the magneto-rotational instability and gravito-turbulence. In modeling the DCO$^+$ emission we also find that it is confined to three concentric rings at 65.7$\pm$0.9 au, 149.9$^{+0.5}_{-0.7}$ au and 259$\pm$1 au, indicative of a complex chemical environment.
[16]  [pdf] - 1536198
Lecture notes on the formation and early evolution of planetary systems
Comments: Based on lectures given at the University of Colorado, Boulder. Last revised February 2017. Published, in modified form, as "Astrophysics of Planet Formation" (Cambridge University Press)
Submitted: 2007-01-16, last modified: 2017-02-20
These notes provide an introduction to the theory of the formation and early evolution of planetary systems. Topics covered include the structure, evolution and dispersal of protoplanetary disks; the formation of planetesimals, terrestrial and gas giant planets; and orbital evolution due to gas disk migration, planetesimal scattering, planet-planet interactions, and tides.
[17]  [pdf] - 1533015
EXor outbursts from disk amplification of stellar magnetic cycles
Comments: ApJL, in press
Submitted: 2016-12-05
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.
[18]  [pdf] - 1531319
Tidal disruption events from supermassive black hole binaries
Comments: 31 pages, 25 figures, MNRAS Accepted; updated to incorporate referee comments. Movies of simulations can be found here:
Submitted: 2016-08-19, last modified: 2016-11-15
We investigate the pre-disruption gravitational dynamics and post-disruption hydrodynamics of the tidal disruption of stars by supermassive black hole (SMBH) binaries. We focus on binaries with relatively low mass primaries ($10^6M_{\odot}$), moderate mass ratios, and separations with reasonably long gravitational wave inspiral times (tens of Myr). First, we generate a large ensemble (between 1 and 10 million) of restricted three-body integrations to quantify the statistical properties of tidal disruptions by circular SMBH binaries of initially-unbound stars. Compared to the reference case of a disruption by a single SMBH, the binary potential induces significant variance into the specific energy and angular momentum of the star at the point of disruption. Second, we use Newtonian numerical hydrodynamics to study the detailed evolution of the fallback debris from 120 disruptions randomly selected from the three-body ensemble (excluding only the most deeply penetrating encounters). We find that the overall morphology of the debris is greatly altered by the presence of the second black hole, and the accretion rate histories display a wide range of behaviors, including order of magnitude dips and excesses relative to control simulations that include only one black hole. Complex evolution persists, in some cases, for many orbital periods of the binary. We find evidence for power in the accretion curves on timescales related to the binary orbital period, though there is no exact periodicity. We discuss our results in the context of future wide-field surveys, and comment on the prospects of identifying and characterizing the subset of events occurring in nuclei with binary SMBHs.
[19]  [pdf] - 1457534
Planet Formation Imager (PFI): science vision and key requirements
Comments: 13 pages, 5 figures, Proceedings of SPIE 2016
Submitted: 2016-08-01, last modified: 2016-08-16
The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.
[20]  [pdf] - 1470820
Prompt planetesimal formation beyond the snow line
Comments: ApJ Letters, in press
Submitted: 2016-08-11
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 $\alpha = 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.
[21]  [pdf] - 1444068
Turbulence, Transport and Waves in Ohmic Dead Zones
Comments: 14 pages, 17 figures
Submitted: 2016-06-09
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.
[22]  [pdf] - 1415082
Strongly magnetized accretion discs require poloidal flux
Comments: 7 pages, 3 figures, 1 table; Accepted to MNRAS
Submitted: 2016-02-15, last modified: 2016-05-18
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.
[23]  [pdf] - 1411224
The Mass and Size Distribution of Planetesimals Formed by the Streaming Instability. I. The Role of Self-Gravity
Comments: 18 pages, 16 figures, accepted to ApJ after minor revisions and title change
Submitted: 2015-11-30, last modified: 2016-03-31
We study the formation of planetesimals in protoplanetary disks from the gravitational collapse of solid over-densities generated via the streaming instability. To carry out these studies, we implement and test a particle-mesh self-gravity module for the Athena code that enables the simulation of aerodynamically coupled systems of gas and collisionless self-gravitating solid particles. Upon employment of our algorithm to planetesimal formation simulations, we find that (when a direct comparison is possible) the Athena simulations yield predicted planetesimal properties that agree well with those found in prior work using different numerical techniques. In particular, the gravitational collapse of streaming-initiated clumps leads to an initial planetesimal mass function that is well-represented by a power-law, dN/dM ~ M^{-p}, with p = 1.6 +/- 0.1, which equates to a differential size distribution dN/dR ~ R^{-q}, with q = 2.8 +/- 0.1. We find no significant trends with resolution from a convergence study of up to 512^3 grid zones and Npar = 1.5x10^8 particles. Likewise, the power-law slope appears indifferent to changes in the relative strength of self-gravity and tidal shear, and to the time when (for reasons of numerical economy) self-gravity is turned on, though the strength of these claims is limited by small number statistics. For a typically assumed radial distribution of minimum mass solar nebula solids (assumed here to have dimensionless stopping time \tau = 0.3), our results support the hypothesis that bodies on the scale of large asteroids or Kuiper Belt Objects could have formed as the high-mass tail of a primordial planetesimal population.
[24]  [pdf] - 1379068
Planet-disc evolution and the formation of Kozai-Lidov planets
Comments: Accepted for publication in MNRAS
Submitted: 2016-03-09
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).
[25]  [pdf] - 1392898
On the structure of tidally-disrupted stellar debris streams
Comments: 17 pages, 10 figures, Resubmitted to MNRAS after first referee report
Submitted: 2016-03-02
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.
[26]  [pdf] - 1351582
Accretion disc dynamo activity in local simulations spanning weak-to-strong net vertical magnetic flux regimes
Comments: 19 pages, 13 figures, 3 tables; Accepted to MNRAS
Submitted: 2015-11-19, last modified: 2016-01-12
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, $\beta_0^{\rm mid} = p_{\rm gas} / p_B$. For $10^5 \geq \beta_0^{\rm mid} \geq 10$ the effective $\alpha$-viscosity parameter scales as a power-law. Dynamo activity persists up to and including $\beta_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 $\beta^{\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.
[27]  [pdf] - 1323550
Post-periapsis pancakes: sustenance for self-gravity in tidal disruption events
Comments: 18 pages, 14 figures. Accepted for publication in MNRAS
Submitted: 2015-10-27
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 $\gamma \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.
[28]  [pdf] - 1288444
Magnetically driven accretion in protoplanetary discs
Comments: accepted to MNRAS after very minor revisions
Submitted: 2015-08-04, last modified: 2015-09-04
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.
[29]  [pdf] - 1264136
Accretion disk dynamo as the trigger for X-ray binary state transitions
Comments: ApJ, in press. Minor corrections to Figure 1 and section 4.4 to match proof version
Submitted: 2015-07-14, last modified: 2015-07-29
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.
[30]  [pdf] - 1263889
Signatures of MRI-Driven Turbulence in Protoplanetary Disks: Predictions for ALMA Observations
Comments: accepted to ApJ after minor revisions
Submitted: 2015-01-12, last modified: 2015-06-22
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.
[31]  [pdf] - 1284824
Accretion kinematics through the warped transition disk in HD142527 from resolved CO(6-5) observations
Comments: Submitted to ApJ
Submitted: 2015-05-28
The finding of residual gas in the large central cavity of the HD142527 disk motivates questions on the origin of its non-Keplerian kinematics, and possible connections with planet formation. We aim to understand the physical structure that underlies the intra-cavity gaseous flows, guided by new molecular-line data in CO(6-5) with unprecedented angular resolutions. Given the warped structure inferred from the identification of scattered-light shadows cast on the outer disk, the kinematics are consistent, to first order, with axisymmetric accretion onto the inner disk occurring at all azimuth. A steady-state accretion profile, fixed at the stellar accretion rate, explains the depth of the cavity as traced in CO isotopologues. The abrupt warp and evidence for near free-fall radial flows in HD 142527 resemble theoretical models for disk tearing, which could be driven by the reported low mass companion, whose orbit may be contained in the plane of the inner disk. The companion's high inclination with respect to the massive outer disk could drive Kozai oscillations over long time-scales; high-eccentricity periods may perhaps account for the large cavity. While shadowing by the tilted disk could imprint an azimuthal modulation in the molecular-line maps, further observations are required to ascertain the significance of azimuthal structure in the density field inside the cavity of HD142527.
[32]  [pdf] - 1216944
Importance of thermal diffusion in the gravo-magnetic limit cycle
Comments: 11 Pages, 11 Figures, Accepted for publication in MNRAS
Submitted: 2014-09-15
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.
[33]  [pdf] - 1216679
The Kozai-Lidov Mechanism in Hydrodynamical Disks
Comments: Accepted for publication in ApJ Letters
Submitted: 2014-09-03
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.
[34]  [pdf] - 1215816
Giant Outbursts in Be/X-ray Binaries
Comments: Accepted for publication in ApJ Letters
Submitted: 2014-07-21
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.
[35]  [pdf] - 1209158
Distinguishing an ejected blob from alternative flare models at the Galactic centre with GRAVITY
Comments: 11 pages, 9 figures, accepted by MNRAS
Submitted: 2014-04-24
The black hole at the Galactic centre exhibits regularly flares of radiation, the origin of which is still not understood. In this article, we study the ability of the near-future GRAVITY infrared instrument to constrain the nature of these events. We develop realistic simulations of GRAVITY astrometric data sets for various flare models. We show that the instrument will be able to distinguish an ejected blob from alternative flare models, provided the blob inclination is >= 45deg, the flare brightest magnitude is 14 <= mK <= 15 and the flare duration is >= 1h30.
[36]  [pdf] - 791968
Efficiency of Particle Trapping in the Outer Regions of Protoplanetary Disks
Comments: accepted for publication in ApJ, 8 pages, 5 figures
Submitted: 2014-02-06
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.
[37]  [pdf] - 1202912
A mechanism for hysteresis in black hole binary state transitions
Comments: 4 pages, 1 figure, accepted for publication in The Astrophysical Journal Letters
Submitted: 2014-01-21
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 $\alpha \sim 1$. The transition to a geometrically thin inner disk occurs when $L \sim \alpha^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, $\alpha$ decreases to $\alpha \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.
[38]  [pdf] - 1201919
Convergence of simulations of self-gravitating accretion discs II: Sensitivity to the implementation of radiative cooling and artificial viscosity
Comments: 11 pages, 7 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2013-11-28, last modified: 2014-01-13
Recently it has been suggested that the fragmentation boundary in Smoothed Particle Hydrodynamic (SPH) and FARGO simulations of self-gravitating accretion discs with beta-cooling do not converge as resolution is increased. Furthermore, this recent work suggests that by carefully optimising the artificial viscosity parameters in these codes it can be shown that fragmentation may occur for much longer cooling times than earlier work suggests. If correct, this result is intriguing as it suggests that gas giant planets could form, via direct gravitational collapse, reasonably close to their parent stars. This result is, however, slightly surprising and there have been a number of recent studies suggesting that the result is likely an indication of a numerical problem with the simulations. One suggestion, in particular, is that the SPH results are influenced by the manner in which the cooling is implemented. We extend this work here and show that if the cooling is implemented in a manner that removes a known numerical artefact in the shock regions, the fragmentation boundary converges to a value consistent with earlier work and that fragmentation is unlikely for the long cooling times suggested by this recent work. We also investigate the optimisation of the artificial viscosity parameters and show that the values that appear optimal are likely introducing numerical problems in both the SPH and FARGO simulations. We therefore conclude that earlier predictions for the cooling times required for fragmentation are likely correct and that, as suggested by this earlier work, fragmentation cannot occur in the inner parts (r < 50 au) of typical protostellar discs.
[39]  [pdf] - 1201881
The Long-Term Dynamical Evolution of Planetary Systems
Comments: 23 pages, 10 figures. Refereed review chapter, accepted for publication in Protostars & Planets VI, University of Arizona Press (2014), eds. H.Beuther, C.Dullemond, Th.Henning, R. Klessen
Submitted: 2013-11-26
This chapter concerns the long-term dynamical evolution of planetary systems from both theoretical and observational perspectives. We begin by discussing the planet-planet interactions that take place within our own Solar System. We then describe such interactions in more tightly-packed planetary systems. As planet-planet interactions build up, some systems become dynamically unstable, leading to strong encounters and ultimately either ejections or collisions of planets. After discussing the basic physical processes involved, we consider how these interactions apply to extrasolar planetary systems and explore the constraints provided by observed systems. The presence of a residual planetesimal disc can lead to planetary migration and hence cause instabilities induced by resonance crossing; however, such discs can also stabilise planetary systems. The crowded birth environment of a planetary system can have a significant impact: close encounters and binary companions can act to destabilise systems, or sculpt their properties. In the case of binaries, the Kozai mechanism can place planets on extremely eccentric orbits which may later circularise to produce hot Jupiters.
[40]  [pdf] - 1180576
The Dispersal of Protoplanetary Disks
Comments: 23 pages, 6 figures. Refereed review chapter, accepted for publication in Protostars & Planets VI, University of Arizona Press (2014), eds. H.Beuther, C.Dullemond, Th.Henning, R.Klessen
Submitted: 2013-11-07
Protoplanetary disks are the sites of planet formation, and the evolution and eventual dispersal of these disks strongly influences the formation of planetary systems. Disk evolution during the planet-forming epoch is driven by accretion and mass-loss due to winds, and in typical environments photoevaporation by high-energy radiation from the central star is likely to dominate final gas disk dispersal. We present a critical review of current theoretical models, and discuss the observations that are used to test these models and inform our understanding of the underlying physics. We also discuss the role disk dispersal plays in shaping planetary systems, considering its influence on both the process(es) of planet formation and the architectures of planetary systems. We conclude by presenting a schematic picture of protoplanetary disk evolution and dispersal, and discussing prospects for future work.
[41]  [pdf] - 1180218
Two timescale dispersal of magnetized protoplanetary disks
Comments: ApJL, in press
Submitted: 2013-10-24
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.
[42]  [pdf] - 1172018
Turbulence In the Outer Regions of Protoplanetary Disks. II. Strong Accretion Driven by a Vertical Magnetic Field
Comments: accepted to ApJ after very minor revisions
Submitted: 2013-06-13, last modified: 2013-08-12
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.
[43]  [pdf] - 1172211
Formation of circumbinary planets in a dead zone
Comments: Accepted for publication in ApJ
Submitted: 2013-06-21
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.
[44]  [pdf] - 1159016
The interplay between X-ray photoevaporation and planet formation
Comments: 11 pages, 10 figures, accepted by MNRAS on 31/12/2012
Submitted: 2013-01-14
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.
[45]  [pdf] - 1152141
Turbulence In the Outer Regions of Protoplanetary Disks. I. Weak Accretion with No Vertical Magnetic Flux
Comments: accepted to ApJ after minor revisions
Submitted: 2012-10-15, last modified: 2013-01-14
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.
[46]  [pdf] - 1157710
Mini-Oort clouds: Compact isotropic planetesimal clouds from planet-planet scattering
Comments: Accepted to MNRAS Letters. Associated movie available at
Submitted: 2012-11-12
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.
[47]  [pdf] - 1152349
A limit on eccentricity growth from global 3-D simulations of disc-planet interactions
Comments: 12 pages, 9 figures, 2 tables. Accepted for publication in MNRAS. Movies of some of the simulations can be found at:
Submitted: 2012-10-22
We present high resolution 3-D simulations of the planet-disc interaction using smoothed particle hydrodynamics, to investigate the possibility of driving eccentricity growth by this mechanism. For models with a given disc viscosity (\alpha = 0.01), we find that for small planet masses (a few Jupiter masses) and canonical surface densities, no significant eccentricity growth is seen over the duration of our simulations. This contrasts with the limiting case of large planet mass (over twenty Jupiter masses) and extremely high surface densities, where we find eccentricity growth in agreement with previously published results. We identify the cause of this as being a threshold surface density for a given planet mass below which eccentricity growth cannot be excited by this method. Further, the radial profile of the disc surface density is found to have a stronger effect on eccentricity growth than previously acknowledged, implying that care must be taken when contrasting results from different disc models. We discuss the implication of this result for real planets embedded in gaseous discs, and suggest that the disc-planet interaction does not contribute significantly to observed exoplanet eccentricities.
[48]  [pdf] - 1150797
Turbulent Linewidths as a Diagnostic of Self-Gravity in Protostellar Discs
Comments: 8 pages, 13 figures, accepted for publication in MNRAS
Submitted: 2012-08-21
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.
[49]  [pdf] - 1117286
Global variation of the dust-to-gas ratio in evolving protoplanetary discs
Comments: Accepted by Monthly Notices of the Royal Astronomical Society on the 8th of March 2012. 20 pages. 16 figures
Submitted: 2012-03-13
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.
[50]  [pdf] - 1117023
Emergent Mesoscale Phenomena in Magnetized Accretion Disc Turbulence
Comments: 18 pages, 20 figures, accepted to MNRAS
Submitted: 2012-03-01
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 \Delta 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 {\alpha} 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 {\alpha} = 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)
[51]  [pdf] - 1093045
Debris disks as signposts of terrestrial planet formation. II Dependence of exoplanet architectures on giant planet and disk properties
Comments: Astronomy and Astrophysics, in press. Movies from simulations are at
Submitted: 2012-01-17, last modified: 2012-02-16
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.
[52]  [pdf] - 1083635
The butterfly effect in the extreme-mass ratio inspiral problem
Comments: Accepted by ApJ Letts after minor modifications
Submitted: 2011-08-25, last modified: 2011-12-14
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.
[53]  [pdf] - 1091622
Convergence of SPH simulations of self-gravitating accretion discs: Sensitivity to the implementation of radiative cooling
Comments: 9 pages, 11 figures, MNRAS in press
Submitted: 2011-11-14
Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much more easily than would be expected from a two-dimensional local model. Subsequently, global and local two-dimensional models have been shown to display similar fragmentation properties, leaving it unclear whether the three-dimensional results reflect a physical effect or a numerical problem associated with the treatment of cooling or artificial viscosity in SPH. Here, we study how fragmentation of self-gravitating disc flows in SPH depends upon the implementation of cooling. We run disc simulations that compare a simple cooling scheme, in which each particle loses energy based upon its internal energy per unit mass, with a method in which the cooling is derived from a smoothed internal energy density field. For the simple per particle cooling scheme, we find a significant increase in the minimum cooling time scale for fragmentation with increasing resolution, matching previous results. Switching to smoothed cooling, however, results in lower critical cooling time scales, and tentative evidence for convergence at the highest spatial resolution tested. We conclude that precision studies of fragmentation using SPH require careful consideration of how cooling (and, probably, artificial viscosity) is implemented, and that the apparent non-convergence of the fragmentation boundary seen in prior simulations is likely a numerical effect. In real discs, where cooling is physically smoothed by radiative transfer effects, the fragmentation boundary is probably displaced from the two-dimensional value by a factor that is only of the order of unity.
[54]  [pdf] - 1078054
Turbulent Linewidths in Protoplanetary Disks: Predictions from Numerical Simulations
Comments: Accepted to ApJ after minor changes. A movie can be found at
Submitted: 2011-07-18, last modified: 2011-09-01
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.
[55]  [pdf] - 1083658
Hydrodynamic outcomes of planet scattering in transitional discs
Comments: MNRAS, in press. Movies from example simulations are available through the Data Conservancy, or at
Submitted: 2011-08-26
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.
[56]  [pdf] - 402403
Dynamics of Protoplanetary Disks
Comments: Annual Review of Astronomy and Astrophysics (2011). Final edited version at .High resolution versions of illustrations at
Submitted: 2010-11-05, last modified: 2011-08-23
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.
[57]  [pdf] - 1083225
Stability of self-gravitating discs under irradiation
Comments: MNRAS, in press
Submitted: 2011-08-04
Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress - which can be represented as an effective viscous alpha - out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.
[58]  [pdf] - 1076501
Turbulence in Global Simulations of Magnetized Thin Accretion Disks
Comments: 24 pages and 25 figures. MNRAS in press. Version with high resolution figures available from
Submitted: 2011-05-09
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 $< < \alpha > > \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.
[59]  [pdf] - 1053384
The debris disk - terrestrial planet connection
Comments: Contribution to proceedings of IAU 276: Astrophysics of Planetary Systems
Submitted: 2011-04-14
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.
[60]  [pdf] - 1053080
Debris disks as signposts of terrestrial planet formation
Comments: Astronomy & Astrophysics, in press
Submitted: 2011-03-31
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.
[61]  [pdf] - 1033640
Particle transport in evolving protoplanetary disks: Implications for results from Stardust
Comments: ApJ, in press
Submitted: 2010-07-12
Samples returned from comet 81P/Wild 2 by Stardust confirm that substantial quantities of crystalline silicates were incorporated into the comet at formation. We investigate the constraints that this observation places upon protoplanetary disk physics, assuming that outward transport of particles processed at high temperatures occurs via advection and turbulent diffusion in an evolving disk. We also look for constraints on particle formation locations. Our results are based upon 1D disk models that evolve with time under the action of viscosity and photoevaporation, and track solid transport using an ensemble of individual particle trajectories. We find that two classes of disk model are consistent with the Stardust findings. One class features a high particle diffusivity (a Schmidt number Sc < 1), which suffices to diffuse particles up to 20 microns in size outward against the mean gas flow. For Sc > 1, such models are unlikely to be viable, and significant outward transport requires that the particles of interest settle into a midplane layer that experiences an outward gas flow. In either class of models, the mass of inner disk material that reaches the outer disk is a strong function of the disk's initial compactness. Hence, models of grain transport within steady-state disks underestimate the efficiency of outward transport. Neither model results in sustained outward transport of very large particles exceeding a mm in size. We show that the transport efficiency generally falls off rapidly with time. Hence, high-temperature material must be rapidly incorporated into icy bodies to avoid fallback, and significant radial transport may only occur during the initial phase of rapid disk evolution. It may also vary substantially between disks depending upon their initial mass distributions. We discuss implications for Spitzer observations of crystalline silicates in T Tauri disks.
[62]  [pdf] - 1025273
Connections Between Local and Global Turbulence in Accretion Disks
Comments: 8 Pages, 7 Figures ApJ, In Press
Submitted: 2010-02-18
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.
[63]  [pdf] - 1024706
Planet-planet scattering in planetesimal disks II: Predictions for outer extrasolar planetary systems
Comments: 27 pages, 21 figures. ApJ in press. First paper in series at arXiv:0905.3741
Submitted: 2010-01-20
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.
[64]  [pdf] - 1018256
The role of disc self-gravity in the formation of protostars and protostellar discs
Comments: 11 pages, 8 figures, Accepted for publication in MNRAS
Submitted: 2009-11-06
We use time-dependent, one-dimensional disc models to investigate the evolution of protostellar discs that form through the collapse of molecular cloud cores and in which the primary transport mechanism is self-gravity. We assume that these discs settle into a state of thermal equilibrium with Q = 2 and that the strength of the angular momentum transport is set by the cooling rate of the disc. The results suggest that these discs will attain a quasi-steady state that persists for a number of free-fall times and in which most of the mass within 100 au is located inside 10-20 au. This pile-up of mass in the inner disc could result in temperatures that are high enough for the growth of MHD turbulence which could rapidly drain the inner disc and lead to FU Orionis-like outbursts. In all our simulations, the inner regions of the discs (r < 40 au) were stable against fragmentation, while fragmentation was possible in the outer regions (r > 40 au) of discs that formed from cores that had enough initial angular momentum to deposit sufficient mass in these outer regions. The large amounts of mass in these outer regions, however, suggests that fragmentation will lead to the formation of sub-stellar and stellar mass companions, rather than planetary mass objects. Although mass accretion rates were largely consistent with observations, the large disc masses suggest that an additional transport mechanism (such as MRI occuring in the upper layers of the disc) must operate in order to drain the remaining disc material within observed disc lifetimes.
[65]  [pdf] - 1017774
Black hole mergers: the first light
Comments: 16 pages, accepted by MNRAS. Animations of the simulations are available at
Submitted: 2009-10-01
The coalescence of supermassive black hole binaries occurs via the emission of gravitational waves, that can impart a substantial recoil to the merged black hole. We consider the energy dissipation, that results if the recoiling black hole is surrounded by a thin circumbinary disc. Our results differ significantly from those of previous investigations. We show analytically that the dominant source of energy is often potential energy, released as gas in the outer disc attempts to circularize at smaller radii. Thus, dimensional estimates, that include only the kinetic energy gained by the disc gas, underestimate the real energy loss. This underestimate can exceed an order of magnitude, if the recoil is directed close to the disc plane. We use three dimensional Smooth Particle Hydrodynamics (SPH) simulations and two dimensional finite difference simulations to verify our analytic estimates. We also compute the bolometric light curve, which is found to vary strongly depending upon the kick angle. A prompt emission signature due to this mechanism may be observable for low mass (10^6 Solar mass) black holes whose recoil velocities exceed about 1000 km/s. Emission at earlier times can mainly result from the response of the disc to the loss of mass, as the black holes merge. We derive analytically the condition for this to happen.
[66]  [pdf] - 901554
Giant planet migration, disk evolution, and the origin of transitional disks
Comments: 13 pages, 9 figures. Accepted for publication in ApJ
Submitted: 2009-08-31
We present models of giant planet migration in evolving protoplanetary disks. Our disks evolve subject to viscous transport of angular momentum and photoevaporation, while planets undergo Type II migration. We use a Monte Carlo approach, running large numbers of models with a range in initial conditions. We find that relatively simple models can reproduce both the observed radial distribution of extra-solar giant planets, and the lifetimes and accretion histories of protoplanetary disks. The use of state-of-the-art photoevaporation models results in a degree of coupling between planet formation and disk clearing, which has not been found previously. Some accretion across planetary orbits is necessary if planets are to survive at radii <~1.5AU, and if planets of Jupiter mass or greater are to survive in our models they must be able to form at late times, when the disk surface density in the formation region is low. Our model forms two different types of "transitional" disks, embedded planets and clearing disks, which show markedly different properties. We find that the observable properties of these systems are broadly consistent with current observations, and highlight useful observational diagnostics. We predict that young transition disks are more likely to contain embedded giant planets, while older transition disks are more likely to be undergoing disk clearing.
[67]  [pdf] - 24559
Planet-planet scattering in planetesimal disks
Comments: ApJ Letters, in press
Submitted: 2009-05-22
We study the final architecture of planetary systems that evolve under the combined effects of planet-planet and planetesimal scattering. Using N-body simulations we investigate the dynamics of marginally unstable systems of gas and ice giants both in isolation and when the planets form interior to a planetesimal belt. The unstable isolated systems evolve under planet-planet scattering to yield an eccentricity distribution that matches that observed for extrasolar planets. When planetesimals are included the outcome depends upon the total mass of the planets. For system masses exceeding about one Jupiter mass the final eccentricity distribution remains broad, whereas for lower mass planetary systems a combination of divergent orbital evolution and recircularization of scattered planets results in a preponderance of nearly circular final orbits. We also study the fate of marginally stable multiple planet systems in the presence of planetesimal disks, and find that for high planet masses the majority of such systems evolve into resonance. A significant fraction lead to resonant chains that are planetary analogs of Jupiter's Galilean satellites. We predict that a transition from eccentric to near-circular orbits will be observed once extrasolar planet surveys detect sub-Jovian mass planets at orbital radii of 5-10 AU.
[68]  [pdf] - 437677
Time-dependent models of the structure and evolution of self-gravitating protoplanetary discs
Comments: 10 pages, 10 figures, accepted for publication in MNRAS
Submitted: 2009-04-13
Angular momentum transport within young massive protoplanetary discs may be dominated by self-gravity at radii where the disk is too weakly ionized to allow the development of the magneto-rotational instability. We use time-dependent one-dimensional disc models, based on a local cooling time calculation of the efficiency of transport, to study the radial structure and stability (against fragmentation) of protoplanetary discs in which self-gravity is the sole transport mechanism. We find that self-gravitating discs rapidly attain a quasi-steady state in which the surface density in the inner disc is high and the strength of turbulence very low (alpha ~ 10^{-3} or less inside 5 au). Temperatures high enough to form crystalline silicates may extend out to several au at early times within these discs. None of our discs spontaneously develop regions that would be unambiguously unstable to fragmentation into substellar objects, though the outer regions (beyond 20 au) of the most massive discs are close enough to the threshold that fragmentation cannot be ruled out. We discuss how the mass accretion rates through such discs may vary with disc mass and with mass of the central star, and note that a determination of the \dot{M}-M_* relation for very young systems may allow a test of the model.
[69]  [pdf] - 315764
Planet-planet scattering leads to tightly packed planetary systems
Comments: Accepted to ApJ Letters. 5 pages, 2 figures, 2 tables,
Submitted: 2009-03-26
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.
[70]  [pdf] - 22210
A Proposal for a Renewed Research Emphasis in Astrophysical and Celestial Dynamics
Comments: 6 pp. White Paper submitted for the 2010 Decadal Survey, to panels PSF, SSE, GAN, & GCT
Submitted: 2009-03-09
Given the impressive investment by the nation in observational Astronomy and Astrophysics facilities coming on line now and in the near future, we advocate for an increased investment in applied and fundamental research on Astrophysical and Celestial Dynamics (ACD). Specifically we call for a) continued and expanded support for applied research in ACD, b) creation of support for fundamental research in ACD and its subfields, and c) the creation of a unified program to help scientists coordinate and collaborate in their research in these fields. The benefits of this proposal are threefold. First, it will enable researchers to interpret and understand the implications of newly observed phenomena that will invariably arise from new facilities and surveys. Second, research on fundamentals will foster connections between specialists, leveraging advances found in one sub-field and making them available to others. Third, a coordinated approach for applied and fundamental research in ACD will help academic institutions in the United States to produce future researchers trained and knowledgeable in essential subfields such as Mathematical Celestial Mechanics and able to continue its advancement in conjunction with the increase in observations.
[71]  [pdf] - 21961
The Galactic Center: A Laboratory for Fundamental Astrophysics and Galactic Nuclei
Comments: An ASTRO2010 White Paper. Animations and high resolution images can be found at
Submitted: 2009-03-02
As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH.
[72]  [pdf] - 15865
Massive black hole binary mergers within sub-pc scale gas discs
Comments: 11 pages, 13 figures. MNRAS, accepted. v2: minor changes
Submitted: 2008-09-01, last modified: 2008-11-03
[ABRIDGED] We study supermassive black hole binary mergers driven by angular momentum loss to small-scale gas discs. Such binaries form after major galaxy mergers, but their fate is unclear since hardening through stellar scattering becomes very inefficient at sub-parsec distances. Gas discs may dominate binary dynamics on these scales, and promote mergers. Using numerical simulations, we investigate the evolution of the orbits of binaries embedded within geometrically thin gas discs. Our simulations directly resolve angular momentum transport within the disc, which at the radii of interest is likely dominated by disc self-gravity. We show that the binary decays at a rate which is in good agreement with analytical estimates, while the eccentricity grows. Saturation of eccentricity growth is not observed up to values e > 0.35. Accretion onto the black holes is variable, and is roughly modulated by the binary orbital frequency. Scaling our results, we analytically estimate the maximum rate of binary decay that is possible without fragmentation occuring within the surrounding gas disc, and compare that rate to an estimate of the stellar dynamical hardening rate. For binary masses in the range 10^5 < M/Msun < 10^8 we find that decay due to gas discs may dominate for separations below a ~ 0.01--0.1 pc, where the disc is optically thick. The minimum merger time scale is shorter than the Hubble time for M < 10^7 Msun. This implies that gas discs could commonly attend relatively low mass black hole mergers, and that a significant population of binaries might exist at separations of a few hundredths of a pc, where the combined decay rate is slowest. For more massive binaries, we suggest that scattering of stars formed within a fragmenting gas disc could act as a significant additional sink of binary angular momentum.
[73]  [pdf] - 16539
Mean motion resonances from planet-planet scattering
Comments: 14 pages, accepted to ApJ Letters
Submitted: 2008-09-19
Planet-planet scattering is the leading mechanism to explain the large eccentricities of the observed exoplanet population. However, scattering has not been considered important to the production of pairs of planets in mean motion resonances (MMRs). We present results from a large number of numerical simulations of dynamical instabilities in 3-planet systems. We show that MMRs arise naturally in about five percent of cases. The most common resonances we populate are the 2:1 and 3:1 MMRs, although a wide variety of MMRs can occur, including high-order MMRs (up to eleventh order). MMRs are generated preferentially in systems with uneven mass distributions: the smallest planet is typically ejected after a series of close encounters, leaving the remaining, more massive planets in resonance. The distribution of resonant planets is consistent with the phase-space density of resonant orbits, meaning that planets are randomly thrown into MMRs rather than being slowly pulled into them. It may be possible to distinguish between MMRs created by scattering vs. convergent migration in a gaseous disk by considering planetary mass ratios: resonant pairs of planets beyond ~1 AU with more massive outer planets are likely to have formed by scattering. In addition, scattering may be responsible for pairs of planets in high-order MMRs (3:1 and higher) that are not easily populated by migration. The frequency of MMRs from scattering is comparable to the expected survival rate of MMRs in turbulent disks. Thus, planet-planet scattering is likely to be a major contributor to the population of resonant planets.
[74]  [pdf] - 14748
Microphysical dissipation, turbulence and magnetic fields in hyper-accreting discs
Comments: 14 pages, 10 figures. MNRAS, accepted. Minor revisions. Comments added in discussion section
Submitted: 2008-07-22, last modified: 2008-09-17
Hyper-accreting discs occur in compact-object mergers and collapsars, and may power gamma-ray bursts (GRBs). We calculate the microscopic viscosity and resistivity of plasma in these discs, and discuss the implications for their global structure and evolution. In the neutrino-cooled innermost regions, the viscosity is provided mainly by mildly degenerate electrons, while the resistivity is modified from the Spitzer value due to the effects of both relativity and degeneracy. The plasma behaves as an almost ideal MHD fluid. Among the non-ideal MHD effects the Hall term is relatively the most important, while the magnetic Prandtl number, Pr (the ratio of viscosity to resistivity), is typically larger than unity: 10 < Pr < 6000. The outer radiatively inefficient regions also display high Pr. Numerical simulations of the magneto-rotational instability indicate that the saturation level and angular momentum transport efficiency may be greatly enhanced at high Pr. If this behaviour persists in the presence of a strong Hall effect we would expect that hyper-accreting discs should be strongly magnetised and highly variable. The expulsion of magnetic field that cannot be dissipated at small scales may also favour a magnetic outflow. We note the similaries between the Prandtl number in hyper-accreting discs and X-ray binary discs, which suggests that a comparison between late-time activity in GRBs and X-ray binary accretion states may be fruitful. Our results imply that the behavior of high Prandtl number MHD flows needs to be considered in studies of hyper-accreting discs.
[75]  [pdf] - 14876
Extrasolar Planet Eccentricities from Scattering in the Presence of Residual Gas Disks
Comments: 7 pages, 7 figures. Accepted to ApJ
Submitted: 2008-07-25
Gravitational scattering between massive planets has been invoked to explain the eccentricity distribution of extrasolar planets. For scattering to occur, the planets must either form in -- or migrate into -- an unstable configuration. In either case, it is likely that a residual gas disk, with a mass comparable to that of the planets, will be present when scattering occurs. Using explicit hydrodynamic simulations, we study the impact of gas disks on the outcome of two-planet scattering. We assume a specific model in which the planets are driven toward instability by gravitational torques from an outer low mass disk. We find that the accretion of mass and angular momentum that occurs when a scattered planet impacts the disk can strongly influence the subsequent dynamics by reducing the number of close encounters. The eccentricity of the innermost surviving planet at the epoch when the system becomes Hill stable is not substantially altered from the gas-free case, but the outer planet is circularized by its interaction with the disk. The signature of scattering initiated by gas disk migration is thus a high fraction of low eccentricity planets at larger radii accompanying known eccentric planets. Subsequent secular evolution of the two planets in the presence of damping can further damp both eccentricities, and tends to push systems away from apsidal alignment and toward anti-alignment. We note that the late burst of accretion when the outer planet impacts the disk is in principle observable, probably via detection of a strong near-IR excess in systems with otherwise weak disk and stellar accretion signatures.
[76]  [pdf] - 14410
Binary formation and mass function variations in fragmenting discs with short cooling times
Comments: 10 pages, 5 figures. Accepted for publication in MNRAS. Figures 1 and 3 degraded to meet arXiv size limits - version with high resolution figures available at
Submitted: 2008-07-10
Accretion discs at sub-pc distances around supermassive black holes are likely to cool rapidly enough that self-gravity results in fragmentation. Here, we use high-resolution hydrodynamic simulations of a simplified disc model to study how the outcome of fragmentation depends upon numerical resolution and cooling time, and to investigate the incidence of binary formation within fragmenting discs. We investigate a range of cooling times, from the relatively long cooling time-scales that are marginally unstable to fragmentation down to highly unstable cooling on a time-scale that is shorter than the local dynamical time. The characteristic mass of fragments decreases with reduced cooling time, though the effect is modest and dependent upon details of how rapidly bound clumps radiate. We observe a high incidence of capture binaries, though we are unable to determine their final orbits or probability of survival. The results suggest that faster cooling in the parent disc results in an increased binary fraction, and that a high primordial binary fraction may result from disc fragmentation. We discuss our results in terms of the young massive stars close to the Galactic Centre, and suggest that observations of some stellar binaries close to the Galactic Centre remain consistent with formation in a fragmenting accretion disc.
[77]  [pdf] - 13598
Interaction of the magnetorotational instability with hydrodynamic turbulence in accretion disks
Comments: ApJ, in press
Submitted: 2008-06-16
Accretion disks in which angular momentum transport is dominated by the magnetorotational instability (MRI) can also possess additional, purely hydrodynamic, drivers of turbulence. Even when the hydrodynamic processes, on their own, generate negligible levels of transport, they may still affect the evolution of the disk via their influence on the MRI. Here, we study the interaction between the MRI and hydrodynamic turbulence using local MRI simulations that include hydrodynamic forcing. As expected, we find that hydrodynamic forcing is generally negligible if it yields a saturated kinetic energy density that is small compared to the value generated by the MRI. For stronger hydrodynamic forcing levels, we find that hydrodynamic turbulence modifies transport, with the effect varying depending upon the spatial scale of hydrodynamic driving. Large scale forcing boosts transport by an amount that is approximately linear in the forcing strength, and leaves the character of the MRI (for example the ratio between Maxwell and Reynolds stresses) unchanged, up to the point at which the forced turbulence is an order of magnitude stronger than that generated by the MRI. Low amplitude small scale forcing may modestly suppress the MRI. We conclude that the impact of hydrodynamic turbulence on the MRI is generically ignorable in cases, such as convection, where the additional turbulence arises due to the accretion energy liberated by the MRI itself. Hydrodynamic turbulence may affect (and either enhance or suppress) the MRI if it is both strong, and driven by independent mechanisms such as self-gravity, supernovae, or solid-gas interactions in multiphase protoplanetary disks.
[78]  [pdf] - 12003
Stellar dynamical evidence against a cold disc origin for stars in the Galactic Centre
Comments: 5 pages, 6 colour figures. MNRAS Letters in press. (v2: very minor changes)
Submitted: 2008-04-22, last modified: 2008-05-30
Observations of massive stars within the central parsec of the Galaxy show that, while most stars orbit within a well-defined disc, a significant fraction have large eccentricities and / or inclinations with respect to the disc plane. Here, we investigate whether this dynamically hot component could have arisen via scattering from an initially cold disc -- the expected initial condition if the stars formed from the fragmentation of an accretion disc. Using N-body methods, we evolve a variety of flat, cold, stellar systems, and study the effects of initial disc eccentricity, primordial binaries, very massive stars and intermediate mass black holes. We find, consistent with previous results, that a circular disc does not become eccentric enough unless there is a significant population of undetected 100--1000 Msun objects. However, since fragmentation of an eccentric disc can readily yield eccentric stellar orbits, the strongest constraints come from inclinations. We show that_none_ of our initial conditions yield the observed large inclinations, regardless of the initial disc eccentricity or the presence of massive objects. These results imply that the orbits of the young massive stars in the Galactic Centre are largely primordial, and that the stars are unlikely to have formed as a dynamically cold disc.
[79]  [pdf] - 7389
Quasistars: Accreting black holes inside massive envelopes
Comments: MNRAS, in press. Minor changes from previous version including additional discussion of stability
Submitted: 2007-11-26, last modified: 2008-04-16
We study the structure and evolution of "quasistars", accreting black holes embedded within massive hydrostatic gaseous envelopes. These configurations may model the early growth of supermassive black hole seeds. The accretion rate onto the black hole adjusts so that the luminosity carried by the convective envelope equals the Eddington limit for the total mass. This greatly exceeds the Eddington limit for the black hole mass alone, leading to rapid growth of the black hole. We use analytic models and numerical stellar structure calculations to study the structure and evolution of quasistars. We derive analytically the scaling of the photospheric temperature with the black hole mass and envelope mass, and show that it decreases with time as the black hole mass increases. Once the photospheric temperature becomes lower than 10,000 K, the photospheric opacity drops precipitously and the photospheric temperature hits a limiting value, analogous to the Hayashi track for red giants and protostars, below which no hydrostatic solution for the convective envelope exists. For metal-free (Population III) opacities this limiting temperature is approximately 4000 K. After a quasistar reaches this limiting temperature, the envelope is rapidly dispersed by radiation pressure. We find that black hole seeds with masses between 1000 and 10000 solar masses could form via this mechanism in less than a few Myr.
[80]  [pdf] - 9953
Eccentricity of masing disks in Active Galactic Nuclei
Comments: ApJ, submitted
Submitted: 2008-02-11
Observations of Keplerian disks of masers in NCG 4258 and other Seyfert galaxies can be used to obtain geometric distance estimates and derive the Hubble constant. The ultimate precision of such measurements could be limited by uncertainties in the disk geometry. Using a time-dependent linear theory model, we study the evolution of a thin initially eccentric disk under conditions appropriate to sub-pc scales in Active Galactic Nuclei. The evolution is controlled by a combination of differential precession driven by the disk potential and propagating eccentricity waves that are damped by viscosity. A simple estimate yields a circularization timescale of approximately 10 Myr at 0.1 pc. Numerical solutions for the eccentricity evolution confirm that damping commences on this timescale, but show that the subsequent decay rate of the eccentricity depends upon the uncertain strength of viscous damping of eccentricity. If eccentricity waves are important further decay of the eccentricity can be slow, with full circularization requiring up to 50 Myr for disks at radii of 0.1 pc to 0.2 pc. Observationally, this implies that it is plausible that enough time has elapsed for the eccentricity of masing disks to have been substantially damped, but that it may not be justified to assume vanishing eccentricity. We predict that during the damping phase the pericenter of the eccentric orbits describes a moderately tightly wound spiral with radius.
[81]  [pdf] - 7785
Why are there so few hot Jupiters?
Comments: 7 pages, 8 figures, Accepted for publication in MNRAS
Submitted: 2007-12-05
We use numerical simulations to model the migration of massive planets at small radii and compare the results with the known properties of 'hot Jupiters' (extrasolar planets with semi-major axes a < 0.1 AU). For planet masses Mp sin i > 0.5 MJup, the evidence for any `pile-up' at small radii is weak (statistically insignificant), and although the mass function of hot Jupiters is deficient in high mass planets as compared to a reference sample located further out, the small sample size precludes definitive conclusions. We suggest that these properties are consistent with disc migration followed by entry into a magnetospheric cavity close to the star. Entry into the cavity results in a slowing of migration, accompanied by a growth in orbital eccentricity. For planet masses in excess of 1 Jupiter mass we find eccentricity growth timescales of a few x 10^5 years, suggesting that these planets may often be rapidly destroyed. Eccentricity growth appears to be faster for more massive planets which may explain changes in the planetary mass function at small radii and may also predict a pile-up of lower mass planets, the sample of which is still incomplete.
[82]  [pdf] - 314899
Self-gravitating fragmentation of eccentric accretion disks
Comments: 9 pages, 7 figures. Accepted for publication in ApJ
Submitted: 2007-11-05
We consider the effects of eccentricity on the fragmentation of gravitationally unstable accretion disks, using numerical hydrodynamics. We find that eccentricity does not affect the overall stability of the disk against fragmentation, but significantly alters the manner in which such fragments accrete gas. Variable tidal forces around an eccentric orbit slow the accretion process, and suppress the formation of weakly-bound clumps. The "stellar" mass function resulting from the fragmentation of an eccentric disk is found to have a significantly higher characteristic mass than that from a corresponding circular disk. We discuss our results in terms of the disk(s) of massive stars at ~0.1pc from the Galactic Center, and find that the fragmentation of an eccentric accretion disk, due to gravitational instability, is a viable mechanism for the formation of these systems.
[83]  [pdf] - 1467
Massive planet migration: Theoretical predictions and comparison with observations
Comments: ApJ, in press. References updated to match published version
Submitted: 2007-05-21, last modified: 2007-07-25
We quantify the utility of large radial velocity surveys for constraining theoretical models of Type II migration and protoplanetary disk physics. We describe a theoretical model for the expected radial distribution of extrasolar planets that combines an analytic description of migration with an empirically calibrated disk model. The disk model includes viscous evolution and mass loss via photoevaporation. Comparing the predicted distribution to a uniformly selected subsample of planets from the Lick / Keck / AAT planet search programs, we find that a simple model in which planets form in the outer disk at a uniform rate, migrate inward according to a standard Type II prescription, and become stranded when the gas disk is dispersed, is consistent with the radial distribution of planets for orbital radii 0.1 AU < a < 2.5 AU and planet masses greater than 1.65 Jupiter masses. Some variant models are disfavored by existing data, but the significance is limited (~95%) due to the small sample of planets suitable for statistical analysis. We show that the favored model predicts that the planetary mass function should be almost independent of orbital radius at distances where migration dominates the massive planet population. We also study how the radial distribution of planets depends upon the adopted disk model. We find that the distribution can constrain not only changes in the power-law index of the disk viscosity, but also sharp jumps in the efficiency of angular momentum transport that might occur at small radii.
[84]  [pdf] - 1420
Role of Dynamical Research in the Detection and Characterization of Exoplanets
Comments: 9 pages, White Paper Submitted to the ExoPlanet Task Force
Submitted: 2007-05-18
(Abridged) The discovery of extrasolar planetary systems revealed an unexpected diversity of planetary systems that has revolutionized planet formation theory. A strong program of theoretical research is essential to maximize both the discovery potential and the scientific returns of future observational programs, so as to achieve a deeper understanding of the formation and evolution of planetary systems. We outline three broad categories of theoretical research: detailed studies of specific planetary systems, testing planet formation models by comparing their predictions to the observed exoplanet population, and detailed modeling of specific physical processes. We describe how such theoretical research plays an important role in analyzing observations for a wide range detection methods and contributes to understanding the Earth's place in the universe and the potential for Earth-like life beyond our solar system. In this white paper, we suggest how to maximize the scientific return of future exoplanet observations. Our recommendations include a strong theory program, support for multiple observational programs that will study a diverse set of planets and stars, significant observing time devoted to follow-up observations, and healthy collaboration between observers and theorists.
[85]  [pdf] - 316906
Extrasolar planetary dynamics with a generalized planar Laplace-Lagrange secular theory
Comments: 14 pages, 12 figures, 1 table, accepted for publication in ApJ
Submitted: 2007-03-05
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.
[86]  [pdf] - 87207
Dust dynamics during protoplanetary disc clearing
Comments: 13 pages, 5 figures. Accepted for publication in MNRAS
Submitted: 2006-11-27
We consider the dynamics of dust and gas during the clearing of protoplanetary discs. We work within the context of a photoevaporation/viscous model for the evolution of the gas disc, and use a two-fluid model to study the dynamics of dust grains as the gas disc is cleared. Small (<~10um) grains remain well-coupled to the gas, but larger (~1mm) grains are subject to inward migration from large radii (~50AU), suggesting that the time-scale for grain growth in the outer disc is ~10^4-10^5yr. We describe in detail the observable appearance of discs during clearing, and find that pressure gradients in the gas disc result in a strong enhancement of the local dust-to-gas ratio in a ring near to the inner disc edge. Lastly, we consider a simple model of the disc-planet interaction, and suggest that observations of disc masses and accretion rates provide a straightforward means of discriminating between different models of disc clearing.
[87]  [pdf] - 85400
Dust filtration at gap edges: Implications for the spectral energy distributions of discs with embedded planets
Comments: 9 pages, 7 figures, Accepted fir publication in MNRAS
Submitted: 2006-09-29
The spectral energy distributions (SEDs) of some T Tauri stars display a deficit of near-IR flux that could be a consequence of an embedded Jupiter-mass planet partially clearing an inner hole in the circumstellar disc. Here, we use two-dimensional numerical simulations of the planet-disc interaction, in concert with simple models for the dust dynamics, to quantify how a planet influences the dust at different radii within the disc. We show that pressure gradients at the outer edge of the gap cleared by the planet act as a filter - letting particles smaller than a critical size through to the inner disc while holding back larger particles in the outer disc. The critical particle size depends upon the disc properties, but is typically of the order of 10 microns. This filtration process will lead to discontinuous grain populations across the planet's orbital radius, with small grains in the inner disc and an outer population of larger grains. We show that this type of dust population is qualitatively consistent with SED modelling of systems that have optically thin inner holes in their circumstellar discs. This process can also produce a very large gas-to-dust ratio in the inner disc, potentially explaining those systems with optically thin inner cavities that still have relatively high accretion rates.
[88]  [pdf] - 316560
Constraints on the Stellar Mass Function from Stellar Dynamics at the Galactic Center
Comments: 9 pages, 5 figures. Accepted for publication in ApJ
Submitted: 2006-09-29
We consider the dynamical evolution of a disk of stars orbiting a central black hole. In particular, we focus on the effect of the stellar mass function on the evolution of the disk, using both analytic arguments and numerical simulations. We apply our model to the ring of massive stars at ~0.1pc from the Galactic Center, assuming that the stars formed in a cold, circular disk, and find that our model requires the presence of a significant population of massive (>100Msun) stars in order to explain the the observed eccentricities of 0.2-0.3. Moreover, in order to limit the damping of the heavier stars' eccentricities, we also require fewer low-mass stars than expected from a Salpeter mass function, giving strong evidence for a significantly ``top-heavy'' mass function in the rings of stars seen near to the Galactic Center. We also note that the maximum possible eccentricities attainable from circular initial conditions at ages of <10Myr are around 0.4-0.5, and suggest that any rings of stars found with higher eccentricities were probably not formed from circular disks.
[89]  [pdf] - 83471
Planetesimal formation via fragmentation in self-gravitating protoplanetary discs
Comments: 5 pages, Accepted for publication in MNRAS
Submitted: 2006-07-12
An unsolved issue in the standard core accretion model for gaseous planet formation is how kilometre-sized planetesimals form from, initially, micron-sized dust grains. Solid growth beyond metre sizes can be difficult both because the sticking efficiency becomes very small, and because these particles should rapidly migrate into the central star. We consider here how metre-sized particles evolve in self-gravitating accretion discs using simulations in which the gravitational influence of the solid particles is also included. Metre-sized particles become strongly concentrated in the spiral structures present in the disc and, if the solid to gas density ratio is sufficiently high, can fragment due to their own self-gravity to form planetesimals directly. This result suggests that planetesimal formation may occur very early in the star formation process while discs are still massive enough to be self-gravitating. The dependence of this process on the surface density of the solids is also consistent with the observation that extrasolar planets are preferentially found around high metallicity stars.
[90]  [pdf] - 81146
Predictions for the correlation between giant and terrestrial extrasolar planets in dynamically evolved systems
Comments: 8 pages, 8 figures, 1 table, accepted for publication in ApJ
Submitted: 2006-04-04
The large eccentricities of many giant extrasolar planets may represent the endpoint of gravitational scattering in initially more crowded systems. If so, the early evolution of the giant planets is likely to be more restrictive of terrestrial planet formation than would be inferred from the current, dynamically quiescent, configurations. Here, we study statistically the extent of the anti-correlation between giant planets and terrestrial planets expected in a scattering model. We use marginally stable systems of three giant planets, with a realistic range of planetary masses, as a simple model for the initial conditions prior to scattering, and show that after scattering the surviving planets reproduce well the known extrasolar planet eccentricities beyond a > 0.5 AU. By tracking the minimum periastron values of all planets during the evolution, we derive the distribution of orbital radii across which strong perturbations (from crossing orbits) are likely to affect low mass planet formation. We find that scattering affects inner planet formation at orbital separations less than 50% of the final periastron distance of the innermost massive planet in approximately 30% of the realizations, and can occasionally influence planet formation at orbital separations less than 20% of the final periastron distance of the innermost massive planet. The domain of influence of the scattering massive planets increases as the mass differential between the massive planets decreases. Observational study of the correlation between massive and terrestrial extrasolar planets in the same system has the potential to constrain the origin of planetary eccentricity.
[91]  [pdf] - 79635
The stellar mass-accretion rate relation in T Tauri stars and brown dwarfs
Comments: 4 pages, 2 figures. Accepted for publication in ApJ Letters
Submitted: 2006-02-02
Recent observations show a strong correlation between stellar mass and accretion rate in young stellar and sub-stellar objects, with the scaling $\dot{M}_{acc} \propto M_*^2$ holding over more than four orders of magnitude in accretion rate. We explore the consequences of this correlation in the context of disk evolution models. We note that such a correlation is not expected to arise from variations in disk angular momentum transport efficiency with stellar mass, and suggest that it may reflect a systematic trend in disk initial conditions. In this case we find that brown dwarf disks initially have rather larger radii than those around more massive objects. By considering disk evolution, and invoking a simple parametrization for a shut-off in accretion at the end of the disk lifetime, we show that such models predict that the scatter in the stellar mass-accretion rate relationship should increase with increasing stellar mass, in rough agreement with current observations.
[92]  [pdf] - 77864
Flares in long and short Gamma-Ray Bursts: a common origin in a hyperaccreting accretion disk
Comments: Accepted to the Astrophysical Journal Letters
Submitted: 2005-11-16
Early-time X-ray observations of GRBs with the Swift satellite have revealed a more complicated phenomenology than was known before. In particular, the presence of flaring activity on a wide range of time scales probably requires late-time energy production within the GRB engine. Since the flaring activity is observed in both long and short GRBs, its origin must be within what is in common for the two likely progenitors of the two classes of bursts: a hyperaccreting accretion disk around a black hole of a few solar masses. Here, we show that some of the observational properties of the flares, such as the duration-time scale correlation, and the duration-peak luminosity anticorrelation displayed by most flares within a given burst, are qualitatively consistent with viscous disk evolution, provided that the disk at large radii either fragments or otherwise suffers large amplitude variability. We discuss the physical conditions in the outer parts of the disk, and conclude that gravitational instability, possibly followed by fragmentation, is the most likely candidate for this variability.
[93]  [pdf] - 75995
Investigating fragmentation conditions in self-gravitating accretion discs
Comments: 5 pages, MNRAS, accepted
Submitted: 2005-09-14
The issue of fragmentation in self-gravitating gaseous accretion discs has implications both for the formation of stars in discs in the nuclei of active galaxies, and for the formation of gaseous planets or brown dwarfs in circumstellar discs. It is now well established that fragmentation occurs if the disc is cooled on a timescale smaller than the local dynamical timescale, while for longer cooling times the disc reaches a quasi-steady state in thermal equilibrium, with the cooling rate balanced by the heating due to gravitational stresses. We investigate here how the fragmentation boundary depends on the assumed equation of state. We find that the cooling time required for fragmentation increases as the specific heat ratio, gamma, decreases, exceeding the local dynamical timescale for gamma = 7/5. This result can be easily interpreted as a consequence of there being a maximum stress (in units of the local disc pressure) that can be sustained by a self-gravitating disc in quasi-equilibrium. Fragmentation occurs if the cooling time is such that the stress required to reach thermal equilibrium exceeds this value, independent of gamma. This result suggest that a quasi-steady, self-gravitating disc can never produce a stress that results in the viscous alpha parameter exceeding ~ 0.06.
[94]  [pdf] - 74661
Planetary migration
Comments: Revised version corrects summary of simulation results on co-orbital torques. To appear in, "A Decade Of Extrasolar Planets Around Normal Stars", STScI May Symposium 2005. Accompanying animation of migration regimes at
Submitted: 2005-07-20, last modified: 2005-09-07
Gravitational torques between a planet and gas in the protoplanetary disk result in orbital migration of the planet, and are likely to play an important role in the formation and early evolution of planetary systems. For masses comparable to those of observed giant extrasolar planets, the interaction with the disk is strong enough to form a gap, leading to coupled evolution of the planet and disk on a viscous time scale (Type II migration). Both the existence of hot Jupiters, and the statistical distribution of observed orbital radii, are consistent with an important role for Type II migration in the history of currently observed systems. We discuss the possibility of improving constraints on migration by including information on the host stars' metallicity, and note that migration could also form a population of massive planets at large orbital radii that may be indirectly detected via their influence on debris disks. For lower mass planets with masses of the order of that of the Earth, surface density perturbations created by the planet are small, and migration in a laminar disk is driven by an intrinsic and apparently robust asymmetry between interior and exterior torques. Analytic and numerical calculations of this Type I migration are in reasonable accord, and predict rapid orbital decay during the final stages of the formation of giant planet cores. The difficulty of reconciling Type I migration with giant planet formation may signal basic errors in our understanding of protoplanetary disks, core accretion, or both. We discuss physical effects that might alter Type I behavior, in particular the possibility that for sufficiently low masses turbulent fluctuations in the gas surface density dominate the torque, leading to random walk migration of very low mass bodies.
[95]  [pdf] - 75382
Eccentricity of supermassive black hole binaries coalescing from gas rich mergers
Comments: ApJ, in press
Submitted: 2005-08-23
Angular momentum loss to circumbinary gas provides a possible mechanism for overcoming the `last parsec' problem and allowing the most massive black hole binaries formed from galactic mergers to coalesce. Here, we show that if gas disks also catalyze the merger of the somewhat lower mass binaries detectable with the Laser Interferometer Space Antenna (LISA), then there may be a purely gravitational wave signature of the role of gas in the form of a small but finite eccentricity just prior to merger. Numerical simulations suggest that eccentricity, excited by the interaction between the binary and surrounding gas disk, is only partially damped during the final phase of gravitational radiation-driven inspiral. We estimate that the typical eccentricity at one week prior to coalescence is of the order of 0.01. Higher terminal eccentricities, which can approach 0.1, are possible if the binary has an extreme mass ratio. The detection of even a small eccentricity prior to merger by LISA provides a possible discriminant between gas-driven inspirals and those effected by stellar processes.
[96]  [pdf] - 73379
Quantifying orbital migration from exoplanet statistics and host metallicities
Comments: ApJ, in press
Submitted: 2005-05-30
We investigate how the statistical distribution of extrasolar planets may be combined with knowledge of the host stars' metallicity to yield constraints on the migration histories of gas giant planets. At any radius, planets that barely manage to form around the lowest metallicity stars accrete their envelopes just as the gas disk is being dissipated, so the lower envelope of planets in a plot of metallicity vs semi-major axis defines a sample of non-migratory planets that will have suffered less than average migration subsequent to gap opening. Under the assumption that metallicity largely controls the initial surface density of planetesimals, we use simplified core accretion models to calculate how the minimum metallicity needed for planet formation varies as a function of semi-major axis. Models that do not include core migration prior to gap opening (Type I migration) predict that the critical metallicity is largely flat between the snow line and a semimajor axis of about 6 AU, with a weak dependence on the initial surface density profile of planetesimals. When slow Type I migration is included, the critical metallicity is found to increase steadily from 1-10 AU. Large planet samples, that include planets at modestly greater orbital radii than present surveys, therefore have the potential to quantify the extent of migration in both Type I and Type II regimes.
[97]  [pdf] - 70468
Influence of massive planet scattering on nascent terrestrial planets
Comments: 6 pages, 2 figures, 2 tables, accepted for publication in ApJL
Submitted: 2005-01-17
In most extrasolar planetary systems, the present orbits of known giant planets admit the existence of stable terrestrial planets. Those same giant planets, however, have typically eccentric orbits that hint at violent early dynamics less benign for low mass planet formation. Under the assumption that massive planet eccentricities are the end point of gravitational scattering in multiple planet systems, we study the evolution of the building blocks of terrestrial planets during the scattering process. We find that typically, evolutionary sequences that result in a moderately eccentric giant planet orbiting at a ~ 2.5 AU eject over 95% of the material initially present within the habitable zone. Crossing orbits largely trigger the ejection, and leave the surviving material with a wide dispersion in semi-major axis, eccentricity and inclination. Based on these results, we predict that radial velocity follow-up of terrestrial planet systems found by Kepler will find that these are anti-correlated with the presence of eccentric giant planets orbiting at a few AU.
[98]  [pdf] - 66881
Accelerated planetesimal growth in self-gravitating protoplanetary discs
Comments: 11 pages, 9 figures, Accepted for publication in MNRAS
Submitted: 2004-08-21
In this paper we consider the evolution of small planetesimals in marginally stable, self-gravitating protoplanetary discs. The drag force between the disc gas and the embedded planetesimals generally causes the planetesimals to drift inwards through the disc at a rate that depends on the particle size. In a marginally stable, self-gravitating disc, however, the planetesimals are significantly influenced by the non-axisymmetric spiral structures resulting from the growth of the gravitational instability. For small particles, that are strongly coupled to the disc gas, and for large particles, that have essentially decoupled from the disc gas, the effect is not particularly significant. Intermediate sized particles, which would generally have the largest radial drift rates, do, however, become significantly concentrated at the peaks of the spiral arms. These high density regions may persist for, of order, an orbital period and may attain densities comparable to that of the disc gas. Although at the end of the simulation only ~ 25 % of the planetesimal particles lie in regions of enhanced density, during the course of the simulation at least 75 % of the planetesimal particles have at some stage been in a such a region. We find that the concentration of particles in the spiral arms results in an increased collision rate, an effect that could significantly accelerate planetesimal growth. The density enhancements may also be sufficient for the growth of planetesimals through direct gravitational collapse. The interaction between small planetesimals and self-gravitating spiral structures may therefore play an important role in the formation of large planetesimals that will ultimately coagulate to form terrestrial planets or the cores of gas/ice giant planets.
[99]  [pdf] - 65844
The Dynamics of Two Massive Planets on Inclined Orbits
Comments: 25 pages, 15 figures, 2 tables, accepted for publication in Icarus
Submitted: 2004-06-30
The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.
[100]  [pdf] - 1348135
Theory of disk accretion onto supermassive black holes
Comments: Mostly introductory review, to appear in "Supermassive black holes in the distant Universe", ed. A.J. Barger, Kluwer Academic Publishers, in press
Submitted: 2004-05-09
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.
[101]  [pdf] - 61532
MHD simulations of the collapsar model for GRBs
Comments: 4 pages, to appear in proceedings of "2003 GRB Conference" (Santa Fe, NM, September 8-12, 2003), needs aipprocs LaTeX class, movies are available at
Submitted: 2003-12-11
We present results from axisymmetric, time-dependent magnetohydrodynamic (MHD) simulations of the collapsar model for gamma-ray bursts. Our main conclusion is that, within the collapsar model, MHD effects alone are able to launch, accelerate and sustain a strong polar outflow. We also find that the outflow is Poynting flux-dominated, and note that this provides favorable initial conditions for the subsequent production of a baryon-poor fireball.
[102]  [pdf] - 60342
Substellar companions and isolated planetary mass objects from protostellar disc fragmentation
Comments: 5 pages, accepted for publication in MNRAS
Submitted: 2003-10-23
Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time scale that is short compared to the orbital period. We use a combination of hydrodynamic simulations and N-body orbit integrations to study the long term evolution of a fragmenting disc with an initial mass ratio to the star of M_disc/M_star = 0.1. For a disc which is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar mass star) up to ~0.01 M_sun. Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 -- in which the companions are close to or beyond the deuterium burning limit -- appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ~1 au should differ from that of stars with lower mass planetary companions.
[103]  [pdf] - 59854
On the formation time scale and core masses of gas giant planets
Comments: ApJL, in press
Submitted: 2003-10-07
Numerical simulations show that the migration of growing planetary cores may be dominated by turbulent fluctuations in the protoplanetary disk, rather than by any mean property of the flow. We quantify the impact of this stochastic core migration on the formation time scale and core mass of giant planets at the onset of runaway gas accretion. For standard Solar Nebula conditions, the formation of Jupiter can be accelerated by almost an order of magnitude if the growing core executes a random walk with an amplitude of a few tenths of an au. A modestly reduced surface density of planetesimals allows Jupiter to form within 10 Myr, with an initial core mass below 10 Earth masses, in better agreement with observational constraints. For extrasolar planetary systems, the results suggest that core accretion could form massive planets in disks with lower metallicities, and shorter lifetimes, than the Solar Nebula.
[104]  [pdf] - 59824
Migration of extrasolar planets to large orbital radii
Comments: MNRAS, in press
Submitted: 2003-10-06
Observations of structure in circumstellar debris discs provide circumstantial evidence for the presence of massive planets at large (several tens of au) orbital radii, where the timescale for planet formation via core accretion is prohibitively long. Here, we investigate whether a population of distant planets can be produced via outward migration subsequent to formation in the inner disc. Two possibilities for significant outward migration are identified. First, cores that form early at radii of around 10 au can be carried to larger radii via gravitational interaction with the gaseous disc. This process is efficient if there is strong mass loss from the disc - either within a cluster or due to photoevaporation from a star more massive than the Sun - but does not require the extremely destructive environment found, for example, in the core of the Orion Nebula. We find that, depending upon the disc model, gas disc migration can yield massive planets (several Jupiter masses) at radii of around 20-50 au. Second, interactions within multiple planet systems can drive the outer planet into a large, normally highly eccentric orbit. A series of scattering experiments suggests that this process is most efficient for lower mass planets within systems of unequal mass ratio. This mechanism is a good candidate for explaining the origin of relatively low mass giant planets in eccentric orbits at large radii.
[105]  [pdf] - 59665
Axisymmetric MHD simulations of the collapsar model for gamma-ray bursts
Comments: submitted in ApJ Letters
Submitted: 2003-09-30
We present results from axisymmetric, time-dependent magnetohydrodynamic (MHD) simulations of the collapsar model for gamma-ray bursts. We begin the simulations after the 1.7 MSUN iron core of a 25 MSUN presupernova star has collapsed and study the ensuing accretion of the 7 MSUN helium envelope onto the central black hole formed by the collapsed iron core. We consider a spherically symmetric progenitor model, but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum and a weak radial magnetic field. Our MHD simulations include a realistic equation of state, neutrino cooling, photodisintegration of helium, and resistive heating. Our main conclusion is that, within the collapsar model, MHD effects alone are able to launch, accelerate and sustain a strong polar outflow. We also find that the outflow is Poynting flux-dominated, and note that this provides favorable initial conditions for the subsequent production of a baryon-poor fireball.
[106]  [pdf] - 58486
On the stability of self-gravitating protoplanetary discs
Comments: 6 pages, 6 figures, to appear in Proceedings of 'DARWIN / TPF and the search for extrasolar terrestrial planets'
Submitted: 2003-08-12
It has already been shown, using a local model, that accretion discs with cooling times t_cool <= 3 Omega^-1 fragment into gravitationally bound objects, while those with cooling times t_cool > 3 Omega^-1 evolve into a quasi-steady state. We present results of three-dimensional simulations that test if the local result still holds globally. We find that the fragmentation boundary is close to that determined using the local model, but that fragmentation may occur for longer cooling times when the disc is more massive or when the mass is distributed in such a way as to make a particular region of the disc more susceptible to the growth of the gravitational instability. These results have significant implications for the formation of of gaseous planets in protoplanetary discs and also for the redistribution of angular momentum which could be driven by the presence of relatively massive, bound objects within the disc.
[107]  [pdf] - 57816
The spectroscopic signature of hot Jupiters in FU Orionis objects
Comments: MNRAS, in press
Submitted: 2003-07-04
We show that if FU Orionis objects harbour hot Jupiters embedded in their discs, the resulting non-axisymmetric dissipation profile in the disc would be manifest as time-dependent distortions in the absorption line profiles of these objects. In order to affect the infrared line profiles, planets must lie within approximately 0.5 au of the central star, whereas only planets within 0.1 au would influence the optical line profiles. The timescale for modulation of the line profiles is relatively short (months) in each case, so that the effect could not have been discovered from published spectra (which combine data taken in different observing seasons). The detection of hot Jupiters in FU Orionis objects would be in line with the expectations of tidal migration theories (which predict a high incidence of close planets around young stars) and would also lend support to models which link the triggering of rapid rise FU Orionis events to the existence of a close massive planet.
[108]  [pdf] - 55582
Dispersion in the lifetime and accretion rate of T Tauri discs
Comments: MNRAS, in press
Submitted: 2003-03-14
We compare evolutionary models for protoplanetary discs that include disc winds with observational determinations of the disc lifetime and accretion rate in Taurus. Using updated estimates for stellar ages in Taurus, together with published classifications, we show that the evolution of the disc fraction with stellar age is similar to that derived for ensembles of stars within young clusters. Around 30 percent of stars lose their discs within 1 Myr, while the remainder have disc lifetimes that are typically in the 1-10 Myr range. We show that the latter range of ages is consistent with theoretical models for disc evolution, provided that there is a dispersion of around 0.5 in the log of the initial disc mass. The same range of initial conditions brackets the observed variation in the accretion rate of Classical T Tauri stars at a given age. We discuss the expected lifetime of discs in close binary systems, and show that our models predict that the disc lifetime is almost constant for separations exceeding 10 au. This implies a low predicted fraction of binaries that pair a Classical T Tauri star with a Weak-lined T Tauri star, and is in better agreement with observations of the disc lifetime in binaries than disc models that do not include disc mass loss in a wind.
[109]  [pdf] - 54898
The variability of accretion onto Schwarzschild black holes from turbulent magnetized discs
Comments: MNRAS, in press. For accompanying movies see
Submitted: 2003-02-13
We use global magnetohydrodynamic simulations, in a pseudo-Newtonian potential, to investigate the temporal variability of accretion discs around Schwarzschild black holes. We use the vertically-averaged magnetic stress in the simulated disc as a proxy for the rest-frame dissipation, and compute the observed emission by folding this through the transfer function describing the relativistic beaming, light bending and time delays near a non-rotating black hole. The temporal power spectrum of the predicted emission from individual annuli in the disc is described by a broken power law, with indices of -3.5 at high frequency and 0 to -1 at low frequency. Integrated over the disc, the power spectrum is approximated by a single power law with an index of -2. Increasing inclination boosts the relative power at frequencies of around 0.3 times the orbital frequency at the marginally stable orbit, but no evidence is found for sharp quasi-periodic oscillations in the lightcurve. Assuming that fluorescent iron line emission locally tracks the continuum flux, we compute simulated broad iron line profiles. We find that relativistic beaming of the non-axisymmetric emission profile, induced by turbulence, produces high-amplitude variability in the iron line profile. We show that this substructure within the broad iron line profile can survive averaging over a number of orbital periods, and discuss the origin of the anomalous X-ray spectral features, recently reported by Turner et al. (2002) for the Seyfert galaxy NGC3516, in the context of turbulent disc models.
[110]  [pdf] - 144689
Constraints on a planetary origin for the gap in GM Aurigae's protoplanetary disc
Comments: 8 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2003-02-11
The unusual spectral energy distribution (SED) of the classical T Tauri star GM Aurigae provides evidence for the presence of an inner disc hole extending to several au. Using a combination of hydrodynamical simulations and Monte Carlo radiative transport, we investigate whether the observed SED is consistent with the inner hole being created and maintained by an orbiting planet. We show that an ~ 2 M_Jupiter planet, orbiting at 2.5 au in a disc with mass 0.047 M_sun and radius 300 au, provides a good match both to the SED and to CO observations which constrain the velocity field in the disc. A range of planet masses is allowed by current data, but could in principle be distinguished with further observations between 3 and ~ 20 microns. Future high precision astrometric instruments should also be able to detect the motion of the central star due to an orbiting Jupiter mass planet. We argue that the small number of T Tauri stars with SEDs resembling that of GM Aur is broadly consistent with the expected statistics of embedded migrating planets.
[111]  [pdf] - 53223
A reduced efficiency of terrestrial planet formation following giant planet migration
Comments: ApJL, in press
Submitted: 2002-11-21
Substantial orbital migration of massive planets may occur in most extrasolar planetary systems. Since migration is likely to occur after a significant fraction of the dust has been locked up into planetesimals, ubiquitous migration could reduce the probability of forming terrestrial planets at radii of the order of 1 au. Using a simple time dependent model for the evolution of gas and solids in the disk, I show that replenishment of solid material in the inner disk, following the inward passage of a giant planet, is generally inefficient. Unless the timescale for diffusion of dust is much shorter than the viscous timescale, or planetesimal formation is surprisingly slow, the surface density of planetesimals at 1 au will typically be depleted by one to two orders of magnitude following giant planet migration. Conceivably, terrestrial planets may exist only in a modest fraction of systems where a single generation of massive planets formed and did not migrate significantly.
[112]  [pdf] - 52823
The effect of cooling on the global stability of self-gravitating protoplanetary discs
Comments: 6 pages, 7 figures, accepted for publication in MNRAS
Submitted: 2002-11-05
Using a local model Gammie (2001) has shown that accretion discs with cooling times t_cool <= 3/Omega fragment into gravitationally bound objects, while those with cooling times t_cool > 3/Omega evolve into a quasi-steady state. We use three-dimensional smoothed particle hydrodynamic simulations of protoplanetary accretion discs to test if the local results hold globally. We find that for disc masses appropriate for T Tauri discs, the fragmentation boundary still occurs at a cooling time close to t_cool = 3/Omega. For more massive discs, which are likely to be present at an earlier stage of the star formation process, fragmentation occurs for longer cooling times, but still within a factor of two of that predicted using a local model. These results have implications not only for planet formation in protoplanetary discs and star formation in AGN discs, but also for the redistribution of angular momentum which could be driven by the presence of relatively massive objects within the accretion disc.
[113]  [pdf] - 52233
Orbital migration and the brown dwarf desert
Comments: Proceedings, IAU Symp. 211 on "Brown Dwarfs", ed. E. L. Martin
Submitted: 2002-10-08
The orbital elements of extreme mass ratio binaries will be modified by interactions with surrounding circumstellar disks. For brown dwarf companions to Solar-type stars the resulting orbital migration is sufficient to drive short period systems to merger, creating a brown dwarf desert at small separations. We highlight the similarities and the differences between the migration of brown dwarfs and massive extrasolar planets, and discuss how observations can test a migration model for the brown dwarf desert.
[114]  [pdf] - 51556
Astrometric signatures of self-gravitating protoplanetary discs
Comments: 6 pages, 9 figures, accepted for publication in MNRAS
Submitted: 2002-09-09
We use high resolution numerical simulations to study whether gravitational instabilities within circumstellar discs can produce astrometrically detectable motion of the central star. For discs with masses of M_disc = 0.1 M_star, which are permanantly stable against fragmentation, we find that the magnitude of the astrometric signal depends upon the efficiency of disc cooling. Short cooling times produce prominent filamentary spiral structures in the disc, and lead to stellar motions that are potentially observable with future high precision astrometric experiments. For a disc that is marginally unstable within radii of \~10 au, we estimate astrometric displacements of 10-100 microarcsec on decade timescales for a star at a distance of 100 pc. The predicted displacement is suppressed by a factor of several in more stable discs in which the cooling time exceeds the local dynamical time by an order of magnitude. We find that the largest contribution comes from material in the outer regions of the disc and hence, in the most pessimistic scenario, the stellar motions caused by the disc could confuse astrometric searches for low mass planets orbiting at large radii. They are, however, unlikely to present any complications in searches for embedded planets orbiting at small radii, relative to the disc size, or Jupiter mass planets or greater orbiting at large radii.
[115]  [pdf] - 48546
Predictions for the frequency and orbital radii of massive extrasolar planets
Comments: MNRAS, in press
Submitted: 2002-03-29
We investigate the migration of massive extrasolar planets due to gravitational interaction with a viscous protoplanetary disc. We show that a model in which planets form at 5 AU at a constant rate, before migrating, leads to a predicted distribution of planets that is a steeply rising function of log (a), where a is the orbital radius. Between 1 AU and 3 AU, the expected number of planets per logarithmic interval in orbital radius roughly doubles. We demonstrate that, once selection effects are accounted for, this is consistent with current data, and then extrapolate the observed planet fraction to masses and radii that are inaccessible to current observations. In total, about 15 percent of stars targeted by existing radial velocity searches are predicted to possess planets with masses 0.3 M_Jupiter < M_p sin (i) < 10 M_Jupiter, and radii 0.1 AU < a < 5 AU. A third of these planets (around 5 percent of the target stars) lie at the radii most amenable to detection via microlensing. A further 5-10 percent of stars could have planets at radii of 5 AU < a < 8 AU that have migrated outwards. We discuss the probability of forming a system (akin to the Solar System) in which significant radial migration of the most massive planet does not occur. About 10-15 percent of systems with a surviving massive planet are estimated to fall into this class. Finally, we note that a smaller fraction of low mass planets than high mass planets is expected to survive without being consumed by the star. The initial mass function for planets is thus predicted to rise more steeply towards small masses than the observed mass function.
[116]  [pdf] - 354071
Accretion during the merger of supermassive black holes
Comments: ApJL, in press
Submitted: 2002-01-18
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.
[117]  [pdf] - 46374
The brown dwarf desert as a consequence of orbital migration
Comments: MNRAS (Letters), in press
Submitted: 2001-11-30
We show that the dearth of brown dwarfs in short-period orbits around Solar-mass stars - the brown dwarf desert - can be understood as a consequence of inward migration within an evolving protoplanetary disc. Brown dwarf secondaries forming at the same time as the primary star have masses which are comparable to the initial mass of the protoplanetary disc. Subsequent disc evolution leads to inward migration, and destruction of the brown dwarf, via merger with the star. This is in contrast with massive planets, which avoid this fate by forming at a later epoch when the disc is close to being dispersed. Within this model, a brown dwarf desert arises because the mass at the hydrogen burning limit is coincidentally comparable to the initial disc mass for a Solar mass star. Brown dwarfs should be found in close binaries around very low mass stars, around other brown dwarfs, and around Solar-type stars during the earliest phases of star formation.
[118]  [pdf] - 45727
Magnetic activity in accretion disc boundary layers
Comments: MNRAS, in press. Movies from the simulations are available at
Submitted: 2001-10-30
We use three dimensional magnetohydrodynamic simulations to study the structure of the boundary layer between an accretion disc and a non-rotating, unmagnetized star. Under the assumption that cooling is efficient, we obtain a narrow but highly variable transition region in which the radial velocity is only a small fraction of the sound speed. A large fraction of the energy dissipation occurs in high density gas adjacent to the hydrostatic stellar envelope, and may therefore be reprocessed and largely hidden from view of the observer. As suggested by Pringle (1989), the magnetic field energy in the boundary layer is strongly amplified by shear, and exceeds that in the disc by an order of magnitude. These fields may play a role in generating the magnetic activity, X-ray emission, and outflows in disc systems where the accretion rate is high enough to overwhelm the stellar magnetosphere.
[119]  [pdf] - 45307
Ab initio simulations of accretion disc boundary layers
Comments: To appear in "The Physics of Cataclysmic Variables and Related Objects", eds. B. Gaensicke, K. Beuermann & K. Reinsch. Animations from boundary layer and black hole plunging region simulations at
Submitted: 2001-10-10
I discuss the results of simplified three dimensional magnetohydrodynamic simulations of the boundary layer between a disc and a non-rotating, unmagnetized star. Strong magnetic fields, possibly approaching equipartition with the thermal energy, occur in the boundary layer due to the shearing of disc-generated fields. The mean boundary layer magnetic field, which is highly variable on an orbital timescale, is estimated to exceed 50 kG for a CV with an accretion rate of 10^(-9) solar masses per year. However, these fields do not drive efficient angular momentum transport within the boundary layer. As a consequence the radial velocity in the boundary layer is low, and the density high.
[120]  [pdf] - 45086
A variable efficiency for thin disk black hole accretion
Comments: 4 pages, 4 postscript figures. Accepted for publication in Astrophysical Journal Letters
Submitted: 2001-10-01
We explore the presence of torques at the inner edges of geometrically-thin black hole accretion disks using 3-dimensional magnetohydrodynamic (MHD) simulations in a pseudo-Newtonian potential. By varying the saturation level of the magnetorotational instability that leads to angular momentum transport, we show that the dynamics of gas inside the radius of marginal stability varies depending upon the magnetic field strength just outside that radius. Weak fields are unable to causally connect material within the plunging region to the rest of the disk, and zero torque is an approximately correct boundary condition at the radius of marginal stability. Stronger fields, which we obtain artificially but which may occur physically within more complete disk models, are able to couple at least some parts of the plunging region to the rest of the disk. In this case, angular momentum (and implicitly energy) is extracted from the material in the plunging region. Furthermore, the magnetic coupling to the plunging region can be highly time dependent with large fluctuations in the torque at the radius of marginal stability. This implies varying accretion efficiencies, both across systems and within a given system at different times. The results suggest a possible link between changes in X-ray and outflow activity, with both being driven by transitions between weak and strong field states.
[121]  [pdf] - 42139
Eccentric extrasolar planets from resonant Saturn mass companions
Comments: Withdrawn. Mechanism works convincingly only for the higher mass companions observed in the known resonant systems
Submitted: 2001-04-24, last modified: 2001-07-18
Withdrawn in current form.
[122]  [pdf] - 40756
Magnetohydrodynamic simulations of black hole accretion
Comments: 6 pages, to appear in proceedings, 20th Texas Symposium on Relativistic Astrophysics, eds J.C. Wheeler and H. Martel
Submitted: 2001-02-02
We discuss the results of three-dimensional magnetohydrodynamic simulations, using a pseudo-Newtonian potential, of thin disk (h/r ~ 0.1) accretion onto black holes. We find (i) that magnetic stresses persist within the marginally stable orbit, and (ii) that the importance of those stresses for the dynamics of the flow depends upon the strength of magnetic fields in the disk outside the last stable orbit. Strong disk magnetic fields (alpha > 0.1) lead to a gross violation of the zero-torque boundary condition at the last stable orbit, while weaker fields (alpha ~ 0.01) produce results more akin to traditional models for thin disk accretion onto black holes. Fluctuations in the magnetic field strength in the disk could lead to changes in the radiative efficiency of the flow on short timescales.
[123]  [pdf] - 40391
Episodic accretion in magnetically layered protoplanetary discs
Comments: MNRAS, in press
Submitted: 2001-01-15
We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 AU such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion onto the protostar occurs in bursts with an accretion rate of 10^{-5} solar masses / yr, separated by quiescent intervals where the accretion rate is 10^{-8} solar masses / yr. Such bursts could lead to repeated episodes of strong mass outflow in Young Stellar Objects. The transition to this episodic mode of accretion occurs at an early epoch (t < 1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at around 1 AU. The predicted rate of low mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1-3 AU.
[124]  [pdf] - 1468133
Simulations of accretion flows crossing the last stable orbit
Comments: ApJ, in press. Revised to match accepted version. References corrected
Submitted: 2000-07-04, last modified: 2000-10-11
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.
[125]  [pdf] - 36882
Suppression of giant planet formation in stellar clusters
Comments: A&A in press. Minor revisions to match accepted version
Submitted: 2000-07-04, last modified: 2000-09-02
Photoevaporation driven by the ultraviolet radiation from massive stars severely limits the lifetime of protoplanetary discs around stars formed within stellar clusters. I investigate the resulting influence of clustered environments on the probability of giant planet formation, and show that for clusters as rich, or richer than, Orion, the time available for planet formation is likely to be limited to the length of any delay between low mass and high mass star formation. Under popular models for the formation of massive planets, the fraction of stars with giant planets in rich clusters is expected to be substantially suppressed as compared to less clustered star formation environments.
[126]  [pdf] - 109810
Early planet formation as a trigger for further planet formation
Comments: Nature in press, short
Submitted: 1999-12-07
Recent discoveries of extrasolar planets at small orbital radii, or with significant eccentricities, indicate that interactions between massive planets and the disks of gas and dust from which they formed are vital for determining the final shape of planetary systems. We show that if this interaction occurs at an early epoch, when the protoplanetary disc was still massive, then rapid planet growth through accretion causes an otherwise stable disc to fragment into additional planetary mass bodies when the planetary mass reaches 4-5 Jupiter masses. We suggest that such catastrophic planet formation could account for apparent differences in the mass function of massive planets and brown dwarfs, and the existence of young stars that appear to have dissipated their discs at an early epoch. Subsequent gravitational interactions will lead to planetary systems comprising a small number of massive planets in eccentric orbits.
[127]  [pdf] - 107523
The Blandford-Znajek mechanism and emission from isolated accreting black holes
Comments: ApJL, in press
Submitted: 1999-07-21
In the presence of a magnetic field, rotational energy can be extracted from black holes via the Blandford-Znajek mechanism. We use self-similar advection dominated accretion (ADAF) models to estimate the efficiency of this mechanism for black holes accreting from geometrically thick disks, in the light of recent magnetohydrodynamic disk simulations, and show that the power from electromagnetic energy extraction exceeds the accretion luminosity for ADAFs at sufficiently low accretion rates. We consider the detectability of isolated stellar mass black holes accreting from the ISM, and show that for any rapidly rotating holes the efficiency of energy extraction could reach 0.01. The estimated total luminosity would be consistent with the tentative identification of some EGRET sources as accreting isolated black holes, if that energy is radiated primarily as gamma rays. We discuss the importance of emission from the Blandford-Znajek mechanism for the spectra of other advection dominated accretion flows, especially those in low luminosity galactic nuclei.
[128]  [pdf] - 104108
Warped discs and the directional stability of jets in Active Galactic Nuclei
Comments: MNRAS, in press. Revised to match accepted version, with one new figure showing alignment timescale as a function of black hole mass
Submitted: 1998-11-30, last modified: 1999-06-26
Warped accretion discs in Active Galactic Nuclei (AGN) exert a torque on the black hole that tends to align the rotation axis with the angular momentum of the outer disc. We compute the magnitude of this torque by solving numerically for the steady state shape of the warped disc, and verify that the analytic solution of Scheuer and Feiler (1996) provides an excellent approximation. We generalise these results for discs with strong warps and arbitrary surface density profiles, and calculate the timescale on which the black hole becomes aligned with the angular momentum in the outer disc. For massive holes and accretion rates of the order of the Eddington limit the alignment timescale is always short (less than a Myr), so that jets accelerated from the inner disc region provide a prompt tracer of the angular momentum of gas at large radii in the disc. Longer timescales are predicted for low luminosity systems, depending on the degree of anisotropy in the disc's hydrodynamic response to shear and warp, and for the final decay of modest warps at large radii in the disc that are potentially observable via VLBI. We discuss the implications of this for the inferred accretion history of those Active Galactic Nuclei whose jet directions appear to be stable over long timescales. The large energy deposition rate at modest disc radii during rapid realignment episodes should make such objects transiently bright at optical and infrared wavelengths.
[129]  [pdf] - 105293
Lense-Thirring precession of accretion disks around compact objects
Comments: ApJ, in press. Minor revisions to match accepted version. Animations showing warp evolution are available at
Submitted: 1999-02-18, last modified: 1999-06-26
Misaligned accretion disks surrounding rotating compact objects experience a torque due to the Lense-Thirring effect, which leads to precession of the inner disk. It has been suggested that this effect could be responsible for some low frequency Quasi-Periodic Oscillations observed in the X-ray lightcurves of neutron star and galactic black hole systems. We investigate this possibility via time-dependent calculations of the response of the inner disk to impulsive perturbations for both Newtonian point mass and Paczynski-Wiita potentials, and compare the results to the predictions of the linearized twisted accretion disk equations. For most of a wide range of disk models that we have considered, the combination of differential precession and viscosity causes the warps to decay extremely rapidly. Moreover, at least for relatively slowly rotating objects, linear calculations in a Newtonian point mass potential provide a good measure of the damping rate, provided only that the timescale for precession is much shorter than the viscous time in the inner disk. The typically rapid decay rates suggest that coherent precession of a fluid disk would not be observable, though it remains possible that the damping rate of warp in the disk could be low enough to permit weakly coherent signals from Lense-Thirring precession.
[130]  [pdf] - 106749
Black hole formation via hypercritical accretion during common envelope evolution
Comments: ApJ, submitted
Submitted: 1999-06-01
Neutron stars inspiralling into a stellar envelope can accrete at rates vastly exceeding the Eddington limit if the flow develops pressures high enough to allow neutrinos to radiate the released gravitational energy. It has been suggested that this hypercritical mode of accretion leads inevitably to the formation of stellar mass black holes during common envelope evolution. We study the hydrodynamics of this flow at large radii (R >> R_ns), and show that for low Mach number flows, in two dimensions, modest density gradients in the stellar envelope suffice to produce a hot, advection dominated accretion disk around the accreting object. The formation of outflows from such a disk is highly probable, and we discuss the impact of the resultant mass loss and feedback of energy into the envelope for the survival of the neutron star. Unless outflows are weaker than those inferred for well observed accreting systems, we argue that in most cases insufficient accretion occurs to force collapse to a black hole before the envelope has been ejected. This conclusions is of interest for black hole formation in general, for some models of gamma ray bursts, and for predictions of the event rate in future LIGO observations.
[131]  [pdf] - 103606
Accretion Disc Evolution in Single and Binary T Tauri Stars
Comments: 11 pages, MNRAS in press
Submitted: 1998-10-29
We present theoretical models for the evolution of T Tauri stars surrounded by circumstellar discs. The models include the effects of pre-main-sequence stellar and time dependent disc evolution, and incorporate the effects of stellar magnetic fields acting on the inner disc. For single stars, consistency with observations in Taurus-Auriga demands that disc dispersal occurs rapidly, on much less than the viscous timescale of the disc, at roughly the epoch when heating by stellar radiation first dominates over internal viscous dissipation. Applying the models to close binaries, we find that because the initial conditions for discs in binaries are uncertain, studies of extreme mass ratio systems are required to provide a stringent test of theoretical disc evolution models. We also note that no correlation of the infra-red colours of T Tauri stars with their rotation rate is observed, in apparent contradiction to the predictions of simple magnetospheric accretion models.
[132]  [pdf] - 102829
Tilted accretion discs in cataclysmic variables: tidal instabilities and superhumps
Comments: 7 pages, 7 figures, accepted for publication in MNRAS
Submitted: 1998-09-10
We investigate the growth of tidal instabilities in accretion discs in a binary star potential, using three dimensional numerical simulations. As expected from analytic work, the disc is prone to an eccentric instability provided that it is large enough to extend to the 3:1 resonance. The eccentric disc leads to positive superhumps in the light curve. It has been proposed that negative superhumps might arise from a tilted disc, but we find no evidence that the companion gravitational tilt instability can grow fast enough in a fluid disc to create a measurable inclination. The origin of negative superhumps in the light curves of cataclysmic variables remains a puzzle.
[133]  [pdf] - 102828
The effects of tidally induced disc structure on white dwarf accretion in intermediate polars
Comments: 9 pages, 5 figures, accepted for publication in MNRAS
Submitted: 1998-09-10
We investigate the effects of tidally induced asymmetric disc structure on accretion onto the white dwarf in intermediate polars. Using numerical simulation, we show that it is possible for tidally induced spiral waves to propagate sufficiently far into the disc of an intermediate polar that accretion onto the central white dwarf could be modulated as a result. We suggest that accretion from the resulting asymmetric inner disc may contribute to the observed X-ray and optical periodicities in the light curves of these systems. In contrast to the stream-fed accretion model for these periodicities, the tidal picture predicts that modulation can exist even for systems with weaker magnetic fields where the magnetospheric radius is smaller than the radius of periastron of the mass transfer stream. We also predict that additional periodic components should exist in the emission from low mass ratio intermediate polars displaying superhumps.
[134]  [pdf] - 101354
Turbulence and angular momentum transport in a global accretion disk simulation
Comments: ApJ Letters, in press. The paper and additional visualizations are available at
Submitted: 1998-05-11
The global development of magnetohydrodynamic turbulence in an accretion disk is studied within a simplified disk model that omits vertical stratification. Starting with a weak vertical seed field, a saturated state is obtained after a few tens of orbits in which the energy in the predominantly toroidal magnetic field is still subthermal. The efficiency of angular momentum transport, parameterized by the Shakura-Sunyaev alpha parameter, is of the order of 0.1. The dominant contribution to alpha comes from magnetic stresses, which are enhanced by the presence of weak net vertical fields. The power spectra of the magnetic fields are flat or decline only slowly towards the largest scales accessible in the calculation, suggesting that the viscosity arising from MHD turbulence may not be a locally determined quantity. I discuss how these results compare with observationally inferred values of alpha, and possible implications for models of jet formation.
[135]  [pdf] - 101189
Simulations of spiral structure in the accretion disc of IP Pegasi during outburst
Comments: 6 pages, MNRAS in press. Color figures at
Submitted: 1998-04-28
We consider the implications of the detection of spiral structure in the accretion disc of the binary IP Pegasi. We use numerical simulations of the development of a disc outburst to construct predicted Doppler tomograms, which are found to be in close agreement with the observations if the spiral pattern arises as a transient feature when the disc expands viscously at the start of the outburst. The good agreement of such viscous disc simulations with the data is consistent with models in which most of the angular momentum transport in the disc originates in internal stresses rather than globally excited waves or shocks. Future detailed observations of the development of transient spiral features offer the potential to measure the dependence of the disc viscosity on the local physical conditions in the disc.
[136]  [pdf] - 1469557
Hydrodynamics of the stream-disk impact in interacting binaries
Comments: 16 pages, including 8 figures. 1 color figure as a jpeg. ApJ, in press
Submitted: 1997-09-03
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.
[137]  [pdf] - 98325
Radiation induced warping of protostellar accretion disks
Comments: 12 pages, including 3 figures. ApJ Letters, in press
Submitted: 1997-08-15
We examine the consequences of radiatively driven warping of accretion disks surrounding pre-main-sequence stars. These disks are stable against warping if the luminosity arises from a steady accretion flow, but are unstable at late times when the intrinsic luminosity of the star overwhelms that provided by the disk. Warps can be excited for stars with luminosities of around 10 solar luminosities or greater, with larger and more severe warps in the more luminous systems. A twisted inner disk may lead to high extinction towards stars often viewed through their disks. After the disk at all radii becomes optically thin, the warp decays gradually on the local viscous timescale, which is likely to be long. We suggest that radiation induced warping may account for the origin of the warped dust disk seen in Beta Pictoris, if the star is only around 10-20 Myr old, and could lead to non-coplanar planetary systems around higher mass stars.
[138]  [pdf] - 98178
The stream-disk interaction
Comments: 8 pages, including 3 figures. To appear in 13th North American workshop on CVs. Paper and color visualizations also available at
Submitted: 1997-07-31
I review theoretical aspects of the interaction between the accretion stream and the disk in interacting binary systems, concentrating on recent hydrodynamic calculations. At low accretion rates, cooling is expected to be efficient, and the interaction leads to a nearly ballistic stream overflowing the disk rim towards smaller radii. If cooling is ineffectual, the shocked gas produces a bulge on the disk rim, and there is no coherent stream inward of the disk edge. Results are presented for the mass fraction and velocity structure of the overflowing component, and the implications for X-ray observations of `dips', and doppler tomography of cataclysmic variables are briefly discussed.
[139]  [pdf] - 1234807
The ejection of T Tauri stars from molecular clouds and the fate of circumstellar discs
Comments: 8 pages, 3 figures, LaTeX using mn.sty. Accepted by MNRAS. Also available in HTML format at
Submitted: 1996-10-21
We investigate the evolution of circumstellar discs around T Tauri stars that are ejected from small stellar clusters within molecular clouds. In particular, we study how the interaction that leads to ejection may hasten the transition between Classical and Weak-lined T Tauri status. In our models, ejections of T Tauri stars at velocities of 3-10 km/s truncate the accretion disc at radii between 1 and 10 a.u., reducing the viscous evolution time of the disc so that accretion rapidly ceases. The observational appearance of the resulting systems is then dependent on the presence or absence of a stellar magnetic field. For non-magnetic stars we find that a near-infra red excess should persist due to reprocessing of stellar radiation, but that this is greatly diminished for magnetic T Tauri stars by the presence of a magnetosphere extending to corotation. In either case, there is a period when ejected stars should appear as non-accreting systems with detectable circumstellar material at wavelengths of 5 microns and beyond. We discuss the implications of these results for models in which ejected stars contribute to the halo of pre-main-sequence objects discovered from ROSAT observations of star forming regions and the All-Sky Survey.
[140]  [pdf] - 94798
Accretion disks in interacting binaries: Simulations of the stream-disk impact
Comments: LaTeX (aaspp.sty), 19 pages including 7 of 8 figures. The paper and remaining figure are also available at . ApJ in press
Submitted: 1996-06-06
We investigate the impact between the gas stream from the inner Lagrangian point and the accretion disk in interacting binaries, using three dimensional Smooth Particle Hydrodynamics simulations. We find that a significant fraction of the stream material can ricochet off the disk edge and overflow towards smaller radii, and that this generates pronounced non-axisymmetric structure in the absorption column towards the central object. We discuss the implications of our results for observations and time-dependent models of low-mass X-ray binaries, cataclysmic variables and supersoft X-ray sources.
[141]  [pdf] - 94693
Red giant-disk encounters: Food for quasars?
Comments: 24 pages of LaTeX, including 10 of 11 figures. Uses aaspp.sty. Colour versions of Fig. 4 (not included) are available at . ApJ in press
Submitted: 1996-05-22
We explore the role that red giants might play in the central regions of Active Galactic Nuclei. Due to their large radii and the low binding energy of the stellar envelope, giants are vulnerable to envelope stripping from collisions with the accretion disk. Using hydrodynamic simulations we show that such collisions will typically deposit a substantial fraction of the envelope mass into the disk on each passage. Repeated encounters will then lead to the complete destruction of the star save for the dense core. We estimate the rate of fuel supply by this mechanism using simple models for the AGN disk and central stellar cluster. If the central stellar density is of order 10^7 solar masses per cubic pc, then stripping of giants could account for the activity of typical AGN provided that the accretion disk extends beyond 0.1 pc. For AGN with smaller disks, or clusters of lower central density, giant stripping could be an important source of gas enriched via stellar nucleosynthesis.
[142]  [pdf] - 1234526
Magnetic braking of T Tauri stars
Comments: 10 pages of uu-encoded gz-compressed Postscript, including 9 figures. MNRAS in press
Submitted: 1995-12-05
We construct models for the rotation rates of T Tauri stars whose spin is regulated by magnetic linkage between the star and a surrounding accretion disc. Our models utilise a time-dependent disc code to follow the accretion process and include the effects of pre-main-sequence stellar evolution. We find that the initial disc mass controls the evolution of the star-disc system. For sufficiently massive discs, a stellar field of $\sim$ 1 kG is able to regulate the spin rate to the observed values during the classical T Tauri phase. The field then acts to expel the disc and the star spins up at constant angular momentum as a weak-line system. Lower mass discs are ejected at an early epoch and fail to brake the star significantly. We extend the model to close binary systems, and find that the removal of angular momentum from the disc by the secondary significantly prolongs the inner disc lifetime. Such systems should therefore be relatively slow rotators. We also discuss the implications of our model for the spectral energy distributions and variability of T Tauri stars.
[143]  [pdf] - 1348172
Dynamo driven accretion discs and dwarf nova eruptions
Comments: uu-encoded gz-compressed Postscript file, 18 pages including 6 figures. ApJ in press
Submitted: 1995-11-20
We explore the consequences of a magnetic dynamo origin for the viscosity in accretion discs, for the structure and evolution of discs in dwarf nova systems. We propose that the rapid cooling that sets in at the end of a dwarf nova eruption acts to inhibit the Balbus-Hawley instability, and thereby to quench dynamo action and so reduce disc viscosity. We demonstrate that a modified disc instability model can reproduce the basic properties of dwarf nova eruptions, as well as some properties of quiescent discs. We also discuss some observational consequences of our model.
[144]  [pdf] - 1234336
Magnetic cycles and photometric variability of T Tauri stars
Comments: uuencoded compressed postscript file. 8 pages including 3 figures. To appear in MNRAS.
Submitted: 1995-03-13
We consider the effects of stellar magnetic cycles on disc accretion in T Tauri stars where the inner region of the disc is disrupted by the stellar field. Using time-dependent disc models, we show that order unity variations in field strength can strongly modulate the accretion flow for cycle time-scales of years to decades. Photometric variability is expected to be greatest in the ultraviolet, with smaller but still significant variations probable in the near-infrared. We suggest that variable stellar magnetic fields may contribute to the long-time-scale photometric variability of T Tauri stars. The long-term variability of T Tauri stars may therefore provide a clue to the nature of magnetic dynamo cycles in these stars.
[145]  [pdf] - 1234320
Magnetically modulated accretion in T Tauri stars
Comments: LaTeX file (requires mn.sty), 4 pages, no figures or tables. To appear in MNRAS.
Submitted: 1995-02-11
We examine how accretion on to T Tauri stars may be modulated by a time-dependent `magnetic gate' where the inner edge of the accretion disc is disrupted by a varying stellar field. We show that magnetic field variations on time-scales shorter than 10^5 yr can modulate the accretion flow, thus providing a possible mechanism both for the marked photometric variability of T Tauri stars and for the possible conversion of T Tauri stars between classical and weak line status. We thus suggest that archival data relating to the spectrophotometric variability of T Tauri stars may provide an indirect record of magnetic activity cycles in low-mass pre-main-sequence stars.