Advances in High Energy Physics

Volume 2017 (2017), Article ID 7638204, 10 pages

https://doi.org/10.1155/2017/7638204

## Strong Electroweak Phase Transition in a Model with Extended Scalar Sector

Department of Physics, Faculty of Basic Sciences, University of Mazandaran, P.O. Box 47416-95447, Babolsar, Iran

Correspondence should be addressed to Ali Tofighi; ri.ca.zmu@ihgifot.a

Received 4 January 2017; Revised 16 March 2017; Accepted 30 March 2017; Published 30 April 2017

Academic Editor: Enrico Lunghi

Copyright © 2017 Mina Saeedhosseini and Ali Tofighi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.

#### Abstract

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 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 . 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 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], real gauge singlets are coupled to Higgs boson via trilinear interactions. Previous studies of multiscalar singlet extension of SM impose separate 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 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 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 has been considered. At low energy the scalar sector of their model is composed of two scalars, one with an unbroken symmetry and another scalar with a broken 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 real scalar gauge singlets with trilinear coupling to Higgs. We present the phenomenology of the model for 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 where , , and are the number densities of baryons, antibaryons, and photons. In a symmetric universe one expects 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 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] where is the number of degrees of freedom of the particle, is the field dependent particle mass, and and and 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 in a power series of . 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 , the effective potential had an absolute minimum which was located at the origin. At 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 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 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 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 is a quartic function of the temperature; hence in this phase .

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, . But the rate of sphaleron induced baryon violating process is suppressed by a Boltzmann factor but at phase transition where is the broken phase minimum at the critical temperature .

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 where denotes the Hubble parameter at the time of phase transition.

The density of latent heat released into the plasma is where 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 where the radiation energy density of the plasma . And the parameter is the effective degrees of freedom in the thermal bath at the phase transition. In this work we assume , where 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 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 isSince 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 () are positive definite (see the Appendix for details).

At high temperature we have The parameters of (11) are given byThe critical temperature of model is given by The strength of the phase transition is denoted by and it is given by Moreover, is the of the scalar field at the second minimum of the effective potential at :

##### 3.1. Phenomenology of the Models with Trilinear Interactions

The structure of the scalar mass matrix isThe physical mass squares of the scalar sector of our model are the eigenvalues of this matrix. The characteristic equation isHence in general one should solve a polynomial of degree of . From the above equation we see that the coefficient of is unity. If we designate the coefficient of by , then is equal to the sum of physical masses of the model; namely, where the physical mass of the th scalar is denoted by . From this relation we find an important relation about the Higgs self-coupling; namely, where 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 , then the deviation of the parameter from the Higgs self-coupling of the Standard Model will be very small.

##### 3.2. The Special Case

Here we consider the case . The scalar mass matrix of the model at zero temperature for this case isThe physical mass squares of the model can be obtained from where Now one of the eigenvalues is equal to ; hence we obtainMoreover, by minimizing the potential with respect to variables , , we obtain By minimizing the potential with respect to variable we obtain The parameters of the model must satisfy (23) and (25). Hence the parameter space of the model in this case contains four independent parameters ().

If we subtract (23) from (21) we get 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 where is the Higgs triple coupling of the and is the Higgs triple coupling of the multisinglet extension of the . 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 between the CP-even component of the doublet () and the physical Higgs (whose 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 .

The results are presented in Table 1 for the onset of a strong EWPT; namely, . 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.