text:
# 1 Introduction

## 1 Introduction

Supersymmetry (SUSY) is one of the most promising ideas to solve the hierarchy problem, if broken at a low scale $`M_{SUSY}\mathrm{TeV}`$. It plays an even more important role in a quantum theory of gravity, such as string theory, since it provides a classical and quantum stabilization mechanism for string/M-theory vacua. The construction of absolutely stable string models, realistic or not, where SUSY is broken by means of any symmetry breaking mechanism is a fundamental and outstanding problem in string theory. In order to try to shed some light into this crucial question, it is essential to analyze the form and nature of instabilities that affect the known non-SUSY string models.

One of the most interesting mechanisms of SUSY breaking and essentially the unique available for oriented closed strings at the string level, is the Scherk-Schwarz (SS) mechanism  in which SUSY is broken by twisting the periodicity conditions along some compact dimensions. This mechanism allows to build classically stable models without tachyons if the compactification radius along which the SS twist is performed is chosen to be large enough. String models with SS SUSY breaking mechanism have received some attention in both heterotic  and open strings . On the contrary, not much analysis of the quantum stability of these vacua have been performed, showing that the models evolve either to smaller radii reaching eventually the tachyonic regime , or they recover unbroken SUSY by a decompactification limit 111See  for a similar analysis performed in an M-theory context and  for a nice general analysis in non-SUSY heterotic models. See also  for an analysis of the stability of a certain class of non-SUSY 6D orientifolds..

In this paper, along the lines of , we compute the one-loop induced vacuum energy density (cosmological constant) as a function of the radius of the twisted direction, for a certain class of Type IIB orbifold and orientifold string models with SS SUSY breaking. This is obtained by computing the one-loop partition function on the relevant world-sheet surfaces: the torus for IIB orbifolds, and the torus, Klein bottle, annulus and Möbius strip for IIB orientifolds.

The first model we consider is a simple 9-dimensional (9D) $`𝐙_2`$ IIB orientifold, whose vacuum energy density has already been considered in . In agreement with , we show how the cosmological constant of this model crucially depends on the choice made for the $`𝐙_2`$ Chan-Paton twist matrix $`\gamma _g`$. Depending on $`\gamma _g`$, the model evolves either towards the tachyonic regime or to the decompactification limit (see figure 1).

Orbifolds and orientifolds on twisted Asymptotically Locally Euclidean (ALE) spaces of the form $`(\times S^1)/𝐙_N`$, where $`𝐙_N`$ ($`N`$ odd) acts simultaneously as a rotation on $``$ and a translation on $`S^1`$, are studied next. Such spaces upon reduction on $`S^1`$ give rise to Melvin backgrounds , and have recently received renewed interest. We show that for any $`N`$ odd, with or without open strings, the vacuum energy density $`\rho ^N(R)`$ is a negative monotonic function, reaching zero as $`R\mathrm{}`$. In the orbifold models we find that $`\rho ^N(R)>\rho ^M(R)`$ for $`N<M`$, for any value of $`R`$, and for the lower $`N,M=3,5,7,9`$ considered (see figure 2). In the orientifold models, with a suitable choice of twist matrices, we find the same behavior (see figure 3). In both cases, then, it is resonable to assume that $`\rho ^N(R)>\rho ^M(R)`$ for $`N<M`$, $`N,M`$. The perturbative quantum fate of these models is then clear. For any initial value of $`R`$ and $`N`$, $`R`$ will shrink until $`R=R_T`$, the critical radius where some twisted string mode becomes tachyonic. The dynamics is now governed by the classical tachyonic instability and, according to the analysis of , the model will undergo a phase transition towards another $`(\times S^1)/𝐙_N^{}`$ model, with $`N^{}<N`$, with $`R`$ increasing along the transition . The $`𝐙_N^{}`$ model will also undergo a phase transition and so on, until the twist vanishes and a flat SUSY space-time is recovered. Quite interestingly, we get exactly the same energy density $`\rho ^N(R)`$ for non-compact $`(\times S^1)/𝐙_N`$ or compact $`(T^2\times S^1)/𝐙_N`$ models222In the compact case, $`N`$ has to be odd and has to preserve the lattice structure of the torus., although the fate of the two twisted directions seem to crucially depend on whether they are compact or not .

The last model we consider is a 4D chiral orientifold model, discussed in . This is based on a $`T^6/(𝐙_6^{}\times 𝐙_2^{})`$ orbifold, where $`𝐙_6^{}`$ is the usual rotation defining a SUSY orientifold  and $`𝐙_2^{}=\sigma ()^F`$, with $`\sigma `$ an order two translation along one of the radii ($`R`$) of a torus and $`F`$ the space-time fermion number operator. The model has $`D9`$-branes, as well as $`D5`$ and $`\overline{D}5`$-branes. We first study $`\rho `$ as a function of the single modulus $`R`$, fixing the remaining 5 other compactified directions. For a proper choice of the Chan-Paton twist matrix associated to the $`𝐙_2^{}`$ element in the $`D9`$ sector, we find a minimum of $`\rho (R)`$ at a finite value $`R_0`$, where $`\rho (R_0)<0`$, providing in this way an interesting non-trivial stabilization mechanism for the SS direction $`R`$ (see figure 4). In order to decide whether this minimum is actually an absolute minimum or not, we also study the dependence of $`R`$ on some other compact directions, taken to be equal for simplicity. We are not able to give a definite answer to whether $`R_0`$ is an absolute minimum or not, but our analysis seems to suggest that for finite values of the remaining compact directions, the minimum is always present, but it could run-away to infinity or towards the tachyonic instability due to the dynamics along the other directions (see figure 5).

An interesting by-product of our studies is provided by understanding which string states contribute to $`\rho (R)`$. Generalizing the well-known technique of unfolding the fundamental domain of the torus , we show that the closed string contribution to $`\rho (R)`$, in both orbifold and orientifolds, can be analytically computed and is given by untwisted closed strings only, where no winding modes along the SS direction appear. As far as the $`R`$ dependence is concerned, the whole one-loop string partition function looks effectively like that of a purely quantum field theory. This provides a generalization of what is well-known to happen to strings at finite temperature  and to $`𝐙_2`$ SS twists . For open strings, it is important to distinguish between longitudinal and transverse SS breaking, depending on whether open strings can propagate or not along the SS direction .

The paper is organized as follows. In section 2 the general form of $`\rho (R)`$ for both closed and open strings is shown, for both longitudinal and transverse SS breaking. The 9D IIB orientifold is discussed in section 3 and in section 4 the class of orbifold/orientifold models on twisted ALE spaces are considered. Section 5 is devoted to the 4D chiral orientifold. Finally, some conclusions are given in section 6. In an appendix, we collect some details on the technique of unfolding the torus fundamental domain to provide an analytic integration of the torus one-loop partition function.

## 2 General form of the vacuum energy

The vacuum energy at one-loop level in any quantum field theory depends simply on the mass spectrum of the theory and its degeneracy. In $`D`$ space-time dimensions, in a Schwinger proper time parametrization, it can be written as

$$E_D=V_{D1}_0^{\mathrm{}}\frac{dt}{2t}\underset{i}{}d_i()^{2J_i}\frac{d^Dk}{(2\pi )^D}e^{t(k^2+m_i^2)},$$
(1)
with $`J_i`$ and $`m_i`$ the spin and mass of the state $`i`$, $`d_i`$ its degeneracy, $`k`$ its euclidean momentum and $`V_{D1}`$ the volume of the $`D1`$ spatial dimensions. More succinctly,

$$E_D=_0^{\mathrm{}}\frac{dt}{2t}\mathrm{Tr}_{}e^{tH},$$
(2)
where $``$ is the full Hilbert space of our system (including $`V_{D1}`$) and $`H=k^2+m^2`$. We are interested in the explicit form of (1) for orbifold and orientifold-derived theories, with the usual tower of massive string states, where supersymmetry is broken by a SS mechanism. We focus our attention to the case in which the twisted periodicity conditions defining the SS breaking are taken only along one direction, a circle of radius $`R`$, henceforth denoted SS direction. As we will see shortly, similar considerations apply also when the SS direction is an $`S^1/𝐙_2`$ orbifold, such as in the chiral 4D model considered in section 5. More general configurations with two (or more) SS directions are possible and have been explicitly constructed , but will not be considered here.

It is convenient to treat separately the contribution to the vacuum energy of states propagating or not along the SS direction. We denote the symmetry breaking in the two cases respectively as longitudinal and transverse SS breaking.

### 2.1 Longitudinal SS breaking

This is the situation that will concern us mostly, since it applies to all closed strings and to open strings on $`D9`$-branes. In both cases, the vacuum energy can be written as in (1). For open strings this is clear, $`t`$ being the modulus of the annulus and $`d_i`$ the string degeneracy of the state. For closed strings the situation is more complicated because the $`PSL(2,)`$ modular invariance of the torus restricts the integration over the modulus $`\tau `$ to the fundamental domain. Nevertheless, generalizing standard techniques (see e.g.  and the appendix) one can unfold the fundamental domain to the strip and rewrite the whole closed string contribution, including the Klein bottle term, in the form (1), where $`d_i`$ is now the string degeneracy of the state, with the level-matching conditions imposed by means of the $`\tau _1`$ integration. Quite interestingly, in this way only untwisted closed string states will explicitly appear in the computation. Similarly, winding modes along the SS direction will not be present, so that as far as the $`R`$ dependence is concerned, the whole amplitude looks effectively like that of a purely quantum field theory with an infinite number of states.

It is useful to distinguish closed and open string contributions to the vacuum energy. Recall that in string theory $`H^{(closed)}=L_0+\overline{L}_0=\alpha ^{}p^2/2+(N+\overline{N})/\alpha ^{}`$ and $`H^{(open)}=L_0=\alpha ^{}p^2+N/\alpha ^{}`$. Therefore, rescaling $`t\pi \alpha ^{}t`$ and $`t2\pi \alpha ^{}t`$ respectively in the two cases, one gets

$`\rho _D^{(C)}`$ $`=`$ $`{\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{dt}{2t^{\frac{D+2}{2}}}}{\displaystyle \underset{i}{}}d_i^{(C)}()^{2J_i}e^{\pi t\alpha ^{}m_i^2},`$
$`\rho _D^{(O)}`$ $`=`$ $`2^{D/2}{\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{dt}{2t^{\frac{D+2}{2}}}}{\displaystyle \underset{i}{}}d_i^{(O)}()^{2J_i}e^{2\pi t\alpha ^{}m_i^2},`$ (3)
where $`(C)`$ and $`(O)`$ stand respectively for closed and open, and we have defined the energy densities

$$\rho _D\frac{E_D(4\pi ^2\alpha ^{})^{\frac{D}{2}}}{V_{D1}},$$
(4)
with $`V_{D1}`$ the $`(D1)`$-dimensional spatial volume of the non-compact dimensions. The generic mass of a given closed and open string Kaluza-Klein (KK) state at level $`n`$ along the SS direction is respectively

$`m_{i,n}^2={\displaystyle \frac{2(N+\overline{N})}{\alpha ^{}}}+\mathrm{}+{\displaystyle \frac{(n+q)^2}{R^2}}m_i^2+{\displaystyle \frac{(n+q)^2}{R^2}},`$
$`m_{i,n}^2={\displaystyle \frac{N}{\alpha ^{}}}+\mathrm{}+{\displaystyle \frac{[n+q(G)]^2}{R^2}}m_i^2+{\displaystyle \frac{[n+q(G)]^2}{R^2}},`$ (5)
where $`q`$ is the twisted charge given by the SS breaking, $`N`$ and $`\overline{N}`$ are the string oscillator numbers and the dots stand for the KK and winding mode contributions along the other compact directions. The index $`i`$ in (3) includes thus a sum over $`N`$, $`\overline{N}`$, KK and winding modes over all the compact directions (but the SS direction) of states of given charge $`q`$ and then a sum over all possible twists $`q`$. Since the SS breaking can (and must) be implemented in the gauge sector as well, for open strings $`q`$ depends also on the gauge degrees of freedom, $`q=q(G)`$. The $`m_i^2`$ mass terms are typically functions of the geometric moduli of the compactification, except $`R`$, but for simplicity of notation this dependence will be left implicit in the following. Eq.(3) can be rewritten as

$`\rho _D^{(C)}`$ $`=`$ $`{\displaystyle \underset{n,i}{}}{\displaystyle \underset{F=0,1}{}}{\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{dt}{2t^{\frac{D+2}{2}}}}d_i(q,F)e^{i\pi F}e^{\pi t\alpha ^{}\left(m_i^2+\frac{(n+q)^2}{R^2}\right)},`$
$`\rho _D^{(O)}`$ $`=`$ $`2^{D/2}{\displaystyle \underset{n,i,G}{}}{\displaystyle \underset{F=0,1}{}}{\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{dt}{2t^{\frac{D+2}{2}}}}d_i(q,F)e^{i\pi F}e^{2\pi t\alpha ^{}\left(m_i^2+\frac{[n+q(G)]^2}{R^2}\right)},`$ (6)
where $`F`$ is the space-time fermion number operator, $`G`$ denotes a sum over the gauge indices and $`d_i(q,F)`$ are the string degeneracy factors, in general depending on $`q`$ and $`F`$, that include also the degeneracy arising from the expansion of the modular functions. By a Poisson resummation on the index $`n`$ and some algebra, it is not difficult to explicitly compute $`\rho _D^{(C)}`$ and $`\rho _D^{(O)}`$. It is convenient to separate the $`m_i^2=0`$ contribution, denoted by $`\rho _{D,0}^{(C,O)}`$, from the remaining ones. One gets for both closed and open strings:

$$\rho _{D,0}^{(C,O)}=\left(\frac{\sqrt{\alpha ^{}}}{R}\right)^D\frac{\mathrm{\Gamma }(\frac{D+1}{2})}{\pi ^{\frac{D+1}{2}}}\underset{q;F=0,1}{}e^{i\pi F}d_0(q,F)\frac{\mathrm{Li}_{D+1}(e^{2i\pi q})+\mathrm{Li}_{D+1}(e^{2i\pi q})}{2},$$
(7)
where

$$\mathrm{Li}_p(z)=\underset{w=1}{\overset{\mathrm{}}{}}\frac{z^w}{w^p}$$
(8)
are the polylogarithm functions, the sum over the gauge degrees of freedom for open strings is implicit, and we denoted by $`d_0`$ the degeneracy of states with $`m_i^2=0`$. As far as the $`m_i^20`$ contributions are concerned, we get

$$\rho _{D,i}^{(C,O)}=2\left(\frac{\sqrt{\alpha ^{}}}{R}\right)^{\frac{D1}{2}}(\alpha ^{}m_i^2)^{\frac{D+1}{4}}\underset{q;F=0,1}{}e^{i\pi F}d_i(q,F)\underset{w=1}{\overset{\mathrm{}}{}}\frac{\mathrm{cos}(2i\pi qw)}{w^{\frac{D+1}{2}}}K_{\frac{D+1}{2}}(2\pi Rwm_i),$$
(9)
where $`K_n`$ are the modified Bessel functions and again the sum over gauge indices has been omitted. The full vacuum energy $`\rho _D`$ is obtained by summing the total closed and open string energy density contributions, $`\rho _D=\rho _D^{(C)}+\rho _D^{(O)}`$, where

$$\rho _D^{(C,O)}=\rho _{D,0}^{(C,O)}+\underset{i0}{}\rho _{D,i}^{(C,O)}$$
(10)
and $`i0`$ indicates that states with $`m_i^2=0`$ should not be included in the sum. Both $`\rho _D^{(C)}`$ and $`\rho _D^{(O)}`$ are finite, as expected by the non-local nature of the SS breaking. Potential divergences should be local in space-time, but locally SUSY is preserved and hence no divergences at all are present. Indeed, in both eqs.(7) and (9) there would be a potential $`R`$-independent UV divergence arising from the $`w=0`$ term, where the index $`w`$, entering in (8), is obtained by a Poisson resummation on the index $`n`$ of eq.(6). This term vanishes because at each mass level $`i`$ the total number of bosons, summed over all possible twists $`q`$, equal the total number of fermions:

$$\underset{q}{}d_i(q,0)=\underset{q}{}d_i(q,1)i.$$
(11)
String models with SS SUSY breaking generically have winding modes that become tachyonic below certain values of $`R`$, where the vacuum energy diverges. Though the $`m_i^2`$ mass terms defined in (5) are always positive, this divergence appears in (9) from the sum over all massive string states. As is well known, the degeneracy of massive string states, for large masses, has a leading exponential behavior $`d_i\mathrm{exp}(2\pi c\sqrt{\alpha ^{}}m)`$, with $`c`$ a given constant. On the other hand, for large values of its argument, the modified Bessel function $`K_n(z)`$ admits an asymptotic expansion whose leading term is $`\mathrm{exp}(z)`$. Hence, we see that the infinite sum over $`i`$ in (9) converges only for $`R>R_T=\sqrt{\alpha ^{}}c`$. When the SS twist $`q=F`$, one can easily recognize eqs.(7) and (9) to be closely related to the free energy of string/field theory-derived models, with $`1/T=2\pi R`$ (see e.g. ) and $`T_H=1/(2\pi R_T)`$ being the Hagedorn temperature.

The general form of the vacuum energy for $`m_iR>>1`$ is easily extracted. We see from (7) that the $`m_i^2=0`$ contributions is power-like in $`R`$, whereas for large $`R`$ (9) is exponentially suppressed in $`R`$. More precisely we get

$$\rho _D\frac{1}{R^D}\underset{q}{}\left[d_0(q,1)d_0(q,0)\right]C_D(q)+O\left(\frac{e^{m_iR}}{R^{\frac{D}{2}}}\right),$$
(12)
where $`d_0(q,1)`$ and $`d_0(q,0)`$ are the total (closed + open) number of fermionic and bosonic massless states in $`D+1`$ dimensions (before the SS compactification) with charge $`q`$ and $`C_D(q)`$ are certain functions easily obtained from (7). For $`𝐙_2`$ or $`𝐙_3`$ twists, in which we get two independent twists $`q=0,1/2`$ or $`q=0,1/3`$, the constraint (11) implies that for large $`R`$ the vacuum energy is dominated by the difference between the total number of fermionic and bosonic $`D`$-dimensional, rather than $`D+1`$-dimensional, massless states $`d_0(0,1)`$ and $`d_0(0,0)`$ ;

$$\rho _D\frac{d_0(0,1)d_0(0,0)}{R^D}+O\left(\frac{e^{m_iR}}{R^{\frac{D}{2}}}\right).$$
(13)
In these cases, an exponentially small one-loop cosmological constant requires $`d_0(0,0)=d_0(0,1)`$ . However, for more general twists, we notice that the leading power-like behavior can be vanishing in a non-trivial way, thanks to a compensation between bosonic and fermionic contributions with different twists, even if $`d_0(0,0)d_0(0,1)`$. It would be quite interesting to fully exploit this observation and see if there exist string models with a spectrum satisfying this property.

All the above considerations are easily generalized to the case in which the SS direction is an $`S^1/𝐙_2`$ orbifold. Bulk states propagating along the orbifold are now classified according to their $`𝐙_2`$ parities. The massive spectra of $`𝐙_2`$-even and $`𝐙_2`$-odd states differ only by the presence or not of a zero mode along the orbifold. Both contributions can be summed together in the form (6), where the KK level runs over all integers. Possible left-over $`n=0`$ terms in the process of recombining $`𝐙_2`$-even and $`𝐙_2`$-odd contributions must vanish, since SUSY is broken by the compactification and they do not depend on $`R`$. Equation (6) and all the analysis that follows is then still valid for the $`S^1/𝐙_2`$ orbifold. The remaining compact directions can be instead arbitrary, as far as SUSY is broken only by the twist on $`R`$. Their structure will affect the explicit form of $`m_i^2`$ as well as the degeneracy factors $`d_i(q,F)`$.

### 2.2 Transverse SS breaking

States that do not propagate along the SS direction do not have a KK decomposition along that direction and hence their contribution to the vacuum energy requires a separate analysis. In string theory, states of this kind can arise either as twisted closed strings located at fixed points orthogonal to the SS direction or as open strings on $`D`$-branes transverse to the SS direction. In our class of models, the first kind of states appears always with unbroken tree-level SUSY and hence will never contribute to the one-loop vacuum energy. The same applies also to open strings on $`D`$-branes transverse to the SS direction, that present unbroken SUSY at the classical level. The only exception arises for open strings stretched between $`D5`$-branes/$`O5`$-planes and $`\overline{D}5`$-branes/$`\overline{O}5`$-planes, where SUSY is broken at tree level, for the 4D model discussed in section 5. In this case, the one-loop open string amplitude is more conveniently expressed as a tree-level exchange of closed string states, propagating from one object to the other. This contribution can be thus summarized as follows:

$$\rho _DV_n\underset{i}{}\underset{\widehat{w}}{}Q^{(i)}\overline{Q}^{(i)}𝒢_d^{(i)}[\mathrm{\Delta }x(\widehat{w})],$$
(14)
where $`𝒢_d^{(i)}`$ is the $`d`$-dimensional propagator of a particle of mass $`m_i`$ and spin $`J_i`$ and $`V_n`$ is the volume of the compact longitudinal directions along which the states can propagate. $`\mathrm{\Delta }x(\widehat{w})`$ is the transverse $`D`$-brane/$`O`$-plane-$`\overline{D}`$-brane/$`\overline{O}`$-plane distance modulo windings, since a given closed string state can wind $`w`$ times along all transverse directions before ending on a $`D`$-brane/$`O`$-plane, and $`Q^{(i)}=Q_D^{(i)}/Q_O^{(i)}`$ is the $`D`$-brane/$`O`$-plane charge for each state $`i`$. The same applies to the anti-brane and anti-plane charges $`\overline{O}^{(i)}`$. As in last subsection, the sum over $`i`$ runs over the string and winding modes along the longitudinal directions, whereas the sum over $`\widehat{w}`$ runs over all the possible windings along the transverse $`d`$ dimensions. One should recall that the subscripts $`d`$ and $`n`$ in (14) represent the number of space-time directions in which closed strings propagate and thus it can be different for untwisted and various twisted closed string states. As in the last subsection, the massless $`i=0`$ contributions in (14) are power-like in $`R`$, whereas massive ones give an infinite sum over modified Bessel functions $`K_n`$ and thus are exponentially suppresed in $`R`$, for large $`R`$. The sign of the vacuum energy contribution is now given by the brane/plane charge and is always negative. This is intuitively clear, since objects with opposite charges feel an attractive potential between each other. The divergence due to the open string tachyon arising below a given radius is determined by looking at the asymptotic form of the modified Bessel functions and at their degeneracies for large masses.

## 3 A nine-dimensional model

We compute in this section the energy density $`\rho _9`$ of a simple 9D Type IIB orientifold compactified on $`S^1/(\mathrm{\Omega }\times 𝐙_2)`$, where $`\mathrm{\Omega }`$ is the world-sheet parity operator and $`𝐙_2`$ is generated by $`g`$, the product of a $`\pi R`$ translation $`\sigma `$ along the circle and $`()^F`$, with $`F`$ the space-time fermion number operator; $`g=\sigma ()^F`$. Such a computation has already been done in ; we review it in the following as a simple example of the kind of computation we are going to perform in the next sections. We are interested in the radius dependence of $`\rho _9=\rho _9(R)`$; for simplicity we do not include continuous Wilson lines but discuss the dependence of $`\rho _9(R)`$ on the twist matrices $`\gamma _g`$ that embed the $`𝐙_2`$ group in the Chan-Paton degrees of freedom. They can alternatively be considered as $`𝐙_2`$ discrete Wilson lines.

The construction of the model is straightforward and will not be presented in detail (see e.g. ). The only massless tadpoles are those for the dilaton, graviton (NSNS) and for the untwisted 10-form (RR), whose cancellation requires the presence of 32 $`D9`$-branes. All neutral fermions are anti-periodic along the circle and thus massive, with a mass $`1/R`$. In the twisted sector we get a tower of real would-be tachyons starting from $`\alpha ^{}m^2=R^2/(4\alpha ^{})2`$. The twist matrix $`\gamma _g`$ is arbitrary and must only satisfy the group algebra $`\gamma _g^2=\pm I`$. The gauge group is $`SO(n)\times SO(32n)`$, with massless fermions in the bifundamentals $`(𝐧,\mathrm{𝟑𝟐}𝐧)`$, for $`\gamma _g=\mathrm{diag}(I_n,I_{32n})`$. If $`\gamma _g`$ is taken traceless, the twisted tadpoles associated to the above would-be tachyons cancel as well. From this perspective, such a choice is different from the others. If $`\gamma _g`$ is chosen antisymmetric, one gets a $`U(16)`$ group, or subgroups thereof, with massless fermions in antisymmetric representations333In this model, local tadpole cancellation for both untwisted NSNS and RR tadpoles cannot be achieved. By choosing respectively $`\gamma _g`$ symmetric or antisymmetric one can cancel respectively the massive NSNS or RR tadpoles.. We focus in the following on symmetric twist matrices for definiteness, but since the energy density $`\rho _9(R)`$ depends on $`\gamma _g`$ only through its trace, the analysis applies equally well for other more general choices.

The full energy density of the model is obtained by summing the closed string contribution $`\rho _9^{(C)}`$ (torus+Klein bottle) and the open one, $`\rho _9^{(O)}`$ (annulus+Möbius strip).

### 3.1 Closed string contribution

Since the Klein bottle amplitude vanishes identically (the $`\mathrm{\Omega }`$-projection acts in a supersymmetric manner on the closed string spectrum), the whole contribution is given by the torus amplitude:

$$\rho _9^{(C)}=_{}\frac{d^2\tau }{2\tau _2^{11/2}}\frac{1}{2|\eta |^{24}}\underset{n,w}{}\left\{|\theta _2^4|^2\mathrm{\Lambda }_{n,w}()^n+|\theta _4^4|^2\mathrm{\Lambda }_{n,w+\frac{1}{2}}+|\theta _3^4|^2\mathrm{\Lambda }_{n,w+\frac{1}{2}}()^n\right\},$$
(15)
where

$$\mathrm{\Lambda }_{n,w}=q^{\frac{\alpha ^{}}{4}\left(\frac{n}{R}+\frac{wR}{\alpha ^{}}\right)^2}\overline{q}^{\frac{\alpha ^{}}{4}\left(\frac{n}{R}\frac{wR}{\alpha ^{}}\right)^2},$$
(16)
and $`\theta _i=\theta _i(0|\tau )`$ are the usual theta functions, defined as in . In eq.(15) and throughout all the paper, we often omit the modular dependence of $`\theta _i`$ on $`\tau `$ and always leave implicit its vanishing argument $`z=0`$. Using the unfolding technique (see the appendix), (15) can be written as

$$\rho _9^{(C)}=\frac{R}{\sqrt{\alpha ^{}}}_0^{\mathrm{}}\frac{d\tau _2}{2\tau _2^6}\underset{N}{}d_N^2e^{4\pi \tau _2N}\underset{p=}{}\left[\frac{1()^p}{2}\right]e^{p^2{\scriptscriptstyle \frac{\pi R^2}{4\tau _2\alpha ^{}}}},$$
(17)
where

$$\frac{1}{2}_{1/2}^{1/2}𝑑\tau _1\frac{|\theta _2^4|^2}{|\eta |^{24}}\underset{N}{}d_N^2e^{4\pi \tau _2N}.$$
(18)
Eq.(17) can also be written by using a Poisson resummation as

$$\rho _9^{(C)}=_0^{\mathrm{}}\frac{d\tau _2}{2\tau _2^{11/2}}\underset{N}{}d_N^2e^{4\pi \tau _2N}\underset{n}{}\left[e^{\pi \tau _2\frac{\alpha ^{}n^2}{(R/2)^2}}e^{\pi \tau _2\frac{\alpha ^{}(n+1/2)^2}{(R/2)^2}}\right].$$
(19)
Eq.(19) is precisely of the general form (6) with $`D=9`$, $`q=F/2`$, $`i=N`$, $`\alpha ^{}m_N^2=4N`$, $`d_N(0,0)=d_N(1/2,1)=d_N^2`$, $`d_N(0,1)=d_N(1/2,0)=0`$, and $`R/2R`$. The rescaling of the radius is standard and is due to the identification $`xx+\pi R`$ induced by the freely acting action. We thus notice, as mentioned in the last section, that the full string theory contribution to the cosmological constant (all KK, winding, untwisted and twisted string states) is automatically encoded in a field theory contribution with only untwisted states, whose KK modes are shifted by the SS mechanism (bosons/fermions periodic/anti-periodic along the SS circle R) and a reduced massive string spectrum, where only the diagonal $`N=\overline{N}`$ states contribute . The integration on $`\tau _2`$ in (19) can thus be read off directly from (7) for $`N=0`$ and (9) for $`N0`$.

### 3.2 Open string contribution

The annulus and Möbius strip contributions, in the absence of Wilson lines, are

$`\rho _9^A=`$ $`\left(\mathrm{Tr}\gamma _g\right)^2{\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{dt}{4(2t)^{11/2}}}{\displaystyle \frac{\theta _2^4}{\eta ^{12}}}(it){\displaystyle \underset{m}{}}()^me^{2\pi tm^2\frac{\alpha ^{}}{R^2}},`$
$`\rho _9^M=`$ $`+32{\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{dt}{4(2t)^{11/2}}}{\displaystyle \frac{\theta _2^4}{\eta ^{12}}}(it{\displaystyle \frac{1}{2}}){\displaystyle \underset{m}{}}()^me^{2\pi tm^2\frac{\alpha ^{}}{R^2}}.`$ (20)
The integrations of (20) are easily performed and one gets

$`\rho _{9,0}^A+\rho _{9,0}^M=`$ $`\left[32(\mathrm{Tr}\gamma _g)^2\right]\left({\displaystyle \frac{\sqrt{\alpha ^{}}}{R}}\right)^92^{\frac{15}{2}}d_0{\displaystyle \frac{\mathrm{\Gamma }(5)}{\pi ^5}}\xi (10)(22^9),`$ (21)
$`\rho _{9,N}^A+\rho _{9,N}^M=`$ $`\left[32()^N(\mathrm{Tr}\gamma _g)^2\right]\left({\displaystyle \frac{\sqrt{\alpha ^{}}}{R}}\right)^42^{\frac{7}{2}}d_NN^{\frac{5}{2}}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1()^n}{n^5}}K_5\left[\pi Rn\sqrt{{\displaystyle \frac{N}{\alpha ^{}}}}\right],`$
with $`d_N`$ as in (18). It is not difficult to realize that these expressions are of the general form (6) with $`D=9`$, $`i=N`$, $`\alpha ^{}m_N^2=N`$ and $`q_{ij}(G)=F/2+\widehat{q}_i+\widehat{q}_j`$, where $`i,j`$ run over the fundamental representation of $`G=SO(32)`$. The parameters $`\widehat{q}_i=0,1/2`$ are related to the eigenvalues of $`\gamma _g`$ by

$$\gamma _g=\mathrm{diag}(e^{2i\pi \widehat{q}_1},\mathrm{},e^{2i\pi \widehat{q}_{32}}).$$
(22)
As mentioned, the resulting 9D gauge group is $`SO(n)\times SO(32n)`$ if $`\gamma _g=\mathrm{diag}(I_n,I_{32n})`$, with bosons and fermions in antisymmetric, symmetric and bifundamental representations, depending on the KK and string mass level $`N`$. Notice that the leading open string contribution to the energy density is positive whenever $`14n18`$, with a maximum for $`n=16`$.

The full energy density $`\rho _9=\rho _9^{(C)}+\rho _9^{(O)}=\rho _9^{(C)}+\rho _9^A+\rho _9^M`$ can then be numerically evaluated as a function of $`R`$ by truncating the infinite sums of modified Bessel functions.

The $`R`$-dependence of $`\rho _9`$ crucially depends on $`\gamma _g`$. For $`n14`$ or $`n18`$, $`\rho _9<0`$ monotonically and leads the system to a tachyonic instability, like in the early work of . For $`15n17`$, $`\rho _9`$ can be positive and a maximum close to $`R_T`$ is obtained for $`n=15`$ or $`n=17`$.

In figure (1) we show these different behaviors plotting $`\rho _9(R)`$ for $`n=16`$ (blue/solid line), $`n=15n=17`$ (green/dotted-dashed line) and $`n=14n=18`$ (red/dashed line).

## 4 Strings on twisted ALE spaces

An interesting class of models whose energy density can be studied are orbifold or orientifold models on Asymptotically Locally Euclidean (ALE) spaces, non-trivially fibered along an $`S^1`$, of the form $`(\times S^1)/𝐙_N`$, or their compact versions $`(T^2\times S^1)/𝐙_N`$. The $`𝐙_N`$ generator is a product of a $`4\pi /N`$ rotation along $``$ and of a translation of $`2\pi R/N`$ along the circle. The rotation is taken to be of angle $`4\pi /N`$ so that $`g^N=1`$ on spinors. A non-trivial fibration is necessary to implement the Scherk-Schwarz supersymmetry breaking and to lift the mass of the would-be tachyons. Upon compactification on $`S^1`$, such spaces give rise to a Melvin background  (see  and  for strings and $`D`$-branes on Melvin backgrounds).

Consider then Type IIB string theory on $`(\times S^1)/𝐙_N`$ or $`(T^2\times S^1)/𝐙_N`$ (as we will see, our results are the same for the non-compact or compact version). In order to keep our analysis as simple as possible, we focus on $`N`$ odd. Other values of $`N`$ would require the introduction of $`O7`$-planes and $`D7`$-branes, in addition to $`D9`$-branes and the $`O9`$-plane, when considering orientifolds. Moreover, for $`N`$ odd, would-be tachyons appear localized in space-time and an understanding of the possible tachyonic instabilities can be obtained .

Uncharged fermions in all these models are massive for $`N3`$, with a mass $`1/R`$. In each $`k`$-twisted sector we get a tower of complex would-be tachyons starting from $`\alpha ^{}m^2=4\frac{k}{N}+\left(\frac{k}{N}\right)^2\frac{R^2}{\alpha ^{}}`$.

### 4.1 Orbifold models

The only relevant one-loop world-sheet surface in this case is the torus. Its contribution to the vacuum energy density $`\rho _7^{N,T}(R)`$ can be written, using the unfolding technique discussed in the appendix (see  for a more detailed discussion), as an integral over the strip of the untwisted sector only

$$\rho _7^{N,T}=\frac{1}{N}\frac{R}{\sqrt{\alpha ^{}}}_0^{\mathrm{}}\frac{d\tau _2}{2\tau _2^5}\underset{k=1}{\overset{N1}{}}\underset{M}{}\left[d_M(B)^kd_M(F)^k\right]e^{4\pi \tau _2M}\underset{n}{}e^{\frac{\pi R^2}{\alpha ^{}\tau _2}\left(n+\frac{k}{N}\right)^2},$$
(23)
where we have defined for later convenience the coefficients $`d_M(B,F)^k`$ as

$$\underset{M}{}\left[d_M(B)^kd_M(F)^k\right]e^{4\pi \tau _2M}=\frac{d\tau _1}{4}\left|2\mathrm{sin}(\frac{2\pi k}{N})\frac{_{\alpha ,\beta }\eta _{\alpha \beta }\theta \left[\genfrac{}{}{0pt}{}{\alpha }{\beta }\right]^3\theta \left[\genfrac{}{}{0pt}{}{\alpha }{\beta +\frac{2k}{N}}\right]}{\eta ^9\theta \left[\genfrac{}{}{0pt}{}{\frac{1}{2}}{\frac{1}{2}+\frac{2k}{N}}\right]}\right|^2.$$
(24)
If, similarly to the 9D case of last section, the torus contribution to the vacuum energy is exactly encoded in a field theory-like contribution with only KK and untwisted states, (23) should admit a rewriting such as (6), where the twist $`q`$ will be given by the Lorentz $`SO(2)`$ charge of the states along the twisted directions, with the correct degeneracies. It is useful to compute the latter in some detail for the massless states of IIB string theory.

The above orbifold breaks the Lorentz group $`SO(10)SO(8)\times SO(2)`$. A generic field $`\mathrm{\Phi }`$ will have the following periodicity conditions along $`S^1`$;444Notice that due to the freely acting action, the radius $`R`$ entering in (25) is $`N`$ times smaller than the $`R`$ appearing in eq.(23).

$$\mathrm{\Phi }(y+2\pi R)=e^{2i\pi \widehat{q}\alpha }\mathrm{\Phi }(y),$$
(25)
where $`\widehat{q}`$ is its charge under the $`SO(2)U(1)`$ internal Lorentz group and $`\alpha =2/N`$ is the twist induced by the $`𝐙_N`$ action. Upon reduction on $`S^1`$, a field with charge $`\widehat{q}`$ will have a tower of KK states with masses $`M_n^{(q)}=(n+q)/R`$, where $`q=\widehat{q}\alpha `$. Massless states are not present whenever $`q`$.

The $`U(1)`$ charges of the massless states of IIB string theory are easily obtained. The Ramond-Ramond (RR) four, two and zero-forms have the following $`SO(8)\times U(1)`$ decomposition:

$`\mathrm{𝟑𝟓}`$ $`=`$ $`\mathrm{𝟏𝟓}_0\mathrm{𝟏𝟎}_1\mathrm{𝟏𝟎}_1,`$
$`\mathrm{𝟐𝟖}`$ $`=`$ $`\mathrm{𝟏𝟓}_0\mathrm{𝟔}_1\mathrm{𝟔}_1\mathrm{𝟏}_0,`$
$`\mathrm{𝟏}`$ $`=`$ $`\mathrm{𝟏}_0,`$ (26)
where the subscript denotes the $`U(1)`$ Lorentz charge $`\widehat{q}`$. One gets for the Neveu-Schwarz/Neveu-Schwarz (NSNS) graviton, $`B`$-field and dilaton:

$`\mathrm{𝟑𝟓}`$ $`=`$ $`\mathrm{𝟐𝟏}_0\mathrm{𝟔}_1\mathrm{𝟔}_1\mathrm{𝟏}_2\mathrm{𝟏}_2,`$
$`\mathrm{𝟐𝟖}`$ $`=`$ $`\mathrm{𝟏𝟓}_0\mathrm{𝟔}_1\mathrm{𝟔}_1\mathrm{𝟏}_0,`$
$`\mathrm{𝟏}`$ $`=`$ $`\mathrm{𝟏}_0.`$ (27)
For generic $`N2,4`$ we thus get 70 bosonic massless states. For fermions one has two copies of:

$`\mathrm{𝟓𝟔}`$ $`=`$ $`\mathrm{𝟐𝟒}_{1/2}\mathrm{𝟐𝟒}_{1/2}\mathrm{𝟒}_{3/2}\mathrm{𝟒}_{3/2},`$
$`\mathrm{𝟖}`$ $`=`$ $`\mathrm{𝟒}_{1/2}\mathrm{𝟒}_{1/2}.`$ (28)
Notice that all fermions are always massive except for the case $`N=3`$, where we get 16 massless states from the decomposition of the two gravitinos.

We are now ready to explicitly show how these massless states and the corresponding twists arise from (23). We focus on the $`N=3`$ case, but the analysis can be generalized to the other cases. Reintroducing also the vanishing $`k=0`$ term in (23) we have ($`\tau _2t`$ and $`RNR`$):

$`\rho _7^T={\displaystyle \frac{R}{\sqrt{\alpha ^{}}}}{\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{dt}{2t^5}}{\displaystyle \underset{M}{}}{\displaystyle \underset{n}{}}e^{4\pi tM}`$ $`\{[d_M^0(B)d_M^0(F)]e^{(3n)^2\frac{\pi R^2}{\alpha ^{}t}}+`$ (29)
$`\left[d_M^1(B)d_M^1(F)\right]e^{(3n+1)^2\frac{\pi R^2}{\alpha ^{}t}}+`$
$`[d_M^2(B)d_M^2(F)]e^{(3n+2)^2\frac{\pi R^2}{\alpha ^{}t}}\}.`$
Notice that the above sums can be rewritten as follows:

$`{\displaystyle \underset{n}{}}e^{(3n)^2\frac{\pi R^2}{\alpha ^{}t}}`$ $`=`$ $`{\displaystyle \underset{n}{}}{\displaystyle \frac{1+e^{\frac{2i\pi n}{3}}+e^{\frac{2i\pi n}{3}}}{3}}e^{n^2\frac{\pi R^2}{\alpha ^{}t}},`$ (30)
$`{\displaystyle \underset{n}{}}\left[e^{(3n+1)^2\frac{\pi R^2}{\alpha ^{}t}}+e^{(3n+2)^2\frac{\pi R^2}{\alpha ^{}t}}\right]`$ $`=`$ $`{\displaystyle \underset{n}{}}{\displaystyle \frac{2e^{\frac{2i\pi n}{3}}e^{\frac{2i\pi n}{3}}}{3}}e^{n^2\frac{\pi R^2}{\alpha ^{}t}}.`$
Since $`d_M(B,F)^1=d_M(B,F)^2`$ $`M`$, eq.(23) can be rewritten, by performing a Poisson resummation in $`n`$, as

$`\rho _7^T={\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{dt}{2t^{\frac{9}{2}}}}{\displaystyle \underset{n}{}}{\displaystyle \underset{M}{}}e^{4\pi tM}\{`$ $`\left[d_M(0,0)d_M(0,1)\right]e^{\frac{\pi t\alpha ^{}n^2}{R^2}}`$
$`+`$ $`[d_M({\displaystyle \frac{1}{3}},0)d_M({\displaystyle \frac{1}{3}},1)]e^{\frac{\pi t\alpha ^{}(n+1/3)^2}{R^2}}\},`$ (31)
where

$`d_M(0,0/1)`$ $`=`$ $`{\displaystyle \frac{1}{3}}\left[d_M(B/F)^0+2d_M(B/F)^1\right],`$
$`d_M({\displaystyle \frac{1}{3}},0/1)`$ $`=`$ $`{\displaystyle \frac{2}{3}}\left[d_M(B/F)^0d_M(B/F)^1\right]`$ (32)
are the total number of bosonic and fermionic states at level $`M`$ with twist 0 and 1/3 respectively. One can easily check that for $`M=0`$ the above coefficients precisely coincide with the field theory results above. Indeed, from (24) one finds $`d_0(B)^0=d_0(F)^0=128`$, $`d_0(B)^1=d_0(B)^2=41`$, $`d_0(F)^1=d_0(F)^2=40`$ and hence $`d_0(0,0)=70`$, $`d_0(0,1)=16`$, $`d_0(1/3,0)=58`$, $`d_0(1/3,1)=112`$, as expected. Eq.(4.1) is hence precisely of the form (6) with $`D=7`$, $`i=M`$, $`\alpha ^{}m_M^2=4M`$ and $`q=0,1/3`$.

Interestingly, even for $`𝐙_3`$ (and most likely for all $`𝐙_N`$, at least with $`N`$ odd) the full string theory computation is encoded in a field theory computation where only untwisted states with the uncompactified level matching conditions $`M=\overline{M}`$ enter.

The form of $`\rho _7^N(R)`$ for various $`𝐙_N`$ orbifold models is reported in figure 2. As anticipated in the introduction, we find that $`\rho _7^N(R)>\rho _7^M(R)`$ for $`N<M`$ and for all the values of $`R`$ for which both energy densities are well defined.

### 4.2 Orientifold models

Orientifold models are obtained as usual by modding out the above orbifolds by the world-sheet parity operator $`\mathrm{\Omega }`$. As in the 9D model discussed in last section, massless tadpoles are cancelled by introducing 32 $`D9`$-branes. We do not include continuous Wilson lines but consider the dependence of $`\rho _7^N(R)`$ on the twist matrices $`\gamma _g`$, whose only constraint comes from the group algebra: $`\gamma _g^N=\pm I`$. The resulting gauge group depends on the precise form of $`\gamma _g`$, although the energy density is sensitive only on its trace. In addition to the torus amplitude (23), we have now to consider also the Klein bottle, annulus and Möbius strip surfaces.

The Klein bottle contribution, $`\rho _7^{N,K}`$, is easily obtained and has a form similar to (23). Defining

$`{\displaystyle \underset{M}{}}[D_M(B)^{2k}D_M(F)^{2k}]e^{4\pi Mt}=2\mathrm{sin}({\displaystyle \frac{4\pi k}{N}}){\displaystyle \frac{_{\alpha ,\beta }\eta _{\alpha \beta }\theta \left[\genfrac{}{}{0pt}{}{\alpha }{\beta }\right]^3\theta \left[\genfrac{}{}{0pt}{}{\alpha }{\beta +\frac{4k}{N}}\right]}{2\eta ^9\theta \left[\genfrac{}{}{0pt}{}{\frac{1}{2}}{\frac{1}{2}+\frac{4k}{N}}\right]}}(2it),`$ (33)
$`\rho _7^{N,K}`$ can be in fact nicely combined with $`\rho _7^{N,T}`$ above, so that the full closed string contribution reads exactly as (23) with $`d_M(B,F)^k[d_M(B,F)^k+D_M(B,F)^k]/2`$. As expected, the degeneracies of the massless states again agree with the field theory expectations. This can be easily checked by noting that from (33) one gets $`D_0(B)^2=D_0(B)^4=9`$, $`D_0(F)^2=D_0(F)^4=0`$ and that the RR four and zero-form, the NSNS $`B`$ field and half of the fermions are now projected out, leaving the decomposition of the remaining states as above. The open string contribution $`\rho _7^{N,(O)}`$ is easily computed:

$`\rho _7^{N,(O)}=2^{7/2}{\displaystyle _0^{\mathrm{}}}{\displaystyle \frac{dt}{2t^{9/2}}}{\displaystyle \underset{M,n,k=1}{\overset{N1}{}}}\left[\left(\mathrm{Tr}\gamma _\mathrm{k}\right)^2(1)^M\mathrm{Tr}\gamma _{2\mathrm{k}}\right]e^{2\pi in\frac{k}{N}}e^{2\pi t\alpha ^{}\frac{n^2}{R^2}}D_M^ke^{2\pi tM},`$ (34)
where the coefficients $`D_N^k`$ are the same as those defined in (33) and $`\gamma _k=(\gamma _g)^k`$. Eq.(34) is not manifestly in the form (6) because the embedding of the SS breaking in the gauge sector through the twist matrices $`\gamma _g`$ requires a little bit of algebra. However, we do not need to work out these details, since the vacuum energy depends only on the trace of $`\gamma _g`$. The integration in $`t`$ is by now standard and leads to a power-like behavior in $`R`$ for the $`M=0`$ terms and a sum over modified Bessel functions for $`M0`$.

The form of $`\rho _7^N(R)`$ for various $`𝐙_N`$ orientifold models is reported in figure 3 for a proper choice of $`\gamma _g`$ that maximizes $`\rho _7^N(R)`$. As can be seen from the figure, the qualitative structure of $`\rho _7^N`$ is not modified by the orientifold projection. In particular, we still find that $`\rho _7^N(R)>\rho _7^M(R)`$ for $`N<M`$. For different choices of $`\gamma _g`$, $`\rho _7^N`$ has always the same form as in figure 3, but the above ordering of $`\rho _7^N`$ can be lost.

The vacuum energy density $`\rho _7^N(R)`$ receives a non-vanishing contribution only for states that do not propagate along the two twisted directions for both orbifold and orientifold models. This implies that the above computation apply equally well for non-compact $`(\times S^1)/𝐙_N`$ or compact $`(T^2\times S^1)/𝐙_N`$ models. The energy density is located at the tip of the cone or at the fixed points of $`T^2`$, in the two cases. This is expected since these are precisely the loci where would-be tachyons are localized.

It is interesting to notice that if we plot $`\rho _7^N`$ as a function of $`R/N`$, the effective radius of the SS direction, we get again the behavior as in figures 2 and 3, but now $`\rho _7^N(R/N)<\rho _7^M(R/M)`$ for $`N<M`$, $`N,M`$ and $`R`$. It would be interesting to have some dynamical understanding of this ordering of the energy densities for this class of models.

The $`𝐙_2`$ orbifold and orientifold models can also be considered. In this case space-time is flat with fermions antiperiodic along $`R`$ and thus we recover exactly the previous 9D model of section 2 or its compactified version on a $`T^2`$ torus.

## 5 A four dimensional model

In this section we study the vacuum energy density of a 4D IIB orientifold model compactified on $`T^6/(𝐙_6^{}\times 𝐙_2^{})`$. The $`𝐙_6^{}`$ group is generated by the element $`\theta `$, acting as SUSY rotations of angles $`2\pi v_i^\theta `$ in the three tori $`T_i^2(i=1,2,3)`$, with $`v_i^\theta =1/6(1,2,3)`$ . The complex structure of the first and second torus is fixed by the $`𝐙_6^{}`$ action, whereas the third torus can be taken rectangular, since the action is a $`𝐙_2`$ reflection. We define $`P_{1,2}`$ such that the volume of the $`i`$-th torus equals to $`4\pi ^2P_i^2/[3(1+\delta _{i,1})]`$, and $`R`$ and $`S`$ to be the radii of the two circles of $`T_3^2`$. The $`𝐙_2^{}`$ group is generated by $`\beta =\sigma ()^F`$, the product of a translation $`\sigma `$ of length $`\pi R`$ along the Scherk-Schwarz direction, the circle with radius $`R`$, and $`()^F`$, where $`F`$ is the 4D space-time fermion number. The model has $`D9`$, $`D5`$ and $`\overline{D}5`$-branes. All $`D5`$-branes are located at the same point in space-time and at $`y=0`$ and $`y=\pi R`$ along the SS direction. Similarly, all $`\overline{D}5`$-branes are at the same point in space-time and at $`y=\pi R/2`$ and $`y=3\pi R/2`$ along the SS direction.

We compute the cosmological constant $`\rho _4`$ as a function of the SS radius and the other moduli. Contrary to the previous models, $`\rho _4`$ gets a non-vanishing contribution both from states with longitudinal and transverse SS SUSY breaking. For simplicity we do not include continuous Wilson lines. A detailed construction of the model is presented in .

### 5.1 Longitudinal SS breaking contribution

There are two main longitudinal contributions coming from the torus amplitude and the open sector of $`D9`$-branes.

#### Closed string contribution.

Due to the unfolding technique, the longitudinal closed string contribution to the energy density, denoted by $`\rho _4^{(C,l)}`$, can be written as

$$\rho _4^{(C,l)}=_0^{\mathrm{}}\frac{dt}{2t^3}\underset{i}{}\underset{j=U,T}{}d_i^{(C,j)}e^{\pi t\alpha ^{}m_{i}^{2}{}_{}{}^{(C,j)}}\underset{n}{}[e^{\pi t\alpha ^{}\frac{(2n)^2}{R^2}}e^{\pi t\alpha ^{}\frac{(2n+1)^2}{R^2}}],$$
(35)
where $`U`$ and $`T`$ stand for the untwisted and $`\theta ^2/\theta ^4`$-twisted string sectors—the only sectors with a non-vanishing contribution—and the index $`i`$ includes the sum over the string, KK and winding modes over the non-SS directions with the level matching conditions imposed by means of the $`\tau _1`$ integration555This equals $`\overline{N}=N+_{a=1}^5m_an_a`$ in the untwisted sector and $`N=\overline{N}+w_1n_1`$ in the $`\theta ^2/\theta ^4`$ twisted ones..

The masses of the $`i`$-th level for the untwisted ($`j=U`$) and twisted ($`j=T`$) sectors are given by the following expression:

$`m_{i}^{2}{}_{}{}^{(C,j)}`$ $`=`$ $`{\displaystyle \frac{2(N+\overline{N})}{\alpha ^{}}}+{\displaystyle \frac{\alpha ^{}n_1^2}{S^2}}+{\displaystyle \frac{w_1^2S^2}{\alpha ^{}}}+\delta _{j,U}{\displaystyle \underset{a=1}{\overset{2}{}}}\left[{\displaystyle \frac{2\alpha ^{}N_a^2}{\sqrt{3}P_a^2}}+{\displaystyle \frac{\sqrt{3}W_a^2P_a^2}{2\alpha ^{}}}\right],`$ (36)
where we have defined the combinations of KK and winding modes $`N_a^2`$ and $`W_a^2`$ as

$`N_a^2`$ $`=`$ $`n_{2a}^2+n_{2a+1}^2+n_{2a}n_{2a+1},`$
$`W_a^2`$ $`=`$ $`w_{2a}^2+w_{2a+1}^2w_{2a}w_{2a+1}.`$ (37)
The degeneracy of the $`i`$-th state is given by $`d_i^{(C,j)}=d_N^jd_{\overline{N}}^j/24`$, with

$`{\displaystyle \underset{N,\overline{N}}{}}d_N^jd_{\overline{N}}^jq^N\overline{q}^{\overline{N}}=\{\begin{array}{cc}|\theta _2^4|^2|\eta |^{24},\hfill & j=U,\\ 9|\theta _2^2|^2|\eta |^{12}\left|\theta \left[\genfrac{}{}{0pt}{}{1/2}{1/3}\right]\theta \left[\genfrac{}{}{0pt}{}{1/2}{5/6}\right]^1\right|^4,\hfill & j=T.\end{array}`$ (40)
Once again, eq.(35) is of the same form as (3).

#### Open string contribution.

The longitudinal open string contribution to the energy density, denoted by $`\rho _4^{(O,l)}`$, is

$$\rho _4^{(O,l)}=_0^{\mathrm{}}\frac{dt}{4(2t)^3}\underset{i}{}\underset{j=U,T}{}d_i^{(O,j)}e^{2\pi t\alpha ^{}m_{i}^{2}{}_{}{}^{(O,j)}}\underset{n}{}[e^{\pi t\alpha ^{}\frac{(2n)^2}{R^2}}e^{\pi t\alpha ^{}\frac{(2n+1)^2}{R^2}}],$$
(41)
where

$`m_{i}^{2}{}_{}{}^{(O,j)}`$ $`=`$ $`{\displaystyle \frac{N}{\alpha ^{}}}+{\displaystyle \frac{\alpha ^{}n_1^2}{S^2}}+\delta _{j,U}{\displaystyle \underset{a=1}{\overset{2}{}}}{\displaystyle \frac{2\alpha ^{}N_a^2}{\sqrt{3}P_a^2}},`$ (42)
with $`N_a^2`$ as in (37) and:

$`d_i^{(O,U)}`$ $`=`$ $`\left[\left(\mathrm{Tr}\gamma _\beta \right)^232(1)^N\right]{\displaystyle \frac{d_N^U}{6}},`$ (43)
$`d_i^{(O,T)}`$ $`=`$ $`\left[\left(\mathrm{Tr}\gamma _\theta ^2\gamma _\beta \right)^2\left(\mathrm{Tr}\gamma _\theta ^4\right)(1)^N\right]{\displaystyle \frac{d_N^T}{6}},`$ (44)
with $`\gamma _\theta `$ and $`\gamma _\beta `$ the twist matrices associated respectively to the SUSY and non-SUSY twists $`\theta `$ and $`\beta `$, and $`d_N^{U/T}`$ defined as in (40).

The untwisted contribution is essentially the same as the open contribution (20) of the 9D model, the only difference being an overall factor and the Kaluza-Klein lattices from the compactified directions. As in the 9D case, $`\gamma _\beta `$ is unconstrained (provided that $`\gamma _\beta ^2=\pm I`$). The value of $`\mathrm{Tr}\gamma _\theta ^4`$ is fixed to be $`8`$ by tadpole cancellation.

### 5.2 Transverse SS breaking contribution

This contribution arises only when considering $`D5`$, $`\overline{D}5`$-branes, as well as $`O5`$ and $`\overline{O}5`$-planes. As shown in section 2, these terms can be written as a sum of propagators of closed string states with mass $`m_i`$, propagating along the number of compact dimensions that are orthogonal to the brane.

Tadpole cancellations are automatically encoded when we add up the Klein bottle, annulus and Möbius strip contributions in the closed string channel. Notice that the cancellation is obtained for all the KK modes along the SS direction due to the choice of the brane positions, ensuring the local tadpole cancellation in that dimension.

In terms of closed string states, the full amplitude is a sum of two contributions, one coming from the untwisted and one from the $`\theta ^2/\theta ^4`$ twisted string sector.

#### Untwisted string contribution

This is given by the sum of the relevant Klein bottle, annulus and Möbius strip amplitudes. Denoting by $`\rho _4^{(U,t)}`$ the contribution to the vacuum energy of this sector, we get:

$`\rho _4^{(U,t)}={\displaystyle \frac{8P_2^2\pi ^2}{3\alpha ^{}}}{\displaystyle \underset{i,\widehat{w}}{}}d_i^U\left\{8G_4^{(i)}\left[\mathrm{\Delta }x(\widehat{w})^{(O,U)}\right]+\left[12(1)^N\right]G_4^{(i)}\left[\mathrm{\Delta }x(\widehat{w})^{(C,U)}\right]\right\},`$ (45)
where $`G_d^{(i)}(x)`$ is the $`d`$-dimensional propagator of a scalar particle of mass ($`a=O/C`$):

$`m_i^{2(a,U)}`$ $`=`$ $`{\displaystyle \frac{4N}{\alpha ^{}}}+(1+3\delta _{a,C}){\displaystyle \frac{2P_2^2W_2^2}{\sqrt{3}}},`$ (46)
and

$`\left(\mathrm{\Delta }x(\widehat{w})^{(a,U)}\right)^2={\displaystyle \frac{\pi ^2R^2}{\alpha ^{}}}\left(w+{\displaystyle \frac{1}{2}}\right)^2+(1+3\delta _{a,O})\pi ^2\left[{\displaystyle \frac{S^2}{\alpha ^{}}}n^2+{\displaystyle \frac{2P_1^2W_1^2}{\sqrt{3}}}\right].`$ (47)
In eq.(45), the three terms in curly brackets are given by the annulus, Klein bottle and Möbius strip surfaces, respectively. Notice that the two labels $`O`$ and $`C`$ are needed because the winding and KK modes exchanged between two $`D`$-branes (annulus contribution) are not equal to those exchanged between $`O`$-planes (Klein bottle) or an $`O`$-plane and a $`D`$-brane (Möbius strip). More precisely, $`O`$-planes couple only to even winding modes along the longitudinal directions (second torus) and this explains the factor $`(1+3\delta _{a,C})`$ in (46). Similarly, with our choice of $`D/\overline{D}5`$-brane positions, closed strings exchanged between two branes have an integer winding mode along the first torus and the $`S`$ direction, whereas half-windings appear between $`O`$-planes. This explains the factor $`(1+3\delta _{a,O})`$ in (47). The degeneracy is given by $`d_i^U=d_N^U`$, with $`d_N^U`$ defined in (40).

#### Twisted string contribution

As before, this is given by the sum of the Klein bottle, annulus and Möbius strip amplitudes. It is given by

$`\rho _4^{(T,t)}={\displaystyle \frac{4\pi }{3}}{\displaystyle \underset{i,\widehat{w}}{}}d_i^T\left\{2G_2^{(i)}\left[\mathrm{\Delta }x(\widehat{w})^{(O,T)}\right]+\left[12(1)^N\right]G_2^{(i)}\left[\mathrm{\Delta }x(\widehat{w})^{(C,T)}\right]\right\},`$ (48)
where the index $`i`$ runs over $`N`$, the string oscillator number, and hence

$`m_i^{2(a,T)}`$ $`=`$ $`{\displaystyle \frac{4N}{\alpha ^{}}},`$ (49)
independent of $`a`$, whereas $`w`$ includes the winding modes in the second torus:

$`\left(\mathrm{\Delta }x(\widehat{w})^{(a,T)}\right)^2={\displaystyle \frac{\pi ^2R^2}{\alpha ^{}}}\left(w+{\displaystyle \frac{1}{2}}\right)^2+(1+3\delta _{a,O}){\displaystyle \frac{\pi ^2S^2}{\alpha ^{}}}w^{\mathrm{\hspace{0.17em}2}}.`$ (50)
The degeneracy is given by

$`{\displaystyle \underset{i}{}}d_i^Te^{4\pi il}=3{\displaystyle \frac{\theta _2^2}{\eta ^6}}{\displaystyle \frac{\theta \left[\genfrac{}{}{0pt}{}{1/6}{0}\right]}{\theta \left[\genfrac{}{}{0pt}{}{1/6}{1/2}\right]}}{\displaystyle \frac{\theta \left[\genfrac{}{}{0pt}{}{1/3}{0}\right]}{\theta \left[\genfrac{}{}{0pt}{}{1/3}{1/2}\right]}}.`$ (51)
All the considerations performed after (47) apply also here, in order to understand the form of (48) and (50).

It is interesting to note that the transverse contribution is always negative, due to the fact that the Klein bottle and annulus amplitudes are negative and dominant over the positive Möbius strip amplitude.

### 5.3 Behavior of the cosmological constant

The full 4D cosmological constant is given by adding together all the contributions, coming from both the longitudinal and the transverse SS breaking sectors:

$$\rho _4=\rho _4^{(C,l)}+\rho _4^{(O,l)}+\rho _4^{(U,t)}+\rho _4^{(T,t)}.$$
(52)
For generic values of $`\mathrm{Tr}\gamma _\beta 0`$, the longitudinal and transverse contributions, both negative, result in a monotonic $`\rho _4<0`$, leading the system towards the tachyonic instability. The same instability affects the IIB orbifold model before the $`\mathrm{\Omega }`$-projection, as can be seen by studying the torus contribution. Things are more interesting when $`\mathrm{Tr}\gamma _\beta =0=\mathrm{Tr}\gamma _\theta ^2\gamma _\beta `$. In this case, the $`D9`$ annulus contribution vanishes and the longitudinal contribution to $`\rho _4`$ is greater than zero and can thus compensate the always negative transverse contribution. This feature is characteristic of a longitudinal breaking when the SUSY-breaking operator is of order even. In presence of operators of order odd, as in section 4, it is not possible to take such a choice for the twist matrices. For this reason, although we have not performed a detailed analysis, we expect that 4D models with odd SS SUSY twists, such as the model constructed in , will unavoidably have $`\rho _4<0`$ $`R`$, and end up in the tachyonic regime. In the following we thus focus to the case in which both $`\gamma _\beta `$ and $`\gamma _\theta ^2\gamma _\beta `$ are taken to be traceless.

In Fig. 4 we present the behavior of the one-loop cosmological constant as a function of the SS radius where the other radii has been fixed to $`S=P_1=P_210`$ in units of $`\alpha ^{}`$. We found a minimum for $`\rho _4`$ which is essentially due to a compensation between the longitudinal (positive) and transverse (negative) monotonic contributions.

In order to get a better understanding of the fate of this minimum as the other moduli are varied, we study the behavior of $`\rho _4`$ as a function of the two parameters $`R`$ and $`S=P_1=P_2`$. The numeric result (figure 5) shows that the structure of the minimum in $`R`$ is still preserved as long as $`S>8`$ but $`\rho _4(R,S)`$ drives both moduli to larger values and hence to a decompactification limit.

Other cases similar to the latest one have been considered by studying $`\rho _4`$ as a function of two moduli, one of them being the SS circle radius $`R`$ and the other one the value of one or more of the extra moduli. In all the studies the scenario with a minimum in the $`R`$ direction is preserved, whereas the runaway behavior depends on how the extra moduli are fixed. As an example, fixing all the moduli in the first and second torus, $`\rho _4(R,S)`$ still develops a minimum in $`R`$ but the behavior along $`S`$ drives the system to smaller radii and thus towards the tachyonic instability.

## 6 Conclusions

In this paper, the one-loop instabilities of a certain class of IIB orbifold and orientifold string models have been analyzed. We have shown that typically in these models the SS direction tends to shrink and to reach the tachyon instability, as in the early work of . Only for orientifolds with a $`𝐙_2`$ SS twist and with a proper choice of Chan-Paton twist matrices this situation is modified. In this case the tachyonic instability is avoided, but the radius increases with a runaway behavior towards the decompactification limit in which SUSY is restored. When more geometric moduli are involved, such as in the 4D model of section 5, the situation is more interesting but less clear.

Our results are preliminary in various respects. First of all the issue of moduli stabilization should be considered for all moduli at the same time, both geometrical and not. This approach is clearly tremendously hard unless some other arguments, like that in , can predict the behaviour of the vacuum energy density as moduli are varied666See also the recent interesting proposal to stabilize moduli in string theory .. In addition, all our results hold at one-loop level only and thus can be spoiled by higher loop corrections. It has been shown in , for instance, that higher loop corrections to the vacuum energy density may lead, for large $`R`$, to $`\mathrm{log}R`$ dependencies, in addition to the usual $`1/R^n`$ terms777Such terms seem to appear also at one-loop level in certain string models with uncancelled local tadpoles for $`R\mathrm{}`$ .. Nevertheless, we think that our results indicate that the issue of moduli stabilization is somehow more interesting for models admitting a non-SUSY $`𝐙_2`$ action and for which the corresponding $`\rho (R)>0`$ for large $`R`$. Interestingly, models with a $`𝐙_2`$ non-SUSY action and hence anti-periodic fermions in the bulk are affected by possible semi-classical instabilities at strong coupling, where space-time is eaten by a bubble of nothing .

In the context of closed string tachyon condensation, it would be interesting to see if the jump in the vacuum energy $`\rho _7`$ when the $`𝐙_N𝐙_N^{}`$ ($`N^{}<N`$) transition takes place is exactly accounted for the tachyon condensation. In this case, assuming that no energy in form of radiation is released in the transition (like in the above case of the semi-classical false vacuum decay) and knowing that the final stage is flat space-time vacuum, one could have some hint on the value of the potential of some twisted closed string tachyon.

Acknowledgements

We would like to thank M. Blau, D. Perini and C.A. Scrucca for useful discussions. This work was partially supported by the EC through the RTN network “The quantum structure of space-time and the geometric nature of fundamental interactions”, contract HPRN-CT-2000-00131. M. B. thanks CONACyT (Mexico).

## Appendix

In the following appendix we briefly discuss the unfolding technique (UT) used in this paper to calculate the torus contribution to the vacuum energy density. This technique was first introduced in string theory in the context of strings at finite temperature . It has also been efficiently used in the context of threshold corrections to gauge couplings in heterotic string theories (see e.g. , and more recently ).

One loop closed string amplitudes on the torus always involve an integral over the fundamental region $``$ which is given by the complex upper half plane $`\tau _2>0`$ modded out by the modular group $`PSL(2,)`$. Analytic integration over this region is difficult to perform even for zero-point amplitudes, i.e. in computing the partition function of a model. The UT gives a systematic procedure that allows one to unfold $``$ into the strip $`S`$: $`\tau _2>0`$, $`1/2\tau _11/2`$. In many instances, such as in the models studied in this paper, such technique is essential to be able to compute the torus partition function.

The integrand $`Z_T`$ of a torus partition function is generically given by a sum of terms. Although $`Z_T`$ is always modular invariant in any consistent string theory, each term in the sum generically is not, being mapped to another term under the modular group. It is always possible, however, to choose a set of representative terms $`Z_T^i`$ having as an orbit under the action of the modular group the whole $`Z_T`$. We can thus integrate a given set of representatives of the orbit over the unfolding of $``$, namely the strip $`S`$, instead of integrating the whole $`Z_T`$ over $``$:

$`{\displaystyle _{}}Z_T={\displaystyle _{S/G}}{\displaystyle \underset{gG}{}}{\displaystyle \underset{i}{}}g\left[Z_T^i\right]={\displaystyle _S}{\displaystyle \underset{i}{}}Z_T^i.`$ (53)
As an example, we consider the torus amplitude of the 9D model of section 3, and show in detail how the UT allows to rewrite eq.(15) as (17). As can be easily shown, each of the three terms in curly brackets in (15) is not modular invariant. The first term is invariant under $`\mathrm{\Gamma }_0[2]/𝐙_2`$ which is given by the set of unimodular transformations

$`\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)`$ (56)
with $`c`$ even and $`d`$ odd, where $`𝐙_2`$ acts as $`a,b,c,da,b,c,d`$. The two remaining terms can be reached from the second one by the action of $`S`$ and $`ST`$, where $`S(\tau 1/\tau )`$ and $`T(\tau \tau +1)`$ are the usual $`SL(2,)`$ generators. Moreover, $`PSL(2,)=\{I,S,ST\}\times \mathrm{\Gamma }_0[2]/𝐙_2`$. The lattice contribution of the first term in (15) can be written as:

$`{\displaystyle \frac{\sqrt{\alpha ^{}\tau _2}}{R}}\mathrm{\Lambda }_{n,w}()^n`$ $`=`$ $`{\displaystyle \underset{n,w}{}}e^{\frac{\pi R^2}{\alpha ^{}\tau _2}\left|n+1/2+w\tau \right|^2}={\displaystyle \underset{n\mathrm{odd},w\mathrm{even}}{}}e^{\frac{\pi R^2}{4\alpha ^{}\tau _2}\left|n+w\tau \right|^2}`$ (59)
$`=`$ $`{\displaystyle \underset{p\mathrm{odd}}{}}{\displaystyle \underset{\begin{array}{c}c\mathrm{odd},\\ d\mathrm{even}\end{array}}{}}e^{\frac{\pi R^2}{4\tau _2\alpha ^{}}p^2\left|c+d\tau \right|^2},`$
where we have performed a Poisson resummation and then the change of variable $`n=pc`$, $`w=pd`$ with $`c0`$ and $`p`$ the minimum common divisor between $`n`$ and $`w`$.

By applying $`\stackrel{~}{\mathrm{\Gamma }}_0[2]=\mathrm{\Gamma }_0[2]/(𝐙_2\times T)`$ to $`\mathrm{exp}(\frac{R^2}{4\alpha ^{}\tau _2}p^2)`$, we obtain $`\mathrm{exp}(\frac{R^2}{4\alpha ^{}\tau _2}p^2\left|c+d\tau \right|^2)`$. Therefore (15) can finally be written as (17):

$`\rho _9^{(C)}`$ $`=`$ $`R{\displaystyle _{}}{\displaystyle \frac{d^2\tau }{4\tau _2^6}}(1+S+ST){\displaystyle \underset{g\stackrel{~}{\mathrm{\Gamma }}_0[2]}{}}g\left[{\displaystyle \frac{\left|\theta _2^4\right|^2}{\left|\eta \right|^{24}}}{\displaystyle \underset{p}{}}{\displaystyle \frac{1()^p}{2}}e^{p^2\frac{\pi R^2}{4\alpha ^{}\tau _2}}\right]`$ (60)
$`=`$ $`{\displaystyle \frac{R}{\sqrt{\alpha ^{}}}}{\displaystyle _S}{\displaystyle \frac{d^2\tau }{4\tau _2^6}}{\displaystyle \frac{\left|\theta _2^4\right|^2}{\left|\eta \right|^{24}}}{\displaystyle \underset{p}{}}{\displaystyle \frac{1()^p}{2}}e^{p^2\frac{\pi R^2}{4\alpha ^{}\tau _2}},`$
where $`S`$ is the strip $`\tau _2>0`$, $`1/2\tau _11/2`$, the unfolded region of $``$.