Source author record

Antony Lewis

Antony Lewis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2025arXiv

A flexible parameterization to test early physics solutions to the Hubble tension with future CMB data

One approach to reconciling local measurements of a high expansion rate with observations of acoustic oscillations in the CMB and galaxy clustering (the "Hubble tension") is to introduce additional contributions to the $Λ$CDM model that are relevant before recombination. While numerous possibilities exist, none are currently well-motivated or preferred by data. However, future CMB experiments, which will measure acoustic peaks to much smaller scales and resolve polarization signals with higher signal-to-noise over large sky areas, should detect almost any such modification at high significance. We propose a model-agnostic method to capture most relevant possible deviations from $Λ$CDM due to additional non-interacting components, while remaining sufficiently constraining to enable detection across various scenarios. The phenomenological model uses a fluid model with four parameters governing additional density contributions that peak at different redshifts, and two sound speed parameters. We forecast possible constraints with Simons Observatory, explore parameter degeneracies that arise in $Λ$CDM, and demonstrate that this method could detect a range of specific models. Which of the new parameters gets excited can indicate the nature of any new physics, while the generality of the model allows for testing with future data in a way that should not be plagued by a posteriori choices or publication bias. When testing our model with Planck data, we find good consistency with the $Λ$CDM model, but the data also allows for large Hubble parameter, especially if the sound speed of an additional component is not too different to that of radiation.

preprint2025arXiv

Improving CMB constraints on early Universe physics with LSS: A multi-probe forecast including cross-covariance

Extensions to the $Λ\textrm{CDM}$ model prior to recombination can modify the growth of perturbations around radiation-matter equality, leaving a distinct signature in the matter power spectrum. Upcoming large-scale structure surveys will be sensitive to these features, allowing tests of early physics that are complementary to the CMB observations. In this paper, we forecast how well the combination of galaxy clustering, weak lensing and CMB lensing two point statistics, also known as $6\times2$pt analysis, will tighten constraints on extensions to the $Λ\textrm{CDM}$ model in the early Universe. We find significant improvements, in particular in the case of early dark energy, where the uncertainty on its density parameter could be divided by a factor of $3$ to $4$ when combining Euclid observables with Simons Observatory or CMB-S4, compared to using CMB observations alone. Testing for different scale cuts, we find that much of the constraining power comes from the largest scales which are less prone to systematic uncertainties. We take into account the most significant terms in the cross-covariance between large-scale structure tracers and CMB power spectra, which arises from gravitational lensing. Assessing the impact of this additional cross-covariance on the constraints, we find small corrections for most parameters, except for $A_s$ and $τ$ where the lensing induced covariance leads to a more significant degradation of constraints. This forecast analysis highlights the potential of combining CMB and galaxy survey data to test the cosmological model. In particular, early Universe physics, relevant before recombination, stands out as a promising area that benefits substantially from this approach.

preprint2025arXiv

Unified dark fluid with null sound speed as an alternative to phantom dark energy

Recent BAO measurements from DESI, when combined with CMB and supernovae data, suggest evolving dark energy and in particular point to a possible phantom regime, with an equation of state parameter $w<-1$. We explore an alternative phenomenological way to model dark matter and dark energy based on a unified dark fluid (UDF). By construction, our model reproduces the same background expansion history as DESI's best-fit using the CPL parametrization, but assumes a vanishing rest-frame sound speed and no anisotropic stress. This simple prescription ensures a consistent and physical treatment of perturbations and, in our case, the use of a unified dark sector avoids phantom behaviour. We model CMB, LSS, and redshift-space distortion observables, and find mostly small differences with CPL, suggesting that while stage IV CMB and galaxy surveys will be able to test these models, achieving a decisive distinction between them may prove challenging on linear scales. At the non-linear level, we study spherical collapse in the UDF and show that within this framework, structure formation proceeds very similarly to standard scenarios. Using Planck, DESI BAO DR2, and DES Y5 supernovae data, we demonstrate that this simple UDF model fits current observations nearly as well as CPL, while treating perturbations consistently. Because most cosmological observations are not sensitive to how the dark sector is split, the unified framework can also approximate the phenomenology of interacting dark energy-dark matter scenarios or evolving dark matter, making it a general way to model the data, at least as long as the dark components have a vanishing sound speed, which is the most distinctive feature of our analysis. Our results highlight that a unified dark fluid with evolving equation of state and null sound speed is sufficient to pass current constraints without invoking a phantom component.

preprint2022arXiv

Constraints on $τ_\mathrm{NL}$ from Planck temperature and polarization

We update constraints on the amplitude of the primordial trispectrum, using the final Planck mission temperature and polarization data. In the squeezed limit, a cosmological local trispectrum would be observed as a spatial modulation of small-scale power on the CMB sky. We reconstruct this signal as a source of statistical anisotropy via quadratic estimator techniques. We systematically demonstrate how the estimated power spectrum of a reconstructed modulation field can be translated into a constraint on $τ_\mathrm{NL}$ via likelihood methods, demonstrating the procedures effectiveness by inferring known $τ_\mathrm{NL}$ signal(s) from simulations. Our baseline results constrain $τ_\mathrm{NL} < 1700$ at the 95\% confidence level, providing the most stringent constraints to date.

preprint2014arXiv

Accuracy of cosmological parameters using the baryon acoustic scale

Percent-level measurements of the comoving baryon acoustic scale standard ruler can be used to break degeneracies in parameter constraints from the CMB alone. The sound horizon at the epoch of baryon drag is often used as a proxy for the scale of the peak in the matter density correlation function, and can conveniently be calculated quickly for different cosmological models. However, the measurements are not directly constraining this scale, but rather a measurement of the full correlation function, which depends on the detailed evolution through decoupling. We assess the level of reliability of parameter constraints based on a simple approximation of the acoustic scale compared to a more direct determination from the full numerical two-point correlation function. Using a five-parameter fitting technique similar to recent BAO data analyses, we find that for standard ΛCDM models and extensions with massive neutrinos and additional relativistic degrees of freedom, the approximation is at better than 0.15% for most parameter combinations varying over reasonable ranges.

preprint2014arXiv

CMB lensing reconstruction using cut sky polarization maps and pure-$B$ modes

Detailed measurements of the CMB lensing signal are an important scientific goal of ongoing ground-based CMB polarization experiments, which are mapping the CMB at high resolution over small patches of the sky. In this work we simulate CMB polarization lensing reconstruction for the $EE$ and $EB$ quadratic estimators with current-generation noise levels and resolution, and show that without boundary effects the known and expected zeroth and first order $N^{(0)}$ and $N^{(1)}$ biases provide an adequate model for non-signal contributions to the lensing power spectrum estimators. Small sky areas present a number of additional challenges for polarization lensing reconstruction, including leakage of $E$ modes into $B$ modes. We show how simple windowed estimators using filtered pure-$B$ modes can greatly reduce the mask-induced mean-field lensing signal and reduce variance in the estimators. This provides a simple method (used with recent observations) that gives an alternative to more optimal but expensive inverse-variance filtering.

preprint2014arXiv

Impact of polarisation on the intrinsic CMB bispectrum

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.

preprint2014arXiv

Rayleigh scattering: blue sky thinking for future CMB observations

Rayleigh scattering from neutral hydrogen during and shortly after recombination causes the CMB anisotropies to be significantly frequency dependent at high frequencies. This may be detectable with Planck, and would be a strong signal in any future space-based CMB missions. The later peak of the Rayleigh visibility compared to Thomson scattering gives an increased large-scale CMB polarization signal that is a greater than 4% effect for observed frequencies greater than 500GHz. There is a similar magnitude suppression on small scales from additional damping. Due to strong correlation between the Rayleigh and primary signal, measurement of the Rayleigh component is limited by noise and foregrounds, not cosmic variance of the primary CMB, and should observable over a wide range of angular scales at frequencies between roughly 200GHz and 800GHz. I give new numerical calculations of the temperature and polarization power spectra, and show that future CMB missions could measure the temperature Rayleigh cross-spectrum at high precision, detect the polarization from Rayleigh scattering, and also accurately determine the cross-spectra between the Rayleigh temperature signal and primary polarization. The Rayleigh scattering signal may provide a powerful consistency check on recombination physics. In principle it can be used to measure additional horizon-scale primordial perturbation modes at recombination, and distinguish a significant tensor mode B-polarization signal from gravitational lensing at the power spectrum level.

preprint2013arXiv

Efficient sampling of fast and slow cosmological parameters

Physical parameters are often constrained from the data likelihoods using sampling methods. Changing some parameters can be much more computationally expensive (`slow') than changing other parameters (`fast parameters'). I describe a method for decorrelating fast and slow parameters so that parameter sampling in the full space becomes almost as efficient as sampling in the slow subspace when the covariance is well known and the distributions are simple. This gives a large reduction in computational cost when there are many fast parameters. The method can also be combined with a fast 'dragging' method proposed by Neal (2005) that can be more robust and efficient when parameters cannot be fully decorrelated a priori or have more complicated dependencies. I illustrate these methods for the case of cosmological parameter estimation using data likelihoods from the Planck satellite observations with dozens of fast nuisance parameters, and demonstrate a speed up by a factor of five or more. In more complicated cases, especially where the fast subspace is very fast but complex or highly correlated, the fast-slow sampling methods can in principle give arbitrarily large performance gains. The new samplers are implemented in the latest version of the publicly available CosmoMC code.

preprint2013arXiv

On the joint analysis of CMB temperature and lensing-reconstruction power spectra

Gravitational lensing provides a significant source of cosmological information in modern CMB parameter analyses. It is measured in both the power spectrum and trispectrum of the temperature fluctuations. These observables are often treated as independent, although as they are both determined from the same map this is impossible. In this paper, we perform a rigorous analysis of the covariance between lensing power spectrum and trispectrum analyses. We find two dominant contributions coming from: (i) correlations between the disconnected noise bias in the trispectrum measurement and sample variance in the temperature power spectrum; and (ii) sample variance of the lenses themselves. The former is naturally removed when the dominant N0 Gaussian bias in the reconstructed deflection spectrum is dealt with via a partially data-dependent correction, as advocated elsewhere for other reasons. The remaining lens-cosmic-variance contribution is easily modeled but can safely be ignored for a Planck-like experiment, justifying treating the two observable spectra as independent. We also test simple likelihood approximations for the deflection power spectrum, finding that a Gaussian with a parameter-independent covariance performs well.

preprint2012arXiv

CMB lensing and primordial squeezed non-Gaussianity

Squeezed primordial non-Gaussianity can strongly constrain early-universe physics, but it can only be observed on the CMB after it has been gravitationally lensed. We give a new simple non-perturbative prescription for accurately calculating the effect of lensing on any squeezed primordial bispectrum shape, and test it with simulations. We give the generalization to polarization bispectra, and discuss the effect of lensing on the trispectrum. We explain why neglecting the lensing smoothing effect does not significantly bias estimators of local primordial non-Gaussianity, even though the change in shape can be >~10%. We also show how tau_NL trispectrum estimators can be well approximated by much simpler CMB temperature modulation estimators, and hence that there is potentially a ~10-30% bias due to very large-scale lensing modes, depending on the range of modulation scales included. Including dipole sky modulations can halve the tau_NL error bar if kinematic effects can be subtracted using known properties of the CMB temperature dipole. Lensing effects on the g_NL trispectrum are small compared to the error bar. In appendices we give the general result for lensing of any primordial bispectrum, and show how any full-sky squeezed bispectrum can be decomposed into orthogonal modes of distinct angular dependence.

preprint2012arXiv

CMB power spectrum parameter degeneracies in the era of precision cosmology

Cosmological parameter constraints from the CMB power spectra alone suffer several well-known degeneracies. These degeneracies can be broken by numerical artefacts and also a variety of physical effects that become quantitatively important with high-accuracy data e.g. from the Planck satellite. We study degeneracies in models with flat and non-flat spatial sections, non-trivial dark energy and massive neutrinos, and investigate the importance of various physical degeneracy-breaking effects. We test the CAMB power spectrum code for numerical accuracy, and demonstrate that the numerical calculations are accurate enough for degeneracies to be broken mainly by true physical effects (the integrated Sachs-Wolfe effect, CMB lensing and geometrical and other effects through recombination) rather than numerical artefacts. We quantify the impact of CMB lensing on the power spectra, which inevitably provides degeneracy-breaking information even without using information in the non-Gaussianity. Finally we check the numerical accuracy of sample-based parameter constraints using CAMB and CosmoMC. In an appendix we document recent changes to CAMB's numerical treatment of massive neutrino perturbations, which are tested along with other recent improvements by our degeneracy exploration results.

preprint2012arXiv

Constraining Primordial Magnetism

Primordial magnetic fields could provide an explanation for the galactic magnetic fields observed today, in which case they may also leave interesting signals in the CMB and the small-scale matter power spectrum. We discuss how to approximately calculate the important non-linear magnetic effects within the guise of linear perturbation theory, and calculate the matter and CMB power spectra including the SZ contribution. We then use various cosmological datasets to constrain the form of the magnetic field power spectrum. Using solely large-scale CMB data (WMAP7, QUaD and ACBAR) we find a 95% CL on the variance of the magnetic field at 1 Mpc of B_λ< 6.4 nG. When we include SPT data to constrain the SZ effect, we find a revised limit of B_λ< 4.1 nG. The addition of SDSS Lyman-alpha data lowers this limit even further, roughly constraining the magnetic field to B_λ< 1.3 nG.

preprint2012arXiv

Estimators for CMB Statistical Anisotropy

We use quadratic maximum-likelihood (QML) estimators to constrain models with Gaussian but statistically anisotropic Cosmic Microwave Background (CMB) fluctuations, using CMB maps with realistic sky-coverage and instrumental noise. This approach is optimal when the anisotropy is small, or when checking for consistency with isotropy. We demonstrate the power of the QML approach by applying it to the WMAP data to constrain several models which modulate the observed CMB fluctuations to produce a statistically anisotropic sky. We first constrain an empirically motivated spatial modulation of the observed CMB fluctuations, reproducing marginal evidence for a dipolar modulation pattern with amplitude 7% at L < 60, but demonstrate that the effect decreases at higher multipoles and is 1% at L~500. We also look for evidence of a direction-dependent primordial power spectrum, finding a very statistically significant quadrupole signal nearly aligned with the ecliptic plane; however we argue this anisotropy is largely contaminated by observational systematics. Finally, we constrain the anisotropy due to a spatial modulation of adiabatic and isocurvature primordial perturbations, and discuss the close relationship between anisotropy and non-Gaussianity estimators.

preprint2012arXiv

Likelihood Analysis of CMB Temperature and Polarization Power Spectra

Microwave background temperature and polarization observations are a powerful way to constrain cosmological parameters if the likelihood function can be calculated accurately. The temperature and polarization fields are correlated, partial sky coverage correlates power spectrum estimators at different ell, and the likelihood function for a theory spectrum given a set of observed estimators is non-Gaussian. An accurate analysis must model all these properties. Most existing likelihood approximations are good enough for a temperature-only analysis, however they cannot reliably handle a temperature-polarization correlations. We give a new general approximation applicable for correlated Gaussian fields observed on part of the sky. The approximation models the non-Gaussian form exactly in the ideal full-sky limit and is fast to evaluate using a pre-computed covariance matrix and set of power spectrum estimators. We show with simulations that it is good enough to obtain correct results at ell >~ 30 where an exact calculation becomes impossible. We also show that some Gaussian approximations give reliable parameter constraints even though they do not capture the shape of the likelihood function at each ell accurately. Finally we test the approximations on simulations with realistically anisotropic noise and asymmetric foreground mask.

preprint2012arXiv

Recovering 3D clustering information with angular correlations

We study how to recover the full 3D clustering information of P(\vec{k},z), including redshift space distortions (RSD), from 2D tomography using the angular auto and cross spectra of different redshift bins C_\ell(z,z'). We focus on quasilinear scales where the minimum scale λ_{min} or corresponding maximum wavenumber k_{max}= 2π/λ_{min} is targeted to be between k_{max}={0.05-0.2} h/Mpc. For spectroscopic surveys, we find that we can recover the full 3D clustering information when the redshift bin width Δz used in the 2D tomography is similar to the targeted minimum scale, i.e. Δz ~ {0.6-0.8} λ_{min} H(z)/c which corresponds to Δz ~ 0.01-0.05 for z<1. This value of Δz is optimal in the sense that larger values of Δz lose information, while smaller values violate our minimum scale requirement. For a narrow-band photometric survey, with photo-z error σ_z=0.004, we find almost identical results to the spectroscopic survey because the photo-z error is smaller than the optimal bin width σ_z<Δz. For a typical broad-band photometric survey with σ_z=0.1, we have that σ_z>Δz and most radial information is intrinsically lost. The remaining information can be recovered from the 2D tomography if we use Δz ~ 2σ_z. While 3D and 2D analysis are shown here to be equivalent, the advantage of using angular positions and redshifts is that we do not need a fiducial cosmology to convert to 3D coordinates. This avoids assumptions and marginalization over the fiducial model. In addition, it becomes straight forward to combine RSD, clustering and weak lensing in 2D space.

preprint2012arXiv

The full squeezed CMB bispectrum from inflation

The small-scale CMB temperature we observe on the sky is modulated by perturbations that were super-horizon at recombination, giving differential focussing and lensing that generate a non-zero bispectrum even for single-field inflation where local physics is identical. Understanding this signal is important for primordial non-Gaussianity studies and also parameter constraints from the CMB lensing bispectrum signal. Because of cancellations individual effects can appear larger or smaller than they are in total, so a full analysis may be required to avoid biases. I relate angular scales on the sky to physical scales at recombination using the optical equations, and give full-sky results for the large-scale adiabatic temperature bispectrum from Ricci focussing (expansion of the ray bundle), Weyl lensing (convergence and shear), and temperature redshift modulations of small-scale power. The delta N expansion of the beam is described by the constant temperature 3-curvature, and gives a nearly-observable version of the consistency relation prediction from single-field inflation. I give approximate arguments to quantify the likely importance of dynamical effects, and argue that they can be neglected for modulation scales l <~ 100, which is sufficient for lensing studies and also allows robust tests of local primordial non-Gaussianity using only the large-scale modulation modes. For accurate numerical results early and late-time ISW effects must be accounted for, though I confirm that the late-time non-linear Rees-Sciama contribution is negligible compared to other more important complications. The total corresponds to f_NL ~ 7 for Planck-like temperature constraints and f_NL ~ 11 for cosmic-variance limited data to lmax=2000. Temperature lensing bispectrum estimates are affected at the 0.2 sigma level by Ricci focussing, and up to 0.5 sigma with polarization.

preprint2012arXiv

The real shape of non-Gaussianities

I review what bispectra and trispectra look like in real space, in terms of the sign of particular shaped triangles and tetrahedrons. Having an equilateral density bispectrum of positive sign corresponds to having concentrated overdensities surrounded by larger weaker underdensities. In 3D these are concentrated density filaments, as expected in large-scale structure. As the shape changes from equilateral to flattened the concentrated overdensities flatten into lines (3D planes). I then focus on squeezed bispectra, which can be thought of as correlations of changes in small-scale power with large-scale fields, and discuss the general non-perturbative form of the squeezed bispectrum and its angular dependence. A general trispectrum has tetrahedral form and I show examples of what this can look like in real space. Squeezed trispectra are of particular interest and come in two forms, corresponding to large-scale variance of small-scale power, and correlated modulations of an equilateral-form bispectrum. Flattened trispectra can be produced by line-like features in 2D, for example from cosmic strings, and randomly located features also give a non-Gaussian signal. There are relationships between the squeezed types of non-Gaussianity, and also a useful interpretation in terms of statistical anisotropy. I discuss the various possible physical origins of cosmological non-Gaussianities, both in terms of primordial perturbations and late-time dynamical and geometric effects.

preprint2011arXiv

The linear power spectrum of observed source number counts

We relate the observable number of sources per solid angle and redshift to the underlying proper source density and velocity, background evolution and line-of-sight potentials. We give an exact result in the case of linearized perturbations assuming general relativity. This consistently includes contributions of the source density perturbations and redshift distortions, magnification, radial displacement, and various additional linear terms that are small on sub-horizon scales. In addition we calculate the effect on observed luminosities, and hence the result for sources observed as a function of flux, including magnification bias and radial-displacement effects. We give the corresponding linear result for a magnitude-limited survey at low redshift, and discuss the angular power spectrum of the total count distribution. We also calculate the cross-correlation with the CMB polarization and temperature including Doppler source terms, magnification, redshift distortions and other velocity effects for the sources, and discuss why the contribution of redshift distortions is generally small. Finally we relate the result for source number counts to that for the brightness of line radiation, for example 21-cm radiation, from the sources.

preprint2011arXiv

The shape of the CMB lensing bispectrum

Lensing of the CMB generates a significant bispectrum, which should be detected by the Planck satellite at the 5-sigma level and is potentially a non-negligible source of bias for f_NL estimators of local non-Gaussianity. We extend current understanding of the lensing bispectrum in several directions: (1) we perform a non-perturbative calculation of the lensing bispectrum which is ~10% more accurate than previous, first-order calculations; (2) we demonstrate how to incorporate the signal variance of the lensing bispectrum into estimates of its amplitude, providing a good analytical explanation for previous Monte-Carlo results; and (3) we discover the existence of a significant lensing bispectrum in polarization, due to a previously-unnoticed correlation between the lensing potential and E-polarization as large as 30% at low multipoles. We use this improved understanding of the lensing bispectra to re-evaluate Fisher-matrix predictions, both for Planck and cosmic variance limited data. We confirm that the non-negligible lensing-induced bias for estimation of local non-Gaussianity should be robustly treatable, and will only inflate f_NL error bars by a few percent over predictions where lensing effects are completely ignored (but note that lensing must still be accounted for to obtain unbiased constraints). We also show that the detection significance for the lensing bispectrum itself is ultimately limited to 9 sigma by cosmic variance. The tools that we develop for non-perturbative calculation of the lensing bispectrum are directly relevant to other calculations, and we give an explicit construction of a simple non-perturbative quadratic estimator for the lensing potential and relate its cross-correlation power spectrum to the bispectrum. Our numerical codes are publicly available as part of CAMB and LensPix.

preprint2010arXiv

Asymmetric Beams and CMB Statistical Anisotropy

Beam asymmetries result in statistically-anisotropic cosmic microwave background (CMB) maps. Typically, they are studied for their effects on the CMB power spectrum, however they more closely mimic anisotropic effects such as gravitational lensing and primordial power asymmetry. We discuss tools for studying the effects of beam asymmetry on general quadratic estimators of anisotropy, analytically for full-sky observations as well as in the analysis of realistic data. We demonstrate this methodology in application to a recently-detected 9 sigma quadrupolar modulation effect in the WMAP data, showing that beams provide a complete and sufficient explanation for the anomaly.

preprint2010arXiv

Improving CMB non-Gaussianity estimators using tracers of local structure

Local non-Gaussianity causes correlations between large scale perturbation modes and the small scale power. The large-scale CMB signal has contributions from the integrated Sachs Wolfe (ISW) effect, which does not correlate with the small scale power. If this ISW contribution can be removed, the sensitivity to local non-Gaussianity is improved. Gravitational lensing and galaxy counts can be used to trace the ISW contribution; in particular we show that the CMB lensing potential is highly correlated with the ISW signal. We construct a nearly-optimal estimator for the local non-Gaussianity parameter $\fnl$ and investigate to what extent we can use this to decrease the variance on ${\fnl}$. We show that the variance can be decreased by up to $20\%$ at Planck sensitivity using galaxy counts. CMB lensing is a good bias-independent ISW tracer for future more sensitive observations, though the fractional decrease in variance is small if good polarization data is also available.

preprint2010arXiv

Lensed CMB simulation and parameter estimation

Modelling of the weak lensing of the CMB will be crucial to obtain correct cosmological parameter constraints from forthcoming precision CMB anisotropy observations. The lensing affects the power spectrum as well as inducing non-Gaussianities. We discuss the simulation of full sky CMB maps in the weak lensing approximation and describe a fast numerical code. The series expansion in the deflection angle cannot be used to simulate accurate CMB maps, so a pixel remapping must be used. For parameter estimation accounting for the change in the power spectrum but assuming Gaussianity is sufficient to obtain accurate results up to Planck sensitivity using current tools. A fuller analysis may be required to obtain accurate error estimates and for more sensitive observations. We demonstrate a simple full sky simulation and subsequent parameter estimation at Planck-like sensitivity. The lensed CMB simulation and parameter estimation codes are publicly available.

preprint2010arXiv

Massive Neutrinos and Magnetic Fields in the Early Universe

Primordial magnetic fields and massive neutrinos can leave an interesting signal in the CMB temperature and polarization. We perform a systematic analysis of general perturbations in the radiation-dominated universe, accounting for any primordial magnetic field and including leading- order effects of the neutrino mass. We show that massive neutrinos qualitatively change the large- scale perturbations sourced by magnetic fields, but that the effect is much smaller than previously claimed. We calculate the CMB power spectra sourced by inhomogeneous primordial magnetic fields, from before and after neutrino decoupling, including scalar, vector and tensor modes, and consistently modelling the correlation between the density and anisotropic stress sources. In an appendix we present general series solutions for the possible regular primordial perturbations.

preprint2009arXiv

Results of the GREAT08 Challenge: An image analysis competition for cosmological lensing

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

preprint2009arXiv

Weak lensing of the CMB

The cosmic microwave background (CMB) represents a unique source for the study of gravitational lensing. It is extended across the entire sky, partially polarized, located at the extreme distance of z=1100, and is thought to have the simple, underlying statistics of a Gaussian random field. Here we review the weak lensing of the CMB, highlighting the aspects which differentiate it from the weak lensing of other sources, such as galaxies. We discuss the statistics of the lensing deflection field which remaps the CMB, and the corresponding effect on the power spectra. We then focus on methods for reconstructing the lensing deflections, describing efficient quadratic maximum-likelihood estimators and delensing. We end by reviewing recent detections and observational prospects.

preprint1999arXiv

Quadratic Lagrangians and Topology in Gauge Theory Gravity

We consider topological contributions to the action integral in a gauge theory formulation of gravity. Two topological invariants are found and are shown to arise from the scalar and pseudoscalar parts of a single integral. Neither of these action integrals contribute to the classical field equations. An identity is found for the invariants that is valid for non-symmetric Riemann tensors, generalizing the usual GR expression for the topological invariants. The link with Yang-Mills instantons in Euclidean gravity is also explored. Ten independent quadratic terms are constructed from the Riemann tensor, and the topological invariants reduce these to eight possible independent terms for a quadratic Lagrangian. The resulting field equations for the parity non-violating terms are presented. Our derivations of these results are considerably simpler that those found in the literature.