Is the magnitude of the Peccei-Quinn scale set by the landscape?

TL;DR

The paper investigates whether the string landscape influences the Peccei-Quinn scale, suggesting it is set around 10^{11}-10^{12} GeV due to SUSY breaking considerations, with implications for dark matter abundance.

Contribution

It proposes a landscape-based explanation for the PQ scale within a SUSY axion model, linking it to soft SUSY breaking terms and dark matter constraints.

Findings

PQ scale likely around 10^{11}-10^{12} GeV due to SUSY considerations

Anthropic setting of PQ scale is unlikely because of dark matter overproduction

Preferred f_a value aligns with intermediate scale dark matter abundance

Abstract

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.

Figures

Tables

Table 1. Table 1 : Input parameters (TeV) and masses (GeV) for a landscape SUSY benchmark point from the NUHM3 model with m t = 173.2 subscript 𝑚 𝑡 173.2 m_{t}=173.2 GeV using Isajet 7.88 [ 51 ] .

parameter	landSUSY
$m_0(1,2)$	$16$ TeV
$m_0\left(3\right)$	$5$ TeV
$m_{1/2}$	1.5 TeV
$A_0$	-7 TeV
$\tan \beta$	10
$\mu$	0.2 TeV
$m_A$	3 TeV
$m_{\tilde{g}}$	3619 GeV
$m_{{\tilde{u}}_L}$	16211 GeV
$m_{{\tilde{u}}_R}$	16264 GeV
$m_{{\tilde{e}}_R}$	15956 GeV
$m_{{\tilde{t}}_1}$	1294 GeV
$m_{{\tilde{t}}_2}$	3561 GeV
$m_{{\tilde{b}}_1}$	3605 GeV
$m_{{\tilde{b}}_2}$	4999 GeV
$m_{{\tilde{\tau }}_1}$	4749 GeV
$m_{{\tilde{\tau }}_2}$	4982 GeV
$m_{{\tilde{\nu }}_{\tau }}$	4951 GeV
$m_{{\tilde{\chi }}_1^{\pm }}$	210 GeV
$m_{{\tilde{\chi }}_2^{\pm }}$	1312 GeV
$m_{{\tilde{\chi }}_1^0}$	200 GeV
$m_{{\tilde{\chi }}_2^0}$	207 GeV
$m_{{\tilde{\chi }}_3^0}$	688 GeV
$m_{{\tilde{\chi }}_4^0}$	1320 GeV
$m_h$	125 GeV
${\Omega }_{{\tilde{\chi }}_1^0}^{std}{h}^2$	0.01
$BF\left(b\to s\gamma \right)\times {10}^4$	$3.0$
$BF\left(B_s\to {\mu }^{+}{\mu }^{-}\right)\times {10}^9$	$3.8$
${\sigma }^{SI}({\tilde{\chi }}_1^0,p)$ (pb)	$1.0\times {10}^{-9}$
${\sigma }^{SD}\left({\tilde{\chi }}_1^{0}p\right)$ (pb)	$2.0\times {10}^{-5}$
${⟨\sigma v⟩|}_{v\to 0}$ (cm3/sec)	$2\times {10}^{-25}$
${\Delta }_{\mathrm{EW}}$	23.3

Table 2. Table 2 : ℤ 24 R superscript subscript ℤ 24 𝑅 \mathbb{Z}_{24}^{R} and PQ charge assignments for various superfields of the hyCCK model.

multiplet	$H_u$	$H_d$	$Q_i$	$L_i$	$U_i^c$	$D_i^c$	$E_i^c$	$N_i^c$	X	Y
${\mathbb{Z}}_{24}^R$ charge	16	12	5	9	5	9	5	1	-1	5
PQ charge	-1	-1	1	1	0	0	0	0	1	-3

Equations

$$W_{hyCCK}$$

$$V_{hyCCK}\ni \frac{f^2}{m_P^2}\left[9|{\varphi }_X{|}^4|{\varphi }_Y{|}^2+|{\varphi }_X{|}^6\right]+m_X^2|{\varphi }_X{|}^2+m_Y^2|{\varphi }_Y{|}^2+(f{A}_f{\varphi }_X^3{\varphi }_Y/m_P+h.c.).$$

$${\Omega }_a^{\mathrm{CO}}{h}^2\simeq 0.23f({\theta }_i){\theta }_i^2{\left(\frac{f_a/N_{\mathrm{DW}}}{10^{12}\mathrm{GeV}}\right)}^{7/6}$$

$$\rho$$

$$\mathrm{with}\Omega$$

Full text

OU-HEP-190501

**Is the magnitude of the Peccei-Quinn scale

set by the landscape?

** Howard Baer1111Email: baer@ou.edu , Vernon Barger2222Email: barger@pheno.wisc.edu, Dibyashree Sengupta1333Email: Dibyashree.Sengupta-1@ou.edu

Hasan Serce2444Email: serce@ou.edu, Kuver Sinha1555Email: kuver.sinha@ou.edu and Robert Wiley Deal1666Email: rwileydeal@ou.edu

1*Dept. of Physics and Astronomy, University of Oklahoma, Norman, OK 73019, USA

2*Dept. of Physics, University of Wisconsin, Madison, WI 53706 USA

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}\sim 100$ GeV. Such a picture leads to the prediction that $m_{h}\simeq 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 $\mathbb{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}\sim 10-100$ TeV places the mixed axion-neutralino dark matter abundance into the intermediate scale sweet zone where $f_{a}\sim 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}\sim 10^{12}$ GeV.

1 Introduction

The Standard Model (SM) is beset with several fine-tuning problems that call for new physics beyond the SM. These include:

The gauge hierarchy problem[1] wherein the weak scale $m_{weak}\simeq m_{W,Z,h}\simeq 100$ GeV $\sim 10^{-14}m_{\rm GUT}$ (where $m_{\rm GUT}\simeq 2\times 10^{16}$ GeV). 2. 2.

The cosmological constant (CC) problem[2] wherein the measured value of the cosmological constant $\Lambda\simeq 10^{-120}m_{P}^{4}$ (where naively it is expected that $\Lambda\simeq m_{P}^{4}$ with $m_{P}$ being the reduced Planck mass). 3. 3.

The strong CP problem[3] wherein the $\bar{\theta}$ coefficient of the CP-violating $\bar{\theta}(g_{s}^{2}/32\pi)G_{\mu\nu A}\tilde{G}_{A}^{\mu\nu}$ QCD Lagrangian term is $\lesssim 10^{-10}$ times its expected value from the ’t Hooft theta vacuum solution to the $U(1)_{A}$ problem[4].

The most elegant solution to problem #1 is to extend the Poincare’ group of spacetime symmetries to include weak scale broken supersymmetry (WSS)[5, 6] wherein all quadratic divergences to $m_{h}^{2}$ necessarily cancel. WSS is in fact supported by four sets of measurements:

• •

The measured values of the gauge couplings actually unify under Minimal Supersymmetric Standard Model (MSSM) evolution whereas they do not under SM evolution[7].

• •

The measured value of the top quark mass is just right to generate a radiative breakdown of EW symmetry in the MSSM (assumed valid up to scales $Q\simeq m_{\rm GUT}$)[8].

• •

The measured value of $m_{h}\simeq 125$ GeV falls squarely within the narrow window of MSSM required values where $m_{h}\lesssim 135$ GeV[9] and

• •

the measured value of $m_{W}\simeq 80.4$ GeV favors heavy WSS over the SM assuming the measured value of $m_{t}\simeq 173.2$ GeV[10].

Early concern over the non-appearance of SUSY particles at LHC Run 2 has been allayed by renewed scrutiny of naturalness measures[11, 12]: updated analyses now require only $m_{\tilde{g}}\lesssim 6$ TeV and $m_{\tilde{t}_{1}}\lesssim 3$ TeV[12, 13, 14] as compared to present LHC limits that $m_{\tilde{g}}\gtrsim 2$ TeV and $m_{\tilde{t}_{1}}\gtrsim 1$ TeV. In the MSSM, only higgsinos are required by naturalness to have weak scale masses, and these particles are very difficult (but not impossible) to see at LHC[16, 15].

Some understanding of the magnitude of the cosmological constant has emerged from the string theory landscape of vacua solutions[17]. In a discretuum[18] of vacua states with all possible values for $\Lambda:-m_{P}^{4}\rightarrow+m_{P}^{4}$, then one expects $\Lambda$ as large as possible subject to the requirement that galaxies be able to condense, which is a seeming precondition for a pocket universe containing sentient observers. Indeed, Weinberg[19] used such reasoning to predict the value of $\Lambda$ to within a factor of several of its value which was measured more than a decade later. Such an anthropic solution to the CC problem emerges naturally from string theory incuding a vast landscape of flux vacua, estimated as $\sim 10^{500}$ such possibilities[20].

The most elegant solution to the strong CP problem involves the introduction of a new global $U(1)$ PQ symmetry[21] which is spontaneously broken (SSB) at some scale $f_{a}\sim 10^{9}-10^{16}$ GeV.111In accord with the PDG[22], we take $f_{A}\equiv f_{a}/N_{\rm DW}$ where $N_{\rm DW}$ is the domain-wall number which is $N_{\rm DW}=6$ for the DFSZ axion model assumed here. The (pseudo-)Goldstone boson which emerges from SSB, the axion $a$[23], allows for a dynamical relaxation of the $G\tilde{G}$ QCD Lagrangian terms to zero thus solving the strong CP problem. A remnant of the PQ procedure is the existence of a physical axion particle which also turns out to be a solid candidate for cold dark matter (CDM) in the universe[24]. Remnant axion CDM is being searched for at a variety of experiments, the most sensitive of which are the microwave cavity searches[25]. While the PQ axion solution to the strong CP problem is indeed compelling, it is beset by two problems of its own.

• •

Global symmetries are not respected by gravitational interactions and thus the PQ symmetry is not expected to be fundamental[26, 27]. Instead, PQ symmetry is expected to emerge as an accidental, approximate symmetry from some more fundamental gravity-safe symmetry. To be gravity-safe, the resulting PQ symmetry must be of exceptionally high quality: the PQ breaking contributions to the axionic potential must be suppressed by at least eight powers of $m_{P}$[26] in order for $\bar{\theta}\lesssim 10^{-10}$: $V_{PQB}\sim\phi^{10}/m_{P}^{8}$ (where $\phi$ stands for generic scalar fields). This is a very high bar to hurdle!

• •

In string theory, many candidate axions can emerge, but with a PQ scale $f_{a}\sim m_{\rm GUT}$ to $m_{\rm string}$[28, 29]. Meanwhile, cosmological (dark matter) constraints seem to require $f_{a}\sim 10^{11}-10^{12}$ GeV[24]. A further problem then is: what accounts for the apparent suppression of the PQ breaking scale?

The goal of this paper is to examine the second of these axionic concerns in the context of the string theory landscape: can the magnitude of the PQ breaking scale be understood from landscape considerations within a well-motivated model for axion (and WIMP) dark matter? A number of previous works have also addressed this question and these will be briefly reviewed in Subsec. 1.1. The first of the PQ concerns–gravity safety– was recently addressed within the context of anomaly-free (up to a Green-Schwarz term) discrete $R$ symmetries which forbid the SUSY $\mu$ term while respecting grand unification conditions[30]. It was found in Ref. [31] that two closely related SUSY axion models– dubbed hybrid CCK[32] (hyCCK) and hybrid SPM[33] (hySPM)– were found to be gravity-safe under a $\mathbb{Z}_{24}^{R}$ discrete $R$-symmetry which also led to

• •

suppression of the SUSY $\mu$ term,

• •

suppression of the various renormalizable $R$-parity violating terms and

• •

suppression of the dangerous dimension-5 proton decay operators,

all while allowing for see-saw neutrinos. In such a setting, the $\mathbb{Z}_{24}^{R}$ symmetry (and resulting accidental approximate $U(1)_{PQ}$ symmetry) are broken as a consequence of SUSY breaking. Thus, the PQ scale emerges as a derived value depending on the soft SUSY breaking terms. We will assume the hyCCK model in Sec. 2 and derive a probability distribution for the resulting magnitude of the PQ scale $f_{a}$ for an assumed upper bound on the allowed dark matter abundance in pocket universes. Requiring a (pocket) universe without overproduction of mixed axion-WIMP dark matter by a (arbitrary) factor four excess beyond its measured value then leads to a most probable value of $f_{a}\sim 10^{14}$ GeV. We also derive a probability distribution for the initial axion misalignment angle $\theta_{i}$ which tends to favor smaller values of its expected range. Typically this leads to a large overproduction of dark matter beyond its observed value. If one assumes axion-only dark matter, then one may compensate for the large value of $f_{a}$ by allowing for a small value of the initial axion misalignment angle $\theta_{i}\sim 0$[34]. In our approach, where SUSY stabilizes the weak scale and the $\mathbb{Z}_{24}^{R}$ symmetry yields gravity-safety, then WIMPs are also overproduced and a tiny value of $\theta_{i}$ cannot save the day for $f_{a}\gtrsim 10^{14}$ GeV.

Some motivation for our approach comes from earlier analyses by Douglas which explored a statistical approach to the magnitude of the SUSY breaking scale. In Ref’s [35] and [36], it is assumed that all real-valued SUSY breaking scales are equally likely in a fertile patch of the landscape which contains the MSSM as the low energy effective theory. In the case of $F$-term SUSY breaking, since the $F$ terms are complex-valued fields and the magnitude of SUSY breaking depends on the modulus of $\langle F\rangle$, then one expects the magnitude of soft terms $m_{soft}\sim m_{3/2}\sim m_{hidden}^{2}/m_{P}$ to enjoy a linearly increasing statistical distribution in the landscape. For a variety of hidden sector $F$ and $D$-term fields contributing to SUSY breaking, then one expects instead $m_{soft}^{n}$ where $n=2n_{F}+n_{D}-1$ and where $n_{F}$ is the number of contributing $F$-term SUSY breaking fields and $n_{D}$ is the number of $D$-term SUSY breaking fields.

Naively, the increasing soft term prior would suggest soft terms most probably at energy scales far beyond the weak scale. However, Agrawal et al.[37] have computed that if the weak scale is increasing by a factor of 2-5 beyond its measured value, then nuclear physics is modified in ways that are unlikely to lead to a livable universe (as we understand it). Requiring SUSY contributions to the weak scale to be within a factor of a few (four) of its measured value is the same as requiring the naturalness measure[11, 12] $\Delta_{\rm EW}\lesssim 30$[38]. By requiring that the weak scalar potential is properly broken to $SU(3)_{C}\times U(1)_{EM}$ and no contribution to the weak scale exceeds a factor four beyond the weak scale (in the presence of a natural value of $\mu\sim 200$ GeV which emerges from our assumed solution to the SUSY $\mu$ problem), then the distribution of soft terms becomes bounded from above. For the case of a mild $n=1$ statistical draw on soft terms, then there is large mixing in the top-squark mass matrix leading to enhanced Higgs mass radiative corrections and hence $m_{h}\simeq 125$ GeV whilst the gluinos and squarks are pulled beyond LHC search limits[39, 40]. An exception is the SUSY-preserving $\mu$ term which directly contributes to the magnitude of the weak scale and thus must be of order $\mu\sim 100-350$ GeV. The $\mu$ parameter leads to a rather light set of four higgsinos whose parameter space is only now beginning to be explored at LHC via the soft dilepton plus jet plus $\not\!\!\!{E_{T}}$ signature[16]. Thus, in the SUSY landscape picture with a mild statistical draw towards large soft terms, the prediction is that the LHC will see exactly that which it does see: a light Higgs of mass $m_{h}\simeq 125$ GeV with as yet no sign of sparticles[40].

1.1 Review of some previous work

Here, we briefly review some previous related works to give context for our directions.

• •

In Ref. [34], Linde assumes a scenario where inflation continues past the PQ phase transition (which alleviates the axion domain wall problem) giving rise to a uniform axion field strength $a\equiv\theta_{i}f_{a}$ throughout the observable universe. He argues that an increased value of initial axion field strength by a factor ten in pocket universes leads to an increase in matter content of galaxies by a factor $\sim 10^{8}$ likely leading to conditions hostile to life as we know it. Since $\Omega_{a}h^{2}\sim\theta_{i}^{2}f_{a}^{7/6}$, then small values of $\theta_{i}$ are favored and hence values of $f_{a}\gg 10^{12}$ GeV can be accommodated.

• •

In Hellerman and Walcher Ref. [41], the authors examine axion CDM and galaxy formation in a multiverse setting where all other parameters of $\Lambda$CDM models are fixed at the values of our pocket universe. By requiring structure formation before the onset of cosmological constant domination, the authors derive upper bounds on the ratio of dark matter-to-baryons $\rho_{DM}/\rho_{B}$ which are typically of order $10^{4}-10^{5}$ while our universe sits at the lower edge of allowed values $\sim 5$. They conclude that anthropic constraints are unlikely to explain our observed value of $\rho_{DM}/\rho_{B}\sim 5$.

• •

In Wilczek and in Tegmark et al. Ref’s. [42, 43], a model with axion-only CDM is assumed where PQ breaking occurs before the end of inflation. For uniform sampling of $\theta_{i}$, then the prior distribution of axion CDM is the mild distribution $f(\rho_{\rm CDM})\sim 1/\sqrt{\rho_{\rm CDM}}$. Imposing anthropic limits from galaxy formation, star and black hole formation, solar system stability and stiff upper limits on the density of habitable halos, then it is found that the measured DM density is comparable to expectations from the multiverse.

• •

Freivogel in Ref. [44] adopts the same axion-only CDM model as in Ref [43] but then uses Bousso’s causal diamond measure to provide selection constraints. He fixes $\rho_{B}/\rho_{\gamma}$ to its observed value, fixes $f_{a}\gg 10^{12}$ GeV and allows only the misalignment angle $\theta_{i}$ to vary uniformly and calculates the probability of observing a particular value of $\xi=\rho_{\rm DM}/\rho_{B}$ which is allowed to vary. He finds 68% of observers see $\xi\leq 15$ and 95% of observers see $\xi\leq 65$. In our pocket universe with $\xi\sim 5$, he concludes that the observed value of $\xi$ is then reasonable. However, it is not completely improbable to observe $\xi$ values different from our universe.

Some related papers[45, 46, 47, 48] argue that comparable values for $\Omega_{\Lambda}$ and $\Omega_{\rm DM}$ should emerge from the multiverse.

In these works, a particularly parsimonious approach– of the SM amended with a simple axionic extension– is adopted. While admittedly simple, it contains two problems which make them likely unrealistic. First, there is nothing in them to stabilize the weak scale, and thus we would expect a value for the weak scale far beyond its measured value. The inclusion of WSS solves this problem. Second, it is well known that the global $U(1)_{PQ}$ needed for an axionic solution to the strong CP problem is intrinsically incompatible with the inclusion of gravitation and with embedding in string theory in particular. In the following, we adopt a particular SUSY axion model based upon a more fundamental discrete $R$-symmetry which may arise from compactification of extra space dimensions in string theory. For a strong enough discrete symmetry, in this case $\mathbb{Z}_{24}^{R}$, then the emerging accidental, approximate global PQ symmetry is sharp enough such as to allow for the PQ solution to the strong CP problem.

2 Probability distributions for the PQ scale and

$\theta_{i}$ from a SUSY axion model based upon a gravity-safe $\mathbb{Z}_{24}^{R}$ symmetry

In this Section, we first introduce the gravity-safe hyCCK SUSY axion model where the global PQ symmetry emerges as an accidental, approximate global symmetry from a more fundamental $\mathbb{Z}_{24}^{R}$ symmetry. Next, we review calculation of the mixed axion-WIMP dark matter abundance using eight coupled Boltzmann equations. Then, we combine statistical selection of the PQ scale with the requirement against overproduction of dark matter within pocket universes which populate an eternally-inflating[49] multiverse to derive probability distributions for the magnitude of the PQ energy scale and also for the initial axion misalignment angle $\theta_{i}$.

2.1 MSSM augmented by gravity-safe PQ sector based on $\mathbb{Z}_{24}^{R}$ discrete $R$-symmetry

In Ref. [39], probability distributions for Higgs and sparticle masses were derived from the string landscape assuming Douglas’[35, 36] $m_{soft}^{n}$ (with $n=0,1,2$) soft term prior coupled with a veto $\Theta(30-\Delta_{\rm EW})$ on non-standard vacua or too large of SUSY contributions to the weak scale. The $n=1$ sampling generates a probability distribution for $m_{h}$ which peaked at $m_{h}\simeq 125$ GeV while sparticle mass predictions were characterized by

$m_{\tilde{g}}\sim 4\pm 2$ TeV, 2. 2.

$m_{\tilde{t}_{1}}\sim 1.5\pm 0.5$ TeV, 3. 3.

$m_{A}\sim 3\pm 2$ TeV, 4. 4.

$m_{\widetilde{\chi}^{\pm}_{1},\widetilde{\chi}^{0}_{1,2}}\sim 200\pm 100$ GeV and 5. 5.

$m_{\tilde{q},\tilde{\ell}}\sim 20\pm 10$ TeV.

In accord with these results, we will adopt for illustrative purposes a SUSY benchmark point from the NUHM3 model[50] with parameters $m_{0}(1,2)=16$ TeV, $m_{0}(3)=5$ TeV, $m_{1/2}=1.5$ TeV, $A_{0}=-7$ TeV, $\tan\beta=10$ with $\mu=200$ GeV and $m_{A}=3$ TeV. We generate sparticle mass spectra using Isajet 7.88[51] and find the spectra provided in Table 1 and labelled as landSUSY (landscape SUSY). The spectra lie well within the landscape SUSY predictions for an $n=1$ mild draw to large soft terms[39]. Our final results will hardly depend on significant deviations from the landSUSY benchmark model which mainly sets a natural value for the $\mu$ parameter and the saxion and axino branching fractions which enter into the relic density calculation.

We will augment the landscape SUSY spectra with a PQ sector from the hybrid CCK[32] model where the superpotential is given by

[TABLE]

where we have introduced MSSM singlet fields $X$ and $Y$ with PQ charges listed in Table 2. The hyCCK model is thus a particular example of a SUSY DFSZ axion model where PQ field $X$ couples to the two Higgs doublets thus providing a solution to the SUSY $\mu$ problem222When the PQ field $Y$ couples to the two Higgs doublets, the resulting model is hySPM model which is also an example of a gravity-safe model that solves strong CP problem and the SUSY $\mu$ problem. This model gives $f_{a}$ values similar to those obtained in the hyCCK model[31].. The model assumes an underlying $\mathbb{Z}_{24}^{R}$ discrete $R$ symmetry which forbids

1. renormalizable RPV terms, 2. the usual $\mu$ term and
2. dimension-5 proton decay operators while allowing the required Yukawa couplings and see-saw neutrino terms[30].

The lowest order PQ breaking superpotential terms are $X^{8}Y^{2}/m_{P}^{7}$, $Y^{10}/m_{P}^{7}$ and $X^{4}Y^{6}/m_{P}^{7}$ leading to lowest order scalar potential terms suppressed by powers of $m_{P}^{8}$: thus, the underlying $\mathbb{Z}_{24}^{R}$ symmetry renders the model gravity-safe according to the KM-R requirements[26]. The required $U(1)_{PQ}$ symmetry arises as an accidental approximate global symmetry which arises from an underlying more fundamental discrete $\mathbb{Z}_{24}^{R}$ symmetry.

The scalar potential for hyCCK, augmented with the corresponding soft SUSY breaking terms, is given by

[TABLE]

Minimization conditions for the hyCCK model can be found in Ref. [52]. The scalar potential develops a non-zero minimum at $\langle\phi_{X}\rangle\equiv v_{X}$ and $\langle\phi_{Y}\rangle\equiv v_{Y}$ for a sufficiently large soft term $-A_{f}$, thus breaking the underlying $\mathbb{Z}_{24}^{R}$ and accidental, approximate $PQ$ symmetries. The PQ breaking vev is given by $v_{PQ}=\sqrt{v_{X}^{2}+9v_{Y}^{2}}$ where then $f_{a}=\sqrt{2}v_{PQ}$. In accord with expectations from supergravity models, we will assume $m_{X}=m_{Y}=m_{\tilde{a}}=m_{s}\equiv m_{3/2}$[53]. Thus, in this model, the PQ scale is a derived consequence of SUSY breaking.

The calculated value of $f_{a}$ is given in Fig. 1 as a function of the $-A_{f}$ soft term assuming various values of $m_{X}=m_{Y}=m_{0}(1,2)\equiv m_{3/2}$ and three different values of $f$. From Fig. 1, we see that $f_{a}$ has a monotonically increasing value with increasing $|-A_{f}|$. For a particular value of $-A_{f}$ and $m_{X}=m_{Y}=m_{0}(1,2)\equiv m_{3/2}$ if the value of $f$ is reduced by a factor of 2, then $f_{a}$ increases by approximately $41\%$ and if the value of $f$ is increased by a factor of 2, then $f_{a}$ decreases by approximately $41\%$. Since $-A_{f}$ doesn’t contribute directly to the determination of the weak scale, then there is no (anthropic) upper bound on its value and one might expect $-A_{f}$ and hence $f_{a}$ to lie far beyond the well-known cosmological sweet spot where $f_{a}\sim 10^{11}-10^{12}$ GeV[24].

2.2 Relic density of mixed axion-WIMP dark matter

The evaluation of mixed axion-WIMP dark matter from SUSY axion models is more complicated than simply adding the WIMP thermal abundance to the coherent-oscillation-produced axions. To evaluate the mixed neutralino-axion relic density, we apply the eight-coupled-Boltzmann equation computer code developed in Ref’s [54, 55]. For brevity, we will not reproduce here the eight coupled Boltzmann equations which can instead be found in Ref. [54]. The code relies on the IsaReD[56] calculation of $\langle\sigma v\rangle(T)$ which is a crucial input to the coupled Boltzmann calculation. Starting from the time of re-heat with temperature $T_{R}$ at the end of the inflationary epoch, the computer code tracks the coupled abundances of radiation (i.e. SM particles), neutralinos, axinos, gravitinos, saxions and axions (the latter two consists of both thermal/decay-produced and coherent oscillation-produced (CO) components).

The CO-produced abundance of axions is determined in part by the axion field initial misalignment angle $\theta_{i}$ [24, 57]. For numerical analyses, we adopt a simple formula

[TABLE]

where $f(\theta_{i})=\left[\log\left(e/(1-\theta_{i}^{2}/\pi^{2})\right)\right]^{7/6}$ is the anharmonicity factor [57] and $N_{\rm DW}$ is the domain wall number ($=6$ for the DFSZ model). In previous work the initial misalignment angle $\theta_{i}$ is adjusted to gain the measured value of the relic abundance. In the current work, we will allow for a uniform distribution of $\theta_{i}:0-\pi$ values since we are scanning over many pocket universes which arise as subuniverses of the more vast multiverse.

In Fig. 2 we show the energy densities of various species vs. scale factor $R/R_{0}$ that influence the ultimate dark matter abundance for the landscape SUSY benchmark point landSUSY in Table 1. Here, $R_{0}$ is the reference scale factor at the beginning of re-heat and the corresponding temperature $T$ is shown by the dashed green line (where instead the $y$-axis is interpreted as temperature in GeV). We take $T_{R}=10^{7}$ GeV 333This value of $T_{R}$ is in accord with well-motivated baryogenesis mechanisms such as non-thermal or Affleck-Dine leptogenesis [58]. and $f_{a}=s_{0}=10^{12}$ GeV and where $s_{0}$ denotes the initial saxion field value. We also take $m_{\tilde{a}}=m_{s}=m_{3/2}=16$ TeV. The blue curve denotes the neutralino abundance which freezes out at $R\sim 10^{6}R_{0}$ or $T\sim 10$ GeV. The saxion and axion contributions are split into their thermally- and decay-produced components and their coherent-oscillation (CO) produced components. Saxions decay around $R\sim 10^{5}R_{0}$ ($T\sim 10$ GeV) whilst axinos decay around $R\sim 10^{6}R_{0}$ ( or $T\sim 1$ GeV). The saxion decays depend on a model dependent coupling $\xi_{s}$ which governs the saxion decay rate $s\rightarrow aa$ and $s\rightarrow\tilde{a}\tilde{a}$[54]. We take $\xi_{s}=1$ so these decays are turned on. (Of course, for our case the $s\rightarrow\tilde{a}\tilde{a}$ decay is not kinematically open so $s$ decays mainly to $aa$ but also to other MSSM particles). CO-produced axions (brown curve) start to oscillate around $T\sim 1$ GeV and become the dominant component of dark matter as one enters the era of entropy conservation on the right-hand-side of the plot. Due to late decays of axinos, which occur after neutralino freeze-out, the neutralino abundance increases to $\Omega_{\widetilde{\chi}^{0}_{1}}h^{2}\simeq 0.02.$

To gain some perspective on the expected relative abundances of mixed axion-WIMP dark matter, in Fig. 3 we show the relic density of mixed axion-WIMP dark matter vs. $f_{a}$ for the landSUSY benchmark point and with $T_{R}=10^{7}$ GeV and $m_{s}=m_{\tilde{a}}=m_{3/2}=16$ TeV and where $\theta_{i}=\theta_{s}=1$444Here, the saxion field strength $s=\theta_{s}.f_{a}$.. The green curve corresponds to the axion relic density while the blue curve corresponds to the WIMP relic density. The red curve shows the total relic density. We see that for low values of $f_{a}$, the axion relic density– arising here from coherent oscillations corresponding to Eq. 3– is highly suppressed. Also, the thermally-produced WIMP dark matter is highly suppressed due to the higgsino-like nature of the LSP which enhances its annihilation rate. The WIMP relic density is also highly suppressed by co-annihilations with the slightly heavier higgsinos $\widetilde{\chi}^{\pm}_{1}$ and $\widetilde{\chi}^{0}_{2}$. Thus, for $f_{a}\sim 10^{10}-10^{12}$ GeV, we expect typically an under-production of mixed axion-higgsino DM. As $f_{a}$ increases, the CO-produced axions steadily increase while WIMPs remain at their thermally-produced level. By $f_{a}\sim 10^{12}$ GeV, the axino and saxion decay rates are sufficiently suppressed (by $\Gamma_{\tilde{a},s}\sim 1/f_{a}^{2}$) that they begin decaying into higgsinos after WIMP freeze-out, thus augmenting the WIMP abundance with a non-thermal, decay-produced component. By $f_{a}\sim 3\times 10^{12}$ GeV, then the mixed axion-WIMP abundance saturates the measured value $\Omega_{\rm CDM}h^{2}\simeq 0.12$, and where at this point CDM consists nearly equally of axions along with a comparable thermal and non-thermal WIMP component. In this region, the non-thermal WIMP component arises mainly from thermal axino production followed by late $\tilde{a}$ decays in the early universe. As $f_{a}$ increases further, the thermal axino production rate falls off rapidly so that the WIMP abundance levels off. For even higher values of $f_{a}\gtrsim 10^{14}$ GeV, saxion production via COs becomes large and so saxion-decay produced WIMP production rapidly rises. In addition, the $s\rightarrow aa$ decays sharply increase the already over-produced axions. These relativistic axions also lead to violation of limits on relativistic species present in the early universe characterized in terms of the effective number of neutrinos parameter $\Delta N_{\nu}^{eff}$ which we take (very conservatively) to be $\lesssim 1$ (brown curve).555 In the Particle Data Book[22], it is tabulated that $N_{eff}=3.13\pm 0.32$. For $f_{a}\gtrsim 2\times 10^{15}$ GeV, then entropy dilution of all relics from CO-produced saxions can suppress the mixed axion-neutralino relic abundance.

In terms of the string theory landscape, we see that allowing values of $f_{a}\sim m_{\rm GUT}$ could lead to dark matter overproduction by a factor of $\sim 10^{4}$ compared to its measured value. As noted by Linde and others, it might be hard to visualize the existence of observers in a universe with such an overabundance of dark matter. Precisely how much of an overabundance of dark matter is anthropically too much is an open question. But clearly, if such a limit exists, then it would place an upper limit on the value of $f_{a}$. Even requiring a modest factor of four overabundance, indicated by the dashed gray horizontal line, would already require a value $f_{a}\lesssim 10^{13}$ GeV. This upper bound is well below the expected magnitude for $f_{a}$ from string theory where instead $f_{a}\sim 10^{16}-10^{18}$ GeV is typically expected[59]. The bound on $f_{a}$ from the axion abundance may be considered a softer bound since it is possible to lower the axion abundance with a smaller value of $\theta_{i}\sim 0$ (although if $\theta_{i}$ scans on the landscape, then $\theta_{i}\sim 1$ is to be expected). However, we see that a bound on $f_{a}$ still obtains from the WIMP contribution to $\Omega_{a\widetilde{\chi}^{0}_{1}}h^{2}$, although this bound on WIMP overproduction occurs at over an order of magnitude higher values: in Fig. 3, $f_{a}\lesssim 10^{14}$ GeV occurs from just overproduction of the WIMP component of dark matter.

In Fig. 4, we show the total mixed WIMP plus axion dark matter abundance but this time assuming $T_{R}=10^{6}$ GeV and $10^{8}$ GeV with $\theta_{i}=0.1,\ 0.5,\ 1$ and 2. For $T_{R}\gtrsim 10^{9}$ GeV, thermal production and late decay of gravitinos can lead to conflict with bounds from late-decaying neutral particles in the early universe: in this case, the gravitino problem[60, 61]. We see that for different $\theta_{i}$ values the upper limit on $f_{a}$ can move around by typically an order of magnitude: nonetheless, an upper bound on $f_{a}$ from overproduction of dark matter should obtain which is still much less than the the string/GUT scale. We also show variation in the $a-\widetilde{\chi}^{0}_{1}$ dark matter relic density versus varying $T_{R}$. For the DFSZ axion model, the axino and saxion production rates in the early universe hardly depend on $T_{R}$[62] (unlike the case of the KSVZ axion model[63]). Some variation in relic density is seen for $f_{a}\gtrsim 10^{14}$ GeV where gravitino production, which does depend on $T_{R}$[64], becomes important and augments the non-thermal WIMP abundance.

2.3 PQ scale from the landscape

In this Subsection, we investigate whether landscape considerations can determine the magnitude of the PQ scale $f_{a}$. We assume an $n=1$ statistical draw towards large soft terms $-A_{f}$ which in turn leads to large PQ scales along the lines of Fig. 1 where the PQ scale is related to the breakdown of supersymmetry. For our landscape benchmark point landSUSY, the magnitude of $f_{a}$ is determined by the quartic soft term $A_{f}$. However, since $-A_{f}$ is not connected with EWSB, then it need not be susceptible to the same bounds on MSSM soft terms that emerge from requiring an appropriate breakdown of electroweak symmetry with independent contributions to $m_{weak}$ not more than a factor of a few from its value $m_{weak}\simeq 100$ GeV. Instead, the PQ scale $f_{a}$ is intimately related to the production of both axion dark matter and (natural) higgsino-like WIMP dark matter.

Since we are working within a multiverse scenario wherein each pocket universe may have different laws of physics, and the multiverse is an expression of the universe emerging from a spacetime continuum characterized by eternal inflation, then of course inflationary cosmology is an essential component of our overall scheme. In inflationary cosmology, the universe has an early exponential expansion phase which drives the universe to flatness, which requires an overall energy density teetering on the boundary between an open or a closed universe. Such a universe is characterized by the overall energy density lying at its critical closure density:

[TABLE]

and where $\Omega_{curv}=0$ for an inflationary universe which gives rise to a flat geometry. For our pocket universe, the measured value of the Hubble constant is $H_{0}=100h$ km/s/Mpc with $h=0.678\pm 0.009$ but for other pocket universes then $H_{0}$ will be different depending on the various constituencies. We will adopt as usual $\rho_{B}/\rho_{\gamma}$ equal to the value of our universe since we are assuming a “friendly” fertile patch of the multiverse where the SM remains as the low energy effective theory.666 Anthropic arguments usually depend on a so-called “friendly” landscape wherein one focuses on most parameters asuming their SM values so as to retain predictivity[65]. Sometimes these are called fertile patches of the landscape of vacua since they should lead to the standard cosmological and particle physics models aside from just the few mass scales which may scan in the multiverse. Thus, in the fertile patch of multiverse assumed here, only $\Lambda$, $m_{soft}$ and $\theta_{i}$ are assumed to scan. The scanning of the soft term $A_{f}$ sets the value of $f_{a}$ for otherwise fixed values of scalar masses as expected for our landSUSY benchmark point: i.e. we assume a common value of all scalar masses $m_{0}(1,2)=m_{X}=m_{Y}\equiv m_{3/2}$. We allow smaller values of $m_{0}(3)$ as occurs in the mini-landscape picture wherein third generation fields lie on the bulk of the compactified orbifold whilst first/second generation fields lie near orbifold fixed points[66].

The generated probability distribution for $-A_{f}$ is shown in Fig. 5a), which is seen to rise linearly as expected. For a given value of $A_{f}$, then the value of $f_{a}$ is determined by the minimization conditions arising from Eq. 2. In frame b), we show the derived distribution $dP/df_{a}$. Here, the probability distribution is seen to favor the highest values of $f_{a}$ possible, which would be generated from very large values of $-A_{f}$.

At this point, our prior distribution for $f_{a}$ is set, but we will also need some selection criterion to avoid $f_{a}$ exploding up to huge values, leading to perhaps a gross overproduction of dark matter. Thus, the question now is: how much dark matter is too much dark matter for our fertile patch of pocket universes within the greater multiverse? Some of the papers of Subsec. 1.1 have entertained values of $\rho_{DM}/\rho_{B}$ as high as 25-100.

For illustrative purposes, we will consider the effect of limiting pocket universes to a modest bound of four times greater dark matter density than in our universe: suppose $\Omega_{\rm DM}h^{2}\lesssim 0.48$. Such a bound would saturate the case where we maintain our measured value of $\rho_{c}$ but allow the dark matter abundance to nearly saturate $\rho_{c}$ at the expense of a dark energy component. Such models were commonly contemplated before the discovery of a non-zero dark energy component.

In Fig. 6(a)), we show the resulting probability distribution $dP/df_{a}$ which results from an $n=1$ draw on $-A_{f}$ coupled to an anthropic/cosmological selection bound $\Omega_{a\widetilde{\chi}^{0}_{1}}h^{2}<0.48$ (green curve). Even with our proposed modest selection bound, we see that the value of $f_{a}$ is driven to its nearly maximal value such as to avoid overproduction of dark matter. From the plot, we would expect that a value of $f_{a}\sim 10^{14}$ GeV or only somewhat lower, with a rather sharp cutoff $f_{a}\lesssim 8\times 10^{13}$ GeV. For comparison, we also show the black histogram where we instead require that the upper bound on dark matter abundance is only slightly beyond our measured value: $\Omega_{\rm DM}<0.15$. This case would prefer $f_{a}\sim 5\times 10^{12}$ GeV.

Let us compare the results of Fig. 6(a)) with those of Fig. 7 which shows the allowed mixed axion-WIMP dark matter abundance for our landSUSY benchmark point in the generic SUSY DFSZ axion model while scanning uniformly over $\theta_{i}$ and uniformly over $\log(f_{a})$. From Fig. 7, we see that for $f_{a}\sim 10^{13}-10^{14}$ GeV, we are already overproducing dark matter compared to our universe with $\Omega_{\rm DM}h^{2}=0.12$. There is only a miniscule probability to obtain from Fig. 6(a)) $f_{a}$ values low enough to match the measured value, which occurs for $f_{a}\sim 10^{11}$ – $\sim 4\times 10^{12}$ GeV. In the references from Subsec. 1.1, large values of $f_{a}\sim 10^{14}-10^{16}$ GeV could be compensated for by selecting on small values of $\theta_{i}$. For our case of natural mixed axion-WIMP dark matter, this compensation is not permitted because large $f_{a}$ also leads to large (non-thermal) overproduction of WIMP dark matter via delayed axino and saxion decays in the early universe. From Fig. 6(a), we would expect that if the landscape is involved in determining the PQ scale $f_{a}$, then its value should be very near the maximally allowed abundance of DM in pocket universes such as to allow observers to exist. But it is hard to believe that our pocket universe’s value of dark matter abundance is nearly anthropically maximal (as depicted by the black curve of Fig. 6(a)).

In Fig. 6(b)), we show the corresponding distribution $dP/d\theta_{i}$ from the $n=1$ pull on soft terms coupled with our modest anthropic veto that $\Omega_{\rm DM}h^{2}<0.48$. The plot shows a probability that $\theta_{i}$ is peaked around its smallest allowed values. This is easy to understand in that while the landscape prior strongly favors large values of $f_{a}$, from Eq. 3 we see that overproduction of axions can be avoided by selecting only those vacua with correspondingly tiny values of $\theta_{i}$. This effect is easily understood from Fig. 8 where we show regions of the $\theta_{i}$ vs. $f_{a}$ plane for our landSUSY benchmark point which lead to $\Omega_{\rm DM}h^{2}<0.48$ (green points) or $\Omega_{\rm DM}h^{2}>0.48$ (yellow points). The brown points denote where also $\Delta N_{eff}>1$. From the figure, we see that for large $f_{a}\sim 8\times 10^{13}$ GeV, only a small range of $\theta_{i}$ allows for non-overproduction of dark matter. And once $f_{a}\gtrsim 8\times 10^{13}$ GeV, then no value of $\theta_{i}$ is possible which allows one to avoid DM overproduction.

3 Prediction of PQ scale from generic SUSY DFSZ

axion model with uniform scan on $\theta_{i}$

In Sec. 2, we adopted a particular gravity-safe SUSY axion model based on a $\mathbb{Z}_{24}^{R}$ discrete $R$-symmetry and the hyCCK superpotential Eq. 1 to show that a statistical draw towards large soft terms also yields a draw to large PQ breaking scale $f_{a}$. The value of $f_{a}$ gains an upper bound by requiring no overproduction of dark matter. For the modest assumption of less than a factor four times the measured abundance of dark matter, then we found $f_{a}\sim 10^{14}$ GeV which is well below the values expected from pre-landscape string theory but which typically leads to much more dark matter production than we observe in our universe.

In this Section, we try to be more general by eschewing a particular SUSY axion model and instead assume a generic SUSY DFSZ axion model[54] where $f_{a}$ is an input instead of an output parameter. In this case, we will adopt a uniform distribution in $\theta_{i}$ in accord with expectations from the landscape, but then require that the dark matter abundance lie at its measured value: $\Omega_{a\widetilde{\chi}^{0}_{1}}h^{2}=0.12$. From this, we can then determine the necessary value of $f_{a}$ such that, for scanned values of $m_{\tilde{a}}$, $m_{s}$ and $m_{3/2}$, the measured abundance of mixed axion-neutralino dark matter is obtained. We will scan uniformly over each of $m_{\tilde{a}},\ m_{s}$ and $m_{3/2}:1-50$ TeV.

For a SUSY benchmark point within a two-component dark matter framework, direct and indirect WIMP dark matter searches can put a stringent upper limits upon the neutralino density which are more severe than the measured value, $\Omega_{\rm DM}h^{2}=0.12$. In many cases, indirect DM detection (IDD) offers the most contraining limits on the non-thermal, decay-produced neutralinos for models with thermally underproduced higgsino-like neutralinos. In natural SUSY models from the $n=1$ landscape, thermally produced neutralinos typically make up 5-20% of the total CDM density which renders them safe from Fermi-LAT+MAGIC[67] limits on overproduction of gamma rays in dwarf spheroidal galaxies[68]. In Fig. 9, we show the allowed percentage of WIMP dark matter compared to $m_{\widetilde{\chi}^{0}_{1}}$ along with the Fermi-LAT+MAGIC IDD limit for our landSUSY benchmark point. If we increase $\mu\sim 340$ GeV, then $m_{\widetilde{\chi}^{0}_{1}}\sim 340$ GeV and all generated points would be Fermi-LAT+MAGIC allowed. The gray-shaded region shows the excluded WIMP composition for all landSUSY points within a good approximation.

In Fig. 10, we work within the SUSY DFSZ axion model using again our landSUSY benchmark point but with input parameters $m_{s},\ \theta_{s},\ \theta_{i},\ f_{a},\ T_{R}$ and $m_{3/2}$. Here, we fix $T_{R}=10^{7}$ GeV and $\theta_{s}=1$ but allow $m_{s}$, $m_{\tilde{a}}$ and $m_{3/2}$ to scan over the range given above with a uniform scan on $\theta_{i}$ and a log prior scan on $f_{a}$. We only accept solutions with $\Omega_{a\widetilde{\chi}^{0}_{1}}h^{2}=0.12$. We show the parameter space with augmented neutralino densities $\Omega_{\widetilde{\chi}^{0}_{1}}h^{2}<0.12,0.06$ and 0.03 with black, orange and purple colors respectively. We impose an upper limit on $\theta_{i}$ ($\theta_{i}<3.14$) so that the highly fine-tuned region $\theta_{i}\simeq\pi$ is not present in our analysis.

In Fig. 10 frame a), the resulting abundance of neutralino dark matter is shown while the remainder of DM is made of DFSZ axions. The horizontal line around $f_{a}\lesssim 10^{11}$ GeV is just the expected thermal abundance of 200 GeV higgsino-like WIMP dark matter. For higher $f_{a}>10^{11}$ GeV, then non-thermal LSP production begins to occur where axinos can be produced in the early universe and decay to LSPs after neutralino freeze-out. There is a gap around $f_{a}\sim 10^{13}$ GeV where axino decays are still contributing to the neutralino density. For $f_{a}\sim 10^{14}$ GeV, points again become allowed due to diminished thermal production of axinos in the early universe. The WIMP abundance increases for $f_{a}\gtrsim 10^{14}$ GeV due to increasing CO-production of saxions which then decay (in part) to WIMPs[54]. An upper limit of $f_{a}\lesssim 2\times 10^{14}$ GeV ensues in this case since for large $f_{a}$, the DM is always overproduced. There is no conflict here with Fig. 7 since in this case with random values of $m_{s}$ and $m_{\tilde{a}}$, then relative axino and saxion production and decay rates can vary which leads to allowed points for $f_{a}\sim 10^{14}$ GeV.

The corresponding correlation of $\theta_{i}$ with the required value of $f_{a}$ to make $\Omega_{a\widetilde{\chi}^{0}_{1}}=0.12$ is shown in frame b). Here, large values of $\theta_{i}$ are correlated with low values of $f_{a}$ to boost the axion production to gain accord with the measured relic abundance. For very large $f_{a}$, then consequently small values of $\theta_{i}$ are required to allow for $\Omega_{a\widetilde{\chi}^{0}_{1}}h^{2}=0.12$.

In frame c), we show the resulting probability distribution $dP/df_{a}$ versus $f_{a}$. Including all points with the measured abundance, then one obtains the black histogram which peaks around $f_{a}\sim 2\times 10^{12}$ GeV but with a tail extending to over $10^{14}$ GeV. For the cases in which neutralino makes less than half of the measured DM density, the peak shifts to lower values of $f_{a}$ (orange and purple histograms) with a small probability at high $f_{a}$.

In frame d), we show the probability distribution $dP/d\theta_{i}$ for the three cases considered. Here, the black histogram is almost uniform across its range whilst the orange histogram displays a gap at small $\theta_{i}$ where no allowed solutions occur. The high $f_{a}$ region does not show up for $\Omega_{\widetilde{\chi}^{0}_{1}}h^{2}\lesssim 0.04$ and $\theta_{i}$ can only take values greater than $\sim$1 when the neutralino makes less than 25% of the total DM abundance (purple).

4 Conclusions

In this paper we have sought to answer the question: is the magnitude of the PQ scale $f_{a}$ set by the landscape, or by something else? To address this question, we have adopted the scenario advocated by Douglas wherein the soft terms are statistically favored by a prior distribution $m_{soft}^{2n_{F}+n_{D}-1}$ and where we take the value $n=2n_{F}+n_{D}-1=1$ (i.e. a linear distribution favoring large soft SUSY breaking terms). Along with this prior distribution, we invoke a selection criteria that vetos models with inappropriate EW breaking (CCB minima or no EWSB) and vetos models with contributions to the weak scale $\gtrsim 4$ (corresponding to $\Delta_{\rm EW}>30$) in accord with nuclear physics constraints derived by Agrawal et al. on anthropically allowed values for the weak scale. Such an approach receives support in that previously it has been shown that $n=1$ (or 2) has a most probable Higgs mass of $m_{h}\simeq 125$ GeV whilst lifting sparticle masses beyond the reach of Run 2 of the LHC.

We implement this approach within a highly motivated SUSY axion model labeled as hyCCK. We feel this is an improvement upon previous work in that the model

1. stabilizes the weak scale via SUSY, 2. solves the strong CP problem via a gravity-safe $\mathbb{Z}_{24}^{R}$ discrete symmetry which gives rise to an accidental, approximate global $U(1)_{PQ}$, 3. solves the SUSY $\mu$ problem[69] via the presence of a Kim-Nilles superpotential[70] and 4. explains $R$-parity conservation and proton stability as consequences of the more fundamental $\mathbb{Z}_{24}^{R}$ which could arise from compactification of extra space dimensions in string theory. By choosing a MSSM benchmark point in accord with $n=1$ landscape predictions and allowing for the PQ soft term $-A_{f}$ to scan linearly and to set the magnitude of the PQ scale $f_{a}$, then we find that as large as possible values of $f_{a}$ are statistically prefered. In this approach, typically both WIMP and axion dark matter are overproduced. Axions are overproduced via vacuum misalignment at $f_{a}\gg 10^{12}$ GeV unless the large value of $f_{a}$ is compensated for by a small value of $\theta_{i}$. However, since WIMPs are also overproduced at large $f_{a}$, due to axino and saxion production coupled with delayed decays to SUSY LSPs after neutralino freeze-out (i.e. non-thermal WIMP production), then even with small $\theta_{i}$ one cannot avoid overproduction of mixed axion-WIMP dark matter. From this approach, using even a modest DM overproduction bound of a factor four, then there is only a tiny probability to gain the measured value of dark matter density in our own pocket universe. Thus, the answer to the question posed in the title is: No, in our well-motivated landscape SUSY model based upon gravity-safe, electroweak natural hyCCK SUSY axion model, the magnitude of the PQ scale is highly unlikely to be set by the landscape.

Instead, an alternative but perhaps underappreciated mechanism is available to set the magnitude of the PQ scale. This is that in generic supergravity models with hidden sector SUGRA breaking via the superHiggs mechanism, then soft terms arise from SUGRA breaking with magnitudes of order the gravitino mass $m_{3/2}$. For a well-specified hidden sector, then the soft terms are all calculable and correlated. For our landscape SUSY model with $m_{0}(1,2)\sim m_{3/2}$, we would also expect $-A_{f}\sim m_{3/2}\sim 10-100$ TeV. This places us from Fig. 1 into the zone where $f_{a}\sim 10^{11}-10^{12}$ GeV which is the sweet spot for generating a thermal underabundance of higgsino-like WIMP dark matter but with mainly SUSY DFSZ axion dark matter.

In Sec. 3, we implemented instead a uniform scan over $\theta_{i}$ and a scan over independent $m_{\tilde{a}}$, $m_{s}$ and $m_{3/2}$ values with the corresponding $f_{a}$ value such that $\Omega_{a\widetilde{\chi}^{0}_{1}}h^{2}=0.12$. In this case, we are pushed to the cosmological sweet spot where $f_{a}\simeq 10^{11}-10^{13}$ GeV.

From our overall approach, we understand why WIMPs haven’t yet been detected: it is because they make up only a small portion of the total dark matter. In addition, even though the non-thermal decay-production of WIMPs is allowed for, more than 40% of the points in Fig. 10a) (Sec. 3) still have the same neutralino abundance as the thermally-produced value. Nonetheless, WIMP discovery should be possible at multi-ton noble liquid direct WIMP detection experiments[71]. Axions– while likely to make up the bulk of dark matter– are much more difficult to detect since the presence of higgsinos in their $g_{a\gamma\gamma}$ coupling diagram suppresses their detection rate compared to KSVZ or non-SUSY DFSZ axion models[68]. Meanwhile, sparticle detection may be possible at HL-LHC via the soft OS dilepton+jet+MET channel which arises from direct higgsino pair production[16]. Detection of gluinos, top squarks or winos might be possible at HL-LHC if we are lucky, but otherwise may have to await construction of higher energy $pp$ colliders[14].

Acknowledgements: This work was supported in part by the US Department of Energy, Office of High Energy Physics. The computing for this project was performed at the OU Supercomputing Center for Education & Research (OSCER) at the University of Oklahoma (OU).