-
Fundamental Physics Opportunities with the Next-Generation Event Horizon Telescope
Authors:
Dimitry Ayzenberg,
Lindy Blackburn,
Richard Brito,
Silke Britzen,
Avery E. Broderick,
Raúl Carballo-Rubio,
Vitor Cardoso,
Andrew Chael,
Koushik Chatterjee,
Yifan Chen,
Pedro V. P. Cunha,
Hooman Davoudiasl,
Peter B. Denton,
Sheperd S. Doeleman,
Astrid Eichhorn,
Marshall Eubanks,
Yun Fang,
Arianna Foschi,
Christian M. Fromm,
Peter Galison,
Sushant G. Ghosh,
Roman Gold,
Leonid I. Gurvits,
Shahar Hadar,
Aaron Held
, et al. (23 additional authors not shown)
Abstract:
The Event Horizon Telescope (EHT) Collaboration recently published the first images of the supermassive black holes in the cores of the Messier 87 and Milky Way galaxies. These observations have provided a new means to study supermassive black holes and probe physical processes occurring in the strong-field regime. We review the prospects of future observations and theoretical studies of supermass…
▽ More
The Event Horizon Telescope (EHT) Collaboration recently published the first images of the supermassive black holes in the cores of the Messier 87 and Milky Way galaxies. These observations have provided a new means to study supermassive black holes and probe physical processes occurring in the strong-field regime. We review the prospects of future observations and theoretical studies of supermassive black hole systems with the next-generation Event Horizon Telescope (ngEHT), which will greatly enhance the capabilities of the existing EHT array. These enhancements will open up several previously inaccessible avenues of investigation, thereby providing important new insights into the properties of supermassive black holes and their environments. This review describes the current state of knowledge for five key science cases, summarising the unique challenges and opportunities for fundamental physics investigations that the ngEHT will enable.
△ Less
Submitted 4 December, 2023;
originally announced December 2023.
-
Testing the Existence of Event Horizons against Rotating Reflecting Surfaces
Authors:
Joost de Kleuver,
Thomas Bronzwaer,
Heino Falcke,
Ramesh Narayan,
Yosuke Mizuno,
Oliver Porth,
Hector Olivares
Abstract:
Recently the Event Horizon Telescope observed black holes at event horizon scales for the first time, enabling us to now test the existence of event horizons. Although event horizons have by definition no observable features, one can look for their non-existence. In that case, it is likely that there is some kind of surface, which like any other surface could absorb (and thermally emit) and/or ref…
▽ More
Recently the Event Horizon Telescope observed black holes at event horizon scales for the first time, enabling us to now test the existence of event horizons. Although event horizons have by definition no observable features, one can look for their non-existence. In that case, it is likely that there is some kind of surface, which like any other surface could absorb (and thermally emit) and/or reflect radiation. In this paper, we study the potential observable features of such rotating reflecting surfaces. We construct a general description of reflecting surfaces in arbitrary spacetimes. This is used to define specific models for static and rotating reflecting surfaces, of which we study the corresponding light paths and synthetic images. This is done by numerical integration of the geodesic equation and by the use of the general relativistic radiative transfer code RAPTOR. The reflecting surface creates an infinite set of ring-like features in synthetic images inside the photon ring. There is a central ring in the middle and higher order rings subsequently lie exterior to each other converging to the photon ring. The shape and size of the ring features change only slightly with the radius of the surface R, spin a and inclination i, resulting in all cases in features inside the 'shadow region'. We conclude that rotating reflecting surfaces have clear observable features and that the Event Horizon Telescope is able to observe the difference between reflecting surfaces and an event horizon for high reflectivities. Such reflecting surface models can be excluded, which strengthens the conclusion that the black hole shadow indeed indicates the existence of an event horizon.
△ Less
Submitted 9 November, 2023;
originally announced November 2023.
-
Magnetic flux eruptions at the root of time-lags in low-luminosity AGN
Authors:
Jesse Vos,
Jordy Davelaar,
Hector Olivares,
Christiaan Brinkerink,
Heino Falcke
Abstract:
Sagittarius A$^\ast$ is a compact radio source at the center of the Milky Way that has not conclusively shown evidence for the presence of a relativistic jet. Nevertheless, indirect methods at radio frequencies do indicate consistent outflow signatures. Temporal shifts between features in frequency bands are known as time lags, associated with flares or outflows of the accretion system. It is poss…
▽ More
Sagittarius A$^\ast$ is a compact radio source at the center of the Milky Way that has not conclusively shown evidence for the presence of a relativistic jet. Nevertheless, indirect methods at radio frequencies do indicate consistent outflow signatures. Temporal shifts between features in frequency bands are known as time lags, associated with flares or outflows of the accretion system. It is possible to gain information on the emission and outflow mechanics by interpreting these time lags. By means of a combined general-relativistic magnetrohydrodynamical and radiative transfer modeling, we study the origin of the time-lags for magnetically arrested disc models at three black hole spins ($a_\ast \in \{ -0.9375, 0, 0.9375 \}$). We exclusively modeled the emission from the source across a frequency range of $ν$ = 19-47 GHz. Our study also includes a targeted `slow light' study for one of the best-fitting `fast light' windows. We recovered the observational time-lag relations in various windows of our simulated light curves. The theoretical interpretation of these most promising time-lag windows is threefold; i) a magnetic flux eruption perturbs the jet-disc boundary and creates a flux tube, ii) the flux tube orbits and creates a clear emission feature, and iii) the flux tube interacts with the jet-disc boundary. The best-fitting windows have an intermediate (i=30$^\circ$/50$^\circ$) inclination and zero-BH-spin. The targeted `slow light' study did not yield better-fitting time-lag results, which indicates that the fast vs. slow light paradigm is often not intuitively understood and is likely influential in timing-sensitive studies. Sophisticated general-relativistic magnetrohydrodynamical models consistently capture the observational time-lag behavior, which is rooted in the complex dynamic interplay between the flux tube and coupled disk-jet system.
△ Less
Submitted 7 June, 2024; v1 submitted 25 October, 2023;
originally announced October 2023.
-
Plasmoid identification and statistics in two-dimensional Harris sheet and GRMHD simulations
Authors:
Jesse Vos,
Hector Olivares,
Benoit Cerutti,
Monika Moscibrodzka
Abstract:
Magnetic reconnection is a ubiquitous phenomenon for magnetized plasmas and leads to the rapid reconfiguration of magnetic field lines. During reconnection events, plasma is heated and accelerated until the magnetic field lines enclose and capture the plasma within a circular configuration. These plasmoids could therefore observationally manifest themselves as hot spots that are associated with fl…
▽ More
Magnetic reconnection is a ubiquitous phenomenon for magnetized plasmas and leads to the rapid reconfiguration of magnetic field lines. During reconnection events, plasma is heated and accelerated until the magnetic field lines enclose and capture the plasma within a circular configuration. These plasmoids could therefore observationally manifest themselves as hot spots that are associated with flaring behavior in supermassive black hole systems, such as Sagittarius A$^\ast$. We have developed a novel algorithm for identifying plasmoid structures, which incorporates watershed and custom closed contouring steps. From the identified plasmoids, we determine the plasma characteristics and energetics in magnetohydrodynamical simulations. The algorithm's performance is showcased for a high-resolution suite of axisymmetric ideal and resistive magnetohydrodynamical simulations of turbulent accretion discs surrounding a supermassive black hole. For validation purposes, we also evaluate several Harris current sheets that are well-investigated in the literature. Interestingly, we recover the characteristic power-law distribution of plasmoid sizes for both the black hole and Harris sheet simulations. This indicates that while the dynamics are vastly different, with different dominant plasma instabilities, the plasmoid creation behavior is similar. Plasmoid occurrence rates for resistive general relativistic magnetohydrodynamical simulations are significantly higher than for the ideal counterpart. Moreover, the largest identified plasmoids are consistent with sizes typically assumed for semi-analytical interpretation of observations. We recover a positive correlation between the plasmoid formation rate and a decrease in black-hole-horizon-penetrating magnetic flux. The developed algorithm has enabled an extensive quantitative analysis of plasmoid formation in black hole accretion simulations.
△ Less
Submitted 14 November, 2023; v1 submitted 6 September, 2023;
originally announced September 2023.
-
First-order hyperbolic formulation of the pure tetrad teleparallel gravity theory
Authors:
Ilya Peshkov,
Héctor Olivares,
Evgeniy Romenski
Abstract:
This paper presents a derivation of a first-order reduction and 3+1 decomposition of the teleparallel equivalent of general relativity (TEGR) in the pure-tetrad formulation (no spin connection). Our analysis demonstrates that in vacuum spacetimes, our 3+1 TEGR equations has the principal part of the differential operator equivalent to the one of tetrad reformulation of general relativity by Estabr…
▽ More
This paper presents a derivation of a first-order reduction and 3+1 decomposition of the teleparallel equivalent of general relativity (TEGR) in the pure-tetrad formulation (no spin connection). Our analysis demonstrates that in vacuum spacetimes, our 3+1 TEGR equations has the principal part of the differential operator equivalent to the one of tetrad reformulation of general relativity by Estabrook, Robinson, Wahlquist, and Buchman and Bardeen, and therefore the presented 3+1 decomposition of TEGR also admits a symmetric hyperbolic formulation, a desirable property for ensuring well-posedness of the initial value problem. Furthermore, the structure of the 3+1 equations possess a lot of similarities with the equations of relativistic electrodynamics and the recently proposed dGREM tetrad-reformulation of general relativity.
△ Less
Submitted 12 October, 2025; v1 submitted 24 November, 2022;
originally announced November 2022.
-
A new first-order formulation of the Einstein equations exploiting analogies with electrodynamics
Authors:
H. Olivares,
I. M. Peshkov,
E. R. Most,
F. M. Guercilena,
L. J. Papenfort
Abstract:
The Einstein and Maxwell equations are both systems of hyperbolic equations which need to satisfy a set of elliptic constraints throughout evolution. However, while electrodynamics (EM) and magnetohydrodynamics (MHD) have benefited from a large number of evolution schemes that are able to enforce these constraints and are easily applicable to curvilinear coordinates, unstructured meshes, or N-body…
▽ More
The Einstein and Maxwell equations are both systems of hyperbolic equations which need to satisfy a set of elliptic constraints throughout evolution. However, while electrodynamics (EM) and magnetohydrodynamics (MHD) have benefited from a large number of evolution schemes that are able to enforce these constraints and are easily applicable to curvilinear coordinates, unstructured meshes, or N-body simulations, many of these techniques cannot be straightforwardly applied to existing formulations of the Einstein equations. We develop a 3+1 a formulation of the Einstein equations which shows a striking resemblance to the equations of relativistic MHD and to EM in material media. The fundamental variables of this formulation are the frame fields, their exterior derivatives, and the Nester-Witten and Sparling forms. These mirror the roles of the electromagnetic 4-potential, the electromagnetic field strengths, the field excitations and the electric current. The role of the lapse function and shift vector, corresponds exactly to that of the scalar electric potential. The formulation, that we name dGREM (for differential forms, General Relativity and Electro-Magnetism), is manifestly first order and flux-conservative, which makes it suitable for high-resolution shock capturing schemes and finite-element methods. Being derived as a system of equations in exterior derivatives, it is directly applicable to any coordinate system and to unstructured meshes, and leads to a natural discretization potentially suitable for the use of machine-precision constraint propagation techniques such as the Yee algorithm and constrained transport. Due to these properties, we expect this new formulation to be beneficial in simulations of many astrophysical systems, such as binary compact objects and core-collapse supernovae as well as cosmological simulations of the early universe.
△ Less
Submitted 9 November, 2021;
originally announced November 2021.
-
Impact of non-thermal particles on the spectral and structural properties of M87
Authors:
Christian M. Fromm,
Alejandro Cruz-Osorio,
Yosuke Mizuno,
Antonios Nathanail,
Ziri Younsi,
Oliver Porth,
Hector Olivares,
Jordy Davelaar,
Heino Falcke,
Michael Kramer,
Luciano Rezzolla
Abstract:
The recent 230 GHz observations of the Event Horizon Telescope (EHT) are able to image the innermost structure of the M87 and show a ring-like structure which is in agreement with thermal synchrotron emission generated in a torus surrounding a supermassive black hole. However, at lower frequencies M87 is characterised by a large-scale and edge-brightened jet with clear signatures of non-thermal em…
▽ More
The recent 230 GHz observations of the Event Horizon Telescope (EHT) are able to image the innermost structure of the M87 and show a ring-like structure which is in agreement with thermal synchrotron emission generated in a torus surrounding a supermassive black hole. However, at lower frequencies M87 is characterised by a large-scale and edge-brightened jet with clear signatures of non-thermal emission. In order to bridge the gap between these scales and to provide a theoretical interpretation of these observations we perform general relativistic magnetohydrodynamic simulations of accretion on to black holes and jet launching.
M87 has been the target for multiple observations across the entire electromagnetic spectrum. Among these VLBI observations provide unique details on the collimation profile of the jet down to several gravitational radii. In this work we aim to model the observed broad-band spectrum of M87 from the radio to the NIR regime and at the same time fit the jet structure as observed with Global mm-VLBI at 86 GHz. We use general relativistic magnetohydrodynamics and simulate the accretion of the magnetised plasma onto Kerr-black holes in 3D. The radiative signatures of these simulations are computed taking different electron distribution functions into account and a detailed parameter survey is performed in order to match the observations.
The results of our simulations show that magnetically arrested disks around fast spinning black holes ($a_\star\geq0.5$) together with a mixture of thermal and non-thermal particle distributions are able to model simultaneously the broad-band spectrum and the innermost jet structure of M87
△ Less
Submitted 3 November, 2021;
originally announced November 2021.
-
Comparison of the ion-to-electron temperature ratio prescription: GRMHD simulations with electron thermodynamics
Authors:
Yosuke Mizuno,
Christian M. Fromm,
Ziri Younsi,
Oliver Porth,
Hector Olivares,
Luciano Rezzolla
Abstract:
The Event Horizon Telescope (EHT) collaboration, an Earth-size sub-millimetre radio interferometer, recently captured the first images of the central supermassive black hole in M87. These images were interpreted as gravitationally-lensed synchrotron emission from hot plasma orbiting around the black hole. In the accretion flows around low-luminosity active galactic nuclei such as M87, electrons an…
▽ More
The Event Horizon Telescope (EHT) collaboration, an Earth-size sub-millimetre radio interferometer, recently captured the first images of the central supermassive black hole in M87. These images were interpreted as gravitationally-lensed synchrotron emission from hot plasma orbiting around the black hole. In the accretion flows around low-luminosity active galactic nuclei such as M87, electrons and ions are not in thermal equilibrium. Therefore, the electron temperature, which is important for the thermal synchrotron radiation at EHT frequencies of 230 GHz, is not independently determined. In this work, we investigate the commonly used parameterised ion-to-electron temperature ratio prescription, the so-called R-$β$ model, considering images at 230 GHz by comparing with electron-heating prescriptions obtained from general-relativistic magnetohydrodynamical (GRMHD) simulations of magnetised accretion flows in a Magnetically Arrested Disc (MAD) regime with different recipes for the electron thermodynamics. When comparing images at 230 GHz, we find a very good match between images produced with the R-$β$ prescription and those produced with the turbulent- and magnetic reconnection- heating prescriptions. Indeed, this match is on average even better than that obtained when comparing the set of images built with the R-$β$ prescription with either a randomly chosen image or with a time-averaged one. From this comparative study of different physical aspects, which include the image, visibilities, broadband spectra, and light curves, we conclude that, within the context of images at 230 GHz relative to MAD accretion flows around supermassive black holes, the commonly-used and simple R-$β$ model is able to reproduce well the various and more complex electron-heating prescriptions considered here.
△ Less
Submitted 17 June, 2021;
originally announced June 2021.
-
Constraints on black-hole charges with the 2017 EHT observations of M87*
Authors:
Prashant Kocherlakota,
Luciano Rezzolla,
Heino Falcke,
Christian M. Fromm,
Michael Kramer,
Yosuke Mizuno,
Antonios Nathanail,
Hector Olivares,
Ziri Younsi,
Kazunori Akiyama,
Antxon Alberdi,
Walter Alef,
Juan Carlos Algaba,
Richard Anantua,
Keiichi Asada,
Rebecca Azulay,
Anne-Kathrin Baczko,
David Ball,
Mislav Balokovic,
John Barrett,
Bradford A. Benson,
Dan Bintley,
Lindy Blackburn,
Raymond Blundell,
Wilfred Boland
, et al. (212 additional authors not shown)
Abstract:
Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87*…
▽ More
Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes.
△ Less
Submitted 19 May, 2021;
originally announced May 2021.
-
Fuzzball Shadows: Emergent Horizons from Microstructure
Authors:
Fabio Bacchini,
Daniel R. Mayerson,
Bart Ripperda,
Jordy Davelaar,
Héctor Olivares,
Thomas Hertog,
Bert Vercnocke
Abstract:
We study the physical properties of four-dimensional, string-theoretical, horizonless "fuzzball" geometries by imaging their shadows. Their microstructure traps light rays straying near the would-be horizon on long-lived, highly redshifted chaotic orbits. In fuzzballs sufficiently near the scaling limit this creates a shadow much like that of a black hole, while avoiding the paradoxes associated w…
▽ More
We study the physical properties of four-dimensional, string-theoretical, horizonless "fuzzball" geometries by imaging their shadows. Their microstructure traps light rays straying near the would-be horizon on long-lived, highly redshifted chaotic orbits. In fuzzballs sufficiently near the scaling limit this creates a shadow much like that of a black hole, while avoiding the paradoxes associated with an event horizon. Observations of the shadow size and residual glow can potentially discriminate between fuzzballs away from the scaling limit and alternative models of black compact objects.
△ Less
Submitted 20 October, 2021; v1 submitted 22 March, 2021;
originally announced March 2021.
-
Using space-VLBI to probe gravity around Sgr A*
Authors:
C. M. Fromm,
Y. Mizuno,
Z. Younsi,
H. Olivares,
O. Porth,
M. De Laurentis,
H. Falcke,
M. Kramer,
L. Rezzolla
Abstract:
The Event Horizon Telescope (EHT) will soon provide the first high-resolution images of the Galactic Centre supermassive black hole (SMBH) candidate Sagittarius A* (Sgr A*), enabling us to probe gravity in the strong-field regime. Besides studying the accretion process in extreme environments, the obtained data and reconstructed images could be used to investigate the underlying spacetime structur…
▽ More
The Event Horizon Telescope (EHT) will soon provide the first high-resolution images of the Galactic Centre supermassive black hole (SMBH) candidate Sagittarius A* (Sgr A*), enabling us to probe gravity in the strong-field regime. Besides studying the accretion process in extreme environments, the obtained data and reconstructed images could be used to investigate the underlying spacetime structure. In its current configuration, the EHT is able to distinguish between a rotating Kerr black hole and a horizon-less object like a boson star. Future developments can increase the ability of the EHT to tell different spacetimes apart. We investigate the capability of an advanced EHT concept, including an orbiting space antenna, to image and distinguish different spacetimes around Sgr A*. We use GRMHD simulations of accreting compact objects (Kerr and dilaton black holes, as well as boson stars) and compute their radiative signatures via general relativistic radiative transfer calculations. To facilitate comparison with upcoming and future EHT observations we produce realistic synthetic data including the source variability, diffractive and refractive scattering while incorporating the observing array, including a space antenna. From the generated synthetic observations we dynamically reconstructed black hole shadow images using regularised Maximum Entropy methods. We employ a genetic algorithm to optimise the orbit of the space antenna with respect to improved imaging capabilities and u-v-plane coverage of the combined array (ground array and space antenna and developed a new method to probe the source variability in Fourier space. The inclusion of an orbiting space antenna improves the capability of the EHT to distinguish the spin of Kerr black holes and dilaton black holes based on reconstructed radio images and complex visibilities.
△ Less
Submitted 21 January, 2021;
originally announced January 2021.
-
Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole
Authors:
Dimitrios Psaltis,
Lia Medeiros,
Pierre Christian,
Feryal Ozel,
Kazunori Akiyama,
Antxon Alberdi,
Walter Alef,
Keiichi Asada,
Rebecca Azulay,
David Ball,
Mislav Balokovic,
John Barrett,
Dan Bintley,
Lindy Blackburn,
Wilfred Boland,
Geoffrey C. Bower,
Michael Bremer,
Christiaan D. Brinkerink,
Roger Brissenden,
Silke Britzen,
Dominique Broguiere,
Thomas Bronzwaer,
Do-Young Byun,
John E. Carlstrom,
Andrew Chael
, et al. (163 additional authors not shown)
Abstract:
The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p…
▽ More
The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources.
△ Less
Submitted 2 October, 2020;
originally announced October 2020.
-
Two-moment scheme for general-relativistic radiation hydrodynamics: a systematic description and new applications
Authors:
Lukas R. Weih,
Hector Olivares,
Luciano Rezzolla
Abstract:
We provide a systematic description of the steps necessary -- and of the potential pitfalls to be encountered -- when implementing a two-moment scheme within an Implicit-Explicit (IMEX) scheme to include radiative-transfer contributions in numerical simulations of general-relativistic (magneto-)hydrodynamics. We make use of the M1 closure, which provides an exact solution for the optically thin an…
▽ More
We provide a systematic description of the steps necessary -- and of the potential pitfalls to be encountered -- when implementing a two-moment scheme within an Implicit-Explicit (IMEX) scheme to include radiative-transfer contributions in numerical simulations of general-relativistic (magneto-)hydrodynamics. We make use of the M1 closure, which provides an exact solution for the optically thin and thick limit, and an interpolation between these limits. Special attention is paid to the efficient solution of the emerging set of implicit conservation equations. In particular, we present an efficient method for solving these equations via the inversion of a $4\times 4$-matrix within an IMEX scheme. While this method relies on a few approximations, it offers a very good compromise between accuracy and computational efficiency. After a large number of tests in special relativity, we couple our new radiation code, \texttt{FRAC}, with the general-relativistic magnetohydrodynamics code \texttt{BHAC} to investigate the radiative Michel solution, namely, the problem of spherical accretion onto a black hole in the presence of a radiative field. By performing the most extensive exploration of the parameter space for this problem, we find that the accretion's efficiency can be expressed in terms of physical quantities such as temperature, $T$, luminosity, $L$, and black-hole mass, $M$, via the expression $\varepsilon=(L/L_{\rm Edd})/(\dot{M}/\dot{M}_{\rm Edd})= 7.41\times 10^{-7}\left(T/10^6\,\mathrm{K}\right)^{0.22} \left(L/L_\odot\right)^{0.48} \left(M/M_\odot\right)^{0.48}$, where $L_{\mathrm{Edd}}$ and $\dot{M}_{\mathrm{Edd}}$ are the Eddington luminosity and accretion rate, respectively. Finally, we also consider the accretion problem away from spherical symmetry, finding that the solution is stable under perturbations in the radiation field.
△ Less
Submitted 19 May, 2020; v1 submitted 30 March, 2020;
originally announced March 2020.
-
Plasmoid formation in global GRMHD simulations and AGN flares
Authors:
Antonios Nathanail,
Christian M. Fromm,
Oliver Porth,
Hector Olivares,
Ziri Younsi,
Yosuke Mizuno,
Luciano Rezzolla
Abstract:
One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes,…
▽ More
One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes, reconnection layers can be developed and destroyed rapidly during the turbulent evolution of the flow. We present a series of two-dimensional general-relativistic magnetohydrodynamic simulations of tori accreting onto rotating black holes focusing our attention on the formation and evolution of current sheets. Initially, the tori are endowed with a poloidal magnetic field having a multi-loop structure along the radial direction and with an alternating polarity. During reconnection processes, plasmoids and plasmoid chains are developed leading to a flaring activity and hence to a variable electromagnetic luminosity. We describe the methods developed to track automatically the plasmoids that are generated and ejected during the simulation, contrasting the behaviour of multi-loop initial data with that encountered in typical simulations of accreting black holes having initial dipolar field composed of one loop only. Finally, we discuss the implications that our results have on the variability to be expected in accreting supermassive black holes.
△ Less
Submitted 5 February, 2020;
originally announced February 2020.
-
Deep Horizon; a machine learning network that recovers accreting black hole parameters
Authors:
Jeffrey van der Gucht,
Jordy Davelaar,
Luc Hendriks,
Oliver Porth,
Hector Olivares,
Yosuke Mizuno,
Christian M. Fromm,
Heino Falcke
Abstract:
The Event Horizon Telescope recently observed the first shadow of a black hole. Images like this can potentially be used to test or constrain theories of gravity and deepen the understanding in plasma physics at event horizon scales, which requires accurate parameter estimations. In this work, we present Deep Horizon, two convolutional deep neural networks that recover the physical parameters from…
▽ More
The Event Horizon Telescope recently observed the first shadow of a black hole. Images like this can potentially be used to test or constrain theories of gravity and deepen the understanding in plasma physics at event horizon scales, which requires accurate parameter estimations. In this work, we present Deep Horizon, two convolutional deep neural networks that recover the physical parameters from images of black hole shadows. We investigate the effects of a limited telescope resolution and observations at higher frequencies. We trained two convolutional deep neural networks on a large image library of simulated mock data. The first network is a Bayesian deep neural regression network and is used to recover the viewing angle $i$, and position angle, mass accretion rate $\dot{M}$, electron heating prescription $R_{\rm high}$ and the black hole mass $M_{\rm BH}$. The second network is a classification network that recovers the black hole spin $a$. We find that with the current resolution of the Event Horizon Telescope, it is only possible to accurately recover a limited number of parameters of a static image, namely the mass and mass accretion rate. Since potential future space-based observing missions will operate at frequencies above 230 GHz, we also investigated the applicability of our network at a frequency of 690 GHz. The expected resolution of space-based missions is higher than the current resolution of the Event Horizon Telescope, and we show that Deep Horizon can accurately recover the parameters of simulated observations with a comparable resolution to such missions.
△ Less
Submitted 27 May, 2020; v1 submitted 29 October, 2019;
originally announced October 2019.
-
General relativistic resistive magnetohydrodynamics with robust primitive variable recovery for accretion disk simulations
Authors:
Bart Ripperda,
Fabio Bacchini,
Oliver Porth,
Elias R. Most,
Hector Olivares,
Antonios Nathanail,
Luciano Rezzolla,
Jannis Teunissen,
Rony Keppens
Abstract:
Recent advances in black hole astrophysics, particularly the first visual evidence of a supermassive black hole at the center of the galaxy M87 by the Event Horizon Telescope (EHT), and the detection of an orbiting "hot spot" nearby the event horizon of Sgr A* in the Galactic center by the Gravity Collaboration, require the development of novel numerical methods to understand the underlying plasma…
▽ More
Recent advances in black hole astrophysics, particularly the first visual evidence of a supermassive black hole at the center of the galaxy M87 by the Event Horizon Telescope (EHT), and the detection of an orbiting "hot spot" nearby the event horizon of Sgr A* in the Galactic center by the Gravity Collaboration, require the development of novel numerical methods to understand the underlying plasma microphysics. Non-thermal emission related to such hot spots is conjectured to originate from plasmoids that form due to magnetic reconnection in thin current layers in the innermost accretion zone. Resistivity plays a crucial role in current sheet formation, magnetic reconnection, and plasmoid growth in black hole accretion disks and jets. We included resistivity in the three-dimensional general-relativistic magnetohydrodynamics (GRMHD) code BHAC and present the implementation of an Implicit-Explicit scheme to treat the stiff resistive source terms of the GRMHD equations. The algorithm is tested in combination with adaptive mesh refinement to resolve the resistive scales and a constrained transport method to keep the magnetic field solenoidal. Several novel methods for primitive variable recovery, a key part in relativistic magnetohydrodynamics codes, are presented and compared for accuracy, robustness, and efficiency. We propose a new inversion strategy that allows for resistive-GRMHD simulations of low gas-to-magnetic pressure ratio and highly magnetized regimes as applicable for black hole accretion disks, jets, and neutron star magnetospheres. We apply the new scheme to study the effect of resistivity on accreting black holes, accounting for dissipative effects as reconnection.
△ Less
Submitted 6 August, 2019; v1 submitted 16 July, 2019;
originally announced July 2019.
-
Constrained transport and adaptive mesh refinement in the Black Hole Accretion Code
Authors:
Hector Olivares,
Oliver Porth,
Jordy Davelaar,
Elias R. Most,
Christian M. Fromm,
Yosuke Mizuno,
Ziri Younsi,
Luciano Rezzolla
Abstract:
Worldwide very long baseline radio interferometry arrays are expected to obtain horizon-scale images of supermassive black hole candidates as well as of relativistic jets in several nearby active galactic nuclei. This motivates the development of models for magnetohydrodynamic flows in strong gravitational fields. The Black Hole Accretion Code (BHAC) intends to aid with the modelling of such sourc…
▽ More
Worldwide very long baseline radio interferometry arrays are expected to obtain horizon-scale images of supermassive black hole candidates as well as of relativistic jets in several nearby active galactic nuclei. This motivates the development of models for magnetohydrodynamic flows in strong gravitational fields. The Black Hole Accretion Code (BHAC) intends to aid with the modelling of such sources by means of general relativistic magnetohydrodynamical (GRMHD) simulations in arbitrary stationary spacetimes. New additions were required to guarantee an accurate evolution of the magnetic field when small and large scales are captured simultaneously. We discuss the adaptive mesh refinement (AMR) techniques employed in BHAC, essential to keep several problems computationally tractable, as well as staggered-mesh-based constrained transport (CT) algorithms to preserve the divergence-free constraint of the magnetic field, including a general class of prolongation operators for face-allocated variables compatible with them. Through several standard tests, we show that the choice of divergence-control method can produce qualitative differences in simulations of scientifically relevant accretion problems. We demonstrate the ability of AMR to reduce the computational costs of accretion simulations while sufficiently resolving turbulence from the magnetorotational instability. In particular, we describe a simulation of an accreting Kerr black hole in Cartesian coordinates using AMR to follow the propagation of a relativistic jet while self-consistently including the jet engine, a problem set up-for which the new AMR implementation is particularly advantageous. The CT methods and AMR strategies discussed here are being employed in the simulations performed with BHAC used in the generation of theoretical models for the Event Horizon Telescope Collaboration.
△ Less
Submitted 25 June, 2019;
originally announced June 2019.
-
Modeling non-thermal emission from the jet-launching region of M 87 with adaptive mesh refinement
Authors:
J. Davelaar,
H. Olivares,
O. Porth,
T. Bronzwaer,
M. Janssen,
F. Roelofs,
Y. Mizuno,
C. M. Fromm,
H. Falcke,
L. Rezzolla
Abstract:
The galaxy M 87 harbors a kiloparsec-scale relativistic jet, whose origin coincides with a supermassive black hole. Observational mm-VLBI campaigns are capable of resolving the jet-launching region at the scale of the event horizon. In order to provide a context for interpreting these observations, realistic general-relativistic magnetohydrodynamical (GRMHD) models of the accretion flow are constr…
▽ More
The galaxy M 87 harbors a kiloparsec-scale relativistic jet, whose origin coincides with a supermassive black hole. Observational mm-VLBI campaigns are capable of resolving the jet-launching region at the scale of the event horizon. In order to provide a context for interpreting these observations, realistic general-relativistic magnetohydrodynamical (GRMHD) models of the accretion flow are constructed. The characteristics of the observed spectral-energy distribution (SED) depend on the shape of the electrons' energy-distribution function (eDF). The dependency on the eDF is omitted in the modeling of the first Event Horizon Telescope results. In this work, we aim to model the M 87 SED from radio up to NIR/optical frequencies using a thermal-relativistic Maxwell- Jüttner distribution, as well as a relativistic $κ$-distribution function. The electrons are injected based on sub-grid, particle-in-cell parametrizations for sub-relativistic reconnection. A GRMHD simulation in Cartesian-Kerr-Schild coordinates, using eight levels of adaptive mesh refinement (AMR), forms the basis of our model. To obtain spectra and images, the GRMHD data is post-processed with the ray-tracing code RAPTOR, which is capable of ray tracing through AMR GRMHD simulation data. We obtain radio spectra in both the thermal-jet and $κ$-jet models consistent with radio observations. Additionally, the $κ$-jet models also recover the NIR/optical emission. The models recover the observed source sizes and core shifts and obtain a jet power of $\approx 10^{43}$ ergs/s. In the $κ$-jet models, both the accretion rates and jet powers are approximately two times lower than the thermal-jet model. The frequency cut-off observed at $ν\approx 10^{15}$ Hz is recovered when the accelerator size is $10^6$ - $10^8$ cm, this could potentially point to an upper limit for plasmoid sizes in the jet of M 87.
△ Less
Submitted 24 June, 2019;
originally announced June 2019.
-
The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project
Authors:
Oliver Porth,
Koushik Chatterjee,
Ramesh Narayan,
Charles F. Gammie,
Yosuke Mizuno,
Peter Anninos,
John G. Baker,
Matteo Bugli,
Chi-kwan Chan,
Jordy Davelaar,
Luca Del Zanna,
Zachariah B. Etienne,
P. Chris Fragile,
Bernard J. Kelly,
Matthew Liska,
Sera Markoff,
Jonathan C. McKinney,
Bhupendra Mishra,
Scott C. Noble,
Héctor Olivares,
Ben Prather,
Luciano Rezzolla,
Benjamin R. Ryan,
James M. Stone,
Niccolò Tomei
, et al. (3 additional authors not shown)
Abstract:
Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic…
▽ More
Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.
△ Less
Submitted 5 August, 2019; v1 submitted 9 April, 2019;
originally announced April 2019.
-
How to tell an accreting boson star from a black hole
Authors:
Hector Olivares,
Ziri Younsi,
Christian M. Fromm,
Mariafelicia De Laurentis,
Oliver Porth,
Yosuke Mizuno,
Heino Falcke,
Michael Kramer,
Luciano Rezzolla
Abstract:
The capability of the Event Horizon Telescope (EHT) to image the nearest supermassive black hole candidates at horizon-scale resolutions offers a novel means to study gravity in its strongest regimes and to test different models for these objects. Here, we study the observational appearance at 230 GHz of a surfaceless black hole mimicker, namely a non-rotating boson star, in a scenario consistent…
▽ More
The capability of the Event Horizon Telescope (EHT) to image the nearest supermassive black hole candidates at horizon-scale resolutions offers a novel means to study gravity in its strongest regimes and to test different models for these objects. Here, we study the observational appearance at 230 GHz of a surfaceless black hole mimicker, namely a non-rotating boson star, in a scenario consistent with the properties of the accretion flow onto Sgr A*. To this end, we perform general relativistic magnetohydrodynamic simulations followed by general relativistic radiative transfer calculations in the boson star space-time. Synthetic reconstructed images considering realistic astronomical observing conditions show that, despite qualitative similarities, the differences in the appearance of a black hole -- either rotating or not -- and a boson star of the type considered here are large enough to be detectable. These differences arise from dynamical effects directly related to the absence of an event horizon, in particular, the accumulation of matter in the form of a small torus or a spheroidal cloud in the interior of the boson star, and the absence of an evacuated high-magnetization funnel in the polar regions. The mechanism behind these effects is general enough to apply to other horizonless and surfaceless black hole mimickers, strengthening confidence in the ability of the EHT to identify such objects via radio observations.
△ Less
Submitted 6 August, 2020; v1 submitted 23 September, 2018;
originally announced September 2018.
-
The Current Ability to Test Theories of Gravity with Black Hole Shadows
Authors:
Yosuke Mizuno,
Ziri Younsi,
Christian M. Fromm,
Oliver Porth,
Mariafelicia De Laurentis,
Hector Olivares,
Heino Falcke,
Michael Kramer,
Luciano Rezzolla
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 April, 2018;
originally announced April 2018.
-
The Black Hole Accretion Code: adaptive mesh refinement and constrained transport
Authors:
Hector Olivares,
Oliver Porth,
Yosuke Mizuno
Abstract:
With the forthcoming VLBI images of Sgr A* and M87, simulations of accretion flows onto black holes acquire a special importance to aid with the interpretation of the observations and to test the predictions of different accretion scenarios, including those coming from alternative theories of gravity. The Black Hole Accretion Code (BHAC) is a new multidimensional general-relativistic magnetohydron…
▽ More
With the forthcoming VLBI images of Sgr A* and M87, simulations of accretion flows onto black holes acquire a special importance to aid with the interpretation of the observations and to test the predictions of different accretion scenarios, including those coming from alternative theories of gravity. The Black Hole Accretion Code (BHAC) is a new multidimensional general-relativistic magnetohydrondynamics (GRMHD) module for the MPI-AMRVAC framework. It exploits its adaptive mesh refinement techniques (AMR) to solve the equations of ideal magnetohydrodynamics in arbitrary curved spacetimes with a significant speedup and saving in computational cost. In a previous work, this was shown using a Generalized Lagrange Multiplier (GLM) to enforce the solenoidal constraint of the magnetic field. While GLM is fully compatible with MPI-AMRVAC's AMR infrastructure, we found that simulations were sensible to the divergence control technique employed, resulting in an improved behavior for those using Constrained Transport (CT). However, cell-centered CT is incompatible with AMR, and several modifications were required to make AMR compatible with staggered CT. We present here preliminary results of these new additions, which achieved machine precision fulfillment of the solenoidal constraint and a significant speedup in a problem close to the intended scientific application.
△ Less
Submitted 2 February, 2018;
originally announced February 2018.
-
The Black Hole Accretion Code
Authors:
Oliver Porth,
Hector Olivares,
Yosuke Mizuno,
Ziri Younsi,
Luciano Rezzolla,
Monika Moscibrodzka,
Heino Falcke,
Michael Kramer
Abstract:
We present the black hole accretion code (BHAC), a new multidimensional general-relativistic magnetohydrodynamics module for the MPI-AMRVAC framework. BHAC has been designed to solve the equations of ideal general-relativistic magnetohydrodynamics in arbitrary spacetimes and exploits adaptive mesh refinement techniques with an efficient block-based approach. Several spacetimes have already been im…
▽ More
We present the black hole accretion code (BHAC), a new multidimensional general-relativistic magnetohydrodynamics module for the MPI-AMRVAC framework. BHAC has been designed to solve the equations of ideal general-relativistic magnetohydrodynamics in arbitrary spacetimes and exploits adaptive mesh refinement techniques with an efficient block-based approach. Several spacetimes have already been implemented and tested. We demonstrate the validity of BHAC by means of various one-, two-, and three-dimensional test problems, as well as through a close comparison with the HARM3D code in the case of a torus accreting onto a black hole. The convergence of a turbulent accretion scenario is investigated with several diagnostics and we find accretion rates and horizon-penetrating fluxes to be convergent to within a few percent when the problem is run in three dimensions. Our analysis also involves the study of the corresponding thermal synchrotron emission, which is performed by means of a new general-relativistic radiative transfer code, BHOSS. The resulting synthetic intensity maps of accretion onto black holes are found to be convergent with increasing resolution and are anticipated to play a crucial role in the interpretation of horizon-scale images resulting from upcoming radio observations of the source at the Galactic Center.
△ Less
Submitted 18 June, 2017; v1 submitted 29 November, 2016;
originally announced November 2016.
-
BlackHoleCam: fundamental physics of the Galactic center
Authors:
C. Goddi,
H. Falcke,
M. Kramer,
L. Rezzolla,
C. Brinkerink,
T. Bronzwaer,
R. Deane,
M. De Laurentis,
G. Desvignes,
J. R. J. Davelaar,
F. Eisenhauer,
R. Eatough,
R. Fraga-Encinas,
C. M. Fromm,
S. Gillessen,
A. Grenzebach,
S. Issaoun,
M. Janßen,
R. Konoplya,
T. P. Krichbaum,
R. Laing,
K. Liu,
R. -S. Lu,
Y. Mizuno,
M. Moscibrodzka
, et al. (14 additional authors not shown)
Abstract:
Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme lim…
▽ More
Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.
△ Less
Submitted 7 February, 2017; v1 submitted 28 June, 2016;
originally announced June 2016.
-
Simulations of recoiling black holes: adaptive mesh refinement and radiative transfer
Authors:
Zakaria Meliani,
Yosuke Mizuno,
Hector Olivares,
Oliver Porth,
Luciano Rezzolla,
Ziri Younsi
Abstract:
(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 accur…
▽ More
(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.
△ Less
Submitted 10 November, 2016; v1 submitted 27 June, 2016;
originally announced June 2016.