-
The Santa Barbara Binary-Disk Code Comparison
Authors:
Paul C. Duffell,
Alexander J. Dittmann,
Daniel J. D'Orazio,
Alessia Franchini,
Kaitlin M. Kratter,
Anna B. T. Penzlin,
Enrico Ragusa,
Magdalena Siwek,
Christopher Tiede,
Haiyang Wang,
Jonathan Zrake,
Adam M. Dempsey,
Zoltan Haiman,
Alessandro Lupi,
Michal Pirog,
Geoffrey Ryan
Abstract:
We have performed numerical calculations of a binary interacting with a gas disk, using eleven 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…
▽ More
We have performed numerical calculations of a binary interacting with a gas disk, using eleven 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 (a^2 Ω_B /ν)$ binary orbits, where $ν$ is the kinematic viscosity. We did not investigate dependence on binary mass ratio, eccentricity, disk temperature, or disk viscosity; therefore, these benchmarks may act as guides towards expanding converged solutions to the wider parameter space but might need to be updated in a future study that investigates dependence on system parameters. We find the most major discrepancies between codes resulted from the dimensionality of the setup (3D vs 2D disks). Beyond this, 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.
△ Less
Submitted 20 June, 2024; v1 submitted 20 February, 2024;
originally announced February 2024.
-
Analytical and Numerical Analysis of Circumbinary Disk Dynamics -- I: Coplanar Systems
Authors:
Siddharth Mahesh,
Sean T. McWilliams,
Michal Pirog
Abstract:
We present an analytical and numerical study of a system composed of a stellar binary pair and a massless, locally isothermal viscous accretion disk that is coplanar to the binary orbital plane. Analytically, we study the effect of the binary's gravitational potential over short timescales through the study of stability for epicyclic orbits, and over long timescales by revisiting the concept of re…
▽ More
We present an analytical and numerical study of a system composed of a stellar binary pair and a massless, locally isothermal viscous accretion disk that is coplanar to the binary orbital plane. Analytically, we study the effect of the binary's gravitational potential over short timescales through the study of stability for epicyclic orbits, and over long timescales by revisiting the concept of resonant torques. Numerically, we perform two-dimensional Newtonian numerical simulations of the disk-binary system over a range of binary mass ratios. We find that the results of our simulations are consistent with previous numerical studies. We additionally show, by comparison of the analytical and numerical results, that the circumbinary gap is maintained on the orbital timescale through the driving of epicyclic instabilities, and does not depend on resonant torquing, contrary to standard lore. While our results are applicable to any disk-binary system, we highlight the importance of this result in the search for electromagnetic and gravitational-wave signatures from supermassive black-hole binaries.
△ Less
Submitted 10 December, 2023; v1 submitted 2 May, 2023;
originally announced May 2023.
-
Self-gravitating axially symmetric disks in general-relativistic rotation
Authors:
Janusz Karkowski,
Wojciech Kulczycki,
Patryk Mach,
Edward Malec,
Andrzej Odrzywolek,
Michal Pirog
Abstract:
We integrate numerically axially symmetric stationary Einstein equations describing self-gravitating disks around spinless black holes. The numerical scheme is based on a method developed by Shibata, but contains important new ingredients. We derive a new general-relativistic Keplerian rotation law for self-gravitating disks around spinning black holes. Former results concerning rotation around sp…
▽ More
We integrate numerically axially symmetric stationary Einstein equations describing self-gravitating disks around spinless black holes. The numerical scheme is based on a method developed by Shibata, but contains important new ingredients. We derive a new general-relativistic Keplerian rotation law for self-gravitating disks around spinning black holes. Former results concerning rotation around spin-less black holes emerge in the limit of a vanishing spin parameter. These rotation curves might be used for the description of rotating stars, after appropriate modification around the symmetry axis. They can be applied to the description of compact torus--black hole configurations, including active galactic nuclei or products of coalescences of two neutron stars.
△ Less
Submitted 30 May, 2018; v1 submitted 8 February, 2018;
originally announced February 2018.
-
Geometry of Keplerian disk systems and bounds on masses of their components
Authors:
E. Malec,
P. Mach,
M. Pirog
Abstract:
We investigate accreting disk systems with polytropic gas in Keplerian motion. Numerical data and partial analytic results show that the self-gravitation of the disk speeds up its rotation -- its rotational frequency is larger than that given by the well known strictly Keplerian formula that takes into account the central mass only. Thus determination of central mass in systems with massive disks…
▽ More
We investigate accreting disk systems with polytropic gas in Keplerian motion. Numerical data and partial analytic results show that the self-gravitation of the disk speeds up its rotation -- its rotational frequency is larger than that given by the well known strictly Keplerian formula that takes into account the central mass only. Thus determination of central mass in systems with massive disks requires great care -- the strictly Keplerian formula yields only an upper bound. The effect of self-gravity depends on geometric aspects of disk configurations. Disk systems with a small (circa $10^{-4}$) ratio of the innermost radius to the outermost disk radius have the central mass close to the upper limit, but if this ratio is of the order of unity then the central mass can be smaller by many orders of magnitude from this bound.
△ Less
Submitted 22 November, 2012;
originally announced November 2012.
-
Estimating masses of Keplerian disk systems: the case of AGN in NGC 4258
Authors:
Patryk Mach,
Michal Pirog,
Edward Malec
Abstract:
The Keplerian motion of accretion disks in Active Galactic Nuclei (AGN) is usually believed to be generated by a heavy central mass. We investigate accreting disk systems with polytropic gas in Keplerian rotation and obtain a phenomenological formula that relates the Keplerian angular frequency to the ratio of disk and central masses. Central mass approaches the Keplerian value, if the inner bound…
▽ More
The Keplerian motion of accretion disks in Active Galactic Nuclei (AGN) is usually believed to be generated by a heavy central mass. We investigate accreting disk systems with polytropic gas in Keplerian rotation and obtain a phenomenological formula that relates the Keplerian angular frequency to the ratio of disk and central masses. Central mass approaches the Keplerian value, if the inner boundary of a disk is close to the minimal stable orbit of a black hole. These results are applied to NGC 4258, the unique AGN with a finely measured Keplerian rotation curve of the central disk, with the conclusion that its rotation curve is, in fact, determined by the central black hole. The mass of the accretion disk exceeds 100 solar masses.
△ Less
Submitted 30 October, 2012;
originally announced October 2012.