Context. The supermassive black hole in the center of the Milky Way, Sgr A*, displays a nearly flat radio spectrum that is typical for jets in active galactic nuclei. Indeed, time-dependent ...magnetized models of radiatively inefficient accretion flows (RIAFs), which are commonly used to explain the millimeter, near-infrared, and X-ray emission of Sgr A*, often also produce jet-like outflows. However, the emission from these models has so far failed to reproduce the flat radio spectrum. Aims. We investigate whether current accretion simulations can produce the compact flat spectrum emission by simply using a different prescription for the heating of the radiating particles in the jet. Methods. We studied the radiative properties of accretion flows onto a black hole produced in time-dependent general-relativistic magnetohydrodynamic (GRMHD) simulations. A crucial free parameter in these models has always been the electron temperature, and here we allowed for variations in the proton-to-electron temperature ratios in the jet and disk. Results. We found that the flat spectrum is readily reproduced by a standard GRMHD model if one has an almost isothermal jet coupled to a two-temperature accretion flow. The low-frequency radio emission comes from the outflowing sheath of matter surrounding the strongly magnetized nearly empty jet. The model is consistent with the radio sizes and spectrum of Sgr A*. Conclusions. Hence, GRMHD models of accreting black holes can in principle naturally reproduce jets that match observed characteristics. For Sgr A* the model fit to the spectrum predicts higher mass-accretion rates when a jet is included than without a jet. Hence, the impact of the recently discovered G2 cloud that is expected to be accreted onto Sgr A* might be less severe than currently thought.
Abstract Synchrotron emissivities, absorptivities, and Faraday rotation and conversion coefficients are needed in modeling a variety of astrophysical sources, including Event Horizon Telescope (EHT) ...sources. We develop a method for estimating transfer coefficients that exploits their linear dependence on the electron distribution function, decomposing the distribution function into a sum of parts each of whose emissivity can be calculated easily. We refer to this procedure as stochastic averaging and apply it in two contexts. First, we use it to estimate the emissivity of an isotropic κ distribution function with a high-energy cutoff. The resulting coefficients can be evaluated efficiently enough to be used directly in ray-tracing calculations, and we provide an example calculation. Second, we use stochastic averaging to assess the effect of subgrid turbulence on the volume-averaged emissivity and along the way provide a prescription for a turbulent emissivity. We find that for parameters appropriate to EHT sources turbulence reduces the emissivity slightly. In the infrared, turbulence can dramatically increase the emissivity.
Abstract
The Event Horizon Telescope (EHT) has recently released high-resolution images of accretion flows onto two supermassive black holes. Our physical understanding of these images depends on the ...accuracy and precision of numerical models of plasma and radiation around compact objects. The goal of this work is to speed up radiative-transfer simulations used to create mock images of black holes for comparison with the EHT observations. A ray-tracing code for general relativistic and fully polarized radiative transfer through plasma in strong gravity is ported onto a graphics processing unit (GPU). We describe our GPU implementation and carry out speedup tests using models of optically thin advection-dominated accretion flow onto a black hole realized semianalytically and in 3D general relativistic magnetohydrodynamic simulations, low and very high image pixel resolutions, and two different sets of CPU+GPUs. We show that a GPU with high double precision computing capability can significantly reduce the image production computational time, with a speedup factor of up to approximately 1200. The significant speedup facilitates, e.g., dynamic model fitting to the EHT data, including polarimetric data. The method extension may enable studies of emission from plasma with nonthermal particle distribution functions for which accurate approximate synchrotron emissivities are not available. The significant speedup reduces the carbon footprint of the generation of the EHT image libraries by at least an order of magnitude.
Context. Radiatively inefficient accretion flows (RIAFs) are believed to power supermassive black holes in the underluminous cores of galaxies. Such black holes are typically accompanied by ...flat-spectrum radio cores indicating the presence of moderately relativistic jets. One of the best constrained RIAFs is associated with the supermassive black hole in the Galactic center, Sgr A*. Aims. Since the plasma in RIAFs is only weakly collisional, the dynamics and the radiative properties of these systems are very uncertain. Here we want to study the impact of varying electron temperature on the appearance of accretion flows and jets. Methods. Using three-dimensional general relativistic magnetohydrodynamics accretion flow simulations, we use ray tracing methods to predict spectra and radio images of RIAFs allowing for different electron heating mechanisms in the in- and outflowing parts of the simulations. Results. We find that small changes in the electron temperature can result in dramatic differences in the relative dominance of jets and accretion flows. Application to Sgr A* shows that radio spectrum and size of this source can be well reproduced with a model where electrons are more efficiently heated in the jet. The X-ray emission is sensitive to the electron heating mechanism in the jets and disk and therefore X-ray observations put strong constraints on electron temperatures and geometry of the accretion flow and jet. For Sgr A*, the jet model also predicts a significant frequency-dependent core shift which could place independent constraints on the model once measured accurately. Conclusions. We conclude that more sophisticated models for electron distribution functions are crucial for constraining GRMHD simulations with actual observations. For Sgr A*, the radio appearance may well be dominated by the outflowing plasma. Nonetheless, at the highest radio frequencies, the shadow of the event horizon should still be detectable with future Very Long Baseline Interferometric observations.
Context. It has been proposed that Very Long Baseline Interferometry (VLBI) at submillimeter waves will allow us to image the shadow of the black hole in the center of our Milky Way, Sagittarius A* ...(Sgr A*), and thereby test basic predictions of the theory of general relativity. Aims. This paper presents imaging simulations of a new Space VLBI (SVLBI) mission concept. An initial design study of the concept has been presented in the form of the Event Horizon Imager (EHI). The EHI may be suitable for imaging Sgr A* at high frequencies (up to ∼690 GHz), which has significant advantages over performing ground-based VLBI at 230 GHz. The concept EHI design consists of two or three satellites in polar or equatorial circular medium-Earth orbits (MEOs) with slightly different radii. Due to the relative drift of the satellites along the individual orbits over the course of several weeks, this setup will result in a dense spiral-shaped uv-coverage with long baselines (up to ∼60 Gλ), allowing for extremely high-resolution and high-fidelity imaging of radio sources. Methods. We simulated observations of general relativistic magnetohydrodynamics (GRMHD) models of Sgr A* for the proposed configuration and calculate the expected noise based on preliminary system parameters. On long baselines, where the signal-to-noise ratio may be low, fringes could be detected assuming that the system is sufficiently phase stable and the satellite orbits can be reconstructed with sufficient accuracy. Averaging visibilities accumulated over multiple epochs of observations could then help improving the image quality. With three satellites instead of two, closure phases could be used for imaging. Results. Our simulations show that the EHI could be capable of imaging the black hole shadow of Sgr A* with a resolution of 4 μas (about 8% of the shadow diameter) within several months of observing time. Conclusion. Our preliminary study of the EHI concept shows that it is potentially of high scientific value. It could be used to measure black hole shadows much more precisely than with ground-based VLBI, allowing for stronger tests of general relativity and accretion models.
ABSTRACT
Accreting black holes tend to display a characteristic dark central region called the black hole shadow, which depends only on space–time/observer geometry and which conveys information ...about the black hole’s mass and spin. Conversely, the observed central brightness depression, or image shadow, additionally depends on the morphology of the emission region. In this paper, we investigate the astrophysical requirements for observing a meaningful black hole shadow in GRMHD-based models of accreting black holes. In particular, we identify two processes by which the image shadow can differ from the black hole shadow: evacuation of the innermost region of the accretion flow, which can render the image shadow larger than the black hole shadow, and obscuration of the black hole shadow by optically thick regions of the accretion flow, which can render the image shadow smaller than the black hole shadow, or eliminate it altogether. We investigate in which models the image shadows of our models match their corresponding black hole shadows, and in which models the two deviate from each other. We find that, given a compact and optically thin emission region, our models allow for measurement of the black hole shadow size to an accuracy of 5 per cent. We show that these conditions are generally met for all MAD simulations we considered, as well as some of the SANE simulations.
We propose a new imaging technique for radio and optical/infrared interferometry. The proposed technique reconstructs the image from the visibility amplitude and closure phase, which are standard ...data products of short-millimeter very long baseline interferometers such as the Event Horizon Telescope (EHT) and optical/infrared interferometers, by utilizing two regularization functions: the 1-norm and total variation (TV) of the brightness distribution. In the proposed method, optimal regularization parameters, which represent the sparseness and effective spatial resolution of the image, are derived from data themselves using cross-validation (CV). As an application of this technique, we present simulated observations of M87 with the EHT based on four physically motivated models. We confirm that 1 + TV regularization can achieve an optimal resolution of ∼20%-30% of the diffraction limit λ/Dmax, which is the nominal spatial resolution of a radio interferometer. With the proposed technique, the EHT can robustly and reasonably achieve super-resolution sufficient to clearly resolve the black hole shadow. These results make it promising for the EHT to provide an unprecedented view of the event-horizon-scale structure in the vicinity of the supermassive black hole in M87 and also the Galactic center Sgr A*.
We propose a new technique for radio interferometry to obtain superresolution full-polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the ...image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the 1 norm and the total variation (TV) of the brightness distribution. As an application of this technique, we present simulated linear polarization observations of two physically motivated models of M87 with the Event Horizon Telescope. We confirm that 1+TV regularization can achieve an optimal resolution of ∼25%-30% of the diffraction limit , which is the nominal spatial resolution of a radio interferometer for both the total intensity (i.e., Stokes I) and linear polarizations (i.e., Stokes Q and U). This optimal resolution is better than that obtained from the widely used Cotton-Schwab CLEAN algorithm or from using 1 or TV regularizations alone. Furthermore, we find that 1+TV regularization can achieve much better image fidelity in linear polarization than other techniques over a wide range of spatial scales, not only in the superresolution regime, but also on scales larger than the diffraction limit. Our results clearly demonstrate that sparse reconstruction is a useful choice for high-fidelity full-polarimetric interferometric imaging.
Millimeter polarimetry of Sgr A* probes the linearly polarized emission region on a scale of ∼10 Schwarzschild radii (RS), as well as the dense, magnetized accretion flow on scales out to the Bondi ...radius (∼105RS) through Faraday rotation. We present here multi-epoch ALMA Band 6 (230 GHz) polarimetry of Sgr A*. The results confirm a mean rotation measure, , consistent with measurements over the past 20 yr, and support an interpretation of the RM as originating from a radiatively inefficient accretion flow with . Variability is observed for the first time in the RM on timescales that range from hours to months. The long-term variations may be the result of changes in the line-of-sight properties in a turbulent accretion flow. Short-term variations in the apparent RM are not necessarily the result of Faraday rotation and may be the result of complex emission and propagatation effects close to the black hole, some of which have been predicted in numerical modeling. We also confirm the detection of circular polarization at a mean value of −1.1% 0.2%. It is variable in amplitude on timescales from hours to months, but the handedness remains unchanged from that observed in past centimeter- and millimeter-wavelength detections. These results provide critical constraints for the analysis and interpretation of Event Horizon Telescope data of Sgr A*, M87, and similar sources.
The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the ...analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe the PATOKA pipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use, iharm, ipole, and igrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced by PATOKA and conclude with a discussion of limitations and future directions.