Calculation of Atmospheric Neutrino Flux Based on AMS02 Observation - - PowerPoint PPT Presentation

calculation of atmospheric neutrino flux based on ams02
SMART_READER_LITE
LIVE PREVIEW

Calculation of Atmospheric Neutrino Flux Based on AMS02 Observation - - PowerPoint PPT Presentation

Calculation of Atmospheric Neutrino Flux Based on AMS02 Observation July 13, 2017 M. Honda @ ICRC2017 1. Cosmic Ray Observation of AMS02 and Some Others 2. Cosmic Spectra Model 3. Muon Calibration of Hadronic Interaction Model 4. Calculation


slide-1
SLIDE 1

Calculation of Atmospheric Neutrino Flux Based

  • n AMS02 Observation

July 13, 2017 M. Honda @ ICRC2017

  • 1. Cosmic Ray Observation of AMS02 and Some Others
  • 2. Cosmic Spectra Model
  • 3. Muon Calibration of Hadronic Interaction Model
  • 4. Calculation of Atmospheric Neutrino Flux
  • 5. Comparison with Previous Calculation and Estimation of Errors.
  • 6. Summary
slide-2
SLIDE 2

Primary Cosmic Ray Model and referred data (2004)

Other chemical compositions are also considered in the calculation, but they give small contributions.

slide-3
SLIDE 3

Balloon Borne (BESS) Satellite (ISS, AMS02)

Direct Observation

slide-4
SLIDE 4

Recent Cosmic Ray observation and Available High Energy data

slide-5
SLIDE 5

Proton Closeup

slide-6
SLIDE 6

Helium Closeup

slide-7
SLIDE 7

Cosmic Ray Spectra Model Based on AMS02 Observation

Looking forward to hearing from CALET and ISS-CREAM

slide-8
SLIDE 8

Solar Modulation and Neutron Monitor

slide-9
SLIDE 9

Solar Modulated Flux at Fixed Energy

slide-10
SLIDE 10

Practical Formula for Solar Modulation

Where a is from right figure, and N is the Count of Newark Neutron Monitor.

slide-11
SLIDE 11

Muon Calibration

Primary cosmic ray flux Interaction model Calculation scheme (include rigidity cutoff) Atmosphere model

slide-12
SLIDE 12

Muon Calibration of Interaction Model

Quick 3D Calculation for Muon flux. As the muon flux is a “local quantity” (γ ct 〜 60km at10 GeV ), We can calculate it in a quick calculation method:

  • 1. Inject cosmic rays just above the observation point,
  • 2. Analyze all muons reach the surface of Earth.
slide-13
SLIDE 13

Comparison of Quick 3D and Full 3D calculations

Full 3D Quick 3D

μ

+

μ

This method works above 0.2 GeV/c.

slide-14
SLIDE 14

Muon Calibration of inclusive DPMJET-III Data are larger by ~15% Data are larger by ~0.05 Data are smaller by ~0.05 ==> DPMJET-III Should be Modifjed ~15% scatter ?

slide-15
SLIDE 15

Comparison AFTER the modification

slide-16
SLIDE 16

JAM + Modified DPMJET-II vs Muons at the Balloon altitude ( HKKM2011) Good agreement ! Use DPMJET-III above 32 GeV and JAM below 32 GeV

slide-17
SLIDE 17
slide-18
SLIDE 18

With Cosmic ray spectra model based

  • n AMS02 observation

After muon re-calibration Before muon re-calibration

slide-19
SLIDE 19

Simulation Sphere (Rs = 10 x Re)

Full 3D-Calculation

Injection Sphere (Re +100lm)

Virtual Detector

The neutrino flux is calculated from the number of neutrinos path through with virtual detector correction.

Cosmic Rays are sampled and injected here

Re = 6378km

Cosmic ray go out this sphere are discarded. Cosmic rays go beyond are pass the rigidity cutofg test

slide-20
SLIDE 20

Atmosphere Model

Air density comparison with MSISE90 US-starndard'76 may be used as the global approximation of the Atmosphere.

slide-21
SLIDE 21

Kamioka INO site South Pole Near North Pole (Physalmi)

Atmosphere model (NRLMSISE-00) and seasonal variations

slide-22
SLIDE 22

IGRF10 Geomagnetic Horizontal Field Strength

slide-23
SLIDE 23

Based On AMS02 Obervation (Preliminary)

slide-24
SLIDE 24

Estimated Error in Atmospheric ν-flux Calculation (HKKMS07)   -observation error + Residual of reconstruction K  air

Kaon production uncertainty Mean free path (interaction crossection) uncertainty Atmosphere density profule uncertainty Possible Error with JAM (HKKM11)

slide-25
SLIDE 25

Summary

  • According to the new cosmic ray data provided by AMS02 and some
  • thers, We have renewed the Cosmic Ray Spectra Model including the

practical formula for the Solar Modulations.

  • The Muon Calibration of Hadronic Interaction Model was repeated

again, with the newly constructed Cosmic Ray Spectra Model.

  • The atmospheric neutrino flux calcuoated with New Cosmic Ray

Spectra Model and the Muon-Calibrated Hadronic Interaction Model wit it is generally similar to the previously calculated one, considering the difference of the Muon-Calibration Policy.

  • However, a ~ 5% deficit of the flux at ~ 10GeV is seen in the newly

calculated atmospheric neutrino flux irrespective of arriving directions

  • r kind of neutrinos.
slide-26
SLIDE 26

Comparison with Neutrino Data

slide-27
SLIDE 27

From K.Okumura in ICRC2015

slide-28
SLIDE 28

Best fit corresponds to 62 %

  • f the predicted variations

Solar Modulation of Atmospheric Neutrinos

From PHD thesis of

  • E. Richard

Poor Statistics

Picked up mainly the forbush decrease ?

slide-29
SLIDE 29

Observed Azimuthal Variation of flux (from PHD thesis of E.Richard)

νe

Zenith Angle Binned All Energies Energy Binned All Azimuth angles

slide-30
SLIDE 30

νμ

Observed Azimuthal Variation of flux (from PHD thesis of E.Richard)

Zenith Angle Binned All Energies Energy Binned All Azimuth angles

air

slide-31
SLIDE 31

Comparison in [Flux/depth]

DPMJET-III show the best agreement

slide-32
SLIDE 32

Pν=0.32

At GeV/c

π

± → νμ(̄

νμ)+ μ

±

p+ N(air)→ π

±+ X

μ

± → ̄

νμ(νμ)+ νe(̄ νe)+ e

±

Phase Space Study

slide-33
SLIDE 33

Pμ=0.32

At GeV/c

 pN air X

Size of Phase space to give the variation

slide-34
SLIDE 34

e  e   

Vertical neutrino flux

slide-35
SLIDE 35

e  e   

Horizontal neutrino flux

slide-36
SLIDE 36

Horizontal enhancement of neutrino fmux (Battistoni et al. Astropart. Phys 1999)

slide-37
SLIDE 37

Interpretation of horizontal enhancement Longer integration length in the neutrino production zone for horizontal directions