ILLUMINATING THE UNIVERSE: USING THE COSMIC MICROWAVE BACKGROUND TO PROBE STRUCTURE AT INTERMEDIATE REDSHIFT
A DISSERTATION SUBMITTED TO THE DEPARTMENT OF PHYSICS AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY
Stephen John Osborne June 2013
© 2013 by Stephen John Osborne. All Rights Reserved. Re-distributed by Stanford University under license with the author.
This dissertation is online at: http://purl.stanford.edu/sz245wk1516
ii I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy.
Sarah Church, Primary Adviser
I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy.
Steven Allen
I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy.
Chao-Lin Kuo
Approved for the Stanford University Committee on Graduate Studies. Patricia J. Gumport, Vice Provost Graduate Education
This signature page was generated electronically upon submission of this dissertation in electronic format. An original signed hard copy of the signature page is on file in University Archives.
iii iv Abstract
The cosmic microwave background (CMB) provides a backlight that allows us to probe structure out to the last scattering surface. We exploit observations of the sky at microwave and sub-mm wavelengths to measure properties of galaxies and galaxy clusters, as well as to search for possible pre-inflationary signals. In Chapters 2–4 we measure the correla- tion between the dark matter distribution and the microwave and sub-mm emission from galaxies to probe the connection between dark and luminous matter at redshifts 1 3. ⇠ CMB photons are gravitationally deflected by dark matter overdensities, with the majority of the 3 arcminute RMS deflection occurring between redshift 2 and 3. The dark matter ⇠ structures that lens the CMB are traced by dusty star-forming galaxies that emit strongly in the infrared, and have a redshift distribution that peaks between redshift 1 and 3. We use observations of the CMB from the Planck satellite to reconstruct the deflection angles with statistical estimators, and we correlate the deflections with observations of the infrared background light at 100-850 GHz. We find that the two signals are strongly correlated, with a correlation coe cient of approximately 0.8, and we use the measured cross spectrum to estimate the minimum mass scale at which dark matter halos host a CIB source, as well as the star formation rate density in three redshift bins between redshift 1 and 7. In Chapter 5 we use the Doppler shift of CMB light scattered by moving galaxy clusters, known as the kinetic Sunyaev-Zeldovich (kSZ) e↵ect, to put a limit on the large-scale velocity distri- bution of a sample of galaxy clusters observed in WMAP CMB data. On 100 Mpc scales cluster velocities relative to the CMB are expected to be small, originating from gravita- tional instabilities. Larger motions could be generated by pre-inflationary inhomogeneities that leave a “tilt” across our horizon, resulting in a uniform matter flow across the horizon. The kSZ e↵ect is sensitive to such a flow, and we use it to constrain the radial and dipole
v velocity of a sample of 736 clusters with mean redshift 0.12, finding no evidence for either. In Chapters 6 and 7 we search for a possible pre-inflationary signal in CMB data. Mod- els of inflation suggest that our current patch of the universe could have been created as a nucleation bubble from a phase of false vacuum eternal inflation. If additional bubbles are produced, then it is possible that one of them intersected our past lightcone at the time of decoupling, imprinting a disk-shaped signal in the CMB. We have searched for this signal in the WMAP data using optimal algorithms that evaluate the exact posterior likelihood in an e cient and computationally fast way. We find no evidence for the signal, and place limits on the curvature perturbation generated by a collision intersecting the last scattering surface.
vi Preface
Each chapter in this thesis, with the exception of the introduction and conclusion, is a complete paper, and was written in collaboration with others. Some of the papers are in the process of journal submission, and the final versions may di↵er from those presented here. The formatting has been changed for consistency between chapters, and additional sections have been added to Chapter 2 to provide background information. The work was done in collaboration, and I will now outline the contributions that I made. Several of the results use data from the Planck satellite, which is a European Space Agency experiment, with a significant NASA contribution, designed to measure CMB tem- perature and polarization anisotropies. The project grew from two proposed missions, the cosmic background radiation anisotropy satellite (COBRAS) and the satellite for measure- ment of background anisotropies (SAMBA), which were combined and selected as the 3rd Medium-sized mission in the Horizon 2000 Scientific Programme. The success of the mis- sion and the high quality of the data are due to the combined e↵ort of a large team of scientists and engineers. I am a member of the Core Team of scientists on the High Fre- quency Instrument (HFI), one of two instruments onboard Planck (the other being the Low Frequency Instrument) and I joined the project after the satellite had been constructed. I was involved in the pre-flight testing of the instrument, analyzing data from tests done at Centre Spatial de Liege` with Jean-Michel Lamarre. I worked with Andrew Lange’s group at Caltech and the Planck group at JPL on several aspects of the low-level data processing. My main contribution was to better determine the HFI detector time response, by analyzing data from tests where the detector bias current was stepped, and by comparing scans of the Galaxy made six months apart. These projects were done in collaboration with Brendan Crill and Guillaume Patanchon.
vii The analysis presented in Chapter 2 was developed over many iterations between my- self and Duncan Hanson, with Olivier Dore´ providing valuable comments, suggestions, and guidance. The paper presented in Chapter 3 is one of the Planck Collaboration papers, and as such represents the work of many people. I made significant contributions to all aspects of the analysis directly related to the work in Chapter 3, including the measure- ment, null tests, and modeling sections. The project was initiated, proposed, and lead by Olivier, who contributed to all aspects of the paper. Ultimately, the lens reconstructions used in the analysis were produced by Duncan; although I made my own reconstructions, we used Duncan’s to be consistent with the Planck lensing power spectrum paper. In addi- tion, Duncan generated the simulations necessary to calculate the estimator normalization 1 and mean-field, and developed and implemented the C pipeline. I performed the cal- culations in the modeling section, however the code to calculate the CIB auto and cross power spectra with the halo model was given to me by Olivier. The work in Chapter 4 is an extension of the modeling work in Chapter 3 and was done with Olivier. Chapter 5 is largely my own work. I wrote the code and performed the analysis, while Sarah Church, Elena Pierpaoli, and Daisy Mak provided valuable comments and sugges- tions. Many of the results were checked by Daisy who developed an independent analysis pipeline. Chapters 6 and 7 were written in collaboration with Leonardo Senatore and Kendrick Smith. The algorithms were developed in collaboration, and it is di cult to separate in- dividual contributions. The theoretical background was largely written by Leonardo, and the frequentist analysis was primarily developed by Kendrick. I wrote the code, with the 1 exception of the C function and the code to calculate the bubble profiles smoothed by the CMB transfer function, both of which were written by Kendrick.
viii Acknowledgments
Without the help, advice, and encouragement of many people this work would not have been successfully completed. The main person I have to thank is my advisor, Sarah Church, who has supported me throughout, providing advice and assistance, and ensuring that I always had funding to pursue interesting research ideas. Thanks to Sarah, I was able to participate in several observing runs, first to Mauna Kea with the SuZIE experiment, and later to Chile with the QUIET telescope. I am grateful to Charles Lawrence for letting me join the Planck Collaboration and allowing me to work at JPL for extended periods of time, and to Andrew Lange who let me work with his group at Caltech when I first joined the Planck team. I owe a huge amount to Olivier Dore,´ who advised me on all aspects of research, helping me to fully participate in the Planck scientific analysis and introducing me to several of his collaborators. It was a pleasure writing papers, proposals, and presentations with Olivier, and I learned a lot about galaxies, star formation, and the infrared background. I visited Caltech and JPL many times to meet with Olivier, and every trip was engaging. I enjoyed working with Duncan Hanson on the Planck lensing analysis and learning about statistical estimators, and am grateful for the significant amount of time he has spent answering my— often trivial—questions. It was always entertaining visiting Duncan at Caltech and taking his dog on walks. After I started working with Duncan and Olivier, analyzing Planck data became much easier and more fun! I enjoyed collaborating with Elena Pierpaoli and Daisy Mak on the bulk flow project. They provided invaluable help and guidance, especially when trying to get the paper pub- lished. I would like to thank both Leonardo Senatore and Kendrick Smith for making the bubbles project so much fun. I learned a lot from Leonardo about how to approach
ix problems, and I am in awe of Kendrick’s knowledge of statistics. In addition, Leonardo, Kendrick, and Olivier organized a superb conference in Aspen! In my first three years at Stanford I shared an o ce with Ed Wu, and I want to thank him for teaching me about computer science and data analysis—I learned a lot discussing research ideas with Ed. I have had the privilege of working with a talented group of ex- perimentalists in Sarah’s group, and I have benefited from working in the lab with Matt Sieth, Patricia Voll, Kiruthika Devaraj, and Judy Lau, learning how to design experiments and build instruments. I want to thank Keith Thompson for working with me to design and build the QUIET star-tracker, and for teaching me about telescope design, as well as how to use the telescopes at the student observatory. Dana Volponi and Maria Frank were friendly and supportive, and ensured that I successfully navigated the Stanford bureaucracy. Finally, I want to thank my beautiful and brilliant girlfriend Becky Du↵ett, who has put up with many late nights working, and many months travelling. In addition, I am grateful to Becky’s family. The Du↵etts welcomed me for many Christmases and holidays in Portland when I couldn’t make it home.
x Contents
Abstract v
Preface vii
Acknowledgments ix
1 Introduction 1 1.1 Our Cosmology ...... 1 1.1.1 The Friedmann Equations ...... 5 1.2 The Early Universe ...... 6 1.2.1 Inflation ...... 6 1.2.2 Eternal Inflation ...... 8 1.3 Growth of Structure ...... 10 1.4 The Cosmic Microwave Background ...... 14 1.4.1 The Sunyaev Zeldovich E↵ect ...... 15 1.5 The Cosmic Infrared Background ...... 16 1.6 Lensing of the CMB ...... 19 1.7 The Planck Mission ...... 20
2 Foreground Contamination in CMB Lens Reconstruction 39 2.1 Introduction ...... 40 2.2 Lens Reconstruction ...... 43 2.2.1 Lensing Potential ...... 43 2.2.2 Lensing Power Spectrum ...... 45
xi 2.3 Source Trispectrum ...... 46 2.4 Source Estimators ...... 49 2.5 Foreground Simulations ...... 52 2.6 Source Modeling ...... 54 2.6.1 Radio Point Sources ...... 58 2.6.2 SZ ...... 62 2.6.3 Dusty Galaxies ...... 65 2.7 Additional Bias Terms ...... 67 2.7.1 Inhomogeneous Flux Sensitivity ...... 67 2.7.2 Source Masking ...... 69 2.8 Results and Discussion ...... 71 2 A Flat Sky Expressions ...... 83 2 B Supplementary Plots ...... 85 2 B 1 Estimator Sensitivity ...... 85 2 B 2 Frequency Dependence of the Foreground Bias ...... 85 2 B 3 Bispectrum- and Trispectrum-Related Spectra ...... 87 2 B 4 Results on Simulations ...... 88
3 Gravitational Lensing by Dusty Star-Forming Galaxy Halos 93 3.1 Introduction ...... 94 3.2 Data Sets ...... 96 3.2.1 Planck Maps ...... 96 3.2.2 External Data Sets ...... 99 3.3 Cross-Correlation Formalism and Implementation ...... 101 3.3.1 Reconstructing the CMB Lensing Potential ...... 101 3.3.2 Decreasing the Foreground Noise ...... 102 3.3.3 Measuring Cross-Correlations ...... 105 3.3.4 Simulating Cross-Correlations ...... 106 3.4 A Strong Signal Using Planck HFI Data ...... 107 3.5 Statistical and Systematic Error Budget ...... 109 3.5.1 Statistical Error Budget ...... 109
xii 3.5.2 Instrumental and Observational Systematic E↵ects ...... 111 3.5.3 Astrophysical Contamination ...... 115 3.5.4 Final Statistical and Systematic Error Budget ...... 129 3.6 Interpretation and Discussion ...... 129 3.6.1 Model Comparison Methodology ...... 130 3.6.2 Two Modelling Approaches ...... 131 3.6.3 Interpreting the Reconstructed Emissivities ...... 140 3.6.4 Discussion and Outlook ...... 141 3 A Statistical Errors ...... 163
4 Reconstructing The Redshift Distribution Of The CIB Emissivity 167 4.1 Introduction ...... 168 4.2 Reconstructing the CIB Emissivity ...... 169 4.3 Fisher Matrix Analysis ...... 172 4.3.1 Noise Covariance Matrix ...... 173 4.3.2 Experimental Model ...... 174 4.3.3 Principal Components ...... 176 4.4 Results ...... 177
5 Measuring The Galaxy Cluster Bulk Flow From WMAP Data 184 5.1 Introduction ...... 185 5.2 Theory ...... 188 5.2.1 SZ Clusters as Tracers of the Velocity Field ...... 188 5.2.2 Expected Signal for a Generic Velocity Tracer ...... 190 5.3 Data ...... 193 5.3.1 The X-ray Selected Cluster Sample ...... 194 5.3.2 External Simulations ...... 195 5.4 Method ...... 195 5.4.1 Outline of the Method ...... 195 5.4.2 Cluster Modeling ...... 196 5.4.3 Filters ...... 199 5.4.4 Dipole Fitting Procedure ...... 210
xiii 5.4.5 Conversion to Velocity Dipole ...... 212 5.4.6 Tests of the Method ...... 213 5.5 Systematic E↵ects ...... 214 5.5.1 Thermal SZ ...... 214 5.5.2 Radio Point Sources ...... 219 5.5.3 Galactic Emission ...... 220 5.6 Results ...... 224 5.6.1 Dipole ...... 224 5.6.2 Monopole and Higher Moments ...... 227 5.6.3 Comparison with SuZIE Measurements ...... 228 5.7 Conclusions ...... 228 5 A Expected Cluster Velocity Dipole ...... 241 5 B Comparison with KAKE Results ...... 242 5 B 1 Filter Pipeline Comparison ...... 242 5 B 2 Cluster Dipole Comparison ...... 244
6 Optimal Analysis of Azimuthal Features in the CMB 248 6.1 Introduction ...... 249 6.2 Mini-Summary of the Theory of Bubble Collisions ...... 252 6.3 Bubble Signal ...... 255 6.3.1 Ramp and Step Models ...... 255 6.3.2 CMB Temperature Profiles ...... 256 6.3.3 Can CMB Transfer Functions be Neglected? ...... 258 6.3.4 Discretizing the Bubble Radius ...... 260 6.4 Data Analysis ...... 262 6.4.1 Definitions and Notation ...... 262 6.4.2 An Algorithm for Fast Calculation of 2 ...... 264 6.4.3 Bayesian Analysis ...... 267 6.4.4 Frequentist Analysis ...... 270 6.4.5 Computational Cost ...... 274 6.5 Large Numbers of Bubbles ...... 275
xiv 6.6 Conclusions ...... 278 1 6 A Details of the C Filter ...... 284 6 B Likelihood Model ...... 286
7 Collisions with Other Universes: the Optimal Analysis of the WMAP Data 290 7.1 Introduction ...... 291 7.2 Bubble Signal ...... 292 7.3 Method ...... 293 7.3.1 Bayesian Analysis ...... 294 7.3.2 Frequentist Analysis ...... 297 7.4 Many Bubbles ...... 299 7.5 Conclusions ...... 299
8 Summary 302 8.1 Lensing and the Cosmic Infrared Background ...... 302 8.2 Kinetic SZ Velocity Constraints ...... 304 8.3 Cosmic Bubble Collisions ...... 305
xv List of Tables
1.1 Evolution of a single mode of the potential , baryon density b, dark mat-
ter density c, and radiation temperature ⇥r, based on Ma [19]. Each mode begins evolving when it enters the horizon, which may not occur before decoupling. D(a) is the growth function...... 13
2.1 Primary source terms discussed in Sec. 2.3. The first row gives the trispec- trum due to lensing for comparison purposes...... 48 2.2 One-point amplitudes calculated at 150 GHz, with a 200 mJy (6 mJy in brackets) maximum flux density. For radio point sources S is not shown h i since Eq. (2.45) gives an unphysical negative flux...... 62
3.1 Point source estimator. The measured quantity S 3 , as defined in Eq. 9, is given as a function of frequency...... D . .E ...... 121 3.2 Cross-spectrum detection band-powers. All values are in units of µK.sr. fore The extragalactic foreground contribution, C` has been removed from
C`. Both statistical and systematic errors are given (see Sect. 3.5 for details). 128 3.3 Reconstructed emissivity as a function of redshift and associated star for- mation rate. At each frequency and for each of the three redshift bins the first quantity corresponds to the mean emissivity in the corresponding red- 1 1 shift bin, ¯j(z), in Jy Mpc sr , while the second corresponds to the SFR 3 1 density, ⇢SFR, in M Mpc yr ...... 143
5.1 Simulated thermal SZ monopole and measured monopole...... 200 5.2 Monopole and dipole in the PSM simulations for the z = 0 1 shell. . . . . 223 xvi 5.7 95% Confidence Limits on the Dipole...... 227
5.3 Results from the KAKE filter.⇤ ...... 239
5.4 Results from our Wiener filter.⇤ ...... 239 5.5 Results from the matched filter...... 240 5.6 Results from the tSZ bias removing filter...... 240 5.8 Number of clusters satisfying the luminosity and redshift cut in our sample and the KAKE sample...... 245 5.9 Comparison of the cluster monopole and dipole in our cluster sample with the KAKE sample...... 245 5.10 Same as table 5.9 for the W band channels only...... 246 5.11 Result from Kashlinsky et al. [54]...... 246 5.12 Mean optical depth of our cluster sample...... 247 5.13 Result from our filter pipeline with the luminosity cut...... 247
xvii List of Figures
1.1 Simulation of the spatial distribution of bubbles at a late time in an eternally inflating false vacuum. From [16]...... 9
1.2 Map of CMB temperature anisotropies measured by the WMAP satellite (WMAP Collaboration)...... 15
1.3 Power spectrum of CMB temperature anisotropies measured by the WMAP satellite (WMAP Collaboration)...... 16
2 1.4 Maps of the 26 deg N1 field (galactic coordinates l, b = 85.33 , 44.28 )
at 217, 353, 545, and 857 GHz, smoothed with a 100 Gaussian, with the CMB, bright point sources, and galactic cirrus subtracted. From Planck Collaboration et al. [39]...... 18
1.5 Power spectrum of the CMB lensing potential in the ⇤CDM model. The solid line is the linear theory model, with the dashed line also including non-linear corrections from HALOFIT [44]. From [43]...... 20
1.6 Theoretical temperature and polarization power spectra. The primordial B-mode signal is shown in blue with two di↵erent values for the tensor to scalar ratio. The dashed curves are expected noise levels at 143 GHz (red) and 217 GHz (blue) for two sky surveys...... 22
xviii 2.1 Trispectrum-related power spectra biases ˆab calculated for an L = CL cd max 1500 experiment. The amplitudes of the four trispectrum shapes are arbi- trary (such that ˆab 1 at its peak.), but are used consistently in each of CL ab ⇡ the four plots. Note that the relative size of the bias terms in each plot has no physical significance. The colored solid lines give the total contributions to ˆab from the indicted trispectra. Corresponding dashed lines give only CL the contribution from the primary contraction of the trispectrum (which in the upper left panel is equal to the lensing potential spectrum)...... 52
2.2 Same as Fig. 2.1, but for an experiment which is cosmic variance (C.V.)
limited to Lmax = 3000...... 53
2.3 Temperature power spectra of the simulated sky components for the Lmax = 1500 and 3000 experimental models, including unlensed and lensed CMB (black and orange), radio sources (red), infrared sources (blue), thermal SZ (green), and kinetic SZ (purple) at 148 GHz. Bright sources have been removed before taking the spectra (see Sec. 2.5 for details)...... 54
2.4 Power spectra calculated using our foreground model. The solid lines show the CMB (unlensed in black, lensed in orange), radio point source (red), SZ (green), and CIB (blue) spectra. Radio and CIB point sources with flux density S > 200 mJy are assumed to be masked. The dashed lines show the point source signals for a 6 mJy flux density cut. For comparison the blue dotted lines show the CIB one and two-halo terms calculated using the Ade et al. [27] CIB model...... 59
2.5 Redshift distribution of power in the radio point source (red), SZ (green), and CIB (blue) spectra calculated using our foreground model at 150 GHz for ` = 50 (solid) and ` = 2000 (dotted). The power in the lensing potential spectrum is shown in black. The curves are normalized by area...... 61
2.6 As in Fig. 2.5, except that the foreground spectra show the redshift distri- S 2S 2 bution of power in C` ...... 61
xix 2.7 Power spectrum of S 2 from simulation (colored lines) and from our model (solid black lines) for radio point sources (upper panel), infrared sources (middle panel), and SZ (lower panel). The infrared model has been scaled by a factor of 30 to better fit the simulation results. The additional dashed and dotted lines are results from models with di↵erent parameters, dis- cussed in the text...... 63
2.8 Cross power spectrum of S 2 and calculated from the simulated S and maps (colored boxes) and the model predictions (black lines). To better fit the simulations the CIB model spectrum has been multiplied by 25 for the B12 model and by 10 for the B10 model. The SZ spectrum has been multiplied by 2. The empirical model from Sec 2.4 is almost identical to the radio point source model up to an amplitude...... 64
2.9 Contribution to the reconstructed lensing spectrum from simulations for
each foreground component when Lmax = 1500 (red line). The blue line gives the result after subtracting the disconnected term, which is the green line. The residual foreground contribution calculated using the results from Sec. 2.3 is the orange line. The blue bars give the lensing spectrum recon- structed from the full-sky simulated CMB map...... 69
2.10 Same as Fig. 2.9, but for Lmax = 3000. The di↵erence between lensing the reconstructed lensing potential and the simulated lensing potential is due to the N(1) noise bias [17], which has not been subtracted...... 70
2.11 Detection significance as a function of the minimum (black) and maximum (blue) multipole used in the estimator. The black curve is calculated with
`max = 3000 and the blue curve with `min = 2...... 86
ˆ 2.12 Bias to C` calculated using the Lmax = 1500 (upper panel) or Lmax = 3000 (lower panel) estimator and masking scheme from the simulated radio point source (red), infrared point source (blue), thermal SZ (green), and kinetic SZ (purple) maps. The gray bars show the error bars from a Fisher matrix calculation...... 90
xx 2.13 The bispectrum and trispectrum-related spectra calculated from the 148 GHz simulated map. The red line shows the best-fit spectrum, which depends on a single parameter: S 3 for the bispectrum, and S 4 for the trispectrum. The noise bias has beenD E subtracted from the trispectrumD E estimate, and the bias terms caused by CMB lensing are shown in blue...... 91
2.14 The green line shows the bias from the unresolved foreground components in the masked 148 GHz simulated map. The solid blue line is the mea- sured bias, the dashed blue line is the disconnected term, and the red line is the sum of the blue lines. The gray error bars are from a Fisher matrix
calculation for the Lmax = 1500 estimator...... 92
3.1 Redshift- and mass- integrand for the CIB and CMB lensing potential power spectra at ` = 500, calculated using the CIB halo model of [17], evaluated at 217 GHz. The good match between the redshift and halo mass distri- butions leads to an expected correlation up to 80 %. The sharper features in the CIB kernel are artefacts from the [27] model. We note that the low mass, high z behavior of our calculation is limited by the accuracy of the 5 mass function we use [28]. All of our mass integrals use Mmin = 10 M .. . 97
3.2 Combined Galactic, point-source and Hi mask with sky fractions 16, 30 and 43 %...... 100
3.3 Angular cross-spectra between the reconstructed lensing map and the tem- perature map at the six HFI frequencies. The error bars correspond to the scatter within each band. The solid line is the expected result based on the PER model and is not a fit to these data (see Fig. 3.16 for an adjusted model), although it is already a satisfying model. In each panel we also show the correlation between the lens reconstruction at 143 GHz and the 143 GHz temperature map in grey. This is a simple illustration of the fre- quency scaling of our measured signal and also the strength of our signal as compared to possible intra-frequency systematic errors...... 104
xxi 3.4 Temperature maps of size 1 deg2 at 545 and 857 GHz stacked on the 20,000 brightest peaks (left column), troughs (centre column) and random map lo- cations (right column). The stacked (averaged) temperature maps is in K. The arrows indicate the lensing deflection angle deduced from the gradient of the band-pass filtered lensing potential map stacked on the same peaks.
The longest arrow corresponds to a deflection of 6.300, which is only a frac- tion of the total deflection angle because of our filtering. This stacking allows us to visualize in real space the lensing of the CMB by the galaxies
that generate the CIB. A small and expected o↵set ( 10) was corrected by ' hand when displaying the deflection field...... 107
T 3.5 Naive analytical estimates of the contribution to the C` variance as a func- tion of multipole and frequency as given in Eq. 4. We assume the same bin sizes as in Fig. 3.3. The di↵erent lines are the contribution to the ana- 2 CIB CIB lytical error from the signal only: C` C` + C` (green), noise only: ˆ ,N ˆCIB,N ˆCIB,N C` C` (blue), and the mixed signal and⇣ noise terms:⌘ C` C` (yel- ˆ ,N CIB low) and C` C` (orange). The total contribution is the solid black line, CIB 2 and the theory spectrum, C` , is the dashed line...... 111 ⇣ ⌘
3.6 Null tests at 545 GHz. Left: di↵erence spectra obtained by nulling the sig- nal in the HR temperature map before correlating it with our nominal reconstruction. Middle: temperature signal nulled using di↵erent detectors at 545 GHz. Right: temperature signal nulled using the first and second sur- vey maps. The black error bars correspond to the scatter measured within an `-bin, while the coloured bands correspond to the analytical estimate. Except for the survey null test (see text for details), these tests are passed 2 satisfactorily except, as illustrated by the quoted and Ndof, thus strength- ening confidence in our signal...... 115
xxii 3.7 Left: di↵erence between the cross-spectra measured using the 20 % Galac- tic mask (20 % is the unmasked sky fraction) from that measured with our default 40 % Galactic mask. Middle: spectra obtained when di↵erencing the 60 % and 40 % Galaxy mask measurements. For both left and mid- dle panels and all Galactic masks, the same point source and Hi masks are used, which removes an additional fraction of the sky. Right: di↵erence between the cross-spectra calculated with the Hi cleaned temperature maps from those with no Hi cleaning. This cross-spectrum is thus the correlation between the Hi template and the reconstruction. The error bars are cal- culated in the same way as in Fig. 3.6. Again, the null tests are passed with an acceptable 2...... 116
3.8 Left: di↵erence between cross-spectra calculated using the lens reconstruc- tion at 100 GHz with the nominal 143 GHz reconstruction. We see an over- all shift, which leads to a high reduced 2. This shift can be explained by the expected overall normalization uncertainties of the 100 GHz and 143 GHz reconstructions. While this uncertainty is not included in the 2 or the solid bars, it is included later in our analysis in Sect. 3.6.1. Middle: same as the left panel, but the 217 GHz reconstruction is used instead of the 100 GHz reconstruction. Right: di↵erence between cross-spectra when we consider the 143 GHz lens reconstruction calculated with a less restric- tive Galaxy mask (that excludes only 20 % of the sky) and the nominal reconstruction mask that excludes 40 % of the sky...... 117
3.9 Results from the point source contamination estimator of Eq. (5). The best- fit cross-spectra associated with shot noise are plotted in black. We do not show the best-fit at 545 and 857 GHz since the signal-to-noise ratio is low. The grey line is a prediction for the bias from the CMB lensing - infrared correlation, and has been subtracted from the spectra (plotted as black points). We see that with the subtraction of the bias from CMB lensing, the measured bispectrum-related spectrum is generally consistent either with zero, or with the shape expected for shot noise...... 120
xxiii 3.10 Frequency spectrum of our cross-spectra averaged within `-bins (black points with associated error bars). The light shaded regions correspond to the HFI frequency bands. The solid black curve corresponds to the best- fit CIB assuming a [64] spectrum, while the dot-dashed line assumes a [9] spectrum. The dashed black line corresponds to the best-fit model when allowing for an SZ component in addition to the [64] CIB shape. The blue dots correspond to the associated absolute value of the best-fit SZ com- ponent. We conclude from this plot that the SZ e↵ect is not an important contaminant...... 122
3.11 Foreground components at each frequency. The data points and error bars show our results. The dashed line is an estimated upper limit on the mag- nitude of the SZ contamination derived in Sect. 3.5.3. We show the abso- lute value of this contribution, which is negative at frequencies less than 217 GHz. The dot-dashed line is the extragalactic point source contribu- tion, with an amplitude measured from our data as derived in Sect. 3.5.3. Again we show the absolute value, with the signal being negative below ` 1200. The oscillating solid line corresponds to the calculated ISW ⇠ contamination...... 125
3.12 Cross-spectrum of the 545 GHz lens reconstruction correlated with the 545 GHz temperature map with di↵erent Galactic masks. The legend gives the vis- ible sky fractions. The solid line represents the analytic unclustered shot-
noise contribution fit to the fsky = 0.09 points above ` = 1300...... 126
3.13 Cross-correlation coe cients calculated from the model spectrum and best-fit halo model at each frequency. The CIB is a spectacular tracer of CMB lensing, and vice-versa. The data points represent the measured cross-correlation divided by the best-fit auto power spectra models at 545 GHz.132
3.14 Marginalized 2-D distributions of zc and z for the linear bias model, fit to all frequencies simultaneously. The orange dots indicate the parameter values at the minimum 2...... 133
xxiv 3.15 Measured cross-spectra with the best-fit j reconstruction model fit to both the CIB auto- and CIB-lensing cross-spectra (solid coloured), and the best- fit linear bias model (dashed coloured). The 2 values quoted in each panel are the contribution to the global 2 from the data in the panel for the halo model, and loosely indicate the goodness of fit (see text for details). The one and two-halo contributions are shown as the dashed and solid black lines, respectively. A light dashed black horizontal line is indicating the zero level...... 136 3.16 PIR auto-spectra with the best-fit mean emissivity j reconstruction model fit for the CIB auto and CIB-lensing cross-spectra (solid coloured). The 2 values are defined as in Fig. 3.15. The one and two-halo contributions are shown as the dashed and solid black lines, respectively, while shot noise is the dot-dashed black line...... 137
3.17 Marginalized 2-D distribution of log10 (Mmin/M ) and log10 (Msat/M ) for our overall HOD model when the CIB-lensing cross-spectra are combined with the CIB auto-spectra and FIRAS measurements. The orange dot shows the best-fit value. The contours correspond to 68%, 95% and 99.7% confi- dence intervals...... 139 3.18 Reconstructed mean emissivity, ¯j, for each frequency as a function of red- shift. The solid line at low z and the dashed line at higher z correspond to the B11 model. The B11 emissivity model at z > 1 is not used, and is shown only for reference. The black error bars correspond to the 68% C.L. while the color shading display the full posterior distribution...... 142 3.19 Marginalized 1-D distribution of the emissivity in the high redshift bin at 353 GHz with (black line) or without (blue line) including the CIB-lensing correlation. Its inclusion helps to constrain the emissivity at high redshift, transforming an upper limit into a detection...... 143 3.20 Correlation between the lensing potential and the IRIS map at 100 µm using our nominal lens reconstruction. We clearly see a correlation and estimate the significance to be 9 , ignoring possible systematic e↵ects. The solid line represents a simple reasonnable prediction for this signal...... 144
xxv 3.21 Ratio of various error estimation procedures to the errors obtained with the data-based analytical estimate. At each frequency the numerator is given by: (i) the scatter within an `-bin in simulations (solid black line); (ii) the scatter within an `-bin in the data (solid dashed black line); (iii) the scatter of bins across simulated realizations (solid coloured line); (iv) the analytical errors calculated from the simulations (dashed coloured line); (v) the scatter across realizations for the cross-correlation between the sim- ulated temperature map and the lensing potential reconstructed from the data (coloured dot-dashed line). The grey envelope is the precision of the simulated errors expected from 100 simulations (shown as a spread around unity)...... 163
4.1 Derivatives of the CIB auto spectrum (top and middle panels) and CIB lens- ing potential cross spectrum (bottom panel) with respect to the amplitudes,
ji, of the binned emissivity...... 171
4.2 CIB auto and CIB lensing potential cross spectra relative to our fiducial model at 217 GHz assuming 10 µK-arcmin noise, 5’ beam FWHM and
Lmax = 3000 (dark gray bars), and 100 µK-arcmin noise, 5’ beam FWHM
and Lmax = 1500 (light gray bars). The colored lines show the spectra if the emissivity in a single redshift bin is doubled while holding the mean CIB intensity constant...... 174
4.3 The five (normalized) principal components of the CIB emissivity with the largest eigenvalues, obtained with experimental model 1. In the upper panel we assume no knowledge of the frequency dependence of j, while in the lower panel we assume that the frequency dependence is known exactly. Note the di↵erent y-scale in the two panels...... 175
xxvi 4.4 Reconstructed emissivity relative to the true emissivity for experimental model 1 (solid) and model 2 (dashed), assuming that the frequency depen- dence of j is precisely known. The colored lines give the reconstruction calculated using only the principal components measured with signal-to- noise ratio greater than 3, as determined from a Fisher matrix calculation, using only the CIB auto spectra, auto and cross frequency spectra, or CIB lensing potential cross spectra...... 177
5.1 Redshift distribution for the cluster sample used in this paper...... 191 5.2 Expected bulk flow velocity in ⇤CDM cosmology with a selection function calculated from our cluster sample in redshift shells extending from z = 0 to the specified x-axis value. The dashed line shows the result with a uniform selection function. The dark and light shaded areas are the 68% and 95% confidence limits from cosmic variance...... 192 5.3 Multi-frequency matched filters for the WMAP channels...... 202 5.4 Unbiased multi-frequency matched filters. A discussion of the filter shape is given in the text...... 202 5.5 Thermal SZ bias removing filters for an experiment with four frequencies:
100 GHz, 143 GHz, 217 GHz and 353 GHz, beam FWHM of 100,70,50
and 50, and white noise levels 25, 15, 25 and 75 µK/K arcmin...... 204 5.6 Spectra of a simulated kSZ map (with signal only at the locations of our cluster sample) convolved with the WMAP beams. The simulations were performed both with (top 3 spectra) and without (bottom 3 spectra) a bulk flow component of amplitude 1000 km/s. The dashed lines are the WMAP beam functions, scaled to the signal amplitude...... 205 5.7 Sum of the eight WMAP maps in the Q, V and W bands filtered by the matched filters in Fig. 5.3. The map is noise dominated...... 206 5.8 Sum of the eight WMAP maps in the Q, V and W bands filtered by the tSZ removing filters in Fig. 5.4. The values in this map are larger than in Fig. 5.7 since the map is noise dominated and the unbiased filters give a lower signal to noise measurement of the kSZ signal...... 206
xxvii 5.9 Wiener filter used by KAKE...... 208
5.10 Wiener filter for each of the WMAP channels...... 208
5.11 Suppression of the kSZ signal by the beam and filter in the V band...... 209
5.12 Left: Error in the recovered dipole direction from simulated kSZ maps with a bulk flow velocity of 0 km/s (top), 500 km/s (middle) and 1000 km/s
(bottom) in the directions: (latitude, longitude) = (0 ,0 ) black lines, (0 ,
90 ) blue lines and (90 ,0 ) red lines, using the MF. Right: Same for the UF. The arrows are the 95% confidence limits...... 215
5.13 Recovered bulk flow velocity in simulated maps containing CMB, noise and kSZ using the MF (black) and the UF (blue). All clusters are included in the fit. The x-axis is the bulk velocity input into the simulated maps, the y-axis is the recovered velocity. The red line indicates perfect recovery. When the input bulk flow is small, a dipole is not detected and the data points provide an estimate of the scatter in the recovered amplitude. . . . . 216
5.14 Error in the direction of the bulk flow velocity in simulated maps containing CMB, noise and kSZ with the MF (black) and UF (blue), when all clusters are included in the fit...... 217
5.15 Error in the recovered monopole velocity in simulated maps containing CMB, noise and kSZ with the MF (black) and UF (blue)...... 218
5.16 Sum of the PSM simulated maps convolved with the WMAP beams and filtered with the MF. This map shows galactic emission around the edges of the mask that is absent from the WMAP foreground reduced maps due to WMAP’s foreground reduction method [6]. The color scale is altered to better show the galactic emission by mimicking the e↵ect of replacing the 1 data by sinh (data) [29]...... 221
5.17 Sum of the PSM simulated maps convolved with the WMAP beams and 1 filtered with the UF. The sinh color scale that was used in Fig. 5.16 is used here...... 221
xxviii 5.18 Left: Cluster dipole amplitude in the maps filtered by the KAKE filter. Points with a plus sign have redshift shells with minimum redshift of 0, those with a cross have shells with minimum redshift of 0.05 and 0.12. The green line is the noise bias, the red line is the 95% confidence limit that there is no bulk flow and the blue line is the 99.7% confidence limit. Right: Cluster dipole amplitude in km/s...... 225
5.19 Cluster dipole amplitude in the maps filtered by the WF...... 225
5.20 Cluster dipole amplitude in the maps filtered by the MF...... 226
5.21 Cluster dipole amplitude in the maps filtered by the UF...... 226
5.22 95% confidence upper limit to the bulk flow in the redshift 0 1 shell. We find a bias in simulated thermal SZ maps equivalent to a bulk velocity of 2500 km/s...... 227 ⇠ 5.23 Spectrum of the W1 channel map filtered by the KAKE filter (black) and the di↵erence between the filtered W band map spectra from our KAKE fil- ter pipeline and the spectra of the publicly available maps used in the Kash- linsky et al. [54] analysis...... 243
5.24 Filtered W1 channel map spectrum from the KAKE pipeline and our pipeline.243
5.25 Integrated signal from a beam smoothed and filtered cluster in our simu- lated maps...... 244
6.1 A spacetime diagram showing the causal structure of a single bubble col- lision, based on Fig. 3 from [15]. Coordinates are chosen so that light
propagating in the plane of the diagram moves along 45 lines...... 253
6.2 The Earth’s last scattering surface at the time of decoupling, based on Fig. 4 from [15]. The color of the shaded region indicates the magnitude of the curvature perturbation, assumed to be of the “ramp” form given below in Eq. (4)...... 253
xxix 6.3 Angular bubble profiles b(✓) defined in 6.3.2, for bubbles at comoving dis- § tance r = 13886.6 Mpc (corresponding to angular size ✓bubble = 11.39 ) and arbitrary normalization. The ramp and step profiles are obtained by evolv- ing an initial feature in the adiabatic curvature perturbation (Eqs. (4), (5)) forward to obtain a CMB temperature profile, using the CMB transfer func- tion. The cosine and disc profiles are approximations to the ramp and step profiles in which the CMB transfer function is omitted (Eq. (14)). Qual- itatively, the e↵ect of including the CMB transfer function is to smooth
sharp features in the profile, with smoothing length 1 given by the CMB ⇡ acoustic scale...... 261
6.4 The WMAP 7-year best fit ⇤CDM CMB spectrum (black line), the spec- trum of bubbles with the ramp profile normalized to the CMB spectrum at ` = 2 (blue), and the spectrum of bubbles with the step profile normalized in the same way (red)...... 275
6.5 Left panel: Histogrammed maximum likelihood bubble amplitude aML, for
many Monte Carlo simulations of the data d. Counterintuitively, aML is within 0.001 of zero 11% of the time. Right panel: Posterior likelihoods (d a) for three randomly simulated data realizations d. It is seen that the L | derivative (@ /@a) is always zero at a = 0, but this point can be either the L global maximum likelihood (blue, solid curve), a local maximum which is not the global maximum (magneta, dotted curve), or a local minimum (red, dashed curve). Likelihoods in this figure were generated using a toy model
with 1000 µK-arcmin isotropic noise, 5 Gaussian beam, no sky cut, and
assuming fixed bubble size ✓bubble = 30 ...... 286
xxx 7.1 Bayesian analysis of the bubble parameter space, assuming the “ramp model” for the bubble profile. Top panel: Posterior likelihood (aramp d) for the L | amplitude parameter aramp, defined to be the slope of the initial curva- 1 ture perturbation in Mpc , given the WMAP data d (solid black), after marginalizing the bubble radius. As explained in the text, the tails of the likelihood are slow to decay, due to a poorly constrained region of parame- ter space with small bubble radius. We illustrate this by showing the likeli-
hood calculated using bubbles with a subset of angular sizes: ✓bubble < 20
(blue short-dashed), and ✓bubble > 20 (red long-dashed). Bottom panel: Posterior likelihood (↵ramp d), obtained from the top panel by changing L | ramp ramp ramp 4/3 variables from a to ↵ = a (sin ✓bubble) . After this change of variables, the likelihood is narrower and less sensitive to marginalization over the bubble radius. Vertical lines are 95% confidence limits on the amplitude parameter ↵ramp. The likelihood is consistent with no bubbles (↵ramp = 0)...... 295 7.2 Bayesian analysis of the “step” bubble model. We show the posterior likeli- hood (↵step d) after changing variables from astep to ↵step = astep(sin ✓ )1/3 L | bubble to remove degeneracies. The dashed coloured lines have the same meaning as in Fig. 7.1. The likelihood is consistent with no detection of the signal. . 297 7.3 Distribution of the likelihood ratio from simulations. An x-axis value of zero means that the likelihood peaks at a = 0. The red up-shaded distri- bution is consistent with the likelihood peaking at a = 0 within numerical precision. The WMAP value is consistent with a peak at a = 0, which occurs in 51.4% of the simulations...... 298
8.1 Comparison of the Planck and SPT-SPIRE CIB-CMB lensing potential cross spectrum. The SPT-SPIRE data is from Holder et al. [2]...... 304
xxxi xxxii Chapter 1
Introduction
“I dont know if you have had the same experience, but the snag I always come up against when I’m telling a story is this dashed di cult problem of where to begin it.”
Right Ho, Jeeves P.G. Wodehouse
1.1 Our Cosmology
Observations over the past few decades have allowed us to build a relatively simple cos- mological model of the universe that can account for many of the observed structures we see. The ⇤ Cold Dark Matter (⇤CDM) model allows us to explain the structure that we see in cosmic microwave background (CMB) images as well as in galaxy surveys, and it enables us to calculate the age and size of the universe. The model assumes that our place in the universe is not special (the Copernican principle) and that we observe a sample of the universe that is representative of the whole, with the same physical laws at all locations. The following properties are assumed:
Isotropy — The large scale properties of the universe look the same in all directions. • This means that global properties of the large scale structure of the universe can be inferred from observations made in one direction. In addition, isotropy implies that a
1 2 CHAPTER 1. INTRODUCTION
preferred direction or axis of the universe cannot to be chosen, for example on large enough scales objects do not preferentially move towards us in one direction and away from us in another.
Homogeneity — Any property of the universe when averaged over large scales is the • same at all locations. This is observed to be true even for regions of the universe that have, based on the ⇤CDM model, never been in causal contact. The explanation for this is currently unknown, with one possibility being a period of rapid expansion early in the history of the universe.
These assumptions are found to be in excellent agreement with all measurements, and we present one such test in Chapter 5, where we measure the bulk motion of a sample of galaxy clusters. As well as these assumptions, there are several other important features of the universe that must be accounted for by any cosmological model:
Expansion of space — Evidence for the expansion of space was obtained in 1929 • when Hubble observed that galaxies are receding from us with a velocity propor- tional to their distance from us. This can be accounted for in General Relativ- ity by including a cosmological constant in the Einstein equations (which can ei- ther be interpreted as a constant energy density of space, or equivalently as being a new type of matter/energy with negative pressure). The expansion of space implies that the universe starts in a singularity, the Big Bang, and allows us to define an age for the universe. The current expansion rate of the universe is measured to be H = 74.2 3.6 km/s [1] and the dark energy density (defined relative to the critical 0 ± +0.015 density) is ⌦⇤ = 0.728 0.016 [2]. Presence of dark matter — First postulated in 1932 by Jan Oort to account for higher • than expected velocities of stars in the Milky Way given the amount of visible matter present in the Galaxy, and in 1933 by Fritz Zwicky to account for the discrepancy between the Coma cluster mass inferred from its luminosity and the mass inferred from the velocities of galaxies within it. Both observations can be explained by the presence of matter which does not interact electromagnetically with baryonic matter. Any explanation using modified gravity theories must explain the spatial separation 1.1. OUR COSMOLOGY 3
of visible and dark matter observed in colliding galaxy clusters (for example, the Bullet Cluster, Clowe et al. 3). To form structures with the density contrasts observed today, N-body simulations suggest that the dark matter is non-relativistic (or “cold”). The dark matter density has been measured to be ⌦ h2 = 0.1123 0.0035 [2, using c ± the WMAP7+BAO+H0 combination].
Spatial flatness — The evolution of the universe is described by the Friedmann equa- • tions (which we briefly describe below in Section 1.1.1), which derive from Ein-
stein’s equations. They allow us to define a critical density, ⇢c, which is the energy
density of a flat universe. If the energy density is greater than ⇢c, then the universe
is closed, if it is below ⇢c, then the universe is open. Observations of the CMB combined with measurements of the Hubble Constant and Type Ia supernovae show that the curvature is flat to within experimental error, with ⌦ ⌦ + ⌦ = tot ⌘ i i ⇤ 1.002 0.011 [4] where ⌦ is the energy density relative to the critical density. A ± P major focus of this work is to measure gravitational lensing of the CMB, which will ultimately allow more precise measurements of curvature using the CMB alone [e.g. 5].
Combining the Standard Model of particle physics with the ⇤CDM model allows the universe to be described by 31 free parameters [6], with all current cosmological observa- tions being explained with 5 parameters, although more are likely to be needed to accom- modate more precise observations. In addition to the dark matter and dark energy densities above, the additional parameters are the baryonic matter density, ⌦ h2 = 0.02260 0.00053, b ± +0.088 9 the amplitude of scalar fluctuations, As = 2.441 0.092 10 , and the spectral index of ⇥ scalars n = 0.963 0.012 [2]. Additional parameters may be required, such as the tensor s ± ⇣ ⌘ to scalar ratio, r < 0.72 (95% confidence) [7], the tensor spectral index, and the neutrino 2 density, ⌦⌫h = ⌫ m⌫/93 eV with ⌫ m⌫ < 0.17 eV (95% confidence) [8], with a con- straint obtained from the CMB, large scale structure, Lyman ↵ forest, and SNIa data, and P P a lower limit of ⌫ m⌫ & 0.05 from neutrino oscillation experiments, as well as parameters needed for moreP complex models. In Chapter 2 we will present methodology designed to aid measurements of the sum of neutrino masses. 4 CHAPTER 1. INTRODUCTION
While this model has been hugely successful in allowing us to explain observations of structure in the universe it does not tell us the nature of dark energy or dark matter. Several unsolved problems remain:
Horizon problem — The universe appears statistically homogeneous and isotropic on • large scales, and there is no observational evidence to the contrary. We observe the CMB to have the same temperature on opposite sides of the sky, despite the fact that the regions are separated by a distance large enough that, according to the ⇤CDM model, they have never been in causal contact.
Flatness problem — The energy density of the universe is (to within observational • limits) equal to the critical density required to make the universe flat, despite there being no clear reason why this should be the case. The problem is even more serious given that a flat universe is in an unstable equilibrium if the energy density is slightly over the critical value in the early universe, then the universe will become denser relative to the critical density as it evolves. We will show in Section 1.1.1 that the e↵ective curvature density has a (1 + z)2 dependence, whereas the matter density has a (1 + z)3 dependence and the radiation density has a (1 + z)4 dependence. In the past the curvature contribution therefore dominates, and for the universe to be flat today the curvature contribution must have been exponentially small in the early universe.
Absence of defects — If the early universe is hot, then a large number of stable • particles, such as magnetic monopoles, should be produced [e.g. 9]. If they survived until the present day, then they would be the primary constituent of the universe, however, searches for them have not been successful1
We now overview the equations that underpin the ⇤CDM model before discussing a possible solution to the above problems, known as inflation. This will motivate work done to search for evidence of eternal inflation, which we describe in Chapter 7. We will then describe how density perturbations grow in the universe and discuss specific astrophysical signals that we can use to learn about structure growth, namely the CMB, the Sunyaev- Zeldovich e↵ect, and the cosmic infrared background. We will end this chapter with an 1When assessing the importance of these objects in the universe it should be remembered that their exis- tence is not motivated by any observational evidence. 1.1. OUR COSMOLOGY 5
overview of the Planck experiment and briefly mention work we have done to mitigate systematic uncertainties in the Planck data.
1.1.1 The Friedmann Equations
Starting from the Einstein equations of general relativity we can derive the Friedmann equations governing the evolution of the scale factor. The Einstein equations are:
1 Rµ⌫ gµ⌫R + gµ⌫⇤ = 8⇡GTµ⌫ , (1.1) 2 where Tµ⌫ is the stress-energy tensor, Rµ⌫ is the Ricci curvature tensor, ⇤ is the cosmological constant, gµ⌫ is the metric tensor, and we are using units with c = 1. On large scales the universe is found to be homogeneous and isotropic, and the most general metric satisfying these properties is the Friedmann-Robertson-Walker (FRW) metric:
2 2 µ ⌫ 2 2 dr 2 2 2 2 ds = gµ⌫dx dx = dt + a (t) + r (d✓ + sin ✓d ) , (1.2) 1 kr2 " # where a is the scale factor, r is a distance coordinate, and k is the spatial curvature with values of -1, 0, or +1 for an open, flat, or closed universe. To proceed further we require information about the energy density of the universe, which can be encoded in a stress- energy tensor. The di↵erent components of the stress-energy tensor are the energy density
T00, energy flux density T0i, momentum density Ti0, and the spatial stress tensor Tij. For a perfect fluid with no shear stresses, viscosity, or heat conduction, that is in thermodynamic equilibrium the stress-energy tensor is:
Tµ⌫ = (⇢ + p) uµu⌫ + pgµ⌫ , (1.3)
where ⇢ is the energy density, p is the hydrostatic pressure, and uµ is the velocity of the fluid, with gµ⌫u u = 1. If the fluid is at rest then the only non-zero entries are T = ⇢, µ ⌫ 00 and Tii = p. Solving the Einstein equations with the FRW metric gives the Friedmann Equations: a¨ 4⇡G ⇤ = (⇢ + 3p) + (1.4) a 3 3 6 CHAPTER 1. INTRODUCTION
a˙ 2 8⇡G⇢ + ⇤ k = . (1.5) a 3 a2 ✓ ◆ The dimensionless density parameters are defined by:
8⇡G⇢ ⇢ ⌦ = = , 2 (1.6) 3H ⇢c where H = a˙/a is the Hubble constant, and ⇢c is the critical density.
1.2 The Early Universe
1.2.1 Inflation
Inflation provides a solution to all three of the horizon, flatness, and defect problems. In the 1970s Guth discovered inflation while trying to understand why magnetic monopoles are not observed today [10]. He realized that an expansion of space resulting from a constant energy density would also resolve the horizon and flatness problems [11]. Inflation would cause the universe to rapidly expand making regions causally connected before expansion appear causally disconnected today and, depending on the number of e-folds (multiplica- tion by e), driving the universe to flatness. This expansion is thought to occur between 36 32 10 and 10 seconds after the Big Bang. After this period of rapid expansion ⇠ ⇠ the universe continues expanding but at a slower rate. Since the universe is known to be expanding now such a period of expansion is attractive theoretically since the two mecha- nisms could be related [e.g. 12]. In addition, inflation can provide a mechanism for generating the observed density fluc- tuations in the universe today. Any particle that causes inflation (known as the inflaton) would undergo quantum fluctuations in density, which would be magnified in size during expansion. Since space expands rapidly the fluctuations cannot be “undone”, as would be the case in a non-inflating universe, and they become the seeds for structure to grow. In a similar way quantum fluctuations in the graviton (the hypothetical particle that transmits the gravitational force) would lead to gravitational waves in the universe today that can potentially be detected. We note that inflation is only a hypothesis, and that models of the universe can be created to overcome the horizon problem without invoking inflation. 1.2. THE EARLY UNIVERSE 7
From equation 1.6 we see that the e↵ective curvature density is ⌦ = ⌦ 1 = /H2a2. k tot The flatness problem can therefore be solved if:
d 1 < 0 , (1.7) dt aH ! and the curvature density would decrease with time until it is close enough to zero that deviations from flatness after inflation would appear at much later times. Equation 1.7 also solves the horizon problem. The size of the comoving horizon (the radius of the region causally connected to an observer2) is:
a da0 1 ⌘ = . (1.8) a a H(a ) Z0 0 0 0 If 1/aH decreases with time then the size of the comoving horizon will shrink and re- gions previously in causal contact will become disconnected. If the Hubble constant is positive and approximately constant during inflation, then since H = a˙/a we find a(t) = a exp (H(t t )), where a and t are the values at the end of inflation. 0 0 0 0 One hypothesis is that inflation is caused by a scalar field, (t). For inflation to generate accelerating expansion d2a/dt2 must be positive, and so from the Friedmann equations p < ⇢/3. Since the energy density is positive the pressure is negative. Following Dodelson [14] we can find the pressure and energy density of the scalar field from the stress-energy tensor:
1 2 ⇢ = ˙ + V( ) 2 (1.9) 1 p = ˙ V( ) , 2 where V( ) is the potential for the field. From these equations we see how the pressure can be negative: if the field has more potential energy than kinetic energy. If the scalar field is in a local minimum of the potential but not a global minimum (known as a false vacuum) then it will have little kinetic energy but non-zero potential energy, and hence negative pressure.
2The number of “universe sized” patches with the same volume as ours must be at least 21 (95% confi- dence) based on current flatness limits Scott and Zibin [13]. 8 CHAPTER 1. INTRODUCTION
The Friedmann equations for the scalar field are:
dV ¨ + 3H ˙ = d (1.10) 8⇡G 1 H2 = ˙2 + V( ) . 3 2 ! Inflation ends when the scalar field rolls down the potential and the potential energy of the field is converted into radiation and particles in a process known as “reheating”. The universe then enters a period where the energy density is dominated by radiation.
1.2.2 Eternal Inflation
Quantum fluctuations in the inflaton field give some regions of space a larger inflaton po- tential. From equation 1.10 we see that regions with larger potential have greater expansion rates, and so at a later time they will occupy more space than regions with a lower potential. Over time these regions dominate the volume, and so if the decay rate to non-inflating uni- verses is su ciently low then the universe will end up “eternally inflating”. Regions with lower vacuum energy density appear as bubble universes. It is therefore possible that our universe is a bubble with a relatively low vacuum energy density, and hence low expansion rate, compared with the surrounding universe. An alternative possibility is that the inflaton potential could contain many false vacuum states. In standard inflation expansion occurs while the universe is in a false vacuum state and ends when there is a phase transition to the true vacuum. If there are many false vacua then the inflaton in each region of space will be in one of many possible false vacuum states and can undergo a phase transition to a region with a lower vacuum energy density and di↵erent physical constants. The universe would then have a fractal like structure with bubbles in one vacuum state containing bubbles with lower vacuum energy densities. That the universe could contain a large number of false vacua is a prediction of string theory [e.g. 15]. In Fig. 1.1 we show a simulation of the spatial distribution of bubbles at a late time in an eternally inflating false vacuum from Kleban [16]. At late times the physical volume of the false vacuum is larger, resulting in a larger number of small bubbles. It is di cult to test this hypothesis experimentally since we can only observe events 1.2. THE EARLY UNIVERSE 9
Figure 1.1 Simulation of the spatial distribution of bubbles at a late time in an eternally inflating false vacuum. From [16]. within our horizon. However, if other bubble universes have collided with our own uni- verse, and the collision boundary is within our horizon (or has passed entirely through it) then there would be observational e↵ects that we could detect. There are several ways to search for a collision:
A collision leaves an imprint on the temperature and polarization pattern of the CMB • if the collision intersects the last scattering surface. We discuss this possibility further in Chapter 7, where we present results from a search for this signal.
If the collision wall has passed entirely through the last scattering surface then there • is a dipole induced in the CMB. A dipole asymmetry would be di cult to detect 10 CHAPTER 1. INTRODUCTION
above the motion of the Earth and Solar System, and depending on its amplitude may be impossible to measure.
In principle any measurement of large scale density perturbations is sensitive to the • e↵ects of a bubble collision. For example a collision would introduce a coherent flow of objects.
In addition to the observational challenges associated with providing evidence for eter- nal inflation, falsifying the theory would be di cult since a detection of any of these signa- tures could be explained by other phenomena. However, the polarization pattern introduced during a collision has a signature that would be di cult to explain in other models [17].
1.3 Growth of Structure
After inflation density perturbations begin to grow due to gravity. We now overview the equations governing the linear growth structure so that we can connect our observations with underlying physical processes. While we do not know the details of the evolution of the universe immediately after inflation (or even if inflation occurred) we know that at some point the universe had cooled su ciently to form an ionized plasma of protons and electrons. Coulomb scattering strongly couples the electrons and protons in the plasma and Compton scattering couples the electrons and photons, with all three components a↵ected by gravity. Over-dense regions of the plasma therefore collapse under gravity, resisted by the pressure of the fluid. The relative strength of these two forces at a given time deter- mines whether the fluid collapses or expands, and it undergoes acoustic oscillations with the compressed regions becoming hotter and the rarefied regions colder. The universe had cooled su ciently 380,000 years after the Big Bang for electrons and protons to condense into atoms, and as the number of free electrons becomes smaller the Coulomb length gets larger, decoupling the photons from the electrons and protons. Once the Coulomb length is greater than the Hubble length photons travel through the universe without scattering, and the baryons fall into the dark matter potential wells. Dark matter and neutrinos do not interact electromagnetically (with neutrino decoupling occurring 1 second after the Big ⇠ Bang), and their dynamics are only a↵ected by gravity. 1.3. GROWTH OF STRUCTURE 11
For small density fluctuations, = ⇢/⇢ 1, the evolution of each component is accu- ⌧ rately described by perturbing the Einstein equations and neglecting second order terms. At 1 the present time, linear perturbation theory applies on scales larger than k & 10 Mpc. On smaller scales gravitational collapse must be treated using a non-linear theory, for example with the halo model, which we use in Chapters 3 and 4. We will only state the main results here since a full description requires lengthy treatment. We closely follow Dodelson [14], which should be referred to for more details.
For each component of the fluid a Boltzmann equation can be written that describes its evolution from the initial inhomogeneities and accounts for its interaction with all of the other components. The Boltzmann equation describes the time evolution of the phase space density ⇢(x, p, t), where x and p are position and momentum. It formalizes the statement that the rate of change of the number of particles is equal to the di↵erence between the rate of production and the rate of destruction. While analytical solutions to the equations can be found in limiting cases when either a single component dominates or on large or small scales, on intermediate scales when the matter and radiation energy densities are similar numerical solutions must be found. This can be done using a numerical code such as CAMB [18].
We start by writing down the perturbed metric. While the unperturbed metric depends only on the scale factor a(t), the perturbed metric depends on two additional functions: and :
g00(x, t) = 1 2 (x, t) g0i(x, t) = 0 (1.11) 2 gij(x, t) = a ij (1 + 2 (x, t)) , where is the perturbation to the spatial curvature and is the Newtonian potential. For linear theory it is easiest to work in Fourier space where the linearity ensures that each Fourier mode of the density and velocity fields evolves independently. Solving the Boltz- mann equations to first order in and v and assuming that the velocity field is irrotational 12 CHAPTER 1. INTRODUCTION
we find: