arxiv:hep-ph/ v1 22 Jun 2005

Transcription

1 IPPP/05/9 DCPT/05/58 March 18, 008 On the Neutral Scalar Sector of the General R-parity Violating MSSM arxiv:hep-ph/050609v1 Jun 005 A. Dedes (1, S. Rimmer (1, J. Rosiek (, M. Schmidt-Sommerfeld (1 (1 Institute for Particle Physics Phenomenology (IPPP, Durham DH1 3LE, UK ( Institute of Theoretical Physics, Warsaw University, Hoza 69, Warsaw, Poland ABSTRACT Starting out from the most general, gauge invariant and renormalizable scalar potential of the R- parity violating MSSM and performing a calculable rotation to the scalar fields we arrive at the so called Grossman-Haber basis where the sneutrino VEVs are zero. The advantage of our rotation is that, in addition, we obtain diagonal soft supersymmetry breaking sneutrino masses and all potential parameters and VEVs real, proving that the MSSM scalar potential does not exhibit spontaneous or explicit CP-violation at tree level. The model has five CP-even and four CPodd physical neutral scalars, with at least one CP-even scalar lighter than M Z. We parametrise the neutral scalar sector in a way that resembles the parametrisation of the R-parity conserving MSSM, analyze its mass spectrum, the coupling to the gauge sector and the stability of the potential. 1

2 1 Introduction That none of the terms in the Standard Model (SM violate lepton number (L is not due to an imposed symmetry, but merely reflects the fact that all such combinations of SM fields are ruled out by consideration of gauge invariance and renormalisability [1]. For supersymmetric extensions of the SM this is no longer true. In the Minimal Supersymmetric Standard Model (MSSM [], lepton number violating terms (and baryon number (B violating terms appear naturally, giving rise to tree-level processes, proton decay for example, which are already strongly constrained by experiment. Either, bounds can be set on Lagrangian parameters, or a further discrete symmetry can be imposed on the Lagrangian, such that these processes are absent. The discrete symmetry most commonly imposed is known as R-parity (R P [3, 4]. Under R-parity the particles of the Standard Model including the scalar Higgs fields are even, while all their superpartners are odd. Imposing this symmetry has a number of effects. Firstly, any interaction terms which violate lepton number or baryon number will not appear. Secondly, the decay of the lightest supersymmetric particle (LSP into SM particles would violate R P ; the LSP is therefore stable. The sneutrino vacuum expectation values (VEVs are zero; without extending the MSSM field content, spontaneous generation of R P violating terms is phenomenologically discounted [5]. If R P conservation is not imposed, fields with different R P mix [6 1]. In particular, the neutrinos will mix with the neutralinos and the sneutrinos will mix with the neutral scalar Higgs fields; all five complex neutral scalar fields can acquire vacuum expectation values. Minimising this ten parameter potential in general is not straightforward, it is more convenient to simplify the system by choosing an appropriate basis in the neutral scalar sector. As one of the Higgs doublets carries the same quantum numbers as the lepton doublets (apart from the non-conserved lepton number, it is convenient to introduce the notation L α = (H 1, L i where H 1 and L i are the chiral superfields containing one Higgs doublet and the leptons, respectively (α = 0,..., 3 and i = 1,...,3. Furthermore, starting from the interaction basis, we are free to rotate the fields and choose the direction corresponding to that of the Higgs field. Assuming that the system defining the five complex vacuum expectation values of the fields was solved, four complex VEVs v α would define a direction in the four dimensional (H 1, L i space. One can then choose the basis vector which defines the Higgs fields to point in the direction defined by the vacuum expectation values. We refer to this basis, in which the sneutrino (as we call the fields perpendicular to the Higgs field VEVs are zero, as the Grossman-Haber (GH basis [13]. This basis has the

3 virtue of simplifying the mass matrices and vertices of the theory and thus is better suited for calculations. Basis independent parameterisations can be chosen which explicitly show the amount of physical lepton number violation [14 16]. Values for physical observables such as sneutrino masses and mass splitting between CP-even and CP-odd sneutrinos have been derived in the literature in terms of these combinations but usually under some approximations (for example the number of generations or CP-conservation. We find this procedure in general complicated for practical applications and we shall not adopt it here. Instead, we present in the next section a calculable procedure for framing the most general MSSM scalar potential in the GH basis. An advantage of our procedure is to obtain a diagonal slepton mass matrix, two real non-zero vacuum expectation values and real parameters of the neutral scalar potential in the rotated basis. The latter proves that the neutral scalar sector of the most general R-parity violating MSSM exhibits neither spontaneous nor explicit CP-violation in agreement with [17]. In Section 3, the tree-level masses and mixing of the neutral scalar sector is investigated. Using the Courant-Fischer theorem for the interlaced eigenvalues, we prove that there is always at least one neutral scalar which is lighter than the Z-gauge boson. We present approximate formulae which relate the Higgs masses, mixing and Higgs-gauge boson vertices of the R-parity conserving (RPC case with the R-parity violating (RPV one. In Section 4, the positiveness of the scalar mass matrices and stability of the vacuum is discussed. Basis choice in the neutral scalar sector In this section we develop a procedure connecting a general neutral scalar basis with the Grossman-Haber (GH basis with three vanishing sneutrino VEVs, the latter being more convenient for practical applications. The most general, renormalizable, gauge invariant superpotential that contains the minimal content of fields, is given by W = ǫ ab [ 1 λ αβk L a α Lb β Ēk + λ αjk La α Qb x j D k x µ α L a α Hb + (Y U ij Q a x i H b Ūj x ] + 1 ǫ xyz λ ijk Ūx i D y j D z k, (.1 where Q a i x, Dx i, Ūi x, La i, Ē i, H1 a, Ha are the chiral superfield particle content, i = 1,, 3 is a generation index, x = 1,, 3 and a = 1, are SU(3 and SU( gauge indices, respectively. The simple form of (.1 results when combining the chiral doublet superfields with common 3

4 hypercharge Y = 1 into La,...,3 = (H a 1, L a i=1,,3. µ α is the generalized dimensionful µ- parameter, and λ αβk, λ αjk, λ ijk, (Y U ij are Yukawa matrices with ǫ ab and ǫ xyz being the totally anti-symmetric tensors, with ǫ 1 = ǫ 13 = +1. Then the five neutral scalar fields, ν Lα, h 0 from the SU( doublets, L α = ( ν Lα, ẽ Lα T and H = (h +, h 0 T, form the most general neutral scalar potential of the MSSM, V neutral = ( m L αβ ν Lα ν Lβ + µ αµ β ν Lα ν Lβ + µ αµ α h 0 h 0 + m H h 0 h 0 b α ν Lα h 0 b α ν Lαh (g + g [h 0 h 0 ν Lα ν Lα ], (. where general complex parameters b α, an hermitian matrix ( m L and αβ m H all arise from the supersymmetry breaking sector of the theory. The last term in (. originates from the D-term contributions of the superfields L α, H. Defining ( M L αβ ( m L one can rewrite the potential in (. in a compact form as V neutral = ( M L αβ + µ αµ β, and m m H + µ αµ α, (.3 αβ ν Lα ν Lβ + m h0 h0 (b α ν Lα h 0 + H.c (g + g [h0 h0 ν Lα ν Lα ]. (.4 In order to go to the GH basis, we redefine the Higgs-sneutrino fields where U is a 4 4 unitary matrix ν Lα = U αβ ν Lβ, (.5 U = V diag(e iφα Z, (.6 being composed of three other matrices which we define below, V unitary and Z real orthogonal. The potential in the primed basis becomes, [ ( ] V neutral = Z T ˆM L Z ν αβ Lα ν Lβ + m h0 h0 ] [(b Z α ν Lα h0 + H.c + 1 ( 8 (g + g h 0 h0 Lα ν Lα ν, (.7 where ( ˆM L = diag(e iφα V ( M L V diag(e iφ α, b T = b T V diag(e iφα. (.8 4

5 ( The unitary matrix, V, is chosen such that ˆM L is real and diagonal - the hat (ˆ is used to denote a diagonal matrix. The phases φ α are chosen such that b α is real and positive [they are equal to the phases of (b T V α]. The minimisation conditions for the scalar fields are now derived, to obtain conditions for the vacuum expectation values, V [ ( ] ν = Z T L ˆM Z ν Lβ (b Z α h 0 1 ( Lα αβ 4 (g + g h 0 h 0 ν Lγ ν Lγ ν Lα = 0, vac vac V = m h0 (b Z α ν Lα + 1 ( 4 (g + g h 0 h0 Lγ ν ν Lγ h 0 = 0, (.9 h 0 vac where vac indicates that the fields have to be replaced by their VEVs, vac ν Lα = v α, h 0 = v u. (.10 The U(1 Y symmetry of the unbroken Lagrangian was used to set the phase of v u to zero, however, at this stage all other vacuum expectation values will be treated as complex variables. By combining Eqs. (.9,.10 we obtain [ ( ] Z T L ˆM Z v β (b Z α v u 1 αβ 8 (g + g (v u v γ v γ v α = 0, (.11 m v u (b Z α v α (g + g (v u v γv γ v u = 0. (.1 In a general basis, it is difficult to solve the above system with respect to the VEVs without making some approximations, for example assuming small sneutrino VEVs [11]. In order to simplify calculations we would like to find a basis where the sneutrino VEVs vanish, v 1 = v = v 3 = 0. In other words, we are seeking an orthogonal matrix Z, such that the following equation, [ ( ] Z T L ˆM Z v 0 (b Z α v u 1 α0 M Z v u v 0 v u + v 0 v 0 δ 0α = 0, (.13 holds. If the above system is satisfied, then a solution with zero sneutrino VEVs exists. The other solutions, with non-vanishing sneutrino VEVs will be discussed later. In Eq. (.13, M Z = 1 4 ( ( g + g v u + v0, (.14 is the Z-gauge boson mass squared. It is obvious that when v i = 0, v 0 is real. It is now useful to define tanβ v u v 0. (.15 5

6 To determine Z, multiplying (.13 by Z γα, summing over α and solving for Z α0, yields, Z α0 = ( ˆM L b α tanβ αα 1 M tan β 1 Z tan β+1. (.16 For given set of model parameters, Z α0 depends only on tanβ which we can now fix by solving the orthonormality condition, Z α0 Z α0 = [( L ˆM b α tan β 1 M Z αα ] = 1. (.17 tan β 1 tan β+1 This equation can be easily be solved numerically for any given set of model parameters. It is worth ( noting that when b i = 0 and using notation more typical for this case, b 0 m 1, L ˆM m 1, Eq. (.17 reduces to one of the standard RPC MSSM equations 00 for the Higgs VEVs: m 1 v d = v u [ m (g + g (v u v d ]. (.18 For some parameter choices Eq. (.17 may admit multiple solutions for tan β. Each of the possible tan β specify a different basis, and each of these bases has one solution of the minimisation conditions with vanishing sneutrino VEVs. The subtlety highlighted earlier is the following: all possible solutions of the minimisation conditions can be found in each basis, so, in general, each basis contains a number of extrema equal to the number of possible solutions for tanβ. Hence, a solution with v i = 0 in one basis, is a solution with v i 0 in another basis. The important point to note is that by considering all possible values of tanβ, and selecting the value which corresponds to the deepest minima for the solution with vanishing sneutrino VEVs, all the solutions will have been accounted for, and the GH basis will have been determined correctly. The value of the potential at the vacuum, in terms of tanβ is given by M4 Z V (tan β = (g + g ( tan β 1. (.19 tan β + 1 The obvious conclusion from the equation above is that the deepest minimum of the potential is given by the solution for tanβ or cot β which is greatest. Knowing tan β, one should fix m using Eqs. (.1, (again in the analogy with RPC MSSM where m is usually given in terms of M A, tan β. Namely m = Z α0b α cot β 1 M Z 6 tan β 1 tan β + 1. (.0

7 In this way m is chosen to give the correct value of the the Z-boson mass. Only the first column of the Z matrix, Z α0, is defined by Eq. (.16. The remaining elements of Z must still be determined. Having fixed the first column of the matrix, the other three columns can be chosen to be orthogonal to the first column and to each other. This leaves us with an O(3 invariant subspace, such that the matrix Z is given by ( 1 0 Z = O, (.1 0 X 3 3 where O = Z 00 Z 10 + Z 0 + Z Z 10 Z 00 Z 10 Z10 +Z 0 +Z3 30 Z 0 Z 00 Z 0 Z10 +Z 0 +Z3 30 Z 30 Z 00 Z 30 Z10 +Z 0 +Z3 30 Z 0 +Z30 3 Z10 +Z 0 +Z3 30 Z 0 +Z3 30 Z 0 +Z3 30 Z 10 Z 0 Z10 +Z 0 +Z3 30 Z 10 Z 30 Z10 +Z 0 +Z Z 30 Z0 +Z3 30 Z 0 Z 0 +Z3 30, (. and X is an, as yet, undetermined 3 3 orthogonal [ matrix ( determined ] by three angles. This remaining freedom can be used to diagonalise Z T L ˆM Z, i.e. the (real symmetric ( ij sneutrino part of the Z T L ˆM Z matrix, with entries ( ˆM L i. We have now accomplished our aim of finding the matrices V and Z which, after inserting into potential of Eq. (.7 and dropping the primes, reduce the scalar potential to the form V neutral = (M L [ ] αβ ν Lα ν Lβ + m h0 h0 B α ν Lα h 0 + H.c (g + g ( h 0 h0 ν Lα ν Lα, (.3 where ( L with ˆM simple form (M L αβ [ ( ] Z T L ˆM Z αβ and B α (b Z α, (.4 and b given by Eq. (.8. In this basis the matrix M L adopts a particularly (M L αβ = B 0 tanβ 1 M Z cos β B j tanβ B i tanβ ( ˆM L i δ ij, (.5 where there is no sum over i in the down-right part of the matrix. Notice that we did not only succeed to self consistently go to a basis where the sneutrino VEVs are zero, but also 7

8 we managed to have the sneutrino masses ( ˆM L i diagonal and all the parameters of the scalar potential in Eq. (.3 real. As a byproduct of our procedure, we denote here that the potential of Eq. (.3 exhibits neither spontaneous nor explicit CP-violation at the tree level. The latter is in agreement with the results of Ref. [17] following a different method. Of course, the parameters µ α of the superpotential and the soft supersymmetry breaking couplings stay in general complex. A question arises when we include high order corrections to the potential. Then the vanishing sneutrino VEVs will be shifted to non-zero values by tadpoles originating, for example, from the LQD contribution in the superpotential (.1. The sneutrino VEVs maybe set back to zero by a renormalization condition such that a counterterm for these VEVs set their one particle irreducible (1PI tadpole corrections to zero. To conclude, it is worth to ( make a remark about the sign of B 0. As is clear from the form of Eqs. (.4,.8,.16, if L ˆM 1 M tan β 1 Z > 0 for all α, B tan αα β+1 0 is always positive in GH basis. 3 Parametrising the neutral scalar mass matrices The neutral scalar sector of the R-parity violating MSSM is in general very complicated. This is due to the fact that the scalars mix through the lepton number violating terms proportional to B i, ML and unless all of these parameters and VEVs are real one has a matrix to consider. However, for any given set of model parameters, one can always perform the basis change described in the previous section and arrive to the potential defined by Eq. (.3, with only real parameters. Consequently, the physical neutral scalars are, at the tree level, exact CP-eigenstates. This implies that the neutral scalar mass matrix decouples into two 5 5 matrices, one for the CP-odd particles and one for CP-even. Ultimately, one would like to parametrise the scalar sector resulting from the potential in (.3 with as few parameters as possible in order to make contact with phenomenology. These parameters in the case of the R-parity conserving MSSM are: the physical mass of the CP-odd Higgs boson M A = B 0 sin β, (3.1 and tanβ. An advantage of the form of potential in Eq. (.3,.4,.5 is that, M A and tan β can still be used for parametrising the general Higgs sector in the R-parity violating MSSM. M A is the mass of the lightest CP-odd Higgs boson in the R-parity conserving 8

9 MSSM; as such, it is used here as a parameter. m A is used to denote the physical treelevel mass of the lightest CP-odd Higgs in the R-parity violating MSSM (the convention adopted is that masses in the RPC case, parameters in this model, are denoted by M, and the masses in the RPV model are denoted by m. 3.1 CP-even neutral scalar masses and couplings The Lagrangian after spontaneous gauge symmetry breaking contains the terms ( Re h 0 L Re h 0 Re ν L0 Re ν Li M EVEN Re ν L0. (3. Re ν Lj As such, the scalar CP-even Higgs squared mass matrix becomes cos βma + sin βmz 1 sin β(m A + M Z B j M EVEN = 1 sin β(m A + M Z sin βma + cos βmz B j tanβ B i B i tanβ Mi δ ij, (3.3 where M i ( ˆM L i + 1 cos βm Z, (3.4 are in fact the sneutrino physical masses of the RPC case. It is important here to notice that the top-left sub-matrix is identical to the RPC case, for which the Higgs masses are given by M h,h = 1 (M Z + MA ± (MZ + M A 4MA M Z cos β and will be used as parameters in the RPV model., (3.5 The matrix (3.3 always has one eigenvalue which is smaller than MZ. This may be proved as follows: one first observes that the upper left submatrix of (3.3, call it A, has at least one eigenvalue smaller than or equal to MZ. Then using the Courant-Fischer theorem [18] of the linear matrix algebra one proves that, for one flavour, the eigenvalues of the 3 3 matrix M EVEN, are interlaced with those of A. This means that the matrix M EVEN with i = 1 has at least one eigenvalue smaller or equal than M Z. Repeating this procedure twice, proves our statement. Furthermore, it is interesting to notice that in the region where tan β 1, the eigenvector (sin β, cosβ, 0, 0, 0 T corresponds to the eigenvalue 9

10 with mass approximately MZ. Notice that this is the same eigenvector as in the RPC case which corresponds to the Higgs boson which couples almost maximally to the Z-gauge boson. Lepton flavour violating processes have not been observed as yet and therefore, bearing in mind cancellations, the parameters B i tan β have to be much smaller than min(m A, M i. To get a rough estimate, consider the dominant contribution from neutral scalars and neutralinos in the loop [19], m ν a ew B tan β 1 ev, (3.6 16π m 3 with m = max(m A, M i and B O(B i. This shows that B i tan β m m 0.1%. (3.7 With this approximation, it is not hard to find a matrix Z R which rotates the fields into the mass basis, such that ] Z T R M EVEN Z R = diag [m h 0, mh 0, (m ν+ i, (3.8 with m h being the lightest neutral scalar mass and 0 Z R = cos α sin α cos(β α cos αb j + sin(β α sinαb j cos β(mj M h cos β(mj M H sin α cos βp h i B i cos(β α cosα cos βp H i B i sin(β α cos(β α sin αb j + sin(β α cos αb j cos β(mj M h cos β(mj M H δ ij, (3.9 where there is no sum over i and (Mj, Mh, M H are defined in (3.4,3.5. In addition, tanα = tan β M A + M Z M A M Z and P h,h i = M Z cos β M h,h cos β(m H M h (M i M h,h, (3.10 (the common convention is to choose 0 β π/ and π/ α 0. The mass eigenstates of the RPV model are therefore given by ( cosβp h 0 cosαreh h 0 sin α Re ν L0 + i B i Re ν Li, (3.11 cos(β α ( cosβp H 0 sin α Reh H 0 + cosαre ν L0 + ( ν + i + ( cos(β α cosαb j i B i sin(β α Re ν Li, cosβ(mj M h + sin(β αb j sin αb j cosβ(mj M H ( cos(β α sin αbj cos β(mj M h + sin(β α cosαb j cos β(mj M H 10 Re h 0 Re ν L0 + Re ν Li,

11 with corresponding masses, m h 0 M h M Z cos β M h (M H M h cos β m H 0 M H + M Z cos β M H (M H M h cos β i=1 i=1 B i M i M h B i M i M H B 4 + O( M 6 cos 4 β, (3.1 B 4 + O( M 6 cos 4 β, (3.13 (m ν + i ( ˆM ν i + B i cos β B 4 Mi MZ cos β [ Mi 4 M i (M A + M Z + M A M Z cos β + O( M 6 cos 4 β. (3.14 The above expressions, are useful in relating the masses of the neutral scalars in the RPC and RPV case in the valid approximation B tanβ << min(m A, M i. They are presented here for the first time except the mass in (3.14 which agrees with Ref. [13]. We note here that these formulae are not valid if some of the diagonal entries in the mass matrix are closely degenerated - in such case even small B i terms lead to the strong mixing of respective fields. However in many types of calculations (e.g. various loop calculations one can still formally use such expansion - in the final result one often gets expressions of the type f(m 1 f(m m 1 m which have a well defined and correct limit also for degenerate masses, even if the expansion used in the intermediate steps was, in principle, wrong. It is interesting to note that the rotation matrix U defined in (.6, although explicitly calculated in this article, does not appear to all the neutral scalar vertices. For example, the vertices of the CP-even neutral scalars with the gauge bosons read as 1, L VVH = 1 g M Z cosθ w (cosβz Rs + sin βz R1s Z µ Z µ H 0 s + 1 g M W (cos βz Rs + sin βz R1s W +µ W µ H0 s, (3.15 where H 0 s=1,...,5 are the Higgs boson fields, h 0, H 0, ( ν + 1, ( ν +, ( ν + 3 respectively. From (3.9 and L VVH above, it is easy to see that the light Higgs boson coupling to the vector bosons (V = Z, W, is proportional to sin(β α as in the RPC case. In fact, the coupling 1 Note that the matrix Z defined in (.1 has nothing to do with neither Z R nor Z A defined in this section. We follow the conventions of Ref. [0]. ] 11

12 sum rule, 5 g V V H = s 0 g V V φ, (3.16 s=1 valid in the RPC case for s = 1,, persists also here, where g V V H 0 s are the couplings appearing in (3.15 and g V V φ the corresponding coupling appearing in the Standard Model. 3. CP-odd neutral scalar masses and couplings For the CP-odd case one finds, ( L Im h 0 Im ν L0 Im ν Li M ODD Im h 0 Im ν L0, (3.17 Im ν Lj where the CP-odd mass matrix reads, cos βma + ξ sin βm 1 Z sin β(m A ξm Z B j M ODD = 1 sin β(m A ξm Z sin βma + ξ cos βmz B j tan β, (3.18 B i B i tan β Mi δ ij and ξ is the gauge fixing parameter in R ξ gauge. In fact, by using an orthogonal rotation V = sin β cosβ 0 cos sin β , (3.19 we can always project out the would-be Goldstone mode, of the CP-odd scalar matrix and thus ξmz 0 0 V T M ODDV = 0 MA B j cos β. (3.0 0 Mi δ ij B i cos β Under the approximation of small bilinear RPV couplings [see Eq. (3.7], a solution is determined for the matrix Z A which rotates the fields into the mass basis, such that ] Z T A M ODD Z A = diag [m G 0, ma 0, (m ν i, (3.1 1

13 sin β Z A = cosβ 0 with the mass eigenstates given by cosβ sin β B i cos β(m i M A B j Mj M A B j tan β Mj M A δ ij, (3. G 0 sin β Im h 0 cosβ Im ν L0, A 0 cosβ Im h 0 + sin β Im ν L0 + ( ν i with corresponding masses, B i M j M A m A M A 1 cos β B i cosβ(m i M A Im ν Li, Im h 0 + B j tan β Mi Im ν L0 + Im ν Li, (3.3 M A i=1 B i M i M A B 4 + O( M 6 cos 4 β, (3.4 (m ν i M i B i (M A M i cos β + O( B 4 M 6 cos 4 β. (3.5 The coupling of the Z gauge boson to the CP-even and CP-odd neutral scalar fields is given by L ZHA = ig c W [ ( ] (p H 0 s p A 0 p µ Z R (+αs Z A (+αp Z R 1s Z A 1p Z µ Hs 0 A 0 p, (3.6 where the four momenta p µ H, p µ s 0 A are incoming and the fields A 0 0 p=1,...5 correspond to p G 0, A 0, ( ν 1, ( ν, ( ν 3 respectively. One may check that the coupling Z G 0 h 0 derived from (3.6 is proportional to sin(α β as it should be. 4 Positiveness and stability of the scalar potential 4.1 Positiveness In general, one should inspect whether all physical masses in the CP-odd and CP-even sector are positive. For that, all diagonal square subdeterminants of mass matrices should be positive. One can easily check that both CP-odd and CP-even mass matrices in (3.3,3.18 respectively, lead, in the rotated basis, to the same set of conditions, Mi > 0 with i = 1,, 3 and MA > 1 cos β 13 i=1 Bi Mi. (4.1

14 Using the form of MA in (3.1, the last equation can be rewritten in the form B 0 > tanβ i=1 Bi Mi. (4. Excluding some very singular mass configurations, the above conditions are rather trivially fulfilled if one takes into account the bound of Eq. ( Stability The question of whether the potential is stable, i.e. bounded from below, is far more complicated. In most cases the quartic (D term dominates and there is no problem. The only exception being when the fields follow the direction h 0 = 4 i=0 ν L i. In such a case, one should check whether the remaining part of the potential is positive along this direction. 3 Denoting R i=0 ν i and h 0 = Re iφ, where φ is a free phase, and using Eqs. (.0,.5,3.4, one can write down the scalar potential along this direction in the GH basis as B 0 V neutral = sin β cosβ ν L0 ν L0 + [ Mi + B 0 cot β ] ν Li ν Li + B i tanβ ( ν L0 ν Li + ν L0 ν Li B α ( νlα h 0 + H.c. ν L Q ν L ( B T ν L Re iφ + H.c.. (4.3 where the real symmetric matrix Q is ( Q = M A B i tan β B j tanβ [M i + B 0 cot β]δ ij. (4.4 Finding the stability conditions for the potential (4.3 is difficult, it depends on nine real variables (4 moduli and five phases of the fields. To simplify the problem, we perform one more field rotation to the basis in which the matrix Q is diagonal. This can be done, in general, by numerical routines (routines where already used in calculating the GH basis, and therefore, finding the stability conditions for the general scalar potential always has to involve some numerical analysis. We thus define the matrix P, ν L P ν L, as P QP = diag(x 0, X 1, X, X 3. (4.5 In fact, Q is real, so we can choose P to be real orthogonal. Also, we denote D β B α P αβ. Obviously, the rotation P preserves the value of R = h 0. 14

15 The potential becomes: V neutral = [ Xα ν Lα D α R ( ν Lα e iφ + H.c. ], (4.6 where X 0 has to be positive, otherwise for φ = 0 along the direction ν Li = Im ν L0 = 0 the potential V neutral = Re ν L0 [X 0 D 0 sign(re ν L0 ] falls to at least for one direction along the Re ν L0 axis. In fact the condition on X α is X α D α. Thus our first conclusion is that the matrix Q has to be positively defined. One can write down appropriate conditions in the same manner as for the scalar mass matrices; comparing with Eq. (4.1, it can be observed that this condition is automatically fulfilled if relation (4.1 holds. With X α positive, one can write down the potential as: V neutral = Xα ν Lα D α Re iφ R Xα D α X α. (4.7 To further simplify the problem, denote ν Lα = u α e i(φ φα, where u α 0 are field moduli and φ α are free phases. Then 3 ( V neutral = R u α Xα R D α e iφα Xα 3 D α X α, (4.8 where R = i=0 ν i = i=0 u i. Phases φ α can be adjusted independently of u α. The worst case from the point of view of potential stability, the smallest first term inside the parenthesis, occurs for D α e iφα = D α. Denoting further ǫ α = u α /R, 0 ǫ α 1, one can reduce our initial problem to the question whether the function g(ǫ α = Xα ǫ α Dα Xα D α X α = (X α ǫ α D α ǫ α, (4.9 depending now on four real positive parameters, is non-negative on the unit sphere 3 ǫ α = 1. In general such problem can be solved numerically using the method of Lagrange multipliers. For X i > X 0 D 0, the minimum occurs for ǫ α = D α X α + λ, (4.10 where λ can be found numerically as a root of the following equation: D α (X α + λ = 1. (

16 For smaller X i, the minimum is realized for ǫ i = 0 for one or more values of i and requires analysis of various special cases. Having found the correct minimum, to prove the stability of the potential one needs to show that the function g at the minimum is non-negative. As shown in Eq. (3.7, B i terms and thus also D i terms are usually very small. In this case one can set approximate, sufficient conditions for the stability of the potential, without resorting to solving Eq. (4.11, numerically. Denote D = 3 i=1 D i and X min = 3 3 min(x 1, X, X 3. Then, using the inequality D i ǫ i i=1 D i i=1 ǫ i = D 1 ǫ 0, one has g(ǫ α X 0 ǫ 0 + X min (1 ǫ 0 + (X i X min ǫ i D 0 ǫ 0 D 1 ǫ 0. (4.1 Terms (X i X min ǫ i are always non-negative. The worst case being when the vector (ǫ 1, ǫ, ǫ 3 is along the minimal X i axis, where these terms vanish. Other terms are rotation invariant in the 3-dimensional space (ǫ 1, ǫ, ǫ 3, so Eq. (4.1 is equivalent to finding parameters X 0, X min, D 0, D for which the expression (4.13, depending on just one real variable, is positive: g (ǫ 0 = X 0 ǫ 0 + X min(1 ǫ 0 D 0 ǫ 0 D 1 ǫ 0 0. (4.13 Analysis of (4.13 is further simplified by one more approximation, justified for small D: g (ǫ 0 X 0 ǫ 0 + X min(1 ǫ 0 D 0 ǫ 0 D. (4.14 The rhs of Eq. (4.14 is now trivial. Following approximate conditions for the stability of the potential can be summarized as follows: X min range X min X 0 D 0 Stability requires X 0 D 0 + D 0 < X min < X 0 D 0 (X 0 X min (X min D D0 Both conditions are sufficient, but not minimal - we have made some approximations and there may be parameters which do not fall into either of the categories above, and yet still give a stable potential. For example, if X 0 = X 1 = X = X 3 X, one can easily derive the exact necessary and sufficient condition for potential stability as X D0 + D, less strict than X ( D 0 + D which would be given by the table above. For complementary work the reader is referred to Ref. [1]. 16

17 5 Conclusions In this letter we present a procedure for calculating the rotation matrix which brings the neutral scalar fields of the general R-parity violating MSSM onto the GH basis where they develop n zero VEVs, with n being the number of flavour generations. In doing so, we have made no assumption about the complexity of the parameters. We consider the case of n = 3 generations, but our approach immediately applies to other cases, apart from obvious modifications of the form of Z matrix defined in (.1... As a byproduct of basis change, we prove that the tree level MSSM potential does not exhibit any form of CP-violation, neither explicit nor spontaneous. Consequently, the neutral scalar fields can be divided into CP-even and CP-odd sectors with the 5 5 neutral scalar squared mass matrices, taking a very simple form with only small RPV masses sitting on their off diagonal elements. We can thus expand along small RPV masses and find analytic approximate formulae which relate the RPC and the RPV neutral scalar masses. Furthermore we also find, that in general there is always at least one neutral scalar field with mass lighter than M Z which couples maximally to the Z-gauge boson in the case of large tanβ and large M A. Our procedure for finding the rotation matrix U has been coded 3 and is numerically stable. In the end, we are aiming to construct the most general MSSM quantum field theory structure resorting neither to R-parity violation nor to other approximations. This will be useful for examining the phenomenology of the MSSM as a whole. The convenient choice of the basis for the neutral sector found in this paper is a first step towards this direction. Acknowledgments A.D. and M.S.-S. would like to thank The Nuffield Foundation for financial support. J.R. would like to thank IPPP for the hospitality during his visit. His work was supported in part by the KBN grant P03B ( S.R. acknowledges the award of a UK PPARC studentship. References [1] G. Feinberg, P. Kabir and S. Weinberg, Phys. Rev. Lett. 3 ( See also S. Weinberg, The Quantum Theory Of Fields. Vol. I: Foundations Cambridge, UK: Univ. Pr. ( The code will be available from rosiek/physics/rpv/scalar.html 17

18 [] H. P. Nilles, Phys. Rept. 110 (1984 1; H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985; S. P. Martin, arxiv:hep-ph/ [3] P. Fayet, Phys. Lett. B 69 ( ; G. R. Farrar and P. Fayet, Phys. Lett. B 76, 575 (1978. [4] For reviews, see H. K. Dreiner, arxiv:hep-ph/ ; G. Bhattacharyya, Nucl. Phys. Proc. Suppl. 5A ( [arxiv:hep-ph/ ]; M. Chemtob, Prog. Part. Nucl. Phys. 54 ( [arxiv:hep-ph/040609]. [5] Y. Grossman and H. E. Haber, Phys. Rev. D 67 ( [arxiv:hep-ph/01073]. [6] S. Weinberg, Phys. Rev. D 6 ( Phenomenology of the L = 1 terms LLĒ, LQ D and LH were first considered in L. J. Hall and M. Suzuki, Nucl. Phys. B 31 ( Phenomenology of the B = 1 term Ū D D in the superpotential was first considered in F. Zwirner, Phys. Lett. B 13 ( [7] T. Banks, Y. Grossman, E. Nardi and Y. Nir, Phys. Rev. D 5 ( [arxiv:hep-ph/950548]. [8] R. Hempfling, Nucl. Phys. B 478 ( [arxiv:hep-ph/951188]. [9] M. Nowakowski and A. Pilaftsis, Nucl. Phys. B 461 ( [arxiv:hep-ph/950871]. [10] F. de Campos, M. A. Garcia-Jareno, A. S. Joshipura, J. Rosiek and J. W. F. Valle, Nucl. Phys. B 451, 3 (1995 [arxiv:hep-ph/95037]; M. Hirsch, M. A. Diaz, W. Porod, J. C. Romao and J. W. F. Valle, Phys. Rev. D 6 ( [Erratum-ibid. D 65 ( ] [arxiv:hep-ph/ ]. [11] B. C. Allanach, A. Dedes and H. K. Dreiner, Phys. Rev. D 69 ( [arxiv:hep-ph/ ]. [1] D. W. Jung, S. K. Kang, J. D. Park and E. J. Chun, JHEP 0408 ( [arxiv:hep-ph/ ]. [13] Y. Grossman and H. E. Haber, Phys. Rev. D 59, (1999 [arxiv:hep-ph/ ]. [14] S. Davidson, M. Losada and N. Rius, Nucl. Phys. B 587 ( [arxiv:hep-ph/ ]. 18

19 [15] Y. Grossman and H. E. Haber, Phys. Rev. D 63 ( [arxiv:hep-ph/000576]. [16] S. Davidson and H. E. Haber, arxiv:hep-ph/ [17] M. Masip and A. Rasin, Phys. Rev. D 58 ( [arxiv:hep-ph/980371]. [18] See for instance, Carl D. Meyer, Matrix Analysis and Applied Linear Algebra, Siam, 000; R. A. Horn and C. R. Johnson, Topics in Matrix analysis, Cambridge University Press, [19] S. Davidson and M. Losada, JHEP 0005 (000 01[arXiv:hep-ph/ ]; A. Abada, S. Davidson and M. Losada, Phys. Rev. D 65 ( [arxiv:hep-ph/011133]; M. A. Diaz, M. Hirsch, W. Porod, J. C. Romao and J. W. F. Valle, Phys. Rev. D 68 ( [Erratum-ibid. D 71 ( ] [arxiv:hep-ph/03001]; Y. Grossman and S. Rakshit, Phys. Rev. D 69 ( [arxiv:hep-ph/ ]. [0] J. Rosiek, Phys. Rev. D 41 ( ; Erratum arxiv:hep-ph/ [1] S. A. Abel and C. A. Savoy, Nucl. Phys. B 53 ( [arxiv:hep-ph/980318]. 19