SciELO - Scientific Electronic Library Online

 
vol.62 número4Estudio teórico sobre el cambio de índice de refracción y la absorción óptica en un punto cuántico en presencia de un campo magnético uniformeDiffusion of chemically reactive species in stagnation point flow of a micropolar viscoelastic fluid in a porous medium over a stretching/shrinking sheet índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.62 no.4 México jul./ago. 2016

 

Research

Analytical calculation of radiative corrections of a THDM potential

E. Díaz Méndez* 

1*Facultad de Ciencias Físico-Matemáticas, Benemérita Universidad Autónoma de Puebla, 72570, Puebla, Pue., México.


Abstract:

We obtain a closed form effective potential at the one-loop level of a Two Higgs Doublet Model. Through the loop expansion we reproduce the expression presented by Weinberg and Coleman, showing explicitly every step involved in the calculation. The formalism is then extended to include interaction terms between two scalar doublets and a final expression of the one-loop level contributions is presented.

Keywords: Two Higgs doublet model; effective potential

PACS: 12.15.Lk; 12.60.Fr

1. Introduction

In the Standard Model (SM) spontaneous symmetry breaking is performed making use of tree level terms in the potential, but this is only an approximate method in the full quantum theory, since it disregards the quantum corrections coming from virtual processes at the loop level, which are particularly important when considering zero and finite temperature studies. In Ref. 1 a first effort to extend the potential considered taking into account the one-loop level corrections, in order to study the phase transition between SSB and non-SSB scenarios for different models: the authors used explicitly the Landau gauge for the scalar electrodynamics theory. Later on, Jackiw 2 made the calculation in an arbitrary gauge up to two-loop level for a set of scalars with an O(n) internal symmetry. Subsequently, Arnold et al. 3 extended the result for the whole SM with a single Higgs doublet including the thermal contributions.

However, as far as we know, no one has determined from first principles the one-loop level vacuum contributions in which more than one scalar field are present and have interaction terms, such as in the Two Higgs Doublet Model (THDM). Although different models have been presented in which several non interacting scalar fields are present. In the cosmological arena the scalar sector effective potential is useful to study the electroweak phase transition ,4 on which both vacuum and finite-temperature contributions must be taken into account.

In a recent paper 5 a first attempt was made to have an effective potential in which two scalar fields have coupling terms present in the potential. Unfortunately the way in which multiple scalar potentials have been constructed makes use of a generic expression coming from the single scalar one-loop contribution, without an explicit and detailed derivation. In the present work, we obtain an effective potential for the scalar sector of a generic THDM. We first perform a detailed analysis for the single massless and massive field theories and then work out the THDM extension, we do this only for the scalar sector, since all other sectors are equal to those of the Minimal Standard Model, for which the effective potential has been obtained even up to the three loop level in some scenarios 6.

The rest of our paper is structured as follows: In Sec. 2 we present the general features used in the derivation of the effective potential. And then proceed to calculate the massless self interacting real scalar field potential in Sec. 3. We then work out the case of a massive self interacting real scalar field potential in Sec. 4, and in Sec. 5 we obtain the expression for the THDM effective potential. Our conclusions are presented in Sec. 6.

2. Derivation of the Effective Potential

For a single real massless self interactive scalar, whose Lagrangian is given by

L=12(μϕ)2-λ4!ϕ4, (1)

the one loop level corrections to the potential are given by an infinite series of Feynman diagrams; the first four of which are shown in Fig. 1.

FIGURE 1 First four one-loop Feynman diagrams that contribute to the effective potential. 

The first diagram contains only one vertex, so we only need one factor of λ (the factor of i cancels with the one of the propagator); and it has a single propagator inside the loop. The second diagram, has two vertices and two propagators. In general, the n-th diagram contains n vertices and n propagators, each one corresponding to a side of the polygon formed in the center of each diagram.

So, from the first diagram we obtain:

λ1k2+iϵ, (2)

where k represents the momentum of the particle going around the loop. It is a virtual particle, so we actually have to integrate over all of its possible values, that is, from minus infinity to plus infinity (for each component):

d4k(2π)4λk2+iϵ. (3)

The denominator (2π)4 is conventional, and the integration limits are usually left implicit when they are ±∞. Now we must add combinatoric factors: the diagram stays the same if we interchange the two external legs attached to the vertex, so we add a 1/2; and it is also invariant under a reflection through a vertical axis passing through the center of the loop, so we have another 1/2:

d4k(2π)41212λk2+iϵ. (4)

Finally, because of the definition of the connected generating functional, we have a factor of i and because of the form of the Taylor series expansion, we must add a ɸ 2 for each vertex in the diagram, so the final expression for the first diagram is

id4k(2π)41212λϕ2k2+iϵ, (5)

where the ɸ is defined as the functional derivative of the connected generating functional with respect to the external source.

Repeating the above reasoning, we arrive at the conclusion that the n-th diagram is represented by the expression

id4k(2π)412n12λϕ2k2+iϵn. (6)

Summing all of them, we obtain

id4k(2π)4n=112n12λϕ2k2+iϵn. (7)

Performing a Wick rotation, that is, making the change of variable k0 ik0 of the time component of the momentum; we eliminate the pre-factor of i and leave a minus sign. We also identify the infinite sum with the following Taylor series expansion of the logarithm:

ln(1-x)=-n=1xnn, (8)

and dropping the ,i€

id4k(2π)4n=112n(12λϕ2k2+iϵ)n=12d4k(2π)4ln(1+λϕ22k2)=Λ464π2ln(1+λϕ22Λ2)+λϕ2Λ2128π2-λ2ϕ4256π2ln(1+2Λ2λϕ2). (9)

In general, we set the second derivative of the potential, evaluated at ɸ=0, equal to the (squared) renormalized mass μ of the particle:

d2Vdϕ2|ϕ=0=μ2. (10)

and for the the particular case of a massless scalar (which we will treat in the following section), we have

d2Vdϕ2|ϕ=0=0. (11)

Another renormalization condition that is employed concerns the fourth derivative of the potential which is equated to the coupling constant λ

d4Vdϕ4|ϕ=M=λ, (12)

where we evaluate ɸ at M to avoid IR divergences. M is called the renormalization mass which is an arbitrary quantity.

3. Massless Self Interacting Real Scalar Field

The Lagrangian for a real massless self interacting scalar field model is of the form

L=12(μϕ)2+12A(μϕ)2-λ4!ϕ4-12Bϕ2-14!Cϕ4, (13)

where the terms containing A, B, and C are the counterterms. Hence the one-loop-level effective potential is

V(ϕ)=λ4!ϕ4+12Bϕ2+14!Cϕ4+id4k(2π)4n=112n12λϕ2k2+iϵn, (14)

the integral over the sum of the diagrams can be expressed as the logarithm

12d4k(2π)4ln1+λϕ22k2=Λ464π2ln1+λϕ22Λ2+λϕ2Λ2128π2-λ2ϕ4256π2ln1+2Λ2λϕ2. (15)

The detailed derivation of the integral is found in Appendix A. We can drop the first term because it vanishes in the limit of large Λ; also in this limit we have

ln1+2Λ2λϕ2=-lnλϕ22Λ2, (16)

substituting the analytical solution of the logarithmic integral that we solve we now arrive at the following expression for the potential

V(ϕ)=λ4!ϕ4+12Bϕ2+14!Cϕ4+λϕ2Λ2128π2+λ2ϕ4256π2ln(λϕ22Λ2), (17)

the counterterm coefficient B is obtained through Eq. 11

B=-λΛ264π2, (18)

and the counterterm coefficient C, from the condition given by Eq. 12

C=-λ225+6lnλM22Λ264π2, (19)

after substituting them into the effective potential, and a little bit of work, we get

V(ϕ)=λ4!ϕ4+λ2ϕ4256lnϕ2M2-256. (20)

It is important to notice that although the coefficients of the counter terms B and C that we obtained are different to those obtained in Ref. 1, our final expression for the potential is the same.

4. Massive Self Interacting Real Scalar Field

The Lagrangian for a massive self interacting scalar field theory is given by

L=12(μϕ)2+12A(μϕ)2-μ22ϕ2-λ4!ϕ4-12Bϕ2-14!Cϕ4, (21)

so we identify the potential as

V(ϕ)=μ22ϕ2+λ4!ϕ4+12Bϕ2+14!Cϕ4+id4k(2π)4n=112n12λϕ2k2-μ2+iϵn. (22)

Once again we are left with an integral over the sum of the one loop level self interacting contributions which we write as a logarithmic integral (the detailed derivation is found in Appendix B)

12d4k(2π)4ln1+λϕ22(k2+μ2-iϵ)=-2μ4lnμ2Λ2+12(2μ2+λϕ2)2ln2μ2+λϕ22Λ2+λΛ2ϕ2, (23)

from which we arrive at the effective one loop level potential for a massive interacting scalar field

V(ϕ)=μ22ϕ2+λ4!ϕ4+164π2μ2+λϕ222ln1+λϕ22μ2-12λμ2ϕ2-2524λ2ϕ4+14λ2ϕ4ln2μ2λM2. (24)

The methodology used in obtaining the analytical solution for this case will now be extended to calculate the potential of interacting scalar fields in the THDM.

5. Effective Potential for a Two Higgs Doublet Model

Because we are interested in the one-loop vacuum corrections of the scalar sector of the THDM (a review on the THMD can be found in Ref. 7), we will only work with the contributions arising from the scalar potential

V(Φ1,Φ2)=m12Φ1Φ1+m22Φ2Φ2-m122Φ1Φ2+Φ2Φ1+λ1Φ1Φ12+λ2Φ2Φ22+λ3Φ1Φ1Φ2Φ2+λ4Φ1Φ2Φ2Φ1+λ5Φ1Φ22+Φ2Φ12, (25)

where the Higgs doublets can be written as

Φ1=ϕ1+ϕ10,Φ2=ϕ2+ϕ20, (26)

but the physical fields of the charged sector are defined through the transformation

G±H±=cccosβsinβ-sinβcosβϕ1±ϕ2±, (27)

while the imaginary terms are transformed via

H0h0=cccosαsinα-sinαcosαReϕ10Reϕ20, (28)

and the real neutral sector is obtained with the transformation

G0A0=cccosβsinβ-sinβcosβImϕ10Imϕ20. (29)

Through the previously defined transformations we are able to obtain all of the self interaction vertices between each of the fields h0 ,H0 , A0 , H± , G0 and G± found in Ref. 8. But since we choose to work in the unitary gauge, the Goldstone bosons will be “eaten” by the SU(2) x U(1) gauge fields. Therefore, we will only be dealing with the physical fields.

Furthermore, we only consider quartic vertices. Diagrams containing triple vertices such as those shown in Fig.2 give loop corrections to the propagator, but they are not relevant to our present work, so we do not consider them. The effective potential in integral form is given in Appendix C, from which we obtain the effective potential as a function of the physical fields:

FIGURE 2 Scalar propagator correction diagrams. 

V(Φ1,Φ2)=Vtree+VH0H0H0H0+Vh0h0h0h0+VA0A0A0A0+VH0h0h0h0+Vh0H0H0H0+VH0H0h0h0+VH0H0A0A0+Vh0h0A0A0+VH+H+H-H-+Vh0h0H-H++VH0H0H-H++VA0A0H-H++VA0A0H0h0+VH0h0H-H+, (30)

the first three terms corresponding to neutral Higgs of the same type VSSSS where S=h0, H0, A0 are given by

VSSSS=-λμS2S2-2512λ2S4+λ2S4(16λμS2M2-64μS4)12(λM2+2μS2)2+2λμS2S2ln1+λS22μS2+2μS4ln1+λS22μS2+12λ2S4lnμS2+λS22μS2+λM22, (31)

next we have the coupling of three fields of the same type (S1 ) with another one of a different kind (S2 ) given by

VS1S1S1S2=-λμS12S2S11-2ln1+λS2S12μS12+λ24S22S122ln1+λS2S12μS12-3+2μS14ln1+λS2S12μS12, (32)

during the calculation for this case, we set the one of the mass counter term coefficients B=0, because if one includes all the mass counter terms you end up with a divergent logarithm.

The subsequent four terms of the effective potential which have two equal pairs S1, S2 are given by

VS1S1S2S2=-λμS2(S12+S22)-34λ2(S14+S24)+4(S14+S24)M2λ33(M2λ+2μS2)2(μS2-M2λ)+2λμS2S12ln1+λS122μS2+2λμS2S22ln1+λS222μS2+λ22S14lnμS2+λS122μS2+λM22+λ22S24lnμS2+λS222μS2+λM22. (33)

Finally we have the terms in which we have a single pair of scalars of the same kind and the remaining two are of a different type

VS1S1S2S3=-λμS12S2S31-2ln1+λS2S32μS12+λ24S22S322ln1+λS2S32μS12-3+2μS14ln1+λS2S32μS12. (34)

the last term in the potential involves the coupling of fields which are all of a different kind, such as H0h0H+H- , this type of coupling does not permit a polygon diagram whose series is constructible.

6. Conclusions

We calculated the contributions of the THDM effective potential that come from one-loop diagrams with quartic coupling vertices which are dimensionally consistent and renormalizable. We hope this analytical scheme will be of use to model builders interested in theories that must take in to consideration vacuum corrections to two or more scalars.

Acknowledgments

We acknowledge support from CONACYT (México), and appreciate the advice and discussions from Gerardo F. Torres del Castillo and O. Meza-Aldama.

References

1. S.R. Coleman, and E. Weinberg, J. Phys. Rev. D 7 (1973) 1888-1910. [ Links ]

2. R. Jackiw, Phys. Rev. D 9 (1974) 1686. [ Links ]

3. P.B. Arnold and O. Espinosa, Phys. Rev. D 47 (1993) 3546. Erratum-ibid. D 50 (1994) 6662. arXiv:hep-ph/9212235. [ Links ]

4. R.R. Caldwell, and S.S. Gubser, 87 (2013) 063523. arXiv:1302.1201. [ Links ]

5. Indrani Chakraborty, and Anirban Kundu, (2015), arXiv: 1508.00702. [ Links ]

6. P. Stephen Martin, Phys. Rev., D 89 (2014) 013003. arXiv: 1310.7553. [ Links ]

7. G.C. Branco, P.M. Ferreira, L. Lavoura, M.N. Rebelo, M. Sher, and J.P. Silva, Phys. Rept. 516 (2012) 1. arXiv:1106.0034. [ Links ]

8. J. F. Gunion, H. E. Haber, G. L. Kane and S. Dawson, Front. Phys. 80 (2000) 1. [ Links ]

Appendix A. Solution to the integral of the massless effective potential

We proceed to solve the integral by first changing from an integral over 4-momentum to an integral over hyper spherical coordinates

k0=ρcosψ,k1=ρsinψcosθ,k2=ρsinψsinθcosϕ,k3=ρsinψsinθsinϕ,d4k=ρ3sin2ψsinθdρdψdθdϕ,

Where θ ,ψ (0,π) and ɸ € (0,2π).

From which we obtain

12d4k(2π)4ln(1+λϕ22k2)=12(2π)40πsin2ψdψ×0πsinθdθ02πdϕ0Λρ3ln(1+λϕ22ρ2)dρ. (A.1)

We now solve the integral by

(A.2)

the first term cancels because when Λ→∞ the logarithm goes to zero, as for the last term

I1=λϕ232π20Λρ3dρ2ρ2+λϕ2, (A.3)

making x =2ρ2dx =4ρdρxdx =8ρ3dρ,

(A.4)

Now defining y=x+λϕ2 dy=dx

(A.5)

(A.6)

In the large cutoff limit, Λ1, we have

ln(1+2Λ2λϕ2)=-ln(λϕ22Λ2), (A.7)

so the effective potential is

(A.8)

B. Solution to the integral of the massive effective potential

In this case we have that the sum of 1-loop diagrams is given by the following integral:

I12d4k(2π)4ln(1+λϕ22(k2+μ2-iϵ))=12(2π2)1(2π)4ρ3ln(1+λϕ22(ρ2+μ2))dρ, (B.1)

we now make the change of variable x=2ρ2 dx =4ρdρ xdx =8ρ3dρ. Then dropping the iϵ,

I=116π218xln(1+λϕ2x+μ2)dx=116π218ln(1+λϕ2x+μ2)d(12x2)=1128π212x2ln(1+λϕ2x+2μ2)-1128π212x211+λϕ2x+2μ2-λϕ2(x+2μ2)dx=1256π2x2ln(1+λϕ2x+2μ2)+λϕ2256π2x2dx(x+2μ2+λϕ2)(x+2μ2). (B.2)

Let’s now solve the rightmost integral (redefining some x-independent quantities as a and b):

I1=x2dx(x+a)(x+b), (B.3)

setting u = x2 and dv=dx(x+a)(x+b), to get v we integrate:

(B.4)

We now use the following known integral

(B.5)

where in our case B=1, A=-a+b, and u=y, thus

(B.6)

For each term on the right we have

(B.7)

(B.8)

putting both together, and doing some algebra,

I1=b2a-bln(x+b)+a2b-aln(x+a)+x+32(b+a).

Finally, placing a cutoff at k 22, or equivalently x 2 = 2Λ2,

(B.9)

We neglect the first term since when Λ→∞ it tends to zero. Then, evaluating, we’re left only with

(B.10)

In the limit of large Λ, we get

(B.11)

C. THDM intergral form effective Potential

There are essentially three types of vertices involved in the calculation: those with four equal fields, with two pair of equal fields, and with only one pair of equal fields. We only consider the case in which only one particle goes around the loop. Thus, when we analyze a vertex with more than one particle interacting, we write only the propagators of the particle that is “repeated”. In the first case, the integral form of the potential is

VSSSS=Vtree+d4k(2π)4ln(1+λS22(k2+μS2)), (C.1)

where λ represents genericallly the self-coupling constant and μs is the mass of the particle S, and we’ve omitted the i iϵ. In the second case,

VS1S1S2S2=Vtree+d4k(2π)4ln(1+λS122(k2+μS22)), (C.2)

where there’s an obvious symmetry between S1 and S2 . Finally, for the last case,

VS1S1S2S3=Vtree+d4k(2π)4ln(1+λS2S32(k2+μS12)). (C.3)

Received: February 24, 2016; Accepted: March 29, 2016

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License