The "Space VLBI 2020: Science and Technology Futures" meeting was the second in The Future of High-Resolution Radio Interferometry in Space series. The first meeting (2018 September 5--6; Noordwijk, the Netherlands) focused on the full range of science applications possible for very long baseline interferometry (VLBI) with space-based antennas. Accordingly, the observing frequencies (wavelengths) considered ranged from below 1~MHz (> 300 m) to above 300~GHz (< 1 mm). For this second meeting, the focus was narrowed to mission concepts and the supporting technologies to enable the highest angular resolution observations at frequencies of 30~GHz and higher (< 1 cm). This narrowing of focus was driven by both scientific and technical considerations. First, results from the RadioAstron mission and the Event Horizon Telescope (EHT) have generated considerable excitement for studying the inner portions of black hole (BH) accretion disks and jets and testing elements of the General Theory of Relativity (GR). Second, the technologies and requirements involved in space-based VLBI differ considerably between 100~MHz and 100~GHz; a related consideration is that there are a number of existing instruments or mission concepts for frequencies of approximately 100~MHz and below, while it has been some time since attention has been devoted to space VLBI at frequencies above 10~GHz. This conference summary attempts to capture elements of presentations and discussions that occurred.

Interferometric visibilities, reflecting the complex correlations between signals recorded at antennas in an interferometric array, carry information about the angular structure of a distant source. While unknown antenna gains in both amplitude and phase can prevent direct interpretation of these measurements, certain combinations of visibilities called closure phases and closure amplitudes are independent of antenna gains and provide a convenient set of robust observables. However, these closure quantities have subtle noise properties and are generally both linearly and statistically dependent. These complications have obstructed the proper use of closure quantities in interferometric analysis, and they have obscured the relationship between analysis with closure quantities and other analysis techniques such as self calibration. We review the statistics of closure quantities, noting common pitfalls that arise when approaching low signal-to-noise due to the nonlinear propagation of statistical errors. We then develop a strategy for isolating and fitting to the independent degrees of freedom captured by the closure quantities through explicit construction of linearly independent sets of quantities along with their noise covariance in the Gaussian limit, valid for moderate signal-to-noise, and we demonstrate that model fits have biased posteriors when this covariance is ignored. Finally, we introduce a unified procedure for fitting to both closure information and partially calibrated visibilities, and we demonstrate both analytically and numerically the direct equivalence of inference based on closure quantities to that based on self calibration of complex visibilities with unconstrained antenna gains.

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 Event Horizon Telescope image of the supermassive black hole in the galaxy M87 is dominated by a bright, unresolved ring. General relativity predicts that embedded within this image lies a thin "photon ring," which is composed of an infinite sequence of self-similar subrings that are indexed by the number of photon orbits around the black hole. The subrings approach the edge of the black hole "shadow," becoming exponentially narrower but weaker with increasing orbit number, with seemingly negligible contributions from high order subrings. Here, we show that these subrings produce strong and universal signatures on long interferometric baselines. These signatures offer the possibility of precise measurements of black hole mass and spin, as well as tests of general relativity, using only a sparse interferometric array.

We report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.

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

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

The Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.

Operating 1.5 million km from Earth at the Sun-Earth L2 Lagrange point, the Origins Space Telescope equipped with a slightly modified version of its HERO heterodyne instrument could function as a uniquely valuable node in a VLBI network. The unprecedented angular resolution resulting from the combination of Origins with existing ground-based millimeter/submillimeter telescope arrays would increase the number of spatially resolvable black holes by a factor of a million, permit the study of these black holes across all of cosmic history, and enable new tests of general relativity by unveiling the photon ring substructure in the nearest black holes.

We present the design and implementation of an automated data calibration and reduction pipeline for very-long-baseline interferometric (VLBI) observations taken at millimeter wavelengths. These short radio-wavelengths provide the best imaging resolution available from ground-based VLBI networks such as the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA), but require specialized processing due to the strong effects from atmospheric opacity and turbulence as well as the heterogeneous nature of existing global arrays. The pipeline builds upon a calibration suite (HOPS) originally designed for precision geodetic VLBI. To support the reduction of data for astronomical observations, we have developed an additional framework for global phase and amplitude calibration which provides output in a standard data format for astronomical imaging and analysis. We test the pipeline on observations taken at 3.5 mm (86 GHz) by the GMVA joined by the phased Atacama Large Millimeter/submillimeter Array in April 2017, and demonstrate the benefits from the specialized processing of high frequency VLBI data with respect to classical analysis techniques.

This paper proposes a cloud architecture for the correlation of wide bandwidth VLBI data. Cloud correlation facilitates processing of entire experiments in parallel using flexibly allocated and practically unlimited compute resources. This approach offers a potential improvement over dedicated correlation clusters, which are constrained by a fixed number of installed processor nodes and playback units. Additionally, cloud storage offers an alternative to maintaining a fleet of hard-disk drives that might be utilized intermittently. We describe benchmarks of VLBI correlation using the DiFX-2.5.2 software on the Google Cloud Platform to assess cloud-based correlation performance. The number of virtual CPUs per Virtual Machine was varied to determine the optimum configuration of cloud resources. The number of stations was varied to determine the scaling of correlation time with VLBI arrays of different sizes. Data transfer rates from Google Cloud Storage to the Virtual Machines performing the correlation were also measured. We also present an example cloud correlation configuration. Current cloud service and equipment pricing data is used to compile cost estimates allowing an approximate economic comparison between cloud and cluster processing. The economic comparisons are based on cost figures which are a moving target, and are highly dependent on factors such as the utilization of cluster and media, which are a challenge to estimate. Our model suggests that shifting to the cloud is an alternative path for high data rate, low duty cycle wideband VLBI correlation that should continue to be explored. In the production phase of VLBI correlation, the cloud has the potential to significantly reduce data processing times and allow the processing of more science experiments in a given year for the petabyte-scale data sets increasingly common in both astronomy and geodesy VLBI applications.

Very-long-baseline interferometry (VLBI) at frequencies above 230 GHz with Earth-diameter baselines gives spatial resolution finer than the ${\sim}50 \mu$as "shadow" of the supermassive black hole at the Galactic Center, Sagittarius A* (Sgr A*). Imaging static and dynamical structure near the "shadow" provides a test of general relativity and may allow measurement of black hole parameters. However, traditional Earth-rotation synthesis is inapplicable for sources (such as Sgr A*) with intra-day variability. Expansions of ground-based arrays to include space-VLBI stations may enable imaging capability on time scales comparable to the prograde innermost stable circular orbit (ISCO) of Sgr A*, which is predicted to be 4-30 minutes, depending on black hole spin. We examine the basic requirements for space-VLBI, and we develop tools for simulating observations with orbiting stations. We also develop a metric to quantify the imaging capabilities of an array irrespective of detailed image morphology or reconstruction method. We validate this metric on example reconstructions of simulations of Sgr A* at 230 and 345 GHz, and use these results to motivate expanding the Event Horizon Telescope (EHT) to include small dishes in Low Earth Orbit (LEO). We demonstrate that high-sensitivity sites such as the Atacama Large Millimeter/Submillimeter Array (ALMA) make it viable to add small orbiters to existing ground arrays, as space-ALMA baselines would have sensitivity comparable to ground-based non-ALMA baselines. We show that LEO-enhanced arrays sample half of the diffraction-limited Fourier plane of Sgr A* in less than 30 minutes, enabling reconstructions of near-horizon structure with normalized root-mean-square error $\lesssim0.3$ on sub-ISCO timescales.

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

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

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.

The Greenland Telescope project has recently participated in an experiment to image the supermassive black hole shadow at the center of M87 using Very Long Baseline Interferometry technique in April of 2018. The antenna consists of the 12-m ALMA North American prototype antenna that was modified to support two auxiliary side containers and to withstand an extremely cold environment. The telescope is currently at Thule Air Base in Greenland with the long-term goal to move the telescope over the Greenland ice sheet to Summit Station. The GLT currently has a single cryostat which houses three dual polarization receivers that cover 84-96 GHz, 213-243 GHz and 271-377 GHz bands. A hydrogen maser frequency source in conjunction with high frequency synthesizers are used to generate the local oscillator references for the receivers. The intermediate frequency outputs of each receiver cover 4-8 GHz and are heterodyned to baseband for digitization within a set of ROACH-2 units then formatted for recording onto Mark-6 data recorders. A separate set of ROACH-2 units operating in parallel provides the function of auto-correlation for real-time spectral analysis. Due to the stringent instrumental stability requirements for interferometry a diagnostic test system was incorporated into the design. Tying all of the above equipment together is the fiber optic system designed to operate in a low temperature environment and scalable to accommodate a larger distance between the control module and telescope for Summit Station. A report on the progress of the above electronics instrumentation system will be provided.

The Event Horizon Telescope (EHT) is a very-long-baseline interferometry (VLBI) experiment that aims to observe supermassive black holes with an angular resolution that is comparable to the event horizon scale. The South Pole occupies an important position in the array, greatly increasing its north-south extent and therefore its resolution. The South Pole Telescope (SPT) is a 10-meter diameter, millimeter-wavelength telescope equipped for bolometric observations of the cosmic microwave background. To enable VLBI observations with the SPT we have constructed a coherent signal chain suitable for the South Pole environment. The dual-frequency receiver incorporates state-of-the-art SIS mixers and is installed in the SPT receiver cabin. The VLBI signal chain also includes a recording system and reference frequency generator tied to a hydrogen maser. Here we describe the SPT VLBI system design in detail and present both the lab measurements and on-sky results.

Centaurus A (Cen A) is a bright radio source associated with the nearby galaxy NGC 5128 where high-resolution radio observations can probe the jet at scales of less than a light-day. The South Pole Telescope (SPT) and the Atacama Pathfinder Experiment (APEX) performed a single-baseline very-long-baseline interferometry (VLBI) observation of Cen A in January 2015 as part of VLBI receiver deployment for the SPT. We measure the correlated flux density of Cen A at a wavelength of 1.4 mm on a $\sim$7000 km (5 G$\lambda$) baseline. Ascribing this correlated flux density to the core, and with the use of a contemporaneous short-baseline flux density from a Submillimeter Array observation, we infer a core brightness temperature of $1.4 \times 10^{11}$ K. This is close to the equipartition brightness temperature, where the magnetic and relativistic particle energy densities are equal. Under the assumption of a circular Gaussian core component, we derive an upper limit to the core size $\phi = 34.0 \pm 1.8~\mu\textrm{as}$, corresponding to 120 Schwarzschild radii for a black hole mass of $5.5 \times 10^7 M_{\odot}$.

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.

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

The Atacama Millimeter/submillimeter Array (ALMA) Phasing Project (APP) has developed and deployed the hardware and software necessary to coherently sum the signals of individual ALMA antennas and record the aggregate sum in Very Long Baseline Interferometry (VLBI) Data Exchange Format. These beamforming capabilities allow the ALMA array to collectively function as the equivalent of a single large aperture and participate in global VLBI arrays. The inclusion of phased ALMA in current VLBI networks operating at (sub)millimeter wavelengths provides an order of magnitude improvement in sensitivity, as well as enhancements in u-v coverage and north-south angular resolution. The availability of a phased ALMA enables a wide range of new ultra-high angular resolution science applications, including the resolution of supermassive black holes on event horizon scales and studies of the launch and collimation of astrophysical jets. It also provides a high-sensitivity aperture that may be used for investigations such as pulsar searches at high frequencies. This paper provides an overview of the ALMA Phasing System design, implementation, and performance characteristics.

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

The Event Horizon Telescope, an Earth-sized interferometer, aims to capture an image of a black hole's event horizon to test the theory of general relativity and probe accretion processes.

The goal of the Event Horizon Telescope (EHT) is to provide spatially resolved images of Sgr A*, the source associated with the Galactic Center black hole. Because Sgr A* varies on timescales short compared to an EHT observing campaign, it is interesting to ask whether variability contains information about the structure and dynamics of the accretion flow. In this paper, we introduce "time-domain filtering", a technique to filter time fluctuating images with specific temporal frequency ranges, and demonstrate the power and usage of the technique by applying it to mock millimeter wavelength images of Sgr A*. The mock image data is generated from General Relativistic Magnetohydrodynamic (GRMHD) simulation and general relativistic ray-tracing method. We show that the variability on each line of sight is tightly correlated with a typical radius of emission. This is because disk emissivity fluctuates on a timescale of order the local orbital period. Time-domain filtered images therefore reflect the model dependent emission radius distribution, which is not accessible in time-averaged images. We show that, in principle, filtered data have the power to distinguish between models with different black hole spins, different disk viewing angles, and different disk orientations in the sky.

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

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

Very long baseline interferometry (VLBI) is a technique for imaging celestial radio emissions by simultaneously observing a source from telescopes distributed across Earth. The challenges in reconstructing images from fine angular resolution VLBI data are immense. The data is extremely sparse and noisy, thus requiring statistical image models such as those designed in the computer vision community. In this paper we present a novel Bayesian approach for VLBI image reconstruction. While other methods often require careful tuning and parameter selection for different types of data, our method (CHIRP) produces good results under different settings such as low SNR or extended emission. The success of our method is demonstrated on realistic synthetic experiments as well as publicly available real data. We present this problem in a way that is accessible to members of the community, and provide a dataset website (vlbiimaging.csail.mit.edu) that facilitates controlled comparisons across algorithms.

The Event Horizon Telescope is a global very-long baseline interferometer capable of probing potential deviations from the Kerr metric, which is believed to provide the unique description of astrophysical black holes. Here we report an updated constraint on the quadrupolar deviation of Sagittarius A* within the context of a radiatively inefficient accretion flow model in a quasi-Kerr background. We also simulate near-future constraints obtainable by the forthcoming eight-station array and show that in this model already a one-day observation can measure the spin magnitude to within $0.005$, the inclination to within $0.09^\circ$, the position angle to within $0.04^\circ$, and the quadrupolar deviation to within $0.005$ at $3\sigma$ confidence. Thus, we are entering an era of high-precision strong gravity measurements.

General Relativity predicts that the emission close to a black hole must be lensed by its strong gravitational field, illuminating the last photon orbit. This results in a dark circular area known as the black hole 'shadow'. The Event Horizon Telescope (EHT) is a (sub)mm VLBI network capable of Schwarzschild-radius resolution on Sagittarius A* (or Sgr A*), the 4 million solar mass black hole at the Galactic Center. The goals of the Sgr A* observations include resolving and measuring the details of its morphology. However, EHT data are sparse in the visibility domain, complicating reliable detailed image reconstruction. Therefore, direct pixel imaging should be complemented by other approaches. Using simulated EHT data from a black hole emission model we consider an approach to Sgr A* image reconstruction based on a simple and computationally efficient analytical model that produces images similar to the synthetic ones. The model consists of an eccentric ring with a brightness gradient and a two-dimensional Gaussian. These elemental forms have closed functional representations in the visibility domain, which lowers the computational overhead of fitting the model to the EHT observations. For model fitting we use a version of the Markov chain Monte-Carlo (MCMC) algorithm based on the Metropolis-Hastings sampler with replica exchange. Over a series of simulations we demonstrate that our model can be used for determining geometric measures of a black hole, thus providing information on the shadow size, linking General Relativity with accretion theory.

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

Images of the linear polarization of synchrotron radiation around Active Galactic Nuclei (AGN) identify their projected magnetic field lines and provide key data for understanding the physics of accretion and outflow from supermassive black holes. The highest resolution polarimetric images of AGN are produced with Very Long Baseline Interferometry (VLBI). Because VLBI incompletely samples the Fourier transform of the source image, any image reconstruction that fills in unmeasured spatial frequencies will not be unique and reconstruction algorithms are required. In this paper, we explore extensions of the Maximum Entropy Method (MEM) to linear polarimetric VLBI imaging. In contrast to previous work, our polarimetric MEM algorithm combines a Stokes I imager that uses only bispectrum measurements that are immune to atmospheric phase corruption with a joint Stokes Q and U imager that operates on robust polarimetric ratios. We demonstrate the effectiveness of our technique on 7- and 3-mm wavelength quasar observations from the VLBA and simulated 1.3-mm Event Horizon Telescope observations of Sgr A* and M87. Consistent with past studies, we find that polarimetric MEM can produce superior resolution compared to the standard CLEAN algorithm when imaging smooth and compact source distributions. As an imaging framework, MEM is highly adaptable, allowing a range of constraints on polarization structure. Polarimetric MEM is thus an attractive choice for image reconstruction with the EHT.

Originally developed to image the shadow region of the central black hole in Sagittarius A* and in the nearby galaxy M87, the Event Horizon Telescope (EHT) provides deep, very high angular resolution data on other AGN sources too. The challenges of working with EHT data have spurred the development of new image reconstruction algorithms. This work briefly reviews the status of the EHT and its utility for observing AGN sources, with emphasis on novel imaging techniques that offer the promise of better reconstructions at 1.3 mm and other wavelengths.

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

The radio emission from Sgr A$^\ast$ is thought to be powered by accretion onto a supermassive black hole of $\sim\! 4\times10^6~ \rm{M}_\odot$ at the Galactic Center. At millimeter wavelengths, Very Long Baseline Interferometry (VLBI) observations can directly resolve the bright innermost accretion region of Sgr A$^\ast$. Motivated by the addition of many sensitive, long baselines in the north-south direction, we developed a full VLBI capability at the Large Millimeter Telescope Alfonso Serrano (LMT). We successfully detected Sgr A$^\ast$ at 3.5~mm with an array consisting of 6 Very Long Baseline Array telescopes and the LMT. We model the source as an elliptical Gaussian brightness distribution and estimate the scattered size and orientation of the source from closure amplitude and self-calibration analysis, obtaining consistent results between methods and epochs. We then use the known scattering kernel to determine the intrinsic two dimensional source size at 3.5 mm: $(147\pm7~\mu\rm{as}) \times (120\pm12~\mu\rm{as})$, at position angle $88^\circ\pm7^\circ$ east of north. Finally, we detect non-zero closure phases on some baseline triangles, but we show that these are consistent with being introduced by refractive scattering in the interstellar medium and do not require intrinsic source asymmetry to explain.

An initial three-station version of the Event Horizon Telescope, a millimeter-wavelength very-long baseline interferometer, has observed Sagittarius A* (Sgr A*) repeatedly from 2007 to 2013, resulting in the measurement of a variety of interferometric quantities. Of particular importance, there is now a large set of closure phases, measured over a number of independent observing epochs. We analyze these observations within the context of a realization of semi-analytic radiatively inefficient disk models, implicated by the low luminosity of Sgr A*. We find a broad consistency among the various observing epochs and between different interferometric data types, with the latter providing significant support for this class of models of Sgr A*. The new data significantly tighten existing constraints on the spin magnitude and its orientation within this model context, finding a spin magnitude of $a=0.10^{+0.30+0.56}_{-0.10-0.10}$, an inclination with respect to the line of sight of $\theta={60^\circ}^{+5^\circ+10^\circ}_{-8^\circ-13^\circ}$, and a position angle of $\xi={156^\circ}^{+10^\circ+14^\circ}_{-17^\circ-27^\circ}$ east of north. These are in good agreement with previous analyses. Notably, the previous $180^\circ$ degeneracy in the position angle has now been conclusively broken by the inclusion of the closure phase measurements. A reflection degeneracy in the inclination remains, permitting two localizations of the spin vector orientation, one of which is in agreement with the orbital angular momentum of the infrared gas cloud G2 and the clockwise disk of young stars. This possibly supports a relationship between Sgr A*'s accretion flow and these larger-scale features.

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.

Magnetic fields are believed to drive accretion and relativistic jets in black hole accretion systems, but the magnetic-field structure that controls these phenomena remains uncertain. We perform general relativistic (GR) polarized radiative transfer of time-dependent three-dimensional GR magnetohydrodynamical (MHD) simulations to model thermal synchrotron emission from the Galactic Center source Sagittarius A$^\ast$ (Sgr A*). We compare our results to new polarimetry measurements by the Event Horizon Telescope (EHT) and show how polarization in the visibility (Fourier) domain distinguishes and constrains accretion flow models with different magnetic field structures. These include models with small-scale fields in disks driven by the magnetorotational instability (MRI) as well as models with large-scale ordered fields in magnetically-arrested disks (MAD). We also consider different electron temperature and jet mass-loading prescriptions that control the brightness of the disk, funnel-wall jet, and Blandford-Znajek-driven funnel jet. Our comparisons between the simulations and observations favor models with ordered magnetic fields near the black hole event horizon in Sgr A*, although both disk- and jet-dominated emission can satisfactorily explain most of the current EHT data. We show that stronger model constraints should be possible with upcoming circular polarization and higher frequency ($349~{\rm GHz}$) measurements.

In general relativity, the angular radius of the shadow of a black hole is primarily determined by its mass-to-distance ratio and depends only weakly on its spin and inclination. If general relativity is violated, however, the shadow size may also depend strongly on parametric deviations from the Kerr metric. Based on a reconstructed image of Sagittarius A* (Sgr A*) from a simulated one-day observing run of a seven-station Event Horizon Telescope (EHT) array, we employ a Markov chain Monte Carlo algorithm to demonstrate that such an observation can measure the angular radius of the shadow of Sgr A* with an uncertainty of ~1.5 uas (6%). We show that existing mass and distance measurements can be improved significantly when combined with upcoming EHT measurements of the shadow size and that tight constraints on potential deviations from the Kerr metric can be obtained.

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

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.

We show that interferometry can be applied to study irregular, rapidly rotating structures, as are expected in the turbulent accretion flow near a black hole. Specifically, we analyze the lagged covariance between interferometric baselines of similar lengths but slightly different orientations. For a flow viewed close to face-on, we demonstrate that the peak in the lagged covariance indicates the direction and angular velocity of the emission pattern from the flow. Even for moderately inclined flows, the covariance robustly estimates the flow direction, although the estimated angular velocity can be significantly biased. Importantly, measuring the direction of the flow as clockwise or counterclockwise on the sky breaks a degeneracy in accretion disk inclinations when analyzing time-averaged images alone. We explore the potential efficacy of our technique using three-dimensional, general relativistic magnetohydrodynamic (GRMHD) simulations, and we highlight several baseline pairs for the Event Horizon Telescope (EHT) that are well-suited to this application. These results indicate that the EHT may be capable of estimating the direction and angular velocity of the emitting material near Sagittarius A*, and they suggest that a rotating flow may even be utilized to improve imaging capabilities.

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}}$.

We report about a 230 GHz very long baseline interferometry (VLBI) fringe finder observation of blazar 3C 279 with the APEX telescope in Chile, the phased submillimeter array (SMA), and the SMT of the Arizona Radio Observatory (ARO). We installed VLBI equipment and measured the APEX station position to 1 cm accuracy (1 sigma). We then observed 3C 279 on 2012 May 7 in a 5 hour 230 GHz VLBI track with baseline lengths of 2800 M$\lambda$ to 7200 M$\lambda$ and a finest fringe spacing of 28.6 micro-arcseconds. Fringes were detected on all baselines with SNRs of 12 to 55 in 420 s. The correlated flux density on the longest baseline was ~0.3 Jy/beam, out of a total flux density of 19.8 Jy. Visibility data suggest an emission region <38 uas in size, and at least two components, possibly polarized. We find a lower limit of the brightness temperature of the inner jet region of about 10^10 K. Lastly, we find an upper limit of 20% on the linear polarization fraction at a fringe spacing of ~38 uas. With APEX the angular resolution of 230 GHz VLBI improves to 28.6 uas. This allows one to resolve the last-photon ring around the Galactic Center black hole event horizon, expected to be 40 uas in diameter, and probe radio jet launching at unprecedented resolution, down to a few gravitational radii in galaxies like M 87. To probe the structure in the inner parsecs of 3C 279 in detail, follow-up observations with APEX and five other mm-VLBI stations have been conducted (March 2013) and are being analyzed.

The 6 billion solar mass supermassive black hole at the center of the giant elliptical galaxy M87 powers a relativistic jet. Observations at millimeter wavelengths with the Event Horizon Telescope have localized the emission from the base of this jet to angular scales comparable to the putative black hole horizon. The jet might be powered directly by an accretion disk or by electromagnetic extraction of the rotational energy of the black hole. However, even the latter mechanism requires a confining thick accretion disk to maintain the required magnetic flux near the black hole. Therefore, regardless of the jet mechanism, the observed jet power in M87 implies a certain minimum mass accretion rate. If the central compact object in M87 were not a black hole but had a surface, this accretion would result in considerable thermal near-infrared and optical emission from the surface. Current flux limits on the nucleus of M87 strongly constrain any such surface emission. This rules out the presence of a surface and thereby provides indirect evidence for an event horizon.

We report the discovery of Faraday rotation toward radio source 3C 84, the active galactic nucleus in NGC1275 at the core of the Perseus Cluster. The rotation measure (RM), determined from polarization observations at wavelengths of 1.3 and 0.9 mm, is (8.7 +/- 2.3) x 10^5 radians/m^2, among the largest ever measured. The RM remained relatively constant over a 2 year period even as the intrinsic polarization position angle wrapped through a span of 300 degrees. The Faraday rotation is likely to originate either in the boundary layer of the radio jet from the nucleus, or in the accretion flow onto the central black hole. The accretion flow probably is disk-like rather than spherical on scales of less than a parsec, otherwise the RM would be even larger.

Observations of the black hole in the center of the Milky Way with the Event Horizon Telescope at 1.3 mm have revealed a size of the emitting region that is smaller than the size of the black-hole shadow. This can be reconciled with the spectral properties of the source, if the accretion flow is seen at a relatively high inclination (50-60 degrees). Such an inclination makes the angular momentum of the flow, and perhaps of the black hole, nearly aligned with the angular momenta of the orbits of stars that lie within 3 arcsec from the black hole. We discuss the implications of such an alignment for the properties of the black hole and of its accretion flow. We argue that future Event-Horizon-Telescope observations will not only refine the inclination of Sgr A* but also measure precisely its orientation on the plane of the sky.

The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ~50 microarcseconds. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry (VLBI). However, strong-field GR features of interest will be blurred at lambda >= 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.

We demonstrate that polarimetric interferometry can be used to extract precise spatial information about compact polarized flares of Sgr A*. We show that, for a faint dynamical component, a single interferometric baseline suffices to determine both its polarization and projected displacement from the quiescent intensity centroid. A second baseline enables two-dimensional reconstruction of the displacement, and additional baselines can self-calibrate using the flare, enhancing synthesis imaging of the quiescent emission. We apply this technique to simulated 1.3-mm wavelength observations of a "hot spot" embedded in a radiatively inefficient accretion disk around Sgr A*. Our results indicate that, even with current sensitivities, polarimetric interferometry with the Event Horizon Telescope can achieve ~5 microarcsecond relative astrometry of compact flaring structures near Sgr A* on timescales of minutes.

A 12-m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline > 9,000 km to achieve an exceptional angular resolution of 20 micro arc sec at 350 GHz, which will enable us to resolve the shadow size of ~40 micro arc sec. The triangle with the longest baselines formed by the GLT, the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile, and the Submillimeter Array (SMA) in Hawaii will play a key role for the M87 observations. We have been working on the image simulations based on realistic conditions for a better understanding of the possible observed images. In parallel, retrofitting of the telescope and the site developments are in progress. Based on three years of opacity monitoring at 225 GHz, our measurements indicate that the site is excellent for submm observations, comparable to the ALMA site. The GLT is also expected to make single-dish observations up to 1.5 THz.

The Event Horizon Telescope (EHT) is a project to assemble a Very Long Baseline Interferometry (VLBI) network of mm wavelength dishes that can resolve strong field General Relativistic signatures near a supermassive black hole. As planned, the EHT will include enough dishes to enable imaging of the predicted black hole "shadow", a feature caused by severe light bending at the black hole boundary. The center of M87, a giant elliptical galaxy, presents one of the most interesting EHT targets as it exhibits a relativistic jet, offering the additional possibility of studying jet genesis on Schwarzschild radius scales. Fully relativistic models of the M87 jet that fit all existing observational constraints now allow horizon-scale images to be generated. We perform realistic VLBI simulations of M87 model images to examine detectability of the black shadow with the EHT, focusing on a sequence of model images with a changing jet mass load radius. When the jet is launched close to the black hole, the shadow is clearly visible both at 230 and 345 GHz. The EHT array with a resolution of 20-30$\mu$as resolution ($\sim$2-4 Schwarzschild radii) is able to image this feature independent of any theoretical models and we show that imaging methods used to process data from optical interferometers are applicable and effective for EHT data sets. We demonstrate that the EHT is also capable of tracing real-time structural changes on a few Schwarzschild radii scales, such as those implicated by VHE flaring activity of M87. While inclusion of ALMA in the EHT is critical for shadow imaging, generally the array is robust against loss of a station.

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.

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 report results from 5-day VLBI observations of the well-known quasar 3C 279 at 1.3 mm (230 GHz) in 2011. The measured nonzero closure phases on triangles including stations in Arizona, California and Hawaii indicate that the source structure is spatially resolved. We find an unusual inner jet direction at scales of $\sim$1 parsec extending along the northwest-southeast direction (PA = $127^{\circ}\pm3^{\circ}$), as opposed to other (previously) reported measurements on scales of a few parsecs showing inner jet direction extending to the southwest. The 1.3 mm structure corresponds closely with that observed in the central region of quasi-simultaneous super-resolution VLBA images at 7 mm. The closure phase changed significantly on the last day when compared with the rest of observations, indicating that the inner jet structure may be variable on daily timescales. The observed new direction of the inner jet shows inconsistency with the prediction of a class of jet precession models. Our observations indicate a brightness temperature of $\sim 8\times10^{10}$ K in the 1.3 mm core, much lower than that at centimeter wavelengths. Observations with better uv coverage and sensitivity in the coming years will allow the discrimination between different structure models and will provide direct images of the inner regions of the jet with 20--30 $\mu$as (5--7 light months) resolution.

Global VLBI imaging at millimeter and sub-millimeter wavelength overcomes the opacity barrier of synchrotron self-absorption in AGN and opens the direct view into sub-pc scale regions not accessible before. Since AGN variability is more pronounced at short millimeter wavelength, mm-VLBI can reveal structural changes in very early stages after outbursts. When combined with observations at longer wavelength, global 3mm and 1mm VLBI adds very detailed information. This helps to determine fundamental physical properties at the jet base, and in the vicinity of super-massive black holes at the center of AGN. Here we present new results from multi-frequency mm-VLBI imaging of OJ287 during a major outburst. We also report on a successful 1.3mm VLBI experiment with the APEX telescope in Chile. This observation sets a new record in angular resolution. It also opens the path towards future mm-VLBI with ALMA, which aims at the mapping of the black hole event horizon in nearby galaxies, and the study of the roots of jets in AGN.

Approximately 10% of active galactic nuclei exhibit relativistic jets, which are powered by accretion of matter onto super massive black holes. While the measured width profiles of such jets on large scales agree with theories of magnetic collimation, predicted structure on accretion disk scales at the jet launch point has not been detected. We report radio interferometry observations at 1.3mm wavelength of the elliptical galaxy M87 that spatially resolve the base of the jet in this source. The derived size of 5.5 +/- 0.4 Schwarzschild radii is significantly smaller than the innermost edge of a retrograde accretion disk, suggesting that the M87 jet is powered by an accretion disk in a prograde orbit around a spinning black hole.

The recent development of a relatively inexpensive 16-Gbps data-recording system based on commercial off-the-shelf technology and open-source software, along with parallel development in broadband Very Long Baseline Interferometry (VLBI) techniques, is enabling dramatically improved sensitivity for both astronomical and geodetic VLBI. The system is described, including the results of a demonstration VLBI experiment that illustrates a number of cutting-edge technologies that can be deployed in the near future to significantly enhance the power of the VLBI technique.

We present the first 1.3 mm (230 GHz) very long baseline interferometry model image of an AGN jet using closure phase techniques with a four-element array. The model image of the quasar 1924-292 was obtained with four telescopes at three observatories: the James Clerk Maxwell Telescope (JCMT) on Mauna Kea in Hawaii, the Arizona Radio Observatory's Submillimeter Telescope (SMT) in Arizona, and two telescopes of the Combined Array for Research in Millimeterwave Astronomy (CARMA) in California in April 2009. With the greatly improved resolution compared with previous observations and robust closure phase measurement, the inner jet structure of 1924-292 was spatially resolved. The inner jet extends to the northwest along a position angle of $-53^\circ$ at a distance of 0.38\,mas from the tentatively identified core, in agreement with the inner jet structure inferred from lower frequencies, and making a position angle difference of $\sim 80^{\circ}$ with respect to the cm-jet. The size of the compact core is 0.15\,pc with a brightness temperature of $1.2\times10^{11}$\,K. Compared with those measured at lower frequencies, the low brightness temperature may argue in favor of the decelerating jet model or particle-cascade models. The successful measurement of closure phase paves the way for imaging and time resolving Sgr A* and nearby AGN with the Event Horizon Telescope.

Dynamical mass measurements to date have allowed determinations of the mass M and the distance D of a number of nearby supermassive black holes. In the case of Sgr A*, these measurements are limited by a strong correlation between the mass and distance scaling roughly as M ~ D^2. Future very-long baseline interferometric (VLBI) observations will image a bright and narrow ring surrounding the shadow of a supermassive black hole, if its accretion flow is optically thin. In this paper, we explore the prospects of reducing the correlation between mass and distance with the combination of dynamical measurements and VLBI imaging of the ring of Sgr A*. We estimate the signal to noise ratio of near-future VLBI arrays that consist of five to six stations, and we simulate measurements of the mass and distance of Sgr A* using the expected size of the ring image and existing stellar ephemerides. We demonstrate that, in this best-case scenario, VLBI observations at 1 mm can improve the error on the mass by a factor of about two compared to the results from the monitoring of stellar orbits alone. We identify the additional sources of uncertainty that such imaging observations have to take into account. In addition, we calculate the angular diameters of the bright rings of other nearby supermassive black holes and identify the optimal targets besides Sgr A* that could be imaged by a ground-based VLBI array or future space-VLBI missions allowing for refined mass measurements.

Millimeter wave Very Long Baseline Interferometry (mm-VLBI) provides access to the emission region surrounding Sagittarius A*, the supermassive black hole at the center of the Milky Way, on sub-horizon scales. Recently, a closure phase of 0+-40 degrees was reported on a triangle of Earth-sized baselines (SMT-CARMA-JCMT) representing a new constraint upon the structure and orientation of the emission region, independent from those provided by the previously measured 1.3mm-VLBI visibility amplitudes alone. Here, we compare this to the closure phases associated with a class of physically motivated, radiatively inefficient accretion flow models, and present predictions for future mm-VLBI experiments with the developing Event Horizon Telescope (EHT). We find that the accretion flow models are capable of producing a wide variety of closure phases on the SMT-CARMA-JCMT triangle, and thus not all models are consistent with the recent observations. However, those models that reproduce the 1.3mm-VLBI visibility amplitudes overwhelmingly have SMT-CARMA-JCMT closure phases between +-30 degrees, and are therefore broadly consistent with all current mm-VLBI observations. Improving station sensitivity by factors of a few, achievable by increases in bandwidth and phasing together multiple antennas at individual sites, should result in physically relevant additional constraints upon the model parameters and eliminate the current 180 degree ambiguity on the source orientation. When additional stations are included, closure phases of order 45--90 degrees are typical. In all cases the EHT will be able to measure these with sufficient precision to produce dramatic improvements in the constraints upon the spin of Sgr A*.

Extension of very long baseline interferometry (VLBI) to observing wavelengths shorter than 1.3mm provides exceptional angular resolution (~20 micro arcsec) and access to new spectral regimes for the study of astrophysical phenomena. To maintain phase coherence across a global VLBI array at these wavelengths requires that ultrastable frequency references be used for the heterodyne receivers at all participating telescopes. Hydrogen masers have traditionally been used as VLBI references, but atmospheric turbulence typically limits (sub) millimeter VLBI coherence times to ~1-30 s. Cryogenic Sapphire Oscillators (CSO) have better stability than Hydrogen masers on these time scale and are potential alternatives to masers as VLBI references. Here, We describe the design, implementation and tests of a system to produce a 10 MHz VLBI frequency standard from the microwave (11.2 GHz) output of a CSO. To improve long-term stability of the new reference, the CSO was locked to the timing signal from the Global Positioning System satellites and corrected for the oscillator aging. The long-term performance of the CSO was measured by comparison against a hydrogen maser in the same laboratory. The superb short-term performance, along with the improved long-term performance achieved by conditioning, makes the CSO a suitable reference for VLBI at wavelengths less than 1.3mm.

Millimeter very-long baseline interferometry (mm-VLBI) provides the novel capacity to probe the emission region of a handful of supermassive black holes on sub-horizon scales. For Sagittarius A* (Sgr A*), the supermassive black hole at the center of the Milky Way, this provides access to the region in the immediate vicinity of the horizon. Broderick et al. (2009) have already shown that by leveraging spectral and polarization information as well as accretion theory, it is possible to extract accretion-model parameters (including black hole spin) from mm-VLBI experiments containing only a handful of telescopes. Here we repeat this analysis with the most recent mm-VLBI data, considering a class of aligned, radiatively inefficient accretion flow (RIAF) models. We find that the combined data set rules out symmetric models for Sgr A*'s flux distribution at the 3.9-sigma level, strongly favoring length-to-width ratios of roughly 2.4:1. More importantly, we find that physically motivated accretion flow models provide a significantly better fit to the mm-VLBI observations than phenomenological models, at the 2.9-sigma level. This implies that not only is mm-VLBI presently capable of distinguishing between potential physical models for Sgr A*'s emission, but further that it is sensitive to the strong gravitational lensing associated with the propagation of photons near the black hole. Based upon this analysis we find that the most probable magnitude, viewing angle, and position angle for the black hole spin are a=0.0(+0.64+0.86), theta=68(+5+9)(-20-28) degrees, and xi=-52(+17+33)(-15-24) east of north, where the errors quoted are the 1-sigma and 2-sigma uncertainties.

Sagittarius A*, the ~4 x 10^6 solar mass black hole candidate at the Galactic Center, can be studied on Schwarzschild radius scales with (sub)millimeter wavelength Very Long Baseline Interferometry (VLBI). We report on 1.3 mm wavelength observations of Sgr A* using a VLBI array consisting of the JCMT on Mauna Kea, the ARO/SMT on Mt. Graham in Arizona, and two telescopes of the CARMA array at Cedar Flat in California. Both Sgr A* and the quasar calibrator 1924-292 were observed over three consecutive nights, and both sources were clearly detected on all baselines. For the first time, we are able to extract 1.3 mm VLBI interferometer phase information on Sgr A* through measurement of closure phase on the triangle of baselines. On the third night of observing, the correlated flux density of Sgr A* on all VLBI baselines increased relative to the first two nights, providing strong evidence for time-variable change on scales of a few Schwarzschild radii. These results suggest that future VLBI observations with greater sensitivity and additional baselines will play a valuable role in determining the structure of emission near the event horizon of Sgr A*.

The Murchison Wide-field Array (MWA) is a low frequency radio telescope, currently under construction, intended to search for the spectral signature of the epoch of re-ionisation (EOR) and to probe the structure of the solar corona. Sited in Western Australia, the full MWA will comprise 8192 dipoles grouped into 512 tiles, and be capable of imaging the sky south of 40 degree declination, from 80 MHz to 300 MHz with an instantaneous field of view that is tens of degrees wide and a resolution of a few arcminutes. A 32-station prototype of the MWA has been recently commissioned and a set of observations taken that exercise the whole acquisition and processing pipeline. We present Stokes I, Q, and U images from two ~4 hour integrations of a field 20 degrees wide centered on Pictoris A. These images demonstrate the capacity and stability of a real-time calibration and imaging technique employing the weighted addition of warped snapshots to counter extreme wide field imaging distortions.

It is shown that the excellent Murchison Radio-astronomy Observatory site allows the Murchison Widefield Array to employ a simple RFI blanking scheme and still calibrate visibilities and form images in the FM radio band. The techniques described are running autonomously in our calibration and imaging software, which is currently being used to process an FM-band survey of the entire southern sky.

Sagittarius A* is the source of near infrared, X-ray, radio, and (sub)millimeter emission associated with the supermassive black hole at the Galactic Center. In the submillimeter regime, Sgr A* exhibits time-variable linear polarization on timescales corresponding to <10 Schwarzschild radii of the presumed 4 million solar mass black hole. In previous work, we demonstrated the potential for total-intensity (sub)millimeter-wavelength VLBI to detect time-variable -- and periodic -- source structure changes in the Sgr A* black hole system using nonimaging analyses. Here we extend this work to include full polarimetric VLBI observations. We simulate full-polarization (sub)millimeter VLBI data of Sgr A* using a hot-spot model that is embedded within an accretion disk, with emphasis on nonimaging polarimetric data products that are robust against calibration errors. Although the source-integrated linear polarization fraction in the models is typically only a few percent, the linear polarization fraction on small angular scales can be much higher, enabling the detection of changes in the polarimetric structure of Sgr A* on a wide variety of baselines. The shortest baselines track the source-integrated linear polarization fraction, while longer baselines are sensitive to polarization substructures that are beam-diluted by connected-element interferometry. The detection of periodic variability in source polarization should not be significantly affected even if instrumental polarization terms cannot be calibrated out. As more antennas are included in the (sub)mm-VLBI array, observations with full polarization will provide important new diagnostics to help disentangle intrinsic source polarization from Faraday rotation effects in the accretion and outflow region close to the black hole event horizon.

Very strong evidence suggests that Sagittarius A*, a compact radio source at the center of the Milky Way, marks the position of a super massive black hole. The proximity of Sgr A* in combination with its mass makes its apparent event horizon the largest of any black hole candidate in the universe and presents us with a unique opportunity to observe strong-field GR effects. Recent millimeter very long baseline interferometric observations of Sgr A* have demonstrated the existence of structures on scales comparable to the Schwarzschild radius. These observations already provide strong evidence in support of the existence of an event horizon. (Sub)Millimeter VLBI observations in the near future will combine the angular resolution necessary to identify the overall morphology of quiescent emission, such as an accretion disk or outflow, with a fine enough time resolution to detect possible periodicity in the variable component of emission. In the next few years, it may be possible to identify the spin of the black hole in Sgr A*, either by detecting the periodic signature of hot spots at the innermost stable circular orbit or parameter estimation in models of the quiescent emission. Longer term, a (sub)millimeter VLBI "Event Horizon Telescope" will be able to produce images of the Galactic center emission to the see the silhouette predicted by general relativistic lensing. These techniques are also applicable to the black hole in M87, where black hole spin may be key to understanding the jet-launching region.

A long standing goal in astrophysics is to directly observe the immediate environment of a black hole with angular resolution comparable to the event horizon. Realizing this goal would open a new window on the study of General Relativity in the strong field regime, accretion and outflow processes at the edge of a black hole, the existence of an event horizon, and fundamental black hole physics (e.g., spin). Steady long-term progress on improving the capability of Very Long Baseline Interferometry (VLBI) at short wavelengths has now made it extremely likely that this goal will be achieved within the next decade. The most compelling evidence for this is the recent observation by 1.3mm VLBI of Schwarzschild radius scale structure in SgrA*, the compact source of radio, submm, NIR and xrays at the center of the Milky Way. SgrA* is thought to mark the position of a ~4 million solar mass black hole, and because of its proximity and estimated mass presents the largest apparent event horizon size of any black hole candidate in the Universe. Over the next decade, existing and planned mm/submm facilities will be combined into a high sensitivity, high angular resolution "Event Horizon Telescope" that will bring us as close to the edge of black hole as we will come for decades. This white paper describes the science case for mm/submm VLBI observations of both SgrA* and M87 (a radio loud AGN of a much more luminous class that SgrA*). We emphasize that while there is development and procurement involved, the technical path forward is clear, and the recent successful observations have removed much of the risk that would normally be associated with such an ambitious project.

The Murchison Widefield Array (MWA) is a dipole-based aperture array synthesis telescope designed to operate in the 80-300 MHz frequency range. It is capable of a wide range of science investigations, but is initially focused on three key science projects. These are detection and characterization of 3-dimensional brightness temperature fluctuations in the 21cm line of neutral hydrogen during the Epoch of Reionization (EoR) at redshifts from 6 to 10, solar imaging and remote sensing of the inner heliosphere via propagation effects on signals from distant background sources,and high-sensitivity exploration of the variable radio sky. The array design features 8192 dual-polarization broad-band active dipoles, arranged into 512 tiles comprising 16 dipoles each. The tiles are quasi-randomly distributed over an aperture 1.5km in diameter, with a small number of outliers extending to 3km. All tile-tile baselines are correlated in custom FPGA-based hardware, yielding a Nyquist-sampled instantaneous monochromatic uv coverage and unprecedented point spread function (PSF) quality. The correlated data are calibrated in real time using novel position-dependent self-calibration algorithms. The array is located in the Murchison region of outback Western Australia. This region is characterized by extremely low population density and a superbly radio-quiet environment,allowing full exploitation of the instrumental capabilities.

Recent millimeter-VLBI observations of Sagittarius A* (Sgr A*) have, for the first time, directly probed distances comparable to the horizon scale of a black hole. This provides unprecedented access to the environment immediately around the horizon of an accreting black hole. We leverage both existing spectral and polarization measurements and our present understanding of accretion theory to produce a suite of generic radiatively inefficient accretion flow (RIAF) models of Sgr A*, which we then fit to these recent millimeter-VLBI observations. We find that if the accretion flow onto Sgr A* is well described by a RIAF model, the orientation and magnitude of the black hole's spin is constrained to a two-dimensional surface in the spin, inclination, position angle parameter space. For each of these we find the likeliest values and their 1-sigma & 2-sigma errors to be a=0(+0.4+0.7), inclination=50(+10+30)(-10-10) degrees, and position angle=-20(+31+107)(-16-29) degrees, when the resulting probability distribution is marginalized over the others. The most probable combination is a=0(+0.2+0.4), inclination=90(-40-50) degrees and position angle=-14(+7+11)(-7-11) degrees, though the uncertainties on these are very strongly correlated, and high probability configurations exist for a variety of inclination angles above 30 degrees and spins below 0.99. Nevertheless, this demonstrates the ability millimeter-VLBI observations, even with only a few stations, to significantly constrain the properties of Sgr A*.

The recent detection of Sagittarius A* at lambda = 1.3 mm on a baseline from Hawaii to Arizona demonstrates that millimeter wavelength very long baseline interferometry (VLBI) can now spatially resolve emission from the innermost accretion flow of the Galactic center region. Here, we investigate the ability of future millimeter VLBI arrays to constrain the spin and inclination of the putative black hole and the orientation of the accretion disk major axis within the context of radiatively inefficient accretion flow (RIAF) models. We examine the range of baseline visibility and closure amplitudes predicted by RIAF models to identify critical telescopes for determining the spin, inclination, and disk orientation of the Sgr A* black hole and accretion disk system. We find that baseline lengths near 3 gigalambda have the greatest power to distinguish amongst RIAF model parameters, and that it will be important to include new telescopes that will form north-south baselines with a range of lengths. If a RIAF model describes the emission from Sgr A*, it is likely that the orientation of the accretion disk can be determined with the addition of a Chilean telescope to the array. Some likely disk orientations predict detectable fluxes on baselines between the continental United States and even a single 10-12 m dish in Chile. The extra information provided from closure amplitudes by a four-antenna array enhances the ability of VLBI to discriminate amongst model parameters.

Advances in VLBI instrumentation now allow wideband recording that significantly increases the sensitivity of short wavelength VLBI observations. Observations of the super-massive black hole candidate at the center of the Milky Way, SgrA*, with short wavelength VLBI reduces the scattering effects of the intervening interstellar medium, allowing observations with angular resolution comparable to the apparent size of the event horizon of the putative black hole. Observations in April 2007 at a wavelength of 1.3mm on a three station VLBI array have now confirmed structure in SgrA* on scales of just a few Schwarzschild radii. When modeled as a circular Gaussian, the fitted diameter of SgrA* is 37 micro arcsec (+16,-10; 3-sigma), which is smaller than the expected apparent size of the event horizon of the Galactic Center black hole. These observations demonstrate that mm/sub-mm VLBI is poised to open a new window onto the study of black hole physics via high angular resolution observations of the Galactic Center.

The super massive black hole candidate, Sagittarius A*, exhibits variability from radio to X-ray wavelengths on time scales that correspond to < 10 Schwarzschild radii. We survey the potential of millimeter-wavelength VLBI to detect and constrain time variable structures that could give rise to such variations, focusing on a model in which an orbiting hot spot is embedded in an accretion disk. Non-imaging algorithms are developed that use interferometric closure quantities to test for periodicity, and applied to an ensemble of hot-spot models that sample a range of parameter space. We find that structural periodicity in a wide range of cases can be detected on most potential VLBI arrays using modern VLBI instrumentation. Future enhancements of mm/sub-mm VLBI arrays including phased array processors to aggregate VLBI station collecting area, increased bandwidth recording, and addition of new VLBI sites all significantly aid periodicity detection. The methods described herein can be applied to other models of Sagittarius A*, including jet outflows and Magneto-Hydrodynamic accretion simulations.

The cores of most galaxies are thought to harbour supermassive black holes, which power galactic nuclei by converting the gravitational energy of accreting matter into radiation (ref 1). Sagittarius A*, the compact source of radio, infrared and X-ray emission at the centre of the Milky Way, is the closest example of this phenomenon, with an estimated black hole mass that is 4 million times that of the Sun (refs. 2,3). A long-standing astronomical goal is to resolve structures in the innermost accretion flow surrounding Sgr A* where strong gravitational fields will distort the appearance of radiation emitted near the black hole. Radio observations at wavelengths of 3.5 mm and 7 mm have detected intrinsic structure in Sgr A*, but the spatial resolution of observations at these wavelengths is limited by interstellar scattering (refs. 4-7). Here we report observations at a wavelength of 1.3 mm that set a size of 37 (+16, -10; 3-sigma) microarcseconds on the intrinsic diameter of Sgr A*. This is less than the expected apparent size of the event horizon of the presumed black hole, suggesting that the bulk of SgrA* emission may not be not centred on the black hole, but arises in the surrounding accretion flow.

Multiwavelength monitoring observations of Sagittarius A* exhibit variability on timescales of minutes to hours, indicating emission regions localized near the event horizon. (Sub)Millimeter-wavelength VLBI is uniquely suited to probe the environment of the assumed black hole on these scales. We consider a range of orbiting hot-spot and accretion-disk models and find that periodicity in Sgr A* flares is detectable using closure quantities. Our methods are applicable to any model producing source structure changes near the black hole, including jets and magnetohydrodynamic disk instabilities, and suggest that (sub)millimeter VLBI will play a prominent role in investigating Sgr A* near the event horizon.

We present high sensitivity Very Long Baseline Interferometry (VLBI) observations 806 days after the gamma-ray burst of 2003 March 29 (GRB030329). The angular diameter of the radio afterglow is measured to be 0.347 +- 0.09 mas, corresponding to 0.99 +- 0.26 pc at the redshift of GRB030329 (z = 0.1685). The evolution of the image size favors a uniform external density over an R^-2 wind-like density profile (at distances of R >~10^18 cm from the source), although the latter cannot be ruled out yet. The current apparent expansion velocity of the image size is only mildly relativistic, suggesting a non-relativistic transition time of t_NR ~ 1 yr. A rebrightening, or at least a significant flattening in the flux decay, is expected within the next several years as the counter-jet becomes visible (this has not yet been observed). An upper limit of <1.9c is set on the proper motion of the flux centroid.

Experiments were performed with prototype antenna tiles for the Mileura Widefield Array--Low Frequency Demonstrator (MWA-LFD) to better understand the widefield, wideband properties of their design and to characterize the radio frequency interference (RFI) between 80 and 300 MHz at the site in Western Australia. Observations acquired during the six month deployment confirmed the predicted sensitivity of the antennas, sky-noise dominated system temperatures, and phase-coherent interferometric measurements. The radio spectrum is remarkably free of strong terrestrial signals, with the exception of two narrow frequency bands allocated to satellite downlinks and rare bursts due to ground-based transmissions being scattered from aircraft and meteor trails. Results indicate the potential of the MWA-LFD to make significant achievements in its three key science objectives: epoch of reionziation science, heliospheric science, and radio transient detection.

We report on new developments in VLBI, with emphasis on experiments performed at the highest frequencies possible to date (so called mm-VLBI). We have observed the nucleus of M 87 (Virgo A) with global VLBI at 3 mm. We show a new image of the inner-most jet region with an angular resolution of approx. 300 x 60 micro-arcseconds. In terms of Schwarzschild radii, this leads to an upper limit of the jet base of approx. 100 x 20 Schwarzschild radii. We also report on two VLBI pilot-experiments, which demonstrate the technical feasibility of global VLBI at 150 and 230 GHz (2 mm and 1.3 mm). The experiments lead to upper limits to the size of the unresolved AGN-cores in the 25 - 30 micro-arcsecond range. The participation of new and near-future mm-telescopes (like APEX, CARMA, SMA, LMT, ALMA, etc.) in global mm-VLBI will provide the necessary sensitivity for the imaging of black holes and their immediate environment.

We summarize the present status of VLBI experiments at 3 mm (86 GHz), 2 mm (129-150 GHz) and 1.3 mm (215-230 GHz). We present and discuss a new 3 mm VLBI map of M87 (Virgo A), which has a spatial resolution of only approx. 20 Schwarzschild radii. We discuss recent VLBI results for SgrA* and argue in favor of new observations within an extended European mm-VLBI network, in order to search for variability. We discuss the possibilities to image the `event horizon' of a super-massive black hole at wavelengths < 2mm, and conclude that the addition of large and sensitive millimetre telescopes such as CARMA, the SMA, the LMT and ALMA will be crucial for this.

We present the first contemporaneous 43GHz and 86GHz VLBI images of the v=1 J=2-1 and J=1-0 SiO masers in the Orion-KL nebula. Both maser species exhibit the same general morphology of earlier J=1-0 maser images which appear to trace the edges of a bi-polar conical outflow. Surprisingly, the J=2-1 masers form further from the central protostar than the J=1-0 masers, a fact not readily explained by current SiO maser pumping models. The average magnitude of offsets between corresponding regions of the two masing transitions is approximately 14% of the total radial extent of the SiO maser emission. This offset indicates that each transition must trace different physical conditions.

In April 2002 an array of antennas operating at 129GHz successfully detected VLBI fringes on both continuum AGN and SiO maser sources. The 129GHz fringes on maser sources represent the highest frequency spectral line VLBI detections to date. The AGN 3C279 was detected on long baselines at both 129GHz (and at 147GHz, see Krichbaum et al in these proceedings) yielding fringe spacings of 50-56 micro arc seconds, an angular resolution record. The array consisted of the University of Arizona Kittpeak 12m antenna, the Heinrich Hertz 10m Telescope (HHT), and the IRAM 30m dish on Pico Veleta. At 129GHz, a number of evolved stars and several young stellar objects exhibit strong SiO maser emission in the v=1 J=3-2 transition. Preliminary cross power spectra of SiO masers around the red hypergiant VYCMa on the HHT-KittPeak baseline ~190km are consistent with multiple spatially separate maser spots associated with the star. Future observations will include continuum observations of the radio source at the Galactic Center, SgrA*, and higher frequency maser lines including HCN and methanol.

At 147 GHz (2mm wavelength), we detected three prominent AGN (NRAO150, 3C279, 1633+382) with Very Long Baseline Interferometry (VLBI) with an angular resolution of only 18 micro-arcseconds on the baseline between two antennas in Arizona (10m HHT and 12m KittPeak) and the IRAM 30m antenna on Pico Veleta in Spain. This is a new world record in radio interferometry and astronomical imaging and opens fascinating future possibilities to directly image and study the innermost regions in Quasars and other Active Galactic Nuclei.

At radio wavelengths, images of the compact radio source Sagittarius A* (Sgr A*) in the Galactic Center are scatter broadened with a lambda^2 dependence due to an intervening ionized medium. We present VLBI observations of Sgr A* at 86 GHz using a six station array including the VLBA antennas at Pie Town, Fort Davis and Los Alamos, the 12m antenna at Kitt Peak and the millimeter arrays at Hat Creek and Owens Valley. To avoid systematic errors due to imperfect antenna calibration, the data were modeled using interferometric closure information. The data are best modeled by a circular Gaussian brightness distribution of FWHM 0.18 +- 0.02 mas. The data are also shown to be consistent with an elliptical model corresponding to the scattering of a point source. The source structure in the N-S direction, which is less well determined than in the E-W direction due to the limited N-S (u,v) coverage of the array, is constrained to be less than 0.27 mas by these measurements. These results are consistent with extrapolations of intrinsic structure estimates obtained with VLBI at 7mm wavelength assuming the intrinsic size of Sgr A* has a greater dependence than lambda^0.9 with wavelength.

While VLBI experiments at 3mm are nowadays fairly standard, VLBI at 1mm still is technically difficult. Here we summarize the results from previous VLBI pilot experiments at 1.3mm. These included antennas at IRAM (Pico Veleta, Plateau de Bure) and more recently the Heinrich Hertz telescope on Mt. Graham, Arizona.