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.


Title: Grid-based simulations of polar circumbinary discs: polar alignment and vortex formation
ABSTRACT 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
PAR ID:
10396215
Author(s) / Creator(s):
; ; ;
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
  1. Abstract 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 codephantom. 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
  2. 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
  3. Abstract Mutually misaligned circumbinary planets may form in a warped or broken gas disk or from later planet–planet interactions. With numerical simulations and analytic estimates we explore the dynamics of two circumbinary planets with a large mutual inclination. A coplanar inner planet causes prograde apsidal precession of the binary and the stationary inclination for the outer planet is higher for larger outer planet orbital radius. In this case a coplanar outer planet always remains coplanar. On the other hand, a polar inner planet causes retrograde apsidal precession of the binary orbit and the stationary inclination is smaller for larger outer planet orbital radius. For a range of outer planet semimajor axes, an initially coplanar orbit is librating meaning that the outer planet undergoes large tilt oscillations. Circumbinary planets that are highly inclined to the binary are difficult to detect—it is unlikely for a planet to have an inclination below the transit detection limit in the presence of a polar inner planet. These results suggest that there could be a population of circumbinary planets that are undergoing large tilt oscillations. 
    more » « less
  4. Vortex crystals are quasiregular arrays of like-signed vortices in solid-body rotation embedded within a uniform background of weaker vorticity. Vortex crystals are observed at the poles of Jupiter and in laboratory experiments with magnetized electron plasmas in axisymmetric geometries. We show that vortex crystals form from the free evolution of randomly excited two-dimensional turbulence on an idealized polar cap. Once formed, the crystals are long lived and survive until the end of the simulations (300 crystal-rotation periods). We identify a fundamental length scale, L γ = ( U / γ ) 1 / 3 , characterizing the size of the crystal in terms of the mean-square velocity U of the fluid and the polar parameter γ = f p / a p 2 , with f p the Coriolis parameter at the pole and a p the polar radius of the planet. 
    more » « less
  5. null (Ed.)
    Abstract Among the potential milliHz gravitational wave (GW) sources for the upcoming space-based interferometer LISA are extreme- or intermediate-mass ratio inspirals (EMRI/IMRIs). These events involve the coalescence of supermassive black holes in the mass range 105M⊙ ≲ M ≲ 107M⊙ with companion BHs of much lower masses. A subset of E/IMRIs are expected to occur in the accretion discs of active galactic nuclei (AGN), where torques exerted by the disc can interfere with the inspiral and cause a phase shift in the GW waveform. Here we use a suite of two-dimensional hydrodynamical simulations with the moving-mesh code DISCO to present a systematic study of disc torques. We measure torques on an inspiraling BH and compute the corresponding waveform deviations as a function of the binary mass ratio q ≡ M2/M1, the disc viscosity (α), and gas temperature (or equivalently Mach number; $$\mathcal {M}$$). We find that the absolute value of the gas torques is within an order of magnitude of previously determined planetary migration torques, but their precise value and sign depends non-trivially on the combination of these parameters. The gas imprint is detectable by LISA for binaries embedded in AGN discs with surface densities above $$\Sigma _0\ge 10^{4-6} \rm \, g cm^{-2}$$, depending on q, α and $$\mathcal {M}$$. Deviations are most pronounced in discs with higher viscosities, and for E/IMRIs detected at frequencies where LISA is most sensitive. Torques in colder discs exhibit a noticeable dependence on the GW-driven inspiral rate as well as strong fluctuations at late stages of the inspiral. Our results further suggest that LISA may be able to place constraints on AGN disc parameters and the physics of disc-satellite interaction. 
    more » « less