All-in-one Relaxion, a unified solution to five BSM puzzles

We present a unified relaxion solution to the five major outstanding issues in particle physics: the hierarchy problem, dark matter, matter-antimatter asymmetry, neutrino masses and the strong CP problem. The only additional field content in our construction with respect to standard relaxion models is an up-type vector-like fermion pair and three right-handed neutrinos charged under the relaxion shift symmetry. The observed dark matter abundance is generated automatically by oscillations of the relaxion field that begin once it is misaligned from its original stopping point after reheating. The matter-antimatter asymmetry arises from spontaneous baryogenesis induced by the CPT violation due to the rolling of the relaxion after reheating. The CPT violation is communicated to the baryons and leptons via an operator, $\partial_\mu \phi J^\mu$, where $J^\mu$ consists of right-handed neutrino currents arising naturally from a simple neutrino mass model. Finally, the strong CP problem is solved via the Nelson-Barr mechanism, i.e. by imposing CP as a symmetry of the Lagrangian that is broken only spontaneously by the relaxion. The CP breaking is such that although an ${\cal O}(1)$ strong CKM phase is generated, the induced strong CP phase is much smaller, i.e., within experimental bounds.


I. INTRODUCTION
In recent times particle physics research has been driven to a large extent by the expectation of physics beyond the Standard Model (BSM) at the TeV scale. While there are many theoretical and observational reasons to extend the Standard Model (SM) -such as the hierarchy problem, dark matter, matter-antimatter asymmetry, neutrino masses and the strong CP problem -only the first of these issues is believed to necessarily require TeV scale new physics.
In fact if the hierarchy problem is ignored and new physics scales much beyond the TeV scale are allowed, the other issues can be solved by very minimal extensions of the SM. For instance, in the so-called 'SMASH' model of Ref. [1,2], a simultaneous solution to neutrino masses, baryogenesis, the strong CP problem, dark matter and a model for inflation was achieved with the addition of only three right-handed SM-singlet neutrinos, a new vectorlike color triplet fermion and a complex SM singlet scalar. In another set of papers (the flaxion/axiflavon constructions of Ref. [3,4]), the Peccei-Quinn solution [5] to the strong CP and dark matter puzzles was unified with the Froggatt-Nielsen solution [6] to the SM flavour puzzle by identifying the Froggatt-Nielsen and Peccei-Quinn symmetries. Furthermore in Ref. [3], it was also shown that the radial component of the flavon field in such models can act as the inflaton and the baryon asymmetry can be explained by leptogenesis.
It is arguably much more challenging to find an explanation (apart from tuning or anthropics) for a light Higgs mass with a high new physics scale. While conventional wisdom says this is impossible, the recently-proposed cosmological relaxation (or relaxion) models [7] aim to find just such an explanation. In these models the rolling of the so-called relaxion field during inflation leads to a scanning of the the Higgs mass squared from positive to negative values. Once the Higgs mass squared becomes negative it triggers a backreaction potential that stops the scanning soon after, at a value much smaller than the cut-off, thus dynamically solving the hierarchy problem.
Building on previous work [8][9][10], we show in this paper that the relaxion construction already has ingredients that can provide solutions to multiple BSM puzzles in addition to the hierarchy problem. Our starting point is Ref. [10], in which it was shown that Spontaneous Baryogenesis [11,12] can be successfully achieved in scenarios where the relaxion rolls after reheating (due to the vanishing of the backreaction at high temperatures). The spontaneous CPT violation due to the relaxion rolling biases thermal equilibrium such that sphalerons processes lead to a net generation of B + L number. The CPT breaking is communicated to the baryons/leptons via an operator, O SB = ∂ µ φJ µ , where J µ is a current containing the B + L current. Once the temperature drops sufficiently the backreaction potential reappears and the relaxion, being misaligned from the original minima, starts oscillating [13]. Remarkably, in the finite parameter space where baryon asymmetry can be successfully explained, relaxion oscillations can also explain the observed dark matter abundance [10].
In this paper we improve the above spontaneous relaxion baryogenesis (SRB) set-up of Ref. [10] by clarifying the origin of its two non-renormalisable parts, namely the operator O SB and the rolling potential. The first improvement is a much needed explanation for the operator O SB . In this work, the current J µ consists of a combination of right-handed neutrino currents that arises naturally from a neutrino mass model; the operator O SB would thus be practically unconstrained, leading to an opening up of the parameter space with respect to Ref. [10].
Finally, we provide a renormalisable sector that can ra-diatively generate the relaxion rolling potential. Borrowing from the Nelson-Barr relaxion construction of Ref. [8], we present a model that generates the rolling potential while also solving the strong CP problem via the Nelson-Barr mechanism [14,15]. Following the basic philosophy of Nelson-Barr models, CP is assumed to be an exact symmetry in the UV and is only broken spontaneously by the relaxion in a way that does not generate a strong CP phase greater than allowed limits. On the other hand, the phase of the relaxion at its stopping point gets mapped to the CKM phase. Thus, with a well-motivated completion of the minimal SRB scenario of Ref. [10] we achieve a unified solution to five BSM puzzles, namely the hierarchy problem, dark matter, matter-antimatter asymmetry, neutrino masses and the strong CP problem. We describe our set-up in detail in the following sections.

II. REVIEW AND BASIC SET-UP
In relaxion models, the Higgs mass squared parameter, µ 2 , is promoted to a dynamical quantity µ 2 (φ), which varies due to its couplings to the new relaxion field, φ, with, Here, H is the SM Higgs doublet, λ H is its quartic coupling, and M is the UV cut-off of the Higgs effective theory [33]. The coefficient r roll is model-dependent and κ < ∼ 1. During inflation the relaxion field slow-rolls due to the last term in Eq. 1 such that the value µ 2 parameter is slowly scanned starting from positive values. Starting from a relaxion field value φ < φ c = −|F cos −1 κ|, the relaxion field slow-rolls down the slope, increasing the value of φ and decreasing the value of |µ 2 |. After crossing the point φ = φ c , µ 2 becomes negative prompting electroweak symmetry breaking. This in turn activates the backreaction potential which induces periodic 'wiggles' on top of the linear envelope, were Λ 4 c = m n v(φ) 4−n is an increasing function of the Higgs vacuum expectation value (VEV). These wiggles cause the relaxion field to come to a halt at φ = φ 0 satisfying, where we have taken sin Hence, if F/f k 1, a large hierarchy can be achieved between the Higgs VEV, v, and the cut-off M .  Figure 1: Schematic representation of our set-up with the vertical line representing the clockwork system. The three right-handed neutrinos at the 0th site, ni, couple to the SM in the usual way, generating neutrino masses; they also provide the current in the all-important operator for spontaneous baryogenesis, OSB = ∂µφJ µ . At the kth site we introduce a new strong sector which couples to SM via its fermionic matter content, (L, L c , N, N c ). This sector generates the backreaction wiggles and relaxion oscillations inside these wiggles generate the observed dark matter abundance. Finally, at site N there is a Nelson-Barr sector that radiatively generates the rolling potential while also providing a solution to the strong CP problem. This sector couples to the SM up sector via a new vector-like quark pair (ψ, ψ c ).
As discussed in Ref. [7], the cut-off, M , cannot be raised to an arbitrarily high value because of the following cosmological requirements: (1) the relaxion dynamics should be classical, i.e., quantum fluctuations should be sub-dominant, which implies that H I < V /H 2 I , H I being the Hubble scale during inflation and (2) the relaxion vacuum energy should not drive inflation and should thus be negligible compared to the total vacuum energy, i.e. M 4 < H 2 M 2 P l . Together these conditions give the following upper bound on the cut-off (where we have used Eq. 4),

M < ∼
Let us now discuss what happens after inflation. First note that for the backreaction sector we will adopt the non-QCD model of Ref. [7] where φ is the axion of a new strong sector. This sector couples to SM via vector-like electroweak doublets and singlets (L, L c ) and (N, N c ) that are charged under the new strong group and are given Yukawa couplings to the Higgs. If the reheating temperature is greater than the critical temperature of the chiral phase transition of the new sector, i.e. if T r > T c ∼ √ 4πf π , the wiggles disappear and the relaxion starts rolling again. Here f π is the analog of the 'pion decay constant' of the strong sector. When the universe cools below the temperature T c again, the backreaction potential reappears and the rolling eventually stops provided, This condition is obtained by demanding that the relaxion does not pick up enough kinetic energy to overshoot the barriers once the backreaction potential reappears [13,16,17]. We note for future reference that the backreaction scale Λ c depends not just on f π but also on other technically-natural couplings [7], such that Λ c and f π ∼ T c are two independent scales. Also note that in this model we must have Λ 4 c < ∼ 16π 2 v 4 so that Higgsindependent contributions are sub-dominant compared to Higgs-dependant ones.
It is this second phase of rolling that can lead to a generation of both the observed dark matter abundance as well as the baryon asymmetry. The explanation for dark matter requires no additional ingredient. This is due to the fact that during the second phase of rolling, the relaxion gets misaligned from its original stopping point by an angle [13], As shown in Ref. [13], this sets off relaxion oscillations that can give rise to the observed dark matter relic abundance, Note that the correct relic density can always be reproduced by choosing an appropriate value of tan φ0 f . While there is some room for this in the relaxion mechanism, as the relaxion is spread across multiple vacua at the end of its rolling, the probability distribution of the relaxion field peaks for O(1) values of tan φ0 f [18]. Thus the extent to which tan φ0 f deviates from unity can be interpreted as a measure of the tuning required to get the correct relic abundance.
It was shown in Ref. [10] that with just one additional requirement, this second phase of rolling can also explain the observed baryon asymmetry. One additionally requires that some fermions with B + L charge must be charged under the relaxion shift symmetry. This leads to the presence of the operator, ∂ µ φJ µ , where J µ contains the B + L current. This operator can generate a chemical potential for B + L violation once the second phase of relaxion rolling results in a CPT breaking expectation value for ∂ µ φ, giving rise to a baryon asymmetry due to the presence of (B + L)-violating sphaleron transitions.
As we will show later, to generate a baryon asymmetry of the observed size one requires a hierarchy f f k . This and the fact that the relaxion, in any case, requires a large hierarchy between f k and its field excursion during rolling, f k F , are problematic as pointed out in Ref. [19]. The solution to generating the latter hierarchy is embedding the relaxion construction in a so-called clockwork model [20][21][22]; this can easily be extended to also generate the former hierarchy, giving the overall hierarchy f f k F . In clockwork models there is a system of interacting complex scalars, Φ i , all of which get a VEV such that Φ i = f √ 2 e iπi/f . There is an approximate abelian symmetry, U(1) i , at each site which is spontaneously broken to give rise to a corresponding pseudo-goldstone mode π i . Explicit breaking effects give the angular fields, π i , a mass matrix such that the lightest state is a massless (Goldstone) mode given by, Given that the mixing angle π k |φ ∼ 3 −k , any Lagrangian term where the angular field π k couples with a decay constant f translates to an interaction of φ with an exponentially-enhanced effective decay constant, 3 k f , in the mass basis. Coming back to our set-up, the hierarchy f f k F can be obtained by having the operator, O SB , at the 0-th site, the backreaction sector at the intermediate k-th site and the rolling potential at the N th site, such that f k = 3 k f and F = 3 N f . This is schematically shown in Fig. 1. The rolling and backreaction potentials eventually lift the flat direction in Eq. 9, giving a mass to the relaxion.
While the SRB set-up reviewed above is an elegant solution to multiple BSM puzzles with only a minimal modification of standard relaxion models, it is incomplete in two respects: the operator, O SB , and the rolling potential are non-renormalisable and introduced in a somewhat ad hoc way. While O SB arises naturally if baryons and/or lepton are charged under the abelian symmetry of which the relaxion is a goldstone boson, charging the SM fermions seems to have no purpose other than generating O SB . Furthermore, the charge assignments have to be carefully chosen such that they are anomaly free with respect to QCD (a necessary condition to avoid generating a strong CP phase) and preferably also with respect to electromagnetism (to avoid the generation of a φγγ coupling that rules out most of the parameter space in this set-up [10]). Here we complete the SRB set-up as follows: • Instead of introducing the operator, O SB , by hand, we propose a simple neutrino mass model at site 0, which generates this operator with a current containing only three new right-handed (RH) neutrino fields. The operator arises because the RH fields are charged under the relaxion shift symmetry, which in turn is identified with a Froggatt-Nielsen symmetry. This also explains the observed smallness of neutrino masses. As only the SM singlet RH neutrinos are charged under it, the relaxion shift symmetry is automatically anomaly-free with respect to both QCD and electromagnetism.
• We show that the rolling potential can be generated by the addition of a single up-type vector-like pair at the site N . In addition, this modification can solve the strong CP problem via the Nelson-Barr mechanism [14,15]. In Nelson-Barr models, CP is a good symmetry in the UV and is spontaneously broken at an intermediate scale to generate an O(1) CKM phase but a much smaller strong CP phase (within allowed constraints). We borrow the Nelson-Barr relaxion sector from Ref. [8], where the relaxion phase upon stopping results in an O(1) CKM phase.
From the point of view of the Nelson-Barr relaxion model, baryogenesis is not straightforward to achieve due to the lack of any explicit CP violation. Spontaneous baryogenesis is thus an attractive potential feature for the Nelson-Barr relaxion, as it does not require explicit CP violation.

III.
NEUTRINO MASSES AND SPONTANEOUS BARYOGENESIS

Getting the operator OSB
At site 0 we introduce a sector that simultaneously generates small neutrino masses and an operator suitable for spontaneous baryogenesis. We introduce three RH neutrinos, n i , that are charged under the abelian symmetry at the 0-th site of the clockwork chain, U(1) 0 . We fixed the charge of Φ 0 to be −1 under this symmetry and take all SM fields to be neutral. The Lagrangian for the couplings of these right-handed neutrinos is given by, where l i are the SM lepton doublets, q nj are the abelian charges for the sterile neutrinos and M ij n are the associated Majorana mass matrix elements. Given that we will eventually use the Nelson-Barr solution to the strong CP problem, we impose that CP is an exact symmetry of the Lagrangian so that all the couplings above are real. Following spontaneous symmetry breaking, the first term becomes, The effective Yukawa couplings, Y ij n = y ij n ( F N ) qn j can thus be exponentially small, which, as we will soon see, provides a natural explanation for neutrino masses with a smaller-than-usual seesaw scale. This is just an implementation of the Froggatt-Nielsen mechanism [6] to naturally obtain small Yukawa couplings. While the first term in Eq. 10 is still non-renormalisable, we note that such operators can be UV completed in standard ways involving chains of vector-like fermions (see for eg. Ref. [9]).
The factor e iπ0qn j /f in front of the Yukawa couplings can be rotated away by the field redefinition n j → n j e −iπ0qn j /f , which, through the redefinition of the kinetic terms for the RH neutrinos, yields the desired op- where we ignore an O(1) factor corresponding to π 0 |φ .

Getting Baryon Asymmetry from OSB
Now we will demonstrate how the presence of O SB can lead to spontaneous baryogenesis, an idea developed in Refs. [11,12]. The essential feature of this mechanism is the presence of a rolling field that breaks CPT. It was shown in Ref. [10] that this role can be played by the rolling of the relaxion after reheating.
During the second phase of the rolling of the relaxion, the operator O SB causes equal and opposite shifts in the energy of particles versus antiparticles, implying, where µ i (μ i ) is the chemical potential for (anti-) particles of the ith specie; q i is its charges under U(1) 0 (which is non-zero only for the RH neutrinos); Q i is the electromagnetic charge; T 3i is the charge corresponding to the diagonal generator of SU(2) L ; and the chemical potentials µ Q,T3,B−L have been introduced to enforce conservation of Q, T 3 and B − L. In the presence of (B + L)violating sphaleron processes, we find, for fermions and bosons respectively, where f (T, µ) is the Fermi-Dirac (Bose-Einstein) distribution for fermions (bosons). We have taken µ T and the factor g i denotes the number of degrees of freedom for each species. The quantities µ T3,Q,B−L can be obtained by imposing n T3 = n Q = n B−L = 0. For temperatures above the critical temperature for the electroweak phase transition, we obtain the following chemical potentials, taking all the q ni = Q n . We subsequently obtain a baryon number density, where g SB = 3 Q n /4 and finally for its ratio with the entropy density, The equilibrium distribution changes after electroweak symmetry breaking, when there is no longer a need to conserve T 3 . This gives µ Q = −(4/11)µ B−L = −Q nφ /12f and, once again, g SB = 3Q n /4. However, species such as the RH fermions, which are coupled very weakly to the thermal plasma, would not be able to reequilibrate on the time scale of the electroweak phase transition. The precise value of g SB is thus hard to compute without considering the full dynamics of the process and may be different from the value obtained above. Given that these subtleties only lead to an O(1) ambiguity in g SB , for definiteness we stick to the value derived in Eq. 15. The value of η is frozen at T = T sph = 130 GeV, the temperature at which the sphaleron processes decouple [23]. Requiring that this reproduces the observed baryon asymmetry, η 0 = 8.7 × 10 −11 , we obtain, where we have used V ∼ 5Hφ ∼ Λ 4 c /f k via Eq. 6. To obtain the numerical value above we have taken the value of g SB in Eq. 15 with Q n = 6, where the latter choice will be justified in the next subsection. It is crucial that the relaxion keeps rolling with a non-zeroφ when the value of η gets frozen at T = T sph . To ensure this, we require that the critical temperature for the phase transition of the strong sector, T c , is lower than T sph .
Note that in the above analysis we have assumed that the RH neutrinos are relativistic and in equilibrium at the temperatures relevant for the calculation above. First of all, this requires that the seesaw scale, where M ij n = M n . The second requirement -that the RH neutrinos are in equilibrium -implies that the interaction rate of n i with SM particles satisfies, As discussed in detail in Ref. [24], taking into account all the processes that contribute to the equilibration of the RH neutrinos, we dimensionally expect, Γ(n) ∼ g 2 Y 2 N T , where g is the weak coupling and Y ij n ∼ Y n . Requiring Γ > H(T sph ) gives,

Neutrino masses
The Lagrangian in Eq. 10 generates masses for the SM neutrinos, Given that spontaneous baryogenesis demands M n < ∼ T sph , we require small effective Yukawas Y n < ∼ 10 −6 , which can be naturally obtained due to the Froggatt-Nielsen mechanism, as explained above. Note that without the Majorana mass term the neutrinos would only have a Dirac mass, which would require a much smaller Yukawa coupling. Our explanation for the smallness of neutrino masses thus utilises a combination of the Froggatt-Nielsen and seesaw mechanisms. Eq. 18 and Eq. 20 imply that our model requires a finite range for both the Yukawa coupling and Majorana mass scale, While sterile neutrinos with masses below 500 MeV are in tension with big-bang nucleosynthesis (BBN), for masses around a few GeV they might not be too far from the reach of future experiments such as SHiP [28]. For F N = 0.1 the above range of the Yukawa couplings can be obtained for 6 ≤ Q n ≤ 8, where we have taken all q ni = Q n . We have chosen equal charges for all the RH neutrinos because our objective is just to illustrate how the smallness of the neutrino mass scale can be explained. While this will result in an anarchic neutrino mass matrix (for an anarchic y ij n ) more structure in the masses and mixings can be obtained, for instance, by also charging the lepton doublets under the Froggatt-Nielsen symmetry (for an example see Ref. [3]). Finally note that as all the couplings in Eq. 11 are real, a concrete prediction of our model is a CP-conserving phase in the neutrino mass matrix which is still well within the range allowed by neutrino experiments [25].

IV. NELSON-BARR SECTOR AND THE
ROLLING POTENTIAL

Generating the Rolling Potential
To generate the rolling potential the abelian symmetry at the last site of the clockwork, U(1) N , needs to be broken. In Ref. [8] it was shown that this can be achieved with a minimal modification of the SM up sector, namely the addition of a vector-like pair, (ψ, ψ c ), where ψ has the same quantum numbers as an up-type singlet. Furthermore, as shown in Ref. [26], the same modification can also give a Nelson-Barr solution [14,15] to the strong CP problem, provided we impose an additional Z 2 symmetry. The Lagrangian terms for the relevant interactions are, where the ψ, ψ c and Φ N are odd under the Z 2 symmetry. The Z 2 symmetry forbids the term QHψ c . Recall that an exact CP symmetry has been imposed and all the couplings appearing in the Lagrangian, here and elsewhere, are real. It is clear that the U(1) N symmetry is collectively broken by y i ψ andỹ i ψ . This leads to the breaking of the relaxion shift symmetry and the generation of the rolling potential in Eq. (1),(2) with, The 1-loop Φ N → Φ N diagram gives the first term whereas the Φ 2 N H † H box diagram gives the second term. For the loop diagram generating the first term, we have taken the the cut-off for the ψu c loop to be the mass of the clockwork radial modes, m ρ ∼ f .

Nelson-Barr solution to the strong CP problem
We now show that the Lagrangian in Eq. (23) also provides a solution to the strong CP problem. Once the relaxion stops, the phase, θ N = π N /f ∼ φ/F , enters the 4 × 4 matrix for the up sector, where, The phase, θ N , is nothing but the phase of the cosine of the rolling potential at the relaxion stopping point. Note that the bottom-left element in the above mass matrix is zero due to the absence of the QHψ c term in the Lagrangian, which in turn is a direct consequence of the Z 2 symmetry. This ensures that at tree level there is no contribution toθ QCD from the phase, θ N , as, Arg(det(M u )) = Arg(µ ψ · det(vY u )) = 0, where we use the fact that µ ψ is real. On the other hand, for µ 2 + B i B * i v 2 we can integrate out the vector-like pair to give the following effective 3 × 3 mass squared matrix of the SM up quark sector: For | y ψ × ỹ ψ |/| y ψ + ỹ ψ | 2 1 there is an O(1) phase in the second term above and thus in the unitary matrix, V u L , diagonalising the above matrix. This is translated into an O(1) phase in the CKM matrix V CKM = V † u L V d L . In general, the delicate structure of the mass matrix, M u , is spoiled by radiative effects. For instance, these effects might induce a non-zero QHψ c term, or a phase in the diagonal terms, inducing a finiteθ QCD . Such effects For each value of f , taking the maximal value y ψ = 10 −2 , the value of the Higgs mass cut-off, M , is fixed as shown on the right-hand side of the frame. The red band shows the region where the cut-off exceeds the bounds imposed by cosmological requirements; blue denotes the region where the Higgs-independent contributions to the backreaction are no longer subdominant; the orange region shows the fifth force exclusions due to mixing of φ with the Higgs; green denotes the region in which the relaxion overshoots the backreaction barriers after reheating and the purple dashed lines show the value of tan(φ0/f ) (which can also be interpreted as the required amount of tuning) needed to reproduce the correct relic density. Finally, the grey band at the bottom shows the region where the relaxion mechanism is unable to raise the Higgs cut-off beyond 2 TeV. For more details see Sec. V.
were carefully analysed in Ref. [8] where it was shown thatθ QCD is within experimental bounds provided, whereỹ i ψ ∼ y i ψ ∼ y ψ . Before going to the next section we would like to comment that there appears to be no obvious difficulty in extending our model along the lines of Ref. [9] to also address the SM flavour puzzle for the charged leptons and quarks. At the cost of complicating our model, this can be achieved by identifying one of the intermediate sites of the clockwork chain with the flavon for the charged fermions and the abelian symmetry at this site with a Froggatt-Nielsen flavour symmetry. In order not to generate aθ QCD , the charge assignment of the SM fermions must be anomaly-free with respect to QCD as emphasised in Ref. [9] where an example charge assignment was also presented. We do not explore this direction further and stick to our more minimal set-up here.

V. PARAMETER SPACE
In this section we impose the constraints derived throughout the previous sections on the (m φ , f ) plane. The results are shown in Fig. 2. We fix the f k /f ratio according to Eq. 17 such that each point in the plot gives the correct baryon asymmetry. First, let us consider the constraints that arise on the SRB scenario for the right-handed neutrino current introduced in Sec. III. The vertical green band shows the region that is ruled out by requiring, m φ < 5H(T c ), the condition in Eq. 6 that the relaxion does not overshoot the barriers of the backreaction potential once they reappear after reheating. Here we have taken the maximal value T c = T sph (see Sec. III). The blue shaded region corresponds to the region Λ 2 c > 16π 2 v 2 ruled out by the requirement (derived in Ref. [7]) that the Higgs-dependent parts of the backreaction potential dominate over any Higgs-independent contribution. The dashed lines show the required value of tan(φ 0 /f k ) to reproduce the correct dark matter density in Eq. 8. As explained below Eq. 8, the extent to which tan φ0 f deviates from unity can be interpreted as a measure of the required tuning. The orange region shows fifth force constraints that arise due to the fact that the relaxion mixes with the Higgs boson with a mixing angle (see Ref. [27]), The rest of the constraints arise from the implementation of the Nelson-Barr mechanism in Sec. IV. First of all, from Eq. 24 we see that for a given value of f and y ψ one can fix the value of the Higgs mass cut-off , M , giving us the scale on the right hand side of the frame. The red band at the top shows the region where the value of the cut-off, M , exceeds the upper bound imposed in Eq. 5. Finally the grey band at the bottom shows the region where the relaxion mechanism is unable to raise the Higgs cut-off beyond 2 TeV.
We see from Fig. 2 that after all the constraints are imposed, a finite allowed region remains that remarkably contains the region in which tuning to obtain the correct relic density for dark matter is minimal. We will comment in the next section on how this allowed region can be probed further by future experiments.

VI. DISCUSSION AND CONCLUSION
We have presented a simultaneous solution to five BSM puzzles: the hierarchy problem, dark matter, matterantimatter asymmetry, neutrino masses and the strong CP problem. While our construction is admittedly more involved than some other attempts to solve BSM puzzles in a unified way [1][2][3][4], this is because we also use cosmological relaxation to achieve the challenging task of solving the hierarchy problem without adding any TeV scale states that cut-off the Higgs mass divergence. Our construction has all the ingredients of a standard relaxion model -such as a chain of clockwork scalars and a TeV-scale strong sector -but beyond this we only make minimal modifications by adding three RH neutrinos and an up-type SU(2) L singlet vector-like quark pair.
We utilise the fact that the standard relaxion mechanism has many interesting features built-in, such as spontaneous CPT violation during its rolling, spontaneous CP violation when it stops and oscillations about its stopping point after reheating. Motivated by previous work, we use these features to solve other BSM puzzles. The spontaneous CPT violation leads to spontaneous baryogenesis during the rolling of the relaxion after reheating; the spontaneous CP violation leads to a Nelson-Barr solution of the strong CP problem; and the relaxion oscillations generate the observed dark matter abundance. The spontaneous baryogenesis mechanism requires that baryons and/or leptons are charged under the relaxion shift symmetry. In this work the relaxion shift symmetry is identified with a Froggatt-Nielsen symmetry, under which three new RH neutrino states (but no SM states) are charged. This satisfies the requirement of spontaneous baryogenesis while also giving an explanation for the smallness of neutrino masses for a seesaw scale below the weak scale.
Our all-in-one relaxion set-up gives a diverse set of observational predictions. Our construction predicts the absence of a CP violating phase in the neutrino mass matrix and GeV scale sterile neutrinos that might not be too far from the reach of future experiments such as SHiP [28]. The strong CP phase in our model is non-zero and may be detectable in future experiments. The finite allowed parameter space in Fig. 2 can also be probed by future improvements in fifth force experiments. Finally, we would like to point out the interesting trade-off that exists in Fig. 2 between the Higgs mass cut-off scale up to which the hierarchy problem can be solved and the tuning required to reproduce the correct relic abundance of dark matter. The least tuned regions correspond to cut-off values smaller than 100 TeV. Top partners in a full solution to the hierarchy problem can thus be expected to be seen at these mass scales in a future 100 TeV collider.
We conclude by mentioning an interesting future direction. Our construction involves the standard relaxion mechanism, which utilises Hubble friction to stop the relaxion. It will be interesting to see if some of our ideas can be implemented in alternative models involving particle production [29][30][31][32], which are attractive because they decouple the relaxion mechanism from inflation.