We study the effects of general relativity (GR) on the evolution and alignment of circumbinary disks around binaries on all scales. We implement relativistic apsidal precession of the binary into the hydrodynamics code
We describe the first grid-based simulations of the polar alignment of a circumbinary disc. We simulate the evolution of an inclined disc around an eccentric binary using the grid-based code athena++ . The use of a grid-based numerical code allows us to explore lower disc viscosities than have been examined in previous studies. We find that the disc aligns to a polar orientation when the α viscosity is high, while discs with lower viscosity nodally precess with little alignment over 1000 binary orbital periods. The time-scales for polar alignment and disc precession are compared as a function of disc viscosity, and are found to be in agreement with previous studies. At very low disc viscosities (e.g. α = 10−5), anticyclonic vortices are observed along the inner edge of the disc. These vortices can persist for thousands of binary orbits, creating azimuthally localized overdensities and multiple pairs of spiral arms. The vortex is formed at ∼3–4 times the binary semimajor axis, close to the inner edge of the disc, and orbits at roughly the local Keplerian speed. The presence of a vortex in the disc may play an important role in the evolution of circumbinary systems, such as driving episodic accretion and accelerating the formation of polar circumbinary planets.
more » « less- Award ID(s):
- 1753168
- NSF-PAR ID:
- 10396215
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 520
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 2138-2147
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract phantom . We find that the effects of GR can suppress the stable polar alignment of a circumbinary disk, depending on how the relativistic binary apsidal precession timescale compares to the disk nodal precession timescale. Studies of circumbinary disk evolution typically ignore the effects of GR, which is an appropriate simplification for low-mass or widely separated binary systems. In this case, polar alignment occurs, provided that the disks initial misalignment is sufficiently large. However, systems with a very short relativistic precession timescale cannot polar align and instead move toward coplanar alignment. In the intermediate regime where the timescales are similar, the outcome depends upon the properties of the disk. Polar alignment is more likely in the wavelike disk regime (where the disk viscosity parameter is less than the aspect ratio,α <H /r ), since the disk is in good radial communication. In the viscous disk regime, disk breaking is more likely. Multiple rings can destructively interact with one another, resulting in short disk lifetimes and the disk moving toward coplanar alignment. Around main-sequence star or stellar mass black hole binaries, polar alignment may be suppressed far from the binary, but in general, the inner parts of the disk can align to polar. Polar alignment may be completely suppressed for disks around supermassive black holes for close binary separations. -
more » « less
ABSTRACT Hydrodynamical interaction in circumbinary discs (CBDs) plays a crucial role in various astrophysical systems, ranging from young stellar binaries to supermassive black hole binaries in galactic centres. Most previous simulations of binary-disc systems have adopted locally isothermal equation of state. In this study, we use the grid-based code Athena++ to conduct a suite of two-dimensional viscous hydrodynamical simulations of circumbinary accretion on a Cartesian grid, resolving the central cavity of the binary. The gas thermodynamics is treated by thermal relaxation towards an equilibrium temperature (based on the constant − β cooling ansatz, where β is the cooling time in units of the local Keplerian time). Focusing on equal mass, circular binaries in CBDs with (equilibrium) disc aspect ratio H/R = 0.1, we find that the cooling of the disc gas significantly influences the binary orbital evolution, accretion variability, and CBD morphology, and the effect depends sensitively on the disc viscosity prescriptions. When adopting a constant kinematic viscosity, a finite cooling time (β ≳ 0.1) leads to a binary inspiral as opposed to an outspiral and the CBD cavity becomes more symmetric. When adopting a dynamically varying α-viscosity, binary inspiral only occurs within a narrow range of cooling time (corresponding to β around 0.5).
-
more » « less
ABSTRACT A test particle orbit around an eccentric binary has two stationary states in which there is no nodal precession: coplanar and polar. Nodal precession of a misaligned test particle orbit centres on one of these stationary states. A low-mass circumbinary disc undergoes the same precession and moves towards one of these states through dissipation within the disc. For a massive particle orbit, the stationary polar alignment occurs at an inclination less than 90°, which is the prograde-polar stationary inclination. A sufficiently high angular momentum particle has an additional higher inclination stationary state, the retrograde-polar stationary inclination. Misaligned particle orbits close to the retrograde-polar stationary inclination are not nested like the orbits close to the other stationary points. We investigate the evolution of a gas disc that begins close to the retrograde-polar stationary inclination. With hydrodynamical disc simulations, we find that the disc moves through the unnested crescent shape precession orbits and eventually moves towards the prograde-polar stationary inclination, thus increasing the parameter space over which circumbinary discs move towards polar alignment. If protoplanetary discs form with an isotropic orientation relative to the binary orbit, then polar discs may be more common than coplanar discs around eccentric binaries, even for massive discs. This has implications for the alignment of circumbinary planets.
-
more » « less
ABSTRACT We present hydrodynamical simulations to model the accretion flow from a polar circumbinary disc on to a high eccentricity (e = 0.78) binary star system with near unity mass ratio (q = 0.83), as a model for binary HD 98800 BaBb. We compare the polar circumbinary disc accretion flow with the previously studied coplanar case. In the coplanar case, the circumbinary disc becomes eccentric and the accretion alternates from being dominant on to one binary member to the other. For the polar disc case involving a highly eccentric binary, we find that the circumbinary disc retains its initially low eccentricity and that the primary star accretion rate is always about the same as the secondary star accretion rate. Recent observations of the binary HD 98800 BaBb, which has a polar circumbinary disc, have been used to determine the value of the $\rm H\,\alpha$ flux from the brighter component. From this value, we infer that the accretion rate is much lower than for typical T Tauri stars. The eccentric orbit of the outer companion HD 98800 A increases the accretion rate on to HD 98800 B by ∼20 per cent after each periastron passage. Our hydrodynamical simulations are unable to explain such a low accretion rate unless the disc viscosity parameter is very small, α < 10−5. Additional observations of this system would be useful to check on this low accretion rate.
-
more » « less
Abstract Hydrodynamical interactions between binaries and circumbinary disks (CBDs) play an important role in a variety of astrophysical systems, from young stellar binaries to supermassive black hole binaries. Previous simulations of CBDs have mostly employed locally isothermal equations of state. We carry out 2D viscous hydrodynamic simulations of CBDs around equal-mass, circular binaries, treating the gas thermodynamics by thermal relaxation toward equilibrium temperature (the constant-
β cooling ansatz, whereβ is the cooling time in units of the local Keplerian time). As an initial study, we use the grid-based codeAthena++ on a polar grid, covering an extended disk outside the binary co-orbital region. We find that with a longer cooling time, the accretion variability is gradually suppressed, and the morphology of the CBD becomes more symmetric. The disk also shows evidence of hysteresis behavior depending on the initial conditions. Gas cooling also affects the rate of angular momentum transfer between the binary and the CBD, where given our adopted disk thickness and viscosity (H /r ∼ 0.1 andα ∼ 0.1), the binary orbit expands while undergoing accretion for mostβ values between 0 and 4.0 except over a narrow range of intermediateβ values. The validity of using a polar grid excising the central domain is also discussed.
