Showing 1–50 of 363 results for author: Narayan, R

Accessing a New Population of Supermassive Black Holes with Extensions to the Event Horizon Telescope

Authors: Xinyue Alice Zhang, Angelo Ricarte, Dominic W. Pesce, Michael D. Johnson, Neil Nagar, Ramesh Narayan, Venkatessh Ramakrishnan, Sheperd Doeleman, Daniel C. M. Palumbo

Abstract: The Event Horizon Telescope has produced resolved images of the supermassive black holes Sgr A* and M87*, which present the largest shadows on the sky. In the next decade, technological improvements and extensions to the array will enable access to a greater number of sources, unlocking studies of a larger population of supermassive black holes through direct imaging. In this paper, we identify 12… ▽ More The Event Horizon Telescope has produced resolved images of the supermassive black holes Sgr A* and M87*, which present the largest shadows on the sky. In the next decade, technological improvements and extensions to the array will enable access to a greater number of sources, unlocking studies of a larger population of supermassive black holes through direct imaging. In this paper, we identify 12 of the most promising sources beyond Sgr A* and M87* based on their angular size and millimeter flux density. For each of these sources, we make theoretical predictions for their observable properties by ray tracing general relativistic magnetohydrodynamic models appropriately scaled to each target's mass, distance, and flux density. We predict that these sources would have somewhat higher Eddington ratios than M87*, which may result in larger optical and Faraday depths than previous EHT targets. Despite this, we find that visibility amplitude size constraints can plausibly recover masses within a factor of 2, although the unknown jet contribution remains a significant uncertainty. We find that the linearly polarized structure evolves substantially with Eddington ratio, with greater evolution at larger inclinations, complicating potential spin inferences for inclined sources. We discuss the importance of 345 GHz observations, milli-Jansky baseline sensitivity, and independent inclination constraints for future observations with upgrades to the Event Horizon Telescope (EHT) through ground updates with the next-generation EHT (ngEHT) program and extensions to space through the Black Hole Explorer (BHEX). △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 21 pages, 9 figures, 1 table, submitted to ApJ

The Black Hole Explorer: Motivation and Vision

Authors: Michael D. Johnson, Kazunori Akiyama, Rebecca Baturin, Bryan Bilyeu, Lindy Blackburn, Don Boroson, Alejandro Cardenas-Avendano, Andrew Chael, Chi-kwan Chan, Dominic Chang, Peter Cheimets, Cathy Chou, Sheperd S. Doeleman, Joseph Farah, Peter Galison, Ronald Gamble, Charles F. Gammie, Zachary Gelles, Jose L. Gomez, Samuel E. Gralla, Paul Grimes, Leonid I. Gurvits, Shahar Hadar, Kari Haworth, Kazuhiro Hada , et al. (43 additional authors not shown)

Abstract: We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery… ▽ More We present the Black Hole Explorer (BHEX), a mission that will produce the sharpest images in the history of astronomy by extending submillimeter Very-Long-Baseline Interferometry (VLBI) to space. BHEX will discover and measure the bright and narrow "photon ring" that is predicted to exist in images of black holes, produced from light that has orbited the black hole before escaping. This discovery will expose universal features of a black hole's spacetime that are distinct from the complex astrophysics of the emitting plasma, allowing the first direct measurements of a supermassive black hole's spin. In addition to studying the properties of the nearby supermassive black holes M87* and Sgr A*, BHEX will measure the properties of dozens of additional supermassive black holes, providing crucial insights into the processes that drive their creation and growth. BHEX will also connect these supermassive black holes to their relativistic jets, elucidating the power source for the brightest and most efficient engines in the universe. BHEX will address fundamental open questions in the physics and astrophysics of black holes that cannot be answered without submillimeter space VLBI. The mission is enabled by recent technological breakthroughs, including the development of ultra-high-speed downlink using laser communications, and it leverages billions of dollars of existing ground infrastructure. We present the motivation for BHEX, its science goals and associated requirements, and the pathway to launch within the next decade. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: Proceedings for SPIE Astronomical Telescopes and Instrumentation

Machine learning evaluation in the Global Event Processor FPGA for the ATLAS trigger upgrade

Authors: Zhixing Jiang, Scott Hauck, Dennis Yin, Bowen Zuo, Ben Carlson, Shih-Chieh Hsu, Allison Deiana, Rohin Narayan, Santosh Parajuli, Jeff Eastlack

Abstract: The Global Event Processor (GEP) FPGA is an area-constrained, performance-critical element of the Large Hadron Collider's (LHC) ATLAS experiment. It needs to very quickly determine which small fraction of detected events should be retained for further processing, and which other events will be discarded. This system involves a large number of individual processing tasks, brought together within th… ▽ More The Global Event Processor (GEP) FPGA is an area-constrained, performance-critical element of the Large Hadron Collider's (LHC) ATLAS experiment. It needs to very quickly determine which small fraction of detected events should be retained for further processing, and which other events will be discarded. This system involves a large number of individual processing tasks, brought together within the overall Algorithm Processing Platform (APP), to make filtering decisions at an overall latency of no more than 8ms. Currently, such filtering tasks are hand-coded implementations of standard deterministic signal processing tasks. In this paper we present methods to automatically create machine learning based algorithms for use within the APP framework, and demonstrate several successful such deployments. We leverage existing machine learning to FPGA flows such as hls4ml and fwX to significantly reduce the complexity of algorithm design. These have resulted in implementations of various machine learning algorithms with latencies of 1.2us and less than 5% resource utilization on an Xilinx XCVU9P FPGA. Finally, we implement these algorithms into the GEP system and present their actual performance. Our work shows the potential of using machine learning in the GEP for high-energy physics applications. This can significantly improve the performance of the trigger system and enable the ATLAS experiment to collect more data and make more discoveries. The architecture and approach presented in this paper can also be applied to other applications that require real-time processing of large volumes of data. △ Less

Submitted 7 May, 2024; originally announced June 2024.

Comments: 14 pages, 4 figures, 6 tables. Accepted by JINST on April 3, 2024

Multi-Zone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Authors: Hyerin Cho, Ben S. Prather, Kung-Yi Su, Ramesh Narayan, Priyamvada Natarajan

Abstract: Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. We describe and test a "multi-zone" technique which is designed to tackle this difficult problem in 3D general relativistic magnetohydrodynamic (GRMHD) simulations. We simulate accretion on a non-spinning SMBH ($a_*=0$) usi… ▽ More Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. We describe and test a "multi-zone" technique which is designed to tackle this difficult problem in 3D general relativistic magnetohydrodynamic (GRMHD) simulations. We simulate accretion on a non-spinning SMBH ($a_*=0$) using initial conditions from a large scale galaxy simulation, and achieve steady state over 8 decades in radius. The density scales with radius as $ρ\propto r^{-1}$ inside the Bondi radius $R_B$, which is located at $R_B=2\times 10^5 \,r_g$ ($\approx 60\,{\rm pc}$ for M87) where $r_g$ is the gravitational radius of the SMBH; the plasma-$β\sim$ unity, indicating an extended magnetically arrested state; the mass accretion rate $\dot{M}$ is $\approx 1\%$ of the analytical Bondi accretion rate $\dot{M}_B$; and there is continuous energy feedback out to $\approx 100R_B$ (or beyond $>\,{\rm kpc}$) at a rate $\approx 0.02 \dot{M}c^2$. Surprisingly, any ordered rotation in the external medium does not survive as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multi-zone method, we simulate GRMHD accretion for a wide range of Bondi radii, $R_{\rm B} \sim 10^2 - 10^7\,r_{\rm g}$, and find that $\dot{M}/\dot{M}_B\approx (R_B/6\, r_g)^{-0.5}$. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 19 pages, 9 figures, submitted to ApJ

Self-Gravitating Matter in Stationary and Axisymmetric Black Hole Spacetimes

Authors: Prashant Kocherlakota, Ramesh Narayan

Abstract: All black holes (BHs) in nature are expected to be described by the Kerr vacuum solution of general relativity. However, the Kerr solution comes with several difficulties such as the existence of Cauchy horizons, curvature singularities, and causality-violating regions. Attempts to resolve some of these issues include phenomenological BH models, which typically contain nontrivial matter content. W… ▽ More All black holes (BHs) in nature are expected to be described by the Kerr vacuum solution of general relativity. However, the Kerr solution comes with several difficulties such as the existence of Cauchy horizons, curvature singularities, and causality-violating regions. Attempts to resolve some of these issues include phenomenological BH models, which typically contain nontrivial matter content. We introduce a simple framework here to examine the properties of matter in such phenomenological models for a broad class of stationary and axisymmetric spinning BH spacetimes, generated from nonspinning seed solutions via a metric ansatz. We apply this framework to a representative set of spinning BH spacetimes and the non-spinning seeds from which they are derived. The models span different types of matter - fluids, scalar fields, electromagnetic fields. For each model, we calculate the timelike four-velocity of the matter and thereby identify the rest frame of the matter, both outside and inside the horizon. We then examine the spatial distribution of the matter rest-frame energy density $ε$ and the principal pressures. This provides a complete picture of how the matter moves, what its material properties are, and whether it obeys the classical energy conditions. Notably, at a horizon, the normal component of the pressure always satisfies $p_n = -ε$. We also investigate the expansions of the principal null congruences and explore the Hawking mass profiles of these spacetimes. These provide glimpses into the geometry of the stationary BH exterior as well as the nonstationary interior cosmology. The axisymmetric metric ansatz we work with can be used to generate new spinning solutions from a variety of nonspinning seeds. The matter in these models often satisfies the weak energy condition, at least in the BH exterior, and some models exhibit non-rigid, differential rotation. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 25 pages, 9 figures, 2 tables

Density jump for oblique collisionless shocks in pair plasmas: physical solutions

Authors: Antoine Bret, Colby C. Haggerty, Ramesh Narayan

Abstract: Collisionless shocks are frequently analyzed using the magnetohydrodynamics (MHD) formalism, even though MHD assumes a small mean free path. Yet, isotropy of pressure, fruit of binary collisions and assumed in MHD, may not apply in collisionless shocks. This is especially true within a magnetized plasma, where the field can stabilize an anisotropy. In a previous article \citep{BretJPP2022b}, a m… ▽ More Collisionless shocks are frequently analyzed using the magnetohydrodynamics (MHD) formalism, even though MHD assumes a small mean free path. Yet, isotropy of pressure, fruit of binary collisions and assumed in MHD, may not apply in collisionless shocks. This is especially true within a magnetized plasma, where the field can stabilize an anisotropy. In a previous article \citep{BretJPP2022b}, a model was presented capable of dealing with the anisotropies that may arise at the front crossing. It was solved for any orientation of the field with respect to the shock front. Yet, for some values of the upstream parameters, several downstream solutions were found. Here, we complete the work started in \cite{BretJPP2022b} by showing how to pick the physical solution out of the ones offered by the algebra. This is achieved by 2 means: 1) selecting the solution that has the downstream field obliquity closest to the upstream one. This criterion is exemplified on the parallel case and backed up by Particle-in-Cell simulations. 2) Filtering out solutions which do not satisfy a criteria already invoked to trim multiple solutions in MHD: the evolutionarity criterion, that we assume valid in the collisionless case. The end result is a model in which a given upstream configuration results in a unique, or none (like in MHD), downstream configuration. The largest departure from MHD is found for the case of a parallel shock. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 22 pages, 9 figures, to appear in Journal of Plasma Physics

Observational Signatures of AGN Feedback in the Morphology and the Ionization States of Milky Way-like Galaxies

Authors: Nadia Qutob, Razieh Emami, Kung-Yi Su, Randall Smith, Lars Hernquist, Dian P. Triani, Cameron Hummels, Drummond Fielding, Philip F. Hopkins, Rachel S. Somerville, David R. Ballantyne, Mark Vogelsberger, Grant Tremblay, James F. Steiner, Douglas Finkbeiner, Ramesh Narayan, Minjung Park, Josh Grindlay, Priyamvada Natarajan, Christopher C. Hayward, Dušan Kereš, Sam B. Ponnada, Sirio Belli, Rebecca Davies, Gabriel Maheson , et al. (2 additional authors not shown)

Abstract: We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. Three jet models, including cosmic ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. We examine the distribution of Mg II, O VI, and O VIII ions, along… ▽ More We make an in-depth analysis of different AGN jet models' signatures, inducing quiescence in galaxies with a halo mass of $10^{12} M_\odot$. Three jet models, including cosmic ray-dominant, hot thermal, and precessing kinetic jets, are studied at two energy flux levels each, compared to a jet-free, stellar feedback-only simulation. We examine the distribution of Mg II, O VI, and O VIII ions, alongside gas temperature and density profiles. Low-energy ions, like Mg II, concentrate in the ISM, while higher energy ions, e.g., O VIII, prevail at the AGN jet cocoon's edge. High-energy flux jets display an isotropic ion distribution with lower overall density. High-energy thermal or cosmic ray jets pressurize at smaller radii, significantly suppressing core density. The cosmic ray jet provides extra pressure support, extending cool and warm gas distribution. A break in the ion-to-mass ratio slope in O VI and O VIII is demonstrated in the ISM-to-CGM transition (between 10-30 kpc), growing smoothly towards the CGM at greater distances. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Comments: 17 pages, 13 figures

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.

Comments: 14 pages, 10 figures; submitted to Astronomy and Astrophysics

On the Universality of Energy Extraction from Black Hole Spacetimes

Authors: Koushik Chatterjee, Ziri Younsi, Prashant Kocherlakota, Ramesh Narayan

Abstract: The launching of astrophysical jets provides the most compelling observational evidence for direct extraction of black hole (BH) spin energy via the Blandford-Znajek (BZ) mechanism. Whilst it is know that spinning Kerr BHs within general relativity (GR) follow the BZ jet power relation, the nature of BH energy extraction in general theories of gravity has not been adequately addressed. This study… ▽ More The launching of astrophysical jets provides the most compelling observational evidence for direct extraction of black hole (BH) spin energy via the Blandford-Znajek (BZ) mechanism. Whilst it is know that spinning Kerr BHs within general relativity (GR) follow the BZ jet power relation, the nature of BH energy extraction in general theories of gravity has not been adequately addressed. This study performs the first comprehensive investigation of the BZ jet power relation by utilising a generalised BH spacetime geometry which describes parametric deviations from the Kerr metric of GR, yet recovers the Kerr metric in the limit that all deviation parameters vanish. Through performing and analysing an extensive suite of three-dimensional covariant magnetohydrodynamics (MHD) simulations of magnetised gas accretion onto these generalised BH spacetimes we find that the BZ jet power relation still holds, in some instances yielding jet powers far in excess of what can be produced by even extremal Kerr BHs. It is shown that the variation of the quadrupole moment of the BH can enhance or suppress the effects of BH spin, and by extension of frame-dragging. This variation greatly enhances or suppresses the observed jet power and underlying photon ring image asymmetry, introducing a previously unexplored yet important degeneracy in BH parameter inference. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: To be submitted to journal. Comments are welcome

Energy Extraction from Spinning Stringy Black Holes

Authors: Koushik Chatterjee, Prashant Kocherlakota, Ziri Younsi, Ramesh Narayan

Abstract: We perform the first numerical simulations modeling the inflow and outflow of magnetized plasma in the Kerr-Sen spacetime, which describes classical spinning black holes (BHs) in string theory. We find that the Blandford-Znajek (BZ) mechanism, which is believed to power astrophysical relativistic outflows or ``jets'', is valid even for BHs in an alternate theory of gravity, including near the extr… ▽ More We perform the first numerical simulations modeling the inflow and outflow of magnetized plasma in the Kerr-Sen spacetime, which describes classical spinning black holes (BHs) in string theory. We find that the Blandford-Znajek (BZ) mechanism, which is believed to power astrophysical relativistic outflows or ``jets'', is valid even for BHs in an alternate theory of gravity, including near the extremal limit. The BZ mechanism releases outward Poynting-flux-dominated plasma as frame-dragging forces magnetic field lines to twist. However, for nonspinning BHs, where the frame-dragging is absent, we find an alternate powering mechanism through the release of gravitational potential energy during accretion. Outflows from non-spinning stringy BHs can be approximately $250\%$ more powerful as compared to Schwarzschild BHs, due to their relatively smaller event horizon sizes and, thus, higher curvatures. Finally, by constructing the first synthetic images of near-extremal non-Kerr BHs from time-dependent simulations, we find that these can be ruled out by horizon-scale interferometric images of accreting supermassive BHs. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: To be submitted to journal. Comments are welcome

Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD

Authors: Hyerin Cho, Ben S. Prather, Ramesh Narayan, Priyamvada Natarajan, Kung-Yi Su, Angelo Ricarte, Koushik Chatterjee

Abstract: Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magneto-hydrodynamic (GRMHD) code KHARMA that allows us to span $7$ orders of magnitude in spatial scale, to simulate ac… ▽ More Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magneto-hydrodynamic (GRMHD) code KHARMA that allows us to span $7$ orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with Bondi radius $R_B\approx 2\times 10^5 \,GM_\bullet/c^2$, where $M_\bullet$ is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as $\sim r^{-1}$ rather than $r^{-3/2}$ and the accretion rate $\dot{M}$ is consequently suppressed by over 2 orders of magnitude below the Bondi rate $\dot{M}_B$. We find continuous energy feedback from the accretion flow to the external medium at a level of $\sim10^{-2}\dot{M}c^2 \sim 5\times 10^{-5} \dot{M}_B c^2$. Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback. △ Less

Submitted 1 December, 2023; v1 submitted 29 October, 2023; originally announced October 2023.

Comments: 8 pages with 5 pages of Appendices, 4 figures, accepted to ApJ Letters

Galaxy formation with Wave/Fuzzy Dark Matter: The core-halo structure

Authors: Alvaro Pozo, Razieh Emami, Philip Mocz, Tom Broadhurst, Lars Hernquist, Mark Vogelsberger, Randall Smith, Grant Tremblay, Ramesh Narayan, James Steiner, Josh Grindlay, George Smoot

Abstract: Dark matter-dominated cores have long been claimed for the well-studied local group dwarf galaxies. More recently, extended stellar halos have been uncovered around several of these dwarfs through deeper imaging and spectroscopy. Such core-halo structures are not a feature of conventional cold dark matter (CDM), based on collisionless particles where smooth, scale-free profiles are predicted. In c… ▽ More Dark matter-dominated cores have long been claimed for the well-studied local group dwarf galaxies. More recently, extended stellar halos have been uncovered around several of these dwarfs through deeper imaging and spectroscopy. Such core-halo structures are not a feature of conventional cold dark matter (CDM), based on collisionless particles where smooth, scale-free profiles are predicted. In contrast, smooth and prominent dark matter cores are predicted for Warm and Fuzzy/Wave Dark Matter (WDM/$ψ$DM) respectively. The question arises to what extent the visible stellar profiles should reflect this dark matter core structure. Here we compare cosmological hydrodynamical simulations of CDM, WDM $\&$ $ψ$DM, aiming to predict the stellar profiles for these three DM scenarios. We show that cores surrounded by extended halos are distinguishable for WDM and $ψ$DM, with the most prominent cores in the case of $ψ$DM, where the stellar density is enhanced in the core due to the presence of the relatively dense soliton. Our analysis demonstrates that such behavior does not appear in CDM, implying that the small-scale cut-off in the power spectrum present for WDM and $ψ$DM provides a core-halo transition. Consequently, we estimate the mass of the $ψ$DM particle at this core-halo transition point. Furthermore, we observe the anticipated asymmetry for $ψ$DM due to the soliton's random walk, a distinctive characteristic not found in the symmetric distributions of stars in Warm and CDM models. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 17 pages and 15 figures

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

IIDS: Design of Intelligent Intrusion Detection System for Internet-of-Things Applications

Authors: KG Raghavendra Narayan, Srijanee Mookherji, Vanga Odelu, Rajendra Prasath, Anish Chand Turlapaty, Ashok Kumar Das

Abstract: With rapid technological growth, security attacks are drastically increasing. In many crucial Internet-of-Things (IoT) applications such as healthcare and defense, the early detection of security attacks plays a significant role in protecting huge resources. An intrusion detection system is used to address this problem. The signature-based approaches fail to detect zero-day attacks. So anomaly-bas… ▽ More With rapid technological growth, security attacks are drastically increasing. In many crucial Internet-of-Things (IoT) applications such as healthcare and defense, the early detection of security attacks plays a significant role in protecting huge resources. An intrusion detection system is used to address this problem. The signature-based approaches fail to detect zero-day attacks. So anomaly-based detection particularly AI tools, are becoming popular. In addition, the imbalanced dataset leads to biased results. In Machine Learning (ML) models, F1 score is an important metric to measure the accuracy of class-level correct predictions. The model may fail to detect the target samples if the F1 is considerably low. It will lead to unrecoverable consequences in sensitive applications such as healthcare and defense. So, any improvement in the F1 score has significant impact on the resource protection. In this paper, we present a framework for ML-based intrusion detection system for an imbalanced dataset. In this study, the most recent dataset, namely CICIoT2023 is considered. The random forest (RF) algorithm is used in the proposed framework. The proposed approach improves 3.72%, 3.75% and 4.69% in precision, recall and F1 score, respectively, with the existing method. Additionally, for unsaturated classes (i.e., classes with F1 score < 0.99), F1 score improved significantly by 7.9%. As a result, the proposed approach is more suitable for IoT security applications for efficient detection of intrusion and is useful in further studies. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Toward General-Relativistic Magnetohydrodynamics Simulations in Stationary Non-Vacuum Spacetimes

Authors: Prashant Kocherlakota, Ramesh Narayan, Koushik Chatterjee, Alejandro Cruz-Osorio, Yosuke Mizuno

Abstract: Accretion of magnetized gas on compact astrophysical objects such as black holes has been successfully modeled using general relativistic magnetohydrodynamic (GRMHD) simulations. These simulations have largely been performed in the Kerr metric, which describes the spacetime of a vacuum and stationary spinning black hole (BH) in general relativity (GR). The simulations have revealed important clues… ▽ More Accretion of magnetized gas on compact astrophysical objects such as black holes has been successfully modeled using general relativistic magnetohydrodynamic (GRMHD) simulations. These simulations have largely been performed in the Kerr metric, which describes the spacetime of a vacuum and stationary spinning black hole (BH) in general relativity (GR). The simulations have revealed important clues on the physics of accretion and jets near the BH event horizon, and have been used to interpret recent Event Horizon Telescope images of the supermassive BHs, M87$^*$ and Sgr A$^*$. GRMHD simulations require the spacetime metric in horizon-penetrating coordinates such that all metric coefficients are regular at the event horizon. The Kerr metric and its electrically charged spinning analog, the Kerr-Newman metric, are currently the only metrics available in such coordinates. We report here horizon-penetrating forms of a large class of stationary, axisymmetric, spinning metrics. These can be used to carry out GRMHD simulations of accretion on spinning, nonvacuum BHs and non-BHs within GR, as well as accretion on spinning objects described by non-GR metric theories of gravity. △ Less

Submitted 5 March, 2024; v1 submitted 27 July, 2023; originally announced July 2023.

Comments: Corrected typo in eqs. 23, C20

Journal ref: ApJL 956, L11 (2023)

Recipes for Jet Feedback and Spin Evolution of Black Holes with Strongly-Magnetized Super-Eddington Accretion Disks

Authors: Angelo Ricarte, Ramesh Narayan, Brandon Curd

Abstract: A spinning black hole accreting from a disk of strongly magnetized plasma via a magnetically arrested disk is known to produce an efficient electromagnetic jet powered by the black hole's spin energy. We present general relativistic radiative magnetohydrodynamic simulations of magnetically arrested systems covering a range of sub- to super-Eddington accretion rates. Using the numerical results fro… ▽ More A spinning black hole accreting from a disk of strongly magnetized plasma via a magnetically arrested disk is known to produce an efficient electromagnetic jet powered by the black hole's spin energy. We present general relativistic radiative magnetohydrodynamic simulations of magnetically arrested systems covering a range of sub- to super-Eddington accretion rates. Using the numerical results from these simulations, we develop formulae to describe the magnetization, jet efficiency, and spin evolution of an accreting black hole as a function of its spin and accretion rate. A black hole with near-Eddington accretion experiences a mild degree of spin-down because of angular momentum loss through the jet, leading to an equilibrium spin of 0.8 rather than 1.0 at the Eddington limit. As the accretion rate increases above Eddington, the spin-down effect becomes progressively stronger, ultimately converging on previous predictions based on non-radiative simulations. In particular, spin evolution drives highly super-Eddington systems toward a black hole spin near zero. The formulae developed in this letter may be applied to galaxy and cosmological scale simulations that include black holes. If magnetically arrested disk accretion is common among supermassive black holes, the present results have broad implications for active galactic nucleus feedback and cosmological spin evolution. △ Less

Submitted 9 August, 2023; v1 submitted 10 July, 2023; originally announced July 2023.

Comments: 21 pages, 7 figures, 1 table. Accepted by ApJL. See the journal for downloadable data files as well as a helpful Python script referenced in the data availability statement

The EB-correlation in Resolved Polarized Images: Connections to Astrophysics of Black Holes

Authors: Razieh Emami, Sheperd S. Doeleman, Maciek Wielgus, Dominic Chang, Koushik Chatterjee, Randall Smith, Matthew Liska, James F. Steiner, Angelo Ricarte, Ramesh Narayan, Grant Tremblay, Douglas Finkbeiner, Lars Hernquist, Chi-Kwan Chan, Lindy Blackburn, Ben S. Prather, Paul Tiede, Avery E. Broderick, Mark Vogelsberger, Charles Alcock, Freek Roelofs

Abstract: We present an in-depth analysis of a newly proposed correlation function in visibility space, between the E and B modes of the linear polarization, hereafter the EB-correlation, for a set of time-averaged GRMHD simulations compared with the phase map from different semi-analytic models as well as the Event Horizon Telescope (EHT) 2017 data for M87* source. We demonstrate that the phase map of the… ▽ More We present an in-depth analysis of a newly proposed correlation function in visibility space, between the E and B modes of the linear polarization, hereafter the EB-correlation, for a set of time-averaged GRMHD simulations compared with the phase map from different semi-analytic models as well as the Event Horizon Telescope (EHT) 2017 data for M87* source. We demonstrate that the phase map of the time-averaged EB-correlation contains novel information that might be linked to the BH spin, accretion state and the electron temperature. A detailed comparison with a semi-analytic approach with different azimuthal expansion modes shows that to recover the morphology of the real/imaginary part of the correlation function and its phase, we require higher orders of these azimuthal modes. To extract the phase features, we propose to use the Zernike polynomial reconstruction developing an empirical metric to break degeneracies between models with different BH spins that are qualitatively similar. We use a set of different geometrical ring models with various magnetic and velocity field morphologies and show that both the image space and visibility based EB-correlation morphologies in MAD simulations can be explained with simple fluid and magnetic field geometries as used in ring models. SANEs by contrast are harder to model, demonstrating that the simple fluid and magnetic field geometries of ring models are not sufficient to describe them owing to higher Faraday Rotation depths. A qualitative comparison with the EHT data demonstrates that some of the features in the phase of EB-correlation might be well explained by the current models for BH spins as well as electron temperatures, while others may require a larger theoretical surveys. △ Less

Submitted 3 May, 2023; v1 submitted 30 April, 2023; originally announced May 2023.

Comments: 26 pages, 21 Figures

Linear analysis of the Kelvin-Helmholtz instability in relativistic magnetized symmetric flows

Authors: Anthony Chow, Michael E. Rowan, Lorenzo Sironi, Jordy Davelaar, Gianluigi Bodo, Ramesh Narayan

Abstract: We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed $c_{\rm s}$, Alfvén speed $v_{\rm A}$, and magnetic field orientation. For the instability asso… ▽ More We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed $c_{\rm s}$, Alfvén speed $v_{\rm A}$, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity onto the fluid-frame wavevector is larger than the projection of the Alfvén speed onto the same direction, i.e., shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed $v$, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{re} \equiv v/v_{\rm f\perp} \sqrt{(1-v_{\rm f\perp}^2)/(1-v^2)} \cosθ=\sqrt{2}$, where $v_{rm f\perp}=[v_A^2+c_{\rm s}^2(1-v_A^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and $θ$ is the laboratory-frame angle between the flow velocity and the wavevector projection onto the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes -- both for single objects and for merging binaries -- where the Alfvén speed may approach the speed of light. △ Less

Submitted 15 June, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: 11 pages, 7 figures, 1 table, Accepted for publication in Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 524, Issue 1, September 2023, Pages 90

Key Science Goals for the Next-Generation Event Horizon Telescope

Authors: Michael D. Johnson, Kazunori Akiyama, Lindy Blackburn, Katherine L. Bouman, Avery E. Broderick, Vitor Cardoso, R. P. Fender, Christian M. Fromm, Peter Galison, José L. Gómez, Daryl Haggard, Matthew L. Lister, Andrei P. Lobanov, Sera Markoff, Ramesh Narayan, Priyamvada Natarajan, Tiffany Nichols, Dominic W. Pesce, Ziri Younsi, Andrew Chael, Koushik Chatterjee, Ryan Chaves, Juliusz Doboszewski, Richard Dodson, Sheperd S. Doeleman , et al. (20 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) has led to the first images of a supermassive black hole, revealing the central compact objects in the elliptical galaxy M87 and the Milky Way. Proposed upgrades to this array through the next-generation EHT (ngEHT) program would sharply improve the angular resolution, dynamic range, and temporal coverage of the existing EHT observations. These improvements will u… ▽ More The Event Horizon Telescope (EHT) has led to the first images of a supermassive black hole, revealing the central compact objects in the elliptical galaxy M87 and the Milky Way. Proposed upgrades to this array through the next-generation EHT (ngEHT) program would sharply improve the angular resolution, dynamic range, and temporal coverage of the existing EHT observations. These improvements will uniquely enable a wealth of transformative new discoveries related to black hole science, extending from event-horizon-scale studies of strong gravity to studies of explosive transients to the cosmological growth and influence of supermassive black holes. Here, we present the key science goals for the ngEHT and their associated instrument requirements, both of which have been formulated through a multi-year international effort involving hundreds of scientists worldwide. △ Less

Submitted 21 April, 2023; originally announced April 2023.

Comments: 32 pages, 11 figures, accepted for publication in a special issue of Galaxies on the ngEHT (https://www.mdpi.com/journal/galaxies/special_issues/ngEHT_blackholes)

Journal ref: Galaxies 2023, 11(3), 61

BenchMD: A Benchmark for Unified Learning on Medical Images and Sensors

Authors: Kathryn Wantlin, Chenwei Wu, Shih-Cheng Huang, Oishi Banerjee, Farah Dadabhoy, Veeral Vipin Mehta, Ryan Wonhee Han, Fang Cao, Raja R. Narayan, Errol Colak, Adewole Adamson, Laura Heacock, Geoffrey H. Tison, Alex Tamkin, Pranav Rajpurkar

Abstract: Medical data poses a daunting challenge for AI algorithms: it exists in many different modalities, experiences frequent distribution shifts, and suffers from a scarcity of examples and labels. Recent advances, including transformers and self-supervised learning, promise a more universal approach that can be applied flexibly across these diverse conditions. To measure and drive progress in this dir… ▽ More Medical data poses a daunting challenge for AI algorithms: it exists in many different modalities, experiences frequent distribution shifts, and suffers from a scarcity of examples and labels. Recent advances, including transformers and self-supervised learning, promise a more universal approach that can be applied flexibly across these diverse conditions. To measure and drive progress in this direction, we present BenchMD: a benchmark that tests how well unified, modality-agnostic methods, including architectures and training techniques (e.g. self-supervised learning, ImageNet pretraining),perform on a diverse array of clinically-relevant medical tasks. BenchMD combines 19 publicly available datasets for 7 medical modalities, including 1D sensor data, 2D images, and 3D volumetric scans. Our benchmark reflects real-world data constraints by evaluating methods across a range of dataset sizes, including challenging few-shot settings that incentivize the use of pretraining. Finally, we evaluate performance on out-of-distribution data collected at different hospitals than the training data, representing naturally-occurring distribution shifts that frequently degrade the performance of medical AI models. Our baseline results demonstrate that no unified learning technique achieves strong performance across all modalities, leaving ample room for improvement on the benchmark. Code is released at https://github.com/rajpurkarlab/BenchMD. △ Less

Submitted 26 June, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Black Holes Up Close

Authors: Ramesh Narayan, Eliot Quataert

Abstract: Recent developments have ushered in a new era in the field of black hole astrophysics, providing our first direct view of the remarkable environment near black hole event horizons. These observations have enabled astronomers to confirm long-standing ideas on the physics of gas flowing into black holes with temperatures that are hundreds of times greater than at the center of the Sun. At the same t… ▽ More Recent developments have ushered in a new era in the field of black hole astrophysics, providing our first direct view of the remarkable environment near black hole event horizons. These observations have enabled astronomers to confirm long-standing ideas on the physics of gas flowing into black holes with temperatures that are hundreds of times greater than at the center of the Sun. At the same time, the observations have conclusively shown that light rays near a black hole experience large deflections which cause a dark shadow in the center of the image, an effect predicted by Einstein's theory of General Relativity. With further investment, this field is poised to deliver decades of advances in our understanding of gravity and black holes through new and stringent tests of General Relativity, as well as new insights into the role of black holes as the central engines powering a wide range of astronomical phenomena. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 15 pages, 4 Figs; Nature Perspective

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Using Machine Learning to Link Black Hole Accretion Flows with Spatially Resolved Polarimetric Observables

Authors: Richard Qiu, Angelo Ricarte, Ramesh Narayan, George N. Wong, Andrew Chael, Daniel Palumbo

Abstract: We introduce a new library of 535,194 model images of the supermassive black holes and Event Horizon Telescope (EHT) targets Sgr A* and M87*, computed by performing general relativistic radiative transfer calculations on general relativistic magnetohydrodynamics simulations. Then, to infer underlying black hole and accretion flow parameters (spin, inclination, ion-to-electron temperature ratio, an… ▽ More We introduce a new library of 535,194 model images of the supermassive black holes and Event Horizon Telescope (EHT) targets Sgr A* and M87*, computed by performing general relativistic radiative transfer calculations on general relativistic magnetohydrodynamics simulations. Then, to infer underlying black hole and accretion flow parameters (spin, inclination, ion-to-electron temperature ratio, and magnetic field polarity), we train a random forest machine learning model on various hand-picked polarimetric observables computed from each image. Our random forest is capable of making meaningful predictions of spin, inclination, and the ion-to-electron temperature ratio, but has more difficulty inferring magnetic field polarity. To disentangle how physical parameters are encoded in different observables, we apply two different metrics to rank the importance of each observable at inferring each physical parameter. Details of the spatially resolved linear polarization morphology stand out as important discriminators between models. Bearing in mind the theoretical limitations and incompleteness of our image library, for the real M87* data, our machinery favours high-spin retrograde models with large ion-to-electron temperature ratios. Due to the time-variable nature of these targets, repeated polarimetric imaging will further improve model inference as the EHT and next-generation (EHT) continue to develop and monitor their targets. △ Less

Submitted 9 February, 2023; v1 submitted 9 December, 2022; originally announced December 2022.

Comments: 24 pages, 27 figures

Density jump for oblique collisionless shocks in pair plasmas: allowed solutions

Authors: Antoine Bret, Ramesh Narayan

Abstract: Shockwaves in plasma are usually dealt with using Magnetohydrodynamics (MHD). Yet, MHD entails the assumption of a short mean free path, which is not fulfilled in a collisionless plasma. Recently, for pair plasmas, we devised a model allowing to account for kinetic effects within an MHD-like formalism. Its relies on an estimate of the anisotropy generated when crossing the front, with a subsequent… ▽ More Shockwaves in plasma are usually dealt with using Magnetohydrodynamics (MHD). Yet, MHD entails the assumption of a short mean free path, which is not fulfilled in a collisionless plasma. Recently, for pair plasmas, we devised a model allowing to account for kinetic effects within an MHD-like formalism. Its relies on an estimate of the anisotropy generated when crossing the front, with a subsequent assessment of the stability of this anisotropy in the downstream. We solved our model for parallel, perpendicular and switch-on shocks. Here we bridge between all these cases by treating the problem of an arbitrarily, but coplanar, oriented magnetic field. Even though the formalism presented is valid for anisotropic upstream temperatures, only the case of a cold upstream is solved. We find extra solutions which are not part of the MHD catalog, and a density jump that is notably less in the quasi parallel, highly magnetized, regime. Given the complexity of the calculations, this work is mainly devoted to the presentation of the mathematical aspect of our model. A forthcoming article will be devoted to the physics of the shocks here defined. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 16 pages, 8 figures, to appear in Journal of Plasma Physics

Probing plasma composition with the next generation Event Horizon Telescope (ngEHT)

Authors: Razieh Emami, Richard Anantua, Angelo Ricarte, Sheperd S. Doeleman, Avery Broderick, George Wong, Lindy Blackburn, Maciek Wielgus, Ramesh Narayan, Grant Tremblay, Charles Alcock, Lars Hernquist, Randall Smith, Matthew Liska, Priyamvada Natarajan, Mark Vogelsberger, Brandon Curd, Joana A. Kramer

Abstract: We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models including a constant electron to magnetic pressure (constant $β_e$ model) with a population… ▽ More We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models including a constant electron to magnetic pressure (constant $β_e$ model) with a population of non-thermal electrons as well as a R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations get enhanced per adding pairs. On the contrary, we show that at higher frequencies a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 14 pages, 9 figures

Tracing the hot spot motion using the next generation Event Horizon Telescope (ngEHT)

Authors: Razieh Emami, Paul Tiede, Sheperd S. Doeleman, Freek Roelofs, Maciek Wielgus, Lindy Blackburn, Matthew Liska, Koushik Chatterjee, Bart Ripperda, Antonio Fuentes, Avery Broderick, Lars Hernquist, Charles Alcock, Ramesh Narayan, Randall Smith, Grant Tremblay, Angelo Ricarte, He Sun, Richard Anantua, Yuri Y. Kovalev, Priyamvada Natarajan, Mark Vogelsberger

Abstract: We propose to trace the dynamical motion of a shearing hot spot near the SgrA* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its bo… ▽ More We propose to trace the dynamical motion of a shearing hot spot near the SgrA* source through a dynamical image reconstruction algorithm, StarWarps. Such a hot spot may form as the exhaust of magnetic reconnection in a current sheet near the black hole horizon. A hot spot that is ejected from the current sheet into an orbit in the accretion disk may shear and diffuse due to instabilities at its boundary during its orbit, resulting in a distinct signature. We subdivide the motion to two distinct phases; the first phase refers to the appearance of the hot spot modelled as a bright blob, followed by a subsequent shearing phase simulated as a stretched ellipse. We employ different observational arrays, including EHT(2017,2022) and the next generation event horizon telescope (ngEHTp1, ngEHT) arrays, in which few new additional sites are added to the observational array. We make dynamical image reconstructions for each of these arrays. Subsequently, we infer the hot spot phase in the first phase followed by the axes ratio and the ellipse area in the second phase. We focus on the direct observability of the orbiting hot spot in the sub-mm wavelength. Our analysis demonstrates that newly added dishes may easily trace the first phase as well as part of the second phase, before the flux is reduced substantially. The algorithm used in this work can be extended to any other types of the dynamical motion. Consequently, we conclude that the ngEHT is a key to directly observe the dynamical motions near variable sources, such as SgrA*. △ Less

Submitted 15 November, 2022; v1 submitted 12 November, 2022; originally announced November 2022.

Comments: 16 pages, 8 pages

Observational Signatures of Frame Dragging in Strong Gravity

Authors: Angelo Ricarte, Daniel C. M. Palumbo, Ramesh Narayan, Freek Roelofs, Razieh Emami

Abstract: Objects orbiting in the presence of a rotating massive body experience a gravitomagnetic frame-dragging effect, known as the Lense-Thirring effect, that has been experimentally confirmed in the weak-field limit. In the strong-field limit, near the horizon of a rotating black hole, frame dragging becomes so extreme that all objects must co-rotate with the black hole's angular momentum. In this work… ▽ More Objects orbiting in the presence of a rotating massive body experience a gravitomagnetic frame-dragging effect, known as the Lense-Thirring effect, that has been experimentally confirmed in the weak-field limit. In the strong-field limit, near the horizon of a rotating black hole, frame dragging becomes so extreme that all objects must co-rotate with the black hole's angular momentum. In this work, we perform general relativistic numerical simulations to identify observable signatures of frame dragging in the strong-field limit that appear when infalling gas is forced to flip its direction of rotation as it is being accreted. In total intensity images, infalling streams exhibit "S"-shaped features due to the switch in the tangential velocity. In linear polarization, a flip in the handedness of spatially resolved polarization ticks as a function of radius encodes a transition in the magnetic field geometry that occurs due to magnetic flux freezing in the dragged plasma. Using a network of telescopes around the world, the Event Horizon Telescope collaboration has demonstrated that it is now possible to directly image black holes on event horizon scales. We show that the phenomena described in this work would be accessible to the next-generation Event Horizon Telescope (ngEHT) and extensions of the array into space, which would produce spatially resolved images on event horizon scales with higher spatial resolution and dynamic range. △ Less

Submitted 2 November, 2022; originally announced November 2022.

Comments: Submitted to ApJL. 15 pages, 8 figures

Flux eruption events drive angular momentum transport in magnetically arrested accretion flows

Authors: Koushik Chatterjee, Ramesh Narayan

Abstract: We evolve two high-resolution general relativistic magnetohydrodynamic (GRMHD) simulations of advection-dominated accretion flows around non-spinning black holes (BHs), each over a duration $\sim 3\times 10^5\,GM_{\rm BH}/c^3$. One model captures the evolution of a weakly magnetized (SANE) disk and the other a magnetically arrested disk (MAD). Magnetic flux eruptions in the MAD model push out gas… ▽ More We evolve two high-resolution general relativistic magnetohydrodynamic (GRMHD) simulations of advection-dominated accretion flows around non-spinning black holes (BHs), each over a duration $\sim 3\times 10^5\,GM_{\rm BH}/c^3$. One model captures the evolution of a weakly magnetized (SANE) disk and the other a magnetically arrested disk (MAD). Magnetic flux eruptions in the MAD model push out gas from the disk and launch strong winds with outflow efficiencies at times reaching $10\%$ of the incoming accretion power. Despite the substantial power in these winds, average mass outflow rates remain small out to a radius $\sim100\,GM_{\rm BH}/c^2$, only reaching $\sim 60-80\%$ of the horizon accretion rate. The average outward angular momentum transport is primarily radial in both modes of accretion, but with a clear distinction: magnetic flux eruption-driven disk winds cause a strong vertical flow of angular momentum in the MAD model, while for the SANE model, the magnetorotational instability (MRI) moves angular momentum mostly equatorially through the disk. Further, we find that the MAD state is highly transitory and non-axisymmetric, with the accretion mode often changing to a SANE-like state following an eruption before reattaining magnetic flux saturation with time. The Reynolds stress changes direction during such transitions, with the MAD (SANE) state showing an inward (outward) stress, possibly pointing to intermittent MRI-driven accretion in MADs. Pinning down the nature of flux eruptions using next-generation telescopes will be crucial in understanding the flow of mass, magnetic flux and angular momentum in sub-Eddington accreting BHs like M87$^*$ and Sagittarius A$^*$. △ Less

Submitted 14 October, 2022; originally announced October 2022.

Comments: Submitted to ApJ. Comments are welcome

Observational Constraints on Direct Electron Heating in the Hot Accretion Flows in Sgr A* and M87*

Authors: Fu-Guo Xie, Ramesh Narayan, Feng Yuan

Abstract: An important parameter in the theory of hot accretion flows around black holes is $δ$, which describes the fraction of ``viscously'' dissipated energy in the accretion flow that goes directly into heating electrons. For a given mass accretion rate, the radiative efficiency of a hot accretion flow is determined by $δ$. Unfortunately, the value of $δ$ is hard to determine from first principles. The… ▽ More An important parameter in the theory of hot accretion flows around black holes is $δ$, which describes the fraction of ``viscously'' dissipated energy in the accretion flow that goes directly into heating electrons. For a given mass accretion rate, the radiative efficiency of a hot accretion flow is determined by $δ$. Unfortunately, the value of $δ$ is hard to determine from first principles. The recent Event Horizon Telescope Collaboration (EHTC) results on M87* and Sgr A* provide us with a different way of constraining $δ$. By combining the mass accretion rates in M87* and Sgr A* estimated by the EHTC with the measured bolometric luminosities of the two sources, we derive good constraints on the radiative efficiencies of the respective accretion flows. In parallel, we use a theoretical model of hot magnetically arrested disks (MAD) to calculate the expected radiative efficiency as a function of $δ$ (and accretion rate). By comparing the EHTC-derived radiative efficiencies with the theoretical results from MAD models, we find that Sgr A* requires $δ\ga 0.3$. %with the most likely value being $δ\sim 0.5$. A similar comparison in the case of M87* gives inconclusive results as there is still a large uncertainty in the accretion rate in this source. △ Less

Submitted 29 November, 2022; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 9 pages, 2 figures. To appear in ApJ

Unraveling Twisty Linear Polarization Morphologies in Black Hole Images

Authors: Razieh Emami, Angelo Ricarte, George N. Wong, Daniel Palumbo, Dominic Chang, Sheperd S. Doeleman, Avery Broaderick, Ramesh Narayan, Maciek Wielgus, Lindy Blackburn, Ben S. Prather, Andrew A. Chael, Richard Anantua, Koushik Chatterjee, Ivan Marti-Vidal, Jose L. Gomez, Kazunori Akiyama, Matthew Liska, Lars Hernquist, Grant Tremblay, Mark Vogelsberger, Charles Alcock, Randall Smith, James Steiner, Paul Tiede , et al. (1 additional authors not shown)

Abstract: We investigate general relativistic magnetohydrodynamic simulations (GRMHD) to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterising the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnet… ▽ More We investigate general relativistic magnetohydrodynamic simulations (GRMHD) to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterising the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnetic field structure in the emitting region. Moreover, the dependence of this twisty pattern on spin can be attributed to changes in the magnetic field geometry that occur due to the frame dragging. By studying an analytic ring model, we find that the roles of Doppler boosting and lensing are subdominant. Faraday rotation may cause a systematic shift in the linear polarization pattern, but we find that its impact is subdominant for models with strong magnetic fields and modest ion-to-electron temperature ratios. Models with weaker magnetic fields are much more strongly affected by Faraday rotation and have more complicated emission geometries than can be captured by a ring model. However, these models are currently disfavoured by the recent EHT observations of M87*. Our results suggest that linear polarization maps can provide a probe of the underlying magnetic field structure around a black hole, which may then be usable to indirectly infer black hole spins. The generality of these results should be tested with alternative codes, initial conditions, and plasma physics prescriptions. △ Less

Submitted 28 March, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 26 pages, 20 figures, accepted for publication in ApJ

GRRMHD Simulations of MAD Accretion Disks Declining from Super-Eddington to Sub-Eddington Accretion Rates

Authors: Brandon Curd, Ramesh Narayan

Abstract: We present two general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of magnetically arrested disks (MADs) around non-spinning ($a_*=0$) and spinning ($a_*=0.9$) supermassive black holes (BHs). In each simulation, the mass accretion rate is decreased with time such that we sample Eddington-scaled rates over the range $3 \gtrsim \dot{M}/\dot{M}_{\rm{Edd}}\gtrsim 0.3$. For the non… ▽ More We present two general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of magnetically arrested disks (MADs) around non-spinning ($a_*=0$) and spinning ($a_*=0.9$) supermassive black holes (BHs). In each simulation, the mass accretion rate is decreased with time such that we sample Eddington-scaled rates over the range $3 \gtrsim \dot{M}/\dot{M}_{\rm{Edd}}\gtrsim 0.3$. For the non-spinning BH model, the total and radiative efficiencies increase as the accretion rate decreases, varying over the range $η_{\rm{tot}}\sim9-16\%$ and $η_{\rm{rad}}\sim6-12\%$, respectively. This model shows very little jet activity. In contrast, the spinning BH model has a strong relativistic jet powered by spin energy extracted from the BH. The jet power declines with accretion rate such that $η_{\rm{jet}}\sim 18-39\%$ while the total and radiative efficiencies are $η_{\rm{tot}}\sim 64-100\%$ and $η_{\rm{rad}}\sim 45-79\%$, respectively. We confirm that mildly sub-Eddington disks can extract substantial power from a spinning BH, provided they are in the MAD state. The jet profile out to $100\, GM/c^2$ is roughly parabolic with a power-law index of $k\approx0.43-0.53$ during the sub-Eddington evolution. Both models show significant variability in the outgoing radiation which is likely associated with episodes of magnetic flux eruptions. The $a_*=0.9$ model shows semi-regular variations with a period of $\sim2000\, GM/c^3$ over the final $\sim10,000\, GM/c^3$ of the simulation, which suggests that magnetic flux eruptions may be an important source of quasi-periodic variability. For the simulated accretion rates, the $a_*=0$ model is spinning up while the $a_*=0.9$ model is spinning down. Spinup-spindown equilibrium of the BH will likely be achieved at $0.5 < a_{*,{\rm{eq}}} < 0.6$, assuming continuous accretion in the MAD state. △ Less

Submitted 24 September, 2022; originally announced September 2022.

Comments: 23 pages, 18 figures

Density jump as a function of magnetic field for switch-on collisionless shocks in pair plasmas

Authors: Antoine Bret, Ramesh Narayan

Abstract: The properties of collisionless shocks, like the density jump, are usually derived from magnetohydrodynamics (MHD), where isotropic pressures are assumed. Yet, in a collisionless plasma, an external magnetic field can sustain a stable anisotropy. In \cite{BretJPP2018}, we devised a model for the kinetic history of the plasma through the shock front, allowing to self-consistently compute the downst… ▽ More The properties of collisionless shocks, like the density jump, are usually derived from magnetohydrodynamics (MHD), where isotropic pressures are assumed. Yet, in a collisionless plasma, an external magnetic field can sustain a stable anisotropy. In \cite{BretJPP2018}, we devised a model for the kinetic history of the plasma through the shock front, allowing to self-consistently compute the downstream anisotropy, hence the density jump, in terms of the upstream parameters. This model dealt with the case of a parallel shock, where the magnetic field is normal to the front both in the upstream and the downstream. Yet, MHD also allows for shock solutions, the so-called switch-on solutions, where the field is normal to the front only in the upstream. This article consists in applying our model to these switch-on shocks. While MHD offers only 1 switch-on solution within a limited range of Alfvén Mach numbers, our model offers 2 kinds of solutions within a slightly different range of Alfvén Mach numbers. These 2 solutions are most likely the outcome of the intermediate and fast MHD shocks under our model. While the intermediate and fast shocks merge in MHD for the parallel case, they do not within our model. For simplicity, the formalism is restricted to non-relativistic shocks in pair plasmas where the upstream is cold. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 18 pages. 6 figures. To appear in Journal of Plasma Physics

Jets from SANE Super-Eddington Accretion Disks: Morphology, Spectra, and Their Potential as Targets for ngEHT

Authors: Brandon Curd, Razieh Emami, Richard Anantua, Daniel Palumbo, Sheperd Doeleman, Ramesh Narayan

Abstract: We present general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of super-Eddington accretion flows around supermassive black holes (SMBHs) which may apply to tidal disruption events (TDEs). We perform long duration ($t\geq81,200\, GM/c^3$) simulations which achieve mass accretion rates $\gtrsim 11$ times the Eddington rate and produce thermal synchrotron spectra and images of t… ▽ More We present general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of super-Eddington accretion flows around supermassive black holes (SMBHs) which may apply to tidal disruption events (TDEs). We perform long duration ($t\geq81,200\, GM/c^3$) simulations which achieve mass accretion rates $\gtrsim 11$ times the Eddington rate and produce thermal synchrotron spectra and images of their jets. The jet reaches a maximum velocity of $v/c \approx 0.5-0.9$, but the density weighted outflow velocity is $\sim0.2-0.35c$. Gas flowing beyond the funnel wall expands conically and drives a strong shock at the jet head while variable mass ejection along the jet axis results in internal shocks and dissipation. For a $T_i/T_e=1$ model, the radio/submillimeter spectra peak at $>100$ GHz and the luminosity increases with BH spin, exceeding $\sim 10^{41} \, \rm{erg\, s^{-1}}$ in the brightest models. The emission is extremely sensitive to $T_i/T_e$ as some models show an order of magnitude decrease in the peak frequency and up to four orders of magnitude decline in their radio/submillimeter luminosity as $T_i/T_e$ approaches 20. Assuming a maximum VLBI baseline distance of $10 \ {\rm{G}}λ$, 230 GHz images of $T_i/T_e=1$ models shows that the jet head may be bright enough for its motion to be captured with the EHT (ngEHT) at $D\lesssim110$ (180) Mpc at the $5σ$ significance level. Resolving emission from internal shocks requires $D\lesssim45$ Mpc for both the EHT or ngEHT. The 5 GHz emission in each model is dimmer ($\lesssim10^{36} \ {\rm{erg\, s^{-1}}}$) than upper limits placed on TDEs with no radio emission which suggests jets similar to our models may have gone undetected in previous observations. Our models suggest that the ngEHT may be utilized for $>230$ GHz radio/submillimeter followup of future TDEs. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Detectability of Wandering Intermediate-Mass Black Holes in the Milky Way Galaxy from Radio to X-rays

Authors: Bryan Seepaul, Fabio Pacucci, Ramesh Narayan

Abstract: Intermediate-mass black holes (IMBHs, $10^{3-6} \, \rm{M_\odot}$), are typically found at the center of dwarf galaxies and might be wandering, thus far undetected, in the Milky Way (MW). We use model spectra for advection-dominated accretion flows to compute the typical fluxes, in a range of frequencies spanning from radio to X-rays, emitted by a putative population of $10^5 \, \rm{M_\odot}$ IMBHs… ▽ More Intermediate-mass black holes (IMBHs, $10^{3-6} \, \rm{M_\odot}$), are typically found at the center of dwarf galaxies and might be wandering, thus far undetected, in the Milky Way (MW). We use model spectra for advection-dominated accretion flows to compute the typical fluxes, in a range of frequencies spanning from radio to X-rays, emitted by a putative population of $10^5 \, \rm{M_\odot}$ IMBHs wandering in five realistic, volume-weighted, MW environments. We predict that $\sim 27\%$ of the wandering IMBHs can be detected in the X-ray with Chandra, $\sim 37\%$ in the near-infrared with the Roman Space Telescope, $\sim 49\%$ in the sub-mm with CMB-S4 and $\sim 57\%$ in the radio with ngVLA. We find that the brightest fluxes are emitted by IMBHs passing through molecular clouds or cold neutral medium, where they are always detectable. We propose criteria to facilitate the selection of candidates in multi-wavelength surveys. Specifically, we compute the X-ray to optical ratio ($α_{\rm ox}$) and the optical to sub-mm ratio, as a function of the accretion rate of the IMBH. We show that at low rates the sub-mm emission of IMBHs is significantly higher than the optical, UV and X-ray emission. Finally, we place upper limits on the number $N_\bullet$ of these objects in the MW: $N_\bullet<2000$ and $N_\bullet<100$, based on our detectability expectations and current lack of detections in molecular clouds and cold neutral medium, respectively. These predictions will guide future searches of IMBHs in the MW, which will be instrumental to understanding their demographics and evolution. △ Less

Submitted 19 July, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: Accepted for publication in MNRAS. This is the final version of the manuscript. 11 pages, 7 figures. Includes predictions for Chandra, eROSITA, Planck, SCUBA, AXIS, Roman, CMB-S4, ngVLA

Analytical Model of Disk Evaporation and State Transitions in Accreting Black Holes

Authors: Hyerin Cho, Ramesh Narayan

Abstract: State transitions in black hole X-ray binaries are likely caused by gas evaporation from a thin accretion disk into a hot corona. We present a height-integrated version of this process which is suitable for analytical and numerical studies. With radius $r$ scaled to Schwarzschild units and coronal mass accretion rate $\dot{m}_c$ to Eddington units, the results of the model are independent of black… ▽ More State transitions in black hole X-ray binaries are likely caused by gas evaporation from a thin accretion disk into a hot corona. We present a height-integrated version of this process which is suitable for analytical and numerical studies. With radius $r$ scaled to Schwarzschild units and coronal mass accretion rate $\dot{m}_c$ to Eddington units, the results of the model are independent of black hole mass. State transitions should thus be similar in X-ray binaries and AGN. The corona solution consists of two power-law segments separated at a break radius $r_b \sim10^3 \,(α/0.3)^{-2}$, where $α$ is the viscosity parameter. Gas evaporates from the disk to the corona for $r>r_b$, and condenses back for $r<r_b$. At $r_b$, $\dot{m}_c$ reaches its maximum, $\dot{m}_{c,{\rm max}} \approx 0.02\, (α/0.3)^3$. If at $r\gg r_b$ the thin disk accretes with $\dot{m}_d < \dot{m}_{c,{\rm max}} $, then the disk evaporates fully before reaching $r_b$, giving the hard state. Otherwise, the disk survives at all radii, giving the thermal state. While the basic model considers only bremsstrahlung cooling and viscous heating, we also discuss a more realistic model which includes Compton cooling and direct coronal heating by energy transport from the disk. Solutions are again independent of black hole mass, and $r_b$ remains unchanged. This model predicts strong coronal winds for $r>r_b$, and a $T\sim 5\times 10^8\,{\rm K}$ Compton-cooled corona for $r < r_b$. Two-temperature effects are ignored, but may be important at small radii. △ Less

Submitted 23 June, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Comments: 40 pages, 12 figures, published in ApJ

Observational properties of puffy disks: radiative GRMHD spectra of mildly sub-Eddington accretion

Authors: Maciek Wielgus, Debora Lancova, Odele Straub, Wlodek Kluzniak, Ramesh Narayan, David Abarca, Agata Rozanska, Frederic Vincent, Gabriel Torok, Marek Abramowicz

Abstract: Numerical general relativistic radiative magnetohydrodynamic simulations of accretion disks around a stellar mass black hole with a luminosity above 0.5 of the Eddington value reveal their stratified, elevated vertical structure. We refer to these thermally stable numerical solutions as puffy disks. Above a dense and geometrically thin core of dimensionless thickness $h/r \sim 0.1$, crudely resemb… ▽ More Numerical general relativistic radiative magnetohydrodynamic simulations of accretion disks around a stellar mass black hole with a luminosity above 0.5 of the Eddington value reveal their stratified, elevated vertical structure. We refer to these thermally stable numerical solutions as puffy disks. Above a dense and geometrically thin core of dimensionless thickness $h/r \sim 0.1$, crudely resembling a classic thin accretion disk, a puffed-up, geometrically thick layer of lower density and $h/r \sim 1.0$ is formed. We discuss the observational properties of puffy disks, in particular the geometrical obscuration of the inner disk by the elevated puffy region at higher observing inclinations, and collimation of the radiation along the accretion disk spin axis, which may explain the apparent super-Eddington luminosity of some X-ray objects. We also present synthetic spectra of puffy disks, and show that they are qualitatively similar to those of a Comptonized thin disk. We demonstrate that the existing xspec spectral fitting models provide good fits to synthetic observations of puffy disks, but cannot correctly recover the input black hole spin. The puffy region remains optically thick to scattering; in its spectral properties the puffy disk roughly resembles that of a warm corona sandwiching the disk core. We suggest that puffy disks may correspond to X-ray binary systems of luminosities above 0.3 of the Eddington luminosity in the intermediate spectral states. △ Less

Submitted 14 July, 2022; v1 submitted 17 February, 2022; originally announced February 2022.

Journal ref: MNRAS, Volume 514, Issue 1, pp.780-789 (2022)

The Phase-I Trigger Readout Electronics Upgrade of the ATLAS Liquid Argon Calorimeters

Authors: G. Aad, A. V. Akimov, K. Al Khoury, M. Aleksa, T. Andeen, C. Anelli, N. Aranzabal, C. Armijo, A. Bagulia, J. Ban, T. Barillari, F. Bellachia, M. Benoit, F. Bernon, A. Berthold, H. Bervas, D. Besin, A. Betti, Y. Bianga, M. Biaut, D. Boline, J. Boudreau, T. Bouedo, N. Braam, M. Cano Bret , et al. (173 additional authors not shown)

Abstract: The Phase-I trigger readout electronics upgrade of the ATLAS Liquid Argon calorimeters enhances the physics reach of the experiment during the upcoming operation at increasing Large Hadron Collider luminosities. The new system, installed during the second Large Hadron Collider Long Shutdown, increases the trigger readout granularity by up to a factor of ten as well as its precision and range. Cons… ▽ More The Phase-I trigger readout electronics upgrade of the ATLAS Liquid Argon calorimeters enhances the physics reach of the experiment during the upcoming operation at increasing Large Hadron Collider luminosities. The new system, installed during the second Large Hadron Collider Long Shutdown, increases the trigger readout granularity by up to a factor of ten as well as its precision and range. Consequently, the background rejection at trigger level is improved through enhanced filtering algorithms utilizing the additional information for topological discrimination of electromagnetic and hadronic shower shapes. This paper presents the final designs of the new electronic elements, their custom electronic devices, the procedures used to validate their proper functioning, and the performance achieved during the commissioning of this system. △ Less

Submitted 16 May, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 56 pages, 41 figures, 6 tables

Journal ref: 2022 JINST 17 P05024

Probing Plasma Physics with Spectral Index Maps of Accreting Black Holes on Event Horizon Scales

Authors: Angelo Ricarte, Charles Gammie, Ramesh Narayan, Ben S. Prather

Abstract: The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of the supermassive black holes at the centre of our galaxy and at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. In this work, we predict spectral index… ▽ More The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of the supermassive black holes at the centre of our galaxy and at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. In this work, we predict spectral index maps of both M87* and Sgr A* by applying the general relativistic radiative transfer (GRRT) code IPOLE to a suite of general relativistic magnetohydrodynamic (GRMHD) simulations. We analytically show that the spectral index increases with increasing magnetic field strength, electron temperature, and optical depth. Consequently, spectral index maps grow more negative with increasing radius in almost all models, since all of these quantities tend to be maximised near the event horizon. Additionally, photon ring geodesics exhibit more positive spectral indices, since they sample the innermost regions of the accretion flow with the most extreme plasma conditions. Spectral index maps are sensitive to highly uncertain plasma heating prescriptions (the electron temperature and distribution function). However, if our understanding of these aspects of plasma physics can be tightened, even the spatially unresolved spectral index around 230 GHz can be used to discriminate between models. In particular, Standard and Normal Evolution (SANE) flows tend to exhibit more negative spectral indices than Magnetically Arrested Disk (MAD) flows due to differences in the characteristic magnetic field strength and temperature of emitting plasma. △ Less

Submitted 21 December, 2022; v1 submitted 4 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS

Building a weak shockwave from linear modes

Authors: Antoine Bret, Ramesh Narayan

Abstract: In shockwave theory, the density, velocity and pressure jumps are derived from the conservation equations. Here, we address the physics of a weak shock the other way around. We first show that the density profile of a weak shockwave in a fluid can be expressed as a sum of linear acoustic modes. The shock so built propagates at the speed of sound and matter is exactly conserved at the front crossin… ▽ More In shockwave theory, the density, velocity and pressure jumps are derived from the conservation equations. Here, we address the physics of a weak shock the other way around. We first show that the density profile of a weak shockwave in a fluid can be expressed as a sum of linear acoustic modes. The shock so built propagates at the speed of sound and matter is exactly conserved at the front crossing. Yet, momentum and energy are only conserved up to order 0 in powers of the shock amplitude. The density, velocity and pressure jumps are similar to those of a fluid shock, and an equivalent Mach number can be defined. A similar process is possible in magnetohydrodynamic. Yet, such a decomposition is found impossible for collisionless shocks due to the dispersive nature of ion acoustic waves. Weakly nonlinear corrections to their frequency do not solve the problem. Weak collisionless shocks could be inherently nonlinear, non-amenable to any linear superposition. Or they could be nonexistent, as hinted by recent works. △ Less

Submitted 22 November, 2021; v1 submitted 19 November, 2021; originally announced November 2021.

Comments: 8 pages, 2 figures, to appear in Journal of Plasma Physics

Journal ref: J. Plasma Phys. 88 (2022) 905880108

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

Applications and Techniques for Fast Machine Learning in Science

Authors: Allison McCarn Deiana, Nhan Tran, Joshua Agar, Michaela Blott, Giuseppe Di Guglielmo, Javier Duarte, Philip Harris, Scott Hauck, Mia Liu, Mark S. Neubauer, Jennifer Ngadiuba, Seda Ogrenci-Memik, Maurizio Pierini, Thea Aarrestad, Steffen Bahr, Jurgen Becker, Anne-Sophie Berthold, Richard J. Bonventre, Tomas E. Muller Bravo, Markus Diefenthaler, Zhen Dong, Nick Fritzsche, Amir Gholami, Ekaterina Govorkova, Kyle J Hazelwood , et al. (62 additional authors not shown)

Abstract: In this community review report, we discuss applications and techniques for fast machine learning (ML) in science -- the concept of integrating power ML methods into the real-time experimental data processing loop to accelerate scientific discovery. The material for the report builds on two workshops held by the Fast ML for Science community and covers three main areas: applications for fast ML ac… ▽ More In this community review report, we discuss applications and techniques for fast machine learning (ML) in science -- the concept of integrating power ML methods into the real-time experimental data processing loop to accelerate scientific discovery. The material for the report builds on two workshops held by the Fast ML for Science community and covers three main areas: applications for fast ML across a number of scientific domains; techniques for training and implementing performant and resource-efficient ML algorithms; and computing architectures, platforms, and technologies for deploying these algorithms. We also present overlapping challenges across the multiple scientific domains where common solutions can be found. This community report is intended to give plenty of examples and inspiration for scientific discovery through integrated and accelerated ML solutions. This is followed by a high-level overview and organization of technical advances, including an abundance of pointers to source material, which can enable these breakthroughs. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 66 pages, 13 figures, 5 tables

Report number: FERMILAB-PUB-21-502-AD-E-SCD

Journal ref: Front. Big Data 5, 787421 (2022)

Faraday depolarization and induced circular polarization by multi-path propagation with application to FRBs

Authors: Paz Beniamini, Pawan Kumar, Ramesh Narayan

Abstract: We describe how the observed polarization properties of an astronomical object are related to its intrinsic polarization properties and the finite temporal and spectral resolutions of the observing device. Moreover, we discuss the effect that a scattering screen, with non-zero magnetic field, between the source and observer has on the observed polarization properties. We show that the polarization… ▽ More We describe how the observed polarization properties of an astronomical object are related to its intrinsic polarization properties and the finite temporal and spectral resolutions of the observing device. Moreover, we discuss the effect that a scattering screen, with non-zero magnetic field, between the source and observer has on the observed polarization properties. We show that the polarization properties are determined by the ratio of observing bandwidth and coherence bandwidth of the scattering screen and the ratio of temporal resolution of the instrument and the variability time of screen, as long as the length over which the Faraday rotation induced by the screen changes by $\simπ$ is smaller than the size of the screen visible to the observer. We describe the conditions under which a source that is 100\% linearly polarized intrinsically might be observed as partially depolarized, and how the source's temporal variability can be distinguished from the temporal variability induced by the scattering screen. In general, linearly polarized waves passing through a magnetized scattering screen can develop a significant circular polarization. We apply the work to the observed polarization properties of a few fast radio bursts (FRBs), and outline potential applications to pulsars. △ Less

Submitted 28 December, 2021; v1 submitted 30 September, 2021; originally announced October 2021.

Comments: 15 pages, 9 figures. Accepted for publication in MNRAS

Jets in Magnetically Arrested Hot Accretion Flows: Geometry, Power and Black Hole Spindown

Authors: Ramesh Narayan, Andrew Chael, Koushik Chatterjee, Angelo Ricarte, Brandon Curd

Abstract: We present the results of nine simulations of radiatively-inefficient magnetically arrested disks (MADs) across different values of the black hole spin parameter $a_*$: $-0.9$, $-0.7$, $-0.5$, $-0.3$, 0, 0.3, 0.5, 0.7, and 0.9. Each simulation was run up to $t \gtrsim 100,000\,GM/c^3$ to ensure disk inflow equilibrium out to large radii. We find that the saturated magnetic flux level, and conseque… ▽ More We present the results of nine simulations of radiatively-inefficient magnetically arrested disks (MADs) across different values of the black hole spin parameter $a_*$: $-0.9$, $-0.7$, $-0.5$, $-0.3$, 0, 0.3, 0.5, 0.7, and 0.9. Each simulation was run up to $t \gtrsim 100,000\,GM/c^3$ to ensure disk inflow equilibrium out to large radii. We find that the saturated magnetic flux level, and consequently also jet power, of MAD disks depends strongly on the black hole spin, confirming previous results. Prograde disks saturate at a much higher relative magnetic flux and have more powerful jets than their retrograde counterparts. MADs with spinning black holes naturally launch jets with generalized parabolic profiles whose widths vary as a power of distance from the black hole. For distances up to $100\;GM/c^2$, the power-law index is $k \approx 0.27-0.42$. There is a strong correlation between the disk-jet geometry and the dimensionless magnetic flux, resulting in prograde systems displaying thinner equatorial accretion flows near the black hole and wider jets, compared to retrograde systems. Prograde and retrograde MADs also exhibit different trends in disk variability: accretion rate variability increases with increasing spin for $a_*>0$ and remains almost constant for $a_*\lesssim 0$, while magnetic flux variability shows the opposite trend. Jets in the MAD state remove more angular momentum from black holes than is accreted, effectively spinning down the black hole. If powerful jets from MAD systems in Nature are persistent, this loss of angular momentum will notably reduce the black hole spin over cosmic time. △ Less

Submitted 27 January, 2022; v1 submitted 27 August, 2021; originally announced August 2021.

Comments: Accepted by MNRAS 2022 January 26. Received 2022 January 26; in original form 2021 August 27

Toward determining the number of observable supermassive black hole shadows

Authors: Dominic W. Pesce, Daniel C. M. Palumbo, Ramesh Narayan, Lindy Blackburn, Sheperd S. Doeleman, Michael D. Johnson, Chung-Pei Ma, Neil M. Nagar, Priyamvada Natarajan, Angelo Ricarte

Abstract: We present estimates for the number of shadow-resolved supermassive black hole (SMBH) systems that can be detected using radio interferometers, as a function of angular resolution, flux density sensitivity, and observing frequency. Accounting for the distribution of SMBHs across mass, redshift, and accretion rate, we use a new semi-analytic spectral energy distribution model to derive the number o… ▽ More We present estimates for the number of shadow-resolved supermassive black hole (SMBH) systems that can be detected using radio interferometers, as a function of angular resolution, flux density sensitivity, and observing frequency. Accounting for the distribution of SMBHs across mass, redshift, and accretion rate, we use a new semi-analytic spectral energy distribution model to derive the number of SMBHs with detectable and optically thin horizon-scale emission. We demonstrate that (sub)millimeter interferometric observations with ${\sim}0.1$ $μ$as resolution and ${\sim}1$ $μ$Jy sensitivity could access ${>}10^6$ SMBH shadows. We then further decompose the shadow source counts into the number of black holes for which we could expect to observe the first- and second-order lensed photon rings. Accessing the bulk population of first-order photon rings requires ${\lesssim}2$ $μ$as resolution and ${\lesssim}0.5$ mJy sensitivity, while doing the same for second-order photon rings requires ${\lesssim}0.1$ $μ$as resolution and ${\lesssim}5$ $μ$Jy sensitivity. Our model predicts that with modest improvements to sensitivity, as many as $\sim$5 additional horizon-resolved sources should become accessible to the current Event Horizon Telescope (EHT), while a next-generation EHT observing at 345 GHz should have access to ${\sim}$3 times as many sources. More generally, our results can help guide enhancements of current arrays and specifications for future interferometric experiments that aim to spatially resolve a large population of SMBH shadows or higher-order photon rings. △ Less

Submitted 10 January, 2022; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 37 pages, 18 figures, published in ApJ

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.

Comments: 15 pages, 4 figures, published in PRD on May 19

Brightness Asymmetry of Black Hole Images as a Probe of Observer Inclination

Authors: Lia Medeiros, Chi-Kwan Chan, Ramesh Narayan, Feryal Ozel, Dimitrios Psaltis

Abstract: The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which show a ring-like emission structure with the South side only slightly brighter than the North side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17 deg for M87), whi… ▽ More The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which show a ring-like emission structure with the South side only slightly brighter than the North side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17 deg for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counter examples, we show that low-inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black-hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 10 pages, 11 figures, submitted to ApJ

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, 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 , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

Black Hole Magnetic Fields and Their Imprint on Circular Polarization Images

Authors: Angelo Ricarte, Richard Qiu, Ramesh Narayan

Abstract: The circular polarization of black hole accretion flows can encode properties of the underlying magnetic field structure. Using general relativistic magnetohydrodynamics (GRMHD) simulations, we study the imprint of magnetic field geometry on circular polarization images potentially observable by the Event Horizon Telescope (EHT). We decompose images into the different mechanisms that generate circ… ▽ More The circular polarization of black hole accretion flows can encode properties of the underlying magnetic field structure. Using general relativistic magnetohydrodynamics (GRMHD) simulations, we study the imprint of magnetic field geometry on circular polarization images potentially observable by the Event Horizon Telescope (EHT). We decompose images into the different mechanisms that generate circular polarization in these models, which are sensitive to both the line of sight direction and twist of the magnetic field. In these models, a stable sign of the circular polarization over time, as observed for several sources, can be attributed to a stability of these properties. We illustrate how different aspects of a generic helical magnetic field geometry become imprinted on a circular polarization image. We also identify novel effects of light bending that affect the circular polarization image on event horizon scales. One consequence is the sign flipping of successive photon rings in face-on systems, which if observable and uncorrupted by Faraday rotation, can directly encode the handedness of the approaching magnetic field. △ Less

Submitted 22 April, 2021; originally announced April 2021.

Comments: 20 pages, 19 figures, in review with MNRAS

Light echos and coherent autocorrelations in a black hole spacetime

Authors: Paul M. Chesler, Lindy Blackburn, Sheperd S. Doeleman, Michael D. Johnson, James M. Moran, Ramesh Narayan, Maciek Wielgus

Abstract: The Event Horizon Telescope recently produced the first images of a black hole. These images were synthesized by measuring the coherent correlation function of the complex electric field measured at telescopes located across the Earth. This correlation function corresponds to the Fourier transform of the image under the assumption that the source emits spatially incoherent radiation. However, blac… ▽ More The Event Horizon Telescope recently produced the first images of a black hole. These images were synthesized by measuring the coherent correlation function of the complex electric field measured at telescopes located across the Earth. This correlation function corresponds to the Fourier transform of the image under the assumption that the source emits spatially incoherent radiation. However, black holes differ from standard astrophysical objects: in the absence of absorption and scattering, an observer sees a series of increasingly demagnified echos of each emitting location. These echos correspond to rays that orbit the black hole one or more times before reaching the observer. This multi-path propagation introduces spatial and temporal correlations into the electric field that encode properties of the black hole, irrespective of intrinsic variability. We explore the coherent temporal autocorrelation function measured at a single telescope. Specifically, we study the simplified toy problem of scalar field correlation functions $\langle Ψ(t) Ψ(0) \rangle$ sourced by fluctuating matter located near a Schwarzschild black hole. We find that the correlation function is peaked at times equal to integer multiples of the photon orbit period; the corresponding power spectral density vanishes like $λ/r_{\rm g}$ where $r_{\rm g} = G M / c^{2}$ is the gravitational radius of the black hole and $λ$ is the wavelength of radiation observed. For supermassive black holes observed at millimeter wavelengths, the power in echos is suppressed relative to direct emission by $\sim 10^{-13} λ_{\rm mm}/M_{6}$, where $λ_{\rm mm} = λ/(1\,{\rm mm})$ and $M_6 = M/(10^6 M_\odot)$. Consequently, detecting multi-path propagation near a black hole using the coherent electric field autocorrelation is infeasible with current technology. △ Less

Submitted 21 December, 2020; originally announced December 2020.

Comments: 15 pages, 5 figures