Showing 1–50 of 166 results for author: Armitage, P J

Decretion disc size in Be/X-ray binaries depends upon the disc aspect ratio

Authors: Rebecca G. Martin, Stephen H. Lubow, Philip J. Armitage, Daniel J. Price

Abstract: With three-dimensional hydrodynamical simulations we show that the size of the decretion disc and the structure of the accretion flow onto the neutron star in a Be/X-ray binary strongly depends upon the disc aspect ratio, $H/R$. We simulate a Be star disc that is coplanar to the orbit of a circularly or moderately eccentric neutron star companion, thereby maximising the effects of tidal truncation… ▽ More With three-dimensional hydrodynamical simulations we show that the size of the decretion disc and the structure of the accretion flow onto the neutron star in a Be/X-ray binary strongly depends upon the disc aspect ratio, $H/R$. We simulate a Be star disc that is coplanar to the orbit of a circularly or moderately eccentric neutron star companion, thereby maximising the effects of tidal truncation. For low disc aspect ratio, $H/R\lesssim 0.1$, the disc is efficiently tidally truncated by the neutron star. Most material that escapes the Roche lobe of the Be star is accreted by the neutron star through tidal streams. For larger disc aspect ratio, the outflow rate through the Be star disc is higher, tidal truncation becomes inefficient, the disc fills the Roche lobe and extends to the orbit of the companion. Some material escapes the binary as a gas stream that begins near the L2 point. While the accretion rate onto the neutron star is higher, the fraction of the outflow that is accreted by the neutron star is smaller. Low density Be star discs are expected to be approximately isothermal, such that $H/R$ increases with radius. Tidal truncation is therefore weaker for larger separation binaries, and lower mass primaries. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS

A thermodynamic criterion for the formation of Circumplanetary Disks

Authors: Leonardo Krapp, Kaitlin M. Kratter, Andrew N. Youdin, Pablo Benítez-Llambay, Frédéric Masset, Philip J. Armitage

Abstract: The formation of circumplanetary disks is central to our understanding of giant planet formation, influencing their growth rate during the post-runaway phase and observability while embedded in protoplanetary disks. We use 3D global multifluid radiation hydrodynamics simulations with the FARGO3D code to define the thermodynamic conditions that enable circumplanetary disk formation around Jovian pl… ▽ More The formation of circumplanetary disks is central to our understanding of giant planet formation, influencing their growth rate during the post-runaway phase and observability while embedded in protoplanetary disks. We use 3D global multifluid radiation hydrodynamics simulations with the FARGO3D code to define the thermodynamic conditions that enable circumplanetary disk formation around Jovian planets on wide orbits. Our simulations include stellar irradiation, viscous heating, static mesh refinement, and active calculation of opacity based on evolving dust fluids. We find a necessary condition for the formation of circumplanetary disks in terms of a mean cooling time: when the cooling time is at least one order of magnitude shorter than the orbital time scale, the specific angular momentum of the gas is nearly Keplerian at scales of $R_{\rm{Hill}}/3$. We show that the inclusion of multifluid dust dynamics favors rotational support because dust settling produces an anisotropic opacity distribution that favors rapid cooling. In all our models with radiation hydrodynamics, specific angular momentum decreases as time evolves in agreement with the formation of an inner isentropic envelope due to compressional heating. The isentropic envelope can extend up to $R_{\rm{Hill}}/3$ and shows negligible rotational support. Thus, our results imply that young gas giant planets may host spherical isentropic envelopes, rather than circumplanetary disks. △ Less

Submitted 26 February, 2024; v1 submitted 22 February, 2024; originally announced February 2024.

Comments: 17 pages, 10 figures, submitted to ApJ. Comments are welcome

Streaming Instability and Turbulence: Conditions for Planetesimal Formation

Authors: Jeonghoon Lim, Jacob B. Simon, Rixin Li, Philip J. Armitage, Daniel Carrera, Wladimir Lyra, David G. Rea, Chao-Chin Yang, Andrew N. Youdin

Abstract: The streaming instability (SI) is a leading candidate for planetesimal formation, which can concentrate solids through two-way aerodynamic interactions with the gas. The resulting concentrations can become sufficiently dense to collapse under particle self-gravity, forming planetesimals. Previous studies have carried out large parameter surveys to establish the critical particle to gas surface den… ▽ More The streaming instability (SI) is a leading candidate for planetesimal formation, which can concentrate solids through two-way aerodynamic interactions with the gas. The resulting concentrations can become sufficiently dense to collapse under particle self-gravity, forming planetesimals. Previous studies have carried out large parameter surveys to establish the critical particle to gas surface density ratio ($Z$), above which SI-induced concentration triggers planetesimal formation. The threshold $Z$ depends on the dimensionless stopping time ($τ_s$, a proxy for dust size). However, these studies neglected both particle self-gravity and external turbulence. Here, we perform 3D stratified shearing box simulations with both particle self-gravity and turbulent forcing, which we characterize via $α_D$ that measures turbulent diffusion. We find that forced turbulence, at amplitudes plausibly present in some protoplanetary disks, can increase the threshold $Z$ by up to an order of magnitude. For example, for $τ_s = 0.01$, planetesimal formation occurs when $Z \gtrsim 0.06$, $\gtrsim 0.1$, and $\gtrsim 0.2$ at $α_D = 10^{-4}$, $10^{-3.5}$, and $10^{-3}$, respectively. We provide a single fit to the critical $Z$ as a function of $α_D$ and $τ_s$ required for the SI to work (though limited to the range $τ_s = 0.01$--0.1). Our simulations also show that planetesimal formation requires a mid-plane particle-to-gas density ratio that exceeds unity, with the critical value being independent of $α_D$. Finally, we provide the estimation of particle scale height that accounts for both particle feedback and external turbulence. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: 27 pages, 13 figures, submitted to The Astrophysical Journal

Tidal truncation of circumplanetary disks fails above a critical disk aspect ratio

Authors: Rebecca G. Martin, Philip J. Armitage, Stephen H. Lubow, Daniel J. Price

Abstract: We use numerical simulations of circumplanetary disks to determine the boundary between disks that are radially truncated by the tidal potential, and those where gas escapes the Hill sphere. We consider a model problem, in which a coplanar circumplanetary disk is resupplied with gas at an injection radius smaller than the Hill radius. We evolve the disk using the PHANTOM Smoothed Particle Hydrodyn… ▽ More We use numerical simulations of circumplanetary disks to determine the boundary between disks that are radially truncated by the tidal potential, and those where gas escapes the Hill sphere. We consider a model problem, in which a coplanar circumplanetary disk is resupplied with gas at an injection radius smaller than the Hill radius. We evolve the disk using the PHANTOM Smoothed Particle Hydrodynamics code until a steady-state is reached. We find that the most significant dependence of the truncation boundary is on the disk aspect ratio $H/R$. Circumplanetary disks are efficiently truncated for $H/R \lesssim 0.2$. For $H/R \simeq 0.3$, up to about half of the injected mass, depending on the injection radius, flows outwards through the decretion disk and escapes. As expected from analytic arguments, the conditions ($H/R$ and Shakura-Sunyaev $α$) required for tidal truncation are independent of planet mass. A simulation with larger $α=0.1$ shows stronger outflow than one with $α=0.01$, but the dependence on transport efficiency is less important than variations of $H/R$. Our results suggest two distinct classes of circumplanetary disks: tidally truncated thin disks with dust-poor outer regions, and thicker actively decreting disks with enhanced dust-to-gas ratios. Applying our results to the PDS 70c system, we predict a largely truncated circumplanetary disk, but it is possible that enough mass escapes to support an outward flow of dust that could explain the observed disk size. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

The distribution of accretion rates as a diagnostic of protoplanetary disc evolution

Authors: R. Alexander, G. Rosotti, P. J. Armitage, G. J. Herczeg, C. F. Manara, B. Tabone

Abstract: We show that the distribution of observed accretion rates is a powerful diagnostic of protoplanetary disc physics. Accretion due to turbulent ("viscous") transport of angular momentum results in a fundamentally different distribution of accretion rates than accretion driven by magnetised disc winds. We find that a homogeneous sample of $\gtrsim$300 observed accretion rates would be sufficient to d… ▽ More We show that the distribution of observed accretion rates is a powerful diagnostic of protoplanetary disc physics. Accretion due to turbulent ("viscous") transport of angular momentum results in a fundamentally different distribution of accretion rates than accretion driven by magnetised disc winds. We find that a homogeneous sample of $\gtrsim$300 observed accretion rates would be sufficient to distinguish between these two mechanisms of disc accretion at high confidence, even for pessimistic assumptions. Current samples of T Tauri star accretion rates are not this large, and also suffer from significant inhomogeneity, so both viscous and wind-driven models are broadly consistent with the existing observations. If accretion is viscous, the observed accretion rates require low rates of disc photoevaporation ($\lesssim$$10^{-9}$M$_{\odot}$yr$^{-1}$). Uniform, homogeneous surveys of stellar accretion rates can therefore provide a clear answer to the long-standing question of how protoplanetary discs accrete. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 10 pages, 8 figures. Accepted for publication in MNRAS

The role of the drag force in the gravitational stability of dusty planet-forming disc -- II. Numerical simulations

Authors: Cristiano Longarini, Philip J. Armitage, Giuseppe Lodato, Daniel J. Price, Simone Ceppi

Abstract: Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed p… ▽ More Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed particle hydrodynamics simulations of two-fluid discs, in the regime where the Stokes number of the particles St>1, to study how the formation of solid clumps depends on the disc-to-star mass ratio, the strength of gravitational instability, and the Stokes number. Gravitational instability of the simulated discs is sustained by local cooling. We find that the ability of the spiral structures to concentrate solids increases with the cooling time, and decreases with the Stokes number, while the relative dynamical temperature between gas and dust of the particles decreases with the cooling time and the disc-to-star mass ratio, and increases with the Stokes number. Dust collapse occurs in a subset of high disc mass simulations, yielding clumps whose mass is close to linear theory estimates, namely 1-10 Earth masses. Our results suggest that if planet formation occurs via this mechanism, the best conditions correspond to near the end of the self-gravitating phase, when the cooling time is long and the Stokes number close to unity. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS, 20 pages

Revisiting collisional dust growth in Class 0/I protostellar disks: Sweep-up can convert a few $10 M_\oplus$ of dust into kg pebbles in 0.1 Myr

Authors: Wenrui Xu, Philip J. Armitage

Abstract: Recent observations suggest that the first stages of planet formation likely take place in the Class 0/I phase of Young Stellar Object evolution, when the star and the disk are still embedded in an infalling envelope. In this study we perform grain coagulation calculations to investigate the very first stage of planet formation, the collisional growth of dust grains, in Class 0/I disks. We find th… ▽ More Recent observations suggest that the first stages of planet formation likely take place in the Class 0/I phase of Young Stellar Object evolution, when the star and the disk are still embedded in an infalling envelope. In this study we perform grain coagulation calculations to investigate the very first stage of planet formation, the collisional growth of dust grains, in Class 0/I disks. We find that the slow increase in grain mass by high-velocity collision with much smaller grains ("sweep-up") allows $\sim 50 M_\oplus$ of grains to grow well beyond the fragmentation barrier into $\sim$kg pebbles by the end of Class 0/I (0.1 Myr). We analyze the linear growth and saturation of sweep-up to understand our results quantitatively, and test whether the sweep-up outcome is sensitive to disk parameters and details of the grain coagulation model. The sweep-up pebble population could be important for planet formation, because they are less well-coupled to the gas (compared to the main population below the fragmentation barrier) and therefore more favorable to known mechanisms of dust clump formation (which initiate planetesimal formation). It also contains enough mass to form all planet cores, based on observational estimates of the planet mass budget. Our findings motivate future studies of grain growth and planetesimal formation in Class 0/I disks, including the subsequent evolution of this sweep-up population. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 16 pages, 10 figures. Accepted for publication in ApJ

The role of the drag force in the gravitational stability of dusty planet forming disc -- I. Analytical theory

Authors: Cristiano Longarini, Giuseppe Lodato, Giuseppe Bertin, Philip J. Armitage

Abstract: Recent observations show that planet formation is already underway in young systems, when the protostar is still embedded into the molecular cloud and the accretion disc is massive. In such environments, the role of self gravity (SG) and gravitational instability (GI) is crucial in determining the dynamical evolution of the disc. In this work, we study the dynamical role of drag force in self-grav… ▽ More Recent observations show that planet formation is already underway in young systems, when the protostar is still embedded into the molecular cloud and the accretion disc is massive. In such environments, the role of self gravity (SG) and gravitational instability (GI) is crucial in determining the dynamical evolution of the disc. In this work, we study the dynamical role of drag force in self-gravitating discs as a way to form planetesimals in early protoplanetary stages. We obtain the dispersion relation for density-wave perturbations on a fluid composed of two phases (gas and dust) interacting through the common gravitation field and the mutual drag force, and we find that the stability threshold is determined by three parameters: the local dust-to-gas density ratio, the dust relative temperature and the relevant Stokes number. In a region of parameters space, where young protoplanetary discs are likely to be found, the instability can be dust driven, occurring at small wavelengths. In this regime, the Jeans mass is much smaller than the one predicted by the standard gravitational instability model. This mechanism can be a viable way to form planetary cores in protostellar discs, since their predicted mass is about 10 Earth masses. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: 14 pages, 8 figures, 1 appendix, accepted for publication on MNRAS

Saturation of the magnetorotational instability and the origin of magnetically elevated accretion discs

Authors: Mitchell C. Begelman, Philip J. Armitage

Abstract: We propose that the strength of angular momentum transport in accretion discs threaded by net vertical magnetic field is determined by a self-regulation mechanism: the magnetorotational instability (MRI) grows until its own turbulent resistivity damps the fastest growing mode on the scale of the disc thickness. Given weak assumptions as to the structure of MRI-derived turbulence, supported by prio… ▽ More We propose that the strength of angular momentum transport in accretion discs threaded by net vertical magnetic field is determined by a self-regulation mechanism: the magnetorotational instability (MRI) grows until its own turbulent resistivity damps the fastest growing mode on the scale of the disc thickness. Given weak assumptions as to the structure of MRI-derived turbulence, supported by prior simulation evidence, the proposed mechanism reproduces the known scaling of the viscous $α$-parameter, $α\propto β_z^{-1/2}$. Here, $β_z = 8πp_g/B_{z0}^2$ is the initial plasma $β$-parameter on the disc midplane, $B_{z0}$ is the net field, and $p_g $ is the midplane gas pressure. We generalize the argument to discs with strong suprathermal toroidal magnetic fields, where the MRI growth rate is modified from the weak-field limit. Additional sources of turbulence are required if such discs are to become magnetically elevated, with the increased scale heights near the midplane that are seen in simulations. We speculate that tearing modes, associated with current sheets broadened by the effective resistivity, are a possible source of enhanced turbulence in elevated discs. △ Less

Submitted 16 April, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: 9 pages, 2 figures, to be published in Monthly Notices of the Royal Astronomical Society. Replaces arXiv:2211.09261, eliminating section 2.2 and with moderate revisions and clarifications in sections 3.1, 3.3, 4.2 and 4.3

A Neural Network Subgrid Model of the Early Stages of Planet Formation

Authors: Thomas Pfeil, Miles Cranmer, Shirley Ho, Philip J. Armitage, Tilman Birnstiel, Hubert Klahr

Abstract: Planet formation is a multi-scale process in which the coagulation of $\mathrm{μm}$-sized dust grains in protoplanetary disks is strongly influenced by the hydrodynamic processes on scales of astronomical units ($\approx 1.5\times 10^8 \,\mathrm{km}$). Studies are therefore dependent on subgrid models to emulate the micro physics of dust coagulation on top of a large scale hydrodynamic simulation.… ▽ More Planet formation is a multi-scale process in which the coagulation of $\mathrm{μm}$-sized dust grains in protoplanetary disks is strongly influenced by the hydrodynamic processes on scales of astronomical units ($\approx 1.5\times 10^8 \,\mathrm{km}$). Studies are therefore dependent on subgrid models to emulate the micro physics of dust coagulation on top of a large scale hydrodynamic simulation. Numerical simulations which include the relevant physical effects are complex and computationally expensive. Here, we present a fast and accurate learned effective model for dust coagulation, trained on data from high resolution numerical coagulation simulations. Our model captures details of the dust coagulation process that were so far not tractable with other dust coagulation prescriptions with similar computational efficiency. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 6 pages, 4 figures, accepted at the Machine Learning and the Physical Sciences workshop, NeurIPS 2022

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

Authors: Yan-Mong Chan, Natascha Manger, Yin Li, Chao-Chin Yang, Zhaohuan Zhu, Philip J. Armitage, Shirley Ho

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

Submitted 6 January, 2024; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 20 pages, 13 figures; Accepted to ApJ on 2023 October 30

Journal ref: ApJ 960 19 (2024)

Lifetime of the Outer Solar System Nebula From Carbonaceous Chondrites

Authors: C. S. Borlina, B. P. Weiss, J. F. J. Bryson, P. J. Armitage

Abstract: The evolution and lifetime of protoplanetary disks (PPDs) play a central role in the formation and architecture of planetary systems. Astronomical observations suggest that PPDs evolve in two timescales, accreting onto the star for up to several million years (Myr) followed by gas dissipation within <1 Myr. Because solar nebula magnetic fields are sustained by the gas of the protoplanetary disk, w… ▽ More The evolution and lifetime of protoplanetary disks (PPDs) play a central role in the formation and architecture of planetary systems. Astronomical observations suggest that PPDs evolve in two timescales, accreting onto the star for up to several million years (Myr) followed by gas dissipation within <1 Myr. Because solar nebula magnetic fields are sustained by the gas of the protoplanetary disk, we can use paleomagnetic measurements to infer the lifetime of the solar nebula. Here, we use paleomagnetic measurements of meteorites to constrain this lifetime and investigate whether the solar nebula had a two-timescale evolution. We report on paleomagnetic measurements of bulk subsamples of two CO carbonaceous chondrites: Allan Hills A77307 and Dominion Range 08006. If magnetite in these meteorites can acquire a crystallization remanent magnetization that recorded the ambient field during aqueous alteration, our measurements suggest that the local magnetic field strength at the CO parent body location was <0.9 \muT at some time between 2.7 and 5.1 Myr after the formation of calcium-aluminum-rich inclusions. Coupled with previous paleomagnetic studies, we conclude that the dissipation of the solar nebula in the 3-7 AU region occurred <1.5 Myr after the dissipation of the nebula in the 1-3 AU region, suggesting that protoplanetary disks go through a two-timescale evolution in their lifetime, consistent with dissipation by photoevaporation and/or magnetohydrodynamic winds. We also discuss future directions necessary to obtain robust records of solar nebula fields using bulk chondrites, including obtaining ages from meteorites and experimental work to determine how magnetite acquires magnetization during chondrite parent body alteration. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Comments: 30 pages, 7 figures, 4 tables

Journal ref: JGR: Planets (2022)

Demographics of young stars and their protoplanetary disks: lessons learned on disk evolution and its connection to planet formation

Authors: C. F. Manara, M. Ansdell, G. P. Rosotti, A. M. Hughes, P. J. Armitage, G. Lodato, J. P. Williams

Abstract: Since Protostars and Planets VI (PPVI), our knowledge of the global properties of protoplanetary and debris disks, as well as of young stars, has dramatically improved. At the time of PPVI, mm-observations and optical to near-infrared spectroscopic surveys were largely limited to the Taurus star-forming region, especially of its most massive disk and stellar population. Now, near-complete surveys… ▽ More Since Protostars and Planets VI (PPVI), our knowledge of the global properties of protoplanetary and debris disks, as well as of young stars, has dramatically improved. At the time of PPVI, mm-observations and optical to near-infrared spectroscopic surveys were largely limited to the Taurus star-forming region, especially of its most massive disk and stellar population. Now, near-complete surveys of multiple star-forming regions cover both spectroscopy of young stars and mm interferometry of their protoplanetary disks. This provides an unprecedented statistical sample of stellar masses and mass accretion rates, as well as disk masses and radii, for almost 1000 young stellar objects within 300 pc from us, while also sampling different evolutionary stages, ages, and environments. At the same time, surveys of debris disks are revealing the bulk properties of this class of more evolved objects. This chapter reviews the statistics of these measured global star and disk properties and discusses their constraints on theoretical models describing global disk evolution. Our comparisons of observations to theoretical model predictions extends beyond the traditional viscous evolution framework to include analytical descriptions of magnetic wind effects. Finally, we discuss how recent observational results can provide a framework for models of planet population synthesis and planet formation. △ Less

Submitted 1 March, 2023; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: Review Chapter for Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Final accepted version after interaction with the referees and community feedback. 21 pages (24 with references), 8 figures. Data table available at http://ppvii.org/chapter/15/

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

Authors: Matthew M. Murphy, Philip J. Armitage

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

Submitted 16 March, 2022; originally announced March 2022.

Comments: 9 pages, 9 figures, Accepted to MNRAS

Lecture notes on accretion disk physics

Authors: Philip J. Armitage

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

Submitted 18 January, 2022; originally announced January 2022.

Comments: Expanded write-up of introductory graduate level lectures

MHD disc winds can reproduce fast disc dispersal and the correlation between accretion rate and disc mass in Lupus

Authors: Benoît Tabone, Giovanni P. Rosotti, Giuseppe Lodato, Philip J. Armitage, Alexander J. Cridland, Ewine F. van Dishoeck

Abstract: The final architecture of planetary systems depends on the extraction of angular momentum and mass-loss processes of the discs in which they form. Theoretical studies proposed that magnetized winds launched from the discs (MHD disc winds) could govern accretion and disc dispersal. In this work, we revisit the observed disc demographics in the framework of MHD disc winds, combining analytical solut… ▽ More The final architecture of planetary systems depends on the extraction of angular momentum and mass-loss processes of the discs in which they form. Theoretical studies proposed that magnetized winds launched from the discs (MHD disc winds) could govern accretion and disc dispersal. In this work, we revisit the observed disc demographics in the framework of MHD disc winds, combining analytical solutions of disc evolution and a disc population synthesis approach. We show that MHD disc winds alone can account for both disc dispersal and accretion properties. The decline of disc fraction over time is reproduced by assuming that the initial accretion timescale (a generalization of the viscous timescale) varies from disc to disc and that the decline of the magnetic field strength is slower than that of the gas. The correlation between accretion rate and disc mass, and the dispersion of the data around the mean trend as observed in Lupus is then naturally reproduced. The model also accounts for the rapidity of the disc dispersal. This paves the way for planet formation models in the paradigm of wind-driven accretion. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS Letters

Secular evolution of MHD wind-driven discs: analytical solutions in the expanded α-framework

Authors: Benoît Tabone, Giovanni P. Rosotti, Alexander J. Cridland, Philip J. Armitage, Giuseppe Lodato

Abstract: The evolution of protoplanetary discs and the related process of planet formation is regulated by angular momentum transport and mass-loss processes. Over the past decade, the paradigm of viscosity has been challenged and MHD disc winds appear as a compelling scenario to account for disc accretion. In this work, we aim to construct the equivalent of the widely used analytical description of viscou… ▽ More The evolution of protoplanetary discs and the related process of planet formation is regulated by angular momentum transport and mass-loss processes. Over the past decade, the paradigm of viscosity has been challenged and MHD disc winds appear as a compelling scenario to account for disc accretion. In this work, we aim to construct the equivalent of the widely used analytical description of viscous evolution for the MHD wind case. The transport of angular momentum and mass induced by the wind is parameterized by an $α$-like parameter and by the magnetic lever arm parameter $λ$. Extensions of the paradigmatic Lynden-Bell and Pringle similarity solutions to the wind case are presented. We show that wind-driven accretion leads to a steeper decrease in the disc mass and accretion rate than in viscous models due to the absence of disc spreading. If the decline of the magnetic field strength is slower than that of the gas surface density, the disc is dispersed after a finite time. The evolution of the disc in the $\dot{M}_*-M_D$ plane is sensitive to the wind and turbulence parameters. A disc population evolving under the action of winds can exhibit a correlation between $\dot{M}_*$ and $M_D$ depending on the initial conditions. The simplified framework proposed in this work opens to a new avenue to test the effectiveness of wind-driven accretion from the observed disc demographics and constitutes an important step to include wind-driven accretion in planet population synthesis models. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

Dynamical stability of giant planets: the critical adiabatic index in the presence of a solid core

Authors: Suman Kumar Kundu, Eric R. Coughlin, Andrew N. Youdin, Philip J. Armitage

Abstract: The dissociation and ionization of hydrogen, during the formation of giant planets via core accretion, reduces the effective adiabatic index $γ$ of the gas and could trigger dynamical instability. We generalize the analysis of Chandrasekhar, who determined that the threshold for instability of a self-gravitating hydrostatic body lies at $γ=4/3$, to account for the presence of a planetary core, whi… ▽ More The dissociation and ionization of hydrogen, during the formation of giant planets via core accretion, reduces the effective adiabatic index $γ$ of the gas and could trigger dynamical instability. We generalize the analysis of Chandrasekhar, who determined that the threshold for instability of a self-gravitating hydrostatic body lies at $γ=4/3$, to account for the presence of a planetary core, which we model as an incompressible fluid. We show that the dominant effect of the core is to stabilize the envelope to radial perturbations, in some cases completely (i.e. for all $γ> 1$). When instability is possible, unstable planetary configurations occupy a strip of $γ$ values whose upper boundary falls below $γ=4/3$. Fiducial evolutionary tracks of giant planets forming through core accretion appear unlikely to cross the dynamical instability strip that we define. △ Less

Submitted 6 August, 2021; originally announced August 2021.

Global 3D Radiation Hydrodynamic Simulations of Proto-Jupiter's Convective Envelope

Authors: Zhaohuan Zhu, Yan-Fei Jiang, Hans Baehr, Andrew N. Youdin, Philip J. Armitage, Rebecca G. Martin

Abstract: The core accretion model of giant planet formation has been challenged by the discovery of recycling flows between the planetary envelope and the disc that can slow or stall envelope accretion. We carry out 3D radiation hydrodynamic simulations with an updated opacity compilation to model the proto-Jupiter's envelope. To isolate the 3D effects of convection and recycling, we simulate both isolated… ▽ More The core accretion model of giant planet formation has been challenged by the discovery of recycling flows between the planetary envelope and the disc that can slow or stall envelope accretion. We carry out 3D radiation hydrodynamic simulations with an updated opacity compilation to model the proto-Jupiter's envelope. To isolate the 3D effects of convection and recycling, we simulate both isolated spherical envelopes and envelopes embedded in discs. The envelopes are heated at given rates to achieve steady states, enabling comparisons with 1D models. We vary envelope properties to obtain both radiative and convective solutions. Using a passive scalar, we observe significant mass recycling on the orbital timescale. For a radiative envelope, recycling can only penetrate from the disc surface until $\sim$0.1-0.2 planetary Hill radii, while for a convective envelope, the convective motion can "dredge up" the deeper part of the envelope so that the entire convective envelope is recycled efficiently. This recycling, however, has only limited effects on the envelopes' thermal structure. The radiative envelope embedded in the disc has identical structure as the isolated envelope. The convective envelope has a slightly higher density when it is embedded in the disc. We introduce a modified 1D approach which can fully reproduce our 3D simulations. With our updated opacity and 1D model, we recompute Jupiter's envelope accretion with a 10 $M_{\oplus}$ core, and the timescale to runaway accretion is shorter than the disc lifetime as in prior studies. Finally, we discuss the implications of the efficient recycling on the observed chemical abundances of the planetary atmosphere (especially for super-Earths and mini-Neptunes). △ Less

Submitted 1 September, 2021; v1 submitted 22 June, 2021; originally announced June 2021.

Comments: 24 pages, 17 figures, Accepted for publication in MNRAS

Primordial giant planet obliquity driven by a circumplanetary disk

Authors: Rebecca G. Martin, Philip J. Armitage

Abstract: Detached circumplanetary disks are unstable to tilting as a result of the stellar tidal potential. We examine how a tilted circumplanetary disk affects the evolution of the spin axis of an oblate planet. The disk is evolved using time-dependent equations for linear wave-like warp evolution, including terms representing the effect of the tidal potential and planetary oblateness. For a disk with a s… ▽ More Detached circumplanetary disks are unstable to tilting as a result of the stellar tidal potential. We examine how a tilted circumplanetary disk affects the evolution of the spin axis of an oblate planet. The disk is evolved using time-dependent equations for linear wave-like warp evolution, including terms representing the effect of the tidal potential and planetary oblateness. For a disk with a sufficiently large mass, we find that the planet spin quickly aligns to the misaligned disk. The tilt of the planetary spin axis then increases on the same timescale as the disk. This can be an efficient mechanism for generating primordial obliquity in giant planets. We suggest that directly imaged exoplanets at large orbital radii, where the disk mass criterion is more likely to be satisfied, could have significant obliquities due to the tilt instability of their circumplanetary disks. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: Accepted for publication in ApJL

Partial tidal disruption events by stellar mass black holes: gravitational instability of stream and impact from remnant core

Authors: Yi-Han Wang, Rosalba Perna, Philip J. Armitage

Abstract: In dense star clusters, such as globular and open clusters, dynamical interactions between stars and black holes (BHs) can be extremely frequent, leading to various astrophysical transients. Close encounters between a star and a stellar mass BH make it possible for the star to be tidally disrupted by the BH. Due to the relative low mass of the BH and the small cross section of the tidal disruption… ▽ More In dense star clusters, such as globular and open clusters, dynamical interactions between stars and black holes (BHs) can be extremely frequent, leading to various astrophysical transients. Close encounters between a star and a stellar mass BH make it possible for the star to be tidally disrupted by the BH. Due to the relative low mass of the BH and the small cross section of the tidal disruption event (TDE) for cases with high penetration, disruptions caused by close encounters are usually partial disruptions. The existence of the remnant stellar core and its non-negligible mass compared to the stellar mass BH alters the accretion process significantly. We study this problem with SPH simulations using the code {\tt Phantom}, with the inclusion of radiation pressure, which is important for small mass BHs. Additionally, we develop a new, more general method of computing the fallback rate which does not rely on any approximation. Our study shows that the powerlaw slope of the fallback rate has a strong dependence on the mass of the BH in the stellar mass BH regime. Furthermore, in this regime, self-gravity of the fallback stream and local instabilities become more significant, and cause the disrupted material to collapse into small clumps before returning to the BH. This results in an abrupt increase of the fallback rate, which can significantly deviate from a powerlaw. Our results will help in the identification of TDEs by stellar mass BHs in dense clusters. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS

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

Authors: Rebecca G. Martin, Zhaohuan Zhu, Philip J. Armitage, Chao-Chin Yang, Hans Baehr

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

Submitted 22 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS

A Bayesian neural network predicts the dissolution of compact planetary systems

Authors: Miles Cranmer, Daniel Tamayo, Hanno Rein, Peter Battaglia, Samuel Hadden, Philip J. Armitage, Shirley Ho, David N. Spergel

Abstract: Despite over three hundred years of effort, no solutions exist for predicting when a general planetary configuration will become unstable. We introduce a deep learning architecture to push forward this problem for compact systems. While current machine learning algorithms in this area rely on scientist-derived instability metrics, our new technique learns its own metrics from scratch, enabled by a… ▽ More Despite over three hundred years of effort, no solutions exist for predicting when a general planetary configuration will become unstable. We introduce a deep learning architecture to push forward this problem for compact systems. While current machine learning algorithms in this area rely on scientist-derived instability metrics, our new technique learns its own metrics from scratch, enabled by a novel internal structure inspired from dynamics theory. Our Bayesian neural network model can accurately predict not only if, but also when a compact planetary system with three or more planets will go unstable. Our model, trained directly from short N-body time series of raw orbital elements, is more than two orders of magnitude more accurate at predicting instability times than analytical estimators, while also reducing the bias of existing machine learning algorithms by nearly a factor of three. Despite being trained on compact resonant and near-resonant three-planet configurations, the model demonstrates robust generalization to both non-resonant and higher multiplicity configurations, in the latter case outperforming models fit to that specific set of integrations. The model computes instability estimates up to five orders of magnitude faster than a numerical integrator, and unlike previous efforts provides confidence intervals on its predictions. Our inference model is publicly available in the SPOCK package, with training code open-sourced. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: 8 content pages, 7 appendix and references. 8 figures. Source code at: https://github.com/MilesCranmer/bnn_chaos_model inference code at https://github.com/dtamayo/spock

Survivor bias: divergent fates of the Solar System's ejected vs. persisting planetesimals

Authors: Sean N. Raymond, Nathan A. Kaib, Philip J. Armitage, Jonathan J. Fortney

Abstract: The orbital architecture of the Solar System is thought to have been sculpted by a dynamical instability among the giant planets. During the instability a primordial outer disk of planetesimals was destabilized and ended up on planet-crossing orbits. Most planetesimals were ejected into interstellar space but a fraction were trapped on stable orbits in the Kuiper belt and Oort cloud. We use a suit… ▽ More The orbital architecture of the Solar System is thought to have been sculpted by a dynamical instability among the giant planets. During the instability a primordial outer disk of planetesimals was destabilized and ended up on planet-crossing orbits. Most planetesimals were ejected into interstellar space but a fraction were trapped on stable orbits in the Kuiper belt and Oort cloud. We use a suite of N-body simulations to map out the diversity of planetesimals' dynamical pathways. We focus on two processes: tidal disruption from very close encounters with a giant planet, and loss of surface volatiles from repeated passages close to the Sun. We show that the rate of tidal disruption is more than a factor of two higher for ejected planetesimals than for surviving objects in the Kuiper belt or Oort cloud. Ejected planetesimals are preferentially disrupted by Jupiter and surviving ones by Neptune. Given that the gas giants contracted significantly as they cooled but the ice giants did not, taking into account the thermal evolution of the giant planets decreases the disruption rate of ejected planetesimals. The frequency of volatile loss and extinction is far higher for ejected planetesimals than for surviving ones and is not affected by the giant planets' contraction. Even if all interstellar objects were ejected from Solar System-like systems, our analysis suggests that their physical properties should be more diverse than those of Solar System small bodies as a result of their divergent dynamical histories. This is consistent with the characteristics of the two currently-known interstellar objects. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: ApJ Letters, in press. 6 pages, 4 figures

Visualizing the Kinematics of Planet Formation

Authors: Disk Dynamics Collaboration, Philip J. Armitage, Jaehan Bae, Myriam Benisty, Edwin A. Bergin, Simon Casassus, Ian Czekala, Stefano Facchini, Jeffrey Fung, Cassandra Hall, John D. Ilee, Miriam Keppler, Aleksandra Kuznetsova, Romane Le Gal, Ryan A. Loomis, Wladimir Lyra, Natascha Manger, Sebastian Perez, Christophe Pinte, Daniel J. Price, Giovanni Rosotti, Judit Szulagyi, Kamber Schwarz, Jacob B. Simon, Richard Teague , et al. (1 additional authors not shown)

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

Submitted 15 September, 2020; v1 submitted 9 September, 2020; originally announced September 2020.

Comments: To be submitted to PASA. Comments welcome

Predicting the long-term stability of compact multiplanet systems

Authors: Daniel Tamayo, Miles Cranmer, Samuel Hadden, Hanno Rein, Peter Battaglia, Alysa Obertas, Philip J. Armitage, Shirley Ho, David Spergel, Christian Gilbertson, Naireen Hussain, Ari Silburt, Daniel Jontof-Hutter, Kristen Menou

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

Submitted 14 July, 2020; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: Published week of July 13th in Proceedings of the National Academy of Sciences: https://www.pnas.org/cgi/doi/10.1073/pnas.2001258117. Check out simple usage of package (and regenerate paper figures) at: https://github.com/dtamayo/spock

A fast-growing tilt instability of detached circumplanetary disks

Authors: Rebecca G. Martin, Zhaohuan Zhu, Philip J. Armitage

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

Submitted 9 July, 2020; originally announced July 2020.

Comments: Accepted for publication in ApJL

Future Simulations of Tidal Disruption Events

Authors: Julian H. Krolik, Philip J. Armitage, Yanfei Jiang, Giuseppe Lodato

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

Submitted 5 June, 2020; originally announced June 2020.

Comments: A chapter in the ISSI review book "The Tidal Disruption of Stars by Massive Black Holes", to be published in Space Science Reviews

Bondi on spherically symmetric accretion

Authors: Philip J. Armitage

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

Submitted 31 March, 2020; originally announced April 2020.

Comments: Part of a series, marking the Royal Astronomical Society's bicentennial, revisiting influential papers published by the Society

Journal ref: Astronomy & Geophysics, Volume 61, 2.40-2.42 (2020)

Turbulence Regulates the Rate of Planetesimal Formation via Gravitational Collapse

Authors: Daniel A. Gole, Jacob B. Simon, Rixin Li, Andrew N. Youdin, Philip J. Armitage

Abstract: We study how the interaction between the streaming instability and intrinsic gas-phase turbulence affects planetesimal formation via gravitational collapse in protoplanetary disks. Turbulence impedes the formation of 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 ef… ▽ More We study how the interaction between the streaming instability and intrinsic gas-phase turbulence affects planetesimal formation via gravitational collapse in protoplanetary disks. Turbulence impedes the formation of 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 turbulence. We find that planetesimal formation is suppressed by turbulence once velocity fluctuations exceed $\langle δv^2 \rangle \simeq 10^{-3.5} - 10^{-3} c_s^2$. Turbulence whose strength is just below the threshold reduces the rate at which solids are bound into clumps. Our results suggest that the well-established turbulent thickening of the mid-plane solid layer is the primary mechanism by which turbulence influences planetesimal formation and that planetesimal formation requires a mid-plane solid-to-gas ratio $ε\gtrsim 0.5$. We also quantify the initial planetesimal mass function using a new clump-tracking method to determine each planetesimal mass shortly after collapse. For models in which planetesimals form, we show that the 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. Planetesimal formation may therefore be unviable in the inner disk out to 2-3 times the dust sublimation radius. △ Less

Submitted 21 October, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: Accepted to ApJ after revisions

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

Authors: Xinyu Li, Philip Chang, Yuri Levin, Christopher D. Matzner, Philip J. Armitage

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

Submitted 30 March, 2020; v1 submitted 14 December, 2019; originally announced December 2019.

Comments: Accepted by MNRAS

The Influence of Black Hole Binarity on Tidal Disruption Events

Authors: Eric R. Coughlin, Philip J. Armitage, Giuseppe Lodato, C. J. Nixon

Abstract: 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 companio… ▽ More 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. △ Less

Submitted 27 September, 2019; originally announced September 2019.

Comments: 22 pages, 5 figures. Accepted chapter for Springer Space Science Reviews book on tidal disruption events

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

Authors: Bhupendra Mishra, Mitchell C. Begelman, Philip J. Armitage, Jacob B. Simon

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

Submitted 21 January, 2020; v1 submitted 21 July, 2019; originally announced July 2019.

Envisioning the next decade of Galactic Center science: a laboratory for the study of the physics and astrophysics of supermassive black holes

Authors: Tuan Do, Andrea Ghez, Jessica R. Lu, Mark Morris, Matthew Hosek Jr., Aurelien Hees, Smadar Naoz, Anna Ciurlo, Philip J. Armitage, Rachael L Beaton, Eric Becklin, Andrea Bellini, Rory O. Bentley, Joss Bland-Hawthorn, Sukanya Chakrabarti, Zhuo Chen, Devin S. Chu, Arezu Dehghanfar, Charles F. Gammie, Abhimat K. Gautam, Reinhard Genzel, Jenny Greene, Daryl Haggard, Joseph Hora, Wolfgang E. Kerzendorf , et al. (16 additional authors not shown)

Abstract: 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 desc… ▽ More 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. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: 9 pages, 4 figures, submitted for Astro2020 White Paper

The Mass and Size Distribution of Planetesimals Formed by the Streaming Instability. II. The Effect of the Radial Gas Pressure Gradient

Authors: Charles P. Abod, Jacob B. Simon, Rixin Li, Philip J. Armitage, Andrew N. Youdin, Katherine A. Kretke

Abstract: 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 ins… ▽ More 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 π^5 G^2 Σ_p^3 / Ω^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.. △ Less

Submitted 29 August, 2019; v1 submitted 23 October, 2018; originally announced October 2018.

Comments: 18 pages, 10 figures, 3 tables, accepted to ApJ

A brief overview of planet formation

Authors: Philip J. Armitage

Abstract: 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 perva… ▽ More 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. △ Less

Submitted 28 March, 2018; originally announced March 2018.

Comments: Short introduction, aimed at new students, for a forthcoming survey volume

Interstellar object 'Oumuamua as an extinct fragment of an ejected cometary planetesimal

Authors: Sean N. Raymond, Philip J. Armitage, Dimitri Veras

Abstract: '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 'Oum… ▽ More '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. △ Less

Submitted 7 March, 2018; originally announced March 2018.

Comments: ApJ Letters, in press

Circumbinary, not transitional: On the spiral arms, cavity, shadows, fast radial flows, streamers and horseshoe in the HD142527 disc

Authors: Daniel J. Price, Nicolas Cuello, Christophe Pinte, Daniel Mentiplay, Simon Casassus, Valentin Christiaens, Grant M. Kennedy, Jorge Cuadra, Sebastian Perez M., Sebastian Marino, Philip J. Armitage, Alice Zurlo, Attila Juhasz, Enrico Ragusa, Guillaume Laibe, Giuseppe Lodato

Abstract: 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 sh… ▽ More 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. △ Less

Submitted 6 March, 2018; originally announced March 2018.

Comments: 16 pages, 12 figures, accepted to MNRAS. Movies at http://users.monash.edu.au/~dprice/pubs/HD142527/

Implications of the interstellar object 1I/'Oumuamua for planetary dynamics and planetesimal formation

Authors: Sean N. Raymond, Philip J. Armitage, Dimitri Veras, Elisa V. Quintana, Thomas Barclay

Abstract: '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 Gala… ▽ More '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. △ Less

Submitted 26 February, 2018; v1 submitted 27 November, 2017; originally announced November 2017.

Comments: MNRAS, in press

Circumbinary discs from tidal disruption events

Authors: Eric R. Coughlin, Philip J. Armitage

Abstract: 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.… ▽ More 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. △ Less

Submitted 18 July, 2017; originally announced July 2017.

Comments: 5 pages, 2 figures, MNRAS Letters accepted

A Three-Dimensional View of Turbulence: Constraints on Turbulent Motions in the HD 163296 Protoplanetary Disk using DCO$^+$

Authors: Kevin M. Flaherty, A. Meredith Hughes, Sanaea C. Rose, Jacob B. Simon, Chunhua Qi, Sean M. Andrews, Agnes Kospal, David J. Wilner, Eugene Chiang, Philip J. Armitage, Xue-ning Bai

Abstract: 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 combine… ▽ More 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 $α\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. △ Less

Submitted 14 June, 2017; originally announced June 2017.

Comments: Accepted to ApJ

Tidal disruption by extreme mass ratio binaries and application to ASASSN-15lh

Authors: Eric R. Coughlin, Philip J. Armitage

Abstract: 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… ▽ More 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. △ Less

Submitted 1 December, 2017; v1 submitted 12 May, 2017; originally announced May 2017.

Comments: 10 pages, 3 figures; updated to include revisions from referee; fixed figures

Evidence for universality in the initial planetesimal mass function

Authors: Jacob B. Simon, Philip J. Armitage, Andrew N. Youdin, Rixin Li

Abstract: 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 o… ▽ More 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 τ\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 $Δ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. △ Less

Submitted 13 September, 2017; v1 submitted 10 May, 2017; originally announced May 2017.

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

EXor outbursts from disk amplification of stellar magnetic cycles

Authors: Philip J. Armitage

Abstract: 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 thr… ▽ More 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. △ Less

Submitted 5 December, 2016; originally announced December 2016.

Comments: ApJL, in press

Tidal disruption events from supermassive black hole binaries

Authors: Eric R. Coughlin, Philip J. Armitage, Chris Nixon, Mitchell C. Begelman

Abstract: 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 an… ▽ More 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. △ Less

Submitted 15 November, 2016; v1 submitted 19 August, 2016; originally announced August 2016.

Comments: 31 pages, 25 figures, MNRAS Accepted; updated to incorporate referee comments. Movies of simulations can be found here: http://w.astro.berkeley.edu/~eric_coughlin/movies.html

Prompt planetesimal formation beyond the snow line

Authors: Philip J. Armitage, Josh A. Eisner, Jacob B. Simon

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

Submitted 11 August, 2016; originally announced August 2016.

Comments: ApJ Letters, in press

Turbulence, Transport and Waves in Ohmic Dead Zones

Authors: Daniel Gole, Jacob B. Simon, Stephen H. Lubow, Philip J. Armitage

Abstract: 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… ▽ More 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. △ Less

Submitted 9 June, 2016; originally announced June 2016.

Comments: 14 pages, 17 figures

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

Authors: Rebecca G. Martin, Stephen H. Lubow, Chris Nixon, Philip J. Armitage

Abstract: 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 sta… ▽ More 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). △ Less

Submitted 9 March, 2016; originally announced March 2016.

Comments: Accepted for publication in MNRAS

On the structure of tidally-disrupted stellar debris streams

Authors: Eric R. Coughlin, Chris Nixon, Mitchell C. Begelman, Philip J. Armitage

Abstract: 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… ▽ More 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. △ Less

Submitted 2 March, 2016; originally announced March 2016.

Comments: 17 pages, 10 figures, Resubmitted to MNRAS after first referee report

Strongly magnetized accretion discs require poloidal flux

Authors: Greg Salvesen, Philip J. Armitage, Jacob B. Simon, Mitchell C. Begelman

Abstract: 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 wit… ▽ More 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. △ Less

Submitted 18 May, 2016; v1 submitted 15 February, 2016; originally announced February 2016.

Comments: 7 pages, 3 figures, 1 table; Accepted to MNRAS