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 Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

The Probe of Inflation and Cosmic Origins (PICO) is a NASA-funded study of a Probe-class mission concept. The top-level science objectives are to probe the physics of the Big Bang by measuring or constraining the energy scale of inflation, probe fundamental physics by measuring the number of light particles in the Universe and the sum of neutrino masses, to measure the reionization history of the Universe, and to understand the mechanisms driving the cosmic star formation history, and the physics of the galactic magnetic field. PICO would have multiple frequency bands between 21 and 799 GHz, and would survey the entire sky, producing maps of the polarization of the cosmic microwave background radiation, of galactic dust, of synchrotron radiation, and of various populations of point sources. Several instrument configurations, optical systems, cooling architectures, and detector and readout technologies have been and continue to be considered in the development of the mission concept. We will present a snapshot of the baseline mission concept currently under development.

The Probe of Inflation and Cosmic Origins (PICO) is a probe-class mission concept currently under study by NASA. PICO will probe the physics of the Big Bang and the energy scale of inflation, constrain the sum of neutrino masses, measure the growth of structures in the universe, and constrain its reionization history by making full sky maps of the cosmic microwave background with sensitivity 80 times higher than the Planck space mission. With bands at 21-799 GHz and arcmin resolution at the highest frequencies, PICO will make polarization maps of Galactic synchrotron and dust emission to observe the role of magnetic fields in Milky Way's evolution and star formation. We discuss PICO's optical system, focal plane, and give current best case noise estimates. The optical design is a two-reflector optimized open-Dragone design with a cold aperture stop. It gives a diffraction limited field of view (DLFOV) with throughput of 910 square cm sr at 21 GHz. The large 82 square degree DLFOV hosts 12,996 transition edge sensor bolometers distributed in 21 frequency bands and maintained at 0.1 K. We use focal plane technologies that are currently implemented on operating CMB instruments including three-color multi-chroic pixels and multiplexed readouts. To our knowledge, this is the first use of an open-Dragone design for mm-wave astrophysical observations, and the only monolithic CMB instrument to have such a broad frequency coverage. With current best case estimate polarization depth of 0.65 microK(CMB}-arcmin over the entire sky, PICO is the most sensitive CMB instrument designed to date.

EBEX was a long-duration balloon-borne experiment to measure the polarization of the cosmic microwave background. The experiment had three frequency bands centered at 150, 250, and 410 GHz and was the first to use a kilo-pixel array of transition edge sensor (TES) bolometers aboard a balloon platform; shortly after reaching float we operated 504, 342, and 109 TESs at each of the bands, respectively. We describe the design and characterization of the array and the readout system. We give the distributions of measured thermal conductances, normal resistances, and transition temperatures. With the exception of the thermal conductance at 150 GHz. We measured median low-loop-gain time constants $\tau_{0}=$ 88, 46, and 57 ms and compare them to predictions. Two measurements of bolometer absorption efficiency show high ($\sim$0.9) efficiency at 150 GHz and medium ($\sim$0.35, and $\sim$0.25) at the two higher bands, respectively. We measure a median total optical load of 3.6, 5.3 and 5.0 pW absorbed at the three bands, respectively. EBEX pioneered the use of the digital version of the frequency domain multiplexing (FDM) system which multiplexed the bias and readout of 16 bolometers onto two wires. We present accounting of the measured noise equivalent power. The median per-detector noise equivalent temperatures referred to a black body with a temperature of 2.725 K are 400, 920, and 14500 $\mu$K$\sqrt{s}$ for the three bands, respectively. We compare these values to our pre-flight predictions and to a previous balloon payload, discuss the sources of excess noise, and the path for a future payload to make full use of the balloon environment.

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.

Red novae are optical transients erupting at luminosities typically higher than those of classical novae. Their outbursts are believed to be caused by stellar mergers. We present millimeter/submillimeter-wave observations with ALMA and SMA of the three best known Galactic red novae, V4332 Sgr, V1309 Sco, and V838 Mon. The observations were taken 22, 8, and 14 yr after their respective eruptions and reveal the presence of molecular gas at excitation temperatures of 35-200 K. The gas displays molecular emission in rotational transitions with very broad lines (full width $\sim$400 km\s). We found emission of CO, SiO, SO, SO$_2$ (in all three red novae), H$_2$S (covered only in V838 Mon) and AlO (present in V4332 Sgr and V1309 Sco). No anomalies were found in the isotopic composition of the molecular material and the chemical (molecular) compositions of the three red novae appear similar to those of oxygen-rich envelopes of classical evolved stars (RSGs, AGBs, post-AGBs). The minimum masses of the molecular material that most likely was dispersed in the red-nova eruptions are 0.1, 0.01, and 10$^{-4}$ M$_{\odot}$ for V838 Mon, V4332 Sgr, and V1309 Sco, respectively. The molecular outflows in V4332 Sgr and V1309 Sco are spatially resolved and appear bipolar. The kinematic distances to V1309 Sco and V4332 Sgr are 2.1 and 4.2 kpc, respectively. The kinetic energy stored in the ejecta of the two older red-nova remnants of V838 Mon and V4332 Sgr is of order $10^{46}$ erg, similar to values found for some post-AGB (pre-PN) objects whose bipolar ejecta were also formed in a short-duration eruption. Our observations strengthen the link between these post-AGB objects and red novae and support the hypothesis that some of the post-AGB objects were formed in a common-envelope ejection event or its most catastrophic outcome, a merger.

We discuss a systematic effect associated with measuring polarization with a continuously rotating half-wave plate. The effect was identified with the data from the E and B Experiment (EBEX), which was a balloon-borne instrument designed to measure the polarization of the CMB as well as that from Galactic dust. The data show polarization fraction larger than 10\% while less than 3\% were expected from instrumental polarization. We give evidence that the excess polarization is due to detector non-linearity in the presence of a continuously rotating HWP. The non-linearity couples intensity signals into polarization. We develop a map-based method to remove the excess polarization. Applying this method for the 150 (250) GHz bands data we find that 81\% (92\%) of the excess polarization was removed. Characterization and mitigation of this effect is important for future experiments aiming to measure the CMB B-modes with a continuously rotating HWP.

We present our work to develop and characterize low thermal conductance bolometers that are part of sinuous antenna multi-chroic pixels (SAMP). We use longer, thinner and meandered bolometer legs to achieve 9 pW/K thermal conductance bolometers. We also discuss the development of inductor-capacitor chips operated at 4 K to extend the multiplexing factor of the frequency domain multiplexing to 105, an increase of 60% compared to the factor currently demonstrated for this readout system. This technology development is motivated by EBEX-IDS, a balloon-borne polarimeter designed to characterize the polarization of foregrounds and to detect the primordial gravity waves through their B-mode signature on the polarization of the cosmic microwave background. EBEX-IDS will operate 20,562 transition edge sensor bolometers spread over 7 frequency bands between 150 and 360 GHz. Balloon and satellite platforms enable observations at frequencies inaccessible from the ground and with higher instantaneous sensitivity. This development improves the readiness of the SAMP and frequency domain readout technologies for future satellite applications.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]

CMB-S4 is a proposed experiment to map the polarization of the Cosmic Microwave Background (CMB) to nearly the cosmic variance limit for angular scales that are accessible from the ground. The science goals and capabilities of CMB-S4 in illuminating cosmic inflation, measuring the sum of neutrino masses, searching for relativistic relics in the early universe, characterizing dark energy and dark matter, and mapping the matter distribution in the universe have been described in the CMB-S4 Science Book. This Technology Book is a companion volume to the Science Book. The ambitious science goals of CMB-S4, a "Stage-4" experiment, require a step forward in experimental capability from the current Stage=II experiments. To guide this process, we summarize the current state of CMB instrumentation technology, and identify R&D efforts necessary to advance it for use in CMB-S4. The book focuses on technical challenges in four broad areas: Telescope Design; Receiver Optics; Focal-Plane Optical Coupling; and Focal-Plane Sensor and Readout.

We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $\tau=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4\sigma$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged]

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We used two novel approaches to produce sub-wavelength structure (SWS) anti-reflection coatings (ARC) on silicon for the millimeter and sub-millimeter (MSM) wave band: picosecond laser ablation and dicing with beveled saws. We produced pyramidal structures with both techniques. The diced sample, machined on only one side, had pitch and height of 350 $\mu$m and 972 $\mu$m. The two laser ablated samples had pitch of 180 $\mu$m and heights of 720 $\mu$m and 580 $\mu$m; only one of these samples was ablated on both sides. We present measurements of shape and optical performance as well as comparisons to the optical performance predicted using finite element analysis and rigorous coupled wave analysis. By extending the measured performance of the one-sided diced sample to the two-sided case, we demonstrate 25 % band averaged reflectance of less than 5 % over a bandwidth of 97 % centered on 170 GHz. Using the two-sided laser ablation sample, we demonstrate reflectance less than 5 % over 83 % bandwidth centered on 346 GHz.

The E and B Experiment (EBEX) was a long-duration balloon-borne cosmic microwave background polarimeter that flew over Antarctica in 2013. We describe the experiment's optical system, receiver, and polarimetric approach, and report on their in-flight performance. EBEX had three frequency bands centered on 150, 250, and 410 GHz. To make efficient use of limited mass and space we designed a 115 cm$^{2}$sr high throughput optical system that had two ambient temperature mirrors and four anti-reflection coated polyethylene lenses per focal plane. All frequency bands shared the same optical train. Polarimetry was achieved with a continuously rotating achromatic half-wave plate (AHWP) that was levitated with a superconducting magnetic bearing (SMB). Rotation stability was 0.45 % over a period of 10 hours, and angular position accuracy was 0.01 degrees. This is the first use of a SMB in astrophysics. The measured modulation efficiency was above 90 % for all bands. To our knowledge the 109 % fractional bandwidth of the AHWP was the broadest implemented to date. The receiver that contained one lens and the AHWP at a temperature of 4 K, the polarizing grid and other lenses at 1 K, and the two focal planes at 0.25 K performed according to specifications giving focal plane temperature stability with fluctuation power spectrum that had $1/f$ knee at 2 mHz. EBEX was the first balloon-borne instrument to implement technologies characteristic of modern CMB polarimeters including high throughput optical systems, and large arrays of transition edge sensor bolometric detectors with mutiplexed readouts.

We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 {\mu}K$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 {\mu}K$\cdot$arcmin.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

The E and B Experiment (EBEX) was a long-duration balloon-borne instrument designed to measure the polarization of the cosmic microwave background (CMB) radiation. EBEX was the first balloon-borne instrument to implement a kilo-pixel array of transition edge sensor (TES) bolometric detectors and the first CMB experiment to use the digital version of the frequency domain multiplexing system for readout of the TES array. The scan strategy relied on 40 s peak-to-peak constant velocity azimuthal scans. We discuss the unique demands on the design and operation of the payload that resulted from these new technologies and the scan strategy. We describe the solutions implemented including the development of a power system designed to provide a total of at least 2.3 kW, a cooling system to dissipate 590 W consumed by the detectors' readout system, software to manage and handle the data of the kilo-pixel array, and specialized attitude reconstruction software. We present flight performance data showing faultless management of the TES array, adequate powering and cooling of the readout electronics, and constraint of attitude reconstruction errors such that the spurious B-modes they induced were less than 10% of CMB B-mode power spectrum with $r=0.05$.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

Based on sub-arcsecond Atacama Large Millimeter/submillimeter Array (ALMA) and Submillimeter Array (SMA) 1.3 mm continuum images of the massive protocluster NGC 6334I obtained in 2015 and 2008, we find that the dust emission from MM1 has increased by a factor of 4.0$\pm$0.3 during the intervening years, and undergone a significant change in morphology. The continuum emission from the other cluster members (MM2, MM4 and the UCHII region MM3 = NGC 6334F) has remained constant. Long term single-dish maser monitoring at HartRAO finds that multiple maser species toward NGC 6334I flared beginning in early 2015, a few months before our ALMA observation, and some persist in that state. New ALMA images obtained in 2016 July-August at 1.1 and 0.87 mm confirm the changes with respect to SMA 0.87 mm images from 2008, and indicate that the (sub)millimeter flaring has continued for at least a year. The excess continuum emission, centered on the hypercompact HII region MM1B, is extended and elongated ($1.6" \times 1.0" \approx 2100 \times 1300$ au) with multiple peaks, suggestive of general heating of the surrounding subcomponents of MM1, some of which may trace clumps in a fragmented disk rather than separate protostars. In either case, these remarkable increases in maser and dust emission provide direct observational evidence of a sudden accretion event in the growth of a massive protostar yielding a sustained luminosity surge by a factor of $70\pm20$, analogous to the largest events in simulations by Meyer et al. (2017). This target provides an excellent opportunity to assess the impact of such a rare event on a protocluster over many years.

A 32 GHz bandwidth VLBI capable correlator and phased array has been designed and deployed at the Smithsonian Astrophysical Observatory's Submillimeter Array (SMA). The SMA Wideband Astronomical ROACH2 Machine (SWARM) integrates two instruments: a correlator with 140 kHz spectral resolution across its full 32 GHz band, used for connected interferometric observations, and a phased array summer used when the SMA participates as a station in the Event Horizon Telescope (EHT) Very Long Baseline Interferometry (VLBI) array. For each SWARM quadrant, Reconfigurable Open Architecture Computing Hardware (ROACH2) units shared under open source from the Collaboration for Astronomy Signal Processing and Electronics Research (CASPER) are equipped with a pair of ultra-fast Analog-to- Digital Converters (ADCs), a Field Programmable Gate Array (FPGA) processor, and eight 10 Gigabit Ethernet ports. A VLBI data recorder interface designated the SWARM Digital Back End, or SDBE, is implemented with a ninth ROACH2 per quadrant, feeding four Mark6 VLBI recorders with an aggregate recording rate of 64 Gbps. This paper describes the design and implementation of SWARM, as well as its deployment at SMA with reference to verification and science data.

We report on the first use of laser ablation to make sub-millimeter, broad-band, anti-reflection coatings (ARC) based on sub-wavelength structures (SWS) on alumina and sapphire. We used a 515 nm laser to produce pyramid-shaped structures with pitch of about 320 $\mu$m and total height of near 800 $\mu$m. Transmission measurements between 70 and 140 GHz are in agreement with simulations using electromagnetic propagation software. The simulations indicate that SWS ARC with the fabricated shape should have a fractional bandwidth response of $\Delta \nu / \nu_{center}\ = 0.55$ centered on 235 GHz for which reflections are below 3%. Extension of the bandwidth to both lower and higher frequencies, between few tens of GHz and few THz, should be straightforward with appropriate adjustment of laser ablation parameters.

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.

The E and B Experiment (EBEX) is a balloon-borne polarimeter designed to measure the polarization of the cosmic microwave background radiation and to characterize the polarization of galactic dust. EBEX was launched December 29, 2012 and circumnavigated Antarctica observing $\sim$6,000 square degrees of sky during 11 days at three frequency bands centered around 150, 250 and 410 GHz. EBEX was the first experiment to operate a kilo-pixel array of transition-edge sensor bolometers and a continuously rotating achromatic half-wave plate aboard a balloon platform. It also pioneered the use of detector readout based on digital frequency domain multiplexing. We describe the temperature calibration of the experiment. The gain response of the experiment is calibrated using a two-step iterative process. We use signals measured on passes across the Galactic plane to convert from readout-system counts to power. The effective smoothing scale of the EBEX optics and the star camera-to-detector offset angles are determined through \c{hi}2 minimization using the compact HII region RCW 38. This two-step process is initially performed with parameters measured before the EBEX 2013 flight and then repeated until the calibration factor and parameters converge.

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 present the full catalog of Young Stellar Objects (YSOs) identified in the 18 molecular clouds surveyed by the Spitzer Space Telescope "cores to disks" (c2d) and "Gould Belt" (GB) Legacy surveys. Using standard techniques developed by the c2d project, we identify 3239 candidate YSOs in the 18 clouds, 2966 of which survive visual inspection and form our final catalog of YSOs in the Gould Belt. We compile extinction corrected SEDs for all 2966 YSOs and calculate and tabulate the infrared spectral index, bolometric luminosity, and bolometric temperature for each object. We find that 326 (11%), 210 (7%), 1248 (42%), and 1182 (40%) are classified as Class 0+I, Flat-spectrum, Class II, and Class III, respectively, and show that the Class III sample suffers from an overall contamination rate by background AGB stars between 25% and 90%. Adopting standard assumptions, we derive durations of 0.40-0.78 Myr for Class 0+I YSOs and 0.26-0.50 Myr for Flat-spectrum YSOs, where the ranges encompass uncertainties in the adopted assumptions. Including information from (sub)millimeter wavelengths, one-third of the Class 0+I sample is classified as Class 0, leading to durations of 0.13-0.26 Myr (Class 0) and 0.27-0.52 Myr (Class I). We revisit infrared color-color diagrams used in the literature to classify YSOs and propose minor revisions to classification boundaries in these diagrams. Finally, we show that the bolometric temperature is a poor discriminator between Class II and Class III YSOs.

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 Spitzer Space Telescope mapped the Perseus molecular cloud complex with IRAC and MIPS as part of the c2d Spitzer Legacy project. This paper combines the observations from both instruments giving an overview of low-mass star formation across Perseus from 3.6 to 70 micron. We provide an updated list of young stellar objects with new classifications and source fluxes from previous works, identifying 369 YSOs in Perseus with the Spitzer dataset. By synthesizing the IRAC and MIPS maps of Perseus and building on the work of previous papers in this series (Jorgensen et al. 2006, Rebull et al. 2007), we present a current census of star formation across the cloud and within smaller regions. 67% of the YSOs are associated with the young clusters NGC 1333 and IC 348. The majority of the star formation activity in Perseus occurs in the regions around the clusters, to the eastern and western ends of the cloud complex. The middle of the cloud is nearly empty of YSOs despite containing regions of high visual extinction. The western half of Perseus contains three-quarters of the total number of embedded YSOs (Class 0+I and Flat SED sources) in the cloud and nearly as many embedded YSOs as Class II and III sources. Class II and III greatly outnumber Class 0+I objects in eastern Perseus and IC 348. These results are consistent with previous age estimates for the clusters. Across the cloud, 56% of YSOs and 91% of the Class 0+I and Flat sources are in areas where Av > 5 mag, indicating a possible extinction threshold for star formation.

We report the discovery in space of a disilicon species, SiCSi, from observations between 80 and 350 GHz with the IRAM 30m radio telescope. Owing to the close coordination between laboratory experiments and astrophysics, 112 lines have now been detected in the carbon-rich star CWLeo. The derived frequencies yield improved rotational and centrifugal distortion constants up to sixth order. From the line profiles and interferometric maps with the Submillimeter Array, the bulk of the SiCSi emis- sion arises from a region of 6 arcseconds in radius. The derived abundance is comparable to that of SiC2. As expected from chemical equilibrium calculations, SiCSi and SiC2 are the most abundant species harboring a SiC bond in the dust formation zone and certainly both play a key role in the formation of SiC dust grains.

Rotational spectra in four new excited vibrational levels of the linear carbon chain radical C$_4$H radical were observed in the millimeter band between 69 and 364 GHz in a low pressure glow discharge, and two of these were observed in a supersonic molecular beam between 19 and 38 GHz. All have rotational constants within 0.4% of the $^2\Sigma^+$ ground vibrational state of C$_4$H and were assigned to new bending vibrational levels, two each with $^2\Sigma$ and $^2\Pi$ vibrational symmetry. The new levels are tentatively assigned to the $1\nu_6$ and $1\nu_5$ bending vibrational modes (both with $^2\Pi$ symmetry), and the $1\nu_6 + 1\nu_7$ and $1\nu_5 + 1\nu_6$ combination levels ($^2\Sigma$ symmetry) on the basis of the derived spectroscopic constants, relative intensities in our discharge source, and published laser spectroscopic and quantum calculations. Prior spectroscopic constants in the $1\nu_7$ and $2\nu_7$ levels were refined. Also presented are interferometric maps of the ground state and the $1\nu_7$ level obtained with the SMA near 257 GHz which show that C$_4$H is present near the central star in IRC+10216. We found no evidence with the SMA for the new vibrationally excited levels of C$_4$H at a peak flux density averaged over a $3^{\prime\prime}$ synthesized beam of $\ge 0.15$ Jy/beam in the 294-296 and 304-306 GHz range, but it is anticipated that rotational lines in the new levels might be observed in IRC+10216 when ALMA attains its full design capability.

Using the SMA and VLA, we have imaged the massive protocluster NGC6334I(N) at high angular resolution (0.5"~650AU) from 6cm to 0.87mm, detecting 18 new compact continuum sources. Three of the new sources are coincident with previously-identified water masers. Together with the previously-known sources, these data bring the number of likely protocluster members to 25 for a protostellar density of ~700 pc^-3. Our preliminary measurement of the Q-parameter of the minimum spanning tree is 0.82 -- close to the value for a uniform volume distribution. All of the (nine) sources with detections at multiple frequencies have SEDs consistent with dust emission, and two (SMA1b and SMA4) also have long wavelength emission consistent with a central hypercompact HII region. Thermal spectral line emission, including CH3CN, is detected in six sources: LTE model fitting of CH3CN(J=12-11) yields temperatures of 72-373K, confirming the presence of multiple hot cores. The fitted LSR velocities range from -3.3 to -7.0 km/s, with an unbiased mean square deviation of 2.05 km/s, implying a dynamical mass of 410+-260 Msun for the protocluster. From analysis of a wide range of hot core molecules, the kinematics of SMA1b are consistent with a rotating, infalling Keplerian disk of diameter 800AU and enclosed mass of 10-30 Msun that is perpendicular (within 1 degree) to the large-scale bipolar outflow axis. A companion to SMA1b at a projected separation of 0.45" (590AU; SMA1d), which shows no evidence of spectral line emission, is also confirmed. Finally, we detect one 218.440GHz and several 229.7588GHz Class-I methanol masers.

A spectral line survey of the oxygen-rich red supergiant VY Canis Majoris was made between 279 and 355 GHz with the Submillimeter Array. Two hundred twenty three spectral features from 19 molecules (not counting isotopic species of some of them) were observed, including the rotational spectra of TiO, TiO2, and AlCl for the first time in this source. The parameters and an atlas of all spectral features is presented. Observations of each line with a synthesized beam of ~0.9 arcsec, reveal the complex kinematics and morphology of the nebula surrounding VY CMa. Many of the molecules are observed in high lying rotational levels or in excited vibrational levels. From these, it was established that the main source of the submillimeter-wave continuum (dust) and the high excitation molecular gas (the star) are separated by about 0.15 arcsec. Apparent coincidences between the molecular gas observed with the SMA, and some of the arcs and knots observed at infrared wavelengths and in the optical scattered light by the Hubble Space Telescope are identified. The observations presented here provide important constraints on the molecular chemistry in oxygen-dominated circumstellar environments and a deeper picture of the complex circumstellar environment of VY CMa.

Phosphorus-bearing compounds have only been studied in the circumstellar environments (CSEs) of the asymptotic giant branch (AGB) star IRC +10216 and the protoplanetary nebula CRL 2688, both C-rich objects, and the O-rich red supergiant VY CMa. The current chemical models cannot reproduce the high abundances of PO and PN derived from observations of VY CMa. No observations have been reported of phosphorus in the CSEs of O-rich AGB stars. We aim to set observational constraints on the phosphorous chemistry in the CSEs of O-rich AGB stars, by focussing on the Mira-type variable star IK Tau. Using the IRAM 30m telescope and the Submillimeter Array (SMA), we observed four rotational transitions of PN (J=2-1,3-2,6-5,7-6) and four of PO (J=5/2-3/2,7/2-5/2,13/2-11/2,15/2-13/2). The IRAM 30m observations were dedicated line observations, while the SMA data come from an unbiased spectral survey in the frequency range 279-355 GHz. We present the first detections of PN and PO in an O-rich AGB star and estimate abundances X(PN/H2) of about 3x10^-7 and X(PO/H2) in the range 0.5-6.0x10^-7. This is several orders of magnitude higher than what is found for the C-rich AGB star IRC +10216. The diameter (<=0.7") of the PN and PO emission distributions measured in the interferometric data corresponds to a maximum radial extent of about 40 stellar radii. The abundances and the spatial occurrence of the molecules are in very good agreement with the results reported for VY CMa. We did not detect PS or PH3 in the survey. We suggest that PN and PO are the main carriers of phosphorus in the gas phase, with abundances possibly up to several 10^-7. The current chemical models cannot account for this, underlining the strong need for updated chemical models that include phosphorous compounds.

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.

We describe the design and performance of the digital servo and motion control system for the 6-meter diameter parabolic antennas of the Submillimeter Array (SMA) on Mauna Kea, Hawaii. The system is divided into three nested layers operating at a different, appropriate bandwidth. (1) A rack-mounted, real-time Unix system runs the position loop which reads the high resolution azimuth and elevation encoders and sends velocity and acceleration commands at 100 Hz to a custom-designed servo control board (SCB). (2) The microcontroller-based SCB reads the motor axis tachometers and implements the velocity loop by sending torque commands to the motor amplifiers at 558 Hz. (3) The motor amplifiers implement the torque loop by monitoring and sending current to the three-phase brushless drive motors at 20 kHz. The velocity loop uses a traditional proportional-integral-derivative (PID) control algorithm, while the position loop uses only a proportional term and implements a command shaper based on the Gauss error function. Calibration factors and software filters are applied to the tachometer feedback prior to the application of the servo gains in the torque computations. All of these parameters are remotely adjustable in software. The three layers of the control system monitor each other and are capable of shutting down the system safely if a failure or anomaly occurs. The Unix system continuously relays antenna status to the central observatory computer via reflective memory. In each antenna, a Palm Vx hand controller displays system status and allows full local control of the drives in an intuitive touchscreen user interface. It can also be connected outside the cabin for convenience during antenna reconfigurations. Excellent tracking performance (0.3 arcsec rms) is achieved with this system. It has been in reliable operation on 8 antennas for over 10 years and has required minimal maintenance.

The formation of inorganic dust in circumstellar environments of evolved stars is poorly understood. Spectra of molecules thought to be most important for the nucleation, i.e. AlO, TiO, and TiO2, have been recently detected in the red supergiant VY CMa. These molecules are effectively formed in VY CMa and the observations suggest that non-equilibrium chemistry must be involved in their formation and nucleation into dust. In addition to exploring the recent observations of VY CMa, we briefly discuss the possibility of detecting these molecules in the dust-poor circumstellar environment of Betelgeuse.

We report the first detection of pure rotational transitions of TiO and TiO_2 at (sub-)millimeter wavelengths towards the red supergiant VY CMa. A rotational temperature, T_rot, of about 250 K was derived for TiO_2. Although T_rot was not well constrained for TiO, it is likely somewhat higher than that of TiO_2. The detection of the Ti oxides confirms that they are formed in the circumstellar envelopes of cool oxygen-rich stars and may be the "seeds" of inorganic-dust formation, but alternative explanations for our observation of TiO and TiO_2 in the cooler regions of the envelope cannot be ruled out at this time. The observations suggest that a significant fraction of the oxides is not converted to dust, but instead remains in the gas phase throughout the outflow.

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.

We have constructed a five station 12 GHz atmospheric phase interferometer (API) for the Submillimeter Array (SMA) located near the summit of Mauna Kea, Hawaii. Operating at the base of unoccupied SMA antenna pads, each station employs a commercial low noise mixing block coupled to a 0.7 m off-axis satellite dish which receives a broadband, white noise-like signal from a geostationary satellite. The signals are processed by an analog correlator to produce the phase delays between all pairs of stations with projected baselines ranging from 33 to 261 m. Each baseline's amplitude and phase is measured continuously at a rate of 8 kHz, processed, averaged and output at 10 Hz. Further signal processing and data reduction is accomplished with a Linux computer, including the removal of the diurnal motion of the target satellite. The placement of the stations below ground level with an environmental shield combined with the use of low temperature coefficient, buried fiber optic cables provides excellent system stability. The sensitivity in terms of rms path length is 1.3 microns which corresponds to phase deviations of about 1 degree of phase at the highest operating frequency of the SMA. The two primary data products are: (1) standard deviations of observed phase over various time scales, and (2) phase structure functions. These real-time statistical data measured by the API in the direction of the satellite provide an estimate of the phase front distortion experienced by the concurrent SMA astronomical observations. The API data also play an important role, along with the local opacity measurements and weather predictions, in helping to plan the scheduling of science observations on the telescope.

We report Submillimeter Array observations of SO2 emission in the circumstellar envelope of the red supergiant VY CMa, with an angular resolution of ~1". SO2 emission appears in three distinct outflow regions surrounding the central continuum peak emission that is spatially unresolved. No bipolar structure is noted in the sources. A fourth source of SO2 is identified as a spherical wind centered at the systemic velocity. We estimate the SO2 column density and rotational temperature assuming local thermal equilibrium (LTE) as well as perform non-LTE radiative transfer analysis using RADEX. Column densities of SO2 are found to be ~10^16 cm^-2 in the outflows and in the spherical wind. Comparison with existing maps of the two parent species OH and SO shows the SO2 distribution to be consistent with that of OH. The abundance ratio f_SO2/f_SO is greater than unity for all radii greater than at least 3x10^16 cm. SO2 is distributed in fragmented clumps compared to SO, PN, and SiS molecules. These observations lend support to specific models of circumstellar chemistry that predict f_SO2/f_SO>1 and may suggest the role of localized effects such as shocks in the production of SO2 in the circumstellar envelope.

We report a spectral-line survey of the extreme carbon star IRC+10216 carried out between 293.9 and 354.8 GHz with the Submillimeter Array. A total of 442 lines were detected, more than 200 for the first time; 149 are unassigned. Maps at an angular resolution of ~3" were obtained for each line. A substantial new population of narrow lines with an expansion velocity of ~4 km/s (i.e. ~ 30% of the terminal velocity) was detected. Most of these are attributed to rotational transitions within vibrationally excited states, emitted from energy levels above the v=0, J=0 ground state with excitation energy of 1000-3000 K. Emission from these lines appears to be centered on the star with an angular extent of <1". We use multiple transitions detected in several molecules to derive physical conditions in this inner envelope of IRC+10216.

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

Observations of Lynds Dark Nebula 1221 from the Spitzer Space Telescope are presented. These data show three candidate protostars towards L1221, only two of which were previously known. The infrared observations also show signatures of outflowing material, an interpretation which is also supported by radio observations with the Very Large Array. In addition, molecular line maps from the Five College Radio Astronomy Observatory are shown. One-dimensional dust continuum modelling of two of these protostars, IRS1 and IRS3, is described. These models show two distinctly different protostars forming in very similar environments. IRS1 shows a higher luminosity and larger inner radius of the envelope than IRS3. The disparity could be caused by a difference in age or mass, orientation of outflow cavities, or the impact of a binary in the IRS1 core.

We used the Extended Submillimeter Array (eSMA) in its most extended configuration to investigate the innermost (within a radius of 290 R* from the star) circumstellar envelope (CSE) of IRC+10216. We imaged the CSE using HCN and other molecular lines with a beam size of 0."22 x 0."46, deeply into the very inner edge (15 R*) of the envelope where the expansion velocity is only 3 km/s. The excitation mechanism of hot HCN and KCl maser lines is discussed. HCN maser components are spatially resolved for the first time on an astronomical object. We identified two discrete regions in the envelope: a region with a radius of . 15 R*, where molecular species have just formed and the gas has begun to be accelerated (region I) and a shell region (region II) with a radius of 23 R* and a thickness of 15 R*, whose expansion velocity has reached up to 13 km/s, nearly the terminal velocity of 15 km/s. The Si$^{34}$S line detected in region I shows a large expansion velocity of 16 km/s due to strong wing components, indicating that the emission may arise from a shock region in the innermost envelope. In region II, the P.A. of the most copious mass loss direction was found to be 120 +/- 10 degrees, which may correspond to the equatorial direction of the star. Region II contains a torus-like feature. These two regions may have emerged due to significant differences in the size distributions of the dust particles in the two regions.

We aim to study the structure, dynamics and physical conditions of Gomez's Hamburger (IRAS 18059-3211; GoHam). We confirm that GoHam essentially consists of a flaring disk in keplerian rotation around a young, probably pre-MS star. We present high resolution SMA maps of 12CO J=2-1, 13CO J=2-1, 12CO J=3-2, and C17O J=3-2, as well as data on 12CO J=6-5 and the continuum flux at these wavelengths. Spatial resolutions up to 1" are obtained. Except for the C17O data, the dynamical ranges are larger than 10. The maps are compared to a numerical model, which simulates the emission of a rotating disk with the expected general properties of such objects; a very satisfactory fitting of our maps is obtained. The meaning and reliability of our results are thoroughly discussed. Our observations allow measurement of the main properties of GoHam at scales between ~ 1" (~ 5 10^15 cm, for the assumed distance, 300 pc) and the total extent of the nebula, 14". We are able to measure the global structure of the gas-rich disk, which is found to be flaring, and its dynamics, which is clearly dominated by keplerian rotation, with a very small degree of turbulence. The combination of different lines, particularly showing different opacities, allows us to reasonably estimate the distributions of the gas temperature and density. We clearly find a significant and sharp increase in temperature at large distances from the equator, accompanied by a decrease in density of the same order. Finally, we identify and study a condensation in the southern part of the disk that has no counterparts in the rest of the nebula. This condensation is quite extended (about 5 10^15 cm), contains a significant amount of mass (roughly, ~ 6 10^-3 Mo), and seems to be associated with a detectable distortion of the global rotation kinematics.

Using the Submillimeter Array we have detected the J=3-2 and 2-1 rotational transitions from within the first vibrationally excited state of CO toward the extreme carbon star IRC+10216 (CW Leo). The emission remains spatially unresolved with an angular resolution of ~2" and, given that the lines originate from energy levels that are ~3100 K above the ground state, almost certainly originates from a much smaller (~10^{14} cm) sized region close to the stellar photosphere. Thermal excitation of the lines requires a gas density of ~10^{9} cm^{-3}, about an order of magnitude higher than the expected gas density based previous infrared observations and models of the inner dust shell of IRC+10216.

We report the first science observations and results obtained with the "extended" SMA (eSMA), which is composed of the SMA (Submillimeter Array), JCMT (James Clerk Maxwell Telescope) and CSO (Caltech Submillimeter Observatory). Redshifted absorptions at z=0.886 of CI (^3P_1 - ^3P_0) were observed with the eSMA with an angular resolution of 0.55"x0.22" at 1.1 mm toward the southwestern image of the remarkable lensed quasar PKS 1830-211, but not toward the northeastern component at a separation of ~1". Additionally, SMA observations of CO, 13CO and C18O (all J=4-3) were obtained toward this object: CO was also detected toward the SW component, but none of the isotopologues were. This is the first time [CI] is detected in this object, allowing the first direct determination of relative abundances of neutral atomic carbon to CO in the molecular clouds of a spiral galaxy at z>0.1. The [CI] and CO profiles can be decomposed into two and three velocity components respectively. We derive C/CO column density ratios ranging from <0.5 (representative of dense cores) to ~2.5 (close to translucent clouds values). This could indicate that we are seeing environments with different physical conditions or that we are witnessing chemical evolution of regions where C has not completely been converted into CO.

A spectral-line survey of IRC+10216 in the 345 GHz band has been undertaken with the Submillimeter Array. Although not yet completed, it has already yielded a fairly large sample of narrow molecular emission lines with line-widths indicating expansion velocities of ~4 km/s, less than 3 times the well-known value of the terminal expansion velocity (14.5 km/s) of the outer envelope. Five of these narrow lines have now been identified as rotational transitions in vibrationally excited states of previously detected molecules: the v=1, J=17--16 and J=19--18 lines of Si34S and 29SiS and the v=2, J=7--6 line of CS. Maps of these lines show that the emission is confined to a region within ~60 AU of the star, indicating that the narrow-line emission is probing the region of dust-formation where the stellar wind is still being accelerated.

The cores of most galaxies are thought to harbour supermassive black holes, which power galactic nuclei by converting the gravitational energy of accreting matter into radiation (ref 1). Sagittarius A*, the compact source of radio, infrared and X-ray emission at the centre of the Milky Way, is the closest example of this phenomenon, with an estimated black hole mass that is 4 million times that of the Sun (refs. 2,3). A long-standing astronomical goal is to resolve structures in the innermost accretion flow surrounding Sgr A* where strong gravitational fields will distort the appearance of radiation emitted near the black hole. Radio observations at wavelengths of 3.5 mm and 7 mm have detected intrinsic structure in Sgr A*, but the spatial resolution of observations at these wavelengths is limited by interstellar scattering (refs. 4-7). Here we report observations at a wavelength of 1.3 mm that set a size of 37 (+16, -10; 3-sigma) microarcseconds on the intrinsic diameter of Sgr A*. This is less than the expected apparent size of the event horizon of the presumed black hole, suggesting that the bulk of SgrA* emission may not be not centred on the black hole, but arises in the surrounding accretion flow.

The eSMA ("extended SMA") combines the SMA, JCMT and CSO into a single facility, providing enhanced sensitivity and spatial resolution owing to the increased collecting area at the longest baselines. Until ALMA early science observing (2011), the eSMA will be the facility capable of the highest angular resolution observations at 345 GHz. The gain in sensitivity and resolution will bring new insights in a variety of fields, such as protoplanetary/transition disks, high-mass star formation, solar system bodies, nearby and high-z galaxies. Therefore the eSMA is an important facility to prepare the grounds for ALMA and train scientists in the techniques. Over the last two years, and especially since November 2006, there has been substantial progress toward making the eSMA into a working interferometer. In particular, (i) new 345-GHz receivers, that match the capabilities of the SMA system, were installed at the JCMT and CSO; (ii) numerous tests have been performed for receiver, correlator and baseline calibrations in order to determine and take into account the effects arising from the differences between the three types of antennas; (iii) first fringes at 345 GHz were obtained on August 30 2007, and the array has entered the science-verification stage. We report on the characteristics of the eSMA and its measured performance at 230 GHz and that expected at 345 GHz. We also present the results of the commissioning and some initial science-verification observations, including the first absorption measurement of the C/CO ratio in a galaxy at z=0.89, located along the line of sight to the lensed quasar PKS1830-211, and on the imaging of the vibrationally excited HCN line towards IRC+10216.

We study the nature of Gomez's Hamburger (IRAS 18059-3211), a nebula that has been proposed to be a post-AGB object. Such a classification is not confirmed; instead, we argue that it will be a key object in the study of disks rotating around young stars. We present high resolution SMA maps of CO J=2--1 in Gomez's Hamburger. The data are analyzed by means of a code that simulates the emission of a nebula showing a variety of physical conditions and kinematics. Our observations clearly show that the CO emitting gas in Gomez's Hamburger forms a spectacular disk in keplerian rotation. Model calculations undoubtly confirm this result. The central (mainly stellar) mass is found to be high, ~ 4 Mo for a distance of 500 pc. The mass and (relatively low) luminosity of the source are, independent of the assumed distance, very different from those possible in evolved stars. Gomez's Hamburger is probably a transitional object between the pre-MS and MS phases, still showing interstellar material around the central star or stellar system.

We used the Submillimeter Array to map the angular distribution of the H30$\alpha$ recombination line (231.9 GHz) in the circumstellar region of the peculiar star MWC349A. The resolution was $1\farcs2$, but because of high signal-to-noise ratio we measured the positions of all maser components to accuracies better than $0\farcs01$, at a velocity resolution of $1 kms$. The two strongest maser components (called high velocity components) at velocities near -14 and $32 kms$ are separated by $0\farcs048 \pm 0\farcs001$ (60 AU) along a position angle of $102 \pm 1\arcdeg$. The distribution of maser emission at velocities between and beyond these two strongest components were also provided. The continuum emission lies at the center of the maser distribution to within 10 mas. The masers appear to trace a nearly edge-on rotating disk structure, reminiscent of the water masers in Keplerian rotation in the nuclear accretion disk of the galaxy NGC4258. However, the maser components in MWC349A do not follow a simple Keplerian kinematic prescription with $v \sim r^{-1/2}$, but have a larger power law index. We explore the possibility that the high velocity masers trace spiral density or shock waves. We also emphasize caution in the interpretation of relative centroid maser positions where the maser is not clearly resolved in position or velocity, and we present simulations that illustrate the range of applicability of the centroiding method.

We present maps of 14.4 deg^2 of the Ophiuchus dark clouds observed by the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). These high quality maps depict both numerous point sources as well as extended dust emission within the star-forming and non-star-forming portions of these clouds. Using PSF-fitting photometry, we detect 5779 sources at 24 um and 81 sources at 70 um at the 10 sigma level of significance. Three hundred twenty-three candidate young stellar objects (YSOs) were identified according to their positions on the MIPS/2MASS K versus K$-$[24] color-magnitude diagrams as compared to 24 um detections in the SWIRE extragalactic survey. We find that more than half of the YSO candidates, and almost all the ones with protostellar Class I spectral energy distributions, are confined to the known cluster and aggregates.

We report on the results of a Submillimeter Array interferometric observation of the proto-planetary nebula CRL 618 in the 12CO J=6-5 line. With the new capability of SMA enabling us to use two receivers at a time, we also observed simultaneously in the 12CO J=2-1 and 13CO J=2-1 lines. The 12CO J=6-5 and 13CO J=2-1 lines were first interferometrically observed toward CRL 618. The flux of the high velocity component of the 12CO J=6-5 line is almost fully recovered, while roughly 80% of the flux of the low velocity component is resolved out. The low recovery rate suggests that the emission region of the low velocity component of the 12CO J=6-5 line is largely extended. Continuum emission is detected both at 230 and 690 GHz. The flux of the 690 GHz continuum emission seems to be partially resolved out, suggesting dust emission partly contaminates the 690 GHz continuum flux. The cavity structure, which has been confirmed in a previous observation in the 12CO J=2-1 line, is not clearly detected in the 12CO J=6-5 line, and only the south wall of the cavity is detected. This result suggests that the physical condition of the molecular envelope of CRL 618 is not exactly axial symmetric.

H13CN J=8-7 sub-millimetre line emission produced in the circumstellar envelope around the extreme carbon star IRC+10216 has been imaged at sub-arcsecond angular resolution using the SMA. Supplemented by a detailed excitation analysis the average fractional abundance of H13CN in the inner wind (< 5E15 cm) is estimated to be about 4E-7, translating into a total HCN fractional abundance of 2E-5 using the isotopic ratio 12C/13C=50. Multi-transitional single-dish observations further requires the H13CN fractional abundance to remain more or less constant in the envelope out to a radius of about 4E16 cm, where the HCN molecules are effectively destroyed, most probably, by photodissociation. The large amount of HCN present in the inner wind provides effective line cooling that can dominate over that generated from CO line emission. It is also shown that great care needs to be taken in the radiative transfer modelling where non-local, and non-LTE, effects are important and where the radiation field from thermal dust grains plays a major role in exciting the HCN molecules. The amount of HCN present in the circumstellar envelope around IRC+10216 is consistent with predicted photospheric values based on equilibrium chemical models and indicates that any non-equilibrium chemistry occurring in the extended pulsating atmosphere has no drastic net effect on the fractional abundance of HCN molecules that enters the outer envelope. It further suggests that few HCN molecules are incorporated into dust grains.

We present maps of 7.78 square degrees of the Lupus molecular cloud complex at 24, 70, and $160\:\mu$m. They were made with the Spitzer Space Telescope's Multiband Imaging Photometer for Spitzer (MIPS) instrument as part of the Spitzer Legacy Program, ``From Molecular Cores to Planet-Forming Disks'' (c2d). The maps cover three separate regions in Lupus, denoted I, III, and IV. We discuss the c2d pipeline and how our data processing differs from it. We compare source counts in the three regions with two other data sets and predicted star counts from the Wainscoat model. This comparison shows the contribution from background galaxies in Lupus I. We also create two color magnitude diagrams using the 2MASS and MIPS data. From these results, we can identify background galaxies and distinguish them from probable young stellar objects. The sources in our catalogs are classified based on their spectral energy distribution (SED) from 2MASS and Spitzer wavelengths to create a sample of young stellar object candidates. From 2MASS data, we create extinction maps for each region and note a strong corresponence between the extinction and the $160\:\mu$m emission. The masses we derived in each Lupus cloud from our extinction maps are compared to masses estimated from $^{13}$CO and C$^{18}$O and found to be similar to our extinction masses in some regions, but significantly different in others. Finally, based on our color-magnitude diagrams, we selected 12 of our reddest candidate young stellar objects for individual discussion. Five of the 12 appear to be newly-discovered YSOs.

We have undertaken an unprecedentedly large 1.1 millimeter continuum survey of three nearby star forming clouds using Bolocam at the Caltech Submillimeter Observatory. We mapped the largest areas in each cloud at millimeter or submillimeter wavelengths to date: 7.5 sq. deg in Perseus (Paper I), 10.8 sq. deg in Ophiuchus (Paper II), and 1.5 sq. deg in Serpens with a resolution of 31", detecting 122, 44, and 35 cores, respectively. Here we report on results of the Serpens survey and compare the three clouds. Average measured angular core sizes and their dependence on resolution suggest that many of the observed sources are consistent with power-law density profiles. Tests of the effects of cloud distance reveal that linear resolution strongly affects measured source sizes and densities, but not the shape of the mass distribution. Core mass distribution slopes in Perseus and Ophiuchus (alpha=2.1+/-0.1 and alpha=2.1+/-0.3) are consistent with recent measurements of the stellar IMF, whereas the Serpens distribution is flatter (alpha=1.6+/-0.2). We also compare the relative mass distribution shapes to predictions from turbulent fragmentation simulations. Dense cores constitute less than 10% of the total cloud mass in all three clouds, consistent with other measurements of low star-formation efficiencies. Furthermore, most cores are found at high column densities; more than 75% of 1.1 mm cores are associated with Av>8 mag in Perseus, 15 mag in Serpens, and 20-23 mag in Ophiuchus.

Discovered in 1995 at the Caltech Submillimeter Observatory (CSO), the vibrationally-excited water maser line at 658 GHz (455 micron) is seen in oxygen-rich giant and supergiant stars. Because this maser can be so strong (up to thousands of Janskys), it was very helpful during the commissioning phase of the highest frequency band (620-700 GHz) of the Submillimeter Array (SMA) interferometer. From late 2002 to early 2006, brief attempts were made to search for emission from additional sources beyond the original CSO survey. These efforts have expanded the source count from 10 to 16. The maser emission appears to be quite compact spatially, as expected from theoretical considerations; thus these objects can potentially be used as atmospheric phase calibrators. Many of these objects also exhibit maser emission in the vibrationally-excited SiO maser at 215 GHz. Because both maser lines likely originate from a similar physical region, these objects can be used to test techniques of phase transfer calibration between millimeter and submillimeter bands. The 658 GHz masers will be important beacons to assess the performance of the Atacama Large Millimeter Array (ALMA) in this challenging high-frequency band.

Phase closure at 682 GHz and 691 GHz was first achieved using three antennas of the Submillimeter Array (SMA) interferometer located on Mauna Kea, Hawaii. Initially, phase closure was demonstrated at 682.5 GHz on Sept. 19, 2002 using an artificial ground-based "beacon" signal. Subsequently, astronomical detections of both Saturn and Uranus were made at the frequency of the CO(6-5) transition (691.473 GHz) on all three baselines on Sept. 22, 2002. While the larger planets such as Saturn are heavily resolved even on these short baselines (25.2m, 25.2m and 16.4m), phase closure was achieved on Uranus and Callisto. This was the first successful experiment to obtain phase closure in this frequency band. The CO(6-5) line was also detected towards Orion BN/KL and other Galactic sources, as was the vibrationally-excited 658 GHz water maser line toward evolved stars. We present these historic detections, as well as the first arcsecond-scale images obtained in this frequency band.

We present observations of 10.6 square degrees of the Perseus molecular cloud at 24, 70, and 160 microns with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). The image mosaics show prominent, complex extended emission dominated by illuminating B stars on the East side of the cloud, and by cold filaments of 160 micron emission on the West side. Of 3950 point sources identified at 24 microns, 1141 have 2MASS counterparts. A quarter of these populate regions of the Ks vs. Ks-[24] diagram that are distinct from stellar photospheres and background galaxies, and thus are likely to be cloud members with infrared excess. Nearly half (46%) of these 24 micron excess sources are distributed outside the IC 348 and NGC 1333 clusters. NGC 1333 shows the highest fraction of stars with flat or rising spectral energy distributions (28%), while Class II SEDs are most common in IC 348. These results are consistent with previous relative age determinations for the two clusters. The intercluster region contains several tightly clumped (r~0.1 pc) young stellar aggregates whose members exhibit a wide variety of infrared spectral energy distributions characteristic of different circumstellar environments. One possible explanation is a significant age spread among the aggregate members, such that some have had time to evolve more than others. Alternatively, if the aggregate members all formed at roughly the same time, then remarkably rapid circumstellar evolution would be required to account for the association of Class I and Class III sources at ages <~1 Myr. We highlight important results for several other objects as well (full abstract in the paper).

We present IRAC (3.6, 4.5, 5.8, and 8.0 micron) observations of the Chamaeleon II molecular cloud. The observed area covers about 1 square degree defined by $A_V >2$. Analysis of the data in the 2005 c2d catalogs reveals a small number of sources (40) with properties similar to those of young stellaror substellar objects (YSOs). The surface density of these YSO candidates is low, and contamination by background galaxies appears to be substantial, especially for sources classified as Class I or flat SED. We discuss this problem in some detail and conclude that very few of the candidate YSOs in early evolutionary stages are actually in the Cha II cloud. Using a refined set of criteria, we define a smaller, but more reliable, set of 24 YSO candidates.

We report high angular-resolution (~1") CO J=3--2 interferometric mapping, using the Submillimeter Array (SMA), of IRAS22036+5306 (I22036), a bipolar pre-planetary nebula (PPN) with knotty jets discovered in our HST SNAPshot survey of young PPNs. In addition, we have obtained supporting lower-resolution (~10") CO and 13CO J=1-0 observations with the Owens Valley Radio Observatory (OVRO) interferometer, as well as optical long-slit echelle spectra at the Palomar Observatory. The CO J=3-2 observations show the presence of a very fast (~220 km/s), highly collimated, massive (0.03 Msun) bipolar outflow with a very large scalar momentum (about 10^{39} g cm s^{-1}), and the characteristic spatio-kinematic structure of bow-shocks at the tips of this outflow. The Halpha line shows an absorption feature blue-shifted from the systemic velocity by ~100 km/s, which most likely arises in neutral interface material between the fast outflow and the dense walls of the bipolar lobes at low latitudes. The fast outflow in I22036, as in most PPNs, cannot be driven by radiation pressure. We find an unresolved source of submillimeter (and millimeter-wave) continuum emission in I22036, implying a very substantial mass (0.02-0.04 Msun) of large (radius >~1 mm), cold (< ~50 K) dust grains associated with I22036's toroidal waist. We also find that the 13C/12C ratio in I22036 is very high (0.16), close to the maximum value achieved in equilibrium CNO-nucleosynthesis (0.33). The combination of the high circumstellar mass (i.e., in the extended dust shell and the torus) and the high 13C/12C ratio in I22036 provides strong support for this object having evolved from a massive (>~4 Msun) progenitor in which hot-bottom-burning has occurred.

We present submillimeter observations of dark clouds that are part of the Spitzer Legacy Program, From Molecular Cores to Planet-Forming Disks (c2d). We used the Submillimetre Common User's Bolometer Array to map the regions observed by Spitzer by the c2d program to create a census of dense molecular cores including data from the infrared to the submillimeter. In this paper, we present the basic data from these observations: maps, fluxes, and source attributes. We also show data for an object just outside the Perseus cloud that was serendipitously observed in our program. We propose that this object is a newly discovered, evolved protostar.

We present observations of 3.86 sq. deg. of the Perseus molecular cloud complex with the Spitzer Space Telescope Infrared Array Camera (IRAC). The maps show strong extended emission arising from shocked H2 in outflows in the region and from polycyclic aromatic hydrocarbon features. More than 120,000 sources are extracted toward the cloud. Based on their IRAC colors and comparison to off-cloud and extragalactic fields, we identify 400 candidate young stellar objects. About two thirds of these are associated with the young clusters IC348 and NGC1333, while the last third is distributed over the remaining cloud. We classify the young stellar objects according to the traditional scheme based on the slope of their spectral energy distributions. Significant differences are found for the numbers of embedded Class I objects relative to more evolved Class II objects in IC348, NGC1333 and the remaining cloud with the embedded Class I and "flat spectrum" YSOs constituting 14%, 36% and 47% of the total number of YSOs identified in each of these regions. The high number of Class I objects in the extended cloud (61% of the Class I objects in the entire cloud) suggests that a significant fraction of the current star formation is occuring outside the two main clusters. Finally we discuss a number of outflows and identify their driving sources, including the known deeply embedded Class 0 sources outside the two major clusters. The Class 0 objects are found to be detected by Spitzer and have very red [3.6]-[4.5] colors but do not show similarly red [5.8]-[8.0] colors. The Class 0 objects are easily identifiable in color-color diagrams plotting these two colors but are problematic to extract automatically due to the extended emission from shocked gas or scattered light in cavities related to the associated outflows.

We present a large-scale millimeter continuum map of the Ophiuchus molecular cloud. Nearly 11 square degrees, including all of the area in the cloud with visual extinction more than 3 magnitudes, was mapped at 1.1 mm with Bolocam on the Caltech Submillimeter Observatory (CSO). By design, the map also covers the region mapped in the infrared with the Spitzer Space Telescope. We detect 44 definite sources, and a few likely sources are also seen along a filament in the eastern streamer. The map indicates that dense cores in Ophiuchus are very clustered and often found in filaments within the cloud. Most sources are round, as measured at the half power point, but elongated when measured at lower contour levels, suggesting spherical sources lying within filaments. The masses, for an assumed dust temperature of 10 K, range from 0.24 to 3.9 solar masses, with a mean value of 0.96 solar masses. The total mass in distinct cores is 42 solar masses, 0.5 to 2% of the total cloud mass, and the total mass above 4 sigma is about 80 solar masses. The mean densities in the cores are quite high, with an average of 1.6 x 10^6 per cc, suggesting short free-fall times. The core mass distribution can be fitted with a power law with slope of 2.1 plus or minus 0.3 for M>0.5 solar masses, similar to that found in other regions, but slightly shallower than that of some determinations of the local IMF. In agreement with previous studies, our survey shows that dense cores account for a very small fraction of the cloud volume and total mass. They are nearly all confined to regions with visual extinction at least 9 mag, a lower threshold than found previously.

We have completed a 1.1 mm continuum survey of 7.5 sq deg of the Perseus Molecular Cloud using Bolocam at the Caltech Submillimeter Observatory. This represents the largest millimeter or submillimeter continuum map of Perseus to date. Our map covers more than 30,000 31" (FWHM) resolution elements to a 1 sigma RMS of 15 mJy/beam. We detect a total of 122 cores above a 5 sigma point source mass detection limit of 0.18 M_sun, assuming a dust temperature of 10 K, 60 of which are new millimeter or submillimeter detections. The 1.1 mm mass function is consistent with a broken power law of slope -1.3 (0.5 M_sun<M<2.5 M_sun) and -2.6 (M>2.5 M_sun), similar to the local initial mass function slope. No more than 5% of the total cloud mass is contained in discrete 1.1 mm cores, which account for a total mass of 285 M_sun. We suggest an extinction threshold for millimeter cores of Av~5 mag, based on our calculation of the probability of finding a 1.1 mm core as a function of Av. Much of the cloud is devoid of compact millimeter emission; despite the significantly greater area covered compared to previous surveys, only 5-10 of the newly identified sources lie outside previously observed areas. The two-point correlation function confirms that dense cores in the cloud are highly structured, with significant clustering on scales as large as 2e5 AU. These 1.1 mm results, especially when combined with recently acquired c2d Spitzer Legacy data, will provide a census of dense cores and protostars in Perseus and improve our understanding of the earliest stages of star formation in molecular clouds.

The Submillimeter Array (SMA), a collaboration between the Smithsonian Astrophysical Observatory and the Academica Sinica Institute for Astronomy and Astrophysics of Taiwan, is an eight-element radio-interferometer designed to operate throughout the major atmospheric windows from about 180 to 900 GHz. In an effort to mitigate the effects of atmospheric instabilities which limit the phase coherence of the array especially in the higher frequency bands, the array was designed to allow simultaneous operation of a low frequency receiver (<350 GHz) with a high frequency receiver (>330 GHz). The overlap region of 330-350 GHz was included to facilitate dual polarization measurements in the frequency range considered to offer the highest sensitivity for continuum observations with the array. So far, the array is equipped with working SIS receivers covering the frequency ranges 176-256 GHz, 260-350 GHz, and 600-700 GHz, and single frequency operation has been routine in the lower two frequency bands for the past year. More recently, with the completion of IF hardware required to make full use of the SMA cross-correlator, dual receiver operation became possible. We have since made a number of Galactic and extra-galactic astronomical observations in dual-band mode with the hopes of using the 230 GHz receiver as a phase reference to enable improved interferometry in the 650 GHz band. We will present the current antenna and receiver performance, some of the first interferometric images in the 650 GHz receiver band, and our initial attempts at phase transfer.

We present maps of over 1.5 square degrees in Chamaeleon (Cha) II at 24, 70, and 160 micron observed with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS) and a 1.2 square degree millimeter map from SIMBA on the Swedish-ESO Submillimetre Telescope (SEST). The c2d Spitzer Legacy Team's data reduction pipeline is described in detail. Over 1500 24 micron sources and 41 70 micron sources were detected by MIPS with fluxes greater than 10-sigma. More than 40 potential YSOs are identified with a MIPS and 2MASS color-color diagram and by their spectral indices, including two previously unknown sources with 24 micron excesses. Our new SIMBA millimeter map of Cha II shows that only a small fraction of the gas is in compact structures with high column densities. The extended emission seen by MIPS is compared with previous CO observations. Some selected interesting sources, including two detected at 1 mm, associated with Cha II are discussed in detail and their SEDs presented. The classification of these sources using MIPS data is found to be consistent with previous studies.

We used the Submillimeter Array to image the SiO maser emission in the $v=1$, \$J=5-4$ transition associated with the peculiar red supergiant VY Canis Majoris. We identified seven maser components and measured their relative positions and linear polarization properties. Five of the maser components are coincident to within about 150 mas ($\sim$ 200 AU at the distance of 1.5 kpc); most of them may originate in the circumstellar envelope at a radius of about 50 mas from the star along with the SiO masers in the lowest rotational transitions. Our measurements show that two of the maser components may be offset from the inner stellar envelope (at the 3$\sigma$ level of significance) and may be part of a larger bipolar outflow associated with VY CMa identified by Shinnaga et al. The strongest maser feature at a velocity of 35.9 kms$^{-1}$ has a 60 percent linear polarization, and its polarization direction is aligned with the bipolar axis. Such a high degree of polarization suggests that maser inversion is due to radiative pumping. Five of the other maser features have significant linear polarization.

Using the partially completed Submillimeter Array with five antennas, we have observed the CO J=2-1 and 3-2 emission from the envelope surrounding the carbon star V Hya. The high-angular resolution (2"-4") maps show that V Hya is powering a bipolar molecular jet having an extreme velocity of 70-185 km/s. The axis of this high velocity jet is perpendicular to the major axis of the flattened disk-like envelope, which is expanding with a velocity of ~16 km/s.There is a third kinematic component, a medium-velocity wind having a deprojected velocity of 40-120 km/s moving along the disk plane. Both the high velocity jet and the medium velocity wind have a dynamical time scale of a few hundred years. The flattened structure and the collimated jet observed in V Hya suggests that the formation of asymmetrical structure proceeds while the central star is still in the AGB phase

We present maps of the 6 cm and 1.3 mm transitions of formaldehyde toward three cold, dense pre-protostellar cores: L1498, L1512, and L1544. The 6 cm transition is a unique probe of high density gas at low temperature. However, our models unequivocally indicate that H2CO is depleted in the interiors of PPCs, and depletion significantly affects how H2CO probes the earliest stages of star formation. Multi-stage, self-consistent models, including gas--dust energetics, of both H2CO transitions are presented, and the implications of the results are discussed.

We present imaging observations of the evolved star IRC+10216 in the CS J=14--13 line at 685.4 GHz and associated submillimeter continuum at about 2'' resolution made with the partially constructed Submillimeter Array. The CS J=14--13 line emission from the stellar envelope is well resolved both spatially and spectrally. The strong central concentration of the line emission provides direct evidence that CS is a parent molecule that forms close to the stellar photosphere, in accord with previous images of the lower excitation CS J=2--1 line and inferences from unresolved observations of vibrationally excited transitions. The continuum emission is dominated by a compact, unresolved component, consistent with the photospheric emission, that accounts for about 20% of the broadband 450 micron flux. These are the first interferometer imaging observations made in the semi-transparent 450 micron atmospheric window.

We have mapped 63 regions forming high-mass stars in CS J=5-4 using the CSO. The CS peak position was observed in C34S J=5-4 towards 57 cores and in 13CS J=5-4 towards the 9 brightest cores. The sample is a subset of a sample originally selected toward water masers; the selection on maser sources should favor sources in an early stage of evolution. The integrated intensity of CS J=5-4 correlates very well with the dust continuum emission at 350 microns. The distributions of size, virial mass, surface density, and luminosity are all peaked with a few cores skewed towards much larger values than the mean. We find a weak correlation between C34S linewidth and size, consistent with Dv ~ R^{0.3}. The linewidths are much higher than would be predicted by the usual relations between linewidth and size determined from regions of lower mass. These regions are very turbulent. The derived virial mass agrees within a factor of 2 to 3 with mass estimates from dust emission at 350 microns after corrections for the density structure are accounted for. The resulting cumulative mass spectrum of cores above 1000 solar masses can be approximated by a power law with a slope of about -0.9, steeper than that of clouds measured with tracers of lower density gas and close to that for the total masses of stars in OB associations. The median turbulent pressures are comparable to those in UCHII regions, and the pressures at small radii are similar to those in hypercompact-HII regions (P/k ~ 10^{10} K cm^{-3}). The filling factors for dense gas are substantial, and the median abundance of CS is about 10^{-9}. The ratio of bolometric luminosity to virial mass is much higher than the value found for molecular clouds as a whole, and the correlation of luminosity with mass is tighter. (Abridged).

We report the results from recent observations of Sgr A* at short-/sub-millimeter wavelengths made with the partially finished Sub-Millimeter Array (SMA) on Mauna Kea. A total of 25 epochs of observations were carried out over the past 15 months in 2001 March to 2002 May. Noticeable variations in flux density at 1.3 mm were observed showing three ``flares''. The SMA observations suggest that Sgr A* highly increases towards submillimeter wavelengths during a flare suggesting the presence of a break wavelength in spectral index around 3 mm. A cross-correlation of the SMA data at 1 mm with the VLA data at 1 cm show a global delay of $t_{delay} > 3d$, suggesting that sub-millimeter wavelengths tend to peak first. Only marginal day-to-day variations in flux density (2-3 $\sigma$) have been detected at 1.3 mm. No significant flares on a short time scale ($\sim1$ hr) have been observed at 1.3 mm. We also failed to detect significant periodic signals at a level of 5% (3$\sigma$) from Sgr A* in a periodic searching window ranging from 10 min to 2.5 hr. The flares observed at the wavelengths between short-centimeter and sub-millimeter might be a result of collective mass ejections associated with X-ray flares that originate from the inner region of the accretion disk near the supermassive black hole.

The GLAST Balloon Flight Engineering Model (BFEM) represents one of 16 towers that constitute the Large Area Telescope (LAT), a high-energy (>20 MeV) gamma-ray pair-production telescope being built by an international partnership of astrophysicists and particle physicists for a satellite launch in 2006. The prototype tower consists of a Pb/Si pair-conversion tracker (TKR), a CsI hodoscopic calorimeter (CAL), an anti-coincidence detector (ACD) and an autonomous data acquisition system (DAQ). The self-triggering capabilities and performance of the detector elements have been previously characterized using positron, photon and hadron beams. External target scintillators were placed above the instrument to act as sources of hadronic showers. This paper provides a comprehensive description of the BFEM data-reduction process, from receipt of the flight data from telemetry through event reconstruction and background rejection cuts. The goals of the ground analysis presented here are to verify the functioning of the instrument and to validate the reconstruction software and the background-rejection scheme.

Atomic carbon has been detected in the envelopes of three carbon-rich evolved stars: HD 44179 (=AFGL 915, the `Red Rectangle'); HD 56126; and, tentatively, the carbon star V Hya. This brings to seven the number of evolved star envelopes in which CI has been detected. Upper limits were found for several other stars, including R CrB. CI was not detected in several oxygen rich post asymptotic giant branch (AGB) stars (OH231.8+4.2, for example), although it is detected in their carbon-rich analogues. Two trends are evident in the data. First, circumstellar envelopes with detectable CI are overwhelmingly carbon rich, suggesting that much of the CI is produced by the dissociation of molecules other than CO. Second, the more evolved the envelope away from the AGB, the higher is the CI/CO ratio. The oxygen-rich supergiant star Alpha Ori remains the only oxygen rich star with a wind containing detectable CI. These data suggest an evolutionary sequence for the CI/CO ratio in cool circumstellar envelopes. This ratio is small (a few %) while the star is on the AGB, and the CI is located in the outer envelope and produced by photodissociation. The ratio increases to about 0.5 as the star evolves away from the AGB because of the dissociation of CO and other carbon-bearing molecules by shocks caused by the fast winds which appear at the end of evolution on the AGB. Finally, the ratio becomes >> 1 as the central star becomes hot enough to photodissociate CO.

V Hya, an evolved carbon star with a complex circumstellar envelope, has two variability periods, 530d and 6000d (17 years). We analyze recent light curve data and show that both variations have benn present for at least 100 years and have been regular over this time. The 530d period and its 1.5 - 2 magnitude amplitude show that V Hya is a Mira variable. We suggest that the star is in a binary system (as also suspected from the structure of the circumstellar envelope) and that the 17-year variation is due to extinction by circumstellar dust orbiting with the companion. The properties of the envelope found from molecular line observations: the fast molecular wind, the relatively small size of the dense circumstellar envelope, and the high mass loss rate, all suggest that V Hya has entered its `superwind' phase. However, its spectral type, period, colors, and lack of ionizing radiation show that the star is still on the AGB. These properties add to the evidence that the complex structures of many planetary planetary nebulae, including fast stellar winds, originate during the final phases of mass loss on the AGB.

CO(J = 2-1) and SiO(J = 5-4) emission has been observed from the molecular envelope around the nearby S star pi Gru. The CO line profile differs from the usual parabolic shape seen in uniformly expanding envelopes; it has a Voigt-like profile and two horns. A model for line formation in the envelope shows that a tilted, expanding disk reproduces the observations well. The star also has a fast molecular wind, with a projected outflow speed of at least 70 and perhaps as high as 90 km/s. The fast wind is presumably ejected from the poles of the disk. These observations show that the complex structure seen in many planetary nebulae, including quadrupolar structure and fast winds, may largely evolve from structure formed while the progenitor star is in the last stages of evolution on the AGB.

This paper describes observations of a new phenomenon in evolved mass-losing AGB stars: the presence of two winds with different expansion velocities. CO(2-1) and CO(3-2) line emission was observed for 45 AGB stars at high velocity resolution and double winds found in 20% of the sample. Highly asymmetric lines were found in six other stars. The data tentatively suggest that double winds occur when the star undergoes a change (pulsational mode, chemical composition) and that the very narrow components represent the onset of a new phase of mass loss.

We have detected the J=4-3 rotational transition of 12CO in absorption at z=0.89 towards the quasar PKS1830-211, but not the 12CO 5-4 or the 3P1-3P0 fine structure line of neutral carbon. The intervening molecular medium thus has a total 12CO column density of N(CO)~2x10^18 cm^-2, which corresponds to the large column density of molecular hydrogen of N(H2)=2.5x10^22 cm^-2 and a reddening of Av=25 magnitudes. The 12CO excitation temperature is low, below 15 K. Comparison with the molecular absorption results of Wiklind and Combes (1996) shows that the absorbing material has similar molecular abundances to Galactic dark clouds. We find an upper limit for atomic carbon of N(CI)<10^18 cm^-2, which again would be the case for most Galactic dark clouds. We also report new observations of the absorbing system towards B0218+357 at z=0.68. We have tentatively detected the 13CO 4-3 line, but for H2O, although a feature is seen at the correct velocity, due to the inadequate signal to noise ratio we report only an upper limit for the fundamental line of ortho water vapor. The tentative detection of the 13CO J=4-3 line implies that the 13CO excitation temperature is lower than 20K and the column density is fairly large, N(13CO)~10^17 cm^-2, giving rise to saturated absorption in the J=2-1 transition. The total column density of molecular gas is again large in this source, N(H2)>2x10^22 cm^-2, corresponding to a reddening >20 magnitudes.

The 609u 3P1 -> 3P0 line of neutral atomic carbon has been detected in the protoplanetary nebulae CRL 618 and CRL 2688. CI appears to be about 7% as abundant as CO in CRL 2688, and 70% as abundant as CO in CRL 618. CI emission arises primarily from the slow component of the wind from these objects. The data presented here indicate that CI gradually becomes more abundant as the star leaves the AGB, and CI is significantly enhanced relative to CO before the star develops a large HII region.

We report the results of a search for neutrino-induced particle cascades using a deep ocean water Cherenkov detector. The effective mass of the detector, a string of seven 40 cm diameter photomultipliers at 5.2 m spacing, is found through simulation analysis to be surprisingly large: greater than 1 megaton of water at incident neutrino energies of 1 PeV. We find no evidence for neutrino-induced cascades in 18.6 hours of observation. Although the limit implied by this observation is the strongest yet for predictions of active galatic nuclei (AGN) neutrinos at energies above 100 TeV, perhaps the more intriguing result is that the power of these techniques can be exploited to test these AGN models in a relatively short time.

We report detection of the 609u line of neutral atomic carbon in globules of the Helix nebula. The measurements were made towards the position of peak CO emission. At the same position, we obtained high-quality CO(2-1) and 13CO(2-1) spectra and a 135" x 135" map in CO(2-1). The velocity distribution of CI shows six narrow (1 -> 2 km/sec) components which are associated with individual globules traced in CO. The CI column densities are 0.5 -> 1.2 x 10^16/cm^2. CI is found to be a factor of ~6 more abundant than CO. Our estimate for the mass of the neutral envelope is an order of magnitude larger than previous estimates. The large abundance of CI in the Helix can be understood as a result of the gradual photoionisation of the molecular envelope by the central star's radiation field.

A survey of CO(3-2) emission from optically visible oxygen-rich Mira variable stars within 500 pc of the sun was conducted. A molecular envelope was detected surrounding 36 of the 66 stars examined. Some of these stars have lower outflow velocities than any Miras previously detected in CO. The average terminal velocity of the ejected material was 7.0 km/sec, about half the value found in Miras selected by infrared criteria. None of the stars with spectral types earlier than M 5.5 were detected. The terminal velocity increases as the temperature of the stellar photosphere decreases, as would be expected for a radiation driven wind. Mass loss rates for the detected objects were calculated, and it was found that there is no correlation between the infrared color of a Mira variable, and its mass loss rate. The mass loss rate is correlated with the far infrared luminosity, although a few stars appear to have extensive dust envelopes without any detectable molecular wind. A power-law relationship is found to hold between the mass loss rate and the terminal velocity of the ejected material. This relationship indicates that the dust envelope should be optically thick in the near infrared and visible regions of the spectrum when the outflow velocity is > 17 km/sec. At the low end of the range of outflow velocities seen, the dust drift velocity may be high enough to lead to the destruction of the grains via sputtering. Half of the stars which were detected were re-observed in the CO(4-3) transition. A comparison of the outflow velocities obtained from these observations with those obtained by other investigators at lower frequencies shows no evidence for gradual acceleration of the outer molecular envelope.

The DUMAND-II detector will search for astronomical sources of high energy neutrinos. Successful deployment of the basic infrastructure, including the shore cable, the underwater junction box, and an environmental module was accomplished in December, 1993. One optical module string was also deployed and operated, logging data for about 10 hours. The underwater cable was connected to the shore station where we were able to successfully exercise system controls and log further environmental data. After this time, water leaking into the electronics control module for the deployed string disabled the string electrical system. The acquired data are consistent with the expected rate of downgoing muons, and our ability to reconstruct muons was demonstrated. The measured acoustical backgrounds are consistent with expectation, which should allow acoustical detection of nearby PeV particle cascades. The disabled string has been recovered and is undergoing repairs ashore. We have identified the source of the water leak and implemented additional testing and QC procedures to ensure no repetition in our next deployment. We will be ready to deploy three strings and begin continuous data taking in 1995.

We present maps at $20''$ resolution of the molecular emission around 18 evolved stars (14~asymptotic giant branch stars, one supergiant, two proto-planetary nebulae and one planetary nebula), mostly in the $^{12}$CO(3-2) line. Almost all molecular envelopes appear to be at least marginally resolved at this resolution. A substantial fraction of the molecular envelopes show clear deviations from spherical symmetry in the form of elliptical or bipolar envelopes. This indicates that there is a need to implement non-spherical mass loss in current scenarios of the late stages of stellar evolution, in particular on the asymptotic giant branch.