Positive energy and non-SUSY flows in ISO(7) gauged supergravity

We consider maximal gauged supergravity in 4D with ISO(7) gauge group, which arises from a consistent truncation of massive IIA supergravity on a six-sphere. Within its G$_2$ invariant sector, the theory is known to possess a supersymmetric AdS extremum as well as two non-supersymmetric ones. In this context we provide a first order formulation of the theory by making use of the Hamilton-Jacobi (HJ) formalism. This allows us to derive a positive energy theorem for both non-supersymmetric extrema. Subsequently, we also find novel non-supersymmetric domain walls (DW) interpolating between the supersymmetric extremum and each of the other two. Finally, we discuss a perturbative HJ technique that may be used in order to solve for curved DW geometries.


Introduction
One of the most controversial aspects of string phenomenology is the possible existence of a huge landscape of vacua of the theory, covering a wide range of lower dimensional features. In the original argument of [1], the string landscape was estimated to contain O(10 500 ) different vacua, while a more recent F-theory analysis extends it even further to O(10 272000 ) [2]. However, these estimations are usually thought of as very rough, since they try to count all distinct sets of topological data of the compact manifold, together with all possible internal flux configurations, without though really taking the dynamics into account.
On the other hand, regardless of what the exact order of magnitude may be, the concept itself of such a wide landscape seems to spoil predictivity of string theory and moreover, appears to be in tension with the idea of a theory of everything with no free parameters. By adopting a proper top-down approach, one might say that the only lower dimensional constructions that have all rights to be in the landscape are those that can be directly obtained by a string compactification within the right perturbative corner, i.e. where quantum loop, higher derivative and instanton corrections are all under control.
This way, one would naturally exclude any other kind of lower dimensional EFT that seems perfectly consistent from a low energy viewpoint and yet it has no known string theory origin. All such constructions should then be relegated to the so-called swampland. Nevertheless, since there might be constructions that do admit a stringy embedding which is simply very difficult to find or just not yet fully understood, it may be worth while looking for a set of purely bottom-up criteria that one can use to distinguish the landscape from the swampland. This was exactly the philosophy proposed by the authors of [3] when they formulated the weak gravity conjecture (WGC) as the first swampland criterion. This criterion states that any consistent EFT must retain gravity as the weakest force in the game, at any energy scale. As a consequence, within an EFT with gauge interactions controlled by a coupling g, the existence of microscopic particles (with mass smaller than the charge in Planck units) is needed in order to obstruct taking the ungauged limit of the theory g → 0.
In more recent years, a whole String Swampland Program [4] has been developing with the aim of defining an entire set of swampland criteria based on quantum gravity principles, that should help us select good stringy EFT's from the rest. Of course, to start with, these criteria are formulated as conjectures, which are supported by stringy insights, and only checked explicitly in some specific string theory setup's. This way, we came to conceive a web of different swampland criteria, which get more and more refined and better understood as time goes by, and interesting relations among all of them arise. For a complete review of the topic we address the reader to [5] and references therein.
In this paper we want to focus our attention on the so-called AdS Swampland Conjecture proposed by [6], according to which all non-supersymmetric AdS vacua must be ultimately unstable in string theory. The intuition behind this relies on a stronger version of the WGC applied to membranes (see also [7] for a similar logic). According to this, if the vacuum is supported by flux (which is generically the case), it will be unstable against spontaneous nucleation of microscopic membranes (tension smaller than the charge) that will eventually discharge the flux and hence destroy the associated AdS vacuum. Note that such a process is intrinsically non-perturbative and it does not affect supersymmetric vacua, since the corresponding spectrum only has BPS membranes (tension equal to the charge).
Though the intuitive picture is very clear, there are still very few explicit examples of such stringy constructions within a controlled setup (see e.g. the one in [8], where spherical D3 brane nucleation is found to destroy the non-supersymmetric and yet tachyon-free orbifolds of AdS 5 × S 5 ). In this paper we study non-supersymmetric AdS 4 × S 6 vacua of massive IIA supergravity. Such solutions are captured by a 4D gauged supergravity theory with gauge group ISO(7) [9], and in particular its G 2 invariant subsector. Remarkably, its G 2 preserving non-supersymmetric vacuum has already passed all perturbative checks, including the computation of the full KK spectrum [10], as well as brane jet instabilities of different sorts [11]. One really needs a non-perurbative analysis to explicitly check validity of the AdS Swampland Conjecture in this case.
In this paper we try to take a complementary approach and check for the possible existence of a dynamical protection mechanism for the aforementioned non-supersymmetric vacuum, by trying to derive a positive energy theorem, analogous to the one protecting supersymmetric AdS [12]. We will study the problem within a 4D N = 1 supergravity model capturing both the supersymmetric extremum as well as the non-supersymmetric one. The technical tool we will use is the Hamilton-Jacobi (HJ) formalism, which will allow us to derive the aforementioned positive energy theorem. It is worth stressing that this has no direct implications concerning full non-perturbative stability of the above vacuum. The reason for this is that what we conventionally refer to as 4D reduced model is not a good EFT in a Wilsonian sense, it simply captures a self-contained sector of the theory in a mathematical sense. Therefore, the existence of a positive energy theorem will simply tell us that potential non-perturbative instabilities, if any, must be looked for elsewhere. Indeed, while the vacuum is found to be non-perturbatively stable by virtue of the positive energy theorem within this sector, a non-perturbative instability completely analogous to the one in [8] is found in the companion paper [13], where the role of D3 branes is taken there by D2 branes.
The paper is organized as follows. In section 2 we introduce our setup, i.e. a particular 4D minimal supergravity model admitting SUSY as well as non-SUSY AdS extrema. In section 3 we discuss the central role of first order flows when studying bubble nucleation and positive energy theorems. In section 4 we focus on radial flows featuring a codimension one flat slicing of the metric and introduce their HJ formulation. In section 5 we review a general positive energy theorem and give a constructive proof thereof in our reduced 4D model. This will imply a dynamical protection from decay for the non-SUSY vacua within our model. In section 6 we briefly discuss radial flows involving a curved slicing of the 4D metric. Within this context, we present a viable algorithm that allows one to perturbatively integrate the HJ problem at any desired order. Finally, we collect in appendix A some details concerning the massive IIA supergravity origin of the 4D reduced model studied here.
In these gaugings the gauge group is embedded in the full E 7(7) global symmetry of the ungauged theory through its maximal SL(8, R) subgroup. The corresponding embedding tensor of the theory lies within the 36 ⊕ 36 of SL(8, R): The standard ISO(7) theory as arising from the reduction of massless IIA supergravity on a six-sphere [17], corresponds to the choice However, one may turn on an extra parameter within this gauging, which essentially also takes a possible symplectic deformation [18] into account. This gives rise to which of course solves the quadratic constraint (2.2) as well as the standard ISO(7) theory did. In type IIA language, the new p parameter controlling the symplectic deformation turns out to correspond to the Romans' mass, and hence the theory can be obtained as a consistent truncation of massive IIA supergravity on a six-sphere [19]. Some details concerning the truncation to its G 2 invariant sector are collected in appendix. As we will see later on, contrary to its massless counterpart, this theory contains non-trivial AdS 4 vacua.
The 32 real supercharges of the theory are organized into the fundamental representation of the R-symmetry group SU(8) R . When restricting to the G 2 invariant sector, these decompose as Hence, the truncated theory enjoys N = 1 supersymmetry. The 70 physical scalars of maximal supergravity decompose as where the G 2 singlet scalars are organized into one complex field named S. The corresponding embedding tensor deformations contain four real G 2 singlets in total, two coming from the A 1 piece of the T-tensor and the remaining two coming from the A 2 piece. However, out of these four parameters, only two remain independent when specifying to the ISO (7) gauging. In terms of minimal supergravity language, these two parameters generate the following superpotential deformation The Kähler potential determining the scalar kinetic coupling reads The corresponding action is given by [19] where κ 4 denotes the 4D gravitational coupling, K SS ≡ ∂ S ∂SK, and the scalar potential V is determined by the superpotential through 1 where D denotes the Kähler covariant derivative, e.g. D S W ≡ ∂ S W + ∂ S K W. After fixing the specific scalar parametrization given by the explicit expression of the scalar potential reads 2 This potential admits three different critical points, one of which is supersymmetric, while the remaining ones have sponteneously broken SUSY. With this choice of embedding tensor normalization, the three critical points are illustrated in table 1, while in figure 1, we show the level curves of the scalar potential around the three critical points. In [11] ID SUSY G res.  (7) gauged supergravity preserving at least G 2 as residual symmetry. The squared masses of the G 2 singlet modes are normalized to the absolute value of the cosmological constant. In these units the BF bound [20] lies the complete mass spectra for the entire set of scalar fields within maximal supergravity was computed and, while the SO(7) invariant critical point exhibits the presence of modes below the BF bound, the non-SUSY G 2 extremum remains completely tachyon free. In addition to this perturbative analysis, still in [11], a different stability analysis was performed, i.e. the one w.r.t. brane jets by making use of the probe approximation. As a result, the G 2 point was found to be also protected from this type of instability. Finally in [10], the full KK spectrum was checked by using ExFT techniques and even this calculation showed no signs of instabilities for this vacuum.
As a consequence, one might argue that this special example of non-SUSY AdS vacuum of massive IIA supergravity poses a challenge to the common lore brought up by the Swampland Program, which generically expects instabilities to arise, both at a perturbative [21,22] and non-perturbative level [6], whenever SUSY is broken. In the following sections, we aim at studying possible dynamical protection mechanisms preventing nonperturbative instabilities, at least within controlled EFT setups such as the 4D N = 1 description at hand.
We will indeed be able to derive a positive energy theorem for all vacua in this truncated theory that will exclude the occurence of non-perturbative decays within this universal sector. As already stressed in the introduction though, this will of course not imply the full stability of the corresponding vacua, but will only tell us that possible non-perturbative decay channels should be searched for within a more general setup. In analogy with [23], Figure 1: The level curves of the scalar potential for the G 2 invariant sector of ISO (7) gauged supergravity. The horizontal axis represents the ϕ direction, while the vertical axis represents χ. Both the SUSY & non-SUSY G 2 extrema are local minima, while the SO (7) critical point is actually a saddle.
the technical tool used here to our scope will be the HJ formalism, which will allow us to translate the problem of (non-)existence of certain radial flows into the (non-)existence of HJ generating functions with certain local and global properties.
solutions started to be considered in situations where the singularity can be resolved by a brane in string theory. In these cases, the corresponding non-perturbative decay process described by the associated geometry is the spontaneous nucleation of a charged spherical brane. Therefore, in a more general sense, non-perturbative instabilities of a given AdS vacuum may be seen as solutions to the reduced dynamics of gravity coupled to scalars, which describe a foliation of the lower dimensional metric (4D, in our case) with an FRW slice in one dimension lower (3D, in our case) where A, B and the scalar fields are suitable functions of the radial coordinate r.
In order for this to actually represent a viable decay channel, one must of course make sure that the following requirements are met: • the effective FRW cosmology on the bubble wall is of an expanding type, • the asymptotic geometry is our (putative) AdS vacuum, First of all, one can easily check that the dynamics of the system in the absence of coupling with extra matter, implies the slicing to be maximally symmetric, hence AdS 3 , Mkw 3 or dS 3 , depending on the value of the curvature on the 3D slices. As a consequence, expanding bubbles will necessarily correspond to dS slicings of our AdS 4 . Then the specific type of instability will crucially depend on the global properties of the solution, and especially on its asymptotic behavior at the other boundary of the r coordinate. If it asymptotes to a different AdS geometry, what we have is a CDL bubble [24] describing decay to a true vacuum through bubble neacleation. If it approaches a singular metric, we need to look at its UV description in terms of string/M-theory in order to understand the physics behind. It could be a bubble of nothing (BoN) [25] if the singularity gets resolved in higher dimensions into a smoothly shrinking Euclidean geometry, or conversely it could correspond to a spherical brane nucleation process if there truns out to be a stringy source placed at the singularity, or again it could instead correspond to a decompactification limit.
Secondly, the last condition given in terms of the variation of the Euclidean action may be viewed as a zero energy condition for the bubble geometry. In other words, the Euclidean action is extremized at a finite value only if its asymptotically AdS geometry has zero energy w.r.t. empty AdS itself. For a CDL bubble in the thin wall limit, i.e. where the scalar fields are approximated by step functions and hence carry no kinetic energy along the radial flow, this zero energy condition is understood as a bound on the tension of the (thin) bubble wall, which is usually referred to as the CDL bound. This condition was directly interpreted in [26] as a consequence of energy conservation in a bubble nucleation process. Later, in the classification of [27], bubble walls respecting the CDL bound are called ultra-extremal walls and are indeed found to lead to instabilities [28].
In this context, flat domain wall (DW) geometries represent a limit of the aforementioned bubble walls when the Euclidean action can only be extremized at an infinite value of the bubble radius, and this happens in turn when the CDL bound is exactly saturated. However, these do not correspond to real instabilities but rather as walls of marginal stability separating two phases that would require infinite time to decay into each other.
In some sense, one might already expect that the existence of interpolating solutions with flat slicing (Mkw 3 ) could in itself exclude the existence of actual bubbles with the same initial and final points.
More precisely, it's within the same context of flows with flat slices that one can find the objects proving a positive energy theorem for a given AdS vacuum. By this we mean that the existence of some special flat radial flows automatically implies that any bulk perturbations of the AdS asymptotic geometry will necessarily have positive energy, and hence no finite Euclidean action and therefore no significant contribution to the path integral. In the next sections, we will first formulate HJ for flat flows and show its relation to fake supergravity, and subsequently use it to give a constructive proof of the positive energy theorem for the G 2 invariant vacua of ISO (7) gauged supergravity, as well as to find interpolating flat DW's.

Flat first order flows: separable HJ treatment
Flat radial flows within the G 2 invariant sector of ISO(7) gauged supergravity may be found by exploiting the first order formulation through the HJ formalism. These flows are solutions of the form where A & B, as well as the scalar fields ϕ & χ, are functions of the radial coordinate r.
It is worth mentioning that B is actually pure gauge. However, as often in these cases, it may be useful to make use of this gauge fixing freedom in order to better integrate the first order flow.
When plugging the Ansatz (4.1) into the action (2.9), one obtains the following 1D reduced Lagrangian 3 where denotes differentiation w.r.t. r, and V is the scalar potential introduced in (2.12). It is worth stressing that evaluating the supergravity action on a specific Ansatz as we do here does not lead to a consistent set of field equations in general. Nevertheless, in this particular case, it has been checked explicitly that the 1D Lagrangian (4.2) yields exactly the same set of field equations as those obtained from the 4D equations of motion evaluated on (4.1). Now we introduce the conjugate momenta for our dynamical fields According to the HJ formalism, the second order dynamics associated with the Lagrangian (4.2) turns out be equivalent to the following first order flow defined by a suitable gener- when supplemented with the HJ constraint H eff | π I →F I ! = 0, which fixes the form of the generating function F (A, ϕ, χ). In this case, the HJ constraint takes the following explicit form Since the effective potential has a factorized dependence of the warp factor A, the corresponding generating functional may be found by casting a separable Ansatz F (A, ϕ, χ) = e 3A f (ϕ, χ), where the function f of the scalar fields satisfies the following non-linear PDE where K ij ≡ diag 2 7 , 2 7 e −2ϕ is nothing but the inverse kinetic metric. This identity manifestly shows the interpretation of these first order flows as a realization of fake supergravity [29].
Indeed, one can show that is a global solution to the PDE (4.7), corresponding to the actual superpotential of the theory. Thanks to the existence of the superpotential, one may construct a globally bounding function as −3f 2 SUSY , which dynamically protects the SUSY vacuum from nonperturbative decay by virtue of the positive energy theorem for supersymmetry [12] (see also [30,31]).

Positive energy theorems for non-SUSY vacua
Now, because the PDE (4.7) characterizing the function f is non-linear, f SUSY may not be the only global solution. In [32] it was studied how solutions of the aforementioned PDE can imply dynamical protection from non-perturbative decay. The claim can be summarized into the following theorem [33].
Theorem: If a scalar potential V (φ i ) in a given EFT can be written as for a suitable globally defined function f of the scalars such that: then other points in the scalar manifold M scalar have higher energy than φ 0 itself, whence φ 0 is stable against non-perturbative decay within the given EFT.
As remarked earlier, due to the non-linearity of the defining PDE, multiple global f 's might exist. Moreover, since any critical point would represent a singular initial condition, multiple solutions are possible even at a local level. By following the analysis in [23], one may construct all of the local solutions of (4.7) by perturbatively expanding around each of the critical points of the scalar potential. This way, we may look for explicit solutions of the form 5 where the numerical coefficients f (k) (φ 0 ) of the above expansion are in fact rank k tensors. After fixing the 0th and 1st order coefficients to be the second order coefficients satisfy second degree algebraic equations admitting 2 2 = 4 different branches of solutions. After fixing this discrete choice at 2nd order, the whole tower of higher order coefficients will be completely and uniquely determined by just solving degree 1 algebraic equations.
As required by the above argument, the particular choice of local f having a global extremum at a given φ 0 is relevant for discussing positive energy theorems and nonperturbative stability of the related vacuum. Different choices of local branches turn out to define flat DW solutions interpolating between different vacua. The correct discrete choice for the SUSY extremum labeled by '2' proving the positive energy theorem turns out be 6 the above values exactly matching the entries of the Hessian matrix of f SUSY , evaluated at φ 2 . For the SO(7) and the G 2 invariant non-SUSY critical points, these are respectively given by

Interpolating flat DW solutions
As far as interpolating DW's are concerned, there turn out to exist solutions connecting the SO(7) vacuum with both G 2 preserving ones, while the only way of connecting these 6 We have fixed m = 2     need not be monotic. What has to satisfy monotonicity is the function f when projected along the flow, and indeed, in both situations, this turns out to be the case. where we explicitly pulled out the dS radius L in such a way that the exponential warp factors appearing in the above metric are dimensionless.
The whole Hamiltonian formulation goes through almost identically to the previous flat case studied in section 4, and the resulting Hamiltonian precisely coincides with the one in (4.4), except that the effective potential now has an extra term induced by the dS curvature, namely 7 It may be worth commenting that, even though it seems a very marginal complication, it is in fact rather a drastic one. The reason for this is that adding such a potential term spoils the factorized structure of the A dependence which allowed us to decouple A from the system, when studying the flat case. Thus, we will now have to face a genuine 3-field situation. 7 We still stick to the B = 0 gauge for the sake of simplicity here.
We propose here an appropriate generalization of the perturbative procedure introduced in [23] and applied here in section 4, which in principle allows one to determine F (A, φ i ) at arbitrarily high orders. Due to the non-separability of our effective potential, we have to cast a more general Ansatz for F where the coefficients of our perturbative expansion are now promoted to functions of A.
The HJ constraint (4.6) specified for the Ansatz (6.3) will now impose a set of ODE's (each for any perturbative order) for the functions F (k) (A) rather than algebraic conditions as it was in the separable case. Nevertheless, the condition F (1) (A) = 0 still makes sure that the whole perturbative tower remain decoupled. At 0th order one finds where the dot denotes differentiation w.r.t. A, and we have set V (φ 0 ) ≡ − 3 2 . The above ODE is then easily solved by 8 Subsequently, in analogy with the flat case, the only crucial step in the whole algorithm is determining F (2) , as it satisfies a non-linear ODE. Generically, this will admit a continuous family of solutions depending on a choice of the integration constants c ij (each for every independent component of F (2) ). Once that one is fixed, higher perturbative orders will give rise to linear ODE's fixing each F (k) (A) in terms of the previously determined ones.
For our specific case, , (6.6) whose entries are fixed by the following system of 1st order ODE's where V ij (ϕ 0 , χ 0 ) and ϕ 0 are known data concerning the AdS vacuum, and the function F (0) (A) was determined in (6.5) by solving the HJ constraint at 0th order. After choosing a particular solution to the above 1st order system, all higher order coefficients F (k) (A) for k > 2 will be fixed by solving linear ODE's. This concrete procedure offers a viable numerical method to integrate the HJ PDE defining F even in the case of curved flows, i.e. where the separable Ansatz no longer works.

Concluding Remarks
In this paper we discussed positive energy theorems for non-SUSY AdS extrema within the G 2 invariant sector of maximal supergravity with ISO(7) gauge group. Given the massive IIA origin of the theory, this analysis is directly relevant to the issue of nonperturbative stability of non-SUSY AdS. Our results show the existence of a dynamical protection mechanism against non-perturbative decay of the non-SUSY vacua, enforced by the validity of a positive energy theorem. We were able to prove this by using HJ techniques.
However, it is worth remarking that this result by no means guarantees full nonperturbative stability of the non-SUSY vacua at hand. If nothing else, this is because the N = 1 model we adopted is not a good effective description in a physical sense, since it does not capture the lightest modes. Of course it is still a valid description in a mathematical sense, meaning that it describes a consitent set of 10D deformations of the metric and the rest of the IIA fields. In this context then, a fair conclusion to our analysis is the statement that possible non-perturbative instabilities of the non-SUSY G 2 vacuum have to be searched for within a different setup, becuase the positive energy theorem proven in this paper prevents the existence of any expanding bubble solutions within this truncation.
Thinking along the lines of [8], we actually know that, if an instability exists, it will be of a brane nucleation type, due to the presence of a Freund-Rubin flux F (6) that fills internal space completely. The source for this flux singularity has to be a D2 brane. The corresponding metric and dilaton would then be of the form where H D2 is a function of r. Now, for a fully localized D2 within the transverse R 7 , this function must be given by (1 + Q D2 r 5 ), which behaves as r −5 close to the source. It is straightforward to check that this behavior respects the truncation Ansatz in (A.1). In particular, a near horizon D2 solution is given by the following 4D scalar configuration As a consequence, any brane nucleation process involving a fully localized D2 should be captured by our minimal supergravity model. Because we have just proven that such instabilities cannot occur in our model, we have to conclude that an instability of this vacuum must necessarily involve non-spherically symmetric D2 charge distributions.
Indeed, in [13], the non-perturbative instability that the non-SUSY vacua are found to suffer from corresponds to D2 brane charge which is smeared over S 6 . Note that this requires scaling the internal components of the metric along S 6 and the 10D dilaton independently, and this goes manifestly beyond the truncation Ansatz in (A.1).

Acknowledgments
We would like to thank Pieter Bomans, Davide Cassani and Nicolò Petri for collaborating in a related project. We further thank Nicolò Petri for some valuable comments on a draft version of this manuscript. The work of GD is supported by the STARS at unipd grant named THEsPIAN.
A Consistent truncations of massive type IIA supergravity on S 6 In this appendix we include some details concerning the 10D origin of G 2 invariant ISO (7) gauged supergravity in four dimensions. Massive type IIA supergravity admits a consistent truncation on the six-sphere, for which the full non-linear KK Ansatz was determined in [19]. When restricting to the G 2 invariant sector thereof, this reduction precisely turns out to yield the minimal supergravity theory discussed in section 2. Its field content is simply given by the 4D metric g µν and the complex scalar field S = χ + ie −ϕ .
The consistency of the truncation automatically guarantees the equivalence between the 4D equations of motion obtained from the reduced Lagrangian (2.9) of ISO(7) gauged supergravity and the 10D equations of motion of massive IIA supergravity. In particular, by virtue of the uplift formulae given in (A.1), the three vacuum solution collected in table 1 may be directly interpreted as AdS 4 × S 6 solutions of massive IIA supergravity. More specifically, the SO(7) preserving vacuum is simply of a Freund-Rubin type, i.e. where the only non-vanishing flux (besides of course F (0) = m) is the (magnetic) 6-form flux filling internal space completely. In the G 2 invariant vacua, on the other hand, the S 6 retains its round metric, while non-trivial lower-form internal fluxes such as F (2) , F (4) and H (3) are turned on by exploiting the fundamental 2-and 3-forms defined by its NK structure.