Search for: All records

Abstract We introduceMahakala, aPython-based, modular, radiative ray-tracing code for curved spacetimes. We employ Google’sJAXframework for accelerated automatic differentiation, which can efficiently compute Christoffel symbols directly from the metric, allowing the user to easily and quickly simulate photon trajectories through non-Kerr spacetimes.JAXalso enablesMahakalato run in parallel on both CPUs and GPUs.Mahakalanatively uses the Cartesian Kerr–Schild coordinate system, which avoids numerical issues caused by the pole in spherical coordinate systems. We demonstrateMahakala’s capabilities by simulating 1.3 mm wavelength images (the wavelength of Event Horizon Telescope observations) of general relativistic magnetohydrodynamic simulations of low-accretion rate supermassive black holes. The modular nature ofMahakalaallows us to quantitatively explore how different regions of the flow influence different image features. We show that most of the emission seen in 1.3 mm images originates close to the black hole and peaks near the photon orbit. We also quantify the relative contribution of the disk, forward jet, and counterjet to 1.3 mm images. 

Abstract Models of the resolved Event Horizon Telescope (EHT) sources Sgr A* and M87* are constrained by observations at multiple wavelengths, resolutions, polarizations, and time cadences. In this paper, we compare unresolved circular polarization (CP) measurements to a library of models, where each model is characterized by a distribution of CP over time. In the library, we vary the spin of the black hole, the magnetic field strength at the horizon (i.e., both SANE and magnetically arrested disk or MAD models), the observer inclination, a parameter for the maximum ion–electron temperature ratio assuming a thermal plasma, and the direction of the magnetic field dipole moment. We find that Atacama Large Millimeter/submillimeter Array (ALMA) observations of Sgr A* are inconsistent with all edge-on (i= 90°) models. Restricting attention to the MAD models favored by earlier EHT studies of Sgr A*, we find that only models with magnetic dipole moment pointing away from the observer are consistent with ALMA data. We also note that in 26 of the 27 passing MAD models, the accretion flow rotates clockwise on the sky. We provide a table of the means and standard deviations of the CP distributions for all model parameters, along with their trends. 

Models for the observational appearance of astrophysical black holes rely critically on accurate general-relativistic ray tracing and radiation transport to compute the intensity measured by a distant observer. In this paper, we illustrate how the choice of coordinates and initial conditions affect this process. In particular, we show that propagating rays from the camera to the source leads to different solutions if the spatial part of the momentum of the photon points towards the horizon or away from it. In doing this, we also show that coordinates that are well suited for numerical general-relativistic magnetohydrodynamic (GRMHD) simulations are typically not optimal for generic ray tracing. We discuss the implications for black hole images and show that radiation transport in optimal and nonoptimal spacetime coordinates lead to the same images up to numerical errors and algorithmic choices. 

Abstract Direct imaging of supermassive black holes (SMBHs) at event horizon-scale resolutions, as recently done by the Event Horizon Telescope, allows for testing alternative models to SMBHs such as Kerr naked singularities (KNSs). We demonstrate that the KNS shadow can be closed, open, or vanishing, depending on the spins and observational inclination angles. We study the critical parameters where the KNS shadow opens a gap, a distinctive phenomenon that does not happen with the black hole shadow. We show that the KNS shadow can only be closed for dimensionless spina≲ 1.18 and vanishing fora≳ 1.18 for certain ranges of inclination angles. We further analyze the effective angular momentum of photon orbits to demonstrate the fundamental connections between light geodesics and the KNS shadow geometry. We also perform numerical general relativistic ray-tracing calculations, which reproduce the analytical topological change in the KNS shadow, and illustrate other observational features within the shadow due to the lack of an event horizon. By comparing the geometric features of the KNS shadow with black hole shadow observations, the topological change in the shadow of KNSs can be used to test the cosmic censorship hypothesis and KNSs as alternative models to SMBHs. 

Abstract We present an in-depth analysis of the newly proposed correlation function in visibility space, between the E and B modes of linear polarization, hereafter the EB correlation, for a set of time-averaged general relativistic magnetohydrodynamical simulations compared with the phase map from different semianalytic models and the Event Horizon Telescope (EHT) 2017 data for M87*. We demonstrate that the phase map of time-averaged EB correlation contains novel information that might be linked to black hole (BH) spin, accretion state, and electron temperature. A detailed comparison with a semianalytic approach with different azimuthal expansion modes shows that to recover the morphology of real/imaginary part of the correlation function and its phase, we require higher orders of azimuthal modes. To extract the phase features, we use 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 geometrical ring models with various magnetic and velocity field morphologies, showing that both the image space and visibility-based EB -correlation morphologies in magnetically arrested disk  simulations can be explained with simple fluid and magnetic field geometries as used in ring models. Standard and normal evolutions 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 and electron temperatures, while others require larger theoretical surveys. 

Abstract The Event Horizon Telescope (EHT) has produced images of M87* and Sagittarius A*, and will soon produce time sequences of images, or movies. In anticipation of this, we describe a technique to measure the rotation rate, or pattern speed Ω_p, from movies using an autocorrelation technique. We validate the technique on Gaussian random field models with a known rotation rate and apply it to a library of synthetic images of Sgr A* based on general relativistic magnetohydrodynamics simulations. We predict that EHT movies will have Ω_p≈ 1° perGMc⁻³, which is of order 15% of the Keplerian orbital frequency in the emitting region. We can plausibly attribute the slow rotation seen in our models to the pattern speed of inward-propagating spiral shocks. We also find that Ω_pdepends strongly on inclination. Application of this technique will enable us to compare future EHT movies with the clockwise rotation of Sgr A* seen in near-infrared flares by GRAVITY. Pattern speed analysis of future EHT observations of M87* and Sgr A* may also provide novel constraints on black hole inclination and spin, as well as an independent measurement of black hole mass. 

Abstract The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which shows 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° 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 counterexamples, 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.