Cosmogenic neutrino fluxes under the effect of active-sterile secret interactions Damiano Fiorillo, Gennaro Miele, Stefano Morisi Dipartimento di Fisica "Ettore Pancini", Universit`adegli studi di Napoli Federico II, Complesso Univ. Monte S. Angelo, I-80126 Napoli, Italy and INFN - Sezione di Napoli, Complesso Univ. Monte S. Angelo, I-80126 Napoli, Italy Ninetta Saviano INFN - Sezione di Napoli, Complesso Univ. Monte S. Angelo, I-80126 Napoli, Italy (Dated: May 13, 2020) Ultra High Energy cosmogenic neutrinos may represent a unique opportunity to unveil possible new physics interactions once restricted to the neutrino sector only. In the present paper we study the observable effects of a secret active-sterile interactions, mediated by a pseudoscalar, on the expected flux of cosmogenic neutrinos. The results show that for masses of sterile neutrinos and pseudoscalars of hundreds MeV, necessary to evade cosmological, astrophysical and elementary particle constraints, the presence of such new interactions can significantly change the energy spectrum of cosmogenic neutrinos at Earth in the energy range from PeV to ZeV. Interestingly, the distortion of the spectrum results to be detectable at apparatus like GRAND or ARIANNA if the mediator mass is around 250 MeV and the UHECRs are dominated by the proton component. Larger mediator masses or a chemical composition of UHECRs dominated by heavier nuclei would require much larger cosmic rays apparatus which might be available in future. Keywords: Ultra High Energy neutrinos, cosmogenic neutrinos, secret interactions I. INTRODUCTION GRAND [34], ARIANNA [35], ARA [36] and JEM-EUSO [37]. Unfortunately, a possible additional difficulty lies in The neutrino sector still represents a partially un- the fact that cosmogenic neutrino fluxes have a strong de- known territory. Fundamental questions like the nature pendence on the chemical composition of cosmic rays, as of neutrinos (Dirac or Majorana) or the possible con- shown for instance in [38]. This is particularly important nection of their small masses with physics Beyond the since recent results suggest that the chemical composi- Standard Model (BSM) can represent possible windows tion is actually mixed, containing significant amounts of on new physics. In the last decade increasing attention heavier nuclei rather than simply protons [39], and unfor- has been devoted to high energy astrophysical neutrinos, tunately for heavier nuclei one expects a suppression in after the observation of the first events at the IceCube the neutrino production. For this reason, in this work we detector [1]. These astrophysical fluxes are of extreme have analyzed two benchmark scenarios for the chemical relevance for the neutrino sector since they provide a composition of the UHECRs, namely the cases in which powerful tool of investigation for Beyond Standard Model the dominant component is either protons or Helium nu- physics, such as sterile neutrinos, Lorentz violations and clei. non standard model interactions. In this context, cosmo- In view of the above considerations a future detection genic neutrinos, that are mainly expected to have ener- of cosmogenic neutrinos would allow us to infer crucial gies up to 1012 GeV, could represent a unique opportu- astrophysical informations concerning UHECRs chem- nity. ical composition, interactions and origin. However, at Cosmogenic neutrinos are produced by the photo- the same time the shape of the cosmogenic neutrino hadronic interactions of Ultra High Energy Cosmic rays spectrum, that could be strongly distorted by the (UHECRs), whose precise chemical composition and ori- presence of new physics, could unveil the presence of gin is still unknown, with the Cosmic Microwave Back- these new interactions at least in the neutrino sector. ground (CMB) [2]. The same process is also responsible In particular, any kind of BSM interactions involving for a depletion in the flux of UHECRs, which is known in active neutrinos ν could modify the expected spectrum literature as the GZK cutoff [3, 4]. For this reason cos- and rate of cosmogenic neutrinos at Earth. Examples of arXiv:2002.10125v2 [hep-ph] 12 May 2020 mogenic neutrinos are also known as GZK-neutrinos and such new couplings discussed in literature are: they have been extensely studied in a number of works, i) Non Standard Interactions (NSI) for active neutrinos see for instance [5{30]. νν ! ff [40{51], where f denotes quarks or charged The desirable observation of high-energy cosmogenic leptons; neutrinos would be particularly relevant in order to de- ii)Secret interactions (SI) mediated by some new boson termine the origin of UHECRs. Several cosmic rays appa- (scalar or vector), and just restricted to the active ratus, like HiRes [31] and the Pierre Auger Observatory neutrino sector νν ! νν [24, 52{61]; (PAO) [32, 33] for example, have already tried to per- iii)Secret interactions just restricted to the sterile form such a measurement, and others are planned to do neutrino sector νsνs ! νsνs for different sterile neutrino it in future with much better chances, like for instance mass scales, [62{74]. 2 iv) Secret interactions involving active and sterile The quantity F [z0;E(1 + z0)] depends of course on neutrinos simultaneously νν ! νsνs[75, 76] the proton spectrum, which is itself the solution of a Boltzmann equation (see [38]), which takes into account In this paper we consider a scheme of SI similar to the proton energy losses due to Bethe-Heitler processes the point iv, wherein the new interaction, mediated by a and their depletion due to pγ processes. This calcula- new pseudoscalar boson, intervene both active and sterile tion has been performed for example in [38], which pro- neutrinos. While [75] explores the effects of this interac- vides the neutrino spectrum expected at Earth. It is tion on primordial nucleosynthesis, and [76] studies the important to notice that they assume a cosmic ray spec- effects of the interaction on neutrinos in the IceCube en- trum purely made of protons and emitted by sources ergy range of interest, we have analyzed the effects of whose density follows the star formation evolution. In this interaction on cosmogenic neutrino fluxes to ques- other words, the proton luminosity can be written as tion their observability. We will assume throughout that Lp(z; E) = H(z)Qp(E), where Qp(E) is the proton in- both the active and the sterile neutrinos are Majorana jection spectrum from the single sources and H(z) is the particles: this is also the reason why we have to choose Star Forming rate [77] a pseudoscalar mediator. In fact the scalar contraction 8 3:4 νν¯ s is antihermitean and therefore not admissible as a >(1 + z) z ≤ 1; < −0:3 possible interaction operator: therefore the only possi- H(z) = N1(1 + z) 1 < z ≤ 4; (3) ble contraction isνγ ¯ 5νs. To preserve parity, we take the > −3:5 :N1N4(1 + z) z > 4; mediator to be a pseudoscalar. In particular we study the distortion implied by such new coupling on the ex- 3:7 3:2 where N1 = 2 and N4 = 5 . The proton luminosity pected cosmogenic neutrino flux estimating the possibil- Lp(z; E) works as an input to the Boltzmann equation ity to measure this effect in apparatus like GRAND [34]. which provides the function F [z0;E(1 + z0)] in Eq. (2). For the sake of simplicity we will assume just one sterile Since Ref. [38] only furnishes the final neutrino spec- neutrino (hereafter denoted by νs) coupling with the ac- trum at Earth, and does not provide the proton spectrum tive sector via this new interaction. Moreover, with the Qp(E) at each redshift, in principle one would not be aim to catch the main implications of such scenario on able to reproduce the function F [z0;E(1 + z0)] in Eq. (2). cosmogenic neutrino flux we consider an active neutrino However, at sufficiently high energies the mean free path only (hereafter denoted by ν), hence focussing on a more for pγ interaction becomes so small that neutrinos can simple 1 + 1 framework. However it would be straight- be assumed to be produced exactly at the same place in forward to extend our analysis to three active neutrinos which the emission of the protons occurs. If we make even though we do not expect that the main results, here the further hypothesis that at each redshift the injection obtained, would be drastically changed in a more realistic spectrum of the protons has exactly the same form, with 3 + 1 framework. The interaction term then becomes a redshifted energy, the function F [z0;E(1 + z0)] takes on the form LSI = λ νγ5νs'; (1) F [z0;E(1 + z0)] = ρ(z0)f[E(1 + z0)] ; (4) where λ is a dimensionless free coupling. where ρ(z0) is proportional to the Star Forming Rate given in Eq. (3). II. COSMOGENIC NEUTRINO FLUX AT After this simplification, there is a unique function EARTH WITHOUT SECRET INTERACTIONS f[E(1 + z0)] which reproduces the spectrum obtained by [38]. While the effects of these simplifications may be As described above, cosmogenic neutrinos are pro- relevant at lower energies, they should be almost irrele- duced by the scattering of high energy protons from the vant in the high energy part of the cosmogenic spectrum, cosmic rays with the CMB photons. The cosmogenic where the hypothesis of a small mean free path is natu- neutrino flux φν , expected to be isotropic, can be pa- ral. In particular, the mean free path for pγ interactions rameterized in the form is typically of the order of 50 Mpc, which is small com- pared to the cosmological distances. dφ Z dz0 ν 0 0 The inversion of Eq. (2) under the ansatz of Eq. (4) is = 0 F [z ;E(1 + z )] ; (2) dEdΩ H(z ) in principle possible in an exact way through the use of a Mellin transform.