skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 10:00 PM ET on Friday, February 6 until 10:00 AM ET on Saturday, February 7 due to maintenance. We apologize for the inconvenience.


Title: The Santa Barbara Binary−disk Code Comparison
We have performed numerical calculations of a binary interacting with a gas disk, using 11 different numerical methods and a standard binary−disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary−disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, 0.2(a2ΩB/ν) binary orbits, whereνis the kinematic viscosity. The largest discrepancies between codes resulted from the dimensionality of the setup (3D vs. 2D disks). We find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.  more » « less
Award ID(s):
2206299 2006176
PAR ID:
10537734
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
970
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
156
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We examine the geometry of the post–asymptotic giant branch (AGB) star binary AC Her and its circumbinary disk. We show that the observations describe a binary orbit that is perpendicular to the disk with an angular momentum vector that is within 9° of the binary eccentricity vector, meaning that the disk is close to a stable polar alignment. The most likely explanation for the very large inner radius of the dust is a planet within the circumbinary disk. This is therefore both the first reported detection of a polar circumbinary disk around a post-AGB binary and the first evidence of a polar circumbinary planet. We consider the dynamical constraints on the circumbinary disk size and mass. The polar circumbinary disk feeds circumstellar disks with gas on orbits that are highly inclined with respect to the binary orbit plane. The resulting circumstellar disk inclination could be anywhere from coplanar to polar depending upon the competition between the mass accretion and binary torques. 
    more » « less
  2. Abstract Pebble accretion is recognized as a significant accelerator of planet formation. Yet only formulae for single-sized (monodisperse) distribution have been derived in the literature. These can lead to significant underestimates for Bondi accretion, for which the best accreted pebble size may not be the one that dominates the mass distribution. We derive in this paper the polydisperse theory of pebble accretion. We consider a power-law distribution in pebble radius, and we find the resulting surface and volume number density distribution functions. We derive also the exact monodisperse analytical pebble accretion rate for which 3D accretion and 2D accretion are limits. In addition, we find analytical solutions to the polydisperse 2D Hill and 3D Bondi limits. We integrate the polydisperse pebble accretion numerically for the MRN distribution, finding a slight decrease (by an exact factor 3/7) in the Hill regime compared to the monodisperse case. In contrast, in the Bondi regime, we find accretion rates 1–2 orders of magnitude higher compared to monodisperse, also extending the onset of pebble accretion to 1–2 orders of magnitude lower in mass. We find megayear timescales, within the disk lifetime, for Bondi accretion on top of planetary seeds of masses 10−6to 10−4M, over a significant range of the parameter space. This mass range overlaps with the high-mass end of the planetesimal initial mass function, and thus pebble accretion is possible directly following formation by streaming instability. This alleviates the need for mutual planetesimal collisions as a major contribution to planetary growth. 
    more » « less
  3. Abstract The architecture of planetary systems depends on the evolution of the disks in which they form. In this work, we develop a population synthesis approach to interpret the Atacama Large Millimeter/submillimeter Array survey of Gas Evolution of PROtoplanetary Disks (AGE-PRO) measurements of disk gas mass and size considering two scenarios: turbulence-driven evolution with photoevaporative winds and MHD wind-driven evolution. A systematic method is proposed to constrain the distribution of disk parameters from the disk fractions, accretion rates, disk gas masses, and CO gas sizes. We find that turbulence-driven accretion with initially compact disks (R0 ≃ 5–20 au), low mass-loss rates, and relatively long viscous timescales (tν,0 ≃ 0.4–3 Myr orαSS ≃ 2–4 × 10−4) can reproduce the disk fractions and gas sizes. However, the distribution of apparent disk lifetimes defined as the M D / M ̇ * ratio is severely overestimated by turbulence-driven models. On the other hand, MHD wind-driven accretion can reproduce the bulk properties of disk populations from Ophiuchus to Upper Scorpius assuming compact disks with an initial magnetization of aboutβ ≃ 105DW ≃ 0.5–1 × 10−3) and a magnetic field that declines with time. More studies are needed to confirm the low masses found by AGE-PRO, notably for compact disks that question turbulence-driven accretion. The constrained synthetic disk populations can now be used for realistic planet population models to interpret the properties of planetary systems on a statistical basis. 
    more » « less
  4. Abstract We analyze accretion-rate time series for equal-mass binaries in coplanar gaseous disks spanning a continuous range of orbital eccentricities up to 0.8 for both prograde and retrograde systems. The dominant variability timescales match those of previous investigations; the binary orbital period is dominant for prograde binaries withe≳ 0.1, with a 5 × longer “lump” period taking over fore≲ 0.1. This lump period fades and drops from 5 × to 4.5 × the binary period aseapproaches 0.1, where it vanishes. For retrograde orbits, the binary orbital period dominates ate≲ 0.55 and is accompanied by a 2 × longer timescale periodicity at higher eccentricities. The shape of the accretion-rate time series varies with binary eccentricity. For prograde systems, the orientation of an eccentric disk causes periodic trading of accretion between the binary components in a ratio that we report as a function of binary eccentricity. We present a publicly available tool,binlite, that can rapidly (≲0.01 s) generate templates for the accretion-rate time series onto either binary component for choice of binary eccentricity below 0.8. As an example use case, we build lightcurve models where the accretion rate through the circumbinary disk and onto each binary component sets contributions to the emitted specific flux. We combine these rest-frame, accretion-variability lightcurves with observer-dependent Doppler boosting and binary self-lensing. This allows a flexible approach to generating lightcurves over a wide range of binary and observer parameter space. We envisionbinliteas the access point to a living database that will be updated with state-of-the-art hydrodynamical calculations as they advance. 
    more » « less
  5. Abstract We present an intensive multiwavelength monitoring campaign of the quasar PG 1302−102 with Swift and the Las Cumbres Observatory network telescopes. Atz∼ 0.3, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time the RM technique has been applied to test disk structures predicted in the SMBH binary model that has been suggested for this source. PG 1302−102 was observed at a ∼daily cadence for ∼9 months in 14 bands spanning from X-ray to UV and optical wavelengths, and it shows moderate to significant levels of variability correlated between wavelengths. We measure the interband time lags, which are consistent with aτ∝λ4/3relation as expected from standard disk reprocessing, albeit with large uncertainties. The disk size implied by the lag spectrum is consistent with the expected disk size for its black hole mass within uncertainties. While the source resembles other reverberation-mapped active galactic nuclei in many respects, and we do not find evidence supporting the prevalent hypothesis that it hosts an SMBH binary, we demonstrate the feasibility of studying SMBH binaries from this novel angle and suggest possibilities for the LSST Deep Drilling Fields. 
    more » « less