Neutrino decay modifies neutrino propagation in a unique way; not only is there flavor changing as there is in neutrino oscillations, there is also energy transport from initial to final neutrinos. The most sensitive direct probe of neutrino decay is currently IceCube which can measure the energy and flavor of neutrinos traveling over extragalactic distances. For the first time we calculate the flavor transition probability for the cases of visible and invisible neutrino decay, including the effects of the expansion of the universe, and consider the implications for IceCube. As an example, we demonstrate how neutrino decay addresses a tension in the IceCube data. We also provide a publicly available code to calculate the effect of visible decay.

This report summarizes the present status of neutrino non-standard interactions (NSI). After a brief overview, several aspects of NSIs are discussed, including connection to neutrino mass models, model-building and phenomenology of large NSI with both light and heavy mediators, NSI phenomenology in both short- and long-baseline neutrino oscillation experiments, neutrino cross-sections, complementarity of NSI with other low- and high-energy experiments, fits with neutrino oscillation and scattering data, DUNE sensitivity to NSI, effective field theory of NSI, as well as the relevance of NSI to dark matter and cosmology. We also discuss the open questions and interesting future directions that can be pursued by the community at large. This report is based on talks and discussions during the Neutrino Theory Network NSI workshop held at Washington University in St. Louis from May 29-31, 2019 (https://indico.cern.ch/event/812851/)

The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.

The sources of cosmic rays with energies above 55 EeV are still mysterious. A guaranteed associated flux of ultra high energy neutrinos known as the cosmogenic neutrino flux will be measured by next generation radio facilities, such as the proposed Giant Radio Array for Neutrino Detection (GRAND). By using the orthogonal information provided by the cosmogenic neutrino flux, we here determine the prospects of GRAND to constrain the source redshift evolution and the chemical composition of the cosmic ray sources. If the redshift evolution is known, independently on GRAND's energy resolution, GRAND with 200,000 antennas will constrain the proton/iron fraction to the $\sim5-10\%$ level after one year of data taking; on the other hand, if hints on the average source composition are given, GRAND will measure the redshift evolution of the sources to a $\sim 10\%$ uncertainty. However, the foreseen configuration of GRAND alone will not be able to break the degeneracy between redshift evolution of the sources and their composition. Our findings underline the discriminating potential of next generation radio array detectors and motivate further efforts in this direction.

The initial data from the Event Horizon Telescope (EHT) on M87$^*$, the supermassive black hole at the center of the M87 galaxy, provide direct observational information on its mass, spin, and accretion disk properties. A combination of the EHT data and other constraints provide evidence that M87$^*$ has a mass $\sim 6.5 \times 10^9\,M_\odot$ and dimensionless spin parameter $|a^*|\gtrsim 0.5$. These determinations disfavor ultra light bosons of mass $\mu_b\sim 10^{-21}$ eV, within the range considered for fuzzy dark matter, invoked to explain dark matter distribution on $\sim$ kpc scales. Future observations of M87$^*$ could be expected to strengthen our conclusions.

Constraints on inflationary models typically assume only the standard models of cosmology and particle physics. By extending the neutrino sector to include a new interaction with a light scalar mediator ($m_{\phi}\sim$MeV), it is possible to relax these constraints, in particular via opening up regions of the parameter space of the spectral index $n_s$. These new interactions can be probed at IceCube via interactions of astrophysical neutrinos with the Cosmic Neutrino Background for nearly all of the relevant parameter space.

The IceCube Neutrino Observatory detects high energy astrophysical neutrinos in two event topologies: tracks and cascades. Since the flavor composition of each event topology differs, tracks and cascades can be used to test the neutrino properties and the mechanisms behind the neutrino production in astrophysical sources. Assuming a conventional model for the neutrino production, the IceCube data sets related to the two channels are in $>3\sigma$ tension with each other. Invisible neutrino decay with lifetime $\tau/m=10^2$ s/eV solves this tension. Noticeably, it leads to an improvement over the standard non-decay scenario of more than $3\sigma$ while remaining consistent with all other multi-messenger observations. In addition, our invisible neutrino decay model predicts a reduction of $59\%$ in the number of observed $\nu_\tau$ events which is consistent with the current observational deficit.

The detection of the diffuse supernova neutrino background (DSNB) will preciously contribute to gauge the properties of the core-collapse supernova population. We estimate the DSNB event rate in the next-generation neutrino detectors, Hyper-Kamiokande enriched with Gadolinium, JUNO, and DUNE. The determination of the supernova unknowns through the DSNB will be heavily driven by Hyper-Kamiokande, given its higher expected event rate, and complemented by DUNE that will help in reducing the parameters uncertainties. Meanwhile, JUNO will be sensitive to the DSNB signal over the largest energy range. A joint statistical analysis of the expected rates in 20 years of data taking from the above detectors suggests that we will be sensitive to the local supernova rate at most at a 20-33% level. A non-zero fraction of supernovae forming black holes will be confirmed at a 90% CL, if the true value of that fraction is larger than 20%. On the other hand, the DSNB events show extremely poor statistical sensitivity to the nuclear equation of state and mass accretion rate of the progenitors forming black holes.

The increasing statistics of the high-energy neutrino flux observed by the IceCube Observatory points towards an excess of events above the atmospheric neutrino background in the 30--400 TeV energy range. Such an excess is compatible with the findings of the ANTARES Telescope and it would naturally imply the possibility that more than one source class contributes to the observed flux. Electromagnetically hidden sources have been invoked to interpret this excess of events at low energies. By adopting a unified model for the electromagnetically bright and choked gamma-ray bursts and taking into account particle acceleration at the internal and collimation shock radii, we discuss whether bright and choked bursts are viable candidates. Our findings suggest that, although producing a copious neutrino flux, choked and bright astrophysical jets cannot be the dominant sources of the excess of neutrino events. A fine tuning of the model parameters or distinct scenarios for choked jets should be invoked in order to explain the low-energy neutrino data of IceCube and ANTARES.

Long duration gamma-ray bursts (GRBs) have often been considered as the natural evolution of some core-collapse supernova (CCSN) progenitors. However, the fraction of CCSNe linked to astrophysical jets and their properties are still poorly constrained. While any successful astrophysical jet harbored in a CCSN should produce high energy neutrinos, photons may be able to successfully escape the stellar envelope only for a fraction of progenitors, possibly leading to the existence of high-luminosity, low-luminosity and not-electromagnetically bright ("choked") GRBs. By postulating a CCSN-GRB connection, we accurately model the jet physics within the internal-shock GRB model and assume scaling relations for the GRB parameters that depend on the Lorentz boost factor $\Gamma$. The IceCube high energy neutrino flux is then employed as an upper limit of the neutrino background from electromagnetically bright and choked GRBs to constrain the jet and the progenitor properties. The current IceCube data set is compatible with up to $1\%$ of all CCSNe harboring astrophysical jets. Interestingly, those jets are predominantly choked. Our findings suggest that neutrinos can be powerful probes of the burst physics and can provide major insights on the CCSN-GRB connection.

High energy neutrinos have been detected by IceCube, but their origin remains a mystery. Determining the sources of this flux is a crucial first step towards multi-messenger studies. In this work we systematically compare two classes of sources with the data: Galactic and extragalactic. We assume that the neutrino sources are distributed according to a class of Galactic models. We build a likelihood function on an event by event basis including energy, event topology, absorption, and direction information. We present the probability that each high energy event with deposited energy $E_{\rm dep}>60$ TeV in the HESE sample is Galactic, extragalactic, or background. For Galactic models considered the Galactic fraction of the astrophysical flux has a best fit value of $1.3\%$ and is $<9.5\%$ at 90\% CL. A zero Galactic flux is allowed at $<1\sigma$.

Determining anisotropies in the arrival directions of cosmic rays at the highest energy is an important task in astrophysics. It is common and useful to partition the sky into spherical harmonics as a measure of anisotropy. The two lowest nontrivial spherical harmonics, the dipole and the quadrupole, are of particular interest, since these distributions encapsulate a dominant single source and a plane of sources, as well as offering relatively high statistics. The best experiments for the detection of ultra high energy cosmic rays currently are all ground-based, with highly nonuniform exposures on the sky resulting from the fixed experimental locations on the earth. This nonuniform exposure increases the complexity and error in inferring anisotropies. It turns out that there is an optimal latitude for an experiment at which nonuniform exposure does not diminish the inference of the quadrupole moment. We derive the optimal latitude, and find that (presumably by a fortuitous coincidence) this optimal latitude runs through the largest cosmic ray experiment, Pierre Auger Observatory (PAO) in the Southern Hemisphere, and close to the largest cosmic ray experiment in the Northern Hemisphere, Telescope Array (TA). Consequently, assuming a quadrupole distribution, PAO and TA can reconstruct the cosmic ray quadrupole distribution to a high precision, without concern for their partial sky exposure.

The two main advantages of space-based observation of extreme energy ($\gtrsim5\times10^{19}$ eV) cosmic rays (EECRs) over ground based observatories are the increased field of view and the full-sky coverage with nearly uniform systematics across the entire sky. The former guarantees increased statistics, whereas the latter enables a clean partitioning of the sky into spherical harmonics. The discovery of anisotropies would help to identify the long sought origin of EECRs. We begin an investigation of the reach of a full-sky space-based experiment such as EUSO to detect anisotropies in the extreme-energy cosmic-ray sky compared to ground based partial-sky experiments such as the Pierre Auger Observatory and Telescope Array. The technique is explained here, and simulations for a Universe with just two nonzero multipoles, monopole plus either dipole or quadrupole, are presented. These simulations quantify the advantages of space-based, all-sky coverage.

We investigate the possibility to find a characteristic TeV scale quantum black holes decay signature in the data recorded by cosmic rays experiments. TeV black holes can be produced via the collisions of ultra high energetic protons (E > $10^18$ eV) with nucleons the from atmosphere. We focus on the case when the black holes decay into two particles moving in the forward direction in the Earth reference frame (back-to-back in the center of mass reference frame) and induce two overlapping showers. When reconstructing both the energy and the shape of the resultant air shower, there is a significant difference between showers induced only via standard model interactions and showers produced via the back-to-back decay of black holes as intermediate states.

We discuss a number of experimental constraints on Weinberg's Higgs portal model. In this framework, the standard model (SM) particle spectrum is extended to include one complex scalar field S and one Dirac fermion \psi. These new fields are singlets under the SM gauge group and are charged under a global U(1) symmetry. Breaking of this U(1) symmetry results in a massless Goldstone boson \alpha and a massive CP-even scalar r, and splits the Dirac fermion into two new mass-eigenstates \psi_\pm, corresponding to Majorana fermions. The interest on such a minimal SM extension is twofold. On the one hand, if the Goldstone bosons are in thermal equilibrium with SM particles until the era of muon annihilation their contribution to the effective number of neutrino species can explain the hints from cosmological observations of extra relativistic degrees of freedom at the epoch of last scattering. On the other hand, the lightest Majorana fermion \psi_- provides a plausible dark matter candidate. Mixing of r with the Higgs doublet \phi is characterized by the mass of hidden scalar m_h and the mixing angle \theta. We constrain this parameter space using a variety of experimental data, including heavy meson decays with missing energy, the invisible Higgs width, and direct dark matter searches. We show that different experimental results compress the allowed parameter space in complementary ways, covering a large range of \psi_- masses (5 GeV \alt m_- \alt 100 GeV). Though current results narrow the parameter space significantly (for the mass range of interest, \theta \alt 10^{-3} to 10^{-4}), there is still room for discovery (\alpha decoupling at the muon annihilation era requires \theta \agt 10^{-5} to 10^{-4}). In the near future, measurements from ATLAS, CMS, LHCb, NA62, XENON1T, LUX, and CDMSlite will probe nearly the full parameter space.

The two main advantages of space-based observation of extreme-energy ($\gtrsim 10^{19}$~eV) cosmic-rays (EECRs) over ground-based observatories are the increased field of view, and the all-sky coverage with nearly uniform systematics of an orbiting observatory. The former guarantees increased statistics, whereas the latter enables a partitioning of the sky into spherical harmonics. We have begun an investigation, using the spherical harmonic technique, of the reach of \J\ into potential anisotropies in the extreme-energy cosmic-ray sky-map. The technique is explained here, and simulations are presented. The discovery of anisotropies would help to identify the long-sought origin of EECRs.