Magnetars are highly magnetized young neutron stars that occasionally produce enormous bursts and flares of X-rays and gamma-rays. Of the approximately thirty magnetars currently known in our Galaxy and Magellanic Clouds, five have exhibited transient radio pulsations. Fast radio bursts (FRBs) are millisecond-duration bursts of radio waves arriving from cosmological distances. Some have been seen to repeat. A leading model for repeating FRBs is that they are extragalactic magnetars, powered by their intense magnetic fields. However, a challenge to this model has been that FRBs must have radio luminosities many orders of magnitude larger than those seen from known Galactic magnetars. Here we report the detection of an extremely intense radio burst from the Galactic magnetar SGR 1935+2154 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project. The fluence of this two-component bright radio burst and the estimated distance to SGR 1935+2154 together imply a 400-800 MHz burst energy of $\sim 3 \times 10^{34}$ erg, which is three orders of magnitude brighter than those of any radio-emitting magnetar detected thus far. Such a burst coming from a nearby galaxy would be indistinguishable from a typical FRB. This event thus bridges a large fraction of the radio energy gap between the population of Galactic magnetars and FRBs, strongly supporting the notion that magnetars are the origin of at least some FRBs.

We evaluate the ability of Convolutional Neural Networks (CNNs) to predict galaxy cluster masses in the BAHAMAS hydrodynamical simulations. We train four separate single-channel networks using: stellar mass, soft X-ray flux, bolometric X-ray flux, and the Compton $y$ parameter as observational tracers, respectively. Our training set consists of $\sim$6400 synthetic cluster images generated from the simulation, while an additional $\sim$1600 images form a test set. We also train a `multi-channel' CNN by combining the four observational tracers. We utilize $\texttt{Keras}$ with a $\texttt{Tensorflow}$ backend to train the network, and all four converge within 2000 epochs. The cluster masses predicted from these networks are evaluated using the average fractional difference between predicted cluster mass and true cluster mass. The resulting predictions are especially precise for halo masses in the range $10^{13.25}M_{\odot}<M<10^{14.5}M_{\odot}$, where all five networks produce mean mass biases of order $\approx 1\%$ with a scatter on the mean bias of $\approx 0.5\%$. The network trained with Compton $y$ parameter maps yields the most precise predictions. We interpret the network's behaviour using two diagnostic tests to determine which features are used to predict cluster mass. The CNN trained with stellar mass images detect galaxies (not surprisingly), while CNNs trained with gas-based tracers utilise the shape of the signal to estimate cluster mass.

We report on the discovery and analysis of bursts from nine new repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 195 to 1380 pc cm$^{-3}$. We detect two bursts from three of the new sources, three bursts from four of the new sources, four bursts from one new source, and five bursts from one new source. We determine sky coordinates of all sources with uncertainties of $\sim$10$^\prime$. We detect Faraday rotation measures for two sources, with values $-20(1)$ and $-499.8(7)$ rad m$^{-2}$, that are substantially lower than the RM derived from bursts emitted by FRB 121102. We find that the DM distribution of our events, combined with the nine other repeaters discovered by CHIME/FRB, is indistinguishable from that of thus far non-repeating CHIME/FRB events. However, as previously reported, the burst widths appear statistically significantly larger than the thus far non-repeating CHIME/FRB events, further supporting the notion of inherently different emission mechanisms and/or local environments. These results are consistent with previous work, though are now derived from 18 repeating sources discovered by CHIME/FRB during its first year of operation. We identify candidate galaxies that may contain FRB 190303.J1353+48 (DM = 222.4 pc cm$^{-3}$).

The location of a galaxy cluster's centroid is typically derived from observations of the galactic and/or gas component of the cluster, but these typically deviate from the true centre. This can produce bias when observations are combined to study average cluster properties. Using data from the BAHAMAS cosmological hydrodynamic simulations we study this bias in both two and three dimensions for 2000 clusters over the $10^{13} - 10^{15} ~\mathrm{M_{\odot}}$ mass range. We quantify and model the offset distributions between observationally-motivated centres and the `true' centre of the cluster, which is taken to be the most gravitationally bound particle measured in the simulation. We fit the cumulative distribution function of offsets with an exponential distribution and a Gamma distribution fit well with most of the centroid definitions. The galaxy-based centres can be seen to be divided into a mis-centred group and a well-centred group, with the well-centred group making up about $60\%$ of all the clusters. Gas-based centres are overall less scattered than galaxy-based centres. We also find a cluster-mass dependence of the offset distribution of gas-based centres, with generally larger offsets for smaller mass clusters. We then measure cluster density profiles centred at each choice of the centres and fit them with empirical models. Stacked, mis-centred density profiles fit to the Navarro-Frenk-White dark-matter profile and Komatsu-Seljak gas profile show that recovered shape and size parameters can significantly deviate from the true values. For the galaxy-based centres, this can lead to cluster masses being underestimated by up to $10\%$.

We construct the mean thermal Sunyaev-Zel'dovich (tSZ) Comptonization y profile around Luminous Red Galaxies (LRGs) in the redshift range 0.16 < z < 0.47 from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) using the Planck y map. The mean central tSZ signal for the full sample is y ~ 1.8 * 10^(-7) and we detect tSZ emission out to ~30 arcmin, which is well beyond the 10 arcmin angular resolution of the y map and well beyond the virial radii of the LRGs. We compare the measured profile with predictions from the cosmo-OWLS suite of cosmological hydrodynamical simulations. This comparison agrees well for models that include feedback from active galactic nuclei (AGN), but not with hydrodynamic models without this energetic feedback mechanism. This suggests that an additional heating mechanism is required over SNe feedback and star formation to explain the y data profile. We also compare our results with predictions based on the halo model with a universal pressure profile (UPP) giving the y signal. The predicted profile is consistent with the data, but only if we account for the clustering of haloes via a two-halo term and if halo masses are estimated using the mean stellar-to-halo mass (SHM) relation of Coupon et al. (2015) or Wang et al.(2016) estimated from gravitational lensing measurements. We also discuss the importance of scatter in the SHM relation on the model predictions.

We search the Planck data for a thermal Sunyaev-Zel'dovich (tSZ) signal due to gas filaments between pairs of Luminous Red Galaxies (LRG's) taken from the Sloan Digital Sky Survey Data Release 12 (SDSS/DR12). We identify $\sim$260,000 LRG pairs in the DR12 catalog that lie within 6-10 $h^{-1} \mathrm{Mpc}$ of each other in tangential direction and within 6 $h^{-1} \mathrm{Mpc}$ in radial direction. We stack pairs by rotating and scaling the angular positions of each LRG so they lie on a common reference frame, then we subtract a circularly symmetric halo from each member of the pair to search for a residual signal between the pair members. We find a statistically significant (5.3$\sigma$) signal between LRG pairs in the stacked data with a magnitude $\Delta y = (1.31 \pm 0.25) \times 10^{-8}$. The uncertainty is estimated from two Monte Carlo null tests which also establish the reliability of our analysis. Assuming a simple, isothermal, cylindrical filament model of electron over-density with a radial density profile proportional to $r_c/r$ (as determined from simulations), where $r$ is the perpendicular distance from the cylinder axis and $r_c$ is the core radius of the density profile, we constrain the product of over-density and filament temperature to be $\delta_c \times (T_{\rm e}/10^7 \, {\rm K}) \times (r_c/0.5h^{-1} \, {\rm Mpc}) = 2.7 \pm 0.5$. To our knowledge, this is the first detection of filamentary gas at over-densities typical of cosmological large-scale structure. We compare our result to the BAHAMAS suite of cosmological hydrodynamic simulations (McCarthy et al. 2017) and find a slightly lower, but marginally consistent Comptonization excess, $\Delta y = (0.84 \pm 0.24) \times 10^{-8}$.

The Canadian Hydrogen Observatory and Radio-transient Detector (CHORD) is a next-generation radio telescope, proposed for construction to start immediately. CHORD is a pan-Canadian project, designed to work with and build on the success of the Canadian Hydrogen Intensity Mapping Experiment (CHIME). It is an ultra-wideband, "large-N, small-D" telescope, consisting of a central array of 512x6-m dishes, supported by a pair of distant outrigger stations, each equipped with CHIME-like cylinders and a 64-dish array. CHORD will measure the distribution of matter over a huge swath of the Universe, detect and localize tens of thousands of Fast RadioBursts (FRBs), and undertake cutting-edge measurements of fundamental physics.

We report on the discovery of eight repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 103.5 to 1281 pc cm$^{-3}$. They display varying degrees of activity: six sources were detected twice, another three times, and one ten times. These eight repeating FRBs likely represent the bright and/or high-rate end of a distribution of infrequently repeating sources. For all sources, we determine sky coordinates with uncertainties of $\sim$10$^\prime$. FRB 180916.J0158+65 has a burst-averaged DM = $349.2 \pm 0.3$ pc cm$^{-3}$ and a low DM excess over the modelled Galactic maximum (as low as $\sim$20 pc cm$^{-3}$); this source also has a Faraday rotation measure (RM) of $-114.6 \pm 0.6$ rad m$^{-2}$, much lower than the RM measured for FRB 121102. FRB 181030.J1054+73 has the lowest DM for a repeater, $103.5 \pm 0.3$ pc cm$^{-3}$, with a DM excess of $\sim$ 70 pc cm$^{-3}$. Both sources are interesting targets for multi-wavelength follow-up due to their apparent proximity. The DM distribution of our repeater sample is statistically indistinguishable from that of the first 12 CHIME/FRB sources that have not repeated. We find, with 4$\sigma$ significance, that repeater bursts are generally wider than those of CHIME/FRB bursts that have not repeated, suggesting different emission mechanisms. Our repeater events show complex morphologies that are reminiscent of the first two discovered repeating FRBs. The repetitive behavior of these sources will enable interferometric localizations and subsequent host galaxy identifications.

Canada has thriving communities in CMB (cosmic microwave background) studies, cosmology and submillimetre (submm) astronomy, with involvement in many facilities that featured prominently in previous Astronomy Long Range Plans. The standard cosmological model continues to be well fit using a small number of parameters. No one expects this model to be complete and so we need to continue to challenge it with data; moreover, it does not explain how galaxies and other structures form. So, how do we improve the precision of our understanding of structure formation within this model? Wavelengths from the microwave to the submm will be particularly fruitful for answering this question. That's because, in addition to the CMB anisotropies, there are other signals that can be extracted from large maps at these wavelengths - particularly the cosmic infrared and submm backgrounds, the thermal and kinetic Sunyaev-Zeldovich effects, and CMB lensing. Such signals carry a wealth of information about the cosmological model, as well as how dust, gas and star-formation evolve within dark-matter halos. Cross-correlations between these signals and those coming from the radio, optical and X-ray surveys, will provide even more information. Canadians are already members of teams for several related facilities and are working to be involved in others. In order for Canada to be fully engaged in exploiting the detailed information coming from these cosmological signatures, it is crucial that we find the resources to participate competitively in a combination of projects currently being planned. Examples include CMB-S4, CCAT-prime, AtLAST, a new camera for JCMT, balloon projects such as BFORE and a future ambitious CMB satellite.

We search for the potential contamination in the Planck thermal Sunyaev-Zeldovich (tSZ) map by calculating the cross-correlation between the tSZ signal and weak lensing by large scale structure and the Cosmic Microwave Background (CMB). The lensing data we use is the convergence map from the Red Sequence Cluster Lensing Survey (RCSLens) and the Planck CMB lensing map. We reconstruct the tSZ $y$ map with a Needlet Internal Linear Combination method using the HFI sky maps from the Planck satellite. We remove the CMB signal while minimizing the residual noise. The cross correlation signal from our reconstructed $y$ map is consistent with that from the Planck team's NILC $y$ map. The CIB and galactic dust emission are two potential sources of contamination in the reconstructed $y$ map. We remove the CIB signal by subtracting the CIB maps reconstructed by Planck collaboration from the raw temperature maps. We find that cross-correlation between the CIB and galactic lensing contributes to $(7.5\pm6.0)\%$ in the Planck NILC tSZ cross galactic lensing signal within $100<\ell<2000$, which implies that previous detections of the tSZ cross galactic lensing is robust to CIB contaminations. In contrast, the Planck NILC tSZ cross CMB lensing is biased by $(18.4\pm2.8)\%$ in the same $\ell$ range. Galactic dust contamination is tested by projecting out a grey-body dust models with different dust spectral indices. Galactic dust does not affect galactic lensing cross tSZ signal significantly.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

Angular power spectra computed from Planck HFI 353 GHz intensity and polarization maps produce a TB correlation that can be approximated by a power-law (Planck Collaboration XI 2018). Whether the observed TB correlation is an induced systematic feature or a physical property of Galactic dust emission is of interest both for cosmological and Galactic studies. We investigate the large angular scale E and B-mode morphology of microwave polarized thermal dust emission, and relate it to physical quantities of polarization angle and polarization fraction. We use empirical models of polarized dust to show that dust polarization angle is a key factor in producing the TB correlation. A small sample of both simulated and observed polarization angle maps are combined with 353 GHz intensity and dust polarization fraction to produce a suite of maps from which we compute TB and EB. Model realizations that produce a positive TB correlation are common and can result from large-scale (>5 degree) structure in the polarization angle. The TB correlation appears robust to introduction of individual intensity, polarization angle and polarization fraction model components that are independent of the 353 GHz observations. We conclude that the observed TB correlation is likely the result of large-scale Galactic dust polarization properties.

We report the detection of a single burst from the first-discovered repeating Fast Radio Burst source, FRB 121102, with CHIME/FRB, which operates in the frequency band 400-800 MHz. The detected burst occurred on 2018 November 19 and its emission extends down to at least 600 MHz, the lowest frequency detection of this source yet. The burst, detected with a significance of 23.7$\sigma$, has fluence 12$\pm$3 Jy ms and shows complex time and frequency morphology. The 34 ms width of the burst is the largest seen for this object at any frequency. We find evidence of sub-burst structure that drifts downward in frequency at a rate of -3.9$\pm$0.2 MHz ms$^{-1}$. Our best fit tentatively suggests a dispersion measure of 563.6$\pm$0.5 pc cm$^{-3}$, which is ${\approx}$1% higher than previously measured values. We set an upper limit on the scattering time at 500 MHz of 9.6 ms, which is consistent with expectations from the extrapolation from higher frequency data. We have exposure to the position of FRB 121102 for a total of 11.3 hrs within the FWHM of the synthesized beams at 600 MHz from 2018 July 25 to 2019 February 25. We estimate on the basis of this single event an average burst rate for FRB 121102 of 0.1-10 per day in the 400-800 MHz band for a median fluence threshold of 7 Jy ms in the stated time interval.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

The discovery of a repeating Fast Radio Burst (FRB) source, FRB 121102, eliminated models involving cataclysmic events for this source. No other repeating FRB has yet been detected in spite of many recent FRB discoveries and follow-ups, suggesting repeaters may be rare in the FRB population. Here we report the detection of six repeat bursts from FRB 180814.J0422+73, one of the 13 FRBs detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project during its pre-commissioning phase in July and August 2018. These repeat bursts are consistent with originating from a single position on the sky, with the same dispersion measure (DM), ~189 pc cm-3. This DM is approximately twice the expected Milky Way column density, and implies an upper limit on the source redshift of 0.1, at least a factor of ~2 closer than FRB 121102. In some of the repeat bursts, we observe sub-pulse frequency structure, drifting, and spectral variation reminiscent of that seen in FRB 121102, suggesting similar emission mechanisms and/or propagation effects. This second repeater, found among the first few CHIME/FRB discoveries, suggests that there exists -- and that CHIME/FRB and other wide-field, sensitive radio telescopes will find -- a substantial population of repeating FRBs.

Fast radio bursts (FRBs) are highly dispersed millisecond-duration radio flashes likely arriving from far outside the Milky Way galaxy. This phenomenon was discovered at radio frequencies near 1.4 GHz and to date has been observed in one case at as high as 8 GHz, but not below 700 MHz in spite of significant searches at low frequencies. Here we report detections of FRBs at radio frequencies as low as 400 MHz, on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) using the CHIME/FRB instrument. We present 13 FRBs detected during a telescope pre-commissioning phase, when our sensitivity and field-of-view were not yet at design specifications. Emission in multiple events is seen down to 400 MHz, the lowest radio frequency to which we are sensitive. The FRBs show a variety of temporal scattering behaviours, with the majority significantly scattered, and some apparently unscattered to within measurement uncertainty even at our lowest frequencies. Of the 13 reported here, one event has the lowest dispersion measure yet reported, implying it is among the closest yet known, and another has shown multiple repeat bursts, as described in a companion paper. Our low-scattering events suggest that efforts to detect FRBs at radio frequencies below 400 MHz will eventually be successful. The overall scattering properties of our sample suggest that FRBs as a class are preferentially located in environments that scatter radio waves more strongly than the diffuse interstellar medium (ISM) in the Milky Way.

The reionization optical depth is the most poorly determined of the six $\Lambda$CDM parameters fit to CMB anisotropy data. Instrumental noise and systematics have prevented uncertainties from reaching their cosmic variance limit. At present, the datasets providing the most statistical constraining power are the WMAP, Planck LFI, and Planck HFI full-sky polarization maps. As the reprocessed HFI data with reduced systematics are not yet publicly unavailable, we examine determinations of $\tau$ using 9-year WMAP and 2015 Planck LFI data, with an emphasis on characterizing potential systematic bias resulting from foreground template and masking choices. We find evidence for a low-level systematic in the LFI polarization data with a roughly common-mode morphology across the LFI frequencies and a spectrum consistent with leakage of intensity signal into the polarization channels. We demonstrate significant bias in the optical depth derived when using the LFI 30 GHz map as a template to clean synchrotron from WMAP data, and recommend against use of the 2015 LFI 30 GHz polarization data as a foreground template for non-LFI datasets. We find an inconsistency between versions of the 2015 polarized 353 GHz dust templates reconstructed from the Planck likelihood and those from delivered maps, which can affect $\tau$ at the 1$\sigma$ level. The spread in $\tau$ values over the ensemble of data combinations we study suggests that systematic uncertainties still contribute significantly to the current uncertainty in $\tau$, but all values are consistent with the range of $\tau$ = 0.07 +/- 0.02.

We analyze simulated maps of the Cosmology Large Angular Scale Surveyor (CLASS) experiment and recover a nearly cosmic-variance limited estimate of the reionization optical depth $\tau$. We use a power spectrum-based likelihood to simultaneously clean foregrounds and estimate cosmological parameters in multipole space. Using software specifically designed to constrain $\tau$, the amplitude of scalar fluctuations $A_s$, and the tensor-to-scalar ratio $r$, we demonstrate that the CLASS experiment will be able to estimate $\tau$ within a factor of two of the full-sky cosmic variance limit allowed by cosmic microwave background polarization measurements. Additionally, we discuss the role of CLASS's $\tau$ constraint in conjunction with gravitational lensing of the CMB on obtaining a $\gtrsim4\sigma$ measurement of the sum of the neutrino masses.

The Cosmology Large Angular Scale Surveyor consists of four instruments performing a CMB polarization survey. Currently, the 40 GHz and first 90 GHz instruments are deployed and observing, with the second 90 GHz and a multichroic 150/220 GHz instrument to follow. The receiver is a central component of each instrument's design and functionality. This paper describes the CLASS receiver design, using the first 90 GHz receiver as a primary reference. Cryogenic cooling and filters maintain a cold, low-noise environment for the detectors. We have achieved receiver detector temperatures below 50 mK in the 40 GHz instrument for 85% of the initial 1.5 years of operation, and observed in-band efficiency that is consistent with pre-deployment estimates. At 90 GHz, less than 26% of in-band power is lost to the filters and lenses in the receiver, allowing for high optical efficiency. We discuss the mounting scheme for the filters and lenses, the alignment of the cold optics and detectors, stray light control, and magnetic shielding.

The search for inflationary primordial gravitational waves and the measurement of the optical depth to reionization, both through their imprint on the large angular scale correlations in the polarization of the cosmic microwave background (CMB), has created the need for high sensitivity measurements of polarization across large fractions of the sky at millimeter wavelengths. These measurements are subject to instrumental and atmospheric $1/f$ noise, which has motivated the development of polarization modulators to facilitate the rejection of these large systematic effects. Variable-delay polarization modulators (VPMs) are used in the Cosmology Large Angular Scale Surveyor (CLASS) telescopes as the first element in the optical chain to rapidly modulate the incoming polarization. VPMs consist of a linearly polarizing wire grid in front of a movable flat mirror. Varying the distance between the grid and the mirror produces a changing phase shift between polarization states parallel and perpendicular to the grid which modulates Stokes U (linear polarization at $45^\circ$) and Stokes V (circular polarization). The CLASS telescopes have VPMs as the first optical element from the sky; this simultaneously allows a lock-in style polarization measurement and the separation of sky polarization from any instrumental polarization further along in the optical path. The Q-band CLASS VPM was the first VPM to begin observing the CMB full time, starting in the Spring of 2016. The first W-band CLASS VPM was installed in the Spring of 2018.

The Cosmology Large Angular Scale Surveyor (CLASS) aims to detect and characterize the primordial B-mode signal and make a sample-variance-limited measurement of the optical depth to reionization. CLASS is a ground-based, multi-frequency microwave polarimeter that surveys 70% of the microwave sky every day from the Atacama Desert. The focal plane detector arrays of all CLASS telescopes contain smooth-walled feedhorns that couple to transition-edge sensor (TES) bolometers through symmetric planar orthomode transducer (OMT) antennas. These low noise polarization-sensitive detector arrays are fabricated on mono-crystalline silicon wafers to maintain TES uniformity and optimize optical efficiency throughout the wafer. In this paper, we discuss the design and characterization of the first CLASS 93 GHz detector array. We measure the dark parameters, bandpass, and noise spectra of the detectors and report that the detectors are photon-noise limited. With current array yield of 82%, we estimate the total array noise-equivalent power (NEP) to be 2.1 aW$\sqrt[]{\mathrm{s}}$.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a novel transit radio telescope operating across the 400-800-MHz band. CHIME is comprised of four 20-m x 100-m semi-cylindrical paraboloid reflectors, each of which has 256 dual-polarization feeds suspended along its axis, giving it a >200 square degree field-of-view. This, combined with wide bandwidth, high sensitivity, and a powerful correlator makes CHIME an excellent instrument for the detection of Fast Radio Bursts (FRBs). The CHIME Fast Radio Burst Project (CHIME/FRB) will search beam-formed, high time-and frequency-resolution data in real time for FRBs in the CHIME field-of-view. Here we describe the CHIME/FRB backend, including the real-time FRB search and detection software pipeline as well as the planned offline analyses. We estimate a CHIME/FRB detection rate of 2-42 FRBs/sky/day normalizing to the rate estimated at 1.4-GHz by Vander Wiel et al. (2016). Likely science outcomes of CHIME/FRB are also discussed. CHIME/FRB is currently operational in a commissioning phase, with science operations expected to commence in the latter half of 2018.

We examine the impact of baryon acoustic oscillation (BAO) scale measurements on the discrepancy between the value of the Hubble constant ($H_0$) inferred from the local distance ladder and from Planck cosmic microwave background (CMB) data. While the BAO data alone cannot constrain $H_0$, we show that combining the latest BAO results with WMAP, Atacama Cosmology Telescope (ACT), or South Pole Telescope (SPT) CMB data produces values of $H_0$ that are $2.4-3.1\sigma$ lower than the distance ladder, independent of Planck, and that this downward pull was less apparent in some earlier analyses that used only angle-averaged BAO scale constraints rather than full anisotropic information. At the same time, the combination of BAO and CMB data also disfavors the lower values of $H_0$ preferred by the Planck high-multipole temperature power spectrum. Combining galaxy and Lyman-$\alpha$ forest (Ly$\alpha$) BAO with a precise estimate of the primordial deuterium abundance produces $H_0=66.98\pm1.18$ km s$^{-1}$ Mpc$^{-1}$ for the flat $\Lambda$CDM model. This value is completely independent of CMB anisotropy constraints and is $3.0\sigma$ lower than the latest distance ladder constraint, although $2.4\sigma$ tension also exists between the galaxy BAO and Ly$\alpha$ BAO. These results show that it is not possible to explain the $H_0$ disagreement solely with a systematic error specific to the Planck data. The fact that tensions remain even after the removal of any single data set makes this intriguing puzzle all the more challenging to resolve.

We study the scaling relations between the baryonic content and total mass of groups of galaxies, as these systems provide a unique way to examine the role of non-gravitational processes in structure formation. Using Planck and ROSAT data, we conduct detailed comparisons of the stacked thermal Sunyaev-Zel'dovich (tSZ) effect and X-ray scaling relations of galaxy groups found in the the Galaxy And Mass Assembly (GAMA) survey and the BAHAMAS hydrodynamical simulation. We use weak gravitational lensing data from the Kilo Degree Survey (KiDS) to determine the average halo mass of the studied systems. We analyse the simulation in the same way, using realistic weak lensing, X-ray, and tSZ synthetic observations. Furthermore, to keep selection biases under control, we employ exactly the same galaxy selection and group identification procedures to the observations and simulation. Applying this careful comparison, we find that the simulations are in agreement with the observations, particularly with regards to the scaling relations of the lensing and tSZ results. This finding demonstrates that hydrodynamical simulation have reached the level of realism that is required to interpret observational survey data and study the baryon physics within dark matter haloes, where analytical modelling is challenging. Finally, using simulated data, we demonstrate that our observational processing of the X-ray and tSZ signals is free of significant biases. We find that our optical group selection procedure has, however, some room for improvement.

We present results from a new incoherent-beam Fast Radio Burst (FRB) search on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder. Its large instantaneous field of view (FoV) and relative thermal insensitivity allow us to probe the ultra-bright tail of the FRB distribution, and to test a recent claim that this distribution's slope, $\alpha\equiv-\frac{\partial \log N}{\partial \log S}$, is quite small. A 256-input incoherent beamformer was deployed on the CHIME Pathfinder for this purpose. If the FRB distribution were described by a single power-law with $\alpha=0.7$, we would expect an FRB detection every few days, making this the fastest survey on sky at present. We collected 1268 hours of data, amounting to one of the largest exposures of any FRB survey, with over 2.4\,$\times$\,10$^5$\,deg$^2$\,hrs. Having seen no bursts, we have constrained the rate of extremely bright events to $<\!13$\,sky$^{-1}$\,day$^{-1}$ above $\sim$\,220$\sqrt{(\tau/\rm ms)}$ Jy\,ms for $\tau$ between 1.3 and 100\,ms, at 400--800\,MHz. The non-detection also allows us to rule out $\alpha\lesssim0.9$ with 95$\%$ confidence, after marginalizing over uncertainties in the GBT rate at 700--900\,MHz, though we show that for a cosmological population and a large dynamic range in flux density, $\alpha$ is brightness-dependent. Since FRBs now extend to large enough distances that non-Euclidean effects are significant, there is still expected to be a dearth of faint events and relative excess of bright events. Nevertheless we have constrained the allowed number of ultra-intense FRBs. While this does not have significant implications for deeper, large-FoV surveys like full CHIME and APERTIF, it does have important consequences for other wide-field, small dish experiments.

The Cosmology Large Angular Scale Surveyor (CLASS) is a four telescope array designed to characterize relic primordial gravitational waves from inflation and the optical depth to reionization through a measurement of the polarized cosmic microwave background (CMB) on the largest angular scales. The frequencies of the four CLASS telescopes, one at 38 GHz, two at 93 GHz, and one dichroic system at 145/217 GHz, are chosen to avoid spectral regions of high atmospheric emission and span the minimum of the polarized Galactic foregrounds: synchrotron emission at lower frequencies and dust emission at higher frequencies. Low-noise transition edge sensor detectors and a rapid front-end polarization modulator provide a unique combination of high sensitivity, stability, and control of systematics. The CLASS site, at 5200 m in the Chilean Atacama desert, allows for daily mapping of up to 70\% of the sky and enables the characterization of CMB polarization at the largest angular scales. Using this combination of a broad frequency range, large sky coverage, control over systematics, and high sensitivity, CLASS will observe the reionization and recombination peaks of the CMB E- and B-mode power spectra. CLASS will make a cosmic variance limited measurement of the optical depth to reionization and will measure or place upper limits on the tensor-to-scalar ratio, $r$, down to a level of 0.01 (95\% C.L.).

We present measurements of the spatial mapping between (hot) baryons and the total matter in the Universe, via the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) map from Planck and the weak gravitational lensing maps from the Red Sequence Cluster Survey (RCSLenS). The cross-correlations are performed on the map level where all the sources (including diffuse intergalactic gas) contribute to the signal. We consider two configuration-space correlation function estimators, $\xi^{ y-\kappa}$ and $\xi^ {y-\gamma_{t}}$, and a Fourier space estimator, $C_{\ell}^{y-\kappa}$, in our analysis. We detect a significant correlation out to three degrees of angular separation on the sky. Based on statistical noise only, we can report 13$\sigma$ and 17$\sigma$ detections of the cross-correlation using the configuration-space $y-\kappa$ and $y-\gamma_{t}$ estimators, respectively. Including a heuristic estimate of the sampling variance yields a detection significance of 6$\sigma$ and 8$\sigma$, respectively. A similar level of detection is obtained from the Fourier-space estimator, $C_{\ell}^{y-\kappa}$. As each estimator probes different dynamical ranges, their combination improves the significance of the detection. We compare our measurements with predictions from the cosmo-OWLS suite of cosmological hydrodynamical simulations, where different galactic feedback models are implemented. We find that a model with considerable AGN feedback that removes large quantities of hot gas from galaxy groups and WMAP-7yr best-fit cosmological parameters provides the best match to the measurements. All baryonic models in the context of a Planck cosmology over-predict the observed signal. Similar cosmological conclusions are drawn when we employ a halo model with the observed `universal' pressure profile.

The Primordial Inflation Polarization ExploreR (PIPER) is a balloon-borne telescope designed to measure the polarization of the Cosmic Microwave Background on large angular scales. PIPER will map 85% of the sky at 200, 270, 350, and 600 GHz over a series of 8 conventional balloon flights from the northern and southern hemispheres. The first science flight will use two 32x40 arrays of backshort-under-grid transition edge sensors, multiplexed in the time domain, and maintained at 100 mK by a Continuous Adiabatic Demagnetization Refrigerator. Front-end cryogenic Variable-delay Polarization Modulators provide systematic control by rotating linear to circular polarization at 3 Hz. Twin telescopes allow PIPER to measure Stokes I, Q, U, and V simultaneously. The telescope is maintained at 1.5 K in an LHe bucket dewar. Cold optics and the lack of a warm window permit sensitivity at the sky-background limit. The ultimate science target is a limit on the tensor-to-scalar ratio of r ~ 0.007, from the reionization bump to l ~ 300. PIPER's first flight will be from the Northern hemisphere, and overlap with the CLASS survey at lower frequencies. We describe the current status of the PIPER instrument.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder radio telescope is currently surveying the northern hemisphere between 400 and 800 MHz. By mapping the large scale structure of neutral hydrogen through its redshifted 21 cm line emission between $z \sim 0.8-2.5$ CHIME will contribute to our understanding of Dark Energy. Bright astrophysical foregrounds must be separated from the neutral hydrogen signal, a task which requires precise characterization of the polarized telescope beams. Using the DRAO John A. Galt 26 m telescope, we have developed a holography instrument and technique for mapping the CHIME Pathfinder beams. We report the status of the instrument and initial results of this effort.

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations, which includes different galactic feedback models, to predict the cross-correlation signal between weak gravitational lensing and the thermal Sunyaev-Zeldovich (tSZ) $y$-parameter. The predictions are compared to the recent detection reported by van Waerbeke and collaborators. The simulations reproduce the weak lensing-tSZ cross-correlation, $\xi_{y\kappa}(\theta)$, well. The uncertainty arising from different possible feedback models appears to be important on small scales only ($\theta \lesssim 10$ arcmin), while the amplitude of the correlation on all scales is sensitive to cosmological parameters that control the growth rate of structure (such as $\sigma_8$, $\Omega_m$ and $\Omega_b$). This study confirms our previous claim (in Ma et al.) that a significant proportion of the signal originates from the diffuse gas component in low-mass ($M_{\rm{halo}} \lesssim 10^{14} M_{\odot}$) clusters as well as from the region beyond the virial radius. We estimate that approximately 20$\%$ of the detected signal comes from low-mass clusters, which corresponds to about 30$\%$ of the baryon density of the Universe. The simulations also suggest that more than half of the baryons in the Universe are in the form of diffuse gas outside halos ($\gtrsim 5$ times the virial radius) which is not hot or dense enough to produce a significant tSZ signal or be observed by X-ray experiments. Finally, we show that future high-resolution tSZ-lensing cross-correlation observations will serve as a powerful tool for discriminating between different galactic feedback models.

We examine the internal consistency of the Planck 2015 cosmic microwave background (CMB) temperature anisotropy power spectrum. We show that tension exists between cosmological constant cold dark matter (LCDM) model parameters inferred from multipoles l<1000 (roughly those accessible to Wilkinson Microwave Anisotropy Probe), and from l>=1000, particularly the CDM density, Omega_ch^2, which is discrepant at 2.5 sigma for a Planck-motivated prior on the optical depth, tau=0.07+/-0.02. We find some parameter tensions to be larger than previously reported because of inaccuracy in the code used by the Planck Collaboration to generate model spectra. The Planck l>=1000 constraints are also in tension with low-redshift data sets, including Planck's own measurement of the CMB lensing power spectrum (2.4 sigma), and the most precise baryon acoustic oscillation (BAO) scale determination (2.5 sigma). The Hubble constant predicted by Planck from l>=1000, H_0=64.1+/-1.7 km/s/Mpc, disagrees with the most precise local distance ladder measurement of 73.0+/-2.4 km/s/Mpc at the 3.0 sigma level, while the Planck value from l<1000, 69.7+/-1.7 km/s/Mpc, is consistent within 1 sigma. A discrepancy between the Planck and South Pole Telescope (SPT) high-multipole CMB spectra disfavors interpreting these tensions as evidence for new physics. We conclude that the parameters from the Planck high-multipole spectrum probably differ from the underlying values due to either an unlikely statistical fluctuation or unaccounted-for systematics persisting in the Planck data.

The Cosmology Large Angular Scale Surveyor (CLASS) will measure the polarization of the Cosmic Microwave Background to search for and characterize the polarized signature of inflation. CLASS will operate from the Atacama Desert and observe $\sim$70% of the sky. A variable-delay polarization modulator (VPM) modulates the polarization at $\sim$10 Hz to suppress the 1/f noise of the atmosphere and enable the measurement of the large angular scale polarization modes. The measurement of the inflationary signal across angular scales that span both the recombination and reionization features allows a test of the predicted shape of the polarized angular power spectra in addition to a measurement of the energy scale of inflation. CLASS is an array of telescopes covering frequencies of 38, 93, 148, and 217 GHz. These frequencies straddle the foreground minimum and thus allow the extraction of foregrounds from the primordial signal. Each focal plane contains feedhorn-coupled transition-edge sensors that simultaneously detect two orthogonal linear polarizations. The use of single-crystal silicon as the dielectric for the on-chip transmission lines enables both high efficiency and uniformity in fabrication. Integrated band definition has been implemented that both controls the bandpass of the single mode transmission on the chip and prevents stray light from coupling to the detectors.

Approximately half of the Universe's baryons are in a form that has been hard to detect directly. However, the missing component can be traced through the cross-correlation of the thermal Sunyaev-Zeldovich (tSZ) effect with weak gravitational lensing. We build a model for this correlation and use it to constrain the extended baryon component, employing data from the Canada France Hawaii Lensing Survey and the {\it Planck\/} satellite. The measured correlation function is consistent with an isothermal $\beta$-model for the halo gas pressure profile, and the 1- and 2-halo terms are both detected at the 4$\sigma$ level. In addition, we measure the hydrostatic mass bias $(1-b)=0.79^{+0.07}_{-0.10}$, which is consistent with numerical simulation results and the constraints from X-ray observations. The effective temperature of the gas is found to be in the range ($7\times10^{5}$--$3 \times10^{8}$)\,K, with approximately $50\%$ of the baryons appearing to lie beyond the virial radius of the halos, consistent with current expectations for the warm-hot intergalactic medium.

We present a procedure for efficiently compressing astronomical radio data for high performance applications. Integrated, post-correlation data are first passed through a nearly lossless rounding step which compares the precision of the data to a generalized and calibration-independent form of the radiometer equation. This allows the precision of the data to be reduced in a way that has an insignificant impact on the data. The newly developed Bitshuffle lossless compression algorithm is subsequently applied. When the algorithm is used in conjunction with the HDF5 library and data format, data produced by the CHIME Pathfinder telescope is compressed to 28% of its original size and decompression throughputs in excess of 1 GB/s are obtained on a single core.

The CHIME Pathfinder is a new interferometric radio telescope that uses a hybrid FPGA/GPU FX correlator. The GPU-based X-engine of this correlator processes over 819 Gb/s of 4+4-bit complex astronomical data from N=256 inputs across a 400 MHz radio band. A software framework is presented to manage this real-time data flow, which allows each of 16 processing servers to handle 51.2 Gb/s of astronomical data, plus 8 Gb/s of ancillary data. Each server receives data in the form of UDP packets from an FPGA F-engine over the eight 10 GbE links, combines data from these packets into large (32MB-256MB) buffered frames, and transfers them to multiple GPU co-processors for correlation. The results from the GPUs are combined and normalized, then transmitted to a collection server, where they are merged into a single file. Aggressive optimizations enable each server to handle this high rate of data; allowing the efficient correlation of 25 MHz of radio bandwidth per server. The solution scales well to larger values of N by adding additional servers.

We present the design and implementation of a custom GPU-based compute cluster that provides the correlation X-engine of the CHIME Pathfinder radio telescope. It is among the largest such systems in operation, correlating 32,896 baselines (256 inputs) over 400MHz of radio bandwidth. Making heavy use of consumer-grade parts and a custom software stack, the system was developed at a small fraction of the cost of comparable installations. Unlike existing GPU backends, this system is built around OpenCL kernels running on consumer-level AMD GPUs, taking advantage of low-cost hardware and leveraging packed integer operations to double algorithmic efficiency. The system achieves the required 105TOPS in a 10kW power envelope, making it among the most power-efficient X-engines in use today.

We examine the consistency of WMAP9 and Planck data. We compare sky maps, power spectra, and inferred LCDM cosmological parameters. Residual dipoles are seen in the WMAP and Planck sky map differences, but are consistent within the uncertainties and are not large enough to explain the widely-noted differences in angular power spectra at higher l. After removing residual dipoles and galactic foregrounds, the residual difference maps exhibit a quadrupole and other large-scale systematic structure. We identify this structure as possibly originating from Planck's beam sidelobe pick-up, but note that it appears to have insignificant cosmological impact. We develop an extension of the internal linear combination technique and find features that plausibly originate in the Planck data. We examine LCDM model fits to the angular power spectra and conclude that the ~2.5% difference in the spectra at multipoles greater than l~100 are significant at the 3-5 sigma level. We revisit the analysis of WMAP's beam data and conclude that previously-derived uncertainties are robust and cannot explain the power spectrum differences. Finally, we examine the consistency of the LCDM parameters inferred from each data set taking into account the fact that both experiments observe the same sky, but cover different multipole ranges, apply different sky masks, and have different noise. While individual parameter values agree within the uncertainties, the 6 parameters taken together are discrepant at the ~6 sigma level, with chi2=56 for 6 dof (PTE = 3e-10). Of the 6 parameters, chi2 is best improved by marginalizing over Omega_c h^2, giving chi2=5.2 for 5 degrees of freedom. We find that perturbing the WMAP window function by its dominant beam error profile has little effect on Omega_c h^2, while perturbing the Planck window function by its corresponding error profile has a much greater effect on Omega_c h^2.

The determination of the Hubble constant has been a central goal in observational astrophysics for nearly 100 years. Extraordinary progress has occurred in recent years on two fronts: the cosmic distance ladder measurements at low redshift and cosmic microwave background (CMB) measurements at high redshift. The CMB is used to predict the current expansion rate through a best-fit cosmological model. Complementary progress has been made with baryon acoustic oscillation (BAO) measurements at relatively low redshifts. While BAO data do not independently determine a Hubble constant, they are important for constraints on possible solutions and checks on cosmic consistency. A precise determination of the Hubble constant is of great value, but it is more important to compare the high and low redshift measurements to test our cosmological model. Significant tension would suggest either uncertainties not accounted for in the experimental estimates, or the discovery of new physics beyond the standard model of cosmology. In this paper we examine in detail the tension between the CMB, BAO, and cosmic distance ladder data sets. We find that these measurements are consistent within reasonable statistical expectations, and we combine them to determine a best-fit Hubble constant of 69.6+/-0.7 km/s/Mpc. This value is based upon WMAP9+SPT+ACT+6dFGS+BOSS/DR11+H_0/Riess; we explore alternate data combinations in the text. The combined data constrain the Hubble constant to 1%, with no compelling evidence for new physics.

The Cosmology Large Angular Scale Surveyor (CLASS) experiment aims to map the polarization of the Cosmic Microwave Background (CMB) at angular scales larger than a few degrees. Operating from Cerro Toco in the Atacama Desert of Chile, it will observe over 65% of the sky at 38, 93, 148, and 217 GHz. In this paper we discuss the design, construction, and characterization of the CLASS 38 GHz detector focal plane, the first ever Q-band bolometric polarimeter array.

The Cosmology Large Angular Scale Surveyor (CLASS) is an experiment to measure the signature of a gravita-tional-wave background from inflation in the polarization of the cosmic microwave background (CMB). CLASS is a multi-frequency array of four telescopes operating from a high-altitude site in the Atacama Desert in Chile. CLASS will survey 70\% of the sky in four frequency bands centered at 38, 93, 148, and 217 GHz, which are chosen to straddle the Galactic-foreground minimum while avoiding strong atmospheric emission lines. This broad frequency coverage ensures that CLASS can distinguish Galactic emission from the CMB. The sky fraction of the CLASS survey will allow the full shape of the primordial B-mode power spectrum to be characterized, including the signal from reionization at low $\ell$. Its unique combination of large sky coverage, control of systematic errors, and high sensitivity will allow CLASS to measure or place upper limits on the tensor-to-scalar ratio at a level of $r=0.01$ and make a cosmic-variance-limited measurement of the optical depth to the surface of last scattering, $\tau$.

The recent BICEP2 measurement of primordial gravity waves (r = 0.2^{+0.07}_{-0.05}) appears to be in tension with the upper limit from WMAP (r<0.13 at 95% CL) and Planck (r<0.11 at 95% CL). We carefully quantify the level of tension and show that it is very significant (around 0.1% unlikely) when the observed deficit of large-scale temperature power is taken into account. We show that measurements of TE and EE power spectra in the near future will discriminate between the hypotheses that this tension is either a statistical fluke, or a sign of new physics. We also discuss extensions of the standard cosmological model that relieve the tension, and some novel ways to constrain them.

The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne cosmic microwave background (CMB) polarimeter designed to search for evidence of inflation by measuring the large-angular scale CMB polarization signal. BICEP2 recently reported a detection of B-mode power corresponding to the tensor-to-scalar ratio r = 0.2 on ~2 degree scales. If the BICEP2 signal is caused by inflationary gravitational waves (IGWs), then there should be a corresponding increase in B-mode power on angular scales larger than 18 degrees. PIPER is currently the only suborbital instrument capable of fully testing and extending the BICEP2 results by measuring the B-mode power spectrum on angular scales $\theta$ = ~0.6 deg to 90 deg, covering both the reionization bump and recombination peak, with sensitivity to measure the tensor-to-scalar ratio down to r = 0.007, and four frequency bands to distinguish foregrounds. PIPER will accomplish this by mapping 85% of the sky in four frequency bands (200, 270, 350, 600 GHz) over a series of 8 conventional balloon flights from the northern and southern hemispheres. The instrument has background-limited sensitivity provided by fully cryogenic (1.5 K) optics focusing the sky signal onto four 32x40-pixel arrays of time-domain multiplexed Transition-Edge Sensor (TES) bolometers held at 140 mK. Polarization sensitivity and systematic control are provided by front-end Variable-delay Polarization Modulators (VPMs), which rapidly modulate only the polarized sky signal at 3 Hz and allow PIPER to instantaneously measure the full Stokes vector (I, Q, U, V) for each pointing. We describe the PIPER instrument and progress towards its first flight.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark energy. We have deployed a pathfinder version of CHIME that will yield constraints on the BAO power spectrum and provide a test-bed for our calibration scheme. I will discuss the CHIME calibration requirements and describe instrumentation we are developing to meet these requirements.

A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryon acoustic oscillation (BAO) scale across this poorly probed redshift range where dark energy becomes a significant contributor to the evolution of the Universe. The instrument revives the cylinder design in radio astronomy with a wide field survey as a primary goal. Modern low-noise amplifiers and digital processing remove the necessity for the analog beamforming that characterized previous designs. The Pathfinder consists of two cylinders 37\,m long by 20\,m wide oriented north-south for a total collecting area of 1,500 square meters. The cylinders are stationary with no moving parts, and form a transit instrument with an instantaneous field of view of $\sim$100\,degrees by 1-2\,degrees. Each CHIME Pathfinder cylinder has a feedline with 64 dual polarization feeds placed every $\sim$30\,cm which Nyquist sample the north-south sky over much of the frequency band. The signals from each dual-polarization feed are independently amplified, filtered to 400-800\,MHz, and directly sampled at 800\,MSps using 8 bits. The correlator is an FX design, where the Fourier transform channelization is performed in FPGAs, which are interfaced to a set of GPUs that compute the correlation matrix. The CHIME Pathfinder is a 1/10th scale prototype version of CHIME and is designed to detect the BAO feature and constrain the distance-redshift relation.

We examine the Sunyaev--Zeldovich (SZ) effect in the 7-yr Wilkinson Microwave Anisotropy Probe (WMAP) data by cross-correlating it with the Planck Early-release Sunyaev--Zeldovich catalog. Our analysis proceeds in two parts. We first perform a stacking analysis in which the filtered WMAP data are averaged at the locations of the 175 Planck clusters. We then perform a regression analysis to compare the mean amplitude of the SZ signal, Y_{500}, in the WMAP data to the corresponding amplitude in the Planck data. The aggregate Planck clusters are detected in the 7-yr WMAP data with a signal-to-noise ratio of 16.3. In the regression analysis we find that the SZ amplitude measurements agree to better than 25%: $a=1.23 \pm 0.18$ for the fit Y^{wmap}_{500} = aY^{planck}_{500}$.

We report the first detection of a correlation between gravitational lensing by large scale structure and the thermal Sunyaev-Zeldovich (tSZ) effect. Using the mass map from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS) and a newly-constructed tSZ map from Planck, we measure a non-zero correlation between the two maps out to one degree angular separation on the sky, with an overall significance of 6 sigma. The tSZ maps are formed in a manner that removes primary cosmic microwave background fluctuations and minimizes residual contamination by galactic and extragalactic dust emission, and by CO line emission. We perform numerous tests to show that our measurement is immune to these residual contaminants. The resulting correlation function is consistent with the existence of a warm baryonic gas tracing the large scale structure with a bias b_gas. Given the shape of the lensing kernel, our signal sensitivity peaks at a redshift z~0.4, where half a degree separation on the sky corresponds to a physical scale of ~10 Mpc. The amplitude of the signal constrains the product (b_gas/1)(T_e / 0.1 keV)(n_e / 1 m^-3)=2.01\pm 0.31\pm 0.21, at redshift zero. Our study suggests that a substantial fraction of the "missing" baryons in the universe may reside in a low density warm plasma that traces dark matter.

We constrain cosmological parameters using combined measurements of the baryon acoustic oscillation (BAO) feature in the correlation function of galaxies and Ly-\alpha absorbers that together cover 0.1 < z < 2.4. The BAO position measurements alone -- without fixing the absolute sound horizon `standard ruler' length with cosmic microwave background (CMB) data -- constrain \Omega_m = 0.303 +/- 0.040 (68 per cent confidence) for a flat \Lambda CDM model, and w = -1.06^{+0.33}_{-0.32}, \Omega_m = 0.292^{+0.045}_{-0.040} for a flat wCDM model. Adding other large-scale structure (LSS) clustering constraints -- correlation function shape, the Alcock-Paczynski test and growth rate information -- to the BAO considerably tightens constraints (\Omega_m = 0.290 +/- 0.019, H_0 = 67.5 +/- 2.8 km s^{-1} Mpc^{-1}, \sigma_8 = 0.80 +/- 0.05 for \Lambda CDM, and w = -1.14 +/- 0.19 for wCDM). The LSS data mildly prefer a lower value of H_0, and a higher value of \Omega_m, than local distance ladder and type IA supernovae (SNe) measurements, respectively. While tension in the combined CMB, SNe and distance ladder data appears to be relieved by allowing w < -1, this freedom introduces tension with the LSS \sigma_8 constraint from the growth rate of matter fluctuations. The combined constraint on w from CMB, BAO and LSS clustering for a flat wCDM model is w = -1.03 +/- 0.06.

We present cosmological parameter constraints based on the final nine-year WMAP data, in conjunction with additional cosmological data sets. The WMAP data alone, and in combination, continue to be remarkably well fit by a six-parameter LCDM model. When WMAP data are combined with measurements of the high-l CMB anisotropy, the BAO scale, and the Hubble constant, the densities, Omegabh2, Omegach2, and Omega_L, are each determined to a precision of ~1.5%. The amplitude of the primordial spectrum is measured to within 3%, and there is now evidence for a tilt in the primordial spectrum at the 5sigma level, confirming the first detection of tilt based on the five-year WMAP data. At the end of the WMAP mission, the nine-year data decrease the allowable volume of the six-dimensional LCDM parameter space by a factor of 68,000 relative to pre-WMAP measurements. We investigate a number of data combinations and show that their LCDM parameter fits are consistent. New limits on deviations from the six-parameter model are presented, for example: the fractional contribution of tensor modes is limited to r<0.13 (95% CL); the spatial curvature parameter is limited to -0.0027 (+0.0039/-0.0038); the summed mass of neutrinos is <0.44 eV (95% CL); and the number of relativistic species is found to be 3.84+/-0.40 when the full data are analyzed. The joint constraint on Neff and the primordial helium abundance agrees with the prediction of standard Big Bang nucleosynthesis. We compare recent PLANCK measurements of the Sunyaev-Zel'dovich effect with our seven-year measurements, and show their mutual agreement. Our analysis of the polarization pattern around temperature extrema is updated. This confirms a fundamental prediction of the standard cosmological model and provides a striking illustration of acoustic oscillations and adiabatic initial conditions in the early universe.

We present the final nine-year maps and basic results from the WMAP mission. We provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate CMB anisotropy from foreground emission, and both types of signals are analyzed in detail. The WMAP mission has resulted in a highly constrained LCDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that Big Bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (3.84+/-0.40). The model fit also implies that the age of the universe is 13.772+/-0.059 Gyr, and the fit Hubble constant is H0 = 69.32+/-0.80 km/s/Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity reported earlier by WMAP now has high statistical significance (n_s = 0.9608+/-0.0080); and the universe is close to flat/Euclidean, Omega_k = -0.0027 (+0.0039/-0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter LCDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further. With no significant anomalies and an adequate goodness-of-fit, the inflationary flat LCDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

The Primordial Inflation Explorer (PIXIE) is an Explorer-class mission to measure the gravity-wave signature of primordial inflation through its distinctive imprint on the linear polarization of the cosmic microwave background. The instrument consists of a polarizing Michelson interferometer configured as a nulling polarimeter to measure the difference spectrum between orthogonal linear polarizations from two co-aligned beams. Either input can view the sky or a temperature-controlled absolute reference blackbody calibrator. PIXIE will map the absolute intensity and linear polarization (Stokes I, Q, and U parameters) over the full sky in 400 spectral channels spanning 2.5 decades in frequency from 30 GHz to 6 THz (1 cm to 50 um wavelength). Multi-moded optics provide background-limited sensitivity using only 4 detectors, while the highly symmetric design and multiple signal modulations provide robust rejection of potential systematic errors. The principal science goal is the detection and characterization of linear polarization from an inflationary epoch in the early universe, with tensor-to-scalar ratio r < 10^{-3} at 5 standard deviations. The rich PIXIE data set will also constrain physical processes ranging from Big Bang cosmology to the nature of the first stars to physical conditions within the interstellar medium of the Galaxy.

(Abridged) A simple six-parameter LCDM model provides a successful fit to WMAP data, both when the data are analyzed alone and in combination with other cosmological data. Even so, it is appropriate to search for any hints of deviations from the now standard model of cosmology, which includes inflation, dark energy, dark matter, baryons, and neutrinos. The cosmological community has subjected the WMAP data to extensive and varied analyses. While there is widespread agreement as to the overall success of the six-parameter LCDM model, various "anomalies" have been reported relative to that model. In this paper we examine potential anomalies and present analyses and assessments of their significance. In most cases we find that claimed anomalies depend on posterior selection of some aspect or subset of the data. Compared with sky simulations based on the best fit model, one can select for low probability features of the WMAP data. Low probability features are expected, but it is not usually straightforward to determine whether any particular low probability feature is the result of the a posteriori selection or of non-standard cosmology. We examine in detail the properties of the power spectrum with respect to the LCDM model. We examine several potential or previously claimed anomalies in the sky maps and power spectra, including cold spots, low quadrupole power, quadropole-octupole alignment, hemispherical or dipole power asymmetry, and quadrupole power asymmetry. We conclude that there is no compelling evidence for deviations from the LCDM model, which is generally an acceptable statistical fit to WMAP and other cosmological data.

(Abridged) We present WMAP seven-year observations of bright sources which are often used as calibrators at microwave frequencies. Ten objects are studied in five frequency bands (23 - 94 GHz): the outer planets (Mars, Jupiter, Saturn, Uranus and Neptune) and five fixed celestial sources (Cas A, Tau A, Cyg A, 3C274 and 3C58). The seven-year analysis of Jupiter provides temperatures which are within 1-sigma of the previously published WMAP five-year values, with slightly tighter constraints on variability with orbital phase, and limits (but no detections) on linear polarization. Scaling factors are provided which, when multiplied by the Wright Mars thermal model predictions at 350 micron, reproduce WMAP seasonally averaged observations of Mars within ~2%. An empirical model is described which fits brightness variations of Saturn due to geometrical effects and can be used to predict the WMAP observations to within 3%. Seven-year mean temperatures for Uranus and Neptune are also tabulated. Uncertainties in Uranus temperatures are 3%-4% in the 41, 61 and 94 GHz bands; the smallest uncertainty for Neptune is ~8% for the 94 GHz band. Intriguingly, the spectrum of Uranus appears to show a dip at ~30 GHz of unidentified origin, although the feature is not of high statistical significance. Flux densities for the five selected fixed celestial sources are derived from the seven-year WMAP sky maps, and are tabulated for Stokes I, Q and U, along with polarization fraction and position angle. Fractional uncertainties for the Stokes I fluxes are typically 1% to 3%. Source variability over the seven-year baseline is also estimated. Significant secular decrease is seen for Cas A and Tau A: our results are consistent with a frequency independent decrease of about 0.53% per year for Cas A and 0.22% per year for Tau A.

[Abridged] We present updated estimates of Galactic foreground emission using seven years of WMAP data. Using the power spectrum of differences between multi-frequency template-cleaned maps, we find no evidence for foreground contamination outside of the updated (KQ85y7) foreground mask. We place a 15 microKelvin upper bound on rms foreground contamination in the cleaned maps used for cosmological analysis. We find no indication in the polarization data of an extra "haze" of hard synchrotron emission from energetic electrons near the Galactic center. We provide an updated map of the cosmic microwave background (CMB) using the internal linear combination (ILC) method, updated foreground masks, and updates to point source catalogs with 62 newly detected sources. Also new are tests of the Markov chain Monte Carlo (MCMC) foreground fitting procedure against systematics in the time-stream data, and tests against the observed beam asymmetry. Within a few degrees of the Galactic plane, WMAP total intensity data show a rapidly steepening spectrum from 20-40 GHz, which may be due to emission from spinning dust grains, steepening synchrotron, or other effects. Comparisons are made to a 1-degree 408 MHz map (Haslam et al.) and the 11-degree ARCADE 2 data (Singal et al.). We find that spinning dust or steepening synchrotron models fit the combination of WMAP and 408 MHz data equally well. ARCADE data appear inconsistent with the steepening synchrotron model, and consistent with the spinning dust model, though some discrepancies remain regarding the relative strength of spinning dust emission. More high-resolution data in the 10-40 GHz range would shed much light on these issues.

(Abridged) The 7-year WMAP data and improved astrophysical data rigorously test the standard cosmological model and its extensions. By combining WMAP with the latest distance measurements from BAO and H0 measurement, we determine the parameters of the simplest LCDM model. The power-law index of the primordial power spectrum is n_s=0.968+-0.012, a measurement that excludes the scale-invariant spectrum by 99.5%CL. The other parameters are also improved from the 5-year results. Notable examples of improved parameters are the total mass of neutrinos, sum(m_nu)<0.58eV, and the effective number of neutrino species, N_eff=4.34+0.86-0.88. We detect the effect of primordial helium on the temperature power spectrum and provide a new test of big bang nucleosynthesis. We detect, and show on the map for the first time, the tangential and radial polarization patterns around hot and cold spots of temperature fluctuations, an important test of physical processes at z=1090 and the dominance of adiabatic scalar fluctuations. With the 7-year TB power spectrum, the limit on a rotation of the polarization plane due to potential parity-violating effects has improved to Delta(alpha)=-1.1+-1.4(stat)+-1.5(syst) degrees. We report significant detections of the SZ effect at the locations of known clusters of galaxies. The measured SZ signal agrees well with the expected signal from the X-ray data. However, it is a factor of 0.5 to 0.7 times the predictions from "universal profile" of Arnaud et al., analytical models, and hydrodynamical simulations. We find, for the first time in the SZ effect, a significant difference between the cooling-flow and non-cooling-flow clusters (or relaxed and non-relaxed clusters), which can explain some of the discrepancy. This lower amplitude is consistent with the lower-than-theoretically-expected SZ power spectrum recently measured by the South Pole Telescope collaboration.

(Abridged) We present the angular power spectra derived from the 7-year maps and discuss the cosmological conclusions that can be inferred from WMAP data alone. The third acoustic peak in the TT spectrum is now well measured by WMAP. In the context of a flat LambdaCDM model, this improvement allows us to place tighter constraints on the matter density from WMAP data alone, and on the epoch of matter-radiation equality, The temperature-polarization (TE) spectrum is detected in the 7-year data with a significance of 20 sigma, compared to 13 sigma with the 5-year data. The low-l EE spectrum, a measure of the optical depth due to reionization, is detected at 5.5 sigma significance when averaged over l = 2-7. The BB spectrum, an important probe of gravitational waves from inflation, remains consistent with zero. The upper limit on tensor modes from polarization data alone is a factor of 2 lower with the 7-year data than it was using the 5-year data (Komatsu et al. 2010). We test the parameter recovery process for bias and find that the scalar spectral index, ns, is biased high, but only by 0.09 sigma, while the remaining parameters are biased by < 0.15 sigma. The improvement in the third peak measurement leads to tighter lower limits from WMAP on the number of relativistic degrees of freedom (e.g., neutrinos) in the early universe: Neff > 2.7 (95% CL). Also, using WMAP data alone, the primordial helium mass fraction is found to be YHe = 0.28+0.14-0.15, and with data from higher-resolution CMB experiments included, we now establish the existence of pre-stellar helium at > 3 sigma (Komatsu et al. 2010).

(Abridged) New full sky temperature and polarization maps based on seven years of data from WMAP are presented. The new results are consistent with previous results, but have improved due to reduced noise from the additional integration time, improved knowledge of the instrument performance, and improved data analysis procedures. The improvements are described in detail. The seven year data set is well fit by a minimal six-parameter flat Lambda-CDM model. The parameters for this model, using the WMAP data in conjunction with baryon acoustic oscillation data from the Sloan Digital Sky Survey and priors on H_0 from Hubble Space Telescope observations, are: Omega_bh^2 = 0.02260 +-0.00053, Omega_ch^2 = 0.1123 +-0.0035, Omega_Lambda = 0.728 +0.015 -0.016, n_s = 0.963 +-0.012, tau = 0.087 +-0.014 and sigma_8 = 0.809 +-0.024 (68 % CL uncertainties). The temperature power spectrum signal-to-noise ratio per multipole is greater that unity for multipoles < 919, allowing a robust measurement of the third acoustic peak. This measurement results in improved constraints on the matter density, Omega_mh^2 = 0.1334 +0.0056 -0.0055, and the epoch of matter- radiation equality, z_eq = 3196 +134 -133, using WMAP data alone. The new WMAP data, when combined with smaller angular scale microwave background anisotropy data, results in a 3 sigma detection of the abundance of primordial Helium, Y_He = 0.326 +-0.075.The power-law index of the primordial power spectrum is now determined to be n_s = 0.963 +-0.012, excluding the Harrison-Zel'dovich-Peebles spectrum by >3 sigma. These new WMAP measurements provide important tests of Big Bang cosmology.

We describe a sampling method to estimate the polarized CMB signal from observed maps of the sky. We use a Metropolis-within-Gibbs algorithm to estimate the polarized CMB map, containing Q and U Stokes parameters at each pixel, and its covariance matrix. These can be used as inputs for cosmological analyses. The polarized sky signal is parameterized as the sum of three components: CMB, synchrotron emission, and thermal dust emission. The polarized Galactic components are modeled with spatially varying power law spectral indices for the synchrotron, and a fixed power law for the dust, and their component maps are estimated as by-products. We apply the method to simulated low resolution maps with pixels of side 7.2 degrees, using diagonal and full noise realizations drawn from the WMAP noise matrices. The CMB maps are recovered with goodness of fit consistent with errors. Computing the likelihood of the E-mode power in the maps as a function of optical depth to reionization, tau, for fixed temperature anisotropy power, we recover tau=0.091+-0.019 for a simulation with input tau=0.1, and mean tau=0.098 averaged over 10 simulations. A `null' simulation with no polarized CMB signal has maximum likelihood consistent with tau=0. The method is applied to the five-year WMAP data, using the K, Ka, Q and V channels. We find tau=0.090+-0.019, compared to tau=0.086+-0.016 from the template-cleaned maps used in the primary WMAP analysis. The synchrotron spectral index, beta, averaged over high signal-to-noise pixels with standard deviation sigma(beta)<0.25, but excluding ~6% of the sky masked in the Galactic plane, is -3.03+-0.04. This estimate does not vary significantly with Galactic latitude, although includes an informative prior.

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.

We present a new estimate of foreground emission in the WMAP data, using a Markov chain Monte Carlo (MCMC) method. The new technique delivers maps of each foreground component for a variety of foreground models, error estimates of the uncertainty of each foreground component, and provides an overall goodness-of-fit measurement. The resulting foreground maps are in broad agreement with those from previous techniques used both within the collaboration and by other authors. We find that for WMAP data, a simple model with power-law synchrotron, free-free, and thermal dust components fits 90% of the sky with a reduced chi-squared of 1.14. However, the model does not work well inside the Galactic plane. The addition of either synchrotron steepening or a modified spinning dust model improves the fit. This component may account for up to 14% of the total flux at Ka-band (33 GHz). We find no evidence for foreground contamination of the CMB temperature map in the 85% of the sky used for cosmological analysis.

This paper focuses on cosmological constraints derived from analysis of WMAP data alone. A simple LCDM cosmological model fits the five-year WMAP temperature and polarization data. The basic parameters of the model are consistent with the three-year data and now better constrained: Omega_b h^2 = 0.02273+-0.00062, Omega_c h^2 = 0.1099+-0.0062, Omega_L = 0.742+-0.030, n_s = 0.963+0.014- 0.015, tau = 0.087+-0.017, sigma_8 = 0.796+-0.036. With five years of polarization data, we have measured the optical depth to reionization, tau>0, at 5 sigma significance. The redshift of an instantaneous reionization is constrained to be z_reion = 11.0+-1.4 with 68% confidence. This excludes a sudden reionization of the universe at z=6 at more than 3.5 sigma significance, suggesting that reionization was an extended process. Using two methods for polarized foreground cleaning we get consistent estimates for the optical depth, indicating an error due to foreground treatment of tau~0.01. This cosmological model also fits small-scale CMB data, and a range of astronomical data measuring the expansion rate and clustering of matter in the universe. We find evidence for the first time in the CMB power spectrum for a non-zero cosmic neutrino background, or a background of relativistic species, with the standard three light neutrino species preferred over the best-fit LCDM model with N_eff=0 at >99.5% confidence, and N_eff > 2.3 (95% CL) when varied. The five-year WMAP data improve the upper limit on the tensor-to-scalar ratio to r < 0.43 (95% CL), for power-law models. With longer integration we find no evidence for a running spectral index, with dn_s/dlnk = -0.037+-0.028.

Cosmology and other scientific results from the WMAP mission require an accurate knowledge of the beam patterns in flight. While the degree of beam knowledge for the WMAP one-year and three-year results was unprecedented for a CMB experiment, we have significantly improved the beam determination as part of the five-year data release. Physical optics fits are done on both the A and the B sides for the first time. The cutoff scale of the fitted distortions on the primary mirror is reduced by a factor of ~2 from previous analyses. These changes enable an improvement in the hybridization of Jupiter data with beam models, which is optimized with respect to error in the main beam solid angle. An increase in main-beam solid angle of ~1% is found for the V2 and W1-W4 differencing assemblies. Although the five-year results are statistically consistent with previous ones, the errors in the five-year beam transfer functions are reduced by a factor of ~2 as compared to the three-year analysis. We present radiometry of the planet Jupiter as a test of the beam consistency and as a calibration standard; for an individual differencing assembly, errors in the measured disk temperature are ~0.5%.

We present new full-sky temperature and polarization maps in five frequency bands from 23 to 94 GHz, based on data from the first five years of the WMAP sky survey. The five-year maps incorporate several improvements in data processing made possible by the additional years of data and by a more complete analysis of the instrument calibration and in-flight beam response. We present several new tests for systematic errors in the polarization data and conclude that Ka band data (33 GHz) is suitable for use in cosmological analysis, after foreground cleaning. This significantly reduces the overall polarization uncertainty. With the 5 year WMAP data, we detect no convincing deviations from the minimal 6-parameter LCDM model: a flat universe dominated by a cosmological constant, with adiabatic and nearly scale-invariant Gaussian fluctuations. Using WMAP data combined with measurements of Type Ia supernovae and Baryon Acoustic Oscillations, we find (68% CL uncertainties): Omega_bh^2 = 0.02267 \pm 0.00059, Omega_ch^2 = 0.1131 \pm 0.0034, Omega_Lambda = 0.726 \pm 0.015, n_s = 0.960 \pm 0.013, tau = 0.084 \pm 0.016, and Delta_R^2 = (2.445 \pm 0.096) x 10^-9. From these we derive: sigma_8 = 0.812 \pm 0.026, H_0 = 70.5 \pm 1.3 km/s/Mpc, z_{reion} = 10.9 \pm 1.4, and t_0 = 13.72 \pm 0.12 Gyr. The new limit on the tensor-to-scalar ratio is r < 0.22 (95% CL). We obtain tight, simultaneous limits on the (constant) dark energy equation of state and spatial curvature: -0.14 < 1+w < 0.12 and -0.0179 < Omega_k < 0.0081 (both 95% CL). The number of relativistic degrees of freedom (e.g. neutrinos) is found to be N_{eff} = 4.4 \pm 1.5, consistent with the standard value of 3.04. Models with N_{eff} = 0 are disfavored at >99.5% confidence.

We present the list of point sources found in the WMAP 5-year maps. The technique used in the first-year and three-year analysis now finds 390 point sources, and the five-year source catalog is complete for regions of the sky away from the galactic plane to a 2 Jy limit, with SNR > 4.7 in all bands in the least covered parts of the sky. The noise at high frequencies is still mainly radiometer noise, but at low frequencies the CMB anisotropy is the largest uncertainty. A separate search of CMB-free V-W maps finds 99 sources of which all but one can be identified with known radio sources. The sources seen by WMAP are not strongly polarized. Many of the WMAP sources show significant variability from year to year, with more than a 2:1 range between the minimum and maximum fluxes.

We present the temperature and polarization angular power spectra of the cosmic microwave background (CMB) derived from the first 5 years of WMAP data. The 5-year temperature (TT) spectrum is cosmic variance limited up to multipole l=530, and individual l-modes have S/N>1 for l<920. The best fitting six-parameter LambdaCDM model has a reduced chi^2 for l=33-1000 of chi^2/nu=1.06, with a probability to exceed of 9.3%. There is now significantly improved data near the third peak which leads to improved cosmological constraints. The temperature-polarization correlation (TE) is seen with high significance. After accounting for foreground emission, the low-l reionization feature in the EE power spectrum is preferred by \Delta\chi^2=19.6 for optical depth tau=0.089 by the EE data alone, and is now largely cosmic variance limited for l=2-6. There is no evidence for cosmic signal in the BB, TB, or EB spectra after accounting for foreground emission. We find that, when averaged over l=2-6, l(l+1)C^{BB}_l/2\pi < 0.15 uK^2 (95% CL).

(Abridged) The WMAP 5-year data strongly limit deviations from the minimal LCDM model. We constrain the physics of inflation via Gaussianity, adiabaticity, the power spectrum shape, gravitational waves, and spatial curvature. We also constrain the properties of dark energy, parity-violation, and neutrinos. We detect no convincing deviations from the minimal model. The parameters of the LCDM model, derived from WMAP combined with the distance measurements from the Type Ia supernovae (SN) and the Baryon Acoustic Oscillations (BAO), are: Omega_b=0.0456+-0.0015, Omega_c=0.228+-0.013, Omega_Lambda=0.726+-0.015, H_0=70.5+-1.3 km/s/Mpc, n_s=0.960+-0.013, tau=0.084+-0.016, and sigma_8=0.812+-0.026. With WMAP+BAO+SN, we find the tensor-to-scalar ratio r<0.22 (95% CL), and n_s>1 is disfavored regardless of r. We obtain tight, simultaneous limits on the (constant) equation of state of dark energy and curvature. We provide a set of "WMAP distance priors," to test a variety of dark energy models. We test a time-dependent w with a present value constrained as -0.33<1+w_0<0.21 (95% CL). Temperature and matter fluctuations obey the adiabatic relation to within 8.9% and 2.1% for the axion and curvaton-type dark matter, respectively. The TE and EB spectra constrain cosmic parity-violation. We find the limit on the total mass of neutrinos, sum(m_nu)<0.67 eV (95% CL), which is free from the uncertainty in the normalization of the large-scale structure data. The effective number of neutrino species is constrained as N_{eff} = 4.4+-1.5 (68%), consistent with the standard value of 3.04. Finally, limits on primordial non-Gaussianity are -9<f_{NL}^{local}<111 and -151<f_{NL}^{equil}<253 (95% CL) for the local and equilateral models, respectively.

We present a full-sky model of polarized Galactic microwave emission based on three years of observations by the Wilkinson Microwave Anisotropy Probe (WMAP) at frequencies from 23 to 94 GHz. The model compares maps of the Stokes Q and U components from each of the 5 WMAP frequency bands in order to separate synchrotron from dust emission, taking into account the spatial and frequency dependence of the synchrotron and dust components. This simple two-component model of the interstellar medium accounts for at least 97% of the polarized emission in the WMAP maps of the microwave sky. Synchrotron emission dominates the polarized foregrounds at frequencies below 50 GHz, and is comparable to the dust contribution at 65 GHz. The spectral index of the synchrotron component, derived solely from polarization data, is -3.2 averaged over the full sky, with a modestly flatter index on the Galactic plane. The synchrotron emission has mean polarization fraction 2--4% in the Galactic plane and rising to over 20% at high latitude, with prominent features such as the North Galactic Spur more polarized than the diffuse component. Thermal dust emission has polarization fraction 1% near the Galactic center, rising to 6% at the anti-center. Diffuse emission from high-latitude dust is also polarized with mean fractional polarization 0.036 +/- 0.011.

A simple cosmological model with only six parameters (matter density, Omega_m h^2, baryon density, Omega_b h^2, Hubble Constant, H_0, amplitude of fluctuations, sigma_8, optical depth, tau, and a slope for the scalar perturbation spectrum, n_s) fits not only the three year WMAP temperature and polarization data, but also small scale CMB data, light element abundances, large-scale structure observations, and the supernova luminosity/distance relationship. Using WMAP data only, the best fit values for cosmological parameters for the power-law flat LCDM model are (Omega_m h^2, Omega_b h^2, h, n_s, tau, sigma_8) = 0.1277+0.0080-0.0079, 0.02229+-0.00073, 0.732+0.031-0.032, 0.958+-0.016, 0.089+-0.030, 0.761+0.049-0.048). The three year data dramatically shrink the allowed volume in this six dimensional parameter space. Assuming that the primordial fluctuations are adiabatic with a power law spectrum, the WMAP data_alone_ require dark matter, and favor a spectral index that is significantly less than the Harrison-Zel'dovich-Peebles scale-invariant spectrum (n_s=1, r=0). Models that suppress large-scale power through a running spectral index or a large-scale cut-off in the power spectrum are a better fit to the WMAP and small scale CMB data than the power-law LCDM model; however, the improvement in the fit to the WMAP data is only Delta chi^2 = 3 for 1 extra degree of freedom. The combination of WMAP and other astronomical data yields significant constraints on the geometry of the universe, the equation of state of the dark energy, the gravitational wave energy density, and neutrino properties. Consistent with the predictions of simple inflationary theories, we detect no significant deviations from Gaussianity in the CMB maps.

The WMAP satellite has completed 3 years of observations of the cosmic microwave background radiation. The 3-year data products include several sets of full sky maps of the Stokes I, Q and U parameters in 5 frequency bands, spanning 23 to 94 GHz, and supporting items, such as beam window functions and noise covariance matrices. The processing used to produce the current sky maps and supporting products represents a significant advancement over the first year analysis, and is described herein. Improvements to the pointing reconstruction, radiometer gain modeling, window function determination and radiometer spectral noise parametrization are presented. A detailed description of the updated data processing that produces maximum likelihood sky map estimates is presented, along with the methods used to produce reduced resolution maps and corresponding noise covariance matrices. Finally two methods used to evaluate the noise of the full resolution sky maps are presented along with several representative year-to-year null tests, demonstrating that sky maps produced from data from different observational epochs are consistent.

We present new full-sky temperature maps in five frequency bands from 23 to 94 GHz, based on the first three years of the WMAP sky survey. The new maps, which are consistent with the first-year maps and more sensitive, incorporate improvements in data processing made possible by the additional years of data and by a more complete analysis of the polarization signal. These include refinements in the gain calibration and beam response models. We employ two forms of multi-frequency analysis to separate astrophysical foreground signals from the CMB, each of which improves on our first-year analyses. First, we form an improved 'Internal Linear Combination' map, based solely on WMAP data, by adding a bias correction step and by quantifying residual uncertainties in the resulting map. Second, we fit and subtract new spatial templates that trace Galactic emission; in particular, we now use low-frequency WMAP data to trace synchrotron emission. The WMAP point source catalog is updated to include 115 new sources. We derive the angular power spectrum of the temperature anisotropy using a hybrid approach that combines a maximum likelihood estimate at low l (large angular scales) with a quadratic cross-power estimate for l>30. Our best estimate of the CMB power spectrum is derived by averaging cross-power spectra from 153 statistically independent channel pairs. The combined spectrum is cosmic variance limited to l=400, and the signal-to-noise ratio per l-mode exceeds unity up to l=850. The first two acoustic peaks are seen at l=220.8 +- 0.7 and l=530.9 +- 3.8, respectively, while the first two troughs are seen at l=412.4 +- 1.9 and l=675.1 +- 11.1, respectively. The rise to the third peak is unambiguous; when the WMAP data are combined with higher resolution CMB measurements, the existence of a third acoustic peak is well established.

The Wilkinson Microwave Anisotropy Probe WMAP has mapped the entire sky in five frequency bands between 23 and 94 GHz with polarization sensitive radiometers. We present three-year full-sky maps of the polarization and analyze them for foreground emission and cosmological implications. These observations open up a new window for understanding the universe. WMAP observes significant levels of polarized foreground emission due to both Galactic synchrotron radiation and thermal dust emission. The least contaminated channel is at 61 GHz. Informed by a model of the Galactic foreground emission, we subtract the foreground emission from the maps. In the foreground corrected maps, for l=2-6, we detect l(l+1) C_l^{EE} / (2 pi) = 0.086 +-0.029 microkelvin^2. This is interpreted as the result of rescattering of the CMB by free electrons released during reionization and corresponds to an optical depth of tau = 0.10 +- 0.03. We see no evidence for B-modes, limiting them to l(l+1) C_l^{BB} / (2 pi) = -0.04 +- 0.03 microkelvin^2. We find that the limit from the polarization signals alone is r<2.2 (95% CL) corresponding to a limit on the cosmic density of gravitational waves of Omega_{GW}h^2 < 5 times 10^{-12}. From the full WMAP analysis, we find r<0.55 (95% CL) corresponding to a limit of Omega_{GW}h^2 < 10^{-12} (95% CL).

The Primordial Anisotropy Polarization Pathfinder Array (PAPPA) is a balloon-based instrument to measure the polarization of the cosmic microwave background and search for the signal from gravity waves excited during an inflationary epoch in the early universe. PAPPA will survey a 20 x 20 deg patch at the North Celestial Pole using 32 pixels in 3 passbands centered at 89, 212, and 302 GHz. Each pixel uses MEMS switches in a superconducting microstrip transmission line to combine the phase modulation techniques used in radio astronomy with the sensitivity of transition-edge superconducting bolometers. Each switched circuit modulates the incident polarization on a single detector, allowing nearly instantaneous characterization of the Stokes I, Q, and U parameters. We describe the instrument design and status.

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.

We describe our methodology for comparing the WMAP measurements of the cosmic microwave background (CMB) and other complementary data sets to theoretical models. The unprecedented quality of the WMAP data, and the tight constraints on cosmological parameters that are derived, require a rigorous analysis so that the approximations made in the modeling do not lead to significant biases. We describe our use of the likelihood function to characterize the statistical properties of the microwave background sky. We outline the use of the Monte Carlo Markov Chains to explore the likelihood of the data given a model to determine the best fit cosmological parameters and their uncertainties. We add to the WMAP data the l>~700 CBI and ACBAR measurements of the CMB, the galaxy power spectrum at z~0 obtained from the 2dF galaxy redshift survey (2dFGRS), and the matter power spectrum at z~3 as measured with the Ly alpha forest. These last two data sets complement the CMB measurements by probing the matter power spectrum of the nearby universe. Combining CMB and 2dFGRS requires that we include in our analysis a model for galaxy bias, redshift distortions, and the non-linear growth of structure. We show how the statistical and systematic uncertainties in the model and the data are propagated through the full analysis.

The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps significant at more than 10 standard deviations. The correlations are present in all WMAP frequency bands with similar amplitude from 23 to 94 GHz, and are consistent with a superposition of a CMB signal with a weak foreground. The fitted CMB component is robust against different data combinations and fitting techniques. On small angular scales theta < 5 deg, the WMAP data show the temperature-polarization correlation expected from adiabatic perturbations in the temperature power spectrum. The data for l > 20 agree well with the signal predicted solely from the temperature power spectra, with no additional free parameters. We detect excess power on large angular scales (theta > 10 deg) compared to predictions based on the temperature power spectra alone. The excess power is well described by reionization at redshift 11 < z_r < 30 at 95% confidence, depending on the ionization history. A model-independent fit to reionization optical depth yields results consistent with the best-fit LambdaCDM model, with best fit value tau = 0.17 +- 0.04 at 68% confidence, including systematic and foreground uncertainties. This value is larger than expected given the detection of a Gunn-Peterson trough in the absorption spectra of distant quasars, and implies that the universe has a complex ionization history: WMAP has detected the signal from an early epoch of reionization.

The CMB has distinct peaks in both its temperature angular power spectrum (TT) and temperature-polarization cross-power spectrum (TE). From the WMAP data we find the first peak in the temperature spectrum at l = 220.1 +- 0.8 with an amplitude of 74.7 +- 0.5 microK; the first trough at l = 411.7 +- 3.5 with an amplitude of 41.0 +- 0.5 microK; and the second peak at l = 546 +- 10 with an amplitude of 48.8 +- 0.9 microK. The TE spectrum has an antipeak at l = 137 +- 9 with a cross-power of -35 +- 9 microK^2, and a peak at l = 329 +- 19 with cross-power 105 +- 18 microK^2. All uncertainties are 1 sigma and include calibration and beam errors. An intuition for how the data determine the cosmological parameters may be gained by limiting one's attention to a subset of parameters and their effects on the peak characteristics. We interpret the peaks in the context of a flat adiabatic LambdaCDM model with the goal of showing how the cosmic baryon density, Omega_b h^2, matter density, Omega_m h^2, scalar index, n_s, and age of the universe are encoded in their positions and amplitudes. To this end, we introduce a new scaling relation for the TE antipeak-to-peak amplitude ratio and recompute known related scaling relations for the TT spectrum in light of the WMAP data. From the scaling relations, we show that WMAP's tight bound on Omega_b h^2 is intimately linked to its robust detection of the first and second peaks of the TT spectrum.

Knowledge of the beam profiles is of critical importance for interpreting data from cosmic microwave background experiments. In this paper, we present the characterization of the in-flight optical response of the WMAP satellite. The main beam intensities have been mapped to < -30 dB of their peak values by observing Jupiter with the satellite in the same observing mode as for CMB observations. The beam patterns closely follow the pre-launch expectations. The full width at half maximum is a function of frequency and ranges from 0.82 degrees at 23 GHz to 0.21 degrees at 94 GHz; however, the beams are not Gaussian. We present: (a) the beam patterns for all ten differential radiometers and show that the patterns are substantially independent of polarization in all but the 23 GHz channel; (b) the effective symmetrized beam patterns that result from WMAP's compound spin observing pattern; (c) the effective window functions for all radiometers and the formalism for propagating the window function uncertainty; and (d) the conversion factor from point source flux to antenna temperature. A summary of the systematic uncertainties, which currently dominate our knowledge of the beams, is also presented. The constancy of Jupiter's temperature within a frequency band is an essential check of the optical system. The tests enable us to report a calibration of Jupiter to 1-3% accuracy relative to the CMB dipole.

We present limits to the amplitude of non-Gaussian primordial fluctuations in the WMAP 1-year cosmic microwave background sky maps. A non-linear coupling parameter, f_NL, characterizes the amplitude of a quadratic term in the primordial potential. We use two statistics: one is a cubic statistic which measures phase correlations of temperature fluctuations after combining all configurations of the angular bispectrum. The other uses the Minkowski functionals to measure the morphology of the sky maps. Both methods find the WMAP data consistent with Gaussian primordial fluctuations and establish limits, -58<f_NL<134, at 95% confidence. There is no significant frequency or scale dependence of f_NL. The WMAP limit is 30 times better than COBE, and validates that the power spectrum can fully characterize statistical properties of CMB anisotropy in the WMAP data to high degree of accuracy. Our results also validate the use of a Gaussian theory for predicting the abundance of clusters in the local universe. We detect a point-source contribution to the bispectrum at 41 GHz, b_src = (9.5+-4.4) X 1e-5 uK^3 sr^2, which gives a power spectrum from point sources of c_src = (15+-6) X 1e-3 uK^2 sr in thermodynamic temperature units. This value agrees well with independent estimates of source number counts and the power spectrum at 41 GHz, indicating that b_src directly measures residual source contributions.

WMAP precision data enables accurate testing of cosmological models. We find that the emerging standard model of cosmology, a flat Lambda-dominated universe seeded by nearly scale-invariant adiabatic Gaussian fluctuations, fits the WMAP data. With parameters fixed only by WMAP data, we can fit finer scale CMB measurements and measurements of large scle structure (galaxy surveys and the Lyman alpha forest). This simple model is also consistent with a host of other astronomical measurements. We then fit the model parameters to a combination of WMAP data with other finer scale CMB experiments (ACBAR and CBI), 2dFGRS measurements and Lyman alpha forest data to find the model's best fit cosmological parameters: h=0.71+0.04-0.03, Omega_b h^2=0.0224+-0.0009, Omega_m h^2=0.135+0.008-0.009, tau=0.17+-0.06, n_s(0.05/Mpc)=0.93+-0.03, and sigma_8=0.84+-0.04. WMAP's best determination of tau=0.17+-0.04 arises directly from the TE data and not from this model fit, but they are consistent. These parameters imply that the age of the universe is 13.7+-0.2 Gyr. The data favors but does not require a slowly varying spectral index. By combining WMAP data with other astronomical data sets, we constrain the geometry of the universe, Omega_tot = 1.02 +- 0.02, the equation of state of the dark energy w < -0.78 (95% confidence limit assuming w >= -1), and the energy density in stable neutrinos, Omega_nu h^2 < 0.0076 (95% confidence limit). For 3 degenerate neutrino species, this limit implies that their mass is less than 0.23 eV (95% confidence limit). The WMAP detection of early reionization rules out warm dark matter.

We present full sky microwave maps in five bands (23 to 94 GHz) from the WMAP first year sky survey. Calibration errors are <0.5% and the low systematic error level is well specified. The 2<l<900 anisotropy power spectrum is cosmic variance limited for l<354 with a signal-to-noise ratio >1 per mode to l=658. The temperature-polarization cross-power spectrum reveals both acoustic features and a large angle correlation from reionization. The optical depth of reionization is 0.17 +/- 0.04, which implies a reionization epoch of 180+220-80 Myr (95% CL) after the Big Bang at a redshift of 20+10-9 (95% CL) for a range of ionization scenarios. This early reionization is incompatible with the presence of a significant warm dark matter density. The age of the best-fit universe is 13.7 +/- 0.2 Gyr old. Decoupling was 379+8-7 kyr after the Big Bang at a redshift of 1089 +/- 1. The thickness of the decoupling surface was dz=195 +/- 2. The matter density is Omega_m h^2 = 0.135 +0.008 -0.009, the baryon density is Omega_b h^2 = 0.0224 +/- 0.0009, and the total mass-energy of the universe is Omega_tot = 1.02 +/- 0.02. The spectral index of scalar fluctuations is fit as n_s = 0.93 +/- 0.03 at wavenumber k_0 = 0.05 Mpc^-1, with a running index slope of dn_s/d ln k = -0.031 +0.016 -0.018 in the best-fit model. This flat universe model is composed of 4.4% baryons, 22% dark matter and 73% dark energy. The dark energy equation of state is limited to w<-0.78 (95% CL). Inflation theory is supported with n_s~1, Omega_tot~1, Gaussian random phases of the CMB anisotropy, and superhorizon fluctuations. An admixture of isocurvature modes does not improve the fit. The tensor-to-scalar ratio is r(k_0=0.002 Mpc^-1)<0.90 (95% CL).

Full sky maps are made in five microwave frequency bands to separate the temperature anisotropy of the CMB from foreground emission. We define masks that excise regions of high foreground emission. The effectiveness of template fits to remove foreground emission from the WMAP data is examined. These efforts result in a CMB map with minimal contamination and a demonstration that the WMAP CMB power spectrum is insensitive to residual foreground emission. We construct a model of the Galactic emission components. We find that the Milky Way resembles other normal spiral galaxies between 408 MHz and 23 GHz, with a synchrotron spectral index that is flattest (beta ~ -2.5) near star-forming regions, especially in the plane, and steepest (beta ~ -3) in the halo. The significant synchrotron index steepening out of the plane suggests a diffusion process in which the halo electrons are trapped in the Galactic potential long enough to suffer synchrotron and inverse Compton energy losses and hence a spectral steepening. The synchrotron index is steeper in the WMAP bands than in lower frequency radio surveys, with a spectral break near 20 GHz to beta < -3. The modeled thermal dust spectral index is also steep in the WMAP bands, with beta ~ 2.2. Microwave and H alpha measurements of the ionized gas agree. Spinning dust emission is limited to < ~5% of the Ka-band foreground emission. A catalog of 208 point sources is presented. Derived source counts suggest a contribution to the anisotropy power from unresolved sources of (15.0 +- 1.4) 10^{-3} microK^2 sr at Q-band and negligible levels at V-band and W-band.

We confront predictions of inflationary scenarios with the WMAP data, in combination with complementary small-scale CMB measurements and large-scale structure data. The WMAP detection of a large-angle anti-correlation in the temperature--polarization cross-power spectrum is the signature of adiabatic superhorizon fluctuations at the time of decoupling. The WMAP data are described by pure adiabatic fluctuations: we place an upper limit on a correlated CDM isocurvature component. Using WMAP constraints on the shape of the scalar power spectrum and the amplitude of gravity waves, we explore the parameter space of inflationary models that is consistent with the data. We place limits on inflationary models; for example, a minimally-coupled lambda phi^4 is disfavored at more than 3-sigma using WMAP data in combination with smaller scale CMB and large scale structure survey data. The limits on the primordial parameters using WMAP data alone are: n_s(k_0=0.002 Mpc^{-1})=1.20_{-0.11}^{+0.12}, dn/dlnk=-0.077^{+0.050}_{- 0.052}, A(k_0=0.002 Mpc}^{-1})=0.71^{+0.10}_{-0.11} (68% CL), and r(k_0=0.002 Mpc^{-1})<1.28 (95% CL).

The first-year WMAP data, in combination with any one of a number of other cosmic probes, show that we live in a flat \Lambda-dominated CDM universe with \Omega_m ~ 0.27 and \Omega_\Lambda ~ 0.73. In this model the late-time action of the dark energy, through the integrated Sachs-Wolfe effect, should produce CMB anisotropies correlated with matter density fluctuations at z<2 (Crittenden & Turok 1996). The measurement of such a signal is an important independent check of the model. We cross-correlate the NRAO VLA Sky Survey radio source catalog (Condon et al. 1998) with the WMAP data in search of this signal, and see indications of the expected correlation. Assuming a flat \Lambda-CDM cosmology, we find \Omega_\Lambda>0 (95% CL, statistical errors only) with the peak of the likelihood at \Omega_\Lambda=0.68, consistent with the preferred WMAP value. A closed model with \Omega_m=1.28, h=0.33, and no dark energy component (\Omega_\Lambda=0), marginally consistent with the WMAP CMB TT angular power spectrum, would produce an anti-correlation between the matter distribution and the CMB. Our analysis of the cross-correlation of the WMAP data with the NVSS catalog rejects this cosmology at the 3\sigma level.

The Microwave Anisotropy Probe (MAP) satellite, launched June 30, 2001, will produce full sky maps of the cosmic microwave background radiation in 5 frequency bands spanning 20 - 106 GHz. MAP contains 20 differential radiometers built with High Electron Mobility Transistor (HEMT) amplifiers with passively cooled input stages. The design and test techniques used to evaluate and minimize systematic errors and the pre-launch performance of the radiometers for all five bands are presented.

The WMAP satellite has completed one year of measurements of the Cosmic Microwave Background (CMB) radiation using 20 differential high-electron-mobility-transistor (HEMT) based radiometers. All the radiometers are functioning nominally, and characterizations of the on-orbit radiometer performance are presented, with an emphasis on properties that are required for the production of sky maps from the time ordered data. A radiometer gain model, used to smooth and interpolate the CMB dipole gain measurements is also presented. No degradation in the sensitivity of any of the radiometers has been observed during the first year of observations.

Since the Galactic center is ~1000 times brighter than fluctuations in the Cosmic Microwave Background (CMB), CMB experiments must carefully account for stray Galactic pickup. We present the level of contamination due to sidelobes for the year one CMB maps produced by the WMAP observatory. For each radiometer, full 4 pi sr antenna gain patterns are determined from a combination of numerical prediction, ground-based and space-based measurements. These patterns are convolved with the WMAP year one sky maps and observatory scan pattern to generate expected sidelobe signal contamination, for both intensity and polarized microwave sky maps. Outside of the Galactic plane, we find rms values for the expected sidelobe pickup of 15, 2.1, 2.0, 0.3, 0.5 uK for K, Ka, Q, V, and W bands respectively. Except at K band, the rms polarized contamination is <<1 uK. Angular power spectra of the Galactic pickup are presented.

We present the angular power spectrum derived from the first-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We study a variety of power spectrum estimation methods and data combinations and demonstrate that the results are robust. The data are modestly contaminated by diffuse Galactic foreground emission, but we show that a simple Galactic template model is sufficient to remove the signal. Point sources produce a modest contamination in the low frequency data. After masking ~700 known bright sources from the maps, we estimate residual sources contribute ~3500 uK^2 at 41 GHz, and ~130 uK^2 at 94 GHz, to the power spectrum l*(l+1)*C_l/(2*pi) at l=1000. Systematic errors are negligible compared to the (modest) level of foreground emission. Our best estimate of the power spectrum is derived from 28 cross-power spectra of statistically independent channels. The final spectrum is essentially independent of the noise properties of an individual radiometer. The resulting spectrum provides a definitive measurement of the CMB power spectrum, with uncertainties limited by cosmic variance, up to l~350. The spectrum clearly exhibits a first acoustic peak at l=220 and a second acoustic peak at l~540 and it provides strong support for adiabatic initial conditions. Kogut et al. (2003) analyze the C_l^TE power spectrum, and present evidence for a relatively high optical depth, and an early period of cosmic reionization. Among other things, this implies that the temperature power spectrum has been suppressed by \~30% on degree angular scales, due to secondary scattering.

We describe the calibration and data processing methods used to generate full-sky maps of the cosmic microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy Probe (WMAP) observations. Detailed limits on residual systematic errors are assigned based largely on analyses of the flight data supplemented, where necessary, with results from ground tests. The data are calibrated in flight using the dipole modulation of the CMB due to the observatory's motion around the Sun. This constitutes a full-beam calibration source. An iterative algorithm simultaneously fits the time-ordered data to obtain calibration parameters and pixelized sky map temperatures. The noise properties are determined by analyzing the time-ordered data with this sky signal estimate subtracted. Based on this, we apply a pre-whitening filter to the time-ordered data to remove a low level of 1/f noise. We infer and correct for a small ~1% transmission imbalance between the two sky inputs to each differential radiometer, and we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. No other systematic error corrections are applied to the data. Calibration and baseline artifacts, including the response to environmental perturbations, are negligible. Systematic uncertainties are comparable to statistical uncertainties in the characterization of the beam response. Both are accounted for in the covariance matrix of the window function and are propagated to uncertainties in the final power spectrum. We characterize the combined upper limits to residual systematic uncertainties through the pixel covariance matrix.

We present the design, manufacturing methods, and characterization of 20 microwave feed horns currently in use on the Microwave Anisotropy Probe (MAP) satellite. The nature of the cosmic microwave background (CMB) anisotropy requires a detailed understanding of the properties of every optical component of a microwave telescope. In particular, the properties of the feeds must be known so that the forward gain and sidelobe response of the telescope can be modeled and so that potential systematic effects may be computed. MAP requires low emissivity, azimuthally symmetric, low-sidelobe feeds in five microwave bands (K, Ka, Q, V, and W) that fit within a constrained geometry. The beam pattern of each feed is modeled and compared with measurements; the agreement is generally excellent to the -60 dB level (80 degrees from the beam peak). This agreement verifies the beam-predicting software and the manufacturing process. The feeds also affect the properties and modeling of the microwave receivers. To this end, we show that the reflection from the feeds is less than -25 dB over most of each band and that their emissivity is acceptable. The feeds meet their multiple requirements.

The primary goal of the MAP satellite, now in orbit, is to make high fidelity polarization sensitive maps of the full sky in five frequency bands between 20 and 100 GHz. From these maps we will characterize the properties of the cosmic microwave background (CMB) anisotropy and Galactic and extragalactic emission on angular scales ranging from the effective beam size, <0.23 degree, to the full sky. MAP is a differential microwave radiometer. Two back-to-back shaped offset Gregorian telescopes feed two mirror symmetric arrays of ten corrugated feeds. We describe the prelaunch design and characterization of the optical system, compare the optical models to the measurements, and consider multiple possible sources of systematic error.

The purpose of the MAP mission is to determine the geometry, content, and evolution of the universe via a 13 arcmin full-width-half-max (FWHM) resolution full sky map of the temperature anisotropy of the cosmic microwave background radiation with uncorrelated pixel noise, minimal systematic errors, multifrequency observations, and accurate calibration. These attributes were key factors in the success of NASA's Cosmic Background Explorer (COBE) mission, which made a 7 degree FWHM resolution full sky map, discovered temperature anisotropy, and characterized the fluctuations with two parameters, a power spectral index and a primordial amplitude. Following COBE considerable progress has been made in higher resolution measurements of the temperature anisotropy. With 45 times the sensitivity and 33 times the angular resolution of the COBE mission, MAP will vastly extend our knowledge of cosmology. MAP will measure the physics of the photon-baryon fluid at recombination. From this, MAP measurements will constrain models of structure formation, the geometry of the universe, and inflation. In this paper we present a pre-launch overview of the design and characteristics of the MAP mission. This information will be necessary for a full understanding of the MAP data and results, and will also be of interest to scientists involved in the design of future cosmic microwave background experiments and/or space science missions.

The spherical Mexican Hat wavelet is introduced in this paper, with the aim of testing the Gaussianity of the Cosmic Microwave Background temperature fluctuations. Using the information given by the wavelet coefficients at several scales, we have performed several statistical tests on the COBE-DMR maps to search for evidence of non-Gaussianity. Skewness, kurtosis, scale-scale correlations (for two and three scales) and Kolmogorov-Smirnov tests indicate that the COBE-DMR data are consistent with a Gaussian distribution. We have extended the analysis to compare temperature values provided by COBE-DMR data with distributions (obtained from Gaussian simulations) at each pixel and at each scale. The number of pixels with temperature values outside the 95% and the 99% is consistent with that obtained from Gaussian simulations, at all scales. Moreover, the extrema values for COBE-DMR data (maximum and minimum temperatures in the map) are also consistent with those obtained from Gaussian simulations.

We present an overview of the upcoming Microwave Anisotropy Probe (MAP) mission, with an emphasis on those aspects of the mission that simplify the data analysis. The method used to make sky maps from the differential temperature data is reviewed and we present some of the noise properties expected from these maps. An overview of the method we plan to use to mine the angular power spectrum from the mega-pixel sky maps closes the paper.

We use a simple model to investigate the effect of polarized Galactic foreground emission on the ability of planned CMB missions to detect and model CMB polarization. Emission from likely polarized sources (synchrotron and spinning dust) would dominate the polarization of the microwave sky at frequencies below 90 GHz if known Galactic foregrounds are at least 10% polarized at high latitude (|b| > 30 deg). Maps of polarization at frequencies below 90 GHz will likely require correction for foreground emission to enable statistical analysis of the individual Stokes Q or U components. The temperature-polarization cross-correlation is less affected by foreground emission, even if the foreground polarization is highly correlated with the foreground intensity. Polarized foregrounds, even if uncorrected, do not dominate the uncertainty in the temperature-polarization cross-correlation for instrument noise levels typical of the MAP experiment. Methods which remove galactic signals at the cost of signal to noise ratio should carefully balance the value of rejecting faint foregrounds with the cost of increased instrument noise.

We present an analysis of the Gaussianity of the 4-year COBE-DMR data (in HEALPix pixelisation) based on spherical wavelets. The skewness, kurtosis and scale-scale correlation spectra are computed from the detail wavelet coefficients at each scale. The sensitivity of the method to the orientation of the data is also taken into account. We find a single detection of non-Gaussianity at the $>99%$ confidence level in one of our statistics. We use Monte-Carlo simulations to assess the statistical significance of this detection and find that the probability of obtaining such a detection by chance for an underlying Gaussian field is as high as 0.69. Therefore, our analysis does not show evidence of non-Gaussianity in the COBE-DMR data.

We investigate the Gaussianity of the 4-year COBE-DMR data (in HEALPix pixelisation) using an analysis based on spherical Haar wavelets. We use all the pixels lying outside the Galactic cut and compute the skewness, kurtosis and scale-scale correlation spectra for the wavelet coefficients at each scale. We also take into account the sensitivity of the method to the orientation of the input signal. We find a detection of non-Gaussianity at $> 99$ per cent level in just one of our statistics. Taking into account the total number of statistics computed, we estimate that the probability of obtaining such a detection by chance for an underlying Gaussian field is 0.69. Therefore, we conclude that the spherical wavelet technique shows no strong evidence of non-Gaussianity in the COBE-DMR data.

There is enormous potential to advance cosmology from statistical characterizations of cosmic microwave background sky maps. The angular power spectrum of the microwave anisotropy is a particularly important statistic. Existing algorithms for computing the angular power spectrum of a pixelized map typically require O(N^3) operations and O(N^2) storage, where N is the number of independent pixels in the map. The MAP and Planck satellites will produce megapixel maps of the cosmic microwave background temperature at multiple frequencies; thus, existing algorithms are not computationally feasible. In this article, we introduce an algorithm that requires O(N^2) operations and O(N^1.5) storage that can find the minimum variance power spectrum from sky map data roughly one million times faster than was previously possible. This makes feasible an analysis that was hitherto intractable.

Inflationary cosmology predicts that the particle horizon should be generically much bigger than the present-day Hubble radius, $1/H_0$. This implies a special regime of super-Hubble scale energy density fluctuations imprinted on the cosmic microwave background radiation (CMBR), which from present theory could only be explained by inflation Causality constraints are used to determine models for the power spectrum that accommodate a suppression scale. A three parameter likelihood analysis is performed of the COBE-DMR 4-year data with respect to the amplitude, spectral index, and suppression scale. It is found that all suppression length scales larger than $1/H_0$ are consistent with the data, but that scales of order $1/H_0$ are slightly preferred, at roughly the one-sigma level. Many non-inflation models would be consistent with a small suppression length scale, whereas for standard inflation models, the duration of the inflation epoch would have to be bounded by a fairly small upper limit. Suppression scales smaller than $1/H_0$ are strongly excluded by the anisotrophy data.

We analyze the Cosmic Microwave Background (CMB) anisotropy data from the independent COBE FIRAS and DMR observations. We extract the frequency spectrum of the FIRAS signal that has the spatial distribution seen by DMR and show that it is consistent with CMB temperature fluctuations in the radiation well into the Wien region of the spectrum. Conversely, we form a map of the Planckian component of the sky temperature from FIRAS and show that it correlates with the DMR anisotropy map. The rms fluctuations at angular scales of 7 degrees are 48 \pm 14 uK for the FIRAS data compared with 35 \pm 2 uK for the DMR data and 31 \pm 6 uK for the correlated combination (1 sigma uncertainties). The consistency of these data, from very different instruments with very different observing strategies, provide compelling support for the interpretation that the signal seen by DMR is, in fact, temperature anisotropy of cosmological origin. The data also limit rms fluctuations in the Compton y parameter, observable via the Sunyaev- Zel'dovich effect, to Delta_y < 3 x 10^{-6} (95% CL) on 7 degree angular scales.

Small departures from a homogeneous isotropic spacetime create observable features in the large-scale anisotropy of the cosmic microwave background. We cross-correlate the maps of the cosmic microwave background anisotropy from the Cosmic Background Explorer (COBE) Differential Microwave Radiometers (DMR) 4-year data set with template maps from Bianchi VII_h cosmological models to limit global rotation or shear in the early universe. On the largest scales, spacetime is well described by the Friedmann-Robertson-Walker metric, with departures from isotropy about each spatial point limited to shear sigma/H_0 < 10^{-9} and rotation omega/H_0 < 6 x 10^{-8} for 0.1 < Omega_0 < 1.

Recent detections of cosmic microwave background (CMB) anisotropy at half-degree angular scales show considerable scatter in the reported amplitude even at similar angular resolution. We use Monte Carlo techniques to simulate the current set of medium-scale CMB observations, including all relevant aspects of sky coverage and measurement technique. The scatter in the reported amplitudes is well within the range expected for the standard cold dark matter (CDM) cosmological model, and results primarily from the restricted sky coverage of each experiment. Within the context of standard CDM current observations of CMB anisotropy support the detection of a ``Doppler peak'' in the CMB power spectrum consistent with baryon density 0.01 < Omega_b < 0.13 (95% confidence) for Hubble constant H_0 = 50 km/s/Mpc. The uncertainties are approximately evenly divided between instrument noise and cosmic variance arising from the limited sky coverage.

We use the COBE Differential Microwave Radiometers (DMR) 4-year sky maps to model Galactic microwave emission at high latitudes (|b| > 20 deg). Cross-correlation of the DMR maps with Galactic template maps detects fluctuations in the high-latitude microwave sky brightness with the angular variation of the DIRBE far-infrared dust maps and a frequency dependence consistent with a superposition of dust and free-free emission. We find no significant correlations between the DMR maps and various synchrotron templates. On the largest angular scales (e.g., quadrupole), Galactic emission is comparable in amplitude to the anisotropy in the cosmic microwave background (CMB). The CMB quadrupole amplitude, after correction for Galactic emission, has amplitude $Q_{rms}$ = 10.7 uK with random uncertainty 3.6 uK and systematic uncertainty 7.1 uK from uncertainty in our knowledge of Galactic microwave emission.

We employ a pixel-based likelihood technique to estimate the angular power spectrum of the COBE Differential Microwave Radiometer (DMR) 4-year sky maps. The spectrum is consistent with a scale-invariant power-law form with a normalization, expressed in terms of the expected quadrupole anisotropy, of Q_{rms-PS|n=1} = 18 +/- 1.4 uK, and a best-fit spectral index of 1.2 +/- 0.3. The normalization is somewhat smaller than we concluded from the 2-year data, mainly due to additional Galactic modeling. We extend the analysis to investigate the extent to which the ``small" quadrupole observed in our sky is statistically consistent with a power-law spectrum. The most likely quadrupole amplitude is somewhat dependent on the details of Galactic foreground subtraction and data selection, ranging between 7 and 10 uK, but in no case is there compelling evidence that the quadrupole is too small to be consistent with a power-law spectrum. We conclude with a likelihood analysis of the band power amplitude in each of four spectral bands between l = 2 and 40, and find no evidence for deviations from a simple power-law spectrum.

The angular power spectrum estimator developed by Peebles (1973) and Hauser & Peebles (1973) has been modified and applied to the 4 year maps produced by the COBE DMR. The power spectrum of the observed sky has been compared to the power spectra of a large number of simulated random skies produced with noise equal to the observed noise and primordial density fluctuation power spectra of power law form, with $P(k) \propto k^n$. The best fitting value of the spectral index in the range of spatial scales corresponding to spherical harmonic indices $3 \leq \ell \lesssim 30$ is an apparent spectral index $n_{app}$ = 1.13 (+0.3) (-0.4) which is consistent with the Harrison-Zel'dovich primordial spectral index $n_{pri} = 1$ The best fitting amplitude for $n_{app} = 1$ is $\langle Q_{RMS}^2\rangle^{0.5}$ = 18 uK.

The 2-point temperature correlation function is evaluated from the 4-year COBE DMR microwave anisotropy maps. We examine the 2-point function, which is the Legendre transform of the angular power spectrum, and show that the data are statistically consistent from channel to channel and frequency to frequency. The most likely quadrupole normalization is computed for a scale-invariant power-law spectrum of CMB anisotropy, using a variety of data combinations. For a given data set, the normalization inferred from the 2-point data is consistent with that inferred by other methods. The smallest and largest normalization deduced from any data combination are 16.4 and 19.6 uK respectively, with a value ~18 uK generally preferred.

We search the high-latitude portion of the COBE Differential Microwave Radiometers (DMR) 4-year sky maps for evidence of a non-Gaussian temperature distribution in the cosmic microwave background. The genus, 3-point correlation function, and 2-point correlation function of temperature maxima and minima are all in excellent agreement with the hypothesis that the CMB anisotropy on angular scales of 7 degrees or larger represents a random-phase Gaussian field. A likelihood comparison of the DMR sky maps to a set of random-phase non-Gaussian toy models selects the exact Gaussian model as most likely. Monte Carlo simulations show that the 2-point correlation of the peaks and valleys in the maps provides the greatest discrimination among the class of models tested.

Fourier analysis and power spectrum estimation of the cosmic microwave background anisotropy on an incompletely sampled sky developed by Gorski (1994) has been applied to the high-latitude portion of the 4-year COBE DMR 31.5, 53 and 90 GHz sky maps. Likelihood analysis using newly constructed Galaxy cuts (extended beyond |b| = 20deg to excise the known foreground emission) and simultaneously correcting for the faint high latitude galactic foreground emission is conducted on the DMR sky maps pixelized in both ecliptic and galactic coordinates. The Bayesian power spectrum estimation from the foreground corrected 4-year COBE DMR data renders n ~ 1.2 +/- 0.3, and Q_{rms-PS} ~ 15.3^{+3.7}_{-2.8} microK (projections of the two-parameter likelihood). These results are consistent with the Harrison-Zel'dovich n=1 model of amplitude Q_{rms-PS} ~ 18 microK detected with significance exceeding 14sigma (dQ/Q < 0.07). (A small power spectrum amplitude drop below the published 2-year results is predominantly due to the application of the new, extended Galaxy cuts.)

We limit the possible contributions from non-cosmological sources to the COBE-DMR four-year sky maps. The DMR data are cross-correlated with maps of rich clusters, extragalactic IRAS sources, HEAO-1 A-2 X-ray emission and 5 GHz radio sources using a Fourier space technique. There is no evidence of significant contamination by such sources at an rms level of ~8 uK (95% confidence level at 7 degree resolution) in the most sensitive 53 GHz sky map. This level is consistent with previous limits set by analysis of earlier DMR data and by simple extrapolations from existing source models. We place a 95% confidence level rms limit on the Comptonization parameter averaged over the high-latitude sky of $\delta y\ <\ 1\ \times\ 10^{-6}$. Extragalactic sources have an insignificant effect on CMB power spectrum parameterizations determined from the DMR data.

The sky-RMS is the simplest model-independent characterization of a cosmological anisotropy signal. The RMS temperature fluctuations determined from the $COBE$-DMR four-year sky maps are frequency independent, consistent with the hypothesis that they are cosmological in origin, with a typical amplitude at 7 degrees of ~ 35 \pm\ 2 uK and at 10 degrees of ~29 \pm\ 1 uK. A joint analysis of the 7 deg and 10 deg "cross"-RMS derived from the data in both Galactic and Ecliptic coordinates is used to determine the RMS quadrupole normalization, $Q_{rms-PS}$, for a scale-invariant Harrison-Zel'dovich power law model. The difference in the inferred normalizations either including or excluding the quadrupole is a consequence of the sensitivity of the method to the Galaxy contaminated observed sky quadrupole. Whilst there are variations depending on the data selection, all results are consistent with an inferred \qrms\ normalization of ~18 \pm\ 2 uK.

The Differential Microwave Radiometers (DMR) instrument aboard the Cosmic Background Explorer (COBE) has mapped the full microwave sky to mean sensitivity 26 uK per 7 deg field of view. The absolute calibration is determined to 0.7% with drifts smaller than 0.2% per year. We have analyzed both the raw differential data and the pixelized sky maps for evidence of contaminating sources such as solar system foregrounds, instrumental susceptibilities, and artifacts from data recovery and processing. Most systematic effects couple only weakly to the sky maps. The largest uncertainties in the maps result from the instrument susceptibility to the Earth's magnetic field, microwave emission from the Earth, and upper limits to potential effects at the spacecraft spin period. Systematic effects in the maps are small compared to either the noise or the celestial signal: the 95% confidence upper limit for the pixel-pixel rms from all identified systematics is less than 6 uK in the worst channel. A power spectrum analysis of the (A-B)/2 difference maps shows no evidence for additional undetected systematic effects.

The cosmic microwave background radiation provides unique constraints on cosmological models. In this Letter we present a summary of the spatial properties of the cosmic microwave background radiation based on the full 4 years of COBE DMR observations, as detailed in a set of companion Letters. The anisotropy is consistent with a scale-invariant power law model and Gaussian statistics. With full use of the multi-frequency 4-year DMR data, including our estimate of the effects of Galactic emission, we find a power-law spectral index of $n=1.2\pm 0.3$ and a quadrupole normalization $Q_{rms-PS}=15.3^{+3.8}_{-2.8}$ $\mu$K. For $n=1$ the best-fit normalization is $Q_{rms-PS}\vert_{n=1}=18\pm 1.6$ $\mu$K. These values are consistent with both our previous 1-year and 2-year results. The results include use of the $\ell=2$ quadrupole term; exclusion of this term gives consistent results, but with larger uncertainties. The 4-year sky maps, presented in this Letter, portray an accurate overall visual impression of the anisotropy since the signal-to-noise ratio is ~2 per 10 degree sky map patch. The improved signal-to-noise ratio of the 4-year maps also allows for improvements in Galactic modeling and limits on non-Gaussian statistics.

A major goal of cosmology is to obtain sensitive, high resolution maps of the Cosmic Microwave Background (CMB) anisotropy. Such maps, as would be produced by the recently proposed Microwave Anisotropy Probe (MAP), will contain a wealth of primary information about conditions in the early universe. To mitigate systematic effects when observing the microwave background, it is desirable for the raw data to be collected in differential form: as a set of temperature differences between points in the sky. However, the production of large (mega-pixel) maps from a set of temperature differences is a potentially severe computational challenge. We present a new technique for producing maps from differential radiometer data that has a computational cost that grows in the slowest possible way with increasing angular resolution and number of map pixels. The required central processor (CPU) time is proportional to the number of differential data points and the required random access memory (RAM) is proportional to the number of map pixels. We test our technique, and demonstrate its feasibility, by simulating one year of a space-borne anisotropy mission.

We cross-correlate the COBE-DMR 2-year sky maps with spatial templates from long-wavelength radio surveys and the far-infrared COBE DIRBE maps. We place an upper limit on the spectral index of synchrotron radiation beta_{synch} < -2.9 between 408 MHz and 31.5 GHz. We obtain a statistically significant cross-correlation with the DIRBE maps whose dependence on the DMR frequencies indicates a superposition of dust and free-free emission. The high-latitude dust emission (|b| > 30 deg) is well fitted by a single dust component with temperature T = 18 (+3, -7) K and emissivity epsilon ~ nu^beta with beta = 1.9 (+3.0, -0.5). The free-free emission is spatially correlated with the dust on angular scales larger than the 7 degree DMR beam, with rms variations 5.3 +/- 1.8 uK at 53 GHz and angular power spectrum P ~ ell^{-3}. If this correlation persists to smaller angular scales, free-free emission should not be a significant contaminant to measurements of the cosmic microwave anisotropy at degree angular scales for frequencies above 20 GHz.

We compute the three-point temperature correlation function of the COBE Differential Microwave Radiometer (DMR) two-year sky maps to search for evidence of non-Gaussian temperature fluctuations. We detect three-point correlations in our sky with a substantially higher signal-to-noise ratio than from the first year data. However, the magnitude of the signal is consistent with the level of cosmic variance expected from Gaussian fluctuations, even when the low order multipole moments, up to l = 9, are filtered from the data. These results do not strongly constrain most existing models of structure formation, but the absence of intrinsic three-point correlations on large angular scales is an important consistency test for such models.

The frequency-independent RMS temperature fluctuations determined from the COBE-DMR two year sky maps are used to infer the parameter Q_{rms-PS}, which characterizes the normalization of power law models of primordial cosmological temperature anisotropy. In particular, a 'cross'-RMS statistic is used to determine Q_{rms-PS} for a forced fit to a scale-invariant Harrison-Zel'dovich (n = 1) spectral model. Using a joint analysis of the 7 degree and 10 degree RMS temperature derived from both the 53 and 90 GHz sky maps, we find Q_{rms-PS} = 17.0^{+2.5}_{-2.1} uK when the low quadrupole is included, and Q_{rms-PS} = 19.4^{+2.3}_{-2.1} uK excluding the quadrupole. These results are consistent with the n = 1 fits from more sensitive methods (e.g. power spectrum, correlation function). The effect of the low quadrupole derived from the COBE-DMR data on the inferred Q_{rms-PS} normalization is investigated. A bias to lower Q_{rms-PS} is found when the quadrupole is included. The higher normalization for a forced n = 1 fit is then favored by the cross-RMS technique. As initially pointed out in Wright et al. (1994a) and further discussed here, analytic formulae for the RMS sky temperature fluctuations will NOT provide the correct normalization amplitude.

We use the two-point correlation function of the extrema points (peaks and valleys) in the COBE Differential Microwave Radiometers (DMR) 2-year sky maps as a test for non-Gaussian temperature distribution in the cosmic microwave background anisotropy. A maximum likelihood analysis compares the DMR data to n=1 toy models whose random-phase spherical harmonic components $a_{\ell m}$ are drawn from either Gaussian, $\chi^2$, or log-normal parent populations. The likelihood of the 53 GHz (A+B)/2 data is greatest for the exact Gaussian model. All non-Gaussian models tested are ruled out at 90\% confidence, limited by type II errors in the statistical inference. The extrema correlation function is a stronger test for this class of non-Gaussian models than topological statistics such as the genus.

High-contrast peaks in the cosmic microwave background (CMB) anisotropy can appear as unresolved sources to observers. We fit simulated CMB maps generated with a cold dark matter model to a set of unresolved features at instrumental resolution 0.5 degrees to 1.5 degrees and present the integral density per steradian n(>|T|) of unresolved features brighter than threshold temperature |T|. A typical medium-scale experiment observing 0.001 sr at 0.5 degree resolution would expect to observe one unresolved feature brighter than 85 \muK after convolution with the beam profile, with less than 5% probability to observe a source brighter than 150 \muK. Increasing the power-law index of primordial density perturbations n from 1 to 1.5 raises these temperature limits |T| by a factor of 2. The results are only weakly dependent on the assumed values of the Hubble constant and baryon density.

We have simulated full-sky maps of the Cosmic Microwave Background (CMB) anisotropy expected from Cold Dark Matter (CDM) models at 0.5 and 1.0 degree angular resolution. Statistical properties of the maps are presented as a function of sky coverage, angular resolution, and instrument noise, and the implications of these results for observability of the Doppler peak are discussed. The rms fluctuations in a map are not a particularly robust probe of the existence of a Doppler peak, however, a full correlation analysis can provide reasonable sensitivity. We find that sensitivity to the Doppler peak depends primarily on the fraction of sky covered, and only secondarily on the angular resolution and noise level. Color plates and one-dimensional scans of the maps are presented to visually illustrate the anisotropies.

The first two years of COBE DMR observations of the CMB anisotropy are analyzed and compared with our previously published first year results. The results are consistent, but the addition of the second year of data increases the precision and accuracy of the detected CMB temperature fluctuations. The two-year 53 GHz data are characterized by RMS temperature fluctuations of DT=44+/-7 uK at 7 degrees and DT=30.5+/-2.7 uK at 10 degrees angular resolution. The 53X90 GHz cross-correlation amplitude at zero lag is C(0)^{1/2}=36+/-5 uK (68%CL) for the unsmoothed 7 degree DMR data. A likelihood analysis of the cross correlation function, including the quadrupole anisotropy, gives a most likely quadrupole-normalized amplitude Q_{rms-PS}=12.4^{+5.2}_{-3.3} uK (68% CL) and a spectral index n=1.59^{+0.49}_{-0.55} for a power law model of initial density fluctuations, P(k)~k^n. With n fixed to 1.0 the most likely amplitude is 17.4 +/-1.5 uK (68% CL). Excluding the quadrupole anisotropy we find Q_{rms-PS}= 16.0^{+7.5}_{-5.2} uK (68% CL), n=1.21^{+0.60}_{-0.55}, and, with n fixed to 1.0 the most likely amplitude is 18.2+/-1.6 uK (68% CL). Monte Carlo simulations indicate that these derived estimates of n may be biased by ~+0.3 (with the observed low value of the quadrupole included in the analysis) and {}~+0.1 (with the quadrupole excluded). Thus the most likely bias-corrected estimate of n is between 1.1 and 1.3. Our best estimate of the dipole from the two-year DMR data is 3.363+/-0.024 mK towards Galactic coordinates (l,b)= (264.4+/-0.2 degrees, +48.1+/-0.4 degrees), and our best estimate of the RMS quadrupole amplitude in our sky is 6+/-3 uK.

A new technique of Fourier analysis on a cut sky (Gorski, 1994) has been applied to the two year COBE DMR sky maps. The Bayesian power spectrum estimation results are consistent with the Harrison-Zel'dovich n=1 model. The maximum likelihood estimates of the parameters of the power spectrum of primordial perturbations are n=1.22 (1.02) and Q_{rms-PS}=17 (20) uK including (excluding) the quadrupole anisotropy. The marginal likelihood function on n renders n=1.10 \pm 0.32 (0.87 \pm 0.36).

The {\it COBE} DMR sky maps contain low-level correlated noise. We obtain estimates of the amplitude and pattern of the correlated noise from three techniques: angular averages of the covariance matrix, Monte Carlo simulations of two-point correlation functions, and direct analysis of the DMR maps. The results from the three methods are mutually consistent. The noise covariance matrix of a DMR sky map is diagonal to an accuracy of better than 1\%. For a given sky pixel, the dominant noise covariance occurs with the ring of pixels at an angular separation of $60 \deg$ due to the $60 \deg$ separation of the DMR horns. The mean covariance at $60 \deg$ is $0.45\% ^{+0.18}_{-0.14}$ of the mean variance. Additionally, the variance in a given pixel is $0.7\%$ greater than would be expected from a single beam experiment with the same noise properties. Auto-correlation functions suffer from a $\sim 1.5\; \sigma$ positive bias at $60 \deg$ while cross-correlations have no bias. Published {\it COBE} DMR results are not significantly affected by correlated noise. COBE pre-print 94-

We have searched the temperature maps from the COBE Differential Microwave Radiometers (DMR) first two years of data for evidence of unresolved sources. The high-latitude sky (|b| > 30\deg) contains no sources brighter than 192 uK thermodynamic temperature (322 Jy at 53 GHz). The cumulative count of sources brighter than threshold T, N(> T), is consistent with a superposition of instrument noise plus a scale-invariant spectrum of cosmic temperature fluctuations normalized to Qrms-PS = 17 uK. We examine the temperature maps toward nearby clusters and find no evidence for any Sunyaev-Zel'dovich effect, \Delta y < 7.3 x 10^{-6} (95% CL) averaged over the DMR beam. We examine the temperature maps near the brightest expected radio sources and detect no evidence of significant emission. The lack of bright unresolved sources in the DMR maps, taken with anisotropy measurements on smaller angular scales, places a weak constraint on the integral number density of any unresolved Planck-spectrum sources brighter than flux density S, n(> S) < 2 x 10^4 (S/1 Jy)^{-2} sr^{-1}.

We use statistical and topological quantities to test the COBE-DMR first year sky maps against the hypothesis that the observed temperature fluctuations reflect Gaussian initial density perturbations with random phases. Recent papers discuss specific quantities as discriminators between Gaussian and non-Gaussian behavior, but the treatment of instrumental noise on the data is largely ignored. The presence of noise in the data biases many statistical quantities in a manner dependent on both the noise properties and the unknown CMB temperature field. Appropriate weighting schemes can minimize this effect, but it cannot be completely eliminated. Analytic expressions are presented for these biases, and Monte Carlo simulations used to assess the best strategy for determining cosmologically interesting information from noisy data. The genus is a robust discriminator that can be used to estimate the power law quadrupole-normalized amplitude independently of the 2-point correlation function. The genus of the DMR data are consistent with Gaussian initial fluctuations with Q_rms = 15.7 +/- 2.2 - (6.6 +/- 0.3)(n - 1) uK where n is the power law index. Fitting the rms temperature variations at various smoothing angles gives Q_rms = 13.2 +/- 2.5 uK and n = 1.7 +0.3 -0.6. While consistent with Gaussian fluctuations, the first year data are only sufficient to rule out strongly non-Gaussian distributions of fluctuations.

The FIRAS instrument on the COBE satellite has conducted an unbiased survey of the far-infrared emission from our Galaxy. The first results of this survey were reported by Wright et al. (1991). We report the results of new analyses of this spectral survey, which includes emission lines from 158 um C+, 122 um and 205 um N+, 370 um and 609 um C, and CO J=2-1 through 5-4. We report the morphological distribution along the galactic plane (b=0) of the spectral line emission, and the high galactic latitude intensities of the C+ and 205 um N+ emission. The high galactic latitude intensity cosecant of the 158 um fine structure transition from C+ is presented, and C+ is seen to decrease more rapidly than the far infrared intensity with increasing galactic latitude. C+ and H I emission are closely correlated with a C+ cooling rate of (2.65 +/- 0.15)x10^{-26} erg/s/H-atom. We conclude that this emission arises almost entirely from the Cold Neutral Medium. The high galactic latitude intensity of the 205 um fine structure transition from N+ arises entirely from the Warm Ionized Medium, and its cosecant is presented. We estimate the total ionizing photon rate in the Galaxy to be 3.5x10^{53} ionizing photons per second, based on the 205 um N+ transition.

We compute the three-point temperature correlation function of the {\it COBE} Differential Microwave Radiometer (DMR) first-year sky maps to search for non-Gaussian temperature fluctuations. The level of fluctuations seen in the computed correlation function are too large to be attributable solely to instrument noise. However the fluctuations are consistent with the level expected to result from a superposition of instrument noise and sky signal arising from a Gaussian power law model of initial fluctuations, with a quadrupole normalized amplitude of 17 $\mu$K and a power law spectral index $n = 1$. We place limits on the amplitude of intrinsic three-point correlations with a variety of predicted functional forms.