We constrain the spectral index of polarized synchrotron emission, $\beta_s$, by correlating the recently released 2.3 GHz S-Band Polarization All Sky Survey (S-PASS) data with the 23 GHz 9-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We sub-divide the S-PASS field, which covers the Southern Ecliptic hemisphere, into $15^{\circ}\times 15^{\circ}$ regions, and estimate the spectral index of polarized synchrotron emission within each region using a simple but robust T-T plot technique. Three different versions of the S-PASS data are considered, corresponding to either no correction for Faraday rotation; Faraday correction based on the rotation measure model presented by the S-PASS team; or Faraday correction based on a rotation measure model presented by Hutschenreuter and En{\ss}lin. We find that the correlation between S-PASS and WMAP is strongest when applying the S-PASS model. Adopting this correction model, we find that the mean spectral index of polarized synchrotron emission gradually steepens from $\beta_s\approx-2.8$ at low Galactic latitudes to $\beta_s\approx-3.2$ at high Galactic latitudes, in good agreement with previously published results. Finally, we consider two special cases defined by the BICEP2 and SPIDER fields, and obtain mean estimates of $\beta_{BICEP2}=-3.22\pm0.06$ and $\beta_{SPIDER}=-3.21\pm0.03$, respectively. Adopting the WMAP 23 GHz sky map bandpass filtered to including angular scales only between $2^{\circ}$ and $10^{\circ}$ as a spatial template, we constrain the root-mean-square synchrotron polarization amplitude to be less than $0.03\mu K$ ($0.009\mu K$) at 90 GHz (150 GHz) for the BICEP2 field, corresponding roughly to a tensor-to-scalar ratio of $r\lesssim0.02$ ($r\lesssim0.005$), respectively. Very similar constraints are obtained for the SPIDER field.

We develop a time-ordered data simulator and map-maker for the proposed PIXIE Fourier transform spectrometer and use them to investigate the impact of polarization leakage, imperfect collimation, elliptical beams, sub-pixel effects, correlated noise and spectrometer mirror jitter on the PIXIE data analysis. We find that PIXIE is robust to all of these effects, with the exception of mirror jitter which could become the dominant source of noise in the experiment if the jitter is not kept significantly below $0.1\mu m\sqrt{s}$. Source code is available at https://github.com/amaurea/pixie.

We apply a messenger field method to solve the linear minimum-variance mapmaking equation in the context of Cosmic Microwave Background (CMB) observations. In simulations, the method produces sky maps that converge significantly faster than those from a conjugate gradient descent algorithm with a diagonal preconditioner, even though the computational cost per iteration is similar. The messenger method recovers large scales in the map better than conjugate gradient descent, and yields a lower overall $\chi^2$. In the single, pencil beam approximation, each iteration of the messenger mapmaking procedure produces an unbiased map, and the iterations become more optimal as they proceed. A variant of the method can handle differential data or perform deconvolution mapmaking. The messenger method requires no preconditioner, but a high-quality solution needs a cooling parameter to control the convergence. We study the convergence properties of this new method, and discuss how the algorithm is feasible for the large data sets of current and future CMB experiments.

We present polarization measurements of extragalactic radio sources observed during the Cosmic Microwave Background polarization survey of the Q/U Imaging Experiment (QUIET), operating at 43 GHz (Q-band) and 95 GHz (W-band). We examine sources selected at 20 GHz from the public, $>$40 mJy catalog of the Australia Telescope (AT20G) survey. There are $\sim$480 such sources within QUIET's four low-foreground survey patches, including the nearby radio galaxies Centaurus A and Pictor A. The median error on our polarized flux density measurements is 30--40 mJy per Stokes parameter. At S/N $> 3$ significance, we detect linear polarization for seven sources in Q-band and six in W-band; only $1.3 \pm 1.1$ detections per frequency band are expected by chance. For sources without a detection of polarized emission, we find that half of the sources have polarization amplitudes below 90 mJy (Q-band) and 106 mJy (W-band), at 95% confidence. Finally, we compare our polarization measurements to intensity and polarization measurements of the same sources from the literature. For the four sources with WMAP and Planck intensity measurements $>1$ Jy, the polarization fraction are above 1% in both QUIET bands. At high significance, we compute polarization fractions as much as 10--20% for some sources, but the effects of source variability may cut that level in half for contemporaneous comparisons. Our results indicate that simple models---ones that scale a fixed polarization fraction with frequency---are inadequate to model the behavior of these sources and their contributions to polarization maps.

We estimate the spectral index, beta, of polarized synchrotron emission as observed in the 9 yr WMAP sky maps using two methods, linear regression ("T-T plot") and maximum likelihood. We partition the sky into 24 disjoint sky regions, and evaluate the spectral index for all polarization angles between 0 deg and 85 deg in steps of 5. Averaging over polarization angles, we derive a mean spectral index of beta_all-sky=-2.99+-0.01 in the frequency range of 23-33 GHz. We find that the synchrotron spectral index steepens by 0.14 from low to high Galactic latitudes, in agreement with previous studies, with mean spectral indices of beta_plane=-2.98+-0.01 and beta_high-lat=-3.12+-0.04. In addition, we find a significant longitudinal variation along the Galactic plane with a steeper spectral index toward the Galactic center and anticenter than toward the Galactic spiral arms. This can be well modeled by an offset sinusoidal, beta(l)=-2.85+0.17sin(2l-90). Finally, we study synchrotron emission in the BICEP2 field, in an attempt to understand whether the claimed detection of large-scale B-mode polarization could be explained in terms of synchrotron contamination. Adopting a spectral index of beta=-3.12, typical for high Galactic latitudes, we find that the most likely bias corresponds to about 2% of the reported signal (r=0.003). The flattest index allowed by the data in this region is beta=-2.5, and under the assumption of a straight power-law frequency spectrum, we find that synchrotron emission can account for at most 20% of the reported BICEP2 signal.

High-resolution Maximum Likelihood map-making of the Cosmic Microwave Background is usually performed using Conjugate Gradients with a preconditioner that ignores noise correlations. We here present a new preconditioner that approximates the map noise covariance as circulant, and show that this results in a speedup of up to 400% for a realistic scanning pattern from the Atacama Cosmology Telescope. The improvement is especially large for polarized maps.

We present the temperature power spectrum of the Cosmic Microwave Background obtained by cross-correlating maps from the Atacama Cosmology Telescope (ACT) at 148 and 218 GHz with maps from the Planck satellite at 143 and 217 GHz, in two overlapping regions covering 592 square degrees. We find excellent agreement between the two datasets at both frequencies, quantified using the variance of the residuals between the ACT power spectra and the ACTxPlanck cross-spectra. We use these cross-correlations to calibrate the ACT data at 148 and 218 GHz, to 0.7% and 2% precision respectively. We find no evidence for anisotropy in the calibration parameter. We compare the Planck 353 GHz power spectrum with the measured amplitudes of dust and cosmic infrared background (CIB) of ACT data at 148 and 218 GHz. We also compare planet and point source measurements from the two experiments.

Cosmic Microwave Background lensing simulation by Taylor expansion has long been considered impractical due to slow convergence, but a recent flat-sky implementation shows that a simple trick eliminates this problem, making Taylor lensing a fast and simple lensing algorithm for the flat sky. Here we generalize the method to the full sky, and study its convergence and performance relative to a commonly used numerical code, Lenspix, with extensive benchmarks of both. Compared to the flat sky case, the method takes a speed hit due to the slow speed of spherical harmonic transforms compared to fast Fourier transforms, resulting in speeds of 1/3 to 2/3 of Lenspix for similar accuracy.

We present efficient algorithms for CMB lensing simulation and power spectrum es- timation for flat-sky CMB polarization maps. We build a pure B-mode estimator to remedy E to B leakage due to partial sky coverage. We show that our estimators are unbiased, and consistent with the projected errors. We demonstrate our algorithm using simulated observations of small sky patches with realistic noise and weights for upcoming CMB polarization experiments.

We investigate sets of random variables that can be arranged sequentially such that a given variable only depends conditionally on its immediate predecessor. For such sets, we show that the full joint probability distribution may be expressed exclusively in terms of uni- and bivariate marginals. Under the assumption that the CMB power spectrum likelihood only exhibits correlations within a banded multipole range, \Delta l, we apply this expression to two outstanding problems in CMB likelihood analysis. First, we derive a statistically well-defined hybrid likelihood estimator, merging two independent (e.g., low- and high-l) likelihoods into a single expression that properly accounts for correlations between the two. Applying this expression to the WMAP likelihood, we verify that the effect of correlations on cosmological parameters in the transition region is negligible in terms of cosmological parameters for WMAP; the largest relative shift seen for any parameter is 0.06\sigma. However, because this may not hold for other experimental setups (e.g., for different instrumental noise properties or analysis masks), but must rather be verified on a case-by-case basis, we recommend our new hybridization scheme for future experiments for statistical self-consistency reasons. Second, we use the same expression to improve the convergence rate of the Blackwell-Rao likelihood estimator, reducing the required number of Monte Carlo samples by several orders of magnitude, and thereby extend it to high-l applications.

The Q/U Imaging ExperimenT (QUIET) has observed the cosmic microwave background (CMB) at 43 and 95GHz. The 43-GHz results have been published in QUIET Collaboration et al. (2011), and here we report the measurement of CMB polarization power spectra using the 95-GHz data. This data set comprises 5337 hours of observations recorded by an array of 84 polarized coherent receivers with a total array sensitivity of 87 uK sqrt(s). Four low-foreground fields were observed, covering a total of ~1000 square degrees with an effective angular resolution of 12.8', allowing for constraints on primordial gravitational waves and high-signal-to-noise measurements of the E-modes across three acoustic peaks. The data reduction was performed using two independent analysis pipelines, one based on a pseudo-Cl (PCL) cross-correlation approach, and the other on a maximum-likelihood (ML) approach. All data selection criteria and filters were modified until a predefined set of null tests had been satisfied before inspecting any non-null power spectrum. The results derived by the two pipelines are in good agreement. We characterize the EE, EB and BB power spectra between l=25 and 975 and find that the EE spectrum is consistent with LCDM, while the BB power spectrum is consistent with zero. Based on these measurements, we constrain the tensor-to-scalar ratio to r=1.1+0.9-0.8 (r<2.8 at 95% C.L.) as derived by the ML pipeline, and r=1.2+0.9-0.8 (r<2.7 at 95% C.L.) as derived by the PCL pipeline. In one of the fields, we find a correlation with the dust component of the Planck Sky Model, though the corresponding excess power is small compared to statistical errors. Finally, we derive limits on all known systematic errors, and demonstrate that these correspond to a tensor-to-scalar ratio smaller than r=0.01, the lowest level yet reported in the literature.

We present a fully parallelized grid-based parameter estimation algorithm for investigating multidimensional likelihoods called Snake, and apply it to cosmological parameter estimation. The basic idea is to map out the likelihood grid-cell by grid-cell according to decreasing likelihood, and stop when a certain threshold has been reached. This approach improves vastly on the "curse of dimensionality" problem plaguing standard grid-based parameter estimation simply by disregarding grid-cells with negligible likelihood. The main advantages of this method compared to standard Metropolis-Hastings MCMC methods include 1) trivial extraction of arbitrary conditional distributions; 2) direct access to Bayesian evidences; 3) better sampling of the tails of the distribution; and 4) nearly perfect parallelization scaling. The main disadvantage is, as in the case of brute-force grid-based evaluation, a dependency on the number of parameters, N_par. One of the main goals of the present paper is to determine how large N_par can be, while still maintaining reasonable computational efficiency; we find that N_par = 12$ is well within the capabilities of the method. The performance of the code is tested by comparing cosmological parameters estimated using Snake and the WMAP-7 data with those obtained using CosmoMC, the current standard code in the field. We find fully consistent results, with similar computational expenses, but shorter wall time due to the perfect parallelization scheme.

The Q/U Imaging ExperimenT (QUIET) employs coherent receivers at 43GHz and 95GHz, operating on the Chajnantor plateau in the Atacama Desert in Chile, to measure the anisotropy in the polarization of the CMB. QUIET primarily targets the B modes from primordial gravitational waves. The combination of these frequencies gives sensitivity to foreground contributions from diffuse Galactic synchrotron radiation. Between 2008 October and 2010 December, >10,000hours of data were collected, first with the 19-element 43GHz array (3458hours) and then with the 90-element 95GHz array. Each array observes the same four fields, selected for low foregrounds, together covering ~1000deg^2. This paper reports initial results from the 43GHz receiver which has an array sensitivity to CMB fluctuations of 69uK sqrt(s). The data were extensively studied with a large suite of null tests before the power spectra, determined with two independent pipelines, were examined. Analysis choices, including data selection, were modified until the null tests passed. Cross correlating maps with different telescope pointings is used to eliminate a bias. This paper reports the EE, BB and EB power spectra in the multipole range ell=25-475. With the exception of the lowest multipole bin for one of the fields, where a polarized foreground, consistent with Galactic synchrotron radiation, is detected with 3sigma significance, the E-mode spectrum is consistent with the LCDM model, confirming the only previous detection of the first acoustic peak. The B-mode spectrum is consistent with zero, leading to a measurement of the tensor-to-scalar ratio of r=0.35+1.06-0.87. The combination of a new time-stream double-demodulation technique, Mizuguchi-Dragone optics, natural sky rotation, and frequent boresight rotation leads to the lowest level of systematic contamination in the B-mode power so far reported, below the level of r=0.1

We describe a Bayesian framework for estimating the time-domain noise covariance of CMB observations, typically parametrized in terms of a 1/f frequency profile. This framework is based on the Gibbs sampling algorithm, which allows for exact marginalization over nuisance parameters through conditional probability distributions. In this paper we implement support for gaps in the data streams and marginalization over fixed time-domain templates, and also outline how to marginalize over confusion from CMB fluctuations, which may be important for high signal-to-noise experiments. As a by-product of the method, we obtain proper constrained realizations, which themselves can be useful for map making. To validate the algorithm, we demonstrate that the reconstructed noise parameters and corresponding uncertainties are unbiased using simulated data. The CPU time required to process a single data stream of 100 000 samples with 1000 samples removed by gaps is 3 seconds if only the maximum posterior parameters are required, and 21 seconds if one also want to obtain the corresponding uncertainties by Gibbs sampling.

A recent application of the Kolmogorov-Smirnov test to the WMAP 7 year W-band maps claims evidence that the CMB is "weakly random", and that only 20% of the signal can be explained as a random Gaussian field. I here repeat this analysis, and in contrast to the original result find no evidence for deviation from the standard {\Lambda}CDM model. Instead, the results of the original analysis are consistent with not properly taking into account the correlations of the {\Lambda}CDM power spectrum.