Cosmology with Large-scale Structure of the Universe Eiichiro - - PowerPoint PPT Presentation

cosmology with large scale structure of the universe
SMART_READER_LITE
LIVE PREVIEW

Cosmology with Large-scale Structure of the Universe Eiichiro - - PowerPoint PPT Presentation

Cosmology with Large-scale Structure of the Universe Eiichiro Komatsu (Texas Cosmology Center, UT Austin) Korean Young Cosmologists Workshop, June 27, 2011 Cosmology Update: WMAP 7-year+ Standard Model H&He = 4.58% (0.16%)


slide-1
SLIDE 1

Cosmology with Large-scale Structure of the Universe

Eiichiro Komatsu (Texas Cosmology Center, UT Austin) Korean Young Cosmologists Workshop, June 27, 2011

slide-2
SLIDE 2

Cosmology Update: WMAP 7-year+

  • Standard Model
  • H&He = 4.58% (±0.16%)
  • Dark Matter = 22.9% (±1.5%)
  • Dark Energy = 72.5% (±1.6%)
  • H0=70.2±1.4 km/s/Mpc
  • Age of the Universe = 13.76 billion

years (±0.11 billion years)

“ScienceNews” article on the WMAP 7-year results

2

slide-3
SLIDE 3

Cosmology: Next Decade?

  • Astro2010: Astronomy & Astrophysics Decadal Survey
  • Report from Cosmology and Fundamental Physics Panel

(Panel Report, Page T

  • 3):

3

slide-4
SLIDE 4

Cosmology: Next Decade?

  • Astro2010: Astronomy & Astrophysics Decadal Survey
  • Report from Cosmology and Fundamental Physics Panel

(Panel Report, Page T

  • 3): Translation

Inflation Dark Energy Dark Matter Neutrino Mass

4

slide-5
SLIDE 5

Cosmology: Next Decade?

  • Astro2010: Astronomy & Astrophysics Decadal Survey
  • Report from Cosmology and Fundamental Physics Panel

(Panel Report, Page T

  • 3): Translation

Inflation Dark Energy Dark Matter Neutrino Mass

5

Large-scale structure of the universe has a potential to give us valuable information on all of these items.

slide-6
SLIDE 6

What to measure?

  • Inflation
  • Shape of the initial power spectrum (ns; dns/dlnk; etc)
  • Non-Gaussianity (3pt fNLlocal; 4pt τNLlocal; etc)
  • Dark Energy
  • Angular diameter distances over a wide redshift range
  • Hubble expansion rates over a wide redshift range
  • Growth of linear density fluctuations over a wide

redshift range

  • Shape of the matter power spectrum (modified grav)

6

slide-7
SLIDE 7

What to measure?

  • Neutrino Mass
  • Shape of the matter power spectrum
  • Dark Matter
  • Shape of the matter power spectrum (warm/hot DM)
  • Large-scale structure traced by γ-ray photons

7

slide-8
SLIDE 8

Shape of the Power Spectrum, P(k)

Hlozek et al., arXiv:1105.4887 non-linear P(k) at z=0 linear P(k)

Matter density fluctuations measured by various tracers, extrapolated to z=0 CMB, z=1090 (l=2–3000) Galaxy, z=0.3 Gas, z=3

8

slide-9
SLIDE 9

Shape of the Power Spectrum, P(k)

non-linear P(k) at z=0 linear P(k)

Matter density fluctuations measured by various tracers, extrapolated to z=0 CMB, z=1090 (l=2–3000) Galaxy, z=0.3 Gas, z=3

Primordial spectrum, Pprim(k)~kns

9

slide-10
SLIDE 10

non-linear P(k) at z=0 linear P(k) asymptotes to kns(lnk)2/k4

T(k): Suppression of power during the radiation- dominated era.

Primordial spectrum, Pprim(k)~kns

The suppression depends

  • n Ωcdmh2 and Ωbaryonh2

P(k)=A x kns x T2(k)

10

slide-11
SLIDE 11

Current Limit on ns

  • Planck’s CMB data are expected to improve the error

bar by a factor of ~4.

  • Limit on the tilt of the power spectrum:
  • ns=0.968±0.012 (68%CL; Komatsu et al. 2011)
  • Precision is dominated by the WMAP 7-year data

11

slide-12
SLIDE 12

Probing Inflation (2-point Function)

  • Joint constraint on the

primordial tilt, ns, and the tensor-to-scalar ratio, r.

  • Not so different from the

5-year limit.

  • r < 0.24 (95%CL)
  • Limit on the tilt of the

power spectrum: ns=0.968±0.012 (68%CL)

12

Komatsu et al. (2011)

r = (gravitational waves)2 / (gravitational potential)2 Planck?

slide-13
SLIDE 13

Role of the Large-scale Structure of the Universe

  • However, CMB data can’t go much beyond k=0.2 Mpc–1

(l=3000).

  • Large-scale structure data are required to go to

smaller scales.

13

slide-14
SLIDE 14

Shape of the Power Spectrum, P(k)

non-linear P(k) at z=0 linear P(k)

Matter density fluctuations measured by various tracers, extrapolated to z=0 CMB, z=1090 (l=2–3000) Galaxy, high-z Gas, z=3

14

slide-15
SLIDE 15

Measuring a scale- dependence of ns(k)

  • As far as the value of ns is concerned, CMB is probably

enough.

  • However, if we want to measure the scale-dependence of

ns, i.e., deviation of Pprim(k) from a pure power-law, then we need the small-scale data.

  • This is where the large-scale structure data become

quite powerful (Takada, Komatsu & Futamase 2006)

  • Schematically:
  • dns/dlnk = [ns(CMB) - ns(LSS)]/(lnkCMB - lnkLSS)

15

slide-16
SLIDE 16

Probing Inflation (3-point Function)

  • Inflation models predict that primordial fluctuations are very

close to Gaussian.

  • In fact, ALL SINGLE-FIELD models predict a particular form
  • f 3-point function to have the amplitude of fNLlocal=0.02.
  • Detection of fNL>1 would rule out ALL single-field models!

16

Can We Rule Out Inflation?

slide-17
SLIDE 17

Bispectrum

  • Three-point function!
  • Bζ(k1,k2,k3)

= <ζk1ζk2ζk3> = (amplitude) x (2π)3δ(k1+k2+k3)F(k1,k2,k3)

17

model-dependent function

k1 k2 k3 Primordial fluctuation

slide-18
SLIDE 18

MOST IMPORTANT

slide-19
SLIDE 19

Single-field Theorem (Consistency Relation)

  • For ANY single-field models*, the bispectrum in the

squeezed limit is given by

  • Bζ(k1~k2<<k3)≈(1–ns) x (2π)3δ(k1+k2+k3) x Pζ(k1)Pζ(k3)
  • Therefore, all single-field models predict fNL≈(5/12)(1–ns).
  • With the current limit ns=0.968, fNL is predicted to be 0.01.

Maldacena (2003); Seery & Lidsey (2005); Creminelli & Zaldarriaga (2004)

* for which the single field is solely responsible for driving inflation and generating observed fluctuations.

19

slide-20
SLIDE 20

Probing Inflation (3-point Function)

  • No detection of 3-point functions of primordial curvature
  • perturbations. The 95% CL limit is:
  • –10 < fNLlocal < 74
  • The 68% CL limit: fNLlocal = 32 ± 21
  • The WMAP data are consistent with the prediction of

simple single-field inflation models: 1–ns≈r≈fNL

  • The Planck’s expected 68% CL uncertainty: ΔfNLlocal = 5

20

Komatsu et al. (2011)

slide-21
SLIDE 21

Trispectrum

  • Tζ(k1,k2,k3,k4)=(2π)3δ(k1+k2+k3+k4) {gNL[(54/25)

Pζ(k1)Pζ(k2)Pζ(k3)+cyc.] +τNL[Pζ(k1)Pζ(k2)(Pζ(| k1+k3|)+Pζ(|k1+k4|))+cyc.]} k3 k4 k2 k1

gNL

k2 k1 k3 k4

τNL

21

slide-22
SLIDE 22

τNLlocal–fNLlocal Diagram

  • The current limits

from WMAP 7-year are consistent with single-field or multi- field models.

  • So, let’s play around

with the future.

22

ln(fNL) ln(τNL) 74 3.3x104

(Smidt et

  • al. 2010)

x0.5

slide-23
SLIDE 23

Case A: Single-field Happiness

  • No detection of

anything after

  • Planck. Single-field

survived the test (for the moment: the future galaxy surveys can improve the limits by a factor of ten). ln(fNL) ln(τNL) 10 600

23

x0.5

slide-24
SLIDE 24

Case B: Multi-field Happiness

  • fNL is detected. Single-

field is dead.

  • But, τNL is also

detected, in accordance with multi- field models: τNL>0.5 (6fNL/5)2 [Sugiyama, Komatsu & Futamase (2011)] ln(fNL) ln(τNL) 600

24

30 x0.5

slide-25
SLIDE 25

Case C: Madness

  • fNL is detected. Single-

field is dead.

  • But, τNL is not

detected, inconsistent with the multi-field bound.

  • (With the caveat that

this bound may not be completely general) BOTH the single-field and multi-field are gone. ln(fNL) ln(τNL) 30 600

25

x0.5

slide-26
SLIDE 26

Beyond CMB: Large-scale Structure!

  • In principle, the large-scale structure of the universe
  • ffers a lot more statistical power, because we can get

3D information. (CMB is 2D, so the number of Fourier modes is limited.)

26

slide-27
SLIDE 27

Beyond CMB: Large-scale Structure?

  • Statistics is great, but the large-scale structure is non-

linear, so perhaps it is less clean?

  • Not necessarily.

27

slide-28
SLIDE 28

MOST IMPORTANT

slide-29
SLIDE 29

Non-linear Gravity

  • For a given k1, vary k2 and k3, with k3≤k2≤k1
  • F2(k2,k3) vanishes in the squeezed limit, and peaks at the

elongated triangles.

29

slide-30
SLIDE 30

Non-linear Galaxy Bias

  • There is no F2: less suppression at the squeezed, and

less enhancement along the elongated triangles.

  • Still peaks at the equilateral or elongated forms.

30

slide-31
SLIDE 31

Primordial Non-Gaussianity

  • This gives the peaks at the squeezed configurations,

clearly distinguishable from other non-linear/ astrophysical effects.

31

Sefusatti & Komatsu (2007); Jeong & Komatsu (2010)

slide-32
SLIDE 32

Bispectrum is powerful

  • fNLlocal ~ O(1) is quite possible with the bispectrum
  • method. (See Donghui Jeong’s talk)
  • This needs to be demonstrated by the real data! (e.g.,

SDSS-LRG)

32

slide-33
SLIDE 33

Need For Dark “Energy”

  • First of all, DE does not even need to be an energy.
  • At present, anything that can explain the observed

(1) Luminosity Distances (Type Ia supernovae) (2) Angular Diameter Distances (BAO, CMB) simultaneously is qualified for being called “Dark Energy.”

  • The candidates in the literature include: (a) energy, (b)

modified gravity, and (c) extreme inhomogeneity.

  • Measurements of the (3) growth of structure break
  • degeneracy. (The best data right now is the X-ray clusters.)

33

slide-34
SLIDE 34

H(z): Current Knowledge

  • H2(z) = H2(0)[Ωr(1+z)4+Ωm(1+z)3+Ωk(1+z)2+Ωde(1+z)3(1+w)]
  • (expansion rate) H(0) = 70.2 ± 1.4 km/s/Mpc
  • (radiation) Ωr = (8.4±0.3)x10-5
  • (matter) Ωm = 0.275±0.016
  • (curvature) Ωk < 0.008 (95%CL)
  • (dark energy) Ωde = 0.725±0.015
  • (DE equation of state) w = –1.00±0.06

WMAP7+

34

slide-35
SLIDE 35

H(z) to Distances

  • Comoving Distance
  • χ(z) = c∫z[dz’/H(z’)]
  • Luminosity Distance
  • DL(z) = (1+z)χ(z)[1–(k/6)χ2(z)/R2+...]
  • R=(curvature radius of the universe); k=(sign of curvature)
  • WMAP 7-year limit: R>2χ(∞); justify the Taylor expansion
  • Angular Diameter Distance
  • DA(z) = [χ(z)/(1+z)][1–(k/6)χ2(z)/R2+...]

35

slide-36
SLIDE 36

DA(z) = (1+z)–2 DL(z)

  • To measure DA(z), we need to know the intrinsic size.
  • What can we use as the standard ruler?

Redshift, z

0.2 2 6 1090

Type 1a Supernovae Galaxies (BAO) CMB

DL(z) DA(z)

0.02

36

slide-37
SLIDE 37

How Do We Measure DA(z)?

  • If we know the intrinsic physical sizes, d, we can

measure DA. What determines d?

Redshift, z

0.2 2 6 1090

Galaxies CMB

0.02

DA(galaxies)=dBAO/θ

dBAO dCMB

DA(CMB)=dCMB/θ

θ θ

37

slide-38
SLIDE 38

CMB as a Standard Ruler

  • The existence of typical spot size in image space yields
  • scillations in harmonic (Fourier) space. What

determines the physical size of typical spots, dCMB?

θ θ~the typical size of hot/cold spots θ θ θ θ θ θ θ

38

slide-39
SLIDE 39

Sound Horizon

  • The typical spot size, dCMB, is determined by the

physical distance traveled by the sound wave from the Big Bang to the decoupling of photons at zCMB~1090 (tCMB~380,000 years).

  • The causal horizon (photon horizon) at tCMB is given by
  • dH(tCMB) = a(tCMB)*Integrate[ c dt/a(t), {t,0,tCMB}].
  • The sound horizon at tCMB is given by
  • ds(tCMB) = a(tCMB)*Integrate[ cs(t) dt/a(t), {t,0,tCMB}],

where cs(t) is the time-dependent speed of sound

  • f photon-baryon fluid.

39

slide-40
SLIDE 40
  • The WMAP 7-year values:
  • lCMB = π/θ = πDA(zCMB)/ds(zCMB) = 302.69±0.76
  • CMB data constrain the ratio, DA(zCMB)/ds(zCMB).
  • rs(zCMB)=(1+zCMB)ds(zCMB)=146.6±1.6 Mpc (comoving)

lCMB=302.69±0.76

40

slide-41
SLIDE 41
  • Color: constraint from

lCMB=πDA(zCMB)/ds(zCMB) with zEQ & Ωbh2.

  • Black contours: Markov

Chain from WMAP 3yr (Spergel et al. 2007)

What DA(zCMB)/ds(zCMB) Gives You (3-year example)

lCMB=301.8±1.2

1-Ωm-ΩΛ = 0.3040Ωm +0.4067ΩΛ

41

slide-42
SLIDE 42

0.0 0.5 1.0 1.5 2.0 M 0.0 0.5 1.0 1.5 2.0

  • ESSENCE+SNLS+gold

(M,) = (0.27,0.73) Total=1

42

slide-43
SLIDE 43

BAO in Galaxy Distribution

  • The same acoustic oscillations should be hidden in this

galaxy distribution... 2dFGRS

43

slide-44
SLIDE 44

BAO as a Standard Ruler

  • The existence of a localized clustering scale in the 2-point

function yields oscillations in Fourier space. (1+z)dBAO Percival et al. (2006) Okumura et al. (2007)

Position Space Fourier Space

44

slide-45
SLIDE 45

Sound Horizon Again

  • The clustering scale, dBAO, is given by the physical distance

traveled by the sound wave from the Big Bang to the decoupling of baryons at zBAO=1020.5±1.6 (c.f., zCMB=1091±1).

  • The baryons decoupled slightly later than CMB.
  • By the way, this is not universal in cosmology, but

accidentally happens to be the case for our Universe.

  • If 3ρbaryon/(4ρphoton) =0.64(Ωbh2/0.022)(1090/(1+zCMB)) is

greater than unity, zBAO>zCMB. Since our Universe happens to have Ωbh2=0.022, zBAO<zCMB. (ie, dBAO>dCMB)

45

slide-46
SLIDE 46

Standard Rulers in CMB & Matter

  • For flat LCDM, but very similar results for w≠–1 and

curvature≠0!

46

Komatsu et al. (2009)

slide-47
SLIDE 47

Not Just DA(z)...

  • A really nice thing about BAO at a given redshift is that

it can be used to measure not only DA(z), but also the expansion rate, H(z), directly, at that redshift.

  • BAO perpendicular to l.o.s

=> DA(z) = ds(zBAO)/θ

  • BAO parallel to l.o.s

=> H(z) = cΔz/[(1+z)ds(zBAO)]

47

slide-48
SLIDE 48

Transverse=DA(z); Radial=H(z)

Two-point correlation function measured from the SDSS Luminous Red Galaxies (Gaztanaga, Cabre & Hui 2008) (1+z)ds(zBAO)

θ = ds(zBAO)/DA(z) cΔz/(1+z) = ds(zBAO)H(z)

Linear Theory SDSS Data

48

slide-49
SLIDE 49

Percival et al. (2010)

Redshift, z

2dFGRS and SDSS main samples SDSS LRG samples

(1+zBAO)ds(zBAO)/DV(z)

Ωm=0.278, ΩΛ=0.722

49

0.2 0.3 0.4

DV(z) = {(1+z)2DA2(z)[cz/H(z)]}1/3

Since the current data are not good enough to constrain DA(z) and H(z) separately, a combination distance, DV(z), has been constrained.

slide-50
SLIDE 50

WMAP7+BAO+...

  • At the moment, BAO is

great for fixing curvature, but not good for fixing w

  • We still need supernovae

for fixing w, but this would change as more BAO data (especially at higher redshifts) become available.

50

Komatsu et al. (2011)

slide-51
SLIDE 51

w(z)=w0+wa*z/(1+z)

Komatsu et al. (2011)

  • Cosmological constant,

w0=–1 and wa=0, are perfectly consistent with data.

  • Of course we all want this

to change at some point...

51

slide-52
SLIDE 52

Hobby-Eberly Telescope Dark Energy Experiment (HETDEX)

52

Use 9.2-m HET to map the universe using 0.8M Lyman-alpha emitting galaxies in z=1.9–3.5

slide-53
SLIDE 53

HETDEX Foot-print (in RA-DEC coordinates)

53

slide-54
SLIDE 54

HETDEX: Sound Waves in the Distribution of Galaxies

  • 1000
  • 500

500 1000

  • 1000
  • 500

500 1000

Sloan Digital Sky Survey

54

Small Scale Large Scale

slide-55
SLIDE 55

HETDEX: Sound Waves in the Distribution of Galaxies

  • 1000
  • 500

500 1000

  • 1000
  • 500

500 1000

HETDEX

HETDEX vs SDSS

10x more galaxies observed 3x larger volume surveyed Will survey the previously unexplored discovery space

55

Small Scale Large Scale

slide-56
SLIDE 56

Beyond BAO

  • BAOs capture only a fraction of the information

contained in the galaxy power spectrum!

  • The full usage of the 2-dimensional power spectrum

leads to a substantial improvement in the precision of distance and expansion rate measurements.

56

slide-57
SLIDE 57

BAO vs Full Modeling

  • Full modeling improves upon

the determinations of DA & H by more than a factor of two.

  • On the DA-H plane, the size
  • f the ellipse shrinks by more

than a factor of four. Shoji, Jeong & Komatsu (2008)

57

slide-58
SLIDE 58

Alcock-Paczynski: The Most Important Thing For HETDEX

  • Where does the improvement

come from?

  • The Alcock-Paczynski test is the key.

This is the most important component for the success of the HETDEX survey.

58

slide-59
SLIDE 59

The AP Test: How That Works

  • The key idea: (in the absence of the redshift-space

distortion - we will include this for the full analysis; we ignore it here for simplicity), the distribution of the power should be isotropic in Fourier space.

59

slide-60
SLIDE 60
  • DA: (RA,Dec) to the transverse separation, rperp, to the

transverse wavenumber

  • kperp = (2π)/rperp = (2π)[Angle on the sky]/DA
  • H: redshifts to the parallel separation, rpara, to the

parallel wavenumber

  • kpara = (2π)/rpara = (2π)H/(cΔz)

The AP Test: How That Works

If DA and H are correct: kpara kperp If DA is wrong: kperp If H is wrong: kperp

60

slide-61
SLIDE 61
  • DA: (RA,Dec) to the transverse separation, rperp, to the

transverse wavenumber

  • kperp = (2π)/rperp = (2π)[Angle on the sky]/DA
  • H: redshifts to the parallel separation, rpara, to the

parallel wavenumber

  • kpara = (2π)/rpara = (2π)H/(cΔz)

The AP Test: How That Works

If DA and H are correct: kpara kperp If DA is wrong: kperp If H is wrong: kperp kperp If DA and H are wrong:

slide-62
SLIDE 62

DAH from the AP test

  • So, the AP test can’t be used

to determine DA and H separately; however, it gives a measurement of DAH.

  • Combining this with the BAO

information, and marginalizing

  • ver the redshift space

distortion, we get the solid contours in the figure.

62

slide-63
SLIDE 63

Redshift Space Distortion

  • Both the AP test and the redshift space distortion make

the distribution of the power anisotropic. Would it spoil the utility of this method?

  • Some, but not all!

63

f is marginalized over. f is fixed.

slide-64
SLIDE 64

WMAP Amplitude Prior

  • WMAP measures the amplitude of curvature

perturbations at z~1090. Let’s call that Rk. The relation to the density fluctuation is

  • Variance of Rk has been constrained as:

64

where kWMAP=0.027 Mpc–1

slide-65
SLIDE 65

Then Solve This Diff. Equation...

g(z)=(1+z)D(z)

65

w w

slide-66
SLIDE 66

Degeneracy Between Amplitude at z=0 (σ8) and w

Flat Universe Non-flat Univ.

66

slide-67
SLIDE 67

Alexey Vikhlinin, from a slide presented at the IPMU Dark Energy Conference in Japan, June 2009

g(z)=

67

slide-68
SLIDE 68

HETDEX and Neutrino Mass

  • Neutrinos suppress

the matter power spectrum on small scales (k>0.1 h Mpc–1).

  • A useful number to

remember:

  • For ∑mν=0.1 eV, the

power spectrum at k>0.1 h Mpc–1 is suppressed by ~7%.

  • We can measure this

easily!

For 10x the number density of HETDEX

68

slide-69
SLIDE 69

Neutrino Mass and P(k)

  • Total neutrino mass: coming from the small scale
  • ΔP/P ~ –8Ων/Ωm = –[8/(Ωmh2)]∑mν/(
  • Where the suppression begins depends on individual

masses!

  • For 10x the number density of HETDEX

69

slide-70
SLIDE 70

Expectation for HETDEX

  • CV limited: error goes as 1/sqrt(volume)
  • SN limited: error goes as 1/(number density)/sqrt(volume)

cosmic variance limited regime shot noise limited regime

70

slide-71
SLIDE 71

Expected HETDEX Limit

  • ~6x better than WMAP 7-year+H0

71

slide-72
SLIDE 72

Summary

  • Three (out of four) questions:
  • What is the physics of inflation?
  • P(k) shape (esp, dn/dlnk) and non-Gaussianity
  • What is the nature of dark energy?
  • DA(z), H(z), growth of structure
  • What is the mass of neutrinos?
  • P(k) shape
  • CMB and large-scale structure observations can lead to

major breakthroughs in any of the above questions.

72