Source author record

Peter Coles

Peter Coles 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

35works
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

35 published item(s)

preprint2022arXiv

A Revolution in Science: the Eclipse Expeditions of 1919

The first direct experimental test of Einstein's theory of general relativity involved a pair of expeditions to measure the bending of light at a total solar eclipse that took place one hundred years ago, on 29 May 1919. So famous is this experiment, and so dramatic was the impact on Einstein himself, that history tends not to recognise the controversy that surrounded the results at the time. In this article, I discuss the experiment in its scientific and historical background context and explain why it was, and is, such an important episode in the development of modern physics.

preprint2022arXiv

Gravitational Stability of Vortices in Bose-Einstein Condensate Dark Matter

We investigate a simple model for a galactic halo under the assumption that it is dominated by a dark matter component in the form of a Bose-Einstein condensate involving an ultra-light scalar particle. In particular we discuss the possibility if the dark matter is in superfluid state then a rotating galactic halo might contain quantised vortices which would be low-energy analogues of cosmic strings. Using known solutions for the density profiles of such vortices we compute the self-gravitational interactions in such halos and place bounds on the parameters describing such models, such as the mass of the particles involved.

preprint2016arXiv

Dynamical Analysis of Scalar Field Cosmologies with Spatial Curvature

We explore the dynamical behaviour of cosmological models involving a scalar field (with an exponential potential and a canonical kinetic term) and a matter fluid with spatial curvature included in the equations of motion. Using appropriately defined parameters to describe the evolution of the scalar field energy in this situation, we find that there are two extra fixed points that are not present in the case without curvature. We also analyse the evolution of the effective equation-of-state parameter for different initial values of the curvature.

preprint2015arXiv

Galaxy Clustering in 3D and Modified Gravity Theories

We study Modified Gravity (MG) theories by modelling the redshifted matter power spectrum in a spherical Fourier-Bessel (sFB) basis. We use a fully non-linear description of the real-space matter power-spectrum and include the lowest-order redshift-space correction (Kaiser effect), taking into account some additional non-linear contributions. Ignoring relativistic corrections, which are not expected to play an important role for a shallow survey, we analyse two different modified gravity scenarios, namely the generalised Dilaton scalar-tensor theories and the $f({R})$ models in the large curvature regime. We compute the 3D power spectrum ${\cal C}^s_{\ell}(k_1,k_2)$ for various such MG theories with and without redshift space distortions, assuming precise knowledge of background cosmological parameters. Using an all-sky spectroscopic survey with Gaussian selection function $φ(r)\propto \exp(-{r^2 / r^2_0})$, $r_0 = 150 \, h^{-1} \, {\textrm{Mpc}}$, and number density of galaxies $\bar {\textrm{N}} =10^{-4}\;{\textrm{Mpc}}^{-3}$, we use a $χ^2$ analysis, and find that the lower-order $(\ell \leq 25)$ multipoles of ${\cal C}^s_\ell(k,k')$ (with radial modes restricted to $k < 0.2 \, h \,{\textrm{Mpc}}^{-1}$) can constraint the parameter $f_{R_0}$ at a level of $2\times 10^{-5} (3\times 10^{-5})$ with $3 σ$ confidence for $n=1(2)$. Combining constraints from higher $\ell > 25$ modes can further reduce the error bars and thus in principle make cosmological gravity constraints competitive with solar system tests. However this will require an accurate modelling of non-linear redshift space distortions. Using a tomographic $β(a)$-$m(a)$ parameterization we also derive constraints on specific parameters describing the Dilaton models of modified gravity.

preprint2015arXiv

Testing Modified Gravity with Cosmic Shear

We use the cosmic shear data from the Canada-France-Hawaii Telescope Lensing Survey to place constraints on $f(R)$ and {\it Generalized Dilaton} models of modified gravity. This is highly complimentary to other probes since the constraints mainly come from the non-linear scales: maximal deviations with respects to the General-Relativity + $Λ$CDM scenario occurs at $k\sim1 h \mbox{Mpc}^{-1}$. At these scales, it becomes necessary to account for known degeneracies with baryon feedback and massive neutrinos, hence we place constraints jointly on these three physical effects. To achieve this, we formulate these modified gravity theories within a common tomographic parameterization, we compute their impact on the clustering properties relative to a GR universe, and propagate the observed modifications into the weak lensing $ξ_{\pm}$ quantity. Confronted against the cosmic shear data, we reject the $f(R)$ $\{ |f_{R_0}|=10^{-4}, n=1\}$ model with more than 99.9% confidence interval (CI) when assuming a $Λ$CDM dark matter only model. In the presence of baryonic feedback processes and massive neutrinos with total mass up to 0.2eV, the model is disfavoured with at least 94% CI in all different combinations studied. Constraints on the $\{ |f_{R_0}|=10^{-4}, n=2\}$ model are weaker, but nevertheless disfavoured with at least 89% CI. We identify several specific combinations of neutrino mass, baryon feedback and $f(R)$ or Dilaton gravity models that are excluded by the current cosmic shear data. Notably, universes with three massless neutrinos and no baryon feedback are strongly disfavoured in all modified gravity scenarios studied. These results indicate that competitive constraints may be achieved with future cosmic shear data.

preprint2014arXiv

Cross-correlating Sunyaev-Zel'dovich and Weak Lensing Maps

We present novel statistical tools to cross-correlate frequency cleaned thermal Sunyaev-Zel'dovich (tSZ) maps and tomographic weak lensing (wl) convergence maps. Moving beyond the lowest order cross-correlation, we introduce a hierarchy of mixed higher-order statistics, the cumulants and cumulant correlators, to analyze non-Gaussianity in real space, as well as corresponding polyspectra in the harmonic domain. Using these moments, we derive analytical expressions for the joint two-point probability distribution function (2PDF) for smoothed tSZ (y_s) and convergence (κ_s) maps. The presence of tomographic information allows us to study the evolution of higher order {\em mixed} tSZ-weak lensing statistics with redshift. We express the joint PDFs p_{κy}(κ_s,y_s) in terms of individual one-point PDFs (p_κ(κ_s), p_y(y_s)) and the relevant bias functions (b_κ(κ_s), b_y(y_s)). Analytical results for two different regimes are presented that correspond to the small and large angular smoothing scales. Results are also obtained for corresponding {\em hot spots} in the tSZ and convergence maps. In addition to results based on hierarchical techniques and perturbative methods, we present results of calculations based on the lognormal approximation. The analytical expressions derived here are generic and applicable to cross-correlation studies of arbitrary tracers of large scale structure including e.g. that of tSZ and soft X-ray background.

preprint2013arXiv

New Approaches to Probing Minkowski Functionals

We generalize the concept of the ordinary skew-spectrum to probe the effect of non-Gaussianity on the morphology of Cosmic Microwave Background (CMB) maps in several domains: in real-space (where they are commonly known as cumulant-correlators), and in harmonic and needlet bases. The essential aim is to retain more information than normally contained in these statistics, in order to assist in determining the source of any measured non-Gaussianity, in the same spirit as Munshi & Heavens (2010) skew-spectra were used to identify foreground contaminants to the CMB bispectrum in Planck data. Using a perturbative series to construct the Minkowski Functionals (MFs), we provide a pseudo-Cl based approach in both harmonic and needlet representations to estimate these spectra in the presence of a mask and inhomogeneous noise. Assuming homogeneous noise we present approx- imate expressions for error covariance for the purpose of joint estimation of these spectra. We present specific results for four different models of primordial non-Gaussianity local, equilateral, orthogonal and enfolded models, as well as non-Gaussianity caused by unsubtracted point sources. Closed form results of next-order corrections to MFs too are obtained in terms of a quadruplet of kurt-spectra. We also use the method of modal decomposition of the bispectrum and trispectrum to reconstruct the MFs as an alternative method of reconstruction of morphological properties of CMB maps. Finally, we introduce the odd-parity skew-spectra to probe the odd-parity bispectrum and its impact on the morphology of the CMB sky. Although developed for the CMB, the generic results obtained here can be useful in other areas of cosmology.

preprint2013arXiv

Testing Cosmology with Extreme Galaxy Clusters

Motivated by recent suggestions that a number of observed galaxy clusters have masses which are too high for their given redshift to occur naturally in a standard model cosmology, we use Extreme Value Statistics to construct confidence regions in the mass-redshift plane for the most extreme objects expected in the universe. We show how such a diagram not only provides a way of potentially ruling out the concordance cosmology, but also allows us to differentiate between alternative models of enhanced structure formation. We compare our theoretical prediction with observations, placing currently observed high and low redshift clusters on a mass-redshift diagram and find -- provided we consider the full sky to avoid a posteriori selection effects -- that none are in significant tension with concordance cosmology.

preprint2012arXiv

Probing CMB Secondary Anisotropies through Minkowski Functionals

Secondary contributions to the anisotropy of the Cosmic Microwave Background (CMB), such as the integrated Sachs-Wolfe (ISW) effect, the thermal Sunyaev-Zel'dovich effect (tSZ), and the effect of gravitational lensing, have distinctive non-Gaussian signatures, and full descriptions therefore require information beyond that contained in their power spectra. In this paper we use the recently introduced skew-spectra associated with the Minkowski Functionals (MF) to probe the topology of CMB maps to probe the secondary non-Gaussianity as a function of beam-smoothing in order to separate various contributions. We devise estimators for these spectra in the presence of a realistic observational masks and present expressions for their covariance as a function of instrumental noise. Specific results are derived for the mixed ISW-lensing and tSZ-lensing bispectra as well as contamination due to point sources for noise levels that correspond to the Planck (143 GHz channel) and EPIC (150 GHz channel) experiments. The cumulative signal to noise ration $S/N$ for one-point generalized skewness-parameters can reach an order of ${\cal O}(10)$ for Planck and two orders of magnitude higher for EPIC, i.e. ${\cal O}(10^3)$. We also find that these three spectra skew-spectra are correlated, having correlation coefficients $r \sim 0.5-1.0$; higher $l$ modes are more strongly correlated. Though the values of $S/N$ increase with decreasing noise, the triplets of skew-spectra that determine the MFs bcome more correlated; the $S/N$ ratios of lensing-induced skew-spectra are smaller compared to that of a frequency-cleaned tSZ map.

preprint2012arXiv

Statistical Properties of Thermal Sunyaev-Zel'dovich Maps

At high angular frequencies, beyond the damping tail of the primary power spectrum, the dominant contribution to the power spectrum of cosmic microwave background (CMB) temperature fluctuations is the thermal Sunyaev-Zel'dovich (tSZ) effect. We investigate various important statistical properties of the Sunyaev-Zel'dovich maps, using well-motivated models for dark matter clustering to construct statistical descriptions of the tSZ effect to all orders enabling us to determine the entire probability distribution function (PDF). Any generic deterministic biasing scheme can be incorporated in our analysis and the effects of projection, biasing and the underlying density distribution can be analysed separately and transparently in this approach. We introduce the cumulant correlators as tools to analyse tSZ catalogs and relate them to corresponding statistical descriptors of the underlying density distribution. The statistics of hot spots in frequency-cleaned tSZ maps are also developed in a self-consistent way to an arbitrary order, to obtain results complementary to those found using the halo model. We also consider different beam sizes, to check the extent to which the PDF can be extracted from various observational configurations. The formalism is presented with two specific models for underlying matter clustering: (1) the hierarchical ansatz; and (2) the lognormal distribution. We find both models to be in very good agreement with the simulation results, though the lognormal model has an edge over the hierarchical model.

preprint2012arXiv

Strong Gravitational Lens Statistics using the Herschel-ATLAS

The identification of strong gravitational lenses in large surveys has historically been a rather time consuming exercise. Early data from the Herschel Astrophysical Terahertz Large Area Survey (Herschel-ATLAS) demonstrate that lenses can be identified efficiently at submillimetre wavelengths using a simple flux criteria. Motivated by that development, this work considers the statistical properties of strong gravitational lens systems which have been, and will be, found by the Herschel-ATLAS. Analytical models of lens statistics are tested with the current best estimates for the various model ingredients. These include the cosmological parameters, the mass function and the lens density profile, for which we consider the singular isothermal sphere (SIS) and the Navarro, Frenk & White (NFW) approximations. The five lenses identified in the Herschel-ATLAS Science Demonstration Phase suggest a SIS density profile is preferred, but cannot yet constrain Ω_Λto an accuracy comparable with other methods. The complete Herschel-ATLAS data set should be sufficient to yield competitive constraints on Ω_Λ. Whilst this huge number of lenses has great potential for constraining cosmological parameters, they will be most powerful in constraining uncertainty in astrophysical processes. Further investigation is needed to fully exploit this unprecedented data set.

preprint2011arXiv

Exact Extreme Value Statistics and the Halo Mass Function

Motivated by observations that suggest the presence of extremely massive clusters at uncomfortably high redshifts for the standard cosmological model to explain, we develop a theoretical framework for the study of the most massive haloes, e.g. the most massive cluster found in a given volume, based on Extreme Value Statistics (EVS). We proceed from the exact distribution of the extreme values drawn from a known underlying distribution, rather than relying on asymptotic theory (which is independent of the underlying form), arguing that the former is much more likely to furnish robust statistical results. We illustrate this argument with a discussion of the use of extreme value statistics as a probe of primordial non-Gaussianity.

preprint2011arXiv

From Weak Lensing to non-Gaussianity via Minkowski Functionals

We present a new harmonic-domain approach for extracting morphological information, in the form of Minkowski Functionals (MFs), from weak lensing (WL) convergence maps. Using a perturbative expansion of the MFs, which is expected to be valid for the range of angular scales probed by most current weak-lensing surveys, we show that the study of three generalized skewness parameters is equivalent to the study of the three MFs defined in two dimensions. We then extend these skewness parameters to three associated skew-spectra which carry more information about the convergence bispectrum than their one-point counterparts. We discuss various issues such as noise and incomplete sky coverage in the context of estimation of these skew-spectra from realistic data. Our technique provides an alternative to the pixel-space approaches typically used in the estimation of MFs, and it can be particularly useful in the presence of masks with non-trivial topology. Analytical modeling of weak lensing statistics relies on an accurate modeling of the statistics of underlying density distribution. We apply three different formalisms to model the underlying dark-matter bispectrum: the hierarchical ansatz, halo model and a fitting function based on numerical simulations; MFs resulting from each of these formalisms are computed and compared. We investigate the extent to witch late-time gravity-induced non-Gaussianity (to which weak lensing is primarily sensitive) can be separated from primordial non-Gaussianity and how this separation depends on source redshift and angular scale.

preprint2011arXiv

Models for the Clustering of Far-Infrared and Sub-millimetre selected Galaxies

We discuss and compare two alternative models for the two-point angular correlation function of galaxies detected through the sub-millimetre emission using the Herschel Space Observatory. The first, now-standard Halo Model, which represents the angular correlations as arising from one-halo and two-halo contributions, is flexible but complex and rather unwieldy. The second model is based on a much simpler approach: we incorporate a fitting function method to estimate the matter correlation function with approximate model of the bias inferred from the estimated redshift distribution to find the galaxy angular correlation function. We find that both models give a good account of the shape of the correlation functions obtained from published preliminary studies of the HerMES and H-ATLAS surveys performed using Herschel, and yield consistent estimates of the minimum halo mass within which the sub-millimetre galaxies must reside. We note also that both models predict an inflection in the correlation function at intermediate angular scales, so the presence of the feature in the measured correlation function does not unambiguously indicate the presence of intra-halo correlations. The primary barrier to more detailed interpretation of these clustering measurements lies in the substantial uncertainty surrounding the redshift distribution of the sources.

preprint2011arXiv

Statistical Characterization of Temperature Patterns in Anisotropic Cosmologies

We consider the issue of characterizing the coherent large-scale patterns from CMB temperature maps in globally anisotropic cosmologies. The methods we investigate are reasonably general; the particular models we test them on are the homogeneous but anisotropic relativistic cosmologies described by the Bianchi classification. Although the temperature variations produced in these models are not stochastic, they give rise to a "non-Gaussian" distribution of temperature fluctuations over the sky that is a partial diagnostic of the model. We explore two methods for quantifying non-Gaussian and/or non-stationary fluctuation fields in order to see how they respond to the Bianchi models.We first investigate the behavior of phase correlations between the spherical harmonic modes of the maps. Then we examine the behavior of the multipole vectors of the temperature distribution which, though defined in harmonic space, can indicate the presence of a preferred direction in real space, i.e. on the 2-sphere. These methods give extremely clear signals of the presence of anisotropy when applied to the models we discuss, suggesting that they have some promise as diagnostics of the presence of global asymmetry in the Universe.

preprint2011arXiv

The Morphology of the Thermal Sunyaev-Zel'dovich Sky

At high angular frequencies the thermal Sunyaev-Zel'dovich (tSZ) effect constitutes the dominant signal in the CMB sky. The tSZ effect is caused by large scale pressure fluctuations in the baryonic distribution in the Universe so its statistical properties provide estimates of corresponding properties of the projected 3D pressure fluctuations. It's power spectrum is a sensitive probe of the density fluctuations, and the bispectrum can be used to separate the bias associated with pressure. The bispectrum is often probed with a one-point real-space analogue, the skewness. In addition to the skewness the morphological properties, as probed by the well known Minkowski Functionals (MFs), also require the generalized one-point statistics, which at the lowest order are identical to the skewness parameters. The concept of generalized skewness parameters can be extended to define a set of three associated generalized skew-spectra. We use these skew-spectra to probe the morphology of the tSZ sky or the y-sky. We show how these power spectra can be recovered from the data in the presence of arbitrary mask and noise templates using the well known Pseudo-Cl (PCL) approach for arbitrary beam shape. We also employ an approach based on the halo model to compute the tSZ bispectrum. The bispectrum from each of these models is then used to construct the generalized skew-spectra. We consider the performance of an all-sky survey with Planck-type noise and compare the results against a noise-free ideal experiment using a range of smoothing angles. We find that the skew-spectra can be estimated with very high signal-to-noise ratio from future frequency cleaned tSZ maps that will be available from experiments such as Planck. This will allow their mode by mode estimation for a wide range of angular frequencies and will help us to differentiate them from various other sources of non-Gaussianity.

preprint2011arXiv

The Statistics of Cosmological Lyman-alpha Absorption

We study the effect of the non-Gaussianity induced by gravitational evolution upon the statistical properties of absorption in quasar (QSO) spectra. Using the generic hierarchical ansatz and the lognormal approximation we derive the analytical expressions for the one-point PDF as well as for the joint two-point probability distribution (2PDF) of transmitted fluxes in two neighbouring QSOs. These flux PDFs are constructed in 3D as well as in projection (i.e. in 2D). The PDFs are constructed by relating the lower-order moments, i.e. cumulants and cumulant correlators, of the fluxes to the 3D neutral hydrogen distribution which is, in turn, expressed as a function of the underlying dark matter distribution. The lower-order moments are next modelled using a generating function formalism in the context of a {\em minimal tree-model} for the higher-order correlation hierarchy. These different approximations give nearly identical results for the range of redshifts probed, and we also find a very good agreement between our predictions and outputs of hydrodynamical simulations. The formalism developed here for the joint statistics of flux-decrements concerning two lines of sight can be extended to multiple lines of sight, which could be particularly important for the 3D reconstruction of the cosmic web from QSO spectra (e.g. in the BOSS survey). These statistics probe the underlying projected neutral hydrogen field and are thus linked to "hot-spots" of absorption. The results for the PDF and the bias presented here use the same functional forms of scaling functions that have previously been employed for the modelling of other cosmological observation such as the Sunyaev-Zel'dovich effect.

preprint2011arXiv

Tomography and Weak lensing Statistics

Extending previous studies, we derive generic predictions for lower order cumulants and their correlators for individual tomographic bins as well as between two different bins. We derive the corresponding one- and two-point joint probability distribution function for the tomographic convergence maps from different bins as a function of angular smoothing scale. The modelling of weak lensing statistics is obtained by adopting a detailed prescription for the underlying density contrast. In this paper we concentrate on the convergence field $κ$ and use top-hat filter; though the techniques presented can readily be extended to model the PDF of shear components or to include other windows such as the compensated filter. The functional form for the underlying PDF and bias is modelled in terms of the non-linear or the quasilinear form depending on the smoothing angular scale. Results from other semi-analytical models e.g. the lognormal distribution are also presented. Introducing a reduced convergence for individual bins, we are able to show that the tomographic PDFs and bias for each bin sample the same functional form of the underlying PDF of density contrast but with varying variance. The joint probability distribution of the convergence maps that correspond to two different tomographic bins can be constructed from individual tomographic PDF and bias. We study their dependence on cosmological parameters for source distributions corresponding to the realistic surveys such as LSST and DES. We briefly outline how photometric redshift information can be incorporated in our computation of cumulants, cumulant correlators and the PDFs. Connection of our results to the full 3D calculations is elucidated. Analytical results for inclusion of realistic noise and finite survey size are presented in detail.

preprint2011arXiv

Understanding the WMAP Cold Spot mystery

The first and third year data releases from the WMAP provide evidence of an anomalous Cold Spot (CS) at galactic latitude b=-57deg and longitude l=209deg. We have examined the properties of the CS in some detail in order to assess its cosmological significance. We have performed a cluster analysis of the local extrema in the CMB signal to show that the CS is actually associated with a large group of extrema rather than just one. In the light of this we have re-examined the properties of the WMAP ILC and co-added "cleaned" WCM maps, which have previously been used for the analysis of the properties of the signal in the vicinity of the CS. These two maps have remarkably similar properties on equal latitude rings for |b|>30deg, as well as in the vicinity of the CS. We have also checked the idea that the CMB signal has a non-Gaussian tail, localized in the low multipole components of the signal. For each ring we apply a linear filter with characteristic scale R, dividing the CMB signal in two parts: the filtered part, with characteristic scale above that of the filter R, and the difference between the initial and filtered signal. Using the filter scale as a variable, we can maximize the skewness and kurtosis of the smoothed signal and minimize these statistics for the difference between initial and filtered signal. We have discovered that the shape of the CS is formed primarily by the components of the CMB signal represented by multipoles between 10<=L<=20, with a corresponding angular scale about 5-10 degs. This signal leads to modulation of the whole CMB sky, clearly seen at |b|>30deg in both the ILC and WCM maps, rather than a single localized feature. After subtraction of this modulation, the remaining part of the CMB signal appears to be consistent with statistical homogeneity and Gaussianity.

preprint2010arXiv

Higher Order Statistics for Three-dimensional Shear and Flexion

We introduce a collection of statistics appropriate for the study of spinorial quantities defined in three dimensions, focussing on applications to cosmological weak gravitational lensing studies in 3D. In particular, we concentrate on power spectra associated with three- and four-point statistics, which have the advantage of compressing a large number of typically very noisy modes into a convenient data set. It has been shown previously by \cite{MuHe09} that, for non--Gaussianity studies in the microwave background, such compression can be lossless for certain purposes, so we expect the statistics we define here to capture the bulk of the cosmological information available in these higher-order statistics. We consider the effects of a sky mask and noise, and use Limber's approximation to show how, for high-frequency angular modes, confrontation of the statistics with theory can be achieved efficiently and accurately. We focus on scalar and spinorial fields including convergence, shear and flexion of 3D weak lensing, but many of the results apply for general spin fields.

preprint2010arXiv

Higher-order Convergence Statistics for Three-dimensional Weak Gravitational Lensing

Weak gravitational lensing on a cosmological scales can provide strong constraints both on the nature of dark matter and the dark energy equation of state. Most current weak lensing studies are restricted to (two-dimensional) projections, but tomographic studies with photometric redshifts have started, and future surveys offer the possibility of probing the evolution of structure with redshift. In future we will be able to probe the growth of structure in 3D and put tighter constraints on cosmological models than can be achieved by the use of galaxy redshift surveys alone. Earlier studies in this direction focused mainly on evolution of the 3D power spectrum, but extension to higher-order statistics can lift degeneracies as well as providing information on primordial non-gaussianity. We present analytical results for specific higher-order descriptors, the bispectrum and trispectrum, as well as collapsed multi-point statistics derived from them, i.e. cumulant correlators. We also compute quantities we call the power spectra associated with the bispectrum and trispectrum, the Fourier transforms of the well-known cumulant correlators. We compute the redshift dependence of these objects and study their performance in the presence of realistic noise and photometric redshift errors.

preprint2010arXiv

Higher-order Statistics of Weak Lensing Shear and Flexion

Owing to their more extensive sky coverage and tighter control on systematic errors, future deep weak lensing surveys should provide a better statistical picture of the dark matter clustering beyond the level of the power spectrum. In this context, the study of non-Gaussianity induced by gravity can help tighten constraints on the background cosmology by breaking parameter degeneracies, as well as throwing light on the nature of dark matter, dark energy or alternative gravity theories. Analysis of the shear or flexion properties of such maps is more complicated than the simpler case of the convergence due to the spinorial nature of the fields involved. Here we develop analytical tools for the study of higher-order statistics such as the bispectrum (or trispectrum) directly using such maps at different source redshift. The statistics we introduce can be constructed from cumulants of the shear or flexions, involving the cross-correlation of squared and cubic maps at different redshifts. Typically, the low signal-to-noise ratio prevents recovery of the bispectrum or trispectrum mode by mode. We define power spectra associated with each multi- spectra which compresses some of the available information of higher order multispectra. We show how these can be recovered from a noisy observational data even in the presence of arbitrary mask, which introduces mixing between Electric (E-type) and Magnetic (B-type) polarization, in an unbiased way. We also introduce higher order cross-correlators which can cross-correlate lensing shear with different tracers of large scale structures.

preprint2010arXiv

Primordial Non-Gaussianity from a Joint Analysis of Cosmic Microwave Background Temperature and Polarization

We explore a systematic approach to the analysis of primordial non-Gaussianity using fluctuations in temperature and polarization of the Cosmic Microwave Background (CMB). Following Munshi & Heavens (2009), we define a set of power-spectra as compressed forms of the bispectrum and trispectrum derived from CMB temperature and polarization maps; these spectra compress the information content of the corresponding full multispectra and can be useful in constraining early Universe theories. We generalize the standard pseudo-C_l estimators in such a way that they apply to these spectra involving both spin-0 and spin-2 fields, developing explicit expressions which can be used in the practical implementation of these estimators. While these estimators are suboptimal, they are nevertheless unbiased and robust hence can provide useful diagnostic tests at a relatively small computational cost. We next consider approximate inverse-covariance weighting of the data and construct a set of near-optimal estimators based on that approach. Instead of combining all available information from the entire set of mixed bi- or trispectra, i.e multispectra describing both temperature and polarization information, we provide analytical constructions for individual estimators, associated with particular multispectra. The bias and scatter of these estimators can be computed using Monte-Carlo techniques. Finally, we provide estimators which are completely optimal for arbitrary scan strategies and involve inverse covariance weighting; we present the results of an error analysis performed using a Fisher-matrix formalism at both the one-point and two-point level.

preprint2010arXiv

Temperature and Polarization Patterns in Anisotropic Cosmologies

We study the coherent temperature and polarization patterns produced in homogeneous but anisotropic cosmological models. We show results for all Bianchi types with a Friedman-Robertson-Walker limit (i.e. Types I, V, VII$_{0}$, VII$_{h}$ and IX) to illustrate the range of possible behaviour. We discuss the role of spatial curvature, shear and rotation in the geodesic equations for each model and establish some basic results concerning the symmetries of the patterns produced. We also give examples of the time-evolution of these patterns in terms of the Stokes parameters $I$, $Q$ and $U$.

preprint2010arXiv

The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

preprint2010arXiv

Zonal Modes of Cosmic Microwave Background Temperature Maps

All-sky maps of the cosmic microwave background temperature fluctuations are usually represented by a spherical harmonic decomposition involving modes labelled by their degree l and order m (where -l < m < +l). The zonal modes (i.e those with m = 0) are of particular interest because they vary only with galactic latitude; any anomalous behaviour in them might therefore be an indication of erroneous foreground substraction. We perform a simple statistical analysis of the modes with low l for sky maps derived via different cleaning procedures from the Wilkinson Microwave Anisotropy Probe (WMAP) and show that the zonal modes provide a useful diagnostic of possible systematics.

preprint2009arXiv

Adiabatic versus Isocurvature Non--Gaussianity

We study the extent to which one can distinguish primordial non--Gaussianity (NG) arising from adiabatic and isocurvature perturbations. We make a joint analysis of different NG models based on various inflationary scenarios: local-type and equilateral-type NG from adiabatic perturbations and local-type and quadratic-type NG from isocurvature perturbations together with a foreground contamination by point sources. We separate the Fisher information of the bispectrum of CMB temperature and polarization maps by l for the skew spectrum estimator introduced by Munshi & Heavens (2009) to study the scale dependence of the signal-to-noise ratio of different NG components and their correlations. We find that the adiabatic and the isocurvature modes are strongly correlated, though the phase difference of acoustic oscillations helps to distinguish them. The correlation between local-and equilateral-type is weak, but the two isocurvature modes are too strongly correlated to be discriminated. Point source contamination, to the extent to which it can be regarded as white noise, can be almost completely separated from the primordial components for l>100. Including correlations among the different components, we find that the errors of the NG parameters increase by 20-30% for the WMAP 5-year observation, but 5% for Planck observations.

preprint2009arXiv

New Optimised Estimators for the Primordial Trispectrum

Cosmic microwave background studies of non-Gaussianity involving higher-order multispectra can distinguish between early universe theories that predict nearly identical power spectra. However, the recovery of higher-order multispectra is difficult from realistic data due to their complex response to inhomogeneous noise and partial sky coverage, which are often difficult to model analytically. A traditional alternative is to use one-point cumulants of various orders, which collapse the information present in a multispectrum to one number. The disadvantage of such a radical compression of the data is a loss of information as to the source of the statistical behaviour. A recent study by Munshi & Heavens (2009) has shown how to define the skew spectrum (the power spectra of a certain cubic field, related to the bispectrum) in an optimal way and how to estimate it from realistic data. The skew spectrum retains some of the information from the full configuration-dependence of the bispectrum, and can contain all the information on non-Gaussianity. In the present study, we extend the results of the skew spectrum to the case of two degenerate power-spectra related to the trispectrum. We also explore the relationship of these power-spectra and cumulant correlators previously used to study non-Gaussianity in projected galaxy surveys or weak lensing surveys. We construct nearly optimal estimators for quick tests and generalise them to estimators which can handle realistic data with all their complexity in a completely optimal manner. We show how these higher-order statistics and the related power spectra are related to the Taylor expansion coefficients of the potential in inflation models, and demonstrate how the trispectrum can constrain both the quadratic and cubic terms.

preprint2008arXiv

The Effect of Primordial Non--Gaussianity on the Topology of Large-Scale Structure

We study the effect of primordial non-Gaussianity on the development of large-scale cosmic structure using high-resolution N-body simulations. In particular, we focus on the topological properties of the "cosmic web", quantitatively characterized by the Minkowski Functionals, for models with quadratic non-linearities with different values of the usual non-Gaussianity parameter fNL. In the weakly non-linear regime, we find that analytic formulae derived from perturbation theory agree with the numerical results within a few percent of the amplitude of each MF when |fNL|<1000. In the non-linear regime, the detailed behavior of the MFs as functions of threshold density deviates more strongly from the analytical curves, while the overall amplitude of the primordial non-Gaussian effect remains comparable to the perturbative prediction. When smaller-scale information is included, the influence of primordial non-Gaussianity becomes increasingly significant statistically due to decreasing sample variance. We find that the effect of the primordial non-Gaussianity with |fNL|=50 is comparable to the sample variance of mass density fields with a volume of 0.125(Gpc/h)^3 when they are smoothed by Gaussian filter at a scale of 5Mpc/h. The detectability of this effect in actual galaxy surveys will strongly depend upon residual uncertainties in cosmological parameters and galaxy biasing.

preprint2007arXiv

The One-dimensional Fourier Representation and Large Angular Scale Foreground Contamination in the 3-year Wilkinson Microwave Anisotropy Probe data

We employ the one-dimensional Fourier representation (1DFR) to analyze the 3-year WMAP de-biased internal linear combination (DILC) map and its possible contamination by galactic foregrounds. The 1DFR is a representation of the spherical harmonic coefficients for each l mode using an inverse Fourier transform into one-dimensional curves. Based on the a priori assumption that the CMB signal should be statistically independent of, and consequently have no significant correlation with, any foregrounds, we cross-correlate the 1DFR curves of 2 <= l <= 10 modes, which are claimed by the WMAP team to be free of contamination and suitable for whole sky analysis. We find that 8 out of the 9 modes are negatively cross correlated with the foreground maps, an event which has a probability of only 9/512 ~ 0.0176 for uncorrelated signals. Furthermore, the local extrema of the 1DFR curves between the DILC and those of the foregrounds for l=2 and 6 are correlated with significance level below 0.04. We also discuss the minimum variance optimization method and use the properties of the measured cross-correlation to estimate the possible level of contamination present in the DILC map.

preprint2006arXiv

Testing Gaussian random hypothesis with the cosmic microwave background temperature anisotropies in the three-year WMAP data

We test the hypothesis that the temperature of the cosmic microwave background is consistent with a Gaussian random field defined on the celestial sphere, using de-biased internal linear combination (DILC) map produced from the 3-year WMAP data. We test the phases for spherical harmonic modes with l <= 10 (which should be the cleanest) for their uniformity, randomness, and correlation with those of the foreground templates. The phases themselves are consistent with a uniform distribution, but not for l <= 5, and the differences between phases are not consistent with uniformity. For l=3 and l=6, the phases of the CMB maps cross-correlate with the foregrounds, suggestion the presence of residual contamination in the DLC map even on these large scales. We also use a one-dimensional Fourier representation to assemble a_lm into the ΔT_l(ϕ) for each l mode, and test the positions of the resulting maxima and minima for consistency with uniformity randomness on the unit circle. The results show significant departures at the 0.5% level, with the one-dimensional peaks being concentrated around ϕ=180 degs. This strongly significant alignment with the Galactic meridian, together with the cross-correlation of DILC phases with the foreground maps, strongly suggests that even the lowest spherical harmonic modes in the map are significantly contaminated with foreground radiation.

preprint1996arXiv

Quantifying the Topology of Large-Scale Structure

We propose and investigate a new algorithm for quantifying the topological properties of cosmological density fluctuations. We first motivate this algorithm by drawing a formal distinction between two definitions of relevant topological characteristics, based on concepts, on the one hand, from differential topology and, on the other, from integral geometry. The former approach leads one to concentrate on properties of the contour surfaces which, in turn, leads to the algorithms CONTOUR2D and CONTOUR3D familiar to cosmologists. The other approach, which we adopt here, actually leads to much simpler algorithms in both two and three dimensions. (The 2D algorithm has already been introduced to the astronomical literature.) We discuss the 3D case in some detail and compare results obtained with it to analogous results using the CONTOUR3D algorithm.

preprint1996arXiv

Redshift evolution of clustering

We discuss how the redshift dependence of the observed two-point correlation function of various classes of objects can be related to theoretical predictions. This relation involves first a calculation of the redshift evolution of the underlying matter correlations. The next step is to relate fluctuations in mass to those of any particular class of cosmic objects; in general terms, this means a model for the bias and how it evolves with cosmic epoch. Only after these two effects have been quantified can one perform an appropriate convolution of the non-linearly evolved two-point correlation function of the objects with their redshift distribution to obtain the `observed' correlation function for a given sample. This convolution in itself tends to mask the effect of evolution by mixing amplitudes at different redshifts. We develop a formalism which incorporates these requirements and, in particular, a set of plausible models for the evolution of the bias factor. We apply this formalism to the spatial, angular and projected correlation functions from different samples of high-redshift objects, assuming a simple phenomenological model for the initial power-spectrum and an Einstein-de Sitter cosmological model. We find that our model is roughly consistent with data on the evolution of QSO and galaxy clustering, but only if the effective degree of biasing is small. We discuss the differences between our analysis and other theoretical studies of clustering evolution and argue that the dominant barrier to making definitive predictions is uncertainty about the appropriate form of the bias and its evolution with cosmic epoch.

preprint1995arXiv

Getting the Measure of the Flatness Problem

The problem of estimating cosmological parameters such as $Ω$ from noisy or incomplete data is an example of an inverse problem and, as such, generally requires a probablistic approach. We adopt the Bayesian interpretation of probability for such problems and stress the connection between probability and information which this approach makes explicit. This connection is important even when information is ``minimal'' or, in other words, when we need to argue from a state of maximum ignorance. We use the transformation group method of Jaynes to assign minimally--informative prior probability measure for cosmological parameters in the simple example of a dust Friedman model, showing that the usual statements of the cosmological flatness problem are based on an inappropriate choice of prior. We further demonstrate that, in the framework of a classical cosmological model, there is no flatness problem.

preprint1995arXiv

The Cluster Distribution as a Test of Dark Matter Models. III. the Cluster Velocity Field

We study the large-scale velocity fields traced by galaxy clusters in numerical simulations of a box of side 960 $h^{-1}$ Mpc, and compare them with available data on real clusters. In order to test the reliability of the simulations, which are based on an optimized version of the Zel'dovich approximation, we compare their cluster velocities with those of "exact" N-body simulations, and find a remarkable agreement between the two according to a variety of statistical tests. We analyse Cold Dark Matter models with density parameter in the range $0.2\le Ω_0\le 1$, both with and without the cosmological constant term to provide a flat geometry. We also simulate a Cold+Hot Dark Matter model, with 30% provided by the hot component. Comparison with real data is performed by applying tests based on the cumulative velocity frequency distribution (CVFD) and bulk flow statistics. For the CVFD, we use observational velocity data from different authors. By merging all available data in a combined sample, the CVFD test is able to exclude only open models with $Ω_0\le 0.4$ and flat models with $Ω_0=0.2$. However, the analysis of individual observational samples gives contradictory results; in particular the recent Giovanelli data, which imply much lower velocities (with significantly smaller uncertainties) than the other data, are actually only consistent with the open $Ω_0=0.4$ and the flat $Ω_0=0.2$ models, i.e. the models which are excluded by the other data. Furthermore, CVFD and bulk flow analyses of the Branchini et al. reconstructed velocity data disfavour precisely those models accepted on the grounds of Giovanelli's sample. Finally, we confirm that the Lauer & Postman bulk flow determination is an extremely rare event in the cosmological models we have analysed.