The Higgs Trilinear Coupling and the Scale of New Physics for the SM-Axion-Seesaw-Higgs Portal Inﬂation (SMASH) Model

: In the extended scalar sector of the SMASH (Standard Model - Axion-Seesaw-Higgs portal inﬂation) framework, we conduct a phenomenological investigation of the observable effects. In a suitable region of the SMASH scalar parameter spaces, we solve the vacuum metastability problem and discuss the one-loop correction to the triple Higgs coupling, λ HHH . The λ HHH and SM Higgs quartic coupling λ H corrections are found to be proportional to the threshold correction. A large λ HHH correction ( (cid:38) 5%) implies vacuum instability in the model and thus limits the general class of theories that use threshold correction. We performed a full two-loop renormalization group analysis of the SMASH model. The SMASH framework has also been used to estimate the evolution of lepton asymmetry in the universe.


Introduction
After the discovery of the Standard Model (SM) Higgs boson [1,2], every elementary particle of the SM has been confirmed to exist.Even though the past forty years have been a spectacular triumph for the SM, the mass of the Higgs boson (m H = 125.25 ± 0.17 GeV) [3] poses a serious problem for the SM [4].It is well-known that the SM Higgs potential is metastable [5], as the sign of the quartic coupling, λ H , turns negative at instability scale Λ IS ∼ 10 11 GeV.On the other hand, the SM is devoid of non-perturbative problems since the non-perturbative scale Λ NS M Pl , where M Pl = 1.22 × 10 19 GeV is the Planck scale, but still there are studies on non-perturbative effects of the SM [6][7][8][9][10].In the post-Planckian regime, effects of quantum gravity are expected to dominate, and the non-perturbative scale is therefore well beyond the validity region of the SM, unlike the instability scale.The largest uncertainties in SM vacuum stability are driven by top quark pole mass and the mass of the SM Higgs boson [11].The current data are in significant tension with the stability hypothesis, making it more likely that the universe is in a false vacuum state [12][13][14][15].The expected lifetime of vacuum decay to a true vacuum is extraordinarily long, and it is unlikely to affect the evolution of the universe [16,17].However, it is unclear why the vacuum state entered into a false vacuum to begin with during the early universe.In this post-SM era, the emergence of vacuum stability problems (among many others) forces the particle theorists to expand the SM in such a way that the λ H will stay positive during the run all the way up to the Planck scale.
It is possible that at or below the instability scale, heavy degrees of freedom originating from a theory beyond the SM start to alter the running of the SM parameters of renormalization group equations (RGE).It has been shown that incorporating the Type-I seesaw mechanism [18][19][20][21][22][23][24][25][26][27][28] will have a large destabilizing effect if the neutrino Yukawa couplings are large [29], and an insignificantly small effect if they are small.Thus, to solve the vacuum stability problem simultaneously with neutrino mass, a larger theory extension is required.Embedding the invisible axion model [30][31][32] together with the Type-I seesaw was considered in [33,34].The axion appears as a phase of a complex singlet scalar field.This approach aims to solve the vacuum stability problem by proving that the universe is currently in a true vacuum.The scalar sector of such a theory may stabilize the vacuum with a threshold mechanism [35,36].The effective SM Higgs coupling gains a positive correction δ ≡ λ 2  Hσ /λ σ at m ρ , where λ Hσ is the Higgs doublet-singlet portal coupling and λ σ is the quartic coupling of the new scalar.
Corrections altering λ H in such a model would also induce corrections to the triple Higgs coupling, λ tree HHH = 3m 2 H /v, where v = 246.22 GeV is the SM Higgs vacuum expectation value (VEV) [37][38][39].The triple Higgs coupling is uniquely determined by the SM but is unmeasured.In fact, the Run 2 data from the Large Hadron Collider (LHC) has only been able to determine the upper limit of the coupling to be 12 times the SM prediction [3].Therefore, future prospects of measuring a deviation of triple Higgs coupling by the high-luminosity upgrade of the LHC (HL-LHC) [40,41] or by a planned next-generation Future Circular Collider (FCC) [42][43][44][45][46][47][48][49] give us hints of the structure of the scalar sector of a beyond-the-SM theory.Previous work has shown that large corrections to triple Higgs coupling might originate from a theory with one extra Dirac neutrino [50,51], inverse seesaw model [52], two Higgs doublet model [38,39,53,54], one extra scalar singlet [37,55,56] or in the Type II seesaw model [57].
The complex singlet scalar, and consequently, the corresponding threshold mechanism, is embedded in a recent SMASH [58][59][60] theory, which utilizes it at λ Hσ ∼ −10 −6 and λ σ ∼ 10 −10 .The mechanism turns out to be dominant unless the new Yukawa couplings of SMASH are O(1).In addition to its simple scalar sector extension, SMASH includes electroweak singlet quarks Q and Q and three heavy right-handed Majorana neutrinos N 1 , N 2 and N 3 to generate masses for neutrinos.
The structure of this paper is as follows: In Section 2, we summarize the SMASH model and cover the relevant details of its scalar sector.We also establish the connection between the threshold correction and the leading order λ HHH correction.In Section 3, we discuss the methods, numerical details, RGE running, and our choice of benchmark points.Our results are presented in Section 4, where the viable parameter space is constrained by various current experimental limits.In SMASH, one can obtain at most ∼5% correction to λ HHH while simultaneously stabilizing the vacuum.We give our short conclusions on Section 5.

Theory
The SMASH framework [58][59][60] expands the scalar sector of the SM by introducing a complex singlet field where ρ and A (the axion) are real scalar fields, and v σ v is the VEV of the complex singlet.The scalar potential of SMASH is then Defining φ 1 = H and φ 2 = σ, in basis (H, σ), the scalar mass matrix of this potential is which has eigenvalues and At the heavy singlet limit and Defining threshold correction δ ≡ λ 2 Hσ /λ σ in Equation ( 13), and The first term in the Equation ( 9) is the leading component.
The SMASH framework also includes a new quark-like field, Q, which has color but is an electro-weak singlet.It gains its mass via the Higgs mechanism, through a complex singlet σ.It arises from the Yukawa term We will show later that Y Q = O(1) is forbidden by the vacuum stability requirement.The hypercharge of Q is chosen to be q = −1/3, even though q = 2/3 is possible.Our analysis is almost independent of the hypercharge assignment.Threshold correction: Consider an energy scale below m ρ < Λ IS , where the heavy scalar ρ is integrated out.The low-energy Higgs potential should match the SM Higgs potential It turns out that the quartic coupling we measure has an additional term Since the SM Higgs quartic coupling will be approximately λ H (M Pl ) ≈ −0.02, the threshold correction Hσ λ σ (13) should have a minimum value close to |λ H (M Pl )| or slightly larger to push the highenergy counterpart λ H to positive value all the way up to M Pl .A too-large correction will, however, increase λ H too rapidly, exceeding the perturbativity limit √ 4π.We demonstrate the conditions for δ in Section 4. Similar to λ H , the SM Higgs quadratic parameter µ H gains a threshold correction In the literature [35,36], there are two possible ways of implementing this threshold mechanism.One may start by solving the SM RGE's up to m ρ , where the new singlet effects kick in, and the quadratic and quartic couplings gain sudden increments.Continuation of RGE analysis to even higher scales then requires utilizing the new RGE's up to the Planck scale.
Another approach is to only solve the new RGEs on the SM scale while ignoring the low-scale SM RGEs entirely.We will use the former approach.
One-loop correction to triple Higgs coupling: The portal term of the Higgs potential contains the trilinear couplings for HHρ and Hρρ vertices.The vertex factors for HHρ and Hρρ vertices are introduced in Figure 1.The one-loop diagrams contributing to SM triple Higgs coupling are in Figure 2. We denote the SM tree-level triple Higgs coupling as λ HHH .The correction is gained by adding all the triangle diagrams (taking into account the symmetry factors) Here, p and q are the external momenta and the loop integral is defined as The contribution from diagram Figure 3 is subleading, since it is proportional to They can be derived from Equation (2).We denote ρ and its propagator by red color.The process H * → HH is disallowed for on-shell external momenta, so at least one of them must be off-shell.Specifically, the momentum-dependent correction to the triple coupling at the tree-level is an effective coupling that enters the specific process with one off-shell Higgs decaying into two real Higgses.Note that the correction is dependent on the Higgs off-shell momentum q ≡ q * , which we assume to be at O(1) TeV at the LHC and HL-LHC.The first diagram is dominant due to the heaviness of the ρ scalar.Therefore, we may ignore the subleading contributions of diagrams involving two or more ρ propagators.We integrate out the heavy scalar, causing the finite integral in Equation ( 16) to be logarithmically divergent.We calculate the finite part of it using dimensional regularization and obtain where z ≡ 1 + (4m 2 H /q 2 ) and µ = m ρ is the regularization scale 1 .We have used the modified minimal subtraction scheme (MS), where the terms ln 4π and Euler-Mascheroni constant γ E ≈ 0.57722 emerging in the calculation are absorbed to the regularization scale µ.For calculations, we use the value q * = 1 TeV.It is especially interesting to see that at the leading order, the triple Higgs coupling correction is proportional to the threshold corrections.This intimate connection forbids a too-large correction.In fact, the bound from vacuum stability turns out to constrain the triple Higgs coupling correction to 5%, as we shall see in Section 4. Consequently, if LHC or HL-LHC manages to measure a correction to λ HHH , this will rule out theories that utilize exclusively threshold correction mechanisms as a viable solution to the vacuum stability problem.Indeed, there are alternate ways to produce large ∆λ HHH without expanding the scalar sector [50,52].
It should be noted that loop corrections contributing to the final to-be-observed value are included in the SM.Indeed, experiments are measuring λ SM HHH = λ SM(tree) , where the SM one-loop correction depends on the Higgs off-shell momentum.At the O(1) TeV scale we are considering, the SM 1-loop correction amounts to approximately −7% [50].
Light neutrino masses: The neutrino sector of SMASH is able to generate correct neutrino masses and observe the baryon asymmetry of the universe with suitable benchmarks.The relevant Yukawa terms for neutrinos in the model are We take a simplified approach: Dirac and Majorana Yukawa matrices (Y ν and Y n , respectively) are assumed to be diagonal.
To generate baryonic asymmetry in the universe, SMASH utilizes the thermal leptogenesis scenario [61], which generates lepton asymmetry in the early universe and leads to baryon asymmetry.In the scenario, heavy neutrinos require a sufficient mass hierarchy [62,63] and one or more Yukawa couplings must have complex CP phase factors.We assume the CP phases are O(1) radians to near-maximize the CP asymmetry [64][65][66] If the CP violation is maximal, the largest value is obtained.To produce matter-antimatter asymmetry in the universe, a large asymmetry is required.Following [58], we set the heavy neutrino mass hierarchy at These choices give the full 6 × 6 neutrino mass matrix which is in block form, and contains two free parameters: √ 2 is the Majorana mass term.Light neutrino masses are then generated via the well-known Type I seesaw mechanism [18][19][20][21][22][23][24][25][26][27][28], by block diagonalizing the full neutrino mass matrix M ν .
The light neutrino mass matrix is After removing the irrelevant sign via field redefinition, where we have denoted C = v 2 /( √ 2v σ ) and assumed normal mass ordering m 1 < m 2 < m 3 .This gives the neutrino masses m i = Cy 2 i /Y i .We do not know the absolute masses, but the mass squared differences have been measured by various neutrino oscillation experiments [67,79].Nevertheless, their values provide two constraints, leaving three free parameters.However, the heavy neutrino Yukawa couplings Y i must be no larger than O(10 −3 ) to avoid vacuum instability [59].
In addition, an order-of-magnitude estimate of the generated matter-antimatter asymmetry (baryon-to-photon ratio) is directly proportional to the CP asymmetry where κ ∼ 0.01-0.1 is an efficiency factor.We arrive at which, in principle, can be consistent with the observed η.To achieve successful resonance leptogenesis, v σ should be between 10 10 and 10 12 GeV (Table 1).We will provide suitable benchmark points in the next section.The estimation of lepton asymmetry, which is one of the crucial implications of SMASH as the framework claims to solve the matter-asymmetry issue since the scenario only consists of the decay and inverse decay of N 2 or N 3 to N 1 .The leptogenesis evolution for the benchmark values shown in Table 2 is in Figure 4.   12) and Figure 1 of [74] with mass bound from [68])  The evolution of the "abundance" of N 2 in blue, the "abundance" of N 2 in thermal equilibrium in red, and the lepton asymmetry generated by the CP violating decays and inverse decays of N 2 divided by the CP asymmetry parameter ε CP in black.T is the temperature of the universe in GeV, as well as the mass of heavy neutrinos M 2 in GeV.∆L is the number of changes in lepton number over entropy density s.
We will investigate the influence of N 1 , N 2 , and N 3 oscillations (i.e., right-handed neutrino oscillations) on leptogenesis evolutions, predict baryon-to-photon ratios for different set masses of light active left-handed neutrinos, and evaluate a more precise value of κ by solving complicated Boltzmann equations in the future course of analysis in the SMASH framework.

Methods
We generate the suitable benchmark points demonstrating different physics aspects of the model in the neutrino sector by fitting in the known neutrino mass squared differences ∆m 2 ij , assuming normal mass ordering (m 1 < m 2 < m 3 ).This leaves three free neutrino parameters, the values of which we generate by logarithmically distributed random sampling.These are the candidates for benchmark points.We then require that the candidate points be consistent with the bounds for the sum of light neutrino masses [68][69][70][71][72][73][74][75][76][77][78].The next step is to choose suitable values for other unknown parameters, using the stability of the vacuum as a requirement.
The authors of [58] have generated the corrections to the two-loop β functions of SMASH.We solve numerically the full two-loop 14 coupled renormalization group differential equations with SMASH corrections with respect to Yukawa (Y u , Y d , Y e , Y ν , Y n , Y Q ), gauge (g 1 , g 2 , g 3 ) and scalar couplings (µ 2 H , µ 2 S , λ H , λ σ , λ Hσ ), ignoring the light SM degrees of freedom, from M Z to Planck scale.We assume Yukawa matrices are on a diagonal basis, with the exception of Y ν .We use the MS scheme for the running of the RGE's.Since the top quark MS mass is different from its pole mass, the difference is taken into account via the relation [80] where α 3 ≡ g 2 3 /4π ≈ 0.1085 at µ = m Z .We define the Higgs quadratic coupling as µ H = m H / √ 2 and quartic coupling as λ H = m 2 H /2v 2 .We use MATLAB R2019's ode45-solver.See Table 1 for the used SMASH benchmark points and Table 3 for our SM input [3].Our scale convention is t ≡ log 10 µ/ GeV.In some papers, the running of SM parameters (Y t , Y b , Y τ , g 1 , g 2 , g 3 , µ 2 H , λ H ) obeys the SM RGE's without corrections from a more effective theory until some intermediate scale Λ BSM [35], after which the SM parameters gain threshold correction (where it is relevant) and the running of all SM parameters follows the new RGE's from that point onwards.We choose to utilize this approach while acknowledging an alternative approach, where the threshold correction is applied at the beginning (µ = m Z ) [36], and both approaches give almost the same results.As previously stated, SM Higgs quadratic and quartic couplings will gain the threshold correction.
Our aim is to find suitable benchmark points, which: • Allow the quartic and Yukawa couplings of the theory to remain positive and perturbative up to the Planck scale; • Utilize a threshold correction mechanism to λ H via δ 0.1; • Avoid the overproduction of dark radiation via the cosmic axion background (requiring λ Hσ < 0); • Produce a significant contribution matter-antimatter asymmetry via leptogenesis (requiring hierarchy between the heavy neutrinos); • Produce a ∼5% correction to triple Higgs coupling λ HHH .

Stability of vacuum:
We have plotted how the running of the SM quartic coupling, λ H changes with each benchmark point in Figure 5.Note that all the threshold corrections are utilized well before the SM instability scale Λ IS .One can choose v σ > Λ IS if m ρ < Λ IS is ensured.This is the case with BP3.Our results for the chosen benchmarks are in Figure 6, where the SM best fit is denoted by a red star.Clearly, the electroweak vacuum is stable with our benchmark points, and it is assigned to m pole t 172.69 ± 0.3 GeV and m H 125.25 ± 0.17 GeV [3].For every case, we investigated the running of the quartic couplings of the scalar potential.We used the following stability conditions and for λ Hσ < 0 [35] −λ Hσ (µ) If one or more conditions are not met on the scale µ ∈ [m Z , M Pl ], we denote this point as unstable.If any of the quartic couplings rises above √ 4π, we denote this point non-perturbative.) plane with benchmark points BP1-BP3.The red star corresponds to the SM best-fit value.The height and width of the star correspond to the present uncertainties.The vacuum is stable in the yellow region.The contour numbers n correspond to the vacuum instability scale 10 n GeV.
We chose the new scalar parameters in such a way that the threshold correction is large but allowed, 0.1 < δ < λ H .This changes the behavior of the coupling's running so that after the correction, the λ H increases in energy instead of decreasing, the opposite of the coupling's running in a pure SM scenario.A too-large threshold correction will have an undesired effect, lowering the non-perturbative scale to energies lower than the Planck scale.These effects are visualized in Figure 7, where for each benchmark point kept λ σ at its designated value in Table 1.Instead, we let the portal coupling, λ Hσ , vary between 0 and √ 0.6λ σ .This demonstrates the small range of viable parameters space.We have also investigated the significance of v σ on the bounds of threshold correction δ.A choice of δ is available as long as v σ 2 × 10 13 GeV.This can be seen clearly from Figure 8.Given a fixed δ, the result is independent of λ Hσ and λ σ .The lower and higher bound for δ increases as a function of v σ .Instability bound increases, since the needed vacuum-stabilizing threshold effect increases as one approaches the SM instability scale Λ IS .At v σ 2 × 10 13 GeV, the m ρ > Λ IS , so the quartic coupling λ H will turn negative before threshold correction is utilized.On the other hand, the non-perturbative scale increases, since, as the cutoff point m ρ increases, the quartic coupling λ H decreases and correspondingly, the largest possible threshold correction increases.The available parameter space is consistent with a stable vacuum in (v σ , δ) space.λ σ is fixed, while λ Hσ is determined by Equation ( 13) and m ρ by Equation ( 9).We have denoted our benchmark points with a red star.
Our next scan was over the new quartic couplings, log 10 (−λ Hσ ) ∈ [−7, 0] and log 10 λ σ ∈ [−10, 0].The scalar potential is stable and the couplings remain perturbative at only a narrow band, where δ ∼ 0.01-0.1,see Figure 9.If one considers small δ, the SM Higgs quartic coupling will decrease to near zero at µ = M Pl .This corresponds to a region near the left side of the stability band.In contrast, we chose our benchmarks with large δ, placing it near the right side of the stability band, corresponding to the large value of λ H at µ = M Pl .This was a deliberate choice to maximize the correction to λ HHH .In addition, we have scanned the Dirac neutrino and new quark-like particle Yukawa couplings (y 1 and Y Q , respectively) over y 1 ∈ [0, 2] and Y Q ∈ [0, 0.04], keeping y 2 and y 3 small, real 2 and positive but non-zero.See Figure 10 for details corresponding to each benchmark point.There, we point to an area producing a stable vacuum.The Dirac neutrino Yukawa couplings may have a maximum value of O( 1), but a more stringent constraint is found for Y Q .It should be noted that even though, from the vacuum instability point of view, Y max Q < y max 1 , this does not imply Y Q < y 1 , since both are in principle free parameters.See Table 2 for computed values for neutrino masses for normal hierarchy (m 1 < m 2 < m 3 ) corresponding to each benchmark.Note that all BP1-BP3 produce a value of baryon-to-photon ratio comparable to experimental values and a mass of axion consistent with axion dark matter scenario, because it requires axion decay constant f A ≡ v σ to be O (10 11 ) GeV [30][31][32].In Figure 4, we show the evolution of N 2 or N 3 abundance, as well as the lepton asymmetry generated by the CP violating decays and inverse decays of N 2 or N 3 , divided by the CP asymmetry parameter ε CP as per [62].The resulting lepton asymmetry is translated to baryon asymmetry via the sphaleron process with a c s fraction.We have also shown the N 2 or N 3 abundance in thermal equilibrium.The number density n of particles decreases in an expanding universe if there are now particle number-changing interactions.However, the ratio of number density n to entropy density s, that is, "abundance" = n/s, is constant.Changing "abundance" during the early universe thus indicates particle interactions, or in our case, N 2 or N 3 decays and inverse decays.A corresponding mass hierarchy for right-handed neutrinos implies an upper bound of ε CP ∼ 10 −5 to 10 −6 [63,81].
Correction to SM triple Higgs coupling: According to PDG [3], the largest possible experimental value for λ HHH is 12 times the SM prediction 3 , from Run 2 data for the b bγγ channel alone.The real singlet scalar ρ mixes with the SM Higgs, providing a one-loop correction to SM triple Higgs coupling λ HHH .We scanned the parameter space with log 10 (−λ Hσ ) ∈ [−7, 0] and log 10 λ σ ∈ [−10, 0].At each point, we calculated the correction to λ HHH .See Figure 11 for details.We identified a section of parameter space excluded by triple Higgs coupling searches from LHC run 2 and determined the area sensitive to future experiments, namely, HL-LHC and FCC-hh.We assume HL-LHC uses 14 TeV center-ofmass energy and integrated luminosity L = 3 ab −1 , for FCC-hh we assume center-of-mass energy 100 TeV and integrated luminosity L = 3 ab −1 .The relative correction in Table 4 is calculated with respect to the SM tree-level prediction.We chose our benchmark points in a way that their correction to triple Higgs coupling will be borderline observable at FCC-hh [82], that is, the correction will be ∼5%.Thus, η in BP3 for a factor of 10 larger is necessary for stable vacuum and FCC-hh better detection shown in Figure 11.Future FCC-hh accelerator is sensitive to ∼5% deviation of the Standard Model prediction.This is demonstrated by the benchmark points we have chosen.Although the model's stable region allows for even smaller deviations, part of the region is still accessible by FCC-hh.This has implications for a general class of BSM theories that utilize complex singlet scalars and other new non-scalar fields.If the corrections from non-scalar contributions to SM triple Higgs and quartic couplings are tiny, any large correction to λ HHH (such as a discrepancy from a SM value measured by the HL-LHC) would rule out such a class of theories, including SMASH.It will be up to the HL-LHC experiment to determine whether this is the case.

Conclusions
We have investigated suitable benchmark scenarios for the simplest SMASH model regarding the scalars and neutrinos, constraining the new Yukawa couplings and scalar couplings via the vacuum stability and theory perturbativity requirements.The model can easily account for the neutrino sector, predicting the correct light neutrino mass spectrum while evading the experimental bounds for right-handed heavy sterile Majorana neutrinos.In [58], the authors of the SMASH model performed a one-loop RGE analysis of the model and presented the two-loop RGE's.We have extended the analysis to two-loop to gain the increased precision needed for the combined achievement of a stabilized electroweak vacuum and a large-enough triple Higgs coupling correction to be sensitive at FCC-hh.To the best of the authors' knowledge, this is the first report on the connection between threshold correction to λ H and one-loop correction to λ HHH .
We found an interesting interplay between the triple Higgs coupling correction and the SM Higgs quartic coupling correction.A successful vacuum stabilization mechanism (threshold mechanism) in SMASH is consistent with small triple Higgs coupling corrections, requiring it to be at most ∼5%.Since the ∆λ HHH is proportional to the threshold correction δ, a large correction to ∆λ HHH inevitably leads to a large threshold correction.Detecting a λ HHH correction larger than ∼35% is within the sensitivity of a future high-luminosity upgrade of the LHC [40,41].If detected, it would, therefore, rule out the simplest scalar sector of the model completely.This would force the model to develop non-minimal alternatives, such as an additional scalar doublet or triplet instead of a singlet.These alternatives were considered by the authors of the SMASH model in their recently updated study [60].The lepton asymmetry |∆L/ε CP | is around 6 × 10 −7 to 10 −5 for the present-day scenario of the universe, which can be verified experimentally [83] at the FCC [47], LHC [84], and by the Circular Electron Positron Collider (CEPC) [85] for the SMASH framework.

Figure 1 .
Figure 1.Vertex factors on trilinear vertices involving the SM Higgs boson as well as a real singlet ρ.They can be derived from Equation(2).We denote ρ and its propagator by red color.

Figure 2 .
Figure 2. One-loop corrections to SM triple Higgs coupling induced by the existence of an extra scalar singlet.In Equation (15), the correction ∆λ HHH is derived.

Figure 3 .
Figure 3. One-loop SM triple Higgs coupling correction diagram with a cubic vertex and a quartic vertex.

Figure 4 .
Figure 4.The evolution of the "abundance" of N 2 in blue, the "abundance" of N 2 in thermal equilibrium in red, and the lepton asymmetry generated by the CP violating decays and inverse decays of N 2 divided by the CP asymmetry parameter ε CP in black.T is the temperature of the universe in GeV, as well as the mass of heavy neutrinos M 2 in GeV.∆L is the number of changes in lepton number over entropy density s.

Figure 5 .
Figure 5.Running of SM Higgs quartic coupling in Standard Model (dashed line) and in SMASH with benchmark points BP1-BP3 (solid line).Threshold correction is utilized at m ρ .

Figure 6 .
Figure 6.Vacuum stability of SMASH in (m H , m pole t

Figure 8 .
Figure 8.The available parameter space is consistent with a stable vacuum in (v σ , δ) space.λ σ is fixed, while λ Hσ is determined by Equation (13) and m ρ by Equation(9).We have denoted our benchmark points with a red star.

Figure 9 .
Figure 9. (Above): Different regions in the logarithmic (−λ Hσ , λ σ ) plane.The contour numbers n above the yellow band correspond to the vacuum instability scale 10 n GeV.Below the yellow band, the contour numbers m correspond to the non-perturbative scale 10 m GeV.The color coding is interpreted as in Figure 6.For non-perturbative scale calculations, we used BP1.(Below): Zoomed-in detail of the figure above, showing in addition our chosen benchmarks.

Figure 11 .
Figure 11.(Above): Different regions in the logarithmic (−λ Hσ , λ σ ) plane.The yellow band corresponds to a stable vacuum configuration.The red area is excluded from the second run of the Large Hadron Collider since the triple Higgs coupling corrections to SMASH would be too large.The dashed line corresponds to the expected sensitivity of the high-luminosity LHC, and the dotted line to the expected sensitivity of the Future Circular Collider in hadronic collision mode.(Below): Zoomed-in detail of the figure above, showing in addition our chosen benchmarks.

Table 1 .
Used benchmark points (BP) in our analysis.Note that we assume specific texture to right-handed neutrino Yukawa matrix Y n .

Table 3 .
[3]d SM inputs in our analysis, at µ = m Z = 91.18GeV,with the exception of top mass, which is evaluated at µ = m t .Masses and vacuum expectation values are in GeV units[3].