We present a study on the possibility of searching for long-lived supersymmetric partners with the MoEDAL experiment at the LHC. MoEDAL is sensitive to highly ionising objects such as magnetic monopoles or massive (meta)stable electrically charged particles. We focus on prospects of directly detecting long-lived sleptons in a phenomenologically realistic model which involves an intermediate neutral long-lived particle in the decay chain. This scenario is not yet excluded by the current data from ATLAS or CMS, and is compatible with astrophysical constraints. Using Monte Carlo simulation, we compare the sensitivities of MoEDAL versus ATLAS in scenarios where MoEDAL could provide discovery reach complementary to ATLAS and CMS, thanks to looser selection criteria combined with the virtual absence of background. It is also interesting to point out that, in such scenarios, in which charged staus are the main long-lived candidates, the relevant mass range for MoEDAL is compatible with a potential role of Supersymmetry in providing an explanation for the anomalous events observed by the ANITA detector.

Motivated by the recent galactic center gamma-ray excess identified in the Fermi-LAT data, we perform a detailed study of QCD fragmentation uncertainties in the modeling of the energy spectra of gamma-rays from Dark-Matter (DM) annihilation. When Dark-Matter particles annihilate to coloured final states, either directly or via decays such as $W^{(*)}\to q\bar{q}'$, photons are produced from a complex sequence of shower, hadronisation and hadron decays. In phenomenological studies, their energy spectra are typically computed using Monte Carlo event generators. These results have however intrinsic uncertainties due to the specific model used and the choice of model parameters, which are difficult to asses and which are typically neglected. We derive a new set of hadronisation parameters (tunes) for the \textsc{Pythia~8.2} Monte Carlo generator from a fit to LEP and SLD data at the $Z$ peak. For the first time, we also derive a conservative set of uncertainties on the shower and hadronisation model parameters. Their impact on the gamma-ray energy spectra is evaluated and discussed for a range of DM masses and annihilation channels. The spectra and their uncertainties are also provided in tabulated form for future use. The fragmentation-parameter uncertainties may be useful for collider studies as well.

We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.

We show that the Galactic Center Excess (GCE) emission, as recently updated by the Fermi-LAT Collaboration, could be explained by the sum of Fermi-bubbles-like emission plus dark matter (DM) annihilation, in the context of a scalar-singlet Higgs portal scenario (SHP). In fact, the standard SHP, where the DM particle, $S$, only has renormalizable interactions with the Higgs, is non-operational due to strong constraints, specially from DM direct detection limits. Thus we consider a most economical extension, called ESHP (for extended SHP), which simply consists in the addition of a second (heavier) scalar singlet in the dark sector. The second scalar can be integrated-out, leaving a standard SHP plus a dimension-6 operator. Essentially, this model has only two relevant parameters (the DM mass and the coupling of the dim-6 operator). DM annihilation occurs mainly into two Higgs bosons, $SS\rightarrow hh$. We demonstrate that, despite its economy, the ESHP model provides excellent fits to the GCE (with p-value $\sim 0.6-0.7$) for very reasonable values of the parameters, in particular $m_S \simeq 130$ GeV. This is achieved without conflicting with other observables and keeping the $S-$particle relic density at the accepted value for the DM content in the universe.

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.

We present a new method to interpret the $\gamma$-ray data of our inner Galaxy as measured by the Fermi Large Area Telescope (Fermi LAT). We train and test convolutional neural networks with simulated Fermi-LAT images based on models tuned to real data. We use this method to investigate the origin of an excess emission of GeV $\gamma$-rays seen in previous studies. Interpretations of this excess include $\gamma$ rays created by the annihilation of dark matter particles and $\gamma$ rays originating from a collection of unresolved point sources, such as millisecond pulsars. Our new method allows precise measurements of the contribution and properties of an unresolved population of $\gamma$-ray point sources in the interstellar diffuse emission model.

One of the most promising strategies to identify the nature of dark matter consists in the search for new particles at accelerators and with so-called direct detection experiments. Working within the framework of simplified models, and making use of machine learning tools to speed up statistical inference, we address the question of what we can learn about dark matter from a detection at the LHC and a forthcoming direct detection experiment. We show that with a combination of accelerator and direct detection data, it is possible to identify newly discovered particles as dark matter, by reconstructing their relic density assuming they are weakly interacting massive particles (WIMPs) thermally produced in the early Universe, and demonstrating that it is consistent with the measured dark matter abundance. An inconsistency between these two quantities would instead point either towards additional physics in the dark sector, or towards a non-standard cosmology, with a thermal history substantially different from that of the standard cosmological model.

The Fermi Collaboration has recently updated their analysis of gamma rays from the center of the Galaxy. They reconfirm the presence of an unexplained emission feature which is most prominent in the region of $1-10$ GeV, known as the Galactic Center GeV excess (GCE). Although the GCE is now firmly detected, an interpretation of this emission as a signal of self-annihilating dark matter (DM) particles is not unambiguously possible due to systematic effects in the gamma-ray modeling estimated in the Galactic Plane. In this paper we build a covariance matrix, collecting different systematic uncertainties investigated in the Fermi Collaboration's paper that affect the GCE spectrum. We show that models where part of the GCE is due to annihilating DM can still be consistent with the new data. We also re-evaluate the parameter space regions of the minimal supersymmetric Standard Model (MSSM) that can contribute dominantly to the GCE via neutralino DM annihilation. All recent constraints from DM direct detection experiments such as PICO, LUX, PandaX and Xenon1T, limits on the annihilation cross section from dwarf spheroidal galaxies and Large Hadron Collider limits are considered in this analysis. Due to a slight shift in the energy spectrum of the GC excess with respect to the previous Fermi analysis, and the recent limits from direct detection experiments, we find a slightly shifted parameter region of the MSSM compared to our previous analysis that is consistent with the GCE. Neutralinos with a mass between $85-220$ GeV can describe the excess via annihilation into a pair of $W$-bosons or top quarks. Remarkably, there are low fine-tuning models among the regions that we have found. The complete set of solutions will be probed by upcoming direct detection experiments and with dedicated searches in the upcoming data of the Large Hadron Collider.

The $\mu\nu$SSM solves the $\mu$ problem of supersymmetric models and reproduces neutrino data, simply using couplings with right-handed neutrinos $\nu$'s. Given that these couplings break explicitly $R$ parity, the gravitino is a natural candidate for decaying dark matter in the $\mu \nu$SSM. In this work we carry out a complete analysis of the detection of $\mu \nu$SSM gravitino dark matter through $\gamma$-ray observations. In addition to the two-body decay producing a sharp line, we include in the analysis the three-body decays producing a smooth spectral signature. We perform first a deep exploration of the low-energy parameter space of the $\mu \nu$SSM taking into account that neutrino data must be reproduced. Then, we compare the $\gamma$-ray fluxes predicted by the model with Fermi-LAT observations. In particular, with the 95$\%$ CL upper limits on the total diffuse extragalactic $\gamma$-ray background using 50 months of data, together with the upper limits on line emission from an updated analysis using 69.9 months of data. For standard values of bino and wino masses, gravitinos with masses larger than about 4 GeV, or lifetimes smaller than $10^{28}$ s, produce too large fluxes and are excluded as dark matter candidates. However, when limiting scenarios with large and close values of the gaugino masses are considered, the constraints turn out to be less stringent, excluding masses larger than 17 GeV and lifetimes smaller than $4\times 10^{25}$ s.

Global fit studies performed in the pMSSM and the photon excess signal originating from the Galactic Center seem to suggest compressed electroweak supersymmetric spectra with a $\sim$100 GeV bino-like dark matter particle. We find that these scenarios are not probed by traditional electroweak supersymmetry searches at the LHC. We propose to extend the ATLAS and CMS electroweak supersymmetry searches with an improved strategy for bino-like dark matter, focusing on chargino plus next-to-lightest neutralino production, with a subsequent decay into a tri-lepton final state. We explore the sensitivity for pMSSM scenarios with $\Delta m = m_{\rm NLSP} - m_{\rm LSP} \sim (5 - 50)$ GeV in the $\sqrt{s} = 14$ TeV run of the LHC. Counterintuitively, we find that the requirement of low missing transverse energy increases the sensitivity compared to the current ATLAS and CMS searches. With 300 fb$^{-1}$ of data we expect the LHC experiments to be able to discover these supersymmetric spectra with mass gaps down to $\Delta m \sim 9$ GeV for DM masses between 40 and 140 GeV. We stress the importance of a dedicated search strategy that targets precisely these favored pMSSM spectra.

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.

We present global fits of an effective field theory description of real, and complex scalar dark matter candidates. We simultaneously take into account all possible dimension 6 operators consisting of dark matter bilinears and gauge invariant combinations of quark and gluon fields. We derive constraints on the free model parameters for both the real (five parameters) and complex (seven) scalar dark matter models obtained by combining Planck data on the cosmic microwave background, direct detection limits from LUX, and indirect detection limits from the Fermi Large Area Telescope. We find that for real scalars indirect dark matter searches disfavour a dark matter particle mass below 100 GeV. For the complex scalar dark matter particle current data have a limited impact due to the presence of operators that lead to p-wave annihilation, and also do not contribute to the spin-independent scattering cross- section. Although current data are not informative enough to strongly constrain the theory parameter space, we demonstrate the power of our formalism to reconstruct the theoretical parameters compatible with an actual dark matter detection, by assuming that the excess of gamma rays observed by the Fermi Large Area Telescope towards the Galactic centre is entirely due to dark matter annihilations. Please note that the excess can very well be due to astrophysical sources such as millisecond pulsars. We find that scalar dark matter interacting via effective field theory operators can in principle explain the Galactic centre excess, but that such interpretation is in strong tension with the non-detection of gamma rays from dwarf galaxies in the real scalar case. In the complex scalar case there is enough freedom to relieve the tension.

Observations with the Fermi Large Area Telescope (LAT) indicate an excess in gamma rays originating from the center of our Galaxy. A possible explanation for this excess is the annihilation of Dark Matter particles. We have investigated the annihilation of neutralinos as Dark Matter candidates within the phenomenological Minimal Supersymmetric Standard Model (pMSSM). An iterative particle filter approach was used to search for solutions within the pMSSM. We found solutions that are consistent with astroparticle physics and collider experiments, and provide a fit to the energy spectrum of the excess. The neutralino is a Bino/Higgsino or Bino/Wino/Higgsino mixture with a mass in the range $84-92$~GeV or $87-97$~GeV annihilating into W bosons. A third solutions is found for a neutralino of mass $174-187$~GeV annihilating into top quarks. The best solutions yield a Dark Matter relic density $0.06 < \Omega h^2 <0.13$. These pMSSM solutions make clear forecasts for LHC, direct and indirect DM detection experiments. If the MSSM explanation of the excess seen by Fermi-LAT is correct, a DM signal might be discovered soon.

We present statistically convergent profile likelihood maps obtained via global fits of a phenomenological Minimal Supersymmetric Standard Model with 15 free parameters (the MSSM-15), based on over 250M points. We derive constraints on the model parameters from direct detection limits on dark matter, the Planck relic density measurement and data from accelerator searches. We provide a detailed analysis of the rich phenomenology of this model, and determine the SUSY mass spectrum and dark matter properties that are preferred by current experimental constraints. We evaluate the impact of the measurement of the anomalous magnetic moment of the muon ($g-2$) on our results, and provide an analysis of scenarios in which the lightest neutralino is a subdominant component of the dark matter. The MSSM-15 parameters are relatively weakly constrained by current data sets, with the exception of the parameters related to dark matter phenomenology ($M_1$, $M_2$, $\mu$), which are restricted to the sub-TeV regime, mainly due to the relic density constraint. The mass of the lightest neutralino is found to be < 1.5 TeV at 99% C.L., but can extend up to 3 TeV when excluding the $g - 2$ constraint from the analysis. Low-mass bino-like neutralinos are strongly favoured, with spin-independent scattering cross-sections extending to very small values, $\sim 10^{-20}$ pb. ATLAS SUSY null searches strongly impact on this mass range, and thus rule out a region of parameter space that is outside the reach of any current or future direct detection experiment. The best-fit point obtained after inclusion of all data corresponds to a squark mass of 2.3 TeV, a gluino mass of 2.1 TeV and a 130 GeV neutralino with a spin-independent cross-section of $2.4 \times 10^{-10}$ pb, which is within the reach of future multi-ton scale direct detection experiments and of the upcoming LHC run at increased centre-of-mass energy.

We study in detail the impact of the current uncertainty in nucleon matrix elements on the sensitivity of direct and indirect experimental techniques for dark matter detection. We perform two scans in the framework of the cMSSM: one using recent values of the pion-sigma term obtained from Lattice QCD, and the other using values derived from experimental measurements. The two choices correspond to extreme values quoted in the literature and reflect the current tension between different ways of obtaining information about the structure of the nucleon. All other inputs in the scans, astrophysical and from particle physics, are kept unchanged. We use two experiments, XENON100 and IceCube, as benchmark cases to illustrate our case. We find that the interpretation of dark matter search results from direct detection experiments is more sensitive to the choice of the central values of the hadronic inputs than the results of indirect search experiments. The allowed regions of cMSSM parameter space after including XENON100 constrains strongly differ depending on the assumptions on the hadronic matrix elements used. On the other hand, the constraining potential of IceCube is almost independent of the choice of these values.

We present global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM) and the Non-Universal Higgs Model (NUHM), including the most recent CMS constraint on the Higgs boson mass, 5.8/fb integrated luminosity null Supersymmetry searches by ATLAS, the new LHCb measurement of the Bs to mu+mu- branching ratio and the 7-year WMAP dark matter relic abundance determination. We include the latest dark matter constraints from the XENON100 experiment, marginalising over astrophysical and particle physics uncertainties. We present Bayesian posterior and profile likelihood maps of the highest resolution available today, obtained from up to 350M points. We find that the new constraint on the Higgs boson mass has a dramatic impact, ruling out large regions of previously favoured cMSSM and NUHM parameter space. In the cMSSM, light sparticles and predominantly gaugino-like dark matter with a mass of a few hundred GeV are favoured. The NUHM exhibits a strong preference for heavier sparticle masses and a Higgsino-like neutralino with a mass of 1 TeV. The future ton-scale XENON1T direct detection experiment will probe large portions of the currently favoured cMSSM and NUHM parameter space. The LHC operating at 14 TeV collision energy will explore the favoured regions in the cMSSM, while most of the regions favoured in the NUHM will remain inaccessible. Our best-fit points achieve a satisfactory quality-of-fit, with p-values ranging from 0.21 to 0.35, so that none of the two models studied can be presently excluded at any meaningful significance level.

We present new global fits of the cMSSM, including LHC 1/fb integrated luminosity SUSY exclusion limits, recent LHC 5/fb constraints on the mass of the Higgs boson and XENON100 direct detection data. Our analysis fully takes into account astrophysical and hadronic uncertainties that enter the analysis when translating direct detection limits into constraints on the cMSSM parameter space. We provide results for both a Bayesian and a Frequentist statistical analysis. We find that LHC 2011 constraints in combination with XENON100 data can rule out a significant portion of the cMSSM parameter space. Our results further emphasise the complementarity of collider experiments and direct detection searches in constraining extensions of Standard Model physics. The LHC 2011 exclusion limit strongly impacts on low-mass regions of cMSSM parameter space, such as the stau co-annihilation region, while direct detection data can rule out regions of high SUSY masses, such as the Focus-Point region, which is unreachable for the LHC in the near future. We show that, in addition to XENON100 data, the experimental constraint on the anomalous magnetic moment of the muon plays a dominant role in disfavouring large scalar and gaugino masses. We find that, should the LHC 2011 excess hinting towards a Higgs boson at 126 GeV be confirmed, currently favoured regions of the cMSSM parameter space will be robustly ruled out from both a Bayesian and a profile likelihood statistical perspective.

We present updated global fits of the constrained Minimal Supersymmetric Standard Model (cMSSM), including the most recent constraints from the ATLAS and CMS detectors at the LHC, as well as the most recent results of the XENON100 experiment. Our robust analysis takes into account both astrophysical and hadronic uncertainties that enter in the calculation of the rate of WIMP-induced recoils in direct detection experiment. We study the consequences for neutralino Dark Matter, and show that current direct detection data already allow to robustly rule out the so-called Focus Point region, therefore demonstrating the importance of particle astrophysics experiments in constraining extensions of the Standard Model of Particle Physics. We also observe an increased compatibility between results obtained from a Bayesian and a Frequentist statistical perspective. We find that upcoming ton-scale direct detection experiments will probe essentially the entire currently favoured region (at the 99% level), almost independently of the statistical approach used. Prospects for indirect detection of the cMSSM are further reduced.

The unfortunate case where the Large Hadron Collider (LHC) fails to discover physics Beyond the Standard Model (BSM) is sometimes referred to as the "Nightmare scenario" of particle physics. We study the consequences of this hypothetical scenario for Dark Matter (DM), in the framework of the constrained Minimal Supersymmetric Standard Model (cMSSM). We evaluate the surviving regions of the cMSSM parameter space after null searches at the LHC, using several different LHC configurations, and study the consequences for DM searches with ton-scale direct detectors and the IceCube neutrino telescope. We demonstrate that ton-scale direct detection experiments will be able to conclusively probe the cMSSM parameter space that would survive null searches at the LHC with 100fb$^{-1}$ of integrated luminosity at 14TeV. We also demonstrate that IceCube (80 strings plus DeepCore) will be able to probe as much as 17% of the currently favoured parameter space after 5 years of observation.

Even if Supersymmetric particles are found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the bulk of the Dark Matter (DM) in the Universe using LHC data alone. We study the complementarity of LHC and DM indirect searches, working out explicitly the reconstruction of the DM properties for a specific benchmark model in the coannihilation region of a 24-parameters supersymmetric model. Combining mock high-luminosity LHC data with present-day null searches for gamma-rays from dwarf galaxies with the Fermi LAT, we show that current Fermi LAT limits already have the capability of ruling out a spurious Wino-like solution that would survive using LHC data only, thus leading to the correct identification of the cosmological solution. We also demonstrate that upcoming Planck constraints on the reionization history will have a similar constraining power, and discuss the impact of a possible detection of gamma-rays from DM annihilation in Draco with a CTA-like experiment. Our results indicate that indirect searches can be strongly complementary to the LHC in identifying the DM particles, even when astrophysical uncertainties are taken into account.

Motivated by the recent re-confirmation by CoGENT of the low-energy excess of events observed last year, and the recent improved limits from the XENON-100 experiment that are in contention with the CoGENT data, we re-examine the low mass neutralino region of the Minimal Supersymmetric Standard Model and of the Next-to-Minimal Supersymmetric Standard Model, both without assuming gaugino mass unification. We make several focused scans for each model, determining conservative constraints on input parameters. We then determine how these constraints are made increasingly stringent as we re-invoke our experimental constraints involving the dark matter relic abundance, collider constraints from LEP and the Tevatron, and then from flavour physics, as a series of successive 2 sigma hard cuts. We find that for both models, when all relevant constraints are applied in this fashion, we do not generate neutralino LSPs that possess a spin-independent scattering cross section in excess of 10^-5 pb and a mass ~7 GeV < m_chi < ~9 GeV that is necessary in order to explain the CoGENT observations.

We investigate the reconstruction capabilities of Dark Matter mass and spin-independent cross-section from future ton-scale direct detection experiments using germanium, xenon or argon as targets. Adopting realistic values for the exposure, energy threshold and resolution of Dark Matter experiments which will come online within 5 to 10 years, the degree of complementarity between different targets is quantified. We investigate how the uncertainty in the astrophysical parameters controlling the local Dark Matter density and velocity distribution affects the reconstruction. For a 50 GeV WIMP, astrophysical uncertainties degrade the accuracy in the mass reconstruction by up to a factor of $\sim 4$ for xenon and germanium, compared to the case when astrophysical quantities are fixed. However, combination of argon, germanium and xenon data increases the constraining power by a factor of $\sim 2$ compared to germanium or xenon alone. We show that future direct detection experiments can achieve self-calibration of some astrophysical parameters, and they will be able to constrain the WIMP mass with only very weak external astrophysical constraints.

Research in many areas of modern physics such as, e.g., indirect searches for dark matter and particle acceleration in SNR shocks, rely heavily on studies of cosmic rays (CRs) and associated diffuse emissions (radio, microwave, X-rays, gamma rays). While very detailed numerical models of CR propagation exist, a quantitative statistical analysis of such models has been so far hampered by the large computational effort that those models require. Although statistical analyses have been carried out before using semi-analytical models (where the computation is much faster), the evaluation of the results obtained from such models is difficult, as they necessarily suffer from many simplifying assumptions, The main objective of this paper is to present a working method for a full Bayesian parameter estimation for a numerical CR propagation model. For this study, we use the GALPROP code, the most advanced of its kind, that uses astrophysical information, nuclear and particle data as input to self-consistently predict CRs, gamma rays, synchrotron and other observables. We demonstrate that a full Bayesian analysis is possible using nested sampling and Markov Chain Monte Carlo methods (implemented in the SuperBayeS code) despite the heavy computational demands of a numerical propagation code. The best-fit values of parameters found in this analysis are in agreement with previous, significantly simpler, studies also based on GALPROP.

In theories with Universal Extra-Dimensions (UED), the gamma_1 particle, first excited state of the hypercharge gauge boson, provides an excellent Dark Matter (DM) candidate. Here we use a modified version of the SuperBayeS code to perform a Bayesian analysis of the minimal UED scenario, in order to assess its detectability at accelerators and with DM experiments. We derive in particular the most probable range of mass and scattering cross sections off nucleons, keeping into account cosmological and electroweak precision constraints. The consequences for the detectability of the gamma_1 with direct and indirect experiments are dramatic. The spin-independent cross section probability distribution peaks at ~ 10^{-11} pb, i.e. below the sensitivity of ton-scale experiments. The spin-dependent cross-section drives the predicted neutrino flux from the center of the Sun below the reach of present and upcoming experiments. The only strategy that remains open appears to be direct detection with ton-scale experiments sensitive to spin-dependent cross-sections. On the other hand, the LHC with 1 1/fb of data should be able to probe the current best-fit UED parameters.

Dark matter (DM) is currently searched for with a variety of detection strategies. Accelerator searches are particularly promising, but even if Weakly Interacting Massive Particles (WIMPs) are found at the Large Hadron Collider (LHC), it will be difficult to prove that they constitute the bulk of the DM in the Universe. We show that a significantly better reconstruction of the DM properties can be obtained with a combined analysis of LHC and direct detection (DD) data, by making a simple Ansatz on the WIMP local density, i.e. by assuming that the local density scales with the cosmological relic abundance. We demonstrate this method in an explicit example in the context of a 24-parameter supersymmetric model, with a neutralino LSP in the stau co-annihilation region. Our results show that future ton-scale DD experiments will allow to break degeneracies in the SUSY parameter space and achieve a significantly better reconstruction of the neutralino composition and its relic density than with LHC data alone.

In this work, we explore the "degenerate gravitino" scenario where the mass difference between the gravitino and the lightest MSSM particle is much smaller than the gravitino mass itself. In this case, the energy released in the decay of the next to lightest sypersymmetric particle (NLSP) is reduced. Consequently the cosmological and astrophysical constraints on the gravitino abundance, and hence on the reheating temperature, become softer than in the usual case. On the other hand, such small mass splittings generically imply a much longer lifetime for the NLSP. We find that, in the constrained MSSM (CMSSM), for neutralino LSP or NLSP, reheating temperatures compatible with thermal leptogenesis are reached for small splittings of order 10^{-2} GeV. While for stau NLSP, temperatures of 4x10^9 GeV can be obtained even for splittings of order of tens of GeVs. This "degenerate gravitino" scenario offers a possible way out to the gravitino problem for thermal leptogenesis in supersymmetric theories.

The $\mu\nu$SSM provides a solution to the $\mu$-problem of the MSSM and explains the origin of neutrino masses by simply using right-handed neutrino superfields. Given that R-parity is broken in this model, the gravitino is a natural candidate for dark matter since its lifetime becomes much longer than the age of the Universe. We consider the implications of gravitino dark matter in the $\mu\nu$SSM, analyzing in particular the prospects for detecting gamma rays from decaying gravitinos. If the gravitino explains the whole dark matter component, a gravitino mass larger than 20 GeV is disfavored by the isotropic diffuse photon background measurements. On the other hand, a gravitino with a mass range between 0.1-20 GeV gives rise to a signal that might be observed by the FERMI satellite. In this way important regions of the parameter space of the $\mu\nu$SSM can be checked.

We study in detail the ability of the nominal configuration of the IceCube neutrino telescope (with 80 strings) to probe the parameter space of the Constrained MSSM (CMSSM) favoured by current collider and cosmological data. Adopting conservative assumptions about the galactic halo model and the expected experiment performance, we find that IceCube has a probability between 2% and 12% of achieving a 5sigma detection of dark matter annihilation in the Sun, depending on the choice of priors for the scalar and gaugino masses and on the astrophysical assumptions. We identify the most important annihilation channels in the CMSSM parameter space favoured by current constraints, and we demonstrate that assuming that the signal is dominated by a single annihilation channel canlead to large systematic errors in the inferred WIMP annihilation cross section. We demonstrate that ~ 66% of the CMSSM parameter space violates the equilibrium condition between capture and annihilation in the center of the Sun. By cross-correlating our predictions with direct detection methods, we conclude that if IceCube does detect a neutrino flux from the Sun at high significance while direct detection experiments do not find a signal above a spin-independent cross section sigma_SI^p larger than 5x10^{-9} pb, the CMSSM will be strongly disfavoured, given standard astrophysical assumptions for the WIMP distribution. This result is robust with respect to a change of priors. We argue that the proposed low-energy DeepCore extension of IceCube will be an ideal instrument to focus on relevant CMSSM areas of parameter space.

In the framework of the Constrained MSSM we derive the most probable ranges of the diffuse gamma radiation flux from the direction of the Galactic center and of the positron flux from the Galactic halo due to neutralino dark matter annihilation. We find that, for a given halo model, and assuming flat priors, the 68% probability range of the integrated gamma-ray flux spans about one order of magnitude, while the 95% probability range can be much larger and extend over four orders of magnitude (even exceeding five for a tiny region at small neutralino mass). The detectability of the signal by GLAST depending primarily on the cuspiness of the halo profile. The positron flux, on the other hand, appears to be too small to be detectable by PAMELA, unless the boost factor is at least of order ten and/or the halo profile is extremely cuspy. We also briefly discuss the sensitivity of our results to the choice of priors.

Standard Big Bang nucleosynthesis at baryonic density as inferred by WMAP implies a primordial Li7 abundance factor of two to three larger than that inferred by observations of low--metallicity halo stars. Recent observations of Li6 in halo stars suggest a pre--galactic origin of this isotope, which is exceedingly difficult to explain by putative high redshift cosmic ray populations. We investigate if one or both of these lithium problems may be solved by late decaying relic particles in the Constrained Minimal Supersymmetric Standard Model (CMSSM) coupled to gravity. Assuming that the gravitino is the lightest supersymmetric particle (LSP) we find that in large parts of the CMSSM parameter space either of these problems can be potentially solved. In particular, we also find solutions where both lithium problems may be solved simultaneously. These entail the hadronic decay of relic 1 TeV staus into 100 GeV gravitinos at 1000 sec after the Big Bang, as proposed by one of us before. Moreover, the produced gravitinos naturally contribute a large fraction, or all, to the by WMAP required dark matter density. A prediction of this model is the dark matter to be lukewarm with free--streaming lengths of the order of a 4 keV early freezing--out relic particle. Such lukewarm gravitino dark matter may have beneficial effects to the formation of galactic structures and may be constrained by the reionisation history of the Universe. The current version of the paper presents results for both cases (a) when catalytic nuclear reations are included and (b) when they are neglected.

After a period of inflationary expansion, the Universe reheated and reached full thermal equilibrium at the reheating temperature T_R. In this work we point out that, in the context of effective low-energy supersymmetric models, LHC measurements may allow one to determine T_R as a function of the mass of the dark matter particle assumed to be either an axino or a gravitino. An upper bound on their mass may also be derived.

We re-examine the properties of the Constrained MSSM in light of updated constraints, paying particular attention to the impact of the recent substantial shift in the Standard Model prediction for BR(B to X_s gamma). With the help of a Markov Chain Monte Carlo scanning technique, we vary all relevant parameters simultaneously and derive Bayesian posterior probability maps. We find that the case of \mu>0 remains favored, and that for \mu<0 it is considerably more difficult to find a good global fit to current constraints. In both cases we find a strong preference for a focus point region. This leads to improved prospects for detecting neutralino dark matter in direct searches, while superpartner searches at the LHC become more problematic, especially when \mu<0. In contrast, prospects for exploring the whole mass range of the lightest Higgs boson at the Tevatron and the LHC remain very good, which should, along with dark matter searches, allow one to gain access to the otherwise experimentally challenging focus point region. An alternative measure of the mean quality-of-fit which we also employ implies that present data are not yet constraining enough to draw more definite conclusions. We also comment on the dependence of our results on the choice of priors and on some other assumptions.

We use cosmic microwave background and large scale structure data to test a broad and physically well-motivated class of inflationary models: those with flat tree-level potentials (typical in supersymmetry). The non-trivial features of the potential arise from radiative corrections which give a simple logarithmic dependence on the inflaton field, making the models very predictive. We also consider a modified scenario with new physics beyond a certain high-energy cut-off showing up as non-renormalizable operators (NRO) in the inflaton field. We find that both kinds of models fit remarkably well CMB and LSS data, with very few free parameters. Besides, a large part of these models naturally predict a reasonable number of e-folds. A robust feature of these scenarios is the smallness of tensor perturbations (r < 10^{-3}). The NRO case can give a sizeable running of the spectral index while achieving a sufficient number of e-folds. We use Bayesian model comparison tools to assess the relative performance of the models. We believe that these scenarios can be considered as a standard physical class of inflationary models, on a similar footing with monomial potentials.

In the framework of the Constrained MSSM we re--examine the gravitino as the lightest superpartner and a candidate for cold dark matter in the Universe. Unlike in other recent studies, we include both a thermal contribution to its relic population from scatterings in the plasma and a non--thermal one from neutralino or stau decays after freeze--out. Relative to a previous analysis [1] we update, extend and considerably improve our treatment of constraints from observed light element abundances on additional energy released during BBN in association with late gravitino production. Assuming the gravitino mass in the GeV to TeV range, and for natural ranges of other supersymmetric parameters, the neutralino region is excluded, while for smaller values of the gravitino mass it becomes allowed again. The gravitino relic abundance is consistent with observational constraints on cold dark matter from BBN and CMB in some well defined domains of the stau region but, in most cases, only due to a dominant contribution of the thermal population. This implies, depending on the gravitino mass, a large enough reheating temperature. If $\mgravitino>1$ GeV then $T_R>10^7$ GeV, if allowed by BBN and other constraints but, for light gravitinos, if $\mgravitino>100$ keV then $T_R>3\times 10^3$ GeV. On the other hand, constraints mostly from BBN imply an upper bound $T_R \lsim {a few}x 10^8\times10^9$ GeV which appears inconsistent with thermal leptogenesis. Finally, most of the preferred stau region corresponds to the physical vacuum being a false vacuum. The scenario can be partially probed at the LHC.

We perform a comprehensive exploration of the Constrained MSSM parameters using a Markov Chain Monte Carlo technique and a Bayesian analysis. We compute superpartner masses and other collider observables as well as a cold dark matter abundance, and compare them with experimental data. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant SM parameters. We delineate probability distributions of the CMSSM parameters, the collider and cosmological observables as well as a dark matter direct detection cross section. The 68% probability intervals of the CMSSM parameters are: 0.52 TeV < m_{1/2} < 1.26 TeV, m_0 <2.10 TeV, -0.34 TeV < A_0 < 2.41 TeV and 38.5< tan(beta) <54.6. Generally, large fractions of high probability ranges of the superpartner masses will be probed at the LHC. We highlight a complementarity between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. We further expose a number of correlations among the observables, in particular between BR(B_s \to \mu^+ \mu^-) and BR({\bar B}\to X_s\gamma) or sigma_p^{SI}. Once SUSY is discovered, this and other correlations may prove helpful in distinguishing the CMSSM from other supersymmetric models. The robustness of our results is investigated in terms of the assumed ranges of CMSSM parameters and the effect of the (g-2)_mu anomaly which shows some tension with the other observables. We find that the results for m_0, and the observables which strongly depend on it, are sensitive to our assumptions, while our conclusions for the other variables are robust.

We outline the WIMP dark matter parameter space in the Constrained MSSM by performing a comprehensive statistical analysis that compares with experimental data predicted superpartner masses and other collider observables as well as a cold dark matter abundance. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant Standard Model parameters. We present high-probability regions for neutralino dark matter direct detection cross section, and we find that 10^-10 pb < sigma_p^SI < 10^-8 for direct WIMP detection (with details slightly dependent on the assumptions made). We highlight a complementarity between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. We conclude that most of the 95% probability region for the cross section will be explored by future one-tonne detectors, that will therefore cover most of the currently favoured region of parameter space.

In the framework of the CMSSM we study the gravitino as the lightest supersymmetric particle and the dominant component of cold dark matter in the Universe. We include both a thermal contribution to its relic abundance from scatterings in the plasma and a non--thermal one from neutralino or stau decays after freeze--out. In general both contributions can be important, although in different regions of the parameter space. We further include constraints from BBN on electromagnetic and hadronic showers, from the CMB blackbody spectrum and from collider and non--collider SUSY searches. The region where the neutralino is the next--to--lightest superpartner is severely constrained by a conservative bound from excessive electromagnetic showers and probably basically excluded by the bound from hadronic showers, while the stau case remains mostly allowed. In both regions the constraint from CMB is often important or even dominant. In the stau case, for the assumed reasonable ranges of soft SUSY breaking parameters, we find regions where the gravitino abundance is in agreement with the range inferred from CMB studies, provided that, in many cases, a reheating temperature $\treh$ is large, $\treh\sim10^{9}\gev$. On the other side, we find an upper bound $\treh\lsim 5\times 10^{9}\gev$. Less conservative bounds from BBN or an improvement in measuring the CMB spectrum would provide a dramatic squeeze on the whole scenario, in particular it would strongly disfavor the largest values of $\treh\sim 10^{9}\gev$. The regions favored by the gravitino dark matter scenario are very different from standard regions corresponding to the neutralino dark matter, and will be partly probed at the LHC.

We consider a supersymmetric grand unified model which naturally solves the strong CP and mu problems via a Peccei-Quinn symmetry and leads to the standard realization of hybrid inflation. We show that the Peccei-Quinn field of this model can act as curvaton. In contrast to the standard curvaton hypothesis, both the inflaton and the curvaton contribute to the total curvature perturbation. The model predicts an isocurvature perturbation too which has mixed correlation with the adiabatic one. The cold dark matter of the universe is mostly constituted by axions plus a small amount of lightest sparticles. The predictions of the model are confronted with the Wilkinson microwave anisotropy probe and other cosmic microwave background radiation data. We analyze two representative choices of parameters and derive bounds on the curvaton contribution to the adiabatic perturbation. We find that, for the choice which provides the best fitting of the data, the curvaton contribution to the adiabatic amplitude must be smaller than about 67% (at 95% confidence level). The best-fit power spectra are dominated by the adiabatic part of the inflaton contribution. We use Bayesian model comparison to show that this choice of parameters is disfavored with respect to the pure inflaton scale-invariant case with odds of 50 to 1. For the second choice of parameters, the adiabatic mode is dominated by the curvaton, but this choice is strongly disfavored relative to the pure inflaton scale-invariant case (with odds of 10^7 to 1). We conclude that in the present framework the perturbations must be dominated by the adiabatic component from the inflaton.

If the axino is the lightest superpartner and satisfies cosmological bounds, including a preferred range of the relic abundance of cold dark matter, then the usual stringent constraints on the parameter space of the CMSSM become greatly relaxed. The lightest superpartner of the usual CMSSM spectrum will appear to be stable in collider experiments but will not necessarily obey relic abundance constraints. It may be either neutral (lightest neutralino) or charged (typically a stau). With the axino as cold dark matter, large regions of the CMSSM, often corresponding to heavy superpartners, become allowed, depending on the axino mass and the reheating temperature.

In contrast to the inflaton's case, the curvature perturbations due to the curvaton field depend strongly on the evolution of the curvaton before its decay. We study in detail the dynamics of the curvaton evolution during and after inflation. We consider that the flatness of the curvaton potential may be affected by supergravity corrections, which introduce an effective mass proportional to the Hubble parameter. We also consider that the curvaton potential may be dominated by a quartic or by a non-renormalizable term. We find analytic solutions for the curvaton's evolution for all these possibilities. In particular, we show that, in all the above cases, the curvaton's density ratio with respect to the background density of the Universe decreases. Therefore, it is necessary that the curvaton decays only after its potential becomes dominated by the quadratic term, which results in (Hubble damped) sinusoidal oscillations. In the case when a non-renormalizable term dominates the potential, we find a possible non-oscillatory attractor solution that threatens to erase the curvature perturbation spectrum. Finally, we study the effects of thermal corrections to the curvaton's potential and show that, if they ever dominate the effective mass, they lead to premature thermalization of the curvaton condensate. To avoid this danger, a stringent bound has to be imposed on the coupling of the curvaton to the thermal bath.

A simple extension of the minimal supersymmetric standard model which naturally and simultaneously solves the strong CP and mu problems via a Peccei-Quinn and a continuous R symmetry is considered. This model is supplemented with hybrid inflation and leptogenesis, but without taking the specific details of these scenarios. It is shown that the Peccei-Quinn field can successfully act as a curvaton generating the total curvature perturbation in the universe in accord with the cosmic background explorer measurements. A crucial phenomenon, which assists us to achieve this, is the "tachyonic amplification" of the perturbation acquired by this field during inflation if the field, in its subsequent evolution, happens to be stabilized for a while near a maximum of the potential. In this case, the contribution of the field to the total energy density is also enhanced ("tachyonic effect"), which helps too. The cold dark matter in the universe consists, in this model, mainly of axions which carry an isocurvature perturbation uncorrelated with the total curvature perturbation. There are also lightest sparticles (neutralinos) which, like the baryons, originate from the inflationary reheating and, thus, acquire an isocurvature perturbation fully correlated with the curvature perturbation. So, the overall isocurvature perturbation has a mixed correlation with the adiabatic one. It is shown that the presently available bound on such an isocurvature perturbation from cosmic microwave background radiation and other data is satisfied. Also, the constraint on the non-Gaussianity of the curvature perturbation obtained from the recent Wilkinson microwave anisotropy probe data is fulfilled thanks to the "tachyonic effect".

We derive the allowed ranges of the spin--independent interaction cross section $\sigsip$ for the elastic scattering of neutralinos on proton for wide ranges of parameters of the general Minimal Supersymmetric Standard Model. We investigate the effects of the lower limits on Higgs and superpartner masses from colliders, as well as the impact of constraints from $\bsgamma$ and the new measurement of $\gmtwo$ on the upper and lower limits on $\sigsip$. We further explore the impact of the neutralino relic density, including coannihilation, and of theoretical assumptions about the largest allowed values of the supersymmetric parameters. For $\mu>0$, requiring the latter to lie below $1\tev$ leads to $\sigsip\gsim 10^{-11}\pb$ at $\mchi\sim100\gev$ and $\sigsip\gsim 10^{-8}\pb$ at $\mchi\sim1\tev$. When the supersymmetric parameters are allowed above $1\tev$, for $440\gev \lsim \mchi\lsim 1020 \gev$ we derive a {\em parameter--independent lower limit} of $\sigsip \gsim 2\times 10^{-12}\pb$. (No similar lower limits can be set for $\mu<0$ nor for $1020\gev\lsim\mchi\lsim2.6\tev$.) Requiring $\abundchi<0.3$ implies a {\em parameter--independent upper limit} $\mchi\lsim2.6\tev$. The new $\epem$--based measurement of $(g-2)_{\mu}$ restricts $\mchi\lsim 350\gev$ at $1 \sigma$ CL and $\mchi\lsim515\gev$ at $2 \sigma$ CL, and implies $\mu>0$. The largest allowed values of $\sigsip$ have already become accessible to recent experimental searches.

Coannihilation processes provide an important additional mechanism for reducing the density of stable relics in the Universe. In the case of the stable lightest neutralino of the MSSM, and in particular the Constrained MSSM (CMSSM), the coannihilation with sleptons plays a major role in opening up otherwise cosmologically excluded ranges of supersymmetric parameters. In this paper, we derive a full set of exact, analytic expressions for the coannihilation of the lightest neutralino with the sleptons into all two--body tree--level final states in the framework of minimal supersymmetry. We make no simplifying assumptions about the neutralino nor about sfermion masses and mixings other than the absence of explicit CP--violating terms and inter--family mixings. The expressions should be particularly useful in computing the neutralino WIMP relic abundance without the approximation of partial wave expansion. We illustrate the effect of our analytic results with numerical examples and demonstrate a sizeable difference with approximate expressions available in the literature.

We compute the neutralino relic density in the minimal supersymmetric standard model by using exact expressions for the neutralino annihilation cross section into all tree-level final states, including all contributions and interference terms. We find that several final states may give comparable contributions to the relic density, which illustrates the importance of performing a complete calculation. We compare the exact results with those of the usual expansion method and demonstrate a sizeable discrepancy (of more than 10%) over a significant range of the neutralino mass of up to several tens of GeV which is caused by the presence of resonances and new final-state thresholds. We perform several related checks and comparisons. In particular, we find that the often employed approximate iterative procedure of computing the neutralino freeze-out temperature gives generally very accurate results, except when the expansion method is used near resonances and thresholds.