Search

In dense neutrino environments, such as provided by core-collapse supernovae or neutron-star mergers, neutrino angular distributions may be unstable to collective flavor conversions, whose outcome remains to be fully understood. These conversions are much faster than hydrodynamical scales, suggesting that self-consistent configurations may never be strongly unstable. With this motivation in mind, we study weakly unstable modes, i.e., those with small growth rates. We show that our newly developed dispersion relation (Paper~I of this series) allows for an expansion in powers of the small growth rate. For weakly unstable distributions, we show that the unstable modes must either move with subluminal phase velocity, or very close to the speed of light. The instability is fed from neutrinos moving resonantly with the waves, allowing us to derive explicit expressions for the growth rate. For axisymmetric distributions, often assumed in the literature, numerical examples show the accuracy of these expressions. We also note that for the often-studied one-dimensional systems one should not forget the axial-symmetry-breaking modes, and we provide explicit expressions for the range of wavenumbers that exhibit instabilities., Comment: Title changed to match Part I format

The cores of active galactic nuclei (AGN) are potential accelerators of 10-100 TeV cosmic rays, in turn producing high-energy neutrinos. This picture was confirmed by the compelling evidence of a TeV neutrino signal from the nearby active galaxy NGC 1068, leaving open the question of which is the site and mechanism of cosmic ray acceleration. One candidate is the magnetized turbulence surrounding the central supermassive black hole. Recent particle-in-cell simulations of magnetized turbulence indicate that stochastic cosmic ray acceleration is non-resonant, in contrast to the assumptions of previous studies. We show that this has important consequences on a self-consistent theory of neutrino production in the corona, leading to a more rapid cosmic ray acceleration than previously considered. The turbulent magnetic field fluctuations needed to explain the neutrino signal are consistent with a magnetically powered corona. We find that strong turbulence, with turbulent magnetic energy density higher than $1\%$ of the rest mass energy density, naturally explains the normalization of the IceCube neutrino flux, in addition to the neutrino spectral shape. Only a fraction of the protons in the corona, which can be directly inferred from the neutrino signal, are accelerated to high energies. Thus, in this framework, the neutrino signal from NGC 1068 provides a testbed for particle acceleration in magnetized turbulence., Comment: 10 pages, 5 figures

Many particles predicted by extensions of the Standard Model feature interactions with neutrinos, e.g., Majoron-like bosons $\phi$. If the mass of $\phi$ is larger than about 10 keV, they can be produced abundantly in the core of the next galactic core-collapse supernova through neutrino coalescence, and leave it with energies of around 100 MeV. Their subsequent decay to high-energy neutrinos and anti-neutrinos provides a distinctive signature at Earth. Ongoing and planned neutrino and dark matter experiments allow us to reconstruct the energy, flavor, and time of arrival of these high-energy neutrinos. For the first time, we show that these measurements can help pinpointing the mass of $\phi$ and its couplings to neutrinos of different flavor. Our results can be generalized in a straightforward manner to other hypothetical feebly interacting particles, like novel gauge bosons or heavy neutral leptons, that decay into neutrinos., Comment: 38 pages, 7 figures

Neutrino-neutrino refraction leads to collective flavor evolution that can include fast flavor conversion, an ingredient still missing in numerical simulations of core-collapse supernovae. We provide a theoretical framework for the linear regime of this phenomenon using the language of response theory. In analogy to electromagnetic waves, we introduce a flavor susceptibility as the linear response to an external flavor field. By requiring self-consistency, this approach leads to the usual dispersion relation for growing modes, but differs from the traditional treatment in that it predicts Landau damping of subluminal collective modes. The new dispersion relation has definite analyticity properties and can be expanded for small growth rates. This approach simplifies and intuitively explains Morinaga's proof of sufficiency for the occurrence of growing modes. We show that weakly growing modes arise as soon as an angular crossing is formed, due to their resonant interaction with individual neutrino modes. For longitudinal plasma waves, a similar resonance causes Landau damping or conversely, the two-stream instability., Comment: Title matched to published version

Since neutrinos have mass differences, they could decay into one another. But their lifetimes are likely long, even when shortened by new physics, so decay likely impacts neutrinos only during long trips. This makes high-energy astrophysical neutrinos, traveling for up to billions of light-years, sensitive probes of decay. However, their sensitivity must be tempered by reality. We derive from them thorough bounds on the neutrino lifetimes accounting for critical astrophysical unknowns and the nuances of neutrino detection. Using the diffuse neutrino flux, we disfavor lifetimes $\tau \lesssim 20$-450 s $(m/{\rm eV})$, based on present IceCube data, and forecast factor-of-10 improvements by upcoming detectors. Using, for the first time, neutrinos from the active galaxy NGC 1068, extant unknowns preclude placing lifetime bounds today, but upcoming detectors could disfavor $\tau \sim 100$-5000 s $(m/{\rm eV})$., Comment: 20 pages, 6 figures, 2 tables, plus appendices

We revisit supernova (SN) bounds on a hidden sector consisting of millicharged particles $\chi$ and a massless dark photon. Unless the self-coupling is fine-tuned to be small, rather than exiting the SN core as a gas, the particles form a relativistic fluid and subsequent dark QED fireball, streaming out against the drag due to the interaction with matter. Novel bounds due to excessive energy deposition in the mantle of low-energy SNe can be obtained. The cooling bounds from SN~1987A are unexpectedly not affected in the free-streaming regime. The inclusion of $\chi\bar{\chi}\rightarrow \rm SM$ substantially modifies the constraints in the trapping regime, which can only be treated hydrodynamically. Our results can be adapted to generic sub-GeV self-interacting dark sectors., Comment: 5 pages, 2 figures, with Supplemental Material

A dense neutrino gas exhibiting angular crossings in the electron lepton number is unstable and develops fast flavor conversions. Instead of assuming an unstable configuration from the onset, we imagine that the system is externally driven toward instability. We use the simplest model of two neutrino beams initially of different flavor that either suddenly appear or one or both slowly build up. Flavor conversions commence well before the putative unstable state is fully attained, and the final outcome depends on how the system is driven. Our results suggest that in an astrophysical setting, one should focus less on flavor instabilities in the neutrino radiation field and more on the external dynamics that leads to the formation of the unstable state., Comment: 5 pages, 1 figure

Flavor-dependent neutrino transport is described by a well-known kinetic equation for occupation-number matrices in flavor space. However, as an overlooked theoretical problem, we show that in the inhomogeneous case, neutrino-neutrino refractive energy is not conserved. We derive the missing gradient terms in the fast flavor limit (vanishing neutrino masses), and prove that the missing refractive energy is traded with the huge reservoir of neutrino kinetic energy through gradients of the weak interaction potential. Even small changes of the kinetic energy accommodate the refractive energy gained or lost. Flavor evolution alone is negligibly affected by the new terms., Comment: Small changes to the text, version published on Physical Review Letters; 5 pages, 1 figures, plus Supplemental Material

The paradigm-changing possibility of collective neutrino-antineutrino oscillations was recently advanced in analogy to collective flavor oscillations. However, the amplitude for the backward scattering process $\nu_{\mathbf{p}_1}\overline\nu_{\mathbf{p}_2}\to\nu_{\mathbf{p}_2}\overline\nu_{\mathbf{p}_1}$ is helicity-suppressed and vanishes for massless neutrinos, implying that there is no off-diagonal refractive index between $\nu$ and $\overline\nu$ of a single flavor of massless neutrinos. For a nonvanishing mass, collective helicity oscillations are possible, representing de-facto $\nu$--$\overline\nu$ oscillations in the Majorana case. However, such phenomena are suppressed by the smallness of neutrino masses as discussed in the previous literature., Comment: 3 pages, 2 figures, with appendices

The flavor composition of TeV--PeV astrophysical neutrinos, i.e., the proportion of neutrinos of different flavors in their flux, is a versatile probe of high-energy astrophysics and fundamental physics. Because flavor identification is challenging and the number of detected high-energy astrophysical neutrinos is limited, so far measurements of the flavor composition have represented an average over the range of observed neutrino energies. Yet, this washes out the potential existence of changes in the flavor composition with energy and weakens our sensitivity to the many models that posit them. For the first time, we measure the energy dependence of the flavor composition, looking for a transition from low to high energies. Our present-day measurements, based on the 7.5-year public sample of IceCube High-Energy Starting Events (HESE), find no evidence of a flavor transition. The observation of HESE and through-going muons jointly by next-generation neutrino telescopes Baikal-GVD, IceCube-Gen2, KM3NeT, P-ONE, TAMBO, and TRIDENT may identify a flavor transition around 200TeV by 2030. By 2040, we could infer the flavor composition with which neutrinos are produced with enough precision to establish the transition from neutrino production via the full pion decay chain at low energies to muon-damped pion decay at high energies., Comment: 23+7 pages, 7+6 figures

In dense astrophysical environments, notably core-collapse supernovae and neutron star mergers, neutrino-neutrino forward scattering can spawn flavor conversion on very short scales. Scattering with the background medium can impact collective flavor conversion in various ways, either damping oscillations or possibly setting off novel collisional flavor instabilities (CFIs). A key feature in this process is the slowness of collisions compared to the much faster dynamics of neutrino-neutrino refraction. Assuming spatial homogeneity, we leverage this hierarchy of scales to simplify the description accounting only for the slow dynamics driven by collisions. We illustrate our new approach both in the case of CFIs and in the case of fast instabilities damped by collisions. In both cases, our strategy provides new equations, the slow-dynamics equations, that simplify the description of flavor conversion and allow us to qualitatively understand the final state of the system after the instability, either collisional or fast, has saturated., Comment: Added footnote 2 to correct an algebraic error; version corresponding to Erratum in PRD

The recent discovery of astrophysical neutrinos from the Seyfert galaxy NGC 1068 suggests the presence of non-thermal protons within a compact "coronal" region close to the central black hole. The acceleration mechanism of these non-thermal protons remains elusive. We show that a large-scale magnetic reconnection layer, of the order of a few gravitational radii, may provide such a mechanism. In such a scenario, rough energy equipartition between magnetic fields, X-ray photons, and non-thermal protons is established in the reconnection region. Motivated by recent three-dimensional particle-in-cell simulations of relativistic reconnection, we assume that the spectrum of accelerated protons is a broken power law, with the break energy being constrained by energy conservation (i.e., the energy density of accelerated protons is at most comparable to the magnetic energy density). The proton spectrum is $dn_p/dE_p\propto E_p^{-1}$ below the break, and $dn_p/dE_p\propto E_p^{-s}$ above the break, with IceCube neutrino observations suggesting $s \simeq 3$. Protons above the break lose most of their energy within the reconnection layer via photohadronic collisions with the coronal X-rays, producing a neutrino signal in good agreement with the recent observations. Gamma-rays injected in photohadronic collisions are cascaded to lower energies, sustaining the population of electron-positron pairs that makes the corona moderately Compton thick., Comment: Updated discussion, 8 pages, 3 figures, with Appendices

Ultra-high-energy (UHE) cosmic neutrinos, with energies above 100 PeV, could be finally discovered in the near future. Measuring their flavor composition would reveal information about their production and propagation, but new techniques are needed for UHE neutrino telescopes to have the capabilities to do it. We address this by proposing a new way to measure the UHE neutrino flavor composition that does not rely on individual telescopes having flavor-identification capabilities. We manufacture flavor sensitivity from the joint detection by one telescope sensitive to all flavors -- the radio array of IceCube-Gen2 -- and one mostly sensitive to $\nu_\tau$ -- GRAND. From this limited flavor sensitivity, predominantly to $\nu_\tau$, and even under conservative choices of neutrino flux and detector size, we extract important insight. For astrophysics, we forecast meaningful constraints on the neutrino production mechanism; for fundamental physics, vastly improved constraints on Lorentz-invariance violation. These are the first measurement forecasts of the UHE $\nu_\tau$ content., Comment: Version published in PRD

We return to interpreting the historical SN~1987A neutrino data from a modern perspective. To this end, we construct a suite of spherically symmetric supernova models with the Prometheus-Vertex code, using four different equations of state and five choices of final baryonic neutron-star (NS) mass in the 1.36-1.93 M$_\odot$ range. Our models include muons and proto-neutron star (PNS) convection by a mixing-length approximation. The time-integrated signals of our 1.44 M$_\odot$ models agree reasonably well with the combined data of the four relevant experiments, IMB, Kam-II, BUST, and LSD, but the high-threshold IMB detector alone favors a NS mass of 1.7-1.8 M$_\odot$, whereas Kam-II alone prefers a mass around 1.4 M$_\odot$. The cumulative energy distributions in these two detectors are well matched by models for such NS masses, and the previous tension between predicted mean neutrino energies and the combined measurements is gone, with and without flavor swap. Generally, our predicted signals do not strongly depend on assumptions about flavor mixing, because the PNS flux spectra depend only weakly on antineutrino flavor. While our models show compatibility with the events detected during the first seconds, PNS convection and nucleon correlations in the neutrino opacities lead to short PNS cooling times of 5-9 s, in conflict with the late event bunches in Kam-II and BUST after 8-9 s, which are also difficult to explain by background. Speculative interpretations include the onset of fallback of transiently ejected material onto the NS, a late phase transition in the nuclear medium, e.g., from hadronic to quark matter, or other effects that add to the standard PNS cooling emission and either stretch the signal or provide a late source of energy. More research, including systematic 3D simulations, is needed to assess these open issues., Comment: Added new Appendix D, additional references included; 53 pages, 31 figures, 8 tables

Neutrino-neutrino scattering could have a large secret component that would turn neutrinos within a supernova (SN) core into a self-coupled fluid. Neutrino transport within the SN core, emission from its surface, expansion into space, and the flux spectrum and time structure at Earth might all be affected. We examine these questions from first principles. First, diffusive transport differs only by a modified spectral average of the interaction rate. We next study the fluid energy transfer between a hot and a cold blackbody surface in plane-parallel and spherical geometry. The key element is the decoupling process within the radiating bodies, which themselves are taken to be isothermal. For a zero-temperature cold plate, mimicking radiation into free space by the hot plate, the energy flux is 3--4\% smaller than the usual Stefan-Boltzmann Law. The fluid energy density just outside the hot plate is numerically 0.70 of the standard case, the outflow velocity is the speed of sound $v_s=c/\sqrt{3}$, conspiring to a nearly unchanged energy flux. Our results provide the crucial boundary condition for the expansion of the self-interacting fluid into space, assuming an isothermal neutrino sphere. We also derive a dynamical solution, assuming the emission suddenly begins at some instant. A neutrino front expands in space with luminal speed, whereas the outflow velocity at the radiating surface asymptotically approaches $v_s$ from above. Asymptotically, one thus recovers the steady-state emission found in the two-plate model. A sudden end to neutrino emission leads to a fireball with constant thickness equal to the duration of neutrino emission., Comment: Version accepted at PRD, with revised notation; 20 pages, 8 figures

When hypothetical neutrino secret interactions ($\nu$SI) are large, they form a fluid in a supernova (SN) core, flow out with sonic speed, and stream away as a fireball. For the first time, we tackle the complete dynamical problem and solve all steps, systematically using relativistic hydrodynamics. The impact on SN physics and the neutrino signal is remarkably small. For complete thermalization within the fireball, the observable spectrum changes in a way that is independent of the coupling strength. One potentially large effect beyond our study is quick deleptonization if $\nu$SI violate lepton number. By present evidence, however, SN physics leaves open a large region in parameter space, where laboratory searches and future high-energy neutrino telescopes will probe $\nu$SI., Comment: Version accepted in PRL, added Supplemental Material with derivation of fluid equations and comparison with previous literature; 5 pages, 3 figures, plus Supplemental Material

The origin of the TeV--PeV astrophysical neutrinos seen by the IceCube telescope is unknown. If they are made in proton-photon interactions in astrophysical sources, their spectrum may show bump-like features. We search for such features in the 7.5-years High-Energy Starting Events (HESE), and forecast the power of such searches using larger data samples expected from upcoming telescopes. Present-day data reveals no evidence of bump-like features, which allows us to constrain candidate populations of photohadronic neutrino sources. Near-future forecasts show promising potential for stringent constraints or decisive discovery of bump-like features. Our results provide new insight into the origins of high-energy astrophysical neutrinos, complementing those from point-source searches., Comment: Submitted as a proceeding for ICRC 2023. arXiv admin note: substantial text overlap with arXiv:2301.00024

Next-generation ultra-high-energy (UHE) neutrino telescopes, presently under planning, will have the potential to probe the decay of heavy dark matter (DM) into UHE neutrinos, with energies in excess of $10^7$~GeV. Yet, this potential may be deteriorated by the presence of an unknown background of UHE neutrinos, cosmogenic or from astrophysical sources, not of DM origin and seemingly large enough to obscure the DM signature. We show that leveraging the angular and energy distributions of detected events safeguards future searches for DM decay against such backgrounds. We focus on the radio-detection of UHE neutrinos in the planned IceCube-Gen2 neutrino telescope, which we model in state-of-the-art detail. We report promising prospects for the discovery potential of DM decay into UHE neutrinos, the measurement of DM mass and lifetime, and limits on the DM lifetime, despite the presence of a large background, without prior knowledge of its size and shape., Comment: Version accepted for publication in PRD; 19 pages, 13 figures, plus 4 figures in appendix

The metastable hypermassive neutron star produced in the coalescence of two neutron stars can copiously produce axions that radiatively decay into $\mathcal{O}(100)$~MeV photons. These photons can form a fireball with characteristic temperature smaller than $1\rm\, MeV$. By relying on X-ray observations of GW170817/GRB 170817A with CALET CGBM, Konus-Wind, and Insight-HXMT/HE, we present new bounds on the axion-photon coupling for axion masses in the range $1$-$400\,\rm MeV$. We exclude couplings down to $5\times 10^{-11}\,\rm GeV^{-1}$, complementing and surpassing existing constraints. Our approach can be extended to any feebly-interacting particle decaying into photons., Comment: Version accepted at Phys. Rev. Lett.; 5 pages, 3 figures + 5 pages, 3 figures

We consider a dense neutrino gas in the "fast-flavor limit" (vanishing neutrino masses). For the first time, we identify exact solutions of the nonlinear wave equation in the form of solitons. They can propagate with both sub- or superluminal speed, the latter not violating causality. The soliton with infinite speed is a homogeneous solution and coincides with the usual fast-flavor pendulum except that it swings only once instead of being periodic. The subluminal soliton in the static limit corresponds to a one-swing "spatial pendulum". A necessary condition for such solutions to exist is a ``crossed'' neutrino angle distribution. Based on the Nyquist criterion, we derive a new sufficient condition without solving the dispersion relation. The solitons are very fragile: they are as unstable as the homogeneous neutrino gas alone. Moreover, in the presence of matter, only the solution survives that is homogeneous in a frame comoving with the matter current. Generally, the matter effect cannot be eliminated by transformations in flavor space, but instead has a real physical impact., Comment: 12 pages, 8 figures, plus appendices

New feebly interacting particles would emerge from a supernova core with 100-MeV-range energies and produce $\gamma$-rays by subsequent decays. These would contribute to the diffuse cosmic $\gamma$-ray background or would have shown up in the Solar Maximum Mission (SMM) satellite from SN~1987A. However, we show for the example of axion-like particles (ALPs) that, even at distances beyond the progenitor star, the decay photons may not escape, and can instead form a fireball, a plasma shell with $T\lesssim 1$ MeV. Thus, existing arguments do not exclude ALPs with few 10 MeV masses and a two-photon coupling of a few $10^{-10}~{\rm GeV}^{-1}$. However, the energy would have showed up in sub-MeV photons, which were not seen from SN 1987A in the Pioneer Venus Orbiter (PVO), closing again this new window. A careful re-assessment is required for other particles that were constrained in similar ways., Comment: 12 pages, 7 figures. Corrected spectrum for the photons, according to Erratum submitted to Phys. Rev. D. Conclusions are unchanged

The flavor evolution of a neutrino gas can show ''slow'' or ''fast'' collective motion. In terms of the usual Bloch vectors to describe the mean-field density matrices of a homogeneous neutrino gas, the slow two-flavor equations of motion (EOMs) are $\dot{\mathbf{P}}_\omega=(\omega\mathbf{B}+\mu\mathbf{P})\times\mathbf{P}_\omega$, where $\omega=\Delta m^2/2E$, $\mu=\sqrt{2} G_{\mathrm{F}} (n_\nu+n_{\bar\nu})$, $\mathbf{B}$ is a unit vector in the mass direction in flavor space, and $\mathbf{P}=\int d\omega\,\mathbf{P}_\omega$. For an axisymmetric angle distribution, the fast EOMs are $\dot{\mathbf{D}}_v=\mu(\mathbf{D}_0-v\mathbf{D}_1)\times\mathbf{D}_v$, where $\mathbf{D}_v$ is the Bloch vector for lepton number, $v=\cos\theta$ is the velocity along the symmetry axis, $\mathbf{D}_0=\int dv\,\mathbf{D}_v$, and $\mathbf{D}_1=\int dv\,v\mathbf{D}_v$. We discuss similarities and differences between these generic cases. Both systems can have pendulum-like instabilities (soliton solutions), both have similar Gaudin invariants, and both are integrable in the classical and quantum case. Describing fast oscillations in a frame comoving with $\mathbf{D}_1$ (which itself may execute pendulum-like motions) leads to transformed EOMs that are equivalent to an abstract slow system. These conclusions carry over to three flavors., Comment: 16 pages; typo below Eq. (5) and in Eq. (47) corrected

The origin of the bulk of the high-energy astrophysical neutrinos seen by IceCube, with TeV--PeV energies, is unknown. If they are made in photohadronic, i.e., proton-photon, interactions in astrophysical sources, this may manifest as a bump-like feature in their diffuse flux, centered around a characteristic energy. We search for evidence of this feature, allowing for variety in its shape and size, in 7.5 years of High-Energy Starting Events (HESE) collected by the IceCube neutrino telescope, and make forecasts using larger data samples from upcoming neutrino telescopes. Present-day data reveals no evidence of bump-like features, which allows us to constrain candidate populations of photohadronic neutrino sources. Near-future forecasts show promising potential for stringent constraints or decisive discovery of bump-like features. Our results provide new insight into the origins of high-energy astrophysical neutrinos, complementing those from point-source searches., Comment: 29 pages, 13 figures

Starburst galaxies are well-motivated astrophysical emitters of high-energy gamma-rays. They are well-known cosmic-ray "reservoirs", thanks to their large magnetic fields which confine high-energy protons for $\sim 10^5$ years. Over such long times, cosmic-ray transport can be significantly affected by scatterings with sub-GeV dark matter. Here we point out that this scattering distorts the cosmic-ray spectrum, and the distortion can be indirectly observed by measuring the gamma-rays produced by cosmic-rays via hadronic collisions. Present gamma-ray data show no sign of such a distortion, leading to stringent bounds on the cross section between protons and dark matter. These are highly complementary with current bounds and have large room for improvement with the future gamma-ray measurements in the 0.1-10 TeV range from the Cherenkov Telescope Array, which can strengthen the limits by as much as two orders of magnitude., Comment: 17 pages, 7 figures. v2: version published in Physical Review Letters

Majoron-like bosons would emerge from a supernova (SN) core by neutrino coalescence of the form $\nu\nu\to\phi$ and $\bar\nu\bar\nu\to\phi$ with 100 MeV-range energies. Subsequent decays to (anti)neutrinos of all flavors provide a flux component with energies much larger than the usual flux from the "neutrino sphere." The absence of 100 MeV-range events in the Kamiokande-II and Irvine-Michigan-Brookhaven signal of SN 1987A implies that less than 1% of the total energy was thus emitted and provides the strongest constraint on the Majoron-neutrino coupling of $g\lesssim 10^{-9}\,{\rm MeV}/m_\phi$ for $100~{\rm eV}\lesssim m_\phi\lesssim 100~{\rm MeV}$. It is straightforward to extend our new argument to other hypothetical feebly interacting particles., Comment: 5+8 pages, 2+6 figures. Matches published version

The discovery of ultra-high-energy neutrinos, with energies above 100 PeV, may soon be within reach of upcoming neutrino telescopes. We present a robust framework to compute the statistical significance of point-source discovery via the detection of neutrino multiplets. We apply it to the radio array component of IceCube-Gen2. To identify a source with $3\sigma$ significance, IceCube-Gen2 will need to detect a triplet, at best, and an octuplet, at worst, depending on whether the source is steady-state or transient, and on its position in the sky. The discovery, or absence, of sources significantly constrains the properties of the source population., Comment: Submission to SciPost Phys. Proc., presented at the 21st International Symposium on Very High Energy Cosmic Ray Interactions (ISVHECRI 2022)

The detection of very high-energy neutrinos by IceCube experiment supports the existence of a comparable gamma-ray counterpart from the same cosmic accelerators. Under the likely assumption that the sources of these particles are of extragalactic origin, the emitted photon flux would be significantly absorbed during its propagation over cosmic distances. However, in the presence of photon mixing with ultra-light axion-like-particles (ALPs), this expectation would be strongly modified. Notably, photon-ALP conversions in the host galaxy would produce an ALP flux which propagates unimpeded in the extragalactic space. Then, the back-conversion of ALPs in the Galactic magnetic field leads to a diffuse high-energy photon flux. In this context, the recent detection of the diffuse high-energy photon flux by the Large High Altitude Air Shower Observatory (LHAASO) allows us to exclude at the $95\%$ CL an ALP-photon coupling $g_{a\gamma}\gtrsim 3.9-7.8 \times 10^{-11}~\mathrm{GeV^{-1}}$ for $m_{a}\lesssim 4\times10^{-7}~\mathrm{eV}$, depending on the assumptions on the magnetic fields and on the original gamma-ray spectrum. This new bound is complementary with other ALP constraints from very-high-energy gamma-ray experiments and sensitivities of future experiments., Comment: 17 pages, 14 figures

Upcoming neutrino telescopes may discover ultra-high-energy (UHE) cosmic neutrinos, with energies beyond 100 PeV, in the next 10-20 years. Finding their sources would identify guaranteed sites of interaction of UHE cosmic rays, whose origin is unknown. We search for sources by looking for multiplets of UHE neutrinos arriving from similar directions. Our forecasts are state-of-the-art, geared at neutrino radio-detection in IceCube-Gen2. They account for detector energy and angular response, and for critical, but uncertain backgrounds. Sources at declination of $-45^\circ$ to $0^\circ$ will be easiest to discover. Discovering even one steady-state source in 10 years would imply that the source has an UHE neutrino luminosity at least larger than about $10^{43}$ erg/s (depending on the source redshift evolution). Discovering no transient source would disfavor transient sources brighter than $10^{53}$ erg as dominant. Our results aim to inform the design of upcoming detectors., Comment: 12 pages, 3 figures, plus references and technical appendices

Primordial black holes are a possible component of dark matter, and a most promising way of investigating them is through the product of their Hawking evaporation. As a result of this process, any species lighter than the Hawking temperature is emitted, including possible new particles beyond the Standard Model. These can then be detected in lab-based experiments via their interaction with the Standard Model particles. In a previous work, we have first proposed and studied this scenario in the presence of an interaction between the light new species and nucleons. Here we extend this discussion to include the case of interaction with electrons. We show that the simultaneous presence of primordial black holes and species lighter than about $100$ MeV can be constrained by the measurements of direct detection experiments, such as XENON1T, and water Cherenkov neutrino detectors, such as Super-Kamiokande. Our results provide a complementary and alternative way of investigation with respect to cosmological and collider searches., Comment: 11 pages, 3 figures v2: version published in Physical Review D text typo corrects

Neutrinos are the Standard Model (SM) particles which we understand the least, often due to how weakly they interact with the other SM particles. Beyond this, very little is known about interactions among the neutrinos, i.e., their self-interactions. The SM predicts neutrino self-interactions at a level beyond any current experimental capabilities, leaving open the possibility for beyond-the-SM interactions across many energy scales. In this white paper, we review the current knowledge of neutrino self-interactions from a vast array of probes, from cosmology, to astrophysics, to the laboratory. We also discuss theoretical motivations for such self-interactions, including neutrino masses and possible connections to dark matter. Looking forward, we discuss the capabilities of searches in the next generation and beyond, highlighting the possibility of future discovery of this beyond-the-SM physics., Comment: Editors: Nikita Blinov, Mauricio Bustamante, Kevin J. Kelly, Yue Zhang. 29 pages, 16 figures, plus references. Contribution to Snowmass 2021

Axion-like particles (ALP) can in principle be produced in very hot and dense astrophysical environments, escape from the extreme object where such conditions are met, and then be converted in gamma--rays in the magnetic fields intervening between the event and the Earth. This process potentially offers a new window on both the physics of the axions, and the inner working of the astrophysical objects where they are produced. Interestingly, while this process has been studied for core--collapse supernovae and other extreme astrophysical events, no estimate exists for Neutron Star Mergers, objects recently identified through the detection of gravitational waves. In this work we study the production of ALPs in neutron star mergers, finding that for a large region of the ALP parameter space its magnitude at the source is such to produce a sizable gamma-ray signal at Earth. We show detection forecasts for such events placed in nearby galaxies, finding that they are potentially observable with the Fermi--LAT, thus opening a new window into both the astrophysics of these cataclysmic events, and of new particles beyond the standard model., Comment: 5 pages, 3 figures

Among the several strategies for indirect searches of dark matter, one very promising one is to look for the gamma-rays from decaying dark matter. Here we use the most up-to-date upper bounds on the gamma-ray flux from $10^5$ to $10^{11}$ GeV, obtained from CASA-MIA, KASCADE, KASCADE-Grande, Pierre Auger Observatory, and Telescope Array. We obtain global limits on dark matter lifetime in the range of masses $m_\mathrm{DM}=[10^7-10^{15}]~\mathrm{GeV}$. We provide the bounds for a set of decay channels chosen as representatives. The constraints derived here are new and cover a region of the parameter space not yet explored. We compare our results with the projected constraints from future neutrino telescopes, in order to quantify the improvement that will be obtained by the complementary high-energy neutrino searches., Comment: 8 pages, 5 figures

Astrophysical neutrino fluxes are often modeled as power laws of the energy. This is reasonable in the case of hadronic sources, but it does not capture the behavior in photohadronic sources, where the spectrum depends on the properties of the target photons on which protons collide. This limits the possibility of a unified treatment of different sources. In order to overcome this difficulty, we model the target photons by a blackbody spectrum. This model is sufficiently flexible to reproduce neutrino fluxes from known photohadronic sources; we apply it to study the sensitivity of Dense Neutrino Arrays, Neutrino Telescopes and Neutrino Radio Arrays to photohadronic sources. We also classify the flavor composition of the neutrino spectrum in terms of the parameter space. We discuss the interplay with the experiments, studying the changes in the track-to-shower ratio induced by different flavor compositions, both within and outside the region of the Glashow resonance., Comment: Presented at the 37th International Cosmic Ray Conference (ICRC 2021)

Among the information provided by high energy neutrinos, a promising possibility is to analyze the effects of a Violation of Equivalence Principle (VEP) on neutrino oscillations. We analyze the IceCube data on atmospheric neutrino fluxes under the assumption of a VEP and obtain updated constraints on the parameter space with the benchmark choice that neutrinos with different masses couple with different strengths to the gravitational field. In this case we find that the VEP parameters times the local gravitational potential at Earth can be constrained at the level of $10^{-27}$. We show that the constraints from atmospheric neutrinos strongly depend on the assumption that the neutrino eigenstates interacting diagonally with the gravitational field coincide with the mass eigenstates, which is not a priori justified: this is particularly clear in the case that the basis of diagonal gravitational interaction coincide with the flavor basis, which cannot be constrained by the observation of atmospheric neutrinos. Finally, we quantitatively study the effect of a VEP on the flavor composition of the astrophysical neutrinos, stressing again the interplay with the basis in which the VEP is diagonal: we find that for some choices of such basis the flavor ratio measured by IceCube can significantly change., Comment: Presented at the 37th International Cosmic Ray Conference (ICRC 2021)

The direct detection of sub-GeV dark matter interacting with nucleons is hampered by the low recoil energies induced by scatterings in the detectors. This experimental difficulty is avoided in the scenario of boosted dark matter where a component of dark matter particles is endowed with large kinetic energies. In this Letter, we point out that the current evaporation of primordial black holes with masses from $10^{14}$ to $10^{16}$ g is a source of boosted light dark matter with energies of tens to hundreds of MeV. Focusing on the XENON1T experiment, we show that these relativistic dark matter particles could give rise to a signal orders of magnitude larger than the present upper bounds. Therefore, we are able to significantly constrain the combined parameter space of primordial black holes and sub-GeV dark matter. In the presence of primordial black holes with a mass of $10^{15}~\mathrm{g}$ and an abundance compatible with present bounds, the limits on DM-nucleon cross-section are improved by four orders of magnitude., Comment: 8 pages, 3 figures. v2: updated to published version

The symmetry of the theory of relativity under diffeomorphisms strongly depends on the equivalence principle. Violation of Equivalence Principle (VEP) can be tested by looking for deviations from the standard framework of neutrino oscillations. In recent works, it has been shown that strong constraints on the VEP parameter space can be placed by means of the atmospheric neutrinos observed by the IceCube neutrino telescope. In this paper, we focus on the KM3NeT neutrino telescope and perform a forecast analysis to assess its capacity to probe VEP. Most importantly, we examine the crucial role played by systematic uncertainties affecting the neutrino observations. We find that KM3NeT will constrain VEP parameters times the local gravitational potential at the level of $10^{-27}$. Due to the systematic-dominated regime, independent analyses from different neutrino telescopes are fundamental for robustly testing the equivalence principle., Comment: 7 pages, 2 figures

Star-forming and starburst galaxies, which are well-known cosmic-rays reservoirs, are expected to emit gamma-rays and neutrinos predominantly via hadronic collisions. In this Letter, we analyze the 10-year Fermi-LAT spectral energy distributions of 13 nearby galaxies by means of a physical model which accounts for high-energy proton transport in starburst nuclei and includes the contribution of primary and secondary electrons. In particular, we test the hypothesis that the observed gamma-ray fluxes are mostly due to star-forming activity, in agreement with the available star formation rates coming from IR and UV observations. Through this observation-based approach, we determine the most-likely neutrino counterpart from star-forming and starburst galaxies and quantitatively assess the ability of current and upcoming neutrino telescopes to detect them as point-like sources. Remarkably, we find that the cores of the Small Magellanic Cloud and the Circinus galaxy are potentially observable by KM3NeT/ARCA with 6 years of observation. Moreover, most of the nearby galaxies are likely to be just a factor of a few below the KM3NeT and IceCube-Gen2 point-like sensitivities. After investigating the prospects for detection of gamma-rays above TeV energies from these sources, we conclude that the joint observations of high-energy neutrinos and gamma-rays with upcoming telescopes will be an objective test for our emission model and may provide compelling evidence of star-forming activity as a tracer of neutrino production., Comment: 7 pages, 2 figures, version 2, accepted and published in APJ letters

Recently Dolgov and Pozdnyakov proposed a new baryogenesis mechanism in which baryon asymmetry is produced without violating baryon number at the Lagrangian level. In this scenario, baryon asymmetry is generated by absorption of a new particle X carrying baryon number onto Primordial Black Holes (PBHs). Assuming CP-violation, the particle X is absorbed at a different rate than the antiparticle $\bar{X}$, producing an asymmetry in the baryonic number. We independently test this scenario, finding that it suffers from two fundamental issues.\\ At the phenomenological level, strong absorption by PBHs initially increases the baryon asymmetry. However, at later times such asymmetry is completely absorbed by PBHs. In order to overcome this issue, we account for PBH evaporation, which provides a natural way of halting the absorption while keeping a finite baryon asymmetry. We provide a systematic study of the parameter space, identifying the regions leading to the production of the baryon asymmetry without violating the known constraints on PBHs concentration. At the theoretical level, a model realizing the CP-violation postulated in this scenario is difficult to realize. We show, by implementing a minimal model, that the framework proposed in the original work in order to produce CP-violation, even if qualitatively correct, is quantitatively in disagreement with the observed baryon asymmetry, namely this mechanism produces only a fraction of the total baryon asymmetry., Comment: 8 pages, 3 figures

Primordial black holes (PBHs) hypothetically generated in the first instants of life of the Universe are potential dark matter (DM) candidates. Focusing on PBHs masses in the range $[5 \times10^{14} - 5 \times 10^{15}]$g, we point out that the neutrinos emitted by PBHs evaporation can interact through the coherent elastic neutrino nucleus scattering (CE$\nu$NS) producing an observable signal in multi-ton DM direct detection experiments. We show that with the high exposures envisaged for the next-generation facilities, it will be possible to set bounds on the fraction of DM composed by PBHs improving the existing neutrino limits obtained with Super-Kamiokande. We also quantify to what extent a signal originating from a small fraction of DM in the form of PBHs would modify the so-called "neutrino floor", the well-known barrier towards detection of weakly interacting massive particles (WIMPs) as the dominant DM component., Comment: 8 pages; 4 figures

The sources of IceCube neutrinos are as yet unknown. The multi-messenger observation of their emission in $\gamma$-rays can be a guide to their identification, as exemplified by the case of TXS 0506+056. We suggest a new method of searching for $\gamma$-rays with Imaging Air Cherenkov Telescopes from sources in coincidence with possible astrophysical neutrinos. We propose that searches of $\gamma$-rays are extended, from the current practice of only a few days, to up to one month from a neutrino alert. We test this strategy on simulated sources modeled after the blazar \emph{TXS 0506+056-like}, emitting neutrinos and $\gamma$-rays via photohadronic interactions: the $\gamma$-rays are subsequently reprocessed in the VHE range. Using MAGIC as a benchmark example, we show that current Cherenkov Telescopes should be able to detect$\gamma$-ray counterparts to neutrino alerts with a rate of approximately one per year. It has been proposed that the high-energy diffuse neutrino flux can be explained by $\sim$ 5\% of all blazars flaring in neutrinos once every 10 years, with a neutrino luminosity similar to that of TXS 0506+056 during the 2014-2015 neutrino flare. The implementation of our strategy could lead, over a timescale of one or few years, either to the detection of this subclass of blazars contributing to the diffuse neutrino flux, or to a constraint on this model., Comment: 10 pages, 2 figures

Many extensions of the Standard Model predict the existence of new charged or neutral gauge bosons, with a wide variety of phenomenological implications depending on the model adopted. The search for such particles is extensively carried through at the Large Hadron Collider (LHC), and it is therefore of crucial importance to have for each proposed scenario quantitative predictions that can be matched to experiments. In this work we focus on the implications of one of these models, the TopFlavor Model, proposing a charged $\text{W}^\prime$ boson that has preferential couplings to the third generation fermions. We compare such predictions to the ones from the so called Sequential Standard Model (SSM), that is used as benchmark, being one of the simplest and most commonly considered models for searches at the LHC. We identify the parameter space still open for searches at the LHC, and in particular we show that the cross section for the processes $pp \to \text{W}^\prime \to \tau \nu$ and $pp \to \text{W}^\prime \to tb$ can be up to two orders of magnitude smaller with respect to the SSM, depending on the free parameters of the model, like the particle mass and its width. This study makes the case for further searches at the LHC, and shows how a complete and systematic model independent analysis of $\text{W}^\prime$ boson phenomenology at colliders is essential to provide guidance for future searches., Comment: 11 pages, 3 figures

In the next decades, ultra-high-energy neutrinos in the EeV energy range will be potentially detected by next-generation neutrino telescopes. Although their primary goals are to observe cosmogenic neutrinos and to gain insight into extreme astrophysical environments, they can also indirectly probe the nature of dark matter. In this paper, we study the projected sensitivity of up-coming neutrino radio telescopes, such as RNO-G, GRAND and IceCube-gen2 radio array, to decaying dark matter scenarios. We investigate different dark matter decaying channels and masses, from $10^7$ to $10^{15}$ GeV. By assuming the observation of cosmogenic or newborn pulsar neutrinos, we forecast conservative constraints on the lifetime of heavy dark matter particles. We find that these limits are competitive with and highly complementary to previous multi-messenger analyses., Comment: 11 pages, 4 figures. v3: Fig.2 has been fixed, results are unchanged

Among the information provided by high energy neutrinos, a promising possibility is to analyze the effects of a Violation of Equivalence Principle (VEP) on neutrino oscillations. We analyze the recently released IceCube data on atmospheric neutrino fluxes under the assumption of a VEP and obtain updated constraints on the parameter space with the benchmark choice that neutrinos with different masses couple with different strengths to the gravitational field. In this case we find that the VEP parameters times the local gravitational potential at Earth can be constrained at the level of $10^{-27}$. We show that the constraints from atmospheric neutrinos strongly depend on the assumption that the neutrino eigenstates interacting diagonally with the gravitational field coincide with the mass eigenstates, which is not \textit{a priori} justified: this is particularly clear in the case that the basis of diagonal gravitational interaction coincide with the flavor basis, which cannot be constrained by the observation of atmospheric neutrinos. Finally, we quantitatively study the effect of a VEP on the flavor composition of the astrophysical neutrinos, stressing again the interplay with the basis in which the VEP is diagonal: we find that for some choices of such basis the flavor ratio measured by IceCube can significantly change., Comment: 9 pages, 3 figures

Starburst galaxies, which are known as "reservoirs" of high-energy cosmic-rays, can represent an important high-energy neutrino "factory" contributing to the diffuse neutrino flux observed by IceCube. In this paper, we revisit the constraints affecting the neutrino and gamma-ray hadronuclear emissions from this class of astrophysical objects. In particular, we go beyond the standard prototype-based approach leading to a simple power-law neutrino flux, and investigate a more realistic model based on a data-driven blending of spectral indexes, thereby capturing the observed changes in the properties of individual emitters. We then perform a multi-messenger analysis considering the extragalactic gamma-ray background (EGB) measured by Fermi-LAT and different IceCube data samples: the 7.5-year High-Energy Starting Events (HESE) and the 6-year high-energy cascade data. Along with starburst galaxies, we take into account the contributions from blazars and radio galaxies as well as the secondary gamma-rays from electromagnetic cascades. Remarkably, we find that, differently from the highly-constrained prototype scenario, the spectral index blending allows starburst galaxies to account for up to $40\%$ of the HESE events at $95.4\%$ CL, while satisfying the limit on the non-blazar EGB component. Moreover, values of $\mathcal{O}(100~\mathrm{PeV})$ for the maximal energy of accelerated cosmic-rays by supernovae remnants inside the starburst are disfavoured in our scenario. In broad terms, our analysis points out that a better modeling of astrophysical sources could alleviate the tension between neutrino and gamma-ray data interpretation., Comment: 20 pages, 15 figures. v2: updated to published version

The next generation of Imaging Atmospheric Cherenkov telescope, like CTA, is going to strongly improve the detection capability of high-energy cosmic rays. In our paper we discuss the possibility to use such apparatus to detect Earth-skimming tau neutrinos. Interestingly the analysis shows that order few events per year can be detected for energies above $10^8$ GeV in the optimistic case of larger neutrino fluxes produced by Flat Spectrum Radio Quasars. However, even for more conservative cosmogenic neutrino fluxes such rate will be also obtained, but for a decade of running. This estimate seems to candidate a set up like CTA for performing high energy neutrino astronomy as well., Comment: 9 pages, 6 figures

We consider the effects of active-sterile secret neutrino interactions, mediated by a new pseudoscalar particle, on high- and ultra high-energy neutrino fluxes. In particular, we focus on the case of 3 active and 1 sterile neutrino coupled by a flavor dependent interaction, extending the case of 1 active and 1 sterile neutrino we have recently examined. We find that, depending on the kind of interaction of sterile neutrino with the active sector, new regions of the parameter space for secret interactions are now allowed leading to interesting phenomenological implications on two benchmark fluxes we consider, namely an astrophysical power law flux, in the range below 100 PeV, and a cosmogenic flux, in the Ultrahigh energy range. First of all, the final active fluxes can present a measurable depletion observable in future experiments. Especially, in the case of only tau neutrino interacting, we find that the effects on the astrophysical power law flux can be so large to be already probed by the IceCube experiment. Moreover, we find intriguing features in the energy dependence of the flavor ratio., Comment: 11 pages, 8 figures

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 GRAND apparatus if the scalar 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., Comment: 10 pages, 3 figures