Covariant Compton Scattering Kernel in General Relativistic Radiative Transfer

A covariant scattering kernel is a core component in any self-consistent general relativistic radiative transfer formulation in scattering media. An explicit closed-form expression for a covariant Compton scattering kernel with a good dynamical energy range has unfortunately not been available thus far. Such an expression is essential to obtain numerical solutions to the general relativistic radiative transfer equations in complicated astrophysical settings where strong scattering effects are coupled with highly relativistic flows and steep gravitational gradients. Moreover, this must be performed in an efficient manner. With a self-consistent covariant approach, we have derived a closed-form expression for the Compton scattering kernel for arbitrary energy range. The scattering kernel and its angular moments are expressed in terms of hypergeometric functions, and their derivations are shown explicitly in this paper. We also evaluate the kernel and its moments numerically, assessing various techniques for their calculation. Finally, we demonstrate that our closed-form expression produces the same results as previous calculations, which employ fully numerical computation methods and are applicable only in more restrictive settings.Comment: 29 pages, 10 figures, 2 tables; Accepted for publication in MNRA

Variations in emission from episodic plasmoid ejecta around black holes

The X-ray and radio flares observed in X-ray binaries and active galactic nuclei (AGN) are attributed to energetic electrons in the plasma ejecta from the accretion flows near the black hole in these systems. It is argued that magnetic reconnection could occur in the coronae above the accretion disk around the black hole, and that this drives plasmoid outflows resembling the solar coronal mass ejection (CME) phenomenon. The X-ray and radio flares are emission from energetic electrons produced in the process. As the emission region is located near the black hole event horizon, the flare emission would be subject to special- and general-relativistic effects. We present calculations of the flaring emission from plasmoids orbiting around a black hole and plasmoid ejecta launched from the inner accretion disk when general-relativistic effects are crucial in determining the observed time-dependent properties of the emission. We consider fully general-relativistic radiative transfer calculations of the emission from evolving ejecta from black hole systems, with proper accounting for differential arrival times of photons emitted from the plasmoids, and determine the emission lightcurves of plasmoids when they are in orbit and when they break free from their magnetic confinement. The implications for interpreting time-dependent spectroscopic observations of flaring emission from accreting black holes are discussed.Comment: 18 pages, 15 figures; Accepted for publication in MNRA

Dark matter concentrations in galactic nuclei according to polytropic models

We calculate the radial profiles of galaxies where the nuclear region is self-gravitating, consisting of self-interacting dark matter (SIDM) with $F$ degrees of freedom. For sufficiently high density this dark matter becomes collisional, regardless of its behaviour on galaxy scales. Our calculations show a spike in the central density profile, with properties determined by the dark matter microphysics, and the densities can reach the `mean density' of a black hole (from dividing the black-hole mass by the volume enclosed by the Schwarzschild radius). For a galaxy halo of given compactness ($\chi=2GM/Rc^2$), certain values for the dark matter entropy yield a dense central object lacking an event horizon. For some soft equations of state of the SIDM (e.g. $F\ge6$), there are multiple horizonless solutions at given compactness. Although light propagates around and through a sphere composed of dark matter, it is gravitationally lensed and redshifted. While some calculations give non-singular solutions, others yield solutions with a central singularity. In all cases the density transitions smoothly from the central body to the dark-matter envelope around it, and to the galaxy's dark matter halo. We propose that pulsar timing observations will be able to distinguish between systems with a centrally dense dark matter sphere (for different equations of state) and conventional galactic nuclei that harbour a supermassive black hole.Comment: MNRAS accepted, 24 pages, 12 figure

Simulations of recoiling black holes: adaptive mesh refinement and radiative transfer

(Abridged) We here continue our effort to model the behaviour of matter when orbiting or accreting onto a generic black hole by developing a new numerical code employing advanced techniques geared solve the equations of in general-relativistic hydrodynamics. The new code employs a number of high-resolution shock-capturing Riemann-solvers and reconstruction algorithms, exploiting the enhanced accuracy and the reduced computational cost of AMR techniques. In addition, the code makes use of sophisticated ray-tracing libraries that, coupled with general-relativistic radiation-transfer calculations, allow us to compute accurately the electromagnetic emissions from such accretion flows. We validate the new code by presenting an extensive series of stationary accretion flows either in spherical or axial symmetry and performed either in 2D or 3D. In addition, we consider the highly nonlinear scenario of a recoiling black hole produced in the merger of a supermassive black hole binary interacting with the surrounding circumbinary disc. In this way we can present, for the first time, ray-traced images of the shocked fluid and the light-curve resulting from consistent general-relativistic radiation-transport calculations from this process. The work presented here lays the ground for the development of a generic computational infrastructure employing AMR techniques to deal accurately and self-consistently with accretion flows onto compact objects. In addition to the accurate handling of the matter, we provide a self-consistent electromagnetic emission from these scenarios by solving the associated radiative-transfer problem. While magnetic fields are presently excluded from our analysis, the tools presented here can have a number of applications to study accretion flows onto black holes or neutron stars.Comment: 20 pages, 20 figures, accepted for publication in A&

Observable Emission Features of Black Hole GRMHD Jets on Event Horizon Scales

The general-relativistic magnetohydrodynamical (GRMHD) formulation for black hole-powered jets naturally gives rise to a stagnation surface, wherefrom inflows and outflows along magnetic field lines that thread the black hole event horizon originate. We derive a conservative formulation for the transport of energetic electrons which are initially injected at the stagnation surface and subsequently transported along flow streamlines. With this formulation the energy spectra evolution of the electrons along the flow in the presence of radiative and adiabatic cooling is determined. For flows regulated by synchrotron radiative losses and adiabatic cooling, the effective radio emission region is found to be finite, and geometrically it is more extended along the jet central axis. Moreover, the emission from regions adjacent to the stagnation surface is expected to be the most luminous as this is where the freshly injected energetic electrons concentrate. An observable stagnation surface is thus a strong prediction of the GRMHD jet model with the prescribed non-thermal electron injection. Future millimeter/sub-millimeter (mm/sub-mm) very-long-baseline interferometric (VLBI) observations of supermassive black hole candidates, such as the one at the center of M87, can verify this GRMHD jet model and its associated non-thermal electron injection mechanism.Comment: 19 pages, 12 figure; accepted for publication in Ap

The Current Ability to Test Theories of Gravity with Black Hole Shadows

Our Galactic Center, Sagittarius A* (Sgr A*), is believed to harbour a supermassive black hole (BH), as suggested by observations tracking individual orbiting stars. Upcoming sub-millimetre very-long-baseline-interferometry (VLBI) images of Sgr A* carried out by the Event-Horizon-Telescope Collaboration (EHTC) are expected to provide critical evidence for the existence of this supermassive BH. We assess our present ability to use EHTC images to determine if they correspond to a Kerr BH as predicted by Einstein's theory of general relativity (GR) or to a BH in alternative theories of gravity. To this end, we perform general-relativistic magnetohydrodynamical (GRMHD) simulations and use general-relativistic radiative transfer (GRRT) calculations to generate synthetic shadow images of a magnetised accretion flow onto a Kerr BH. In addition, and for the first time, we perform GRMHD simulations and GRRT calculations for a dilaton BH, which we take as a representative solution of an alternative theory of gravity. Adopting the VLBI configuration from the 2017 EHTC campaign, we find that it could be extremely difficult to distinguish between BHs from different theories of gravity, thus highlighting that great caution is needed when interpreting BH images as tests of GR.Comment: Published in Nature Astronomy on 16.04.18 (including supplementary information); simulations at https://blackholecam.org/telling_bhs_apart