We consider an extension of the Standard Model (SM) with additional gauge singlets which exhibits a strong first-order phase transition. Due to this first-order phase transition in the early universe gravitational waves are produced. We estimate the contributions such as the sound wave, the bubble wall collision, and the plasma turbulence to the stochastic gravitational wave background, and we find that the strength at the peak frequency is large enough to be detected at future gravitational interferometers such as eLISA. Deviations in the various Higgs boson self-couplings are also evaluated.
1. Introduction
The discovery of a narrow resonance, with a mass near 125 GeV, at the Large Hadron Collider (LHC) with properties similar to those of the Higgs boson predicted by the Standard Model (SM) [1, 2], sparked a lot of excitements among high energy physicists. But in spite of this important discovery the SM is considered to be incomplete. For instance, in the Standard Model of particle physics a strong first-order phase transition (SFPT) does not occur [3]. However SFPT is needed to justify the baryon asymmetry of our universe [4]; moreover, the SM does not have a candidate for dark matter (DM).
Therefore, some new models are required to address these issues. A popular model is to couple a singlet scalar to Higgs boson. In [5] it is shown that it is possible to modify the standard theory by adding a scalar which possesses a discrete Z2 symmetry and to address the issue of the dark matter of the universe, within the frame work of singlet extended SM; the issue of dark matter has been studied in [6–13], while electroweak phase transition was studied in [14–23] and in [24, 25]; the authors attempt to explain electroweak phase transition EWPT and dark matter by singlet extended SM.
Another class of models are the multisinglet extensions of the SM [26–36]. These models have a larger parameter space in comparison to the singlet extended models; hence they can address several issues; in [26, 27] cosmological implications of such models with classical conformal invariance are presented. Electroweak phase transitions in two-Higgs doublet model are analyzed in [37, 38] and within supersymmetric models in [39–45]. A comprehensive review of EWPT within various models has been given in [46].
In order to investigate the dynamics of the electroweak phase transition EWPT one has to utilize techniques from the domain of thermal field theory [47–50]. The occurrence of a first-order phase transition requires that the electroweak breaking and preserving minima be degenerate, an event which happens at a critical temperature Tc. Moreover, to prevent the washout of any baryon asymmetry by electroweak sphalerons, the electroweak phase transition must be strongly first order. Namely, the ratio of vacuum expectation value of the Higgs field to the critical temperature needs to be greater than unity. As we describe in the next section if a first-order phase transition occurs in the early universe, the dynamics of bubble collision and subsequent turbulence of the plasma are expected to generate gravitational wave (GW). If we detect these GW, then we can obtain information about symmetry breaking in early universe. GW signals from phase transitions have been discussed in [51–59].
In [60] the authors study EWPT within several exotic models and in [61] an analysis of the EWPT of a large number of minimal extensions of the SM and their classically conformal limits is presented. Complex conformal singlet extension of SM with emphasis on the issue of dark matter and Higgs phenomenology is studied in [62]. Recently baryogenesis within a φ6 model is addressed in [63]. An investigation of electroweak phase transitions in a singlet extended model in the 100 (TeV) range is given in [64]. The authors of [65] study strong first-order EWPT in a singlet scalar extension of the SM where the singlet scalar is coupled nonminimally to gravity. In this scheme the singlet field first derives inflation and at a later time causes a strong EWPT; in a new study a first-order EWPT in the SM is obtained by varying Yukawas during phase transition [66].
In this work we propose a new model and we investigate the strength of EWPT within this model. In this model which is a generalization of [19], N real gauge singlets are coupled to Higgs boson via trilinear interactions. Previous studies of multiscalar singlet extension of SM impose separate Z2 symmetries on the singlets [26, 28, 29, 31, 33, 35, 36]; however in our model we do not require such symmetry. In spite of its simple form, this model has a very rich phenomenology. The main feature of the singlet extended SM (without Z2 symmetry) is that the potential barrier between the true and the false vacua necessary for first-order EWPT can be formed mainly by tree-level interactions. But in the singlet extended SM (with Z2 symmetry) nondecoupling loop effects are needed for the occurrence of a strong first-order EWPT; however, these models have DM candidate. In [67] a nonminimal composite model based on the coset SO7/SO(6) has been considered. At low energy the scalar sector of their model is composed of two scalars, one with an unbroken Z2 symmetry and another scalar with a broken Z2 symmetry.
The plan of this paper is as follows.
In Section 2 we summarize the basic notions of EWPT. We describe the finite temperature effective potential at one loop. Then by emphasizing the underlying physical mechanisms, we describe the basic quantities of interest such as the strength of a phase transition, the rate of variation of bubble nucleation rate per volume, and the ratio of the latent heat released at the phase transition to the radiation energy density. In Section 3 we consider a simple extension of the SM by the addition of N real scalar gauge singlets with trilinear coupling to Higgs. We present the phenomenology of the model for N=2 and due to lack of protective symmetry these gauge singlets are not candidates for DM and we discuss the issue of observability of gravitational wave of our model. And finally in Section 4 we present our conclusions. Technical details are explained in the Appendix.
2. Electroweak Phase Transition
In this section we summarize basic notions and definitions of the electroweak phase transitions.
The observable universe consists predominantly of matter. The asymmetry between the matter and antimatter content of the universe is expressed by the baryon to photon ratio (1)ρ=nb-nb-nγ~10-9,where nb, nb-, and nγ are the number densities of baryons, antibaryons, and photons. In a symmetric universe one expects ρ=0 but experiments reveal that ρ has a tiny but nonzero value. This paradox can be resolved by requiring baryon number violation, C and CP violation, and departure from thermal equilibrium [4].
Baryogenesis is the physical process which is responsible for this observed baryon asymmetry of universe (BAU). In electroweak baryogenesis one assumes that the physical mechanism is the occurrence of a strong first-order electroweak phase transition (EWPT); namely, a smooth or a weak phase transition can not explain BAU [68, 69].
A convenient tool for investigation of EWPT is effective potential. For any quantum field theory if we replace the quantum field by its vacuum expectation value in the presence of a source the result for the potential energy to lowest order in perturbation theory is called the effective potential. The physical meaning of the effective potential is that it represents an energy density. In general effective potential contains other terms (loop corrections) [70, 71]. By using path integrals it is possible to find effective potential. In an Euclidean space-time at one-loop order the result is (2)Veffφc=V0φc+i2∫d4p2π4logp2+m2φcp2+Oħ2+⋯,which is known as the Coleman-Weinberg potential. The first term in (2) is the classical tree-level potential.
While considering physical processes in a hot environment such as early universe a proper method is finite temperature field theory [47–50]. The expression for the one-loop effective potential at finite temperature is [72] (3)V1effφc,T=∑iniT42π2J∓miφcT,where ni is the number of degrees of freedom of the particle, mi is the field dependent particle mass, and (4)J∓=±∫0∞dyy2log1∓exp-x2+y2,and J-(x) and J+(x) denote the contribution from bosons and fermions. In addition there is another contribution to thermal effective potential from the daisy subtraction [73, 74].
The best way for the evaluation of the finite temperature effective potential is to use special packages (codes) developed for a specific model.
But when the temperature is much greater than the masses of the various particles of the system under consideration it is possible to expand J∓ in a power series of m/T. In this work we use this high-temperature approximation. Recently, a new scheme for computation and resummation of thermal masses beyond the high-temperature approximation in general beyond SM scenarios has been proposed [75].
Let us consider the shape of the effective potential. For a generic model, as the universe cooled down at temperature above a critical temperature Tc, the effective potential had an absolute minimum which was located at the origin. At Tc there were two degenerate minima, which were separated by an energy barrier. At temperature below the critical temperature the second minimum became the global one and presently there is no energy barrier. Moreover the rate of expansion of the universe slowed, as the Hubble parameter H which characterizes the rate of expansion of the universe depends quadratically on the temperature.
In 1976 ’t Hooft discovered that baryon number is violated in the Standard Model and it is due to transition between two topologically distinct SU(2)L ground states [76, 77]. While, at zero temperature, the probability for barrier penetration is vanishingly small, at nonzero temperature the transition between two ground states differing by a unit of topological charge can be achieved by a classical motion over the barrier. Unstable static solutions of the field equations with energy equal to the height of the barrier separating two topologically distinct SU(2)L ground states (sphalerons) have been reported in [78]. Hence in the electroweak baryogenesis, the baryon asymmetry of universe is generated through the sphaleron process in the symmetric phase during the EWPT. At the symmetric phase the rate of baryon violating process which we denote by Γ~sph is a quartic function of the temperature; hence in this phase Γ~sph≫H.
But the third conditions to generate BAU is the departure from thermal equilibrium. Hence, the baryon number changing sphaleron interaction must quickly decouple in the broken phase; that is, Γ~sph<H. But the rate of sphaleron induced baryon violating process is suppressed by a Boltzmann factor (5)Γ~sph∝exp-EsphTT,but at phase transition Esph(T)∝φc(6)φcTc≥1,where φc is the broken phase minimum at the critical temperature Tc.
A first-order phase transition proceeds by nucleation of bubbles of the broken symmetry phase within the symmetric phase. The underlying mechanisms for bubble nucleation are quantum tunneling and thermal fluctuations. These bubbles then expand, merge, and collide. Gravitational waves are produced due to collision of bubble walls and turbulence in the plasma after the collisions. In addition while these bubbles pass through the plasma sound waves are created. These sound waves can provide additional sources of gravitational waves.
A crucial parameter for the calculation of the gravitational wave spectrum is the rate of variation of the bubble nucleation rate per volume, called β. It is common to use a normalized dimensionless parameter which is defined as (7)β~=βH∗,where H∗ denotes the Hubble parameter at the time of phase transition.
The density of latent heat released into the plasma is (8)ϵ∗=-VmineffT+TddTVmineffTT=T∗,where Vmineff(T) is the temperature-dependent true minimum of the effective potential of the scalar fields which causes the phase transition; moreover, its value must be set to zero by adding a constant at each time. Another dimensionless parameter for characterizing the spectrum of gravitational wave is (9)α=ϵ∗ρrad,where the radiation energy density of the plasma ρrad=π2/30g∗Tc4. And the parameter g∗ is the effective degrees of freedom in the thermal bath at the phase transition. In this work we assume g∗=106.75+NS, where NS denotes the number of singlet scalars that facilitates the electroweak phase transition.
3. The Model
In [19] the effects of a light scalar on the electroweak phase transition have been considered. In their model the scalar sector has been extended by an addition of a singlet, which has a trilinear interaction with the Higgs boson. Recently an extension of the SM with addition of N isospin-singlet, which has a quartic interaction with the Higgs boson, has been considered [57]. Here we consider a generalization of model of [19]. At zero temperature the effective potential of the scalar sector of our model is(10)V0=-DT02φ2+λ4φ4+12∑i=1Nmi2si2+φ2∑i=1Nκisi.Since we have not included quartic self-coupling for the extra scalars, in order to have a stable potential, we assume that the squares of the mass parameters of the extra scalars (mi2) are positive definite (see the Appendix for details).
At high temperature we have (11)VT=DT2-T02φ2-ETφ3+λT4φ4+12∑i=1Nmi2si2+φ2∑i=1Nκisi+T212∑i=1Nκisi.The parameters of (11) are given by(12)D=18v22mW2+mZ2+2mt2+2λv2,E=18πv34mW3+2mZ3,λT=λ-116π2v46mW4lnmw2aBT2+3mZ4lnmZ2aBT2-12mt4lnmt2aFT2,lnaB=3.91,lnaF=1.14.The critical temperature of model is given by (13)Tc=T01-E2/DλT-2ζ-ζ/12D,whereζ=∑i=1Nκi2mi2.The strength of the phase transition is denoted by ξ and it is given by (14)ξ=φcTc=2EλT-2ζ.Moreover, sic is the vev of the scalar field si at the second minimum of the effective potential at Tc:(15)sic=-κimi2φc2+Tc212.
3.1. Phenomenology of the Models with N Trilinear Interactions
The structure of the scalar mass matrix is(16)2λv22κ1v2κ2v2κ3v⋯⋯2κ1vm1200002κ2v0m220002κ3v0m32⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯.The physical mass squares of the scalar sector of our model are the eigenvalues of this matrix. The characteristic equation is(17)detω-2λv2-2κ1v-2κ2v-2κ3v⋯⋯-2κ1vω-m120000-2κ2v0ω-m22000-2κ3v0ω-m32⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯⋯=0.Hence in general one should solve a polynomial of degree (N+1) of ω. From the above equation we see that the coefficient of ωN+1 is unity. If we designate the coefficient of ωN by α, then -α is equal to the sum of physical masses of the model; namely, (18)2λv2+∑i=1Nmi2=mH2+∑i=1Nχi2,where the physical mass of the ith scalar is denoted by χi. From this relation we find an important relation about the Higgs self-coupling; namely, (19)λ-λSM=12v2∑i=1Nχi2-mi2,where λSM is the Higgs quartic self-coupling of the Standard Model.
We see that if the mass parameters of the scalars as well as the physical masses of the scalars are much smaller than v, then the deviation of the parameter λ from the Higgs self-coupling of the Standard Model will be very small.
3.2. The Special Case N=2
Here we consider the case N=2. The scalar mass matrix of the model at zero temperature for this case is(20)M2=2λv22κ1v2κ2v2κ1vm1202κ2v0m22.The physical mass squares of the model can be obtained from (21)ω3+Aω2+Bω+C=0,where (22)A=-2λv2+m12+m22,B=2v2λm12+m22-2κ12+κ22+m12m22,C=4v2κ12m22+κ22m12-2λ2v2m12m22.Now one of the eigenvalues is equal to mH2; hence we obtain(23)mH6+AmH4+BmH2+C=0.Moreover, by minimizing the potential with respect to variables s1, s2, we obtain (24)si=-κiv2mi2,i=1,2.By minimizing the potential with respect to variable φ we obtain (25)-2DT02+λv2-2∑i=12κi2v2mi2=0.The parameters of the model must satisfy (23) and (25). Hence the parameter space of the model in this case contains four independent parameters (κ1,κ2,m1,m2).
If we subtract (23) from (21) we get (26)χ4+A+mH2χ2+mH4+mH2A+B=0.Therefore, the physical mass squares of the singlets are determined.
But models with extended Higgs sectors predicting strongly first-order phase transition simultaneously predict a significant deviation in the triple Higgs boson coupling as well [79]. This deviation at the tree level in [17] and at loop level in [33, 79] has been studied, with (27)Δhhh=λhhhMSM-λhhhSMλhhhSM,where λhhhSM is the Higgs triple coupling of the SM and λhhhMSM is the Higgs triple coupling of the multisinglet extension of the SM. Collider experiments could measure the Higgs triple coupling. Here we want to explore the region of intermediate mass of the singlets. But there are bounds on Higgs-Portal models from the LHC Higgs data [80–82]. For instance, the Higgs doublet is mixed with the extra singlet scalars and the mixing element cos(φ) between the CP-even component of the doublet (φ) and the physical Higgs (whose mH=125.09 GeV) is not arbitrary and it is subject of a constraint coming from the Higgs coupling to the W gauge bosons. Current data [81] suggests cos(φ)>0.86.
The results are presented in Table 1 for the onset of a strong EWPT; namely, ξ=1. For each configuration in the table we present our results for the deviation of Higgs triple coupling as well. Hence by adding one scalar to the model we can have singlets in the intermediate mass region.
Different configurations associated with the onset of a strong first-order phase transition (ξ=1).
Set
κ1
κ2
m1 (GeV)
m2 (GeV)
χ1 (GeV)
χ2 (GeV)
cos(φ)
Δλhhh
Δλhhhh
I
0.2
3.7
69.2
16.0
6.2
69.2
0.993
13.9%
4.1%
II
0.5
3.1
18.7
13.5
5.2
18.6
0.995
14.5%
3.0%
III
0.8
2.3
11.2
10.4
4.0
11.1
0.997
15.1%
1.8%
IV
1.2
0.7
18.1
3.2
17.5
1.3
0.999
5.7%
0.6%
For completeness the predictions of the model for the deviation of Higgs boson quartic coupling Δλhhhh is given in Table 1.
Moreover, the existence of the extra singlet scalars could affect the total Higgs decay if they are light enough, which becomes (28)Γtotalh=cos2φΓtotalhSM+∑Γh⟶si+sk,where s denote all the scalars. The deviation from the SM value ΔΓtotal<1.4(MeV).
In the SM, a total Higgs decay width around 4 MeV is predicted. In this work we assume the decay h→si+sk is an invisible decay. However, current analysis [82] suggests that the Higgs invisible decay branching ratio should be less than 17%. In Table 2 we present Higgs invisible decay branching ratio in various decay modes, as well as the total Higgs invisible branching ratio for all of Higgs invisible decay modes. The invisible Higgs width, total Higgs width, and the deviation of total Higgs width of our model from that of the SM are also shown. The unit for the various width in this table is MeV; moreover in our calculation we have assumed Γtotal(hSM)=4(MeV).
Values of Higgs invisible branching ratios for various decay modes, total value of Higgs invisible branching ratio, Higgs invisible width, Higgs total width, and deviation of the Higgs width from the SM value are shown for various configurations. The unit for width in this table is MeV. All four configurations are consistent with current experimental data.
Set
B(h→s1s1)
B(h→s1s2)
B(h→s2s2)
Btot(h→inv)
Γinv(h)
Γtot(h)
ΔΓtot
I
1.45%
0.31%
0.0%
1.49%
0.06
4.006
0.006
II
0.75%
0.04%
3.6×10-5%
0.8%
0.03
3.992
0.008
III
0.24%
0.004%
1.3×10-6%
.24%
0.01
3.99
0.01
IV
0.02%
4.6×10-7%
6.3×10-12%
0.02%
5.8×10-4
3.993
0.007
It turns out that the critical temperature for this model Tc<100(GeV). Hence it is a good approximation to use λ instead of λT and in [19] this approximation is used to study SFPT for the case N=1 but the mass of the light scalar is up to 12 GeV; in [20] a one-loop study of the same model has been presented but the mass of the light scalar to catalyze a SFPT is up to 20 GeV. But for the model presented in this work and using this approximation for the case N=2 the mass of the scalar to catalyze a SFPT is 69 GeV. By adding more scalars we expect to have heavy singlets.
3.3. Detection of Gravitational Waves
The direct detection of gravitational waves by LIGO [83] generated a lot of interest among researchers in cosmology, astrophysics, and particle physics. Major sources of gravitational waves are inflation, compact binary systems, or cosmological phase transitions.
In this section under phenomenological description we aim to make an estimate for the observability of the gravitational waves produced by the model presented in this work.
In all of previous multisinglet models a Z2 symmetry (either broken or unbroken) is imposed on the fields. However, neither of the two fields of our model has this symmetry. This is the main difference between our model and previous studies.
The main parameters that are of value for obtaining the spectrum of gravitational waves are the parameters α and β, which we discussed in Section 2. And they are obtained from the thermal effective potential.
For the set I of Table 1 and from (8) and (9) we obtain α=0.16. But the velocity of the bubble wall and the efficiency factor, the fraction of the latent heat which is converted to the kinetic energy of the plasma, are determined from (29)vb=1/3+α2+2α/31+α,κ=11+0.715α0.715α+4273α2.Thus for this configuration vb=0.81 and κ=0.17. The peak frequency for bubble collision contribution is given by [84] (30)fcol=16.5×10-60.62vb2-0.1vb+1.8βH∗T∗100g∗1001/6.Hz.And the energy density at this frequency is given by (31)Ωhcol2=1.67×10-5βH∗20.11vb30.42+vb2κα1+α2g∗100-1/3,where h is the reduced Hubble constant at present.
Another source of gravitational waves is the compression waves in the plasma (sound waves) and the peak frequency is (32)fsw=1.9×10-5βH∗vb-1T∗100g∗1001/6,Hz;the energy density at this peak frequency is obtained from (33)Ωhsw2=2.65×10-6βH∗-1κα1+α2g∗100-1/3vb.And finally the peak frequency for gravitational waves which are caused by the turbulence of the plasma is (34)fturbo=2.7×10-5βH∗vb-1T∗100g∗1001/6Hz,and the peak energy density of this part is given by (35)Ωhturbo2=3.35×10-4βH∗-1εκα1+α3/2g∗100-1/3vb1211/31+8πfturbo/H∗,where ε denotes the fraction of latent heat that is transformed into turbulent motion of the plasma. We choose ε=0.05.
Hence by knowing the values of the parameters α, β~, and T∗, we can obtain the spectra of the GW. Even though the nucleation temperature T∗ is lower than the critical temperature in this work we assume T∗≈Tc.
The standard method of calculation of the parameter β~ from the effective potential is to compute the Euclidean action of the model and it is explained in [36, 85], but in an approximate scheme it is found that [59] (36)β~≈170-4lnT∗1GeV-2lng∗,and in our case when the phase transition happens at weak scale, β~≈144.
Hence, in order to assess the implications of the model on the spectra of GW we present results by varying this parameter in the interval 50≤β~≤250. In Table 3 we present our results. The unit for various frequencies is mHz.
The spectra of gravitational wave as predicted by our model. The emitted GW are within the reach of eLISA C1.
GW spectra
fsw
fcol
fturb
Ωhsw2
Ωhcol2
Ωhturb2
Ωhtotal2
β~=50
0.65
0.12
1.85
2.34×10-11
1.98×10-13
1.67×10-15
2.36×10-11
β~=100
1.30
0.24
3.69
1.17×10-11
9.9×10-14
8.35×10-16
1.18×10-11
β~=250
3.25
0.60
9.23
4.67×10-12
3.96×10-14
3.34×10-16
4.71×10-12
The results of Table 3 show that in our model the contribution from the turbulent motion of always a few orders of magnitude is smaller than the previous two. Moreover, the contribution from the sound waves is the dominant source of the total GW spectrum. By studying the frequency dependent spectra [84], we find that the peak of energy density of the sound wave contribution and the peak of energy density due to collision are well separated, a desirable feature while detecting these waves. These waves are within the reach of future gravitational wave interferometers (eLISA C1).
4. Conclusions
In this work we have presented a new extension of the SM. In this model we amend the SM by N gauge singlets. In order to avoid proliferation of the parameters we considered the most economical model. For each scalar we allowed a mass term with a positive squared mass parameter to insure vacuum stability in the direction of that scalar and a triple coupling with the Higgs field to facilitate a strong electroweak phase transition, and for the special case N=2 we find SFPT with gauge singlets in the intermediate mass range. We also investigated the deviations of the Higgs coupling constants from the SM values. And we find that the deviation of the triple Higgs boson coupling can be as large as 15%.
Finally, we have obtained the gravitational wave spectrum from the electroweak phase transition. We have shown that the gravitational wave signal can be detected by eLISA. The present model has a large parameter space and as a result a richer phenomenology in comparison to [19]. It would be of interest to amend the model by inclusion of quartic self-coupling of the scalars as well as the mixing term between scalars (see the Appendix for details).
It would be of interest to consider higher values of N for the model proposed in this work. These and other related issues are presently under consideration.
AppendixVacuum Stability Conditions
For the special case N=2 of the model the weak condition for vacuum stability is that if the value of the fields tends to infinity, then the tree-level potential is not unbounded-from-below directions. This leads to λ>0. But if only the field s1 tends to infinity, then the condition vacuum stability will be m12>0 and by similar argument for the field s2 the restriction will be m22>0.
However, the strong condition for vacuum stability as stated in [35, 86, 87] is that the masses of the particles of a particular model are positive.
The necessary conditions for a symmetric matrix A of order 3 to have real eigenvalues are (A.1)a11>0,a22>0,a33>0,a-12=a12+a11a22>0,a-13=a13+a11a33>0,a-23=a23+a22a33>0,(A.2)a11a22a33+a12a33+a13a22+a23a11+2a-12a-13a-23>0.The constraints stated in (A.1) apply to quantities with dimension of squared mass. From our mass matrix of Section 3.2, the first constraint is 2λv2>0, which leads to λ>0. Hence, for the model presented in this work, the criteria which provide the necessary and sufficient vacuum stability conditions are given by (A.3)λ>0,m12>0,m22>0,κ1>-λm122,κ2>-λm222,(A.4)λm12m22+κ12m22+κ22m12+2κ1+2λm122κ2+2λm22m12m22>0,where the conditions stated in the first line of (A.3) have been obtained by using weak conditions for vacuum stability.
In our model we have not included quartic terms for the gauge singles such as (A.5)Vquartic=∑i=12λisi44+δs12s22in the effective potential, where λ1, λ2, and δ are dimensionless coupling parameters; in this case the vacuum stability in the direction of the extra scalars is maintained if (A.6)λ1>0,λ2>0,λ1λ2>-2δ,and the squares of the mass parameters m12 and m22 in principle can assume any value (positive, zero, or negative); for instance, in the multiscalar model described in [26] all of the extra scalars do not have a mass term, while in the absence of quartic terms the squared mass parameters of the model presented in Section 3 must be positive as required by vacuum stability.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
AadG.AbajyanT.AbbottB.Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC2012716112910.1016/j.physletb.2012.08.020ChatrchyanS.KhachatryanV.SirunyanA. M.Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC20127161306110.1016/j.physletb.2012.08.021AndersonG. W.HallL. J.Electroweak phase transition and baryogenesis19924582685269810.1103/PhysRevD.45.26852-s2.0-0000198042SakharovA.Violation of CP invariance, c asymmetry, and baryon asymmetry of the universe19675242710.1070/PU1991v034n05ABEH002497SilveiraV.ZeeA.Scalar phantoms19851611–313614010.1016/0370-2693(85)90624-02-s2.0-0000831993McDonaldJ.Gauge singlet scalars as cold dark matter19945063637364910.1103/PhysRevD.50.36372-s2.0-0000802841BurgessC. P.PospelovM.Ter VeldhuisT.The minimal model of nonbaryonic dark matter: a singlet scalar20016191–370972810.1016/s0550-3213(01)00513-22-s2.0-0035944988van der BijJ.The minimal non-minimal standard model20066361565910.1016/j.physletb.2006.03.018BargerV.LangackerP.McCaskeyM.Ramsey-MusolfM. J.ShaughnessyG.CERN LHC phenomenology of an extended standard model with a real scalar singlet200877303500510.1103/physrevd.77.0350052-s2.0-41049116164GuoW.-L.WuY.-L.The real singlet scalar dark matter model20101010, article 083YagunaC. E.The singlet scalar as FIMP dark matter20111108, article 6010.1007/JHEP08(2011)060CoimbraR.SampaioM. O.SantosR.ScannerS: constraining the phase diagram of a complex scalar singlet at the LHC2013735242810.1140/epjc/s10052-013-2428-4FengL.ProfumoA. S.UbaldiB. L.Closing in on singlet scalar dark matter: LUX, invisible Higgs decays and gamma-ray lines201520153, article 4510.1007/JHEP03(2015)045EspinosaJ.QuirósM.The electroweak phase transition with a singlet19933051-29810510.1016/0370-2693(93)91111-YChoiJ.VolkasR.Real Higgs singlet and the electroweak phase transition in the standard model1993317338539110.1016/0370-2693(93)91013-DProfumoS.Ramsey-MusolfM. J.ShaughnessyG.Singlet Higgs phenomenology and the electroweak phase transition20070708, article 01010.1088/1126-6708/2007/08/010NobleA.PerelsteinM.Higgs self-coupling as a probe of the electroweak phase transition200878606351810.1103/physrevd.78.063518AshoorioonA.KonstandinT.Strong electroweak phase transition without collider traces20090907, article 086DasS.FoxP. J.KumarA.WeinerN.The dark side of the electroweak phase transition2010201011, article 108EspinosaJ. R.KonstandinT.RivaF.Strong electroweak phase transitions in the standard model with a singlet2012854359263010.1016/j.nuclphysb.2011.09.0102-s2.0-80054058868AlanneT.TuominenK.VaskonenV.Strong phase transition, dark matter and vacuum stability from simple hidden sectors201488969271110.1016/j.nuclphysb.2014.11.001ZBL1326.812562-s2.0-84909619162CurtinD.MeadeP.YuC.Testing electroweak baryogenesis with future colliders20142014, article 12710.1007/jhep11(2014)127ProfumoS.Ramsey-MusolfM. J.WainwrightC. L.WinslowP.Singlet-catalyzed electroweak phase transitions and precision Higgs boson studies201591303501810.1103/physrevd.91.0350182-s2.0-84923338406ClineJ. M.KainulainenK.ScottP.WenigerC.Update on scalar singlet dark matter201388055025ClineJ. M.KainulainenK.Electroweak baryogenesis and dark matter from a singlet Higgs20132013, article 01210.1088/1475-7516/2013/01/012EspinosaJ. R.QuirósM.Novel effects in electroweak breaking from a hidden sector200776707600410.1103/physrevd.76.076004EspinosaJ. R.KonstandinT.NoJ. M.QuirósM.Some cosmological implications of hidden sectors2008781212352810.1103/physrevd.78.123528DrozdA.GrzadkowskiB.WudkaJ.Cosmology of multi-singlet-scalar extensions of the standard model201142112255226210.5506/aphyspolb.42.22552-s2.0-82455219286AbadaA.GhafforD.NasriS.Two-singlet model for light cold dark matter201183909502110.1103/physrevd.83.0950212-s2.0-79960771895AbadaA.NasriS.Phenomenology of a light cold dark matter two-singlet model201285707500910.1103/physrevd.85.075009DrozdA.GrzadkowskiB.WudkaJ.Multi-scalar-singlet extension of the standard model—the case for dark matter and invisible Higgs boson201220124, article 006AhricheA.NasriS.Light dark matter, light Higgs boson, and the electroweak phase transition201285909300710.1103/physrevd.85.093007AhricheA.ArhribA.NasriS.Higgs phenomenology in two-singlet model201420142, article 042TofighiA.GhodsiO. N.SaeedhoseiniM.Phase transition in multi-scalar-singlet extensions of the Standard Model201574820821110.1016/j.physletb.2015.07.0092-s2.0-84937893011ModakK. P.MajumdarD.RakshitS.A possible explanation of low energy γ-ray excess from galactic center and fermi bubble by a dark matter model with two real scalars201520153, article 011JinnoR.NakayamaK.TakimotoM.Gravitational waves from the first order phase transition of the Higgs field at high energy scales201693404502410.1103/PhysRevD.93.045024MR3499412ClineJ. M.LemieuxP.-A.Electroweak phase transition in two Higgs doublet models19975563873388110.1103/physrevd.55.38732-s2.0-0009043025FrommeL.HuberS. J.SeniuchM.Baryogenesis in the two-Higgs doublet model20060611, article 03810.1088/1126-6708/2006/11/038HuberS. J.SchmidtM. G.SUSY variants of the electroweak phase transition199910347348110.1007/s100520050770FunakuboK.SenahaE.Electroweak phase transition, critical bubbles, and sphaleron decoupling condition in the MSSM2009791111502410.1103/physrevd.79.115024ChungD. J. H.LongA. J.Electroweak phase transition in the μν SSM20108112353110.1103/PhysRevD.81.123531KozaczukJ.ProfumoS.HaskinsL. S.WainwrightC. L.Cosmological phase transitions and their properties in the NMSSM20152015, article 14410.1007/jhep01(2015)144BalázsC.MazumdarA.PukartasE.WhiteG.Baryogenesis, dark matter and inflation in the next-to-minimal supersymmetric standard model20142014, article 07310.1007/jhep01(2014)073HuangW.KangZ.ShuJ.WuP.YangJ. M.New insights in the electroweak phase transition in the NMSSM201591202500610.1103/PhysRevD.91.025006ChiangC.-W.SenahaE.Electroweak phase transitions in the secluded U(1)-prime-extended MSSM201020106, article 030ChungD. J.LongA. J.WangL.-T.125 GeV Higgs boson and electroweak phase transition model classes201387202350910.1103/physrevd.87.023509DolanL.JackiwR.Symmetry behavior at finite temperature19749123320334110.1103/PhysRevD.9.33202-s2.0-26544432182KolbE. W.TurnerM. S.The early universe199069, article 1547KapustaJ. I.1989Cambridge, UKCambridge University PressCambridge Monographs on Mathematical PhysicsMR1110583Le-BellacM.2000Cambridge University PressKosowskyA.MackA.KahniashviliT.Gravitational radiation from cosmological turbulence20026621702403010.1103/physrevd.66.024030MR19366322-s2.0-0037101305KosowskyA.TurnerM. S.WatkinsR.Gravitational radiation from colliding vacuum bubbles199245124514453510.1103/PhysRevD.45.45142-s2.0-26344449570KamionkowskiM.KosowskyA.TurnerM. S.Gravitational radiation from first-order phase transitions19944962837285110.1103/PhysRevD.49.28372-s2.0-24244452960CapriniC.DurrerR.Gravitational waves from stochastic relativistic sources: primordial turbulence and magnetic fields200674606352110.1103/PhysRevD.74.063521LeitaoL.MegevandA.Gravitational waves from a very strong electroweak phase transition2016201610.1088/1475-7516/2016/05/037ArtymowskiM.LewickiM.WellsJ. D.Gravitational wave and collider implications of electroweak baryogenesis aided by non-standard cosmology20172017, article 06610.1007/jhep03(2017)066KakizakiM.KanemuraS.MatsuiT.Gravitational waves as a probe of extended scalar sectors with the first order electroweak phase transition2015921111500710.1103/physrevd.92.115007HuangP.LongA. J.WangL.-T.Probing the electroweak phase transition with Higgs factories and gravitational waves201695075008BalázsC.FowlieA.MazumdarA.WhiteG. A.Gravitational waves at aLIGO and vacuum stability with a scalar singlet extension of the standard model201795404350510.1103/physrevd.95.043505KatzA.PerelsteinM.Higgs couplings and electroweak phase transotion201420147, article 108KrogJ.MojazaM.SanninoF.Four-fermion limit of gauge-Yukawa theories201592804501510.1103/physrevd.92.085043WangZ.SteeleT. G.HanifT.MannR. B.Conformal complex singlet extension of the Standard Model: scenario for dark matter and a second Higgs boson201620168, article 06510.1007/jhep08(2016)065LewickiM.Rindler-DallerT.WellsJ. D.Enabling electroweak baryogenesis through dark matter201620166, article 055KotwalA. V.Ramsey-MusolfM. J.NoJ. M.WinslowP.Singlet-catalyzed electroweak phase transitions in the 100 TeV frontier201694303502210.1103/physrevd.94.035022BaldesI.KonstandinT.ServantJ.A First-Order Electroweak Phase Transition in the Standard Model from Varying Yukawashttps://arxiv.org/abs/1604.04526TenkanenT.TuominenK.VaskonenV.A strong electroweak phase transition from the inflaton field20162016, article 037ChalaM.NardinibG.SobolevcI.Unified explanation for dark matter and electroweak baryogenesis with direct detection and gravitational wave signatures20169405500610.1103/physrevd.94.055006KuzminV. A.RubakovV. A.ShaposhnikovM. E.On anomalous electroweak baryon-number non-conservation in the early universe19851551-2364210.1016/0370-2693(85)91028-72-s2.0-0003017602MorrisseyD. E.Ramsey-MusolfM. J.Electroweak baryogenesis20121412500310.1088/1367-2630/14/12/1250032-s2.0-84871910908ColemanS.1988Cambridge, UKCambridge University PressZeeA.20102ndPrinceton, NJ, USAPrinceton University PressMR2797348CarenaM.MegevandA.QuirósM.WagnerC. E. M.Electroweak baryogenesis and new TeV fermions20057161-231935110.1016/j.nuclphysb.2005.03.0252-s2.0-18344381984ArnoldP. B.EspinosaO.Effective potential and first-order phase transitions: beyond leading order19934783546357910.1103/physrevd.47.35462-s2.0-0001645345FodorZ.HebeckerA.Finite temperature effective potential to order g4, λ2 and the electroweak phase transition19944321-212714610.1016/0550-3213(94)90596-72-s2.0-0001457711CurtinD.MeadeP.RamaniH.Thermal Resummation and Phase Transitionshttps://arxiv.org/abs/1612.00466't HooftG.Symmetry breaking through bell-jackiw anomalies197637181110.1103/physrevlett.37.8't HooftG.Computation of the quantum effects due to a four-dimensional pseudoparticle197614123432345010.1103/PhysRevD.14.34322-s2.0-4243603598MantonN. S.Topology in the weinberg-salam theory19832882019202610.1103/physrevd.28.2019MR7209932-s2.0-0001571293KanemuraS.KikuchiM.YagyuK.One-loop corrections to the Higgs self-couplings in the singlet extension201791715417710.1016/j.nuclphysb.2017.02.004MR3627324RobensT.StefaniakT.LHC benchmark scenarios for the real Higgs singlet extension of the standard model201676, article 26810.1140/epjc/s10052-016-4115-8CheungK.KoP.LeeJ. S.TsengP.Bounds on Higgs-portal models from the LHC Higgs data2015201510, article 05710.1007/jhep10(2015)057BechtleP.HeinemeyerS.StålO.StefaniakT.WeigleinG.Probing the Standard Model with Higgs signal rates from the Tevatron, the LHC and a future ILC2014201411, article 3910.1007/jhep11(2014)0392-s2.0-84919933842CapriniC.HindmarshM.HuberS.KonstandinT.KozaczukJ.NardiniG.NardiniG.NoJ. M.PetiteauA.SchwallerP.ServantG.WeirD. J.Observation of gravitational waves from a binary black hole merger2016116CapriniC.HindmarshM.HuberS.Science with the space-based interferometer eLISA. II: gravitational waves from cosmological phase transitions20162016, article 001VaskonenV.Electroweak baryogenesis and gravitational waves from a real scalar singlethttps://arxiv.org/abs/1611.02073KannikeK.Vacuum stability conditions from copositivity criteria2012727, article 209310.1140/epjc/s10052-012-2093-zChakraborttyJ.KonarP.MondalT.Copositive criteria and boundedness of the scalar potential201489909500810.1103/physrevd.89.095008