sign in sign up
<<
Home , Baryon asymmetry

arXiv:1806.00597v2 [hep-ph] 9 Sep 2018 atta ed ob xlie nalcsooia oes hsasymm This models. ( density cosmological all or in density explained to be density to number needs that fact W where where oe feeetr atce l nw n eie Pvoaini enc mixin is quark violation the to CP . the referring verified in processes and relevant known the significan all that special particles of expects which elementary among fie of [1], and Sakharov model particles by of proposed theory first quantum as the This on . based electroweak scale is electroweak asymmetry explaining for aynccag is: baryonic n[] hsfia siaesoswtotdutta ti mosbei h s the in impossible is experimen the it by that requested analysis doubt asymmetry comprehensive without baryon and shows desired detailed estimate the A final reach [7]. This in [8]. disputed in further was [6] in h ueia eut hra eto sddctdt Conclusions to dedicated poten effective is V thermal section calculation. loop whereas our one results the for d numerical of set-up has the framework the that the contains and in II [8] analytically in Section estimates theory. quantu theoretical simple model the a with present agreement will we good Here baryogenesis. electroweak the † ˜ mi:[email protected] Email: ic h nvrei otycmoe fmte,tematter-antim the , of composed mostly is Universe the Since h iegneo h ayncreti h tnadmdli give is model standard the in current baryon the of divergence The h acltoso h aynaymtyaecmlxadintricate and complex are asymmetry baryon the of calculations The baryon the computing for methods developed authors the [2]-[6] In aµν sisda.Here, dual. its is n q f saeteqak n hr sa mle u vrflvr.Te ythe by Then flavors. over sum implied an is there and quarks the are ’s stenme ffamilies, of number the is a ainlIsiueo hsc n ula niern OBo PO Engineering Nuclear and of Institute National ihteol oreo Pvoaini h K arx u res lit Our c the standard in matrix. in found CKM asymmetry baryon the asymmetry required in baryon experimentally violation for the CP limits of theoretical source the only the with xeddt eodsadr oe theories. model standard beyond to extended ecmuetebro smer nteUies nteframew the in Universe the in asymmetry baryon the compute We aynaymtyi h tnadmdlrevisited model standard the in asymmetry Baryon ∂ µ F j µν B µ = steeetoantctensor, electromagnetic the is n f Dtd etme 1 2018) 11, September (Dated: " B .INTRODUCTION I. 32 g ( I H SET-UP THE II. t π 2 eaaJora Renata = ) 2 j B W µ Z µν = a η W d 1 2 ˜ 3 = qγ ¯ aµν xj µ n B 0 b q, n − ( − γ ,t x, a n ¯ em r fmxmmiprac o h baryon the for importance maximum of are terms g 32 a data. tal b d.Te h hoyrqie he anproperties main three requires theory the Then lds. g † ) nscinIItebro smer scomputed is asymmetry baryon the III section In ei h rsneo Pvoain ntestandard the In violation. CP of presence the is ce edter ehdta ed orslsi very in results to leads that method theory field m or ′ π . rastegnrto fbro smer tthe at asymmetry baryon of generation the treats 2 . rc plcblt oaybyn ftestandard the of beyond any to applicability irect 2 soois h ehdcnb easily be can method The osmologies. ilo h tnadmdl eto Vincludes IV Section model. standard the of tial rtr n o ml oacutfor account to small too and erature y[9]: by n F η n r ae ncranasmtosregarding assumptions certain on based are and G6 uhrs-auee Romania Bucharest-Magurele, MG-6, x smer ntesadr oe.Teresult The model. standard the in asymmetry F ˜ µν te smer sa xeietlybased experimentally an is asymmetry atter = µν F ntesadr oe a hnperformed then was model standard the in ˜ puae nteCMmti.Tu one Thus matrix. CKM the in apsulated n sisdual, its is µν tyi esrdb h ai fbaryon of ratio the by measured is etry l swti n eycoeto close very and within is ult b − ml otx ftesadr oe to model standard the of context imple s # n ¯ r ftesadr model standard the of ork , b .Oeo h otrlal theories reliable most the of One ). ote hoe h conserved the theorem Noether W µν a sthe is SU (2) L esrand tensor (1) (3) (2) 2

The change in is then,

tf tf 3 0 3 µ B(tf ) − B(0) = dt d x∂0jB = dt d x∂µjB . (4) Z0 Z Z0 Z Since the electromagnetic field in Eq. (1) turns out irrelevant one may also write:

△B = B(tf ) − B(0) = nf [NCs(tf ) − NCs(0)] = g2 tf n dt d3xW a W˜ aµν . (5) f 32π2 µν Z0 Z a Here NCS is the Chern-Simons number associated to the field Wµ . Calculating the expression on the right hand side of Eq. (5) is a non-trivial task. In what follow we will consider a simpler and equivalent approach. The quantum associated to the baryon divergence of current is:

µ a ¯ µ 4 h∂µjBi = dAµdΨdΨdh∂µjB exp[i d xL], (6) Z Z ¯ a where Ψ and Ψ stand for all the fermion species, Aµ for all gauge bosons and h for all scalars in the theory. Eq. (6) may be rewritten as:

µ h∂µjB i = 1 ∂S ∂S dAa dΨ¯ dΨdh[iΨ¯ + i Ψ] exp[i d4xL]= 2 µ ∂Ψ¯ ∂Ψ Z Z 1 = 48δ(0) dAa dΨ¯ dΨdh exp[i d4x[L]. (7) 2 µ Z Z Here the equation of motion for the fermion fields and also the Dyson Schwinger equations were applied. Moreover we summed up over the fermions flavors. If one considers an approach in which the effective Higgs potential is determined at some relevant temperature then Eq. (4) together with Eq. (7) lead to:

B(tf ) − B(0) ≈ 24Vtf δ(0) exp[iΓ(T )], (8) where Γ(T ) is the effective action at the temperature T . Here V is the space volume and tf is the time interval. The baryon asymmetry in the Universe is defined as [9]: n η = B , (9) nγ where nB = nb − n¯b is the difference between the baryon and antibaryons numbers per unite volume and nγ is the photon density. One can further write Eq. (9) as [10]:

△BBnγ η ≈ ≈△BB, (10) nγ 1 ¯ where this time △BB = 2 (△B −△B). According to our result in Eq. (8) then: 1 △ B = (△B −△B¯)=12Vt δ(0)[exp[iΓ(T )] − exp[iΓ(¯ T )], (11) B 2 f where Γ(¯ T ) is the CP conjugate effective action at temperature T . It is more convenient to use the following definition for the baryon asymmetry:

nB △BB η1 = ≈ , (12) s g∗

∗ where s ≈ nγ g is the entropy density and g∗ = 106.75 is the number of degrees of freedom in the [11]. Our next task is to determine the difference that appears on the right hand side of Eq. (11). 3

III. ASYMMETRY OF THE EXPONENTIAL OF THE EFFECTIVE ACTION

The effective action of the standard model at finite temperature is obtained from the regular partition function in 1 terms of the classical Higgs field calculated for a finite time of order T [12]. Let us start with the definition of interest and leave aside for the moment the temperature:

a ¯ 4 exp[iΓ] = dAµdhdΨdΨ exp[i d xL(Φ)]. (13) Z Z Since we are actually interested in the difference between the quantity in Eq. (13) and its CP conjugate we need to determine the term in the Lagrangian that is changed by CP conjugation. We will refer here only to the quark contribution. It is helpful to remove the CP violating part from the quark mixing term to the down quark mass matrix (and this will not alter the results). We denote by Md the diagonal mass of the down quarks. We perform the change of variables:

Ud = d′ (14) where U is the CKM matrix and d represents all down quark states. Then the mass term of the down quarks becomes:

′ † ′ d¯ UMdU d . (15)

We denote:

† t t t t A + iB = UMdU = (U1MdU1 + U2MdU2)+ i(U2MdU1 − U1MdU2), (16)

† where A and B are the real and imaginary parts of UMdU and U1 and U2 are the real and imaginary parts of the CKM matrix. We separate the relevant term in the Lagrangian according to:

L = L1 − d¯(A + iB)d, (17) for further convenience. Then one can write Eq. (13) as:

∞ 1 exp[iΓ] = dAa dhdΨ¯dΨ [ d4xdBd¯ ]n exp[i d4x[L(Φ) − dAd¯ ]] = µ n! Z n=0 " Z # Z ∞ X 1 ∂ dAa dhdΨ¯ dΨ [iB ]n exp[i d4xL(Φ) − dAd¯ ]= µ n! ij ∂A Z n=0 " ij # Z X∞ (i)n ∂ dAa dhdΨ¯ dΨ [V B V t ]n exp[i d4xL(Φ) − dAd¯ ]= µ n! ki ij kj ∂m′ Z n=0 " dk # Z ∞ X (i)n ∂ [V B V t ]n exp[iΓ(A)], (18) n! ki ij kj ∂m′ n=0 " dk # X where Γ(A) denotes the effective action with the mass matrix A for the down quarks. Here V is the unitary matrix t ′ that diagonalizes A according to V AV = Ad where Ad is the matrix of eigenvalues (mdk)of the matrix A. We also used:

∂ ∂VkiAij Vjk ∂ ∂ = ′ = VkiVjk ′ . (19) ∂Aij ∂Aij ∂mdk ∂mdk It is straightforward to deduce that:

∞ (−i)n d exp[iΓ]¯ = [V B V t ]n exp[iΓ(A)]. (20) n! ki ij kj dm′ n=0 " dk # X since:

L(CP )= L1 − d¯(A − iB)d. (21) 4

Then Eqs. (18) and (20) lead to:

∞ (i)n − (−i)n d exp[iΓ] − exp[iΓ]¯ = [V B V t ]n exp[iΓ(A)] = n! ki ij kj dm′ n=0 " k # X (i)n ∂ 2 [V B V t ]n exp[iΓ(A)]. (22) n! ki ij kj ∂m′ n=2k+1 " dk # X

IV. BARYON ASYMMETRY IN THE STANDARD MODEL

The one loop temperature dependent effective potential for the standard model has the expression [12], [13]:

T Veff (Φ,T )= V0(Φ) + V1(Φ) + V1 (Φ,T ) (23) where,

2 m 2 λ 4 V0(Φ) = − Φ + Φ 2 4 2 2s 1 4 mi (φ) V1(Φ) = n (−1) i m (Φ) ln[ ] − C , (24) i 64π2 i µ2 i i " # X Here i represents the particle, si its spin, ni the number of degrees of freedom associated to the specific particle, µ is the renormalization scale and Ci are constants depending on the renormalization scheme. For example ng = 3 for the 3 gauge bosons, nq = 12 for the quarks and nl = 4 for the leptons. For our calculations relevant is only Cq = 2 . The one loop thermal corrections are given by:

T 4 m2 T 4 m V T = n J ( i ) − n J ( j ), (25) 1 i 2π2 b T 2 j 2π2 f T 2 i=boson j=fermion X X where Jb and Jb are loop functions with the expression:

∞ 2 2 2 2 Jb[y ]= dxx ln 1 + exp[− x + y ] 0 " # Z p ∞ 2 2 2 2 Jf [y ]= dxx ln 1 − exp[− x + y ] . (26) 0 " # Z p The following approximation for y ≪ 1 will be useful in what follows:

4 2 2 2 7π π 2 1 4 y 3 Jf (y ) ≈ − y − y ln( )+ O(y ), (27) 360 24 32 af where ln(af ) ≈ 2.6351. In practice we will use only the first two terms in the expression in Eq. (27) as the other corrections turn out too small. Then the effective action is: 1 Γ(Φ) = −Vt V = − V . (28) f eff T 4 eff 1 Here we considered a thermal bath characterized by the volume Vtf = T 4 where T is the temperature and tf is the corresponding time [11]. We apply Eq. (18) to the potential and effective action in Eqs. (23) to obtain:

∞ (i)n ∂ [V B V t ]n exp[iΓ (A)] ≈ n! ki ij kj ∂m′ eff n=0 " dk X ∞ n in 3 m′2 1 (−i)n V B V [− m′3 [ln[ dk ] − 1]+ m′ T 2] exp[iΓ (A)]. (29) n! ki ij kj 4π2 dk µ2 2 dk eff n=0 " # X 5

Here the effective potential is considered for the masses of the down quarks as extracted from the matrix A and this does not correspond to the full down quark masses as in the standard model.. Then one can derive immediately:

exp[−Γ] − exp[−Γ]¯ = n (i)n 3 m′2 1 2 V B V (− m′3 [ln[ dk − 1]+ m′ T 2] exp[−Γ (A)] = n! ki ij kj 4π2 dk µ2 2 dk eff n=2k+1 " # X ′2 3 ′3 mdk 1 ′ 2 2i sin VkiBij Vkj [− 2 mdk[ln[ 2 ] − 1]+ mdkT ] exp[−Γeff (A)], (30) " 4π µ 2 # where this time the calculations are done in the euclidean space. Then from Eqs. (11) and (30) the baryon asymmetry in the Universe reads:

△BB = 12δ(0)Vtf [exp[−Γ(T )] − exp[−Γ(¯ T )] = ′2 3 mdk 1 ′2 2 i24δ(0)Vtf sin VkiBij Vkj (− 2 [ln[ 2 − 1]+ mdkT ] exp[−Γeff (A)] = " 4π µ 2 # ′2 1 3 mdk 1 ′2 2 i24δ(0) 4 sin VkiBij Vkj (− 2 [ln[ 2 − 1]+ mdkT ] exp[−Γeff (A)]. (31) T " 4π µ 2 #

Here all the masses are computed at the minimum of the thermal corrected potential.

V. NUMERICAL CALCULATIONS

For the electroweak baryogenesis the temperature should be at the electroweak scale and it is standard to take it T ≈ 100 GeV [8]. Moreover the perturbation theory works at best when µ = max[mi(Φ)] so we will take µ = mt = 173.3 GeV. All the masses for the standard model fermions, Higgs or gauge bosons are taken at the electroweak scale and have the usual values indicated in [14]. An exception are the down quarks for which the mass matrix eigenvalues are calculated according to Eqs. (16). Thus in the effective action the mass A is diagonalized by the matrix V such that:

t V AV = Ad t t A = (U1MdU1 + U2MdU2). (32)

The matrices U1, U2 are the real and imaginary parts of the CKM matrix and Md is the diagonal mass matrix of the down quarks in the standard model computed at the minimum of the thermal potential. Then the eigenvalues of the down quarks that appear in the effective action are slightly different than those in the standard model:

′ −3 −3 md1 = z5.42 × 10 as compared to md1 =4.7 × 10 GeV ′ md2 = z0.096 ascomparedto md2 =0.096 GeV ′ md3 = z4.179 ascomparedto md3 =4.18 GeV, (33)

Φmin where mdk are the down quark masses in the standard model. Here z = v where Φmin is the minimum of the full temperature dependent effective potential. Moreover for completeness we give here also the matrix B:

t t B = (U2MdU1 − U1MdU2). (34) where again Md is the diagonal mass of the down quarks computed at Φmin. 1 The factor Vtf is replaced by T 4 where T is the temperature for the thermal bath. Next we need to regularize the δ(0) function for which we adopt the simplest procedure [15].

d4k −[i∂ ]2 δ(0) = lim δ(x − y) = lim exp[ µ ] exp[i(x − y)] = x→y x→y (2π)4 M 2 Z d4k −k2 µ4 lim exp[ ] exp[i(x − y)] = i , (35) x→y (2π)4 µ2 16π2 Z 6 where µ is the renormalization scale and practically the highest scale in the theory besides the Higgs . Next step is to determine the value of the effective action. The effective potential for the standard model at the finite temperature T has been calculated retaining only the relevant terms in [12]. It has the expression:

λ(T ) V (Φ ,T )= D(T 2 − T 2) − ET Φ3 + Φ4, (36) c 0 c 4 c where the coefficients are: 2m2 + m2 +2m2 D = W Z t 8v2 2m2 + m2 E = W Z 4πv3 m2 − 8Bv2 T 2 = h 0 4D 3 B = (2m4 + m4 − 4m4) 64π2v4 W Z t 2 2 2 3 4 mW 4 mZ 4 mt λ(T )= λ − 2 2 2mW ln[ 2 ]+ mZ ln[ 2 ] − 4mt ln[ 2 . (37) 16π v " ABT ABT AF T #

3 3 Here ln AB = ln ab − 2 where ln ab =5.4076 and ln AF = ln af − 2 where ln af =2.6351. ∂V The minimum of the potential is compute by considering = 0for T = 100 GeV. and it is given by Φmin = 208.77 ∂ΦC GeV. Then the effective action at the minimum can be computed as:

V (Φ ) Γ = − eff min . (38) eff T 4

As we mentioned before all the masses for the down quarks in Eq. (31) must be computed at Φmin such that they Φmin gain a scale factor v =0.848 = z. Considering all the evaluations in the present section Eq. (31) can be calculated and leads to: ¯ △BB (△B −△B) −27 η1 = = ≈−5.145 × 10 , (39) g∗ 2g∗ where g∗ = 106.75.

VI. CONCLUSIONS

In this work we computed the baryon asymmetry in the standard model assuming that the only source of CP viola- tion is the CKM matrix and without considering the exact processes that take place at the electroweak baryogenesis. Our method is clear and straightforward and can be easily generalized to any kind of standard model extensions. The result of the present work is in good agreement with other estimates in the literature calculated in a more complicated manner. Thus the estimate,

nb − n¯b −27 η1 = ≈−5.145 × 10 , (40) s is very close and within the limits calculated in [8]: n − n¯ |η| = | b b |≤ 6 × 10−27. (41) s It would be of great interest to perform a similar calculation for the most common supersymmetric extensions of the standard model. This will be done in further work.

[1] A. D. Sakharov, JETP Lett. 5, 24 (1967). 7

[2] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 155, 36 (1985). [3] M. Dine, P. Huet and R. Singleton Jr.,NUcl. Phys. B 375, 625 (1992). [4] M. Dine, P. Huet, R. G. Leigh, A. Linde and D. Linde , Phys. Lett. B 283, 319 (1992); Phys. Rev. D 46, 550 (1992). [5] A. G. Cohen, D. B. Kaplan and A. E. Nelson, Ann. REv. NUcl. Part. Sci. 43, 27 (1993). [6] G. R. Farrar and M. E. Shaposhnikov, Phys. Rev. Lett. 70, 2833 (1993). [7] M. B. Gavela, P. Hernandez, J. Orloff and O. Pene, preprint 1993, CERN 93/7081, LPTHE Orsay-93/48, HUTP-93/48, HD-THEP-93-45. [8] P. Huet and E. Sather, Phys. Rev. D 51, 379-394 (1995). [9] M. Trodden, Rev. Mod. Phys. 71, 1463-1500 (1999). [10] A. Riotto, CERN-TH-98-204 Conference: C98-06-29.8, p. 326-436 Proceedings, arXIv:hep-ph/9807454 (1998). [11] J. M. Cline, Lectures given at Conference: C06-07-31.1, arXiv:hep-ph/0609145 (2006). [12] M. Quiros, IEM-FT-187-99 Conference C-98-06-29.8, p. 187-259 Proceedings, arXiv:hep-ph/9901312 (1999). [13] D. E. Morrissey and M. J. Ramsey-Musolf, New J. Phys. 14, 125003 (2012) [14] C. Patrignani et al. (Particle Data Group), Chin. Phys. C 40, 100001 (2016). [15] M. Srednicki, ”Quantum Field Theory”, Cambridge University Press, New York 2010 (pg.472-477).