Recent observations of SgrA* by the GRAVITY instrument have astrometrically tracked infrared flares (IR) at distances of $\sim 10$ gravitational radii ($r_g$). In this paper, we study a model for the flares based on 3D general relativistic magnetohydrodynamic (GRMHD) simulations of magnetically arrested accretion disks (MADs) which exhibit violent episodes of flux escape from the black hole magnetosphere. These events are attractive for flare modeling for several reasons: i) the magnetically dominant regions can resist being disrupted via magneto-rotational turbulence and shear, ii) the orientation of the magnetic field is predominantly vertical as suggested by the GRAVITY data, iii) magnetic reconnection associated with the flux eruptions could yield a self-consistent means of particle heating/acceleration during the flare events. In this analysis we track erupted flux bundles and provide distributions of sizes, energies and plasma parameter. In our simulations, the orbits tend to circularize at a range of radii from $\sim 5-40 r_g$. The magnetic energy contained within the flux bundles ranges up to $\sim10^{40}$ erg, enough to power IR and X-ray flares. We find that the motion within the magnetically supported flow is substantially sub-Keplerian, in tension with the inferred period-radius relation of the three GRAVITY flares.

The Event Horizon Telescope recently observed the first shadow of a black hole. Images like this can potentially be used to test or constrain theories of gravity and deepen the understanding in plasma physics at event horizon scales, which requires accurate parameter estimations. In this work, we present Deep Horizon, two convolutional deep neural networks that recover the physical parameters from images of black hole shadows. We investigate the effects of a limited telescope resolution and observations at higher frequencies. We trained two convolutional deep neural networks on a large image library of simulated mock data. The first network is a Bayesian deep neural regression network and is used to recover the viewing angle $i$, and position angle, mass accretion rate $\dot{M}$, electron heating prescription $R_{\rm high}$ and the black hole mass $M_{\rm BH}$. The second network is a classification network that recovers the black hole spin $a$. We find that with the current resolution of the Event Horizon Telescope, it is only possible to accurately recover a limited number of parameters of a static image, namely the mass and mass accretion rate. Since potential future space-based observing missions will operate at frequencies above 230 GHz, we also investigated the applicability of our network at a frequency of 690 GHz. The expected resolution of space-based missions is higher than the current resolution of the Event Horizon Telescope, and we show that Deep Horizon can accurately recover the parameters of simulated observations with a comparable resolution to such missions.

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.

The observations of GW170817/GRB170817A have confirmed that the coalescence of a neutron-star binary is the progenitor of a short gamma-ray burst. In the standard picture of a short gamma-ray burst, a collimated highly relativistic outflow is launched after merger and it successfully breaks out from the surrounding ejected matter. Using initial conditions inspired from numerical-relativity binary neutron-star merger simulations, we have performed general-relativistic hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations in which the jet is launched and propagates self-consistently. The complete set of simulations suggests that: (i) MHD jets have an intrinsic energy and velocity polar structure with a ``hollow core'' subtending an angle $\theta_{\rm core}\approx4^{\circ}-5^{\circ}$ and an opening angle of $\theta_{\rm jet}\gtrsim10^{\circ}$; (ii) MHD jets eject significant amounts of matter and two orders of magnitude more than HD jets; (iii) the energy stratification in MHD jets naturally yields the power-law energy scaling $E(>\Gamma\beta)\propto(\Gamma\beta)^{-4.5}$; (iv) MHD jets provide fits to the afterglow data from GRB170817A that are comparatively better than those of the HD jets and without free parameters; (v) finally, both of the best-fit HD/MHD models suggest an observation angle $\theta_{\rm obs} \simeq 21^{\circ}$ for GRB170817A.

One of the main dissipation processes acting on all scales in relativistic jets is thought to be governed by magnetic reconnection. Such dissipation processes have been studied in idealized environments, such as reconnection layers, which evolve in merging islands and lead to the production of plasmoids, ultimately resulting in efficient particle acceleration. In accretion flows onto black holes, reconnection layers can be developed and destroyed rapidly during the turbulent evolution of the flow. We present a series of two-dimensional general-relativistic magnetohydrodynamic simulations of tori accreting onto rotating black holes focusing our attention on the formation and evolution of current sheets. Initially, the tori are endowed with a poloidal magnetic field having a multi-loop structure along the radial direction and with an alternating polarity. During reconnection processes, plasmoids and plasmoid chains are developed leading to a flaring activity and hence to a variable electromagnetic luminosity. We describe the methods developed to track automatically the plasmoids that are generated and ejected during the simulation, contrasting the behaviour of multi-loop initial data with that encountered in typical simulations of accreting black holes having initial dipolar field composed of one loop only. Finally, we discuss the implications that our results have on the variability to be expected in accreting supermassive black holes.

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.

Upcoming VLBI observations will resolve nearby supermassive black holes, most notably Sagittarius A* and M87, on event horizon-scales. Recent observations of Sagittarius A* with the Event Horizon Telescope have revealed horizon-scale structure. Accordingly, the detection and measurement of the back hole ``shadow" is expected to enable the existence of astrophysical black holes to be verified directly. Although the theoretical description of the shadow is straightforward, its observational appearance is largely determined by the properties of the surrounding accretion flow, which is highly turbulent. We introduce a new polarised general-relativistic radiative transfer code, \texttt{BHOSS}, which accurately solves the equations of polarised radiative transfer in arbitrary strong-gravity environments, providing physically-realistic images of astrophysical black holes on event horizon-scales, as well as also providing insight into the fundamental properties and nature of the surrounding accretion flow environment.

Worldwide very long baseline radio interferometry arrays are expected to obtain horizon-scale images of supermassive black hole candidates as well as of relativistic jets in several nearby active galactic nuclei. This motivates the development of models for magnetohydrodynamic flows in strong gravitational fields. The Black Hole Accretion Code (BHAC) intends to aid with the modelling of such sources by means of general relativistic magnetohydrodynamical (GRMHD) simulations in arbitrary stationary spacetimes. New additions were required to guarantee an accurate evolution of the magnetic field when small and large scales are captured simultaneously. We discuss the adaptive mesh refinement (AMR) techniques employed in BHAC, essential to keep several problems computationally tractable, as well as staggered-mesh-based constrained transport (CT) algorithms to preserve the divergence-free constraint of the magnetic field, including a general class of prolongation operators for face-allocated variables compatible with them. Through several standard tests, we show that the choice of divergence-control method can produce qualitative differences in simulations of scientifically relevant accretion problems. We demonstrate the ability of AMR to reduce the computational costs of accretion simulations while sufficiently resolving turbulence from the magnetorotational instability. In particular, we describe a simulation of an accreting Kerr black hole in Cartesian coordinates using AMR to follow the propagation of a relativistic jet while self-consistently including the jet engine, a problem set up-for which the new AMR implementation is particularly advantageous. The CT methods and AMR strategies discussed here are being employed in the simulations performed with BHAC used in the generation of theoretical models for the Event Horizon Telescope Collaboration.

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.

Radio-astronomical observations of the supermassive black-hole candidate in the galactic center will soon offer the possibility to study gravity in its strongest regimes and to test different models for these compact objects. Studies based on semi-analytic models and strong-field images of stationary plasma configurations around boson stars have stressed the difficulty to distinguish them from black holes. For the first time, we present general-relativistic magnetohydrodynamic simulations followed by consistent general-relativistic radiative-transfer calculations in a boson-star spacetime to revisit the appearance of these objects in a realistic accretion scenario. We find that the absence of an event horizon in a boson star leads to important differences in the dynamics of the accretion and results in both the formation of a small torus in the interior of the boson star and in the absence of an evacuated high-magnetization funnel in the polar regions. Synthetic reconstructed images considering realistic astronomical observing conditions show that differences in the appearance of the two compact objects are large enough to be detectable. These results, which also apply to other horizonless compact objects, strengthen confidence in the ability to determine the presence of an event horizon via radio observations and highlight the importance of self-consistent multidimensional simulations to study the compact object at the galactic center.

High-resolution Very-Long-Baseline Interferometry observations of NGC 1052 show a two sided jet with several regions of enhanced emission and a clear emission gap between the two jets.This gap shrinks with increasing frequency and vanishes around $\nu\sim43$ GHz. The observed structures are due to both the macroscopic fluid dynamics interacting with the surrounding ambient medium including an obscuring torus and the radiation microphysics. In this paper we investigate the possible physical conditions in relativistic jets of NGC 1052 by directly modelling the observed emission and spectra via state-of-the-art special-relativistic hydrodynamic (SRHD) simulations and radiative transfer calculations. To investigate the physical conditions in the relativistic jet we coupled our radiative transfer code to evolutionary algorithms and performed simultaneous modelling of the observed jet structure and the broadband radio spectrum. During the calculation of the radiation we consider both thermal and non-thermal emission. In order to compare our model to VLBI observations we take into account the sparse sampling of the u-v plane, the array properties and the imaging algorithm. We present for the first time an end-to-end pipeline for fitting numerical simulations to VLBI observations of relativistic jets taking into account the macrophysics including fluid dynamics and ambient medium configurations together with thermal/non-thermal emission and the properties of the observing array. The detailed analysis of our simulations shows that the structure and properties of the observed relativistic jets in NGC 1052 can be reconstructed by a slightly over-pressured jet ($d_k\sim1.5$) embedded in a decreasing pressure ambient medium

Context: Few Active Galactic Nuclei (AGN) reveal a double sided jet system. However, these are crucial to understand basic physical properties of extragalactic jets. Aims: We address the questions whether jets in AGN are symmetric in nature, how well they are collimated on small scales, and how do they evolve with time. Methods: We monitored the sub-parsec scale morphology of NGC1052 with the VLBA at 43GHz from 2005 to 2009. Results: A detailed study of 29 epochs show an asymmetry between both jets. A kinematic analysis of the outflows reveals higher apparent velocities for the eastern jet as compared to the western one: beta_ej = 0.529+-0.038 and beta_wj = 0.343+-0.037, respectively. Contradictory to previous studies, we find higher flux densities for the western jet. The distribution of brightness temperature and jet width features well collimated jets up to 1mas distance to the dynamic center and a nearly conical outflow further outwards. Combining flux density ratios and velocities of the jet flow, we were not able to find a combination of intrinsic velocity and inclination angle of the jets that is consistent for all observations, which is contradictory to symmetrically evolving jets. Spectral index maps between quasi-simultaneous 22GHz and 43GHz observations support the existence of an optically thick absorber covering the innermost 1.6 mas around the 43GHz central feature and an optically thin jet emission with a spectral index of <= -1. Conclusions: Our results fit into a picture in which we expect larger internal energy and/or magnetic flux in the western jet and higher kinetic energy in the eastern jet. Previous observations at lower frequencies found slower velocities of the moving jet features as compared to this work. Considering the different velocities at different areas we suggest a spine-sheath structure with a faster inner layer and a slower outer layer.

Our Galactic Center, Sagittarius A* (Sgr A*), is believed to harbour a supermassive black hole (BH), as suggested by observations tracking individual orbiting stars. Upcoming sub-millimetre very-long-baseline-interferometry (VLBI) images of Sgr A* carried out by the Event-Horizon-Telescope Collaboration (EHTC) are expected to provide critical evidence for the existence of this supermassive BH. We assess our present ability to use EHTC images to determine if they correspond to a Kerr BH as predicted by Einstein's theory of general relativity (GR) or to a BH in alternative theories of gravity. To this end, we perform general-relativistic magnetohydrodynamical (GRMHD) simulations and use general-relativistic radiative transfer (GRRT) calculations to generate synthetic shadow images of a magnetised accretion flow onto a Kerr BH. In addition, and for the first time, we perform GRMHD simulations and GRRT calculations for a dilaton BH, which we take as a representative solution of an alternative theory of gravity. Adopting the VLBI configuration from the 2017 EHTC campaign, we find that it could be extremely difficult to distinguish between BHs from different theories of gravity, thus highlighting that great caution is needed when interpreting BH images as tests of GR.

High-resolution Very-Long-Baseline Interferometry observations of active galactic nuclei have revealed asymmetric structures in the jets of radio galaxies. These asymmetric structures may be due to internal asymmetries in the jet, could be induced by the different conditions in the surrounding ambient medium including the obscuring torus, or a combination of the two. In this paper we investigate the influence of the ambient medium (including the obscuring torus) on the observed properties of jets from radio galaxies. We performed special-relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets using the special-relativistic hydrodynamics code \texttt{Ratpenat}, which is based on a second-order accurate finite-volume method and an approximate Riemann solver. Using a newly developed emission code to compute the electromagnetic emission, we have investigated the influence of different ambient medium and torus configurations on the jet structure and subsequently computed the non-thermal emission produced by the jet and the thermal absorption due to the torus. To better compare the emission simulations with observations we produced synthetic radio maps, taking into account the properties of the observatory. The detailed analysis of our simulations shows that the observed asymmetries can be produced by the interaction of the jet with the ambient medium and by the absorption properties of the obscuring torus.

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

High resolution Very Long Baseline Interferometry (VLBI) observations of Active Galactic Nuclei (AGN) revealed traveling and stationary or quasi-stationary radio-components in several blazar jets. The traveling ones are in general interpreted as shock waves generated by pressure perturbations injected at the jet nozzle. The stationary features can be interpreted as recollimation shocks in non-pressure matched jets if they show a quasi-symmetric bump in the spectral index distribution. In some jets there may be interactions between the two kinds of shocks. These shock--shock interactions are observable with VLBI techniques, and their signature should also be imprinted on the single--dish light curves. We performed relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets. To simulate the shock interaction we injected a perturbation at the jet nozzle once a steady-state was reached. We computed the non-thermal emission (including adiabatic and synchotron losses) resulting from the simulation. We show that the injection of perturbations in a jet can produce a bump in emission at GHz frequencies previous to the main flare, which is produced when the perturbation fills the jet in the observer's frame. The detailed analysis of our simulations and the non-thermal emission calculations show that interaction between a recollimation shock and traveling shock produce a typical and clear signature in both the single--dish light curves and in the VLBI observations: the flaring peaks are higher and delayed with respect to the evolution of a perturbation through a conical jet. This fact can allow to detect such interactions for stationary components lying outside of the region in where the losses are dominated by inverse Compton scattering.

Blazars are among the most powerful extragalactic objects, as a sub-class of active galactic nuclei. They launch relativistic jets and their emitted radiation shows strong variability across the entire electro-magnetic spectrum. The mechanisms producing the variability are still controversial and different models have been proposed to explain the observed variations in multi-frequency blazar light curves.We investigate the capabilities of the classical shock-in-jet model to explain and reconstruct the observed evolution of flares in the turnover frequency turnover flux density plane and their frequency-dependent light curve parameters. With a detailed parameter space study we provide the framework for future, detailed comparisons of observed flare signatures with the shock-in-jet scenario. Based on the shock model we compute synthetic single-dish light curves at different radio frequencies (2.6 to 345 GHz) and for different physical conditions in a conical jet (e.g. magnetic field geometry and Doppler factor). From those we extract the slopes of the different energy loss stages within the $\nu_\mathrm{m}$-$S_\mathrm{m}$ plane and deduce the frequency-dependence of different light curve parameters such as flare amplitude, time scale and cross-band delays. The evolution of the Doppler factor along the jet has the largest influence on the evolution of the flare and on the frequency-dependent light curve parameters. The synchrotron stage can be hidden in the Compton or in the adiabatic stage, depending mainly on the evolution of the Doppler factor, which makes it difficult to detect its signature in observations. In addition, we show that the time lags between different frequencies can be used as an efficient tool to better constrain the physical properties of these objects.

Relativistic jets in active galactic nuclei represent one of the most powerful phenomena in the Universe. They form in the surroundings of the supermassive black holes as a by-product of accretion onto the central black hole in active galaxies. The flow in the jets propagates at velocities close to the speed of light. The distance between the first part of the jet that is visible in radio images (core) and the black hole is still a matter of debate. Only very-long-baseline interferometry observations resolve the innermost compact regions of the radio jet. Those can access the jet base, and combining data at different wavelenghts, address the physical parameters of the outflow from its emission. We have performed an accurate analysis of the frequency-dependent shift of the VLBI core location for a multi-wavelength set of images of the blazar CTA 102 including data from 6 cm down to 3 mm. The measure of the position of the central black hole, with mass $\sim 10^{8.93}\,M_\odot$, in the blazar CTA 102 reveals a distance of $\sim 8\times10^4$ gravitational radii to the 86 GHz core, in agreement with similar measures obtained for other blazars and distant radio galaxies, and in contrast with recent results for the case of nearby radio galaxies, which show distances between the black hole and the radio core that can be two orders of magnitude smaller.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

To obtain a better understanding of the location and mechanisms for the production of the gamma-ray emission in jets of AGN we present a detailed study of the HST-1 structure, 0.8 arcsec downstream the jet of M87, previously identified as a possible candidate for TeV emission. HST-1 shows a very peculiar structure with superluminal as well as possible stationary sub-components, and appears to be located in the transition from a parabolic to a conical jet shape, presumably leading to the formation of a recollimation shock. This scenario is supported by our new RHD simulations in which the interaction of a moving component with a recollimation shock leads to the appearance of a new superluminal component. To discern whether HST-1 is produced by a recollimation shock or some other MHD instability, we present new polarimetric 2.2 and 5 GHz VLBA, as well as 15, 22 and 43 GHz JVLA observations obtained between November 2012 and March 2013.

If we assume that nuclei of extragalactic radio sources contain binary black hole systems, the two black holes can eject VLBI components in which case two families of different VLBI trajectories will be observed. Another important consequence of a binary black hole system is that the VLBI core is associated with one black hole, and if a VLBI component is ejected by the second black hole, one expects to be able to detect the offset of the origin of the VLBI component ejected by the black hole that is not associated with the VLBI core. The ejection of VLBI components is perturbed by the precession of the accretion disk and the motion of the black holes around the center of gravity of the binary black hole system. We modeled the ejection of the component taking into account the two pertubations and present a method to fit the coordinates of a VLBI component and to deduce the characteristics of the binary black hole system. Specifically, this is the ratio Tp/Tb where Tp is the precession period of the accretion disk and Tb is the orbital period of the binary black hole system, the mass ratio M1/M2, and the radius of the Binary Black Hole system Rbin. From the variations of the coordinates as a function of time of the ejected VLBI component, we estimated the inclination angle io and the bulk Lorentz factor gamma of the modeled component. We applied the method to component S1 of 1823+568 and to component C5 of 3C 279.

The temporal and spatial spectral evolution of the jets of AGN can be studied with multi-frequency, multi-epoch VLBI observations. The combination of both, morphological and spectral parameters can be used to derive source intrinsic physical properties such as the magnetic field and the non-thermal particle density. In the first two papers of this series, we analyzed the single-dish light curves and the VLBI kinematics of the blazar CTA 102 and suggested a shock-shock interaction between a traveling and a standing shock wave as a possible scenario to explain the observed evolution of the component associated to the 2006 flare. In this paper we investigate the core-shift and spectral evolution to test our hypothesis of a shock-shock interaction. We used 8 multi-frequency VLBA observations to analyze the temporal and spatial evolution of the spectral parameters during the flare. We observed CTA 102 between May 2005 and April 2007 using the VLBA at six different frequencies spanning from 2 - 86 GHz. After the calibrated VLBA images were corrected for opacity, we performed a detailed spectral analysis. From the derived values we estimated the magnetic field and the density of the relativistic particles. The detailed analysis of the opacity shift reveals that the position of the jet core is proportional to nu^-1 with some temporal variations. The value suggests possible equipartition between magnetic field energy and particle kinetic energy densities at the most compact regions. From the variation of the physical parameters we deduced that the 2006 flare in CTA 102 is connected to the ejection of a new traveling feature (t=2005.9) and the interaction between this shock wave and a stationary structure around 0.1 mas from the core. The source kinematics together with the spectral and structural variations can be described by helical motions in an over-pressured jet.

The 2006 radio flare in the blazar CTA102 is the largest ever reported for this source. The analysis of the single--dish light curves revealed a possible shock-shock interaction as a scenario for the 2006 outburst. In order to confirm this hypothesis we analyzed eight multi-frequency (2 GHz - 86 GHz) VLBI observations covering the 2006 flare. We detected a radio component ejected around 2005.9 and a stationary feature around 0.1 mas away from the core. We also performed a core-shift and spectral analysis on the fully calibrated VLBI maps. The frequency dependent position of the core slightly deviates from the typical $\nu^{-1}$-dependence during the flaring state. This behaviour may reflect a mismatch between magnetic energy density and energy density of the relativistic particles. From the core-shift corrected maps we obtained the spectral parameters and computed the temporal and spatial evolution of the jet intrinsic parameters such as the magnetic field and the particle density. The overall picture of the source could be best interpreted as an over-pressured jet perturbed by traveling shock waves.

We performed multifrequency multiepoch Very Long Baseline Array (VLBA) observations of the blazar CTA 102 during its 2006 radio flare, the strongest ever reported for this source. These observations provide an excellent opportunity to investigate the evolution of the physical properties of blazars, especially during these flaring events. We want to study the kinematic changes in the source during the strong radio outburst in April 2006 and test the assumption of a shock-shock interaction. This assumption is based on the analysis and modeling of the single-dish observations of CTA\,102 (Paper I). In this paper we study the kinematics of CTA 102 at several frequencies using VLBI observations. From the modeled jet features we derived estimates for the evolution of the physical parameters, such as the particle density and the magnetic field. Furthermore, we combined our observations during the 2006 flare with long-term VLBA monitoring of the source at 15 GHz and 43 GHz. We cross-identified seven features throughout our entire multifrequency observations and find evidence of two possible recollimation shocks around 0.1 mas (deprojected 18 pc at a viewing angle of 2.6 degrees) and 6.0 mas (deprojected 1 kpc) from the core. The 43 GHz observations reveal a feature ejected at epoch $t_\mathrm{ej}=2005.9\pm0.2$, which could be connected to the 2006 April radio flare. Furthermore, this feature might be associated with the traveling component involved in the possible shock-shock interaction, which gives rise to the observed double peak structure in the single-dish light curves reported in Paper I.

The F-GAMMA program is among the most comprehensive programs that aim at understanding the physics in active galactic nuclei through the multi-frequency monitoring of Fermi blazars. Here we discuss monthly sampled broad-band radio spectra (2.6 - 142 GHz). Two different studies are presented. (a) We discuss that the variability patterns traced can be classified into two classes: (1) to those showing intense spectral-evolution and (2) those showing a self-similar quasi-achromatic behaviour. We show that a simple two-component model can very well reproduce the observed phenomenologies. (b) We present the cm-to-mm behaviour of three gamma-ray bright Narrow Line Seyfert 1 galaxies over time spans varying between ~1.5 and 3 years and compare their variability characteristics with typical blazars.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of Fermi blazars. In the current study we show and discuss the evolution of broad-band radio spectra, which are measured at ten frequencies between 2.64 and 142 GHz using the Effelsberg 100-m and the IRAM 30-m telescopes. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. It is argued that these can be attributed to only two classes of variability mechanisms. The first four types are dominated by spectral evolution and can be described by a simple two-component system composed of: (a) a steep quiescent spectral component from a large scale jet and (b) a time evolving flare component following the "Shock-in-Jet" evolutionary path. The fifth type is characterised by an achromatic change of the broad band spectrum, which could be attributed to a different mechanism, likely involving differential Doppler boosting caused by geometrical effects. Here we present the classification, the assumed physical scenario and the results of calculations that have been performed for the spectral evolution of flares.

The F-GAMMA program is a coordinated effort to investigate the physics of Active Galactic Nuclei (AGNs) via multi-frequency monitoring of {\em Fermi} blazars. The current study is concerned with the broad-band radio spectra composed of measurement at ten frequencies between 2.64 and 142 GHz. It is shown that any of the 78 sources studied can be classified in terms of their variability characteristics in merely 5 types of variability. The first four types are dominated by spectral evolution and can be reproduced by a simple two-component system made of the quiescent spectrum of a large scale jet populated with a flaring event evolving according to Marscher & Gear (1985). The last type is characterized by an achromatic change of the broad-band spectrum which must be attributed to a completely different mechanism. Here are presented, the classification, the assumed physical system and the results of simulations that have been conducted.

The analysis of the radio light curves of the blazar CTA102 during its 2006 flare revealed a possible interaction between a standing shock wave and a traveling one. In order to better understand this highly non-linear process, we used a relativistic hydrodynamic code to simulate the high energy interaction and its related emission. The calculated synchrotron emission from these simulations showed an increase in turnover flux density, $S_{m}$, and turnover frequency, $\nu_{m}$, during the interaction and decrease to its initial values after the passage of the traveling shock wave.

Context: The blazar CTA 102 (z=1.037) underwent a historical radio outburst in April 2006. This event offered a unique chance to study the physical properties of the jet. Aims: We used multifrequency radio and mm observations to analyze the evolution of the spectral parameters during the flare as a test of the shock-in-jet model under these extreme conditions. Methods: For the analysis of the flare we took into account that the flaring spectrum is superimposed on a quiescent spectrum. We reconstructed the latter from archival data and fitted a synchrotron self-absorbed distribution of emission. The uncertainties of the derived spectral parameters were calculated using Monte Carlo simulations. The spectral evolution is modeled by the shock-in-jet model, and the derived results are discussed in the context of a geometrical model (varying viewing angle) and shock-shock interaction. Results: The evolution of the flare in the turnover frequency-turnover flux density plane shows a double peak structure. The nature of this evolution is dicussed in the frame of shock-in-jet models. We discard the generation of the double peak structure in the turnover frequency-turnover flux density plane purely based on geometrical changes (variation of the Doppler factor). The detailed modeling of the spectral evolution favors a shock-shock interaction as a possible physical mechanism behind the deviations from the standard shock-in-jet model.

Centaurus A is the closest active galactic nucleus. High resolution imaging using Very Long Baseline Interferometry (VLBI) enables us to study the spectral and kinematic behavior of the radio jet-counterjet system on sub-parsec scales, providing essential information for jet emission and formation models. Our aim is to study the structure and spectral shape of the emission from the central-parsec region of Cen A. As a target of the Southern Hemisphere VLBI monitoring program TANAMI (Tracking Active Galactic Nuclei with Milliarcsecond Interferometry), VLBI observations of Cen A are made regularly at 8.4 and 22.3 GHz with the Australian Long Baseline Array (LBA) and associated telescopes in Antarctica, Chile, and South Africa. The first dual-frequency images of this source are presented along with the resulting spectral index map. An angular resolution of 0.4 mas x 0.7 mas is achieved at 8.4 GHz, corresponding to a linear scale of less than 0.013 pc. Hence, we obtain the highest resolution VLBI image of Cen A, comparable to previous space-VLBI observations. By combining with the 22.3 GHz image, which has been taken without contributing transoceanic baselines at somewhat lower resolution, we present the corresponding dual-frequency spectral index distribution along the sub-parsec scale jet revealing the putative emission regions for recently detected gamma-rays from the core region by Fermi/LAT. We resolve the innermost structure of the milliarcsecond scale jet and counterjet system of Cen A into discrete components. The simultaneous observations at two frequencies provide the highest resolved spectral index map of an AGN jet allowing us to identify multiple possible sites as the origin of the high energy emission.

We have found evidence for interaction between a standing and a traveling shock in the jet of the blazar CTA 102. Our result is based in the study of the spectral evolution of the turnover frequency-turnover flux density plane. The radio/mm light curves were taken during a major radio outburst in April 2006.

The blazar CTA 102 underwent a major radio flare in April 2006. We used several 15 GHz VLBI observations from the MOJAVE program to investigate the influence of this extreme event on jet kinematics. The result of modeling and analysis lead to the suggestion of an interaction between traveling and standing shocks 0.2 mas away from the VLBI core.

The radio light curve and spectral evolution of the blazar CTA 102 during its 2006 outburst can be rather well explained by the standard shock-in-jet model. The results of a pixel-to-pixel spectral analysis of multi-frequency VLBI images, together with kinematics derived from the MOJAVE survey lead to the picture of an over-pressured jet with respect to the ambient medium. The interaction of a traveling shock wave with a standing one is a possible scenario which could explain the observed spectral behaviour

At a distance of about 3.8 Mpc, the radio galaxy Centaurus A is the closest active galaxy. Therefore it is a key target for studying the innermost regions of active galactic nuclei (AGN). VLBI observations conducted within the framework of the TANAMI program enable us to study the central region of the Cen A jet with some of the highest linear resolutions ever achieved in an AGN. This region is the likely origin of the gamma-ray emission recently detected by the Fermi Large Area Telescope (LAT). TANAMI monitors a sample of radio and gamma-ray selected extragalactic jets south of -30 degrees declination at 8.4 GHz and 22.3 GHz with the Australian Long Baseline Array (LBA) and the transoceanic antennas Hartebeesthoek in South Africa, the 6 m Transportable Integrated Geodetic Observatory (TIGO) in Chile and the 9 m German Antarctic Receiving Station (GARS) in O'Higgins, Antarctica. The highest angular resolution achieved at 8.4 GHz in the case of Cen A is 0.59mas x 0.978mas (natural weighting) corresponding to a linear scale of less than 18 milliparsec. We show images of the first three TANAMI 8.4 GHz observation epochs of the sub-parsec scale jet-counterjet system of Cen A. With a simultaneous 22.3 GHz observation in 2008 November, we present a high resolution spectral index map of the inner few milliarcseconds of the jet probing the putative emission region of gamma-ray-photons.