We report the discovery of a unique gravitational lens system, SDSSJ2222+2745, producing five spectroscopically confirmed images of a z_s=2.82 quasar lensed by a foreground galaxy cluster at z_l=0.49. We also present photometric and spectroscopic evidence for a sixth lensed image of the same quasar. The maximum separation between the quasar images is 15.1". Both the large image separations and the high image multiplicity of the lensed quasar are in themselves exceptionally rare, and observing the combination of these two factors is an exceptionally unlikely occurrence in present datasets. This is only the third known case of a quasar lensed by a cluster, and the only one with six images. The lens system was discovered in the course of the Sloan Giant Arcs Survey, in which we identify candidate lenses in the Sloan Digital Sky Survey and target these for follow up and verification with the 2.56m Nordic Optical Telescope. Multi-band photometry obtained over multiple epochs from September 2011 to September 2012 reveal significant variability at the ~10-30% level in some of the quasar images, indicating that measurements of the relative time delay between quasar images will be feasible. In this lens system we also identify a bright (g = 21.5) giant arc corresponding to a strongly lensed background galaxy at z_s=2.30. We fit parametric models of the lens system, constrained by the redshift and positions of the quasar images and the redshift and position of the giant arc. The predicted time delays between different pairs of quasar images range from ~100 days to ~6 years.

We generate simulations of the CMB temperature field as observed by the WMAP satellite, taking into account the detailed shape of the asymmetric beams and scanning strategy of the experiment, and use these to re-estimate the WMAP beam transfer functions. This method avoids the need of artificially symmetrizing the beams, as done in the baseline WMAP approach, and instead measures the total convolution effect by direct simulation. We find only small differences with respect to the nominal transfer functions, typically less than 1% everywhere, and less than 0.5% at l<400. The net effect on the CMB power spectrum is less than 0.6%. The effect on all considered cosmological parameters are negligible. For instance, we find that the spectral index of scalar perturbations after taking into account the beam asymmetries is n_s = 0.964 +- 0.014, corresponding to a negative shift of -0.1 sigma compared to the previously released WMAP results. Our CMB sky simulations are made publicly available, and can be used for general studies of asymmetric beam effects in the WMAP data.

Non-commutative geometry at inflation can give arise to parity violating modulations of the primordial power spectrum. We develop the statistical tools needed for investigating whether these modulations are evident in the Cosmic Microwave Background (CMB). The free parameters of the models are two directional parameters (theta,phi), the signal amplitude A*, and a tilt parameter n* that modulates correlation power on different scales. The signature of the model corresponds to a kind of hemispherical power asymmetry. When analyzing the 7-year WMAP data we find a weak signature for a preferred direction in the Q-, V-, and W bands with direction (l,b) = (-225 deg,-25 deg) +- (20 deg, 20 deg), which is close to another previously discovered hemispherical power asymmetry. Although these results are intriguing, the significance of the detection in the W-, V- and Q-bands are nonzero at about 2 sigma, suggesting that the simplest parameterization of the leading correction represents only partially the effects of the space-time non-commutativity possibly responsible for the hemispherical asymmetry. Our constraints on the presence of a dipole are independent of its physical origin and prefer a blue-tilted spectral index n* ~ 0 with the amplitude A* ~ 0.18.

Using the 5 year WMAP data, we re-investigate claims of non-Gaussianities and asymmetries detected in local curvature statistics of the 1 year WMAP data. In Hansen et al 2004, it was found that the northern ecliptic hemisphere was non-Gaussian at the ~1% level testing the densities of hill-, lake and saddle points based on the second derivatives of the CMB temperature map. The 5 year WMAP data has a much lower noise level and better control of systematics. Using these, we find that the anomalies are still present at a consistent level. Also the direction of maximum non-Gaussianity remains. Due to limited availability of computer resources, Hansen et al. 2004 were unable to calculate the full covariance matrix for the chi^2 test used. Here we apply the full covariance matrix instead of the diagonal approximation and find that the non-Gaussianities disappear and there is no preferred non-Gaussian direction. We compare with simulations of weak lensing to see if this may cause the observed non-Gaussianity when using diagonal covariance matrix. We conclude that weak lensing does not produce non-Gaussianity in the local curvature statistics at the scales investigated in this paper. The cause of the non-Gaussian detection in the case of a diagonal matrix remains unclear.

We revisit the anisotropic universe model previously developed by Ackerman, Carroll and Wise (ACW), and generalize both the theoretical and computational framework to include polarization and various forms of systematic effects. We apply our new tools to simulated WMAP data in order to understand the potential impact of asymmetric beams, noise mis-estimation and potential Zodiacal light emission. We find that neither has any significant impact on the results. We next show that the previously reported ACW signal is also present in the 1-year WMAP temperature sky map presented by Liu & Li, where data cuts are more aggressive. Finally, we reanalyze the 5-year WMAP data taking into account a previously neglected (-i)^{l-l'}-term in the signal covariance matrix. We still find a strong detection of a preferred direction in the temperature map. Including multipoles up to l=400, the anisotropy amplitude for the W-band is found to be g = 0.29 +- 0.031, nonzero at 9 sigma. However, the corresponding preferred direction is also shifted very close to the ecliptic poles at (l,b)= (96,30), in agreement with the analysis of Hanson & Lewis, indicating that the signal is aligned along the plane of the solar system. This strongly suggests that the signal is not of cosmological origin, but most likely is a product of an unknown systematic effect. Determining the nature of the systematic effect is of vital importance, as it might affect other cosmological conclusions from the WMAP experiment. Finally, we provide a forecast for the Planck experiment including polarization.

We present local extrema studies of two models that introduce a preferred direction into the observed cosmic microwave background (CMB) temperature field. In particular, we make a frequentist comparison of the one- and two-point statistics for the dipole modulation and ACW models with data from the five-year Wilkinson Microwave Anisotropy Probe (WMAP). This analysis is motivated by previously revealed anomalies in the WMAP data, and particularly the difference in the statistical nature of the temperature anisotropies when analysed in hemispherical partitions. The analysis of the one-point statistics indicates that the previously determined hemispherical variance difficulties can be apparently overcome by a dipole modulation field, but new inconsistencies arise if the mean and the l-dependence of the statistics are considered. The two-point correlation functions of the local extrema, the temperature pair product and the point-point spatial pair-count, demonstrate that the impact of such a modulation is to over-`asymmetrise' the temperature field on smaller scales than the wave-length of the dipole or quadrupole, and this is disfavored by the observed data.The results from the ACW model predictions, however, are consistent with the standard isotropic hypothesis. The two-point analysis confirms that the impact of this type of violation of isotropy on the temperature extrema statistics is relatively weak. From this work, we conclude that a model with more spatial structure than the dipole modulated or rotational-invariance breaking models are required to fully explain the observed large-scale anomalies in the WMAP data.

We develop a new Bayesian method for estimating white noise levels in CMB sky maps, and apply this algorithm to the 5-year WMAP data. We assume that the amplitude of the noise RMS is scaled by a constant value, alpha, relative to a pre-specified noise level. We then derive the corresponding conditional density, P(alpha | s, Cl, d), which is subsequently integrated into a general CMB Gibbs sampler. We first verify our code by analyzing simulated data sets, and then apply the framework to the WMAP data. For the foreground-reduced 5-year WMAP sky maps and the nominal noise levels initially provided in the 5-year data release, we find that the posterior means typically range between alpha=1.005 +- 0.001 and alpha=1.010 +- 0.001 depending on differencing assembly, indicating that the noise level of these maps are biased low by 0.5-1.0%. The same problem is not observed for the uncorrected WMAP sky maps. After the preprint version of this letter appeared on astro-ph., the WMAP team has corrected the values presented on their web page, noting that the initially provided values were in fact estimates from the 3-year data release, not from the 5-year estimates. However, internally in their 5-year analysis the correct noise values were used, and no cosmological results are therefore compromised by this error. Thus, our method has already been demonstrated in practice to be both useful and accurate.

We present a consistent self-contained and pedagogical review of the CMB Gibbs sampler, focusing on computational methods and code design. We provide an easy-to-use CMB Gibbs sampler named SLAVE developed in C++ using object-oriented design. While discussing why the need for a Gibbs sampler is evident and what the Gibbs sampler can be used for in a cosmological context, we review in detail the analytical expressions for the conditional probability densities and discuss the problems of galactic foreground removal and anisotropic noise. Having demonstrated that SLAVE is a working, usable CMB Gibbs sampler, we present the algorithm for white noise level estimation. We then give a short guide on operating SLAVE before introducing the post-processing utilities for obtaining the best-fit power spectrum using the Blackwell-Rao estimator.

We introduce a new CMB temperature likelihood approximation called the Gaussianized Blackwell-Rao (GBR) estimator. This estimator is derived by transforming the observed marginal power spectrum distributions obtained by the CMB Gibbs sampler into standard univariate Gaussians, and then approximate their joint transformed distribution by a multivariate Gaussian. The method is exact for full-sky coverage and uniform noise, and an excellent approximation for sky cuts and scanning patterns relevant for modern satellite experiments such as WMAP and Planck. A single evaluation of this estimator between l=2 and 200 takes ~0.2 CPU milliseconds, while for comparison, a single pixel space likelihood evaluation between l=2 and 30 for a map with ~2500 pixels requires ~20 seconds. We apply this tool to the 5-year WMAP temperature data, and re-estimate the angular temperature power spectrum, $C_{\ell}$, and likelihood, L(C_l), for l<=200, and derive new cosmological parameters for the standard six-parameter LambdaCDM model. Our spectrum is in excellent agreement with the official WMAP spectrum, but we find slight differences in the derived cosmological parameters. Most importantly, the spectral index of scalar perturbations is n_s=0.973 +/- 0.014, 1.9 sigma away from unity and 0.6 sigma higher than the official WMAP result, n_s = 0.965 +/- 0.014. This suggests that an exact likelihood treatment is required to higher l's than previously believed, reinforcing and extending our conclusions from the 3-year WMAP analysis. In that case, we found that the sub-optimal likelihood approximation adopted between l=12 and 30 by the WMAP team biased n_s low by 0.4 sigma, while here we find that the same approximation between l=30 and 200 introduces a bias of 0.6 sigma in n_s.

We extend the previously described CMB Gibbs sampling framework to allow for exact Bayesian analysis of anisotropic universe models, and apply this method to the 5-year WMAP temperature observations. This involves adding support for non-diagonal signal covariance matrices, and implementing a general spectral parameter MCMC sampler. As a worked example we apply these techniques to the model recently introduced by Ackerman et al., describing for instance violations of rotational invariance during the inflationary epoch. After verifying the code with simulated data, we analyze the foreground-reduced 5-year WMAP temperature sky maps. For l <= 400 and the W-band data, we find tentative evidence for a preferred direction pointing towards (l,b) = (110 deg, 10 deg) with an anisotropy amplitude of g_* = 0.15 +/- 0.039. Similar results are obtained from the V-band data [g_* = 0.10 +/- 0.04; (l,b) = (130 deg, 20 deg)]. Further, the preferred direction is stable with respect to multipole range, seen independently in both l=[2,100] and [100,400], although at lower statistical significance. We have not yet been able to establish a fully satisfactory explanation for the observations in terms of known systematics, such as non-cosmological foregrounds, correlated noise or asymmetric beams, but stress that further study of all these issues is warranted before a cosmological interpretation can be supported.

Quantum gravity effects are expected to modify the primordial density fluctuations produced during inflation and leave their imprint on the cosmic microwave background observed today. We present a new analysis discussing whether these effects are detectable, considering both currently available data and simulated results from an optimal CMB experiment. We find that the WMAP (Wilkinson Microwave Anisotropy Probe) data show no evidence for the particular signature considered in this work but give an upper bound on the parameters of the model. However, a hypothetical experiment shows that with proper data, the trans-Planckian effects should be detectable through alternate sampling methods. This fuzzy conclusion is a result of the nature of the oscillations, since they give rise to a likelihood hypersurface riddled with local maxima. A simple Bayesian analysis shows no significant evidence for the simulated data to prefer a trans-Planckian model. Conventional Markov chain Monte Carlo (MCMC) methods are not suitable for exploring this complicated landscape, but alternative methods are required to solve the problem. This, however, requires extremely high-precision data.

We investigate the possibility of using cosmological observations to probe and constrain an imperfect dark energy fluid. We consider a general parameterization of the dark energy component accounting for an equation of state, speed of sound and viscosity. We use present and future data from the cosmic microwave background radiation (CMB), large scale structures and supernovae type Ia. We find that both the speed of sound and viscosity parameters are difficult to nail down with the present cosmological data. Also, we argue that it will be hard to improve the constraints significantly with future CMB data sets. The implication is that a perfect fluid description might ultimately turn out to be a phenomenologically sufficient description of all the observational consequences of dark energy. The fundamental lesson is however that even then one cannot exclude, by appealing to observational evidence alone, the possibility of imperfectness in dark energy.