Issue |
A&A
Volume 636, April 2020
|
|
---|---|---|
Article Number | C3 | |
Number of page(s) | 3 | |
Section | Astrophysical processes | |
DOI | https://doi.org/10.1051/0004-6361/201732392e | |
Published online | 23 April 2020 |
Ultra-high-energy cosmic rays and neutrinos from tidal disruptions by massive black holes (Corrigendum)
1
Sorbonne Université, CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France
e-mail: claire.guepin@iap.fr
2
Department of Astronomy, University of Maryland, College Park, MD 20742, USA
3
Joint Space-Science Institute, University of Maryland, College Park, MD 20742, USA
4
SISSA, Via Bonomea 265, 34136 Trieste, Italy
5
INFN Sezione di Trieste, Trieste, Italy
6
Institute for Fundamental Physics of the Universe, Via Beirut 2, 34014 Trieste, Italy
7
Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA
8
Department of Physics, The Pennsylvania State University, University Park, PA 16802, USA
9
Department of Astronomy & Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA
10
Center for Particle and Gravitational Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA
11
Yukawa Institute for Theoretical Physics, Kyoto 606-8502, Japan
Key words: astroparticle physics / neutrinos / errata, addenda
We noticed a normalization error in our calculations and a bug in our intergalactic propagation code, which impacted our predictions of the ultra-high-energy cosmic ray (UHECR) and high-energy (HE) neutrino fluxes from jetted tidal disruption events (TDE; Guépin et al. 2018). Correcting these errors led to sizable changes in the cosmic rays and neutrino fluxes, but our conclusions remain the same.
The corrected UHECR and HE neutrino spectra from single TDEs, in the high and medium states, are illustrated in Figs. 1 and 2. The high state of the UHECR spectrum is slightly modified as the maximum energy of injected particles is determined by the competition between the acceleration and the dynamical timescales, and not between the acceleration and the cooling timescales. As stated above, the HE neutrino fluxes are one order of magnitude higher that the ones derived previously. Thus we show that the HE neutrino flux from one source at dL, 2 = 50 Mpc reaches the IceCube detection limit. We note that considering the observed event rate density of jetted TDE is ∼0.03 Gpc−3 yr−1, the rate of TDEs at distances smaller than 50 Mpc is 4 × 10−6 yr−1, which is extremely low. The sensitivity limit of GRAND has also been updated (GRAND Collaboration 2020) as GRAND cannot observe a TDE during its entire duration: The source moves out of the field of view of the detector because of the Earth’s rotation. Thus the GRAND sensitivity limit for long duration bursts (lasting for more than a few hours) is reduced when compared to the sensitivity limit for short duration bursts, which was previously shown.
![]() |
Fig. 1. Cosmic-ray spectra for one source with pure iron injection with spectral index α = 1.8, photon field with |
![]() |
Fig. 2. Neutrino spectra for three flavors for one source with same characteristics as in Fig. 1. We show the total spectra (in eV cm−2) for a high state SED with Lpk = 1047.5 erg s−1 (left) and a medium state SED with Lpk = 1046 erg s−1 (right). We consider two different distances dL, 1 = 1.88 Gpc (dark blue) and dL, 2 = 50 Mpc (light blue). The IceCube and projected GRAND (GRAND Collaboration 2020) sensitivities are also shown (dashed black and green lines). For the IceCube sensitivities, we show two cases depending on the declination: 0° < δ < 30° (most favorable case, thin line) and 30° < δ < 60° (Swift J1644+57 case, thick line) (Aartsen et al. 2014). |
The corrected diffuse UHECR and HE neutrino spectra, as well as the corresponding UHECR composition, are illustrated in Figs. 3 and 4. The bug in the intergalactic propagation code only affected the spectrum of high-energy protons. A large number of protons pile up below 1019 eV. It affects the composition of UHECR, which better reproduces, than our precious results, the standard deviation of the atmospheric depth of air shower maximum development σ(Xmax) observed by the Pierre Auger observatory below 1019 eV.
![]() |
Fig. 3. Diffuse cosmic-ray (left) and neutrino (right) spectra from a population of jetted TDEs obtained for an injection of 70% Si and 30% Fe, with spectral index α = 1.5, ξCRfs = 2.6 × 10−3 and source evolution derived in this work, with maximum bolometric luminosity in the population Lmax = 6.8 × 1046 erg s−1 (in the medium state). Left: total spectrum (black) and its composition. We superimposed data from the Auger experiment (black dots, The Pierre Auger Collaboration 2015) and from the Telescope Array experiment (gray dots, Fukushima 2015) for which only statistical uncertainties are shown. Right: diffuse neutrino flux for three flavors (black), the diffuse neutrino flux measured by the IceCube experiment (IceCube Collaboration 2017), and the projected limits for GRAND (GRAND Collaboration 2020), ARA/ARIANNA (Allison et al. 2015; Barwick et al. 2015), and POEMMA (Neronov et al. 2017). |
![]() |
Fig. 4. Mean and standard deviation of Xmax for the spectrum shown in Fig. 3 (pale gray band). We also show Auger measurements (Aab et al. 2014) with uncertainties (black dots) and simulation results for pure proton injection (red band) and pure iron injection (blue band). The bands are obtained by accounting for hadronic model uncertainties (EPOS-LHC, SIBYLL 2.1 or QGSJET II-04). |
The fraction fs of the total population of TDE required to reproduce the data from the Pierre Auger observatory is modified. The diffuse UHECR spectrum is calculated for an injection of 70% Si and 30% Fe, a spectral index α = 1.5, an acceleration efficiency ηacc = 0.1, a TDE source evolution, and a fraction ξCRfs = 2.6 × 10−3 of the local event rate density ṅtde,0 = 1.5 × 102 Gpc−3 yr−1. This rate is computed from the TDE rate per galaxy obtained in the case of core galaxies. The population fraction approximately corresponds to the rate density ∼0.4 Gpc−3 yr−1. Moreover, the predicted diffuse HE neutrino spectrum is lower than in our previous estimates, as our estimate of fs has decreased, and is about two orders of magnitude below the spectrum of HE neutrinos detected by IceCube. The detection of very-high-energy neutrinos from TDE by detectors, such as GRAND, also seems disfavored.
From these results, accounting for time variability of multi-messenger emissions seems to be important as the bulk of the different messengers are not typically emitted during the same phase. In this example, most of the UHECR contributing to the diffuse flux, allowing the Auger data to be reproduced, are produced during the medium state, whereas most of the HE neutrinos are produced during the high state. In order to reproduce the Auger data, a model accounting only for the high state would require a higher baryon loading or a higher TDE event rate, which would subsequently increase the associated HE neutrino flux.
References
- Aab, A., Abreu, P., Aglietta, M., et al. 2014, Phys. Rev. D, 90, 122005 [NASA ADS] [CrossRef] [Google Scholar]
- Aartsen, M. G., Ackermann, M., Adams, J., et al. 2014, ApJ, 796, 109 [Google Scholar]
- Allison, P., Auffenberg, J., Bard, R., et al. 2015, Astropart. Phys., 70, 62 [NASA ADS] [CrossRef] [Google Scholar]
- Barwick, S. W., Berg, E. C., Besson, D. Z., et al. 2015, Astropart. Phys., 70, 12 [NASA ADS] [CrossRef] [Google Scholar]
- Fukushima, M. 2015, Eur. Phys. J. Web Conf., 99, 04004 [CrossRef] [Google Scholar]
- GRAND Collaboration (Alvarez-Muniz, J., et al.) 2020, Sci. China Phys. Mech. Astron., 63, 219501 [CrossRef] [Google Scholar]
- Guépin, C., Kotera, K., Barausse, E., Fang, K., & Murase, K. 2018, A&A, 616, A179 [CrossRef] [EDP Sciences] [Google Scholar]
- IceCube Collaboration (Aartsen, M. G., et al.) 2017, ArXiv e-prints [arXiv:1710.01191] [Google Scholar]
- Neronov, A., Semikoz, D. V., Anchordoqui, L. A., Adams, J. H., & Olinto, A. V. 2017, Phys. Rev. D, 95, 023004 [NASA ADS] [CrossRef] [Google Scholar]
- The Pierre Auger Collaboration (Aab, A., et al.) 2015, ArXiv e-prints [arXiv:1509.03732] [Google Scholar]
© ESO 2020
All Figures
![]() |
Fig. 1. Cosmic-ray spectra for one source with pure iron injection with spectral index α = 1.8, photon field with |
In the text |
![]() |
Fig. 2. Neutrino spectra for three flavors for one source with same characteristics as in Fig. 1. We show the total spectra (in eV cm−2) for a high state SED with Lpk = 1047.5 erg s−1 (left) and a medium state SED with Lpk = 1046 erg s−1 (right). We consider two different distances dL, 1 = 1.88 Gpc (dark blue) and dL, 2 = 50 Mpc (light blue). The IceCube and projected GRAND (GRAND Collaboration 2020) sensitivities are also shown (dashed black and green lines). For the IceCube sensitivities, we show two cases depending on the declination: 0° < δ < 30° (most favorable case, thin line) and 30° < δ < 60° (Swift J1644+57 case, thick line) (Aartsen et al. 2014). |
In the text |
![]() |
Fig. 3. Diffuse cosmic-ray (left) and neutrino (right) spectra from a population of jetted TDEs obtained for an injection of 70% Si and 30% Fe, with spectral index α = 1.5, ξCRfs = 2.6 × 10−3 and source evolution derived in this work, with maximum bolometric luminosity in the population Lmax = 6.8 × 1046 erg s−1 (in the medium state). Left: total spectrum (black) and its composition. We superimposed data from the Auger experiment (black dots, The Pierre Auger Collaboration 2015) and from the Telescope Array experiment (gray dots, Fukushima 2015) for which only statistical uncertainties are shown. Right: diffuse neutrino flux for three flavors (black), the diffuse neutrino flux measured by the IceCube experiment (IceCube Collaboration 2017), and the projected limits for GRAND (GRAND Collaboration 2020), ARA/ARIANNA (Allison et al. 2015; Barwick et al. 2015), and POEMMA (Neronov et al. 2017). |
In the text |
![]() |
Fig. 4. Mean and standard deviation of Xmax for the spectrum shown in Fig. 3 (pale gray band). We also show Auger measurements (Aab et al. 2014) with uncertainties (black dots) and simulation results for pure proton injection (red band) and pure iron injection (blue band). The bands are obtained by accounting for hadronic model uncertainties (EPOS-LHC, SIBYLL 2.1 or QGSJET II-04). |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.