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 $10^{10\text{--}12}$ 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 $10^{16\text{--}17}$ 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] 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 [3]. This feature is a great advantage compared to 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 [4]. 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][11][12] GeV [5,6]; see also [7,8]. This model also gets renewed attention in the context of the Twin Higgs model as a solution to the hierarchy problem as well as the parity solution to the strong CP problem [9,10]. 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. [5,6] 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 considerably 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][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 Sec. 2, we summarize the SO(10) model which has the minimal LR-symmetric gauge group at the intermediate stage. In Sec. 3, we discuss the gauge coupling unification in the minimal LR symmetric model. In Sec. 4, we study the proton lifetime. In Sec. 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. 2 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 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 4 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 represents 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 . 5 . 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 Sec. 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 Here and hereafter, the overline on the extra fields denotes a charge conjugation rather than a Dirac adjoint. By using the G LR representations, the Eq. (6) is decomposed as 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 follow (see also [9,10]).
• 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 [12].
• 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 [3] tend to be much heavier than the observed ones. This is because the Dirac neutrino Yukawa coupling for 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. [13]). 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 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. Above M R , those coefficients are given by Below M R , on the other hand, they are given by which come only from the SM particle contribution. 6 To calculate the RG flow for the gauge couplings, we consider the one-loop matching condition at the renormalization scale, Here, α 2R is the gauge coupling of the SU (2) R gauge interaction, and in this model, the value of gauge coupling for SU (2) R group is the same as that for SU (2) L group above M R , i.e. α 2R = α 2L = α 2 . As we are taking the MS renormalization scheme, there is a mass independent threshold correction in the right-hand side [14]. 7 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. 8 In Figure 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 [15]: The gauge couplings for M R = 10 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 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 As a measure of the quality of the unification, we define, where we take µ = Λ = M X . The definition of∆ is different from the unification measure ∆ defined in [9,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 , while ∆ in [9,10] 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.∆ 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 . 9 The parameters ∆1,2,3 also get contributions from higher dimensional operators H45 /M Pl , although we assume that H45 /M Pl O(1).  Figure 2: 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.

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 anti-symmetric representation 10 of the minimal SU (5) GUT gauge symmetry, respectively. In general setup of the SO(10) GUT, they have different masses (see e.g. [17]), 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 4. Those operators are reduced to in terms of the G SM fields [18] (see also [19]). 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 4. In Eq. (20), we do not take account of the effects of the quark mixing angles [20]. 10 The partial decay widths 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 as I=1, 3 In this calculation, W 0 for p → π 0 e + decay mode is −0.131 GeV 2 which have been obtained by lattice simulation [21].
To calculate the coefficients of the proton decay operators at the proton mass scale m p , we have to consider the renormalization factor A. In this paper, we consider the one-loop 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 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 a-th gauge interaction for an each proton decay operator: 10 In the present model, an SM fermion is a linear combination of the spinor F16 and the extra particles E10 and E45 [9,10], 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 contribution from extra particles in this paper. The proton decay operators which come from the gauge interaction of the extra particle E10 are summarized in Ref. [16].  Figure 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 grey 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.
The coefficient C a is summarized in Ref. [22], with which the renormalization factors are given by 11 In Eq. (26), we double 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, 12 In Figure 3, we overlay the current limit and the future prospects on the proton lifetime for p → π 0 e + decay mode on Figure 2. The current limit is the 90%CL exclusion limit by Superkamiokande (SK) experiment, 1.6 × 10 34 years [23], which is shown as the black solid line. The future prospects is the expected exclusion limit at 90%C.L. of the Hyper-K (HK) experiment, 1.3 × 10 35 years [24], which is shown as the black dashed line. The figure shows that some of the parameter region with moderate coupling unification have 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 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.

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 solution to the Strong CP problem [25,26]. 13 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 [27,28]. The axion model not only solve the strong CP problem, but also provides a good candidate for cold dark matter [29,30,31,32]. In fact, the axion dark matter model is successful when the PQ breaking scale is of 10 11-12 GeV, which is close to the LR-breaking scale discussed in this paper; see [33] 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 that the PQ charges of F 16 are 1/2, while that of H 16 is vanishing. 14 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 13 Alternatively, the strong CP problem can be solved in the LR symmetric model by imposing space-time parity appropriately; see [9,10] and references therein. 14 We may consider a model in which E10 and E45 have the opposite PQ charges. In this case the PQ charge of F16 is vanishing, although the domain wall number is again NDW = −21.
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. 15 As both the extra fermions as well as the SM fermions possesses the PQ charges, this model is in between the KSVZ [34,35] and DFSZ [36,37] invisible axion models, and is in principle distinguishable from these models.
The coherent oscillation of the axion turns into the dark matter density [38], where we have defined F eff = f a /N DW . ∆a i /F eff ∈ [−π, π) 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 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 higher dimensional operator to generate the SM Yukawa interactions of the first generation explicitly break the PQ symmetry. 16 We argue that the axion in our setup is within the reach of future detection. Due to the non-vanishing axion potential, the axion get a mass given by [39] m a 5.7 µeV 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 axionphoton coupling, which is parameterized as with [40] 17 Note that g aγγ in our model is equivalent to that in the DFSZ axion model [36,37], which is already excluded by the current ADMX experiment for m a 2.7 -3.3 µeV [41,42]. 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 [41], CULTASK [43] and MADMAX [44]; see also Ref. [45]. 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 PQsymmetry 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 [46] which requires a rather high reheating temperature, T R 10 9-10 GeV [47,48,49]. 18 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 [53]. 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. Ref. [33,54]. 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. Ref. [55,57].
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 [56,57,58,59,60,61]. When explicit breaking terms exist, the effective θ angle of QCD is non-vanishing 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 [62]. 19 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. [63]. This argument strengthens the PQ mechanism in the present model. 20

Summary
In this paper, we have investigated the proton lifetime in the SO(10) GUT which is broken down by the VEV of H 45 to the minimal LR-symmetric gauge group SU  10 17 GeV from that expected in Refs. [5,6] 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-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. axion domains in a/F eff = N DW × [−π, π) are gauge equivalent with each other, and hence, the axion domain wall configurations which connects different domains are not completely stable. As we will see, however, the axion domain wall problem remains even in the model with the discrete gauge symmetry.
To make our discussion concrete, let us assume that the discrete Z N DW gauge symmetry originates from a U (1) gauge symmetry broken by the VEV of a complex scalar Φ whose gauge charge is large, N DW 1. Note that this U (1) gauge symmetry is different from the global U (1) PQ symmetry. The U (1) gauge charge of the PQ breaking field P is 1 as in Sec. 5. 21 The VEV of the PQ breaking field P eventually breaks the Z N DW symmetry.
In this model, the stable topological defect is not the domain wall but the local strings which are associated with the spontaneous U (1) gauge symmetry breaking. For example, a cosmic local string around which the phase of Φ winds from 0-2π are expected to be formed when Φ obtains a VEV at a very high energy scale. The phase of the PQ breaking field P is changed by 2π/N DW under the parallel transport around this local string, which corresponds to the Aharanov-Bohem effect.
Now let us assume that the spontaneous symmetry breaking of the U (1) gauge symmetry by the VEV of Φ takes place well before inflation, while the approximate global PQ symmetry breaking occurs after inflation. In this case, the cosmic local strings that are formed when Φ obtains a VEV have been diluted away by cosmic inflation. After inflation, the cosmic temperature decreases below the PQ breaking scale. Then, associated with the spontaneous breaking of the approximate global U (1) PQ symmetry, a few cosmic global strings are expected to be formed in each Hubble volume. Note that these global strings are different from the ones diluted away during the inflation. When we turn around the global string, the phase of the PQ field P takes values from 0 to 2π when the winding number is one, and hence, the axion field takes values from 0 to f a × 2π = N DW × F eff × 2π.
Around the global string, the [0, N DW × F eff × 2π) region has N DW domains that are gauge equivalent under the Z N DW . Since the approximate U (1) PQ symmetry is highly protected by the Z N DW symmetry, the tension of the domain walls connecting the N DW domains is negligibly small. Therefore, we have no domain wall problem associated with the Z N DW symmetry breaking by P = 0. When the cosmic temperature decreases further, the cosmic global string networks follow the so-called scaling solution where the number of the cosmic global strings in each Hubble volume at that time remains of O(1); see e.g. [68]. When the axion potential is generated at around the QCD scale Λ QCD , potential barriers appear around each global string which result in N DW domain walls whose boundary is the global string.
As mentioned earlier, each domain wall attached to the global string connects different domains which are gauge equivalent under the discrete Z N DW symmetry. Therefore, this domain wall is not completely stable. In fact, each wall can be punctured by a loop of the earlier mentioned local string, around which the phase of Φ winds from 0-2π, since this local string connects the different axion domains without 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 . 22 As a result, the domain wall is virtually stable below the QCD scale and they immediately dominate over the energy density of the universe, which causes the domain wall problem. 23