Welcome to visit Communications in Theoretical Physics,
Particle Physics and Quantum Field Theory

Four-family ${ \mathcal N }=1$ supersymmetric Pati–Salam models from intersecting D6-branes

  • Tianjun Li 1, 2 ,
  • Rui Sun , 3 ,
  • Chi Zhang 1, 2
Expand
  • 1CAS Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China
  • 2School of Physical Sciences, University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing 100049, China
  • 3 Korea Institute for Advanced Study, 85 Hoegiro, Dongdaemun-Gu, Seoul 02455, Republic of Korea

Received date: 2022-03-21

  Revised date: 2022-04-13

  Accepted date: 2022-04-14

  Online published: 2022-06-13

Copyright

© 2022 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing

Abstract

We investigate the construction of four-family ${ \mathcal N }=1$ supersymmetric Pati–Salam models from Type IIA ${{\mathbb{T}}}^{6}/{{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$ orientifold with intersecting D6-branes. Utilizing the deterministic algorithm introduced in He W, Li T and Sun R (2021 arXiv: hep-th/2112.09632), we obtain 274 types of models with three rectangular tori and distinct gauge coupling relations at string scale, while 6 types of models with two rectangular tori and one tilted torus. In both cases, there exists a class of models with gauge coupling unification at string scale. In particular, for the models with two rectangular tori, one tilted torus and gauge coupling unification, the gaugino condensations are allowed, and thus supersymmetry breaking and moduli stabilization are possible for further phenomenological study.

Cite this article

Tianjun Li , Rui Sun , Chi Zhang . Four-family ${ \mathcal N }=1$ supersymmetric Pati–Salam models from intersecting D6-branes[J]. Communications in Theoretical Physics, 2022 , 74(6) : 065201 . DOI: 10.1088/1572-9494/ac6747

1. Introduction

One of the main motivations of string phenomenology is to find a unifying ${ \mathcal N }=1$ supersymmetric quantum field theory, a competent framework among various extensions of the four-dimensional Standard Model (SM). D-branes play an important role in constructing interesting models at the phenomenological level, especially in Type I, Type IIA and Type IIB string theories. Chiral fermions appear at

worldvolume singularities of D-branes [28]; and

intersecting loci of D-branes in the internal space [9].

Intersecting D6-branes on Type IIA orientifolds have been used to construct three-family non-supersymmetric models and grand unified models [1030]. Even though these models satisfy the Ramond–Ramond (RR) tadpole cancellation conditions, there are Neveu–Schwarz–Neveu–Schwarz (NS-NS) tadpoles remaining due to their non-supersymmetric nature. Moreover, the string scale is close to the Planck scale since the intersecting D6-branes are not transversal in the internal space. As a result, there are large Planck scale corrections at the loop level, leading to the gauge hierarchy problem. As a remedy, a large number of supersymmetric standard-like models and grand unified models [3149] have been constructed, with the gauge hierarchy problem solved. We refer to [50] for a comprehensive review of such non-supersymmetric and supersymmetric models.
Among these supersymmetric models, Pati–Salam models have been a prominent road to the Standard Model, without adding any extra U(1) symmetry around the electroweak level. In [40], Cvetič Liu and one of us (TL) showed how to systematically construct ${ \mathcal N }=1$ supersymmetric Pati–Salam models from intersecting D6-branes on Type IIA ${{\mathbb{T}}}^{6}/{{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$-orientifold. After D-brane splitting and supersymmetric preserving Higgs mechanism applied, the Pati–Salam gauge symmetry SU(4)C × SU(2)L × SU(2)R eventually breaks down to the SM gauge symmetry. Due to the supersymmetry breaking and moduli stabilization triggered by their two confining groups in the hidden sectors, these models do have realistic and phenomenological consequences, as shown in [43, 51, 52]. Until very recently, intriguingly all possible 202752 three-family ${ \mathcal N }=1$ supersymmetric Pati–Salam models on ${{\mathbb{T}}}^{6}/{{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$ have been found and classified with 33 types of independent models according to gauge coupling relations [1].
Most of these above achievements are based on three-family model buildings. In fact, the study of the SM with four-families is also worthy of attention due to the flavor democracy hypothesis. In [53], a democratic mass matrix model was introduced to mainly fix the mass gap problem between three families of fermions of the SM, as well as the hierarchy problem of the Yukawa couplings. Allowing three families of fermions, one gets typical SM predictions, such as a low mass of the top quark and the inequality between three neutrino masses [54]. If one allows four families of fermions in the democratic mass matrix model, three families of precisely massless neutrinos and a massive neutrino can be realized without the assumption of a larger hierarchy of the Yukawa couplings. Moreover, via a slight breaking of democracy, the three massless neutrinos obtain small masses. Then the flavor problem of the SM can be solved naturally by putting the flavor democracy hypothesis due to the important role the democratic mass matrix model plays.
As a matter of fact, the SM does not make a theoretical prediction on the number of families. The only restriction for the number of SM families comes from the requirement made by the quantum chromodynamics (QCD). The asymptotic freedom of QCD provides an upper bound of 8 for the number of families as discussed in [55]. And if we allow four families of fermions for SM, many open issues of SM can be solved in a natural way. For example, introducing a fourth massive generation to SM can alter the cross-section and decay channels of the Higgs particles [56]. When the Yukawa couplings of the fourth generation particles are large enough, these particles are natural candidates for electroweak symmetry breaking [57]. Motivated by constructing realistic four-family SM, in this paper, we concentrate on building four-family ${ \mathcal N }=1$ supersymmetric models with gauge symmetry SU(4)C × SU(2)L × SU(2)R on ${{\mathbb{T}}}^{6}/{{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$ orientifold with intersecting D6-branes as an extension to [1, 40]. We obtain four-family Pati–Salam models with gauge coupling unification at string scale or near string scale. In particular, there are models without any filler brane required as well.
The paper is organized as follows. In section 2 we will briefly review the basics of model building from intersecting D6-branes on ${{\mathbb{T}}}^{6}/{{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$ orientifold. Constraints on D6-brane configuration, such as RR tadpole cancellation condition and supersymmetric condition are also reviewed. In section 3, we present the symmetry-breaking mechanism for U(4)C × U(2)L × U(2)R, i.e. D6-brane splitting and the Higgs mechanism. Four-family of chiral fermion conditions and various symmetry relations, such as T-dualities are also discussed. Section 4 is devoted to the phenomenological features of the four-family models. For each class of model, we list its chiral spectrum for the open string sector. In section 5, we draw conclusions and briefly discuss the limitations of our work. The four-family Pati–Salam models are presented in the appendix.

2. Basics of ${{\mathbb{T}}}^{6}/({{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2})$-orientifolds model construction

To construct realistic and four-family supersymmetric models, we recall the basics of model construction from Type IIA string theory compactified on ${{\mathbb{T}}}^{6}/({{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2})$-orientifold, with D6-branes intersecting at general tilted angles, under similar settings as in [31] and [33].
We begin with the orientifold ${{\mathbb{T}}}^{6}/({{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2})$, where D6-branes can be naturally viewed as general 3-cycles. Using the canonical isomorphism ${{\mathbb{T}}}^{6}\simeq {{\mathbb{T}}}^{2}\times {{\mathbb{T}}}^{2}\times {{\mathbb{T}}}^{2}$, we can easily write down the coordinate chart of the orientifold. Let zi, i = 1, 2, 3 be the complex coordinates of the ith torus, respectively. Let θ and ω be the two generators of the abelian group ${{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$. In coordinate (z1, z2, z3), we define orientifold actions θ and ω of ${{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$ on ${{\mathbb{T}}}^{6}$ by
$\begin{eqnarray}\begin{array}{l}\theta ({z}_{1},{z}_{2},{z}_{3}):= (-{z}_{1},-{z}_{2},{z}_{3}),\\ \omega ({z}_{1},{z}_{2},{z}_{3}):= ({z}_{1},-{z}_{2},-{z}_{3}).\end{array}\end{eqnarray}$
In addition, we define actions Ω and R of ${{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$ on ${{\mathbb{T}}}^{6}$, by Ω the parity-reversion on the world-sheet, and R the complex conjugate of ${{\mathbb{T}}}^{6}$ as a complex manifold. By the very definition of an orientifold, where the world-sheet parity and complex conjugate make no difference, there are exactly four components of the ${{\mathbb{T}}}^{6}/({{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2})$-orientifold, namely the image of ${{\mathbb{T}}}^{6}$ under the action of ΩR, ΩRθ, ΩRω and ΩRθω. These components are objects bearing RR charges, and thus D6-branes are introduced to cancel their RR charges.
Generally, Dp-branes are (p + 1)-dimensional objects in the spacetime, where strings start from and land on. In our case, viewed from the internal space ${{\mathbb{T}}}^{6}/{{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$, the D6-branes are three-dimensional objects. It is sufficient to identify D6-branes and 3-cycles for physical considerations. The powerful Eilenberg–Zilber theorem tells us that
$\begin{eqnarray*}\begin{array}{l}{H}_{3}({{\mathbb{T}}}^{6};{\mathbb{Z}})\simeq {H}_{3}({{\mathbb{T}}}^{2}\times {{\mathbb{T}}}^{2}\times {{\mathbb{T}}}^{2};{\mathbb{Z}})\\ \,\,\,\simeq {H}_{1}({{\mathbb{T}}}^{2};{\mathbb{Z}})\times {H}_{1}({{\mathbb{T}}}^{2};{\mathbb{Z}})\times {H}_{1}({{\mathbb{T}}}^{2};{\mathbb{Z}})\\ \,\,\,\simeq {{\mathbb{Z}}}^{2}\times {{\mathbb{Z}}}^{2}\times {{\mathbb{Z}}}^{2}.\end{array}\end{eqnarray*}$
In the following discussion, we denote [ai], [bi], i = 1, 2, 3 as generators of ${H}_{3}({{\mathbb{T}}}^{2};{\mathbb{Z}})$ for the ith torus respectively. Under the orientifold action ${{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$ and taking world-sheet into consideration, we find that there are only two patterns for the lattice ${H}_{1}({{\mathbb{T}}}^{2};{\mathbb{Z}})\simeq {{\mathbb{Z}}}^{2}$: the rectangular one and the tilted one [11, 33, 35, 45]. In the basis [ai], [bi], i = 1, 2, 3, the former is generated by [ai], [bi] over ${\mathbb{Z}}$, while the latter is generated by $[{\tilde{a}}_{i}],[{b}_{i}]$ over ${\mathbb{Z}}$ with $[{\tilde{a}}_{i}]:= [{a}_{i}]+\tfrac{1}{2}[{b}_{i}]$. Under the basis [a1], [b1], [a2], [b2], [a3], [b3], we can represent a 3-cycle [Πa] of the orientifold ${{\mathbb{T}}}^{6}/({{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2})$ by the coordinate $({n}_{a}^{1},{2}^{-{\beta }_{1}}{l}_{a}^{1})\times ({n}_{a}^{2},{2}^{-{\beta }_{2}}{l}_{a}^{2})\times ({n}_{a}^{3},{2}^{-{\beta }_{3}}{l}_{a}^{3})$, where ${n}_{a}^{i},{l}_{a}^{i}$ are all integers, βi takes 0 when the ith torus is rectangular and takes value 1 for the tilted case. Under the action of ΩR, a D6-brane $[{{\rm{\Pi }}}_{a}]=({n}_{a}^{1},{l}_{a}^{1})\times ({n}_{a}^{2},{l}_{a}^{2})\times ({n}_{a}^{3},{l}_{a}^{3})$ becomes $[{{\rm{\Pi }}}_{a^{\prime} }]=({n}_{a}^{1},-{l}_{a}^{1})\times ({n}_{a}^{2},-{l}_{a}^{2})\times ({n}_{a}^{3},-{l}_{a}^{3})$, with $[{{\rm{\Pi }}}_{a^{\prime} }]$ a short-handed notation for the image of [Πa]. To sum up, a D6-brane [Πa] and its orientifold image $[{{\rm{\Pi }}}_{a^{\prime} }]$ take the form
$\begin{eqnarray}\begin{array}{rcl}\left[{{\rm{\Pi }}}_{a}\right] & = & ({n}_{a}^{1},{2}^{-{\beta }_{1}}{l}_{a}^{1})\times ({n}_{a}^{2},{2}^{-{\beta }_{2}}{l}_{a}^{2})\\ & & \times ({n}_{a}^{3},{2}^{-{\beta }_{3}}{l}_{a}^{3})\\ \left[{{\rm{\Pi }}}_{a^{\prime} }\right] & = & ({n}_{a}^{1},-{2}^{-{\beta }_{1}}{l}_{a}^{1})\times ({n}_{a}^{2},-{2}^{-{\beta }_{2}}{l}_{a}^{2})\\ & & \times ({n}_{a}^{3},-{2}^{-{\beta }_{3}}{l}_{a}^{3}).\end{array}\end{eqnarray}$
Denote [ΠΩR] = 23(1, 0) × (1, 0) × (1, 0), then we have
$\begin{eqnarray}\begin{array}{rcl}\left[{{\rm{\Pi }}}_{{\rm{\Omega }}R}\right] & = & {2}^{3}(1,0)\times (1,0)\times (1,0),\\ \left[{{\rm{\Pi }}}_{{\rm{\Omega }}R\omega }\right] & = & -{2}^{3-{\beta }_{2}-{\beta }_{3}}(1,0)\times (0,1)\times (0,1),\\ \left[{{\rm{\Pi }}}_{{\rm{\Omega }}R\theta \omega }\right] & = & -{2}^{3-{\beta }_{1}-{\beta }_{3}}(0,1)\times (1,0)\times (0,1),\\ \left[{{\rm{\Pi }}}_{{\rm{\Omega }}R\omega \theta }\right] & = & -{2}^{3-{\beta }_{1}-{\beta }_{2}}(0,1)\times (0,1)\times (0,1),\end{array}\end{eqnarray}$
where [ΠΩRω], [ΠΩRθ], [ΠΩRθω] are the images of [ΠΩR] under the action of ω, θ and θω respectively. The coefficient 23 comes from identifying 8 O6-branes (±1, 0) × (±1, 0) × (±1, 0), under the quotient of Ω and R actions. The intersection number of D6-branes can be easily computed as
$\begin{eqnarray}\begin{array}{l}{I}_{{ab}}=[{{\rm{\Pi }}}_{a}]\,[{{\rm{\Pi }}}_{b}]={2}^{-k}({n}_{a}^{1}\,{l}_{b}^{1}-{n}_{b}^{1}\,{l}_{a}^{1})\\ \,\times \,({n}_{a}^{2}\,{l}_{b}^{2}-{n}_{b}^{2}\,{l}_{a}^{2})({n}_{a}^{3}\,{l}_{b}^{3}-{n}_{b}^{3}\,{l}_{a}^{3}),\\ {I}_{{ab}^{\prime} }=[{{\rm{\Pi }}}_{a}]\,[{{\rm{\Pi }}}_{b^{\prime} }]=-{2}^{-k}({n}_{a}^{1}\,{l}_{b}^{1}+{n}_{b}^{1}\,{l}_{a}^{1})\\ \,\times \,({n}_{a}^{2}\,{l}_{b}^{2}+{n}_{b}^{2}\,{l}_{a}^{2})({n}_{a}^{3}\,{l}_{b}^{3}+{n}_{b}^{3}\,{l}_{a}^{3}),\\ {I}_{{aa}^{\prime} }=[{{\rm{\Pi }}}_{a}]\,[{{\rm{\Pi }}}_{a^{\prime} }]=-{2}^{3-k}{n}_{a}^{1}\,{l}_{a}^{1}\,{n}_{a}^{2}\,{l}_{a}^{2}\,{n}_{a}^{3}\,{l}_{a}^{3},\end{array}\end{eqnarray}$
where k = β1 + β2 + β3. If we denote [ΠO6] = [ΠΩR] +[ΠΩRω] + [ΠΩRθ] + [ΠΩRθω], we have
$\begin{eqnarray}\begin{array}{rcl}{I}_{a\,O6} & = & [{{\rm{\Pi }}}_{a}]\,[{{\rm{\Pi }}}_{O6}]={2}^{3-k}(-{l}_{a}^{1}\,{l}_{a}^{2}\,{l}_{a}^{3}\\ & & +{l}_{a}^{1}\,{n}_{a}^{2}\,{n}_{a}^{3}+{n}_{a}^{1}\,{l}_{a}^{2}\,{n}_{a}^{3}+{n}_{a}^{1}\,{n}_{a}^{2}\,{l}_{a}^{3}).\end{array}\end{eqnarray}$
The spectrum of intersecting D6-branes is given in table 1.
Table 1. Spectrum of intersecting D6-branes.
Sectors Representations
aa U(Na/2) vector multiplets
3 adjoint chiral multiplets
ab + ba ${I}_{{ab}}({\square }_{a},\overline{{\square }_{b}})$ fermions
${ab}^{\prime} +b^{\prime} a$ ${I}_{{ab}^{\prime} }({\square }_{a},{\square }_{b})$ fermions
${aa}^{\prime} +a^{\prime} a$ $\tfrac{1}{2}({I}_{{aa}^{\prime} }-\tfrac{1}{2}{I}_{a,O6})\,\square \square $ fermions
$\tfrac{1}{2}({I}_{{aa}^{\prime} }+\tfrac{1}{2}{I}_{a,O6})$ fermions
For three-family ${ \mathcal N }=1$ supersymmetric Pati–Salam model building in Type IIA orientifolds on ${{\mathbb{T}}}^{6}/({{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2})$ with intersecting D6-branes in which the SU(4)C × SU(2)L × SU(2)R gauge symmetries arise from U(n) branes. To get four families of fermions, we require
$\begin{eqnarray}\begin{array}{l}{I}_{{ab}}+{I}_{{ab}^{\prime} }=4,\\ {I}_{{ac}}=-4,{I}_{{ac}^{\prime} }=0.\end{array}\end{eqnarray}$
The conditions ${I}_{{ab}}+{I}_{{ab}^{\prime} }=4$ and Iac = −4 ensure that there are four families of SM fermions, and the condition ${I}_{{ac}^{\prime} }=0$ means that stack a of D6-branes is parallel to the orientifold image of stack c of D6-branes. Thus there should be open strings stretching between the two stack of D6-branes. The light scalar from the NS scalar will obtain mass ${Z}_{{ac}^{\prime} }^{2}/4\pi \alpha ^{\prime} $, with ${Z}_{{ac}^{\prime} }^{2}$ the minimal squared length of the stretching string. Similarly, the light fermions from the R sector acquire the same masses [13, 14, 39]. These light scalars and fermions form the Higgs fields needed to break the Pati–Salam gauge symmetry to the SM gauge symmetry.
In addition, there are two main constraints for D6-brane configurations and O6-plane configurations, namely the RR Tadpole Cancellation Condition and Supersymmetry Condition. As the sources of RR charges, D6-branes and O6-planes should satisfy the Gauss'law, for the flux of RR fields through the compact space ${{\mathbb{T}}}^{6}/({{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2})$ without boundary should be conserved. This form of Gauss'law is the so-called RR tadpole cancellation condition. To satisfy this condition, stacks of Na, a = 1, 2, 3 D6-branes are often needed to be introduced as the so-called filler brane, with a running through three families of gauge groups. Then the RR tadpole cancellation condition reads
$\begin{eqnarray}\sum _{a=1}^{3}{N}_{a}[{{\rm{\Pi }}}_{a}]+\sum _{a=1}^{3}{N}_{a}[{{\rm{\Pi }}}_{a^{\prime} }]-4[{{\rm{\Pi }}}_{O6}]=0,\end{eqnarray}$
and the coefficient 4 before $[{{\rm{\Pi }}}_{O6}]$ comes from the −4 RR charges in the D6-brane charge unit. Other than O6-planes in equation (2.7), we can also introduce D6-branes between the four O6-planes to cancel out the total RR charges, and rewrite equation (2.7) free of O6-planes. To do this, we note
$\begin{eqnarray}\begin{array}{l}{A}_{a}=-{n}_{a}^{1}\,{n}_{a}^{2}\,{n}_{a}^{3},{B}_{a}={n}_{a}^{1}\,{l}_{a}^{2}\,{l}_{a}^{3},\\ {C}_{a}={l}_{a}^{1}\,{n}_{a}^{2}\,{l}_{a}^{3},{D}_{a}={l}_{a}^{1}\,{l}_{a}^{2}\,{n}_{a}^{3};\\ {\tilde{A}}_{a}=-{l}_{a}^{1}\,{l}_{a}^{2}\,{l}_{a}^{3},{\tilde{B}}_{a}={l}_{a}^{1}\,{n}_{a}^{2}\,{n}_{a}^{3},\\ {\tilde{C}}_{a}={n}_{a}^{1}\,{l}_{a}^{2}\,{n}_{a}^{3},{\tilde{D}}_{a}={n}_{a}^{1}\,{n}_{a}^{2}\,{l}_{a}^{3}.\end{array}\end{eqnarray}$
Therefore the tadpole cancellation conditions 2.7 can be expressed in the form
$\begin{eqnarray}\begin{array}{l}-{2}^{k}{N}^{(1)}+\sum _{a=1}^{3}{N}_{a}{A}_{a}=-16,\\ -{2}^{k}{N}^{(2)}+\sum _{a=1}^{3}{N}_{a}{B}_{a}=-16,\\ -{2}^{k}{N}^{(3)}+\sum _{a=1}^{3}{N}_{a}{C}_{a}=-16,\\ -{2}^{k}{N}^{(4)}+\sum _{a=1}^{3}{N}_{a}{D}_{a}=-16,\end{array}\end{eqnarray}$
where 2N(i), i = 1, 2, 3, 4 are the number of filler branes between the four O6-branes, as shown in table 2.
Table 2. Configuration of four O6-planes.
Orientifold actions O6-planes (n1, l1) × (n2, l2) × (n3, l3)
ΩR 1 $({2}^{{\beta }_{1}},0)\times ({2}^{{\beta }_{2}},0)\times ({2}^{{\beta }_{3}},0)$
ΩRω 2 $({2}^{{\beta }_{1}},0)\times (0,-{2}^{{\beta }_{2}})\times (0,{2}^{{\beta }_{3}})$
ΩRθω 3 $(0,-{2}^{{\beta }_{1}})\times ({2}^{{\beta }_{2}},0)\times (0,{2}^{{\beta }_{3}})$
ΩRθ 4 $(0,-{2}^{{\beta }_{1}})\times (0,{2}^{{\beta }_{2}})\times ({2}^{{\beta }_{3}},0)$
Furthermore, the four-dimensional ${ \mathcal N }=1$ supersymmetric conditions require 1/4 supercharges conserved under (i) orientation reversion of D6-branes; and (ii) group action of ${{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$. As shown in [9], the four-dimensional ${ \mathcal N }=1$ supersymmetry is preserved under the orientation reversion if and only if rotations of D6-branes with respect to O6-planes are elements of SU(3), while their total rotation angles equal to 0. When the four-dimensional ${ \mathcal N }=1$ supersymmetry is preserved under orientation reversion, it will be preserved under the ${{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$-action manifestly. The supersymmetric condition can be written as [33]
$\begin{eqnarray}\begin{array}{l}{x}_{A}{\tilde{A}}_{a}+{x}_{B}{\tilde{B}}_{a}+{x}_{C}{\tilde{C}}_{a}+{x}_{D}{\tilde{D}}_{a}=0,\\ {A}_{a}/{x}_{A}+{B}_{a}/{x}_{B}+{C}_{a}/{x}_{C}\\ +\,{D}_{a}/{x}_{D}\lt 0,a\,=\,1,2,3,\end{array}\end{eqnarray}$
with ${x}_{A}=\lambda ,{x}_{B}=\lambda {2}^{{\beta }_{2}+{\beta }_{3}}{/\chi }_{2}{\chi }_{3},{x}_{D}=\lambda {2}^{{\beta }_{1}+{\beta }_{2}}/{\chi }_{1}{\chi }_{2}$. The positive parameter λ is introduced to put equation (2.10) on an equal scaling, and ${\chi }_{i}={R}_{i}^{2}/{R}_{i}^{1},i=1,2,3$ is the complex structure modulus parameter for the ith torus respectively.

3. Gauge symmetry breaking via brane splittings

To obtain SM or standard-like models via the mechanism of intersecting D6-branes, there should be at least two extra U(1) gauge symmetries for either supersymmetric models or non-supersymmetric models, as a result of the constraints on the quantum number of the right handed electron [14, 3133]. Among these two U(1) gauge symmetries, one is lepton number symmetry U(1)L and another ${\rm{U}}{(1)}_{{I}_{3R}}$ is an analogy for right-hand weak isospin. We have the hypercharge QY expressed in the form
$\begin{eqnarray}{Q}_{Y}={Q}_{{I}_{3R}}+\displaystyle \frac{{Q}_{B}-{Q}_{L}}{2}.\end{eqnarray}$
The baryonic charge QB arises from U(1)B, via the decomposition U(3)C ≃ SU(3)C × U(1)B. On the other hand, since the U(1)I3R gauge field should be massless, the gauge group U(1)I3R must come from the non-abelian component of U(2)R or USp symmetry, otherwise the U(1)I3R will acquire mass from the BF couplings. To get an anomaly-free U(1)BL, the U(1)L symmetry should come from some non-abelian group for similar reasons. In previous studies of supersymmetric model building [31, 32], U(1)I3R comes from USp groups. These models indeed have two extra anomaly-free U(1) symmetries, and have at least 8 Higgs doublets. One could in principle break their symmetry groups down to the SM symmetry, but cannot do this without violating the D-flatness and F-flatness, thus the supersymmetry.
In this paper, as introduced in [52], we study a generalized version of the four-family MSSM models. In these models, the b or c stacks are with twice the number of D6-branes. Then the gauge symmetries of these generalized four-family models can be broken to the standard four-family gauge symmetries SU(4)C × SU(2)L × SU(2)R via the Higgs mechanism. Taking gauge symmetries ${\rm{U}}{(4)\times {\rm{U}}{(4)}_{L}\times {\rm{U}}(2)}_{R}$ as example, we consider a U(4) gauge theory with a scalar field in the adjoint representation. By choosing appropriate rotations commuting with the generators of the Lie algebra of SU(4), one can break U(4) to U(2) × U(2), and finally to U(2). We choose the rotation for U(4) scalar field acting on the vacuum expectation value Φ0 as
$\begin{eqnarray}{{\rm{\Phi }}}_{1}=\left(\begin{array}{cccc}{V}_{1} & 0 & 0 & 0\\ 0 & {V}_{1} & 0 & 0\\ 0 & 0 & -{V}_{1} & 0\\ 0 & 0 & 0 & -{V}_{1}\end{array}\right).\end{eqnarray}$
We find that the U(4) gauge symmetry breaks spontaneously to U(2) × U(2), as the matrix in (3.2) lies in the center of the Lie algebra of U(4). What is left to us is to break down U(2) × U(2) to U(2). For each U(2) component in U(2) × U(2), the generators of its Lie algebra are the standard Pauli matrices. If we choose the rotation matrix for U(2) × U(2) as
$\begin{eqnarray}{{\rm{\Phi }}}_{2}=\left(\begin{array}{cccc}0 & 0 & {V}_{2} & 0\\ 0 & 0 & 0 & {V}_{2}\\ {V}_{2} & 0 & 0 & 0\\ 0 & {V}_{2} & 0 & 0\end{array}\right),\end{eqnarray}$
we find that U(2) × U(2) breaks down to U(2). We also note that a mass of
$\begin{eqnarray}{m}^{2}=8{g}^{2}{V}_{2}^{2}\end{eqnarray}$
is acquired through the above process. For models with gauge symmetry ${\rm{U}}{(4)\times {\rm{U}}{(2)}_{L}\times {\rm{U}}(4)}_{R}$ and ${\rm{U}}{(4)\times {\rm{U}}{(4)}_{L}\times {\rm{U}}(4)}_{R}$, one can similarly break the symmetry down to ${\rm{U}}{(4)\times {\rm{U}}{(2)}_{L}\times {\rm{U}}(2)}_{R}$ following the above procedure. The anomalies of the overall U(1) symmetries are canceled by the generalized Green-Schwarz mechanism [13, 14, 31], while their fields get massive from the linear BF couplings.
The gauge symmetry SU(4)C × SU(2)L × SU(2)R can be further broken down to SM gauge symmetry by D6-brane splitting and Higgs mechanism. Firstly, one can split stack a of Na = 8 D6-branes into stack a1 of ${N}_{{a}_{1}}=6$ D6-branes and stack a2 ${N}_{{a}_{2}}=2$ D6-branes. Then the U(4)C symmetry breaks down to U(3) × U(1). We denote the numbers of symmetric and anti-symmetric representations for SU(4)C , SU(2)L and SU(2)R by ${n}_{\square \square }^{a}$ and , ${n}_{\square \square }^{b}$, and ${n}_{\square \square }^{c}$, . After splitting, the symmetric and anti-symmetric representations of SU(4)C descend to symmetric representations of SU(3)C and U(1)BL, and anti-symmetric representations of SU(3)C. Note that there are ${I}_{{a}_{1}a{{\prime} }_{2}}$ new fields, arising from the intersection of a1 stack and a2 stack of D6-branes. The anomaly-free gauge symmetry is SU(3)C × U(1)BL, a subgroup of SU(4)C.
Similarly, the stack c of Nc = 4 D6-branes can be broken into stack c1 of ${N}_{{c}_{1}}=2$ D6-branes and stack c2 ${N}_{{c}_{2}}=2$ D6-branes. Then the U(2)R symmetry breaks down to U(1)I3R. The symmetric representations of SU(2)R descend to the symmetric representations of U(1)I3R only. Also, there are ${I}_{{c}_{1}c{{\prime} }_{2}}$ new fields, arising from the intersection of c1 stack and c2 stack of D6-branes. The anomaly-free gauge symmetry is U(1)I3R, a subgroup of SU(2)R. After splitting, the gauge symmetry of our model breaks down to SU(3)C ×SU(2)L × U(1)BL × U(1)I3R .
To get just the SM gauge symmetry, we assume the minimal squared distance ${Z}_{{a}_{2}c{{\prime} }_{1}}^{2}$ between a2 stack and the orientifold image of c1 stack of D6-branes to be very small. Then there are ${I}_{{a}_{2}c{{\prime} }_{1}}$ chiral multiplets of light fermions, arising from the open string stretching between a2 stack of D6-branes and the orientifold image of c1 stack of D6-branes. These particles break down SU(3)C × SU(2)L × U(1)BL × U(1)I3R to the SM gauge symmetry, playing the same role as the right-handed neutrinos and their complex conjugates. Meanwhile, they preserve the D-flatness and F-flatness, thus the supersymmetry. In conclusion, the whole symmetry breaking chain is
$\begin{eqnarray}\begin{array}{l}\mathrm{SU}{(4)}_{C}\times \mathrm{SU}{(2)}_{L}\times \mathrm{SU}{(2)}_{R}\\ \mathop{\longrightarrow }\limits^{a\to {a}_{1}+{a}_{2}}\mathrm{SU}{(3)}_{C}\times \mathrm{SU}{(2)}_{L}\times \mathrm{SU}{(2)}_{R}\times U{(1)}_{B-L}\\ \mathop{\longrightarrow }\limits^{c\to {c}_{1}+{c}_{2}}\mathrm{SU}{(3)}_{C}\times \mathrm{SU}{(2)}_{L}\\ \times {\rm{U}}{(1)}_{I3R}\times {\rm{U}}{(1)}_{B-L}\\ \mathop{\longrightarrow }\limits^{{\rm{Higgs}}\,{\rm{Mechanism}}}\mathrm{SU}{(3)}_{C}\\ \times \mathrm{SU}{(2)}_{L}\times {\rm{U}}{(1)}_{Y}.\end{array}\end{eqnarray}$
The process of dynamical supersymmetry breaking has been studied in [36] for D6-brane models from Type IIA orientifolds. The kinetic function for a stack a of D6-branes is of the form [43]
$\begin{eqnarray}\begin{array}{rcl}{f}_{a} & = & \displaystyle \frac{1}{4{\kappa }_{a}}\left({n}_{a}^{1}\,{n}_{a}^{2}\,{n}_{a}^{3}s-\displaystyle \frac{{n}_{a}^{1}\,{l}_{a}^{2}\,{l}_{a}^{3}\,{u}^{1}}{{2}^{{\beta }_{2}+{\beta }_{3}}}\right.\\ & & \left.-\displaystyle \frac{{l}_{a}^{1}\,{n}_{a}^{2}\,{l}_{a}^{3}\,{u}^{2}}{{2}^{{\beta }_{1}+{\beta }_{3}}}-\displaystyle \frac{{l}_{a}^{1}\,{l}_{a}^{2}\,{n}_{a}^{3}{u}^{3}}{{2}^{{\beta }_{1}+{\beta }_{2}}}\right),\end{array}\end{eqnarray}$
where κa is a constant with respect to the gauge groups, for instance, κa = 1 for SU(Na). We use moduli parameters s and ui, i = 1, 2, 3 in a supergravity basis, which are related to four-dimensional dilation parameter φ4 and complex structure moduli parameters Ui, i = 1, 2, 3 as follows
$\begin{eqnarray}\begin{array}{rcl}\mathrm{Re}(s) & = & \displaystyle \frac{{{\rm{e}}}^{-{\phi }_{4}}}{2\pi }\displaystyle \frac{\sqrt{\mathrm{Im}({U}^{1})\,\mathrm{Im}({U}^{2})\,\mathrm{Im}({U}^{3})}}{\left|{U}^{1}{U}^{2}{U}^{3}\right|},\\ \mathrm{Re}({u}^{1}) & = & \displaystyle \frac{{{\rm{e}}}^{-{\phi }_{4}}}{2\pi }\sqrt{\displaystyle \frac{\mathrm{Im}({U}^{1})}{\mathrm{Im}({U}^{2})\,\mathrm{Im}({U}^{3})}}\left|\displaystyle \frac{{U}^{2}{U}^{3}}{{U}^{1}}\right|,\\ \mathrm{Re}({u}^{2}) & = & \displaystyle \frac{{{\rm{e}}}^{-{\phi }_{4}}}{2\pi }\sqrt{\displaystyle \frac{\mathrm{Im}({U}^{2})}{\mathrm{Im}({U}^{1})\,\mathrm{Im}({U}^{3})}}\left|\displaystyle \frac{{U}^{1}{U}^{3}}{{U}^{2}}\right|,\\ \mathrm{Re}({u}^{3}) & = & \displaystyle \frac{{{\rm{e}}}^{-{\phi }_{4}}}{2\pi }\sqrt{\displaystyle \frac{\mathrm{Im}({U}^{3})}{\mathrm{Im}({U}^{1})\,\mathrm{Im}({U}^{2})}}\left|\displaystyle \frac{{U}^{1}{U}^{2}}{{U}^{3}}\right|.\end{array}\end{eqnarray}$
In our present models, the Ui, i = 1, 2, 3 can be computed as in [45]
$\begin{eqnarray}\begin{array}{rcl}{U}^{1} & = & {\rm{i}}\,{\chi }_{1},{U}^{2}={\rm{i}}{\chi }_{2},\\ {U}^{3} & = & \displaystyle \frac{2\,{\chi }_{3}^{2}+4\,{\rm{i}}{\chi }_{3}}{4+{\chi }_{3}^{2}}.\end{array}\end{eqnarray}$
Moreover, the Kähler potential takes the form of
$\begin{eqnarray}K=-\mathrm{ln}(S+\overline{S})-\sum _{i=1}^{3}\mathrm{ln}({U}^{i}+{\overline{U}}^{i}).\end{eqnarray}$
Note that the three moduli parameters χ1, χ2, χ3 are not independent, as they can be expressed in terms of xA, xB, xC, xD and the latter parameters are related by the supersymmetric condition (2.10). Actually, one can determine χ1, χ2, χ3 up to an overall coefficient, namely an action of dilation on these parameters. So one has to stabilize this dilation to determine all the moduli parameters. Previous studies [5860] employ mechanisms like gaugino condensation to stabilize this overall coefficient, dictating that there should be at least two USp groups in the hidden sectors. Moreover, the one-loop beta functions [40]
$\begin{eqnarray}\begin{array}{rcl}{\beta }_{i}^{g} & = & -3\left(\displaystyle \frac{{N}^{(i)}}{2}+1\right)+2\left|{I}_{{ai}}\right|+\left|{I}_{{bi}}\right|\\ & & +\left|{I}_{{ci}}\right|+3\left(\displaystyle \frac{{N}^{(i)}}{2}-1\right)\\ & = & -6+2\left|{I}_{{ai}}\right|+\left|{I}_{{bi}}\right|+\left|{I}_{{ci}}\right|\end{array}\end{eqnarray}$
for each USp(N(i)) arising from 2N(i) filler branes are required to be negative. However, in this paper, to include other potential mechanisms, we do not restrict ourselves only to models with at least two USp groups in the hidden sectors. The gauge coupling constant related to stack a of D6-branes is
$\begin{eqnarray}{g}_{a}^{-2}=\left|\mathrm{Re}({f}_{a})\right|,\end{eqnarray}$
and the coupling constant of stack b and stack c of D6-branes are determined in the same way. The kinetic function for U(1)Y is a linear combination of those for SU(4)C and SU(2)R, as shown in [11, 43]
$\begin{eqnarray}{f}_{Y}=\displaystyle \frac{3}{5}\left(\displaystyle \frac{2}{3}{f}_{a}+{f}_{c}\right).\end{eqnarray}$
The coupling constant gY is determined by
$\begin{eqnarray}{g}_{Y}^{-2}=\left|\mathrm{Re}({f}_{Y})\right|.\end{eqnarray}$
At tree-level, the gauge couplings have the relation
$\begin{eqnarray}{g}_{a}^{2}=\alpha {g}_{b}^{2}=\beta \displaystyle \frac{5}{3}{g}_{Y}^{2}=\gamma (\pi {{\rm{e}}}^{-{\phi }_{4}}),\end{eqnarray}$
where α, β and γ are ratios between the strong coupling and the weak coupling, and hypercharge coupling, respectively.

3.1. T-duality and its variations

In string theory, if two models are related by T-duality, these models are considered equivalent. By applying T-duality, one can tremendously simplify the process of searching inequivalent models. Before the discussion of T-duality, we first point out two obvious symmetries that relate equivalent models.
i

(i)Two models are equivalent if they are related by a permutation of three ${{\mathbb{T}}}^{2};$ and

ii

(ii)Two D6-models are equivalent if their wrapping numbers on any two ${{\mathbb{T}}}^{2}$ are in opposite signs, while are the same on the third ${{\mathbb{T}}}^{2}$.

The above two symmetries are known as the D6-brane Sign Equivalent Principle (DSEP). Then, follow the convention of [40], we introduce Type I and Type II T-dualities. Type I duality transformation acts on arbitrary two ${{\mathbb{T}}}^{2}$, say the jth and kth ${{\mathbb{T}}}^{2}$. The wrapping numbers on these tori transform as follows
$\begin{eqnarray}\begin{array}{l}({n}_{a}^{j},{l}_{a}^{j})\mapsto (-{l}_{a}^{j},{n}_{a}^{j}),\\ ({n}_{a}^{k},{l}_{a}^{k})\mapsto ({l}_{a}^{k},-{n}_{a}^{k}),\end{array}\end{eqnarray}$
when Type I T-duality applies. Recall the definitions in equation (2.8), Type I T-duality only makes an exchange between $({A}_{a},{\tilde{A}}_{a})$, $({B}_{a},{\tilde{B}}_{a})$, $({C}_{a},{\tilde{C}}_{a})$, $({D}_{a},{\tilde{D}}_{a})$. Moreover, Type I duality transformation is often combined with the trivial two ${{\mathbb{T}}}^{2}$ exchange, and we call the combination an extended Type I T-duality.
As for Type II T-duality, it acts on all three different ${{\mathbb{T}}}^{2}$. For instance, if we pick the ith, jth and kth ${{\mathbb{T}}}^{2}$, the wrapping numbers on these tori transform as
$\begin{eqnarray}\begin{array}{l}({n}_{a}^{i},{l}_{a}^{i})\mapsto (-{n}_{a}^{i},{l}_{a}^{i}),\\ ({n}_{a}^{j},{l}_{a}^{j})\mapsto ({l}_{a}^{j},{n}_{a}^{j}),\\ ({n}_{a}^{k},{l}_{a}^{k})\mapsto ({l}_{a}^{k},{n}_{a}^{k}).\end{array}\end{eqnarray}$
In [40], Type II T-duality often combines with the interchange between b and c stacks of D6-branes
$\begin{eqnarray}b\leftrightarrow c\end{eqnarray}$
associated to SU(2)L and SU(2)R gauge groups.
If we composite Type I T-duality and DSEP, we will get a variation of Type II T-duality. Under this symmetry transformation, the wrapping numbers on all three ${{\mathbb{T}}}^{2}$ change as
$\begin{eqnarray}\begin{array}{rcl}({n}_{a}^{1},{l}_{a}^{1}) & \mapsto & ({n}_{a}^{1},-{l}_{a}^{1}),\\ ({n}_{a}^{2},{l}_{a}^{2}) & \mapsto & ({n}_{a}^{2},-{l}_{a}^{2}),\\ ({n}_{a}^{3},{l}_{a}^{3}) & \mapsto & ({n}_{a}^{3},-{l}_{a}^{3}),\\ b & \leftrightarrow & c.\end{array}\end{eqnarray}$
It is easy to see that under variation of Type II T-duality, only the signs of ${\tilde{A}}_{a},{\tilde{B}}_{a},{\tilde{C}}_{a},{\tilde{D}}_{a}$ in equation (2.8) change.
However, it is worth mentioning that, the variation (3.18) of Type II T-duality is not an equivalence in our construction of four-family supersymmetric models, if the model is not invariant under SU(2)L and SU(2)R interchange. This observation makes sense at least phenomenologically. For a four-family supersymmetric model, one can obtain a new model by exchanging the b and c stacks of D6-branes associated to the SU(2)L and SU(2)R groups, as the quantum numbers for SU(2)L and SU(2)R in the particle spectrum will exchange, so will the gauge couplings for these two groups at the string level. We will present with examples of this inequivalence in the next section.

3.2. Supersymmetric 4-family models

Employing the deterministic algorithm in [1], we do not restrict the number of USp groups, and consider two cases, one without any titled torus, and the other with the third torus to be tilted, without loss of generality. Note that four-family models are completely different from the three-family models, thus the argument in [40] to exclude which torus is tilted or cannot be directly applied to our case, due to the even number of generations.

3.2.1. Models without tilted torus

We obtain six classes of 274 supersymmetric four-family models with the deterministic algorithm introduced in [1] with representative models presented in appendix A. The classification is based on the gauge groups, T-equivalences and phenomenological considerations such as gauge coupling relations.
Model 14 is a class of its own, as it is the representative model without tilted torus achieving exact gauge coupling unification at the string scale. The Higgs-like particles in this model arise from the intersection at b and $c^{\prime} $ stacks of D6-branes, while the Higgs doublets arise from the massless open string states in a N = 2 subsector and form vector-like pairs.
The second class of models has no USp group, which means the tadpole cancellation conditions are satisfied without any filler brane as a rare case. These models are represented by Model 15 and Model 16, which are independent models of T-equivalence. These models are the first four-family examples having no confining USp groups to achieve approximate gauge coupling unification at the string scale.
The third class of models includes Model 17 and Model 18 with negative β function and positive β function respectively.
The fourth class of models includes Model 1921 with two USp groups. These models are independent under T-duality.
As discussed in [59] there are at least two confining USp groups needed, with negative β function, and thus allow for gaugino condensations, these models would need alternative mechanisms to break the supersymmetry. Furthermore, we can observe from the spectrum tables that Model 15, 17 and 20 do not have the proper Higgs doublets with quantum number $(1,\,2,\,\overline{2})$ under U(4)C × U(2)L × U(2)R gauge symmetry, from neither the b and c or $c^{\prime} $ stacks of brane intersection nor the massless open string states in a N = 2 subsector. Thus, we do not have the SM fermion Yukawa couplings at the renormalizable level which are invariant under the global U(1)C × U(1)L × U(1)R symmetry in these models.
The fifth class of models no less than three USp groups, and includes Models 1925. Among these models, Model 24 and Model 25 are related by T-dualities with b and c stacks of branes swapping. To see this, we show how Model 24 and Model 25 are related. To begin with, a stacks of D6-branes in Model 24 and Model 25 are related by the DSEP:
$\begin{eqnarray*}\begin{array}{l}(-1,0)\times (-1,1)\times (1,1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}\,{\rm{(i)}}}(-1,0)\times (1,1)\times (-1,1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}\,{\rm{(i)}}}(1,0)\times (1,1)\times (1,-1).\end{array}\end{eqnarray*}$
And the b stack of D6-branes in Model 24 are related to the c stack of D6-branes in Model 25, by Type I T-duality (3.15), the DSEP, and the interchange (3.17) of b and c stacks:
$\begin{eqnarray*}\begin{array}{l}(0,1)\times (-1,1)\times (-3,1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}\,{\rm{(i)}}}(0,1)\times (-3,1)\times (-1,1)\\ \mathop{\longrightarrow }\limits^{{\rm{Type}}\,{\rm{I}}\,{\rm{T-duality}}}(0,1)\times (-1,-3)\times (1,1)\\ \mathop{\longrightarrow }\limits^{\text{DSEP (ii) and b \unicode{x02194} c}}(0,1)\times (1,3)\times (-1,-1).\end{array}\end{eqnarray*}$
While the c stack of D6-branes in Model 24 and the b stack of D6-branes in Model 25 are related by Type II T-duality, DSEP and the bc exchange:
$\begin{eqnarray*}\begin{array}{l}(1,2)\times (-1,0)\times (-1,1)\\ \mathop{\longrightarrow }\limits^{{\rm{Type}}\,{\rm{II}}\,{\rm{T-duality}}}(-1,2)\times (0,-1)\times (1,-1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}\,{\rm{(i)}}}(-1,2)\times (1,-1)\times (0,-1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}\,{\rm{(ii)}}\,{\rm{and}}\,{\rm{b\leftrightarrow c}}}(-1,2)\times (-1,1)\times (0,1).\end{array}\end{eqnarray*}$
Even though these two models are related by the generalized T-duality as above, they are not phenomenologically equivalent. Model 24 achieves U(4) and U(2)R unification, while Model 25 has U(4) and U(2)L unification due to b and c stacks swapping. The Higgs particles of Model 24 come from the intersection of b and $c^{\prime} $ stacks of D6-branes. Since all the beta functions of models within this class are negative, we may break the supersymmetry and stabilize the moduli via gaugino condensations.
The sixth class of models are with large wrapping numbers 5, 6, 7, 8, 9, 10, 11, 13, 15, 17, represented by Models 2635 which did not appear in the former search. Three-family models with large wrapping number 5 have been found in [51], but it is the first time to find four-family models with wrapping numbers at this scale.

3.2.2. Models with one tilted torus

Employing the deterministic algorithm, we obtain in total 6 types of gauge coupling relations with represented models presented in appendix B.
Model 36 is a class of its own, as it is the only type of model with one tilted torus achieving exact gauge coupling unification at the string scale. The Higgs particles in this model arise from the massless open string states in a N = 2 subsector and form vector-like pairs because the b stack branes for these models are parallel to c stack brane images on the third two-tori. Since all the beta functions are negative in this model, we can break the supersymmetry and stabilize the moduli via gaugino condensations.
The second class of models includes Model 37 and Model 38, and has no USp group. These models are related by type II T-duality. More specifically, we will below show how they are related by T-dualities explicitly. It is easy to find that the a stacks of both models are related by DSEP and Type II T-duality
$\begin{eqnarray*}\begin{array}{l}(1,-1)\times (-1,0)\times (-1,-1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}}(-1,0)\times (1,-1)\times (-1,-1)\\ \mathop{\longrightarrow }\limits^{{\rm{Type}}\,{\rm{II}}\,{\rm{T-duality}}}(1,0)\times (-1,1)\times (-1,-1).\end{array}\end{eqnarray*}$
Note that the b stack of wrapping numbers of Model 38 are obtained by applying the variation of Type II T-duality and DSEP on the b stack of wrapping numbers of Model 37:
$\begin{eqnarray*}\begin{array}{l}(0,1)\times (-1,-2)\times (1,1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}}(-1,-2)\times (0,1)\times (1,1)\\ \mathop{\longrightarrow }\limits^{{\rm{variation}}\,{\rm{of}}\,{\rm{Type}}\,{\rm{II}}\,{\rm{T-duality}}}(-1,2)\times (0,-1)\times (1,-1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}}(1,-2)\times (0,-1)\times (-1,1).\\ \end{array}\end{eqnarray*}$
The c stacks of Model 37 and Model 38 are related by applying DSEP twice:
$\begin{eqnarray*}\begin{array}{l}(1,0)\times (1,-2)\times (1,1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}}(1,-2)\times (1,0)\times (1,1)\\ \mathop{\longrightarrow }\limits^{{\rm{DSEP}}}(1,-2)\times (-1,0)\times (-1,-1).\end{array}\end{eqnarray*}$
Models 3942 with one USp group in the hidden sector are related by T-dualities in a similar way. The Higgs particles in Model 37 again come from the massless open string states in a N = 2 subsector and form vector-like pairs because the b stack branes for these models are parallel to c stack brane on the third two-tori. Since there are no USp groups in the hidden sectors, gaugino condensations do not work in this case. One needs to stabilize the modulus and break the supersymmetry via a different mechanism. For Model 39, there are four Higgs doublets arising from N = 2 subsectors due to the parallel of b stack branes and c stack brane on the third two-tori. In addition, there are models 43 and 44 with distinct gauge coupling relations at string scale.

4. Phenomenological analysis

4.1. Models without tilted torus

We begin with Model 14. The gauge group of Model 14 is ${\rm{U}}{(4)\times {\rm{U}}{(2)}_{L}\times {\rm{U}}(2)}_{R}\times \mathrm{USp}(2)\times \mathrm{USp}(4)$. We tabulate the full spectrum of chiral particles of Models 14 and 15 in tables 3 and 4 respectively. Interestingly both Model 15 and Model 16 do not have any USp group, and then they do not have any exotic particles charged under USp groups as well. We tabulate the full spectrum of chiral particles of Model 15 below as a representative for this class.
Table 3. Spectrum of chiral particles of Model 14.
Model 14 SU(4)C × SU(2)L × SU(2)R × USp(2) × USp(4) Q4C Q2L Q2R Qem BL Field
ab $4\times (4,\overline{2},1,1,1)$ 1 −1 0 $-\tfrac{1}{3},\tfrac{2}{3},-1,0$ $\tfrac{1}{3},-1$ QL, LL
ac $4\times (\overline{4},1,2,1,1)$ −1 0 1 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QR, LR
${bc}^{\prime} $ $8\times (1,\overline{2},\overline{2},1,1)$ 0 −1 −1 0, ±1 0 $H^{\prime} $
b1 $4\times (1,\overline{2},1,2,1)$ 0 −1 0 $\mp \tfrac{1}{2}$ 0
a4 $2\times (4,1,1,1,\overline{4})$ 1 0 0 $\tfrac{1}{6},-\tfrac{1}{2}$ $\tfrac{1}{3},-1$
c4 $2\times (1,1,\overline{2},1,4)$ 0 0 −1 $\pm \tfrac{1}{2}$ 0
${a}_{\overline{\square \square }}$ $1\times (\overline{10},1,1,1,1)$ −2 0 0 $-\tfrac{1}{3},1$ $-\tfrac{2}{3},2$
1 × (6, 1, 1, 1, 1) 2 0 0 $\tfrac{1}{3},-\tfrac{1}{3},-1$ $\tfrac{2}{3},-2$
b□□ 3 × (1, 3, 1, 1, 1) 0 2 0 0, ±1 0
3 × (1, 1, 1, 1, 1) 0 2 0 0 0
c□□ 1 × (1, 1, 3, 1, 1) 0 0 2 0, ±1 0
1 × (1, 1, 1, 1, 1) 0 0 2 0 0
bc $2\times (1,\overline{2},2,1,1)$ 0 −1 1 0, ±1 0 ${H}_{u}^{i},{H}_{d}^{i}$
$2\times (1,2,\overline{2},1,1)$ 0 1 −1
Table 4. Spectrum of chiral particles of Model 15.
Model 15 SU(4)C × SU(2)L × SU(2)R Q4C Q2L Q2R Qem BL Field
ab $8\times (4,\overline{2},1)$ 1 −1 0 $-\tfrac{1}{3},\tfrac{2}{3},-1,0$ $\tfrac{1}{3},-1$ QL, LL
${ab}^{\prime} $ $4\times (\overline{4},\overline{2},1)$ −1 −1 0 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QL, LL
ac $4\times (\overline{4},1,2)$ −1 0 1 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QR, LR
bc $16\times (1,\overline{2},2)$ 0 −1 1 0 ± 1 0 $H^{\prime} $
b□□ 5 × (1, 3, 1) 0 2 0 0, ±1 0
5 × (1, 1, 1) 0 2 0 0 0
${c}_{\overline{\square \square }}$ $7\times (1,1,\overline{3})$ 0 0 −2 0, ±1 0
9 × (1, 1, 1) 0 0 2 0, ±1 0
${bc}^{\prime} $ $8\times (1,\overline{2},\overline{2})$ 0 −1 −1 0, ±1 0 $H^{\prime} $
8 × (1, 2, 2) 0 1 1
Model 17 and Model 18 are constructed with one USp group. We tabulate the full spectrum of chiral particles of Model 17 in table 5 as a representative of this class of models.
Table 5. Spectrum of chiral particles of Model 17.
Model 17 SU(4)C × SU(2)L × SU(2)R × USp(2) Q4C Q2L Q2R Qem BL Field
ab $6\times (4,\overline{2},1,1)$ 1 −1 0 $-\tfrac{1}{3},\tfrac{2}{3},-1,0$ $\tfrac{1}{3},-1$ QL, LL
${ab}^{\prime} $ $2\times (\overline{4},\overline{2},1,1)$ −1 −1 0 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QL, LL
ac $4\times (\overline{4},1,2,1)$ −1 0 1 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QR, LR
bc $8\times (1,\overline{2},2,1)$ 0 −1 1 0, ±1 0 $H^{\prime} $
c2 $4\times (1,1,2,\overline{2})$ 0 0 1 $\pm \tfrac{1}{2}$ 0
b□□ 3 × (1, 3, 1, 1) 0 2 0 0, ±1 0
3 × (1, 1, 1, 1) 0 2 0 0 0
${c}_{\overline{\square \square }}$ $7\times (1,1,\overline{3},1)$ 0 0 −2 0, ±1 0
9 × (1, 1, 1, 1) 0 0 2 0 0
${bc}^{\prime} $ $6\times (1,\overline{2},\overline{2},1)$ 0 −1 −1 0, ±1 0 $H^{\prime} $
6 × (1, 2, 2, 1) 0 1 1
Models 1921 are built with two USp groups. The gauge groups for Model 19, Model 20 and Model 21 are ${\rm{U}}(4)\,\times {\rm{U}}{(4)}_{L}\times {{\rm{U}}(2)}_{R}\times \mathrm{USp}(2)\times \mathrm{USp}(4)$, ${\rm{U}}(4)\times {\rm{U}}{(2)}_{L}\times {{\rm{U}}(2)}_{R}\times \mathrm{USp}{\left(4\right)}^{2}$ and ${\rm{U}}{(4)\times {\rm{U}}{(4)}_{L}\times {\rm{U}}(2)}_{R}\times \mathrm{USp}(4)\times \mathrm{USp}(12)$, respectively. We tabulate the full spectrum of chiral particles in Model 20 in table 6.
Table 6. Spectrum of chiral particles of Model 20.
Model 20 $\mathrm{SU}{(4)}_{C}\times \mathrm{SU}{(2)}_{L}\times \mathrm{SU}{(2)}_{R}\times \mathrm{USp}{\left(4\right)}^{2}$ Q4C Q2L Q2R Qem BL Field
ab $6\times (4,\overline{2},1,1,1)$ 1 −1 0 $-\tfrac{1}{3},\tfrac{2}{3},-1,0$ $\tfrac{1}{3},-1$ QL, LL
${ab}^{\prime} $ $2\times (\overline{4},\overline{2},1,1,1)$ −1 −1 0 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QL, LL
ac $4\times (\overline{4},1,2,1,1)$ −1 0 1 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QR, LR
bc $4\times (1,\overline{2},2,1,1)$ 0 −1 1 0, ±1 0 $H^{\prime} $
c2 $1\times (1,1,\overline{2},4,1)$ 0 0 −1 $\mp \tfrac{1}{2}$ 0
a4 $2\times (\overline{4},1,1,1,4)$ −1 0 0 $-\tfrac{1}{6},\tfrac{1}{2}$ $-\tfrac{1}{3},1$
b4 $3\times (1,2,1,1,\overline{4})$ 0 1 0 $\pm \tfrac{1}{2}$ 0
c4 $4\times (1,1,2,1,\overline{4})$ 0 0 1 $\pm \tfrac{1}{2}$ 0
a□□ 1 × (10, 1, 1, 1, 1) 2 0 0 $\tfrac{1}{3},-\tfrac{1}{3},-1$ $\tfrac{2}{3},-2$
$1\times (\overline{6},1,1,1,1)$ −2 0 0 $-\tfrac{1}{3},\tfrac{1}{3},1$ $-\tfrac{2}{3},2$
b□□ 1 × (1, 3, 1, 1, 1) 0 2 0 0, ±1 0
1 × (1, 1, 1, 1, 1) 0 2 0 0 0
c□□ 5 × (1, 1, 3, 1, 1) 0 0 2 0, ±1 0
27 × (1, 1, 1, 1, 1) 0 0 2 0 0
${bc}^{\prime} $ $7\times (1,\overline{2},\overline{2},1,1)$ 0 −1 −1 0, ±1 0 $H^{\prime} $
7 × (1, 2, 2, 1, 1) 0 1 1
Models 2225 are built with at least three USp groups. Their gauge groups are ${\rm{U}}(4)\times {\rm{U}}{(2)}_{L}\times {{\rm{U}}(2)}_{R}\times \mathrm{USp}{\left(2\right)}^{2}\,\times \mathrm{USp}(4)$, ${\rm{U}}{(4)\times {\rm{U}}{(2)}_{L}\times {\rm{U}}(2)}_{R}\times \mathrm{USp}{\left(2\right)}^{3}\times \mathrm{USp}(8)$, ${\rm{U}}(4)\,\times {{\rm{U}}{(2)}_{L}\times {\rm{U}}(2)}_{R}\times \mathrm{USp}{\left(2\right)}^{3}\times \mathrm{USp}(4)$ and ${\rm{U}}(4)\times {\rm{U}}{(2)}_{L}\,\times {{\rm{U}}(2)}_{R}\,\times \mathrm{USp}{\left(2\right)}^{3}\times \mathrm{USp}(4)$, respectively. Note that Model 24 and Model 25 are related by T-duality, but are not phenomenologically equivalent. This can be easily seen from the fact that Model 24 has U(4) and U(2)R gauge coupling unification at the string scale, while Model 25 has U(4) and U(2)L gauge coupling unification. We represent the full spectrum of chiral particles in Model 24 in table 7.
Table 7. Spectrum of chiral particles of Model 24.
Model 24 $\mathrm{SU}{(4)}_{C}\times \mathrm{SU}{(2)}_{L}\times \mathrm{SU}{(2)}_{R}\times \mathrm{USp}{\left(2\right)}^{3}\times \mathrm{USp}(4)$ Q4C Q2L Q2R Qem BL Field
${ab}^{\prime} $ 4 × (4, 2, 1, 1, 1, 1, 1) 1 1 0 $-\tfrac{1}{3},\tfrac{2}{3},-1,0$ $\tfrac{1}{3},-1$ QL, LL
ac $4\times (\overline{4},1,2,1,1,1,1)$ −1 0 1 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QR, LR
bc $2\times (1,2,\overline{2},1,1,1,1)$ 0 1 −1 0, ±1 0 $H^{\prime} $
${bc}^{\prime} $ $4\times (1,\overline{2},\overline{2},1,1,1,1)$ 0 −1 −1 0, ±1 0 H
a3 $1\times (\overline{4},1,1,1,1,2,1)$ −1 0 0 $-\tfrac{1}{6},\tfrac{1}{2}$ $-\tfrac{1}{3},1$
a4 $2\times (4,1,1,1,1,1,\overline{4})$ 1 0 0 $\tfrac{1}{6},-\tfrac{1}{2}$ $\tfrac{1}{3},-1$
b1 $1\times (1,\overline{2},1,2,1,1,1)$ 0 −1 0 $\mp \tfrac{1}{2}$ 0
b2 $3\times (1,2,1,1,\overline{2},1,1)$ 0 1 0 $\pm \tfrac{1}{2}$ 0
c2 $2\times (1,1,2,1,\overline{2},1,1)$ 0 0 1 $\pm \tfrac{1}{2}$ 0
c4 $1\times (1,1,\overline{2},1,1,1,4)$ 0 0 −1 $\mp \tfrac{1}{2}$ 0
${b}_{\overline{\square \square }}$ $2\times (1,\overline{3},1,1,1,1,1)$ 0 −2 0 0, ±1 0
2 × (1, 1, 1, 1, 1, 1, 1) 0 2 0 0 0
${c}_{\overline{\square \square }}$ $1\times (1,1,\overline{3},1,1,1,1)$ 0 0 −2 0, ±1 0
1 × (1, 1, 1, 1, 1, 1, 1) 0 0 2 0 0
The exotic particles charged by USp groups may form bound states and composite particles at some intermediate energy scale, as the strong coupling dynamics of the USp groups require. The composite particles are consistent with anomaly cancellation conditions, as in the QCD case. These composite particles thus are charged only under the SM gauge symmetry [34]. There are essentially two kinds of neutral bound states. The first one comes from decomposing the rank 2 anti-symmetric representation of the USp groups into two fundamental representations and then taking the pseudo inner product of the fundamental representations. The second one comes from the rank 2N anti-symmetric representation of USp(2N) for N ≥ 2. The first bound state is similar to a meson that is the inner product of a fundamental representation and an anti-fundamental representation of SU(3)C in QCD. The second bound state is a USp(2N) singlet, which is an analog to a baryon being a rank 3 anti-symmetric representation of SU(3)C. Models 14, 20 and 25 contain the second kind of bound state. Now we take Models 14, 20 and 25 as examples to show explicitly how bound states are formed.
The composite particle spectrum for Model 14 is listed in table 8. The confining group is USp(4), with two charged intersections. The mixing of intersection a4 and c4 results in the chiral supermultiplets (4, 1, 2, 1, 1).
Table 8. The composite particle spectrum for Model 14.
Model 14 SU(4)C × SU(2)L × SU(2)R × USp(2) × USp(4)
Confining force Intersection Exotic particle spectrum Confined particle spectrum
$\mathrm{USp}{\left(4\right)}_{4}$ a4 $2\times (4,1,1,1,\overline{4})$ 3 × (6, 1, 1, 1, 1), 3 × (10, 1, 1, 1, 1), 4 × (4, 1, 2, 1, 1)
c4 2 × (1, 1, 2, 1, 4) 3 × (1, 1, 1, 1, 1), 3 × (1, 1, 3, 1, 1)
$\mathrm{USp}{\left(2\right)}_{1}$ b1 4 × (1, 2, 1, 2, 1) 10 × (1, 1, 1, 1, 1), 10 × (1, 1, 3, 1, 1)
Moreover, the confined particle spectrum for Model 20 in table 9. The confining group is USp(4), and has three charged intersections. Besides self-confinement, it is also viable to form mixed confinement between sections within the same confining group. The chiral supermultiplets $(\overline{4},2,1,1,1),(1,2,2,1,1)$ and $(\overline{4},1,2,1,1)$ are yielded by the mixed-confinement between intersections a4, b4 and c4.
Table 9. The composite particle spectrum for Model 20.
Model 20 $\mathrm{SU}{(4)}_{C}\times \mathrm{SU}{(2)}_{L}\times \mathrm{SU}{(2)}_{R}\times \mathrm{USp}{\left(4\right)}^{2}$
Confining force Intersection Exotic particle spectrum Confined particle spectrum
$\mathrm{USp}{\left(4\right)}_{4}$ a4 $2\times (\overline{4},1,1,1,4)$ $3\times (\overline{6},1,1,1,1),3\times (\overline{10},1,1,1,1),6\times (\overline{4},2,1,1,1)$
b4 $3\times (1,2,1,1,\overline{4})$ 6 × (1, 1, 1, 1, 1), 6 × (1, 3, 1, 1, 1), 12 × (1, 2, 2, 1, 1)
c4 $4\times (1,1,2,1,\overline{4})$ $10\times (1,1,1,1,1),10\times (1,1,3,1,1),8\times (\overline{4},1,2,1,1)$
$\mathrm{USp}{\left(4\right)}_{2}$ c2 $1\times (1,1,\overline{2},4,1)$ $1\times (1,1,\overline{3},1,1),1\times (1,1,1,1,1)$
For Model 25, the composite particle spectrum is given in table 10. There are two confining groups USp(4) and USp(2), each with two charged intersections. The mixed-confinement between intersections b2, c2 yields the chiral supermultiplet (1, 2, 2, 1, 1, 1, 1), while the mixed-confinement between intersections a4, c4 yields the chiral supermultiplet (4, 1, 2, 1, 1, 1, 1). Note that when there is only one charged intersection, mixed-confinement will not be formed, thus only the tensor representations from self-confinement are left. Checking from the composite particle spectra, one finds that no new anomaly is introduced to the remaining gauge symmetry. Thus our models are free of anomalies. The above analysis for composite particles applies to all our models except Model 15 and Model 16 without any confining group.
Table 10. The composite particle spectrum for Model 25.
Model 25 $\mathrm{SU}{(4)}_{C}\times \mathrm{SU}{(2)}_{L}\times \mathrm{SU}{(2)}_{R}\times \mathrm{USp}{\left(2\right)}^{3}\times \mathrm{USp}(4)$
Confining force Intersection Exotic particle spectrum Confined particle spectrum
$\mathrm{USp}{\left(4\right)}_{2}$ b2 $3\times (1,2,1,1,\overline{4},1,1)$ 6 × (1, 1, 1, 1, 1, 1, 1), 6 × (1, 3, 1, 1, 1, 1, 1), 6 × (1, 2, 2, 1, 1, 1, 1)
c2 $2\times (1,1,2,1,\overline{4},1,1)$ 3 × (1, 1, 1, 1, 1, 1, 1, 1), 3 × (1, 1, 3, 1, 1, 1, 1)
$\mathrm{USp}{\left(2\right)}_{4}$ a4 $2\times (4,1,1,1,1,1,\overline{2})$ 3 × (6, 1, 1, 1, 1, 1, 1), 3 × (10, 1, 1, 1, 1, 1, 1), 2 × (4, 1, 2, 1, 1, 1, 1)
c4 1 × (1, 1, 2, 1, 1, 1, 2) 1 × (1, 1, 1, 1, 1, 1, 1), 1 × (1, 1, 3, 1, 1, 1, 1)
$\mathrm{USp}{\left(2\right)}_{1}$ b1 $1\times (1,2,1,\overline{2},1,1,1)$ 1 × (1, 1, 1, 1, 1, 1, 1), 1 × (1, 3, 1, 1, 1, 1, 1)
$\mathrm{USp}{\left(2\right)}_{3}$ a3 1 × (4, 1, 1, 1, 1, 2, 1) 1 × (10, 1, 1, 1, 1, 1, 1), 1 × (6, 1, 1, 1, 1, 1, 1)

4.2. Models with one tilted torus

In this section, we show basic phenomenological properties of models with one tilted torus. Model 36 represents the models with exact gauge coupling unification at the string level so far. The gauge symmetries therein are ${\rm{U}}{(4)\times {\rm{U}}{(2)}_{L}\times {\rm{U}}(2)}_{R}\times \mathrm{USp}{\left(2\right)}^{2}$ with two confining groups in the hidden sector. The full spectrum of this model is shown in table 11.
Table 11. Spectrum of chiral particles of Model 36.
Model 36 $\mathrm{SU}{(4)}_{C}\times \mathrm{SU}{(2)}_{L}\times \mathrm{SU}{(2)}_{R}\times \mathrm{USp}{\left(2\right)}^{2}$ Q4C Q2L Q2R Qem BL Field
${ab}^{\prime} $ 4 × (4, 2, 1, 1, 1) 1 1 0 $-\tfrac{1}{3},\tfrac{2}{3},-1,0$ $\tfrac{1}{3},-1$ QL, LL
${ac}^{\prime} $ $4\times (\overline{4},1,\overline{2},1,1)$ −1 0 −1 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QR, LR
b4 $4\times (1,2,1,1,\overline{2})$ 0 1 0 $\pm \tfrac{1}{2}$ 0
c2 $4\times (1,1,2,\overline{2},1)$ 0 0 1 $\pm \tfrac{1}{2}$ 0
${b}_{\overline{\square \square }}$ $3\times (1,\overline{3},1,1,1)$ 0 −2 0 0, ±1 0
2 × (1, 1, 1, 1, 1) 0 2 0 0 0
c□□ 3 × (1, 1, 3, 1, 1) 0 0 2 0, ±1 0
3 × (1, 1, 1, 1, 1) 0 0 2 0 0
bc $8\times (1,2,\overline{2},1,1)$ 0 1 −1 0, ±1 0 ${H}_{u}^{i},{H}_{d}^{i}$
$8\times (1,\overline{2},2,1,1)$ 0 −1 1
Model 37 and Model 38 have no USp group. The gauge group for the two models is U(4)C × U(2)L × U(2)R. Since there is no USp group, the gaugino condensation mechanism will not work. Thus, one needs to find other mechanisms for supersymmetry breaking. Also, there are no exotic particles in these two models, as exotic particles are charged under USp groups. We show the full chiral spectrum in the open string sectors for Model 37 in table 12.
Table 12. Spectrum of chiral particles of Model 37.
Model 37 SU(4)C × SU(4)L × SU(4)R Q4C Q2L Q2R Qem BL Field
${ab}^{\prime} $ $4\times (\overline{4},\overline{4},1)$ −1 −1 0 $\tfrac{1}{3},-\tfrac{2}{3},1,0$ $-\tfrac{1}{3},1$ QL, LL
${ac}^{\prime} $ 4 × (4, 1, 4) 1 0 1 $-\tfrac{1}{3},\tfrac{2}{3},-1,0$ $\tfrac{1}{3},-1$ QR, LR
${b}_{\overline{\square \square }}$ $2\times (1,\overline{10},1)$ 0 −2 0 0, ±1 0
2 × (1, 6, 1) 0 2 0 0 0
c□□ 2 × (1, 1, 10) 0 0 2 0, ±1 0
2 × (1, 1, 6) 0 0 2 0 0
bc $4\times (1,\overline{4},4)$ 0 −1 1 0, ±1 0 ${H}_{u}^{i},{H}_{d}^{i}$
$4\times (1,4,\overline{4})$ 0 1 −1
Models 39, 40, 41, 42 have only one USp(2) group in the hidden sector. Their gauge symmetry is all ${\rm{U}}{(4)\times {\rm{U}}{(4)}_{L}\times {\rm{U}}(2)}_{R}\times \mathrm{USp}(2)$. We note that Model 39 and Model 40 are T-dual to each other, as well as Model 41 and Model 42. But they are clearly not equivalent models at the phenomenological level due to b and c stacks of brane swapping. In Model 39, we have SU(3)C and U(1)Y gauge coupling unification at the string level, while in Model 40 this gauge unification gets swapped to SU(3)C and SU(2)L gauge coupling unification at the string level. Similarly, this b and c stacks of brane swapping appear between Model 41 and Model 42 as well. The whole spectrum of chiral particles of Model 39 in table 13.
Table 13. Spectrum of chiral particles of Model 39.
Model 39 SU(4)C × SU(4)L × SU(2)R × USp(2) Q4C Q4L Q2R Qem BL Field
ab $4\times (4,\overline{4},1,1)$ 1 −1 0 $-\tfrac{1}{3},\tfrac{2}{3},-1,0$ $\tfrac{1}{3},-1$ QL, LL
ac $4\times (\overline{4},1,2,1)$ −1 0 1 $-\tfrac{1}{3},\tfrac{2}{3},-1,0$ $\tfrac{1}{3},-1$ QR, LR
${bc}^{\prime} $ 4 × (1, 4, 2, 1) 0 1 1 0, ±1 0 $H^{\prime} $
b1 $4\times (1,2,1,\overline{2})$ 0 1 0 $\pm \tfrac{1}{2}$ 0
b□□ 2 × (1, 10, 1, 1) 0 2 0 $\tfrac{1}{3},-\tfrac{1}{3},-1$ $\tfrac{2}{3},-2$
$2\times (1,\overline{6},1,1)$ 0 −2 0 $-\tfrac{1}{3},1$ $-\tfrac{2}{3},2$
${c}_{\overline{\square \square }}$ $3\times (1,1,\overline{3},1)$ 0 0 2 0, ±1 0
3 × (1, 1, 1, 1) 0 0 2 0 0
bc $6\times (1,2,\overline{2},1)$ 0 1 −1 0, ±1 0 ${H}_{u}^{i},{H}_{d}^{i}$
$6\times (1,\overline{2},2,1)$ 0 −1 1
For the models with one tilted torus, the models represented by Model 36 are the only class with gauge coupling unification and carry two confining USp groups. Thus gaugino condensation can trigger supersymmetry breaking and moduli stabilization [36].

5. Discussions and conclusions

Utilizing the deterministic algorithm, we obtain various classes of four-family supersymmetric models from intersecting D6-branes on ${{\mathbb{T}}}^{6}/{{\mathbb{Z}}}_{2}\times {{\mathbb{Z}}}_{2}$ orientifold, with and without tilted torus. In total, there are 274 physical independent four-family supersymmetric models without tilted torus, and 6 physical independent four-family supersymmetric models with the third torus tilted, without loss of generality.
For models without tilted torus, Model 14 represents the model with gauge coupling unification at the string scale. Models 15 and 16 are the rare models without any USp group in the hidden sectors, with tadpole cancellation conditions satisfied. Models 2226 are with at least two confining USp groups. Thus gaugino condensation can be triggered to break the supersymmetry and stabilize the moduli. Moreover, there are Models 2735 with wrapping numbers absolutely larger than 5 which was not reached for three-family supersymmetric Pati–Salam models as discussed in [1].
For models with one tilted torus, Models 37 and 38 satisfy the tadpole cancellation conditions without any filler branes. Models such as Models 39 and 40 are related with b and c stacks of branes swapping. Model 36 represents the models with exact gauge coupling unification at the string scale, with two confining USp groups allowing gaugino condensation as well. This class of models would be ideal for further phenomenology model buildings (such as in [61]) as gaugino condensation can be triggered to break the supersymmetry and stabilize the moduli, while gauge coupling unified at string scale.

Appendix A. Four-family standard models from intersecting D6-branes without tilted toriAppendix B.Four-family standard models from intersecting D6-branes with one tilted torus

In the appendix, we list all representative four-family models obtained from our random scanning method. In the first columns for each table, a, b, c represent three stacks of D6-branes, respectively. Also in the first columns, 1, 2, 3, 4 is a short-handed notation for the filler branes along the ΩR, ΩRω, ΩRθω and ΩRθ O6-planes, respectively. The second columns for each table list the numbers of D6-branes in every stack, respectively. In the third column of each table, we record wrapping numbers of each D6-brane configuration and designate the third ${{\mathbb{T}}}^{2}$ to be tilted.
Table 14. D6-brane configurations and intersection numbers of Model 14, and its gauge coupling relation is ${g}_{a}^{2}={g}_{b}^{2}={g}_{c}^{2}=\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{8}{3}\sqrt[4]{2}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 14 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(2)\times {USp}(4)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 1 4
a 8 (−1, 0) × (−1, 1) × (1, 2) −1 1 4 0 −4 0 0 2
b 4 (0, 1) × (−1, 2) × (−1, 2) 3 −3 0 −8 −4 0
c 4 (1, 1) × (−1, 0) × (−1, 2) 1 −1 0 −2
1 2 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=\tfrac{1}{4}{x}_{B}=\tfrac{1}{4}{x}_{C}=\tfrac{1}{2}{x}_{D}$
4 4 (0, 1) × (0, 1) × (1, 0) ${\beta }_{1}^{g}=-2$, ${\beta }_{4}^{g}=0$
${\chi }_{1}=\tfrac{1}{\sqrt{2}}$, ${\chi }_{2}=\tfrac{1}{\sqrt{2}}$, ${\chi }_{3}=\tfrac{1}{\sqrt{2}}$
Table 15. D6-brane configurations and intersection numbers of Model 15, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{7}{9}{g}_{b}^{2}=\tfrac{7}{3}{g}_{c}^{2}=\tfrac{35}{23}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{4\ {5}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}}{3\sqrt{3}}$.
Model 15 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $
a 8 (−1, 0) × (−1, 1) × (1, 1) 0 0 8 −4 −4 0
b 4 (−1, 2) × (0, 1) × (−1, 3) 5 −5 −16 0
c 4 (1, 2) × (−2, 1) × (−1, 1) −7 −9
${x}_{A}=\tfrac{1}{5}{x}_{B}=\tfrac{1}{6}{x}_{C}=\tfrac{1}{6}{x}_{D}$
${\chi }_{1}=\tfrac{\sqrt{5}}{6}$, ${\chi }_{2}=\tfrac{1}{\sqrt{5}}$, ${\chi }_{3}=\tfrac{2}{\sqrt{5}}$
Table 16. D6-brane configurations and intersection numbers of Model 16, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{10}{9}{g}_{b}^{2}=\tfrac{10}{3}{g}_{c}^{2}=\tfrac{50}{29}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{2}{3}\sqrt{\tfrac{2}{3}}{11}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 16 $U(4)\times U{\left(4\right)}_{L}\times U{\left(4\right)}_{R}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $
a 8 (−1, 0) × (1, 1) × (−1, 1) 0 0 8 −4 −4 0
b 8 (−1, 1) × (−1, 3) × (0, 1) 4 −4 0 −8
c 8 (−1, 1) × (1, 2) × (−1, −1) 8 24
${x}_{A}=\tfrac{1}{11}{x}_{B}=\tfrac{1}{3}{x}_{C}=\tfrac{1}{3}{x}_{D}$
${\chi }_{1}=\tfrac{\sqrt{11}}{3}$, ${\chi }_{2}=\tfrac{1}{\sqrt{11}}$, ${\chi }_{3}=\tfrac{2}{\sqrt{11}}$
Table 17. D6-brane configurations and intersection numbers of Model 17, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{5}{6}{g}_{b}^{2}=\tfrac{5}{2}{g}_{c}^{2}=\tfrac{25}{16}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{4}{9}{14}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 17 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(2)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 2
a 8 (−1, 0) × (−1, 1) × (1, 1) 0 0 6 −2 −4 0 0
b 4 (−1, 2)× (0, 1) × (−1, 2) 3 −3 −8 0 0
c 4 (1, 2) × (−2, 1) × (−1, 1) −7 −9 4
2 2 (1, 0) × (0, 1) × (0, 1) ${x}_{A}=\tfrac{2}{7}{x}_{B}=\tfrac{1}{4}{x}_{C}=\tfrac{1}{4}{x}_{D}$
${\beta }_{2}^{g}=-2$
${\chi }_{1}=\tfrac{\sqrt{\tfrac{7}{2}}}{4}$, ${\chi }_{2}=\sqrt{\tfrac{2}{7}}$, ${\chi }_{3}=2\sqrt{\tfrac{2}{7}}$
Table 18. D6-brane configurations and intersection numbers of Model 18, and its gauge coupling relation is ${g}_{a}^{2}={g}_{b}^{2}=3{g}_{c}^{2}=\tfrac{5}{3}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{8}{3}{2}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 18 $U(4)\times U{\left(4\right)}_{L}\times U{\left(4\right)}_{R}\times {USp}(4)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 3
a 8 (−1, 0) × (−1, 2) × (1, 1) 1 −1 4 0 −4 0 −2
b 8 (−1, 1) × (0, 1) × (−1, 1) 0 0 0 0 2
c 8 (1, 1) × (−1, 1) × (−1, 1) −4 −12 −2
3 4 (0, 1) × (1, 0) × (0, 1) ${x}_{A}=\tfrac{1}{4}{x}_{B}={x}_{C}=\tfrac{1}{2}{x}_{D}$
${\beta }_{3}^{g}=2$
${\chi }_{1}=\sqrt{2}$, ${\chi }_{2}=\tfrac{1}{2\sqrt{2}}$, ${\chi }_{3}=\sqrt{2}$
Table 19. D6-brane configurations and intersection numbers of Model 19, and its gauge coupling relation is ${g}_{a}^{2}=2{g}_{b}^{2}={g}_{c}^{2}=\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{8}{3}{2}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 19 $U(4)\times U{\left(4\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(2)\times {USp}(4)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 1 4
a 8 (−1, 0) × (−1, 1) × (1, 2) −1 1 0 4 −4 0 0 2
b 8 (0, 1) × (−1, 1)×(−1, 1) 0 0 2 −6 −2 0
c 4 (1, 1) × (−1, 0) × (−1, 2) 1 −1 0 −2
1 2 (1, 0) × (1, 0) × (1, 0) xA = xB = xC = 2xD
4 4 (0, 1) × (0, 1) × (1, 0) ${\beta }_{1}^{g}=-4$, ${\beta }_{4}^{g}=0$
${\chi }_{1}=\sqrt{2}$, ${\chi }_{2}=\sqrt{2}$, ${\chi }_{3}=\sqrt{2}$
Table 20. D6-brane configurations and intersection numbers of Model 20, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{21}{22}{g}_{b}^{2}=\tfrac{7}{2}{g}_{c}^{2}=\tfrac{7}{4}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{4}{11}\sqrt{3}{10}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 20 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}{\left(4\right)}^{2}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 2 4
a 8 (−1, 0) × (1, 1) × (−1, 2) 1 −1 6 −2 −4 0 0 −2
b 4 (−3, 2) × (−1, 2) × (0, 1) 1 −1 −4 0 0 3
c 4 (−2, 1) × (1, 2) × (−1, −2) 5 27 −1 4
2 4 (1, 0) × (0, 1) × (0, 1) ${x}_{A}=\tfrac{1}{20}{x}_{B}=\tfrac{3}{8}{x}_{C}=\tfrac{3}{4}{x}_{D}$
4 4 (0, 1) × (0, 1) × (1, 0) ${\beta }_{2}^{g}=-5$, ${\beta }_{4}^{g}=5$
${\chi }_{1}=\tfrac{3\sqrt{\tfrac{5}{2}}}{2}$, ${\chi }_{2}=\tfrac{1}{\sqrt{10}}$, ${\chi }_{3}=\tfrac{1}{\sqrt{10}}$
Table 21. D6-brane configurations and intersection numbers of Model 21, and its gauge coupling relation is ${g}_{a}^{2}={g}_{b}^{2}=4{g}_{c}^{2}=\tfrac{20}{11}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{32}{9}\sqrt[4]{2}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 21 $U(4)\times U{\left(4\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(4)\times {USp}(12)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 2 4
a 8 (−1, 0) × (1, 1) × (−1, 2) 1 −1 4 0 −4 0 0 −2
b 8 (−1, 1) × (−1, 1) × (0, 1) 0 0 −6 6 0 2
c 4 (−2, 1) × (1, 2) × (−1, −2) 5 27 −1 4
2 12 (1, 0) × (0, 1) × (0, 1) ${x}_{A}=\tfrac{1}{16}{x}_{B}=\tfrac{1}{2}{x}_{C}={x}_{D}$
4 4 (0, 1) × (0, 1) × (1, 0) ${\beta }_{2}^{g}=-5$, ${\beta }_{4}^{g}=4$
${\chi }_{1}=2\sqrt{2}$, ${\chi }_{2}=\tfrac{1}{2\sqrt{2}}$, ${\chi }_{3}=\tfrac{1}{2\sqrt{2}}$
Table 22. D6-brane configurations and intersection numbers of Model 22, and its gauge coupling relation is ${g}_{a}^{2}={g}_{b}^{2}=2{g}_{c}^{2}=\tfrac{10}{7}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=2\sqrt{2}\sqrt[4]{3}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 22 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}{\left(2\right)}^{2}\times {USp}(4)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 2 3 4
a 8 (−1, 0) × (−1, 1) × (1, 1) 0 0 4 0 −4 0 0 −1 1
b 4 (0, 1) × ( − 1, 3) × ( − 1, 1) 2 −2 0 −12 1 0 0
c 4 (1, 1) × ( − 1, 0) × ( − 2, 2) 0 0 2 0 −2
2 4 (1, 0) × (0, 1) × (0, 1) ${x}_{A}=\tfrac{1}{3}{x}_{B}=\tfrac{1}{3}{x}_{C}=\tfrac{1}{3}{x}_{D}$
3 2 (0, 1) × (1, 0) × (0, 1) ${\beta }_{2}^{g}=-3$, ${\beta }_{3}^{g}=-4$, ${\beta }_{4}^{g}=-2$
4 2 (0, 1) × (0, 1) × (1, 0) ${\chi }_{1}=\tfrac{1}{\sqrt{3}}$, ${\chi }_{2}=\tfrac{1}{\sqrt{3}}$, ${\chi }_{3}=\tfrac{2}{\sqrt{3}}$
Table 23. D6-brane configurations and intersection numbers of Model 23, and its gauge coupling relation is ${g}_{a}^{2}={g}_{b}^{2}=\tfrac{5}{2}{g}_{c}^{2}=\tfrac{25}{16}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{2}{3}{11}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 23 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}{\left(2\right)}^{3}\times {USp}(8)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 1 2 3 4
a 8 ( − 1, 0) × (1, 1) × ( − 1, 1) 0 0 4 0 −4 0 0 0 1 −1
b 4 ( − 1, 2) × ( − 1, 1) × (0, 1) 1 −1 −4 6 −2 0 0 1
c 4 ( − 1, 1) × (1, 3) × ( − 1, − 1) 0 12 3 −1 3 1
1 2 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=\tfrac{1}{11}{x}_{B}=\tfrac{1}{2}{x}_{C}=\tfrac{1}{2}{x}_{D}$
2 8 (1, 0) × (0, 1) × (0, 1) ${\beta }_{1}^{g}=-1$, ${\beta }_{2}^{g}=-5$, ${\beta }_{3}^{g}=-1$, ${\beta }_{4}^{g}=-2$
3 2 (0, 1) × (1, 0) × (0, 1) ${\chi }_{1}=\tfrac{\sqrt{11}}{2}$, ${\chi }_{2}=\tfrac{1}{\sqrt{11}}$, ${\chi }_{3}=\tfrac{2}{\sqrt{11}}$
4 2 (0, 1) × (0, 1) × (1, 0)
Table 24. D6-brane configurations and intersection numbers of Model 24, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{3}{2}{g}_{b}^{2}={g}_{c}^{2}=\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=2\ {3}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 24 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}{\left(2\right)}^{3}\times {USp}(4)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 1 2 3 4
a 8 (−1, 0) × (−1, 1) × (1, 1) 0 0 0 4 −4 0 0 0 −1 1
b 4 (0, 1) × (−1, 1) × (−3, 1) −2 2 2 −4 −1 3 0 0
c 4 (1, 2) × (−1, 0) × (−1, 1) −1 1 0 2 0 −1
1 2 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=3{x}_{B}=\tfrac{3}{2}{x}_{C}=\tfrac{3}{2}{x}_{D}$
2 4 (1, 0) × (0, 1) × (0, 1) ${\beta }_{1}^{g}=-5$, ${\beta }_{2}^{g}=-1$, ${\beta }_{3}^{g}=-4$, ${\beta }_{4}^{g}=-3$
3 2 (0, 1) × (1, 0) × (0, 1) ${\chi }_{1}=\tfrac{\sqrt{3}}{2}$, ${\chi }_{2}=\sqrt{3}$, ${\chi }_{3}=2\sqrt{3}$
4 2 (0, 1) × (0, 1) × (1, 0)
Table 25. D6-brane configurations and intersection numbers of Model 25, and its gauge coupling relation is ${g}_{a}^{2}={g}_{b}^{2}=\tfrac{3}{2}{g}_{c}^{2}=\tfrac{5}{4}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=2\ {3}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 25 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}{\left(2\right)}^{3}\times {USp}(4)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b $b^{\prime} $ c $c^{\prime} $ 1 2 3 4
a 8 (1, 0) × (1, 1) × (1, −1) 0 0 4 0 −4 0 0 0 1 −1
b 4 (−1, 2) × (−1, 1) × (0, 1) 1 −1 4 2 −2 0 0 1
c 4 (0, 1) × (1, 3) × (−1, −1) −2 2 3 −1 0 0
1 4 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=\tfrac{1}{3}{x}_{B}=\tfrac{1}{2}{x}_{C}=\tfrac{1}{2}{x}_{D}$
2 2 (1, 0) × (0, 1) × (0, 1) ${\beta }_{1}^{g}=-1$, ${\beta }_{2}^{g}=-5$, ${\beta }_{3}^{g}=-4$, ${\beta }_{4}^{g}=-3$
3 2 (0, 1) × (1, 0) × (0, 1) ${\chi }_{1}=\tfrac{\sqrt{3}}{2}$, ${\chi }_{2}=\tfrac{1}{\sqrt{3}}$, ${\chi }_{3}=\tfrac{2}{\sqrt{3}}$
4 2 (0, 1) × (0, 1) × (1, 0)
Table 26. D6-brane configurations and intersection numbers of Model 26, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{11}{10}{g}_{b}^{2}=\tfrac{11}{5}{g}_{c}^{2}=\tfrac{55}{37}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{8\ {13}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}}{7\sqrt{5}}$.
Model 26 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}{\left(2\right)}^{3}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 2 3 4
a 8 (1, −1) × (1, 0) × (1, 1) 0 0 −3 7 −4 0 −1 0 1
b 4 (−5, 2) × (1, 1) × (−1,0) −3 3 −2 0 −2 5 0
c 4 (−3, 1) × (−1, 1) × (−1, 1) 0 12 1 3 3
2 2 (1, 0) × (0, 1) × (0, 1) ${x}_{A}=\tfrac{26}{5}{x}_{B}=13{x}_{C}=\tfrac{26}{5}{x}_{D}$
3 2 (0, 1) × (1, 0) × (0, 1) ${\beta }_{2}^{g}$ = −1, ${\beta }_{3}^{g}=2$, βg4 = −1
4 2 (0, 1) × (0, 1) × (1, 0) ${\chi }_{1}=\sqrt{13}$, ${\chi }_{2}=\tfrac{2\sqrt{13}}{5}$, ${\chi }_{3}=2\sqrt{13}$
Table 27. D6-brane configurations and intersection numbers of Model 27, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{6}{7}{g}_{b}^{2}=\tfrac{2}{7}{g}_{c}^{2}=\tfrac{2}{5}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{6}{7}\sqrt{2}\sqrt[4]{3}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 27 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c
a 8 (−2, −1) × (1, 1) × (1, 1) 0 −8 12 −8 −4 0
b 4 (−1, 1) × (6, 2) × (−1,0) 4 −4 8 0
c 4 (1, 1) × (−1,0) × (−2, 2) 0 0
xA = 9xB = 3xC = 9xD
${\chi }_{1}=\sqrt{3}$, ${\chi }_{2}=3\sqrt{3}$, ${\chi }_{3}=2\sqrt{3}$
Table 28. D6-brane configurations and intersection numbers of Model 28, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{47}{112}{g}_{b}^{2}=\tfrac{18}{7}{g}_{c}^{2}=\tfrac{30}{19}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{1}{3}\sqrt{\tfrac{2}{7}}{23}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 28 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(2)\times {USp}(4)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 1 2
a 8 (1, −1) × (1, 0) × (2, 1) 1 −1 −5 9 −4 0 0 −2
b 4 (−7, 2) × (1, 1) × (−1,0) −5 5 9 11 0 −2
c 4 (−2, 1) × (−2, 1) × (−2, 1) 5 27 −1 4
1 2 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=\tfrac{92}{7}{x}_{B}=46{x}_{C}=\tfrac{46}{7}{x}_{D}$
2 4 (1, 0) × (0, 1) × (0, 1) ${\beta }_{1}^{g}=-5$, βg2 = 4
${\chi }_{1}=\sqrt{23}$, ${\chi }_{2}=\tfrac{2\sqrt{23}}{7}$, ${\chi }_{3}=4\sqrt{23}$
Table 29. D6-brane configurations and intersection numbers of Model 29, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{5}{24}{g}_{b}^{2}=\tfrac{19}{8}{g}_{c}^{2}=\tfrac{95}{62}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{2}{45}{86}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 29 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}{\left(2\right)}^{2}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 2 4
a 8 (2, −1) × (1, 0) × (2, 1) 0 0 −6 10 −4 0 −2 2
b 4 (−8, 1) × (1, 1) × (−1,0) −7 7 15 11 −1 0
c 4 (−3, 1) × (−2, 1) × (−2, 1) 9 39 4 6
2 2 (1, 0) × (0, 1) × (0, 1) ${x}_{A}=\tfrac{43}{4}{x}_{B}=86{x}_{C}=\tfrac{43}{4}{x}_{D}$
4 2 (0, 1) × (0, 1) × (1, 0) ${\beta }_{2}^{g}=3$, ${\beta }_{4}^{g}=4$
${\chi }_{1}=\sqrt{86}$, ${\chi }_{2}=\tfrac{\sqrt{\tfrac{43}{2}}}{4}$, ${\chi }_{3}=2\sqrt{86}$
Table 30. D6-brane configurations and intersection numbers of Model 30, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{760}{31}{g}_{b}^{2}=9{g}_{c}^{2}=\tfrac{15}{7}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{48}{31}\sqrt{2}\sqrt[4]{3}{7}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 30 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}{\left(2\right)}^{2}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 1 3
a 8 (−1, 1) × (−1,0) × (1, 1) 0 0 −5 9 −4 0 0 0
b 4 (−2, 1) × (−1, 1) × (−4, 1) 3 29 33 −35 −1 8
c 4 (1, 0) × (−9, −2) × (−1, 1) 7 −7 0 2
1 2 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=42{x}_{B}=\tfrac{28}{3}{x}_{C}=42{x}_{D}$
3 2 (0, 1) × (1, 0) × (0, 1) ${\beta }_{1}^{g}=-5$, ${\beta }_{3}^{g}=4$
${\chi }_{1}=2\sqrt{\tfrac{7}{3}}$, ${\chi }_{2}=3\sqrt{21}$, ${\chi }_{3}=4\sqrt{\tfrac{7}{3}}$
Table 31. D6-brane configurations and intersection numbers of Model 31, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{1309}{47}{g}_{b}^{2}=5{g}_{c}^{2}=\tfrac{25}{13}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{8}{47}\sqrt[4]{2}{185}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 31 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c
a 8 (0, 1) × (−1, −2) × (2, 1) 0 0 10 −6 −4 0
b 4 (1, −2) × (1, −3) × (4, 1) −23 −73 48 −24
c 4 (1, −10) × (0, −1) × (−2, 1) 8 −8
${x}_{A}={x}_{B}=\tfrac{1}{5}{x}_{C}=\tfrac{1}{74}{x}_{D}$
${\chi }_{1}=\tfrac{1}{\sqrt{370}}$, ${\chi }_{2}=\sqrt{\tfrac{5}{74}}$, ${\chi }_{3}=2\sqrt{\tfrac{74}{5}}$
Table 32. D6-brane configurations and intersection numbers of Model 32, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{1248}{41}{g}_{b}^{2}=11{g}_{c}^{2}=\tfrac{11}{5}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{64}{123}{77}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 32 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(14)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 1
a 8 ( − 1, 1) × ( − 1, 0) × (1, 1) 0 0 −5 9 −4 0 0
b 4 (−2, 1) × (−2, 1) × (−4, 1) 13 51 45 −35 −1
c 4 (1, 0) × (−11, −2) × (−1, 1) 9 −9 0
1 14 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=56{x}_{B}=\tfrac{112}{11}{x}_{C}=56{x}_{D}$
${\beta }_{1}^{g}=-5$
${\chi }_{1}=4\sqrt{\tfrac{7}{11}}$, ${\chi }_{2}=2\sqrt{77}$, ${\chi }_{3}=8\sqrt{\tfrac{7}{11}}$
Table 33. D6-brane configurations and intersection numbers of Model 33, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{1680}{47}{g}_{b}^{2}=13{g}_{c}^{2}=\tfrac{65}{29}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{256}{141}\sqrt[4]{2}{13}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 33 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(10)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 1
a 8 (−1, 1) × (−1,0) × (1, 1) 0 0 −5 9 −4 0 0
b 4 (−2, 1) × (−2, 1) × (−4, 1) 13 51 51 −45 −1
c 4 (1, 0) × (−13, −2) × (−1, 1) 11 −11 0
1 10 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=64{x}_{B}=\tfrac{128}{13}{x}_{C}=64{x}_{D}$
βg1 = −5
${\chi }_{1}=8\sqrt{\tfrac{2}{13}}$, ${\chi }_{2}=4\sqrt{26}$, ${\chi }_{3}=16\sqrt{\tfrac{2}{13}}$
Table 34. D6-brane configurations and intersection numbers of Model 34, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{2176}{53}{g}_{b}^{2}=15{g}_{c}^{2}=\tfrac{25}{11}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{192}{53}\sqrt[4]{3}{5}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 34 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(6)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 1
a 8 (−1, 1) × (−1,0) × (1, 1) 0 0 −5 9 −4 0 0
b 4 (−2, 1) × (−2, 1) × (−4, 1) 13 51 57 −55 −1
c 4 (1, 0) × (−15, −2) × (−1, 1) 13 −13 0
1 6 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=72{x}_{B}=\tfrac{48}{5}{x}_{C}=72{x}_{D}$
βg1 = −5
${\chi }_{1}=4\sqrt{\tfrac{3}{5}}$, ${\chi }_{2}=6\sqrt{15}$, ${\chi }_{3}=8\sqrt{\tfrac{3}{5}}$
Table 35. D6-brane configurations and intersection numbers of Model 35, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{2736}{59}{g}_{b}^{2}=17{g}_{c}^{2}=\tfrac{85}{37}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{64}{177}{170}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 35 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(2)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 1
a 8 (−1, 1) × (−1,0) × (1, 1) 0 0 −5 9 −4 0 0
b 4 (−2, 1) × (−2, 1) × (−4, 1) 13 51 63 −65 −1
c 4 (1, 0) × (−17, −2) × (−1, 1) 15 −15 0
1 2 (1, 0) × (1, 0) × (1, 0) ${x}_{A}=80{x}_{B}=\tfrac{160}{17}{x}_{C}=80{x}_{D}$
βg1 = −5
${\chi }_{1}=4\sqrt{\tfrac{10}{17}}$, ${\chi }_{2}=2\sqrt{170}$, ${\chi }_{3}=8\sqrt{\tfrac{10}{17}}$
The rest columns of each table record intersection numbers between stacks. For instance, in the b column of table 1444, from top to bottom, the numbers represent intersection numbers Iab, Ibb, Icb, etc.. As usual, $b^{\prime} $ and $c^{\prime} $ are the orientifold ΩR image of b and c stacks of D6-branes. We also list the relation between xA, xB, xC, xD, which are determined by the supersymmetry condition equation (2.10), as well as the relation between the moduli parameter χ1, χ2, χ3. The one loop β functions ${\beta }_{i}^{g}$ are also listed. To have a clearer sight of gauge couplings, we list them up in the caption of each table, which makes it easier to check whether they are unified.
Table 36. D6-brane configurations and intersection numbers of Model 36, and its gauge coupling relation is ${g}_{a}^{2}={g}_{b}^{2}={g}_{c}^{2}=\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=2\sqrt{2}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 36 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}{\left(2\right)}^{2}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 2 4
a 8 (1, 1) × (0, −1) × (1, 1) 0 0 0 4 0 −4 0 0
b 4 (1, 0) × (−4, 1) × (−1, −1) −3 3 0 0 0 4
c 4 (0, −1) × (4, 1) × (1, 1) 3 −3 4 0
2 2 (1, 0) × (0, 1) × (0, −2) xA = 4xB = xC = 4xD
4 2 (0, 1) × (0, 1) × (−2, 0) βg2 = −2, βg4 = −2
χ1 = 1, χ2 = 4, χ3 = 2
Table 37. D6-brane configurations and intersection numbers of Model 37, and its gauge coupling relation is ${g}_{a}^{2}=2{g}_{b}^{2}=2{g}_{c}^{2}=\tfrac{10}{7}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=4\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 37 $U(4)\times U{\left(4\right)}_{L}\times U{\left(4\right)}_{R}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c
a 8 (1, −1) × (−1,0) × (−1, −1) 0 0 0 4 0 −4
b 8 (0, 1) × (−1, −2) × (1, 1) −2 2 0 0
c 8 (1, 0) × (1, −2) × (1, 1) 2 −2
${x}_{A}=\tfrac{1}{2}{x}_{B}={x}_{C}=\tfrac{1}{2}{x}_{D}$
χ1 = 1, ${\chi }_{2}=\tfrac{1}{2}$, χ3 = 2
Table 38. D6-brane configurations and intersection numbers of Model 38, and its gauge coupling relation is ${g}_{a}^{2}=2{g}_{b}^{2}=2{g}_{c}^{2}=\tfrac{10}{7}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=4\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 38 $U(4)\times U{\left(4\right)}_{L}\times U{\left(4\right)}_{R}$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c
a 8 (1, 0) × (−1, 1) × (−1, −1) 0 0 4 0 0 −4
b 8 (1, −2) × (0, −1) × (−1, 1) 2 −2 0 0
c 8 (1, −2) × (−1,0) × (−1, −1) 2 −2
${x}_{A}={x}_{B}=\tfrac{1}{2}{x}_{C}=\tfrac{1}{2}{x}_{D}$
${\chi }_{1}=\tfrac{1}{2}$, χ2 = 1, χ3 = 2
Table 39. D6-brane configurations and intersection numbers of Model 39, and its gauge coupling relation is ${g}_{a}^{2}=2{g}_{b}^{2}={g}_{c}^{2}=\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{8}{3}\sqrt[4]{2}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 39 $U(4)\times U{\left(4\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(2)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 1
a 8 (1, −1) × (1, 0) × (1, 1) 0 0 4 0 −4 0 0
b 8 (0, 1) × (−1, 2) × (−1, 1) 2 −2 0 4 4
c 4 (−1,0) × (−1, −4) × (1, −1) −3 3 0
1 2 (1, 0) × (1, 0) × (−2, 0) ${x}_{A}=\tfrac{1}{2}{x}_{B}=2{x}_{C}=\tfrac{1}{2}{x}_{D}$
βg1 = −2
${\chi }_{1}=\sqrt{2}$, ${\chi }_{2}=\tfrac{1}{2\sqrt{2}}$, ${\chi }_{3}=2\sqrt{2}$
Table 40. D6-brane configurations and intersection numbers of Model 40, and its gauge coupling relation is ${g}_{a}^{2}={g}_{b}^{2}=2{g}_{c}^{2}=\tfrac{10}{7}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{8}{3}\sqrt[4]{2}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 40 $U(4)\times U{\left(2\right)}_{L}\times U{\left(4\right)}_{R}\times {USp}(2)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 2
a 8 (−1, −1) × (0, −1) × (−1, −1) 0 0 4 0 −4 0 0
b 4 (−1,0) × (4, 1) × (−1, 1) 3 −3 0 −4 0
c 8 (0, −1) × (2, −1) × (−1, 1) −2 2 4
2 2 (1, 0) × (0, 1) × (0, −2) xA = 2xB = xC = 4xD
βg2 = −2
${\chi }_{1}=\sqrt{2}$, ${\chi }_{2}=2\sqrt{2}$, ${\chi }_{3}=\sqrt{2}$
Table 41. D6-brane configurations and intersection numbers of Model 41, and its gauge coupling relation is ${g}_{a}^{2}=2{g}_{b}^{2}={g}_{c}^{2}=\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{8}{3}\sqrt[4]{2}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 41 $U(4)\times U{\left(4\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(2)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 2
a 8 (−1,0) × (1, 1) × (−1, 1) 0 0 0 4 0 −4 0
b 8 (1, 2) × (1, 0) × (1, −1) −2 2 0 −4 4
c 4 (−1,4) × (0, −1) × (1, −1) 3 −3 0
2 2 (1, 0) × (0, 1) × (0, −2) ${x}_{A}=\tfrac{1}{2}{x}_{B}=\tfrac{1}{4}{x}_{C}=\tfrac{1}{4}{x}_{D}$
${\beta }_{2}^{g}=-2$
${\chi }_{1}=\tfrac{1}{2\sqrt{2}}$, ${\chi }_{2}=\tfrac{1}{\sqrt{2}}$, ${\chi }_{3}=\sqrt{2}$
Table 42. D6-brane configurations and intersection numbers of Model 42, and its gauge coupling relation is ${g}_{a}^{2}={g}_{b}^{2}=2{g}_{c}^{2}=\tfrac{10}{7}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{8}{3}\sqrt[4]{2}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 42 $U(4)\times U{\left(2\right)}_{L}\times U{\left(4\right)}_{R}\times {USp}(2)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 1
a 8 (−1,0) × (−1, −1) × (1, −1) 0 0 0 4 −4 0 0
b 4 (−1, −4) × (−1,0) × (1, −1) −3 3 4 0 0
c 8 (−1, −2) × (0, 1) × (1, 1) −2 2 4
1 2 (1, 0) × (1, 0) × (−2, 0) ${x}_{A}=2{x}_{B}=\tfrac{1}{2}{x}_{C}=\tfrac{1}{2}{x}_{D}$
βg1 = −2
${\chi }_{1}=\tfrac{1}{2\sqrt{2}}$, ${\chi }_{2}=\sqrt{2}$, ${\chi }_{3}=2\sqrt{2}$
Table 43. D6-brane configurations and intersection numbers of Model 43, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{8}{3}{g}_{b}^{2}={g}_{c}^{2}=\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{16}{9}{2}^{3/4}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 43 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times {USp}(2)$
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c 3
a 8 (0, 1) × (1, −1) × (1, −1) 0 0 1 3 −4 0 0
b 4 (−1, 1) × (0, 1) × (−2,4) 2 −2 −9 −5 2
c 4 (4, −1) × (1, 0) × (1, 1) −3 3 0
3 2 (0, 1) × (1, 0) × (0, 2) ${x}_{A}={x}_{B}=\tfrac{1}{2}{x}_{C}=4{x}_{D}$
βg3 = −4
${\chi }_{1}=\sqrt{2}$, ${\chi }_{2}=2\sqrt{2}$, ${\chi }_{3}=\tfrac{1}{\sqrt{2}}$
Table 44. D6-brane configurations and intersection numbers of Model 44, and its gauge coupling relation is ${g}_{a}^{2}=\tfrac{16}{5}{g}_{b}^{2}=2{g}_{c}^{2}=\tfrac{10}{7}\left(\tfrac{5}{3}{g}_{Y}^{2}\right)=\tfrac{16}{5}\sqrt{2}\pi {{\rm{e}}}^{{\phi }^{4}}$.
Model 44 $U(4)\times U{\left(2\right)}_{L}\times U{\left(2\right)}_{R}\times $
stack N (n1, l1) × (n2, l2) × (n3, l3) n□□ b b c c
a 8 (0, 1) × (1, −1) × (1, −1) 0 0 1 3 −4 0
b 4 (−1, 1) × (0, 1) × (−2,4) 2 −2 −6 −6
c 4 (2, −1) × (1, 0) × (2, 2) −2 2
${x}_{A}={x}_{B}=\tfrac{1}{2}{x}_{C}=2{x}_{D}$
χ1 = 1, χ2 = 2, χ3 = 1

TL and CZ are supported by the National Key Research and Development Program of China Grant No. 2020YFC2201504, by the Projects No. 11875062, No. 11947302, and No. 12047 503 supported by the National Natural Science Foundation of China, as well as by the Key Research Program of the Chinese Academy of Sciences, Grant NO. XDPB15. RS is supported by KIAS Individual Grant PG080701. RS would like to thank Weikun He for their useful discussions. CZ would like to thank Lina Wu for their helpful discussions.

1
He W Li T Sun R 2021 The complete search for the supersymmetric Pati–Salam models from intersecting D6-branes arXiv: hep-th/2112.09632

2
Berkooz M Leigh R G 1997 A D = 4 N = 1 orbifold of type I strings Nucl. Phys. B 483 187

DOI

3
Shiu G Tye S-H H 1998 TeV scale superstring and extra dimensions Phys. Rev. D 58 106007

DOI

4
Lykken J Poppitz E Trivedi S P 1999 Branes with GUTs and supersymmetry breaking Nucl. Phys. B 543 105

DOI

5
Cvetic M Plumacher M Wang J 2000 Three family type IIB orientifold string vacua with non-abelian wilson lines J. High Energy Phys. JHEP00(2000)004

DOI

6
Aldazabal G Font A Ibanez L E Violero G 1998 D = 4, N = 1, type iib orientifolds Nucl. Phys. B 536 29

DOI

7
Aldazabal G Ibanez L E Quevedo F Uranga A M 2000 D-branes at singularities : a bottom-up approach to the string embedding of the standard model J. High Energy Phys. JHEP00(2000)002

DOI

8
Klein M Rabadan R 2000 ${{\mathbb{Z}}}_{N}\times {{\mathbb{Z}}}_{M}$ orientifolds with and without discrete torsion J. High Energy Phys. JHEP00(2000)049

DOI

9
Berkooz M Douglas M R Leigh R G 1996 Branes intersecting at angles Nucl. Phys. B 480 265

DOI

10
Blumenhagen R Goerlich L Kors B Lust D 2000 Noncommutative compactifications of type i strings on tori with magnetic background flux J. High Energy Phys. JHEP00(2000)006

DOI

11
Blumenhagen R Kors B Lust D 2001 Type I strings with F- and B-Flux J. High Energy Phys. JHEP01(2001)030

DOI

12
Aldazabal G Franco S Ibanez L E Rabadan R Uranga A M 2001 Intersecting brane worlds J. High Energy Phys. JHEP01(2001)047

DOI

13
Aldazabal G Franco S Ibanez L E Rabadan R Uranga A M 2001 D = 4 chiral string compactifications from intersecting branes J. Math. Phys. 42 3103

DOI

14
Ibanez L E Marchesano F Rabadan R 2001 Getting just the standard model at intersecting branes J. High Energy Phys. JHEP01(2001)002

15
Angelantonj C Antoniadis I Dudas E Sagnotti A 2000 Type-I strings on magnetised orbifolds and brane transmutation Phys. Lett. B 489 223

DOI

16
Forste S Honecker G Schreyer R 2001 Orientifolds with branes at angles J. High Energy Phys. JHEP01(2001)004

DOI

17
Forste S Honecker G Schreyer R 2001 Supersymmetric ${{\mathbb{Z}}}_{N}\times {{\mathbb{Z}}}_{M}$ orientifolds in 4D with D-Branes at angles Nucl. Phys. B 593 127

DOI

18
Blumenhagen R Kors B Lust D Ott T 2001 The standard model from stable intersecting brane world orbifolds Nucl. Phys. B 616 3

DOI

19
Cremades D Ibanez L E Marchesano F 2002 Intersecting brane models of particle physics and the higgs mechanism J. High Energy Phys. JHEP02(2002)022

DOI

20
Cremades D Ibanez L E Marchesano F 2002 Standard model at intersecting d5-branes: lowering the string scale Nucl. Phys. B 643 93

DOI

21
Cremades D Ibanez L E Marchesano F 2002 SUSY quivers, intersecting branes and the modest hierarchy problem J. High Energy Phys. JHEP02(2002)009

DOI

22
Bailin D Kraniotis G V Love A 2003 Intersecting D5-brane models with massive vector-like leptons J. High Energy Phys. JHEP03(2003)052

DOI

23
Bailin D Kraniotis G V Love A 2002 New standard-like models from intersecting D4-branes Phys. Lett. B 547 43

DOI

24
Bailin D Kraniotis G V Love A 2002 Standard-like models from intersecting D4-branes Phys. Lett. B 530 202

DOI

25
Bailin D Kraniotis G V Love A 2003 Standard-like models from intersecting D5-branes Phys. Lett. B 553 79

DOI

26
Kokorelis C 2002 Deformed intersecting D6-brane GUTS I J. High Energy Phys. JHEP02(2002)027

DOI

27
Kokorelis C 2002 Deformed intersecting D6-brane GUTS II arXiv: hep-th/0210200

28
Kokorelis C 2002 Exact standard model compactifications from intersecting branes J. High Energy Phys. JHEP02(2002)036

DOI

29
Kokorelis C 2004 Exact standard model structures from intersecting d5-branes Nucl. Phys. B 677 115

DOI

30
Kokorelis C 2002 New standard model vacua from intersecting branes J. High Energy Phys. JHEP02(2002)029

DOI

31
Cvetic M Shiu G Uranga A M 2001 Chiral four-dimensional N = 1 supersymmetric type IIA orientifolds from intersecting D6-branes Nucl. Phys. B 615 3

DOI

32
Cvetic M Papadimitriou I 2003 More supersymmetric standard-like models from intersecting D6-branes on type IIA orientifolds Phys. Rev. D 67 126006

DOI

33
Cvetic M Papadimitriou I Shiu G 2004 Supersymmetric three family SU(5) grand unified models from type IIA orientifolds with intersecting D6-branes Nucl. Phys. B 696 298

DOI

34
Cvetic M Langacker P Shiu G 2002 Phenomenology of a three-family standard-like string model Phys. Rev. D 66 066004

DOI

35
Cvetic M Shiu G Uranga A M 2001 Three-family supersymmetric standard-like models from intersecting brane worlds Phys. Rev. Lett. 87 201801

DOI

36
Cvetic M Langacker P Wang J 2003 Dynamical supersymmetry breaking in standard-like models with intersecting D6-branes Phys. Rev. D 68 046002

DOI

37
Blumenhagen R Goerlich L Ott T 2003 Supersymmetric Intersecting Branes on the type iia ${{\mathbb{T}}}^{6}/{{\mathbb{Z}}}_{4}$ orientifold J. High Energy Phys. JHEP03(2003)021

DOI

38
Honecker G 2003 Chiral supersymmetric models on an orientifold of ${{\mathbb{Z}}}_{4}\times {{\mathbb{Z}}}_{2}$ with intersecting D6-branes Nucl. Phys. B 666 175

DOI

39
Li T Liu T 2003 Quasi-supersymmetric G3 unification from intersecting D6-branes on type IIA orientifolds Phys. Lett. B 573 193

DOI

40
Cvetic M Li T Liu T 2004 Supersymmetric pati-salam models from intersecting d6-branes: a road to the standard model Nucl. Phys. B 698 163

DOI

41
Cvetic M Langacker P Li T Liu T 2005 D6-brane splitting on type IIA orientifolds Nucl. Phys. B 709 241

DOI

42
Chen C-M Kraniotis G V Mayes V E Nanopoulos D V Walker J W 2005 A K-theory anomaly free supersymmetric flipped SU(5) model from intersecting branes Phys. Lett. B 625 96

DOI

43
Chen C-M Li T Mayes V E Nanopoulos D V 2008 Towards realistic supersymmetric spectra and Yukawa textures from intersecting branes Phys. Rev. D 77 125023

DOI

44
Chen C-M Li T Mayes V E Nanopoulos D V 2008 A realistic world from intersecting D6-branes Phys. Lett. B 665 267

DOI

45
Chen C-M Li T Mayes V E Nanopoulos D V 2008 Realistic yukawa textures and SUSY spectra from intersecting branes Phys. Rev. D 77 125023

DOI

46
Chen C-M Li T Nanopoulos D V 2006 Standard-like model building on type II orientifolds Nucl. Phys. B 732 224

DOI

47
Chen C-M Kraniotis G V Mayes V E Nanopoulos D V Walker J W 2005 A supersymmetric flipped SU(5) intersecting brane world Phys. Lett. B 611 156

DOI

48
Chen C-M Li T Mayes V E Nanopoulos D V 2008 Variations of the hidden sector in a realistic intersecting brane model J. Phys. G: Nucl. Part. Phys. 35 095008

DOI

49
Chen C-M Li T Mayes V E Nanopoulos D V 2008 Yukawa corrections from four-point functions in intersecting D6-brane models Phys. Rev. D 78 105015

DOI

50
Blumenhagen R Cvetic M Langacker P Shiu G 2005 Toward realistic intersecting D-brane models Ann. Review Nucl. Particle Sci. 55 71

DOI

51
Li T Mansha A Sun R 2021 Revisiting the supersymmetric Pati–Salam models from intersecting D6-branes Eur. Phys. J. C 81 82

DOI

52
Li T Mansha A Sun R Wu L He W 2021 N = 1 supersymmetric ${SU}{\left(12\right)}_{C}\times {SU}{\left(2\right)}_{L}\times {SU}{\left(2\right)}_{R}$ models, ${SU}{\left(4\right)}_{C}\times {SU}{\left(6\right)}_{L}\times {SU}{\left(2\right)}_{R}$ models, and ${SU}{\left(4\right)}_{C}\times {SU}{\left(2\right)}_{L}\times {SU}{\left(6\right)}_{R}$ models from intersecting D6-branes Phys. Rev. D 104 046018

DOI

53
Harari H Haut H Weyers J 1978 Quark masses and cabibbo angles Phys. Lett. B 78 459

DOI

54
Sultansoy S 2007 Flavor democracy in particle physics AIP Conf. Proc. 899 49

55
Gross D J Wilczek F 1973 Ultraviolet behavior of non-abelian gauge theories Phys. Rev. Lett. 30 1343

DOI

56
Gronqvist H Particle physics models with four generations PhD Thesis

57
Hung P Q Xiong C 2011 Dynamical electroweak symmetry breaking with a heavy fourth generation Nucl. Phys. B 848 288

DOI

58
Taylor T R 1990 Dilaton, gaugino condensation and supersymmetry breaking Phys. Lett. B 252 59

DOI

59
Brustein R Steinhardt P 1993 Challenges for superstring cosmology Phys. Lett. B 302 196

DOI

60
de Carlos B Casas J A Muñoz C 1993 supersymmetry breaking and determination of the unification gauge coupling constant in string theories Nucl. Phys. B 399 623

DOI

61
He W Li T Sun R Wu L The final model building for the supersymmetric pati-salam models from intersecting D6-branes arXiv: abs/2112.09630

Outlines

/