General relativistic radiative transfer in black hole systems

Accretion onto compact objects plays a central role in high-energy astrophysics. The presence of a compact object considerably alters the structure and dynamics of the accreting plasma, as well as its radiative emissions. For accreting black holes in Active Galactic Nuclei (AGN) this is especially true. A significant fraction of the emission may originate or pass near the event horizon. Strong gravity modifies the radiation from an emission source. Photons no longer propagate in straight lines and experience frequency shifts. Gravitational lensing creates multiple images of an emission source, further modifying its temporal and spectral properties. Addressing these effects, the first part of this thesis formulates the equations of radiative transfer for particles with and without mass in a manifestly covariant form. Using ray-tracing, the observed images and line emission from accretion disks and tori are calculated. The effects of absorption, emission and optical depth gradients are investigated. The second part of this thesis examines scattering in general relativity. The general relativistic Compton scattering kernel and its angular moments are expressed in closed-form for the first time, in terms of hypergeometric functions. This has the advantage of being fast, accurate and not restricted by specific energy ranges. The results are in perfect agreement with semi-analytic calculations and Monte-Carlo simulations of Compton scattering of monochromatic emission lines. Finally, I investigate the effects of variability in the accretion flow. Two models are considered: a plasmoid on a Keplerian orbit around a black hole and a magnetically-driven plasmoid ejection from the disk corona. Deriving a new time-dependent radiative transfer formulation, I calculate this variable emission, presenting the results in the form of spectrograms and lightcurves

Flares in the Galactic Centre – I. Orbiting flux tubes in magnetically arrested black hole accretion discs

Recent observations of Sgr A* by the GRAVITY instrument have astrometrically tracked infrared (IR) flares at distances of ∼10 gravitational radii (rg). In this paper, we study a model for the flares based on 3D general relativistic magnetohydrodynamic (GRMHD) simulations of magnetically arrested accretion discs (MADs) that exhibit violent episodes of flux escape from the black hole magnetosphere. These events are attractive for flare modelling for several reasons: (i) the magnetically dominant regions can resist being disrupted via magnetorotational turbulence and shear; (ii) the orientation of the magnetic field is predominantly vertical as suggested by the GRAVITY data; and (iii) the magnetic reconnection associated with the flux eruptions could yield a self-consistent means of particle heating/acceleration during the flare events. In this analysis, we track erupted flux bundles and provide distributions of sizes, energies, and plasma parameter. In our simulations, the orbits tend to circularize at a range of radii from ∼5 to 40rg⁠. The magnetic energy contained within the flux bundles ranges up to ∼1040erg⁠, enough to power IR and X-ray flares. We find that the motion within the magnetically supported flow is substantially sub-Keplerian, in tension with the inferred period–radius relation of the three GRAVITY flares

Flares in the Galactic Centre - I:Orbiting flux tubes in magnetically arrested black hole accretion discs

Recent observations of SgrA* by the GRAVITY instrument have astrometrically tracked infrared flares (IR) at distances of $\sim 10$ gravitational radii ($r_g$). In this paper, we study a model for the flares based on 3D general relativistic magnetohydrodynamic (GRMHD) simulations of magnetically arrested accretion disks (MADs) which exhibit violent episodes of flux escape from the black hole magnetosphere. These events are attractive for flare modeling for several reasons: i) the magnetically dominant regions can resist being disrupted via magneto-rotational turbulence and shear, ii) the orientation of the magnetic field is predominantly vertical as suggested by the GRAVITY data, iii) magnetic reconnection associated with the flux eruptions could yield a self-consistent means of particle heating/acceleration during the flare events. In this analysis we track erupted flux bundles and provide distributions of sizes, energies and plasma parameter. In our simulations, the orbits tend to circularize at a range of radii from $\sim 5-40 r_g$. The magnetic energy contained within the flux bundles ranges up to $\sim10^{40}$ erg, enough to power IR and X-ray flares. We find that the motion within the magnetically supported flow is substantially sub-Keplerian, in tension with the inferred period-radius relation of the three GRAVITY flares.Comment: accepted for publication by MNRAS, 18-Jan-202

Jet-torus connection in radio galaxies: Relativistic hydrodynamics and synthetic emission

High-resolution Very-Long-Baseline Interferometry observations of active galactic nuclei have revealed asymmetric structures in the jets of radio galaxies. These asymmetric structures may be due to internal asymmetries in the jet, could be induced by the different conditions in the surrounding ambient medium including the obscuring torus, or a combination of the two. In this paper we investigate the influence of the ambient medium (including the obscuring torus) on the observed properties of jets from radio galaxies. We performed special-relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets using the special-relativistic hydrodynamics code \texttt{Ratpenat}, which is based on a second-order accurate finite-volume method and an approximate Riemann solver. Using a newly developed emission code to compute the electromagnetic emission, we have investigated the influence of different ambient medium and torus configurations on the jet structure and subsequently computed the non-thermal emission produced by the jet and the thermal absorption due to the torus. To better compare the emission simulations with observations we produced synthetic radio maps, taking into account the properties of the observatory. The detailed analysis of our simulations shows that the observed asymmetries can be produced by the interaction of the jet with the ambient medium and by the absorption properties of the obscuring torus.Comment: 14 pages, 17 figures, submitted to A&

Aberrational Effects for Shadows of Black Holes

In this paper, we discuss how the shadow of a Kerr black hole depends on the motion of the observer. In particular, we derive an analytical formula for the boundary curve of the shadow for an observer moving with given four-velocity at given Boyer--Lindquist coordinates. We visualize the shadow for various values of parameters.Comment: 12 pages, 3 figures; Proceedings of the 524. WE-Heraeus-Seminar held at the Physikzentrum, Bad Honnef, Germany, 17.--23.2.201

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 (χ ≡ 2GM/Rc2), 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 6), 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

Visibility of black hole shadows in low-luminosity AGN

Accreting black holes tend to display a characteristic dark central region called the black hole shadow, which depends only on space–time/observer geometry and which conveys information about the black hole’s mass and spin. Conversely, the observed central brightness depression, or image shadow, additionally depends on the morphology of the emission region. In this paper, we investigate the astrophysical requirements for observing a meaningful black hole shadow in GRMHD-based models of accreting black holes. In particular, we identify two processes by which the image shadow can differ from the black hole shadow: evacuation of the innermost region of the accretion flow, which can render the image shadow larger than the black hole shadow, and obscuration of the black hole shadow by optically thick regions of the accretion flow, which can render the image shadow smaller than the black hole shadow, or eliminate it altogether. We investigate in which models the image shadows of our models match their corresponding black hole shadows, and in which models the two deviate from each other. We find that, given a compact and optically thin emission region, our models allow for measurement of the black hole shadow size to an accuracy of 5 per cent. We show that these conditions are generally met for all MAD simulations we considered, as well as some of the SANE simulations

Black hole parameter estimation with synthetic very long baseline interferometry data from the ground and from space

Context. The Event Horizon Telescope (EHT) has imaged the shadow of the supermassive black hole in M 87. A library of general relativistic magnetohydrodynamics (GMRHD) models was fit to the observational data, providing constraints on black hole parameters. Aims. We investigate how much better future experiments can realistically constrain these parameters and test theories of gravity. Methods. We generated realistic synthetic 230 GHz data from representative input models taken from a GRMHD image library for M 87, using the 2017, 2021, and an expanded EHT array. The synthetic data were run through an automated data reduction pipeline used by the EHT. Additionally, we simulated observations at 230, 557, and 690 GHz with the Event Horizon Imager (EHI) Space VLBI concept. Using one of the EHT parameter estimation pipelines, we fit the GRMHD library images to the synthetic data and investigated how the black hole parameter estimations are affected by different arrays and repeated observations. Results. Repeated observations play an important role in constraining black hole and accretion parameters as the varying source structure is averaged out. A modest expansion of the EHT already leads to stronger parameter constraints in our simulations. High-frequency observations from space with the EHI rule out all but ∼15% of the GRMHD models in our library, strongly constraining the magnetic flux and black hole spin. The 1σ constraints on the black hole mass improve by a factor of five with repeated high-frequency space array observations as compared to observations with the current ground array. If the black hole spin, magnetization, and electron temperature distribution can be independently constrained, the shadow size for a given black hole mass can be tested to ∼0.5% with the EHI space array, which allows tests of deviations from general relativity. With such a measurement, high-precision tests of the Kerr metric become within reach from observations of the Galactic Center black hole Sagittarius A*

Black Hole Flares: Ejection of Accreted Magnetic Flux through 3D Plasmoid-mediated Reconnection

Magnetic reconnection can power bright, rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high-resolution (5376 × 2304 × 2304 cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current sheet forms in a transient, nonaxisymmetric, low-density magnetosphere within the inner few Schwarzschild radii. Magnetic flux bundles escape from the event horizon through reconnection at the universal plasmoid-mediated rate in this current sheet. The reconnection feeds on the highly magnetized plasma in the jets and heats the plasma that ends up trapped in flux bundles to temperatures proportional to the jet's magnetization. The escaped flux bundles can complete a full orbit as low-density hot spots, consistent with Sgr A* observations by the GRAVITY interferometer. Reconnection near the horizon produces sufficiently energetic plasma to explain flares from accreting black holes, such as the TeV emission observed from M87. The drop in the mass accretion rate during the flare and the resulting low-density magnetosphere make it easier for very-high-energy photons produced by reconnection-accelerated particles to escape. The extreme-resolution results in a converged plasmoid-mediated reconnection rate that directly determines the timescales and properties of the flare