No Arabic abstract
Accreting supermassive binary black holes (SMBBHs) are potential multi-messenger sources because they emit both gravitational wave and electromagnetic (EM) radiation. Past work has shown that their EM output may be periodically modulated by an asymmetric density distribution in the circumbinary disk, often called an overdensity or lump; this modulation could possibly be used to identify a source as a binary. We explore the sensitivity of the overdensity to SMBBH mass ratio and magnetic flux through the accretion disk. We find that the relative amplitude of the overdensity and its associated EM periodic signal both degrade with diminishing mass ratio, vanishing altogether somewhere between 1:2 and 1:5. Greater magnetization also weakens the lump and any modulation of the light output. We develop a model to describe how lump formation results from internal stress degrading faster in the lump region than it can be rejuvenated through accretion inflow, and predicts a threshold value in specific internal stress below which lump formation should occur and which all our lump-forming simulations satisfy. Thus, detection of such a modulation would provide a constraint on both mass-ratio and magnetic flux piercing the accretion flow.
We carry out 2D viscous hydrodynamics simulations of circumbinary disk (CBD) accretion using {footnotesize AREPO}. We resolve the accretion flow from a large-scale CBD down to the streamers and disks around individual binary components. Extending our recent studies citep{mun19}, we consider circular binaries with various mass ratios ($0.1leq q_{rm{b}}leq1$) and study accretion from ``infinite, steady-supply disks and from finite-sized, viscously spreading tori. For ``infinite disks, a global steady state can be reached, and the accretion variability has a dominant frequency ${sim}0.2Omega_{rm{b}}$ for $q_{rm{b}}>0.5$ and $Omega_{rm{b}}$ for $q_{rm{b}}<0.5$, ($Omega_{rm{b}}$ is the binary angular frequency). We find that the accretion ``eigenvalue $l_0$ -- the net angular momentum transfer from the disk to the binary per unit accreted mass -- is always positive and falls in the range ($0.65$-$0.85)a_{rm b}^2Omega_{rm{b}}$ (with $a_{rm{b}}$ the binary separation), depending weakly on the mass ratio and viscosity. This leads to binary expansion when $q_{rm{b}}gtrsim0.3$. Accretion from a finite torus can be separated into two phases: an initial transient phase, corresponding to the filling of the binary cavity, followed by a viscous pseudo-stationary phase, during which the torus viscously spreads and accretes onto the binary. In the viscous phase, the net torque on the binary per unit accreted mass is close to $l_0$, the value derived for ``infinite disks. We conclude that similar-mass binaries accreting from CBDs gain angular momentum and expand over long time scales. This result significantly impacts the coalescence of supermassive binary black holes and newly formed binary stars. We offer a word of caution against conclusions drawn from simulations of transient accretion onto empty circumbinary cavities.
We investigate sustenance and dependence on magnetic Prandtl number (${rm Pm}$) for magnetorotational instability (MRI)-driven turbulence in astrophysical Keplerian disks with zero net magnetic flux using standard shearing box simulations. We focus on the turbulence dynamics in Fourier space, capturing specific/noncanonical anisotropy of nonlinear processes due to disk flow shear. This is a new type of nonlinear redistribution of modes over wavevector orientations in Fourier space -- the nonlinear transverse cascade -- which is generic to shear flows and fundamentally different from usual direct/inverse cascade. The zero flux MRI has no exponentially growing modes, so its growth is transient, or nonmodal. Turbulence self-sustenance is governed by constructive cooperation of the transient growth of MRI and the nonlinear transverse cascade. This cooperation takes place at small wavenumbers (on the flow size scales) referred to as the vital area in Fourier space. The direct cascade transfers mode energy from the vital area to larger wavenumbers. At large ${rm Pm}$, the transverse cascade prevails over the direct one, keeping most of modes energy contained in small wavenumbers. With decreasing ${rm Pm}$, however, the action of the transverse cascade weakens and can no longer oppose the action of direct cascade which more efficiently transfers energy to higher wavenumbers, leading to increased resistive dissipation. This undermines the sustenance scheme, resulting in the turbulence decay. Thus, the decay of zero net flux MRI-turbulence with decreasing ${rm Pm}$ is attributed to topological rearrangement of the nonlinear processes when the direct cascade begins to prevail over the transverse cascade.
High-mass multiples might form via fragmentation of self-gravitational disks or alternative scenarios such as disk-assisted capture. However, only few observational constraints exist on the architecture and disk structure of high-mass protobinaries and their accretion properties. Here we report the discovery of a close ($57.9pm0.2$mas=170au) high-mass protobinary, IRAS17216-3801, where our VLTI/GRAVITY+AMBER near-infrared interferometry allows us to image the circumstellar disks around the individual components with 3 milliarcsecond resolution. We estimate the component masses to $sim20$ and $sim18 M_{odot}$ and find that the radial intensity profiles can be reproduced with an irradiated disk model, where the inner regions are excavated of dust, likely tracing the dust sublimation region in these disks. The circumstellar disks are strongly misaligned with respect to the binary separation vector, which indicates that the tidal forces did not have time to realign the disks, pointing towards a young dynamical age of the system. We constrain the distribution of the Br$gamma$ and CO-emitting gas using VLTI/GRAVITY spectro-interferometry and VLT/CRIRES spectro-astrometry and find that the secondary is accreting at a higher rate than the primary. VLT/NACO imaging shows $L$-band emission on 3-4 times larger scales than the binary separation, matching the expected dynamical truncation radius for the circumbinary disk. The IRAS17216-3801 system is $sim3times$ more massive and $sim5times$ more compact than other high-mass multiplies imaged at infrared wavelengths and the first high-mass protobinary system where circumstellar and circumbinary dust disks could be spatially resolved. This opens exciting new opportunities for studying star-disk interactions and the role of multiplicity in high-mass star formation.
As 2 black holes bound to each other in a close binary approach merger their inspiral time becomes shorter than the characteristic inflow time of surrounding orbiting matter. Using an innovative technique in which we represent the changing spacetime in the region occupied by the orbiting matter with a 2.5PN approximation and the binary orbital evolution with 3.5PN, we have simulated the MHD evolution of a circumbinary disk surrounding an equal-mass non-spinning binary. Prior to the beginning of the inspiral, the structure of the circumbinary disk is predicted well by extrapolation from Newtonian results. The binary opens a low-density gap whose radius is roughly two binary separations, and matter piles up at the outer edge of this gap as inflow is retarded by torques exerted by the binary; nonetheless, the accretion rate is diminished relative to its value at larger radius by only about a factor of 2. During inspiral, the inner edge of the disk at first moves inward in coordination with the shrinking binary, but as the orbital evolution accelerates, the rate at which the inner edge moves toward smaller radii falls behind the rate of binary compression. In this stage, the rate of angular momentum transfer from the binary to the disk slows substantially, but the net accretion rate decreases by only 10-20%. When the binary separation is tens of gravitational radii, the rest-mass efficiency of disk radiation is a few percent, suggesting that supermassive binary black holes in galactic nuclei could be very luminous at this stage of their evolution. If the luminosity were optically thin, it would be modulated at a frequency that is a beat between the orbital frequency of the disks surface density maximum and the binary orbital frequency. However, a disk with sufficient surface density to be luminous should also be optically thick; as a result, the periodic modulation may be suppressed.
The magnetic fields of accretion disks play an important role in studying their evolution. We may assume that its generation is connected to the dynamo mechanism, which is similar with that in the galactic disks. Here, we propose a model of the magnetic field of the accretion disk that uses the same approaches that have been used for galaxies. It is necessary to obtain the field, which is expected to be less than the equipartition value, and without destroying the disk. To do so, it is necessary to formulate the basic properties of the ionized medium and to estimate the parameters governing the dynamo. We used the no-z approximation that has been developed for thin disks. We also take different boundary conditions that can change the value of the field significantly. We show that the magnetic field strictly depends on the boundary conditions. Taking zero conditions and the fixed magnetic field condition on the inner boundary, which are connected to the physical properties of the accretion disk, we can avoid solutions that are greater than the equipartition field.