Using HST, we identify circumnuclear ($100$-$500$ pc scale) structures in nine new H$_2$O megamaser host galaxies to understand the flow of matter from kpc-scale galactic structures down to the supermassive black holes (SMBHs) at galactic centers. We double the sample analyzed in a similar way by Greene et al. (2013) and consider the properties of the combined sample of 18 sources. We find that disk-like structure is virtually ubiquitous when we can resolve $<200$ pc scales, in support of the notion that non-axisymmetries on these scales are a necessary condition for SMBH fueling. We perform an analysis of the orientation of our identified nuclear regions and compare it with the orientation of megamaser disks and the kpc-scale disks of the hosts. We find marginal evidence that the disk-like nuclear structures show increasing misalignment from the kpc-scale host galaxy disk as the scale of the structure decreases. In turn, we find that the orientation of both the $\sim100$ pc scale nuclear structures and their host galaxy large-scale disks is consistent with random with respect to the orientation of their respective megamaser disks.

We use new precision measurements of black hole masses from water megamaser disks to investigate scaling relations between macroscopic galaxy properties and supermassive black hole (BH) mass. The megamaser-derived BH masses span 10^6-10^8 M_sun, while all the galaxy properties that we examine (including stellar mass, central mass density, central velocity dispersion) lie within a narrow range. Thus, no galaxy property correlates tightly with M_BH in ~L* spiral galaxies. Of them all, stellar velocity dispersion provides the tightest relation, but at fixed sigma* the mean megamaser M_BH are offset by -0.6+/-0.1 dex relative to early-type galaxies. Spiral galaxies with non-maser dynamical BH masses do not show this offset. At low mass, we do not yet know the full distribution of BH mass at fixed galaxy property; the non-maser dynamical measurements may miss the low-mass end of the BH distribution due to inability to resolve the spheres of influence and/or megamasers may preferentially occur in lower-mass BHs.

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 the direct measurement of the Hubble constant, yielding the direct measurement of the angular-diameter distance to NGC 6264 using the H$_{2}$O megamaser technique. Our measurement is based on sensitive observations of the circumnuclear megamaser disk from four observations with the Very Long Baseline Array, the Green Bank Telescope and the Effelsberg Telescope. We also monitored the maser spectral profile for 2.3 years using the Green Bank Telescope to measure accelerations of maser lines by tracking their line-of-sight velocities as they change with time. The measured accelerations suggest that the systemic maser spots have a significantly wider radial distribution than in the archetypal megamaser in NGC 4258. We model the maser emission as arising from a circumnuclear disk with orbits dominated by the central black hole. The best fit of the data gives a Hubble constant of $H_{0} =$ 68$\pm$9 km s$^{-1}$ Mpc$^{-1}$, which corresponds to an angular-diameter distance of 144$\pm$19 Mpc. In addition, the fit also gives a mass of the central black hole of (3.09$\pm$0.42)$\times10^{7}$ $M_{\odot}$. The result demonstrates the feasibility of measuring distances to galaxies located well into the Hubble flow by using circumnuclear megamaser disks.

As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH.