We explore the possibility that the electroweak phase transition happens at a scale much higher than the electroweak scale today. In this context, high scale CP-violating sources for electroweak baryogenesis are not constrained by low-energy experiments. We propose a scenario of high-scale electroweak baryogenesis linked to flavour physics. This scenario allows for a period of enhanced Yukawa couplings during the evolution of the universe, which source time-dependent CP violation. The electroweak symmetry is never restored after the high-scale phase transition due to negative contributions to the Higgs thermal mass squared from a large number of additional electroweak-scale neutral scalars coupling to the Higgs. As a result, the washout of the baryon asymmetry is avoided.
July 16, 2019
High Scale Electroweak Phase Transition:
Baryogenesis & Symmetry Non-Restoration
Iason Baldes and Géraldine Servant
DESY, Notkestraße 85, D-22607 Hamburg, Germany
II. Institute of Theoretical Physics, University of Hamburg, D-22761 Hamburg, Germany
- 1 Introduction
- 2 Toy Example
- 3 Complete model
- 4 Phase transition and evolution of the potential
- 5 Earlier cosmological history
- 6 Phenomenology
- 7 Conclusions
- A Coefficients of the mass matrix
- B Phase transition with constant Yukawas
Two leading theories for explaining the matter-antimatter asymmetry are leptogenesis  and electroweak baryogenesis (EWBG) [2, 3]. The former is intrinsically tied to the mass scale of the Majorana neutrinos, , typically taken to be at a very high scale. In standard, non-flavoured, thermal leptogenesis with a hierarchy between the lightest and heavier sterile states, the Davidson-Ibarra bound requires GeV . Nevertheless, detailed work has shown that the mass scale can be reduced, even down to GeV, when flavour effects and the possibility of resonant enhancements with quasi-degenerate spectra are taken into account .
In contrast, it is usually assumed EWBG is intrinsically tied to the scale of electroweak symmetry breaking GeV. This is positive from the viewpoint of testability, as the required CP violation can be constrained from measurements of electric dipole moments (EDMs) , and the requirement for a strong first order phase transition typically leads to sufficient modifications of the zero-temperature Higgs potential for deviations to be measurable at colliders in the near future [7, 8, 9, 10]. This has led to a healthy tightening of the constraints on the scenario in the past few decades. Recently, the idea of linking EWBG to the flavour sector has been advocated. This has the advantage of possibly: (i) providing the CP violation required for EWBG while making the CP violation time dependent [11, 12, 13, 14, 15, 16, 17, 18, 19, 20] and hence naturally evading the EDM bounds, (ii) making the phase transition strong through the varying Yukawa couplings [14, 15]. Bringing flavour constraints into the game, however, makes model building in this framework challenging. It would therefore be helpful to raise the scale of EWBG, so we can in turn also raise the flavour scale and hence more easily satisfy the flavour constraints.
More broadly, raising the scale of EW symmetry breaking is anyway an exciting theoretical possibility, not limited to the context of the flavour model considered below. The aim of this paper therefore is to study the possibility of high scale EWBG, in which the Higgs first obtains a large vacuum expectation value (VEV), which is later gradually decreased to GeV while in the broken electroweak phase. The VEV can be gradually decreased using a symmetry non-restoration effect, in which the Higgs — through the coupling to other scalar fields — gains a negative thermal mass squared and hence a VEV proportional to the temperature [21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32].111For brevity, we omit “squared” when discussing the thermal masses of scalar fields from now on. In the models of symmetry non-restoration considered so far, the broken symmetry is not restored at any temperature. For electroweak baryogenesis, however, we want the Higgs to start in the symmetric phase and undergo a phase transition into the broken phase. Here, we will first show the two conditions can be realised together generically, through a simple toy model example, sketched in Fig. 1.
Motivated by our findings, we then return to flavour considerations in a more complete model, in which the Yukawa couplings are field-dependent and large at early times. The flavor sector contains extra fermions whose mass is controlled by the VEV of a scalar field that sets the flavour scale, (10) TeV, today. The broken EW phase minimum develops at large Higgs values once the temperature drops to the flavour scale. The Higgs then undergoes a strong first order phase transition from a point in field space in which the Yukawa couplings are (1) — allowing for enhanced CP violation compared to the SM — into the broken phase minimum [11, 12, 13, 14, 15, 16, 17, 18, 19]. This is when baryogenesis takes place.
Through another phase transition the Yukawa couplings are suppressed to their present day values. The Higgs also obtains a negative thermal mass and the VEV of the Higgs gradually decreases to GeV as the temperature drops. The washout avoidance condition, , is maintained throughout the evolution of the potential following the first phase transition. The sequence of phase transitions is summarised in Fig. 2. The scenario is a novel realisation of high scale electroweak baryogenesis linked to flavour physics.
The working model we present here is intended as an initial exploration of such ideas, through which model building problems can be identified and hopefully serve as inspiration for future work in more realistic contexts. The scenario shares some aspects of the leptogenesis scenario in Ref. . Though in that model — as the symmetry non-restoration effect is not used — EW symmetry is restored in an intermediate step. Hence, as was noted by the authors, it would require an additional source of number violation if it were to be modified for EWBG.
The paper is organised as follows. In the next section we illustrate the idea of a high scale phase transition or crossover using the symmetry non-restoration effect in a toy model. In Sec. 3 we describe the full model, starting from the fermionic sector and flavour structure and then moving onto the scalar potential and the symmetry non-restoration effect. In Sec. 4 we discuss the details of the phase transition. In Sec. 5, we briefly comment on a possible earlier cosmological history to justify our choice of potential parameters. We then discuss related phenomenology in Sec. 6, namely possible gravitational signals and the constraints on the low mass scalar states required by the symmetry non-restoration effect, before concluding our discussion.
2 Toy Example
High temperature symmetry non-restoration was studied some time ago [21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32], mainly in the context of GUT theories or in the context of SUSY flat directions . The phenomenon has been confirmed by lattice simulations [35, 36] and non-perturbative methods . For the electroweak symmetry, it was considered only a few times. The possible existence of a broken phase of electroweak symmetry at high temperature in Little Higgs extensions of the Standard Model was investigated in [38, 39]. The theory, however, exhibits a restoration of electroweak symmetry as long as temperatures are not pushed beyond the range of validity of the EFT for a finite temperature calculation . This conclusion is generalised to Twin Higgs models in  and confirms earlier findings in . The case of composite Higgs models with partial fermion compositeness in which the Higgs is a PNGB has been studied recently in and these models also lead to EW symmetry restoration [18, 19].
Here we will implement the ideas illustrated in Fig. 1, and show how a phase transition or crossover can occur at a high scale, i.e. above the zero-temperature minimum of the scalar potential, using an extension of the symmetry non-restoration effect. Unlike in earlier realisations of the symmetry non-restoration effect, the symmetry is actually restored at a sufficiently large temperature, i.e. above some mass threshold. Here, by symmetry non-restoration, we mean that at temperatures below the phase transition one of the scalar fields obtains a VEV proportional to the temperature.
The main idea is to induce a negative thermal mass for the Higgs through a negative cross-quartic coupling between the Higgs and a large number of additional scalar fields. Consider a toy model of scalar fields, , , and , where is a generational index (the reason for considering multiple generations will be made clear below). We denote the degrees of freedom with , , and (the sector therefore has in total degrees of freedom). In this section is acting as a placeholder for the EW Higgs, though we switch off the usual SM Yukawa and gauge interactions for the discussion in this section. For the purposes of our example, the relevant terms in the tree level potential are given by
where for simplicity we assume degenerate masses and couplings for the generations and that the cross quartic is negligible. As we shall be choosing , stability of the tree level potential requires
At high temperatures, , the thermal masses of the fields are 
Now consider a judicious choice of parameters so that: (i) and always have positive thermal masses, (ii) is positive at high temperature, (iii) becomes negative when the contribution of to its thermal mass becomes negligible, i.e. once . The effective potential in the direction, when can be approximated as . Positive returns a minimum at , but for negative we will find a minimum at . The latter solution is the usual symmetry non-restoration effect [21, 22, 23, 24, 25, 26, 27, 28, 29, 32]. What is new here is the presence of the additional field which can switch the sign of when reaches a mass threshold, leading to a phase transition or crossover. (Similarly, the symmetry non-restoration effect disappears if falls sufficiently below .) Eventually, for , the VEV is set by the usual zero-temperature minimization conditions.
We numerically evaluate the effective potential including the tree-level terms, zero and finite-temperature one-loop terms, and the daisy resummation.222We use the Arnold-Espinosa method of implementing the daisy resummation . We cut off the contribution of to the thermal masses with an exponential factor, , in order to avoid spurious contributions to the daisy resummation. We checked that the thermal mass estimated using the high-temperature expansion is consistent with the second derivative of the one loop thermal terms. In fact, the phase transition is stronger when using the numerical value rather than the high-temperature expansion value. The latter is crucial and weakens the phase transition. To give a concrete example, consider the choice of parameters333Motivated by flavour bounds, we take a characteristic scale TeV for illustration. The scale of the transition, however, can be taken much larger. The main limit for baryogenesis is around GeV when the sphalerons become out-of-equilibrium in the symmetric phase.
In Fig. 3 we show the resulting cross over, together with the thermal contributions from the and scalars and the daisy resummation. In Fig. 4 we plot the evolution of the VEV and effective mass of as a function of , showing the various stages discussed above. As mentioned previously, the mass threshold is naively at , however, additional factors which enter the full expressions lead to the non-zero VEV only developing at in our example. We have checked the VEVs remain zero throughout due to positive thermal contributions in the field directions.
The reason for requiring multiple generations of is revealed by considering the thermal mass of the , Eq. (3). A large thermal mass spoils the symmetry non-restoration effect once it enters the effective potential through the daisy resummation . This is because a large thermal mass can make the vacuum contribution, , which leads to the symmetry non-restoration effect, negligible in the effective potential. (This is not captured in the naive Eq. (5) which is simply based on a high- expansion.) Assuming, as we do, that , the use of multiple generations means the thermal mass of the can be reduced, assuming the inter-generational interactions are negligible. Thus allowing for the symmetry non-restoration phase to proceed even once the daisy resummation is included. Furthermore, the use of multiple generations allows us entertain the possibility that the are singlet fields, i.e. , in our full model below, which leads to simpler low energy phenomenology.
It is interesting that the stability constraint implies
which reveals that a negative thermal mass can be achieved for a sufficiently large , while keeping small enough, and the potential stable.
In our example here, we do not have a first order phase transition required for EWBG. Nevertheless, we shall see below that in our full model a strong enough phase transition can be achieved. What is important here is that we can start in the symmetric phase at high temperature and make a transition to a period in which obtains a large temperature-dependent VEV. We can then use the additional freedom gained, e.g. by introducing additional field directions, to arrange for a strong first order electroweak phase transition at a high scale followed by the use of the symmetry non-restoration effect to avoid washout.
3 Complete model
To realise the sequence of events described in the introduction, we need two main ingredients: First, there should be negative contributions to the Higgs thermal mass to enable the Higgs potential to exhibit a period of symmetry non-restoration. This can be achieved by adding extra scalar fields coupling to the Higgs as we have just seen in the previous section. Second, for a successful implementation of the EWBG mechanism, we need to start in the symmetric phase, at the highest temperatures. As we will see, we realise this through thermal effects from additional fermions together with higher dimensional terms in the potential, rather than scalar degrees-of-freedom as in the toy example.
In the following, we start discussing the new fermionic degrees of freedom connected to the flavour sector, as motivated in Froggatt-Nielsen models . These will provide positive thermal contributions to the effective potential and will also be responsible for creating a minimum in the Higgs potential at large Higgs values, through zero temperature one-loop effects (hence related to step 1 and 2 in Fig. 2).444Although the fermions create a minimum at large field values, we cannot trap the EW Higgs at this point to temperatures below the EW scale today without: (i) huge fine tuning of the polynomial potential, (ii) diluting the baryon asymmetry due to the false vacuum energy becoming dominant. Hence the symmetry non-restoration effect is still required to avoid washout. Furthermore, these fermions help us achieve a strong first order phase transition and are also our source of CP violation [14, 15].
We then discuss the extended scalar sector consisting of the EW higgs together with a scalar which controls the masses of the exotic fermions. The phase transition in the direction is essential in reducing the effective Yukawa couplings to their SM values (step 2). Finally we discuss the new scalars leading to EW symmetry non-restoration, which provide a negative thermal contribution for the first phase transition (step 1), and also enable the final stage of the mechanism (step 3). Combining all these effects, we obtain the more complicated sequence of phase tranisitions sketched in Fig. 2, in contrast to the transition in a single field direction as in the toy model.
3.1 The fermionic sector
To illustrate our scenario we focus on the top and charm quarks using the Froggatt-Nielsen (FN) mechanism. The mass matrix follows the pattern,
where . In the FN picture where is the FN scale set by the mass of vector-like quarks and is the flavon VEV or an explicit soft breaking of the FN symmetry by one unit. It is necessary to explicitly break the FN symmetry, unless the FN symmetry is gauged in an extended model, in order to avoid the appearance of a massless Goldstone boson. Hence, we shall assume below that arises from an explicit breaking, in order to avoid having to study the dynamics of the flavon field.555The flavon will eventually gain a VEV and — if it is of the same order as the other dimensionful terms in the flavour sector — it will also be of the same order as the explicit breaking scale. Hence it is not expected to change our overall picture. This means the scalar and pseudoscalar components of the flavon will also end up with masses at a similar scale. To be safe from limits from mixing this mass scale should be at least several TeV if the flavon couples to all quark flavours . To generate the above mass matrix we assign the following FN charges to the SM quarks:
In the UV completion we add vector-like quarks which transform as under the SM gauge group. We require three such vector-like quarks,
each with twelve degrees-of-freedom, where and denote the usual chiral components and the number in the superscript denotes the negative FN charge. The full mass matrix is then given by
where we have suppressed factors of and indicate bare mass terms allowed by the FN and electroweak symmetries by . The entries proportional to break the FN symmetry by one unit. In principle there may be even smaller entries in the terms which break the FN symmetry by more than one unit, which may be generated by renormalization group flow. For simplicity we assume these are negligible and set the corresponding entries to zero. We next imagine the mass terms as arising from a bare contribution, which we take to be , and through the Yukawa coupling to another scalar in the form . The full mass matrix is therefore given by
The eventually obtains a large VEV, , giving . It is useful to define effective Yukawa couplings
The effective Yukawa couplings to the Higgs and the fermion masses along the axis are shown in Fig. 5. At the zero temperature minimum (where , not shown in in Fig. 5), we obtain three super heavy mass eigenstates, , corresponding to the FN fermions, a mass eigenstate corresponding to the top, , and one corresponding to the charm, . Below we shall study the temperature evolution of together with the Higgs field . The numerical values of the coefficients in the fermionic mass matrix used in our analysis are given in Appendix A.
3.2 The scalar potential
In this section we describe the two-field scalar potential consisting of the real scalar and the electroweak Higgs . We write the tree-level potential as
We fix the electroweak Higgs mass and VEV to the observed values, GeV and GeV. The higher dimensional terms are required to stabilise the potential due to the large number of fermions with sizable Yukawa couplings. In particular should not be too far above the scale of the FN fermions or there would be an instability in the direction. Note that generically, to achieve GeV and GeV requires a large degree of fine-tuning already in the tree level terms of the potential as soon as , which is required from flavour constraints. This is seen in the relation for the physical Higgs mass
where the large terms on the right-hand-side must be tuned to return the much smaller . Here our philosophy is to assume these large tree-level contributions cancel, possibly due to a Higgs relaxation mechanism operating during an earlier period, as discussed in Sec. 5. In our example we choose the parameters to be
As we shall see below, the dimensionless couplings have been chosen as to obtain the required pattern of symmetry breaking in the thermal evolution of the potential. With the parameters chosen above we find . The effective electroweak quartic, i.e. , which enters the triple Higgs cross section for collider searches, remains close to its SM value. The potential is shown on the left in Fig. 6. The one-loop potential, including the effect of the fermions after diagonalising (12), is shown on the right in Fig. 6. As can be seen in the figure, although the tree-level potential has a barrier in the direction, this is almost completely erased at loop-level once the effect of the fermions is taken into account. The appropriate strong first order phase transition can be achieved by an interplay of the fermionic degrees-of-freedom together with the symmetry non-restoration sector which we discuss next.
3.3 Symmetry non-restoration sector
As in the toy example, we obtain a negative contribution to the thermal mass of through negative cross quartics. Let us introduce generations of singlet scalars , i.e. each with degree of freedom. The relevant terms in the tree level potential are
where we again assume for simplicity that the couplings are universal and that any interactions between the and are negligible (also in order to keep small enough to not spoil the symmetry non-restoration effect). At high temperatures the thermal masses of the fields are
where we have introduced the sum over the effective, field-dependent, Yukawa couplings of the Higgs. As the SM contributions to already amount to , we require a large in order to obtain a negative thermal mass for while remaining consistent with the stability constraint.666The stability constraint is relaxed in the presence of the higher dimensional operators. Here we choose parameters sufficient for stability of the potential, i.e. consistent with stability in the limit . Here, for illustration, we choose
Note that the size of the symmetry non-restoration effect depends on the field values and , through the effective Yukawas . This is illustrated in Fig. 7. Nevertheless, this does not mean there is necessarily a minimum for non-zero and at all temperatures, because of (i) the higher dimensional terms and because (ii) the FN fermions also couple to raising the potential due to finite temperature effects. Because some of the fermions already have masses of the order of the critical temperature in the symmetric phase, , a simple use of Eq. (19) is not possible here, and a numerical evaluation of the effective potential is required.
4 Phase transition and evolution of the potential
We now describe the various effects contributing to the sequence of phase transitions in this model. As advertised above, temperature effects from the additional FN fermions play a crucial role in maintaining the global minimum of the potential at and at high temperatures. Eventually drops, a non-trivial interplay between the one-loop terms for the fermions, assisted by those for the EW gauge bosons, and the symmetry non-restoration effect allows for a broken phase minimum to develop. Depending on the parameters chosen, we find we can obtain a first order phase transition along the direction, followed by another first order phase transition in the direction. Afterwards there is a slow evolution to the zero temperature minimum, at large and small , all the while maintaining .
In order to avoid an early transition along the direction, which leads eventually to a cross over in the direction due to the negative thermal mass at large values, the quartic should be small enough. The mass parameter must be at most EW scale, in order to maintain a negative down to , for simplicity we have set it to zero for our plots in this section. We discuss relaxing this assumption in Sec. 6.2. In order to obtain a strong first order phase transition in the direction, we require to be below some value, otherwise the phase transition occurs too early. On the other hand, to maintain a large enough during the subsequent evolution, we require a large enough . Keeping all other parameters fixed, we find the correct evolution of the potential for , when we set GeV (which is relevant for the low analysis).
A detailed plot of the potential at the critical temperature is shown in Fig. 8. Note the interplay between the symmetry non-restoration effect — arising from loop effects of on the potential — and the positive fermionic, gauge and daisy terms leads to the barrier. We wish to emphasise the particular importance of the varying Yukawas in achieving a strong first order phase transition in our example . The phase transition, in contrast, is much weaker if we switch off the Yukawa variation effects. This is discussed in further detail in Appendix B. The overall evolution of the potential for our choice of parameters is shown in Fig. 9. We have calculated the symmetric bounce action for the bubble, denoted , using the AnyBubble code . The probability of nucleating a bubble in a Hubble volume reaches in a radiation dominated Universe when 
In our example we find the step 1 phase transition occurs at TeV. Here the washout parameter reads
After remaining in the TeV minimum and supercooling to TeV, we find the step 2 phase transition occurs. The path of the two-field bounce solution is shown in Fig. 9. Note the initial bubble is thick-walled, meaning the centre of the bubble is away from the true minimum in field space. Nevertheless, as the bubble expands the fields will quickly relax down to the minimum of the potential. We have checked that with the given parameters the minimum does indeed respect until the EW minimum is reached. We have also checked that the positive thermal contributions in the field direction are sufficient to keep the VEVs at zero throughout.
5 Earlier cosmological history
We now comment on the hierarchy problem in this framework. As shown in Eq. (15), we need to tune parameters to keep the Higgs mass parameter small. One way to address this is to stipulate that a relaxion mechanism took place before the EW phase transition . The relaxation of the Higgs parameter would have to take place during inflation. We would then have the following cosmological history:
Inflation begins, drastically lowering the temperature of the thermal bath. At this stage, the Higgs mass parameter and the mass parameter are both large.
Relaxation starts for the Higgs and its mass parameter, Eq. (15), is relaxed to the usual low weak scale value. Relaxation ends when the Higgs obtains a small VEV, leading to the backreaction on the relaxion potential. The VEV of does not contribute to the relaxion potential barriers and hence it can naturally be of a larger scale.
Reheating in the visible sector: EW symmetry is restored and is also reheated so its VEV goes back to zero.
The evolution described in Sec. 4 takes place.
In this context, we would assume that the scalar sector is also relaxed during stage 3. So the Higgs and are part of a common sector (perhaps composite) and they are relaxed together, while from the flavour sector is not subject to relaxation.
6.1 Gravitational wave signal
During the cosmological evolution, after the EW phase transition and baryogenesis, between step 1 and step 2, the scalar fields become stuck at a false minimum and there is some super cooling. The timescale of the transition is
The ratio of energy released compared to radiation bath, however, is rather suppressed
due to the contribution of to . Thus the resulting stochastic gravitational wave background generated during the first-order phase transition  is suppressed. It is too small to be detected by LISA  but it is within the BBO sensitivity , as illustrated in Fig. 10. Due to the unusual situation of a relatively low combined with a suppressed , the bubble wall collisions (envelope contribution), gives the dominant effect. The step-1 phase transition is characterised by and , and returns a completely negligible gravitational wave background.
On the other hand, if another cosmological gravitational wave background due to cosmic strings exists, the rapid and huge change in due to the at the EW scale leads to some feature at a characteristic frequency in the spectrum of gravitational waves emitted in the radiation dominated era . A similar feature can be expected in suitable gravitational wave backgrounds coming from inflation. We leave this study for future investigation.
6.2 Scalar sector in the IR
The scenario relies on the scalar degrees-of-freedom to guide the electroweak minimum to its present value. Hence, it is necessary for the mass to be at or below the EW scale otherwise, once , the symmetry non-restoration effect disappears and will become small. This is shown in Fig. 11. The experimental constraint on such a scenario comes from searches for these light scalars. Note while we have considered universal mass and coupling terms for the , we can imagine that in a more realistic scenario the masses are split in a spectrum of states with masses . The partial width of the SM Higgs to the is given by
where denotes the number of generations with mass below the threshold . Demanding at most an modification to the SM Higgs signal strength requires for . Hence the states must lie above this threshold. In summary, we obtain
by combining the EW Higgs decay constraint with the washout avoidance condition shown in Fig. 11.
The states will become thermally populated and should not over-produce DM. The cross quartic is too small for annihilation solely through the Higgs portal and anyway, at these masses, is ruled out by direct detection [54, 55, 56, 57, 58]. Hence we need to arrange for the to decay.777Alternatively, provided the additional interaction does not lead to a too large thermal mass, the could annihilate into dark radiation , or a dark mediator which subsequently decays . This can be achieved if the obtain VEVs and can mix with the Higgs. Here we assume the obtain VEVs by introducing a small explicit breaking of the symmetry in Eq. (17). Explicitly this may be introduced through a linear term in the scalar potential
which, in the limit induces VEVs
The mixing angle for the mixing of a singlet state with the Higgs is given by
The can then mix with the SM Higgs and decay into light SM degrees of freedom. We demand that the decay before their energy density grows to dominate the universe, as otherwise they would dilute the baryon asymmetry . This can be achieved provided the decay rate of the states, MeV , is larger than the Hubble rate, , when . This implies for GeV.
Further limits come from EW precision observables, Higgs signal strength measurements and direct searches [63, 64]. We may derive an approximate constraint by considering a degenerate spectrum, keeping in mind direct search limits will not apply once the masses are split in a more realistic model. Given a mass of the singlet states GeV, the limit on the sum of the mixing angles reads , depending on the mass [63, 64]. Together with the rapid decay condition, this bounds the mixing angle to lie in the range
Translated into a bound on the VEVs this reads
which shows Eq. (28) can be applied consistently. By introducing a mixing with the SM Higgs, we also open up a decay channels of the form , , and . Nevertheless, a calculation reveals that these are completely negligible. For example, from dimensional analysis,
A similar calculation reveals
We have confirmed these with a more detailed calculation, which also displays the additional suppression expected as and the available phase space is reduced. The four-body decay channels are even more suppressed.
It is usually thought that the EW phase transition occurs when the Universe cools to temperatures GeV. In this paper we have instead speculated on the possibility of high scale EW phase transition and EW baryogenesis. This requires additional field content in order to break the EW symmetry at a high scale and to also suppress the sphalerons to avoid washout while the EW VEV is lowered to its present day value. We first showed the generic ingredients required for a transition to occur at a temperature far above the scale of the zero temperature minimum of a theory. Such findings may be of more general interest.
We then moved onto our specific scenario. In our example we have demonstrated the combination of: (i) a flavour model, (ii) the symmetry non-restoration effect can give us a novel scenario of high-scale EW baryogenesis. Both the Froggatt-Nielsen fermions responsible for field dependent Yukawas and the non-restoring scalar degrees-of-freedom combine to give us a strong first order phase transition. The fermions also help to control the symmetry non-restoration effect. Furthermore, the large Yukawas during the phase transition can act as the source of CP violation required to obtain the baryon asymmetry. This naturally allows for an absence of measurable EDMs.
The generic prediction of the scenario is a large number of light scalars, around the EW scale, with a small coupling to the EW Higgs. In our scenario these scalars mix with the SM-like Higgs, although a more complete construction with an extended hidden sector may eventually show that this is not generically necessary. The model presented here may well not be the simplest or most elegant realisation of these ideas, it is presented as a proof-of-principle, we hope it facilitates further exploration of this intriguing possibility.
G. Servant thanks R. Rattazzi for raising a question after her talk at the 2017 Johns Hopkins workshop that stimulated this work, and is grateful to C. Csaki and D. Kaplan for organising the workshop. We thank N. Suresh for earlier collaboration on phase transitions in the context of FN models and J.R. Espinosa and T. Konstandin for discussions.
While this paper was completed, Ref.  appeared, which deals with similar ideas, although is not motivated by baryogenesis as the EW symmetry is never restored in their case. We also learnt a high scale EW phase transition is being considered by Glioti, Rattazzi, and Vecchi (to appear).
Appendix A Coefficients of the mass matrix
For completeness, we provide the coefficients of our generalised mass matrix, Eq. (12). These were found by generating uniformly-distributed pseudorandom numbers with magnitudes in the range and phases in the range . We found an initial seed returning approximately the top and charm masses after attempts. Some entries were then further adjusted by in order to return the mass eigenvalues,
at GeV, TeV and TeV, where the top and charm correspond to and respectively. The coefficients of the entries proportional to in matrix (12) are given by
The coefficients of the entries are given by
Finally, the coefficients of the entries are given by
The mass matrix , in Eq. (12), is then formed by summing the three matrices above, multiplied by the relevant field values or factor, and dividing by . We then diagonalised at discrete points in field space and then interpolated over these points for reasons of efficiency in the numerical work.
Appendix B Phase transition with constant Yukawas
To contrast with our analysis above, we now consider the Step 1 phase transition with the mixing terms in the fermionic mass matrix switched off, i.e. we set . The Yukawa couplings and masses of the fermions along the, , axis are shown in Fig. 12. There are three zero mass eigenstates and two with couplings to . We calculate the critical temperature of the phase transition and find a very weak first order transition at TeV. This is shown in Fig. 13, along with the potential at zero temperature, showing that the qualitative difference to the case is apparent once the finite effects are taken into account. For other proposals of using fermions to achieve a strong first order phase transition see [14, 15, 66, 67].
-  M. Fukugita and T. Yanagida, Baryogenesis Without Grand Unification, Phys. Lett. B 174 (1986) 45–47.
-  M. E. Shaposhnikov, Baryon Asymmetry of the Universe in Standard Electroweak Theory, Nucl. Phys. B 287 (1987) 757–775.
-  A. G. Cohen, D. B. Kaplan and A. E. Nelson, Baryogenesis at the weak phase transition, Nucl. Phys. B 349 (1991) 727–742.
-  S. Davidson and A. Ibarra, A Lower bound on the right-handed neutrino mass from leptogenesis, Phys. Lett. B 535 (2002) 25–32, [hep-ph/0202239].
-  E. K. Akhmedov, V. A. Rubakov and A. Yu. Smirnov, Baryogenesis via neutrino oscillations, Phys. Rev. Lett. 81 (1998) 1359–1362, [hep-ph/9803255].
-  ACME collaboration, J. Baron et al., Order of Magnitude Smaller Limit on the Electric Dipole Moment of the Electron, Science 343 (2014) 269–272, [1310.7534].
-  C. Grojean, G. Servant and J. D. Wells, First-order electroweak phase transition in the standard model with a low cutoff, Phys. Rev. D 71 (2005) 036001, [hep-ph/0407019].
-  D. Curtin, P. Meade and C.-T. Yu, Testing Electroweak Baryogenesis with Future Colliders, JHEP 11 (2014) 127, [1409.0005].
-  S. Di Vita, C. Grojean, G. Panico, M. Riembau and T. Vantalon, A global view on the Higgs self-coupling, JHEP 09 (2017) 069, [1704.01953].
-  S. Di Vita, G. Durieux, C. Grojean, J. Gu, Z. Liu, G. Panico et al., A global view on the Higgs self-coupling at lepton colliders, JHEP 02 (2018) 178, [1711.03978].
-  M. Berkooz, Y. Nir and T. Volansky, Baryogenesis from the Kobayashi-Maskawa phase, Phys. Rev. Lett. 93 (2004) 051301, [hep-ph/0401012].
-  G. Perez and T. Volansky, Split fermions baryogenesis from the Kobayashi-Maskawa phase, Phys. Rev. D 72 (2005) 103522, [hep-ph/0505222].
-  L. J. Hall, H. Murayama and G. Perez, Electroweak baryogenesis from late neutrino masses, Phys. Rev. Lett. 95 (2005) 111301, [hep-ph/0504248].
-  I. Baldes, T. Konstandin and G. Servant, A First-Order Electroweak Phase Transition in the Standard Model from Varying Yukawas, 1604.04526.
-  I. Baldes, T. Konstandin and G. Servant, Flavor Cosmology: Dynamical Yukawas in the Froggatt-Nielsen Mechanism, JHEP 12 (2016) 073, [1608.03254].
-  B. von Harling and G. Servant, Cosmological evolution of Yukawa couplings: the 5D perspective, JHEP 05 (2017) 077, [1612.02447].
-  S. Bruggisser, T. Konstandin and G. Servant, CP-violation for Electroweak Baryogenesis from Dynamical CKM Matrix, JCAP 1711 (2017) 034, [1706.08534].
-  S. Bruggisser, B. Von Harling, O. Matsedonskyi and G. Servant, The Baryon Asymmetry from a Composite Higgs, 1803.08546.
-  S. Bruggisser, B. Von Harling, O. Matsedonskyi and G. Servant, Electroweak Phase Transition and Baryogenesis in Composite Higgs Models, 1804.07314.
-  G. Servant, The serendipity of electroweak baryogenesis, Phil. Trans. Roy. Soc. Lond. A376 (2018) 20170124.
-  S. Weinberg, Gauge and Global Symmetries at High Temperature, Phys. Rev. D 9 (1974) 3357–3378.
-  R. N. Mohapatra and G. Senjanovic, Soft CP Violation at High Temperature, Phys. Rev. Lett. 42 (1979) 1651.
-  Y. Fujimoto and S. Sakakibara, On symmetry nonrestoration at high temperature, Phys. Lett. B 151 (1985) 260–262.
-  P. Salomonson, B. S. Skagerstam and A. Stern, On the Primordial Monopole Problem in Grand Unified Theories, Phys. Lett. B 151 (1985) 243–246.
-  P. Salomonson and B.-S. K. Skagerstam, High Temperature Phases in an O() X O() Symmetric Four Dimensional Vector Model, Phys. Lett. B 155 (1985) 100–102.
-  G. R. Dvali, A. Melfo and G. Senjanovic, Is There a monopole problem?, Phys. Rev. Lett. 75 (1995) 4559–4562, [hep-ph/9507230].
-  G. Bimonte and G. Lozano, Can symmetry nonrestoration solve the monopole problem?, Nucl. Phys. B 460 (1996) 155–166, [hep-th/9509060].
-  G. Bimonte and G. Lozano, On Symmetry nonrestoration at high temperature, Phys. Lett. B 366 (1996) 248–252, [hep-th/9507079].
-  G. R. Dvali, A. Melfo and G. Senjanovic, Nonrestoration of spontaneously broken P and CP at high temperature, Phys. Rev. D 54 (1996) 7857–7866, [hep-ph/9601376].
-  J. Orloff, The UV price for symmetry nonrestoration, Phys. Lett. B403 (1997) 309–315, [hep-ph/9611398].
-  M. Pietroni, N. Rius and N. Tetradis, Inverse symmetry breaking and the exact renormalization group, Phys. Lett. B397 (1997) 119–125, [hep-ph/9612205].
-  M. B. Gavela, O. Pene, N. Rius and S. Vargas-Castrillon, The Fading of symmetry nonrestoration at finite temperature, Phys. Rev. D 59 (1999) 025008, [hep-ph/9801244].
-  Y. Hamada and M. Yamada, Baryogenesis in false vacuum, Eur. Phys. J. C 77 (2017) 643, [1605.06897].
-  B. Bajc and G. Senjanovic, High temperature symmetry breaking via flat directions, Phys. Rev. D61 (2000) 103506, [hep-ph/9811321].
-  K. Jansen and M. Laine, Inverse symmetry breaking with 4-D lattice simulations, Phys. Lett. B435 (1998) 166–174, [hep-lat/9805024].
-  G. Bimonte, D. Iniguez, A. Tarancon and C. L. Ullod, Inverse symmetry breaking on the lattice: An Accurate MC study, Nucl. Phys. B559 (1999) 103–122, [hep-lat/9903027].
-  M. B. Pinto and R. O. Ramos, A Nonperturbative study of inverse symmetry breaking at high temperatures, Phys. Rev. D61 (2000) 125016, [hep-ph/9912273].
-  J. R. Espinosa, M. Losada and A. Riotto, Symmetry nonrestoration at high temperature in little Higgs models, Phys. Rev. D 72 (2005) 043520, [hep-ph/0409070].
-  S. Aziz, B. Ghosh and G. Dey, Broken electroweak phase at high temperature in the Littlest Higgs model with T-parity, Phys. Rev. D 79 (2009) 075001, [0901.3442].
-  A. Ahriche, The Restoration of the Electroweak Symmetry at High Temperature for Little Higgs, 1003.5045.
-  C. Kilic and S. Swaminathan, Can A Pseudo-Nambu-Goldstone Higgs Lead To Symmetry Non-Restoration?, JHEP 01 (2016) 002, [1508.05121].
-  A. Katz and M. Perelstein, Higgs Couplings and Electroweak Phase Transition, JHEP 07 (2014) 108, [1401.1827].
-  P. B. Arnold and O. Espinosa, The Effective potential and first order phase transitions: Beyond leading-order, Phys. Rev. D 47 (1993) 3546, [hep-ph/9212235].
-  C. D. Froggatt and H. B. Nielsen, Hierarchy of Quark Masses, Cabibbo Angles and CP Violation, Nucl. Phys. B 147 (1979) 277.
-  A. Masoumi, K. D. Olum and B. Shlaer, Efficient numerical solution to vacuum decay with many fields, JCAP 1701 (2017) 051, [1610.06594].
-  J. M. Moreno, M. Quiros and M. Seco, Bubbles in the supersymmetric standard model, Nucl. Phys. B 526 (1998) 489–500, [hep-ph/9801272].
-  M. Drees, F. Hajkarim and E. R. Schmitz, The Effects of QCD Equation of State on the Relic Density of WIMP Dark Matter, JCAP 1506 (2015) 025, [1503.03513].
-  S. Borsanyi et al., Calculation of the axion mass based on high-temperature lattice quantum chromodynamics, Nature 539 (2016) 69–71, [1606.07494].
-  P. W. Graham, D. E. Kaplan and S. Rajendran, Cosmological Relaxation of the Electroweak Scale, Phys. Rev. Lett. 115 (2015) 221801, [1504.07551].
-  C. Grojean and G. Servant, Gravitational Waves from Phase Transitions at the Electroweak Scale and Beyond, Phys. Rev. D75 (2007) 043507, [hep-ph/0607107].
-  C. Caprini et al., Science with the space-based interferometer eLISA. II: Gravitational waves from cosmological phase transitions, JCAP 1604 (2016) 001, [1512.06239].
-  E. Thrane and J. D. Romano, Sensitivity curves for searches for gravitational-wave backgrounds, Phys. Rev. D 88 (2013) 124032, [1310.5300].
-  Y. Cui, M. Lewicki, D. E. Morrissey and J. D. Wells, Cosmic Archaeology with Gravitational Waves from Cosmic Strings, Phys. Rev. D 97 (2018) 123505, [1711.03104].
-  XENON collaboration, E. Aprile et al., First Dark Matter Search Results from the XENON1T Experiment, Phys. Rev. Lett. 119 (2017) 181301, [1705.06655].
-  PandaX-II collaboration, X. Cui et al., Dark Matter Results From 54-Ton-Day Exposure of PandaX-II Experiment, Phys. Rev. Lett. 119 (2017) 181302, [1708.06917].
-  LUX collaboration, D. S. Akerib et al., Results from a search for dark matter in the complete LUX exposure, Phys. Rev. Lett. 118 (2017) 021303, [1608.07648].
-  I. Baldes, N. F. Bell, K. Petraki and R. R. Volkas, Two radiative inverse seesaw models, dark matter, and baryogenesis, JCAP 1307 (2013) 029, [1304.6162].
-  M. Duerr, P. Fileviez PÃ©rez and J. Smirnov, Scalar Dark Matter: Direct vs. Indirect Detection, JHEP 06 (2016) 152, [1509.04282].
-  I. Baldes and K. Petraki, Asymmetric thermal-relic dark matter: Sommerfeld-enhanced freeze-out, annihilation signals and unitarity bounds, JCAP 1709 (2017) 028, [1703.00478].
-  I. Baldes, M. Cirelli, P. Panci, K. Petraki, F. Sala and M. Taoso, Asymmetric dark matter: residual annihilations and self-interactions, SciPost Phys. 4 (2018) 041, [1712.07489].
-  R. J. Scherrer and M. S. Turner, Decaying Particles Do Not Heat Up the Universe, Phys. Rev. D 31 (1985) 681.
-  J. R. Ellis, M. K. Gaillard and D. V. Nanopoulos, A Phenomenological Profile of the Higgs Boson, Nucl. Phys. B 106 (1976) 292.
-  T. Robens and T. Stefaniak, Status of the Higgs Singlet Extension of the Standard Model after LHC Run 1, Eur. Phys. J. C 75 (2015) 104, [1501.02234].
-  A. Falkowski, C. Gross and O. Lebedev, A second Higgs from the Higgs portal, JHEP 05 (2015) 057, [1502.01361].
-  P. Meade and H. Ramani, Unrestored Electroweak Symmetry, 1807.07578.
-  M. Carena, A. Megevand, M. Quiros and C. E. M. Wagner, Electroweak baryogenesis and new TeV fermions, Nucl. Phys. B 716 (2005) 319–351, [hep-ph/0410352].
-  D. Egana-Ugrinovic, The minimal fermionic model of electroweak baryogenesis, JHEP 12 (2017) 064, [1707.02306].