The imaging fidelity of the Event Horizon Telescope (EHT) is currently determined by its sparse baseline coverage. In particular, EHT coverage is dominated by long baselines, and is highly sensitive to atmospheric conditions and loss of sites between experiments. The limited short/mid-range baselines especially affect the imaging process, hindering the recovery of more extended features in the image. We present an algorithmic contingency for the absence of well-constrained short baselines in the imaging of compact sources, such as the supermassive black holes observed with the EHT. This technique enforces a specific second moment on the reconstructed image in the form of a size constraint, which corresponds to the curvature of the measured visibility function at zero baseline. The method enables the recovery of information lost in gaps of the baseline coverage on short baselines and enables corrections of any systematic amplitude offsets for the stations giving short-baseline measurements present in the observation. The regularization can use historical source size measurements to constrain the second moment of the reconstructed image to match the observed size. We additionally show that a characteristic size can be derived from available short-baseline measurements, extrapolated from other wavelengths, or estimated without complementary size constraints with parameter searches. We demonstrate the capabilities of this method for both static and movie reconstructions of variable sources.

It has been proposed that Very Long Baseline Interferometry (VLBI) at sub-millimeter 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 general relativity. This paper presents imaging simulations of a new Space VLBI mission concept. An initial design study of the concept has been presented as 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 with slightly different radii. Due to the relative drift of the satellites along the individual orbits, this setup will result in a dense spiral-shaped uv-coverage with long baselines (up to ~60 Glambda), allowing for extremely high-resolution and high-fidelity imaging of radio sources. We simulate EHI observations of general relativistic magnetohydrodynamics models of Sgr A* 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 if 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, closure phases could be used for imaging. Our simulations show that the EHI could be capable of imaging the black hole shadow of Sgr A* with a resolution of 4 uas (about 8% of the shadow diameter) within several months of observing time. The EHI concept could thus 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.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

The compact radio source Sagittarius~A$^*$ (Sgr~A$^*$)in the Galactic Center is the primary supermassive black hole candidate. General relativistic magnetohydrodynamical (GRMHD) simulations of the accretion flow around Sgr\,A$^*$ predict the presence of sub-structure at observing wavelengths of $\sim 3$\,mm and below (frequencies of 86\,GHz and above). For very long baseline interferometry (VLBI) observations of Sgr\,A$^*$ at this frequency the blurring effect of interstellar scattering becomes subdominant, and arrays such as the High Sensitivity Array (HSA) and the global mm-VLBI Array (GMVA) are now capable of resolving potential sub-structure in the source. Such investigations improve our understanding of the emission geometry of the mm-wave emission of Sgr\,A$^*$, which is crucial for constraining theoretical models and for providing a background to interpret 1\,mm VLBI data from the Event Horizon Telescope (EHT). We performed high-sensitivity very long baseline interferometry (VLBI) observations of Sgr\,A$^*$ at 3\,mm using the Very Long Baseline Array (VLBA) and the Large Millimeter Telescope (LMT) in Mexico on two consecutive days in May 2015, with the second epoch including the Green Bank Telescope (GBT). We find an overall source geometry that matches previous findings very closely, showing a deviation in fitted model parameters less than 3\% over a time scale of weeks and suggesting a highly stable global source geometry over time. The reported sub-structure in the 3\,mm emission of Sgr\,A$^*$ is consistent with theoretical expectations of refractive noise on long baselines. However, comparing our findings with recent results from 1\,mm and 7\,mm VLBI observations, which also show evidence for east-west asymmetry, an intrinsic origin cannot be excluded. Confirmation of persistent intrinsic substructure will require further VLBI observations spread out over multiple epochs.

We present a full 360 degree (i.e., 4$\pi$ steradian) general-relativistic ray-tracing and radiative transfer calculations of accreting supermassive black holes. We perform state-of-the-art three-dimensional general relativistic magnetohydrodynamical simulations using the BHAC code, subsequently post-processing this data with the radiative transfer code RAPTOR. All relativistic and general-relativistic effects, such as Doppler boosting and gravitational redshift, as well as geometrical effects due to the local gravitational field and the observer's changing position and state of motion, are therefore calculated self-consistently. Synthetic images at four astronomically-relevant observing frequencies are generated from the perspective of an observer with a full 360-degree view inside the accretion flow, who is advected with the flow as it evolves. As an example, we calculated images based on recent best-fit models of observations of Sagittarius A*. These images are combined to generate a complete 360-degree Virtual Reality movie of the surrounding environment of the black hole and its event horizon. Our approach also enables the calculation of the local luminosity received at a given fluid element in the accretion flow, providing important applications in, e.g., radiation feedback calculations onto black hole accretion flows. In addition to scientific applications, the 360-degree Virtual Reality movies we present also represent a new medium through which to communicate black hole physics to a wider audience, serving as a powerful educational tool.

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. 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. 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. 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%. 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.

We propose a new technique for radio interferometry to obtain super-resolution 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 $\ell _1$-norm and 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 (EHT). We confirm that $\ell _1$+TV regularization can achieve an optimal resolution of $\sim 25-30$\% of the diffraction limit $\lambda/D_{\rm max}$, 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 $\ell _1$ or TV regularizations alone. Furthermore, we find that $\ell _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 super-resolution 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.

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 $\ell_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 $\ell_1$+TV regularization can achieve an optimal resolution of $\sim 20-30$% of the diffraction limit $\lambda/D_{\rm max}$, 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 super-massive black hole in M87 and also the Galactic center Sgr A*.

We present the results of a closure phase analysis of 3 mm very long baseline interferometry (VLBI) measurements performed on Sagittarius A* (Sgr A*). We have analyzed observations made in May 2015 using the Very Long Baseline Array, the Robert C. Byrd Green Bank Telescope and the Large Millimeter Telescope Alfonso Serrano and obtained non-zero closure phase measurements on several station triangles - indicative of a non-point-symmetric source structure. The data are fitted with an asymmetric source structure model in Sgr A*, represented by a simple two-component model, which favours a fainter component due East of the main source. This result is discussed in light of a scattering screen with substructure or an intrinsically asymmetric source.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

We explore the conditions for the thermal instability to operate in the mini-spiral region of the Galactic centre (Sgr A*), where both the hot and cold media are known to coexist. The photoionisation Cloudy calculations are performed for different physical states of plasma. We neglect the dynamics of the material and concentrate on the study of the parameter ranges where the thermal instability may operate, taking into account the past history of Sgr A* bolometric luminosity. We show that the thermal instability does not operate at the present very low level of the Sgr A* activity. However, Sgr A* was much more luminous in the past. For the highest luminosity states the two-phase medium can be created up to 1.4 pc from the centre. The presence of dust grains tends to suppress the instability, but the dust is destroyed in the presence of strong radiation field and hot plasma. The clumpiness is thus induced in the high activity period, and the cooling/heating timescales are long enough to preserve later the past multi-phase structure. The instability enhances the clumpiness of the mini-spiral medium and creates a possibility of episodes of enhanced accretion of cold clumps towards Sgr A*. The mechanism determines the range of masses and sizes of clouds; under the conditions of Sgr A*, the likely values come out $1$ - $10^2M_{\oplus}$ for the cloud typical mass.