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.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

The Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.

This paper outlines the science case for line-intensity mapping with a space-borne instrument targeting the sub-millimeter (microwaves) to the far-infrared (FIR) wavelength range. Our goal is to observe and characterize the large-scale structure in the Universe from present times to the high redshift Epoch of Reionization. This is essential to constrain the cosmology of our Universe and form a better understanding of various mechanisms that drive galaxy formation and evolution. We argue that the proposed frequency range would make it possible to probe important metal cooling lines such as [CII] up to very high redshift as well as a large number of rotational lines of the CO molecule. These can be used to trace molecular gas and dust evolution and constrain the buildup in both the cosmic star formation rate density and the cosmic infrared background (CIB). Moreover, surveys at the highest frequencies will detect FIR lines which are used as diagnostics of galaxies and AGN. Tomography of these lines over a wide redshift range will enable invaluable measurements of the cosmic expansion history at epochs inaccessible to other methods, competitive constraints on the parameters of the standard model of cosmology, and numerous tests of dark matter, dark energy, modified gravity and inflation. To reach these goals, large-scale structure must be mapped over a wide range in frequency to trace its time evolution over a reasonable fraction of the volume of the observable Universe. In addition, the surveyed area needs to be very large to beat cosmic variance and to probe the largest scales where its easier to separate the astrophysical and cosmological contributions to the observed signal. Only, a space-borne mission can properly meet these requirements.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

Constraints on primordial non-Gaussianity (PNG) will shed light on the origin of primordial fluctuations and the physics of the early universe. The intensity mapping technique is a promising probe of structure formation on large scales; at high redshifts, it can provide complementary information to other cosmological probes. We explore the potential of future wide-field [CII] and CO intensity mapping surveys in constraining PNG of the local shape, which induces a distinct, scale-dependent correction to the line bias. We explore the parameter space of CO and [CII] survey designs that can achieve the nominal target of $\sigma(f_{\rm NL}^{\rm loc}) =1$ and further calculate what constraints such surveys can place on PNG of the equilateral, orthogonal, and quasi-single field shapes. We further test the dependence of these constraints on various model assumptions; namely the halo mass function, modeling of line bias, and correlation between mass and [CII]/CO luminosity. We find that the ability of CO and [CII] intensity mapping surveys to constrain $f_{\rm NL}^{\rm loc}$ relies heavily on the spectral (redshift) coverage, requiring at least a coverage of an octave in frequency, to produce significant results, with the optimized surveys generally covering between $z\sim[2-8]$. As this redshift window partially overlaps with 21-cm EoR experiments like the Hydrogen Epoch of Reionization Array (HERA), we additionally explore the prospects for a hypothetical set of surveys to place constraints on $\sigma(f_{\rm NL}^{\rm loc})$ utilizing multi-tracer analysis between the lines of HI, CO, and [CII].

Primordial non-Gaussianity is an invaluable window into the physical processes that gave rise to the cosmological structure. The presence of local shape PNG imprints a distinct scale-dependent correction to the bias of dark matter tracers on large scales, which can be effectively probed via the technique of intensity mapping. Considering an upcoming generation of experiments, we demonstrate that intensity mapping of CO and [CII] emission can improve upon the current best constraints from the {\it Planck} satellite. We show that measurement of the CO intensity power spectrum by a hypothetical next stage of the ground-based COMAP experiment can achieve $\sigma(f_{\rm NL}^{\rm loc}) = 3.4$, and that the proposed CMB satellite mission PIXIE can achieve $\sigma(f_{\rm NL}^{\rm loc}) = 3.9$ via measurement of [CII] intensity power spectrum.

We present a new set of measurements obtained with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) of Jupiter's microwave thermal emission near the 1.3 cm ammonia (NH$_{3}$) absorption band. We use these observations to investigate the ammonia mole fraction in the upper troposphere, near $0.3 < P < 2$ bar, based on radiative transfer modeling. We find that the ammonia mole fraction must be $\sim2.4\times 10^{-4}$ below the NH$_{3}$ ice cloud, i.e., at $0.8 < P < 8$ bar, in agreement with results by de Pater et al. (2001, 2016a). We find the NH$_{3}$ cloud-forming region between $0.3 < P < 0.8$ bar to be globally sub-saturated by 55% on average, in accordance with the result in Gibson et al. (2005). Although our data are not very sensitive to the region above the cloud layer, we are able to set an upper limit of $2.4\times 10^{-7}$ on the mole fraction here, a factor of $\sim$10 above the saturated vapor curve.

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

We present a measurement of the abundance of carbon monoxide in the early Universe, utilizing the final results from the CO Power Spectrum Survey (COPSS). Between 2013 and 2015, we performed observations with the Sunyaev-Zel'dovich Array (SZA) to measure the aggregate CO emission from $z\sim3$ galaxies with the intensity mapping technique. Data were collected on 19 fields, covering an area of 0.7 square degrees, over a frequency range of 27 to 35 GHz. With these data, along with data analyzed in COPSS I, we are able to observe the CO(1-0) transition in the redshift range of $z=2.3-3.3$ for spatial frequencies between $k=0.5-10\ h\,\textrm{Mpc}^{-1}$, spanning a comoving volume of $4.9\times10^{6}\ h^{-3}\,\textrm{Mpc}^{3}$. We present estimates of contributions from continuum sources and ground illumination within our measurement. We constrain the amplitude of the CO power spectrum to be $P_{\textrm{CO}}=3.0_{-1.3}^{+1.3}\times10^{3}\ \mu\textrm{K}^{2} (h^{-1}\,\textrm{Mpc})^{3}$, or $\Delta^{2}_{\textrm{CO}}(k\!=\!1\ h\,\textrm{Mpc}^{-1})=1.5^{+0.7}_{-0.7} \times10^{3}\ \mu\textrm{K}^{2}$, at 68% confidence, and $P_{\textrm{CO}}>0$ at 98.9% confidence. These results are a factor of 10 improvement in sensitivity compared to the previous limits set in COPSS I. We use this measurement to place constraints on the CO(1-0) galaxy luminosity function at $z\sim3$. We constrain the ratio of CO(1-0) luminosity to host halo mass to $A_{\textrm{CO}}=6.3_{-2.1}^{+1.4}\times10^{-7}\ L_{\odot}\ M_{\odot}^{-1}$, and estimate a mass fraction of molecular gas of $f_{\textrm{H}_{2}}=5.5^{+3.4}_{-2.2}\times10^{-2}$ for halos with masses of order $10^{12}M_{\odot}$. Using simple theoretical estimates for the scaling of molecular gas mass fraction and halo mass, we estimate the global density of molecular gas to be $\rho_{z\sim3}(\textrm{H}_{2})=1.1_{-0.4}^{+0.7}\times10^{8}\ M_{\odot}\ \textrm{Mpc}^{-3}$.

We present constraints on the abundance of carbon-monoxide in the early Universe from the CO Power Spectrum Survey (COPSS). We utilize a data set collected between 2005 and 2008 using the Sunyaev-Zel'dovich Array (SZA), which were previously used to measure arcminute-scale fluctuations of the CMB. This data set features observations of 44 fields, covering an effective area of 1.7 square degrees, over a frequency range of 27 to 35 GHz. Using the technique of intensity mapping, we are able to probe the CO(1-0) transition, with sensitivity to spatial modes between $k=0.5{-}2\ h\,\textrm{Mpc}^{-1}$ over a range in redshift of $z=2.3{-}3.3$, spanning a comoving volume of $3.6\times10^{6}\ h^{-3}\,\textrm{Mpc}^{3}$. We demonstrate our ability to mitigate foregrounds, and present estimates of the impact of continuum sources on our measurement. We constrain the CO power spectrum to $P_{\textrm{CO}}<2.6\times10^{4}\ \mu\textrm{K}^{2} (h^{-1}\,\textrm{Mpc})^{3}$, or $\Delta^{2}_{\textrm{CO}}(k\! = \! 1 \ h\,\textrm{Mpc}^{-1})<1.3 \times10^{3}\ \mu\textrm{K}^{2}$, at $95\%$ confidence. This limit resides near optimistic predictions for the CO power spectrum. Under the assumption that CO emission is proportional to halo mass during bursts of active star formation, this corresponds to a limit on the ratio of $\textrm{CO}(1{-}0)$ luminosity to host halo mass of $A_{\textrm{CO}}<1.2\times10^{-5}\ L_{\odot}\ M_{\odot}^{-1}$. Further assuming a Milky Way-like conversion factor between CO luminosity and molecular gas mass ($\alpha_{\textrm{CO}}=4.3\ M_{\odot}\ (\textrm{K}\ \textrm{km}\ \textrm{s}^{-1}\ \textrm{pc}^{-2})^{-1}$), we constrain the global density of molecular gas to $\rho_{z\sim3}(M_{\textrm{H}_{2}})\leq 2.8 \times10^{8}\ M_{\odot}\ \textrm{Mpc}^{-3}$.

Searches for slow radio transients and variables have generally focused on extragalactic populations, and the basic parameters of Galactic populations remain poorly characterized. We present a large 3 GHz survey performed with the Allen Telescope Array (ATA) that aims to improve this situation: ASGARD, the ATA Survey of Galactic Radio Dynamism. ASGARD observations spanned 2 years with weekly visits to 23 deg^2 in two fields in the Galactic Plane, totaling 900 hr of integration time on science fields and making it significantly larger than previous efforts. The typical blind unresolved source detection limit was 10 mJy. We describe the observations and data analysis techniques in detail, demonstrating our ability to create accurate wide-field images while effectively modeling and subtracting large-scale radio emission, allowing standard transient-and-variability analysis techniques to be used. We present early results from the analysis of two pointings: one centered on the microquasar Cygnus X-3 and one overlapping the Kepler field of view (l = 76{\deg}, b = +13.5{\deg}). Our results include images, catalog statistics, completeness functions, variability measurements, and a transient search. Out of 134 sources detected in these pointings, the only compellingly variable one is Cygnus X-3, and no transients are detected. We estimate number counts for potential Galactic radio transients and compare our current limits to previous work and our projection for the fully-analyzed ASGARD dataset.

The Allen Telescope Array (ATA) is a cm-wave interferometer in California, comprising 42 antenna elements with 6-m diameter dishes. We characterize the antenna optical accuracy using two-antenna interferometry and radio holography. The distortion of each telescope relative to the average is small, with RMS differences of 1 percent of beam peak value. Holography provides images of dish illumination pattern, allowing characterization of as-built mirror surfaces. The ATA dishes can experience mm-scale distortions across -2 meter lengths due to mounting stresses or solar radiation. Experimental RMS errors are 0.7 mm at night and 3 mm under worst case solar illumination. For frequencies 4, 10, and 15 GHz, the nighttime values indicate sensitivity losses of 1, 10 and 20 percent, respectively. The ATA.s exceptional wide-bandwidth permits observations over a continuous range 0.5 to 11.2 GHz, and future retrofits may increase this range to 15 GHz. Beam patterns show a slowly varying focus frequency dependence. We probe the antenna optical gain and beam pattern stability as a function of focus and observation frequency, concluding that ATA can produce high fidelity images over a decade of simultaneous observation frequencies. In the day, the antenna sensitivity and pointing accuracy are affected. We find that at frequencies greater than 5 GHz, daytime observations greater than 5 GHz will suffer some sensitivity loss and it may be necessary to make antenna pointing corrections on a 1 to 2 hourly basis.

We present the first observations of the Sun with the Allen Telescope Array (ATA). We used up to six frequencies, from 1.43 to 6 GHz, and baselines from 6 to 300 m. To our knowledge, these are the first simultaneous multifrequency full-Sun maps obtained at microwave frequencies without mosaicing. The observations took place when the Sun was relatively quiet, although at least one active region was present each time. We present multi-frequency flux budgets for each sources on the Sun. Outside of active regions, assuming optically thin bremsstrahlung (free--free) coronal emission on top of an optically thick ~10 000 K chromosphere, the multi-frequency information can be condensed into a single, frequency-independent, "coronal bremsstrahlung contribution function" [EM/sqrt(T)] map. This technique allows the separation of the physics of emission as well as a measurement of the density structure of the corona. Deviations from this simple relationship usually indicate the presence of an additional gyroresonance-emission component, as is typical in active regions.

We present results from daily radio continuum observations of the Bootes field as part of the Pi GHz Sky Survey (PiGSS). These results are part of a systematic and unbiased campaign to characterize variable and transient sources in the radio sky. The observations include 78 individual epochs distributed over 5 months at a radio frequency of 3.1 GHz with a median RMS image noise in each epoch of 2.8 mJy. We produce 5 monthly images with a median RMS of 0.6 mJy. No transient radio sources are detected in the daily or monthly images. At 15 mJy, we set an upper limit (2 sigma) to the surface density of 1-day radio transients at 0.025 deg^-2. At 5 mJy, we set an upper limit (2\sigma) to the surface density of 1-month radio transients at 0.18 deg^-2. We also produce light curves for 425 sources and explore the variability properties of these sources. Approximately 20% of the sources exhibit some variability on daily and monthly time scales. The maximum RMS fractional modulations on the 1 day and 1 month time scales for sources brighter than 10 mJy are 2 and 0.5, respectively. The probability of a daily fluctuation for all sources and all epochs by a factor of 10 is less than 10^-4. We compare the radio to mid-infrared variability for sources in the field and find no correlation. Finally, we apply the statistics of transient and variable populations to constrain models for a variety of source classes.

We present our second paper on the Allen Telescope Array Twenty-centimeter Survey (ATATS), a multi-epoch, ~700 sq. deg. radio image and catalog at 1.4 GHz. The survey is designed to detect rare, bright transients as well as to commission the ATA's wide-field survey capabilities. ATATS explores the challenges of multi-epoch transient and variable source surveys in the domain of dynamic range limits and changing (u,v) coverage. Here we present images made using data from the individual epochs, as well as a revised image combining data from all ATATS epochs. The combined image has RMS noise 3.96 mJy / beam, with a circular beam of 150 arcsec FWHM. The catalog, generated using a false detection rate algorithm, contains 4984 sources, and is >90% complete to 37.9 mJy. The catalogs generated from snapshot images of the individual epochs contain between 1170 and 2019 sources over the 564 sq. deg. area in common to all epochs. The 90% completeness limits of the single epoch catalogs range from 98.6 to 232 mJy. We compare the catalog generated from the combined image to those from individual epochs, and from the NRAO VLA Sky Survey (NVSS), a legacy survey at the same frequency. We are able to place new constraints on the transient population: fewer than 6e-4 transients / sq. deg., for transients brighter than 350 mJy with characteristic timescales of minutes to days. This strongly rules out an astronomical origin for the ~1 Jy sources reported by Matsumura et al. (2009), based on their stated rate of 3.1e-3 / sq. deg.

We have observed 37 bright, polarized radio sources with the Allen Telescope Array (ATA) to present a novel analysis of their Faraday rotation properties. Each source was observed during the commissioning phase with 2 to 4 100-MHz bands at frequencies ranging from 1 to 2 GHz. These observations demonstrate how the continuous frequency coverage of the ATA's log-periodic receiver can be applied to the study of Faraday rotation measures (RMs). We use RM synthesis to show that wide-bandwidth data can find multiple RM components toward a single source. Roughly a quarter of the sources studied have extra RM components with high confidence (brighter than ~40 mJy), when observing with a RM resolution of roughly 100 rad/m2. These extra components contribute 10%-70% of the total polarized flux. This is the first time multiple RM components have been identified in a large sample of point sources. For our observing configuration, these extra RM components bias the measurement of the peak RM by 10-15 rad/m2 ; more generally, the peak RM cannot be determined more precisely than the RM beam size. Comparing our 1-2 GHz RM spectra to VLBA polarimetric maps shows both techniques can identify complicated Faraday structures in the sources. However, the RM values and fractional polarization are generally smaller at lower frequencies than in the higher-frequency VLBA maps. With a few exceptions, the RMs from this work are consistent with that of earlier, narrow-bandwidth, all-sky surveys. This work also describes the polarimetry calibration procedure and that on-axis ATA observations of linear polarization can be calibrated to an accuracy of 0.2% of Stokes I. Future research directions include studying the time-dependent RM structure in Active Galactic Nuclei (AGNs) and enabling accurate, wide-area RM surveys to test models of Galactic and extragalactic magnetic fields.

The Allen Telescope Array (ATA) is a Large-Number-Small-Diameter radio telescope array currently with 42 individual antennas and 5 independent back-end science systems (2 imaging FX correlators and 3 time domain beam formers) located at the Hat Creek Radio Observatory (HCRO). The goal of the ATA is to run multiple back-ends simultaneously, supporting multiple science projects commensally. The primary software control systems are based on a combination of Java, JRuby and Ruby on Rails. The primary control API is simplified to provide easy integration with new back-end systems while the lower layers of the software stack are handled by a master observing system. Scheduling observations for the ATA is based on finding a union between the science needs of multiple projects and automatically determining an efficient path to operating the various sub-components to meet those needs. When completed, the ATA is expected to be a world-class radio telescope, combining dedicated SETI projects with numerous radio astronomy science projects.

The Pi GHz Sky Survey (PiGSS) is a key project of the Allen Telescope Array. PiGSS is a 3.1 GHz survey of radio continuum emission in the extragalactic sky with an emphasis on synoptic observations that measure the static and time-variable properties of the sky. During the 2.5-year campaign, PiGSS will twice observe ~250,000 radio sources in the 10,000 deg^2 region of the sky with b > 30 deg to an rms sensitivity of ~1 mJy. Additionally, sub-regions of the sky will be observed multiple times to characterize variability on time scales of days to years. We present here observations of a 10 deg^2 region in the Bootes constellation overlapping the NOAO Deep Wide Field Survey field. The PiGSS image was constructed from 75 daily observations distributed over a 4-month period and has an rms flux density between 200 and 250 microJy. This represents a deeper image by a factor of 4 to 8 than we will achieve over the entire 10,000 deg^2. We provide flux densities, source sizes, and spectral indices for the 425 sources detected in the image. We identify ~100$ new flat spectrum radio sources; we project that when completed PiGSS will identify 10^4 flat spectrum sources. We identify one source that is a possible transient radio source. This survey provides new limits on faint radio transients and variables with characteristic durations of months.

We present the Allen Telescope Array Twenty-centimeter Survey (ATATS), a multi-epoch (12 visits), 690 square degree radio image and catalog at 1.4GHz. The survey is designed to detect rare, very bright transients as well as to verify the capabilities of the ATA to form large mosaics. The combined image using data from all 12 ATATS epochs has RMS noise sigma = 3.94mJy / beam and dynamic range 180, with a circular beam of 150 arcsec FWHM. It contains 4408 sources to a limiting sensitivity of S = 20 mJy / beam. We compare the catalog generated from this 12-epoch combined image to the NRAO VLA Sky Survey (NVSS), a legacy survey at the same frequency, and find that we can measure source positions to better than ~20 arcsec. For sources above the ATATS completeness limit, the median flux density is 97% of the median value for matched NVSS sources, indicative of an accurate overall flux calibration. We examine the effects of source confusion due to the effects of differing resolution between ATATS and NVSS on our ability to compare flux densities. We detect no transients at flux densities greater than 40 mJy in comparison with NVSS, and place a 2-sigma upper limit on the transient rate for such sources of 0.004 per square degree. These results suggest that the > 1 Jy transients reported by Matsumura et al. (2009) may not be true transients, but rather variable sources at their flux density threshold.

The first 42 elements of the Allen Telescope Array (ATA-42) are beginning to deliver data at the Hat Creek Radio Observatory in Northern California. Scientists and engineers are actively exploiting all of the flexibility designed into this innovative instrument for simultaneously conducting surveys of the astrophysical sky and conducting searches for distant technological civilizations. This paper summarizes the design elements of the ATA, the cost savings made possible by the use of COTS components, and the cost/performance trades that eventually enabled this first snapshot radio camera. The fundamental scientific program of this new telescope is varied and exciting; some of the first astronomical results will be discussed.