Blob formation and ejection from the radiative inefficient accretion flow around massive black hole

Aims. It is speculated that a merger of two massive stellar-origin black holes in a dense stellar environment may lead to the formation of a massive black hole in the pair-instability mass gap (∼50−135 M⊙). Such a merger-formed black hole is expected to typically have a high spin (a ∼ 0.7). If such a massive black hole acquires another black hole it may lead to another merger detectable by LIGO/Virgo in gravitational waves. Acquiring a companion may be hindered by gravitational-wave kick/recoil, which accompanies the first merger and may quickly remove the massive black hole from its parent globular or nuclear cluster. We test whether it is possible for a massive merger-formed black hole in the pair-instability gap to be retained in its parent cluster and have low spin. Such a black hole would be indistinguishable from a primordial black hole. Methods. We employed results from numerical relativity calculations of black hole mergers to explore the range of gravitational-wave recoil velocities for various combinations of merging black hole masses and spins. We compared merger-formed massive black hole speeds with typical escape velocities from globular and nuclear clusters. Results. We show that a globular cluster is highly unlikely to form and retain a ∼100 M⊙ black hole if the spin of the black hole is low (a ≲ 0.3). Massive merger-formed black holes with low spins acquire high recoil speeds (≳ 200 km s−1) from gravitational-wave kick during formation that exceed typical escape speeds from globular clusters (∼ 50 km s−1). However, a very low-spinning (a ∼ 0.1) and massive (∼100 M⊙) black hole could be formed and retained in a galactic nuclear star cluster. Even though such massive merger-formed black holes with such low spins acquire high speeds during formation (∼ 400 km s−1), they may avoid ejection since massive nuclear clusters have high escape velocities (∼ 300−500 km s−1). A future detection of a massive black hole in the pair-instability mass gap with low spin would therefore not be proof of the existence of primordial black holes, which are sometimes claimed to have low spins and arbitrarily high masses.

Black holes that accrete far below the Eddington limit are believed to accrete through a geometrically thick, optically thin, rotationally supported plasma that we will refer to as a radiatively inefficient accretion flow (RIAF). RIAFs are typically collisionless in the sense that the Coulomb mean free path is large compared to $GM/c^2$, and relativistically hot near the event horizon. In this paper we develop a phenomenological model for the plasma in RIAFs, motivated by the application to sources such as Sgr A* and M87. The model is derived using Israel-Stewart theory, which considers deviations up to second order from thermal equilibrium, but modified for a magnetized plasma. This leads to thermal conduction along magnetic field lines and a difference in pressure, parallel and perpendicular to the field lines (which is equivalent to anisotrotropic viscosity). In the non-relativistic limit, our model reduces to the widely used Braginskii theory of magnetized, weakly collisional plasmas. We compare our model to the existing literature on dissipative relativistic fluids, describe the linear theory of the plasma, and elucidate the physical meaning of the free parameters in the model. We also describe limits of the model when the conduction is saturated and when the viscosity implies a large pressure anisotropy. In future work, the formalism developed in this paper will be used in numerical models of RIAFs to assess the importance of non-ideal processes for the dynamics and radiative properties of slowly accreting black holes.

Using general relativistic radiative transfer (GRRT) simulations, we investigate the bright ring features and polarization structures in images of the Kerr-Sen black hole associated with Sgr A*, as illuminated by 230 GHz thermal synchrotron emission from radiatively inefficient accretion flows (RIAF). Our findings reveal that an increase in the dilaton parameter leads to a shrinking of the bright ring, accompanied by enhancements in both its width and brightness. As the disk thickness grows, the bright ring's diameter and width both decrease. The brightness enhancement induced by the disk thickness is less prominent than that driven by the dilaton parameter. Comparing with the Event Horizon Telescope (EHT) observational data of SgrA*, we present the allowed ranges of black hole parameters, and find that effects of the disk thickness on the allowed parameter space are stronger than those of the observer's inclination. Furthermore, we analyze the coefficient β 2, associated with the two-fold rotational symmetry of the electric vector position angles (EVPA), to probe the polarization structure of the black hole images, and reveal that effects of the disk thickness on β 2 are much weaker than those from the dilaton parameter.

In this pilot study, we investigate the use of a deep learning (DL) model to temporally evolve the dynamics of gas accreting on to a black hole in the form of a radiatively inefficient accretion flow (RIAF). We have trained a convolutional neural network (CNN) on a data set that consists of numerical solutions of the hydrodynamical equations for a range of initial conditions. We find that deep neural networks trained on one simulation seem to learn reasonably well the spatiotemporal distribution of densities and mass continuity of a black hole accretion flow over a duration of 8 × 104GM/c3, comparable to the viscous time-scale at r = 400GM/c2; after that duration, the model drifts from the ground truth suffering from excessive artificial mass injection. Models trained on simulations with different initial conditions show some promise of generalizing to configurations not present in the training set, but also suffer from mass continuity issues. We discuss the caveats behind this method and the potential benefits that DL models offer. For instance, once trained the model evolves an RIAF on a single GPU four orders of magnitude faster than usual fluid dynamics integrators running in parallel on 200 CPU cores. We speculate that a data-driven machine learning approach should be very promising for accelerating simulations of accreting black holes.

Recent work has suggested that black holes (BHs) could be cosmologically coupled to the accelerated expansion of the Universe, potentially becoming a candidate for origin of dark energy. This would imply that BH mass growth follows cosmological expansion, with the masses of individual BHs growing as MBH ∝ (1 + z)3. In this Letter, we discuss the binary systems Gaia BH1 and Gaia BH2, which contain ∼9 M⊙ BHs orbited by ∼1 M⊙ stars in widely separated orbits. The ages of both systems can be constrained by the properties of the luminous stars. If BH masses are indeed growing as (1 + z)3, the masses of both BHs at formation would have been significantly smaller than today. We find a 77% probability that the mass of the BH in Gaia BH2 would have been below 2.2 M⊙ at formation. This is below the classical Tolman–Oppenheimer–Volkov limit, though it is not yet clear whether or not BHs subject to cosmological coupling should obey this limit. For Gaia BH1, the same probability is 70%. This analysis is consistent with results from two BHs in the globular cluster NGC 3201, but unlike the NGC 3201 BHs, the Gaia BHs have well-constrained inclinations and therefore firm upper mass limits. The discovery of more BHs in binary systems with Gaia astrometry in the coming years will allow us to test the cosmological coupling hypothesis decisively.

The cosmic-ray (CR) accelerator at the Galactic center (GC) is not yet established by current observations. Here we investigate the radiative-inefficient accretion flow (RIAF) of Sagittarius A* (SgrA*) as a CR accelerator assuming acceleration by turbulent magnetic reconnection, and derive possible emission fluxes of CRs interacting within the RIAF (the central ∼1013 cm). The target environment of the RIAF is modeled with numerical, general relativistic magnetohydrodynamics together with leptonic radiative transfer simulations. The acceleration of the CRs is not computed here. Instead, we inject CRs constrained by the magnetic reconnection power of the accretion flow and compute the emission/absorption of γ-rays due to these CRs interacting with the RIAF, through Monte Carlo simulations employing the CRPropa 3 code. The resulting very-high-energy (VHE) fluxes are not expected to reproduce the point source HESS J1745-290 as the emission of this source is most likely produced at parsec scales. The emission profiles derived here intend to trace the VHE signatures of the RIAF as a CR accelerator and provide predictions for observations of the GC with improved angular resolution and differential flux sensitivity as those of the forthcoming Cerenkov Telescope Array (CTA). Within the scenario presented here, we find that for mass accretion rates ≳10−7 M ⊙yr−1, the RIAF of SgrA* produces VHE fluxes that are consistent with the High Energy Stereoscopic System (H.E.S.S.) upper limits for the GC and potentially observable by the future CTA. The associated neutrino fluxes are negligible compared with the diffuse neutrino emission measured by the IceCube.

Most massive galaxies have supermassive black holes at their centres, and the masses of the black holes are believed to correlate with properties of the host-galaxy bulge component. Several explanations have been proposed for the existence of these locally established empirical relationships, including the non-causal, statistical process of galaxy-galaxy merging, direct feedback between the black hole and its host galaxy, and galaxy-galaxy merging and the subsequent violent relaxation and dissipation. The empirical scaling relations are therefore important for distinguishing between various theoretical models of galaxy evolution, and they furthermore form the basis for all black-hole mass measurements at large distances. Observations have shown that the mass of the black hole is typically 0.1 per cent of the mass of the stellar bulge of the galaxy. Until now, the galaxy with the largest known fraction of its mass in its central black hole (11 per cent) was the small galaxy NGC 4486B. Here we report observations of the stellar kinematics of NGC 1277, which is a compact, lenticular galaxy with a mass of 1.2 × 10(11) solar masses. From the data, we determine that the mass of the central black hole is 1.7 × 10(10) solar masses, or 59 per cent of its bulge mass. We also show observations of five other compact galaxies that have properties similar to NGC 1277 and therefore may also contain over-massive black holes. It is not yet known if these galaxies represent a tail of a distribution, or if disk-dominated galaxies fail to follow the usual black-hole mass scaling relations.

We present results on the evolution of the intrinsic scatter of black hole masses considering different implementations of a model in which black holes only grow via mergers. We demonstrate how merger driven growth affects the correlations between black hole mass and host bulge mass. The simple case of an initially lognormal distributed scatter in black hole and bulge masses combined with random merging within the galaxy population results in a decreasing scatter with merging generation/number as predicted by the central-limit theorem. In general we find that the decrease in scatter σ is well approximated by σ<inf>merg</inf>(<it>m</it>) ∼σ<inf>ini</inf>× (<it>m</it>+ 1)−<it>a</it>/2 with <it>a</it>= 0.42 for a range of mean number of mergers <it>m</it> < 50. For a large mean number of mergers (<it>m</it> > 100) we find a convergence to <it>a</it>= 0.61. This is valid for a wide range of different initial distributions, refill-scenarios or merger mass ratios. Growth scenarios based on halo merger trees of a (100 Mpc)3ΛCDM-simulation show a similar behaviour with a scatter decrease of <it>a</it>= 0.30 with typical number of mergers <it>m</it> < 50 consistent with random merging (best matching model: <it>a</it>= 0.34). Assuming a present-day scatter of 0.3 dex in black hole mass and a mean number of mergers not exceeding <it>m</it>= 50 our results imply a scatter of 0.6 dex at <it>z</it>= 3 and thus a possible scenario in which overmassive (and undermassive) black holes at high redshift are a consequence of a larger intrinsic scatter in black hole mass. A simple toy model connecting the growth of black holes to the growth of ΛCDM haloes via mergers, neglecting any contribution from accretion, yields a consistent <it>M</it><inf>•</inf>–<it>M</it><inf>Bulge</inf> relation at <it>z</it>= 0– if we assume the correct initial relation.

QUANTUM gravitational effects are usually ignored in calculations of the formation and evolution of black holes. The justification for this is that the radius of curvature of space-time outside the event horizon is very large compared to the Planck length (Għ/c3)1/2 ≈ 10−33 cm, the length scale on which quantum fluctuations of the metric are expected to be of order unity. This means that the energy density of particles created by the gravitational field is small compared to the space-time curvature. Even though quantum effects may be small locally, they may still, however, add up to produce a significant effect over the lifetime of the Universe ≈ 1017 s which is very long compared to the Planck time ≈ 10−43 s. The purpose of this letter is to show that this indeed may be the case: it seems that any black hole will create and emit particles such as neutrinos or photons at just the rate that one would expect if the black hole was a body with a temperature of (κ/2π) (ħ/2k) ≈ 10−6 (M/M)K where κ is the surface gravity of the black hole1. As a black hole emits this thermal radiation one would expect it to lose mass. This in turn would increase the surface gravity and so increase the rate of emission. The black hole would therefore have a finite life of the order of 1071 (M/M)−3 s. For a black hole of solar mass this is much longer than the age of the Universe. There might, however, be much smaller black holes which were formed by fluctuations in the early Universe2. Any such black hole of mass less than 1015 g would have evaporated by now. Near the end of its life the rate of emission would be very high and about 1030 erg would be released in the last 0.1 s. This is a fairly small explosion by astronomical standards but it is equivalent to about 1 million 1 Mton hydrogen bombs. It is often said that nothing can escape from a black hole. But in 1974, Stephen Hawking realized that, owing to quantum effects, black holes should emit particles with a thermal distribution of energies — as if the black hole had a temperature inversely proportional to its mass. In addition to putting black-hole thermodynamics on a firmer footing, this discovery led Hawking to postulate 'black hole explosions', as primordial black holes end their lives in an accelerating release of energy.

The low-luminosity black hole Sgr A* provides a test bed for models of radiatively inefficient accretion flows (RIAFs). Recent submillimeter linear polarization measurements of Sgr A* have provided evidence that the electrons in the accretion flow are relativistic over a large range of radii. Here we show that these high temperatures result in elliptical plasma normal modes. Thus, polarized millimeter and submillimeter radiation emitted within RIAFs will undergo generalized Faraday rotation, a cyclic conversion between linear and circular polarization. This effect will not depolarize the radiation even if the rotation measure is extremely high. Rather, the beam will take on the linear and circular polarization properties of the plasma normal modes. As a result, polarization measurements of Sgr A* in this frequency regime will constrain the temperature, density, and magnetic profiles of RIAF models.

A significant fraction of the stellar population in the cusp around central\nblack holes of galaxies consists of compact remnants of evolved stars, such as\nwhite dwarfs, neutron stars and stellar mass black holes. We estimate the rate\nof capture of compact objects by massive central black holes, assuming most\nspiral galaxies have a central black hole of modest mass ($\\sim 10^6 \\msun$),\nand a cuspy spheroid. It is likely that the total capture rate is dominated by\nnucleated spirals. We estimate the flux of gravitational wave radiation from\nsuch coalescences, and the estimated detectable source count for proposed\nspace--based gravitational wave observatories such as LISA. About one event per\nyear should be detectable within $1 \\, {\\rm Gpc}$, given very conservative\nestimates of the black hole masses and central galactic densities. We expect\n$10^2$--$10^3$ detectable sources at lower frequencies ($10^{-4}$ Hz) ``en\nroute'' to capture. If stellar mass black holes are ubiquitous, the signal may\nbe dominated by stellar mass black holes coalescing with massive black holes.\nThe rate of white dwarf--white dwarf mergers in the cores of nucleated spirals\nis estimated at $\\sim10^{-6}$ per year per galaxy.\n

We study high-energy neutrino and cosmic-ray (CR) emission from the cores of low-luminosity active galactic nuclei (LLAGN). In LLAGN, the thermalization of particles is expected to be incomplete in radiatively inefficient accretion flows (RIAFs), allowing the existence of non-thermal particles. In this work, assuming stochastic particle acceleration due to turbulence in RIAFs, we solve the Fokker-Planck equation and calculate spectra of escaping neutrinos and CRs. The RIAF in LLAGN can emit CR protons with $\gtrsim10$ PeV energies and TeV-PeV neutrinos generated via $pp$ and/or $p\gamma$ reactions. We find that, if $\sim1$\% of the accretion luminosity is carried away by non-thermal ions, the diffuse neutrino intensity from the cores of LLAGN may be as high as $E_\nu^2\Phi_\nu\sim3\times{10}^{-8} {\rm GeV} {\rm cm}^{-2} {\rm s}^{-1}{\rm sr}^{-1}$, which can be compatible with the observed IceCube data. This result does not contradict either of the diffuse gamma-ray background observed by {\it Fermi} or observed diffuse cosmic-ray flux. Our model suggests that, although very-high-energy gamma rays may not escape, radio-quiet AGN with RIAFs can emit GeV gamma-rays, which could be used for testing the model. We also calculate the neutron luminosity from RIAFs of LLAGN, and discuss a strong constraint on the model of jet mass loading mediated by neutrons from the diffuse neutrino observation.

We investigate a mechanism for a super-massive black hole at the center of a galaxy to wander in the nucleus region. A situation is supposed in which the central black hole tends to move by the gravitational attractions from the nearby molecular clouds in a nuclear bulge but is braked via the dynamical frictions from the ambient stars there. We estimate the approximate kinetic energy of the black hole in an equilibrium between the energy gain rate through the gravitational attractions and the energy loss rate through the dynamical frictions in a nuclear bulge composed of a nuclear stellar disk and a nuclear stellar cluster as observed from our Galaxy. The wandering distance of the black hole in the gravitational potential of the nuclear bulge is evaluated to get as large as several 10 pc, when the black hole mass is relatively small. The distance, however, shrinks as the black hole mass increases, and the equilibrium solution between the energy gain and loss disappears when the black hole mass exceeds an upper limit. As a result, we can expect the following scenario for the evolution of the black hole mass: When the black hole mass is smaller than the upper limit, mass accretion of the interstellar matter in the circumnuclear region, causing the AGN activities, makes the black hole mass larger. However, when the mass gets to the upper limit, the black hole loses the balancing force against the dynamical friction and starts spiraling downward to the gravity center. From simple parameter scaling, the upper mass limit of the black hole is found to be proportional to the bulge mass, and this could explain the observed correlation of the black hole mass with the bulge mass.

The discovery of gravitational-wave radiation from merging black holes (BHs) also uncovered BHs with masses in the range of ≈20–160 M ⊙. In contrast, the most massive Galactic stellar-mass BH currently known has a mass of ≈21 M ⊙. While low-mass X-ray binaries (LMXBs) will never independently evolve into a binary BH system, and binary evolution effects can play an important role in explaining the different BH masses found through studies of X-ray binaries and gravitational-wave events, (electromagnetic) selection effects may also play a role in this discrepancy. Assuming BH LMXBs originate in the Galactic plane, we show that the spatial distributions of the current samples of confirmed and candidate BH LMXBs are both biased to sources that lie at a large distance from the plane. Specifically, most of the confirmed and candidate BH LMXBs are found at a Galactic height larger than three times the scale height for massive star formation. In addition, the confirmed BH LMXBs are found at larger distances to the Galactic center than the candidate BH LMXBs. Interstellar absorption makes candidate BH LMXBs in the plane and bulge too faint for a dynamical mass measurement using current instrumentation. Given the observed and theoretical evidence for BH natal and/or Blaauw kicks, their relation with BH mass and binary orbital period, and the relation between outburst recurrence time and BH mass, the observational selection effects imply that the current sample of confirmed BH LMXBs is biased against the most massive BHs.

We compile black hole (BH) masses for ~60,000 quasars in the redshift range 0.1 z 4.5 included in the Fifth Data Release of the Sloan Digital Sky Survey, using virial BH mass estimators based on the Hβ, Mg II, and C IV emission lines. Within our sample, the widths of the three lines follow lognormal distributions, with means and dispersions that do not depend strongly on luminosity or redshift. The Mg II- and Hβ-estimated BH masses are consistent with one another, but there is a positive bias between the C IV- and Mg II-estimated BH masses correlated with the C IV-Mg II blueshift, suggesting that the C IV estimator is more severely affected by a disk wind. If the underlying BH mass distribution decreases with mass and the Eddington ratio distribution at fixed true BH mass has nonzero width, we show that the measured virial BH mass and Eddington ratio distributions within finite luminosity bins are subject to Malmquist bias. We present a model that reproduces the observed virial mass distribution, quasar luminosity function, and line width distribution of our sample; it has an underlying BH mass distribution dN/dlog M ∝ M−2.6 and a lognormal true Eddington ratio distribution at fixed true mass with dispersion 0.4 dex and mean dependent on BH mass. In this model, the observed virial mass (Eddington ratio) distribution for the SDSS sample is biased high (low) by ~0.6 dex within finite luminosity bins. Finally, we compare virial BH masses of radio and broad absorption line quasars with ordinary quasars matched in redshift and luminosity.