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

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

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

We present a measurement of the Hubble constant made using geometric distance measurements to megamaser-hosting galaxies. We have applied an improved approach for fitting maser data and obtained better distance estimates for four galaxies previously published by the Megamaser Cosmology Project: UGC 3789, NGC 6264, NGC 6323, and NGC 5765b. Combining these updated distance measurements with those for the maser galaxies CGCG 074-064 and NGC 4258, and assuming a fixed velocity uncertainty of 250 km s$^{-1}$ associated with peculiar motions, we constrain the Hubble constant to be $H_0 = 73.9 \pm 3.0$ km s$^{-1}$ Mpc$^{-1}$ independent of distance ladders and the cosmic microwave background. This best value relies solely on maser-based distance and velocity measurements, and it does not use any peculiar velocity corrections. Different approaches for correcting peculiar velocities do not modify $H_0$ by more than ${\pm}1{\sigma}$, with the full range of best-fit Hubble constant values spanning 71.8-76.9 km s$^{-1}$ Mpc$^{-1}$. We corroborate prior indications that the local value of $H_0$ exceeds the early-Universe value, with a confidence level varying from 95-99% for different treatments of the peculiar velocities.

We present a new investigation of the dependence of H2O maser detection rates and properties on the mid-IR AGN luminosity, L_AGN, and the obscuring column density, N_H, based on mid-IR and hard X-ray photometry. Based on spectral energy distribution fitting that allows for decomposition of the black hole accretion and star-formation components in the mid-infrared, we show that the megamaser (disk maser) detection rate increases sharply for galaxies with 12 micron AGN luminosity L^{AGN}_{12 micron} greater than 10^42 erg/s, from ~<3%(~<2%) to ~12%(~5%). By using the ratio of the observed X-ray to mid-IR AGN luminosity as an indicator of N_H, we also find that the megamaser (disk maser) detection rates are boosted to 15%(7%) and 20%(9%) for galaxies with N_H >= 10^23 cm^{-2} and N_H >= 10^{24} cm^{-2}, respectively. Combining these column density cuts with a constraint for high L^{AGN}_{12 micron} (>=10^42 erg/s) predicts further increases in the megamaser (disk maser) detection rates to 19%(8%) and 27%(14%), revealing unprecedented potential boosts of the megamaser and disk maser detection rates by a factor of 7-15 relative to the current rates, depending on the chosen sample selection criteria. A noteworthy aspect of these new predictions is that the completeness rates are only compromised mildly, with the rates remaining at the level of ~95%(~50%) for sources with N_H >= 10^{23} cm^{-2} (N_H >= 10^24 cm^-2). Applying these selection methods to the current X-ray AGN surveys predicts the detection of >~15 new megamaser disks.

As part of the survey component of the Megamaser Cosmology Project, we have discovered a disk megamaser system in the galaxy CGCG 074-064. Using the GBT and the VLA, we have obtained spectral monitoring observations of this maser system at a monthly cadence over the course of two years. We find that the systemic maser features display line-of-sight accelerations of ~4.4 km s$^{-1}$ yr$^{-1}$ that are nearly constant with velocity, while the high-velocity maser features show accelerations that are consistent with zero. We have also used the HSA to make a high-sensitivity VLBI map of the maser system in CGCG 074-064, which reveals that the masers reside in a thin, edge-on disk with a diameter of ~1.5 mas (0.6 pc). Fitting a three-dimensional warped disk model to the data, we measure a black hole mass of $2.42^{+0.22}_{-0.20} \times 10^7$ M$_{\odot}$ and a geometric distance to the system of $87.6^{+7.9}_{-7.2}$ Mpc. Assuming a CMB-frame recession velocity of $7308 \pm 150$ km s$^{-1}$, we constrain the Hubble constant to $H_0 = 81.0^{+7.4}_{-6.9}$ (stat.) $\pm 1.4$ (sys.) km s$^{-1}$ Mpc$^{-1}$.

NGC 4258 is a critical galaxy for establishing the extragalactic distance scale and estimating the Hubble constant (Ho). Water masers in the nucleus of the galaxy orbit about its supermassive black hole, and very long baseline interferometric observations of their positions, velocities, and accelerations can be modeled to give a geometric estimate of the angular-diameter distance to the galaxy. We have improved the technique to obtain model parameter values, reducing both statistical and systematic uncertainties compared to previous analyses. We find the distance to NGC 4258 to be 7.576 +/- 0.082 (stat.) +/- 0.076 (sys.) Mpc. Using this as the sole source of calibration of the Cepheid-SN Ia distance ladder results in Ho = 72.0 +/- 1.9 km/s/Mpc, and in concert with geometric distances from Milky Way parallaxes and detached eclipsing binaries in the LMC we find Ho = 73.5 +/- 1.4 km/s/Mpc. The improved distance to NGC 4258 also provides a new calibration of the tip of the red giant branch of M_{F814W} = -4.01 +/- 0.04$ mag, with reduced systematic errors for the determination of Ho compared to the LMC-based calibration, because it is measured on the same Hubble Space Telescope photometric system and through similarly low extinction as SN Ia host halos. The result is Ho = 71.1 +/- 1.9 km/s/Mpc, in good agreement with the result from the Cepheid route, and there is no difference in Ho when using the same calibration from NGC 4258 and the same SN Ia Hubble diagram intercept to start and end both distance ladders.

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.

We present VLBI images and kinematics of water maser emission in three active galaxies: NGC 5728, Mrk 1, and IRAS 08452-0011. IRAS 08452-0011 is a triple-peaked H2O megamaser, consistent with a Keplerian rotating disk, indicating a black hole mass of (3.2+/-0.2)x10^7 M_sun. NGC 5728 and Mrk 1 display double-peaked spectra and VLBI imaging reveal complicated gas kinematics, which do not allow for a robust determination of black hole mass. The two systems are either gas disks perturbed by AGN winds or part of outflows. We find that disturbed morphology and kinematics are a ubiquitous feature of all double-peaked maser systems, implying that these maser sources may reside in environments where AGN winds are prominent at ~1 pc scale and have significant impact on the masing gas. Such AGN tend to have black hole mass M_BH < 8x10^6 M_sun and Eddington ratios lambda_Edd >~ 0.1, while the triple-peaked megamasers show an opposite trend.

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

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.

In combination with observations of the Cosmic Microwave Background, a measurement of the Hubble Constant provides a direct test of the standard $\Lambda$CDM cosmological model and a powerful constraint on the equation of state of dark energy. Observations of circumnuclear water vapor megamasers at 22 GHz in nearby active galaxies can be used to measure their distances, geometrically, and thereby provide a direct, one step measurement of the Hubble Constant. The measurement is independent of distance ladders and standard candles. With present-day instrumentation, the Megamaser Cosmology Project is expected to reach a $\sim$4% measurement using the megamaser technique, based on distances to fewer than 10 megamaser systems. A goal of the observational cosmology community is to reach a percent-level measurement of H$_0$ from several independent astrophysical measurements to minimize the systematics. The ngVLA will provide the sensitivity at 22 GHz required to discover and measure the additional megamaser disks that will enable an H$_0$ measurement at the $\sim$1% level.

The Next-Generation Very Large Array (ngVLA) has the potential to be a workhorse for the discovery and study of paired supermassive black holes either at large separations (dual) or in tightly bound systems (binary). In this chapter, we outline the science case for the study of these supermassive pairs, and summarize discovery methods that can be used at radio wavelengths to discover them: including morphological, spectral, and time-domain identifications. One critical aspect of this work is that multi-messenger binary black hole studies may be possible with the ngVLA when combined with gravitational-wave searches using pulsar timing array techniques. However, long-baseline interferometery (>>1000 km) will make this possibility more likely by expanding the redshift range at which radio emission arising from two separate black holes may be resolved and studied.

H$_2$O megamasers residing in the accretion disks of active galactic nuclei (AGN) exhibit Keplerian rotation about the central supermassive black hole (SMBH). Such disk maser systems are excellent tools for diagnosing the kinematic status of the SMBH, and they currently provide the only direct and unambiguous measure of SMBH velocities outside of the Milky Way. We have measured the galaxy recession velocities for a sample of 10 maser disk systems using a combination of spatially resolved HI disk modeling, spatially integrated HI profile fitting, and optical spectral line and continuum fitting. In comparing the SMBH velocities to those of their host galaxies, we find two (out of 10) systems -- J0437+2456 and NGC 6264 -- for which the SMBH and galaxy velocities show a statistically significant ($>$3$\sigma$) difference. For NGC 6264 the apparent velocity offset can likely be explained by ionized gas motion within the host galaxy (e.g., from AGN-driven shocks). The velocity measurements for J0437+2456, however, imply a SMBH peculiar velocity of $69.6 \pm 12.7$ km s$^{-1}$ (5.5$\sigma$). We thus consider J0437+2456 to be a promising candidate for hosting either a recoiling or binary SMBH, though additional observations are necessary to exclude the possibility of a systematic offset between the galactic recession velocity and that measured using the optical spectrum.

Water megamasers from circumnuclear disks in galaxy centers provide the most accurate measurements of supermassive black hole masses and uniquely probe the sub-parsec accretion processes. At the same time, these systems offer independent crucial constraints of the Hubble Constant in the nearby universe, and thus, the arguably best single constraint on the nature of dark energy. The chances of finding these golden standards are however abysmally low, at an overall =< 3% for any level of water maser emission detected at 22 GHz, and =< 1% for those exhibiting the disk-like configuration. We provide here a thorough summary of the current state of the detection of water megamaser disks, along with a novel investigation of the likelihood of increasing their detection rates based on a multivariate parameter analysis of optical and mid-infrared photometric properties of the largest database of galaxies surveyed for the 22 GHz emission. We find that galaxies with water megamaser emission tend to associate with strong emission in all WISE mid-infrared wavelengths, with the strongest enhancement in the W4 band, at 22 micron, as well as with previously proposed and newly found indicators of AGN strength in the mid-infrared, such as red W1-W2 and W1-W4 colors, and the integrated mid-infrared luminosity of the host galaxy. These trends offer a potential boost of the megamaser detection rates to 6 - 15%, or a factor of 2-8 relative to the current rates, depending on the chosen sample selection criteria, while fostering real chances for discovering >= 20 new megamaser disks.

We present high resolution (sub-mas) VLBI maps of nuclear H2O megamasers for seven galaxies. In UGC6093, the well-aligned systemic masers and high-velocity masers originate in an edge-on, flat disk and we determine the mass of the central SMBH to be M_SMBH = 2.58*10^7Msun(+-7%). For J1346+5228, the distribution of masers is consistent with a disk, but the faint high-velocity masers are only marginally detected, and we constrain the mass of the SMBH to be in the range 1.5-2.0*10^7Msun. The origin of the masers in Mrk1210 is less clear, as the systemic and high-velocity masers are misaligned and show a disorganized velocity structure. We present one possible model in which the masers originate in a tilted, warped disk, but we do not rule out the possibility of other explanations including outflow masers. In NGC6926, we detect a set of redshifted masers, clustered within a pc of each other, and a single blueshifted maser about 4.4pc away, an offset that would be unusually large for a maser disk system. Nevertheless, if it is a disk system, we estimate the enclosed mass to be M_SMBH<4.8*10^7 Msun . For NGC5793, we detect redshifted masers spaced about 1.4pc from a clustered set of blueshifted features. The orientation of the structure supports a disk scenario as suggested by Hagiwara et al.(2001). We estimate the enclosed mass to be M SMBH<1.3*10^7 Msun. For NGC2824 and J0350-0127, the masers may be associated with pc or sub-pc scale jets or outflows.

As part of the Megamaser Cosmology Project (MCP), we present VLBI maps of nuclear water masers toward five galaxies. The masers originate in sub-parsec circumnuclear disks. For three of the galaxies, we fit Keplerian rotation curves to estimate their supermassive black hole (SMBH) masses, and determine (2.9 $\pm$ 0.3) $\times~10^{6}M_\odot$ for J0437+2456, (1.7 $\pm$ 0.1) $\times~10^{7}M_\odot$ for ESO 558$-$G009, and (1.1 $\pm$ 0.2) $\times~10^{7}M_\odot$ for NGC 5495. In the other two galaxies, Mrk 1029 and NGC 1320, the geometry and dynamics are more complicated and preclude robust black hole mass estimates. Including our new results, we compiled a list of 15 VLBI-confirmed disk maser galaxies with robust SMBH mass measurements. With this sample, we confirm the empirical relation of $R_{out} \propto 0.3 M_{SMBH}$ reported in Wardle & Yusef-Zadeh (2012). We also find a tentative correlation between maser disk outer radii and WISE luminosity. We find no correlations of maser disk size with X-ray 2-10 keV luminosity or [O III] luminosity.

We present 321 GHz observations of five AGN from ALMA Cycle 0 archival data: NGC 5793, NGC 1068, NGC 1386, NGC 4945, and the Circinus galaxy. Submillimeter maser emission is detected for the first time towards NGC 4945, and we present a new analysis of the submillimeter maser system in Circinus. None of the other three galaxies show maser emission, though we have detected and imaged the continuum from every galaxy. Both NGC 4945 and Circinus are known to host strong ($\gtrsim 10$ Jy) 22 GHz megamaser emission, and VLBI observations have shown that the masers reside in the innermost $\sim 1$ parsec of the galaxies. The peak flux densities of the 321 GHz masers in both systems are substantially weaker (by a factor of $\sim$100) than what is observed at 22 GHz, though the corresponding isotropic luminosities are more closely matched (within a factor of $\sim$10) between the two transitions. We compare the submillimeter spectra presented here to the known 22 GHz spectra in both galaxies, and we argue that while both transitions originate from the gaseous environment near the AGN, not all sites are necessarily in common. In Circinus, the spectral structure of the 321 GHz masers indicates that they may trace the accretion disk at radii interior to the 22 GHz masers. The continuum emission in NGC 4945 and NGC 5793 shows a spatial distribution indicative of an origin in the galactic disks (likely thermal dust emission), while for the other three galaxies the emission is centrally concentrated and likely originates from the nucleus.

As part of the Megamaser Cosmology Project (MCP), here we present a new geometric distance measurement to the megamaser galaxy NGC 5765b. Through a series of VLBI observations, we have confirmed the water masers trace a thin, sub-parsec Keplerian disk around the nucleus, implying an enclosed mass of 4.55 $\pm$ 0.40 $\times~10^{7}M_\odot$. Meanwhile, from single dish monitoring of the maser spectra over two years, we measured the secular drifts of maser features near the systemic velocity of the galaxy with rates between 0.5 and 1.2 km s$^{-1}$ yr$^{-1}$. Fitting a warped, thin disk model to these measurements, we determine a Hubble Constant $H_{0}$ of 66.0 $\pm$ 6.0 km s$^{-1}$ Mpc$^{-1}$ with the angular-diameter distance to NGC 5765b of 126.3 $\pm$ 11.6 Mpc. Apart from the distance measurement, we also investigate some physical properties related to the maser disk in NGC 5765b. The high-velocity features are spatially distributed into several clumps, which may indicate the existence of a spiral density wave associated with the accretion disk. For the red-shifted features, the envelope defined by the peak maser intensities increases with radius. The profile of the systemic masers in NGC 5765b is smooth and shows almost no structural changes over the two years of monitoring time, which differs from the more variable case of NGC 4258.

We use single-dish radio spectra of known 22 GHz H$_2$O megamasers, primarily gathered from the large dataset observed by the Megamaser Cosmology Project, to identify Keplerian accretion disks and to investigate several aspects of the disk physics. We test a mechanism for maser excitation proposed by Maoz & McKee (1998), whereby population inversion arises in gas behind spiral shocks traveling through the disk. Though the flux of redshifted features is larger on average than that of blueshifted features, in support of the model, the high-velocity features show none of the predicted systematic velocity drifts. We find rapid intra-day variability in the maser spectrum of ESO 558-G009 that is likely the result of interstellar scintillation, for which we favor a nearby ($D \approx 70$ pc) scattering screen. In a search for reverberation in six well-sampled sources, we find that any radially-propagating signal must be contributing $\lesssim$10% of the total variability. We also set limits on the magnetic field strengths in seven sources, using strong flaring events to check for the presence of Zeeman splitting. These limits are typically 200--300 mG ($1\sigma$), but our most stringent limits reach down to 73 mG for the galaxy NGC 1194.

We present imaging and spectroscopic observations with HST and VLT of the ring of SN 1987A from 1994 to 2014. After an almost exponential increase of the shocked emission from the hotspots up to day ~8,000 (~2009), both this and the unshocked emission are now fading. From the radial positions of the hotspots we see an acceleration of these up to 500-1000 km/s, consistent with the highest spectroscopic shock velocities from the radiative shocks. In the most recent observations (2013 and 2014), we find several new hotspots outside the inner ring, excited by either X-rays from the shocks or by direct shock interaction. All of these observations indicate that the interaction with the supernova ejecta is now gradually dissolving the hotspots. We predict, based on the observed decay, that the inner ring will be destroyed by ~2025.

We present observations of the H2O megamasers in the accretion disk of NGC 6323. By combining interferometric and spectral monitoring data, we estimate H$_{0} = 73^{+26}_{-22}$ km/s/Mpc, where the low strength of the systemic masers (<15 mJy) limits the accuracy of this estimate. The methods developed here for dealing with weak maser emission provide guidance for observations of similar sources, until significant increases in radio telescope sensitivity, such as anticipated from the next generation Very Large Array, are realized.