Large-amplitude Sgr A* near-infrared flares result from energy injection into electrons near the black hole event horizon. Astrometry data show continuous rotation of the emission region during bright flares, and corresponding rotation of the linear polarization angle. One broad class of physical flare models invokes magnetic reconnection. Here we show that such a scenario can arise in a general relativistic magnetohydrodynamic simulation of a magnetically arrested disc. Saturation of magnetic flux triggers eruption events, where magnetically dominated plasma is expelled from near the horizon and forms a rotating, spiral structure. Dissipation occurs via reconnection at the interface of the magnetically dominated plasma and surrounding fluid. This dissipation is associated with large increases in near-infrared emission in models of Sgr A*, with durations and amplitudes consistent with the observed flares. Such events occur at roughly the timescale to re-accumulate the magnetic flux from the inner accretion disc, 10h for Sgr A*. We study near-infrared observables from one sample event to show that the emission morphology tracks the boundary of the magnetically dominated region. As the region rotates around the black hole, the near-infrared centroid and linear polarization angle both undergo continuous rotation, similar to the behavior seen in Sgr A* flares.

Infrared interferometry is a new frontier for precision ground based observing, with new instrumentation achieving milliarcsecond (mas) spatial resolutions for faint sources, along with astrometry on the order of 10 microarcseconds. This technique has already led to breakthroughs in the observations of the supermassive black hole at the Galactic centre and its orbiting stars, AGN, and exo-planets, and can be employed for studying X-ray binaries (XRBs), microquasars in particular. Beyond constraining the orbital parameters of the system using the centroid wobble and spatially resolving jet discrete ejections on mas scales, we also propose a novel method to discern between the various components contributing to the infrared bands: accretion disk, jets and companion star. We demonstrate that the GRAVITY instrument on the Very Large Telescope Interferometer (VLTI) should be able to detect a centroid shift in a number of sources, opening a new avenue of exploration for the myriad of transients expected to be discovered in the coming decade of radio all-sky surveys. We also present the first proof-of-concept GRAVITY observation of a low-mass X-ray binary transient, MAXI J1820+070, to search for extended jets on mas scales. We place the tightest constraints yet via direct imaging on the size of the infrared emitting region of the compact jet in a hard state XRB.

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

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 Galactic center black hole candidate Sgr A* is the best target for studies of low-luminosity accretion physics, including with near-infrared and submillimeter wavelength long baseline interferometry experiments. Here we compare images and spectra generated from a parameter survey of general relativistic MHD simulations to a set of radio to near-infrared observations of Sgr A*. Our models span the limits of weak and strong magnetization and use a range of sub-grid prescriptions for electron heating. We find two classes of scenarios can explain the broad shape of the submillimeter spectral peak and the highly variable near-infrared flaring emission. Weakly magnetized "disk-jet" models where most of the emission is produced near the jet wall, consistent with past work, as well as strongly magnetized (magnetically arrested disk) models where hot electrons are present everywhere. Disk-jet models are strongly depolarized at submillimeter wavelengths as a result of strong Faraday rotation, inconsistent with observations of Sgr A*. We instead favor the strongly magnetized models, which provide a good description of the median and highly variable linear polarization signal. The same models can also explain the observed mean Faraday rotation measure and potentially the polarization signals seen recently in Sgr A* near-infrared flares.

We analyze two 3D general-relativistic magnetohydrodynamic accretion simulations in the context of how they would manifest in Event Horizon Telescope (EHT) observations of supermassive black holes. The two simulations differ only in whether the initial angular momentum of the plasma is aligned with the rapid (a = 0.9) spin of the black hole. Both have low net magnetic flux. Ray tracing is employed to generate resolved images of the synchrotron emission. When using parameters appropriate for Sgr A* and assuming a viewing angle aligned with the black hole spin, we find the most prominent difference is that the central shadow in the image is noticeably eccentric in tilted models, with the ring of emission relatively unchanged. Applying this procedure to M87 with a viewing angle based on the large-scale jet, we find that adding tilt increases the angular size of the ring for fixed black hole mass and distance, while at the same time increasing the number of bright spots in the image. Our findings illustrate observable features that can distinguish tilted from aligned flows. They also show that tilted models can be viable for M87, and that not accounting for tilt can bias inferences of physical parameters. Future modeling of horizon-scale observations should account for potential angular momentum misalignment, which is likely generic at the low accretion rates appropriate for EHT targets.

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

We analyze the two brightest Chandra X-ray flares detected from Sagittarius A*, with peak luminosities more than 600 x and 245 x greater than the quiescent X-ray emission. The brightest flare has a distinctive double-peaked morphology --- it lasts 5.7 ksec ($\sim 2$ hours), with a rapid rise time of 1500 sec and a decay time of 2500 sec. The second flare lasts 3.4 ksec, with rise and decay times of 1700 sec and 1400 sec. These luminous flares are significantly harder than quiescence: the first has a power law spectral index $\Gamma = 2.06\pm 0.14$ and the second has $\Gamma = 2.03\pm 0.27$, compared to $\Gamma = 3.0\pm0.2$ for the quiescent accretion flow. These spectral indices (as well as the flare hardness ratios) are consistent with previously-detected Sgr A* flares, suggesting that bright and faint flares arise from similar physical processes. Leveraging the brightest flare's long duration and high signal-to-noise, we search for intraflare variability and detect excess X-ray power at a frequency of $\nu \approx 3$ mHz, but show that it is an instrumental artifact and not of astrophysical origin. We find no other evidence (at the 95% confidence level) for periodic or quasi-periodic variability in either flares' time series. We also search for non-periodic excess power but do not find compelling evidence in the power spectrum. Bright flares like these remain our most promising avenue for identifying Sgr A*'s short timescale variability in the X-ray, which may probe the characteristic size scale for the X-ray emission region.

We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region.

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.

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation.

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 use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity.

Precise measurements of the S-stars orbiting SgrA* have set strong constraints on the nature of the compact object at the centre of the Milky Way. The presence of a black hole in that region is well established, but its neighboring environment is still an open debate. In that respect, the existence of dark matter in that central region may be detectable due to its strong signatures on the orbits of stars: the main effect is a Newtonian precession which will affect the overall pericentre shift of S2, the latter being a target measurement of the GRAVITY instrument. The exact nature of this dark matter (e.g., stellar dark remnants or diffuse dark matter) is unknown. This article assumes it to be an scalar field of toroidal distribution, associated with ultra-light dark matter particles, surrounding the Kerr black hole. Such a field is a form of "hair" expected in the context of superradiance, a mechanism that extracts rotational energy from the black hole. Orbital signatures for the S2 star are computed and shown to be detectable by GRAVITY. The scalar field can be constrained because the variation of orbital elements depends both on the relative mass of the scalar field to the black hole and on the field mass coupling parameter.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

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.

We analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical "changing look" quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|\beta_{He}-\beta_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.

In a dynamically relaxed cluster around a massive black hole a dense stellar cusp of old stars is expected to form. Previous observations showed a relative paucity of red giant stars within the central 0.5 pc in the Galactic Center. By co-adding spectroscopic observations taken over a decade, we identify new late-type stars, including the first five warm giants (G2-G8III), within the central 1 arcsec 2 (0.04 {\times} 0.04 pc^2) of the Galaxy. Our findings increase the number of late-type stars to 21, of which we present deep spectra for 16. The updated star count, based on individual spectral classification, is used to reconstruct the surface density profile of giant stars. Our study, for the first time, finds a cusp in the surface number density of the spectroscopically identified old (>3 Gyr) giants population (m K<17) within 0.02-0.4 pc described by a single power law with an exponent {\Gamma}= 0.34 {\pm} 0.04.

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

The Galactic Center black hole Sgr A* is the archetypical example of an underfed massive black hole. The extremely low accretion rate can be understood in radiatively inefficient accretion flow models. Testing those models has proven to be difficult due to the lack of suitable probes. Radio and submm polarization measurements constrain the flow very close to the event horizon. X-ray observations resolving the Bondi radius yield an estimate roughly four orders of magnitude further out. Here, we present a new, indirect measurement of the accretion flow density at intermediate radii. We use the dynamics of the gas cloud G2 to probe the ambient density. We detect the presence of a drag force slowing down G2 with a statistical significance of approx 9 {\sigma}. This probes the accretion flow density at around 1000 Schwarzschild radii and yields a number density of approx. 4 x 10^3 cm^-3. Self-similar accretion models where the density follows a power law radial profile between the inner zone and the Bondi radius have predicted similar values.

The microquasar SS 433 is well-known for its precessing, relativistic baryonic jets. Depending on their heating mechanism, the optical jet bullets may serve as a probe of the collimated radiation coming from the inner region close to the compact object. The optical interferometer VLTI/GRAVITY has allowed to spatially resolved the optical jets in SS 433 for the first time. We present here the second such observation taken over three nights in July 2017. In addition, we use the multi-wavelength XSHOOTER spectrograph at VLT to study the optical bullets in SS 433 in detail. GRAVITY reveals elongated exponential-like spatial profiles for the optical jets, suggestive of a heating mechanism acting throughout a long portion of the jet and naturally explained by photoionization by the collimated radiation. We also spatially resolve the movement of the optical bullets for the first time, detecting extended jet components corresponding to previous ejections. \texttt{Cloudy} photoionization models explain both the spatial intensity profiles measured with GRAVITY and the line ratios from XSHOOTER, and constrain the properties of the optical bullets and the ionizing radiation. We find that the latter must peak in the UV with an isotropic luminosity (as inferred by a face on observer) $\approx 10^{41}$ erg/s. Provided that the X-ray SED is sufficiently hard, the collimated X-ray luminosity could still be high enough so that the face on observer would see SS 433 as an ULX ($L_X \lesssim 10^{40}$ erg/s) and it would still be compatible with the H/He/He+ ionization balance of the optical bullets. The kinetic power in the optical jets is constrained to $3-20 \times 10^{38}$ erg/s, and the extinction in the optical jets to $A_V = 6.7 \pm 0.1$. We suggest there may be substantial $A_V \gtrsim 1$ and structured circumstellar extinction in SS 433, likely arising from dust formed in equatorial outflows.

The microquasar SS 433 is the only known steady supercritical accretor in the Galaxy. It is well-known for its relativistic baryonic jets, but the system also drives equatorial outflows. These have been routinely detected in radio images, and components associated with a circumbinary disk have also been suggested in optical emission lines. We wish to spatially resolve the regions producing the stationary emission lines of SS 433 to shed light on its circumbinary structure and outflows. With an estimated binary orbit size $\lesssim 0.1 \text{ mas}$, this requires optical interferometry. We use the optical interferometer VLTI+GRAVITY to spatially resolve SS 433 in the near-infrared K band at high spectral resolution ($R\approx 4000$) on three nights in July 2017. The stationary Br$\gamma$ line is clearly dominated by an extended $\sim 1 \text{ mas} \sim 5 \text{ AU}$ circumbinary structure perpendicular to the jets and with a strong rotation component. The rotation direction is retrograde relative to the jet precession, in accordance with the slaved disk precession model. The structure has a very high specific angular momentum and is too extended to be a stable circumbinary disk in Keplerian rotation; interpreting it as such leads to a very high enclosed mass $M \gtrsim 400 M_{\odot}$. We instead interpret it as the centrifugal ejection of the circumbinary disk, with the implication that there must be an efficient transfer of specific angular momentum from the binary to the disk. We suggest that the equatorial outflows sometimes seen in radio images result from similar episodes of circumbinary disk centrifugal ejection. In addition to the equatorial structure, we find a very extended $\sim 6 \text{ mas} \sim 30 \text{ AU}$ spherical wind component to the Br$\gamma$ line: the entire binary is engulfed in an optically thin spherical line emission envelope.

The broadening of atomic emission lines by high-velocity motion of gas near accreting supermassive black holes is an observational hallmark of quasars. Observations of broad emission lines could potentially constrain the mechanism for transporting gas inwards through accretion disks or outwards through winds. The size of this broad-line region has been estimated by measuring the light travel time delay between the variable nuclear continuum and the emission lines - a method known as reverberation mapping. In some models the emission lines arise from a continuous outflow, whereas in others they are produced by orbiting gas clouds. Directly imaging such regions has not hitherto been possible because of their small angular sizes (< 0.1 milli-arcseconds). Here we report a spatial offset (with a spatial resolution of ten micro-arcseconds or about 0.03 parsecs for a distance of 550 million parsecs) between the red and blue photo-centres of the broad Paschen-{\alpha} line of the quasar 3C 273 perpendicular to the direction of its radio jet. This spatial offset corresponds to a gradient in the velocity of the gas and thus implies that the gas is orbiting the central supermassive black hole. The data are well fitted by a broad-line-region model of a thick disk of gravitationally bound material orbiting a black hole of 300 million solar masses. We infer a disk radius of 150 light days; a radius of 100-400 light days was found previously using reverberation mapping. The rotation axis of the disk aligns in inclination and position angle with the radio jet. Our results support the methods that are often used to estimate the masses of accreting supermassive black holes and to study their evolution over cosmic time.

We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.

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.

Pulsars in the Galactic Center (GC) are important probes of General Relativity, star formation, stellar dynamics, stellar evolution, and the interstellar medium. Despite years of searching, only a handful of pulsars in the central 0.5 deg are known. The high-frequency sensitivity of ngVLA will open a new window for discovery and characterization of pulsars in the GC. A pulsar in orbit around the GC black hole, Sgr A*, will provide an unprecedented probe of black hole physics and General Relativity.

This work presents an interferometric study of the massive-binary fraction in the Orion Trapezium Cluster with the recently comissioned GRAVITY instrument. We observe a total of 16 stars of mainly OB spectral type. We find three previously unknown companions for $\theta ^1$ Ori B, $\theta ^2$ Ori B, and $\theta ^2$ Ori C. We determine a separation for the previously suspected companion of NU Ori. We confirm four companions for $\theta ^1$ Ori A, $\theta ^1$ Ori C, $\theta ^1$ Ori D, and $\theta ^2$ Ori A, all with substantially improved astrometry and photometric mass estimates. We refine the orbit of the eccentric high-mass binary $\theta ^1$ Ori C and we are able to derive a new orbit for $\theta ^1$ Ori D. We find a system mass of 21.7 $M_{\odot}$ and a period of $53$ days. Together with other previously detected companions seen in spectroscopy or direct imaging, eleven of the 16 high-mass stars are multiple systems. We obtain a total number of 22 companions with separations up to 600 AU. The companion fraction of the early B and O stars in our sample is about 2, significantly higher than in earlier studies of mostly OB associations. The separation distribution hints towards a bimodality. Such a bimodality has been previously found in A stars, but rarely in OB binaries, which up to this point have been assumed to be mostly compact with a tail of wider companions. We also do not find a substantial population of equal-mass binaries. The observed distribution of mass ratios declines steeply with mass, and like the direct star counts, indicates that our companions follow a standard power law initial mass function. Again, this is in contrast to earlier findings of flat mass ratio distributions in OB associations. We exclude collision as a dominant formation mechanism but find no clear preference for core accretion or competitive accretion.

Radio images of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), are dominated by interstellar scattering. Previous studies of Sgr A* have adopted an anisotropic Gaussian model for both the intrinsic source and the scattering, and they have extrapolated the scattering using a purely $\lambda^2$ scaling to estimate intrinsic properties. However, physically motivated source and scattering models break all three of these assumptions. They also predict that refractive scattering effects will be significant, which have been ignored in standard model fitting procedures. We analyze radio observations of Sgr A* using a physically motivated scattering model, and we develop a prescription to incorporate refractive scattering uncertainties when model fitting. We show that an anisotropic Gaussian scattering kernel is an excellent approximation for Sgr A* at wavelengths longer than 1cm, with an angular size of $(1.380 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the major axis, $(0.703 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the minor axis, and a position angle of $81.9^\circ \pm 0.2^\circ$. We estimate that the turbulent dissipation scale is at least $600\,{\rm km}$, with tentative support for $r_{\rm in} = 800 \pm 200\,{\rm km}$, suggesting that the ion Larmor radius defines the dissipation scale. We find that the power-law index for density fluctuations in the scattering material is $\beta < 3.47$, shallower than expected for a Kolmogorov spectrum ($\beta=11/3$), and we estimate $\beta = 3.38^{+0.08}_{-0.04}$ in the case of $r_{\rm in} = 800\,{\rm km}$. We find that the intrinsic structure of Sgr A* is nearly isotropic over wavelengths from 1.3mm to 1.3cm, with a size that is roughly proportional to wavelength. We discuss implications for models of Sgr A*, for theories of interstellar turbulence, and for imaging Sgr A* with the Event Horizon Telescope.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.

Rapid, large amplitude variability at optical to X-ray wavelengths is now seen in an increasing number of Seyfert galaxies and luminous quasars. The variations imply a global change in accretion power, but are too rapid to be communicated by inflow through a standard thin accretion disc. Such discs are long known to have difficulty explaining the observed optical/UV emission from active galactic nuclei. Here we show that alternative models developed to explain these observations have larger scale heights and shorter inflow times. Accretion discs supported by magnetic pressure in particular are geometrically thick at all luminosities, with inflow times as short as the observed few year timescales in extreme variability events to date. Future time-resolved, multi-wavelength observations can distinguish between inflow through a geometrically thick disc as proposed here, and alternative scenarios of extreme reprocessing of a central source or instability-driven limit cycles.

We report the first detection of the Galactic Centre massive black hole, Sgr~A*, in the far infrared. Our measurements were obtained with PACS on board the \emph{Herschel} satellite at $100~\mathrm{\mu m}$ and $160~\mathrm{\mu m}$. While the warm dust in the Galactic Centre is too bright to allow for a direct detection of Sgr~A*, we measure a significant and simultaneous variation of its flux of $\Delta F_{\nu\widehat{=}160 ~\mathrm{\mu m}} = (0.27\pm0.06)~\mathrm{Jy}$ and $\Delta F_{\nu\widehat{=}100 ~\mathrm{\mu m}}= (0.16\pm0.10)~\mathrm{Jy}$ during one observation. The significance level of the $160 ~\mathrm{\mu m}$ band variability is $4.5\sigma$ and the corresponding $100 ~\mathrm{\mu m}$ band variability is significant at $1.6\sigma$. We find no example of an equally significant false positive detection. Conservatively assuming a variability of $25\%$ in the FIR, we can provide upper limits to the flux. Comparing the latter with theoretical models we find that 1D RIAF models have difficulties explaining the observed faintness. However, the upper limits are consistent with modern ALMA and VLA observations. Our upper limits provide further evidence for a spectral peak at $\sim 10^{12} ~ \mathrm{Hz}$ and constrain the number density of $\gamma \sim 100$ electrons in the accretion disk and or outflow.

We report results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional {\it uv} coverage in the N-S direction, and leads to a spatial resolution of $\sim$30 $\mu$as ($\sim$3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of $\sim$4-13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of $\sim$3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow a more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

We study model images and polarization maps of Sagittarius A* at 230 GHz. We post-process GRMHD simulations and perform a fully relativistic radiative transfer calculation of the emitted synchrotron radiation to obtain polarized images for a range of mass accretion rates and electron temperatures. At low accretion rates, the polarization map traces the underlying toroidal magnetic field geometry. At high accretion rates, we find that Faraday rotation internal to the emission region can depolarize and scramble the map. We measure the net linear polarization fraction and find that high accretion rate "jet-disc" models are heavily depolarized and are therefore disfavoured. We show how Event Horizon Telescope measurements of the polarized "correlation length" over the image provide a model-independent upper limit on the strength of these Faraday effects, and constrain plasma properties like the electron temperature and magnetic field strength.

In this research note, we present a set of distortion solutions that may be used to correct geometric optical distortion in images taken with the S13 camera of the NACO adaptive optics imager.

Astrometric and spectroscopic monitoring of individual stars orbiting the supermassive black hole in the Galactic Center offer a promising way to detect general relativistic effects. While low-order effects are expected to be detected following the periastron passage of S2 in Spring 2018, detecting higher-order effects due to black hole spin will require the discovery of closer stars. In this paper, we set out to determine the requirements such a star would have to satisfy to allow the detection of black hole spin. We focus on the instrument GRAVITY, which saw first light in 2016 and which is expected to achieve astrometric accuracies $10-100 \mu$as. For an observing campaign with duration $T$ years, $N_{obs}$ total observations, astrometric precision $\sigma_x$ and normalized black hole spin $\chi$, we find that $a_{orb}(1-e^2)^{3/4} \lesssim 300 R_S \sqrt{\frac{T}{4 \text{years}}} \left(\frac{N_{obs}}{120}\right)^{0.25} \sqrt{\frac{10 \mu as}{\sigma_x}} \sqrt{\frac{\chi}{0.9}}$ is needed. For $\chi=0.9$ and a potential observing campaign with $\sigma_x = 10 \mu$as, 30 observations/year and duration 4-10 years, we expect $\sim 0.1$ star with $K<19$ satisfying this constraint based on the current knowledge about the stellar population in the central 1". We also propose a method through which GRAVITY could potentially measure radial velocities with precision $\sim 50$ km/s. If the astrometric precision can be maintained, adding radial velocity information increases the expected number of stars by roughly a factor of two. While we focus on GRAVITY, the results can also be scaled to parameters relevant for future extremely large telescopes.

Millimeter wavelength polarimetry of accreting black hole systems can provide a tomographic probe of the accretion flow on a wide range of linear scales. We searched for linear polarization in two low luminosity active galactic nuclei (LLAGN), M81 and M84, using the Combined Array for Millimeter Astronomy (CARMA) and the Submillimeter Array (SMA). We find upper limits of $\sim 1 - 2\%$ averaging over the full bandwidth and with a rotation measure (RM) synthesis technique. These low polarization fractions, along with similar low values for LLAGN M87 and 3C84, suggest that LLAGN have qualitatively different polarization properties than radio-loud sources and Sgr A*. If the sources are intrinsically polarized and then depolarized by Faraday rotation then we place lower limits on the RM of a few times $10^7\, {\rm rad\, m^{-2}}$ for the full bandwidth case and $\sim 10^9\, {\rm rad\, m^{-2}}$ for the RM synthesis analysis. These limits are inconsistent with or marginally consistent with expected accretion flow properties. Alternatively, the sources may be depolarized by cold electrons within a few Schwarzschild radii from the black hole, as suggested by numerical models.

After changing optical AGN type from 1.9 to 1 in 1984, the AGN Mrk 1018 recently reverted back to its type 1.9 state. Our ongoing monitoring now reveals that the AGN has halted its dramatic dimming, reaching a minimum around October 2016. The minimum was followed by an outburst rising with $\sim$0.25 U-band mag/month. The rebrightening lasted at least until February 2017, as confirmed by joint Chandra and Hubble observations. Monitoring was resumed in July 2017 after the source emerged from sunblock, at which point the AGN was found only $\sim$0.4 mag brighter than its minimum. The intermittent outburst was accompanied by the appearance of a red wing asymmetry in broad-line shape, indicative of an inhomogeneous broad-line region. The current flickering brightness of Mrk 1018 following its rapid fading either suggests that the source has reignited, remains variable at a low level, or may continue dimming over the next few years. Distinguishing between these possibilities requires continuous multiwavelength monitoring.

The G2 object has recently passed its pericenter passage in our Galactic Center. While the $Br_\gamma$ emission shows clear signs of tidal interaction, the change in the observed luminosity is only of about a factor of 2, in contention with all previous predictions. We present high resolution simulations performed with the moving mesh code, RICH, together with simple analytical arguments that reproduce the observed $Br_\gamma$ emission. In our model, G2 is a gas cloud that undergoes tidal disruption in a dilute ambient medium. We find that during pericenter passage, the efficient cooling of the cloud results in a vertical collapse, compressing the cloud by a factor of $\sim5000$. By properly taking into account the ionization state of the gas, we find that the cloud is UV starved and are able to reproduce the observed $Br_\gamma$ luminosity. For densities larger than $\approx500\;\mathrm{cm}^{-3}$ at pericenter, the cloud fragments, due to cooling instabilities and the emitted radiation is inconsistent with observations. For lower densities, the cloud survives the pericenter passage intact and its emitted radiation matches the observed lightcurve. From the duration of $Br_\gamma$ emission which contains both redshifted and blueshifted components, we show that the cloud is not spherical but rather elongated with a size ratio of 4 at year 2001. The simulated cloud's elongation grows as it travels towards pericenter and is consistent with observations, due to viewing angles. The simulation is also consistent with having a spherical shape at apocenter.

To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Our observations in the near-infrared K-band continuum reveal a disk around each binary component, with similar half-flux radii of about 0.1 au at d$\sim$130 pc, inclinations ($i=$28$\pm$3$^o$\ and $i=$22$\pm$6$^o$), and position angles (PA=0$^o\pm$6$^o$ and PA=-2$^o\pm$12$^o$), suggesting that they formed from the fragmentation of a common disk. The S CrA N spectrum shows bright HeI and Br$\gamma$ line emission exhibiting inverse P-Cygni profiles, typically associated with infalling gas. The continuum-compensated Br$\gamma$ line visibilities of S CrA N show the presence of a compact Br$\gamma$ emitting region the radius of which is about $\sim$0.06 au, which is twice as big as the truncation radius. This component is mostly tracing a wind. Moreover, a slight radius change between the blue- and red-shifted Br$\gamma$ line components is marginally detected. The presence of an inverse P-Cygni profile in the HeI and Br$\gamma$ lines, along with the tentative detection of a slightly larger size of the blue-shifted Br$\gamma$ line component, hint at the simultaneous presence of a wind and magnetospheric accretion in S CrA N.

We study the young S-stars within a distance of 0.04 pc from the supermassive black hole in the center of our Galaxy. Given how inhospitable the region is for star formation, their presence is more puzzling the younger we estimate their ages. In this study, we analyse the result of 12 years of high resolution spectroscopy within the central arcsecond of the Galactic Center (GC). By co-adding between 55 and 105 hours of spectra we have obtained high signal to noise H- and K-band spectra of eight stars orbiting the central supermassive black hole. Using deep H-band spectra, we show that these stars must be high surface gravity (dwarf) stars. We compare these deep spectra to detailed model atmospheres and stellar evolution models to infer the stellar parameters. Our analysis reveals an effective temperature of 21000-28500 K, a rotational velocity of 60-170 km/s, and a surface gravity of 4.1-4.2. These parameters imply a spectral type of B0-B3V for these stars. The inferred masses lie within 8-14 Msun. We derive an age of 6.6^{+3.4}{-4.7} Myr for the star S2, which is compatible with the age of the clockwise rotating young stellar disk in the GC. We estimate the age of all other studied S-stars to be less than 15 Myr, which are compatible with the age of S2 within the uncertainties. The relatively low ages for these S-stars favor a scenario in which the stars formed in a local disk rather than the field-binary-disruption scenario throughout a longer period of time.

We use VLBA+VLA observations to measure the sizes of the scatter-broadened images of 6 of the most heavily scattered known pulsars: 3 within the Galactic Centre (GC) and 3 elsewhere in the inner Galactic plane. By combining the measured sizes with temporal pulse broadening data from the literature and using the thin-screen approximation, we locate the scattering medium along the line of sight to these 6 pulsars. At least two scattering screens are needed to explain the observations of the GC sample. We show that the screen inferred by previous observations of SGR J1745-2900 and Sgr A*, which must be located far from the GC, falls off in strength on scales < 0.2 degree. A second scattering component closer to (< 2 kpc) or even (tentatively) within (< 700 pc) the GC produces most or all of the temporal broadening observed in the other GC pulsars. Outside the GC, the scattering locations for all three pulsars are ~2 kpc from Earth, consistent with the distance of the Carina-Sagittarius or Scutum spiral arm. For each object the 3D scattering origin coincides with a known HII region (and in one case also a supernova remnant), suggesting that such objects preferentially cause the intense interstellar scattering seen towards the Galactic plane. We show that the HII regions should contribute > 25% of the total dispersion measure (DM) towards these pulsars, and calculate reduced DM distances. Those distances for other pulsars lying behind HII regions may be similarly overestimated.

Monitoring of Sagittarius A* from X-ray to radio wavelengths has revealed structured variability --- including X-ray flares --- but it is challenging to establish correlations between them. Most studies have focused on variability in the X-ray and infrared, where variations are often simultaneous, and because long time series at sub-millimeter and radio wavelengths are limited. Previous work on sub-mm and radio variability hints at a lag between X-ray flares and their candidate sub-millimeter or radio counterparts, with the long wavelength data lagging the X-ray. However, there is only one published time lag between an X-ray flare and a possible radio counterpart. Here we report 9 contemporaneous X-ray and radio observations of Sgr A*. We detect significant radio variability peaking $\gtrsim$176 minutes after the brightest X-ray flare ever detected from Sgr A*. We also report other potentially associated X-ray and radio variability, with the radio peaks appearing $\lesssim$80 minutes after these weaker X-ray flares. Taken at face value, these results suggest that stronger X-ray flares lead to longer time lags in the radio. However, we also test the possibility that the variability at X-ray and radio wavelengths is not temporally correlated. We cross-correlate data from mismatched X-ray and radio epochs and obtain comparable correlations to the matched data. Hence, we find no overall statistical evidence that X-ray flares and radio variability are correlated, underscoring a need for more simultaneous, long duration X-ray--radio monitoring of Sgr A*.

Understanding the origin of the flaring activity from the Galactic center supermassive black hole, Sagittarius A*, is a major scientific goal of the NuSTAR Galactic plane survey campaign. We report on the data obtained between July 2012 and April 2015, including 27 observations on Sgr A* with a total exposure of ~ 1 Ms. We found a total of ten X-ray flares detected in the NuSTAR observation window, with luminosities in the range of $L_{3-79~keV}$~$(0.2-4.0) \times 10^{35}~erg~s^{-1}$. With this largest hard X-ray Sgr A* flare dataset to date, we studied the flare spectral properties. Seven flares are detected above 5{\sigma} significance, showing a range of photon indices ({\Gamma} ~ 2.0-2.8) with typical uncertainties of +/-0.5 (90% confidence level). We found no significant spectral hardening for brighter flares as indicated by a smaller sample. The accumulation of all the flare spectra in 1-79 keV can be well fit with an absorbed power-law model with {\Gamma}=2.2+/-0.1, and does not require the existence of a spectral break. The lack of variation in X-ray spectral index with luminosity would point to a single mechanism for the flares and is consistent with the synchrotron scenario. Lastly, we present the quiescent state spectrum of Sgr A*, and derived an upper limit on the quiescent luminosity of Sgr A* above 10 keV to be $L_{Xq, 10-79 keV}$ < $(2.9{\pm}0.2) \times 10^{34}~erg~s^{-1}$.

GRAVITY is a new instrument to coherently combine the light of the European Southern Observatory Very Large Telescope Interferometer to form a telescope with an equivalent 130 m diameter angular resolution and a collecting area of 200 m$^2$. The instrument comprises fiber fed integrated optics beam combination, high resolution spectroscopy, built-in beam analysis and control, near-infrared wavefront sensing, phase-tracking, dual beam operation and laser metrology [...]. This article gives an overview of GRAVITY and reports on the performance and the first astronomical observations during commissioning in 2015/16. We demonstrate phase tracking on stars as faint as m$_K$ ~ 10 mag, phase-referenced interferometry of objects fainter than m$_K$ ~ 15 mag with a limiting magnitude of m$_K$ ~ 17 mag, minute long coherent integrations, a visibility accuracy of better than 0.25 %, and spectro-differential phase and closure phase accuracy better than 0.5{\deg}, corresponding to a differential astrometric precision of better than 10 microarcseconds ({\mu}as). The dual-beam astrometry, measuring the phase difference of two objects with laser metrology, is still under commissioning. First observations show residuals as low as 50 {\mu}as when following objects over several months. We illustrate the instrument performance with the observations of archetypical objects for the different instrument modes. Examples include the Galactic Center supermassive black hole and its fast orbiting star S2 for phase referenced dual beam observations and infrared wavefront sensing, the High Mass X-Ray Binary BP Cru and the Active Galactic Nucleus of PDS 456 for few {\mu}as spectro-differential astrometry, the T Tauri star S CrA for a spectro-differential visibility analysis, {\xi} Tel and 24 Cap for high accuracy visibility observations, and {\eta} Car for interferometric imaging with GRAVITY.

We present the first optical observation at sub-milliarcsecond (mas) scale of the microquasar SS 433 obtained with the GRAVITY instrument on the VLT interferometer. The 3.5 hour exposure reveals a rich K-band spectrum dominated by hydrogen Br$\gamma\ $ and \ion{He}{i} lines, as well as (red-shifted) emission lines coming from the jets. The K-band continuum emitting region is dominated by a marginally resolved point source ($<$ 1 mas) embedded inside a diffuse background accounting for 10\% of the total flux. The jet line positions agree well with the ones expected from the jet kinematic model, an interpretation also supported by the consistent sign (i.e. negative/positive for the receding/approaching jet component) of the phase shifts observed in the lines. The significant visibility drop across the jet lines, together with the small and nearly identical phases for all baselines, point toward a jet that is offset by less than 0.5 mas from the continuum source and resolved in the direction of propagation, with a typical size of 2 mas. The jet position angle of $\sim$80$^{\circ}$ is consistent with the expected one at the observation date. Jet emission so close to the central binary system would suggest that line locking, if relevant to explain the amplitude and stability of the 0.26c jet velocity, operates on elements heavier than hydrogen. The Br$\gamma\ $ profile is broad and double peaked. It is better resolved than the continuum and the change of the phase signal sign across the line on all baselines suggests an East-West oriented geometry alike the jet direction and supporting a (polar) disk wind origin.

We observe the HMXB BP Cru using interferometry in the near-infrared K band with VLTI/GRAVITY. Continuum visibilities are at most partially resolved, consistent with the predicted size of the hypergiant. Differential visibility amplitude ($\Delta |V| \sim 5\%$) and phase ($\Delta \phi \sim 2 \deg$) signatures are observed across the HeI $2.059 \mu$m and Br$\gamma$ lines, the latter seen strongly in emission, unusual for the donor star's spectral type. For a baseline $B \sim 100$m, the differential phase RMS $\sim 0.2 \deg$ corresponds to an astrometric precision of $\sim 2 \mu$as. A model-independent analysis in the marginally resolved limit of interferometry reveals asymmetric and extended emission with a strong wavelength dependence. We propose geometric models based on an extended and distorted wind and/or a high density gas stream, which has long been predicted to be present in this system. The observations show that optical interferometry is now able to resolve HMXBs at the spatial scale at which accretion takes place, and therefore probe the effects of the gravitational and radiation fields of the compact object on its environment.

In early 2014 the fast-moving near-infrared source G2 reached its closest approach to the supermassive black hole Sgr A* in the Galactic Center. We report on the evolution of the ionized gaseous component and the dusty component of G2 immediately after this event, revealed by new observations obtained in 2015 and 2016 with the SINFONI integral field spectrograph and the NACO imager at the ESO VLT. The spatially resolved dynamics of the Br$\gamma$ line emission can be accounted for by the ballistic motion and tidal shearing of a test-particle cloud that has followed a highly eccentric Keplerian orbit around the black hole for the last 12 years. The non-detection of a drag force or any strong hydrodynamic interaction with the hot gas in the inner accretion zone limits the ambient density to less than a few 10$^3$ cm$^{-3}$ at the distance of closest approach (1500 $R_s$), assuming G2 is a spherical cloud moving through a stationary and homogeneous atmosphere. The dust continuum emission is unresolved in L'-band, but stays consistent with the location of the Br$\gamma$ emission. The total luminosity of the Br$\gamma$ and L' emission has remained constant to within the measurement uncertainty. The nature and origin of G2 are likely related to that of the precursor source G1, since their orbital evolution is similar, though not identical. Both object are also likely related to a trailing tail structure, which is continuously connected to G2 over a large range in position and radial velocity.

We present the first fully simultaneous fits to the NIR and X-ray spectral slope (and its evolution) during a very bright flare from Sgr A*, the supermassive black hole at the Milky Way's center. Our study arises from ambitious multi-wavelength monitoring campaigns with XMM-Newton, NuSTAR and SINFONI. The average multi-wavelength spectrum is well reproduced by a broken power-law with $\Gamma_{NIR}=1.7\pm0.1$ and $\Gamma_X=2.27\pm0.12$. The difference in spectral slopes ($\Delta\Gamma=0.57\pm0.09$) strongly supports synchrotron emission with a cooling break. The flare starts first in the NIR with a flat and bright NIR spectrum, while X-ray radiation is detected only after about $10^3$ s, when a very steep X-ray spectrum ($\Delta\Gamma=1.8\pm0.4$) is observed. These measurements are consistent with synchrotron emission with a cooling break and they suggest that the high energy cut-off in the electron distribution ($\gamma_{max}$) induces an initial cut-off in the optical-UV band that evolves slowly into the X-ray band. The temporal and spectral evolution observed in all bright X-ray flares are also in line with a slow evolution of $\gamma_{max}$. We also observe hints for a variation of the cooling break that might be induced by an evolution of the magnetic field (from $B\sim30\pm8$ G to $B\sim4.8\pm1.7$ G at the X-ray peak). Such drop of the magnetic field at the flare peak would be expected if the acceleration mechanism is tapping energy from the magnetic field, such as in magnetic reconnection. We conclude that synchrotron emission with a cooling break is a viable process for Sgr A*'s flaring emission.

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.

Recent studies have shown possible connections between highly magnetized neutron stars ("magnetars"), whose X-ray emission is too bright to be powered by rotational energy, and ordinary radio pulsars. In addition to the magnetar SGR J1745-2900, one of the radio pulsars in the Galactic centre (GC) region, PSR J1746-2850, had timing properties implying a large magnetic field strength and young age, as well as a flat spectrum. All characteristics are similar to those of rare, transient, radio-loud magnetars. Using several deep non-detections from the literature and two new detections, we show that this pulsar is also transient in the radio. Both the flat spectrum and large amplitude variability are inconsistent with the light curves and spectral indices of 3 radio pulsars with high magnetic field strengths. We further use frequent, deep archival imaging observations of the GC in the past 15 years to rule out a possible X-ray outburst with a luminosity exceeding the rotational spin down rate. This source, either a transient magnetar without any detected X-ray counterpart or a young, strongly magnetized radio pulsar producing magnetar-like radio emission, further blurs the line between the two categories. We discuss the implications of this object for the radio emission mechanism in magnetars and for star and compact object formation in the GC.

We propose a new technique for radio interferometry to obtain super-resolution full polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the $\ell _1$-norm and total variation (TV) of the brightness distribution. As an application of this technique, we present simulated linear polarization observations of two physically motivated models of M87 with the Event Horizon Telescope (EHT). We confirm that $\ell _1$+TV regularization can achieve an optimal resolution of $\sim 25-30$\% of the diffraction limit $\lambda/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer for both the total intensity (i.e. Stokes $I$) and linear polarizations (i.e. Stokes $Q$ and $U$). This optimal resolution is better than that obtained from the widely used Cotton-Schwab CLEAN algorithm or from using $\ell _1$ or TV regularizations alone. Furthermore, we find that $\ell _1$+TV regularization can achieve much better image fidelity in linear polarization than other techniques over a wide range of spatial scales, not only in the super-resolution regime, but also on scales larger than the diffraction limit. Our results clearly demonstrate that sparse reconstruction is a useful choice for high-fidelity full-polarimetric interferometric imaging.

We propose a new imaging technique for radio and optical/infrared interferometry. The proposed technique reconstructs the image from the visibility amplitude and closure phase, which are standard data products of short-millimeter very long baseline interferometers such as the Event Horizon Telescope (EHT) and optical/infrared interferometers, by utilizing two regularization functions: the $\ell_1$-norm and total variation (TV) of the brightness distribution. In the proposed method, optimal regularization parameters, which represent the sparseness and effective spatial resolution of the image, are derived from data themselves using cross validation (CV). As an application of this technique, we present simulated observations of M87 with the EHT based on four physically motivated models. We confirm that $\ell_1$+TV regularization can achieve an optimal resolution of $\sim 20-30$% of the diffraction limit $\lambda/D_{\rm max}$, which is the nominal spatial resolution of a radio interferometer. With the proposed technique, the EHT can robustly and reasonably achieve super-resolution sufficient to clearly resolve the black hole shadow. These results make it promising for the EHT to provide an unprecedented view of the event-horizon-scale structure in the vicinity of the super-massive black hole in M87 and also the Galactic center Sgr A*.

Until recently, the strong interstellar scattering observed towards the Galactic center (GC) black hole, Sgr A*, was thought to come from dense gas within the GC region. The pulse broadening towards the transient magnetar SGR J1745-2900 near Sgr A* has shown that the source of the scattering is instead located much closer to Earth, possibly in a nearby spiral arm. We show that a single HII region along the line of sight, 1.5-4.8 kpc away from Earth with density $n_e$ of a few 100 cm$^{-3}$ and radius R = 1.8-3.2 pc can explain the observed angular broadening of Sgr A*. Clouds closer to the GC overproduce the observed DM, providing an independent location constraint that agrees with that from the magnetar pulse broadening. Our model predicts that sources within 10 pc should show the same scattering origin as the magnetar and Sgr A*, while the nearest known pulsars with separations > 20 pc should not. The radio spectrum of Sgr A* should show a cutoff from free-free absorption at 0.2-1 GHz. For a magnetic field strength $B = 15-70 \mu$G, the HII region could produce the rotation measure of the magnetar, the largest of any known pulsar, without requiring the gas near Sgr A* to be strongly magnetised.

Using 25 years of data from uninterrupted monitoring of stellar orbits in the Galactic Center, we present an update of the main results from this unique data set: A measurement of mass of and distance to SgrA*. Our progress is not only due to the eight year increase in time base, but also due to the improved definition of the coordinate system. The star S2 continues to yield the best constraints on the mass of and distance to SgrA*; the statistical errors of 0.13 x 10^6 M_sun and 0.12 kpc have halved compared to the previous study. The S2 orbit fit is robust and does not need any prior information. Using coordinate system priors, also the star S1 yields tight constraints on mass and distance. For a combined orbit fit, we use 17 stars, which yields our current best estimates for mass and distance: M = 4.28 +/- 0.10|stat. +/. 0.21|sys. x 10^6 M_sun and R_0 = 8.32 +/- 0.07|stat. +/- 0.14|sys. kpc. These numbers are in agreement with the recent determination of R_0 from the statistical cluster parallax. The positions of the mass, of the near-infrared flares from SgrA* and of the radio source SgrA* agree to within 1mas. In total, we have determined orbits for 40 stars so far, a sample which consists of 32 stars with randomly oriented orbits and a thermal eccentricity distribution, plus eight stars for which we can explicitly show that they are members of the clockwise disk of young stars, and which have lower eccentricity orbits.

Decomposing an arbitrary electron energy distribution into sums of Maxwellian and power law components is an efficient method to calculate synchrotron emission and absorption. We use this method to study the effect of non-thermal electrons on submm images and spectra of the Galactic center black hole, Sgr A*. We assume a spatially uniform functional form for the electron distribution function and use a semi-analytic radiatively inefficient accretion flow and a 2D general relativistic MHD snapshot as example models of the underlying accretion flow structure. We develop simple analytic models which allow us to generalize from the numerical examples. A high energy electron component containing a small fraction (few per cent) of the total internal energy (e.g. a "power law tail") can produce a diffuse halo of emission, which modifies the observed image size and structure. A population of hot electrons with a larger energy fraction (e.g. resulting from a diffusion in electron energy space) can dominate the emission, so that the observed images and spectra are well approximated by considering only a single thermal component for a suitable choice of the electron temperature. We discuss the implications of these results for estimating accretion flow or black hole parameters from images and spectra, and for the identification of the black hole "shadow" in future mm-VLBI data. In particular, the location of the first minimum in visibility profiles does not necessarily correspond to the shadow size as sometimes assumed.

We recently discovered that the active galactic nucleus (AGN) of Mrk 1018 has changed optical type again after 30 years as a type 1 AGN. Here we combine Chandra, NuStar, Swift, Hubble Space Telescope and ground-based observations to explore the cause of this change. The 2-10keV flux declines by a factor of ~8 between 2010 and 2016. We show with our X-ray observation that this is not caused by varying neutral hydrogen absorption along the line-of-sight up to the Compton-thick level. The optical-UV spectral energy distributions are well fit with a standard geometrically thin optically thick accretion disc model that seems to obey the expected $L\sim T^4$ relation. It confirms that a decline in accretion disc luminosity is the primary origin for the type change. We detect a new narrow-line absorber in Lya blue-shifted by ~700km/s with respect to the systemic velocity of the galaxy. This new Lya absorber could be evidence for the onset of an outflow or a companion black hole with associated gas that could be related to the accretion rate change. However, the low column density of the absorber means that it is not the direct cause for Mrk 1018's changing-look nature.

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

Ray tracing radiative transfer is a powerful method for comparing theoretical models of black hole accretion flows and jets with observations. We present a public code, grtrans, for carrying out such calculations in the Kerr metric, including the full treatment of polarised radiative transfer and parallel transport along geodesics. The code is written in Fortran 90 and efficiently parallelises with OpenMP, and the full code and several components have Python interfaces. We describe several tests which are used for verifiying the code, and we compare the results for polarised thin accretion disc and semi-analytic jet problems with those from the literature as examples of its use. Along the way, we provide accurate fitting functions for polarised synchrotron emission and transfer coefficients from thermal and power law distribution functions, and compare results from numerical integration and quadrature solutions of the polarised radiative transfer equations. We also show that all transfer coefficients can play an important role in predicted images and polarisation maps of the Galactic center black hole, Sgr A*, at submillimetre wavelengths.

Gamma-ray observations have shown pulsars to be efficient converters of rotational energy into GeV photons and it is of wide-ranging interest to determine their contribution to the gamma-ray background. We arrive at flux predictions from both the young (<~ Myr) and millisecond (~Gyr) Galactic pulsar populations. We find that unresolved pulsars can yield both a significant fraction of the excess GeV gamma rays near the Galactic Center and an inverse Compton flux in the inner kpc similar to that inferred by Fermi. We compare models of the young pulsar population and millisecond pulsar population to constraints from gamma-ray and radio observations. Overall, we find that the young pulsars should outnumber millisecond pulsars as unassociated gamma-ray point sources in this region. The number of young radio pulsars discovered near the Galactic Center is in agreement with our model of the young pulsar population. Deeper radio observations at higher latitudes can constrain the total gamma-ray emission from both young and millisecond pulsars from the inner galaxy. While this is a step towards better understanding of pulsars, cosmic rays in the Milky Way, and searches for dark matter, we also discuss a few interesting puzzles that arise from the underlying physics of pulsar emission and evolution.

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

Near-infrared observations of stellar orbits at the Galactic Center provide conclusive evidence for a massive black hole associated with the compact radio source Sgr A*. The astrometric reference frame for these observations is tied to a set of red giant stars, which are also detectable at radio wavelengths through SiO maser emission in their envelopes. We have improved the precision and long-term stability of this reference frame, in which Sgr A* is localized to within a factor 5 better than previously: ~0.17 mas in position (in 2009) and ~0.07 mas/yr in velocity. This improvement is the result of modeling and correcting optical distortion in the VLT/NACO imager to a sub-mas level and including new infrared and radio measurements, which now both span more than a decade in time. A further improvement will follow future observations and facilitate the detection of relativistic orbital effects.

We analyze the light curves of 413 radio sources at submillimeter wavelengths using data from the Submillimeter Array calibrator database. The database includes more than 20,000 observations at 1.3 and 0.8 mm that span 13 years. We model the light curves as a damped random walk and determine a characteristic time scale $\tau$ at which the variability amplitude saturates. For the vast majority of sources, primarily blazars and BL Lac objects, we find only lower limits on $\tau$. For two nearby low luminosity active galactic nuclei, M81 and M87, however, we measure $\tau=1.6^{+3.0}_{-0.9}$ days and $\tau=45^{+61}_{-24}$ days, respectively ($2\sigma$ errors). Including the previously measured $\tau=0.33\pm 0.16$ days for Sgr A*, we show an approximately linear correlation between $\tau$ and black hole mass for these nearby LLAGN. Other LLAGN with spectra that peak in the submm are expected to follow this correlation. These characteristic time scales are comparable to the minimum time scale for emission processes close to an event horizon, and suggest that the underlying physics may be independent of black hole mass, accretion rate, and jet luminosity.

X-ray flares have routinely been observed from the supermassive black hole, Sagittarius A$^\star$ (Sgr A$^\star$), at our Galactic center. The nature of these flares remains largely unclear, despite of many theoretical models. In this paper, we study the statistical properties of the Sgr A$^\star$ X-ray flares, by fitting the count rate (CR) distribution and the structure function (SF) of the light curve with a Markov Chain Monte Carlo (MCMC) method. With the 3 million second \textit{Chandra} observations accumulated in the Sgr A$^\star$ X-ray Visionary Project, we construct the theoretical light curves through Monte Carlo simulations. We find that the $2-8$ keV X-ray light curve can be decomposed into a quiescent component with a constant count rate of $\sim6\times10^{-3}~$count s$^{-1}$ and a flare component with a power-law fluence distribution $dN/dE\propto E^{-\alpha_{\rm E}}$ with $\alpha_{\rm E}=1.65\pm0.17$. The duration-fluence correlation can also be modelled as a power-law $T\propto E^{\alpha_{\rm ET}}$ with $\alpha_{\rm ET} < 0.55$ ($95\%$ confidence). These statistical properties are consistent with the theoretical prediction of the self-organized criticality (SOC) system with the spatial dimension $S = 3$. We suggest that the X-ray flares represent plasmoid ejections driven by magnetic reconnection (similar to solar flares) in the accretion flow onto the black hole.

We report on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm \mu$as $\sim$ 5.5 $R_{{\rm s}}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{{\rm s}}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm \gamma$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{{\rm s}}$.

Studies of Fermi data indicate an excess of GeV gamma rays around the Galactic center (GC), possibly due to dark matter. We show that young gamma-ray pulsars can yield a similar signal. First, a high concentration of GC supernovae naturally leads to a population of kicked pulsars symmetric about the GC. Second, while very-young pulsars with soft spectra reside near the Galactic plane, pulsars with spectra that have hardened with age accumulate at larger angles. This combination, including unresolved foreground pulsars, traces the morphology and spectrum of the Excess.

We report new observations with the Very Large Array, Atacama Large Millimeter Array, and Submillimeter Array at frequencies from 1.0 to 355 GHz of the Galactic Center black hole, Sagittarius A*. These observations were conducted between October 2012 and November 2014. While we see variability over the whole spectrum with an amplitude as large as a factor of 2 at millimeter wavelengths, we find no evidence for a change in the mean flux density or spectrum of Sgr A* that can be attributed to interaction with the G2 source. The absence of a bow shock at low frequencies is consistent with a cross-sectional area for G2 that is less than $2 \times 10^{29}$ cm$^2$. This result fits with several model predictions including a magnetically arrested cloud, a pressure-confined stellar wind, and a stellar photosphere of a binary merger. There is no evidence for enhanced accretion onto the black hole driving greater jet and/or accretion flow emission. Finally, we measure the millimeter wavelength spectral index of Sgr A* to be flat; combined with previous measurements, this suggests that there is no spectral break between 230 and 690 GHz. The emission region is thus likely in a transition between optically thick and thin at these frequencies and requires a mix of lepton distributions with varying temperatures consistent with stratification.

We present a statistical analysis of the X-ray flux distribution of Sgr A* from the Chandra X-ray Observatory's 3 Ms Sgr A* X-ray Visionary Project (XVP) in 2012. Our analysis indicates that the observed X-ray flux distribution can be decomposed into a steady quiescent component, represented by a Poisson process with rate $Q=(5.24\pm0.08)\times10^{-3}$ cts s$^{-1},$ and a variable component, represented by a power law process ($dN/dF\propto F^{-\xi},$ $\xi=1.92_{-0.02}^{+0.03}$). This slope matches our recently-reported distribution of flare luminosities. The variability may also be described by a log-normal process with a median unabsorbed 2-8 keV flux of $1.8^{+0.9}_{-0.6}\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ and a shape parameter $\sigma=2.4\pm0.2,$ but the power law provides a superior description of the data. In this decomposition of the flux distribution, all of the intrinsic X-ray variability of Sgr A* (spanning at least three orders of magnitude in flux) can be attributed to flaring activity, likely in the inner accretion flow. We confirm that at the faint end, the variable component contributes ~10% of the apparent quiescent flux, as previously indicated by our statistical analysis of X-ray flares in these Chandra observations. Our flux distribution provides a new and important observational constraint on theoretical models of Sgr A*, and we use simple radiation models to explore the extent to which a statistical comparison of the X-ray and infrared can provide insights into the physics of the X-ray emission mechanism.

The Event Horizon Telescope (EHT), a global sub-millimeter wavelength very long baseline interferometry array, is now resolving the innermost regions around the supermassive black holes Sgr A* and M87. Using black hole images from both simple geometric models and relativistic magnetohydrodynamical accretion flow simulations, we perform a variety of experiments to assess the promise of the EHT for studying strong gravity and accretion physics during the stages of its development. We find that (1) the addition of the LMT and ALMA along with upgraded instrumentation in the "Complete" stage of the EHT allow detection of the photon ring, a signature of Kerr strong gravity, for predicted values of its total flux; (2) the inclusion of coherently averaged closure phases in our analysis dramatically improves the precision of even the current array, allowing (3) significantly tighter constraints on plausible accretion models and (4) detections of structural variability at the levels predicted by the models. While observations at 345 GHz circumvent problems due to interstellar electron scattering in line-of-sight to the galactic center, short baselines provided by CARMA and/or the LMT could be required in order to constrain the overall shape of the accretion flow. Given the systematic uncertainties in the underlying models, using the full complement of two observing frequencies (230 and 345 GHz) and sources (Sgr A* and M87) may be critical for achieving transformative science with the EHT experiment.

Sagittarius A* harbors the supermassive black hole that lies at the dynamical center of our Galaxy. Sagittarius A* spends most of its time in a low luminosity emission state but flares frequently in the infrared and X-ray, increasing up to a few hundred fold in brightness for up to a few hours at a time. The physical processes giving rise to the X-ray flares are uncertain. Here we report the detection with the NuSTAR observatory in Summer and Fall 2012 of four low to medium amplitude X-ray flares to energies up to 79 keV. For the first time, we clearly see that the power-law spectrum of Sagittarius A* X-ray flares extends to high energy, with no evidence for a cut off. Although the photon index of the absorbed power-law fits are in agreement with past observations, we find a difference between the photon index of two of the flares (significant at the 95% confidence level). The spectra of the two brightest flares (~55 times quiescence in the 2-10 keV band) are compared to simple physical models in an attempt to identify the main X-ray emission mechanism, but the data do not allow us to significantly discriminate between them. However, we confirm the previous finding that the parameters obtained with synchrotron models are, for the X-ray emission, physically more reasonable than those obtained with inverse-Compton models. One flare exhibits large and rapid (< 100 s) variability, which, considering the total energy radiated, constrains the location of the flaring region to be within ~10 Schwarzschild radii of the black hole.

Pulsars orbiting the Galactic center black hole, Sgr A*, would be potential probes of its mass, distance and spin, and may even be used to test general relativity. Despite predictions of large populations of both ordinary and millisecond pulsars in the Galactic center, none have been detected within 25 pc by deep radio surveys. One explanation has been that hyperstrong temporal scattering prevents pulsar detections, but the recent discovery of radio pulsations from a highly magnetized neutron star (magnetar) within 0.1 pc shows that the temporal scattering is much weaker than predicted. We argue that an intrinsic deficit in the ordinary pulsar population is the most likely reason for the lack of detections to date: a "missing pulsar problem" in the Galactic center. In contrast, we show that the discovery of a single magnetar implies efficient magnetar formation in the region. If the massive stars in the central parsec form magnetars rather than ordinary pulsars, their short lifetimes could explain the missing pulsars. Efficient magnetar formation could be caused by strongly magnetized progenitors, or could be further evidence of a top-heavy initial mass function. Furthermore, current high-frequency surveys should already be able to detect bright millisecond pulsars, given the measured degree of temporal scattering.

The well-known bluer-when-brighter trend observed in quasar variability is a signature of the complex processes in the accretion disk, and can be a probe of the quasar variability mechanism. Using a sample of 604 variable quasars with repeat spectra in SDSS-I/II, we construct difference spectra to investigate the physical causes of this bluer-when-brighter trend. The continuum of our composite difference spectrum is well-fit by a power-law, with a spectral index in excellent agreement with previous results. We measure the spectral variability relative to the underlying spectra of the quasars, which is independent of any extinction, and compare to model predictions. We show that our SDSS spectral variability results cannot be produced by global accretion rate fluctuations in a thin disk alone. However, we find that a simple model of a inhomogeneous disk with localized temperature fluctuations will produce power-law spectral variability over optical wavelengths. We show that the inhomogeneous disk will provide good fits to our observed spectral variability if the disk has large temperature fluctuations in many independently varying zones, in excellent agreement with independent constraints from quasar microlensing disk sizes, their strong UV spectral continuum, and single-band variability amplitudes. Our results provide an independent constraint on quasar variability models, and add to the mounting evidence that quasar accretion disks have large localized temperature fluctuations.

Magnetically arrested accretion discs (MADs), where the magnetic pressure in the inner disc is dynamically important, provide an alternative mechanism for regulating accretion to what is commonly assumed in black hole systems. We show that a global magnetic field inversion in the MAD state can destroy the jet, significantly increase the accretion rate, and move the effective inner disc edge in to the marginally stable orbit. Reconnection of the MAD field in the inner radii launches a new type of transient outflow containing hot plasma generated by magnetic dissipation. This transient outflow can be as powerful as the steady magnetically-dominated Blandford-Znajek jet in the MAD state. The field inversion qualitatively describes many of the observational features associated with the high luminosity hard to soft state transition in black hole X-ray binaries: the jet line, the transient ballistic jet, and the drop in rms variability. These results demonstrate that the magnetic field configuration can influence the accretion state directly, and hence the magnetic field structure is an important second parameter in explaining observations of accreting black holes across the mass and luminosity scales.

We propose a new model of the steep power law state of luminous black hole X-ray binaries. The model uses the fact that at high luminosities, the inner radii of radiation pressure dominated accretion discs are expected to i) become effectively optically thin and ii) produce significant luminosities. The gas temperature therefore rises sharply inwards, producing local saturated Compton spectra with rapidly increasing peak energy. These spectra sum together to form a steep power law tail to the spectrum. A given photon energy on this tail corresponds to a narrow range in radius, so that local vertical oscillations of the disc naturally produce high quality high-frequency quasi-periodic oscillations (HFQPOs) in the hard X-ray band. The two lowest order modes have a robust frequency ratio of ~1.53. This model explains the appearance of steep power law spectra and HFQPOs at high luminosity, the 3:2 HFQPO frequency ratios, and their association with the power law spectral component. We predict an increase in QPO quality factor when the power spectrum is restricted to a narrower photon energy band, and an increase in HFQPO frequency at higher X-ray energies or lower luminosities. Future X-ray telescopes could detect additional HFQPOs from higher order modes. We demonstrate how this model could be used to measure black hole spin from HFQPOs, and qualitatively estimate the spin of GRO J1655-40 as a/M ~ 0.4 - 0.7.

Over the last decade, X-ray observations of Sgr A* have revealed a black hole in a deep sleep, punctuated roughly once per day by brief flares. The extreme X-ray faintness of this supermassive black hole has been a long-standing puzzle in black hole accretion. To study the accretion processes in the Galactic Center, Chandra (in concert with numerous ground- and space-based observatories) undertook a 3 Ms campaign on Sgr A* in 2012. With its excellent observing cadence, sensitivity, and spectral resolution, this Chandra X-ray Visionary Project (XVP) provides an unprecedented opportunity to study the behavior of the closest supermassive black hole. We present a progress report from our ongoing study of X-ray flares, including the brightest flare ever seen from Sgr A*. Focusing on the statistics of the flares and the quiescent emission, we discuss the physical implications of X-ray variability in the Galactic Center.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

Most supermassive black holes (SMBHs) are accreting at very low levels and are difficult to distinguish from the galaxy centers where they reside. Our own Galaxy's SMBH provides a uniquely instructive exception, and we present a close-up view of its quiescent X-ray emission based on 3 mega-second of Chandra observations. Although the X-ray emission is elongated and aligns well with a surrounding disk of massive stars, we can rule out a concentration of low-mass coronally active stars as the origin of the emission based on the lack of predicted Fe Kalpha emission. The extremely weak H-like Fe Kalpha line further suggests the presence of an outflow from the accretion flow onto the SMBH. These results provide important constraints for models of the prevalent radiatively inefficient accretion state.

We compile and analyse long term (~10 year) submillimetre (1.3, 0.87, 0.43 mm, submm) wavelength light curves of the Galactic centre black hole, Sagittarius A*. The 0.87 and 0.43 mm data are taken from the literature, while the majority of the 1.3 mm light curve is from previously unpublished SMA and CARMA data. We use Monte Carlo simulations to show that on minute to few hour time-scales the variability is consistent with a red noise process with a 230 GHz power spectrum slope of 2.3+0.8-0.6 at 95% confidence. The light curve is de-correlated (white noise) on very long (month to year) times. In order to identify the transition time between red and white noise, we model the light curves as a stochastic damped random walk process. The models allow a quantitative estimate of this physical characteristic time-scale of 8-4+3 hours at 230 GHz at 95% confidence, with consistent results at 345 and 690 GHz. This corresponds to ~10 orbital times or ~1 inflow (viscous) time at R = 3 Rs, a typical radius producing the 230 GHz emission as measured by very long baseline interferometry and found in theoretical accretion flow and jet models. This time-scale is significantly shorter (longer) than those measured at radio (near-infrared, NIR) wavelengths, and is marginally inconsistent with the same variability mechanism operating in the submm and NIR for the expected t ~ R^3/2 scaling. It is crudely consistent with the analogous time-scale inferred in studies of quasar optical light curves after accounting for the difference in emission radius. We find evidence that the submm variability persists at least down to the ISCO, if not the event horizon. These results can be compared quantitatively with similar analyses at different wavebands to test for connections between the variability mechanisms, and with light curves from theoretical models of accreting black holes.

We present the first systematic analysis of the X-ray variability of Sgr A* during the Chandra X-ray Observatory's 2012 Sgr A* X-ray Visionary Project (XVP). With 38 High Energy Transmission Grating Spectrometer (HETGS) observations spaced an average of 7 days apart, this unprecedented campaign enables detailed study of the X-ray emission from this supermassive black hole at high spatial, spectral and timing resolution. In 3 Ms of observations, we detect 39 X-ray flares from Sgr A*, lasting from a few hundred seconds to approximately 8 ks, and ranging in 2-10 keV luminosity from ~1e34 erg/s to 2e35 erg/s. Despite tentative evidence for a gap in the distribution of flare peak count rates, there is no evidence for X-ray color differences between faint and bright flares. Our preliminary X-ray flare luminosity distribution dN/dL is consistent with a power law with index -1.9 (+0.3 -0.4); this is similar to some estimates of Sgr A*'s NIR flux distribution. The observed flares contribute one-third of the total X-ray output of Sgr A* during the campaign, and as much as 10% of the quiescent X-ray emission could be comprised of weak, undetected flares, which may also contribute high-frequency variability. We argue that flares may be the only source of X-ray emission from the inner accretion flow.

The Event Horizon Telescope (EHT), a global very long baseline interferometry array operating at millimetre wavelengths, is spatially resolving the immediate environments of black holes for the first time. The current observations of the Galactic center black hole, Sagittarius A* (Sgr A*), and M87 have been interpreted in terms of either geometric models (e.g., a symmetric Gaussian) or detailed calculations of the appearance of black hole accretion flows. The former are not physically motivated, while the latter are subject to large systematic uncertainties. Motivated by the dominant relativistic effects of Doppler beaming and gravitational lensing in many calculations, we propose a geometric crescent model for black hole images. We show that this simple model provides an excellent statistical description of the existing EHT data of Sgr A* and M87, superior to other geometric models for Sgr A*. It also qualitatively matches physically predicted models, bridging accretion theory and observation. Based on our results, we make predictions for the detectability of the black hole shadow, a signature of strong gravity, in future observations.

Some fraction of the material ejected in a core collapse supernova explosion may remain bound to the compact remnant, and eventually turn around and fall back. We show that the late time (> days) power associated with the accretion of this "fallback" material may significantly affect the optical light curve, in some cases producing super-luminous or otherwise peculiar supernovae. We use spherically symmetric hydrodynamical models to estimate the accretion rate at late times for a range of progenitor masses and radii and explosion energies. The accretion rate onto the proto-neutron star or black hole decreases as Mdot ~ t^-5/3 at late times, but its normalization can be significantly enhanced at low explosion energies, in very massive stars, or if a strong reverse shock wave forms at the helium/hydrogen interface in the progenitor. If the resulting super-Eddington accretion drives an outflow which thermalizes in the outgoing ejecta, the supernova debris will be re-energized at a time when photons can diffuse out efficiently. The resulting light curves are different and more diverse than previous fallback supernova models which ignored the input of accretion power and produced short-lived, dim transients. The possible outcomes when fallback accretion power is significant include super-luminous (> 10^44 ergs / s) Type II events of both short and long durations, as well as luminous Type I events from compact stars that may have experienced significant mass loss. Accretion power may unbind the remaining infalling material, causing a sudden decrease in the brightness of some long duration Type II events. This scenario may be relevant for explaining some of the recently discovered classes of peculiar and rare supernovae.

We present the first spectral energy distributions produced self-consistently by 2.5D general relativistic magneto-hydrodynamical (GRMHD) numerical simulations, where radiative cooling is included in the dynamical calculation. As a case study, we focus on the accretion flow around the supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), which has the best constrained physical parameters. We compare the simulated spectra to the observational data of Sgr A* and explore the parameter space of our model to determine the effect of changing the initial magnetic field configuration, ion to electron temperature ratio T_i/T_e and the target accretion rate. We find the best description of the data for a mass accretion rate of ~ 1e-9 Msun/yr, and rapid spin (0.7 < a_* < 0.9). The submillimeter peak flux seems largely independent of initial conditions, while the higher energies can be very sensitive to the initial magnetic field configuration. Finally, we also discuss flaring features observed in some simulations, that may be due to artifacts of the 2D configuration.

Based on Bremer et al. (2011) and Eckart et al. (2012) we report on simultaneous observations and modeling of the millimeter, near-infrared, and X-ray flare emission of the source Sagittarius A* (SgrA*) associated with the super-massive black hole at the Galactic Center. We study physical processes giving rise to the variable emission of SgrA* from the radio to the X-ray domain. To explain the statistics of the observed variability of the (sub-)mm spectrum of SgrA*, we use a sample of simultaneous NIR/X-ray flare peaks and model the flares using a synchrotron and SSC mechanism. The observations reveal flaring activity in all wavelength bands that can be modeled as the signal from adiabatically expanding synchrotron self-Compton (SSC) components. The model parameters suggest that either the adiabatically expanding source components have a bulk motion larger than v_exp or the expanding material contributes to a corona or disk, confined to the immediate surroundings of SgrA*. For the bulk of the synchrotron and SSC models, we find synchrotron turnover frequencies in the range 300-400 GHz. For the pure synchrotron models this results in densities of relativistic particles of the order of 10^6.5 cm^-3 and for the SSC models, the median densities are about one order of magnitude higher. However, to obtain a realistic description of the frequency-dependent variability amplitude of SgrA*, models with higher turnover frequencies and even higher densities are required. We discuss the results in the framework of possible deviations from equilibrium between particle and magnetic field energy. We also summarize alternative models to explain the broad-band variability of SgrA*.

We present general relativistic magnetohydrodynamic (GRMHD) numerical simulations of the accretion flow around the supermassive black hole in the Galactic centre, Sagittarius A* (Sgr A*). The simulations include for the first time radiative cooling processes (synchrotron, bremsstrahlung, and inverse Compton) self-consistently in the dynamics, allowing us to test the common simplification of ignoring all cooling losses in the modeling of Sgr A*. We confirm that for Sgr A*, neglecting the cooling losses is a reasonable approximation if the Galactic centre is accreting below ~10^{-8} Msun/yr i.e. Mdot < 10^{-7} Mdot_Edd. But above this limit, we show that radiative losses should be taken into account as significant differences appear in the dynamics and the resulting spectra when comparing simulations with and without cooling. This limit implies that most nearby low-luminosity active galactic nuclei are in the regime where cooling should be taken into account. We further make a parameter study of axisymmetric gas accretion around the supermassive black hole at the Galactic centre. This approach allows us to investigate the physics of gas accretion in general, while confronting our results with the well studied and observed source, Sgr A*, as a test case. We confirm that the nature of the accretion flow and outflow is strongly dependent on the initial geometry of the magnetic field. For example, we find it difficult, even with very high spins, to generate powerful outflows from discs threaded with multiple, separate poloidal field loops.

Observations of black hole binaries (BHBs) have established a rich phenomenology of X-ray states. The soft states range from the low variability, accretion disc dominated thermal state (TD) to the higher variability, non-thermal steep power law state (SPL). The disc component in all states is typically modeled with standard thin disc accretion theory. However, this theory is inconsistent with optical/UV spectral, variability, and gravitational microlensing observations of active galactic nuclei (AGNs), the supermassive analogs of BHBs. An inhomogeneous disc (ID) model with large (~0.4 dex) temperature fluctuations in each radial annulus can qualitatively explain all of these AGN observations. The inhomogeneity may be a consequence of instabilities in radiation dominated discs, and therefore may be present in BHBs as well. We show that ID models can explain many features of the TD and SPL states of BHBs. The observed relationships between spectral hardness, disc fraction, and rms variability amplitude in BHBs are reproduced with temperature fluctuations similar to those inferred in AGNs, suggesting a unified picture of luminous accretion discs across orders of magnitude in black hole mass. This picture can be tested with spectral fitting of ID models, X-ray polarization observations, and radiation MHD simulations. If BHB accretion discs are indeed inhomogeneous, only the most disc dominated states (disc fraction > 0.95) can be used to robustly infer black hole spin using current continuum fitting methods.

High-resolution, multi-wavelength, and time-domain observations of the Galactic centre black hole candidate, Sgr A*, allow for a direct test of contemporary accretion theory. To date, all models have assumed alignment between the accretion disc and black hole angular momentum axes, but this is unjustified for geometrically thick accretion flows like that onto Sgr A*. Instead, we calculate images and spectra from a set of simulations of accretion flows misaligned ('tilted') by 15 degrees from the black hole spin axis and compare them with millimetre (mm) to near-infrared (NIR) observations. Non-axisymmetric standing shocks from eccentric fluid orbits dominate the emission, leading to a wide range of possible image morphologies. These effects invalidate previous parameter estimates from model fitting, including estimates of the dimensionless black hole spin, except possibly at low values of spin or tilt. At 1.3mm, the images have crescent morphologies, and the black hole shadow may still be accessible to future mm-VLBI observations. Shock heating leads to high energy electrons (T > 10^12 K), which can naturally produce the observed NIR flux, spectral index, and rapid variability ('flaring'). This NIR emission is uncorrelated with that in the mm, which also agrees with observations. These are the first models to self-consistently explain the time-variable mm to NIR emission of Sgr A*. Predictions of the model include significant structural changes observable with mm-VLBI on both the dynamical (hour) and Lense-Thirring precession (day-year) timescales; and ~30-50 microarcsecond changes in centroid position from extreme gravitational lensing events during NIR flares, detectable with the future VLT instrument GRAVITY. If the observed NIR emission is caused by shock heating in a tilted accretion disc, then the Galactic centre black hole has a positive, non-zero spin parameter (a > 0).

Ongoing millimeter VLBI observations with the Event Horizon Telescope allow unprecedented study of the innermost portion of black hole accretion flows. Interpreting the observations requires relativistic, time-dependent physical modeling. We discuss the comparison of radiative transfer calculations from general relativistic MHD simulations of Sagittarius A* and M87 with current and future mm-VLBI observations. This comparison allows estimates of the viewing geometry and physical conditions of the Sgr A* accretion flow. The viewing geometry for M87 is already constrained from observations of its large-scale jet, but, unlike Sgr A*, there is no consensus for its millimeter emission geometry or electron population. Despite this uncertainty, as long as the emission region is compact, robust predictions for the size of its jet launching region can be made. For both sources, the black hole shadow may be detected with future observations including ALMA and/or the LMT, which would constitute the first direct evidence for a black hole event horizon.

The supermassive black hole candidate at the center of M87 drives an ultra-relativistic jet visible on kiloparsec scales, and its large mass and relative proximity allow for event horizon scale imaging with very long baseline interferometry at millimeter wavelengths (mm-VLBI). Recently, relativistic magneto-hydrodynamic (MHD) simulations of black hole accretion flows have proven capable of launching magnetically-dominated jets. We construct time-dependent disc/jet models of the innermost portion of the M87 nucleus by performing relativistic radiative transfer calculations from one such simulation. We identify two types of models, jet-dominated or disc/jet, that can explain the spectral properties of M87, and use them to make predictions for current and future mm-VLBI observations. The Gaussian source size for the favored sky orientation and inclination from observations of the large-scale jet is 33-44 microarcseconds (~4-6 Schwarzschild radii) on current mm-VLBI telescopes, very similar to existing observations of Sgr A*. The black hole shadow, direct evidence of an event horizon, should be visible in future measurements using baselines between Hawaii and Mexico. Both models exhibit variability at millimeter wavelengths with factor of ~2 amplitudes on year timescales. For the low inclination of M87, the counter-jet dominates the event horizon scale millimeter wavelength emission from the jet-forming region.

Geometrically thick accretion flows may be present in black hole X-ray binaries observed in the low/hard state and in low-luminosity active galactic nuclei. Unlike in geometrically thin disks, the angular momentum axis in these sources is not expected to align with the black hole spin axis. We compute images from three-dimensional general relativistic magnetohydrodynamic simulations of misaligned (tilted) accretion flows using relativistic radiative transfer, and compare the estimated locations of the radiation edge with expectations from their aligned (untilted) counterparts. The radiation edge in the tilted simulations is independent of black hole spin for a tilt of 15 degrees, in stark contrast to the results for untilted simulations, which agree with the monotonic dependence on spin expected from thin accretion disk theory. Synthetic emission line profiles from the tilted simulations depend strongly on the observer's azimuth, and exhibit unique features such as broad "blue wings." Coupled with precession, the azimuthal variation could generate time fluctuations in observed emission lines, which would be a clear "signature" of a tilted accretion flow. Finally, we evaluate the possibility that the observed low- and high-frequency quasi-periodic oscillations (QPOs) from black hole binaries could be produced by misaligned accretion flows. Although low-frequency QPOs from precessing, tilted disks remains a viable option, we find little evidence for significant power in our light curves in the frequency range of high-frequency QPOs.

Active galactic nuclei (AGN) have been observed to vary stochastically with 10-20 rms amplitudes over a range of optical wavelengths where the emission arises in an accretion disk. Since the accretion disk is unlikely to vary coherently, local fluctuations may be significantly larger than the global rms variability. We investigate toy models of quasar accretion disks consisting of a number of regions, n, whose temperatures vary independently with an amplitude of \sigma_T in dex. Models with large fluctuations (\sigma_T=0.35-0.50) in 100-1000 independently fluctuating zones for every factor of two in radius can explain the observed discrepancy between thin accretion disk sizes inferred from microlensing events and optical luminosity while matching the observed optical variability. For the same range of \sigma_T, inhomogeneous disk spectra provide excellent fits to the HST quasar composite without invoking global Compton scattering atmospheres to explain the high levels of observed UV emission. Simulated microlensing light curves for the Einstein cross from our time-varying toy models are well fit using a time-steady power-law temperature disk, and produce magnification light curves that are consistent with current microlensing observations. Deviations due to the inhomogeneous, time-dependent disk structure should occur above the 1% level in the light curves, detectable in future microlensing observations with millimag sensitivity.

Recent high resolution observations of the Galactic center black hole allow for direct comparison with accretion disk simulations. We compare two-temperature synchrotron emission models from three dimensional, general relativistic magnetohydrodynamic simulations to millimeter observations of Sgr A*. Fits to very long baseline interferometry and spectral index measurements disfavor the monochromatic face-on black hole shadow models from our previous work. Inclination angles \le 20 degrees are ruled out to 3 \sigma. We estimate the inclination and position angles of the black hole, as well as the electron temperature of the accretion flow and the accretion rate, to be i=50+35-15 degrees, \xi=-23+97-22 degrees, T_e=(5.4 +/- 3.0)x10^10 K and Mdot=(5+15-2)x10^-9 M_sun / yr respectively, with 90% confidence. The black hole shadow is unobscured in all best fit models, and may be detected by observations on baselines between Chile and California, Arizona or Mexico at 1.3mm or .87mm either through direct sampling of the visibility amplitude or using closure phase information. Millimeter flaring behavior consistent with the observations is present in all viable models, and is caused by magnetic turbulence in the inner radii of the accretion flow. The variability at optically thin frequencies is strongly correlated with that in the accretion rate. The simulations provide a universal picture of the 1.3mm emission region as a small region near the midplane in the inner radii of the accretion flow, which is roughly isothermal and has \nu/\nu_c ~ 1-20, where \nu_c is the critical frequency for thermal synchrotron emission.

The recent VLBI observation of the Galactic center black hole candidate Sgr A* at 1.3mm shows source structure on event-horizon scales. This detection enables a direct comparison of the emission region with models of the accretion flow onto the black hole. We present the first results from time-dependent radiative transfer of general relativistic MHD simulation data, and compare simulated synchrotron images at black hole spin a=0.9 with the VLBI measurements. After tuning the accretion rate to match the millimeter flux, we find excellent agreement between predicted and observed visibilities, even when viewed face-on (i < 30 degrees). VLBI measurements on 2000-3000km baselines should constrain the inclination. The data constrain the accretion rate to be (1.0-2.3)x10^-9 M_sun / yr with 99% confidence, consistent with but independent of prior estimates derived from spectroscopic and polarimetric measurements. Finally, we compute light curves, which show that magnetic turbulence can directly produce flaring events with .5 hour rise times, 2-3.5 hour durations and 40-50% flux modulation, in agreement with observations of Sgr A* at millimeter wavelengths.

Relativistic radiative transfer problems require the calculation of photon trajectories in curved spacetime. We present a novel technique for rapid and accurate calculation of null geodesics in the Kerr metric. The equations of motion from the Hamilton-Jacobi equation are reduced directly to Carlson's elliptic integrals, simplifying algebraic manipulations and allowing all coordinates to be computed semi-analytically for the first time. We discuss the method, its implementation in a freely available FORTRAN code, and its application to toy problems from the literature.