Next Article in Journal
Thermofractals, Non-Additive Entropy, and q-Calculus
Next Article in Special Issue
The Merger Rate of Black Holes in a Primordial Black Hole Cluster
Previous Article in Journal
Light-Induced Effects in Amorphous Chalcogenide Glasses: Femtoseconds to Seconds
Previous Article in Special Issue
Scattering on Quasi-Spherical Black-Holes: Features and Beyond
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy Production Due to Electroweak Phase Transition in the Framework of Two Higgs Doublet Model

1
Department of Physics and Astronomy, Novosibirsk State University, Lyapunova 2, 630090 Novosibirsk, Russia
2
Institute of Physics, Southern Federal University, 344090 Rostov on Don, Russia
3
CNRS, Astroparticule et Cosmologie, Université de Paris, F-75013 Paris, France
4
Center for Cosmopartilce Physics “Cosmion”, National Research Nuclear University “MEPhI” (Moscow Engineering Physics Institute), 115409 Moscow, Russia
*
Author to whom correspondence should be addressed.
Submission received: 4 March 2021 / Revised: 8 April 2021 / Accepted: 19 April 2021 / Published: 29 April 2021
(This article belongs to the Special Issue Beyond the Standard Models of Physics and Cosmology)

Abstract

:
We revisit the possibility of first order electroweak phase transition (EWPT) in one of the simplest extensions of the Standard Model scalar sector, namely the two-Higgs-doublet model (2HDM). We take into account the ensuing constraints from the electroweak precision tests, Higgs signal strengths and the recent LHC bounds from direct scalar searches. By studying the vacuum transition in 2HDM, we discuss in detail the entropy released in the first order EWPT in various parameter planes of a 2HDM.

1. Introduction

It is a well-established fact that electroweak phase transition (EWPT) is either a second order or a smooth crossover in the Standard Model (SM) of particle physics. It is also well-established that the entropy density in the early universe plasma is conserved in the course of the cosmological expansion if the plasma is in thermal equilibrium state with negligible chemical potential of every species [1,2]. The entropy conservation law is given by
s = P + ρ T a 3 = const .
where s is the entropy density, a ( t ) is the scale factor, T ( t ) is the temperature of the fluid (or plasma) are function of time t, ρ and P are the energy density and pressure of the plasma, respectively.
In the early universe, the state of matter is quite close to the equilibrium as the reaction rate Γ n σ v , n = 1 , 2 , represents decays of one and two-body reactions, respectively, σ v is the product of annihilation cross section, is much faster than the cosmological expansion rate, i.e., the Hubble parameter H = a ˙ / a T 2 / m P l , where a ˙ is the time derivative of the scale factor. The equilibrium condition Γ > H is always satisfied for at temperature T < α m Pl . Here α is the coupling constant of the particle interaction of the order of 10 2 and m Pl is the Planck mass. Due to the large value of m Pl , thermal equilibrium exists in most of the history of the universe, if α 1 .
As mentioned above, during thermal equilibrium, the entropy density in the comoving volume is conserved. However, there are scenarios where the entropy density is not conserved. For example, if the universe at a certain stage was dominated by primordial black holes [3], the entropy production can be very high, high enough to delete the pre-existing baryon asymmetry [4]. In the context of the modern cosmological paradigm of inflationary Universe with baryosynthesis and dark matter/energy, physics beyond the Standard Model (BSM) underlying these necessary elements of the modern cosmological model can provide many examples of various mechanisms of high entropy production (see e.g., [5] for review and references). Taking apart the wide range of various possibilities we consider here the problem of entropy production by minimal extension of SM and start the discussion from the SM predictions for the cosmological entropy production.
A large entropy production could take place during during quantum chromodynamical (QCD) phase phase transition at T 100–200 MeV. However, due to difficulty in numerical computation, QCD phase transition in early universe cosmology is not known in detail. For reference, see Ref. [6].
Few mechanisms of realistic though very weak entropy production could take place during the freeze-out of dark matter (DM) particles. However, usually, the fraction of DM density was quite low at the freezing out temperature and the effect is tiny.
An interesting effect, not covered in this paper, is the formation of bubble walls that can take place in the early universe. Their collision can lead to the formation of primordial black holes due to first order phase transition with background gravitational waves [7,8]. The entropy released in the process can wipe out the dark matter density that was present before EWPT. Recent results on muon ( g 2 ) experiment, see Refs. [9,10,11], have opened up a new window to study dark matter in the framework of extended SM physics which shall be done in our next work.
Most probably, the largest entropy release in the standard model took place in the process of the electroweak transition from symmetric to asymmetric electroweak phase in the course of the cosmological cooling down. In the SM with one Higgs field, the process is a mild crossover and the entropy production is about 13 % [12].
According to the electroweak (EW) theory at the temperature higher than a critical one, T > T c , the expectation value of the Higgs field, ϕ , in the fluid (plasma) is zero and the universe is in electroweak symmetric phase [13].When the temperature drops below T c , a nonzero expectation value is created, which gradually rises, with decreasing temperature, up to the vacuum expectation value η . Such a state does not satisfy the conditions necessary for the entropy conservation and an entropy production is expected.
A huge amount of entropy is released if EWPT is first order, which is the case even with the minimalist extension of standard model known as two-Higgs-doublet Model (2HDM). In what follows, we have considered a real 2HDM of type-I and are scanned over certain parameter spaces and used numerical analysis to calculate the entropy production for some interesting and unique benchmark points.
The paper is arranged as follows: in the next section details about 2HDM are given along with some some Large Hadron Collider (LHC) constraints, followed by the theoretical framework of the process. Due to cumbersome and very difficult analytical calculations, we did numerical analysis of certain parameters using the BSMPT package [14]. A generic discussion and conclusion is given there after. The paper closes by two appendices that give details about the metric being used here and also the masses of the scalar bosons generated by 2HDM.

2. 2HDM: A Brief Review and Current Constraints

There are two scalar doublets in the framework and they are defined as:
φ I = ϕ I + 1 2 ( v I + ρ I ) + i η I ,
with I = 1 , 2 . Here ϕ I ± , ρ I , η I , and v I indicate the charged, neutral CP-even and neutral CP-odd degrees of freedom (d.o.f.) and the vacuum expectation value (VEV) of the I-th doublet, respectively.
Prior to spontaneous symmetry breaking (SSB), the tree-level 2HDM Lagrangian, comprising of 6-dimensional operators, assumes the form
L = L kin + L Yuk V ( φ 1 , φ 2 ) + L 6 ,
where,
L kin = 1 4 X = G a , W i , B X μ ν X μ ν + I = 1 , 2 | D μ φ I | 2 + ψ = Q , L , u , d , l ψ ¯ i D ψ , L Yuk = I = 1 , 2 Y I e l ¯ e φ I + I = 1 , 2 Y I d q ¯ d φ I + I = 1 , 2 Y I u q ¯ u φ ˜ I , V ( φ 1 , φ 2 ) = m 11 2 | φ 1 | 2 + m 22 2 | φ 2 | 2 ( μ 2 φ 1 φ 2 + h . c . ) + λ 1 | φ 1 | 4 + λ 2 | φ 2 | 4 + λ 3 | φ 1 | 2 | φ 2 | 2 + λ 4 | φ 1 φ 2 | 2 + λ 5 2 φ 1 φ 2 + λ 6 | φ 1 | 2 + λ 7 | φ 2 | 2 φ 1 φ 2 + h . c . , L 6 = i c i O i / f 2 .
L kin is the kinetic term of the Lagrangian, L Yuk is the Lagrangian originating from the Yukawa interaction. The indices μ , ν = 0 , 1 , 2 , 3 are the indexes for the time-space components. W i , i = 1 , 2 , 3 , and B are the four gauge bosons, G a ( a = 1 , 2 ) are the Goldstone bosons and Y e , Y d and Y u are the Yukawa coupling constants. D is the covariant derivative, see Ref. [15], and q , u , d , e , l , q , m 11 , m 22 are defined in Appendix B. The term ‘h.c.’ and the symbol † stay for the Hermitian conjugate. c i is the Wilson coefficient of the 6-dimensinal operator O i and f. The terms proportional to λ 6 , 7 are known as ‘hard- Z 2 violating’ because, not only do they lead to a quadratically divergent amplitude for φ 1 φ 2 transition [16] but also to CP-violation in the scalar sector for complex values [17]. However, it is possible to realize the CP-conserving limit with nonzero values of λ 6 , 7 as well [18]. In this paper, we contain our discussion to the CP-conserving 2HDM, and hence λ 6 , 7 = 0 . The electroweak symmetry is broken by the VEVs, namely v 1 and v 2 corresponding to the two doublets φ 1 , 2 respectively. This leads to the mixing of same types of degrees of freedom of φ 1 , 2 . For the case of CP-conservation, the mass matrices of the neutral CP-even and odd scalars and the charged scalars are diagonalized by the following operation:
H h = R ( α ) ρ 1 ρ 2 , W L ± H ± = R ( β ) ϕ 1 ± ϕ 2 ± , Z L A = R ( β ) η 1 η 2 .
Here,
R ( θ ) = cos θ sin θ sin θ cos θ ,
h , H are the neutral CP-even physical d.o.f., and A and H ± are the neutral CP-odd and charged d.o.f., respectively. Z L and A are the other scalars. From Equation (5) it can be seen that β is the mixing angle of the charged and CP-odd sectors and is given by β = tan 1 ( v 2 / v 1 ) . Δ is the mixing angle of the CP-even neutral scalars and is expressed as
Δ = sin 1 M ρ 12 2 ( M ρ 12 2 ) 2 + ( M ρ 11 2 m h 2 ) 2 ,
with M ρ 2 is the mass-squared matrix in the neutral CP-even sector and the repeating and non-repeating indices refers to the diagonal and off diagonal elements, respectively. Please see Appendix B for detailed values of the masses. Here, we have assumed that h is the SM-like Higgs with mass of m h 125.09 GeV and m H > m h . It has been shown that the tree-level Higgs-mediated flavor-changing neutral currents (FCNC) appear in models where more than one scalar doublet gives mass to the same kind of SM fermions [19,20]. Such a situation can be avoided under the framework of various discrete symmetries, for example, a Z 2 symmetry [19,20].
Due to the rotation in the scalar sector following Equation (5), the couplings of the SM gauge bosons and fermions to the SM-like Higgs boson are required to be re-scaled compared to the corresponding SM values. After SSB, the Yukawa sector of the 2HDM can be written as:
L Yuk = 1 2 ( κ D s β α + ρ D c β α ) D ¯ D h + 1 2 ( κ D c β α ρ D s β α ) D ¯ D H + 1 2 ( κ U s β α + ρ U c β α ) U ¯ U h + 1 2 ( κ U c β α ρ U s β α ) U ¯ U H + 1 2 ( κ L s β α ρ L c β α ) L ¯ L h + 1 2 ( κ L c β α ρ L s β α ) L ¯ L H i 2 U ¯ γ 5 ρ U U A + i 2 D ¯ γ 5 ρ D D A + i 2 L ¯ γ 5 ρ L L A + U ¯ ( V CKM ρ D P R ρ U V CKM P L ) D H + + ν ¯ ρ L P R L H + + h . c . .
Here, κ f = 2 M f / v for f = U , D and L . L and R are defined after Equation (11). For details of other parameters see Ref. [21]. The generation indices of the fermionic fields have been suppressed in Equation (8). As mentioned earlier, the measurement of the signal strengths of the SM-like Higgs at LHC demands that the properties of one of the neutral CP-even neutral scalars, here h, should closely resemble that of the SM Higgs. As Equation (8) indicates, this is satisfied at the vicinity of the so-called ‘alignment limit’, i.e., cos ( β α ) 0 . The current measurement of Higgs signal strengths have pushed the 2HDMs close to the alignment limit [22,23,24,25]. The measurement of the Higgs signal strengths dictate that for type-II 2HDM, at tan β 1 , the constraint on cos ( β α ) is given by 0.05 cos ( β α ) 0.15 at 95 % CL. The allowed region becomes even smaller for higher values of tan β . The situation for type-III and -IV are quite similar to that of type-II 2HDM.
This constraint is comparably relaxed in type-I 2HDM, where the allowed range is | cos ( β α ) | 0.4 . Among the tree-level scalar-gauge couplings which are important for the cascade decays of the new scalars, A Z h and H ± h W are proportional to cos ( β α ) , whereas A Z H and H ± H W are proportional to sin ( β α ) . The tripple letters denote tripple coupling and Z is one of the SM guage bosons. It is possible to realize an exact alignment in the multi-Higgs-doublet models in the framework of certain additional symmetries of the 2HDM potential [26,27,28,29,30].
The impact of the measurement of the Higgs signal strengths in each individual search channels on the cos ( β α ) vs. tan β plane has been discussed in Ref. [31]. It should be mentioned that the coupling multipliers of the SM-like Higgs also becomes close to unity when sin ( β + α ) = 1 , i.e., at the so-called ‘wrong-sign Yukawa’ limit [32] for type-II, type-III, and type-IV 2HDM. Though, with better measurement of the processes like V h b b ¯ , h γ γ , Υ γ , (where b is the bottom quark), [33,34] the fate of the wrong-sign Yukawa region will be decided in near future.
It is possible to constrain the parameter space of the type-II 2HDM even at the alignment limit from the nonobservation of the heavier scalars [18] in processes like g g / b b ¯ H / A τ τ ¯ , γ γ , g g H h h , etc. Here g and τ represents gluon and tau and their respective interactions. Both ATLAS and CMS experiments are involved in numerous dedicated searches of these kinds, for instance Refs. [35,36,37,38,39], resulting in significant constraints on the 2HDM parameter space.
The constraints from the decay of the new scalars into SM particles are significantly relaxed in the hierarchical scenario compared to the degenerate case [25]. However. for the hierarchical spectrum of new scalars, the channels like H ( A ) Z A ( H ) dominates the total decay width of such states, leading to new bounds on the parameter space which are not applicable for the degenerate case. A hierarchical spectrum such as m A > m H m H ± v can lead to a first order electroweak phase transition providing an explanation for the matter-antimatter asymmetry, with A Z H being its smoking gun signature at LHC [40,41]. In general, the importance of Higgs cascade decays as the possible probes of an extended scalar sector have been discussed in the literature [42,43,44,45,46,47] and A Z H decay is dubbed as a ‘golden channel’ in this context [48].

3. EWPT Theory in 2HDM

The Lagrangian density of the electroweak theory (discussed in details in the previous section) in 2HDM can be expressed as [15]
L = L f + L Yuk + L gauge , kin + L Higgs .
The first term on the right hand side, L f , is the kinetic term for the fermion fields:
L f = j i Ψ ¯ L ( j ) D Ψ L ( j ) + Ψ ¯ R ( j ) D Ψ R ( j )
= i Ψ ¯ L γ μ ( μ + i g W μ + i g Y L B μ ) Ψ L + i Ψ ¯ R γ μ ( μ + i g W μ + i g Y R B μ ) Ψ R ,
where L and R represents the left and right chiral field of that fermion and D is the Feynman operator applied on the covariant derivative (it is defined as A = γ μ A μ , where γ are the gamma matrices) [15] and j runs over all fermionic species (the field Ψ j ) listed in Table A1. And g is the coupling constant. The partial derivative denotes derivative over space time, e.g., 0 = d / d t and μ runs from 0 to 3.
The second term of Equation (9), Yukawa interaction term (for details, see previous section), L Yuk is [49]
L Yuk = y e e R ¯ Φ a L L + y e L L ¯ Φ a e R + ,
where y e is a complex dimensionless constant, Φ a ( a = 1 , 2 ) is a S U ( 2 ) L doublet and for the Lagrangian to be gauge invariant it is coupled with another S U ( 2 ) L fermion L L . e R is the right chiral electron field and the same goes for other fermions like quarks, neutrinos, etc.
The third term L gauge , kin represents U ( 1 ) invariant kinetic term of four gauge bosons ( W i , i = 1 , 2 , 3 and B). It can be written as
L gauge , kin = 1 4 G μ ν i G i μ ν 1 4 F μ ν B F B μ ν ,
where G μ ν i = μ W ν i ν W μ i g ϵ i j k W μ j W ν k and F μ ν B = μ B ν ν B μ .
The Lagragian density for the doublet Higgs bosons is given by
L Higgs = ( D μ Φ 1 ) ( D μ Φ 2 ) + ( D μ Φ 1 ) ( D μ Φ 2 ) V tot ( Φ 1 , Φ 2 ) = { ( μ + i g T i W μ i + i g Y B μ ) Φ 1 } { ( μ + i g T i W μ i + i g Y B μ ) Φ 1 } + { ( μ + i g T i W μ i + i g Y B μ ) Φ 2 } { ( μ + i g T i W μ i + i g Y B μ ) Φ 2 } V tot ( Φ 1 , Φ 2 , T ) .
We define
W μ = g T i W μ i + g Y B μ .
Thus, from Equation (14), we get:
L Higgs , kin = ( μ Φ a ) ( μ Φ a ) i ( W μ Φ a ) ( μ Φ a ) + i ( μ Φ a ) W μ Φ a + ( W μ Φ a ) W μ Φ a .
The standard CP-conserving 2HDM potential V tot ( Φ 1 , Φ 2 , T ) consists of tree level potential V tree ( Φ 1 , Φ 2 ) :
V tree ( Φ 1 , Φ 2 ) = m 11 2 Φ 1 Φ 1 + m 22 2 Φ 2 Φ 2 m 12 2 Φ 1 Φ 2 + m 12 Φ 2 Φ 1 + 1 2 λ 1 Φ 1 Φ 1 2 + 1 2 λ 2 Φ 2 Φ 2 2 + λ 3 Φ 1 Φ 1 Φ 2 Φ 2 + λ 4 Φ 1 Φ 2 Φ 2 Φ 1 + 1 2 λ 5 Φ 1 Φ 2 2 + 1 2 λ 5 Φ 2 Φ 1 2
and other correction terms V CW v 1 , v 2 and V T . The correction terms are defined in Refs. [50,51]:
V CW v 1 + v 2 = i n i 64 π 2 ( 1 ) 2 s i m i 4 v 1 , v 2 log m i 2 v 1 , v 2 μ 2 c i ,
V T = T 4 2 π 2 i = bosons n i J B m i 2 ( v 1 , v 2 ) T 2 + i = fermions n i J F m i 2 ( v 1 , v 2 ) T 2 ,
where μ is the renormalisation scale which we take to be 246 GeV. For a detailed overview of 2HDM scenarios and other relevant aspects, see Ref. [21].
The potential-dependent mass of fermions and bosons m i v 1 + v 2 and the corresponding n i , s i , and c i are discussed in details in Appendix B.
J B and J F are approximated Landau gauge up to leading orders:
T 4 J B m 2 T = π 4 T 4 45 + π 2 12 T 2 m 2 π 6 T ( m 2 ) 3 / 2 1 32 m 4 ln m 2 a b T 2 + ,
T 4 J F m 2 T = 7 π 4 T 4 360 π 2 24 T 2 m 2 1 32 m 4 ln m 2 a f T 2 + ,
where a b = 16 a f = 16 π 2 exp ( 3 / 2 2 γ E ) with γ E being the Euler–Mascheroni constant.
When the temperature of the universe drops down to the critical temperature T c , a second local minimum appears with the same height of the global minimum situated at Φ 1 = Φ 2 = 0 [52]. The critical temperature can be obtained using the following expression:
V tot Φ 1 = 0 , Φ 2 = 0 , T c = V tot Φ 1 = v 1 , Φ 2 = v 2 , T c .
During EWPT, if ρ and P are, respectively, the energy density and pressure of the fluid determining the course of evolution of the early universe, then one gets [9]:
ρ = ρ f + ρ gauge , kin + ρ Higgs g 00 L Yuk ,
P = P f + P gauge , kin + P Higgs 1 3 g i i L Yuk .
where ρ and P are the energy and the pressure density of the plasma and the the first three terms represents in the above equations the contribution from the fermionic, gauge kinetic and the Higgs sectors, respectively. We have assumed that dark matter and other components might have been present but they did not contribute much to the energy density of the universe during the particular epoch of EWPT which happened in radiation domination. The expressions for ρ f , P f and ρ gauge , kin , P gauge , kin appear solely from fermionic and gauge sectors and their interactions (the stress–energy tensor for the above quantities is mentioned in Appendix A):
ρ H + F + G = { 0 Φ a i ( W 0 Φ a ) } 0 Φ a + { 0 Φ a + i W 0 Φ a } 0 Φ a , g 00 ( α Φ a ) ( α Φ a ) i ( W α Φ a ) ( α Φ a ) + i ( α Φ a ) W α Φ a + ( W α Φ a ) W α Φ a g 00 V tot ( Φ 1 , Φ 2 , T ) ,
P H + F + G = { q Φ a i ( W q Φ a ) } q Φ a + { q Φ a + i W q Φ a } q Φ a g q q ( α Φ a ) ( α Φ a ) i ( W α Φ a ) ( α Φ a ) + i ( α Φ a ) W α Φ a + ( W α Φ a ) W α Φ a g q q V tot ( Φ 1 , Φ 2 , T ) .
The subscript H+F+G stands for the total contribution that comes from the Higgs sector, the fermionic sector and the gauge sector. The index 00 means the first element of the metric tensor the repeated indices q q means the diagonal elements.
The early universe was flat, hence the metric g μ ν = ( + , , , ) . Furthermore, hence ρ H + F + G + P H + F + G become:
ρ H + F + G + P H + F + G = 0 Φ a 0 Φ a + 2 ( W α Φ a ) W α Φ a , + ( 0 Φ a ) ( 0 Φ a ) i ( W 0 Φ a ) ( 0 Φ a ) + i ( 0 Φ a ) W 0 Φ a
where the explicit expression for ρ H + F + G is given as
ρ H + F + G = 0 Φ a 0 Φ a + ( W α Φ a ) W α Φ a V tot ( Φ 1 , Φ 2 , T ) + L Yuk .
The oscillations of the Higgs fields around minimum after it appeared in the course of the phase transition, are damped due to particle production by the oscillating field. The characteristic time is equal to the decay width of the Higgses and it is large in comparison with the expansion rate and the universe cooling rate. So we may assume that Higgses essentially live in the minimum, of the potential. In principle, it can be calculated numerically by the solution of the corresponding Klein–Gordon equation with damping induced by the particle production [12].
With the above assumption,
ρ = Φ ˙ a , min 2 + V tot ( Φ 1 , Φ 2 , T ) + g π 2 30 T 4 .
where g is the effective degrees of freedom.
The last term in Equation (29) arises from the Yukawa interaction between fermions and Higgs bosons and from the energy density of the fermions, the gauge bosons and the interaction between the Higgs and gauge bosons. This is the energy density of the relativistic particles which have not gained mass till the moment of EWPT. g is the effective degrees of freedom. We have assumed that the energy density consists of the two parts: the energy density of the fields ϕ ( t ) s sitting at the minima of the potential and the contribution from relativistic matter.
Since for relativistic species P = ( 1 / 3 ) ρ , we can write
P = Φ ˙ a , min 2 2 + 1 3 g π 2 30 T 4 .
The oscillations of scalar fields around their minima are quickly damped, so we take the time derivative of the fields equivalent to the their derivative around the minima and neglect higher order terms of their time derivative ϕ ˙ 2 and so on. Furthermore, as a result the evolution of the minima induced by the expansion of the universe is very slow.
The entropy conservation law holds when the plasma (assumed to be an ideal fluid) was in thermal equilibrium with negligible chemical potential. However, as the temperature went below T c , EWPT happened and the universe went into a thermally nonequilibrium state. It is to be noted that one of the main consequences of EWPT is electroweak baryogenesis and, following Sakharov’s principle [53], the out-of-equilibrium process is a necessary condition for successful baryogenesis.
As a result of this deviation from thermal equilibrium, the entropy conservation law is no longer valid during EWPT and, hence, a rise in the entropy production can be noticed significantly during this process.
To calculate this production, it is necessary to solve the evolution equation for energy density conservation,
ρ ˙ = 3 H ( ρ + P ) .
From henceforward computational analysis was used for further calculations which are discussed in the next section.

4. Entropy Release in 2HDM Scenarios

In the early universe when the temperature of universe T T c , the universe was in thermal equilibrium and also was dominated by relativistic species. Almost all of the fermions and bosons were massless and the contribution from those who were already massive (e.g., DM) to the total energy density of the universe was insignificant. During that epoch, the chemical potential of the massless bosons was zero, and with the assumption that chemical potential of the fermions was negligible, the entropy density per comoving volume was conserved and given by
s ρ r + P r T a 3 = const . ,
where the subscript r is used to indicate relativistic components. For our scenario, generalising Equation (29),
ρ r + P r g T 4 ,
g 🟉 is not constant over time; it depends on the components of the hot primordial hot soup. These two equations, Equation (32) and Equation (33), imply
T a 1 .
As long as the thermal equilibrium was maintained, s remained constant. If the thermal equilibrium was not maintained at some epoch at later time, the value of s and thus g 🟉 T a 3 T 3 , might have increased as entropy can only either increase or remain constant. The contribution to entropy is dominated by a single particle species, namely by that with the largest mass in the temperature range m ( T ) < T . So for these temperatures g ( T ) T = c o n s t . and the relative entropy rise is given just by a 3 T 3 . Since the final temperature T f = m ( T f ) , below which new particle species start to dominate, does not depend upon g , the relatives entropy rise is determined by T 3 a 3 . In other words, we calculate below the entropy rise over background of the constant g a 3 T 2 .
As the universe dropped down to T c the Higgs potential gets degenerate minima. If the temperature drops to the mass of any component of the relativistic plasma, that component gained mass and became non-relativistic and decoupled from relativistic fluid. We are assuming this process was instantaneous and the universe was not in thermal equilibrium.There is change in g 🟉 of the relativistic plasma. This led to increase of s.
Suppose at T c , the scale factor was a c and g 🟉 g 🟉 , tot = 110.75 for our 2HDM model and thus at that moment s c g 🟉 , tot a c 3 T c 3 . As the temperature drops down to T T x , the component ‘x’ would decoupled and thus g 🟉 would decrease. If g 🟉 , + and g 🟉 , be the g 🟉 -factor before and after the decoupling of the ‘x’, then change in s relative to the time of critical temperature is:
δ s s c = g 🟉 , + a x T x 3 g 🟉 , tot a c T c 3 g 🟉 , tot a c T c 3 .
In what follows, the BSMPT C++ package, which deals with various properties and features related to 2HDM and baryon asymmetry is applied. The package was used to calculate the critical temperature for T c and the vacuum expectation value VEV and V eff ( T ) for each benchmark points. The calculation was repeated for 4 parameter space, the first one being the benchmark points provided in the BSMPT manual. The parameters are chosen in such a way that they satisfy the limiting conditions for type-I real 2HDM. It should also be noted that we have assumed that VEV/ T c > 1 . If this limiting assumption is varied, more benchmark points would arise, which is beyond the scope of this paper.
The differential equation Equation (31) was solved numerically by interpolating the data for V eff ( T ) for all the benchmark points and the entropy release was calculated for the same. For four different benchmark points, as mentioned in Table 1, the entropy release has been calculated with the assumption of a c T c 1 and and the results are presented in Figure 1.
As it is seen from Figure 1, the amount of entropy released increases as the critical temperature for EWPT increases. For example, the entropy production for T c = 139.5 GeV is 41 % , for T c = 151 GeV is ∼52 %, for T c = 173.5 GeV is ∼63 % and for T c = 255.5 GeV is ∼73 %. All these results are way higher than the entropy release by EWPT in the SM which is ∼13% [4].
The main reason for this excess in the production of entropy is the extra scalar bosons produced in the 2HDM which contributes the most to the process. It is to be noted that the contributions from lighter particles like electrons and neutrinos are similar to that of the SM.

5. Conclusions

In this paper, it is shown that the total release of entropy due to electroweak phase transition (EWPT) is very large even in the framework of minimal extension of the Standard Model (SM) of particle physics namely type-I two-Higgs-doublet model (2HDM) compared to the SM. It is a proven fact that unlike the SM where electroweak phase transition EWPT is of second order, in the mere extension of the SM, EWPT becomes a first order phase transition. An interesting fact is that as g , which is the effective degrees of freedom, decreases as the temperature falls down. However, as we go to a very low temperature scale, the minimum temperature takes the value of the particle mass and their contribution remains the same as that of SM.
There are two points which should be noted. First, the benchmark points are unique and they were calculated using the BSMPT package. The limiting condition of the BSMPT was set so that the vacuum expectation value (VEV) exceeds the critical temperature, T c : VEV/ T c > 1 . All the benchmark points used here satisfy this condition. Second, in this paper, we have considered only the real sector of 2HDM. If other extensions of 2HDM such as the complex 2HDM is considered, there might be considerable change in the entropy production.
In addition, one needs to mention two effects even though they are beyond the scope of this paper. First, the entropy release due to the EWPT can considerably reduce the abundance of dark matter present in the universe before EWPT. Detailed calculation of this dilution factor for SM scenario is done in Ref. [4]. Second, bubble walls that were formed might collide and may produce primordial black holes and might lead to a sufficient entropy production. The bubble collisions are also the source of primordial gravitational wave background. These will be studied in the subsequent papers.

Author Contributions

Article by A.C. and M.Y.K. The authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

The work of AC is funded by RSF Grant 19-42-02004. The work by M.K. has been performed with a support of the Ministry of Science and Higher Education of the Russian Federation, Project “Fundamental problems of cosmic rays and dark matter”, No 0723-2020-0040.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Alexander Dolgov for carefully reading the manuscript and for several helpful discussions and comments. We acknowledge private communication with Phillipp Basler regarding BSMPT. We also thank Indrani Chakraborty, Anirban Kundu, Madhumita Santra and Soubhik Mondal for discussions.

Conflicts of Interest

There has been no conflict of interest among the authors or with any published literature.

Appendix A. Energy Momentum Tensor

T f μ ν = j i Ψ ¯ L ( j ) γ μ ν Ψ ¯ L ( j ) + Ψ ¯ R ( j ) γ μ ν Ψ ¯ R ( j ) g μ ν L f ,
T gauge , kin μ ν = + F B μ α ν B α 1 4 g μ ν F α β B F B α β + G i μ α ν W α 1 4 g μ ν G α β i G i α β g ϵ i j k W μ j W α k ν W α W α j W ν k μ W α .
These are the stress-energy tensors for the fermioinc and the gauge sectors, respectively.

Appendix B. Masses of New Scalars

c i = 5 6 , i = W ± , Z , γ , 3 2 , otherwise .
m W 2 = g 2 4 v 2 ,
m Z 2 = g 2 + g 2 4 v 2 ,
m γ 2 = 0 ,
m ¯ H ± 2 = 1 2 M 11 C + M 22 C + 1 2 4 M 12 C 2 + M 13 C 2 + M 11 C M 22 C 2 ,
m ¯ G ± 2 = 1 2 M 11 C + M 22 C + 1 2 4 M 12 C 2 + M 13 C 2 + M 11 C M 22 C 2 ,
where
c 1 = 1 48 12 λ 1 + 8 λ 3 + 4 λ 4 + 3 3 g 2 + g 2 ,
c 2 = 1 48 12 λ 2 + 8 λ 3 + 4 λ 4 + 3 3 g 2 + g 2 + 24 v 2 2 m t 2 ( T = 0 ) + 1 2 v 2 2 m b 2 ( T = 0 ) ,
where m t ( T = 0 ) = 172.5 GeV and m b ( T = 0 ) = 4.92 GeV . For our case ( v 3 = 0 ) ,
M 11 C = m 11 2 + λ 1 v 1 2 2 + λ 3 v 2 2 2 ,
M 22 C = m 22 2 + λ 2 v 2 2 2 + λ 3 v 1 2 2 ,
M 12 C = v 1 v 2 2 λ 4 + λ 5 m 12 2 ,
M 13 C = 0 .
Masses of h, H and A are the eigen values of the matrix
M ¯ N = M N .
For our case ( v 3 = 0 ) ,
M 11 N = m 11 2 + 3 λ 1 2 v 1 2 + λ 3 + λ 4 2 v 2 2 + 1 2 λ 5 v 2 2 ,
M 22 N = m 11 2 + λ 1 2 v 1 2 + λ 3 + λ 4 2 v 2 2 1 2 λ 5 v 2 2 ,
M 33 N = m 22 2 + 3 λ 2 2 v 2 2 + 1 2 λ 3 + λ 4 + λ 5 v 1 2 ,
M 44 N = m 22 2 + λ 2 2 v 2 2 + 1 2 λ 3 + λ 4 λ 5 v 1 2 ,
M 12 N = 0 ,
M 13 N = m 12 2 + λ 3 + λ 4 + λ 5 v 1 v 2 ,
M 14 N = 0 ,
M 23 N = 0 ,
M 24 N = m 12 2 + λ 5 v 1 v 2 ,
M 34 N = 0 .
Table A1. Field-dependent mass of all fermions.
Table A1. Field-dependent mass of all fermions.
Fermions n i s i m f ( T = 0 )
e4 1 2 y e 2 v k lepton
μ 4 1 2 y μ 2 v k lepton
τ 4 1 2 y τ 2 v k lepton
u12 1 2 y u 2 v k quark
c12 1 2 y c 2 v k quark
t12 1 2 y t 2 v k quark
d12 1 2 y d 2 v k quark
s12 1 2 y s 2 v k quark
b12 1 2 y b 2 v k quark
Bosons n i s i m ( v ) 2
h11eigenvalues of (A15)Higgs
H11eigenvalues of (A15)Higgs
A11eigenvalues of (A15)Higgs
G 0 11eigenvalues of (A15)Goldstone
H ± 21Equation (A7)charged Higgs
G ± 21Equation (A8)charged Goldstone
Z L 11Equation (A5)Higgs
Z T 22Equation (A5)Higgs
W L 21Equation (A4)Higgs
W T 42Equation (A4)Higgs
γ L 12Equation (A6)
γ T 22Equation (A6)

References

  1. Gorbunov, D.S.; Rubakov, V.A. Introduction to the Theory of the Early Universe: Cosmological Perturbations and Inflationary Theory; World Scientific Publisher: Singapore, 2011. [Google Scholar]
  2. Bambi, C.; Dolgov, A.D. Introduction to Particle Cosmology; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar]
  3. Dolgov, A.D.; Naselsky, P.D.; Novikov, I.D. Gravitational waves, baryogenesis, and dark matter from primordial black holes. arXiv 2000, arXiv:astro-ph/0009407. [Google Scholar]
  4. Chaudhuri, A.; Dolgov, A. PBH evaporation, baryon asymmetry, and dark matter. arXiv 2020, arXiv:2001.11219. [Google Scholar]
  5. Khlopov, M. What comes after the Standard model? Prog. Part. Nucl. Phys. 2021, 116, 103824. [Google Scholar] [CrossRef]
  6. Boeckel, T.; Schettler, S.; Schaffner-Bielich, J. The Cosmological QCD phase transition revisited. Prog. Part. Nucl. Phys. 2011, 66, 266–270. [Google Scholar] [CrossRef] [Green Version]
  7. Konoplich, R.V.; Rubin, S.G.; Khlopov, M.Y. Formation of black holes in first-order phase transitions in the Universe. Astron. Lett. 1998, 24, 413–417. [Google Scholar]
  8. Jedamzik, K.; Niemeyer, J.C. Primordial black hole formation during first-order phase transitions. Phys. Rev. D 1999, 59, 124014. [Google Scholar] [CrossRef] [Green Version]
  9. Chen, C.H.; Chiang, C.W.; Nomura, T. Muon g-2 in two-Higgs-doublet model with type-II seesaw mechanism. arXiv 2021, arXiv:2104.03275. [Google Scholar]
  10. Han, X.F.; Li, T.; Wang, H.X.; Wang, L.; Zhang, Y. Lepton-specific inert two-Higgs-doublet model confronted with the new results for muon and electron g-2 anomaly and multi-lepton searches at the LHC. arXiv 2021, arXiv:2104.03227. [Google Scholar]
  11. Hernández, A.E.C. Fermion masses and mixings, dark matter, leptogenesis and g-2 muon anomaly in an extended 2HDM with inverse seesaw. arXiv 2021, arXiv:2104.02730. [Google Scholar]
  12. Chaudhuri, A.; Dolgov, A. Electroweak phase transition and entropy release in the early universe. JCAP 2018, 1801, 032. [Google Scholar] [CrossRef] [Green Version]
  13. Bochkarev, A.I.; Kuzmin, S.V.; Shaposhnikov, M.E. Electroweak baryogenesis and the Higgs boson mass problem. Phys. Lett. B 1990, 244, 275–278. [Google Scholar] [CrossRef]
  14. Basler, P.; Muhlleitner, M.; Müller, J. BSMPT v2. A Tool for the electroweak phase transition and the baryon asymmetry of the Universe in extended Higgs sectors. arXiv 2020, arXiv:2007.01725. [Google Scholar]
  15. Langacker, P. Introduction to the Standard Model and electroweak physics. In The Dawn of the LHC Era (TASI 2008); Han, T., Ed.; World Scientific Publisher: Singapore, 2010; pp. 3–48. [Google Scholar]
  16. Ginzburg, I.F.; Krawczyk, M. Symmetries of two Higgs doublet model and CP violation. Phys. Rev. D 2005, 72, 115013. [Google Scholar] [CrossRef] [Green Version]
  17. El Kaffas, A.W.; Osl, P.; Ogreid, O.M. CP violation, stability and unitarity of the two Higgs doublet model. Nonlin. Phenom. Complex Syst. 2007, 10, 347. [Google Scholar]
  18. Grzadkowski, B.; Haber, H.E.; Ogreid, O.M.; Osl, P. Heavy Higgs boson decays in the alignment limit of the 2HDM. J. High Energy Phys. 2018, 12, 056. [Google Scholar] [CrossRef] [Green Version]
  19. Glashow, S.L.; Weinberg, S. Natural conservation laws for neutral currents. Phys. Rev. D 1977, 15, 1958. [Google Scholar] [CrossRef]
  20. Paschos, E.A. Diagonal Neutral Currents. Phys. Rev. D 1977, 15, 1966. [Google Scholar] [CrossRef]
  21. Branco, G.C.; Ferreira, P.M.; Lavoura, L.; Rebelo, M.N.; Sher, M.; Silva, J.P. Theory and phenomenology of two-Higgs-doublet models. Phys. Rep. 2012, 516. [Google Scholar] [CrossRef] [Green Version]
  22. Gu, J. Learning from Higgs Physics at Future Higgs Factories. J. High Energy Phys. 2017, 1712, 153. [Google Scholar] [CrossRef] [Green Version]
  23. Bernon, J. Scrutinizing the alignment limit in two-Higgs-doublet models. II. mh = 125 GeV. Phys. Rev. D 2016, 93, 035027. [Google Scholar] [CrossRef] [Green Version]
  24. Bernon, J. Scrutinizing the alignment limit in two-Higgs-doublet models: mh = 125 GeV. Phys. Rev. D 2015, 92, 075004. [Google Scholar] [CrossRef] [Green Version]
  25. Dorsch, G.C.; Huber, J.S.; Mimasu, K.; No, J.M. Hierarchical versus degenerate 2HDM: The LHC run 1 legacy at the onset of run 2. Phys. Rev. D 2016, 93, 115033. [Google Scholar] [CrossRef] [Green Version]
  26. Bhupal Dev, P.S.; Pilaftsis, A. Erratum to: Maximally Symmetric Two Higgs Doublet Model with Natural Standard Model Alignment. J. High Energy Phys. 2015, 1511, 147. First published: J. High Energy Phys. 2014, 1412, 24. [Google Scholar] [CrossRef] [Green Version]
  27. Pilaftsis, A. Symmetries for standard model alignment in multi-Higgs doublet models. Phys. Rev. D 2016, 93, 075012. [Google Scholar] [CrossRef] [Green Version]
  28. Das, D.; Dey, U.K.; Pal, P.B. Quark mixing in an S3 symmetric model with two Higgs doublets. Phys. Rev. D 2017, 96, 031701. [Google Scholar] [CrossRef] [Green Version]
  29. Benakli, K.; Chen, Y.; Lafforgue-Marmet, G. R-symmetry for Higgs alignment without decoupling. Eur. Phys. J. C 2019, 79, 172. [Google Scholar] [CrossRef]
  30. Pramanick, S.; Raychaudhuri, A. Three-Higgs-doublet model under A4 symmetry implies alignment. J. High Energy Phys. 2018, 1801, 11. [Google Scholar] [CrossRef] [Green Version]
  31. Craig, N.; Galloway, J.; Thomas, S. Searching for signs of the second Higgs doublet. arXiv 2013, arXiv:1305.2424. [Google Scholar]
  32. Ferreira, P.M. The Wrong Sign limit in the 2HDM. arXiv 2014, arXiv:1410.1926. [Google Scholar]
  33. Modak, T. Constraining wrong-sign hbb couplings with hΥγ. Phys. Rev. D 2016, 94, 075017. [Google Scholar] [CrossRef] [Green Version]
  34. Ferreira, P.M. Probing wrong-sign Yukawa couplings at the LHC and a future linear collider. Phys. Rev. D 2014, 89, 115003. [Google Scholar] [CrossRef] [Green Version]
  35. Sirunyan, A.M. [CMS Collaboration] Search for Higgs boson pair production in the γ γ b b ¯ final state in pp collisions at s = 13 TeV. arXiv 2018, arXiv:1806.00408. [Google Scholar]
  36. The ATLAS Collaboration. Search for a CP-odd Higgs Boson Decaying to Zh in pp Collisions at s = 13 TeV with the ATLAS Detector. ATLAS-CONF-2016-015; CERN: Geneva, Swizerland, 2016. Available online: https://cds.cern.ch/record/2141003/files/ATLAS-CONF-2016-015.pdf (accessed on 27 April 2021).
  37. The ATLAS Collaboration. Search for Additional Heavy Neutral Higgs and Gauge Bosons in the Ditau Final State Produced in 36.1 fb−1 of pp Collisions at s = 13 TeV with the ATLAS Detector. ATLAS-CONF-2017-050; CERN: Geneva, Swizerland, 2017. Available online: http://cds.cern.ch/record/2273866/files/ATLAS-CONF-2017-050.pdf (accessed on 27 April 2021).
  38. The CMS Collaboration. Search for Additional Neutral Higgs Bosons Decaying to a Pair of Tau Leptons in pp Collisions at s = 7 and 8 TeV. CMS-PAS-HIG-14-029. CERN: Geneva, Swizerland, 2015. Available online: https://cds.cern.ch/record/2041463/files/HIG-14-029-pas.pdf (accessed on 27 April 2021).
  39. The ATLAS Collaboration. Search for Heavy Resonances Decaying to a W or Z Boson and a Higgs Boson in Final States with Leptons and b-Jets in 36.1~fb−1 of pp Collision Data at s = 13 TeV with the ATLAS Detector. ATLAS-CONF-2017-055; CERN: Geneva, Switzerland, 2017. Available online: https://cds.cern.ch/record/2273871/files/ATLAS-CONF-2017-055.pdf (accessed on 27 April 2021).
  40. Dorsch, G.C.; Huber, S.J.; Mimasu, K.; No, J.M. Echoes of the electroweak phase transition: Discovering a second Higgs doublet through A0ZH0. Phys. Rev. Lett. 2014, 113, 211802. [Google Scholar] [CrossRef] [Green Version]
  41. Dorsch, G.C.; Huber, S.J.; No, J.M. A strong electroweak phase transition in the 2HDM after LHC8. J. High Energy Phys. 2013, 1310, 29. [Google Scholar] [CrossRef] [Green Version]
  42. Gao, C.; Luty, M.A.; Mulhearn, M.; Neill, N.A.; Wang, Z. Searching for additional Higgs bosons via Higgs cascades. Phys. Rev. D 2018, 97, 075040. [Google Scholar] [CrossRef] [Green Version]
  43. Bauer, M.; Carena, M.; Gemmler, K. Flavor from the electroweak scale. J. High Energy Phys. 2015, 1511, 016. [Google Scholar] [CrossRef] [Green Version]
  44. Kling, F.; Pyarelal, A.; Su, S. Light charged Higgs bosons to AW/HW via top decay. J. High Energy Phys. 2015, 1511, 051. [Google Scholar] [CrossRef] [Green Version]
  45. Coleppa, B.; Kling, F.; Su, S. Charged Higgs search via AW±/HW± channel. J. High Energy Phys. 2014, 1412, 148. [Google Scholar] [CrossRef]
  46. Kling, F.; No, J.M.; Su, S. Anatomy of exotic Higgs decays in 2HDM. J. High Energy Phys. 2016, 1609, 093. [Google Scholar] [CrossRef] [Green Version]
  47. Adhikary, A.; Banerjee, S.; Barman, R.K.; Bhattacherjee, B. Resonant heavy Higgs searches at the HL-LHC. J. High Energy Phys. 2019, 9, 068. [Google Scholar] [CrossRef] [Green Version]
  48. Kling, F.; Li, H.; Pyarelal, A.; Song, H.; Su, S. Exotic Higgs decays in type-II 2HDMs at the LHC and future 100 TeV hadron colliders. J. High Energy Phys. 2019, 6, 31. [Google Scholar] [CrossRef] [Green Version]
  49. Logan, H.E. TASI 2013 Lectures on Higgs physics within and beyond the Standard Model. arXiv 2013, arXiv:1406.1786. [Google Scholar]
  50. Blinov, N.; Profumo, S.; Stefaniak, T. The electroweak phase transition in the inert doublet model. JCAP 2015, 7, 028. [Google Scholar] [CrossRef]
  51. Basler, P.; Krause, M.; Muhlleitner, M.; Wittbrodt, J.; Wlotzka, A. Strong first order electroweak phase transition in the CP-conserving 2HDM revisited. J. High Energy Phys. 2017, 2, 121. [Google Scholar] [CrossRef] [Green Version]
  52. Katz, A.; Perelstein, M. Higgs couplings and electroweak phase transition. J. High Energy Phys. 2014, 7, 108. [Google Scholar] [CrossRef] [Green Version]
  53. Sakharov, A.D. Violation of CP Invariance, C Asymmetry, and Baryon Asymmetry of the Universe. J. Exp. Theor. Phys. Lett. 1967, 5, 24–27. [Google Scholar]
Figure 1. The entropy release for four different scenarios are shown above. The left panel shows the entropy production for T c = 139.5 GeV and 151 GeV (Benchmark-IV and III, respectively). The right panel shows the entropy production for T c = 173.5 GeV (blue line, Benchmark-II) and the bottom panel shows the entropy production for T c = 255.5 GeV (black line, Benchmark-I).
Figure 1. The entropy release for four different scenarios are shown above. The left panel shows the entropy production for T c = 139.5 GeV and 151 GeV (Benchmark-IV and III, respectively). The right panel shows the entropy production for T c = 173.5 GeV (blue line, Benchmark-II) and the bottom panel shows the entropy production for T c = 255.5 GeV (black line, Benchmark-I).
Physics 03 00020 g001
Table 1. 2HDM benchmark points for entropy production.
Table 1. 2HDM benchmark points for entropy production.
m h [GeV] m H [GeV] m H ± = m A [GeV] tan β
Benchmark-I12550050010
Benchmark-II12540050010
Benchmark-III1259040010
Benchmark-IV125.09228.172336.94
λ 1 λ 2 λ 3 λ 4 λ 5 T c [GeV]
Benchmark-I4.130.294.1500255.5
Benchmark-II0.250.2512.65−1.48−1.48173.5
Benchmark-III0.1330.2595.02−2.51−2.51151
Benchmark-IV1.220.29−0.514.07−3.86139.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chaudhuri, A.; Khlopov, M.Y. Entropy Production Due to Electroweak Phase Transition in the Framework of Two Higgs Doublet Model. Physics 2021, 3, 275-289. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3020020

AMA Style

Chaudhuri A, Khlopov MY. Entropy Production Due to Electroweak Phase Transition in the Framework of Two Higgs Doublet Model. Physics. 2021; 3(2):275-289. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3020020

Chicago/Turabian Style

Chaudhuri, Arnab, and Maxim Yu. Khlopov. 2021. "Entropy Production Due to Electroweak Phase Transition in the Framework of Two Higgs Doublet Model" Physics 3, no. 2: 275-289. https://0-doi-org.brum.beds.ac.uk/10.3390/physics3020020

Article Metrics

Back to TopTop