We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code
To better understand the orbital dynamics of exoplanets around close binary stars, i.e., circumbinary planets (CBPs), we applied techniques from dynamical systems theory to a physically motivated set of solutions in the Circular Restricted Three-Body Problem (CR3BP). We applied Floquet theory to characterize the linear dynamical behavior—static, oscillatory, or exponential—surrounding planar circumbinary periodic trajectories (limit cycles). We computed prograde and retrograde limit cycles and analyzed their geometries, stability bifurcations, and dynamical structures. Orbit and stability calculations are exact computations in the CR3BP and reproducible through the open-source Python package
- NSF-PAR ID:
- 10399921
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astronomical Journal
- Volume:
- 165
- Issue:
- 4
- ISSN:
- 0004-6256
- Format(s):
- Medium: X Size: Article No. 140
- Size(s):
- ["Article No. 140"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract Bothros to post-process data from 3D general relativistic magnetohydrodynamic simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller innermost stable circular orbit angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems. -
more » « less
Abstract Stellar-mass black holes (sBHs) embedded in gaseous disks of active galactic nuclei (AGN) can be important sources of detectable gravitational radiation for LIGO/Virgo when they form binaries and coalesce due to orbital decay. In this paper, we study the effect of dynamical friction (DF) on the formation of BH binaries in AGN disks using
N -body simulations. We employ two simplified models of DF, with the force on the BH depending on Δv v τ satisfies ΩKτ ≲ 20–30 (where ΩKis the Keplerian frequency around the SMBH): prograde binaries are formed when the DF is stronger (smallerτ ), while retrograde binaries dominate when the DF is weaker (largerτ ). We determine the distribution of both prograde and retrograde binaries as a function of initial orbital separation and the DF strength. Using our models of DF, we show that for a given sBH number density in the AGN disk, the formation rate of sBH binaries increases with decreasingτ and can reach a moderate value with a sufficiently strong DF. -
We present an approach based upon binary tree tensor network (BTTN) states for computing steady-state current statistics for a many-particle 1D ratchet subject to volume exclusion interactions. The ratcheted particles, which move on a lattice with periodic boundary conditions subject to a time-periodic drive, can be stochastically evolved in time to sample representative trajectories via a Gillespie method. In lieu of generating realizations of trajectories, a BTTN state can variationally approximate a distribution over the vast number of many-body configurations. We apply the density matrix renormalization group algorithm to initialize BTTN states, which are then propagated in time via the time-dependent variational principle (TDVP) algorithm to yield the steady-state behavior, including the effects of both typical and rare trajectories. The application of the methods to ratchet currents is highlighted, but the approach extends naturally to other interacting lattice models with time-dependent driving. Although trajectory sampling is conceptually and computationally simpler, we discuss situations for which the BTTN TDVP strategy can be beneficial.more » « less
-
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. -
more » « less
Abstract The HH 24 complex harbors five collimated jets emanating from a small protostellar multiple system. We have carried out a multiwavelength study of the jets, their driving sources, and the cloud core hosting the embedded stellar system, based on data from the Hubble Space Telescope, Gemini, Subaru, Apache Point Observatory 3.5 m, Karl G. Jansky Very Large Array, and Atacama Large Millimeter/submillimeter Array (ALMA) telescopes. The data show that the multiple system, SSV 63, contains at least 7 sources, ranging in mass from the hydrogen-burning limit to proto-Herbig Ae stars. The stars are in an unstable nonhierarchical configuration, and one member, a borderline brown dwarf, is moving away from the protostellar system with 25 km s−1, after being ejected ∼5800 yr ago as an orphaned protostar. Five of the embedded sources are surrounded by small, possibly truncated, disks resolved at 1.3 mm with ALMA. Proper motions and radial velocities imply jet speeds of 200–300 km s−1. The two main HH 24 jets, E and C, form a bipolar jet system that traces the innermost portions of parsec-scale chains of Herbig–Haro and H2shocks with a total extent of at least 3 pc. H2CO and C18O observations show that the core has been churned and continuously fed by an infalling streamer.13CO and12CO trace compact, low-velocity, cavity walls carved by the jets and an ultracompact molecular outflow from the most embedded object. Chaotic
N -body dynamics likely will eject several more of these objects. The ejection of stars from their feeding zones sets their masses. Dynamical decay of nonhierarchical systems can thus be a major contributor to establishing the initial mass function.
