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 provide an updated and improved study of the prospects of the H.E.S.S. and Cherenkov Telescope Array (CTA) experiments in testing neutralino dark matter in the Minimal Supersymmetric Standard Model with nine free parameters (p9MSSM). We include all relevant experimental constraints and theoretical developments, in particular the calculation of the Sommerfeld enhancement for both present-day annihilation and the relic abundance. We perform a state-of-the-art analysis of the CTA sensitivity with a log-likelihood test ratio statistics and apply it to a numerical scan of the p9MSSM parameter space focusing on a TeV scale dark matter. We find that, assuming Einasto profile of dark matter halo in the Milky Way, H.E.S.S. has already been able to nearly reach the so-called thermal WIMP value, while CTA will go below it by providing a further improvement of at least an order of magnitude. Both H.E.S.S. and CTA are sensitive to several cases for which direct detection cross section will be below the so-called neutrino floor, with H.E.S.S. being sensitive to most of the wino region, while CTA also covering a large fraction of the ~1 TeV higgsino region. We show that CTA sensitivity will be further improved in the monochromatic photon search mode for both single-component and underabundant dark matter.

Motivated by the recent anomalies in $b\rightarrow s$ transitions that emerged at LHCb, we consider a model with an $L_\mu - L_\tau$ gauge symmetry and additional vector-like fermions. We find that by introducing supersymmetry the model can be made consistent with the long-standing deviation in the measured value of the anomalous magnetic moment of the muon, $g-2$, and neutralino dark matter of broad mass ranges and properties. In particular, dark matter candidates include the well-known 1 TeV higgsino, which in the MSSM is typically not compatible with solutions to the $g-2$ puzzle. Moreover, its spin-independent cross section could be at the origin of the recent small excess in XENON-1T data. We apply to the model constraints arising from flavor precision measurements and direct searches at the Large Hadron Collider and show that they do not currently exclude the relevant parameter space regions.

We review several current aspects of dark matter theory and experiment. We overview the present experimental status, which includes current bounds and recent claims and hints of a possible signal in a wide range of experiments: direct detection in underground laboratories, gamma-ray, cosmic ray, X-ray, neutrino telescopes, and the LHC. We briefly review several possible particle candidates for a Weakly Interactive Massive Particle (WIMP) and dark matter that have recently been considered in the literature. We pay particular attention to the lightest neutralino of supersymmetry as it remains the best motivated candidate for dark matter and also shows excellent detection prospects. Finally we briefly review some alternative scenarios that can considerably alter properties and prospects for the detection of dark matter obtained within the standard thermal WIMP paradigm.

We examine the question to what extent prospective detection of dark matter by direct and indirect- detection experiments could shed light on what fraction of dark matter was generated thermally via the freeze-out process in the early Universe. By simulating putative signals that could be seen in the near future and using them to reconstruct WIMP dark matter properties, we show that, in a model- independent approach this could only be achieved in a thin sliver of the parameter space. However, with additional theoretical input the hypothesis about the thermal freeze-out as the dominant mechanism for generating dark matter can potentially be verified. We illustrate this with two examples: an effective field theory of dark matter with a vector messenger and a higgsino or wino dark matter within the MSSM.

Using the existing simplified model framework, we build several dark matter models which have suppressed spin-independent scattering cross section. We show that the scattering cross section can vanish due to interference effects with models obtained by simple combinations of simplified models. For weakly interacting massive particle (WIMP) masses $\gtrsim$10 GeV, collider limits are usually much weaker than the direct detection limits coming from LUX or XENON100. However, for our model combinations, LHC analyses are more competitive for some parts of the parameter space. The regions with direct detection blind spots can be strongly constrained from the complementary use of several Large Hadron Collider (LHC) searches like mono-jet, jets + missing transverse energy, heavy vector resonance searches, etc. We evaluate the strongest limits for combinations of scalar + vector, "squark" + vector, and scalar + "squark" mediator, and present the LHC 14 TeV projections.

We examine the projected ability to reconstruct the mass, scattering, and annihilation cross section of dark matter in the new generation of large underground detectors, XENON-1T, SuperCDMS, and DarkSide-G2, in combination with diffuse gamma radiation from expected 15 years of data from Fermi-LAT observation of 46 local spiral dwarf galaxies and projected CTA sensitivity to a signal from the Galactic Center. To this end we consider several benchmark points spanning a wide range of WIMP mass, different annihilation final states, and large enough event rates to warrant detection in one or more experiments. As previously shown, below some 100 GeV only direct detection experiments will in principle be able to reconstruct the WIMP mass well. This may, in case a signal at Fermi-LAT is also detected, additionally help restricting \sigma v and the allowed decay branching rates. In the intermediate range between some 100 GeV and up a few hundred GeV, direct and indirect detection experiments can be used in complementarity to ameliorate the respective determinations, which in individual experiments can at best be rather poor, thus making the WIMP reconstruction in this mass range very challenging. At large WIMP mass, ~1 TeV, CTA will have the ability to reconstruct mass, annihilation cross section, and the allowed decay branching rates to very good precision for the $\tau^+\tau^-$ or purely leptonic final state, good for the $W^+W^-$ case, and rather poor for $b\bar{b}$. A substantial improvement can potentially be achieved by reducing the systematic uncertainties, increasing exposure, or by an additional measurement at Fermi-LAT that would help reconstruct the annihilation cross section and the allowed branching fractions to different final states.

We construct models of dark matter with suppressed spin-independent scattering cross section utilizing the existing simplified model framework. Even simple combinations of simplified models can exhibit interference effects that cause the tree level contribution to the scattering cross section to vanish, thus demonstrating that direct detection limits on simplified models are not robust when embedded in a more complicated and realistic framework. In general for fermionic WIMP masses >10 GeV direct detection limits on the spin-independent scattering cross section are much stronger than those coming from the LHC. However these model combinations, which we call less-simplified models, represent situations where LHC searches become more competitive than direct detection experiments even for moderate dark matter mass. We show that a complementary use of several searches at the LHC can strongly constrain the direct detection blind spots by setting limits on the coupling constants and mediators' mass. We derive the strongest limits for combinations of vector + scalar, vector + "squark", and "squark" + scalar mediator, and present the corresponding projections for the LHC 14 TeV for a number of searches: mono-jet, jets + missing energy, and searches for heavy vector resonances.

We examine axino dark matter in the regime of a low reheating temperature T_R after inflation and taking into account that reheating is a non-instantaneous process. This can have a significant effect on the dark matter abundance, mainly due to entropy production in inflaton decays. We study both thermal and non-thermal production of axinos in the context of the MSSM with ten free parameters. We identify the ranges of the axino mass and the reheating temperature allowed by the LHC and other particle physics data in different models of axino interactions. We confront these limits with cosmological constraints coming the observed dark matter density, large structures formation and big bang nucleosynthesis. We find a number of differences in the phenomenologically acceptable values of the axino mass and the reheating temperature relative to previous studies. In particular, an upper bound on the axino mass becomes dependent on T_R, reaching a maximum value at T_R~10^2 GeV. If the lightest ordinary supersymmetric particle is a wino or a higgsino, we obtain lower a limit of approximately 10 GeV for the reheating temperature. We demonstrate also that entropy production during reheating affects the maximum allowed axino mass and lowest values of the reheating temperature.

We present an updated analysis of the CMSSM and the NUHM using the latest experimental data and numerical tools. We map out favored regions of Bayesian posterior probability in light of data from the LHC, flavor observables, the relic density and dark matter searches. We present some updated features with respect to our previous analyses: we include the effects of corrections to the light Higgs mass beyond the 2-loop order using FeynHiggs v2.10.0; we include in the likelihood the latest limits from direct searches for squarks and gluinos at ATLAS with ~20/fb; the latest constraints on the spin-independent scattering cross section of the neutralino from LUX are applied taking into account uncertainties in the nuclear form factors. We find that in the CMSSM the posterior distribution now tends to favor smaller values of Msusy than in the previous analyses. As a consequence, the statistical weight of the A-resonance region increases to about 30% of the total probability, with interesting new prospects for the 14 TeV run at the LHC. The most favored region, on the other hand, still features multi-TeV squarks and gluinos, and ~1TeV higgsino dark matter whose detection prospects by current and one-tonne detectors look very promising. The same region is predominant in the NUHM, although the A-resonance region is also present there as well as a new solution, of neutralino-stau coannihilation through the channel stau stau -> hh at very large \mu. We derive the expected sensitivity of the future CTA experiment to ~1 TeV higgsino dark matter for both models and show that the prospects for probing both models are realistically good. We comment on the complementarity of this search to planned direct detection one-tonne experiments.

We investigate the prospects for detection of neutralino dark matter in the 19-parameter phenomenological MSSM (pMSSM). We explore very wide ranges of the pMSSM parameters but pay particular attention to the higgsino-like neutralino at the ~ 1 TeV scale, which has been shown to be a well motivated solution in many constrained supersymmetric models, as well as to a wino-dominated solution with the mass in the range of 2-3 TeV. After summarising the present bounds on the parameter space from direct and indirect detection experiments, we focus on prospects for detection of the Cherenkov Telescope Array (CTA). To this end, we derive a realistic assessment of the sensitivity of CTA to photon fluxes from dark matter annihilation by means of a binned likelihood analysis for the Einasto and Navarro-Frenk-White halo profiles. We use the most up to date instrument response functions and background simulation model provided by the CTA Collaboration. We find that, with 500 hours of observation, under the Einasto profile CTA is bound to exclude at the 95% C.L. almost all of the ~ 1 TeV higgsino region of the pMSSM, effectively closing the window for heavy supersymmetric dark matter in many realistic models. CTA will be able to probe the vast majority of cases corresponding to a spin-independent scattering cross section below the reach of 1-tonne underground detector searches for dark matter, in fact even well below the irreducible neutrino background for direct detection. On the other hand, many points lying beyond the sensitivity of CTA will be within the reach of 1-tonne detectors, and some within collider reach. Altogether, CTA will provide a highly sensitive way of searching for dark matter that will be partially overlapping and partially complementary with 1-tonne detector and collider searches, thus being instrumental to effectively explore the nearly full parameter space of the pMSSM.

Increasingly stringent limits from LHC searches for new physics, coupled with lack of convincing signals of weakly interacting massive particle (WIMP) in dark matter searches, have tightly constrained many realizations of the standard paradigm of thermally produced WIMPs as cold dark matter. In this article, we review more generally both thermally and non-thermally produced dark matter (DM). One may classify DM models into two broad categories: one involving bosonic coherent motion (BCM) and the other involving WIMPs. BCM and WIMP candidates need, respectively, some approximate global symmetries and almost exact discrete symmetries. Supersymmetric axion models are highly motivated since they emerge from compelling and elegant solutions to the two fine-tuning problems of the Standard Model: the strong CP problem and the gauge hierarchy problem. We review here non-thermal relics in a general setup, but we also pay particular attention to the rich cosmological properties of various aspects of mixed SUSY/axion dark matter candidates which can involve both WIMPs and BCM in an interwoven manner. We also review briefly a panoply of alternative thermal and non-thermal DM candidates.

We examine a scenario in which the reheating temperature $T_R$ after inflation is so low that it is comparable to, or lower than, the freeze out temperature of ordinary WIMPs. In this case the dark matter relic abundance is reduced, thus relaxing the impact of the usually strong constraint coming from the requirement that the universe does not overclose. We first re-examine the dynamics of freezeout during reheating. Next we apply a Bayesian approach to study the parameter space of the MSSM with ten free parameters, the CMSSM and the singlino-dominated regions of the NMSSM. In each case we find dramatic departures from the usually considered regime of high $T_R$, with important implications for direct detection dark matter searches. In the MSSM we examine WIMP mass range up to ~5 TeV, and find regions of bino dark matter over the whole mass range, and of higgsino dark matter with mass over a similar range but starting from the ~1 TeV value of the standard high $T_R$ scenario. We show that the prospects for bino detection strongly depend on $T_R$, while the higgsino is for the most part detectable by future one-tonne detectors. The wino, which is excluded in the standard scenario, becomes allowed again if its mass is roughly above 3.5 TeV, and can be detectable. In the CMSSM, the bino and higgsino mass ranges become more constrained although detection prospects remain similar. In the Next-to-MSSM at low enough $T_R$ wide ranges of singlino-dominated parameter space of the model become cosmologically allowed. We also study the contribution to the DM relic density from direct and cascade decays of the inflaton. Finally, we consider the case of a gravitino as dark matter. We find strong bounds from overclosure and Big Bang Nucleosynthesis, and derive lower limits on $T_R$ which depend on the gravitino mass and on the nature of the lightest ordinary superpartner.

We explore the MSSM with 9 free parameters (p9MSSM) that have been selected as a minimum set that allows an investigation of neutralino dark matter and collider signatures while maintaining consistency with several constraints. These include measurement of the dark matter relic density from PLANCK, main properties of the discovered Higgs boson, LHC direct SUSY searches, recent evidence for a Standard Model-like BR(Bs -> \mu+ \mu-), and the measurement of \delta(g-2), plus a number of other electroweak and flavor physics constraints. We perform a simulation of two LHC direct SUSY searches at sqrt(s)=8 TeV: the CMS inclusive \alpha_T search for squarks and gluinos and the CMS electroweak production search with 3l+E_T^miss in the final state. We use the latter to identify the regions of the parameter space, consistent at 2\sigma\ with \delta(g-2), that are not excluded by the direct limits from the electroweak production. We find that they correspond to a neutralino mass in the window 200 GeV<m_\chi<500 GeV. We also implement the likelihood for the XENON100 exclusion bound, in which we consider for the first time the impact of a recent determination of the \Sigma_{\pi N} term from CHAOS data, \Sigma_{\pi N}=43\pm12 MeV. We show that in light of this measurement, the present statistical impact of the XENON100 bound is greatly reduced, although future sensitivities of the LUX and XENON1T experiments will have decisive impact on the mixed bino/higgsino composition of the neutralino. We point out some tension between the constraints from \delta(g-2) and XENON100. Finally, we present prospects for various indirect searches of dark matter, namely \gamma-ray fluxes from dSphs and the Galactic Center at Fermi-LAT, and the positron flux at AMS02. We also show the 5-year sensitivity on the spin-dependent neutralino-proton cross section due to neutrino fluxes from the Sun at IceCube.

We review the status of axino dark matter. Two hierarchy problems, the strong CP problem and the gauge hierarchy problem, have led to introducing into particle physics a spontaneously broken global Peccei-Quinn symmetry and a softly broken supersymmetry, respectively. Combining them implies the presence of not only an axion, but also of its scalar component, saxion, and their fermionic partner, axino. Among these, the axion and the axino are attractive dark matter candidates. Various possibilities for the axino as dark matter are discussed.

We investigate the impact of recent limits from LHC searches for supersymmetry and from direct and indirect searches for dark matter on global Bayesian inferences of the parameter space of the Constrained Minimal Supersymmetric Standard Model (CMSSM). In particular we apply recent exclusion limits from the CMS \alpha_T analysis of 1.1/fb of integrated luminosity, the current direct detection dark matter limit from XENON100, as well as recent experimental constraints on \gamma-ray fluxes from dwarf spheroidal satellite galaxies of the Milky Way from the FermiLAT telescope, in addition to updating values for other non-LHC experimental constraints. We extend the range of scanned parameters to include a significant fraction of the focus point/hyperbolic branch region. While we confirm earlier conclusions that at present LHC limits provide the strongest constraints on the model's parameters, we also find that when the uncertainties are not treated in an excessively conservative way, the new bounds from dwarf spheroidal have the power to significantly constrain the focus point/hyperbolic branch region. Their effect is then comparable, if not stronger, to that from XENON100. We further analyze the effects of one-year projected sensitivities on the neutrino flux from the Sun in the 86-string IceCube+DeepCore configuration at the South Pole. We show that data on neutrinos from the Sun, expected for the next few months at IceCube and DeepCore, have the potential to further constrain the same region of parameter space independently of the LHC and can yield additional investigating power for the model.

Axino arises in supersymmetric versions of axion models and is a natural candidate for cold or warm dark matter. Here we revisit axino dark matter produced thermally and non-thermally in light of recent developments. First we discuss the definition of axino relative to low energy axion one for several KSVZ and DFSZ models of the axion. Then we review and refine the computation of the dominant QCD production in order to avoid unphysical cross-sections and, depending on the model, to include production via SU(2) and U(1) interactions and Yukawa couplings.

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 extent to which the uncertainties associated with the propagation of Galactic cosmic rays impact upon estimates for the gamma-ray flux from the mid-latitude region. We consider contributions from both standard astrophysical background (SAB) processes as well as resolved point sources. We have found that the uncertainties in the total gamma-ray flux from the mid-latitude region relating to propagation parameter values consistent with local B/C and Be10/Be9 data dominate by 1-2 orders of magnitude. These uncertainties are reduced to less than an order of magnitude when the normalisations of the SAB spectral components are fitted to the corresponding Fermi LAT data. We have found that for many propagation parameter configurations (PPCs) our fits improve when an extragalactic background (EGB) component is simultaneously fitted to the data. We also investigate the improvement in our fits when a flux contribution from neutralino dark matter (DM), described by the Minimal Supersymmetric Standard Model, was simultaneously fitted to the data. We consider three representative cases of neutralino DM for both Burkert and Einasto DM density profiles, in each case simultaneously fitting a boost factor of the DM contribution together with the SAB and EGB components. We have found that for several PPCs there are significant improvements in our fits, yielding both substantial EGB and DM components, where for a few of these PPCs the best-fit EGB component is consistent with recent estimates by the Fermi Collaboration.

We re-consider the gravitino as dark matter in the framework of the Constrained MSSM. We include several recently suggested improvements on: (i) the thermal production of gravitinos, (ii) the calculation of the hadronic spectrum from NLSP decay and (iii) the BBN calculation including stau bound-state effects. In most cases we find an upper bound on the reheating temperature $\treh\lsim {\rm a few}\times 10^7\gev$ from over-production of $^6{\rm Li}$ from bound state effects. We also find an upper limit on the stau lifetime of $3\times10^4 \sec$, which is nearly an order of magnitude larger than the simple limit $5\times 10^3 \sec$ often used to avoid the effect of bound-state catalysis. The bound on $\treh$ is relaxed to $\lsim {\rm a few}\times 10^8\gev$ when we use a more conservative bound on ${^6}{\rm Li}/{^7}{\rm Li}$, in which case a new region at small stau mass at $\sim 100\gev$ and much longer lifetimes opens up. Such a low stau mass region can be easily tested at the LHC.

The diffuse Galactic gamma-ray data from the region of the Galactic Center has been collected by the LAT instrument on the Fermi Gamma-Ray Space Telescope. In this paper we argue that it may be able to provide an unambiguous evidence of originating, in addition to known astrophysical sources, from dark matter annihilations in the halo, independently of the mass and other properties of the dark matter particle. We also show that the recently released high precision data from mid-latitudes is already providing an upper bound, albeit still a weak one, on the cuspiness of the dark matter density profile as a function of the mass of the dark matter assumed to be a stable neutralino of minimal supersymmetry.

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.

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 show that the $\Omega_b - \Omega_{DM}$ coincidence can naturally be explained in a framework where axino cold dark matter is predominantly produced in non-thermal processes involving decays of Q-balls formed in Affleck-Dinebaryogenesis. In this approach the similarity between $\Omega_b$ and $\Omega_{DM}$ is a direct consequence of the (sub-)GeV scale of the mass of the axino.

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.

The axino and the gravitino are well-motivated candidates for the lightest supersymmetric particle (LSP) and also for cold dark matter in the Universe. Assuming that a charged slepton is the next-to-lightest supersymmetric particle (NLSP), we show how the NLSP decays can be used to probe the axino LSP scenario in hadronic axion models as well as the gravitino LSP scenario at the Large Hadron Collider and the International Linear Collider. We show how one can identify experimentally the scenario realized in nature. In the case of the axino LSP, the NLSP decays will allow one to estimate the value of the axino mass and the Peccei-Quinn scale.

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.

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 investigate the role of an effective dimension-4 axino-quark-squark coupling in the thermal processes producing stable cold axino relics in the early Universe. We find that, while the induced squark and quark scattering processes are always negligible, squark decays become important in the case of low reheat temperature and large gluino mass. The effect can tighten the bounds on the scenario from the requirement that cold dark matter axinos do not overclose the Universe.

Supersymmetric extensions of the Standard Model that incorporate the axion solution to the strong CP problem necessarily contain also the axino, the fermionic partner of the axion. In contrast to the neutralino and the gravitino, the axino mass is generically not of the order of the supersymmetry-breaking scale and can be much smaller. The axino is therefore an intriguing candidate for a stable superpartner. In a previous Letter [1] it was shown that axinos are a natural candidate for cold dark matter in the Universe when they are generated non-thermally through out-of-equilibrium neutralino decays. Here, we extend the study of non-thermal production and include a competing thermal production mechanism through scatterings and decays of particles in the plasma. We identify axino masses in the range of tens of MeV to several GeV (depending on the scenario) as corresponding to cold axino relics if the reheating temperature $\treh$ is less than about $5\times10^4\gev$. At higher $\treh$ and lower mass, axinos could constitute warm dark matter. In the scenario with axinos as relics, the gravitino problem finds a natural solution. The lightest superpartner of the Standard Model spectrum remains effectively stable in high-energy detectors but may be either neutral or charged. The usual highly restrictive constraint $\abundchi\lsim1$ on the relic abundance of the lightest neutralino becomes void.

As experimental sensitivity increases, one is approaching the range of WIMP-nucleon interaction strengths characteristic of neutralinos. But this continuing progress also unearths new experimental challenges and uncertainties.

There exist several well-motivated candidates for non-baryonic cold dark matter, including neutralinos, axions, axinos, gravitinos, Wimpzillas. I review the dark matter properties of the neutralino and the current status of its detection. I also discuss the axino as a new interesting alternative.

Supersymmetric extensions of the Standard Model when combined with the Peccei-Quinn solution to the strong CP problem necessarily contain also the axino, the fermionic partner of the axion. In contrast to the neutralino and the gravitino, the axino mass is generically not of the order of the supersymmetry breaking scale and can be much smaller. The axino is therefore an intriguing candidate for a stable superpartner. The axinos are a natural candidate for cold dark matter in the Universe when they are generated non-thermally through out-of-equilibrium neutralino decays or via a competing thermal production mechanism through scatterings and decays of particles in the plasma. We identify axino masses in the range of tens of MeV to several GeV (depending on the scenario) as corresponding to cold axino relics if the reheating temperature T_R is less than about 5 \times 10^4 GeV. At higher T_R and lower mass, axinos could constitute warm dark matter. In the scenario with axinos as stable relics the gravitino problem finds a natural solution. The lightest superpartner of the Standard Model spectrum will remain stable in high-energy detectors but may be either neutral or charged. The usual constraint Omega*h^2<1 on the relic abundance of the lightest neutralino does not hold.

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.

We show that axinos produced in the early Universe in the decay of the lightest neutralinos are a natural candidate for cold dark matter. We argue that axinos may well provide the main component of the missing mass in the Universe because their relic density is often comparable with the critical density.

I review axions, neutralinos, axinos, gravitinos and super-massive Wimpzillas as dark matter candidates.

We examine the effect of some uncertainties in the input astrophysical parameters on direct detection searches for WIMPs in the Galactic halo. We concentrate on the possible WIMP annual modulation signal recently reported by the DAMA Collaboration. We find that allowing for a reasonable uncertainty in a WIMP Maxwellian velocity distribution leads to significantly relaxed constraints on the WIMP mass as compared to the original DAMA analysis.

We explore a model of inflation where the inflaton mass-squared is generated at a high scale by gravity-mediated soft supersymmetry breaking, and runs at lower scales to the small value required for slow-roll inflation. The running is supposed to come from the coupling of the inflaton to a non-Abelian gauge field. In contrast with earlier work, we do not constrain the magnitude of the supersymmetry breaking scale, and we find that the model might work even if squark and slepton masses come from gauge-mediated supersymmetry breaking. With the inflaton and gaugino masses in the expected range, and $\alpha = g^2/4\pi $ in the range $10^{-2}$ to $10^{-3}$ (all at the high scale) the model can give the observed cosmic microwave anisotropy, and a spectral index in the observed range. The latter has significant variation with scale, which can confirm or rule out the model in the forseeable future.

In constrained minimal supersymmetry the lightest neutralino of bino-type is the only neutral candidate for dark matter. As a result, one is typically able to restrict all the supersymmetric mass spectra below roughly 1-2\tev\ {\em without} imposing an ill-defined fine-tuning constraint.

In the constrained MSSM one is typically able to restrict the supersymmetric mass spectra below roughly 1-2\tev\ {\em without} resorting to the ambiguous fine-tuning constraint.

A simple procedure is presented which leads to a dramatic simplification in the calculation of the relic density of stable particles in the Universe.

I address the question of whether supersymmetry provides a viable candidate for the dark matter in the Universe. I review the properties of the lightest neutralino as a candidate for solving the dark matter problem. I discuss the neutralino's phenomenological and cosmological properties, and constraints from present and future experiments. In the minimal supersymmetric model, the neutralino mass has been experimentally excluded below some $20$ GeV, and is not expected to be significantly larger than about $150$ GeV. I identify a gaugino-like neutralino as the most natural dark matter candidate for a plausible range of parameters. The requirement that the lightest neutralino be the dominant matter component in the flat Universe provides non-trivial restrictions on other parameters of the model, in particular on the masses of the sfermions. Next, I study the consequences of adopting further grand unification assumptions. In both scenarios I find sfermion masses most likely beyond the reach of LEP 200 and the Tevatron but well within the discovery potential of the SSC and the LHC. I also comment on the effects of relaxing grand unification assumptions. Finally, I briefly outline prospects for the neutralino dark matter searches.

A recent discovery by the COBE satellite is presented, and its importance discussed in the context of the physics of the early Universe. The various implications for our understanding of the structure of the Universe are outlined.