Integrated Cosmological Probe Combination
DISS. ETH NO. 25195
Integrated cosmological probe combination
A thesis submitted to attain the degree of
DOCTOR OF SCIENCES of ETH ZURICH
(Dr. sc. ETH Zurich)
presented by
Andrina Nicola
Master of Science in Physics
ETH Zürich
born on February 5th, 1988
citizen of Roveredo (GR), Switzerland
accepted on the recommendation of
Prof. Dr. Alexandre Réfrégier
Dr. Adam Amara
Prof. Dr. Eiichiro Komatsu
2018
De sui ipsius et multorum ignorantia. — Francesco Petrarca
Abstract
Our Universe is an exciting laboratory, which we can study by observing many of its properties, called cosmological probes. Recent observational progress has led to the establishment of the ΛCDM cosmological model, which relies on the general theory of relativity, dark matter, dark energy and the existence of primordial perturbations. In spite of this progress, several questions remain unanswered, such as the nature of dark matter and dark energy and the details of the initial conditions. Constraining the key components of the standard model of cosmology is thus one of the main goals of observational cosmology today. In this thesis, we aim to improve our understanding of ΛCDM by combining the information from different cosmological probes, taking into account both the auto- and the cross-correlations. This allows us to go beyond the notion of isolated probes, towards an integrated view of the Universe and cosmology. In a first step, we compare the constraining power of a future spectroscopic galaxy redshift survey when analyzed using the spherical harmonic tomography power spectrum and the spherical Fourier-Bessel power spectrum in a Fisher matrix analysis. We then develop a framework for integrated analysis of cosmological probes in which the probes are combined starting at the map level. In a first implementation, we apply this framework to a combination of CMB temperature anisotropy measurements from the Planck satellite as well as galaxy clustering and weak lensing measured by the Sloan Digital Sky Survey. We then extend the analysis to also include CMB lensing from Planck, weak lensing from Dark Energy Survey Science Verification data as well as background probes. These analyses yield self-consistent and competitive constraints on cosmological parameters and provide a confirmation of ΛCDM through the consistency of different probes. As a complement, we quantify possible tensions between these constraints and other existing results through the relative entropy and find all data sets considered to be consistent in the framework of ΛCDM. Finally, we revisit the relative entropy and propose a novel model selection method, which combines relative entropy and posterior predictive distributions. In a series of toy models and applications to cosmological data, we show that this algorithm gives results consistent with expectations and thus appears promising for model selection in cosmology.
i
Introduzione - Zusammenfassung - Zämefassig
L’universo é un appassionante laboratorio, che si lascia decifrare studiando le cosiddette osservabili cosmologiche. La cosmologia ha fatto dei forti progressi negli ultimi anni e questi esperimenti hanno stabilito l’attuale modello standard cosmologico, tecnicamente chiamato ΛCDM. Questo modello si basa sulla teoria della relatività generale, la materia oscura, l’energia oscura e l’esistenza di fluttuazioni primordiali. Malgrado questo progresso, molte questioni rimangono aperte: per esempio la natura della materia e dell’energia oscure, e il meccanismo delle fluttuazioni primordiali. Quindi, uno degli scopi principali della cosmologia sperimen- tale é di studiare queste componenti. In questa tesi usiamo l’informazione contenuta in diverse osservabili per migliorare la nostra conoscenza del modello standard della cosmologia. Allontanandoci quindi dallo studio di singole osservabili, arriviamo ad una visione globale dell’universo e della cosmologia. In einem ersten Projekt vergleichen wir Analysemethoden für eine zukünftige spektroskopi- sche Galaxienkartierung. Wir vergleichen deren Informationsgehalt für eine Analyse mittels Kugelfunktionen und sphärischen Fourier-Bessel Funktionen. Wir entwickeln dann eine Methode um verschiedene kosmologische Observablen gemeinsam zu analysieren. In einer ersten Arbeit verwenden wir diese Methode um Daten der kosmischen Hintergrundstrahlung (aufgenommen mit dem Planck Satelliten) mit Karten von Galaxie- positionen und -scherungen der Sloan Digital Sky Survey zu kombinieren. In einer zweiten Arbeit erweitern wir diese Analyse: (i) auf die Scherung der kosmischen Hintergrundstrahlung, gemessen mit Planck, (ii) auf Karten von Galaxiescherungen der Dark Energy Survey und (iii) auf Messungen des homogenen Universums. Diese Studien erlauben uns die Werte kosmolo- gischer Parameter einzuschränken. Zusätzlich liefert die Übereinstimmung der verschiedenen Datensätze eine Bestätigung des kosmologischen Standardmodells. Wir ergänzen diese Arbeiten indem wir unsere Resultate mit Hilfe der Kullback-Leibler (KL) Divergenz statistisch mit anderen Resultaten vergleichen. Diese Studie zeigt, dass alle be- trachteten Datensätze sowohl miteinander als auch mit dem kosmologischen Standardmodell übereinstimmen. In ere letschte Arbet chömmer nomol uf d KL Divergänz zrugg und schlönd en neui Methode vor, wo chan zum Modäll teschte brucht wärde. Die Methode verbindet d KL Divergänz mit de sogenannte posterior predictive distribution. Mir teschte d Methode i mehrere vereifachte Awändige und zeiged, dass sie s Potenzial hett, um ide Zuekunft i de Kosmologie zum Modäll teschte agwändet z wärde.
iii
Zusammenfassung
Das Universum ist ein aufregendes Labor, das wir durch kosmologische Observablen studie- ren können. In den letzten Jahren gab es grosse Fortschritte in der Kosmologie, und diese Beobachtungen haben zur Aufstellung des kosmologischen Standardmodells, ΛCDM, geführt. Dieses basiert auf der allgemeinen Relativitätstheorie, der Existenz der dunklen Materie, der dunklen Energie und anfänglicher Dichteschwankungen. Trotz dieses Fortschritts bleiben viele Fragen bis heute ungeklärt, wie beispielsweise die physikalische Natur der dunklen Mate- rie und Energie. Eine der Hauptmotivationen der heutigen Kosmologie ist es deshalb, diese Hauptbestandteile des kosmologischen Standardmodells besser zu verstehen. In dieser Arbeit kombinieren wir mehrere kosmologische Observablen, was uns eine umfassende Sicht auf das Universums und die Kosmologie erlaubt. In einem ersten Projekt vergleichen wir Analysemethoden für eine zukünftige spektroskopi- sche Galaxienkartierung. Wir vergleichen deren Informationsgehalt für eine Analyse mittels Kugelfunktionen und sphärischen Fourier-Bessel Funktionen. Wir entwickeln dann eine Methode um verschiedene kosmologische Observablen gemeinsam zu analysieren. In einer ersten Arbeit verwenden wir diese Methode um Daten der kosmischen Hintergrundstrahlung (aufgenommen mit dem Planck Satelliten) mit Karten von Galaxie- positionen und -scherungen der Sloan Digital Sky Survey zu kombinieren. In einer zweiten Arbeit erweitern wir diese Analyse: (i) auf die Scherung der kosmischen Hintergrundstrahlung, gemessen mit Planck, (ii) auf Karten von Galaxiescherungen der Dark Energy Survey und (iii) auf Messungen des homogenen Universums. Diese Studien erlauben uns die Werte kosmolo- gischer Parameter einzuschränken. Zusätzlich liefert die Übereinstimmung der verschiedenen Datensätze eine Bestätigung des kosmologischen Standardmodells. Wir ergänzen diese Arbeiten indem wir unsere Resultate mit Hilfe der Kullback-Leibler (KL) Divergenz statistisch mit anderen Resultaten vergleichen. Diese Studie zeigt, dass alle be- trachteten Datensätze sowohl miteinander als auch mit dem kosmologischen Standardmodell übereinstimmen. In einer letzten Arbeit kehren wir zur KL Divergenz zurück und stellen eine neue Methode zur Modellüberprüfung vor. Diese basiert auf der KL Divergenz und der sogenannten posterior predictive distribution. Wir zeigen anhand verschiedener, vereinfachter Modelle, dass diese Methode das Potential hat, um in Zukunft in der Kosmologie angewandt zu werden.
v
Contents
Abstract i
Introduzione - Zusammenfassung - Zämefassig iii
Zusammenfassungv
1 Introduction 1 1.1 The homogeneous Universe...... 3 1.2 The inhomogeneous Universe...... 7 1.2.1 Linear perturbation theory...... 7 1.2.2 Cosmological structure formation...... 9 1.3 Observational probes...... 12 1.3.1 Statistics of random fields...... 12 1.3.2 Cosmic Microwave Background...... 14 1.3.3 Galaxy clustering...... 15 1.3.4 Cosmic shear...... 18 1.3.5 Background probes...... 21
2 Three-dimensional spherical analyses of cosmological spectroscopic surveys 23 2.1 Introduction...... 23 2.2 Comparison baseline...... 24 2.3 Three-dimensional spherical power spectra...... 26 2.3.1 The Cartesian power spectrum...... 26 2.3.2 The spherical harmonic tomography power spectrum...... 29 2.3.3 The spherical Fourier Bessel power spectrum...... 31 2.4 Fisher matrices for 3D spherical power spectra...... 34 2.4.1 The Fisher matrix for the spherical harmonic tomography power spectrum 34 2.4.2 The Fisher matrix for the SFB power spectrum...... 35 2.5 Results...... 36 2.5.1 Spherical harmonic tomography power spectrum...... 36 2.5.2 Spherical Fourier Bessel power spectrum...... 40 2.5.3 Comparison between spherical harmonic tomography and SFB power spectrum...... 43
vii Contents
2.5.4 Implementation effects on estimated survey constraining power.... 44 2.6 Conclusions...... 45 2.A The radialization of the SFB power spectrum in the presence of RSDs...... 46 2.B The SFB power spectrum for a generic distance-redshift relation...... 46 2.C Derivation of the Fisher matrix for the SFB power spectrum...... 48 2.C.1 Full Fisher matrix...... 48 2.C.2 Diagonal Fisher matrix...... 50 2.D Comparison of FM computation techniques for the SFB power spectrum: Results 50
3 Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing 53 3.1 Introduction...... 53 3.2 Framework...... 55 3.3 Theoretical predictions...... 56 3.4 Maps...... 57 3.4.1 Cosmic Microwave Background...... 57 3.4.2 Galaxy overdensity...... 60 3.4.3 Weak lensing...... 61 3.5 Spherical harmonic power spectra...... 63 3.5.1 CMB...... 64 3.5.2 Galaxy clustering...... 66 3.5.3 Cosmic shear...... 66 3.5.4 CMB and galaxy overdensity cross-correlation...... 68 3.5.5 CMB and weak lensing shear cross-correlation...... 68 3.5.6 Galaxy overdensity and weak lensing shear cross-correlation...... 69 3.6 Covariance matrix...... 69 3.6.1 Theoretical covariance estimate...... 70 3.6.2 Covariance estimate from Gaussian simulations...... 70 3.7 Cosmological constraints...... 71 3.8 Conclusions...... 77 3.A Theoretical prediction for CMB and weak lensing shear cross-correlation.... 79 3.B Treatment of systematic uncertainties in galaxy clustering data...... 81 3.C PSF correction and construction of weak lensing shear maps...... 81 3.D Transformation of weak lensing shear under rotation...... 85 3.E Choice of PolSpice parameter settings...... 85 3.F Correlated maps of spin-0 and spin-2 fields...... 86 3.G Validation of spherical harmonic power spectrum measurements...... 89 3.H Spherical harmonic power spectrum robustness tests...... 89 3.H.1 Comparison between spherical harmonic power spectra in equatorial and Galactic coordinates...... 91 3.H.2 Comparison between spherical harmonic power spectra derived from different foreground-reduced CMB maps...... 92 viii Contents
3.H.3 Impact of systematics correction on galaxy clustering power spectrum. 92
4 Integrated cosmological probes: extended analysis 95 4.1 Introduction...... 95 4.2 Framework for integrated probe combination...... 96 4.3 Data...... 96 4.3.1 DES weak lensing...... 97 4.3.2 CMB lensing convergence...... 99 4.3.3 Type Ia supernovae...... 100 4.3.4 Hubble parameter...... 100 4.4 Model predictions...... 100 4.5 Spherical harmonic power spectra...... 103 4.5.1 DES SV cosmic shear...... 104 4.5.2 CMB temperature and DES SV weak lensing shear cross-correlation... 105 4.5.3 CMB lensing convergence and galaxy overdensity cross-correlation... 105 4.5.4 CMB lensing convergence and CMB temperature cross-correlation... 106 4.5.5 CMB lensing convergence and SDSS Stripe 82 weak lensing shear cross- correlation...... 106 4.5.6 CMB lensing convergence and DES SV weak lensing shear cross-correlation106 4.6 Systematics...... 108 4.6.1 CMB temperature anisotropies...... 108 4.6.2 Galaxy overdensity...... 108 4.6.3 Weak lensing...... 108 4.6.4 CMB lensing convergence...... 109 4.6.5 SNe Ia...... 110 4.7 Covariance matrix...... 110 4.8 Parameter inference...... 111 4.9 Cosmological constraint results...... 114 4.10 Conclusions...... 115 4.A Impact of intrinsic alignments...... 123 4.B Impact of baryonic processes on the dark matter power spectrum...... 125 4.C Correlated spin-0 and spin-2 fields...... 127 4.D Validation of spherical harmonic power spectrum measurement...... 131 4.E Spherical harmonic power spectrum robustness tests...... 131 4.E.1 Comparison between spherical harmonic power spectra in equatorial and Galactic coordinates...... 131 4.E.2 Comparison between spherical harmonic power spectra measured from different foreground-reduced CMB temperature maps...... 134 4.F Impact of unresolved foregrounds on CMB temperature anisotropies...... 134
ix Contents
5 Integrated cosmological probes: concordance quantified 139 5.1 Introduction...... 139 5.2 Assessing concordance...... 140 5.3 Data...... 143 5.4 Results...... 144 5.5 Conclusions...... 147 5.A Gaussianity tests...... 152
6 Consistency tests in cosmology using relative entropy 155 6.1 Introduction...... 155 6.2 Data and parameter space...... 156 6.3 Data set consistency: fixed model and different data...... 157 6.3.1 Sequential experiments...... 157 6.3.2 Toy model...... 162 6.4 Model rejection: fixed data and different models...... 164 6.4.1 Implementation...... 166 6.4.2 Toy model...... 167 6.4.3 Application to cosmological data...... 169 6.5 Conclusions...... 174 6.A Assessing consistency for arbitrary updates...... 175 6.B Gaussianity tests...... 175 6.C Alternative method for computing the relative entropy...... 176 6.D Toy model implementation details...... 176 6.D.1 Implementation choices...... 176 6.D.2 Choice of fiducial polynomial coefficients...... 177 6.E SNe Ia and CMB posterior means...... 177
7 Conclusions 181
Bibliography 185
List of symbols 207
List of abbreviations 209
Acknowledgements 213
x 1 Introduction
I say, there is no darkness but ignorance. — William Shakespeare, Twelfth night (IV.II)
The last three decades have seen cosmology transit from a mostly theoretical to an experi- mental science. Recent observations have yielded precise measurements of the Large-Scale Structure (LSS) of the Universe and the Cosmic Microwave Background (CMB) radiation, leading to the establishment of the ΛCDM cosmological model.
ΛCDM is based on General Relativity (GR) as the theory of gravity and the assumption of isotropy and homogeneity of the Universe on large spatial scales (the cosmological principle). Early observations led to the discovery that the Universe is expanding and cooling down in the process, as described by the Big Bang model of the Universe. The evolution of the Universe in time is determined by its constituents: ordinary matter, as included in the standard model of particle physics; a cold dark matter component that interacts gravitationally, but has as yet not been found to interact through any other means; and a dark energy component, which is responsible for the late-time accelerated expansion of the Universe. The inhomogeneities and anisotropies detected in the CMB temperature and the galaxy distribution finally, are assumed to have been seeded by a period of exponential expansion shortly after the Big Bang, called inflation.
This simple model is described by a set of six parameters, θ {h, Ω , Ω , n , σ , τ }, where = m b s 8 reion h is the dimensionless Hubble parameter, Ωm is the fractional matter density today, Ωb is the fractional baryon density today, ns denotes the scalar spectral index of initial perturbations, σ8 1 is the r.m.s. of matter fluctuations in spheres of comoving radius 8h− Mpc and τreion denotes the optical depth to reionization.
The ΛCDM cosmological model currently provides a good fit to observations on large scales.
1 Chapter 1. Introduction
Despite its observational success, this model poses theoretical challenges. Some of the most pressing ones are that no physical theory of cosmology can explain the Universe at times earlier than the Planck time and that current theories of inflation face conceptual problems. Furthermore, despite being required by the standard cosmological model, as yet neither dark matter nor dark energy have been detected experimentally and we are far from an understanding of the latter, given the discrepancy between expectation and measurement of the value of the cosmological constant.
Several current and upcoming experiments have set out to answer these fundamental ques- tions, making the current period very exciting for cosmology. These experiments include the Dark Energy Survey (DES), the Dark Energy Spectroscopic Instrument (DESI), the Large Synop- tic Survey Telescope (LSST), Euclid, the Wide Field Infrared Survey Telescope (WFIRST) as well as Planck, the Simons Observatory and CMB Stage 4, which will cover large, overlapping parts of the sky and provide measurements at an unprecedented precision. The tightest constraints on our cosmological model will come from combining the data from all of these surveys, as these cosmological probes are complementary and correlated. It is therefore crucial to develop analysis techniques that allow us to combine the information from different cosmological probes, taking both the auto- and the cross-correlations into account. These methods will be useful for improving our fundamental understanding of cosmology and for understanding possible tensions between data sets.
In the main part of this thesis, we present a framework for integrated probe combination and describe its application to several different data sets.
This thesis is organized as follows. In Chapter1 we give a brief overview of cosmology, starting with a description of the homogeneous and the inhomogeneous Universe. We then review cosmological observations, with a particular focus on CMB observations, galaxy clustering, weak gravitational lensing and background probes. In Chapter2 we present the comparison of the constraining power of a future spectroscopic galaxy redshift survey when analyzed using the spherical harmonic tomography power spectrum and the spherical Fourier-Bessel power spectrum in a Fisher matrix analysis. In Chapters3 and4 we present an integrated framework for cosmological probe combination in which the probes are combined into a common framework starting at the map level. We describe the application of this framework to a combination of CMB temperature anisotropies, CMB lensing, galaxy clustering, weak lensing,
Supernovae Ia data and H0 measurements. In Chapter5 we investigate possible tensions between the derived constraints and other existing results within the framework of ΛCDM. Finally, in Chapter6 we revisit the relative entropy as a consistency measure and propose a novel framework for model selection based on relative entropy and posterior predictive distributions. We conclude in Chapter7.
2 1.1. The homogeneous Universe
1.1 The homogeneous Universe
On large spatial scales, the Universe appears isotropic and homogeneous, implying that perturbations around the average properties of the Universe are small. This notion of statistical isotropy and homogeneity is formalized in the cosmological principle. The discovery of the linear relation between galaxy distance and recession velocity by Edwin Hubble in 1929 [134], showed that the Universe is expanding and is therefore evolving in time. Thus the Universe can be described as homogeneous and isotropic in space, but not in time. The theoretical framework for describing the Universe as a whole is the general theory of relativity developed by Albert Einstein in 1915 [77]. In this introduction we will give a short summary of physical cosmology, mainly following Refs. [71, 27, 242] and we will give additional references where needed.
The cosmological principle within GR determines the three-dimensional spatial geome- try of the Universe. There exist only three spacetime geometries that can be foliated into isotropic and homogeneous spatial slices and these are described by the Friedmann-Lemaître-
Robertson-Walker (FLRW) metric gµν. If we choose spherical coordinates for the spatial part of the metric, the line element is given by [47]
ds2 g dxµdxν c2dt 2 a2(t)£dχ2 S2 (χ)dΩ2¤, (1.1) = µν = − + + K where c is the speed of light, a(t) is the scale factor accounting for the expansion of the Universe, χ is the comoving distance, K denotes the spatial curvature, µ,ν {0,...,3} and ∈ we assume the Einstein summation convention. The three different spatial geometries are parametrized by the quantity S and these correspond to flat (K 0), positively (K 1) or K = = negatively (K 1) curved space. We have = − 1 c p χ H sinh( K c/H ), K 1, pK 0 0 = − SK (χ) χ, K 0, (1.2) = = χ 1 c sin(p K ), K 1. p K H0 c/H0 | | | | =
The motion of particles in an FLRW spacetime is fully determined by the metric through the geodesic equation, which for a particle trajectory xµ(λ) is given by [71,27]
2 µ α β d x µ dx dx Γ , (1.3) dλ2 = − αβ dλ dλ
µ where λ parametrizes the particle path. The Christoffel symbols Γαβ are defined as
µν · ¸ µ g ∂gαν ∂gβν ∂gαβ Γ . (1.4) αβ = 2 ∂xβ + ∂xα − ∂xν
We can study the Universe by measuring the electromagnetic radiation and gravitational
3 Chapter 1. Introduction waves emitted by astronomical sources. An important consequence of the expansion of the
Universe is that radiation emitted by a distant source at time temit at a wavelength λemit will be observed at time tobs at a modified wavelength λobs. Using the FLRW metric we find [71,27]
λ a(t ) obs obs . (1.5) λemit = a(temit)
We can define the redshift z as
λ λemit a(t ) z obs − obs 1. (1.6) = λemit = a(temit) −
Setting t to today, i.e. t t , and letting a(t ) 1, a(t ) a gives obs obs = 0 obs = emit = 1 1 z . (1.7) + = a
Therefore, we can determine the value of the scale factor at emission by measuring the observed wavelength shift in known atomic transition lines.
Defining and measuring distances to astronomical objects in an expanding Universe is non- trivial. In general there exist two types of distances: comoving distances, which remain fixed through the expansion of the Universe and physical distances, which get stretched with time. The comoving coordinate distance is defined as [71, 27]
cdt Z t0 dt dχ χ c , (1.8) = a(t) =⇒ = te a(t) and it only depends on the spatial curvature through a(t) if we choose the form of the FLRW metric as in Eq. 1.1. Since photons travel on null-geodesics (ds2 0), this is equivalent to the = coordinate distance traveled by a radially-moving photon from its emission at time te to today. For t 0, Eq. 1.8 is the largest distance traveled by photons from the Big bang to today, which e = corresponds to today’s comoving horizon, i.e.
Z t0 dt χH c . (1.9) = 0 a(t)
dt The conformal time η is defined as dη /a(t) and is closely related to the comoving distance as = cdη dχ. (1.10) = If we set c 1, the conformal time and the comoving distance both have units of time or length. = The conformal time can thus be used as an alternative time coordinate, where the conformal time from t 0 to a given time t is equivalent to the comoving particle horizon. 0 = The comoving distance defined in Eq. 1.8 is not observable. If we know the intrinsic physical length of a source d or its absolute luminosity L, we can however define two indirect distance measures. The angular diameter distance is defined as the distance that relates observed
4 1.1. The homogeneous Universe angular size ∆θ to known physical extent d, i.e. [71,27]
d SK (χ(z)) dA(z) . (1.11) = ∆θ = (1 z) + The luminosity distance of a distant emitter is defined as the distance relating observed flux F to intrinsic luminosity L of the source, i.e. [71, 27]
s L dL(z) SK (χ(z))(1 z). (1.12) = 4πF = +
The dynamics of the Universe is determined by GR. The Einstein field equations (EFE) of GR relate the curvature of four-dimensional spacetime to the matter content of a given physical system and are given by [47]
8πG Gµν Λgµν Tµν, (1.13) + = c4 where G denotes the gravitational constant. Gµν is a measure for spacetime curvature and is called Einstein tensor; Tµν is the stress-energy tensor, which is a measure for the Universe’s energy density content, and Λ is the cosmological constant. For an ideal fluid with density ρ, pressure p and four-velocity U µ, it is given by ³ p ´ Tµν ρ UµUν pgµν (1.14) = c2 + +
The evolution equation for the scale factor a(t) can be obtained by requiring the FLRW metric to be a solution to the EFE for an ideal fluid with energy density ρ and pressure p. This gives us two evolution equations for a(t), called the Friedmann equations [71,27, 242,47]
µ a˙ ¶2 8πG K c2 H 2(t) ρ , (1.15) = a = 3 − a2 a¨ 4πG µ 3p ¶ ρ , (1.16) a = − 3 + c2 where denotes derivative with respect to time t and we have defined the Hubble parameter · a˙ H(t) /a. = Equations 1.15 and 1.16 show that the dynamics of the Universe’s expansion depends on its constituents. Cosmological observations to date suggest that the Universe’s energy density is distributed between several different species, which can be modeled as perfect fluids. Due to energy-momentum conservation, the energy density of these species evolves differently as the Universe expands. We can introduce an equation of state parameter w, which relates energy density and pressure of an ideal fluid, i.e. [47]
p wρc2. (1.17) =
5 Chapter 1. Introduction
The conservation of energy-momentum then implies that [47]
3(1 w) ρ a− + . (1.18) ∝ The most important component for the early-time evolution of the Universe are relativistic species such as photons and massless neutrinos. These have an equation of state parameter 1 4 w and therefore their energy density falls as ρ a− as the Universe expands, making = 3 r ∝ them dynamically important at early times. At intermediate times in the Universe’s history, the energy density is dominated by non-relativistic species, such as baryons. Starting as early as the 1930s, observations showed evidence for an additional matter component in the Universe, which was subsequently called dark matter [286]. As yet, dark matter has only been detected through its gravitational interactions. It is thought to be collisionless and to interact at most very weakly with ordinary matter (e.g. [59, 55, 73, 185, 230]). Both ordinary and dark matter have an equation of state parameter w 0 and thus only dilute 3 = as ρ a− . The discovery of the accelerated expansion of the Universe in the late 1990s m ∝ [225, 205] has shown that the observations can only be explained by invoking an additional component. This component is responsible for the late-time acceleration of the Universe and 1 is called dark energy. From Equations 1.16 and 1.17 we see that a¨ 0 requires w /3, i.e. a > < so-called dark energy component with negative pressure. The case w 1 is special, as it = − gives constant energy density, i.e. ρ const. (c.f. Eq. 1.18). The energy-momentum tensor for = such a component follows from Eq. 1.14 as T pg , which has the same form as both the µν = µν energy-momentum tensor of the vacuum and the cosmological constant Λ. Therefore, the cosmological constant can be interpreted as due to the energy density of the vacuum [47].
The energy density scale factor scalings derived above allow us to rewrite the Friedmann equation. Let us define the critical density of the Universe as [47]
3H 2 ρcrit . (1.19) = 8πG
The ratio between the total density of the Universe and the critical density determines the spatial geometry of the Universe; we have [47]:
ρ ρ : K 0 (1.20) > crit > ρ ρ : K 0 (1.21) = crit = ρ ρ : K 0. (1.22) < crit <
We can finally define the density parameter Ωi for each species including the curvature as [47]
ρi Ωi . (1.23) = ρcrit
6 1.2. The inhomogeneous Universe
Using these definitions we can finally rewrite Eq 1.15 as [47]
2 H 4 3 2 Ω a− Ω a− Ω Ω a− , (1.24) 2 = r + m + Λ + K H0 where H0 denotes the present day Hubble parameter. Observations to date show that the Universe is spatially very close to flat [212], in the following discussion we will therefore set the spatial curvature to zero and assume a flat Universe.
1.2 The inhomogeneous Universe
1.2.1 Linear perturbation theory
The observed Universe is not homogeneous. At small scales we see structures formed through the action of gravity such as galaxy clusters, galaxies, stars and planets. In order to understand how these formed, we need to study the generation and evolution of perturbations to these deviations from isotropy and homogeneity. As the Universe is smooth on large scales, we can treat these deviations as perturbations to an otherwise smooth FLRW background g¯µν, i.e. [279, 27, 242]
g g¯ h , (1.25) µν = µν + µν where h is a symmetric four-tensor and h g¯ . Of the ten independent components µν | µν| ¿ µν of hµν, only six constitute physical degrees of freedom, as the rest is eliminated by gauge freedom.
The perturbations hµν can be decomposed in a scalar, vector and tensor part. This decompo- sition is helpful, as the equations for scalar, vector and tensor perturbations decouple in linear theory. Focusing only on the scalar perturbations, which are the relevant ones for our work, and adopting the Newtonian gauge we can write the perturbed line element as [71]
µ 2Ψ(x,t)¶ µ 2Φ(x,t)¶ ds2 c2 1 dt 2 a2(t) 1 £dχ2 χ2dΩ2¤, (1.26) = − + c2 + + c2 + where the functions Ψ(x,t) and Φ(x,t) are the only two scalar metric perturbations, called cur- vature and Newtonian potential respectively. Furthermore, the quantity dΩ dθ2 sin2 θdφ2 = + denotes the angular part of the metric. Analogously, we can define small perturbations around the mean properties of the Universe’s components; for the particle density for example, we define:
ρi (x,t) ρ¯i (x,t) δi (x,t) − . (1.27) = ρ¯i (x,t)
The metric perturbations are coupled to the perturbations in the particle distributions through gravity. In order to obtain evolution equations for these quantities we need to split the problem
7 Chapter 1. Introduction into two different parts: First, we need to consider the effect of metric perturbations on the distribution of matter and radiation and second, the effect of perturbations to the energy density of the Universe on the metric perturbations. The first part amounts to deriving the Boltzmann equations for the particle distributions, while for the second part we need to derive the EFE for the perturbed FLRW metric. This gives rise to a set of coupled differential equations, called Einstein-Boltzmann equations. In deriving these equations, we can make several simplifications: Firstly, we can work at linear order and drop all higher order terms, as we assume all perturbations to be small. Secondly, we can transform from real to Fourier space in order to remove derivative terms and thus significantly simplify the equations. Working in Fourier space has the additional advantage that different Fourier modes evolve independently in the linear regime.
We adopt the following Fourier convention in our derivation:
Z d3k f (x) eikx f (k), (1.28) = (2π)3 where f is an arbitrary function on R3. With this convention, the Boltzmann equations for radiation, dark matter and baryons finally become [71]
· ¸ Φ˙ ikµΨ µvb 1 Θ˙ ikcµΘ τ˙ Θ0 Θ P2(µ)Π , (1.29) + = −c2 − c − − + c − 2 Π Θ Θ Θ , (1.30) = 2 + P2 + P0 · 1 ¸ Θ˙ P ikµΘP τ˙ ΘP (1 P2(µ))Π , (1.31) + = − − + 2 − 3Φ˙ δ˙ ikv , (1.32) + = − c2 a˙ v˙ v ikΨ, (1.33) + a = − 3Φ˙ δ˙ ikv , (1.34) b + b = − c2 a˙ τ˙ v˙ ikΨ [v 3icΘ1], (1.35) b + a = − + R b + Φ˙ ikµΨ Θ˙ ν ikcµΘν , (1.36) + = −c2 − c where we have expressed the perturbation variables in terms of conformal time and thus · denotes derivative w.r.t. η. Furthermore, k is the wave vector, µ is the cosine of the angle between the photon momentum and the wave vector, P` is the Legendre polynomial of order `, Θ denotes the perturbations to the photon temperature (assuming an Einstein-Boltzmann distribution at zeroth order), Θi are the multipole moments of the temperature perturbation field, τ is the optical depth, ΘP denotes the perturbations to the radiation polarization field and ΘPi its multipole moments, δ is the dark matter overdensity field, v its velocity, δb is the
3ρb,0 baryon overdensity field, v its velocity, R is the baryon-to-photon ratio, R /4ργ,0, and finally b = Θν denotes the perturbations to the neutrino temperature.
8 1.2. The inhomogeneous Universe
Eq. 1.29 governs the evolution of the perturbations to the photon temperature and includes effects such as Compton scattering between photons and baryons. Eq. 1.31 describes the generation of photon polarization, mainly from the quadrupole moment of the temperature field. Equations 1.32 to 1.35 describe the evolution of dark and baryonic matter perturbations and the associated velocity fields. The dark matter velocity is only affected by gravitational forces, while the baryon velocity field is additionally modified by Compton scattering with photons. Finally Eq. 1.36 describes the evolution of neutrino temperature perturbations, assuming a vanishing neutrino mass.
This set of equations needs to be closed by the Einstein equations describing the effect of perturbations in the Universe’s energy density on the metric perturbations. The linearized EFE give us [71] µ ¶ 2 a˙ a˙ 2 £ ¤ k Φ 3 Φ˙ Ψ 4πGa ρmδm 4ρrΘr,0 , (1.37) + a − a = + k2 (Φ Ψ) 32πGa2ρ Θ , (1.38) + = − r r,2 where the subscript m includes all matter components and r all radiation components. From Eq. 1.38 we see that within GR and in absence of quadrupole moments in the radiation distribution (the so-called anisotropic stress), we have Φ Ψ. This is a very good assumption = − during matter domination [33] and we therefore set Φ Ψ for the remainder of this thesis. = − In general, the set of coupled Einstein-Boltzmann equations needs to be solved numerically and several codes exist and are publicly available (e.g. [243, 167, 36, 221]). To get more intuition about the solutions of these equations, we will briefly discuss the solutions of these equations qualitatively for matter and radiation perturbations.
1.2.2 Cosmological structure formation
The physics of cosmological structure formation can be pictured as an interplay between gravity, the expansion of the Universe and thermal motion/pressure of the particles. As matter is essentially pressureless, while radiation has high pressure, perturbations to these two components evolve differently, as discussed below.
Matter perturbations
The evolution of matter perturbations is mainly determined by the properties of dark matter, as it makes up most of the matter in the Universe. In general, the evolution of a perturbation at scale k can be divided into different phases: At early times, all perturbations will lie outside the horizon. As long as the scale k of a perturbation is larger than the horizon, no causal physics can affect the mode. At intermediate times, two things happen: (i) As the size of the horizon grows, it will reach the size of the perturbation and the perturbation will cross the horizon. When this happens the mode can be affected by causal physics. (ii) The Universe
9 Chapter 1. Introduction transitions from radiation to matter domination. Since gravitational potential perturbations evolve differently during radiation and matter domination, the properties of a given mode depend on horizon-crossing taking place during radiation or matter domination.
On super-horizon scales, matter perturbations do not evolve and stay constant, except for a small change in amplitude during the transition from radiation to matter domination. As the size of the horizon grows, small-scale perturbations enter the horizon, followed by large- scale perturbations. Small-scale perturbations tend to cross the horizon during radiation domination, while perturbation on large-scales only enter during matter domination. The causal growth of matter perturbations during radiation domination is not very efficient as the modes only grow logarithmically with time. This is due to the oscillation of gravitational potentials during radiation domination, as we will see below. When the Universe enters matter domination, the gravitational potentials stay constant. On scales significantly smaller than the horizon, perturbations to the matter distribution therefore start growing linearly with scale factor a(t), independently of their scale k. This suggests that the amplitude of matter perturbations at scale k and late times a, δ(k,a), can be written as
δ(k,a) k2Φ (k)T (k)D (a). (1.39) ∝ p 1
The quantity Φp denotes the primordial potential perturbation amplitude, which we will discuss below. The transfer function T (k) accounts for all scale-dependent processes during horizon crossing and radiation domination. Finally the linear growth factor D1(a) accounts for scale-independent growth during matter domination.
Even though baryons make up only a small fraction of the total matter in the Universe, they have a characteristic impact on the matter distribution, which is due to electromagnetic interactions. We will discuss these Baryonic Acoustic Oscillations (BAO) in Chap. 1.3.3.
Radiation perturbations
The evolution of photon perturbations can be divided into two different phases: the evolution before recombination at z 1100 and after recombination. Before recombination, photons ? ≈ and baryons are tightly coupled through Compton scattering and can be treated as a single photon-baryon fluid [132, 71]. During recombination, electrons and protons combine to form neutral hydrogen and the fraction of free electrons falls significantly. Photons thus decouple from baryons and essentially start free-streaming through the Universe. Even though decoupling and last scattering take place during recombination, which lasts a finite time, here we will assume that recombination occurs instantaneously and thus use last scattering and recombination interchangeably.
In order to better understand the physics of radiation anisotropies, it is instructive to study the perturbations in configuration rather than Fourier space [26]. After horizon-crossing overdense regions will start to gravitationally collapse. As the photon density increases, so
10 1.2. The inhomogeneous Universe does its pressure and as soon as the pressure dominates, it will send out a spherical acoustic sound wave in the photon-baryon fluid. This sound wave will travel until recombination. At recombination, photons and baryons decouple and the photons will thus start free-streaming through the Universe, leaving the baryons behind. This process imprints the distance trav- eled by the photon-baryon fluid until recombination rs(η?) as a characteristic scale on the distribution of radiation perturbations. In configuration space, this gives rise to enhanced correlations between photon properties separated by a distance rs(η?). In Fourier space on the other hand, this leads to oscillations in the radiation distribution. Denoting the sound speed of the fluid as cs(η), we have [71]
Z η? rs(η?) cs(η0)dη0, (1.40) = 0 where the subscript ? denotes evaluation at recombination.
There are several additions to the simplified picture given above. Firstly, even during recombi- nation, photons and baryons are not perfectly coupled, which means that photons travel a finite distance between different scattering events. The total distance traveled by photons in a Hubble time, called diffusion length, is determined by their mean free path and the average number of scattering events. Any perturbation on scales smaller than the photon diffusion length will therefore be damped. Secondly the photon distribution is affected by several late-time phenomena after recombination, such as rescattering with free electrons after reion- ization, temporal changes in gravitational potentials due to dark energy, and gravitational lensing.
Initial conditions
The study of cosmological structure formation immediately raises the question about the origin of primordial fluctuations. The theory of inflation, which was developed in the early 1980s [104, 172, 16], can both explain the flatness and the curvature problem of cosmology as well as provide a mechanism for the generation of primordial perturbations. However, this theory suffers from several, currently unresolved conceptual problems, such as eternal inflation (multiverse) or the trans-Planckian problem (see e.g. [278, 40]).
The theory of inflation postulates that the Universe underwent an early phase of exponential expansion driven by the potential energy of a scalar field with negative equation of state parameter. Such an inflationary phase brings the Universe to spatial flatness and causes the observable Universe to arise from a causally connected region. Furthermore the exponential expansion will stretch the ubiquitous quantum mechanical fluctuations in the scalar field to observable scales, thus providing a natural mechanism for the generation of primordial per- turbations and the initial conditions for structure formation. Due to the quantum mechanical nature of the initial conditions, we cannot theoretically predict the value of the perturbation field but only its statistical properties. The primordial fluctuations generated by single-field inflationary models are Gaussian to a very good approximation [179, 10] and can therefore be
11 Chapter 1. Introduction fully described by their two-point function or power spectrum1. Focusing on such models, the predicted primordial potential power spectrum can be written as [71]
2 µ ¶ns 1 µ ¶2 50π ck − 2 Ωm PΦ (k) δ , (1.41) p = 9k3 H H D (a 1) 0 1 = where δH denotes the scalar power spectrum amplitude. The primordial power spectrum thus forms the seeds of all cosmological structures we see today.
1.3 Observational probes
In cosmology, our laboratory is the Universe, which we can observe without any control on the experimental conditions. We can make inference about the Universe by observing its properties, the cosmological probes, and comparing them to theoretical predictions. Since we can only make theoretical predictions for statistical properties, most cosmological analyses focus on the two-point function or power spectrum of the observables. If the observable is a Gaussian random field, these statistics capture all of the information. If the observable has non-Gaussian features, it can still be analyzed in terms of the two-point function but a significant amount of information will only be accessible through higher-order statistics.
There is a fundamental measurement uncertainty associated with cosmological observations that is due to the fact that we can only observe our particular realization of the Universe and therefore cannot repeat an experiment. This uncertainty is called the cosmic variance [71].
1.3.1 Statistics of random fields
As a consequence of the cosmological principle (i.e. the symmetry of the background cos- mology), the Gaussian primordial fluctuations are statistically isotropic and homogeneous. Linear evolution of perturbation amplitudes preserves these properties, whereas non-linear evolution creates significant non-Gaussianity and can break the isotropy and homogeneity of the Universe. In the work presented in Chapters2-5 we mostly focus on Gaussian fields that can be fully described in terms of two-point correlators. In this section, we thus review the statistical description of random fields in terms of two-point functions, mostly following Ref. [51].
Let f (x) be a random field on R3. The three-dimensional correlation function (two-point function) is defined as the ensemble average of the product of the value of the field at two different points x and x0 x r, i.e. [51] = −
ξ(x,x0) f (x)f (x0) , (1.42) = 〈 〉 where denotes the average over all realizations. The quantity ξ is thus a measure for the 〈···〉
1These two quantities will be discussed in more detail in Sec. 1.3.1.
12 1.3. Observational probes correlation between the field values at different positions. If f is statistically isotropic and homogeneous, i.e. its statistical properties are invariant under translations and rotations, then the correlation function only depends on the norm of the separation of the two points i.e. [51]
ξ(x,x0) ξ( x x0 ) ξ(r ) f (x)f (x0) , (1.43) = | − | = = 〈 〉 where r r . Since cosmological perturbation theory generally predicts Fourier-space quanti- = | | ties, it is instructive to consider the analogous relations in Fourier space. Using Eq. 1.28, the three-dimensional power spectrum of f is defined as [51]
P(k,k0) f (k)f ∗(k0) , (1.44) = 〈 〉 where the superscript denotes complex conjugation. For statistically isotropic and homoge- ∗ neous fields, this relation simplifies and we obtain [51]
f (k)f ∗(k0) P(k)δ(k k0), (1.45) 〈 〉 = − i.e. different Fourier modes are uncorrelated and the power spectrum only depends on the modulus of the wave vector k.
In order to take into account the spherical sky geometry, we can perform analogous steps for fields f (xˆ) defined on the two-sphere S2. The two-point correlation function of f is defined as [51]
ξ(θ,θ0) f (θ)f (θ0) , (1.46) = 〈 〉 2 where θ,θ0 S . As before, the correlation function of an isotropic field only depends on the ∈ angular separation of the points θ,θ0, i.e. we have [51]
ξ(∆θ) f (θ)f (θ0) , (1.47) = 〈 〉 where ∆θ θ θ0 . Any function defined on the sphere can be decomposed into spherical = | − | harmonics as
X X f (θ) f`mY`m(θ), (1.48) = ` m ` | |≤ where Y (θ) denotes the spherical harmonic of degree ` and order m, m `. The quantities `m | | ≤ f`m are the spherical harmonic coefficients of f . For a statistically isotropic field, different harmonic modes are uncorrelated and the spherical harmonic power spectrum (or angular power spectrum) C` of f is then defined by [51]
f`m f ∗ C`δ`` δmm . (1.49) 〈 `0m0 〉 = 0 0 This is the Fourier-space equivalent of Eq. 1.47 for random fields defined on the sphere.
13 Chapter 1. Introduction
1.3.2 Cosmic Microwave Background
The Cosmic Microwave Background is the most constraining cosmological probe to date and it consists of the photons that last scattered with electrons before recombination and have been free-streaming through the Universe since. This relic radiation has been discovered by Penzias and Wilson in 1965 [203] and found to be nearly isotropic. In 1992, the COBE DMR satellite measured its blackbody spectrum to high precision and detected fluctuations in the ∆T 5 CMB temperature at the level of /T 10− [252]. More recently, the WMAP [30, 124] and later ∼ the Planck [210, 215, 212] satellite produced high-resolution maps of the distribution of CMB temperature and polarization anisotropies on the sky and measured its statistical properties.
The CMB is an inherently two-dimensional quantity and can thus be analyzed in terms of the spherical harmonic power spectrum C`. Theoretically, the CMB power spectrum is well described within linear theory, as the amplitude of radiation perturbations is small at all times. These perturbations thus evolve linearly from the initial fluctuations, preserve Gaussianity and can thus be characterized through the angular power spectrum, which is mostly sensitive to early Universe physics. We can separate the direction- and initial condition dependence of the CMB temperature anisotropies T and relate their spherical harmonic power spectrum to the primordial power spectrum generated during inflation and the evolution of matter perturbations as [71]
Z ¯ ¯2 2 ¯∆T`(k)¯ C TT dk k2P lin(k)¯ ¯ , (1.50) ` = π δδ ¯ δ(k) ¯ where ∆T` denotes the transfer function of the temperature anisotropies and δ is the dark matter overdensity. The CMB temperature anisotropy power spectrum is mainly sensitive to the spatial curvature density of the Universe, Ω , the physical baryon density ω Ω h2, the K b = b physical matter density ω Ω h2 and the scalar spectral index n . m = m s The observed CMB temperature anisotropies receive contributions not only from early uni- verse (i.e. pre-recombination) physics but also from late-time phenomena. These secondary CMB anisotropies largely fall into two classes, gravitational and scattering phenomena [132]. Scattering secondaries are due to the fact that galaxy formation results in free electrons that rescatter with CMB photons at late times. An example is the Sunyaev-Zel’dovich effect [260, 261], which is due to scattering of CMB photons with free electrons with large bulk or thermal velocities. An example for a gravitational effect is the integrated Sachs-Wolfe (ISW) effect. During matter-domination, the amplitude of gravitational potential perturbations stays constant while it starts to decay as soon as dark energy becomes significant for the expansion. The ISW is due to the blue- or redshifting of CMB photons as they travel through these decaying large-scale gravitational potential perturbations and thus provides a direct evidence for accelerated expansion.
The most prominent gravitational effect on the CMB however, is the lensing of CMB photons by the LSS. Analogously to galaxy weak lensing, the CMB photons get deflected by the cumulative
14 1.3. Observational probes lensing effect of the LSS. Therefore the observed CMB temperature at position θˆ can be obtained from the unlensed temperature as [206, 166]
T (θˆ) T (θˆ φ(θˆ)) T (θˆ) i φ(θˆ) T (θˆ) O(φ2), (1.51) obs = unl + ∇ = unl + ∇ ∇i obs + where denotes the two-dimensional gradient on the sphere. The quantity φ(θˆ), which gives ∇ rise the deflection, is the CMB lensing potential and is given by [206, 166]
Z 2 dz χ? χ(z) φ(θˆ) − Ψ(χ(z)θˆ,z), (1.52) = −c H(z) χ?χ(z) where χ? denotes the comoving distance to the last scattering surface. The CMB lensing convergence is defined as κ(θˆ) 1 2φ(θˆ)[35]. As we know that the statistical properties of = − 2 ∇ the unlensed CMB are those of a Gaussian random field, this means that [206]
(i) if we average over realizations of the CMB temperature for fixed φ(θˆ), CMB lensing generates anisotropies in the CMB temperature covariance matrix,
(ii) if we average over the CMB temperature anisotropies and φ(θˆ), CMB lensing leads to a non-vanishing four-point function in the observed CMB data.
These two signatures allow us to reconstruct the CMB lensing potential from CMB temperature anisotropy maps and obtain information from the lensing convergence. This is in contrast to galaxy weak lensing. Since we do not have any information about the intrinsic distribution of galaxies and their initial positions, we cannot reconstruct the deflection field but need to focus on cosmic shear, which correlates the shapes of nearby galaxies [166].
The CMB lensing potential can be reconstructed from the CMB temperature and polarization field using quadratic estimators (e.g. [195]). It is sensitive to the integrated matter distribution from the last scattering surface to today and the CMB lensing kernel peaks at a redshift z 2. ∼ In order to obtain cosmological information from the CMB lensing potential, we can measure its angular power spectrum or the cross-correlation between the CMB lensing potential and any other tracer of the LSS.
CMB lensing together with secondary CMB temperature anisotropies in general, make the CMB a powerful probe of not only the high- but also the low redshift Universe. Finally, the recent progress in CMB polarization measurements opens up a new window for CMB science.
1.3.3 Galaxy clustering
While the primary CMB anisotropies provide a picture of the matter distribution during the epoch of last scattering (z 1100), a wealth of complementary cosmological information is ? ∼ contained in the three-dimensional matter distribution in the Universe at lower redshift. This information can be extracted by measuring the statistical properties of the matter overdensity
15 Chapter 1. Introduction
field δm, e.g. through the correlation function ξmm(r ) in real space or the power spectrum Pmm(k,z) in Fourier space.
Observations have shown that dark matter (DM) makes up approximately 80% of the Universe’s matter content, whereas baryons only constitute the remaining 20% (e.g. [124, 212]). However, as DM does not or interacts very weakly, we need to infer the matter distribution in the Universe using the distribution of the luminous matter, such as galaxies or diffuse gas. In this introduction, we focus on galaxies as a tracer of the underlying DM distribution.
Galaxies form at the density maxima of the matter field and are therefore biased tracers of the overall matter distribution [24]. Furthermore, as galaxies are discrete tracers of an underlying continuous density field, the measured galaxy density receives a contribution due to shot noise. We can thus write the observed galaxy density as [201]
δ (k,z) b(k,z)δ(k,z) n(k,z), (1.53) g = + where b(k,z) is called galaxy bias parameter and parametrizes the relation between galaxy position and the underlying dark matter field. The quantity n(k,z) accounts for shot noise.
The galaxy overdensity field is three-dimensional and optimal extraction of cosmological information from galaxy clustering thus requires an analysis in 3D. Measuring galaxy redshifts to high precision however is challenging and requires spectroscopic galaxy redshift surveys, which infer galaxy redshifts by measuring the electromagnetic spectrum of each galaxy and identifying known atomic transition lines. These measurements are time-consuming and spectroscopic surveys thus typically contain significantly less galaxies than photometric sur- veys, which infer galaxy redshifts by measuring galaxy fluxes in a given number of photometric passbands. Since photometrically determined galaxy redshifts have large uncertainties, galax- ies in these surveys are typically analyzed by splitting the sample in redshift and measuring the two-dimensional clustering in each of these bins. Spectroscopic galaxy redshift surveys on the other hand are typically analyzed using the power spectrum or the three-dimensional correlation function.
We can compute the three-dimensional galaxy clustering power spectrum from Eq. 1.53. Noting that δ(k,z) and n(k,z) are uncorrelated, we obtain the power spectrum at a given redshift z as
2 δ (k)δ∗(k0) P (k,z)δ(k k0) b (k,z)P(k,z)δ(k k0) P (k,z)δ(k k0), (1.54) 〈 g g 〉 = gg − = − + SN − where PSN(k,z) denotes the shot noise power spectrum. The galaxy power spectrum and the correlation function can be interpreted as the probability of detecting two galaxies separated by a distance r [201]. In contrast to the CMB, the matter distribution, including galaxies and clusters, shows large perturbations to the underlying FLRW background and therefore, its statistical properties cannot generally be modeled using linear theory. As large scales have not had the time to collapse yet, these non-linear effects are however mostly confined to small
16 1.3. Observational probes scales. Non-linear evolution has generally two effects: first the non-linear power spectrum deviates from the linear one and second a significant amount of cosmological information is contained in non-Gaussian statistics of the matter-density field, such as the bispectrum. Several approaches to modeling the non-linear matter power spectrum have been developed in the literature. These range from analytic perturbation theory approaches (for a review see [31]) to fitting functions calibrated through numerical simulations (e.g. [251, 267, 186, 187]).
Returning to Eq. 1.54, the shot noise power spectrum for Poisson samples of the underlying density field is given by
1 PSN(k,z) , (1.55) = n¯ where n¯ is the mean galaxy number density of the sample. In practice however we need to take into account modifications to the shot noise power spectrum due to survey masks or non-Poissonian nature of the tracers (see e.g. Chapter3,[199]).
In our discussion so far we have ignored the baryons since they are subdominant with re- spect to the dark matter. However, the sound waves in the photon-baryon fluid prior to recombination leave characteristic imprints on the matter distribution. Let us recall that after recombination, the matter distribution consists of a superposition of spherical dark matter overdensities surrounded by rings of baryons at a distance of rs(η?). As structures grow, baryons will fall into the DM potential wells and equally the DM will fall into the potentials generated by the baryon distribution. This results in an enhanced probability of finding two galaxies separated by a distance rs(η?), which manifests itself as a peak in the galaxy correla- tion function and corresponding oscillations in the power spectrum, called baryonic acoustic oscillations (BAO).
As cosmological constraints from measurements of the shape of Pgg(k,z) are generally weaker than those from the CMB, measurements of the 3D power spectrum are mainly used to extract specific features in the power spectrum, such as the BAOs or redshift space distortions (RSD). The BAO scale, which is imprinted on the galaxy distribution, can be used as a standard ruler to infer the angular diameter distance dA(z) (Eq. 1.11) to an effective redshift. RSDs arise due to peculiar galaxy velocities. When estimating the distance of galaxies from their redshift we make the assumption that it is solely due to the Hubble flow. However, galaxies tend to fall into overdense regions and these peculiar velocities give rise to an anisotropy in the galaxy clustering signal along the line-of-sight. This allows for the measurement of the linear growth rate f (z)σ8(z), which is given by:
dlogD1(a) f (z)σ8(z) a (1 z) 1 σ8(z). (1.56) = dloga | = + −
In photometric data, these features are damped due to galaxy redshift uncertainties and thus the 2D power spectrum is usually used to derive constraints from the power spectrum shape.
17 Chapter 1. Introduction
1.3.4 Cosmic shear
Clustering studies treat galaxies as point-particles and measure the cosmological information contained in the statistical properties of their relative positions. Galaxies however are not point-like but rather extended objects with a variety of shapes and forms. As photons emitted by distant galaxies travel towards us, their paths get deflected by the gravitational effect of the intervening matter, causing galaxy shapes to appear distorted. This effect is called gravitational lensing. If the distortions are large, i.e. when photons are deflected by a large overdensity such as a galaxy cluster or a massive galaxy, we speak of strong gravitational lensing. Weak gravitational lensing on the other hand, describes the tiny but coherent distortions caused by the cumulative lensing effect of the LSS between the distant galaxy and the observer. These two regimes of gravitational lensing carry complementary information: strong lensing predominantly probes the matter distribution in the lensing galaxy and therefore the small- scale distribution of matter. Weak lensing or cosmic shear is sensitive to the integrated distribution of gravitationally interacting matter and it thus probes the large-scale distribution of both luminous and dark matter in the Universe. In the following, we will focus on weak gravitational lensing only, as it is one of the cosmological probes used in the work presented later in this thesis.
Weak lensing formalism
In order to derive the light deflections caused by the LSS, let us look at the situation of a source galaxy located at the true angular position θS as seen from the observer but observed under the angle θ due to the deflection of light (see Fig. 1.1). Since gravitational lensing conserves surface brightness we have [71]
I (θ) I (θ ), (1.57) obs = true S where Iobs(θ) is the observed specific intensity distribution and Itrue(θS) denotes the surface brightness at the source plane. To derive the relation between the observed θ and true angular positions θS, we need to solve the geodesic equation for photons (Eq. 1.3) in the presence of inhomogeneities along the photon path. This gives us [71]
Z χ µ ¶ i i 2 ∂Φ(x(χ0)) χ0 θS θ 2 dχ0 1 , (1.58) = + c 0 ∂xi − χ where i 1,2 and Φ denotes the gravitational potential perturbation. In the weak-field regime, = we can linearize this equation and define the transformation matrix Ai j , which defines a linear mapping between the angular positions on the source plane and the observed positions, as
i à ! ∂θS 1 κ γ1 γ2 Ai j − − − , (1.59) = ∂θ j = γ 1 κ γ − 2 − + 1
18 1.3. Observational probes
Iobs(✓)
Itrue(✓S)
✓ ✓ O S
Figure 1.1: Schematic illustration of weak gravitational lensing.
where κ denotes the scalar convergence and γ (γ ,γ ) denotes the spin-two shear. The = 1 2 convergence affects both components of the source’s angular positions equally and it thus describes isotropic magnification or demagnification of the source image. The shear on the other hand describes anisotropic stretching of the galaxy shape. As we have no information on intrinsic galaxy sizes, weak gravitational lensing studies focus on measuring cosmic shear γ, which is given by
A11 A22 γ1 − , (1.60) = − 2 γ A . (1.61) 2 = − 12
We can relate the shear to the gravitational potential perturbations along the line of sight by combining Equations 1.58 and 1.59: From Eq. 1.58 we see that the deflection is sensitive to the gravitational potential integrated along the perturbed photon path x(χ0). We can simplify this by making use of the fact that perturbations in the weak-field regime are small and set the perturbed photon path equal to the unperturbed one, i.e. x(χ0) χ0θ [155]. Using this = so-called Born approximation, we finally obtain [71]
à ! Z χ 2 µ ¶ κ γ1 γ2 2 ∂ Φ(x(χ0)) χ0 Ai j δi j − − − 2 dχ0 χ0 1 . (1.62) − = γ κ γ = c 0 ∂xi ∂x j − χ − 2 − + 1 Measuring the statistical properties of the cosmic shear field thus allows us to infer the statisti- cal properties of the gravitational potential perturbations.
The shear is a spin-2 field, i.e. it is invariant under rotations of 180◦ and can be decomposed into spin-2 spherical harmonics. In order to compute its spin-0 angular power spectrum, we
19 Chapter 1. Introduction need to express the information contained in γ in terms of scalar variables. Any spin-2 field on S2 can be decomposed into two scalar fields, the curl-free E-mode and the B-mode field. Since cosmic shear is sourced by a single scalar field in the Born approximation, the B-mode vanishes and we have [25, 155]:
C EE C γγ C κκ, (1.63) ` ≡ ` = ` C BB 0. (1.64) ` = Using the Poisson equation (late-time limit of Eq. 1.37) to relate the potential to the matter overdensity we obtain the cosmic shear E-mode power spectrum for a source at redshift z as ∗ [71]
à !2 2 Z 2 µ ¶2 µ 1 ¶ γγ 3 ΩmH0 c χ χ χ nl ` /2 C dz ∗ − P k + ,z , (1.65) ` = 2 c2 H(z) a2 χ δδ = χ(z) ∗ where we have used the Limber approximation to simplify the integral [169, 147, 148]. In weak gravitational lensing surveys, we usually measure the power spectrum by averaging over a large number of source galaxies with redshift distribution n(z), defined as Z dz n(z) 1. (1.66) =
For such a galaxy population, Eq. 1.65 becomes [71]
γ ¡ ¢ γ ¡ ¢ µ 1 ¶ Z c W χ(z) W χ(z) ` /2 C γγ dz P nl k + ,z , (1.67) ` = H(z) χ2(z) δδ = χ(z) where we have defined the weak lensing shear window function as
2 Z χh 3 ΩmH χ(z) χ(z0) χ(z) γ ¡ ¢ 0 − W χ(z) 2 dz0n(z0) . (1.68) = 2 c a χ(z) χ(z0)
We can see from Eq. 1.67 that the cosmic shear power spectrum is indeed sensitive to the matter power spectrum and it thus too becomes increasingly non-linear at small scales. In addition, numerical simulations have shown that the dark matter power spectrum on small scales is significantly affected by baryonic processes such as gas cooling and feedback from Active Galactic Nuclei (AGN) and supernovae (e.g. [275]). Both these effects need to be taken into account when modeling the cosmic shear power spectrum.
Galaxy ellipticities
Following the discussion above, we expect weak gravitational lensing to cause the shapes of galaxies to appear distorted and the shear can thus be inferred from the observed galaxy ellipticities i.e. shapes. This is challenging since galaxies show a variety of intrinsic ellipticities 1 with an r.m.s. of the order of ²2 /2 0.3 whereas the cosmic shear distortions to these 〈 〉 ∼
20 1.3. Observational probes intrinsic shapes are on the level of a few percent, i.e. γ 0.03 [155]. Therefore, cosmic shear ∼ cannot be detected on single galaxies and weak lensing measurements need to be performed on large galaxy samples in order to overcome the statistical noise. The intrinsic shapes of galaxies have additionally been found to exhibit significant correlations, which are called intrinsic alignments [120, 180, 128, 80, 196, 140, 247] and are due to the tidal gravitational field surrounding the galaxy. Weak gravitational lensing assumes all correlations between galaxy shapes to be due to gravitational lensing and therefore intrinsic alignments need to be taken into account in cosmic shear analyses in order to not bias cosmological parameter constraints (see e.g. [159]).
Several different ellipticity estimators have been developed in the literature (e.g. [149, 188, 285, 245, 246]) and they mainly fall into two categories: moment-based or model-based methods. Moment-based methods infer the galaxy ellipticity from the moments of the galaxy light profile, while model-based methods infer ellipticities by fitting a model galaxy to the image. All these methods need to operate in the low signal-to-noise regime and most of them are affected by biases, which need to be calibrated using simulations or estimated theoretically.
1.3.5 Background probes
All the cosmological probes discussed so far measure the cosmological information contained in the inhomogeneities of the Universe. In addition to these, there exist several cosmological probes that measure the homogeneous expansion history of the Universe and thus derive cosmological constraints. Examples include BAOs, which have been discussed in Chap. 1.3.3, Type Ia Supernovae and direct Hubble constant measurements. In the following, we focus on the latter two as these are the most important ones for this work.
Type Ia Supernovae
Type Ia Supernovae are thought to be the result of thermonuclear explosions triggered by mass accretion onto white dwarf stars. Since the electron-degeneracy pressure inside white dwarfs cannot support masses larger than the Chandrasekhar mass [53], any mass accretion leading to the star exceeding this limit will cause it to gravitationally collapse [97]. Although the exact processes leading to supernova explosions are still under debate, SNe Ia have similar light profiles and peak magnitudes, and can thus be used as standardizable candles. By measuring the distance moduli µ for SNe Ia at redshifts zSNe we can thus infer the luminosity distance dL(z), i.e. µ ¶ dL(zSNe) µ(zSNe) mSNe M 5log , (1.69) = − = 10 10[pc] where mSNe is the observed SNe Ia magnitude and M is the absolute magnitude. This relation follows directly from Eq. 1.12 upon taking the logarithm. In recent years the absolute mag- nitudes of Type Ia Supernovae have been found to exhibit correlations with their color, the
21 Chapter 1. Introduction duration of the explosion and host galaxy properties, which need to be taken into account in cosmological analyses.
Hubble constant measurements
The Hubble constant H0 measures the local expansion rate in the Universe and it is the only cosmological parameter that can be measured both assuming a cosmological model and without any assumptions on cosmology. The idea behind these measurements is that the intercept of the logarithmic magnitude-redshift relation of SNe Ia is a measure for H0 if the distance modulus of the Supernovae is known [137]. Since SNe with recession velocities dominated by the Hubble flow typically lie at distances of a few tens of Mpc, their distances cannot be measured directly. Instead, direct Hubble constant measurements rely on distance ladder measurements. The latest measurements have employed a distance ladder build on local parallax measurements, the maser-based distance to the NGC 4258 galaxy and Cepheids to calibrate SNe Ia distances and measure H0 [228]. These analyses result in measurements complementary to the CMB and are derived without the assumption of a specific cosmological model.
22 2 Three-dimensional spherical analyses of cosmological spectroscopic surveys
We are all in the gutter, but some of us are looking at the stars. — Oscar Wilde, Lady Windermere’s Fan
This chapter appeared in a similar form as a part of Nicola, Réfrégier, Amara & Paranjape, 2014 [190].
2.1 Introduction
Constraining the nature and properties of dark energy and dark matter are amongst the most intriguing tasks of current physics. Spectroscopic galaxy redshift surveys offer a way to probe the matter distribution at low redshift, which is strongly affected by the properties of the dark sector. Upcoming spectroscopic clustering surveys like DESI [165], HETDEX [123] and PFS [266], are therefore amongst the most promising tools to achieve these tasks. As discussed in Chapter1, galaxy surveys are inherently three-dimensional, making their analysis more complex as opposed to the CMB, which can be analyzed through two-dimensional maps on the sky. Depending on the galaxy survey geometry, different analysis methods have thus been proposed. For surveys with limited angular sky coverage, the sky can be approximated as flat. Therefore, the clustering of galaxies can be analyzed in three-dimensional Cartesian coordi- nates by means of the spatial correlation function ξ(r ) or through its Fourier counterpart, the Cartesian power spectrum P(k,r ) (for details on these quantities, see Chapter1 and references therein).
In recent years, galaxy redshift surveys have become both wider and deeper, leading us to in- vestigate analysis methods other than P(k,r ), which do not rely on the flat-sky approximation and which facilitate combination of galaxy clustering data with other cosmological probes. A statistic which naturally incorporates the curvature of the sky is the spherical harmonic tomog-
23 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys
i j raphy (SHT) power spectrum C` , the spherical harmonic transform of the angular correlation function at redshifts zi (for theoretical studies see e.g. [38, 52, 70, 21] and for application to data see e.g. [111, 129]). The three-dimensional information can partly be retrieved from this tomographic analysis by performing the spherical harmonics decomposition at a number of different redshifts. Tomographic analyses of the matter overdensity field require subdivision of data into bins, since a finite redshift resolution is needed to compute angular correlations in practice.
Another common way to analyze the three-dimensional matter overdensity field in spherical geometry, which has been applied to galaxy redshift surveys (e.g. [115, 219]), weak lensing (e.g. [112,48]) and the integrated Sachs-Wolfe effect [244], is to measure its spherical Fourier transform. The result is the three-dimensional spherical Fourier Bessel (SFB) power spec- trum C`(k,k0) where the angular dependence is encoded in the multipole ` and the radial dependence in the wave vector k. This statistic allows us to retrieve the clustering information without having to adopt the flat sky approximation or the need for redshift binning.
Recently, the spherical harmonic tomography power spectrum has been compared to P(k,r ), showing that the two methods yield consistent results [21] and both these methods have been employed to investigate the complementarity of weak lensing and galaxy redshift surveys (see e.g. [156, 67]).
With the aforementioned galaxy redshift surveys under development, it becomes increas- ingly important to further test and compare the applicability of these statistics to survey requirements. In this chapter, we compare the two spherical-sky statistics, i.e. the SHT power spectrum and the SFB power spectrum using a Fisher analysis. We study the sensitivity of these statistics to the detail of their implementation, placing particular emphasis on the advantages and disadvantages of each method, some of which we illustrate with simplified toy models.
This chapter is organized as follows. In Section 2.2 we summarize our comparison baseline model. In Sections 2.3 and 2.4 we review three-dimensional spherical analyses of the matter overdensity field as well as Fisher matrix forecasting techniques and present applications to the SFB power spectrum. In Section 2.5 we present a comparison of the spherical harmonic tomography and the SFB power spectrum. We conclude in Section 2.6. Derivations and discussion of employed toy models are deferred to the Appendix.
2.2 Comparison baseline
In this work, we consider a wCDM cosmological model in the framework of general relativity specified by the set of 7 cosmological parameters θ (h,Ω ,Ω ,w ,w ,n ,σ ), where we fix = m Λ 0 a s 8 the baryon density Ω 0.045. This model allows for a dynamical evolution of dark energy as b = well as curvature and is characterized by 7 parameters: the mean fractional matter density
Ωm, the fractional density of dark energy ΩΛ, the Hubble constant H0 100h km/s/Mpc, the 1 = r.m.s. of matter fluctuations σ8 in spheres of comoving radius 8h− Mpc, the scalar spectral
24 2.2. Comparison baseline
106 4.0 ×
3.2
2.4 ) z ( W 1.6
0.8
0.0 0.5 1.0 1.5 2.0 2.5 3.0 z
Figure 2.1: Survey window function with mean redshift z¯ 1.2 and constant galaxy vol- 3 ¡ 1 ¢ 3 = ume density n0 10− h− Mpc − . This corresponds to a galaxy surface density of nA 2 = = 0.49 arcmin− .
index ns and two parameters w0 and wa that characterize the equation of state of dark energy ([173, 56])
w(a) w (1 a)w . (2.1) = 0 + − a We choose fiducial values θ (0.7,0.3,0.69, 0.95,0,1.0,0.8), which are consistent with the fid = − recent results by WMAP 9 [124]. In all calculations we fix w 0. a = We define the baseline survey in terms of total surveyed volume V and number of detected galaxies Ngal by expressing the galaxy number density as [282]
Ngal n (r ) n0 φ(r ) φ(r ), (2.2) = = V where r is the comoving distance and φ(r ) defines the radial survey selection function with normalization R d3r φ(r ) V as in [219, 282]. The normalization condition defines the window = function in redshift W (z)[282]. For incomplete sky-coverage, parametrized by the fraction of sky covered in the survey fsky, it becomes
Z Z c Z V dr r 2φ(r ) dz r 2φ(r ) dz W (z) . (2.3) = H(z) = = 4πfsky
For all surveys considered, we assume fractional sky coverage f 0.125, which represents a sky = lower limit to the sky fraction covered by future surveys, and a radial selection function given
25 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys by
2 ³ r ´ φ(r ) e− r0 , (2.4) = 3 3 which yields a survey volume of V π 2 r . This choice is motivated by the fact that it allows for = 0 the analytical computation of SFB power spectra in the absence of evolution, which is useful to test the full results. In order to match window functions characteristic of upcoming galaxy 1 redshift surveys, we set r r (z 1) 2354 h− Mpc, as computed in our fiducial cosmological 0 = = = model.
Fig. 2.1 shows the angular galaxy density as a function of redshift for a volume density of 3 ¡ 1 ¢ 3 n 10− h− Mpc − ; it further illustrates the binning scheme adopted in the tomographic 0 = analysis. Each bin configuration is chosen by requiring the same number of galaxies in each redshift bin. Our baseline choices are summarized in Table 2.1. For comparison and future reference we show recent choices made in the literature in Table 2.2. This table highlights the breadth of possible choices to make when analyzing a galaxy redshift survey. The most important ones include:
1. Statistic used: spherical harmonic tomography power spectrum, spherical Fourier Bessel power spectrum or Cartesian power spectrum
2. Physical effects included: redshift space distortions, relativistic corrections (see e.g. [38, 282])
3. Implementation scheme: examples include simplifying assumptions like the Limber ap- proximation, Fisher matrix computation method or number of cosmological parameters considered
From Table 2.2 we see that the choices made differ considerably. One of the aims of this chapter is therefore to investigate how much parameter forecasts and constraints are influenced by some of these choices.
2.3 Three-dimensional spherical power spectra
2.3.1 The Cartesian power spectrum
A common approach to analyze the observed matter overdensity field δ(x,r ) is to expand it into its Cartesian Fourier components δ(k,r ), where we use the comoving distance r as a measure of time t. The real-space overdensity field is related to its Fourier counterpart through
Z 1 3 ik x δ(x) d k δ(k,r )e · . (2.5) = (2π)3
26 2.3. Three-dimensional spherical power spectra
Table 2.1: Baseline specification
Sky coverage: fsky 0.125 = 1 Selection function: φ(r ) as in Eq. 2.4 with r0 2354 h− Mpc 3 ¡ 1= ¢ 3 Survey Galaxy volume density: n0 10− h− Mpc − = 2 Galaxy surface density: n 0.49 arcmin− A = Angular scales covered: ` [2,50] ∈ Cosmological parameters: θ (h,Ω ,Ω ,w ,n ,σ )a = m Λ 0 s 8 Fiducial values: θ (0.7,0.3,0.69, 0.95,1.0,0.8) fid = − Prior: none Model Galaxy bias: b (k,r ) 1 = Redshift space distortions: yes Relativistic corrections: no SFB SHT Implementation k range: 1 Number of bins: 7 k (0.0007,0.2) hMpc− ∈ Redshift range: 0.05 z 3.0 Fisher matrix: diagonal ≤ ≤ a We further fix Ω 0.045 and wa 0. b = =
The Fourier space correlation function is the Cartesian matter power spectrum P (k,r ) defined by
3 δ(k,r )δ(k0,r ) (2π) δ (k k0)P (k,r ). (2.6) 〈 〉 = D − We only focus on the linear matter power spectrum which factorizes into a time and scale dependent part as
P (k,r ) D2(r )P (k), (2.7) = 0 where D(r ) is the linear growth factor and P (k) is the power spectrum at redshift z 0. We 0 = assume a transfer function as summarized in [256] and neglect both the effects of baryon oscillations (BAOs) and neutrinos in our baseline configurations. In Section 2.5 we investigate the impact of BAOs on our results.
The measurement of this statistic from galaxy redshift surveys bears one complication: the observables in these surveys are the galaxy redshifts z and their angular positions on the sky ¡θ,φ¢. Therefore, in order to compute any three-dimensional power spectrum from data, the redshift needs to be related to a wave number k through the assumption of a radial distance. This transformation depends on the choice of a fiducial cosmological model. Any three- dimensional analysis of the matter overdensity field therefore requires the assumption of a cosmological model [223], prior to testing it.
Another consequence of the fact that radial galaxy distances are only accessible through their redshift, is that the distance estimates will be affected by peculiar galaxy velocities. The
27 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys
Table 2.2: Compilation of different implementations used in the literature.
Rel. Limber a Paper Statistic RSD Bias Np Prior corr. approx. Gaztanaga et i j C ,P(k) b 8 Planck + al., 2012 [90] ` SN-II de Putter et al., i j C ,P(k) b 9 Planck Fisher 2013 [67] ` analysis Cai et al., 2012 i j C ,P(k) b 6 CMB p. s.c [44] ` Kirk et al., 2015 i j C 7 flat [156] ` Font-Ribera et i j C ,P(k) b 8 Planck al., 2013 [84] ` Di Dio et al., i j C 5 none 2014 [70] ` This work i j C` ,C`(k,k0) 6 none Bonvin et al., i j C -- 2011 [38] ` General Challinor et al., i j C -- analysis 2011 [52] ` Rassat et al., C (k,k ) -- 2012 [219] ` 0 Yoo et al., 2013 C (k,k ) -- [282] ` 0 Pratten et al., C (k,k ) -- 2013 [218] ` 0
a Number of parameters used in analysis. i j b Transverse modes with C . RSDs taken into account for the radial modes in each redshift bin using P(k). ` c Assume primordial CMB power spectrum known.
28 2.3. Three-dimensional spherical power spectra comoving galaxy distances s inferred from their redshifts are given by [146] v n s r · , (2.8) = + a H(a) where v is the galaxy velocity due to the linear collapse of overdensities, H(a) is the Hubble parameter and n denotes the line of sight direction. These redshift space distortions (RSDs) lead to an enhancement of the Cartesian matter power spectrum given by [146]
2 P (k,r ) P (k,r )¡1 βµ2 ¢ , (2.9) s = + k
f d logD(a) where β /b. The quantity f /d loga denotes the linear growth rate, b is the galaxy bias, = = discussed below, and µk is the cosine of the angle between the line of sight and the wave vector k. Measuring the power spectrum in redshift space therefore allows us to also estimate the growth rate f of matter perturbations.
Irrespective of the analysis method, galaxy redshift surveys pose an additional complication. Since galaxies are only expected to form inside the peaks of the overdensity field [24] and galaxy formation is not completely understood yet, the galaxy overdensity field δg(k,r ) is expected to constitute a biased tracer of the underlying dark matter distribution δ(k,r ), i.e. δ (k,r ) b(k,r ) δ(k,r ). In this chapter we assume that galaxies perfectly trace dark matter, g = which amounts to setting the bias parameter b(k,r ) 1. Since in this work we focus on = clustering on large scales, where the scale-dependence of galaxy bias is negligible (e.g. [182]), we believe that this simplified assumption is appropriate because the statistics we compare will all equally suffer from the problem of bias. An investigation of the effects of scale-dependent bias on our results would be interesting for future work.
2.3.2 The spherical harmonic tomography power spectrum
The need for assuming a cosmological model, which arises in three-dimensional analyses of galaxy clustering, can be circumvented with a tomographic analysis. This amounts to discretizing the redshift and analyzing the angular dependence of galaxy clustering through i j the spherical harmonic tomography power spectrum C` at a number of different redshifts in order to partly recover the three-dimensional information. In practice, a galaxy catalogue is analyzed by subdividing the galaxies into redshift bins and computing both the auto- and cross-power spectra for all the bins.
Assuming the overdensity field to be statistically isotropic and homogeneous, the spherical harmonic tomography power spectrum including RSDs between redshift bins i and j, with radial selection functions φi (r ) and φj (r ) respectively, is given by [198] Z i j 2 2 ³ i i,r ´³ j j,r ´ C dkk P0 (k) W (k) βW (k) W (k) βW (k) , (2.10) ` = π ` + ` ` + ` where the auto power spectra are obtained for i j and the cross power spectra for i j . The = 6=
29 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys
5 10−
5 10− z¯ = 0.30 6 10−
z¯ = 0.66
7 10− z¯ = 0.88
6
10− | l z¯ = 1.10 l 10 8 C
C − |
z¯ = 1.40 9 10− z¯ = 1.80 7 10− 10 10− z¯ = 2.50
11 10− 101 102 101 102 l l
i j Figure 2.2: Spherical harmonic tomography power spectra C` for the 7 redshift bins and survey specified in Sec. 2.2, where z¯ denotes the mean redshift of each bin. The left panel shows the auto-power spectra, while the absolute value of the cross-power spectra are shown on the right hand side.
R selection functions are normalized i.e. dr φi (r ) 1. Their unnormalized counterparts, the = i redshift distributions for each bin, are shown in Fig. 2.1. W` is the real-space window function i,r whereas W` accounts for the corrections due to RSDs; both window functions are defined in terms of the spherical Bessel functions j` as [198] Z W (k) dr D(r )φ (r )j (kr ), (2.11) ` = i `
Z " ¡ 2 ¢ r 2` 2` 1 W (k) dr D(r )φi (r ) + − j` (kr ) ` = (2` 3)(2` 1) + − `(` 1) (` 1)(` 2) ¸ − j` 2 (kr ) + + j` 2 (kr ) . (2.12) − (2` 1)(2` 1) − − (2` 1)(2` 3) + − + + + i j Fig. 2.2 shows both the auto and the cross (neighboring redshift bins) SHT power spectra C` as a function of angular scale ` for the baseline configuration defined in Sec. 2.2.
The computation of the spherical harmonic tomography power spectrum through Eq. 2.10 can be computationally expensive and it is therefore common to resort to the small angle and wide selection function approximation. At large `, the spherical harmonic tomography power spectrum can be approximated through Limber’s approximation as [168]
Z Ã 1 ! i j φi (r )φj (r ) ` 2 C dr D(r )P0 k + . (2.13) ` ' r 2 = r
30 2.3. Three-dimensional spherical power spectra
2.3.3 The spherical Fourier Bessel power spectrum
An alternative way for analyzing the galaxy overdensity field is through the spherical Fourier Bessel transform (e.g. [115, 219, 218, 282]). The galaxy overdensity field δ(r) can be expanded into its translationally invariant parts i.e. the eigenfunctions of the Laplacian in spherical coordinates. In flat-space these are given by products of spherical Bessel functions and spherical harmonics Y`m(θ,φ) which leads to the expansion r 2 Z X δ(r) dk δ`m(k)k j` (kr )Y`m(θ,φ). (2.14) = π `,m
The coefficients δ`m(k) are given by
r Z Z 1 2 ` 2 3 ¡ ¢ ¡ ¢ δ`m(k) ki dr r d k0 δ k0,r j` (kr ) j` k0r Y ∗ (θk ,φk ). (2.15) = 2π2 π `m 0 0
The SFB power spectrum is defined as the variance of these coefficients as given below, where the last equality holds if the overdensity field δ(r) is statistically isotropic and homogeneous (SIH) [219]
δ`m(k)δ∗ (k0) C`(k,k0)δ`` δmm C`(k)δD(k k0)δ`` δmm . (2.16) 〈 `0m0 〉 = 0 0 = − 0 0 denotes an ensemble average. Under the SIH assumption we further obtain [48] 〈〉 C (k) P (k). (2.17) ` = A similar identity can be obtained in the presence of RSDs in the flat-sky limit and for high radial wavevectors k. In this case we approximately obtain
2 C (k) P (k)¡1 β¢ . (2.18) ` ' + A derivation of this radialization in the presence of RSDs can be found in Appendix 2.A. In cosmology, the SIH condition will usually be violated for two reasons [219] (i) the observed fields are generally confined to finite regions of space defined by the survey selection function and (ii) the field δ(r) and the power spectrum P (k,r ) evolve.
In order to compute the observed SFB power spectrum, constraints on radial survey geometry can be imposed through a radial selection function φ(r ), as defined in Sec. 2.2[219]
δ (r) φ(r )δ(r). (2.19) obs = Accounting for time-evolution as well as the effects of RSDs, the observed SFB power spectrum of the overdensity field δ becomes [219]
µ ¶2 Z 2 2 ¡ r ¢¡ r ¢ C`(k,k0) dk00k00 P0(k00) W`(k,k00) W (k,k00) W`(k0,k00) W (k0,k00) . (2.20) = π + ` + `
31 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys
r W`(k,k00) is the real-space window function whereas W` (k,k00) accounts for the corrections due to RSDs; they are defined as ([219, 218]) Z 2 ¡ ¢ W (k,k00) dr r D(r )φ(r )k j (kr ) j k00r , (2.21) ` = ` `
Z 2 · 2 r 2 k ` ¡ ¢ W` (k,k00) dr r β D(r )φ(r ) 2 j` 1 (kr ) j` 1 k00r = k00 (2` 1) − − + 2 ¸ `(` 1) © ¡ ¢ ¡ ¢ª (` 1) ¡ ¢ + j` 1 (kr ) j` 1 k00r j` 1 (kr ) j` 1 k00r + j` 1 (kr ) j` 1 k00r − (2` 1)2 − + + + − + (2` 1)2 + + + + (2.22) Z µ ¶ 2 k dφ(r ) ` ¡ ¢ (` 1) ¡ ¢ dr r β D(r ) j` (kr ) j` 1 k00r + j` (kr ) j` 1 k00r . + k dr (2` 1) − − (2` 1) + 00 + + Eq. 2.22 allows for a time-dependence of the overdensity field δ(k,r ), since the survey selection functions tend to be broad in redshift as opposed to the redshift bins in Section 2.3.2.
r 1 Fig. 2.3 shows the SFB power spectrum C 0 (k,k) C (k,k)( 0/2p2π)− , both as a function of angu- ` = ` lar scale ` and wave vector k for the selection function defined in Eq. 2.4. The normalization follows [282] and facilitates comparison of the SFB with the Cartesian power spectrum P (k,r ).
Just as for the spherical harmonic tomography power spectrum it is useful to obtain approxi- mations to Eq. 2.20. In [282] it is shown that in the limit of large angular multipoles `, the SFB power spectrum, neglecting RSDs, can be approximated by [282]
à 1 ! à 1 ! 2 ` 2 ` 2 ¡ ¢ C (k,k0) φ + D + P0(k)δD k k0 . (2.23) ` ' k k −
In the cases considered, there always remains a significant difference between Eq. 2.20 and Eq. 2.23, which is the SFB analogue of Limber’s approximation. Nonetheless it proves useful to test the accuracy of the full equation.
The SFB power spectrum for a generic distance-redshift relation
The spherical Fourier Bessel coefficients are functions of the wave vector k, which in turn depends on the measure of separation in real space. It is customary to choose k conjugate to the comoving separation r , making the assumption of a cosmological model inevitable when computing the SFB transform. As an alternative, in this chapter we additionally compute the SFB power spectrum for two distance-redshift relations which can directly be computed from observable quantities.
In order to derive an expression for the SFB power spectrum, we assume a generic distance- redshift relation defined as r˜(z) where r˜ is an arbitrary monotonic function of the redshift z.
32 2.3. Three-dimensional spherical power spectra
Cl0(k = k = 0.0048)
Cl0(k = k = 0.0078)
Cl0(k = k = 0.1700)
] 4
] 4
3 10
3 10 Mpc) Mpc) 1 1 − − h h ) [( ) [( k, k k, k
( 3 ( 0 3 l 0 10 l 10 C C
Cl0(k, k, l = 2)
Cl0(k, k, l = 10)
Cl0(k, k, l = 50) 2 2 10 3 2 1 10 10− 10− 10− 101 1 k [hMpc− ] l
Figure 2.3: The SFB auto power spectrum C`0 (k,k) as a function of wave vector k and angular scale `.
With this choice of distance-redshift relation the SFB power spectrum reduces to
µ ¶2 Z 2 2 ¡ r ¢¡ r ¢ C`(ν,ν0) dk00k00 P0(k00) W`(ν,k00) W (ν,k00) W`(ν0,k00) W (ν0,k00) , (2.24) = π + ` + ` where ν denotes the wave vector conjugate to the new separation measure and φ0(s) is the selection function in this coordinate system. W`(ν,k00) is the real-space window function r whereas W` (ν,k00) accounts for the corrections due to RSDs; they are given by Z 2 ¡ ¢ W (ν,k00) dr r D(r )φ0(r˜)νj (νr˜) j k00r , (2.25) ` = ` `
Z ν2 dr˜ H(z) · `2 r 2 ˜ ˜ ¡ ¢ W` (ν,k00) dr r β D(r )φ0(r ) 2 j` 1 (νr ) j` 1 k00r = k00 dz c (2` 1) − − + 2 ¸ `(` 1) © ¡ ¢ ¡ ¢ª (` 1) ¡ ¢ + j` 1 (νr˜) j` 1 k00r j` 1 (νr˜) j` 1 k00r + j` 1 (νr˜) j` 1 k00r − (2` 1)2 − + + + − + (2` 1)2 + + + + (2.26) Z µ ¶ 2 ν dφ0 dr˜ H(z) ` ¡ ¢ (` 1) ¡ ¢ dr r β D(r ) j` (νr˜) j` 1 k00r + j` (νr˜) j` 1 k00r . + k dr˜ dz c (2` 1) − − (2` 1) + 00 + +
Since the selection function transforms as a scalar quantity, φ0 (r˜) is related to the selection function in terms of the comoving distance through φ0 (r˜(z)) φ(r (z)). For a derivation of = these identities, the reader is referred to Appendix 2.B.
As an illustration of the impact of the choice of distance-redshift relation, we consider two
33 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys alternatives to the comoving distance
c r˜(z) z, (2.27) = H0 c r˜(z) ln(1 z). = H0 +
The first is the linear approximation to the comoving distance valid for low redshifts;, while the second is a logarithmic approximation to r chosen to reproduce both its behavior at low and intermediate redshift. The resulting SFB power spectra are shown in Fig. 2.6. The normalization again follows [282].
2.4 Fisher matrices for 3D spherical power spectra
The Fisher matrix (FM) allows us to forecast the constraints on cosmological parameters obtainable with future surveys under the approximation of Gaussianity (for an overview of Fisher forecasting see e.g. [270, 113] on which this summary is based). This method can be applied to survey optimization or, as done in this chapter, it can be used to assess the constraining power of different data analysis methods. The FM allows for the propagation of uncertainties in the measurement to uncertainties on the model parameters, which here are the parameters of the wCDM cosmological model. Bayes’ theorem allows us to relate the posterior probability distribution p (θ x) around the maximum likelihood (ML) estimator | to the data likelihood L (x;θ). The inverse covariance matrix of the posterior distribution is called the Fisher matrix and given by
∂2 logL Fαβ . (2.28) = 〈− ∂θα∂θβ 〉
When several parameters are simultaneously estimated from the data, the marginalized uncer- q tainty on each parameter θ is bounded by ∆θ F 1 [152]. The fixed uncertainty, obtained α α ≥ αα− when keeping all parameters except one fixed, is smaller or equal to the former and given by 1 ∆θ /pFαα [152]. α ≥
2.4.1 The Fisher matrix for the spherical harmonic tomography power spectrum
The FM for a tomographic survey employing N redshift bins can be derived from Eq. 2.28 assuming a Gaussian likelihood for the spherical harmonics coefficients. The result is [133]
2 1 h i X ( ` )∆` ˜ 1 ˜ 1 Fαβ fsky + Tr D`αC`− D`βC`− , (2.29) = ` 2
34 2.4. Fisher matrices for 3D spherical power spectra where the sum is over bands of width ∆` in the power spectrum and we set ∆` 1. The data = covariance is an N N matrix given by × x x x x [C ]i j C i j N i j , (2.30) ` = ` + ` where the xi ,x j denote the respective bins. The first term in Eq. 2.30 represents the innate xi x j 1 cosmic variance, while the second term is due to shot noise and given by N /nA, where ` = nA is the galaxy surface density of the survey.
The matrix D`α contains the dependence of the observables on the parameters θα and has elements given by
x x ∂C i j i j ` [D`α] . (2.31) = ∂θα
The simple scaling with fsky accounts for the fact that angular modes become coupled for incomplete sky coverage. This reduces the number of independent modes at a given angular scale ` and therefore increases the uncertainties due to cosmic variance [238].
2.4.2 The Fisher matrix for the SFB power spectrum
The computation of the FM for the SFB power spectrum from the Gaussian likelihood for the SFB coefficients δ`m(k) is challenging due to the correlations between different k modes, which are due to time-evolution of the overdensity field and finite survey effects. The compli- cation arising from the non-diagonal data covariance matrix can be dealt with in two different ways: (i) by choosing a finite grid in k space and computing the FM on this discrete grid or (ii) by approximating the full FM by assuming a diagonal data covariance matrix. Drawing from previous work ([268, 116, 112, 282]) we can find expressions for the FM in both cases.
The FM for a measurement of the SFB power spectrum for n discrete wave vectors ki can be written as [112]
2 1 · Cˆ Cˆ ¸ X ( ` )∆` ˆ 1 ∂ ` ˆ 1 ∂ ` Fαβ fsky + Tr C`− C`− , (2.32) = ` 2 ∂θα ∂θβ where the sum is over bands of width ∆` in the power spectrum and we set ∆` 1 and the = scaling with fsky accounts for incomplete sky coverage. Cˆ ` is the non-diagonal covariance matrix for given angular multipole `
C˜ (k ,k ) C˜ (k ,k ) C˜ (k ,k ) ` 1 1 ` 1 2 ··· ` 1 n C˜ (k2,k1) C˜ (k2,k2) C˜ (k2,kn) ˆ ` ` ··· ` C` . . . , (2.33) = . . .. . . . . . C˜ (k ,k ) C˜ (k ,k ) C˜ (k ,k ) ` n 1 ` n 2 ··· ` n n
35 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys and C˜ (k ,k ) C (k ,k ) N (k ,k ). The first term is again the cosmic variance and the ` i j = ` i j + ` i j second is the shot noise in a survey with galaxy volume density n¯(r ) given by [282]
µ ¶Z 2ki k j 2 ¡ ¢ 1 N`(ki ,k j ) dr r φ(r )j` (ki r ) j` k j r . (2.34) = π n¯(r )
If we assume a broad window function, such that mode coupling can be neglected [82], we can approximate C (k,k0) 0 for k k0. This allows us to obtain a simplification of Eq. 2.32 ` = 6= given by
k Zmax X (2` 1)∆` Ldk 1 ∂C`(k,k) ∂C`(k,k) + Fαβ fsky 2 , (2.35) = ` 2 2π (C`(k,k) N`(k,k)) ∂θα ∂θβ kmin + where the sum is over bands of width ∆` in the power spectrum, L denotes the maximal length scale probed in the survey and kmin, kmax denote the wave vector limits of the survey. For our calculations we set L to the characteristic survey depth i.e. L r 1 and ∆` 1. = 0 = For a detailed derivation of Equations 2.32 and 2.35 the reader is referred to Appendix 2.C. We note that we do not include any optimal weighting of the data [115], a subject which will be interesting for future work.
2.5 Results
As a means for assessing the applicability of both the spherical harmonic tomography and the spherical Fourier Bessel power spectrum to upcoming galaxy redshift surveys, we compare their forecasted performance in a Fisher matrix analysis. From the numerous possible combi- nations discussed in Section 2.2, we have chosen to place our emphasis on two topics: We first focus on each statistic separately and address the main complication associated with it; then we compare the constraining power of both statistics for the baseline survey (Section 2.2).
2.5.1 Spherical harmonic tomography power spectrum
The constraints on cosmological parameters obtained when analyzing the baseline survey (2.2) using the SHT power spectrum and only taking into account auto-power spectra 2 are highlighted in Table 2.3. For the baseline configuration we obtain uncertainties of the order of the parameter value, which is due to the fact that we only consider large-scale informa- tion from angular multipoles ` [2,50] 3. Increasing the maximal angular scale probed to ∈
1The maximal length scale probed L is not a well-defined quantity, but parameter constraints seem stable 1 against changing specification, since setting L V 3 changes results by at most 10%. = survey 2We find that including the cross-correlations does not affect results significantly and we thus neglect them in order to match the specifications used for the SFB power spectrum more closely. 3We believe that this reduced range does not affect our results because we are mainly concerned with comparing two different statistics.
36 2.5. Results
Table 2.3: Parameter constraints obtained for different implementations of the SHT and SFB power spectrum. Results for the baseline configuration are marked in bold. Note that in the case of the SFB the full cov. results neglect contributions due to shot noise, whereas all other constraints assume shot noise as specified in Sec. 2.2.
a Statistic Implementation Radial res. σh σΩm σΩΛ σw0 σns σσ8 Limber ` 50 7 bins 5.5 1.8 1.2 2.8 0.79 0.22 max = no ` 50 7 bins 2.2 0.32 0.50 1.2 1.4 0.15 RSD max = SHT RSD ` 200 7 bins 0.082 0.040 0.10 0.23 0.10 0.030 max = ` 50 7 bins 0.42 0.17 0.48 1.0 0.33 0.15 max = ` 50 max = 7 bins 0.38 0.12 0.45 0.97 0.28 0.11 RSD w/ BAOs ` 50 10 bins 0.078 0.098 0.26 0.59 0.30 0.084 max = ` 50 20 bins 0.15 0.073 0.19 0.38 0.20 0.074 max = `max 50 30 bins 0.12 0.058 0.17 0.34 0.18 0.060 = 1 Log. 0.2 h Mpc− 0.30 0.048 0.32 0.46 0.31 0.67 1 full cov. Lin. 0.2 h Mpc− 0.22 0.047 0.29 0.40 0.25 0.49 1 Com. 0.2 h Mpc− 0.38 0.043 0.55 1.3 0.38 0.86 1 Log. 0.2 h Mpc− 0.33 0.059 0.52 0.71 0.30 0.74 1 SFB Linear 0.2 h Mpc− 0.53 0.086 0.82 1.0 0.39 1.2 Com. 0.2 h Mpc 1 0.37 0.091 0.72 2.7 0.32 0.81 diag cov. − Com. w/ 0.2 h Mpc 1 0.11 0.028 0.31 1.2 0.15 0.26 BAOs − 1 Com. 0.15 h Mpc− 0.48 0.097 0.79 2.8 0.51 1.1 1 Com. 0.1 h Mpc− 0.55 0.10 0.92 2.9 0.62 1.2 a For the SFB power spectrum this corresponds to kmax.
37 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys
` 200, which corresponds to the non-linear cutoff for the lowest redshift bin, consider- max = ably improves parameter constraints. The restriction to large-scale angular perturbations is due to the calculation of the SFB power spectrum, which becomes slow for smaller scales. Nonetheless, the matter density of the universe and the power spectrum amplitude are already sensibly constrained whereas the dark energy sector is poorly constrained due to the signif- icant degeneracies present. Before comparing these results to the constraints obtained for the SFB power spectrum, we first discuss the main complication associated with tomographic analyses.
The SHT power spectrum necessitates tomographic analyses of galaxy catalogues, which amounts to splitting the data into redshift bins. This additional freedom raises the question of how to optimally perform this subdivision. For a fixed baseline survey and therefore data, we expect to see small differences between binning schemes. As we will show below, on the contrary, we can identify instabilities when implementing different bin configurations for a given survey, when we estimate their respective constraining power in a Fisher analysis using Eq. 2.29. Not only is this behavior unexpected but it also implies that these instabilities need to be kept in mind when e.g. comparing the forecasted performance of future surveys.
As an example for studying the effects of redshift binning on parameter constraints, we investigate configurations that differ in the amount of redshift bin overlap. In general there are two causes for bin overlap in galaxy redshift surveys: (i) in spectroscopic surveys, bins can purposely be tailored to have overlap while (ii) in photometric surveys, redshift bins will overlap due to inaccurate redshift measurements. We focus on spectroscopic surveys and therefore only consider case (i). A nonzero overlap between redshift bins will cause them to be correlated, if we assume that they are both located in the same part of the sky. Investigating the impact of bin overlap/correlation on the constraining power of galaxy redshift surveys therefore not only highlights instabilities with data binning but it also addresses the core of the same sky-different sky issue (see e.g. [156, 67]), which is the question of how much correlations between data sets can affect parameter constraints.
We investigate the effects of bin configuration on constraining power using a series of highly simplified toy models, which are based on the Limber approximation and ignore shot noise contributions. For a detailed description of these, the reader is referred to [189]. As shown in [189], we find that increasing the amount of overlap between bins, while keeping their mean redshifts fixed and taking into account correlations, can result in an improvement of cosmological parameter constraints by as much as a factor of 2. This behavior is only found when constraining parameters that exhibit a high level of redshift degeneracy between each other i.e. parameters which can only be distinguished with information at separate redshifts; an example from cosmology is the redshift degeneracy between parameters which control the growth of structure as a function of time and the overall clustering amplitude. On the other hand, we find that constraints on non-redshift degenerate parameters as well as fixed errors are insensitive to changes in bin overlap. It is important to point out that these conclusions do not apply to bin overlap caused by redshift uncertainties (case (ii)). Redshift errors cannot be
38 2.5. Results modeled solely as a redshift bin broadening, since this approach does not take into account the uncertainty introduced in the redshift distribution. If redshift uncertainties are implemented as in [17], we find that increased bin overlap due to larger redshift uncertainties worsens parameter constraints as intuitively expected.
These results suggest that the main effect of overlap between redshift bins on spectroscopic surveys is to break redshift degeneracies between parameters. This seems counterintuitive but as we show in Ref. [189], the dependence of parameter constraints on correlation is a generic feature of such data sets. This suggests that the observed sensitivity of parameter constraints on binning scheme is due to the fact that the amount of correlation between redshift bins, which has an effect on parameter constraints, is scheme dependent.
The results presented so far have been based on tomographic analyses consisting of only two redshift bins; as the number of redshift bins is increased, the effect of bin overlap becomes negligible. The more available cosmological information is recovered from a survey, the less sensitive parameter constraints become to binning schemes. In order to obtain stable parameter constraints from a tomographic analysis of galaxy redshift surveys it is therefore essential to ensure that the available information is well recovered by the survey. We will review the limitations imposed on a tomographic analysis returning back to our baseline survey.
1.2 σh
σΩm 1.0 σΩΛ
σw0
0.8 σns
σσ8
σ 0.6
0.4
0.2
7 10 20 30 NBin
Figure 2.4: Uncertainties (1 σ) on cosmological parameters obtained with SHT power spec- trum as a function of the number of redshift bins NBin.
The parameter that mainly controls the constraining power of a tomographic survey is its radial resolution; in practice this is the number of redshift bins. Starting from the baseline survey, we increase the number of redshift bins from N 7 to N 30 as shown in Table Bin = Bin = 2.3 and Fig. 2.4. This reduces uncertainties by almost a factor of three, which is in agreement 1 with the rough /pNBin scaling of parameter constraints with bin number when shot noise is not
39 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys yet dominant [70]. Therefore the maximal cosmological information retrievable analyzing a survey using the SHT power spectrum is limited by redshift accuracy and shot noise.
As seen from Table 2.3 we reanalyze the baseline survey neglecting redshift space distortions; once performing the full calculation and once assuming the Limber approximation. The results indicate that both changes deteriorate parameter constraints, showing that including RSDs in tomographic analyses increases the amount of cosmological information. Furthermore, the results obtained with the Limber approximation and the exact calculation deviate significantly (differences of up to a factor of 5), which suggests that approximations in power spectrum calculations should be used carefully when computing Fisher matrices.
2.5.2 Spherical Fourier Bessel power spectrum
The constraints on cosmological parameters obtained when analyzing the baseline survey (2.2) using the spherical Fourier Bessel power spectrum are emphasized in Table 2.3. In agreement with the results for the SHT power spectrum, we obtain constraints of the same order of mag- nitude as the parameters themselves, an effect which we again attribute to the small multipole range considered. The best-constrained parameter is the matter density Ωm, whereas the SFB analysis mostly yields larger uncertainties on the remaining cosmological parameters than its tomographic counterpart. As in the previous section, we address particular issues associated with this analysis before turning to the comparison between the two methods.
Comparison of Fisher matrix computation techniques
The computation of Fisher matrices for the SFB power spectrum through Eq. 2.32 4 is time- consuming and numerically challenging because it requires inverting the covariance matrix defined in Eq. 2.33. This step is tricky for two reasons: Firstly, as seen from Fig. 2.3, the SFB
` power spectrum falls off sharply for large scales with k /r0 for given angular scale ` and ¿ 1 survey depth r . This is because in a cone of angular extent θ /`, large radial modes cannot 0 ∼ be measured when the survey depth is finite. The rapid decrease in large scale power results in a considerable dynamic range in the covariance matrix eigenvalues, making it almost singular. Secondly, neighboring wave vectors k are strongly correlated, which further complicates the inversion of the covariance matrices.
In order to estimate the SFB Fisher matrix through Eq. 2.32 it is therefore inevitable to restrict
` calculations to separated wave vectors with k /r0 to overcome numerical instabilities. In ≥ practice we cut off all large-scale information for each angular scale ` as soon as it causes the covariance matrix condition number, which is a measure for accuracy loss in matrix inversion, to exceed κ 102. This is possible but it seems desirable to investigate alternatives to this crit = 4There is one subtlety involved with Fisher matrix calculations for the SFB power spectrum: as can be seen from Eq. 2.34, the shot noise contribution to the power spectrum is essentially a galaxy number count and therefore cosmology dependent. Since we ignore information from non-linear scales, we won’t include this information when estimating the constraining power of a survey.
40 2.5. Results
“brute-force” approach.
Even though neighboring wave vectors are strongly correlated, the correlations tend to rapidly decrease as we move away from the diagonal. This suggests resorting to the approximation ¡ ¢ C k,k0 0 for k k0 in order to obtain useful approximations to Eq. 2.32. The most straight- ` ' 6= forward implementation of these ideas is given in Eq. 2.35. Despite being an approximation to the full Fisher matrix given in Eq. 2.32, we find that Eq. 2.35 yields no-shot noise as well as shot noise constraints which are mostly accurate to better than a factor of 2 for the baseline survey (for detailed results, see Table 2.4 in Appendix 2.D). These results agree with those obtained for the Cartesian matter power spectrum [268] and encourage the use of Eq. 2.35 for fast calculations which allow errors of up to a factor of a few.
Making use of this simplifying approximation, we investigate the impact of the non-linearity wave-vector cut on parameter constraints. As can be seen from Table 2.3 and Fig. 2.5 decreasing the maximal wave vector kmax by a factor of two increases parameter uncertainties by almost the same amount. The increase is larger than theoretically expected, since decreasing the cutoff scale by a factor of two will halve the number of available modes and should thus lead to an increase in uncertainty by a factor p2. We believe that this is due to the fact that increasing the cutoff scale additionally helps breaking parameter degeneracies, since we approximately observe the theoretical scaling for the fixed parameter constraints.
The choice of distance-redshift relation
As discussed in Sec. 2.3.3, the need for assuming a fiducial cosmological model can be miti- gated by analyzing surveys using distance-redshift relations that can directly be computed from observable quantities. Using the expressions derived in Sec. 2.3.3 we can investigate the impact of the distance-redshift relation on the obtained power spectra as well as survey constraining power.
We focus on two simple alternatives to the comoving distance as given in Eq. 2.27. Both approximations are fairly accurate at low redshift; at very high redshifts on the other hand, both approximations break down because they considerably overestimate the comoving separation. The SFB power spectra obtained with these two distance-redshift relations are shown in Fig. 2.6. The choice of a different fiducial distance causes a shift in the observed SFB power spectra, because the window functions are offset from those in comoving distance.
These simple distance-redshift relations are viable alternatives to the comoving distance only if analyzing a survey in terms of them does not significantly reduce its constraining power. To test their performance, we compare their constraints on cosmological parameters for a survey as defined in 2.2 in two different ways: Since the diagonal Fisher matrix is an acceptable approximation to the full calculation, we will employ it to compare the forecasted parameter constraints obtained with all three distance measures taking shot noise into account. As a mere illustration, we additionally compare the constraints obtained from the full Fisher matrix,
41 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys
1.6 σh
σΩm
σΩΛ
σns 1.2 σσ8
σ 0.8
0.4
0.0 0.08 0.12 0.16 0.20 1 kmax [hMpc− ]
Figure 2.5: Uncertainties (1 σ) on cosmological parameters obtained using the SFB power spectrum as a function of maximal wave number kmax. Constraints on w0 are not shown for clarity.
neglecting any shot noise contributions.
Both these results are shown in Table 2.3. The constraints are similar, irrespective of the distance-redshift relation chosen. The only parameter exhibiting a significant dependence on the way distance is related to redshift is the dark energy equation of state parameter w0. An analogous behavior is seen in the no shot-noise constraints obtained with the full Fisher matrix: we obtain different constraints especially for those parameters, which the comoving distance depends on, whereas the remaining constraints are largely insensitive to the distance- redshift relation of choice. Choosing a distance-redshift relation other than the comoving distance therefore appears to cause the SFB power spectrum to exhibit a stronger dependence on these former parameters because the volume element and the distance in Eq. 2.25 and 2.26 do not change in the same way.
As we include contributions due to shot noise, this potential constraining power is considerably reduced because alternative distance-redshift relations tend to overestimate the comoving separations for large redshift. The shot noise is therefore enhanced, which largely removes the gain from increased sensitivity.
The above considerations illustrate that constraints obtained from an SFB analysis of a galaxy redshift survey seem to be mostly stable against changes in the assumed distance-redshift relation. This suggests that it could be possible to analyze galaxy clustering using distance- redshift relations which only depend on observable quantities, without too large a loss in constraining power.
42 2.5. Results
2.5.3 Comparison between spherical harmonic tomography and SFB power spec- trum
After focusing on each of the two statistics separately we can compare their baseline con- straints shown in Table 2.3. Unexpectedly, the survey constraining power is significantly affected by the choice of analysis method: When baryon acoustic oscillations are neglected, the SFB power spectrum yields weaker constraints, particularly on those cosmological pa- rameters that need redshift leverage in order to be distinguished (i.e. growth and amplitude parameters). This is already evident for our baseline but the effect is enhanced if we consider larger bin numbers in the SHT analysis.
The weakness of constraints on growth as well as amplitude of matter fluctuations seems to be a generic feature of SFB analyses: When the underlying field is statistically isotropic and homogeneous (SIH), the SFB coefficients are given by Eq. 2.36 i.e. they correspond to an angular average of the Cartesian Fourier coefficient. In practice, the SIH condition is not met due to time evolution of the overdensity field and finite survey effects and the SFB coefficients are related to their Cartesian counterparts through Eq. 2.15. Any violation of the SIH condition will therefore introduce a coupling between the considered radial and angular scales k, ` and the redshift at which δ(k) will mostly be measured. Around r0, the decrease in the selection function amplitude breaks the SIH condition which causes modes to add incoherently and leads to cancellations. The contributions to the SFB power spectrum will therefore be preferentially weighted toward lower redshift, which results in smaller redshift leverage and weaker combined constraints.
These results therefore suggest that, even though both analysis methods are equivalent for infinite survey extent and recovery of all available modes, because the information content of the overdensity field does not depend on the basis set in which it is analyzed, they appear not to be equivalent for finite surveys. Analyses of galaxy redshift surveys through the SFB power spectrum are more affected by finite survey effects, which means that some high redshift information will be downweighted. The SHT power spectrum analysis, on the other hand, allows us to probe more efficiently the complete high redshift range of the survey.
We test the impact of baryon oscillations on our results by reanalyzing our baseline config- urations using the transfer function as specified in [78]. From Table 2.3 we see that adding BAOs improves the SFB constraints, while leaving the SHT power spectrum constraints mostly unchanged. This suggests that a significant fraction of the information lost due to finite survey effects in SFB analyses can be compensated by the fact that its three-dimensional nature allows us to recover the information contained in the BAOs, while an SHT analysis tends to dilute these features 5. The inclusion of BAOs results in comparable constraints for both methods. This suggests that SFB is better suited for capturing information beyond the smooth
5The change of transfer function from [256] to [78] without BAOs (the so-called “no-wiggles” approximation obtained with fractional baryon density Ω 0 but ignoring the oscillatory contribution to the power spectrum) b 6= does not significantly affect parameter constraints, which suggests that the improvement in constraining power can be mainly attributed to the presence of BAOs.
43 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys shape of the power spectrum.
2.5.4 Implementation effects on estimated survey constraining power
The constraining power of a particular survey, as estimated from Fisher matrix calculations, is clearly determined by survey specifications and included physics. Nevertheless, details in the particular implementation can also affect parameter constraints and we can investigate the magnitude of this effect using our simplified toy model (see Ref. [189]). Our findings suggest that the magnitude can potentially equal that of changing survey specifications, when one of the following two conditions is fulfilled: (i) As illustrated in Ref. [189], the choice of binning scheme and thus implementation ceases to be relevant as more information is retrieved from a particular survey and parameter constraints become tighter. This therefore suggests that the choice of prior can significantly affect the stability of parameter constraints: applying a tight prior reduces the susceptibility of parameter constraints to implementation. (ii) A second essential choice is the set of constrained parameters: as highlighted by the toy model, parameters which are redshift degenerate with each other are particularly affected by changes in implementation. On the other hand, constraints on non-degenerate parameters are expected to be more stable. It is thus important to be aware of these instabilities whenever constraining a set of degenerate and loosely constrained parameters.
These susceptibilities can further be amplified due to numerical instabilities in Fisher matrix calculations. Fisher matrices can have large condition numbers κ(F ) i.e. be close to numeri- cally singular, if one or more parameters are not well constrained by the data. Therefore the error introduced by the Fisher matrix inversion can be of order 100%, if the accuracy in the 1 Fisher matrix elements is smaller than κ(F )− [274].
These findings further suggest that care has to be taken when comparing Fisher matrix results. Whenever two different results need to be compared it is essential to make sure that not only the survey specifications are similar but also that priors and set of constrained parameters agree with each other. This is relevant in light of the recent discussion regarding benefits of performing spectroscopic and photometric galaxy surveys in the same part of the sky (see e.g. [156, 67]): a reliable comparison between the results obtained by different groups seems difficult due to the differing choices of priors and constrained parameters. Exactly matching the baseline survey is especially important in this case since it investigates the impact of cross-correlations on parameter constraints, which, as indicated by the toy model, only has an effect when parameters are loosely constrained and degenerate with each other. This suggests that changes in implementation have the potential to even affect qualitative results in this particular case.
44 2.6. Conclusions
2.6 Conclusions
Using a Fisher analysis, we have investigated three-dimensional analyses of galaxy redshift surveys on a spherical sky. In the course of our analysis, it has become evident that Fisher matrix results need to be carefully analyzed and compared. Especially when Fisher matrix methods are used to forecast constraints on large parameter sets, exhibiting degeneracies among one another, the obtained constraints are susceptible to details in implementation.
In particular we have compared the SHT and the SFB power spectrum, two statistics that are designed for the analysis of galaxy redshift surveys in a spherical geometry. By comparing their forecasted constraints on cosmological parameters, we have shown the applicability of approximations, such as the Limber approximation, and the numerical issues associated with computing these statistics. We have also studied the sensitivity of these statistics to the detail of their implementation. Our analysis is based on several simplifying assumptions; in particular, we restrict ourselves to linear, unbiased galaxy clustering and only focus on large scale power spectrum modes. For future work it would be interesting to include a treatment of these effects.
Using toy models, we find that constraints obtained from a tomographic analysis of galaxy redshift surveys can be susceptible to implementation effects like redshift bin overlap, if only a limited amount of the total available information is retrieved. This suggests that in order to be stable against changes in implementation, it is important to retrieve as much information as possible from the tomographic analyses, e.g. by using a large number of redshift bins.
The computation of the SFB power spectrum from data relies on the assumption of a distance- redshift relation, usually the cosmology-dependent comoving distance. Analyses of galaxy redshift surveys by means of the SFB power spectrum therefore require the assumption of a cosmological model, prior to testing it. Comparing the SFB parameter constraints obtained using alternative distance-redshift relations, we find them to be largely stable against changes in the assumed distance. This suggests that future surveys could in principle be analyzed using distance-redshift relations only relying on observable quantities without too large a degradation in parameter constraints.
For the baseline survey configuration we considered, we find that the SHT power spectrum yields somewhat tighter constraints than the SFB power spectrum. When we add baryon oscillations, on the other hand, the two methods yield comparable constraints. We attribute the former to the fact that the SFB power spectrum is less sensitive to modes at high redshift near the survey boundary, while the SHT power spectrum can be tailored to probe these modes. In the presence of BAOs, this effect can be compensated by the fact that the three- dimensional nature of the SFB transform allows us to resolve the baryonic oscillations, which tend to be diluted in SHT analyses. This fact would make SHT analyses advantageous for future spectroscopic galaxy redshift surveys mainly focusing on the power spectrum shape, while the SFB power spectrum may be well suited for specific applications like BAOs.
45 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys
2.A The radialization of the SFB power spectrum in the presence of RSDs
The overdensity field in the absence of RSDs at a constant time r is isotropic and homogeneous and can be directly related to the Cartesian matter power spectrum as in Eq. 2.17[ 48]. The isotropy is broken in the presence of RSDs but an approximate relation between these two quantities still holds. The SFB coefficients of the overdensity field are related to their Cartesian Fourier counterpart through Z 1 ` δ`m(k,r ) ki dΩk δ(k,r )Y`∗m(θk ,φk ). (2.36) = p8π3
The contribution to the Cartesian Fourier coefficient due to RSDs is given by
δ (k,r ) βµ2 δ(k,r ), (2.37) RSD = k ¡ ¢ 2 4 with a power spectrum δ (k,r )δ∗ k0,r β µ P(k,r ) δ (k0 k). The quantity µ is 〈 RSD RSD 〉 = k D − k the cosine of the angle between the line of sight direction and the wave vector direction. In the flat-sky limit we can assume that the line of sight direction is constant. For a fixed angular multipole `, the SFB power spectrum will obtain contributions from increasingly radial wave vectors for larger k. In the flat-sky limit and large radial wavevectors k we can thus approximate µ 1. Inserting this into Eq. 2.36 using the identities k ' 1 Z 3 (k0 k) x δD(k0 k) d xe − · , (2.38) − = (2π)3 and [9]
ik r X ` ˆ e · 4π i j`(kr )Y`∗m(k)Y`m (nˆ), (2.39) = `,m gives the contribution to the SFB power spectrum due to RSDs
2 δ`m(k)δ∗ (k0) β P(k,r ) δD(k0 k)δ`` δmm . (2.40) 〈 `0m0 〉 = − 0 0 In the flat-sky and large wave vector limit even the RSD contribution approximately radializes in absence of a selection function and time-dependence of the overdensity field. This behavior is perceivable in Fig. 4 of [282], illustrating that the curvature of the sky is negligible for small scale perturbations.
2.B The SFB power spectrum for a generic distance-redshift rela- tion
The need for assuming a cosmological model before testing it can be avoided choosing a distance-redshift relation which does not depend on cosmology. In order to derive an expres-
46 2.B. The SFB power spectrum for a generic distance-redshift relation sion for the SFB power spectrum, we assume a generic distance-redshift relation defined as r˜(z) where r˜ is an arbitrary monotonic function of the redshift z. The measured redshift will be affected by peculiar galaxy velocities v along the line of sight n [200] i.e. v n z ztrue · , (2.41) obs ' + ac where c is the speed of light. The distance s inferred from the galaxy redshifts therefore becomes
v n dr˜ s r˜(z ) r˜(ztrue) · . (2.42) = obs ' + ac dz
The overdensity field can be decomposed in the SFB basis set with coefficients given by
r Z 2 3 δ`m(ν) d r˜φ0(r˜)δ(r˜)νj` (νr˜)Y ∗ (θ,φ), (2.43) = π `m where ν denotes the wave vector conjugate to r˜ and φ0(r˜) is the selection function in the new coordinate system. Since the overdensity field is independent of the distance measure we have d3r˜ δ(r˜) d3r δ(r) where r is the comoving distance. Eq. 2.43 therefore reduces to = r Z 2 3 δ`m(ν) d rφ0(r˜)δ(r)νj` (νr˜)Y ∗ (θ,φ). (2.44) = π `m
Following [116], the functions of r˜ can be expanded as
µ ¶ dφ0 v n dr˜ φ0(r˜) φ0(r˜) · , (2.45) ' + dr˜ ac dz dj (νr ) µv n dr˜ ¶ j (νr˜) j (νr˜) ` · , ` ' ` + dr˜ ac dz which can be inserted into Eq. 2.44 to yield to first order
r ½Z Z 2 1 3 3 ik0 r δ`m(ν) d r d k’φ0(r˜)δ(k0)e · νj` (νr˜)Y ∗ (θ,φ) = π (2π)3 `m Z Z ¾ 3 3 v(k0) n ik0 r dr˜ d d r d k’ · e · ν [φ0(r˜)j (νr˜)]Y ∗ (θ,φ) . (2.46) + ac dz dr˜ ` `m
The linear continuity equation allows us to relate the Fourier transform of the galaxy velocity field to the overdensity through [71]
aH(z)δ(k) v(k) iβ k. (2.47) = k2
Thus the SFB power spectrum for the distance-redshift relation r˜(z) reduces to
µ ¶2 Z 2 2 r 2 C`(ν,ν0) dk00k00 P0(k00) W`(ν,k00) W (ν,k00) δ`` δmm . (2.48) = π | + ` | 0 0
47 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys
r W`(ν,k00) is the real-space window function whereas W` (ν,k00) accounts for the corrections due to RSDs; they are given by Z 2 ¡ ¢ W (ν,k00) dr r D(r )φ0(r˜)νj (νr˜) j k00r , (2.49) ` = ` `
Z ν2 dr˜ H(z) · `2 r 2 ˜ ˜ ¡ ¢ W` (ν,k00) dr r β D(r )φ0(r ) 2 j` 1 (νr ) j` 1 k00r = k00 dz c (2` 1) − − + 2 ¸ `(` 1) © ¡ ¢ ¡ ¢ª (` 1) ¡ ¢ + j` 1 (νr˜) j` 1 k00r j` 1 (νr˜) j` 1 k00r + j` 1 (νr˜) j` 1 k00r − (2` 1)2 − + + + − + (2` 1)2 + + + + (2.50) Z µ ¶ 2 ν dφ0 dr˜ H(z) ` ¡ ¢ (` 1) ¡ ¢ dr r β D(r ) j` (νr˜) j` 1 k00r + j` (νr˜) j` 1 k00r . + k dr˜ dz c (2` 1) − − (2` 1) + 00 + +
Since the selection function transforms as a scalar, φ0 (r˜) is related to the selection function in comoving distance through φ0 (r˜(z)) φ(r (z)). = Any measurement of the SFB power spectrum will be affected by shot noise. The number of galaxies Ngal observed in a given survey is independent of the distance-redshift relation of choice i.e. Z Z 3 3 N d r˜ φ0 (r˜)n0 (r˜) d r φ(r )n (r ), (2.51) gal = = where n (r ) is the galaxy volume density in comoving coordinates and n0 (r˜) is the volume density in terms of r˜. Together with the identity φ0 (r˜(z)) φ(r (z)), this implies = ¯ 3 ¯ ¯d r¯ n0 (r˜) ¯ ¯n (r ). (2.52) = ¯d3r˜ ¯
For a generic distance-redshift relation the shot-noise therefore reduces to
Z Z ¯ ¯ ¡ ¢ 2νν0 2 1 ¡ ¢ 2νν0 2 ¯dr˜ ¯ 1 ¡ ¢ N` ν,ν0 dr r φ0 (r˜) j` (νr˜) j` ν0r˜ dr r˜ ¯ ¯φ(r ) j` (νr˜) j` ν0r˜ . = π n0 (r˜) = π ¯dr ¯ n (r ) (2.53)
2.C Derivation of the Fisher matrix for the SFB power spectrum
2.C.1 Full Fisher matrix
The FM obtained from a data likelihood with covariance matrix C and mean µ is given by [270]
· ¸ T 1 1 ∂C 1 ∂C ∂µ 1 ∂µ Fαβ Tr C− C− C− . (2.54) = 2 ∂θα ∂θβ + 〈 ∂θα 〉 〈∂θβ 〉
48 2.C. Derivation of the Fisher matrix for the SFB power spectrum
] 4 ] 4
3 10 3 10 Mpc) Mpc) 1 1 − − h h ) [( ) [( k, k k, k
( 3 ( 3 0 0 l 10 l 10 C C
Cl0(k, k, l = 2) Cl0(k, k, l = 2)
Cl0(k, k, l = 10) Cl0(k, k, l = 10)
Cl0(k, k, l = 50) Cl0(k, k, l = 50) 2 2 10 3 2 1 10 3 2 1 10− 10− 10− 10− 10− 10− 1 1 k [hMpc− ] k [hMpc− ]
(a) Logarithmic distance-redshift relation (b) Linear distance-redshift relation
Figure 2.6: The SFB auto power spectrum C`0 (k,k) for two different distance redshift relation as a function of angular scale `.
Assuming a measurement of the SFB coefficients of the matter overdensity field δ`m(k) for a discrete set of radial wave vectors k denoted (k ,k , ,k ), the data covariance matrix is 1 2 ··· n given by
Cˆ 0 0 0 0 `1 ··· ...... ...... . . . . ˆ 0 C`1 0 0 Cˆ ··· ··· , (2.55) = 0 0 0 Cˆ 0 `2 ··· ...... ...... . . . . . 0 0 0 0 Cˆ ··· `max where the sub-covariance matrices are defined as
C˜ (k ,k ) C˜ (k ,k ) C˜ (k ,k ) `i 1 1 `i 1 2 ··· `i 1 n C˜ (k2,k1) C˜ (k2,k2) C˜ (k2,kn) ˆ `i `i ··· `i C`i . . . . (2.56) = . . .. . . . . . C˜ (k ,k ) C˜ (k ,k ) C˜ (k ,k ) `i n 1 `i n 2 ··· `i n n Since Cˆ is blockdiagonal and µ δ (k ) 0, Eq. 2.54 yields = 〈 `m i 〉 = 2 1 · Cˆ Cˆ ¸ X ( ` )∆` ˆ 1 ∂ ` ˆ 1 ∂ ` Fαβ fsky + Tr C`− C`− , (2.57) = ` 2 ∂θα ∂θβ where the sum is over bands of width ∆` in the power spectrum.
49 Chapter 2. Three-dimensional spherical analyses of cosmological spectroscopic surveys
2.C.2 Diagonal Fisher matrix
The computation of the SFB FM through Eq. 2.32 is numerically challenging and it is therefore desirable to investigate possible approximations. Although the amount of cross-correlation
C`(k,k0) between neighboring k vectors can be considerable, it tends rapidly to zero for separated wave vectors. If we assume a broad window function, such that mode coupling can 6 be neglected [82], we can approximate C (k,k0) 0 for k k0 . Assuming a measurement of a ` = 6= set of discrete SFB modes δ (k ) up to ` ` which satisfy δ (k ) 0 and defining `m i ≤ max 〈 `m i 〉 =
δ (k ) δ , (2.58) `m i = `m,i δ (k )δ (k ) C (k ,k ) N (k ,k ) ∆2 , 〈 `m i `m i 〉 = ` i i + ` i i = `,i the data likelihood can be written as
δ2 1 P `m,i ,m,i 1 2 ` ∆2 − `,i L (x;θ) n e . (2.59) = (`max 1)(`max 3) 2 Q 2` 1 (2π) − + `,i ∆`,i+
Applying Eq. 2.28, the FM becomes
X (2` 1)∆` 1 ∂C`(ki ,ki ) ∂C`(ki ,ki ) Fαβ fsky + , (2.60) = 2 (C (k ,k ) N (k ,k ))2 ∂θ ∂θ `,i ` i i + ` i i α β where the sum is over bands of width ∆` in the power spectrum and wave vectors ki . To proceed, we assume that the maximal length scale probed by the survey is given by L. Therefore 2π the minimal measurable mode is k /L which also defines the k-space resolution. The min = 2π maximal measurable mode is determined by the smallest distance ∆L and given by k /∆L. max = Turning the Riemann sum in Eq 2.60 into a continuous integral yields
k Zmax X (2` 1)∆` Ldk 1 ∂C`(k,k) ∂C`(k,k) + Fαβ fsky 2 . (2.61) = ` 2 2π (C`(k,k) N`(k,k)) ∂θα ∂θβ kmin +
2.D Comparison of FM computation techniques for the SFB power spectrum: Results
The parameter constraints obtained for the SFB power spectrum using both Fisher matrix (FM) computation methods and applying the same wave vector cuts, both neglecting and including shot noise contributions, are shown in Table 2.4. Apart from the constraints on w0 the results agree reasonably well, as mentioned in 2.5.2.
6This approximation is equally justified when we assume that the SFB power spectrum is computed for radial wave vector bins which are broader than the correlation scale due to finite survey effects.
50 2.D. Comparison of FM computation techniques for the SFB power spectrum: Results
Table 2.4: Comparison of parameter constraints for the SFB power spectrum obtained using discrete (Eq. 2.32) and continuous (Eq. 2.35) Fisher matrix
Statistic Implementation Shot noise σh σΩm σΩΛ σw0 σns σσ8 Comoving no 0.38 0.043 0.55 1.3 0.38 0.86 full cov. Comoving yes 0.48 0.056 0.68 1.5 0.46 1.1 SFB Comoving no 0.30 0.071 0.57 2.2 0.27 0.66 diag cov. Comoving yes 0.39 0.10 0.77 3.0 0.33 0.85
51
3 Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
Zwei Dinge erfüllen das Gemüt mit immer neuer und zunehmender Bewunderung und Ehrfurcht, je öfter und anhaltender sich das Nachdenken damit beschäftigt: Der bestirnte Himmel über mir, und das moralische Gesetz in mir. — Immanuel Kant, Kritik der praktischen Vernunft, Beschluss
This chapter appeared in a similar form as Nicola, Réfrégier & Amara, 2016 [191].
3.1 Introduction
The past two decades have seen immense progress in observational cosmology that has lead to the establishment of the ΛCDM model for cosmology (for details, see Chapter1 and references therein). This development is mainly based on the combination of different cosmological probes such as the CMB temperature anisotropies, galaxy clustering, weak gravitational lensing, supernovae and galaxy clusters. Until now, these probes have been, for the most part, measured and analyzed separately using different techniques and combined at late stages of the analysis, i.e. when deriving constraints on cosmological parameters. However, this approach is not ideal for current and future surveys such as the Dark Energy Survey (DES1), the Dark Energy Spectroscopic Instrument (DESI2), the Large Synoptic Survey Telescope (LSST3), Euclid4 and the Wide Field Infrared Survey Telescope (WFIRST5) for several reasons. First,
1http://www.darkenergysurvey.org. 2http://desi.lbl.gov. 3http://www.lsst.org. 4http://sci.esa.int/euclid/. 5http://wfirst.gsfc.nasa.gov.
53 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing these surveys will cover large, overlapping regions of the observable universe and are therefore not statistically independent. In addition, the analysis of these surveys requires tight control of systematic effects, which might be identified by a direct cross-correlation of the probes statistics. Moreover, each probe provides a measurement of the cosmic structures through a different physical field, such as density, velocity, gravitational potentials, and temperature. A promising way to test for new physics, such as modified gravity, is to look directly for deviations from the expected relationships of the statistics of the different fields. The integrated treatment of the probes from the early stages of the analysis will thus provide the cross-checks and the redundancy needed not only to achieve high-precision but also to challenge the different sectors of the cosmological model.
Several earlier studies have considered joint analyses of various cosmological probes. Mandel- baum et al. [181], Cacciato et al. [43] and Kwan et al. [162] for example derived cosmological constraints from a joint analysis of galaxy-galaxy lensing and galaxy clustering while Liu et al. [175] used the cross-correlation between the galaxy shear field and the overdensity field to- gether with the cross-correlation of the galaxy overdensity with CMB lensing to constrain multiplicative bias in the weak lensing shear measurement in CFHTLenS. Recently, Singh et al. [248] performed a joint analysis of CMB lensing as well as galaxy clustering and weak lensing. Furthermore, Eifler et al. [76] and Krause et al. [159] have theoretically investigated joint analyses for photometric galaxy surveys by modeling the full non-Gaussian covariance matrix between cosmic shear, galaxy-galaxy lensing, galaxy clustering, photometric baryon acoustic oscillations (BAO), galaxy cluster number counts and galaxy cluster weak lensing.
Extending beyond this, we present and implement an integrated approach to probe combi- nation. In this first implementation we combine data from CMB temperature anisotropies, galaxy overdensities and weak lensing. We use data from Planck 2015 [209] for the CMB, for galaxy clustering we use photometric data from the 8th data release of the Sloan Digital Sky Survey (SDSS DR 8) [13] and the weak lensing shear data comes from SDSS Stripe 82 [20]. We combine these probes into a common framework at the map level by creating projected 2-dimensional maps of CMB temperature, galaxy overdensity and the weak lensing shear field. In order to jointly analyze this set of maps we consider the spherical harmonic power spectra of the probes including their cross-correlations. This leads to a spherical harmonic power spectrum matrix that combines CMB temperature anisotropies, galaxy clustering, cosmic shear, galaxy-galaxy lensing and the ISW [232] effect with galaxy and weak lensing shear tracers. We combine this power spectrum matrix together with the full Gaussian covariance matrix and derive constraints on the parameters of the ΛCDM cosmological model, marginalizing over a constant linear galaxy bias and a parameter accounting for possible multiplicative bias in the weak lensing shear measurement. In this first implementation, we use some conservative and simplifying assumptions. For instance we include a limited range of angular scales for the different probes to reduce our sensitivity to systematics, nuisance parameters and nonlinear corrections. With this, we work under the assumption of Gaussian covariance matrices and with a reduced set of nuisance parameters.
54 3.2. Framework
Data
Maps
M1 M2 MM3 3
Power Spectra ij C`
Cov[Cij,Ci0j0 ] P(M|D) Theory ` `0
Ωm h
σ8 ΩΛ
Figure 3.1: Synopsis of the framework for integrated probe combination employed in this work.
This chapter is organized as follows. In Section 3.2 we describe the framework for integrated probe combination employed in this work. The theoretical modeling of the cosmological observables is summarized in Section 3.3. Section 3.4 describes the data analysis for each probe, especially the map-making procedure. The computation of the spherical harmonic auto- and cross-power spectra is discussed in Section 3.5 and the estimation of the covariance matrix is detailed in Section 3.6. In Section 3.7 we present the cosmological constraints derived from the joint analysis and we conclude in Section 3.8. More detailed descriptions of data analysis as well as robustness tests are deferred to the Appendix.
3.2 Framework
The framework for integrated probe combination employed in this work is illustrated in Fig. 3.1. In a first step we collect data for different cosmological probes as taken by either separate surveys or by the same survey. For our first implementation described below we use cosmological data from the CMB temperature anisotropies, the galaxy overdensity field and
55 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing the weak lensing shear field. After data collection, we perform probe specific data analysis which involves data selection and systematics removal. We then homogenize the data format by creating projected 2-dimensional maps for all probes considered. The common data format allows us to combine the cosmological probes into a common framework at the map level. We compute both the spherical harmonic auto- and cross-power spectra of this set of maps i j and combine them into the spherical harmonic power spectrum matrix C` . This matrix captures the cosmological information contained in the two-point statistics of the maps. In a last step we compute the power spectrum covariance matrix and combine it with theoretical predictions to derive constraints on cosmological parameters from a joint fit to the measured spherical harmonic power spectra. The details of the implementation for CMB temperature anisotropies, galaxy overdensities and weak lensing are described below.
3.3 Theoretical predictions
The statistical properties of both galaxy overdensity δg and weak lensing shear γ, as well as their cross-correlation can be measured from their spherical harmonic power spectra. These generally take the form of weighted integrals of the nonlinear matter power spectrum nl Pδδ(k,z) multiplied with spherical Bessel functions j`(kχ(z)). Their computation is time- consuming and we therefore resort to the the Limber approximation [169, 147, 148] to speed up calculations. This is a valid approximation for small angular scales, typically ` O(10), > and broad redshift bins [200]. For simplicity, we further focus on flat cosmological models, i j i.e. Ω 0, for the theoretical predictions. The spherical harmonic power spectrum C at k = ` multipole ` between cosmological probes i, j {δ ,γ} can then be expressed as: ∈ g
Z i ¡ ¢ j ¡ ¢ µ 1 ¶ i j c W χ(z) W χ(z) ` /2 C dz P nl k + ,z , (3.1) ` = H(z) χ2(z) δδ = χ(z) where c is the speed of light, χ(z) the comoving distance, H(z) the Hubble parameter and i 0 ¡ ¢ W χ(z) denotes the window function for probe i 0.
For galaxy clustering the window function is given by
H(z) W δg ¡χ(z)¢ b(z)n(z), (3.2) = c where b(z) denotes a linear galaxy bias and n(z) is the normalized redshift selection function of the survey i.e. R dz n(z) 1. We focus on scale-independent galaxy bias since we restrict = the analysis to large scales, which are well-described by linear theory.
The window function for weak lensing shear is
2 Z χh 3 ΩmH χ(z) χ(z0) χ(z) γ ¡ ¢ 0 − W χ(z) 2 dz0n(z0) , (3.3) = 2 c a χ(z) χ(z0) where Ωm denotes the matter density parameter today, H0 is the present-day Hubble parame-
56 3.4. Maps
ter, χh is the comoving distance to the horizon and a denotes the scale factor.
Similarly to the spherical harmonic power spectra of galaxy clustering and weak lensing the spherical harmonic power spectrum of CMB temperature anisotropies T can be related to the primordial matter power spectrum generated during inflation as [71]
Z ¯ ¯2 2 ¯∆T`(k)¯ C TT dk k2P lin(k)¯ ¯ , (3.4) ` = π δδ ¯ δ(k) ¯ where ∆T` denotes the transfer function of the temperature anisotropies and δ is the matter overdensity.
The CMB temperature anisotropies are correlated to tracers of the large-scale structure (LSS) such as galaxy overdensity and weak lensing shear primarily through the integrated Sachs- Wolfe effect [232]. On large enough scales where linear theory holds, the spherical harmonic power spectra between these probes can be computed from expressions similar to those above. In the Limber approximation [169, 147, 148], the spherical harmonic power spectrum between CMB temperature anisotropies and a tracer i of the LSS becomes: [62]
2 µ 1 ¶ ΩmH TCMB 1 Z d ` /2 C iT 3 0 dz [D(z)(1 z)]D(z)W i ¡χ(z)¢P lin k + ,0 , (3.5) ` 2 1 2 δδ = c (` /2) dz + = χ(z) + where T denotes the mean temperature of the CMB today, i {δ ,γ} and W i ¡χ(z)¢ repre- CMB ∈ g sents the window functions defined in Equations 3.2 and 3.3. We have further split the linear lin matter power spectrum Pδδ (k,z) into its time-dependent part parametrized by the growth lin factor D(z) and the scale-dependent part Pδδ (k,0). For a derivation of Eq. 3.5 for the galaxy γT overdensity field as tracer of the LSS see e.g. Padmanabhan et al. [197]. The derivation for C` is similar and is detailed in Appendix 3.A.
To compute the auto-power spectrum of the CMB temperature anisotropies we use the publicly available Boltzmann code CLASS6 [164]. For the other power spectra we use PYCOSMO [221]. We calculate the linear matter power spectrum from the transfer function derived by Eisenstein and Hu [78]. To compute the nonlinear matter power spectrum we use the HALOFIT fitting function [251] with the revisions of Takahashi et al. [267].
3.4 Maps
3.4.1 Cosmic Microwave Background
We use the foreground-reduced CMB anisotropy maps provided by the Planck collabora- tion [210] in their 2015 data release. We choose these over the uncleaned single-frequency maps because they allow to perform the foreground correction on the maps rather than the power spectrum level. This is important when considering probe combination. The Planck
6http : //class-code.net.
57 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
Table 3.1: Summary of used data.
CMB Survey: Planck 2015 [210] temperature Fiducial foreground-reduced map: Commander anisotropies Sky coverage: f 0.776 sky = Survey: SDSS DR8 [13] Sky coverage: f 0.27 galaxy sky = Galaxy sample: CMASS1-4 overdensity Number of galaxies: N 854063 gal = Photometric redshift range 0.45 z 0.65 ≤ phot < Survey: SDSS Stripe 82 co-add [20] Sky coverage: f 0.0069 sky = weak lensing Number of galaxies: N 3322915 gal = Photometric redshift range: 0.1 . zphot . 1.1 r.m.s. ellipticity per component: σ 0.43 e ∼
foreground-reduced CMB anisotropy maps have been derived using four different algorithms: Commander, NILC, SEVEM and SMICA. The maps are given in HEALPix7 [99] format and are provided in Galactic coordinates at two different resolutions of NSIDE 1024 and NSIDE = 2048. These correspond to pixel areas of 11.8 and 2.95 arcmin2 respectively. Different data = configurations are available [210]; we use both the half-mission half-sum (HMHS) maps, which contain both signal and noise, and the half-mission half-difference maps (HMHD), which contain only noise and potential residual systematic uncertainties. All four maps yield consistent estimates of both the spherical harmonic power spectrum of the CMB temperature anisotropies as well as the spherical harmonic cross-power spectrum between CMB tempera- ture anisotropies and tracers of the LSS [210, 215, 213]. Since the Planck collaboration found the Commander approach to be the preferred solution for studying the CMB anisotropies at large and intermediate angular scales, we also choose it for our analysis. Each of the four fore- ground reduction methods also provides a confidence mask inside which the CMB solution is trusted. Following the Planck collaboration [210], we adopt the union of the confidence masks for Commander, SEVEM and SMICA. This is referred to as the UT78 mask and covers 77.6% of the sky at a resolution of NSIDE 2048. To downgrade the mask to NSIDE 1024, we follow = = the description outlined in Planck Collaboration et al. [210]. The HMHS CMB anisotropy map derived using Commander is shown in the top panel of Fig. 3.2 for resolution NSIDE 1024 and = the corresponding HMHD map is shown in Fig. 3.18 in the Appendix.
7http : //healpix.sourceforge.net.
58 3.4. Maps
31.0 CMB temperature 4 10 ⇥ 4.5 60
30
120 60 0 60 120
b 33.5 0
30
60 4.5 36.0 55.5 53.0 50.5 T [K] l
31.0 Galaxy density 10 60
30
120 60 0 60 120
b 33.5 0
30
60 1 36.0 55.5 53.0 50.5 g [1] l
31.0 Weak lensing 0.1000 60
30
120 60 0 60 120
b 33.5 0
30
60 0.0003 36.0 ˆ [1] 55.5 53.0 50.5 | | l
Figure 3.2: Summary of the three maps in Galactic coordinates used in this analysis. The all-sky maps are in Mollweide projection while the zoom-in versions are in Gnomonic projection. The HMHS map of CMB temperature anisotropies as derived using Commander is shown in the top panel. It is masked using the UT78 mask. The middle panel shows the systematics-corrected (see text) galaxy overdensity map for CMASS1-4 galaxies. Grey areas have been masked either because they lie outside the survey footprint or are potentially contaminated by systematics. The lower panel shows the map of the SDSS Stripe 82 shear modulus γˆ . Grey areas have been | | masked because they are either unobserved or do not contain galaxies for shear measurement. The zoom-in figures (left) are enlarged versions of the 5 5 deg2 region centered on (l,b) × (53°, 33.5°) shown in the maps. The zoom-in for the galaxy shear map is overlaid with a = − whisker plot of the galaxy shears. All three maps have resolution NSIDE 1024. =
59 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
3.4.2 Galaxy overdensity
The SDSS [283, 79, 102, 103] obtained wide-field images of 14555 deg2 of the sky in 5 photo- metric passbands (u,g,r,i,z [89, 250, 72]) up to a limiting r-band magnitude of r 22.5. The ' photometric data is complemented with spectroscopic data from the Baryonic Oscillations Spectroscopic Survey (BOSS) [79, 65, 249]. BOSS was conducted as part of SDSS III [79] and obtained spectra of approximately 1.5 million luminous galaxies distributed over 10000 deg2 of the sky. The SDSS photometric redshifts for DR8 [13] are estimated using a local regression model trained on a spectroscopic training set consisting of 850000 SDSS DR8 spectra and spectroscopic data from other surveys8. The algorithm is outlined in Beck et al. [28].
In our analysis, which is described in the following, we largely follow Ho et al. [129]. We select objects classified as galaxies from the PHOTOPRIMARY table in the Catalog Archive Server (CAS9). To obtain a homogeneous galaxy sample we further select CMASS galaxies using the color-magnitude cuts used for BOSS target selection [79] and outlined in Ho et al. [129]. This selection isolates luminous, high-redshift galaxies that are approximately stellar mass limited [280, 229]. We further restrict the sample to CMASS galaxies with SDSS photometric redshifts between 0.45 z 0.65, i.e. we consider the photometric redshift slices CMASS1-4. This ≤ < selection yields a total of N 1096455 galaxies. gal = To compute the galaxy overdensity field, we need to characterize the full area observed by the survey and mask regions heavily affected by foregrounds or potential systematics. The area imaged by the SDSS is divided into units called fields. Several such fields have been observed multiple times in the SDSS imaging runs. The survey footprint is the union of the best observed (primary) fields at each position and is described in terms of MANGLE [106, 107, 263] spherical polygons. Each of these polygons is matched to the SDSS field fully covering it10. In order to select the survey area least affected by foregrounds and potential systematics we follow Ho et al. [129] and Ross et al. [229] and restrict the analysis to polygons covered by fields with score11 0.6, full width at half maximum (FWHM) of the point spread function (PSF) ≥ PSF-FWHM 2.0 arcsec in the r-band and Galactic extinction E(B V) 0.08 as determined < − ≤ from the extinction maps from Schlegel et al. [235].
To facilitate a joint analysis between the LSS probes and the CMB, which is given as a map in Galactic coordinates, we transform both the galaxy positions as well as the survey mask from equatorial (RA, DEC) to Galactic (l, b) coordinates. We construct the continuous galaxy overdensity field by pixelizing the galaxy overdensities δ δn/n¯ onto a HEALPix pixelization g = of the sphere with resolution NSIDE 1024. We mask the galaxy overdensity map with a = HEALPix version of the SDSS survey mask, which is obtained by random sampling of the MANGLE mask. To account for the effect of bright stars, we use the Tycho astrometric catalog [131] and define magnitude-dependent stellar masks as defined in Padmanabhan et al. [198].
8More details can be found on http://www.sdss3.org/dr8/algorithms/photo-z.php. 9The SDSS Catalog Archive Server can be accessed through http://skyserver.sdss.org/CasJobs/SubmitJob.aspx. 10This information is found in the files window_unified.fits and window_flist.fits. 11http://www.sdss3.org/dr10/algorithms/resolve.php.
60 3.4. Maps
2.0 W γ(χ(z)) γ 1.5 n (z) nδg (z) ) z
( 1.0 n
0.5
0.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 z
Figure 3.3: Redshift distribution for the LSS probes. The figure shows the redshift selection function of SDSS CMASS1-4 galaxies, the redshift selection function for the SDSS Stripe 82 galaxies as well as the weak lensing shear window function defined in Eq. 3.3. The redshift selection function for CMASS1-4 galaxies as well as the weak lensing shear window function have been rescaled relative to the Stripe 82 redshift selection function.
We remove galaxies inside the bright star masks and correct for the area covered by the bright stars by removing the area covered by the star masks from the pixel area A pix,corr = A A when computing the galaxy overdensity. The final map covers a fraction pix,uncorr − stars f 0.27 of the sky and contains N 854063 galaxies. sky ≈ gal = Even after masking and removal of high contamination regions, there are still systematics left in the galaxy overdensity map. The correction for residual systematic uncertainties in the maps follows Ross et al. [229] and Hernández-Monteagudo et al. [119] and is described in Appendix 3.B. The final map is shown in the middle panel of Fig. 3.2.
As well as the maps we need an estimate for the redshift distribution of the galaxies in our sample. To this end we follow Ho et al. [129] and match photometrically detected galaxies to galaxies observed spectroscopically in SDSS DR9 [12]. We then estimate the redshift distribu- tion of the photometric galaxies from the spectroscopic redshift distribution of the matching galaxies. The selected CMASS1-4 galaxies have spectroscopic redshifts 0.4 . z . 0.7 as can be seen from the redshift distribution shown in Fig. 3.3.
3.4.3 Weak lensing
We take weak lensing data from the SDSS Stripe 82 co-add [20], which comprises 275 deg2 of co-added SDSS imaging data with a limiting r-band magnitude r 23.5 and r-band median ≈ seeing of 1.1 arcsec. The shapes of objects detected in the SDSS were measured from the
61 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing adaptive moments [32] by the PHOTO pipeline [176] and are available on the CAS12. Photo- metric redshifts for all detected galaxies were computed using a neural network approach as described in Reis et al. [224] and are available as a DR7 value added catalog13.
In the following analysis we closely follow the work by Lin et al. [170]. We select objects identified as galaxies in the co-add data (i.e. run 106 or run 206) from the CAS and = = we restrict the sample to galaxies with extinction corrected i-band magnitudes in the range 18 i 24. Further we select only objects that pass the clean photometry cuts as defined < < by the SDSS14 and do not have flags indicating problems with the measurement of adaptive moments as well as negative errors on those. The former cuts especially exclude galaxies containing saturated pixels. We use shapes measured in the i-band since it has the smallest seeing (1.05 arcsec) [20, 170] and further consider only galaxies with observed sizes at least 50% larger than the PSF. This requirement is quantified by requiring the resolution factor R 1 mRrCcPSF/mRrCc [32] to satisfy R 0.33, where mRrCc and mRrCcPSF denote the sum = − > of the second order moments in the CCD column and row direction for both the object and the PSF.
For the above galaxy sample we compute PSF-corrected galaxy ellipticities using the linear PSF correction algorithm as described in Hirata and Seljak [125]. For weak lensing shear measurement we follow Lin et al. [170] and restrict the sample to galaxies with PSF-corrected ellipticity components e ,e satisfying e 1.4 as well as e 1.4 and photometric redshift 1 2 | 1| < | 2| < uncertainties σ 0.15. This sample has an r.m.s. ellipticity per component of σ 0.43. We z < e ∼ then turn the PSF-corrected ellipticities for this sample into shear estimates. The details of the analysis are described in Appendix 3.C.
After computing weak lensing shear estimates from the ellipticities we apply a rotation to both the galaxy positions and shears from equatorial to Galactic coordinates15 to allow for combination with the CMB. We pixelize both weak lensing shear components onto separate HEALPix pixelizations of the sphere choosing a resolution of NSIDE 1024 as for the galaxy = overdensity map. At this resolution the mean number of galaxies per pixel is about 38, which 2 corresponds to n 3.2 arcmin− . We apply a mask to both maps, which accounts for both gal ' unobserved and empty pixels. The final maps are constructed using N 3322915 galaxies gal = and cover a sky fraction f 0.0069. The map of the shear modulus γˆ is shown in the sky ≈ | | bottom panel of Fig. 3.2 together with a zoom-in region with overlaid whisker plot illustrating the magnitude and direction of the weak lensing shear.
We follow Lin et al. [170] and estimate the redshift distribution of the galaxies from their photometric redshift distribution. The redshift distribution is shown in Fig. 3.3 together with the window function for weak lensing shear defined in Eq. 3.3. We see that the selected galaxies have photometric redshifts z . 1.0.
12See footnote9. 13http : //classic.sdss.org/dr7/products/value_added/, http : //das.sdss.org/va/coadd_galaxies/. 14http://www.sdss.org/dr12/tutorials/flags/. 15The exact rotation of the shears is described in Appendix 3.D.
62 3.5. Spherical harmonic power spectra
3.5 Spherical harmonic power spectra
We calculate the spherical harmonic power spectra of the maps presented in the previous section using the publicly available code PolSpice16 [265, 57]. The PolSpice code is designed to combine both real and Fourier space in order to correct spherical harmonic power spectra measured on a cut-sky from the effect of the mask. The algorithm can be summarized as follows: starting from a masked HEALPix map, PolSpice first computes the so-called pseudo power spectrum, which is then Fourier transformed to the real space correlation function. In order to correct for the effects of the mask, the latter is divided by the mask correlation function. In a last step, the demasked correlation function is Fourier transformed back to the spherical harmonic power spectrum. This approach ensures that PolSpice can exploit the advantages of real space while still performing the computationally expensive calculations in Fourier space.
Demasking can only be performed on angular scales on which information is available, which translates to a maximal angular scale θmax for which a demasked correlation function can be computed. This maximal scale leads to ringing when transforming back from real to Fourier space, which can be reduced by apodizing the correlation function prior to inversion. Both these steps lead to biases in the power spectrum recovered by PolSpice. The kernels relating the average PolSpice estimates to the true power spectra can be computed theoretically for a given maximal angular scale and apodization prescription and need to be corrected for when comparing theoretical predictions to observed power spectra.
An additional difficulty arises in the computation of spherical harmonic power spectra of spin- 2 fields. Finite sky coverage tends to cause mixing between E-and B-modes. The polarization version of PolSpice is designed to remove E- to B-mode leakage in the mean [57]. Details on our earlier application of PolSpice to LSS data are described in Appendix A of Becker et al. [29].
In order to calculate both the auto- and cross-power spectra for all probes, we need to estimate the maximal angular scale θmax. This is not a well-defined quantity but we can separately estimate it for each probe from the real space correlation function of its mask. The real space correlation function of the survey mask will fall off significantly or vanish for scales larger than
θmax. We therefore estimate θmax as the scale around which the mask correlation function significantly decreases in amplitude. Appendix 3.E illustrates this analysis for the example of the SDSS Stripe 82 weak lensing shear mask. In order to reduce Fourier ringing we apodize the correlation function using a Gaussian window function; following Chon et al. [57] we choose
θ the FWHM of the Gaussian window as θ max/2. Survey masks with complicated angular FWHM = dependence might not exhibit a clear fall-off, which complicates the choice of θmax. We therefore validate our choices of θmax and θFWHM with the Gaussian simulations as described in Appendix 3.F and 3.G. We find our choices to allow the recovery of the input power spectra for all the probes and settings.
16http://www2.iap.fr/users/hivon/software/PolSpice/.
63 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
Table 3.2: Spherical harmonic power spectrum parameters and angular multipole ranges.
Power spectrum θmax [deg] θFWHM [deg] `-range ∆` TT C` 40 20 [10, 610] 30 δg δg C` 80 40 [30, 210] 30 γγ C` 10 5 [70, 610] 60 δg T C` 40 20 [30, 210] 30 γT C` 10 5 [70, 610] 60 γδg C` 10 5 [30, 210] 60
All spherical harmonic power spectra are corrected for the effect of the HEALPix pixel window function and the power spectra involving the CMB map are further corrected for the Planck effective beam window function, which complements the CMB maps.
We now separately describe the measurement of all the six spherical harmonic power spec- tra. To compute the power spectra, we use the maps and masks described in Section 3.4 at resolution NSIDE 1024, except for the CMB temperature power spectrum. For the latter = we use the maps at resolution NSIDE 2048, but we do not expect this to make a significant = difference. The PolSpice parameter settings used to compute the power spectra are sum- marized in Tab. 3.2. This table further gives the angular multipole range as well as binning scheme employed for the cosmological analysis. For all probes considered, the uncertainties are derived from the Gaussian simulations described in Section 3.6.2 and Appendix 3.F.
3.5.1 CMB
We use the half-mission half-sum (HMHS) map to estimate the CMB signal power spectrum and the half-mission half-difference (HMHD) map to estimate the noise in the power spectrum of the HMHS map.
The minimal angular multipole used in the cosmological analysis is chosen such as to min- imize demasking effects and the cut at ` 610 ensures that we are not biased by residual = foregrounds in the maps as discussed in Section 3.7. The resulting power spectrum is shown in the top panel of Fig. 3.4. In Appendix 3.H we compare the CMB auto-power spectrum computed from the different foreground-reduced maps. As illustrated in Fig. 3.16 in Appendix 3.H we find that the measured CMB auto-power spectrum is unaffected by the choice of foreground-reduced map.
64 3.5. Spherical harmonic power spectra
6000 ] 2 TT K 5000 µ )[
π 4000 (2 / 3000 TT ` C 2000 + 1)
` 1000 ( ` 0 0 100 200 300 400 500 600 7 10− 4 1.0 × 10− 0.8 δgT δgδg 0.6 g [K]
0.4 δ g T δ ` g δ ` 0.2 C `C 0.0
0.2 5 − 10− 0.4 − 50 100 150 200 3 101 102 9 ×5 10− 10− 2.5 × 1.5 × 2.0 γT γδg γγ 1.0 ) 5 1.5 π 10− (2 /
g 0.5 [K] γγ 1.0 ` γδ ` T C γ ` 0.5 `C 0.0 6 `C 10− + 1)
0.0 ` ( 0.5 ` 0.5 − − 7 1.0 1.0 10− − 100 200 300 400 500 600 − 50 100 150 200 102 6 102 × ` ` `
Figure 3.4: Spherical harmonic power spectra for all probes used in the cosmological analy- sis. The top left panel shows the power spectrum of CMB anisotropies computed from the Commander CMB temperature map at resolution of NSIDE 2048. The middle left panel shows = the cross-power spectrum between CMB temperature anisotropies and galaxy overdensity computed from the systematics-reduced SDSS CMASS1-4 map and the Commander map at res- olution NSIDE 1024. The middle right panel shows the spherical harmonic power spectrum = of the galaxy overdensity computed from the systematics-reduced SDSS CMASS1-4 map at NSIDE 1024. The bottom left panel shows the spherical harmonic power spectrum between = CMB temperature anisotropies and weak lensing shear measured from the Commander CMB map and the SDSS Stripe 82 weak lensing maps at resolution NSIDE 1024. The bottom- = middle panel shows the spherical harmonic power spectrum between galaxy overdensity and galaxy weak lensing shear computed from the systematics-reduced SDSS CMASS1-4 map and the SDSS Stripe 82 galaxy weak lensing shear map at resolution NSIDE 1024. The bottom = right panel shows the spherical harmonic power spectrum of cosmic shear E-modes computed from the SDSS Stripe 82 weak lensing shear maps. The angular multipole ranges and binning schemes for all power spectra are summarized in Table 3.2. All power spectra are derived from the maps in Galactic coordinates. The solid lines show the theoretical predictions for the best-fit cosmological model determined from the joint analysis which is summarized in Tab. 3.3. The theoretical predictions have been convolved with the PolSpice kernels as described in Section 3.5. The error bars are derived from the Gaussian simulations described in Section 3.6.2 and Appendix 3.F.
65 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
3.5.2 Galaxy clustering
The galaxy overdensity maps described in Section 3.4 are estimated from discrete galaxy tracers. Therefore, their spherical harmonic power spectrum receives contributions from the galaxy clustering signal and Poisson shot noise. To estimate the noise power spectrum, we resort to simulations. We generate noise maps by randomizing the positions of all the galaxies in the sample inside the mask. Since this procedure removes all correlations between galaxy positions, the power spectra of these maps will give an estimate of the level of Poisson shot noise present in the data. In order to obtain a robust noise power spectrum, we generate 100 noise maps and estimate the noise power spectrum from the mean of these power spectra.
The spherical harmonic galaxy clustering power spectrum contains significant contributions from nonlinear structure formation at small angular scales. The effects of nonlinear galaxy bias are difficult to model and we therefore restrict our analysis to angular scales for which nonlinear corrections are small. We can estimate the significance of nonlinear effects by comparing the spherical harmonic galaxy clustering power spectrum computed using the nonlinear matter power spectrum as well as the linear matter power spectrum. Since galaxies are more clustered than dark matter this is likely to underestimate the effect. We find that the difference between the two reaches 5% of the power spectrum uncertainties and thus becomes mildly significant at around ` 250. This difference is smaller than the difference max ∼ derived in Ho et al. [129] and de Putter et al. [66] which is likely due to the fact that we consider a single redshift bin and do not split the data into low and high redshifts. In order not to bias our results we choose ` 210 which is comparable to the limit used in Ho et al. [129] and max = de Putter et al. [66]. To determine the minimal angular multipole we follow Ho et al. [129], who determined that the Limber approximation becomes accurate for scales larger than ` 30. = The middle right panel in Fig. 3.4 shows the spherical harmonic galaxy clustering power spectrum computed from the systematics-corrected map in Galactic coordinates. In Appendix 3.H, we compare the spherical harmonic power spectrum derived from the systematics- corrected maps in Galactic and equatorial coordinates. We find small differences at large angular scales, but the effect on the bandpowers considered in this analysis is negligible, as can be seen from Appendix 3.H (Fig. 3.14). To test the procedure for removing systematic uncertainties, we compare the spherical harmonic power spectra before and after correcting the maps for residual systematics. We find that the removal of systematics marginally reduces the clustering amplitude on large scales, which is expected since Galactic foregrounds exhibit significant large scale clustering. Small angular scales on the other hand, are mostly unaffected by the corrections applied. These results are shown in Appendix 3.H (Fig. 3.17).
3.5.3 Cosmic shear
The power spectrum computed from the weak lensing shear maps described in Section 3.4.3 contains contributions from both the cosmic shear signal and the shape noise of the galaxies, which is due to intrinsic galaxy ellipticities. In order to estimate the shape noise power
66 3.5. Spherical harmonic power spectra spectrum we follow the same methodology as for galaxy clustering and resort to simulations. We generate noise-only maps by rotating the shears of all the galaxies in our sample by a random angle. This procedure removes spatial correlations between galaxy shapes. Since the weak lensing shear signal is at least an order of magnitude smaller than the intrinsic galaxy ellipticities, the power spectrum of the randomized map gives an estimate of the shape noise power spectrum. As for galaxy clustering, we compute 100 noise maps and estimate the shape noise power spectrum from the mean of these 100 noise power spectra.
For the cosmological analysis we choose broader multipole bins than for the CMB temperature anisotropies and galaxy clustering since the small sky fraction covered by SDSS Stripe 82 causes the cosmic shear power spectrum to be correlated across a significantly larger multipole range. The low and high ` limits are chosen to minimize demasking uncertainties and the impact of nonlinearities in the cosmic shear power spectrum.
The spherical harmonic power spectrum of the weak lensing shear E-mode is displayed in the bottom right panel of Fig. 3.4 and the B-mode power spectrum is shown in the Appendix (Fig. 3.19). We see that the E-mode power spectrum is intrinsically low as compared to the best-fit theory power spectrum. These results are similar to those derived by Lin et al. [170], 0.7 who found a low value of Ωm σ8 for Stripe 82 cosmic shear. As can be seen, we do not detect a significant B-mode signal.
When comparing the weak lensing shear E-mode power spectra computed from the maps in Galactic and equatorial coordinates, we find discrepancies. These are mainly caused by the correction for additive bias in the weak lensing shears. As described in Appendix 3.C, the PSF-corrected galaxy shears are affected by an additive bias. Following Lin et al. [170], we correct for this bias by subtracting the mean shear of each CCD camera column from the galaxy shears. This correction is performed in equatorial coordinates and ensures that the mean shear vanishes in this coordinate system. When the galaxy positions and shears are rotated from equatorial to Galactic coordinates, this ceases to be true. Therefore the correction for additive bias is coordinate-dependent and it is this effect that causes the main discrepancies between the measured power spectra. Further descriptions of the impact of the additive shear bias correction can be found in Appendix 3.H.1.
The discrepancies between the cosmic shear power spectra measured from maps in Galactic and equatorial coordinates are still within the experimental uncertainties. We therefore choose to correct for the additive shear bias in equatorial coordinates, apply the rotation to the corrected shears and compute the cosmic shear power spectrum from the maps in Galactic coordinates. We note however, that these differences will become significant for surveys measuring cosmic shear with higher precision. It is therefore important to develop coordinate-independent methods for shear bias correction when performing a joint analysis of different cosmological probes.
67 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
3.5.4 CMB and galaxy overdensity cross-correlation
To compute the spherical harmonic cross-power spectrum between CMB temperature anisotropies and the galaxy overdensity, we use the maps and masks described in Sections 3.4.1 and 3.4.2.
We generally have two possibilities to compute cross-correlations between two maps with different angular masks. We can either compute the cross-correlation by keeping the respective mask for each probe, or we can compute a combined mask, which is the union of all pixels masked in at least one of the maps. When testing both these cases on Gaussian simulations, we observed a better recovery of the input power spectra when applying the combined mask to both maps. We therefore mask both maps with the combined mask, which covers a fraction of sky f 0.26. sky ∼ The spherical harmonic cross-power spectrum between CMB temperature anisotropies and galaxy overdensity is shown in the middle left panel of Fig. 3.4. We see that the ISW power spectrum is very noisy, which makes its detection significance small. Since the power spectrum uncertainties for the considered angular scales are mainly due to cosmic variance, we suspect that the low signal-to-noise is mainly due to the fraction of sky covered by the SDSS CMASS1-4 galaxies. Despite its low significance, we include the ISW power spectrum in our analysis, because we expect it to help break degeneracies between cosmological parameters. We check that the ISW power spectrum does not depend on the choice of foreground-reduced CMB map. We find that the results using the maps provided by the NILC, SEVEM and SMICA algorithms are virtually the same, as illustrated in Appendix 3.H (Fig. 3.16).
3.5.5 CMB and weak lensing shear cross-correlation
We estimate the spherical harmonic cross-power spectrum between CMB temperature anisotropies and the weak lensing shear E-mode field from the maps and masks described in Sections 3.4.1 and 3.4.3. Both maps are masked with the combination of the masks, which covers a fraction of sky f 0.0065. sky ∼ The bottom left panel in Fig. 3.4 shows the spherical harmonic power spectrum between CMB temperature anisotropies and the weak lensing shear E-mode field. As can be seen, the noise level is too high to allow for a detection of the ISW correlation between CMB temperature anisotropies and weak lensing shear. This is to be expected due to the small sky fraction covered by the SDSS Stripe 82 galaxies and the intrinsically low signal-to-noise of this cross- correlation. Nevertheless, we include the power spectrum into the joint analysis to provide an upper limit on the ISW from weak lensing. The measured power spectrum is unaffected by the choice of CMB mapmaking method, as illustrated in Fig. 3.16 in Appendix 3.H.
68 3.6. Covariance matrix
3.5.6 Galaxy overdensity and weak lensing shear cross-correlation
We compute the spherical harmonic cross-power spectrum between the galaxy overdensity and weak lensing shear E-mode field from the maps and masks described in Sections 3.4.2 and 3.4.3. We mask both maps with the combination of the two masks. The combined mask covers a sky fraction f 0.0053. sky ∼ The spherical harmonic cross-power spectrum between galaxy overdensity and weak lensing shear E-mode is shown in the bottom-middle panel of Fig. 3.4. We see that the signal-to-noise of the power spectrum is low at the angular scales considered. This is probably due to the small sky fraction covered by Stripe 82 galaxies. We nevertheless include this cross-correlation in our analysis to serve as an upper limit. In Appendix 3.H we show the comparison between the power spectra measured from the maps in Galactic and in equatorial coordinates. We find reasonable agreement between the two, even though the discrepancies are significantly enhanced compared to the effects on the galaxy overdensity power spectrum. As discussed in Section 3.5.3 this is probably due to the coordinate-dependence of the additive shear bias correction.
3.6 Covariance matrix
In order to obtain cosmological constraints from a joint analysis of CMB temperature anisotropies, galaxy clustering and weak lensing we need to estimate the joint covariance matrix of these cosmological probes. In this work we assume all the fields to be Gaussian random fields, i.e. we assume the covariance between all probes to be Gaussian and neglect any non-Gaussian contribution. This is appropriate for the CMB temperature field as well as the galaxy over- density field at the scales considered but it is only an approximation for the weak lensing shear field [233]. For example, for a survey with source redshifts z 0.6, Sato et al. [233] s = found that neglecting non-Gaussian contributions leads to an underestimation of the diagonal terms in the cosmic shear covariance matrix by a factor of approximately 5 at multipoles ` 600. In our case the discrepancy may be more pronounced since our sample contains a ∼ significant number of galaxies with z 0.6. On the other hand we will be less sensitive to s < the non-Gaussian nature of the covariance matrix since the covariance for our galaxy sample is dominated by shape noise especially at the highest multipoles considered. We therefore decide to leave the introduction of non-Gaussian covariance matrices to future work.
In this work, we employ two different models for the joint Gaussian covariance matrix CG : the first is a theoretical model and the second is based on simulations of correlated Gaussian realizations of the three cosmological probes. We use the theoretical covariance matrix to validate the covariance matrix obtained from the simulations.
69 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
3.6.1 Theoretical covariance estimate
The covariance between cosmological spherical harmonic power spectra is composed of two parts: cosmic variance and noise. For spherical harmonic power spectra computed over the full sky, different ` modes are uncorrelated and the covariance matrix is diagonal. Partial sky coverage, i.e. f 1, has the effect to couple different ` modes and thus leads to a non- sky < diagonal covariance matrix. This covariance becomes approximately diagonal if it is binned into approximately uncorrelated bandpowers of width ∆` [42]. Cabré et al. [42] found the empirical relation ∆`fsky 2 to be a good approximation. In this case the covariance matrix ∼ i j i j between binned power spectra C and C 0 0 can be approximated as [133,42,76] ` `0 h i j i 0 j 0 i j i 0 j 0 δ``0 ii 0 ii 0 j j 0 j j 0 CovG (C ,C ) ∆C ∆C (C N )(C N ) ` `0 = 〈 ` `0 〉' (2` 1)∆`f ` + ` + + sky (3.6) i j ji i (C 0 N i j 0 )(C 0 N ji 0 ) , + ` + ` + where i, j, i 0, j 0 denote different cosmological probes; in our case i, j, i 0, j 0 {T,δ ,γ}. The ∈ g quantities N i j are the noise power spectra of the different probes, which vanish unless i j. = Given a cosmological model and survey specifications such as fractional sky coverage and noise level, we can approximate CG using Eq. 3.6 for each block covariance matrix. We choose a hybrid approach: we adopt a cosmological model to calculate the signal power spectra whereas we approximate N i j with the measured noise power spectra used to remove the noise bias in the data as described in Section 3.5.
3.6.2 Covariance estimate from Gaussian simulations
The theoretical covariance matrix estimate described above is expected to only yield accurate results for uncorrelated binned power spectra, since in this approximation the covariance matrix is fully diagonal. For this reason we also estimate the covariance matrix in an alternative way that does not rely on this approximation: we estimate an empirical covariance matrix from the sample variance of Gaussian simulations of the three cosmological probes. To this end, we simulate correlated realizations of both the two spin-0 fields, CMB temperature and galaxy overdensity, as well as the spin-2 weak lensing shear field. We follow the approach outlined in Giannantonio et al. [93] for simulating correlated maps of spin-0 fields and we make use of the polarization version of the HEALPix routine synfast to additionally simulate correlated maps of the spin-2 field. We estimate noise maps from the data and add these to the correlated signal maps. The details of the algorithm are outlined in Appendix 3.F.
In order to compute the power spectrum covariance matrix, we apply the masks used on the data to the simulated maps and calculate both the auto- and the cross-power spectra of all the probes using the same methodology and PolSpice settings as described in Section 3.5. We
70 3.7. Cosmological constraints
TT Tδg δgδg Tγ δgγ γγ 610
0.90 γγ 0.45 0.15 70 210 δ γ 30 g 610
Tγ
70 210 δgδg 0
` 10 210 Tδg 0.00 10 610
TT
0.15 10 − 10 610 210 210 610 210 610 10 10 70 30 70 `
Figure 3.5: Correlation matrix for the spherical harmonic power spectra derived from the sample variance of the Gaussian simulations. The binning scheme and angular multipole range for each probe follow those outlined in Tab. 3.2.
generate Nsim random realizations and estimate the covariance matrix as
Nsim i j i 0 j 0 1 X h i j i j ih i 0 j 0 i 0 j 0 i CovG (C ,C ) C (`) C¯ (`) C (`0) C¯ (`0) , (3.7) ` `0 k k k k = Nsim 1 k 1 − − − = ¯i j where Ck (`) denotes the mean over all realizations.
The accuracy of the sample covariance estimate depends on the number of simulations. As described in Cabré et al. [42], N 1000 achieves better than 5% accuracy for estimating sim = the covariance matrix for the ISW effect from Gaussian simulations. We therefore follow Cabré et al. [42] and compute the covariance matrix from the sample variance of N 1000 sim = Gaussian realizations of the 4 maps or 6 spherical harmonic power spectra respectively.
The correlation matrix for the spherical harmonic power spectra derived from the Gaussian simulations for binning schemes and angular multipole ranges described in Section 3.5 is shown in Fig. 3.5. We see that the survey masks lead to significant correlations between bandpowers.
3.7 Cosmological constraints
Each of the power spectra presented in Section 3.5 carries cosmological information with probe-specific sensitivities and degeneracies. An integrated combination of these cosmologi-
71 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing cal probes therefore helps break these parameter degeneracies. It further provides robust cos- mological constraints since it is derived from a joint fit to the auto- as well as cross-correlations of three cosmological probes.
In order to derive cosmological constraints from a joint fit to the six spherical harmonic power spectra discussed in Section 3.5, we assume the joint likelihood to be Gaussian, i.e.
1 1 obs theor T 1 obs theor 2 (C` C` ) CG− (C` C` ) L (D θ) 1 e− − − , (3.8) d /2 | = [(2π) detCG ]
theor where CG denotes the Gaussian covariance matrix. C` denotes the theoretical prediction obs for the spherical harmonic power spectrum vector of dimension d and C` is the observed power spectrum vector, defined as
³ δ T δ δ T γδ ´ Cobs C TT C g C g g C γ C g C γγ . (3.9) ` = ` ` ` ` ` ` obs A Gaussian likelihood is a justified assumption for both the CMB temperature anisotropy and galaxy clustering power spectra due to the central limit theorem. Since the weak lensing shear power spectrum receives significant contribution from non-linear structure formation, its likelihood will deviate from being purely Gaussian [110]. It has been shown however, that a Gaussian likelihood is a sensible approximation, especially when CMB data is added to weak lensing [234]. In our first implementation we will thus assume both a joint Gaussian likelihood and Gaussian single probe likelihoods.
We estimate the covariance matrix using both methods outlined in Section 3.6. In both cases we compute the covariance for a ΛCDM cosmological model, which we keep fixed in the joint fit. Note that the covariance matrices depend on the cosmological model and should therefore vary in the fitting procedure [75]. Following standard practice, (e.g. [8]), we approximate the covariance matrix to be constant and compute it for a ΛCDM cosmological model with parameter values {h, Ω , Ω , n , σ , τ , T } {0.7, 0.3, 0.049, 1.0, 0.88, 0.078, 2.275K}, m b s 8 reion CMB = where h is the dimensionless Hubble parameter, Ωm is the fractional matter density today, Ωb is the fractional baryon density today, ns denotes the scalar spectral index, σ8 is the r.m.s. of 1 matter fluctuations in spheres of comoving radius 8h− Mpc and τreion denotes the optical depth to reionization. We further set the linear, redshift-independent galaxy bias parameter to b 2. To obtain an unbiased estimate of the inverse of the covariance matrix derived from the = Gaussian simulations, we apply the correction derived in Kaufman [150], Hartlap et al. [109] and Anderson [19], i.e. we multiply the inverse covariance matrix by (N d 2)/(N 1). sim − − sim − The theoretical covariance matrix estimate does not suffer from this bias and is thus left unchanged.
From the likelihood given in Eq. 3.8, we derive constraints in the framework of a flat ΛCDM cos- mological model, where our fiducial model includes one massive neutrino eigenstate of mass
0.06 eV as in [212]. Our parameter set consists of the six ΛCDM parameters {h, Ωm, Ωb, ns, σ8, τreion}. We further marginalize over two additional parameters: a redshift independent, linear galaxy
72 3.7. Cosmological constraints
Table 3.3: Parameters varied in the MCMC with their respective priors and posterior means. The uncertainties denote the 68% c.l..
Parameter Prior Posterior mean h flat [0.2, 1.2] 0.699 0.018 ∈ ± Ω flat [0.1, 0.7] 0.278 0.019 m ∈ +0.020 Ω flat [0.01, 0.09] 0.0455 −0.0018 b ∈ ± n flat [0.1, 1.8] 0.975 0.019 s ∈ +0.018 σ flat [0.4, 1.5] 0.799 −0.029 8 ∈ ± τ Gaussian with µ 0.089, σ 0.02a 0.0792 0.0196 reion = = ± b flat [1., 3.] 2.13 0.06 ∈ ± 0.080 m Gaussian with µ 0.0, σ 0.1 0.142+0.081 = = − − a This corresponds to a WMAP9 [124] prior with increased variance to accom- modate the Planck results. bias parameter b and a multiplicative bias parameter m for the weak lensing shear. The multiplicative bias parametrizes unaccounted calibration uncertainties affecting the weak lensing shear estimator γˆ and is defined as [120]
γˆ (1 m)γ. (3.10) = + We note that we do not include additional nuisance parameters such as additive weak lensing shear bias, stochastic and scale-dependent galaxy bias [269, 202, 68], photometric redshift uncertainties, intrinsic galaxy alignments (for reviews, see e.g. [273, 141]) or parameters describing the effect of unresolved point sources on the CMB temperature anisotropy power spectrum [208]. In this present work we restrict the analysis to angular scales where these effects are expected to be subdominant.
We sample the parameter space with a Monte Carlo Markov Chain (MCMC) using CosmoHammer [15]. The parameters sampled are summarized in Table 3.3 along with their priors. We choose
flat, uniform priors for all parameters except for τreion and m. The optical depth to reion- ization can only be constrained with CMB polarization data. Since we do not include CMB polarization in this analysis, we apply a Gaussian prior with µ 0.089 and σ 0.02 on τ . = = reion This corresponds to a WMAP9 [124] prior with increased variance to accommodate the Planck 2015 results [212]. We further apply a Gaussian prior on the multiplicative bias m with mean µ 0 and σ 0.1. This is motivated by Hirata and Seljak [125], who found the multiplicative = = bias for the linear PSF correction method to lie in the range m [ 0.08,0.13] for the sample ∈ − considered in this analysis.
In our fiducial configuration presented below we use the covariance matrix derived from the Gaussian simulations as described in Section 3.6.2. We find that this choice does not influence our results since the constraints derived using the theoretical covariance are consistent. In order to further assess the impact of a cosmology-dependent covariance matrix, we perform
73 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing the equivalent analysis using a covariance matrix computed with a cosmological model with 5% lower σ . We find that the derived parameter values change by at most 0.5σ. The width ∼ 8 of the contours is only marginally changed.
In addition to the joint analysis, we also derive parameter constraints from separate analyses of TT δg δg γγ the three auto-power spectra C` ,C` and C` . In all three cases we assume a Gaussian like- lihood as in Eq. 3.8 and derive constraints on the base ΛCDM parameters {h, Ωm, Ωb, ns, σ8} as well as additional parameters constrained by each probe. These are τreion for the CMB temperature anisotropies, b for galaxy clustering and m for the cosmic shear.
Fig. 3.6 shows the constraints on the ΛCDM parameters {h, Ωm, Ωb, ns, σ8} derived from the joint analysis using the spherical harmonic power spectrum vector and likelihood defined in
Equations 3.8 and 3.9. These have been marginalized over τreion, b and m. Also shown are TT δg δg the constraints derived from separate analyses of the three auto-power spectra C` ,C` γγ and C` , each of them marginalized over the respective nuisance parameter. As expected, we find that the constraints derived from the CMB anisotropies are the strongest, followed by the galaxy clustering and cosmic shear constraints, which both constrain the full ΛCDM model rather weakly. The constraints from the CMB temperature anisotropies are broader and have central values which differ from those derived in Planck Collaboration et al. [212]. The reason for these discrepancies is the limited angular multipole range ` [10, 610] employed in the ∈ CMB temperature analysis. This causes the CMB posterior to become broader, asymmetric and results in a shift of the parameter means. We have verified that the Planck likelihood and our analysis give consistent results when the latter is restricted to a similar `-range. If on the other hand, we increase the high multipole limit to ` 1000, we find significant differences max = between our analysis and the Planck likelihood. We therefore choose to be conservative and use ` 610 throughout this work. Comparing the single probe constraints to one another max = we see that they agree reasonably well, the only slight discrepancy being the low value of both
Ωm and σ8 derived from the cosmic shear analysis. This is similar to the results derived in Lin et al. [170] even though the values for Ωm and σ8 are even lower in our analysis. However, care must be taken since the amplitude of the cosmic shear auto-power spectrum appears to have a small dependence on the choice of the coordinate system as discussed in Appendix 3.H.
The potential of the joint analysis emerges when the three auto-power spectra are combined together with their three cross-power spectra. Due to the complementarity of the differ- ent probes the constraints tighten and the allowed parameter space volume is significantly reduced. This is especially true in our case, since the constraints from CMB temperature anisotropies are broadened due to the restricted multipole range that we employed. Including more CMB data would significantly reduce the impact of adding additional cosmological probes. The numerical values of the best-fit parameters and their 68% confidence limits (c.l.) derived from the joint analysis are given in Tab. 3.3.
Fig. 3.7 compares the constraints derived from the joint analysis to the constraints derived by the Planck Collaboration [212]. We show two versions of the Planck constraints: the
74 3.7. Cosmological constraints
SDSS Stripe 82 weak lensing SDSS DR8 galaxy clustering Planck 2015 CMB, this work joint analysis
60 0.
.45 m 0 Ω 30 0.
15 0.
08 0.
.06
b 0 Ω 04 0.
02 0.
6 1.
2 1. s n 8 0.
4 0.
25 1.
8 .00
σ 1 75 0. 50 0.
.4 .6 .8 .0 15 30 45 60 02 04 06 08 .4 .8 .2 .6 50 75 00 25 0 0 0 1 0. 0. 0. 0. 0. 0. 0. 0. 0 0 1 1 0. 0. 1. 1. h Ωm Ωb ns σ8
Figure 3.6: Cosmological parameter constraints derived from the joint analysis, marginalized over τreion, b and m and from the single probes. The single probe constraints have been marginalized over the respective nuisance parameters i.e. τreion for the CMB temperature anisotropies, b for galaxy clustering and m for the weak lensing shear. In each case the inner (outer) contour shows the 68% c.l. (95% c.l.). For clarity the cosmic shear 68% c.l. are solid while the 95% c.l. are dashed.
75 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
joint analysis Planck Collaboration 2015, TT + lowP Planck Collaboration 2015, TT, TE, EE + lowP + lensing + ext
36 0.
32 0.
m 28 Ω 0.
24 0.
051 0.
.048 b 0 Ω 045 0. 042 0.
04 1.
00 1. s n 96 0.
92 0. 90 0.
85 0.
8 80 σ 0.
75 0.
70 0. 12 0. 09 0. reion τ 06 0. 03 0.
.60 .65 .70 .75 .80 .24 .28 .32 .36 042 045 048 051 .92 .96 .00 .04 .70 .75 .80 .85 .90 .03 .06 .09 .12 0 0 0 0 0 0 0 0 0 0. 0. 0. 0. 0 0 1 1 0 0 0 0 0 0 0 0 0 h Ωm Ωb ns σ8 τreion
Figure 3.7: Comparison between the parameter constraints derived from the joint analysis, marginalized over b and m and the constraints from Planck Collaboration et al. [212] using only CMB data (TT+lowP) or adding external data (TT,TE,EE+lowP+lensing+ext). The Planck constraints are marginalized over all nuisance parameters. In each case the inner (outer) contour shows the 68% c.l. (95% c.l.).
76 3.8. Conclusions constraints derived from the combination of CMB temperature anisotropies with the Planck low-` polarization likelihood (TT+lowP) and the ones derived from a combination of the latter with the Planck polarization power spectra, CMB lensing and external data sets (TT, TE, EE + lowP + lensing + BAO + JLA + H0). We see that the joint analysis prefers slightly lower values of the parameters Ωm and Ωb and a higher Hubble parameter h, but these differences are not significant. Despite this fact we find sensible overall agreement between the constraints derived in this work with both versions of the Planck constraints. While the constraints we derived in this analysis are broadened by the restricted multipole range we used, the results already demonstrate the power of integrated probe combination: the complementarity of different cosmological probes and their cross-correlations allows us to obtain reasonable constraints.
The measured power spectra together with the theoretical predictions for the best-fitting cos- mological model derived from the joint analysis are shown in Fig. 3.4. The best-fit cosmology γδg γγ provides a rather good fit to all power spectra except C` and C` , whose measured values are generally lower than our best-fit model. This is mainly due to the assumed Gaussian prior on the multiplicative shear bias m, which does not allow for more negative values of m as would be preferred by the data. If we relax the prior to a Gaussian with standard deviation σ 0.2, we find a best-fit value for the multiplicative bias parameter of m 0.276 0.108. = γδg γγ = − ± This results in an improved fit to both C` and C` , but is in tension with the values derived for the multiplicative bias by Hirata and Seljak [125]. We therefore find evidence for a slight tension between CMB temperature anisotropy data and weak gravitational lensing, as already seen by e.g. [177, 101].
3.8 Conclusions
To further constrain our cosmological model and gain more information about the dark sector, it will be essential to combine the constraining power of different cosmological probes. This work presents a first implementation of an integrated approach to combine cosmological probes into a common framework at the map level. In our first implementation we combine CMB temperature anisotropies, galaxy clustering and weak lensing shear. We use CMB data from Planck 2015 [209], photometric galaxy data from the SDSS DR8 [13] and weak lensing data from SDSS Stripe 82 [20]. We take into account both the information contained in the separate maps as well as the information contained in the cross-correlation between the maps by measuring their spherical harmonic power spectra. This leads to a power spectrum matrix with associated covariance, which combines CMB temperature anisotropies, galaxy clustering, cosmic shear, galaxy-galaxy lensing and the ISW [232] effect with galaxy and weak lensing shear tracers.
From the power spectrum matrix we derive constraints in the framework of a ΛCDM cosmo- logical model assuming both a Gaussian covariance as well as a Gaussian likelihood. We find that the constraints derived from the combination of all probes are significantly tightened
77 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing compared to the constraints derived from each of the three separate auto-power spectra. This is due to the complementary information carried by different cosmological probes. We further compare these constraints to existing ones derived by the Planck collaboration and find reasonable agreement, even though the joint analysis slightly prefers lower values of both
Ωm and Ωb and a higher Hubble parameter h. For a joint analysis of three cosmological probes, the constraints derived are still relatively weak, which is mainly due to our conservative cuts in angular scales. Nevertheless this analysis already demonstrates the potential of integrated probe combination: the complementarity of different data sets, that alone yield rather weak constraints on the full ΛCDM parameter space, allows us to obtain robust constraints which are significantly tighter than those obtained from probes taken individually. In addition, our analysis reveals challenges intrinsic to probe combination. Examples are the need for foreground-correction at the map as opposed to the power spectrum level and the need for coordinate-independent bias corrections.
In this first implementation we have made simplifying assumptions. We assume a Gaussian covariance matrix for all cosmological probes considered. This is justified for the CMB temper- ature anisotropies and the galaxy overdensity at large scales. The galaxy shears on the other hand exhibit non-linearities already at large scales and their covariance therefore receives significant non-Gaussian contributions [233]. Furthermore, we do not take into account the cosmology-dependence of the covariance matrix [75]. In addition we only include systematic uncertainties from a potential multiplicative bias in the weak lensing shear measurement and neglect effects from other sources. Finally we also used the Limber approximation for the theoretical predictions. We leave these extensions to future work but we do not expect them to have a significant impact on our results since we restrict the analysis to scales where the above effects are minimized.
In order to fully exploit the wealth of cosmological information contained in upcoming surveys, it will be essential to investigate ways in which to combine these experiments. It will be thus interesting to extend the framework presented here to include additional cosmological probes, 3-dimensional tomographic information and tests of cosmological models beyond ΛCDM.
78 3.A. Theoretical prediction for CMB and weak lensing shear cross-correlation
3.A Theoretical prediction for CMB and weak lensing shear cross- correlation
The CMB temperature anisotropies are correlated with the weak lensing shear due to the ISW effect. The anisotropies in the temperature field generated by time-varying gravitational potentials Φ are given by (see e.g. [197]):
Z η0 ∂Φ ∆TISW(θ) TCMB δTISW 2TCMB dη , (3.11) = = ηr ∂η where η0 denotes the conformal time today and ηr is the conformal time at recombination. Note that we follow the conventions for the gravitational potential Φ as in Bartelmann [25]. These anisotropies can be decomposed into spherical harmonics with multipole coefficients
Z η0 Z 3 ` d k d ∆TISW,`m 4πi 2TCMB dη 3 [Φ(k,z)]j`(kχ(z))Y`∗m(θk ). (3.12) = ηr (2π) dη
The multipole coefficients of the weak lensing shear E-modes can be expressed through the lensing potential ψ and are given by [25]
s s 1 (` 2)! (` 2)! ` aE,`m + ψ`,m + 4πi = −2 (` 2)! = − (` 2)! − − Z Z 3 d k0 dχg(χ) Φ(k0,z)j`(k0χ(z))Y ∗ (θk ), (3.13) × (2π)3 `m 0 where
Z χh 1 χ(z0) χ(z) g(χ) dz0 − n(z0). (3.14) = χ(z) χ(z) χ(z0)
γT The spherical harmonic power spectrum C` between CMB temperature anisotropies and the weak lensing shear is defined as
γT ∆TISW,`m a∗ C δ`` δmm . (3.15) 〈 E,`0m0 〉 = ` 0 0 Expressing the integrals in terms of redshift and interchanging the integration boundaries gives
s (` 2)! Z z Z d3k d 2 + ∗ ∆TISW,`m aE,∗ ` m (4π) 2TCMB dz 3 [D(z)(1 z)]Φ(k,z 0) 〈 0 0 〉 = (` 2)! 〈 0 (2π) dz + = − Z Z 3 c d k0 j (kχ(z))Y (θ ) dz g ¡χ(z )¢ Φ(k ,z )j (k χ(z ))Y (θ ) , (3.16) ` `∗m k 0 0 3 0 0 `0 0 0 `0m0 k0 × H(z0) (2π) 〉 where z? denotes the redshift at recombination. In order to derive Eq. 3.16 we have used that in linear perturbation theory the time- and scale-dependence of the gravitational potentials
79 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing can be separated i.e.:
Φ(k,z) Φ(k,z 0)D(z)(1 z), (3.17) = = + where D(z) denotes the linear growth factor. We further have that
3 lin Φ(k,z 0)Φ(k0,z0 0) (2π) P (k,z 0)δ(k k0), (3.18) 〈 = = 〉 = ΦΦ = − and therefore Eq. 3.16 reduces to
s (` 2)! Z z Z k2dk d γT 2 + ∗ C` (4π) 2TCMB dz 3 [D(z)(1 z)] = (` 2)! 0 (2π) dz + Z − c ¡ ¢ lin dz0 g χ(z0) D(z0)(1 z0)PΦΦ(k,z 0) j`(kχ(z)) j`(kχ(z0)). (3.19) × H(z0) + = Eq. 3.19 is the exact expression for the spherical harmonic cross-power spectrum between CMB temperature anisotropies and weak lensing shear. In order to speed up computations, it can be simplified by resorting to the Limber approximation [169, 147, 148] which gives
s (` 2)! Z z d γT + ∗ C` 2TCMB dz [D(z)(1 z)] = (` 2)! 0 dz + − ¡ ¢ µ 1 ¶ g χ(z) ` /2 D(z)(1 z) P lin k + ,z 0 . (3.20) × + χ2(z) ΦΦ = χ(z) =
The power spectrum of the gravitational potential at late times is related to the matter power spectrum through Poisson’s equation
µ ¶2 2 4 lin 3 ΩmH P (k,z 0) P lin (k,z 0) 0 δδ = . (3.21) ΦΦ = = 2 c4 k4
For large ` we can make the approximations
s (` 2)! + `2, (` 2)! ∼ (3.22) − 1 2 2 (` /2) ` . + ∼ Using Equations 3.21 and 3.22 we can write Eq. 3.20 as
2 Z µ 1 ¶ γT ΩmH TCMB 1 d ` /2 C 3 0 dz [D(z)(1 z)]D(z)W γ ¡χ(z)¢P lin k + ,0 , (3.23) ` 2 1 2 δδ = c (` /2) dz + = χ(z) + which is the expression given in Eq. 3.5.
80 3.B. Treatment of systematic uncertainties in galaxy clustering data
3.B Treatment of systematic uncertainties in galaxy clustering data
The number density of galaxies observed in SDSS DR8 photometric data is affected by various systematic uncertainties such as stellar density, Galactic extinction and PSF size variation [229, 129]. These effects remain even after masking and removal of the highest contamination regions. In order to obtain an unbiased galaxy overdensity map, we need to correct for the number density variation due to systematics. The SDSS recorded the values of several potential systematic uncertainties for the observed fields: airmass, Galactic extinction and seeing (as measured by the PSF FWHM) in all 5 SDSS bands for the field each galaxy has been observed in as well as sky emission at the position of the galaxy for all the 5 bands. These quantities can be queried for each galaxy position on the CAS17. In this work, we consider four different observational systematics: Galactic extinction in the r-band as well as FWHM of the PSF, airmass and sky emission in the i-band. A further potential systematic uncertainty is the presence of foreground stars. Ross et al. [229] show that the effects of foreground stars on galaxy number density are largely independent of the magnitude of the stars. We therefore follow Ho et al. [129] and investigate how the number density of stars with i-band magnitudes in the range 18.0 i 18.5 affects the number density of detected galaxies. ≤ < We pixelize all quantities onto HEALPix maps of resolution NSIDE 1024 and compute the = number density of galaxies relative to their mean number density as a function of the value of the systematic in the pixel. In order to correct for these systematic uncertainties, we fit a 3rd-order polynomial to the functional dependence of the relative galaxy number density on the systematic. Then we multiply the uncorrected number densities by the inverse of this function. Various potential systematics such as Galactic extinction and stellar density are spatially correlated to one another. When correcting for various systematics simultaneously, the order in which the corrections are applied could influence results [229]. In our sample we find that the corrections are both independent of ordering and SDSS band and correcting for the effect in one band simultaneously corrects for all the other bands. The results are shown in Fig. 3.8 and we use those to correct the galaxy maps from residual systematic uncertainties. We clip the systematics maps at the minimum and maximum systematics value shown in the figure and apply the fitted corrections to the galaxy number density. The galaxy clustering spherical harmonic power spectra before and after correcting for systematic uncertainties are discussed in Section 3.5.2 and shown in Appendix 3.H.
3.C PSF correction and construction of weak lensing shear maps
The galaxy shapes measured from images represent a convolution of the intrinsic galaxy shapes with the PSF of the telescope and the atmosphere. We therefore need to correct this effect using PSF estimates measured from the shapes of stars observed in the survey. As described in Annis et al. [20] and Lin et al. [170], the PSF model for SDSS Stripe 82 data is derived from weighted sums of shapes measured in the individual runs as opposed to the co-adds. Lin et al.
17See footnote9.
81 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
1.03 1.02 1.010 1.02 1.01 1.005 1.01 1.00 1.00 0.99 1.000 gal 0.99 gal gal ¯ ¯ ¯ n n n
/ / 0.98 /
gal 0.98 gal gal n n 0.97 n 0.995 0.97 0.96 0.96 0.990 0.95 0.95 0.94 0.94 0.985 0 2 4 6 8 10 0.8 1.0 1.2 1.4 1.6 1.8 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1 stars [pixel− ] psffwhm i [arcsec] airmass i 1.02 1.06
1.01 1.04 1.02 1.00 fit 1.00 uncorr. gal 0.99 gal ¯ ¯ n n
/ / 0.98 corr.
gal 0.98 gal n n 0.96 0.97 0.94 0.96 0.92 0.95 0.90 6 8 10 12 14 0.01 0.02 0.03 0.04 0.05 0.06 0.07 sky i [nanomaggies/arcsec2] extinction r [mag]
Figure 3.8: Galaxy number density ngal relative to mean galaxy density n¯gal as a function of potential systematic value. The figures show both the uncorrected data and the data corrected using a third order polynomial fit to the uncorrected relation. The error bars assume Poisson noise and are thus likely underestimated due to the correlations between galaxy positions.
[170] found that this leads to biases that need to be removed prior to PSF correction. In order to correct for these effects, we follow the steps outlined in Lin et al. [170]. We select bright stars with i-band magnitudes in 16 i 17, which pass the clean photometry cuts18, and fit < < polynomials to the residuals between their shapes measured from the co-adds and the PSF model for these stars. The residuals before and after subtraction of the polynomial fit are shown in Fig. 3.9. We see that the correction introduced in Lin et al. [170] has considerably removed both an overall bias as well as discontinuities at the CCD camera column (camcol) edges.
Using the revised PSF model, we correct the measured shapes for the effect of the PSF.We use the linear PSF correction algorithm derived in Hirata and Seljak [125], which can be applied to the adaptive moment measurements from the SDSS PHOTO pipeline19.
In order to obtain a galaxy sample for reliable weak lensing shear measurement, we follow Lin et al. [170] and perform two additional selection cuts on the galaxies after PSF correction: we select galaxies with ellipticity components e ,e satisfying e 1.4 as well as e 1.4 and 1 2 | 1| < | 2| < photometric redshift uncertainties σ 0.15. This additional selection leaves a galaxy sample z < consisting of N 3322915 galaxies. gal = Lin et al. [170] found a camcol dependent additive bias in the PSF-corrected ellipticities. The 5 mean ellipticities for each camcol lie in the range e¯ [6 10− ,0.02] and e¯ [0.002,0.009], | 1| = × | 2| = which is larger than expected for a mean zero field [170]. We therefore follow Lin et al. [170] and correct for the additive bias by subtracting the mean ellipticity for each camcol. We choose
18See footnote14. 19 Note that there is a typo in Hirata and Seljak [125]: The quantities Cg ,C f ,Dg ,D f in Eq. (B9) should be squared.
82 3.C. PSF correction and construction of weak lensing shear maps
0.10 0.10
0.05 0.05
mE1PSF 0.00 mE1PSF 0.00 − − 1 1
mE 0.05 mE 0.05 − −
0.10 0.10 − 1.5 1.0 0.5 0.0 0.5 1.0 1.5 − 1.5 1.0 0.5 0.0 0.5 1.0 1.5 − − − − − − DEC [deg] DEC [deg]
Figure 3.9: Residuals between ellipticity component e1 of a random sample of 10000 bright stars measured on the co-add images and PSF models for these objects as a function of declination (DEC) before (left panel) and after (right panel) applying the correction described in the text.
to perform this step prior to coordinate transformation (i.e. for ellipticities defined relative to equatorial coordinates) as opposed to after rotation. We find that removing the mean camcol ellipticity reduces PSF leakage to a level which is subdominant in our analysis.
Fig. 3.10 shows the distributions of the ellipticity components e1 and e2 defined relative to equatorial coordinates. They are averaged over HEALPix pixels of resolution NSIDE 512, = which corresponds to a pixel area of A 0.013 deg2. The figure displays the ellipticity pix ≈ histograms both prior to PSF correction and subtraction of additive bias as well as the final dis- tributions obtained after applying both corrections. We see that the corrections have removed the effects of the PSF and the final histograms can be described by Gaussian distributions.
In the final step, these ellipticities need to be transformed to shear estimates by correcting for the shear resolution factor R, which is defined as ¿ À ∂γˆi R . (3.24) = ∂γi
The shear resolution factor R quantifies the response of the estimated mean ellipticity to an applied shear. For the adaptive moment method described in Bernstein and Jarvis [32] it is given by R 2(1 e2 ), where e denotes the intrinsic r.m.s. ellipticity per component. We = − int int follow Lin et al. [170] and use e 0.37 as measured by Hirata et al. [127]. int = In order to construct the final weak lensing shear maps we thus apply the resolution correction to the ellipticities and transform them from equatorial to Galactic coordinates.
83 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
103 103 before PSF after PSF corr. corr. 102 102 ) ) 1 1 e e ( ( n n 101 101
100 100 0.25 0.20 0.15 0.10 0.05 0.00 0.05 0.10 0.15 0.20 0.20 0.15 0.10 0.05 0.00 0.05 0.10 0.15 0.20 − − − − − − − − − e1 e1 103 103 before PSF after PSF corr. corr. 102 102 ) ) 2 2 e e ( ( n n 101 101
100 100 0.20 0.15 0.10 0.05 0.00 0.05 0.10 0.15 0.20 0.20 0.15 0.10 0.05 0.00 0.05 0.10 0.15 0.20 − − − − − − − − e2 e2
Figure 3.10: Histograms of the ellipticity components e1 and e2 averaged over HEALPix pixels of resolution NSIDE 512. The left panels show the distributions before correction for the = PSF and subtraction of a camcol dependent additive bias, while the right panels show the distributions after the application of these corrections.
84 3.D. Transformation of weak lensing shear under rotation
3.D Transformation of weak lensing shear under rotation
We rotate the weak lensing galaxy shears from equatorial to Galactic coordinates following the implementation in HEALPix. The method is briefly summarized below.
The rotation angle of the shears under a coordinate rotation as described by the rotation matrix R is equal to twice the rotation angle ψ of the coordinate axes with respect to which they are defined. In HEALPix the x-axis is in the direction of eφ and the y-axis in the direction of eθ.
In order to derive ψ we define the following quantities: The position before rotation is denoted ³ ´ ³ ´ as r x y z and the position after rotation is r0 x y z . We further define the = = 0 0 0 vector towards the north pole in the unrotated coordinate system, which is given by p ³ ´ ³ ´ 1 = x0 y0 z0 0 0 1 . Under the inverse rotation R− the north polar vector is mapped to ³ = ´ p00 x y z . At the position r the unit vectors in θ- and φ-direction are given by = 000 000 000 p r eφ × , = p r | × | (3.25) (p r) r e × × . θ = (p r) r | × × | We have the following identities
1 R eφ\θ p eφ\θ R− p, · = · (3.26) 1 R− p r p R r. · = · Taking into account the left-handedness of the HEALPix coordinate system and inserting the explicit expressions, it follows that
c cosψ (zz0 z000), = p1 z2 − c− (3.27) sinψ (x y000 yx000), = p1 z2 − − where c is a constant, which we can remove by ensuring that sin2 ψ cos2 ψ 1. Under this + = rotation the weak lensing shear transforms as
γ10 cos2ψγ1 sin2ψγ2, = + (3.28) γ0 sin2ψγ cos2ψγ . 2 = − 1 + 2
3.E Choice of PolSpice parameter settings
In this section we illustrate the determination of the maximal angular scale θmax used to compute spherical harmonic power spectra on the example of the SDSS Stripe 82 mask. Fig. 3.11 shows the real space correlation function of this mask. It is non-zero for small angular
85 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
0.007 0.006 0.005 0.004 ) θ
( 0.003 C 0.002 0.001 0.000 0.001 − 0 2 4 6 8 10 12 14 θ [deg]
Figure 3.11: Real space correlation function of the SDSS Stripe 82 mask. The dashed line denotes the value chosen for θmax.
scales, then starts to fall-off and approximately vanishes for large angular separations. From this figure we see that θmax is not a well-defined quantity. Our approach is thus to choose a maximal angular scale by eye and validate it on the Gaussian simulations. For the SDSS Stripe 82 mask, we choose θ 10 degrees and θ 5 degrees. When testing these PolSpice max = FWHM = settings on the simulations, we find a reasonable agreement between input and recovered power spectra.
3.F Correlated maps of spin-0 and spin-2 fields
Our analysis relies on Gaussian simulations both for validation of the data analysis pipeline and covariance matrix estimation. We thus need to generate correlated HEALPix maps of both spin-0 and spin-2 fields from input auto- and cross-power spectra. Cabré et al. [42] and Giannantonio et al. [93] describe an algorithm for generating correlated HEALPix maps of spin-0 fields. In order to consistently simulate the weak lensing shear field, we extend this algorithm to also include correlations between spin-0 and spin-2 fields.
These algorithms are all based on the HEALPix routine synfast, which generates HEALPix ii maps of realizations of input spherical harmonic power spectra C` . If the fields are addition- ally mean-subtracted, this is equivalent to requiring that the spherical harmonic coefficients a`m of the maps satisfy
D i E a`m 0, = (3.29) D i i E ii a a ∗ C δ`` δmm . `m `0m0 = ` 0 0 In synfast these conditions are imposed by assigning a random phase ξ with mean 0, ξ 0, 〈 〉 = 86 3.F. Correlated maps of spin-0 and spin-2 fields
and unit variance, ξξ∗ 1, to each spherical harmonic mode ` and setting 〈 〉 = q ai C ii ξ. (3.30) `m = ` As derived in Giannantonio et al. [93], this method can be extended to correlated maps using more random phases. The simplest case is to create two correlated spin-0 zero maps with ii j j i j power spectra C` , C` and cross-power spectrum C` . This is the only case relevant for our work and it is achieved by choosing the amplitudes of the maps of the two probes i, j as [93]
q ai C ii ξ , `m = ` 1 v i j u i j 2 (3.31) j C u j j (C ) a ` tC ` . `m q ξ1 ` ξ2 = ii + − C ii C` `
As described in [93] this algorithm can be implemented using synfast by first creating a ii map with power spectrum C` and a second map using the same seed with power spectrum i j (C )2 ii ` /C` . Finally the second map needs to be added to a third map, generated with a different
j j i j (C )2 ii random seed and with power spectrum C ` /C . This ensures the desired auto- and ` − ` cross-correlations.
To extend this algorithm to spin-2 fields, we make use of the polarization version of synfast, which allows us to generate correlated spin-0 and spin-2 maps consistent with input auto- 00 and cross-power spectra. Let 0 denote the spin-0 field. Then C` denotes the auto-power EE BB spectrum of the spin-0 field, C` , C` are the E- and B-mode power spectra of the spin-2 field 0E and C` is the cross-power spectrum between the spin-0 field and the spin-2 E mode. Given these input power spectra, the polarization mode of synfast generates a map of the spin-0 field and two maps of the spin-2 field with the desired auto- and cross-power spectra.
In order to obtain correlated maps mT,mδg ,mγ1 ,mγ2 of CMB temperature anisotropies, galaxy TT δg T overdensity and galaxy weak lensing shear with auto- and cross-power spectra C` ,C` , δg δg γT γδg γγ C` ,C` ,C` ,C` we therefore proceed as follows:
(i) We first create three correlated HEALPix maps using synfast in polarization mode with the power spectra
C 00 C TT, ` = ` EE C γγ C ` /2, ` = C BB 0, ` = T C 0E C γ . ` = `
These maps are denoted m1, where i {T,γ ,γ }. i ∈ 1 2 87 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
(ii) Following Eq. 3.31, we then create three maps with a new random seed and the power spectra
δ δ 00 g g δg T C C (C )2/C TT, ` = ` − ` ` EE C γγ C ` /2, ` = C BB 0, ` = γδ C 0E C g . ` = `
These maps are denoted m2, where i {δ ,γ ,γ }. i ∈ g 1 2 (iii) We create a spin-0 map generated with the same seed as used for m1 with the power spectrum
00 δg T C (C )2/C TT, ` = ` ` which is called m3.
(iv) Finally we combine the maps i.e.
m m1, T = T 2 3 mδ m m , g = δg + m m1 m2 , γ1 = γ1 + γ1 m m1 m2 . γ2 = γ2 + γ2
TT δg T This procedure yields four correlated maps with auto- and cross-power spectra C` ,C` , δg δg γT γδg γγ C` ,C` ,C` ,C` . The algorithm described above introduces an unwanted, additional correlation between mδg and mγ1 ,mγ2 . It can in principle be corrected for by adding counter- terms to the respective maps. Since the additional correlation is subdominant in the present case, we neglect these counter-terms.
In order to obtain realistic maps we need to account for the effects of HEALPix pixel and beam window function. The signal measured in each HEALPix pixel is a convolution of the underlying signal with the HEALPix window function. If further experimental beams are present, the signal is additionally convolved with the beam window function. Since a convolution in real space is equivalent to a multiplication in Fourier space, we account for these effects by multiplying the input power spectra by the power spectra of the respective window functions prior to generating the HEALPix maps.
To compute the covariance matrix as well as to validate the analysis pipeline we need to add realistic noise to the correlated Gaussian simulations. We choose to add the noise on the map level. For the CMB temperature anisotropies we add the Commander HMHD map provided by the Planck collaboration to each simulated temperature map. We do not randomize the
88 3.G. Validation of spherical harmonic power spectrum measurements noise map for each new realization since the HMHD map features significant correlations which would be lost by randomizing. Since we are adding the same noise map to each random realization we expect to slightly underestimate the noise using our simulations. However, we do not expect this to have a significant effect on our results, since the noise in the CMB temperature power spectrum is dominated by cosmic variance at the scales considered. For the galaxy overdensity field we create noise maps by randomizing the positions of the galaxies in our data inside the survey mask. We then pixelize those on a HEALPix map and add the noise map to the simulated map. The galaxy shear noise maps are created by rotating each galaxy shear by a random angle and repixelizing the rotated shears onto HEALPix maps. As before these noise maps are added to the signal maps to produce the Gaussian simulations including both signal and noise.
3.G Validation of spherical harmonic power spectrum measurements
We validate the spherical harmonic power spectrum measurement outlined in Section 3.5 using the correlated Gaussian simulations described in Appendix 3.F. We compute theo- retical predictions for the six spherical harmonic power spectra considered in this work, i.e. TT δg T δg δg γT γδg γγ C` , C` , C` , C` , C` , C` for a ΛCDM cosmological model with parameters {h, Ωm, Ωb, Ω , n , σ , τ , T } {0.7, 0.3, 0.049, 1.0, 0.88, 0.078, 2.275K}. We further set the linear, Λ s 8 reion CMB = redshift-independent galaxy bias parameter to b 2. = Using the algorithm described above, we generate N 1000 Gaussian realizations of these sim = power spectra and add the noise maps determined from the data. We then apply angular masks equivalent to those in the data and compute the spherical harmonic power spectra from the masked maps using the exact same methodology and PolSpice settings as applied on the data. To estimate the noise bias we follow the same randomization approaches as outlined in Section 3.5.
Figures 3.12 and 3.13 show a comparison between the input power spectra for all the six probes and the means of the recovered realizations. The uncertainties are derived from the sample covariance of the Gaussian realizations. We see that the input power spectra are recovered rather well. Also shown are the χ2 values between the reconstruction and the theory. These are not rigorous measures for the goodness of the recovery since they strongly depend on binning and angular multipole range chosen.
3.H Spherical harmonic power spectrum robustness tests
This section summarizes the robustness tests performed for the spherical harmonic power spectra.
89 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
6000 χ2/dof = 1.82
) 5000 π (2
/ 4000 TT ` theor.
C 3000 reconst.
+ 1) 2000 ` ( ` 1000 0 0 100 200 300 400 500 600 9 10− 4 2.0 × 10− χ2/dof = 1.15 χ2/dof = 0.64 1.5
1.0 g T δ g g δ δ ` ` C 0.5 C
5 0.0 10−
0.5 − 50 100 150 200 3 101 102 × 11 4 10− 10− 12 χ2/dof = 1.33 χ2/dof = 0.23 χ2/dof = 2.91 10−
13 10− 7
g 10− T γγ γ 14 ` 5 γδ ` 10− ` 10− C C C 15 10−
16 10−
17 8 6 10− 10− 10− 102 6 102 3 101 102 102 6 102 × × × ` ` `
Figure 3.12: Comparison between input power spectra and mean recovered power spectra as estimated from N 1000 Gaussian realizations generated using the algorithm outlined in sim = Appendix 3.F. The noise level of the Gaussian realizations is tuned to the data and the spherical harmonic power spectra have been computed using the same methodology and PolSpice settings as applied on the data. The angular multipole ranges and binning schemes for all power spectra are summarized in Table 3.2. Dashed lines denote negative spherical harmonic power spectrum values.
6 10− 0.4 × 0.2 0.0 0.2 BB ` − C 0.4 − 0.6 reconst. − 2 0.8 χ /dof = 1.33 − 1.0 − 100 200 300 400 500 600 `
Figure 3.13: The cosmic shear B-mode power spectrum reconstructed from N 1000 sim = Gaussian realizations generated using the algorithm outlined in Appendix 3.F. The angular multipole range and binning scheme is summarized in Table 3.2.
90 3.H. Spherical harmonic power spectrum robustness tests
4 10− δgδg g δ
g gal. δ `
C eq.
5 10− 3 101 102 ×5 10− 1.5 × γδg γγ
1.0 ) 5 π 10− (2 /
g 0.5 γγ ` γδ ` C
`C 0.0 6 10− + 1) ` ( 0.5 ` − 7 1.0 10− − 50 100 150 200 102 6 102 × ` `
Figure 3.14: Comparison between spherical harmonic power spectra computed from the maps in Galactic and equatorial coordinates.
3.H.1 Comparison between spherical harmonic power spectra in equatorial and Galactic coordinates
We test that the spherical harmonic power spectra involving maps which can be transformed between coordinate systems, i.e. galaxy overdensity and weak lensing shear maps, are unaf- fected by the rotation. The comparison between spherical harmonic power spectra computed from maps in Galactic and equatorial coordinates are shown in Fig. 3.14. We find good agree- δg δg γδg ment between the two power spectra for both C` and C` , while we find discrepancies for γγ C` . We attribute this to the additive bias correction applied to the galaxy shears as outlined in Section 3.5.3. The additive bias correction described in Appendix 3.C, causes an asymmetry between the galaxy shears in different coordinate systems, which is the cause for the large discrepancies detected. This can be seen from Fig. 3.15, which shows a comparison between the cosmic shear power spectra prior to noise removal as estimated from maps in Galactic and equatorial coordinates. The left panel shows the comparison when the additive bias correction is applied while in the right panel we do not apply any correction. As can be seen, we find discrepancies when we apply the additive bias correction in equatorial coordinates and then rotate the corrected shears to Galactic coordinates. Not applying any additive bias correction on the other hand, removes most of these effects.
91 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
9 9 10− 10− × × gal. gal. 5 eq. eq. 4 γγ ` 3 C 2 1 0 100 200 300 400 500 600 100 200 300 400 500 600 ` `
Figure 3.15: Comparison between cosmic shear spherical harmonic power spectra prior to noise removal measured from the maps in Galactic and equatorial coordinates. The left hand panel shows the results when applying the correction for additive bias in equatorial coordinates and then rotating the shears to Galactic coordinates. The right hand panel shows the results when no additive bias correction is applied.
3.H.2 Comparison between spherical harmonic power spectra derived from dif- ferent foreground-reduced CMB maps
We test that the spherical harmonic power spectra involving CMB data are unaffected by our choice of foreground-reduced map. The power spectra involving CMB data are shown in Fig. 3.16 for the foreground-reduced CMB maps derived using the component separation methods Commander, NILC, SEVEM and SMICA. As can be seen, the power spectra are virtually the same.
3.H.3 Impact of systematics correction on galaxy clustering power spectrum
We further investigate the effect of systematics correction on the galaxy clustering power spectrum. The galaxy clustering spherical harmonic power spectra before and after correcting for systematic uncertainties are shown in Fig. 3.17. Our systematics removal method slightly reduces the clustering amplitude at large angular scales, while leaving small angular scales al- most unaffected. This is to be expected since Galactic foregrounds typically exhibit significant large scale clustering.
92 3.H. Spherical harmonic power spectrum robustness tests
7 10− 6000 × ] 2 TT 0.8 δgT K 5000 µ
)[ 0.6
π 4000 [K] (2 0.4 / T
3000 g δ ` TT ` 0.2 C `C 2000 0.0 + 1) `
( 1000 0.2 ` − 0 0.4 0 100 200 300 400 500 600 − 50 100 150 200 9 10− ` 2.5 × 2.0 γT
1.5 com. smica
[K] 1.0 sevem nilc T γ ` 0.5 `C 0.0 0.5 − 1.0 − 100 200 300 400 500 600 `
TT δg T γT Figure 3.16: Comparison between spherical harmonic power spectra C` , C` , C` derived using the four different foreground-reduced CMB maps from Commander, NILC, SEVEM and SMICA.
4 10− gal. corr. gal. uncorr. g δ g δ ` C
5 10− 3 101 102 × `
Figure 3.17: Comparison between galaxy overdensity power spectra computed before and after systematics removal.
93 Chapter 3. Integrated approach to cosmology: combining CMB, large-scale structure and weak lensing
60◦
30◦
120◦ 60◦ 0◦ 60◦ 120◦ 0◦ − −
30◦ −
60◦ −
7.5 7.5 − 5 ∆T [K] 10− ×
Figure 3.18: Half-mission half-difference (HMHD) Commander CMB temperature anisotropy map. This map contains only noise and potential residual systematics.
5 10− 1.5 × 1.0 ) π
(2 0.5 / BB `
C 0.0
+ 1) 0.5 `
( − ` 1.0 − 1.5 − 100 200 300 400 500 600 `
Figure 3.19: Spherical harmonic power spectrum of cosmic shear B-modes computed from the SDSS Stripe 82 maps in Galactic coordinates. The angular multipole range and binning scheme is summarized in Table 3.2.
94 4 Integrated cosmological probes: ex- tended analysis
E quindi uscimmo a riveder le stelle. — Dante, Inferno, Canto XXXIV
This chapter appeared in a similar form as Nicola, Réfrégier & Amara, 2017 [193].
4.1 Introduction
In Chapter3, we implemented an integrated approach to cosmology in which the cosmological probes are combined into a common framework at the map level. This has the advantage of taking full account of the correlations between the different probes which generally probe common survey volumes, to provide a stringent test of systematics through the test of the consistency between the probes and to yield a test of the validity of the cosmological model. We applied this framework to a combination of the CMB temperature from the Planck mission [209], galaxy clustering from the eighth data release of the Sloan Digital Sky Survey (SDSS DR8) [13] and weak lensing from SDSS Stripe 82 [20], making simplifying approximations but also conservative cuts on the data.
In the present work, we extend the integrated analysis of Chapter3 to also include CMB lensing maps from the Planck mission [211], the recent weak lensing measurement with the publicly available Dark Energy Survey (DES) Science Verification (SV) data [138], SNe Ia data from the joint light curve analysis (JLA) [34] and constraints on the Hubble parameter from the Hubble Space Telescope (HST) [227, 74]. This yields 12 auto and cross power spectra which include the CMB temperature power spectrum, cosmic shear, galaxy clustering, galaxy- galaxy lensing CMB lensing cross-correlation along with other cross-correlations as well as background probes. Furthermore, we extend the treatment of systematic uncertainties and relax some of the approximations as compared to Chapter3. In particular, we study the impact of intrinsic alignments, baryonic corrections, residual foregrounds in the CMB temperature,
95 Chapter 4. Integrated cosmological probes: extended analysis and calibration factors for the different power spectra. This extended analysis allows us to derive more robust constraints on the ΛCDM cosmological model and a more thorough test of the consistency between the different probes. Other joint analyses of different sets of cosmological probes have been performed (see references in Chapter3 and Refs. [ 94, 92, 253]).
This chapter is organized as follows. We review the framework for cosmological probe com- bination employed in this work in Section 4.2. In Section 4.3, we describe the data used in this work and Section 4.4 describes the theoretical modeling of the cosmological observables. We detail the computation of spherical harmonic power spectra in Section 4.5, while Section 4.6 summarizes the systematic uncertainties considered in this work. The computation of the covariance matrix is described in Section 4.7. The method for parameter inference is described in Sec. 4.8 and our results on cosmological constraints are presented in Section 4.9. We conclude in Section 4.10. Robustness tests as well as implementation details are deferred to the Appendices.
4.2 Framework for integrated probe combination
The integrated probe combination framework is described in detail in Chapter3, but here we summarize the most important steps: first, we create projected two-dimensional maps for the large-scale structure (LSS) and CMB probes. We then compute both the spherical harmonic auto-power spectra of these probes as well as the cross-power spectra for physically overlapping surveys. This yields a set of 12 spherical harmonic power spectra, which does not include the auto-power spectrum of the CMB lensing convergence but only its cross- correlations. We complement the observed power spectra with theoretical predictions and an estimate of their covariance matrix and combine these into a Gaussian likelihood. We combine the power spectrum likelihood with the likelihood of SNe Ia distance moduli and a constraint on the Hubble parameter, assuming these probes to be independent. In a last step we compute cosmological parameter constraints in a joint fit to these data. The implementation details for the CMB lensing convergence, weak lensing data from DES SV, SNe Ia and the Hubble constant measurement are described below. For a description of the remaining data the reader is referred to Chapter3.
4.3 Data
The data used in this analysis is summarized in Table 4.1 and the footprints of the different surveys are illustrated in Figure 4.1 together with the background probes. We consider the data used in Chapter3 namely, the Planck 2015 foreground-reduced CMB temperature anisotropy map derived using the Commander algorithm [210], a map of the galaxy overdensity field derived using the CMASS1-4 sample from SDSS DR8 [13, 129, 229] and a weak lensing shear map derived using SDSS Stripe 82 co-add data [20, 170]. In addition to the cosmological maps we also employ the binary survey masks presented in Chapter3. In this work, we complement
96 4.3. Data
46 44 42
) 40 z (
µ 38 36 34 32 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 z 0.14 0.12 0.10 )
0 0.08 H (
p 0.06 0.04 0.02 0.00 60 65 70 75 80 H0
Planck temp. Planck lensing SDSS DR8 SDSS Stripe 82 DES SV
Figure 4.1: Summary of the data used in this work. The left hand side shows an overlay of the footprints of the surveys used in this work: the CMB temperature and CMB lensing convergence from Planck, the galaxy density from SDSS DR8 and weak lensing from SDSS Stripe 82 and DES SV. The right hand side shows the background probes: SNe Ia from JLA and Hubble parameter data from HST. All footprints are shown in Galactic coordinates. The map is shown in Mollweide projection at a HEALPix resolution of NSIDE = 1024. (See Tab. 4.1 for references for these different surveys.)
these three maps with several data sets as described below.
4.3.1 DES weak lensing
We use publicly available data from the DES SV.1 The DES is an ongoing survey, imaging the sky in five photometric bandpasses (g,r,i,z,Y ) using DECam [83]. After its five-year duration the DES will have covered approximately 5000 deg2 of the southern sky to a limiting magnitude of about 24. The SV data were taken before the start of the main survey and they consist of more than 250 deg2 [63]. In our analysis we use the largest contiguous area in the DES SV data, which is part of the South Pole Telescope East (SPT-E) field and covers an area of approximately 139 deg2. The weak lensing shear for galaxies in the SPT-E region has been measured using two independent shape measurement codes, NGMIX [245] and IM3SHAPE [285]. Both are model-fitting shear measurement codes and are described in Ref. [138]. Photometric redshifts (photo-z) have been obtained using four different methods as described in Ref. [37]. The photometric redshift catalogs both provide the full photo-z probability distribution function (pdf) as well as an estimate of the mean of the pdf for each galaxy. We follow the choice of
1The data is available at: https : //des.ncsa.illinois.edu/releases/sva1.
97 Chapter 4. Integrated cosmological probes: extended analysis
Table 4.1: Summary of the data sets used in our analysis.
Survey: Planck 2015 [210] CMB temperature Fiducial foreground-reduced map: Commander Sky coverage: f 0.776 sky = Survey: SDSS DR8 [13] Sky coverage: f 0.27 sky = Galaxy density Galaxy sample: CMASS1-4 Number of galaxies: N 854063 gal = Chapter3 Photometric redshift range 0.45 z 0.65 ≤ phot < Survey: SDSS Stripe 82 co-add [20] Sky coverage: f 0.0069 SDSS sky = Number of galaxies: N 3322915 Stripe 82 gal = Photometric redshift range: 0.1 . zphot . 1.1 r.m.s. ellipticity per component: σ 0.43 Weak e ∼ lensing Survey: DES SV [138] Sky coverage: f 0.0039 sky = Number of galaxies: N 3279967 DES gal = Sec. 4.3.1 Photometric redshift range: 0.3 z 1.3 < phot < r.m.s. weighted ellipticity per component: σ 0.24 e ∼ Survey: Planck 2015 [211] CMB lensing Sec. 4.3.2 Sky coverage: f 0.67 sky = Compilation: JLA [34] SNe Type Ia Number of SNe: N 740 Sec. 4.3.3 SNe = Redshift range: 0.01 z 1.3 < < Distance anchor: NGC 4258 [135] Number of Cepheids: N 600 [227] Hubble parameter Ceph. = Sec. 4.3.4 Number of SNe: N 8 [227] SNe = Analysis: Efstathiou [74]
98 4.3. Data
fiducial catalog of Refs. [29, 8] and perform our analysis using the galaxy shapes measured by NGMIX and the photometric redshifts determined using SkyNet.
Our analysis closely follows the spherical harmonic power spectrum measurement described in Appendix A of Ref. [29]. We select objects passing the SVA1 and the NGMIX cuts defined in Ref. [138] which fall into any of the three tomographic redshift bins described in Refs. [37, 138]. This selection yields N 3279967 galaxies. gal = In order to construct the weak lensing shear maps we weight each galaxy’s shear by its inverse variance weight described in Ref. [138]. The galaxy shapes given in the DES SV shear catalogues are biased estimators of the galaxy shears and the NGMIX shape estimates therefore need to be corrected for the sensitivity as described in Ref. [138]. Since this correction factor is a noisy estimate of the true correction it cannot be applied on single galaxies. In order to avoid introducing a bias caused by the noisy estimators of the sensitivity we therefore follow Ref. [29] and estimate the weighted average of the galaxy sensitivities in our sample and correct each galaxy shape with this mean correction.
We then rotate the galaxy shears from equatorial to Galactic coordinates2 and pixelize them onto HEALPix3 [99] pixelizations of the sphere choosing a resolution of NSIDE = 1024. This resolution corresponds to a pixel area of 11.8 arcmin2. We apply a binary mask constructed from the union of unobserved and empty pixels to both shear maps. The final maps cover a fraction of sky f 0.0039. The mean number of galaxies per pixel is approximately given by sky = n 2 gal/pix 67, which corresponds to n 5.73 arcmin− . Figure 4.2 shows the map of the shear = gal = modulus together with a zoom-in region with overlaid whisker plot illustrating the direction of the shear.
We follow Refs. [37, 29] and estimate the redshift distribution of the galaxies from the sum of the individual galaxy pdfs, weighted by their weak lensing shear weights. The resulting redshift distribution together with the weak lensing window function is shown in the Appendix (Fig. 4.19).
4.3.2 CMB lensing convergence
CMB lensing causes statistical anisotropies in CMB maps and the lensing potential can be reconstructed from these maps using a quadratic estimator [195]. We use the CMB lensing potential estimate φˆCMB provided by the Planck Collaboration in their second data release [211]. This estimator has been derived from the foreground-reduced CMB temperature and polarization maps computed using the SMICA algorithm [211, 210]. The use of both CMB temperature as well as polarization data allows for several CMB lensing potential estimators
(φˆTT,φˆTE,φˆEE,φˆEB,φˆTB), which can be combined into a minimal-variance estimator. This estimate is given in the form of spherical harmonic coefficients of the CMB lensing conver-
2The exact rotation applied is given in Chapter3. 3http : //healpix.sourceforge.net.
99 Chapter 4. Integrated cosmological probes: extended analysis gence κ in the angular multipole range 8 ` 2048. These are related to the CMB lensing CMB ≤ ≤ potential φCMB through:
`(` 1) κ + φ . (4.1) CMB,`,m = 2 CMB,`,m
We use these spherical harmonic coefficients to create a HEALPix map of resolution NSIDE = 1024 using the HEALPix routine alm2map. The analysis mask derived by the Planck Collabora- tion is provided as a HEALPix map of NSIDE = 2048. We downgrade this map to a resolution of NSIDE = 1024 following the procedure outlined in Ref. [210], which yields a binary analysis mask. We choose the CMB lensing convergence over the CMB lensing potential map since the lensing convergence is more local and should thus be less affected by masking effects arising when computing angular power spectra. The CMB lensing convergence map covers a fraction of sky f 0.67 and is shown in Fig. 4.2. sky =
4.3.3 Type Ia supernovae
We complement the CMB and LSS data with geometrical constraints on the homogeneous Universe from the distance-redshift relation measured from Type Ia supernovae. We use data from the JLA [34], which is a compilation of 740 SNe Ia comprising data from SDSS-II [88, 154, 254, 163, 46], the Supernova Legacy Survey (SNLS) [22, 258], the HST [226, 262] and several low-redshift experiments [34].4 The JLA data consist of SNe Ia light curve parameters which can be used to calculate observed distance moduli.
4.3.4 Hubble parameter
We also add a local H0 measurement from HST [227] to our analysis. We use the Hubble parameter estimate by Ref. [74], which is a revision of the measurement presented in Ref. [227]. Both measurements are derived from Cepheid-calibrated SNe Ia distance moduli but the former uses a revised distance to the anchor NGC 4258 [135] to calibrate the Cepheid distances. 1 1 This analysis constrains the Hubble parameter to be given by H 70.6 3.3 km s− Mpc− , 0 = ± where the uncertainties are 1σ and assumed to be Gaussian.
4.4 Model predictions
The auto- and cross-correlations of the CMB and LSS cosmological probes can be computed theoretically from the primordial power spectrum. In order to compute the power spectra of the cosmological fields δg , γ and κCMB we employ the Limber approximation [169, 147, 148] as in Chapter3. We further assume a flat cosmological model, i.e. Ωk 0. With these i j = approximations the spherical harmonic power spectrum C` between cosmological probes
4The data can be found at: http : //supernovae.in2p3.fr/sdss_snls_jla/ReadMe.html.
100 4.4. Model predictions
31.0 CMB convergence 3.5 60
30
120 60 0 60 120
b 33.5 0
30
60 3.5 36.0 55.5 53.0 50.5 CMB [1] l
42.5 DES weak lensing 0.1000 60
30
120 60 0 60 120
b 40.0 0
30
60 0.0003 37.5 92.5 95.0 97.5 ˆ [1] l | |
Figure 4.2: New maps used in this analysis in addition to the CMB, galaxy clustering and SDSS weak lensing maps of Chapter3. The full-sky maps are in Galactic coordinates and are shown in Mollweide projection while the zoom-in versions are in Gnomonic projection. The CMB lensing convergence map derived from the foreground-removed CMB temperature and polarization anisotropy maps from SMICA is shown in the top panel. It is masked using the analysis mask provided by the Planck Collaboration. The zoom-in shows an enlarged version of the 5 5 deg2 region centered on (l,b) (53°, 33.5°) shown in the map. The bottom panel × = − shows the map of the weak lensing shear modulus γˆ derived from DES SV. Gray regions | | are masked because they are either unobserved or they do not contain any galaxies at our resolution. The zoom-in shows an enlarged version of the 5 5 deg2 region centered on (l,b) × ( 95°, 40°) shown in the map. It is overlaid with a whisker plot illustrating the direction of = − − the shear. Both maps are shown at a HEALPix resolution of NSIDE = 1024.
101 Chapter 4. Integrated cosmological probes: extended analysis i, j [δ ,γ,κ ] at angular multipole ` can be expressed as ∈ g CMB
Z i ¡ ¢ j ¡ ¢ µ 1 ¶ i j c W χ(z) W χ(z) ` /2 C dz P nl k + ,z , (4.2) ` = H(z) χ2(z) δδ = χ(z) where H(z) is the Hubble parameter, χ(z) the comoving distance and c denotes the speed of nl light. Furthermore, Pδδ (k,z) denotes the nonlinear matter power spectrum at redshift z and i 0 ¡ ¢ wave vector k and W χ(z) is the window function for probe i 0.
The window functions for δg and γ are given in Chapter3. Since the CMB lensing convergence is approximately sourced by a single-lens plane located at the last scattering surface with redshift z its window function can be expressed as the single-plane limit of the weak lensing ∗ shear window function. We therefore have
2 3 ΩmH χ(z) χ(z?) χ(z) κCMB ¡ ¢ 0 − W χ(z) 2 , (4.3) = 2 c a χ(z?) where Ωm is the fractional matter density today and a is the scale factor. In our calculations we set z 1090. ? = The power spectra involving CMB temperature anisotropies can also be related to the primor- dial density fluctuations. The expression for the CMB temperature power spectrum is given in Chapter3. The observed CMB temperature anisotropies are further correlated to tracers of the LSS. For the galaxy overdensity and weak lensing shear this cross-correlation is mainly due to the integrated Sachs-Wolfe (ISW) [232] effect and the resulting cross-power spectra are given in Chapter3. The cross-correlation between the CMB temperature anisotropies and the CMB lensing convergence is dominated by the ISW but receives further contributions from Doppler effects arising from bulk velocities of electrons scattering the CMB photons and from the Sunyaev-Zel’dovich (SZ) [259] effect (for a description of these effects see e.g. Refs. [96, 60]). The cross-correlation due to the SZ effect is not observable using the foreground-reduced CMB temperature anisotropy maps from Ref. [210] but the remaining effects are observable. The cross-power spectrum between CMB temperature anisotropies and CMB lensing convergence can be computed from
a a C TκCMB δ δ , (4.4) 〈 T,`m κCMB,`0m0 〉 = ` ``0 mm0 where aT,`m denotes the spherical harmonic coefficients of the CMB temperature anisotropies
∆T(θ) and aκCMB,`0m0 denotes the spherical harmonic coefficients of the CMB lensing conver- gence defined as
Z c κCMB ¡ ¢ κCMB(θ) dz W χ(z) δ(χ(z)θ,z). (4.5) = H(z)
While the observables discussed so far probe cosmic structure formation, SNe Ia mainly probe the background evolution through their distance moduli. The distance modulus µ of a Type Ia
102 4.5. Spherical harmonic power spectra
Table 4.2: Summary of spherical harmonic power spectrum parameters and angular multipole ranges used in this analysis. The first six power spectra are described in Chapter3.
Power spectrum θmax [deg] θFWHM [deg] `-range ∆` TT C` 40 20 [10, 610] 30 δg δg C` 80 40 [30, 210] 30 γ1γ1 C` 10 5 [70, 610] 60 δg T C` 40 20 [30, 210] 30 γ1T C` 10 5 [70, 610] 60 γ1δg C` 10 5 [30, 210] 60
κT C` 40 20 [40, 400] 60 δg κ C` 80 40 [40, 190] 30 γ1κ C` 10 5 [70, 370] 60 γ2T C` 15 7.5 [70, 610] 60 γ2κ C` 15 7.5 [70, 370] 60 γ2γ2 C` 15 7.5 [70, 610] 60
supernova at redshift zSNe is given by µ ¶ dL(zSNe) µ(zSNe) 5log , (4.6) = 10 10[pc] where dL(zSNe) is the luminosity distance to redshift zSNe.
To compute theoretical predictions for all cosmological observables we follow Chapter3. We use the publicly available Boltzmann code CLASS5 [164] to compute the CMB temper- ature anisotropy power spectra and the cross-correlation between the CMB temperature anisotropies and the CMB lensing convergence. For the other observables we use PYCOSMO [221]. As in Chapter3, we calculate the linear matter power spectrum from the transfer func- tion derived by Ref. [78]. In order to compute nonlinear matter power spectra we use the HALOFIT fitting function [251] with the revisions of Ref. [267].
4.5 Spherical harmonic power spectra
Following Chapter3 we use PolSpice6 [265, 264, 57] to measure the demasked spherical harmonic power spectra from the maps. We calculate the auto-power spectrum of the DES SV weak lensing shear map as well as the cross-correlations between the maps discussed in Section 4.3 and in Chapter3 which have overlaps.
5http : //class-code.net. 6http://www2.iap.fr/users/hivon/software/PolSpice/.
103 Chapter 4. Integrated cosmological probes: extended analysis
We do not include the auto-power spectrum of the CMB lensing convergence in this analysis. This is due to the fact that the auto-power spectrum of the CMB lensing convergence estimator is a biased estimate of the CMB lensing convergence auto-power spectrum because it probes both the connected and the disconnected part of the 4-point function of the CMB temperature anisotropies [206]. In order to obtain the power spectrum of the CMB lensing convergence, the auto-power spectrum of the estimator thus needs to be corrected for this disconnected bias [206], which is beyond the scope of this chapter.
In order to compute the power spectra, we follow the method outlined in Chapter3 to estimate the values of the maximal angular scale used by PolSpice θmax and the apodization parameter θFWHM. We validate these settings using the Gaussian simulations described in Appendix 4.C. The demasking procedure used by PolSpice leads to biases in the recovered power spectra, as discussed in Chapter3. The kernels that relate average PolSpice estimates to the true power spectra can be computed analytically for each choice of θmax,θFWHM and we take them into account by convolving all theoretical predictions with these kernels. The choice of angular multipole ranges follows that described in Chapter3 for all power spectra already included in that analysis. The angular multipole range for power spectra involving the CMB lensing convergence follows the conservative choice described in Ref. [211]. The low-` limit for the power spectra is chosen to minimize the impact of mean field corrections; the high-` limit is chosen because of mild evidence for systematic errors at higher multipoles [211]. The chosen bin widths largely follow the conservative binning outlined in Ref. [211](∆` 45) and the = binning scheme in Chapter3, which is chosen to roughly correspond to the width of the PolSpice kernels. Where we choose angular multipole bins broader than necessary this is mainly done to reduce the size of the spherical harmonic power spectrum vector. The binning schemes and PolSpice parameters used for all power spectra are summarized in Table 4.2.
All the spherical harmonic power spectra are computed from the maps of resolution NSIDE = 1024. They are further corrected for the effect of the HEALPix pixel window function and the power spectra involving the CMB temperature anisotropy map are further corrected for the Planck effective beam window function. The uncertainties are derived from the Gaussian simulations described in Appendix 4.C.
The power spectra computed in this work are described in more detail below; for a description of the remaining power spectra the reader is referred to Chapter3.
4.5.1 DES SV cosmic shear
We compute the cosmic shear power spectrum for DES SV using the map and mask described in Section 4.3.1. In order to estimate the contribution of shape noise, we follow Chapter3 and resort to simulations. We generate 100 noise maps by rotating the galaxy shears by a random angle. We then calculate the power spectra of these maps and our estimator of the shape noise power spectrum is given by the mean power spectrum of the noise maps.
104 4.5. Spherical harmonic power spectra
The weak lensing shear E-mode power spectrum for DES SV is shown in the 4,4-panel of Fig. 4.3 and the B-mode power spectrum is shown in the Appendix (Fig. 4.14). The noise level of DES SV is lower than the one from SDSS Stripe 82 as can be seen by comparing panels 4,4 and 3,3. This is to be expected from the higher galaxy number density and smaller measurement noise of DES SV data.
In Appendix 4.E we compare the DES SV cosmic shear power spectra computed from the maps in Galactic and equatorial coordinates. We find discrepancies similar to those found for SDSS Stripe 82 (Chapter3), especially at small angular scales. Since the differences detected are within the uncertainties of the measurement, we use the cosmic shear power spectrum calculated from the maps in Galactic coordinates in our integrated analysis.
4.5.2 CMB temperature and DES SV weak lensing shear cross-correlation
We compute the cross-power spectrum between the DES SV weak lensing shear and the CMB temperature anisotropies using the maps and masks presented in Section 4.3.1 and Chapter 3. As discussed in Chapter3, we choose to compute cross-correlations using the combined masks of the respective probes rather than the single-probe masks. This is due to the fact that the former approach results in a better recovery of the input cross-power spectra in the Gaussian simulations described in Appendix 4.C. We therefore mask both maps with the combination of the single-probe masks, which covers a fraction of sky f 0.0035. sky ∼ The resulting spherical harmonic power spectrum is shown in the 4,0-panel of Fig. 4.3. As can be seen, the noise level is too high to allow for a detection of the ISW from DES SV weak lensing shear. This is in agreement with the results found for SDSS Stripe 82 in Chapter3. We nevertheless include the power spectrum in our analysis since it provides an upper limit to the ISW signal from weak lensing.
In Appendix 4.E we investigate the impact of our choice of fiducial foreground-reduced CMB temperature map by comparing the power spectra obtained using the four different foreground-reduction algorithms employed by Ref. [210]. As can be seen from Fig. 4.16 we find the measured power spectra to be virtually the same.
4.5.3 CMB lensing convergence and galaxy overdensity cross-correlation
To compute the cross-power spectrum between the CMB lensing convergence and the SDSS DR8 galaxy overdensity we use the maps and masks described in Section 4.3.2 and Chapter3. We mask both maps with their combined mask, which covers a fraction of sky f 0.26. sky ∼ The spherical harmonic cross-power spectrum between the CMB lensing convergence and the galaxy overdensity is shown in the 2,1-panel in Fig. 4.3. We see that we clearly detect a nonzero correlation between the CMB lensing convergence and the galaxy overdensity.
105 Chapter 4. Integrated cosmological probes: extended analysis
4.5.4 CMB lensing convergence and CMB temperature cross-correlation
To compute the cross-power spectrum between the CMB lensing convergence and the CMB temperature anisotropies we use the maps and masks presented in Section 4.3.2 and Chapter 3. The combined mask of both probes covers a fraction of sky f 0.65 and we apply this sky ∼ mask to both maps.
The resulting spherical harmonic power spectrum is shown in the 2,0-panel of Fig. 4.3. Com- paring to the results derived in Ref. [211], we find good overall agreement.
In Appendix 4.E, we again compare the power spectra obtained from the different foreground- reduced CMB temperature anisotropy maps and find them to agree rather well.
4.5.5 CMB lensing convergence and SDSS Stripe 82 weak lensing shear cross-correlation
We estimate the cross-power spectrum between the CMB lensing convergence and the SDSS Stripe 82 weak lensing shear map using the maps and masks described in Section 4.3.2 and Chapter3. We mask both maps with their combined mask, which covers a fraction of sky f 0.0064. sky ∼ The spherical harmonic power spectrum is illustrated in the 3,2-panel in Fig. 4.3. We see that the obtained cross-power spectrum is rather noisy and it does not allow for a detection of the correlation between CMB lensing convergence and SDSS Stripe 82 weak lensing shear. This is probably due to the combined effect of a low fractional sky coverage of SDSS Stripe 82 data and significant noise in both the CMB lensing convergence and SDSS Stripe 82 weak lensing shear. We nevertheless include this cross-correlation into our analysis to serve as an upper limit.
4.5.6 CMB lensing convergence and DES SV weak lensing shear cross-correlation
To compute the cross-power spectrum between the CMB lensing convergence and the DES SV weak lensing shear map we use the maps and masks presented in Sections 4.3.1 and 4.3.2. The combined mask of both maps covers a fraction of sky f 0.0037 and we apply it to both sky ∼ maps.
The 4,2-panel in Fig. 4.3 shows the resulting spherical harmonic power spectrum. We see that the signal-to-noise of the cross-correlation is low for the angular scales considered, which we attribute to both a small sky coverage of DES SV data and the noise level in both maps. Nevertheless we include the power spectrum in our analysis since it provides an upper limit to the cross-correlation of the CMB lensing convergence and the DES SV weak lensing shear field.
106 4.5. Spherical harmonic power spectra
] 6000 2
K TT
µ 5000 )[
⇡ 4000 (2 / 3000 TT `
C 2000
+ 1) 1000 ` ( ` 0 46 0.14 0 100 200 300 400 500 600 44 0.12 7 10 4 42 0.10
1.0 10 )
⇥ ) 40 0 0.08 z