We present the design and performance of broadband and tunable infrared-blocking filters for millimeter and sub-millimeter astronomy composed of small scattering particles embedded in an aerogel substrate. The ultra-low-density (typically < 150 mg/cm^3) aerogel substrate provides an index of refraction as low as 1.05, removing the need for anti-reflection coatings and allowing for broadband operation from DC to above 1 THz. The size distribution of the scattering particles can be tuned to provide a variable cutoff frequency. Aerogel filters with embedded high-resistivity silicon powder are being produced at 40-cm diameter to enable large-aperture cryogenic receivers for cosmic microwave background polarimeters, which require large arrays of sub-Kelvin detectors in their search for the signature of an inflationary gravitational-wave background.

The Wide-Field Infrared Survey Telescope (WFIRST) will answer fundamental questions about the evolution of dark energy over time and expand the catalog of known exoplanets into new regions of parameter space. Using a Hubble-sized mirror and 18 newly developed HgCdTe 4K x 4K photodiode arrays (H4RG-10), WFIRST will measure the positions and shapes of hundreds of millions of galaxies, the light curves of thousands of supernovae, and the microlensing signals of over a thousand exoplanets toward the bulge of the Galaxy. These measurements require unprecedented sensitivity and characterization of the Wide Field Instrument (WFI), particularly its detectors. The WFIRST project undertook an extensive detector development program to create focal plane arrays that meet these science requirements. These prototype detectors have been characterized and their performance demonstrated in a relevant space-like environment (thermal vacuum, vibration, acoustic, and radiation testing), advancing the H4RG-10's technology readiness level (TRL) to TRL-6. We present the performance characteristics of these TRL-6 demonstration devices.

In order to extract cosmological information from observations of the millimeter and submillimeter sky, foreground components must first be removed to produce an estimate of the cosmic microwave background (CMB). We developed a machine learning approach for doing so for full-sky temperature maps of the millimeter and submillimeter sky. We constructed a Bayesian spherical convolutional neural network architecture to produce a model that captures both spectral and morphological aspects of the foregrounds. Additionally, the model outputs a per-pixel error estimate that incorporates both statistical and model uncertainties. The model was then trained using simulations that incorporated knowledge of these foreground components that was available at the time of the launch of the Planck satellite. On simulated maps, the CMB is recovered with a mean absolute difference of $<4\mu$K over the full sky after masking map pixels with a predicted standard error of $>50\mu$K; the angular power spectrum is also accurately recovered. Once validated with the simulations, this model was applied to Planck temperature observations from its 70GHz through 857GHz channels to produce a foreground-cleaned CMB map at a Healpix map resolution of NSIDE=512. Furthermore, we demonstrate the utility of the technique for evaluating how well different simulations match observations, particularly in regard to the modeling of thermal dust.

The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array that observes the cosmic microwave background (CMB) over 75% of the sky from the Atacama Desert, Chile, at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the large angular scale ($1^\circ\lesssim\theta\leqslant 90^\circ$) CMB polarization to constrain the tensor-to-scalar ratio at the $r\sim0.01$ level and the optical depth to last scattering to the sample variance limit. This paper presents the optical characterization of the 40 GHz telescope during its first observation era, from 2016 September to 2018 February. High signal-to-noise observations of the Moon establish the pointing and beam calibration. The telescope boresight pointing variation is $<0.023^\circ$ ($<1.6$% of the beam's full width at half maximum (FWHM)). We estimate beam parameters per detector and in aggregate, as in the CMB survey maps. The aggregate beam has an FWHM of $1.579^\circ\pm.001^\circ$ and a solid angle of $838 \pm 6\ \mu{\rm sr}$, consistent with physical optics simulations. The corresponding beam window function has a sub-percent error per multipole at $\ell < 200$. An extended $90^\circ$ beam map reveals no significant far sidelobes. The observed Moon polarization shows that the instrument polarization angles are consistent with the optical model and that the temperature-to-polarization leakage fraction is $<10^{-4}$ (95% C.L.). We find that the Moon-based results are consistent with measurements of M42, RCW 38, and Tau A from CLASS's CMB survey data. In particular, Tau A measurements establish degree-level precision for instrument polarization angles.

Weak gravitational lensing has emerged as a leading probe of the growth of cosmic structure. However, the shear signal is very small and accurate measurement depends critically on our ability to understand how non-ideal instrumental effects affect astronomical images. WFIRST will fly a focal plane containing 18 Teledyne H4RG-10 near infrared detector arrays, which present different instrument calibration challenges from previous weak lensing observations. Previous work has shown that correlation functions of flat field images are effective tools for disentangling linear and non-linear inter-pixel capacitance (IPC) and the brighter-fatter effect (BFE). Here we present a Fourier-domain treatment of the flat field correlations, which allows us to expand the previous formalism to all orders in IPC, BFE, and classical non-linearity. We show that biases in simulated flat field analyses in Paper I are greatly reduced through the use of this formalism. We then apply this updated formalism to flat field data from three WFIRST flight candidate detectors, and explore the robustness to variations in the analysis. We find that the BFE is present in all three detectors, and that its contribution to the flat field correlations dominates over the non-linear IPC. The magnitude of the BFE is such that the effective area of a pixel is increased by $(3.54\pm0.03)\times 10^{-7}$ for every electron deposited in a neighboring pixel. We compare IPC maps from flat field autocorrelation measurements to those obtained from the single pixel reset method and find a median difference of 0.113%. After further diagnosis of this difference, we ascribe it largely to an additional source of cross-talk, the vertical trailing pixel effect, and recommend further work to develop a model for this effect. These results represent a significant step toward calibration of the non-ideal effects in WFIRST detectors.

The Earth's magnetic field induces Zeeman splitting of the magnetic dipole transitions of molecular oxygen in the atmosphere, which produces polarized emission in the millimeter-wave regime. This polarized emission is primarily circularly polarized and manifests as a foreground with a dipole-shaped sky pattern for polarization-sensitive ground-based cosmic microwave background experiments, such as the Cosmology Large Angular Scale Surveyor (CLASS), which is capable of measuring large angular scale circular polarization. Using atmospheric emission theory and radiative transfer formalisms, we model the expected amplitude and spatial distribution of this signal and evaluate the model for the CLASS observing site in the Atacama Desert of northern Chile. Then, using two years of observations at 32.3 GHz to 43.7 GHz from the CLASS Q-band telescope, we present a detection of this signal and compare the observed signal to that predicted by the model. We recover an angle between magnetic north and true north of $(-5.5 \pm 0.6)^\circ$, which is consistent with the expectation of $-5.9^\circ$ for the CLASS observing site. When comparing dipole sky patterns fit to both simulated and data-derived sky maps, the dipole directions match to within a degree, and the measured amplitudes match to within ${\sim}20\%$.

We explore the fundamental limits to which reionization histories can be constrained using only large-scale cosmic microwave background (CMB) anisotropy measurements. The redshift distribution of the fractional ionization $x_e(z)$ affects the angular distribution of CMB polarization. We project constraints on the reionization history of the universe using low-noise full-sky temperature and E-mode measurements of the CMB. We show that the measured TE power spectrum, $\hat C_\ell^\mathrm{TE}$, has roughly one quarter of the constraining power of $\hat C_\ell^\mathrm{EE}$ on the reionization optical depth $\tau$, and its addition improves the precision on $\tau$ by 20% over using $\hat C_\ell^\mathrm{EE}$ only. We also use a two-step reionization model with an additional high redshift step, parametrized by an early ionization fraction $x_e^\mathrm{min}$, and a late reionization step at $z_\mathrm{re}$. We find that future high signal-to-noise measurements of the multipoles $10\leqslant\ell<20$ are especially important for breaking the degeneracy between $x_e^\mathrm{min}$ and $z_\mathrm{re}$. In addition, we show that the uncertainties on these parameters determined from a map with sensitivity $10\,\mathrm{\mu K\,arcmin}$ are less than 5% larger than the uncertainties in the noiseless case, making this noise level a natural target for future large sky area E-mode measurements.

We report on the development of a polarization-sensitive dichroic (150/220 GHz) detector array for the Cosmology Large Angular Scale Surveyor (CLASS) delivered to the telescope site in June 2019. In concert with existing 40 and 90 GHz telescopes, the 150/220 GHz telescope will make observations of the cosmic microwave background over large angular scales aimed at measuring the primordial B-mode signal, the optical depth to reionization, and other fundamental physics and cosmology. The 150/220 GHz focal plane array consists of three detector modules with 1020 transition edge sensor (TES) bolometers in total. Each dual-polarization pixel on the focal plane contains four bolometers to measure the two linear polarization states at 150 and 220 GHz. Light is coupled through a planar orthomode transducer (OMT) fed by a smooth-walled feedhorn array made from an aluminum-silicon alloy (CE7). In this work, we discuss the design, assembly, and in-lab characterization of the 150/220 GHz detector array. The detectors are photon-noise limited, and we estimate the total array noise-equivalent power (NEP) to be 2.5 and 4 aW$\sqrt{\mathrm{s}}$ for 150 and 220 GHz arrays, respectively.

We report circular polarization measurements from the first two years of observation with the 40 GHz polarimeter of the Cosmology Large Angular Scale Surveyor (CLASS). CLASS is conducting a multi-frequency survey covering 75% of the sky from the Atacama Desert designed to measure the cosmic microwave background (CMB) linear E and B polarization on angular scales $1^\circ \lesssim \theta \leq 90^\circ$, corresponding to a multipole range of $2 \leq \ell \lesssim 200$. The modulation technology enabling measurements of linear polarization at the largest angular scales from the ground, the Variable-delay Polarization Modulator, is uniquely designed to provide explicit sensitivity to circular polarization (Stokes $V$). We present a first detection of circularly polarized atmospheric emission at 40 GHz that is well described by a dipole with an amplitude of $124\pm4\,\mathrm{\mu K}$ when observed at an elevation of $45^\circ$, and discuss its potential impact as a foreground to CMB experiments. Filtering the atmospheric component, CLASS places a 95% C.L. upper limit of $0.4\,\mathrm{\mu K}^2$ to $13.5\,\mathrm{\mu K}^2$ on $\ell(\ell+1)C_\ell^{VV}/(2\pi)$ between $1 \leq \ell \leq 120$, representing a two-orders-of-magnitude improvement over previous limits.

The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne telescope mission to search for inflationary gravitational waves from the early universe. PIPER employs two 32x40 arrays of superconducting transition-edge sensors, which operate at 100 mK. An open bucket dewar of liquid helium maintains the receiver and telescope optics at 1.7 K. We describe the thermal design of the receiver and sub-kelvin cooling with a continuous adiabatic demagnetization refrigerator (CADR). The CADR operates between 70-130 mK and provides ~10 uW cooling power at 100 mK, nearly five times the loading of the two detector assemblies. We describe electronics and software to robustly control the CADR, overall CADR performance in flight-like integrated receiver testing, and practical considerations for implementation in the balloon float environment.

Consistency checks of cosmological data sets are an important tool because they may suggest systematic errors or the type of modifications to \lcdm\ necessary to resolve current tensions. In this work, we derive an analytic method for calculating the level of correlations between model parameters from two correlated cosmological data sets, which complements more computationally expensive simulations. This method is an extension of the Fisher analysis that assumes a Gaussian likelihood and a known data covariance matrix. We apply this method to the SPTpol temperature and polarization CMB spectra (TE and EE). We find weak correlations between \lcdm\ parameters with a 9$\%$ correlation between the TE-only and EE-only constraints on $H_0$ and a 25$\%$ and 32$\%$ correlation for log($A_s$) and $n_s$ respectively. Despite the negative correlations between the TE and EE power spectra, the correlations in the parameters are positive. The TE-EE parameter differences are consistent with zero, with a PTE of 0.53, in contrast to the PTE of 0.017 reported by SPTpol for the consistency of the TE and EE power spectra with \lcdm. Using simulations we find that the results of these two tests are independent and that this difference can arise simply from statistical fluctuations. Ignoring correlations in the TT-TE and TE-EE comparisons biases the $\chi^2$ low, artificially making parameters look more consistent. Therefore, we conclude that these correlations need to be accounted for when performing internal consistency checks of the TT vs TE vs EE power spectra for future CMB analyses.

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 perform forecasts for how baryon acoustic oscillation (BAO) scale and redshift-space distortion (RSD) measurements from future spectroscopic emission line galaxy (ELG) surveys such as Euclid are degraded in the presence of spectral line misidentification. Using analytic calculations verified with mock galaxy catalogs from log-normal simulations we find that constraints are degraded in two ways, even when the interloper power spectrum is modeled correctly in the likelihood. Firstly, there is a loss of signal-to-noise ratio for the power spectrum of the target galaxies, which propagates to all cosmological constraints and increases with contamination fraction, $f_c$. Secondly, degeneracies can open up between $f_c$ and cosmological parameters. In our calculations this typically increases BAO scale uncertainties at the 10-20% level when marginalizing over parameters determining the broadband power spectrum shape. External constraints on $f_c$, or parameters determining the shape of the power spectrum, for example from cosmic microwave background (CMB) measurements, can remove this effect. There is a near-perfect degeneracy between $f_c$ and the power spectrum amplitude for low $f_c$ values, where $f_c$ is not well determined from the contaminated sample alone. This has the potential to strongly degrade RSD constraints. The degeneracy can be broken with an external constraint on $f_c$, for example from cross-correlation with a separate galaxy sample containing the misidentified line, or deeper sub-surveys.

The Cosmology Large Angular Scale Surveyor (CLASS) is mapping the polarization of the Cosmic Microwave Background (CMB) at large angular scales ($2<\ell\lesssim200$) in search of a primordial gravitational wave B-mode signal down to a tensor-to-scalar ratio of $r \approx 0.01$. The same data set will provide a near sample-variance-limited measurement of the optical depth to reionization. Between June 2016 and March 2018, CLASS completed the largest ground-based Q-band CMB survey to date, covering over 31 000~square-degrees (75% of the sky), with an instantaneous array noise-equivalent temperature (NET) sensitivity of $32~\mu \mbox{K}_{cmb}\sqrt{\mbox{s}}$. We demonstrate that the detector optical loading ($1.6~\mbox{pW}$) and noise-equivalent power ($19~\mbox{aW}\sqrt{\mbox{s}}$) match the expected noise model dominated by photon bunching noise. We derive a $13.1\pm0.3~\mbox{K/pW}$ calibration to antenna temperature based on Moon observations, which translates to an optical efficiency of $0.48\pm0.04$ and a $27~\mbox{K}$ system noise temperature. Finally, we report a Tau A flux density of $308\pm11~\mbox{Jy}$ at $38.4\pm0.2~\mbox{GHz}$, consistent with the WMAP Tau A time-dependent spectral flux density model.

Current cosmological tensions motivate investigating extensions to the standard $\Lambda$CDM model. Additional model parameters are typically varied one or two at a time, in a series of separate tests. The purpose of this paper is to highlight that information is lost by not also examining the correlations between these additional parameters, which arise when their effects on model predictions are similar, even if the parameters are not varied simultaneously. We show how these correlations can be quantified with simulations and Markov Chain Monte Carlo (MCMC) methods. As an example, we assume that $\Lambda$CDM is the true underlying model, and calculate the correlations expected between the phenomenological lensing amplitude parameter, $A_L$, the running of the spectral index, $n_{\rm run}$, and the primordial helium mass fraction, $Y_P$, when these parameters are varied one at a time along with the $\Lambda$CDM parameters in fits to the $\textit{Planck}$ 2015 temperature power spectrum. These correlations are not small, ranging from 0.31 ($A_L-n_{\rm run}$) to $-0.93$ ($n_{\rm run}-Y_P$). We find that the values of these three parameters from the $\textit{Planck}$ data are consistent with $\Lambda$CDM expectations within $0.9\sigma$ when the correlations are accounted for. This does not explain the 1.8-2.7$\sigma$ $\textit{Planck}$ preference for $A_L>1$, but provides an additional $\Lambda$CDM consistency test. For example, if $A_L>1$ was a symptom of an underlying systematic error or some real but unknown physical effect that also produced spurious correlations with $n_{\rm run}$ or $Y_P$ our test might have revealed this. We recommend that future cosmological analyses examine correlations between additional model parameters in addition to investigating them separately, one a time.

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

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

We present the design, manufacturing technique, and characterization of a 3D-printed broadband graded index millimeter wave absorber. The absorber is additively manufactured using a fused filament fabrication (FFF) 3D printer out of a carbon-loaded high impact polystyrene (HIPS) filament and is designed using a space-filling curve to optimize manufacturability using said process. The absorber's reflectivity is measured from 63 GHz to 115 GHz and from 140 GHz to 215 GHz and is compared to electromagnetic simulations. The intended application is for terminating stray light in Cosmic Microwave Background (CMB) telescopes, and the absorber has been shown to survive cryogenic thermal cycling.

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 perform a comparison of WMAP 9-year (WMAP9) and Planck 2015 cosmic microwave background (CMB) temperature power spectra across multipoles $30\leq\ell\leq1200$. We generate simulations to estimate the correlation between the two datasets due to cosmic variance from observing the same sky. We find that their spectra are consistent within $1\sigma$. While we do not implement the optimal "$C^{-1}$" estimator on WMAP maps as in the WMAP9 analysis, we demonstrate that the change of pixel weighting only shifts our results at most at the $0.66\sigma$ level. We also show that changing the fiducial power spectrum for simulations only impacts the comparison at around $0.1\sigma$ level. We exclude $\ell<30$ both because WMAP9 data were included in the Planck 2015 $\ell<30$ analysis, and because the cosmic variance uncertainty on these scales is large enough that any remaining systematic difference between the experiments is extremely unlikely to affect cosmological constraints. The consistency shown in our analysis provides high confidence in both the WMAP9 temperature power spectrum and the overlapping multipole region of Planck 2015's, virtually independent of any assumed cosmological model. Our results indicate that cosmological model differences between Planck and WMAP do not arise from measurement differences, but from the high multipoles not measured by WMAP.

The OSIRIS-REx Asteroid Sample Return Mission is the third mission in NASA's New Frontiers Program and is the first U.S. mission to return samples from an asteroid to Earth. The most important decision ahead of the OSIRIS-REx team is the selection of a prime sample-site on the surface of asteroid (101955) Bennu. Mission success hinges on identifying a site that is safe and has regolith that can readily be ingested by the spacecraft's sampling mechanism. To inform this mission-critical decision, the surface of Bennu is mapped using the OSIRIS-REx Camera Suite and the images are used to develop several foundational data products. Acquiring the necessary inputs to these data products requires observational strategies that are defined specifically to overcome the challenges associated with mapping a small irregular body. We present these strategies in the context of assessing candidate sample-sites at Bennu according to a framework of decisions regarding the relative safety, sampleability, and scientific value across the asteroid's surface. To create data products that aid these assessments, we describe the best practices developed by the OSIRIS-REx team for image-based mapping of irregular small bodies. We emphasize the importance of using 3D shape models and the ability to work in body-fixed rectangular coordinates when dealing with planetary surfaces that cannot be uniquely addressed by body-fixed latitude and longitude.

We revisit the degeneracy between the Hubble constant, $H_0$, and matter density, $\Omega_m$, for current cosmic microwave background (CMB) observations within the standard $\Lambda CDM$ model. We show that Planck, Wilkinson Microwave Anisotropy Probe (WMAP), South Pole Telescope (SPT), and Atacama Cosmology Telescope Polarimeter (ACTPol) temperature power spectra produce different values of the exponent $x$ from minimizing the variance of the product $\Omega_mH_0^x$. The distribution of $x$ from the different data sets does not follow the Markov Chain Monte Carlo (MCMC) best-fit values for $H_0$ or $\Omega_m$. Particularly striking is the difference between Planck multipoles $\ell\leq800$ ($x=2.81$), and WMAP ($x = 2.94$), despite very similar best-fit cosmologies. We use a Fisher matrix analysis to show that, in fact, this range in exponent values is exactly as expected in $\Lambda CDM$ given the multipole coverage and power spectrum uncertainties for each experiment. We show that the difference in $x$ from the Planck $\ell \leq 800$ and WMAP data is explained by a turning point in the relationship between $x$ and the maximum effective multipole, at around $\ell=700$. The value of $x$ is determined by several physical effects, and we highlight the significant impact of gravitational lensing for the high-multipole measurements. Despite the spread of $H_0$ values from different CMB experiments, the experiments are consistent with their sampling of the $\Omega_m-H_0$ degeneracy and do not show evidence for the need for new physics or for the presence of significant underestimated systematics according to these tests. The Fisher calculations can be used to predict the $\Omega_m-H_0$ degeneracy of future experiments.

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}}$.

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 describe feedhorn-coupled polarization-sensitive detector arrays that utilize monocrystalline silicon as the dielectric substrate material. Monocrystalline silicon has a low-loss tangent and repeatable dielectric constant, characteristics that are critical for realizing efficient and uniform superconducting microwave circuits. An additional advantage of this material is its low specific heat. In a detector pixel, two Transition-Edge Sensor (TES) bolometers are antenna-coupled to in-band radiation via a symmetric planar orthomode transducer (OMT). Each orthogonal linear polarization is coupled to a separate superconducting microstrip transmission line circuit. On-chip filtering is employed to both reject out-of-band radiation from the upper band edge to the gap frequency of the niobium superconductor, and to flexibly define the bandwidth for each TES to meet the requirements of the application. The microwave circuit is compatible with multi-chroic operation. Metalized silicon platelets are used to define the backshort for the waveguide probes. This micro-machined structure is also used to mitigate the coupling of out-of-band radiation to the microwave circuit. At 40 GHz, the detectors have a measured efficiency of 90%. In this paper, we describe the development of the 90 GHz detector arrays that will be demonstrated using the Cosmology Large Angular Scale Surveyor (CLASS) ground-based telescope.

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

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.

Project Matsu is a collaboration between the Open Commons Consortium and NASA focused on developing open source technology for the cloud-based processing of Earth satellite imagery. A particular focus is the development of applications for detecting fires and floods to help support natural disaster detection and relief. Project Matsu has developed an open source cloud-based infrastructure to process, analyze, and reanalyze large collections of hyperspectral satellite image data using OpenStack, Hadoop, MapReduce, Storm and related technologies. We describe a framework for efficient analysis of large amounts of data called the Matsu "Wheel." The Matsu Wheel is currently used to process incoming hyperspectral satellite data produced daily by NASA's Earth Observing-1 (EO-1) satellite. The framework is designed to be able to support scanning queries using cloud computing applications, such as Hadoop and Accumulo. A scanning query processes all, or most of the data, in a database or data repository. We also describe our preliminary Wheel analytics, including an anomaly detector for rare spectral signatures or thermal anomalies in hyperspectral data and a land cover classifier that can be used for water and flood detection. Each of these analytics can generate visual reports accessible via the web for the public and interested decision makers. The resultant products of the analytics are also made accessible through an Open Geospatial Compliant (OGC)-compliant Web Map Service (WMS) for further distribution. The Matsu Wheel allows many shared data services to be performed together to efficiently use resources for processing hyperspectral satellite image data and other, e.g., large environmental datasets that may be analyzed for many purposes.

Variable-delay Polarization Modulators (VPMs) are currently being implemented in experiments designed to measure the polarization of the cosmic microwave background on large angular scales because of their capability for providing rapid, front-end polarization modulation and control over systematic errors. Despite the advantages provided by the VPM, it is important to identify and mitigate any time-varying effects that leak into the synchronously modulated component of the signal. In this paper, the effect of emission from a $300$ K VPM on the system performance is considered and addressed. Though instrument design can greatly reduce the influence of modulated VPM emission, some residual modulated signal is expected. VPM emission is treated in the presence of rotational misalignments and temperature variation. Simulations of time-ordered data are used to evaluate the effect of these residual errors on the power spectrum. The analysis and modeling in this paper guides experimentalists on the critical aspects of observations using VPMs as front-end modulators. By implementing the characterizations and controls as described, front-end VPM modulation can be very powerful for mitigating $1/f$ noise in large angular scale polarimetric surveys. None of the systematic errors studied fundamentally limit the detection and characterization of B-modes on large scales for a tensor-to-scalar ratio of $r=0.01$. Indeed, $r<0.01$ is achievable with commensurately improved characterizations and controls.

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.

We consider the effectiveness of foreground cleaning in the recovery of Cosmic Microwave Background (CMB) polarization sourced by gravitational waves for tensor-to-scalar ratios in the range $0<r<0.1$. Using the planned survey area, frequency bands, and sensitivity of the Cosmology Large Angular Scale Surveyor (CLASS), we simulate maps of Stokes $Q$ and $U$ parameters at 40, 90, 150, and 220 GHz, including realistic models of the CMB, diffuse Galactic thermal dust and synchrotron foregrounds, and Gaussian white noise. We use linear combinations (LCs) of the simulated multifrequency data to obtain maximum likelihood estimates of $r$, the relative scalar amplitude $s$, and LC coefficients. We find that for 10,000 simulations of a CLASS-like experiment using only measurements of the reionization peak ($\ell\leq23$), there is a 95% C.L. upper limit of $r<0.017$ in the case of no primordial gravitational waves. For simulations with $r=0.01$, we recover at 68% C.L. $r=0.012^{+0.011}_{-0.006}$. The reionization peak corresponds to a fraction of the multipole moments probed by CLASS, and simulations including $30\leq\ell\leq100$ further improve our upper limits to $r<0.008$ at 95% C.L. ($r=0.01^{+0.004}_{-0.004}$ for primordial gravitational waves with $r=0.01$). In addition to decreasing the current upper bound on $r$ by an order of magnitude, these foreground-cleaned low multipole data will achieve a cosmic variance limited measurement of the E-mode polarization's reionization peak.

Characterization of the minute cosmic microwave background polarization signature requires multi-frequency, high-throughput precision instrument systems. We have previously described the detector fabrication of a 40 GHz focal plane and now describe the fabrication of detector modules for measurement of the CMB at 90 GHz. The 90 GHz detectors are a scaled version of the 40 GHz architecture where, due to smaller size detectors, we have implemented a modular (wafer level) rather than the chip-level architecture. The new fabrication process utilizes the same design rules with the added challenge of increased wiring density to the 74 TES's as well as a new wafer level hybridization procedure. The hexagonally shaped modules are tile-able, and as such, can be used to form the large focal planes required for a space-based CMB polarimeter. The detectors described here will be deployed in two focal planes with 7 modules each in the Johns Hopkins University led ground-based Cosmology Large Angular Scale Surveyor (CLASS) telescope.

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.

The Prime Focus Spectrograph (PFS) is an optical/near-infrared multifiber spectrograph with 2394 science fibers distributed across a 1.3-deg diameter field of view at the Subaru 8.2-m telescope. The wide wavelength coverage from 0.38 {\mu}m to 1.26 {\mu}m, with a resolving power of 3000, simultaneously strengthens its ability to target three main survey programs: cosmology, galactic archaeology and galaxy/AGN evolution. A medium resolution mode with a resolving power of 5000 for 0.71 {\mu}m to 0.89 {\mu}m will also be available by simply exchanging dispersers. We highlight some of the technological aspects of the design. To transform the telescope focal ratio, a broad-band coated microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of the cable system, optimizing overall throughput; a fiber with low focal ratio degradation is selected for the fiber-positioner and fiber-slit components, minimizing the effects of fiber movements and fiber bending. Fiber positioning will be performed by a positioner consisting of two stages of piezo-electric rotary motors. The positions of these motors are measured by taking an image of artificially back-illuminated fibers with the metrology camera located in the Cassegrain container; the fibers are placed in the proper location by iteratively measuring and then adjusting the positions of the motors. Target light reaches one of the four identical fast-Schmidt spectrograph modules, each with three arms. The PFS project has passed several project-wide design reviews and is now in the construction phase.

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) 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$.

We report on the status and development of polarization-sensitive detectors for millimeter-wave applications. The detectors are fabricated on single-crystal silicon, which functions as a low-loss dielectric substrate for the microwave circuitry as well as the supporting membrane for the Transition-Edge Sensor (TES) bolometers. The orthomode transducer (OMT) is realized as a symmetric structure and on-chip filters are employed to define the detection bandwidth. A hybridized integrated enclosure reduces the high-frequency THz mode set that can couple to the TES bolometers. An implementation of the detector architecture at Q-band achieves 90% efficiency in each polarization. The design is scalable in both frequency coverage, 30-300 GHz, and in number of detectors with uniform characteristics. Hence, the detectors are desirable for ground-based or space-borne instruments that require large arrays of efficient background-limited cryogenic detectors.

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 Prime Focus Spectrograph (PFS) is an optical/near-infrared multi-fiber spectrograph with 2394 science fibers, which are distributed in 1.3 degree diameter field of view at Subaru 8.2-meter telescope. The simultaneous wide wavelength coverage from 0.38 um to 1.26 um, with the resolving power of 3000, strengthens its ability to target three main survey programs: cosmology, Galactic archaeology, and galaxy/AGN evolution. A medium resolution mode with resolving power of 5000 for 0.71 um to 0.89 um also will be available by simply exchanging dispersers. PFS takes the role for the spectroscopic part of the Subaru Measurement of Images and Redshifts project, while Hyper Suprime-Cam works on the imaging part. To transform the telescope plus WFC focal ratio, a 3-mm thick broad-band coated glass-molded microlens is glued to each fiber tip. A higher transmission fiber is selected for the longest part of cable system, while one with a better FRD performance is selected for the fiber-positioner and fiber-slit components, given the more frequent fiber movements and tightly curved structure. Each Fiber positioner consists of two stages of piezo-electric rotary motors. Its engineering model has been produced and tested. Fiber positioning will be performed iteratively by taking an image of artificially back-illuminated fibers with the Metrology camera located in the Cassegrain container. The camera is carefully designed so that fiber position measurements are unaffected by small amounts of high special-frequency inaccuracies in WFC lens surface shapes. Target light carried through the fiber system reaches one of four identical fast-Schmidt spectrograph modules, each with three arms. Prototype VPH gratings have been optically tested. CCD production is complete, with standard fully-depleted CCDs for red arms and more-challenging thinner fully-depleted CCDs with blue-optimized coating for blue arms.

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 Wilkinson Microwave Anisotropy Probe (WMAP) mapped the distribution of temperature and polarization over the entire sky in five microwave frequency bands. These full-sky maps were used to obtain measurements of temperature and polarization anisotropy of the cosmic microwave background with the unprecedented accuracy and precision. The analysis of two-point correlation functions of temperature and polarization data gives determinations of the fundamental cosmological parameters such as the age and composition of the universe, as well as the key parameters describing the physics of inflation, which is further constrained by three-point correlation functions. WMAP observations alone reduced the flat $\Lambda$ cold dark matter ($\Lambda$CDM) cosmological model (six) parameter volume by a factor of >68,000 compared with pre-WMAP measurements. The WMAP observations (sometimes in combination with other astrophysical probes) convincingly show the existence of non-baryonic dark matter, the cosmic neutrino background, flatness of spatial geometry of the universe, a deviation from a scale-invariant spectrum of initial scalar fluctuations, and that the current universe is undergoing an accelerated expansion. The WMAP observations provide the strongest ever support for inflation; namely, the structures we see in the universe originate from quantum fluctuations generated during inflation.

We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ~ 4 Jupiter masses hosting a sub-Earth mass moon. The data are well fit by this exomoon model, but an alternate star+planet model fits the data almost as well. Nevertheless, these results indicate the potential of microlensing to detect exomoons, albeit ones that are different from the giant planet moons in our solar system. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M pi_rel, where M is the lens system mass and pi_rel is the lens-source relative parallax. If the lens system is nearby (large pi_rel), then M is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, mu_rel = 19.6 +- 1.6 mas/yr, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, the alternate solution has a lower proper motion, which is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M_host = 0.12 (+0.19 -0.06) M_solar and m_comp = 18 (+28 -100 M_earth, at a projected separation of a_perp ~ 0.84 AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the low mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.

We analyze MOA-2010-BLG-311, a high magnification (A_max>600) microlensing event with complete data coverage over the peak, making it very sensitive to planetary signals. We fit this event with both a point lens and a 2-body lens model and find that the 2-body lens model is a better fit but with only Delta chi^2~80. The preferred mass ratio between the lens star and its companion is $q=10^(-3.7+/-0.1), placing the candidate companion in the planetary regime. Despite the formal significance of the planet, we show that because of systematics in the data the evidence for a planetary companion to the lens is too tenuous to claim a secure detection. When combined with analyses of other high-magnification events, this event helps empirically define the threshold for reliable planet detection in high-magnification events, which remains an open question.

We analyze the planetary microlensing event MOA-2010-BLG-328. The best fit yields host and planetary masses of Mh = 0.11+/-0.01 M_{sun} and Mp = 9.2+/-2.2M_Earth, corresponding to a very late M dwarf and sub-Neptune-mass planet, respectively. The system lies at DL = 0.81 +/- 0.10 kpc with projected separation r = 0.92 +/- 0.16 AU. Because of the host's a-priori-unlikely close distance, as well as the unusual nature of the system, we consider the possibility that the microlens parallax signal, which determines the host mass and distance, is actually due to xallarap (source orbital motion) that is being misinterpreted as parallax. We show a result that favors the parallax solution, even given its close host distance. We show that future high-resolution astrometric measurements could decisively resolve the remaining ambiguity of these solutions.

We present 450 {\mu}m polarimetric observations of the M17 molecular cloud obtained with the SHARP polarimeter at the Caltech Submillimeter Observatory. Across the observed region, the magnetic field orientation is consistent with previous submillimeter and far-infrared polarization measurements. Our observations are centered on a region of the molecular cloud that has been compressed by stellar winds from a cluster of OB stars. We have compared these new data with previous 350 {\mu}m polarimetry and find an anti-correlation between the 450 to 350 {\mu}m polarization magnitude ratio and the ratio of 21 cm to 450 {\mu}m intensity. The polarization ratio is lower near the east end of the studied region where the cloud is exposed to stellar winds and radiation. At the west end of the region, the polarization ratio is higher. We interpret the varying polarization spectrum as evidence supporting the radiative alignment torque (RAT) model for grain alignment, implying higher alignment efficiency in the region that is exposed to a higher anisotropic radiation field.

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 Cosmology Large Angular Scale Surveyor (CLASS) instrument will measure the polarization of the cosmic microwave background at 40, 90, and 150 GHz from Cerro Toco in the Atacama desert of northern Chile. In this paper, we describe the optical design of the 40 GHz telescope system. The telescope is a diffraction limited catadioptric design consisting of a front-end Variable-delay Polarization Modulator (VPM), two ambient temperature mirrors, two cryogenic dielectric lenses, thermal blocking filters, and an array of 36 smooth-wall scalar feedhorn antennas. The feed horns guide the signal to antenna-coupled transition-edge sensor (TES) bolometers. Polarization diplexing and bandpass definition are handled on the same microchip as the TES. The feed horn beams are truncated with 10 dB edge taper by a 4 K Lyot-stop to limit detector loading from stray light and control the edge illumination of the front-end VPM. The field-of-view is 19deg x 14deg with a resolution for each beam on the sky of 1.5deg FWHM.

The Galactic bulge source MOA-2010-BLG-523S exhibited short-term deviations from a standard microlensing lightcurve near the peak of an Amax ~ 265 high-magnification microlensing event. The deviations originally seemed consistent with expectations for a planetary companion to the principal lens. We combine long-term photometric monitoring with a previously published high-resolution spectrum taken near peak to demonstrate that this is an RS CVn variable, so that planetary microlensing is not required to explain the lightcurve deviations. This is the first spectroscopically confirmed RS CVn star discovered in the Galactic bulge.

The Prime Focus Spectrograph (PFS) of the Subaru Measurement of Images and Redshifts (SuMIRe) project has been endorsed by Japanese community as one of the main future instruments of the Subaru 8.2-meter telescope at Mauna Kea, Hawaii. This optical/near-infrared multi-fiber spectrograph targets cosmology with galaxy surveys, Galactic archaeology, and studies of galaxy/AGN evolution. Taking advantage of Subaru's wide field of view, which is further extended with the recently completed Wide Field Corrector, PFS will enable us to carry out multi-fiber spectroscopy of 2400 targets within 1.3 degree diameter. A microlens is attached at each fiber entrance for F-ratio transformation into a larger one so that difficulties of spectrograph design are eased. Fibers are accurately placed onto target positions by positioners, each of which consists of two stages of piezo-electric rotary motors, through iterations by using back-illuminated fiber position measurements with a wide-field metrology camera. Fibers then carry light to a set of four identical fast-Schmidt spectrographs with three color arms each: the wavelength ranges from 0.38 {\mu}m to 1.3 {\mu}m will be simultaneously observed with an average resolving power of 3000. Before and during the era of extremely large telescopes, PFS will provide the unique capability of obtaining spectra of 2400 cosmological/astrophysical targets simultaneously with an 8-10 meter class telescope. The PFS collaboration, led by IPMU, consists of USP/LNA in Brazil, Caltech/JPL, Princeton, & JHU in USA, LAM in France, ASIAA in Taiwan, and NAOJ/Subaru.

We describe the conceptual design of the camera cryostats, detectors, and detector readout electronics for the SuMIRe Prime Focus Spectrograph (PFS) being developed for the Subaru telescope. The SuMIRe PFS will consist of four identical spectrographs, each receiving 600 fibers from a 2400 fiber robotic positioner at the prime focus. Each spectrograph will have three channels covering wavelength ranges 3800 {\AA} - 6700 {\AA}, 6500 {\AA} - 10000 {\AA}, and 9700 {\AA} - 13000 {\AA}, with the dispersed light being imaged in each channel by a f/1.10 vacuum Schmidt camera. In the blue and red channels a pair of Hamamatsu 2K x 4K edge-buttable CCDs with 15 um pixels are used to form a 4K x 4K array. For the IR channel, the new Teledyne 4K x 4K, 15 um pixel, mercury-cadmium-telluride sensor with substrate removed for short-wavelength response and a 1.7 um cutoff will be used. Identical detector geometry and a nearly identical optical design allow for a common cryostat design with the only notable difference being the need for a cold radiation shield in the IR camera to mitigate thermal background. This paper describes the details of the cryostat design and cooling scheme, relevant thermal considerations and analysis, and discusses the detectors and detector readout electronics.

We present a method to extract the redshift-space distortions beta parameter in configuration space with a minimal set of cosmological assumptions. We show that a novel combination of the observed monopole and quadrupole correlation functions can remove efficiently the impact of mild non linearities and redshift errors. The method offers a series of convenient properties: it does not depend on the theoretical linear correlation function, the mean galaxy density is irrelevant, only convolutions are used, there is no explicit dependence on linear bias. Analyses based on dark matter N-body simulations and Fisher matrix demonstrate that errors of a few percent on beta are possible with a full sky, 1(Gpc/h)^3 survey centered at a redshift of unity and with negligible shot noise. We also find a baryonic feature in the normalized quadrupole in configuration space that should complicate the extraction of the growth parameter from the linear theory asymptote, but that does not have a major impact with our method.

Microlensing detections of cool planets are important for the construction of an unbiased sample to estimate the frequency of planets beyond the snow line, which is where giant planets are thought to form according to the core accretion theory of planet formation. In this paper, we report the discovery of a giant planet detected from the analysis of the light curve of a high-magnification microlensing event MOA-2010-BLG-477. The measured planet-star mass ratio is $q=(2.181\pm0.004)\times 10^{-3}$ and the projected separation is $s=1.1228\pm0.0006$ in units of the Einstein radius. The angular Einstein radius is unusually large $\theta_{\rm E}=1.38\pm 0.11$ mas. Combining this measurement with constraints on the "microlens parallax" and the lens flux, we can only limit the host mass to the range $0.13<M/M_\odot<1.0$. In this particular case, the strong degeneracy between microlensing parallax and planet orbital motion prevents us from measuring more accurate host and planet masses. However, we find that adding Bayesian priors from two effects (Galactic model and Keplerian orbit) each independently favors the upper end of this mass range, yielding star and planet masses of $M_*=0.67^{+0.33}_{-0.13}\ M_\odot$ and $m_p=1.5^{+0.8}_{-0.3}\ M_{\rm JUP}$ at a distance of $D=2.3\pm0.6$ kpc, and with a semi-major axis of $a=2^{+3}_{-1}$ AU. Finally, we show that the lens mass can be determined from future high-resolution near-IR adaptive optics observations independently from two effects, photometric and astrometric.

We present the analysis of the light curves of 9 high-magnification single-lens gravitational microlensing events with lenses passing over source stars, including OGLE-2004-BLG-254, MOA-2007-BLG-176, MOA-2007-BLG-233/OGLE-2007-BLG-302, MOA-2009-BLG-174, MOA-2010-BLG-436, MOA-2011-BLG-093, MOA-2011-BLG-274, OGLE-2011-BLG-0990/MOA-2011-BLG-300, and OGLE-2011-BLG-1101/MOA-2011-BLG-325. For all events, we measure the linear limb-darkening coefficients of the surface brightness profile of source stars by measuring the deviation of the light curves near the peak affected by the finite-source effect. For 7 events, we measure the Einstein radii and the lens-source relative proper motions. Among them, 5 events are found to have Einstein radii less than 0.2 mas, making the lenses candidates of very low-mass stars or brown dwarfs. For MOA-2011-BLG-274, especially, the small Einstein radius of $\theta_{\rm E}\sim 0.08$ mas combined with the short time scale of $t_{\rm E}\sim 2.7$ days suggests the possibility that the lens is a free-floating planet. For MOA-2009-BLG-174, we measure the lens parallax and thus uniquely determine the physical parameters of the lens. We also find that the measured lens mass of $\sim 0.84\ M_\odot$ is consistent with that of a star blended with the source, suggesting that the blend is likely to be the lens. Although we find planetary signals for none of events, we provide exclusion diagrams showing the confidence levels excluding the existence of a planet as a function of the separation and mass ratio.

We present a one-port calibration technique for characterization of beam waveguide components with a vector network analyzer. This technique involves using a set of known delays to separate the responses of the instrument and the device under test. We demonstrate this technique by measuring the reflected performance of a millimeter-wave variable-delay polarization modulator.

We present the discovery and mass measurement of the cold, low-mass planet MOA-2009-BLG-266Lb, made with the gravitational microlensing method. This planet has a mass of m_p = 10.4 +- 1.7 Earth masses and orbits a star of mass M_* = 0.56 +- 0.09 Solar masses at a semi-major axis of a = 3.2 (+1.9 -0.5) AU and an orbital period of P = 7.6 (+7.7 -1.5} yrs. The planet and host star mass measurements are enabled by the measurement of the microlensing parallax effect, which is seen primarily in the light curve distortion due to the orbital motion of the Earth. But, the analysis also demonstrates the capability to measure microlensing parallax with the Deep Impact (or EPOXI) spacecraft in a Heliocentric orbit. The planet mass and orbital distance are similar to predictions for the critical core mass needed to accrete a substantial gaseous envelope, and thus may indicate that this planet is a "failed" gas giant. This and future microlensing detections will test planet formation theory predictions regarding the prevalence and masses of such planets.

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

A gravitational lens system can be perturbed by "rogue systems" in angular proximities but at different distances. A point mass perturbed by another point mass can be considered as a large separation approximation of the double scattering two point mass (DSTP) lens. The resulting effective lens depends on whether the perturber is closer to or farther from the observer than the main lens system. The caustic is smaller than that of the large separation binary lens when the perturber is the first scatterer; the caustic is similar in size with the large separation binary lens when the perturber is the last scatterer. Modelling of a gravitational lensing by a galaxy requires extra terms other than constant shear for the perturbers at different redshifts. Double scattering two distributed mass (DSTD) lens is considered. The perturbing galaxy behaves as a monopole -- or a point mass -- because the dipole moment of the elliptic mass distribution is zero.

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.

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.

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

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

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

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

Albert Einstein postulated the equivalence of energy and mass, developed the theory of special relativity, explained the photoelectric effect, and described Brownian motion in five papers, all published in 1905, 100 years ago. With these papers, Einstein provided the framework for understanding modern astrophysical phenomena. Conversely, astrophysical observations provide one of the most effective means for testing Einstein's theories. Here, I review astrophysical advances precipitated by Einstein's insights, including gravitational redshifts, gravitational lensing, gravitational waves, the Lense-Thirring effect, and modern cosmology. A complete understanding of cosmology, from the earliest moments to the ultimate fate of the universe, will require developments in physics beyond Einstein, to a unified theory of gravity and quantum physics.

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.

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.

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.

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.

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

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.

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

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.

We are surrounded by radiation that originated from the big bang. It has traveled to us from the farthest reaches of the Universe, carrying with it an unaltered record of the beginning of time and space. The radiation is in the form of gravitational waves - propagating ripples in the curvature of spacetime. We describe a mission to detect these Gravitational Echos Across Time (GREAT) that would open up a new window on the very early universe. By studying the gravitational echoes of the big bang we will gain insight into the fundamental structure of matter, gravity, and how the Universe formed.

Currently three imaging spectrometer architectures, tunable filter, dispersive, and Fourier transform, are viable for imaging the universe in three dimensions. There are domains of greatest utility for each of these architectures. The optimum choice among the various alternative architectures is dependent on the nature of the desired observations, the maturity of the relevant technology, and the character of the backgrounds. The domain appropriate for each of the alternatives is delineated; both for instruments having ideal performance as well as for instrumentation based on currently available technology. The environment and science objectives for the Next Generation Space Telescope will be used as a specific representative case to provide a basis for comparison of the various alternatives.

The COBE Far InfraRed Absolute Spectrophotometer (FIRAS) observations constitute an unbiased survey over the wavelength range from 100 um to 1 cm over 99% of the sky. Improved calibration of the FIRAS instrument and the inclusion of all of the FIRAS data allow an improved signal to noise determination of the spectral lines by a factor of ~2 over our previous results. The resolution is low (0.45 cm^{-1}) so only the strongest lines are observable. The CO chain from J=1-0 to J=8-7 is observed towards the Galactic center. The line ratios are roughly consistent with a 40 K excitation temperature. The 157.7 um C II and 205.3 um N II lines are observable over most of the sky. The 370.4 um and 609.1 um lines of C I, and the 121.9 um line of N II are observed in the Galactic plane. The line ratios at the Galactic center are consistent with a density of n_0 ~30 cm^{-3} and a UV flux of G_0 ~15 uW m^{-2} sr^{-1} (10 Habing units). The 269 um H2O line is observed toward the Galactic center in absorption.

We describe a concept for an imaging spectrograph for a large orbiting observatory such as NASA's proposed Next Generation Space Telescope (NGST) based on an imaging Fourier transform spectrograph (IFTS). An IFTS has several important advantages which make it an ideal instrument to pursue the scientific objectives of NGST. We review the operation of an IFTS and make a quantitative evaluation of the signal-to-noise performance of such an instrument in the context of NGST. We consider the relationship between pixel size, spectral resolution, and diameter of the beamsplitter for imaging and non-imaging Fourier transform spectrographs and give the condition required to maintain spectral modulation efficiency over the entire field of view. We give examples of scientific programs that could be performed with this facility.

The COBE FIRAS data contain foreground emission from interplanetary, Galactic interstellar dust and extragalactic background emission. We use three different methods to separate the various emission components, and derive the spectrum of the extragalactic Far InfraRed Background (FIRB). Each method relies on a different set of assumptions, which affect the FIRB spectrum in different ways. Despite this, the FIRB spectra derived by these different methods are remarkably similar. The average spectrum that we derive in the \nu = 5 - 80 cm^{-1} (2000-125 \um) frequency interval is: I(\nu) = (1.3 \pm 0.4)\times10^{-5} (\nu /\nu_0)^{0.64 \pm 0.12} P_{\nu}(18.5 \pm 1.2 K), where \nu_0=100 cm^{-1} (\lambda_0=100 \um) and P is the Planck function. The derived FIRB spectrum is consistent with the DIRBE 140 and 240 \um detections. The total intensity received in the 5 - 80 cm^{-1} frequency interval is 14 nW m^{-2} sr^{-1}, and comprises about 20% of the total intensity expected from the energy release from nucleosynthesis throughout the history of the universe.

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.

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

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.

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.

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.

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 derive Galactic continuum spectra from 5-96/cm from COBE/FIRAS observations. The spectra are dominated by warm dust emission, which may be fit with a single temperature in the range 16-21 K (for nu^2 emissivity) along each line of sight. Dust heated by the attenuated radiation field in molecular clouds gives rise to intermediate temperature (10-14 K) emission in the inner Galaxy only. A widespread, very cold component (4-7 K) with optical depth that is spatially correlated with the warm component is also detected. The cold component is unlikely to be due to very cold dust shielded from starlight, because it is present at high latitude. We consider hypotheses that the cold component is due to enhanced submillimeter emissivity of the dust that gives rise to the warm component, or that it may be due to very small, large, or fractal particles. Lack of substantial power above the emission from warm dust places strong constraints on the amount of cold gas in the Galaxy. The microwave sky brightness due to interstellar dust is dominated by the cold component, and its angular variation could limit our ability to discern primordial fluctuations in the cosmic microwave background radiation.

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 three papers summarize the methods employed by the COBE-DMR group to analyze the power spectrum of the two-year DMR anisotropy maps. Results for the slope and normalization of the power spectrum are presented and compared. The implications of these results for various dark matter-based structure formation models are also discussed.

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.

The properties of the cosmic microwave background radiation provide unique constraints on cosmological models, i.e. on the content, history, and evolution of the Universe. I discuss the latest measurements of the spectral and spatial properties of the cosmic microwave background radiation. Recent measurements from NASA's Cosmic Background Explorer (COBE) satellite, and from balloon-borne and ground-based platforms, are summarized and their cosmological implications are discussed.

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.

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.

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.

The angular power spectrum estimator developed by Peebles (1973) and Hauser & Peebles (1973) has been modified and applied to the 2 year maps produced by the COBE DMR. The power spectrum of the real 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$. Within the limited range of spatial scales covered by the COBE DMR, corresponding to spherical harmonic indices $3 \leq \ell \lsim 30$, the best fitting value of the spectral index is $n = 1.25^{+0.4}_{-0.45}$ with the Harrison-Zeldovich value $n = 1$ approximately 0.5$\sigma$ below the best fit. For $3 \leq \ell \lsim 19$, the best fit is $n = 1.46^{+0.39}_{-0.44}$. Comparing the COBE DMR $\Delta T/T$ at small $\ell$ to the $\Delta T/T$ at $\ell \approx 50$ from degree scale anisotropy experiments gives a smaller range of acceptable spectral indices which includes $n = 1$.

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 present a determination of the cosmic microwave background dipole amplitude and direction from the COBE Differential Microwave Radiometers (DMR) first year of data. Data from the six DMR channels are consistent with a Doppler-shifted Planck function of dipole amplitude Delta T = 3.365 +/-0.027 mK toward direction (l,b) = (264.4 +/- 0.3 deg, 48.4 +/- 0.5 deg). The implied velocity of the Local Group with respect to the CMB rest frame is 627 +/- 22 km/s toward (l,b) = (276 +/- 3 deg, 30 +/- 3 deg). DMR has also mapped the dipole anisotropy resulting from the Earth's orbital motion about the Solar system barycenter, yielding a measurement of the monopole CMB temperature at 31.5, 53, and 90 GHz, to be 2.75 +/- 0.05 K.

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.