BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Roman R. Rafikov

Roman R. Rafikov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GAastro-ph.HEastro-phastro-ph.CO

Catalog footprint

What is connected

42works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

42 published item(s)

preprint2023arXiv

Vortex weighing and dating of planets in protoplanetary discs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Doppler Boosting of the S-stars in the Galactic Center

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Eclipse timing the Milky Way's gravitational potential

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Planet-disk interaction in disks with cooling: basic theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

O(N^2) fragmentation algorithm

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Generalized Similarity for Accretion/Decretion Disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Eccentricity Excitation in Post-Main Sequence Binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Protoplanetary Disk Heating and Evolution Driven by the Spiral Density Waves

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Birth Locations of the Kepler Circumbinary Planets

The Kepler mission has discovered about a dozen circumbinary planetary systems, all containing planets on ~ 1 AU orbits. We place bounds on the locations in the circumbinary protoplanetary disk, where these planets could have formed through collisional agglomeration starting from small (km-sized or less) planetesimals. We first present a model of secular planetesimal dynamics that accounts for the (1) perturbation due to the eccentric precessing binary, as well as the (2) gravity and (3) gas drag from a precessing eccentric disk. Their simultaneous action leads to rich dynamics, with (multiple) secular resonances emerging in the disk. We derive analytic results for size-dependent planetesimal eccentricity, and demonstrate the key role of the disk gravity for circumbinary dynamics. We then combine these results with a simple model for collisional outcomes and find that in systems like Kepler 16, planetesimal growth starting with 10-100 m planetesimals is possible outside a few AU. The exact location exterior to which this happens is sensitive to disk eccentricity, density and precession rate, as well as to the size of the first generation of planetesimals. Strong perturbations from the binary in the inner part of the disk, combined with a secular resonance at a few AU inhibit the growth of km-sized planetesimals within 2 - 4 AU of the binary. In situ planetesimal growth in the Kepler circumbinary systems is possible only starting from large (few-km-sized) bodies in a low-mass disk with surface density less than 500 g/cm^2 at 1 AU.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tatooine Nurseries: Structure and Evolution of Circumbinary Protoplanetary Disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Structure of Spiral Shocks Excited by Planetary-mass Companions

Direct imaging observations have revealed spiral structures in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain some of these spiral patterns, due to the large pitch angle and high contrast of the spiral arms in observations. We have carried out three dimensional (3-D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the non-linear density wave theory. A secondary (or even a tertiary) spiral arm, especially for inner arms, is also excited by a massive planet. With a more massive planet in the disk, the excited spiral arms have larger pitch angle and the separation between the primary and secondary arms in the azimuthal direction is also larger. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3-D hydrodynamical models than images based on 2-D models assuming vertical hydrostatic equilibrium, indicating the need to model observations with full 3-D hydrodynamics. Overall, companion-induced spiral arms not only pinpoint the companion's position but also provide three independent ways (pitch angle, separation between two arms, and contrast of arms) to constrain the companion's mass.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Viscosity prescription for gravitationally unstable accretion disks

Gravitationally unstable accretion disks emerge in a variety of astrophysical contexts - giant planet formation, FU Orioni outbursts, feeding of AGNs, and the origin of Pop III stars. When a gravitationally unstable disk is unable to cool rapidly it settles into a quasi-stationary, fluctuating gravitoturbulent state, in which its Toomre Q remains close to a constant value Q_0~1. Here we develop an analytical formalism describing the evolution of such a disk, which is based on the assumptions of Q=Q_0 and local thermal equilibrium. Our approach works in the presence of additional sources of angular momentum transport (e.g. MRI), as well as external irradiation. Thermal balance dictates a unique value of the gravitoturbulent stress α_{gt} driving disk evolution, which is a function of the local surface density and angular frequency. We compare this approach with other commonly used gravitoturbulent viscosity prescriptions, which specify the explicit dependence of stress α_{gt} on Toomre Q in an ad hoc fashion, and identify the ones that provide consistent results. We nevertheless argue that our Q=Q_0 approach is more flexible, robust, and straightforward, and should be given preference in applications. We illustrate this with a couple of analytical calculations - locations of the snow line and of the outer edge of the dead zone in a gravitoturbulent protoplanetary disk - which clearly show the simplicity and versatility of the Q=Q_0 approach.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Particle Concentration At Planet Induced Gap Edges and Vortices: I. Inviscid 3-D Hydro Disks

We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global 2-D and 3-D inviscid hydrodynamic simulations. Lagrangian particles have been implemented into magnetohydrodynamic code Athena with cylindrical coordinates. We find two distinct outcomes depending on the mass of the embedded planet. In the presence of a low mass planet ($8 M_{\oplus}$), two narrow gaps start to open in the gas on each side of the planet where the density waves shock. These shallow gaps can dramatically affect particle drift speed and cause significant, roughly axisymmetric dust depletion. On the other hand, a more massive planet ($>0.1 M_{J}$) carves out a deeper gap with sharp edges, which are unstable to the vortex formation. Particles with a wide range of sizes ($0.02<Ωt_{s}<20$) are trapped and settle to the midplane in the vortex, with the strongest concentration for particles with $Ωt_{s}\sim 1$. The dust concentration is highly elongated in the $ϕ$ direction, and can be as wide as 4 disk scale heights in the radial direction. Dust surface density inside the vortex can be increased by more than a factor of 10$^2$ in a very non-axisymmetric fashion. For very big particles ($Ωt_{s}\gg 1$) we find strong eccentricity excitation, in particular around the planet and in the vicinity of the mean motion resonances, facilitating gap opening there. Our results imply that in weakly turbulent protoplanetary disk regions (e.g. the "dead zone") dust particles with a very wide range of sizes can be trapped at gap edges and inside vortices induced by planets with $M_{p}<M_{J}$, potentially accelerating planetesimal and planet formation there, and giving rise to distinctive features that can be probed by ALMA and EVLA.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Planet formation in binaries: dynamics of planetesimals perturbed by the eccentric protoplanetary disk and the secondary

Detections of planets in eccentric, close (separations of ~20 AU) binary systems such as αCen or γCep provide an important test of planet formation theories. Gravitational perturbations from the companion are expected to excite high planetesimal eccentricities resulting in destruction, rather than growth, of objects with sizes of up to several hundred km in collisions of similar-size bodies. It was recently suggested that gravity of a massive axisymmetric gaseous disk in which planetesimals are embedded drives rapid precession of their orbits, suppressing eccentricity excitation. However, disks in binaries are themselves expected to be eccentric, leading to additional planetesimal excitation. Here we develop secular theory of eccentricity evolution for planetesimals perturbed by the gravity of an elliptical protoplanetary disk (neglecting gas drag) and the companion. For the first time we derive an expression for the disturbing function due to an eccentric disk, which can be used for a variety of other astrophysical problems. We obtain explicit analytical solutions for planetesimal eccentricity evolution and delineate four different regimes of dynamical excitation. We show that in systems with massive (>10^{-2}M_Sun) disks, planetesimal eccentricity is usually determined by the gravity of the eccentric disk alone, and is comparable to the disk eccentricity. As a result, the latter imposes a lower limit on collisional velocities of solids, making their growth problematic. This fragmentation barrier can be removed if the gaseous disk rapidly precesses or if its own self-gravity is efficient at lowering disk eccentricity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Planet formation in stellar binaries I: planetesimal dynamics in massive protoplanetary disks

About $20\%$ of exoplanets discovered by radial velocity surveys reside in stellar binaries. To clarify their origin one has to understand the dynamics of planetesimals in protoplanetary disks within binaries. The standard description, accounting for only gas drag and gravity of the companion star has been challenged recently, as the gravity of the protoplanetary disk was shown to play a crucial role in planetesimal dynamics. An added complication is the tendency of protoplanetary disks in binaries to become eccentric, giving rise to additional excitation of planetesimal eccentricity. Here, for the first time, we analytically explore secular dynamics of planetesimals in binaries such as $α$ Cen and $γ$ Cep under the combined action of (1) gravity of the eccentric protoplanetary disk, (2) perturbations due to the (coplanar) eccentric companion, and (3) gas drag. We derive explicit solutions for the behavior of planetesimal eccentricity ${\bf e}_p$ in non-precessing disks (and in precessing disks in certain limits). We obtain the analytical form of the distribution of relative velocities of planetesimals, which is a key input for understanding their collisional evolution. Disk gravity strongly influences relative velocities and tends to push sizes of planetesimals colliding with comparable objects at the highest speed to small values, $\sim 1$ km. We also find that planetesimals in eccentric protoplanetary disks apsidally aligned with the binary orbit collide at lower relative velocities than in mis-aligned disks. Our results highlight a decisive role that disk gravity plays in planetesimal dynamics in binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Planet formation in stellar binaries II: overcoming the fragmentation barrier in alpha Centauri and gamma Cephei-like systems

Planet formation in small-separation (~20 AU) eccentric binaries such as gamma Cephei or alpha Centauri is believed to be adversely affected by the presence of the stellar companion. Strong dynamical excitation of planetesimals by the eccentric companion can result in collisional destruction (rather than growth) of 1-100 km objects, giving rise to the "fragmentation barrier" for planet formation. We revise this issue using a novel description of secular dynamics of planetesimals in binaries, which accounts for the gravity of the eccentric, coplanar protoplanetary disk, as well as gas drag. By studying planetesimal collision outcomes we show, in contrast to many previous studies, that planetesimal growth and subsequent formation of planets (including gas giants) in AU-scale orbits within ~20 AU separation binaries may be possible, provided that the protoplanetary disks are massive (>10^{-2}M_\odot) and only weakly eccentric (disk eccentricity <0.01). These requirements are compatible with both the existence of massive (several M_J) planets in gamma Cep-like systems and the results of recent simulations of gaseous disks in eccentric binaries. Terrestrial and Neptune-like planets can also form in lower-mass disks at small (sub-AU) radii. We find that fragmentation barrier is less of a problem in eccentric disks which are apsidally aligned with the binary orbit. Alignment gives rise to special locations, where (1) relative planetesimal velocities are low and (2) the timescale of their drag-induced radial drift is long. This causes planetesimal pileup at such locations in the disk and promotes their growth.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probing Oort Cloud and Local Interstellar Medium Properties via Dust Produced in Cometary Collisions

The Oort Cloud remains one of the most poorly explored regions of the Solar System. We propose that its properties can be constrained by studying a population of dust grains produced in collisions of comets in the outer Solar System. We explore the dynamics of micron-size grains outside the heliosphere (beyond ~250 AU), which are affected predominantly by the magnetic field of the interstellar medium (ISM) flow past the Sun. We derive analytic models for the production and motion of small particles as a function of their birth location in the Cloud and calculate the particle flux and velocity distribution in the inner Solar System. These models are verified by direct numerical simulations. We show that grains originating in the Oort Cloud have a unique distribution of arrival directions, which should easily distinguish them from both interplanetary and interstellar dust populations. We also demonstrate that the distribution of particle arrival velocities is uniquely determined the mass distribution and dust production rate in the Cloud. Cometary collisions within the Cloud produce a flux of micron-size grains in the inner Solar System of up to several per square meter per year. The next-generation dust detectors may be sensitive enough to detect and constrain this dust population, which will illuminate us about the Oort Cloud's properties. We also show that the recently-detected mysterious population of large (micron-size) unbound particles, which seems to arrive with the ISM flow is unlikely to be generated by collisions of comets in the Oort Cloud.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scattering outcomes of close-in planets: constraints on planet migration

Many exoplanets in close-in orbits are observed to have relatively high eccentricities and large stellar obliquities. We explore the possibility that these result from planet-planet scattering by studying the dynamical outcomes from a large number of orbit integrations in systems with two and three gas-giant planets in close-in orbits (0.05 AU < a < 0.15 AU). We find that at these orbital separations, unstable systems starting with low eccentricities and mutual inclinations ($e\lesssim0.1$, $i\lesssim0.1$) generally lead to planet-planet collisions in which the collision product is a planet on a low-eccentricity, low-inclination orbit. This result is inconsistent with the observations. We conclude that eccentricity and inclination excitation from planet-planet scattering must precede migration of planets into short-period orbits. This result constrains theories of planet migration: the semi-major axis must shrink by 1-2 orders of magnitude without damping the eccentricity and inclination.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Secular dynamics of the triple system harbouring PSR J0337+1715 and implications for the origin of its orbital configuration

We explore secular dynamics of a recently discovered hierarchical triple system consisting of the radio pulsar PSR J0337+1715 and two white dwarfs (WDs). We show that three body interactions endow the inner binary with a large forced eccentricity and suppress its apsidal precession, to about 24% of the rate due to the general relativity. However, precession rate is still quite sensitive to the non-Newtonian effects and may be used to constrain gravity theories if measured accurately. Small value of the free eccentricity of the inner binary $e_{i}^{free}\approx 2.6\times 10^{-5}$ and vanishing forced eccentricity of the outer, relatively eccentric binary naturally result in their apsidal near-alignment. In addition, this triple system provides a unique opportunity to explore excitation of both eccentricity and inclination in neutron star-WD binaries, e.g. due to random torques caused by convective eddies in the WD progenitor. We show this process to be highly anisotropic and more effective at driving eccentricity rather than inclination. The outer binary eccentricity as well as $e_{i}^{free}$ exceed by more than an order of magnitude the predictions of the eccentricity-period relation of Phinney (1992), which is not uncommon. We also argue that the non-zero mutual inclination of the two binaries emerges at the end of the Roche lobe overflow of the outer (rather than the inner) binary.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analysis of Spin-Orbit Misalignment in Eclipsing Binary DI Herculis

Eclipsing binary DI Herculis (DI Her) is known to exhibit anomalously slow apsidal precession, below the rate predicted by the general relativity. Recent measurements of the Rossiter-McLauglin effect indicate that stellar spins in DI Her are almost orthogonal to the orbital angular momentum, which explains the anomalous precession in agreement with the earlier theoretical suggestion by Shakura. However, these measurements yield only the projections of the spin-orbit angles onto the sky plane, leaving the spin projection onto our line of sight unconstrained. Here we describe a method of determining the full three-dimensional spin orientation of the binary components relying on the use of the gravity darkening effect, which is significant for the rapidly rotating stars in DI Her. Gravity darkening gives rise to nonuniform brightness distribution over the stellar surface, the pattern of which depends on the stellar spin orientation. Using archival photometric data obtained during multiple eclipses spread over several decades we are able to constrain the unknown spin angles in DI Her with this method, finding that spin axes of both stars lie close to the plane of the sky. Our procedure fully accounts for the precession of stellar spins over the long time span of observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Angular Momentum Transport by Acoustic Modes Generated in the Boundary Layer II: MHD Simulations

We perform global unstratified 3D magnetohydrodynamic simulations of an astrophysical boundary layer (BL) -- an interface region between an accretion disk and a weakly magnetized accreting object such as a white dwarf -- with the goal of understanding the effects of magnetic field on the BL. We use cylindrical coordinates with an isothermal equation of state and investigate a number of initial field geometries including toroidal, vertical, and vertical with zero net flux. Our initial setup consists of a Keplerian disk attached to a non-rotating star. In a previous work, we found that in hydrodynamical simulations, sound waves excited by shear in the BL were able to efficiently transport angular momentum and drive mass accretion onto the star. Here we confirm that in MHD simulations, waves serve as an efficient means of angular momentum transport in the vicinity of the BL, despite the magnetorotational instability (MRI) operating in the disk. In particular, the angular momentum current due to waves is at times larger than the angular momentum current due to MRI. Our results suggest that angular momentum transport in the BL and its vicinity is a global phenomenon occurring through dissipation of waves and shocks. This point of view is quite different from the standard picture of transport by a local anomalous turbulent viscosity. In addition to angular momentum transport, we also study magnetic field amplification within the BL. We find that the field is indeed amplified in the BL, but only by a factor of a few and remains subthermal.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Low mass planets in protoplanetary disks with net vertical magnetic fields: the Planetary Wake and Gap Opening

We study wakes and gap opening by low mass planets in gaseous protoplanetary disks threaded by net vertical magnetic fields which drive magnetohydrodynamical (MHD) turbulence through the magnetorotational instabilty (MRI), using three dimensional simulations in the unstratified local shearing box approximation. The wakes, which are excited by the planets, are damped by shocks similar to the wake damping in inviscid hydrodynamic (HD) disks. Angular momentum deposition by shock damping opens gaps in both MHD turbulent disks and inviscid HD disks even for low mass planets, in contradiction to the "thermal criterion" for gap opening. To test the "viscous criterion", we compared gap properties in MRI-turbulent disks to those in viscous HD disks having the same stress, and found that the same mass planet opens a significantly deeper and wider gap in net vertical flux MHD disks than in viscous HD disks. This difference arises due to the efficient magnetic field transport into the gap region in MRI disks, leading to a larger effective αwithin the gap. Thus, across the gap, the Maxwell stress profile is smoother than the gap density profile, and a deeper gap is needed for the Maxwell stress gradient to balance the planetary torque density. We also confirmed the large excess torque close to the planet in MHD disks, and found that long-lived density features (termed zonal flows) produced by the MRI can affect planet migration. The comparison with previous results from net toroidal flux/zero flux MHD simulations indicates that the magnetic field geometry plays an important role in the gap opening process. Overall, our results suggest that gaps can be commonly produced by low mass planets in realistic protoplanetary disks, and caution the use of a constant α-viscosity to model gaps in protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Search for Supermassive Black Hole Binaries in the Sloan Digital Sky Survey Spectroscopic Sample

Supermassive black hole (SMBH) binaries are expected in a Lambda CDM cosmology given that most (if not all) massive galaxies contain a massive black hole at their center. So far, however, direct evidence for such binaries has been elusive. We use cross-correlation to search for temporal velocity shifts in the MgII broad emission lines of 0.36 < z < 2 quasars with multiple observations in the Sloan Digital Sky Survey. For ~ 10^9 Msun BHs in SMBH binaries, we are sensitive to velocity drifts for binary separations of ~ 0.1 pc with orbital periods of ~100 years. We find seven candidate sub-pc--scale binaries with velocity shifts > 3.4 sigma ~ 280 km/s, where sigma is our systematic error. Comparing the detectability of SMBH binaries with the number of candidates (N < 7), we can rule out that most 10^9 Msun BHs exist in ~ 0.03-0.2 pc scale binaries, in a scenario where binaries stall at sub-pc scales for a Hubble time. We further constrain that < one-third of quasars host SMBH binaries after considering gas-assisted sub-pc evolution of SMBH binaries, although this result is very sensitive to the assumed size of the broad line region. We estimate the detectability of SMBH binaries with ongoing or next-generation surveys (e.g., BOSS, Subaru Prime Focus Spectrograph), taking into account the evolution of the sub-parsec binary in circumbinary gas disks. These future observations will provide longer time baselines for searches similar to ours and may in turn constrain the evolutionary scenarios of SMBH binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Angular Momentum Transport and Variability in Boundary Layers of Accretion Disks Driven by Global Acoustic Modes

Disk accretion onto a weakly magnetized central object, e.g. a star, is inevitably accompanied by the formation of a boundary layer near the surface, in which matter slows down from the highly supersonic orbital velocity of the disk to the rotational velocity of the star. We perform high resolution 2D hydrodynamical simulations in the equatorial plane of an astrophysical boundary layer with the goal of exploring the dynamics of non-axisymmetric structures that form there. We generically find that the supersonic shear in the boundary layer excites non-axisymmetric quasi-stationary acoustic modes that are trapped between the surface of the star and a Lindblad resonance in the disk. These modes rotate in a prograde fashion, are stable for hundreds of orbital periods, and have a pattern speed that is less than and of order the rotational velocity at the inner edge of the disk. The origin of these intrinsically global modes is intimately related to the operation of a corotation amplifier in the system. Dissipation of acoustic modes in weak shocks provides a universal mechanism for angular momentum and mass transport even in purely hydrodynamic (i.e. non-magnetized) boundary layers. We discuss the possible implications of these trapped modes for explaining the variability seen in accreting compact objects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Disk-satellite interaction in disks with density gaps

Gravitational coupling between a gaseous disk and an orbiting perturber leads to angular momentum exchange between them which can result in gap opening by planets in protoplanetary disks and clearing of gas by binary supermassive black holes (SMBHs) embedded in accretion disks. Understanding the co-evolution of the disk and the orbit of the perturber in these circumstances requires knowledge of the spatial distribution of the torque exerted by the latter on a highly nonuniform disk. Here we explore disk-satellite interaction in disks with gaps in linear approximation both in Fourier and in physical space, explicitly incorporating the disk non-uniformity in the fluid equations. Density gradients strongly displace the positions of Lindblad resonances in the disk (which often occur at multiple locations), and the waveforms of modes excited close to the gap edge get modified compared to the uniform disk case. The spatial distribution of the excitation torque density is found to be quite different from the existing prescriptions: most of the torque is exerted in a rather narrow region near the gap edge where Lindblad resonances accumulate, followed by an exponential fall-off with the distance from the perturber. Despite these differences, for a given gap profile the full integrated torque exerted on the disk agrees with the conventional uniform disk theory prediction at the level of ~10%. The nonlinearity of the density wave excited by the perturber is shown to decrease as the wave travels out of the gap, slowing down its nonlinear evolution and damping. Our results suggest that gap opening in protoplanetary disks and gas clearing around SMBH binaries can be more efficient than the existing theories predict. They pave the way for self-consistent calculations of the gap structure and the orbital evolution of the perturber using accurate prescription for the torque density behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Global Models of Runaway Accretion in White Dwarf Debris Disks

A growing sample of white dwarfs (WDs) with metal-enriched atmospheres are accompanied by excess infrared emission, indicating that they are encircled by a compact dusty disk of solid debris. Such `WD debris disks' are thought to originate from the tidal disruption of asteroids or other minor bodies, but the precise mechanism(s) responsible for transporting matter to the WD surface remains unclear, especially in those systems with the highest inferred metal accretion rates dM_Z/dt ~ 1e8-1e10 g/s. Here we present global time-dependent calculations of the coupled evolution of the gaseous and solid components of WD debris disks. Solids transported inwards (initially due to PR drag) sublimate at tens of WD radii, producing a source of gas that accretes onto the WD surface and viscously spreads outwards in radius, where it overlaps with the solid disk. If the aerodynamic coupling between the solids and gaseous disks is sufficiently strong (and/or the gas viscosity sufficiently weak), then gas builds up near the sublimation radius faster than it can viscously spread away. Since the rate of drag-induced solid accretion increases with gas density, this results in a runaway accretion process, during which the WD accretion rate reaches values orders of magnitude higher than can be achieved by PR drag alone. We explore the evolution of WD debris disks across a wide range of physical conditions and calculate the predicted distribution of observed accretion rates dM_Z/dt, finding reasonable agreement with the current sample. Although the conditions necessary for runaway accretion are at best marginally satisfied given the minimal level of aerodynamic drag between circular gaseous and solid disks, the presence of other stronger forms of solid-gas coupling---such as would result if the gaseous disk is only mildly eccentric---substantially increase the likelihood of runaway accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How to build Tatooine: reducing secular excitation in Kepler circumbinary planet formation

Circumbinary planetary systems recently discovered by Kepler represent an important testbed for planet formation theories. Planetesimal growth in disks around binaries has been expected to be inhibited interior to ~10 AU by secular excitation of high relative velocities between planetesimals, leading to their collisional destruction (rather than agglomeration). Here we show that gravity of the gaseous circumbinary disk in which planets form drives fast precession of both the planetesimal and binary orbits, resulting in strong suppression of planetesimal eccentricities beyond 2-3 AU and making possible growth of 1-100 km objects in this region. The precise location of the boundary of accretion-friendly region depends on the size of the inner disk cavity cleared by the binary torques and on the disk mass (even 0.01 M_Sun disk strongly suppresses planetesimal excitation), among other things. Precession of the orbit of the central binary, enhanced by the mass concentration naturally present at the inner edge of a circumbinary disk, plays key role in this suppression, which is a feature specific to the circumbinary planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Inner edges of compact debris disks around metal-rich white dwarfs

A number of metal-rich white dwarfs (WDs) are known to host compact, dense particle disks, which are thought to be responsible for metal pollution of these stars. In many such systems the inner radii of disks inferred from their spectra are so close to the WD that particles directly exposed to starlight must be heated above 1500 K and are expected to be unstable against sublimation. To reconcile this expectation with observations we explore particle sublimation in H-poor debris disks around WDs. We show that because of the high metal vapor pressure the characteristic sublimation temperature in these disks is 300-400 K higher than in their protoplanetary analogues, allowing particles to survive at higher temperatures. We then look at the structure of the inner edges of debris disks and show that they should generically feature superheated inner rims directly exposed to starlight with temperatures reaching 2500-3500 K. Particles migrating through the rim towards the WD (and rapidly sublimating) shield the disk behind them from strong stellar heating, making the survival of solids possible close to the WD. Our model agrees well with observations of WD+disk systems provided that disk particles are composed of Si-rich material such as olivine, and have sizes in the range ~(0.03-30) cm.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planet-Disk interaction in 3D: the importance of buoyancy waves

We carry out local three dimensional (3D) hydrodynamic simulations of planet-disk interaction in stratified disks with varied thermodynamic properties. We find that whenever the Brunt-Vaisala frequency (N) in the disk is nonzero, the planet exerts a strong torque on the disk in the vicinity of the planet, with a reduction in the traditional "torque cutoff". In particular, this is true for adiabatic perturbations in disks with isothermal density structure, as should be typical for centrally irradiated protoplanetary disks. We identify this torque with buoyancy waves, which are excited (when N is non-zero) close to the planet, within one disk scale height from its orbit. These waves give rise to density perturbations with a characteristic 3D spatial pattern which is in close agreement with the linear dispersion relation for buoyancy waves. The torque due to these waves can amount to as much as several tens of per cent of the total planetary torque, which is not expected based on analytical calculations limited to axisymmetric or low-m modes. Buoyancy waves should be ubiquitous around planets in the inner, dense regions of protoplanetary disks, where they might possibly affect planet migration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structure and evolution of circumbinary disks around supermassive black hole (SMBH) binaries

It is generally believed that gaseous disks around supermassive black hole (SMBH) binaries in centers of galaxies can facilitate binary merger and give rise to observational signatures both in electromagnetic and gravitational wave domains. We explore general properties of circumbinary disks by reformulating standard equations for the viscous disk evolution in terms of the viscous angular momentum flux F_J. In steady state F_J is a linear function of the specific angular momentum, which is a generalization of (but is not equivalent to) the standard constant \dot M disk solution. If the torque produced by the central binary is effective at stopping gas inflow and opening a gap (or cavity) in the disk, then the inner part of the circumbinary disk can be approximated as a constant F_J disk. We compute properties of such disks in different physical regimes relevant for SMBH binaries and use these results to understand the gas-assisted evolution of SMBH pairs starting at separations 10^{-4}-10^{-2} pc. We find the following. (1) Pile-up of matter at the inner edge of the disk leads to continuous growth of the torque acting on the binary and can considerably accelerate its orbital evolution compared to the gravitational wave-driven decay. (2) Torque on the binary is determined non-locally and does not in general reflect the disk properties in the vicinity of the binary. (3) Binary evolution depends on the past history of the disk evolution. (4) Eddington limit can be important in circumbinary disks even if they accrete at sub-Eddington rates at late stages of binary evolution. (5) Circumbinary disk self-consistently evolved under the action of the binary torque emits more power and has spectrum different from the spectrum of constant \dot M disk - it is steeper (νF_ν\propto ν^{12/7}) and extends to shorter wavelength, facilitating its detection.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Density Waves Excited by Low-Mass Planets in Protoplanetary Disks I: Linear Regime

Density waves excited by planets embedded in protoplanetary disks play a central role in planetary migration and gap opening processes. We carry out 2D shearing sheet simulations to study the linear regime of wave evolution with the grid-based code Athena, and provide detailed comparisons with the theoretical predictions. Low mass planets (down to ~0.03 Earth mass at 1 AU) and high spatial resolution (256 grid points per scale height) are chosen to mitigate the effects of wave nonlinearity. To complement the existing numerical studies, we focus on the primary physical variables such as the spatial profile of the wave, torque density, and the angular momentum flux carried by the wave, instead of secondary quantities such as the planetary migration rate. Our results show percent level agreement with theory in both physical and Fourier space. New phenomena such as the change of the toque density sign far from the planet are discovered and discussed. Also, we explore the effect of the numerical algorithms, and find that a high order of accuracy, high resolution, and an accurate planetary potential are crucial to achieve good agreement with the theory. We find that the use of a too large time-step without properly resolving the dynamical time scale around the planet produces incorrect results, and may lead to spurious gap opening. Global simulations of planet migration and gap opening violating this requirement may be affected by spurious effects resulting in e.g. the incorrect planetary migration rate and gap opening mass.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Density Waves Excited by Low-Mass Planets in Protoplanetary Disks II: High-Resolution Simulations of the Nonlinear Regime

We investigate numerically the propagation of density waves excited by a low-mass planet in a protoplanetary disk in the nonlinear regime, using 2D local shearing box simulations with the grid-based code Athena at high spatial resolution (256 grid points per scale height h). The nonlinear evolution results in the wave steepening into a shock, causing damping and angular momentum transfer to the disk. On long timescales this leads to spatial redistribution of the disk density, causing migration feedback and potentially resulting in gap opening. Previous numerical studies concentrated on exploring these secondary phenomena as probes of the nonlinear wave evolution. Here we focus on exploring the evolution of the basic wave properties, such as its density profile evolution, shock formation, post-shock wave behavior, and provide comparison with analytical theory. The generation of potential vorticity at the shock is computed analytically and is subsequently verified by simulations and used to pinpoint the shock location. We confirm the theoretical relation between the shocking length and the planet mass (including the effect of the equation of state), and the post-shock decay of the angular momentum flux carried by the wave. The post-shock evolution of the wave profile is explored, and we quantitatively confirm its convergence to the theoretically expected N-wave shape. The accuracy of various numerical algorithms used to compute the nonlinear wave evolution is also investigated: we find that higher order spatial reconstruction and high resolution are crucial for capturing the shock formation correctly.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Global modeling of radiatively driven accretion of metals from compact debris disks onto the white dwarfs

Recent infrared observations have revealed presence of compact (radii < R_Sun) debris disks around more than a dozen of metal-rich white dwarfs (WD), likely produced by tidal disruption of asteroids. Accretion of high-Z material from these disks may account for the metal contamination of these WDs. It was previously shown using local calculations that the Poynting-Robertson (PR) drag acting on the dense, optically thick disk naturally drives metal accretion onto the WD at the typical rate \dot M_{PR} \approx 10^8 g/s. Here we extend this local analysis by exploring global evolution of the debris disk under the action of the PR drag for a variety of assumptions about the disk properties. We find that massive disks (mass > 10^{20} g), which are optically thick to incident stellar radiation inevitably give rise to metal accretion at rates \dot M > 0.2\dot M_{PR}. The magnitude of \dot M and its time evolution are determined predominantly by the initial pattern of the radial distribution of the debris (i.e. ring-like vs. disk-like) but not by the total mass of the disk. The latter determines only the disk lifetime, which can be several Myr or longer. Evolution of an optically thick disk generically results in the development of a sharp outer edge of the disk. We also find that the low mass (< 10^{20} g), optically thin disks exhibit \dot M << \dot M_{PR} and evolve on characteristic timescale \sim 10^5-10^6 yr, independent of their total mass.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Metal accretion onto white dwarfs caused by Poynting-Robertson drag on their debris disks

Recent discoveries of compact (sizes $<R_\odot$) debris disks around more than a dozen of metal-rich white dwarfs (WDs) suggest that pollution of these stars with metals may be caused by accretion of high-Z material from the disk. But the mechanism responsible for efficient transfer of mass from a particulate disk to the WD atmosphere has not yet been identified. Here we demonstrate that radiation of the WD can effectively drive accretion of matter through the disk towards the sublimation radius (located at several tens of WD radii), where particles evaporate, feeding a disk of metal gas accreting onto the WD. We show that, contrary to some previous claims, Poynting-Robertson (PR) drag on the debris disk is effective at providing metal accretion rate $\dot M_{PR}\sim 10^8$ g/s and higher, scaling quadratically with WD effective temperature. We compare our results with observations and show that, as expected, no WD hosting a particulate debris disk shows evidence of metal accretion rate below that produced by the PR drag. Existence of WDs accreting metals at rates significantly higher than $\dot M_{PR}$ suggests that another mechanism in addition to PR drag drives accretion of high-Z elements in these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Runaway accretion of metals from compact debris disks onto white dwarfs

It was recently proposed that metal-rich white dwarfs (WDs) accrete their metals from compact debris disks found to exist around more than a dozen of them. At the same time, elemental abundances measured in atmospheres of some WDs imply vigorous metal accretion at rates up to $10^{11}$ g/s, far in excess of what can be supplied solely by Poynting-Robertson drag acting on such debris disks. To explain this observation we propose a model, in which rapid transport of metals from the disk onto the WD naturally results from interaction between this particulate disk and spatially coexisting disk of metallic gas. The latter is fed by evaporation of debris particles at the sublimation radius located at several tens of WD radii. Because of pressure support gaseous disk orbits WD slower than particulate disk. Resultant azimuthal drift between them at speed ~1 m/s causes aerodynamic drag on the disk of solids and drives inward migration of its constituent particles. Upon reaching the sublimation radius particles evaporate, enhancing the density of metallic gaseous disk and leading to positive feedback. Under favorable circumstances (low viscosity in the disk of metallic gas and efficient aerodynamic coupling between the disks) system evolves in a runaway fashion, destroying debris disk on time scale of $\sim 10^5$ yr, and giving rise to high metal accretion rates up to $10^{10}-10^{11}$ g/s, in agreement with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The origin of the negative torque density in disk-satellite interaction

Tidal interaction between a gaseous disk and a massive orbiting perturber is known to result in angular momentum exchange between them. Understanding astrophysical manifestations of this coupling such as gap opening by planets in protoplanetary disks or clearing of gas by binary supermassive black holes (SMBHs) embedded in accretion disks requires knowledge of the spatial distribution of the torque exerted on the disk by a perturber. Recent hydrodynamical simulations by Dong et al (2011) have shown evidence for the tidal torque density produced in a uniform disk to change sign at the radial separation of $\approx 3.2$ scale heights from the perturber's orbit, in clear conflict with the previous studies. To clarify this issue we carry out a linear calculation of the disk-satellite interaction putting special emphasis on understanding the behavior of the perturbed fluid variables in physical space. Using analytical as well as numerical methods we confirm the reality of the negative torque density phenomenon and trace its origin to the overlap of Lindblad resonances in the vicinity of the perturber's orbit - an effect not accounted for in previous studies. These results suggest that calculations of the gap and cavity opening in disks by planets and binary SMBHs should rely on more realistic torque density prescriptions than the ones used at present.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamical evolution of thin dispersion-dominated planetesimal disks

We study the dynamics of a vertically thin, dispersion-dominated disk of planetesimals with eccentricities $e$ and inclinations $i$ (normalized in Hill units) satisfying $e >> 1$, $i << e^{-2} << 1$. This situation may be typical for e.g. a population of protoplanetary cores in the end of the oligarchic phase of planet formation. In this regime of orbital parameters planetesimal scattering has an anisotropic character and strongly differs from scattering in thick ($i ~ e$) disks. We derive analytical expressions for the planetesimal scattering coefficients and compare them with numerical calculations. We find significant discrepancies in the inclination scattering coefficients obtained by the two approaches and ascribe this difference to the effects not accounted for in the analytical calculation: multiple scattering events (temporary captures, which may be relevant for the production of distant planetary satellites outside the Hill sphere) and distant interaction of planetesimals prior to their close encounter. Our calculations show that the inclination of a thin, dispersion-dominated planetesimal disk grows exponentially on a very short time scale implying that (1) such disks must be very short-lived and (2) planetesimal accretion in this dynamical phase is insignificant. Our results are also applicable to the dynamics of shear-dominated disks switching to the dispersion-dominated regime.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Convective cooling and fragmentation of gravitationally unstable disks

Gravitationally unstable disks can fragment and form bound objects provided that their cooling time is short. In protoplanetary disks radiative cooling is likely to be too slow to permit formation of planets by fragmentation within several tens of AU from the star. Recently, convection has been suggested as a faster means of heat loss from the disk but here we demonstrate that it is only marginally more efficient than the radiative cooling. The crucial factor is the rate at which energy can be radiated from the disk photosphere, which is robustly limited from above in the convective case by the adiabatic temperature gradient (given a certain midplane temperature). Thus, although vigorous convection is definitely possible in disks, the inefficiency of radiative loss from the photosphere may create a bottleneck limiting the ability of the disk to form self-gravitating objects. Based on this argument we derive a set of analytical constraints which diagnose the susceptibility of an unstable disk to fragmentation and show that the formation of giant planets by fragmentation of protoplanetary disks is unlikely to occur at distances of tens of AU. At the same time these constraints do not preclude the possibility of fragmentation and star formation in accretion disks around supermassive black holes.

0 citations0 reviewsNo code linkedOpen access