This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

11institutetext: STAG Research Centre, Physics and Astronomy, University of Southampton,
Highfield, Southampton SO17 1BJ, United Kingdom
22institutetext: Department of Chemistry and Physics, SUNY Old Westbury, Old Westbury, NY, United States33institutetext: INFN, sezione di Milano-Bicocca, Piazza della Scienza 3, I-20126 Milano, Italy

From Large to Small 𝓝=(𝟒,𝟒)\mathcal{N}=(4,4) Superconformal Surface Defects in Holographic 6d SCFTs

Pietro Capuozzo 2    ​​, John Estes 3    ​​, Brandon Robinson 1    ​​, and Benjamin Suzzoni [email protected] [email protected] [email protected] [email protected]
Abstract

Two-dimensional (2d) 𝒩=(4,4)\mathcal{N}=(4,4) Lie superalgebras can be either “small” or “large”, meaning their R-symmetry is either 𝔰𝔬(4)\mathfrak{so}(4) or 𝔰𝔬(4)𝔰𝔬(4)\mathfrak{so}(4)\oplus\mathfrak{so}(4), respectively. Both cases admit a superconformal extension and fit into the one-parameter family 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}\left(2,1;\gamma\right)\oplus\mathfrak{d}\left(2,1;\gamma\right), with parameter γ(,)\gamma\in(-\infty,\infty). The large algebra corresponds to generic values of γ\gamma, while the small case corresponds to a degeneration limit with γ\gamma\to-\infty. In 11d supergravity, we study known solutions with superisometry algebra 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}\left(2,1;\gamma\right)\oplus\mathfrak{d}\left(2,1;\gamma\right) that are asymptotically locally AdS×7𝕊4{}_{7}\times\mathds{S}^{4}. These solutions are holographically dual to the 6d maximally superconformal field theory with 2d superconformal defects invariant under 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}\left(2,1;\gamma\right)\oplus\mathfrak{d}\left(2,1;\gamma\right). We show that a limit of these solutions, in which γ\gamma\to-\infty, reproduces another known class of solutions, holographically dual to small 𝒩=(4,4)\mathcal{N}=(4,4) superconformal defects. We then use this limit to generate new small 𝒩=(4,4)\mathcal{N}=(4,4) solutions with finite Ricci scalar, in contrast to the known small 𝒩=(4,4)\mathcal{N}=(4,4) solutions. We then use holography to compute the entanglement entropy of a spherical region centered on these small 𝒩=(4,4)\mathcal{N}=(4,4) defects, which provides a linear combination of defect Weyl anomaly coefficients that characterizes the number of defect-localized degrees of freedom. We also comment on the generalization of our results to include 𝒩=(0,4)\mathcal{N}=(0,4) surface defects through orbifolding.

arxiv: 2402.11745

1 Introduction

Six dimensional (6d) superconformal field theories (SCFTs) hold a special place among quantum field theories (QFTs). Owing to the classification discovered in the seminal work by Nahm Nahm:1977tg , superconformal symmetry is only possible in six and fewer spacetime dimensions. Moreover, 𝒩=(2,0)\mathcal{N}=(2,0) is the maximal amount of supersymmetry (SUSY) that a 6d theory can have. Combining this amount of SUSY with conformal symmetry constrains a 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT to such a degree that the only additional information that is necessary to completely determine the theory is the choice of a gauge algebra.

The study of the 6d 𝒩=(2,0)\mathcal{N}=(2,0) theory is thus of fundamental importance in QFT, for many reasons. For example, the 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT determines the physics of many other QFTs in 6d, via SUSY-breaking deformations Antoniadis:1998ki ; Antoniadis:1998ep ; Cordova:2016xhm such as orbifolds Blum:1997fw ; Blum:1997mm ; Brunner:1997gf ; Intriligator:1997dh . By suitable (partial) topological twisting, the 6d theory compactified on, e.g., a Riemann surface Gaiotto:2009we or a 3-manifold Dimofte:2011ju can also determine the physics of infinite families of QFTs in d<6d<6.

The 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT is also of fundamental importance in quantum gravity. Currently, the leading candidate for an ultra-violet (UV)-complete theory of quantum gravity is M-theory. M-theory’s fundamental objects are M2-branes Duff:1990xz and M5-branes Gueven:1992hh , and the low-energy worldvolume theory on MM coincident M5-branes is the 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT with gauge algebra AM1A_{M-1} Strominger:1995ac . Understanding the 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT is thus essential to understanding M-theory in general. In particular, via the Anti-de Sitter/CFT (AdS/CFT) correspondence, the 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT can provide a fully non-perturbative definition of M-theory on an asymptotically 7d AdS spacetime, AdS7, times a four-sphere, 𝕊4\mathds{S}^{4} Witten:1996hc ; Maldacena:1997re ; Witten:1998wy .

Strongly interacting SCFTs constructed in string- and M-theory, including the non-Abelian 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT, are prohibitively difficult to study, for many reasons, of which we will mention only three. First, the 𝒩=(1,0)\mathcal{N}=(1,0) and 𝒩=(2,0)\mathcal{N}=(2,0) SUSY multiplets include a chiral two-form gauge field, and writing a local, gauge-invariant Lagrangian for a non-Abelian higher-form gauge field remains a major open problem. These SCFTs thus have no known Lagrangian descriptions. Second, in the space of renormalization group (RG) flows, these SCFTs are isolated fixed points, and in particular they cannot be reached as infra-red (IR) fixed points of RG flows from free ultra-violet (UV) fixed points. Third, these SCFTs are intrinsically strongly interacting. For example, the 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT has no dimensionless parameter besides MM that can be tuned to allow a perturbative expansion.

As a result, practically all of our direct knowledge111Indirect methods such as dimensional reduction to 5d 𝒩2\mathcal{N}\leq 2 SUSY QFTs have also been used to great effect to study these theories, e.g. by using the resulting lower dimensional Lagrangian description together with supersymmetric localization techniques Bullimore:2014upa . of interacting 6d SCFTs comes from non-perturbative methods, such as the superconformal bootstrap Beem:2015aoa , F-theory Heckman:2015bfa , and especially AdS/CFT Apruzzi:2013yva , where holographic computations of quantities like Weyl anomalies Henningson:1998gx and entanglement entropy (EE) Hung:2011xb are used to great effect to characterize 6d SCFTs at large MM.

An aspect of 6d SCFTs, and generally of QFTs in three and higher dimensions, that requires particularly careful treatment to characterize is the spectrum of 2d, string-like or surface, defects. In the co-dimension one case, 2d defects in 3d QFTs arise as boundaries or interfaces, and so are more easily studied and, thus, more familiar than their higher co-dimension realizations. Despite being somewhat more exotic in standard treatments of QFTs, 2d defects of co-dimension two and greater show up in a number of settings222This is by no means an exhaustive list of the work done on 2d defects. For a recent review of boundaries and defects in QFTs and further references on the topic, see Andrei:2018die .: from free field theories Soderberg:2017oaa ; Lauria:2020emq ; Giombi:2021uae ; Bianchi:2021snj to strongly interacting and non-Lagrangian 4d QFTs, e.g. Gukov:2006jk ; Gukov:2008sn ; Alday:2009fs , to being fundamental objects in 6d SCFTs Ganor:1996nf and in the study of EE and Renyi entropies Hung:2014npa ; Lewkowycz:2014jia . As such, the last few decades have seen tremendous advancements in characterizing Billo:2016cpy and constraining Jensen:2015swa ; Kobayashi:2018lil the properties of 2d defects, and it is vitally important in the study of QFTs, generally, to continue this effort by finding novel constructions of surface defects and examining their unique physics.

Of interest to us in the current work are the holographic descriptions, afforded by AdS/CFT, of 6d SCFTs and the defects that they support. In particular, we will primarily focus our attention on solutions to 11d supergravity (SUGRA) that are contained in a one-parameter family of solutions with superisometry given by the exceptional Lie superalgebra 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}(2,1;\gamma)\oplus\mathfrak{d}(2,1;\gamma) Bachas:2013vza . Crucially, an asymptotically AdS×7𝕊4{}_{7}\times\mathds{S}^{4} solution is possible only at certain values of γ\gamma. Indeed, as a historical note, prior to the full classification given in Bachas:2013vza , evidence of 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}(2,1;\gamma)\oplus\mathfrak{d}(2,1;\gamma) invariant SUGRA solutions that exist for general γ\gamma were found in superconformal Janus solutions in 4d4d gauged 𝒩=8\mathcal{N}=8 SUGRA Bobev:2013yra , which extended beyond the known γ=1\gamma=1 AdS×4𝕊7{}_{4}\times\mathds{S}^{7} Janus solution of 11d SUGRADHoker:2009lky .

More pertitent to the cases that we will study below, the most well studied case among the values of γ\gamma that admit asymptotically locally AdS×7𝕊4{}_{7}\times\mathds{S}^{4} solutions is γ=1/2\gamma=-1/2, which holographically describes 1/2-BPS Wilson surface operators in the 6d AM1A_{M-1} 𝒩=(2,0)\mathcal{N}=(2,0) SCFT that preserve a large 𝒩=(4,4)\mathcal{N}=(4,4) 2d SUSY. These BPS Wilson surface operators have a long history in M-theory descriptions of 6d SCFTs Strominger:1995ac ; Ganor:1996nf ; Howe:1997ue ; Maldacena:1998im ; Berenstein:1998ij , and there has been a recent resurgence of interest in these defects where holographic Gentle:2015jma ; Estes:2018tnu ; Jensen:2018rxu and field theoretic Chalabi:2020iie ; Drukker:2020atp ; Drukker:2020dcz computations have characterized these defect CFTs through their EE and Weyl anomalies.

Recently, new solutions in 11d SUGRA have been constructed that are proposed to be holographically dual to 2d BPS surface defects in 6d 𝒩=(1,0)\mathcal{N}=(1,0) SCFTs preserving “small” 𝒩=(4,4)\mathcal{N}=(4,4) or 𝒩=(0,4)\mathcal{N}=(0,4) 2d SUSY Faedo:2020nol . In the sections below, we will clearly demonstrate that these new solutions also fit into the one-parameter family of solutions in Bachas:2013vza in the limit where γ\gamma\to-\infty. It will turn out that the solutions in Faedo:2020nol are, in fact, a particular case within a broader class of solutions that can be realized in the γ\gamma\to-\infty limit, which we will characterize by computing their contributions to the EE of a spherical region.

A crucial point that we emphasize in our construction of the γ\gamma\to-\infty solutions is that the superisometry algebra of the near-horizon limit of a stack of M5-branes, 𝔬𝔰𝔭(8|4)\mathfrak{osp}(8^{*}|4), does not contain 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}\left(2,1;\gamma\right)\oplus\mathfrak{d}\left(2,1;\gamma\right) as a sub-superalgebra frappat1996dictionary ; D_Hoker_2008 ; D_Hoker_2008b . In the supergravity solution, this is manifested by the appearance of additional terms in the four-form flux as compared to that due to pure M5-branes. Nevertheless, the geometry is still locally asymptotically AdS×7𝕊4{}_{7}\times\mathds{S}^{4}, which is in line with the fact that the bosonic subalgebra of 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}\left(2,1;\gamma\right)\oplus\mathfrak{d}\left(2,1;\gamma\right) is a subalgebra of 𝔬𝔰𝔭(8|4)\mathfrak{osp}(8^{*}|4).

One upshot of our analysis that follows from the identification of the global symmetries in the γ\gamma\to-\infty limit is that it allows for a suitable regulation scheme, which we will employ when we compute the defect sphere EE. Lacking a generalized program of holographic renormalization for SUGRA solutions dual to defects on the field theory side, we will use a background subtraction scheme in order to remove the ambient degrees of freedom and isolate contributions to physical quantities coming from the defect. The key step in this background subtraction scheme is the correct identification of the vacuum solution, and as will be made clear, the ambient theory with a “trivial defect” is a deformed 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT preserving the bosonic superconformal subalgebra 𝔰𝔬(2,2)𝔰𝔬(4)𝔰𝔬(5)𝔬𝔰𝔭(8|4){\mathfrak{so}}(2,2)\oplus{\mathfrak{so}}(4)\oplus{\mathfrak{so}}(5)\subset{\mathfrak{osp}}(8*|4). Since the solutions in Faedo:2020nol belong to the class of solutions that we study below, in finding a vacuum solution that manifests the corresponding isometries and carefully carrying out background subtraction, we will also resolve a puzzle in Faedo:2020nol , where physical quantities like the “defect central charge” were divergent.

Ultimately, we will see that that universal part of the defect sphere EE, SEE(univ)S_{\rm EE}^{\tiny\rm(univ)}, i.e. the coefficient of its log-divergent part, is determined in terms of the highest weight vector, ϖ\varpi, of an AM1A_{M-1} irreducible representation determined by a Young diagram that encodes the partition of M5-branes that specifies the defect. Explicitly, we will show that

SEE(univ)=(ϖ,ϖ)5,\displaystyle S_{\rm EE}^{\tiny\rm(univ)}=-\frac{(\varpi,\varpi)}{5}~{}, (1)

where (,)(\cdot,\cdot) is the scalar product on the weight space. This result is similar to the Wilson surface sphere EE Estes:2018tnu in that both are expressible in terms of scalar quantities derived from representation data but differ in that eq. (1) is negative definite333Unlike in an ordinary 2d CFT, eq. (1) being strictly negative does not necessarily signal that the theory may be non-unitary. Indeed, the 2d defect sphere EE is expressible as a signed linear combination of defect Weyl anomaly coefficients Jensen:2018rxu , which is not bounded from below. and is completely determined by the highest weight vector.

In section 2, we begin by reviewing the 11d supergravity solutions dual to 2d small 𝒩=(4,4)\mathcal{N}=(4,4) defects in 6d 𝒩=(1,0)\mathcal{N}=(1,0) SCFTs found in Faedo:2020nol . In section 3 we briefly review the 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}(2,1;\gamma)\oplus\mathfrak{d}(2,1;\gamma)-invariant solutions to 11d supergravity found in Bachas:2013vza , and we show that by orbifolding the solutions in the γ\gamma\to-\infty limit, we can recover the solutions in Faedo:2020nol . We then use the γ\gamma\to-\infty limit to construct new 2d small 𝒩=(4,4)\mathcal{N}=(4,4) defects with finite Ricci scalar in section 4. Further in section 4, we demonstrate that the naïve AdS×7𝕊4{}_{7}\times\mathds{S}^{4} vacuum is inappropriate to use in a background subtraction scheme for regulating holographic computations in the γ\gamma\to-\infty limit, and we identify the correct background to use in this scheme. In section 5, we utilize the new γ\gamma\to-\infty supergravity solutions and correct regulating scheme in a computation of the contribution of a flat 2d small 𝒩=(4,4)\mathcal{N}=(4,4) superconformal defect to the EE of a spherical region in a 6d SCFT. We then summarize our findings and discuss remaining issues and open questions surrounding these new defect solutions in section 6.

In addition, there are two appendices that detail technical aspects of the computations in the main text. First, in appendix A, we analyze the asymptotic expansion of the supergravity data that specify the new solutions in the γ\gamma\to-\infty limit and construct the map to Fefferman-Graham (FG) gauge. Lastly, in appendix B, we carefully treat the integral in the area functional of the Ryu-Takayanagi (RT) surface in the computation of the holographic EE of the defect in the dual field theory.

2 Review: small 𝓝=(𝟒,𝟒){\mathcal{N}=(4,4)} surface defects

1,1\mathbb{R}^{1,1} rr θ1\theta^{1} θ2\theta^{2} χ\chi zz ρ\rho φ1\varphi^{1} φ2\varphi^{2} ϕ\phi
KK
\cdot \cdot \cdot ISO
M5
\cdot \cdot \cdot \cdot \cdot
M2
\cdot \cdot \cdot \cdot \sim \sim \sim \sim
M5
\cdot \cdot \cdot \cdot \sim
KK
\cdot \cdot \cdot ISO
Table 1: The 1/8-BPS brane setup of Faedo:2020nol , with M2-M5 defect branes intersecting orthogonal M5-branes, and with both stacks of 5-branes probing A-type singularities. In our conventions, –, \cdot and \sim denote directions along which a brane is extended, localized, and smeared, respectively, while ISO(metric) denotes the compact direction of the KK-monopoles.

We begin with a brief summary of the results of Faedo:2020nol . The particular 11d supergravity metric constructed therein is the uplift of a 7d charged AdS×3𝕊3{}_{3}\times\mathds{S}^{3} domain wall initially found in Dibitetto:2017tve , and is given by

ds2=4kQM5HM51/3(dsAdS32+ds𝕊3/k2)+HM52/3(dz2+dρ2+ρ2ds𝕊~3/k2),\displaystyle\text{d}s^{2}=4kQ_{\text{M5}}H_{\text{M5}^{\prime}}^{-1/3}\left(\text{d}s^{2}_{\text{AdS}_{3}}+\text{d}s^{2}_{\mathds{S}^{3}/\mathbb{Z}_{k}}\right)+H_{\text{M5}^{\prime}}^{2/3}\left(\text{d}z^{2}+\text{d}\rho^{2}+\rho^{2}\text{d}s^{2}_{\tilde{\mathds{S}}^{3}/\mathbb{Z}_{k^{\prime}}}\right), (2)

for some parameter QM5Q_{\text{M5}} and a function HM5H_{\text{M5}^{\prime}} defined over a 4d space parametrized by the coordinates {z,ρ,φ1,φ2}\{z,\rho,\varphi^{1},\varphi^{2}\}. The solution above captures the near-horizon geometry of the brane intersection depicted in table 1. Namely, a “bound state” (in the sense discussed in footnote 4) of M2- and M5-branes, with charges QM2Q_{\text{M2}} and QM5Q_{\text{M5}}, intersects an orthogonal stack of M5-branes, thus forming a 1/4-BPS brane setup. In 2d notation, this corresponds to 𝒩=(4,4)\mathcal{N}=(4,4) supersymmetry, with the large R-symmetry realized geometrically as the isometry of the two 3-spheres 𝕊3\mathds{S}^{3} and 𝕊~3\tilde{\mathds{S}}^{3} with coordinates {χ,θ1,θ2}\{\chi,\theta^{1},\theta^{2}\} and {ϕ,φ1,φ2}\{\phi,\varphi^{1},\varphi^{2}\}, respectively. In addition, the two stacks of 5-branes can be made to probe ALE singularities by introducing two Kaluza-Klein (KK) monopoles, with charges kk and kk^{\prime} and Taub-NUT directions χ\partial_{\chi} and ϕ\partial_{\phi}, respectively. The inclusion of any one of the two KK monopoles results in a further breaking of the preserved supersymmetries and a degeneration to an 1/8-BPS setup, which in 2d language corresponds to (large) 𝒩=(0,4)\mathcal{N}=(0,4) supersymmetry. The presence of the second KK monopole does not incur a further loss of supersymmetry, so the final brane configuration is always at least a 1/81/8-BPS supergravity solution.

Furthermore, the defect M2- and M5-branes are fully localized within a 2d submanifold of the worldvolume of the M5-branes. This is to be expected from the holographic realization of a surface defect in a 6d SCFT; this interpretation was first attached to the 7d domain wall in Dibitetto:2017klx and to the full 11d SUGRA background in Faedo:2020nol . On the other hand, the defect branes are smeared in the directions transverse to the worldvolume of the M5-branes444It is this property – that the M2- and M5-branes do not share transverse directions with the M5-branes, other than those along which the former are smeared – which we are implicitly referring to when we describe the M2- and M5-branes as forming a “bound state”. This aligns with the terminology used in Faedo:2020nol , and should not to be confused with the dyonic supermembrane Izquierdo:1995ms , which is a rather different multimembrane solution of 11d supergravity. Indeed, in the setup described by table 1, there is no M2-brane charge dissolved within the M5-brane worldvolume, nor do the M2-branes polarize into a fuzzy 3-sphere via the Myers effect., so that their charge is localized within 𝕊3/k\mathds{S}^{3}/\mathbb{Z}_{k}, but not 𝕊~3/k\tilde{\mathds{S}}^{3}/\mathbb{Z}_{k^{\prime}}. Therefore, while the metric in eq. (2) manifests the isometry groups of both (orbifolded) 3-spheres, the R-symmetry is partially broken, giving rise to small 𝒩=(0,4)\mathcal{N}=(0,4) supersymmetry. For k=1k^{\prime}=1, the solution above fits into the classification of 𝒩=(0,4)\mathcal{N}=(0,4) AdS×3𝕊3/k×CY2{}_{3}\times\mathds{S}^{3}/\mathbb{Z}_{k}\times\operatorname{CY}_{2} backgrounds foliated over an interval performed in Lozano:2020bxo , where CY2=4\operatorname{CY}_{2}=\mathbb{R}^{4} contains the (round) 𝕊~3\tilde{\mathds{S}}^{3}. In particular, the solution above corresponds to taking the M5-branes to be completely localized in their transverse space.

On shell, the defect brane charges are constrained to be equal, QM2=QM5Q_{\text{M2}}=Q_{\text{M5}}, while the function HM5H_{\text{M5}^{\prime}} satisfies Faedo:2020nol

ρ^32HM5(z,ρ^)+kz2HM5(z,ρ^)ρ^=0,\displaystyle\nabla^{2}_{\mathbb{R}^{3}_{\hat{\rho}}}H_{\text{M5}^{\prime}}(z,\hat{\rho})+\frac{k^{\prime}\partial_{z}^{2}H_{\text{M5}^{\prime}}(z,\hat{\rho})}{\hat{\rho}}=0, (3)

where we rescaled ρ^=ρ2/(4k)\hat{\rho}=\rho^{2}/(4k^{\prime}), and denoted by ρ^3\mathbb{R}^{3}_{\hat{\rho}} the three-dimensional subspace which is transverse to the M5-branes, parallel to the M5-branes, and along which the M2-branes are smeared. Following the parametrization adopted in table 1, ρ^3\mathbb{R}^{3}_{\hat{\rho}} is then the space spanned by {ρ^,φ1,φ2}\{\partial_{\hat{\rho}},\partial_{\varphi^{1}},\partial_{\varphi^{2}}\}, with ρ^\hat{\rho} being the radial coordinate and {φ1,φ2}\{\varphi^{1},\varphi^{2}\} parametrizing a 2-sphere.

In terms of the brane setup described above, then, the spacetime in eq. (2) is reached by approaching the brane intersection locus from within the worldvolume of the M5-branes in a radial fashion, i.e. by taking r0r\to 0. In this limit, the ISO(1,1)ISO(1,1) isometry group gets promoted to SO(2,2)SO(2,2), and the M5-brane worldvolume becomes AdS×3𝕊3/k{}_{3}\times\mathds{S}^{3}/\mathbb{Z}_{k}.

A particular solution to eq. (3) is, for any α\alpha\in\mathbb{R},

HM5(z,ρ)=42g31P+PP+2+P24α2+2P+PP+2+P2+2P+P,H_{\text{M5}^{\prime}}(z,\rho)=\frac{4\sqrt{2}}{g^{3}}\frac{1}{P_{+}P_{-}}\frac{\sqrt{P_{+}^{2}+P_{-}^{2}-4\alpha^{2}+2P_{+}P_{-}}}{P_{+}^{2}+P_{-}^{2}+2P_{+}P_{-}}, (4)

where P±=z2+(ρ±α)2P_{\pm}=\sqrt{z^{2}+(\rho\pm\alpha)^{2}}. By redefining

ρ\displaystyle\rho =αcosξ1μ5\displaystyle=\alpha\frac{\cos\xi}{\sqrt{1-\mu^{5}}} (5a)
z\displaystyle z =αμ52sinξ1μ5\displaystyle=\alpha\mu^{\frac{5}{2}}\frac{\sin\xi}{\sqrt{1-\mu^{5}}} (5b)

and setting α=(27/4g3/2kQM5)1\alpha=(2^{7/4}g^{3/2}kQ_{\text{M5}})^{-1}, the particular solution in eq. (4) can be matched to the one found in eq. (2.17) of Faedo:2020nol , which we now reproduce for clarity555To see this, it is convenient to first rationalize the product P+PP_{+}P_{-} as P+P=α2(1+μ5(1μ5)cos(2ξ))/2(1μ5)P_{+}P_{-}=\alpha^{2}\big{(}1+\mu^{5}-(1-\mu^{5})\cos(2\xi)\big{)}/2(1-\mu^{5}).:

HM5(μ,ξ)=227/4(gkQM5)3μ5/2(1μ5)3/2μ5cos2ξ+sin2ξ.H_{\text{M5}^{\prime}}(\mu,\xi)=2^{27/4}(\sqrt{g}kQ_{\text{M5}})^{3}\frac{\mu^{5/2}(1-\mu^{5})^{3/2}}{\mu^{5}\cos^{2}\xi+\sin^{2}\xi}~{}. (6)

Furthermore, in Faedo:2020nol it was argued that, as μ1\mu\to 1 (which corresponds to a non-linear limit in the original coordinates ρ^\hat{\rho} and zz), the near-horizon geometry in eq. (2) locally recovers the AdS/7k×𝕊4/k{}_{7}/\mathbb{Z}_{k}\times\mathds{S}^{4}/\mathbb{Z}_{k^{\prime}} vacuum of M-theory. This was shown by realizing the 11d line element as the uplift of the domain wall solution to 𝒩=1\mathcal{N}=1, d=7d=7 supergravity (whose gauge coupling constant gg appears in eq. (6) above) found in Dibitetto:2017tve , which interpolates between AdS7 (as μ1\mu\to 1) and an infrared singularity (at μ=0\mu=0).

While the singular nature of the solution is not immediately obvious, it can be made manifest by studying the Ricci scalar, \mathcal{R}, in the z0z\to 0 limit. For metrics generally of the form of eq. (2), and following D_Hoker_2008b , we can write the expression for the Ricci scalar for an arbitrary harmonic function HM5H_{\text{M5}^{\prime}} as

\displaystyle\mathcal{R} =HM52/3[16(zHM5)2+(ρHM5)2HM5223z2HM5+ρ2HM5HM52ρHM5ρHM5].\displaystyle=H_{\rm M5^{\prime}}^{-2/3}\left[\frac{1}{6}\frac{(\partial_{z}H_{\rm M5^{\prime}})^{2}+(\partial_{\rho}H_{\rm M5^{\prime}})^{2}}{H_{\rm M5^{\prime}}^{2}}-\frac{2}{3}\frac{\partial_{z}^{2}H_{\rm M5^{\prime}}+\partial_{\rho}^{2}H_{\rm M5^{\prime}}}{H_{\rm M5^{\prime}}}-2\frac{\partial_{\rho}H_{\rm M5^{\prime}}}{\rho H_{\rm M5^{\prime}}}\right]. (7)

The function HM5H_{\rm M5^{\prime}} has a branch point located at z=0z=0 and ρ=α\rho=\alpha and correspondingly admits two different expansions as z0z\rightarrow 0, depending on whether ρ\rho is larger or smaller than α\alpha. The choice of sign for the branch cut can be determined by the requirement P+P0P_{+}P_{-}\geq 0. For ρ>α\rho>\alpha, this leads to

=g2(2α23ρ2)212ρ2/3(ρ2α2)5/3+𝒪(z2)\displaystyle\mathcal{R}=\frac{g^{2}(2\alpha^{2}-3\rho^{2})^{2}}{12\rho^{2/3}(\rho^{2}-\alpha^{2})^{5/3}}+\mathcal{O}(z^{2}) (8)

which has a pole as ρα\rho\rightarrow\alpha. This corresponds to setting ξ=0\xi=0, with the pole appearing as μ0\mu\rightarrow 0. For ρ<α\rho<\alpha, we find

=g2α2/3(α2ρ2)12z8/3+𝒪(1z2/3)\displaystyle\mathcal{R}=\frac{g^{2}\alpha^{2/3}(\alpha^{2}-\rho^{2})}{12z^{8/3}}+\mathcal{O}\bigg{(}\frac{1}{z^{2/3}}\bigg{)} (9)

which corresponds to taking μ0\mu\rightarrow 0 with ξ0\xi\neq 0. Thus for all values of ξ\xi, we find that the Ricci scalar diverges as μ0\mu\rightarrow 0.

3 From large to small 𝓝=(𝟒,𝟒)\mathcal{N}=(4,4) surface defects

To begin, we will briefly review the classification due to D_Hoker_2008b ; Bachas:2013vza of the d=11d=11 supergravity solutions with superisometry given by two copies of the exceptional Lie superalgebra 𝔡(2,1;γ){\mathfrak{d}}(2,1;\gamma). This classification represents the foundation for the new defect solutions we present below.

In particular, we will be interested in the real form 𝔡(2,1;γ;0)\mathfrak{d}(2,1;\gamma;0) which arises as a real subsuperalgebra of the fixed points of an involutive semimorphism of 𝔡(2,1;γ)\mathfrak{d}(2,1;\gamma) frappat1996dictionary . Recall that the 9-dimensional maximal bosonic subalgebra of the real form 𝔡(2,1;γ;0){\mathfrak{d}}(2,1;\gamma;0) is 𝔰𝔬(2,1)𝔰𝔬(3)𝔰𝔬(3){\mathfrak{so}}(2,1)\oplus{\mathfrak{so}}(3)\oplus{\mathfrak{so}}(3). Labelling each factor in this subalgebra by an index a{1,2,3}a\in\{1,2,3\}, the generators T(a)T^{(a)} satisfy

[Ti(a),Tj(b)]=iδabεijkηaklTl(a)for i,j{1,2,3},\displaystyle[T^{(a)}_{i},\,T^{(b)}_{j}]=i\delta^{ab}\varepsilon_{ijk}\eta^{kl}_{a}T^{(a)}_{l}\quad\text{for }i,j\in\{1,2,3\}~{}, (10)

where εijk\varepsilon_{ijk} and ηakl\eta^{kl}_{a} are the totally anti-symmetric tensor (with ε123=1\varepsilon_{123}=1) and the canonical metric induced by the Killing form, respectively. In addition to the bosonic sector, 𝔡(2,1;γ){\mathfrak{d}}(2,1;\gamma) contains an 8-dimensional fermionic generator with components FA1A2A3F_{A_{1}A_{2}A_{3}}, where the indices Aa{±}A_{a}\in\{\pm\} transform in the spinorial representation 𝟐\boldsymbol{2} of the atha^{\rm th} factor in the even subalgebra. Furthermore, the components FA1A2A3F_{A_{1}A_{2}A_{3}} obey the following anti-commutation relation

{FA1A2A3,FB1B2B3}=\displaystyle\{F_{A_{1}A_{2}A_{3}},F_{B_{1}B_{2}B_{3}}\}=\ β1CA2B2CA3B3(Cσi)A1B1Ti(1)\displaystyle\beta_{1}C_{A_{2}B_{2}}C_{A_{3}B_{3}}(C\sigma^{i})_{A_{1}B_{1}}T_{i}^{(1)} (11)
+β2CA1B1CA3B3(Cσi)A2B2Ti(2)\displaystyle+\beta_{2}C_{A_{1}B_{1}}C_{A_{3}B_{3}}(C\sigma^{i})_{A_{2}B_{2}}T_{i}^{(2)}
+β3CA1B1CA2B2(Cσi)A3B3Ti(3),\displaystyle+\beta_{3}C_{A_{1}B_{1}}C_{A_{2}B_{2}}(C\sigma^{i})_{A_{3}B_{3}}T_{i}^{(3)}~{},

where Ciσ2C\equiv i\sigma^{2} is the charge conjugation matrix, {σi}i=13\{\sigma^{i}\}_{i=1}^{3} are the Pauli matrices, and {βi}i=13\{\beta_{i}\}_{i=1}^{3} are real parameters satisfying a=13βa=0\sum_{a=1}^{3}\beta_{a}=0, which follows from the (generalized) Jacobi identity. This last constraint, together with the possibility of absorbing any rescaling {λβ1,λβ2,λβ3}\{\lambda\beta_{1},\lambda\beta_{2},\lambda\beta_{3}\} for λ\lambda\in\mathbb{C} into a redefinition of the normalization of the fermionic generator, implies that 𝔡(2,1;γ){\mathfrak{d}}(2,1;\gamma) is entirely specified by the choice of a ratio of any two of the three β\beta parameters; here, we take γβ2/β3\gamma\equiv\beta_{2}/\beta_{3}. Note that 𝔡(2,1;γ){\mathfrak{d}}(2,1;\gamma) is the only (finite-dimensional) Lie superalgebra admitting a continuous parametrization.

Amongst the possible values that γ\gamma can take, there are clearly three special values corresponding to the vanishing of any one of the β\beta parameters. Specifically, choosing β1=0\beta_{1}=0 fixes γ=1\gamma=-1. The more interesting case, and the one relevant to our analysis, is β3=0\beta_{3}=0, which corresponds to γ±\gamma\to\pm\infty. The case β2=0\beta_{2}=0 corresponds to γ=0\gamma=0 and is equivalent under a group involution as discussed at the end of this section. In the limit β3=0\beta_{3}=0, the anticommutator in eq. (11) degenerates to

{FA1A2A3,FB1B2B3}=β1(CA2B2CA3B3(Cσi)A1B1Ti(1)CA1B1CA3B3(Cσi)A2B2Ti(2)).\{F_{A_{1}A_{2}A_{3}},F_{B_{1}B_{2}B_{3}}\}=\beta_{1}\left(C_{A_{2}B_{2}}C_{A_{3}B_{3}}(C\sigma^{i})_{A_{1}B_{1}}T_{i}^{(1)}-C_{A_{1}B_{1}}C_{A_{3}B_{3}}(C\sigma^{i})_{A_{2}B_{2}}T_{i}^{(2)}\right)~{}. (12)

Consequently, the 𝔰𝔬(4)R𝔰𝔬(3)𝔰𝔬(3)\mathfrak{so}(4)_{R}\cong\mathfrak{so}(3)\oplus\mathfrak{so}(3) R-symmetry of the large superalgebra is broken into the single 𝔰𝔬(3)R\mathfrak{so}(3)_{R} factor which constitutes the R-symmetry of a small superalgebra. Note that the other 𝔰𝔬(3)\mathfrak{so}(3) factor remains a bosonic symmetry of the supergravity solution; however, it is now realized as a flavor symmetry, rather than as an outer automorphism of the supersymmetry algebra.

In addition, there are isolated points of interest in the γ\gamma-parameter space where the real form 𝔡(2,1;γ;0)\mathfrak{d}(2,1;\gamma;0) reduces to classical Lie superalgebras:

𝔡(2,1;γ;0)\displaystyle\mathfrak{d}(2,1;\gamma;0) =𝔬𝔰𝔭(4|2)\displaystyle=\mathfrak{osp}(4^{*}|2)\quad for γ{2,1/2}\displaystyle\text{for }\gamma\in\{-2,-1/2\} (13a)
𝔡(2,1;γ;0)\displaystyle\mathfrak{d}(2,1;\gamma;0) =𝔬𝔰𝔭(4|2;)\displaystyle=\mathfrak{osp}(4|2;\mathbb{R})\quad for γ=1.\displaystyle\text{for }\gamma=1~{}. (13b)

In particular, γ=1/2\gamma=-1/2 is the only value (up to algebra involutility, as we will discuss shortly) for which 𝔡(2,1;γ;0)𝔡(2,1;γ;0)\mathfrak{d}(2,1;\gamma;0)\oplus\mathfrak{d}(2,1;\gamma;0) admits a canonical inclusion into the superisometry algebra 𝔬𝔰𝔭(8|4)\mathfrak{osp}(8^{*}|4) of AdS×7𝕊4{}_{7}\times\mathds{S}^{4}. This case was studied extensively in Estes:2018tnu ; it is the holographic realization of Wilson surfaces, 1/2-BPS codimension-4 superconformal solitons, within the 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT. In this case, the ambient 6d SCFT is undeformed, in the sense that the supergroup of symmetries preserved by the defect is a subgroup of the 6d superconformal symmetry group. This is a common feature throughout the study of defects embedded in higher-dimensional theories.

For generic values of γ\gamma, including the limit γ±\gamma\to\pm\infty, the superisometry algebra 𝔡(2,1;γ;0)𝔡(2,1;γ;0)\mathfrak{d}(2,1;\gamma;0)\oplus\mathfrak{d}(2,1;\gamma;0) is not a subalgebra of the 𝔬𝔰𝔭(8|4)\mathfrak{osp}(8^{*}|4) superconformal Lie superalgebra Estes:2018tnu . As a result, we expect the 6d ambient theory to be deformed as we tune γ\gamma away from the special value γ=1/2\gamma=-1/2. In the supergravity solutions we construct below, this deformation will appear at leading order in the asymptotic expansion of the four-form flux. Additionally, the warp factors of the symmetric spaces corresponding to T(1)T^{(1)} and T(2)T^{(2)} have the correct normalization to create an AdS7 space only when |β1|=|β2||\beta_{1}|=|\beta_{2}|, which can happen only when γ=1/2\gamma=-1/2 or γ±\gamma\to\pm\infty.

Finally, we note that the complex Lie superalgebra 𝔡(2,1;γ)\mathfrak{d}(2,1;\gamma) enjoys a triality symmetry generated by γγ1\gamma\mapsto\gamma^{-1}, γ(γ+1)\gamma\mapsto-(\gamma+1), and γγ/(γ+1)\gamma\mapsto-\gamma/(\gamma+1), any two of which are linearly independent. However, only the γγ1\gamma\mapsto\gamma^{-1} involutility survives at the level of the real form 𝔡(2,1;γ;0)\mathfrak{d}(2,1;\gamma;0), due to the distinguished nature of the T(1)T^{(1)} generator. Therefore, our analysis below can be equivalently carried out in the γ±\gamma\to\pm\infty limits, or even in the γ0\gamma\to 0 limit (albeit with a permutation of the two 𝔰𝔬(3)𝔰𝔬(3){\mathfrak{so}}(3)\oplus{\mathfrak{so}}(3) factors contained in two copies of 𝔡(2,1;γ;0){\mathfrak{d}}(2,1;\gamma;0))666In spite of the involution relating the γ\gamma\to-\infty and γ0\gamma\to 0 limits, these two descriptions cannot be smoothly deformed into one another by varying γ\gamma. Indeed, the two regimes are separated by the decompactification point γ=1\gamma=-1.. The physical interpretation of the choice of either limit corresponds to fixing the orientation of the M5-branes that engineer the ambient 6d SCFT.

3.1 Supergravity solutions for general 𝜸\gamma

In this section, we will review the structure of supergravity solutions with 𝔡(2,1;γ;0)𝔡(2,1;γ;0){\mathfrak{d}}(2,1;\gamma;0)\oplus{\mathfrak{d}}(2,1;\gamma;0) superisometry algebra, for generic γ\gamma, as first described in D_Hoker_2008b ; Estes:2012vm ; Bachas:2013vza . The odd subspace of this superalgebra is 16-dimensional for all γ\gamma, so that all supergravity backgrounds discussed below are 1/21/2-BPS. Furthermore, the maximal bosonic subalgebra of 𝔡(2,1;γ;0)𝔡(2,1;γ;0){\mathfrak{d}}(2,1;\gamma;0)\oplus{\mathfrak{d}}(2,1;\gamma;0) is

𝔰𝔬(2,2)𝔰𝔬(4)𝔰𝔬(4).\mathfrak{so}(2,2)\oplus\mathfrak{so}(4)\oplus\mathfrak{so}(4)~{}. (14)

To realize this superisometry, the supergravity metric must be of the form777For the previously mentioned special values γ=1\gamma=-1 and γ=0\gamma=0, AdS3 Wigner-İnönü contracts to 2,1\mathbb{R}^{2,1} and one of the 3-spheres decompactifies into 3\mathbb{R}^{3}, respectively.

(AdS3×𝕊3×𝕊~3)Σ2,(\text{AdS}_{3}\times\mathds{S}^{3}\times\tilde{\mathds{S}}^{3})\ltimes\Sigma_{2}~{}, (15)

for a Riemann surface Σ2\Sigma_{2}.

The Killing spinor equations on such manifolds can be recast into conditions on two auxiliary functions, hh and GG, defined over the Riemann surface Σ2\Sigma_{2} D_Hoker_2008b ; Estes:2012vm ; Bachas:2013vza . Specifically, if we employ local complex888The SUGRA orientation and Riemannian metric automatically endow the Riemann surface Σ2\Sigma_{2} with a complex structure. coordinates {w,w¯}\{w,\bar{w}\} on Σ2\Sigma_{2}, so that the metric on the Riemann surface is given by dsΣ22=dwdw¯\text{d}s^{2}_{\Sigma_{2}}=\text{d}w\text{d}\bar{w}, the function hh is constrained to be \mathbb{R}-valued and harmonic

ww¯h=0,\partial_{w}\partial_{\bar{w}}h=0~{}, (16)

while GG is \mathbb{C}-valued and satisfies the conformally covariant equation

wG=12(G+G¯)wlnh.\partial_{w}G=\frac{1}{2}\left(G+\bar{G}\right)\partial_{w}\ln h~{}. (17)

Note that both constraints above hold for generic γ\gamma.

Regularity conditions on the auxiliary functions hh and GG have been discussed in Bachas:2013vza . Here, we will consider Riemann surfaces with a boundary, Σ2\partial\Sigma_{2}\neq\emptyset. In short, regularity of the supergravity solution constrains G=±iG=\pm i on Σ2\partial\Sigma_{2}, while the (holomorphic part of the) function hh must either vanish or have a simple pole at any point of Σ2\partial\Sigma_{2}.

The triple (h,G,γ)(h,G,\gamma) uniquely specifies a bosonic background within the space of 11d supergravity solutions with superisometry algebra 𝔡(2,1;γ;0)𝔡(2,1;γ;0){\mathfrak{d}}(2,1;\gamma;0)\oplus{\mathfrak{d}}(2,1;\gamma;0). In particular, the metric field in the supergravity solution can be written in terms of this data as

ds2=fAdS32dsAdS32+f𝕊32ds𝕊32+f𝕊~32ds𝕊~32+fΣ22dsΣ22,\text{d}s^{2}=f_{\text{AdS}_{3}}^{2}\text{d}s^{2}_{\text{AdS}_{3}}+f_{\mathds{S}^{3}}^{2}\text{d}s^{2}_{\mathds{S}^{3}}+f_{\tilde{\mathds{S}}^{3}}^{2}\text{d}s^{2}_{\tilde{\mathds{S}}^{3}}+f_{\Sigma_{2}}^{2}\text{d}s^{2}_{\Sigma_{2}}~{}, (18)

where the metric factors are functions over Σ2\Sigma_{2} and can be expressed in terms of the auxiliary functions

W±:=|G±i|2+γ±1(GG¯1)\displaystyle W_{\pm}:=|G\pm i|^{2}+\gamma^{\pm 1}(G\bar{G}-1) (19)

as Bachas:2013vza

fAdS36\displaystyle f_{\text{AdS}_{3}}^{6} =h2W+Wβ16(GG¯1)2,\displaystyle=\frac{h^{2}W_{+}W_{-}}{\beta_{1}^{6}(G\bar{G}-1)^{2}}~{}, (20a)
f𝕊36\displaystyle f_{{\mathds{S}}^{3}}^{6} =h2(GG¯1)Wβ23β33W+2,\displaystyle=\frac{h^{2}(G\bar{G}-1)W_{-}}{\beta_{2}^{3}\beta_{3}^{3}W_{+}^{2}}~{}, (20b)
f𝕊~36\displaystyle f_{\tilde{\mathds{S}}^{3}}^{6} =h2(GG¯1)W+β23β33W2,\displaystyle=\frac{h^{2}(G\bar{G}-1)W_{+}}{\beta_{2}^{3}\beta_{3}^{3}W_{-}^{2}}~{}, (20c)
fΣ26\displaystyle f_{\Sigma_{2}}^{6} =|wh|6β23β33h4(GG¯1)W+W.\displaystyle=\frac{|\partial_{w}h|^{6}}{\beta_{2}^{3}\beta_{3}^{3}h^{4}}(G\bar{G}-1)W_{+}W_{-}~{}. (20d)

The bosonic sector of the 11d supergravity solution is completed by a three-form gauge potential 𝒞(3)\mathcal{C}_{(3)}, which can be written in terms of the (h,G,γ)(h,G,\gamma) data as

𝒞(3)=i=13bivoli\displaystyle\mathcal{C}_{(3)}=\sum_{i=1}^{3}b_{\mathcal{M}_{i}}\operatorname{vol}_{\mathcal{M}_{i}} (21)

where voli\operatorname{vol}_{\mathcal{M}_{i}} denotes the volume form on the manifold i={AdS3,𝕊3,𝕊~3}i\mathcal{M}_{i}=\{\text{AdS}_{3},\mathds{S}^{3},\tilde{\mathds{S}}^{3}\}_{i}. We also introduced the gauge potentials

bAdS3\displaystyle b_{\text{AdS}_{3}} :=τ1β13[h(G+G¯)1GG¯+(2+γ+γ1)Φ(γγ1)h~+b10],\displaystyle:=\frac{\tau_{1}}{\beta_{1}^{3}}\left[-\frac{h(G+\bar{G})}{1-G\bar{G}}+(2+\gamma+\gamma^{-1})\Phi-(\gamma-\gamma^{-1})\tilde{h}+b_{1}^{0}\right]~{}, (22a)
b𝕊3\displaystyle b_{\mathds{S}^{3}} :=τ2β23[γh(G+G¯)W++γ(Φh~)+b20],\displaystyle:=\frac{\tau_{2}}{\beta_{2}^{3}}\left[-\frac{\gamma h(G+\bar{G})}{W_{+}}+\gamma(\Phi-\tilde{h})+b_{2}^{0}\right]~{}, (22b)
b𝕊~3\displaystyle b_{\tilde{\mathds{S}}^{3}} :=τ3β33[h(G+G¯)γWΦ+h~γ+b30],\displaystyle:=\frac{\tau_{3}}{\beta_{3}^{3}}\left[\frac{h(G+\bar{G})}{\gamma W_{-}}-\frac{\Phi+\tilde{h}}{\gamma}+b_{3}^{0}\right]~{}, (22c)

where {bi0}i=13\{b^{0}_{i}\}_{i=1}^{3} are integration constants, while the dual harmonic function h~\tilde{h} and the real auxiliary function Φ\Phi are defined in terms of hh and GG via

wh~\displaystyle\partial_{w}\tilde{h} =iwh,\displaystyle=-i\partial_{w}h~{}, (23a)
wΦ\displaystyle\partial_{w}\Phi =G¯wh.\displaystyle=\bar{G}\partial_{w}h~{}. (23b)

Finally, in eq. (22) we made use of the signs τi=±1\tau_{i}=\pm 1 for i{1,2,3}i\in\{1,2,3\}, which are subject to the constraint

i=13βifi+hi=13τi=0.\prod_{i=1}^{3}\beta_{i}f_{\mathcal{M}_{i}}+h\prod_{i=1}^{3}\tau_{i}=0~{}. (24)

As discussed above, regularity conditions force hh and GG to have prescribed behavior on Σ2\partial\Sigma_{2}. The solutions that we are interested in studying are asymptotically AdS×7𝕊4{}_{7}\times\mathds{S}^{4}, which are indeed described by G=±iG=\pm i and hh having a single simple pole, which corresponds to having a single asymptotic AdS×7𝕊4{}_{7}\times\mathds{S}^{4} region. Generically at any point on Σ2\partial\Sigma_{2}, the regularity constraint on hh implies that its Laurent expansion reduces to

h=ih0w+c.c.=2h0sinϑϱ,h=-ih_{0}w+\operatorname{c.c.}=\frac{2h_{0}\sin\vartheta}{\varrho}~{}, (25)

for some real constant h0h_{0}. In the second equality, we have introduced new convenient set of coordinates {ϱ,ϑ}\{\varrho,\vartheta\} on Σ2\Sigma_{2} defined by w=eiϑ/ϱw=e^{i\vartheta}/\varrho. Adopting these new polar coordinates and using eq. (25) allows us to perturbatively solve eq. (17) in small ϱ\varrho to find

G=i+a1ϱeiϑsinϑ+𝒪(ϱ2)G=-i+a_{1}\varrho e^{i\vartheta}\sin\vartheta+\mathcal{O}(\varrho^{2}) (26)

for some constant a1a_{1}. If we then insert eq. (26) into eqs. (19) and (27), we can perturbatively expand eq. (18) in small ϱ\varrho to find

ds2=L2dϱ2ϱ22γL2a1(1+γ)2ϱ(dsAdS32+(1+γ)2γ2ds𝕊32)+L2(dϑ2+sin2ϑds𝕊~32)+\text{d}s^{2}=L^{2}\frac{\text{d}\varrho^{2}}{\varrho^{2}}-\frac{2\gamma L^{2}}{a_{1}(1+\gamma)^{2}\varrho}\left(\text{d}s^{2}_{\text{AdS}_{3}}+\frac{(1+\gamma)^{2}}{\gamma^{2}}\text{d}s^{2}_{\mathds{S}^{3}}\right)+L^{2}\left(\text{d}\vartheta^{2}+\sin^{2}\vartheta\text{d}s^{2}_{\tilde{\mathds{S}}^{3}}\right)+\dots (27)

with

L6=a12h02(1+γ)6β16γ2.L^{6}=\frac{a_{1}^{2}h_{0}^{2}(1+\gamma)^{6}}{\beta_{1}^{6}\gamma^{2}}~{}. (28)

We can see the AdS×7𝕊4{}_{7}\times\mathds{S}^{4} asymptotic geometry in eq. (27) clearly for certain values of γ\gamma. The obvious case is if we choose γ=1/2\gamma=-1/2, which was studied extensively in Bachas:2013vza ; Gentle:2015jma ; Estes:2018tnu . The other possibility, namely the limits γ±\gamma\to\pm\infty, will be discussed in the next section.

For γ<0\gamma<0 and for the function hh given in eq. (25), the general solution to eq. (17) for the function GG with an even number 2n+22n+2 of branch points {ξi}i=12n+2Σ2\{\xi_{i}\}_{i=1}^{2n+2}\subset\partial\Sigma_{2} was found in DHoker:2008rje to be

G=i(1+j=12n+2(1)jwξj|wξj|),G=-i\left(1+\sum_{j=1}^{2n+2}(-1)^{j}\frac{w-\xi_{j}}{|w-\xi_{j}|}\right)~{}, (29)

where GG flips sign ±ii\pm i\to\mp i upon crossing each branch point ξi\xi_{i}.

3.2 The 𝜸\gamma\to-\infty limit

Of particular interest to us in the following sections are solutions that are constructed by taking the scaling limit

γ±,\displaystyle\gamma\to\pm\infty~{}, γa1constant,\displaystyle\gamma a_{1}\to\text{constant}~{}, Lconstant.\displaystyle L\to\text{constant}~{}. (30)

From eq. (27) and the requirement γa1constant\gamma a_{1}\to\text{constant}, we can see that this scaling limit realizes an AdS×7𝕊4{}_{7}\times\mathds{S}^{4} asymptotic geometry. In this and the following subsections, we will consider the γ\gamma\to-\infty limit in greater detail and show that the solutions engineered in this limit contain the class of solutions found in Faedo:2020nol which we reviewed in section 2. Below, we will expand upon these solutions and construct new surface defects with small 𝒩=(4,4)\mathcal{N}=(4,4) SUSY.

To begin, we introduce a complex function FF via

G=i(1+γ1F),G=-i\left(1+\gamma^{-1}F\right)~{}, (31)

and we rescale β1=β^(γ)1/3\beta_{1}=\hat{\beta}(-\gamma)^{1/3}. This ensures that LL remains finite as required by the limiting procedure of eq. (30). In this limit, the metric in eq. (27) becomes

ds2\displaystyle\text{d}s^{2} =[4h2β^6(F+F¯)]13(dsAdS32+ds𝕊32)+[h2(F+F¯)216β^6]13(ds𝕊~32+4|wh|2h2dwdw¯).\displaystyle=\left[\frac{4h^{2}}{\hat{\beta}^{6}(F+\bar{F})}\right]^{\frac{1}{3}}\left(\text{d}s^{2}_{\text{AdS}_{3}}+\text{d}s^{2}_{\mathds{S}^{3}}\right)+\left[\frac{h^{2}(F+\bar{F})^{2}}{16\hat{\beta}^{6}}\right]^{\frac{1}{3}}\left(\text{d}s^{2}_{\tilde{\mathds{S}}^{3}}+\frac{4|\partial_{w}h|^{2}}{h^{2}}\text{d}w\text{d}\bar{w}\right)~{}. (32)

In addition, the gauge potentials are given by

bAdS3=b𝕊3=2h~β^3,b𝕊~3=h2β^3i(FF¯)2Φ^,\displaystyle b_{\text{AdS}_{3}}=b_{\mathds{S}^{3}}=-\frac{2\tilde{h}}{\hat{\beta}^{3}},\qquad b_{\tilde{\mathds{S}}^{3}}=\frac{h}{2\hat{\beta}^{3}}\frac{-i(F-\bar{F})}{2}-\hat{\Phi}~{}, (33)

where the function Φ^\hat{\Phi} is defined via

wΦ^\displaystyle\partial_{w}\hat{\Phi} =F¯wh~β^3.\displaystyle=\bar{F}\frac{\partial_{w}\tilde{h}}{\hat{\beta}^{3}}~{}. (34)

Let us introduce a new set of coordinates {z,ρ}\{z,\rho\} on Σ2\Sigma_{2} defined via w=z+iρw=z+i\rho, and choose h0h_{0} in eq. (25) such that the harmonic function h=β^3h1ρh=\hat{\beta}^{3}h_{1}\rho for some constant h1h_{1}. Let us also parametrize the real and imaginary parts of the complex function FF such that

F=2ρ2h1H+iFI,F=\frac{2\rho^{2}}{h_{1}}H+iF_{I}~{}, (35)

for real functions HH and FIF_{I}. Note that the latter does not enter the metric field; indeed, we can rewrite the line element in eq. (32) entirely in terms of HH as

ds2=h1H13(dsAdS32+ds𝕊32)+H23(dz2+dρ2+ρ2ds𝕊~32).\text{d}s^{2}=h_{1}H^{-\frac{1}{3}}\left(\text{d}s^{2}_{\text{AdS}_{3}}+\text{d}s^{2}_{\mathds{S}^{3}}\right)+H^{\frac{2}{3}}\left(\text{d}z^{2}+\text{d}\rho^{2}+\rho^{2}\text{d}s^{2}_{\tilde{\mathds{S}}^{3}}\right)~{}. (36)

The condition in eq. (17) implies that the complex function FF satisfies the first order equation

wF=12(FF¯)wlnh,\partial_{w}F=\frac{1}{2}(F-\bar{F})\partial_{w}\ln h~{}, (37)

which in turn can be recast into a second order equation for its real part,

ρ2H+3ρρH+z2H=0.\partial_{\rho}^{2}H+\frac{3}{\rho}\partial_{\rho}H+\partial_{z}^{2}H=0~{}. (38)

The supergravity solution in the {z,ρ}\{z,\rho\} parametrization of Σ2\Sigma_{2} is completed by the following expressions for the gauge potentials,

bAdS3=b𝕊3=2h1z,b𝕊~3=h1ρ2FI+Φ^.\displaystyle b_{\text{AdS}_{3}}=b_{\mathds{S}^{3}}=2h_{1}z,\qquad b_{\tilde{\mathds{S}}^{3}}=\frac{h_{1}\rho}{2}F_{I}+\hat{\Phi}~{}. (39)

In particular, a more explicit form can be given for the derivatives of b𝕊~3b_{\tilde{\mathds{S}}^{3}} as

ρb𝕊~3=ρ3zH,zb𝕊~3=ρ3ρH.\displaystyle\partial_{\rho}b_{\tilde{\mathds{S}}^{3}}=\rho^{3}\partial_{z}H,\qquad\partial_{z}b_{\tilde{\mathds{S}}^{3}}=-\rho^{3}\partial_{\rho}H~{}. (40)

These solutions have a scaling symmetry under the transformation

zλz,\displaystyle z\rightarrow\lambda z, ρλρ,\displaystyle\rho\rightarrow\lambda\rho, Hλ3H,\displaystyle H\rightarrow\lambda^{-3}H, h1λ1h1,\displaystyle h_{1}\rightarrow\lambda^{-1}h_{1}, (41)

for which the metric and flux are invariant. This transformation could be used to fix the value of h1h_{1} so that the solution is uniquely given by choice of function HH.

3.3 Inclusion of KK-monopoles and recovering known solutions

At the level of the local geometry described by eq. (36), KK-monopoles can be included in a straightforward fashion by replacing 𝕊3\mathds{S}^{3} and 𝕊~3\tilde{\mathds{S}}^{3} with the lens spaces 𝕊3/k\mathds{S}^{3}/\mathbb{Z}_{k} and 𝕊~3/k\tilde{\mathds{S}}^{3}/\mathbb{Z}_{k^{\prime}}, respectively, where kk and kk^{\prime} are the orbifold charges. The lens spaces can be realized as the total spaces of circle bundles over 2-spheres, where the orbifold acts on the Hopf fiber; this amounts to the substitutions

ds𝕊32\displaystyle\text{d}s^{2}_{\mathds{S}^{3}}\to ds𝕊3/k2=14[(dχk+ω)2+ds𝕊22],\displaystyle\ \text{d}s^{2}_{\mathds{S}^{3}/\mathbb{Z}_{k}}\ =\frac{1}{4}\left[\left(\frac{\text{d}\chi}{k}+\omega\right)^{2}+\text{d}s^{2}_{\mathds{S}^{2}}\right]~{}, (42a)
ds𝕊~32\displaystyle\text{d}s^{2}_{\tilde{\mathds{S}}^{3}}\to ds𝕊~3/k2=14[(dϕk+η)2+ds𝕊~22],\displaystyle\ \text{d}s^{2}_{\tilde{\mathds{S}}^{3}/\mathbb{Z}_{k^{\prime}}}=\frac{1}{4}\left[\left(\frac{\text{d}\phi}{k^{\prime}}+\eta\right)^{2}+\text{d}s^{2}_{\tilde{\mathds{S}}^{2}}\right]~{}, (42b)

where dω=vol𝕊2\text{d}\omega=\operatorname{vol}_{\mathds{S}^{2}} and dη=vol𝕊~2\text{d}\eta=\operatorname{vol}_{\tilde{\mathds{S}}^{2}}. The inclusion of either KK-monopole incurs the loss of 1/21/2 of the existing supersymmetry generators, resulting in a small 𝒩=(0,4)\mathcal{N}=(0,4) supersymmetry algebra. Indeed, dimensional reduction along the Taub-NUT direction produces a D6-brane which breaks half of the supersymmetries of the massless type IIA AdS7 vacuum. The second KK-monopole can be added without bringing about any further breaking of the supersymmetries, and provides a second isometric direction upon which the background can be dimensionally reduced to massless type IIA supergravity. Therefore, the final solution is an 1/81/8-BPS configuration. At the level of the superconformal symmetry algebra, this corresponds to a reduction from 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}\left(2,1;\gamma\right)\oplus\mathfrak{d}\left(2,1;\gamma\right) to 𝔡(2,1;γ)\mathfrak{d}\left(2,1;\gamma\right).

The inclusion of orbifolds allows us to match the γ\gamma\to-\infty solutions of section 3.2 to those found in Faedo:2020nol and reviewed in section 2. In particular, we see that the γ\gamma\to-\infty metric in eq. (36), subject to the inclusion of KK-monopoles as in eq. (42), matches the M2-M5-KK-M5-KK near-horizon in eq. (2) under the identifications h1=4kQM5h_{1}=4kQ_{\text{M5}} and H=HM5H=H_{\text{M5}^{\prime}}. In particular, the functions HH and HM5H_{\text{M5}^{\prime}} solve the same equation, since eq. (38) maps to eq. (3) under the rescaling ρ^=ρ2/4\hat{\rho}=\rho^{2}/4. We have thus managed to fully recover the solutions described in Faedo:2020nol as a limiting case of those in Bachas:2013vza .

4 New small 𝓝=(𝟒,𝟒)\mathcal{N}=(4,4) surface defects

In this section, we will construct a new explicit family of solutions H(z,ρ)H(z,\rho) to eq. (38). In particular, without loss of generality we will specialize to the γ\gamma\to-\infty limit. We recall that for γ<0\gamma<0, a general solution with an even number of branch points {ξi}i=12n+2Σ2\{\xi_{i}\}_{i=1}^{2n+2}\subset\partial\Sigma_{2} was given by eq. (29). We rescale the singular points ξj\xi_{j} as

ξj=νjγ1ξ^jΣ2forj{1,2,,2n+2},\xi_{j}=\nu_{j}-\gamma^{-1}\hat{\xi}_{j}\in\partial\Sigma_{2}\quad\text{for}\quad j\in\{1,2,\dots,2n+2\}~{}, (43)

where the collapse points satisfy νjνj+1\nu_{j}\leq\nu_{j+1} for all jj, so as to preserve the total order of the set {ξj}\{\xi_{j}\} following the γ\gamma\to-\infty limit. Furthermore, to ensure finiteness of the complex function FF in this limit, we must also demand that all points νj\nu_{j}, for 1j2n+21\leq j\leq 2n+2, correspond to the collapse of even-dimensional clusters {ξk,ξk+1,,ξk+2m+1}\{\xi_{k},\xi_{k+1},\ldots,\xi_{k+2m+1}\} of singular points, where kjk+2m+1k\leq j\leq k+2m+1. This can be implemented by identifying

νkνk+1νk+2m+1\nu_{k}\equiv\nu_{k+1}\equiv\ldots\equiv\nu_{k+2m+1} (44)

for each cluster. Finally, in order to preserve the ordering of the singular points throughout the γ\gamma\to-\infty limit, we must also arrange the corresponding collapse parameters within each cluster so that ξ^k<ξ^k+1<<ξ^k+2m+1\hat{\xi}_{k}<\hat{\xi}_{k+1}<\ldots<\hat{\xi}_{k+2m+1}. Without loss of generality, then, it suffices to consider a pairwise collapse of neighboring singular points, which can be realized by identifying

ν2i1ν2ifori{1,2,,n+1}.\nu_{2i-1}\equiv\nu_{2i}\quad\text{for}\quad i\in\{1,2,\dots,n+1\}~{}. (45)

However, we will later comment on certain phenomena which appear only when the singular points collapse in clusters of 4 or more.

Σ2\Sigma_{2}zzξ1\xi_{1}ν1\nu_{1}ξ2\xi_{2}ξ3\xi_{3}ν3\nu_{3}ξ4\xi_{4}\cdotsξ2n+1\xi_{2n+1}ν2n+1\nu_{2n+1}ξ2n+2\xi_{2n+2}γ\gamma\rightarrow-\infty
Figure 1: A particular choice of collapse points νjΣ2\nu_{j}\in\partial\Sigma_{2} and the behavior of the singular points ξiΣ2\xi_{i}\in\partial\Sigma_{2} as γ\gamma\to-\infty under pairwise collapse. The collapse points νj\nu_{j} can be chosen arbitrarily, while maintaining νj<νj+1\nu_{j}<\nu_{j+1}, but for clarity in the figure have been depicted at the midpoint in the region [ξj,ξj+1][\xi_{j},\,\xi_{j+1}] to demonstrate pairwise collapse.

An illustration of the pairwise collapse described by eq. (43) and eq. (45) is shown in figure 1. Under this collapse dynamic, the complex function FF becomes, in the γ\gamma\to-\infty limit,

F(w,w¯)=j=12n+2(1)jξ^jw¯w2(w¯νj)|wνj|.F(w,\bar{w})=\sum_{j=1}^{2n+2}(-1)^{j}\hat{\xi}_{j}\frac{\bar{w}-w}{2(\bar{w}-\nu_{j})|w-\nu_{j}|}~{}. (46)

The corresponding function HH is given by

H(z,ρ)=h12j=12n+2(1)jξ^j(ρ2+(zνj)2)3/2.H(z,\rho)=\frac{h_{1}}{2}\sum_{j=1}^{2n+2}\frac{(-1)^{j}\hat{\xi}_{j}}{(\rho^{2}+(z-\nu_{j})^{2})^{3/2}}~{}. (47)

For pure AdS×7𝕊4{}_{7}\times\mathds{S}^{4}, this procedure produces what we refer to as the single-pole vacuum (1PV) solution with

F1PV(w,w¯)\displaystyle F_{\text{1PV}}(w,\bar{w}) =ξ^w¯ww¯|w|,\displaystyle=\hat{\xi}\frac{\bar{w}-w}{\bar{w}|w|}~{}, (48a)
H1PV(z,ρ)\displaystyle H_{\text{1PV}}(z,\rho) =h1ξ^(z2+ρ2)32,\displaystyle=\frac{h_{1}\hat{\xi}}{(z^{2}+\rho^{2})^{\frac{3}{2}}}~{}, (48b)

which can be obtained by taking n=0n=0, ξ2=ξ1=ξ\xi_{2}=-\xi_{1}=\xi, and ν1=ν2=0\nu_{1}=\nu_{2}=0. Alternatively, this solution can also be obtained by collapsing all singular points to the origin, i.e. νj=0\nu_{j}=0 for all jj. Note that eq. (38) is linear in HH and admits a translation symmetry under shifts of zz. Using these two properties, the general solutions given in eq. (47) can be reconstructed from the 1PV solution by taking linear combinations and making use of the translation symmetry.

One may worry that since the solution for the potential HH in eq. (47) is generically singular at the points (z,ρ)=(νi,0)(z,\rho)=(\nu_{i},0) for all i{1,2,,2n+2}i\in\{1,2,\dots,2n+2\}, the resulting supergravity metric may have a singularity or these points may simply correspond to horizons. In order to investigate the regularity of the spacetime geometry, we compute the scalar curvature at these points. Without loss of generality we can choose to evaluate the Ricci scalar at the point (z,ρ)=(ν2j,0)(z,\rho)=(\nu_{2j},0). Recall that for metrics generally of the form in eq. (36), the scalar curvature is given by eq. (7) with HM5H_{\rm M5^{\prime}} replaced by HH. Note that since all collapse points ν2k1\nu_{2k-1} with odd labels are identified with even-labelled ones ν2k\nu_{2k} via eq. (45), the point (z,ρ)=(ν2j,0)(z,\rho)=(\nu_{2j},0) is indeed generic within the set of collapsed points. Ultimately, we find

|(z,ρ)=(ν2j,0)=32L𝕊42(ξ^2jξ^2j1m^1)2/3,\mathcal{R}\bigr{|}_{(z,\rho)=(\nu_{2j},0)}=\frac{3}{2L_{\mathds{S}^{4}}^{2}}\left(\frac{\hat{\xi}_{2j}-\hat{\xi}_{2j-1}}{\hat{m}_{1}}\right)^{-2/3}~{}, (49)

which is indeed finite, where L𝕊4L_{\mathbb{S}^{4}} and m^1\hat{m}_{1} are strictly positive constants introduced below. This suggests that the geometry is regular at these points.

We also note that the 1PV solution corresponds to a spacetime with a constant scalar curvature given by

=32L𝕊42,\mathcal{R}=\frac{3}{2L_{\mathds{S}^{4}}^{2}}~{}, (50)

with L𝕊4=(h1ξ^)1/3L_{\mathds{S}^{4}}=(h_{1}\hat{\xi})^{1/3}.

4.1 Asymptotic local behavior

We will now show explicitly that the solutions described by eq. (47) are asymptotically locally AdS×7𝕊4{}_{7}\times\mathds{S}^{4}. As described in greater detail in appendix A, the Riemann surface Σ2\Sigma_{2} admits a parametrization by Fefferman-Graham (FG) coordinates {v,ϕ}\{v,\phi\} in terms of which the line element takes the following asymptotic form for small vv,

ds2=4L𝕊42v2[dv2+α1(dsAdS32+ds𝕊32)]+L𝕊42[α3dϕ2+α4sin2ϕds𝕊~32],\text{d}s^{2}=\frac{4L_{\mathds{S}^{4}}^{2}}{v^{2}}\left[\text{d}v^{2}+\alpha_{1}\left(\text{d}s^{2}_{\text{AdS}_{3}}+\text{d}s^{2}_{\mathds{S}^{3}}\right)\right]+L_{\mathds{S}^{4}}^{2}\left[\alpha_{3}\text{d}\phi^{2}+\alpha_{4}\sin^{2}\phi\ \text{d}s^{2}_{\tilde{\mathds{S}}^{3}}\right], (51)

where the metric factors are of the form αi=1+𝒪(v4)\alpha_{i}=1+\mathcal{O}(v^{4}) for i{1,3,4}i\in\{1,3,4\}. They are given explicitly, together with the asymptotic mapping to FG coordinates on Σ2\Sigma_{2}, in appendix A. The asymptotic 𝕊4\mathds{S}^{4} radius L𝕊4L_{\mathds{S}^{4}} can be expressed in terms of the moments

m^k:=j=12n+2(1)jξ^jk\hat{m}_{k}:=\sum_{j=1}^{2n+2}(-1)^{j}\hat{\xi}_{j}^{k} (52)

as

L𝕊43=h1m^12.L_{\mathds{S}^{4}}^{3}=\frac{h_{1}\hat{m}_{1}}{2}~{}. (53)

As claimed above, we may recognize within eq. (51), at leading order in vv, the large-xx limit of the line element of AdS7 (with radius 2L𝕊42L_{\mathds{S}^{4}}) written in AdS3 slicing,

dsAdS72=4L𝕊42[dx2+cosh2xdsAdS32+sinh2xds𝕊32],\text{d}s^{2}_{\text{AdS}_{7}}=4L^{2}_{\mathds{S}^{4}}\left[\text{d}x^{2}+\cosh^{2}x\ \text{d}s^{2}_{\text{AdS}_{3}}+\sinh^{2}x\ \text{d}s^{2}_{\mathds{S}^{3}}\right], (54)

where the coordinate normal to the AdS3 foliation is related to the FG coordinate via x=log(2/v)x=\log(2/v). The large-xx limit trivializes the relative warping coth2x\coth^{2}x between the AdS3 and the 𝕊3\mathds{S}^{3} subspaces, matching the behavior seen in eq. (51).

In the γ\gamma\to-\infty limit, the auxiliary functions Φ\Phi and h~\tilde{h} admit the following FG expansions,

Φ=h~\displaystyle\Phi=-\tilde{h} =4m^1cosϕv2+2n^1m^1+m^1n^2n^12m^13cosϕv2+𝒪(v4),\displaystyle=\frac{4\hat{m}_{1}\cos\phi}{v^{2}}+\frac{2\hat{n}_{1}}{\hat{m}_{1}}+\frac{\hat{m}_{1}\hat{n}_{2}-\hat{n}_{1}^{2}}{\hat{m}_{1}^{3}}\cos\phi~{}v^{2}+\mathcal{O}(v^{4})~{}, (55)

where the moments n^i\hat{n}_{i} are defined in eq. (97). Using these expansions, we find that the asymptotic geometry in eq. (51) is supported by a four-form flux999In the following, we have made a choice on the signs τ1,2,3\tau_{1,2,3} in line with the constraint i=13τi=1\prod_{i=1}^{3}\tau_{i}=1 which follows from evaluating eq. (24) on the γ\gamma\to-\infty solutions.

(4)L𝕊43\displaystyle\frac{\mathcal{F}_{(4)}}{L_{\mathds{S}^{4}}^{3}} =16cosϕv3dv(vol𝕊3+volAdS3)8sinϕv2dϕ(vol𝕊3+volAdS3)\displaystyle=-\frac{16\cos\phi}{v^{3}}\text{d}v\wedge(\operatorname{vol}_{\mathds{S}^{3}}+\operatorname{vol}_{\text{AdS}_{3}})-\frac{8\sin\phi}{v^{2}}\text{d}\phi\wedge(\operatorname{vol}_{\mathds{S}^{3}}+\operatorname{vol}_{\text{AdS}_{3}}) (56)
+3sin3ϕdϕvol𝕊~3+ 4cosϕm^1n^2n^12m^14vdv(vol𝕊3+volAdS3)+𝒪(v2).\displaystyle\quad+3\sin^{3}\phi\ \text{d}\phi\wedge\operatorname{vol}_{\tilde{\mathds{S}}^{3}}+\ 4\cos\phi\ \frac{\hat{m}_{1}\hat{n}_{2}-\hat{n}_{1}^{2}}{\hat{m}_{1}^{4}}v\ \text{d}v\wedge(\operatorname{vol}_{\mathds{S}^{3}}+\operatorname{vol}_{\text{AdS}_{3}})+\mathcal{O}(v^{2}).

The first line in the equation above manifests the deformation of the 6d ambient SCFT by the insertion of a source, as already hinted at by the superalgebra structure discussed in section 3. This deformation takes the form of an S-wave over the internal 𝕊~3\tilde{\mathds{S}}^{3}. It can be compared to the undeformed theory, which corresponds to γ=1/2\gamma=-1/2. In that case, the four-form field strength at the conformal boundary v=0v=0 is (4)=3L𝕊43vol𝕊4\mathcal{F}_{(4)}=3L^{3}_{\mathds{S}^{4}}\operatorname{vol}_{\mathds{S}^{4}}, where 𝕊4\mathds{S}^{4} is the internal 4-sphere spanned by ϕ\phi and 𝕊~3\tilde{\mathds{S}}^{3}. We can recognize this term as the first term in the second line of eq. (56).

4.2 Single-pole vacuum

Σ2\Sigma_{2}zzν\nuξ1\xi_{1}ξ2\xi_{2}a\mathfrak{C}_{a}a\mathfrak{C}_{a}Σ2\Sigma_{2}ξ1\xi_{1}ν\nuξ2\xi_{2}
Figure 2: The 1PV solution described in section 4.2 is characterized by two singular points ξ{1,2}\xi_{\{1,2\}}, which collapse to the same νΣ2\nu\in\partial\Sigma_{2} in the γ\gamma\to-\infty limit. As discussed in section 4.3, the basis of non-contractible four-cycles for this solution is one-dimensional. The profile of a representative cycle a\mathfrak{C}_{a} along Σ2\Sigma_{2} is shown; note that the same 3-sphere collapses at both of its endpoints on Σ2\partial\Sigma_{2}. On the right, the 1PV on the upper half plane is mapped to a semi-infinite strip.

For later reference, we now identify a vacuum solution within the family of supergravity backgrounds presented above. An appropriate choice of vacuum is necessary for computing properties associated to a defect embedded in a holographic CFT. Indeed, in order to isolate quantities which are intrinsic to a defect, one must ensure that the contributions due to the ambient degrees of freedom are taken into account. The gravitational analogue of this operation, upon recasting a field theory computation into a bulk one, is vacuum subtraction. In particular we must use a vacuum solution which is characterized by the same bulk deformation as the general solutions discussed above. As discussed above, pure AdS×7𝕊4{}_{7}\times\mathds{S}^{4} with γ=1/2\gamma=-1/2 does not satisfy this criteria, as the solutions we consider here with γ\gamma\rightarrow-\infty contain a bulk deformation as can be seen in the asymptotic expression for the flux given in eq. (56).

We take the vacuum to be the 1PV we identified in eq. (48), which is the solution corresponding to the γ\gamma\rightarrow-\infty of pure AdS×7𝕊4{}_{7}\times\mathds{S}^{4}, i.e. two singular points ξ{1,2}\xi_{\{1,2\}} collapsing to a single point ν\nu, as shown in figure 2. More precisely, for general γ\gamma, the 1PV metric is given in FG form by

ds1PV2(γ)L𝕊42=4v2\displaystyle\allowdisplaybreaks\frac{\text{d}s^{2}_{\text{1PV}}(\gamma)}{L_{\mathds{S}^{4}}^{2}}=\frac{4}{v^{2}} [dv2+(1+2γ+3(2γ+1)c2ϕ16(γ+1)2v2)dsAdS32\displaystyle\left[\text{d}v^{2}+\left(1+\frac{2\gamma+3-(2\gamma+1)c_{2\phi}}{16(\gamma+1)^{2}}v^{2}\right)\text{d}s^{2}_{\text{AdS}_{3}}\right. (57)
+((γ+1)2γ2+2γ1(2γ+1)c2ϕ16γ2v2)ds𝕊32+𝒪(v4)]\displaystyle\ +\left.\left(\frac{(\gamma+1)^{2}}{\gamma^{2}}+\frac{2\gamma-1-(2\gamma+1)c_{2\phi}}{16\gamma^{2}}v^{2}\right)\text{d}s^{2}_{\mathds{S}^{3}}+\mathcal{O}(v^{4})\right]
+\displaystyle+ [(1+(2γ+1)(2c2ϕ+1)12(γ+1)2v2)sϕ2ds𝕊~32+(1+(2γ+1)cϕ24(γ+1)2v2)dϕ2+𝒪(v4)],\displaystyle\left[\left(1+\frac{(2\gamma+1)(2c_{2\phi}+1)}{12(\gamma+1)^{2}}v^{2}\right)s^{2}_{\phi}\text{d}s_{\tilde{\mathds{S}}^{3}}^{2}+\left(1+\frac{(2\gamma+1)c^{2}_{\phi}}{4(\gamma+1)^{2}}v^{2}\right)\text{d}\phi^{2}+\mathcal{O}(v^{4})\right],

where in order to keep the expression manageable, we have adopted the abusive notation

cosxcx,andsinxsx,\displaystyle\cos x\equiv c_{x},\quad\text{and}\quad\sin x\equiv s_{x}, (58)

which will be employed from this point forward. For γ=1/2\gamma=-1/2, the 1PV recovers the AdS×7𝕊4{}_{7}\times\mathds{S}^{4} vacuum in FG gauge. In the γ\gamma\to-\infty limit which is of relevance here, the line element of the 1PV becomes instead

ds1PV2(γ)=4L𝕊4v2[dv2+dsAdS32+ds𝕊32+𝒪(v4)]+L𝕊42[sϕ2ds𝕊~32+dϕ2+𝒪(v4)].\text{d}s^{2}_{\text{1PV}}(\gamma\to-\infty)=\frac{4L_{\mathds{S}^{4}}}{v^{2}}\left[\text{d}v^{2}+\text{d}s^{2}_{\text{AdS}_{3}}+\text{d}s^{2}_{\mathds{S}^{3}}+\mathcal{O}(v^{4})\right]+L_{\mathds{S}^{4}}^{2}\left[s^{2}_{\phi}\text{d}s_{\tilde{\mathds{S}}^{3}}^{2}+\text{d}\phi^{2}+\mathcal{O}(v^{4})\right]. (59)

Furthermore, as a quick sanity check, we can take the 1PV limit of eq. (49), which recovers |1PV=3/(2L𝕊42)\mathcal{R}|_{\rm 1PV}=3/(2L_{\mathds{S}^{4}}^{2}) as expected.

As remarked above, the 1PV contains no explicit defect data – as we will see later, no Young Tableau can be associated to it. However, it does fully capture the bulk S-wave deformation discussed previously. This follows from the fact that all terms in eq. (56) which do not vanish at the conformal boundary v=0v=0 are independent of the moments {m^i,n^j}\{\hat{m}_{i},\hat{n}_{j}\}. Therefore, the 6d ambient theory dual to the 1PV is deformed by the same sources as the theories dual to completely generic γ\gamma\to-\infty bulk solutions, and so, in this limit, there is no smooth deformation of the field theory parameters that restores the ambient conformal symmetry in full. This is to be contrasted with the global AdS×7𝕊4{}_{7}\times\mathds{S}^{4} solution, which enjoys the full SO(6,2)SO(6,2) conformal symmetry. Lastly, in taking the 1PV limit, the solution exhibits a flavor symmetry enhancement SO(4)SO(5)SO(4)\to SO(5).

In the construction of Faedo:2020nol reviewed in section 2, the 1PV solution can be obtained from eq. (3) by taking α0\alpha\rightarrow 0 and setting g3=22/h1ξ^g^{3}=2\sqrt{2}/h_{1}\hat{\xi}. Alternatively, it can also be obtained in the limit g0g\rightarrow 0. To see this, first make a scale transformation using the scaling symmetry given by eq. (41), take gg to scale with λ\lambda as g3=22/h1ξ^λ3g^{3}=2\sqrt{2}/h_{1}\hat{\xi}\lambda^{3}, and then take λ\lambda\rightarrow\infty.

4.3 Partition data

In this subsection, our aim is to identify within the γ\gamma\to-\infty solutions of eq. (47) a basis of independent, non-contractible cycles threaded by four-form flux. Integrating the flux along these cycles will allow us to compute the integral M-brane charges that label a supergravity solution. In turn, this characterization will enable us to recast the specification of a supergravity solution in the form of a partition containing the representation data associated to the defect string in the dual gauge theory, in analogy with the Wilson surfaces of the γ=1/2\gamma=-1/2 solutions Estes:2018tnu . We begin by searching for independent, non-contractible four-cycles through the 11d geometry of the pairwise-collapsed γ\gamma\to-\infty solutions. Later, we will comment on how the cycles are modified if less generic collapse dynamics are considered.

It is straightforward to see that, in the γ\gamma\rightarrow-\infty limit, the volume of 𝕊3\mathds{S}^{3} vanishes on Σ2{νj}j=12n+2\partial\Sigma_{2}-\{\nu_{j}\}_{j=1}^{2n+2}, i.e. all along the boundary of the Riemann surface, except at the locations of the collapse points νj\nu_{j}. In turn, any open curve on Σ2\Sigma_{2} with end points on Σ2{νj}j=12n+2\partial\Sigma_{2}-\{\nu_{j}\}_{j=1}^{2n+2} will be a closed curved in the full 11d geometry. Therefore, any curve encircling at least one of the distinct collapse midpoints νj\nu_{j} will not be contractible to a point. This observation allows us to build a basis for non-contractible four-cycles a\mathfrak{C}_{a}, by taking the product of such curves on Σ2\Sigma_{2} with 𝕊~3\tilde{\mathds{S}}^{3}, i.e.

a{Reiθν2a1| 0θπ}×𝕊~3,\mathfrak{C}_{a}\equiv\{Re^{i\theta}-\nu_{2a-1}\,|\,0\leq\theta\leq\pi\}\times\tilde{\mathds{S}}^{3}\,, (60)

for 1an+11\leq a\leq n+1. Note that any curve enveloping multiple νa\nu_{a}’s can be decomposed into a linear combination of curves enclosing a single collapse point; hence, to build an irreducible basis of cycles, we take the radius RR above such that a\mathfrak{C}_{a} encircles only a single collapse point νa\nu_{a}.

The a\mathfrak{C}_{a} cycles alone do not exhaust the set of all non-contractible four-cycles. Indeed, we can consider a four-cycle constructed from a curve on Σ2\Sigma_{2} connecting singular points νa\nu_{a}. For irreducibility, we only consider curves connecting neighboring collapse points. While the singular nature of their endpoints might make such curves appear problematic at first glance, building a regular four-cycle from such a curve is possible as the volume of 𝕊3\mathds{S}^{3} vanishes at both endpoints, while the volume of 𝕊~3\tilde{\mathds{S}}^{3} remains finite. Indeed, the metric factor of 𝕊3\mathds{S}^{3} contains h1H1/3h_{1}H^{-1/3}, while the metric factor of 𝕊~3\tilde{\mathds{S}}^{3} has ρ2H2/3\rho^{2}H^{2/3}, and from eq. (47), the function HH goes as ρ3\rho^{-3} near any collapse point νa\nu_{a}. Thus, the cycle

a{12(ν2a+1ν2a1)eiθ+12(ν2a+1+ν2a1)| 0θπ}×𝕊3,\mathfrak{C}^{\prime}_{a}\equiv\left\{\frac{1}{2}(\nu_{2a+1}-\nu_{2a-1})e^{i\theta}+\frac{1}{2}(\nu_{2a+1}+\nu_{2a-1})\,\Big{|}\,0\leq\theta\leq\pi\right\}\times\mathds{S}^{3}~{}, (61)

has the desired behavior of a non-contractible four-cycle. The distinction between the a\mathfrak{C}_{a} and a\mathfrak{C}_{a}^{\prime} cycles is illustrated in figure 3. Given that, by construction, they connect neighboring collapse points, there are nn distinct such cycles a\mathfrak{C}_{a}^{\prime}. Combining them with the n+1n+1 cycles a\mathfrak{C}_{a}, we can thus build a basis consisting of 2n+12n+1 non-contractible four-cycles in the solutions defined by eq. (47).

Σ2\Sigma_{2}zzρ\rhoν2a3\nu_{2a-3}ν2a1\nu_{2a-1}ν2a+1\nu_{2a+1}a\mathfrak{C}^{\prime}_{a}a\mathfrak{C}_{a}
Figure 3: The profile of the non-contractible four-cycles a\mathfrak{C}_{a} and a\mathfrak{C}^{\prime}_{a} in Σ2\Sigma_{2}. Every point on the red and blue curves is a 3-sphere.

Having identified a basis for non-contractible four-cycles, we are now in a position to derive the M-brane charges by integrating the four-form flux over the a\mathfrak{C}_{a}. To ease the computation, we abstractly write the flux as Faedo:2020nol

(4)=2h1volAdS3dz+2h1vol𝕊3dz+zHρ3dρvol𝕊~3ρHρ3dzvol𝕊~3.\begin{split}\mathcal{F}_{(4)}&=2h_{1}\operatorname{vol}_{\text{AdS}_{3}}\wedge\text{d}z+2h_{1}\operatorname{vol}_{\mathds{S}^{3}}\wedge\text{d}z\\ &+\partial_{z}H\rho^{3}\text{d}\rho\wedge\operatorname{vol}_{\tilde{\mathds{S}}^{3}}-\partial_{\rho}H\rho^{3}\text{d}z\wedge\operatorname{vol}_{\tilde{\mathds{S}}^{3}}~{}.\end{split} (62)

The pull-back of the four-form field strength (4)\mathcal{F}_{(4)} onto any of the a\mathfrak{C}_{a} cycles eliminates the first two terms in eq. (62). Integrating along a given a\mathfrak{C}_{a} then yields

aPa[(4)]=2h1Vol(𝕊~3)(ξ^2aξ^2a1).\int_{\mathfrak{C}_{a}}P_{\mathfrak{C}_{a}}[\mathcal{F}_{(4)}]=2h_{1}\operatorname{Vol}({\tilde{\mathds{S}}^{3}})(\hat{\xi}_{2a}-\hat{\xi}_{2a-1})\,. (63)

The choice of orientation we made while defining the cycles a\mathfrak{C}_{a} ensures that the integral above is positive. This enables us to define the following charge

Ma\displaystyle M_{a} =12(4π2GN)1/3aPa[(4)],\displaystyle=\frac{1}{2(4\pi^{2}G_{N})^{1/3}}\int_{\mathfrak{C}_{a}}P_{\mathfrak{C}_{a}}[\mathcal{F}_{(4)}]\,, (64)

which can be interpreted as the number of M5-branes in the atha^{\rm th} stack Bachas:2013vza , which obey aMa=M\sum_{a}M_{a}=M with MM being the total number of M5-branes.

It is also possible to generalize the construction above. As we mentioned previously, we can in fact build four-cycles surrounding more than one collapse point νa\nu_{a}. The four-cycle defined by bca=bca\mathfrak{C}_{bc}\equiv\sum_{a=b}^{c}\mathfrak{C}_{a}, with 1bcn+11\leq b\leq c\leq n+1, is also non-contractible by construction, and is characterized by a charge

bcPbc[(4)]=2h1Vol(𝕊~3)a=bc(ξ^2aξ^2a1),\int_{\mathfrak{C}_{bc}}P_{\mathfrak{C}_{bc}}[\mathcal{F}_{(4)}]=2h_{1}\operatorname{Vol}({\tilde{\mathds{S}}^{3}})\sum_{a=b}^{c}(\hat{\xi}_{2a}-\hat{\xi}_{2a-1})\,, (65)

under the four-form field strength.

We can follow a similar analysis for the flux threading the other set of non-contractible four-cycles, which we labelled a\mathfrak{C}^{\prime}_{a} earlier. In this case, only the second term in eq. (62) provides a non-vanishing contribution after pulling back the four-form field strength (4)\mathcal{F}_{(4)} to the cycle a\mathfrak{C}^{\prime}_{a}. Integrating the flux through a\mathfrak{C}^{\prime}_{a} gives

aPa[(4)]=2h1Vol(𝕊3)(ν2a+1ν2a1).\int_{\mathfrak{C}^{\prime}_{a}}P_{\mathfrak{C}^{\prime}_{a}}[\mathcal{F}_{(4)}]=2h_{1}\operatorname{Vol}({\mathds{S}^{3}})(\nu_{2a+1}-\nu_{2a-1})~{}. (66)

Similar to the integrated fluxes through the a\mathfrak{C}_{a} cycles, we define the charge

Ma\displaystyle M^{\prime}_{a} =12(4π2GN)1/3aPa[(4)],\displaystyle=\frac{1}{2(4\pi^{2}G_{N})^{1/3}}\int_{\mathfrak{C}^{\prime}_{a}}P_{\mathfrak{C}^{\prime}_{a}}[\mathcal{F}_{(4)}]\,, (67)

which is read as the number of M5-branes in the atha^{\rm th} stack (see Table 1).

Again, we can easily generalize this analysis to four-cycles that connect non-neighboring collapse points νa\nu_{a}. Denoting the sum of four-cycles as bca=bca\mathfrak{C}^{\prime}_{bc}\equiv\sum_{a=b}^{c}\mathfrak{C}^{\prime}_{a}, where given the construction of the a\mathfrak{C}_{a}^{\prime} in eq. (66) 1bcn1\leq b\leq c\leq n, the integral of the flux through this cycle is simply

P[(4)]=2h1Vol(𝕊3)a(ν2a+1ν2a1).\int_{\mathfrak{C}^{\prime}}P_{\mathfrak{C}^{\prime}}[\mathcal{F}_{(4)}]=2h_{1}\operatorname{Vol}({\mathds{S}^{3}})\sum_{a}(\nu_{2a+1}-\nu_{2a-1})\,. (68)

In addition, following Bachas:2013vza we can deduce the number of M2-branes ending on the atha^{\text{th}} stack of M5-branes, which we denote

Na=b=anMb=h1Vol(𝕊3)(4π2GN)1/3b=an(ν2b+1ν2b1).N_{a}=\sum_{b=a}^{n}M^{\prime}_{b}=\frac{h_{1}\operatorname{Vol}({\mathds{S}^{3}})}{(4\pi^{2}G_{N})^{1/3}}\sum_{b=a}^{n}(\nu_{2b+1}-\nu_{2b-1})\,. (69)

One may notice how the definition of the collapse points νa\nu_{a} influences the properties of NaN_{a}. Indeed, since {νa}\{\nu_{a}\} is an ordered set, we see that NaNbN_{a}\geq N_{b} for aba\leq b. In other words, the set {Na}\{N_{a}\} forms a partition of the total number of M2-branes N=aNaN=\sum_{a}N_{a}, and so it is possible to define a Young diagram to encode the brane charges as illustrated in figure 4(a).

N1N_{1}M1M_{1}M2M_{2}M3M_{3}M4M_{4}N4N_{4}
(a)
N1N_{1}M1M_{1}M2M_{2}M3M_{3}N3N_{3}
(b)
Figure 4: (a) Young diagram corresponding to the partition specifying a γ\gamma\to-\infty solution with 55 distinct νa\nu_{a} constructed by pairwise collapse of 10 different ξj\xi_{j}. (b) Young diagram obtained by “multiwise” collapse of 10 different ξj\xi_{j} to 44 distinct νa\nu_{a}. Another way to realize the construction of (b) starts from the partition in (a) and collapses ν2=ν3\nu_{2}=\nu_{3}.

We can now recast the various moments defined in eqs. 52 and 97 in terms of this partition data. Since only m^1\hat{m}_{1} appears in the asymptotic expansions in appendix A, and so also in the physical quantities computed using those expansions, we will not treat the higher m^j\hat{m}_{j} moments and simply write

m^1=j=12n+2(1)jξ^j=(4π2GN)1/3h1Vol(𝕊~3)M,\hat{m}_{1}=\sum_{j=1}^{2n+2}(-1)^{j}\hat{\xi}_{j}=\frac{(4\pi^{2}G_{N})^{1/3}}{h_{1}\operatorname{Vol}({\tilde{\mathds{S}}^{3}})}M\,, (70)

which gives back the usual relation between MM and the length scale on the 𝕊4\mathds{S}^{4},

L𝕊43=h1m^12=(GN)1/3(2π)4/3M.L^{3}_{\mathds{S}^{4}}=\frac{h_{1}\hat{m}_{1}}{2}=\frac{(G_{N})^{1/3}}{(2\pi)^{4/3}}M\,. (71)

The n^j\hat{n}_{j} moments are a bit trickier to re-express in terms of MaM_{a} and NaN_{a}, but one can show that the following relation holds for any jj

k=0j(1)j+kj!k!(jk)!ν2n+1kn^jk=(GN1/3h1π(2π)1/3)j+1a=1n+1MaNaj,\displaystyle\sum_{k=0}^{j}(-1)^{j+k}\frac{j!}{k!(j-k)!}\nu_{2n+1}^{k}~{}\hat{n}_{j-k}=\left(\frac{G_{N}^{1/3}}{h_{1}\pi(2\pi)^{1/3}}\right)^{j+1}~{}\sum_{a=1}^{n+1}M_{a}N_{a}^{j}~{}, (72)

where it is understood that n^0=m^1\hat{n}_{0}=\hat{m}_{1}. However, we will only require relations up to j=2j=2 moving forward. Using the expressions above, we can write

n^12m^1n^2=GN4/3h14π4(2π)4/3[(a=1n+1MaNa)2Ma=1n+1MaNa2],\displaystyle\hat{n}_{1}^{2}-\hat{m}_{1}\hat{n}_{2}=\frac{G_{N}^{4/3}}{h_{1}^{4}\pi^{4}(2\pi)^{4/3}}\left[\left(\sum_{a=1}^{n+1}M_{a}N_{a}\right)^{2}-M\sum_{a=1}^{n+1}M_{a}N_{a}^{2}\right]~{}, (73)

which will be useful in the following section.

Before moving on, we recall from our previous discussion that it is also possible to consider solutions where more than two branch points ξj\xi_{j} collapse to a single point νa\nu_{a} and, of course, build non-contractible four-cycles around or connecting them. Let us again index the pp distinct loci of collapse as νa\nu_{a}, where now 1ap+11\leq a\leq p+1 with pnp\leq n. The limiting case p=np=n recovers the pairwise collapse described above. It will be useful to label IaI_{a} and KaK_{a} respectively as the smallest and largest jj-indices of the ξj\xi_{j} branch points which collapse to a given νa\nu_{a}, i.e. νa:=νIa=νIa+1==νKa\nu_{a}:=\nu_{I_{a}}=\nu_{I_{a}+1}=\cdots=\nu_{K_{a}}. The ordering of the collapse points νa\nu_{a} and of the branch points ξj\xi_{j} was discussed previously in section 4; in particular, we recall that the parameters ξ^j\hat{\xi}_{j}’s associated to the branch points collapsing to the same νa\nu_{a} are ordered amongst themselves. In the end, the analysis for these “multiwise” collapse solutions is identical to the one presented above for the pairwise collapse scenario, up to the replacements

ξ^2aξ^2a1j=IaKa(1)jξ^j\displaystyle\hat{\xi}_{2a}-\hat{\xi}_{2a-1}\quad\longrightarrow\quad\sum_{j=I_{a}}^{K_{a}}(-1)^{j}\hat{\xi}_{j}~{} (74)

and npn\to p throughout the expressions above. The net effect is that the partition data yields a differently shaped Young diagram, illustrated in figure 4(b). In particular, since the multiwise collapse can greatly reduce the number of singular points on the boundary of the Riemann surface, we see that the Young tableaux specifying the γ\gamma\to-\infty solutions have at most pp rows. This is a striking difference compared to the Young tableau construction of Estes:2018tnu for the Wilson surface solutions with γ=1/2\gamma=-1/2, whose associated Young tableaux always have nn rows. This difference is maximal in the 1PV solution described in section 4.2, to which one cannot attach any Young tableau interpretation at all. This is due to the fact that the 1PV solution features a single collapse point ν\nu, so that the construction of the a{\mathfrak{C}}_{a}^{\prime} cycle fails, thus preventing the definition of the MaM_{a}^{\prime} and NaN_{a} charges.

5 Entanglement entropy of small 𝓝=(𝟒,𝟒)\mathcal{N}=(4,4) surface defects

The entanglement entropy SEES_{\text{EE}} of a spatial subregion \mathcal{B} within a QFT is defined as the von Neumann entropy of the reduced density matrix obtained by tracing out the states in the complementary region ¯\overline{\mathcal{B}} of the QFT. For CFTs with weakly coupled gravity duals, the RT prescription Ryu:2006bv ; Ryu:2006ef ; Nishioka:2021uef holographically recasts the computation of SEES_{\text{EE}} into the following Plateau problem in the asymptotically AdS bulk,

SEE=minζ𝒜[ζ]4GN,S_{\text{EE}}=\min_{\zeta}\frac{\mathcal{A}[\zeta]}{4G_{\text{N}}}~{}, (75)

where 𝒜[ζ]\mathcal{A}[\zeta] is the area functional evaluated on a bulk hypersurface ζ\zeta which is homologous to the chosen spatial subregion in the dual CFT: ζ=b\zeta\cup\mathcal{B}=\partial b for some static bulk subregion bb. We will denote this extremal bulk hypersurface as ζRT\zeta_{\text{RT}}. In the computations below, we choose the spatial subregion to be an Euclidean 5-ball, =𝔹R55\mathcal{B}=\mathbb{B}_{R}^{5}\hookrightarrow\mathbb{R}^{5}, with radius RR. We take \mathcal{B} to be centered on the spatial extent of the surface defect, which has a Lorentzian worldvolume Υ2=1,1\Upsilon_{2}=\mathbb{R}^{1,1}.

In an ordinary QFT, the presence of highly entangled UV degrees of freedom induces short-distance divergences near the surface =𝕊4\partial\mathcal{B}=\mathds{S}^{4}. Most of the divergences in the EE of a general QFT are non-universal and shape-dependent. However, in even dimensional theories, there are universal log-divergent contributions to the EE, which are generically related at conformal fixed points to the Weyl anomalies of the CFT.

This is true in the presence of a defect as well, but we see additional divergences that arise from the defect degrees of freedom near Υ2\Upsilon_{2}\cap\partial\mathcal{B}. In order to isolate these defect-localized contributions, we will adopt a scheme where we subtract off the EE of the deformed, vacuum ambient CFT, SEE[]S_{\rm EE}[\emptyset], from the EE computed in the presence of the defect SEE[Υ2]S_{\rm EE}[\Upsilon_{2}]. We can then extract the coefficient of the universal, log-divergent part of the defect contribution to the sphere EE by

SEE(univ)\displaystyle S_{\rm EE}^{\rm(univ)} =RddR(SEE[Υ2]SEE[])|R0.\displaystyle=R\frac{\text{d}}{\text{d}R}\left(S_{\text{EE}}[\Upsilon_{2}]-S_{\text{EE}}[\emptyset]\right)\big{|}_{R\to 0}~{}. (76)

In Jensen:2018rxu , it was shown that contribution to the log-divergent part of the EE of a spherical region coming from a flat, 22d conformal defect embedded in a dd-dimensional flat-space ambient CFT takes the form of a linear combination of defect localized Weyl anomalies. Explicitly, for a 66d ambient CFT

SEE(univ)=13(aΥ35d2).\displaystyle S_{\rm EE}^{\rm(univ)}=\frac{1}{3}\left(a_{\Upsilon}-\frac{3}{5}d_{2}\right)~{}. (77)

where aΥa_{\Upsilon} is the A-type defect Weyl anomaly coefficient–i.e. appearing with the intrinsic Euler density–and d2d_{2} is the B-type anomaly that enters with the trace of the pullback of the ambient Weyl tensor. Importantly, while Jensen:2018rxu demonstrated that d20d_{2}\geq 0 based on energy conditions, aΥa_{\Upsilon}, though obeying a defect “c-theorem”, has no positivity constraints. This means that as opposed to an ordinary, unitary 22d CFT where the universal part of the EE is is proportional to the central charge Holzhey:1994we and, hence, is non-negative, it is clear from eq. (77) that SEE(univ)S_{\rm EE}^{\rm(univ)} is not similarly bounded nor is it RG monotonic.

In completing the holographic computation, we also need to contend with the fact that the FG expansion is not globally defined as it typically breaks down in a region near the AdS submanifold that is dual to the insertion of the defect in the field theory. However, we have a full analysis of the asymptotic expansions of the data specifying the γ\gamma\to-\infty solutions in appendix A, and we have the general prescription for the FG transformation suitable for defect EE in Estes:2014hka . Together, we will be able to unambiguously holographically compute SEE(univ)S_{\rm EE}^{\rm(univ)}.

To begin the holographic computation of the defect EE, we choose the following parametrization for the AdS3 subspace in eq. (51),

dsAdS32=1u2(du2dt2+dx2).\text{d}s^{2}_{\text{AdS}_{3}}=\frac{1}{u^{2}}\left(\text{d}u^{2}-\text{d}t^{2}+\text{d}x^{2}_{\parallel}\right). (78)

Furthermore, we take the RT hypersurface ζ\zeta to wrap both 𝕊3\mathds{S}^{3} and 𝕊~3\tilde{\mathds{S}}^{3}, and its profile in the remaining subspace to be described by x(u,ρ,z)x_{\parallel}(u,\rho,z). The area of ζ\zeta as measured against the metric in eq. (51) is then

A[ζ]=Vol(𝕊3)Vol(𝕊~3)duΣ2dρdz,A[\zeta]=\operatorname{Vol}(\mathds{S}^{3})\operatorname{Vol}(\tilde{\mathds{S}}^{3})\int\text{d}u\int_{\Sigma_{2}}\ \text{d}\rho\text{d}z\ \mathcal{L}~{}, (79)

where the Lagrangian is

\displaystyle\mathcal{L} =h12ρ3u2[h1((ρx)2+(zx)2)H(z,ρ)+u2(1+(ux)2)H(z,ρ)2]1/2.\displaystyle=\frac{h_{1}^{2}\rho^{3}}{u^{2}}\left[h_{1}\left((\partial_{\rho}x_{\parallel})^{2}+(\partial_{z}x_{\parallel})^{2}\right)H(z,\rho)+u^{2}\left(1+(\partial_{u}x_{\parallel})^{2}\right)H(z,\rho)^{2}\right]^{1/2}. (80)

As shown in Jensen:2013lxa , the minimal area surface wraps the Riemann surface Σ2\Sigma_{2} too, so that ρx=zx=0\partial_{\rho}x_{\parallel}=\partial_{z}x_{\parallel}=0. The Lagrangian is thus minimized by

x2+u2=R2,x_{\parallel}^{2}+u^{2}=R^{2}, (81)

for a constant RR, so that the area of the extremal hypersurface ζRT\zeta_{\text{RT}} is

A[ζRT]=h13Vol(𝕊3)Vol(𝕊~3)log(2Rϵu)Σ2dρdzρ3j=12n+2(1)jξ^j(ρ2+(zνj)2)3/2+𝒪(ϵu2),A[\zeta_{\text{RT}}]=h_{1}^{3}\operatorname{Vol}(\mathds{S}^{3})\operatorname{Vol}(\tilde{\mathds{S}}^{3})\log\left(\frac{2R}{\epsilon_{u}}\right)\int_{\Sigma_{2}}\text{d}\rho\text{d}z\ \rho^{3}\sum_{j=1}^{2n+2}\frac{(-1)^{j}\hat{\xi}_{j}}{\left(\rho^{2}+(z-\nu_{j})^{2}\right)^{3/2}}+\mathcal{O}(\epsilon_{u}^{2}), (82)

where we introduced a small-uu cut-off ϵu>0\epsilon_{u}>0.

The evaluation of the integral above is performed in detail in appendix B. The resulting holographic entanglement entropy is

SEE[Υ2]=π4L𝕊49GNlog(2Rϵu)[6431ϵv4+165n^12m^14165n^2m^13+𝒪(ϵv2)]+𝒪(ϵu2),S_{\text{EE}}[\Upsilon_{2}]=\frac{\pi^{4}L_{\mathds{S}^{4}}^{9}}{G_{\text{N}}}\log\left(\frac{2R}{\epsilon_{u}}\right)\left[\frac{64}{3}\frac{1}{\epsilon_{v}^{4}}+\frac{16}{5}\frac{\hat{n}_{1}^{2}}{\hat{m}_{1}^{4}}-\frac{16}{5}\frac{\hat{n}_{2}}{\hat{m}_{1}^{3}}+\mathcal{O}(\epsilon_{v}^{2})\right]+\mathcal{O}(\epsilon_{u}^{2}), (83)

where ϵv>0\epsilon_{v}>0 is a small-vv cutoff in the FG parametrization.

Subtracting off the 1PV contribution to the entanglement entropy, SEE1PVS_{\text{EE}}^{\text{1PV}}, precisely removes the ϵv4\epsilon_{v}^{-4} divergence from eq. (83), and leaves the 𝒪(ϵv0)\mathcal{O}(\epsilon_{v}^{0}) term unchanged. To see this, we recall that the 1PV limit takes n^k0\hat{n}_{k}\to 0, which in eq. (83) gives

SEE1PV=π4L𝕊49GNlog(2Rϵu)[6431ϵv4+𝒪(ϵv2)]+𝒪(ϵu2).\displaystyle S_{\text{EE}}^{\rm 1PV}=\frac{\pi^{4}L_{\mathds{S}^{4}}^{9}}{G_{\text{N}}}\log\left(\frac{2R}{\epsilon_{u}}\right)\left[\frac{64}{3}\frac{1}{\epsilon_{v}^{4}}+\mathcal{O}(\epsilon_{v}^{2})\right]+\mathcal{O}(\epsilon_{u}^{2}). (84)

Therefore, we can at once compute the coefficient of the universal part of the defect sphere EE using eq. (76) and plugging in eqs. (83) and (84) with SEE[]=SEE1PVS_{\rm EE}[\emptyset]=S_{\rm{EE}}^{\rm 1PV} to find

SEE(univ)\displaystyle S_{\rm EE}^{\rm(univ)} =165π4L𝕊49GNn^12m^1n^2m^14\displaystyle=\frac{16}{5}\frac{\pi^{4}L_{\mathds{S}^{4}}^{9}}{G_{\text{N}}}\frac{\hat{n}_{1}^{2}-\hat{m}_{1}\hat{n}_{2}}{\hat{m}_{1}^{4}} (85a)
=15M[(a=1n+1MaNa)2Ma=1n+1MaNa2],\displaystyle=\frac{1}{5M}\left[\left(\sum_{a=1}^{n+1}M_{a}N_{a}\right)^{2}-M\sum_{a=1}^{n+1}M_{a}N_{a}^{2}\right], (85b)

where we have mapped to field theory quantities using L𝕊43=GN1/3(2π)4/3ML_{\mathds{S}^{4}}^{3}=\frac{G_{N}^{1/3}}{(2\pi)^{4/3}}M and used the definitions of moments n^j\hat{n}_{j}, m^j\hat{m}_{j} in terms of the numbers of branes in eqs. (70) and (73). In terms of the highest weight ϖ\varpi of the AM1A_{M-1} irreducible representation encoded in the Young diagrams that specify the defect discussed in the previous section, we can re-express the defect sphere EE as Estes:2018tnu

SEE(univ)=(ϖ,ϖ)5,\displaystyle S_{\rm EE}^{\rm(univ)}=-\frac{(\varpi,\varpi)}{5}~{}, (86)

where (,)(\cdot,\cdot) is the scalar product on the weight space induced by the Killing form.

We also note that the contribution of the aforementioned bulk deformation to the coefficient of the universal, log-divergent component of the vacuum-subtracted entanglement entropy is

SEE(univ,bulkdef.)\displaystyle S_{\rm EE}^{\rm(univ,bulk-def.)} =83π4L𝕊49GN=M36.\displaystyle=\frac{8}{3}\frac{\pi^{4}L_{\mathds{S}^{4}}^{9}}{G_{\text{N}}}~{}=\frac{M^{3}}{6}~{}. (87)

This is independent of the moments (m^i,n^j)(\hat{m}_{i},\hat{n}_{j}), in line with the lack of Young Tableau data associated to the bulk deformation. Had we subtracted in eq. (76) the AdS×7𝕊4{}_{7}\times\mathds{S}^{4} vacuum, rather than the 1PV, the resulting entanglement entropy would have received both the defect and bulk deformation contributions above.

Finally, we note that the same quantity can be trivially computed in the orbifolded theory described in section 3.3 simply by rescaling

SEE(univ)Vol(𝕊3/k)Vol(𝕊~3/k)Vol(𝕊3)Vol(𝕊~3)SEE(univ)=SEE(univ)kk.\displaystyle S_{\rm EE}^{\rm(univ)}\longrightarrow\frac{\operatorname{Vol}(\mathds{S}^{3}/\mathbb{Z}_{k})\operatorname{Vol}(\tilde{\mathds{S}}^{3}/\mathbb{Z}_{k^{\prime}})}{\operatorname{Vol}(\mathds{S}^{3})\operatorname{Vol}(\tilde{\mathds{S}}^{3})}S_{\rm EE}^{\rm(univ)}=\frac{S_{\rm EE}^{\rm(univ)}}{kk^{\prime}}~{}. (88)

6 Summary and Outlook

In this work, we have constructed a novel class of solutions in 11d SUGRA that are holographically dual to 2d superconformal defects preserving small 𝒩=(4,4)\mathcal{N}=(4,4) and 𝒩=(0,4)\mathcal{N}=(0,4) SUSY in 6d SCFTs at large MM. These solutions fit into the one-parameter family organized in a general classification scheme of 11d SUGRA solutions with superisometry 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}(2,1;\gamma)\oplus\mathfrak{d}(2,1;\gamma) Bachas:2013vza ; specifically, they are obtained in the γ\gamma\to-\infty limit. There are several features of these new solutions that separate them from the more familiar γ=1/2\gamma=-1/2 case that holographically corresponds to 1/21/2-BPS Wilson surface type defects in the 6d 𝒩=(2,0)\mathcal{N}=(2,0) AM1A_{M-1} SCFT.

Within the one-parameter family of solutions labelled by γ\gamma, the γ\gamma\to-\infty limit is slightly unusual from the superalgebra perspective. Despite producing an AdS×7𝕊4{}_{7}\times\mathds{S}^{4} asymptotic geometry as shown in section 3, taking the γ\gamma\to-\infty limit means that 𝔡(2,1;γ)𝔡(2,1;γ)\mathfrak{d}(2,1;\gamma)\oplus\mathfrak{d}(2,1;\gamma) is not a subalgebra of the 𝔬𝔰𝔭(8|4)\mathfrak{osp}(8^{*}|4) superisometry of AdS×7𝕊4{}_{7}\times\mathds{S}^{4}. On the field theory side of the holographic duality, this means that the ambient theory into which the defects are inserted is some deformation of the 6d 𝒩=(2,0)\mathcal{N}=(2,0) SCFT.

We have seen that choosing all of the singular loci in the internal space to collapse to a single point – a configuration which we label 1PV – destroys the data that specifies the defect, i.e. the Young diagram corresponding to the arrangement of M5-branes. However, as is clear from the discussion in section 4.2, the vacuum that we arrive at has an isometry group of SO(2,2)×SO(3)×SO(5)SO(2,2)\times SO(3)\times SO(5), as opposed to the vacuum solution at γ=1/2\gamma=-1/2, which instead enjoys the full SO(6,2)SO(6,2). In the latter case, the trivial defect corresponds to a Wilson surface transforming in the 𝟏\bf{1} of AM1A_{M-1}, and the conformal symmetry of the ambient 6d 𝒩=(2,0)\mathcal{N}=(2,0) theory is restored. For γ\gamma\to-\infty, the trivial defect dual to the 1PV still possesses what looks like the ‘defect’ conformal symmetry despite being the vacuum solution, which renders giving a precise definition for and interpretation of the defect CFT difficult. The 1PV does, however, enable us to correctly employ a background subtraction scheme101010Here, “correctly” refers to a background subtraction which also removes any contributions from the trivial defect. As we have demonstrated, the same cannot be said of a subtraction scheme which utilizes vacuum AdS×7𝕊4{}_{7}\times\mathds{S}^{4}. and arrive at a finite result for SEE(univ)S_{\rm EE}^{\rm(univ)} and, we believe, resolves the puzzling appearance of divergences in the “defect central charge” computed in Faedo:2020nol .

On a more fundamental level, the small 𝒩=(4,4)\mathcal{N}=(4,4) defects at γ\gamma\to-\infty cannot be viewed as a smooth deformation of the Wilson surface defects at γ=1/2\gamma=-1/2. Indeed, the γ\gamma\to-\infty solutions cannot even be smoothly deformed into the solution at γ=0\gamma=0, to which they are related by the involution γ1/γ\gamma\mapsto 1/\gamma with an exchange of the 𝔰𝔬(3)𝔰𝔬(3)\mathfrak{so}(3)\oplus\mathfrak{so}(3) factors in 𝔡(2,1;γ;0)𝔡(2,1;γ;0)\mathfrak{d}(2,1;\gamma;0)\oplus\mathfrak{d}(2,1;\gamma;0). The reason is that there is a special point at γ=1\gamma=-1 where the real form 𝔡(2,1;γ;0)\mathfrak{d}(2,1;\gamma;0) becomes 𝔬𝔰𝔭(4|2;)\mathfrak{osp}(4|2;\mathbb{R}). At this value of γ\gamma, SO(2,2)SO(2,2) Wigner-İnönü contracts to ISO(1,2)ISO(1,2), and AdS3 becomes 2,1\mathbb{R}^{2,1}. Therefore, the γ\gamma\to-\infty solutions are isolated from the other class of asymptotically AdS×7𝕊4{}_{7}\times\mathds{S}^{4} geometries.

In light of the new small 𝒩=(4,4)\mathcal{N}=(4,4) solutions that we have constructed and holographically studied, there are a number of open questions that remain to be answered.

Firstly, as we discussed at the start of section 5 , the contribution from a flat 2d conformal defect to the log-divergent, universal part of the EE of a spherical region in a d4d\geq 4 ambient CFT is built from a linear combination of two defect Weyl anomaly coefficients, aΥa_{\Upsilon} and d2d_{2} that characterize the defect theory. In order to disentangle these two fundamental defect quantities, we would compute a second holographic quantity that contains either aΥa_{\Upsilon} or d2d_{2}. For instance, d2d_{2} controls the normalization of the one-point function Tμν\left<T_{\mu\nu}\right> of the stress-energy tensor, which can be readily computed for most 10d or 11d supergravity solutions by dimensional reduction on the internal space deHaro:2000vlm ; Skenderis:2006uy . Therefore, it is natural to try to compute d2d_{2} and SEE(univ)S_{\rm EE}^{\rm(univ)} in order to isolate the independent defect Weyl anomalies111111In fact, aΥa_{\Upsilon} and d2d_{2} are the only independent defect Weyl anomaly coefficients for superconformal defects preserving at least 2d 𝒩=(0,2)\mathcal{N}=(0,2) supersymmetry. This was shown for co-dimension four defects in 6d SCFTs in Drukker:2020atp and co-dimension two defects in 4d SCFT in Bianchi:2019sxz . . This was successfully done for the Wilson surfaces at γ=1/2\gamma=-1/2 in Estes:2018tnu and for codimension-2 defects in Capuozzo:2023fll . However, the γ\gamma\to-\infty solutions are more subtle, and a naïve application of dimensional reduction techniques would be inappropriate. Namely, in reducing the 11d solutions to 7d, the presence of the non-trivial four-form flux modifies the gravitational equations of motion at the conformal boundary of AdS7, which violates the assumptions in deHaro:2000vlm . Therefore, computing d2d_{2} holographically from Tμν\left<T_{\mu\nu}\right> requires a generalization to account for flux contributions, which is the subject of ongoing work.

Furthermore, it is natural to look for physical observables which can be reliably computed on the field theory side and employed to test the holographic predictions made above. For 2d BPS conformal defects in 6d AM1A_{M-1} and DMD_{M} 𝒩=(2,0)\mathcal{N}=(2,0) SCFTs at large MM, recent developments in analytic bootstrap methods have enabled the computation of correlations functions in the presence of 2d defects that are controlled by anomalies Drukker:2020atp ; Meneghelli:2022gps . Further, despite the lack of a Lagrangian description and of supersymmetric localization methods for 6d SCFTs at large MM, chiral algebra methods have also been shown to give exact results for defect correlators Meneghelli:2022gps and the defect SUSY Casimir energy Bullimore:2014upa ; Chalabi:2020iie . Currently, only Wilson surface type defects, i.e. the holographically dual theories to the γ=1/2\gamma=-1/2 solutions, have been studied using these field theory techniques. It is reasonable to wonder whether any of these methods are applicable to the types of defects in the deformed 6d theory that we have constructed in this work.

The biggest hurdle to clear in trying to generalize bootstrap or chiral algebra methods for use in the dual to the 1PV of the γ\gamma\to-\infty solutions is clarifying the precise role of the deformation parameter γ\gamma. As we have explained in the γ\gamma\to-\infty limit, the ambient 6d theory has reduced global and conformal symmetries, and there is no smooth path in field theory space as γ\gamma is varied through γ=1\gamma=-1 to get from the 6d AM1A_{M-1} 𝒩=(2,0)\mathcal{N}=(2,0) theory to the dual of the 1PV. It is unclear at the moment precisely what symmetry breaking operators are sourced on the field theory side in the γ\gamma\to-\infty limit.

Acknowledgments

The authors would like to thank Yolanda Lozano and Nicolò Petri for useful discussions during the completion of this work and for their comments on a draft of this manuscript. The authors would also like to thank Andy O’Bannon for his contributions during the early phase of this research. The work of PC is supported by a Mayflower studentship from the University of Southampton. The work of BR is supported by the INFN. The work of BS is supported in part by the STFC consolidated grant ST/T000775/1. This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of High Energy Physics under Award Number DE-SC0024557.

Disclaimer: “This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.”

Appendix A Fefferman-Graham parametrization

In this Appendix, we will be concerned with finding a reparametrization {w,w¯}{v,ϕ}\{w,\bar{w}\}\to\{v,\phi\} of the Riemann surface Σ2\Sigma_{2} under which the 11d SUGRA solution in eq. (18) can be asymptotically (in particular, for small vv) recast into the following form,

ds2=4L𝕊42v2[dv2+α1dsAdS32+α2ds𝕊32]+L𝕊42[α3dϕ2+α4sϕ2ds𝕊~32],\text{d}s^{2}=\frac{4L_{\mathds{S}^{4}}^{2}}{v^{2}}\left[\text{d}v^{2}+\alpha_{1}\text{d}s^{2}_{\text{AdS}_{3}}+\alpha_{2}\text{d}s^{2}_{\mathds{S}^{3}}\right]+L_{\mathds{S}^{4}}^{2}\left[\alpha_{3}\text{d}\phi^{2}+\alpha_{4}s^{2}_{\phi}\text{d}s^{2}_{\tilde{\mathds{S}}^{3}}\right]~{}, (89)

for some 𝒪(v0)\mathcal{O}(v^{0}) metric factors {αi}i=14\{\alpha_{i}\}_{i=1}^{4}. By imposing that the metric factor in the vv direction be exact in vv, as above, we can reconstruct the asymptotic reparametrization order by order in vv and in terms of polar coordinates

r=z2+ρ2andθ=arctan(ρ/z)\displaystyle r=\sqrt{z^{2}+\rho^{2}}\quad\text{and}\quad\theta=\arctan(\rho/z) (90)

on Σ2\Sigma_{2} to be the following,

r(v,ϕ)\displaystyle\allowdisplaybreaks r(v,\phi) =2(γ+1)2m1γv2+(2γ+1)m1(c2ϕ3)24γ+m2cϕ2m1+[3γ2m22(7c2ϕ+1)4m13\displaystyle=-\frac{2(\gamma+1)^{2}m_{1}}{\gamma v^{2}}+\frac{(2\gamma+1)m_{1}\left(c_{2\phi}-3\right)}{24\gamma}+\frac{m_{2}c_{\phi}}{2m_{1}}+\Biggl{[}\frac{3\gamma^{2}m_{2}^{2}\left(7c_{2\phi}+1\right)}{4m_{1}^{3}} (91a)
γ2m3(5c2ϕ+3)m122γ(2γ+1)m2cϕsϕ2m1116(8γ(γ+1)+3)m1c2ϕ\displaystyle\qquad-\frac{\gamma^{2}m_{3}\left(5c_{2\phi}+3\right)}{m_{1}^{2}}-\frac{2\gamma(2\gamma+1)m_{2}c_{\phi}s_{\phi}^{2}}{m_{1}}-\frac{1}{16}(8\gamma(\gamma+1)+3)m_{1}c_{2\phi}
+3748γ(γ+1)m1+73m1192+19192(2γ+1)2m1c4ϕ]v248γ(γ+1)2+𝒪(v4),\displaystyle\qquad+\frac{37}{48}\gamma(\gamma+1)m_{1}+\frac{73m_{1}}{192}+\frac{19}{192}(2\gamma+1)^{2}m_{1}c_{4\phi}\Biggr{]}\frac{v^{2}}{48\gamma(\gamma+1)^{2}}+\mathcal{O}(v^{4}),
θ(v,ϕ)\displaystyle\theta(v,\phi) =ϕ+(2γ+1)m12cϕ+3γm212(γ+1)2m12sϕv2+[9γ2m22s2ϕ8m145γ2m3s2ϕ6m13\displaystyle=\phi+\frac{(2\gamma+1)m_{1}^{2}c_{\phi}+3\gamma m_{2}}{12(\gamma+1)^{2}m_{1}^{2}}s_{\phi}v^{2}+\Biggl{[}\frac{9\gamma^{2}m_{2}^{2}s_{2\phi}}{8m_{1}^{4}}-\frac{5\gamma^{2}m_{3}s_{2\phi}}{6m_{1}^{3}} (91b)
+cϕsϕ(5(2γ+1)2c2ϕ24γ(γ+1)7)48\displaystyle\qquad+\frac{c_{\phi}s_{\phi}\left(5(2\gamma+1)^{2}c_{2\phi}-24\gamma(\gamma+1)-7\right)}{48}
+γ(2γ+1)m2(3c2ϕ1)sϕ12m12]v416(γ+1)4+𝒪(v6),\displaystyle\qquad+\frac{\gamma(2\gamma+1)m_{2}\left(3c_{2\phi}-1\right)s_{\phi}}{12m_{1}^{2}}\Biggr{]}\frac{v^{4}}{16(\gamma+1)^{4}}+\mathcal{O}(v^{6})~{},

where we introduced the moments

mk:=j=12n+2(1)jξjk.m_{k}:=\sum_{j=1}^{2n+2}(-1)^{j}\xi_{j}^{k}~{}. (92)

Under this mapping, the metric takes the desired form of eq. (89) with the following metric factors

α1(γ)\displaystyle\allowdisplaybreaks\alpha_{1}(\gamma) =1+2γ+3(1+2γ)c2ϕ16(γ+1)2v2+[9(4(3γ13)γ+16γ2κ17)\displaystyle=1+\frac{2\gamma+3-(1+2\gamma)c_{2\phi}}{16(\gamma+1)^{2}}v^{2}+\left[9\left(4(3\gamma-13)\gamma+16\gamma^{2}\kappa-17\right)\right. (93a)
67(2γ+1)2c4ϕ+12(8(3γ+1)γ+20γ2κ+3)c2ϕ]v418432(γ+1)4+𝒪(v6),\displaystyle\quad\left.-67(2\gamma+1)^{2}c_{4\phi}+12\left(8(3\gamma+1)\gamma+20\gamma^{2}\kappa+3\right)c_{2\phi}\right]\frac{v^{4}}{18432(\gamma+1)^{4}}+\mathcal{O}(v^{6})~{},
α2(γ)\displaystyle\alpha_{2}(\gamma) =(γ+1)2γ2+2γ1(1+2γ)c2ϕ16γ2v2+[9(4(3γ+19)γ+16γ2κ+47)\displaystyle=\frac{(\gamma+1)^{2}}{\gamma^{2}}+\frac{2\gamma-1-(1+2\gamma)c_{2\phi}}{16\gamma^{2}}v^{2}+\left[9\left(4(3\gamma+19)\gamma+16\gamma^{2}\kappa+47\right)\right. (93b)
67(2γ+1)2c4ϕ+12(8(3γ+5)γ+20γ2κ+19)c2ϕ]v418432γ2(γ+1)2+𝒪(v6),\displaystyle\quad\left.-67(2\gamma+1)^{2}c_{4\phi}+12\left(8(3\gamma+5)\gamma+20\gamma^{2}\kappa+19\right)c_{2\phi}\right]\frac{v^{4}}{18432\gamma^{2}(\gamma+1)^{2}}+\mathcal{O}(v^{6})~{},
α3(γ)\displaystyle\alpha_{3}(\gamma) =1+(2γ+1)cϕ24(γ+1)2v2+[12(γ+1)γ24γ2κ9+6(2γ+1)2c2ϕ\displaystyle=1+\frac{(2\gamma+1)c_{\phi}^{2}}{4(\gamma+1)^{2}}v^{2}+\left[-12(\gamma+1)\gamma-24\gamma^{2}\kappa-9+6(2\gamma+1)^{2}c_{2\phi}\right. (93c)
+7(2γ+1)2c4ϕ]v4768(γ+1)4+𝒪(v6),\displaystyle\quad\left.+7(2\gamma+1)^{2}c_{4\phi}\right]\frac{v^{4}}{768(\gamma+1)^{4}}+\mathcal{O}(v^{6})~{},
α4(γ)\displaystyle\alpha_{4}(\gamma) =1+(2γ+1)(2c2ϕ+1)12(γ+1)2v2+[52(γ+1)γ24γ2κ19\displaystyle=1+\frac{(2\gamma+1)(2c_{2\phi}+1)}{12(\gamma+1)^{2}}v^{2}+\left[-52(\gamma+1)\gamma-24\gamma^{2}\kappa-19\right. (93d)
+79(2γ+1)2c4ϕ+12(2(γ+1)γ10γ2κ3)c2ϕ]v44608(γ+1)4+𝒪(v6),\displaystyle\quad\left.+79(2\gamma+1)^{2}c_{4\phi}+12\left(-2(\gamma+1)\gamma-10\gamma^{2}\kappa-3\right)c_{2\phi}\right]\frac{v^{4}}{4608(\gamma+1)^{4}}+\mathcal{O}(v^{6})~{},

where, to the order shown, the moments {mi}\{m_{i}\} enter the metric factors only in the combination

κ3m224m1m3m14.\kappa\equiv\frac{3m_{2}^{2}-4m_{1}m_{3}}{m_{1}^{4}}~{}. (94)

Note that κ\kappa is invariant under coordinate transformations {ξiξi+λ}\{\xi_{i}\mapsto\xi_{i}+\lambda\}, even though m2m_{2} and m3m_{3} are not individually. From this mapping, we also deduce the curvature scale of the asymptotic local 𝕊4\mathds{S}^{4} in FG gauge,

L𝕊43=|1+γ|3γ2h1m12.L_{\mathds{S}^{4}}^{3}=\frac{|1+\gamma|^{3}}{\gamma^{2}}\frac{h_{1}m_{1}}{2}~{}. (95)

We can now take the γ\gamma\to-\infty limit, accompanied by appropriate rescalings as described in section 4, to determine the asymptotic local behavior of the metric in eq. (36). In this limit, we find that

κγ212(n^12m^1n^2)m^14,\frac{\kappa}{\gamma^{2}}\to\frac{12(\hat{n}_{1}^{2}-\hat{m}_{1}\hat{n}_{2})}{\hat{m}_{1}^{4}}~{}, (96)

where we introduced the additional moments

n^k:=j=12n+2(1)jνjkξ^j.\hat{n}_{k}:=\sum_{j=1}^{2n+2}(-1)^{j}\nu_{j}^{k}\hat{\xi}_{j}~{}. (97)

Furthermore, we find that this limiting procedure trivializes the relative warping between AdS3 and 𝕊3\mathds{S}^{3}; that is, α1(γ)=α2(γ)\alpha_{1}(\gamma\to-\infty)=\alpha_{2}(\gamma\to-\infty). Overall, the limit γ\gamma\to-\infty produces the FG line element advertized in eq. (51), with the following metric factors,

α1(γ)\displaystyle\alpha_{1}(\gamma\to-\infty) =1+(3+5c2ϕ)(n^12m^1n^2)32m^14v4\displaystyle=1+\frac{(3+5c_{2\phi})(\hat{n}_{1}^{2}-\hat{m}_{1}\hat{n}_{2})}{32\hat{m}_{1}^{4}}v^{4} (98a)
5(1+7c2ϕ)(2n^133m^1n^1n^2+m^12n^3)cϕ144m^16v6+𝒪(v8),\displaystyle\quad-\frac{5(1+7c_{2\phi})(2\hat{n}_{1}^{3}-3\hat{m}_{1}\hat{n}_{1}\hat{n}_{2}+\hat{m}_{1}^{2}\hat{n}_{3})c_{\phi}}{144\hat{m}_{1}^{6}}v^{6}+\mathcal{O}(v^{8}),
α3(γ)\displaystyle\alpha_{3}(\gamma\to-\infty) =1+3(m^1n^2n^12)8m^14v4+5(2n^133m^1n^1n^2+m^12n^3)cϕ12m^16v6+𝒪(v8),\displaystyle=1+\frac{3(\hat{m}_{1}\hat{n}_{2}-\hat{n}_{1}^{2})}{8\hat{m}_{1}^{4}}v^{4}+\frac{5(2\hat{n}_{1}^{3}-3\hat{m}_{1}\hat{n}_{1}\hat{n}_{2}+\hat{m}_{1}^{2}\hat{n}_{3})c_{\phi}}{12\hat{m}_{1}^{6}}v^{6}+\mathcal{O}(v^{8}), (98b)
α4(γ)\displaystyle\alpha_{4}(\gamma\to-\infty) =1+(1+5c2ϕ)(m^1n^2n^12)16m^14v4\displaystyle=1+\frac{(1+5c_{2\phi})(\hat{m}_{1}\hat{n}_{2}-\hat{n}_{1}^{2})}{16\hat{m}_{1}^{4}}v^{4} (98c)
+5(5cϕ+7c3ϕ)(2n^133m^1n^1n^2+m^12n^3)144m^16v6+𝒪(v8),\displaystyle\quad+\frac{5(5c_{\phi}+7c_{3\phi})(2\hat{n}_{1}^{3}-3\hat{m}_{1}\hat{n}_{1}\hat{n}_{2}+\hat{m}_{1}^{2}\hat{n}_{3})}{144\hat{m}_{1}^{6}}v^{6}+\mathcal{O}(v^{8})~{},

and with 𝕊4\mathds{S}^{4} radius given by eq. (53).

Finally, we note that, if all points are taken to collapse to the origin, νj=0\nu_{j}=0 for j{1,2,,2n+2}j\in\{1,2,\ldots,2n+2\}, the asymptotic reparametrization in eq. (91) becomes exact,

r(v,ϕ)\displaystyle r(v,\phi) =2m^1v2\displaystyle=\frac{2\hat{m}_{1}}{v^{2}} (99a)
θ(v,ϕ)\displaystyle\theta(v,\phi) =ϕ.\displaystyle=\phi~{}. (99b)

Appendix B Area of the Ryu-Takayanagi hypersurface

In this appendix, we fill in the technical details of the computation of the holographic entanglement entropy in eq. (83).

We begin by treating the integral in eq. (82). In order to exploit the geometry of the internal space and facilitate an easier path to evaluating the remaining integral, we first transform to polar coordinates {r,θ}\{r,\,\theta\}, introduced in eq. (90), with θ[0,π]\theta\in[0,\pi]. We can then expand the summands in eq. (82) on the complete basis of L2(0,π)L^{2}(0,\pi) functions spanned by the Legendre polynomials, Pk(cθ)P_{k}(c_{\theta}). The harmonic function HH can then be written in two equivalent representations

H(r,θ)={h12j=12n+2(1)jξ^j|νj|3(k=0rkνjkPk(cθ))3,r[0,|νj|),h12j=12n+2(1)jξ^jr3(k=0rkνjkPk(cθ))3,r(|νj|,Λr(ϵv,0)],H(r,\theta)=\begin{cases}\frac{h_{1}}{2}\sum_{j=1}^{2n+2}(-1)^{j}\hat{\xi}_{j}\left|\nu_{j}\right|^{-3}\left(\sum_{k=0}^{\infty}r^{k}\nu_{j}^{-k}P_{k}(c_{\theta})\right)^{3},\quad r\in[0,|\nu_{j}|),\vspace{0.2cm}\\ \frac{h_{1}}{2}\sum_{j=1}^{2n+2}(-1)^{j}\hat{\xi}_{j}r^{-3}\left(\sum_{k=0}^{\infty}r^{-k}\nu_{j}^{k}P_{k}(c_{\theta})\right)^{3},\quad r\in(|\nu_{j}|,\Lambda_{r}(\epsilon_{v},0)],\end{cases} (100)

which converge when integrated in rr over their respective domains. We have also introduced in the second line of eq. (100) a cutoff scale Λr\Lambda_{r} at large rr, which from the transformation to FG gauge in eq. (91) can be expressed in a small ϵv\epsilon_{v} expansion as

Λr(ϵv,θ)=2m^1ϵv2+n^1cθm^1+(3+5c2θ)m^1n^2(5+3c2θ)n^1216m^13ϵv2+𝒪(ϵv4).\Lambda_{r}(\epsilon_{v},\theta)=\frac{2\hat{m}_{1}}{\epsilon_{v}^{2}}+\frac{\hat{n}_{1}c_{\theta}}{\hat{m}_{1}}+\frac{(3+5c_{2\theta})\hat{m}_{1}\hat{n}_{2}-(5+3c_{2\theta})\hat{n}_{1}^{2}}{16\hat{m}_{1}^{3}}\ \epsilon_{v}^{2}+\mathcal{O}(\epsilon_{v}^{4}). (101)

Thus, the remaining integral in the area functional can be partitioned into two separate contributions

A[ζRT]=32π4L𝕊49m^13log(2Rϵu)j=12n+2(1)jξ^j(j(1)+j(2))+𝒪(ϵu2),A[\zeta_{\text{RT}}]=\frac{32\pi^{4}L_{\mathds{S}^{4}}^{9}}{\hat{m}_{1}^{3}}\log\left(\frac{2R}{\epsilon_{u}}\right)\sum_{j=1}^{2n+2}(-1)^{j}\hat{\xi}_{j}\left(\mathcal{I}_{j}^{(1)}+\mathcal{I}_{j}^{(2)}\right)+\mathcal{O}(\epsilon_{u}^{2}), (102)

where both integrals

j(1)\displaystyle\mathcal{I}_{j}^{(1)} =0πdθ0|νj|drr4sθ3(1|νj|k=0(rνj)kPk(cθ))3,\displaystyle=\int_{0}^{\pi}\text{d}\theta\int_{0}^{\left|\nu_{j}\right|}\text{d}r\ r^{4}s^{3}_{\theta}\left(\frac{1}{\left|\nu_{j}\right|}\sum_{k=0}^{\infty}\left(\frac{r}{\nu_{j}}\right)^{k}P_{k}(c_{\theta})\right)^{3}, (103a)
j(2)\displaystyle\mathcal{I}_{j}^{(2)} =0πdθ|νj|Λr(ϵv,θ)drrsθ3(k=0(νjr)kPk(cθ))3,\displaystyle=\int_{0}^{\pi}\text{d}\theta\int_{|\nu_{j}|}^{\Lambda_{r}(\epsilon_{v},\theta)}\text{d}r\ r\ s^{3}_{\theta}\left(\sum_{k=0}^{\infty}\left(\frac{\nu_{j}}{r}\right)^{k}P_{k}(c_{\theta})\right)^{3}, (103b)

are convergent.

At first glance, eq. (103) may not seem to put us in a better position to evaluate the integral, but we can exploit the properties of Pk(cθ)P_{k}(c_{\theta}) over the interval θ[0,π]\theta\in[0,\pi]. In particular, we can make use of the orthogonality relation of the triple-product of Legendre polynomials

0πdθsθP1(cθ)P2(cθ)P3(cθ)=2(123000)2,\int_{0}^{\pi}\text{d}\theta\ s_{\theta}\ P_{\ell_{1}}(c_{\theta})P_{\ell_{2}}(c_{\theta})P_{\ell_{3}}(c_{\theta})=2\begin{pmatrix}\ell_{1}&\ell_{2}&\ell_{3}\\ 0&0&0\end{pmatrix}^{2}, (104)

where the right-hand side is the Wigner 3j3j-symbol. Since this particular 3j3j-symbol has vanishing magnetic quantum numbers, if 1+2+3\ell\equiv\ell_{1}+\ell_{2}+\ell_{3} is even and together the i\ell_{i} satisfy the triangle inequality, then it can neatly be expressed as

(123000)=(1)/2(/2)!(+1)!i=13(2i)!(/2i)!.\begin{pmatrix}\ell_{1}&\ell_{2}&\ell_{3}\\ 0&0&0\end{pmatrix}=\frac{(-1)^{\ell/2}(\ell/2)!}{\sqrt{(\ell+1)!}}\prod_{i=1}^{3}\frac{\sqrt{(\ell-2\ell_{i})!}}{(\ell/2-\ell_{i})!}~{}. (105)

Otherwise, if \ell is not even or the triangle inequality is violated, the 3j3j-symbol vanishes. Additionally, the following identity

P1P2=3=|12|1+2(123000)2(23+1)P3,P_{\ell_{1}}P_{\ell_{2}}=\sum_{\ell_{3}=\lvert\ell_{1}-\ell_{2}\rvert}^{\ell_{1}+\ell_{2}}\begin{pmatrix}\ell_{1}&\ell_{2}&\ell_{3}\\ 0&0&0\end{pmatrix}^{2}(2\ell_{3}+1)P_{\ell_{3}}~{}, (106)

is particularly useful in the evaluation of the θ\theta-integrals in eq. (103), where for brevity we denoted PP(cθ)P_{\ell}\equiv P_{\ell}(c_{\theta}). Eventually, after applying both eqs. (104) and (106), we find that the θ\theta-integrals in j(1)\mathcal{I}^{(1)}_{j} can be handled with the help of

0πdθsθ3P1P2P3=43(123000)243k=|21|2+1(2k+1)(21k000)2(k23000)2.\int_{0}^{\pi}\text{d}\theta\ s^{3}_{\theta}P_{\ell_{1}}P_{\ell_{2}}P_{\ell_{3}}=\frac{4}{3}\begin{pmatrix}\ell_{1}&\ell_{2}&\ell_{3}\\ 0&0&0\end{pmatrix}^{2}-\frac{4}{3}\sum_{k=\lvert 2-\ell_{1}\rvert}^{2+\ell_{1}}(2k+1)\begin{pmatrix}2&\ell_{1}&k\\ 0&0&0\end{pmatrix}^{2}\begin{pmatrix}k&\ell_{2}&\ell_{3}\\ 0&0&0\end{pmatrix}^{2}. (107)

The integrals j(1)\mathcal{I}_{j}^{(1)} and j(2)\mathcal{I}_{j}^{(2)} may be further simplified by considering the convergence properties of the sum, which allow us to integrate each term in rr separately. Doing so, we rapidly see the usefulness of the previous θ\theta-integral formulae, which appear explicitly as

j(1)=\displaystyle\mathcal{I}_{j}^{(1)}= 1,2,315+νj+2|νj|0πdθsθ3P1P2P3,\displaystyle\sum_{\ell_{1},\ell_{2},\ell_{3}}\frac{1}{5+\ell}\frac{\nu_{j}^{\ell+2}}{|\nu_{j}|^{\ell}}\int_{0}^{\pi}\text{d}\theta\ s^{3}_{\theta}P_{\ell_{1}}P_{\ell_{2}}P_{\ell_{3}}, (108a)
j(2)=\displaystyle\mathcal{I}_{j}^{(2)}= 1,2,32νj20πdθsθ3P1P2P3[Λr2(ϵv,θ)νj2|νj|]\displaystyle\sum_{\begin{subarray}{c}\ell_{1},\ell_{2},\ell_{3}\\ \ell\neq 2\end{subarray}}\frac{\nu_{j}^{\ell}}{2-\ell}\int_{0}^{\pi}\text{d}\theta\ s^{3}_{\theta}P_{\ell_{1}}P_{\ell_{2}}P_{\ell_{3}}\left[\Lambda_{r}^{2-\ell}(\epsilon_{v},\theta)-\frac{\nu_{j}^{2}}{|\nu_{j}|^{\ell}}\right]
+1,2,3=2νj20πdθsθ3P1P2P3ln(Λr(ϵv,θ)|νj|),\displaystyle+\sum_{\begin{subarray}{c}\ell_{1},\ell_{2},\ell_{3}\\ \ell=2\end{subarray}}\nu_{j}^{2}\int_{0}^{\pi}\text{d}\theta\ s^{3}_{\theta}P_{\ell_{1}}P_{\ell_{2}}P_{\ell_{3}}\ln\left(\frac{\Lambda_{r}(\epsilon_{v},\theta)}{\lvert\nu_{j}\rvert}\right)\,, (108b)

where we have isolated the =2\ell=2 mode in j(2)\mathcal{I}^{(2)}_{j} in order to handle the potential log divergence as a separate case.

Starting with the sum on the second line of eq. (108b), we first expand the integral in small ϵv\epsilon_{v} using eq. (101). Using the integral formulae for the Legendre polynomials above, we find that the leading ln(2m^1/ϵv2)\ln\left(2\hat{m}_{1}/\epsilon_{v}^{2}\right) divergence contains no additional θ\theta-dependence, and so due to the constraint that =2\ell=2, is weighted with P12+P2P_{1}^{2}+P_{2}, and vanishes upon integration. Thus, we find that

j=12n+2(1)jξ^j1,2,3=2νj20πdθsθ3P1P2P3ln(Λr(ϵv,θ)|νj|)=𝒪(ϵv4).\displaystyle\sum_{j=1}^{2n+2}(-1)^{j}\hat{\xi}_{j}\sum_{\begin{subarray}{c}\ell_{1},\ell_{2},\ell_{3}\\ \ell=2\end{subarray}}\nu_{j}^{2}\int_{0}^{\pi}\text{d}\theta\ s^{3}_{\theta}P_{\ell_{1}}P_{\ell_{2}}P_{\ell_{3}}\ln\left(\frac{\Lambda_{r}(\epsilon_{v},\theta)}{\lvert\nu_{j}\rvert}\right)=\mathcal{O}(\epsilon_{v}^{4})~{}. (109)

Hence, the only meaningful contributions to A[ζRT]A[\zeta_{\text{RT}}] from j(2)\mathcal{I}^{(2)}_{j} come from the first line in eq. (108b).

Moving on to the cutoff-dependent integrand on the first line of eq. (108b), we can again utilize the small ϵv\epsilon_{v} expansion in eq. (101). Since the leading divergence in Λr(ϵv,θ)\Lambda_{r}(\epsilon_{v},\theta) is 𝒪(1/ϵv2)\mathcal{O}(1/\epsilon_{v}^{2}), we can neglect any integral for >2\ell>2 as it will vanish as 𝒪(ϵv2)\mathcal{O}(\epsilon_{v}^{2}). Truncating to the sum to <2\ell<2, expanding in small ϵv\epsilon_{v}, and evaluating the sum over jj, we find that the total contribution to A[ζRT]A[\zeta_{\text{RT}}] from the first sum in eq. (108b) is

j=12n+2(1)jξ^j1,2,32νj20πdθsθ3P1P2P3Λr2(ϵv,θ)=8m^133ϵv4+2n^125m^1+𝒪(ϵv2).\displaystyle\sum_{j=1}^{2n+2}(-1)^{j}\hat{\xi}_{j}\sum_{\begin{subarray}{c}\ell_{1},\ell_{2},\ell_{3}\\ \ell\neq 2\end{subarray}}\frac{\nu_{j}^{\ell}}{2-\ell}\int_{0}^{\pi}\text{d}\theta\ s^{3}_{\theta}P_{\ell_{1}}P_{\ell_{2}}P_{\ell_{3}}\Lambda_{r}^{2-\ell}(\epsilon_{v},\theta)=\frac{8\hat{m}_{1}^{3}}{3\epsilon_{v}^{4}}+\frac{2\hat{n}_{1}^{2}}{5\hat{m}_{1}}+\mathcal{O}(\epsilon_{v}^{2}). (110)

Finally, we treat the remaining sums in eq. (108a) and the second term on the first line of eq. (108b) together. Firstly, we note that the integral at =2\ell=2 in j(1)\mathcal{I}_{j}^{(1)} vanishes due the integrand being of the form P12+P2P_{1}^{2}+P_{2}. Secondly, we make use of eq. (107) explicitly and observe that only the =0\ell=0 term contributes. That is, if we decompose the sum as partial sums in \ell,

1,2,32(15+12)0πdθsθ3P1P2P3\displaystyle\sum_{\begin{subarray}{c}\ell_{1},\ell_{2},\ell_{3}\\ \ell\neq 2\end{subarray}}\left(\frac{1}{5+\ell}-\frac{1}{2-\ell}\right)\int_{0}^{\pi}\text{d}\theta\ s^{3}_{\theta}P_{\ell_{1}}P_{\ell_{2}}P_{\ell_{3}} =a=0a2(15+a12a)1,2,3=a0πdθsθ3P1P2P3\displaystyle=\sum_{\begin{subarray}{c}a=0\\ a\neq 2\end{subarray}}^{\infty}\left(\frac{1}{5+a}-\frac{1}{2-a}\right)\sum_{\begin{subarray}{c}\ell_{1},\ell_{2},\ell_{3}\\ \ell=a\end{subarray}}\int_{0}^{\pi}\text{d}\theta\ s^{3}_{\theta}P_{\ell_{1}}P_{\ell_{2}}P_{\ell_{3}}
=25561,2,3=10πdθsθ3P1P2P3+,\displaystyle=-\frac{2}{5}-\frac{5}{6}\sum_{\begin{subarray}{c}\ell_{1},\ell_{2},\ell_{3}\\ \ell=1\end{subarray}}\int_{0}^{\pi}\text{d}\theta\ s^{3}_{\theta}P_{\ell_{1}}P_{\ell_{2}}P_{\ell_{3}}+\ldots, (111)

we find that all the partial sums with >0\ell>0 vanish and 2/5-2/5 is the exact result.

Putting all of the results above together and taking the sum over jj, we find

j=12n+2(1)jξ^j(j(1)+j(2))=8m^133ϵv4+25n^12n^2m^1m^1+𝒪(ϵv2).\displaystyle\sum_{j=1}^{2n+2}(-1)^{j}\hat{\xi}_{j}(\mathcal{I}_{j}^{(1)}+\mathcal{I}_{j}^{(2)})=\frac{8\hat{m}_{1}^{3}}{3\epsilon_{v}^{4}}+\frac{2}{5}\frac{\hat{n}_{1}^{2}-\hat{n}_{2}\hat{m}_{1}}{\hat{m}_{1}}+\mathcal{O}(\epsilon_{v}^{2}). (112)

Plugging in to eq. (102), the unregulated area of the RT hypersurface is

A[ζRT]=π4L𝕊49log(2Rϵu)[25631ϵv4+645n^12m^14645n^2m^13+𝒪(ϵv2)]+𝒪(ϵu2).\displaystyle A[\zeta_{\text{RT}}]=\pi^{4}L_{\mathbb{S}^{4}}^{9}\log\left(\frac{2R}{\epsilon_{u}}\right)\left[\frac{256}{3}\frac{1}{\epsilon_{v}^{4}}+\frac{64}{5}\frac{\hat{n}_{1}^{2}}{\hat{m}_{1}^{4}}-\frac{64}{5}\frac{\hat{n}_{2}}{\hat{m}_{1}^{3}}+\mathcal{O}(\epsilon_{v}^{2})\right]+\mathcal{O}(\epsilon_{u}^{2}). (113)

It is then straightforward to see that SEES_{\rm EE} is given by eq. (83).

References