Lecture 3 - Temperature anisotropy from sound waves (continued) - - - PowerPoint PPT Presentation

lecture 3
SMART_READER_LITE
LIVE PREVIEW

Lecture 3 - Temperature anisotropy from sound waves (continued) - - - PowerPoint PPT Presentation

Lecture 3 - Temperature anisotropy from sound waves (continued) - Cosmological parameter dependence of the temperature power spectrum Stone: Fluctuations entering the horizon This is a tricky concept, but it is important Suppose


slide-1
SLIDE 1

Lecture 3

  • Temperature anisotropy from sound waves

(continued)

  • Cosmological parameter dependence of the

temperature power spectrum

slide-2
SLIDE 2
slide-3
SLIDE 3
slide-4
SLIDE 4

Stone: Fluctuations “entering the horizon”

  • This is a tricky concept, but it is important
  • Suppose that there are fluctuations at all wavelengths,

including the ones that exceed the Hubble length (which we loosely call our “horizon”)

  • Let’s not ask the origin of these “super-horizon

fluctuations”, but just assume their existence

  • As the Universe expands, our horizon grows and we can see

longer and longer wavelengths

  • Fluctuations “entering the horizon”
slide-5
SLIDE 5

10 Gpc/h today 1 Gpc/h today 100 Mpc/h today 10 Mpc/h today 1 Mpc/h today “enter the horizon” Radiation Era Last scattering Matter Era

slide-6
SLIDE 6

Three Regimes

  • Super-horizon scales [q < aH]
  • Only gravity is important
  • Evolution differs from Newtonian
  • Sub-horizon but super-sound-horizon [aH < q < aH/cs]
  • Only gravity is important
  • Evolution similar to Newtonian
  • Sub-sound-horizon scales [q > aH/cs]
  • Hydrodynamics important -> Sound waves
slide-7
SLIDE 7

qEQ

  • Which fluctuation entered the horizon before the matter-

radiation equality?

  • qEQ = aEQHEQ ~ 0.01 (ΩMh2/0.14) Mpc–1
  • At the last scattering surface, this subtends the multipole
  • f lEQ = qEQrL ~ 140
slide-8
SLIDE 8

Entered the horizon during the radiation era

slide-9
SLIDE 9

What determines the locations and heights of the peaks? Does the sound-wave solution explain it?

slide-10
SLIDE 10

Peak Locations?

  • VERY roughly speaking, the angular power spectrum Cl is given

by [ ]2 with q -> l/rL

  • Question: What are the integration constants, A and B?
  • Answer: They depend on the initial conditions; namely,

adiabatic or not?

  • For adiabatic initial condition, A >> B when q is large

High-frequency solution, for q >> aH

[We will show it later.]

slide-11
SLIDE 11

Peak Locations?

  • VERY roughly speaking, the angular power spectrum Cl is given

by [ ]2 with q -> l/rL

  • If A>>B, the locations of peaks are

High-frequency solution, for q >> aH

slide-12
SLIDE 12
slide-13
SLIDE 13

The simple estimates do not match! This is simply because these angular scales do not satisfy q >> aH, i.e, the

  • scillations are not pure

cosine even for the adiabatic initial condition. We need a better solution!

slide-14
SLIDE 14

Better Solution in Radiation-dominated Era

  • In the radiation-dominated era, R << 1
  • Change the independent variable from the time (t) to

Going back to the original tight-coupling equation..

slide-15
SLIDE 15

Better Solution in Radiation-dominated Era

  • In the radiation-dominated era, R << 1
  • Change the independent variable from the time (t) to

Then the equation simplifies to

where

slide-16
SLIDE 16

Better Solution in Radiation-dominated Era

Then the equation simplifies to

where

The solution is

slide-17
SLIDE 17

Better Solution in Radiation-dominated Era

Then the equation simplifies to

where

The solution is

where

slide-18
SLIDE 18

Einstein’s Equations

  • Now we need to know Newton’s gravitational potential, φ,

and the scalar curvature perturbation, ψ.

  • Einstein’s equations - let’s look up any text books:
slide-19
SLIDE 19

Einstein’s Equations

  • Now we need to know Newton’s gravitational potential, φ,

and the scalar curvature perturbation, ψ.

  • Einstein’s equations - let’s look up any text books:
slide-20
SLIDE 20

Einstein’s Equations

  • Now we need to know Newton’s gravitational potential, φ,

and the scalar curvature perturbation, ψ.

  • Einstein’s equations - let’s look up any text books:

Will come back to this later. For now, let’s ignore any viscosity.

slide-21
SLIDE 21

Einstein’s Equations

  • Now we need to know Newton’s gravitational potential, φ,

and the scalar curvature perturbation, ψ.

  • Einstein’s equations - let’s look up any text books:

Will come back to this later. For now, let’s ignore any viscosity.

slide-22
SLIDE 22

Einstein’s Equations in Radiation-dominated Era

  • Now we need to know Newton’s gravitational potential, φ,

and the scalar curvature perturbation, ψ.

  • Einstein’s equations - let’s look up any text books:

“non-adiabatic” pressure

slide-23
SLIDE 23

Einstein’s Equations in Radiation-dominated Era

  • Now we need to know Newton’s gravitational potential, φ,

and the scalar curvature perturbation, ψ.

  • Einstein’s equations - let’s look up any text books:

“non-adiabatic” pressure We shall ignore this

slide-24
SLIDE 24

Solution (Adiabatic) in Radiation-dominated Era

  • Low-frequency limit (super-sound-horizon scales, qrs << 1)
  • ΦADI -> –2ζ/3 = constant
  • High-frequency limit (sub-sound-horizon scales, qrs >> 1)
  • ΦADI -> 2ζ

ADI

where

damp

Kodama & Sasaki (1986, 1987)

slide-25
SLIDE 25

Solution (Adiabatic) in Radiation-dominated Era

  • Low-frequency limit (super-sound-horizon scales, qrs << 1)
  • ΦADI -> –2ζ/3 = constant
  • High-frequency limit (sub-sound-horizon scales, qrs >> 1)
  • ΦADI -> 2ζ

ADI

where

damp

Poisson Equation & oscillation solution for radiation

slide-26
SLIDE 26

Solution (Adiabatic) in Radiation-dominated Era

  • Low-frequency limit (super-sound-horizon scales, qrs << 1)
  • ΦADI -> –2ζ/3 = constant
  • High-frequency limit (sub-sound-horizon scales, qrs >> 1)
  • ΦADI -> 2ζ

ADI

where

damp

slide-27
SLIDE 27

ζ: Conserved on large scales

  • For the adiabatic initial condition, there exists a useful quantity,

ζ, which remains constant on large scales

(super-horizon scales, q << aH) regardless of the contents of the Universe

  • ζ is conserved regardless of whether the Universe is

radiation-dominated, matter-dominated, or whatever

  • Energy conservation for q << aH:

Bardeen, Steinhardt & Turner (1983); Weinberg (2003); Lyth, Malik & Sasaki (2005)

slide-28
SLIDE 28

ζ: Conserved on large scales

  • If pressure is a function of the energy density only, i.e.,

Integrate

, then

integration constant Bardeen, Steinhardt & Turner (1983); Weinberg (2003); Lyth, Malik & Sasaki (2005)

slide-29
SLIDE 29

ζ: Conserved on large scales

  • If pressure is a function of the energy density only, i.e.,

, then

integration constant

For the adiabatic initial condition, all species share the same value of ζα, i.e., ζα=ζ

Bardeen, Steinhardt & Turner (1983); Weinberg (2003); Lyth, Malik & Sasaki (2005)

slide-30
SLIDE 30

Sound Wave Solution in the Radiation-dominated Era

The solution is

where Kodama & Sasaki (1986, 1987); Baumann, Green, Meyers & Wallisch (2016)

slide-31
SLIDE 31

Sound Wave Solution in the Radiation-dominated Era

The solution is

where i.e.,

ADI ADI

Kodama & Sasaki (1986, 1987); Baumann, Green, Meyers & Wallisch (2016)

slide-32
SLIDE 32

Sound Wave Solution in the Radiation-dominated Era

The adiabatic solution is with

Therefore, the solution is a pure cosine

  • nly in the high-frequency limit!

Kodama & Sasaki (1986, 1987); Baumann, Green, Meyers & Wallisch (2016)

slide-33
SLIDE 33
slide-34
SLIDE 34

Roles of viscosity

  • Neutrino viscosity
  • Modify potentials:
  • Photon viscosity
  • Viscous photon-baryon fluid: damping of sound waves

Silk (1968) “Silk damping”

slide-35
SLIDE 35

High-frequency solution without neutrino viscosity

The solution is

where

slide-36
SLIDE 36

High-frequency solution with neutrino viscosity

The solution is

where Chluba & Grin (2013) non-zero value!

slide-37
SLIDE 37

High-frequency solution with neutrino viscosity

The solution is

where Hu & Sugiyama (1996) Bashinsky & Seljak (2004) Phase shift!

slide-38
SLIDE 38

High-frequency solution with neutrino viscosity

The solution is

where Hu & Sugiyama (1996) Bashinsky & Seljak (2004) Phase shift!

Thus, the neutrino viscosity will: (1) Reduce the amplitude

  • f sound waves at large multipoles

(2) Shift the peak positions

  • f the temperature power spectrum
slide-39
SLIDE 39

Photon Viscosity

  • In the tight-coupling approximation, the photon viscosity

damps exponentially

  • To take into account a non-zero photon viscosity, we go

to a higher order in the tight-coupling approximation

slide-40
SLIDE 40

Tight-coupling Approximation (1st-order)

  • When Thomson scattering is efficient, the relative velocity

between photons and baryons is small. We write

[d is an arbitrary dimensionless variable]

  • And take *. We obtain

*In this limit, viscosity πγ is exponentially suppressed. This result comes from the Boltzmann equation but we do not derive it here. It makes sense physically.

slide-41
SLIDE 41

Tight-coupling Approximation (2nd-order)

  • When Thomson scattering is efficient, the relative velocity

between photons and baryons is small. We write

[d2 is an arbitrary dimensionless variables]

  • And take .. We obtain

where

slide-42
SLIDE 42

Tight-coupling Approximation (2nd-order)

  • Eliminating d2 and using the fact that R is proportional to

the scale factor, we obtain

  • Getting πγ requires an approximate solution of the Boltzmann

equation in the 2nd-order tight coupling. We do not derive it

  • here. The answer is

Kaiser (1983)

slide-43
SLIDE 43

Tight-coupling Approximation (2nd-order)

  • Eliminating d2 and using the fact that R is proportional to

the scale factor, we obtain

  • Getting πγ requires an approximate solution of the Boltzmann

equation in the 2nd-order tight coupling. We do not derive it

  • here. The answer is

Kaiser (1983) given by the velocity potential

  • a well-known result in fluid

dynamics

slide-44
SLIDE 44

Damped Oscillator

  • Using the energy conservation to replace δuγ with δργ/ργ,

we obtain, for q >> aH,

New term, giving damping! where

slide-45
SLIDE 45

Damped Oscillator

  • Using the energy conservation to replace δuγ with δργ/ργ,

we obtain, for q >> aH,

New term, giving damping! where Important for high frequencies (large multipoles)

slide-46
SLIDE 46

Damped Oscillator

  • Using the energy conservation to replace δuγ with δργ/ργ,

we obtain, for q >> aH,

New term, giving damping! Exponential dampling! SOLUTION:

slide-47
SLIDE 47

Damped Oscillator

  • Using the energy conservation to replace δuγ with δργ/ργ,

we obtain, for q >> aH,

New term, giving damping! Exponential dampling! SOLUTION: Silk Silk “diffusion length” = length traveled by photon’s random walks

slide-48
SLIDE 48

Planck Collaboration (2016)

Sachs-Wolfe Sound Wave Silk Damping?

slide-49
SLIDE 49

Additional Damping

Landau

( )

  • The power spectrum is [ ]2 with q -> l/rL. The damping factor

is thus exp(–2q2/qsilk2)

  • qsilk(tL) = 0.139 Mpc–1. This corresponds to a multipole of lsilk ~ qsilk

rL/√2 = 1370. Seems too large, compared to the exact calculation

  • There is an additional damping due to a finite width of the last

scattering surface, σ~250 K

  • “Fuzziness damping” – Bond (1996)
  • “Landau damping” - Weinberg (2001)
slide-50
SLIDE 50

Planck Collaboration (2016)

Sachs-Wolfe Sound Wave Silk+Landau Damping Total damping: qD–2 = qsilk–2 + qlandau–2 qD ~ 0.11 Mpc–1, giving lD ~ qDrL/√2 ~ 1125

slide-51
SLIDE 51

Planck Collaboration (2016)

Sachs-Wolfe Sound Wave Next: Let’s understand the peak heights Silk+Landau Damping

slide-52
SLIDE 52

Matching Solutions

  • We have a very good analytical solution valid at low and

high frequencies during the radiation era:

  • Now, match this to a high-frequency solution valid at the

last-scattering surface (when R is no longer small)

slide-53
SLIDE 53

Matching Solutions

  • We have a very good analytical solution valid at low and

high frequencies during the radiation era:

  • Now, match this to a high-frequency solution valid at the

last-scattering surface (when R is no longer small)

Slightly improved solution, with a weak time dependence of R using the WKB method [Peebles & Yu (1970)]

slide-54
SLIDE 54

High-frequency Solution(*) at the Last Scattering Surface

  • (*) To a good approximation, the low-frequency solution is

given by setting R=0 because sound waves are not important at large scales

Weinberg “Cosmology”, Eq. (6.5.7) where T(q), S(q), θ(q) are “transfer functions” that smoothly interpolate two limits as

q q q

q << qEQ: q >> qEQ:

with qEQ = aEQHEQ ~ 0.01 Mpc–1, giving lEQ=qEQrL ~ 140

“EQ” for “matter-radiation Equality epoch”

q q

slide-55
SLIDE 55

High-frequency Solution(*) at the Last Scattering Surface

  • (*) To a good approximation, the low-frequency solution is

given by setting R=0 because sound waves are not important at large scales

Weinberg “Cosmology”, Eq. (6.5.7) where T(q), S(q), θ(q) are “transfer functions” that smoothly interpolate two limits as

q q q

q << qEQ: q >> qEQ:

with qEQ = aEQHEQ ~ 0.01 Mpc–1, giving lEQ=qEQrL ~ 140

“EQ” for “matter-radiation Equality epoch”

Due to the decay of gravitational potential during the radiation dominated era

q q

slide-56
SLIDE 56

High-frequency Solution(*) at the Last Scattering Surface

  • (*) To a good approximation, the low-frequency solution is

given by setting R=0 because sound waves are not important at large scales

Weinberg “Cosmology”, Eq. (6.5.7) where T(q), S(q), θ(q) are “transfer functions” that smoothly interpolate two limits as

q q q

q << qEQ: q >> qEQ:

with qEQ = aEQHEQ ~ 0.01 Mpc–1, giving lEQ=qEQrL ~ 140

“EQ” for “matter-radiation Equality epoch”

Due to the neutrino anisotropic stress

q q

slide-57
SLIDE 57
  • (*) To a good approximation, the low-frequency solution is

given by setting R=0 because sound waves are not important at large scales

High-frequency Solution(*) at the Last Scattering Surface

Weinberg “Cosmology”, Eq. (6.5.7)

q q q

q -> 0(*)

−ζ 5

This should agree with the Sachs-Wolfe result: Φ/3; thus,

Φ = −3ζ/5 in the matter-dominated era

slide-58
SLIDE 58
  • (*) To a good approximation, the low-frequency solution is

given by setting R=0 because sound waves are not important at large scales

Effect of Baryons

Weinberg “Cosmology”, Eq. (6.5.7)

q q q

Shift the zero-point of

  • scillations

Reduce the amplitude of

  • scillations
slide-59
SLIDE 59

` ≈ 302 × qrs/⇡

No Baryon [R=0]

slide-60
SLIDE 60

` ≈ 302 × qrs/⇡

No Baryon [R=0]

B

  • s

t d u e t

  • d

e c a y i n g p

  • t

e n t i a l d u r i n g t h e r a d i a t i

  • n

e r a

slide-61
SLIDE 61

` ≈ 302 × qrs/⇡

No Baryon [R=0]

Silk damping

slide-62
SLIDE 62

` ≈ 302 × qrs/⇡

Effect of baryons

slide-63
SLIDE 63

` ≈ 302 × qrs/⇡

Zero-point shift of the

  • scillations

Effect of baryons

slide-64
SLIDE 64

` ≈ 302 × qrs/⇡

WKB factor (1+R)-1/4 and Silk damping compensate the zero- point shift

Effect of baryons

slide-65
SLIDE 65

Effect of Total Matter

Weinberg “Cosmology”, Eq. (6.5.7) where T(q), S(q), θ(q) are “transfer functions” that smoothly interpolate two limits as

q q q

q << qEQ: q >> qEQ: with qEQ = aEQHEQ ~ 0.01 (ΩMh2/0.14) Mpc–1

“EQ” for “matter-radiation Equality epoch”

q q

slide-66
SLIDE 66

` ≈ 302 × qrs/⇡

[ΩMh2=0.07]

Smaller matter density

  • > More potential decay
  • > Larger boost
slide-67
SLIDE 67

Recap

  • The basic structure of the temperature power spectrum is
  • The Sachs-Wolfe “plateau” at low multipoles
  • Sound waves at intermediate multipoles
  • 1st-order tight-coupling
  • Silk damping and Landau damping at high multipoles
  • 2nd-order tight-coupling
slide-68
SLIDE 68

In more details…

  • Decay of gravitational potentials boosts the temperature

power at high multipoles by a factor of 5 compared to

the Sachs-Wolfe plateau

  • Where this boost starts depends on the total matter density
  • Baryon density shifts the zero-point of the oscillation, boosting

the odd peaks relative to the even peaks

  • However, the WKB factor (1+R)–1/4 and damping make the

boosting of the 3rd and 5th peaks not so

  • bvious
slide-69
SLIDE 69

Not quite there yet…

  • The first peak is too low
  • We need to include the “integrated Sachs-Wolfe effect”
  • How to fill zeros between the

peaks?

  • We need to include the Doppler shift of light