Unstoppable brane-flux decay of D6¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\mathrm{D}6} $$\end{document} branes

We investigate pD6¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\mathrm{D}6} $$\end{document} branes inside a flux throat that carries K × M D6 charges with K the 3-form flux quantum and M the Romans mass. In such a setup brane-flux annihilation can proceed through the nucleation of KK5 branes. We find that within the calculable supergravity regime where gsp is large, the D6¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\mathrm{D}6} $$\end{document} branes annihilate immediately against the fluxes despite the existence of a metastable state at small p/M in the probe approximation. The crucial property that causes this naive conflict with effective field theory is a singularity in the 3-form flux, which we cut off at string scale. Our result explains the absence of regular solutions at finite temperature and suggests there should be a smooth time-dependent solution. We also discuss the qualitative differences between D6¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\mathrm{D}6} $$\end{document} branes and D3¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\mathrm{D}3} $$\end{document} branes, which makes it a priori not obvious to conclude the same instability for D3¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{\mathrm{D}3} $$\end{document} branes.


Introduction
A useful way to break supersymmetry in flux backgrounds is to insert branes that preserve different supercharges than the background fluxes and branes. Supersymmetry breaking by combining branes and antibranes is perturbatively unstable since the branes are mobile and can move towards each other and annihilate. The virtue of combining branes with fluxes that have opposite orientation is that there is no direct analog of brane/antibrane annihilation. Instead there can be brane-flux annihilation [1] which proceeds via first nucleating branes out JHEP03(2017)141 of the fluxes which then annihilate with the antibranes. The question then is whether this type of supersymmetry breaking can be metastable. If so, then it would be of great practical use in constructing string models for dS vacua [2], inflation [3], near-extremal black hole micro-states [4,5] and holographic duals to dynamical supersymmetry-breaking [1,[6][7][8][9][10][11] (and see [12] for related work).
In [1] it was indeed argued that metastable states do exist for p D3 branes in the Klebanov-Strassler throat [13] with three-form fluxes carrying M × K D3 charges, where K and M are respectively the H 3 and F 3 flux quanta. Brane-flux annihilation proceeds via the Myers effect [14]: the p D3 branes polarise into an NS5 brane wrapping a two-cycle inside the 3-cycle (i.e. A-cycle) at the tip of the throat. The NS5 brane carries D3-charge that is a function of its position on the A-cycle. When the NS5 brane pinches off at the opposite side of the 3-cycle it induces M − p D3 charges instead of p D3 charges. Effectively, the motion over the 3-cycle corresponds to brane-flux annihilation: out of the background fluxes M D3 branes have materialised and at the same time K dropped by one unit. The M D3 branes then annihilate with p D3 branes to give a supersymmetric vacuum with M − p D3-branes. Kachru, Pearson and Verlinde (KPV) argued [1] that this supersymmetry breaking is metastable in the limit of small p/M . This limit is the same limit in which the characteristic size of 3-sphere at the tip R 2 S 3 ∼ g S M 2 s is much bigger then the characteristic size of the D3 horizon R 4 D3 ∼ g s p 4 s . Hence one might expect the D3 to only modify the geometry locally. KPV argued the metastability by computing the potential energy for the NS5 brane and finding it has a local minimum at a radius r ∼ (p/M ) √ g s M s . In the past years concerns have been raised in the literature about the validity of this result. The first concern is that the NS5 potential was computed by using an action obtained from S-dualising the D5 action. The D5 action is valid at weak coupling and hence the obtained NS5 action can only be argued to be valid at strong coupling. An alternative analysis using the non-Abelian D3 action also suffers from the same problem since an S-dual version of the standard D3 action is required. A second concern is that the KPV computation works at the probe level and backreaction effects should therefore be subleading and tunably small. However many recent works have shown that the backreaction leads to a singularity which is difficult to interpret.
In this paper we study whether backreaction destroys the metastable state of D6 branes. For that we must first establish the metastability of D6-branes at the probe level which has not appeared in the literature before (see however [15]). We then indeed find that backreaction destabilises the D6 branes. We are unable to address the D3-brane stability conclusively but we comment on it.
In section 2 we recall the existing results on the backreaction of antibranes. Section 3 contains a description of the background flux solution in which we place the D6-branes. We then show, in section 4, that D6-branes give rise to metastable states at the probe level when p/M is small (M is the Romans mass, p is the antibrane charge). This meta-stable state is however washed away when backreaction corrections are added (assuming g s p is large enough to allow for a supergravity description) as shown in section 5. In section 6 we comment on the case of branes of lower dimensionality, and discuss a possible loophole that could allow the interesting case of D3 branes to escape this instability. We conclude in section 7. All explicit computations can be found in the appendices.

JHEP03(2017)141 2 Backreaction effects and instabilities
An explicit study of antibrane supersymmetry breaking is hard due to the lack of explicit solutions in the cases of interest. The best understood solutions at the moment are based on antibranes smeared over the A-cycle which has been most studied for D3 branes in the KS throat [16] (see also [17]). The only case for which a quantitative analysis is possible that does not rely on smearing are D6-branes since those solutions are described by ODE's [18]. 1 The qualitative feature that these studies have shown is the presence of a singularity in the H 3 flux density. This property can be shown to be unavoidable even for localised branes [21,22]. Two prominent interpretations of this singularity have appeared in the literature in the past years.
• Singularities are actually to be expected [23] since near the D3-brane we have an AdS 5 ×S 5 perturbed by three-form fluxes. Analogous, albeit supersymmetric, examples of this type are known to give rise to singularities that are resolved by letting the 3branes polarise to (p, q)-5 branes [24]. For smeared D3-branes there is a growing body of work that indicates that certain expected polarisation channels seem absent [25][26][27] and others lead to tachyonic modes [28][29][30][31]. Analogous results for fully localized branes are not as strong and a recent paper found that polarised D3-branes can be made regular by fixing boundary values of some fields [32]. It still remains to show that a full solution with these boundary values exists.
• The 3-form singularity of Dk solutions has a simple interpretation [33] which can be found in the sign of the divergent (but integrable) charge density F 6−k ∧ H 3 . This sign is opposite to the sign of the antibrane charge. The reason is that the flux is attracted gravitationally and electromagnetically towards the antibranes. If the antibranes are replaced instead with branes, the electromagnetic repulsion will exactly counterbalance the gravitational pull [34]. This entails an obvious instability since a clumping of fluxes carrying Dk charges will enhance the brane-flux annihilation. This can be seen at the level of the probe actions describing the Myers effect [27,35,36] or from the point of view of bubble nucleation [36]. If the flux clumps in an unbound way one necessarily crosses a critical value that causes immediate brane-flux decay. If this picture is correct there should exist a smooth time-dependent solution describing the clumping of flux, which reaches a critical value at which the antibranes decay and then the flux clumping should stop and the BPS state should be reached after the closed string radiation decouples. This picture is not inconsistent with the expectation that backreaction is small when p/M is small, but the essential property that makes this possible is that the local backreaction, which can be big, determines stability.
In this paper we demonstrate that the second viewpoint is correct for D6 branes at large g s p but still arbitrarily small p/M .

D6 brane solutions
In order to analyse brane-flux annihilation, we first require a background with fluxes that carry D6 charges without the presence of localized D6 branes. This is possible in massive IIA supergravity as can be seen from the Bianchi identity for F 2 : Here F 2 is the RR 2-form field strength, F 0 is Romans mass and H 3 is the NSNS 3-form field strength. Notice that the term F 0 H 3 appears on the same footing as the D6 source term N δ 6 . Clearly when F 0 = 0 this term acts as a smooth source for 6-brane charge.
A general study of backgrounds with these ingredients was carried out in [18]. It includes the "massive D6" brane with flat worldvolume and non-compact transverse space [37] and D6 branes with AdS 7 worldvolume and compact transverse space [18,38]. We start with brane-flux decay in the first example but before doing so we emphasize why the D6-brane is special.

Why the D6 is special
D6 branes stand out against antibranes of other dimensionality for their simplicity. This is most obvious in the description of the backreaction of D6 branes, which is described by ODE's [18,33,38] without having to smear the branes. A second feature-crucial for this paper-is that the annihilation of the H 3 -flux does not proceed via a polarisation into a higher-dimensional object. All Dk branes with k < 6 polarise into an NS5 branes [15]. D6 branes instead polarise into KK5 branes, which are smaller in dimension. However the KK5 assumes a circular isometry transverse to its worldvolume, one can therefore think of KK5s as 5-branes that are smeared over a circle. This circular isometry direction will live inside the D6 worldvolume such that the KK5 brane looks like a 6-brane and its backreaction is identical to that of the D6 brane since it carries D6 charge and tension. This is related to the fact that D6 solution is the only antibrane solution for which the singularity in H 3 is consistent with the singularity of the metric. In other words: the backreaction of the diverging H 3 flux does not destroy the local D6-metric, whereas it would for k < 6. For instance for k = 3 one finds that the H 2 3 scalar blows up near the horizon, but the horizon of a 3-brane should be smooth AdS 5 × S 5 . Hence the 3-form fluxes necessarily destroy the local AdS 5 throat, which is well known from the Polchinski-Strassler (PS) model [24]. In the PS model the singularity at the would-be horizon can be turned into a physical singularity corresponding to (p, q) 5-branes. Something similar can be expected for the D3 solution and indeed reference [32] has shown that a local 5-brane background with good singularities it is at least not inconsistent with having ISD fluxes at the UV and a non-zero (generalised) ADM mass.

D6 branes in flat space
We start with a flux background geometry with a flat D6 worldvolume. The metric and dilaton 2 is that of a D6 brane To have D6 charge dissolved in fluxes we include a non-trivial H 3 flux and Romans mass given by Here˜ 3 is the unwarped Hodge operator on the transverse space such that˜ 3 1 = r 2 dr ∧Ω 2 = vol 3 and Ω 2 is the volume form on the unit 2-sphere. We denote the line element of 7D Minkowski space by ds 2 7 . Further on, once we discuss the KK5 brane it will turn out to be useful to pick out a specific direction ψ that plays the role of the isometry direction of the KK5. We then write ds 2 7 = ds 2 6 + dψ 2 with ψ a circle direction. The fluxes obey a Hodge-duality relation: where the Hodge-star includes warp factors. This condition is identical to the famous Imaginary Self-Dual (ISD) condition for fractional D3 brane backgrounds after performing three T-dualities along the Minkowski directions. One can furthermore verify that (3.4) ensures a no-force condition for probe D6 branes. D6 branes on the other hand do feel a force. The Bianchi identity for F 2 (with N D6-brane sources located at the origin r = 0) implies˜ to which a spherically symmetric solution is where N is the number of D6-branes sitting at r = 0 and v is an integration constant. We will put N to zero for now and consider a background without any explicit D6 sources. This solution is then regular around r = 0. The solution has a peculiar singularity for

JHEP03(2017)141
In order to examine the local behaviour of the singularity we expand around it: and find This is precisely the local singularity structure of an O6-plane (see appendix A). A discussion of this singularity can be found in [39] where it was argued that strings are well behaved in the singular background. Our interpretation of the singularity as an O6-plane confirms that analysis. We will not go into more details on the global analysis of the solution since the computation we perform is local to r = 0. This solution is not supersymmetric. Supersymmetric solutions for D6 branes in massive IIA do not preserve the Poincaré symmetries of the D6 worldvolume as shown in [37,40]. The flux background is however extremal in the following sense • The expression of the metric and form fields is the one typical for extremal p-branes.
• There is a no-force condition for inserting D6 branes into the background whereas D6 branes are pulled towards the tip of the throat.
• It is T-dual to supersymmetry breaking using (3, 0)-fluxes in fractional D3 backgrounds which is rather well understood and does not influence any of the details of brane-flux decay in that context [1]. 3 Inserting D6 branes. To describe the same throat geometry but with a stack of D6's at the tip, we use the following general Ansatz preserving the worldvolume symmetries (possibly broken at finite temperature) and the transverse rotational symmetries [33] (in string frame): From combining the form field equations of motion one can deduce that It was shown in [18,33,41] that insisting that the metric has an D6 singularity at r = 0 and asymptotes to the ISD solution, we find that α goes to a non-zero constant near the D6-horizon 4 which we denote α 0 . This implies: since the combination e −7A+φ−f blows up at zero or non-zero temperature. At non-zero temperature this is solely due to e f going to zero near the horizon whereas the other fields remain finite. Whereas at zero temperature e f = 1 but e −7A+φ is combination that necessarily diverges near a 6-brane.
This singularity is not a coordinate artefact since it appears in scalar quantities like the H 3 -flux density: (3.12) One readily verifies that the singularity in the associated charge density M H is still integrable, since H 3 ∼ r −1 Vol 3 = rdrdΩ 2 near the source at r = 0. Another peculiar feature is the fact that the backreaction of the fluxes does not modify the local 6-brane singularity of the metric to some unphysical naked singularity. This is different for the Dk solutions with k < 6 as explained in the Introduction.
The sign of the singularity in λ is such that it corresponds to a diverging D6 charge density dissolved in the fluxes as illustrated in figure 1.
Given the interpretation of this singularity as a diverging charge density one should be worried about the stability of this background against brane-flux annihilation. Heuristically one expects a large charge density dissolved in fluxes to increase the probability of localised and mobile D6 branes materialising out of the flux cloud that subsequently annihilate against the D6 branes. This is the topic of section 4.

Warped AdS 7 with D6 branes
If one allows the D6 worldvolume to be AdS 7 instead of Mink 7 the space transversal to the brane can be made compact and conformal to S 3 [34]. The resulting solution was described qualitatively in [18,33]. A much deeper understanding of the solution as well as a better numerical control was reached in [38], where it was understood that the solution can even be supersymmetric, rectifying some statements in [33]. 5 These backgrounds are relevant as holographic descriptions of a large class of strongly coupled (1,0) CFT's in D = 6. The details of these holographic descriptions were uncovered in [44] (see also [45]).
What counts for this paper, all details aside, is that there are various AdS 7 solutions. They can be classified according to the way the six-brane charges are divided over the compact manifold. The D6 charges are either due to explicit D6 branes sitting at the poles of the S 3 or due to spherical D8 branes wrapping contractible S 2 's inside the S 3 as in figure 2. Those D8 branes can induce D6 charges when they have the right kind of worldvolume fluxes on them.
The presence of D8 branes carrying the D6 charges is a crucial difference with the Minkowski solutions treated above. In the latter case it can be shown that the D8 branes would never sit at a stable position and pinch off to become pure D6 branes [25]. Only in AdS 7 can they reach a stable position [42]. Crucially, when the charges are carried by spherical D8 branes they do not cause flux-clumping singularities [25].
The role of the D8 brane can be confusing so we spend a few more words on this. Fluxes can decay against the p D6-charges in two ways [15]: 1) either by having a D8 brane polarising and move all the way over the S 3 which decreases the Romans mass or 2) by having a KK5 brane moving inside the D6 worldvolume which decreases the H 3 -flux. In this paper we focus on the second channel. Concerning the first channel the results have been computed recently: for AdS 7 the D8 will be stuck at a stable position [42] but for Minkowski worldvolume polarisation into D8 branes is impossible [25].
The presence of the 3-form singularity strongly affects the KK5 channel since the density of 3-form fluxes drives the KK5 polarisation. Hence, we expect that KK5 polarisation does JHEP03(2017)141 not occur at all for an AdS worldvolume since the D8 polarisation regularises the 3-form singularity. This is stricktly speaking an assumption for our setup but has been verified explicitly in the T-dual picture [27]. In contrast with AdS, we argue in this paper that the KK5 branes will polarise without finding a meta-stable state for flat worldvolumes. This drastic qualitative difference is caused by having the D8 brane polarised in one case and not in the other [27]. Because if the D8 polarises it removes the 3-form singularity, which lies at the heart of the immediate brane-flux decay.
4 Brane-flux annihilation at probe level

Brane-flux decay and T-duality
To analyse the stability of D6-branes in the background in question we first study the probe approximation where the branes do not backreact on the geometry. We must identify branes that can carry the D6-charge away, similar to the spherical NS5-brane in KPV [1]. These turn out to be KK5-branes that source D6-charge on their worldvolume [15]. This may seem strange as the KK5 is a 5-brane whereas the D6 is a 6-brane, but the KK5 has one special transverse direction, the NUT direction, which assumes circular isometry. For all intents and purposes we can think of the KK5 as a brane that is smeared along this NUT direction and hence can carry D6 direction. In our setup the NUT direction of the KK5-branes is a circle in the D6 worldvolume parametrized by ψ. The relevant coupling in the KK5 action is where n is the 'worldvolume flux' and k is the Taub-NUT vector. The induced 6-brane charge depends on the position ψ and equals −(p−ψF 0 ). In order to not introduce monopole KK5 charges, we consider a pair of KK5 and KK5 that share a NUT direction ψ but are otherwise separated on that circle. We illustrate the D6 decay process in figure 3. The spacetime domain wall that mediates the decay of H-flux and D6-charge is a NS5-brane wrapping the S 1 but traveling radially in Mink 6 .

The probe potential
To compute the probe potential for KK5 branes in massive IIA we could use the actions derived in [46]. It is however technically easier to work in the T-dual frame, and use that the probe potentials are invariant under T-duality [15]. The T-dual setup involves a pair of probe NS5-branes, with opposite charge in the geometry Mink 6 × S 1 × M 3 where M 3 is the space transversal to the original 6-brane and is conformal to R 3 . The branes are localised on the circle and on M 3 . The pair of NS5-branes induce D5 charge between them as explained below. This calculation was already sketched out in [15] but we will go through the details here for completeness. We start by writing down the background T-dual to (3.2) along the ψ-direction, The function S takes takes the same form as above After T-duality the coordinate ψ remains periodic with ψ ∼ ψ + s . Before continuing on to the probe computation we review how a NS5-NS5 pair can induce a D5-brane charge. Let the NS5-brane be located at ψ = ψ 0 and the NS5-brane be located at ψ = −ψ 0 . The probe action relevant for these branes is The upper sign is used for the NS5-brane and the lower is for the antibrane. The field G ± is a worldvolume field where p is a constant and C 0 is the local gauge potential for F 1 . The D5 charge induced by the brane pair can be inferred from the sum of the above WZ couplings, where we have used C 0 = (M/ s )ψ and that C 6 is independent of ψ. The induced D5-charge is therefore controlled by the combination p − 2M ψ 0 . Since ∆ψ ≡ 2ψ 0 is the separation of the brane pair we see that when pulled together to the point ψ = 0 they induce the charge of p D5 whereas on the other pole ψ 0 = s /2 the charge is p − M . The D5-brane charge carried by the NS5-NS5 pair depends on the relative seperation of the branes. The computation of the probe potential is spelled out in appendices B and C. Including the interaction term in the full probe potential we find: The potential is displayed in figure 4 and shows that the hight of the barrier is controlled by M g s . One can verify that there is a metastable state for M g s large whenever Concerning the validity of our probe potential we suffer from the same pathologies as the original work of [1], i.e. the NS5 action is obtained by S-duality of the D5 action, which implies we can stricktly speaking only use it at strong coupling. The size of the S 1 is a modulus and has no effect in our computation so we choose it large enough for supergravity to hold. of vanishing p/M . For the special case p = 1 some initial steps towards a direct effective field theory approach were taken in [47] based on field theory renormalisation [48]. In what follows we will use an alternative description (see for instance [24]). This description is valid when p/M is small but g s p is nevertheless big. An effective field theory description of any kind will have to deal with the singularities typical to brane solutions, such as the singular warpfactors and diverging gauge fields sourced by the brane. Just as in standard electrodynamics one does not expect diverging field strengths to cause any harm as they correspond to an infinite self-energy that is effectively subtracted off in the probe action. For diverging warpfactor, the same story applies.
The diverging H 3 field on the other hand is not as clear since it is not sourced by the 6-brane directly. It is tempting to think of it as the diverging field strength corresponding to the NS5 (KK5)/D6 boundstate since NS5-branes do source singular H 3 flux. 6 However it has been found in [32,49] that this cannot be, at least not within the Ansatz presented in this paper, since then the divergence in H 3 is not of the right kind to correspond to a NS5 (KK5) source. As mentioned earlier the singularity really corresponds to a divergence in the flux clumping, which makes the D6 charge density dissolved in flux diverge near the D6 brane. Our approach to deal with this singularity is to peal off a small amount n p of D6 branes and consider them to be a probes in the backround of p D6 backreacting branes. Clearly this is only possible in the regime p 1, which we anyway require in our supergravity approximation. In that way the singularities in the warpfactor and the gauge fields, generated by the backreacting stack (see picture 5), cancel each other. This cancellation occurs because the leading contribution to the WZ term is the same as the leading contribution in the DBI term with a different sign. So we rely on that fact that branes, even when non-supersymmetric, locally behave like supersymmetric branes for which tension equals charge. Although this effectively eliminates the diverging contributions of the geometry and the gauge field that couples to the brane, it does not remove the singularity in H 3 (λ). The singularity in λ occurs at a place where the space-time metric is not to be JHEP03(2017)141 trusted due to high curvature. This implies one has to be careful when concluding that the infinite flux-clumping destabilizes the system due to direct brane-flux decay. We therefore cut off the solution a string-length away from the singularities and investigate how large the flux-clumping is.
In this section we find that the flux-clumping λ is proportional to g s p at the cut-off distance s away from the backreacting branes. This then implies that for large g s p, where we can trust supergravity at that scale, the system is destabilised independent of the ratio p/M being very small. The argument for using the string scale to cut-off can be found in [47,48]. Supergravity should not receive large corrections beyond that scale since for D6 branes in flat space the curvature scales are small enough already at that length scale.
Another systematic way to approach the effective field theory would be to consider high-order corrections in the blackfold approach [50,51], which we hope to adress in the future.

Flux-clumping at the cut-off
To evaluate the flux-clumping at a finite distance r = s from the antibranes, we need to evaluate λ there. Here we work directly with the D6 solution but later when we consider the T-dual background where the expression for λ remains the same. Near a stack of extremal antibranes we have Since λ = αe φ−7A , we can calculate λ once we know the value of α near the cut-of. The latter we obtain from the Smarr-like relations found in [21,32,49] (see also [52]), that relate the energy E (above the supersymmetric vacuum), to the values of the gauge fields at the sources: such that E = µ 6 pα(r = 0). To find α(0) we therefore need to know E. Following the reasoning of [6] we have (the background has e A = 1 at the tip) From the Smarr relation (5.3) we then find Concerning λ we have λ(l s ) = g 1/4 s π 2 pα(l s ). We take α(l s ) = α(0) and find λ ∼ g s p, consistent with the estimates in [47]. One can only consistently take α(l s ) = α(0) if α (0) is sufficiently small but this is easy to see. At the brane source we must have F 2 = Q Ω 2 and this provides an expression for α (0). Near the source we have: such that JHEP03(2017)141

The corrected probe result at large p
We now compute the probe potential for n D6-branes near the tip geometry with the backreacting of p D6-branes included. By cutting of the divergent flux-clumping by the value found above: λ ∼ g s p we find a finite and well-defined probe potential. To carry out this computation we need the solution T-dual to the D6 solution. A straightforward computation gives the following metric for the smeared D5-branes at r = 0: ds 2 = e 2A ds 2 6 + e −2A dψ 2 + e 2B (dr 2 + r 2 dΩ 2 2 ) , (5.8) Where as before, the equations of motion demand α = λe 7A−φ and the fields have the standard D-brane singularities describing smeared 5-branes at r → 0. We now follow the exact same procedure of appendix B to compute the probe potential by placing a single NS5/NS5-pair in the above background. This calculation can be found in appendix D. The result is a run-away potential: We note that the interaction term has a different dependence on the dilaton and warpfactor such that the warpfactor and dilaton do not simply factor to the front. The factor e 3A−φ in front of the interaction term reduces to g −1 s close to the brane, consistent with the uncorrected probe potential.
Note that this potential receives large corrections as ψ grows. The reason is that the decay of the branes causes at the same time a materialization of M − n D6-branes which quickly dominate over the p D6-branes such that the original background geometry cannot be trusted. Essentially, corrections are subleading until the value of ψ is such that the drop in energy is twice the tension of the n D6-branes. It is satisfactory that at the value for ψ where the induced 6-brane charge in the KK5/KK5 pair becomes zero, the drop in energy is two times the tension of the n D6-branes.
For a stable configuration we would expect that by dialling down n/M and p/M a metastable state would appear for the n mobile antibranes. This is however not the case which indicates that the whole stack of branes may be unstable. Of course when p keeps decreasing and eventually reaches small values for which our analysis is not valid. This is the limit of our approach and a different analysis must be carried out to determine the stability of small number of antibranes.
6 Dk-branes with k < 6 In this section we discuss how simple dimensional analysis can be employed to study flux clumping around antibranes of different dimensions.

JHEP03(2017)141
The flux-clumping for unpolarised Dk-branes can be inferred either via the supergravity techniques of [21,22] or the brane effective field theory techniques of [47]. Both give the same result, namely with r a local radial variable. Hence the flux clumping at the cut-off scale r = s always equals In the supergravity regime this is too large to ensure stability of the probe. Now the question arises how λ c changes when the branes polarise via the Myers effect as predicted by KPV [1]. The essential observation of this paper is that λ c does not change for D6-branes polarising into a KK5/KK5 pair, because the pair has the same dimensionality as the original branes. We expect similar conclusion for D5-branes and hence they should be unstable as well in the supergravity regime. However when k < 5 the Dk charge spreads out on a spherical NS5-brane of radius R such that one expects, from dimensional analysis, The largest possible value for R is the size of the full A-cycle, which we take proportional to R 2 max ∼ g s M 2 s like in the KS throat. For D4-branes we then have that this minimal value of λ c is still too large. However for D3-branes the situation improves. The radius at which the clumping is only of order 1 is found to be: This is far off the radius one finds from the probe computation, which is linear in p [1]: In fact at the probe radius the clumping can be estimated to be too high and should destabilise the brane. This already shows that the probe computation at best only gives qualitative information but certainly not quantitative results, very much in contrast with supersymmetric vacua [53]. Interestingly in reference [32] a √ p dependence of the polarisation radius was inferred on very different grounds: namely by using the aforementioned boundary condition for polarised branes that can regularise λ without the need to cut-off. This potential match seems interesting and we leave it for future research.

Conclusions
We now summarise our results and then give an outlook on the metastability of general antibranes.

Summary
We have shown that at probe level D6-branes experience a classical barrier against annihilation with fluxes of the opposite charge on the condition that the Romans mass quantum M is much bigger than the antibrane charge p, thereby extending the results of [15]. This condition is (qualitatively) T-dual to the condition for the metastability of D3-branes [1], where Romans mass gets replaced with RR 3-form flux. We have then shown that backreaction washes away the classical barrier, confirming the arguments of [35] and [36].
The reason this could be firmly shown rests on 2 facts: 1. D6-branes are special in the sense that they brane-flux decay via mutating into a KK5/KK5 pair carrying D6-charge. The pair has effectively the same dimensionality and backreaction as localised D6-branes.
2. The backreaction of localised D6-branes is described by ODE's. Near the source the backreacted geometries are T-dual to Dk-branes with k < 6 smeared over the k − 6 directions inside the A-cycle that is filled with RR-flux. For such branes there was already strong evidence that they are unstable [27].
The vanishing of the classical barrier upon backreaction is due to the attraction of the fluxes towards the D6-branes, something that is not taken into account in the probe approximation. This increase in flux diverges close to the antibranes at the classical level. Hence we followed the procedure of [47,48], which amounts to cutting off the singularity at string scale. The main observation is that the cut-off value for the flux density is still too high and causes immediate brane-flux decay. Effectively the increased flux density causes an increase in the probability that actual D6-branes materialise out of the flux cloud which subsequently annihilate with the D6-branes. Our result explains the lack of smooth D6 solutions with finite temperature [41].
The reason for this apparent failure of the probe approximation could be due to the fact that the probe calculation is outside of the regime of validity of the NS5 (or KK5) brane action.
We have also investigated whether the same picture could still be valid for Dk-branes with k < 6. Using dimensional analysis, similar to [47] we argue that the same physics should hold for D4 and D5 branes but that D3-branes cannot be shown to be unstable using the method described in this paper. We do however find that, if a metastable state exists for D3-branes, the radius of the spherical NS5-brane must scale different with the charge p from what is predicted by the probe calculation [1].

Outlook
It has been argued that the instability we find can be circumvented by considering small brane charges (e.g. p = 1) [47,54]. For small brane charges the cut-off value of the flux clumping is small and furthermore there is no polarisation potential that can be computed within the supergravity regime. Given that flux clumping is not severe for small p, an instability cannot be argued solely based on a singular flux density. Instead brane effective field theory should be employed to argue for stability [47]. Stability would most likely JHEP03(2017)141 imply a finite temperature solution with small p in supergravity. However for D6-branes the absence of a smooth finite temperature solution was demonstrated in [41]. It could of course be that the metastable state at p = 1 has a very small gap such that the minimal temperature needed for a classical horizon already destabilises the vacuum [54]. This potential loophole might become visible by computing finite temperature corrections to probe actions, something that can be done in the blackfold approach [50,51].
Despite the possible room for the existence of metastable D3-branes, we want to be open-minded and contemplate about the possibility that also D3-branes are unstable against the flux-clumping effect [35,36]. As long as numerical evidence showing that the boundary conditions for regular solutions described in [32] and [49] can be taken, does not exist, the option of an instability has to be taken seriously. There exists a heuristic interpretation of the "Smarr relation" for flux throats that supports this interpretation.
The generalised ADM mass M ADM of a flux throat with sources, measures the energy above the supersymmetric vacuum. It was shown that M ADM obeys a Smarr-like relation [32,49,55] where Q M (Q D ) stands for monopole (dipole) charge and Φ M (Φ D ) denotes the gauge field sources by the monopole (dipole) charge evaluated at the horizon in a specific gauge. The exact expression, on the nose, corresponds to the on-shell Wess-Zumino action for the would-be polarised brane (see [21] for a similar derivation using on-shell brane actions). This almost fits the standard expectation that the energy of the flux throat equals twice the antibrane tension since where the DBI-term and WZ-term exactly equal each other in the probe approximation. In the backreacted case it seems all energy is therefore carried by the WZ term alone. This comes about since the warpfactor redshifts the on-shell DBI to zero, and alternatively this can be seen as a flux clumping effect, which makes the WZ dominate over the DBI. If this interpretation is correct then a metastable state can never exist since metastability truly requires a competition between the DBI and WZ term. We do want to point out that this is not a proof since our interpretation of the Smarr relation as on-shell probe potential is heuristic. Finally, we have not touched upon the 3-form singularities in those (supersymmetric) AdS 7 solutions of [38] with pure D6-branes (ie not all 6-brane charges come from spherical D8-branes). The logic of this paper suggests that those AdS 7 backgrounds, despite supersymmetry, are unstable and not suitable candidates for holographic descriptions of certain 6D CFT's. The typical 6D CFT's classified by these massive IIA solutions are however described mainly in terms of spherical 8-branes [44], but it would be interesting to understand this better. 7 We have mentioned several times that a boundary condition for the gauge fields exists that allows a solution without infinite flux clumping. This boundary condition is simply ΦM = 0 and cannot exist for D6-branes as shown in [33].

A O6 planes in flat space
The O6-plane metric is given by We expand this metric around the critical radius The result is

B Computing the probe potential
We now turn to calculating the probe potential. In order to obtain the expression for B 6 we make use of the Bianchi identity d(H 7 − F 1 C 6 ) = 0, from which we infer The RR gauge potential C 6 can be calculated from the expression for F 7 = dC 6 as

JHEP03(2017)141
C The force between an NS5/KK5 pair We can compute the interaction potential between the two probe NS5-branes by placing one NS5-brane in the geometry, let it backreact, and evaluate that backreaction in the probe brane action of the other NS5-brane. Hence, to find the force between an NS5/NS5 pair amounts to solving a differential equation equation for the harmonic function, denoted K, describing a NS5-brane in our throat geometry. We are only interested in the probe potential and thereby the interaction term close to r = 0, which simplifies our task immensely. The background metric at r = 0 simplifies to a cylinder: ds 2 → v −1 ds 2 6 + v(dψ 2 + dr 2 + r 2 dΩ 2 2 ) . where we have put the NS5-brane at the location r = 0 and ∆ψ measures the distance from the brane along the isometry direction. Furthermore we have put the number of NS5-branes to one and used that in our conventions µ 5 = 2π. In the region around r = 0 we have B 6 = 1 v 2 g 2 s sin 2 (π∆ψ/ s ) 4π 4 + sin 2 (π∆ψ/ s ) . (C.8) We obtain the interaction term by inserting this into the action of the other NS5-brane so that we get − µ 5 e −2φ √ −g + B 6 = −2 µ 5 v 2 g 2 s sin 2 (π∆ψ/ s ) 4π 4 + sin 2 (π∆ψ/ s ) vol 6 . (C.9) The potential describing the force starts rising until it reaches a maximum and then goes back in an identical matter. In reference [15] this potential was approximated to be piecewise linear, which is qualitatively similar. The reason the potential grows roughly linear with ψ is due to the throat geometry; the brane/antibrane forces are effectively confined to the circle dimension such that it mimics electrodynamics in 1 + 1 dimensions.

JHEP03(2017)141
D Computing the corrected probe potential To evaluate the probe potential of a single NS5/NS5-pair in the background described in section 5.3 we require the gauge field expressions: To calculate the interaction term we again follow the same procedure as in appendix C. The metric is approximated by the cylinder metric evaluated at a given radius r = s ds 2 = e 2A ds 2 6 + e 2B (e −2B−2A dψ 2 + dr 2 + r 2 dΩ 2 2 ) , (D.7) We rescale r to absorb the factors of e B−A asr = e B−A r. Using the method described in appendix C we find the interaction potential where we threw away an additive constant. If we were to taker to zero we would find that λ blows up in exactly the right way, see (3.11), to compensate for the vanishing prefactor, and hence we obtain finite result for the terms that include λ. All other terms, including the interaction term, vanish in this limit, which physically can be understood as an effect due to redshift. Since we cannot trust these expressions below the string scale, we need to evaluate them atr = s to investagate whether the redshift is large enough within the regime under our control. In paricular, we find that λ at its cut-off value is λ = g s p, which, indeed, is large.

JHEP03(2017)141
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.