Showing 1–50 of 118 results for author: Coles, P

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

The MillenniumTNG Project: The hydrodynamical full physics simulation and a first look at its galaxy clusters

Authors: Ruediger Pakmor, Volker Springel, Jonathan P. Coles, Thomas Guillet, Christoph Pfrommer, Sownak Bose, Monica Barrera, Ana Maria Delgado, Fulvio Ferlito, Carlos Frenk, Boryana Hadzhiyska, César Hernández-Aguayo, Lars Hernquist, Rahul Kannan, Simon D. M. White

Abstract: Cosmological simulations are an important theoretical pillar for understanding nonlinear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large volume simulations of cosmic structure formation aiming to better understand physical processes o… ▽ More Cosmological simulations are an important theoretical pillar for understanding nonlinear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large volume simulations of cosmic structure formation aiming to better understand physical processes on large scales and to help interpreting upcoming large-scale galaxy surveys. We here focus on the full physics box MTNG740 that computes a volume of $(740\,\mathrm{Mpc})^3$ with a baryonic mass resolution of $3.1\times~10^7\,\mathrm{M_\odot}$ using \textsc{arepo} with $80.6$~billion cells and the IllustrisTNG galaxy formation model. We verify that the galaxy properties produced by MTNG740 are consistent with the TNG simulations, including more recent observations. We focus on galaxy clusters and analyse cluster scaling relations and radial profiles. We show that both are broadly consistent with various observational constraints. We demonstrate that the SZ-signal on a deep lightcone is consistent with Planck limits. Finally, we compare MTNG740 clusters with galaxy clusters found in Planck and the SDSS-8 RedMaPPer richness catalogue in observational space, finding very good agreement as well. However, {\it simultaneously} matching cluster masses, richness, and Compton-$y$ requires us to assume that the SZ mass estimates for Planck clusters are underestimated by $0.2$~dex on average. Thanks to its unprecedented volume for a high-resolution hydrodynamical calculation, the MTNG740 simulation offers rich possibilities to study baryons in galaxies, galaxy clusters, and in large scale structure, and in particular their impact on upcoming large cosmological surveys. △ Less

Submitted 6 December, 2022; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 18 pages, 14 figures, accepted for publication by MNRAS, comments welcome

The SATIN project I: Turbulent multi-phase ISM in Milky Way simulations with SNe feedback from stellar clusters

Authors: Rebekka Bieri, Thorsten Naab, Sam Geen, Jonathan P. Coles, Rüdiger Pakmor, Stefanie Walch

Abstract: We introduce the star formation and Supernova (SN) feedback model of the SATIN (Simulating AGNs Through ISM with Non-Equilibrium Effects) project to simulate the evolution of the star forming multi-phase interstellar medium (ISM) of entire disk galaxies. This galaxy-wide implementation of a successful ISM feedback model naturally covers an order of magnitude in gas surface density, shear and radia… ▽ More We introduce the star formation and Supernova (SN) feedback model of the SATIN (Simulating AGNs Through ISM with Non-Equilibrium Effects) project to simulate the evolution of the star forming multi-phase interstellar medium (ISM) of entire disk galaxies. This galaxy-wide implementation of a successful ISM feedback model naturally covers an order of magnitude in gas surface density, shear and radial motions. It is implemented in the adaptive mesh refinement code RAMSES at a peak resolution of 9 pc. New stars are represented by star cluster (sink) particles with individual SN delay times for massive stars. With SN feedback, cooling and gravity, the galactic ISM develops a realistic three-phase structure. The star formation rates naturally follow observed scaling relations for the local Milky Way gas surface density. SNe drive additional turbulence in the warm (300 K < $T$ < 10$^4$ K) gas and increase the kinetic energy of the cold gas, cooling out of the warm phase. The majority of the gas leaving the galactic ISM is warm and hot with mass loading factors of $3 \le η\le 10$. While the hot gas is leaving the system, the warm and cold gas falls back onto the disc in a galactic fountain flow. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: Submitted to MNRAS

Evolution of Cosmic Voids in the Schrodinger-Poisson Formalism

Authors: Aoibhinn Gallagher, Peter Coles

Abstract: We investigate the evolution of cosmic voids in the Schrodinger Poisson formalism, finding wave mechanical solutions for the dynamics in a standard cosmological background with appropriate boundary conditions. We compare the results in this model to those obtained using the Zel'dovich approximation. We discuss the advantages of studying voids in general and the advantages of the Schrodinger Poisso… ▽ More We investigate the evolution of cosmic voids in the Schrodinger Poisson formalism, finding wave mechanical solutions for the dynamics in a standard cosmological background with appropriate boundary conditions. We compare the results in this model to those obtained using the Zel'dovich approximation. We discuss the advantages of studying voids in general and the advantages of the Schrodinger Poisson description over other approaches. In particular, emphasizing the utility of the free particle approximation. We also discuss a dimensionless number, similar to the Reynolds number, for this system which allows our void solutions to be scaled to systems of different physical dimensions. △ Less

Submitted 17 November, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: 9 pages, 10 figures

Bridging the Gap Between Simply Parametrized and Free-Form Pixelated Models of Galaxy Lenses: The Case of WFI 2033-4723 Quad

Authors: Bernardo Barrera, Liliya L. R. Williams, Jonathan P. Coles, Philipp Denzel

Abstract: We study the radial and azimuthal mass distribution of the lensing galaxy in WFI2033-4723. Mindful of the fact that modeling results depend on modeling assumptions, we examine two very different recent models: simply parametrized (SP) models from the H0LiCOW collaboration, and pixelated free-form (FF) GLASS models. In addition, we fit our own models which are a compromise between the astrophysical… ▽ More We study the radial and azimuthal mass distribution of the lensing galaxy in WFI2033-4723. Mindful of the fact that modeling results depend on modeling assumptions, we examine two very different recent models: simply parametrized (SP) models from the H0LiCOW collaboration, and pixelated free-form (FF) GLASS models. In addition, we fit our own models which are a compromise between the astrophysical grounding of SP, and the flexibility of FF approaches. Our models consist of two offset parametric mass components, and generate many solutions, all fitting the quasar point image data. Among other results, we show that to reproduce point image properties the lensing mass must be lopsided, but the origin of this asymmetry can reside in the main lens plane or along the line of sight. We also show that there is a degeneracy between the slope of the density profile and the magnitude of external shear, and that the models from various modeling approaches are connected not by the mass sheet degeneracy, but by a more generalized transformation. Finally, we discuss interpretation degeneracy which afflicts all mass modeling: inability to correctly assign mass to the main lensing galaxy vs. nearby galaxies or line of sight structures. While this may not be a problem for the determination of $H_0$, interpretation degeneracy may become a major issue for the detailed study of galaxy structure. △ Less

Submitted 27 September, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 13 pages, 7 figures, submitted to OJA

Journal ref: The Open Journal of Astrophysics, Vol. 4, September 29, 2021

The lens SW05 J143454.4+522850: a fossil group at redshift 0.6?

Authors: Philipp Denzel, Onur Çatmabacak, Jonathan P. Coles, Claude Cornen, Robert Feldmann, Ignacio Ferreras, Xanthe Gwyn Palmer, Rafael Küng, Dominik Leier, Prasenjit Saha, Aprajita Verma

Abstract: Fossil groups are considered the end product of natural galaxy group evolution in which group members sink towards the centre of the gravitational potential due to dynamical friction, merging into a single, massive, and X-ray bright elliptical. Since gravitational lensing depends on the mass of a foreground object, its mass concentration, and distance to the observer, we can expect lensing effects… ▽ More Fossil groups are considered the end product of natural galaxy group evolution in which group members sink towards the centre of the gravitational potential due to dynamical friction, merging into a single, massive, and X-ray bright elliptical. Since gravitational lensing depends on the mass of a foreground object, its mass concentration, and distance to the observer, we can expect lensing effects of such fossil groups to be particularly strong. This paper explores the exceptional system $\mathrm{J}143454.4+522850$. We combine gravitational lensing with stellar population-synthesis to separate the total mass of the lens into stars and dark matter. The enclosed mass profiles are contrasted with state-of-the-art galaxy formation simulations, to conclude that SW05 is likely a fossil group with a high stellar to dark matter mass fraction $0.027\pm0.003$ with respect to expectations from abundance matching $0.012\pm0.004$, indicative of a more efficient conversion of gas into stars in fossil groups. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 8 pages, 9 figures, submitted to MNRAS

The Hubble constant from eight time-delay galaxy lenses

Authors: Philipp Denzel, Jonathan P. Coles, Prasenjit Saha, Liliya L. R. Williams

Abstract: We present a determination of the Hubble constant from the joint, free-form analysis of 8 strongly, quadruply lensing systems. In the concordance cosmology, we find $H_0 = 71.8^{+3.9}_{-3.3}\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ with a precision of $4.97\%$. This is in agreement with the latest measurements from Supernovae Type Ia and Planck observations of the cosmic microwave backgro… ▽ More We present a determination of the Hubble constant from the joint, free-form analysis of 8 strongly, quadruply lensing systems. In the concordance cosmology, we find $H_0 = 71.8^{+3.9}_{-3.3}\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ with a precision of $4.97\%$. This is in agreement with the latest measurements from Supernovae Type Ia and Planck observations of the cosmic microwave background. Our precision is lower compared to these and other recent time-delay cosmography determinations, because our modelling strategies reflect the systematic uncertainties of lensing degeneracies. We furthermore are able to find reasonable lensed image reconstructions by constraining to either value of $H_0$ from local and early Universe measurements. This leads us to conclude that current lensing constraints on $H_0$ are not strong enough to break the "Hubble tension" problem of cosmology. △ Less

Submitted 19 February, 2021; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 19 pages, 11 figures, published in MNRAS

Journal ref: MNRAS, Volume 501, Issue 1, February 2021, Pages 784-801

Molecular Cross Sections for High Resolution Spectroscopy of Super Earths, Warm Neptunes and Hot Jupiters

Authors: Siddharth Gandhi, Matteo Brogi, Sergei N. Yurchenko, Jonathan Tennyson, Phillip A. Coles, Rebecca K. Webb, Jayne L. Birkby, Gloria Guilluy, George A. Hawker, Nikku Madhusudhan, Aldo S. Bonomo, Alessandro Sozzetti

Abstract: High resolution spectroscopy (HRS) has been used to detect a number of species in the atmospheres of hot Jupiters. Key to such detections is accurately and precisely modelled spectra for cross-correlation against the R$\gtrsim$20,000 observations. There is a need for the latest generation of opacities which form the basis for high signal-to-noise detections using such spectra. In this study we pre… ▽ More High resolution spectroscopy (HRS) has been used to detect a number of species in the atmospheres of hot Jupiters. Key to such detections is accurately and precisely modelled spectra for cross-correlation against the R$\gtrsim$20,000 observations. There is a need for the latest generation of opacities which form the basis for high signal-to-noise detections using such spectra. In this study we present and make publicly available cross sections for six molecular species, H$_2$O, CO, HCN, CH$_4$, NH$_3$ and CO$_2$ using the latest line lists most suitable for low- and high-resolution spectroscopy. We focus on the infrared (0.95-5~$μ$m) and between 500-1500~K where these species have strong spectral signatures. We generate these cross sections on a grid of pressures and temperatures typical for the photospheres of super Earth, warm Neptunes and hot Jupiters using the latest H$_2$ and He pressure broadening. We highlight the most prominent infrared spectral features by modelling three representative exoplanets, GJ~1214~b, GJ~3470~b and HD~189733~b, which encompass a wide range in temperature, mass and radii. In addition, we verify the line lists for H$_2$O, CO and HCN with previous high resolution observations of hot Jupiters. However, we are unable to detect CH$_4$ with our new cross sections from HRS observations of HD~102195~b. These high accuracy opacities are critical for atmospheric detections with HRS and will be continually updated as new data becomes available. △ Less

Submitted 8 April, 2020; originally announced April 2020.

Comments: 16 pages, 11 figures, accepted for publication in MNRAS

ExoMol molecular line lists XXXV: a rotation-vibration line list for hot ammonia

Authors: Phillip A. Coles, Sergei N. Yurchenko, Jonathan Tennyson

Abstract: A new hot line list for $^{14}$NH$_3$ is presented. The line list CoYuTe was constructed using an accurate, empirically refined potential energy surface and a CCSD(T)/aug-cc-pVQZ ab initio dipole moment surface of ammonia, previously reported. The line list is an improvement of the ammonia line list BYTe [Yurchenko et al., Mon. Not. R. Astron. Soc., 413, 1828 (2011)]. The CoYuTe line list covers w… ▽ More A new hot line list for $^{14}$NH$_3$ is presented. The line list CoYuTe was constructed using an accurate, empirically refined potential energy surface and a CCSD(T)/aug-cc-pVQZ ab initio dipole moment surface of ammonia, previously reported. The line list is an improvement of the ammonia line list BYTe [Yurchenko et al., Mon. Not. R. Astron. Soc., 413, 1828 (2011)]. The CoYuTe line list covers wavenumbers up to 20000 cm$^{-1}$, i.e. wavelengths beyond 0.5 $μ$m for temperatures up to 1500 K. Comparisons with the high temperature experimental data from the literature show excellent agrement for wavenumbers below 6000 cm$^{-1}$. The CoYuTe line list contains 16.9 billion transitions and is available from the ExoMol website (www.exomol.com) and the CDS database. △ Less

Submitted 23 November, 2019; originally announced November 2019.

Journal ref: Monthly Notices of the Royal Astronomical Society, 490, 4638-4647 (2019)

Lessons from a blind study of simulated lenses: image reconstructions do not always reproduce true convergence

Authors: Philipp Denzel, Sampath Mukherjee, Jonathan P. Coles, Prasenjit Saha

Abstract: In the coming years, strong gravitational lens discoveries are expected to increase in frequency by two orders of magnitude. Lens-modelling techniques are being developed to prepare for the coming massive influx of new lens data, and blind tests of lens reconstruction with simulated data are needed for validation. In this paper we present a systematic blind study of a sample of 15 simulated strong… ▽ More In the coming years, strong gravitational lens discoveries are expected to increase in frequency by two orders of magnitude. Lens-modelling techniques are being developed to prepare for the coming massive influx of new lens data, and blind tests of lens reconstruction with simulated data are needed for validation. In this paper we present a systematic blind study of a sample of 15 simulated strong gravitational lenses from the EAGLE suite of hydrodynamic simulations. We model these lenses with a free-form technique and evaluate reconstructed mass distributions using criteria based on shape, orientation, and lensed image reconstruction. Especially useful is a lensing analogue of the Roche potential in binary star systems, which we call the $\textit{lensing Roche potential}$. This we introduce in order to factor out the well-known problem of steepness or mass-sheet degeneracy. Einstein radii are on average well recovered with a relative error of ${\sim}5\%$ for quads and ${\sim}25\%$ for doubles; the position angle of ellipticity is on average also reproduced well up to $\pm10^{\circ}$, but the reconstructed mass maps tend to be too round and too shallow. It is also easy to reproduce the lensed images, but optimising on this criterion does not guarantee better reconstruction of the mass distribution. △ Less

Submitted 17 June, 2020; v1 submitted 22 October, 2019; originally announced October 2019.

Comments: 20 pages, 12 figures. Published in MNRAS. Agrees with published version

A variationally computed room temperature line list for AsH$_{3}$

Authors: Phillip A. Coles, Sergei N. Yurchenko, Richard P. Kovacich, James Hobby, Jonathan Tennyson

Abstract: Calculations are reported on the rotation-vibration energy levels of the arsine molecule with associated transition intensities. A potential energy surface (PES) obtained from \textit{ab initio} electronic structure calculations is refined to experimental data, and the resulting energy levels display sub-wavenumber accuracy for all reliably known $J=0$ term values under 6500 cm$^{-1}$. After a sma… ▽ More Calculations are reported on the rotation-vibration energy levels of the arsine molecule with associated transition intensities. A potential energy surface (PES) obtained from \textit{ab initio} electronic structure calculations is refined to experimental data, and the resulting energy levels display sub-wavenumber accuracy for all reliably known $J=0$ term values under 6500 cm$^{-1}$. After a small empirical adjustment of the band centres, our calculated ($J=1-6$) rovibrational states reproduce 578 experimentally derived energies with a root-mean-square error of 0.122 cm$^{-1}$. Absolute line intensities are computed using the refined PES and a new dipole moment surface (DMS) for transitions between states with energies up to 10~500 cm$^{-1}$ and rotational quantum number $J=30$. The computed DMS reproduces experimental line intensities to within 10\% uncertainty for the $ν_1$ and $ν_3$ bands. Furthermore, our calculated absorption cross-sections display good agreement with the main absorption features recorded in Pacific Northwest National Laboratory (PNNL) for the complete range of $600-6500$ cm$^{-1}$. △ Less

Submitted 6 February, 2019; originally announced February 2019.

Comments: Eq. (16) was corrected comparing to the published version

Journal ref: Phys. Chem. Chem. Phys., 2019, 21, 3264-3277

Analysis of gaseous ammonia (NH$_3$) absorption in the visible spectrum of Jupiter - Update

Authors: Patrick G. J. Irwin, Neil Bowles, Ashwin S. Braude, Ryan Garland, Simon Calcutt, Phillip A. Coles, Sergey N. Yurchenko, Jonathan Tennyson

Abstract: An analysis of currently available ammonia (NH$_3$) visible-to-near-infrared gas absorption data was recently undertaken by Irwin et al. (Icarus, 302 (2018) 426) to help interpret Very Large Telescope (VLT) MUSE observations of Jupiter from 0.48 - 0.93 $μ$m, made in support of the NASA/Juno mission. Since this analysis a newly revised set of ammonia line data, covering the previously poorly constr… ▽ More An analysis of currently available ammonia (NH$_3$) visible-to-near-infrared gas absorption data was recently undertaken by Irwin et al. (Icarus, 302 (2018) 426) to help interpret Very Large Telescope (VLT) MUSE observations of Jupiter from 0.48 - 0.93 $μ$m, made in support of the NASA/Juno mission. Since this analysis a newly revised set of ammonia line data, covering the previously poorly constrained range 0.5 - 0.833 $μ$m, has been released by the ExoMol project, "C2018" (Coles et al., JQSRT 219, 199 - 122, 2018), which demonstrates significant advantages over previously available data sets, and providing for the first time complete line data for the previously poorly constrained 5520- and 6475-Å bands of NH$_3$. In this paper we compare spectra calculated using the ExoMol-C2018 data set (Coles et al., JQSRT 219, 199 - 122, 2018) with spectra calculated from previous sources to demonstrate its advantages. We conclude that at the present time the ExoMol-C2018 dataset provides the most reliable ammonia absorption source for analysing low- to medium-resolution spectra of Jupiter in the visible/near-IR spectral range, but note that the data are less able to model high-resolution spectra owing to small, but significant inaccuracies in the line wavenumber estimates. This work is of significance not only for solar system planetary physics, but for future proposed observations of Jupiter-like planets orbiting other stars, such as with NASA's planned Wide-Field Infrared Survey Telescope (WFIRST). △ Less

Submitted 13 December, 2018; originally announced December 2018.

Comments: 12 Figures

Journal ref: Icarus 321 (2019) 572 - 582

An Optimizing Symbolic Algebra Approach for Generating Fast Multipole Method Operators

Authors: Jonathan P. Coles, Rebekka Bieri

Abstract: We have developed a symbolic algebra approach to automatically produce, verify, and optimize computer code for the Fast Multipole Method (FMM) operators. This approach allows for flexibility in choosing a basis set and kernel, and can generate computer code for any expansion order in multiple languages. The procedure is implemented in the publicly available Python program Mosaic. Optimizations per… ▽ More We have developed a symbolic algebra approach to automatically produce, verify, and optimize computer code for the Fast Multipole Method (FMM) operators. This approach allows for flexibility in choosing a basis set and kernel, and can generate computer code for any expansion order in multiple languages. The procedure is implemented in the publicly available Python program Mosaic. Optimizations performed at the symbolic level through algebraic manipulations significantly reduce the number of mathematical operations compared with a straightforward implementation of the equations. We find that the optimizer is able to eliminate 20-80% of the floating-point operations and for the expansion orders $p \le 10$ it changes the observed scaling properties. We present our approach using three variants of the operators with the Cartesian basis set for the harmonic potential kernel $1/r$, including the use of totally symmetric and traceless multipole tensors. △ Less

Submitted 28 May, 2020; v1 submitted 15 November, 2018; originally announced November 2018.

Comments: Updated to final version submitted to Computer Physics Communications. Accepted on 20 November 2019

Journal ref: Computer Physics Communications 251 (2020) 107081

Hyperfine-resolved rotation-vibration line list of ammonia (NH$_3$)

Authors: Phillip A. Coles, Alec Owens, Jochen Küpper, Andrey Yachmenev

Abstract: A comprehensive, hyperfine-resolved rotation-vibration line list for the ammonia molecule ($^{14}$NH$_3$) is presented. The line list, which considers hyperfine nuclear quadrupole coupling effects, has been computed using robust, first principles methodologies based on a highly accurate empirically refined potential energy surface. Transitions between levels with energies below $8000$~cm$^{-1}$ an… ▽ More A comprehensive, hyperfine-resolved rotation-vibration line list for the ammonia molecule ($^{14}$NH$_3$) is presented. The line list, which considers hyperfine nuclear quadrupole coupling effects, has been computed using robust, first principles methodologies based on a highly accurate empirically refined potential energy surface. Transitions between levels with energies below $8000$~cm$^{-1}$ and total angular momentum $F\leq14$ are considered. The line list shows excellent agreement with a range of experimental data and will significantly assist future high-resolution measurements of NH$_3$, both astronomically and in the laboratory. △ Less

Submitted 18 September, 2018; originally announced September 2018.

The Integrated Bispectrum and Beyond

Authors: Dipak Munshi, Peter Coles

Abstract: The position-dependent power spectrum has been recently proposed as a descriptor of gravitationally induced non-Gaussianity in galaxy clustering, as it is sensitive to the "soft limit" of the bispectrum (i.e. when one of the wave number tends to zero). We generalise this concept to higher order and clarify their relationship to other known statistics such as the skew-spectrum, the kurt-spectra and… ▽ More The position-dependent power spectrum has been recently proposed as a descriptor of gravitationally induced non-Gaussianity in galaxy clustering, as it is sensitive to the "soft limit" of the bispectrum (i.e. when one of the wave number tends to zero). We generalise this concept to higher order and clarify their relationship to other known statistics such as the skew-spectrum, the kurt-spectra and their real-space counterparts the cumulants correlators. Using the {\em Hierarchical Ansatz} (HA) as a toy model for the higher order correlation hierarchy, we show how in the soft limit, polyspectra at a given order can be identified with lower order polyspectra with the same geometrical dependence but with {\em renormalised} amplitudes expressed in terms of amplitudes of the original polyspectra. We extend the concept of position-dependent bispectrum to bispectrum of the divergence of the velocity field $Θ$ and mixed multispectra involving $δ$ and $Θ$ in the 3D perturbative regime. To quantify the effects of transients in numerical simulations, we also present results for lowest order in Lagrangian perturbation theory (LPT) or the Zel'dovich approximation (ZA). Finally, we discuss how to extend the position-dependent spectrum concept to encompass cross-spectra. And finally study the application of this concept to two dimensions (2D), for projected galaxy maps, convergence $κ$ maps from weak-lensing surveys or maps of CMB secondaries e.g. the frequency cleaned $y$ - parameter maps of thermal Sunyaev-Zel'dovich (tSZ) effect from CMB surveys. △ Less

Submitted 15 August, 2016; originally announced August 2016.

Comments: 30 pages, 4 figures. arXiv admin note: text overlap with arXiv:1407.1836 by other authors

Journal ref: JCAP02(2017)010

The ExoMol database: molecular line lists for exoplanet and other hot atmospheres

Authors: Jonathan Tennyson, Sergei N. Yurchenko, Ahmed F. Al-Refaie, Emma J. Barton, Katy L. Chubb, Phillip A. Coles, S. Diamantopoulou, Maire N. Gorman, Christian Hill, Aden Z. Lam, Lorenzo Lodi, Laura K. McKemmish, Yueqi Na, Alec Owens, Oleg L. Polyansky, Clara Sousa-Silva, Daniel S. Underwood, Andrey Yachmenev, Emil Zak

Abstract: The ExoMol database (www.exomol.com) provides extensive line lists of molecular transitions which are valid over extended temperatures ranges. The status of the current release of the database is reviewed and a new data structure is specified. This structure augments the provision of energy levels (and hence transition frequencies) and Einstein $A$ coefficients with other key properties, including… ▽ More The ExoMol database (www.exomol.com) provides extensive line lists of molecular transitions which are valid over extended temperatures ranges. The status of the current release of the database is reviewed and a new data structure is specified. This structure augments the provision of energy levels (and hence transition frequencies) and Einstein $A$ coefficients with other key properties, including lifetimes of individual states, temperature-dependent cooling functions, Landé $g$-factors, partition functions, cross sections, $k$-coefficients and transition dipoles with phase relations. Particular attention is paid to the treatment of pressure broadening parameters. The new data structure includes a definition file which provides the necessary information for utilities accessing ExoMol through its application programming interface (API). Prospects for the inclusion of new species into the database are discussed. △ Less

Submitted 2 October, 2018; v1 submitted 18 March, 2016; originally announced March 2016.

Comments: 68 pages, 24 Tables, 1 figure. Updated corrects minor errors in format statements

Journal ref: J. Molec. Spectrosc., 327, 73-94 (2016)

Lensing-induced morphology changes in CMB temperature maps in modified gravity theories

Authors: D. Munshi, B. Hu, T. Matsubara, P. Coles, A. Heavens

Abstract: Lensing of the Cosmic Microwave Background (CMB) changes the morphology of pattern of temperature fluctuations, so topological descriptors such as Minkowski Functionals can probe the gravity model responsible for the lensing. We show how the recently introduced two-to-two and three-to-one kurt-spectra (and their associated correlation functions), which depend on the power spectrum of the lensing p… ▽ More Lensing of the Cosmic Microwave Background (CMB) changes the morphology of pattern of temperature fluctuations, so topological descriptors such as Minkowski Functionals can probe the gravity model responsible for the lensing. We show how the recently introduced two-to-two and three-to-one kurt-spectra (and their associated correlation functions), which depend on the power spectrum of the lensing potential, can be used to probe modified gravity theories such as $f({R})$ theories of gravity and quintessence models. We also investigate models based on effective field theory, which include the constant-$Ω$ model, and low-energy Ho\vrava theories. Estimates of the cumulative signal-to-noise for detection of lensing-induced morphology changes, reaches ${\cal O}(10^3)$ for the future planned CMB polarization mission COrE$^{+}$. Assuming foreground removal is possible to $\ell_{max}=3000$, we show that many modified gravity theories can be rejected with a high level of significance, making this technique comparable in power to galaxy weak lensing or redshift surveys. These topological estimators are also useful in distinguishing {\em lensing} from other scattering secondaries at the level of the four-point function or trispectrum. Examples include the kinetic Sunyaev-Zel'dovich (kSZ) effect which shares, with lensing, a lack of spectral distortion. We also discuss the complication of foreground contamination from unsubtracted point sources. △ Less

Submitted 2 February, 2016; originally announced February 2016.

Comments: 32 pages, 7 figures

Journal ref: 2016, JCAP, 04, 056

Extracting the late-time kinetic Sunyaev-Zel'dovich effect

Authors: D. Munshi, I. T. Iliev, K. L. Dixon, P. Coles

Abstract: We propose a novel technique to separate the late-time, post-reionization component of the kinetic Sunyaev-Zeldovich (kSZ) effect from the contribution to it from a (poorly understood and probably patchy) reionization history. The kSZ effect is one of the most promising probe of the {\em missing baryons} in the Universe. We study the possibility of reconstructing it in three dimensions (3D), using… ▽ More We propose a novel technique to separate the late-time, post-reionization component of the kinetic Sunyaev-Zeldovich (kSZ) effect from the contribution to it from a (poorly understood and probably patchy) reionization history. The kSZ effect is one of the most promising probe of the {\em missing baryons} in the Universe. We study the possibility of reconstructing it in three dimensions (3D), using future spectroscopic surveys such as the Euclid survey. By reconstructing a 3D template from galaxy density and peculiar velocity fields from spectroscopic surveys we cross-correlate the estimator against CMB maps. The resulting cross-correlation can help us to map out the kSZ contribution to CMB in 3D as a function of redshift thereby extending previous results which use tomographic reconstruction. This allows the separation of the late time effect from the contribution owing to reionization. By construction, it avoids contamination from foregrounds, primary CMB, tSZ effect as well as from star forming galaxies. Due to a high number density of galaxies the signal-to-noise (S/N) for such cross-correlational studies are higher, compared to the studies involving CMB power spectrum analysis. Using a spherical Bessel-Fourier (sFB) transform we introduce a pair of 3D power-spectra: ${\cal C}^{\parallel}_\ell(k)$ and ${\cal C}^{\perp}_\ell(k)$ that can be used for this purpose. We find that in a future spectroscopic survey with near all-sky coverage and a survey depth of $z\approx 1$, reconstruction of ${\cal C}^{\perp}_\ell(k)$ can be achieved in a few radial wave bands $k\approx(0.01-0.5 h^{-1}\rm Mpc)$ with a S/N of upto ${\cal O}(10)$ for angular harmonics in the range $\ell=(200-2000)$ (abrdiged). △ Less

Submitted 11 November, 2015; originally announced November 2015.

Comments: 19 pages, 6 figures

Journal ref: 2016, MNRAS, 463, 2425

Light versus dark in strong-lens galaxies: Dark matter haloes that are rounder than their stars

Authors: Claudio Bruderer, Justin I. Read, Jonathan P. Coles, Dominik Leier, Emilio E. Falco, Ignacio Ferreras, Prasenjit Saha

Abstract: We measure the projected density profile, shape and alignment of the stellar and dark matter mass distribution in 11 strong-lens galaxies. We find that the projected dark matter density profile - under the assumption of a Chabrier stellar initial mass function - shows significant variation from galaxy to galaxy. Those with an outermost image beyond $\sim 10$ kpc are very well fit by a projected NF… ▽ More We measure the projected density profile, shape and alignment of the stellar and dark matter mass distribution in 11 strong-lens galaxies. We find that the projected dark matter density profile - under the assumption of a Chabrier stellar initial mass function - shows significant variation from galaxy to galaxy. Those with an outermost image beyond $\sim 10$ kpc are very well fit by a projected NFW profile; those with images within 10 kpc appear to be more concentrated than NFW, as expected if their dark haloes contract due to baryonic cooling. We find that over several half-light radii, the dark matter haloes of these lenses are rounder than their stellar mass distributions. While the haloes are never more elliptical than $e_{dm} = 0.2$, their stars can extend to $e_* > 0.2$. Galaxies with high dark matter ellipticity and weak external shear show strong alignment between light and dark; those with strong shear ($γ\gtrsim 0.1$) can be highly misaligned. This is reassuring since isolated misaligned galaxies are expected to be unstable. Our results provide a new constraint on galaxy formation models. For a given cosmology, these must explain the origin of both very round dark matter haloes and misaligned strong-lens systems. △ Less

Submitted 3 November, 2015; originally announced November 2015.

Comments: 16 pages, 7 figures, 4 tables. Accepted for publication by MNRAS

Galaxy Clustering in 3D and Modified Gravity Theories

Authors: Dipak Munshi, Geraint Pratten, Patrick Valageas, Peter Coles, Philippe Brax

Abstract: 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 expect… ▽ More 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 $\varphi(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. △ Less

Submitted 3 August, 2015; originally announced August 2015.

Comments: 21 pages, 15 Figures, 1 Table, Submitted to MNRAS

Testing Modified Gravity with Cosmic Shear

Authors: Joachim Harnois-Déraps, Dipak Munshi, Patrick Valageas, Ludovic van Waerbeke, Philippe Brax, Peter Coles, Luca Rizzo

Abstract: 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 th… ▽ More 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. △ Less

Submitted 20 June, 2015; originally announced June 2015.

Comments: 13 pages, 7 figures, 2 tables

Dynamical Analysis of Scalar Field Cosmologies with Spatial Curvature

Authors: Mateja Gosenca, Peter Coles

Abstract: 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… ▽ More 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. △ Less

Submitted 8 February, 2016; v1 submitted 13 February, 2015; originally announced February 2015.

Comments: 17 pages, 11 figures. Amended in response to peer review in the Open Journal of Astrophysics

Journal ref: The Open Journal of Astrophysics, 2016

Reionization and CMB non-Gaussianity

Authors: D. Munshi, P. S. Corasaniti, P. Coles, A. Heavens, S. Pandolfi

Abstract: We show how cross-correlating a high redshift external tracer field, such as the 21cm neutral hydrogen distribution and product maps involving Cosmic Microwave Background (CMB) temperature and polarisation fields, that probe mixed bispectrum involving these fields, can help to determine the reionization history of the Universe, beyond what can be achieved from cross-spectrum analysis. Taking clues… ▽ More We show how cross-correlating a high redshift external tracer field, such as the 21cm neutral hydrogen distribution and product maps involving Cosmic Microwave Background (CMB) temperature and polarisation fields, that probe mixed bispectrum involving these fields, can help to determine the reionization history of the Universe, beyond what can be achieved from cross-spectrum analysis. Taking clues from recent studies for the detection of primordial non-Gaussianity (Munshi & Heavens 2010), we develop a set of estimators that can study reionization using a power spectrum associated with the bispectrum (or skew-spectrum). We use the matched filtering inherent in this method to investigate different reionization histories. We check to what extent they can be used to rule out various models of reionization and study cross contamination from different sources such as the lensing of the CMB. The estimators can be fine-tuned to optimize study of a specific reionization history. We consider three different types of tracers in our study, namely: proto-galaxies; 21cm maps of neutral hydrogen; quasars. We also consider four alternative models of reionization. We find that the cumulative signal-to-noise (S/N) for detection at $\ell_{max}=2000$ can reach $\mathcal{O}(70)$ for cosmic variance limited all-sky experiments. Combining 100GHz, 143GHz and 217GHz channels of the Planck experiment, we find that the (S/N) lies in the range $\mathcal{O}(5)-\mathcal{O}(35)$. The (S/N) depends on the specific choice of a tracer field, and multiple tracers can be effectively used to map out the entire reionization history with reasonable S/N. Contamination from weak lensing is investigated and found to be negligible, and the effects of Thomson scattering from patchy reionizationare also considered. △ Less

Submitted 4 June, 2014; v1 submitted 6 March, 2014; originally announced March 2014.

Comments: 16 pages, 8 figures; MNRAS in press

Journal ref: Mont.Not.Roy.Astron.Soc.442:3427,2014

Probing Modified Gravity Theories with ISW and CMB Lensing

Authors: D. Munshi, B. Hu, A. Renzi, A. Heavens, P. Coles

Abstract: We use the optimised skew-spectrum as well as the skew-spectra associated with the Minkowski Functionals (MFs) to test the possibility of using the cross-correlation of the Integrated Sachs-Wolfe effect (ISW) and lensing of the cosmic microwave background (CMB) radiation to detect deviations in the theory of gravity away from General Relativity (GR). We find that the although both statistics can p… ▽ More We use the optimised skew-spectrum as well as the skew-spectra associated with the Minkowski Functionals (MFs) to test the possibility of using the cross-correlation of the Integrated Sachs-Wolfe effect (ISW) and lensing of the cosmic microwave background (CMB) radiation to detect deviations in the theory of gravity away from General Relativity (GR). We find that the although both statistics can put constraints on modified gravity, the optimised skew-spectra are especially sensitive to the parameter $\rm B_0$ that denotes the the {\em Compton wavelength} of the scalaron at the present epoch. We investigate three modified gravity theories, namely: the Post-Parametrised Friedmanian (PPF) formalism; the Hu-Sawicki (HS) model; and the Bertschinger-Zukin (BZ) formalism. Employing a likelihood analysis for an experimental setup similar to ESA's Planck mission, we find that, assuming GR to be the correct model, we expect the constraints from the first two skew-spectra, $S_{\ell}^{(0)}$ and $S_{\ell}^{(1)}$, to be the same: $\rm B_0<0.45$ at $95%$ confidence level (CL), and $\rm B_0<0.67$ at $99%$ CL in the BZ model. The third skew-spectrum does not give any meaningful constraint. We find that the optimal skew-spectrum provides much more powerful constraint, giving $\rm B_0<0.071$ at $95%$ CL and $\rm B_0<0.15$ at $99%$ CL, which is essentially identical to what can be achieved using the full bispectrum. △ Less

Submitted 6 May, 2014; v1 submitted 4 March, 2014; originally announced March 2014.

Comments: MNRAS in press

Journal ref: 2014, MNRAS, 442(1), 821

Gravitational Lens Recovery with GLASS: Measuring the mass profile and shape of a lens

Authors: Jonathan P. Coles, Justin I. Read, Prasenjit Saha

Abstract: We use a new non-parametric gravitational modelling tool -- \Glass{} -- to determine what quality of data (strong lensing, stellar kinematics, and/or stellar masses) are required to measure the circularly averaged mass profile of a lens and its shape. \Glass{} uses an under-constrained adaptive grid of mass pixels to model the lens, searching through thousands of models to marginalise over model u… ▽ More We use a new non-parametric gravitational modelling tool -- \Glass{} -- to determine what quality of data (strong lensing, stellar kinematics, and/or stellar masses) are required to measure the circularly averaged mass profile of a lens and its shape. \Glass{} uses an under-constrained adaptive grid of mass pixels to model the lens, searching through thousands of models to marginalise over model uncertainties. Our key findings are as follows: (i) for pure lens data, multiple sources with wide redshift separation give the strongest constraints as this breaks the well-known mass-sheet or steepness degeneracy; (ii) a single quad with time delays also performs well, giving a good recovery of both the mass profile and its shape; (iii) stellar masses -- for lenses where the stars dominate the central potential -- can also break the steepness degeneracy, giving a recovery for doubles almost as good as having a quad with time delay data, or multiple source redshifts; (iv) stellar kinematics provide a robust measure of the mass at the half light radius of the stars $r_{1/2}$ that can also break the steepness degeneracy if the Einstein radius $r_E \neq r_{1/2}$; and (v) if $r_E \sim r_{1/2}$, then stellar kinematic data can be used to probe the stellar velocity anisotropy $β$ -- an interesting quantity in its own right. Where information on the mass distribution from lensing and/or other probes becomes redundant, this opens up the possibility of using strong lensing to constrain cosmological models. △ Less

Submitted 29 August, 2014; v1 submitted 30 January, 2014; originally announced January 2014.

Comments: 18 pages, 11 figures, accepted to MNRAS

Probing CMB Secondary Anisotropies through Minkowski Functionals

Authors: Dipak Munshi, Peter Coles, Alan Heavens

Abstract: 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 introdu… ▽ More 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. △ Less

Submitted 26 July, 2012; originally announced July 2012.

Comments: 17 Pages 11 Figures

Strong Gravitational Lens Statistics using the Herschel-ATLAS

Authors: Jo Short, Elizabeth Pearson, Peter Coles, Steve Eales

Abstract: 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 propertie… ▽ More 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. △ Less

Submitted 21 June, 2012; originally announced June 2012.

Comments: 15 pages, 12 figures, submitted to MNRAS

The Statistics of Cosmological Lyman-alpha Absorption

Authors: Dipak Munshi, Peter Coles, Matteo Viel

Abstract: 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 f… ▽ More 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. △ Less

Submitted 19 December, 2011; originally announced December 2011.

Comments: 17 pages, 8 figures

Tomography and Weak lensing Statistics

Authors: Dipak Munshi, Peter Coles, Martin Kilbinger

Abstract: 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 statist… ▽ More 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. △ Less

Submitted 2 December, 2011; originally announced December 2011.

Comments: 19 pages, 10 Figures

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

Authors: Dipak Munshi, Shahab Joudaki, Peter Coles, Joseph Smidt, Scott T. Kay

Abstract: 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 har… ▽ More 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. △ Less

Submitted 2 May, 2014; v1 submitted 21 November, 2011; originally announced November 2011.

Comments: 25 Pages. 17 Figures. Accepted by MNRAS

Journal ref: MNRAS, 2014, 442, 69

Testing Cosmology with Extreme Galaxy Clusters

Authors: Ian Harrison, Peter Coles

Abstract: 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 o… ▽ More 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. △ Less

Submitted 28 July, 2013; v1 submitted 4 November, 2011; originally announced November 2011.

Comments: 5 pages, 3 figures. Small correction to contours in figure 1. Conclusions unchanged

Journal ref: MNRAS 421, L19-L23 (2012)

Exact Extreme Value Statistics and the Halo Mass Function

Authors: Ian Harrison, Peter Coles

Abstract: 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 draw… ▽ More 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. △ Less

Submitted 25 November, 2011; v1 submitted 5 August, 2011; originally announced August 2011.

Comments: 5 pages, 4 figures, minor changes to match published version

Journal ref: MNRAS 418, L20-L24 (2011)

Statistical Properties of Thermal Sunyaev-Zel'dovich Maps

Authors: Dipak Munshi, Shahab Joudaki, Joseph Smidt, Peter Coles, Scott T. Kay

Abstract: 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 st… ▽ More 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. △ Less

Submitted 21 November, 2012; v1 submitted 3 June, 2011; originally announced June 2011.

Comments: 18 pages, 6 figures, Version accepted by MNRAS

Journal ref: MNRAS, (2013), 429, 1564

The Morphology of the Thermal Sunyaev-Zel'dovich Sky

Authors: Dipak Munshi, Joseph Smidt, Shahab Joudaki, Peter Coles

Abstract: 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 densit… ▽ More 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. △ Less

Submitted 25 May, 2011; originally announced May 2011.

Comments: 18 pages, 10 figures, submitted to MNRAS

Journal ref: MNRAS, (2012), 419, 138

From Weak Lensing to non-Gaussianity via Minkowski Functionals

Authors: Dipak Munshi, Ludovic van Waerbeke, Joseph Smidt, Peter Coles

Abstract: 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… ▽ More 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. △ Less

Submitted 9 March, 2011; originally announced March 2011.

Comments: 22 Pages, 12 Figures. Submitting To MNRAS

Journal ref: MNRAS, (2012), 419, 536

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

Authors: Jo Short, Peter Coles

Abstract: 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… ▽ More 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. △ Less

Submitted 10 January, 2011; originally announced January 2011.

Comments: 5 pages, 6 figures, 1 table, accepted for publication in MNRAS

Higher Order Statistics for Three-dimensional Shear and Flexion

Authors: Dipak Munshi, Thomas Kitching, Alan Heavens, Peter Coles

Abstract: 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 convenie… ▽ More 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. △ Less

Submitted 16 December, 2010; originally announced December 2010.

Journal ref: MNRAS, (2011), 416, 629

New Approaches to Probing Minkowski Functionals

Authors: Dipak Munshi, Joseph Smidt, Asantha Cooray, Alessandro Renzi, Alan Heavens, Peter Coles

Abstract: 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 dete… ▽ More 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. △ Less

Submitted 26 June, 2013; v1 submitted 23 November, 2010; originally announced November 2010.

Comments: Version accepted by MNRAS. 27 Pages, 14 Figures

Journal ref: MNRAS, 434, 2830 (2013)

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

Authors: Asantha Cooray, Steve Eales, Scott Chapman, David L. Clements, Olivier Dore, Duncan Farrah, Matt J. Jarvis, Manoj Kaplinghat, Mattia Negrello, Alessandro Melchiorri, Hiranya Peiris, Alexandra Pope, Mario G. Santos, Stephen Serjeant, Mark Thompson, Glenn White, Alexandre Amblard, Manda Banerji, Pier-Stefano Corasaniti, Sudeep Das, Francesco de_Bernardis, Gianfranco de_Zotti, Tommaso Giannantonio, Joaquin Gonzalez-Nuevo Gonzalez, Ali Ahmad Khostovan , et al. (225 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 22 July, 2010; v1 submitted 20 July, 2010; originally announced July 2010.

Comments: White paper supplement to the proposal submitted by the HSLS science team to ESA for Herschel open-time programs

Statistical Characterization of Temperature Patterns in Anisotropic Cosmologies

Authors: Rockhee Sung, Jo Short, Peter Coles

Abstract: 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 sto… ▽ More 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. △ Less

Submitted 11 January, 2011; v1 submitted 12 April, 2010; originally announced April 2010.

Comments: 14 pages, 10 figures, 4 tables, accepted by MNRAS

Temperature and Polarization Patterns in Anisotropic Cosmologies

Authors: Rockhee Sung, Peter Coles

Abstract: 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 es… ▽ More 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$. △ Less

Submitted 15 April, 2010; v1 submitted 6 April, 2010; originally announced April 2010.

Comments: 24 pages, 7 Figures, submitted to JCAP. Revised version: numerous references added, text rewritten, and errors corrected.

Higher-order Statistics of Weak Lensing Shear and Flexion

Authors: Dipak Munshi, Joseph Smidt, Alan Heavens, Peter Coles, Asantha Cooray

Abstract: 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… ▽ More 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. △ Less

Submitted 25 March, 2010; originally announced March 2010.

Comments: 16 pages, 2 figures

Journal ref: MNRAS, (2011), 411, 2241

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

Authors: Dipak Munshi, Peter Coles, Asantha Cooray, Alan Heavens, Joseph Smidt

Abstract: 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 th… ▽ More 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. △ Less

Submitted 26 February, 2010; originally announced February 2010.

Comments: 25 Pages, 4 Figures

Journal ref: MNRAS, (2011), 410, 1295

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

Authors: Dipak Munshi, Alan Heavens, Peter Coles

Abstract: 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.… ▽ More 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. △ Less

Submitted 10 February, 2010; originally announced February 2010.

Comments: 21 pages

Journal ref: MNRAS, (2011), 411, 2161

The Herschel ATLAS

Authors: S. Eales, L. Dunne, D. Clements, A. Cooray, G. De Zotti, S. Dye, R. Ivison, M. Jarvis, G. Lagache, S. Maddox, M. Negrello, S. Serjeant, M. A. Thompson, E. Van Kampen, A. Amblard, P. Andreani, M. Baes, A. Beelen, G. J. Bendo, D. Benford, F. Bertoldi, J. Bock, D. Bonfield, A. Boselli, C. Bridge , et al. (72 additional authors not shown)

Abstract: The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data… ▽ More The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands. △ Less

Submitted 22 October, 2009; originally announced October 2009.

Comments: submitted to PASP

New Optimised Estimators for the Primordial Trispectrum

Authors: Dipak Munshi, Alan Heavens, Asantha Cooray, Joseph Smidt, Peter Coles, Paolo Serra

Abstract: 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… ▽ More 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. △ Less

Submitted 19 October, 2009; originally announced October 2009.

Comments: 19 pages, 2 figures

Journal ref: Mon.Not.Roy.Astron.Soc.412:1993,2011

Detecting relics of a thermal gravitational wave background in the early Universe

Authors: W. Zhao, D. Baskaran, P. Coles

Abstract: A thermal gravitational wave background can be produced in the early Universe if a radiation dominated epoch precedes the usual inflationary stage. This background provides a unique way to study the initial state of the Universe. We discuss the imprint of this thermal spectra of gravitons on the cosmic microwave background (CMB) power spectra, and its possible detection by CMB observations. Assu… ▽ More A thermal gravitational wave background can be produced in the early Universe if a radiation dominated epoch precedes the usual inflationary stage. This background provides a unique way to study the initial state of the Universe. We discuss the imprint of this thermal spectra of gravitons on the cosmic microwave background (CMB) power spectra, and its possible detection by CMB observations. Assuming the inflationary stage is a pure de Sitter expansion we find that, if the number of e-folds of inflation is smaller than 65, the signal of this thermal spectrum can be detected by the observations of Planck and PolarBear experiments, or the planned EPIC experiments. This bound can be even looser if inflation-like stage is the sub-exponential. △ Less

Submitted 30 September, 2009; v1 submitted 24 July, 2009; originally announced July 2009.

Comments: 9 pages, 4 figures, published version

Journal ref: Phys.Lett.B680:411-416,2009

Zonal Modes of Cosmic Microwave Background Temperature Maps

Authors: Jo Short, Peter Coles

Abstract: 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 erroneou… ▽ More 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. △ Less

Submitted 12 April, 2010; v1 submitted 23 July, 2009; originally announced July 2009.

Comments: 5 pages, 4 tables, 2 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.401:2202-2206,2010

Adiabatic versus Isocurvature Non--Gaussianity

Authors: Chiaki Hikage, Dipak Munshi, Alan Heavens, Peter Coles

Abstract: 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 contaminatio… ▽ More 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. △ Less

Submitted 2 July, 2009; originally announced July 2009.

Comments: 7 pages, 3 figures

Journal ref: Mon. Not. Roy. Astron. Soc.404:1505-1511, 2010

Polarized Spots in Anisotropic Open Universes

Authors: Rockhee Sung, Peter Coles

Abstract: We calculate the temperature and polarization patterns generated in anisotropic cosmological models drawn from the Bianchi classification. We show that localized features in the temperature pattern, perhaps similar to the cold spot observed in the Wilkinson Microwave Anisotropy Probe (WMAP) data, can be generated in models with negative spatial curvature, i.e. Bianchi types V and VII$_{h}$. Both… ▽ More We calculate the temperature and polarization patterns generated in anisotropic cosmological models drawn from the Bianchi classification. We show that localized features in the temperature pattern, perhaps similar to the cold spot observed in the Wilkinson Microwave Anisotropy Probe (WMAP) data, can be generated in models with negative spatial curvature, i.e. Bianchi types V and VII$_{h}$. Both these models also generate coherent polarization patterns. In Bianchi VII$_h$, however, rotation of the polarization angle as light propagates along geodesics can convert E modes into B modes but in Bianchi V this is not necessarily the case. It is in principle possible, therefore, to generate localized temperature features without violating existing observational constraints on the odd-parity component of the polarization of the cosmic microwave background. △ Less

Submitted 17 August, 2009; v1 submitted 14 May, 2009; originally announced May 2009.

Comments: 6 pages, 1 figure, revisions to match version published in CQG

Journal ref: Class.Quant.Grav.26:172001,2009