Rather general considerations of the string theory landscape imply a mild statistical draw towards large soft SUSY breaking terms tempered by the requirement of proper electroweak symmetry breaking where SUSY contributions to the weak scale are not too far from m(weak)~ 100 GeV. Such a picture leads to the prediction that m_h~ 125 GeV while most sparticles are beyond current LHC reach. Here we explore the possibility that the magnitude of the Peccei-Quinn (PQ) scale f_a is also set by string landscape considerations within the framework of a compelling SUSY axion model. First, we examine the case where the PQ symmetry arises as an accidental approximate global symmetry from a more fundamental gravity-safe Z(24)^R symmetry and where the SUSY mu parameter arises from a Kim-Nilles operator. The pull towards large soft terms then also pulls the PQ scale as large as possible. Unless this is tempered by rather severe (unknown) cosmological or anthropic bounds on the density of dark matter, then we would expect a far greater abundance of dark matter than is observed. This conclusion cannot be negated by adopting a tiny axion misalignment angle theta_i because WIMPs are also overproduced at large f_a. Hence, we conclude that setting the PQ scale via anthropics is highly unlikely. Instead, requiring soft SUSY breaking terms of order the gravity-mediation scale m_{3/2}~ 10-100 TeV places the mixed axion-neutralino dark matter abundance into the intermediate scale sweet zone where f_a~ 10^{11}-10^{12} GeV. We compare our analysis to the more general case of a generic SUSY DFSZ axion model with uniform selection on theta_i but leading to the measured dark matter abundance: this approach leads to a preference for f_a~ 10^{12} GeV.

The string theory landscape of vacua solutions provides physicists with some understanding as to the magnitude of the cosmological constant. Similar reasoning can be applied to the magnitude of the soft SUSY breaking terms in supersymmetric models of particle physics: there appears to be a statistical draw towards large soft terms which is tempered by the anthropic requirement of the weak scale lying not too far from ~100 GeV. For a mild statistical draw of m_{soft}^n with n=1 (as expected from SUSY breaking due to a single F term) then the light Higgs mass is preferred at ~125 GeV while sparticles are all pulled beyond LHC bounds. We confront a variety of LHC and WIMP dark matter search limits with the statistical expectations from a fertile patch of string theory landscape. The end result is that LHC and WIMP dark matter detectors see exactly that which is expected from the string theory landscape: a Standard Model-like Higgs boson of mass 125 GeV but as yet no sign of sparticles or WIMP dark matter. SUSY from the n=1 landscape is most likely to emerge at LHC in the soft opposite-sign dilepton plus jet plus MET channel. Multi-ton noble liquid WIMP detectors should be able to completely explore the n=1 landscape parameter space.

Particle physics models with Peccei-Quinn (PQ) symmetry breaking as a consequence of supersymmetry (SUSY) breaking are attractive in that they solve the strong CP problem with a SUSY DFSZ-like axion, link the SUSY breaking and PQ breaking intermediate mass scales and can resolve the SUSY mu problem with a naturalness-required weak scale mu term whilst soft SUSY breaking terms inhabit the multi-TeV regime as required by LHC sparticle mass limits and the Higgs mass measurement. On the negative ledger, models based on global symmetries suffer a generic gravity spoliation problem. We present a model based on the discrete R-symmetry Z_{24}^R-- which may emerge from compactification of 10-d Lorentzian spacetime in string theory-- where the mu term and dangerous proton decay and R-parity violating operators are either suppressed or forbidden while a gravity-safe PQ symmetry emerges as an accidental approximate global symmetry leading to a solution to the strong CP problem and a weak-scale/natural value for the mu term.

The requirement of electroweak naturalness in supersymmetric (SUSY) models of particle physics necessitates light higgsinos not too far from the weak scale characterized by m(weak)~ m(W,Z,h)~100 GeV. On the other hand, LHC Higgs mass measurements and sparticle mass limits point to a SUSY breaking scale in the multi-TeV regime. Under such conditions, the lightest SUSY particle is expected to be a mainly higgsino-like neutralino with non-negligible gaugino components (required by naturalness). The computed thermal WIMP abundance in natural SUSY models is then found to be typically a factor 5-20 below its measured value. To gain concordance with observations, either an additional DM particle (the axion is a well-motivated possibility) must be present or additional non-thermal mechanisms must augment the neutralino abundance. We compare present direct and indirect WIMP detection limits to three natural SUSY models based on gravity-, anomaly- and mirage-mediation. We show that the case of natural higgsino-only dark matter where non-thermal production mechanisms augment its relic density, is essentially excluded by a combination of direct detection constraints from PandaX-II, LUX and Xenon-1t experiments, and by bounds from Fermi-LAT/MAGIC observations of gamma rays from dwarf spheroidal galaxies.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

Under the expectation that nature is natural, we extend the Standard Model to include SUSY to stabilize the electroweak sector and PQ symmetry to stabilize the QCD sector. Then natural SUSY arises from a Kim-Nilles solution to the SUSY mu problem which allows for a little hierarchy where mu ~ f_a^2/M_P ~ 100-300 GeV while the SUSY particle mass scale m(SUSY)~ 1-10 TeV >> mu. Dark matter then consists of two particles: a higgsino-like WIMP and a SUSY DFSZ axion. The range of allowed axion mass values m(axion) depends on the mixed axion-higgsino relic density. The range of m(axion) is actually restricted in this case by limits on WIMPs from direct and indirect detection experiments. We plot the expected axion detection rate at microwave cavity experiments. The axion-photon-photon coupling is severely diminished by charged higgsino contributions to the anomalous coupling. In this case, the axion may retreat, at least temporarily, back into the regime of near invisibility. From our results, we urge new ideas for techniques which probe both deeper and more broadly into axion coupling versus axion mass parameter space.

We examine updated prospects for detecting WIMPs in supersymmetric models via direct and indirect dark matter search experiments. We examine several historical and also still viable scenarios: projections for well-tempered neutralinos (WTN), projections from the MasterCode (MC), BayesFits (BF) and Fittino (FO) collaborations, non-thermal wino dark matter (NThW) and finally mixed axion-higgsino dark matter from SUSY with radiatively-driven naturalness (RNS). The WTN is ruled out by recent limits from XENON and LUX collaborations. The NThW scenario, previously on tenuous ground due to gamma-line searches, appears also ruled out by recent combined Fermi-LAT/MAGIC limits combined with new HESS results from continuum gamma rays. Substantial portions of MC parameter space and 1 TeV higgsino parameter space from BF group are ruled out. The 100-300 GeV higgsino-like WIMP from RNS survives due to its possible depleted local abundance (where the axion may make up the bulk of dark matter). Projections from ton-scale noble liquid detectors should discover or rule out WIMPs from the remaining parameter space of these surviving models.

The strong CP problem of QCD is at heart a problem of naturalness: why is the F\tilde{F} term highly suppressed in the QCD Lagrangian when it seems necessary to explain why there are three and not four light pions? The most elegant solution posits a spontaneously broken Peccei-Quinn (PQ) symmetry which requires the existence of the axion field a. The axion field settles to the minimum of its potential thus removing the offensive term but giving rise to the physical axion whose coherent oscillations can make up the cold dark matter. Only now are experiments such as ADMX beginning to explore QCD axion parameter space. Since a bonafide scalar particle-- the Higgs boson-- has been discovered, one might expect its mass to reside at the axion scale f_a~ 10^{11} GeV. The Higgs mass is elegantly stabilized by supersymmetry: in this case the axion is accompanied by its axino and saxion superpartners. Requiring naturalness also in the electroweak sector implies higgsino-like WIMPs so then we expect mixed axion-WIMP dark matter. Ultimately we would expect detection of both an axion and a WIMP while signals for light higgsinos may show up at LHC and must show up at ILC.

Supersymmetric models with radiatively-driven electroweak naturalness require light higgsinos of mass ~ 100-300 GeV. Naturalness in the QCD sector is invoked via the Peccei-Quinn (PQ) axion leading to mixed axion-higgsino dark matter. The SUSY DFSZ axion model provides a solution to the SUSY mu problem and the Little Hierarchy \mu << m_{3/2} may emerge as a consequence of a mismatch between PQ and hidden sector mass scales. The traditional gravitino problem is now augmented by the axino and saxion problems, since these latter particles can also contribute to overproduction of WIMPs or dark radiation, or violation of BBN constraints. We compute regions of the T_R vs. m_{3/2} plane allowed by BBN, dark matter and dark radiation constraints for various PQ scale choices f_a. These regions are compared to the values needed for thermal leptogenesis, non-thermal leptogenesis, oscillating sneutrino leptogenesis and Affleck-Dine leptogenesis. The latter three are allowed in wide regions of parameter space for PQ scale f_a~ 10^{10}-10^{12} GeV which is also favored by naturalness: f_a~ \sqrt{\mu M_P/\lambda_\mu }\sim 10^{10}-10^{12} GeV. These f_a values correspond to axion masses somewhat above the projected ADMX search regions.

We examine mixed axion/neutralino cold dark matter production in the SUSY DFSZ axion model where an axion superfield couples to Higgs superfields. We calculate a wide array of axino and saxion decay modes along with their decay temperatures, and thermal and non-thermal production rates. For a SUSY benchmark model with a standard underabundance (SUA) of Higgsino-like dark matter (DM), we find for the PQ scale f_a< 10^{12} GeV that the DM abundance is mainly comprised of axions as the saxion/axino decay occurs before the standard neutralino freeze-out and thus its abundance remains suppressed. For 10^{12}<~ f_a<~ 10^{14} GeV, the saxion/axino decays occur after neutralino freeze-out so that the neutralino abundance is enhanced by the production via decay and subsequent re-annihilation. For f_a>~ 10^{14} GeV, both neutralino dark matter and dark radiation are typically overproduced. For judicious parameter choices, these can be suppressed and the combined neutralino/axion abundance brought into accord with measured values. A SUSY benchmark model with a standard overabundance (SOA) of bino DM is also examined and typically remains excluded due at least to too great a neutralino DM abundance for f_a<~ 10^{15} GeV. For f_a>~ 10^{15} GeV and lower saxion masses, large entropy production from saxion decay can dilute all relics and the SOA model can be allowed by all constraints.

By insisting on naturalness in both the electroweak and QCD sectors of the MSSM, the portrait for dark matter production is seriously modified from the usual WIMP miracle picture. In SUSY models with radiatively-driven naturalness (radiative natural SUSY or RNS) which include a DFSZ-like solution to the strong CP and SUSY mu problems, dark matter is expected to be an admixture of both axions and higgsino-like WIMPs. The WIMP/axion abundance calculation requires simultaneous solution of a set of coupled Boltzmann equations which describe quasi-stable axinos and saxions. In most of parameter space, axions make up the dominant contribution of dark matter although regions of WIMP dominance also occur. We show the allowed range of PQ scale f_a and compare to the values expected to be probed by the ADMX axion detector in the near future. We also show WIMP detection rates which are suppressed from usual expectations because now WIMPs comprise only a fraction of the total dark matter. Nonetheless, ton-scale noble liquid detectors should be able to probe the entirety of RNS parameter space. Indirect WIMP detection rates are less propitious since they are reduced by the square of the depleted WIMP abundance.

By eschewing finetuning from the electroweak and QCD sectors of supersymmetry (natural supersymmetry or SUSY), and by invoking the Kim-Nilles solution to the SUSY mu problem, one is lead to models wherein the dark matter is comprised of a mixture of axions and higgsino-like WIMPs. Over a large range of Peccei-Quinn breaking scale f_a~ 10^9-10^{12} GeV, one then expects about 90-95% axion dark matter. In such a scenario, both axion and WIMP direct detection may be expected.

The supersymmetrized DFSZ axion model is highly motivated not only because it offers solutions to both the gauge hierarchy and strong CP problems, but also because it provides a solution to the SUSY mu problem which naturally allows for a Little Hierarchy. We compute the expected mixed axion-neutralino dark matter abundance for the SUSY DFSZ axion model in two benchmark cases-- a natural SUSY model with a standard neutralino underabundance (SUA) and an mSUGRA/CMSSM model with a standard overabundance (SOA). Our computation implements coupled Boltzmann equations which track the radiation density along with neutralino, axion (produced thermally (TH) and via coherent oscillations (CO)), saxion (TH- and CO-produced), axino and gravitino densities. In the SUSY DFSZ model, axions, axinos and saxions go through the process of freeze-in-- in contrast to freeze-out or out-of-equilibrium production as in the SUSY KSVZ model-- resulting in thermal yields which are largely independent of the re-heat temperature. We find the SUA case with suppressed saxion-axion couplings (\xi=0) only admits solutions for PQ breaking scale f_a~< 5\times 10^{12} GeV where the bulk of parameter space tends to be axion-dominated. For SUA with allowed saxion-axion couplings (\xi =1), then f_a values up to ~ 2\times 10^{14} GeV are allowed. For the SOA case, almost all of SUSY DFSZ parameter space is disallowed by a combination of overproduction of dark matter, overproduction of dark radiation or violation of BBN constraints. An exception occurs at very large f_a~ 10^{15}-10^{16} GeV where large entropy dilution from CO-produced saxions leads to allowed models.

A variety of supersymmetric models give rise to a split mass spectrum characterized by very heavy scalars but sub-TeV gauginos, usually with a wino-like LSP. Such models predict a thermally-produced underabundance of wino-like WIMP dark matter so that non-thermal DM production mechanisms are necessary. We examine the case where theories with a wino-like LSP are augmented by a Peccei-Quinn sector including an axion-axino-saxion supermultiplet in either the SUSY KSVZ or SUSY DFSZ models and with/without saxion decays to axions/axinos. We show allowed ranges of PQ breaking scale f_a for various cases which are generated by solving the necessary coupled Boltzmann equations. We also present results for a model with radiatively-driven naturalness but with a wino-like LSP.

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.

Recent null results from LHC8 SUSY searches along with the discovery of a SM-like Higgs boson with mass m(h) ~ 125.5 GeV indicates sparticle masses in the TeV range, causing tension with conventional measures of electroweak fine-tuning. We propose a simple Fine-tuning Rule which should be followed under any credible evaluation of fine-tuning. We believe that overestimates of electroweak fine-tuning by conventional measures all arise from violations of this rule. We show that to gain accord with the Fine-tuning Rule, then both Higgs mass and the traditional \Delta_{BG} fine-tuning measures reduce to the model-independent electroweak fine-tuning measure \Delta_{EW}. This occurs by combining dependent contributions to m(Z) or m(h) into independent units. Then, using \Delta_{EW}, we evaluate EW fine-tuning for a variety of SUSY models including mSUGRA, NUHM1, NUHM2, mGMSB, mAMSB, hyper-charged AMSB and nine cases of mixed moduli-anomaly (mirage) mediated SUSY breaking models (MMAMSB) whilst respecting LHC Higgs mass and B-decay constraints (we do not impose LHC8 sparticle mass constraints due to the possibility of compressed spectra within many of these models). We find mSUGRA, mGMSB, mAMSB and MMAMSB models all to be highly fine-tuned. The NUHM1 model is moderately fine-tuned while NUHM2 which allows for radiatively-driven naturalness (RNS) allows for fine-tuning at a meager 10% level in the case where m(higgsinos) ~ 100-200 GeV and the TeV-scale top squarks are well-mixed. Models with RNS may or may not be detect at LHC14. A \sqrt{s} ~ 500 GeV e^+e^- collider will be required to make a definitive search for the requisite light higgsinos.

In this Report we discuss the four complementary searches for the identity of dark matter: direct detection experiments that look for dark matter interacting in the lab, indirect detection experiments that connect lab signals to dark matter in our own and other galaxies, collider experiments that elucidate the particle properties of dark matter, and astrophysical probes sensitive to non-gravitational interactions of dark matter. The complementarity among the different dark matter searches is discussed qualitatively and illustrated quantitatively in several theoretical scenarios. Our primary conclusion is that the diversity of possible dark matter candidates requires a balanced program based on all four of those approaches.

In the post-LHC8 era, it is perceived that what is left of SUSY model parameter space is highly finetuned in the EW sector (EWFT). We discuss how conventional measures overestimate EWFT in SUSY theory. Radiatively-driven natural SUSY (RNS) models maintain the SUSY GUT paradigm with low EWFT at 10% level, but are characterized by light higgsinos ~100-300 GeV and a thermal underabundance of WIMP dark matter. Implementing the SUSY DFSZ solution to the strong CP problem explains the small \mu parameter but indicates dark matter should be comprised mainly of axions with a small admixture of higgsino-like WIMPs. While RNS might escape LHC14 searches, we would expect ultimately direct detection of both WIMPs and axions. An e^+e^- collider with \sqrt{s} ~ 500-600 GeV should provide a thorough search for the predicted light higgsinos.

Supersymmetric models with low electroweak fine-tuning contain light higgsinos with mass not too far from m_h\simeq 125 GeV, while other sparticles can be much heavier. In the R-parity conserving MSSM, the lightest neutralino is then a higgsino-like WIMP (albeit with non-negligible gaugino components), with thermal relic density well below measured values. This leaves room for axions (or other, perhaps not as well motivated, stable particles) to function as co-dark matter particles. The local WIMP abundance is then expected to be below standard estimates, and direct and indirect detection rates must be accordingly rescaled. We calculate rescaled direct and indirect higgsino-like WIMP detection rates in SUSY models that fulfil the electroweak naturalness condition. In spite of the rescaling, we find that ton-scale noble liquid detectors can probe the entire higgsino-like WIMP parameter space, so that these experiments should either discover WIMPs or exclude the concept of electroweak naturalness in R-parity conserving natural SUSY models. Prospects for spin-dependent or indirect detection are more limited due in part to the rescaling effect.

If one simultaneously invokes the SUSY solution to the gauge hierarchy problem and the PQ solution to the strong CP problem, then one might expect mixed axion/neutralino dark matter (DM), i.e. two dark matter particles. In this case, dark matter production involves oscillation production of axions, thermal and non-thermal neutralino production and thermal and non-thermal production of axinos and saxions. The latter particles may feed additional neutralinos into the cosmic soup, or dilute all relics via additional entropy production, or feed dark radiation into the cosmos. Depending on the SUSY and PQ model, the PQ scale can easily exist from 10^9-10^{16} GeV. Axions may lie well beyond the range currently being probed by experiment. In the case of relic WIMPs, theoretical predictions must be rescaled for the possibility that WIMPs only comprise a fraction of the local DM abundance.

Supersymmetric models which fulfill the conditions of electroweak naturalness generally contain light higgsinos with mass not too far from m_h ~125 GeV, while other sparticles can be much heavier. In R-parity conserving models, the lightest neutralino is then a higgsino-like WIMP (albeit with non-negligible gaugino components), with thermal relic density well below measured values. This leaves room for axions to function as co-dark matter particles. The local WIMP abundance is then expected to be below standard estimates, and direct and indirect detection rates must be accordingly rescaled. We calculate rescaled direct and indirect higgsino-like WIMP detection rates in SUSY models that fulfil the electroweak naturalness condition. In spite of the rescaling, we find that ton-scale noble liquid detectors can probe the entire higgsino-like WIMP parameter space, so that these experiments should either discover WIMPs or exclude the concept of electroweak naturalness in R-parity conserving natural SUSY models. Prospects for spin-dependent or indirect detection are more limited due in part to the rescaling effect.

Recent analyses of WMAP9 data show that dark radiation-- parametrized by the apparent number of additional neutrinos \Delta N_{eff} contributing to the cosmic expansion-- is consistent with the Standard Model and bounded from above by about \Delta N_{eff} ~< 0.5 at 95% CL. We consider the mixed axion/neutralino cold dark matter scenario which arises in R-parity conserving supersymmetric (SUSY) models wherein the strong CP problem is solved by hadronic axions with a concommitant axion(a)/saxion(s)/axino(\ta) supermultiplet. Our new results include improved calculations of thermal axion and saxion production and include effects of saxion decay to axinos and axions. We show that the above bound on \Delta N_{eff} is easily satisfied if saxions are mainly thermally produced and m_{LSP} < m_{\ta} \lesssim m_s. However, if the dominant mechanism of saxion production is through coherent oscillations, the WMAP9 data provides a strong bound on saxion production followed by saxion decays to axions. Furthermore we show that scenarios with mixed neutralino/axion dark matter are highly constrained by combined WMAP9, BBN and Xe-100 constraints. In particular, supersymmetric models with a standard overabundance of neutralino dark matter are excluded for all values of the Peccei-Quinn breaking scale. Next generation WIMP direct detection experiments may be able to discover or exclude mixed axion-neutralino CDM scenarios where $s\to aa$ is the dominant saxion decay mode.

Models of natural supersymmetry seek to solve the little hierarchy problem by positing a spectrum of light higgsinos \lesssim 200 GeV and light top squarks \lesssim 500 GeV along with very heavy squarks and TeV-scale gluinos. Such models have low electroweak finetuning and are safe from LHC searches. However, in the context of the MSSM, they predict too low a value of m(h) and the relic density of thermally produced higgsino-like WIMPs falls well below dark matter (DM) measurements. Allowing for high scale soft SUSY breaking Higgs mass m_{H_u}> m_0 leads to natural cancellations during RG running, and to radiatively induced low finetuning at the electroweak scale. This model of radiative natural SUSY (RNS), with large mixing in the top squark sector, allows for finetuning at the 5-10% level with TeV-scale top squarks and a 125 GeV light Higgs scalar h. If the strong CP problem is solved via the PQ mechanism, then we expect an axion-higgsino admixture of dark matter, where either or both the DM particles might be directly detected.

The ATLAS and CMS collaborations have reported an excess of events in the \gamma\gamma, ZZ^*\to 4\ell and WW^* search channels at an invariant mass m \simeq 125 GeV, which could be the first evidence for the long-awaited Higgs boson. We investigate the consequences of requiring m_h\simeq 125 GeV in both the mSUGRA and NUHM2 SUSY models. In mSUGRA, large values of trilinear soft breaking parameter |A_0| are required, and universal scalar m_0\agt 0.8 TeV is favored so that we expect squark and slepton masses typically in the multi-TeV range. This typically gives rise to an "effective SUSY" type of sparticle mass spectrum. In this case, we expect gluino pair production as the dominant sparticle creation reaction at LHC. For m_0< 5 TeV, the superpotential parameter \mu > 2 TeV and m_A> 0.8 TeV, greatly restricting neutralino annihilation mechanisms. These latter conclusions are softened if m_0\sim 10-20 TeV or if one proceeds to the NUHM2 model. The standard neutralino abundance tends to be far above WMAP-measured values unless the neutralino is higgsino-like. We remark upon possible non-standard (but perhaps more attractive) cosmological scenarios which can bring the predicted dark matter abundance into accord with the measured value, and discuss the implications for direct and indirect detection of neutralino cold dark matter.

The minimal supergravity (mSUGRA or CMSSM) model is an oft-used framework for exhibiting the properties of neutralino (WIMP) cold dark matter (CDM). However, the recent evidence from Atlas and CMS on a light Higgs scalar with mass m_h\simeq 125 GeV highly constrains the superparticle mass spectrum, which in turn constrains the neutralino annihilation mechanisms in the early universe. We find that stau and stop co-annihilation mechanisms -- already highly stressed by the latest Atlas/CMS results on SUSY searches -- are nearly eliminated if indeed the light Higgs scalar has mass m_h\simeq 125 GeV. Furthermore, neutralino annihilation via the A-resonance is essentially ruled out in mSUGRA so that it is exceedingly difficult to generate thermally-produced neutralino-only dark matter at the measured abundance. The remaining possibility lies in the focus-point region which now moves out to m_0\sim 10-20 TeV range due to the required large trilinear soft SUSY breaking term A_0. The remaining HB/FP region is more fine-tuned than before owing to the typically large top squark masses. We present updated direct and indirect detection rates for neutralino dark matter, and show that ton scale noble liquid detectors will either discover mixed higgsino CDM or essentially rule out thermally-produced neutralino-only CDM in the mSUGRA model.

While the SUSY flavor, CP and gravitino problems seem to favor a very heavy spectrum of matter scalars, fine-tuning in the electroweak sector prefers low values of superpotential mass \mu. In the limit of low \mu, the two lightest neutralinos and light chargino are higgsino-like. The light charginos and neutralinos may have large production cross sections at LHC, but since they are nearly mass degenerate, there is only small energy release in three-body sparticle decays. Possible dilepton and trilepton signatures are difficult to observe after mild cuts due to the very soft p_T spectrum of the final state isolated leptons. Thus, the higgsino-world scenario can easily elude standard SUSY searches at the LHC. It should motivate experimental searches to focus on dimuon and trimuon production at the very lowest p_T(\mu) values possible. If the neutralino relic abundance is enhanced via non-standard cosmological dark matter production, then there exist excellent prospects for direct or indirect detection of higgsino-like WIMPs. While the higgsino-world scenario may easily hide from LHC SUSY searches, a linear e^+e^- collider or a muon collider operating in the \sqrt{s}\sim 0.5-1 TeV range would be able to easily access the chargino and neutralino pair production reactions.

One of the main channels which allows for a large rate of neutralino dark matter annihilation in the early Universe is via the pseudoscalar Higgs A-resonance. In this case, the measured dark matter abundance can be obtained in the minimal supergravity (mSUGRA) model when tan(beta)\sim 50 and 2m_{\tz_1}\sim m_A. We investigate the reaction pp\to b\phi\to b\mu^+\mu^- +X (where \phi =A or H) at the CERN LHC where requiring the tag of a single b-jet allows for amplification of the signal-to-background ratio. The rare but observable Higgs decay to muon pairs allows for a precise measurement of the Higgs boson mass and decay width. We evaluate signal and background using CalcHEP, with muon energy smearing according to the CMS detector. We find that the Higgs width (\Gamma_A) can typically be determined with the accuracy up to \sim 8% (\sim 17%) for m_A\sim 400 (600) GeV assuming 10^3 fb^{-1} of integrated luminosity. Therefore, the pp\to b\phi\to b\mu^+\mu^- +X process provides a unique possibility for direct \Gamma_A measurement at the LHC. While the Higgs width is correlated with the parameter \tan\beta for a given value of m_A, extracting \tan\beta is complicated by an overlap of the A and H peaks, radiative corrections to the b and \tau Yukawa couplings, and the possibility that SUSY decay modes of the Higgs may be open. In the case where a dilepton mass edge from \tz_2\to\ell^+\ell^-\tz_1 is visible, it should be possible to test the relation that 2m_{\tz_1}\sim m_A.

We examine some conditions which are needed to allow the Peccei-Quinn scale (f_a) as high as the SUSY GUT scale (M_{GUT} \sim 10^{16} GeV) in the context of the PQMSSM with either an axino or a neutralino LSP. The main problem in non-SUSY models with f_a\sim M_{GUT} is the generation of an overabundance of axion dark matter (\sim f_a^{7/6}) due to vacuum misalignment. We show that once all the components of the axion supermultiplet are included, the upper limit on f_a can be evaded due to large entropy injection from saxion decays. This large entropy injection also dilutes all other quasi-stable relic densities, naturally evading the BBN constraints and solving the gravitino problem. We find that f_a\sim M_{GUT} can be allowed by relic density/BBN constraints provided that the saxion mass m_s> 50 TeV, the initial saxion field value is of order of the PQ scale and the initial axion mis-alignment angle \theta_i\alt 0.05 . These restrictions can be considerably loosened for f_a\sim 10^{15} GeV. The allowed range for the re-heat temperature (T_R) is strongly dependent on the nature of the LSP. For the axino LSP, T_R \gtrsim 10^8 GeV, while for neutralino LSP any value is allowed. In the latter case, f_a\sim M_{GUT} can be more easily accommodated than in the axino LSP scenario. For f_a\sim M_{GUT} in SUSY models, the dark matter abundance should be dominated by axions, albeit with mass m_a\sim 10^{-10} eV, far below the region currently probed by experiment.

We consider supersymmetric (SUSY) models wherein the strong CP problem is solved by the Peccei-Quinn (PQ) mechanism with a concommitant axion/axino supermultiplet. We examine R-parity conserving models where the neutralino is the lightest SUSY particle, so that a mixture of neutralinos and axions serve as cold dark matter. The mixed axion/neutralino CDM scenario can match the measured dark matter abundance for SUSY models which typically give too low a value of the usual thermal neutralino abundance, such as models with wino-like or higgsino-like dark matter. The usual thermal neutralino abundance can be greatly enhanced by the decay of thermally-produced axinos to neutralinos, followed by neutralino re-annihilation at temperatures much lower than freeze-out. In this case, the relic density is usually neutralino dominated, and goes as \sim (f_a/N)/m_{axino}^{3/2}. If axino decay occurs before neutralino freeze-out, then instead the neutralino abundance can be augmented by relic axions to match the measured abundance. Entropy production from late-time axino decays can diminish the axion abundance, but ultimately not the neutralino abundance. In mixed axion/neutralino CDM models, it may be possible to detect both a WIMP and an axion as dark matter relics. We also discuss possible modifications of our results due to production and decay of saxions. In the appendices, we present expressions for the Hubble expansion rate and the axion and neutralino relic densities in radiation, matter and decaying-particle dominated universes.

A successful implementation of thermal leptogenesis requires the re-heat temperature after inflation T_R to exceed ~2\times 10^9 GeV. Such a high T_R value typically leads to an overproduction of gravitinos in the early universe, which will cause conflicts, mainly with BBN constraints. Asaka and Yanagida (AY) have proposed that these two issues can be reconciled in the context of the Peccei-Quinn augmented MSSM (PQMSSM) if one adopts a mass hierarchy m(sparticle)>m(gravitino)>m(axino), with m(axino) keV. We calculate the relic abundance of mixed axion/axino dark matter in the AY scenario, and investigate under what conditions a value of T_R sufficient for thermal leptogenesis can be generated. A high value of PQ breaking scale f_a is needed to suppress overproduction of axinos, while a small vacuum misalignment angle \theta_i is needed to suppress overproduction of axions. The large value of f_a results in late decaying neutralinos. To avoid BBN constraints, the AY scenario requires a low thermal abundance of neutralinos and high values of neutralino mass. We include entropy production from late decaying saxions, and find the saxion needs to be typically at least several times heavier than the gravitino. A viable AY scenario suggests that LHC should discover a spectrum of SUSY particles consistent with weak scale supergravity; that the apparent neutralino abundance is low; that a possible axion detection signal (probably with m_axion in the sub-micro-eV range) should occur, but no direct or indirect signals for WIMP dark matter should be observed.

Successful implementation of thermal leptogenesis requires re-heat temperatures T_R\agt 2\times 10^9 GeV, in apparent conflict with SUSY models with TeV-scale gravitinos, which require much lower T_R in order to avoid Big Bang Nucleosynthesis (BBN) constraints. We show that mixed axion/axino dark matter can reconcile thermal leptogenesis with the gravitino problem in models with m_{\tG}\agt 30 TeV, a rather high Peccei-Quinn breaking scale and an initial mis-alignment angle \theta_i < 1. We calculate axion and axino dark matter production from four sources, and impose BBN constraints on long-lived gravitinos and neutralinos. Moreover, we discuss several SUSY models which naturally have gravitino masses of the order of tens of TeV. We find a reconciliation difficult in Yukawa-unified SUSY and in AMSB with a wino-like lightest neutralino. However, T_R\sim 10^{10}-10^{12} GeV can easily be achieved in effective SUSY and in models based on mixed moduli-anomaly mediation. Consequences of this scenario include: 1. an LHC SUSY discovery should be consistent with SUSY models with a large gravitino mass, 2. an apparent neutralino relic abundance \Omega_{\tz_1}h^2\alt 1, 3. no WIMP direct or indirect detection signals should be found, and 4. the axion mass should be less than \sim 10^{-6} eV, somewhat below the conventional range which is explored by microwave cavity axion detection experiments.

We calculate the relic abundance of thermally produced neutralino cold dark matter in the general 19 parameter supergravity (SUGRA-19) model. A scan over GUT scale parameters reveals that models with a bino-like neutralino typically give rise to a dark matter density \Omega_{\tz_1}h^2\sim 1-1000, i.e. between 1 and 4 orders of magnitude higher than the measured value. Models with higgsino or wino cold dark matter can yield the correct relic density, but mainly for neutralino masses around 700-1300 GeV. Models with mixed bino-wino or bino-higgsino CDM, or models with dominant co-annihilation or A-resonance annihilation can yield the correct abundance, but such cases are extremely hard to generate using a general scan over GUT scale parameters; this is indicative of high fine-tuning of the relic abundance in these cases. Requiring that m_{\tz_1}\alt 500 GeV (as a rough naturalness requirement) gives rise to a minimal probably dip in parameter space at the measured CDM abundance. For comparison, we also scan over mSUGRA space with four free parameters. Finally, we investigate the Peccei-Quinn augmented MSSM with mixed axion/axino cold dark matter. In this case, the relic abundance agrees more naturally with the measured value. In light of our cumulative results, we conclude that future axion searches should probe much more broadly in axion mass, and deeper into the axion coupling.

Supersymmetric models based on anomaly-mediated SUSY breaking (AMSB) generally give rise to a neutral wino as a WIMP cold dark matter (CDM) candidate, whose thermal abundance is well below measured values. Here, we investigate four scenarios to reconcile AMSB dark matter with the measured abundance: 1. non-thermal wino production due to decays of scalar fields ({\it e.g} moduli), 2. non-thermal wino production due to decays of gravitinos, 3. non-thermal wino production due to heavy axino decays, and 4. the case of an axino LSP, where the bulk of CDM is made up of axions and thermally produced axinos. In cases 1 and 2, we expect wino CDM to constitute the entire measured DM abundance, and we investigate wino-like WIMP direct and indirect detection rates. Wino direct detection rates can be large, and more importantly, are bounded from below, so that ton-scale noble liquid detectors should access all of parameter space for m_{\tz_1}\alt 500 GeV. Indirect wino detection rates via neutrino telescopes and space-based cosmic ray detectors can also be large. In case 3, the DM would consist of an axion plus wino admixture, whose exact proportions are very model dependent. In this case, it is possible that both an axion and a wino-like WIMP could be detected experimentally. In case 4., we calculate the re-heat temperature of the universe after inflation. In this case, no direct or indirect WIMP signals should be seen, although direct detection of relic axions may be possible. For each DM scenario, we show results for the minimal AMSB model, as well as for the hypercharged and gaugino AMSB models.

Over almost all of minimal supergravity (mSUGRA or CMSSM) model parameter space, there is a large overabundance of neutralino cold dark matter (CDM). We find that the allowed regions of mSUGRA parameter space which match the measured abundance of CDM in the universe are highly fine-tuned. If instead we invoke the Peccei-Quinn-Weinberg-Wilczek solution to the strong CP problem, then the SUSY CDM may consist of an axion/axino admixture with an axino mass of order the MeV scale, and where mixed axion/axino or mainly axion CDM seems preferred. In this case, fine-tuning of the relic density is typically much lower, showing that axion/axino CDM (a\tilde{a}CDM) is to be preferred in the paradigm model for SUSY phenomenology. For mSUGRA with a\tilde{a}CDM, quite different regions of parameter space are now DM-favored as compared to the case of neutralino DM. Thus, rather different SUSY signatures are expected at the LHC in the case of mSUGRA with a\tilde{a}CDM, as compared to mSUGRA with neutralino CDM.

Supersymmetric grand unified theories based on the gauge group SO(10) are highly motivated. In the simplest models, one expects t-b-\tau Yukawa coupling unification, in addition to gauge, matter and Higgs unification. Yukawa unification only occurs with very special GUT scale boundary conditions, leading to a spectra with ~10 TeV first and second generation scalars, TeV-scale third generation scalars, and light gauginos. The relic density of neutralino cold dark matter is calculated to be 10^2-10^4 times higher than observation. If we extend the theory with the PQWW solution to the strong CP problem, then instead a mixture of axions and axinos comprises the dark matter, with the measured abundance. Such a solution solves several cosmological problems. We predict a rather light gluino with m(gluino)~300-500 GeV that should be visible in either Tevatron or forthcoming LHC run 1 data. We would also expect ultimately a positive result from relic axion search experiments.

We examine the minimal supergravity (mSUGRA) model under the assumption that the strong CP problem is solved by the Peccei-Quinn mechanism. In this case, the relic dark matter (DM) abundance consists of three components: {\it i}). cold axions, {\it ii}). warm axinos from neutralino decay, and {\it iii}). cold or warm thermally produced axinos. To sustain a high enough re-heat temperature (T_R\agt 10^6 GeV) for many baryogenesis mechanisms to function, we find that the bulk of DM should consist of cold axions, while the admixture of cold and warm axinos should be rather slight, with a very light axino of mass \sim 100 keV. For mSUGRA with mainly axion cold DM (CDM), the most DM-preferred parameter space regions are precisely those which are least preferred in the case of neutralino DM. Thus, rather different SUSY signatures are expected at the LHC in the case of mSUGRA with mainly axion CDM, as compared to mSUGRA with neutralino CDM.

Supersymmetric models with t-b-\tau Yukawa unification at M_{GUT} qualitatively predict a sparticle mass spectrum including first and second generation scalars at the 3--15 TeV scale, third generation scalars at the (few) TeV scale and gluinos in the sub-TeV range. The neutralino relic density in these models typically turns out to lie far above the measured dark matter abundance, prompting the suggestion that instead dark matter is composed of an axion/axino mixture. We explore the axion and thermal and non-thermal axino dark matter abundance in Yukawa-unified SUSY models. We find in this scenario that {\it i}). rather large values of Peccei-Quinn symmetry breaking scale f_a\sim 10^{12} GeV are favored and {\it ii}). rather large values of GUT scale scalar masses \sim 10-15 TeV allow for the re-heat temperature T_R of the universe to be T_R\agt 10^6 GeV. This allows in turn a solution to the gravitino/Big Bang Nucleosynthesis problem while also allowing for baryogenesis via non-thermal leptogenesis. The large scalar masses for Yukawa-unified models are also favored by data on b\to s\gamma and B_s\to \mu^+\mu^- decay. Testable consequences from this scenario include a variety of robust LHC signatures, a possible axion detection at axion search experiments, but null results from direct and indirect WIMP search experiments.

Gravity-mediated SUSY breaking models with R-parity conservation give rise to dark matter in the universe. I review neutralino dark matter in the minimal supergravity model (mSUGRA), models with non-universal soft SUSY breaking terms (NUSUGRA) which yield a well-tempered neutralino, and models with unified Yukawa couplings at the GUT scale (as may occur in an SO(10) SUSY GUT theory). These latter models have difficulty accommodating neutralino dark matter, but work very well if the dark matter particles are axions and axinos.

The NSF has chosen the site for the Deep Underground Science and Engineering Laboratory (DUSEL) to be in Lead, South Dakota. In fact, the state of South Dakota has already stepped up to the plate and contributed its own funding for the proposed lab, see http://www.sanfordlaboratoryathomestake.org/index.html. The final decision by NSF for funding the Initial Suite of Experiments for DUSEL will be made early in 2009. At that time the NSF Science Board must make a decision. Of order 200 experimentalists have already expressed an interest in performing experiments at DUSEL. In order to assess the interest of the theoretical community, the Center for Cosmology and Astro-Particle Physics (CCAPP) at The Ohio State University (OSU) organized a 3-day DUSEL Theory Workshop in Columbus, Ohio from April 4 - 6, 2008. The workshop focused on the scientific case for six proposed experiments for DUSEL: long baseline neutrino oscillations, proton decay, dark matter, astrophysical neutrinos, neutrinoless double beta decay and N-Nbar oscillations. The outcome of this workshop is the DUSEL Theory White paper addressing the scientific case at a level which may be useful in the decision making process for policy makers at the NSF and in the U.S. Congress. In order to assess the physics interest in the DUSEL project we have posted the DUSEL Theory White paper on the following CCAPP link http://ccapp.osu.edu/whitepaper.html . Please read the white paper and, if you are interested, use the link to show your support by co-signing the white paper.

Simple SUSY GUT models based on the gauge group SO(10) require t-b-\tau Yukawa coupling unification, in addition to gauge coupling and matter unification. The Yukawa coupling unification places strong constraints on the expected superparticle mass spectrum, with scalar masses \sim 10 TeV while gaugino masses are quite light. A problem generic to all supergravity models comes from overproduction of gravitinos in the early universe: if gravitinos are unstable, then their late decays may destroy the predictions of Big Bang nucleosynthesis. We present a Yukawa-unified SO(10) SUSY GUT scenario which avoids the gravitino problem, gives rise to the correct matter-antimatter asymmetry via non-thermal leptogenesis, and is consistent with the WMAP-measured abundance of cold dark matter due to the presence of an axino LSP. To maintain a consistent cosmology for Yukawa-unified SUSY models, we require a re-heat temperature T_R\sim 10^6-10^7 GeV, an axino mass around \sim 0.1-10 MeV, and a PQ breaking scale f_a\sim 10^{12} GeV.

An abundance of astrophysical evidence indicates that the bulk of matter in the universe is made up of massive, electrically neutral particles that form the dark matter (DM). While the density of DM has been precisely measured, the identity of the DM particle (or particles) is a complete mystery. In fact, within the laws of physics as we know them (the Standard Model, or SM), none of the particles have the right properties to make up DM. Remarkably, many new physics extensions of the SM -- designed to address theoretical issues with the electroweak symmetry breaking sector -- require the introduction of new particles, some of which are excellent DM candidates. As the LHC era begins, there are high hopes that DM particles, along with their associated new matter states, will be produced in pp collisions. We discuss how LHC experiments, along with other DM searches, may serve to determine the identity of DM particles and elucidate the associated physics. Most of our discussion centers around theories with weak-scale supersymmetry, and allows for several different DM candidate particles.

In the post-LEP2 era, and in light of recent measurements of the cosmic abundance of cold dark matter (CDM) in the universe from WMAP, many supersymmetric models tend to predict 1. an overabundance of CDM and 2. pessimistically low rates for direct detection of neutralino dark matter. However, in models with a ``well-tempered neutralino'', where the neutralino composition is adjusted to give the measured abundance of CDM, the neutralino is typically of the mixed bino-wino or mixed bino-higgsino state. Along with the necessary enhancement to neutralino annihilation rates, these models tend to give elevated direct detection scattering rates compared to predictions from SUSY models with universal soft breaking terms. We present neutralino direct detection cross sections from a variety of models containing a well-tempered neutralino, and find cross section asymptotes with detectable scattering rates. These asymptotic rates provide targets that various direct CDM detection experiments should aim for. In contrast, in models where the neutralino mass rather than its composition is varied to give the WMAP relic density via either resonance annihilation or co-annihilation, the neutralino remains essentially bino-like, and direct detection rates may be below the projected reaches of all proposed experiments.

Low energy antideuterons suffer a very low secondary and tertiary astrophysical background, while they can be abundantly synthesized in dark matter pair annihilations, therefore providing a privileged indirect dark matter detection technique. The recent publication of the first upper limit on the low energy antideuteron flux by the BESS collaboration, a new evaluation of the standard astrophysical background, and remarkable progresses in the development of a dedicated experiment, GAPS, motivate a new and accurate analysis of the antideuteron flux expected in particle dark matter models. To this extent, we consider here supersymmetric, universal extra-dimensions (UED) Kaluza-Klein and warped extra-dimensional dark matter models, and assess both the prospects for antideuteron detection as well as the various related sources of uncertainties. The GAPS experiment, even in a preliminary balloon-borne setup, will explore many supersymmetric configurations, and, eventually, in its final space-borne configuration, will be sensitive to primary antideuterons over the whole cosmologically allowed UED parameter space, providing a search technique which is highly complementary with other direct and indirect dark matter detection experiments.

Within the minimal supergravity model (mSUGRA) framework, the expectation for the relic density of neutralinos exceeds the WMAP determination unless neutralinos a) have a significant higgsino component, b) have a mass close to half that of a heavy Higgs boson, or c)can efficiently co-annihilate with a charged or colored particle. Within a 1-parameter extension of the mSUGRA model which includes non-universal Higgs masses, we show that agreement with the WMAP data can be obtained over a wide range of mSUGRA parameters for scenarios a) and b), so that the phenomenological implications may be much more diverse than in mSUGRA. We show that direct and/or indirect detection of neutralino dark matter should be possible at various current and planned facilities.

We examine the prospects for supersymmetry discovery in the minimal supergravity (mSUGRA) model via indirect detection of neutralino dark matter. We investigate rates for muon detection in neutrino telescopes, and detection of photons, positrons and anti-protons by balloon and space based experiments. We compare the discovery reach in these channels with the reach for direct detection of dark matter, and also with the reach of collider experiments such as Fermilab Tevatron, CERN LHC and a $\sqrt{s}=0.5-1$ TeV linear $e^+e^-$ collider. We pay particular attention to regions of model parameter space in accord with recent WMAP results on the dark matter density of the universe. We find that the IceCube neutrino telescope as well as 3rd generation direct dark matter detection experiments should be able to cover the {\it entire} WMAP allowed portion of the hyperbolic branch/focus point (HB/FP) region of parameter space. This is in contrast to the case of the CERN LHC or a linear $e^+e^-$ collider, where only a fraction of the HB/FP region can be accessed. In addition, we show that detection of $\gamma$s, $e^+$s and $\bar{p}$s should occur in much of the HB/FP region, as well as in the low $m_{1/2}$ portion of the $A$ annihilation funnel, and will be complementary to searches via colliders in these regions.

The lightest neutralino of R-parity conserving supersymmetric models serves as a compelling candidate to account for the presence of cold dark matter in the universe. In the minimal supergravity (mSUGRA) model, a relic density can be found in accord with recent WMAP data for large values of the parameter $\tan\beta$, where neutralino annihilation in the early universe occurs via the broad s-channel resonance of the pseudoscalar Higgs boson $A$. We map out rates for indirect detection of neutralinos via 1. detection of neutrinos arising from neutralino annihilation in the core of the earth or sun and 2. detection of gamma rays, antiprotons and positrons arising from neutralino annihilation in the galactic halo. If indeed $A$-resonance annihilation is the main sink for neutralinos in the early universe, then signals may occur in the gamma ray, antiproton and positron channels, while a signal in the neutrino channel would likely be absent. This is in contrast to the hyperbolic branch/focus point (HB/FP) region where {\it all} indirect detection signals are likely to occur, and also in contrast to the stau co-annihilation region, where {\it none} of the indirect signals are likely to occur.

We evaluate neutralino-nucleon scattering rates in several well-motivated supersymmetric models, and compare against constraints on the neutralino relic density, BF( b\to s\gamma ) as well as the muon anomalous magnetic moment a_\mu . In the mSUGRA model, the indirect constraints favor the hyperbolic branch/focus point (HB/FP) region of parameter space, and in fact this region is just where neutralino-nucleon scattering rates are high enough to be detected in direct dark matter search experiments! In Yukawa unified SUSY SO(10) models with scalar mass non-universality, the relic density of neutralinos is almost always above experimental bounds, while the corresponding direct detection rates are below experimental levels. Conversely, in five dimensional SO(10) models where gauge symmetry breaking is the result of compactification of the extra dimension, and supersymmetry breaking is communicated via gaugino mediation, the relic density is quite low, while direct detection rates can be substantial.

Recent results from the WMAP measurements of the cosmic background radiation yield very tight constraints on the relic density of supersymmetric cold dark matter. We combine the WMAP constraint with those from the anomalous magnetic moment of the muon and the b -> s gamma branching fraction in a chi^2 determination over the minimal supergravity model (mSUGRA) parameter space. The most favored region of mSUGRA parameter space for almost all tan(beta) values is the hyperbolic branch/focus point (HB/FP) region, with moderate to small values of superpotential Higgs mass |mu| and large GUT scale scalar mass m_0. These favored regions of mSUGRA parameter space can be probed by direct search experiments for supersymmetric dark matter. An exception to the HB/FP region can occur at very large tan(beta) with positive mu values, where wide regions allow resonance annihilation of neutralinos in the early universe.

We evaluate the relic density of neutralinos in the minimal supergravity (mSUGRA) model. All 2-->2 neutralino annihilation diagrams, as well as all processes involving sleptons, charginos, neutralinos and third generation squarks are included. Relativistic thermal averaging of the velocity times cross sections is performed. We find that co-annihilation effects are only important on the edges of the model parameter space, where some amount of fine-tuning is necessary to obtain a reasonable relic density. Alternatively, at high tan(beta), annihilation through very broad Higgs resonances gives rise to an acceptable neutralino relic density over broad regions of parameter space where little or no fine-tuning is needed. Finally, we compare our results against the reach of various e+e- and hadron colliders for supersymmetric matter.

Working within the framework of the minimal supergravity model with gauge coupling unification and radiative electroweak symmetry breaking, we evaluate the cosmological relic density from lightest neutralinos produced in the early universe. Our numerical calculation is distinct in that it involves direct evaluation of neutralino annihilation cross sections using helicity amplitude techniques, and thus avoids the usual expansion as a power series in terms of neutralino velocity. Thus, our calculation includes relativistic Boltzmann averaging, neutralino annihilation threshold effects, and proper treatment of integration over Breit-Wigner poles. We map out regions of parameter space that give rise to interesting cosmological dark matter relic densities. We compare these regions with recent calculations of the reach for supersymmetry by LEP2 and Tevatron Main Injector era experiments. The cosmologically favored regions overlap considerably with the regions where large trilepton signals are expected to occur at the Tevatron. The CERN LHC $pp$ collider can make a thorough exploration of the cosmologically favored region via gluino and squark searches. In addition, over most of the favored region, sleptons ought to be light enough to be detectable at both LHC and at a $\sqrt{s}=500$ GeV $e^+e^-$ collider.