Search

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a balloon-borne telescope designed to survey star formation over cosmological time scales using intensity mapping in the 420 - 540 GHz frequency range. EXCLAIM uses a fully cryogenic telescope coupled to six on-chip spectrometers featuring kinetic inductance detectors (KIDs) to achieve high sensitivity, allowing for fast integration in dark atmospheric windows. The telescope receiver is cooled to $\approx$ 1.7 K by immersion in a superfluid helium bath and enclosed in a superfluid-tight shell with a meta-material anti-reflection coated silicon window. In addition to the optics and the spectrometer package, the receiver contains the magnetic shielding, the cryogenic segment of the spectrometer readout, and the sub-Kelvin cooling system. A three-stage continuous adiabatic demagnetization refrigerator (CADR) keeps the detectors at 100 mK while a $^4$He sorption cooler provides a 900 mK thermal intercept for mechanical suspensions and coaxial cables. We present the design of the EXCLAIM receiver and report on the flight-like testing of major receiver components, including the superfluid-tight receiver window and the sub-Kelvin coolers., Comment: 13 pages, 7 figures, 1 table

This paper presents the design, calibration, and survey strategy of the Fast Radio Burst (FRB) digital backend and its real-time data processing pipeline employed in the Tianlai Cylinder Pathfinder array. The array, consisting of three parallel cylindrical reflectors and equipped with 96 dual-polarization feeds, is a radio interferometer array designed for conducting drift scans of the northern celestial semi-sphere. The FRB digital backend enables the formation of 96 digital beams, effectively covering an area of approximately 40 square degrees with 3 dB beam. Our pipeline demonstrates the capability to make automatic search of FRBs, detecting at quasi-real-time and classify FRB candidates automatically. The current FRB searching pipeline has an overall recall rate of 88\%. During the commissioning phase, we successfully detected signals emitted by four well-known pulsars: PSR B0329+54, B2021+51, B0823+26, and B2020+28. We report the first discovery of an FRB by our array, designated as FRB 20220414A. We also investigate the optimal arrangement for the digitally formed beams to achieve maximum detection rate by numerical simulation., Comment: 27 pages, 21 figures, 7 tables, RAA accepted

The visibilities measured by radio astronomical interferometers include non-astronomical correlated signals that arise from the local environment of the array. These correlated signals are especially important in compact arrays such as those under development for 21\,cm intensity mapping. The amplitudes of the contaminated visibilities can exceed the expected 21\,cm signal and represent a significant systematic effect. We study the receiver noise radiated by antennas in compact arrays and develop a model for how it couples to other antennas. We apply the model to the Tianlai Dish Pathfinder Array (TDPA), a compact array of 16, 6-m dish antennas. The coupling model includes electromagnetic simulations, measurements with a network analyzer, and measurements of the noise of the receivers. We compare the model to drift-scan observations with the array and set requirements on the level of antenna cross-coupling for 21\,cm intensity mapping instruments. We find that for the TDPA, cross-coupling would have to be reduced by TBD orders of magnitude in order to contribute negligibly to the visibilities., Comment: 20 pages, 22 figures, accepted for publication by JAI

This paper introduces a technique called NKL, which cleans both polarized and unpolarized foregrounds from HI intensity maps by applying a Karhunen-Lo\`eve transform on the needlet coefficients. In NKL, one takes advantage of correlations not only along the line of sight, but also between different angular regions, referred to as ``chunks". This provides a distinct advantage over many of the standard techniques applied to map-space that one finds in the literature, which do not consider such spatial correlations. Moreover, the NKL technique does not require any priors on the nature of the foregrounds, which is important when considering polarized foregrounds. We also introduce a modified version of GNILC, referred to as MGNILC, which incorporates an approximation of the foregrounds to improve performance. The NKL and MGNILC techniques are tested on simulated maps which include polarized foregrounds. Their performance is compared to the GNILC, GMCA, ICA and PCA techniques. Two separate tests were performed. One at $1.84 < z < 2.55$ and the other at $0.31 < z < 0.45$. NKL was found to provide the best performance in both tests, providing a factor of 10 to 50 improvement over GNILC at $k < 0.1\,{\rm hMpc^{-1}}$ in the higher redshift case and $k < 0.03 \,{\rm hMpc^{-1}}$ in the lower redshift case. However, none of the methods were found to recover the power spectrum satisfactorily at all BAO scales., Comment: 21 pages, 18 figures, accepted for publication by MRNAS

Intensity mapping experiments are beginning to measure the spatial distribution of neutral atomic hydrogen (HI) to constrain cosmological parameters and the large-scale distribution of matter. However, models of the behaviour of HI as a tracer of matter are complicated by galaxy evolution. In this work, we examine the clustering of HI in relation to galaxy colour, stellar mass, and HI mass in IllustrisTNG at $z$ = 0, 0.5, and 1. We compare the HI-red and HI-blue galaxy cross-power spectra, finding that HI-red has an amplitude 1.5 times higher than HI-blue at large scales. The cross-power spectra intersect at $\approx 3$ Mpc in real space and $\approx 10$ Mpc in redshift space, consistent with $z \approx 0$ observations. We show that HI clustering increases with galaxy HI mass and depends weakly on detection limits in the range $M_{\mathrm{HI}} \leq 10^8 M_\odot$. In terms of $M_\star$, we find blue galaxies in the greatest stellar mass bin cluster more than blue galaxies in other stellar mass bins. Red galaxies in the greatest stellar mass bin, however, cluster the weakest amongst red galaxies. These trends arise due to central-satellite compositions. Centrals correlate less with HI for increasing stellar mass, whereas satellites correlate more, irrespective of colour. Despite the clustering relationships with stellar mass, we find that the cross-power spectra are largely insensitive to detection limits in HI and galaxy surveys. Counter-intuitively, all auto and cross-power spectra for red and blue galaxies and HI decrease with time at all scales in IllustrisTNG. We demonstrate that processes associated with quenching contribute to this trend. The complex interplay between HI and galaxies underscores the importance of understanding baryonic effects when interpreting the large-scale clustering of HI, blue, and red galaxies at $z \leq 1$., Comment: 16 pages, 11 figures

This work presents the initial results of the drift-scan observation for the neutral hydrogen (HI) intensity mapping survey with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The data analyzed in this work were collected in night observations from 2019 through 2021. The primary findings are based on 28 hours of drift-scan observation carried out over seven nights in 2021, which covers $60\,{\rm deg}^2$ sky area. Our main findings are: (i) Our calibration strategy can successfully correct both the temporal and bandpass gain variation over the $4$-hour drift-scan observation. (ii) The continuum maps of the surveyed region are made with frequency resolution of $28$ kHz and pixel area of $2.95\,{\rm arcmin}^2$. The pixel noise levels of the continuum maps are slightly higher than the forecast assuming $T_{\rm sys}=20\,{\rm K}$, which are $36.0$ mK (for 10.0 s integration time) at the $1050$--$1150$ MHz band, and $25.9$ mK (for 16.7 s integration time) at the $1323$--$1450$ MHz band, respectively. (iii) The flux-weighted differential number count is consistent with the NRAO-VLA Sky Survey (NVSS) catalog down to the confusion limit $\sim7\,{\rm mJy}/{\rm beam}^{-1}$. (iv) The continuum flux measurements of the sources are consistent with that found in the literature. The difference in the flux measurement of $81$ isolated NVSS sources is about $6.3\%$. Our research offers a systematic analysis for the FAST HI intensity mapping drift-scan survey and serves as a helpful resource for further cosmology and associated galaxies sciences with the FAST drift-scan survey., Comment: 26 pages, 26 figures, and 4 tables

We present the science case for surveys with the Tianlai dish array interferometer tuned to the $\left[ 1300, 1400 \right] \mathrm{MHz}$ frequency range. Starting from a realistic generation of mock visibility data according to the survey strategy, we reconstruct a map of the sky and perform a foreground subtraction. We show that a survey of the North Celestial Polar cap during a year of observing time and covering an area of $150 \, \mathrm{deg^2}$ would reach a sensitivity of $ 1.5-2 \, \mathrm{mK} $ per $1 \, \mathrm{MHz} \times 0.25^2 \, \mathrm{deg^2 }$ voxel and be marginally impacted by mode-mixing. Tianlai would be able to detect a handful $(\sim 10)$ of nearby massive \HI clumps as well as a very strong cross-correlation signal of 21\,cm intensity maps with the North Celestial Cap Survey optical galaxies. We have also studied the performance of a mid-latitude survey, covering $\sim 1500 \, \mathrm{deg^2}$ centered on a declination of $\delta=55^\circ$, which overlaps the Sloan Digital Sky Survey footprint. Despite a higher noise level for the mid-latitude survey, as well as significant distortions due to mode mixing, Tianlai would be able to detect a highly significant cross-correlation between the 21\,cm signal and the Sloan spectroscopic galaxy sample. Using the extragalactic signals from either or both of these surveys, it will be possible to assess the impact of calibration uncertainties, antenna pattern uncertainties, sources of noise, and mode mixing for future surveys requiring higher sensitivity., Comment: 20 pages, 22 figures. Submitted to MNRAS

Microwave Kinetic Inductance Detectors (MKIDs) are highly scalable detectors that have demonstrated nearly background-limited sensitivity in the far-infrared from high-altitude balloon-borne telescopes and space-like laboratory environments. In addition, the detectors have a rich design space with many optimizable parameters, allowing highly sensitive measurements over a wide dynamic range. For these reasons, MKIDs were chosen for the Experiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM), a balloon-borne telescope targeting nearly background-limited performance in a high-altitude atmospheric environment from 420-540 GHz. We describe MKID optimization in the specific context of EXCLAIM and provide general results that apply to broader applications. Extending the established approach of tone frequency tracking, we show that readout power optimization enables significant, further improvement in dynamic range.

A great challenge for 21 cm intensity mapping experiments is the strong foreground radiation which is orders of magnitude brighter than the 21cm signal. Removal of the foreground takes advantage of the fact that its frequency spectrum is smooth while the redshifted 21cm signal spectrum is stochastic. However, a complication is the non-smoothness of the instrument response. This paper describes the electromagnetic simulation of the Tianlai cylinder array, a pathfinder for 21 cm intensity mapping experiments. Due to the vast scales involved, a direct simulation requires large amount of computing resources. We have made the simulation practical by using a combination of methods: first simulate a single feed, then an array of feed units, finally with the feed array and a cylindrical reflector together, to obtain the response for a single cylinder. We studied its radiation pattern, bandpass response and the effects of mutual coupling between feed units, and compared the results with observation. Many features seen in the measurement result are well reproduced in the simulation, especially the oscillatory features which are associated with the standing waves on the reflector. The mutual coupling between feed units is quantified with S-parameters, which decrease as the distance between the two feeds increases. Based on the simulated S-parameters, we estimate the correlated noise which has been seen in the visibility data, the results show very good agreement with the data in both magnitude and frequency structures. These results provide useful insights on the problem of 21cm signal extraction for real instruments., Comment: 24 pages, 24 figures

This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan., Comment: Contribution to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1908.01062, arXiv:1907.04473

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array., Comment: 24 pages, 6 figures. Updated to match published version

Hydrogen intensity mapping is a new field in astronomy that promises to make three-dimensional maps of the matter distribution of the Universe using the redshifted $21\,\textrm{cm}$ line of neutral hydrogen gas (HI). Several ongoing and upcoming radio interferometers, such as Tianlai, CHIME, HERA, HIRAX, etc. are using this technique. These instruments are designed to map large swaths of the sky by drift scanning over periods of many months. One of the challenges of the observations is that the daytime data is contaminated by strong radio signals from the Sun. In the case of Tianlai, this results in almost half of the measured data being unusable. We try to address this issue by developing an algorithm for solar contamination removal (AlgoSCR) from the radio data. The algorithm is based on an eigenvalue analysis of the visibility matrix, and hence is applicable only to interferometers. We apply AlgoSCR to simulated visibilities, as well as real daytime data from the Tianlai dish array. The algorithm can remove most of the solar contamination without seriously affecting other sky signals and thus makes the data usable for certain applications.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a balloon-borne far-infrared telescope that will survey star formation history over cosmological time scales to improve our understanding of why the star formation rate declined at redshift z < 2, despite continued clustering of dark matter. Specifically,EXCLAIM will map the emission of redshifted carbon monoxide and singly-ionized carbon lines in windows over a redshift range 0 < z < 3.5, following an innovative approach known as intensity mapping. Intensity mapping measures the statistics of brightness fluctuations of cumulative line emissions instead of detecting individual galaxies, thus enabling a blind, complete census of the emitting gas. To detect this emission unambiguously, EXCLAIM will cross-correlate with a spectroscopic galaxy catalog. The EXCLAIM mission uses a cryogenic design to cool the telescope optics to approximately 1.7 K. The telescope features a 90-cm primary mirror to probe spatial scales on the sky from the linear regime up to shot noise-dominated scales. The telescope optical elements couple to six {\mu}-Spec spectrometer modules, operating over a 420-540 GHz frequency band with a spectral resolution of 512 and featuring microwave kinetic inductance detectors. A Radio Frequency System-on-Chip (RFSoC) reads out the detectors in the baseline design. The cryogenic telescope and the sensitive detectors allow EXCLAIM to reach high sensitivity in spectral windows of low emission in the upper atmosphere. Here, an overview of the mission design and development status since the start of the EXCLAIM project in early 2019 is presented., Comment: SPIE Astronomical Telescopes + Instrumentation. arXiv admin note: substantial text overlap with arXiv:1912.07118

This work describes the design and implementation of optics for EXCLAIM, the EXperiment for Cryogenic Large-Aperture Intensity Mapping. EXCLAIM is a balloon-borne telescope that will measure integrated line emission from carbon monoxide (CO) at redshifts $z<1$ and ionized carbon ([CII]) at redshifts $z = 2.5-3.5$ to probe star formation over cosmic time in cross-correlation with galaxy redshift surveys. The EXCLAIM instrument is designed to observe at frequencies of $420$--$540$ GHz using six microfabricated silicon integrated spectrometers with spectral resolving power $R = 512$ coupled to kinetic inductance detectors (KIDs). A completely cryogenic telescope cooled to a temperature below 5~K provides low-background observations between narrow atmospheric lines in the stratosphere. Off-axis reflective optics use a $90$-cm primary mirror to provide $4.2^\prime$ full-width at half-maximum (FWHM) resolution at the center of the EXCLAIM band over a field of view of $22.5^\prime$. Illumination of the $1.7$ K cold stop combined with blackened baffling at multiple places in the optical system ensure low ($< -40$ dB) edge illumination of the primary to minimize spill onto warmer elements at the top of the dewar., Comment: 15 pages, 8 figures, 3 tables

The Tianlai Dish Pathfinder Array is a radio interferometer designed to test techniques for 21~cm intensity mapping in the post-reionization universe as a means for measuring large-scale cosmic structure. It performs drift scans of the sky at constant declination. We describe the design, calibration, noise level, and stability of this instrument based on the analysis of about $\sim 5 \%$ of 6,200 hours of on-sky observations through October, 2019. Beam pattern determinations using drones and the transit of bright sources are in good agreement, and compatible with electromagnetic simulations. Combining all the baselines, we make maps around bright sources and show that the array behaves as expected. A few hundred hours of observations at different declinations have been used to study the array geometry and pointing imperfections, as well as the instrument noise behaviour. We show that the system temperature is below 80~K for most feed antennas, and that noise fluctuations decrease as expected with integration time, at least up to a few hundred seconds. Analysis of long integrations, from 10 nights of observations of the North Celestial Pole, yielded visibilities with amplitudes of 20-30~mK, consistent with the expected signal from the NCP radio sky with $<10\,$mK precision for $1 ~\mathrm{MHz} \times 1~ \mathrm{min}$ binning. Hi-pass filtering the spectra to remove smooth spectrum signal yields a residual consistent with zero signal at the $0.5\,$mK level., Comment: 30 pages, 38 figures

CMB-S4---the next-generation ground-based cosmic microwave background (CMB) experiment---is set to significantly advance the sensitivity of CMB measurements and enhance our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. Among the science cases pursued with CMB-S4, the quest for detecting primordial gravitational waves is a central driver of the experimental design. This work details the development of a forecasting framework that includes a power-spectrum-based semi-analytic projection tool, targeted explicitly towards optimizing constraints on the tensor-to-scalar ratio, $r$, in the presence of Galactic foregrounds and gravitational lensing of the CMB. This framework is unique in its direct use of information from the achieved performance of current Stage 2--3 CMB experiments to robustly forecast the science reach of upcoming CMB-polarization endeavors. The methodology allows for rapid iteration over experimental configurations and offers a flexible way to optimize the design of future experiments given a desired scientific goal. To form a closed-loop process, we couple this semi-analytic tool with map-based validation studies, which allow for the injection of additional complexity and verification of our forecasts with several independent analysis methods. We document multiple rounds of forecasts for CMB-S4 using this process and the resulting establishment of the current reference design of the primordial gravitational-wave component of the Stage-4 experiment, optimized to achieve our science goals of detecting primordial gravitational waves for $r > 0.003$ at greater than $5\sigma$, or, in the absence of a detection, of reaching an upper limit of $r < 0.001$ at $95\%$ CL., Comment: 24 pages, 8 figures, 9 tables, submitted to ApJ. arXiv admin note: text overlap with arXiv:1907.04473

In 21~cm intensity mapping, the spectral smoothness of the foreground is exploited to separate it from the much weaker 21~cm signal. However, the non-smooth frequency response of the instrument complicates this process. Reflections and standing waves generate modulations on the frequency response. Here we report the analysis of the standing waves in the bandpass of the signal channels of the Tianlai Cylinder Array. By Fourier transforming the bandpass into the delay time domain, we find various standing waves generated on the telescope. A standing wave with time delay at about 142 ns is most clearly identified which is produced in the 15 meter feed cable. We also find a strong peak at a shorter delay of $\tau < 50 \ns$, which may be a mix of the standing wave between the reflector and feed, and the standing wave on the 4 m intermediate frequency (IF) cable. We also show that a smoother frequency response could be partially recovered by removing the reflection-inducted modulations. However, the standing wave on the antenna is direction-dependent, which poses a more difficult challenge for high precision calibration., Comment: 15 pages, 12 figures

Precision measurement of the beam pattern of an antenna is very important for many applications. While traditionally such measurement is often made in a microwave anechoic chamber or at a test range, measurement using an unmanned aerial vehicle offers a number of advantages: the measurement can be made for the assembled antenna on site, thus reflecting the actual characteristics of the antenna of interest, and more importantly, it can be performed for larger antennas which cannot be steered or easily measured using the anechoic chamber and test range. Here we report our beam measurement experiment with UAV for a 6 meter dish used in the Tianlai array, which is a radio astronomy experiment. Due to the dish's small collecting area, calibration with an astronomical source only allows for determining the antenna beam pattern over a very limited angular range. We describe in detail the setup of the experiment, the components of the signal transmitting system, the design of the flight path and the procedure for data processing. We find the UAV measurement of the beam pattern agrees very well with the astronomical source measurement in the main lobe, but the UAV measurement can be extended to the fourth side lobe. The measured position and width of each lobe also shows good agreement with electromagnetic field simulation. This UAV-based approach of beam pattern measurement is flexible and inexpensive, and the technique may also be applied to other experiments.

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array., Comment: 28 pages, 30 figures

The Packed Ultra-wideband Mapping Array (PUMA) is a proposed low-resolution transit interferometric radio telescope operating over the frequency range 200 - 1100MHz. Its rich science portfolio will include measuring structure in the universe from redshift z = 0.3 to 6 using 21cm intensity mapping, detecting one million fast radio bursts, and monitoring thousands of pulsars. It will allow PUMA to advance science in three different areas of physics (the physics of dark energy, the physics of cosmic inflation and time-domain astrophysics). This document is a response to a request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Decadal Survey on Astronomy and Astrophysics 2020. We present the science case of PUMA, the development path and major risks to the project., Comment: 46 pages, 16 figures, 7 tables; response to the request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Astro2020 Decadal Survey regarding PUMA APC submission (arXiv:1907.12559); v2: updated with correct bbl file

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey., Comment: Project White Paper submitted to the 2020 Decadal Survey, 10 pages plus references. arXiv admin note: substantial text overlap with arXiv:1907.04473

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift., Comment: 10 pages + references, 2 figures, 1 table; APC white paper submitted to the Astro2020 decadal survey. arXiv admin note: text overlap with arXiv:1810.09572

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively., Comment: 10 pages + references, 3 figures, 3 tables; project white paper submitted to the Astro2020 decadal survey; further details in updated arXiv:1810.09572

We present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4., Comment: 287 pages, 82 figures

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later., Comment: 5 pages + references; 1 figure; science white paper submitted to the Astro2020 decadal survey

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale., Comment: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

In this paper we describe QUBIC, an experiment that will observe the polarized microwave sky with a novel approach, which combines the sensitivity of state-of-the art bolometric detectors with the systematic effects control typical of interferometers. QUBIC unique features are the so-called "self-calibration", a technique that allows us to clean the measured data from instrumental effects, and its spectral imaging power, i.e. the ability to separate the signal in various sub-bands within each frequency band. QUBIC will observe the sky in two main frequency bands: 150 GHz and 220 GHz. A technological demonstrator is currently under testing and will be deployed in Argentina during 2019, while the final instrument is expected to be installed during 2020., Comment: Proceedings of the 2018 ICNFP conference, Crete. Published by Universe arXiv admin note: text overlap with arXiv:1801.03730

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics., Comment: This white-paper was prepared under auspices of the Department of Energy Cosmic Visions Dark Energy program; v1: initial release in Oct 2018; v2: Jul 2019 update: strawman design matches PUMA Decadal Submission, updated forecasting, new figures and science section; v3: fixed reference, missing author in metadata

The Tianlai Pathfinder is designed to demonstrate the feasibility of using a wide field of view radio interferometers to map the density of neutral hydrogen in the Universe after the Epoch of Reionizaton. This approach, called 21~cm intensity-mapping, promises an inexpensive means for surveying the large-scale structure of the cosmos. The Tianlai Pathfinder presently consists of an array of three, 15~m $\times$ 40~m cylinder telescopes and an array of sixteen, 6~m diameter dish antennas located in a radio-quiet part of western China. The two types of arrays were chosen to determine the advantages and disadvantages of each approach. The primary goal of the Pathfinder is to make 3D maps by surveying neutral hydrogen over large areas of the sky %$20,000 {\rm deg}^2$ in two different redshift ranges: first at $1.03 > z > 0.78$ ($700 - 800$~MHz) and later at $0.21 > z > 0.12$ ($1170 - 1270$~MHz). The most significant challenge to $21$~cm intensity-mapping is the removal of strong foreground radiation that dwarfs the cosmological signal. It requires exquisite knowledge of the instrumental response, i.e. calibration. In this paper, we provide an overview of the status of the Pathfinder and discuss the details of some of the analysis that we have carried out to measure the beam function of both arrays. We compare electromagnetic simulations of the arrays to measurements, discuss measurements of the gain and phase stability of the instrument, and provide a brief overview of the data processing pipeline., Comment: 18 pages, 10 figures

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., Comment: 191 pages

A wide bandwidth, dual polarized, modified four-square antenna is presented as a feed antenna for radio astronomical measurements. A linear array of these antennas is used as a line-feed for cylindrical reflectors for Tianlai, a radio interferometer designed for 21~cm intensity mapping. Simulations of the feed antenna beam patterns and scattering parameters are compared to experimental results at multiple frequencies across the 650 - 1420 MHz range. Simulations of the beam patterns of the combined feed array/reflector are presented as well., Comment: 15 pages, 15 figures

Fast Radio Bursts are bright, unresolved, non-repeating, broadband, millisecond flashes, found primarily at high Galactic latitudes, with dispersion measures much larger than expected for a Galactic source. The inferred all-sky burst rate is comparable to the core-collapse supernova rate out to redshift 0.5. If the observed dispersion measures are assumed to be dominated by the intergalactic medium, the sources are at cosmological distances with redshifts of 0.2 to 1. These parameters are consistent with a wide range of source models. One fast radio burst showed circular polarization [21(7)%] of the radio emission, but no linear polarization was detected, and hence no Faraday rotation measure could be determined. Here we report the examination of archival data revealing Faraday rotation in a newly detected burst - FRB 110523. It has radio flux at least 0.6 Jy and dispersion measure 623.30(5) pc cm$^{-3}$. Using Galactic contribution 45 pc cm$^{-3}$ and a model of intergalactic electron density, we place the source at a maximum redshift of 0.5. The burst has rotation measure -186.1(1.4) rad m$^{-2}$, much higher than expected for this line of sight through the Milky Way and the intergalactic medium, indicating magnetization in the vicinity of the source itself or within a host galaxy. The pulse was scattered by two distinct plasma screens during propagation, which requires either a dense nebula associated with the source or a location within the central region of its host galaxy. Keeping in mind that there may be more than one type of fast radio burst source, the detection in this instance of source-local magnetization and scattering favours models involving young stellar populations such as magnetars over models involving the mergers of older neutron stars, which are more likely to be located in low density regions of the host galaxy., Comment: 9 pages total, 6 figures, 1 table. Published in Nature

We present in this paper a new Bayesian semi-blind approach for foreground removal in observations of the 21-cm signal with interferometers. The technique, which we call HIEMICA (HI Expectation-Maximization Independent Component Analysis), is an extension of the Independent Component Analysis (ICA) technique developed for two-dimensional (2D) CMB maps to three-dimensional (3D) 21-cm cosmological signals measured by interferometers. This technique provides a fully Bayesian inference of power spectra and maps and separates the foregrounds from signal based on the diversity of their power spectra. Only relying on the statistical independence of the components, this approach can jointly estimate the 3D power spectrum of the 21-cm signal and, the 2D angular power spectrum and the frequency dependence of each foreground component, without any prior assumptions about foregrounds. This approach has been tested extensively by applying it to mock data from interferometric 21-cm intensity mapping observations under idealized assumptions of instrumental effects. We also discuss the impact when the noise properties are not known completely. As a first step toward solving the 21 cm power spectrum analysis problem we compare the semi-blind HIEMICA technique with the commonly used Principal Component Analysis (PCA). Under the same idealized circumstances the proposed technique provides significantly improved recovery of the power spectrum. This technique can be applied straightforwardly to all 21-cm interferometric observations, including epoch of reionization measurements, and can be extended to single-dish observations as well., Comment: 18 pages, 7 figures, added some discussions about the impact of noise misspecification

We present results from an end-to-end simulation pipeline interferometric observations of cosmic microwave background polarization. We use both maximum-likelihood and Gibbs sampling techniques to estimate the power spectrum. In addition, we use Gibbs sampling for image reconstruction from interferometric visibilities. The results indicate the level to which various systematic errors (e.g., pointing errors, gain errors, beam shape errors, cross- polarization) must be controlled in order to successfully detect and characterize primordial B modes as well as other scientific goals. In addition, we show that Gibbs sampling is an effective method of image reconstruction for interferometric data in other astrophysical contexts., Comment: To appear in the proceedings of IAU Symposium 306: Statistical Challenges in 21st-Century Cosmology

We present a novel, general-purpose method for deconvolving and denoising images from gridded radio interferometric visibilities using Bayesian inference based on a Gaussian process model. The method automatically takes into account incomplete coverage of the uv-plane, signal mode coupling due to the primary beam, and noise mode coupling due to uv sampling. Our method uses Gibbs sampling to efficiently explore the full posterior distribution of the underlying signal image given the data. We use a set of widely diverse mock images with a realistic interferometer setup and level of noise to assess the method. Compared to results from a proxy for point source- based CLEAN method we find that in terms of RMS error and signal-to-noise ratio our approach performs better than traditional deconvolution techniques, regardless of the structure of the source image in our test suite. Our implementation scales as O(np log np), provides full statistical and uncertainty information of the reconstructed image, requires no supervision, and provides a robust, consistent framework for incorporating noise and parameter marginalizations and foreground removal., Comment: 12 pages, 10 figures, MNRAS accepted. arXiv admin note: text overlap with arXiv:1109.4640

The detection of the primordial $B$-mode spectrum of the polarized cosmic microwave background (CMB) signal may provide a probe of inflation. However, observation of such a faint signal requires excellent control of systematic errors. Interferometry proves to be a promising approach for overcoming such a challenge. In this paper we present a complete simulation pipeline of interferometric observations of CMB polarization, including systematic errors. We employ two different methods for obtaining the power spectra from mock data produced by simulated observations: the maximum likelihood method and the method of Gibbs sampling. We show that the results from both methods are consistent with each other, as well as, within a factor of 6, with analytical estimates. Several categories of systematic errors are considered: instrumental errors, consisting of antenna gain and antenna coupling errors, and beam errors, consisting of antenna pointing errors, beam cross-polarization and beam shape (and size) errors. In order to recover the tensor-to-scalar ratio, $r$, within a 10% tolerance level, which ensures the experiment is sensitive enough to detect the $B$-signal at $r=0.01$ in the multipole range $28 < \ell < 384$, we find that, for a QUBIC-like experiment, Gaussian-distributed systematic errors must be controlled with precisions of $|g_{rms}| = 0.1$ for antenna gain, $|\epsilon_{rms}| = 5 \times 10^{-4}$ for antenna coupling, $\delta_{rms} \approx 0.7^\circ$ for pointing, $\zeta_{rms} \approx 0.7^\circ$ for beam shape, and $\mu_{rms} = 5 \times 10^{-4}$ for beam cross-polarization., Comment: 15 pages, 6 figures, submitted to ApJS

Detection of B-mode polarization of the cosmic microwave background (CMB) radiation is one of the frontiers of observational cosmology. Because they are an order of magnitude fainter than E-modes, it is quite a challenge to detect B-modes. Having more manageable systematics, interferometers prove to have a substantial advantage over imagers in detecting such faint signals. Here, we present a method for Bayesian inference of power spectra and signal reconstruction from interferometric data of the CMB polarization signal by using the technique of Gibbs sampling. We demonstrate the validity of the method in the flat-sky approximation for a simulation of an interferometric observation on a finite patch with incomplete uv-plane coverage, a finite beam size and a realistic noise model. With a computational complexity of O(n^{3/2}), n being the data size, Gibbs sampling provides an efficient method for analyzing upcoming cosmology observations., Comment: 8 pages, 8 figures, expanded discussion and edited to match ApJS approved version, acknowledgments updated

We investigate the impact of instrumental systematic errors in interferometric measurements of the cosmic microwave background (CMB) temperature and polarization power spectra. We simulate interferometric CMB observations to generate mock visibilities and estimate power spectra using the statistically optimal maximum likelihood technique. We define a quadratic error measure to determine allowable levels of systematic error that do not induce power spectrum errors beyond a given tolerance. As an example, in this study we focus on differential pointing errors. The effects of other systematics can be simulated by this pipeline in a straightforward manner. We find that, in order to accurately recover the underlying B-modes for r=0.01 at 28

Measurements of the 2.7 K cosmic microwave background (CMB) radiation now provide the most stringent constraints on cosmological models. The power spectra of the temperature anisotropies and the $E$-mode polarization of the CMB are explained well by the inflationary paradigm. The next generation of CMB experiments aim at providing the most direct evidence for inflation through the detection of $B$-modes in the CMB polarization, presumed to have been caused by gravitational waves generated during the inflationary epoch around $10^{-34}$s. The $B$-mode polarization signals are very small ($\leq$10$^{-8}$K) compared with the temperature anisotropies ($\sim 10^{-4}$K). Systematic effects in CMB telescopes can cause leakage from temperature anisotropy into polarization. Bolometric interferometry (BI) is a novel approach to measuring this small signal with lower leakage. If BI can be made to work over wide bandwidth ($\sim20-30\%$) it can provide similar sensitivity to imagers. Subdividing the frequency passband of a Fizeau interferometer would mitigate the problem of `fringe smearing.' Furthermore, the approach should allow simultaneous measurements in image space and visibility space. For subdividing the frequency passsband (`sub-band splitting' henceforth), we write an expression for the output from every baseline at every detector in the focal plane as a sum of visibilities in different frequency sub-bands. For operating the interferometer simultaneously as an imager, we write the output as two integrals over the sky and the focal plane, with all the phase differences accounted for.}{The sub-band splitting method described here is general and can be applied to broad-band Fizeau interferometers across the electromagnetic spectrum. Applications to CMB measurements and to long-baseline optical interferometry are promising., Comment: 8 pages, 5 figures, submitted to Astronomy and Astrophysics

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages., Comment: Science White Paper submitted to the US Astro2010 Decadal Survey. Full list of 177 author available at http://cmbpol.uchicago.edu

Using the 21 cm line, observed all-sky and across the redshift range from 0 to 5, the large scale structure of the Universe can be mapped in three dimensions. This can be accomplished by studying specific intensity with resolution ~ 10 Mpc, rather than via the usual galaxy redshift survey. The data set can be analyzed to determine Baryon Acoustic Oscillation wavelengths, in order to address the question: 'What is the nature of Dark Energy?' In addition, the study of Large Scale Structure across this range addresses the questions: 'How does Gravity effect very large objects?' and 'What is the composition our Universe?' The same data set can be used to search for and catalog time variable and transient radio sources., Comment: White Paper for the Astro2010 Astronomy Decadal Review

The White Mountain Polarimeter (WMPol) is a dedicated ground-based microwave telescope and receiver system for observing polarization of the Cosmic Microwave Background. WMPol is located at an altitude of 3880 meters on a plateau in the White Mountains of Eastern California, USA, at the Barcroft Facility of the University of California White Mountain Research Station. Presented here is a description of the instrument and the data collected during April through October 2004. We set an upper limit on $E$-mode polarization of 14 $\mu\mathrm{K}$ (95% confidence limit) in the multipole range $170<\ell<240$. This result was obtained with 422 hours of observations of a 3 $\mathrm{deg}^2$ sky area about the North Celestial Pole, using a 42 GHz polarimeter. This upper limit is consistent with $EE$ polarization predicted from a standard $\Lambda$-CDM concordance model., Comment: 35 pages. 12 figures. To appear in ApJS

One of the most spectacular scientific breakthroughs in past decades was using measurements of the fluctuations in the cosmic microwave background (CMB) to test precisely our understanding of the history and composition of the Universe. This report presents a roadmap for leading CMB research to its logical next step, using precision polarization measurements to learn about ultra-high-energy physics and the Big Bang itself., Comment: This is the final report of the DoE/NASA/NSF interagency task force on CMB research chaired by Rai Weiss. 87 pages

COMPASS is an on-axis 2.6 meter telescope coupled to a correlation polarimeter operating at a wavelength of 1 cm. The entire instrument was built specifically for CMB polarization studies. We report here on observations of February 2001 - April 2001 using this system. We set an upper limit on E-mode polarized anisotropies of 33.5 uK (95% confidence limit) in the l-range 200-600., Comment: Submitted to ApJ 8/19/2003

COMPASS is an on-axis 2.6 meter telescope coupled to a correlation polarimeter. The entire instrument was built specifically for CMB polarization studies. Careful attention was given to receiver and optics design, stability of the pointing platform, avoidance of systematic offsets, and development of data analysis techniques. Here we describe the experiment, its strengths and weaknesses, and the various things we have learned that may benefit future efforts to measure the polarization of the CMB., Comment: To be published in the proceedings of "The Cosmic Microwave Background and its Polarization", New Astronomy Reviews, (eds. S. Hanany and K.A. Olive)

Most of the useful information about inflationary gravitational waves and reionization is on large angular scales where Galactic foreground contamination is the worst, so a key challenge is to model, quantify and remove polarized foregrounds. We use the Leiden radio surveys to quantify the polarized synchrotron radiation at large angular scales, which is likely to be the most challenging polarized contaminant for the WMAP satellite. We find that the synchrotron E- and B-contributions are equal to within 10% from 408-820MHz with a hint of E-domination at higher frequencies. We quantify Faraday Rotation & Depolarization effects and show that they cause the synchrotron polarization percentage to drop both towards lower frequencies and towards lower multipoles., Comment: To be published in the proceedings of "The Cosmic Microwave Background and its Polarization", New Astronomy Reviews, (eds. S. Hanany and K.A. Olive). 4 pages. More foreground information at http://www.hep.upenn.edu/~angelica/foreground.html#polar or from angelica@higgs.hep.upenn.edu