1. Introduction
1.2 CP violation
1.2.4 CP violating experimental data and CKM model
CP violation was first discovered in Kaon mixing [7] in 1964. The CP eigenstates for the K0 and K0 system are
K1 = 1
√2(K0− K0);
K2 = 1
√2(K0+ K0). (1.45)
K0 is the pesudoscalar particle with odd parity and in charge conjugatation trans-form CK0 = K0; CK0 = K0. It is obvious that K1 is CP even eigenstate and K2 is
Fig. 1.2: Box diagrams for neutral K-meson mixing
the eigenstate of CP odd. If the CP eigenstates are also the Hamiltonian eigenstates, it means that the CP is conserved under the system.
In general
H = M − i 2Γ =
M11 M12 M21∗ M22
− i 2
Γ11 Γ12 Γ∗21 Γ22
, (1.46)
where the M and Γ are 2 × 2 Hermitian matrices, so H is not Hermitian obvi-ously. This non-Hermitian Hamiltonian make the two state system decay during time evolution. The mass eigenstates are
KS = 1
p1 + |˜²|2(K1− ˜²K2),
KL = 1
p1 + |˜²|2(K2+ ˜²K1). (1.47)
Where |˜²| = (2.44 ± 0.04) × 10−3 [34] is the small value related to the mixing of two CP eigenstates. This formula indicates that the mass eigenstates(energy eigenstates) are not exactly identical to the CP eigenstates. This experimental data is explained by so-called box diagram shown in Fig.1.2 in the SM.
Direct CP violation in Kaon decay into ππ has also been discovered[35]. When treating K meson decay, we usually take the mass eigenstates KS and KL as the CP eigenstates instead of K1 and K2, because here we discuss only the CP violation from decay. It is convenient to define the quantities related to the decay amplitude
Fig. 1.3: Tree and penguin diagrams for K → ππ decays, where q can be u or d
for K meson to study direct CP violation,
η00 = < π0π0|H|KL>
< π0π0|H|KS >; η+− = < π+π−|H|KL>
< π+π−|H|KS >. (1.48) CP violation in K → ππ is measured by ²0, which is defined as
η00= ˜² − 2²0;
η+−= ˜² + ²0. (1.49)
The Particle Data Group [12] gives the fitting for the value
Re(²0/˜²) = 1.65 ± 0.26. (1.50)
²0 is explained in the SM by the tree and penguin diagrams in Fig.1.3.
B decays can provide many tests for CKM model by measuring α, β, and γ in the unitary triangle in Fig 1.1.
β is the relative angle between VtdVtb∗ and −VcdVcb∗ on the complex plane. There are several ways people often take to determine this angle. The most popular process is the b → c¯cs process. This is because the amplitude of the tree level and loop diagram has approximately the same phase. One of the process often been used is B → J/ψKs. The sin 2β is extracted from the relation [12]
Sf = −ηfsin 2β, (1.51)
where Sf is the quantity related to time-dependent CP asymmetry in B decays [36, 37], and ηf is the CP eigenvalue of f . The experimental result from average of the related decay by BaBar has the value [38]
sin 2β = 0.686 ± 0.039 ± 0.015. (1.52) From Fig 1.1 definition the α is the angle between VudVub∗ and −VtdVtb∗. It can be extract from B → ππ process, via the measurement of Sπ+π− and Cπ+π−. The measurement of BABAR gives that [40]
α = 96◦+10−6◦ ◦. (1.53)
The decay B → ρ0π0 are also applied to determined α, and the experimental result from Belle gives 68◦ < α < 95◦ at CL = 68%[41].
γ is the angle between VcdVcb∗ and −VudVub∗. The measurement of γ determination uses the B decay process B → DK. BELLE [42] measured B− → DK−, B− → D∗K− and B−→ DK∗− to obtain
(γ = 53+15−18± 3 ± 9)◦ (1.54)
The process B → K+π−are usually applied to test the CP asymmetry, too. The experimental average is ACP(B0 → K+π−) = −0.097 ± 0.012 by HFAG [43]. This asymmetry can be explained by SM [44, 45, 46].
The global fit for the unitary triangle is summarized in Fig.1.4 which is from CKMfitter [47]. The PDG review[12] provides the fitting values for Wolfenstein parameters, which are
λ = 0.2272 ± 0.0010; A = 0.818+0.007−0.017;
ρ = 0.221+0.064−0.028; η = 0.340+0.017−0.045. (1.55) Compare these values with the PDG parameters, we can derive s12 = 0.227 ± 0.001, s23= 0.0422±0.0004, s13= 0.00399±0.00007, and the phase sin δ13= 0.839±
0.006. From these values, one obtain[12]
J = (3.08+0.16−0.18) × 10−5. (1.56)
-1 -0.5 0 0.5
ρ
1 1.5 2 excluded area has CL > 0.95 excluded at CL > 0.95Summer 2007
CKMf i t t e r
Fig. 1.4: The experimental fit for the ρ and η from CKMfitter Group [47]
From the above discussion for CKM matrix, we can see that the CKM matrix works very well in describing the meson decay, leptonic and semileptonic decay.
The phase in CKM matrix generates CP violation, and which is consistent with the experimental result for the CP violation phenomena like K meson mixing, K and B meson decay. However there are still some problems. One of them is the baryongenesis. That is, the amount of particles is more than that of antiparticles in our world, and one necessary condition for this phenomenon is the existence of CP violation. The quantity to estimate the asymmetry of universe is nB/nγ, with nB, nγ
denoting the baryon number density, and photon number density respectively. In high temperature the CKM model can produce about nB/nγ ≈ O(10−20)[48], which is too small compared with observation nB/nγ ≈ 10−8[49]. There should be another source of CP violation beyond the CKM matrix in our world. Also, the CKM model does not provide the answer where CP violation is originated, but just put in by hand. It is desirable for some understanding of the origin of CP violation. In this
thesis we try to study how to connect the CKM matrix phase with the spontaneous CP violating phase for the explanation of the source of CP violation from CKM mechanism.
In chapter 2 we are going to discuss what the spontaneous CP violation is and treat some of the multi-Higgs models. In chapter 3 we will build a new model with the connection between spontaneous CP violating phase and CKM matrix phase. In chapter 4 we will use this model to discuss some phenomenology. In the last chapter we will summarize what we do in this thesis.
VIOLATION
Although the CKM matrix and its complex relation explain the CP violation of observation very well, it is still possible that CP is violated from other place. The multi-Higgs model is a popular topic in this area. Such models may also answer that CP violation comes from the so-called spontaneous CP violation(SCPV), a mechanics first proposed by T.D. Lee [50, 51].
When there are more than one Higgs, their vacuum expectation values might have the relative phases difference. If these phases are non-vanishing after symmetry breaking and irreducible in Higgs self-interaction or Yukawa terms with fermions, then they also produce CP violation. Because this kind of CP violation comes from the spontaneous symmetry breaking of Higgs, it is called spontaneous CP violation.
2.1 Two Higgs Doublet Model
In 1973, T. D. Lee proposed a model with two Higgs doublets [50, 51]. The most important property of this model is that if there is a phase difference between VEVs of two Higgs doublets, this phase can give the contribution to the CP violation. The two Higgs doublets are written in the following form [50]
φ1 = eiθ1H1 = eiθ1
√1
2(ρ1+ R1+ iA1) h−1
;
φ2 = eiθ2H2 = eiθ2
√1
2(ρ2+ R2+ iA2) h−2
. (2.1)
H1 and H2 are Higgs with real vacuum expectation values, and R1, R2, A1, and A2 are real parts and imaginary parts of them. The phase difference between two
Higgs doublets δ = θ2 − θ1 is the spontaneous CP violating phase. If this phase is non-zero and can’t be eliminated by fermion rotation, then this could produce the spontaneous CP violation. The Higgs potential can be built in the form [50]
V = − λ1φ†1φ1− λ2φ†2φ2+ A(φ†1φ1)2+ B(φ†2φ2)2 + C(φ†1φ1)(φ†2φ2) + C(φ†1φ2)(φ†2φ1)
+ 1
2[(φ†1φ2)(Dφ†1φ2+ Eφ†1φ1 + F φ†2φ2) + h.c.], (2.2) where λ1, λ2, A − F , and C are all real numbers. The minimal condition by differ-entiating with respect to δ can give
cos δ = −(4Dρ1ρ2)−1[Eρ21+ F ρ22]. (2.3)
If the right handed side is not required to be 1 or -1, CP is violated spontaneously.
There are eight real scalar fields in two Higgs doublets, and three of them are eaten by W± and Z0. So there are five physical states in two Higgs doublet model after spontaneous symmetry breaking. The two Higgs model is different from SM because it has charged Higgs bosons. The interaction of charged Higgs and fermions is similar to charged weak interaction. So this model also have more contribution to flavor change process than SM.
Usually the two Higgs doublet models are classified into three types by different Yukawa interactions.
Type I Type I is that one Higgs couples with each fermions, like the Higgs in standard model, and another Higgs does not couple with fermions as below
QLΓuφ1UR+ QLΓdφe1DR+ LLΓeφe1ER+ h.c., (2.4)
where Γu, Γd, and Γe are real coupling matrices. We construct thise inter-actions of type I by introducing the discrete symmetry φ2 → −φ2 and other fields are unchanged, so that φ2 only exists in the Higgs potential with even powers.
Type II Type II is that one Higgs doublet couples with up-type quarks and another one couples with down-type quarks. For example
QLΓuφ1UR+ QLΓdφe2DR. + LLΓeφe2ER+ h.c.. (2.5)
We can construct type II model by introducing this discrete symmetry:
φ1 → φ1; φ2 → −φ2; UR→ UR; DR→ −DR; ER→ −ER;
QL → QL; LL→ LL. (2.6)
Type III The last model type III is the most general Yukawa interactions in which there are two Higgs coupling with each fermion,
QL(Γu1φ1+ Γu2φ2)UR + QL(Γd1φe1+ Γd2φe2)DR
+ LL(Γe1φe1+ Γe2φe2)ER+ h.c. (2.7)
so there could be the FCNC process because one can not diagonalize mass matrix v1Γ1+ v2Γ2 and coupling matrices Γ1,2 simultaneously.
Type I and II can not have spontaneous CP violation, because (φ†1φ2)(φ†1φ1), and (φ†1φ2)(φ†2φ2) are not allowed, and this results in sin δ = 0. The spontaneous CP violating phase δ vanishes. So type I and type II with discrete symmetry have no spontaneous CP violation.
If we hope that the spontaneous CP violation exists, then only type III is allowed.
However, type III has the tree level FCNC contribution which is severely constrained by experimental data. Also, there are also too many unknown parameters in the model.
2.2 Weinberg Model
In 1976, Weinberg [52] proposed that by using some discrete symmetry, the model with three or more Higgs doublets gives the spontaneous CP violation without tree level flavor change neutral current. The three Higgs doublet model is called Weinberg
model. Three Higgs φ1, φ2, and φ3 are written as follows
φk = eiθkHk = eiθk
√1
2(vk+ Rk+ iAk) Hk−
, k = 1 ∼ 3, (2.8)
where √12vkeiθk is the vacuum expectation value of neutral part in φk, and we let Hk0 = vk+ Rk+ iAk. Rk and Ak are the corresponding real part and imaginary part in Hk. Branco extended this idea to arbitrary number of generations with two sets of discrete symmetry [53]
D1 : φ1 → φ1; φ2 → −φ2; φ3 → φ3; QL → QL; dR→ dR; uR → −uR D2 : φ1 → φ1; φ2 → φ2; φ3 → −φ3; QL → QL; dR→ dR; uR → uR. (2.9)
The D1 implies the constraint that φ1 couples to the up-type quarks singlet uR and φ2 couples to the down-type ones dR, which has the same propose as the discrete symmetry for type II of two Higgs doublet model. The D2 can suppress φ3 not to couple with quarks, but it can couple to leptons. Applying those discrete symmetry can inhibit the tree level FCNC process of Higgs exchange, and the CP violation can arise from the Higgs interaction themselves. The Yukawa terms of Weinberg model are written as
QLΓuφ1UR+ QLΓdφe2DR+ LLΓeφe3ER+ h.c.. (2.10)
These interactions are similar to type II of two Higgs doublet model, and the spon-taneous CP violation will be produced in the Higgs potential.
After spontaneous symmetry breaking the Lagrangian is expanded as follows
L = −1
v1ULMˆuURH10 − 1
v2DLMˆdDRH20− 1
v2LLMˆeERH30
−
√2
v1 DLVCKM† MˆuURH1−+
√2
v2 ULVCKMMˆdDRH2+ +
√2 v3
LLVCKMMˆeERH3++ h.c., (2.11)
where Mu = −√12Γuv1, Md = −√12Γdv2, and Me = −√12Γev3. VCKM is assumed to be real matrix here. That is, CP violation does not come from the CKM matrix.
The CP should be arisen from the Higgs self-interaction. The parametrization for three Higgs doublet potential in discussion [54] is
V = − µ21φ†1φ1− µ22φ†2φ2− µ23φ†3φ3+ h1(φ†1φ1)2+ h2(φ†2φ2)2+ h3(φ†3φ3)2 + f12(φ†1φ1)(φ†2φ2) + f23(φ†2φ2)(φ†3φ3) + f31(φ†3φ3)(φ†1φ1)
+ g12(φ†1φ2)(φ†2φ1) + g23(φ†2φ3)(Dφ†3φ2) + g31(φ†3φ1)(φ†1φ3)
+ (k12(φ†1φ2)2+ k23(φ†2φ3)2+ k31(φ†3φ1)2+ h.c.). (2.12) It can be assumed that the coefficients in above formula are all real. In the potential only two phases δ12 = θ2 − θ1 and δ23 = θ3− θ2 exist. Differentiating with respect to the two phases, we get a condition below
k12v22
v32 sin 2δ12= k23v22
v21 sin 2δ23 = −k13sin 2(δ12+ δ23). (2.13) Eq.(2.13) reflects the fact that the phases δ12 and δ23 can be nonzero, and CP is violated here.
In this model, there are four charged and five neutral physical Higgs bosons. The mass matrix of Higgs can give CP properties of this model [54]. The resulting mass matrix of charged Higgs has the off diagonal complex elements. It means that the CP violation will arise from the exchanging of charge Higgs. Also, the neutral mass matrix has the mixing terms between scalar and pseudoscalar Higgs which lead to the CP violation under neutral Higgs exchange.
Although Weinberg model can provide the spontaneous CP violation without the tree level FCNC which is inhibited in Lee’s two Higgs doublet model, it still has some contradiction which had been provided by many authors[55, 56, 57]. The Weinberg model is decisively ruled out by data on sin 2β measurement in B → KsJ/ψ. In Weinberg model, the upper bound for magnitude of sin 2β is | sin 2β| < 0.05[56, 57].
The present experimental data is shown in Eq.(1.52) that sin 2β = 0.686 ± 0.039 ± 0.015, which means that the Weinberg model has been ruled out. The neutron EDM calculation also rules out the Weinberg model, from which the estimation for neutron EDM has order 10−23e cm [58], whereas the experimental upper bound is
|dn| < 0.29 × 10−25e cm [59].
2.3 The strong CP problem and Peccei-Quinn symmetry
For the Lagrangian in QCD, the term (θgs2/32π2)GµνGeµν is allowed, where eGµν =
1
2²µνρσGρσ. This term also violates P and CP. This is a possible CP violation in strong interaction.
Because P and CP violation will cause the electric dipole moment(EDM) of particles, the measurement of EDM of particles is important to test the CP violation in standard model. The neutron EDM test is especially important, and in present the experimental upper bound is given [59]
|dn| < 0.29 × 10−25e cm. (2.14) SM theoretical calculation of CKM matrix CP violation gives the small value con-tribution about order less than 10−31e cm [61, 62, 63] without considering strong CP violation. With non-zero θ, dn can be much larger. Experimental bound for neutron EDM constrains the θ very strongly for |θ| . 10−10 [12]. This is considered to be unnatural since other couplings with strong interaction are much larger. This is the problem.
For the multi-Higgs model with spontaneous CP violating phase, the strong phase θ would be large [60] at loop level. We need a mechanism to make this phase small.
In 1977, Peccei and Quinn proposed a mechanism to solve this problem [64, 65].
They introduced another global symmetry U(1)PQ in the standard model. This symmetry is generated by the chiral transformation defined as follows.
u → eiαuγ5u; d → eiαdγ5d;
φ1 → ei(αu+αd)φ1; φi → e−i(αu+αd)φi; i 6= 1, (2.15) where αu and αd are the chiral rotational phases for up-type quarks and down-type quarks respectively. The φ1 and other φi are the multi-Higgs doublets. After the chiral rotation, the strong phase becomes
θ → θ − 4αu− 4αd. (2.16)
Since αu,d are arbitrary phases, one can choose these phase as θ = 4(αu + αd), therefore there is no strong CP phase and also without large contribution to neutron EDM.
Models with PQ symmetry have an axion resulting from spontaneous breaking down of PQ symmetry. No axion has been detected in experiments. One has to make the axion invisible, by extending the original PQ model[66, 67].
The multi-Higgs model can have the spontaneous CP violation(SCPV). This is a nice feature which provides a understanding of the origin of CP violation. But the two Higgs doublet model has tree level FCNC, leading to too many unknown parameters. To improve the situation, Weinberg proposed a three Higgs doublet model which has no tree level FCNC. However, the prediction of Weinberg model for sin 2β is not consistent with experimental data as mentioned before. Here we take the idea [1] that the CP violation is arisen from spontaneous symmetry breaking, but further make the spontaneous CP violating phase be identical to the CP violation in CKM matrix. In the following we build specific models to realize this.
3.1 Making SCPV phase identical to CKM phase
Model(a)
In our new model, we try to build a model with the spontaneous CP violating phase from Higgs identical to the the phase in CKM matrix. We couple two independent Higgs doublets to the up-type quarks and one Higgs doublet to the down type quarks as below
L = QL(Γu1φ1+ Γu2φ2)UR+ QLΓdφeddR+ h.c., (3.1) where Γu1 and Γu2 are real 3 × 3 coupling matrices, and the tilde sign on Higgs means eφk= −iσ2φ∗k. The φ1, φ2 and φd are Higgs doublets, which are written as
φk = eiθk
√1
2(vk+ Rk+ iAk) h−k
, (3.2)
where k can be 1, 2, and d. It is convenient to redefine these Higgs doublets so that they have real vacuum expectation values. That is, φk = eiθkHk. Here we call
Eq.(3.1) as model(a).
After spontaneous symmetry breaking, the mass terms of model(a) Lagrangian appear as
Lm = −UL(Mu1eiθ1 + Mu2eiθ2)UR− DL(Mde−iθd)DR+ h.c., (3.3)
where Mu1,u2 = −Γu1,u2v1,2/√
2 and Md = −Γdvd/√
2. The phase θ1 and −θd can be absorbed into the UR and DR respectively.
From previous section, Eq.(1.20) shows the relations between flavor eigenstates and mass eigenstates of quarks
UL= VLuULm; UR= VRuURm; DL = VLdDLm; DR= VRdDRm.
We make the Md to be diagonal without loss of generality. That is, DL and DR are already the mass eigenstates DLm and DmR, and the mixing matrix VLd and VRd are unit matrices. The mass terms become
Lm = −UL(Mu1+ Mu2eiδ)UR− DLMˆdDR+ h.c., (3.4)
where the relative phase δ = θ2− θ1 is the spontaneous CP violating phase in the Yukawa terms, and ˆMd is the diagonal mass matrix. If δ is non-zero, it could cause the spontaneous CP violation in the model. The total mass matrix of up-type quarks can be diagonalized by left and right handed matrices VRu and VLu. That is
Mˆu = VLu†MuVRu, (3.5)
where ˆMuis the diagonal up-quark mass matrix and Mu = Mu1+eiδMu2. To simplify the discussion we assume that VRu is a unit matrix. This simplification can help us to make the identity relation between δ and the phase in CKM matrix. Because VLd is a unit matrix, from VCKM = VLu†VLd it is obvious that VLu† is equal to VCKM. Eq(3.5) becomes ˆMu = VCKM(Mu1+ eiδMu2) , and we obtain the relation
VCKM† = (Mu1+ eiδMu2) ˆMu
−1. (3.6)
At this step we need the explicit CKM parametrization with a phase. First we use the Particle Data Group parametrization which is shown in Eq.(1.30)
VCKM =
This parametrization makes more than one phases in the VCKM elements. That is, one phase δ13 in V21, V22, V31, V32 and another phase −δ13 in V13. We solve the problem by pulling the phase −δ13 out and then decomposing the VCKM into two matrices as below,
Absorbing the left handed diagonal matrix by redefining the quark sector UL. The remaining matrix has the uniform phase in each element. By comparing two side of Eq.(3.6), we introduce the identical relation
δ = −δ13. (3.8)
This relation implies that the CKM phase comes from the spontaneous CP violating phase. Also, this relation is related to the phase δ13 which has been measured in experiments, so if the spontaneous CP violating phase δ can be nonzero after solving the minimal condition in Higgs potential, it is independent of the masses of Higgs, and the CP phenomena always exists.
By substituting Eq.(3.8) into Eq.(3.6) we determine the coupling matrices Mu1
and Mu2
Note that these matrices Mu1 and Mu2 depend on quark masses and the angles in CKM parametrization. This is not true for other multi-Higgs models, with which the spontaneous CP violating phase and CKM phase are concerned.
We apply the same idea to another CKM parametrization, the original Kobayashi Maskawa parametrization in Eq.(1.29)
Then the step of previous discussion for PDG parametrization makes the same relation as Eq.(3.8),
δ = −δKM. (3.10)
This relation gives the expression for the mass matrix with respect to KM parametriza-tion,
In general we can express the two coupling matrices Mu1 and Mu2 in terms of the CKM matrix, quark mass matrices, and the spontaneous CP violating phase as
follows,
Mu1 = VCKM† Mˆu− eiδ
sin δIm(VCKM† ) ˆMu; Mu2 = 1
sin δIm(VCKM† ) ˆMu. (3.12) This relation is useful when we treat the Yukawa couplings of Higgs, and it is inde-pendent of the parametrization of VCKM. Choosing a specific CKM representation implies a choice of a model.
Model(b)
Now we treat another kind of Yukawa interactions which is called model(b),
L = ¯QLΓuφuUR+ ¯QL(Γ1φe1+ Γ2φe2)dR. + h.c.. (3.13) In this Lagrangian two Higgs doublets φ1 and φ2 couple with down-type quarks and one Higgs-doublet φu is coupled with up-type quarks. Γd(1,2) and Γu are 3 × 3 real matrices. Here the φ1 and φ2 are defined as those in model(a), and φu is with the same definition as that in Eq.(3.2) and in which k is replaced by u. After symmetry breaking the mass terms is written below
Lm = −UL(Mueiθu)UR− DL(Md1e−iθ1 + Md2e−iθ2)DR+ h.c., (3.14) where Md(1,2) = −Γd(1,2)v1,2/√
2 and Mu = −Γuvu/√
2. Following the same treat-ment for model(a), we absorb the phases −θ1 and θu into DR and UR respectively, and we also treat the Mu as diagonal mass matrix. So VLu = VRu = 1 and the mass terms become
Lm = −ULMˆuUR− DL(Md1+ Md2e−iδ)DR+ h.c.. (3.15) The ˆMu indicates that it is diagonal, and δ is the spontaneous CP violating phase with δ = θ2 − θ1. The diagonal down-type mass matrix ˆMd has the relation to Md= Md1+ e−iδMd2,
Mˆd= VL†dMdVRd, (3.16) We make VRd= 1 and the VLd† is equal to VCKM† . That makes us to express the mass matrices Md1 and Md2 in terms of VCKM as
VCKM= (Md1+ e−iδMd2) ˆMd−1. (3.17)
Using the PDG parametrization with the same argument from Eq.(3.7), and
Using the PDG parametrization with the same argument from Eq.(3.7), and