Vacuum structure of the left-right symmetric model

The left-right symmetric model (LRSM), originally proposed to explain parity violation in low energy processes, has since emerged as an attractive framework for light neutrino masses via the seesaw mechanism. The scalar sector of the minimal LRSM consists of an $SU(2)$ bi-doublet, as well as left- and right-handed weak isospin triplets, thus making the corresponding vacuum structure much more complicated than that of the Standard Model. In particular, the desired ground state of the Higgs potential should be a charge conserving, and preferably global, minimum with parity violation at low scales. We show that this is not a generic feature of the LRSM potential and happens only for a small fraction of the parameter space of the potential. We also analytically study the potential for some simplified cases and obtain useful conditions (though not necessary) to achieve successful symmetry breaking. We then carry out a detailed statistical analysis of the minima of the Higgs potential using numerical minimization and find that for a large fraction of the parameter space, the potential does not have a good vacuum. Imposing the analytically obtained conditions, we can readily find a small part of the parameter space with good vacua. Consequences for some scalar masses are also discussed.


I. INTRODUCTION
The discovery of neutrino masses is a sure sign of new physics beyond the Standard Model (SM). A simple paradigm for neutrino masses is the seesaw mechanism [1][2][3][4][5] which introduces right-handed neutrinos (RHN) with heavy Majorana masses. Two questions then arise: (i) what is the seesaw scale or the mass of the RHNs? and (ii) what is the ultraviolet (UV)-complete theory that leads naturally to the basic ingredients inherent in the seesaw mechanism i.e. RHNs and a B − L symmetry whose breaking gives rise to their Majorana masses? Two classes of theories with this property are: (i) the SO(10) model whose basic spinor representation contains the RHN and which contains a group generator that is the B − L local symmetry [6], and (ii) the left-right symmetric model (LRSM) [7][8][9], which is the simplest extension of the SM that contains three RHNs to cancel the gauge anomalies and B − L symmetry as a natural symmetry [10].
In this paper, we focus on symmetry breaking aspects of the minimal LRSM and carry out an analysis of its vacuum structure. The general procedure to investigate this is to write down the Higgs potential involving the various scalar multiplets of the LR gauge group SU (2) L × SU (2) R × U (1) B−L and look for the minimum of the potential that breaks the gauge symmetry down to U (1) em . Morever, the theory should be parity violating at low scale and generate naturally small neutrino masses. The detailed analysis of the non-supersymmetric LRSM Higgs potential and its minima have been discussed in many works [11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26] (supersymmetric LRSM Higgs sectors have been studied in Refs. [27,28]). It is known that for certain ranges of the parameters (e.g. negative scalar mass squares and positive values for scalar couplings), a desired (good) vacuum is obtained. However, if all the couplings are chosen randomly to start with, it is not clear how often one gets a good vacuum. Secondly, it is not known whether those minima obtained above are global minima of the potential or are simply the local ones. Third, the boundedness-from-below of the potential has been used as the necessary and sufficient condition for vacuum stability [17,18], but as we show in this paper, a bounded-from-below potential is necessary for a good vacuum, but not sufficient.
We use the gauge freedom of the theory to give simple criteria that are to be fulfilled in order to end up in a good vacuum (i.e. charge conserving and parity violating at low scales). Keeping arbitrary values for the parameters of the Higgs potential we furthermore check for what fraction of the parameter space a global minimum with desired properties is obtained. Once a range of the parameters is determined where a global minimum occurs, one can then use them to find the scalar spectrum corresponding to that choice, which is a potential test of the model. We do not carry out an exhaustive analysis of the scalar masses but rather give some simple examples at viable minima of the model.
The paper is organized as follows. In Sec. II, we review the details of the model and its scalar sector. In Sec. III, we write down the full scalar potential and give criteria for obtaining the good vacua of the model. In Sec. IV we present analytical studies focused on some simplified cases, followed by numerical studies to scan the whole parameter space in Sec. V. Finally we summarize and conclude in Sec. VI. Some technical details are delegated to the appendix.
It is conventional to adopt the matrix representation in which φ and ∆ L, R are written as The corresponding transformation rules are for U L ∈ SU (2) L , U R ∈ SU (2) R and e iθ B−L ∈ U (1) B−L . Note thatφ ≡ σ 2 φ * σ 2 = − φ * transforms in the same way as φ.
The model also has a discrete left-right symmetry, which can either be the P parity or the C parity: The scalar potential with P parity is more constrained than that with C parity 1 as the latter allows several complex phases. In this paper, for simplicity, we assume all the couplings in the scalar potential to be real, i.e. there is no explicit CP violation 2 in the potential. Such a potential respects both parities. In general the full potential contains 17 gauge invariant terms [14]. After spontaneous symmetry breaking, some components of φ and ∆ L,R obtain nonzero vacuum expectation values (VEVs) while the others do not, depending on the parameters of the potential. Since the gauge symmetry SU (2) L ⊗ SU (2) R ⊗ U (1) B−L is required to break to U (1) em by these scalar fields, the desired VEV alignment is [14] φ = 1 √ 2 1 See Eqs. (9) and (10) in Ref. [24] for comparison. 2 We do, however, find that a small portion of the randomly generated samples in Sec. V have complex VEVs though the potential parameters are real, which implies that spontaneous CP violation is possible in this model. See Refs. [11][12][13][14][15][16] for more details.
The VEVs should furthermore obey the hierarchy v L κ 1, 2 v R (v L may vanish) to meet the known phenomenology, such as tiny neutrino masses, heavy RH gauge boson masses, the electroweak precision parameter ρ 1, etc.
Although Eq. (10) is what we need to successfully achieve spontaneous symmetry breaking in the LRSM, for general (arbitrary) values of the parameters, the scalar potential does not necessarily lead to this VEV alignment. For example, if we minimize the scalar potential we may obtain a minimum with nonzero diagonal VEVs of ∆ L or ∆ R , which would break U (1) em . It is also possible to get a minimum with ∆ L = ∆ R which would imply unbroken parity symmetry. The various possibilities of symmetry breaking with the full scalar potential of LRSM, due to the considerable complexity, have never been comprehensively studied before.
In this paper, we will therefore address an essential question of the spontaneous symmetry breaking in the LRSM: How can we obtain the VEV alignment in Eq. (10) and how likely is this?
Since the full scalar potential is very complicated, a purely analytical study is difficult and we mainly adopt a numerical approach. However, we provide some illustrative analytical studies for simplified cases, where a lot of terms in the potential are absent. Nevertheless, the analytical results give us some useful insight into the vacuum structure of the scalar potential and serve as a supplement to the full numerical calculations. Our numerical approach has already been established in Refs. [29][30][31] to successfully analyze beyond the SM scalar potentials. In general, given specific values of the potential parameters, we can always use a computer program to numerically minimize the potential and obtain a minimum. With further developed algorithms (see the details presented in Sec. V), we can make the program capable of identifying the zero entries in Eq. (10). In this way we can find out all possible VEV alignments that can be obtained in the LRSM potential. We choose not to use the Vevacious package [34], which, among other things, can also provide the minima of beyond the SM scalar potentials. This package currently has limited capability to the case we are interested in, since our potential contains many parameters and field components. In addition, we will perform a statistical analysis with a large number of random samples. Hence, we use a self-written dedicated minimization program, which we have made publicly available in GitHub [35].
As noted earlier, the significance of such a study is two-fold. First of all, for any given set of potential parameters, we can infer whether it can lead to successful symmetry breaking and whether the minimum of the potential is a global minimum. This in turn can put constraints on the potential parameters and also on the scalar mass spectrum. If the LRSM is taken as a serious theory of particle interactions beyond the SM, then these theoretical constraints from vacuum stability should be taken into consideration, in combination with other theoretical constraints, such as unitarity and perturbativity [22,23,36], as well as experimental constraints from lepton flavor violation, neutrinoless double beta decay, rare meson decays and colliders [37][38][39][40][41][42][43][44][45][46][47][48][49][50][51][52][53][54].
Secondly, for extended LRSMs with modified scalar sectors, e.g. with additional scalar bi-doublets, triplets, singlets, etc. [55][56][57][58][59][60][61][62][63][64][65][66], one would again be concerned about the question of whether the desired VEVs can be obtained. While the potential for such cases may be too complicated to repeat the analytical calculations given here, our numerical method can be easily implemented to analyze the vacuum structures of such models. We leave these studies for future work.

III. THE SCALAR POTENTIAL
The most general gauge invariant scalar potential invariant under Eq. (6) contains 17 independent terms: where, as we mentioned above, all couplings are assumed real. The spontaneous symmetry breaking in such a complicated potential may result in various types of VEV alignments. Some of them may successfully achieve the desired symmetry breaking given by Eq. (10) and have phenomenologically viable consequences whereas some may not. In what follows, for convenience, we will refer to the former as good vacua and the later as bad vacua.

A. Good Vacua
Since the electromagnetic gauge symmetry U (1) em should not be broken in any extension of the SM, only the electric neutral components (namely φ 0 1 , φ 0 2 , δ 0 L and δ 0 R ) can acquire nonzero VEVs. In general, the charge conserving VEVs of φ, ∆ L and ∆ R should be However, two of the phases (let us take θ 1 and θ R ) can be removed by the transformation (6) with Meanwhile, θ 2 and θ L are transformed to θ 2 + θ 1 and θ L + 2θ 1 − θ R , which for simplicity can be redefined as θ 2 and θ L . Therefore, without loss of generality, one can always set which reduces Eq. (12) to Eq. (10). Assuming the potential has a minimum for the VEVs given by Eq. (10), one can replace the fields with their VEVs and derive the minimization conditions: From Eq. (15) one can derive (see Appendix A) the renowned seesaw relation of VEVs in LRSM [67]: The left-hand side is roughly of the order βv 2 where β = β 1, 2, 3 and v = 246 GeV, and the right-hand side is ρ v L v R where ρ ≡ 2ρ 1 − ρ 3 . From βv 2 = ρ v L v R , one can see that for very large v R (correspondingly very heavy W R ), v L will be suppressed by 1/v R , corresponding to very tiny neutrino masses, known as the seesaw relation of the VEVs.
Here we would like to give two comments.
, then the VEVs (κ 1 , κ 2 , v L , and v R ) may violate the relation (16) while the derivatives in Eq. (15) remain zero. This can be seen from the analytical calculations in Appendix A and is also verified in our numerical studies. Despite that v 2 L = v 2 R is phenomenologically not allowed (since parity is not broken in this case), it turns out that this case appears much more frequently in the numerical scan than the case with v 2 L = v 2 R . Therefore, if in numerical studies a minimum is obtained with the same VEV alignments as in Eq. (12), one should carefully check whether v 2 it is a good vacuum with the VEVs satisfiying the seesaw relation (16).
• Eq. (15) is based on the assumption that a minimum in the form of Eq. (12) exists. The six equations in (15) can (in general) always be solved with respect to the six variables (κ 1 , κ 2 , v L , v R , θ 2 , θ L ). However, the existence of solutions of Eq. (15) implies neither the existence of the minimum, nor that the first-order derivatives with respect to the fields vanish, i.e. ∂V /∂ϕ = 0, where ϕ stands for all components of φ and ∆ L,R .

B. Bad Vacua
Next, we shall investigate other vacua that could appear but would lead to unacceptable physical or phenomenological consequences, e.g. U (1) em breaking. Without the requirement of charge conservation, in general, any components of φ and ∆ L,R could acquire nonzero VEVs. But before we put nonzero VEVs arbitrarily, we need to examine the symmetries in the potential to avoid considering redundant cases.
First of all, there is the gauge symmetry SU (2) L ⊗ SU (2) R ⊗ U (1) B−L , which allows one to remove some degrees of freedom (DOF) in the scalar fields by gauge fixing. By analogy to the SM case, where the Higgs doublet in the unitarity gauge has only one DOF (the physical Higgs boson), in the LRSM we can adopt a similar gauge to remove 3 + 3 + 1 = 7 DOFs, equal to the number of gauge bosons. More specifically, we use the transformations U L and U R in Eq. (6) Some phases can be further removed by Eq. (13) and U (1) B−L transformations 3 so that one can further set As one can check, indeed seven DOFs have been removed. Secondly, a vacuum with φ 1 = φ 2 and nonzero b 1, 2 or c 1, 2 does not necessarily break U (1) em due to additional symmetries in the vacuum. Note that if where U L can be any SU (2) matrix. Accordingly, ∆ L and ∆ R will be transformed by the above U L and U R to other forms. For example, if θ 2 = 0 and the following identical textures of ∆ L,R are realized, then one could choose the following U L,R transformations: It is now straightforward to show that one obtains the same VEV alignment as in Eq. (10). This implies that the vacuum in this example is physically equivalent to the charge-conserving vacuum. A necessary and sufficient condition to infer whether ∆ L and ∆ R really break U (1) em in the absence of φ is that U (1) em is not broken if and only if It is straightforward to see that if U (1) em is not broken, then the determinants must be zero. The converse, however, needs a short proof: first, note that in the absence of φ , for any ∆ L and ∆ R , we can always transform them via the Schur decomposition 4 to the following form: The determinants det ∆ L = x 2 L and det ∆ R = x 2 R imply that if they are zero, the diagonal elements in Eq. (23) must be zero, hence However, in the presence of φ such transformations may be partially or fully forbidden, depending on whether φ 1 = φ 2 or not. In this case, it could be that ∆ L and ∆ R with zero determinants break U (1) em . But in our numerical study presented in Sec. V, among a large number of randomly generated samples, we do not find any samples belonging to this exotic category. Therefore, based on a high-statistics numerical study, we can draw the conclusion that in the presence of φ , generally det ∆ L = det ∆ R = 0 is sufficient to ensure the conservation of U (1) em .
Finally, some discrete symmetries may connect the vacuum in Eq. (10) to other vacua. Consider the following two vacua: As one can check, the above two minima can be transformed to each other by U L = U R = iσ 2 combined with the P parity transformation: They can also be transformed to Eq. (10) by either U L = U R = iσ 2 or P. If the potential has a minimum at Eq. (10), then the above minima (24) and (25) Combining the transformations of both U L = U R = iσ 2 and parity, one can get one more vacuum with the same potential depth.
Therefore, the vacuum of Eq. (10), if it exists, is always accompanied by several wrong vacua which have the same potential depth and are connected by discrete symmetries 5 , which are subgroups of the left-right gauge symmetry. On the other hand, if one finds a minimum corresponding to one of these wrong vacua, then it implies the existence of the true vacuum. In this sense, searching for minima of these wrong types is also useful. This is particularly important for the numerical searches to be performed later.
In summary, a vacuum which superficially seems to be bad may actually imply the coexistence of a good vacuum, or may itself be a good vacuum up to some continuous symmetry transformation. Considering these possibilities, we would like to propose the following gauge independent criteria for the good vacuum: good vacuum criteria : For the VEV alignment in Eq. (10), one can straightforwardly check that the above criteria are satisfied. Conversely, if a minimum of the potential satisfies Eq. (27), the VEVs must be of the form in Eq. (10) or can be transformed to Eq. (10) under the previous mentioned symmetries.

IV. ANALYTICAL STUDY OF LR VACUA IN LIMITING CASES
In principle, we can analytically compute the first-order derivatives of the scalar potential to find out the minima. However, for the full potential in Eq. (12) the calculations are too complicated to perform analytically. In this section, we focus on some simplified cases in which several terms in the potential are absent. Although as such it is not a full analysis (and sometimes even unrealistic), the analytical results obtained in this way provide crucial insight into the vacuum structure of the scalar potential. The details are expanded below and the results are summarized at the end of this section.
The first simplification we will make is to set α and β to be zero because only the αand β-terms "lock" the bidoublet with the triplets. With α = β = 0, the potential can be written as two separate parts where V φ and V ∆ contain only φ and ∆ L/R respectively. In this case, the global symmetry of the potential becomes much larger i.e.
(In fact this is true as long as α 3 = 0 and β i = 0) Once the φ and ∆ vevs are switched on, the resulting symmetry will be U (1) This will lead to 10 massless states out of which six will be absorbed as longitudinal modes of the gauge bosons of the theory leaving four massless states. Clearly therefore, this theory is not realistic. But nevertheless we study the vacuum structure and the symmetries of the vacuum in order to understand the same for the full model with all couplings turned on. Expressed in terms of the explicit components defined in Eq. (17), V φ and V ∆ are Here the field components (e.g., φ 1, 2 , a 1, 2 , b 1, 2 , c 1, 2 , · · · ) have been defined in Eq. (17). To proceed, let us further set µ 2 2 , λ 2 , λ 4 , ρ 2 , and ρ 4 to zero so that we are not bothered by the cosines appearing in the above expression. In this very simplified case, we note that V φ and V ∆ have essentially the same form: Here we have extracted 2ρ 1 δ 2 1 δ 2 2 in Eq. (32) so that V ∆0 shows explicitly the same form as V φ0 . In the following discussion, we will focus on V φ0 while the conclusions can be easily transferred to V ∆0 by replacing µ 2 1 → µ 2 3 , λ 1 → ρ 1 and 4λ 3 → (ρ 3 − 2ρ 1 ). A notable feature of V φ0 is that it respects the following dihedral (D 4 ) symmetry which leads to D 4 -symmetric vacuum structures shown in Fig. 1. The bounded from below (BFB) condition for the potential in Eq. (31) is manifest: In the following discussion of minima, by default we assume the BFB condition should be satisfied. From the first-order derivatives we get three possible solutions Here all the denominators such as 4(λ 1 + λ 3 ) and 2λ 1 are positive according to the BFB condition (35), which implies that the solutions (38) and (39) will not exist if µ 2 1 < 0, because this will lead to negative φ 2 1 or φ 2 1 . Therefore, for µ 2 1 < 0, the potential has only one minimum, which is necessarily the global minimum of the potential. By comparing the potential values V φ0 of the three solutions, we can infer which one can be the global minimum without computing the second-order derivatives (Hessian matrix). One can see from Eq. (38) and (39)  2λ1 is a global minimum. In summary, assuming µ 2 1 > 0 and λ 1 > 0, the vacuum structure depends on the ratio λ 3 /λ 1 as follows: 2λ1 , with equal depth; • λ 3 /λ 1 > 0: V φ0 has four minima at (φ 1 , φ 2 ) ∝ (0, ±1) and (±1, 0), with equal depth.
So far, we have not taken µ 2 2 , λ 2 , λ 4 into consideration. In the presence of these terms, the D 4 symmetry will be broken and the complexity of the above analysis will be increased. We do not plan to derive the corresponding analytic expressions. But we can discuss qualitatively on the consequences using the D 4 -symmetric conclusions. If we set µ 2 2 > 0 while keeping λ 2 and λ 4 still at zero, this is the same as adding the term −4µ 2 2 φ 1 φ 2 cos θ 2 to the D 4 -symmetric potential V φ0 . From the potential V φ0 − 4µ 2 2 φ 1 φ 2 cos θ 2 we can immediately see that cos θ 2 should be 1 (if 4µ 2 2 φ 1 φ 2 > 0) or −1 (if 4µ 2 2 φ 1 φ 2 < 0) to reach a minimum, which implies that we would still get real solutions even if we turn on the complex phase. Therefore instead of V φ0 − 4µ 2 2 φ 1 φ 2 cos θ 2 we can focus on V φ0 − 4µ 2 2 φ 1 φ 2 . In  , the four minima has the same depth due to the D 4 symmetry. When the µ 2 2 term is added (for illustration we choose µ 2 2 = 0.15), as shown in the corresponding lower panel, two of the minima become deeper than the other two. Note that in this case (φ 1 , φ 2 ) at the minima still align in the direction (1, 1) or (1, −1). In the middle panels (λ 3 /λ 1 = 0), the SO(2) symmetry is broken when µ 2 is nonzero, leading also to the VEV alignment (φ 1 , φ 2 ) ∝ (1, 1). In the right panels (λ 3 /λ 1 = 1), the four minima still have equal depth after adding the µ 2 2 term, but the VEV alignment is changed from (1, 0) or (0, 1) to (1, r) or (r, 1) where r depends on µ 2 . Therefore, within the simple D 4 -soft-broken potential V φ0 − 4µ 2 2 φ 1 φ 2 cos θ 2 , we can already get an arbitrary VEV alignment (φ 1 , φ 2 ) ∝ (1, r). Further turning on λ 2 and λ 4 couplings may produce more possibilities (e.g. spontaneous CP breaking) which should include the VEV alignments obtained in V φ0 − 4µ 2 2 φ 1 φ 2 cos θ 2 . Now let us discuss on the triplet sector. As previously mentioned, if ρ 2 = ρ 4 = 0, the potential of ∆ L and ∆ R reduces to V ∆0 which has the same form as V φ0 . Using the previous conclusions on V φ0 , we know that when the D 4 -symmetric potential V ∆0 has four minima at (δ 1 , δ 2 ) ∝ (0, 1) or (1, 0), which implies the VEV alignment can be obtained. We would like to point it out that even if ρ 2 and ρ 4 are nonzero, V ∆ is still D 4 -symmetric. The symmetry transformation is similar to Eq. (33) with (φ 1 , φ 2 ) replaced by (∆ L , ∆ R ). Consequently, the VEV alignments of (∆ L , ∆ R ) from Eq. (31) have only three possibilities: (i) one of (∆ L , ∆ R ) is nonzero and the other is zero, like Eq. (42); (ii) both are zero; or (iii) both are nonzero but ∆ L = ∆ R . However, the above arguments treat ∆ L and ∆ R as two singlets and do not take the fact that they have internal components into consideration. Analyses in terms of the components (a i , b i , c i ) will be much more complicated. We adopt a numerical method to check the above conclusions and find that there are indeed only the three cases, except that (iii) should be interpreted as a 2 1 + b 2 1 + c 2 1 = a 2 2 + b 2 2 + c 2 2 . Besides, one should notice that for nonzero ρ 2 and ρ 4 , Eq. (41) is no longer the condition to get Eq. (42). But at least the simple conclusion with ρ 2 = ρ 4 = 0 is enough to show that the VEV alignment in Eq. (42) can be obtained in a part of the parameter space.
Eq. (42) is what we need to achieve spontaneous parity breaking. The other two possible cases, namely ∆ L = ∆ R = 0 and ∆ L = ∆ R = 0, can be modified by D 4 -breaking terms [e.g. the α and β terms in (12)] so that ∆ L = ∆ R and ∆ R = 0. However the experimental constraints require that ∆ L should be much smaller than ∆ R since the former, limited by the electroweak ρ parameter, has to be lower than the GeV scale while the latter should be above the TeV scale (collider bounds on W R ). Such a strong hierarchy would require substantial fine-tuning in the scalar potential. Therefore, the solution ∆ L = 0 as a consequence of symmetry is a more favored option.
Focused on the solution with zero ∆ L and nonzero ∆ R , we proceed to study the VEV alignments of the internal components of ∆ R . For simplicity, let us first set ρ 4 = 0 (µ 2 3 , ρ 1 , ρ 2 , and ρ 3 are nonzero). Note that in the remaining terms of V ∆ , only the ρ 3 term couples ∆ L to ∆ R . This term has no contribution to the first-order derivatives at the minimum with zero ∆ L because In the absence of ρ 3 and ρ 4 , ∆ R decouples with ∆ L in V ∆ . Hence we only need to consider the following part of the potential where in the second line we have simplified the ρ 2 term which can be verified by an explicit computation. To proceed, we need the following useful relation: which can be proven by where ω = α 2 − 2β 2 + γ 2 . Note that for a fixed value of δ 2 2 , | det ∆ R | can reach any value in the above range. Therefore, we can parametrize | det ∆ R | as 1 2 δ 2 2 cos θ and write From Eq. (45), it is obvious to identify the minima. Let us take ρ 1 > 0 and µ 2 3 > 0, which is necessary to satisfy the BFB condition and obtain nonzero VEVs. If ρ 2 > 0, the minimum should be at cos 2 θ = 0; and if ρ 2 < 0, it should be at cos 2 θ = 0. Therefore the minima of V ∆2 should locate at: As long as α 1, 2, 3 = β 1, 2, 3 = 0, we always have the freedom to transform ∆ R → U R ∆ R U † R individually (without the corresponding transformation of the bidoublet φ) within the triplet potential V ∆ . According again to the Schur decomposition, we can always transform ∆ R to a lower triangular matrix. In this form, Eqs. (46) and (47) should be Eq. (48) is straightforward to get, because if det ∆ R = 0, only the 2-1 element can be nonzero. Eq. (49) has zero a 2 because when | det ∆ R | = 1 2 δ 2 2 , |a 2 | and |c 2 | should be equal, according to the derivation of Eq. (44). One should keep in mind that Eqs. (48) and (49) are derived under the assumption that µ 2 3 > 0, ρ 4 = α 1, 2, 3 = β 1, 2, 3 = 0, and that the potential is BFB.
The above analyses implies the following sufficient but not necessary conditions to get a good vacuum: If the potential parameters satisfy the above conditions, then it can be guaranteed that the potential has a global minimum corresponding to a good vacuum. In practical use, however, α i and β i cannot all be set to zero because of additional massless states, as mentioned before. To solve this problem, we can add small perturbations to α i to avoid the massless states. More explicitly, using the above condition, we can easily find a set of potential parameters satisfying Eqs. (50), (51) and (52). Then if we set α i to zero, the potential has a good vacuum, though not realistic. Next we can explore the nearby around this point by tentatively adding some small perturbations to α i . If the perturbations are small enough, the conclusion should hold as well. Sometimes, the perturbations can be very large without changing the conclusion. The exploration starting from Eq. (53) needs numerical assistance, as will be done in Sec. V. And we will show (see Fig. 5), indeed one can find some deviations from Eq. (53) that lead to successful symmetry breaking. Besides, we have other comments on the above conditions: • The above conditions also guarantee BFB; • The vacuum obtained in this way always has ∆ L = 0.
The conditions (50)-(53) obtained by the above analytic study will be very useful in the subsequent numerical study. It helps to quickly find out a viable region in which the potential has a good vacuum. In addition, it is also important for setting some benchmarks when studying the global minimum constraints.

V. NUMERICAL STUDY
With the modern technology of numerical computation, given a set of specific values of the potential parameters one can readily obtain a numerical minimum. There have been various algorithms well developed to find minima of multi-variable functions 6 . Most algorithms are based on iterative searches which means the program starts by given an initial point and iteratively computes the next step according to some principles until the steps converge to a minimum. The convergence can not be guaranteed, so for a single process searching for minima there is a small probability of failure. In case of failure, one can try to start the process again with a different initial point. Repeating the processes will eventually arrive at a minimum which may be local or global.
Once we get a numerical minimum of the potential, we can inspect the field component values at the minimum, checking if some of them are zero or some of them are equal. However, in numerical calculations it is impossible to have infinite accuracy so the would-be zero numbers are generally nonzero but very small (e.g. 10 −6 ∼ 10 −9 for numerical minimization based on a 8-byte real number system). This implies that, due to limited accuracy, there is no absolute equality in the numerical results. The simplest solution to this problem is by setting a cut on the difference of two numbers, below which the two numbers are thought to be equal and above which they are not 7 . But this may   cause misjudgment since the difference of two actually unequal numbers may occasionally be smaller than the cut. To double check if a field component is zero at the minimum, we can invoke the minimization process again with this component fixed to zero. If it gives a slightly better minimization, then we can conclude that it is zero at the minimum. We apply the above numerical method to the LRSM potential in Eq. (12). First, let us arbitrarily set the values of potential parameters to see in general what VEV alignments would be obtained. Both the quadratic and quartic couplings are generated by a uniform distribution in the interval [−4π, 4π]. In this work, the energy scales of all dimensional quantities such as the quadratic couplings and the field values are not relevant. For example, we can use v ≡ 246 GeV as the energy unit, instead of GeV or TeV, then µ 2 1 = 0.5 × (246 GeV) 2 = 0.5v 2 can be simply written as 0.5 in the computer program.
The quartic couplings generated in the above way can not guarantee that the potential is bounded from below. The BFB check can be numerically performed by setting the quadratic couplings to zero and then run the minimization process. If any point with V < 0 is reached during minimization, then the potential is not BFB and the sample is Table I. VEV alignments of various types of minima found in the numerical search. Except for type E, each type is defined by a relation that violates the good vacuum conditions in Eq. (27). The typical VEV alignments in each type are listed in the 3rd-5th columns, with "ci" (i = 1, 2, 3, 4) standing for independent nonzero values. Note that other VEV alignments that can be converted to the ones in this table are not shown (see the main text for more details about this issue), "· · · " indicates that the enumeration is not exhaustive. abandoned. If a sample passes the BFB check, then we further minimize the potential with the nonzero quadratic couplings. When a minimum is successfully obtained in this process, we check if it violates the four good vacuum conditions in Eq. (27), from (a) to (d) sequentially. If any of them is violated, then we stop checking the remaining conditions and tag it as type A, B, C, or D, corresponding to the violation of condition (a), (b), (c), or (d), respectively. If all the conditions are satisfied, it is tagged as type E, a good and successful vacuum.
In Fig. 3, we present the result of the above analysis on 144891 samples (all passing the BFB check). In this randomly generated data set, 70380, 70625, 683, and 3135 of the samples fall into the categories of type A, B, C, and D, respectively. Only 68 samples are of type E, which is about 0.05% of the total number. We further inspect the VEV alignments of all samples of the five types and list them in in Tab. I. As we have discussed in Sec. III B, some symmetry transformations can transform the VEV alignments from one form to another-see, e.g. Eq. (26).
For such cases, in Tab. I we only list the representative forms. More specifically, whenever appropriate U L and U R transformations are allowed (e.g. φ = 0, or κ 1 = κ 2 , or U L = U R = iσ 2 ), we always use them to transform ∆ L and ∆ R to lower triangular forms according to the Schur decomposition.
The low percentage (0.05%) of type E can be understood from the analytical studies in Sec. IV. First, for those simplified cases we have studied, one can see that the quadratic couplings µ 2 1 , µ 2 2 , and µ 2 3 have to be positive to get nonzero φ and ∆ L/R . Let us assume that for more general potentials (e.g. λ 4 , β 1, 2, 3 are no longer zero) this conclusion approximately holds as well. Then requiring the three quadratic couplings to be positive in the random number generation already produces a factor of (1/2) 3 = 1/8 which suppresses the percentage by one order of magnitude. Moreover, in Eq. (51) and Eq. (52) some quartic couplings may also need to be positive to get a good vacuum. If 11 of the 17 parameters in the full potential are required to be positive, the suppression factor can easily reach (1/2) 11 ≈ 0.05%. Some parameters may contribute suppression factors smaller or larger than 1/2, say 1/p. Generally it is possible to get a significant suppression at the order of (1/p) n where n ≤ 17. This explains why the percentage can be suppressed to the level of 0.05%.
Although the suppression is understandable, it is would be better to avoid the suppression or at least to know a part of the parameter space that would lead to the correct symmetry breaking with a much higher probability. According to our analytical study, we are led to simple conditions to enhance the probability of ending in a good vacuum: Requiring the constraints in Eq. (54), we repeat the numerical process used to generate Fig. 3 and obtain the left plot in Fig. 4. As the plot shows, with these constraints, the percentage of type E is enhanced to 16.2%, which is at the same order of magnitude as the other types. Eq. (54) is proposed based on the analytical result summarized in Eqs. (50)- (53), but it allows more general parameter settings, e.g. λ 2 , ρ 2, 4 , α 1, 2, 3 do not have to be fixed to zero. It is a compromise between the generality (also simplicity) and the enhancement of the percentage. Including more constraints from Eqs. (50)- (53) can further enhance the percentage at the cost of loss of generality. In the right panel of Fig. 4, we include the constraints together with Eq. (54) and obtain a much higher percentage (86%) of type E, which means the majority of the samples generated under these constraints have type E vacua.
As we have seen that among those randomly generated samples some may have type E vacua and some may not, if the potential is required to lead to successful symmetry breaking, it must be subjected to a lot of constraints. Below we would like to study such constraints.
Note that these constraints are not fully equivalent to the requirement that the potential has a type E minimum. If the potential has no type E minimum, then of course it can not lead to successful symmetry breaking. But even if it has a type E minimum, the minimum could be a local one which coexists with other much deeper minima. Then the vacuum at the type E minimum is not absolutely stable as it may decay to other deeper vacua via quantum tunneling or thermal fluctuation. There is a possibility that the decay rate is very low so that the lifetime is longer than the age of the Universe, known as the meta-stability. Since the analyses including meta-stability would be too much involved here, for simplicity, we only consider the absolute stability. Therefore in what follows, when we claim that a potential can lead to successful symmetry breaking, we mean the potential has a global minimum of type E. The corresponding constraints will be referred to as the global minimum constraints.
Let us investigate the effect of global vs. local minimum. Since there is no minimization algorithm that can guarantee to find global minima for general cases, the method of global minimum test used in our program is repeating the minimization process with random initial values for many times. If none of these minima is deeper than the one being tested, the more likely it is a global minimum. Obviously the more times the process is repeated, the more likely it is a global minimum. This method still can not guarantee the correctness of the global minimum, but with a large number of repetitions, the result will be very reliable. To quantify the effect of testing for a global minimum, we present now constraints on some parameters as illustration, choosing the six α 1, 2, 3 and β 1, 2, 3 parameters. In Fig. 5 the green region has successful symmetry breaking with a global minimum, while the blue region only is a global minimum. Each plot shows the constraint on a pair of the parameters (α i − α j or β i − β j ) while the other parameters  Figure 5. Global minimum constraints on α and β. Samples in the green region have type E global minima so they will lead to successful symmetry breaking; the blue region violates the global minimum constraints which means that either the potential does not have type E minimum or its type E minimum is local; the black region violates the BFB condition. Other potential parameters, if not indicated by the plots, are fixed at the benchmark values in Eq. (56). The grid interval is 0.2π. are fixed at: benchmark : This benchmark is set in such a way that by default (i.e. no parameters are changed) it has a global minimum of type E. The plots are produced in coarse grids because for each sample the program has to run the global minimum test for many times which is CPU intensive. So currently we cannot compute too many samples with limited computer power and consequently the interval of grid scan cannot be too small. In Fig. 5 we use 40 × 40 grids in the range [−4π, 4π] 2 with a interval of 0.2π. As one can see in the left panels of Fig. 5, the green regions cover the central point α 1, 2, 3 = 0 and the nearby part (within 3 or 4 blocks) is also green. This implies that small α i indeed can lead to absolutely stable type E vacua, which is a conjecture of our analytical study. This is approximately true also for the β parameters. However, the difference is that the α i do not have to be small (in some direction they can reach 4π) while the β i , at least for this benchmark, have to be small.
We also note that Fig. 3 and Fig. 4 establish our claim in the introduction that the BFB conditions only provide a necessary but not sufficient condition for an acceptable vacuum since all the columns in Fig. 3 and Fig. 4 satisfy the BFB condition whereas only the green column satisfies the desired vacuum condition.
Because the scalar mass spectrum is fully determined by the parameters of the scalar potential, the global minimum constraints on the potential parameters can be converted to constraints on the scalar mass spectrum. After symmetry breaking the scalar sector contains (including the Goldstone bosons) eight electrically neutral bosons, four singly charged bosons and two doubly charged bosons, among which the bosons with the same electric charge generally have mass mixing. Therefore the mass matrices of the neutral and singly charged scalar bosons are quite complicated, but the mass matrix of doubly charged bosons is much simpler. For simplicity, we will thus only discuss the mass spectrum of the doubly charged bosons. Their mass matrix is where Note that if ρ 4 and β 1, 2, 3 are set to zero, then the mass matrix is diagonal and one can immediately obtain the eigenvalues (i.e. the mass squares of the two doubly charged bosons): This greatly simplifies the scenario and we would like to take it as an example to show the global minimum constraints on the mass spectrum. We also set other parameters to the following specific values: except for ρ 1 , ρ 2 and ρ 3 . We take ρ 1 and ρ 3 as free parameters ranging from −ρ max to ρ max , and fix ρ 2 at some values indicated in Fig. 6. Here ρ max is set at a small value 10 −3 (for a reason to be explained below) and v 2 X is a floating energy scale which is always tuned to make κ 2 1 + κ 2 2 = 246 GeV. To make the scenario considered here more realistic, we also require that it contains a SM-like Higgs boson with the mass m h ≈ 125 GeV and a large v R so that the mass of W R is above the LHC constraints. The parameters in Eq. (59) have been tuned in such a way that for ρ 1 and ρ 3 varying in [−ρ max , ρ max ], the SM Higgs mass m h ranges within 125 ± 1.5 GeV and v R ranges from 21 to 67 TeV. We have checked that changing ρ max within one order of magnitude leads to very little change of (M ±± 1 , M ±± 2 ) or m h . Changing ρ max , however, has significant impact on v R . Generally larger ρ max leads to smaller v R , which is the reason why we use a small ρ max here.
With the above parameter setting, we scan the parameter space and compute the mass spectrum of the doubly charged Higgs bosons when the sample satisfies the global minimum requirement. The result is shown in Fig. 6, where the yellow, green, and blue regions are the allowed regions by the global minimum requirement for ρ 2 = 0.1ρ max , 0.2ρ max , and 0.5ρ max respectively. Note that the constraints presented in Fig. 6 are only for a very specific parameter setting so they should not be interpreted as universal constraints on the mass spectra. Changing the parameter setting can easily lead to significant changes of the constraints, as illustrated by different values of ρ 2 in Fig. 6. It is still interesting that there are certain mass ranges for the doubly charged scalars which are forbidden by our analysis and can be used to test the model.

VI. CONCLUSION
We have performed a study on the vacuum structure of the left-right symmetric scalar potential. The goal was to investigate whether the usually considered VEV alignment in Eq. (10) can be obtained from a generic scalar potential as a global minimum. General criteria to identify a charge conserving and parity violating vacuum were derived (see Eq. (27), and it was demonstrated that the potential parameters are subject to many constraints in order to achieve this minimum. In general if we do not put any constraints on the potential parameters, as indicated by Fig. 3, the probability to end up within the desired VEV alignment is very low, only 0.05%.
We have also analytically studied the minima of the potential in the absence of some terms and obtained conditions that enable us to obtain the correct VEV alignment more easily, as shown in Fig. 4. By requiring that the corresponding minimum is global in this case, we also illustrate the constraints on the potential parameters in Fig. 5.
Our work suggests that successful generation of the usually considered VEV alignment in standard left-right symmetric theories and keeping the vacuum absolutely stable would produce important constraints on the parameters of the potential. These constraints may have interesting phenomenological consequences such as constraints on the mass spectrum of scalar bosons, or the Higgs self-couplings, etc., which can be used to test the model. We have given the example of doubly charged scalar masses in the model as an example.
The present paper can be a starting point for further and much more involved analyses along these lines, such as analyzing loop-corrected effective potentials, investigating the vacuum lifetime of non-global minima, more phenomenological consequences of the global minima, or studies of alternative left-right symmetric models. 1 +2κ 2 cos θ 2 3κ 2 1 λ 4 + κ 2 2 λ 4 − 2µ 2 2 + α 2 v 2 L + α 2 v 2 R (A1)