Source author record

Dipak Munshi

Dipak Munshi 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

30works
4topics
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

30 published item(s)

preprint2016arXiv

3D Weak Lensing: Modified Theories of Gravity

Weak lensing (WL) promises to be a particularly sensitive probe of both the growth of large scale structure (LSS) as well as the fundamental relation between matter density perturbations and metric perturbations, thus providing a powerful tool with which we may constrain modified theories of gravity (MG) on cosmological scales. Future deep, wide-field WL surveys will provide an unprecedented opportunity to constrain deviations from General Relativity (GR). Employing a three-dimensional (3D) analysis based on the spherical Fourier-Bessel (sFB) expansion, we investigate the extent to which MG theories will be constrained by a typical 3D WL survey configuration including noise from the intrinsic ellipticity distribution $σ_ε$ of source galaxies. Here we focus on two classes of screened theories of gravity: i) $f(R)$ chameleon models and ii) environmentally dependent dilaton models. We use one-loop perturbation theory combined with halo models in order to accurately model the evolution of matter power-spectrum with redshift in these theories. Using a Fisher information matrix based approach, we show that for an all-sky spectroscopic survey, the parameter $f_{R_0}$ can be constrained in the range $f_{R_0}< 5\times 10^{-6}(9\times 10^{-6})$ for $n=1(2)$ with a 3$σ$ confidence level. This can be achieved by using relatively low order angular harmonics $\ell<100$. Including higher order harmonics $\ell>100$ can further tighten the constraints, making them comparable to current solar-system constraints. We also employ a Principal Component Analysis (PCA) in order to study the parameter degeneracies in the MG parameters. Our results can trivially be extended to other MG theories, such as the K-mouflage models. The confusion from intrinsic ellipticity correlation and modification of the matter power-spectrum at small scale due to feedback mechanisms is briefly discussed.

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

Principal Components of CMB non-Gaussianity

The skew-spectrum statistic introduced by Munshi & Heavens (2010) has recently been used in studies of non-Gaussianity from diverse cosmological data sets including the detection of primary and secondary non-Gaussianity of Cosmic Microwave Background (CMB) radiation. Extending previous work, focussed on independent estimation, here we deal with the question of joint estimation of multiple skew-spectra from the same or correlated data sets. We consider the optimum skew-spectra for various models of primordial non-Gaussianity as well as secondary bispectra that originate from the cross-correlation of secondaries and lensing of CMB: coupling of lensing with the Integrated Sachs-Wolfe (ISW) effect, coupling of lensing with thermal Sunyaev-Zeldovich (tSZ), as well as from unresolved point-sources (PS). For joint estimation of various types of non-Gaussianity, we use the PCA to construct the linear combinations of amplitudes of various models of non-Gaussianity, e.g. $f^{\rm loc}_{\rm NL},f^{\rm eq}_{\rm NL},f^{\rm ortho}_{\rm NL}$ that can be estimated from CMB maps. Bias induced in the estimation of primordial non-Gaussianity due to secondary non-Gaussianity is evaluated. The PCA approach allows one to infer approximate (but generally accurate) constraints using CMB data sets on any reasonably smooth model by use of a lookup table and performing a simple computation. This principle is validated by computing constraints on the DBI bispectrum using a PCA analysis of the standard templates.

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.

preprint2014arXiv

Reconstructing the Thermal Sunyaev-Zel'dovich Effect in 3D

The thermal Sunyaev-Zel'dovich (tSZ) effect measures the line-of-sight projection of the thermal pressure of free electrons and lacks any redshift information. By cross correlating the tSZ effect with an external cosmological tracer we can recover a good fraction of this lost information. Weak lensing (WL) is thought to provide an unbiased probe of the dark Universe, with many WL surveys having sky coverage that overlaps with tSZ surveys. Generalising the tomographic approach, we advocate the use of the spherical Fourier-Bessel (sFB) expansion to perform an analysis of the cross-correlation between the projected (2D) tSZ Compton $y$-parameter maps and 3D weak lensing convergence maps. We use redshift dependent linear biasing and the halo model as a tool to investigate the tSZ-WL cross-correlations in 3D. We use the Press-Schechter (PS) and the Sheth-Tormen (ST) mass-functions in our calculations, finding that the results are quite sensitive to detailed modelling. We provide detailed analysis of surveys with photometric and spectroscopic redshifts. The signal-to-noise (S/N) of the cross-spectra $\mathcal{C}_{\ell} (k)$ for individual 3D modes, defined by the radial and tangential wave numbers $(k;\ell)$, remains comparable to, but below, unity though optimal binning is expected to improve this. The results presented can be generalised to analyse other CMB secondaries, such as the kinetic Sunyaev-Zel'dovich (kSZ) effect.

preprint2013arXiv

Effects of Linear Redshift Space Distortions and Perturbation Theory on BAOs: A 3D Spherical Analysis

The Baryon Acoustic Oscillations (BAO) are features in the matter power spectrum on scales of order 100-150 h^{-1} Mpc that promise to be a powerful tool to constrain and test cosmological models. The BAO have attracted such attention that future upcoming surveys have been designed with the BAO at the forefront of the primary science goals. Recent studies have advocated the use of a spherical-Fourier Bessel (sFB) expansion for future wide field surveys that cover both wide and deep regions of the sky necessitating the simultaneous treatment of the spherical sky geometry as well as the extended radial coverage. Ignoring the possible effects of growth, which is not expected to be significant at low redshifts, we present an extended analysis of the BAOs using the sFB formalism by taking into account the role of non-linearities and linear redshift distortions in the oscillations observed in the galaxy power spectrum. The sFB power spectrum has both radial and tangential dependence and it has been shown that in the limit that we approach a deep survey the sFB power spectrum is purely radial and collapses to the Cartesian Fourier power spectrum. This radialisation of information is shown to hold even in the presence of redshift space distortions (RSD) and 1-loop corrections to the galaxy power spectrum albeit with modified tangential and radial dependence. As per previous studies we find that the introduction of non-linearities leads to a damping of the oscillations in the matter power spectrum.

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.

preprint2012arXiv

Non-Gaussianity in Large Scale Structure and Minkowski Functionals

Minkowski Functionals (MFs) are topological statistics that have become one of many standard tools used for investigating the statistical properties of cosmological random fields. They have found regular use in studies of departures from Gaussianity in a number of important cosmological scenarios. Important examples include the Cosmic Microwave Background (CMB), weak lensing studies, 21cm surveys and large scale structure (LSS). To lowest order the MFs depend on three generalised skewness parameters that can be shown to probe the bispectrum with differing weights. Recent studies have advocated the use of a power spectrum associated with the bispectrum, called the skew-spectrum, that has more power to distinguish between various contributions to the bispectrum than the conventional formalism adopted when using the Minkowski Functionals. In this article we review the motivations for studying non-Gaussianity and emphasize the importance of the momentum dependence of higher order correlators in investigating both inflationary and early Universe models as well as analytical models for gravitational instability. We then introduce the skew-spectra, applied to galaxy surveys, as a tool for investigating various models for primordial and gravitationally induced non-Gaussianities. We present analytical expressions for the skew-spectra for the density field and divergence of the velocity field in 3D and for projected surveys as a function of redshift and a smoothing angular scale. A Gaussian window function is assumed throughout this paper. Analytical results are derived for the case of gravitationally induced non-Gaussianity. These results can be generalised to incorporate redshift space effects. This will be useful in probing primordial and gravitationally induced non-Gaussianity from ongoing and future galaxy surveys.

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.

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

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.

preprint2010arXiv

A Constraint On the Integrated Mass Power Spectrum out to z = 1100 from Lensing of the Cosmic Microwave Background

The temperature fluctuations and polarization of the Cosmic Microwave Background (CMB) are now a well-known probe of the Universe at an infant age of 400,000 years. During the transit to us from the surface of last scattering, the CMB photons are expected to undergo modifications induced by the intervening large-scale structure. Among the expected secondary effects is the weak gravitational lensing of the CMB by the foreground dark matter distribution. We derive a quadratic estimator that uses the non-Gaussianities generated by the lensing effect at the four-point function level to extract the power spectrum of lensing potential fluctuations integrated out to z ~ 1100 with peak contributions from potential fluctuations at z of 2 to 3. Using WMAP 7-year temperature maps, we report the first direct constraints of this lensing potential power spectrum and find that it has an amplitude of A_L = 0.96 +/- 0.60, 1.06 +/- 0.69 and 0.97 +/- 0.47 using the W, V and W+V bands, respectively.

preprint2010arXiv

A Measurement of Cubic-Order Primordial Non-Gaussianity (g_{NL} and τ_{NL}) With WMAP 5-Year Data

We measure two higher-order power spectra involving weighted cubic and squared temperature anisotropy maps from WMAP 5-year data to study the trispectrum generated by primordial non-Gaussianity. Using these measurements combined with Gaussian and noise simulations, we constrain the cubic order non-Gaussianity parameters τ_{NL}, and g_{NL}. With V+W-band data out to l_{max}=600, we find -7.4 < g_{\rm NL}/10^5 < 8.2 and -0.6 < τ_{\rm NL}/10^4 < 3.3 improving the previous COBE-based limit on τ_{\rm NL} < 10^8 nearly four orders of magnitude with WMAP. We find that the ratio of trispectrum to bispectrum amplitude as captured by the ratio of τ_{\rm Nl}/(6f_{\rm NL}/5)^2 ranges from -3 to 21 at the 95% confidence level.

preprint2010arXiv

CMB Constraints on Primordial non-Gaussianity from the Bispectrum (f_{NL}) and Trispectrum (g_{NL} and τ_{NL}) and a New Consistency Test of Single-Field Inflation

We outline the expected constraints on non-Gaussianity from the cosmic microwave background (CMB) with current and future experiments, focusing on both the third (f_{NL}) and fourth-order (g_{NL} and τ_{NL}) amplitudes of the local configuration or non-Gaussianity. The experimental focus is the skewness (two-to-one) and kurtosis (two-to-two and three-to-one) power spectra from weighted maps. In adition to a measurement of τ_{NL} and g_{NL} with WMAP 5-year data, our study provides the first forecasts for future constraints on g_{NL}. We describe how these statistics can be corrected for the mask and cut-sky through a window function, bypassing the need to compute linear terms that were introduced for the previous-generation non-Gaussianity statistics, such as the skewness estimator. We discus the ratio A_{NL} = τ_{NL}/(6f_{NL}/5)^2 as an additional test of single-field inflationary models and discuss the physical significance of each statistic. Using these estimators with WMAP 5-Year V+W-band data out to l_{max}=600 we constrain the cubic order non-Gaussianity parameters τ_{NL}, and g_{NL} and find -7.4 < g_{NL}/10^5 < 8.2 and -0.6 < τ_{NL}/10^4 < 3.3 improving the previous COBE-based limit on τ_{NL} < 10^8 nearly four orders of magnitude with WMAP.

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

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.

preprint2009arXiv

A New Approach to Probing Primordial Non-Gaussianity

We address the dual challenge of estimating deviations from Gaussianity arising in models of the Early Universe, whilst retaining information necessary to assess whether a detection of non-Gaussianity is primordial. We do this by constructing a new statistic, the bispectrum-related power spectrum, which is constructed from a map of the Cosmic Microwave Background. The estimator is optimised for primordial non-Gaussianity detection, but can also be useful in distinguishing primordial non-Gaussianity from secondary non-Gaussianity, such as may arise from unsubtracted point sources, or residuals from component separation. Extending earlier studies we present unbiased non-Gaussianity estimators optimised for partial sky coverage and inhomogeneous noise associated with realistic scan strategies, but which retain the ability to assess foreground contamination.

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.

preprint2009arXiv

Non-Gaussianity in WMAP Data Due to the Correlation of CMB Lensing Potential with Secondary Anisotropies

We measure the skewness power spectrum of the Cosmic Microwave Background (CMB) anisotropies optimized for a detection of the secondary bispectrum generated by the correlation of the CMB lensing potential with integrated Sachs-Wolfe effect and the Sunyaev-Zel'dovich effect. The covariance of our measurements is generated by Monte-Carlo simulations of Gaussian CMB fields with noise properties consistent with Wilkinson Microwave Anisotropy Probe (WMAP) 5-year data. When interpreting multi-frequency measurements we also take into account the confusion resulting from unresolved radio point sources. We analyze Q, V and W-band WMAP 5-year raw and foreground-cleaned maps using the KQ75 mask out to l_max=600. We find no significant evidence for a non-zero non-Gaussian signal from the lensing-secondary correlation in all three bands and we constrain the overall amplitude of the cross power spectrum between CMB lensing potential and the sum of SZ and ISW fluctuations to be 0.42 \pm 0.86 and 1.19 \pm 0.86 in combined V and W-band raw and foreground-cleaned maps provided by the WMAP team, respectively. The point source amplitude at the bispectrum level measured with this skewness power spectrum is higher than previous measurements of point source non-Gaussianity. We also consider an analysis where we also account for the primordial non-Gaussianity in addition to lensing-secondary bispectrum and point sources. The focus of this paper is on secondary anisotropies. Consequently the estimator is not optimised for primordial non-Gaussianity and the limit we find on local non-Gaussianity from the foreground-cleaned V+W maps is f_NL = -13 \pm 62, when marginalized over point sources and lensing-ISW/SZ constributions to the total bispectrum.

preprint2009arXiv

Secondary non-Gaussianity and Cross-Correlation Analysis

We develop optimised estimators of two sorts of power spectra for fields defined on the sky, in the presence of partial sky coverage. The first is the cross-power spectrum of two fields on the sky; the second is the skew spectrum of three fields. These can probe the Integrated Sachs Wolfe Effect (ISW) at large angular scales and the Sunyaev Zéldovich (SZ) effect from hot gas in clusters at small angular scales. The skew spectrum, recently introduced by Munshi & Heavens (2009), is an optimised statistic which can be tuned to study a particular form of non-Gaussianity, such as may arise in the early Universe, but which retains information on the nature of non-Gaussianity. In this paper we develop the mathematical formalism for the skew spectrum of 3 different fields. When applied to the CMB, this allows us to explore the contamination of the skew spectrum by secondary sources of CMB fluctuations. Considering the three-point function, the study of the bispectrum provides valuable information regarding cross-correlation of secondaries with lensing of CMB with much higher significance compared to just the study involving CMB sky alone. We use our analytical models to study specific cases of cosmological interest which include cross-correlating CMB with various large scale tracers to probe ISW and SZ effects for cross spectral analysis and use the formalism to study the signal-to-noise ratio for detection of the weak lensing of the CMB by cross-correlating it with different tracers as well as point sources for CMB experiments such as Planck (abridged).

preprint2004arXiv

Galaxy Clustering and Dark Energy

We study the evolution of galaxy clustering in various cosmological models with quintessence. We investigate how the analytical predictions vary with change of dark energy equation of state $w_X$. Comparing these predictions against available data we discuss to what extent the problems of galaxy biasing can be modelled. This will be key in constraining the dark energy equation of state with future galaxy surveys. We use a compilation of various surveys to study the number density and amplitude of galaxy clustering from observations of the local universe at $z \sim 0$ to that of the Lyman break galaxies and Ly-$α$ emitters at $z \sim 4.9$. We find that there is a degeneracy between the dark energy equation of state and the way galaxies populate dark matter haloes; objects are more biased in models with more negative values of dark energy equation of state $w_X$. We conclude that, while future all sky CMB observations will determine cosmological parameters with unprecedented precision, and cross correlation of weak lensing experiments and galaxy surveys will provide a cleaner and accurate picture of bias associated with collapsed objects, the rate of growth of large scale structure in such surveys can potentially constrain the equation of state of dark energy and the potential of the scalar field associated with quintessence. In particular, we show that the abundance and spatial distribution of galaxy clusters at intermediate redshifts strongly depend on the dark energy equation of state. When accurate measurement of galaxy clustering at high-redshit becomes possible, it will provide constraints on dark energy that are independent and complementary to type Ia supernova studies.

preprint1997arXiv

Modelling the evolution of correlation functions in gravitational clustering

Padmanabhan (1996) has suggested a model to relate the nonlinear two - point correlation function to the linear two - point correlation function. In this paper, we extend this model in two directions: (1) By averaging over the initial Gaussian distribution of density contrasts, we estimate the spectral dependence of the scaling between nonlinear and linear correlation functions. (2) By using a physically motivated ansatz, we generalise the model to N-point correlation functions and relate the nonlinear, volume averaged, N-point correlation function $\barξ_N(x,a)$ with linearly extrapolated volume averaged 2-point correlation function $\bar ξ_2(l,a)$ evaluated at a different scale. We compare the point of transition between different regimes obtained from our model with numerical simulations and show that the spectral dependence of the scaling relations seen in the simulations can be easily understood. Comparison of the calculated form of $\barξ_N$ with the simulations show reasonable agreement. We discuss several implications of the results.