No Arabic abstract
The recent gravitational wave transient GW190521 has been interpreted by the LIGO-Virgo collaboration (LVC) as sourced by a binary black hole (BH) merger. According to the LVC parameter estimation, at least one of these progenitors falls into the so-called pair-instability supernova mass gap. This raises the important question of how and when these progenitors formed. In this paper we use an accretion model with super-Eddington mass accretion rate obtained from General Relativity hydrodynamics simulations to analyse the scenario wherein the GW190521 original progenitors (OPs) formed at lower masses (and spins) and grew to their estimated LVC parameters by relativistic accretion. We consider that the environment wherein the binary is immersed has density gradients as well as a dependence on the Mach number of the gas. Taking the LVC parameter estimation at $z=0.82$ as the endpoint of the accretion evolution, we estimate the initial masses and spins of the OPs at three different red-shifts $z=100, 50$, and $20$. We found three distinct possible types of OPs: $(i)$ $10^{-4} M_{odot} - 3 M_{odot}$ almost non-rotating (with Kerr spin parameter $a_{star}< 10^{-2}$) primordial BHs; $(ii)$ $3 M_{odot} - 40M_{odot}$ slowly rotating ($ 10^{-2} < a_{star} < 0.5$) stellar mass BHs; $(iii)$ $40M_{odot} - 70M_{odot}$ BHs with a moderate spin parameter $a_{star}sim 0.5$, which could originate from the collapse of high mass Pop III stars. The mass spread is due to varying the density gradient and the relativistic Mach number of the cosmic plasma; the variation of the masses due to the origin at different red-shifts, on the other hand, is negligible, $sim 2%$ ...
The gravitational-wave (GW) detection of GW190521 has provided new insights on the mass distribution of black holes and new constraints for astrophysical formation channels. With independent claims of GW190521 having significant pre-merger eccentricity, we investigate what this implies for GW190521-like binaries that form dynamically. The Laser Interferometer Space Antenna (LISA) will also be sensitive to GW190521-like binaries if they are circular from an isolated formation channel. We show, however, that GW190521-like binaries that form dynamically may skip the LISA band entirely. To this end, we simulate GW190521 analogues that dynamically form via post-Newtonian binary-single scattering. From these scattering experiments, we find that GW190521-like binaries may enter the LIGO-Virgo band with significant eccentricity as suggested by recent studies, though well below an eccentricity of $e_{rm 10Hz} lesssim 0.7$. Eccentric GW190521-like binaries further motivate the astrophysical science case for a decihertz GW observatory, such as the kilometer-scale version of the Midband Atomic Gravitational-wave Interferometric Sensor (MAGIS). Pre-merger observations of GW190521-like binaries with such a decihertz GW detector would be able to constrain the eccentricity of GW190521-like binaries to greater precision than with just LIGO-Virgo alone. These eccentricity constraints would also provide additional insights into the possible environments that GW190521-like binaries form in.
When an accretion disk falls prey to the runaway instability, a large portion of its mass is devoured by the black hole within a few dynamical times. Despite decades of effort, it is still unclear under what conditions such an instability can occur. The technically most advanced relativistic simulations to date were unable to find a clear sign for the onset of the instability. In this work, we present three-dimensional relativistic hydrodynamics simulations of accretion disks around black holes in dynamical space-time. We focus on the configurations that are expected to be particularly prone to the development of this instability. We demonstrate, for the first time, that the fully self-consistent general relativistic evolution does indeed produce a runaway instability.
We revisit Bondi accretion - steady-state, adiabatic, spherical gas flow onto a Schwarzschild black hole at rest in an asymptotically homogeneous medium - for stiff polytropic equations of state (EOSs) with adiabatic indices $Gamma > 5/3$. A general relativistic treatment is required to determine their accretion rates, for which we provide exact expressions. We discuss several qualitative differences between results for soft and stiff EOSs - including the appearance of a minimum steady-state accretion rate for EOSs with $Gamma geq 5/3$ - and explore limiting cases in order to examine these differences. As an example we highlight results for $Gamma = 2$, which is often used in numerical simulations to model the EOS of neutron stars. We also discuss a special case with this index, the ultra-relativistic `causal EOS, $P = rho$. The latter serves as a useful limit for the still undetermined neutron-star EOS above nuclear density. The results are useful, for example, to estimate the accretion rate onto a mini-black hole residing at the center of a neutron star.
The exact time-dependent solution is obtained for a magnetic field growth during a spherically symmetric accretion into a black hole (BH) with a Schwarzschild metric. Magnetic field is increasing with time, changing from the initially uniform into a quasi-radial field. Equipartition between magnetic and kinetic energies in the falling gas is established in the developed stages of the flow. Estimates of the synchrotron radiation intensity are presented for the stationary flow. The main part of the radiation is formed in the region $r leq 7 r_g$, here $r_g$ is a BH gravitational radius. The two-dimensional stationary self-similar magnetohydrodynamic solution is obtained for the matter accretion into BH, in a presence of a large-scale magnetic field, when the magnetic field far from the BH is homogeneous and does not influence the flow. At the symmetry plane perpendicular to the direction of the distant magnetic field, the quasi-stationary disk is formed around BH, which structure is determined by dissipation processes. Parameters of the shock forming due to matter infall onto the disk are obtained. The radiation spectrum of the disk and the shock are obtained for the $10,, M_odot$ BH. The luminosity of such object is about the solar one, for a characteristic galactic gas density, with possibility of observation at distances less than 1 kpc. The spectra of a laminar and a turbulent disk structure around BH are very different. The turbulent disk emits a large part of its flux in the infrared. It may occur that some of the galactic infrared star-like sources are a single BH in the turbulent accretion state. The radiative efficiency of the magnetized disk is very high, reaching $sim 0.5,dot M,c^2$ so it was called recently as a magnetically arrested disk (MAD). Numerical simulations of MAD, and its appearance during accretion into neutron stars are considered and discussed.
A space-based interferometer such as eLISA could observe few to few thousands progenitors of black hole binaries (BHBs) similar to those recently detected by Advanced LIGO. Gravitational radiation circularizes the orbit during inspiral, but some BHBs retain a measurable eccentricity at the low frequencies where eLISA is most sensitive. The eccentricity of a BHB carries precious information about its formation channel: BHBs formed in the field, in globular clusters, or close to a massive black hole (MBH) have distinct eccentricity distributions in the eLISA band. We generate mock eLISA observations, folding in measurement errors, and using Bayesian model selection we study whether eLISA measurements can identify the BHB formation channel. We find that a handful of observations would suffice to tell whether BHBs were formed in the gravitational field of a MBH. Conversely, several tens of observations are needed to tell apart field formation from globular cluster formation. A five-year eLISA mission with the longest possible armlength is desirable to shed light on BHB formation scenarios.