High resolution imaging of supermassive black holes shadows is a direct way to verify the theory of general relativity at extreme gravity conditions. Very Long Baseline Interferometry (VLBI) observations at millimeter/sub-millimeter wavelengths can provide such angular resolution for supermassive black holes, located in Sgr A* and M87. Recent VLBI observations of M87 with the Event Horizon Telescope (EHT) has shown such capabilities. The maximum obtainable spatial resolution of EHT is limited by Earth diameter and atmospheric phase variations. In order to improve the image resolution longer baselines are required. Radioastron space mission has successfully demonstrated the capabilities of Space-Earth VLBI with baselines much larger than Earth diameter. Millimetron is a next space mission of the Russian Space Agency that will operate at millimeter wavelengths. Nominal orbit of the observatory will be located around Lagrangian L2 point of the Sun-Earth system. In order to optimize the VLBI mode, we consider a possible second stage of the mission that could use near-Earth high elliptical orbit (HEO). In this contribution a set of near Earth orbits is used for the synthetic space-ground VLBI observations of Sgr A* and M87 in joint Millimetron and EHT configuration. General-relativistic magnetohydrodynamic models (GRMHD) for black hole environment of Sgr A* and M87 are used for static and dynamic imaging simulations at 230 GHz. A comparison is made between ground and space-ground baselines. It demonstrates that joined observations with Millimetron and EHT facilities significantly improve the image resolution and allow obtaining the dynamic snapshot images of Sgr A*.

Very Long Baseline Interferometry (VLBI) at sub-millimeter waves has the potential to image the shadow of the black hole in the Galactic Center, Sagittarius A* (Sgr A*), and thereby test basic predictions of the theory of general relativity. We investigate the imaging prospects of a new Space VLBI mission concept. The setup consists of two satellites in polar or equatorial circular Medium-Earth Orbits with slightly different radii, resulting in a dense spiral-shaped uv-coverage with long baselines, allowing for extremely high-resolution and high-fidelity imaging of radio sources. We simulate observations of a general relativistic magnetohydrodynamics model of Sgr A* for this configuration with noise calculated from model system parameters. After gridding the $uv$-plane and averaging visibilities accumulated over multiple months of integration, images of Sgr A* with a resolution of up to 4 $\mu$as could be reconstructed, allowing for stronger tests of general relativity and accretion models than with ground-based VLBI.

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 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 comparison with observational data. 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. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm 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. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results 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 the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, 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.

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

The galaxy M 87 harbors a kiloparsec-scale relativistic jet, whose origin coincides with a supermassive black hole. Observational mm-VLBI campaigns are capable of resolving the jet-launching region at the scale of the event horizon. In order to provide a context for interpreting these observations, realistic general-relativistic magnetohydrodynamical (GRMHD) models of the accretion flow are constructed. The characteristics of the observed spectral-energy distribution (SED) depend on the shape of the electrons' energy-distribution function (eDF). The dependency on the eDF is omitted in the modeling of the first Event Horizon Telescope results. In this work, we aim to model the M 87 SED from radio up to NIR/optical frequencies using a thermal-relativistic Maxwell- J\"uttner distribution, as well as a relativistic $\kappa$-distribution function. The electrons are injected based on sub-grid, particle-in-cell parametrizations for sub-relativistic reconnection. A GRMHD simulation in Cartesian-Kerr-Schild coordinates, using eight levels of adaptive mesh refinement (AMR), forms the basis of our model. To obtain spectra and images, the GRMHD data is post-processed with the ray-tracing code RAPTOR, which is capable of ray tracing through AMR GRMHD simulation data. We obtain radio spectra in both the thermal-jet and $\kappa$-jet models consistent with radio observations. Additionally, the $\kappa$-jet models also recover the NIR/optical emission. The models recover the observed source sizes and core shifts and obtain a jet power of $\approx 10^{43}$ ergs/s. In the $\kappa$-jet models, both the accretion rates and jet powers are approximately two times lower than the thermal-jet model. The frequency cut-off observed at $\nu \approx 10^{15}$ Hz is recovered when the accelerator size is $10^6$ - $10^8$ cm, this could potentially point to an upper limit for plasmoid sizes in the jet of M 87.

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.

Very long baseline interferometry (VLBI) makes it possible to recover images of astronomical sources with extremely high angular resolution. Most recently, the Event Horizon Telescope (EHT) has extended VLBI to short millimeter wavelengths with a goal of achieving angular resolution sufficient for imaging the event horizons of nearby supermassive black holes. VLBI provides measurements related to the underlying source image through a sparse set spatial frequencies. An image can then be recovered from these measurements by making assumptions about the underlying image. One of the most important assumptions made by conventional imaging methods is that over the course of a night's observation the image is static. However, for quickly evolving sources, such as the galactic center's supermassive black hole (Sgr A*) targeted by the EHT, this assumption is violated and these conventional imaging approaches fail. In this work we propose a new way to model VLBI measurements that allows us to recover both the appearance and dynamics of an evolving source by reconstructing a video rather than a static image. By modeling VLBI measurements using a Gaussian Markov Model, we are able to propagate information across observations in time to reconstruct a video, while simultaneously learning about the dynamics of the source's emission region. We demonstrate our proposed Expectation-Maximization (EM) algorithm, StarWarps, on realistic synthetic observations of black holes, and show how it substantially improves results compared to conventional imaging algorithms. Additionally, we demonstrate StarWarps on real VLBI data of the M87 Jet from the VLBA.

By linking widely separated radio dishes, the technique of very long baseline interferometry (VLBI) can greatly enhance angular resolution in radio astronomy. However, at any given moment, a VLBI array only sparsely samples the information necessary to form an image. Conventional imaging techniques partially overcome this limitation by making the assumption that the observed cosmic source structure does not evolve over the duration of an observation, which enables VLBI networks to accumulate information as the Earth rotates and changes the projected array geometry. Although this assumption is appropriate for nearly all VLBI, it is almost certainly violated for submillimeter observations of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), which has a gravitational timescale of only ~20 seconds and exhibits intra-hour variability. To address this challenge, we develop several techniques to reconstruct dynamical images ("movies") from interferometric data. Our techniques are applicable to both single-epoch and multi-epoch variability studies, and they are suitable for exploring many different physical processes including flaring regions, stable images with small time-dependent perturbations, steady accretion dynamics, or kinematics of relativistic jets. Moreover, dynamical imaging can be used to estimate time-averaged images from time-variable data, eliminating many spurious image artifacts that arise when using standard imaging methods. We demonstrate the effectiveness of our techniques using synthetic observations of simulated black hole systems and 7mm Very Long Baseline Array observations of M87, and we show that dynamical imaging is feasible for Event Horizon Telescope observations of Sgr A*.

General relativistic magnetohydrodynamic (GRMHD) simulations of accretion disks and jets associated with supermassive black holes show variability on a wide range of timescales. On timescales comparable to or longer than the gravitational timescale $t_G=GM/c^3$, variation may be dominated by orbital dynamics of the inhomogeneous accretion flow. Turbulent evolution within the accretion disk is expected on timescales comparable to the orbital period, typically an order of magnitude larger than $t_G$. For Sgr A*, $t_G$ is much shorter than the typical duration of a VLBI experiment, enabling us to study this variability within a single observation. Closure phases, the sum of interferometric visibility phases on a triangle of baselines, are particularly useful for studying this variability. In addition to a changing source structure, variations in observed closure phase can also be due to interstellar scattering, thermal noise, and the changing geometry of projected baselines over time due to Earth rotation. We present a metric that is able to distinguish the latter two from intrinsic or scattering variability. This metric is validated using synthetic observations of GRMHD simulations of Sgr A*. When applied to existing multi-epoch EHT data of Sgr A*, this metric shows that the data are most consistent with source models containing intrinsic variability from source dynamics, interstellar scattering, or a combination of those. The effects of black hole inclination, orientation, spin, and morphology (disk or jet) on the expected closure phase variability are also discussed.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

The black hole in the center of the Galaxy, associated with the compact source Sagittarius A* (Sgr A*), is predicted to cast a shadow upon the emission of the surrounding plasma flow, which encodes the influence of general relativity in the strong-field regime. The Event Horizon Telescope (EHT) is a Very Long Baseline Interferometry (VLBI) network with a goal of imaging nearby supermassive black holes (in particular Sgr A* and M87) with angular resolution sufficient to observe strong gravity effects near the event horizon. General relativistic magnetohydrodynamic (GRMHD) simulations show that radio emission from Sgr A* exhibits vari- ability on timescales of minutes, much shorter than the duration of a typical VLBI imaging experiment, which usually takes several hours. A changing source structure during the observations, however, violates one of the basic assumptions needed for aperture synthesis in radio interferometry imaging to work. By simulating realistic EHT observations of a model movie of Sgr A*, we demonstrate that an image of the average quiescent emission, featuring the characteristic black hole shadow and photon ring predicted by general relativity, can nonetheless be obtained by observing over multiple days and subsequent processing of the visibilities (scaling, averaging, and smoothing) before imaging. Moreover, it is shown that this procedure can be combined with an existing method to mitigate the effects of interstellar scattering. Taken together, these techniques allow the black hole shadow in the Galactic center to be recovered on the reconstructed image.