Linear Modeling of Baryon-Dark Matter Interactions with CLASS - - PowerPoint PPT Presentation

linear modeling of baryon dark matter interactions with
SMART_READER_LITE
LIVE PREVIEW

Linear Modeling of Baryon-Dark Matter Interactions with CLASS - - PowerPoint PPT Presentation

Interacting Dark Matter Markus R. Mosbech Linear Modeling of Baryon-Dark Matter Interactions with CLASS Introducing a Novel Resonance Interaction Type Markus Rasmussen Mosbech February 21, 2020 Interacting Dark Matter Markus R. Mosbech


slide-1
SLIDE 1

Interacting Dark Matter Markus R. Mosbech

Linear Modeling of Baryon-Dark Matter Interactions with CLASS

Introducing a Novel Resonance Interaction Type Markus Rasmussen Mosbech February 21, 2020

slide-2
SLIDE 2

Interacting Dark Matter Markus R. Mosbech

Overview

1 Evolution 2 Dark Matter Interactions 3 Numerics 4 Results

slide-3
SLIDE 3

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-4
SLIDE 4

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-5
SLIDE 5

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-6
SLIDE 6

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-7
SLIDE 7

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-8
SLIDE 8

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-9
SLIDE 9

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-10
SLIDE 10

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-11
SLIDE 11

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-12
SLIDE 12

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Structure Formation

Primordial perturbations Gravity Pressure

Figure: simulations were performed

at the National Center for Supercomputer Applications by Andrey Kravtsov (The University of Chicago) and Anatoly Klypin (New Mexico State University). Visualizations by Andrey Kravtsov.

slide-13
SLIDE 13

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

The CMB spectrum

Observable Recombination Linear

Figure: by ESA/Planck collaboration

slide-14
SLIDE 14

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

The CMB spectrum

Observable Recombination Linear

Figure: by ESA/Planck collaboration

slide-15
SLIDE 15

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

The CMB spectrum

Observable Recombination Linear

500 1000 1500 2000 2500 l 1 2 3 4 l(l + 1)ClTT ×10−9

slide-16
SLIDE 16

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Perturbation Theory

Small deviation from average Early universe Advantages in Fourier space

10−6 10−5 10−4 10−3 10−2 a 10−2 100 102 104 |δb| (k = 0.5 Mpc−1)

slide-17
SLIDE 17

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Perturbation Theory

Small deviation from average Early universe Advantages in Fourier space

10−6 10−5 10−4 10−3 10−2 a 10−2 100 102 104 |δb| (k = 0.5 Mpc−1)

slide-18
SLIDE 18

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Perturbation Theory

Small deviation from average Early universe Advantages in Fourier space

10−6 10−5 10−4 10−3 10−2 a 10−2 100 102 104 |δb| (k = 0.5 Mpc−1)

slide-19
SLIDE 19

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

The Boltzmann Equations

Dark Matter Baryons Photons ˙ δcdm = −θcdm + 3 ˙ Φ, ˙ θcdm = − ˙ a aθcdm + k2Ψ.

slide-20
SLIDE 20

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

The Boltzmann Equations

Dark Matter Baryons Photons ˙ δb = −θb − 3 ˙ Φ, ˙ θb = − ˙ a aθb + c2

sk2δb

+ 4¯ ργ 3¯ ρb aneσT (θγ − θb) + k2Ψ.

slide-21
SLIDE 21

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

The Boltzmann Equations

Dark Matter Baryons Photons ˙ δγ = −4 3θγ + 4 ˙ Φ ˙ θγ = k2 1 4δγ − σγ

  • + aneσT (θb − θγ) + k2Ψ
slide-22
SLIDE 22

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Temperature

Photons Baryons Cold dark matter EDGES

10−4 10−3 10−2 10−1 100 a 100 101 102 103 104 T [K] Tb Tγ

slide-23
SLIDE 23

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Temperature

Photons Baryons Cold dark matter EDGES

10−4 10−3 10−2 10−1 100 a 100 101 102 103 104 T [K] Tb Tγ

slide-24
SLIDE 24

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Temperature

Photons Baryons Cold dark matter EDGES

10−4 10−3 10−2 10−1 100 a 100 101 102 103 104 T [K] Tb Tγ

slide-25
SLIDE 25

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Temperature

Photons Baryons Cold dark matter EDGES

10−4 10−3 10−2 10−1 100 a 100 101 102 103 104 T [K] Tb Tγ

slide-26
SLIDE 26

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Basic Premise

Nonzero interaction Velocity dependence Nonzero temperature

∆ p χ p χ

slide-27
SLIDE 27

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Basic Premise

Nonzero interaction Velocity dependence Nonzero temperature

∆ p χ p χ

slide-28
SLIDE 28

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Basic Premise

Nonzero interaction Velocity dependence Nonzero temperature

∆ p χ p χ

slide-29
SLIDE 29

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Modifjed Boltzmann equations

Baryons Dark Matter ˙ δb = −θb − 3 ˙ Φ, ˙ θb = − ˙ a aθb + c2

sk2δb

+ 4¯ ργ 3¯ ρb aneσT (θγ − θb) + k2Ψ. + ρχ ρb Rχ (θχ − θb) Rχ = aρb mχ + mH σv FHe

slide-30
SLIDE 30

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Modifjed Boltzmann equations

Baryons Dark Matter ˙ δcdm = −θcdm + 3 ˙ Φ, ˙ θcdm = − ˙ a aθcdm + k2Ψ + Rχ (θb − θχ) Rχ = aρb mχ + mH σv FHe

slide-31
SLIDE 31

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Power Law Interaction

Coulomb Dipole

10−4 10−3 10−2 10−1 100 a 10−1 101 103 T [K] Tb Tb ref. Tχ Tγ

Rχ = acnρbσ0 mχ + mH Tb mH + Tχ mχ + V 2

rms

3 n+1

2

FHe

slide-32
SLIDE 32

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Power Law Interaction

Coulomb Dipole

10−4 10−3 10−2 10−1 100 a 10−4 10−2 100 102 104 T [K] Tb Tb ref. Tχ Tγ

Rχ = acnρbσ0 mχ + mH Tb mH + Tχ mχ + V 2

rms

3 n+1

2

FHe

slide-33
SLIDE 33

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Resonance interaction

Assumptions Velocity dependence

10−4 10−3 10−2 10−1 a 10−7 10−5 10−3 10−1 101 R/aH mχ = 10000 mχ = 10 mχ = 1 mχ = 1000 mχ = 100

Rχ = aρbσ0 mχ + mH xbxχQ3/2 mHmχ K1 Q 2 √xbxχ

  • FHe
slide-34
SLIDE 34

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Resonance interaction

Assumptions Velocity dependence

10−4 10−3 10−2 10−1 a 10−7 10−5 10−3 10−1 101 R/aH mχ = 10000 mχ = 10 mχ = 1 mχ = 1000 mχ = 100

Rχ = aρbσ0 mχ + mH xbxχQ3/2 mHmχ K1 Q 2 √xbxχ

  • FHe
slide-35
SLIDE 35

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

CLASS

Boltzmann solver Linear Modular

Figure: from class-code.net

slide-36
SLIDE 36

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

CLASS

Boltzmann solver Linear Modular

Figure: from class-code.net

slide-37
SLIDE 37

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

CLASS

Boltzmann solver Linear Modular

Figure: from class-code.net

slide-38
SLIDE 38

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Thermodynamics

Interactions Temperature Recombination + reionisation

−20 ∆xe/xe,ref [%] Interacting ΛCDM 101 102 103 104 z 10−3 10−1 xe

slide-39
SLIDE 39

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Thermodynamics

Interactions Temperature Recombination + reionisation

−20 ∆xe/xe,ref [%] Interacting ΛCDM 101 102 103 104 z 10−3 10−1 xe

slide-40
SLIDE 40

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Thermodynamics

Interactions Temperature Recombination + reionisation

−20 ∆xe/xe,ref [%] Interacting ΛCDM 101 102 103 104 z 10−3 10−1 xe

slide-41
SLIDE 41

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Perturbations

Modifjed Boltzmann equations Tight coupling ap- proximation

10−6 10−5 10−4 10−3 10−2 a 10−2 100 102 104 |δb| (k = 0.5 Mpc−1) Interacting ΛCDM

slide-42
SLIDE 42

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Perturbations

Modifjed Boltzmann equations Tight coupling ap- proximation

10−6 10−5 10−4 10−3 10−2 a 10−2 100 102 104 |δb| (k = 0.5 Mpc−1) Interacting ΛCDM

slide-43
SLIDE 43

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Upper limit ∆Cls

500 1,000 1,500 2,000 2,500 1 2 3 4 ·10−9 l l(l + 1)ClTT ΛCDM Power law Resonance

slide-44
SLIDE 44

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Upper limit ∆Cls

−3 −2 −1 ∆TT/TTref [%] Resonance 10 MeV Resonance 100 MeV Power law 1 MeV Power law 100 MeV 500 1000 1500 2000 2500 l −3 −2 −1 1 ∆EE/EEref [%]

slide-45
SLIDE 45

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Coulomb-like

10−4 10−3 10−2 10−1 100 a 10−3 10−1 101 103 T [K] Tb Power law 1 MeV Tχ Power law 1 MeV Tb Power law 100 MeV Tχ Power law 100 MeV

slide-46
SLIDE 46

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Coulomb-like

mχ 1 MeV 100 MeV σ0 (95% CL.) 2.2 × 10−41 cm2 2.3 × 10−41 cm2 Tb(z = 17.2) 0.01 K 0.14 K

slide-47
SLIDE 47

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Resonance

10−4 10−3 10−2 10−1 100 a 10−5 10−3 10−1 101 103 T [K] Tχ Resonance 10 MeV Tb Resonance 10 MeV Tχ Resonance 100 MeV Tb Resonance 100 MeV

slide-48
SLIDE 48

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Resonance

mχ 10 MeV 100 MeV σ0 (95% CL.) 5.2 × 10−9 1.9 × 10−7 ˜ σ0 (95% CL.) 1.2 × 10−51 cm2 1.4 × 10−52 cm2 Tb(z = 17.2) 6.97 K 6.71 K

slide-49
SLIDE 49

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Further Results

ǫ [MeV] mχ = 10 MeV mχ = 100 MeV 1 × 10−11 σ0 = 1.05 × 10−8 σ0 = 3.46 × 10−7 1 × 10−10 σ0 = 2.52 × 10−9 σ0 = 8.36 × 10−8 1 × 10−9 σ0 = 1.39 × 10−9 σ0 = 2.09 × 10−8 1 × 10−8 σ0 = 2.84 × 10−9 σ0 = 1.61 × 10−8 1 × 10−7 σ0 = 9.88 × 10−8 σ0 = 5.29 × 10−8 ǫ [MeV] mχ = 1000 MeV mχ = 5000 MeV 1 × 10−11 σ0 = 1.05 × 10−5 σ0 = 1.25 × 10−4 1 × 10−10 σ0 = 2.59 × 10−6 σ0 = 2.95 × 10−5 1 × 10−9 σ0 = 6.79 × 10−7 σ0 = 7.09 × 10−6 1 × 10−8 σ0 = 1.93 × 10−7 σ0 = 2.25 × 10−6 1 × 10−7 σ0 = 2.07 × 10−7 σ0 = 2.52 × 10−6

Table:

Inferred 95% confjdence upper limits of σ0. Obtained with MontePython using Planck 2015 TT+TE+EE data.

slide-50
SLIDE 50

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Further Results

ǫ [MeV] mχ = 10 MeV mχ = 100 MeV 1 × 10−11 ˜ σ0 = 4.43 × 10−52 cm2 ˜ σ0 = 5.29 × 10−53 cm2 1 × 10−10 ˜ σ0 = 3.36 × 10−52 cm2 ˜ σ0 = 4.04 × 10−52 cm2 1 × 10−9 ˜ σ0 = 5.86 × 10−50 cm2 ˜ σ0 = 3.19 × 10−51 cm2 1 × 10−8 ˜ σ0 = 3.79 × 10−48 cm2 ˜ σ0 = 7.78 × 10−50 cm2 1 × 10−7 ˜ σ0 = 4.17 × 10−45 cm2 ˜ σ0 = 8.08 × 10−48 cm2 ǫ [MeV] mχ = 1000 MeV mχ = 5000 MeV 1 × 10−11 ˜ σ0 = 1.29 × 10−53 cm2 ˜ σ0 = 1.48 × 10−53 cm2 1 × 10−10 ˜ σ0 = 1.00 × 10−52 cm2 ˜ σ0 = 1.10 × 10−52 cm2 1 × 10−9 ˜ σ0 = 8.36 × 10−52 cm2 ˜ σ0 = 8.38 × 10−51 cm2 1 × 10−8 ˜ σ0 = 7.52 × 10−51 cm2 ˜ σ0 = 8.41 × 10−51 cm2 1 × 10−7 ˜ σ0 = 2.55 × 10−49 cm2 ˜ σ0 = 2.98 × 10−49 cm2

Table:

Inferred 95% confjdence upper limits of ˜ σ0. Obtained with MontePython using Planck 2015 TT+TE+EE data.

slide-51
SLIDE 51

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Further Results

100 101 102 103 104

z

10−7 10−5 10−3 10−1 101

Rχ/aH ǫ = 10−10 MeV mχ = 10 MeV ǫ = 10−10 MeV mχ = 100 MeV ǫ = 10−10 MeV mχ = 1000 MeV ǫ = 10−10 MeV mχ = 5000 MeV

Figure: The evolution of the momentum exchange rate Rχ for

ǫ = 1 × 10−10 MeV, using the inferred CMB upper limits for σ0.

slide-52
SLIDE 52

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Further Results

−2 −1

∆CT T

l

/CT T

l,ref [%]

mχ =10 MeV mχ =100 MeV mχ =1000 MeV mχ =5000 MeV

500 1000 1500 2000 2500

l

−3 −2 −1 1

∆CEE

l

/CEE

l,ref [%]

Figure: Percent residuals of the TT (upper) and EE (lower) spectra with

respect to ΛCDM for ǫ = 1 × 10−10 MeV, using the inferred CMB upper limits for σ0.

slide-53
SLIDE 53

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Further Results

100 101 102 103 104

z

10−7 10−5 10−3 10−1 101

Rχ/aH ǫ = 10−11 MeV mχ = 100 MeV ǫ = 10−10 MeV mχ = 100 MeV ǫ = 10−9 MeV mχ = 100 MeV ǫ = 10−8 MeV mχ = 100 MeV ǫ = 10−7 MeV mχ = 100 MeV

Figure: The evolution of the momentum exchange rate Rχ for

mχ = 100 MeV, using the inferred CMB upper limits for σ0.

slide-54
SLIDE 54

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Further Results

−3 −2 −1

∆CT T

l

/CT T

l,ref [%]

ǫ = 10−11 MeV ǫ = 10−10 MeV ǫ = 10−9 MeV ǫ = 10−8 MeV ǫ = 10−7 MeV

500 1000 1500 2000 2500

l

−2

∆CEE

l

/CEE

l,ref [%]

Figure: Percent residuals of the TT (upper) and EE (lower) spectra with

respect to ΛCDM for mχ = 100 MeV, using the inferred CMB upper limits for σ0.

slide-55
SLIDE 55

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Further Results

10−8 10−6 10−4

σ0

10−10 10−8

ǫ [MeV]

10−51 10−48 10−45

˜ σ0 [cm2] mχ = 10 MeV mχ = 100 MeV mχ = 1000 MeV mχ = 5000 MeV

Figure: Inferred 95% upper limits for σ0 (upper fjgure) and ˜

σ0 (lower fjgure) as a function of ǫ for each tested mχ.

slide-56
SLIDE 56

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Further Results

10−8 10−6 10−4

σ0 ǫ = 10−11 MeV ǫ = 10−10 MeV ǫ = 10−9 MeV ǫ = 10−8 MeV ǫ = 10−7 MeV

101 102 103

Mχ [MeV]

10−51 10−48 10−45

˜ σ0 [cm2]

Figure: Inferred 95% upper limits for σ0 (upper fjgure) and ˜

σ0 (lower fjgure) as a function of mχ for each tested ǫ.

slide-57
SLIDE 57

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Conclusion

Upper Limits EDGES Other constraints

slide-58
SLIDE 58

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Conclusion

Upper Limits EDGES Other constraints

slide-59
SLIDE 59

Interacting Dark Matter Markus R. Mosbech Evolution DM Interactions Numerics Results

Conclusion

Upper Limits EDGES Other constraints