The most efficient axion production mechanism in a supernova (SN) core is the nucleon-nucleon bremsstrahlung. This process has been often modeled at the level of the vacuum one-pion exchange (OPE) approximation. Starting from this naive recipe, we revise the calculation including systematically different effects, namely a non-vanishing mass for the exchanged pion, the contribution from the two-pions exchange, effective in-medium nucleon masses and multiple nucleon scatterings. Moreover, we allow for an arbitrary degree of nucleon degeneracy. A self consistent treatment of the axion emission rate including all these effects is currently missing. The aim of this work is to provide such an analysis. Furthermore, we demonstrate that the OPE potential with all the previous corrections gives rise to similar results as the on-shell T-matrix, and is therefore well justified for our and similar studies. We find that the axion emissivity is reduced by over an order of magnitude with respect to the basic OPE calculation, after all these effects are accounted for. The implications for the axion mass bound and the impact for the next generation experimental axion searches is also discussed.

Dark matter axion condensates may experience stimulated decays into photon pairs. This effect has been often interpreted as a parametric resonance of photons from the axion-photon coupling, leading to an exponential growth of the photon occupation number in a narrow instability band. Most of the previous literature does not consider the possible evolution of the axion field due to the photon growth. We revisit this effect presenting a mean field solution of the axion-photon kinetic equations, in terms of number of photons and pair correlations. We study the limit of no axion depletion, recovering the known instability. Moreover, we extend the results including a possible depletion of the axion field. In this case we find that the axion condensate exhibits the behaviour of an inverted pendulum. We discuss the relevance of these effects for two different cases: an homogeneous axion field at recombination and a localized axion clump and discuss constraints that could result from the induced photon background.

We update the globular cluster bound on massive ($m_a$ up to a few 100 keV) axion-like particles (ALP) interacting with photons. The production of such particles in the stellar core is dominated by the Primakoff $\gamma + Ze\to Ze +a$ and by the photon coalescence process $\gamma+\gamma\to a$. The latter, which is predominant at high masses, was not included in previous estimations. Consequently, our result strengthens the previous constraints considerably, especially for $m_a \gtrsim 50$ keV. The combined constraints from Globular Cluster stars, SN 1987A, and beam-dump experiments leave a small triangularly shaped region open in the parameter space around $m_a \sim 0.5-1\,$ MeV and $g_{a\gamma} \sim 10^{-5}$ GeV$^{-1}$. This is informally known as the ALP "cosmological triangle" since it can be excluded only using standard cosmological arguments. As we shall discuss, however, there are viable cosmological models that are compatible with axion-like particles with parameters in such region. We also discuss possibilities to explore the cosmological triangle experimentally in upcoming accelerator experiments.

Neutrinos from a supernova (SN) might undergo fast flavor conversions near the collapsed stellar core. We perform a detailed study of this intriguing possibility, analyzing time-dependent state-of-the-art 3D SN models of 9 and 20 Msun. Both models were computed with multi-D three-flavor neutrino transport based on a two-moment solver, and both exhibit the presence of the lepton-number emission self-sustained asymmetry (LESA). The transport solution does not provide the angular distributions of the neutrino fluxes, which are crucial to track the fast flavor instability. To overcome this limitation, we use a recently proposed approach based on the angular moments of the energy-integrated electron lepton-number distribution. With this method we find the possibility of fast neutrino flavor instability at radii <~20 km, which is well interior to the neutrinosphere. Our results confirm recent observations in a 2D SN model and in 2D/3D models with fixed matter background, which were computed with Boltzmann neutrino transport. However, the flavor unstable locations are not isolated points as discussed previously, but thin skins surrounding volumes where electron antineutrinos are more abundant than electron neutrinos. These volumes grow with time and appear first in the convective layer of the proto-neutron star (PNS), where a decreasing electron fraction (Ye) and high temperatures favor the occurrence of regions with negative neutrino chemical potential. Since Ye remains higher in the LESA dipole direction, where convective lepton-number transport out from the nonconvective PNS core slows down the deleptonization, flavor unstable conditions become more widespread in the opposite hemisphere. This interesting phenomenon deserves further investigation, since its impact on SN modeling and possible consequences for SN dynamics and neutrino observations are presently unclear. (abridged)

Heavy sterile neutrinos with masses ${\mathcal O}(100)$ MeV mixing with active neutrinos can be produced in the core of a collapsing supernova (SN). In order to avoid an excessive energy loss, shortening the observed duration of the SN 1987A neutrino burst, we show that the active-sterile neutrino mixing angle should satisfy $\sin^2 \theta \lesssim 5 \times 10^{-7}$. For a mixing with tau flavour, this bound is much stronger than the ones from laboratory searches. Moreover, we show that in the viable parameter space the decay of such "heavy" sterile neutrinos in the SN envelope would lead to a very energetic flux of daughter active neutrinos; if not too far below current limits, this would be detectable in large underground neutrino observatories, like Super-Kamiokande, as a (slightly time-delayed) high-energy bump in the spectrum of a forthcoming Galactic SN event.

The last decade witnessed an increasing interest in axions and axion-like particles with many theoretical works published and many new experimental proposals that started a real race towards their discovery. This paper is the Conceptual Design Report of the KLASH (KLoe magnet for Axion SearcH) experiment at the Laboratori Nazionali di Frascati (LNF). The idea of this experiment has been stimulated by the availability of the large volume superconducting magnet, with a moderate magnetic field of 0.6 T, used in the KLOE detector at the DAFNE collider. The main conclusion we draw from this report is the possibility to build and put in operation at LNF in 2-3 years a large haloscope with the sensitivity to KSVZ axions in the low mass range between 0.2 and 1 $\mu$eV, complementary to that of other experiments. Timeline and cost are competitive with respect to other proposals in the same mass region thanks to the availability of most of the infrastructure, in particular the superconducting magnet and the cryogenics plant.

We revise the theoretical initial mass-final luminosity relation for progenitors of type IIP and IIL supernovae. The effects of the major uncertainties, as those due to the treatment of convection, semiconvection, rotation, mass loss, nuclear reaction rates and neutrinos production rates are discussed in some details. The effects of mass transfer between components of close-binary systems are also considered. By comparing the theoretical predictions to a sample of type II supernovae for which the initial mass of the progenitors and the pre-explosive luminosity are available, we conclude that stellar rotation may explain a few progenitors which appear brighter than expected in case of non-rotating models. In the most extreme case, SN2012ec, an initial rotational velocity up to 300 km s$^{-1}$ is required. Alternatively, these objects could be mass-loosing components of close binaries. However, most of the observed progenitors appear fainter than expected. This occurrence seems to indicate that the Compton and pair neutrino energy-loss rates, as predicted by the standard electro-weak theory, are not efficient enough and that an additional negative contribution to the stellar energy balance is required. We show that axions coupled with parameters accessible to currently planned experiments, such as IAXO and, possibly, Baby-IAXO and ALPS II, may account for the missing contribution to the stellar energy-loss.

Fast flavor conversions of supernova neutrinos, possible near the neutrinosphere, depends on an interesting interplay of collisions and neutrino oscillations. Contrary to naive expectations, the rate of self-induced neutrino oscillations, due to neutrino-neutrino forward scattering, comfortably exceeds the rate of collisions even deep inside the supernova core. Consistently accounting for collisions and oscillations, we present the first calculations to show that collisions can create the conditions for fast flavor conversions of neutrinos, but oscillations can continue without significant damping thereafter. This may have interesting consequences for supernova explosions and the nature of its associated neutrino emission.

Neutrinos emitted from a supernova may undergo flavor conversions almost immediately above the core, with possible consequences for supernova dynamics and nucleosynthesis. However, the precise conditions for such fast conversions can be difficult to compute and require knowledge of the full angular distribution of the flavor-dependent neutrino fluxes, that is not available in typical supernova simulations. In this paper, we show that the overall flavor evolution is qualitatively similar to the growth of a so-called `zero mode', determined by the background matter and neutrino densities, which can be reliably predicted using only the second angular moments of the electron lepton number distribution, i.e., the difference in the angular distributions of $\nu_e$ and $\bar{\nu}_e$ fluxes. We propose that this zero mode, which neither requires computing the full Green's function nor a detailed knowledge of the angular distributions, may be useful for a preliminary diagnosis of possible fast flavor conversions in supernova simulations with modestly resolved angular distributions

Self-induced flavor conversions near the supernova (SN) core can make the fluxes for different neutrino species become almost equal, potentially altering the dynamics of the SN explosion and washing out all further neutrino oscillation effects. We present a new model-independent analysis strategy for the next galactic SN signal that will distinguish this flavor equalization scenario from a matter effects only scenario during the SN accretion phase. Our method does not rely on fitting or modelling the energy-dependent fluences of the different species to a known function, but rather uses a model-independent comparison of charged-current and neutral-current events at large next-generation underground detectors. Specifically, we advocate that the events due to elastic scattering on protons in a scintillator detector, which is insensitive to oscillation effects and can be used as a model-independent normalization, should be compared with the events due to inverse beta decay of $\bar\nu_e$ in a water Cherenkov detector and/or the events due to charged-current interactions of $\nu_e$ in an Argon detector. The ratio of events in these different detection channels allow one to distinguish a complete flavor equalization from a pure matter effect, for either of the neutrino mass orderings, as long as the spectral differences among the different species are not too small.

A Galactic supernova (SN) axion signal would be detected in a future neutrino Mton-class water Cherenkov detector, such as the proposed Hyper-Kamiokande in Japan. The main detection channel for axions is absorption on the oxygen nuclei in the water. The subsequent oxygen de-excitation leads to a potentially detectable gamma signal. In this contribution we present a calculation of the SN axion signal and discuss its detectability in Hyper-Kamiokande.

We present a preliminary study of the Globular Cluster RGB devoted to improve the available constraint for the axion-electron coupling. By means of multi-band IR photometry of the cluster M3 we obtain g_{ae}/10^{-13} < 2.57 (95\% C.L.).

In this study we include axions in stellar evolution models adopting the current stringest constraints for their coupling to photons and electrons. We obtain that the minimum stellar mass of Core Collapse Supernova (CCSN) progenitors is shifted up by nearly 2 Mo. This result seems to be in tension with the observationaly derived minimum mass of CCSN progenitors.

Large-scale extragalactic magnetic fields may induce conversions between very-high-energy photons and axionlike particles (ALPs), thereby shielding the photons from absorption on the extragalactic background light. However, in simplified "cell" models, used so far to represent extragalactic magnetic fields, this mechanism would be strongly suppressed by current astrophysical bounds. Here we consider a recent model of extragalactic magnetic fields obtained from large-scale cosmological simulations. Such simulated magnetic fields would have large enhancement in the filaments of matter. As a result, photon-ALP conversions would produce a significant spectral hardening for cosmic TeV photons. This effect would be probed with the upcoming Cherenkov Telescope Array detector. This possible detection would give a unique chance to perform a tomography of the magnetized cosmic web with ALPs.

Secret contact interactions among eV sterile neutrinos, mediated by a massive gauge boson $X$ (with $M_X \ll M_W$), and characterized by a gauge coupling $g_X$, have been proposed as a mean to reconcile cosmological observations and short-baseline laboratory anomalies. We constrain this scenario using the latest Planck data on Cosmic Microwave Background anisotropies, and measurements of baryon acoustic oscillations (BAO). We consistently include the effect of secret interactions on cosmological perturbations, namely the increased density and pressure fluctuations in the neutrino fluid, and still find a severe tension between the secret interaction framework and cosmology. In fact, taking into account neutrino scattering via secret interactions, we derive our own mass bound on sterile neutrinos and find (at 95% CL) $m_s < 0.82$ eV or $m_s < 0.29$ eV from Planck alone or in combination with BAO, respectively. These limits confirm the discrepancy with the laboratory anomalies. Moreover, we constrain, in the limit of contact interaction, the effective strength $G_X$ to be $ < 2.8 (2.0) \times 10^{10}\,G_F$ from Planck (Planck+BAO). This result, together with the mass bound, strongly disfavours the region with $M_X \sim 0.1$ MeV and relatively large coupling $g_X\sim 10^{-1}$, previously indicated as a possible solution to the small scale dark matter problem.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

In a Galactic core-collapse supernova (SN), axionlike particles (ALPs) could be emitted via the Primakoff process and eventually convert into $\gamma$ rays in the magnetic field of the Milky Way. From a data-driven sensitivity estimate, we find that, for a SN exploding in our Galaxy, the Fermi Large Area Telescope (LAT) would be able to explore the photon-ALP coupling down to $g_{a\gamma} \simeq 2 \times 10^{-13}\,$GeV$^{-1}$ for an ALP mass $m_a \lesssim 10^{-9}\,$eV. These values are out of reach of next generation laboratory experiments. In this event, the Fermi LAT would probe large regions of the ALP parameter space invoked to explain the anomalous transparency of the Universe to $\gamma$ rays, stellar cooling anomalies, and cold dark matter. If no $\gamma$-ray emission were to be detected, Fermi-LAT observations would improve current bounds derived from SN1987A by more than one order of magnitude.

We study the impact of axion emission in simulations of massive star explosions, as an additional source of energy loss complementary to the standard neutrino emission. The inclusion of this channel shortens the cooling time of the nascent protoneutron star and hence the duration of the neutrino signal. We treat the axion-matter coupling strength as a free parameter to study its impact on the protoneutron star evolution as well as on the neutrino signal. We furthermore analyze the observability of the enhanced cooling in current and next-generation underground neutrino detectors, showing that values of the axion mass $m_a \gtrsim 8 \times 10^{-3}$ eV can be probed. Therefore a galactic supernova neutrino observation would provide a valuable possibility to probe axion masses in a range within reach of the planned helioscope experiment, the International Axion Observatory (IAXO).

It has been recently pointed out that neutrino fluxes from a supernova can show substantial flavor conversions almost immediately above the core. Using linear stability analyses and numerical solutions of the fully nonlinear equations of motion, we perform a detailed study of these fast conversions, focussing on the region just above the supernova core. We carefully specify the instabilities for evolution in space or time, andfind that neutrinos travelling towards the core make fast conversions more generic, i.e., possible for a wider range of flux ratios and angular asymmetries that produce a crossing between the zenith-angle spectra of $\nu_e$ and ${\bar\nu_e}$. Using fluxes and angular distributions predicted by supernova simulations, we find that fast conversions can occur within tens of nanoseconds, only a few meters away from the putative neutrinospheres. If these fast flavor conversions indeed take place, they would have important implications for the supernova explosion mechanism and nucleosynthesis.

Recently it has been pointed out that a cosmic background of relativistic axion-like particles (ALPs) would be produced by the primordial decays of heavy fields in the post-inflation epoch, contributing to the extra-radiation content in the Universe today. Primordial magnetic fields would trigger conversions of these ALPs into sub-MeV photons during the dark ages. This photon flux would produce an early reionization of the Universe, leaving a significant imprint on the total optical depth to recombination $\tau$. Using the current measurement of $\tau$ and the limit on the extra-radiation content $\Delta N_{\rm eff} $ by the Planck experiment we put a strong bound on the ALP-photon conversions. Namely we obtain upper limits on the product of the photon-ALP coupling constant $g_{a\gamma}$ times the magnetic field strength $B$ down to $g_{a\gamma} B \gtrsim 6 \times 10^{-18} \textrm{GeV}^{-1} \textrm{nG} $ for ultralight ALPs.

It has been recently shown that the flavor composition of a self-interacting neutrino gas can spontaneously acquire a time-dependent pulsating component during its flavor evolution. In this work, we perform a more detailed study of this effect in a model where neutrinos are assumed to be emitted in a two-dimensional plane from an infinite line that acts as a neutrino antenna. We consider several examples with varying matter and neutrino densities and find that temporal instabilities with various frequencies are excited in a cascade. We compare the numerical calculations of the flavor evolution with the predictions of linearized stability analysis of the equations of motion. The results obtained with these two approaches are in good agreement in the linear regime, while a dramatic speed-up of the flavor conversions occurs in the non-linear regime due to the interactions among the different pulsating modes. We show that large flavor conversions can take place if some of the temporal modes are unstable for long enough, and that this can happen even if the matter and neutrino densities are changing, as long as they vary slowly.

In dense neutrino backgrounds present in supernovae and in the early Universe neutrino oscillations may exhibit complex collective phenomena, such as synchronized oscillations, bipolar oscillations and spectral splits and swaps. We consider in detail possible decoherence effects on the simplest of these phenomena -- synchronized neutrino oscillations that can occur in a uniform and isotropic neutrino gas. We develop an exact formalism of spectral moments of the flavour spin vectors describing such a system and then apply it to find analytical approaches that allow one to study decoherence effects on its late-time evolution. This turns out to be possible in part due to the existence of the (previously unknown) exact conservation law satisfied by the quantities describing the considered neutrino system. Interpretation of the decoherence effects in terms of neutrino wave packet separation is also given, both in the adiabatic and non-adiabatic regimes of neutrino flavour evolution.

Neutrinos play a crucial role in the collapse and explosion of massive stars, governing the infall dynamics of the stellar core, triggering and fueling the explosion and driving the cooling and deleptonization of the newly formed neutron star. Due to their role neutrinos carry information from the heart of the explosion and, due to their weakly interacting nature, offer the only direct probe of the dynamics and thermodynamics at the center of a supernova. In this paper, we review the present status of modelling the neutrino physics and signal formation in collapsing and exploding stars. We assess the capability of current and planned large underground neutrino detectors to yield faithful information of the time and flavor dependent neutrino signal from a future Galactic supernova. We show how the observable neutrino burst would provide a benchmark for fundamental supernova physics with unprecedented richness of detail. Exploiting the treasure of the measured neutrino events requires a careful discrimination of source-generated properties from signal features that originate on the way to the detector. As for the latter, we discuss self-induced flavor conversions associated with neutrino-neutrino interactions that occur in the deepest stellar regions; matter effects that modify the pattern of flavor conversions in the dynamical stellar envelope; neutrino-oscillation signatures that result from structural features associated with the shock-wave propagation as well as turbulent mass motions in post-shock layers. Finally, we highlight our current understanding of the formation of the diffuse supernova neutrino background and we analyse the perspectives for a detection of this relic signal that integrates the contributions from all past core-collapse supernovae in the Universe.

We show that a self-interacting neutrino gas can spontaneously acquire a non-stationary pulsating component in its flavor content, with a frequency that can exactly cancel the "multi-angle" refractive effects of dense matter. This can then enable homogeneous and inhomogeneous flavor conversion instabilities to exist even at large neutrino and matter densities, where the system would have been stable if the evolution were strictly stationary. Large flavor conversions, especially close to a supernova core, are possible via this novel mechanism. This may have important consequences for the explosion dynamics, nucleosynthesis, as well as for neutrino observations of supernovae.

Self-induced flavor conversions of supernova (SN) neutrinos have been characterized in the spherically symmetric "bulb" model, reducing the neutrino evolution to a one dimensional problem along a radial direction. We lift this assumption, presenting a two-dimensional toy-model where neutrino beams are launched in many different directions from a ring. We find that self-interacting neutrinos spontaneously break the spatial symmetries of this model. As a result the flavor content and the lepton number of the neutrino gas would acquire seizable direction-dependent variations, breaking the coherent behavior found in the symmetric case. This finding would suggest that the previous results of the self-induced flavor evolution obtained in one-dimensional models should be critically re-examined.

We consider a simplifed model for self-induced flavor conversions of a dense neutrino gas in two dimensions, showing new solutions that spontaneously break the spatial symmetries of the initial conditions. As a result of the symmetry breaking induced by the neutrino-neutrino interactions, the coherent behavior of the neutrino gas becomes unstable. This instability produces large spatial variations in the flavor content of the ensemble. Furthermore, it also leads to the creation of domains of different net lepton number flux. The transition of the neutrino gas from a coherent to incoherent behavior shows an intriguing analogy with a streaming flow changing from laminar to turbulent regime. These finding would be relevant for the self-induced conversions of neutrinos streaming-off a supernova core.

We revise the bound from the supernova SN1987A on the coupling of ultralight axion-like particles (ALPs) to photons. In a core-collapse supernova, ALPs would be emitted via the Primakoff process, and eventually convert into gamma rays in the magnetic field of the Milky Way. The lack of a gamma-ray signal in the GRS instrument of the SMM satellite in coincidence with the observation of the neutrinos emitted from SN1987A therefore provides a strong bound on their coupling to photons. Due to the large uncertainty associated with the current bound, we revise this argument, based on state-of-the-art physical inputs both for the supernova models and for the Milky-Way magnetic field. Furthermore, we provide major amendments, such as the consistent treatment of nucleon-degeneracy effects and of the reduction of the nuclear masses in the hot and dense nuclear medium of the supernova. With these improvements, we obtain a new upper limit on the photon-ALP coupling: g_{a\gamma} < 5.3 x 10^{-12} GeV^{-1}, for m_a < 4.4 x 10^{-10} eV, and we also give its dependence at larger ALP masses. Moreover, we discuss how much the Fermi-LAT satellite experiment could improve this bound, should a close-enough supernova explode in the near future.

Secret interactions among sterile neutrinos have been recently proposed as an escape-route to reconcile eV sterile neutrino hints from short-baseline anomalies with cosmological observations. In particular models with coupling g_X \gtrsim 10^{-2} and gauge boson mediators $X$ with $M_X \lesssim 10$ MeV lead to large matter potential suppressing the sterile neutrino production before the neutrino decoupling. With this choice of parameter ranges, big bang nucleosynthesis is left unchanged and gives no bound on the model. However, we show that at lower temperatures when active-sterile oscillations are no longer matter suppressed, sterile neutrinos are still in a collisional regime, due to their secret self-interactions. The interplay between vacuum oscillations and collisions leads to a scattering-induced decoherent production of sterile neutrinos with a fast rate. This process is responsible for a flavor equilibration among the different neutrino species. We explore the effect of this large sterile neutrino population on cosmological observables. We find that a signature of strong secret interactions would be a reduction of the effective number of neutrinos $N_{\rm eff}$ at matter radiation equality down to 2.7. Moreover, for $M_X \gtrsim g_X$ MeV sterile neutrinos would be free-streaming before becoming non-relativistic and they would affect the large-scale structure power spectrum. As a consequence, for this range of parameters we find a tension of a eV mass sterile state with cosmological neutrino mass bounds.

Short-baseline neutrino anomalies suggest the existence of low-mass ( m \sim O(1)~eV) sterile neutrinos \nu_s. These would be efficiently produced in the early universe by oscillations with active neutrino species, leading to a thermal population of the sterile states seemingly incompatible with cosmological observations. In order to relieve this tension it has been recently speculated that new "secret" interactions among sterile neutrinos, mediated by a massive gauge boson X (with M_X << M_W), can inhibit or suppress the sterile neutrino thermalization, due to the production of a large matter potential term. We note however, that they also generate strong collisional terms in the sterile neutrino sector that induce an efficient sterile neutrino production after a resonance in matter is encountered, increasing their contribution to the number of relativistic particle species N_ eff. Moreover, for values of the parameters of the \nu_s-\nu_s interaction for which the resonance takes place at temperature T\lesssim few MeV, significant distortions are produced in the electron (anti)neutrino spectra, altering the abundance of light element in Big Bang Nucleosynthesis (BBN). Using the present determination of $^4$He and deuterium primordial abundances we determine the BBN constraints on the model parameters. We find that $^2$H/H density ratio exclude much of the parameter space if one assume a baryon density at the best fit value of Planck experiment, \Omega_B h^2= 0.02207, while bounds become weaker for a higher \Omega_B h^2=0.02261, the 95 % C.L. upper bound of Planck. Due to the large error on its experimental determination, the helium mass fraction Y_p gives no significant bounds.

During the shock-wave propagation in a core-collapse supernova (SN), matter turbulence may affect neutrino flavor conversion probabilities. Such effects have been usually studied by adding parametrized small-scale random fluctuations (with arbitrary amplitude) on top of coarse, spherically symmetric matter density profiles. Recently, however, two-dimensional (2D) SN models have reached a space resolution high enough to directly trace anisotropic density profiles, down to scales smaller than the typical neutrino oscillation length. In this context, we analyze the statistical properties of a large set of SN matter density profiles obtained in a high-resolution 2D simulation, focusing on a post-bounce time (2 s) suited to study shock-wave effects on neutrino propagation on scales as small as O(100) km and possibly below. We clearly find the imprint of a broken (Kolmogorov-Kraichnan) power-law structure, as generically expected in 2D turbulence spectra. We then compute the flavor evolution of SN neutrinos along representative realizations of the turbulent matter density profiles, and observe no or modest damping of the neutrino crossing probabilities on their way through the shock wave. In order to check the effect of possibly unresolved fluctuations at scales below O(100) km, we also apply a randomization procedure anchored to the power spectrum calculated from the simulation, and find consistent results within \pm 1 sigma fluctuations. These results show the importance of anchoring turbulence effects on SN neutrinos to realistic, fine-grained SN models.

We derive a strong bound on the axion-photon coupling $g_{a\gamma}$ from the analysis of a sample of 39 Galactic Globular Clusters. As recognized long ago, the R parameter, i.e. the number ratio of stars in horizontal over red giant branch of old stellar clusters, would be reduced by the axion production from photon conversions occurring in stellar cores. In this regard we have compared the measured R with state-of-the-art stellar models obtained under different assumptions for $g_{a\gamma}$. We show that the estimated value of $g_{a\gamma}$ substantially depends on the adopted He mass fraction Y, an effect often neglected in previous investigations. Taking as benchmark for our study the most recent determinations of the He abundance in H II regions with O/H in the same range of the Galactic Globular Clusters, we obtain an upper bound $g_{a\gamma}<0.66\times 10^{-10}$ GeV$^{-1}$ at 95$\%$ confidence level. This result significantly improves the constraints from previous analyses and is currently the strongest limit on the axion-photon coupling in a wide mass range.

It has been recently discovered that removing the axial symmetry in the "multi-angle effects" associated with the neutrino-neutrino interactions for supernova (SN) neutrinos, a new multi-azimuthal-angle (MAA) instability would trigger flavor conversions in addition to the ones caused by the bimodal and multi-zenith-angle (MZA) instabilities. We investigate the dependence of the MAA instability on the original SN neutrino fluxes, performing a stability analysis of the linearized neutrino equations of motion. We compare these results with the numerical evolution of the SN neutrino non-linear equations, looking at a local solution along a specific line of sight, under the assumption that the transverse variations of the global solution are small. We also assume that self-induced conversions are not suppressed by large matter effects. We show that the pattern of the spectral crossings (energies where F_{\nu_e} = F_{\nu_x}, and F_{\bar\nu_e} = F_{\bar\nu_x}) is crucial in determining the impact of MAA effects on the flavor evolution. For neutrino spectra with a strong excess of \nu_e over \bar\nu_e, presenting only a single-crossing, MAA instabilities would trigger new flavor conversions in normal mass hierarchy. In our simplified flavor evolution scheme, these would lead to spectral swaps and splits analogous to what produced in inverted hierarchy by the bimodal instability. Conversely, in the presence of spectra with a moderate flavor hierarchy, having multiple crossing energies, MZA effects would produce a sizable delay in the onset of the flavor conversions, inhibiting the growth of the MAA instability. In this case the splitting features for the oscillated spectra in both the mass hierarchies are the ones induced by the only bimodal and MZA effects.

It has been recently pointed out that removing the axial symmetry in the "multi-angle effects" associated with the neutrino-neutrino interactions for supernova (SN) neutrinos, a new multi-azimuthal-angle (MAA) instability would arise. In particular, for a flux ordering $F_{\nu_e} > F_{\bar\nu_e} > F_{\nu_x}$, as expected during the SN accretion phase, this instability occurs in the normal neutrino mass hierarchy. However, during this phase the ordinary matter density can be larger than the neutrino one, suppressing the self-induced conversions. At this regard, we investigate the matter suppression of the MAA effects, performing a linearized stability analysis of the neutrino equations of motion, in the presence of realistic SN density profiles. We compare these results with the numerical solution of the SN neutrino non-linear evolution equations. Assuming axially symmetric distributions of neutrino momenta we find that the large matter term strongly inhibits the MAA effects. In particular, the hindrance becomes stronger including realistic forward-peaked neutrino angular distributions. As a result, in our model for a $10.8$ $M_{\odot}$ iron-core SNe, MAA instability does not trigger any flavor conversion during the accretion phase. Instead, for a $8.8$ $M_{\odot}$ O-Ne-Mg core SN model, with lower matter density profile and less forward-peaked angular distributions, flavor conversions are possible also at early times.

The most general case of self-induced neutrino flavor evolution is described by a set of kinetic equations for a dense neutrino gas evolving both in space and time. Solutions of these equations have been typically worked out assuming that either the time (in the core-collapse supernova environment) or space (in the early universe) homogeneity in the initial conditions is preserved through the evolution. In these cases one can gauge away the homogeneous variable and reduce the dimensionality of the problem. In this paper we investigate if small deviations from an initial postulated homogeneity can be amplified by the interacting neutrino gas, leading to a new flavor instability. To this end, we consider a simple two flavor isotropic neutrino gas evolving in time, and initially composed by only $\nu_e$ and $\bar\nu_e$ with equal densities. In the homogeneous case, this system shows a bimodal instability in the inverted mass hierarchy scheme, leading to the well studied flavor pendulum behavior. This would lead to periodic pair conversions $\nu_e \bar\nu_e \leftrightarrow \nu_x \bar\nu_x$. To break space homogeneity, we introduce small amplitude space-dependent perturbations in the matter potential. By Fourier transforming the equations of motion with respect to the space coordinate, we then numerically solve a set of coupled equations for the different Fourier modes. We find that even for arbitrarily tiny inhomogeneities, the system evolution runs away from the stable pendulum behavior: the different modes are excited and the space-averaged ensemble evolves towards flavor equilibrium. We finally comment on the role of a time decaying neutrino background density in weakening these results.

We use cosmological observations in the post-Planck era to derive limits on thermally produced cosmological axions. In the early universe such axions contribute to the radiation density and later to the hot dark matter fraction. We find an upper limit m_a < 0.67 eV at 95% C.L. after marginalising over the unknown neutrino masses, using CMB temperature and polarisation data from Planck and WMAP respectively, the halo matter power spectrum extracted from SDSS-DR7, and the local Hubble expansion rate H_0 released by the Carnegie Hubble Program based on a recalibration of the Hubble Space Telescope Key Project sample. Leaving out the local H_0 measurement relaxes the limit somewhat to 0.86 eV, while Planck+WMAP alone constrain the axion mass to 1.01 eV, the first time an upper limit on m_a has been obtained from CMB data alone. Our axion limit is therefore not very sensitive to the tension between the Planck-inferred H_0 and the locally measured value. This is in contrast with the upper limit on the neutrino mass sum, which we find here to range from 0.27 eV at 95% C.L. combining all of the aforementioned observations, to 0.84 eV from CMB data alone.

The flavor evolution of neutrinos emitted by a supernova (SN) core is strongly affected by the refractive effects associated with the neutrino-neutrino interactions in the deepest stellar regions. Till now, all numerical studies have assumed the axial symmetry for the "multi-angle effects" associated with the neutrino-neutrino interactions. Recently, it has been pointed out in \cite{Raffelt:2013rqa} that removing this assumption, a new multi-azimuthal-angle (MAA) instability would emerge in the flavor evolution of the dense SN neutrino gas in addition to the one caused by multi-zenith-angle (MZA) effects. Inspired by this result, for the first time we numerically solve the non-linear neutrino propagation equations in SN, introducing the azimuthal angle as angular variable in addition to the usual zenith angle. We consider simple energy spectra with an excess of \nu_e over anti-\nu_e. In these cases, we find that even starting with a complete axial symmetric neutrino emission, the MAA effects would lead to significant flavor conversions in normal mass hierarchy, in cases otherwise stable under the only MZA effects. The final outcome of the flavor conversions, triggered by the MAA instability, depends on the initial asymmetry between \nu_e and anti-\nu_e spectra. If it is sufficiently large, final spectra would show an ordered behavior with spectral swaps and splits. Conversely, for small flavor asymmetries flavor decoherence among angular modes develops, affecting the flavor evolution also in inverted mass hierarchy.

It has been speculated that Lorentz-invariance violation (LIV) might be generated by quantum-gravity (QG) effects. As a consequence, particles may not travel at the universal speed of light. In particular, superluminal extragalactic neutrinos would rapidly lose energy via the bremssthralung of electron-positron pairs (nu -> nu e+ e-), damping their initial energy into electromagnetic cascades, a figure constrained by Fermi-LAT data. We show that the two cascade neutrino events with energies around 1 PeV recently detected by IceCube -if attributed to extragalactic diffuse events, as it appears likely- can place the strongest bound on LIV in the neutrino sector, namely delta =(v^2-1) < O(10^(-18)), corresponding to a QG scale M_QG ~ 10^5 M_Pl (M_QG >~ 10^(-4) M_Pl) for a linear (quadratic) LIV, at least for models inducing superluminal neutrino effects (delta > 0).

Light sterile neutrinos can be excited by oscillations with active neutrinos in the early universe. Their properties can be constrained by their contribution as extra-radiation, parameterized in terms of the effective number of neutrino species N_ eff, and to the universe energy density today \Omega_\nu h^2. Both these parameters have been measured to quite a good precision by the Planck satellite experiment. We use this result to update the bounds on the parameter space of (3+1) sterile neutrino scenarios, with an active-sterile neutrino mass squared splitting in the range (10^{-5} - 10^2 ) eV^2. We consider both normal and inverted mass orderings for the active and sterile states. For the first time we take into account the possibility of two non-vanishing active-sterile mixing angles. We find that the bounds are more stringent than those obtained in laboratory experiments. This leads to a strong tension with the short-baseline hints of light sterile neutrinos. In order to relieve this disagreement, modifications of the standard cosmological scenario, e.g. large primordial neutrino asymmetries, are required.

We perform a study of the flavour evolution in the early universe of a multi-flavour active-sterile neutrino system with parameters inspired by the short-baseline neutrino anomalies. In a neutrino-symmetric bath a "thermal" population of the sterile state would quickly grow, but in the presence of primordial neutrino asymmetries a self-suppression as well as a resonant sterile neutrino production can take place, depending on temperature and chosen parameters. In order to characterize these effects, we go beyond the usual average momentum and single mixing approximations and consider a multi-momentum and multi-flavour treatment of the kinetic equations. We find that the enhancement obtained in this case with respect to the average momentum approximation is significant, up to \sim 20 % of a degree of freedom. Such detailed and computationally demanding treatment further raises the asymmetry values required to significantly suppress the sterile neutrino production, up to large and preferentially net asymmetries |L_{\nu}| > O(10^{-2}). For such asymmetries, however, the active-sterile flavour conversions happen so late that significant distortions are produced in the electron (anti)neutrino spectra. The larger |L_{\nu}|, the more the impact of these distortions takes over as dominant cosmological effect, notably increasing the 4 He abundance in primordial nucleosynthesis (BBN). The standard expression of the primordial yields in terms of the effective number of neutrinos and asymmetries is also greatly altered. We numerically estimate the magnitude of such effects for a few representative cases and comment on possible implications for forthcoming cosmological measurements.

Quantum-gravity (QG) effects might generate Lorentz invariance violation by the interaction of energetic particles with the foamy structure of the space-time. As a consequence, particles may not travel at the universal speed of light. We propose to constrain Lorentz invariance violation for energetic neutrinos exploiting the $\nu_e$ neutronization burst from the next galactic supernova (SN). This prompt signal is expected to produce a sharp peak in the SN $\nu_e$ light curve with a duration of $\sim 25$ ms. However presence of energy-dependent Lorentz invariance violation would significantly spread out the time structure of this signal. We find that the detection the SN $\nu_e$ burst from a typical galactic explosion at $d=10$ kpc in a Mton-class water Cerenkov detector, would be sensitive to a quantum-gravity mass scale $M_{\rm QG} \sim 10^{12}$ GeV ($2 \times10^{5}$ GeV) for the linear (quadratic) energy dependence of Lorentz invariance violation. These limits are rather independent of the neutrino mass hierarchy and whether the neutrino velocity is super or subluminal.

Very high energy (VHE, energy >~ 100 GeV) {\gamma}-rays originating from extragalactic sources interact with low energy photons of background radiation fields and produce electron-positron pairs. Alternatively, in the presence of ambient magnetic fields, they can convert into hypothetical axion-like particles (ALPs), pseudo-scalar spin-0 bosons, predicted by extensions of the standard model. These particles propagate unimpeded over cosmological distances. Here, the effect of photon-ALP oscillations in magnetic fields of galaxy clusters and the Milky Way on VHE {\gamma}-ray spectra is studied. It is shown that this mechanism can lead to a substantial enhancement of the VHE flux and a spectral hardening, thus effectively reducing the opacity of the Universe to VHE {\gamma}-rays.

A fraction of AGN producing VHE gamma-rays are located in galaxy clusters. The magnetic field present in the intra-cluster medium would lead to conversions of VHE photons into axion-like particles (ALPs), which are a generic prediction of several extensions of the Standard Model. ALPs produced in this way would traverse cosmological distances unaffected by the extragalactic background light at variance with VHE photons which undergo a substantial absorption. Eventually, a nontrivial fraction of ALPs would re-convert into VHE photons in the magnetic field of the Milky Way. This mechanism produces a significant hardening of the VHE spectrum of AGN in galaxy clusters. As a specific example we consider the energy spectra of two observed VHE gamma-ray sources located in galaxy clusters, namely 1ES 0414+009 at redshift z=0.287 and Mkn 501 at z=0.034. We find that the hardening in the observed spectra becomes relevant at E > 1 TeV. The detection of this signature would allow to indirectly probe the existence of ultra-light ALPs with mass m_a < 10^{-8} eV and photon-ALP coupling g_{a gamma} < 10^{-10} GeV^{-1} with the presently operating Imaging Atmospheric Cherenkov Telescopes like H.E.S.S., MAGIC, VERITAS and CANGAROO-III and even more likely with the planned detectors like CTA, HAWC and HiSCORE. An independent laboratory check of ultra-light ALPs invoked in this mechanism can be performed with the planned upgrade of the photon regeneration experiment ALPS at DESY and with the next generation solar axion detector IAXO.

Neutrino oscillations in the Earth matter may introduce peculiar modulations in the supernova (SN) neutrino spectra. The detection of this effect has been proposed as diagnostic tool for the neutrino mass hierarchy at "large" 1-3 leptonic mixing angle theta13. We perform an updated study on the observability of this effect at large next-generation underground detectors (i.e., 0.4 Mton water Cherenkov, 50 kton scintillation and 100 kton liquid Argon detectors) based on neutrino fluxes from state-of-the-art SN simulations and accounting for statistical fluctuations via Montecarlo simulations. Since the average energies predicted by recent simulations are lower than previously expected and a tendency towards the equalization of the neutrino fluxes appears during the SN cooling phase, the detection of the Earth matter effect will be more challenging than expected from previous studies. We find that none of the proposed detectors shall be able to detect the Earth modulation for the neutrino signal of a typical galactic SN at 10 kpc. It should be observable in a 100 kton liquid Argon detector for a SN at few kpc and all three detectors would clearly see the Earth signature for very close-by stars only (d ~ 0.2 kpc). Finally, we show that adopting IceCube as co-detector together with a Mton water Cherenkov detector is not a viable option either.

Numerical simulations of the supernova (SN) neutrino self-induced flavor conversions, associated with the neutrino-neutrino interactions in the deepest stellar regions, have been typically carried out assuming the "bulb-model". In this approximation, neutrinos are taken to be emitted half-isotropically by a common neutrinosphere. In the recent Ref. \cite{Mirizzi:2011tu} we have removed this assumption by introducing flavor-dependent angular distributions for SN neutrinos, as suggested by core-collapse simulations. We have found that in this case a novel multi-angle instability in the self-induced flavor transitions can arise. In this work we perform an extensive study of this effect, carrying out a linearized flavor stability analysis for different SN neutrino energy fluxes and angular distributions, in both normal and inverted neutrino mass hierarchy. We confirm that spectra of different nu species which cross in angular space (where F_{\nu_e}=F_{\nu_x} and F_{\bar\nu_e}=F_{\bar\nu_x}) present a significant enhancement of the flavor instability, and a shift of the onset of the flavor conversions at smaller radii with respect to the case of an isotropic neutrino emission. We also illustrate how a qualitative (and sometimes quantitative) understanding of the dynamics of these systems follows from a stability analysis.

Light sterile neutrinos mixing with the active ones have been recently proposed to solve different anomalies observed in short-baseline oscillation experiments. These neutrinos can also be produced by oscillations of the active neutrinos in the early universe, leaving possible traces on different cosmological observables. Here we perform an updated study of the neutrino kinetic equations in (3+1) and (2+1) oscillation schemes, dynamically evolving primordial asymmetries of active neutrinos and taking into account for the first time CP-violation effects. In the absence of neutrino asymmetries, eV-mass scale sterile neutrinos would be completely thermalized creating a tension with respect to the CMB, LSS and BBN data. In the past literature, active neutrino asymmetries have been invoked as a way to inhibit the sterile neutrino production via the in-medium suppression of the sterile-active mixing angle. However, neutrino asymmetries also permit a resonant sterile neutrino production. We find that if the active species have equal asymmetries L, a value |L|=10^{-3} is required to start suppressing the resonant sterile production, roughly an order of magnitude larger than what previously expected. When active species have opposite asymmetries the sterile abundance is further enhanced, requiring an even larger |L|\simeq 10^{-2} to start suppressing their production. In the latter case, CP-violation (naturally expected) further exacerbates the phenomenon. Some consequences for cosmological observables are briefly discussed: for example, it is likely that moderate suppressions of the sterile species production are associated with significant spectral distortions of the active neutrino species, with potentially interesting phenomenological consequences especially for BBN.

We revisit our previous results on the matter suppression of self-induced neutrino flavor conversions during a supernova (SN) accretion phase, performing a linearized stability analysis of the neutrino equations of motion, in the presence of realistic SN density profiles. In our previous numerical study, we used a simplified model based on an isotropic neutrino emission with a single typical energy. Here, we take into account realistic neutrino energy and angle distributions. We find that multi-energy effects have a sub-leading impact in the flavor stability of the SN neutrino fluxes with respect to our previous single-energy results. Conversely, realistic forward-peaked neutrino angular distributions would enhance the matter suppression of the self-induced oscillations with respect to an isotropic neutrino emission. As a result, in our models for iron-core SNe, collective flavor conversions have a negligible impact on the characterization of the observable neutrino signal during the accretion phase. Instead, for a low-mass O-Ne-Mg core SN model, with lower matter density profile and less forward-peaked angular distributions, collective conversions are possible also at early times.

The usual description of self-induced flavor conversions for neutrinos (\nu's) in supernovae is based on the simplified assumption that all the \nu's of the different species are emitted "half-isotropically" by a common neutrinosphere, in analogy to a blackbody emission. However, realistic supernova simulations show that \nu angular distributions at decoupling are far from being half-isotropic and, above all, are flavor-dependent. We show that flavor-dependent angular distributions may lead to crossing points in the angular spectra of different \nu species (where F_{\nu_e}=F_{\nu_x} and F_{{\bar\nu}_e}=F_{{\bar\nu}_x}) around which a new multi-angle instability can develop. To characterize this effect, we carry out a linearized flavor stability analysis for different SN neutrino angular distributions. We find that this instability can shift the onset of the flavor conversions toward low-radii and produce a smearing of the splitting features found with trivial $\nu$ emission models. As a result the spectral differences among \nu's of different flavors could be strongly reduced.

The rise time of a Galactic supernova (SN) bar-nue lightcurve, observable at a high-statistics experiment such as the IceCube Cherenkov detector, can provide a diagnostic tool for the neutrino mass hierarchy at "large" 1-3 leptonic mixing angle theta_13. Thanks to the combination of matter suppression of collective effects at early postbounce times on one hand and the presence of the ordinary Mikheyev-Smirnov-Wolfenstein effect in the outer layers of the SN on the other hand, a sufficiently fast rise time on O(100) ms scale is indicative of an inverted mass hierarchy. We investigate results from an extensive set of stellar core-collapse simulations, providing a first exploration of the astrophysical robustness of these features. We find that for all the models analyzed (sharing the same weak interaction microphysics) the rise times for the same hierarchy are similar not only qualitatively, but also quantitatively, with the signals for the two classes of hierarchies significantly separated. We show via Monte Carlo simulations that the two cases should be distinguishable at IceCube for SNe at a typical Galactic distance 99% of the times. Finally, a preliminary survey seems to show that the faster rise time for inverted hierarchy as compared to normal hierarchy is a qualitatively robust feature predicted by several simulation groups. Since the viability of this signature ultimately depends on the quantitative assessment of theoretical/numerical uncertainties, our results motivate an extensive campaign of comparison of different code predictions at early accretion times with implementation of microphysics of comparable sophistication, including effects such like nucleon recoils in weak interactions.

The mixing of photons with axion-like particles (ALPs) in the large-scale magnetic field $B$ changes the polarization angle of a linearly polarized photon beam from active galactic nuclei in radio galaxies as it propagates over cosmological distances. Using available ultraviolet polarization data concerning these sources we derive a new bound on the product of the photon-ALP coupling $g_{a\gamma}$ times $B$. We find $g_{a\gamma} B \lesssim 10^{-11}$ GeV$^{-1}$ nG for ultralight ALPs with $m_a \lesssim 10^{-15}$ eV. We compare our new bound with the ones present in the literature and we comment about possible improvements with observations of more sources.

We propose the liquid-scintillator detector LENA (Low Energy Neutrino Astronomy) as a next-generation neutrino observatory on the scale of 50 kt. The outstanding successes of the Borexino and KamLAND experiments demonstrate the large potential of liquid-scintillator detectors in low-energy neutrino physics. LENA's physics objectives comprise the observation of astrophysical and terrestrial neutrino sources as well as the investigation of neutrino oscillations. In the GeV energy range, the search for proton decay and long-baseline neutrino oscillation experiments complement the low-energy program. Based on the considerable expertise present in European and international research groups, the technical design is sufficiently mature to allow for an early start of detector realization.

The large neutrino fluxes emitted with a distinct flavor hierarchy from core-collapse supernovae (SNe) during the post-bounce accretion phase, offer the best opportunity to detect effects from neutrino flavor oscillations. We perform a dedicated study of the SN neutrino flavor evolution during the accretion phase, using results from recent neutrino radiation hydrodynamics simulations. In contrast to what expected in the presence of only neutrino-neutrino interactions, we find that the multi-angle effects associated with the dense ordinary matter suppress collective oscillations. This is related to the high matter densities during the accretion phase in core-collapse SNe of massive iron-core progenitors. The matter suppression implies that neutrino oscillations will start outside the neutrino transport region and therefore will have a negligible impact on the neutrino heating and the explosion dynamics. Furthermore, the possible detection of the next galactic SN neutrino signal from the accretion phase, based on the usual Mikheyev- Smirnov-Wolfenstein effect in the SN mantle and Earth matter effects, can reveal the neutrino mass hierarchy in the case that the mixing angle $\theta_{13}$ is not very small.

The non-isotropic nature of the neutrino emission from a supernova (SN) core might potentially affect the flavor evolution of the neutrino ensemble, via neutrino-neutrino interactions in the deepest SN regions. We investigate the dependence of these "multi-angle effects" on the original SN neutrino fluxes in a three-flavor framework. We show that the pattern of the spectral crossings (energies where F_\nu_e= F_\nu_x, and F_\bar\nu_e= F_\bar\nu_x) is crucial in determining the impact of multi-angle effects on the flavor evolution. For neutrino spectra presenting only a single-crossing, synchronization of different angular modes prevails over multi-angle effects, producing the known "quasi single-angle" evolution. Conversely, in the presence of spectra with multiple crossing energies, synchronization is not stable. In this situation, multi-angle effects would produce a sizable delay in the onset of the flavor conversions, as recently observed. We show that, due to the only partial adiabaticity of the evolution at large radii, the multi-angle suppression can be so strong to dramatically affect the final oscillated neutrino spectra. In particular three-flavor effects, associated with the solar parameters, could be washed-out in multi-angle simulations.

The usual description of self-induced neutrino flavor conversions in core collapse supernovae (SNe) is based on the dominance of the neutrino density n_nu over the net electron density n_e. However, this condition is not met during the post-bounce accretion phase, when the dense matter in a SN is piled up above the neutrinosphere. As recently pointed-out, a dominant matter term in the anisotropic SN environment would dephase the flavor evolution for neutrinos traveling on different trajectories, challenging the occurrence of the collective behavior in the dense neutrino gas. Using the results from recent long term simulations of core-collapse SN explosions, based on three flavor Boltzmann neutrino transport in spherical symmetry, we find that both the situations of complete matter suppression (when n_e >> n_nu) and matter-induced decoherence (when n_e \gtrsim n_nu) of flavor conversions are realized during the accretion phase. The matter suppression at high densities prevents any possible impact of the neutrino oscillations on the neutrino heating and hence on the dynamics of the explosion. Furthermore, it changes the interpretation of the Earth matter effect on the SN neutrino signal during the accretion phase, allowing the possibility of the neutrino mass hierarchy discrimination at not too small values of the leptonic mixing angle \theta_{13} (i.e. \sin^2{\theta}_{13} \gtrsim 10^{-3}).

We calculate the expected galactic supernova neutrino signal at large next-generation underground detectors. At different epochs after the explosion, the primary fluxes can be quite different. For these primary neutrino fluxes, spectral splits induced by collective neutrino flavor transformations can arise for either mass hierarchy in both neutrino and antineutrino channels. We classify flux models according to the nature and number of these splits, and calculate the observable electron-neutrino and electron-antineutrino spectra at Earth, taking into account subsequent matter effects. We find that some of the spectral splits could occur sufficiently close to the peak energies to produce significant distortions in the observable SN neutrino signal. The most striking signature of this effect would be presence of peculiar energy dependent modulations associated with Earth matter crossing, present only in portions of the SN neutrino energy spectra demarcated by spectral splits. These signatures at proposed large water Cherenkov, scintillation, and liquid Argon detectors could give hints about the primary SN neutrino fluxes, as well as on the neutrino mass hierarchy and the mixing angle theta_{13}.

We update cosmological hot dark matter constraints on neutrinos and hadronic axions. Our most restrictive limits use 7-year data from the Wilkinson Microwave Anisotropy Probe for the cosmic microwave background anisotropies, the halo power spectrum (HPS) from the 7th data release of the Sloan Digital Sky Survey, and the Hubble constant from Hubble Space Telescope observations. We find 95% C.L. upper limits of \sum m_\nu<0.44 eV (no axions), m_a<0.91 eV (assuming \sum m_\nu=0), and \sum m_\nu<0.41 eV and m_a<0.72 eV for two hot dark matter components after marginalising over the respective other mass. CMB data alone yield \sum m_\nu<1.19 eV (no axions), while for axions the HPS is crucial for deriving m_a constraints. This difference can be traced to the fact that for a given hot dark matter fraction axions are much more massive than neutrinos.

Predictions of the thermodynamic conditions for phase transitions at high baryon densities and large chemical potentials are currently uncertain and largely phenomenological. Neutrino observations of core-collapse supernovae can be used to constrain the situation. Recent simulations of stellar core collapse that include a description of quark matter predict a sharp burst of anti \nu_e several hundred milliseconds after the prompt \nu_e neutronization burst. We study the observational signatures of that anti \nu_e burst at current neutrino detectors - IceCube and Super-Kamiokande. For a Galactic core-collapse supernova, we find that signatures of the QCD phase transition can be detected, regardless of the neutrino oscillation scenario. The detection would constitute strong evidence of a phase transition in the stellar core, with implications for the equation of state at high matter density and the supernova explosion mechanism.

It was recently realized that three-flavor effects could peculiarly modify the development of spectral splits induced by collective oscillations, for supernova neutrinos emitted during the cooling phase of a protoneutron star. We systematically explore this case, explaining how the impact of these three-flavor effects depends on the ordering of the neutrino masses. In inverted mass hierarchy, the solar mass splitting gives rise to instabilities in regions of the (anti)neutrino energy spectra that were otherwise stable under the leading two-flavor evolution governed by the atmospheric mass splitting and by the 1-3 mixing angle. As a consequence, the high-energy spectral splits found in the electron (anti)neutrino spectra disappear, and are transferred to other flavors. Imperfect adiabaticity leads to smearing of spectral swap features. In normal mass hierarchy, the three-flavor and the two-flavor instabilities act in the same region of the neutrino energy spectrum, leading to only minor departures from the two-flavor treatment.

Various satellite-borne missions are being planned whose goal is to measure the polarization of a large number of gamma-ray bursts (GRBs). We show that the polarization pattern predicted by current models of GRB emission can be drastically modified by the existence of very light axion-like particles (ALPs), which are present in many extensions of the Standard Model of particle physics. Basically, the propagation of photons emitted by a GRB through cosmic magnetic fields with a domain-like structure induces photon-ALP mixing, which is expected to produce a strong modification of the original photon polarization. Because of the random orientation of the magnetic field in each domain, this effect strongly depends on the orientation of the photon line of sight. As a consequence, photon-ALP conversion considerably broadens the original polarization distribution. Searching for such a peculiar feature through future high-statistics polarimetric measurements is therefore a new opportunity to discover very light ALPs.

Very-high energy photons emitted by distant cosmic sources are absorbed on the extragalactic background light (EBL) during their propagation. This effect can be characterized in terms of a photon transfer function at Earth. The presence of extragalactic magnetic fields could also induce conversions between very high-energy photons and hypothetical axion-like particles (ALPs). The turbulent structure of the extragalactic magnetic fields would produce a stochastic behaviour in these conversions, leading to a statistical distribution of the photon transfer functions for the different realizations of the random magnetic fields. To characterize this effect, we derive new equations to calculate the mean and the variance of this distribution. We find that, in presence of ALP conversions, the photon transfer functions on different lines of sight could have relevant deviations with respect to the mean value, producing both an enhancement or a suppression in the observable photon flux with respect to the expectations with only absorption. As a consequence, the most striking signature of the mixing with ALPs would be a reconstructed EBL density from TeV photon observations which appears to vary over different directions of the sky: consistent with standard expectations in some regions, but inconsistent in others.

We discuss current cosmological constraints on axions, as well as future sensitivities. Bounds on axion hot dark matter are discussed first, and subsequently we discuss both current and future sensitivity to models in which axions play the role as cold dark matter, but where the Peccei-Quinn symmetry is not restored during reheating.

In the framework of the Standard Model we calculate the flavour non-universal correction for neutrino refraction in a neutrino background and verify a similar previous result for the case of ordinary-matter background. The dominant term arises at loop level and involves tau leptons circulating in the loop. These O(G_F m_tau^2) corrections to the tree-level potential provide the dominant refractive difference between nu_mu and nu_tau unless the medium contains mu or tau leptons. Our results affect the flavour evolution of dense neutrino gases and may be of interest for collective three-flavour oscillations of supernova neutrinos. We spell out explicitly how these non-universal neutrino-neutrino interactions enter the flavour oscillation equations.

The presence of a primordial magnetic field would have induced resonant conversions between photons and axion-like particles (ALPs) during the thermal history of the Universe. These conversions would have distorted the blackbody spectrum of the cosmic microwave background (CMB). In this context, we derive bounds on the photon-ALP resonant conversions using the high precision CMB spectral data collected by the FIRAS instrument on board of the Cosmic Background Explorer. We obtain upper limits on the product of the photon-ALP coupling constant g times the magnetic field strength B down to g B > 10^{-13} GeV^{-1} nG for ALP masses below the eV scale.

Various extensions of the Standard Model predict the existence of hidden photons kinetically mixing with the ordinary photon. This mixing leads to oscillations between photons and hidden photons, analogous to the observed oscillations between different neutrino flavors. In this context, we derive new bounds on the photon-hidden photon mixing parameters using the high precision cosmic microwave background spectral data collected by the Far Infrared Absolute Spectrophotometer instrument on board of the Cosmic Background Explorer. Requiring the distortions of the CMB induced by the photon-hidden photon mixing to be smaller than experimental upper limits, this leads to a bound on the mixing angle < 10^{-7}-10^{-5} for hidden photon masses between 10^{-14} eV and 10^{-7} eV. This low-mass and low-mixing region of the hidden photon parameter space was previously unconstrained.

We give a very brief overview of collective effects in neutrino oscillations in core collapse supernovae where refractive effects of neutrinos on themselves can considerably modify flavor oscillations, with possible repercussions for future supernova neutrino detection. We discuss synchronized and bipolar oscillations, the role of energy and angular neutrino modes, as well as three-flavor effects. We close with a short summary and some open questions.

We investigate the effects of nonstandard four-fermion neutrino-neutrino interactions on the flavor evolution of dense neutrino gases. We find that in the regions where the neutrino-neutrino refractive index leads to collective flavor oscillations, the presence of new neutrino interactions can produce flavor equilibration in both normal and inverted neutrino mass hierarchy. In realistic supernova environments, these effects are significant if the nonstandard neutrino-neutrino interaction strength is comparable to the one expected in the standard case, dominating the ordinary matter potential. However, very small nonstandard neutrino-neutrino couplings are enough to trigger the usual collective neutrino flavor transformations in the inverted neutrino mass hierarchy, even if the mixing angle vanishes exactly.

For neutrinos streaming from a supernova (SN) core, dense matter suppresses self-induced flavor transformations if the electron density n_e significantly exceeds the neutrino density n_nu in the conversion region. If n_e is comparable to n_nu one finds multi-angle decoherence, whereas the standard self-induced transformation behavior requires that in the transformation region n_nu is safely above n_e. This condition need not be satisfied in the early phase after supernova core bounce. Our new multi-angle effect is a subtle consequence of neutrinos traveling on different trajectories when streaming from a source that is not point-like.

In core-collapse supernovae, neutrinos and antineutrinos are initially subject to significant self-interactions induced by weak neutral currents, which may induce strong-coupling effects on the flavor evolution (collective transitions). The interpretation of the effects is simplified when self-induced collective transitions are decoupled from ordinary matter oscillations, as for the matter density profile that we discuss. In this case, approximate analytical tools can be used (pendulum analogy, swap of energy spectra). For inverted neutrino mass hierarchy, the sequence of effects involves: synchronization, bipolar oscillations, and spectral split. Our simulations shows that the main features of these regimes are not altered when passing from simplified (angle-averaged) treatments to full, multi-angle numerical experiments.

In the dense supernova core, self-interactions may align the flavor polarization vectors of neutrinos and antineutrinos, and induce collective flavor transformations. Different alignment ansatzes are known to describe approximately the phenomena of synchronized or bipolar oscillations, and the split of neutrino energy spectra. We discuss another phenomenon observed in some numerical experiments in inverted hierarchy, showing features akin to a low-energy split of antineutrino spectra. The phenomenon appears to be approximately described by another alignment ansatz which, in the considered scenario, reduces the (nonadiabatic) dynamics of all energy modes to only two neutrino plus two antineutrino modes. The associated spectral features, however, appear to be fragile when passing from single- to multi-angle simulations.

The rich phenomenology of collective neutrino oscillations has been studied only in one-dimensional or spherically symmetric systems. Motivated by the non-spherical example of coalescing neutron stars, presumably the central engines of short gamma-ray bursts, we use the Liouville equation to formulate the problem for general source geometries. Assuming the neutrino ensemble displays self-maintained coherence, the problem once more becomes effectively one-dimensional along the streamlines of the overall neutrino flux. This approach for the first time provides a formal definition of the ``single-angle approximation'' frequently used for supernova neutrinos and allows for a natural generalization to non-spherical geometries. We study the explicit example of a disk-shaped source as a proxy for coalescing neutron stars.

The prompt nu_e burst from a core-collapse supernova (SN) is subject to both matter-induced flavor conversions and strong neutrino-neutrino refractive effects. For the lowest-mass progenitors, leading to O-Ne-Mg core SNe, the matter density profile can be so shallow that the usual MSW matter effects occur within the dense-neutrino region close to the neutrino sphere. In this case a ``split'' occurs in the emerging spectrum, i.e., the nu_e flavor survival probability shows a step-like feature. We explain this feature analytically as a ``MSW prepared spectral split.'' In a three-flavor treatment, the step-like feature actually consists of two narrowly spaced splits. They are determined by two combinations of flavor-lepton numbers that are conserved under collective oscillations.

We update our previous constraints on two-component hot dark matter (axions and neutrinos), including the recent WMAP 5-year data release. Marginalising over sum m_nu provides m_a < 1.02 eV (95% C.L.) for the axion mass. In the absence of axions we find sum m_nu < 0.63 eV (95% C.L.).

Non-linear effects on supernova neutrino oscillations, associated with neutrino self-interactions, are known to induce collective flavor transitions near the supernova core for theta_13 \neq 0. In scenarios with very shallow electron density profiles, these transformations have been shown to couple with ordinary matter effects, jointly producing spectral distortions both in normal and inverted hierarchy. In this work we consider a complementary scenario, characterized by higher electron density, as indicated by post-bounce shock-wave simulations. In this case, early collective flavor transitions are decoupled from later, ordinary matter effects. Moreover, such transitions become more amenable to both numerical computations and analytical interpretations in inverted hierarchy, while they basically vanish in normal hierarchy. We numerically evolve the neutrino density matrix in the region relevant for self-interaction effects. In the approximation of averaged intersection angle between neutrino trajectories, our simulations neatly show the collective phenomena of synchronization, bipolar oscillations, and spectral split, recently discussed in the literature. In the more realistic (but computationally demanding) case of non-averaged neutrino trajectories, our simulations do not show new significant features, apart from the smearing of ``fine structures'' such as bipolar nutations. Our results seem to suggest that, at least for non-shallow matter density profiles, averaging over neutrino trajectories plays a minor role in the final outcome. In this case, the swap of nu_e and nu_{\mu,\tau} spectra above a critical energy may represent an unmistakable signature of the inverted hierarchy, especially for theta_{13} small enough to render further matter effects irrelevant.

Neutrino oscillation experiments and direct bounds on absolute masses constrain neutrino mass differences to fall into the microwave energy range, for most of the allowed parameter space. As a consequence of these recent phenomenological advances, older constraints on radiative neutrino decays based on diffuse background radiations and assuming strongly hierarchical masses in the eV range are now outdated. We thus derive new bounds on the radiative neutrino lifetime using the high precision cosmic microwave background spectral data collected by the Far Infrared Absolute Spectrophotometer instrument on board of Cosmic Background Explorer. The lower bound on the lifetime is between a few x 10^19 s and 5 x 10^20 s, depending on the neutrino mass ordering and on the absolute mass scale. However, due to phase space limitations, the upper bound in terms of the effective magnetic moment mediating the decay is not better than ~ 10^-8 Bohr magnetons. We also comment about possible improvements of these limits, by means of recent diffuse infrared photon background data. We compare these bounds with pre-existing limits coming from laboratory or astrophysical arguments. We emphasize the complementarity of our results with others available in the literature.

One interpretation of the unexplained signature observed in the PVLAS experiment invokes a new axion-like particle (ALP) with a two-photon vertex, allowing for photon-ALP oscillations in the presence of magnetic fields. In the range of masses and couplings suggested by PVLAS, the same effect would lead to a peculiar dimming of high-energy photon sources. For typical parameters of the turbulent magnetic field in the galaxy, the effect sets in at E_gamma >~ 10 TeV, providing an ALP signature in the spectra of TeV gamma sources that can be probed with Cherenkov telescopes. A dedicated search will be strongly motivated if the ongoing photon regeneration experiments confirm the PVLAS particle interpretation.

We use observations of the cosmological large-scale structure to derive limits on two-component hot dark matter consisting of mass-degenerate neutrinos and hadronic axions, both components having velocity dispersions corresponding to their respective decoupling temperatures. We restrict the data samples to the safely linear regime, in particular excluding the Lyman-alpha forest. Using standard Bayesian inference techniques we derive credible regions in the two-parameter space of m_a and sum(m_nu). Marginalising over sum(m_nu) provides m_a < 1.2 eV (95% C.L.). In the absence of axions the same data and methods give sum(m_nu) < 0.65 eV (95% C.L.). We also derive limits on m_a for a range of axion-pion couplings up to one order of magnitude larger or smaller than the hadronic value.

The Sun is a source of high energy neutrinos (E > 10 GeV) produced by cosmic ray interactions in the solar atmosphere. We study the impact of three-flavor oscillations (in vacuum and in matter) on solar atmosphere neutrinos, and calculate their observable fluxes at Earth, as well as their event rates in a kilometer-scale detector in water or ice. We find that peculiar three-flavor oscillation effects in matter, which can occur in the energy range probed by solar atmosphere neutrinos, are significantly suppressed by averaging over the production region and over the neutrino and antineutrino components. In particular, we find that the relation between the neutrino fluxes at the Sun and at the Earth can be approximately expressed in terms of phase-averaged ``vacuum'' oscillations, dominated by a single mixing parameter (the angle theta_23).

Photon-axion conversion induced by intergalactic magnetic fields causes an apparent dimming of distant sources, notably of cosmic standard candles such as supernovae of type Ia (SNe Ia). We review the impact of this mechanism on the luminosity-redshift relation of SNe Ia, on the dispersion of quasar spectra, and on the spectrum of the cosmic microwave background. The original idea of explaining the apparent dimming of distant SNe Ia without cosmic acceleration is strongly constrained by these arguments. However, the cosmic equation of state extracted from the SN Ia luminosity-redshift relation remains sensitive to this mechanism. For example, it can mimic phantom energy.

Supernova neutrino flavor transitions during the shock wave propagation are known to encode relevant information not only about the matter density profile but also about unknown neutrino properties, such as the mass hierarchy (normal or inverted) and the mixing angle theta_13. While previous studies have focussed on "deterministic" density profiles, we investigate the effect of possible stochastic matter density fluctuations in the wake of supernova shock waves. In particular, we study the impact of small-scale fluctuations on the electron (anti)neutrino survival probability, and on the observable spectra of inverse-beta-decay events in future water-Cherenkov detectors. We find that such fluctuations, even with relatively small amplitudes, can have significant damping effects on the flavor transition pattern, and can partly erase the shock-wave imprint on the observable time spectra, especially for sin^2(theta_13) > O(10^-3).

A model-independent experimental signature for flavor oscillations in the neutrino signal from the next Galactic supernova (SN) would be the observation of Earth matter effects. We calculate the probability for observing a Galactic SN shadowed by the Earth as a function of the detector's geographic latitude. This probability depends only mildly on details of the Galactic SN distribution. A location at the North Pole would be optimal with a shadowing probability of about 60%, but a far-northern location such as Pyhasalmi in Finland, the proposed site for a large-volume scintillator detector, is almost equivalent (58%). We also consider several pairs of detector locations and calculate the probability that only one of them is shadowed, allowing a comparison between a shadowed and a direct signal. For the South Pole combined with Kamioka this probability is almost 75%, for the South Pole combined with Pyhasalmi it is almost 90%. One particular scenario consists of a large-volume scintillator detector located in Pyhasalmi to measure the geo-neutrino flux in a continental location and another such detector in Hawaii to measure it in an oceanic location. The probability that only one of them is shadowed exceeds 50% whereas the probability that at least one is shadowed is about 80%. We provide an online tool to calculate different shadowing probabilities for the one- and two-detector cases.

We analyze the neutrino events from the supernova (SN) 1987A detected by the Kamiokande II (KII) and Irvine-Michigan-Brookhaven (IMB) experiments. For the time-integrated flux we assume a quasi-thermal spectrum of the form $(E/E_0)^\alpha e^{-(\alpha+1)E/E_0}$ where $\alpha$ plays the role of a spectral index. This simple representation not only allows one to fit the total energy $E_{\rm tot}$ emitted in $\bar\nu_e$ and the average energy $<E_{\bar\nu_e}>$, but also accommodates a wide range of shapes, notably anti-pinched spectra that are broader than a thermal distribution. We find that the pile-up of low-energy events near threshold in KII forces the best-fit value for $\alpha$ to the lowest value of any assumed prior range. This applies to the KII events alone as well as to a common analysis of the two data sets. The preference of the data for an ``unphysical'' spectral shape implies that one can extract meaningful values for $<E_{\bar\nu_e}>$ and $E_{\rm tot}$ only if one fixes a prior value for $\alpha$. The tension between the KII and IMB data sets and theoretical expectations for $<E_{\bar\nu_e}>$ is not resolved by an anti-pinched spectrum.

Axion-photon conversion induced by intergalactic magnetic fields has been proposed as an explanation for the dimming of distant supernovae of type Ia SNe Ia) without cosmic acceleration. The effect depends on the intergalactic electron density $n_e$ as well as the $B$-field strength and domain size. We show that for $n_e\alt 10^{-9} {\rm cm}^{-3}$ the same mechanism would cause excessive spectral distortion of the cosmic microwave background (CMB). This small-$n_e$ parameter region had been left open by the most restrictive previous constraints based on the dispersion of quasar (QSO) spectra. The combination of CMB and QSO limits suggests that the photon-axion conversion mechanism can only play a subleading role for SN Ia dimming. A combined analysis of all the observables affected by the photon-axion oscillations would be required to give a final verdict on the viability of this model.

Observations of the cosmological large-scale structure provide well-established neutrino mass limits. We extend this argument to thermal relic axions. We calculate the axion thermal freeze-out temperature and thus their cosmological abundance on the basis of their interaction with pions. For hadronic axions we find a new mass limit $m_a<1.05$ eV (95% CL), corresponding to a limit on the axion decay constant of $f_a>5.7\times 10^6$ GeV. For other models this constraint is significantly weakened only if the axion-pion coupling is strongly suppressed. For comparison we note that the same approach leads to $\sum m_\nu<0.65$ eV (95% CL) for neutrinos.

Several current projects aim at building a large water-Cherenkov detector, with a fiducial volume about 20 times larger than in the current Super-Kamiokande experiment. These projects include the Underground nucleon decay and Neutrino Observatory (UNO) in the Henderson Mine (Colorado), the Hyper-Kamiokande (HK) detector in the Tochibora Mine (Japan), and the MEgaton class PHYSics (MEMPHYS) detector in the Frejus site (Europe). We study the physics potential of a reference next-generation detector (0.4 Mton of fiducial mass) in providing information on supernova neutrino flavor transitions with unprecedented statistics. After discussing the ingredients of our calculations, we compute neutrino event rates from inverse beta decay ($\bar\nu_e p\to e^+ n $), elastic scattering on electrons, and scattering on oxygen, with emphasis on their time spectra, which may encode combined information on neutrino oscillation parameters and on supernova forward (and possibly reverse) shock waves. In particular, we show that an appropriate ratio of low-to-high energy events can faithfully monitor the time evolution of the neutrino crossing probability along the shock-wave profile. We also discuss some background issues related to the detection of supernova relic neutrinos, with and without the addition of gadolinium.

If neutrinos have mass, they can also decay. Decay lifetimes of cosmological interest can be probed, in principle, through the detection of the redshifted, diffuse neutrino flux produced by all past supernovae--the so-called supernova relic neutrino (SRN) flux. In this work, we solve the SRN kinetic equations in the general case of three-generation flavor transitions followed by invisible (nonradiative) two-body decays. We then use the general solution to calculate observable SRN spectra in some representative decay scenarios. It is shown that, in the presence of decay, the SRN event rate can basically span the whole range below the current experimental upper bound--a range accessible to future experimental projects. Radiative SRN decays are also briefly discussed.

It has recently been realized that supernova neutrino signals may be affected by shock propagation over a time interval of a few seconds after bounce. In the standard three-neutrino oscillation scenario, such effects crucially depend on the neutrino level crossing probability P_H in the 1-3 sector. By using a simplified parametrization of the time-dependent supernova radial density profile, we explicitly show that simple analytical expressions for P_H accurately reproduce the phase-averaged results of numerical calculations in the relevant parameter space. Such expressions are then used to study the structure of P_H as a function of energy and time, with particular attention to cases involving multiple crossing along the shock profile. Illustrative applications are given in terms of positron spectra generated by supernova electron antineutrinos through inverse beta decay.