No Arabic abstract
Electron-positron pair creation near sub-Eddington accretion rate black holes is believed to be dominated by the Breit-Wheeler process (photon-photon collisions). The interacting high energy photons are produced when unscreened electric fields accelerate leptons either in coherent, macroscopic gaps or in incoherent structures embedded in the turbulent plasma flow. The latter type of acceleration results in a drizzle of pair production sourced by photons from the background radiation field whose energies are near the pair-production threshold. In this work, we use radiation GRMHD simulations to extend an earlier study of pair drizzle by Moscibrodzka et al. We focus on low-magnetization (SANE) accretion onto supermassive Kerr black holes and consider radiation due to synchrotron, bremsstrahlung, and Compton upscattering processes. We confirm that pair drizzle in M87 is sufficient to keep the magnetospheric charge density orders of magnitude above the Goldreich-Julian density. We also find that pair production peaks along the jet-disk boundary.
We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5times 10^8M_{odot}$ black hole with accretion rates varying from $sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases with and without net vertical magnetic flux. The inner regions of all disks have radiation pressure $sim 10^4-10^6$ times the gas pressure. Non-axisymmetric density waves that steepen into spiral shocks form as gas flows towards the black hole. In simulations without net vertical magnetic flux, Reynolds stress generated by the spiral shocks are the dominant mechanism to transfer angular momentum. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when net vertical magnetic flux is sufficiently large. Outflows are formed with speed $sim 0.1-0.4c$. When the accretion rate is smaller than $sim 500 L_{Edd}/c^2$, outflows start around $10$ gravitational radii and the radiative efficiency is $sim 5%-7%$ with both magnetic field configurations. With accretion rate reaching $1500 L_{Edd}/c^2$, most of the funnel region close to the rotation axis becomes optically thick and the outflow only develops beyond $50$ gravitational radii. The radiative efficiency is reduced to $1%$. We always find the kinetic energy luminosity associated with the outflow is only $sim 15%-30%$ of the radiative luminosity. The mass flux lost in the outflow is $sim 15%-50%$ of the net mass accretion rates. We discuss implications of our simulation results on the observational properties of these disks.
We use global three dimensional radiation magneto-hydrodynamic simulations to study the properties of inner regions of accretion disks around a 5times 10^8 solar mass black hole with mass accretion rates reaching 7% and 20% of the Eddington value. This region of the disk is supported by magnetic pressure with surface density significantly smaller than the values predicted by the standard thin disk model but with a much larger disk scale height. The disks do not show any sign of thermal instability over many thermal time scales. More than half of the accretion is driven by radiation viscosity in the optically thin corona region for the lower accretion rate case, while accretion in the optically thick part of the disk is driven by the Maxwell and Reynolds stresses from MRI turbulence. Coronae with gas temperatures > 10^8 K are generated only in the inner approx 10 gravitational radii in both simulations, being more compact in the higher accretion rate case. In contrast to the thin disk model, surface density increases with increasing mass accretion rate, which causes less dissipation in the optically thin region and a relatively weaker corona. The simulation results may explain the formation of X-ray coronae in Active Galactic Nuclei (AGNs), the compact size of such coronae, and the observed trend of optical to X-ray luminosity with Eddington ratio for many AGNs.
Ultra-luminous X-ray sources (ULXs) have been puzzling us with a debate whether they consist of an intermediate mass black hole or super-Eddington accretion by a stellar mass black hole. Here we suggest that in the presence of large scale strong magnetic fields and non-negligible vertical motion, the luminosity of ULXs, particularly in their hard states, can be explained with sub-Eddington accretion by stellar mass black holes. In this framework of 2.5D magnetized advective accretion flows, magnetic tension plays the role of transporting matter (equivalent to viscous shear via turbulent viscosity) and we neither require to invoke an intermediate mass black hole nor super-Eddington accretion. Our model explains the sources, like, NGC 1365 X1/X2, M82 X42.3+59, M99 X1 etc. which are in their hard power-law dominated states.
The merger rate of stellar-mass black hole binaries (sBHBs) inferred by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) suggests the need for an efficient source of sBHB formation. Active galactic nucleus (AGN) disks are a promising location for the formation of these sBHBs, as well as binaries of other compact objects, because of powerful torques exerted by the gas disk. These gas torques cause orbiting compact objects to migrate towards regions in the disk where inward and outward torques cancel, known as migration traps. We simulate the migration of stellar mass black holes in an example of a model AGN disk, using an augmented N-body code that includes analytic approximations to migration torques, stochastic gravitational forces exerted by turbulent density fluctuations in the disk, and inclination and eccentricity dampening produced by passages through the gas disk, in addition to the standard gravitational forces between objects. We find that sBHBs form rapidly in our model disk as stellar-mass black holes migrate towards the migration trap. These sBHBs are likely to subsequently merge on short time-scales. The process continues, leading to the build-up of a population of over-massive stellar-mass black holes. The formation of sBHBs in AGN disks could contribute significantly to the sBHB merger rate inferred by LIGO.
We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$alpha$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $Gamma$, and the X-ray-to-optical power slopes, $alpha_{rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.