Modern radio interferometers enable observations of spectral lines with unprecedented spatial resolution and sensitivity. In spite of these technical advances, many lines of interest are still at best weakly detected and therefore necessitate detection and analysis techniques specialized for the low signal-to-noise ratio (SNR) regime. Matched filters can leverage knowledge of the source structure and kinematics to increase sensitivity of spectral line observations. Application of the filter in the native Fourier domain improves SNR while simultaneously avoiding the computational cost and ambiguities associated with imaging, making matched filtering a fast and robust method for weak spectral line detection. We demonstrate how an approximate matched filter can be constructed from a previously observed line or from a model of the source, and we show how this filter can be used to robustly infer a detection significance for weak spectral lines. When applied to ALMA Cycle 2 observations of CH3OH in the protoplanetary disk around TW Hya, the technique yields a ~53% SNR boost over aperture-based spectral extraction methods, and we show that an even higher boost will be achieved for observations at higher spatial resolution. A Python-based open-source implementation of this technique is available under the MIT license at https://github.com/AstroChem/VISIBLE.

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.

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.

Turbulence can transport angular momentum in protoplanetary disks and influence the growth and evolution of planets. With spatially and spectrally resolved molecular emission line measurements provided by (sub)millimeter interferometric observations, it is possible to directly measure non-thermal motions in the disk gas that can be attributed to this turbulence. We report a new constraint on the turbulence in the disk around HD 163296, a nearby young A star, determined from ALMA Science Verification observations of four CO emission lines (the CO(3-2), CO(2-1), 13CO(2-1), and C18O(2-1) transitions). The different optical depths for these lines permit probes of non-thermal line-widths at a range of physical conditions (temperature and density) and depths into the disk interior. We derive stringent limits on the non-thermal motions in the upper layers of the outer disk such that any contribution to the line-widths from turbulence is <3% of the local sound speed. These limits are approximately an order of magnitude lower than theoretical predictions for full-blown MHD turbulence driven by the magneto-rotational instability, potentially suggesting that this mechanism is less efficient in the outer (R>30AU) disk than has been previously considered.

We present results from a spectroscopic Spitzer and Herschel mid-to-far-infrared study of the circumbinary disk orbiting the evolved (age ~12-23 Myr) close binary T Tauri system V4046 Sgr. Spitzer IRS spectra show emission lines of [Ne II], H_2 S(1), CO_2 and HCN, while Herschel PACS and SPIRE spectra reveal emission from [O I], OH, and tentative detections of H_2O and high-J transitions of CO. We measure [Ne III]/[Ne II] < 0.13, which is comparable to other X-ray/EUV luminous T Tauri stars that lack jets. We use the H_2 S(1) line luminosity to estimate the gas mass in the relatively warm surface layers of the inner disk. The presence of [O I] emission suggests that CO, H_2O, and/or OH is being photodissociated, and the lack of [C I] emission suggests any excess C may be locked up in HCN, CN and other organic molecules. Modeling of silicate dust grain emission features in the mid-infrared indicates that the inner disk is composed mainly of large (r~5 um) amorphous pyroxene and olivine grains (~86% by mass) with a relatively large proportion of crystalline silicates. These results are consistent with other lines of evidence indicating that planet building is ongoing in regions of the disk within ~30 AU of the central, close binary.

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

Protoplanetary "transition" disks have large, mass-depleted central cavities, yet also deliver gas onto their host stars at rates comparable to disks without holes. The paradox of simultaneous transparency and accretion can be explained if gas flows inward at much higher radial speeds inside the cavity than outside the cavity, since surface density (and by extension optical depth) varies inversely with inflow velocity at fixed accretion rate. Radial speeds within the cavity might even have to approach free-fall values to explain the huge surface density contrasts inferred for transition disks. We identify observational diagnostics of fast radial inflow in channel maps made in optically thick spectral lines. Signatures include (1) twisted isophotes in maps made at low systemic velocities and (2) rotation of structures observed between maps made in high-velocity line wings. As a test case, we apply our new diagnostic tools to archival ALMA data on the transition disk HD 142527, and uncover evidence for free-fall radial velocities inside its cavity. Although the observed kinematics are also consistent with a disk warp, the radial inflow scenario is preferred because it predicts low surface densities that appear consistent with recent observations of optically thin CO isotopologues in this disk. How material in the disk cavity sheds its angular momentum wholesale to fall freely onto the star is an unsolved problem; gravitational torques exerted by giant planets or brown dwarfs are briefly discussed as a candidate mechanism.

We present sensitive, sub-arcsecond resolution Submillimeter Array observations of the protoplanetary disk around the nearby, pre-main sequence spectroscopic binary V4046 Sgr. We report for the first time a large inner hole (r=29 AU) spatially resolved in the 1.3 mm continuum emission and study the structure of this disk using radiative transfer calculations to model the spectral energy distribution (SED), continuum visibilities, and spectral line emission of CO and its main isotopologues. Our modeling scheme demonstrates that the majority of the dust mass is distributed in a narrow ring (centered at 37 AU with a FWHM of 16 AU) that is ~5 times more compact than the gas disk. This structure implies that the dust-to-gas mass ratio has a strong spatial variation, ranging from a value much larger than typical of the interstellar medium (ISM) at the ring to much smaller than that of the ISM at larger disk radii. We suggest that these basic structural features are potentially observational signatures of the accumulation of solids at a local gas pressure maximum. These models also require a substantial population of ~micron-sized grains inside the central disk cavity. We suggest that this structure is likely the result of dynamical interactions with a low-mass companion, although photoevaporation may also play a secondary role.

Theory predicts that giant planets and low mass stellar companions shape circumstellar disks by opening annular gaps in the gas and dust spatial distribution. For more than a decade it has been debated whether this is the dominant process that leads to the formation of transitional disks. In this paper, we present millimeter-wave interferometric observations of the transitional disk around the young intermediate mass star LkHa330. These observations reveal a lopsided ring in the 1.3 mm dust thermal emission characterized by a radius of about 100 AU and an azimuthal intensity variation of a factor of 2. By comparing the observations with a Gaussian parametric model, we find that the observed asymmetry is consistent with a circular arc, that extends azimuthally by about 90 deg and emits about 1/3 of the total continuum flux at 1.3 mm. Hydrodynamic simulations show that this structure is similar to the azimuthal asymmetries in the disk surface density that might be produced by the dynamical interaction with unseen low mass companions orbiting within 70 AU from the central star. We argue that such asymmetries might lead to azimuthal variations in the millimeter-wave dust opacity and in the dust temperature, which will also affect the millimeter-wave continuum emission. Alternative explanations for the observed asymmetry that do not require the presence of companions cannot be ruled out with the existing data. Further observations of both the dust and molecular gas emission are required to derive firm conclusions on the origin of the asymmetry observed in the LkHa330 disk.

We analyze sensitive, sub-arcsecond resolution ALMA Science Verification observations of CO emission lines in the protoplanetary disk hosted by the young, isolated Ae star HD 163296. The observed spatial morphology of the 12CO J=3-2 emission line is asymmetric across the major axis of the disk; the 12CO J=2-1 line features a much less pronounced, but similar, asymmetry. The J=2-1 emission from 12CO and its main isotopologues have no resolved spatial asymmetry. We associate this behavior as the direct signature of a vertical temperature gradient and layered molecular structure in the disk. This is demonstrated using both toy models and more sophisticated calculations assuming non-local thermodynamic equilibrium (non LTE) conditions. A model disk structure is developed to reproduce both the distinctive spatial morphology of the 12CO J=3-2 line as well as the J=2-1 emission from the CO isotopologues assuming relative abundances consistent with the interstellar medium. This model disk structure has tau=1 emitting surfaces for the 12CO emission lines that make an angle of about 15 degrees with respect to the disk midplane. Furthermore, we show that the spatial and spectral sensitivity of these data can distinguish between models that have sub-Keplerian gas velocities due to the vertical extent of the disk and its associated radial pressure gradient (a fractional difference in the bulk gas velocity field of approximately greater than 5%).

We present a substantial extension of the mm-wave continuum photometry catalog for Taurus circumstellar dust disks. Combining new Submillimeter Array data with measurements in the literature, we construct a mm-wave luminosity distribution for Class II disks that is statistically complete for stellar hosts with spectral types earlier than M8.5 and has a (3-sigma) depth of ~3 mJy. The resulting census eliminates a longstanding bias against disks with late-type hosts, and thereby reveals a strong correlation between L_mm and the host spectral type. We confirm that this corresponds to a statistically robust relationship between the masses of dust disks and the stars that host them. A Bayesian regression technique is used to characterize these relationships: the results indicate a typical 1.3 mm flux density of 25 mJy for solar mass hosts and a power-law scaling L_mm \propto M_star^1.5-2.0. We suggest that a reasonable treatment of dust temperature in the conversion from L_mm to M_disk favors an inherently linear M_disk \propto M_star scaling, with a typical disk-to-star mass ratio of $\sim$0.2--0.6%. The RMS dispersion around this regression is 0.7 dex, suggesting that the combined effects of diverse evolutionary states, dust opacities, and temperatures in these disks imprint a FWHM range of a factor of 40 on the inferred M_disk (or L_mm) at any given host mass. We argue that this relationship between M_disk and M_star likely represents the origin of the inferred correlation between giant planet frequency and host star mass in the exoplanet population, and provides some basic support for the core accretion model for planet formation. Moreover, we caution that selection bias must be considered in comparative studies of disk evolution, and illustrate that fact with statistical comparisons of L_mm between Taurus and other clusters (abridged).

We report the first detection of c-C3H2 in a circumstellar disk. The c-C3H2 J=6-5 line (217.882 GHz) is detected and imaged through Atacama Large Millimeter Array (ALMA) Science Verification observations toward the disk around the Herbig Ae star HD 163296 at 0.8" resolution. The emission is consistent with that arising from a Keplerian rotating disk. Two additional c-C3H2 transitions are also tentatively detected, bolstering the identification of this species, but with insufficient signal-to-noise ratio to constrain the spatial distribution. Using a previously developed model for the physical structure of this disk, we fit a radial power-law distribution model to the c-C3H2 6-5 emission and find that c-C3H2 is present in a ring structure from an inner radius of about 30 AU to an outer radius of about 165 AU. The column density is estimated to be 1e12-1e13 cm-2. The clear detection and intriguing ring structure suggest that c-C3H2 has the potential to become a useful probe of radiation penetration in disks.

We present 1.3 millimeter ALMA Cycle 0 observations of the edge-on debris disk around the nearby, ~10 Myr-old, M-type star AU Mic. These observations obtain 0.6 arcsec (6 AU) resolution and reveal two distinct emission components: (1) the previously known dust belt that extends to a radius of 40 AU, and (2) a newly recognized central peak that remains unresolved. The cold dust belt of mass about 1 lunar mass is resolved in the radial direction with a rising emission profile that peaks sharply at the location of the outer edge of the "birth ring" of planetesimals hypothesized to explain the midplane scattered light gradients. No significant asymmetries are discerned in the structure or position of this dust belt. The central peak identified in the ALMA image is ~6 times brighter than the stellar photosphere, which indicates an additional emission process in the inner regions of the system. Emission from a stellar corona or activity may contribute, but the observations show no signs of temporal variations characteristic of radio-wave flares. We suggest that this central component may be dominated by dust emission from an inner planetesimal belt of mass about 0.01 lunar mass, consistent with a lack of emission shortward of 25 microns and a location <3 AU from the star. Future millimeter observations can test this assertion, as an inner dust belt should be readily separated from the central star at higher angular resolution.

We present sensitive, arcsecond-resolution Submillimeter Array observations of the 12CO J=2-1 line emission from the circumstellar disk orbiting the double-lined spectroscopic binary star V4046 Sgr. Based on a simple model of the disk structure, we use a novel Monte Carlo Markov Chain technique to extract the Keplerian velocity field of the disk from these data and estimate the total mass of the central binary. Assuming the distance inferred from kinematic parallax measurements in the literature (d is approximately 73 pc), we determine a total stellar mass M_star = 1.75^{+0.09}_{-0.06} solar masses and a disk inclination i_d = 33.5^{+0.7}_{-1.4} degrees from face-on. These measurements are in excellent agreement with independent dynamical constraints made from multi-epoch monitoring of the stellar radial velocities, confirming the absolute accuracy of this precise (~ few percent uncertainties) disk-based method for estimating stellar masses and reaffirming previous assertions that the disk and binary orbital planes are well aligned (with |i_d - i_star| \approx 0.1\pm1 degree). Using these results as a reference, we demonstrate that various pre-main sequence evolution models make consistent and accurate predictions for the masses of the individual components of the binary, and uniformly imply an advanced age of ~5-30 Myr. Taken together, these results verify that V4046 Sgr is one of the precious few nearby and relatively evolved pre-main sequence systems that still hosts a gas-rich accretion disk.

The processes that form transition disks - disks with depleted inner regions - are not well understood; possible scenarios include planet formation, grain growth and photoevaporation. Disks with spatially resolved dust holes are rare, but, in general, even less is known about the gas structure. The disk surrounding A0 star Oph IRS 48 in the nearby Rho Ophiuchus region has a 30 AU radius hole previously detected in the 18.7 micron dust continuum and in warm CO in the 5 micron fundamental ro-vibrational band. We present here Submillimeter Array 880 micron continuum imaging resolving an inner hole. However, the radius of the hole in the millimeter dust is only 13 AU, significantly smaller than measured at other wavelengths. The nesting structure of the disk is counter-intuitive, with increasingly large radii rings of emission seen in the millimeter dust (12.9 +1.7/-3.4 AU), 5 micron CO (30 AU) and 18.7 micron dust (peaking at 55 AU). We discuss possible explanations for this structure, including self-shadowing that cools the disk surface layers, photodissociation of CO, and photoevaporation. However, understanding this unusual disk within the stringent multi-wavelength spatial constraints will require further observations to search for cold atomic and molecular gas.

We present a detailed analysis of the spatially and spectrally resolved 12CO J=2-1 and J=3-2 emission lines from the TW Hya circumstellar disk, based on science verification data from the Atacama Large Millimeter/Submillimeter Array (ALMA). These lines exhibit substantial emission in their high-velocity wings (with projected velocities out to 2.1 km/s, corresponding to intrinsic orbital velocities >20 km/s) that trace molecular gas as close as 2 AU from the central star. However, we are not able to reproduce the intensity of these wings and the general spatio-kinematic pattern of the lines with simple models for the disk structure and kinematics. Using three-dimensional non-local thermodynamic equilibrium molecular excitation and radiative transfer calculations, we construct some alternative models that successfully account for these features by modifying either (1) the temperature structure of the inner disk (inside the dust-depleted disk cavity; r < 4 AU); (2) the intrinsic (Keplerian) disk velocity field; or (3) the distribution of disk inclination angles (a warp). The latter approach is particularly compelling because a representative warped disk model qualitatively reproduces the observed azimuthal modulation of optical light scattered off the disk surface. In any model scenario, the ALMA data clearly require a substantial molecular gas reservoir located inside the region where dust optical depths are known to be substantially diminished in the TW Hya disk, in agreement with previous studies based on infrared spectroscopy. The results from these updated model prescriptions are discussed in terms of their potential physical origins, which might include dynamical perturbations from a low-mass companion with an orbital separation of a few AU.

We present high resolution (0.3" = 16 AU), high signal-to-noise ratio Submillimeter Array observations of the 870 microns (345 GHz) continuum and CO J=3--2 line emission from the protoplanetary disk around TW Hya. Using continuum and line radiative transfer calculations, those data and the multiwavelength spectral energy distribution are analyzed together in the context of simple two-dimensional parametric disk structure models. Under the assumptions of a radially invariant dust population and (vertically integrated) gas-to-dust mass ratio, we are unable to simultaneously reproduce the CO and dust observations with model structures that employ either a single, distinct outer boundary or a smooth (exponential) taper at large radii. Instead, we find that the distribution of millimeter-sized dust grains in the TW Hya disk has a relatively sharp edge near 60 AU, contrary to the CO emission (and optical/infrared scattered light) that extends to a much larger radius of at least 215 AU. We discuss some possible explanations for the observed radial distribution of millimeter-sized dust grains and the apparent CO-dust size discrepancy, and suggest that they may be hallmarks of substructure in the dust disk or natural signatures of the growth and radial drift of solids that might be expected for disks around older pre-main sequence stars like TW Hya.

We present 870 micron observations of dust continuum emission from the LkCa 15 protoplanetary disk at high angular resolution (with a characteristic scale of 0.25" = 35 AU), obtained with the IRAM Plateau de Bure interferometer and supplemented by slightly lower resolution observations from the Submillimeter Array. We fit these data with simple morphological models to characterize the spectacular ring-like emission structure of this disk. Our analysis indicates that a small amount of 870 micron dust emission (~5 mJy) originates inside a large (40-50 AU radius) low optical depth cavity. This result can be interpreted either in the context of an abrupt decrease by a factor of ~5 in the radial distribution of millimeter-sized dust grains or as indirect evidence for a gap in the disk, in agreement with previous inferences from the unresolved infrared spectrum and scattered light images. A preliminary model focused on the latter possibility suggests the presence of a low-mass (planetary) companion, having properties commensurate with those inferred from the recent discovery of LkCa 15b.

We describe a novel approach to accelerating Monte Carlo Markov Chains. Our focus is cosmological parameter estimation, but the algorithm is applicable to any problem for which the likelihood surface is a smooth function of the free parameters and computationally expensive to evaluate. We generate a high-order interpolating polynomial for the log-likelihood using the first points gathered by the Markov chains as a training set. This polynomial then accurately computes the majority of the likelihoods needed in the latter parts of the chains. We implement a simple version of this algorithm as a patch (InterpMC) to CosmoMC and show that it accelerates parameter estimatation by a factor of between two and four for well-converged chains. The current code is primarily intended as a "proof of concept", and we argue that there is considerable room for further performance gains. Unlike other approaches to accelerating parameter fits, we make no use of precomputed training sets or special choices of variables, and InterpMC is almost entirely transparent to the user.