Search

We study the cosmological evolution of a light QCD axion and identify the parameter space to obtain the correct relic dark matter abundance. The axion potential is flattened at the origin, corresponding to the only minimum, while it is unsuppressed at $\pi$. These potential features arise by assuming a mirror sector with the strong CP phase $\bar\theta$ shifted by $\pi$ compared to the SM sector, which allows the mirror axion potential to be tuned against the usual QCD axion potential. Before the QCD phase transition, assuming the mirror sector is decoupled and much colder than the SM thermal bath, the mirror sector potential dominates, causing the axion to initially roll to a temporary minimum at $\pi$. However, after the QCD phase transition, the potential minimum changes, and the axion relaxes from the newly created "hilltop" near $\pi$ to the CP-conserving minimum at the origin. As the axion adiabatically tracks this shift in the potential minimum through the QCD phase transition, with non-adiabatic evolution near $\pi$ and 0, it alters the usual prediction of the dark matter abundance. Consequently, this "hilltop" misalignment mechanism opens new regions of axion parameter space, with the correct relic abundance while still solving the strong CP problem, that could be explored in future experiments., Comment: 28 pages, 9 figures. Match to the published version in JHEP

We show that a rotating axion field that makes a transition from a matter-like equation of state to a kination-like equation of state around the epoch of recombination can significantly ameliorate the Hubble tension, i.e., the discrepancy between the determinations of the present-day expansion rate $H_0$ from observations of the cosmic microwave background on one hand and Type Ia supernovae on the other. We consider a specific, UV-complete model of such a rotating axion and find that it can relax the Hubble tension without exacerbating tensions in determinations of other cosmological parameters, in particular the amplitude of matter fluctuations $S_8$. We subsequently demonstrate how this rotating axion model can also generate the baryon asymmetry of our universe, by introducing a coupling of the axion field to right-handed neutrinos. This baryogenesis model predicts heavy neutral leptons that are most naturally within reach of future lepton colliders, but in finely-tuned regions of parameter space may also be accessible at the high-luminosity LHC and the beam dump experiment SHiP., Comment: 19 pages, 6 figures

We study lepto-axiogenesis in theories where the right-handed neutrino is light enough that its dynamics affect the determination of the baryon asymmetry. When compared with theories of high-scale lepto-axiogenesis where the Majorana neutrino mass may be treated as an effective dimension-five operator, we find that the predicted saxion mass is lower. Two distinct scenarios emerge. In the first, processes that generate the baryon asymmetry are in equilibrium down to the mass of the right-handed neutrino. In the second, the relevant processes never reach equilibrium; the baryon number freezes in. We comment on implications for supersymmetric spectra and discuss constraints on late decays of supersymmetric relics and from dark radiation. In contrast to high-scale lepto-axiogenesis, which predicts superpartners with masses of 10-100 TeV or more, we find this scenario is consistent with a wider range of superpartner masses, all the way down to current direct search bounds., Comment: 59 pages, 6 figures

Axion rotations can simultaneously explain the dark matter abundance and the baryon asymmetry of the Universe by kinetic misalignment and axiogenesis. We consider a scenario in which the Peccei-Quinn symmetry breaking field is as large as the Planck scale during inflation and the axion rotation is initiated by the inflaton-induced potential immediately after the end of inflation. This is a realization of the cogenesis scenario that is free of problems with domain walls and isocurvature perturbations thanks to large explicit Peccei-Quinn symmetry breaking at the Planck scale during inflation. The baryon asymmetry can be more efficiently produced by lepto-axiogenesis, in which case the axion mass is predicted to be larger than $O(0.1)$ meV. We also discuss a UV complete model in supersymmetric theories., Comment: 12 pages, 1 figure; v2: discussions added

Generating axion dark matter through the kinetic misalignment mechanism implies the generation of large asymmetries for Standard Model fermions in the early universe. Even if these asymmetries are washed out at later times, they can trigger a chiral plasma instability in the early universe. Similarly, a direct coupling of the axion with the hypercharge gauge field can trigger a tachyonic instability. These instabilities produce helical magnetic fields, which are preserved until the electroweak phase transition. At the electroweak phase transition, these become a source of baryon asymmetry, which can be much more efficient than the original axiogenesis proposal. We discuss constraints on axion dark matter production from the overproduction of the baryon asymmetry as well as a minimal, albeit fine-tuned setup, where both the correct dark matter abundance and baryon asymmetry can be achieved. For a given axion decay constant, this leads to a sharp prediction for the mass of the radial direction of the Peccei Quinn field, which is a soft mass scale in supersymmetric theories., Comment: 29 pages + references, 5 figures

We show that a minimal scenario, utilizing only the graviton as an intermediate messenger between the inflaton, the dark sector and the Standard Model (SM), is able to generate $simultaneously$ the observed relic density of dark matter (DM), the baryon asymmetry through leptogenesis, as well as a sufficiently hot thermal bath after inflation. We assume an inflaton potential of the form $V(\phi)\propto \phi^k$ about the minimum at the end of inflation. The possibility of reheating via minimal gravitational interactions has been excluded by constraints on dark radiation for excessive gravitational waves produced from inflation. We thus extend the minimal model in several ways: i) we consider non-minimal gravitational couplings--this points to the parameter range of DM masses $M_{N_1} \simeq 2-10 $ PeV, and right-handed neutrino masses $M_{N_2} \simeq (5-20) \times 10^{11}$ GeV, and $T_\text{rh} \lesssim 3 \times 10^5$ GeV (for $k \le 20$); ii) we propose an explanation for the PeV excess observed by IceCube when the DM has a direct but small Yukawa coupling to the SM; and iii) we also propose a novel scenario, where the gravitational production of DM is a two-step process, first through the production of two scalars, which then decay to fermionic DM final states. In this case, the absence of a helicity suppression enhances the production of DM and baryon asymmetry, and allows a great range for the parameters including a dark matter mass below an MeV where dark matter warmness can be observable by cosmic 21-cm lines, even when gravitational interactions are responsible for reheating. We also show that detectable primordial gravitational wave signals provide the opportunity to probe this scenario for $T_\text{rh}\lesssim 5\times 10^6$ GeV in future experiments, such as BBO, DECIGO, CE and ET., Comment: 38 pages, 11 figures, 3 tables, Gravitational wave analysis added, version accepted for publication in JHEP

Heavy QCD axions are well-motivated extensions of the QCD axion that address the quality problem while still solving the strong CP problem. Owing to the gluon coupling, critical for solving the strong CP problem, these axions can be produced in significant numbers in beam dump and collider environments for axion decay constants as large as PeV, relevant for addressing the axion quality problem. In addition, if these axions have leptonic couplings, they can give rise to long-lived decay into lepton pairs, in particular, dominantly into muons above the dimuon threshold and below the GeV scale in a broad class of axion models. Considering existing constraints, primarily from rare meson decays, we demonstrate that current and future neutrino facilities and long-lived particle searches have the potential to probe significant parts of the heavy QCD axion parameter space via dimuon final states., Comment: 19 pages + references, 6 figures

If the Peccei-Quinn field containing the QCD axion undergoes rotations in the early universe, the dimension-five operator responsible for neutrino masses can generate a lepton asymmetry that ultimately gives rise to the observed baryon asymmetry of the Universe. This lepto-axiogenesis scenario requires a flat potential for the radial direction of the Peccei-Quinn field, naturally realized in supersymmetric models. We carefully compute the efficiency of this mechanism for the Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) and Kim-Shifman-Vainshtein-Zakharov (KSVZ) axion models and place lower bounds on the masses of scalar superpartners required to reproduce the observed baryon asymmetry. For the KSVZ model, we find an efficiency for generation of the asymmetry six times larger than the previously extant computation after including scattering channels involving superpartners. In this case, the superpartner scale should be above $\sim$ 30 TeV for a domain wall number of one; the lower bound weakens for larger domain wall numbers. We find that the superpartner mass scale may also be as low as 30 TeV for the DFSZ model. In all cases, the lower bound on the superpartner masses is inversely proportional to the sum of the squares of the neutrino masses and so can strengthen as the upper bound on the neutrino mass improves. We identify the parameter space where the axion rotation can simultaneously produce axion dark matter via kinetic misalignment; in this case it is possible to put an upper bound of order PeV on the masses of scalar superpartners., Comment: 59 pages, 9 figures

Supersymmetric extensions of the Standard Model with scalar superpartners above 10 TeV are well motivated since the Higgs boson mass can be explained by quantum corrections while maintaining gauge coupling unification. If supersymmetry breaking is transmitted to gauginos via anomaly mediation, the gaugino masses are loop suppressed compared to scalar masses, and the lightest supersymmetric particle is the Higgsino or wino, which can be the dark matter. In this setup, we identify the regions of parameter space that reproduce the observed Higgs boson mass and the thermal abundance of dark matter. We analyze the effects of complex phases in the gaugino mass parameters on the electron electric dipole moment (EDM) and the dark matter scattering cross section. We find that, for scalar masses up to 10 PeV and any size of the complex phases, the model with Higgsino dark matter is within reach of planned experiments -- Advanced ACME via electron EDM and LUX-ZEPLIN via dark matter direct detection -- with complementary discovery potentials, and the model with wino dark matter is within reach of future electron EDM experiments., Comment: 9 pages, 6 figures

We demonstrate that the observed cosmological excess of matter over antimatter may originate from a heavy QCD axion that solves the strong CP problem but has a mass much larger than that given by the Standard Model QCD strong dynamics. We investigate a rotation of the heavy QCD axion in field space, which is transferred into a baryon asymmetry through weak and strong sphaleron processes. This provides a strong cosmological motivation for heavy QCD axions, which are of high experimental interest. The viable parameter space has an axion mass $m_a$ between 1~MeV and 10 GeV and a decay constant $f_a < 10^5$ GeV, which can be probed by accelerator-based direct axion searches and observations of the cosmic microwave background., Comment: 37 pages, 6 figures

We consider the generation of the baryon asymmetry in models with right-handed neutrinos produced through gravitational scattering of the inflaton during reheating. The right-handed neutrinos later decay and generate a lepton asymmetry, which is partially converted to a baryon asymmetry by Standard Model sphaleron processes. We find that a sufficient asymmetry can be generated for a wide range of right-handed neutrino masses and reheating temperatures. We also show that the same type of gravitational scattering produces Standard Model Higgs bosons, which can achieve inflationary reheating consistent with the production of a baryon asymmetry., Comment: 6 pages, 3 figures

Theoretical investigations into the evolution of the early universe are an essential part of particle physics that allow us to identify viable extensions to the Standard Model as well as motivated parameter space that can be probed by various experiments and observations. In this white paper, we review particle physics models of the early universe. First, we outline various models that explain two essential ingredients of the early universe (dark matter and baryon asymmetry) and those that seek to address current observational anomalies. We then discuss dynamics of the early universe in models of neutrino masses, axions, and several solutions to the electroweak hierarchy problem. Finally, we review solutions to naturalness problems of the Standard Model that employ cosmological dynamics., Comment: contribution to Snowmass 2021

The Standard Model of Particle Physics cannot explain the observed baryon asymmetry of the Universe. This observation is a clear sign of new physics beyond the Standard Model. There have been many recent theoretical developments to address this question. Critically, many new physics models that generate the baryon asymmetry have a wide range of repercussions for many areas of theoretical and experimental particle physics. This white paper provides an overview of such recent theoretical developments with an emphasis on experimental testability., Comment: Contribution to Snowmass 2021. Solicited white paper from TF08. v2: includes additional contributions

Complex scalar fields charged under approximate $U(1)$ symmetries appear in well-motivated extensions of the Standard Model. One example is the field that contains the QCD axion field associated with the Peccei-Quinn symmetry; others include flat directions in supersymmetric theories with baryon, lepton, or flavor charges. These fields may take on large values and rotate in field space in the early universe. The relevant approximate $U(1)$ symmetry ensures that the angular direction of the complex field is light during inflation and that the rotation is thermodynamically stable and is long-lived. These properties allow rotating complex scalar fields to naturally serve as curvatons and explain the observed perturbations of the universe. The scenario imprints non-Gaussianity in the curvature perturbations, likely at a level detectable in future large scale structure observations. The rotation can also explain the baryon asymmetry of the universe without producing excessive isocurvature perturbations., Comment: 39 pages + references, 4 figures

We show that the rotation of the QCD axion field, aided by $B-L$ violation from supersymmetric $R$-parity violating couplings, can yield the observed baryon abundance. Strong sphaleron processes transfer the angular momentum of the axion field into a quark chiral asymmetry, which $R$-parity violating couplings convert to the baryon asymmetry of the Universe. We focus on the case of dimensionless $R$-parity violating couplings with textures motivated by grand unified theories and comment on more general scenarios. The axion decay constant and mass spectrum of supersymmetric particles are constrained by Big Bang nucleosynthesis, proton decay from the $R$-parity violation, and successful thermalization of the Peccei-Quinn symmetry breaking field. Axion dark matter may be produced by the axion rotation via the kinetic misalignment mechanism for axion decay constants below $10^{10}$ GeV, or by the conventional misalignment mechanism for $10^{11-12}$ GeV. The viable parameter region can be probed by proton decay and axion searches. This scenario may also have connections with collider experiments, including searches for long-lived particles, and observations of gravitational waves., Comment: 40 pages + references, 5 figures

Rotations of an axion field in field space provide a natural origin for an era of kination domination, where the energy density is dominated by the kinetic term of the axion field, preceded by an early era of matter domination. Remarkably, no entropy is produced at the end of matter domination and hence these eras of matter and kination domination may occur even after Big Bang Nucleosynthesis. We derive constraints on these eras from both the cosmic microwave background and Big Bang Nucleosynthesis. We investigate how this cosmological scenario affects the spectrum of possible primordial gravitational waves and find that the spectrum features a triangular peak. We discuss how future observations of gravitational waves can probe the viable parameter space, including regions that produce axion dark matter by the kinetic misalignment mechanism or the baryon asymmetry by axiogenesis. For QCD axion dark matter produced by the kinetic misalignment mechanism, a modification to the inflationary gravitational wave spectrum occurs above 0.01 Hz and, for high values of the energy scale of inflation, the prospects for discovery are good. We briefly comment on implications for structure formation of the universe., Comment: 57 pages, 17 figures; v2: references and discussions added; v3: matches journal version

Split higgsinos are a compelling class of models to explain dark matter and may be on the verge of detection by multiple current experimental avenues. The idea is based on a large split in scales between the electroweak scale and decoupled scalars, with relatively light higgsinos between the two. Such models enjoy the merit of depending on very few parameters while still explaining gauge coupling unification, dark matter, and most of the hierarchy between the Planck and electroweak scales, and they remain undetected by past experiments. We analyze split higgsinos in view of current and next generation experiments. We discuss the direct and indirect detection prospects and further demonstrate promising discovery potentials in the upcoming electron electric dipole moment experiments. The parameter space of this model is analyzed in terms of experiments expected to run in the coming years and where we should be looking for the next potential discoveries., Comment: 10 pages, 7 figures

An axion rotating in field space can produce dark photons in the early universe via tachyonic instability. This explosive particle production creates a background of stochastic gravitational waves that may be visible at pulsar timing arrays or other gravitational wave detectors. This scenario provides a novel history for dark photon dark matter. The dark photons may be warm at a level detectable in future 21-cm line surveys. For a consistent cosmology, the radial direction of the complex field containing the axion must be thermalized. We explore a concrete thermalization mechanism in detail and also demonstrate how this setup can be responsible for the generation of the observed baryon asymmetry., Comment: 47 pages, 2 figures

Reheating is a process where the energy density of a dominant component of the universe other than radiation, such as a matter component, is transferred into radiation. It is usually assumed that the temperature of the universe decreases due to cosmic expansion even during the reheating process, in which case the maximal temperature of the universe is much higher than the reheat temperature. We point out that the temperature of the universe during reheating may in fact increase in well-motivated scenarios. We derive the necessary conditions for the temperature to increase during reheating and discuss concrete examples involving a scalar field. We comment on implications for particle physics and cosmology due to an increasing temperature during reheating., Comment: 24 pages, 7 figures

We propose a baryogenenesis mechanism that uses a rotating condensate of a Peccei-Quinn (PQ) symmetry breaking field and the dimension-five operator that gives Majorana neutrino masses. The rotation induces charge asymmetries for the Higgs boson and for lepton chirality through sphaleron processes and Yukawa interactions. The dimension-five interaction transfers these asymmetries to the lepton asymmetry, which in turn is transferred into the baryon asymmetry through the electroweak sphaleron process. QCD axion dark matter can be simultaneously produced by dynamics of the same PQ field via kinetic misalignment or parametric resonance, favoring an axion decay constant $f_a \lesssim 10^{10}$ GeV, or by conventional misalignment and contributions from strings and domain walls with $f_a \sim 10^{11}$ GeV. The size of the baryon asymmetry is tied to the mass of the PQ field. In simple supersymmetric theories, it is independent of UV parameters and predicts the supersymmtry breaking mass scale to be $\mathcal{O}(10-10^4)$ TeV, depending on the masses of the neutrinos and whether the condensate is thermalized during a radiation or matter dominated era. We also construct a theory where TeV scale supersymmetry is possible. Parametric resonance may give warm axions, and the radial component of the PQ field may give signals in rare kaon decays from mixing with the Higgs and in dark radiation., Comment: 71 pages, 6 Figures

Adding an axion-like particle (ALP) to the Standard Model, with a field velocity in the early universe, simultaneously explains the observed baryon and dark matter densities. This requires one or more couplings between the ALP and photons, nucleons, and/or electrons that are predicted as functions of the ALP mass. These predictions arise because the ratio of dark matter to baryon densities is independent of the ALP field velocity, allowing a correlation between the ALP mass, $m_a$, and decay constant, $f_a$. The predicted couplings are orders of magnitude larger than those for the QCD axion and for dark matter from the conventional ALP misalignment mechanism. As a result, this scheme, ALP cogenesis, is within reach of future experimental ALP searches from the lab and stellar objects, and for dark matter., Comment: 24 pages, 3 figures

Axion cold dark matter from standard misalignment typically requires a decay constant $f_a~\gtrsim~10^{11}$ GeV. Kinetic misalignment and parametric resonance easily allow lower values of $f_a$ when the radial Peccei-Quinn (PQ) symmetry breaking field takes large initial values. Here, we consider the effects of inflation on kinetic misalignment and parametric resonance. We assume that the initial PQ field value is determined by quantum fluctuations, and is set by the Hubble parameter during inflation, $H_I$, and the PQ field mass. PQ field oscillations begin before or after the completion of reheating after inflation at a temperature $T_R$. We determine the range of $f_a$ and the inflationary parameters $(H_I, T_R)$ consistent with axion dark matter for a quartic potential for the PQ field. We find that $4\times 10^8$ GeV $< f_a < 10^{11}$ GeV can consistently produce axion dark matter. A significant portion of the allowed parameter space predicts rare kaon decays, $K_L \rightarrow (\pi^0 + \rm{missing \; energy})$, and/or suppression of structure formation on small scales., Comment: 40 pages, 6 figures

We consider a cosmological scenario where a relativistic particle and a stable massive particle are simultaneously produced from the decay of a late-decaying particle after Big-Bang Nucleosynthesis but before matter-radiation equality. The relativistic and massive particles behave as dark radiation and warm dark matter, respectively. Due to a common origin, the warmness and abundances are closely related. We refer to the models that lead to such a scenario as Common Origin of Warm and Relativistic Decay Products (COWaRD). We show that COWaRD predicts a correlation between the amount of dark radiation and suppression of the large scale structure, which can be tested in future precision cosmology observations. We demonstrate that COWaRD is realized, as an example, in a class of supersymmetric axion models and that future observations by the next generation Cosmic Microwave Background, Large Scale Structure, and 21-cm surveys can reveal the structure of the theory., Comment: Version accepted for publication. 22 pages + 2 appendices and references

In the conventional misalignment mechanism, the axion field has a constant initial field value in the early universe and later begins to oscillate. We present an alternative scenario where the axion field has a nonzero initial velocity, allowing an axion decay constant much below the conventional prediction from axion dark matter. This axion velocity can be generated from explicit breaking of the axion shift symmetry in the early universe, which may occur as this symmetry is approximate., Comment: 7+4 pages, 2+2 figures; v2: Supplemental Material and references added, matches journal version

Axions may be produced thermally inside the cores of neutron stars (NSs), escape the stars due to their feeble interactions with matter, and subsequently convert into X-rays in the magnetic fields surrounding the stars. We show that a recently-discovered excess of hard X-ray emission in the 2 - 8 keV energy range from the nearby Magnificent Seven isolated NSs could be explained by this emission mechanism. These NSs are unique in that they had previously been expected to only produce observable flux in the UV and soft X-ray bands from thermal surface emission at temperatures ~100 eV. No conventional astrophysical explanation of the Magnificent Seven hard X-ray excess exists at present. We show that the hard X-ray excess may be consistently explained by an axion-like particle with mass $m_a \lesssim 2 \times 10^{-5}$ eV and $g_{a\gamma\gamma} \times g_{ann} \in (2 \times 10^{-21}, 10^{-18})$ GeV$^{-1}$ at 95\% confidence, accounting for both statistical and theoretical uncertainties, where $g_{a\gamma\gamma}$ ($g_{ann}$) is the axion-photon (axion-neutron) coupling constant., Comment: 7 + 10 pages, 4 + 8 figures, version published in PRL

We propose a mechanism called axiogenesis where the cosmological excess of baryons over antibaryons is generated from the rotation of the QCD axion. The Peccei-Quinn (PQ) symmetry may be explicitly broken in the early universe, inducing the rotation of a PQ charged scalar field. The rotation corresponds to the asymmetry of the PQ charge, which is converted into the baryon asymmetry via QCD and electroweak sphaleron transitions. In the concrete model we explore, interesting phenomenology arises due to the prediction of a small decay constant and the connections with new physics at the LHC and future colliders and with axion dark matter., Comment: 6+10 pages, 1+2 figures; v2: a simplified model and references added, matches journal version

Several theoretical motivations point to ultralight QCD axions with large decay constants $f_a \simeq \mathcal{O}(10^{16}-10^{17})$ GeV, to which experimental proposals are dedicated. This regime is known to face the problem of overproduction of axion dark matter from the misalignment mechanism unless the misalignment angle $\theta_{\rm mis}$ is as small as $\mathcal{O}(10^{-3}-10^{-4})$, which is generally considered a fine-tuning problem. We investigate a dynamical explanation for a small $\theta_{\rm mis}$. The axion mass arises from strong dynamics and may be sufficiently enhanced by early dynamics so as to overcome Hubble friction and drive the field value to the bottom of the potential long before the QCD phase transition. Together with an approximate CP symmetry in the theory, this minimum is very closely related to today's value and thus $\theta_{\rm mis}$ can automatically be well under unity. Owing to such efficient relaxation, the isocurvature perturbations are essentially damped. As an existence proof, using supersymmetric theories we illustrate that the Higgs coupling with the inflaton energy can successfully achieve this axion damping in a consistent inflationary cosmology., Comment: 16 pages, 2 figures; v2: minor changes and references added, matches journal version

The QCD axion serves as a well-motivated dark matter candidate and the misalignment mechanism is known to reproduce the observed abundance with a decay constant $f_a \simeq \mathcal{O}(10^{12})$ GeV for a misalignment angle $\theta_{\rm mis} \simeq \mathcal{O}(1)$. While $f_a \ll 10^{12}$ GeV is of great experimental interest, the misalignment mechanism requires the axion to be very close to the hilltop, i.e. $\theta_{\rm mis} \simeq \pi$. This particular choice of $\theta_{\rm mis}$ has been understood as fine-tuning the initial condition. We offer a dynamical explanation for $\theta_{\rm mis} \simeq \pi$ in a class of models. The axion dynamically relaxes to the minimum of the potential by virtue of an enhanced mass in the early universe. This minimum is subsequently converted to a hilltop because the CP phase of the theory shifts by $\pi$ when one contribution becomes subdominant to another with an opposite sign. We demonstrate explicit and viable examples in supersymmetric models where the higher dimensional Higgs coupling with the inflaton naturally achieves both criteria. Associated phenomenology includes a strikingly sharp prediction of $3 \times 10^9$ GeV $\lesssim f_a \lesssim 10^{10}$ GeV and the absence of isocurvature perturbation., Comment: 25 pages, 3 figures; v2: discussions and references added, matches journal version

Despite growing interest and extensive effort to search for ultralight dark matter in the form of a hypothetical dark photon, how it fits into a consistent cosmology is unclear. Several dark photon dark matter production mechanisms proposed previously are known to have limitations, at least in certain mass regimes of experimental interest. In this letter, we explore a novel mechanism, where a coherently oscillating axion-like field can efficiently transfer its energy density to a dark photon field via a tachyonic instability. The residual axion relic is subsequently depleted via couplings to the visible sector, leaving only the dark photon as dark matter. We ensure that the cosmologies of both the axion and dark photon are consistent with existing constraints. We find that the mechanism works for a broad range of dark photon masses, including those of interest for ongoing experiments and proposed detection techniques., Comment: 8 pages, 2 figues; v2: minor changes, matches journal version

We examine the theoretical motivations for long-lived particle (LLP) signals at the LHC in a comprehensive survey of Standard Model (SM) extensions. LLPs are a common prediction of a wide range of theories that address unsolved fundamental mysteries such as naturalness, dark matter, baryogenesis and neutrino masses, and represent a natural and generic possibility for physics beyond the SM (BSM). In most cases the LLP lifetime can be treated as a free parameter from the $\mu$m scale up to the Big Bang Nucleosynthesis limit of $\sim 10^7$m. Neutral LLPs with lifetimes above $\sim$ 100m are particularly difficult to probe, as the sensitivity of the LHC main detectors is limited by challenging backgrounds, triggers, and small acceptances. MATHUSLA is a proposal for a minimally instrumented, large-volume surface detector near ATLAS or CMS. It would search for neutral LLPs produced in HL-LHC collisions by reconstructing displaced vertices (DVs) in a low-background environment, extending the sensitivity of the main detectors by orders of magnitude in the long-lifetime regime. In this white paper we study the LLP physics opportunities afforded by a MATHUSLA-like detector at the HL-LHC. We develop a model-independent approach to describe the sensitivity of MATHUSLA to BSM LLP signals, and compare it to DV and missing energy searches at ATLAS or CMS. We then explore the BSM motivations for LLPs in considerable detail, presenting a large number of new sensitivity studies. While our discussion is especially oriented towards the long-lifetime regime at MATHUSLA, this survey underlines the importance of a varied LLP search program at the LHC in general. By synthesizing these results into a general discussion of the top-down and bottom-up motivations for LLP searches, it is our aim to demonstrate the exceptional strength and breadth of the physics case for the construction of the MATHUSLA detector., Comment: 213 pages, 73 figures. Extended Section 2 to add more detailed discussion of LLP reconstruction and analysis, and background rejection. Updated comparison of MATHUSLA RH neutrino sensitivity to other experiments. Updated analysis of long-lived ALPs produced in weak-scale processes and decaying to jets. Various clarifications, fixed typos, and added references. Results and conclusions unchanged

The QCD axion is a good dark matter candidate. The observed dark matter abundance can arise from misalignment or defect mechanisms, which generically require an axion decay constant $f_a \sim \mathcal{O}(10^{11})$ GeV (or higher). We introduce a new cosmological origin for axion dark matter, parametric resonance from oscillations of the Peccei-Quinn symmetry breaking field, that requires $f_a \sim (10^8 -10^{11})$ GeV. The axions may be warm enough to give deviations from cold dark matter in Large Scale Structure., Comment: 6 pages, 2 figures; v2: discussions and references added

In supersymmetric theories, the gravitino is abundantly produced in the early Universe from thermal scattering, resulting in a strong upper bound on the reheat temperature after inflation. We point out that the gravitino problem may be absent or very mild due to the early dynamics of a supersymmetry breaking field, i.e. a sgoldstino. In models of low scale mediation, the field value of the sgoldstino determines the mediation scale and is in general different in the early Universe from the present one. A large initial field value since the era of the inflationary reheating suppresses the gravitino production significantly. We investigate in detail the cosmological evolution of the sgoldstino and show that the reheat temperature may be much higher than the conventional upper bound, restoring the compatibility with thermal leptogenesis., Comment: 23 pages, 3 figures; v2: discussions added and one figure updated, matches version published in JHEP

In all supersymmetric theories, gravitinos, with mass suppressed by the Planck scale, are an obvious candidate for dark matter; but if gravitinos ever reached thermal equilibrium, such dark matter is apparently either too abundant or too hot, and is excluded. However, in theories with an axion, a saxion condensate is generated during an early era of cosmological history and its late decay dilutes dark matter. We show that such dilution allows previously thermalized gravitinos to account for the observed dark matter over very wide ranges of gravitino mass, keV < $m_{3/2}$ < TeV, axion decay constant, $10^9$ GeV < $f_a$ < $10^{16}$ GeV, and saxion mass, 10 MeV < $m_s$ < 100 TeV. Constraints on this parameter space are studied from BBN, supersymmetry breaking, gravitino and axino production from freeze-in and saxion decay, and from axion production from both misalignment and parametric resonance mechanisms. Large allowed regions of $(m_{3/2}, f_a, m_s)$ remain, but differ for DFSZ and KSVZ theories. Superpartner production at colliders may lead to events with displaced vertices and kinks, and may contain saxions decaying to $(WW,ZZ,hh), gg, \gamma \gamma$ or a pair of Standard Model fermions. Freeze-in may lead to a sub-dominant warm component of gravitino dark matter, and saxion decay to axions may lead to dark radiation., Comment: 30 pages, 4 figures

A new scheme for lightest supersymmetric particle (LSP) dark matter is introduced and studied in theories of TeV supersymmetry with a QCD axion, $a$, and a high reheat temperature after inflation, $T_R$. A large overproduction of axinos ($\tilde{a}$) and gravitinos ($\tilde{G}$) from scattering at $T_R$, and from freeze-in at the TeV scale, is diluted by the late decay of a saxion condensate that arises from inflation. The two lightest superpartners are $\tilde{a}$, with mass of order the TeV scale, and $\tilde{G}$ with mass $m_{3/2}$ anywhere between the keV and TeV scales, depending on the mediation scale of supersymmetry breaking. Dark matter contains both warm and cold components: for $\tilde{G}$ LSP the warm component arises from $\tilde{a} \rightarrow \tilde{G}a$, while for $\tilde{a}$ LSP the warm component arises from $\tilde{G} \rightarrow \tilde{a}a$. The free-streaming scale for the warm component is predicted to be of order 1 Mpc (and independent of $m_{3/2}$ in the case of $\tilde{G}$ LSP). $T_R$ can be as high as $10^{16}$ GeV, for any value of $m_{3/2}$, solving the gravitino problem. The PQ symmetry breaking scale $V_{PQ}$ depends on $T_R$ and $m_{3/2}$ and can be anywhere in the range $(10^{10} - 10^{16})$ GeV. Detailed predictions are made for the lifetime of the neutralino LOSP decaying to $\tilde{a}+ h/Z$ and $\tilde{G}+h/Z/\gamma$, which is in the range of $(10^{-1}-10^6)$m over much of parameter space. For an axion misalignment angle of order unity, the axion contribution to dark matter is sub-dominant, except when $V_{PQ}$ approaches $10^{16}$ GeV., Comment: 43 pages, 16 figures

We present a simple and natural dark sector model in which dark matter particles arise as composite states of hidden strong dynamics and their stability is ensured by accidental symmetries. The model has only a few free parameters. In particular, the gauge symmetry of the model forbids the masses of dark quarks, and the confinement scale of the dynamics provides the unique mass scale of the model. The gauge group contains an Abelian symmetry $U(1)_D$, which couples the dark and standard model sectors through kinetic mixing. This model, despite its simple structure, has rich and distinctive phenomenology. In the case where the dark pion becomes massive due to $U(1)_D$ quantum corrections, direct and indirect detection experiments can probe thermal relic dark matter which is generically a mixture of the dark pion and the dark baryon, and the Large Hadron Collider can discover the $U(1)_D$ gauge boson. Alternatively, if the dark pion stays light due to a specific $U(1)_D$ charge assignment of the dark quarks, then the dark pion constitutes dark radiation. The signal of this radiation is highly correlated with that of dark baryons in dark matter direct detection., Comment: 6 pages, 5 figures; references added

We introduce a class of unified supersymmetric axion theories with unified and PQ symmetries broken by the same set of fields at a scale $\sim 2\times 10^{16}$ GeV. A typical domain wall number of order 30 leads to an axion decay constant $f_a$ of order $10^{15}$ GeV. Inflation generates a large saxion condensate giving a reheat temperature $T_R$ below the QCD scale for supersymmetry breaking of order $1-10$ TeV. Axion field oscillations commence in the saxion matter-dominated era near the QCD scale, and recent lattice computations of the temperature dependence of the axion mass in this era allow a controlled calculation of the axion dark matter abundance. A successful prediction of this abundance results for an initial axion misalignment angle of order unity, $\theta_i \sim 1$. A highly correlated set of predictions is discussed for $f_a$, $T_R$, the supersymmetric Higgs mass parameter $\mu$, the amount of dark radiation $\Delta N_{eff}$, the proton decay rate $\Gamma(p \rightarrow e^+ \pi^0)$, isocurvature density perturbations and the $B$-mode of the cosmic microwave background. The last two are particularly interesting when the energy scale of inflation is also of order $10^{16}$ GeV., Comment: 29 pages, 7 figures; v2, version published in Physical Review D

Dark matter, $X$, may be generated by new physics at the TeV scale during an early matter-dominated (MD) era that ends at temperature $T_R \ll {\rm TeV}$. Compared to the conventional radiation-dominated (RD) results, yields from both Freeze-Out and Freeze-In processes are greatly suppressed by dilution from entropy production, making Freeze-Out less plausible while allowing successful Freeze-In with a much larger coupling strength. Freeze-In is typically dominated by the decay of a particle $B$ of the thermal bath, $B \rightarrow X$. For a large fraction of the relevant cosmological parameter space, the decay rate required to produce the observed dark matter abundance leads to displaced signals at LHC and future colliders, for any $m_X$ in the range ${\rm keV} < m_X < m_B$ and for values of $m_B$ accessible to these colliders. This result applies whether the early MD era arises after conventional inflation, when $T_R$ is the usual reheat temperature, or is a generic MD era with an alternative origin. In the former case, if $m_X$ is sufficiently large to be measured from kinematics, the reheat temperature $T_R$ can be extracted. Our result is independent of the particular particle physics implementation of $B \rightarrow X$, and can occur via any operator of dimension less than 8 (4) for a post-inflation (general MD) cosmology. An interesting example is provided by DFS axion theories with TeV-scale supersymmetry and axino dark matter of mass GeV to TeV, which is typically overproduced in a conventional RD cosmology. If $B$ is the higgsino, $\tilde h$, Higgs, W and Z particles appear at the displaced decays, $\tilde h \rightarrow h \tilde a, Z \tilde a$ and $\tilde h^\pm \rightarrow W^\pm \tilde a$. The scale of axion physics, $f$, is predicted to be in the range $(3\times10^8 - 10^{12})$ GeV and, over much of this range, can be extracted from the decay length., Comment: 26 pages + appendices, 13 figures; v2, JCAP version