May 22, 2026
Non-supersymmetric aspects of superstring theory have been attracting growing attention in recent years [1]–[20]. This growing interest is motivated not only by the phenomenological observation that supersymmetry has yet to be discovered in nature, but also by the increasingly prominent role of non-supersymmetric systems in the swampland program on quantum gravity [21]–[26].
Indeed, the construction and classification of examples, the exploration of dualities among them, and even attempts to lift them to M-theory have a long history [27]–[36]. However, the absence of supersymmetry, and hence of protected objects such as BPS states, often makes such attempts difficult.
Why have dualities in supersymmetric string theory worked so remarkably well? One reason is the existence of K3 surfaces and their rich geometry [37]–[44]. Several non-perturbative aspects of Type II strings in higher dimensions come down to K3 geometry: e.g. the realization of non-abelian gauge symmetries in terms of Kodaira singularities [45]. K3 surfaces also shed light on string dualities in lower dimensions through K3-fibered Calabi-Yau threefolds [46]–[48].
There are non-supersymmetric analogues of K3 surfaces, called Enriques surfaces, which can be obtained as fixed-point-free \(\mathbb{Z}_2\) quotients of K3 surfaces. Since this is not Calabi–Yau, it breaks all supersymmetry. Enriques surfaces have often appeared in string-theoretic contexts before, but the constructions were often arranged so as to preserve supersymmetry [49]–[56]. Recently, Type IIA and IIB strings (more precisely Type 0A and Type 0B) on Enriques have been studied in [12], where they have been conjectured to be dual to an asymmetric orbifold of heterotic strings, as non-supersymmetric extensions of known dualities in six and five dimensions.
In the present work, we study heterotic strings on an Enriques surface with 48 inequivalent choices of shift vectors. We construct the worldsheet theory by taking an Enriques orbifold of supersymmetric heterotic strings on a K3 surface, itself realised at an orbifold point. We then show that the resulting models are closely related to ten-dimensional non-supersymmetric heterotic strings [27], [31]. We evaluate the light spectrum of these theories focusing on massless states and moduli-independent tachyons. Although the only tachyon-free parent 10-dimensional theory is the \(D_8 \times D_8\), also some of the theories related to \(E_8 \times D_8, (E_7 \times A_1)^2, A_{15}\times U(1)\) and \(D_4\times D_{12}\) heterotic strings do not inherit the parent’s tachyons. By contrast, all the compactifications descending from the \(D_{16}\) theory retain tachyonic states.
This paper is organized as follows. In section 2 we examine heterotic strings on an Enriques surface. We describe the Enriques involution and use it to define an orbifold action. We construct the partition function and classify inequivalent shift vectors. In section 3 we analyse the light spectrum of the theory, focusing on massless and tachyonic states. In particular, we show that for a class of shift vectors, all moduli-independent tachyons are removed. Additional material is collected into appendices.
We consider the 1-loop partition function of heterotic strings on Enriques surfaces. In our conventions the left-moving sector contains the superstring and the right-moving sector contains the bosonic string.
In this section we review supersymmetric heterotic strings compactified on \(T^4\). Conventions for the Jacobi theta functions can be found in Appendix 4. The Hilbert space of heterotic strings on \(T^4\) is \[\require{physics} \begin{align} \mathcal{H}_{T^4}=&\mathcal{H}_{\text{Boson}}^{4,4}\otimes\mathcal{H}_{\text{Fermion}}^{8,0}\otimes\mathcal{H}_{\Gamma_{4,20}},\\ \mathcal{H}_{\Gamma_{4,20}}=&\mathcal{H}_{\text{Boson}}^{4,20}\otimes \qty(\bigoplus_{p\in\Gamma_{4,20}}\mathbb{C}\ket{p}). \end{align}\] To compute the partition function, let us work at a special point in moduli space, where we can take \[\begin{align} \Gamma_{4,20} = \Gamma_{4,4} \oplus \Gamma_{0,16} \, , \end{align}\] so that the two lattices are separately self-dual and have no mixing. Here, \(\Gamma_{4,4}=H_1(T^4)\oplus H^1(T^4)\) is a \((4,4)\) even self-dual lattice. We will use the same symbol for the lattice and its theta series when no confusion can arise. The partition function is \[\require{physics} \label{eq:hetonT4} \begin{align} &Z_{T^4}=\frac{1}{(\mathrm{Im}\tau)^2\eta^4\bar\eta^4}\frac{\Theta}{\eta^4} \frac{\Gamma_{4,4}}{\eta^4\bar\eta^{4}}\frac{\bar\Gamma_{16}}{\bar\eta^{16}},\\ \Theta(q)=&\qty(\sum_{r\in V}-\sum_{r\in Sp})q^{\frac{1}{2}r^2},\quad \Gamma_{4,4}=\sum_{p \in \Gamma_{4,4}}q^{\frac{1}{2}p_L^2} \bar q^{\frac{1}{2}p^2_R},\quad \bar\Gamma_{16}=\sum_{p\in\Gamma_{0,16}}\bar q^{\frac{1}{2}p^2} \end{align}\tag{1}\] where \(\Theta\) comes from left-moving fermions, and \(V,Sp\) are the vector and spinor lattices of \(SO(8)\): \[\require{physics} \label{eq:V95and95Sp} \begin{align} V=&\left\{(r_1,\cdots,r_4)\in\mathbb{Z}^4\middle|\sum_{i=1}^4 r_i\in2\mathbb{Z}+1\right\},\\ Sp=&\left\{(r_1,\cdots,r_4)\in\qty(\mathbb{Z}+\frac{1}{2})^4\middle|\sum_{i=1}^4 r_i\in2\mathbb{Z}+1\right\}.\\ \end{align}\tag{2}\] The last factor \(\bar\Gamma_{16}\) comes from internal right-moving bosons and should be even and self-dual for modular invariance. Then it is restricted to one of the even self-dual lattices, \(E_8\times E_8\) or \(\mathrm{Spin}(32)/\mathbb{Z}_2=\mathrm{Spin}(32)/\mathbb{Z}_2\): \[\require{physics} \label{eq:32lattices} \begin{align} E_8=&\left\{p\in\mathbb{Z}^8 \text{ or } (\mathbb{Z}+\tfrac{1}{2})^8\middle| \sum_{i=1}^8 p_i\in2\mathbb{Z}\right\} \, ,\\ \mathrm{Spin}(32)/\mathbb{Z}_2=&\left\{p\in\mathbb{Z}^{16}\text{ or }\qty(\mathbb{Z}+\tfrac{1}{2})^{16}\middle|\sum_{i=1}^{16}p_i \in2\mathbb{Z}\right\} \, , \end{align}\tag{3}\] Modular transformations of each factor are given by: \[\begin{align} T\cdot \frac{1}{\eta^{12}\bar\eta^{24}}=&(-1)\frac{1}{\eta^{12}\bar\eta^{24}},&&S\cdot \frac{1}{(\mathrm{Im}\tau)^2\eta^4\bar\eta^4}=\frac{1}{(\mathrm{Im}\tau)^2\eta^4\bar\eta^4}\\ T\cdot \Theta=&(-1)\Theta, &&S\cdot \frac{\Theta}{\eta^4}=\frac{\Theta}{\eta^4},\\ T\cdot\Gamma_{4,4}=&\Gamma_{4,4},&&S\cdot \frac{\Gamma_{4,4}}{\eta^4\bar\eta^4}=\frac{\Gamma_{4,4}}{\eta^4\bar\eta^4}\\ T\cdot \bar\Gamma_{16}=& \bar\Gamma_{16} , &&S\cdot \frac{\bar\Gamma_{16} }{\bar\eta^{16}}=\frac{\bar\Gamma_{16} }{\bar\eta^{16}}, \end{align}\] Therefore \(Z_{T^4}\) is modular invariant.
Starting from heterotic strings on \(T^4\), a singular Enriques surface can be obtained as the orbifold limit of subsequent \(\mathbb{Z}_4\) and \(\mathbb{Z}_2\) actions. This can be thought of as the "orbifold of an orbifold". Let \(g\) denote the orbifold generator of the \(\mathbb{Z}_2\) action and \(h\) the orbifold generator of the \(\mathbb{Z}_4\) action, then \[\begin{align} g: \begin{cases} X^6 \rightarrow-X^6\\ X^7 \rightarrow -X^7\\ X^8 \rightarrow-X^8\\ X^9 \rightarrow-X^9 \end{cases} \, , \quad \begin{cases} \psi^6 \rightarrow-\psi^6\\ \psi^7 \rightarrow -\psi^7\\ \psi^8 \rightarrow -\psi^8\\ \psi^9 \rightarrow -\psi^9 \end{cases} \, ;&\quad \quad h: \begin{cases} X^6 \rightarrow-X^6\\ X^7 \rightarrow -X^7+\pi R_7\\ X^8 \rightarrow X^8\\ X^9 \rightarrow X^9 + \pi R_9, \end{cases} \, , \quad \begin{cases} \psi^6 \rightarrow -\psi^6\\ \psi^7 \rightarrow -\psi^7\\ \psi^8 \rightarrow \psi^8\\ \psi^9 \rightarrow \psi^9 \end{cases} \, , \\ &[g,h]=0 \, , \nonumber \end{align}\] where \(X^i+2\pi R_i=X^i\) with \(i=6,...,9\) are bosonic compact coordinates, \(\psi^{i}\) are worldsheet fermions and all other coordinates are left untouched. It can be easily seen that the action of \(g\) on \(T^4\) provides \(16\) fixed points, and the quotient \(T^4/g\) can be viewed as an orbifold limit of a \(K3\) surface. It is also easily seen that \(h\) acts freely on \(T^4/g\), and it exchanges singular points pairwise. It should be noted that although \(h^2\) is identity on K3 surfaces, it lifts to \(-1\) on the spinor bundle of \(K3\).
The relation between these manifolds is schematically drawn in Figure 1.
For later use, it will be useful to define the same actions on bosonised fermions. Let \(r\) be a spinor or vector weight of \(SO(8)\), then \[\begin{align} g\ket{r} = e^{2 \pi \mathrm{i}r \cdot v_{1f}} \ket{r} \, , \quad h\ket{r}= e^{2 \pi \mathrm{i}r\cdot v_{2f}} \ket{r} \, \\ v_{1f}= \frac{1}{2}(1,1,0,0) \, , \quad v_{2f} = \frac{1}{2} (1,0,0,0) \, . \end{align}\] Here \(v_{1f}\) flips four fermions \(\psi^6,\cdots,\psi^9\), while \(h\) flips just two fermions \(\psi^6,\psi^7\). Our orbifold actions leave the rank 16 lattice completely unrotated, so that the invariant sublattices \(I\subset \Gamma_{4,4}\) with respect to each action have signatures \[\begin{align} \text{sig}\left (I_g \right) = (0,0) \, , \quad \text{sig}\left(I_{h}\right) = (2,2) \, , \end{align}\] where \(\text{vol}(I_h)=2\)1. The actions of \(g,h\) on \(\mathcal{H}_{\Gamma_{4,4}}\) are explained in Appendix 5.
The orbifold action on \(\Gamma_{0,16}\) is defined up to a phase depending on a shift vector. Since the gauge lattice is left geometrically untouched, these shifts can be chosen, without loss of generality, to lie entirely along the \(\Gamma_{0,16}\) directions. We will denote by \(v,w\) the shift vectors related to the \(\mathbb{Z}_2,\mathbb{Z}_4\) action respectively: \[\begin{align} g\ket{p}=&e^{-2\pi ip\cdot v}\ket{p}\,,\\ h\ket{p}=&e^{-2\pi ip\cdot w}\ket{p}\, , \end{align}\] for \(p\in\Gamma_{0,16}\).
The partition function of heterotic strings on an Enriques surface is constructed as \[\begin{align} Z_{\text{Enriques}} =&\frac{1}{(\mathrm{Im}\tau)^2\eta^4\bar\eta^4}\frac{1}{\eta^8 \bar \eta^{20}}\frac{1}{4}\sum_{k,l=0}^3\frac{1}{2}\sum_{c,d=0}^1Z^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}\\ Z^{k,l} \genfrac{[}{]}{0pt}{}{c}{d} =& \tr_{\mathcal{H}^{g^ch^k}} g^d h^l q^{L_0} \bar q^{\bar L_0}\\ =&\Theta^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}\Gamma^{k,l}_{4,4} \genfrac{[}{]}{0pt}{}{c}{d} \overline{\Gamma}_{16}^{k,l} \genfrac{[}{]}{0pt}{}{c}{d} , \label{eq:32op} \end{align}\tag{4}\] where we have separated contributions from fermions, rank 8 and rank 16 lattices. Here \(\mathcal{H}^{g^ch^k}\) denotes the Hilbert space of the \(g^ch^k\)-twisted sector, which means that boundary conditions of fields along spatial \(S^1\) on \(T^4\) are twisted by \(g^ch^k\). The explicit expressions of these building blocks are2 \[\begin{align} \label{zdhtrcgb} \Theta^{k,l} \genfrac{[}{]}{0pt}{}{c}{d} =&e^{-\pi \mathrm{i}(dv_{1f} +lv_{2f})\cdot (cv_{1f}+kv_{2f})} \left(\sum_{r \in V} - \sum_{r \in Sp} \right)q^{\frac{1}{2}(r+cv_{1f}+kv_{2f})^2}e^{2 \pi \mathrm{i}(r+cv_{1f}+kv_{2f}) \cdot (d v_{1f}+ lv_{2f})} \, , \nonumber\\ \overline{\Gamma}^{k,l}_{16} \genfrac{[}{]}{0pt}{}{c}{d} =& e^{\pi \mathrm{i}(dv+lw)\cdot (cv+kw)} \sum_{p \in \Gamma_{0,16}} \bar q^{\frac{1}{2} (p+cv+kw)^2}e^{-2 \pi \mathrm{i}(p+cv + kw)\cdot (dv + lw)} \, , \end{align}\tag{5}\] where the form of \(\Gamma_{4,4}^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}\) is discussed in Appendix 5. In the present paper we focus on \(N=0\). The shift vectors are subject to the following constraints \[\begin{align} \label{eq:32shift32vec32cond} 2 v \in \Gamma_{0,16} \, , \quad 4w \in \Gamma_{0,16} \, , \quad 2v^2 \in 2 \mathbb{Z}+1 \, , \quad 4w^2 \in 2 \mathbb{Z}+1\, , \quad v \cdot 4w \in 4\mathbb{Z}+1 \, . \end{align}\tag{6}\] These conditions follow from imposing invariance of \(Z^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}\) under \(d\to d+2\) and \(l\to l+4\) , which comes from \(g^2=1\) and \(h^4=1\): \[\require{physics} \begin{align} e^{-2\pi \mathrm{i}(2v\cdot p)}e^{-\pi \mathrm{i}c(2v^2+1)} e^{-\pi \mathrm{i}k(2v\cdot w -\frac{1}{2})}=&1,\\ e^{-2\pi \mathrm{i}(4w\cdot p)}e^{-\pi \mathrm{i}c(1+4v\cdot w)}e^{-\pi \mathrm{i}k\qty(1+4w^2)}=&1, \end{align}\] for every \(c,k\in\mathbb{Z}\) and \(p\in\Gamma_{0,16}\). Here we used the invariance of \(\Gamma_{4,4}^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}\) under \(d\to d+2\) and \(l\to l+4\), which can be explicitly checked by using the expressions in Appendix 5.
A comment on the choice of phases in [eq:32enriques32pf32pieces] is in order. In the purely \(g\) or \(h\) orbifolds, when \(k,l=0\) or \(c,d=0\), the reader will recognise the usual orbifold actions described in the previous section. However, in "mixed" terms, new phases appear. For instance, in the \(g\)-twisted sector the action of \(h\) is defined up to an additional phase \(e^{\pi \mathrm{i}v \cdot w}\): \[\begin{align} \overline{\Gamma}^{0,1}_{16} \genfrac{[}{]}{0pt}{}{1}{0}=e^{\pi \mathrm{i}v \cdot w} \sum_{p \in \Gamma_{0,16}} \bar q^{\frac{1}{2} (p+v)^2} e^{-2 \pi \mathrm{i}(p+v) \cdot w} \, . \end{align}\] The phases are chosen so that each piece of the partition function has the following transformation properties under the modular group \[\label{eq:32STtransf} \begin{align} T\cdot \Theta^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}=&(-1)\Theta^{k,k+l} \genfrac{[}{]}{0pt}{}{c}{c+d},&&S\cdot \frac{\Theta^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}}{\eta^4}=\frac{\Theta^{l,-k} \genfrac{[}{]}{0pt}{}{d}{-c}}{\eta^4},\\ T\cdot\Gamma^{k,l}_{4,4} \genfrac{[}{]}{0pt}{}{c}{d}=&\Gamma^{k,k+l}_{4,4} \genfrac{[}{]}{0pt}{}{c}{c+d}, &&S\cdot\frac{\Gamma^{k,l}_{4,4} \genfrac{[}{]}{0pt}{}{c}{d}}{\eta^4\bar\eta^4}=\frac{\Gamma^{l,-k}_{4,4} \genfrac{[}{]}{0pt}{}{d}{-c}}{\eta^4\bar\eta^4},\\ T\cdot \overline{\Gamma}^{k,l}_{0,16} \genfrac{[}{]}{0pt}{}{c}{d}=& \overline{\Gamma}^{k,k+l}_{0,16} \genfrac{[}{]}{0pt}{}{c}{c+d} , &&S\cdot \frac{1}{\bar\eta^{16}}\overline{\Gamma}^{k,l}_{0,16} \genfrac{[}{]}{0pt}{}{c}{d} =\frac{1}{\bar\eta^{16}}\overline{\Gamma}^{l,-k}_{0,16} \genfrac{[}{]}{0pt}{}{d}{-c} \, . \end{align}\tag{7}\] Then \[\begin{align} T\cdot Z^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}=&Z^{k,k+l} \genfrac{[}{]}{0pt}{}{c}{c+d},\\ S\cdot Z^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}=&Z^{l,-k} \genfrac{[}{]}{0pt}{}{d}{-c}, \end{align}\] and the total partition function of heterotic strings on Enriques is modular invariant.
The \(4^22^2=64\) pieces of the partition function can be organised into 8 orbits with respect to \(T,S\) transformations. Denoting the orbit by its generator:
\(Z^{0,0} \genfrac{[}{]}{0pt}{}{0}{0}\): 1 term; this is the partition function of heterotic strings on \(T^4\) 1 ;
\(Z^{0,0} \genfrac{[}{]}{0pt}{}{0}{1}\): 3 terms; these are part of the partition function of heterotic strings on K3 \(=T^4/g\);
\(Z^{0,2} \genfrac{[}{]}{0pt}{}{0}{0}\): 3 terms; these terms correspond to a direct compactification of 10d non-supersymmetric heterotic strings [27] on \(T^4\);
\(Z^{0,1} \genfrac{[}{]}{0pt}{}{0}{0}\): 12 terms; these terms together with the ones in orbit 3 are part of the partition function of heterotic strings on hyperelliptic surface \(T^4/h\) [1];
\(Z^{0,1} \genfrac{[}{]}{0pt}{}{0}{1}\): 12 terms.
\(Z^{0,2} \genfrac{[}{]}{0pt}{}{0}{1}\): 3 terms.
\(Z^{0,2} \genfrac{[}{]}{0pt}{}{1}{0}\): 6 terms.
\(Z^{0,1} \genfrac{[}{]}{0pt}{}{1}{0}\): 24 terms; all of these vanish because \(\Gamma_{4,4}^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}=0\) in this orbit; see Appendix 5.
Using the conditions derived in 6 , it is possible to classify inequivalent shift vectors. The complete derivation is presented in Appendix 7, here we summarise the final result. Throughout this section we fix the \(\mathbb{Z}_2\) shift to be \[\begin{align} v=\frac{1}{2}(0^{14},1,1) \, , \end{align}\] in the canonical basis.
Classification of shift vectors of \(E_8\times E_8\) is summarized in Appendix 7. As a result, we have \(24\) inequivalent shifts:
\[\begin{align} \text{E_8\times D_8}: &\quad \{\mathsf A_{1},\mathsf A_{2},\mathsf A_{3}\}\times\{\mathsf B_{1},\mathsf B_{2}\} ,\quad \{\mathsf A_{4},\mathsf A_{5}\}\times\{\mathsf B_{3},\mathsf B_{4}\},\\ \text{D_8\times D_8}: &\quad \{\mathsf A_{4},\mathsf A_{5}\}\times\{\mathsf B_{5}\} ,\quad \{\mathsf A_{6}\}\times\{\mathsf B_{1},\mathsf B_{2}\},\\ \text{(E_7\times A_1)^2}: &\quad \{\mathsf A_{7},\mathsf A_{8}\}\times\{\mathsf B_{6},\mathsf B_{7},\mathsf B_{8}\} ,\quad \{\mathsf A_{9},\mathsf A_{10}\}\times\{\mathsf B_{9},\mathsf B_{10}\}. \end{align}\] Here \(E_8\times D_8, D_8\times D_8\) and \((E_7\times A_1)^2\) denote the corresponding ten-dimensional non-supersymmetric parent heterotic theories.
Here \(A_i\) and \(B_j\) denote the components of the shift \(w=(A_i:B_j)\in \frac{1}{4}(E_8\times E_8)\) in the first and second \(E_8\) factors, respectively. Their unbroken algebras are summarized as
\[\begin{array}{@{}l@{\qquad}l@{\qquad\qquad}l@{\qquad}l@{}} \mathsf A_1: E_8 & \mathsf B_1: D_4\times A_1^2\times U(1)^2 & \mathsf A_6: A_7\times U(1) & \mathsf B_6: A_5\times U(1)^3 \\ \mathsf A_2: E_7\times A_1 & \mathsf B_2: D_6\times U(1)^2 & \mathsf A_7: E_7\times U(1) & \mathsf B_7: D_6\times A_1\times U(1) \\ \mathsf A_3: D_8 & \mathsf B_3: E_6\times U(1)^2 & \mathsf A_8: D_6\times A_1\times U(1) & \mathsf B_8: E_7\times U(1) \\ \mathsf A_4: D_7\times U(1) & \mathsf B_4: A_7\times U(1) & \mathsf A_9: E_6\times A_1\times U(1) & \mathsf B_9: D_5\times A_1\times U(1)^2 \\ \mathsf A_5: D_5\times A_3 & \mathsf B_5: A_5\times A_1\times U(1)^2 & \mathsf A_{10}: A_7\times A_1 & \mathsf B_{10}: A_3^2\times A_1\times U(1). \end{array} \label{eq:list95of95AB}\tag{8}\]
For given \(v,w\) and some original root system \(R\) then the root system of the unbroken gauge algebra \(R_{v,w}\) is \[R_{v,w}=\{\alpha\in R|\alpha\cdot v,\alpha\cdot w\in\mathbb{Z}\} \, .\] Since roots transform under the orbifold action with phases \(e^{2 \pi \mathrm{i}\alpha \cdot v},e^{2 \pi \mathrm{i}\alpha \cdot w}\), \(R_{v,w}\) is the sub-system of roots left invariant by the shifts \(v,w\).
Similarly, there are \(24\) inequivalent \(w\)-shifts on \(\mathrm{Spin}(32)/\mathbb{Z}_2\). In the case of \(4w\in \mathbb{Z}^{16}\), there are \(10+8+2=20\) shift vectors: \[\require{physics} \begin{align} W_{0,\frac{1}{2}}(n_1,n_2)=& \qty(\frac{1}{4}\qty(0^{14-n_1-n_2},1^{n_1},2^{n_2}),0,\frac{1}{2}):\\ (n_1,n_2)\in\Bigl\{ (0,0),(0,2)&,(0,4),(0,6), (4,1),(4,3),(4,5), (8,0),(8,2),(12,1) \Bigr\},\\ \end{align}\] and \[\require{physics} \begin{align} W_{\frac{1}{4},\frac{1}{4}}(n_1,n_2)=&\qty(\frac{1}{4}\qty(0^{14-n_1-n_2},1^{n_1},2^{n_2}),\frac{1}{4},\frac{1}{4}):\\ (n_1,n_2)\in\Bigl\{ (2,0),(2,2)&,(2,4),(2,6), (6,1),(6,3),(10,0),(10,2) \Bigr\} \, . \end{align}\] and
\[\require{physics} \widetilde{W}_{\frac{1}{4},\frac{1}{4}}(n)=\qty(-\frac{1}{2},\frac{1}{2}^{11-8n},\frac{1}{4}^{2+8n},\frac{1}{4},\frac{1}{4}),\quad n=0,1.\] The unbroken gauge algebras can be identified by using Weyl subgroup of \(D_{14}\times A_1\) which stabilize \(W_{\ast,\ast}(n_1,n_2)\), as follows: \[\begin{align} W_{0,1/2}(n_1,n_2)&:\quad D_{14-n_1-n_2}\times A_{n_1-1}\times D_{n_2}\times \mathrm{U}(1)^{3},\\ W_{1/4,1/4}(n_1,n_2)&:\quad D_{14-n_1-n_2}\times A_{n_1-1}\times A_1\times D_{n_2}\times U(1)^{2},\\ \widetilde{W}_{\frac{1}{4},\frac{1}{4}}(n)&:D_{12-8n}\times A_{1+8n}\times A_1\times U(1)^2 \end{align}\] Here, \(D_1\) means \(U(1)\).
The shifts in the odd class are \[\require{physics} \begin{align}\label{eq:new95shifts} W_{\frac{1}{8},\frac{3}{8}}(n)=&\left(\qty(\frac{1}{8})^{11-4n},\qty(\frac{3}{8})^{3+4n},\frac{1}{8},\frac{3}{8}\right) ,\quad n=0,1,\\ W_{\frac{5}{8},-\frac{1}{8}}(n)=& \left( \left(\frac{1}{8}\right)^{13-4n}, \left(\frac{3}{8}\right)^{1+4n}, \frac{5}{8}, -\frac{1}{8} \right), \quad n=0,1, \end{align}\tag{9}\]
with gauge algebras \[\begin{align} A_{10-4n}\times A_{2+4n}\times U(1)^{4},\\ A_{12-4n}\times A_{4n}\times U(1)^4, \end{align}\] respectively.
These shift vectors for \(\mathrm{Spin}(32)/\mathbb{Z}_2\) with their parent non-supersymmetric heterotic strings and gauge algebra are summarised in table 1.
| shift | 10d parent | gauge algebra |
|---|---|---|
| \(W_{0,\frac12}(0,n)\) | \(D_{16}\) | \(D_{14-n}\times D_n\times \mathrm{U}(1)^2,\quad n=0,2,4,6\) |
| \(W_{0,\frac12}(4,n)\) | \(D_4\times D_{12}\) | \(D_{10-n}\times A_3\times D_n\times \mathrm{U}(1)^3,\quad n=1,3,5\) |
| \(W_{0,\frac12}(12,1)\) | \(D_4\times D_{12}\) | \(A_{11}\times \mathrm{U}(1)^5\) |
| \(W_{\frac14,\frac14}(2,n)\) | \(D_4\times D_{12}\) | \(D_{12-n}\times A_1^2\times D_n\times \mathrm{U}(1)^2,\quad n=0,2,4,6\) |
| \(W_{\frac14,\frac14}(10,n)\) | \(D_4\times D_{12}\) | \(D_{4-n}\times A_9\times D_n\times A_1\times \mathrm{U}(1)^2,\quad n=0,2\) |
| \(W_{0,\frac12}(8,n)\) | \(D_8\times D_8\) | \(D_{6-n}\times A_7\times D_n\times \mathrm{U}(1)^3,\quad n=0,2\) |
| \(W_{\frac14,\frac14}(6,n)\) | \(D_8\times D_8\) | \(D_{8-n}\times A_5\times A_1\times D_n\times \mathrm{U}(1)^2,\quad n=1,3\) |
| \(\widetilde{W}_{\frac14,\frac14}(n)\) | \(D_4\times D_{12}\) | \(D_{12-8n}\times A_{1+8n}\times A_1\times \mathrm{U}(1)^2,\quad n=0,1\) |
| \(W_{\frac18,\frac38}(n)\) | \(A_{15}\times \mathrm{U}(1)\) | \(A_{10-4n}\times A_{2+4n}\times \mathrm{U}(1)^4,\quad n=0,1\) |
| \(W_{\frac58,-\frac18}(n)\) | \(A_{15}\times \mathrm{U}(1)\) | \(A_{12-4n}\times A_{4n}\times \mathrm{U}(1)^4,\quad n=0,1\) |
Heterotic string theories without spacetime supersymmetry in 10 dimensions can be constructed as \(\mathbb{Z}_2\) asymmetric orbifolds of supersymmetric ones [27]. These orbifolds act only on the internal compact degrees of freedom and hence do not reduce the number of spacetime dimensions. In total, there are seven inequivalent theories corresponding to seven inequivalent shift vectors \(\delta \in \frac{1}{2} \Gamma_{0,16}\). The shift vectors fall into two classes: \(\delta^2=1\) giving rise to tachyonic models and \(\delta^2=2\) that lead to tachyon-free theories. For \(\mathrm{Spin}(32)/\mathbb{Z}_2\), the inequivalent shifts are \[\require{physics} \begin{align} &\delta=\qty(0^{15},1): D_{16} \quad &&\delta=\qty(0^{12},\qty(\frac{1}{2})^4):D_4\times D_{12} ,\\ &\delta=\qty(\qty(\frac{1}{4})^{16}):A_{15}\times U(1) , &&\delta=\qty(\qty(\frac{1}{2})^8,0^8):D_8\times D_8 ,\\ \end{align}\] where the last one corresponds to \(\delta^2=2\). We have also listed the unbroken gauge subalgebra of \(\mathrm{Spin}(32)/\mathbb{Z}_2\) after the orbifold. For \(E_8\times E_8\) the shifts are \[\require{physics} \begin{align} &\delta=(0^8;0^7,1): &&D_8\times E_8\\ &\delta=\qty(0^6,\frac{1}{2},\frac{1}{2};0^6,\frac{1}{2},\frac{1}{2}):&&(E_7\times A_1)^2\\ &\delta=(1,0^7;1,0^7): &&D_8\times D_8 \, , \end{align}\] where the last one corresponds to \(\delta^2=2\). For each shift we list the unbroken gauge subalgebra of \(E_8\times E_8\).
Except for \(D_8\times D_8\) strings, the other six heterotic strings have tachyons. This means that the perturbative expansion is around an unstable vacuum. Tachyon condensation of them was argued in [4], [58] and their (meta)-stable vacuum is in lower dimensions. It was recently shown that these non-supersymmetric heterotic strings describe dynamics of near-horizon limit of heterotic branes [22], [25], which are understood in the context of swampland program in quantum gravity [21]. Fascinating relations of (non-supersymmetric) heterotic strings to a kind of generalized cohomology, called topological modular forms [59]–[61], were pointed out in [10], [62].
We can then associate each Enriques shift \(w\) satisfying \(2w=\delta\) modulo vectors in \(\Gamma_{0,16}\) to its "parent" 10-dimensional non-supersymmetric heterotic theory with shift \(\delta\). Notice, however, that the shift vector \(w\) breaks not only the original \(E_8\times E_8,\mathrm{Spin}(32)/\mathbb{Z}_2,\) but also part of gauge symmetry of 10d non-susy heterotic. It would then be natural to expect that only the Enriques shift vectors with unbroken gauge algebra \(D_8 \times D_8\) will give tachyon-free theories. Surprisingly, we will find more moduli-independent-tachyon-free theories than naively expected, as discussed in section 3.
Enriques surfaces do not admit spin structures. Indeed, the free part of the second cohomology lattice of an Enriques surface \(X\) has intersection form \[H^2(X,\mathbb{Z})_{\rm free}\simeq E_8(-1)\oplus \Gamma_{1,1},\] and hence signature \(\sigma=-8\). If \(X\) were spin, Rokhlin’s theorem would imply \(\sigma\in 16\mathbb{Z}\), which is not the case. Therefore Enriques surfaces are non-spin.
The absence of a spin structure on Enriques surfaces may be one of the reasons why string compactifications on them have not been discussed so far in the literature. Nevertheless, superstring theories can be formulated on Enriques surfaces in certain senses. For example, taking Enriques orbifold of Type IIA on a K3 surface yields Type 0A string theory on an Enriques surface [12]. This is a natural argument because there are no fermions in Type 0A, at least perturbatively.
A more sophisticated answer will be provided in the case of heterotic strings on a non-spin manifold [63]. In the present work we define the theories as exactly solvable orbifold CFTs. The question of its interpretation as a smooth heterotic compactification on a non-spin Enriques surface, including possible global anomalies and gauge-bundle/torsion data, will be treated in the future.
In this section we discuss the tachyonic and massless spectrum of heterotic strings on an Enriques surface. We denote the momenta by \(p\) taking values in \(\overline{\Gamma}^{k,l}_{16} \genfrac{[}{]}{0pt}{}{c}{d}\), \(q_{L,R}\) the momenta taking values in \(\Gamma^{k,l}_{4,4} \genfrac{[}{]}{0pt}{}{c}{d}\) and \(r\) the \(SO(8)\) weights. Since we are interested in the light spectrum, we will analyse the mass formulas for the Neveu-Schwarz (NS) and Ramond (R) states with minimal contribution. To this end, the leading powers of \(q\) in \(\Theta^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}\) are summarised in table 2. We will refer to tachyonic states that are present throughout moduli space as moduli-independent, whereas tachyons that occur only on special loci in moduli space will instead be called moduli-dependent. In the following, we focus on moduli-independent tachyons and massless states. Our aim is to show that appropriate choice of shift vector \(w\) can remove moduli-independent tachyons.
| \(g^c h^k\)-twisted sector | NS leading term | R leading term |
|---|---|---|
| \(1\) | \(q^{1/2}\) | \(8\, q^{1/2}\) |
| \(h\) | \(q^{1/8}\) | \(4q^{3/8}\) |
| \(h^2\) | \(q^0\) | \(8\, q^{1/2}\) |
| \(h^3\) | \(q^{1/8}\) | \(4q^{3/8}\) |
| \(g\) | \(q^{1/4}\) | \(2q^{1/4}\) |
| \(gh\) | \(q^{1/8}\) | \(q^{3/8}\) |
| \(gh^2\) | \(q^{1/4}\) | \(2q^{1/4}\) |
| \(gh^3\) | \(q^{1/8}\) | \(q^{3/8}\) |
To make discussions about projection phases more compact in the following, let us define \[\begin{align} \label{eq:32phases} \Phi^{k,l} \genfrac{[}{]}{0pt}{}{c}{d} =& e^{-\pi \mathrm{i}\,(d v_{1f}+l v_{2f})\cdot(c v_{1f}+k v_{2f})} e^{\pi \mathrm{i}\,(d v+l w)\cdot(c v+k w)}\\ \times&e^{2\pi \mathrm{i}\,(r+c v_{1f}+k v_{2f})\cdot(d v_{1f}+l v_{2f})} e^{-2\pi \mathrm{i}\,(p+c v+k w)\cdot(d v+l w)} \, . \end{align}\tag{10}\]
The masses are \[\begin{align} m^2_{L} = \frac{1}{2}q_L^2 +\frac{1}{2}r^2+N_L-\frac{1}{2} \, ; \quad m^2_R= \frac{1}{2}q^2_R + \frac{1}{2}p^2+N_R-1 \, , \end{align}\] where the oscillators take integer values. For the lightest NS and R states one finds \[\begin{align} m^2_L = N_L \, , \end{align}\] which means that there are no tachyons, although massless states are present.
The masses are \[\begin{align} m^2_{L} = \frac{1}{2}(r+v_{1f})^2+N_L-\frac{1}{4} \, ; \quad m^2_R= \frac{1}{2}(p+v)^2+N_R-\frac{3}{4} \, , \end{align}\] where the oscillators take half-integer values. Lightest NS and R states satisfy \[\begin{align} m^2_L = N_L \, , \end{align}\] so this twisted sector is likewise tachyon-free while containing massless states.
The masses are \[\begin{align} m^2_{L} = \frac{1}{2} q^2_L+\frac{1}{2}(r+v_{2f})^2+N_L-\frac{1}{4} \, ; \quad m^2_R= \frac{1}{2} q^2_R+\frac{1}{2}(p+w)^2+N_R-\frac{3}{4} \, , \end{align}\] where the oscillators take half-integer values. The minimal contributions in NS and R sectors have mass \[\begin{align} \label{eq:32h-twisted32mass} m^2_L(\text{NS}) = \frac{1}{2} q^2_L+N_L-\frac{1}{8} \, ; \quad m^2_L(\text{R})=\frac{1}{2} q^2_L + N_L +\frac{1}{8} \, . \end{align}\tag{11}\] Therefore in the R sector there are no massless nor tachyonic states. Moreover, since in the \(h\)-twisted sector the winding along \(R_9\) is shifted (see Appendix 5), \(q^2=0\) is not a point in the lattice. Consequently, the NS sector may contain moduli-dependent tachyons if \[\begin{align} q^2_L < \frac{1}{4} \, , \end{align}\] but not moduli-independent ones. Equivalent results apply to the \(h^3\)-twisted sector.
The masses are \[\begin{align} m^2_{L} = \frac{1}{2} q^2_L+\frac{1}{2}(r+v_{1f}+v_{2f})^2+N_L-\frac{1}{4} \, ; \quad m^2_R= \frac{1}{2} q^2_R+\frac{1}{2}(p+v+w)^2+N_R-\frac{3}{4} \, , \end{align}\] where the oscillators take half-integer values. The lightest states in NS and R sectors have the same mass as the ones in the \(h\)-twisted sector (11 ). In this sector the winding along \(R_7\) is shifted which means that again there are moduli-dependent tachyons and no massless states. The same results hold in the \(gh^3\)-twisted sector.
The masses are \[\begin{align} m^2_{L} = \frac{1}{2} q^2_L+\frac{1}{2}(r+2v_{2f})^2+N_L-\frac{1}{2} \, ; \quad m^2_R= \frac{1}{2} q^2_R+\frac{1}{2}(p+2w)^2+N_R-1 \, , \end{align}\] where the oscillators take integer values. Lightest NS and R states satisfy \[\begin{align} m^2_{L}(\text{NS})=\frac{1}{2}q^2_L+N_L-\frac{1}{2} \, ; \quad m^2_L(\text{R}) = \frac{1}{2} q^2_L + N_L \, . \end{align}\] In the R sector there are massless states and no tachyonic states. In the NS sector, in addition to moduli-dependent tachyons, there is also a moduli-independent tachyon with mass \(m^2=-\frac{1}{2}\). Level-matching imposes \((p+2w)^2=1\).
The masses are \[\begin{align} m^2_{L} = \frac{1}{2}(r+v_{1f}+2v_{2f})^2+N_L-\frac{1}{4} \, ; \quad m^2_R= \frac{1}{2}(p+v+2w)^2+N_R-\frac{3}{4} \, , \end{align}\] where the oscillators take half-integer values. Lightest NS and R states satisfy \[\begin{align} m^2_L = N_L \, , \end{align}\] so this twisted sector is tachyon-free while containing massless states.
In this subsection we determine which Enriques shift vectors \(w\) remove tachyonic states which are present independently of the value of radii \(R_6,\cdots,R_9\). As discussed above, the only sector in which moduli-independent tachyons can occur is the \(h^2\)-twisted sector. In this sector the lightest NS states have
\[m_L^2({\rm NS}) = \frac{1}{2} q_L^2+N_L-\frac{1}{2} , \qquad m_R^2 = \frac{1}{2} q_R^2+\frac{1}{2}(p+2w)^2+N_R-1 .\]
Thus a moduli-independent tachyon has \[q_L=q_R=0,\qquad N_L=N_R=0, \qquad m^2=-\frac{1}{2} ,\]
and level matching imposes \[(p+2w)^2=1, \qquad p\in \Gamma_{0,16}. \label{eq:h2-tachyon-level-matching}\tag{12}\] It remains to impose the orbifold projection.
For a candidate tachyon in the \(h^2\)-twisted sector, namely \(c=0\), \(k=2\), the leading NS contribution of the fermionic block carries the phase \[\exp\!\left[-\frac{\pi i}{2}(d+l)\right].\] Therefore the projection factor for a state labelled by \(p\) is \[\begin{align} & \frac{1}{8} \sum_{d=0}^1\sum_{l=0}^3 \Phi^{2,l} \genfrac{[}{]}{0pt}{}{0}{d}\\ =& \frac{1}{8} \sum_{d=0}^{1} \sum_{l=0}^{3} \exp\!\left[ -2\pi i d\left(p\cdot v+\frac{1}{2}\right) \right] \exp\!\left[ -2\pi i l\left(p\cdot w+w^2+\frac{1}{4}\right) \right], \label{eq:h2-tachyon-projection} \end{align}\tag{13}\] where in the second line we used \(v\cdot w=\frac{1}{4}\), as is the case for the representatives listed in Section 2.4. Hence a level-matched candidate survives the orbifold projection if and only if \[(p+2w)^2=1, \qquad 2p\cdot v \equiv 1 \pmod 2, \qquad p\cdot w+w^2+\frac{1}{4}\in\mathbb{Z} . \label{eq:tachyon-survival-conditions}\tag{14}\] The last condition can be rewritten in a useful \(w\)-independent form. Indeed, using 12 and \(p^2\in 2\mathbb{Z}\), we have \[p\cdot w+w^2+\frac{1}{4} = \frac{2-p^2}{4}.\] Thus the tachyon survives precisely when \[(p+2w)^2=1, \qquad 2p\cdot v \equiv 1 \pmod 2, \qquad p^2\equiv 2 \pmod 4 . \label{eq:tachyon-survival-congruences}\tag{15}\] We now apply this criterion to all admissible Enriques shift vectors. There are two possible mechanisms by which the moduli-independent tachyon is removed. First, the level-matched candidates \[\left\{ p\in\Gamma_{0,16} \;\middle|\; (p+2w)^2=1 \right\},\] may be empty. Second, there are level-matched candidates, but all the states may fail the congruence conditions 15 and are therefore projected out.
\[\renewcommand{\arraystretch}{1.35} \begin{array}{c|l} \text{mechanism} & \text{shift classes} \\ \hline (p+2w)^2=1\;\text{has no solutions} & (\mathsf A_4,\mathsf B_5),\;(\mathsf A_5,\mathsf B_5),\;(\mathsf A_6,\mathsf B_1),\;(\mathsf A_6,\mathsf B_2) \\[1mm] \text{all level-matched candidates are projected out} & (\mathsf A_2,\mathsf B_2),\; \{\mathsf A_7,\mathsf A_8\}\times\{\mathsf B_7,\mathsf B_8\}, (\mathsf A_9,\mathsf B_9),\;(\mathsf A_{10},\mathsf B_{10}) \end{array}\]
For the \(E_8\times E_8\) lattice, the moduli-independent-tachyon-free shift classes are summarized in the table ¿tbl:tab:E8E895shifts?
Thus \(11\) of the \(24\) inequivalent \(E_8\times E_8\) Enriques shift classes are free of moduli-independent tachyons.
For the \(\mathrm{Spin}(32)/\mathbb{Z}_2\) lattice, the moduli-independent-tachyon-free shift classes are summarized in the table ¿tbl:tab:so3295shift?. Therefore \(9\) of the \(24\) inequivalent \(\mathrm{Spin}(32)/\mathbb{Z}_2\) Enriques shift classes are free of moduli-independent tachyons.
\[\renewcommand{\arraystretch}{1.35} \begin{array}{c|l} \text{mechanism} & \text{shift classes} \\ \hline (p+2w)^2=1\;\text{has no solutions} & W_{0,\frac12}(8,0),\; W_{0,\frac12}(8,2),\; W_{\frac14,\frac14}(6,1),\; W_{\frac14,\frac14}(6,3) \\[1mm] \text{all level-matched candidates are projected out} & W_{\frac14,\frac14}(2,2),\; W_{\frac14,\frac14}(2,6),\; W_{\frac14,\frac14}(10,0),\; W_{\frac18,\frac38}(0),\; W_{\frac18,\frac38}(1) \end{array}\]
We emphasise that this analysis removes tachyons which are present throughout the \(\Gamma_{4,4}\) moduli space. It does not exclude moduli-dependent tachyons.
In this section we analyse the massless spectrum of heterotic strings on Enriques. We restrict to solutions with \(q_{L,R}^2=0\). At special loci in the moduli space of admissible \(\Gamma_{4,4}\) lattices the number of massless states might be enhanced from non-vanishing values of \(q_{L,R}\); these states will be disregarded in the following.
Massless states are organised into representations of the six dimensional little group \(\mathrm{Spin}(4) \sim SU(2) \times SU(2)\). Representations of \(SU(2) \times SU(2)\) will be denoted \((\bullet,\bullet)\), where each entry gives the dimension of the corresponding \(SU(2)\) representation.
As an example, consider the states associated to \(N_L=1\) in heterotic compactified on \(T^4/G\). In 10 dimensions the right-moving excitations comprise 16 internal oscillators in the trivial representations of \(SO(8)\) and 8 spacetime oscillators in the vector \(\mathbf{8_v}\) of \(SO(8)\). Upon branching \[\begin{align} SO(8) \longrightarrow SO(4)\simeq SU(2)\times SU(2) \, , \end{align}\] the six-dimensional states organise as \((\mathbf{2},\mathbf{2})+20(\mathbf{1},\mathbf{1})\).
The computation of massless spectrum proceeds in three steps. First, we solve the level-matched conditions derived at the beginning of this section. We then branch the relevant \(SO(8)\) representations to \(SU(2)\times SU(2)\). Finally, multiplicities are determined from the phases appearing in the partition function ([eq:32enriques32pf32pieces]), which have been collected in the function \(\Phi^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}\) 10 .
We will use the following relations: \(v_{1f}\cdot v_{2f}=\frac{1}{4},\,\, v_{1f}^2=\frac{1}{2}, \, \, v_{2f}^2=\frac{1}{4}\) and \(v\cdot w=\frac{1}{4}\). The results in this section are completely general, except for a class of solutions in the untwisted sector (see table 3) which has been analysed as an example on \(E_8 \times E_8\).
The analysis splits into two classes of level-matched solutions \[\begin{align} r^2=1, q^2_{L}=0=N_L, p^2=0=q^2_R, N_R=1 \, ; \quad r^2=1, q^2_{L}=0=N_L, N_R=0=q^2_R, p^2=2 ,\;. \end{align}\] When \(p^2=0,r^2=1\) the phases simplify as \[\begin{align} \frac{1}{8}\sum_{l=0}^3\sum_{d=0}^1\Phi^{0,l} \genfrac{[}{]}{0pt}{}{0}{d} = \frac{1}{8} \left(1 +e^{2 \pi \mathrm{i}r \cdot v_{1f}}+ \sum_{j=1}^3\big(e^{2 \pi \mathrm{i}r \cdot j v_{2f}} + e^{2 \pi \mathrm{i}r \cdot (v_{1f}+j v_{2f})} \big) \right) \, . \end{align}\] All states in the R sector are projected out, while the surviving NS states combine into \((\mathbf{2},\mathbf{2})\). We also have \(N_{R}=1\) yielding \((\mathbf{2},\mathbf{2})+20(\mathbf{1},\mathbf{1})\), so one finds \[\begin{align} \{ (\mathbf{2},\mathbf{2}) \}_{L} \times \{ (\mathbf{2},\mathbf{2})+20(\mathbf{1},\mathbf{1}) \}_{R} = 20 (\mathbf{2},\mathbf{2})+(\mathbf{3},\mathbf{3})+(\mathbf{3},\mathbf{1})+(\mathbf{1},\mathbf{3})+(\mathbf{1},\mathbf{1}) \, . \end{align}\] For the other solution we have to consider all possible values of dot products between the 480 roots of \(E_8\times E_8\) or \(\mathrm{Spin}(32)/\mathbb{Z}_2\) and the shift vectors. As an example, the results for \(E_8\times E_8\) are \[\begin{align} v\cdot\alpha= \left \{ 0,\pm \frac{1}{2}, \pm 1 \right \} \, ; \quad w\cdot\alpha= \left \{ 0, \pm\frac{1}{4},\pm \frac{1}{2}, \pm \frac{3}{4}, \pm 1, \pm \frac{5}{4}, \pm\frac{7}{4} \right \} \, . \end{align}\] The resulting states are summarised in table 3.
| sector | \(p \cdot v\) | \(p \cdot w\) | state |
|---|---|---|---|
| NS | \(0, \pm 1\) | \(0,1\) | \((\mathbf{2},\mathbf{2})\) |
| NS | \(0, \pm 1\) | \(\pm \frac14, \pm \frac12, \pm \frac34, \pm \frac54, \pm \frac74\) | - |
| NS | \(\pm \frac12\) | \(0, \pm \frac12, \pm1\) | \((\mathbf{1},\mathbf{1})\) |
| NS | \(\pm \frac12\) | \(\pm \frac14, \pm \frac34, \pm \frac54, \pm \frac74\) | - |
| R | \(0,\pm1\) | \(\pm \frac14, \pm \frac34, \pm \frac54, \pm \frac74\) | \((\mathbf{1},\mathbf{2})\) |
| R | \(0, \pm 1\) | \(0, \pm\frac12, \pm 1\) | - |
| R | \(\pm \frac12\) | \(\pm \frac14, \pm \frac34, \pm \frac54, \pm \frac74\) | \((\mathbf{2},\mathbf{1})\) |
| R | \(\pm \frac12\) | \(0, \pm\frac12, \pm 1\) | - |
There are two solutions to the mass constraints \[\begin{align} (r+v_{1f})^2=\frac{1}{2}, N_L=0=N_R, (p+v)^2=\frac{3}{2}\, ; \quad (r+v_{1f})^2=\frac{1}{2},N_L=0, (p+v)^2=\frac{1}{2}=N_R \, , \end{align}\] which behave differently under the orbifold phases. While the second is completely projected out, the first one contributes non-trivially. Taking into account the fixed point degeneracy of 16, it yields \[\begin{align} \frac{16}{8}\sum_{l=0}^3\sum_{d=0}^1\Phi^{0,l} \genfrac{[}{]}{0pt}{}{1}{d} =16\cdot \frac{1+h^2}{4}=\begin{cases} 4+4e^{-4 \pi \mathrm{i}(p\cdot w)} \quad \text{NS}\\ 4-4e^{-4 \pi \mathrm{i}(p\cdot w)} \quad \text{R} \end{cases} \, , \end{align}\] where the lightest eigenstates have eigenvalue of \(g\) equal to 1, so the projection factor \(\frac{1+g}{2}\) simplifies. Moreover, the insertion of \(h,h^3\) in \(g\)-twisted sector vanishes, as derived in Appendix 5 : \(\Gamma_{4,4}^{0,l} \genfrac{[}{]}{0pt}{}{1}{d}=0\) for \(l=1,3\). Supposing there are \(k_1\) solutions to \((p+v)^2=\frac{3}{2}\) such that \(p\cdot4w\) is even and \(k_2\) such that it is odd, the resulting spectrum is \[\begin{align} 16 k_1(\mathbf{1},\mathbf{1}) + 8k_2(\mathbf{2},\mathbf{1}) \, . \end{align}\]
There is a single level-matched solution \[\begin{align} (r+2v_{2f})^2 =1, N_L=0=q^2_L, (p+2w)^2=2, N_{R}=0=q^2_R \, , \end{align}\] corresponding to \(2(\mathbf{1},\mathbf{2})+2(\mathbf{2},\mathbf{1})\). Supposing that there are
\(m_1\) solutions to \((p+2w)^2=2\) such that \(2p\cdot w=\frac{1}{2}\) mod 2 and \(2p \cdot v= 2 \mathbb{Z}\);
\(m_2\) solutions to \((p+2w)^2=2\) such that \(2p\cdot w=\frac{1}{2}\) mod 2 and \(2p \cdot v= 2 \mathbb{Z}+1\);
\(m_3\) solutions to \((p+2w)^2=2\) such that \(2p\cdot w=\frac{3}{2}\) mod 2 and \(2p \cdot v= 2 \mathbb{Z}\);
\(m_4\) solutions to \((p+2w)^2=2\) such that \(2p\cdot w=\frac{3}{2}\) mod 2 and \(2p \cdot v= 2 \mathbb{Z}+1\),
the multiplicities of states are \[\begin{align} \frac{1}{8}\sum_{l=0}^3\sum_{d=0}^1\Phi^{2,l} \genfrac{[}{]}{0pt}{}{0}{d} =\begin{cases} \frac{1}{2} e^{-4 \pi \mathrm{i}w^2} (m_1+m_3) \left(-1+e^{4 \pi \mathrm{i}w^2} \right) \quad \text{for} \, \, (\mathbf{1},\mathbf{2})\\ \frac{1}{2} e^{-4 \pi \mathrm{i}w^2} (m_2+m_4) \left(-1+e^{4 \pi \mathrm{i}w^2} \right) \quad \text{for} \, \, (\mathbf{2},\mathbf{1}) \end{cases} \quad . \end{align}\] The shift vectors \(w\) on \(E_8\times E_8\) and \(\mathrm{Spin}(32)/\mathbb{Z}_2\) fall into two categories: \(w^2=\frac{1}{4}\) mod 1 and \(w^2=\frac{3}{4}\) mod 1. Substituting the values, we find \[\begin{align} (m_1+m_3)(\mathbf{1},\mathbf{2})+(m_2+m_4)(\mathbf{2},\mathbf{1}) \, , \end{align}\] for both classes of shift vectors.
Similarly to what happens in the \(g\)-twisted sector we have \[\begin{align} (r+v_{1f}+2v_{2f})^2=\frac{1}{2}, N_L=0=N_R, (p+v+2w)^2=\frac{3}{2}\, ; \nonumber \\ (r+v_{1f}+2v_{2f})^2=\frac{1}{2},N_L=0, (p+v+2w)^2=\frac{1}{2}=N_R \, . \end{align}\] Again the second class of solutions is projected out and the first one is non-trivial, leading to \[\begin{align} \frac{16}{8}\sum_{l=0}^3\sum_{d=0}^1\Phi^{2,l} \genfrac{[}{]}{0pt}{}{1}{d} =\begin{cases} 4-4e^{-2 \pi \mathrm{i}p\cdot v} \quad \text{NS}\\ 4+4e^{-2 \pi \mathrm{i}p\cdot v} \quad \text{R} \end{cases} \, , \end{align}\] where we have taken into account the fixed point degeneracy of 16. Supposing there are \(n_1\) solutions to \((p+v+2w)^2=\frac{3}{2}\) such that \(p\cdot2v\) is odd and \(n_2\) such that it is even, we get \[\begin{align} 16 n_1(\mathbf{1},\mathbf{1}) + 8n_2(\mathbf{2},\mathbf{1}) \, . \end{align}\] As expected, since \(h^2\) acts trivially on bosons, the NS spectrum of this twisted sector is the same as the one in the \(g\)-sector, while states in the R sector have opposite chirality.
We are grateful to Ida Zadeh for helpful discussions during the early stages of this work and for useful comments on an earlier version of the manuscript. We also thank Yuta Hamada, Hiroki Wada, and Masashi Kawahira for useful comments and discussions. A.I. also thanks Syun’ya Mizoguchi for his encouragement throughout this work.
The theta functions are
\[\require{physics} \label{eq:thetas} \begin{align} \vartheta_1=\vartheta\Bigl[\begin{matrix}1\\1\end{matrix}\Bigr]\mathrel{\vcenter{:}}= &i \sum_{n \in \mathbb{Z}}(-1)^n q^{\frac{1}{2}\qty(n-\frac{1}{2})^2}=0, \\ \vartheta_2=\vartheta\Bigl[\begin{matrix}1\\0\end{matrix}\Bigr]\mathrel{\vcenter{:}}=&\sum_{n \in \mathbb{Z}} q^{\frac{1}{2}\qty(n-\frac{1}{2})^2} =2 q^{\frac{1}{8}} \prod_{n=1}^{\infty}\qty(1-q^n)\qty(1+ q^n)\qty(1+ q^n),\\ \vartheta_3=\vartheta\Bigl[\begin{matrix}0\\0\end{matrix}\Bigr]\mathrel{\vcenter{:}}=&\sum_{n \in \mathbb{Z}} q^{\frac{1}{2} n^2} =\prod_{n=1}^{\infty}\qty(1-q^n)\qty(1+ q^{n-\frac{1}{2}})\qty(1+ q^{n-\frac{1}{2}}),\\ \vartheta_4=\vartheta\Bigl[\begin{matrix}0\\1\end{matrix}\Bigr]\mathrel{\vcenter{:}}=&\sum_{n \in \mathbb{Z}}(-1)^n q^{\frac{1}{2} n^2}=\prod_{n=1}^{\infty}\qty(1-q^n)\qty(1- q^{n-\frac{1}{2}})\qty(1-q^{n-\frac{1}{2}}). \end{align}\tag{16}\] A more general definition is given by
\[\require{physics} \vartheta\genfrac{[}{]}{0pt}{}{a}{b}=\sum_{n\in\mathbb{Z}}q^{\frac{1}{2}\qty(n+\frac{a}{2})^2}e^{\pi i b\qty(n+\frac{a}{2})}, a,b\in\mathbb{Z}.\] Sometimes we will use the short-hand notation, e.g. \(\vartheta_{10}\) for \(\vartheta \genfrac{[}{]}{0pt}{}{1}{0}\).
We used the following equations elsewhere:
\[\require{physics} \label{eq:theta95priodicity} \begin{align} \vartheta \genfrac{[}{]}{0pt}{}{a\pm 2}{b}&=\vartheta \genfrac{[}{]}{0pt}{}{a}{b},&&\vartheta \genfrac{[}{]}{0pt}{}{a}{b\pm 2}=e^{\pm\pi i a}\vartheta \genfrac{[}{]}{0pt}{}{a}{b},\\ T\cdot\vartheta \genfrac{[}{]}{0pt}{}{a}{b}=&e^{-\frac{1}{4}\pi i a(a-2)}\vartheta \genfrac{[}{]}{0pt}{}{a}{a+b-1},&&S\cdot\vartheta \genfrac{[}{]}{0pt}{}{a}{b}=\qty(-i\tau)^{\frac{1}{2}}e^{\frac{1}{2}\pi i ab}\vartheta \genfrac{[}{]}{0pt}{}{b}{-a},\\ \end{align}\tag{17}\]
\[\prod_{n=1}^\infty\frac{1}{(1+q^n)^4}=\frac{4\eta^6}{\vartheta_{10}^2}\prod_{n=1}\frac{1}{(1-q^n)^4}.\]
The definition of the Dedekind eta function and its transformation are \[\require{physics} \begin{align} &\eta(\tau)=q^{\frac{1}{24}}\prod_{n=1}^\infty \qty(1-q^n),\\ \eta(\tau+1)&=e^{\frac{\pi i}{12}}\eta(\tau), \quad\eta\qty(-\frac{1}{\tau})=\sqrt{-i\tau}\eta(\tau). \end{align} \label{eq:eta95modular95tr}\tag{18}\]
In this section we construct \(g,h,gh\)-twisted Hilbert spaces for \(\mathcal{H}_{\Gamma_{4,4}}\) explicitly, and calculate \[\begin{align} \Gamma^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}=\eta^4\bar\eta^4\tr_{\mathcal{H}_{\Gamma_{4,4}}^{g^c h^k}}g^d h^lq^{L_0}\bar q^{\bar{L}_0} \end{align}\] as its trace. On the compactification of heterotic strings on \(T^4\), there is a Hilbert space associated with \(\Gamma_{4,4}=H^1(T^4;\mathbb{Z})\oplus H_1(T^4;\mathbb{Z})\): \[\require{physics} \mathcal{H}_{\Gamma_{4,4}}=\mathcal{H}_{\text{Boson}}^{4,4}\otimes\qty(\bigoplus_{p\in\Gamma_{4,4}}\mathbb{C}\ket{p}).\] The actions of \(g,h\) on the fields are \[\begin{align} g: \begin{cases} X^6 \rightarrow -X^6,\\ X^7 \rightarrow -X^7 ,\\ X^8 \rightarrow -X^8 ,\\ X^9 \rightarrow -X^9, \\ \end{cases} \quad \quad\quad h : \begin{cases} X^6 \rightarrow -X^6,\\ X^7 \rightarrow -X^7+\pi R_7 ,\\ X^8 \rightarrow X^8 ,\\ X^9 \rightarrow X^9+ \pi R_9. \\ \end{cases} \end{align}\] The action of \(g,h\) on a state \(\ket{p}\) are given as follows: \[\begin{align} g\ket{p_6,p_7,p_8,p_9}=&\ket{-p_6,-p_7,-p_8,-p_9},\\ h\ket{p_6,p_7,p_8,p_9}=&(-1)^{R_7p_7+R_9 p_9}\ket{-p_6,-p_7,p_8,p_9},\\ gh\ket{p_6,p_7,p_8,p_9}=&(-1)^{R_7p_7+R_9 p_9}\ket{p_6,p_7,-p_8,-p_9}.\\ \end{align}\]
Now we can compute \(\Gamma^{0,l} \genfrac{[}{]}{0pt}{}{0}{d}\) as follows:
\[\begin{align} &\Gamma_{4,4}^{0,0} \genfrac{[}{]}{0pt}{}{0}{0}=\sum_{p\in\Gamma_{4,4}}q^{\frac{1}{2}p_L^2}\bar q^{\frac{1}{2}p_R^2}, &&\Gamma_{4,4}^{0,0} \genfrac{[}{]}{0pt}{}{0}{1}=\frac{16|\eta|^{12}}{|\vartheta_{10}|^4},\\ &\Gamma^{0,1}_{4,4} \genfrac{[}{]}{0pt}{}{0}{0} =\biggl|\frac{2\eta^3}{\vartheta_{10}}\biggr|^2\vartheta(R_8)\vartheta^{0,1}(R_9), &&\Gamma^{0,1}_{4,4} \genfrac{[}{]}{0pt}{}{0}{1} =\biggl|\frac{2\eta^3}{\vartheta_{10}}\biggr|^2\vartheta(R_6)\vartheta^{0,1}(R_7), \end{align}\] where \(\Gamma_{1,1}(R)\) is a \((1,1)\) even self-dual lattice and \(\vartheta(R),\vartheta^{k,l}(R)\) are its (modified) theta function:
\[\require{physics} \begin{align} \Gamma_{1,1}(R)=&\frac{1}{\sqrt{2}}\Bigl\{\left(\frac{n}{R}+mR;\frac{n}{R}-mR\right)|n,m\in\mathbb{Z}\Bigr\},\\ \vartheta(R)=&\sum_{p\in\Gamma_{1,1}(R)}q^{\frac{1}{2}p_L^2}\bar q^{\frac{1}{2}p_R^2},\\ \theta^{k,l}(R)=&\sum_{n\in\mathbb{Z}}\sum_{m\in \mathbb{Z}+\frac{1}{2}k}(-1)^{ln}q^{\frac{1}{4}\qty(\frac{n}{R}+mR)^2}\bar q^{\frac{1}{4}\qty(\frac{n}{R}-mR)^2}. \end{align}\]
The boundary conditions of \(g\)-twisted bosons \(X^i_g:T^2\to T^4\) are
\[\begin{align} X_g^{i}(\tau,\sigma+2\pi)=&g\cdot X_g^i(\tau,\sigma)\\ =&-X_g^i(\tau,\sigma),\quad6\leq i,j\leq 9. \end{align}\]
Let us denote the Fock space of \(X_g^i\) by \(\mathcal{H}^{4,4}_{X_g}\). The \(g\)-twisted Hilbert space \(\mathcal{H}_{\Gamma_{4,4}}^{g}\) is given concretely as follows: \[\require{physics} \begin{align} \mathcal{H}_{\Gamma_{4,4}}^{g}=&\mathcal{H}^{4,4}_{X_g}\otimes\qty(\bigoplus_{\text{fixed points}}\mathbb{C}\ket{x_6,x_7,x_8,x_9})\\ \end{align}\] where all 16 states \(\ket{x_6,x_7,x_8,x_9}\) are \(g\) invariant, while half pairs of them are \(h\) invariant. Canonical quantization gives
\[\begin{align} [\alpha^i_r,\alpha^j_s]=&[\tilde{\alpha}^i_r,\tilde{\alpha}^j_s]=r\delta^{ij}\delta_{r+s,0},\\ [\alpha^i_r,\tilde{\alpha}^j_s]=&0, \quad r,s\in\mathbb{Z}+\frac{1}{2},\quad6\leq i,j\leq 9, \end{align}\] and \[\require{physics} \begin{align} L_0=&\frac{1}{2}\sum_{i=6}^9\sum_{r=1/2}\qty(\alpha_{-r}^i\alpha_{r}^i+\alpha_r^i\alpha_{-r}^i)\\ =&\sum_{i=6}^9\sum_{r=1/2}\alpha_{-r}^i\alpha_{r}^i+\frac{1}{12}. \end{align}\]
Then it holds that
\[\label{eq:gamma01} \begin{align} &\Gamma_{4,4}^{0,0} \genfrac{[}{]}{0pt}{}{1}{0}=16\biggl|\frac{\eta^{12}}{\vartheta_{01}^4}\biggr|, &&\Gamma_{4,4}^{0,0} \genfrac{[}{]}{0pt}{}{1}{1}=16\biggl|\frac{\eta^{12}}{\vartheta_{00}^4}\biggr|,\\ & \Gamma_{4,4}^{0,1} \genfrac{[}{]}{0pt}{}{1}{0}=\Gamma_{4,4}^{0,1} \genfrac{[}{]}{0pt}{}{1}{1}= 0. \end{align}\tag{19}\]
The boundary conditions of \(h\)-twisted bosons \(X^i_h:T^2\to T^4\) are given as follows: \[\begin{align} X_h^{i}(\sigma_1+2\pi,\sigma_2)=&h\cdot X_h^i(\sigma_1,\sigma_2),i=6,7,8,9, \end{align}\] where \[h : \begin{cases} X^6 \rightarrow -X^6,\\ X^7 \rightarrow -X^7+\pi R_7 ,\\ X^8 \rightarrow X^8,\\ X^9 \rightarrow X^9+\pi R_9. \\ \end{cases}\]
Let us denote the Fock space of \(X_h^i\) by \(\mathcal{H}^{4,4}_{X_h}\). The \(h\)-twisted Hilbert space \(\mathcal{H}_{\Gamma_{4,4}}^{h}\) is given concretely as follows:
\[\begin{align} \mathcal{H}_{\Gamma_{4,4}}^{h}=& \mathcal{H}^{4,4}_{X_h}\otimes \left(\bigoplus_{\substack{x^6=0,\pi R_6 \\ x^7=\pm\frac{1}{2}\pi R_7\\(p_8,p_9)\in\Gamma_{2,2}(R_8,R_9)}}\mathbb{C}\ket{x_6,x_7,p_8,p_9}\right).\\ \end{align}\]
The actions of \(g,h\) are given by
\[\begin{align} g\ket{x_6,x_7,p_8,p_9}=&\ket{-x_6,-x_7,-p_8,-p_9},\\ h\ket{x_6,x_7,p_8,p_9}=&(-1)^{R_9p_9}\ket{-x_6,-x_7+\pi R_7,p_8,p_9},\\ gh\ket{x_6,x_7,p_8,p_9}=&(-1)^{R_9p_9}\ket{x_6,x_7+\pi R_7,-p_8,-p_9}. \end{align}\] The winding number \(m_9\) of \(X^9\) takes its value in \(\mathbb{Z}+\frac{1}{2}\) since it holds that \[\begin{align} \pi R_9=&X^9(2\pi)-X^9(0)\\ =&2\pi m_9 R_9. \end{align}\] Canonical quantization gives
\[\begin{align} [\alpha^i_r,\alpha^j_s]=&[\tilde{\alpha}^i_r,\tilde{\alpha}^j_s]=r\delta^{ij}\delta_{r+s,0},\\ [\alpha^i_r,\tilde{\alpha}^j_s]=&0, \quad r,s\in\mathbb{Z}+\frac{1}{2},i,j=6,7. \end{align}\]
\[\begin{align} [\alpha^i_n,\alpha^j_m]=&[\tilde{\alpha}^i_n,\tilde{\alpha}^j_m]=n\delta^{ij}\delta_{n+m,0},\\ [\alpha^i_n,\tilde{\alpha}^j_m]=&0, \quad n,m\in\mathbb{Z},i,j=8,9. \end{align}\]
\(L_0\) is given by
\[\require{physics} \begin{align} L_0=&\frac{1}{2}\sum_{i=6}^7\sum_{r=1/2}\qty(\alpha_{-r}^i\alpha_{r}^i+\alpha_r^i\alpha_{-r}^i)+\frac{1}{2}\sum_{i=8}^9\sum_{n=1}\qty(\alpha_{-n}^i\alpha_{n}^i+\alpha_n^i\alpha_{-n}^i)\\ =&\sum_{i=6}^7\sum_{r=1/2}\alpha_{-r}^i\alpha_{r}^i+\sum_{i=8}^9\sum_{n=1}\alpha_{-n}^i\alpha_{n}^i-\frac{1}{24}. \end{align}\] The result is \[\begin{align} \Gamma_{4,4}^{1,0} \genfrac{[}{]}{0pt}{}{0}{0} =&4\frac{\eta^3\bar\eta^3}{|\vartheta_{01}|^2}\vartheta(R_8)\vartheta^{1,0}(R_9),\\ \Gamma_{4,4}^{1,0} \genfrac{[}{]}{0pt}{}{0}{1}=& \Gamma_{4,4}^{1,1} \genfrac{[}{]}{0pt}{}{0}{1}=0,\\ \Gamma_{4,4}^{1,1} \genfrac{[}{]}{0pt}{}{0}{0} =&4\frac{\eta^3\bar\eta^3}{|\vartheta_{00}|^2}\vartheta(R_8)\vartheta^{1,1}(R_9),\\ \end{align}\]
The boundary conditions of \(gh\)-twisted bosons \(X^i_{gh}:T^2\to T^4\) are
\[\begin{align} &X_{gh}^{i}(\sigma_1+2\pi,\sigma_2)\\=&gh\cdot X_{gh}^i(\sigma_1,\sigma_2),~~i=6,7,8,9, \end{align}\]
where
\[gh : \begin{cases} X^6 \rightarrow X^6,\\ X^7 \rightarrow X^7 +\pi R_7 ,\\ X^8 \rightarrow -X^8,\\ X^9 \rightarrow -X^9- \pi R_9. \\ \end{cases}\]
The winding number \(m_7\) of \(X^7\) takes its value in \(\mathbb{Z}+\frac{1}{2}\) in the \(gh\)-twisted sector: \[\begin{align} \pi R_7=&X^7(2\pi)-X^7(0)\\ =&2\pi m_7 R_7. \end{align}\]
After the canonical quantization, we find
\[\begin{align} [\alpha^i_n,\alpha^j_m]=&[\tilde{\alpha}^i_n,\tilde{\alpha}^j_m]=n\delta^{ij}\delta_{n+m,0},\\ [\alpha^i_n,\tilde{\alpha}^j_m]=&0, \quad n,m\in\mathbb{Z},i,j=6,7, \end{align}\]
\[\begin{align} [\alpha^i_r,\alpha^j_s]=&[\tilde{\alpha}^i_r,\tilde{\alpha}^j_s]=r\delta^{ij}\delta_{r+s,0},\\ [\alpha^i_r,\tilde{\alpha}^j_s]=&0, \quad r,s\in\mathbb{Z}+\frac{1}{2},i,j=8,9. \end{align}\]
and \[\require{physics} \begin{align} L_0=&\frac{1}{2}\sum_{i=6}^7\sum_{n=1}\qty(\alpha_{-n}^i\alpha_{n}^i+\alpha_n^i\alpha_{-n}^i)+\frac{1}{2}\sum_{i=8}^9\sum_{r=1/2}\qty(\alpha_{-r}^i\alpha_{r}^i+\alpha_r^i\alpha_{-r}^i)\\ =&\sum_{i=6}^7\sum_{n=1}\alpha_{-n}^i\alpha_{n}^i+\sum_{i=8}^9\sum_{r=1/2}\alpha_{-r}^i\alpha_{r}^i-\frac{1}{24}. \end{align}\]
Let us denote the Fock space of \(X_{gh}^i\) by \(\mathcal{H}^{4,4}_{X_{gh}}\). The \(gh\)-twisted Hilbert space \(\mathcal{H}_{\Gamma_{4,4}}^{gh}\) is given concretely as follows:
\[\begin{align} \mathcal{H}_{\Gamma_{4,4}}^{gh}=& \mathcal{H}^{4,4}_{X_{gh}}\otimes \left(\bigoplus_{\substack{ x^8=0,\pi R_8\\x^9=\pm\frac{1}{2}\pi R_9 \\(p_6,p_7)\in\Gamma_{2,2}(R_6,R_7)}}\mathbb{C}\ket{p_6,p_7,x_8,x_9}\right).\\ \end{align}\] The actions of \(g,h,gh\) are defined as
\[\begin{align} g\ket{p_6,p_7,x_8,x_9}=&\ket{-p_6,-p_7,-x_8,-x_9},\\ h\ket{p_6,p_7,x_8,x_9}=&(-1)^{R_7 p_7}\ket{-p_6,-p_7,x_8,x_9+\pi R_9},\\ gh\ket{p_6,p_7,x_8,x_9}=&(-1)^{R_7 p_7}\ket{p_6,p_7,-x_8,-x_9+\pi R_9}.\\ \end{align}\] Then the results are
\[\begin{align} \Gamma_{4,4}^{1,0} \genfrac{[}{]}{0pt}{}{1}{0} =&4\frac{\eta^3\bar\eta^3}{|\vartheta_{01}|^2}\vartheta(R_6)\vartheta^{1,0}(R_7),\\ \Gamma_{4,4}^{1,0} \genfrac{[}{]}{0pt}{}{1}{1}=&\Gamma_{4,4}^{1,1} \genfrac{[}{]}{0pt}{}{1}{0}=0,\\ \Gamma_{4,4}^{1,1} \genfrac{[}{]}{0pt}{}{1}{1} =&4\frac{\eta^3\bar\eta^3}{|\vartheta_{00}|^2}\vartheta(R_6)\vartheta^{1,1}(R_7). \end{align}\]
In this section we obtain the following expression: \[\begin{align} \Gamma^{0,0} \genfrac{[}{]}{0pt}{}{c}{d}=&\frac{16|\eta|^{12}}{|\vartheta \genfrac{[}{]}{0pt}{}{1+c}{1+d}|^4},(c,d)\neq (0,0)\\ \Gamma^{k,l} \genfrac{[}{]}{0pt}{}{0}{0}=&\frac{4\eta^3\bar\eta^3}{|\vartheta \genfrac{[}{]}{0pt}{}{k+1}{l+1}|^2}\theta(R_8)\theta^{k,l}(R_9),\\ \Gamma^{c,d} \genfrac{[}{]}{0pt}{}{c}{d}=&\frac{4\eta^3\bar\eta^3}{|\vartheta \genfrac{[}{]}{0pt}{}{1+c}{1+d}|^2}\theta(R_6)\theta(R_7),(c,d)\neq (0,0).\\ \end{align}\] Others are zero: \[\begin{align} \Gamma^{0,1} \genfrac{[}{]}{0pt}{}{1}{d}=&0,\quad\text{(g-twisted sector)}\\ \Gamma^{1,l} \genfrac{[}{]}{0pt}{}{0}{1}=&0,\quad\text{(h-twisted sector)}\\ \Gamma^{1,d+1} \genfrac{[}{]}{0pt}{}{1}{d}=&0,\quad\text{(gh-twisted sector)}\\ \end{align}\]
In this appendix, we summarise features of the root and weight lattices used in this paper. For notation and more details see [64]. We denote the \(i\)-th standard orthonormal basis vector of \(\mathbb{R}^n\) by \(\varepsilon_i\).
The root system of \(A_n\) is \[\begin{align} &\varepsilon_j-\varepsilon_k, &&\text{for}\quad j\neq k \end{align}\] The basis is \[\begin{align} &\alpha_i=\varepsilon_i-\varepsilon_{i+1},&&\text{for}\quad 1\leq i \leq n. \end{align}\]
The root system of \(D_n\): \[\begin{align} &\pm\varepsilon_j\pm\varepsilon_k,\pm\varepsilon_j\mp\varepsilon_k,&&\text{for}\quad 1\leq j<k\leq n. \end{align}\] The basis is \[\begin{align} &\alpha_i=\varepsilon_i-\varepsilon_{i+1},&&\text{for}\quad 1\leq i\leq n-1, \\ &\alpha_n=\varepsilon_{n-1}+\varepsilon_n. \end{align}\]
The root system of \(E_6\) is \[\require{physics} \begin{align} &\pm\varepsilon_i\pm\varepsilon_j, &&\text{for}\quad1\leq i<j\leq 5,\\ &\pm\frac{1}{2}\qty(\varepsilon_8-\varepsilon_7-\varepsilon_6+\sum_{i=1}^5 (-1)^{\nu_i}\varepsilon_i), &&\text{for}\quad\sum_{i=1}^5 \nu_i\in 2\mathbb{Z}. \end{align}\] The basis is \[\begin{align} &\alpha_1=\frac{1}{2}\left(\varepsilon_1+\varepsilon_8\right)-\frac{1}{2}\left(\varepsilon_2+\varepsilon_3+\varepsilon_4+\varepsilon_5+\varepsilon_6+\varepsilon_7\right), &&\alpha_2=\varepsilon_1+\varepsilon_2, \\ & \alpha_3=\varepsilon_2-\varepsilon_1, &&\alpha_4=\varepsilon_3-\varepsilon_2, \\ &\alpha_5=\varepsilon_4-\varepsilon_3, &&\alpha_6=\varepsilon_5-\varepsilon_4. \end{align}\]
The root system of \(E_7\) is \[\require{physics} \begin{align} &\pm\varepsilon_i\pm\varepsilon_j, &&\text{for}\quad1\leq i<j\leq 6, \\ &\pm(\varepsilon_7-\varepsilon_8),\\ &\pm\frac{1}{2}\qty(\varepsilon_8-\varepsilon_7-\varepsilon_6+\sum_{i=1}^6 (-1)^{\nu_i}\varepsilon_i), &&\text{for}\quad\sum_{i=1}^6 \nu_i\in 2\mathbb{Z}. \end{align}\] The basis is \[\begin{align} & \alpha_1=\frac{1}{2}\left(\varepsilon_1+\varepsilon_8\right)-\frac{1}{2}\left(\varepsilon_2+\varepsilon_3+\varepsilon_4+\varepsilon_5+\varepsilon_6+\varepsilon_7\right), \\ & \alpha_2=\varepsilon_1+\varepsilon_2, \quad\quad\quad\quad\alpha_3=\varepsilon_2-\varepsilon_1, \\ & \alpha_4=\varepsilon_3-\varepsilon_2, \quad\quad\quad\quad \alpha_5=\varepsilon_4-\varepsilon_3, \\ & \alpha_6=\varepsilon_5-\varepsilon_4, \quad\quad\quad\quad \alpha_7=\varepsilon_6-\varepsilon_5 . \end{align}\]
The root basis vectors of \(E_8\) type are \[\label{eq:32e8roots} \begin{align} &\alpha_1=\frac{1}{2}\left(\varepsilon_1+\varepsilon_8\right)-\frac{1}{2}\left(\varepsilon_2+\varepsilon_3+\varepsilon_4+\varepsilon_5+\varepsilon_6+\varepsilon_7\right), \\ &\alpha_2=\varepsilon_1+\varepsilon_2, \quad\quad\quad\quad \alpha_3=\varepsilon_2-\varepsilon_1,\\ &\alpha_4=\varepsilon_3-\varepsilon_2, \quad\quad\quad\quad \alpha_5=\varepsilon_4-\varepsilon_3 \\ &\alpha_6=\varepsilon_5-\varepsilon_4, \quad\quad\quad\quad \alpha_7=\varepsilon_6-\varepsilon_5,\\ &\alpha_8=\varepsilon_7-\varepsilon_6 . \end{align}\tag{20}\]
In this appendix we classify inequivalent shift vectors \(w\) for \(E_8 \times E_8\) and \(\mathrm{SemiSpin}(32)\) while fixing \(v=\frac{1}{2}(0^{14},1,1)\). Let us first state what do we mean by equivalent.
Modular invariance requires order 4 shifts \(w \in \frac{1}{4} \Gamma_{0,16}\) altogether with \(4w^2 \in 2\mathbb{Z}+1\) and \(4w \cdot v \in 4\mathbb{Z}+1\). Two shifts \(w,w'\) are equivalent if they differ by the following transformations: \[\begin{align} w^\prime=s(w)+\alpha \, , \quad \alpha \in \Gamma_{0,16} \, \end{align}\] where \(s\) is an element of Weyl subgroup of gauge algebra which stabilises \(v\). We can check the equivalence of two shift vectors explicitly as follows. Let \(A=cv+kw, B=dv+lw\), then the part of the partition function which comes from the lattice \(\Gamma_{0,16}\) is \[\bar\Gamma_{16,w}^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}=e^{\pi \mathrm{i}A\cdot B}\sum_{p\in\Gamma_{0,16}}\bar q^{\frac{1}{2}(p+A)^2}e^{-2\pi \mathrm{i}(p+A)\cdot B} \, .\] Replacing \(w\) by \(w^\prime=s(w)+\alpha\) yields \[\begin{align} \bar\Gamma_{16,w^\prime}^{k,l}=&e^{\pi \mathrm{i}(s(A)+k\alpha)\cdot (s(B)+l\alpha)}\sum_{p\in\Gamma_{0,16}}\bar q^{\frac{1}{2}(p+s(A)+k\alpha)^2}e^{-2\pi i(p+s(A)+k\alpha)\cdot (s(B)+l\alpha)}\\ =&e^{\pi \mathrm{i}(dk-lc)v\cdot\alpha}\bar\Gamma_{16,w}^{k,l} \, , \end{align}\] which leads to \(\bar\Gamma_{16,w+\alpha}^{k,l} \equiv\bar\Gamma_{16,w}^{k,l}\) up to discrete torsion \(e^{ \pi \mathrm{i}(dk-lc)\mathbb{Z}}\) [57]. Here \(v \cdot \alpha \in \mathbb{Z}\) follows from the mixed condition on the shift vectors. Here we used the fact that a Weyl group action preserve inner product: \(s(A)\cdot s(B)=A\cdot B\), for example.
We will classify inequivalent shift vectors on \(E_8 \times E_8\) using Kac’s algorithm (see [65] and references therein). Given the structure of this algebra, the discussion can be carried out for each \(E_8\) factor separately to combine the admissible pairs at the end. For \(\mathrm{Spin}(32)/\mathbb{Z}_2\) it will be easier to use a different method instead. Recalling the definition of the \(\mathrm{SemiSpin}(32)\) lattice \[\require{physics} \Gamma_{\mathrm{Spin}(32)/\mathbb{Z}_2}=\left\{p\in\mathbb{Z}^{16}\text{ or }\qty(\mathbb{Z}+\tfrac{1}{2})^{16}\middle|\sum_{i=1}^{16}p_i \in2\mathbb{Z}\right\} \, ,\] we will take separately into account vectors of the form \(4w \in \mathbb{Z}^{16}\) and \(4w \in \left( \mathbb{Z}+\frac{1}{2} \right)^{16}\).
We describe the \(E_8\times E_8\) shift classes by a convenient change of simple-root basis. The point is to enlarge the list of simple-root candidates by adding the negative of the highest root.
Let \(\alpha_1,\ldots,\alpha_8\) be the simple roots of \(E_8\), with the conventions used above. We define \(\alpha_9\) by \[\alpha_9 = -\left( 2\alpha_1+3\alpha_2+4\alpha_3+6\alpha_4 +5\alpha_5+4\alpha_6+3\alpha_7+2\alpha_8 \right).\] Thus \(-\alpha_9\) is the highest root of \(E_8\). Equivalently, \(-\alpha_9\) is the positive root such that, for every positive root \(\rho\), the difference \[-\alpha_9-\rho\] is a non-negative integral linear combination of \(\alpha_1,\ldots,\alpha_8\).
Of course, \(\alpha_9,\alpha_1,\ldots,\alpha_8\) are not independent, but they satisfy \[\alpha_9 +2\alpha_1+3\alpha_2+4\alpha_3+6\alpha_4 +5\alpha_5+4\alpha_6+3\alpha_7+2\alpha_8 =0.\] We draw a graph whose nodes are the roots \(\alpha_i\), with an edge between two nodes when the corresponding roots have inner product \(-1\). In our convention the graph is \[\begin{array}{ccccccccccccccc} \alpha_1 & - & \alpha_3 & - & \alpha_4 & - & \alpha_5 & - & \alpha_6 & - & \alpha_7 & - & \alpha_8 & - & \alpha_9 \\[-1mm] &&&& | \\ &&&& \alpha_2 . \end{array}\] This is called the extended Dynkin diagram of \(E_8\).
Let \(w\) be a \(\mathbb{Z}_4\) shift in one \(E_8\) factor. Here \(\alpha_i\cdot w\) means the ordinary Euclidean inner product of the root \(\alpha_i\) with the shift vector \(w\). Up to Weyl reflections and lattice shifts, we choose \(w\) so that \[\alpha_i\cdot w\ge 0 \qquad (i=1,\ldots,8), \qquad -\alpha_9\cdot w\le 1.\] If \(-\alpha_9\cdot w>1\), then an equivalent shift vector \(w_n^\prime=w+n\alpha_9\) satisfies \(-\alpha_9\cdot w^\prime=-\alpha_9\cdot w-2n<1\) for some integer \(n\). We then attach non-negative integers to the nodes by \[s_i=4\,\alpha_i\cdot w \qquad (i=1,\ldots,8), \qquad s_9=4(1+\alpha_9\cdot w).\] We record the shift class as \[[s_1s_2\cdots s_8;s_9].\] The relation among the roots implies \[2s_1+3s_2+4s_3+6s_4+5s_5+4s_6+3s_7+2s_8+s_9=4.\]
The unbroken semisimple root system is obtained by keeping precisely the nodes with label \(0\). In other words, the zero-labelled part of the graph is read as the Dynkin diagram of the unbroken semisimple algebra.
For the first \(E_8\) factor, there are ten shift classes are summarised in table ¿tbl:tab:32first32comp32E8?.
\[\begin{array}{c|c|c|c|c} \text{class} & [s_1\cdots s_8;s_9] & \text{zero-labelled subgraph} & 4\mathsf A_i^2 \bmod 2 \\ \hline \mathsf A_1 & [00000000;4] & E_8 & 0 \\ \mathsf A_2 & [00000002;0] & E_7\oplus A_1 & 0 \\ \mathsf A_3 & [20000000;0] & D_8 & 0 \\ \mathsf A_4 & [10000000;2] & D_7 & 1 \\ \mathsf A_5 & [00000100;0] & D_5\oplus A_3 & 1 \\ \mathsf A_6 & [01000000;1] & A_7 & 0 \\ \mathsf A_7 & [00000001;2] & E_7 & 1/2 \\ \mathsf A_8 & [10000001;0] & D_6\oplus A_1 & 1/2 \\ \mathsf A_9 & [00000010;1] & E_6\oplus A_1 & 3/2 \\ \mathsf A_{10} & [00100000;0] & A_7\oplus A_1 & 3/2 \end{array}\]
For the second \(E_8\) factor, the fixed \(\mathbb{Z}_2\) shift \(v\) first leaves an \(E_7\oplus A_1\) root system. The relevant nonabelian part is read from the \(E_7\) factor. In this paragraph, we again denote the simple roots of this \(E_7\) by \(\alpha_1,\ldots,\alpha_7\). We define \[\alpha_0 = -\left( 2\alpha_1+2\alpha_2+3\alpha_3+4\alpha_4 +3\alpha_5+2\alpha_6+\alpha_7 \right).\] Thus \(-\alpha_0\) is the highest root of this \(E_7\). The roots \(\alpha_0,\alpha_1,\ldots,\alpha_7\) satisfy \[\alpha_0 +2\alpha_1+2\alpha_2+3\alpha_3+4\alpha_4 +3\alpha_5+2\alpha_6+\alpha_7 =0.\] Their graph is \[\begin{array}{ccccccccccccc} \alpha_0 & - & \alpha_1 & - & \alpha_3 & - & \alpha_4 & - & \alpha_5 & - & \alpha_6 & - & \alpha_7 \\[-1mm] &&&&&& | \\ &&&&&& \alpha_2 . \end{array}\] This is called the extended Dynkin diagram of \(E_7\).
For an order-four shift in this \(E_7\) factor, we choose \(w\) so that \[\alpha_i\cdot w\ge0 \qquad (i=1,\ldots,7), \qquad -\alpha_0\cdot w\le1,\] up to lattice shift, and define the node labels by \[r_i=4\,\alpha_i\cdot w \qquad (i=1,\ldots,7), \qquad r_0=4(1+\alpha_0\cdot w).\] We record the class as \[[r_0;r_1r_2\cdots r_7].\] The labels obey \[r_0+2r_1+2r_2+3r_3+4r_4+3r_5+2r_6+r_7=4.\] Again, the zero-labelled part of the graph gives the Dynkin diagram of the unbroken semisimple algebra.
The ten \(\mathsf B_i\) classes are summarized in the table ¿tbl:tab:32second32comp32E8?
\[\begin{array}{c|c|c|c|c} \text{class} & [r_0;r_1\cdots r_7] & \text{zero-labelled subgraph} & 4\mathsf B_i^2 \bmod 2 \\ \hline \mathsf B_1 & [0;1000010] & D_4\oplus A_1\oplus A_1 & 1 \\ \mathsf B_2 & [0;0000012] & D_6 & 1 \\ \mathsf B_3 & [2;0000002] & E_6 & 0 \\ \mathsf B_4 & [0;0200000] & A_7 & 0 \\ \mathsf B_5 & [0;0000101] & A_5\oplus A_1 & 0 \\ \mathsf B_6 & [1;0100001] & A_5 & 1/2 \\ \mathsf B_7 & [0;0000020] & D_6\oplus A_1 & 1/2 \\ \mathsf B_8 & [4;0000000] & E_7 & 1/2 \\ \mathsf B_9 & [2;0000010] & D_5\oplus A_1 & 3/2 \\ \mathsf B_{10} & [0;0001000] & A_3\oplus A_3\oplus A_1 & 3/2 \end{array}\]
Finally, writing the full \(E_8\times E_8\) shift as \[w=(\mathsf A;\mathsf B),\] the condition \(4w^2\in 2\mathbb{Z}+1\) is equivalent to \[\label{eq:here1} \mathsf A^2+\mathsf B^2\equiv 1 \pmod 2.\tag{21}\] Hence the allowed pairs are read off from the last columns of the two tables: \[\begin{align} &\{\mathsf A_1,\mathsf A_2,\mathsf A_3,\mathsf A_6\} \times \{\mathsf B_1,\mathsf B_2\}, \\ &\{\mathsf A_4,\mathsf A_5\} \times \{\mathsf B_3,\mathsf B_4,\mathsf B_5\}, \\ &\{\mathsf A_7,\mathsf A_8\} \times \{\mathsf B_6,\mathsf B_7,\mathsf B_8\}, \\ &\{\mathsf A_9,\mathsf A_{10}\} \times \{\mathsf B_9,\mathsf B_{10}\}. \end{align}\] Then there are \(24\) inequivalent \(E_8\times E_8\) shift classes.
In this subsection let us assume that \(4w\in\mathbb{Z}^{16}\). In this case the components of \(w\) in terms of standard basis on \(\mathbb{R}^{16}\), which include \(\mathrm{Spin}(32)/\mathbb{Z}_2\) lattice, would be \[\pm\frac{1}{4}, \pm\frac{2}{4},\pm\frac{3}{4},\pm 1\] up to lattice shifts. Since \(w\) should satisfy \(4v\cdot w\in 4\mathbb{Z}+1\), the last two components of \(w\) should be one of the following: \[\require{physics} \qty(\frac{1}{4},\frac{1}{4}),\qty(0,\frac{1}{2})\] up to lattice shifts. For instance, one might take them to be \(\require{physics} \qty(\frac{3}{4},-\frac{1}{4})\), but these can be shifted to \(\require{physics} \qty(\frac{1}{4},\frac{1}{4})\) by \(\require{physics} \frac{1}{2}\qty(\cdots,-1,1)\in \Gamma_{\mathrm{Spin}(32)/\mathbb{Z}_2}\).
Claim. One can take the first \(14\) components of \(w\) as follows: \[\require{physics} \qty(0^{n_0},\qty(\pm\frac{1}{4})^{n_1},\qty(\pm \frac{1}{2})^{n_2};)\] up to lattice shifts. What differs from the \(E_8\times E_8\) case is that if we keep this expression for a shift vector, the pair of numbers \((n_0,n_1,n_2)\) is invariant under lattice shifts and Weyl reflection.
Let us show this statement. We start from the form of the first \(14\) components: \[\require{physics} \qty(0^{n_0},\qty(\pm\frac{1}{4})^{n_1},\qty(\pm \frac{1}{2})^{n_2},\qty(\pm \frac{3}{4})^{n_3}).\] Here we can assume that \(n_3=0,1\), because if \(n_3\geq 2\), there is at least one pair of two \(\frac{3}{4}\)’s. Then they can be shifted to a pair of \(\frac{1}{4}^2\); e.g, \(\require{physics} \qty(\cdots,\frac{3}{4}^{2n})\to \qty(\cdots,\frac{1}{4}^{2n})\) by \((\cdots,0,1^{2n})\in \mathrm{Spin}(32)/\mathbb{Z}_2\) lattice.
Next, we assume that \(n_3=1\): \[\require{physics} \qty(0^{n_0},\qty(\pm\frac{1}{4})^{n_1},\qty(\pm \frac{1}{2})^{n_2},\qty(\pm \frac{3}{4})^{1}).\] because if \(n_3=0\), we obtain \(\require{physics} \qty(0^{n_0},\qty(\pm \frac{1}{4})^{n_1},\qty(\pm\frac{1}{2})^{n_2})\) and the proof was already finished. Now we can assume that \(n_0\geq 1\) or \(n_2 \geq 1\), because otherwise \(n_1=13\). Then the norm of this vector is \[\require{physics} \qty(\qty(\pm\frac{1}{4})^{13},\pm \frac{3}{4})^2=\frac{11}{8},\] and can never yield \(4w^2\in 2\mathbb{Z}+1\) with the norms of the last two components: \(\require{physics} \qty(\frac{1}{4},\frac{1}{4})^2=\frac{1}{8}\) or \(\require{physics} \qty(0,\frac{1}{2})^2=\frac{1}{4}\). Now proof is finished and remained thing is to identify possible pair of \((n_1,n_2)\) with \(n_0=14-n_1-n_2\) which satisfies \(4w^2\in2\mathbb{Z}+1\). The result is \[\require{physics} \begin{align} W_{0,\frac{1}{2}}(n_1,n_2)=& \qty(\frac{1}{4}\qty(0^{14-n_1-n_2},1^{n_1},2^{n_2}),0,\frac{1}{2}):\\ (n_1,n_2)\in\Bigl\{ (0,0),(0,2)&,(0,4),(0,6), (4,1),(4,3),(4,5), (8,0),(8,2),(12,1) \Bigr\},\\ \end{align}\] and \[\require{physics} \begin{align} W_{\frac{1}{4},\frac{1}{4}}(n_1,n_2)=&\qty(\frac{1}{4}\qty(0^{14-n_1-n_2},1^{n_1},2^{n_2}),\frac{1}{4},\frac{1}{4}):\\ (n_1,n_2)\in\Bigl\{ (2,0),(2,2)&,(2,4),(2,6), (6,1),(6,3),(10,0),(10,2) \Bigr\} \, . \end{align}\] and
\[\require{physics} \widetilde{W}_{\frac{1}{4},\frac{1}{4}}(n)=\qty(-\frac{1}{2},\frac{1}{2}^{11-8n},\frac{1}{4}^{2+8n},\frac{1}{4},\frac{1}{4}),\quad n=0,1.\\\]
In this subsection we take as \(\require{physics} 4w\in\qty(\mathbb{Z}+\frac{1}{2})^{16}\). By a similar argument, we obtain the following four shift vectors \[\require{physics} \begin{align} W_{\frac{1}{8},\frac{3}{8}}(n)=&\left(\qty(\frac{1}{8})^{11-4n},\qty(\frac{3}{8})^{3+4n},\frac{1}{8},\frac{3}{8}\right) ,\quad n=0,1,\\ W_{\frac{5}{8},-\frac{1}{8}}(n)=& \left( \left(\frac{1}{8}\right)^{13-4n}, \left(\frac{3}{8}\right)^{1+4n}, \frac{5}{8}, -\frac{1}{8} \right), \quad n=0,1.\\ \end{align}\]
Of course \(W_{\frac{1}{8},\frac{3}{8}}(2)\) and \(W_{\frac{5}{8},-\frac{1}{8}}(n), n=2,3\) are also allowed, but they are equivalent to some the above four vectors up to Weyl reflection and lattice shifts: \[\require{physics} \begin{align} W^\prime =-W+\qty((\pm\frac{1}{2})^{14},\frac{1}{2},-\frac{1}{2}). \end{align}\] Permutations of the components with even number of sign flip are allowed in the Weyl group of \(D_n\).
The \(\mathbb{Z}_4\) orbifold is part of a class of orbifolds constructed in [1] with an additional spectator \(S^1\). In their notation, this corresponds to the \((r,a,\delta)=(18,2,0)\) point in the Nikulin involution space. These orbifolds act as \(\mathbb{Z}_2\) on bosons and as \(\mathbb{Z}_4\) on fermions.↩︎
We can include an additional phase \(e^{-\pi i(dk-lc)N}\) for \(N=0,1\), in the partition function; this is the discrete torsion phase [57]. In the present case, however, it does not lead to a different spectrum, because the only modular orbit on which this phase is non-trivial has a vanishing \(\Gamma_{4,4}^{k,l} \genfrac{[}{]}{0pt}{}{c}{d}\) contribution.↩︎