Goldstino condensation at large N

We analyze the standard fermionic 4D Volkov–Akulov action with N global non-linear supersymmetries. We find that the stationary points of the system are described by an effective potential, written in terms of two composite real scalar fields, which approaches the exact quantum effective potential in the large N limit. We identify the existence of at least two stationary points, one representing the original supersymmetry breaking configuration and the other one corresponding to goldstino condensation, where supersymmetry is restored in the deep IR.

Fermionic condensation often plays an important role in understanding the vacuum structure of a quantum theory.One of the basic systems where fermionic condensation has been understood is the Nambu-Jona-Lasinio model [1,2], together with its 2D counter-part, the Gross-Neveu model [3].(For textbook discussions see e.g.[4,5]).An aspect that makes these models stand out is that, due to a large N number of fermion species, the quantum effective potential can be evaluated with arbitrary precision in a 1/N expansion, and the stationary points can be analyzed with confidence.
In [6] the existence of the goldstino tachyonic instability was demonstrated within 4D N=1 non-linear supersymmetry by means of the so-called functional renormalization group flow [32].
The adoption of such a method was necessary because a perturbative loop expansion would not be trustable.In general, however, one would like to have at hand a model where the goldstino condensation effect can be understood directly with the use of more conventional methods of quantum field theory.In particular, a large N demonstration of the existence of this phenomenon would be welcome.In the present work we take the first steps towards filling this gap.It is of course the 4D N=1 goldstino condensation that has the most phenomenological value, but it is theoretically important to analyze the generic existence of goldstino condensation from various perspectives 3 .
We will work with a single and well-established model that includes only fermions: the Volkov-Akulov (VA) model with N non-linearly realized supersymmetries [34][35][36].Since these supersymmetries are non-linearly realized, their number can be arbitrarily large, and the same holds for the number of the accompanying Goldstone fermions, the goldstini.This means that we are not restricted to the typical N ≤ 4 of linearly realized supersymmetry (or N ≤ 8 for supergravity).In addition, our aim here is neither to extract nor to study a phenomenological result, but we would like, instead, to answer one and only one question: can a system with N goldstini have a well-controlled stationary point described by goldstino condensation?Our main result is an affirmative answer to this question, thus giving complementary support to the findings of [6].Our results further suggest that such a configuration corresponds to the restoration of supersymmetry, and we also give some arguments in favor of its protection from higher-order corrections.Finally, we analyze a system that is more relevant for string flux compactifications, where, among N fermions, only one is a 4D N=1 goldstino, with all the others becoming pseudo-goldstini, once they acquire a mass.

The effective action and large N
We work with a system that has N non-linearly realized supersymmetries and focus explicitly on the goldstino sector, which is described by the Lagrangian where the goldstino vielbein A m a is defined as and f is the supersymmetry breaking order parameter whose mass dimension is Note that the actual supersymmetry breaking scale is N 1/4 √ f .The system has a global U(N) R-symmetry under which the spinors G I and G I transform in the fundamental and antifundamental representation, respectively.The goldstino vielbein A m a is, instead, a singlet under such U(N).This theory is defined with a cut-off Λ for which we typically assume For later convenience, and as it is typical in large N models, we have already extracted the N coefficient in front of the starting Lagrangian (2.1): the 't Hooft limit [37] thus corresponds to 3).Let us also observe that the number N of non-linear supersymmetries clearly matches the number of goldstini I, J = 1, . . ., N, whose transformations are Making use of the definition (2.2), the leading-order terms of the Lagrangian are ( The constant −Nf 2 /2 can be always removed so that there is no divergent term in the large N limit.From now on we will stop inserting the explicit N I=1 and summation over the same I, J indices will be implied unless otherwise specified. For the sake of completeness, note that we are using the conventions of [5]: therefore, η ab = diag(+1, −1, −1, −1), σ a = (1 2 , σ), and, because it will be useful later on, γ a = 0 σ a σ a 0 , with To pursue our aim, namely to get access to stationary points that possibly describe goldstino condensation, we write the theory (2.1) as Once we integrate out C a m , we get and we recover the model with N goldstini (2.1).Equivalently, the path integration over C a m yields a delta-function at each point of space-time that enforces the constraint.The advantage of (2.6) is that the action becomes Gaussian in the fermions and their path integration can be performed, leaving behind only a bosonic theory.In doing so, we will also be able to explicitly demonstrate the large N behavior of the model under consideration.
We start with the path integral where we can split the Lagrangian as and (2.10) The N 0 stands for the overall normalization of the path integral.To evaluate the fermionic contribution to the path integral we pair the N Weyl goldstini into N/2 Dirac spinors as follows This seems as if we were assuming that N is even.However, we would like to stress that this is just a formality that allows us to easily evaluate the fermion path integral: for odd N the end result would be the same.Thus, the fermionic contribution to (2.8) formally reduces to up to the overall N 0 factor that we will shortly discuss and specify.
If we now bring the determinants into the exponential and then into the effective action we get the full bosonic theory where N crucially appears as a global factor, being at the exponent of the determinant in (2.12).
The power of the large N construction is now manifest: by taking N parametrically large we can make the classical effects dominant with arbitrary precision over the quantum effects because higher-loop contributions are always suppressed by N factors compared to the tree-level term.
This happens because of the overall N in front of the Lagrangian and means that in the large N limit the quantum effective potential and any stationary points are controlled by the scalar potential of (2.13).
Let us notice that the fermion functional determinant contribution is clearly missing a dimensionful normalization inside the logarithm.This is related to the choice of N 0 in the path integral normalization, which we should specify.For instance, one can either introduce a scale to match the dimensions or insert the inverse propagator i∂ /.Here we will choose N 0 to be Other choices of N 0 reflect the different ways that one can use to express the determinant in perturbation theory and should ultimately not matter.

The effective potential and stationary points
Let us now search for stationary points of the bosonic theory (2.13).From the form of its effective action we directly see that this theory has tensor fields that signal the presence of massive higher-spin excitations 4 .However, since we would like to search for translation-invariant and Lorentz-invariant stationary points, we are interested only in the trace parts of these tensors and, consequently, in the resulting scalar potential.We then describe the VEVs of C a m and e m a as follows: (2.15) In order to get the full scalar potential we need to reduce the functional determinant that includes C a m to a convenient form by treating C a m , and so h, as a constant background field.The reader could be concerned that by ignoring the tensor modes' contribution to the scalar potential we may be missing some non-trivial constraint on the trace parts of the tensors.However, this does not happen for the following reason.One can think of splitting both C a m and e m a into traceful and traceless parts as , where X a a ≡ 0 and Y a a ≡ 0.Then, it is easy to see that there can never exist a linear term containing either X a m or Y m a in the scalar potential simply because there is nothing to contract them with: all terms with X a m and Y m a in the potential are directly quadratic in these fields and, therefore, they can be set to vanish consistently when we are searching for a background solution.There can be kinetic mixing of all sorts, of course, but here we are discussing neither the dynamics of the system nor its spectrum.
To proceed, we normalize the fermionic determinant in (2.12) with i∂ /, and treat the bosons as a background.Focusing directly on the relevant contributions from (2.15) we can see that (2.12) becomes (2.16) We then evaluate the trace of the operator as the sum of its eigenvalues, that is where V T is the four-dimensional volume, which we write as d 4 x.The expression for the fermionic path integral can thus be brought to the form which, once it is evaluated in Euclidean momentum space with a cut-off, gives V F constitutes the full correction to the bosonic effective potential from the fermion Gaussian integral.For completeness, in the Appendix we will also discuss such a calculation from the perspective of dimensional regularization.
We conclude that the stationary points of the theory (2.1) are described by the critical points of the scalar potential

.20)
The scalars h and φ on which V eff depends should not be directly thought of as standard scalar fields, because they are physically parts of the tensor fields C a m and e m a : h is related to the trace of C a m and φ to the trace of e m a .In spite of this, the critical values of h and φ do correspond to stationary points of the system and their potential energies at the stationary points can be legitimately identified as actual energy densities [4,5].The equations for the critical points are and they can be combined into a single equation for φ, which in turn directly gives the value of h.The resultant equation for φ is This equation has at least two real solutions that can be easily found if we solve the system numerically, but can also be determined analytically if we solve the equations perturbatively.
Willing to adopt the latter approach, we assume that Under such requirement, we find and . (2.25) The first solution corresponds to the original VA point and the second one is related to goldstino condensation (GC).On these critical configurations h takes the values . The potential energies corresponding to the stationary points of (2.20) that we have just found can then be extracted: they are and respectively.We see that the configuration (φ GC , h GC ), which describes the large condensate, has lower energy than the VA point (φ VA , h VA ) (which seems to correspond to a small condensate; we will shortly come back to this point).We would like to highlight that, since N is assumed to be very large, the scalar potential (2.20) is arbitrarily close to the full quantum effective potential.
Therefore, its stationary points (φ cl , h cl ) (but only those) correspond to actual quantum states of the theory and the corresponding value of the scalar potential captures the energy density of the state with energy E = (V T )V eff (φ cl , h cl ).As an aside, let us also note that, since the overall value of a potential can always be shifted by a constant, these energy densities have to be considered as relative one to the other.We will come to this point in a while.
It is also worth observing that classically, from (2.7), we have This justifies the interpretation of a non-trivial φ background value as a signal of goldstino condensation.We also see that, when φ approaches unit, the VEV of the condensate approaches Nf 2 and possibly jeopardises the control over the possible higher-order terms.Indeed, since we are working with the pure VA term (2.1), we are certainly ignoring higher-order goldstino self-interactions.We will see shortly that, as long as the Lagrangian is written in terms of A m a and its derivatives, we maintain full control of the vacua at large N. Therefore, no matter what power of φ we have in the higher-order terms, the condensate is robust even when φ goes near unit.
For completeness, let us mention that, contrary to what we have done so far, due to the large number N we could also work under the assumption that √ f < Λ.After further assuming , in order to extract an analytic result, we find two almost degenerate solutions: The original VA point is not a part of the stationary configurations any more and only solutions corresponding to large condensation correspond to stationary points.Since the condensate (2.28) takes parametrically large values in this case, we will not pursue this limit further here.

The stationary points in the deep IR
In this subsection we would like to understand what flowing towards the IR means and how f changes while going to lower energies, assuming that we are already in a low energy regime where the VA model is weakly coupled, that is Since the VA model has a single coupling, f , it would be enough to evaluate the flow of any specific term or interaction, and the other interactions would change accordingly.This in principle would require to have a regularization scheme that respects non-linear supersymmetry (as, for example, in [6]).Within such a setup one could deduce the flow of f by considering the 4-Fermi derivative interaction.However, such a calculation would be quite involved for various reasons.
For instance, it would firstly require the identification of a proper regularization scheme, and then the evaluation of both the fermionic wave-function renormalization and the actual running of the 4-Fermi vertex.Nevertheless, because of (2.29), one expects that the loop contributions would be subdominant (i.e.O Λ √ f ) with respect to the classical running due to the mass dimensionality of f .Therefore, without going into the loop calculation we can focus on the vacuum energy, V vac = Nf 2 2 , as a tool to infer the classical flow of f towards the IR.Even though this is a crude analysis, we will not only be able to capture the dominant running, but we will also see that the assumption (2.29) is enforced by the flow itself.
Following [5], let us consider the VA model in the form (2.1) and the Euclidean path integral where Λ representing a momentum cut-off for the quantum field fluctuations (and without specifying the species index, for simplicity).We then distinguish the integration variables into two groups: the "high-momentum" degrees of freedom ( Ĝ, Ĝ) that have bΛ ≤ |k| < Λ, and the "lowmomentum" modes ( G, G) carrying a momentum |k| < bΛ (the parameter b being a fraction b < 1).We thus have and we further split the path integral normalization factor N 0 into the high-momenta and the low-momenta contributions as

.33)
Let us now integrate over the high-momentum modes and focus on the vacuum energy change.
The leading contribution to it comes from the Gaussian kinetic term of Ĝ and Ĝ.Once we integrate over such degrees of freedom, we find Higher-order contributions coming, for instance, from quartic terms of Ĝ or Ĝ are further suppressed by powers of f .Then, we rescale distances and momenta according to so that the variable k ′ is still integrated over the range |k ′ | < Λ, and the fermionic path integral has once more [D GD G] Λ as its measure.Therefore, we have where, in light of (2.34), Because of the rescaling (2.35), while considering (2.37), a decrease of the parameter b represents how much the system flows towards the IR.The expression (2.37) clearly shows that the vacuum energy tends to increase as b decreases.From V vac = Nf 2 2 we deduce that where f 0 is the starting value of f , before we integrate out any high-momentum modes (i.e. for b = 1), and the sub-leading contributions are of order Λ √ f 0 ≪ 1.This in turn implies that the coupling accompanying the higher-order interactions of the VA model becomes more and more irrelevant as one flows to the IR.
We thus conclude that in the deep IR regime b → 0 and gives In such a limit we can check the relative difference between the energy densities of the two critical points: we obtain (2.40) We see that the goldstino condensation point has parametrically lower energy than the VA configuration.Moreover, under the limit (2.39) the VA point recovers its classical energy density and all corrections to it vanish.We see that (2.40) is quite suggestive in favor of interpreting the GC point as a supersymmetry restoring field configuration.Such interpretation is further corroborated by the properties of the kinetic terms of the fermions.Indeed, in the deep IR regime (2.39), where a scalar VEV is properly defined, we observe that which means that one recovers the classical stationary point for the VA model, and Since we know from (2.10) that the kinetic terms of the fermions on a background defined by the stationary points are we conclude that at the VA point the fermions have canonical kinetic terms, whereas the kinetic terms of the fermions vanish on the GC point in the deep IR.This absence of appropriate Goldstone modes when h → −1 is consistent with the restoration of supersymmetry 5 .
As we mentioned earlier, in the Appendix we are going to present the same analysis by using dimensional regularization.We will see that in the deep IR regime that we just studied, defined by the limit b → 0, the results from the two different regularization methods nicely match.
Our findings also connect with the ERG analysis for the 4D N=1 system that used superfields [6].There, the system is driven to an asymptotic supersymmetric point where the derivatives of the superpotential vanish, and so does the vacuum energy.In particular, the asymptotic supersymmetric point of the 4D N=1 system satisfies , where G is the N=1 goldstino and √ f the N=1 supersymmetry breaking scale.We can interestingly observe that (2.28) for φ ∼ −1 corresponds to a similar limit.For completeness, let us also note that the growth of f in [6] is controlled in the IR by its mass dimension.
3 Robustness against higher-order terms In this section we would like to understand how much the stationary points (φ VA , h VA ) and (φ GC , h GC ) are influenced by the higher-order terms that our starting model is ignoring.Even though there are terms that we cannot account for and may change the solutions, especially if they describe R-symmetry breaking, we will provide a simple rule of thumb for the circumstances when higher-order interactions could be dangerous.More precisely: (A) When higher-order terms appear only through the goldstino vielbein A m a and its derivatives, the goldstino condensation is always robust for large N.
(B) When higher-order terms also include explicit (A −1 ) a m ∂ m G I terms, the goldstino condensation may be jeopardized.
We will prove (A) and we will give two different examples of (B).
The reader should keep in mind that, if the goldstino condensation does restore supersymmetry, then higher-order corrections to the treatment above do not threat its existence; only non-perturbative corrections could do that.This is possibly the reason why it is easy to readily control a large class of higher-order corrections of the form (A). Despite of the lack of a proof, we will also see that the corrections of the form described in (B) seem to remain innocuous most of the time.

Corrections from goldstino vielbeins and matter
Let us start by considering the case in which the higher-order terms are expressed only by the goldstino vielbein A m a and its derivatives.Schematically, they have the form for some scale M and some powers R and T .Since we are focused on the stationary points of a scalar potential, such derivative terms do not change the outcome.More importantly, when the higher-order terms take the form (3.1), even if they are not only derivative interactions, they are always parametrically sub-leading in the large N limit simply because they have no N factor in front of them.We can conclude that for large N the goldstino condensation is not spoilt by such higher-order corrections.
Let us now make the discussion a bit more precise by assuming that we have some massive scalars coupled to the system in a way that preserves the existing non-linear supersymmetry.
We assume that such scalars are in their VEVs so that we can restrict ourselves to consider the Gaussian piece of their action.These scalar fields could represent some degrees of freedom that have been removed from the spectrum to deduce the low energy goldstino theory.Their impact can serve as a proxy for the higher-order corrections.We consider n real scalars b i with We are interested in evaluating the contribution of the functional determinant of the scalars b i to the effective potential.To this end, we expand e m a once more as e m a = (1 + φ)δ m a , treating φ as a background field.We get If we now assume that the N 0 for each of the scalar fields b i corresponds to a free massive scalar, after a Wick rotation, we obtain and, explicitly, (3.5) We would like to see how this new contribution to the effective scalar potential changes the stationary points.Even though it is not necessary, we can assume that the scalars b i are heavy, that is M 2 > Λ 2 .Without actually performing any further calculation, but simply exploiting the large N limit, we observe that the contribution (3.5) to the total scalar potential is parametrically subdominant with respect to (2.20) as long as Indeed, the analysis here falls under the general arguments for the robustness of the stationary points under higher-order deformations of the form (3.1).In particular, as far as (2.21) is concerned, the left-hand-side equation does not change, whereas the right-hand-side equation becomes From here it is evident that the deviation of (3.7) from (2.21) is arbitrarily small at large N.
Let us notice that we can also extend the above conclusion to a more general matter-coupled VA system.Consider a Lagrangian that has a matter part (made by scalars, vectors, spinors) of the form [35,45] Ultimately, its induced contribution to the quantum effective potential boils down to some ∆V eff (φ) and, as a consequence, the deviation from the original system is controlled by (3.7), therefore being arbitrarily small at large N.
We can conclude that goldstino condensation is quite a robust prediction of the large N non-linear supersymmetric theory, assuming that matter is coupled to the starting VA system via (3.8) (which inevitably preserves also the R-symmetry).

Explicit goldstini under derivatives
We now discuss terms where the goldstini explicitly appear under derivatives, thus breaking the assumption (3.1).These terms can possibly jeopardize goldstino condensation even at large N.For example, one can consider a term like for some complex dimensionful coupling g.Such a term potentially has a non-trivial impact: not only it is not of the form (3.8), but it also contributes to the large N functional determinant because it contains the N goldstini.However, since we are interested in scalar backgrounds, (3.10) takes the form and it is then clear that it will never contribute to the quantum effective potential of φ.One can in fact perform an integration by parts, treating φ as a constant background field (because we are interested in the properties of the effective potential), to obtain Let us notice that the same manipulation can be done by going to momentum space and assigning zero momentum to φ, as the standard procedure to evaluate the contributions to the quantum effective potential requires.We conclude that terms like (3.10), if present, do not jeopardize the new stationary point associated to goldstino condensation.
As a further example, we can consider the term for some complex dimensionful couplings g ′ IJ .This term will contribute to the quantum effective potential with a large N coefficient.However, it manifestly changes the number of degrees of freedom because it leads to G I ∂ 2 G J terms, which induce an additional massive fermion in the spectrum for each goldstino.For a consistent EFT such terms should be, in any case, independently highly suppressed.Even so, let us analyze the impact of (3.13), assuming that one uses it only as an interaction vertex.The easiest way to handle such a term is to package the goldstini once again into Dirac spinors Ψ A of the form (2.11), where A = 1, . . ., N/2.For our analysis we will also assume that the only non-zero contributions to (3.13) come from with g ′ ∈ R now.Then, considering only the background h and φ contributions from (3.13), the Gaussian fermionic sector is where the first term originates from (2.10).The functional determinant for a single Dirac spinor of (3.15) becomes We readily see why such a deformation changes the degrees of freedom and introduces new massive fermions.However, as far as our purpose is concerned, we simply need to treat the new contribution to the effective potential from the new massive fermionic functional determinant, having in mind that det [i(1 + h)∂ /] is already included in the effective potential and corresponds to the original goldstini.One way to do this calculation is to recast the overall functional determinant in the form The first term is similar to that which we have already calculated, but with 1 + h replaced by −g ′ (1 + φ) 2 .The second factor corresponds, instead, to the contribution of a massive fermion with canonical kinetic term 6 .We conclude that the potential that we have to extremize is (up to constants) where Contrary to the previous case, we cannot use the large N limit any more to eliminate the new terms.However, the theory still has a valid large N limit and the stationary points of (3.18) correspond to stationary points of the full quantum effective potential (in such a limit).
We would like to investigate the existence of a goldstino condensate for the scalar potential (3.18).In order to be able to continue analytically we make the assumption that there is a hierarchy between the scales at play, namely and we furthermore assume that We will check that this condition holds on the solutions.We can already see that it is satisfied for the VA point, if the latter persists.Under such assumption, and up to constants, the potential that we need to extremize takes exactly the form (2.20) at leading order in the Λ 2 /m 2 -expansion.
In particular, the corrections are all of the form Λ 4 ×(Λ 2 /m 2 +Λ 4 /m 4 +. . .).Therefore, both the VA configuration and the GC solution remain intact.Finally, one can check that both solutions satisfy (3.21), as long as This further implies that M ≫ √ f .
Let us observe that we have worked with a parametric separation between the various scales (M , √ f and Λ) that enter the problem so that we can easily deduce analytic results.Clearly, the solutions still exist for weaker assumptions but they have to be found numerically: we provide few numerical solutions for more conservative values of the coefficients in Table 1.We do not know under which conditions the solutions will seize to exist and if they seize to exist at all.
When M becomes smaller than √ f , the extremization problem can not be approached easily by Table 1: Few instances of numerical solutions for stationary points of (3.18) with Λ = 1, but without making any approximation on the effective potential.The numerical solutions approach the analytic ones as we go closer to the parametric limits that allow our approximations.The vacuum energy at the VA point is always in very good agreement with Nf 2 /2 and therefore we do not write it explicitly.Note that, because we are interested in the orders of magnitude and in the possible existence of a solution, we have rounded-up the presented numerical results.
the adoption of analytical methods, and also the numerical analysis seems to require stronger machines or more refined techniques.
We conclude that higher-order terms with explicit derivatives of the goldstini may seem harmful at first sight, but it is not obvious that they actually have an impact on the system after all.As the reader has appreciated, we have analyzed few such terms and we have seen that the properties of the stationary points do not considerably change.Nonetheless, we do not have a general argument to state that the higher-order corrections lying under the circumstance (B) cannot threat the goldstino condensate.As an aside final remark, let us note that other higher-order terms of a similar form can in principle be reduced to the Gaussian terms that we have studied by using Lagrange multipliers.

A single goldstino and N-1 pseudo-goldstini
In this section we wish to take advantage of the large number of fermions in the system to deduce a result for a model that has only N=1 non-linear supersymmetry.To do this, we make all the fermions massive but one, which corresponds to the single goldstino that the theory has in the low energy regime.
We split the goldstini as where G 0 , which we will denote as G from now on, represents the goldstino for the N=1 nonlinear supersymmetry, and the G i s are 2n pseudo-goldstini for a reason that will be clarified in a while.For convenience we pair the 2n pseudo-goldstini into n Dirac spinors Ψ A following (2.11).
We then explicitly break the extended non-linear supersymmetries down to 1 by introducing a Dirac mass term Note that, if each Dirac spinor is split into two Majorana fermions, then one gets a Majorana mass term for each one of them.Since we get back the full N goldstini system when these masses vanish, we call the massive fermions Ψ A pseudo-goldstini.We have with As elsewhere throughout this paper, we are interested in stationary points that are translationinvariant and Lorentz-invariant: we will consider directly the trace parts of C a m and e m a as in (2.15).
We want once more to perform the Gaussian integral over the fermions and derive the contribution to the effective potential for h and φ.To deduce the relevant modifications to it we perform two formal steps that allow us to get the result directly from the formulas that we already have at our disposal.First, we redefine all the fermions as follows treating h, as always, as a constant (because we are only interested in the effective potential critical points).As a consequence, the Lagrangian takes the form (4.6) The last term appears from the fermionic measure in the path integral because of the redefinition of the fermions.It has to be so in light of the fact that, if no field redefinition is performed, such a contribution appears from the Gaussian integral over the massless fermions.After (4.5) the massless fermion G decouples and it can be eliminated without any effect.This implies that the fermion redefinition has to contribute to the Lagrangian through the path integral measure.
Before proceeding, as a reminder for the reader, let us note that N = 2n + 1.
As we were mentioning just above, the goldstino G can now be integrated over without any effect, except of an overall shift in the vacuum energy, which we ignore.After eliminating G right away, the remaining 2n fermions are combined into n massive Dirac spinors with canonical kinetic terms and mass M (1 + φ) 4 /(1 + h).Knowing from the previous section how to evaluate the functional integral of these massive fermions, once we integrate out all the fermions by performing the corresponding Gaussian integral, we explicitly find where We see that the large N limit still gives us a reliable approximation to the full quantum effective potential with arbitrary precision.To derive the stationary points we extremize (4.7) with respect to h and φ.To avoid clutter we also assume that in the large N limit we can have Then, combining the equations for h and φ we see that we can readily eliminate h because it is bound to satisfy The system of the equations of extremization of (4.7) therefore reduces to a single equation for φ.It is possible to search for stationary point solutions without making any assumption on a hierarchy among the various scales at work, but, in this case, one has to proceed numerically.
For completeness, we give some numerical results in Table 2.
However, we can easily proceed analytically by first invoking the typical hierarchy Λ ≪ √ f and consequently observing that, under such requirement, (4.9), for the goldstino condensation point, gives assuming that These equations are in complete agreement with (2.25) and we will also verify them on the solution.Always within such limits the equation for φ GC is From this relation we can recover the goldstino condensation solution (2.25), if we assume that which (for (2.25)) reduces to M 3 ≪ Λf and therefore also to M 2 ≪ f .Consistently with these bounds we can still have M > Λ or M < Λ: the original goldstino condensation solution is intact for arbitrarily light or for quite heavy pseudo-goldstini.The VA point also remains intact when we have √ f ≫ Λ ≫ M .It is of course good that the limit M → 0 can be taken smoothly: in this limit the N-1 pseudo-goldstini become goldstini and we recover the results for the original N goldstini model of Section 2.
This result is relevant for string flux compactifications that include anti-D3/O3 systems as, for example, KKLT does (where supersymmetry is non-linear [46][47][48][49][50][51][52][53][54][55][56][57]).For the sake of the discussion, let us extrapolate our large N results to the case where N = 4.In [50]   making any approximation at the level of the effective potential.The numerical solutions approach the analytic solutions as we go closer to the parametric limits that allow our approximations.Since the vacuum energy at the VA point is always in very good agreement with Nf 2 /2, we do not write it explicitly.
Here again, because we are interested only in the orders of magnitude and in the existence of a solution, we have rounded-up the presented numerical results.
the extra fermions living on the anti-brane world-volume are discussed, in particular for the three massive fermions belonging to the anti-brane, which are pseudo-goldstini.The mass of these fermions is determined by the n (2,1) ISD flux.As we have seen here a small M or, equivalently, a small n (2,1) ISD flux may lead to goldstino condensation and further support the existence of such effects on an anti-brane.Conversely, if n (2,1) is large, then M becomes large as well, and the system has further issues due to large tadpoles [58,59].

Discussion
In this work we have investigated the existence of new stationary points in the standard 4D Volkov-Akulov fermionic system in the presence of N non-linear supersymmetries.An intuitive way to think of such an investigation and of our findings is the following.From the standard Lagrangian det[A m a ] describing the VA system (see (2.1)), one can derive the classical equations of motion and suspect that these equations have two types of vacuum solutions: Clearly, the vacuum solution where the goldstino vanishes corresponds to the original VA point that describes supersymmetry breaking: there, in fact, det[A m a ] = 1.The solution where the goldstino vielbein determinant vanishes, instead, corresponds to a condensation of the goldstini (see (2.2) for the form of A m a ), and implies that supersymmetry is restored, because the vacuum energy is now vanishing.The actual computation is more involved than simply solving det[A m a ] = 0 for the goldstini and proceeds with path integral methods [4,5], which allow to properly treat fermionic condensates.However, the naive intuitive expectation turns out to be correct and a solution of the form det[A m a ] = 0 does actually exist, as the path integral method that we followed here verifies.
As we already mentioned in the Introduction, even though in the present paper we are working directly in the component form and we are exploiting large N methods, our results lend further support to the goldstino condensation analysis of [6] that was performed with the ERG technique for superfields.These two approaches can be considered complementary and it is gratifying to see that they agree.It is also important to bring to the reader's attention the fact that bosonic systems with bosonic Goldstone modes can have a similar behavior where the classically broken symmetry is restored by quantum effects ( [60][61][62][63]; see [4,5] for textbook analysis).The fact that something similar happens for fermionic systems should not come as a big surprise, then.Moreover, this does not mean that supersymmetry cannot be broken, but it signals that the breaking of supersymmetry is more intricate than what one naively expects and has to be studied with care.
It is also worth noting that our results give further evidence that the anti-D3-brane/O3-plane system is inherently unstable on a flux-less Minkowski background: such a system corresponds, in fact, to N=4 [49].It is true that N=4 is not large, but the large N results may still persist.
In this respect, evaluating the leading 1/N corrections to the potential will be illuminating.
We have also discussed a system where all but one goldstini get masses.Such a setup corresponds to placing the anti-D3-brane/O3-plane system on a flux background [50].In this setup we have seen that the goldstino condensation persists.Such a model can also be studied with the use of constrained superfields satisfying X 2 = 0 = XY i [55], and exploiting the ERG technique to analyze the existence of condensates.The resulting (supersymmetric) backgrounds corresponding to the condensation may be ultimately related to some kind of brane-flux annihilation [64,65], but we cannot know if this is indeed the case yet.This is one of the important questions that we leave for future studies.In addition, the impact of including gravitation in our analysis is not necessarily trivial.
Another path that deserves to be investigated is how the goldstino condensation behaves in different dimensions.It is worth performing a similar analysis for example in 2D or in 3D, especially taking into account that spin-2 fields and gauge fields behave differently compared to the 4D case.An analysis of the condensate directly in 10D would also be illuminating, and especially interesting for the BSB models [66].However, such a study seems more challenging compared to that of the lower dimensional systems.
giving further support to what we state in the bulk, that supersymmetry has to be restored at that point.As we already mentioned, let us notice also that the goldstini stop propagating in such a limit at the goldstino condensation point.
Even though, because of the freedom to shift energies, the true value of the energy density in a QFT is a relative matter, we see that the original VA stationary point reaches its original energy density in the limit M → 0. We can therefore define the energy density of the supersymmetry breaking point with respect to that limit as (A.9) and in the same limit we also find thus interpreting the GC point as a supersymmetry restoration point.We see that this analysis agrees exactly with the analysis that we did in Section 2 for the deep IR limit.It is gratifying to see that the results do not change depending on the regularization scheme.For this reason we work only with the cut-off regularization prescription in the bulk of the article.

Table 2 :
Few instances of numerical solutions for stationary points of (4.7) with Λ = 1, but without