Search

The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed: the main science targets, Sgr A* and M87 along with various calibrators. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-K\"onigl (BK) jet model. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Then, we studied the dependences of the VLBI core flux density, size, and brightness temperature on the frequency measured in the AGN host frame. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Our results indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations.

Using the Event Horizon Telescope (EHT), the gravitationally lensed rings around the supermassive black holes (SMBHs) in Messier 87 (M87) and Sagittarius A* (Sgr A*) have now been successfully imaged at a resolution under 10 gravitational radii (R$_{\rm g}$ $ = \rm{GM/c^2}$). To expand studies beyond M87 and Sgr A*, we have constructed the Event Horizon and Environs (ETHER) sample, a comprehensive database encompassing approximately 3.15 million SMBH mass estimates, $\sim$ 20,000 Very-Long Baseline Interferometry (VLBI) radio flux densities, and $\sim$ 36,000 hard X-ray flux densities. This database is designed to identify and optimize target selection for the EHT and its upgrades on the ground and in space. We have identified a Gold Sample (GS) of nearby low-luminosity Active Galactic Nuclei (AGNs) within it that are ideal for studying jet bases and potentially imaging black hole shadows. We observed 27 of these AGNs using the EHT from 2022 to 2024, providing an opportunity to resolve and image accretion flows and jets at resolutions of $\leq$ 100 R$_{\rm g}$. Only a few SMBHs have sufficiently high enough flux density to be imaged at scales of $\leq$ 50 R$_{\rm g}$ with the present EHT. Among these are M87, Sgr A*, NGC4594 (Sombrero/M104), NGC4261, and NGC4374 (Messier 84/M84). Of these, NGC4261, Sombrero, and M84 have been observed and/or are scheduled for deep imaging with EHT+ALMA from 2023 to 2025. Sombrero, NGC4261, M84, NGC4278, and NGC5232 are clearly detected in our EHT+ALMA observations in 2022, indicating that the 230 GHz flux density from the accretion flows is significantly high. Ongoing imaging of the ETHER GS will enable measurements of black hole mass and spin, help constrain General Relativity, and enrich our understanding of jet launching and accretion inflows across a broad multi-parameter space, including black hole mass, spin, accretion rate, and orientation., Comment: 9 pages, 6 figures, 1 table, published in Proceedings of the 16th EVN Symposium, Ed. E. Ros, P. Benke, S.A. Dzib, I. Rottmann, & J.A. Zensus, Bonn: Max-Planck-Institut f\"ur Radioastronomie, 2024, pages 75-84, https://cloud.mpifr-bonn.mpg.de/index.php/s/BkX2CC2Xjn2aKR4

We investigate how magnetic field variations around accreting black holes on event horizon scales affect the morphology of magnetically-driven jet on larger scales. By performing radiative transfer calculations on general relativistic magnetohydrodynamics simulations, we find that temporal variation in the magnetic flux on the event horizon and the jet power are imprinted on the variability of jet width up to several hundred gravitational radii. When the magnetic flux around the black hole drops and then rises, the jet initially narrows or becomes truncated, then widens, creating a thin-thick pattern that propagates down the jet. This suggests that extended jet observations can provide a history record of horizon-scale magnetic field dynamics, and conversely, upcoming changes in the jet image can be predicted from direct observation of the magnetized accreting plasma near the black hole. Furthermore, the pattern of jet width variations shows acceleration up to the relativistic regime as it moves away from the black hole, aligning with plasma bulk motion. We also find in time-averaged images that both the bulk plasma motion and magnetic field configuration in the jet-launching region, which are sensitive to black hole spin, shape diverse features through relativistic beaming and aberration. Higher black hole spins result in more poloidal bulk motion and toroidal magnetic fields, leading to more symmetric jet images and linear polarization patterns. These results suggest a new method for testing the magnetically arrested disk model and the Blandford-Znajek process, and for determining the black hole spin through observations bridging horizon and jet-launching scales., Comment: 19 pages, 16 figures, submitted to ApJ, movie available at https://youtu.be/FwiKspIXjZk

Primordial black holes can be the entirety of the dark matter in a broad, approximately five-orders-of-magnitude-wide mass range, the ``asteroid mass range'', between $10^{-16}\ M_{\rm Sun}$ -- where constraints originate from evaporation -- and $10^{-11}\ M_{\rm Sun}$ -- from microlensing. A direct detection in this mass range is very challenging with any known observational or experimental methods. Here we point out that, unlike previously asserted in the literature, a transient gravitational wave signal from the inspiral phase of light black hole mergers is in principle detectable with current and future high-frequency gravitational wave detectors, including but not limited to ADMX. The largest detection rates are associated with binaries from non-monochromatic mass functions in early-formed three-body systems., Comment: 24 pages, 9 figures

In this manuscript, we obtain sharp and improved regularity estimates for weak solutions of weighted quasilinear elliptic models of Hardy-H\'{e}non-type, featuring an explicit regularity exponent depending only on universal parameters. Our approach is based on geometric tangential methods and uses a refined oscillation mechanism, compactness, and scaling techniques. In some specific scenarios, we establish higher regularity estimates and non-degeneracy properties, providing further geometric insights into such solutions. Our regularity estimates both enhance and, to some extent, extend the results arising from the $C^{p^{\prime}}$ conjecture for the $p$-Laplacian with a bounded source term. As applications of our results, we address some Liouville-type results for our class of equations. Finally, our results are noteworthy, even in the simplest model case governed by the $p$-Laplacian with regular coefficients: $$ \mathrm{div}\left( |\nabla u|^{p-2}\mathfrak{A}(|x|) \nabla u\right) = |x|^{\alpha}u_+^m(x) \quad \text{in} \quad B_1 $$ under suitable assumptions on the data, with possibly singular weight $\mathfrak{h}(|x|) = |x|^{\alpha}$, which includes the Matukuma and Batt-Faltenbacher-Horst's equations as toy models.

Using the semi-analytic model Serotina, we investigate the cosmic spin evolution of supermassive black holes incorporating recent results from general relativistic magnetohydrodynamics simulations of spin-down from relativistic jets. We compare several variations of our model with compiled black hole spin measurements derived from X-ray reflection spectroscopy, correcting for a bias arising from the spin-dependent radiative efficiency of accretion flows. We show that the observed spin distribution is in agreement with a model that includes jet-driven spin-down, a key mechanism that acts to modulate spins across cosmic time at both high and very low specific accretion rates. The data also clearly prefer models with coherent accretion over models in which accretion disks rapidly switch from prograde to retrograde. We further predict spin distributions accessible via spatially resolved event horizons by the next-generation Event Horizon Telescope (ngEHT) and Black Hole Explorer (BHEX), as well as gravitational waves by the Laser Interferometer Space Antenna (LISA), each of which offer unique and distinct windows into the population of spinning black holes. Jet-driven spin-down is most strongly imprinted on the abundance of very highly spinning objects in our model. In addition, we show that the spin distribution sampled by LISA events may contain a signature of the natal spin distribution of heavy seeds, but not of light seeds, offering additional discrimination between these seeding pathways. Spin distributions from these future observed samples can be used to constrain the detailed physical properties of the accretion flow on horizon scales around supermassive black holes., Comment: 23 pages, 11 figures, 4 tables. Submitted to ApJ

The first very long baseline interferometry (VLBI) detections at 870$\mu$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$\lambda$ corresponding to an angular resolution, or fringe spacing, of 19$\mu$as. The Allan deviation of the visibility phase at 870$\mu$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$\mu$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time., Comment: Corresponding author: S. Doeleman

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)., Comment: 21 pages, 9 figures, 1 table, submitted to ApJ

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., Comment: Proceedings for SPIE Astronomical Telescopes and Instrumentation

Although supermassive black holes (SMBHs) reside in the heart of virtually every massive galaxy, it remains debated whether dwarf galaxies also commonly host SMBHs. Because low-mass galaxies may retain a memory of the assembly history of their black holes, probing the black hole occupation fraction of local dwarf galaxies might offer insights into the growth and seeding mechanisms of the first black holes. In this work, we exploit the Western half of the eROSITA all-sky survey (covering $20,000~\rm{deg^2}$) and compile a catalog of accreting SMBHs in local ($D<200$~Mpc) dwarf galaxies. After cleaning our sample from cosmic X-ray background sources, X-ray binaries, and ultraluminous X-ray sources, we identify 74 AGN-dwarf galaxy pairs. Using this large and uniform sample, we derive a luminosity function of dwarf galaxy AGN, fitting it with a power law function and obtaining ${\rm d}N/{\rm d}L_{\rm X} = (15.9\pm2.2)\times L_{\rm X}^{-1.63\pm0.05}$. Measuring the offset between the centroid of dwarf galaxies and the X-ray sources, we find that about $50\%$ of the AGN are likely off-nuclear, in agreement with theoretical predictions. We also compare the black hole-to-stellar mass relation of the AGN in our sample with the local and high-redshift relations, finding that our sources better adhere to the former. This suggests that local AGN across different mass scales underwent a similar growth history. Finally, we compare our sources with semi-analytical models: while our sample is too shallow to distinguish between different seeding models, it favors a growth mechanism linked to the star-formation rate of the host galaxy., Comment: 17 pages, 12 figures, 1 table. Submitted to ApJ. Comments are welcome

This paper studies Schauder theory to transmission problems modelled by fully nonlinear uniformly elliptic equations of second order. We focus on operators F that fails to be concave or convex in the space of symmetric matrices. In a first scenario, it is considered that F enjoys a small ellipticity aperture. In our second case, we study regularity results where the convexity of the superlevel (or sublevel) sets is verified, implying that the operator F is quasiconcave (or quasiconvex).

The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high energy (VHE) gamma-rays, as well as details of the individual observations and light curves. We also conduct phenomenological modelling to investigate the basic source properties. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV has more than doubled within a period of about 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and emphasises the need for combined image and spectral modelling., Comment: 46 pages, 23 figures, accepted by Astronomy & Astrophysics on August. 29, 2024

In this work, we establish universal moduli of continuity for viscosity solutions to fully nonlinear elliptic equations with oblique boundary conditions, whose general model is given by $$ \left\{ \begin{array}{rcl} F(D^2u,x) &=& f(x) \quad \mbox{in} \,\, \Omega\\ \beta(x) \cdot Du(x) + \gamma(x) \, u(x)&=& g(x) \quad \mbox{on} \,\, \partial \Omega. \end{array} \right. $$ Such regularity estimates are achieved by exploring the integrability properties of $f$ based on different scenarios, making a $\text{VMO}$ assumption on the coefficients of $F$, and by considering suitable smoothness properties on the boundary data $\beta, \gamma$ and $g$. Particularly, we derive sharp estimates for borderline cases where $f \in L^n(\Omega)$ and $f\in p-\textrm{BMO}(\Omega)$. Additionally, for source terms in $L^p(\Omega)$, for $p \in (n, \infty)$, we obtain sharp gradient estimates. Finally, we also address Schauder-type estimates for convex/concave operators and suitable H\"{o}lder data.

3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $\nu_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored., Comment: 15 pages, 6 figures, published in A&A

SS 433 is a microquasar, a stellar binary system with collimated relativistic jets. We observed SS 433 in gamma rays using the High Energy Stereoscopic System (H.E.S.S.), finding an energy-dependent shift in the apparent position of the gamma-ray emission of the parsec-scale jets. These observations trace the energetic electron population and indicate the gamma rays are produced by inverse-Compton scattering. Modelling of the energy-dependent gamma-ray morphology constrains the location of particle acceleration and requires an abrupt deceleration of the jet flow. We infer the presence of shocks on either side of the binary system at distances of 25 to 30 parsecs and conclude that self-collimation of the precessing jets forms the shocks, which then efficiently accelerate electrons., Comment: Submitted 20th Apr. 2023, published 25th January 2024 (accepted version)

In this paper, we introduce a novel data augmentation methodology based on Conditional Progressive Generative Adversarial Networks (CPGAN) to generate diverse black hole (BH) images, accounting for variations in spin and electron temperature prescriptions. These generated images are valuable resources for training deep learning algorithms to accurately estimate black hole parameters from observational data. Our model can generate BH images for any spin value within the range of [-1, 1], given an electron temperature distribution. To validate the effectiveness of our approach, we employ a convolutional neural network to predict the BH spin using both the GRMHD images and the images generated by our proposed model. Our results demonstrate a significant performance improvement when training is conducted with the augmented dataset while testing is performed using GRMHD simulated data, as indicated by the high R2 score. Consequently, we propose that GANs can be employed as cost effective models for black hole image generation and reliably augment training datasets for other parameterization algorithms., Comment: 11 pages, 7 figures. Accepted by Monthly Notices of the Royal Astronomical Society Journal

The nature and origin of supermassive black holes (SMBHs) remain an open matter of debate within the scientific community. While various theoretical scenarios have been proposed, each with specific observational signatures, the lack of sufficiently sensitive X-ray observations hinders the progress of observational tests. In this white paper, we present how AXIS will contribute to solving this issue. With an angular resolution of 1.5$^{\prime\prime}$ on-axis and minimal off-axis degradation, we have designed a deep survey capable of reaching flux limits in the [0.5-2] keV range of approximately 2$\times$10$^{-18}$ \fcgs~ over an area of 0.13 deg$^2$ in approximately 7 million seconds (7 Ms). Furthermore, we have planned an intermediate depth survey covering approximately 2 deg$^2$ and reaching flux limits of about 2$\times$10$^{-17}$ \fcgs ~ in order to detect a significant number of SMBHs with X-ray luminosities (L$_X$) of approximately 10$^{42}$ \lx up to z$\sim$10. These observations will enable AXIS to detect SMBHs with masses smaller than 10$^5$ \ms, assuming Eddington-limited accretion and a typical bolometric correction for Type II AGN. AXIS will provide valuable information on the seeding and population synthesis models of SMBH, allowing for more accurate constraints on their initial mass function (IMF) and accretion history from z$\sim$0-10. To accomplish this, AXIS will leverage the unique synergy of survey telescopes such as JWST, Roman, Euclid, LSST, and the new generation of 30m class telescopes. These instruments will provide optical identification and redshift measurements, while AXIS will discover the smoking gun of nuclear activity, particularly in the case of highly obscured AGN or peculiar UV spectra as predicted and recently observed in the early Universe., Comment: This White Paper is part of a series commissioned for the AXIS Probe Concept Mission. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at http://axis.astro.umd.edu with a mission overview at arXiv:2311.00780

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., Comment: 8 pages with 5 pages of Appendices, 4 figures, accepted to ApJ Letters

The objective of this work is to expand upon previous works, considering socially acceptable behaviours within robot navigation and interaction, and allow a robot to closely approach static and dynamic individuals or groups. The space models developed in this dissertation are adaptive, that is, capable of changing over time to accommodate the changing circumstances often existent within a social environment. The space model's parameters' adaptation occurs with the end goal of enabling a close interaction between humans and robots and is thus capable of taking into account not only the arrangement of the groups, but also the basic characteristics of the robot itself. This work also further develops a preexisting approach pose estimation algorithm in order to better guarantee the safety and comfort of the humans involved in the interaction, by taking into account basic human sensibilities. The algorithms are integrated into ROS's navigation system through the use of the $costmap2d$ and the $move\_base$ packages. The space model adaptation is tested via comparative evaluation against previous algorithms through the use of datasets. The entire navigation system is then evaluated through both simulations (static and dynamic) and real life situations (static). These experiments demonstrate that the developed space model and approach pose estimation algorithms are capable of enabling a robot to closely approach individual humans and groups, while maintaining considerations for their comfort and sensibilities.

Horizon-scale observations of the jetted active galactic nucleus M87 are compared with simulations spanning a broad range of dissipation mechanisms and plasma content in three-dimensional general relativistic flows around spinning black holes. Observations of synchrotron radiation from radio to X-ray frequencies can be compared with simulations by adding prescriptions specifying the relativistic electron-plus-positron distribution function and associated radiative transfer coefficients. A suite of time-varying simulations with various spins, plasma magnetizations and turbulent heating and equipartition-based emission prescriptions (and piecewise combinations thereof) is chosen to represent distinct possibilities for the M87 jet/accretion flow/black hole (JAB) system. Simulation jet morphology, polarization and variation are then "observed" and compared with real observations to infer the rules that govern the polarized emissivity. Our models support several possible spin/emission model/plasma composition combinations supplying the jet in M87, whose black hole shadow has been observed down to the photon ring at 230 GHz by the Event Horizon Telescope (EHT). Net linear polarization and circular polarization constraints favor magnetically arrested disk (MAD) models whereas resolved linear polarization favors standard and normal evolution (SANE) in our parameter space. We also show that some MAD cases dominated by intrinsic circular polarization have near-linear V/I dependence on unpaired electron or positron content while SANE polarization exhibits markedly greater positron-dependent Faraday effects - future probes of the SANE/MAD dichotomy and plasma content with the EHT. This is the second work in a series also applying the "observing" simulations methodology to near-horizon regions of supermassive black holes in Sgr A* and 3C 279., Comment: 27 pages, 27 figures

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 ($\lambda$=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., Comment: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

The recent Chandra-JWST discovery of a quasar in the z = 10.1 galaxy UHZ1 reveals that accreting supermassive black holes (SMBHs) were already in place 470 million years after the Big Bang. The Chandra X-ray source detected in UHZ1 is a Compton-thick quasar with a bolometric luminosity of $L_{\rm bol}\sim5\times10^{45}\ \rm{erg\ s^{-1}},$ which corresponds to an estimated BH mass of $\sim4\times10^7 \ \rm{M_{\odot}}$ assuming accretion at the Eddington rate. JWST NIRCAM and NIRSpec data yield a stellar mass estimate for UHZ1 comparable to its BH mass. These characteristics are in excellent agreement with prior theoretical predictions for a unique class of transient, high-redshift objects, Over-massive Black Hole Galaxies [OBGs] by Natarajan et al. that harbor a heavy initial black hole seed that likely formed from the direct collapse of the gas. Based on the excellent agreement between the observed multi-wavelength properties of UHZ1 with theoretical model template predictions, suggests that UHZ1 is the first detected OBG candidate. Our assertion rests on multiple lines of concordant evidence between model predictions and the following observed properties of UHZ1: its X-ray detection and the estimated ratio of the X-ray flux to the IR flux that is consistent with theoretical expectations for a heavy initial BH seed; its high measured redshift of z = 10.1, as predicted for the transient OBG stage (9 < z< 12); the amplitude and shape of the detected JWST Spectral Energy Distribution (SED) between 1 - 5 microns, which is in very good agreement with simulated template SEDs for OBGs; and the extended JWST morphology of UHZ1 that is suggestive of a recent merger, also expected for the formation of transient OBGs. As the first OBG candidate, UHZ1 provides compelling evidence for the formation of heavy initial seeds from direct collapse in the early Universe., Comment: 8 pages, 4 figures, submitted to ApJ, under review

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., Comment: 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

In July 2021, PKS 1510-089 exhibited a significant flux drop in the high-energy gamma-ray (by a factor 10) and optical (by a factor 5) bands and remained in this low state throughout 2022. Similarly, the optical polarization in the source vanished, resulting in the optical spectrum being fully explained through the steady flux of the accretion disk and the broad-line region. Unlike the aforementioned bands, the very-high-energy gamma-ray and X-ray fluxes did not exhibit a significant flux drop from year to year. This suggests that the steady-state very-high-energy gamma-ray and X-ray fluxes originate from a different emission region than the vanished parts of the high-energy gamma-ray and optical jet fluxes. The latter component has disappeared through either a swing of the jet away from the line-of-sight or a significant drop in the photon production efficiency of the jet close to the black hole. Either change could become visible in high-resolution radio images., Comment: 18 pages, 4 figures, 6 tables; accepted for publication in ApJ Letters; corresponding authors: Joleen Barnard, Markus Boettcher, Hester Schutte, Michael Zacharias

The occupation fraction of massive black holes (MBHs) in dwarf galaxies offers interesting insights into initial black hole seeding mechanisms and their mass assembly history, though disentangling these two effects remains challenging. Using the {\sc Romulus} cosmological simulations we examine the impact of environment on the occupation fraction of MBHs in low mass galaxies. Unlike most modern cosmological simulations, {\sc Romulus} seeds MBHs based on local gas properties, selecting dense ($n>3$ cm$^{-3}$), pristine ($Z<3e-4Z_{\odot}$), and rapidly collapsing regions in the early Universe as sites to host MBHs without assuming anything about MBH occupation as a function of galaxy stellar mass, or halo mass, {\it a priori}. The simulations predict that dwarf galaxies with M$_{\star}<10^9$ M$_{\odot}$ in cluster environments are $\sim2$ times more likely to host a MBH compared to those in the field. The predicted occupation fractions are remarkably consistent with those of nuclear star clusters. Across cluster and field environments, dwarf galaxies with earlier formation times are more likely to host a MBH. While the MBH occupation function is similar between cluster and field environments at high redshift ($z>3$), a difference arises as late-forming dwarfs -- which do not exist in the cluster environment -- begin to dominate in the field and pull the MBH occupation fraction down for low mass galaxies. Additionally, prior to in-fall some cluster dwarfs are similar to progenitors of massive, isolated galaxies, indicating that they might have grown to higher masses had they not been impeded by the cluster environment. While the population of MBHs in dwarf galaxies is already widely understood to be important for understanding MBH formation, this work demonstrates that environmental dependence is important to consider as future observations search for low mass black holes in dwarf galaxies., Comment: 16 pages, 7 figures, Accepted to the Open Journal of Astrophysics

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., Comment: 26 pages, 21 Figures

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., Comment: 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)

Geminga is an enigmatic radio-quiet gamma-ray pulsar located at a mere 250 pc distance from Earth. Extended very-high-energy gamma-ray emission around the pulsar was discovered by Milagro and later confirmed by HAWC, which are both water Cherenkov detector-based experiments. However, evidence for the Geminga pulsar wind nebula in gamma rays has long evaded detection by imaging atmospheric Cherenkov telescopes (IACTs) despite targeted observations. The detection of gamma-ray emission on angular scales > 2 deg poses a considerable challenge for the background estimation in IACT data analysis. With recent developments in understanding the complementary background estimation techniques of water Cherenkov and atmospheric Cherenkov instruments, the H.E.S.S. IACT array can now confirm the detection of highly extended gamma-ray emission around the Geminga pulsar with a radius of at least 3 deg in the energy range 0.5-40 TeV. We find no indications for statistically significant asymmetries or energy-dependent morphology. A flux normalisation of $(2.8\pm0.7)\times10^{-12}$ cm$^{-2}$s$^{-1}$TeV$^{-1}$ at 1 TeV is obtained within a 1 deg radius region around the pulsar. To investigate the particle transport within the halo of energetic leptons around the pulsar, we fitted an electron diffusion model to the data. The normalisation of the diffusion coefficient obtained of $D_0 = 7.6^{+1.5}_{-1.2} \times 10^{27}$ cm$^2$s$^{-1}$, at an electron energy of 100 TeV, is compatible with values previously reported for the pulsar halo around Geminga, which is considerably below the Galactic average., Comment: 16 pages, 15 figures, 7 tables. Accepted for publication in Astronomy & Astrophysics

Primordial Black Holes (PBHs) are hypothetical black holes predicted to have been formed from density fluctuations in the early Universe. PBHs with an initial mass around $10^{14}-10^{15}$g are expected to end their evaporation at present times in a burst of particles and very-high-energy (VHE) gamma rays. Those gamma rays may be detectable by the High Energy Stereoscopic System (H.E.S.S.), an array of imaging atmospheric Cherenkov telescopes. This paper reports on the search for evaporation bursts of VHE gamma rays with H.E.S.S., ranging from 10 to 120 seconds, as expected from the final stage of PBH evaporation and using a total of 4816 hours of observations. The most constraining upper limit on the burst rate of local PBHs is $2000$ pc$^{-3}$ yr$^{-1}$ for a burst interval of 120 seconds, at the 95\% confidence level. The implication of these measurements for PBH dark matter are also discussed., Comment: Accepted for publication in JCAP, corresponding authors: F.Brun, J-F. Glicenstein, V. Marandon, T. Tavernier

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., Comment: Accepted for publication in ApJ

We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over three years (from October 2018 to September 2021). We infer the black hole (BH) mass ($M_{\rm BH}$) with host galaxy scaling relations, showing that the sample $M_{\rm BH}$ ranges from $10^{5.1}\,M_\odot$ to $10^{8.2}\,M_\odot$. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frame $g$-band luminosity function (LF) can be well described by a broken power-law of $\phi (L_g)\propto [(L_g / L_{\rm bk})^{0.3} + (L_g / L_{\rm bk})^{2.6}]^{-1}$, with $L_{\rm bk}=10^{43.1}\,{\rm erg\,s^{-1}}$. In the BH mass regime of $10^{5.3}\lesssim (M_{\rm BH}/M_\odot) \lesssim 10^{7.3}$, the TDE mass function follows $\phi(M_{\rm BH})\propto M_{\rm BH}^{-0.25}$, which favors a flat local BH mass function ($dn_{\rm BH}/d{\rm log}M_{\rm BH}\approx{\rm constant}$). We confirm the significant rate suppression at the high-mass end ($M_{\rm BH}\gtrsim 10^{7.5}\,M_\odot$), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass of $M_{\rm gal}\sim 10^{10}\,M_\odot$, the average optical TDE rate is $\approx 3.2\times 10^{-5}\,{\rm galaxy^{-1}\,yr^{-1}}$. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6$]\times 10^{-5}\,{\rm galaxy^{-1}\,yr^{-1}}$ in galaxies with red, green, and blue colors., Comment: Replaced following peer-review process. 38 pages, 23 figures. Accepted for publication in ApJL

In this work, we develop weighted Lorentz-Sobolev estimates for viscosity solutions of fully nonlinear elliptic equations with oblique boundary condition under weakened convexity conditions in the following configuration $F(D^{2}u, Du, u, x) = f(x)$ in $\Omega$ and $\beta\cdot Du+\gamma u=g$ on $\partial \Omega$, where $\Omega$ is a bounded domain in $\mathbb{R}^{n}$ ($n \geq 2$), under suitable assumptions on the source term f, data $\beta$, $\gamma$ and g. In addition, we obtain Lorentz-Sobolev estimates for solutions to the obstacle problem and others applications., Comment: 25 pages

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., Comment: 24 pages, 27 figures

Dwarf galaxies are promising test beds for constraining models of supermassive and intermediate-mass black holes (MBH) via their black hole occupation fraction (BHOF). Disentangling seeding from the confounding effects of mass assembly over a Hubble time is a challenging problem, that we tackle in this study with a suite of semi-analytical models (SAMs). We show how measured BHOF depends on the lowest black hole mass or AGN luminosity achieved by a survey. To tell seeding models apart, we need to detect or model all AGN brighter than $10^{37}\ \rm{erg \ s^{-1}}$ in galaxies of $M_* \sim 10^{8-10} \ \rm{M_{\odot}}$. Shallower surveys, like eRASS, cannot distinguish between seed models even with the compensation of a much larger survey volume. We show that the AMUSE survey, with its inference of the MBH population underlying the observed AGN, strongly favors heavy seed models, growing with either a power-law Eddington Ratio Distribution Function (ERDF) or one in which black hole accretion is tagged to the star-formation rate (AGN-MS). These two growth channels can then be distinguished by the AGN luminosity function at $> 10^{40}\ \rm{erg \ s^{-1}}$, with the AGN-MS model requiring more accretion than observed at z $\sim$ 0. Thus, current X-ray observations favour heavy seeds whose Eddington ratios follow a power-law distribution. The different models also predict different radio scaling relations, which we quantify using the fundamental plane of black hole activity. We close with recommendations for the design of upcoming multi-wavelength campaigns that can optimally detect MBHs in dwarf galaxies., Comment: Accepted to ApJ

In the past few years, the Event Horizon Telescope (EHT) has provided the first-ever event horizon-scale images of the supermassive black holes (BHs) (M87*) and Sagittarius A$^*$ (Sgr A*). The next-generation EHT project is an extension of the EHT array that promises larger angular resolution and higher sensitivity to the dim, extended flux around the central ring-like structure, possibly connecting the accretion flow and the jet. The ngEHT Analysis Challenges aim to understand the science extractability from synthetic images and movies to inform the ngEHT array design and analysis algorithm development. In this work, we compare the accretion flow structure and dynamics in numerical fluid simulations that specifically target M87* and Sgr A*, and were used to construct the source models in the challenge set. We consider (1) a steady-state axisymmetric radiatively inefficient accretion flow model with a time-dependent shearing hotspot, (2) two time-dependent single fluid general relativistic magnetohydrodynamic (GRMHD) simulations from the H-AMR code, (3) a two-temperature GRMHD simulation from the BHAC code, and (4) a two-temperature radiative GRMHD simulation from the KORAL code. We find that the different models exhibit remarkably similar temporal and spatial properties, except for the electron temperature, since radiative losses substantially cool down electrons near the BH and the jet sheath, signaling the importance of radiative cooling even for slowly accreting BHs such as M87*. We restrict ourselves to standard torus accretion flows, and leave larger explorations of alternate accretion models to future work., Comment: Accepted in Galaxies; 23 pages, 7 figures

We present estimates for the number of supermassive black holes (SMBHs) for which the next-generation Event Horizon Telescope (ngEHT) can identify the black hole ``shadow,'' along with estimates for how many black hole masses and spins the ngEHT can expect to constrain using measurements of horizon-resolved emission structure. Building on prior theoretical studies of SMBH accretion flows and analyses carried out by the Event Horizon Telescope (EHT) collaboration, we construct a simple geometric model for the polarized emission structure around a black hole, and we associate parameters of this model with the three physical quantities of interest. We generate a large number of realistic synthetic ngEHT datasets across different assumed source sizes and flux densities, and we estimate the precision with which our defined proxies for physical parameters could be measured from these datasets. Under April weather conditions and using an observing frequency of 230~GHz, we predict that a ``Phase 1'' ngEHT can potentially measure $\sim$50 black hole masses, $\sim$30 black hole spins, and $\sim$7 black hole shadows across the entire sky., Comment: accepted for publication in the ngEHT special issue of Galaxies

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 $\beta_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., Comment: 14 pages, 9 figures

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*., Comment: 16 pages, 8 pages

This work proposes a methodology to test phenomenologically-motivated emission processes that account for the flux and polarization distribution and global structure of the 230 GHz sources imaged by the Event Horizon Telescope (EHT): Messier (M)87* and Sagittarius (Sgr) A*. We introduce to general relativistic magnetohydrodynamic (GRMHD) simulations some novel models to bridge the largely uncertain mechanisms by which high-energy particles in jet/accretion flow/black hole (JAB) system plasmas attain billion degree temperatures and emit synchrotron radiation. The "Observing" JAB Systems methodology then partitions the simulation to apply different parametric models to regions governed by different plasma physics -- an advance over methods where one parametrization is used over simulation regions spanning thousands of gravitational radii from the central supermassive black hole. We present several classes of viewing-angle dependent morphologies, and highlight signatures of piecewise modeling and positron effects -- including a MAD/SANE dichotomy in which polarized maps appear dominated by intrinsic polarization in the MAD case and by Faraday effects in the SANE case. The library of images thus produced spans a wide range of morphologies awaiting discovery by the groundbreaking EHT instrument and its yet more sensitive, higher resolution next-generation counterpart ngEHT., Comment: 21 pages, 9 figures

While supermassive black hole masses have been cataloged across cosmic time, only a few dozen of them have robust spin measurements. By extending and improving the existing Event Horizon Telescope (EHT) array, the next-generation Event Horizon Telescope (ngEHT) will enable multifrequency, polarimetric movies on event horizon scales, which will place new constraints on the space-time and accretion flow. By combining this information, it is anticipated that the ngEHT may be able to measure tens of supermassive black hole masses and spins. In this white paper, we discuss existing spin measurements and many proposed techniques with which the ngEHT could potentially measure spins of target supermassive black holes. Spins measured by the ngEHT would represent a completely new sample of sources that, unlike pre-existing samples, would not be biased towards objects with high accretion rates. Such a sample would provide new insights into the accretion, feedback, and cosmic assembly of supermassive black holes., Comment: Submitted for Galaxies Special Issue "From Vision to Instrument: Creating a Next-Generation Event Horizon Telescope for a New Era of Black Hole Science"

The Event Horizon Telescope (EHT) Collaboration has successfully produced images of two supermassive black holes, enabling novel tests of black holes and their accretion flows on horizon scales. The EHT has so far published total intensity and linear polarization images, while upcoming images may include circular polarization, rotation measure, and spectral index, each of which reveals different aspects of the plasma and space-time. The next-generation EHT (ngEHT) will greatly enhance these studies through wider recorded bandwidths and additional stations, leading to greater signal-to-noise, orders of magnitude improvement in dynamic range, multi-frequency observations, and horizon-scale movies. In this paper, we review how each of these different observables informs us about the underlying properties of the plasma and the spacetime, and we discuss why polarimetric studies are well-suited to measurements with sparse, long-baseline coverage., Comment: Submitted for Galaxies Special Issue "From Vision to Instrument: Creating a Next-Generation Event Horizon Telescope for a New Era of Black Hole Science"

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., Comment: Submitted to ApJL. 15 pages, 8 figures

Massive black holes in the centers of galaxies today must have grown by several orders of magnitude from seed black holes formed at early times. Detecting a population of intermediate mass black holes (IMBHs) can provide constraints on these elusive BH seeds. Here we use the large volume, cosmological hydrodynamical simulation Astrid, which includes IMBH seeds and dynamical friction to investigate the population of IMBH seeds. Dynamical friction is largely inefficient at sinking and merging seed IMBHs at high-z. This leads to an extensive population (several hundred per galaxy) of wandering IMBHs in large halos at z~2. A small fraction of these IMBHs are detectable as HLXs, Hyper Luminous X-ray sources. Importantly, at z ~ 2, IMBHs mergers produce the peak of GW events. We find close to a million GW events in Astrid between z=2-3 involving seed IMBH mergers. These GW events (almost all detectable by LISA) at cosmic noon should provide strong constraints on IMBH seed models and their formation mechanisms. At the center of massive galaxies, where the number of IMBHs can be as high as 10-100, SMBH-IMBH pairs can form. These Intermediate mass ratio inspirals (IMRIs) and extreme mass ratio inspirals (EMRIs), will require the next generation of milli-muHz space-based GW interferometers to be detected. Large populations of IMBHs around massive black holes will probe their environments and MBH causal structure.

The Event Horizon Telescope (EHT) has produced images of the plasma flow around the supermassive black holes in Sgr A* and M87* with a resolution comparable to the projected size of their event horizons. Observations with the next-generation Event Horizon Telescope (ngEHT) will have significantly improved Fourier plane coverage and will be conducted at multiple frequency bands (86, 230, and 345 GHz), each with a wide bandwidth. At these frequencies, both Sgr A* and M87* transition from optically thin to optically thick. Resolved spectral index maps in the near-horizon and jet-launching regions of these supermassive black hole sources can constrain properties of the emitting plasma that are degenerate in single-frequency images. In addition, combining information from data obtained at multiple frequencies is a powerful tool for interferometric image reconstruction, since gaps in spatial scales in single-frequency observations can be filled in with information from other frequencies. Here we present a new method of simultaneously reconstructing interferometric images at multiple frequencies along with their spectral index maps. The method is based on existing Regularized Maximum Likelihood (RML) methods commonly used for EHT imaging and is implemented in the eht-imaging Python software library. We show results of this method on simulated ngEHT data sets as well as on real data from the VLBA and ALMA. These examples demonstrate that simultaneous RML multi-frequency image reconstruction produces higher-quality and more scientifically useful results than is possible from combining independent image reconstructions at each frequency., Comment: 25 pages, 15 figures. Accepted to ApJ

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., Comment: 26 pages, 20 figures, accepted for publication in ApJ

MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ~ 500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential high-energy and very-high-energy gamma-ray emitting region should be located at a distance from the black hole ranging between 10^11 and 10^13 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA., Comment: 15 pages, 5 figures. Accepted for publication in MNRAS

We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~\mu{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*.

The central region of the Milky Way is one of the foremost locations to look for dark matter (DM) signatures. We report the first results on a search for DM particle annihilation signals using new observations from an unprecedented gamma-ray survey of the Galactic Center (GC) region, ${\it i.e.}$, the Inner Galaxy Survey, at very high energies ($\gtrsim$ 100 GeV) performed with the H.E.S.S. array of five ground-based Cherenkov telescopes. No significant gamma-ray excess is found in the search region of the 2014-2020 dataset and a profile likelihood ratio analysis is carried out to set exclusion limits on the annihilation cross section $\langle \sigma v\rangle$. Assuming Einasto and Navarro-Frenk-White (NFW) DM density profiles at the GC, these constraints are the strongest obtained so far in the TeV DM mass range. For the Einasto profile, the constraints reach $\langle \sigma v\rangle$ values of $\rm 3.7\times10^{-26} cm^3s^{-1}$ for 1.5 TeV DM mass in the $W^+W^-$ annihilation channel, and $\rm 1.2 \times 10^{-26} cm^3s^{-1}$ for 0.7 TeV DM mass in the $\tau^+\tau^-$ annihilation channel. With the H.E.S.S. Inner Galaxy Survey, ground-based $\gamma$-ray observations thus probe $\langle \sigma v\rangle$ values expected from thermal-relic annihilating TeV DM particles., Comment: Accepted in Phys. Rev. Lett., includes Supplemental materials. 5+10 pages, 2+8 figures, 4 tables

In this work we derive global estimates for viscosity solutions to fully nonlinear elliptic equations under relaxed structural assumptions on the governing operator which are weaker than convexity and oblique boundary conditions and under suitable assumptions on the dat. Our approach makes use of geometric tangential methods, which consists of importing "fine regularity estimates" from a limiting profile, i.e., the Recession operator, associated with the original second order one via compactness and stability procedures. As a result, we devote a special attention to the borderline scenario. In such a setting, we prove that solutions enjoy BMO type estimates for their second derivatives. In the end, as another application of our findings, we obtain Hessian estimates to obstacle type problems under oblique boundary conditions and no convexity assumptions, which may have their own mathematical interest. A density result in a suitable class of viscosity solutions will be also addressed.

The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.