(Mainly) axion dark matter^{1}^{1}1Plenary talk given at the Particle Physics and Cosmology 2015 (PPC2015) meeting, Deadwood, SD, June 29, 2015
Abstract
The strong CP problem of QCD is at heart a problem of naturalness: why is the 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 PecceiQuinn (PQ) symmetry which requires the existence of the axion field . 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 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 higgsinolike WIMPs so then we expect mixed axionWIMP 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.
aff1]Howard Baer
Introduction: the strong CP problem and axions
The story of the axion begins at the dawn of QCD, circa mid1970s[1]. With two light quarks, one expects QCD to manifest an approximate chiral symmetry which can be recast as . The symmetry gives rise to the wellknown isospin and baryon number symmetries while the axial is spontaneously broken. Since is rank 4, we expect four pseudoGoldstone bosons– the pions– whilst we see only three with . Weinberg dubbed this conundrum the problem and suggested that somehow nature doesn’t respect the symmetry[2]. Shortly thereafter, ’t Hooft discovered the QCD vacuum and the effect of instantons[3]. Indeed, the ground state of QCD did not respect the symmetry and was explained. The price to pay was that the QCD Lagrangian should contain a term of the form
(0) 
In addition, a complex quark mass matrix gives rise to a second contribution so that . This term violates symmetry and leads to a measurable contribution to the neutron EDM. Measurements of the neutron EDM imply . Somehow the term must be very tiny. This is contrary to the old maxim of quantum mechanics: everything is allowed unless explicitly forbidden. How to suppress or get rid of the awkward Lagrangian term Eq. Introduction: the strong CP problem and axions constitutes the famous strong CP problem.
While a variety of solutions to the strong CP problem have been proposed, one stands out for its simplicity and elegance: the PecceiQuinn (PQ) solution[4]. PQ proposed the existence of a new global symmetry which is spontaneously broken at the PQ energy scale . A new Goldstone field, the axion, arises in accord with Goldstone’s theorem. The QCD Lagrangian now includes
(0) 
where the latter term provides a potential for the axion field:
(0) 
The axion field settles dynamically to its minimum at and the offending term goes to zero: the strong problem is solved! As a consequence of this solution, a physical axion particle should exist with mass[5] where . In the original proposal PQ suggested the scale , a value which was soon ruled out by reactor and beam dump experiments. KSVZ[6] and DFSZ[7] proposed variant models with much higher values of GeV, a value which suppressed anomalous decays leading to the socalled invisible axion. Sikivie stripped off the invisibility cloak by suggesting axion detection in microwave cavity experiments[8].
Axion dark matter: relic density and detection
In spite of the fact that axions are extremely weakly coupled– their coupling strength is suppressed by – they still can play the role of cold dark matter. The equation of motion for the axion field in the early universe is given by
(0) 
with . This is the equation of a damped harmonic oscillator. At high temperatures, and the solution is that a constant. The axion mass turns on for GeV and the axion field begins coherent oscillations for which the equation of state is that of cold dark matter. If PQ symmetry breaks after the end of inflation, then must be averaged over disparate domains leading to a definite prediction of shown by the blue line in Fig. 1. If PQ symmetry breaks before the end of inflation, then the axion relic density computed from this “vacuum misalignment” mechanism is given by[9]
(0) 
so that the measured relic density can be achieved for any value by an appropriate choice of . Values of GeV are disallowed by stellar cooling bounds[10] while large values of Hubble constant at the end of inflation are disallowed by axion isocurvature limits[11]: see Fig. 1 for a panoramic view.
The region of Fig. 1 around GeV has been explored by the ADMX experiment. ADMX makes use of a microwave cavity which can be tuned over a range of frequencies. If the frequency is just right, then resonance occurs and a bump should appear in the spectra. ADMX has been upgraded with a new dilution refrigerator and SQUID electronics and hopes to explore more broadly and deeply in vs. parameter space in the near future.
Expectations for an axion signal are model dependent. Two commonly used models include KSVZ which introduced intermediate mass scale PQ charged heavy quark fields or the DFSZ model wherein two PQ charged Higgs doublets are required which couple to a PQ charged but SM singlet field . The ADMX experiment has reached sensitivity to axions in the KSVZ model but at present is short of the DFSZ coupling strength.
Supersymmetric axions: axinos and saxions
While much of the previous discussion is based on longknown results, the biggest recent development in axion physics is…. the discovery of the Higgs boson with GeV[13, 14]. The issue here is that the Higgs boson exists and appears very much SMlike and as a fundamental scalar field. If one then expands the SM to include axions via e.g. the KSVZ or DFSZ model, then one would expect the Higgs mass to blow up to scale. Of course, one can always finetune to but such tuning is generally symptomatic of some missing ingredient in the model.^{2}^{2}2Weinberg states: “The appearance of finetuning in a scientific theory is like a cry of distress from nature, complaining that something needs to be better explained”[15].
This electroweak finetuning crisis is elegantly solved via the introduction of supersymmetry where the offending quadratic divergences to scalar field masses all neatly cancel. Weak scale SUSY (SUSY with weak scale soft breaking) is supported experimentally by 1. the measured strengths of the gauge couplings which neatly unify within the MSSM[16], 2. the large mass of the top quark which is needed for radiative EWSB[17] and 3. the measured value of which squarely sits in the predicted SUSY window GeV[18].
Nonetheless, it has alarmed many theorists that SUSY matter states have not be found at LHC8 searches. However, many other theorists expected a rather high value of soft breaking scale since it offers a decoupling solution to the SUSY flavor and CP problems and is in accord with the gravitino problem for gravitino mass TeV[19]. The enigma is the connection of the SUSY breaking scale to the weak scale: GeV TeV. This is the emerging Little Hierarchy problem. It requires a scrutinization of naturalness measures.
While simple evaluations of large logs or BarbieriGiudice naturalness pointed to high finetuning, the methodology used in these measures has been criticized[20, 21, 22, 23, 24] in that they neglect dependent contributions of oppositesign which lead to large cancellations. A proper evaluation of finetuning allows for TeVscale highly mixed 3rd generation squarks but does require light higgsinos of mass GeV, the lighter the better. The lightest higgsino is the LSP but is thermally underproduced as WIMP dark matter with typically . The light higgsino spectrum is compressed with interhiggsino mass gaps of GeV so that higgsino decay products are soft and are buried under QCD backgrounds at LHC[25]. The only soft SUSY breaking term required to be near is . This term can be driven radiatively to small instead of large negative values at the weak scale in models with nonuniversal Higgs soft masses at the GUT scale. Such models are those with radiativelydriven naturalness or RNS[26].
By requiring SUSY to solve the EW naturalness problem and PQ symmetry to solve the QCD naturalness (strong CP) problem, then the axion becomes but one element of the axion superfield. It is now accompanied by the spin0 parityeven saxion and the spin parityodd axino . In gravitymediated SUSY breaking models, it is expected that [27]. In such a case, then dark matter is composed of two particles: the SUSY WIMP and the axion.
The presence of axions in weak scale SUSY models offers an important and elegant solution to the SUSY mu problem. Since the term occurs in the superpotential (it is supersymmetric and not SUSY breaking) one expects naively its value to be (the Planck scale) while phenomenology/naturalness require GeV. Kim and Nilles recognized[28] that the supersymmetrized DFSZ axion model offers a solution to the mu problem. The mu term is first forbidden because it violates PQ symmetry but then it is regenerated via the DFSZ coupling where the PQ field develops a vev under PQ breaking. Then
(0) 
while . The Little Hierarchy with is just a reflection of a mismatch between PQ breaking scale and hidden sector mass scale . Such a mismatch can arise in models such as the MSY model of radiativelydriven PQ breaking[29]. In the MSY model, SUSY breaking effects drive one of the PQ fields to negative squared mass causing PQ symmetry to break much as EW symmetry is broken radiatively due to the large top quark Yukawa coupling. Minimizing the PQ scalar potential, then one typically finds GeV is induced by TeV[30]: the mu problem is solved and the Little Hierarchy is explained! The PQ breaking scale sets the mass for the axion, the Higgs and the higgsinos! In addition, Majorana masses are induced for righthand neutrinos. While this behavior is exhibited for the MSY model, it seems typical of a much larger class of models.
Signals at LHC and ILC
SUSY models with radiatively driven naturalness and with a supersymmetrized DFSZ axion solution are natural in both the EW and QCD sectors, can accommodate GeV with Tevscale highly mixed stops and can evade LHC8 searches. In contrast to the usual expectation that , in this case only – the other sparticles can be much heavier, TeV. Upper bounds on sparticle masses are computed in Ref. [26, 31]. For natural solutions with , then TeV, within the reach of highluminosity LHC13[32]. In this case, gluino cascade decay events should contain a characteristic dilepton mass edge[32, 33] with GeV in accord with the interhiggsino mass gap. Also, a unique samesign diboson signature sans hard jet activity should arise from wino pair production[34, 32] where and . For (as defined in e.g. Ref. [26]), then can range up to 4 TeV and lie beyond LHC13 reach. Monojet plus dilepton signals offer another LHC detection possibility[35, 36].^{3}^{3}3For more on monojet signals, see also Ref’s [chan, 38, 39].
The smoking gun signature of RNS SUSY will be direct higgsino pair production at an collider such as ILC[40] which would operate with . Built originally as a higgs factory, ILC will turn out to be a higgsino factory. SUSY can be discovered (or an LHC13 discovery can be confirmed) and precision measurements can be made which test both the higgsino and gaugino sectors.
Signals at axion and WIMP detectors
In SUSY with radiativelydriven naturalness, one expects a mixture of axion plus higgsinolike WIMP dark matter. The computation of the dark matter relic abundance requires the solution of eight coupled Boltzmann equations which track the radiation, neutralino, axino,gravitino, saxion and axion (both CO and THproduction) contributions. Axino production and decay can feed into and augment the neutralino abundance. Saxion production and decay can feed the neutralino abundance or dilute it via decay to radiation; it can also inject dark radiation via decays.^{4}^{4}4See also Ref. [41] If too many WIMPs are produced from axino or saxion decays, then they may reannihilate at the particle decay temperature[42]. The calculation for the SUSY DFSZ model[43] is given in Ref. [44] and Fig. 2 for the DFSZ axion in natural SUSY. For low , higgsinos are underproduced and the relic abundance is axiondominated[45]. For higher values, the axino and saxion decay later and increase the WIMP abundance. For too large , then WIMPs become overproduced and the model becomes excluded.
What of WIMP detection? In models with mixed axionWIMP dark matter, then WIMPs make up only part of the relic abundance so the assumed local abundance of WIMPs must be scaled down by a factor . The results of rescaled higgsinolike WIMP SI direct detection rates are shown in Fig. 3. In spite of the rescaling, tonscale noble liquid detectors should make a complete exploration of the expected parameter space. While prospects for direct WIMP detection are excellent, rates for indirect detection must be rescaled by a factor of (for WIMPWIMP annihilation in the galactic halo) or by (for IceCube searches); these factors typically suppress detection rates to very low levels[47, 46].
In addition to WIMP detection, we also expect detection of a DFSZlike axion. The range of that can be explored by ADMX and successor experiments is shown in Fig. 4.
Bullet point summary

Axion dark matter is a byproduct of the elegant PQWW/KSVZ/DFSZ solution to the strong CP problem.

Axions need SUSY or else we would expect .

Electroweak naturalness requires the SUSY mu term GeV, the lighter the better.

The SUSY DFSZ axion model allows for a solution to the SUSY mu problem.

The Little Hierarchy may be a reflection of .

Small GeV can be generated from large in gravitymediated SUSY breaking models with radiative breaking of PQ symmetry.

LHC may discover natural SUSY but discovery is not guaranteed in SUSY with radiativelydriven naturalness.

The ILC would make a guaranteed search for light higgsinos. Upon discovery, precision measurements of their properties are possible.

In such models, we expect mixed axionhiggsinolike WIMP dark matter.

Ultimately, we expect detection of both an axion and a WIMP.
Acknowledgments
I thank CETUP and B. Szczerbinska for their kind hospitality. I thank my collaborators Kyu Jung Bae, Vernon Barger, K. Y. Choi, J. E. Kim, Andre Lessa, Dan Mickelson, Azar Mustafayev, Maren PadeffkeKirkland, Mike Savoy, Hasan Serce and Xerxes Tata. This research was funded in part by the US Department of Energy Office of High Energy Physics.
References
 [1] R. D. Peccei, AIP Conf. Proc. 1274 (2010) 7 [arXiv:1005.0643 [hepph]].
 [2] S. Weinberg, Phys. Rev. D 11, 3583 (1975).
 [3] G. ’t Hooft, Phys. Rev. Lett. 37, 8 (1976).
 [4] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977); S. Weinberg, Phys. Rev. Lett. 40, 223 (1978); F. Wilczek, Phys. Rev. lett. 40, 279 (1978).
 [5] W. A. Bardeen and S.H. H. Tye, Phys. Lett. B 74 (1978) 229.
 [6] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); M. A. Shifman, A. Vainstein and V. I. Zakharov, Nucl. Phys. B 166, 493(1980).
 [7] M. Dine, W. Fischler and M. Srednicki, Phys. Lett. B 104, 199 (1981); A. P. Zhitnitskii, Sov. J. Phys. 31, 260 (1980).
 [8] P. Sikivie, Phys. Rev. Lett. 51 (1983) 1415.
 [9] L. F. Abbott and P. Sikivie, Phys. Lett. B 120, 133 (1983); J. Preskill, M. Wise and F. Wilczek, Phys. Lett. B 120,127 (1983); M. Dine and W. Fischler, Phys. Lett. 120, 137 (1983); M. Turner, Phys. Rev. D 33, 889 (1986); K. J. Bae, J. H. Huh and J. E. Kim, JCAP 0809 (2008) 005; L. Visinelli and P. Gondolo, Phys. Rev. D 80 (2009) 035024.
 [10] For a review, see G. G. Raffelt, J. Phys. A 40 (2007) 6607.
 [11] M. P. Hertzberg, M. Tegmark and F. Wilczek, Phys. Rev. D 78 (2008) 083507.
 [12] H. Baer, K. Y. Choi, J. E. Kim and L. Roszkowski, Phys. Rept. 555 (2014) 1.
 [13] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B716 (2012) 1.
 [14] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B716 (2012) 30.
 [15] S. Weinberg in To Explain the World: The Discovery of Modern Science (Harper, 2015, 416p).
 [16] U. Amaldi, W. de Boer and H. Furstenau, Phys. Lett. B 260, 447 (1991); J. R. Ellis, S. Kelley and D. V. Nanopoulos, Phys. Lett. B 260 (1991) 131; P. Langacker and M. x. Luo, Phys. Rev. D 44 (1991) 817.
 [17] L. E. Ibañez and G. G. Ross, Phys. Lett. B110, 215 (1982); K. Inoue et al. Prog. Theor. Phys. 68, 927 (1982) and 71, 413 (1984); L. Ibañez, Phys. Lett. B118, 73 (1982); H. P. Nilles, M. Srednicki and D. Wyler, Phys. Lett. B 120 (1983) 346; J. Ellis, J. Hagelin, D. Nanopoulos and M. Tamvakis, Phys. Lett. B125, 275 (1983); L. AlvarezGaumé. J. Polchinski and M. Wise, Nucl. Phys. B221, 495 (1983); B. A. Ovrut and S. Raby, Phys. Lett. B 130 (1983) 277; for a review, see L. E. Ibanez and G. G. Ross, Comptes Rendus Physique 8 (2007) 1013.
 [18] M. S. Carena and H. E. Haber, Prog. Part. Nucl. Phys. 50 (2003) 63.
 [19] M. Dine, A. Kagan and S. Samuel, Phys. Lett. B 243 (1990) 250; A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Lett. B 388 (1996) 588; J. Bagger, J. L. Feng and N. Polonsky, Nucl. Phys. B 563 (1999) 3; T. Moroi and M. Nagai, Phys. Lett. B 723 (2013) 107.
 [20] H. Baer, V. Barger and D. Mickelson, Phys. Rev. D 88 (2013) 9, 095013.
 [21] H. Baer, V. Barger, D. Mickelson and M. PadeffkeKirkland, Phys. Rev. D 89 (2014) 115019.
 [22] A. Mustafayev and X. Tata, Indian J. Phys. 88 (2014) 991.
 [23] H. Baer, V. Barger and M. Savoy, Phys. Scripta 90 (2015) 068003.
 [24] X. Tata, Phys. Scripta 90 (2015) 10, 108001.
 [25] H. Baer, V. Barger and P. Huang, JHEP 1111 (2011) 031.
 [26] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev and X. Tata, Phys. Rev. D 87 (2013) 11, 115028.
 [27] E. J. Chun and A. Lukas Phys. Lett. B 357, 43 (1995); J. E. Kim and M. S. Seo, Nucl. Phys. B 864, 296 (2012).
 [28] J. E. Kim and H. P. Nilles, Phys. Lett. B 138, 150 (1984); E. J. Chun, J. E. Kim and H. P. Nilles, Nucl. Phys. B 370, 105 (1992).
 [29] H. Murayama, H. Suzuki and T. Yanagida, Phys. Lett. B 291 (1992) 418.
 [30] K. J. Bae, H. Baer and H. Serce, Phys. Rev. D 91, 015003 (2015) [arXiv:1410.7500 [hepph]].
 [31] H. Baer, V. Barger and M. Savoy, arXiv:1509.02929 [hepph].
 [32] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev, W. Sreethawong and X. Tata, JHEP 1312 (2013) 013.
 [33] B. Altunkaynak, H. Baer, V. Barger and P. Huang, Phys. Rev. D 92 (2015) 3, 035015.
 [34] H. Baer, V. Barger, P. Huang, D. Mickelson, A. Mustafayev, W. Sreethawong and X. Tata, Phys. Rev. Lett. 110 (2013) 15, 151801.
 [35] Z. Han, G. D. Kribs, A. Martin and A. Menon, Phys. Rev. D 89 (2014) 7, 075007.
 [36] H. Baer, A. Mustafayev and X. Tata, Phys. Rev. D 90 (2014) 11, 115007.
 [37] C. Han, A. Kobakhidze, N. Liu, A. Saavedra, L. Wu and J. M. Yang, JHEP 1402 (2014) 049.
 [38] H. Baer, A. Mustafayev and X. Tata, Phys. Rev. D 89 (2014) 5, 055007.
 [39] D. Barducci, A. Belyaev, A. K. M. Bharucha, W. Porod and V. Sanz, JHEP 1507 (2015) 066.
 [40] H. Baer, V. Barger, D. Mickelson, A. Mustafayev and X. Tata, JHEP 1406 (2014) 172.
 [41] F. S. Queiroz, K. Sinha and W. Wester, Phys. Rev. D 90 (2014) 11, 115009.
 [42] K. Y. Choi, J. E. Kim, H. M. Lee and O. Seto, Phys. Rev. D 77 (2008) 123501; H. Baer, A. Lessa, S. Rajagopalan and W. Sreethawong, JCAP 1106 (2011) 031.
 [43] E. J. Chun, Phys. Rev. D 84, 043509 (2011); K. J. Bae, E. J. Chun and S. H. Im, JCAP 1203, 013 (2012); K. J. Bae, H. Baer and E. J. Chun, JCAP 1312 (2013) 028.
 [44] K. J. Bae, H. Baer, A. Lessa and H. Serce, JCAP 1410 (2014) 082.
 [45] K. J. Bae, H. Baer and E. J. Chun, Phys. Rev. D 89 (2014) 031701.
 [46] K. J. Bae, H. Baer, V. Barger, M. R. Savoy and H. Serce, Symmetry 7 (2015) 2, 788.
 [47] H. Baer, V. Barger and D. Mickelson, Phys. Lett. B 726 (2013) 330; H. Baer, V. Barger, D. Mickelson and X. Tata, arXiv:1306.4183 [hepph].