DEAP-3600 is a single-phase liquid argon detector aiming to directly detect Weakly Interacting Massive Particles (WIMPs), located at SNOLAB (Sudbury, Canada). After analyzing data taken during the first year of operation, a null result was used to place an upper bound on the WIMP-nucleon spin-independent, isoscalar cross section. This study reinterprets this result within a Non-Relativistic Effective Field Theory framework, and further examines how various possible substructures in the local dark matter halo may affect these constraints. Such substructures are hinted at by kinematic structures in the local stellar distribution observed by the Gaia satellite and other recent astronomical surveys. These include the Gaia Sausage (or Enceladus), as well as a number of distinct streams identified in recent studies. Limits are presented for the coupling strength of the effective contact interaction operators $\mathcal{O}_1$, $\mathcal{O}_3$, $\mathcal{O}_5$, $\mathcal{O}_8$, and $\mathcal{O}_{11}$, considering isoscalar, isovector, and xenonphobic scenarios, as well as the specific operators corresponding to millicharge, magnetic dipole, electric dipole, and anapole interactions. The effects of halo substructures on each of these operators are explored as well, showing that the $\mathcal{O}_5$ and $\mathcal{O}_8$ operators are particularly sensitive to the velocity distribution, even at dark matter masses above 100 GeV/$c^2$.

We place updated constraints on dark matter annihilation into neutrinos, using the most recently available data from neutrino detectors. We consider Galactic and extragalactic signals of $s$, $p$, and $d$-wave annihilation processes directly into neutrino pairs, and furthermore provide projections for next-generation neutrino experiments. Our constraints span dark matter masses from MeV to ZeV, ranging from $\langle \sigma v \rangle \sim 5\times10^{-26}~{\rm cm}^3 {\rm s}^{-1}$ to $10^{-19}~{\rm cm}^3{\rm s}^{-1}$, respectively. When directional information is available, constraints are strongest. Experiments that report directional and energy information of their events provide much stronger constraints, outlining the importance of making such data public.

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

We 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 describe nuFATE Neutrino Fast Attenuation Through Earth, a very rapid method of accurately computing the attenuation of high-energy neutrinos during their passage through Earth to detectors such as IceCube, ANTARES or KM3Net, including production of secondary neutrinos from $\tau^\pm$ lepton decay. We then use this method to quantify the error on attenuation due to uncertainties in the isotropic neutrino spectrum, the composition of the Earth, and the parton distribution functions. We show that these can be as large as 20%, which can significantly impact reconstructed astrophysical neutrino parameters, as well as searches for new physics. An implementation of this algorithm is provided as a public code.

We investigate the production of circularly polarised X and gamma-ray signals in cosmic accelerators such as supernova remnants and AGN jets. Proton-proton and proton-photon collisions within these sites produce a charge asymmetry in the distribution of mesons and muons that eventually leads to a net circular polarisation signal as these particles decay radiatively. We find that the fraction of circular polarisation thus produced is at the level of $ 5 \times 10^{-4}$, regardless of the exact beam spectrum, as long as it is made predominantly of protons. While this fraction is very small, the detection of circular polarisation signals in conjunction with high-energy neutrinos would provide an unambiguous signature of the presence of high-energy protons in cosmic accelerators. In Supernovae shocks in particular, this would indicate the presence of relativistic protons hitting stationary protons and/or low-energy photons in the intergalactic or interstellar medium.

As hints of physics beyond the standard model become the major driving force behind future large-scale projects, it is increasingly important to consider all sources of evidence and constraints. Here, I illustrate the importance of considering the connection between particle physics, cosmology and astrophysics, mainly via two examples: 1) the sterile neutrino and its impact on cosmology, and 2) the 511 keV line from electron-positron annihilation in the galactic centre.

The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta.

Several next-generation experiments aim to make the first measurement of the neutrino flux from the Carbon-Nitrogen-Oxygen (CNO) solar fusion cycle. We calculate how much time these experiments will need to run for in order to measure this flux with enough precision to tell us the metal content of the Sun's core, and thereby help to solve the solar metallicity problem. For experiments looking at neutrino-electron scattering, we find that SNO+ will measure this CNO neutrino flux with enough precision after five years in its pure scintillator mode, provided its $^{210}$Bi background is measured to 1% accuracy. By comparison, a 100~ton liquid argon experiment such as Argo will take ten years in Gran Sasso lab, or five years in SNOLAB or Jinping. Borexino could obtain this precision in ten years, but this projection is very sensitive to background assumptions. For experiments looking at neutrino-nucleus scattering, the best prospects are obtained for low-threshold solid state detectors (employing either germanium or silicon). These would require new technologies to lower the experimental threshold close to detection of single electron-hole pairs, and exposures beyond those projected for next-generation dark matter detectors.

We use present cosmological observations and forecasts of future experiments to illustrate the power of large-scale structure (LSS) surveys in probing dark matter (DM) microphysics and unveiling potential deviations from the standard $\Lambda$CDM scenario. To quantify this statement, we focus on an extension of $\Lambda$CDM with DM-neutrino scattering, which leaves a distinctive imprint on the angular and matter power spectra. After finding that future CMB experiments (such as COrE+) will not significantly improve the constraints set by the Planck satellite, we show that the next generation of galaxy clustering surveys (such as DESI) could play a leading role in constraining alternative cosmologies and even have the potential to make a discovery. Typically we find that DESI would be an order of magnitude more sensitive to DM interactions than Planck, thus probing effects that until now have only been accessible via $N$-body simulations.

Dark matter with momentum- or velocity-dependent interactions with nuclei has shown significant promise for explaining the so-called Solar Abundance Problem, a longstanding discrepancy between solar spectroscopy and helioseismology. The best-fit models are all rather light, typically with masses in the range of 3-5 GeV. This is exactly the mass range where dark matter evaporation from the Sun can be important, but to date no detailed calculation of the evaporation of such models has been performed. Here we carry out this calculation, for the first time including arbitrary velocity- and momentum-dependent interactions, thermal effects, and a completely general treatment valid from the optically thin limit all the way through to the optically thick regime. We find that depending on the dark matter mass, interaction strength and type, the mass below which evaporation is relevant can vary from 1 to 4 GeV. This has the effect of weakening some of the better-fitting solutions to the Solar Abundance Problem, but also improving a number of others. As a by-product, we also provide an improved derivation of the capture rate that takes into account thermal and optical depth effects, allowing the standard result to be smoothly matched to the well-known saturation limit.

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, suggest 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 1) slightly suppressed fluxes at energies below a PeV and 2) 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 dataset to constrain the strength of DM-neutrino 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.

The study of anomalous electromagnetic emission in the sky is the basis of indirect searches for dark matter. It is also a powerful tool to constrain the radiative decay of active neutrinos. Until now, quantitative analyses have focused on the flux and energy spectrum of such an emission; polarisation has never been considered. Here we show that we could be missing out on an essential piece of information. The radiative decay of neutrinos, as well as the interactions of dark matter and neutrinos with Standard Model particles can generate a circular polarisation signal in X-rays or \gamma-rays. If observed, this could reveal important information about their spatial distribution and particle-antiparticle ratio, and could even reveal the nature of the high-energy particle physics processes taking place in astrophysical sites. The question of the observability of these polarised signatures and their separation from background astrophysical sources is left for future work.

In recent years, a revised set of solar abundances has led to a discrepancy in the sound-speed profile between helioseismology and theoretical solar models. Conventional solutions require additional mechanisms for energy transport within the Sun. Vincent et al. have recently suggested that dark matter with a momentum or velocity dependent cross section could provide a solution. In this work, we consider three models of dark matter with such cross sections and their effect on the stellar structure. In particular, the three models incorporate dark matter particles interacting through an electromagnetic dipole moment: an electric dipole, a magnetic dipole or an anapole. Each model is implemented in the \texttt{DarkStec} stellar evolution program, which incorporates the effects of dark matter capture and heat transport within the solar interior. We show that dark matter with an anapole moment of $\sim1\mathrm{GeV}^{-2}$ or magnetic dipole moment of $\sim10^{-3}\mu_p$ can improve the sound-speed profile, small frequency separations and convective zone radius with respect to the Standard Solar Model. However, the required dipole moments are strongly excluded by direct detection experiments.

We employ data from the recently observed high-energy neutrino events at the IceCube Neutrino Observatory to constrain interactions between the dark matter (DM) in the Milky Way and the neutrino sector. We construct an extended un-binned likelihood in order to explore the parameter space of allowed interactions. We present results in the specific case of a scalar DM candidate interacting via a scalar mediator, and show that due to the energy dependence of the interaction cross section, this approach can constrain the coupling more strongly than traditional cosmological probes for some regions of the parameter space.

Over the past few decades, an anomalous 511 keV gamma-ray line has been observed from the centre of the Milky Way. Dark matter (DM) in the form of light weakly interacting massive particles (WIMPs) annihilating into electron-positron pairs has been one of the leading hypotheses of the observed emission. Given the small required cross section, a further coupling to lighter particles is required to produce the correct relic density. Here, we derive constraints from the Planck satellite on light WIMPs that were in equilibrium with either the neutrino or electron sector in the early universe. For the neutrino sector, we obtain a lower bound on the WIMP mass of 4 MeV for a real scalar and 10 MeV for a Dirac fermion DM particle, at 95% CL. For the electron sector, we find even stronger bounds of 7 and 11 MeV, respectively. Using these results, we show that, in the absence of additional ingredients such as dark radiation, the light thermally produced WIMP explanation of the 511 keV excess is strongly disfavoured by the latest cosmological data. This suggests an unknown astrophysical or more exotic DM source of the signal.

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 20TeV. 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 power-law 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.

We present updated constraints on dark matter models with momentum-dependent or velocity-dependent interactions with nuclei, based on direct detection and solar physics. We improve our previous treatment of energy transport in the solar interior by dark matter scattering, leading to significant changes in fits to many observables. Based on solar physics alone, DM with a spin-independent $q^{4}$ coupling provides the best fit to data, and a statistically satisfactory solution to the solar abundance problem. Once direct detection limits are accounted for however, the best solution is spin-dependent $v^2$ scattering with a reference cross-section of 10$^{-35}$ cm$^2$ (at a reference velocity of $v_0=220$ km s$^{-1}$), and a dark matter mass of about 5 GeV.

The positron annihilation gamma-ray signal in the Milky Way (MW) shows a puzzling morphology: a very bright bulge and a very low surface-brightness disk. A coherent explanation of the positron origin, propagation through the Galaxy and subsequent annihilation in the interstellar medium has not yet been found. Tentative explanations involve positrons from radioactivity, X-ray binaries, and dark matter (DM). Dwarf satellite galaxies (DSGs) are believed to be DM-dominated and hence promising candidates in the search for 511 keV emission as a result of DM annihilation into electron-positron pairs. The goal of this study is to constrain possible 511 keV gamma-ray signals from 39 DSGs of the MW and to test the annihilating DM scenario. We use the spectrometer SPI on INTEGRAL to extract individual spectra for the studied objects. As the diffuse galactic emission dominates the signal, the large scale morphology of the MW has been modelled accordingly and was included in a maximum likelihood analysis. Alternatively, a distance-weighted stacked spectrum has been determined. Only Reticulum II (Ret II) shows a 3.1 sigma signal. Five other sources show tentative 2 sigma signals. The mass-to-511-keV-luminosity-ratio shows a marginal trend towards higher values for intrinsically brighter objects, opposite to the V band mass-to-light-ratio, which is generally used to uncover DM in DSGs. All derived flux values are above the level implied by a DM interpretation of the MW bulge signal. The signal from Ret II is unlikely to be related to a DM origin alone, otherwise, the MW bulge would be about 100 times brighter than what is seen. Ret II is exceptional considering the DSG sample, and rather points to enhanced recent star formation activity, if its origins are similar to processes in the MW. Understanding this emission may provide further clues regarding the origin of the annihilation emission in the MW.

We critically examine recent claims of a high solar metallicity by von Steiger \& Zurbuchen (2016) based on \textit{in situ} measurements of the solar wind, rather than the standard spectroscopically-inferred abundances (Asplund et al. 2009). We test the claim by Vagnozzi et al. (2016) that a composition based on the solar wind enables one to construct a standard solar model in agreement with helioseismological observations and thus solve the decades-old solar modelling problem. We show that, although some helioseismological observables are improved compared to models computed with spectroscopic abundances, most are in fact worse. The high abundance of refractory elements leads to an overproduction of neutrinos, with a predicted $^8$B flux that is nearly twice its observed value, and $^7$Be and CNO fluxes that are experimentally ruled out at high confidence. A combined likelihood analysis shows that models using the vSZ16 abundances fare worse than AGSS09 despite a higher metallicity. We also present astrophysical and spectroscopic arguments showing the vSZ16 composition to be an implausible representation of the solar interior, identifying the first ionisation potential effect in the outer solar atmosphere and wind as the likely culprit.

We study the effects of energy transport in the Sun by asymmetric dark matter with momentum and velocity-dependent interactions, with an eye to solving the decade-old Solar Abundance Problem. We study effective theories where the dark matter-nucleon scattering cross-section goes as $v_{\rm rel}^{2n}$ and $q^{2n}$ with $n = -1, 0, 1 $ or $2$, where $v_{\rm rel}$ is the dark matter-nucleon relative velocity and $q$ is the momentum exchanged in the collision. Such cross-sections can arise generically as leading terms from the most basic nonstandard DM-quark operators. We employ a high-precision solar simulation code to study the impact on solar neutrino rates, the sound speed profile, convective zone depth, surface helium abundance and small frequency separations. We find that the majority of models that improve agreement with the observed sound speed profile and depth of the convection zone also reduce neutrino fluxes beyond the level that can be reasonably accommodated by measurement and theory errors. However, a few specific points in parameter space yield a significant overall improvement. A 3-5 GeV DM particle with $\sigma_{SI} \propto q^2$ is particularly appealing, yielding more than a $6\sigma$ improvement with respect to standard solar models, while being allowed by direct detection and collider limits. We provide full analytical capture expressions for $q$- and $v_{\rm rel}$-dependent scattering, as well as complete likelihood tables for all models.

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 ($p$, $\bar p$ and He) from the usual light elements (Be, B, C, N, O). We find that the propagation parameters that best fit $p$, $\bar p$, He data are significantly different from those that fit light elements, including the B/C and $^{10}$Be/$^9$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.

Future dedicated radio interferometers, including HERA and SKA, are very promising tools that aim to study the epoch of reionization and beyond via measurements of the 21 cm signal from neutral hydrogen. Dark matter (DM) annihilations into charged particles change the thermal history of the Universe and, as a consequence, affect the 21 cm signal. Accurately predicting the effect of DM strongly relies on the modeling of annihilations inside halos. In this work, we use up-to-date computations of the energy deposition rates by the products from DM annihilations, a proper treatment of the contribution from DM annihilations in halos, as well as values of the annihilation cross section allowed by the most recent cosmological measurements from the Planck satellite. Given current uncertainties on the description of the astrophysical processes driving the epochs of reionization, X-ray heating and Lyman-$\alpha$ pumping, we find that disentangling DM signatures from purely astrophysical effects, related to early-time star formation processes or late-time galaxy X-ray emissions, will be a challenging task. We conclude that only annihilations of DM particles with masses of $\sim100$ MeV, could leave an unambiguous imprint on the 21 cm signal and, in particular, on the 21 cm power spectrum. This is in contrast to previous, more optimistic results in the literature, which have claimed that strong signatures might also be present even for much higher DM masses. Additional measurements of the 21 cm signal at different cosmic epochs will be crucial in order to break the strong parameter degeneracies between DM annihilations and astrophysical effects and undoubtedly single out a DM imprint for masses different from $\sim100$ MeV.

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 high-energy 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.

We present an in-depth analysis of the flavour and spectral composition of the 36 high-energy neutrino events observed after three years of observation by the IceCube neutrino telescope. While known astrophysical sources of HE neutrinos are expected to produce a nearly $(1:1:1)$ flavour ratio (electron : muon : tau) of neutrinos at earth, we show that the best fits based on the events detected above $E_\nu \ge 28$ TeV do not necessarily support this hypothesis. Crucially, the energy range that is considered when analysing the HE neutrino data can have a profound impact on the conclusions. We highlight two intriguing puzzles: an apparent deficit of muon neutrinos, seen via a deficit of track-like events; and an absence of $\bar \nu_e$'s at high energy, seen as an absence of events near the Glashow resonance. We discuss possible explanations, including the misidentification of tracks as showers, and a broken power law, in analogy to the observed HE cosmic ray spectrum.

Broad disagreement persists between helioseismological observables and predictions of solar models computed with the latest surface abundances. Here we show that most of these problems can be solved by the presence of asymmetric dark matter coupling to nucleons as the square of the momentum $q$ exchanged in the collision. We compute neutrino fluxes, small frequency separations, surface helium abundances, sound speed profiles and convective zone depths for a number of models, showing more than a $6\sigma$ preference for $q^2$ models over others, and over the Standard Solar Model. The preferred mass (3\,GeV) and reference dark matter-nucleon cross-section ($10^{-37}$\,cm$^2$ at $q_0 = 40$\,MeV) are within the region of parameter space allowed by both direct detection and collider searches.

We employ state-of-the art cosmological observables including supernova surveys and BAO information to provide constraints on the mass and mixing angle of a non-resonantly produced sterile neutrino species, showing that cosmology can effectively rule out sterile neutrinos which decay between BBN and the present day. The decoupling of an additional heavy neutrino species can modify the time dependence of the Universe's expansion between BBN and recombination and, in extreme cases, lead to an additional matter-dominated period; while this could naively lead to a younger Universe with a larger Hubble parameter, it could later be compensated by the extra radiation expected in the form of neutrinos from sterile decay. However, recombination-era observables including the Cosmic Microwave Background (CMB), the shift parameter $R_{CMB}$ and the sound horizon $r_s$ from Baryon Acoustic Oscillations (BAO) severely constrain this scenario. We self-consistently include the full time-evolution of the coupled sterile neutrino and standard model sectors in an MCMC, showing that if decay occurs after BBN, the sterile neutrino is essentially bounded by the constraint $\sin^2\theta \lesssim 0.026 (m_s/\mathrm{eV})^{-2}$.

The IceCube experiment has recently released 3 years of data of the first ever detected high-energy (>30 TeV) neutrinos, which are consistent with an extraterrestrial origin. In this talk, we compute the compatibility of the observed track-to-shower ratio with possible combinations of neutrino flavors with relative proportion (alpha_e:alpha_mu:alpha_tau). Although this observation is naively favored for the canonical (1:1:1) at Earth, once we consider the IceCube expectations for the atmospheric muon and neutrino backgrounds, this flavor combination presents some tension with data. We find that, for an astrophysical neutrino E_nu^{-2} energy spectrum, (1:1:1) at Earth is currently disfavored at 92% C.L. We discuss the trend of this result by comparing the results with the 2-year and 3-year data. We obtain the best-fit for (1:0:0) at Earth, which cannot be achieved from any flavor ratio at sources with averaged oscillations during propagation. Although it is not statistically significant at present, 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.

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 (alpha_e:alpha_mu:alpha_tau). 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 sigma C.L. Indeed, with the recently-released 3-year 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 use the formalism of Gould and Raffelt to compute the dimensionless thermal conduction coefficients for scattering of dark matter particles with standard model nucleons via cross-sections that depend on the relative velocity or momentum exchanged between particles. Motivated by models invoked to reconcile various recent results in direct detection, we explicitly compute the conduction coefficients $\alpha$ and $\kappa$ for cross-sections that go as $v_{\rm rel}^2$, $v_{\rm rel}^4$, $v_{\rm rel}^{-2}$, $q^2$, $q^4$ and $q^{-2}$, where $v_{\rm rel}$ is the relative DM-nucleus velocity and $q$ is the momentum transferred in the collision. We find that a $v_{\rm rel}^{-2}$ dependence can significantly enhance energy transport from the inner solar core to the outer core. The same can true for any $q$-dependent coupling, if the dark matter mass lies within some specific range for each coupling. This effect can complement direct searches for dark matter; combining these results with state-of-the-art Solar simulations should greatly increase sensitivity to certain DM models. It also seems possible that the so-called Solar Abundance Problem could be resolved by enhanced energy transport in the solar core due to such velocity- or momentum-dependent scatterings.

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 \lesssim 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-$\alpha$ measurements of the matter temperature at $z \sim 4$ to set a 95 % confidence level lower bound on the dark matter lifetime of $\sim 4 \times 10^{25}$ s for $m_\chi = 100$ MeV. This bound becomes lower by an order of magnitude at $m_\chi = 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.

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 mu+mu- 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 (<sigma v> = 3 x 10^{-26} cm^3 s^{-1}) for DM masses below 30 GeV and 15 GeV for the e+e- and mu+mu- 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.

The 511 keV electron-positron annihilation line, most recently characterized by the INTEGRAL/SPI experiment, is highly concentrated towards the Galactic centre. Its origin remains unknown despite decades of scrutiny. We propose a novel scenario in which known extragalactic positron sources such as radio jets of active galactic nuclei (AGN) fill the intergalactic medium with MeV e+e- pairs, which are then accreted into the Milky Way. We show that interpreting the diffuse cosmic radio background (CRB) as arising from radio sources with characteristics similar to the observed cores and radio lobes in powerful AGN jets suggests that the intergalactic positron-to-electron ratio could be as high as 10^{-5}, although this can be decreased if the CRB is not all produced by pairs and if not all positrons escape to the intergalactic medium. Assuming an accretion rate of one solar mass per year of matter into the Milky Way, a positron-to-electron ratio of only 10^{-7} is already enough to account for much of the 511 keV emission of the Galaxy. A simple spherical accretion model predicts an emission profile highly peaked in the central bulge, consistent with INTEGRAL observations. However, a realistic model of accretion with angular momentum would likely imply a more extended emission over the disk, with uncertainties depending on the magnetic field structure and turbulence in the galactic halo.

The ARCADE 2 collaboration has reported a significant excess in the isotropic radio background, whose homogeneity cannot be reconciled with clustered sources. This suggests a cosmological origin prior to structure formation. We investigate several potential mechanisms and show that injection of relativistic electrons through late decays of a metastable particle can give rise to the observed excess radio spectrum through synchrotron emission. However, constraints from the cosmic microwave background (CMB) anisotropy, on injection of charged particles and on the primordial magnetic field, present a challenge. The simplest scenario is with a >~ 9 GeV particle decaying into e+e- at a redshift of z ~ 5, in a magnetic field of ~ 5 microgauss, which exceeds the CMB B-field constraints, unless the field was generated after decoupling. Decays into exotic millicharged particles can alleviate this tension, if they emit synchroton radiation in conjunction with a sufficiently large background magnetic field of a dark U(1)' gauge field.

It is well known that current spectroscopic determinations of the chemical composition of the Sun are starkly at odds with the metallicity implied by helioseismology. We investigate whether the discrepancy may be due to conversion of photons to a new light boson in the solar photosphere. We examine the impact of particles with axion-like interactions with the photon on the inferred photospheric abundances, showing that resonant axion-photon conversion is not possible in the region of the solar atmosphere in which line-formation occurs. Although non-resonant conversion in the line-forming regions can in principle impact derived abundances, constraints from axion-photon conversion experiments rule out the couplings necessary for these effects to be detectable. We show that this extends to hidden photons and chameleons (which would exhibit similar phenomenological behaviour), ruling out known theories of new light bosons as photospheric solutions to the solar abundance problem.

We compare the full-sky morphology of the 511 keV gamma ray excess measured by the INTEGRAL/SPI experiment to predictions of models based on dark matter (DM) scatterings that produce low-energy positrons: either MeV-scale DM that annihilates directly into e+e- pairs, or heavy DM that inelastically scatters into an excited state (XDM) followed by decay into e+e- and the ground state.By direct comparison to the data, we find that such explanations are consistent with dark matter halo profiles predicted by numerical many-body simulations for a Milky Way-like galaxy. Our results favor an Einasto profile over the cuspier NFW distribution and exclude decaying dark matter scenarios whose predicted spatial distribution is too broad. We obtain a good fit to the shape of the signal using six fewer degrees of freedom than previous empirical fits to the 511 keV data. We find that the ratio of flux at Earth from the galactic bulge to that of the disk is between 1.9 and 2.4, taking into account that 73% of the disk contribution may be attributed to the beta decay of radioactive 26Al.

We study preheating in models where the inflaton has a non-canonical kinetic term, containing powers of the usual kinetic energy. The inflaton field oscillating about its potential minimum acts as a driving force for particle production through parametric resonance. Non-canonical kinetic terms can impose a speed limit on the motion of the inflaton, modifying the oscillating inflaton profile. This has two important effects: it turns a smooth sinusoidal profile into a sharp saw-tooth, enhancing resonance, and it lengthens the period of oscillations, suppressing resonance. We show that the second effect dominates over the first, so that preheating with a non-canonical inflaton field is less efficient than with canonical kinetic terms, and that the expansion of the Universe suppresses resonance even further.

We reconsider Sommerfeld-enhanced annihilation of dark matter (DM) into leptons to explain PAMELA and Fermi electron and positron observations, in light of possible new effects from substructure. There is strong tension between getting a large enough lepton signal while respecting constraints on the fluxes of associated gamma rays. We first show that these constraints become significantly more stringent than in previous studies when the contributions from background e^+ e^- are taken into account, so much so that even cored DM density profiles are ruled out. We then show how DM annihilations within subhalos can get around these constraints. Specifically, if most of the observed lepton excess comes from annihilations in a nearby (within 1 kpc) subhalo along a line of sight toward the galactic center, it is possible to match both the lepton and gamma ray observations. We demonstrate that this can be achieved in a simple class of particle physics models in which the DM annihilates via a hidden leptophilic U(1) vector boson, with explicitly computed Sommerfeld enhancement factors. Gamma ray constraints on the main halo annihilations (and CMB constraints from the era of decoupling) require the annihilating component of the DM to be subdominant, of order 10^-2 of the total DM density.

Numerical simulations of dark matter collapse and structure formation show that in addition to a large halo surrounding the baryonic component of our galaxy, there also exists a significant number of subhalos that extend hundreds of kiloparsecs beyond the edge of the observable Milky Way. We find that for dark matter (DM) annihilation models, galactic subhalos can significantly modify the spectrum of electrons and positrons as measured at our galactic position. Using data from the recent Via Lactea II simulation we include the subhalo contribution of electrons and positrons as boundary source terms for simulations of high energy cosmic ray propagation with a modified version of the publicly available GALPROP code. Focusing on the DM DM -> 4e annihilation channel, we show that including subhalos leads to a better fit to both the Fermi and PAMELA data. The best fit gives a dark matter particle mass of 1.2 TeV, for boost factors of 90 in the main halo and 1950-3800 in the subhalos (depending on assumptions about the background), in contrast to the 0.85 TeV mass that gives the best fit in the main halo-only scenario. These fits suggest that at least a third of the observed electron cosmic rays from DM annihilation could come from subhalos, opening up the possibility of a relaxation of recent stringent constraints from inverse Compton gamma rays originating from the high-energy leptons.

Using the gradient expansion method of Rigopoulos, Shellard and van Tent which treats cosmological perturbations as gradients on top of a homogeneous and isotropic FRW background, we study the production of nongaussianities in the roulette model of inflation. Investigating a number of trajectories within this two-field model of inflation, we find that while the superhorizon influence of the isocurvature modes on the curvature bispectrum produces nonzero contribution to f_NL, the effect is negligible next to the standard inflationary prediction |f_NL| ~ n_s - 1. This is the case in both the squeezed and equilateral configurations of the bispectrum, although the former is slightly larger in the trajectories under consideration.

Data assimilation techniques, developed in the last two decades mainly for weather prediction, produce better forecasts by taking advantage of both theoretical/numerical models and real-time observations. In this paper, we explore the possibility of applying the data-assimilation techniques known as 4D-VAR to the prediction of solar flares. We do so in the context of a continuous version of the classical cellular-automaton-based self-organized critical avalanche models of solar flares introduced by Lu and Hamilton (Astrophys. J., 380, L89, 1991). Such models, although a priori far removed from the physics of magnetic reconnection and magneto-hydrodynamical evolution of coronal structures, nonetheless reproduce quite well the observed statistical distribution of flare characteristics. We report here on a large set of data assimilation runs on synthetic energy release time series. Our results indicate that, despite the unpredictable (and unobservable) stochastic nature of the driving/triggering mechanism within the avalanche model, 4D-VAR succeeds in producing optimal initial conditions that reproduce adequately the time series of energy released by avalanches/flares. This is an essential first step towards forecasting real flares.