Lecture 2 - Power spectrum of gravitational anisotropy - Temperature - - PowerPoint PPT Presentation

lecture 2
SMART_READER_LITE
LIVE PREVIEW

Lecture 2 - Power spectrum of gravitational anisotropy - Temperature - - PowerPoint PPT Presentation

Lecture 2 - Power spectrum of gravitational anisotropy - Temperature anisotropy from sound waves Bennett et al. (1996) COBE 4-year Power Spectrum The SW formula allows us to determine the 3d power spectrum of at the last scattering surface


slide-1
SLIDE 1

Lecture 2

  • Power spectrum of gravitational anisotropy
  • Temperature anisotropy from sound waves
slide-2
SLIDE 2

COBE 4-year Power Spectrum

Bennett et al. (1996)

The SW formula allows us to determine the 3d power

spectrum of φ at

the last scattering surface from Cl.

But how?

slide-3
SLIDE 3
  • But this is not exactly what we want. We want the

statistical average of this quantity.

gives…

Sachs & Wolfe (1967)

∆T(ˆ n) T0 = 1 3Φ(tL, ˆ rL)

SW Power Spectrum

slide-4
SLIDE 4

Power Spectrum of φ

  • Statistical average of the right hand side contains

two-point correlation function

If does not depend on locations (x) but only on separations between two points (r), then

where we defined consequence of “statistical homogeneity” φ and used

slide-5
SLIDE 5

Power Spectrum of φ

  • In addition, if depends only
  • n the magnitude of the separation r and not on the

directions, then

Power spectrum!

Generic definition of the power spectrum for statistically homogeneous and isotropic fluctuations

slide-6
SLIDE 6

SW Power Spectrum

  • Thus, the power spectrum of the CMB in the SW limit is
  • In the flat-sky approximation,

Perpendicular wavenumber, (qperp)2

slide-7
SLIDE 7

SW Power Spectrum

  • Thus, the power spectrum of the CMB in the SW limit is
  • In the flat-sky approximation,

For a power-law form, , we get

slide-8
SLIDE 8

SW Power Spectrum

  • Thus, the power spectrum of the CMB in the SW limit is
  • In the flat-sky approximation,

For a power-law form, , we get

n=1

full-sky correction

slide-9
SLIDE 9

n=1 n=1.2 ± 0.3

(68%CL)

Bennett et al. (1996)

slide-10
SLIDE 10

COBE 4-year Power Spectrum

Bennett et al. (1996)

slide-11
SLIDE 11

WMAP 9-year Power Spectrum

Bennett et al. (2013)

slide-12
SLIDE 12

Planck 29-mo Power Spectrum

Planck Collaboration (2016)

slide-13
SLIDE 13

Planck 29-mo Power Spectrum

Planck Collaboration (2016)

Clearly, the SW prediction does not fit! Missing physics: Hydrodynamics (sound waves)

slide-14
SLIDE 14
slide-15
SLIDE 15

La Soupe Miso Cosmique

  • When matter and radiation were hotter than 3000 K,

matter was completely ionised. The Universe was filled with plasma, which behaves just like a soup

  • Think about a Miso soup (if you know what it is).

Imagine throwing Tofus into a Miso soup, while changing the density of Miso

  • And imagine watching how ripples are created and

propagate throughout the soup

slide-16
SLIDE 16
slide-17
SLIDE 17

This is a viscous fluid, in which the amplitude of sound waves damps at shorter wavelength

slide-18
SLIDE 18
slide-19
SLIDE 19

When do sound waves become important?

  • In other words, when would the Sachs-Wolfe approximation

(purely gravitational effects) become invalid?

  • The key to the answer: Sound-crossing Time
  • Sound waves cannot alter temperature anisotropy at a

given angular scale if there was not enough time for sound waves to propagate to the corresponding distance at the last-scattering surface

  • The distance traveled by sound waves within a given

time = The Sound Horizon

slide-20
SLIDE 20

Comoving Photon Horizon

  • First, the comoving distance traveled by photons is given

by setting the space-time distance to be null:

ds2 = −c2dt2 + a2(t)dr2 = 0

rphoton = c Z t dt0 a(t0)

slide-21
SLIDE 21

Comoving Sound Horizon

  • Then, we replace the speed of light with a time-

dependent speed of sound:

rs = Z t dt0 a(t0)cs(t0)

  • We cannot ignore the effects of sound waves if qrs > 1
slide-22
SLIDE 22

Sound Speed

  • Sound speed of an adiabatic fluid is given by
  • δP: pressure perturbation
  • δρ: density perturbation
  • For a baryon-photon system:

We can ignore the baryon pressure because it is much smaller than the photon pressure

slide-23
SLIDE 23

Sound Speed

  • Using the adiabatic relationship between photons and baryons:
  • and pressure-density relation of a relativistic fluid, δPγ=δργ/3,

We obtain

[i.e., the ratio of the number densities of baryons and photons is equal everywhere]

  • Or equivalently

where sound speed is reduced!

slide-24
SLIDE 24

Value of R?

  • The baryon mass density goes like a–3, whereas the

photon energy density goes like a–4. Thus, the ratio of the two, R, goes like a.

  • The proportionality constant is:

where we used for

slide-25
SLIDE 25

Value of R?

  • The baryon mass density goes like a–3, whereas the

photon energy density goes like a–4. Thus, the ratio of the two, R, goes like a.

  • The proportionality constant is:

where we used for

For the last-scattering redshift of zL=1090 (or last-scattering temperature of TL=2974 K),

rs = 145.3 Mpc

We cannot ignore the effects of sound waves if qrs>1. Since l~qrL, this means

l > rL/rs = 96

where we used rL=13.95 Gpc

slide-26
SLIDE 26

Creation of Sound Waves: Basic Equations

  • 1. Conservation equations (energy and momentum)
  • 2. Equation of state, relating pressure to energy density
  • 3. General relativistic version of the “Poisson equation”,

relating gravitational potential to energy density

  • 4. Evolution of the “anisotropic stress” (viscosity)

P = P(ρ)

slide-27
SLIDE 27
  • Total energy conservation:
  • C.f., Total energy conservation [unperturbed]

Energy Conservation

( )

velocity potential anisotropic stress: [or, viscosity]

vα = 1 arδuα

slide-28
SLIDE 28

Energy Conservation

  • Total energy conservation:
  • Again, this is the effect of locally-defined inhomogeneous

scale factor, i.e.,

  • The spatial metric is given by
  • Thus, locally we can define a new scale factor:

ds2 = a2(t) exp(−2Ψ)dx2 ˜ a(t, x) = a(t) exp(−Ψ)

slide-29
SLIDE 29

Energy Conservation

  • Total energy conservation:
  • Momentum flux going outward (inward) -> reduction

(increase) in the energy density

C.f., for a non-expanding medium:

˙ ρ + r · (ρv) = 0

( )

slide-30
SLIDE 30

Momentum Conservation

  • Total momentum conservation
  • Cosmological redshift of the momentum
  • Gravitational force given by potential gradient
  • Force given by pressure gradient
  • Force given by gradient of anisotropic stress

vα = 1 arδuα

slide-31
SLIDE 31
  • Pressure of non-relativistic species (i.e., baryons and cold

dark matter) can be ignored relative to the energy density. Thus, we set them to zero: PB=0=PD and δPB=0=δPD

  • Unperturbed pressure of relativistic species (i.e., photons

and relativistic neutrinos) is given by the third of the energy density, i.e., Pγ=ργ/3 and Pν=ρν/3

  • Perturbed pressure involves contributions from the bulk

viscosity:

Equation of State

δPγ =

δPν =

slide-32
SLIDE 32
  • Pressure of non-relativistic species (i.e., baryons and cold

dark matter) can be ignored relative to the energy density. Thus, we set them to zero: PB=0=PD and δPB=0=δPD

  • Unperturbed pressure of relativistic species (i.e., photons

and relativistic neutrinos) is given by the third of the energy density, i.e., Pγ=ργ/3 and Pν=ρν/3

  • Perturbed pressure involves contributions from the bulk

viscosity:

Equation of State

δPγ =

δPν =

The reason for this is that

trace of the stress-energy

  • f relativistic species

vanishes: ∑μ=0,1,2,3 Τμμ = 0

T 0

0 + 3

X

i=1

T i

i = ρ + 3P + r2π = 0

slide-33
SLIDE 33

Two Remarks

  • In the standard scenario:
  • Energy densities are conserved separately; thus we do

not need to sum over all species

  • Momentum densities of photons and baryons are NOT

conserved separately but they are coupled via Thomson scattering. This must be taken into account when writing down separate conservation equations

slide-34
SLIDE 34
  • Fourier transformation replaces

Conservation Equations for Photons and Baryons

r2 ! q2

momentum transfer via scattering

slide-35
SLIDE 35
  • Fourier transformation replaces

Conservation Equations for Photons and Baryons

r2 ! q2

what about photon’s viscosity?

slide-36
SLIDE 36

Formation of a Photon-baryon Fluid

  • Photons are not a fluid. Photons free-stream at

the speed of light

  • The conservation equations are not enough because we

need to specify the evolution of viscosity

  • Solving for viscosity requires information of the phase-space

distribution function of photons: Boltzmann equation

  • However, frequent scattering of photons with baryons* can

make photons behave as a fluid: Photon-baryon fluid

Peebles & Yu (1970); Sunyaev & Zeldovich (1970) *Photons scatter with electrons via Thomson scattering. Protons scatter with electrons via Coulomb scattering. Thus we can say, effectively, photons scatter with baryons

slide-37
SLIDE 37
  • Fourier transformation replaces

Let’s solve them!

r2 ! q2

slide-38
SLIDE 38

Tight-coupling Approximation

  • 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-39
SLIDE 39

Tight-coupling Approximation

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

the scale factor, we obtain

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

we obtain

Wave Equation, with the speed of sound of cs2 = 1/3(1+R)!

slide-40
SLIDE 40

Sound Wave!

  • To simplify the equation, let’s first look at the high-

frequency solution

  • Specifically, we take q >> aH (the wavelength of

fluctuations is much shorter than the Hubble length). Then we can ignore time derivatives of R and Ψ because they evolve in the Hubble time scale:

Peebles & Yu (1970); Sunyaev & Zeldovich (1970) Solution: SOUND WAVE!

slide-41
SLIDE 41

Recap

  • Photons are not a fluid; but Thomson scattering couples

photons to baryons, forming a photon-baryon fluid

  • The reduced sound speed, cs2=1/3(1+R), emerges

automatically

  • δργ/4ργ is the temperature anisotropy at the bottom of the

potential well. Adding gravitational redshift, the observed temperature anisotropy is δργ/4ργ + Φ,

which is given by

slide-42
SLIDE 42
slide-43
SLIDE 43
slide-44
SLIDE 44

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-45
SLIDE 45

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-46
SLIDE 46

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-47
SLIDE 47

ζ: 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-48
SLIDE 48

ζ: 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-49
SLIDE 49

ζ: 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-50
SLIDE 50

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-51
SLIDE 51

Entered the horizon during the radiation era

slide-52
SLIDE 52

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

slide-53
SLIDE 53

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-54
SLIDE 54

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-55
SLIDE 55
slide-56
SLIDE 56

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-57
SLIDE 57

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-58
SLIDE 58

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-59
SLIDE 59

Better Solution in Radiation-dominated Era

Then the equation simplifies to

where

The solution is

We rewrite this using the formula for trigonometry:

sin(ϕ − ϕ0) = sin(ϕ) cos(ϕ0) − cos(ϕ) sin(ϕ0)

<latexit sha1_base64="DurnUGb/bp0JcYqiyolGaYVoUhg=">ACN3icbVDLSgMxFM34rPU16tLNYJF2Fi0zIuhGKLpxJRXsAzpDyaRpG5pJhiRTKEP/yo2/4U43LhRx6x+YtoOMrRdCTs45l5t7gogSqRznxVhZXVvf2Mxt5bd3dvf2zYPDhuSxQLiOWiFUCJKWG4roiuBUJDMOA4mYwvJnqzREWknD2oMYR9kPYZ6RHEFSa6ph3niSs5I2giAaknN5F+ypL2x7isvSrlbNPO2s2h2z4FScWVnLwE1BAaRV65jPXpejOMRMIQqlbLtOpPwECkUQxZO8F0scQTSEfdzWkMEQSz+Z7T2xTjXTtXpc6MOUNWOzHQkMpRyHgXaGUA3kojYl/9Pasepd+glhUawQ/NBvZhailvTEK0uERgpOtYAIkH0Xy0gAIipaPO6xDcxZWXQeOs4joV9/68UL1O48iBY3ACSsAF6AKbkEN1AECj+AVvIMP48l4Mz6Nr7l1xUh7jsCfMr5/APUrq/A=</latexit><latexit sha1_base64="DurnUGb/bp0JcYqiyolGaYVoUhg=">ACN3icbVDLSgMxFM34rPU16tLNYJF2Fi0zIuhGKLpxJRXsAzpDyaRpG5pJhiRTKEP/yo2/4U43LhRx6x+YtoOMrRdCTs45l5t7gogSqRznxVhZXVvf2Mxt5bd3dvf2zYPDhuSxQLiOWiFUCJKWG4roiuBUJDMOA4mYwvJnqzREWknD2oMYR9kPYZ6RHEFSa6ph3niSs5I2giAaknN5F+ypL2x7isvSrlbNPO2s2h2z4FScWVnLwE1BAaRV65jPXpejOMRMIQqlbLtOpPwECkUQxZO8F0scQTSEfdzWkMEQSz+Z7T2xTjXTtXpc6MOUNWOzHQkMpRyHgXaGUA3kojYl/9Pasepd+glhUawQ/NBvZhailvTEK0uERgpOtYAIkH0Xy0gAIipaPO6xDcxZWXQeOs4joV9/68UL1O48iBY3ACSsAF6AKbkEN1AECj+AVvIMP48l4Mz6Nr7l1xUh7jsCfMr5/APUrq/A=</latexit><latexit sha1_base64="DurnUGb/bp0JcYqiyolGaYVoUhg=">ACN3icbVDLSgMxFM34rPU16tLNYJF2Fi0zIuhGKLpxJRXsAzpDyaRpG5pJhiRTKEP/yo2/4U43LhRx6x+YtoOMrRdCTs45l5t7gogSqRznxVhZXVvf2Mxt5bd3dvf2zYPDhuSxQLiOWiFUCJKWG4roiuBUJDMOA4mYwvJnqzREWknD2oMYR9kPYZ6RHEFSa6ph3niSs5I2giAaknN5F+ypL2x7isvSrlbNPO2s2h2z4FScWVnLwE1BAaRV65jPXpejOMRMIQqlbLtOpPwECkUQxZO8F0scQTSEfdzWkMEQSz+Z7T2xTjXTtXpc6MOUNWOzHQkMpRyHgXaGUA3kojYl/9Pasepd+glhUawQ/NBvZhailvTEK0uERgpOtYAIkH0Xy0gAIipaPO6xDcxZWXQeOs4joV9/68UL1O48iBY3ACSsAF6AKbkEN1AECj+AVvIMP48l4Mz6Nr7l1xUh7jsCfMr5/APUrq/A=</latexit><latexit sha1_base64="DurnUGb/bp0JcYqiyolGaYVoUhg=">ACN3icbVDLSgMxFM34rPU16tLNYJF2Fi0zIuhGKLpxJRXsAzpDyaRpG5pJhiRTKEP/yo2/4U43LhRx6x+YtoOMrRdCTs45l5t7gogSqRznxVhZXVvf2Mxt5bd3dvf2zYPDhuSxQLiOWiFUCJKWG4roiuBUJDMOA4mYwvJnqzREWknD2oMYR9kPYZ6RHEFSa6ph3niSs5I2giAaknN5F+ypL2x7isvSrlbNPO2s2h2z4FScWVnLwE1BAaRV65jPXpejOMRMIQqlbLtOpPwECkUQxZO8F0scQTSEfdzWkMEQSz+Z7T2xTjXTtXpc6MOUNWOzHQkMpRyHgXaGUA3kojYl/9Pasepd+glhUawQ/NBvZhailvTEK0uERgpOtYAIkH0Xy0gAIipaPO6xDcxZWXQeOs4joV9/68UL1O48iBY3ACSsAF6AKbkEN1AECj+AVvIMP48l4Mz6Nr7l1xUh7jsCfMr5/APUrq/A=</latexit>
slide-60
SLIDE 60

Better Solution in Radiation-dominated Era

Then the equation simplifies to

where

The solution is

where

slide-61
SLIDE 61

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-62
SLIDE 62

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-63
SLIDE 63

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-64
SLIDE 64

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-65
SLIDE 65

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-66
SLIDE 66

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-67
SLIDE 67

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-68
SLIDE 68

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-69
SLIDE 69

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-70
SLIDE 70

Sound Wave Solution in the Radiation-dominated Era

The solution is

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

slide-71
SLIDE 71

Sound Wave Solution in the Radiation-dominated Era

The solution is

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

slide-72
SLIDE 72

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-73
SLIDE 73

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-74
SLIDE 74
slide-75
SLIDE 75

Roles of viscosity

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

Silk (1968) “Silk damping”

slide-76
SLIDE 76

High-frequency solution without neutrino viscosity

The solution is

where

ϕ 1

<latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit>

ϕ 1

<latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit>
slide-77
SLIDE 77

High-frequency solution with neutrino viscosity

The solution is

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

ϕ 1

<latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit>

ϕ 1

<latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit><latexit sha1_base64="bNhabzp0K7ZXJMGgnjOef5tGyL8=">AB83icbVBNS8NAEJ3Ur1q/qh69LBbBU0lE0GPRi8cK9gOaUDbTbp0s1l2N4US+je8eFDEq3/Gm/GTZuDtj4YeLw3w8y8UHKmjet+O5WNza3tnepubW/4PCofnzS1WmCO2QlKeqH2JNORO0Y5jhtC8VxUnIaS+c3Bd+b0qVZql4MjNJgwTHgkWMYGMl359iJcfMj2PkDesNt+kugNaJV5IGlGgP61/+KCVZQoUhHGs98Fxpghwrwin85qfaSoxmeCYDiwVOKE6yBc3z9GFVUYoSpUtYdBC/T2R40TrWRLazgSbsV71CvE/b5CZ6DbImZCZoYIsF0UZRyZFRQBoxBQlhs8swUQxeysiY6wMTamg3BW315nXSvmp7b9B6vG627Mo4qnME5XIHN9CB2hDBwhIeIZXeHMy58V5dz6WrRWnDmFP3A+fwCEuJFT</latexit>
slide-78
SLIDE 78

High-frequency solution with neutrino viscosity

Using the formula for trigonometry, we write

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

slide-79
SLIDE 79

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-80
SLIDE 80

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-81
SLIDE 81

Yesterday: 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-82
SLIDE 82

Today: 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-83
SLIDE 83

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)

= −32 45 ¯ ργ σT ¯ ne ∂i∂jδu

<latexit sha1_base64="x5ZMo5o+0rukqRjQKVd9YinM=">ACQ3icbZDLSiQxFIZTOjNqO5dWl7MJ0whubKrUQTeC6MalA7bKdJriVPpUd8YkVSQpoSnq3dz4Au58ATcuFHErTKq7GbzMgcCf/z8nly/JpbAuDG+CmdkPHz/NzS80Fj9/+fqtubR8YrPCcOzwTGbmLAGLUmjsOEknuUGQSUST5Pzgzo/vUBjRaP3SjHnoKBFqng4LwVN3/vrPUAC83N6py62c12bAEDPDLGYDUAqklkxUBCXjIOkx1WdUx1jxXIwToCMxT/1h/VROqBF3GyF7XBc9L2IpqJFpnUN69ZP+OFQu24BGu7UZi7XlmfyVWDVZYzIGfwC7XmpQaHvlmEFV73Tp2lm/NKOjt2XEyUoa0cq8Z0K3NC+zWrzf1m3cOlOrxQ6LxqPrkoLSR1Ga2B0r4wyJ0ceQHcCP9WyofgKTqPveEhRG+/F6cbLSjsB392mrt7U9xzJPv5AdZIxHZJnvkByRDuHktySe/IQXAV3wWPwNGmdCaYzK+RVBc9/AUcjsvA=</latexit><latexit sha1_base64="x5ZMo5o+0rukqRjQKVd9YinM=">ACQ3icbZDLSiQxFIZTOjNqO5dWl7MJ0whubKrUQTeC6MalA7bKdJriVPpUd8YkVSQpoSnq3dz4Au58ATcuFHErTKq7GbzMgcCf/z8nly/JpbAuDG+CmdkPHz/NzS80Fj9/+fqtubR8YrPCcOzwTGbmLAGLUmjsOEknuUGQSUST5Pzgzo/vUBjRaP3SjHnoKBFqng4LwVN3/vrPUAC83N6py62c12bAEDPDLGYDUAqklkxUBCXjIOkx1WdUx1jxXIwToCMxT/1h/VROqBF3GyF7XBc9L2IpqJFpnUN69ZP+OFQu24BGu7UZi7XlmfyVWDVZYzIGfwC7XmpQaHvlmEFV73Tp2lm/NKOjt2XEyUoa0cq8Z0K3NC+zWrzf1m3cOlOrxQ6LxqPrkoLSR1Ga2B0r4wyJ0ceQHcCP9WyofgKTqPveEhRG+/F6cbLSjsB392mrt7U9xzJPv5AdZIxHZJnvkByRDuHktySe/IQXAV3wWPwNGmdCaYzK+RVBc9/AUcjsvA=</latexit><latexit sha1_base64="x5ZMo5o+0rukqRjQKVd9YinM=">ACQ3icbZDLSiQxFIZTOjNqO5dWl7MJ0whubKrUQTeC6MalA7bKdJriVPpUd8YkVSQpoSnq3dz4Au58ATcuFHErTKq7GbzMgcCf/z8nly/JpbAuDG+CmdkPHz/NzS80Fj9/+fqtubR8YrPCcOzwTGbmLAGLUmjsOEknuUGQSUST5Pzgzo/vUBjRaP3SjHnoKBFqng4LwVN3/vrPUAC83N6py62c12bAEDPDLGYDUAqklkxUBCXjIOkx1WdUx1jxXIwToCMxT/1h/VROqBF3GyF7XBc9L2IpqJFpnUN69ZP+OFQu24BGu7UZi7XlmfyVWDVZYzIGfwC7XmpQaHvlmEFV73Tp2lm/NKOjt2XEyUoa0cq8Z0K3NC+zWrzf1m3cOlOrxQ6LxqPrkoLSR1Ga2B0r4wyJ0ceQHcCP9WyofgKTqPveEhRG+/F6cbLSjsB392mrt7U9xzJPv5AdZIxHZJnvkByRDuHktySe/IQXAV3wWPwNGmdCaYzK+RVBc9/AUcjsvA=</latexit><latexit sha1_base64="x5ZMo5o+0rukqRjQKVd9YinM=">ACQ3icbZDLSiQxFIZTOjNqO5dWl7MJ0whubKrUQTeC6MalA7bKdJriVPpUd8YkVSQpoSnq3dz4Au58ATcuFHErTKq7GbzMgcCf/z8nly/JpbAuDG+CmdkPHz/NzS80Fj9/+fqtubR8YrPCcOzwTGbmLAGLUmjsOEknuUGQSUST5Pzgzo/vUBjRaP3SjHnoKBFqng4LwVN3/vrPUAC83N6py62c12bAEDPDLGYDUAqklkxUBCXjIOkx1WdUx1jxXIwToCMxT/1h/VROqBF3GyF7XBc9L2IpqJFpnUN69ZP+OFQu24BGu7UZi7XlmfyVWDVZYzIGfwC7XmpQaHvlmEFV73Tp2lm/NKOjt2XEyUoa0cq8Z0K3NC+zWrzf1m3cOlOrxQ6LxqPrkoLSR1Ga2B0r4wyJ0ceQHcCP9WyofgKTqPveEhRG+/F6cbLSjsB392mrt7U9xzJPv5AdZIxHZJnvkByRDuHktySe/IQXAV3wWPwNGmdCaYzK+RVBc9/AUcjsvA=</latexit>

γ γ

slide-84
SLIDE 84

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)

= −32 45 ¯ ργ σT ¯ ne ∂i∂jδu

<latexit sha1_base64="x5ZMo5o+0rukqRjQKVd9YinM=">ACQ3icbZDLSiQxFIZTOjNqO5dWl7MJ0whubKrUQTeC6MalA7bKdJriVPpUd8YkVSQpoSnq3dz4Au58ATcuFHErTKq7GbzMgcCf/z8nly/JpbAuDG+CmdkPHz/NzS80Fj9/+fqtubR8YrPCcOzwTGbmLAGLUmjsOEknuUGQSUST5Pzgzo/vUBjRaP3SjHnoKBFqng4LwVN3/vrPUAC83N6py62c12bAEDPDLGYDUAqklkxUBCXjIOkx1WdUx1jxXIwToCMxT/1h/VROqBF3GyF7XBc9L2IpqJFpnUN69ZP+OFQu24BGu7UZi7XlmfyVWDVZYzIGfwC7XmpQaHvlmEFV73Tp2lm/NKOjt2XEyUoa0cq8Z0K3NC+zWrzf1m3cOlOrxQ6LxqPrkoLSR1Ga2B0r4wyJ0ceQHcCP9WyofgKTqPveEhRG+/F6cbLSjsB392mrt7U9xzJPv5AdZIxHZJnvkByRDuHktySe/IQXAV3wWPwNGmdCaYzK+RVBc9/AUcjsvA=</latexit><latexit sha1_base64="x5ZMo5o+0rukqRjQKVd9YinM=">ACQ3icbZDLSiQxFIZTOjNqO5dWl7MJ0whubKrUQTeC6MalA7bKdJriVPpUd8YkVSQpoSnq3dz4Au58ATcuFHErTKq7GbzMgcCf/z8nly/JpbAuDG+CmdkPHz/NzS80Fj9/+fqtubR8YrPCcOzwTGbmLAGLUmjsOEknuUGQSUST5Pzgzo/vUBjRaP3SjHnoKBFqng4LwVN3/vrPUAC83N6py62c12bAEDPDLGYDUAqklkxUBCXjIOkx1WdUx1jxXIwToCMxT/1h/VROqBF3GyF7XBc9L2IpqJFpnUN69ZP+OFQu24BGu7UZi7XlmfyVWDVZYzIGfwC7XmpQaHvlmEFV73Tp2lm/NKOjt2XEyUoa0cq8Z0K3NC+zWrzf1m3cOlOrxQ6LxqPrkoLSR1Ga2B0r4wyJ0ceQHcCP9WyofgKTqPveEhRG+/F6cbLSjsB392mrt7U9xzJPv5AdZIxHZJnvkByRDuHktySe/IQXAV3wWPwNGmdCaYzK+RVBc9/AUcjsvA=</latexit><latexit sha1_base64="x5ZMo5o+0rukqRjQKVd9YinM=">ACQ3icbZDLSiQxFIZTOjNqO5dWl7MJ0whubKrUQTeC6MalA7bKdJriVPpUd8YkVSQpoSnq3dz4Au58ATcuFHErTKq7GbzMgcCf/z8nly/JpbAuDG+CmdkPHz/NzS80Fj9/+fqtubR8YrPCcOzwTGbmLAGLUmjsOEknuUGQSUST5Pzgzo/vUBjRaP3SjHnoKBFqng4LwVN3/vrPUAC83N6py62c12bAEDPDLGYDUAqklkxUBCXjIOkx1WdUx1jxXIwToCMxT/1h/VROqBF3GyF7XBc9L2IpqJFpnUN69ZP+OFQu24BGu7UZi7XlmfyVWDVZYzIGfwC7XmpQaHvlmEFV73Tp2lm/NKOjt2XEyUoa0cq8Z0K3NC+zWrzf1m3cOlOrxQ6LxqPrkoLSR1Ga2B0r4wyJ0ceQHcCP9WyofgKTqPveEhRG+/F6cbLSjsB392mrt7U9xzJPv5AdZIxHZJnvkByRDuHktySe/IQXAV3wWPwNGmdCaYzK+RVBc9/AUcjsvA=</latexit><latexit sha1_base64="x5ZMo5o+0rukqRjQKVd9YinM=">ACQ3icbZDLSiQxFIZTOjNqO5dWl7MJ0whubKrUQTeC6MalA7bKdJriVPpUd8YkVSQpoSnq3dz4Au58ATcuFHErTKq7GbzMgcCf/z8nly/JpbAuDG+CmdkPHz/NzS80Fj9/+fqtubR8YrPCcOzwTGbmLAGLUmjsOEknuUGQSUST5Pzgzo/vUBjRaP3SjHnoKBFqng4LwVN3/vrPUAC83N6py62c12bAEDPDLGYDUAqklkxUBCXjIOkx1WdUx1jxXIwToCMxT/1h/VROqBF3GyF7XBc9L2IpqJFpnUN69ZP+OFQu24BGu7UZi7XlmfyVWDVZYzIGfwC7XmpQaHvlmEFV73Tp2lm/NKOjt2XEyUoa0cq8Z0K3NC+zWrzf1m3cOlOrxQ6LxqPrkoLSR1Ga2B0r4wyJ0ceQHcCP9WyofgKTqPveEhRG+/F6cbLSjsB392mrt7U9xzJPv5AdZIxHZJnvkByRDuHktySe/IQXAV3wWPwNGmdCaYzK+RVBc9/AUcjsvA=</latexit>

Given by the spatial gradient of

the velocity field

  • a well-known result in fluid dynamics

γ γ

slide-85
SLIDE 85

Damped Oscillator

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

we obtain, for q >> aH,

New term, giving damping! where

slide-86
SLIDE 86

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-87
SLIDE 87

Damped Oscillator

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

we obtain, for q >> aH,

New term, giving damping! Exponential dampling! SOLUTION:

≈ exp

  • −q2/σT ¯

neH

  • <latexit sha1_base64="RUibAPBZ1AdJB7LbQmlyRbMp6vk=">ACIHicbVBNTxsxEPVCy0daINAjF4uoEj07KJKcET0whEkAkhxWM06s4mF1+vas4holZ/CpX+lx5aVfRWfg1OyIECTxrp6b0ZzczLrFae4vhfNDf/5u3C4tJy4937ldW15vrGmS8rJ7EjS126iw8amWwQ4o0XliHUGQaz7OrxP/BqdV6U5pZHFXgEDo3IlgYKUNvcEWOvKGy7wxgqNOW1/na5uyO8GhSQ1kKC5qdjkYHjJkV+JwaDOlT2mzF7XgK/pIkM9JiMxynzb+iX8qQENSg/fdJLbUq8GRkhrHDVF5tCvYIDdQA0U6Hv19MEx/xiUPs9LF8oQn6pPJ2ovB8VWegsgIb+uTcRX/O6FeX7vVoZWxEa+bgorzSnk/S4n3lUJIeBQLSqXArl0NwIClk2ghJM9fknOdtJ3E5OvrQODmdxLFNtsW2WcL2AE7YseswyS7ZT/YL/Y7+h79jP5Ed4+tc9Fs5gP7D9H9A/Slos4=</latexit><latexit sha1_base64="RUibAPBZ1AdJB7LbQmlyRbMp6vk=">ACIHicbVBNTxsxEPVCy0daINAjF4uoEj07KJKcET0whEkAkhxWM06s4mF1+vas4holZ/CpX+lx5aVfRWfg1OyIECTxrp6b0ZzczLrFae4vhfNDf/5u3C4tJy4937ldW15vrGmS8rJ7EjS126iw8amWwQ4o0XliHUGQaz7OrxP/BqdV6U5pZHFXgEDo3IlgYKUNvcEWOvKGy7wxgqNOW1/na5uyO8GhSQ1kKC5qdjkYHjJkV+JwaDOlT2mzF7XgK/pIkM9JiMxynzb+iX8qQENSg/fdJLbUq8GRkhrHDVF5tCvYIDdQA0U6Hv19MEx/xiUPs9LF8oQn6pPJ2ovB8VWegsgIb+uTcRX/O6FeX7vVoZWxEa+bgorzSnk/S4n3lUJIeBQLSqXArl0NwIClk2ghJM9fknOdtJ3E5OvrQODmdxLFNtsW2WcL2AE7YseswyS7ZT/YL/Y7+h79jP5Ed4+tc9Fs5gP7D9H9A/Slos4=</latexit><latexit sha1_base64="RUibAPBZ1AdJB7LbQmlyRbMp6vk=">ACIHicbVBNTxsxEPVCy0daINAjF4uoEj07KJKcET0whEkAkhxWM06s4mF1+vas4holZ/CpX+lx5aVfRWfg1OyIECTxrp6b0ZzczLrFae4vhfNDf/5u3C4tJy4937ldW15vrGmS8rJ7EjS126iw8amWwQ4o0XliHUGQaz7OrxP/BqdV6U5pZHFXgEDo3IlgYKUNvcEWOvKGy7wxgqNOW1/na5uyO8GhSQ1kKC5qdjkYHjJkV+JwaDOlT2mzF7XgK/pIkM9JiMxynzb+iX8qQENSg/fdJLbUq8GRkhrHDVF5tCvYIDdQA0U6Hv19MEx/xiUPs9LF8oQn6pPJ2ovB8VWegsgIb+uTcRX/O6FeX7vVoZWxEa+bgorzSnk/S4n3lUJIeBQLSqXArl0NwIClk2ghJM9fknOdtJ3E5OvrQODmdxLFNtsW2WcL2AE7YseswyS7ZT/YL/Y7+h79jP5Ed4+tc9Fs5gP7D9H9A/Slos4=</latexit><latexit sha1_base64="RUibAPBZ1AdJB7LbQmlyRbMp6vk=">ACIHicbVBNTxsxEPVCy0daINAjF4uoEj07KJKcET0whEkAkhxWM06s4mF1+vas4holZ/CpX+lx5aVfRWfg1OyIECTxrp6b0ZzczLrFae4vhfNDf/5u3C4tJy4937ldW15vrGmS8rJ7EjS126iw8amWwQ4o0XliHUGQaz7OrxP/BqdV6U5pZHFXgEDo3IlgYKUNvcEWOvKGy7wxgqNOW1/na5uyO8GhSQ1kKC5qdjkYHjJkV+JwaDOlT2mzF7XgK/pIkM9JiMxynzb+iX8qQENSg/fdJLbUq8GRkhrHDVF5tCvYIDdQA0U6Hv19MEx/xiUPs9LF8oQn6pPJ2ovB8VWegsgIb+uTcRX/O6FeX7vVoZWxEa+bgorzSnk/S4n3lUJIeBQLSqXArl0NwIClk2ghJM9fknOdtJ3E5OvrQODmdxLFNtsW2WcL2AE7YseswyS7ZT/YL/Y7+h79jP5Ed4+tc9Fs5gP7D9H9A/Slos4=</latexit>
slide-88
SLIDE 88

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-89
SLIDE 89

Diffusion Length

  • The mean free path of the photon between scatterings is

(σTne)–1

  • Below this scale, you do not have a photon-baryon fluid
  • The number of scatterings per Hubble time is

Nscattering=σTne/H

  • Then, the length traveled by photons by random walks

within the Hubble time is (σTne)–1 times √Nscatterings

  • The diffusion length is thus (σTne)–1 times √Nscatterings =

(σTneH)–1/2

slide-90
SLIDE 90

Diffusion Damping

  • Diffusion mixes hot and cold photons -> Damping of anisotropies

by Wayne Hu

slide-91
SLIDE 91

Planck Collaboration (2016)

Sachs-Wolfe Sound Wave Silk Damping?

slide-92
SLIDE 92

Additional Damping

fuzziness

( )

  • 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-93
SLIDE 93

Planck Collaboration (2016)

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

slide-94
SLIDE 94

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 Fuzziness damping at high multipoles
  • 2nd-order tight-coupling