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.

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.

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

The ALMA North America Prototype Antenna was awarded to the Smithsonian Astrophysical Observatory (SAO) in 2011. SAO and the Academia Sinica Institute of Astronomy & Astrophysics (ASIAA), SAO's main partner for this project, are working jointly to relocate the antenna to Greenland to carry out millimeter and submillimeter VLBI observations. This paper presents the work carried out on upgrading the antenna to enable operation in the Arctic climate by the GLT Team to make this challenging project possible, with an emphasis on the unexpected telescope components that had to be either redesigned or changed. Five-years of inactivity, with the antenna laying idle in the desert of New Mexico, coupled with the extreme weather conditions of the selected site in Greenland have it necessary to significantly refurbish the antenna. We found that many components did need to be replaced, such as the antenna support cone, the azimuth bearing, the carbon fiber quadrupod, the hexapod, the HVAC, the tiltmeters, the antenna electronic enclosures housing servo and other drive components, and the cables. We selected Vertex, the original antenna manufacturer, for the main design work, which is in progress. The next coming months will see the major antenna components and subsystems shipped to a site of the US East Coast for test-fitting the major antenna components, which have been retrofitted. The following step will be to ship the components to Greenland to carry out VLBI and single dish observations. Antenna reassembly at Summit Station should take place during the summer of 2018.

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.

We present the first constraint on Faraday rotation measure (RM) at submillimeter wavelengths for the nucleus of M 87. By fitting the polarization position angles ($\chi$) observed with the SMA at four independent frequencies around $\sim$230 GHz and interpreting the change in $\chi$ as a result of \emph{external} Faraday rotation associated with accretion flow, we determine the rotation measure of the M 87 core to be between $-$7.5$\times$10$^{5}$ and 3.4$\times$10$^{5}$ rad/m$^{2}$. Assuming a density profile of the accretion flow that follows a power-law distribution and a magnetic field that is ordered, radial, and has equipartition strength, the limit on the rotation measure constrains the mass accretion rate $\dot{M}$ to be below 9.2$\times$10$^{-4}$ M$_{\odot}$~yr$^{-1}$ at a distance of 21 Schwarzchild radii from the central black hole. This value is at least two orders of magnitude smaller than the Bondi accretion rate, suggesting significant suppression of the accretion rate in the inner region of the accretion flow. Consequently, our result disfavors the classical \emph{advection dominated accretion flow} (ADAF) and prefers the \emph{adiabatic inflow-outflow solution} (ADIOS) or \emph{convection-dominated accretion flow} (CDAF) for the hot accretion flow in M 87.

The GLT project is deploying a new submillimeter (submm) VLBI station in Greenland. Our primary scientific goal is to image a shadow of the supermassive black hole (SMBH) of six billion solar masses in M87 at the center of the Virgo cluster of galaxies. The expected SMBH shadow size of 40-50 $\mu$as requires superbly high angular resolution, suggesting that the submm VLBI would be the only way to obtain the shadow image. The Summit station in Greenland enables us to establish baselines longer than 9,000 km with ALMA in Chile and SMA in Hawaii as well as providing a unique $u$--$v$ coverage for imaging M87. Our VLBI network will achieve a superior angular resolution of about 20 $\mu$as at 350 GHz, corresponding to $\sim2.5$ times of the Schwarzschild radius of the supermassive black hole in M87. We have been monitoring the atmospheric opacity at 230 GHz since August. 2011; we have confirmed the value on site during the winter season is comparable to the ALMA site thanks to high altitude of 3,200 m and low temperature of $-50\degr$C. We will report current status and future plan of the GLT project towards our expected first light on 2015--2016.

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.

The next-generation space VLBI mission, VSOP-2, is expected to provide unprecedented spatial resolutions at 8.4, 22, and 43GHz. In this report, phase referencing with VSOP-2 is examined in detail based on a simulation tool called ARIS. The criterion for successful phase referencing was to keep the phase errors below one radian. Simulations with ARIS reveal that phase referencing achieves good performance at 8.4GHz, even under poor tropospheric conditions. At 22 and 43GHz, it is recommended to conduct phase referencing observations under good or typical tropospheric conditions. The satellite is required to have an attitude-switching capability with a one-minute or shorter cycle, and an orbit determination accuracy higher than about 10cm at apogee; the phase referencing calibrators are required to have a signal-to-noise ratio larger than four for a single scan. The probability to find a suitable phase referencing calibrator was estimated by using VLBI surveys. From the viewpoint of calibrator availability, VSOP-2 phase referencing at 8.4GHz is promising. However, the change of finding suitable calibrators at 22 and 43GHz is significantly reduced; it is important to conduct specific investigations for each target at those frequencies.

We report on a multi-frequency, multi-epoch campaign of Very Long Baseline Interferometry observations of the radio galaxy 1946+708 using the VLBA and a Global VLBI array. From these high-resolution observations we deduce the kinematic age of the radio source to be $\sim$4000 years, comparable with the ages of other Compact Symmetric Objects (CSOs). Ejections of pairs of jet components appears to take place on time scales of 10 years and these components in the jet travel outward at intrinsic velocities between 0.6 and 0.9 c. From the constraint that jet components cannot have intrinsic velocities faster than light, we derive H_0 > 57 km s^-1 Mpc^-1 from the fastest pair of components launched from the core. We provide strong evidence for the ejection of a new pair of components in ~1997. From the trajectories of the jet components we deduce that the jet is most likely to be helically confined, rather than purely ballistic in nature.

We present a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred to the quasar J1147-0724 located at an angular separation of 1.23$^{\circ}$. In observations spanning nearly two years, we have detected the maser spot at the LSR velocity of 34.7 km s$^{-1}$, for which we measured the annual parallax of 2.33$\pm$0.13 mas corresponding to a distance of 430$^{+25}_{-23}$ pc. This measurement has an accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate the stellar radius of S Crt and compare it with those of Mira variables.

Astrometric accuracy of complex modern VLBI arrays cannot be calculated analytically. We study the astrometric accuracy of phase-referenced VLBI observations for the VLBA, EVN and global VLBI array by simulating VLBI data for targets at declinations -25$^\circ$, 0$^\circ$, 25$^\circ$, 50$^\circ$, 75$^\circ$ and 85$^\circ$. The systematic error components considered in this study are calibrator position, station coordinate, Earth orientation and troposphere parameter uncertainties. We provide complete tables of the astrometric accuracies of these arrays for a source separation of 1$^\circ$ either along the right ascension axis or along the declination axis. Astrometric accuracy is 50microas at mid declination and is 300microas at low (-25$^\circ$) and high (85$^\circ$) declinations for the VLBA and EVN. In extending our simulations to source separations of 0.5$^\circ$ and 2$^\circ$, we establish the formula for the astrometric accuracy of the VLBA: Delta = (Delta_1$^\circ$-14)*d+ 14 (microas) where Delta_1$^\circ$ is the astrometric accuracy for a separation d=1$^\circ$ provided in our tables for various declinations and conditions of the wet troposphere. We argue that this formula is also valid for the astrometric accuracy of the EVN and global VLBI array.