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.

We present {\tt radpol} - a numerical scheme for integrating multifrequency polarized radiative transfer equations along rays propagating in a curved spacetime. The scheme includes radiative processes such as synchrotron emission, absorption, Faraday rotation and conversion, and, for the first time, relativistic Compton scatterings including effects of light polarization. The scheme is fully covariant and is applicable to model radio-$\gamma$-ray emission and its polarization from, e.g., relativistic jets and accretion flows onto black holes and other exotic objects described in alternative metric theories and modeled semi-analytically or with time-dependent magnetohydrodynamical simulations. We perform a few tests to validate the implemented numerical algorithms that handle light polarization in curved spacetime. We demonstrate application of the scheme to model broadband emission spectra from a relativistically hot, geometrically thick coronal-like inflow around a supermassive black hole where the disk model is realized in a general relativistic magnetohydrodynamical simulation.

We model non-thermal emission spectrum of the extremely sub-Eddington X-ray binary system A0620-00. It is believed that this non-thermal emission is produced by a radiatively inefficient "quiescent" accretion onto a stellar-mass black hole present in the system. We post-process GRMHD simulations with multiwavelength, fully polarized, relativistic radiative transfer calculations to predict broadband spectra and emission polarization levels for a range of electron models and accretion rates. We find that a model with strong coupling of electrons and ions in the accretion disk and accretion rate of only $\dot{M}=3\times10^{-13} [M_\odot yr^{-1}]$ is able to recover the observed X-ray spectral slope as well as the excess of linear polarization detected in the source in NIR/optical bands. Our models constrain the spectral properties of a putative relativistic jet produced in this system. In addition we show that the magnetized winds from our hot accretion flow carry away a small fraction of the orbital angular momentum of the binary unable to explain the observed rapid orbital decay of the system. Similar to the present GRMHD simulations are often used to explain emission from sub-Eddington supermassvie black holes in Sgr A* or M87; the present simulations allow us to test whether some aspects of the quiescent black hole accretion are scale invariant.

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.

We present ALMA observations at 233, 678, and 870 GHz of the Galactic Center black hole, Sagittarius A*. These observations reveal a flat spectrum over this frequency range with spectral index $\alpha \approx -0.3$, where the flux density $S \propto \nu^\alpha$. We model the submm and far infrared spectrum with a one zone synchrotron model of thermal electrons. We infer electron densities $n = (2-5) \times 10^6$ cm$^{-3}$, electron temperatures $T_e = (1-3) \times 10^{11}$ K, and magnetic field strength $B = 10-50$ G. The parameter range can be further constrained using the observed quiescent X-ray luminosity. The flat submm spectrum results in a high electron temperature and implies that the emitting electrons are efficiently heated. We also find that the emission is most likely optically thin at 233 GHz. These results indicate that millimeter and submillimeter wavelength very long baseline interferometry of Sgr A* including those of the Event Horizon Telescope should see a transparent emission region down to event horizon scales.

Millimeter polarimetry of Sgr A* probes the linearly polarized emission region on a scale of $\sim 10$ Schwarzschild radii ($R_S$) as well as the dense, magnetized accretion flow on scales out to the Bondi radius ($\sim 10^5 R_S$) through Faraday rotation. We present here multi-epoch ALMA Band 6 (230 GHz) polarimetry of Sgr A*. The results confirm a mean rotation measure, ${\rm RM} \approx -5 \times 10^5 {\rm\ rad\ m^{-2}}$, consistent with measurements over the past 20 years and support an interpretation of the RM as originating from a radiatively inefficient accretion flow (RIAF) with $\dot{M} \approx 10^{-8} { \rm\ M_{\odot}\ y^{-1} }$. Variability is observed for the first time in the RM on time scales that range from hours to months. The long-term variations may be the result of changes in the line of sight properties in a turbulent accretion flow. Short-term variations in the apparent RM are not necessarily the result of Faraday rotation and may be the result of complex emission and propagatation effects close to the black hole, some of which have been predicted in numerical modeling. We also confirm the detection of circular polarization at a mean value of $-1.1 \pm 0.2 \%$. It is variable in amplitude on time scales from hours to months but the handedness remains unchanged from that observed in past centimeter- and millimeter-wavelength detections. These results provide critical constraints for the analysis and interpretation of Event Horizon Telescope data of Sgr A*, M87, and similar sources.

We describe ${\tt ipole}$, a new public ray-tracing code for covariant, polarized radiative transport. The code extends the ${\tt ibothros}$ scheme for covariant, unpolarized transport using two representations of the polarized radiation field: in the coordinate frame, it parallel transports the coherency tensor; in the frame of the plasma it evolves the Stokes parameters under emission, absorption, and Faraday conversion. The transport step is implemented to be as spacetime- and coordinate- independent as possible. The emission, absorption, and Faraday conversion step is implemented using an analytic solution to the polarized transport equation with constant coefficients. As a result, ${\tt ipole}$ is stable, efficient, and produces a physically reasonable solution even for a step with high optical depth and Faraday depth. We show that the code matches analytic results in flat space, and that it produces results that converge to those produced by Dexter's ${\tt grtrans}$ polarized transport code on a complicated model problem. We expect ${\tt ipole}$ will mainly find applications in modeling Event Horizon Telescope sources, but it may also be useful in other relativistic transport problems such as modeling for the IXPE mission.

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. In this paper, we study the radiative properties of jets produced in axisymmetric 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. 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 $\kappa$-distribution function. Inside the disk, we assume a thermal distribution for the electrons. We calculate jet spectra and synchrotron maps by using the ray tracing code {\tt RAPTOR}, and compare the synthetic observations to observations of Sgr~A*. Finally, we compare numerical models of jets to semi-analytical ones. We find that in the $\kappa$-jet models, the radio-emitting region size, radio flux, and spectral index in NIR/optical bands increase for decreasing values of the $\kappa$ parameter, which corresponds to a larger amount of accelerated electrons. The model with $\kappa = 3.5$, $\eta_{\rm acc}=5-10\%$ (the percentage of electrons that are accelerated), and observing angle $i = 30^{\rm o}$ fits the observed Sgr~A* emission in the flaring state from the radio to the NIR/optical regimes, while $\kappa = 3.5$, $\eta_{\rm acc}< 1\%$, and observing angle $i = 30^{\rm o}$ fit the upper limits in quiescence.

We present the black hole accretion code (BHAC), a new multidimensional general-relativistic magnetohydrodynamics module for the MPI-AMRVAC framework. BHAC has been designed to solve the equations of ideal general-relativistic magnetohydrodynamics in arbitrary spacetimes and exploits adaptive mesh refinement techniques with an efficient block-based approach. Several spacetimes have already been implemented and tested. We demonstrate the validity of BHAC by means of various one-, two-, and three-dimensional test problems, as well as through a close comparison with the HARM3D code in the case of a torus accreting onto a black hole. The convergence of a turbulent accretion scenario is investigated with several diagnostics and we find accretion rates and horizon-penetrating fluxes to be convergent to within a few percent when the problem is run in three dimensions. Our analysis also involves the study of the corresponding thermal synchrotron emission, which is performed by means of a new general-relativistic radiative transfer code, BHOSS. The resulting synthetic intensity maps of accretion onto black holes are found to be convergent with increasing resolution and are anticipated to play a crucial role in the interpretation of horizon-scale images resulting from upcoming radio observations of the source at the Galactic Center.

Non-VLBI measurements of Faraday rotation at mm wavelengths have been used to constrain mass accretion rates ($\mdot$) onto supermassive black holes in the centre of the Milky Way and in the centre of M87. We constructed general relativistic magnetohydrodynamics models for these sources that qualitatively well describe their spectra and radio/mm images invoking a coupled jet-disk system. Using general relativistic polarized radiative transfer, we now also model the observed mm rotation measure (RM) of M87. The models are tied to the observed radio flux, however, electron temperature and accretion rate are degenerate parameters and are allowed to vary. For the inferred low viewing angles of the M87 jet, the RM is low even as the black hole $\mdot$ increases by a factor of $\simeq100$. In jet-dominated models, the observed linear polarization is produced in the forward-jet, while the dense accretion disk depolarizes the bulk of the near-horizon scale emission which originates in the counter-jet. In the jet-dominated models, with increasing $\mdot$ and increasing Faraday optical depth one is progressively sensitive only to polarized emission in the forward-jet, keeping the measured RM relatively constant. The jet-dominated model reproduces a low net-polarization of $\simeq1$ per cent and RMs in agreement with observed values due to Faraday depolarization, however, with $\mdot$ much larger than the previously inferred limit of $9\times10^{-4}\,\mdotu$. All jet-dominated models produce much higher RMs for inclination angles $i\gtrsim30^\circ$, where the line-of-sight passes through the accretion flow, thereby providing independent constraints on the viewing geometry of the M87 jet.

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

Compact radio cores associated with relativistic jets are often observed in both active galactic nuclei and X-ray binaries. Their radiative properties follow some general scaling laws which primarily depend on their masses and accretion rates. However, it has been suggested that the black hole spin can also strongly influence the power and radio flux of these. Here, we attempt to estimate the dependency of the radio luminosity of steady jets launched by accretion disks on black hole mass, accretion rate and spin using numerical simulations. We make use of 3D GRMHD simulations of accretion disks around low-luminosity black holes in which the jet radio emission is produced by the jet sheath. We find that the radio flux increases roughly by a factor of 6 as the back hole spin increases from a~0 to a=0.98. This is comparable to the increase in accretion power with spin, meaning that the ratio between radio jet and accretion power is hardly changing. Although our jet spine power scales as expected for the Blandford-Znajek process, the dependency of jet radio luminosity on the black hole spin is somewhat weaker. Also weakly rotating black holes can produce visible radio jets. The overall scaling of the radio emission with black hole mass and accretion rate is consistent with the scale-invariant analytical models used to explain the fundamental plane of black hole activity. Spin does not introduce a significant scatter in this model. The jet-sheath model can describe well resolved accreting systems, such as SgrA* and M87, as well as the general scaling behavior of low-luminosity black holes. Hence the model should be applicable to a wide range of radio jets in sub-Eddington black holes. The black hole spin has an effect on the production of visible radio jet, but it may not be the main driver to produce visible radio jets. An extension of our findings to powerful quasars remains speculative.

(abridged) The connection between black hole, accretion disk, and radio jet can be best constrained by fitting models to observations of nearby low luminosity galactic nuclei, in particular the well studied sources Sgr~A* and M87. There has been considerable progress in modeling the central engine of active galactic nuclei by an accreting supermassive black hole coupled to a relativistic plasma jet. However, can a single model be applied to a range of black hole masses and accretion rates? Here we want to compare the latest three-dimensional numerical model, originally developed for Sgr A* in the center of the Milky Way, to radio observations of the much more powerful and more massive black hole in M87. We postprocess three-dimensional GRMHD models of a jet-producing radiatively inefficient accretion flow around a spinning black hole using relativistic radiative transfer and ray-tracing to produce model spectra and images. As a key new ingredient to these models, we allow the proton-electron coupling in these simulations depend on the magnetic properties of the plasma. We find that the radio emission in M87 is well described by a combination of a two-temperature accretion flow and a hot single-temperature jet. The model fits the basic observed characteristics of the M87 radio core. The best fit model has a mass-accretion rate of Mdot approx 9x10^{-3} MSUN/YR and a total jet power of P_j \sim 10^{43} erg/s. Emission at 1.3mm is produced by the counter jet close to the event horizon. Its characteristic crescent shape surrounding the black hole shadow could be resolved by future millimeter-wave VLBI experiments. The model was successfully derived from one for the supermassive black hole in center of the Milky Way by appropriately scaling mass and accretion rate. This suggests the possibility that this model could also apply to a larger range of low-luminosity black holes.

We measure the proper motion of the pulsar PSR J1745-2900 relative to the Galactic Center massive black hole, Sgr A*, using the Very Long Baseline Array (VLBA). The pulsar has a transverse velocity of 236 +/- 11 km s^-1 at position angle 22 +/- 2 deg East of North at a projected separation of 0.097 pc from Sgr A*. Given the unknown radial velocity, this transverse velocity measurement does not conclusively prove that the pulsar is bound to Sgr A*; however, the probability of chance alignment is very small. We do show that the velocity and position is consistent with a bound orbit originating in the clockwise disk of massive stars orbiting Sgr A* and a natal velocity kick of <~ 500 km s^-1. An origin among the isotropic stellar cluster is possible but less probable. If the pulsar remains radio-bright, multi-year astrometry of PSR J1745-2900 can detect its acceleration and determine the full three-dimensional orbit. We also demonstrate that PSR J1745-2900 exhibits the same angular broadening as Sgr A* over a wavelength range of 3.6 cm to 0.7 cm, further confirming that the two sources share the same interstellar scattering properties. Finally, we place the first limits on the presence of a wavelength-dependent shift in the position of Sgr A*, i.e., the core shift, one of the expected properties of optically-thick jet emission. Our results for PSR J1745-2900 support the hypothesis that Galactic Center pulsars will originate from the stellar disk and deepens the mystery regarding the small number of detected Galactic Center pulsars.

The emission region and black hole shadow of Sagittarius A*, the supermassive black hole at the Galactic Center, can be probed with millimeter Very Long Baseline Interferometry. Our goal is to probe the geometry of the emitting plasma around Sgr A* by using modeled mm-VLBI closure phase calculations at 1.3 mm and to constrain the observer's inclination angle and position angle of the black hole spin axis. We have simulated images for three different models of the emission of Sgr A*: an orbiting spot, a disk model, and a jet model. The orbiting spot model was used as a test case scenario, while the disk and jet models are physically driven scenarios based on standard three-dimensional general relativistic magnetohydrodynamic simulations of hot accretion flows. Our results are compared to currently available closure phase observational limits. Our results indicate that more models with closer to edge-on viewing angles are consistent with observational limits. In general, jet and disk geometries can reproduce similar closure phases for different sets of viewing and position angles. Consequently, the favored black hole spin orientation and its magnitude are strongly model dependent. We find that both the jet and the disk models can explain current VLBI limits. We conclude that new observations at 1.3 mm and possibly at longer wavelengths including other triangles of VLBI baselines are necessary to interpret Sgr A* emission and the putative black hole spin parameters.

We analyze 134 ks Chandra ACIS-I observations of the Galactic Centre (GC) performed in July 2011. The X-ray image with the field of view $17' \times 17'$ contains the hot plasma surrounding the Sgr~A*. The obtained surface brightness map allow us to fit Bondi hot accretion flow to the innermost hot plasma around the GC. We have fitted spectra from region up to $5"$ from Sgr~A* using a thermal bremsstrahlung model and four Gaussian profiles responsible for K$_{\alpha}$ emission lines of Fe, S, Ar, and Ca. The X-ray surface brightness profile up to $3"$ from Sgr~A* found in our data image, was successfully fitted with the dynamical model of Bondi spherical accretion. By modelling the surface brightness profile, we derived the temperature and number density profiles in the vicinity of the black hole. The best fitted model of spherical Bondi accretion shows that this type of flow works only up to $3"$ and implies outer plasma density and temperature to be: $n_{\rm e}^{\rm out}=18.3 \pm {0.1}$ cm$^{-3}$ and $T_{\rm e}^{\rm out}= 3.5 \pm {0.3}$ keV respectively. We show that the Bondi flow can reproduce observed surface brightness profile up to $3"$ from Sgr~A* in the Galactic Center. This result strongly suggests the position of stagnation radius in the complicated dynamics around GC. The Faraday rotation computed from our model towards the pulsar PSR J1745-2900 near the GC agrees with the observed one, recently reported.

Radio and mm-wavelength observations of Sagittarius A* (Sgr A*), the radio source associated with the supermassive black hole at the center of our Galaxy, show that it behaves as a partially self-absorbed synchrotron-emitting source. The measured size of Sgr A* shows that the mm-wavelength emission comes from a small region and consists of the inner accretion flow and a possible collimated outflow. Existing observations of Sgr A* have revealed a time lag between light curves at 43 GHz and 22 GHz, which is consistent with a rapidly expanding plasma flow and supports the presence of a collimated outflow from the environment of an accreting black hole. Here we wish to measure simultaneous frequency-dependent time lags in the light curves of Sgr A* across a broad frequency range to constrain direction and speed of the radio-emitting plasma in the vicinity of the black hole. Light curves of Sgr A* were taken in May 2012 using ALMA at 100 GHz using the VLA at 48, 39, 37, 27, 25.5, and 19 GHz. As a result of elevation limits and the longitude difference between the stations, the usable overlap in the light curves is approximately four hours. Although Sgr A* was in a relatively quiet phase, the high sensitivity of ALMA and the VLA allowed us to detect and fit maxima of an observed minor flare where flux density varied by ~10%. The fitted times of flux density maxima at frequencies from 100 GHz to 19 GHz, as well as a cross-correlation analysis, reveal a simple frequency-dependent time lag relation where maxima at higher frequencies lead those at lower frequencies. Taking the observed size-frequency relation of Sgr A* into account, these time lags suggest a moderately relativistic (lower estimates: 0.5c for two-sided, 0.77c for one-sided) collimated outflow.

Radiatively inefficient accretion flows (RIAFs) are believed to power supermassive black holes (SMBH) in the underluminous cores of galaxies. Such black holes are typically accompanied by flat-spectrum radio cores indicating the presence of moderately relativistic jets. One of the best constrained RIAFs is associated with the SMBH in the Galactic center, Sgr A*. Since the plasma in RIAFs is only weakly collisional, the dynamics and the radiative properties of these systems are very uncertain. Here we want to study the impact of varying electron temperature on the appearance of accretion flows and jets. Using 3-D GRMHD accretion flow simulations, we use ray tracing methods to predict spectra and radio images of RIAFs allowing for different electron heating mechanisms in the in- and outflowing parts of the simulations. We find that small changes in the electron temperature can result in dramatic differences in the relative dominance of jets and accretion flows. Application to Sgr A* shows that radio spectrum and size of this source can be well reproduced with a model where electrons are more efficiently heated in the jet. The X-ray emission is sensitive to the electron heating mechanism in the jets and disk and therefore X-ray observations put strong constraints on electron temperatures and geometry of the accretion flow and jet. For Sgr A*, the jet model also predicts a significant frequency-dependent core shift which could place independent constraints on the model once measured accurately. We conclude that more sophisticated models for electron distribution functions are crucial for constraining GRMHD simulations with actual observations. For Sgr A*, the radio appearance may well be dominated by the outflowing plasma. Nonetheless, at the highest radio frequencies, the shadow of the event horizon should still be detectable with future Very Long Baseline Interferometric observations.

The supermassive black hole in the center of the Milky Way, Sgr A*, displays a nearly flat radio spectrum that is typical for jets in active galactic nuclei. Indeed, time-dependent magnetized models of radiatively inefficient accretion flows (RIAFs), which are commonly used to explain the millimeter, near-infrared, and X-ray emission of Sgr A*, often also produce jet-like outflows. However, the emission from these models has so far failed to reproduce the flat radio spectrum. Aims: We investigate whether current accretion simulations can produce the compact flat spectrum emission by simply using a different prescription for the heating of the radiating particles in the jet. Methods: We studied the radiative properties of accretion flows onto a black hole produced in time-dependent general-relativistic magnetohydrodynamic (GRMHD) simulations. A crucial free parameter in these models has always been the electron temperature, and here we allowed for variations in the proton-to-electron temperature ratios in the jet and disk. Results: We found that the flat spectrum is readily reproduced by a standard GRMHD model if one has an almost isothermal jet coupled to a two-temperature accretion flow. The low-frequency radio mission comes from the outflowing sheath of matter surrounding the strongly magnetized nearly empty jet. The model is consistent with the radio sizes and spectrum of Sgr A*. Conclusions: Hence, GRMHD models of accreting black holes can in principle naturally reproduce jets that match observed characteristics. For Sgr A* the model fit to the spectrum predicts higher mass-accretion rates when a jet is included than without a jet. Hence, the impact of the recently discovered G2 cloud that is expected to be accreted onto Sgr A* might be less severe than currently thought.

We calculate the structure and short-term evolution of a gamma ray burst central engine in the form of a turbulent torus accreting onto a stellar mass black hole. Our models apply to the short gamma ray burst events, in which a remnant torus forms after the neutron star-black hole or a double neutron star merger and is subsequently accreted. We study the 2-dimensional, relativistic models and concentrate on the effects of black hole and flow parameters as well as the neutrino cooling. We compare the resulting structure and neutrino emission to the results of our previous 1-dimensional simulations. We find that the neutrino cooled torus launches a powerful mass outflow, which contributes to the total neutrino luminosity and mass loss from the system. The neutrino luminosity may exceed the Blandford-Znajek luminosity of the polar jets and the subsequent annihilation of neutrino-antineutrino pairs will provide an additional source of power to the GRB emission.

(Abridged) We study stability of gas accretion in Active Galactic Nuclei. Our grid based simulations cover a radial range from 0.1 to 200 pc. Here, as in previous studies by our group, we include gas radiative cooling as well as heating by a sub-Eddington X-ray source near the central supermassive black hole of 10^8 M_{\odot}. Our theoretical estimates and simulations show that for the X-ray luminosity L_X \sim 0.008 L_{Edd}, the gas is thermally and convectivelly unstable within the computational domain. In the simulations, we observe that very tiny fluctuations in an initially smooth, spherically symmetric, accretion flow, grow first linearly and then non-linearly. Consequently, an initially one-phase flow relatively quickly transitions into a two-phase/cold-hot accretion flow. For L_X = 0.015 L_{Edd} or higher, the cold clouds continue to accrete but in some regions of the hot phase, the gas starts to move outward. For L_X < 0.015 L_{Edd}, the cold phase contribution to the total mass accretion rate only moderately dominates over the hot phase contribution. This result might have some consequences for cosmological simulations of the so-called AGN feedback problem. Our simulations confirm the previous results of Barai et al. (2012) who used smoothed particle hydrodynamic simulations to tackle the same problem. However here, because we use a grid based code to solve equations in 1-D and 2-D, we are able to follow the gas dynamics at much higher spacial resolution and for longer time in comparison to the 3-D SPH simulations. One of new features revealed by our simulations is that the cold condensations in the accretion flow initially form long filaments, but at the later times, those filaments may break into smaller clouds advected outwards within the hot outflow. These simulations may serve as an attractive model for the so-called Narrow Line Region in AGN.

In accretion-based models for Sgr A* the X-ray, infrared, and millimeter emission arise in a hot, geometrically thick accretion flow close to the black hole. The spectrum and size of the source depend on the black hole mass accretion rate $\dot{M}$. Since Gillessen et al. have recently discovered a cloud moving toward Sgr A* that will arrive in summer 2013, $\dot{M}$ may increase from its present value $\dot{M}_0$. We therefore reconsider the "best-bet" accretion model of Moscibrodzka et al., which is based on a general relativistic MHD flow model and fully relativistic radiative transfer, for a range of $\dot{M}$. We find that for modest increases in $\dot{M}$ the characteristic ring of emission due to the photon orbit becomes brighter, more extended, and easier to detect by the planned Event Horizon Telescope submm VLBI experiment. If $\dot{M} \gtrsim 8 \dot{M}_0$ this "silhouette of the black hole will be hidden beneath the synchrotron photosphere at 230 GHz, and for $\dot{M} \gtrsim 16 \dot{M}_0$ the silhouette is hidden at 345 GHz. We also find that for $\dot{M} > 2 \dot{M}_0$ the near-horizon accretion flow becomes a persistent X-ray and mid-infrared source, and in the near-infrared Sgr A* will acquire a persistent component that is brighter than currently observed flares.

Gamma Ray Bursts (GRB) are the extremely energetic transient events, visible from the most distant parts of the Universe. They are most likely powered by accretion on the hyper-Eddington rates that proceeds onto a newly born stellar mass black hole. This central engine gives rise to the most powerful, high Lorentz factor jets that are responsible for energetic gamma ray emission. We investigate the accretion flow evolution in GRB central engine, using the 2D MHD simulations in General Relativity. We compute the structure and evolution of the extremely hot and dense torus accreting onto the fast spinning black hole, which launches the magnetized jets. We calculate the chemical structure of the disk and account for neutrino cooling. Our preliminary runs apply to the short GRB case (remnant torus accreted after NS-NS or NS-BH merger). We estimate the neutrino luminosity of such an event for chosen disk and central BH mass

Electron-positron pairs may be produced near accreting black holes by a variety of physical processes, and the resulting pair plasma may be accelerated and collimated into a relativistic jet. Here we use a self-consistent dynamical and radiative model to investigate pair production by \gamma\gamma collisions in weakly radiative accretion flows around a black hole of mass M and accretion rate \dot{M}. Our flow model is drawn from general relativistic magnetohydrodynamic simulations, and our radiation field is computed by a Monte Carlo transport scheme assuming the electron distribution function is thermal. We argue that the pair production rate scales as r^{-6} M^{-1} \dot{M}^{6}. We confirm this numerically and calibrate the scaling relation. This relation is self-consistent in a wedge in M, \dot{M} parameter space. If \dot{M} is too low the implied pair density over the poles of the black hole is below the Goldreich-Julian density and \gamma\gamma pair production is relatively unimportant; if \dot{M} is too high the models are radiatively efficient. We also argue that for a power-law spectrum the pair production rate should scale with the observables L_X \equiv X-ray luminosity and M as L_X^2 M^{-4}. We confirm this numerically and argue that this relation likely holds even for radiatively efficient flows. The pair production rates are sensitive to black hole spin and to the ion-electron temperature ratio which are fixed in this exploratory calculation. We finish with a brief discussion of the implications for Sgr A* and M87.

In this paper we review and discuss some of the intriguing properties of the Galactic Center supermassive black hole candidate Sgr A*. Of all possible black hole sources, the event horizon of Sgr A*, subtends the largest angular scale on the sky. It is therefore a prime candidate to study and image plasma processes in strong gravity and it even allows imaging of the shadow cast by the event horizon. Recent mm-wave VLBI and radio timing observations as well as numerical GRMHD simulations now have provided several breakthroughs that put Sgr A* back into the focus. Firstly, VLBI observations have now measured the intrinsic size of Sgr A* at multiple frequencies, where the highest frequency measurements have approached the scale of the black hole shadow. Moreover, measurements of the radio variability show a clear time lag between 22 GHz and 43 GHz. The combination of size and timing measurements, allows one to actually measure the flow speed and direction of magnetized plasma at some tens of Schwarzschild radii. This data strongly support a moderately relativistic outflow, consistent with an accelerating jet model. This is compared to recent GRMHD simulation that show the presence of a moderately relativistic outflow coupled to an accretion flow Sgr A*. Further VLBI and timing observations coupled to simulations have the potential to map out the velocity profile from 5-40 Schwarzschild radii and to provide a first glimpse at the appearance of a jet-disk system near the event horizon. Future submm-VLBI experiments would even be able to directly image those processes in strong gravity and directly confirm the presence of an event horizon.

We present various instability mechanisms in the accreting black hole systems which might indicate at the connection between the accretion disk and jet. The jets observed in microquasars can have a persistent or blobby morphology. Correlated with the accretion luminosity, this might provide a link to the cyclic outbursts of the disk. Such duty-cycle type of behavior on short timescales results from the thermal instability caused by the radiation pressure domination. The same type of instability may explain the cyclic radioactivity of the supermassive black hole systems. The somewhat longer timescales are characteristic for the instability caused by the partial hydrogen ionization. The distortions of the jet direction and complex morphology of the sources can be caused by precession of the disk-jet axis.

We review results from general relativistic axisymmetric magnetohydrodynamic simulations of accretion in Sgr A*. We use general relativistic radiative transfer methods and to produce a broad band (from millimeter to gamma-rays) spectrum. Using a ray tracing scheme we also model images of Sgr A* and compare the size of image to the VLBI observations at 230 GHz. We perform a parameter survey and study radiative properties of the flow models for various black hole spins, ion to electron temperature ratios, and inclinations. We scale our models to reconstruct the flux and the spectral slope around 230 GHz. The combination of Monte Carlo spectral energy distribution calculations and 230 GHz image modeling constrains the parameter space of the numerical models. Our models suggest rather high black hole spin ($a_*\approx 0.9$), electron temperatures close to the ion temperature ($T_i/T_e \sim 3$) and high inclination angles ($i \approx 90 \deg$).

We study a slightly rotating accretion flow onto a black hole, using the fully three dimensional (3-D)numerical simulations. We consider hydrodynamics of an inviscid flow, assuming a spherically symmetric density distribution at the outer boundary and a small, latitude-dependent angular momentum. We investigate the role of the adiabatic index and gas temperature, and the flow behaviour due to non-axisymmetric effects. Our 3-D simulations confirm axisymmetric results: the material that has too much angular momentum to be accreted forms a thick torus near the equator and the mass accretion rate is lower than the Bondi rate. In our previous study of the 3-D accretion flows, for gamma=5/3, we found that the inner torus precessed, even for axisymmetric conditions at large radii. The present study shows that the inner torus precesses also for other values of the adiabatic index: gamma=4/3, 1.2 and 1.01. However, the time for the precession to set increases with decreasing gamma. In particular, for gamma=1.01 we find that depending on the outer boundary conditions, the torus may shrink substantially due to the strong inflow of the non-rotating matter and the precession will have insufficient time to develop. On the other hand, if the torus is supplied by the continuous inflow of the rotating material from the outer radii, its inner parts will eventually tilt and precess, as it was for the larger gamma's.

Using flow models based on axisymmetric general relativistic magnetohydrodynamics (GRMHD) simulations, we construct radiative models for sgr A*. Spectral energy distributions that include the effects of thermal synchrotron emission and absorption, and Compton scattering, are calculated using a Monte Carlo technique. Images are calculated using a ray-tracing scheme. All models are scaled so that the 230 GHz flux density is 3.4 Jy. The key model parameters are the dimensionless black hole spin a*, the inclination i, and the ion-to-electron temperature ratio Ti/Te. We find that: (1) models with Ti/Te=1 are inconsistent with the observed submillimeter spectral slope; (2) the X-ray flux is a strongly increasing function of a*; (3) the X-ray flux is a strongly increasing function of i; (4) 230 GHz image size is a complicated function of i, a*, and Ti/Te, but the Ti/Te = 10 models are generally large and at most marginally consistent with the 230 GHz VLBI data; (5) for models with Ti/Te=10 and i=85 deg the event horizon is cloaked behind a synchrotron photosphere at 230 GHz and will not be seen by VLBI, but these models overproduce NIR and X-ray flux; (6) in all models whose SEDs are consistent with observations the event horizon is uncloaked at 230 GHz; (7) the models that are most consistent with the observations have a* \sim 0.9. We finish with a discussion of the limitations of our model and prospects for future improvements.

We describe a Monte Carlo radiative transport code intended for calculating spectra of hot, optically thin plasmas in full general relativity. The version we describe here is designed to model hot accretion flows in the Kerr metric and therefore incorporates synchrotron emission and absorption, and Compton scattering. The code can be readily generalized, however, to account for other radiative processes and an arbitrary spacetime. We describe a suite of test problems, and demonstrate the expected $N^{-1/2}$ convergence rate, where $N$ is the number of Monte Carlo samples. Finally we illustrate the capabilities of the code with a model calculation, a spectrum of the slowly accreting black hole Sgr A* based on data provided by a numerical general relativistic MHD model of the accreting plasma.

We study how axisymmetric magnetohydrodynamical (MHD) accretion flows depend on gamma adiabatic index in the polytropic equation of state. This work is an extension of Moscibrodzka & Proga (2008), where we investigated the gamma dependence of 2-D Bondi-like accretion flows in the hydrodynamical (HD) limit. Our main goal is to study if simulations for various gamma can give us insights into to the problem of various modes of accretion observed in several types of accretion systems such as black hole binaries (BHB), active galactic nuclei (AGN), and gamma-ray bursts (GRBs). We find that for gamma >~ 4/3, the fast rotating flow forms a thick torus that is supported by rotation and gas pressure. As shown before for gamma=5/3, such a torus produces a strong, persistent bipolar outflow that can significantly reduce the polar funnel accretion of a slowly rotating flow. For low gamma, close to 1, the torus is thin and is supported by rotation. The thin torus produces an unsteady outflow which is too weak to propagate throughout the polar funnel inflow. Compared to their HD counterparts, the MHD simulations show that the magnetized torus can produce an outflow and does not exhibit regular oscillations. Generally, our simulations demonstrate how the torus thickness affects the outflow production. They also support the notion that the geometrical thickness of the torus correlates with the power of the torus outflow. Our results, applied to observations, suggest that the torus ability to radiatively cool and become thin can correspond to a suppression of a jet as observed in the BHB during a transition from a hard/low to soft/high spectral state and a transition from a quiescent to hard/low state in AGN.

Sgr A* is the closest massive black hole and can be observed with the highest angular resolution. Nevertheless, our current understanding of the accretion process in this source is very poor. The inflow is almost certainly of low radiative efficiency and it is accompanied by a strong outflow and the flow is strongly variable but the details of the dynamics are unknown. Even the amount of angular momentum in the flow is an open question. Here we argue that low angular momentum scenario is better suited to explain the flow variability. We present a new hybrid model which describes such a flow and consists of an outer spherically symmetric Bondi flow and an inner axially symmetric flow described through MHD simulations. The assumed angular momentum of the matter is low, i.e. the corresponding circularization radius in the equatorial plane of the flow is just above the innermost stable circular orbit in pseudo-Newtonian potential. We compare the radiation spectrum from such a flow to the broad band observational data for Sgr A*.

We report on next phase of our study of rotating accretion flows onto black holes. We consider hydrodynamical (HD) accretion flows with a spherically symmetric density distribution at the outer boundary but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum. We study accretion flows by means of numerical two-dimensional, axisymmetric, HD simulations for variety of the adiabatic index, $\gamma$ and the gas temperature at infinity, $c_\infty$. Our work is an extension of work done by Proga & Begelman who consider models for only $\gamma=5/3$. Our main result is that the flow properties such as the topology of the sonic surface and time behavior strongly depend on $\gamma$ but little on $c_\infty$. In particular, for $1 < \gamma < 5/3$, the mass accretion rate shows large amplitude, slow time-variability which is a result of mixing between slow and fast rotating gas. This temporal behavior differs significantly from that in models with $\gamma\simless 5/3$ where the accretion rate is relatively constant and from that in models with $\gamma\simgreat 1$ where the accretion exhibits small amplitude quasi-periodic oscillations. The key parameter responsible for the differences is the sound speed of the accretion flow which in turn determines whether the flow is dominated by gas pressure, radiation pressure or rotation. Despite these differences the time-averaged mass accretion rate in units of the corresponding Bondi rate is a weak function of $\gamma$ and $c_\infty$.

Sgr A* is a source of strongly variable emission in several energy bands. It is generally agreed that this emission comes from the material surrounding the black hole which is either falling in or flowing out. The activity must be driven by accretion but the character of accretion flow in this object is an open question. We suggest that the inflow is dominated by the relatively low angular momentum material originating in one of the nearby group of stars. Such material flows in directly towards the black hole up to the distance of order of ten Schwarzschild radii or less, where it hits the angular momentum barrier which leads naturally to a flow variability. We study both the analytical and the numerical solutions for the flow dynamics, and we analyze the radiation spectra in both cases using the Monte Carlo code to simulate the synchrotron, bremsstrahlung and the Compton scattering. Our model roughly reproduces the broad band spectrum of Sgr A* and its variability if we allow for a small fraction of energy to be converted to non-thermal population of electrons. It is also consistent (for a range of viewing angles) with the strong constraints on the amount of circumnuclear material imposed by the measurements of the Faraday rotation.

Numerical simulations of MHD accretion flows in the vicinity of a supermasssive black hole provide important insights to the problem of why and how systems -- such as the Galactic Center -- are underluminous and variable. To access applicability of such flows to real objects, we examine the dynamical MHD studies with computations of the time dependent radiation spectra predicted by the simulations. We apply Monte Carlo methods to calculate spectra predicted by the time-dependent model of an axisymmetic MHD flow accreting onto a black hole presented by Proga and Begelman. Our calculations show that variability in an accretion flow is not always reflected in the corresponding spectra, at least not in all wavelengths. We find no one-to-one correspondence between the accretion state and the predicted spectrum

The material accreting onto Sgr A* most probably comes from the nearby stars. We analyze the pattern of this flow at distances of a fraction of a parsec and we argue that the net angular momentum of this material is low but non-negligible, and the initially supersonic disk accretion changes into subsonic flow with constant angular momentum. Next we estimate the flow parameters at a distance $R_{BHL}$ from the black hole and we argue that for the plausible parameter range the accretion flow is non-stationary. The inflow becomes supersonic at distance of $\sim 10^4 R_g$ but the solution does not continue below the horizon and the material piles up forming a torus, or a ring, at a distance of a few up to tens of Schwarzchild radii. Such a torus is known to be unstable and may explain strong variability of the flow in Sgr A*. Our considerations show that the temporary formation of such a torus seems to be unavoidable. Our best fitting model predicts a rather large accretion rate of around $4 \cdot 10^{-6} M_{\odot}/yr$ directly on Sgr A*. We argue that magnetic fields in the flow are tangled and this allows our model to overcome the disagreement with the Faraday rotation limits.

We consider the sample of weakly active galaxies situated in 'Local Universe' collected in the paper of Pellegrini (2005) with inferred accretion efficiencies from $10^{-2}$ to $10^{-7}$. We apply a model of spherically symmetrical Bondi accretion for given parameters ($M_{BH}$,$T_{\infty}$,$\rho_{\infty}$,) taken from observation. We calculate spectra emitted by the gas accreting onto its central objects using Monte Carlo method including synchrotron and bremsstrahlung photons as seed photons. We compare our results with observed nuclear X-ray luminosities $L_{X,nuc}$ (0.3-10 keV) of the sample. Model is also tested for different external medium parameters ($\rho_{\infty}$ and $T_{\infty}$) and different free parameters of the model. Our model is able to explain most of the observed nuclear luminosities $L_X$ under an assumption that half of the compresion energy is transfered directly to the electrons.