ν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu $$\end{document} Electroweak baryogenesis: the scalar singlet strikes back

We perform a comprehensive scan of the parameter space of a general singlet scalar extension of the Standard Model to identify the regions which can lead to a strong first-order phase transition, as required by the electroweak baryogenesis mechanism. We find that taking into account bubble nucleation is a fundamental constraint on the parameter space and present a conservative and fast estimate for it so as to enable efficient parameter space scanning. The allowed regions turn out to be already significantly probed by constraints on the scalar mixing from Higgs signal strength measurements. We also consider the addition of new neutrino singlet fields with Yukawa couplings to both scalars and forming heavy (pseudo)-Dirac pairs, as in the linear or inverse Seesaw mechanisms for neutrino mass generation. We find that their inclusion does not alter the allowed parameter space from early universe phenomenology in a significant way. Conversely, there are allowed regions of the parameter space where the presence of the neutrino singlets would remarkably modify the collider phenomenology, yielding interesting new signatures in Higgs and singlet scalar decays.


Introduction
The origin of the observed Baryon Asymmetry of the Universe (BAU) is one of the fundamental open problems of the Standard Model (SM) of particle physics and one of the few precious pieces of experimental evidence for physics beyond the SM together with the existence of neutrino masses and dark matter.The generation of the BAU in the early Universe requires satisfying the three Sakharov conditions [1]: baryon number violation, C and C P violation, and departure a e-mail: rosauro@ijclab.in2p3.fr(corresponding author) from thermal equilibrium.In principle the SM itself could address the origin of the BAU via the electroweak baryogenesis (EWBG) mechanism [2][3][4][5][6][7].However, the SM threefamily quark mixing encoded in the Cabibbo-Kobayashi-Maskawa (CKM) matrix does not provide enough C P violation to generate a sufficient asymmetry [8][9][10], and the early Universe transition from the electroweak (EW) symmetric phase to the EW broken phase in the SM is a smooth crossover [11,12], rather than the strongly first-order transition required by the out-of-equilibrium Sakharov condition.
Nevertheless, simple extensions of the SM could solve these issues and make EWBG viable.In particular, extending the scalar sector of the SM by just a real singlet field could allow for a first-order EW phase transition (see e.g.[13][14][15][16][17]).This new scalar singlet may not be alone, but rather be part of an extended dark sector to which it couples.A particularly motivated such scenario is the extension of the SM by (heavy) fermion singlets, i.e. right-handed neutrinos, able to account for the observed pattern of neutrino masses and mixings in Nature.Remarkably, it was shown in Refs.[18,19] that the new sources of CP violation that arise in this extension of the SM, from the simultaneous presence of Yukawa interactions of the heavy neutrinos with the singlet scalar and with the Higgs doublet and SM neutrinos, could lead to successful EWBG depending on the evolution of the scalar sector during the phase transition (a scenario referred to as ν-EWBG in [19]).
In this work we aim to clarify the conditions on the singlet scalar dynamics during the EW phase transition that possibly allows for EWBG in the above setup.Our scope is however more general, and we study the regions of parameter space of the real singlet scalar extension of the SM yielding a strong first-order phase transition (SFOPT), exploring the correlations among different model parameters, and emphasizing those that might arise between measurable quantities.The aim here is not a high-precision computation of the various thermodynamic quantities of the phase transition, which would be numerically challenging if combined with a thorough scan of the model parameter space.Rather, we focus on exploring the parameter space as efficiently as possible, covering broad areas of the multidimensional space via a number of approximations.Even if these approximations are not suitable to obtain highly-accurate results for the SFOPT quantities, they allow to pinpoint the regions of the parameter space with the desired features and test whether they are presently allowed, for subsequent analyses to concentrate in these regions.We pay particular attention to the SFOPT requirement of bubble nucleation for a successful completion of the EW phase transition, for which we provide conservative and fast estimate for assessing if nucleation would take place.While the nucleation dynamics has been studied previously [20] in the context of a Z 2 -symmetric singlet scalar extension of the SM, here we aim at a more general study without the additional constraint of additional symmetries.To the best of our knowledge, this is the first full scan exploring all the different correlations of the parameter space of the scalar singlet extension of the SM aiming to identify the regions where a SFOPT could take place.In addition, we investigate the impact of the heavy neutrinos on the SFOPT dynamics: while previous studies indicate that sizable values of the neutrino Yukawa interactions with the scalar singlet can strengthen the first-order phase transition [21], we find that large values of the neutrino Yukawas, unless compensated by other parameters, can also have the effect of destabilizing the EW broken minimum and are thus generally disfavoured.However, their inclusion does not alter significantly the allowed regions of the parameter space as compared to the singlet-only case.
Finally, we also discuss the phenomenological impact of the existence of such heavy neutrinos as compared to the minimal singlet scalar extension of the SM, finding that the phenomenology can be altered dramatically with respect to the latter model.Specifically, we find that the singlet-like scalar will dominantly decay into right-handed neutrinos (if allowed by phase space), instead of directly decaying into SM particles.These heavy neutrinos may then subsequently decay into SM particles either promptly or via displaced vertexes, depending on the size of their mixing θ with the active SM neutrinos.Since the production of the right-handed neutrinos from the scalar singlet decay is unrelated to the strongly constrained mixing θ , this heavy neutrino production process could well be the dominant one at colliders like the LHC.Besides, the mixing between the scalar singlet and the Higgs would also lead to exotic Higgs decays into right-handed neutrinos (if these are light enough), which can be probed via Higgs signal strength measurements and also in direct searches for such exotic Higgs decays (see e.g.[22]) at the LHC.
This paper is organized as follows.In Sect. 2 we introduce the real singlet scalar extension of the SM with the addition of heavy neutrinos, and discuss the details of the scalar potential in the early Universe relevant for our SFOPT analysis.Then, in Sect. 3 we analyze the experimental constraints on the model, as well as the possible new combined probes of the existence of the singlet and the heavy neutrinos.In Sect. 4 we give details of our model parameter scan, and we discuss its results in Sect. 5. We finally conclude in Sect.6.

The scalar singlet extension of the SM with heavy neutrinos
The simplest extension of the SM scalar sector is the inclusion of a real scalar singlet s that may mix with the Higgs boson.This small addition to the SM may however significantly alter the scalar sector phenomenology.In particular, it can allow for a SFOPT even at tree-level [13,15], re-opening the possibility of explaining the origin of the observed matter-antimatter asymmetry of the Universe in the context of EWBG if new sources of CP violation beyond the SM are also present.More interestingly, the scalar singlet field s could be a window to a dark sector capable of addressing some of the other open problems of the SM.Indeed, given its singlet nature, renormalizable (and therefore less suppressed) couplings are expected between the scalar and both the SM and the extended dark sector.Such scenarios could for example account for the observed dark matter of the Universe [23] (see also [24][25][26][27][28][29][30][31]) or the simultaneous origin of neutrino masses and the BAU [18,19].Motivated by the latter, but easily generalizable, in this work we will consider a dark sector that comprises the real scalar field s and n new Dirac neutrinos, N = (N L , N R ), singlets under the SM gauge group and with lepton number +1.In order to have an accessible extended neutrino sector with masses around the EW scale that may play a non-trivial role in the SFOPT and the baryogenesis process, we consider low-scale realizations of the Seesaw mechanism [32][33][34][35] with an approximate lepton number conservation so as to protect and ensure the lightness of neutrino masses [36][37][38].Indeed in the so-called inverse [39,40] or linear [41,42] Seesaw mechanisms the new heavy neutrinos arrange in Dirac pairs while the SM neutrinos remain massless if the lepton number symmetry is exact.Upon softly breaking this symmetry, the SM neutrinos will acquire small masses and the mass degeneracy of the two members of each Dirac pair will be slightly broken [43].Neglecting the small lepton-number-violating terms (which are suppressed by the tiny SM neutrino masses), the most general lepton-number-conserving interaction Lagrangian among the dark sector fields s, N and the SM fields is: where is the SU (2) L doublet Higgs field, ˜ = iσ 2 * , L L is the SU (2) L lepton doublet and Y ν and Y N are general 3 × n and n × n Yukawa matrices, respectively.Without loss of generality we will work in the basis where Y N is diagonal.
The most general Lagrangian scalar potential for the Higgs doublet and the singlet scalar s is given by (see e.g.[15]) By writing 2), the scalar potential for the neutral fields h and s, relevant for EW symmetry breaking, is found to be In the rest of this work, we will denote the (zero-temperature) vacuum expectation values (vevs) of the Higgs and singlet fields stemming from the potential (2.3) by v EW and ω EW , respectively.We note that all parameters from V (h, s) are real, which means that the only sources of CP violation beyond the SM would arise from the Yukawa couplings in Eq. (2.1), in the absence of further new physics contributions.
In general, to study the early Universe dynamics of the scalar sector and the possibility to have a SFOPT, zerotemperature loop corrections (at 1-loop, this corresponds to the so-called Coleman-Weinberg contribution [44]) as well as finite-temperature contributions to the scalar potential [45,46], should be taken into account.These corrections, however, introduce gauge dependence [47] and renormalization scale dependence in the effective potential of the theory, leading to important theoretical uncertainties [48,49]. 1 Nevertheless, the scalar potential of the singlet scalar extension of the SM may already lead to the generation of a tree-level 1 A possibility to alleviate these problems consists of performing dimensional reduction, working with a 3-dimensional effective theory [50][51][52][53][54].This procedure consists in practice on successively integratingout all the heavy energy scales of the system (see e.g.[48] for a recent discussion on the topic).
barrier between the EW symmetric and broken minima [15] and, in such a case, an analysis based on the tree-level potential (2.3) supplemented by the leading (∼ T 2 ) thermal corrections in a high-temperature approximation, which do not depend on the choice of the gauge, captures the most relevant features needed for the study of the SFOPT.At the same time, working at this level of approximation allows to study the relevant features of the phase transition analytically, as advocated in Ref. [15].This is very advantageous in order to efficiently scan the parameter space of the model.We have verified the generic validity of this approximation concerning the results of our global parameter scan, as we discuss in more detail in Sect. 4. The finite-temperature effective potential V T can in this case be written as: (2.4) The explicit appearance as a free parameter in Eq. (2.4) of the critical temperature T c , at which the EW symmetric and broken minima are degenerate in energy, proves very convenient in a scan of the model parameter space requiring the presence of a SFOPT.Indeed, when imposing that at T = T c the two minima are degenerate, an analytical condition among the other potential parameters in Eq. (2.4) is obtained, effectively trading its freedom for T c and allowing to explore only potentials for which the two-degenerate-minima condition is fulfilled (see Appendix A for details).The parameters in Eq. (2.4), defined at T = T c , are related to those of Eq. (2.3), defined at T = 0, by μ2 The constants c h , c s and m 3 are given by where g and g are respectively the SU (2) L and U (1) Y gauge couplings, Y t is the top Yukawa coupling, and For the study of the temperature evolution of the scalar potential minima and the SFOPT, it is also convenient to rewrite the potential V T from Eq. (2.4) in terms of the temperature-dependent vevs v T ≡ h (T ) and ω T ≡ s (T ) in the broken minimum as [15] V where all dimensionful parameters have an implicit dependence on the temperature T .The mass parameters m2 h , m 2 s , and m 2 sh are defined as evaluated at the EW broken minimum at T .The effective coupling λ 2 and mass m * in Eq. (2.6) are defined as Furthermore, in the parameter scans in Sect. 4 we will eventually trade m 2 sh for the quantity ω p , defined as which has the advantage of being temperature-independent.The parametrization (2.6) explicitly shows that a shift in the field s → s + σ keeps the finite-temperature scalar potential V T invariant with a redefinition of ω T → ω T + σ .The relations between the new parameters in Eq. (2.6) and the coefficients in Eq. (2.3) are found in Ref. [15].The starting point of our analysis of SFOPT scenarios is the finite-temperature potential V T from Eq. (2.4) at T = T c , with two degenerate minima located in general at (0, ω 0 ) and (v, ω) ≡ (v T c , ω T c ) in the two-dimensional field space.We will require that both the Higgs and the singlet field acquire a vev in the EW broken phase, in order to generate masses for the heavy neutrinos after the phase transition.Furthermore, in Ref. [19] it was shown that successful baryogenesis in the present scenario favours the heavy neutrinos to be approximately massless at the onset of the SFOPT.We will thus make use of the shift symmetry of the potential via ω T → ω T + σ discussed above to set ω 0 = 0 at the critical temperature T c (this corresponds to setting μ 1 = 0 in Eq. (2.4)), and mainly focus on phase transitions from ( h , s ) = (0, 0) → (v, ω).For studies on other phase transition scenarios in the singlet scalar extension of the SM, we refer the reader to Refs.[17,20,49,[55][56][57][58][59][60][61][62][63].Following [15], we can use Eq.(2.6) to conveniently parametrize the finite-temperature potential at the critical temperature as The set of parameters of Eq. (2.4) at T = T c can be recovered from the new parametrization given in Eq. (2.10) (together with μ 1 = 0) via the following relations 2 : ) ) ) (2.14) ) (2.16) Then, the set of parameters {ω, ω p , v, m 2 h , m 2 s , λ m } evaluated at T = T c is enough to fully characterize the scalar potential with two degenerate minima at (0, 0) and (v, ω). 3 Together with T c , this set of parameters allows us to specify the finitetemperature potential as a function of T in the singlet scalar extension of the SM, in the high-T approximation.When considering the addition of extra heavy Dirac neutrinos, we also need to add Y ν and Y N as free parameters in our analysis, since they appear in the thermal corrections of the potential V T , see Eq. (2.5).   1 illustrates an example of the type of scalar potentials that would lead to a SFOPT with the characteristics described above, with the specific values of the corresponding potential parameter set at T = 0 given in Table 1.In the left panel of Fig. 1, the potential V T (h, s, T ) is shown at T = T c with the two degenerate minima, represented by the purple diamonds.The red dot indicates the location of the saddle point yielding a potential barrier between both minima.In the right panel we show the potential at T = 0 where the EW breaking minimum (v EW , ω E W ) is now the true vacuum and also the real singlet has an O(100 GeV) vev that generates EW-scale masses for the heavy neutrinos.
In general, the potential V T from Eq. (2.10) characterized by a random set of parameters {ω, ω p , v, m 2 h , m 2 s , λ m } at T c satisfying the conditions in Appendix A from Ref. [15], despite satisfying the desired property of featuring two degenerate minima at (0, 0) and (v, ω), will not reproduce the correct value for the EW symmetry breaking vev at T = 0, v E W = 246.22GeV, obtained from the measurement of the Fermi constant via the muon decay width [64].In addition, upon diagonalization of the scalar mass matrix at T = 0, with m 2 h , m 2 s and m 2 sh defined in Eq. (2.7) and evaluated in the T = 0 EW broken minimum (v EW , ω EW ), the eigenvalue M H for the mostly-doublet mass eigenstate will generally not reproduce the measured value for the Higgs boson mass M H = 125.10GeV.Satisfying these two physical requirements at T = 0 is rather non-trivial in our setup, and considerably reduces the allowed parameter space: given the high accuracy of the v EW and M H measurements, two combinations of the free parameters in the scalar potential are effectively determined.In Sect. 4 we will discuss how these requirements are implemented in our numerical scan of the parameter space of the model.
Finally, we also need to consider the existing constraints on the mixing ξ between the Higgs doublet and the scalar singlet, arising from the diagonalisation of the T = 0 scalar mass matrix, M 0 s .We have where S (H ) is the mass eigenstate corresponding to the mostly-singlet (doublet) scalar combination with a mass M S (M H = 125.10GeV).In the next section we will discuss the present experimental constraints on the (T = 0) model parameters, affecting in particular the possible allowed values of the mixing ξ , which will also be applied to our parameter scan in Sect. 4.

Experimental constraints and phenomenological probes
In this section we discuss the relevant experimental limits on the singlet scalar extension of the SM, making emphasis on how the possible presence of the extra heavy singlet neutrinos can affect them.These experimental constraints will translate into bounds on the parameters of the potential from Eq. ( 2.4) at T = 0. We also discuss the main phenomenological probes of the model, particularly in connection with both the structure of the scalar potential and the presence of the heavy neutrinos in comparison with the minimal singlet scalar extension of the SM.

SM-heavy neutrino mixing
The new Dirac neutrinos with components N R and N L introduced in Eq. (2.1) mix with the SM neutrinos after spontaneous symmetry breaking (SSB) and may participate in the generation of light neutrino masses.The mixing matrix between the SM active and the heavy sterile neutrinos is given by with the Yukawa matrices Y N and Y ν for singlet and SM neutrinos (see Eq. (2.1) for details), respectively.Thus, neglecting all small L-violating parameters that would eventually lead to the generation of the masses of the mostly SMlike light neutrinos, the heavy neutrinos have Dirac masses M N i ω EW Y Ni and the two chiralities of the mass eigenstates N i are given by For sterile neutrinos with masses M N i > M W , the activesterile neutrino mixing is bounded from above by a combination of EW precision tests 4 and flavour observables [67,68] tr θθ † < 0.0048 (2σ ).
(3.3) 4 Notice that the recent anomalous measurement of M W by the CDF II collaboration [65] could be potentially explained through a non-zero neutrino mixing [66].However, this result is in tension with the other observables and we conservatively do not take it into account here.
For lighter sterile neutrinos (M N i < M W ), direct searches at colliders and beam dump experiments as well as searches for peaks and distortions in the decay products of mesons, leptons and beta decays set much more stringent constraints on the active-sterile neutrino mixing, and we refer the reader to Refs.[69,70] for a comprehensive list of these limits.
The bound (3.3) will be used as a χ 2 contribution added to the weight function that we construct to guide our scan of the parameter space (see Appendix A).To constrain since tr (AB) ≤ tr (A) tr (B).Thus, when imposing the present bound on tr θθ † to the right-hand side of Eq (3.4), a conservative bound is implemented since tr θθ † will always be smaller than this quantity.In practice, because we only have access to Y N and not the individual values of the Yukawas, when constraining Y ν we assume a degenerate spectrum such that tr Any other choice would translate into smaller values for Y ν .

LHC Higgs signal strengths
The latest measurements of the 125 GeV Higgs boson signal strength by the ATLAS [71] and CMS [72] collaborations provide an important constraint on deviations of Higgs couplings from their SM values.In the singlet scalar extension of the SM, all SM couplings to the Higgs-like mass eigenstate H become rescaled relative to the SM values by cos ξ , which, in the absence of exotic Higgs decays (see discussion below), yields an overall suppression of Higgs signal strength given by This allows to constrain the singlet-doublet scalar mixing via Higgs measurements (see [73][74][75][76]).We use the latest measurements of the Higgs signal strength from ATLAS, μ = 1.05 ± 0.06 [77] and CMS, μ = 1.002 ± 0.057 [78], and combine them to derive a bound on cos ξ following the Feldman-Cousins [79] prescription. 5We find μ ≥ 0.94 at the 95% C.L. which translates into |sin ξ | < 0.245, as shown in Fig. 2.This is in fairly good agreement with other recent analyses (see e.g.[76]) where Feldman-Cousins is however not applied in general.
In the presence of an exotic Higgs branching fraction into beyond-the-SM (BSM) states BR H →BSM , the Higgs signal strength modifier μ in Eq. (3.5) becomes BR H +X = 0 (dashed) and BR H +X = 0.9 (dash-dotted), and from LEP searches for light scalars (blue) respectively assuming BR S→BSM ≡ BR X = 0 (dashed), BR X = 0.9 with visible BSM decays (dash-dotted) and BR X = 0.9 with invisible BSM decays (dotted) since cos ξ and BR H →BSM now yield a combined dilution of the global Higgs signal strength with respect to the SM.Thus, the presence of exotic Higgs decays yield a tighter bound on cos 2 ξ from Higgs signal strengths, as shown in Fig. 2 for the specific value BR H →BSM ≡ BR X = 0.04 6 .
In particular, the interactions in Eq. (2.1) between the Higgs boson H and the extra heavy neutrino states This occurs via Y N in Eq. (2.1), through the singlet-doublet scalar mixing, or via Y ν , through the active-sterile neutrino mixing.In addition, the interaction Y ν may also mediate H → ν Ni , ν N i decays [80,81].Since the singlet-doublet scalar mixing is much more weakly constrained than the active-sterile neutrino mixing, the leading interaction (assuming θ 2 sin 2 ξ ) after EW symmetry breaking between the scalar states and the N i (i = 1, . . ., n), which is induced by the Lagrangian from Eq. (2.1), would be )), and we concentrate in the following on the leading interaction from Eq. (3.7).The decay H → N i Ni can have a significant impact on the LHC bounds on the Higgs signal strength.From the latest μ measurements performed by ATLAS [77] and CMS [78] and discussed above, we set the bound 71 The total rate of the Higgs-like boson with a mass of 125 GeV decaying into N i states is given by where k ≤ n is the number of kinematically accessible heavy neutrinos.The maximum possible value of H →N N (occur- which is achieved for Y 2 N i = M 2 H /(10 ω 2 EW ) with i = 1, 2, . . .n, and may well be comparable to the SM Higgs boson total width SM = 0.00412 GeV [82] if the mixing sin ξ is not too suppressed.For a given value of , which is the relevant combination of neutrino Yukawa couplings affecting the thermal history of the scalar sector, the following two "Cases" are possible: 2 : Then all neutrinos are kinematically accessible and k = n in Eq. (3.9) so that the 125 GeV Higgs boson decays into all n heavy neutrinos.For a fixed Y 2 N the maximum value of H →N N is achieved when all Yukawa couplings are equal (Y 2 while the minimum value is obtained when the rate is dominated by a single heavy neutrino contribution (Y 2 The decay of the 125 GeV Higgs boson into at least one heavy neutrino may be kinemat-ically forbidden. 8Therefore, for a given value of Y 2  N , the decay rate can be arbitrarily suppressed depending on the value of the individual Yukawas Y N i (e.g. in the limit Y and no lower bound on H →N N exists.Still, sizable Higgs boson branching ratios into sterile neutrinos are also possible (even reaching H →N N = max ) in this case for a fixed value of Y 2 N .
The inclusion of the exotic H → N N decay channel on the Higgs signal strength bound (recall Eq. (3.8)) allows to exclude a significant fraction of the parameter space in which a SFOPT is possible in the present scenario, as we will show explicitly in Sect. 5. Finally, we stress that the heavy neutrinos produced in the decays of the 125 GeV Higgs bosons may themselves decay visibly inside the detector via active-heavy neutrino mixing, leaving a prompt or displaced vertex signal in the detector depending on the value of θ 2 .In particular, if the heavy neutrinos N i are long-lived, they can lead to a two-displaced-vertices signal in the LHC detectors, which would be a very powerful probe of the model [84,85].

Electroweak precision observables
The properties of the singlet field are also constrained by EW precision observables (EWPO), which limit the value of the mixing sin ξ as a function of the scalar mass M S in the singlet scalar extension of the SM (see e.g.[86,87]).This is a result of the shift induced by the presence of the singlet scalar on the EW oblique parameters S, T , U [88] with respect to the SM.A global fit to EWPO measurements yields the respective values of the shifts on the oblique parameters with respect to their SM predictions [89] S = 0.04 ± 0.11, T = 0.09 ± 0.14, with the following correlation coefficients: +0.92 between S and T , −0.68 between S and U and −0.87 between T and U .Explicit expressions for S, T , U in the singlet scalar extension of the SM are given in [87] as a function of sin ξ and M S .Using these, we obtain the 95% C.L. limits on the (M S , |sin ξ |) plane from a χ 2 fit to the S, T , U measurements from Eq. (3.11).These are shown in Fig. 2, highlighting that for values of M S below a TeV, the bound from Higgs signal strength discussed in the previous section is stronger than that of EWPO.Notice, however, that the same EWPO used to constrain S, T and U and, from there, derive constraints on sin ξ are 8 Writing , for R/n > 1 the Higgs decay into at least one heavy neutrino N i must be kinematically forbidden, while for R/n < 1 it is still possible to have k = n in Eq. (3.9).
affected already at tree level and used to derive the bounds on the heavy-active neutrino mixing θ as outlined above [67].In principle, the two contributions should be studied together to derive a consistent set of constraints.The interplay between new physics contributions to the EWPO through S, T , U and the presence of heavy neutrinos was studied in detail in Refs.[90][91][92][93].In particular, it was realized that most observables depend on the same combination of elements of θ and T and that, if a cancellation between these two contributions is present, the bounds on both sources of new physics would weaken significantly.Nevertheless, for this situation to take place, negative and sizable values of T are required [93].The scalar singlet contribution to T does indeed become negative for masses above the mass of the Higgs (see e.g.[87]).For lighter singlet masses, no cancellation is possible and the two effects would rather reinforce each other, leading to slightly stronger constraints.Nevertheless, since the bounds from Higgs signal strength are more stringent, the potential contribution of the singlet is small and does not alter significantly the constraints on heavy-active neutrino mixing derived in [93].Conversely, for a scalar singlet heavier than the Higgs, the bound tr θθ † ≤ 0.0048 would weaken if −2αT ∼ 0.0048.However, given the bounds on sin ξ from the LHC Higgs signal strength measurements (see Sect. 3.2), this is never achieved for sub-Planckian scalar masses.Thus, for the parameter space under study, the possible interplay between the heavy neutrino and scalar singlet contributions to EWPO can be safely neglected.

Searches for singlet-like scalars at LEP and LHC
Under the assumption that the singlet-like scalar decays into SM particles (i.e. its decay is driven by the singlet-doublet mixing), the null results from LEP searches for Higgs-like particles yield strong upper limits on |sin ξ | for singlet-like scalar masses below M S 115 GeV (see e.g.[73]).These limits are at the level of |sin ξ | 0.2 (or below) for masses M S < 100 GeV.At the same time, LHC searches for BSM scalars decaying to W W , Z Z or H H pairs also constrain the doublet admixture of the singlet-like scalar S for M S > M H .For M S > 200 GeV the strongest such limits are obtained by ATLAS in the Z Z → 4 and Z Z → 2 2ν final states [83]. 9In Fig. 2 we show the corresponding bounds on the (M S , |sin ξ |) plane from both LEP and LHC searches for new scalars, under the assumption BR S→H H +BR S→BSM ≡ BR H +X = 0 (dashed lines).
Nevertheless, compared to the minimal singlet extension of the SM, here the presence of the heavy neutrinos may lead to much less stringent bounds on sin ξ from direct scalar searches.Indeed, the interactions of S with the heavy neutrinos N i in Eq. (3.7) will induce the decay S → N i Ni if available by phase space.The corresponding partial width S N N ∝ cos 2 ξ , in contrast to the partial decay widths of S into SM states, S SM ∝ sin 2 ξ .Thus, S → N i Ni will generally be the dominant decay channel for the singlet-like scalar in the limit |sin ξ | 1 (as favoured by LHC Higgs signal strength measurements, see Sect.3.2): • For light singlets (M S 100 GeV), the S → N i Ni decay channel would significantly relax constraints on sin ξ from LEP searches for Higgs bosons decaying visibly (into SM particles), and we show the corresponding dilution of the limits when BR X = BR S→N i Ni = 0.9 in Fig. 2. Nonetheless, if the heavy neutrinos N i are longlived (e.g. for very small neutrino mixing) and would have escaped the LEP detectors, limits from LEP searches for invisibly decaying Higgses [96][97][98][99] would apply.We also depict the bounds from such searches on sin ξ in Fig. 2, showing that they become very strong for rather light scalars.We nevertheless re-stress that these only apply under specific conditions (very long-lived N i , leading to invisible S decays), which depend on the details of the neutrino sector of the model.
• For M S 200 GeV, the presence of the S → N i Ni decay would weaken the LHC limits on sin ξ from pp → S → Z Z searches, as shown explicitly in Fig. 2 for BR H +X = 0.9. 10 At the same time, this BSM decay would open a new avenue to probe the existence of S and N i at the LHC, either when the N i decay products are resolved in the ATLAS/CMS detector or merge into a single reconstructed object (for M S M N i , producing a "neutrino jet" [103,104]).Yet, current LHC searches for heavy neutrinos generally consider N i production modes (e.g.Drell-Yan or W γ fusion, see [105] for a discussion) which yield kinematic properties of the N i rather different from those of S → N i Ni decays, 11 and as such present LHC limits (see [107,108] for reviews) are difficult to extrapolate to our scenario.Moreover, the possibility that the N i yield displaced decays (for θ 2  1) would dramatically reduce the sensitivity of those existing searches, providing at the same time a new avenue for 10 The presence of a non-zero S → H H partial width S H H would also weaken the limits on sin ξ from pp → S → Z Z searches, allowing at the same time to search for S via resonant di-Higgs production (see e.g.[100,101]).Yet, di-Higgs searches are generally less sensitive than Z Z ones for equal branching fractions, and the equivalence theorem [102] naively yields S H H ∼ S Z Z in the M S v limit (since also S H H ∝ sin 2 ξ ).We have thus not considered here the would-be limits from resonant di-Higgs searches for S H H = 0 for simplicity. 11An exception is given by LHC searches for Z gauge bosons decaying to heavy neutrinos [106], which feature similar kinematics and could allow for a reinterpretation in our setup.We defer this for future work.
discovery in long-lived particle searches to be explored in the future.

Higgs self-coupling
Finally, the existence of the singlet scalar would induce a deviation on the Higgs boson trilinear self-coupling λ H H H from its SM value.This is currently being searched for at the LHC [109][110][111] via non-resonant di-Higgs production, albeit with limited precision at present.At tree-level, we find with cos ξ (sin ξ) = c ξ (s ξ ).Additionally, and particularly relevant in the |sin ξ | 1 limit, the one-loop corrections to the trilinear self-coupling coupling should be taken into account, as they contain terms that do not vanish even when sin ξ → 0. The one-loop contribution reads, in the |sin ξ | 1 limit [17] ( We parametrize the deviation with respect to the SM as with λ SM H H H = M 2 H /(2v EW ) the tree-level value of the SM Higgs boson self-coupling. 12The most stringent constraint on κ λ has been recently set by the ATLAS Collaboration [111] to be − 1.0 ≤ κ λ ≤ 6.6 (95% C.L.). (3.15) We note that the measurements of λ H H H at the LHC via nonresonant di-Higgs production can be significantly altered by the presence of a resonant contribution to the di-Higgs signature (see [113,114] for a discussion), appearing in the singlet scalar extension of the SM via the pp → S → H H process.Still, depending on the singlet-like scalar mass M S , it should be possible to exploit the di-Higgs invariant mass distribution m H H to measure the self-coupling λ H H H [113] from the non-resonant part of the distribution (with the resonant part properly identified and subtracted), achieving comparable precision to the scenario with no resonant S → H H contribution.A detailed analysis of this issue is however beyond the scope of the present work.Moreover, as we will see in Sect. 5 this observable barely deviates from its SM value in the interesting regions of the parameter space and is therefore not a relevant probe of the scenario under study.

Parameter scan
In this section we describe our procedure to scan the parameter space of the model and collect the sets of parameters which fulfill the various necessary conditions for a SFOPT.As described in Sect.2, in general a potential characterized by a random set of parameters {ω, Following Ref. [15], we bias our scan towards the parameter sets that satisfy the necessary conditions for a SFOPT. 13 For this purpose, we have defined an ad-hoc weight function to rate how well the selected points satisfy these conditions, in order to prioritize the parameter regions to which the points belong in our scan.We then use this weight function in place of the log-likelihood for a Markov Chain MonteCarlo (MCMC) using the standard Metropolis Hastings algorithm to sample the interesting regions of the parameter space with MonteCUBES [115].The conditions for a SFOPT and the weight function used in the MCMC are explicitly defined in Appendix A. Our procedure of the parameter scan is summarized in Fig. 3.
Finally, we have verified the validity of the high-T approximation adopted in Sect. 2 through a comparison of the value of T c obtained in this approximation with the one obtained from a potential with the same parameters at T = 0 but implementing the temperature evolution with the full 1-loop ther- 13 Among the different conditions, we look for potentials that are bounded from below for which the EW minimum is the global one.While the EW minimum could be metastable, such setups are beyond the scope of this work.

Keep point in parameter space
No Yes Rejected Accepted Fig. 3 Flowchart for the selection of parameter sets generating the correct Higgs vev and mass at T = 0 and satisfying all necessary conditions to potentially have a SFOPT from Ref. [15] mal potential (see e.g.[116]).The comparison has been performed for points passing all viability criteria, as discussed below.Both values of T c agree within 5% accuracy for the vast majority of the points scanned (a few outliers extend to ∼ 10%).Then, even though in some cases the value for the scalar singlet mass M S found in the scan is of the same order as T c , the high-T approximation can remain suitable for our purposes (a fast efficient scan of the parameter space).In addition, we impose a stringent perturbativity condition on the scalar quartic couplings, λ i ≤ 2 (see below), which favours that radiative corrections are not significant.
The parameter sets output of our MCMC scan are further classified according to the following viability criteria: • Points with the scalar potential quartic couplings λ i ≤ 2 to ensure perturbativity.• Points that lead to a sufficiently strong first-order phase transition (if the phase transition occurs).As a rough estimate, we ask for the ratio v/T c > 1.This is required to make EW baryogenesis possible by decoupling sphaleron processes in the EW broken phase.• Points for which the bubbles of the EW broken phase can actually nucleate and the phase transition does take place.
Although the conditions summarized in Fig. 3 are needed to realize a SFOPT, they are not sufficient to guarantee it.
It is important to study whether a nucleation temperature T N < T c exists for which the bubbles of the EW broken phase (the true vacuum for T < T c ) successfully grow [17,60,62,117,118] (and the Universe does not become trapped in the false vacuum).In our scan of parameters, we gauge the nucleation of EW bubbles as follows: The transition probability from the false to the true vacuum is proportional to e −S 3 /T , with S 3 the threedimensional bounce action.At temperatures slightly below T c , in the so-called thinwall regime for which the two minima are almost degenerate, the action S 3 diverges for T → T c [119,120] and thus no transition is possible in this regime.As the Universe cools down from T c to the nucleation temperature T N , away from the thin-wall limit, the computation of the bounce action becomes more involved.It requires to solve the equations of motion for the background fields to find the bounce solution.This is usually done via numerical solvers like CosmoTransitions [121], BubbleProfiler [122] or FindBounce [123].Nucleation will be possible if there exists a temperature T N > 0 at which S 3 /T N ∼ 140 [124][125][126][127], for which the nucleation rate is comparable to the Hubble expansion rate during radiation domination.
In order to estimate the bounce action and thus the nucleation temperature T N , instead of computing the bounce solution along the path that minimizes the tunneling action, we approximate the solution by calculating the bounce action along a straight path in field space, which connects both minima at T N .The action for such a field configuration, S app 3 , will by construction be larger than the tunneling solution [128][129][130], S 3 ≤ S app 3 .We find that there is good agreement between the true action S 3 and our estimation for the cases of interest, and thus successful nucleation for points in parameter space is expected to occur when S app 3 /T N 140, which in turn represents a conservative estimate.In practice, we rewrite the scalar potential in terms of a linear combination of h and s along the straight line connecting the minima at a given T , φ , and the orthogonal one, φ ⊥ , as prescribed in Ref. [131].By taking φ ⊥ → 0 one can quickly find the bounce solution along the straight line using the overshoot-undershoot method in one dimension with, for example, FindBounce.In this manner, the action S app 3 can be computed at the temperature T .The temperature at which S app 3 /T ∼ 140 defines the nucleation temperature T N .As an example of our approximation, we compare in Fig. 4 the approximated and the actual bounce trajectories in field space with the set of parameters given in Table 1.This approximation allows us to efficiently scan the parameter space.In the results shown in the follow- Fig. 4 The same as Fig. 1 but at the nucleation temperature T N .The blue line corresponds to the straight path in the field configuration space which we use as an approximation to estimate the action and judge if nucleation may happen, while the purple curve corresponds to the actual bounce solution passing closer to the saddle point (red dot) between the two minima (purple diamonds) ing section, the boundaries between the nucleating and non-nucleating points should therefore be understood as an approximated result and somewhat conservative.In Sect. 5 we will perform a comparison between the regions of parameter space selected by our nucleation criteria and those found using CosmoTransitions, finding good agreement between both methods.
• Points with the scalar mixing ξ at T = 0 allowed by collider searches as described in the previous section.
In the next section, we will show the impact of each of these conditions in the parameter space to reveal the correlations among the parameters and the preference for particular parameter regions.As we will see, the condition of the bubble nucleation will prove to be the most constraining one [17,60,62,118], which greatly reduces the allowed parameter space.

Results
In this section we present and analyze the results from the parameter scan described in Sect. 4 where the different constraints and conditions described in the previous section have been implemented.While the scan is performed over all parameters (at h fixed so as to reproduce the correct Higgs mass, the constraint in the activeheavy neutrino mixing tr θθ † | exp ≤ 0.0048 [67] implies Y 2 ν will have a negligible impact on the scalar potential.We will therefore not show this parameter in the follow-   5 Results of the parameter scan in the scalar mass, M S , and mixing, sin ξ , plane.In the right (left) panel the scan was performed with(out) the addition of the heavy Dirac neutrinos.The pink stars correspond to points with non-perturbative couplings, the purple crosses to points for which sphaleron transitions would not decouple in the broken phase.The light blue dots (dark blue pluses) additionally do (not) satisfy the nucleation condition as described in the text.The grey-shaded region corresponds to the bounds on the scalar mixing described in Sect.3. In the right panel, we do not show the bounds that depend on the parameters in each point ing.Instead we will mainly focus on parameters which have a direct connection to experimental observables, and thus refer the reader to Figs. 9 and 10 in Sect.6 for results on the complete set of parameters at T c , w, m 2 h , and T = 0 respectively, as a summary of our results.
In Fig. 5 we show the points collected in our parameter scan in the plane of the mass and mixing of the scalar singlet (M S , |sinξ |).We study and compare the scenarios with(out) the addition of the heavy Dirac neutrinos in the right (left) panels.All the points have been selected according to the algorithm summarized in Fig. 3 and, therefore, satisfy the conditions from Ref. [15] for a SFOPT and have the correct Higgs mass and vev at T = 0.The points with different colours and symbols are classified by the conditions listed in Sect. 4. The pink stars are discarded since they have at least one very large scalar coupling (λ i > 2). 14For the purple crosses, this perturbativity condition is satisfied, but the firstorder phase transition is not strong enough to decouple the sphaleron process in the EW broken phase (v/T c < 1), even if the bubbles of the broken phase may nucleate.The dark blue crosses labeled with "No T N " have λ i < 2 and v/T c > 1 but the nucleation condition S app 3 /T 140 is not satisfied at any T < T c (and therefore, there is no T N ).Finally, the light blue dots labeled with "T N " have λ i < 2 and v/T c > 1 and also fulfill the nucleation condition.Grey-shaded areas in the left panel represent the values of the scalar mixing ruled out by LHC Higgs signal strength measurements (assuming BR X = 0) as described in Sect.3.2, or by direct searches 14 Such large couplings will drive the model into a non-perturbative regime at scales very close to the EW scale, and thus we disregard those points in our scan.
for Higgs-like particles at LEP for M S 100 GeV and at ATLAS for M S > 200 GeV (assuming BR S→H H = 0).On the right panel we only display the conservative Higgs signal strength bound in the absence of exotic Higgs decays, since the bounds from direct scalar searches at LEP and LHC may be diluted when heavy Dirac neutrinos are included, depending on the values of neutrino couplings as discussed in Sect.3.4.As can be seen from the plots, these constraints are quite relevant and a big portion of the parameter space is ruled out by them, so that only small values of sin ξ are still allowed.Moreover, we also find that the condition of successful bubble nucleation considerably reduces the size of the viable parameter space, as pointed out in Refs.[117,118] for other scenarios.Thus, only the light blue dots below the grey-shaded regions are successful candidates for a SFOPT satisfying all phenomenological constraints listed in the previous section.From Fig. 5 we can also see that the Universe may undergo an EW SFOPT only if the mass of the singlet scalar S is M S 300 GeV.Generally speaking, higher values of M S would also imply larger ω and hence a significant distance between the two minima, in general too large to allow for bubble nucleation.The apparent exception to this rule by the few points clustered around M S ∼ 500−1000 GeV can be understood through a closer inspection of their thermal evolution.Indeed, in these cases we find a SFOPT only in the singlet direction at T O(100) GeV.After this transition, both the Higgs and the singlet vevs roll towards their values at T = 0, v EW and ω EW , respectively, with ω E W v EW .Through the comparison between the left and right panels of Fig. 5, we can study the impact of the presence of heavy Dirac neutrinos.While we find new nucleating regions char-Fig.6 Correlations between different parameters in the scan with nonzero neutrino Yukawa coupling.The allowed regions are very similar for the scan without neutrino Yukawas (see also Fig. 10).The color coding is the same as in Fig. 5 with the grey diamonds indicating the parame-ter space which can undergo successful nucleation but are excluded by their value of the scalar mixing ξ .In all panels the parameters are represented at T = T c except for the bottom-right panel with the trilinear Higgs coupling at T = 0 acterized by large ω and sizeable and negative λ m (as can be seen in Fig. 10), these are largely ruled out by the constraints from Higgs signal strength measurements, as shown in Fig. 9. Indeed, as shown in the upper-left panel of Fig. 6, the points that pass the criteria for the case including the heavy neutrinos cluster at small values of Y N or small ω.The condition that leads to the (hyperbola-like) correlation shown in this panel is the requirement of the stability of the EW broken minimum imposed in Eq. (A.3), in particular that m 2 s > 0 at T = 0. Since Y N induces a negative evolution of m 2 s from T c down to T = 0, values of Y N > 1 are constrained in the scan unless T c is small and/or λ m is negative to cancel its effect in Eq. (2.5).This can be seen explicitly in the corresponding panels of Fig. 9.Such a cancellation is however prevented by the bounds on the scalar mixing, as indicated by the grey points in Fig. 9. Thus, from now on, we will only present the results of our scan with Y N = 0, i.e. in the presence of heavy Dirac neutrinos.However, the allowed regions should also be considered generally valid for the Y N = 0 scenario without the heavy neutrinos, with the caveat that direct scalar searches at LEP and ATLAS further constrain the parameter space.
We present the distribution of scan points in various other interesting slices of the parameter space in Fig. 6.For information on the distributions under the other parameters of the scan, we refer the reader to Fig. 9.The color coding and symbols is the same as for Fig. 5, but in addition to them, we now have the grey diamonds (labeled with "T N + cos ξ ") indicating the parameter space points which can undergo successful nucleation but are excluded by their value of the scalar mixing sin ξ .They correspond to the blue dots covered by the grey-shaded area in the right panel of Fig. 5.In the upper right panel of Fig. 6 we show the correlations found in our scan between the two scalar vevs at T = T c .As can be seen, the most significant constraint is the requirement of a sufficiently strong EW phase transition, v/T c > 1.When imposing this together with v EW = v exp EW and M H = M exp H , values of v below ∼ 100 GeV are ruled out.Besides this constraint, we find that large values of ω, beyond ∼ 200 GeV, are disfavored by the requirement of successful nucleation.Indeed, we generally find that, if ω were too large, the field distance between the two minima would be too big to allow for bubble nucleation despite satisfying the rest of the requirements.Even though we find regions of parameter space successfully nucleating for singlet vevs as large as ω ∼ 1000 GeV at T c , a detailed study of these regions shows that these transitions at T N occur from (0, ω − ω) → (0, ω), with ω ω, such as the distance travelled in field space is not qualitatively larger than for the region with ω 200 GeV.
In the lower left panel of Fig. 6 we show the distribution of scan points in the ω −λ m plane at T = T c .We find that the light-blue dots, which pass all requirements and in particular the nucleation condition, display an anticorrelation between these two parameters.Additionally, the bounds on sin ξ rule out most of the points with λ m < 0 unless |λ m | 1.These trends can be understood from the hyperbolic shape of the correlation between λ m and ω p , which is found in the corresponding panel of Fig. 9. Indeed, from Eq. (2.9) this behavior is expected if μ m is negative.Analyzing the accepted samples, we find that negative μ m is preferred in order to satisfy our condition for successful nucleation.In fact, negative μ m decreases the barrier between the two degenerate minima and thus we find no nucleating samples for positive μ m .Finally, the area with negative λ m < 0 and negative ω p is ruled out by the constraints on the scalar singlet mixing since, as expected from Eq. (2.9), −m 2 sh would become too large.We have further analyzed the non-trivial correlation found between ω and λ m when imposing our criteria for nucleation, comparing these results with points that successfully nucleate according to the tunneling module from CosmoTransitions.As can be seen in Fig. 7, for a subset of our sample featuring successful nucleation, the areas found by both our approximate estimate (light-blue points) and CosmoTransitions (black octagons for a first-order EW phase transition) generally agree well.The two exceptions we identify are: (i) the region with negative λ m and large ω, where CosmoTransitions finds successfully tunneling points which are not found by our approximation. 15This region corresponds to significantly more curved trajectories 15 For positive λ m and large ω, a few of the light-blue points yielding successful nucleation with our criteria are instead tagged as second-order phase transitions (2 nd OPT, green pentagons) by CosmoTransitions While we have not explicitly discerned the order of the transition in our scan (which is beyond our present scope), we note that no qualitative new regions appear when considering such parameter points, as this region falls within the areas where a SFOPT is found by our method.than those depicted in Fig. 4, not-well approximated by our straight-line assumption.Nevertheless, this whole area of the parameter space also leads to too large scalar singletdoublet mixing (as described above) and is experimentally excluded.(ii) The points for which CosmoTransitions does not find an EW phase transition, whereas our nucleation proxy does, i.e. the red squares and the blue points with no counterpart (neither red square, green pentagon or black octagon) in Fig. 7.This should a priori never happen, since our criterium for nucleation is conservative.A careful investigation of such points reveals that the phasetracking module of CosmoTransitions does not produce numerically reliable results in such cases. 16We thus conclude that neither of these exceptions is meaningful, and our estimate for nucleation agrees well with the results from CosmoTransitions for the values of the parameters for which CosmoTransitions yields a reliable numerical result, thus representing an efficient and fast alternative for scans of the parameter space discriminating in a conservative way if nucleation could happen.
As discussed in Sect.3, the inclusion of the singlet scalar causes a deviation of the Higgs trilinear coupling from its SM value, which can be parametrized as in Eq. (3.14).The distributions of κ λ as a function of ω E W at T = 0 is shown 16 Even if by construction two phases are always present for our model parameter points, in these cases CosmoTransitions fails to find one of them for the default numerical precision in the code.A significant increase in the numerical precision generally leads to CosmoTransitions finding the second phase and identifying a first-order transition, in agreement with our estimate.Nevertheless, this increase in numerical precision makes the computation too slow to allow for an efficient scan of the parameter space.  in the bottom-right panel of Fig. 6, together with the current bounds from collider searches, which are given by Eq. (3.15) and shown by the grey-shaded regions.We find that the lightblue points satisfying all conditions tend to cluster in a narrow range around κ λ ∼ 1.Thus, given the sensitivity to κ λ ∼ 2.2 [132] of future probes such as the HL-LHC, no deviations caused by a singlet scalar responsible for a SFOPT are to be expected in this observable.Finally, in Fig. 8 we show the regions of the parameter space which are constrained by the possible new decay channel of the Higgs-like state H into heavy neutrinos using Eqs.(3.8) and (3.9) as described in Sect.3.2, for points which pass all the constraints (i.e.light-blue in Fig. 6) in our parameter scan.In the left panel we show the contribution to the corresponding branching ratio assuming a degenerate heavy neutrino spectrum (Y 2 for which H →N N can be comparable to or even exceed the SM Higgs boson total width SM in an important part of the parameter space.The solid red line separates "Case 1" and "Case 2" as discussed in Sect.3.2.Note that in the region above the solid red line ("Case 2"), a different combination of Y N i Yukawa couplings (yielding the same value of Y 2 N ) could arbitrarily reduce the value of BR H →N N by making all neutrinos either too heavy for the Higgs to decay into or with negligible couplings.From this panel we can also infer that the heavy Dirac neutrinos are in general lighter than ∼ 300 GeV.In the right panel of Fig. 8 we instead show the minimum possible value of BR H →N N for each parameter point.Notice that for the points corresponding to "Case 1" (region below the solid red line in the left panel) the exclusion limits from Higgs signal strength measurements, shown in grey, are unavoidable and rule out a significant region of the parameter space, while for "Case 2" the BR H →N N can be made arbitrarily small and thus the bound can always be evaded.In Fig. 9 we show in red the points excluded by B R H →N N in the different relevant planes in parameter space.Even though these constraints are important, as seen in Fig. 8, they do not exclude particular regions of parameter space.

Summary and conclusions
In this work we have explored the parameter space of the scalar singlet extension of the SM with the aim of identifying the regions in which a SFOPT, as required to explain the puzzle to the origin of the observed baryon asymmetry through the EWBG mechanism, can take place.The main goal of the study is to contribute to the predictability of the scenario by relating the areas where a SFOPT can happen with potentially testable observables or correlations among them.
Previous studies [15] showed the conditions that need to be met by the extended scalar potential in order to develop two degenerate minima at some critical temperature T c .Together with the requirement of reproducing the correct mass and vacuum expectation value of the SM-like Higgs and of the required strength of the transition (v/T c > 1), these set of constraints already impose stringent and non-trivial conditions of the allowed parameter space.
Nevertheless, as advocated by [17,60,62,117,118], we find that the requirement that bubble nucleation may actually take place between the two minima is the most constraining requirement, reducing drastically the allowed parameter space.Furthermore, testing explicitly this condition is not possible in a fast and analytical way and relying on the numerical solvers available [121][122][123] necessarily limits the speed of the scan hindering the exploration of large parameter spaces.Moreover, given the complexity of the problem, for some points in the parameter space we find that some numerical solvers fail to find one of the phases, and hence the corresponding transition, or are unable to produce a result.For this reason we have adopted a fast and conservative approximation to the three-dimensional action of the bounce solution, S app 3 , that controls the transition rate between the two minima and requires S app 3 /T N ∼ 140 at some nucleation temperature T N .We find that for most of the sampled points in the parameter space that satisfy this criteria, CosmoTransitions does indeed find a first order phase transition (with a small fraction of second order transitions, something we did not explicitly discriminate), thus validating our approach.We also point out that, for many of the points that passed our selection criteria, CosmoTransitions did not provide an output.Thus, larger regions of the parameter space may be explored in a fast and efficient way through the approximation adopted, although it should be taken as a conservative estimate and not as an exact result.
In our scan of the parameter space we find that the regions with the correct mass and vev for the Higgs and successful nucleation are mainly characterized by values of the singlet vev ω EW 300 GeV.Indeed, if ω EW is too large, the two minima tend to be too far apart in field space 17 and nucleation may not happen.This in turn translates into values of the scalar singlet mass that cluster around M S 300 GeV.The exception to this rule is a clustering of allowed points with large ω and values of M S in the 500-1000 GeV range, which in any case do not produce an EW phase transition given that v(T ) smoothly goes from 0 to v E W as the Universe expands.We have verified that for these points the actual jump in ω during the phase transition is also small.Regarding the most constraining observables, we find that the bounds on the singlet-doublet mixing from Higgs signal strength measurements by ATLAS and CMS are already ruling out important regions of the parameter space.Direct searches for the singlet scalar when its decays are SM Higgs-like both at LEP and at LHC are also relevant.
We have also investigated how this picture is affected when the scalar singlet is not alone, but part of larger dark sector it may interact with.As a particularly motivated scenario, we considered as case study the addition of extra sterile neutrino singlets of both chiralities.These new states will have Yukawa couplings Y N to the scalar singlet, which would induce Dirac masses around the EW scale for these heavy neutral leptons.Furthermore, a Yukawa coupling Y ν among the SM Higgs doublet, the SM neutrinos and the heavy neutrinos would generally also be allowed.The simultaneous presence of Y N and Y ν implies a new source of CP-violation that may be enough to induce the BAU via EWBG [18,19] (the so-called ν-EWBG scenario [19]).Furthermore, if a small source of lepton-number violation is introduced, the presence of Y ν would induce small neutrino masses able to explain the neutrino oscillation phenomenon in the manner of the lowscale symmetry-protected seesaws like the inverse or linear seesaw variants.
Previous studies [21] showed that the presence of the heavy neutrinos increases the strength of the transition by enhancing v/T c .We reproduce this result, but find that sizable Y N , unless accompanied by small ω, can also destabilize the broken minimum.Thus, when our criteria for nucleation and stability are imposed, for the allowed values of the Yukawa couplings the regions of the scalar potential parameter space are comparable to the scenario without heavy neutrinos.Hence the early universe phenomenology regarding the possibility of a SFOPT of both scenarios is very similar, as summarized in Figs. 9 and 10.Even if new areas appear when including the neutrinos, we observe in Fig. 9 that they are excluded by Higgs signal strength measurements.The small values of the Yukawa couplings and the scalar singlet vev ω preferred, seem to make the generation of the BAU via EWBG difficult according to the findings of [19], but a dedicated analysis would be required to confirm or rule out its viability.
Conversely, the presence of the heavy sterile neutrinos may significantly affect the collider phenomenology of the scalar singlet extension.Indeed, while very large values of Y N could hinder vacuum stability, values around Y N ∼ O(1) are perfectly allowed.Such a sizable coupling would on the one hand imply that the scalar singlet decays would be overwhelmingly dominated to the heavy sterile neutrino channel, given that the singlet-doublet mixing is more strongly constrained and the smaller SM Yukawa couplings.This would in turn invalidate the bounds on the scalar mixing derived from direct searches of the singlet with SMlike decays at LEP and LHC.On the other hand, dedicated searches for this new decay channel should be considered.
Furthermore, if allowed by phase space, the decay to heavy sterile neutrinos could also be sizable for the SM-like Higgs scalar via its mixing with the singlet.We have found that this in fact strengthens the Higgs signal strength constraints in significant portions of the parameter space, correspond-  9 Triangle plot for non-zero Yukawas, only for points successfully nucleating.The grey diamonds correspond to regions of parameter space excluded by constraints on the scalar mixing while red dots are excluded by the constraint from H → N N .The blue points satisfy all phenomenological bounds.We note the very strong correlation between the singlet vev ω and ω p at T c for blue points, as well as between λ m and ω p ing to the red points in Fig. 8. Interestingly, the possibility that this is the dominant channel to produce and test for the heavy neutrinos at collider searches also remains open in parts of the parameter space.Indeed, the mixing of the heavy neutrinos with their active counterparts induced by Y ν is more strongly constrained from flavour and electroweak precision observables, as well as collider searches via Drell-Yan production.Thus, if Y N is more sizable than Y ν , the heavy neutrinos would be more easily produced via Higgs or singlet decays.For small enough values of Y ν , the decays of the heavy neutrinos themselves would not be prompt and may induce interesting signatures with displaced vertexes.We thus find that the viable parameter space allows for very striking and non-standard collider phenomenology Data Availability Statement This manuscript has no associated data or the data will not be deposited.[Authors' comment: As this is a theoretical study, no data were collected for our study.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.Funded by SCOAP 3 .SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.

A Parameter scan with weight function
We describe in this section the ad-hoc weight function used for the parameter scan and also show the distribution of the parameter points that satisfy the nucleation condition in the full parameter space at Fig. 9.
A set of necessary conditions for a successful SFOPT were identified, analytically, in Ref. [15], and can be found in Table 1 therein.They fall into the following three categories: • Conditions to have degenerate minima at T = T c : In order for the two minima to be stable we impose with m 2 h | 0 ≡ (ω/2)[λ m − m 2 h /ω 2 − 2m 2 sh /(vω)].These expressions are valid under the assumption that the symmetric minimum sits at (0, 0), which we can assume without loss of generality thanks to the shift symmetry present in the potential.
• Matching conditions at T c : Once a viable degenerate potential is found at T c with the reduced set of parameters {ω, ω p , v, m 2 h , m 2 s , λ m }, we need to set the parameters λ 2 , m * (given by Eq. (2.8)) and m 2 sh to particular values so as to have the general potential from Eq. (2.6).In this step we ensure that the running of the potential according to temperature change makes the broken minimum the deepest one by imposing We require that the potential is bounded from below, and the broken minimum is the global minimum and stable, which are translated into the following conditions on the parameters at T = 0.
Notice that, as seems to be the case in the SM, the electroweak minimum could be metastable at T = 0 with a lifetime longer than the age of the Universe.Thus, this condition is more restrictive than strictly necessary.Nevertheless, it is convenient since it allows to efficiently scan the potential without the need of computing the lifetime of the vacuum.
All the conditions listed above take the form "C X > 0" with C X a given function of the parameters in the potential.In order to guide our scan towards the areas where these conditions are met we construct the following weight function with W X defined as The factor 10 6 is a penalty to the points that do not satisfy one of the conditions, with which we can make sure any point accepted in the MCMC satisfies all the necessary conditions for the SFOPT.We also add a Gaussian χ 2 term for the Higgs mass as well as for the constraint on Y ν from Eq. (3.= 0.17 GeV is its uncertainty [64] and tr θθ † exp = 0.0014 with σ tr(θθ † ) = 0.0014 [67,68], and we take n = 3.

8 Fig. 1
Fig.1Shape of the scalar potential for the example point in parameter space specified in Table1at the critical temperature T c (left panel) and at T = 0 (right panel).The color bar denotes the value of the scalar potential, V T .The purple diamonds denote the location of the minima,

Fig. 2
Fig.2Existing 95% C.L. constraints on the singlet-doublet scalar mixing sin ξ as a function of M S from EWPO (solid green), from LHC Higgs signal strength measurements (red) with BR H →BSM ≡ BR X = 0 (solid) and BR X = 0.04 (dash-dotted), from direct searches for pp → S → Z Z by ATLAS[83] (ocher), with BR S→H H + BR S→BSM ≡ BR H +X = 0 (dashed) and BR H +X = 0.9 (dash-dotted), and from LEP searches for light scalars (blue) respectively assuming BR S→BSM ≡ BR X = 0 (dashed), BR X = 0.9 with visible BSM decays (dash-dotted) and BR X = 0.9 with invisible BSM decays (dotted) will not yield the correct values of the Higgs vev and Higgs mass at T = 0, namely v E W = v exp EW ≡ 246.22 GeV and M H = M exp H ≡ 125.10 GeV.The first condition v E W = 246.22GeV can always be imposed starting from any given set of parameters by shifting all the parameters η with dimension of mass (including T c ) as η → (v exp E W /v EW ) η.To satisfy the second condition, M H = M exp H , we solve for m 2 h for each generated set w = {ω, ω p , v, m 2 s , λ m , T c , Y ν , Y N } to find values which reproduce the correct Higgs boson mass: in practice, for a given set w we scan m 2 h imposing v E W /M H = v exp EW /M exp H before the aforementioned shift of the mass dimensionful parameters η, which guarantees M H = M exp H after it.A solution does not always exist depending on the actual values of w.In this way, we obtain the sets of parameters which have two degenerate minima at T = T c and also reproduce the correct Higgs vev and mass at T = 0.

Fig. 7
Fig. 7 Sample comparison between the results obtained with our proxy (light blue circles) and CosmoTransitions (CT) in the ω − λ m plane for successfully nucleating points.According to CT, red squares correspond to non-nucleating points, green pentagons to those giving a 2 nd order phase transition (OPT) and black octagons to successfully nucleating points

Fig. 8
Fig. 8 Bounds on the regions of parameter space giving rise to successful nucleation.The left panel correspond to the plane Y N against ω EW at T = 0, with the color legend giving the size of BR H →N N for n = 3 degenerate heavy neutrinos.The solid red line corresponding to ω EW Y N = (M H /2) separates the cases where the decay is always allowed (below) and where it depends strongly on the assumptions on Fig. 9 Triangle plot for non-zero Yukawas, only for points successfully nucleating.The grey diamonds correspond to regions of parameter space excluded by constraints on the scalar mixing while red dots are excluded by the constraint from H → N N .The blue points satisfy all

Fig. 10
Fig. 10 Triangle plot comparing the successfully nucleating points for the singlet scalar alone (orange pentagons) and including the neutral fermions (blue dots), without imposing any phenomenological bound, d(V T,b − V T,s ) dT 2 T c = c h v 2 + ω(c s ω + 2m 3 ) > 0, (A.2)where V T,b(s) correspond to the potential in the broken (symmetric) minimum.Notice that in this step the neutrino Yukawa couplings Y ν(N ) may play an important role as they enter into the c h and c s constants which set the running of the potential with the temperature.• Conditions for potential at T ≤ T c :

Table 1
Parameter set at T = 0 .7)The H → N i Ni decay channel is then driven by the scalar mixing, H →N i Ni ∝ sin 2 ξ , while in minimal seesaw scenarios the heavy neutrinos are produced via mixing with the SM neutrinos, leading to H →N i Ni ∝ θ 4 and H →ν N i ∝ θ 2 .In our scenario these generally correspond to subleading effects (the corrections for both H and S interactions in Eq. (3.7) are O(θ 2