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.

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.

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.

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.

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.