... new physics beyond the Standard Model (SM) is needed to achieve a complete description of Nature. First of all, there is overwhelming evidence, ranging from the cosmic microwave background (CMB) to the shapes of the rotation curves of spiral galaxies, that nearly 26% of the Universe is made of yet unidentified dark matter (DM) . Moreover, the SM cannot generate the primordial inflation needed to solve the horizon and flatness problems of the Universe, as well as to explain the statistically isotropic, Gaussian and nearly scale invariant fluctuations of the CMB . The SM also lacks enough CP violation to explain why the Universe contains a larger fraction of baryonic matter than of antimatter. Aside from these three problems at the interface between particle physics and cosmology, the SM suffers from a variety of intrinsic naturalness issues. In particular, the neutrino masses are disparagingly smaller than any physics scale in the SM and, similarly, the strong CP problem states that the θ-parameter of quantum chromodynamics (QCD) is constrained from measurements of the neutron electric dipole moment to lie below an unexpectedly small value.
In this Letter we show that these problems may be intertwined in a remarkably simple way, with a solution pointing to a unique new physics scale around 1011 GeV. The SM extension we consider consists just of the KSVZ axion model [3, 4] and three right-handed (RH) heavy SM-singlet neutrinos [One may also chose alternatively the DFSZ axion model. The inflationary predictions in this model stay the same, but the window in the axion mass will move to larger values . Importantly, in this case the PQ symmetry is required to be an accidental rather than an exact symmetry in order to avoid the overclosure of the Universe due to domain walls .]. This extra matter content was recently proposed in , where it was emphasized that in addition to solving the strong CP problem, providing a good dark matter candidate (the axion), explaining the origin of the small SM neutrino masses (through an induced seesaw mechanism) and the baryon asymmetry of the Universe (via thermal leptogenesis), it could also stabilize the effective potential of the SM at high energies thanks to a threshold mechanism [7, 8]. This extension also leads to successful primordial inflation by using the modulus of the KSVZ SM singlet scalar field . Adding a cosmological constant to account for the present acceleration of the Universe, this Standard Model Axion Seesaw Higgs portal inflation (SMASH) model offers a self-contained description of particle physics from the electroweak scale to the Planck scale and of cosmology from inflation until today...
We extend the SM with a new complex singlet scalar field σ and two Weyl fermions Q and Q in the 3 and 3 representations of SU(3)c and with charges −1/3 and 1/3 under U(1)Y . With these charges, Q and Q can decay into SM quarks, which ensures that they will not become too abundant in the early Universe. We also add three RH fermions Ni. The model is endowed with a new Peccei-Quinn (PQ) global U(1) symmetry, which also plays the role of lepton number. The charges under this symmetry are: q(1/2), u(−1/2), d(−1/2), L(1/2), N(−1/2), E(−1/2), Q(−1/2) Q(−1/2), σ(1); and the rest of the SM fields (e.g. the Higgs) are uncharged. The new Yukawa couplings are: L ⊃ −[FijLiHNj + 1/2 YijσNiNj + yQσQ + yQdiσQdi + h.c.]. The two first terms realise the seesaw mechanism once σ acquires a vacuum expectation value (VEV) <σ> = vσ/√2, giving a neutrino mass matrix of the form mv = −FY -1Ft v2/(√2 vσ), with v = 246 GeV. The strong CP problem is solved as in the standard KSVZ scenario, with the role of the axion decay constant, fA, played by vσ = fA. Due to non-perturbative QCD effects, the angular part of σ = (ρ + vσ) exp(iA/fA)/√2, the axion field A, gains a potential with an absolute minimum at A = 0. At energies above the QCD scale, the axion-gluon coupling is L ⊃ −(αS/8π)(A/fA)GG, solving the strong CP problem when <A> relaxes to zero. The latest lattice computation of the axion mass gives mA=(57.2 ± 0.7)(1011GeV/fA) µeV .
The scalar sector of the model has the potential
V(H,σ) = λH ( H†H − v2/2)2 + λσ (|σ|2 − vσ2 /2)2 + 2λHσ (H†H − v2/2)(|σ|2 − vσ2 /2). (1)
In the unitary gauge, there are two scalar fields that could drive inflation: h, the neutral component of the Higgs doublet Ht=(0, h)/√2, and the modulus of the new singlet, ρ = √2|σ|.... inflation in SMASH is mostly driven by ρ, with a non-minimal coupling 2×10-3≲ξσ≲1. The upper bound on ξσ ensures that the scale of perturbative unitarity breaking is at MP (provided that also ξH ≲1), whereas the lower bound on ξσ corresponds to a tensor-to-scalar ratio r ≲ 0.07 (as constrained by the Planck satellite and the BICEP/Keck array [1, 29]). Neglecting ξH, predictive slow-roll inflation in SMASH in the Einstein frame can be described by a single canonically normalized field χ with potential
V(χ) = (λ/4)ρ(χ)4 (1 + ξσ ρ(χ)2/MP2)-2 , (2)
where λ can be either λσ or λσ =λσ −λHσ 2/λH, with the second case being possible only if λHσ < 0, corresponding to an inflationary valley in a mixed direction in the plane (ρ, h). The field χ is the solution of Ω2dχ/dρ≃ (bΩ2 + 6ξσ2 ρ2/M2P)1/2, where Ω≃ 1 + ξσ ρ2/MP2 is the Weyl transformation into the Einstein frame and b=1 for λ=λσ or b=1+|λHσ /λH|∼ 1 for λ=λσ. The value of b determines the angle in field space described by the inflationary trajectory: h2/ρ2≃ b −1. The predictions in the case λ = λσ (or b → 1) for r vs the scalar spectral index ns are shown in FIG. 1 for various values of ξσ. The running of ns is in the ballpark of 10-4–10-3, which may be probed e.g. by future observations of the 21 cm emission line of Hydrogen . These values of the primordial parameters are perfectly compatible with the latest CMB data, and the amount of inflation that is produced solves the horizon and flatness problems. Given the current bounds on r and ns, fully consistent (and predictive) inflation in SMASH occurs if 10-13≲ λ≲10-9.
For the measured central values of the Higgs and top quark masses, the Higgs quartic coupling of the SM becomes negative at h = ΛI ∼ 1011 GeV . If no new physics changes this behaviour, Higgs inflation is not viable, since it requires a positive potential at Planckian field values...Remarkably, the Higgs portal term ∝ λHσ in (1) allows absolute stability (even when the corresponding low-energy SM potential would be negative if extrapolated to large h) via the threshold-stabilisation mechanism of [7, 8, 22]. In SMASH, instabilities could also originate in the ρ direction due to quantum corrections from the RH neutrinos and KSVZ fermions. For λHσ > 0, absolute stability requires
λH, λσ > 0, for h < √2 Λh and λH, λσ > 0, for h > √2 Λh , (3)
where we define Λh2 = λHσ vσ2 /λH, λH = λH − λHσ 2/λσ and λσ = λσ −λHσ 2/λH . Instead, for λHσ < 0, the stability condition is λH, λσ > 0, for all h .
An analysis based on two-loop renormalization group (RG) equations for the SMASH couplings and one-loop matching with the SM  shows that stability can be achieved for δ ≡ λHσ 2/λσ between 10-3 and 10-1, depending on mt ... The Yukawas must satisfy the bound 6y4 + ∑Yii4 ≲ 16π2 2λσ / log(30MP/√2 λσvσ). It will prove convenient to define SMASH benchmark units:
λ10= λσ /10-10 ; δ3 = δ /0.03 ; v11= vσ /1011 GeV. (4)...
For λHσ>0, the PQ symmetry is restored nonthermally after inflation and then spontaneously broken again before reheating. On the other hand, for λHσ<0 and efficient reheating, the restoration and breaking are thermal. In the phase transition, which happens at a critical temperature Tc ≳ λσ 1/4 vσ , a network of cosmic strings is formed. Its evolution leads to a low-momentum population of axions that together with those arising from the realignment mechanism [43–45] constitute the dark matter in SMASH. Requiring that all the DM is made of axions restricts the symmetry breaking scale to the range 3×1010GeV≲vσ≲1.2×1011GeV, which translates into the mass window 50 µeV≲mA≲200 µeV, (6) updating the results of  with the latest axion mass data . The main uncertainty now arises from the string contribution [46, 47], which is expected to be diminished in the near future [48, 49]. Importantly, the SMASH axion mass window (6) will be probed in the upcoming decade by axion dark matter direct detection experiments such as CULTASK , MADMAX , and ORPHEUS , see also [23, 53] and FIG. 3 for our estimates of their future sensitivity...
The origin of the baryon asymmetry of the Universe is explained in SMASH from thermal leptogenesis . This requires massive RH neutrinos acquiring equilibrium abundances and then decaying when production rates become Boltzmann suppressed. If λHσ <0, then Treheating > Tc for stable models in the DM window (5). The RH neutrinos become massive after the PQ phase transition, and those with masses Mi < Tc retain an equilibrium abundance. The stability bound on the Yukawas Yii enforces Tc > M1, so that at least the lightest RH neutrino stays in equilibrium. Moreover, the annihilations of the RH neutrinos tend to be suppressed with respect to their decays. This allows for vanilla leptogenesis from the decays of a single RH neutrino, which demands M1 ≳ 5×108 GeV [58, 59]. However, for vσ as in (5), this is just borderline compatible with stability. Nevertheless, leptogenesis can occur with a mild resonant enhancement  for a less hierarchical RH neutrino spectrum, which relaxes the stability bound and ensures that all the RH neutrinos remain in equilibrium after the phase transitions.