Proton decay and axion dark matter in SO(10) grand unification via minimal left–right symmetry

We study the proton lifetime in the SO(10) Grand Unified Theory (GUT), which has the left–right (LR) symmetric gauge theory below the GUT scale. In particular, we focus on the minimal model without the bi-doublet Higgs field in the LR symmetric model, which predicts the LR-breaking scale at around 1010-12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{10-12}$$\end{document} GeV. The Wilson coefficients of the proton decay operators turn out to be considerably larger than those in the minimal SU(5) GUT model especially when the Standard Model Yukawa interactions are generated by integrating out extra vector-like multiplets. As a result, we find that the proton lifetime can be within the reach of the Hyper-Kamiokande experiment even when the GUT gauge-boson mass is in the 1016\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{16}$$\end{document}–1017\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$10^{17}$$\end{document} GeV range. We also show that the mass of the extra vector-like multiplets can be generated by the Peccei–Quinn symmetry breaking in a consistent way with the axion dark matter scenario.


Introduction
The Grand Unified Theory (GUT) [1] is one of the most attractive candidates for physics beyond the Standard Model (SM), which provides an explanation of the charge quantization. In particular, the SO(10) gauge group [2,3] is one of the most attractive candidates for the unification group as it not only unifies all the gauge interactions in the SM but also unifies a generation of the SM fermions into one representation. Furthermore, it also predicts the existence of the right-handed neutrinos, which naturally explains the light active neutrino masses through the seesaw mechanism [4][5][6][7][8]. This feature is a great advantage over the SU (5) GUT.
Another interesting feature of the SO (10) GUT is that the rank of SO (10) is larger than the SM. Accordingly, the SO(10) GUT allows various symmetry breaking paths to the SM gauge groups, such as the Left-Right (LR) symmetric groups [9][10][11][12][13][14][15]. Among these possibilities, the minimal model based on the SU (3) C × SU (2) L × SU (2) R ×U (1) B−L gauge group without a bi-doublet of SU (2) L × SU (2) R uniquely predicts an intermediate breaking scale of the LR symmetry to be around 10 10−12 GeV [16,17]; see also [18,19]. This model also gets renewed attention as it can explain the small Higgs quartic coupling constant at a high energy scale while solving the strong CP problem simultaneously [10][11][12][13][14][15]20,21]. In this class of models, all the SM Yukawa interactions are generated by integrating out extra vector-like multiplets at around the LR-breaking scale.
In this paper, we discuss the proton lifetime in this scenario with the simplest possibility of the extra matter content. 1 As we will see, the preferred GUT scale 10 17 GeV is lower than expected in Refs. [16,17] by a factor a few or so, due to the effects of the extra matter multiplets on the renormalization group running. 2 We also find that the Wilson coefficients of the proton decay operators are consid- 1 Hall and Harigaya have extensively studied various possibilities of the extra matter contents [20,21]. In this paper we only study the simplest one among them, which has also been partly discussed in a different context of the axion dark matter scenario [22]. 2 Recall that the proton decay rate is inversely proportional to the fourth power of the GUT gauge-boson mass M X , and hence the slight change of M X significantly affects . erably larger than those in the minimal SU (5) GUT model due to the larger gauge coupling below the GUT scale as well as the SU (2) R gauge interaction at the intermediate scale. As a result, the proton decay rate is enhanced and a parameter region consistent with the gauge coupling unification in the 10 16 -10 17 GeV range can be tested by the Hyper-Kamiokande (Hyper-K) experiment. We also discuss a possibility to generate the mass of the extra vector-like multiplet by the Peccei-Quinn (PQ) symmetry breaking in a consistent way with the axion dark matter scenario.
The organization of the paper is as follows. In Sect. 2, we summarize the SO(10) model which has the minimal LR symmetric gauge group at the intermediate stage. In Sect. 3, we discuss the gauge coupling unification in the minimal LR symmetric model. In Sect. 4, we study the proton lifetime. In Sect. 5, we discuss the mass generation of the extra vector-like multiplets by the PQ symmetry breaking. We give a summary of our discussion in the final section.

The minimal setup of the SO(10) GUT model
In this paper, we discuss SO(10) GUT with the following chain of symmetry breaking: To ensure this chain and subsequent SM symmetry breaking, we introduce an SO (10) where T R 3 is the third generator of SU (2) R . As will be discussed, the typical values of the GUT and the LR symmetry breaking scales are M GUT = O(10 16−17 ) GeV and M R = O(10 10−12 ) GeV, respectively. Throughout this paper, we assume these minimal contents for the Higgs sector, and assume that only the doublet Higgs bosons of SU (2) R and SU (2) L remain massless below the GUT scale.
In the minimal SO(10) GUT model, each generation of the quarks and the leptons of the SM forms an SO(10)-spinor F 16 , which is decomposed into the G LR and G SM representations as where the subscript is for the charges of B-L and Y , respectively. To embed U (1) B−L into SO(10), we renormalize the charges so that the U (1) gauge couplings are given by 5 In the LR symmetric model with only SU (2) L ,R doublet Higgs bosons, the Yukawa interactions in the SM are given by the higher-dimensional operators in Table 1. In the SO(10) notation, they correspond to where i, j = 1, 2, 3 is the flavor index. is the cutoff scale. Hereafter, we suppress the gauge and the flavor indices unless otherwise stated. After the LR symmetry breaking, these operators contribute to the Yukawa interactions: y u contributes to the up-type and neutrino ones, while y d to the down-type and charged-lepton ones. In Table 1, the second and the third columns represent the Yukawa interactions from the higher-dimensional operators in Eq. (5) in the representations of G LR and G SM , respectively.
Obviously, these contributions are too small to realize the observed masses of the heavy flavor fermions in the SM for = M GUT , for example. In fact, since the LR symmetry breaking scale M R is around 10 10 -10 11 GeV, while M GUT = 10 16 -10 17 GeV, the coefficient of these operators are ∼ M R / = 10 −7 -10 −5 , and hence we cannot realize the Yukawa couplings for the second and third generations. To reproduce the observed quark and lepton masses, we need to introduce extra vector-like multiplets with masses of M R so that the terms in Eq. (5) are generated by integrating out those extra multiplets.
In this paper, we assume that all the SM Yukawa interactions are generated by integrating out extra vector-like multiplets. In this case, the minimal extra vector-like fermions consist of three flavors of the fundamental representation of SO(10), E 10 , and three flavors of the adjoint representation of SO(10), E 45 . 6 When the Yukawa interactions of the first generation are provided by the M GUT suppressed operators, two flavors of E 10 and E 45 are enough to reproduce the SM Yukawa interactions. As discussed in Sect. 5, however, the three flavor model is advantageous as the masses of the extra vector-like fermions can be interrelated to the PQ symmetry breaking. In what follows, we denote the number of extra particle flavors by N E .
With these extra matter multiplets, the origin of the Yukawa interactions in Eq. (5) are obtained from the renormalizable interactions, where M extra is the extra particle mass. We assume that the mass parameters for E 10 and E 45 are the same for simplicity. E 10 and E 45 are decomposed into the G LR representations as 6 See [20,21] for other possibilities.
Here and hereafter, the overline on the extra fields denotes a charge conjugation rather than a Dirac adjoint. By using the G LR representations, Eq. (6) is decomposed as (45) .
When we integrate out the extra particles, these contributions become the higher-dimensional operators which are summarized in Table 1. The resultant Yukawa coupling constants in the SM are proportional to M R /M extra , and hence, the top Yukawa coupling requires the extra particle masses should be around LR symmetry breaking scale. Several comments of the above minimal setup is summarized as follows; see also [20,21].
• The large difference between the top mass and the other third generation one is the most serious problem in realization of the observed fermion masses in generic SO (10) GUT models. This is not a matter in our model because there are two origins of the Yukawa interactions. • Small difference between down-type quark masses and charged-lepton masses is introduced by higher-dimensional operators that come from SO(10) breaking effects [25]. • The right-handed neutrino masses are around LR-breaking scale M R . As we assume that the LR-breaking scale is around 10 10 -10 12 GeV, the masses of the active neutrinos generated by the seesaw mechanism [4][5][6][7][8] tend to be much heavier than the observed ones. This is because the Dirac neutrino Yukawa coupling for the third generation is O(1) since it is unified with the top Yukawa coupling. • When we assume that large mixings in the MNS matrix are realized, the CKM matrix also should be a large mixing matrix because of the unification. However, this does not satisfy experimental results.
We can solve the above problems by cancellation between the contributions from the operators in Eq. (5) with some other higher-dimensional operators which include the GUT breaking effects. The latter operators are suppressed by a factor of H 45 / . However, in this model, the suppression factor is not so small even if is around the Planck scale, M Pl ∼ 2.4 × 10 18 GeV, as H 45 can be as large as around 10 17 GeV. By the cancellation, the small neutrino Yukawa coupling can be achieved even for the O(1) top Yukawa coupling, and hence the active neutrino masses satisfy the experimental results. The mixing matrices of the quarks and the neutrinos can also be consistent with each other by cancellation.

Gauge coupling unification
In the previous section, we introduce extra fermions to achieve the Yukawa interactions of the SM. In this section, we consider the renormalization group (RG) flow of the gauge couplings including the contributions of those extra matter multiplets as well as the SU (2) R doublet Higgs boson. We assume for simplicity that the masses of the extra fermions and SU (2) R doublet Higgs are M R . As the extra fermions makes the gauge coupling constants become rather strong at around the GUT scale, it is important to take into account the two-loop contributions of the gauge coupling constants to the RG flow; see e.g. Refs.
[26 -28]. The extra Yukawa interactions to the two-loop RGE may slightly affect the precision of the unification and the GUT gauge-boson mass. Since those effects depend on the detailed mass spectrum of the extra fermions, we neglect those contributions in this paper. The β function of the gauge coupling g a is given by where a, b take values 1, 2, 3 which refer to U (1) Y , SU (2) L , and SU (3) C below M R and take values 1, 2L, 2R, and 3 which refer to U (1) B−L , SU (2) L , SU (2) R , and SU (3) C above M R , respectively: • Above M R , the coefficients of the gauge coupling beta functions are where each of a 0 and b 0 contains contributions from the SM particles and the SU (2) R doublet Higgs; a 10 and b 10 from E 10 ; a 45 and b 45 from E 45 ; and N E is the number of extra particle pairs. These coefficients above M R are given by 7 We set the SU (2) L and SU (2) R gauge couplings equal, which come only from the SM particle contribution.
To calculate the RG flow for the gauge couplings, we consider the one-loop matching condition at the renormalization scale, Recall that the value of gauge coupling for SU (2) R group is the same as that for the SU (2) L group above M R : α 2R = As we are taking the MS renormalization scheme, there is a mass independent threshold correction in the righthand side [29]. 8 In the following, we assume that the massive gauge boson of SU (2) R × U (1) B−L and the extra matter multiplets E 10, 45 have the same mass of M R for simplicity. The contributions of the extra matter do not affect the quality of the unification significantly as long as they have SO(10) consistent masses. 9 In Fig. 1, the RG flow of the gauge couplings is shown. The input values for the RG flow are taken to be the central values of the experimental measurements in [30]:  11 GeV, on the other hand, meet well together before they hit the Landau pole. There, we see that the two-loop contributions are not negligible with which the 8 For the DR renormalization scheme, the mass independent threshold correction is absent. 9 The mass scales of the extra vector-like matter affect the size of the GUT gauge coupling constant. The mass splitting within the extra SO(10) multiplets also affect the precision of the gauge coupling unification at the GUT scale.
RG flow becomes non-linear. The results for M R = 10 12 GeV also show that the gauge couplings become close with each other moderately at around M R = 10 15 GeV.
To quantify the quality of the unification, let us consider the matching conditions between the gauge coupling constants in the LR symmetric model and the SO(10) gauge coupling, α G = g 2 G /4π : The  The quality of the unification¯ as a function of (M X , α −1 G ). The upper and the lower panels are for N E = 3 and N E = 2, respectively. The quality of the unification is reasonably high in the blue shaded region (¯ < 5), while it is moderate in the light-blue shaded region (¯ < 10). The parameter¯ gets contribution not only from the mass splittings of the GUT multiplets but also from the mass difference between the GUT particles and M X . The pink shade region is excluded as M X is above the Landau pole of α 1,2,3 (μ, M R ). We confine ourselves to the region with M X M Pl , so that the effective field theory without gravity is valid other than the GUT gauge bosons, although we do not specify them in this paper. 10 See [20,21] for various contributions of the GUT multiplets to 1,2,3 .
As a measure of the quality of the unification, we definē where we take μ = = M X . The definition of¯ is different from the unification measure defined in [20,21]. The parameter¯ gets a contribution not only from the mass splittings of the GUT multiplets but also from the mass difference between the GUT particles and M X , while in [20,21] purely measures the precision of the unification.
In Fig. 2, we show¯ as a function of (M X , α −1 G ) for a given M R . The quality of the unification is reasonably high in the blue shaded region (¯ < 5), while it is moderate in the light-blue shaded region (¯ < 10). The figure shows that a reasonable unification, i.e.¯ < 5, is not possible for  For comparison, we also show¯ for N E = 2. In this case, the Landau pole is at the very high energy scale and does not exclude the parameter region significantly. For N E = 2, more precise unification is achieved for a lower M R and a higher M X than the case of N E = 3. In such a parameter region, however, there is a tension with the possibility to obtain the first generation Yukawa couplings as the higher-dimensional operators suppressed by M GUT .

Proton lifetime
In the present model, the exchanges of the massive gauge boson in the (3, 2, 2) −2/3 representation, i.e. the X -type gauge bosons, induce the proton decay. Incidentally, the each of the SU (2) R doublet component of the X gauge boson belongs to the adjoint representation 24 and the antisymmetric representation 10 of the minimal SU (5) GUT Table 2 The B and L violating operators mediated by the X gauge boson gauge symmetry, respectively. In general setup of the SO(10) GUT, they have different masses (see e.g. [31]), while they are common in the LR symmetric model. The massive gauge boson in the (3, 1, 1) 4/3 representation, on the other hand, does not lead to the proton decay. After integrating out the X gauge boson, the gauge interaction of the matter field F 16 results in the B and L breaking operators O (1,2) in Table 2. Those operators are reduced to (20) in terms of the G SM fields [32]; see also [33]. Below the electroweak symmetry breaking scale, we may decompose it into the proton decay operators in terms of the SU (3) C × U (1) em fields such that L eff = C I O I as in Table 2. In Eq. (20), we do not take account of the effects of the quark mixing angles [34]. 11 The partial decay width for the p → π 0 e + is given by where m p and m π 0 are the proton and the neutral pion masses, respectively, and W I 0 are the proton form factor. We may safely approximate I =1, 3 In this calculation, W 0 for the p → π 0 e + decay mode is −0.131 GeV 2 , which has been obtained by a lattice simulation [37].
To calculate the coefficients of the proton decay operators at the proton mass scale m p , we have to consider the 11 In the present model, an SM fermion is a linear combination of the spinor F 16 and the extra particles E 10 and E 45 [20,21], and therefore we should consider the proton decay operators which come from the gauge interactions of the extra particles too, strictly speaking. However, we have introduced the extra particles to realize the large Yukawa couplings, while the Yukawa couplings of the first two generations are small. Therefore, we expect that the contributions from extra particles are small for the first generation, and thus we do not consider a contribution from extra particles in this paper. The proton decay operators which come from the gauge interaction of the extra particle E 10 are summarized in Refs. [35,36]. renormalization factor A. In this paper, we consider the oneloop level renormalization factor from gauge interactions. Here, we divide the energy region into two parts. The first region is between the GeV scale and the LR-breaking scale M R , where the renormalization factor is written as A long . The second region is between the LR-breaking scale M R and the GUT scale M GUT , where the renormalization factor is written as A short . The total renormalization factor A is given by the product of these factors, A = A long × A short . We calculate this renormalization factor for each of the proton decay operators O (1) and O (2) .
The one-loop level renormalization factor for each gauge group is given by where M end > M start ; a a is the coefficient for β function for each gauge coupling which are shown in Eqs. (11) and (14). C a is the factor appearing in the anomalous dimension γ a of the ath gauge interaction for an each proton decay operator: The coefficient C a is summarized in Ref. [38], with which the renormalization factors are given by 12 (26) 12 We use the six flavor RG equations even below the electroweak scale. If we instead use the three flavor RG below the electroweak scale, the A (1,2) long are slightly enhanced by about 10%.

Fig. 3
The black solid and black dashed lines are proton decay constraints on the p → π 0 e + decay mode from current SK limit and the future HK prospect. The gray shaded region is excluded by the current SK limit and the region between black solid line and black dashed line will be explored by the HK experiment In Eq. (26), we double the SU (2) L contribution to include the contribution from the SU (2) R gauge interaction. For M R 10 11 GeV, M X 10 16.5 GeV and N E = 3, for example, we find that the renormalization factors are given by 13 In Fig. 3, we overlay the current limit and the future prospects on the proton lifetime for p → π 0 e + decay mode on Fig. 2. The current limit is the 90% CL exclusion limit by Super-kamiokande (SK) experiment, 1.6 × 10 34 years [39], which is shown as the black solid line. The future prospects is the expected exclusion limit at 90% CL of the Hyper-K (HK) experiment, 1.3×10 35 years [40], which is shown as the black dashed line. The figure shows that some part of the parameter region with moderate coupling unification has been excluded by the current SK limit for M R 10 11.5 GeV (N E = 3). The figure also shows that the HK experiment has a sufficient sensitivity to test large portion of the parameter space with moderate coupling unification for M R = O(10 11 ) GeV for N E = 2, 3. 13 For this choice, we find α −1 1 (M X , M R ) 6.9, α −1 2 (M X , M R ) 7.3, and α −1 3 (M X , M R ) 6.0, respectively.

Model with Peccei-Quinn symmetry
In the minimal setup with N E = 3, we assume that all the SM Yukawa interactions are generated by integrating out the extra vector-like multiplets with masses around the LR-breaking scale. In this section, we briefly discuss a possibility to generate those masses by the PQ symmetry breaking. The PQ mechanism is one of the most successful solutions to the strong CP problem [41,42]. 14 There, the effective θ -angle of QCD is canceled by the VEV of the pseudo-Nambu-Goldstone boson, axion a, which is associated with the spontaneous breaking of the PQ symmetry [43,44]. The axion model not only solves the strong CP problem, but also provides a good candidate for cold dark matter [45][46][47][48]; see also Refs. [49][50][51][52]. In fact, the axion dark matter model is successful when the PQ breaking scale is of 10 11 -10 12 GeV, which is close to the LR-breaking scale discussed in this paper; see [53] for review. This coincidence motivates us to see how it is successful to the mass scale of the extra vector-like fermions with the PQ breaking scale. For this purpose, let us introduce a gauge singlet complex scalar field, P, which breaks the PQ symmetry at an intermediate scale. The PQ charge of P is defined to be 1. Below the PQ breaking scale, the axion appears as a phase component of P, where f a is the decay constant of the axion. The PQ symmetry is realized by the shift of a, where the domain of the axion is given by a/ f a = [−π, π).
To generate the extra fermion masses at the PQ scale, we assume that P couples to E 10,45 via, where k 10,45 are the coupling constants. Here, we assume that the PQ charges of E 10 and E 45 are −1/2. In this case, the interaction terms in Eq. (6) impose the requirement that the PQ charges of F 16 are 1/2, while that of H 16 is vanishing. 15 With these charge assignments, we find that the anomalous axion coupling to QCD is given by Here, G andG are the QCD field strength and its hodge dual, respectively. N F 16 = 3 is the number of generation of the SM fermions, and N E 10 = N E 45 = N E = 3. The Lorentz and color indices are suppressed. Below the QCD scale the anomalous coupling of the axion to QCD in Eq. (31) leads to a non-vanishing axion potential and the axion settles down to its minimum which solves the strong CP problem. 16 As both the extra fermions and the SM fermions possess the PQ charges, this model is in between the KSVZ [54,55] and DFSZ [56][57][58] invisible axion models, and is in principle distinguishable from these models. The coherent oscillation of the axion turns into the dark matter density [59], a h 2 0.18 where we have defined denotes the initial misalignment angle of the axion from the N DW degenerate CP conserving vacua. Therefore, the axion dark matter scenario is successful for F eff ∼ 10 11−12 GeV for a typical initial misalignment angle. In this present model, the PQ breaking scale is given by f a = N DW F eff , the axion dark 15 We may consider a model in which E 10 and E 45 have the opposite PQ charges. In this case the PQ charge of F 16 is vanishing, although the domain wall number is again N DW = −21. 16 Here, the origin of the axion field space is taken to be the one at which the effective θ-angle of QCD is vanishing without loss of generality.
matter prefers the PQ breaking scale at f a ∼ 10 12−13 GeV. Accordingly, we find that the extra multiplet masses at M R ∼ 10 11 GeV can be provided for k 10,45 ∼ 10 −(1−2) consistently with the axion dark matter scenario. It should be emphasized that this scenario does not work for N E = 2 since the higherdimensional operator to generate the SM Yukawa interactions of the first generation explicitly break the PQ symmetry. 17 We argue that the axion in our setup is within the reach of future detection. Due to the non-vanishing axion potential, the axion gets a mass given by [60] m a 5.7 µeV 10 12 GeV F eff .
The axion also couples to photons through the electromagnetic anomaly N QED and thorough the mixing with neutral mesons. Many on-going and future axion search experiments utilize the axion-photon coupling, which is parameterized as with [61] 18 Note that g aγ γ in our model is equivalent to that in the DFSZ axion model [56][57][58], which is already excluded by the current ADMX experiment for m a 2.7-3.3 µeV [62,63]. The higher mass range of m a up to 400 μeV (corresponding to F eff ∼ 10 11 GeV) is expected to be covered by future cavity haloscopes such as ADMX [62], CULTASK [64] and MAD-MAX [65]; see also Ref. [66,67]. Several comments are in order. The axion potential induced by the anomalous QCD coupling in Eq. (31) possesses Z N DW discrete symmetry in the domain of the axion a/ f a ∈ [−π, π), or equivalently in a/F eff ∈ N DW ×[−π, π). The discrete symmetry is spontaneously broken by the VEV of the axion. Thus, the domain wall formation takes place after the onset of the coherent oscillation of the axion, if the initial misalignment angle in each Hubble volume of the Universe at that time is random. Once the domain walls are formed, they immediately dominate the Universe, which conflicts with the Standard Cosmology. To avoid this problem, we need to assume that the PQ symmetry breaking takes place before inflation and never gets restored after inflation. Under this assumption, the initial misalignment angle of the axion is uniform in the entire Universe, and hence the axion sits in the same sub-domain and evades the formation of the domain wall.
We mention that the large domain wall number, N DW = −21, is advantageous to avoid the PQ symmetry restoration, since the actual PQ breaking scale is an order of magnitude larger than the effective decay constant F eff appropriate for the axion dark matter scenario, i.e. F eff ∼ 10 11−12 GeV. Therefore, the present model can be consistent with a cosmological scenario with higher reheating temperature than in the conventional axion dark matter models. In this sense, the present model can be more easily consistent with the thermal leptogenesis scenario [68] which requires a rather high reheating temperature, T R 10 9-10 GeV [69][70][71]. 19 As another comment, the massless axion fluctuates quantum mechanically during inflation, which leads to the isocurvature fluctuation of the axion dark matter density when the PQ symmetry breaking takes place before inflation. The dark matter isocurvature fluctuation have been severely constrained by the precise measurements of the cosmic microwave background [75]. The amplitude of the isocurvature fluctuation is proportional to the Hubble parameter during inflation, H I . As a result, H I is constrained from above as H I 10 7−8 GeV to avoid the current constraint; see e.g. Refs. [53,76]. Therefore, the present scenario with the axion dark matter can be refuted if the primordial B-mode polarization in the cosmic microwave background is discovered in near future; see, e.g. Refs. [77][78][79].
Finally, let us comment on the origin of the PQ symmetry. By definition, the U (1) PQ symmetry cannot be an exact symmetry as it is explicitly broken by the QCD anomaly. Besides, it is also argued that any global symmetries are broken by quantum gravity effects [78][79][80][81][82][83][84]. When explicit breaking terms exist, the effective θ angle of QCD is nonvanishing even in the presence of the axion, which spoils the PQ mechanism. For example, if the PQ symmetry is completely broken by the quantum gravity effects, it is expected that there should be a PQ breaking term at least, which drastically affects the axion potential and spoils the PQ mechanism. In the present model, however, we may regard that the discrete Z 2N DW symmetry to be a discrete gauge symmetry as it can satisfy the anomaly free conditions [85]. 20 If Z 2N DW symmetry is a gauge symmetry, the lowest-dimensional operator which breaks the U (1) PQ symmetry but is invariant under the Z 2N DW gauge symmetry is given by which is highly suppressed and does not spoil the PQ mechanism; see e.g. Ref. [86]. This argument strengthens the PQ mechanism in the present model. 21

Summary
In this paper, we have investigated the proton lifetime in the SO (10)  Due to the absence of the bi-doublet Higgs boson, the LRbreaking scale is determined to be at around 10 10−12 GeV in order to achieve the gauge coupling unification. As a notable feature of the model, it requires extra vectorlike fermions to generate the SM Yukawa interactions. Such extra multiplets affect the RG flow, and lower the unification scale down to M X 10 17 GeV from that expected in Refs. [16,17] by a factor a few or so. We have also found that the Wilson coefficients of the proton decay operators are considerably larger than those in the minimal SU (5) GUT model. As a result, the proton decay rate is enhanced and we find that some portion of the parameter space consistent with the gauge coupling unification can be tested by the Hyper-K experiment thorough the proton decay search even when the GUT gauge-boson mass is in the range 10 16 -10 17 GeV.
We also discussed a possibility to generate the mass of the extra vector-like multiplets by the PQ symmetry breaking. We found that the axion dark matter scenario and the present model can be successfully combined for the model with N E = 3. This combination can be tested by the proton decay search, the axion search and the search for the primordial B-mode fluctuation in the cosmic microwave background. 20 Here, we normalize the charges of the discrete symmetry so that the extra vector-like multiplets have a charge −1. Accordingly, the PQ breaking field P possesses the discrete charge 2. 21 We discus the domain wall problem in Appendix A. As another possible justification of the U (1) PQ symmetry, we may consider a U (1) gauge symmetry with an accidental global U (1) PQ symmetry [87], where E 10 and E 45 couple to a different PQ charged complex scalars; see Refs. [88][89][90] for details. potential barrier. Once the domain wall is punctured, the loop of local string expands on the domain wall, and the domain wall disappears eventually. The rate of such a puncturing process, however, is highly suppressed, since the formation of the loop of the local string is suppressed by e −| | 4 / 2 QCD F eff T at a temperature below the QCD scale: T QCD . 23 As a result, the domain walls are virtually stable below the QCD scale and they immediately dominate over the energy density of the Universe, which causes the domain wall problem. 24