Showing 1–50 of 105 results for author: Rafikov, R R

Early stages of gap opening by planets in protoplanetary discs

Authors: Amelia J. Cordwell, Roman R. Rafikov

Abstract: Annular substructures in protoplanetary discs, ubiquitous in sub-mm observations, can be caused by gravitational coupling between a disc and its embedded planets. Planetary density waves inject angular momentum into the disc leading to gap opening only after travelling some distance and steepening into shocks (in the absence of linear damping); no angular momentum is deposited in the planetary coo… ▽ More Annular substructures in protoplanetary discs, ubiquitous in sub-mm observations, can be caused by gravitational coupling between a disc and its embedded planets. Planetary density waves inject angular momentum into the disc leading to gap opening only after travelling some distance and steepening into shocks (in the absence of linear damping); no angular momentum is deposited in the planetary coorbital region, where the wave has not shocked yet. Despite that, simulations show mass evacuation from the coorbital region even in inviscid discs, leading to smooth, double-trough gap profiles. Here we consider the early, time-dependent stages of planetary gap opening in inviscid discs. We find that an often-overlooked contribution to the angular momentum balance caused by the time-variability of the specific angular momentum of the disc fluid (caused, in turn, by the time-variability of the radial pressure support) plays a key role in gap opening. Focusing on the regime of shallow gaps with depths of $\lesssim 20\%$, we demonstrate analytically that early gap opening is a self-similar process, with the amplitude of the planet-driven perturbation growing linearly in time and the radial gap profile that can be computed semi-analytically. We show that mass indeed gets evacuated from the coorbital region even in inviscid discs. This evolution pattern holds even in viscous discs over a limited period of time. These results are found to be in excellent agreement with 2D numerical simulations. Our simple gap evolution solutions can be used in studies of dust dynamics near planets and for interpreting protoplanetary disc observations. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 20 pages, 16 figures, submitted to MNRAS, comments welcome

Orbital evolution of LIGO/Virgo binaries in stellar clusters driven by cluster tides, stellar encounters and general relativity

Authors: Alexander Rasskazov, Roman R. Rafikov

Abstract: Origin of LIGO/Virgo gravitational wave events may involve production of binaries with relativistic components in dense stellar systems - globular or nuclear star clusters - and their subsequent evolution towards merger. Orbital parameters of these binaries (the inner orbit) and their motion inside the cluster (the outer orbit) evolve due to both external agents - random encounters with cluster st… ▽ More Origin of LIGO/Virgo gravitational wave events may involve production of binaries with relativistic components in dense stellar systems - globular or nuclear star clusters - and their subsequent evolution towards merger. Orbital parameters of these binaries (the inner orbit) and their motion inside the cluster (the outer orbit) evolve due to both external agents - random encounters with cluster stars and cluster tides due to the smooth cluster potential - and the internal ones - various sources of dissipation and precession within the binary. We present a numerical framework - Binary Evolution in Stellar Clusters (BESC) - that follows the evolution of the binary inner and outer orbits accounting for all these effects simultaneously, enabling efficient Monte Carlo studies. The secular effect of cluster tides is computed in the singly-averaged approximation, without averaging over the outer binary orbit. As to stellar encounters, we include the effects of both close and distant flybys on the inner and outer orbits of the binary, respectively. In particular, this allows us to explicitly account for the dynamical friction sinking the binary towards the cluster centre. Also, given our focus on the LIGO/Virgo sources, we include the general relativistic precession (which suppresses cluster tides at high eccentricities) and the gravitational wave emission (shrinking the binary orbit). We use BESC to illustrate a number of characteristic binary evolutionary outcomes and discuss relative contributions of different physical processes. BESC can also be used to study other objects in clusters, e.g. blue stragglers, hot Jupiters, X-ray binaries, etc. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 20 pages, 16 figures, submitted to MNRAS

Gravitational torque in circumbinary discs: global radial oscillations

Authors: Nicolas P. Cimerman, Roman R. Rafikov

Abstract: Circumbinary discs (CBDs) arise in many astrophysical settings, including young stellar binaries and supermassive black hole binaries. Their structure is mediated by gravitational torques exerted on the disc by the central binary. The spatial distribution of the binary torque density (so-called excitation torque density) in CBDs is known to feature global large-amplitude, quasi-periodic oscillatio… ▽ More Circumbinary discs (CBDs) arise in many astrophysical settings, including young stellar binaries and supermassive black hole binaries. Their structure is mediated by gravitational torques exerted on the disc by the central binary. The spatial distribution of the binary torque density (so-called excitation torque density) in CBDs is known to feature global large-amplitude, quasi-periodic oscillations, which are often interpreted in terms of the local resonant Lindblad torques. Here we investigate the nature of these torque oscillations using 2D, inviscid hydrodynamic simulations and theoretical calculations. We show that torque oscillations arise due to the gravitational coupling of the binary potential to the density waves launched near the inner cavity and freely propagating out in the disc. We provide analytical predictions for the radial periodicity of the torque density oscillations and verify them with simulations, showing that disc sound speed and the multiplicity of the density wave spiral arms are the key factors setting the radial structure of the oscillations. Resonant Lindblad torques play no direct role in determining the radial structure and periodicity of the torque oscillations and manifest themselves only by driving the density waves in the disc. We also find that vortices forming at the inner edge of the disc can provide a non-trivial contribution to the angular momentum transport in the CBD. Our results can be applied to understanding torque behaviour in other settings, e.g. discs in cataclysmic variables and X-ray binaries. △ Less

Submitted 3 August, 2023; originally announced August 2023.

Comments: 20 pages, 15 figures, submitted to MNRAS

Torque wiggles -- a robust feature of the global disc-planet interaction

Authors: Nicolas P. Cimerman, Roman R. Rafikov, Ryan Miranda

Abstract: Gravitational coupling between planets and protoplanetary discs is responsible for many important phenomena such as planet migration and gap formation. The key quantitative characteristics of this coupling is the excitation torque density -- the torque (per unit radius) imparted on the disc by planetary gravity. Recent global simulations and linear calculations found an intricate pattern of low-am… ▽ More Gravitational coupling between planets and protoplanetary discs is responsible for many important phenomena such as planet migration and gap formation. The key quantitative characteristics of this coupling is the excitation torque density -- the torque (per unit radius) imparted on the disc by planetary gravity. Recent global simulations and linear calculations found an intricate pattern of low-amplitude, quasi-periodic oscillations in the global radial distribution of torque density in the outer disc, which we call torque wiggles. Here we show that torque wiggles are a robust outcome of global disc-planet interaction and exist despite the variation of disc parameters and thermodynamic assumptions (including $β$-cooling). They result from coupling of the planetary potential to the planet-driven density wave freely propagating in the disc. We developed analytical theory of this phenomenon based on approximate self-similarity of the planet-driven density waves in the outer disc. We used it, together with linear calculations and simulations, to show that (a) the radial periodicity of the wiggles is determined by the global shape of the planet-driven density wave (its wrapping in the disc) and (b) the sharp features in the torque density distribution result from constructive interference of different azimuthal (Fourier) torque contributions at radii where the planetary wake crosses the star-planet line. In the linear regime the torque wiggles represent a weak effect, affecting the total (integrated) torque by only a few per cent. However, their significance should increase in the non-linear regime, when a gap (or a cavity) forms around the perturber's orbit. △ Less

Submitted 12 June, 2023; originally announced June 2023.

Comments: 19 pages, 15 figures, submitted to MNRAS

Relativistic phase space diffusion of compact object binaries in stellar clusters and hierarchical triples

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: The LIGO/Virgo detections of compact object mergers have posed a challenge for theories of binary evolution and coalescence. One promising avenue for producing mergers dynamically is through secular eccentricity oscillations driven by an external perturber, be it a tertiary companion (as in the Lidov-Kozai (LK) mechanism) or the tidal field of the stellar cluster in which the binary orbits. The si… ▽ More The LIGO/Virgo detections of compact object mergers have posed a challenge for theories of binary evolution and coalescence. One promising avenue for producing mergers dynamically is through secular eccentricity oscillations driven by an external perturber, be it a tertiary companion (as in the Lidov-Kozai (LK) mechanism) or the tidal field of the stellar cluster in which the binary orbits. The simplest theoretical models of these oscillations use a 'doubly-averaged' (DA) approximation, averaging both over the binary's internal Keplerian orbit and its 'outer' barycentric orbit relative to the perturber. However, DA theories do not account for fluctuations of the perturbing torque on the outer orbital timescale, which are known to increase a binary's eccentricity beyond the maximum DA value, potentially accelerating mergers. Here we reconsider the impact of these short-timescale fluctuations in the test-particle quadrupolar limit for binaries perturbed by arbitrary spherical cluster potentials (including LK as a special case), {in particular including 1pN} general relativistic (GR) apsidal precession of the internal orbit. Focusing on the behavior of the binary orbital elements around peak eccentricity, we discover a new effect, relativistic phase space diffusion (RPSD), in which a binary can jump to a completely new dynamical trajectory on an outer orbital timescale, violating the approximate conservation of DA integrals of motion. RPSD arises from an interplay between secular behavior at extremely high eccentricity, short-timescale fluctuations, and rapid GR precession, and can change the subsequent secular evolution dramatically. This effect occurs even in hierarchical triples, but has not been uncovered until now. △ Less

Submitted 4 December, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: Revised version, accepted for publication in ApJ

Modeling planet-induced gaps and rings in ALMA disks: the role of in-plane radiative diffusion

Authors: Alexandros Ziampras, Richard P. Nelson, Roman R. Rafikov

Abstract: ALMA observations of protoplanetary disks in dust continuum emission reveal a variety of annular structures. Attributing the existence of such features to embedded planets is a popular scenario, supported by studies using hydrodynamical models. Recent work has shown that radiative cooling greatly influences the capability of planet-driven spiral density waves to transport angular momentum, ultimat… ▽ More ALMA observations of protoplanetary disks in dust continuum emission reveal a variety of annular structures. Attributing the existence of such features to embedded planets is a popular scenario, supported by studies using hydrodynamical models. Recent work has shown that radiative cooling greatly influences the capability of planet-driven spiral density waves to transport angular momentum, ultimately deciding the number, position, and depth of rings and gaps that a planet can carve in a disk. However, radiation transport has only been treated via local thermal relaxation, not taking into account radiative diffusion along the disk plane. We compare the previous state-of-the-art models of planet-disk interaction with local cooling prescriptions to our new models that include cooling in the vertical direction and radiative diffusion in the plane of the disk, and show that the response of the disk to the induced spiral waves can differ significantly when comparing these two treatments of the disk thermodynamics. We follow up with synthetic emission maps of ALMA systems, and show that our new models reproduce the observations found in the literature better than models with local cooling. We conclude that appropriate treatment of radiation transport is key to constraining the parameter space when interpreting ALMA observations using the planet-disk interaction scenario. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 18 pages, 28 figures, 4 tables; submitted to MNRAS

Formation of Gaps in Self-gravitating Debris Disks by Secular Resonance in a Single-planet System. II. Towards a Self-consistent Model

Authors: Antranik A. Sefilian, Roman R. Rafikov, Mark C. Wyatt

Abstract: High-resolution observations of several debris disks reveal structures such as gaps and spirals, suggestive of gravitational perturbations induced by underlying planets. Most existing studies of planet--debris disk interactions ignore the gravity of the disk, treating it as a reservoir of massless planetesimals. In this paper, we continue our investigation into the long-term interaction between a… ▽ More High-resolution observations of several debris disks reveal structures such as gaps and spirals, suggestive of gravitational perturbations induced by underlying planets. Most existing studies of planet--debris disk interactions ignore the gravity of the disk, treating it as a reservoir of massless planetesimals. In this paper, we continue our investigation into the long-term interaction between a single eccentric planet and an external, massive debris disk. Building upon our previous work, here we consider not only the axisymmetric component of the disk's gravitational potential, but also the non-axisymmetric torque that the disk exerts on the planet (ignoring for now only the non-axisymmetric component of the disk self-gravity). To this goal, we develop and test a semi-analytic `$N$-ring' framework that is based on a generalized (softened) version of the classical Laplace--Lagrange secular theory. Using this tool, we demonstrate that even when the disk is less massive than the planet, not only can a secular resonance be established within the disk that leads to the formation of a wide gap, but that the very same resonance also damps the planetary eccentricity $e_p$ via a process known as resonant friction. The resulting gap is initially non-axisymmetric (akin to those observed in HD 92945 and HD 206893), but evolves to become more axisymmetric (similar to that in HD 107146) as $e_p(t)\rightarrow0$ with time. We also develop analytic understanding of these findings, finding good quantitative agreement with the outcomes of the $N$-ring calculations. Our results may be used to infer both the dynamical masses of (gapped) debris disks and the dynamical history of the planets interior to them, as we exemplify for HD 206893. △ Less

Submitted 14 July, 2023; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: Accepted version (ApJ): 35 Pages (including 8 pages of Appendices), 15 Figures, 1 Table, 3 Animations (see Ancillary files). A copy of the N-ring code is made publicly available at http://doi.org/10.6084/m9.figshare.22598377

Vortex weighing and dating of planets in protoplanetary discs

Authors: Roman R. Rafikov, Nicolas P. Cimerman

Abstract: High-resolution sub-mm observations of some protoplanetary discs reveal non-asixymmetric features, which can often be interpreted as dust concentrations in vortices that form at the edges of gaps carved out by the embedded planets. We use recent results on the timescale for the planet-driven vortex development in low-viscosity discs to set constraints on the mass and age of a planet producing the… ▽ More High-resolution sub-mm observations of some protoplanetary discs reveal non-asixymmetric features, which can often be interpreted as dust concentrations in vortices that form at the edges of gaps carved out by the embedded planets. We use recent results on the timescale for the planet-driven vortex development in low-viscosity discs to set constraints on the mass and age of a planet producing the vortex. Knowledge of the age of the central star in a vortex-bearing protoplanetary disc system allows one to set a lower limit on the planetary mass at the level of several tens of $M_\oplus$. Also, an independent upper limit on the planetary mass would constrain the planetary age, although given the current direct imaging detection limits this constraint is not yet very stringent (it is also sensitively dependent on the disc scale height). These results can be extended to account for the history of planetary mass accretion if it is known. We apply our calculations to several protoplanetary discs harbouring vortex-like features as revealed by ALMA and set limits of $(30-50)M_\oplus$ (for disc aspect ratio of $0.1$) on the minimum masses of putative planets that could be responsible for these vortices. Our vortex-based method provides an independent way of constraining the properties of embedded planets, complementary to other approaches. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: 7 pages, 4 figures, MNRAS, in press

Emergence of vortices at the edges of planet-driven gaps in protoplanetary discs

Authors: Nicolas P. Cimerman, Roman R. Rafikov

Abstract: Young planets embedded in protoplanetary discs (PPDs) excite spiral density waves, which propagate, shock and deposit angular momentum in the disc. This results in gap opening around the planetary orbit, even for low (sub-thermal) mass planets, provided that the effective viscosity in the disc is low. The edges of these planet-induced gaps are known to be prone to emergence of observable vortices… ▽ More Young planets embedded in protoplanetary discs (PPDs) excite spiral density waves, which propagate, shock and deposit angular momentum in the disc. This results in gap opening around the planetary orbit, even for low (sub-thermal) mass planets, provided that the effective viscosity in the disc is low. The edges of these planet-induced gaps are known to be prone to emergence of observable vortices via the Rossby Wave Instability (RWI). We study timescales for the development of vortices driven by low mass planets in inviscid discs. We employ a recently developed semi-analytical theory of vortensity production by the planet-driven shock to predict vortensity evolution near the planet, from which we derive the radial profile of the planet-induced gap as a function of time (this procedure can have multiple other uses, e.g. to study dust trapping, suppression of pebble accretion, etc.). We then analyze the linear stability of the gap edges against the RWI, obtaining the timescales for the first appearance of unstable modes and (later) fully developed vortices at gap edges. We present useful formulae for these timescales as functions of planetary and disc parameters and provide their physical justification. We also thoroughly test our semi-analytical framework against high resolution 2D hydrodynamic simulations, confirming the accuracy of our theoretical predictions. We discuss ways in which our semi-analytical framework can be extended to incorporate additional physics, e.g. planetary accretion, migration, and non-zero disc viscosity. Our results can be used to interpret observations of PPDs and to predict emergence of vortices in simulations. △ Less

Submitted 6 December, 2022; originally announced December 2022.

Comments: 21 pages, 17 figures, accepted in MNRAS

Density waves in protoplanetary discs excited by eccentric planets: linear theory

Authors: Callum W. Fairbairn, Roman R. Rafikov

Abstract: Spiral density waves observed in protoplanetary discs have often been used to infer the presence of embedded planets. This inference relies both on simulations as well as the linear theory of planet-disc interaction developed for planets on circular orbits to predict the morphology of the density wake. In this work we develop and implement a linear framework for calculating the structure of the de… ▽ More Spiral density waves observed in protoplanetary discs have often been used to infer the presence of embedded planets. This inference relies both on simulations as well as the linear theory of planet-disc interaction developed for planets on circular orbits to predict the morphology of the density wake. In this work we develop and implement a linear framework for calculating the structure of the density wave in a gaseous disc driven by an eccentric planet. Our approach takes into account both the essential azimuthal and temporal periodicities of the problem, allowing us to treat any periodic perturbing potential (i.e. not only that of an eccentric planet). We test our framework by calculating the morphology of the density waves excited by an eccentric, low-mass planet embedded in a globally isothermal disc and compare our results to the recent direct numerical simulations (and heuristic wavelet analysis) of the same problem by Zhu and Zhang. We find excellent agreement with the numerical simulations, capturing all the complex eccentric features including spiral bifurcations, wave crossings and planet-wave detachments, with improved accuracy and detail compared with the wavelet method. This illustrates the power of our linear framework in reproducing the morphology of complicated time-dependent density wakes, presenting it as a valuable tool for future studies of eccentric planet-disc interactions. △ Less

Submitted 27 September, 2022; v1 submitted 29 July, 2022; originally announced July 2022.

Comments: 11 pages, 6 figures, accepted in MNRAS

Radial Profiles of Surface Density in Debris Discs

Authors: Roman R. Rafikov

Abstract: Resolved observations of debris discs can be used to derive radial profiles of Azimuthally-averaged Surface Density (ASD), which carries important information about the disc structure even in presence of non-axisymmetric features and has improved signal-to-noise characteristics. We develop a (semi-)analytical formalism allowing one to relate ASD to the underlying semi-major axis and eccentricity d… ▽ More Resolved observations of debris discs can be used to derive radial profiles of Azimuthally-averaged Surface Density (ASD), which carries important information about the disc structure even in presence of non-axisymmetric features and has improved signal-to-noise characteristics. We develop a (semi-)analytical formalism allowing one to relate ASD to the underlying semi-major axis and eccentricity distributions of the debris particles in a straightforward manner. This approach does not involve the distribution of particle apsidal angles, thus simplifying calculations. It is a much faster, more flexible and effective way of calculating ASD than the Monte Carlo sampling of orbital parameters of debris particles. We present explicit analytical results based on this technique for a number of particle eccentricity distributions, including two cases of particular practical importance: a prescribed radial profile of eccentricity, and the Rayleigh distribution of eccentricities. We then show how our framework can be applied to observations of debris discs and rings for retrieving either the semi-major axis distribution or (in some cases) the eccentricity distribution of debris, thus providing direct information about the architecture and dynamical processes operating in debris discs. Our approach also provides a fast and efficient way of forward modeling observations. Applications of this technique to other astrophysical systems, e.g. the nuclear stellar disc in M31 or tenuous planetary rings, are also discussed. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: 16 pages, 13 figures, submitted to MNRAS, comments welcome

Anatomy of a slow merger: dissecting secularly-driven inspirals of LIGO/Virgo gravitational wave sources

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: The dozens of compact object mergers detected by LIGO/Virgo raise a key theoretical question: how do initially wide binaries shrink sufficiently quickly that they are able to merge via gravitational wave (GW) radiation within a Hubble time? One promising class of answers involves secular driving of binary eccentricity by some external tidal perturbation. This perturbation can arise due to the pres… ▽ More The dozens of compact object mergers detected by LIGO/Virgo raise a key theoretical question: how do initially wide binaries shrink sufficiently quickly that they are able to merge via gravitational wave (GW) radiation within a Hubble time? One promising class of answers involves secular driving of binary eccentricity by some external tidal perturbation. This perturbation can arise due to the presence of a tertiary point mass, in which case the system exhibits Lidov-Kozai (LK) dynamics, or it can stem from the tidal field of the stellar cluster in which the binary orbits. While these secular tide-driven mechanisms have been studied exhaustively in the case of no GW emission, when GWs are included the dynamical behavior is still incompletely understood. In this paper we consider compact object binaries driven to merger via high eccentricity excitation by (doubly-averaged, test-particle quadrupole level) cluster tides - which includes LK-driven mergers as a special case - and include the effects of both general relativistic precession and GW emission. We provide for the first time an analytical understanding of the different evolutionary stages of the binary's semimajor axis, secular oscillation timescale, and phase space structure all the way to merger. Our results will inform future population synthesis calculations of compact object binary mergers from hierarchical triples and stellar clusters. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 30 pages, 12 figures. To be submitted to ApJ. Comments welcome

Eclipse timing the Milky Way's gravitational potential

Authors: Sukanya Chakrabarti, Daniel J. Stevens, Jason Wright, Roman R. Rafikov, Philip Chang, Thomas Beatty, Daniel Huber

Abstract: We show that a small, but \textit{measurable} shift in the eclipse mid-point time of eclipsing binary (EBs) stars of $\sim$ 0.1 seconds over a decade baseline can be used to directly measure the Galactic acceleration of stars in the Milky Way at $\sim$ kpc distances from the Sun. We consider contributions to the period drift rate from dynamical mechanisms other than the Galaxy's gravitational fiel… ▽ More We show that a small, but \textit{measurable} shift in the eclipse mid-point time of eclipsing binary (EBs) stars of $\sim$ 0.1 seconds over a decade baseline can be used to directly measure the Galactic acceleration of stars in the Milky Way at $\sim$ kpc distances from the Sun. We consider contributions to the period drift rate from dynamical mechanisms other than the Galaxy's gravitational field, and show that the Galactic acceleration can be reliably measured using a sample of $\textit{Kepler}$ EBs with orbital and stellar parameters from the literature. Given the uncertainties on the formulation of tidal decay, our approach here is necessarily approximate, and the contribution from tidal decay is an upper limit assuming the stars are not tidally synchronized. We also use simple analytic relations to search for well-timed sources in the \textit{Kepler} field, and find $\sim$ 70 additional detached EBs with low eccentricities that have estimated timing precision better than 1 second. We illustrate the method with a prototypical, precisely timed EB using an archival \textit{Kepler} light curve and a modern synthetic \textit{HST} light curve (which provides a decade baseline). This novel method establishes a realistic possibility for obtaining fundamental Galactic parameters using eclipse timing to measure Galactic accelerations, along with other emerging new methods, including pulsar timing and extreme precision radial velocity observations. This acceleration signal grows quadratically with time. Therefore, given baselines established in the near-future for distant EBs, we can expect to measure the period drift in the future with space missions like \textit{JWST} and the \textit{Roman Space Telescope}. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 7 pages, 3 figures, submitted to ApJ Letters

Boundary Layers of Accretion Disks: Wave-Driven Transport and Disk Evolution

Authors: Matthew S. B. Coleman, Roman R. Rafikov, Alexander A. Philippov

Abstract: Astrophysical objects possessing a material surface (white dwarfs, young stars, etc.) may accrete gas from the disc through the so-called surface boundary layer (BL), in which the angular velocity of the accreting gas experiences a sharp drop. Acoustic waves excited by the supersonic shear in the BL play an important role in mediating the angular momentum and mass transport through that region. He… ▽ More Astrophysical objects possessing a material surface (white dwarfs, young stars, etc.) may accrete gas from the disc through the so-called surface boundary layer (BL), in which the angular velocity of the accreting gas experiences a sharp drop. Acoustic waves excited by the supersonic shear in the BL play an important role in mediating the angular momentum and mass transport through that region. Here we examine the characteristics of the angular momentum transport produced by the different types of wave modes emerging in the inner disc, using the results of a large suite of hydrodynamic simulations of the BLs. We provide a comparative analysis of the transport properties of different modes across the range of relevant disc parameters. In particular, we identify the types of modes which are responsible for the mass accretion onto the central object. We find the correlated perturbations of surface density and radial velocity to provide an important contribution to the mass accretion rate. Although the wave-driven transport is intrinsically non-local, we do observe a clear correlation between the angular momentum flux injected into the disc by the waves and the mass accretion rate through the BL. We find the efficiency of angular momentum transport (normalized by thermal pressure) to be a weak function of the flow Mach number. We also quantify the wave-driven evolution of the inner disc, in particular the modification of the angular frequency profile in the disc. Our results pave the way for understanding wave-mediated transport in future three-dimensional, magnetohydrodynamic studies of the BLs. △ Less

Submitted 10 March, 2022; v1 submitted 4 November, 2021; originally announced November 2021.

Comments: 16 pages, 9 figures, submitted to MNRAS

Planet-driven density waves in protoplanetary discs: numerical verification of nonlinear evolution theory

Authors: Nicolas P. Cimerman, Roman R. Rafikov

Abstract: Gravitational coupling between protoplanetary discs and planets embedded in them leads to the emergence of spiral density waves, which evolve into shocks as they propagate through the disc. We explore the performance of a semi-analytical framework for describing the nonlinear evolution of the global planet-driven density waves, focusing on the low planet mass regime (below the so-called thermal ma… ▽ More Gravitational coupling between protoplanetary discs and planets embedded in them leads to the emergence of spiral density waves, which evolve into shocks as they propagate through the disc. We explore the performance of a semi-analytical framework for describing the nonlinear evolution of the global planet-driven density waves, focusing on the low planet mass regime (below the so-called thermal mass). We show that this framework accurately captures the (quasi-)self-similar evolution of the wave properties expressed in terms of properly rescaled variables, provided that certain theoretical inputs are calibrated using numerical simulations (an approximate, first principles calculation of the wave evolution based on the inviscid Burgers equation is in qualitative agreement with simulations but overpredicts wave damping at the quantitative level). We provide fitting formulae for such inputs, in particular, the strength and global shape of the planet-driven shock accounting for nonlinear effects. We use this nonlinear framework to theoretically compute vortensity production in the disc by the global spiral shock and numerically verify the accuracy of this calculation. Our results can be used for interpreting observations of spiral features in discs, kinematic signatures of embedded planets in CO line emission ("kinks"), and for understanding the emergence of planet-driven vortices in protoplanetary discs. △ Less

Submitted 13 September, 2021; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: 21 pages, 18 figures, Accepted for publication at MNRAS

Planet formation in stellar binaries: Global simulations of planetesimal growth

Authors: Kedron Silsbee, Roman R. Rafikov

Abstract: Planet formation around one component of a tight, eccentric binary system such as $γ$ Cephei (with semimajor axis around 20 AU) is theoretically challenging because of destructive high-velocity collisions between planetesimals. Despite this fragmentation barrier, planets are known to exist in such orbital configurations. Here we present a novel numerical framework for carrying out multi-annulus co… ▽ More Planet formation around one component of a tight, eccentric binary system such as $γ$ Cephei (with semimajor axis around 20 AU) is theoretically challenging because of destructive high-velocity collisions between planetesimals. Despite this fragmentation barrier, planets are known to exist in such orbital configurations. Here we present a novel numerical framework for carrying out multi-annulus coagulation-fragmentation calculations of planetesimal growth, which fully accounts for the specifics of planetesimal dynamics in binaries, details of planetesimal collision outcomes, and the radial transport of solids in the disk due to the gas drag-driven inspiral. Our dynamical inputs properly incorporate the gravitational effects of both the eccentric stellar companion and the massive non-axisymmetric protoplanetary disk in which planetesimals reside, as well as gas drag. We identify a set of disk parameters that lead to successful planetesimal growth in systems such as $γ$ Cephei or $α$ Centauri starting from $1-10$ km size objects. We identify the apsidal alignment of a protoplanetary disk with the binary orbit as one of the critical conditions for successful planetesimal growth: It naturally leads to the emergence of a dynamically quiet location in the disk (as long as the disk eccentricity is of order several percent), where favorable conditions for planetesimal growth exist. Accounting for the gravitational effect of a protoplanetary disk plays a key role in arriving at this conclusion, in agreement with our previous results. These findings lend support to the streaming instability as the mechanism of planetesimal formation. They provide important insights for theories of planet formation around both binary and single stars, as well as for the hydrodynamic simulations of protoplanetary disks in binaries (for which we identify a set of key diagnostics to verify). △ Less

Submitted 23 July, 2021; originally announced July 2021.

Comments: Accepted to A&A

Journal ref: A&A 652, A104 (2021)

Boundary Layers of Accretion Disks: Discovery of Vortex-Driven Modes and Other Waves

Authors: Matthew S. B. Coleman, Roman R. Rafikov, Alexander A. Philippov

Abstract: Disk accretion onto weakly magnetized objects possessing a material surface must proceed via the so-called boundary layer (BL) - a region at the inner edge of the disk, in which the velocity of accreting material abruptly decreases from its Keplerian value. Supersonic shear arising in the BL is known to be conducive to excitation of acoustic waves that propagate into both the accretor and the disk… ▽ More Disk accretion onto weakly magnetized objects possessing a material surface must proceed via the so-called boundary layer (BL) - a region at the inner edge of the disk, in which the velocity of accreting material abruptly decreases from its Keplerian value. Supersonic shear arising in the BL is known to be conducive to excitation of acoustic waves that propagate into both the accretor and the disk, enabling angular momentum and mass transport across the BL. We carry out a numerical exploration of different wave modes that operate near the BL, focusing on their morphological characteristics in the innermost parts of accretion disk. Using a large suite of simulations covering a broad range of Mach numbers (of the supersonic shear flow in the BL), we provide accurate characterization of the different types of modes, verifying their properties against analytical results, when available. We discover new types of modes, in particular, global spiral density waves launched by vortices forming in the disk near the BL as a result of the Rossby wave instability; this instability is triggered by the vortensity production in that region caused by the nonlinear damping of acoustic waves. Azimuthal wavenumbers of the dominant modes that we observe appear to increase monotonically with the Mach number of the runs, but a particular mix of modes found in a simulation is mildly stochastic. Our results provide a basis for better understanding of the angular momentum and mass transport across the BL as well as the emission variability in accreting objects. △ Less

Submitted 11 October, 2021; v1 submitted 22 March, 2021; originally announced March 2021.

Comments: 25 pages, 16 figures, changed title and reorganized structure, replaced to match the version accepted by MNRAS

Secular dynamics of binaries in stellar clusters -- III. Doubly-averaged dynamics in the presence of general relativistic precession

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: Secular evolution of binaries driven by an external (tidal) potential is a classic astrophysical problem. Tidal perturbations can arise due to an external point mass, as in the Lidov-Kozai (LK) theory of hierarchical triples, or due to an extended stellar system (e.g. galaxy or globular cluster) in which the binary resides. For many applications, general-relativistic (GR) apsidal precession is imp… ▽ More Secular evolution of binaries driven by an external (tidal) potential is a classic astrophysical problem. Tidal perturbations can arise due to an external point mass, as in the Lidov-Kozai (LK) theory of hierarchical triples, or due to an extended stellar system (e.g. galaxy or globular cluster) in which the binary resides. For many applications, general-relativistic (GR) apsidal precession is important, and has been accounted for in some LK calculations. Here we generalise and extend these studies by exploring in detail the effect of GR precession on (quadrupole-level) tidal evolution of binaries orbiting in arbitrary axisymmetric potentials (which includes LK theory as a special case). We study the (doubly-averaged) orbital dynamics for arbitrary strengths of GR and binary initial conditions and uncover entirely new phase space morphologies with important implications for the binary orbital evolution. We also explore how GR precession affects secular evolution of binary orbital elements when the binary reaches high eccentricity ($e\to 1$) and delineate several different dynamical regimes. Our results are applicable to a variety of astrophysical systems. In particular, they can be used to understand the high-eccentricity behaviour of (cluster) tide-driven compact object mergers -- i.e. LIGO/Virgo gravitational wave sources -- for which GR effects are crucial. △ Less

Submitted 19 May, 2021; v1 submitted 30 November, 2020; originally announced November 2020.

Comments: Final version, accepted for publication in MNRAS

Formation of Gaps in Self-Gravitating Debris Disks by Secular Resonance in a Single-planet System I: A Simplified Model

Authors: Antranik A. Sefilian, Roman R. Rafikov, Mark C. Wyatt

Abstract: Spatially resolved images of debris disks frequently reveal complex morphologies such as gaps, spirals, and warps. Most existing models for explaining such morphologies focus on the role of massive perturbers (i.e. planets, stellar companions), ignoring the gravitational effects of the disk itself. Here we investigate the secular interaction between an eccentric planet and a massive, external debr… ▽ More Spatially resolved images of debris disks frequently reveal complex morphologies such as gaps, spirals, and warps. Most existing models for explaining such morphologies focus on the role of massive perturbers (i.e. planets, stellar companions), ignoring the gravitational effects of the disk itself. Here we investigate the secular interaction between an eccentric planet and a massive, external debris disk using a simple analytical model. Our framework accounts for both the gravitational coupling between the disk and the planet, as well as the disk self-gravity -- with the limitation that it ignores the non-axisymmetric component of the disk (self-)gravity. We find generally that even when the disk is less massive than the planet, the system may feature secular resonances within the disk (contrary to what may be naively expected), where planetesimal eccentricities get significantly excited. Given this outcome we propose that double-ringed debris disks, such as those around HD 107146 and HD 92945, could be the result of secular resonances with a yet-undetected planet interior to the disk. We characterize the dependence of the properties of the secular resonances (i.e. locations, timescales, and widths) on the planet and disk parameters, finding that the mechanism is robust provided the disk is massive enough. As an example, we apply our results to HD 107146 and find that this mechanism readily produces $\sim 20$ au wide non-axisymmetric gaps. Our results may be used to set constraints on the total mass of double-ringed debris disks. We demonstrate this for HD 206893, for which we infer a disk mass of $\approx 170$ Earth masses by considering perturbations from the known brown dwarf companion. △ Less

Submitted 8 January, 2021; v1 submitted 29 October, 2020; originally announced October 2020.

Comments: Accepted for publication in ApJ: 25 Pages, 13 Figures, 1 Table, 3 Animations (see Ancillary files)

Doppler Boosting of the S-stars in the Galactic Center

Authors: Roman R. Rafikov

Abstract: Astrometry and spectroscopy of the S-stars in the Galactic Center provide a unique way to probe the properties of the central supermassive black hole, as well as the post-Newtonian effects caused by its gravity, e.g. gravitational redshift and general relativistic precession. It has also been suggested that the photometry of S-stars can be used for studying the properties of gaseous environment of… ▽ More Astrometry and spectroscopy of the S-stars in the Galactic Center provide a unique way to probe the properties of the central supermassive black hole, as well as the post-Newtonian effects caused by its gravity, e.g. gravitational redshift and general relativistic precession. It has also been suggested that the photometry of S-stars can be used for studying the properties of gaseous environment of Sgr A$^*$. Due to the high velocities of the S-stars, sometimes approaching $~0.1c$, their photometric signal should be considerably affected by the Doppler boosting. We calculate this relativistic effect for several S-stars closely approaching the central black hole (most of them recently announced) and show that the amplitude of the photometric variability due to the Doppler boosting for some of them (S62 and S4714) exceeds $6\%$; for the well studied star S2 it is about $2\%$. Measurement of the Doppler boosting can confirm the existence and help refine orbital parameters of the S-stars with noisy spectroscopy and astrometry. This effect should be explicitly accounted for when the photometry of S-stars is used for probing the medium around the Sgr A$^*$. We discuss the observability of the Doppler boosting given the complications typical for the Galactic Center and conclude, in particular, that the purely photometric detection of the higher order relativistic corrections to the Doppler boosting signal (due to the gravitational redshift and transverse Doppler shift, which we also calculate) is hardly possible for the S-stars. △ Less

Submitted 19 January, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 6 pages, 4 figures, published in ApJL

Journal ref: The Astrophysical Journal Letters, Volume 905, Issue 2, id.L35 (2020)

Gaps and Rings in Protoplanetary Disks with Realistic Thermodynamics: The Critical Role of In-Plane Radiation Transport

Authors: Ryan Miranda, Roman R. Rafikov

Abstract: Many protoplanetary disks exhibit annular gaps in dust emission, which may be produced by planets. Simulations of planet-disk interaction aimed at interpreting these observations often treat the disk thermodynamics in an overly simplified manner, which does not properly capture the dynamics of planet-driven density waves driving gap formation. Here we explore substructure formation in disks using… ▽ More Many protoplanetary disks exhibit annular gaps in dust emission, which may be produced by planets. Simulations of planet-disk interaction aimed at interpreting these observations often treat the disk thermodynamics in an overly simplified manner, which does not properly capture the dynamics of planet-driven density waves driving gap formation. Here we explore substructure formation in disks using analytical calculations and hydrodynamical simulations that include a physically-motivated prescription for radiative effects associated with the planet-induced density waves. For the first time, our treatment accounts not only for cooling from the disk surface, but also for radiation transport along the disk midplane. We show that this in-plane cooling, with a characteristic timescale typically an order of magnitude shorter than the one due to surface cooling, plays a critical role in density wave propagation and dissipation (we provide a simple estimate of this timescale). We also show that viscosity, at the levels expected in protoplanetary disks ($α\lesssim 10^{-3}$), has a negligible effect on density wave dynamics. Using synthetic maps of dust continuum emission, we find that the multiplicity and shape of the gaps produced by planets are sensitive to the physical parameters---disk temperature, mass, and opacity---that determine the damping of density waves. Planets orbiting at $\lesssim 20$ au produce the most diverse variety of gap/ring structures, although significant variation is also found for planets at $\gtrsim 50$ au. By improving the treatment of physics governing planet-disk coupling, our results present new ways of probing the planetary interpretation of annular substructures in disks. △ Less

Submitted 5 October, 2020; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: 31 pages, 16 figures, accepted for publication in ApJ

O(N^2) fragmentation algorithm

Authors: Roman R. Rafikov, Kedron Silsbee, Richard A. Booth

Abstract: Collisional fragmentation is a ubiquitous phenomenon arising in a variety of astrophysical systems, from asteroid belts to debris and protoplanetary disks. Numerical studies of fragmentation typically rely on discretizing the size distribution of colliding objects into a large number N of bins in mass space, usually logarithmically spaced. A standard approach for redistributing the debris produced… ▽ More Collisional fragmentation is a ubiquitous phenomenon arising in a variety of astrophysical systems, from asteroid belts to debris and protoplanetary disks. Numerical studies of fragmentation typically rely on discretizing the size distribution of colliding objects into a large number N of bins in mass space, usually logarithmically spaced. A standard approach for redistributing the debris produced in collisions into the corresponding mass bins results in O(N^3) calculation, which leads to significant computational overhead when N is large. Here we formulate a more efficient explicit O(N^2) fragmentation algorithm, which works when the size spectrum of fragments produced in an individual collision has a self-similar shape with only a single characteristic mass scale (which can have arbitrary dependence on the energy and masses of colliding objects). Fragment size spectra used in existing fragmentation codes typically possess this property. We also show that our O(N^2) approach can be easily extended to work with non-self-similar fragment size distributions, for which we provide a worked example. This algorithm offers a substantial speedup of fragmentation calculations for large N> 100, even over the implicit methods, making it an attractive tool for studying collisionally evolving systems. △ Less

Submitted 5 November, 2019; originally announced November 2019.

Comments: 9 pages, 5 figures, submitted to AAS Journals

Planet-disk interaction in disks with cooling: basic theory

Authors: Ryan Miranda, Roman R. Rafikov

Abstract: Gravitational coupling between young planets and their parent disks is often explored using numerical simulations, which typically treat the disk thermodynamics in a highly simplified manner. In particular, many studies adopt the locally isothermal approximation, in which the disk temperature is a fixed function of the stellocentric distance. We explore the dynamics of planet-driven density waves… ▽ More Gravitational coupling between young planets and their parent disks is often explored using numerical simulations, which typically treat the disk thermodynamics in a highly simplified manner. In particular, many studies adopt the locally isothermal approximation, in which the disk temperature is a fixed function of the stellocentric distance. We explore the dynamics of planet-driven density waves in disks with more general thermodynamics, in which the temperature is relaxed towards an equilibrium profile on a finite cooling timescale $t_{\rm c}$. We use both linear perturbation theory and direct numerical simulations to examine the global structure of density waves launched by planets in such disks. A key diagnostic used in this study is the behavior of the wave angular momentum flux (AMF), which directly determines the evolution of the underlying disk. The AMF of free waves is constant for slowly cooling (adiabatic) disks, but scales with the disk temperature for rapidly cooling (and locally isothermal) disks. However, cooling must be extremely fast, with $β= Ωt_{\rm c} \lesssim 10^{-3}$ for the locally isothermal approximation to provide a good description of density wave dynamics in the linear regime (relaxing to $β\lesssim 10^{-2}$ when nonlinear effects are important). For intermediate cooling timescales, density waves are subject to a strong linear damping. This modifies the appearance of planet-driven spiral arms and the characteristics of axisymmetric structures produced by massive planets: in disks with $β\approx 0.1$ -- $1$, a near-thermal mass planet opens only a single wide gap around its orbit, in contrast to the several narrow gaps produced when cooling is either faster or slower. △ Less

Submitted 24 February, 2020; v1 submitted 4 November, 2019; originally announced November 2019.

Comments: 22 pages, 8 figures, accepted for publication in ApJ

Envelopes of embedded super-Earths II. Three-dimensional isothermal simulations

Authors: William Béthune, Roman R. Rafikov

Abstract: Massive planetary cores embedded in protoplanetary discs are believed to accrete extended atmospheres, providing a pathway to forming gas giants and gas-rich super-Earths. The properties of these atmospheres strongly depend on the nature of the coupling between the atmosphere and the surrounding disc. We examine the formation of gaseous envelopes around massive planetary cores via three-dimensiona… ▽ More Massive planetary cores embedded in protoplanetary discs are believed to accrete extended atmospheres, providing a pathway to forming gas giants and gas-rich super-Earths. The properties of these atmospheres strongly depend on the nature of the coupling between the atmosphere and the surrounding disc. We examine the formation of gaseous envelopes around massive planetary cores via three-dimensional inviscid and isothermal hydrodynamic simulations. We focus the changes in the envelope properties as the core mass varies from low (sub-thermal) to high (super-thermal) values, a regime relevant to close-in super-Earths. We show that global envelope properties such as the amount of rotational support or turbulent mixing are mostly sensitive to the ratio of the Bondi radius of the core to its physical size. High-mass cores are fed by supersonic inflows arriving along the polar axis and shocking on the densest parts of the envelope, driving turbulence and mass accretion. Gas flows out of the core's Hill sphere in the equatorial plane, describing a global mass circulation through the envelope. The shell of shocked gas atop the core surface delimits regions of slow (inside) and fast (outside) material recycling by gas from the surrounding disc. While recycling hinders the runaway growth towards gas giants, the inner regions of protoplanetary atmospheres, more immune to mixing, may remain bound to the planet. △ Less

Submitted 8 July, 2019; v1 submitted 5 July, 2019; originally announced July 2019.

Comments: 16 pages, 15 figures, accepted for publication in MNRAS

Envelopes of embedded super-Earths I. Two-dimensional simulations

Authors: William Béthune, Roman R. Rafikov

Abstract: Measurements of exoplanetary masses and radii have revealed a population of massive super-Earths --- planets sufficiently large that, according to one dimensional models, they should have turned into gas giants. To better understand the origin of these objects, we carry out hydrodynamical simulations of planetary cores embedded in a nascent protoplanetary disk. In this first paper of a series, to… ▽ More Measurements of exoplanetary masses and radii have revealed a population of massive super-Earths --- planets sufficiently large that, according to one dimensional models, they should have turned into gas giants. To better understand the origin of these objects, we carry out hydrodynamical simulations of planetary cores embedded in a nascent protoplanetary disk. In this first paper of a series, to gain intuition as well as to develop useful diagnostics, we focus on two-dimensional simulations of the flow around protoplanetary cores. We use the pluto code to study isothermal and adiabatic envelopes around cores of sub- to super-thermal masses, fully resolving the envelope properties down to the core surface. Owing to the conservation of vortensity, envelopes acquire a substantial degree of rotational support when the core mass increases beyond the thermal mass, suggesting a limited applicability of one-dimensional models for describing the envelope structure. The finite size of the core (relatively large for super-Earths) also controls the amount of rotational support in the entire envelope. Steady non-axisymmetric shocks develop in the supersonic envelopes of high-mass cores, triggering mass accretion and turbulent mixing in their interiors. We also examine the influence of the gas self-gravity on the envelope structure. Although it only weakly alters the properties of the envelopes, the gas gravity has significant effect on the properties of the density waves triggered by the core in the protoplanetary disk. △ Less

Submitted 3 July, 2019; originally announced July 2019.

Comments: 17 pages, 24 figures; accepted for publication in MNRAS Published by Oxford University Press

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 487, Issue 2, August 2019, Pages 2319-2334

Compact object binary mergers driven by cluster tides: a new channel for LIGO/Virgo gravitational wave events

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: The detections of gravitational waves produced in mergers of binary black holes (BH) and neutron stars (NS) by LIGO/Virgo have stimulated interest in the origin of the progenitor binaries. Dense stellar systems - globular and nuclear star clusters - are natural sites of compact object binary formation and evolution towards merger. Here we explore a new channel for the production of binary mergers… ▽ More The detections of gravitational waves produced in mergers of binary black holes (BH) and neutron stars (NS) by LIGO/Virgo have stimulated interest in the origin of the progenitor binaries. Dense stellar systems - globular and nuclear star clusters - are natural sites of compact object binary formation and evolution towards merger. Here we explore a new channel for the production of binary mergers in clusters, in which the tidal field of the cluster secularly drives the binary to high eccentricity (even in the absence of a central massive black hole) until gravitational wave emission becomes important. We employ the recently developed secular theory of cluster tide-driven binary evolution to compute present day merger rates for BH-BH, NS-BH and NS-NS binaries, varying cluster potential and central concentration of the binary population (but ignoring cluster evolution and stellar flybys for now). Unlike other mechanisms, this new dynamical channel can produce a significant number of mergers out to cluster-centric distances of several pc. For NS-NS binaries we find merger rates in the range $0.01-0.07$ Gpc$^{-3}$ yr$^{-1}$ from globular clusters and $0.1-0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters. For NS-BH and BH-BH binaries we find small merger rates from globular clusters, but a rate of $0.1 - 0.2$ Gpc$^{-3}$ yr$^{-1}$ from cusped nuclear clusters, contributing to the observed LIGO/Virgo rate at the level of several per cent. Therefore, cluster tide-driven mergers constitute a new channel that can be further explored with current and future gravitational wave detectors. △ Less

Submitted 16 August, 2019; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: Final version, accepted for publication in ApJ Letters (9 pages)

On the planetary interpretation of multiple gaps and rings in protoplanetary disks seen by ALMA

Authors: Ryan Miranda, Roman R. Rafikov

Abstract: It has been recently suggested that the multiple concentric rings and gaps discovered by ALMA in many protoplanetary disks may be produced by a single planet, as a result of the complex propagation and dissipation of the multiple spiral density waves it excites in the disk. Numerical efforts to verify this idea have largely utilized the so-called locally isothermal approximation with a prescribed… ▽ More It has been recently suggested that the multiple concentric rings and gaps discovered by ALMA in many protoplanetary disks may be produced by a single planet, as a result of the complex propagation and dissipation of the multiple spiral density waves it excites in the disk. Numerical efforts to verify this idea have largely utilized the so-called locally isothermal approximation with a prescribed disk temperature profile. However, in protoplanetary disks this approximation does not provide an accurate description of the density wave dynamics on scales of tens of au. Moreover, we show that locally isothermal simulations tend to overestimate the contrast of ring and gap features, as well as misrepresent their positions, when compared to simulations in which the energy equation is evolved explicitly. This outcome is caused by the non-conservation of the angular momentum flux of linear perturbations in locally isothermal disks. We demonstrate this effect using simulations of locally isothermal and adiabatic disks (with essentially identical temperature profiles) and show how the dust distributions, probed by mm wavelength observations, differ between the two cases. Locally isothermal simulations may thus underestimate the masses of planets responsible for the formation of multiple gaps and rings on scales of tens of au observed by ALMA. We suggest that caution should be exercised in using the locally isothermal simulations to explore planet-disk interaction, as well as in other studies of wave-like phenomena in astrophysical disks. △ Less

Submitted 20 May, 2019; originally announced May 2019.

Comments: 8 pages, 3 figures, accepted for publication in ApJ Letters

Potential softening and eccentricity dynamics in razor-thin, nearly-Keplerian discs

Authors: Antranik A. Sefilian, Roman R. Rafikov

Abstract: In many astrophysical problems involving discs (gaseous or particulate) orbiting a dominant central mass, gravitational potential of the disc plays an important dynamical role. Its impact on the motion of external objects, as well as on the dynamics of the disc itself, can usually be studied using secular approximation. This is often done using softened gravity to avoid singularities arising in ca… ▽ More In many astrophysical problems involving discs (gaseous or particulate) orbiting a dominant central mass, gravitational potential of the disc plays an important dynamical role. Its impact on the motion of external objects, as well as on the dynamics of the disc itself, can usually be studied using secular approximation. This is often done using softened gravity to avoid singularities arising in calculation of the orbit-averaged potential --- disturbing function --- of a razor-thin disc using classical Laplace-Lagrange theory. We explore the performance of several softening formalisms proposed in the literature in reproducing the correct eccentricity dynamics in the disc potential. We identify softening models that, in the limit of zero softening, give results converging to the expected behavior exactly, approximately or not converging at all. We also develop a general framework for computing secular disturbing function given an arbitrary softening prescription for a rather general form of the interaction potential. Our results demonstrate that numerical treatments of the secular disc dynamics, representing the disc as a collection of $N$ gravitationally interacting annuli, are rather demanding: for a given value of the (dimensionless) softening parameter, $ς\ll 1$, accurate representation of eccentricity dynamics requires $N \sim Cς^{-χ}\gg 1$, with $C\sim O(10)$, $1.5\lesssim χ\lesssim 2$. In discs with sharp edges a very small value of the softening parameter $ς$ ($\lesssim 10^{-3}$) is required to correctly reproduce eccentricity dynamics near the disc boundaries; this finding is relevant for modelling planetary rings. △ Less

Submitted 25 August, 2019; v1 submitted 16 April, 2019; originally announced April 2019.

Comments: Final version accepted for publication in MNRAS: 23 Pages, 9 Figures, 2 Tables

The Gemini Planet Imager Exoplanet Survey: Giant Planet and Brown Dwarf Demographics From 10-100 AU

Authors: Eric L. Nielsen, Robert J. De Rosa, Bruce Macintosh, Jason J. Wang, Jean-Baptiste Ruffio, Eugene Chiang, Mark S. Marley, Didier Saumon, Dmitry Savransky, S. Mark Ammons, Vanessa P. Bailey, Travis Barman, Celia Blain, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, Ian Czekala, Rene Doyon, Gaspard Duchene, Thomas M. Esposito, Daniel Fabrycky, Michael P. Fitzgerald, Katherine B. Follette, Jonathan J. Fortney, Benjamin L. Gerard , et al. (40 additional authors not shown)

Abstract: We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation be… ▽ More We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^αa^β$, finding $α$ = -2.4 $\pm$ 0.8 and $β$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability. △ Less

Submitted 10 April, 2019; originally announced April 2019.

Comments: 52 pages, 18 figures. AJ in press

Multi-Messenger Astrophysics with Pulsar Timing Arrays

Authors: Luke Zoltan Kelley, Maria Charisi, Sarah Burke-Spolaor, Joseph Simon, Laura Blecha, Tamara Bogdanovic, Monica Colpi, Julie Comerford, Daniel J. D'Orazio, Massimo Dotti, Michael Eracleous, Matthew Graham, Jenny E. Greene, Zoltán Haiman, Kelly Holley-Bockelmann, Erin Kara, Bernard Kelly, S. Komossa, Shane L. Larson, Xin Liu, Chung-Pei Ma, Scott Noble, Vasileios Paschalidis, Roman R. Rafikov, Vikram Ravi , et al. (7 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will… ▽ More Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 decadal review. This is one of five core white papers written by members of the NANOGrav Collaboration, see also: J.Cordes et al., S.R.Taylor et al., X.Siemens et al., and E.Fonseca et al

Secular dynamics of binaries in stellar clusters II: dynamical evolution

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: Dense stellar clusters are natural sites for the origin and evolution of exotic objects such as relativistic binaries (potential gravitational wave sources), blue stragglers, etc. We investigate the secular dynamics of a binary system driven by the global tidal field of an axisymmetric stellar cluster in which the binary orbits. In a companion paper (Hamilton & Rafikov 2019a) we developed a genera… ▽ More Dense stellar clusters are natural sites for the origin and evolution of exotic objects such as relativistic binaries (potential gravitational wave sources), blue stragglers, etc. We investigate the secular dynamics of a binary system driven by the global tidal field of an axisymmetric stellar cluster in which the binary orbits. In a companion paper (Hamilton & Rafikov 2019a) we developed a general Hamiltonian framework describing such systems. The effective (doubly-averaged) Hamiltonian derived there encapsulates all information about the tidal potential experienced by the binary in its orbit around the cluster in a single parameter $Γ$. Here we provide a thorough exploration of the phase-space of the corresponding secular problem as $Γ$ is varied. We find that for $Γ> 1/5$ the phase-space structure and the evolution of binary orbital element are qualitatively similar to the Lidov-Kozai problem. However, this is only one of four possible regimes, because the dynamics are qualitatively changed by bifurcations at $Γ= 1/5,0,-1/5$. We show how the dynamics are altered in each regime and calculate characteristics such as secular evolution timescale, maximum possible eccentricity, etc. We verify the predictions of our doubly-averaged formalism numerically and find it to be very accurate when its underlying assumptions are fulfilled, typically meaning that the secular timescale should exceed the period of the binary around the cluster by $\gtrsim 10-10^2$ (depending on the cluster potential and binary orbit). Our results may be relevant for understanding the nature of a variety of exotic systems harboured by stellar clusters. △ Less

Submitted 16 August, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 25 pages, 14 figures. Final version accepted for publication in MNRAS

Secular dynamics of binaries in stellar clusters I: general formulation and dependence on cluster potential

Authors: Chris Hamilton, Roman R. Rafikov

Abstract: Orbital evolution of binary systems in dense stellar clusters is important in a variety of contexts: origin of blue stragglers, progenitors of compact object mergers, millisecond pulsars, and so on. Here we consider the general problem of secular evolution of the orbital elements of a binary system driven by the smooth tidal field of an axisymmetric stellar cluster (globular, nuclear, etc.) in whi… ▽ More Orbital evolution of binary systems in dense stellar clusters is important in a variety of contexts: origin of blue stragglers, progenitors of compact object mergers, millisecond pulsars, and so on. Here we consider the general problem of secular evolution of the orbital elements of a binary system driven by the smooth tidal field of an axisymmetric stellar cluster (globular, nuclear, etc.) in which the binary orbits. We derive a secular Hamiltonian (averaged over both the inner Keplerian orbit of the binary and its outer orbit within the cluster) valid to quadrupole order for an arbitrary cluster potential and explore its characteristics. This doubly-averaged 'tidal' Hamiltonian depends on just two parameters, which fully absorb the information about the background cluster potential and the binary's orbit within it: a dimensional parameter $A$ setting the secular timescale, and a dimensionless parameter $Γ$ which determines the phase portrait of the binary's inner orbital evolution. We examine the dependence of $A$ and $Γ$ on cluster potential (both spherical and axisymmetric) and on the binary orbit within the cluster. Our theory reproduces known secular results - such as Lidov-Kozai evolution and the effect of the Galactic tide on Oort Cloud comets - in appropriate limits, but is more general. It provides a universal framework for understanding dynamical evolution of various types of binaries driven by the smooth tidal field of any axisymmetric potential. In a companion paper (Hamilton & Rafikov 2019b) we provide a detailed exploration of the resulting orbital dynamics. △ Less

Submitted 16 August, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 24 pages, 13 figures. Final version accepted for publication in MNRAS

Multiple Spiral Arms in Protoplanetary Disks: Linear Theory

Authors: Ryan Miranda, Roman R. Rafikov

Abstract: Recent observations of protoplanetary disks, as well as simulations of planet-disk interaction, have suggested that a single planet may excite multiple spiral arms in the disk, in contrast to the previous expectations based on linear theory (predicting a one-armed density wave). We re-assess the origin of multiple arms in the framework of linear theory, by solving for the global two-dimensional re… ▽ More Recent observations of protoplanetary disks, as well as simulations of planet-disk interaction, have suggested that a single planet may excite multiple spiral arms in the disk, in contrast to the previous expectations based on linear theory (predicting a one-armed density wave). We re-assess the origin of multiple arms in the framework of linear theory, by solving for the global two-dimensional response of a non-barotropic disk to an orbiting planet. We show that the formation of a secondary arm in the inner disk, at about half of the orbital radius of the planet, is a robust prediction of linear theory. This arm becomes stronger than the primary spiral at several tenths of the orbital radius of the planet. Several additional, weaker spiral arms may also form in the inner disk. On the contrary, a secondary spiral arm is unlikely to form in the outer disk. Our linear calculations, fully accounting for the global behavior of both the phases and amplitudes of perturbations, generally support the recently proposed WKB phase argument for the secondary arm origin (as caused by the intricacy of constructive interference of azimuthal harmonics of the perturbation at different radii). We provide analytical arguments showing that the process of a single spiral wake splitting up into multiple arms is a generic linear outcome of wave propagation in differentially rotating disks. It is not unique to planet-driven waves and occurs also in linear calculations of spiral wakes freely propagating with no external torques. These results are relevant for understanding formation of multiple rings and gaps in protoplanetary disks. △ Less

Submitted 17 March, 2019; v1 submitted 23 November, 2018; originally announced November 2018.

Comments: 22 pages, 16 figures, accepted for publication in ApJ

Spin Evolution and Cometary Interpretation of the Interstellar Minor Object 1I/2017 'Oumuamua

Authors: Roman R. Rafikov

Abstract: Observations of the first interstellar minor object 1I/2017 'Oumuamua did not reveal direct signs of outgassing that would have been natural if it had volatile-rich composition. However, a recent measurement by Micheli et al (2018) of a substantial non-gravitational acceleration affecting the orbit of this object has been interpreted as resulting from its cometary activity, which must be rather vi… ▽ More Observations of the first interstellar minor object 1I/2017 'Oumuamua did not reveal direct signs of outgassing that would have been natural if it had volatile-rich composition. However, a recent measurement by Micheli et al (2018) of a substantial non-gravitational acceleration affecting the orbit of this object has been interpreted as resulting from its cometary activity, which must be rather vigorous. Here we critically re-assess this interpretation by exploring the implications of measured non-gravitational acceleration for the 'Oumuamua's rotational state. We show that outgassing torques should drive rapid evolution of 'Oumuamua's spin (on a timescale of a few days), assuming torque asymmetry typical for the Solar System comets. However, given the highly elongated shape of the object, its torque asymmetry is likely higher, implying even faster evolution. This would have resulted in rapid rotational fission of 'Oumuamua during its journey through the Solar System and is clearly incompatible with the relative stability of its rotational state inferred from photometric variability. Based on these arguments, as well as the lack of direct signs of outgassing, we conclude that the classification of 'Oumuamua as a comet (invoked to explain its claimed anomalous acceleration) is questionable. △ Less

Submitted 5 November, 2018; v1 submitted 17 September, 2018; originally announced September 2018.

Comments: Replaced to match a version published in The Astrophysical Journal Letters, 2018, Volume 867, Issue 1, article id. L17

Non-Gravitational Forces and Spin Evolution of Comets

Authors: Roman R. Rafikov

Abstract: Motion of many comets is affected by non-gravitational forces caused by outgassing from their surfaces. Outgassing also produces reactive torques resulting in cometary spin evolution. We propose that the two processes are correlated and show that the change of cometary spin rate over its heliocentric orbit scales linearly with the amplitude of its non-gravitational acceleration. The proportionalit… ▽ More Motion of many comets is affected by non-gravitational forces caused by outgassing from their surfaces. Outgassing also produces reactive torques resulting in cometary spin evolution. We propose that the two processes are correlated and show that the change of cometary spin rate over its heliocentric orbit scales linearly with the amplitude of its non-gravitational acceleration. The proportionality constant depends on the comet size and orbital elements (semi-major axis and eccentricity) and on the (dimensionless) lever arm parameter $ζ$ that relates the outgassing-induced torque and acceleration. We determine $ζ$ for 7 comets for which both non-gravitational acceleration and change of spin period $ΔP$ were measured and verify this relation. This sample spanning almost 4 decades in $ΔP$ yields $\logζ=-2.21\pm 0.54$, surprisingly small value and spread. We then apply our framework to 209 comets with measured non-gravitational accelerations and determine the objects most likely to exhibit large spin period changes, $ΔP\gtrsim 20$ min per orbit assuming rotation period of 10 hr and $ζ$ comparable to our control sample. These objects should be primary targets for future studies of cometary spin variability, further constraining distribution of $ζ$. Lack of comets with very high expected spin rate changes (which is not equivalent to having the highest non-gravitational acceleration) suggests that (1) cometary fission due to outgassing-driven spin-up must be an important process and (2) the distribution of $ζ$ has a lower limit $\sim 10^{-3}$. △ Less

Submitted 13 September, 2018; originally announced September 2018.

Comments: 10 pages, 4 figures, submitted to AAS journals, comments welcome

Secular Evolution Driven by Massive Eccentric Disks/Rings: An Apsidally Aligned Case

Authors: Irina Davydenkova, Roman R. Rafikov

Abstract: Massive eccentric disks (gaseous or particulate) orbiting a dominant central mass appear in many astrophysical systems, including planetary rings, protoplanetary and accretion disks in binaries, and nuclear stellar disks around supermassive black holes in galactic centers. We present an analytical framework for treating the nearly Keplerian secular dynamics of test particles driven by the gravity… ▽ More Massive eccentric disks (gaseous or particulate) orbiting a dominant central mass appear in many astrophysical systems, including planetary rings, protoplanetary and accretion disks in binaries, and nuclear stellar disks around supermassive black holes in galactic centers. We present an analytical framework for treating the nearly Keplerian secular dynamics of test particles driven by the gravity of an eccentric, apsidally aligned, zero-thickness disk with arbitrary surface density and eccentricity profiles. We derive a disturbing function describing the secular evolution of coplanar objects, which is explicitly related (via one-dimensional, convergent integrals) to the disk surface density and eccentricity profiles without using any ad hoc softening of the potential. Our analytical framework is verified via direct orbit integrations, which show it to be accurate in the low-eccentricity limit for a variety of disk models (for disk eccentricity < 0.1-0.2). We find that free precession in the potential of a disk with a smooth surface density distribution can naturally change from prograde to retrograde within the disk. Sharp disk features - edges and gaps - are the locations where this tendency is naturally enhanced, while the precession becomes very fast. Radii where free precession changes sign are the locations where substantial (formally singular) growth of the forced eccentricity of the orbiting objects occurs. Based on our results, we formulate a self-consistent analytical framework for computing an eccentricity profile for an aligned, eccentric disk (with a prescribed surface density profile) capable of precessing as a solid body under its own self-gravity. △ Less

Submitted 7 September, 2018; originally announced September 2018.

Comments: 22 pages, 17 figures, published in ApJ, 864, 74 (2018)

Fast and Slow Precession of Gaseous Debris Disks Around Planet-Accreting White Dwarfs

Authors: Ryan Miranda, Roman R. Rafikov

Abstract: Spectroscopic observations of some metal-rich white dwarfs (WDs), believed to be polluted by planetary material, reveal the presence of compact gaseous metallic disks orbiting them. The observed variability of asymmetric, double-peaked emission line profiles in about half of such systems could be interpreted as the signature of precession of an eccentric gaseous debris disk. The variability timesc… ▽ More Spectroscopic observations of some metal-rich white dwarfs (WDs), believed to be polluted by planetary material, reveal the presence of compact gaseous metallic disks orbiting them. The observed variability of asymmetric, double-peaked emission line profiles in about half of such systems could be interpreted as the signature of precession of an eccentric gaseous debris disk. The variability timescales --- from decades down to $1.4$ yr (recently inferred for the debris disk around HE 1349--2305) --- are in rough agreement with the rate of general relativistic (GR) precession in the test particle limit. However, it has not been demonstrated that this mechanism can drive such a fast, coherent precession of a radially extended (out to $1 R_\odot$) gaseous disk mediated by internal stresses (pressure). Here we use the linear theory of eccentricity evolution in hydrodynamic disks to determine several key properties of eccentric modes in gaseous debris disks around WDs. We find a critical dependence of both the precession period and radial eccentricity distribution of the modes on the inner disk radius, $r_\mathrm{in}$. For small inner radii, $r_\mathrm{in} \lesssim (0.2 - 0.4) R_\odot$, the modes are GR-driven, with periods of $\approx 1 - 10$ yr. For $r_\mathrm{in} \gtrsim (0.2 - 0.4) R_\odot$, the modes are pressure-dominated, with periods of $\approx 3 - 20$ yr. Correspondence between the variability periods and inferred inner radii of the observed disks is in general agreement with this trend. In particular, the short period of HE 1349--2305 is consistent with its small $r_\mathrm{in}$. Circum-WD debris disks may thus serve as natural laboratories for studying the evolution of eccentric gaseous disks. △ Less

Submitted 4 April, 2018; v1 submitted 1 February, 2018; originally announced February 2018.

Comments: 7 pages, 4 figures, 1 table; data points for observed systems added to figures for comparison with theoretical results; accepted for publication in ApJ

1I/2017 'Oumuamua-like Interstellar Asteroids as Possible Messengers from the Dead Stars

Authors: Roman R. Rafikov

Abstract: Discovery of the first interstellar asteroid (ISA) - 1I/2017 'Oumuamua - raised a number of questions regarding its origin. Many of them relate to its lack of cometary activity, suggesting refractory composition of 'Oumuamua. Here we explore the possibility that 'Oumuamua-like ISAs are produced in tidal disruption events (TDEs) of refractory planetoids (asteroids, dwarf planets, etc.) by the white… ▽ More Discovery of the first interstellar asteroid (ISA) - 1I/2017 'Oumuamua - raised a number of questions regarding its origin. Many of them relate to its lack of cometary activity, suggesting refractory composition of 'Oumuamua. Here we explore the possibility that 'Oumuamua-like ISAs are produced in tidal disruption events (TDEs) of refractory planetoids (asteroids, dwarf planets, etc.) by the white dwarfs (WDs). This idea is supported by existing spectroscopic observations of metal-polluted WDs, hinting at predominantly volatile-poor composition of accreted material. We show that such TDEs sourced by realistic planetary systems (including a population of >1000 km planetoids and massive perturbers - Neptune-to-Saturn mass planets) can eject to interstellar space up to 30% of planetary mass involved in them. Collisional fragmentation, caused by convergent vertical motion of the disrupted planetoid's debris inside the Roche sphere of the WD, channels most of the original mass into 0.1-1 km fragments, similar to 'Oumuamua. Such size spectrum of ISAs (very different from the top-heavy distributions expected in other scenarios) implies that planetary TDEs can account for a significant fraction (up to ~30% under optimistic assumptions) of the ISAs. This figure is based on existing observations of WD metal pollution and accounts for observational biases by using realistic models of circum-WD planetary systems. ISAs should exhibit kinematic characteristics similar to old, dynamically hot Galactic populations; we interpret 'Oumuamua's slow Galactic motion as a statistical fluctuation. ISA ejection in individual planetary TDEs is highly anisotropic, resulting in large fluctuations of their space density. We also show that other ISA production mechanisms involving stellar remnants - direct ejection by massive planets around WDs and SN explosions - have difficulty explaining 'Oumuamua-like ISAs. △ Less

Submitted 16 May, 2018; v1 submitted 8 January, 2018; originally announced January 2018.

Comments: 19 pages, 7 figures, accepted to ApJ, comments welcome. Replaced to match the accepted version

Disk Accretion Driven by Spiral Shocks

Authors: Lev Arzamasskiy, Roman R. Rafikov

Abstract: Spiral density waves are known to exist in many astrophysical disks, potentially affecting disk structure and evolution. We conduct a numerical study of the effects produced by a density wave, evolving into a shock, on the characteristics of the underlying disk. We measure the deposition of angular momentum in the disk by spiral shocks of different strength and verify the analytical prediction of… ▽ More Spiral density waves are known to exist in many astrophysical disks, potentially affecting disk structure and evolution. We conduct a numerical study of the effects produced by a density wave, evolving into a shock, on the characteristics of the underlying disk. We measure the deposition of angular momentum in the disk by spiral shocks of different strength and verify the analytical prediction of Rafikov (2016) for the behavior of this quantity, using shock amplitude (which is potentially observable) as the input variable. Good agreement between the theory and numerics is found as we vary shock amplitude (including highly nonlinear shocks), disk aspect ratio, equation of state, radial profiles of the background density and temperature, and pattern speed of the wave. We show that high numerical resolution is required to properly capture shock-driven transport, especially at low wave amplitudes. We also demonstrate that relating local mass accretion rate to shock dissipation in rapidly evolving disks requires accounting for the time-dependent contribution to the angular momentum budget, caused by the time dependence of the radial pressure support. We provide a simple analytical prescription for the behavior of this contribution and demonstrate its excellent agreement with the simulation results. Using these findings we formulate a theoretical framework for studying one-dimensional (in radius) evolution of the shock-mediated accretion disks, which can be applied to a variety of astrophysical systems. △ Less

Submitted 5 February, 2018; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: 16 pages, 13 figures, submitted to ApJ

Evidence that the Directly-Imaged Planet HD 131399 Ab is a Background Star

Authors: Eric L. Nielsen, Robert J. De Rosa, Julien Rameau, Jason J. Wang, Thomas M. Esposito, Maxwell A. Millar-Blanchaer, Christian Marois, Arthur Vigan, S. Mark Ammons, Etienne Artigau, Vanessa P. Bailey, Sarah Blunt, Joanna Bulger, Jeffrey Chilcote, Tara Cotten, René Doyon, Gaspard Duchêne, Daniel Fabrycky, Michael P. Fitzgerald, Katherine B. Follette, Benjamin L. Gerard, Stephen J. Goodsell, James R. Graham, Alexandra Z. Greenbaum, Pascale Hibon , et al. (33 additional authors not shown)

Abstract: We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength cove… ▽ More We present evidence that the recently discovered, directly-imaged planet HD 131399 Ab is a background star with non-zero proper motion. From new JHK1L' photometry and spectroscopy obtained with the Gemini Planet Imager, VLT/SPHERE, and Keck/NIRC2, and a reanalysis of the discovery data obtained with VLT/SPHERE, we derive colors, spectra, and astrometry for HD 131399 Ab. The broader wavelength coverage and higher data quality allow us to re-investigate its status. Its near-infrared spectral energy distribution excludes spectral types later than L0 and is consistent with a K or M dwarf, which are the most likely candidates for a background object in this direction at the apparent magnitude observed. If it were a physically associated object, the projected velocity of HD 131399 Ab would exceed escape velocity given the mass and distance to HD 131399 A. We show that HD 131399 Ab is also not following the expected track for a stationary background star at infinite distance. Solving for the proper motion and parallax required to explain the relative motion of HD 131399 Ab, we find a proper motion of 12.3 mas/yr. When compared to predicted background objects drawn from a galactic model, we find this proper motion to be high, but consistent with the top 4% fastest-moving background stars. From our analysis we conclude that HD 131399 Ab is a background K or M dwarf. △ Less

Submitted 8 November, 2017; v1 submitted 18 May, 2017; originally announced May 2017.

Comments: 35 pages, 20 figures, Published in AJ

Journal ref: Nielsen et al. 2017 AJ, 154, 6

Protoplanetary Disks as (Possibly) Viscous Disks

Authors: Roman R. Rafikov

Abstract: Protoplanetary disks are believed to evolve on Myr timescales in a diffusive (viscous) manner as a result of angular momentum transport driven by internal stresses. Here we use a sample of 26 protoplanetary disks resolved by ALMA with measured (dust-based) masses and stellar accretion rates to derive the dimensionless $α$-viscosity values for individual objects, with the goal of constraining the a… ▽ More Protoplanetary disks are believed to evolve on Myr timescales in a diffusive (viscous) manner as a result of angular momentum transport driven by internal stresses. Here we use a sample of 26 protoplanetary disks resolved by ALMA with measured (dust-based) masses and stellar accretion rates to derive the dimensionless $α$-viscosity values for individual objects, with the goal of constraining the angular momentum transport mechanism. We find that the inferred values of $α$ do not cluster around a single value, but instead have a broad distribution extending from $10^{-4}$ to $0.04$. Moreover, they correlate with neither the global disk parameters (mass, size, surface density) nor the stellar characteristics (mass, luminosity, radius). However, we do find a strong linear correlation between $α$ and the central mass accretion rate $\dot M$. This correlation is unlikely to result from the direct physical effect of $\dot M$ on disk viscosity on global scales. Instead, we suggest that it is caused by the decoupling of stellar $\dot M$ from the global disk characteristics in one of the following ways. (1) The behavior (and range) of $α$ is controlled by a yet unidentified parameter (e.g. ionization fraction, magnetic field strength, or geometry), ultimately driving the variation of $\dot M$. (2) The central $\dot M$ is decoupled from the global viscous mass accretion rate as a result of an instability or mass accumulation (or loss) in the inner disk. (3) Perhaps the most intriguing possibility is that angular momentum in protoplanetary disks is transported non-viscously, e.g. via magnetohydrodynamic winds or spiral density waves. △ Less

Submitted 20 January, 2017; v1 submitted 9 January, 2017; originally announced January 2017.

Comments: 14 pages, 12 figures, submitted to ApJ

Radial Transport and Meridional Circulation in Accretion Disks

Authors: Alexander A. Philippov, Roman R. Rafikov

Abstract: Radial transport of particles, elements and fluid driven by internal stresses in three-dimensional (3D) astrophysical accretion disks is an important phenomenon, potentially relevant for the outward dust transport in protoplanetary disks, origin of the refractory particles in comets, isotopic equilibration in the Earth-Moon system, etc. To gain better insight into these processes, we explore the d… ▽ More Radial transport of particles, elements and fluid driven by internal stresses in three-dimensional (3D) astrophysical accretion disks is an important phenomenon, potentially relevant for the outward dust transport in protoplanetary disks, origin of the refractory particles in comets, isotopic equilibration in the Earth-Moon system, etc. To gain better insight into these processes, we explore the dependence of meridional circulation in 3D disks with shear viscosity on their thermal stratification, and demonstrate strong effect of the latter on the radial flow. Previous locally isothermal studies have normally found a pattern of the radial outflow near the midplane, switching to inflow higher up. Here we show, both analytically and numerically, that a flow, which is inward at all altitudes, is possible in disks with entropy and temperature steeply increasing with height. Such thermodynamic conditions may be typical in the optically thin, viscously heated accretion disks. Disks in which these conditions do not hold should feature radial outflow near the midplane, as long as their internal stress is provided by the shear viscosity. Our results can also be used for designing hydrodynamical disk simulations with a prescribed pattern of the meridional circulation. △ Less

Submitted 7 February, 2017; v1 submitted 7 January, 2017; originally announced January 2017.

Comments: 11 pages, 5 figures, under review in ApJ, revised version after the referee report

Searching for Binary Supermassive Black Holes via Variable Broad Emission Line Shifts: Low Binary Fraction

Authors: Lile Wang, Jenny E. Greene, Wenhua Ju, Roman R. Rafikov, John J. Ruan, Donald P. Schneider

Abstract: Supermassive black hole binaries (SMBHs) are expected to result from galaxy mergers, and thus are natural byproducts (and probes) of hierarchical structure formation in the Universe. They are also the primary expected source of low-frequency gravitational wave emission. We search for binary BHs using time-variable velocity shifts in broad Mg II emission lines of quasars with multi-epoch observatio… ▽ More Supermassive black hole binaries (SMBHs) are expected to result from galaxy mergers, and thus are natural byproducts (and probes) of hierarchical structure formation in the Universe. They are also the primary expected source of low-frequency gravitational wave emission. We search for binary BHs using time-variable velocity shifts in broad Mg II emission lines of quasars with multi-epoch observations. First, we inspect velocity shifts of the binary SMBH candidates identified in Ju et al. (2013), using SDSS spectra with an additional epoch of data that lengthens the typical baseline to ~10 yr. We find variations in the line-of-sight velocity shifts over 10 years that are comparable to the shifts observed over 1-2 years, ruling out the binary model for the bulk of our candidates. We then analyze 1438 objects with 8 yr median time baselines, from which we would expect to see velocity shifts >1000 km/s from sub-pc binaries. We find only one object with an outlying velocity of 448 km/s, indicating, based on our modeling, that ~< 1 per cent (the value varies with different assumptions) of SMBHs that are active as quasars reside in binaries with ~0.1 pc separations. Binaries either sweep through these small separations rapidly or stall at larger radii. △ Less

Submitted 31 October, 2016; originally announced November 2016.

Comments: 13 pages, 7 figures; submitted to ApJ

Discovery of a Substellar Companion to the Nearby Debris Disk Host HR 2562

Authors: Quinn M. Konopacky, Julien Rameau, Gaspard Duchene, Joseph C. Filippazzo, Paige A. Giorla Godfrey, Christian Marois, Eric L. Nielsen, Laurent Pueyo, Roman R. Rafikov, Emily L. Rice, Jason J. Wang, S. Mark Ammons, Vanessa P. Bailey, Travis S. Barman, Joanna Bulger, Sebastian Bruzzone, Jeffrey K. Chilcote, Tara Cotten, Rebekah I. Dawson, Robert J. De Rosa, Rene Doyon, Thomas M. Esposito, Michael P. Fitzgerald, Katherine B. Follette, Stephen Goodsell , et al. (32 additional authors not shown)

Abstract: We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observation… ▽ More We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5σ$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk. △ Less

Submitted 23 August, 2016; originally announced August 2016.

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

On the Eccentricity Excitation in Post-Main Sequence Binaries

Authors: Roman R. Rafikov

Abstract: Several classes of stellar binaries with post-main sequence (post-MS) components - millisecond pulsars with the white dwarf companions (MSP+WD) and periods P_b~30 d, binaries hosting post-asymptotic giant branch (post-AGB) stars or barium stars with P_b~several yr - feature high eccentricities (up to 0.4) despite the expectation of their efficient tidal circularization during the AGB phase. It was… ▽ More Several classes of stellar binaries with post-main sequence (post-MS) components - millisecond pulsars with the white dwarf companions (MSP+WD) and periods P_b~30 d, binaries hosting post-asymptotic giant branch (post-AGB) stars or barium stars with P_b~several yr - feature high eccentricities (up to 0.4) despite the expectation of their efficient tidal circularization during the AGB phase. It was suggested that the eccentricities of these binaries can be naturally excited by their tidal coupling to the circumbinary disk, formed by the material ejected from the binary. Here we critically reassess this idea using simple arguments rooted in the global angular momentum conservation of the disk+binary system. Compared to previous studies we (1) fully account for the viscous spreading of the circumbinary disk, (2) consider the possibility of reaccretion from the disk onto the binary (in agreement with simulations and empirical evidence), and (3) allow for the reduced viscosity after the disk expands, cools, and forms dust. These ingredients conspire to significantly lower the efficiency of eccentricity excitation by the disk tides. We find that explaining eccentricities of the post-AGB binaries is difficult and requires massive ($> 10^{-2}M_\odot$), long-lived ($> 10^5$ yr) circumbinary disks, that do not reaccrete. Reaccretion is detrimental for the eccentricity growth also in the MSP+WD systems. Reduced efficiency of the disk-driven excitation motivates study of alternative mechanisms for producing the peculiar eccentricities of the post-MS binaries. △ Less

Submitted 18 May, 2016; originally announced May 2016.

Comments: 9 pages, 4 figures, submitted to ApJ

Generalized Similarity for Accretion/Decretion Disks

Authors: Roman R. Rafikov

Abstract: Decretion (or external) disks are gas disks freely expanding to large radii due to their internal stresses. They are expected to naturally arise in tidal disruption events, around Be stars, in mass-losing post main sequence binaries, as a result of supernova fallback, etc. Their evolution is theoretically understood in two regimes: when the central object does not exert torque on the disk (a stand… ▽ More Decretion (or external) disks are gas disks freely expanding to large radii due to their internal stresses. They are expected to naturally arise in tidal disruption events, around Be stars, in mass-losing post main sequence binaries, as a result of supernova fallback, etc. Their evolution is theoretically understood in two regimes: when the central object does not exert torque on the disk (a standard assumption for conventional accretion disks) or when no mass inflow (or outflow) occurs at the disk center. However, many astrophysical objects - circumbinary disks, Be stars, neutron stars accreting in a propeller regime, etc. - feature non-zero torque simultaneously with the non-zero accretion (or ejection of mass) at the disk center. We provide a general description for the evolution of such disks (both linear and non-linear) in the self-similar regime, to which the disk should asymptotically converge with time. We identify a similarity parameter $λ$, which is uniquely related to the degree, to which the central mass accretion is suppressed by the non-zero central torque. The known decretion disk solutions correspond to the two discrete values of $λ$, while our new solutions cover a continuum of its physically allowed values, corresponding to either accretion or mass ejection by the central object. A direct relationship between $λ$ and central $\dot M$ and torque is also established. We describe the time evolution of the various disk characteristics for different $λ$, and show that the observable properties (spectrum and luminosity evolution) of the decretion disks are in general different from the standard accretion disks with no central torque. △ Less

Submitted 25 April, 2016; originally announced April 2016.

Comments: 16 pages, 13 figures, submitted to ApJ

Accretion and Orbital Inspiral in Gas-Assisted Supermassive Black Hole Binary Mergers

Authors: Roman R. Rafikov

Abstract: Many galaxies are expected to harbor binary supermassive black holes (SMBHs) in their centers. Their interaction with the surrounding gas results in accretion and exchange of angular momentum via tidal torques, facilitating binary inspiral. Here we explore the non-trivial coupling between these two processes and analyze how the global properties of externally supplied circumbinary disks depend on… ▽ More Many galaxies are expected to harbor binary supermassive black holes (SMBHs) in their centers. Their interaction with the surrounding gas results in accretion and exchange of angular momentum via tidal torques, facilitating binary inspiral. Here we explore the non-trivial coupling between these two processes and analyze how the global properties of externally supplied circumbinary disks depend on the binary accretion rate. By formulating our results in terms of the angular momentum flux driven by internal stresses, we come up with a very simple classification of the possible global disk structures, which differ from the standard constant $\dot M$ accretion disk solution. Suppression of accretion by the binary tides, leading to a significant mass accumulation in the inner disk, accelerates binary inspiral. We show that once the disk region strongly perturbed by the viscously transmitted tidal torque exceeds the binary semi-major axis, the binary can merge in less than its mass-doubling time due to accretion. Thus, unlike the inspirals driven by stellar scattering, the gas-assisted merger can occur even if the binary is embedded in a relatively low mass disk (lower than its own mass). This is important for resolving the "last parsec" problem for SMBH binaries and understanding powerful gravitational wave sources in the Universe. We argue that the enhancement of accretion by the binary found in some recent simulations cannot persist for a long time and should not affect the long-term orbital inspiral. We also review the existing simulations of the SMBH binary-disk coupling and propose a numerical setup, which is particularly well suited for verifying our theoretical predictions. △ Less

Submitted 16 February, 2016; originally announced February 2016.

Comments: 12 pages, 3 figures, submitted to ApJ

Protoplanetary Disk Heating and Evolution Driven by the Spiral Density Waves

Authors: Roman R. Rafikov

Abstract: High-resolution imaging of some protoplanetary disks in scattered light reveals presence of the global spiral arms of significant amplitude, likely excited by massive planets or stellar companions. Assuming that these arms are density waves, evolving into spiral shocks, we assess their effect on the thermodynamics, accretion, and global evolution of the disk. We derive analytical expressions for t… ▽ More High-resolution imaging of some protoplanetary disks in scattered light reveals presence of the global spiral arms of significant amplitude, likely excited by massive planets or stellar companions. Assuming that these arms are density waves, evolving into spiral shocks, we assess their effect on the thermodynamics, accretion, and global evolution of the disk. We derive analytical expressions for the direct (irreversible) heating, angular momentum transport, and mass accretion rate induced by the disk shocks of arbitrary strength. We find these processes to be very sensitive to the shock amplitude. Focusing on the waves of moderate strength (density jump at the shock $ΔΣ/Σ\sim 1$) we show the associated disk heating to be negligible (contributing at $\sim 1\%$ level to the energy budget) in passive, irradiated protoplanetary disks on $\sim 100$ AU scales, but becoming important within several AU from the star. At the same time, shock heating can be a significant (or even dominant) energy source in disks of cataclysmic variables, stellar X-ray binaries, and supermassive black hole binaries, heated mainly by viscous dissipation. Mass accretion induced by the global spiral shocks is comparable to (or exceeds) the mass inflow due to viscous stresses. Protoplanetary disks featuring prominent global spirals must be evolving rapidly, in $\lesssim 0.5$ Myr at $\sim 100$ AU. A direct upper limit on the disk evolution timescale can be established via the measurement of the gravitational torque due to the spiral arms from the imaging data. Our findings suggest that, regardless of their origin, global spiral waves must be important agents of the protoplanetary disk evolution. They may serve as an effective mechanism of disk dispersal and could be related to the transitional disk phenomenon. △ Less

Submitted 12 January, 2016; originally announced January 2016.

Comments: 11 pages, 2 figures, submitted to ApJ

Spreading layers in accreting objects: role of acoustic waves for angular momentum transport, mixing and thermodynamics

Authors: Alexander A. Philippov, Roman R. Rafikov, James M. Stone

Abstract: Disk accretion at high rate onto a white dwarf or a neutron star has been suggested to result in the formation of a spreading layer (SL) - a belt-like structure on the object's surface, in which the accreted matter steadily spreads in the poleward (meridional) direction while spinning down. To assess its basic characteristics we perform two-dimensional hydrodynamic simulations of supersonic SLs in… ▽ More Disk accretion at high rate onto a white dwarf or a neutron star has been suggested to result in the formation of a spreading layer (SL) - a belt-like structure on the object's surface, in which the accreted matter steadily spreads in the poleward (meridional) direction while spinning down. To assess its basic characteristics we perform two-dimensional hydrodynamic simulations of supersonic SLs in the relevant morphology with a simple prescription for cooling. We demonstrate that supersonic shear naturally present at the base of the SL inevitably drives sonic instability that gives rise to large scale acoustic modes governing the evolution of the SL. These modes dominate the transport of momentum and energy, which is intrinsically global and cannot be characterized via some form of local effective viscosity (e.g. $α$-viscosity). The global nature of the wave-driven transport should have important implications for triggering Type I X-ray bursts in low mass X-ray binaries. The nonlinear evolution of waves into a system of shocks drives effective re-arrangement (sensitively depending on thermodynamical properties of the flow) and deceleration of the SL, which ultimately becomes transonic and susceptible to regular Kelvin-Helmholtz instability. We interpret this evolution in terms of the global structure of the SL and suggest that mixing of the SL material with the underlying stellar fluid should become effective only at intermediate latitudes on the accreting object's surface, where the flow has decelerated appreciably. In the near-equatorial regions the transport is dominated by acoustic waves and mixing is less efficient. We speculate that this latitudinal non-uniformity of mixing in accreting white dwarfs may be linked to the observed bipolar morphology of classical novae ejecta. △ Less

Submitted 29 November, 2015; originally announced November 2015.

Comments: 20 pages, 25 figures, accepted to ApJ

Tatooine Nurseries: Structure and Evolution of Circumbinary Protoplanetary Disks

Authors: David Vartanyan, Jose A. Garmilla, Roman R. Rafikov

Abstract: Recent discoveries of circumbinary planets by Kepler mission provide motivation for understanding their birthplaces - protoplanetary disks around stellar binaries with separations <1 AU. We explore properties and evolution of such circumbinary disks focusing on modification of their structure caused by tidal coupling to the binary. We develop a set of analytical scaling relations describing viscou… ▽ More Recent discoveries of circumbinary planets by Kepler mission provide motivation for understanding their birthplaces - protoplanetary disks around stellar binaries with separations <1 AU. We explore properties and evolution of such circumbinary disks focusing on modification of their structure caused by tidal coupling to the binary. We develop a set of analytical scaling relations describing viscous evolution of the disk properties, which are verified and calibrated using 1D numerical calculations with realistic inputs. Injection of angular momentum by the central binary suppresses mass accretion onto the binary and causes radial distribution of the viscous angular momentum flux F_J to be different from that in a standard accretion disk around a single star with no torque at the center. Disks with no mass accretion at the center develop F_J profile which is flat in radius. Radial profiles of temperature and surface density are also quite different from those in disks around single stars. Damping of the density waves driven by the binary and viscous dissipation dominate heating of the inner disk (within 1-2 AU), pushing the iceline beyond 3-5 AU, depending on disk mass and age. Irradiation by the binary governs disk thermodynamics beyond ~10 AU. However, self-shadowing by the hot inner disk may render central illumination irrelevant out to ~20 AU. Spectral energy distribution of a circumbinary disk exhibits a distinctive bump around 10 micron, which may facilitate identification of such disks around unresolved binaries. Efficient tidal coupling to the disk drives orbital inspiral of the binary and may cause low-mass and compact binaries to merge into a single star within the disk lifetime. We generally find that circumbinary disks present favorable sites for planet formation (despite wider zone of volatile depletion), in agreement with the statistics of Kepler circumbinary planets. △ Less

Submitted 29 September, 2015; v1 submitted 24 September, 2015; originally announced September 2015.

Comments: 19 pages, 16 figures, submitted to ApJ, minor typos corrected