a state of the art calculation of atmospheric lepton
play

A state-of-the-art calculation of atmospheric lepton fluxes 35 th - PowerPoint PPT Presentation

A state-of-the-art calculation of atmospheric lepton fluxes 35 th International Cosmic Ray Conference 2017, Bexco, Busan, South Korea Anatoli Fedynitch In collaboration with: DESY Zeuthen Hans Dembinski (Bartol, MPIK) Felix Riehn (KIT, LIP)


  1. A state-of-the-art calculation of atmospheric lepton fluxes 35 th International Cosmic Ray Conference 2017, Bexco, Busan, South Korea Anatoli Fedynitch In collaboration with: DESY Zeuthen Hans Dembinski (Bartol, MPIK) Felix Riehn (KIT, LIP) Ralph Engel (KIT) Tom K. Gaisser (Bartol) Todor Stanev (Bartol)

  2. Importance of high precision atmospheric lepton flux prediction > Conventional and prompt leptons are backgrounds for astrophysical neutrino measurements > Estimation of spectral index and normalization depends on the shape of the atmospheric expectation > Accuracy of calculations at lower energies limits neutrino properties measurements with atmospheric neutrinos Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 2

  3. Importance of high precision atmospheric lepton flux prediction IceCube Collaboration, ICRC 2015 > Conventional and prompt leptons are backgrounds for astrophysical neutrino measurements > Estimation of spectral index and normalization depends on the shape of the atmospheric expectation > Accuracy of calculations at lower energies limits neutrino properties measurements with atmospheric neutrinos Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 3 IceCube, PRD 91, 072004 (2015)

  4. Transport equations (hadronic cascade equations) System of non-linear PDE for each particle species h (~62 x #E-bins) : Interactions with air Decays Energy losses (radiative) Re-injection from interactions Re-injection from decays Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 4

  5. Transport equations (hadronic cascade equations) System of non-linear PDE for each particle species h (~62 x #E-bins) : cosmic ray physics Interactions with air Decays atmospheric physics Energy losses (radiative) Re-injection from interactions Re-injection from decays particle physics Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 5

  6. Our tools: MCEq - open-source cascade equation solver > Simultaneous integration of several thousands kinetic equations > Energy range 1 (30) GeV – 10 11 GeV > All models included > High optimization: multi- core, GPU, … (BLAS, MKL, CUDA) (~milli-seconds) > MIT licensed @ https://github.com/afedynitch/MCEq CORSIKA: A. Fedynitch, J. Becker Tjus and P. Desiati, PRD 2012 MCEq: A. Fedynitch, R. Engel, T. K. Gaisser, F. Riehn and S. Todor. PoS ICRC 2015, 1129 Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 6

  7. Models of the comic ray spectrum at the top of the atmosphere > None of the shown models extends over entire energy range and provides error estimates > GST-X and HXa are quite extreme assumptions for UHECR HXa : T. K. Gaisser, Astropart. Phys. 35 (2012) 801 – 806 GH : T. K. Gaisser and M. Honda, Ann. Rev. Nucl. Part. Sci. 52 (2002) 153 – 199 GST-X : T. K. Gaisser, T. Stanev, and S. Tilav, Front. Phys.(Beijing) 8 (2013) 748 – 758 Update of GH (not shown) : J. Evans, D. Porzio et al., 1612.03219 Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 7

  8. GSF flux model (talk by H. Dembinski 19 th July 14:45) Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 8

  9. GSF flux model (talk by H. Dembinski 19 th July 14:45) > Wavy feature at lower energies are due to the hardening of proton and He spectra (see also a recent update of GH model by Evans et al. 1612.03219) > Increase of error around 10 TeV because of the gap between direct and indirect exp. > Higher flux at the knee and harder specrtal index between knee and ankle > Latter effect comes from the lighter composition at the knee as in other models > Mainly driven by KASKADE Gr. and latest data from IceTop and TUNKA Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 9

  10. SIBYLL 2.3c with charm and improved kaon distributions More details in F. Riehn et al., PoS(ICRC2017)301 Talk by Ralph Engel right now  Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 10

  11. Implementation of “Barr et al.” in a numerical scheme > Compute partial derivatives wrt. phase-space Original scheme regions, i.e. > No correlations between phase-space regions (as in Barr et al.) or add. correlations Elements of Jacobian (numerical) Extension for charm Charm Error propagation CH_B + CH_A Energy dependence G. D. Barr, S. Robbins, T. K. Gaisser, and T. Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 11 Stanev, Phys. Rev. D 74, 094009 (2006)

  12. Implementation of “Barr et al.” in a numerical scheme > Compute partial derivatives wrt. phase-space Original scheme regions, i.e. > No correlations between phase-space regions (as in Barr et al.) or add. correlations Elements of Jacobian (numerical) Extension for charm Charm Charm Error propagation CH_B + 70 % + max. 25 % CH_A 10 % Energy dependence Energy dependence G. D. Barr, S. Robbins, T. K. Gaisser, and T. Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 12 Stanev, Phys. Rev. D 74, 094009 (2006)

  13. Uncertainties of lepton fluxes Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 13

  14. Uncertainties of lepton fluxes > Similar error band for each interaction model > Error from Barr et al. scheme approximately universal for all current models > Hadronic interaction models produce compatible results Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 14

  15. Uncertainties of lepton ratios Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 15

  16. Uncertainties of lepton ratios > Only very few benchmark measurements available > Ratios are more sensitive to the leading particle effect > SIBYLL 2.1 and 2.3 known to have problems with K+ > Newer models all compatible Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 16

  17. Prompt neutrinos from decays of charmed mesons > SIBYLL 2.3c is the only full MC model > Compatible with LHC data and IceCube limit > New CR flux model (GSF) changes situation a bit > Uncertainties from QCD theory are large and calculations are compatible within bands IceCube: Astrophys.J. 833 (2016) GMS: Garzelli et al., JHEP 1510 (2015) 115 BERSS: Bhattacharya et al. JHEP 2015: 110 Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 17

  18. Prompt neutrinos from decays of charmed mesons > SIBYLL 2.3c is the only full MC model > Compatible with LHC data and IceCube limit > New CR flux model (GSF) changes situation a bit > Uncertainties from QCD theory are large and calculations are compatible within bands IceCube: Astrophys.J. 833 (2016) GMS: Garzelli et al., JHEP 1510 (2015) 115 BERSS: Bhattacharya et al. JHEP 2015: 110 Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 18

  19. Conclusions > We developed a high -performance open-source code for atmospheric lepton calculations (Updated version of MCEq soon @ https://github.com/afedynitch/MCEq) > Uncertainties are dominated by hadronic interactions and cosmic ray flux errors > The CR flux errors are quantified using a new global fit to experimental data (GSF); the model changes the expectation of the very high energy neutrinos > Hadronic errors are quantified using a scheme after Barr et al., which can be used as a starting point for further studies; the error bands weakly depend on the initial choice of the interaction model > The high performance of the code together with the error propagation scheme can be used to constrain hadronic interactions in a phase-space inaccessible to colliders Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 19

  20. Hadronic uncertainties: current state of the art > “ Uncertainties in atmospheric neutrino fluxes ”, G. D. Barr, S. Robbins, T. K. Gaisser, and T. Stanev, Phys. Rev. D 74, 094009 (2006) > Cut phase-space in regions/slices in E lab and x lab and assign uncertainty to each slice (uncorrelated) > Draw-back 1 : Uncorrelated uncertainty assigned by hand and judged only from availability of experimental data (not how well TARGET described it) > Draw-back 2 : The “central value” is assumed to be TARGET. Scheme doesn’t tell anything about ”best estimate”. Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 20

  21. Hadronic uncertainties computed with different interaction models Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 21

  22. Hadronic uncertainties computed with different interaction models 1. If models with known problems are excluded, then errors are ~ universal = approx. independent on the choice of interaction model Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 22 “Bracketing” with different models can be a good estimate, sometimes 2.

  23. Hadronic uncertainties below 200 GeV Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 23

  24. Hadronic uncertainties for high energy leptons Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 24

  25. Relevant particles for lepton production Muon neutrinos Electron neutrinos charm decay pion decay kaon decay charm decay muon decay kaon decay Anatoli Fedynitch | 35 th ICRC, 2017, Busan, South Korea | 2017/07/10 | Page 25

Recommend


More recommend