Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

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

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

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

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.

Very long baseline interferometry (VLBI) from the ground at millimeter wavelengths can resolve the black hole shadow around two supermassive black holes, Sagittarius A* and M87. The addition of modest telescopes in space would allow the combined array to produce higher-resolution, higher-fidelity images of these and other sources. This paper explores the potential benefits of adding orbital elements to the Event Horizon Telescope. We reconstruct model images using simulated data from arrays including telescopes in different orbits. We find that an array including one telescope near geostationary orbit and one in a high-inclination medium Earth or geosynchronous orbit can succesfully produce high-fidelity images capable of resolving shadows as small as 3 microarcseconds in diameter. One such key source, the Sombrero Galaxy, may be important to address questions regarding why some black holes launch powerful jets while others do not. Meanwhile, higher-resolution imaging of the substructure of M87 may clarify how jets are launched in the first place. The extra resolution provided by space VLBI will also improve studies of the collimation of jets from active galactic nuclei.

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.

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.

Many undergraduate radio astronomy courses are unable to give a detailed treatment of aperture synthesis due to time constraints and limited math backgrounds of students. We have taken a laboratory-based approach to teaching radio interferometry using a set of college-level, table-top exercises. These are performed with the Very Small Radio Telescope (VSRT), an interferometer developed at the Haystack Observatory using satellite TV electronics as detectors and compact fluorescent light bulbs as microwave signal sources. The hands-on experience provided by the VSRT in these labs allows students to gain a conceptual understanding of radio interferometry and aperture synthesis without the rigorous mathematical background traditionally required. The data are quickly and easily processed using a user-friendly data analysis Java package, VSRTI\_Plotter.jar. This software can also be used in the absence of the equipment as an interactive computer activity to demonstrate an interferometer's responses to assorted surface brightness distributions. The students also gain some familiarity with Fourier transforms and an appreciation for the Fourier relations in interferometry using another Java package, the Tool for Interactive Fourier Transforms (TIFT). We have successfully used these tools in multiple offerings of our radio astronomy course at Union College

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.

We propose a new superresolution imaging technique for interferometry using sparse modeling, utilizing two regularization terms: the $\ell_1$-norm and a new function named Total Squared Variation (TSV) of the brightness distribution. TSV is an edge-smoothing variant of Total Variation (TV), leading to reducing the sum of squared gradients. First, we demonstrate that our technique may achieve super-resolution of $\sim 30$% compared to the traditional CLEAN beam size using synthetic observations of two point sources. Second, we present simulated observations of three physically motivated static models of Sgr A* with the Event Horizon Telescope (EHT) to show the performance of proposed techniques in greater detail. We find that $\ell_1$+TSV regularization outperforms $\ell_1$+TV regularization with the popular isotropic TV term and the Cotton-Schwab CLEAN algorithm, demonstrating that TSV is well-matched to the expected physical properties of the astronomical images, which are often nebulous. Remarkably, in both the image and gradient domains, the optimal beam size minimizing root-mean-squared errors is $\lesssim 10$% of the traditional CLEAN beam size for $\ell_1$+TSV regularization, and non-convolved reconstructed images have smaller errors than beam-convolved reconstructed images. This indicates that the traditional post-processing technique of Gaussian convolution in interferometric imaging may not be required for the $\ell_1$+TSV regularization. We also propose a feature extraction method to detect circular features from the image of a black hole shadow with the circle Hough transform (CHT) and use it to evaluate the performance of the image reconstruction. With our imaging technique and the CHT, the EHT can constrain the radius of the black hole shadow with an accuracy of $\sim 10-20$% in present simulations for Sgr A*.

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

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.

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.

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.

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 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 measurement of Zeeman splitting in spectral lines---both in emission and absorption---can provide direct estimates of the magnetic field strength and direction in atomic and molecular clouds, both in our own Milky Way and in external galaxies. This method will probe the magnetic field in the warm and cold neutral components of the interstellar medium, providing a complement to the extensive SKA Faraday studies planning to probe the field in the ionized components.

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

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.

We report on the observation of the 36 GHz methanol maser line in the star forming region DR21W to accurately measure the Zeeman effect. The reported Zeeman signature by Fish et al. (2011) became suspicious after an instrumental effect was discovered in the early days of the Very Large Array Wide-band Digital Architecture (WIDAR) correlator commissioning. We conclude that the previously reported magnetic field strength of 58 mG ((1.7 Hz/mG)/z) is instrumental in nature and thus incorrect. With the improved performance of the array, we now deduce a 3 sigma limit of -4.7 to +0.4 mG ((1.7 Hz/mG)/z) for the line-of-sight component of the magnetic field strength in DR21W.

We have developed a set of college level, table-top labs that can be performed with an interferometer using satellite TV electronics and compact fluorescent lamps as microwave signal sources. This interferometer, which was originally developed at the MIT Haystack Observatory as a Very Small Radio Telescope (VSRT) to observe the Sun, provides students with hands-on experience in the fundamentals of radio interferometry. These labs are easily performed and convey an intuitive sense of how combining the signals from an array of antennas reveals information about the structure of a radio source. We have also developed a package of java programs, called "VSRTI Plotter", which is available as a free-download, to facilitate the data processing and analysis of these labs.

The star-forming region W75N hosts bright OH masers that are observed to be variable. We present observations taken in 2008 of the ground-state OH maser transitions with the Very Long Baseline Array (VLBA) and the Multi-Element Radio-Linked Interferometer Network (MERLIN) and with the Nancay Radio Telescope in 2011. Several of the masers in W75N were observed to be flaring, with the brightest 1720-MHz maser in excess of 400 Jy. The 1720-MHz masers appear to be associated with the continuum source VLA 1, unlike the bright flaring 1665- and 1667-MHz masers, which are associated with VLA 2. The 1720-MHz masers are located in an outflow traced by water masers and are indicative of very dense molecular material near the H II region. The magnetic field strengths are larger in the 1720-MHz maser region than in most regions hosting only main-line OH masers. The density falls off along the outflow, and the order of appearance of different transitions of OH masers is consistent with theoretical models. The 1665- and 1667-MHz VLBA data are compared against previous epochs over a time baseline of over 7 years. The median maser motion is 3.5 km/s, with a scatter that is comparable to thermal turbulence. The general pattern of maser proper motions observed in the 1665- and 1667-MHz transitions is consistent with previous observations.

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

We report on the detection of 44 GHz Class I methanol (CH3OH) maser emission in the SgrA complex with the Expanded Very Large Array (EVLA). These EVLA observations show that the SgrA complex harbors at least four different tracers of shocked regions in the radio regime. The 44 GHz masers correlate with the positions and velocities of previously detected 36 GHz CH3OH masers, but less with 1720 MHz OH masers. Our detections agree with theoretical predictions that the densities and temperatures conducive for 1720 MHz OH masers may also produce 36 and 44 GHz CH3OH maser emission. However, many 44 GHz masers do not overlap with 36 GHz methanol masers, suggesting that 44 GHz masers also arise in regions too hot and too dense for 36 GHz masers to form. This agrees with the non-detection of 1720 MHz OH masers in the same area, which are thought to be excited under even cooler and less dense conditions. We speculate that the geometry of the 36 GHz masers outlines the current location of a shock front.

Class I methanol masers are believed to be produced in the shock-excited environment around star-forming regions. Many authors have argued that the appearance of various subsets of class I masers may be indicative of specific evolutionary stages of star formation or excitation conditions. Until recently, however, no major interferometer was capable of imaging the important 36 GHz transition. We report on Expanded Very Large Array observations of the 36 GHz methanol masers and Submillimeter Array observations of the 229 GHz methanol masers in DR21(OH), DR21N, and DR21W. The distribution of 36 GHz masers in the outflow of DR21(OH) is similar to that of the other class I methanol transitions, with numerous multitransition spatial overlaps. At the site of the main continuum source in DR21(OH), class I masers at 36 and 229 GHz are found in virtual overlap with class II 6.7 GHz masers. To the south of the outflow, the 36 GHz masers are scattered over a large region but usually do not appear coincident with 44 GHz masers. In DR21W we detect an "S-curve" signature in Stokes V that implies a large value of the magnetic field strength if interpreted as due to Zeeman splitting, suggesting either that class I masers may exist at higher densities than previously believed or that the direct Zeeman interpretation of S-curve Stokes V profiles in class I masers may be incorrect. We find a diverse variety of different maser phenomena in these sources, suggestive of differing physical conditions among them.

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

We present the first detection of a 6.7 GHz Class II methanol (CH3OH) maser in the Andromeda galaxy (M31). The CH3OH maser was found in a Very Large Array (VLA) survey during the fall of 2009. We have confirmed the methanol maser with the new Expanded VLA (EVLA), in operation since March 2010, but were unsuccessful in detecting a water maser at this location. A direct application for this methanol maser is the determination of the proper motion of M31, such as was obtained with water masers in M33 and IC10 previously. Unraveling the three-dimensional velocity of M31 would solve for the biggest unknown in the modeling of the dynamics and evolution of the Local Group of galaxies.

We present global VLBI observations of the first-excited state OH masers in the massive star-forming region Onsala 1 (ON 1). The 29 masers detected are nearly all from the 6035 MHz transition, and nearly all are identifiable as Zeeman pair components. The 6030 and 6035 MHz masers are coincident with previously published positions of ground-state masers to within a few milliarcseconds, and the magnetic fields deduced from Zeeman splitting are comparable. The 6.0 GHz masers in ON 1 are always found in close spatial association with 1665 MHz OH masers, in contrast to the situation in the massive star-forming region W3(OH), suggesting that extreme high density OH maser sites (excited-state masers with no accompanying ground-state maser, as seen in W3(OH)) are absent from ON 1. The large magnetic field strength among the northern, blueshifted masers is confirmed. The northern masers may trace an outflow or be associated with an exciting source separate from the other masers, or the relative velocities of the northern and southern masers may be indicative of expansion and rotation. High angular resolution observations of nonmasing material will be required to understand the complex maser distribution in ON 1.

We report on the interferometric detection of 36.2 GHz Class I methanol emission with the new 27-40 GHz Ka band receivers available on the Expanded Very Large Array (EVLA). The brightness temperatures of the interferometric 36 GHz detections unambiguously indicate for the first time that the emission is maser emission. The 36 GHz methanol masers are not co-spatial with 1720 MHz OH masers, indicating that the two species trace different shocks. The 36 GHz and 44 GHz methanol masers, which both are collisionally pumped, do not necessarily co-exist and may trace different methanol gas. The methanol masers seem correlated with NH_3(3,3) density peaks. We favor an explanation in which the 36 GHz Class I methanol masers outline regions of cloud-cloud collisions, perhaps just before the onset of the formation of individual massive stars. The transition of the Very Large Array (VLA) to the EVLA is well under way, and these detections demonstrate the bright future of this completely renewed instrument.

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.

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.

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.

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.

The collisionally pumped, ground-state 1720 MHz maser line of OH is widely recognized as a tracer for shocked regions and observed in star forming regions and supernova remnants. Whereas some lines of excited states of OH have been detected and studied in star forming regions, the subject of excited-state OH in supernova remnants -- where high collision rates are to be expected -- is only recently being addressed. Modeling of collisional excitation of OH demonstrates that 1720, 4765 and 6049 MHz masers can occur under similar conditions in regions of shocked gas. In particular, the 6049 and 4765 MHz masers become more significant at increased OH column densities where the 1720 MHz masers begin to be quenched. In supernova remnants, the detection of excited-state OH line maser emission could therefore serve as a probe of regions of higher column densities. Using the Very Large Array, we searched for excited-state OH in the 4.7, 7.8, 8.2 and 23.8 GHz lines in four well studied supernova remnants with strong 1720 MHz maser emission (SgrAEast, W28, W44 and IC443). No detections were made, at typical detection limits of around 10 mJy/beam. The search for the 6 GHz lines were done using Effelsberg since the VLA receivers did not cover those frequencies, and are reported on in an accompanying letter (Fish, Sjouwerman & Pihlstrom 2007). We also cross-correlated the positions of known supernova remnants with the positions of 1612 MHz maser emission obtained from blind surveys. No probable associations were found, perhaps except in the SgrAEast region. The lack of detections of excited-state OH indicates that the OH column densities suffice for 1720 MHz inversion but not for inversion of excited-state transitions, consistent with the expected results for C-type shocks.

While masers in the 1720 MHz transition of OH are detected toward many supernova remnants (SNRs), no other OH transition is seen as a maser in SNRs. We present a search for masers at 6049 MHz, which has recently been predicted to produce masers by pure collisional excitation at conditions similar to that required for 1720 MHz masing. The Effelsberg 100 m telescope was used to observe the excited-state 6016, 6030, 6035, and 6049 MHz lines of OH toward selected SNRs, most of which have previously-detected bright 1720 MHz masers. No excited-state masers are found toward SNRs, consistent with previous observations of the 6049 MHz and other excited-state transitions. We do not see clear evidence of absorption toward SNR target positions, although we do see evidence of absorption in the molecular cloud at +50 km/s near Sgr A East. Weak absorption is detected at 6016 MHz toward W3(OH), while stronger, narrower emission is seen at 6049 MHz, suggesting that the 6049 MHz emission is a low-gain maser. We conclude that conditions in SNRs are not conducive to excited-state maser emission, especially in excited-state satellite lines.

This Letter reports on initial Expanded Very Large Array (EVLA) observations of the 6035 MHz masers in ON 1. The EVLA data are of good quality, lending confidence in the new receiver system. Nineteen maser features, including six Zeeman pairs, are detected. The overall distribution of 6035 MHz OH masers is similar to that of the 1665 MHz OH masers. The spatial resolution is sufficient to unambiguously determine that the magnetic field is strong (~ -10 mG) at the location of the blueshifted masers in the north, consistent with Zeeman splitting detected in 13441 MHz OH masers in the same velocity range. Left and right circularly polarized ground-state features dominate in different regions in the north of the source, which may be due to a combination of magnetic field and velocity gradients. The combined distribution of all OH masers toward the south is suggestive of a shock structure of the sort previously seen in W3(OH).

We present data of proper motions of OH masers in the massive star-forming regions ON 1, K3-50, and W51 Main/South. OH maser motions in ON 1 are consistent with expansion at approximately 5 km/s, likely tracing the expanding ultracompact H II region. Motions in K3-50 are faster and may be indicating the final stages of OH maser emission in the source, before the OH masers turn off as the H II region transitions from the ultracompact to the compact phase. W51 South shows indications of aspherical expansion, while motions in W51 Main are more difficult to interpret. Nevertheless, it appears that the relative projected separation between W51 Main and W51 South is decreasing, corresponding to an estimate of enclosed mass of at least 1500 solar masses, consistent with estimates derived from millimeter-wavelength dust emission. We confirm the ~20 mG magnetic fields previously seen in W51 Main, which may represent the upper end of the density range allowable for 1665 MHz maser emission. Magnetic field strengths and directions, obtained from Zeeman splitting, in each source are consistent with values obtained in the first epoch four to nine years ago.

Excited-state OH maser emission has previously been reported in the circumstellar envelopes of only two evolved stars: the Mira star AU Geminorum and the hypergiant NML Cygni. We present Very Large Array (VLA) observations of the 1665, 1667, and excited-state 4750 MHz mainline OH transitions in AU Gem and Expanded Very Large Array (EVLA) observations of the excited-state 6030 and 6035 MHz OH mainline transitions in NML Cyg. We detect masers in both mainline transitions in AU Gem but no excited-state emission in either star. We conclude that the excited-state OH emission in AU Gem is either a transient phenomenon (such as for NML Cyg outlined below), or possibly an artifact in the data, and that the excited state OH emission in NML Cyg was generated by an episode of enhanced shock between the stellar mass-loss and an outflow of the Cyg OB2 association. With these single exceptions, it therefore appears that excited-state OH emission indeed should not be predicted nor observable in evolved stars as part of their normal structure or evolution.

The massive star-forming region G11.904-0.141 is one of only 11 sources to show maser emission in the highly-excited 13441 MHz transition of OH. VLBA observations of the 1665, 1667, 4765, and 13441 MHz transitions of OH toward G11.904-0.141 are presented. Masers are detected at 1665, 1667, and 4765 MHz, but the 13441 MHz masers are not detected. Consistent magnetic field strengths of approximately +3.5 mG are detected in the ground-state masers, in contrast with a possible -3.0 mG magnetic field previously detected at 13441 MHz. The variable 13441 MHz masers may be associated with an outflow.

Recent studies of methanol and ground-state OH masers at very high spectral resolution have shed new light on small-scale maser processes. The nearby source W3(OH), which contains numerous bright masers in several different transitions, provides an excellent laboratory for high spectral resolution techniques. We present a model of W3(OH) based on EVN observations of the rotationally-excited 6030 and 6035 MHz OH masers taken at 0.024 km/s spectral resolution. The 6.0 GHz masers are becoming brighter with time and show evidence for tangential proper motions. We confirm the existence of a region of magnetic field oriented toward the observer to the southeast and find another such region to the northeast in W3(OH), near the champagne flow. The 6.0 GHz masers trace the inner edge of a counterclockwise rotating torus feature. Masers at 6030 MHz are usually a factor of a few weaker than at 6035 MHz but trace the same material. Velocity gradients of nearby Zeeman components are much more closely correlated than in the ground state, likely due to the smaller spatial separation between Zeeman components. Hydroxyl maser peaks at very long baseline interferometric resolution appear to have structure on scales both smaller than that resolvable as well as on larger scales.

Recent observational and theoretical advances concerning astronomical masers in star forming regions are reviewed. Major masing species are considered individually and in combination. Key results are summarized with emphasis on present science and future prospects.

We report on a second epoch of VLBA observations of the 1665 and 1667 MHz OH masers in the massive star-forming region W75 N. We find evidence to confirm the existence of very strong (~40 mG) magnetic fields near source VLA 2. The masers near VLA 2 are dynamically distinct and include a very bright spot apparently moving at 50 km/s relative to those around VLA 1. This fast-moving spot may be an example of a rare class of OH masers seen in outflows in star-forming regions. Due to the variability of these masers and the rapidity of their motions, tracking these motions will require multiple observations over a significantly shorter time baseline than obtained here. Proper motions of the masers near VLA 1 are more suggestive of streaming along magnetized shocks rather than Keplerian rotation in a disk. The motions of the easternmost cluster of masers in W75 N (B) may be tracing slow expansion around an unseen exciting source.

Recent observational and theoretical advances have called into question traditional OH maser pumping models in evolved (OH/IR) stars. The detection of excited-state OH lines would provide additional constraints to discriminate amongst these theoretical models. In this Letter, we report on VLA observations of the 4750 MHz and 4765 MHz lines of OH toward 45 sources, mostly evolved stars. We detect 4765 MHz emission in the star forming regions Mon R2 and LDN 1084, but we do not detect excited-state emission in any evolved stars. The flux density and velocity of the 4765 MHz detection in Mon R2 suggests that a new flaring event has begun.

We present VLBA observations of 1667 MHz OH maser emission from the massive star formation region G5.89-0.39. The observations were phase referenced allowing the absolute positions of the masers to be obtained. The 1667 MHz masers have radial velocities that span ~50 km/s but show little evidence of tracing the bipolar molecular outflow, as has been claimed in previous studies. We identify 23 Zeeman pairs through comparison of masers in left and right circular polarization. Magnetic field strengths range from -2 mG to +2 mG, and an ordered reversal in magnetic field direction is observed toward the southern region of the UC HII region. We suggest that the velocity and magnetic field structure of the 1667 MHz masers can be explained in the context of a model in which the masers arise in a neutral shell just outside a rapidly exanding ionized shell.

Aims: We have observed the 6030 and 6035 MHz transitions of OH in high-mass star-forming regions to obtain magnetic field estimates in both maser emission and absorption. Methods: Observations were taken with the Effelsberg 100 m telescope. Results: Our observations are consistent with previous results, although we do detect a new 6030 MHz maser feature near -70 km/s in the vicinity of W3(OH). In absorption we obtain a possible estimate of -1.1 +/- 0.3 mG for the average line-of-sight component of the magnetic field in the absorbing OH gas in K3-50 and submilligauss upper limits for the line-of-sight field strength in DR 21 and W3. Conclusions: These results indicate that the magnetic field strength in the vicinity of OH masers is higher than that of the surrounding, non-masing material, which in turn suggests that the density of masing OH regions is higher than that of their surroundings.

A recent submillimeter line survey of Orion KL claimed detection of SiH. This paper reports on GBT observations of the 5.7 GHz Lambda-doubling transitions of SiH in Orion. Many recombination lines, including C164-delta, are seen, but SiH is not detected. The nondetection corresponds to an upper limit of 1.5 x 10^15 cm^-2 (4 sigma) for the beam-averaged column density of SiH. This suggests that the fractional abundance of SiH in the extended ridge is no more than twice that in the hot core.

We present VLBA observations of the ground-state hydroxyl masers in W3(OH) at 0.02 km s-1 spectral resolution. Over 250 masers are detected, including 56 Zeeman pairs. Lineshapes are predominantly Gaussian or combinations of several Gaussians, with normalized deviations typically of the same magnitude as in masers in other species. Typical FWHM maser linewidths are 0.15 to 0.38 km s-1 and are larger in the 1665 MHz transition than in the other three ground-state transitions. The satellite-line 1612 and 1720 MHz masers show no evidence of sigma^+/-2,3 components. The spatial positions of most masers are seen to vary across the line profile, with many spots showing clear, organized positional gradients. Equivalent line-of-sight velocity gradients in the plane of the sky typically range from 0.01 to 1 km s-1 AU-1 (i.e., positional gradients of 1 to 100 AU (km s-1)-1). Small velocity gradients in the 1667 MHz transition support theoretical predictions that 1667 MHz masers appear in regions with small velocity shifts along the amplification length. Deconvolved maser spot sizes appear to be larger in the line wings but do not support a spherical maser geometry.

We analyze full-polarization VLBA data of ground-state, main-line OH masers in 18 massive star-forming regions previously presented in a companion paper. The OH masers often arise in the shocked neutral gas surrounding ultracompact Hii regions. Magnetic fields as deduced from OH maser Zeeman splitting are highly ordered, both on the scale of a source as well as the maser clustering scale of ~10^15 cm. Results from our large sample show that this clustering scale appears to be universal to these masers. OH masers around ultracompact Hii regions live ~10^4 years and then turn off abruptly, rather than weakening gradually with time. These masers have a wide range of polarization properties. At one extreme (e.g., W75 N), pi-components are detected and the polarization position angles of maser spots show some organization. At the other extreme (e.g., W51 e1/e2), almost no linear polarization is detected and total polarization fractions can be substantially less than unity. A typical source has properties intermediate to these two extremes. In contrast to the well ordered magnetic field inferred from Zeeman splitting, there is generally no clear pattern in the distribution of polarization position angles. This can be explained if Faraday rotation in a typical OH maser source is large on a maser amplification length but small on a single (e-folding) gain length. Increasing or decreasing Faraday rotation by a factor of ~5 among different sources can explain the observed variation in polarization properties. We suggest that almost all pi-components acquire a signficant amount of circular polarization from low-gain stimulated emission of a sigma-component from OH appropriately shifted in velocity and lying along the propagation path.

We present full-polarization VLBA maps of the ground-state, main-line, 2 Pi 3/2, J = 3/2 OH masers in 18 Galactic massive star-forming regions. This is the first large polarization survey of interstellar hydroxyl masers at VLBI resolution. A total of 184 Zeeman pairs are identified, and the corresponding magnetic field strengths are indicated. We also present spectra of the NH3 emission or absorption in these star-forming regions. Analysis of these data will be presented in a companion paper.

We have observed six high-mass star-forming regions in the 2 Pi 3/2, J = 7/2 lines of OH using the GBT in order to investigate whether the magnetic field, and hence the density, measured in absorption differs from that implied by maser Zeeman splitting. We detect absorption in both the 13441 and 13434 MHz main lines in all six sources. Zeeman splitting in the F = 3-3 absorption line in W3(OH) implies a line-of-sight magnetic field strength of 3.0 +/- 0.3 mG. This is significantly less than full magnetic field strengths detected from OH maser Zeeman splitting, suggesting that OH maser regions may be denser than the non-masing OH material by a factor of several. Zeeman splitting is not detected in other sources, but we are able to place upper limits on B_parallel of 1.2 mG in G10.624-0.385 and 2.9 mG in K3-50. These results are consistent with a density enhancement of the masers, but other explanations for the lower magnetic field in absorption compared to maser emission are possible for these two sources. Absorption in one or both of the 13442 and 13433 MHz satellite lines is also seen in four sources. This is the very first detection of the 2 Pi 3/2, J = 7/2 satellite lines. Ratios of satellite-line to main-line absorption suggest enhancement of the satellite lines from local thermodynamic equilibrium values. Masers are seen in the F = 4-4 and 3-3 transitions of W3(OH) and the 4-4 transition of ON 1. A previously undetected 4-4 maser is seen near -44.85 km/s in W3(OH).

The DR 21 massive star forming region contains two cometary HII regions, aligned nearly perpendicular to each other on the sky. This offers a unique opportunity to discriminate among models of cometary HII regions. We present hydrogen recombination and ammonia line observations of DR 21 made with the Very Large Array. The velocity of the molecular gas, measured from ammonia emission and absorption lines, is constant to within +/- 1 km/s across the region. However, the radial velocity of the ionized material, measured from hydrogen recombination lines, differs by approximately 9 km/s between the "heads" of the two cometary HII regions and by up to approximately 7 km/s from that of the molecular gas. These findings indicate a supersonic velocity difference between the compact heads of the cometary regions and between each head and the ambient molecular material. This suggests that the observed cometary morphologies are created largely by the motion of wind-blowing, ionizing stars through the molecular cloud, as in a bow shock model.

We have identified and classified Zeeman pairs in the survey by Argon, Reid, & Menten of massive star-forming regions with 18 cm (2 Pi 3/2, J = 3/2) OH maser emission. We have found a total of more than 100 Zeeman pairs in more than 50 massive star-forming regions. The magnetic field deduced from the Zeeman splitting has allowed us to assign an overall line-of-sight magnetic field direction to many of the massive star-forming regions. Combining these data with other data sets obtained from OH Zeeman splitting, we have looked for correlations of magnetic field directions between star-forming regions scattered throughout the Galaxy. Our data do not support a uniform, Galactic-scale field direction, nor do we find any strong evidence of magnetic field correlations within spiral arms. However, our data suggest that in the Solar neighborhood the magnetic field outside the Solar circle is oriented clockwise as viewed from the North Galactic Pole, while inside the Solar circle it is oriented counterclockwise. This pattern, including the magnetic field reversal near the Sun, is in agreement with results obtained from pulsar rotation measures.