We perform a reconstruction of the coupling function between vacuum energy and geodesic cold dark matter using the latest observational data. We bin the interaction in seventeen redshift bins but use a correlation prior to prevent rapid, unphysical oscillations in the coupling function. This prior also serves to eliminate any dependence of the reconstruction on the binning method. We use two different forms of the correlation prior, finding that both give similar results for the reconstruction of the dark matter -- dark energy interaction. Calculating the Bayes factor for each case, we find no meaningful evidence for deviation from the null interacting case, i.e. $\Lambda$CDM, in our reconstruction.

Likelihood fitting to two-point clustering statistics made from galaxy surveys usually assumes a multivariate normal distribution for the measurements, with justification based on the central limit theorem given the large number of overdensity modes. However, this assumption cannot hold on the largest scales where the number of modes is low. Whilst more accurate distributions have previously been developed in idealized cases, we derive a procedure suitable for analysing measured monopole power spectra with window effects, stochastic shot noise, and the dependence of the covariance matrix on the model being fitted all taken into account. A data transformation is proposed to give an approximately Gaussian likelihood, with a variance--correlation decomposition of the covariance matrix to account for its cosmological dependence. By comparing with the modified-$t$ likelihood derived under the usual normality assumption, we find in numerical tests that our new procedure gives more accurate constraints on the local non-Gaussianity parameter $f_\mathrm{NL}$, which is sensitive to the large-scale power. A simple data analysis pipeline is provided for straightforward application of this new approach in preparation for forthcoming large galaxy surveys such as DESI and Euclid.

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while blind to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for $\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino mass density and including the 457 $\times$ 457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions, and from their combination obtain $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m = 0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 = 0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and $w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

We consider a self-consistent and physical approach to interacting dark energy models described by a Lagrangian, and identify a new class of models with variable dark energy sound speed. We show that if the interaction between dark energy in the form of quintessence and cold dark matter is purely momentum exchange this generally leads to a dark energy sound speed that deviates from unity. Choosing a specific sub-case, we study its phenomenology by investigating the effects of the interaction on the cosmic microwave background and linear matter power spectrum. We also perform a global fitting of cosmological parameters using CMB data, and compare our findings to $\Lambda$CDM.

We show how standard Newtonian N-body simulations can be interpreted in terms of the weak-field limit of general relativity by employing the recently developed Newtonian motion gauge. Our framework allows the inclusion of radiation perturbations and the non-linear evolution of matter. We show how to construct the weak-field metric by combining Newtonian simulations with results from Einstein-Boltzmann codes. We discuss observational effects on weak lensing and ray tracing, identifying important relativistic corrections.

A flat Friedman-Roberson-Walker universe dominated by a cosmological constant ($\Lambda$) and cold dark matter (CDM) has been the working model preferred by cosmologists since the discovery of cosmic acceleration. However, tensions of various degrees of significance are known to be present among existing datasets within the $\Lambda$CDM framework. In particular, the Lyman-$\alpha$ forest measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation Spectroscopic Survey (BOSS) prefers a smaller value of the matter density fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave background (CMB). Also, the recently measured value of the Hubble constant, $H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$, is $3.4\sigma$ higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$ inferred from the Planck CMB data. In this work, we investigate if these tensions can be interpreted as evidence for a non-constant dynamical dark energy (DE). Using the Kullback-Leibler (KL) divergence to quantify the tension between datasets, we find that the tensions are relieved by an evolving DE, with the dynamical DE model preferred at a $3.5\sigma$ significance level based on the improvement in the fit alone. While, at present, the Bayesian evidence for the dynamical DE is insufficient to favour it over $\Lambda$CDM, we show that, if the current best fit DE happened to be the true model, it would be decisively detected by the upcoming DESI survey.

We forecast constraints on neutral hydrogen (HI) and cosmological parameters using near-term intensity mapping surveys with instruments such as BINGO, MeerKAT, and the SKA, and Stage III and IV optical galaxy surveys. If foregrounds and systematic effects can be controlled - a problem which becomes much easier in cross-correlation - these surveys will provide exquisite measurements of the HI density and bias, as well as measurements of the growth of structure, the angular diameter distance, and the Hubble rate, over a wide range of redshift. We also investigate the possibility of detecting the late time ISW effect using the Planck satellite and forthcoming intensity mapping surveys, finding that a large sky survey with Phase 1 of the SKA can achieve a near optimal detection.

Initial conditions for (Newtonian) cosmological N-body simulations are usually set by re-scaling the present-day power spectrum obtained from linear (relativistic) Boltzmann codes to the desired initial redshift of the simulation. This back-scaling method can account for the effect of inhomogeneous residual thermal radiation at early times, which is absent in the Newtonian simulations. We analyse this procedure from a fully relativistic perspective, employing the recently-proposed Newtonian motion gauge framework. We find that N-body simulations for LambdaCDM cosmology starting from back-scaled initial conditions can be self-consistently embedded in a relativistic space-time with first-order metric potentials calculated using a linear Boltzmann code. This space-time coincides with a simple "N-body gauge" for z<50 for all observable modes. Care must be taken, however, when simulating non-standard cosmologies. As an example, we analyse the back-scaling method in a cosmology with decaying dark matter, and show that metric perturbations become large at early times in the back-scaling approach, indicating a breakdown of the perturbative description. We suggest a suitable "forwards approach" for such cases.

The properties of large underdensities in the distribution of galaxies in the Universe, known as cosmic voids, are potentially sensitive probes of fundamental physics. We use data from the MultiDark suite of N-body simulations and multiple halo occupation distribution mocks to study the relationship between galaxy voids, identified using a watershed void-finding algorithm, and the gravitational potential $\Phi$. We find that the majority of galaxy voids correspond to local density minima in larger-scale overdensities, and thus lie in potential wells. However, a subset of voids can be identified that closely trace maxima of the gravitational potential and thus stationary points of the velocity field. We identify a new void observable, $\lambda_v$, which depends on a combination of the void size and the average galaxy density contrast within the void, and show that it provides a good proxy indicator of the potential at the void location. A simple linear scaling of $\Phi$ as a function of $\lambda_v$ is found to hold, independent of the redshift and properties of the galaxies used as tracers of voids. We provide an accurate fitting formula to describe the spherically averaged potential profile $\Phi(r)$ about void centre locations. We discuss the importance of these results for the understanding of the evolution history of voids, and for their use in precision measurements of the integrated Sachs-Wolfe effect, gravitational lensing and peculiar velocity distortions in redshift space.

Small temperature anisotropies in the Cosmic Microwave Background can be sourced by density perturbations via the late-time integrated Sachs-Wolfe effect. Large voids and superclusters are excellent environments to make a localized measurement of this tiny imprint. In some cases excess signals have been reported. We probed these claims with an independent data set, using the first year data of the Dark Energy Survey in a different footprint, and using a different super-structure finding strategy. We identified 52 large voids and 102 superclusters at redshifts $0.2 < z < 0.65$. We used the Jubilee simulation to a priori evaluate the optimal ISW measurement configuration for our compensated top-hat filtering technique, and then performed a stacking measurement of the CMB temperature field based on the DES data. For optimal configurations, we detected a cumulative cold imprint of voids with $\Delta T_{f} \approx -5.0\pm3.7~\mu K$ and a hot imprint of superclusters $\Delta T_{f} \approx 5.1\pm3.2~\mu K$ ; this is $\sim1.2\sigma$ higher than the expected $|\Delta T_{f}| \approx 0.6~\mu K$ imprint of such super-structures in $\Lambda$CDM. If we instead use an a posteriori selected filter size ($R/R_{v}=0.6$), we can find a temperature decrement as large as $\Delta T_{f} \approx -9.8\pm4.7~\mu K$ for voids, which is $\sim2\sigma$ above $\Lambda$CDM expectations and is comparable to previous measurements made using SDSS super-structure data.

We present a new method for detection of the integrated Sachs-Wolfe (ISW) imprints of cosmic superstructures on the cosmic microwave background, based on a matched filtering approach. The expected signal-to-noise ratio for this method is comparable to that obtained from the full cross-correlation, and unlike other stacked filtering techniques it is not subject to an a posteriori bias. We apply this method to Planck CMB data using voids and superclusters identified in the CMASS galaxy data from the Sloan Digital Sky Survey Data Release 12, and measure the ISW amplitude to be $A_\mathrm{ISW}=1.64\pm0.53$ relative to the $\Lambda$CDM expectation, corresponding to a $3.1\sigma$ detection. In contrast to some previous measurements of the ISW effect of superstructures, our result is in agreement with the $\Lambda$CDM model.

The standard numerical tools for studying non-linear collapse of matter are Newtonian $N$-body simulations. Previous work has shown that these simulations are in accordance with General Relativity (GR) up to first order in perturbation theory, provided that the effects from radiation can be neglected. In this paper we show that the present day matter density receives more than 1$\%$ corrections from radiation on large scales if Newtonian simulations are initialised before $z=50$. We provide a relativistic framework in which \emph{unmodified} Newtonian simulations are compatible with linear GR even in the presence of radiation. Our idea is to use GR perturbation theory to keep track of the evolution of relativistic species and the relativistic space-time consistent with the Newtonian trajectories computed in $N$-body simulations. If metric potentials are sufficiently small, they can be computed using a first-order Einstein--Boltzmann code such as CLASS. We make this idea rigorous by defining a class of GR gauges, the \emph{Newtonian motion} gauges, which are defined such that matter particles follow Newtonian trajectories. We construct a simple example of a relativistic space-time within which unmodified Newtonian simulations can be interpreted.

Spatial correlations of the observed sizes and luminosities of galaxies can be used to estimate the magnification that arises through weak gravitational lensing. However, the intrinsic prop- erties of galaxies can be similarly correlated through local physical effects, and these present a possible contamination to the weak lensing estimation. In an earlier paper (Ciarlariello et al. 2015) we modelled the intrinsic size correlations using the halo model, assuming the galaxy sizes reflect the mass in the associated halo. Here we extend this work to consider galaxy magnitudes and show that these may be even more affected by intrinsic correlations than galaxy sizes, making this a bigger systematic for measurements of the weak lensing signal. We also quantify how these intrinsic correlations are affected by sample selection criteria based on sizes and magnitudes.

We present a fully relativistic calculation of the matter bispectrum at second order in cosmological perturbation theory assuming a Gaussian primordial curvature perturbation. For the first time we perform a full numerical integration of the bispectrum for both baryons and cold dark matter using the second-order Einstein-Boltzmann code, SONG. We review previous analytical results and provide an improved analytic approximation for the second-order kernel in Poisson gauge which incorporates Newtonian nonlinear evolution, relativistic initial conditions, the effect of radiation at early times and the cosmological constant at late times. Our improved kernel provides a percent level fit to the full numerical result at late times for most configurations, including both equilateral shapes and the squeezed limit. We show that baryon acoustic oscillations leave an imprint in the matter bispectrum, making a significant impact on squeezed shapes.

We explore the potential of using intensity mapping surveys (MeerKAT, SKA) and optical galaxy surveys (DES, LSST) to detect HI clustering and weak gravitational lensing of 21cm emission in auto- and cross-correlation. Our forecasts show that high precision measurements of the clustering and lensing signals can be made in the near future using the intensity mapping technique. Such studies can be used to test the intensity mapping method, and constrain parameters such as the HI density $\Omega_{\rm HI}$, the HI bias $b_{\rm HI}$ and the galaxy-HI correlation coefficient $r_{\rm HI-g}$.

We measure the cross-correlation between the galaxy density in the Dark Energy Survey (DES) Science Verification data and the lensing of the cosmic microwave background (CMB) as reconstructed with the Planck satellite and the South Pole Telescope (SPT). When using the DES main galaxy sample over the full redshift range $0.2 < z < 1.2$, a cross-correlation signal is detected at $6 \sigma$ and $4\sigma$ with SPT and Planck respectively. We then divide the DES galaxies into five photometric redshift bins, finding significant ($>$$2 \sigma$) detections in all bins. Comparing to the fiducial Planck cosmology, we find the redshift evolution of the signal matches expectations, although the amplitude is consistently lower than predicted across redshift bins. We test for possible systematics that could affect our result and find no evidence for significant contamination. Finally, we demonstrate how these measurements can be used to constrain the growth of structure across cosmic time. We find the data are fit by a model in which the amplitude of structure in the $z<1.2$ universe is $0.73 \pm 0.16$ times as large as predicted in the LCDM Planck cosmology, a $1.7\sigma$ deviation.

An interaction between the vacuum energy and dark matter is an intriguing possibility which may offer a way of solving the cosmological constant problem. Adopting a general prescription for momentum exchange between the two dark components, we reconstruct $\alpha(a)$, the temporal evolution of the coupling strength between dark matter and vacuum energy, in a nonparametric Bayesian approach using combined observational data sets from the cosmic microwave background, supernovae and large scale structure. An evolving interaction between the vacuum energy and dark matter removes some of the tensions between different data sets. However, it is not preferred over $\Lambda$CDM in the Bayesian sense, as improvement in the fit is not sufficient to compensate for the increase in the volume of the parameter space.

The initial conditions for Newtonian $N$-body simulations are usually generated by applying the Zel'dovich approximation to the initial displacements of the particles using an initial power spectrum of density fluctuations generated by an Einstein-Boltzmann solver. We show that in most gauges the initial displacements generated in this way receive a first-order relativistic correction. We define a new gauge, the $N$-body gauge, in which this relativistic correction vanishes and show that a conventional Newtonian $N$-body simulation includes all first-order relativistic contributions (in the absence of radiation) if we identify the coordinates in Newtonian simulations with those in the relativistic $N$-body gauge.

The cross-correlation of a foreground density field with two different background convergence fields can be used to measure cosmographic distance ratios and constrain dark energy parameters. We investigate the possibility of performing such measurements using a combination of optical galaxy surveys and HI intensity mapping surveys, with emphasis on the performance of the planned Square Kilometre Array (SKA). Using HI intensity mapping to probe the foreground density tracer field and/or the background source fields has the advantage of excellent redshift resolution and a longer lever arm achieved by using the lensing signal from high redshift background sources. Our results show that, for our best SKA-optical configuration of surveys, a constant equation of state for dark energy can be constrained to $\simeq 8\%$ for a sky coverage $f_{\rm sky}=0.5$ and assuming a $\sigma(\Omega_{\rm DE})=0.03$ prior for the dark energy density parameter. We also show that using the CMB as the second source plane is not competitive, even when considering a COrE-like satellite.

In this work we study the relevance of the cosmic web and substructures on the matter and lensing power spectra measured from halo mock catalogues extracted from the N-body simulations. Since N-body simulations are computationally expensive, it is common to use faster methods that approximate the dark matter field as a set of halos. In this approximation, we replace mass concentrations in N-body simulations by a spherically symmetric Navarro-Frenk-White halo density profile. We also consider the full mass field as the sum of two distinct fields: dark matter halos ($M>9\times 10^{12}~M_{\odot}$/h) and particles not included into halos. Mock halos reproduce well the matter power spectrum, but underestimate the lensing power spectrum on large and small scales. For sources at $z_{\rm s}=1$ the lensing power spectrum is underestimated by up to 40% at $\ell\approx 10^4$ with respect to the simulated halos. The large scale effect can be alleviated by combining the mock catalogue with the dark matter distribution outside the halos. In addition, to evaluate the contribution of substructures we have smeared out the intra-halo substructures in a N-body simulation while keeping the halo density profiles unchanged. For the matter power spectrum the effect of this smoothing is only of the order of 5%, but for lensing substructures are much more important: for $\ell\approx 10^4$ the internal structures contribute 30% of the total spectrum. These findings have important implications in the way mock catalogues have to be created, suggesting that some approximate methods currently used for galaxy surveys will be inadequate for future weak lensing surveys.

We present a simple model for describing intrinsic correlations for galaxy sizes based on the halo model. Studying these correlations is important both to improve our understanding of galaxy properties and because it is an important potential systematic for weak lensing size magnification measurements. Our model assumes that the density field drives these intrinsic correlations and we also model the distribution of satellite galaxies. We calculate the possible contamination to measurements of lensing convergence power spectrum from galaxy sizes, and show that the cross-correlation of intrinsic sizes with convergence is potentially an important systematic. We also explore how these intrinsic size correlations may affect surveys with different redshift depth. We find that, in this simple approach, intrinsic size correlations cannot be neglected in order to estimate lensing convergence power spectrum for constraining cosmological parameters.

Dark matter and dark energy are usually assumed to couple only gravitationally. An extension to this picture is to model dark energy as a scalar field coupled directly to cold dark matter. This coupling leads to new physical effects, such as a fifth-force and a time-dependent dark matter particle mass. In this work we examine the impact that coupling has on weak lensing statistics by constructing realistic simulated weak-lensing maps using raytracing techniques through N-body cosmological simulations. We construct maps for different lensing quantities, covering a range of scales from a few arcminutes to several degrees. The concordance $\Lambda$CDM model is compared to different coupled dark energy models, described either by an exponential scalar field potential (standard coupled dark energy scenario) or by a SUGRA potential (bouncing model). We analyse several statistical quantities and our results, with sources at low redshifts are largely consistent with previous work on CMB lensing by Carbone et al., 2013. The most significant differences from the $\Lambda$CDM model are due to the enhanced growth of the perturbations and to the effective friction term in non-linear dynamics. For the most extreme models, we see differences in the power spectra up to 40% compared to the $\Lambda$CDM model. The different time evolution of the linear matter overdensity can account for most of the differences, but when controlling for this using a $\Lambda$CDM model having the same normalization, the overall signal is smaller due to the effect of the friction term appearing in the equation of motion for dark matter particles.

We compute the bispectrum induced in the cosmic microwave background (CMB) temperature and polarisation by the evolution of the primordial density perturbations using the second-order Boltzmann code SONG. We show that adding polarisation increases the signal-to-noise ratio by a factor four with respect to temperature alone and we estimate the observability of this intrinsic bispectrum and the bias it induces on measurements of primordial non-Gaussianity. When including all physical effects except the late-time non-linear evolution, we find for the intrinsic bispectrum a signal-to-noise of $S/N=3.8,\,2.9,\,1.6$ and $0.5$ for, respectively, an ideal experiment with an angular resolution of $\ell_\text{max}=3000$, the proposed CMB surveys PRISM and COrE, and Planck's polarised data; the bulk of this signal comes from the $E$-polarisation and from squeezed configurations. We discuss how CMB lensing is expected to reduce these estimates as it suppresses the bispectrum for squeezed configurations and contributes to the noise in the estimator. We find that the presence of the intrinsic bispectrum will bias a measurement of primordial non-Gaussianity of local type by $f_{NL}^\text{bias}=0.66$ for an ideal experiment with $\ell_\text{max}=3000$. Finally, we verify the robustness of our results by reproducing the analytical approximation for the squeezed-limit bispectrum in the general polarised case.

We estimate the B-polarisation induced in the Cosmic Microwave Background by the non-linear evolution of density perturbations. Using the second-order Boltzmann code SONG, our analysis incorporates, for the first time, all physical effects at recombination. We also include novel contributions from the redshift part of the Boltzmann equation and from the bolometric definition of the temperature in the presence of polarisation. The remaining line-of-sight terms (lensing and time-delay) have previously been studied and must be calculated non-perturbatively. The intrinsic B-mode polarisation is present independent of the initial conditions and might contaminate the signal from primordial gravitational waves. We find this contamination to be comparable to a primordial tensor-to-scalar ratio of $r\simeq10^{-7}$ at the angular scale $\ell\simeq100\,$, where the primordial signal peaks, and $r\simeq 5 \cdot 10^{-5}$ at $\ell\simeq700\,$, where the intrinsic signal peaks. Therefore, we conclude that the intrinsic B-polarisation from second-order effects is not likely to contaminate future searches of primordial gravitational waves.

We present the strongest robust constraints on primordial non-Gaussianity (PNG) from currently available galaxy surveys, combining large-scale clustering measurements and their cross-correlations with the cosmic microwave background. We update the data sets used by Giannantonio et al. (2012), and broaden that analysis to include the full set of two-point correlation functions between all surveys. In order to obtain the most reliable constraints on PNG, we advocate the use of the cross-correlations between the catalogs as a robust estimator and we perform an extended analysis of the possible systematics to reduce their impact on the results. To minimize the impact of stellar contamination in our luminous red galaxy (LRG) sample, we use the recent Baryon Oscillations Spectroscopic Survey catalog of Ross et al. (2011). We also find evidence for a new systematic in the NVSS radio galaxy survey similar to, but smaller than, the known declination-dependent issue; this is difficult to remove without affecting the inferred PNG signal, and thus we do not include the NVSS auto-correlation function in our analyses. We find no evidence of primordial non-Gaussianity; for the local-type configuration we obtain for the skewness parameter $ -36 < f_{\mathrm{NL}} < 45 $ at 95 % c.l. ($5 \pm 21$ at $1\sigma$) when using the most conservative part of our data set, improving previous results; we also find no evidence for significant kurtosis, parameterized by $g_{\mathrm{NL}}$. In addition to PNG, we simultaneously constrain dark energy and find that it is required with a form consistent with a cosmological constant.

We present improved methodology for including covariance matrices in the error budget of Baryon Oscillation Spectroscopic Survey (BOSS) galaxy clustering measurements, revisiting Data Release 9 (DR9) analyses, and describing a method that is used in DR10/11 analyses presented in companion papers. The precise analysis method adopted is becoming increasingly important, due to the precision that BOSS can now reach: even using as many as 600 mock catalogues to estimate covariance of 2-point clustering measurements can still lead to an increase in the errors of ~20%, depending on how the cosmological parameters of interest are measured. In this paper we extend previous work on this contribution to the error budget, deriving formulae for errors measured by integrating over the likelihood, and to the distribution of recovered best-fit parameters fitting the simulations also used to estimate the covariance matrix. Both are situations that previous analyses of BOSS have considered. We apply the formulae derived to Baryon Acoustic Oscillation (BAO) and Redshift-Space Distortion (RSD) measurements from BOSS in our companion papers. To further aid these analyses, we consider the optimum number of bins to use for 2-point measurements using the monopole power spectrum or correlation function for BAO, and the monopole and quadrupole moments of the correlation function for anisotropic-BAO and RSD measurements.

We study structure formation in non-minimally coupled dark energy models, where there is a coupling in the Lagrangian between a quintessence scalar field and gravity via the Ricci scalar. We consider models with a range of different non-minimal coupling strengths and compare these to minimally coupled quintessence models with time-dependent dark energy densities. The equations of state of the latter are tuned to either reproduce the equation of state of the non-minimally coupled models or their background history. Thereby they provide a reference to study the unique imprints of coupling on structure formation. We show that the coupling between gravity and the scalar field, which effectively results in a time-varying gravitational constant G, is not negligible and its effect can be distinguished from a minimally coupled model. We extend previous work on this subject by showing that major differences appear in the determination of the mass function at high masses, where we observe differences of the order of 40% at z=0. Our new results concern effects on the non-linear matter power spectrum and on the lensing signal (differences of ~10% for both quantities), where we find that non-minimally coupled models could be distinguished from minimally coupled ones.

We develop a new, efficient code for solving the second-order Einstein-Boltzmann equations, and use it to estimate the intrinsic CMB non-Gaussianity arising from the non-linear evolution of density perturbations. The full calculation involves contributions from recombination and less tractable contributions from terms integrated along the line of sight. We investigate the bias that this intrinsic bispectrum implies for searches of primordial non-Gaussianity. We find that the inclusion or omission of certain line of sight terms can make a large impact. When including all physical effects but lensing and time-delay, we find that the local-type f_nl would be biased by f_nl ~ 0.5, below the expected sensitivity of the Planck satellite. The speed of our code allows us to confirm the robustness of our results with respect to a number of numerical parameters.

We apply a new non-parametric Bayesian method for reconstructing the evolution history of the equation-of-state $w$ of dark energy, based on applying a correlated prior for $w(z)$, to a collection of cosmological data. We combine the latest supernova (SNLS 3-year or Union2.1), cosmic microwave background, redshift space distortion and the baryonic acoustic oscillation measurements (including BOSS, WiggleZ and 6dF) and find that the cosmological constant appears consistent with current data, but that a dynamical dark energy model which evolves from $w<-1$ at $z\sim0.25$ to $w > -1$ at higher redshift is mildly favored. Estimates of the Bayesian evidences show little preference between the cosmological constant model and the dynamical model for a range of correlated prior choices. Looking towards future data, we find that the best fit models for current data could be well distinguished from the $\Lambda$CDM model by observations such as Planck and Euclid-like surveys.

We revisit the state of the integrated Sachs-Wolfe (ISW) effect measurements in light of newly available data and address criticisms about the measurements which have recently been raised. We update the data set previously assembled by Giannantonio et al. to include new data releases for both the cosmic microwave background (CMB) and the large-scale structure (LSS) of the Universe. We find that our updated results are consistent with previous measurements. By fitting a single template amplitude, we now obtain a combined significance of the ISW detection at the 4.4 sigma level, which fluctuates by 0.4 sigma when alternative data cuts and analysis assumptions are considered. We also make new tests for systematic contaminations of the data, focusing in particular on the issues raised by Sawangwit et al. Amongst them, we address the rotation test, which aims at checking for possible systematics by correlating pairs of randomly rotated maps. We find results consistent with the expected data covariance, no evidence for enhanced correlation on any preferred axis of rotation, and therefore no indication of any additional systematic contamination. We publicly release the results, the covariance matrix, and the sky maps used to obtain them.

Assuming General Relativity is correct on large-scales, Redshift-Space Distortions (RSDs) and the Integrated Sachs-Wolfe effect (ISW) are both sensitive to the time derivative of the linear growth function. We investigate the extent to which these probes provide complementary or redundant information when they are combined to constrain the evolution of the linear velocity power spectrum, often quantified by the function $f(z)\sigma_8(z)$, where $f$ is the logarithmic derivative of $\sigma_8$ with respect to $(1+z)$. Using a spherical Fourier-Bessel (SFB) expansion for galaxy number counts and a spherical harmonic expansion for the CMB anisotropy, we compute the covariance matrices of the signals for a large galaxy redshift survey combined with a CMB survey like Planck. The SFB basis allows accurate ISW estimates by avoiding the plane-parallel approximation, and it retains RSD information that is otherwise lost when projecting angular clustering onto redshift shells. It also allows straightforward calculations of covariance with the CMB. We find that the correlation between the ISW and RSD signals are low since the probes are sensitive to different modes. For our default surveys, on large scales ($k<0.05 \Mpc/h$), the ISW can improve constraints on $f\sigma_8$ by more than 10% compared to using RSDs alone. In the future, when precision RSD measurements are available on smaller scales, the cosmological constraints from ISW measurements will not be competitive; however, they will remain a useful consistency test for possible systematic contamination and alternative models of gravity.

We use a Quadratic Maximum Likelihood (QML) method to estimate the angular power spectrum of the cross-correlation between cosmic microwave background and large scale structure maps as well as their individual auto-spectra. We describe our implementation of this method and demonstrate its accuracy on simulated maps. We apply this optimal estimator to WMAP 7-year and NRAO VLA Sky Survey (NVSS) data and explore the robustness of the angular power spectrum estimates obtained by the QML method. With the correction of the declination systematics in NVSS, we can safely use most of the information contained in this survey. We then make use of the angular power spectrum estimates obtained by the QML method to derive constraints on the dark energy critical density in a flat $\Lambda$CDM model by different likelihood prescriptions. When using just the cross-correlation between WMAP 7 year and NVSS maps with 1.8$^\circ$ resolution, the best-fit model has a cosmological constant of approximatively 70% of the total energy density, disfavouring an Einstein-de Sitter Universe at more than 2 $\sigma$ CL (confidence level).

Local non-Gaussianity, parametrized by $f_{\rm NL}$, introduces a scale-dependent bias that is strongest at large scales, precisely where General Relativistic (GR) effects also become significant. With future data, it should be possible to constrain $f_{\rm NL} = {\cal O}(1)$ with high redshift surveys. GR corrections to the power spectrum and ambiguities in the gauge used to define bias introduce effects similar to $f_{\rm NL}= {\cal O}(1)$, so it is essential to disentangle these effects. For the first time in studies of primordial non-Gaussianity, we include the consistent GR calculation of galaxy power spectra, highlighting the importance of a proper definition of bias. We present observable power spectra with and without GR corrections, showing that an incorrect definition of bias can mimic non-Gaussianity. However, these effects can be distinguished by their different redshift and scale dependence, so as to extract the true primordial non-Gaussianity.

The next generation of weak lensing surveys will trace the evolution of matter perturbations and gravitational potentials from the matter dominated epoch until today. Along with constraining the dynamics of dark energy, they will probe the relations between matter overdensities, local curvature, and the Newtonian potential. We work with two functions of time and scale to account for any modifications of these relations in the linear regime from those in the LCDM model. We perform a Principal Component Analysis (PCA) to find the eigenmodes and eigenvalues of these functions for surveys like DES and LSST. This paper builds on and significantly extends the PCA analysis of Zhao et al. (2009) in several ways. In particular, we consider the impact of some of the systematic effects expected in weak lensing surveys. We also present the PCA in terms of other choices of the two functions needed to parameterize modified growth on linear scales, and discuss their merits. We analyze the degeneracy between the modified growth functions and other cosmological parameters, paying special attention to the effective equation of state w(z). Finally, we demonstrate the utility of the PCA as an efficient data compression stage which enables one to easily derive constraints on parameters of specific models without recalculating Fisher matrices from scratch.

We develop an efficient, non-parametric Bayesian method for reconstructing the time evolution of the dark energy equation of state w(z) from observational data. Of particular importance is the choice of prior, which must be chosen carefully to minimise variance and bias in the reconstruction. Using a principal component analysis, we show how a correlated prior can be used to create a smooth reconstruction and also avoid bias in the mean behaviour of w(z). We test our method using Wiener reconstructions based on Fisher matrix projections, and also against more realistic MCMC analyses of simulated data sets for Planck and a future space-based dark energy mission. While the accuracy of our reconstruction depends on the smoothness of the assumed w(z), the relative error for typical dark energy models is <10% out to redshift z=1.5.

We investigate the effects of potential sources of systematic error on the angular and photometric redshift, z_phot, distributions of a sample of redshift 0.4 < z < 0.7 massive galaxies whose selection matches that of the Baryon Oscillation Spectroscopic Survey (BOSS) constant mass sample. Utilizing over 112,778 BOSS spectra as a training sample, we produce a photometric redshift catalog for the galaxies in the SDSS DR8 imaging area that, after masking, covers nearly one quarter of the sky (9,913 square degrees). We investigate fluctuations in the number density of objects in this sample as a function of Galactic extinction, seeing, stellar density, sky background, airmass, photometric offset, and North/South Galactic hemisphere. We find that the presence of stars of comparable magnitudes to our galaxies (which are not traditionally masked) effectively remove area. Failing to correct for such stars can produce systematic errors on the measured angular auto-correlation function, w, that are larger than its statistical uncertainty. We describe how one can effectively mask for the presence of the stars, without removing any galaxies from the sample, and minimize the systematic error. Additionally, we apply two separate methods that can be used to correct the systematic errors imparted by any parameter that can be turned into a map on the sky. We find that failing to properly account for varying sky background introduces a systematic error on w. We measure w, in four z_phot slices of width 0.05 between 0.45 < z_phot < 0.65 and find that the measurements, after correcting for the systematic effects of stars and sky background, are generally consistent with a generic LambdaCDM model, at scales up to 60 degrees. At scales greater than 3 degrees and z_phot > 0.5, the magnitude of the corrections we apply are greater than the statistical uncertainty in w.

Burrage, Davis, and Shaw recently suggested exploiting the correlations between high and low energy luminosities of astrophysical objects to probe possible mixing between photons and axion-like particles (ALP) in magnetic field regions. They also presented evidence for the existence of ALP's by analyzing the optical/UV and X-ray monochromatic luminosities of AGNs. We extend their work by using the monochromatic luminosities of 320 unobscured Active Galactic Nuclei from the Sloan Digital Sky Survey/Xmm-Newton Quasar Survey (Young et al., 2009), which allows the exploration of 18 different combinations of optical/UV and X-ray monochromatic luminosities. However, we do not find compelling evidence for the existence of ALPs. Moreover, it appears that the signal reported by Burrage et al. is more likely due to X-ray absorption rather than to photon-ALP oscillation.

We apply the needlet formalism to the Wilkinson Microwave Anisotropy Probe 5-year data, looking for evidence of non-Gaussianity in the bispectrum of the needlet amplitudes. We confirm earlier findings of an asymmetry in the non-Gaussianity between the northern and southern galactic hemispheres. We attempt to isolate which scales and geometrical configurations are most anomalous, and find the bispectrum is most significant on large scales and in the more co-linear configurations, and also in the `squeezed' configurations. However, these anomalies do not appear to affect the estimate of the non-linear parameter $\fnl$, and we see no significant difference between its value measured in the two hemispheres.

We use a principal component approach to contrast different kinds of probes of dark energy, and to emphasize how an array of probes can work together to constrain an arbitrary equation of state history w(z). We pay particular attention to the role of the priors in assessing the information content of experiments and propose using an explicit prior on the degree of smoothness of w(z) that is independent of the binning scheme. We also show how a figure of merit based on the mean squared error probes the number of new modes constrained by a data set, and use it to examine how informative various experiments will be in constraining the evolution of dark energy.

We constrain the amplitude of primordial non-Gaussianity in the CMB data taking into account the presence of foreground residuals in the maps. We generalise the needlet bispectrum estimator marginalizing over the amplitudes of thermal dust, free-free and synchrotron templates. We apply our procedure to WMAP 5 year data, finding fNL= 38\pm 47 (1 \sigma), while the analysis without marginalization provides fNL= 35\pm 42. Splitting the marginalization over each foreground separately, we found that the estimates of fNL are positively cross correlated of 17%, 12% with the dust and synchrotron respectively, while a negative cross correlation of about -10% is found for the free-free component.

Dark energy can be studied by its influence on the expansion of the Universe as well as on the growth history of the large-scale structure. In this paper, we follow the growth of the cosmic density field in early dark energy cosmologies by combining observations of the primary CMB temperature and polarisation power spectra at high redshift, of the CMB lensing deflection field at intermediate redshift and of weak cosmic shear at low redshifts for constraining the allowed amount of early dark energy. We present these forecasts using the Fisher-matrix formalism and consider the combination of Planck-data with the weak lensing survey of Euclid. We find that combining these data sets gives powerful constraints on early dark energy and is able to break degeneracies in the parameter set inherent to the various observational channels. The derived statistical 1-sigma-bound on the early dark energy density parameter is sigma(Omega_d^e)=0.0022 which suggests that early dark energy models can be well examined in our approach. In addition, we derive the dark energy figure of merit for the considered dark energy parameterisation and comment on the applicability of the growth index to early dark energy cosmologies.

We present a global measurement of the integrated Sachs-Wolfe (ISW) effect obtained by cross-correlating all relevant large scale galaxy data sets with the cosmic microwave background radiation map provided by the Wilkinson Microwave Anisotropy Probe. With these measurements, the overall ISW signal is detected at the ~ 4.5 sigma level. We also examine the cosmological implications of these measurements, particularly the dark energy equation of state w, its sound speed, and the overall curvature of the Universe. The flat LCDM model is a good fit to the data and, assuming this model, we find that the ISW data constrain Omega_m = 0.20 +0.19 -0.11 at the 95% confidence level. When we combine our ISW results with the latest baryon oscillation and supernovae measurements, we find that the result is still consistent with a flat LCDM model with w = -1 out to redshifts z > 1.

We analyse the possibility that topological defects can act as a source of magnetic fields through the Harrison mechanism in the radiation era. We give a detailed relativistic derivation of the Harrison mechanism at first order in cosmological perturbations, and show that it is only efficient for temperatures above T ~ 0.2 keV. Our main result is that the vector metric perturbations generated by the defects cannot induce vorticity in the matter fluids at linear order, thereby excluding the production of currents and magnetic fields. We show that anisotropic stress in the matter fluids is required to source vorticity and magnetic fields. Our analysis is relevant for any mechanism whereby vorticity is meant to be transferred purely by gravitational interactions, and thus would also apply to dark matter or neutrinos.

Most parameterizations of the dark energy equation of state do not reflect realistic underlying physical models. Here, we develop a relatively simple description of dark energy based on the dynamics of a scalar field which is exact in the limit that the equation of state approaches a cosmological constant, assuming some degree of smoothness of the potential. By introducing just two parameters defined in the configuration space of the field we are able to reproduce a wide class of quintessence models. We examine the observational constraints on these models as compared to linear evolution models, and show how priors in the field space translate into priors on observational parameters.

In this paper we show how the rescattering of CMB photons after cosmic reionization can give a significant linear contribution to the temperature-matter cross-correlation measurements. These anisotropies, which arise via a late time Doppler effect, are on scales much larger than the typical scale of non-linear effects at reionization; they can contribute to degree scale cross-correlations and could affect the interpretation of similar correlations resulting from the integrated Sachs-Wolfe effect. While expected to be small at low redshifts, these correlations can be large given a probe of the density at high redshift, and so could be a useful probe of the cosmic reionization history.

We present evidence of a large angle correlation between the cosmic microwave background measured by WMAP and a catalog of photometrically detected quasars from the SDSS. The observed cross correlation is (0.30 +- 0.14) microK at zero lag, with a shape consistent with that expected for correlations arising from the integrated Sachs-Wolfe effect. The photometric redshifts of the quasars are centered at z ~ 1.5, making this the deepest survey in which such a correlation has been observed. Assuming this correlation is due to the ISW effect, this constitutes the earliest evidence yet for dark energy and it can be used to constrain exotic dark energy models.

We present an analysis of the constraining power of future measurements of the Integrated Sachs-Wolfe (ISW) effect on models of the equation of state of dark energy as a function of redshift, w(z). To achieve this, we employ a new parameterization of w, which utilizes the mean value of w(z) (<w>) as an explicit parameter. This helps to separate the information contained in the estimation of the distance to the last scattering surface (from the CMB) from the information contained in the ISW effect. We then use Fisher analysis to forecast the expected uncertainties in the measured parameters from future ISW observations for two models of dark energy with very different time evolution properties. For example, we demonstrate that the cross-correlation of Planck CMB data and LSST galaxy catalogs will provide competitive constraints on w(z), compared to a SNAP-like SNe project, for models of dark energy with a rapidly changing equation of state (e.g. 'Kink' models). Our work confirms that, while SNe measurements are more suitable for constraining variations in w(z) at low redshift, the ISW effect can provide important independent constraints on w(z) at high z.

Magnetic fields in the early universe could have played an important role in sourcing cosmological perturbations. While not the dominant source, even a small contribution might be traceable through its intrinsic non-Gaussianity. Here we calculate analytically the one, two and three point statistics of the magnetic stress energy resulting from tangled Gaussian fields, and confirm these with numerical realizations of the fields. We find significant non-Gaussianity, and importantly predict higher order moments that will appear between the scalar, vector and tensor parts of the stress energy (e.g. scalar-tensor-tensor moments). Such higher order cross correlations are a generic feature of non-linear theories and could prove to be an important probe of the early universe.

Recently, we presented evidence for a cross-correlation of the WMAP satellite map of the cosmic microwave background (CMB) and the HEAO1 satellite map of the hard X-ray background (XRB) with a dimensionless amplitude of 0.14 +- 0.05 normalized to the product of the rms fluctuations of the CMB and XRB (Boughn & Crittenden, 2004). Such a correlation is expected in a universe dominated by a cosmological constant via the integrated Sachs-Wolfe (ISW) effect and the level of the correlation observed is consistent with that predicted by the currently favored Lambda cold dark matter model of the universe. Since this offers independent confirmation of the cosmological model, it is important to verify the origin of the correlation. Here we explore in detail some possible foreground sources of the correlation. The present evidence all supports an ISW origin.

Recent deep X-ray surveys combined with spectroscopic identification of the sources have allowed the determination of the rest-frame 2-8 keV luminosity as a function of redshift. In addition, an analysis of the HEAO1 A2 2-10 keV full-sky map of the X-ray background (XRB) reveals clustering on the scale of several degrees. Combining these two results in the context of the currently favored Lambda-CDM cosmological model implies an average X-ray bias factor, b_x, of b_x^2 = 1.12 +- 0.33, i.e., b_x = 1.06 +- 0.16. These error estimates include only statistical error; the systematic error sources, while comparable, appear to be sub-dominant. This result is in contrast to the large biases of some previous estimates and is more in line with current estimates of the optical bias of L* galaxies.

We have detected statistically significant correlations between the cosmic microwave background and two tracers of large-scale structure, the HEAO1 A2 full sky hard X-ray map and the NVSS 1.4 GHz, nearly full sky radio galaxy survey. The level of correlations in these maps is consistent with that predicted for the integrated Sachs-Wolfe (ISW) effect in the context of a Lambda CDM cosmological model and, therefore, provides independent evidence for a cosmological constant. A maximum likelihood fit to the amplitude of the ISW effect relative to the predicted value is 1.13 +- 0.35 (statistical error only).

We cross correlate the large-scale cosmic microwave background (CMB) sky measured by WMAP with two probes of large-scale structure at z ~ 1. The hard X-ray background, measured by the HEAO-1 satellite, is positively correlated with the WMAP data at the 2.5-3.0 sigma level. The number counts of radio galaxies in the NVSS survey are also correlated at a slightly weaker level (2.-2.5 sigma). These correlations appear to arise from both hemispheres on the sky and are resilient to changes in the levels of masking of the Galaxy and point sources, suggesting that foregrounds are not responsible for the signal. The implication is that some of the observed CMB fluctuations arise at low redshifts. The level of the correlations is consistent with that expected for the cosmological constant (Omega_Lambda = 0.72) concordance model resulting from the integrated Sachs-Wolfe effect. Thus, we may be observing dark energy's effect on the growth of structure.

With the increased accuracy and angular scale coverage of the recent CMB experiments it has become important to include calibration and beam uncertainties when estimating cosmological parameters. This requires an integration over possible values of the calibration and beam size, which can be done numerically but increases computation times. We present a fast and simple algorithm for marginalization over beam and calibration errors by analytical integration. We also illustrate the effect of incorporating these uncertainties by calculating the constraints on various cosmological and inflationary parameters including the spectral index n_s and the physical baryon density Omega_b h^2, using the latest CMB data. We find that parameter constraints are significantly changed when calibration/beam uncertainties are taken into account. Typically the best fit parameters are shifted and the errors bars are increased by up to fifty per cent for e.g. n_s and Omega_b h^2, although as expected there is no change for Omega_K, because it is constrained by the positions of the peaks.

A careful analysis of the HEAO1 A2 2-10 keV full-sky map of the X-ray background (XRB) reveals clustering on the scale of several degrees. After removing the contribution due to beam smearing, the intrinsic clustering of the background is found to be consistent with an auto-correlation function of the form (3.6 +- 0.9) x 10^{-4} theta^{-1} where theta is measured in degrees. If current AGN models of the hard XRB are reasonable and the cosmological constant-cold dark matter cosmology is correct, this clustering implies an X-ray bias factor of b_X ~ 2. Combined with the absence of a correlation between the XRB and the cosmic microwave background, this clustering can be used to limit the presence of an integrated Sachs-Wolfe (ISW) effect and thereby to constrain the value of the cosmological constant, Omega_Lambda < 0.60 (95 % C.L.). This constraint is inconsistent with much of the parameter space currently favored by other observations. Finally, we marginally detect the dipole moment of the diffuse XRB and find it to be consistent with the dipole due to our motion with respect to the mean rest frame of the XRB. The limit on the amplitude of any intrinsic dipole is delta I / I < 5 x 10^{-3} at the 95 % C.L. When compared to the local bulk velocity, this limit implies a constraint on the matter density of the universe of Omega_m^{0.6}/b_X(0) > 0.24.

We present a new limit on the cosmological constant based on the absence of correlations between the cosmic microwave background (CMB) and the distribution of distant radio sources. In the cosmological constant-cold dark matter models currently favored, such correlations should have been produced via the integrated Sachs-Wolfe effect, assuming that radio sources trace the local (z=1) matter density. We find no evidence of correlations between the COBE 53Hz microwave map and the NVSS 1.4 GHz radio survey. The implied 95% CL limit on the cosmological constant is Lambda < 0.74, in marginal agreement with the values suggested by recent measurements of the CMB anisotropy and type-IA supernovae observations, 0.6 < Lambda < 0.7. If the cosmological model does lie in this range, then the integrated Sachs-Wolfe effect should be detectable with upcoming CMB maps and radio surveys.

We describe the subject of Cosmic Microwave Background (CMB) analysis - its past, present and future. The theory of Gaussian primary anisotropies, those arising from linear physics operating in the early Universe, is in reasonably good shape so the focus has shifted to the statistical pipeline which confronts the data with the theory: mapping, filtering, comparing, cleaning, compressing, forecasting, estimating. There have been many algorithmic advances in the analysis pipeline in recent years, but still more are needed for the forecasts of high precision cosmic parameter estimation to be realized. For secondary anisotropies, those arising once nonlinearity develops, the computational state of the art currently needs effort in all the areas: the Sunyaev-Zeldovich effect, inhomogeneous reionization, gravitational lensing, the Rees-Sciama effect, dusty galaxies. We use the Sunyaev-Zeldovich example to illustrate the issues. The direct interface with observations for these non-Gaussian signals is much more complex than for Gaussian primary anisotropies, and even more so for the statistically inhomogeneous Galactic foregrounds. Because all the signals are superimposed, the separation of components inevitably complicates primary CMB analyses as well.

Alignments in the angular momentum vectors of galaxies can induce large scale correlations in their projected orientations. Such alignments arise from the tidal torques exerted on neighboring proto-galaxies by the smoothly varying shear field. Weak gravitational lensing can also induce ellipticity correlations since the images of neighboring galaxies will be distorted coherently by the intervening mass distribution. Comparing these two sources of shape correlations, it is found that for current weak lensing surveys with a median redshift of 1, the intrinsic signal is a contaminant on the order of 1-10% of the measured signal. However, for shallower surveys, the intrinsic correlations dominate over the lensing signal. The distortions induced by lensing are curl-free, whereas those resulting from intrinsic alignments are not. This difference can be used to disentangle these two sources of ellipticity correlations. When the distortions are dominated by lensing, as occurs at high redshifts, the decomposition provides a valuable tool for understanding properties of the noise and systematic errors.

The distortion field defined by the ellipticities of galaxy shapes projected on the sky can be uniquely decomposed into a gradient and a curl component. If the observed ellipticities are induced by weak gravitational lensing, then the distortion field is curl free. Here we show that, in contrast, the distortion field resulting from intrinsic spin alignments is not curl free. This provides a powerful discriminant between lensing and intrinsic contributions to observed ellipticity correlations. We also show how these contributions can be disentangled statistically from the ellipticity correlations or computed locally from circular integrals of the ellipticity field. This allows for an unambiguous detection of intrinsic galaxy alignments in the data. When the distortions are dominated by lensing, as occurs at high redshifts, the decomposition provides a valuable tool for understanding properties of the noise and systematic errors. These techniques can be applied equally well to the polarization of the microwave background, where it can be used to separate curl-free scalar perturbations from those produced by gravity waves or defects.

Large scale correlations in the orientations of galaxies can result from alignments in their angular momentum vectors. These alignments arise from the tidal torques exerted on neighboring proto-galaxies by the smoothly varying shear field. We compute the predicted amplitude of such ellipticity correlations using the Zel'dovich approximation for a realistic distribution of galaxy shapes. Weak gravitational lensing can also induce ellipticity correlations since the images of neighboring galaxies will be distorted coherently. On comparing these two effects that induce shape correlations, we find that for current weak lensing surveys with a median redshift of z_m = 1, the intrinsic signal is of order 1 - 10 percent of the measured signal. However, for shallower surveys with z_m < 0.3, the intrinsic correlations dominate over the lensing signal. The distortions induced by lensing are curl-free, whereas those resulting from intrinsic alignments are not. This difference can be used to disentangle these two sources of ellipticity correlations.

Recent measurements of a peak in the angular power spectrum of the cosmic microwave background appear to suggest that geometry of the universe is close to being flat. But if other accepted indicators of cosmological parameters are also correct then the best fit model is marginally closed, with the peak in the spectrum at larger scales than in a flat universe. Such observations can be reconciled with a flat universe if the fine structure constant had a lower value at earlier times, which would delay the recombination of electrons and protons and also act to suppress secondary oscillations as observed. We discuss evidence for a few percent increase in the fine structure constant between the time of recombination and the present.

The Cosmic Microwave Background (CMB) encodes information on the origin and evolution of the universe, buried in a fractional anisotropy of one part in 100,000 on angular scales from arcminutes to tens of degrees. We await the coming onslaught of data from experiments measuring the microwave sky from the ground, from balloons and from space. However, we are faced with the harsh reality that current algorithms for extracting cosmological information cannot handle data sets of the size and complexity expected even in the next few years. Here we review the challenges involved in understanding this data: making maps from time-ordered data, removing the foreground contaminants, and finally estimating the power spectrum and cosmological parameters from the CMB map. If handled naively, the global nature of the analysis problem renders these tasks effectively impossible given the volume of the data. We discuss possible techniques for overcoming these issues and outline the many other challenges that wait to be addressed.

Upcoming microwave background experiments will see an incredible increase in the volume of data to be analyzed, which makes the choice of how it is discretized on the sky a crucial issue. I discuss criteria for evaluating different pixelizations and the advantages of using an exactly azimuthal or `igloo' pixelization of the sky. Talk given at ``The CMB and the Planck Mission" workshop in Santander, Spain, June 1998.

We investigate various pixelizations of the sky which allow for fast spherical transforms, for implementation in full sky CMB experiments such as Planck and MAP. We study the effect of varying pixel shape and area on the extraction of the CMB power spectrum. We argue for the benefits of having a truly azimuthal, or `igloo' pixelization. Such pixelizations are simple and allow for fast, exact simulations of pixelized skies. They also allow for precise correction to be made which accounts for the effects of pixel smoothing on extracted multipole moments.

In universes with significant curvature or cosmological constant, cosmic microwave background (CMB) anisotropies are created very recently via the Rees-Sciama or integrated Sachs-Wolfe effects. This causes the CMB anisotropies to become partially correlated with the local matter density (z < 4). We examine the prospects of using the hard (2-10 keV) X-ray background as a probe of the local density and the measured correlation between the HEAO1 A2 X-ray survey and the 4-year COBE-DMR map to obtain a constraint on the cosmological constant. The 95% confidence level upper limit on the cosmological constant is \Omega_\Lambda \leq 0.5, assuming that the observed fluctuations in the X-ray map result entirely from large scale structure. (This would also imply that the X-rays trace matter with a bias factor of b_x = 5.6 \Omega_m^{0.53}.) This bound is weakened considerably if a large portion of the X-ray fluctuations arise from Poisson noise from unresolved sources. For example, if one assumes that the X-ray bias is b_x = 2., then the 95% confidence level upper limit is weaker, \Omega_\Lambda \leq 0.7. More stringent limits should be attainable with data from the next generation of CMB and X-ray background maps.

In models with a cosmological constant, a significant component of the large scale cosmic microwave background (CMB) anisotropy is produced at rather low redshifts, z < 1. In these models, the gravitational potential perturbations begin to evolve at late times. Photons passing through these time varying potentials aquire frequency shifts, an effect first noted by Rees and Sciama. Since the gravitational potential is associated with matter at observable redshifts, the local matter density (or some tracer thereof) is correlated with the CMB anisotropy on the sky. We outline the optimal method of measuring this correlation, and discuss the prospects for using an X-ray/COBE comparison to detect a cosmological constant.

We compute the angular power spectrum of temperature anisotropies on the microwave sky in the cosmic texture theory, with standard recombination assumed. The spectrum shows `Doppler' peaks analogous to those in scenarios based on primordial adiabatic fluctuations such as `standard CDM', but at quite different angular scales. There appear to be excellent prospects for using this as a discriminant between inflationary and cosmic defect theories.

We study the polarization-temperature correlations on the cosmic microwave sky resulting from an initial scale invariant spectrum of tensor (gravity wave) fluctuations, such as those which might arise during inflation. The correlation function has the opposite sign to that for scalar fluctuations on large scales, raising the possibility of a direct determination of whether the microwave anisotropies have a significant tensor component. We briefly discuss the important problem of estimating the expected foreground contamination.

The polarization of the CMB offers a potentially valuable probe of the physics of the early universe, specifically of the velocity of the photon-baryon fluid on the last scattering surface. Here we discuss qualitatively the correlations that arise between the temperature anisotropies and polarization of the CMB and how they might be observed.

We study the polarization-polarization and polarization-temperature correlations in standard adiabatic scenarios for structure formation. Temperature anisotropies due to gravitational potential wells and oscillations in the photon-baryon-electron fluid on the surface of last scattering are each associated with a correlated polarization pattern. While the `correlated part' of the polarization has an r.m.s. of only a third of the total signal, it may still be measurable by mapping a large area on the sky. We calculate the expected signal to noise ratio for various measures of the polarization in a hypothetical mapping experiment such as those now being planned.

The cosmic microwave background anisotropy is sensitive to the slope and amplitude of primordial energy density and gravitational wave fluctuations, the baryon density, the Hubble constant, the cosmological constant, the ionization history, {\it etc.} In this Letter, we examine the degree to which these factors can be separately resolved from combined small- and large-angular scale anisotropy observations. We isolate directions of degeneracy in this cosmic parameter space, but note that other cosmic observations can break the degeneracy. (Correction to Eq. 4)

The contribution of gravitational wave (tensor metric) and energy density (scalar metric) fluctuations to the cosmic microwave background polarization is computed by numerically solving the relativistic radiation transfer equations. We find that the tensor contribution is significant only at large angular scales (multipoles $\ell \lta 40$). For standard recombination, the tensor contribution can dominate at $\ell \lta 40$; however, the effect would be difficult to detect since the total (scalar plus tensor) polarization is $< 1$\%. For models with late reionization, the total large angular scale polarization is large ($\sim 7-9$\%), but the tensor fraction is negligibly small. Hence, polarization may be useful for discriminating ionization history, but is much less promising as a means for detecting tensor fluctuations.

Long-wavelength gravitational waves can induce significant temperature anisotropy in the cosmic microwave background. Distinguishing this from anisotropy induced by energy density fluctuations is critical for testing inflationary cosmology and theories of large-scale structure formation. We describe full radiative transport calculations of the two contributions and show that they differ dramatically at angular scales below a few degrees. We show how anisotropy experiments probing large- and small-angular scales can combine to distinguish the imprint due to gravitational waves.

Two qualitatively different modes of ending superluminal expansion are possible in extended inflation. One mode, different from the one envoked in most extended models to date, easily avoids making big bubbles that distort the cosmic microwave background radiation (CMBR). In this mode, the spectrum of density fluctuations is found to be scale-free, $P(k) \propto k^n$, where $n$ might lie anywhere between 0.5 and 1.0 (whereas, previously, it appeared that the range $1.0> n \gtsim 0.84$ was disallowed).