We present the results of the most complete ever scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine learning package. This is the first such study to separate out low-mass isotopes Be secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input galdef-files with these new parameters will be included in an upcoming public GALPROP update.

A full energy and flavor-dependent analysis of the three-year high-energy IceCube neutrino events is presented. By means of multidimensional fits, we derive the current preferred values of the highenergy neutrino flavor ratios, the normalization and spectral index of the astrophysical fluxes, and the expected atmospheric background events, including a prompt component. A crucial assumption resides on the choice of the energy interval used for the analyses, which significantly biases the results. When restricting ourselves to the ∼30 TeV−3 PeV energy range, which contains all the observed IceCube events, we find that the inclusion of the spectral information improves the fit to the canonical flavor composition at Earth, (1 : 1 : 1)⊕, with respect to a single-energy bin analysis. Increasing both the minimum and the maximum deposited energies has dramatic effects on the reconstructed flavor ratios as well as on the spectral index. Imposing a higher threshold of 60 TeV yields a slightly harder spectrum by allowing a larger muon neutrino component, since above this energy most atmospheric tracklike events are effectively removed. Extending the high-energy cutoff to fully cover the Glashow resonance region leads to a softer spectrum and a preference for tau neutrino dominance, as none of the expected electron antineutrino induced showers have been observed so far. The lack of showers at energies above 2 PeV may point to a broken power-law neutrino spectrum. Future data may confirm or falsify whether or not the recently discovered high-energy neutrino fluxes and the long-standing detected cosmic rays have a common origin.PACS numbers: 95.85. Ry, 14.60.Pq, 95.55.Vj, 29.40.Ka

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

We compute the bounds on the dark matter (DM) annihilation cross section using the most recent Cosmic Microwave Background measurements from WMAP9, SPT'11 and ACT'10. We consider DM with mass in the MeV-TeV range annihilating 100% into either an e + e − or a µ + µ − pair. We consider a realistic energy deposition model, which includes the dependence on the redshift, DM mass and annihilation channel. We exclude the canonical thermal relic abundance cross section ( σv = 3×10 −26 cm 3 s −1 ) for DM masses below 30 GeV and 15 GeV for the e + e − and µ + µ − channels, respectively. A priori, DM annihilating in halos could also modify the reionization history of the Universe at late times. We implement a realistic halo model taken from results of state-of-the-art N-body simulations and consider a mixed reionization mechanism, consisting on reionization from DM as well as from first stars. We find that the constraints on DM annihilation remain unchanged, even when large uncertainties on the halo model parameters are considered.

We use the latest cosmic microwave background (CMB) observations to provide updated constraints on the dark matter lifetime as well as on p-wave suppressed annihilation cross sections in the 1 MeV to 1 TeV mass range. In contrast to scenarios with an s-wave dominated annihilation cross section, which mainly affect the CMB close to the last scattering surface, signatures associated with these scenarios essentially appear at low redshifts (z 50) when structure began to form, and thus manifest at lower multipoles in the CMB power spectrum. We use data from Planck, WMAP9, SPT and ACT, as well as Lyman-α measurements of the matter temperature at z ∼ 4 to set a 95% confidence level lower bound on the dark matter lifetime of ∼ 4 × 10 25 s for mχ = 100 MeV. This bound becomes lower by an order of magnitude at mχ = 1 TeV due to inefficient energy deposition into the intergalactic medium. We also show that structure formation can enhance the effect of p-wave suppressed annihilation cross sections by many orders of magnitude with respect to the background cosmological rate, although even with this enhancement, CMB constraints are not yet strong enough to reach the thermal relic value of the cross section.

The IceCube experiment has recently reported the observation of 28 high-energy (>30  TeV) neutrino events, separated into 21 showers and 7 muon tracks, consistent with an extraterrestrial origin. In this Letter, we compute the compatibility of such an observation with possible combinations of neutrino flavors with relative proportion (αe:αμ∶ατ)⊕. Although the 7∶21 track-to-shower ratio is naively favored for the canonical (1∶1∶1)⊕ at Earth, this is not true once the atmospheric muon and neutrino backgrounds are properly accounted for. We find that, for an astrophysical neutrino E(-2) energy spectrum, (1∶1∶1)⊕ at Earth is disfavored at 81% C.L. If this proportion does not change, 6 more years of data would be needed to exclude (1∶1∶1)⊕ at Earth at 3σ C.L. Indeed, with the recently released 3-yr data, that flavor composition is excluded at 92% C.L. The best fit is obtained for (1∶0∶0)⊕ at Earth, which cannot be achieved from any flavor ratio at sources with averaged oscillations during propagation. If confirmed, this result would suggest either a misunderstanding of the expected background events or a misidentification of tracks as showers, or even more compellingly, some exotic physics which deviates from the standard scenario.

We show that the high-energy cosmic neutrinos seen by the IceCube Neutrino Observatory can be used to probe interactions between neutrinos and the dark sector that cannot be reached by current cosmological methods. The origin of the observed neutrinos is still unknown, and their arrival directions are compatible with an isotropic distribution. This observation, together with dedicated studies of Galactic plane correlations, suggests a predominantly extragalactic origin. Interactions between this isotropic extragalactic flux and the dense dark matter (DM) bulge of the Milky Way would thus lead to an observable imprint on the distribution, which would be seen by IceCube as (i) slightly suppressed fluxes at energies below a PeV and (ii) a deficit of events in the direction of the Galactic center. We perform an extended unbinned likelihood analysis using the four-year high-energy starting event data set to constrain the strength of DMneutrino interactions for two model classes. We find that, in spite of low statistics, IceCube can probe regions of the parameter space inaccessible to current cosmological methods.

After four years of data taking, the IceCube neutrino telescope has detected 54 high-energy starting events (HESE, or contained-vertex events) with deposited energies above 20 TeV. They represent the first ever detection of high-energy extraterrestrial neutrinos and therefore, the first step in neutrino astronomy. In order to study the energy, flavor and isotropy of the astrophysical neutrino flux arriving at Earth, we perform different analyses of two different deposited energy intervals, [10 TeV − 10 PeV] and [60 TeV − 10 PeV]. We first consider an isotropic unbroken power-law spectrum and constrain its shape, normalization and flavor composition. Our results are in agreement with the preliminary IceCube results, although we obtain a slightly softer spectrum. We also find that current data are not sensitive to a possible neutrino-antineutrino asymmetry in the astrophysical flux. Then, we show that although a two-component power-law model leads to a slightly better fit, it does not represent a significant improvement with respect to a single powerlaw flux. Finally, we analyze the possible existence of a North-South asymmetry, hinted at by the combination of the HESE sample with the through-going muon data. If only using HESE data, the scarce statistics from the northern hemisphere does not allow us to reach any conclusive answer, which indicates that the HESE sample alone is not driving the potential North-South asymmetry.