On exponential suppression of the cosmological constant in non-SUSY strings at two loops and beyond

Two independent criteria are presented that together guarantee exponential suppression of the two-loop cosmological constant in non-supersymmetric heterotic strings. They are derived by performing calculations in both the full string theory and in its underlying effective field theory, and come respectively from contributions that involve only physical untwisted states, and contributions that include orbifold twisted states. The criteria depend purely on the spectrum and charges, so a model that satisfies them will do so with no fine-tuning. An additional consistency condition (emerging from the so-called separating degeneration limit of the two-loop diagram) is that the one-loop cosmological constant must also be suppressed, by Bose-Fermi degeneracy in the massless spectrum. We remark that theories of this kind, that have continued exponential suppression to higher orders, can form the basis for a string implementation of the"naturalness without supersymmetry"idea.


Introduction
There has been interest recently in non-supersymmetric string theories, in which one might build the Standard Model (SM) directly. One particular object of focus has been the partial solution of the instability problems that generally arise in the absence of space-time supersymmetry (SUSY). In refs. [1][2][3] it was pointed out that a natural starting point for non-supersymmetric strings is a certain set of Scherk-Schwarz (SS) string models that have accidental Bose-Fermi degeneracy in their massless spectra. It is important to realise that this is a possibility even if the effective theory is entirely non-supersymmetric, and indeed such models were explicitly constructed in the string theory in [1]. They have a visible spectrum that resembles that of the SM and a hidden sector whose Bose-Fermi non-degeneracy is equal and opposite to that of the SM, as shown in Figure 1. In theories of this kind, successive Kaluza-Klein (KK) levels are unable to contribute to the one-loop cosmological constant, which can only get contributions from heavy winding modes, string excitation modes and also from non-level matched states. As these modes are all short-range, they are unable to explore the whole compact volume. Consequently, even if the compactification scale is only moderately large, their contribution to the cosmological constant (and hence destabilising dilaton tadpoles) is parametrically exponentially suppressed. As an example, a supersymmetry breaking scale of say 1/R ∼ 10 14 GeV requires a string mass of M s ∼ 10 16 GeV to get the correct Planck scale (where R represents the compactification scale and where the gravitino mass goes like 1/(2R)). Even though the radius is then only 10 2 string lengths and the visible spectrum entirely non-supersymmetric, the cosmological constant is suppressed by an astronomical factor, e −4πRMs ∼ 10 −546 . (It is worth adding that more generally such a configuration seems to be the only way to get an effective 4D non-supersymmetric theory. A generic non-supersymmetric KK construction will be unstable and collapse on timescales of order R, thus it is never really four dimensional. Such a cosmological constant, generated entirely by heavy modes, also allows novel separations of finite UV and IR contributions to the potential [4]. ) An open question is what happens at two-loops and beyond in such theories. Does the exponential suppression continue? Field theory intuition says that generic two-loop contributions will start to make their appearance, but it is conceivable that some kind of string "miracle" appears to save the day, or that a further subset of one-loop suppressed theories may have two-loop suppression in the cosmological constant as well. This paper shows by explicit calculation that (while we cannot rule out the former) the latter is highly likely. We derive two criteria that define a sub-class of theories which continue to enjoy exponential suppression at two-loops. Like the one-loop case, this suppression is simply an Figure 1: The spectrum of theories that satisfy Bose-Fermi degeneracy with a Standard Model-like light sector (reproduced from [1]). As the Standard Model does not have Bose-Fermi degeneracy a cancelling hidden sector is inevitable, but note there is no supersymmetry in the spectrum. Models off this kind were constructed in [1]. accidental consequence of their particle content.
We should note that from a field theory point of view the possibility of such cancellation is surprising. At one-loop, some scalars and some fermions inevitably gain supersymmetry breaking masses of order the KK scale even though some others may remain exponentially light [1]: this then feeds back into the cosmological constant which one generically expects to be of order the KK scale suppressed by a two-loop factor. However the two-loop cosmological constant calculation in the full string theory already contains within it the one-loop shifts in the spectrum. Therefore another way to view the criteria we present, is that they precisely determine when the latter conspire to cancel in the two-loop cosmological constant.
Our programme, and this entire approach, is reminiscent of the field theory ideas of refs. [5][6][7][8] which attempt to achieve naturalness without supersymmetry, by essentially extending the Veltman condition of ref. [9] to all orders. Indeed, it is a remarkable fact that, thanks to the theorem of Kutasov and Seiberg [31], non-supersymmetric string theories with D = 4 whose cosmological constant vanishes at one-loop must also satisfy the "field independent" Veltman condition, namely Str(M 2 ) = 0 [10,11]. Hence although the object of study here is the cosmological constant, not the mass of some putative Higgs, there is a direct link. However the string case goes even further than the field theory one: there are no freely adjustable couplings, since couplings are all either zero or one (or themselves exponentially suppressed by the volume), so there is absolutely no fine-tuning involved. Theories either have the correct massless particle content or they do not. At the one-loop level, because of this connection to the Veltman condition, any model with vanishing cosmological constant can be thought of as a stringy UV completion of the scenario outlined in ref. [8]. Although we stress that the operator being considered here is the cosmological constant, the exact same procedure could be carried out for the Higgs mass-squared itself. This is discussed in more detail in ref. [12]. In the models of ref. [1], this is achieved because a Scherk-Schwarz deformation preserves the Bose-Fermi degeneracy of the massless modes in all of their KK levels as well. In the logarithmically running low energy theory, one then assumes that the relevant scale at which such a relation should be applied is the compactification scale, above which the theory becomes extra dimensional. An important difference though is the motivation for imposing the condition at that scale which has nothing to do with SUSY being restored there, but rather the one-loop cosmological constant vanishing 3 . At the two-loop level, we will find as mentioned two rather different looking criteria for vanishing cosmological constant. The criterion for the vanishing of the entirely untwisted contributions (that is diagrams whose propagators contain only the descendants of broken N = 2 supermultiplets) is a complicated combination of parameters (numbers of gauge bosons, gauginos, hypermultiplets and so forth) that essentially counts the two-loop effective field theory divergences. As we will demonstrate, this parameter is most easily extracted from the constant term in the "q-expansion" of the two-loop string partition function. By contrast diagrams that contain twisted loops (that is loops of twisted states that still appear in complete N = 1 chiral supermultiplets) can vanish due to the cancellation of combinations of "field dependent" Veltman conditions. Such diagrams have a different dependence on the volume modulus from the entirely untwisted ones, so to avoid fine-tuning one has to impose a second independent criterion for the twisted states, of the form U (−1) F U Tr|Y U T T | 2 = 0 where U stands for generic untwisted fields in the theory, and the trace is over the pairs of twisted states to which they couple, with tree-level coupling Y U T T . This criterion is quite Veltman-like, but note that it is the sum over the Veltman conditions of all the twisted states that appears; we do not need to apply them individually. Furthermore the couplings are degenerate, so again the vanishing of this quantity is a question of particle content.
An important aspect to bear in mind is that one requires an absence of gravitationally coupled products of one-loop divergences in order to produce the above criteria. This contri-bution would normally come from the so-called separating degeneration limit of the two-loop partition function, which we will discuss in some detail. Such terms are absent only if one has chosen a theory that already satisfies the criterion for the one-loop cosmological constant to vanish, namely massless Bose-Fermi degeneracy, N Indeed, more generally one can see that at each order, a sensible criterion for continued suppression can only be achieved when the criteria for all the orders below are satisfied.
The results we have found are a natural extension of this work, which leads one to speculate on the existence of three-loop and beyond cancellations, and whether there might be a universal condition for string theories that, like the one conjectured for field theory in ref. [5], ensures cancellation to all orders. Conversely, it raises the possibility that imposing the requirement of continued exponential suppression to ever higher order could give interesting predictions for the particle content of the theory.  [2,3,38,[92][93][94][95][96][97][98][99][100][101][102]. However, care is required from the outset as there are possible pitfalls. In particular, one of the major difficulties in calculating string amplitudes beyond one-loop proved to be the integrating out of the supermoduli. If done incorrectly, computations of this type typically give ambiguous results that depend on the choice of gauge. For example, attempts were made in the past to determine the value of the two-loop vacuum amplitude for the non-supersymmetric models presented in refs. [2,3] (the so-called KKS models). The initial claim was that the cosmological constant is vanishing, but contradictory evidence was presented in ref. [38]. In fact both of these results suffered from the aforementioned issue of gauge dependence. A correct gauge-fixing procedure was later introduced in the work of refs. [103][104][105][106], and the computation was re-done in ref. [107] with the conclusion that the two-loop contribution is indeed non-vanishing for the KKS models. It is these later papers that form the basis of our analysis.
For the type of non-supersymmetric model described in ref. [1], one does not actually expect the two-loop contribution to the cosmological constant to be identically zero. As described in the Introduction, the best one can achieve at one-loop is for it to be exponentially suppressed if the massless spectrum contains an equal number of bosons and fermions. Therefore we seek a similar suppression at higher loop order.
Note that as the main source of the cosmological constant (a.k.a. Casimir energy) in large volume Scherk-Schwarz compactifications is the massless spectrum, one might think it is preferable to approach the entire problem from the perspective of the effective field theory. However at two loops, it is not always obvious how the string computation factorises onto the field theory diagrams. In addition one would have to perform an analysis in the effective softly broken supergravity, and there are certain purely string contributions, in particular the separating degeneration limit (of which more later), that one has to check. These issues are exacerbated by the fact that the string models typically have a large rank making it tedious to count states, and by the fact that one would in any case have to determine all the tree-level couplings of the effective field theory. As we shall see, it is by contrast far easier to simply extract the coefficient of the relevant (constant) term from the q-expansion of the two-loop partition function.
The structure of two-loop superstring amplitudes is built upon the representation of the worldsheet by a super Riemann surface of genus two. Let us start with a brief outline of the essential properties of such surfaces, and as a warm-up exercise then perform the computation of the two-loop cosmological constant in an entirely supersymmetric theory.
Consider a super Riemann surface of genus g with a canonical homology basis of A I and B I cycles as shown in Figure 2. The period matrix Ω IJ is given by holomorphic abelian 1-forms w I dual to the A I -cycles such that (2.1) In addition to the period matrix there is the super period matrix,Ω IJ , which can be defined in a similar way, by integrating superholomorphic 1/2 forms over the A I and B I cycles.
The supermoduli space M g of a genus g super Riemann surface contains 3g − 3 even moduli and 2g − 2 odd moduli for g ≥ 2. Specialising to the case where g = 2, the super period matrix gives a natural projection of the supermoduli space of a super Riemann surface onto the moduli space of a Riemann surface, and its 3 independent complex entries provide complex coordinates for the moduli space of even moduli, M 2 . The super period matrix can be expressed in a simple way in terms of the period matrix and, following the procedure of refs. [103][104][105][106], one can work in the so-called split gauge, which has the main advantage that the period matrix and super period matrix are equivalent, and one can simply use Ω IJ to denote both. It can be parametrized by where τ 11 , τ 12 and τ 22 are the complex variables corresponding to the three moduli (i.e. playing the same role as τ in the one-loop diagrams). To make the discussion widely accessible, we present the result (which derives from refs. [103][104][105][106] after some work and carefully accounting for the measure) in terms of two-loop ϑ-functions, the most natural extension of the standard one-loop formalism. For a genus 2 surface there are 16 independent spin structures, labelled by half-integer where κ is a 2-vector of spin structures on the A I -cycles, and κ is a 2-vector of spin-structures on the B I -cycles.
The ϑ-functions with characteristic v are defined by A given spin structure is said to be even or odd depending on whether 4κ · κ is even or odd. For vanishing characteristics, v = 0, all of the 6 odd spin-structure ϑ-functions are identically zero (much like ϑ 11 in the one-loop case), so that (2.5) The even spin structures will be denoted generically with a δ, and the even ones with a ν: for example even ϑ-functions will be written as ϑ[δ].
After integrating over the supermoduli, enforcing the GSO projection and summing over spin structures, the cosmological constant for the supersymmetric heterotic string can be written [103][104][105][106] , and the integration is over the fundamental domain of the moduli, F 2 , typically taken to be [108][109][110] 1.
The modular forms appearing in eq.(2.6) are defined as follows. First it is useful to define This expression uses the fact that any even spin structure can be written as the sum of three odd spin structures, δ = ν 1 + ν 2 + ν 3 ; in the sum, ν 4,5,6 are the remaining three odd spin structures, and In term of Ξ 6 we then have where the product is obviously over even spin structures only. In the end the two-loop cosmological constant in a SUSY theory is of course zero, as it should be; this is due to the genus two version of the abstruse identity, namely Υ 8 = 0.

The Scherk-Schwarzed cosmological constant
Adapting the technology of the previous section, one can now start to build up the twoloop cosmological constant for the non-supersymmetric theories of ref. [1]. These theories are constructed by taking a 6D theory in the free fermionic formulation and compactifying down to 4D on a T 2 /Z 2 orbifold, breaking spacetime supersymmetry through a coordinate dependence in the compactification (CDC). This is the equivalent of the Scherk-Schwarz mechanism in string theory. Sectors that are twisted under the final orbifolding remain supersymmetric under the deformation, and so their spectrum is unchanged. (Whenever we refer to "twisted" or "untwisted" this will always mean with respect to the final orbifolding.) At genus two there can be a twist associated with each loop, but the focus will mainly be on the totally untwisted sectors since twisted states are involved in a very restricted set of diagrams due to their remaining supersymmetric structure. It is worth elaborating on this last particular aspect before we start the calculation of the totally untwisted diagrams in earnest. One can proceed by constructing an extension of the argument of refs. [1,16]. At one-loop the partition function of the N = 0 deformed theory (whose orbifold action we shall denote by g) is decomposed as where the indices represent the orbifold action on the A and B cycle. The Scherk-Schwarz phases on the world-sheet degrees of freedom are denoted by a vector e. The only dependence on them is in the first totally untwisted term. The second term is (up to the factor of 1/2) the partition function of the non-orbifolded and non-deformed N = 2 theory, while the second line is the partition function of an entirely undeformed N = 1 theory; both are zero, and hence only the first term can give a non-zero contribution to the cosmological constant.
(So for example any N = 2 e → N = 0 un-orbifolded theory with Bose-Fermi degeneracy implies the existence of a chiral orbifolded N = 1 e → N = 0 theory that also has Bose-Fermi degeneracy.) Continuing to two loops, a similar decomposition would look like 4Z(e) = Z 00 00 (e) − Z 00 where now of course there are two cycles. The bracket is the undeformed N = 1 theory and must vanish by supersymmetry, and the first term is the partition function for the entirely unorbifolded theory, representing contributions containing the untwisted fields only. Clearly the one loop argument would go through as before, were it not for the additional e-dependent terms on the first line, which represent diagrams that have twisting on one pair of A I , B I cycles, with the other pair of A I , B I cycles remaining entirely untwisted. Such diagrams will be referred to as "mixed" diagrams. What remains is therefore to determine the contributions of the mixed diagrams at leading order, and the contribution from the entirely untwisted first term, Z 00 00 (e). It is these two different kinds of contribution that lead to the two criteria mentioned in the Introduction.
The former will be dealt with explicitly later, but for the moment let us now turn to the calculation for the entirely untwisted contribution which is (up to a factor) the cosmological constant of the un-orbifolded theory. To define the sums over spin structures, the CDC and vector notation is the standard one, summarized in ref. [1]. In particular dot-products are the usual Lorentzian ones, while a separate sum over basis vectors V a is understood; thus explicitly the collection of spin-structures in a particular sector are α I V ≡ α I a V a and β I V ≡ β I a V a , with a labelling the basis vectors and, recall, I = 1, 2 labelling the A I and B I cycles. The right-and left-moving fermions have spin-structures denoted The primes represent the shift due to the CDC deformation, that is where n I = n 1I + n 2I , I = I 1 + I 2 and n iI are the winding numbers and I i are the dual-KK numbers in the Poisson resummed theory. In the present context, there are 16 transverse right-moving real fermions and 40 transverse left-moving real fermions on the heterotic string (so that S R/L are vectors containing 16 and 40 different spin structures respectively).
Let us describe the factors in detail. In addition to the self-evident fermion factors, the compactification from 6D to 4D has introduced an extra factor of the two-loop Narain partition function for the two compact bosonic degrees of freedom, Γ 2,2 . In its original non-Scherk-Schwarzed and un-Poisson resummed format it would look like and where G ij and B ij are the usual metric and antisymmetric tensor respectively. After introducing the CDC shift and performing a Poisson resummation on all of the m's, it takes the form (2.18) We should point out that in the above equations and in what follows, we have lowered the 'i' index on the winding numbers purely to simplify notation; they have not been lowered through the use of the metric G ij . A word of warning is also required concerning the definition of the {α a , β a } summation in eq.(2.14): the partition function Γ 2,2 is of course a function of I i , n I i , but now so are the S L and S R due to the CDC induced shift. Therefore one cannot really factor the summations as we appear to do above: everything to the right of Γ (2) 2,2 is to be correctly included in the sum over I i , n I i . However the case of ultimate interest is when the radii are moderately large, since as described in the Introduction we wish to determine the presence or otherwise of unsuppressed SS contributions to the vacuum energy. These can only correspond to n I = 0 mod (2) as is evident from eq.(2.17), while we require at least one of the I=1,2 to be equal to 1 mod (2) to avoid cancellation by supersymmetry. The Poisson resummation could have been done for different choices of the I separately but it would amount to the same result. The result is leading terms that carry the usual volume dependence but are otherwise not suppressed. Conversely the sub-leading terms coming from the non-zero n I modes would involve a simple generalisation of the saddlepoint approximation used for the one-loop case in ref. [1] leading inevitably to exponential suppression.
The final ingredients in eq.(2.14) are the GSO projection phases,C α β . These can be deduced from the fact that two-loop partition functions factorize onto products of two oneloop partition functions in a certain limit of moduli space, at which point the GSO coefficients must factorize as well [37,111]. Since the GSO coefficients are completely moduli independent, this factorization must hold everywhere. They can therefore be written as a product of the known genus one coefficientsC As described in ref. [112], most generally these are functions of the structure constants k ab , k eb , k ae and k ee , that take the following form with the vector e assuming a projective role, completely analogous to that of the other basis vectors. For the canonical assignment of structure constants for the CDC vector e, there is no sector dependence in the phases, that is However, note that in eq.(2.14) we actually haveC α β rather thanC α β . This primed definition does not include the factors of exp[2πi(α I a s a + β I a s a )] appearing in the above equations, which are effectively contained within Ξ 6 instead.
Eq.(2.14) is the "master equation" that provides our first criterion. It is straightforward to check that it has the correct modular properties under Sp(4, Z) by considering the transformations given in eq.(B.3). As we are about to see, one can also use it to determine the leading contribution to the cosmological constant by deduce the q-expansions, by inserting the explicit expressions for the two loop ϑ-functions, in Appendix A. Writing the cosmological constant as vanishes. Note that ℵ is a product of the measure and the partition function.

The q-expansion of ℵ
Let us proceed to examine the q-expansions for the cosmological constant in certain limits, in particular the large radius limit. The general form of the integrand in the two-loop cosmological constant is where the variables L 1 , L 2 , L 3 can be interpreted as Schwinger time parameters for the three propagators of the two-loop sunset Feynman diagram shown in Figure 3. With this parametrization, det Im(Ω) = L 1 L 2 +L 2 L 3 +L 1 L 3 , and the fundamental domain F 2 restricts the variables so that 0 < L 2 ≤ L 1 ≤ L 3 . By parameterizing the period matrix in this way, the q I -expansion of ℵ is symmetric with respect to the three q I . It can be relatively straightforwardly evaluated. The q-expansion of Ψ −1 10 is given by 5 The rest of ℵ is model dependent and can be determined using the q I -expansions of the ϑ-functions in Appendix A. Figure 3: Generic sunset diagram for the two-point function.
As an example of the whole procedure we will consider an SO(10) model that has massless Bose-Fermi degeneracy, and hence exponentially suppressed cosmological constant at oneloop. The model is presented in Appendix C, where it is shown explicitly that in the SUSY theory (i.e. the theory without any CDC deformation) the two-loop cosmological constant vanishes. It is also shown there that the one-loop cosmological constant in the broken theory is exponentially suppressed because there is Bose-Fermi degeneracy at the massless level, and hence the constant term in the one-loop partition function is absent.
Recall that non-vanishing two-loop contribution to the cosmological constant comes from sectors in which at least one of 1 and 2 is equal to 1 mod (2). For example, if 1 = 2 = 1, the q-expansion of ℵ in the full non-SUSY SO(10) theory is found to be (2.27) The terms with 1 = 1 and 2 = 0, and with 1 = 0 and 2 = 1 have the coefficients of 1/q i permuted but are otherwise identical. In particular the constant term is the same. In total then, we find a non-vanishing constant piece, and conclude that this particular model gets a generic (i.e. not exponentially suppressed) contribution to the cosmological constant starting at two-loops.

Field theory factorization: identifying leading contributions
Note that the constant piece in ℵ includes various field theoretical contributions, not only the ones corresponding to the sunset topology. For reference the contributions in the field theory are displayed in figure 4 in the parent N = 2 formalism. They can in principle be computed in the 6D field theory following ref. [113]. Given the complexity of the theories involved, and the fact that one would have to determine the spectrum and all the effective couplings, this would be an extremely arduous task, and it is actually much easier to simply determine the two-loop partition function directly as above. Nevertheless it is instructive to see how the expression of eq.(2.27) does indeed give the corresponding field theory contributions in the various degeneration limits. First note that for sufficiently large compactification volume the non-zero winding mode contributions are extremely exponentially suppressed compared to those with n I i = 0. In addition the supersymmetric minimum for the CDC deformations is around U 1 = 1 as discussed in ref. [4]. Expanding around this point and using eq.(2.17), the dominant contributions to the cosmological constant are given by In the regions of the fundamental domain in which the real parts of the three moduli are integrated from −1/2 to 1/2, the only non-zero contributions come from the physical states with a i = b i ≥ 0, and are given by the physical coefficients C aa . (This result is also a consequence of the fact that modular invariance requires a i − b i ∈ Z.) The approximation sign is there because, as was also the case for one-loop integrals, there is a small region of the fundamental domain for which the integration over the real parts of the moduli does not extend over the full domain −1/2 < Re(Ω IJ ) ≤ 1/2. In this region, there is no level-matching and so unphysical states contribute to the vacuum amplitude. Nevertheless as in ref. [1], we find that the contributions from these unphysical states are also extremely exponentially suppressed compared to the both the massless contributions and the lowest lying string excitation mode contribution, provided that the compactification radii are sufficiently large. As per the previous subsection we are therefore interested in the value of C 00 , the coefficient of the constant piece giving leading order contributions. The important observation is that for these massless modes (with a 1 = a 2 = a 3 = 0) the expression in eq.(2.28) has simply degenerated to the 4 dimensional field-theory result in the Schwinger formalism, so the coefficient C 00 could also be calculated in the effective 6D→4D Scherk-Schwarz field-theory. The relevant diagrams are shown together with the coefficients of their contribution to C 00 in figure 4, which are deduced from the calculations in ref. [113]. (Note that all coefficients are written for the fields as they decompose into boson or fermionic components of N = 2 multiplets.) Different limits of the integral in eq.(2.28) generate all the field-theory diagrams in figure  4. In particular the "double-bubble" diagrams come from the region where L 1 , L 3 → ∞, while L 2 1. Explicitly in this limit, one still requires a 1 = a 2 = a 3 = 0 to avoid exponential suppression, but can everywhere replace det Im Ω ≈ L 1 L 3 . The L 2 integral then may be trivially performed (with its upper limit L 1 being effectively infinite). Taking for example 2 2 = 1 1 = 1 in this limit results in an integral proportional to which (taking the upper limit L 3 → ∞ on the L 1 integral) has the form of a product of two one-loop Poisson resummed Schwinger integrals in a KK theory with two extra dimensions. A more complete way to reach this conclusion would be to first go to the "non-separating degeneration" limit of ref. [106], i.e. τ 22 → i∞ with τ 11 , τ 12 fixed, and from there take τ 11 → i∞. The field theory recipe for evaluating C 00 for the Scherk-Schwarzed string theories is therefore as follows: retain in the list of two-loop diagrams only those that are exactly massless, meaning that the states on all propagators do not receive any CDC shift. Then C 00 is precisely twice the resulting sum of coefficients.
The reasoning is straightforward and exactly mirrors what happens in the one-loop case. First recall that we are (for this calculation) considering only untwisted states in the diagrams of figure 4. This implies that there is KK and e charge conservation at the vertices, which in turn implies that the CDC shifts pairs of either Fermion-Fermion or Boson-Boson masses on the sunset diagrams. The nett effect of such a shift is that the space-time statistics of an entire loop on the diagram is reversed, and consequently these diagrams contribute with an additional minus sign. Meanwhile the "superpartner" diagram (in which the space-time statistics really is reversed on that loop) is still present: hence a factor of two.
In principle the sum of coefficients can vanish, and the important aspect that makes this possible is the coupling degeneracy, which is due to the underlying supersymmetry of the undeformed theory, and the N = 2 structure of the untwisted (i.e. SUSY breaking) sector. This is a well-known feature of effective string theories, but the crucial point here is that while at the level of the field theory a complete cancellation of contributions may seem like a miraculous tuning, at the level of the string theory it is merely a consequence of the particle content and the corresponding partition function and measure (and indeed there are no independent couplings). It is worth repeating that from this point of view (and in practice), it is far easier simply to work with the q-expansion of the string partition function, than to attempt to evaluate C 00 for the entire field theory.

The separating degeneration limit
There is one limit that would not be covered by the field theoretic treatment described in the previous sub-section, namely the separating degeneration limit. For a two-loop string vacuum amplitude this corresponds to taking the limit τ 12 → 0 keeping τ 11 , τ 22 fixed. This gives a Riemann surface that looks like two one-loop vacuum amplitudes connected by a long thin tube, as shown in Figure 5. The limits of various objects appearing in the two-loop cosmological constant are given by [106] where µ 1,2,3 and ν 0 are the three even and unique odd genus 1 spin structures respectively, while the genus two Narain lattice Γ 2,2 splits into a product of two genus one Narain lattices. Therefore the full two-loop cosmological constant in the separating degeneration limit takes which is essentially two one-loop vacuum amplitudes connected by a divergent propagator. We therefore make the crucial conclusion that the separating degeneration limit contains the divergence due to any uncancelled one-loop dilaton tadpoles. In general, i.e. at higher loop order, one expects such terms to always be present. That is at n-loop order, any uncancelled tadpoles from the (n − 1)-loop theory will contribute to divergences in the cosmological constant. Thus if the one-loop partition function has Bose-Fermi degeneracy, these terms are a divergence multiplied by an exponentially suppressed coefficient. One may confirm that the same conclusion is arrived at using the full q-expansion in the separating degeneration limit. First of all in this limit we have Returning to the non-SUSY SO(10) model with massless Bose-Fermi degeneracy given in Appendix C, for the untwisted sector with 1,2 odd, the leading term in the q-expansion of the partition function after summing over spin structures is given by (2.32) The constant term has vanished as expected in this limit, for this model.

Comments on the effect of the one-loop tadpole
For the class of non-SUSY string models that we are considering in this paper, it is known that at one-loop order there is an exponentially suppressed but non-zero dilaton tadpole. If this tadpole is left uncancelled, then as we saw in the previous section, it can contribute through the separating degeneration as a divergence in the two-loop cosmological constant. It is well known that infrared divergences can appear in this degeneration [114][115][116], however, our experience from QFT is that these divergences typically arise because we are asking the wrong questions. As we have learned from QFT, what one should in principle do is stabilise the theory in the correct one-loop vacuum so that the tadpole is effectively cancelled. The two loop separating degeneration divergence would then be seen to be merely an artifact that disappears if we perform this procedure. It might also be the case that one could live with the tadpole and have a dynamical cosmologically evolving background as in ref. [117]. These issues have also been discussed in refs. [118][119][120].
In generic non-supersymmetric string models the dilaton tadpoles can be large. Any attempt to cancel the tadpole through a background redefinition would require such a large shift that it is highly unlikely that the new vacuum bears any resemblance to the original, thereby negating any positive phenomenological aspects of the originally constructed model.
The key point about the specific types of models we consider here is that the dilaton tadpoles are exponentially suppressed. If one were to employ a background redefinition, the shift to achieve this should be sufficiently small so as not to result in any appreciable alteration in the phenomenological properties, including the spectrum of the massless states. If this were not the case then clearly there would be a problem, since the construction of models with suppressed cosmological constants is dependent on a careful cancellation of bosonic and fermionic massless degrees of freedom at one-loop order. In theory one is able to perform this background shift at the string theory level (see ref. [121]), however in practice this would be rather involved.
An alternative argument is built around balancing the one-loop tadpole itself against another contribution as in ref. [4] where the mechanism is incorporated in the effective supergravity theory, and of course should not itself result in a large cosmological constant. In a framework that is completely stable, where the dilaton tadpole is cancelled, the divergent contribution to the two-loop cosmological constant should then vanish, while crucially the remaining contributions remain unaltered. For the models which contain a bose-fermi degeneracy, the potential can be written as where V UV is computed in the full string theory while V IR arises from non-perturbative effects in the effective field theory. The key point is that because V UV comes from the contribution of heavy modes only, it is independent of the low-energy IR physics. Therefore, we can introduce some stabilising mechanism in the IR to cancel the UV contribution, and provided this does not alter the masses of states in any way that is not exponentially suppressed, then the massless spectrum will remain unchanged. A full treatment of the tadpole is beyond the scope of this paper and so we leave a complete study of the dynamics to future work. With this in mind, we assume it is fact consistent to study the cosmological constant in our naive vacuum, with the knowledge that the conditions on the structure of the massless spectrum that guarantee exponential suppression will still be satisfied after the shift to the correct vacuum. We emphasise that this would not be the case without exponential suppression of the one-loop tadpole. Those theories would undergo large shifts in the metric upon finding their true vacua, and any putative dilaton stabilisation would most likely be completely invalidated in the process, along with any two-loop discussion.

Suppression of the "mixed" diagrams
This completes the derivation and discussion of the first criterion for vanishing two-loop cosmological constant. It remains to consider the contributions with one untwisted propagator and two twisted ones, i.e. the mixed diagrams. In the untwisted sector, the compactification from 6D to 4D resulted in the inclusion of the two-loop Narain partition function for the two compact bosonic degrees of freedom. This term meant that, for sufficiently large compactification radii, contributions to the cosmological constant from non-level matched states (including the proto-graviton) were exponentially suppressed compared to contributions from both massless states and the lowest lying string excitation modes. By contrast, for the twisted sectors, the partition function for the two compact bosonic degrees of freedom is given by [107,122] where τ is the Prym period and where Z(Ω) is the partition function for two bosonic degrees of freedom in the uncompactified theory.
For twisted sectors involving some twist on only one of the two loops we anticipate that the cosmological constant may still receive a non-zero contribution. First we can see that again it is the massless states which provide the dominant contributions to the cosmological constant, while massive states receive exponential suppression after integrating over the real parts of the three moduli as before. The contributions from non-level matched (i.e. unphysical) states are also exponentially suppressed (for sufficiently large compactifiction radii), despite the fact these sectors do not include the two-loop Narain partition function. Instead, in these sectors there is the factor, Γ 2,2 (τ ) = which just has the form of a one-loop Narain partition function involving the Prym period τ . As usual we can perform a Poisson resummation giving In order to show that the unphysical states are suppressed even in the twisted sectors, we make use of the fact that there is a relation between the Prym period τ and the period matrix Ω. The Schottky relations state that for any i, j = 2, 3, 4 (2.38) In the notation above, for any given twist = 0, there are 6 even spin structures δ where δ + is also even. These 6 spin structures are denoted δ + i and δ − i , for i = 2, 3, 4, where δ − i = δ + i + . The region of moduli space where there is no level-matching is when L 1 , L 2 , L 3 are all sufficiently small and are at most O(1). When the imaginary parts of the three moduli are small, the Schottky relations tell us that Im(τ ) is also small (while it is large when both L 1 and L 3 are sufficiently large) and so by considering the Poisson resummed form of Γ 2,2 (τ ) we see that small values of τ result in exponential suppression.
What remains therefore are the diagrams with a twisted loop and an untwisted propagator containing only physical states. (Due to the Z 2 orbifold, there can only be either U U U or T T U vertices in the superpotential of the unbroken theory, and hence no diagrams with a single twisted propagator.) The coefficients of these diagrams can be easily evaluated in the field theory. The integral for a loop of fermions of mass m 1 and m 2 coupling to a scalar are of the form (2.39) We can assume one mass to be zero, and first consider the fermion as the KK states. Thus we have to consider the Euclideanized integrals (2.40) We also have the case where the scalar is the KK state which involve the integral (2.41) These diagrams will come with a coefficient Tr Y 2 U T T where Y U T T is the tree-level UTT Yukawa coupling in the superpotential; it takes the value √ 2g Y M or 0 depending on whether the charges are conserved at the vertex. The double-bubble diagrams (for Yukawas) will have the same coefficient with a minus sign (2.43) In the untwisted sector, it is possible to show that the sunset diagrams can be reduced to the form of scalar double-bubble diagrams by basic manipulation [113]. However, similar manipulations do not produce the same result in the twisted sectors and so we must evaluate the sunset diagrams as they are. Using the Schwinger formula Re(A) > 0, (2.44) and the integrals we find that the sunset diagrams can be written in the following form, where either m s = 0 if the single untwisted propagator is a fermion, or m f = 0 if it is a scalar: (2.46) The above integral has UV divergences when at least two of the Schwinger parameters y 1 , y 2 , y 3 approach zero. Therefore, when we come to evaluate these diagrams later we will introduce a regulator e −N 1 y 2 + 1 y 3 .

Figure 8 diagrams
We may proceed to calculate the relevant integrals in a similar manner to refs. [123,124]. In the untwisted sector the scalar figure 8 diagram is proportional to J(m 2 (2.47) We need to consider the case with two compact dimensions with radii R 1 and R 2 . We will begin by considering the supersymmetric case in order to verify cancellation between all diagrams. For the scalar mass we therefore have where m 1 and m 2 are Kaluza-Klein numbers. Therefore, again using the Schwinger formula and integrating over the momentum p we obtain To proceed with the calculation we introduce a regulator e −N/t , allowing us to interchange the order of summation and integration. From there we can perform a Poisson resummation on the KK numbers and finally obtain where U 2 = R 2 /R 1 and E(U, n) is the real analytic Eisenstein series with For a twisted loop there are no associated KK states and so we only have the contribution from the massless state. In this case we simply have J = 1 16π 2 N and so for the figure 8 diagram with a single twisted loop we find (2.52)

Sunset diagram
When the untwisted propagator in the sunset diagram is a scalar we obtain the result (2.53) where the notation E (0,1) (U, n) ≡ ∂ n E(U, n). On the other hand when the untwisted propagator is a fermion we have (2.54) Therefore the total contribution from the sunset diagrams with unbroken supersymmetry is which exactly cancels the contribution from the figure 8 diagram as expected.
Finally we can obtain the two-loop contribution to the vacuum energy from the twisted diagrams in a theory with supersymmetry broken by the Scherk-Schwarz mechanism. The masses of the twisted states themselves are unaffected by the supersymmetry breaking, but the masses of the untwisted states to which they couple may still be shifted. The result of Scherk-Schwarz supersymmetry breaking amounts to shifting the KK numbers by 1 2 . We may proceed with the calculation in the same way as before, and find the shift in the KK numbers results in a replacement of the real analytic Eisenstein series E(U, n) by E 1 (2.56) Therefore, we find the contribution from the twisted sectors to be whereẼ(U, n) is an Eisenstein series restricted to l 1 + l 2 = odd, N T is the number of twisted degrees of freedom, and N U b and N U f denote the number of untwisted bosons and fermions respectively that couple to the twisted states and whose masses remain unshifted after supersymmetry breaking. Therefore, we see that if the spectrum contains a degeneracy in the number of massless bosons and fermions in the untwisted sector that couple to twisted states, then the leading contribution from the twisted sectors is zero. Noting that the functional form of this term makes it unnatural for it to cancel against the entirely untwisted contribution, this gives us a second criterion for the vanishing of the two-loop cosmological constant: = 0 where in terms of the couplings we have and where for a given U , the coupling Y U T T is considered to be a matrix with indices running over all the twisted states, and includes both gauge and Yukawa couplings. Taking account of the degeneracy in the couplings, we can write a simple operational expression for , namely = U,T,T =massless where the sum is over all massless physical untwisted fields, and pairs of twisted fields. The δ Q -function imposes either simple charge conservation for the charge vectors of the triplet of fields (i.e. representing superpotential φψ L ψ R type couplings), or charge conservation with an extra unit in the non-compact space-time index (representing gauge A µψ L γ µ ψ L type couplings that have an extra Dirac matrix).

Conclusion
In this paper we have derived two criteria for the exponential suppression of the two-loop cosmological constant in string theories with spontaneously broken supersymmetry. These two criteria determine respectively when the leading order entirely untwisted and partially twisted contributions vanish. The untwisted criterion, in eq.(2.23), is most easily determined in any given model from the vanishing of the constant term in the q-expansion of the integrand in the two-loop cosmological constant. Note that this object contains factors from the partition function but also from the measure; the criterion can not be determined from the partition function alone. The twisted criterion can be determined from the effective field theory, but can most easily be evaluated in a very simple operational way simply with the knowledge of the states in the spectrum and all of their charges. The resulting condition, in eq.(2.58), is the vanishing of a "sum of Veltman conditions" for the twisted fields; that is, in terms of the effective field theory, one can imagine that at the one-loop level the twisted states in the spectrum will receive quadratically divergent contributions to their mass from the leading quadratic divergence in the cosmological constant. At the two-loop level, these terms will enter into "sunset" diagrams, but the degenerate nature of the couplings implies that the sum of such contributions may vanish, depending on the spectrum.
For consistency, one should also impose the vanishing of the one-loop leading contribution to the cosmological constant, which is achieved in theories that have Bose-Fermi degeneracy in their massless physical states. Divergences associated with the one-loop dilaton tadpole would appear at two loop level in the so-called separating degeneration limit of the diagrams, a limit that resembles two one-loop torus diagrams connected by a long thin tube. However, their presence does not actually affect the phenomenology of these models since the crucial point is that because the tadpoles are exponentially suppressed, their effect on the physical spectrum is in fact negligible.
The two criteria we have presented here can be thought of as a stringy implementation of the "naturalness without supersymmetry" idea first proposed in ref. [5] up to the two-loop level. The existence or otherwise of models that satisfy these conditions, and their properties should they exist, is a subject of current study, which will be reported elsewhere [125].
It would also be of interest to search for a subset of theories that mimic the supertrace rules in models involving D3-branes, where vanishing one-loop supertraces are known to extend to higher order automatically [126].

(C.3)
Beginning with a simple model defined only by the vectors V 0 and V 1 , one finds a contribution appearing as an overall factor in the expression for the cosmological constant. This factor comes from the components corresponding to i R = 1, 2, 3, 6 and is given by for any a 2 , b 2 , c 2 , d 2 ∈ {0, 1 2 }. By inspection, this identity also guarantees a vanishing contribution to the one-loop vacuum energy of the full non-SUSY SO(10) model above, from the untwisted sectors in which both 1 , 2 = 0 mod (2) (where 1 = 1 1 + 1 2 and similar for 2 ).

C.2 Massless Bose-Fermi degeneracy and the 1-loop q-expansion
The one-loop partition function after the applying the CDC is proportional to (C.7) The source of the exponential suppression of the one loop cosmological constant is then that, in the sectors where = 1 + 2 is odd (so that the contributions does not just vanish by supersymmetry), the q-expansion of the partition function is found to be missing the constant term due to the Bose-Fermi degeneracy among the massless states: (C.8)