The ANITA experiment has recently observed two anomalous events emerging from well below the horizon. Even though they are consistent with tau cascades, a high energy Standard Model or Beyond the Standard Model explanation is challenging and in tension with other experiments. We study under which conditions the reflection of generic radio pulses can reproduce these signals. We propose that these pulses can be resonantly produced in the ionosphere via axion-photon conversion. This naturally explains the direction and polarization of the events and avoids other experimental bounds.

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 astrophysical neutrinos discovered by IceCube have the highest detected neutrino energies --- from TeV to PeV --- and likely travel the longest distances --- up to a few Gpc, the size of the observable Universe. These features make them naturally attractive probes of fundamental particle-physics properties, possibly tiny in size, at energy scales unreachable by any other means. The decades before the IceCube discovery saw many proposals of particle-physics studies in this direction. Today, those proposals have become a reality, in spite of astrophysical unknowns. We will showcase examples of doing fundamental neutrino physics at these scales, including some of the most stringent tests of physics beyond the Standard Model. In the future, larger neutrino energies --- up to tens of EeV --- could be observed with larger detectors and further our reach.

We quantify our present knowledge of the size and flavor structure of non-standard neutrino interactions which affect the matter background in the evolution of solar, atmospheric, reactor and long-baseline accelerator neutrinos as determined by a global analysis of oscillation data - both alone and in combination with the results on coherent neutrino-nucleus scattering from the COHERENT experiment. We consider general neutral current neutrino interactions with quarks whose lepton-flavor structure is independent of the quark type. We study the dependence of the allowed ranges of non-standard interaction coefficients, the status of the LMA-D solution, and the determination of the oscillation parameters on the relative strength of the non-standard couplings to up and down quarks. Generically we find that the conclusions are robust for a broad spectrum of up-to-down strengths, and we identify and quantify the exceptional cases related to couplings whose effect in neutrino propagation in the Earth or in the Sun is severely suppressed. As a result of the study we provide explicit constraints on the effective couplings which parametrize the non-standard Earth matter potential relevant for long-baseline experiments.

We study the impact of non-standard momentum distributions of cosmic neutrinos on the anisotropy spectrum of the cosmic microwave background and the matter power spectrum of the large scale structure. We show that the neutrino distribution has almost no unique observable imprint, as it is almost entirely degenerate with the effective number of neutrino flavours, $N_{\mathrm{eff}}$, and the neutrino mass, $m_{\nu}$. Performing a Markov chain Monte Carlo analysis with current cosmological data, we demonstrate that the neutrino mass bound heavily depends on the assumed momentum distribution of relic neutrinos. The message of this work is simple and has to our knowledge not been pointed out clearly before: Cosmology allows that neutrinos have larger masses if their average momentum is larger than that of a perfectly thermal distribution. Here we provide an example in which the mass limits are relaxed by a factor of two.

In this work we revisit the question of whether Cosmology can be made compatible with scenarios with light sterile neutrinos, as invoked to explain the SBL anomalies, in the presence of self-interaction among sterile neutrinos mediated by massive gauge bosons. We examine this proposal by deriving the cosmological predictions of the model in a wide range of the model parameters including the effective interaction strength $G_X$, sterile neutrino mass $m_{\rm st}$ and active-sterile mixings. With those we perform a statistical analysis of the cosmological data from BBN, CMB, and BAO data to infer the posterior probabilities of the sterile self-interaction model parameters. BBN mostly provides information about the effective interaction strength and we find that $\log_{10}(G_X/G_F)\geq 3.6$ can describe the primordial abundances at 95% CL. Our analysis of CMB and BAO data show that when allowing a wide prior for the sterile neutrino mass its posterior is bounded to $m_{\rm st}\leq 0.91~{\rm eV}$ (95% CL) considering CMB data only and $m_{\rm st}\leq 0.37~{\rm eV}$ (95% CL) when adding the BAO information. So the mass bounds are slighly relaxed compared with that of a non-interacting sterile neutrino model but a sterile neutrino mass of 1 eV is still excluded at more than $2\sigma$ CL. Conversely if fixing the sterile neutrino mass and mixing to the values prefered by short baseline data we find that while that CMB data alone favors the self-interacting scenario, including the BAO information severly degrades the agreement with the model. Altogether we conclude then that adding the self-interaction can alleviate the tension between eV sterile neutrinos and CMB data, but when including also the BAO results the self-interacting sterile neutrino model cannot lead to a satisfactory description of the data.

Cosmic-ray interactions with the nuclei of the Earth's atmosphere produce a flux of neutrinos in all directions with energies extending above the TeV scale. However, the Earth is not a fully transparent medium for neutrinos with energies above a few TeV. At these energies, the charged-current neutrino-nucleon cross section is large enough so that the neutrino mean-free path in a medium with the Earth's density is comparable to the Earth's diameter. Therefore, when neutrinos of these energies cross the Earth, there is a non-negligible probability for them to be absorbed. Since this effect depends on the distance traveled by neutrinos and on their energy, studying the zenith and energy distributions of TeV atmospheric neutrinos passing through the Earth offers an opportunity to infer the Earth's density profile. Here we perform an Earth tomography with neutrinos using actual data, the publicly available one-year through-going muon sample of the atmospheric neutrino data of the IceCube neutrino telescope. We are able to determine the mass of the Earth, its moment of inertia, the mass of the Earth's core and to establish the core is denser than the mantle, using weak interactions only, in a way completely independent from gravitational measurements. Our results confirm that this can be achieved with current neutrino detectors. This method to study the Earth's internal structure, complementary to the traditional one from geophysics based on seismological data, is starting to provide useful information and it could become competitive as soon as more statistics is available thanks to the current and larger future neutrino detectors.

An ultra-light axion field with mass $\sim 10^{-22}\ {\rm eV}$, also known as wave or fuzzy dark matter, has been proposed as a component of the dark matter in the Universe. We study the evolution of the axion dark matter distribution in the central region of a halo, assuming the mass is dominated by this axion field, and that gravity is the only important interaction. We calculate the excited axion states in the spherical gravitational potential generated by the self-gravitating ground-state, also known as soliton. These excited states are similar to the states of the hydrogen atom with quantum numbers $(n,l,m)$, here designating oscillation modes of a classical wave. At fixed $n$, the modes with highest $l$ have the lowest energy because of the extended mass distribution generating the potential. We use an approximate analytical treatment to derive the distribution of mass in these states when a steady-state is reached by dynamical relaxation, and find that a corona with a mass density profile $\rho \propto r^{-5/3}$ should be set up around the central soliton, analogous to the Bahcall-Wolf cusp predicted for the stellar distribution around a central black hole. The central soliton accretes dark matter from the corona as dynamical relaxation proceeds and negative orbital energy flows out. This density profile should remain valid out to the radius where the mass in the corona is comparable to the mass of the central soliton; further than that, the gravitational potential depends on the initial distribution of dark matter and the relaxation time increases rapidly with radius.

Although there is overwhelming evidence of dark matter from its gravitational interaction, we still do not know its precise gravitational interaction strength or whether it obeys the equivalence principle. Using the latest available cosmological data and working within the framework of $\Lambda\mbox{CDM}$, we first update the measurement of the Newton's constant for all matter: $G_N=7.26^{+0.27}_{-0.27}\times 10^{-11}\,\mbox{m}^{3}\mbox{kg}^{-1}\mbox{s}^{-2}$, which differs by $2.2 \sigma$ from the standard laboratory-based value. In general relativity, dark matter equivalence principle breaking can be mimicked by a long-range dark matter force mediated by an ultra light scalar field. Using the Planck three year data, we find that the dark matter "fifth-force" strength is constrained to be weaker than $10^{-4}$ of the gravitational force. We also introduce a phenomenological, post-Newtonian two-fluid description to explicitly break the equivalence principle by introducing a difference between dark matter inertial and gravitational masses. Depending on the decoupling time of the dark matter and ordinary matter fluids, the ratio of the dark matter gravitational mass to inertial mass is constrained to be unity at the $10^{-6}$ level.

In this paper we revisit the question of the information which cosmology provides on the scenarios with sterile neutrinos invoked to describe the SBL anomalies using Bayesian statistical tests. We perform an analysis of the cosmological data in $\Lambda$CDM$+r+\nu_s$ cosmologies for different cosmological data combinations, and obtain the marginalized cosmological likelihood in terms of the two relevant parameters, the sterile neutrino mass $m_s$ and its contribution to the energy density of the early Universe $N_{\rm eff}$. We then present an analysis to quantify at which level a model with one sterile neutrino is (dis)favoured with respect to a model with only three active neutrinos, using results from both short-baseline experiments and cosmology. We study the dependence of the results on the cosmological data considered, in particular on the inclusion of the recent BICEP2 results and the SZ cluster data from the Planck mission. We find that only when the cluster data is included the model with one extra sterile neutrino can become more favoured that the model with only the three active ones provided the sterile neutrino contribution to radiation density is suppressed with respect to the fully thermalized scenario. We have also quantified the level of (in)compatibility between the sterile neutrino masses implied by the cosmological and SBL results.

We investigate whether a subset of high-energy events observed by IceCube may be due to neutrinos from Sagittarius A*. We check both spatial and temporal coincidences of IceCube events with other transient activities of Sagittarius A*. Among the seven IceCube shower events nearest to the galactic center, we have found that event 25 has a time very close to (around three hours after) the brightest X-ray flare of Sagittarius A* observed by the Chandra X-ray Observatory with a p-value of 0.9%. Furthermore, two of the seven events occurred within one day of each other (there is a 1.6% probability that this would occur for a random distribution in time). Thus, the determination that some IceCube events occur at similar times as X-ray flares and others occur in a burst could be the smoking gun that Sagittarius A* is a point source of very high energy neutrinos. We point out that if IceCube galactic center neutrino events originate from charged pion decays, then TeV gamma rays should come from neutral pion decays at a similar rate. We show that the CTA, HAWC, H.E.S.S. and VERITAS experiments should be sensitive enough to test this.

We perform a geometric analysis for the sky map of the IceCube TeV-PeV neutrino excess and test its compatibility with the sky map of decaying dark matter signals in our galaxy. Using both Kolmogorov-Smirnov and the likelihood-ratio tests, we have found that the observed event sky map prefers to have a combination of the galactic dark matter and a homogeneous background contributions, compared to a purely galactic dark matter origin. For the assumption that the galactic dark matter is responsible for all neutrino excess, the current data can also exclude a wide range of dark matter profiles except flatter profiles such as the isothermal one. We also consider several representative decaying dark matter spectra, which can provide a good fit to the observed spectrum at IceCube with a dark matter lifetime of around 12 orders of magnitude longer than the age of the universe.

Recent cosmological measurements favour additional relativistic energy density beyond the one provided by the three active neutrinos and photons of the Standard Model (SM). This is often referred to as "dark radiation", suggesting the need of new light states in the theory beyond those of the SM. In this paper, we study and numerically explore the alternative possibility that this increase comes from the decay of some new form of heavy matter into the SM neutrinos. We study the constraints on the decaying matter density and its lifetime, using data from the Wilkinson Microwave Anisotropy Probe, the South Pole Telescope, measurements of the Hubble constant at present time, the results from high-redshift Type-I supernovae and the information on the Baryon Acoustic Oscillation scale. We, moreover, include in our analysis the information on the presence of additional contributions to the expansion rate of the Universe at the time of Big Bang Nucleosynthesis. We compare the results obtained in this decaying matter scenario with those obtained with the standard analysis in terms of a constant $N_{\rm eff}$.

We perform a global analysis of cosmological observables in generalized cosmologies which depart from $\Lambda$CDM models by allowing non-vanishing curvature $\Omega_k\neq 0$, dark energy with equation of state with $\omega\neq -1$, the presence of additional relativistic degrees of freedom $\Delta N_{\rm rel}$, and neutrino masses $\Omega_\nu\neq 0$. By combining the data from cosmic microwave background (CMB) experiments (in particular the latest results from WMAP-7), the present day Hubble constant (H0) measurement, the high-redshift Type-I supernovae (SN) results and the information from large scale structure (LSS) surveys, we determine the parameters in the 10-dimensional parameter space for such models. We present the results from the analysis when the full shape information from the LSS matter power spectrum (LSSPS) is included versus when only the corresponding distance measurement from the baryon acoustic oscillations (BAO) is accounted for. We compare the bounds on the neutrino mass scale in these generalized scenarios with those obtained for the 6+1 parameter analysis in $\Lambda{\rm CDM}+m_\nu$ models and we also study the dependence of those on the set of observables included in the analysis. Finally we combine these results with the information on neutrino mass differences and mixing from the global analysis of neutrino oscillation experiments and derive the presently allowed ranges for the two laboratory probes of the absolute scale of neutrino mass: the effective electron neutrino mass in single beta decay and the effective Majorana neutrino mass in neutrinoless $\beta\beta$ decay.

We determine the solar neutrino fluxes from a global analysis of the solar and terrestrial neutrino data in the framework of three-neutrino mixing. Using a Bayesian approach we reconstruct the posterior probability distribution function for the eight normalization parameters of the solar neutrino fluxes plus the relevant masses and mixing, with and without imposing the luminosity constraint. This is done by means of a Markov Chain Monte Carlo employing the Metropolis-Hastings algorithm. We also describe how these results can be applied to test the predictions of the Standard Solar Models. Our results show that, at present, both models with low and high metallicity can describe the data with good statistical agreement.

We discuss the prospects of detecting the sources of ultra-high energy (UHE) cosmic ray (CR) nuclei via their emission of cosmogenic gamma-rays in the GeV to TeV energy range. These gamma-rays result from electromagnetic cascades initiated by high energy photons, electrons and positrons that are emitted by CRs during their propagation in the cosmic radiation background and are independent of the simultaneous emission of gamma-rays in the vicinity of the source. The corresponding production power by UHE CR nuclei (with mass number A and charge Z) is dominated by pion photo-production (~ A) and Bethe-Heitler pair production (~ Z^2). We show that the cosmogenic gamma-ray signal from a single steady UHE CR source is typically more robust with respect to variations of the source composition and injection spectrum than the accompanying signal of cosmogenic neutrinos. We study the diffuse emission from the sum of extragalactic CR sources as well as the point source emission of the closest sources.

We quantify our current knowledge of the size and flavor structure of the matter effects in the evolution of atmospheric and long-baseline neutrinos based solely on the analysis of the corresponding neutrino data. To this aim we generalize the matter potential of the Standard Model by rescaling its strength, rotating it away from the e-e sector, and rephasing it with respect to the vacuum term. This phenomenological parametrization can be easily translated in terms of non-standard neutrino interactions in matter. We show that in the most general case, the strength of the potential cannot be determined solely by atmospheric and long-baseline data. However its flavor composition is very much constrained and the present determination of the neutrino masses and mixing is robust under its presence. We also present an update of the constraints arising from this analysis in the particular case in which no potential is present in the e-mu and e-tau sectors. Finally we quantify to what degree in this scenario it is possible to alleviate the tension between the oscillation results for neutrinos and antineutrinos in the MINOS experiment and show the relevance of the high energy part of the spectrum measured at MINOS.