ALBERT

Description: Context. The observed spectral energy distribution of an accreting supermassive black hole typically forms a power-law spectrum in the near infrared (NIR) and optical wavelengths, that may be interpreted as a signature of accelerated electrons along the jet. However, the details of acceleration remain uncertain. Aim. In this paper, we study the radiative properties of jets produced in axisymmetric general relativistic magnetohydrodynamics (GRMHD) simulations of hot accretion flows onto underluminous supermassive black holes both numerically and semi-analytically, with the aim of investigating the differences between models with and without accelerated electrons inside the jet. Methods. We assume that electrons are accelerated in the jet regions of our GRMHD simulation. To model them, we modify the electrons’ distribution function in the jet regions from a purely relativistic thermal distribution to a combination of a relativistic thermal distribution and the κ-distribution function (the κ-distribution function is itself a combination of a relativistic thermal and a non-thermal power-law distribution, and thus it describes accelerated electrons). Inside the disk, we assume a thermal distribution for the electrons. In order to resolve the particle acceleration regions in the GRMHD simulations, we use a coordinate grid that is optimized for modeling jets. We calculate jet spectra and synchrotron maps by using the ray tracing code RAPTOR, and compare the synthetic observations to observations of Sgr A*. Finally, we compare numerical models of jets to semi-analytical ones. Results. We find that in the κ-jet models, the radio-emitting region size, radio flux, and spectral index in NIR/optical bands increase for decreasing values of the κ parameter, which corresponds to a larger amount of accelerated electrons. This is in agreement with analytical predictions. In our models, the size of the emission region depends roughly linearly on the observed wavelength λ, independently of the assumed distribution function. The model with κ = 3.5, ηacc = 5–10% (the percentage of electrons that are accelerated), and observing angle i = 30° fits the observed Sgr A* emission in the flaring state from the radio to the NIR/optical regimes, while κ = 3.5, ηacc 〈 1%, and observing angle i = 30° fit the upper limits in quiescence. At this point, our models (including the purely thermal ones) cannot reproduce the observed source sizes accurately, which is probably due to the assumption of axisymmetry in our GRMHD simulations. The κ-jet models naturally recover the observed nearly-flat radio spectrum of Sgr A* without invoking the somewhat artificial isothermal jet model that was suggested earlier. Conclusions. From our model fits we conclude that between 5% and 10% of the electrons inside the jet of Sgr A* are accelerated into a κ distribution function when Sgr A* is flaring. In quiescence, we match the NIR upper limits when this percentage is

Description: Context. Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are the most important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a direct comparison with observational data. Aims. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. Methods. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a millimetre VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects, compared to the addition of thermal noise only. Using synthetic data based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M 87, we performed case studies to assess the image quality that can be obtained with the current and future EHT array for different weather conditions. Results. Our synthetic observations show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of our GRMHD source models can be recovered robustly with the EHT2017 array after performing calibration steps, which include fringe fitting, a priori amplitude and network calibration, and self-calibration. With the planned addition of new stations to the EHT array in the coming years, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

Description: Context. Observational efforts to image the immediate environment of a black hole at the scale of the event horizon benefit from the development of efficient imaging codes that are capable of producing synthetic data, which may be compared with observational data. Aims. We aim to present RAPTOR, a new public code that produces accurate images, animations, and spectra of relativistic plasmas in strong gravity by numerically integrating the equations of motion of light rays and performing time-dependent radiative transfer calculations along the rays. The code is compatible with any analytical or numerical spacetime. It is hardware-agnostic and may be compiled and run both on GPUs and CPUs. Methods. We describe the algorithms used in RAPTOR and test the code’s performance. We have performed a detailed comparison of RAPTOR output with that of other radiative-transfer codes and demonstrate convergence of the results. We then applied RAPTOR to study accretion models of supermassive black holes, performing time-dependent radiative transfer through general relativistic magneto-hydrodynamical (GRMHD) simulations and investigating the expected observational differences between the so-called fast-light and slow-light paradigms. Results. Using RAPTOR to produce synthetic images and light curves of a GRMHD model of an accreting black hole, we find that the relative difference between fast-light and slow-light light curves is less than 5%. Using two distinct radiative-transfer codes to process the same data, we find integrated flux densities with a relative difference less than 0.01%. Conclusions. For two-dimensional GRMHD models, such as those examined in this paper, the fast-light approximation suffices as long as errors of a few percent are acceptable. The convergence of the results of two different codes demonstrates that they are, at a minimum, consistent.

Description: Context. Accreting supermassive black holes are sources of polarized radiation that propagates through highly curved spacetime before reaching the observer. Accurate and efficient numerical schemes for polarized radiative transfer in curved spacetime are needed to help interpret observations of such polarized emission. Aims. We aim to extend our publicly available radiative transfer code RAPTOR to include polarized radiative transfer, so that it can produce simulated polarized observations of accreting black holes. The RAPTOR code must remain compatible with arbitrary spacetimes and it must be efficient in operation, despite the added complexity of polarized radiative transfer. Methods. We provide a brief review of various codes and methods for covariant polarized radiative transfer available in the literature and existing codes, and we present an efficient new scheme. For the spacetime propagation aspect of the computation, we developed a compact, Lorentz-invariant representation of a polarized ray. For the plasma-propagation aspect of the computation, we performed a formal analysis of the stiffness of the polarized radiative-transfer equation with respect to our explicit integrator. We also developed a hybrid integration scheme that switches to an implicit integrator in case of stiffness in order to solve the equation with optimal speed and accuracy for all possible values of the local optical/Faraday thickness of the plasma. Results. We performed a comprehensive code verification by solving a number of well-known test problems using RAPTOR and comparing its output to exact solutions. We also demonstrate convergence with existing polarized radiative-transfer codes in the context of complex astrophysical problems, where we found that the integrated flux densities for all Stokes parameters converged to excellent agreement. Conclusions. The RAPTOR code is capable of performing polarized radiative transfer in arbitrary, highly curved spacetimes. This capability is crucial for interpreting polarized observations of accreting black holes, which can yield information about the magnetic-field configuration in such accretion flows. The efficient formalism implemented in RAPTOR is computationally light and conceptually simple. The code is publicly available.