We extend a general maximum likelihood foreground estimation for cosmic microwave background polarization data to include estimation of instrumental systematic effects. We focus on two particular effects: frequency band measurement uncertainty, and instrumentally induced frequency dependent polarization rotation. We assess the bias induced on the estimation of the $B$-mode polarization signal by these two systematic effects in the presence of instrumental noise and uncertainties in the polarization and spectral index of Galactic dust. Degeneracies between uncertainties in the band and polarization angle calibration measurements and in the dust spectral index and polarization increase the uncertainty in the extracted CMB $B$-mode power, and may give rise to a biased estimate. We provide a quantitative assessment of the potential bias and increased uncertainty in an example experimental configuration. For example, we find that with 10\% polarized dust, tensor to scalar ratio of $r=0.05$, and the instrumental configuration of the EBEX balloon payload, the estimated CMB $B$-mode power spectrum is recovered without bias when the frequency band measurement has 5% uncertainty or less, and the polarization angle calibration has an uncertainty of up to 4$^{\circ}$.

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

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.

We show that a one-component variable-emissivity-spectral-index model (the free-{\alpha} model) provides more physically motivated estimates of dust temperature at the Galactic polar caps than one- or two-component fixed-emissivity-spectral-index models (fixed-{\alpha} models) for interstellar dust thermal emission at far-infrared and millimeter wavelengths. For the comparison we have fit all-sky one-component dust models with fixed or variable emissivity spectral index to a new and improved version of the 210-channel dust spectra from the COBE-FIRAS, the 100 - 240 {\mu}m maps from the COBE-DIRBE and the 94 GHz dust map from the WMAP. The best model, the free-{\alpha} model, is well constrained by data at 60-3000 GHz over 86 per cent of the total sky area. It predicts dust temperature (Tdust) to be 13.7-22.7 (+/-1.3) K, the emissivity spectral index ({\alpha}) to be 1.2 - 3.1 (+/-0.3) and the optical depth ({\tau}) to range 0.6 - 46 x 10^(-5) at \nu 0 = 900 GHz (\lambda 0 = 333 {\mu}m) with a 23 per cent uncertainty. Using these estimates, we present all-sky evidence for an inverse correlation between the emissivity spectral index and dust temperature, which fits the relation {\alpha} = 1/({\delta}+{\omega} Tdust) with {\delta} = -0.510 +/- 0.011 and {\omega} = 0.059 +/- 0.001. This best model will be useful to cosmic microwave background experiments for removing foreground dust contamination and it can serve as an all-sky extended-frequency reference for future higher resolution dust models.

We study the impact of the spectral dependence of the linear polarization rotation induced by an achromatic half-wave plate on measurements of cosmic microwave background polarization in the presence of astrophysical foregrounds. We focus on the systematic effects induced on the measurement of inflationary gravitational waves by uncertainties in the polarization and spectral index of Galactic dust. We find that for the experimental configuration and noise levels of the balloon-borne EBEX experiment, which has three frequency bands centered at 150, 250, and 410 GHz, a crude dust subtraction process mitigates systematic effects to below detectable levels for 10% polarized dust and tensor to scalar ratio of as low as r = 0.01. We also study the impact of uncertainties in the spectral response of the instrument. With a top-hat model of the spectral response for each band, characterized by band-center and band-width, and with the same crude dust subtraction process, we find that these parameters need to be determined to within 1 and 0.8 GHz at 150 GHz; 9 and 2.0 GHz at 250 GHz; and 20 and 14 GHz at 410 GHz, respectively. The approach presented in this paper is applicable to other optical elements that exhibit polarization rotation as a function of frequency.

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

The diffuse high-latitude H-alpha background is widely believed to be predominantly the result of in-situ recombination of ionized hydrogen in the warm interstellar medium of the Galaxy. Instead, we show that both a substantial fraction of the diffuse high-latitude H-alpha intensity in regions dominated by Galactic cirrus dust and much of the variance in the high-latitude H-alpha background are the result of scattering by interstellar dust of H-alpha photons originating elsewhere in the Galaxy. We provide an empirical relation, which relates the expected scattered H-alpha intensity to the IRAS 100um diffuse background intensity, applicable to about 81% of the entire sky. The assumption commonly made in reductions of CMB observations, namely that the observed all-sky map of diffuse H-alpha light is a suitable template for Galactic free-free foreground emission, is found to be in need of reexamination.

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

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

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

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.

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.

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

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

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

We have discovered 21 new Type Ia supernovae (SNe Ia) with the Hubble Space Telescope (HST) and have used them to trace the history of cosmic expansion over the last 10 billion years. These objects, which include 13 spectroscopically confirmed SNe Ia at z > 1, were discovered during 14 epochs of reimaging of the GOODS fields North and South over two years with the Advanced Camera for Surveys on HST. Together with a recalibration of our previous HST-discovered SNe Ia, the full sample of 23 SNe Ia at z > 1 provides the highest-redshift sample known. Combined with previous SN Ia datasets, we measured H(z) at discrete, uncorrelated epochs, reducing the uncertainty of H(z>1) from 50% to under 20%, strengthening the evidence for a cosmic jerk--the transition from deceleration in the past to acceleration in the present. The unique leverage of the HST high-redshift SNe Ia provides the first meaningful constraint on the dark energy equation-of-state parameter at z >1. The result remains consistent with a cosmological constant (w(z)=-1), and rules out rapidly evolving dark energy (dw/dz >>1). The defining property of dark energy, its negative pressure, appears to be present at z>1, in the epoch preceding acceleration, with ~98% confidence in our primary fit. Moreover, the z>1 sample-averaged spectral energy distribution is consistent with that of the typical SN Ia over the last 10 Gyr, indicating that any spectral evolution of the properties of SNe Ia with redshift is still below our detection threshold.

I discuss several issues that arise when trying to constrain the dark energy equation of state using correlations of the integrated Sachs-Wolfe effect with galaxy counts and lensing of the cosmic microwave background. These techniques are complementary to others such as galaxy shear surveys, and can use data that will already be obtained from currently planned observations. In regimes where cosmic variance and shot noise are the dominant sources of error, constraints could be made on the mean equation of state to within 0.33 and its first derivative to within 1.0. Perhaps more interesting is that the determination of dark energy parameters by these types of experiments depends strongly on the presence or absence of perturbations in the dark energy fluid.

The increasing precision of cosmological datasets is opening up new opportunities to test predictions from cosmic inflation. Here we study the impact of high precision constraints on the primordial power spectrum and show how a new generation of observations can provide impressive new tests of the slow-roll inflation paradigm, as well as produce significant discriminating power among different slow-roll models. In particular, we consider next-generation measurements of the Cosmic Microwave Background (CMB) temperature anisotropies and (especially) polarization, as well as new Lyman-$\alpha$ measurements that could become practical in the near future. We emphasize relationships between the slope of the power spectrum and its first derivative that are nearly universal among existing slow-roll inflationary models, and show how these relationships can be tested on several scales with new observations. Among other things, our results give additional motivation for an all-out effort to measure CMB polarization.