Primordial black holes (PBH), produced through a variety of processes in the early universe, could fill galactic halos accounting for a fraction or the totality of the dark matter. In particular, PBH with sub-stellar masses could be captured by stars, start to swallow their material, and eventually turn them into BH, hence originating a peculiar new type of stellar catastrophic event. Here we revisit this process in the most favorable case of PBH capture by neutron stars. We detail a number of novel features in the capture phase, during the settling within the star and mass growth of the PBH, and illustrate some phenomenological consequences. In particular, we point out that in the subsonic regime the PBH drag takes the form of a Bondi accretion. As a result, the onset of the final transmutation of the NS into a stellar sized BH is expected with the PBH seed in slight off-center position. We also compute the gravitational wave energy-loss and signals associated to different phases of the PBH-stellar interaction. In particular, the emission associated to the accretion phase is periodic with a few kHz frequency; in the rare case of a nearby Galactic event and for light PBH, it could constitute a warning of the forthcoming transmutation.

Unveiling the origin of the coalescing binaries detected via gravitational waves (GW) is challenging, notably if no multi-wavelength counterpart is detected. One important diagnostic tool is the coalescing binary distribution with respect to the large scale structures (LSS) of the universe, which we quantify via the cross-correlation of galaxy catalogs with GW ones. By using both existing and forthcoming galaxy catalogs and using realistic Monte Carlo simulations of GW events, we find that the cross-correlation signal should be marginally detectable in a 10-year data taking of advanced LIGO-Virgo detectors at design sensitivity, at least for binary neutron star mergers. The expected addition of KAGRA and LIGO-India to the GW detector network would allow for a firmer detection of this signal and, in combination with future cosmological surveys, would also permit the detection of cross-correlation for coalescing black holes. Such a measurement may unveil, for instance, a primordial origin of coalescing black holes. To attain this goal, we find that it is crucial to adopt a tomographic approach and to reach a sufficiently accurate localization of GW events. The depth of forthcoming surveys will be fully exploited by third generation GW detectors such as the Einstein Telescope or the Cosmic Explorer, which will allow one to perform precision studies of the coalescing black hole LSS distribution and attain rather advanced model discrimination capabilities.

Even if massive ($10\,M_\odot \lesssim M \lesssim 10^4 M_\odot$) primordial black holes (PBHs) can only account for a small fraction of the dark matter (DM) in the universe, they may still be responsible for a sizable fraction of the coalescence events measured by LIGO/Virgo, and/or act as progenitors of the supermassive black holes (SMBHs) observed already at high redshift ($z\gtrsim 6$). In presence of a dominant, non-PBH DM component, the bounds set by CMB via an altered ionization history are modified. We revisit the cosmological accretion of a DM halo around PBHs via toy models and dedicated numerical simulations, deriving updated CMB bounds which also take into account the last Planck data release. We prove that these constraints dominate over other constraints available in the literature at masses $M\gtrsim 20-50\,M_\odot$ (depending on uncertainty in accretion physics), reaching the level $f_{\rm PBH}<3\times 10^{-9}$ around $M\sim 10^{4}\,M_\odot$. These tight bounds are nonetheless consistent with the hypothesis of a primordial origin of the SMBH massive seeds.

The vertical diffusive halo size of the Galaxy, $L$, is a key parameter for dark matter indirect searches. It can be better determined thanks to recent AMS-02 data. We set constraints on $L$ from Be/B and $^{10}$Be/Be data, and we performed a consistency check with positron data. We detail the dependence of Be/B and $^{10}$Be/Be on $L$ and forecast on which energy range better data would be helpful for future $L$ improvements. We used USINE v3.5 for the propagation of nuclei, and $e^+$ were calculated with the pinching method of Boudaud et al. (2017). The current AMS-02 Be/B ($\sim3\%$ precision) and ACE-CRIS $^{10}$Be/Be ($\sim 10\%$ precision) data bring similar and consistent constraints on $L$. The AMS-02 Be/B data alone constrain $L=5^{+3}_{-2}$~kpc at a 68\% confidence level (spanning different benchmark transport configurations), a range for which most models do not overproduce positrons. Future experiments need to deliver percent-level accuracy on $^{10}$Be/$^9$Be anywhere below 10 GV to further constrain $L$. Forthcoming AMS-02, HELIX, and PAMELA $^{10}$Be/$^9$Be results will further test and possibly tighten the limits derived here. Elemental ratios involving radioactive species with different lifetimes (e.g. Al/Mg and Cl/Ar) are also awaited to provide complementary and robuster constraints.

The AMS-02 experiment has ushered cosmic-ray physics into precision era. In a companion paper, we designed an improved method to calibrate propagation models on B/C data. Here we provide a robust prediction of the $\bar{p}$ flux, accounting for several sources of uncertainties and their correlations. Combined with a correlation matrix for the $\bar{p}$ data, we show that the latter are consistent with a secondary origin. This paper presents key elements relevant to the dark matter search in this channel, notably by pointing out the inherent difficulties in achieving predictions at the percent-level precision.

We present an updated analysis of the gamma-ray flux from the directions of classical dwarf spheroidal galaxies, deriving new constraints on WIMP dark matter (DM) annihilation using a decade of Fermi-LAT data. Among the major novelties, we infer the dwarfs' J-factors by including new observations without imposing any a priori parametric profile for the DM distribution. While statistically compatible with results obtained from more conventional parameterisations, this procedure reduces the theoretical bias imposed on the data. Furthermore, we retain the full data-driven shape of the J-factors' empirical probability distributions when setting limits on DM, without imposing log-normality as is typically done. In conjunction with the data-driven J-factors, we use a new method for estimating the probability distribution function of the astrophysical background at the dwarf position, fully profiling over background uncertainties. We show that, for most classical dwarfs, the background systematic uncertainty dominates over the uncertainty on their J-factors.

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.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

This article aims at establishing new benchmark scenarios for Galactic cosmic-ray propagation in the GV-TV rigidity range, based on fits to the AMS-02 B/C data with the USINE v3.5 propagation code. We employ a new fitting procedure, cautiously taking into account data systematic error correlations in different rigidity bins and considering Solar modulation potential and leading nuclear cross-section as nuisance parameters. We delineate specific low, intermediate, and high-rigidity ranges that can be related to both features in the data and peculiar microphysics mechanisms resulting in spectral breaks. We single out a scenario which yields excellent fits to the data and includes all the presumably relevant complexity, the BIG model. This model has two limiting regimes: (i) the SLIM model, a minimal diffusion-only setup, and (ii) the QUAINT model, a convection-reacceleration model where transport is tuned by non-relativistic effects. All models lead to robust predictions in the high-energy regime ($\gtrsim10$GV), i.e. independent of the propagation scenario: at $1\sigma$, the diffusion slope $\delta$ is $[0.43-0.53]$, whereas $K_{10}$, the diffusion coefficient at 10GV, is $[0.26-0.36]$kpc$^2$Myr$^{-1}$; we confirm the robustness of the high-energy break, with a typical value $\Delta_h\sim 0.2$. We also find a hint for a similar (reversed) feature at low rigidity around the B/C peak ($\sim 4$GV) which might be related to some effective damping scale in the magnetic turbulence.

The variability of the high-frequency peaked BL Lac object PKS 2155-304 is studied using almost 10 years of optical, X-ray and gamma-rays data. Publicly available data have been gathered and analyzed with the aim to characterize the variability and to search for log-normality or periodic behavior. The optical and X-ray range follow a log-normal process; a hint for a periodicity of about 700 days is found in optical and in the high energy (100 MeV <E<300 GeV) range. A one zone, time-dependent, synchrotron self-Compton model is successfully used to reproduce the evolution with energy of the variability and the tentatively reported periodicity.

If the gamma-ray excess towards the inner Galaxy (GCE) detected in Fermi-LAT data is due to millisecond pulsars (MSPs), one expects an associated gravitational wave (GW) signal, whose intensity exceeds the disk MSP population emission by an order of magnitude. We compute the expected GW counterpart of the bulge MSP population based on fits of the GCE, and estimate the sensitivity reach of current and future terrestrial GW detectors. The bounds on the average population ellipticity $\varepsilon$ are competitive with the existing ones derived by LIGO/Virgo towards known MSPs. With a 10-year data taking in current LIGO/Virgo configuration, one would detect a signal at the level $\varepsilon\simeq 10^{-7}$, while $\varepsilon\simeq 10^{-8}$ would be attainable with a similar data taking period with a third generation GW detector. This sensitivity should be sufficient for crucial diagnostics on the GCE interpretation in terms of MSPs.

One of the seemingly strongest constraints on the fraction of dark matter in the form of primordial black holes (PBH) of ${\cal O}$(10)$\,M_\odot$ relies on the merger rate inferred from the binary BH merger events detected by LIGO/Virgo. The robustness of these bounds depends however on the accuracy with which the formation of PBH binaries in the early Universe can be described. We revisit the standard estimate of the merger rate, focusing on a couple of key ingredients: the spatial distribution of nearest neighbours and the initial clustering of PBHs associated to a given primordial power spectrum. Overall, we confirm the robustness of the results presented in the literature in the case of a narrow mass function (which constrain the PBH fraction of dark matter to be $f_{\rm PBH}\lesssim 0.001-0.01$). The initial clustering of PBHs might have an effect tightening the current constraint, but only for very broad mass functions, corresponding to wide bumps in the primordial power spectra extending at least over a couple of decades in $k$-space.

Theoretical predictions for the cosmic antiproton spectrum currently fall short of the corresponding experimental level of accuracy. Among the main sources of uncertainty are the antiproton production cross sections in cosmic ray inelastic interactions. We analyse existing data on antiproton production in $pp$ scattering, including for the first time the measurements performed by the NA49 Collaboration. We compute the antiproton spectrum finding that in the energy range where data are available (antiproton energies of about 4-550 GeV) different approaches lead to almost equivalent results, with an uncertainty of 10-20\%. Extrapolations outside this region lead to different estimates, with the uncertainties reaching the 50\% level around $1$ TeV, degrading the diagnostic power of the antiproton channel at those energies. We also comment on the uncertainties in the antiproton production source term coming from nuclei heavier than protons and from neutrons produced in $pp$ scatterings, and point out the need for dedicated experimental campaigns for all processes involving antiproton production in collisions of light nuclei.

Dwarf galaxies represent a powerful probe of annihilating dark matter particle models, with gamma-ray data setting some of the best bounds available. A major issue in improving over existing constraints consists in the limited knowledge of the astrophysical background (mostly diffuse photons, but also unresolved sources). Perhaps more worrisome, several approaches in the literature suffer of the difficulty of assessing the systematic error due to background mis-modelling. Here we propose a data-driven method to estimate the background at the dwarf position and its uncertainty, relying on an appropriate use of the whole-sky data, via an optimisation procedure of the interpolation weights. While this article is mostly methodologically oriented, we also report the bounds based on latest Fermi-LAT data and updated information for J-factors for both isolated and stacked dwarfs. Our results are very competitive with the Fermi-LAT ones, while being derived with a more general and flexible method. We discuss the impact of profiling over the J-factor as well as over the background probability distribution function, with the latter resulting for instance crucial in drawing conclusions of compatibility with DM interpretations of the so-called Galactic Centre Excess.

A new measurement of a spatially extended gamma-ray signal from the center of the Andromeda galaxy (M31) has been recently published by the Fermi-LAT collaboration, reporting that the emission broadly resembles the so-called Galactic center excess (GCE) of the Milky Way (MW). At the same time, evidence is accumulating on a millisecond pulsar (MSPs) origin for the GCE. These elements prompt us to compare the mentioned observations with what is, perhaps, the simplest model for an MSP population, solely obtained by rescaling of the MSP luminosity function determined in the local MW disk via the respective stellar mass of the systems. It is remarkable that without free fitting parameters, this model can account for both the energetics and the morphology of the GCE within uncertainties. For M31, the estimated luminosity due to primordial MSPs is expected to contribute only about a quarter of the detected emission, although a dominant contribution cannot be excluded given the large uncertainties. If correct, the model predicts that the M31 disk emission due to MSP is not far below the present upper bound. We also discuss a few refinements of this simple model. In particular, we use the correlation between globular cluster gamma-ray luminosity and stellar encounter rate to gauge the dynamical MSP formation in the bulge. This component is expected to contribute to the GCE only at a level $\lesssim 5\%$, but it may be of some importance in explaining the signal's morphology in the inner region of the Galaxy. We also comment on some effects which may lead to violations of the simple scaling used, on alternative models, and on future perspectives for improved diagnostics.

We analyze in details the standard Primordial Nucleosynthesis scenario. In particular we discuss the key theoretical issues which are involved in a detailed prediction of light nuclide abundances, as the weak reaction rates, neutrino decoupling and nuclear rate modeling. We also perform a new analysis of available data on the main nuclear processes entering the nucleosynthesis reaction network, with particular stress on their uncertainties as well as on their role in determining the corresponding uncertainties on light nuclide theoretical estimates. The current status of theoretical versus experimental results for 2H, 3He, 4He and 7Li is then discussed using the determination of the baryon density as obtained from Cosmic Microwave Background anisotropies.

Using cosmic-ray boron to carbon ratio (B/C) data recently released by the AMS-02 experiment, we find tantalizing indications ({\it decisive evidence}, in Bayesian terms) in favor of a diffusive origin for the broken power-law spectra found in protons ($p$) and helium nuclei (He). The result is robust with respect to currently estimated uncertainties in the cross sections, and in the presence of a small component of primary boron, expected because of spallation at the acceleration site. Reduced errors at high energy as well as further cosmic ray nuclei data (as absolute spectra of C, N, O, Li, Be) may definitively confirm this scenario.

Stellar-mass Primordial Black Holes (PBH) have been recently reconsidered as a Dark Matter (DM) candidate, after the aLIGO discovery of several binary BH mergers with masses of tens of $M_\odot$. Matter accretion on such massive objects leads to the emission of high-energy photons, capable of altering the ionization and thermal history of the universe. This in turn affects the statistical properties of the cosmic microwave background (CMB) anisotropies. Previous analyses have assumed spherical accretion. We argue that this approximation likely breaks down and that an accretion disk should form in the dark ages. Using the most up-to-date tools to compute the energy deposition in the medium, we derive constraints on the fraction of DM in PBH. Provided that disks form early on, even under conservative assumptions for accretion, these constraints exclude a monochromatic distribution of PBH with masses above $\sim 2\, M_\odot$ as the dominant form of DM. The bound on the median PBH mass gets more stringent if a broad, log-normal mass function is considered. A deepened understanding of non-linear clustering properties and BH accretion disk physics would permit an improved treatment and possibly lead to more stringent constraints.

If the electroweak sector of the standard model is described by classically conformal dynamics, the early Universe evolution can be substantially altered. It is already known that---contrarily to the standard model case---a first order electroweak phase transition may occur. Here we show that, depending on the model parameters, a dramatically different scenario may happen: A first-order, six massless quark QCD phase transition occurs first, which then triggers the electroweak symmetry breaking. We derive the necessary conditions for this dynamics to occur, using the specific example of the classically conformal B-L model. In particular, relatively light weakly coupled particles are predicted, with implications for collider searches. This scenario is also potentially rich in cosmological consequences, such as renewed possibilities for electroweak baryogenesis, altered dark matter production, and gravitational wave production, as we briefly comment upon.

We compute cosmic microwave background (CMB) anisotropy constraints on exotic forms of energy injection in electromagnetic (e.m.) channels over a large range of timescales. These constraints are very powerful around or just after recombination, although CMB keeps some sensitivity e.g. to decaying species with lifetimes as long as $10^{25}\,$s. We review here complementary with CMB spectral distortions and primordial nucleosynthesis bounds, which dominate at earlier timescales. For the first time, we describe the effects of the e.m. energy injection on the CMB power spectra as a function of the injection epoch, using the lifetime of a decaying particle as proxy. We identify a suitable on-the-spot approximation. Our results are of interest not only for early universe relics constituting (a fraction of) the dark matter, but also for other exotic injection of e.m. radiation. For illustration, we apply our formalism to: i) Primordial black holes of mass $\in [10^{13.5},10^{16.8}]$ g, showing that the constraints are comparable to the ones obtained from gamma-ray background studies and even dominate below $\sim 10^{14}$g. ii) To a peculiar mass-mixing range in the sterile neutrino parameter space, complementary to other astrophysical and laboratory probes. iii) Finally, we provide a first estimate of the room for improvement left for forthcoming 21 cm experiments, comparing it with the reach of proposed CMB spectral distortion (PiXiE) and CMB angular power spectrum (CORE) missions. We show that the best and most realistic opportunity to look for this signal (or to improve over current constraints) in the 21 cm probe is to focus on the Cosmic Dawn epoch, $15\lesssim z\lesssim30$, where the qualitatively unambiguous signature of a spectrum in emission can be expected for models that evade all current constraints.

In the new precision era for cosmic ray astrophysics, theoretical predictions cannot content themselves with average trends, but need to correctly take into account intrinsic uncertainties. The space-time discreteness of the cosmic ray sources, together with a substantial ignorance of their precise epochs and locations (with the possible exception of the most recent and close ones) play an important role in this sense. We elaborate a statistical theory to deal with this problem, relating the composite probability P({\Psi}) to obtain a flux {\Psi} at the Earth and the single-source probability p({\psi}) to contribute with a flux {\psi}. The main difficulty arises from the fact that p({\psi}) is a fat tail distribution, characterized by power-law or broken power-law behavior up to very large fluxes for which central limit theorem does not hold, and leading to well-known stable laws as opposed Gaussian. We find that relatively simple recipes provide a satisfactory description of the probability P({\Psi}). We also find that a naive Gaussian fit to simulation results would underestimate the probability of very large fluxes, that is, several times above the average, while overestimating the probability of relatively milder excursions. At large energies, large flux fluctuations are prevented by causal considerations, while at low energies a partial knowledge on the recent and nearby population of sources plays an important role. A few proposal have been discussed in the literature to account for spectral breaks recently reported in cosmic ray data in terms of local contributions. We apply our newly developed theory to assess their probabilities, finding that they are relatively small.

A hardening of the proton and helium fluxes is observed above a few hundreds of GeV/nuc. The distribution of local sources of primary cosmic rays has been suggested as a potential solution to this puzzling behavior. Some authors even claim that a single source is responsible for the observed anomalies. But how probable these explanations are? To answer that question, our current description of cosmic ray Galactic propagation needs to be replaced by the Myriad model. In the former approach, sources of protons and helium nuclei are treated as a jelly continuously spread over space and time. A more accurate description is provided by the Myriad model where sources are considered as point-like events. This leads to a probabilistic derivation of the fluxes of primary species, and opens the possibility that larger-than-average values may be observed at the Earth. For a long time though, a major obstacle has been the infinite variance associated to the probability distribution function which the fluxes follow. Several suggestions have been made to cure this problem but none is entirely satisfactory. We go a step further here and solve the infinite variance problem of the Myriad model by making use of the generalized central limit theorem. We find that primary fluxes are distributed according to a stable law with heavy tail, well-known to financial analysts. The probability that the proton and helium anomalies are sourced by local SNR can then be calculated. The p-values associated to the CREAM measurements turn out to be small, unless somewhat unrealistic propagation parameters are assumed.

Current data (LHC direct searches, Higgs mass, dark matter-related bounds) severely affect the constrained minimal SUSY standard model (CMSSM) with neutralinos as dark matter candidates. But the evidence for neutrino masses coming from oscillations requires extending the SM with at least right-handed neutrinos with a Dirac mass term. In turn, this implies extending the CMSSM with right-handed sneutrino superpartners, a scenario we dub $\tilde\nu$CMSSM. These additional states constitute alternative dark matter candidates of the superWIMP type, produced via the decay of the long-lived next-to-lightest SUSY particle (NLSP). Here we consider the interesting and likely case where the NLSP is a $\tilde{\tau}$: despite the modest extension with respect to the CMSSM this scenario has the distinctive signatures of heavy, stable charged particles. After taking into account the role played by neutrino mass bounds and the specific cosmological bounds from the big bang nucleosynthesis in selecting the viable parameter space, we discuss the excellent discovery prospects for this model at the future runs of the LHC. We show that it is possible to probe $\tilde{\tau}$ masses up to 600 GeV at the 14 TeV LHC with $\mathcal{L} = 1100$ fb$^{-1}$ when one considers a pair production of staus with two or more hard jets through all SUSY processes. We also show the complementary discovery prospects from a direct $\tilde{\tau}$ pair production, as well as at the new experiment MoEDAL.

We consider a cosmological model in which a fraction $f$ of the Dark Matter (DM) is allowed to decay in an invisible relativistic component, and compute the resulting constraints on both the decay width (or inverse lifetime) $\Gamma$ and $f$ from purely gravitational arguments. We report a full derivation of the Boltzmann hierarchy, correcting a mistake in previous literature, and compute the impact of the decay --as a function of the lifetime-- on the CMB and matter power spectra. From CMB only, we obtain that no more than 3.8 % of the DM could have decayed in the time between recombination and today (all bounds quoted at 95 % CL). We also comment on the important application of this bound to the case where primordial black holes constitute DM, a scenario notoriously difficult to constrain. For lifetimes longer than the age of the Universe, the bounds can be cast as $f\Gamma < 6.3\times10^{-3}$ Gyr$^{-1}$. For the first time, we also checked that degeneracies with massive neutrinos are broken when information from the large scale structure is used. Even secondary effects like CMB lensing suffice to this purpose. Decaying DM models have been invoked to solve a possible tension between low redshift astronomical measurements of $\sigma_8$ and $\Omega_{\rm m}$ and the ones inferred by Planck. We reassess this claim finding that with the most recent BAO, HST and $\sigma_8$ data extracted from the CFHT survey, the tension is only slightly reduced despite the two additional free parameters, loosening the bound to $f\Gamma < 15.9\times10^{-3}$ Gyr$^{-1}$. The bound however improves to $f\Gamma < 5.9\times10^{-3}$ Gyr$^{-1}$ if only data consistent with the CMB are included. This highlights the importance of establishing whether the tension is due to real physical effects or unaccounted systematics, for settling the reach of achievable constraints on decaying DM.

Blazars represent the most abundant class of high-energy extragalactic $\gamma$-ray sources. The subset of blazars known as BL Lac objects is on average closer to Earth and characterized by harder spectra at high energy than the whole sample. The fraction of BL Lacs that is too dim to be detected and resolved by current $\gamma$-ray telescopes is therefore expected to contribute to the high-energy isotropic diffuse $\gamma$-ray background (IGRB). The IGRB has been recently measured over a wide energy range by the Large Area Telescope (LAT) on board the Gamma-ray Space Telescope ({\it Fermi}). We present a new prediction of the diffuse $\gamma$-ray flux due to the unresolved BL Lac blazar population. The model is built upon the spectral energy distribution and the luminosity function derived from the fraction of BL Lacs detected (and spectrally characterized) in the $\gamma$-ray energy range. We focus our attention on the ${\cal O}(100)$ GeV energy range, predicting the emission up to the TeV scale and taking into account the absorption on the extragalactic background light. In order to better shape the BL Lac spectral energy distribution, we combine the {\it Fermi}-LAT data with Imaging Atmospheric Cerenkov Telescopes measurements of the most energetic sources. Our analysis is carried on separately for low- and intermediate-synchrotron-peaked BL Lacs on one hand, and high-synchrotron-peaked BL Lacs on the other one: we find in fact statistically different features for the two. The diffuse emission from the sum of both BL Lac classes increases from about 10$\%$ of the measured IGRB at 100 MeV to $\sim$100$\%$ of the data level at 100 GeV. At energies greater than 100 GeV, our predictions naturally explain the IGRB data, accommodating their softening with increasing energy. Uncertainties are estimated to be within of a factor of two of the best-fit flux up to 500 GeV.

It is well known that annihilations in the homogeneous fluid of dark matter (DM) can leave imprints in the cosmic microwave background (CMB) anisotropy power spectrum. However, the relevance of DM annihilations in halos for cosmological observables is still subject to debate, with previous works reaching different conclusions on this point. Also, all previous studies used a single type of parameterization for the astrophysical reionization, and included no astrophysical source for the heating of the intergalactic medium. In this work, we revisit these problems. When standard approaches are adopted, we find that the ionization fraction does exhibit a very particular (and potentially constraining) pattern, but the currently measurable optical depth to reionization is left almost unchanged: In agreement with the most of the previous literature, for plausible halo models we find that the modification of the signal with respect to the one coming from annihilations in the smooth background is tiny, below cosmic variance within currently allowed parameter space. However, if different and probably more realistic treatments of the astrophysical sources of reionization and heating are adopted, a more pronounced effect of the DM annihilation in halos is possible. We thus conclude that within currently adopted baseline models the impact of the virialised DM structures cannot be uncovered by CMB power spectra measurements, but a larger impact is possible if peculiar models are invoked for the redshift evolution of the DM annihilation signal or different assumptions are made for the astrophysical contributions. A better understanding (both theoretical and observational) of the reionization and temperature history of the universe, notably via the 21 cm signal, seems the most promising way for using halo formation as a tool in DM searches, improving over the sensitivity of current cosmological probes.

The very high energy Galactic $\gamma$-ray sky is partially opaque in the ($0.1-10$) PeV energy range. In the light of the recently detected high energy neutrino flux by IceCube, a comparable very high energy $\gamma$-ray flux is expected in any scenario with a sizable Galactic contribution to the neutrino flux. Here we elaborate on the peculiar energy and anisotropy features imposed upon these very high energy $\gamma$-rays by the absorption on the cosmic microwave background photons and Galactic interstellar light. As a notable application of our considerations, we study the prospects of probing the PeV-scale decaying DM scenario, proposed as a possible source of IceCube neutrinos, by extensive air shower (EAS) cosmic ray experiments. In particular, we show that anisotropy measurements at EAS experiments are already sensitive to $\tau_{\rm DM}\sim {\cal O}(10^{27})$~s and future measurements, using better gamma/hadron separation, can improve the limit significantly.

Direct techniques for cosmic ray observations have reached an unprecedented level of precision, unveiling fine-details of the energy spectra. I will introduce the evidence for new spectral features which has been accumulated by new experiments over the past few years, and review the main ideas invoked in the theoretical explanations of the revealed spectral breaks and elemental spectra non-universality. I will also briefly comment on the complementary situation of antimatter observations.

Using the updated proton and helium fluxes just released by the AMS-02 experiment we reevaluate the secondary astrophysical antiproton to proton ratio and its uncertainties, and compare it with the ratio preliminarly reported by AMS-02. We find no unambiguous evidence for a significant excess with respect to expectations. Yet, some preference for a flatter energy dependence of the diffusion coefficient starts to emerge. Also, we provide a first assessment of the room left for exotic components such as Galactic Dark Matter annihilation or decay, deriving new stringent constraints.

PAMELA and, more recently, AMS-02, are ushering us into a new era of greatly reduced statistical uncertainties in experimental measurements of cosmic-ray fluxes. In particular, new determinations of traditional diagnostic tools such as the boron-to-carbon ratio (B/C) are expected to significantly reduce errors on cosmic-ray diffusion parameters, with important implications for astroparticle physics, ranging from inferring primary source spectra to indirect dark matter searches. It is timely to stress, however, that the conclusions obtained crucially depend on the framework in which the data are interpreted as well as from some nuclear input parameters. We aim at assessing the theoretical uncertainties affecting the outcome, with models as simple as possible while still retaining the key dependencies. We compare different semi-analytical, two-zone model descriptions of cosmic-ray transport in the Galaxy. We test for the effect of a primary source contamination in the boron flux by parametrically altering its flux, as well as for nuclear cross section uncertainties. Our study on preliminary results from AMS-02 suggests that, differently for instance from the leptonic case, realistic modelling of the geometry of the Galaxy and of the source distribution are of minor importance to correctly reproduce B/C data at high energies and thus, to a large extent, for the extraction of diffusion parameters. The Ansatz on the lack of primary injection of boron represents the most serious bias, and requires multi-messenger studies to be addressed. If this uncertainty could be lifted, nuclear uncertainties would still represent a serious concern, which degrade the systematic error on the inferred parameters to the 20% level, or three times the estimated experimental sensitivity. In order to reduce this, a new nuclear cross section measurement campaign is probably required.

Context: Features in the spectra of primary cosmic rays (CRs) provide invaluable information on the propagation of these particles in the Galaxy. In the rigidity region around a few hundred GV, such features have been measured in the proton and helium spectra by the PAMELA experiment and later confirmed with a higher significance by AMS-02. We investigate the implications of these datasets for the scenario in which CRs propagate under the action of self-generated waves. Aims: We show that the recent data on the spectrum of protons and helium nuclei as collected with AMS-02 and Voyager are in very good agreement with the predictions of a model in which the transport of Galactic CRs is regulated by self-generated waves. We also study the implications of the scenario for the boron-to-carbon ratio: although a good overall agreement is found, at high energy we find marginal support for a (quasi) energy independent contribution to the grammage, that we argue may come from the sources themselves Results: A break in the spectra of all nuclei is found at rigidity of a few hundred GV, as a result of a transition from self-generated waves to pre-existing waves with a Kolmogorov power spectrum. Neither the slope of the diffusion coefficient, nor its normalisation are free parameters. Moreover, at rigidities below a few GV, CRs are predicted to be advected with the self-generated waves at the local Alfv\'en speed. This effect, predicted in our previous work, provides an excellent fit to the Voyager data on the proton and helium spectra at low energies, providing additional support to the model.

In Poulin and Serpico [Phys. Rev. Lett. 114, 091101 (2015)] we have recently argued that when the energy of a photon injected in the primordial plasma falls below the pair-production threshold, the universality of the non-thermal photon spectrum from the standard theory of electromagnetic cascades onto a photon background breaks down. We showed that this could reopen or widen the parameter space for an exotic solution to the 'lithium problem'. Here we discuss another application, namely the impact that this has on non-thermal big bang nucleosynthesis constraints from 4He, 3He and 2H, using the parametric example of monochromatic photon injection of different energies. Typically, we find tighter bounds than those existing in the literature, up to more than one order of magnitude. As a consequence of the non-universality of the spectrum, the energy-dependence of the photodissociation cross-sections is important. We also compare the constraints obtained with current level and future reach of cosmic microwave background spectral distortion bounds.

The standard theory of electromagnetic cascades onto a photon background predicts a quasi-universal shape for the resulting non-thermal photon spectrum. This has been applied to very disparate fields, including non-thermal big bang nucleosynthesis (BBN). However, once the energy of the injected photons falls below the pair-production threshold the spectral shape is very different, a fact that has been overlooked in past literature. This loophole may have important phenomenological consequences, since it generically alters the BBN bounds on non-thermal relics: for instance it allows to re-open the possibility of purely electromagnetic solutions to the so-called "cosmological lithium problem", which were thought to be excluded by other cosmological constraints. We show this with a proof-of-principle example and a simple particle physics model, compared with previous literature.

Several groups of authors have analyzed Fermi LAT data in a region around the Galactic Center finding an unaccounted gamma-ray excess over diffuse backgrounds in the GeV energy range. It has been argued that it is difficult or even impossible to explain this diffuse emission by the leading astrophysical candidates - millisecond pulsars (MSPs). Here we provide a new estimate of the contribution to the excess by a population of yet unresolved MSP located in the bulge of the Milky Way. We simulate this population with the GALPLOT package by adopting a parametric approach, with the range of free parameters gauged on the MSP characteristics reported by the second pulsar catalogue (2PC). We find that the conclusions strongly depend on the details of the MSP luminosity function (in particular, its high luminosity end) and other explicit or tacit assumptions on the MSP statistical properties, which we discuss. Notably, for the first time we study the importance of the possible secondary emission of the MSPs in the Galactic Center, {i.e. the emission via inverse Compton losses of electrons injected in the interstellar medium. Differently from a majority of other authors, we find that within current uncertainties a large if not dominant contribution of MSPs to the excess cannot be excluded. We also show that the sensitivities of future instruments or possibly already of the latest LAT data analysis (Pass 8) provide good perspectives to test this scenario by resolving a significant number of MSPs.

Several groups have recently claimed evidence for an unaccounted gamma-ray excess over {the} diffuse backgrounds at few GeV in {the} Fermi-LAT data in a region around the Galactic Center, consistent with a dark matter annihilation origin. We demonstrate that the main spectral and angular features of this excess can be reproduced if they are mostly due to inverse Compton emission from high-energy electrons injected in a burst event of ~10^52 - 10^53erg roughly O(10^6) years ago. We consider this example as a proof of principle that time-dependent phenomena need to be understood and accounted for - together with detailed diffuse foregrounds and unaccounted "steady state" astrophysical sources - before any robust inference can be made about dark matter signals at the Galactic Center. In addition, we point out that the timescale suggested by our study, which controls both the energy cutoff and the angular extension of the signal, intriguingly matches (together with the energy budget) what is indirectly inferred by other evidences suggesting a very active Galactic Center in the past, for instance related to intense star formation and accretion phenomena.

In the light of the new IceCube data on the (yet unidentified) astrophysical neutrino flux in the PeV and sub-PeV range, we present an update on the status of decaying dark matter interpretation of the events. In particular, we develop further the angular distribution analysis and discuss the perspectives for diagnostics. By performing various statistical tests (maximum likelihood, Kolmogorov-Smirnov and Anderson-Darling tests) we conclude that currently the data show a mild preference (below the two sigma level) for the angular distribution expected from dark matter decay vs. the isotropic distribution foreseen for a conventional astrophysical flux of extragalactic origin. Also, we briefly develop some general considerations on heavy dark matter model building and on the compatibility of the expected energy spectrum of decay products with the IceCube data, as well as with existing bounds from gamma-rays. Alternatively, assuming that the IceCube data originate from conventional astrophysical sources, we derive bounds on both decaying and annihilating dark matter for various final states. The lower limits on heavy dark matter lifetime improve by up to an order of magnitude with respect to existing constraints, definitively making these events---even if astrophysical in origin---an important tool for astroparticle physics studies.

Dwarf spheroidal galaxies of the Local Group are close satellites of the Milky Way characterized by a large mass-to-light ratio and are not expected to be the site of non-thermal high-energy gamma-ray emission or intense star formation. Therefore they are amongst the most promising candidates for indirect dark matter searches. During the last years the High Energy Stereoscopic System (H.E.S.S.) of imaging atmospheric Cherenkov telescopes observed five of these dwarf galaxies for more than 140 hours in total, searching for TeV gamma-ray emission from annihilation of dark matter particles. The new results of the deep exposure of the Sagittarius dwarf spheroidal galaxy, the first observations of the Coma Berenices and Fornax dwarves and the re-analysis of two more dwarf spheroidal galaxies already published by the H.E.S.S. Collaboration, Carina and Sculptor, are presented. In the absence of a significant signal new constraints on the annihilation cross-section applicable to Weakly Interacting Massive Particles (WIMPs) are derived by combining the observations of the five dwarf galaxies. The combined exclusion limit depends on the WIMP mass and the best constraint is reached at 1-2 TeV masses with a cross-section upper bound of ~3.9x10-24 cm^3 s-1 at a 95% confidence level.

We study the impact of the existence of an eV-mass scale sterile neutrino---with parameters in the ballpark of what is required to fit the laboratory anomalies---on the early time profile of the electron neutrino and antineutrino fluxes associated to a core-collapse supernova (SN). In particular, we focus on the universal feature of neutronization burst expected in the first tens of ms of the signal: provided that a detector with sufficient sensitivity is available, it is well-known that in the 3 neutrino framework the detection of the neutronization burst in neutrino channel would signal inverted mass hierarchy. This conclusion is dramatically altered in the presence of a sterile neutrino: we study here both analytically and numerically the region in parameter space where this characteristic signal disappears, mimicking normal hierarchy expectations. Conversely, the detection of a peak consistent with expectations for inverted mass hierarchy would exclude the existence of a sterile state over a much wider parameter space than what required by laboratory anomalies fits, or even probed by detectors coming on-line in the near future. Additionally, we show the peculiar alteration in the energy-time double differential flux, with a delayed peak appearing for kinematical reasons, which might offer a remarkable signature in case of favorable parameters and for a high statistics detection of a Galactic SN. We also comment on additional potentially interesting effects in the electron antineutrino channel, if more than one angle in the active-sterile sector is non-vanishing. As an ancillary result that we derived in the technical resolution of the equations, in an appendix we report the Cayley-Hamilton formalism for the evolution of a four neutrino system in matter, generalizing existing results in the literature.

Dark Matter is essential for structure formation in the late Universe so it must be stable on cosmological time scales. But how stable exactly? Only assuming decays into relativistic particles, we report an otherwise model independent bound on the lifetime of Dark Matter using current cosmological data. Since these decays affect only the low-$\ell$ multipoles of the CMB, the Dark Matter lifetime is expected to correlate with the tensor-to-scalar ratio $r$ as well as curvature $\Omega_k$. We consider two models, including $r$ and $r+\Omega_k$ respectively, versus data from Planck, WMAP, WiggleZ and Baryon Acoustic Oscillations, with or without the BICEP2 data (if interpreted in terms of primordial gravitational waves). This results in a lower bound on the lifetime of CDM given by 160Gyr (without BICEP2) or 200Gyr (with BICEP2) at 95% confidence level.

We re-evaluate the extragalactic gamma-ray flux prediction from dark matter annihilation in the approach of integrating over the nonlinear matter power spectrum, extrapolated to the free-streaming scale. We provide an estimate of the uncertainty based entirely on available N-body simulation results and minimal theoretical assumptions. We illustrate how an improvement in the simulation resolution, exemplified by the comparison between the Millennium and Millennium II simulations, affects our estimate of the flux uncertainty and we provide a "best guess" value for the flux multiplier, based on the assumption of stable clustering for the dark matter perturbations described as a collision-less fluid. We achieve results comparable to traditional Halo Model calculations, but with a much simpler procedure and a more general approach, as it relies only on one, directly measurable quantity. In addition we discuss the extension of our calculation to include baryonic effects as modeled in hydrodynamical cosmological simulations and other possible sources of uncertainty that would in turn affect indirect dark matter signals. Upper limit on the integrated power spectrum from supernovae lensing magnification are also derived and compared with theoretical expectations.

Recent observations by IceCube, notably two PeV cascades accompanied by events at energies ~ (30-400) TeV, are clearly in excess over atmospheric background fluxes and beg for an astroparticle physics explanation. Although some models of astrophysical accelerators can account for the observations within current statistics, intriguing features in the energy and possibly angular distributions of the events make worth exploring alternatives. Here, we entertain the possibility of interpreting the data with a few PeV mass scale decaying Dark Matter, with lifetime of the order of 10^27 s. We discuss generic signatures of this scenario, including its unique energy spectrum distortion with respect to the benchmark $E_\nu^{-2}$ expectation for astrophysical sources, as well as peculiar anisotropies. A direct comparison with the data show a good match with the above-mentioned features. We further discuss possible future checks of this scenario.

We discuss the often-neglected role of bremsstrahlung processes on the interstellar gas in computing indirect signatures of Dark Matter (DM) annihilation in the Galaxy, particularly for light DM candidates in the phenomenologically interesting O(10) GeV mass range. Especially from directions close to the Galactic Plane, the expected gamma-ray spectrum is altered via two effects: directly, by the photons emitted in the bremsstrahlung process on the interstellar gas by energetic electrons which are among the DM annihilation byproducts; indirectly, by the modification of the same electron spectrum, due to the additional energy loss process in the diffusion-loss equation (e.g. the resulting inverse Compton emission is altered). We quantify the importance of the bremsstrahlung emission in the GeV energy range, showing that it is the dominant component of the gamma-ray spectrum for some cases. We also find that, in regions in which bremsstrahlung dominates energy losses, the related gamma-ray emission is only moderately sensitive to possible large variations in the gas density. Still, we stress that, for computing precise spectra in the (sub-)GeV range, it is important to obtain a reliable description of the inner Galaxy gas distribution as well as to compute self-consistently the gamma emission and the solution to the diffusion-loss equation. For example, these are crucial issues to quantify and interpret meaningfully gamma-ray map `residuals' in terms of (light) DM annihilations.

New bounds on decaying Dark Matter are derived from the gamma-ray measurements of (i) the isotropic residual (extragalactic) background by Fermi and (ii) the Fornax galaxy cluster by H.E.S.S. We find that those from (i) are among the most stringent constraints currently available, for a large range of dark matter masses and a variety of decay modes, excluding half-lives up to about 10^26 to few 10^27 seconds. In particular, they rule out the interpretation in terms of decaying dark matter of the e+/- spectral features in PAMELA, Fermi and H.E.S.S., unless very conservative choices are adopted. We also discuss future prospects for CTA bounds from Fornax which, contrary to the present H.E.S.S. constraints of (ii), may allow for an interesting improvement and may become better than those from the current or future extragalactic Fermi data.

We consider the propagation and energy losses of protons and anti-protons produced by dark matter annihilation at redshifts 100<z<~2000. In the case of dark matter annihilations into quarks, gluons and weak gauge bosons, protons and anti-protons carry about 20% of the energy injected into e^\pm and \gamma's, but their interactions are normally neglected when deriving cosmic microwave background bounds from altered recombination histories. Here, we follow numerically the energy-loss history of typical protons/antiprotons in the cosmological medium. We show that about half of their energy is channeled into photons and e^\pm, and we present a simple prescription to estimate the corresponding strengthening of the cosmic microwave background bounds on the dark matter annihilation cross section.

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).

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.

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 derive new bounds on decaying Dark Matter from the gamma ray measurements of (i) the isotropic residual (extragalactic) background by Fermi and (ii) the Fornax galaxy cluster by HESS. We find that those from (i) are among the most stringent constraints currently available, for a large range of DM masses and a variety of decay modes, excluding half-lives up to ~10^26 to few 10^27 seconds. In particular, they rule out the interpretation in terms of decaying DM of the e^\pm spectral features in Pamela, Fermi and HESS, unless very conservative choices are adopted. We also discuss future prospects for CTA bounds from Fornax which, contrary to the present HESS constraints of (ii), may allow for an interesting improvement and may become better than those from the current or future extragalactic Fermi data.

We show that the complex shape of the cosmic ray (CR) spectrum, as recently measured by PAMELA and inferred from Fermi-LAT gamma-ray observations of molecular clouds in the Gould belt, can be naturally understood in terms of basic plasma astrophysics phenomena. A break from a harder to a softer spectrum at blue rigidity R\simeq 10 GV follows from a transition from transport dominated by advection of particles with Alfven waves to a regime where diffusion in the turbulence generated by the same CRs is dominant. A second break at R\simeq 200 GV happens when the diffusive propagation is no longer determined by the self-generated turbulence, but rather by the cascading of externally generated turbulence (for instance due to supernova (SN) bubbles) from large spatial scales to smaller scales where CRs can resonate. Implications of this scenario for the cosmic ray spectrum, grammage and anisotropy are discussed.

Over the last three years, several satellite and balloon observatories have suggested intriguing features in the cosmic ray lepton spectra. Most notably, the PAMELA satellite has suggested an "anomalous" rise with energy of the cosmic ray positron fraction. In this article, we summarize the global picture emerging from the data and recapitulate the main features of different types of explanations proposed. The perspectives in testing different scenarios as well as inferring some astrophysical diagnostics from current/near future experiments are also discussed.

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.

In the curvaton scenario, residual isocurvature perturbations can be imprinted in the cosmic neutrino component after the decay of the curvaton field, implying in turn a non-zero chemical potential in the neutrino distribution. We study the constraints that future experiments like Planck, SPIDER or CMBPol will be able to put on the amplitude of isocurvature perturbations in the neutrino component. We express our results in terms of the square root \gamma of the non-adiabaticity parameter \alpha and of the extra relativistic degrees of freedom \Delta N_eff. Assuming a fiducial model with purely adiabatic fluctuations, we find that Planck (SPIDER) will be able to put the following upper limits at the 1sigma level: \gamma < 5.3x10^-3 (1.2x10^-2) and \Delta N_eff < 0.16 (0.40) . CMBPol will further improve these constraints to \gamma < 1.5x10^-3 and \Delta N_eff < 0.043. Finally, we recast these bounds in terms of the background neutrino degeneracy parameter \xi\ and the corresponding perturbation amplitude \sigma_\xi, and compare with the bounds on \xi\ that can be derived from Big Bang Nucleosynthesis.

We revisit the computation of the extragalactic gamma-ray signal from cosmological dark matter annihilations. The prediction of this signal is notoriously model dependent, due to different descriptions of the clumpiness of the dark matter distribution at small scales, responsible for an enhancement with respect to the smoothly distributed case. We show how a direct computation of this "flux multiplier" in terms of the nonlinear power spectrum offers a conceptually simpler approach and may ease some problems, such as the extrapolation issue. In fact very simple analytical recipes to construct the power spectrum yield results similar to the popular Halo Model expectations, with a straightforward alternative estimate of errors. For this specific application, one also obviates to the need of identifying (often literature-dependent) concepts entering the Halo Model, to compare different simulations.

We derive here a robust bound on the effective number of neutrinos from constraints on primordial nucleosynthesis yields of deuterium and helium. In particular, our results are based on very weak assumptions on the astrophysical determination of the helium abundance, namely that the minimum effect of stellar processing is to keep constant (rather than increase, as expected) the helium content of a low-metallicity gas. Using the results of a recent analysis of extragalactic HII regions as upper limit, we find that Delta Neff<= 1 at 95 % C.L., quite independently of measurements on the baryon density from cosmic microwave background anisotropy data and of the neutron lifetime input. In our approach, we also find that primordial nucleosynthesis alone has no significant preference for an effective number of neutrinos larger than the standard value. The ~2 sigma hint sometimes reported in the literature is thus driven by CMB data alone and/or is the result of a questionable regression protocol to infer a measurement of primordial helium abundance.

Recent data from CREAM seem to confirm early suggestions that primary cosmic ray (CR) spectra at few TeV/nucleon are harder than in the 10-100 GeV range. Also, helium and heavier nuclei spectra appear systematically harder than the proton fluxes at corresponding energies. We note here that if the measurements reflect intrinsic features in the interstellar fluxes (as opposed to local effects) appreciable modifications are expected in the sub-TeV range for the secondary yields, such as antiprotons and diffuse gamma-rays. Presently, the ignorance on the origin of the features represents a systematic error in the extraction of astrophysical parameters as well as for background estimates for indirect dark matter searches. We find that the spectral modifications are appreciable above 100 GeV, and can be responsible for ~30% effects for antiprotons at energies close to 1 TeV or for gamma's at energies close to 300 GeV, compared to currently considered predictions based on simple extrapolation of input fluxes from low energy data. Alternatively, if the feature originates from local sources, uncorrelated spectral changes might show up in antiproton and high-energy gamma-rays, with the latter ones likely dependent from the line-of-sight.

A convincing identification of dark matter (DM) particles can probably be achieved only through a combined analysis of different detections strategies, which provides an effective way of removing degeneracies in the parameter space of DM models. In practice, however, this program is made complicated by the fact that different strategies depend on different physical quantities, or on the same quantities but in a different way, making the treatment of systematic errors rather tricky. We discuss here the uncertainties on the recoil rate in direct detection experiments and on the muon rate induced by neutrinos from dark matter annihilations in the Sun, and we show that, contrarily to the local DM density or overall cross section scale, irreducible astrophysical uncertainties affect the two rates in a different fashion, therefore limiting our ability to reconstruct the parameters of the dark matter particle. By varying within their respective errors astrophysical parameters such as the escape velocity and the velocity dispersion of dark matter particles, we show that the uncertainty on the relative strength of the neutrino and direct-detection signal is as large as a factor of two for typical values of the parameters, but can be even larger in some circumstances.

We consider the diffuse gamma ray data from FERMI first year observations and compare them to the gamma ray fluxes predicted by Dark Matter annihilation or decay (both from prompt emission and from Inverse Compton Scattering), for different observation regions of the sky and a range of Dark Matter masses, annihilation/decay channels and Dark Matter galactic profiles. We find that the data exclude large regions of the Dark Matter parameter space not constrained otherwise and discuss possible directions for future improvements. Also, we further constrain Dark Matter interpretations of the e+e- PAMELA/FERMI spectral anomalies, both for the annihilating and the decaying Dark Matter case: under very conservative assumptions, only models producing dominantly mu+mu- and assuming a cored Dark Matter galactic profile can fit the lepton data with masses around 2 TeV.

Interpretations of indirect searches for dark matter (DM) require theoretical predictions for the annihilation or decay rates of DM into stable particles of the standard model. These predictions include usually only final states accessible as lowest order tree-level processes, with electromagnetic bremsstrahlung and the loop-suppressed two gamma-ray line as exceptions. We show that this restriction may lead to severely biased results for DM tailored to produce only leptons in final states and with mass in the TeV range. For such models, unavoidable electroweak bremsstrahlung of Z and W-bosons has a significant influence both on the branching ratio and the spectral shape of the final state particles. We work out the consequences for two situations: Firstly, the idealized case where DM annihilates at tree level with 100% branching ratio into neutrinos. For a given cross section, this leads eventually to "minimal yields" of photons, electrons, positrons and antiprotons. Secondly, the case where the only allowed two-body final states are electrons. The latter case is typical of models aimed at fitting cosmic ray e^- and e^+ data. We find that the multimessenger signatures of such models can be significantly modified with respect to results presented in the literature.

Annihilation of cosmologically distributed dark matter is predicted to produce a potentially observable flux of high energy photons. Neglecting the contribution from local structure, this signal is predicted to be virtually uniform on the sky and, in order to be identified, must compete with various extragalactic backgrounds. We focus here on unresolved blazars and discuss several techniques for separating the dark matter signal from this background. First, the spectral shapes of the signal and background are expected to differ, a feature which can be exploited with the Fisher Matrix formalism. Second, in any given angular pixel, the number of photons from blazars is drawn from a distribution which is far from Poisson. A knowledge of this distribution enhances one's ability to extract the dark matter signal, while ignorance of it can lead to the introduction of a large systematic error.

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.

A recently proposed model (arXiv:0903.2794) explains the rise in energy of the positron fraction measured by the PAMELA satellite in terms of hadronic production of positrons in aged supernova remnants, and acceleration therein. Here we present a preliminary calculation of the anti-proton flux produced by the same mechanism. While the model is consistent with present data, a rise of the antiproton to proton ratio is predicted at high energy, which strikingly distinguishes this scenario from other astrophysical explanations of the positron fraction (like pulsars). We briefly discuss important implications for Dark Matter searches via antimatter.

We perform an autocorrelation study of the Auger data with the aim to constrain the number density n_s of ultrahigh energy cosmic ray (UHECR) sources, estimating at the same time the effect on $n_s$ of the systematic energy scale uncertainty and of the distribution of UHECR. The use of global analysis has the advantage that no biases are introduced, either in n_s or in the related error bar, by the a priori choice of a single angular scale. The case of continuous, uniformly distributed sources is nominally disfavored at 99% C.L. and the fit improves if the sources follow the large-scale structure of matter in the universe. The best fit values obtained for the number density of proton sources are within a factor ~2 around n_s ~ 10^{-4}/Mpc^3 and depend mainly on the Auger energy calibration scale, with lower densities being preferred if the current scale is correct. The data show no significant small-scale clustering on scales smaller than a few degrees. This might be interpreted as a signature of magnetic smearing of comparable size, comparable with the indication of a ~3 degrees magnetic deflection coming from cross-correlation results. The effects on the above results of some approximations done is also discussed.

We present an up-to-date review of Big Bang Nucleosynthesis (BBN). We discuss the main improvements which have been achieved in the past two decades on the overall theoretical framework, summarize the impact of new experimental results on nuclear reaction rates, and critically re-examine the astrophysical determinations of light nuclei abundances. We report then on how BBN can be used as a powerful test of new physics, constraining a wide range of ideas and theoretical models of fundamental interactions beyond the standard model of strong and electroweak forces and Einstein's general relativity.

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.

Recent results from the PAMELA satellite indicate the presence of a large flux of positrons (relative to electrons) in the cosmic ray spectrum between approximately 10 and 100 GeV. As annihilating dark matter particles in many models are predicted to contribute to the cosmic ray positron spectrum in this energy range, a great deal of interest has resulted from this observation. Here, we consider pulsars (rapidly spinning, magnetized neutron stars) as an alternative source of this signal. After calculating the contribution to the cosmic ray positron and electron spectra from pulsars, we find that the spectrum observed by PAMELA could plausibly originate from such sources. In particular, a significant contribution is expected from the sum of all mature pulsars throughout the Milky Way, as well as from the most nearby mature pulsars (such as Geminga and B0656+14). The signal from nearby pulsars is expected to generate a small but significant dipole anisotropy in the cosmic ray electron spectrum, potentially providing a method by which the Fermi gamma-ray space telescope would be capable of discriminating between the pulsar and dark matter origins of the observed high energy positrons.

Based on general considerations rather than model-dependent fits to specific scenarios, we argue that an increase with energy of the positron fraction in cosmic rays, suggested by several experiments at E>~7 GeV, most likely requires a primary source of electron-positron pairs. We discuss the possible alternatives, and find none of them plausible on astrophysical or particle physics grounds. Further observational ways to test different scenarios are discussed.

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.

Although the emission of radiation from dark matter annihilation is expected to be maximized at the Galactic Center, geometric factors and the presence of point-like and diffuse backgrounds make the choice of the angular window size to optimize the chance of a signal detection a non-trivial problem. Contrarily to what is often assumed, we find that the best strategy is to focus on a window size around the Galactic Center of ~ 1 deg to >~ 30 deg, where the optimal size depends on the angular distribution of the signal and the backgrounds. Although our conclusions are general, we illustrate this point in the particular case of annihilation into two monochromatic photons in the phenomenologically most interesting range of energy 45 GeV < E < 80 GeV, which is of great interest for the GLAST satellite. We find for example that Dark Matter models with sufficiently strong line annihilation signals, like the Inert Doublet Model, may be detectable without or with very moderate boost factors.

We describe a program for computing the abundances of light elements produced during Big Bang Nucleosynthesis which is publicly available at http://parthenope.na.infn.it/. Starting from nuclear statistical equilibrium conditions the program solves the set of coupled ordinary differential equations, follows the departure from chemical equilibrium of nuclear species, and determines their asymptotic abundances as function of several input cosmological parameters as the baryon density, the number of effective neutrino, the value of cosmological constant and the neutrino chemical potential. The program requires commercial NAG library routines.

The arrival directions of ultrahigh energy cosmic rays (UHECRs) may show anisotropies on all scales, from just above the experimental angular resolution up to medium scales and dipole anisotropies. We find that a global comparison of the two-point auto-correlation function of the data with the one of catalogues of potential sources is a powerful diagnostic tool. In particular, this method is far less sensitive to unknown deflections in magnetic fields than cross-correlation studies while keeping a strong discrimination power among source candidates. We illustrate these advantages by considering ordinary galaxies, gamma ray bursts and active galactic nuclei as possible sources. Already the sparse publicly available data suggest that the sources of UHECRs may be a strongly clustered sub-sample of galaxies or of active galactic nuclei. We present forecasts for various cases of source distributions which can be checked soon by the Pierre Auger Observatory.

The GLAST satellite mission will study the gamma ray sky with considerably greater exposure than its predecessor EGRET. In addition, it will be capable of measuring the arrival directions of gamma rays with much greater precision. These features each significantly enhance GLAST's potential for identifying gamma rays produced in the annihilations of dark matter particles. The combined use of spectral and angular information, however, is essential if the full sensitivity of GLAST to dark matter is to be exploited. In this paper, we discuss the separation of dark matter annihilation products from astrophysical backgrounds, focusing on the Galactic Center region, and perform a forecast for such an analysis. We consider both point-like and diffuse astrophysical backgrounds and model them using a point-spread-function for GLAST. While the results of our study depend on the specific characteristics of the dark matter signal and astrophysical backgrounds, we find that in many scenarios it is possible to successfully identify dark matter annihilation radiation, even in the presence of significant astrophysical backgrounds.

Very high energy gamma-rays are expected to be absorbed by the extragalactic background light over cosmological distances via the process of electron-positron pair production. Recent observations of cosmologically distant gamma-ray emitters by ground based gamma-ray telescopes have, however, revealed a surprising degree of transparency of the universe to very high energy photons. One possible mechanism to explain this observation is the oscillation between photons and axion-like-particles (ALPs). Here we explore this possibility further, focusing on photon-ALP conversion in the magnetic fields in and around gamma-ray sources and in the magnetic field of the Milky Way, where some fraction of the ALP flux is converted back into photons. We show that this mechanism can be efficient in allowed regions of the ALP parameter space, as well as in typical configurations of the Galactic Magnetic Field. As case examples, we consider the spectrum observed from two HESS sources: 1ES1101-232 at redshift z=0.186 and H 2356-309 at z=0.165. We also discuss features of this scenario which could be used to distinguish it from standard or other exotic models.

We perform a new estimate of the high energy neutrinos expected from GRBs associated with the first generation of stars in light of new models and constraints on the epoch of reionization and a more detailed evaluation of the neutrino emission yields. We also compare the diffuse high energy neutrino background from Population III stars with the one from "ordinary stars" (Population II), as estimated consistently within the same cosmological and astrophysical assumptions. In disagreement with previous literature, we find that high energy neutrinos from Population III stars will not be observable with current or near future neutrino telescopes, falling below both IceCube sensitivity and atmospheric neutrino background under the most extreme assumptions for the GRB rate. This rules them out as a viable diagnostic tool for these still elusive metal-free stars.

We compare the clustering properties of the combined dataset of ultra-high energy cosmic rays events, reported by the AGASA, HiRes, Yakutsk and Sugar collaborations, with a catalogue of galaxies of the local universe (redshift z<~0.06). We find that the data reproduce particularly well the clustering properties of the nearby universe within z <~0.02. There is no statistically significant cross-correlation between data and structures, although intriguingly the nominal cross-correlation chance probability drops from ~50% to ~10% using the catalogue with a smaller horizon. Also, we discuss the impact on the robustness of the results of deflections in some galactic magnetic field models used in the literature. These results suggest a relevant role of magnetic fields (possibly extragalactic ones, too) and/or possibly some heavy nuclei fraction in the UHECRs. The importance of a confirmation of these hints by Auger data is emphasized.

We propose that axion-like particles (ALPs) with a two-photon vertex, consistent with all astrophysical and laboratory bounds, may lead to a detectable signature in the spectra of high-energy gamma ray sources. This occurs as a result of gamma rays being converted into ALPs in the magnetic fields of efficient astrophysical accelerators according to the "Hillas criterion", such as jets of active galactic nuclei or hot spots of radio galaxies. The discovery of such an effect is possible by GLAST in the 1-100 GeV range and by ground based gamma ray telescopes in the TeV range.

The energy--zenith angular event distribution in a neutrino telescope provides a unique tool to determine at the same time the neutrino-nucleon cross section at extreme kinematical regions, and the high energy neutrino flux. By using a simple parametrization for fluxes and cross sections, we present a sensitivity analysis for the case of a km^3 neutrino telescope. In particular, we consider the specific case of an under-water Mediterranean telescope placed at the NEMO site, although most of our results also apply to an under-ice detector such as IceCube. We determine the sensitivity to departures from standard values of the cross sections above 1 PeV which can be probed independently from an a-priori knowledge of the normalization and energy dependence of the flux. We also stress that the capability to tag downgoing neutrino showers in the PeV range against the cosmic ray induced background of penetrating muons appears to be a crucial requirement to derive meaningful constraints on the cross section.

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.

The high energy neutrino detection by a km^3 Neutrino Telescope placed in the Mediterranean sea provides a unique tool to both determine the diffuse astrophysical neutrino flux and the neutrino-nucleon cross section in the extreme kinematical region, which could unveil the presence of new physics. Here is performed a brief analysis of possible NEMO site performances.

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.

Dark matter annihilating in our Galaxy's halo and elsewhere in the universe is expected to generate a diffuse flux of gamma rays, potentially observable with next generation satellite-based experiments, such as GLAST. In this article, we study the signatures of dark matter in the angular distribution of this radiation, in particular the deterministic ones. Pertaining to the extragalactic contribution, we discuss the effect of the motion of the solar system with respect to the cosmological rest frame, and anisotropies due to the structure of our local universe. For the gamma ray flux from dark matter in our own Galactic halo, we discuss the effects of the offset position of the solar system, the Compton-Getting effect, the asphericity of the Milky Way halo, and the signatures of nearby substructure. We explore the prospects for the detection of these features by the GLAST satellite and find that, if ~10% or more of the diffuse gamma ray background observed by EGRET is the result of dark matter annihilations, then GLAST should be sensitive to anisotropies at the sub-percent level. Such precision would be sufficient to detect many, if not all, of the signatures discussed in this paper.

We perform a reanalysis of the production of CNO elements in a standard Big Bang Nucleosynthesis scenario. The CNO yields in BBN are suppressed by the low density of the plasma, Coulomb barrier effects and the short time scales involved. Yet, the inclusion of nuclides and reactions traditionally disregarded may lead to an increase relevant enough to affect the pristine Population III stars. After a critical reanalysis and upgrade of the nuclear network our results show no major discrepancies with the ones obtained using a smaller nuclear network. The robustness of the standard predictions-the early generation of star developed in a metal-free environment-is confirmed.

Future cosmological data may be sensitive to the effects of a finite sum of neutrino masses even as small as ~0.06 eV, the lower limit guaranteed by neutrino oscillation experiments. We show that a cosmological detection of neutrino mass at that level would improve by many orders of magnitude the existing limits on neutrino lifetime, and as a consequence on neutrino secret interactions with (quasi-)massless particles as in majoron models. On the other hand, neutrino decay may provide a way-out to explain a discrepancy <~ 0.1 eV between cosmic neutrino bounds and Lab data.

If the diffuse extragalactic gamma ray emission traces the large scale structures of the universe, peculiar anisotropy patterns are expected in the gamma ray sky. In particular, because of the cutoff distance introduced by the absorption of 0.1-10 TeV photons on the infrared/optical background, prominent correlations with the local structures within a range of few hundreds Mpc should be present. We provide detailed predictions of the signal based on the PSCz map of the local universe. We also use mock N-body catalogues complemented with the halo model of structures to study some statistical features of the expected signatures. The results are largely independent from cosmological details, and depend mostly on the index of correlation (or bias) of the sources with respect to the large scale distribution of galaxies. For instance, the predicted signal in the case of a quadratic correlation (as it may happen for a dark matter annihilation contribution to the diffuse gamma flux) differs substantially from a linear correlation case, providing a complementary tool to unveil the nature of the sources of the diffuse gamma ray emission. The chances of the present and future space and ground based observatories to measure these features are discussed.

We perform a study of the ultra high energy neutrino detection performances of a km^3 Neutrino Telescope sitting at the three proposed sites for ANTARES, NEMO and NESTOR in the Mediterranean sea. We focus on the effect of the underwater surface profile on the total amount of yearly expected tau and mu crossing the fiducial volume in the limit of full detection efficiency and energy resolution. We also emphasize the possible enhancement of matter effect by a suitable choice of the geometry of the Telescope.

The result of the recent publication [M. Kachelriess and D. V. Semikoz Astropart. Phys. 26, 10 (2006)] of a broad maximum around 25 degrees in the two-point autocorrelation function of ultra-high energy cosmic ray arrival directions has been intriguingly interpreted as the first imprint of the large scale structures (LSS) of baryonic matter in the near universe. We analyze this suggestion in light of the clustering properties expected from the PSCz astronomical catalogue of LSS. The chance probability of the signal is consistent within 2 sigma with the predictions based on the catalogue. No evidence for a significant cross-correlation of the observed events with known overdensities in the LSS is found, which may be due to the role of the galactic and extragalactic magnetic fields, and is however consistent with the limited statistics. The larger statistics to be collected by the Pierre Auger Observatory is needed to answer definitely the question.

We reconsider constraints on Brans-Dicke theories arising from the requirement of successful Big Bang Nucleosynthesis. Such constraints typically arise by imposing that the universe be radiation-dominated at early times, and therefore restricting the contribution that a Brans-Dicke scalar could make to the energy budget of the universe. However, in this paper we show how the dynamics of the Brans-Dicke scalar itself can mimic a portion of the radiation contribution, thereby allowing successful nucleosynthesis with a sizable contribution to the total cosmic energy density. This possibility significantly relaxes the existing bounds on Brans-Dicke fields, and opens the door to new possibilities for early universe cosmology. The necessary fine tunings required by such an arrangement are identified and discussed.

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).

We consider the decoupling of neutrinos in the early Universe in presence of non-standard neutral current neutrino-electron interactions (NSI). We first discuss a semi-analytical approach to solve the relevant kinetic equations and then present the results of fully numerical and momentum-dependent calculations, including flavor neutrino oscillations. We present our results in terms of both the effective number of neutrino species (N_eff) and the impact on the abundance of He-4 produced during Big Bang Nucleosynthesis. We find that, for NSI parameters within the ranges allowed by present laboratory data, non-standard neutrino-electron interactions do not essentially modify the density of relic neutrinos nor the bounds on neutrino properties from cosmological observables, such as their mass. Nonetheless, the presence of neutrino-electron NSI may enhance the entropy transfer from electron-positron pairs into neutrinos instead of photons, up to a value of N_eff=3.12. This is almost three times the correction to N_eff=3 that appears for standard weak interactions.

Deviations from isotropy have been a key tool to identify the origin and the primary type of cosmic rays at low energies. We suggest that the Compton-Getting effect can play a similar role at ultra-high energies: If at these energies the cosmic ray flux is dominated by sources at cosmological distances, then the movement of the Sun relative to the cosmic microwave background frame induces a dipole anisotropy at the 0.6% level. The energy dependence and the orientation of this anisotropy provide important information about the transition between galactic and extragalactic cosmic rays, the charge of the cosmic ray primaries, the galactic magnetic field and, at the highest energies, the energy-loss horizon of cosmic rays. A 3-sigma detection of this effect requires around 10^6 events in the considered energy range and is thus challenging but not impossible with present detectors. As a corollary we note that the Compton-Getting effect allows one also to constrain the fraction of the diffuse gamma-ray background emitted by sources at cosmological distance, with promising detection possibilities for the GLAST satellite.

We study the influence of the regular component of the Galactic magnetic field (GMF) on the arrival directions of ultra-high energy cosmic rays (UHECRs). We find that, if the angular resolution of current experiments has to be fully exploited, deflections in the GMF cannot be neglected even for E=10^20 eV protons, especially for trajectories along the Galactic plane or crossing the Galactic center region. On the other hand, the GMF could be used as a spectrograph to discriminate among different source models and/or primaries of UHECRs, if its structure would be known with sufficient precision. We compare several GMF models introduced in the literature and discuss for the example of the AGASA data set how the significance of small-scale clustering or correlations with given astrophysical sources are affected by the GMF. We point out that the non-uniform exposure to the extragalactic sky induced by the GMF should be taken into account estimating the significance of potential (auto-)correlation signals.

The standard cosmological model predicts the existence of a cosmic neutrino background with a present density of about 110 cm^{-3} per flavour, which affects big-bang nucleosynthesis, cosmic microwave background anisotropies, and the evolution of large scale structures. We report on a precision calculation of the cosmic neutrino background properties including the modification introduced by neutrino oscillations. The role of a possible neutrino-antineutrino asymmetry and the impact of non-standard neutrino-electron interactions on the relic neutrinos are also briefly discussed.

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.

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 discuss the possibility to probe leptonic mixing parameters at high-energy neutrino telescopes in a model-independent way, using astrophysical neutron and pion sources. In particular we show how the octant of the 2-3 mixing angle might be determined independently of prior knowledge of the source, even when current uncertainties on the other mixing parameters are included. We also argue that non-trivial neutrino oscillation effects should be taken into account when using high-energy flavor ratios for astrophysical diagnostics.

Current experiments collecting high statistics in ultra-high energy cosmic rays (UHECRs) are opening a new window on the universe. In this work we discuss a large scale structure model for the UHECR origin which evaluates the expected anisotropy in the UHECR arrival distribution starting from a given astronomical catalogue of the local universe. The model takes into account the main selection effects in the catalogue and the UHECR propagation effects. By applying this method to the IRAS PSCz catalogue, we derive the minimum statistics needed to significatively reject the hypothesis that UHECRs trace the baryonic distribution in the universe, in particular providing a forecast for the Auger experiment.

We review the current cosmological status of neutrinos, with particular emphasis on their effects on Big Bang Nucleosynthesis, Large Scale Structure of the universe and Cosmic Microwave Background Radiation measurements.

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.

We calculate and display the primordial light-element abundances as a function of a neutrino degeneracy parameter \xi common to all flavors. It is the only unknown parameter characterizing the thermal medium at the primordial nucleosynthesis epoch. The observed primordial helium abundance Y_p is the most sensitive cosmic ``leptometer.'' Adopting the conservative Y_p error analysis of Olive and Skillman implies -0.04 \alt \xi \alt 0.07 whereas the errors stated by Izotov and Thuan imply \xi=0.0245+-0.0092 (1 sigma). Improved determinations of the baryon abundance have no significant impact on this situation. A determination of Y_p that reliably distinguishes between a vanishing or nonvanishing \xi is a crucial test of the cosmological standard assumption that sphaleron effects equilibrate the cosmic lepton and baryon asymmetries.

In the early universe, neutrinos are slightly coupled when electron-positron pairs annihilate transferring their entropy to photons. This process originates non-thermal distortions on the neutrino spectra which depend on neutrino flavour, larger for nu_e than for nu_mu or nu_tau. We study the effect of three-neutrino flavour oscillations on the process of neutrino decoupling by solving the momentum-dependent kinetic equations for the neutrino spectra. We find that oscillations do not essentially modify the total change in the neutrino energy density, giving N_eff=3.046 in terms of the effective number of neutrinos, while the small effect over the production of primordial 4He is increased by O(20%), up to 2.1 x 10^{-4}. These results are stable within the presently favoured region of neutrino mixing parameters.

The observed excess of high-energy cosmic rays from the Galactic plane in the energy range \sim 10^18 eV may be naturally explained by neutron primaries generated in the photo-dissociation of heavy nuclei. In this scenario, neutrons with lower energy decay before reaching the Earth and produce a detectable flux in a 1 km^3 neutrino telescope. The initial flavor composition of these neutrinos, \phi(\bar\nu_e):\phi(\bar\nu_\mu):\phi(\bar\nu_\tau)=1:0:0, offers the opportunity to perform a combined \bar\nu_\mu/\bar\nu_\tau appearance and \bar\nu_e disappearance experiment. The observable ratio \phi(\bar\nu_\mu)/\phi(\bar\nu_e+\bar\nu_\tau) of fluxes arriving on Earth depends appreciably on the 13-mixing angle \theta_13 and the leptonic CP phase \delta_CP, opening thus a new experimental avenue to measure these two quantities.

The importance of the Big Bang Nucleosynthesis (BBN) as a unique tool for studying neutrino properties is discussed, and the recent steps towards a self-consistent and robust handling of the weak reaction decoupling from the thermal bath as well as of the neutrino reheating following the e+e- annihilation are summarized. We also emphasize the important role of the Cosmic Microwave Background (CMB) anisotropy in providing an accurate and independent determination of the baryon density parameter omegab. The BBN is presently a powerful parameter-free theory that can test the standard scenario of the neutrino decoupling in the early Universe. Moreover it can constrain new physics in the neutrino sector. The perspectives for improvements in the next years are outlined.

We study the expected diffuse cosmic neutrino flux produced by Population III (PopIII) stars during their nuclear burning phases as well as from their final stages of evolution (core collapse). Assuming a fraction f_III=10^(-3) of all baryons forms PopIII stars, our flux estimate is comparable to the diffuse neutrino flux produced by the ordinary stars and by the ordinary core-collapse supernovae in the universe, i.e. of order 1-10 cm^(-2) s^(-1). Due to the large cosmic redshift, however, the typical energies are in the MeV and sub-MeV range where the solar and geophysical neutrino fluxes are much larger. A direct detection of the diffuse cosmic flux is out of the question with presently known experimental techniques.

Nowadays, the Cosmic Microwave Background (CMB) anisotropies studies accurately determine the baryon fraction omega_b, showing an overall and striking agreement with previous determinations of omega_b obtained from Big Bang Nucleosynthesis (BBN). However, a deeper comparison of BBN predictions with the determinations of the primordial light nuclides abundances shows slight tensions, motivating an effort to further improve the accuracy of theoretical predictions, as well as to better evaluate systematics in both observations and nuclear reactions measurements. We present some results of an important step towards an increasing precision of BBN predictions, namely an updated and critical review of the nuclear network, and a new protocol to perform the nuclear data regression.

The annihilation of new dark matter candidates with masses $m_X$ in the MeV range may account for the galactic positrons that are required to explain the 511 keV $\gamma$-ray flux from the galactic bulge. We study the impact of MeV-mass thermal relic particles on the primordial synthesis of $^2$H, $^4$He, and $^7$Li. If the new particles are in thermal equilibrium with neutrinos during the nucleosynthesis epoch they increase the helium mass fraction for $m_X\alt 10$ MeV and are thus disfavored. If they couple primarily to the electromagnetic plasma they can have the opposite effect of lowering both helium and deuterium. For $m_X=4$--10 MeV they can even improve the overall agreement between the predicted and observed $^2$H and $^4$He abundances.

We report on the status of primordial nucleosynthesis in light of recent results on CMB anisotropies from WMAP experiment. Theoretical estimates for nuclei abundances, along with the corresponding uncertainties, are evaluated using a new numerical code, where all nuclear rates usually considered have been updated using the most recent available data. Moreover, additional processes, neglected in previous calculations, have been included. The combined analysis of CMB and primordial nucleosynthesis prediction for Deuterium gives an effective number of relativistic degrees of freedom in good agreement with the simplest scenario of three non degenerate neutrinos. Our findings seem to point out possible systematics affecting 4He mass fraction measurements, or the effect of exotic physics, like a slightly degenerate relic neutrino background.

The key Standard-Physics inputs of the Big Bang Nucleosynthesis (BBN) are the light nuclei reaction rates. Both the network and the nuclear rates have been recently reanalyzed and updated, and cosmological and New-Physics constraints (taking into account the WMAP Cosmic Microwave Background anisotropies measurement) obtained with a new code are presented.