銀河外ガンマ線による宇宙暗黒物質探索
Kanazawa University Kyoto, August 6, 2018 Koji Ishiwata Based on
- JCAP 1505 (2015) 05, 024 (with S. Ando)
- JCAP 1606 (2016) 06, 045 (with S. Ando)
1. Introduction After the Higgs discovery in 2012, The Standard - - PowerPoint PPT Presentation
Koji Ishiwata Kanazawa University Based on JCAP 1505 (2015) 05, 024 (with S. Ando) JCAP 1606 (2016) 06, 045 (with S. Ando) Kyoto, August 6, 2018 1. Introduction After the Higgs
銀河外ガンマ線による宇宙暗黒物質探索
Kanazawa University Kyoto, August 6, 2018 Koji Ishiwata Based on
After the Higgs discovery in 2012,
theory below a TeV scale
side + something
Cosmological issues:
Cosmological issues:
+ something (DM, …)
XENON1T Fermi-LAT
DM searches
LHC CAST
XENON1T Fermi-LAT
DM searches
LHC CAST
XENON1T Fermi-LAT
DM searches
LHC CAST
NLO calculation@QCD
XENON1T Fermi-LAT
DM searches
LHC CAST
Motivations for indirect DM search (theoretical side)
particles very weakly, i.e., to open the discussion of its stability
AMS-02 ’13
|Rigidity| [GV] 100 200 300 400 500AMS-02 ’16
Excesses over 100 GeV Motivations for indirect DM search (experimental side)
AMS-02 ’13
|Rigidity| [GV] 100 200 300 400 500AMS-02 ’16 Ibe, Matsumoto, Shirai, Yagagida ’14
e.g.,
0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 10 100 1000 Positron Fraction E [GeV]Hamaguchi, Moroi, Nakayama ’15
Decaying DM Decaying/Annihilating DM Motivations for indirect DM search (experimental side)
AMS-02 ’16
indication of a DM signal for DM masses near 80 GeV
Cuoco, Krämer, Korsmeier ’17 Cui, Yuan, Tsai, Fan ’17
4.5σ Motivations for indirect DM search (experimental side)
Are those really DM signals?
Are those really DM signals? We may check with other observables Today’s topic DM search using extragalactic gamma rays (including local galaxy distributions)
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the extragalactic region a). Inverse-Compton (IC) -rays in the extragalactic region
γ
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the a). Inverse-Compton (IC) -rays in the extragalactic region
γ
KI, Matsumoto, Moroi ’09 Profumo, Jeltema ’09
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the extragalactic region a). Inverse-Compton (IC) -rays in the extragalactic region
γ
Ando, KI ’15
Part I
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the extragalactic region a). Inverse-Compton (IC) -rays in the
γ
Cuoco, Xia, Regis, Branchini, Fornengo, Viel ’15
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the extragalactic region a). Inverse-Compton (IC) -rays in the extragalactic region
γ
Part II
Ando, KI ’16
Outline
DM and the extragalactic gamma rays
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the extragalactic region a). Inverse-Compton (IC) -rays in the extragalactic region
γ
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the a). Inverse-Compton (IC) -rays in the extragalactic region
γ
e±
a). Inverse-Compton (IC) -rays in the extragalactic region
γ
γ
e± e±
a). Inverse-Compton (IC) -rays in the extragalactic region
γ
γ
e± γ e±
γ
IC scattering
a). Inverse-Compton (IC) -rays in the extragalactic region
γ
exactly especially for decaying DM
anomalous positron or antiproton excess
KI, Matsumoto, Moroi ’09 Profumo, Jeltema ’09
a). Inverse-Compton (IC) -rays in the extragalactic region
γ
e±
γ
DM model
section
IC scattering
a). Inverse-Compton (IC) -rays in the extragalactic region
γ
flux
mass c- nM
/
lifetime
.Decay
nodes €zy
'*
→e±
→n¥§yn
CMBphotons
decay / annihilation modes$
rytiafnetwintiiefionoossseotiona
iain
→ e±→n¥§yn CMBphotons
Gamma-ray spectrum in various DM models
Ando, KI ’15
a). Inverse-Compton (IC) -rays in the extragalactic region
γ
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the extragalactic region a). Inverse-Compton (IC) -rays in the
γ
Fermi-LAT ’12
Tomborra, Ando, Murase ’14 Inoue ’11 Mauro, Calore, Donato, Ajello, Latronico ’14
They are determined due to recent updates in multi- frequency measurements of gamma rays
b). Astrophysical sources in the extragalactic region
Blazars
SFG
Fermi-LAT ’12
Tomborra, Ando, Murase ’14
b). Astrophysical sources in the extragalactic region
mAGN
Inoue ’11 Mauro, Calore, Donato, Ajello, Latronico ’14
b). Astrophysical sources in the extragalactic region
Ando, KI ’15
Blazars and SFGs well explain the observed gamma rays
b). Astrophysical sources in the extragalactic region
Ando, KI ’15
Constraints on DM scenarios
b). Astrophysical sources in the extragalactic region
Blazars and SFGs well explain the observed gamma rays
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the extragalactic region a). Inverse-Compton (IC) -rays in the extragalactic region
γ
Ando, KI ’15
For the anomalous positron For the TeV anomalous antiproton
Ando, KI ’15
Decaying DM scenarios to explain the anomalous positron or antiproton are partly excluded Decaying DM
Ando, KI ’15
Without astrophysical components
For the anomalous positron For the TeV anomalous antiproton
In the study, we considered that the gamma rays from the extragalactic region is
In the study, we considered that the gamma rays from the extragalactic region is
But due to the recent observational developments,
Important ingredients for our study: c). Tomographic cross-correlation using local galaxy distribution b). Astrophysical sources in the a). Inverse-Compton (IC) -rays in the
γ
DM and local galaxy distributions
Ingredients for further analysis:
Gamma rays are almost isotropic, but ..
c). Tomographic cross-correlation using local galaxy distribution
There’re anisotropies
Fermi-LAT ’12
c). Tomographic cross-correlation using local galaxy distribution
Ingredients for further analysis:
Fermi-LAT ’12
Ingredients for further analysis:
c). Tomographic cross-correlation using local galaxy distribution
2MRS ’11
Galaxy distribution
c). Tomographic cross-correlation using local galaxy distribution
c). Tomographic cross-correlation using local galaxy distribution
8000<v<9000km/s 7000<v<8000km/s 6000<v<7000km/s
We know the distance from each galaxy by its redshift
Ingredients for further analysis:
c). Tomographic cross-correlation using local galaxy distribution
Ingredients for further analysis:
c). Tomographic cross-correlation using local galaxy distribution
c). Tomographic cross-correlation using local galaxy distribution
Gamma rays are expected to trace galaxy distribution
c). Tomographic cross-correlation using local galaxy distribution
Cross correlation
c). Tomographic cross-correlation using local galaxy distribution
✖
c). Tomographic cross-correlation using local galaxy distribution
✖
2MRS,QSO, 2MASS,NVSS,MG,LRG
galaxy catalog:
c). Tomographic cross-correlation using local galaxy distribution
✖
Redshift distribution
c). Tomographic cross-correlation using local galaxy distribution
✖
Redshift distribution
Selecting a galaxy catalog You can get cross-correlation for corresponding redshift region
c). Tomographic cross-correlation using local galaxy distribution
✖
Catalog A
c). Tomographic cross-correlation using local galaxy distribution
✖
redshift Catalog A Catalog B Catalog C
Tomographic cross-correlation
c). Tomographic cross-correlation using local galaxy distribution
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i ˆ n ˆ n + θ
γ-ray flux galaxy distribution
δIγ = Iγ hIγi δΣg = Σg hΣgi
c). Tomographic cross-correlation using local galaxy distribution
Xia, Cuoco, Branchini, Viel ’15
Cross-correlation signal for < 1
✖
th.
✖
Compare both, then exclude the theory which deviates from
Xia, Cuoco, Branchini, Viel ’15
c). Tomographic cross-correlation using local galaxy distribution
✖
Catalog A Catalog B Catalog Cc). Tomographic cross-correlation using local galaxy distribution Theoretical calculation
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i
γ-ray flux galaxy distribution
δIγ = Iγ hIγi δΣg = Σg hΣgi
c). Tomographic cross-correlation using local galaxy distribution Theoretical calculation
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i
γ-ray flux galaxy distribution Σg = Z dχWg(z)ng(χˆ n, z) hngi
Wg(z) = d log Ng dz dz dχ
δIγ = Iγ hIγi δΣg = Σg hΣgi
hΣgi = 1
window function [dimensionless]
c). Tomographic cross-correlation using local galaxy distribution Theoretical calculation = DM + astro. sources
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i
γ-ray flux galaxy distribution
δIγ = Iγ hIγi δΣg = Σg hΣgi
Decaying DM
dΦdm
γdχ (Eγ, z) = 1 4π Ωdmρc mdmτdm 1 1 + z Qdm
γ (E0 γ, z) eτ(E0 γ,z)Idm
γ
= Z dχW dm
γ
(z) ρdm(χˆ n, z) hρdmi
(z) = Z dEγ dΦdm
γdχ (Eγ, z)
[ ] [ ]
cm−2s−1str−1
[ ]
cm−3s−1str−1 GeV−1cm−3s−1str−1
c). Tomographic cross-correlation using local galaxy distribution Theoretical calculation = DM + astro. sources
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i
γ-ray flux galaxy distribution
δIγ = Iγ hIγi δΣg = Σg hΣgi
Annihilating DM
dΦdm
γdχ (Eγ, z) = hσvi 8π ✓Ωdmρc mdm ◆2 (1 + z)3Qdm
γ (E0 γ, z) eτ(E0 γ,z)Idm
γ
= Z dχW dm
γ
(z) ρdm(χˆ n, z) hρdmi 2
W dm
γ(z) = Z dEγ dΦdm
γdχ (Eγ, z)
[ ]
cm−2s−1str−1
[ ]
cm−3s−1str−1
[ ]
GeV−1cm−3s−1str−1
c). Tomographic cross-correlation using local galaxy distribution Theoretical calculation = DM + astro. sources
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i
γ-ray flux galaxy distribution
δIγ = Iγ hIγi δΣg = Σg hΣgi
Decaying/Annihilating DM (source term)
Qdm
γ (E0 γ, z) = Qdm γpr(E0 γ, z) + Qdm γfsr(E0 γ, z) + Qdm γic (E0 γ, z)Qdm
γpr(E0 γ, z) + Qdm γfsr(E0 γ, z) = (1 + z)dNγdE (E0
γ)Qdm
γic(E0 γ, z) = cZ dEe dEγBG(1 + z)dσIC dE0
γ(E0
γ, Ee, EγBG)f BG γ(EγBG, z) Ye(Ee) bIC(Ee, z)
f BG
γ(EγBG, z) = f CMB
γ(EγBG, z) + f EBL
γ(EγBG, z) Ye(Ee) = X
I=e±Z ∞
EedE dNI dE (E)
GeV−1
[ ]
bIC(Ee, z) = Z dE0
γdEγBG(E0 γ − EγBG)dσICdE0
γ(E0
γ, Ee, EγBG)f BG γ(EγBG, z)
[ ]
GeV−1cm−3
[ ]
GeV s−1
c). Tomographic cross-correlation using local galaxy distribution Theoretical calculation = DM + astro. sources
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i
γ-ray flux galaxy distribution
δIγ = Iγ hIγi δΣg = Σg hΣgi
IX
γ =
Z dχW X
γ (z)
nX(χˆ n, z) hnXi
[ ]
cm−3s−1str−1 W X
γ (z) = χ2Z dLγ dnX
γ (Lγ, z)dLγ Fγ(Lγ, z) Gruppioni et al. ’13 Tamborra, Ando, Murase ’14 Ajello et al. ’15 Ackermann et al. ’12 Acero et al. ’15
: number flux of photons from a source with luminosity and redshift
[ ] [ ] [ ] c). Tomographic cross-correlation using local galaxy distribution Theoretical calculation = DM + astro. sources
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i
γ-ray flux galaxy distribution
δIγ = Iγ hIγi δΣg = Σg hΣgi
: luminosity function : luminosity
erg s−1 erg−1s cm−3
[ ]
cm−3s−1str−1 W X
γ (z) = χ2Z dLγ dnX
γ (Lγ, z)dLγ Fγ(Lγ, z)
dnX
γ (Lγ, z)dLγ
Fγ(Lγ, z) Lγ Lγ z cm−2s−1str−1
c). Tomographic cross-correlation using local galaxy distribution Theoretical calculation
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i
γ-ray flux galaxy distribution
δIγ = Iγ hIγi δΣg = Σg hΣgi
Energy [MeV] 3 10 4 10 5 10 6 10 ]Fornasa, Sánchez-Conde ’15
c). Tomographic cross-correlation using local galaxy distribution Advantages of tomography
クロス相関
暗黒物質は、 の領域からの寄与が相対 的に大きい の銀河とクロス相関をとることで、 の領域を選択的に取り出すことができる!
2MRS
ガンマ線も銀河も宇宙の大規模構造を トレースしている つのマップはある程度似ているはず
Ando ’14 Xia, Cuoco, Branchini, Viel ’15
c). Tomographic cross-correlation using local galaxy distribution Advantages of tomography
クロス相関
暗黒物質は、 の領域からの寄与が相対 的に大きい の銀河とクロス相関をとることで、 の領域を選択的に取り出すことができる!
2MRS
ガンマ線も銀河も宇宙の大規模構造を トレースしている つのマップはある程度似ているはず
Ando ’14 Xia, Cuoco, Branchini, Viel ’15
Blazars and SFGs are dominant in z > 0.1
c). Tomographic cross-correlation using local galaxy distribution Advantages of tomography
クロス相関
暗黒物質は、 の領域からの寄与が相対 的に大きい の銀河とクロス相関をとることで、 の領域を選択的に取り出すことができる!
2MRS
ガンマ線も銀河も宇宙の大規模構造を トレースしている つのマップはある程度似ているはず
Ando ’14 Xia, Cuoco, Branchini, Viel ’15
c). Tomographic cross-correlation using local galaxy distribution Theoretical calculation
Cγg(θ) = hδIγ(ˆ n) δΣg(ˆ n + θ)i
γ-ray flux galaxy distribution
δIγ = Iγ hIγi δΣg = Σg hΣgi
= X
`
2` + 1 4⇡ Cg
` P`(cos ✓)
Cg
`
= Z d 2 W(z)Wg(z)Pg ✓ k = ` , z ◆ cross-power spectrum between
γ
c). Tomographic cross-correlation using local galaxy distribution Cg
`
= Z d 2 W(z)Wg(z)Pg ✓ k = ` , z ◆
✖
Catalog A Catalog B Catalog Cc). Tomographic cross-correlation using local galaxy distribution Cg
`
= Z d 2 W(z)Wg(z)Pg ✓ k = ` , z ◆
✖
Catalog A Catalog B Catalog Ccross-correlation
c). Tomographic cross-correlation using local galaxy distribution Cg
`
= Z d 2 W(z)Wg(z)Pg ✓ k = ` , z ◆
✖
Catalog A Catalog B Catalog Ctomography
✖
th.
✖
Compare both, then exclude the theory which deviates from
Xia, Cuoco, Branchini, Viel ’15
c). Tomographic cross-correlation using local galaxy distribution
✖
Catalog A Catalog B Catalog Cトモグラフィー
使用可能な情報 より多くの赤方偏移の情報を使う → トモグラフィー
Catalog Redshift boundaries Ng per bin 2MRS (0.003, 0.1) 43500 2MRS-N2 (0.003, 0.027, 0.1) 21750 2MRS-N3 (0.003, 0.021, 0.035, 0.1) 14500 2MXSC (0.003, 0.3) 770000 2MXSC-N2 (0.003, 0.083, 0.3) 385000 2MXSC-N3 (0.003, 0.066, 0.10, 0.3) 257000 2MXSC-N4 (0.003, 0.058, 0.083, 0.11, 0.3) 193000 2MXSC-N5 (0.003, 0.052, 0.073, 0.093, 0.12, 0.3) 154000 2MXSC-N10 (0.003, 0.039, 0.052, 0.063, 0.073, 77000 0.083, 0.093, 0.10, 0.12, 0.14, 0.3)
Galaxy catalogs
c). Tomographic cross-correlation using local galaxy distribution
Ando ’14
The reported anomalous cosmic rays:
Decaying DM Annihilating DM The reported anomalous cosmic rays:
Decaying DM Annihilating DM
e+ Decaying DM (for the anomalous ) DM → νl±l⌥ (a). νµ±e⌥&νe±e⌥ (b). νµ±µ⌥&νe±µ⌥ Here we focus on three-body leptonic decay: (mainly ) e± (mainly ) µ±
e+
Ando, KI ’16
Decaying DM (for the anomalous )
(a). (b).
e+ Including astrophysical sources give ~10 times stronger constraints
Ando, KI ’16
Decaying DM (for the anomalous )
(a). (b).
The preferred regions are excluded e+
Best fit regions taken from Ibe et al.’14 Ando, KI ’16
Decaying DM (for the anomalous )
(a). (b).
Impacts of IC gamma rays
Ando, KI ’16
(Results without astro. comp.)
Results without IC (consistent with Regis et al. ’15)
(a). (b).
Ando, KI ’16
(a). (b).
IC gamma gives 1-2 orders of magnitude stronger constraints over TeV region Impacts of IC gamma rays
(Results without astro. comp.)
Ando, KI ’16
(a). (b).
IC gamma rays are crucial to constrain over TeV DM Impacts of IC gamma rays
(Results without astro. comp.)
IC gamma gives 1-2 orders of magnitude stronger constraints over TeV region
Ando, KI ’16
Decaying DM (for the anomalous ) The preferred regions are excluded (only by DM component) preferred region DM → W ±µ⌥ O(100) GeV ¯ p
Ando, KI ’16
Annihilating DM (for the anomalous ) The preferred regions are excluded (by including astro components) preferred region O(100) GeV ¯ p DM DM → W +W −
Annihilating DM (for the anomalous )
Ando, KI ’16
DM DM → b¯ b (consistent with
Cuoco et al. ’15)
Obtained constraints are similar to those given by dwarf galaxy O(1) GeV ¯ p
Best fit regions given by Cuoco et al. ’17 Ando, KI ’16
The motivated region is partly excluded DM DM → b¯ b Annihilating DM (for the anomalous ) O(1) GeV ¯ p
Although we’ve shown the constraints on DM models, Our goal is to find the DM! For the goal, a naive step to take next would be to consider
Although we’ve shown the constraints on DM models, Our goal is to find the DM! For the goal, a naive step to take next would be to consider
DM signal might be in PeV neutrino data
IceCube ’15
Although we’ve shown the constraints on DM models, Our goal is to find the DM! For the goal, a naive step to take next would be to consider
Cohen, Murase, Rodd, Safdi, Soreq ’17 Kalashev, Kuznetsov ’16 Kachelriess, Kalashev, Kuznetsov ’18 Dudas, Gherghetta, Kaneta, Mambrini, Olive ’18
People are getting interested in heavier DM
Beyond High Energy Cosmic Rays
High Energy Cosmic Rays (HECRs) Very High Energy Cosmic Rays (VHECRs) Ultra High Energy Cosmic Rays (UHECRs) Extremely High Energy Cosmic Rays (EHECRs)
TeV PeV EeV ZeV
Fonseca ’03
Beyond High Energy Cosmic Rays
High Energy Cosmic Rays (HECRs) Very High Energy Cosmic Rays (VHECRs) Ultra High Energy Cosmic Rays (UHECRs) Extremely High Energy Cosmic Rays (EHECRs)
TeV PeV EeV ZeV
Fonseca ’03 IceCube ’15
Beyond High Energy Cosmic Rays
High Energy Cosmic Rays (HECRs) Very High Energy Cosmic Rays (VHECRs) Ultra High Energy Cosmic Rays (UHECRs) Extremely High Energy Cosmic Rays (EHECRs)
TeV PeV EeV ZeV
Fonseca ’03 AGASA ’03
6 10 3 2 3 10 10 10 19 20 10 10 10 23 24 25 26 J(E) E [m sec sr eV ] 3 −2 −1 −1 2 Energy [eV] Uniform sources C1 C2 C3 C4 C5 C6 C7More will be by Pierre Auger Observatory
Propagation of UHECR
Initial state Target field Process Secondaries Nuclei CBR Pair production (Bethe-Heitler) e± Nuclei CBR Photo-pion production p, n, ⌫, e±, Nuclei CBR Photodisintegration p, n, d, t, 3He, ↵, * Nuclei CBR Elastic scattering*
– Nuclear decay p, n, ⌫, e±, * Photons CBR Pair production* (Breit-Wheeler) e± Photons CBR Double pair production* e± Electrons CBR Triplet pair production* e± Electrons CBR Inverse Compton scattering*
B-field Synchrotron radiation*
Propagation of UHECR
Initial state Target field Process Secondaries Nuclei CBR Pair production (Bethe-Heitler) e± Nuclei CBR Photo-pion production p, n, ⌫, e±, Nuclei CBR Photodisintegration p, n, d, t, 3He, ↵, * Nuclei CBR Elastic scattering*
– Nuclear decay p, n, ⌫, e±, * Photons CBR Pair production* (Breit-Wheeler) e± Photons CBR Double pair production* e± Electrons CBR Triplet pair production* e± Electrons CBR Inverse Compton scattering*
B-field Synchrotron radiation*
Propagation of UHECR nuclei
N + γBG → N + π
A ZX + γBG → A ZX + e+ + e−
Pair production Photo-pion production
xloss(E) = E dE/dx
Stanev, Engel, Mücke, Protheroe, Rachen ’00
Proton
Propagation of UHECR
Initial state Target field Process Secondaries Nuclei CBR Pair production (Bethe-Heitler) e± Nuclei CBR Photo-pion production p, n, ⌫, e±, Nuclei CBR Photodisintegration p, n, d, t, 3He, ↵, * Nuclei CBR Elastic scattering*
– Nuclear decay p, n, ⌫, e±, * Photons CBR Pair production* (Breit-Wheeler) e± Photons CBR Double pair production* e± Electrons CBR Triplet pair production* e± Electrons CBR Inverse Compton scattering*
B-field Synchrotron radiation*
EW cascade
Propagation of UHECR EM particles
e + γBG → e + e+ + e− γ + γBG → e+ + e− γ + γBG → e+ + e− + e+ + e− e + γBG → e + γ
Heiter, Kuempel, Walz, Erdmann ’17
CRpropa 3
A public astrophysical simulation framework for propagating extraterrestrial ultra-high energy particles
Batista, Dundovic, Erdmann, Kampert, Kuempel, Müller, Sigl, Vliet, Walz, Winchen ’16
Energy [eV]
13 10 14 10 15 10 16 10 17 10 18 10 19 10 20 10 21 10 22 10[a.u.]
2dN/dE E
16 10 17 10 18 10 19 10 20 10 21 10 22 10 23 10 24 10 25 10 DINT only CRPropa 3 + DINT Electron Pair Production Photo Pion Production Nuclear DecayProton injection
dN/dE ∝ E−1
at 3-1000 Mpc
CRpropa 3
A public astrophysical simulation framework for propagating extraterrestrial ultra-high energy particles
Batista, Dundovic, Erdmann, Kampert, Kuempel, Müller, Sigl, Vliet, Walz, Winchen ’16
Energy [eV]
13 10 14 10 15 10 16 10 17 10 18 10 19 10 20 10 21 10 22 10[a.u.]
2dN/dE E
16 10 17 10 18 10 19 10 20 10 21 10 22 10 23 10 24 10 25 10 DINT only CRPropa 3 + DINT Electron Pair Production Photo Pion Production Nuclear DecayProton injection
dN/dE ∝ E−1
at 3-1000 Mpc CRs from DM can be simulated similarly
in preparation with
CRpropa 3
A public astrophysical simulation framework for propagating extraterrestrial ultra-high energy particles
Batista, Dundovic, Erdmann, Kampert, Kuempel, Müller, Sigl, Vliet, Walz, Winchen ’16
Energy [eV]
13 10 14 10 15 10 16 10 17 10 18 10 19 10 20 10 21 10 22 10[a.u.]
2dN/dE E
16 10 17 10 18 10 19 10 20 10 21 10 22 10 23 10 24 10 25 10 DINT only CRPropa 3 + DINT Electron Pair Production Photo Pion Production Nuclear DecayProton injection
dN/dE ∝ E−1
at 3-1000 Mpc CRs from DM can be simulated similarly
We have studied DM using extragalactic gamma rays and local galaxy distribution
excluded
____ _ is partly excluded e+ O(100 GeV) ¯ p O(1 GeV) ¯ p