Surface operators in the 6d $\mathcal{N} = (2,0)$ theory

The 6d $\mathcal{N}=(2,0)$ theory has natural surface operator observables, which are akin in many ways to Wilson loops in gauge theories. We propose a definition of a"locally BPS"surface operator and study its conformal anomalies, the analog of the conformal dimension of local operators. We study the abelian theory and the holographic dual of the large $N$ theory refining previously used techniques. Introducing non-constant couplings to the scalar fields allows for an extra anomaly coefficient, which we find in both cases to be related to one of the geometrical anomaly coefficients, suggesting a general relation due to supersymmetry. We also comment on surfaces with conical singularities.


Introduction
Understanding the six dimensional N = (2, 0) superconformal field theory is one of the most intriguing problems in theoretical physics. In this paper we revisit the most natural observables in this theory, surface operators [1]. If we define the theory as arising from N coincident M5-branes, the simplest surface operators correspond to the endpoints of M2branes [2].
In some ways the surface operators in six dimensions are analogous to Wilson loops in lower dimensional gauge theories. Wilson loops are the boundaries of fundamental strings, which are the dimensional reduction of M2-branes, and indeed one obtains Wilson loops in compactifications of the 6d theory with surface operators. Wilson loops are not only interesting due to their physical importance, they are also accessible to many perturbative and non-perturbative calculational tools in supersymmetric field theories: Feynman diagrams, holographic descriptions [3][4][5], localization [6], the defect CFT framework and associated OPE techniques [7,8], integrability [9,10], duality to scattering amplitudes [11] and more. See for instance a recent survey of these techniques, as applied to supersymmetric Wilson loops in ABJM theory [12].
We do not expect all these techniques to extend to surface operators in six dimensions, but it is worthwhile to examine which of them may work, and we hope that some calculations may lead to exact results applicable for all N . Here we take the first step in such an examination, defining the notion of a "locally BPS surface operator" and studying basic properties of their anomalies. This is mainly based on previous work [13][14][15][16][17], which we modify and refine in several ways.
As reviewed in the next section, the evaluation of generic surface operators leads to logarithmic divergences. The anomaly depends on the geometry of the surface, as well as intrinsic properties of the operator which are captured by three numbers, known as anomaly coefficients [18].
The "locally BPS" operator couples to the scalar fields via a unit 5-vector n i . This can be viewed as a coupling to an R-symmetry background, and for non-constant n i we find a new anomaly, proportional to (∂n) 2 , with its own anomaly coefficient.
We perform explicit calculations of the three geometrical and one background coefficients in both the free theory at N = 1 and the holographic description valid at large N . An examination of our results reveals that the new anomaly coefficient matches (up to a sign) one of the geometric ones in both regimes. We present here a simple argument, relying on supersymmetry, why we expect this relation to hold for all N . A more rigorous proof of this relation based on the application of defect CFT techniques to surface operators will be presented in [19].
Beyond the study of N = (2, 0) superconformal symmetry, surface operators in conformal field theories have drawn interest within a number of different contexts. Recent work on entangling surfaces in 4d [20][21][22][23] and theories with boundaries [24][25][26] uses some techniques which apply in our case as well. In particular, the classification of local conformal invariants of surfaces is independent of the codimension and translates to the 6d case [27].
Surface operators in the N = (2, 0) theory have been studied both from a field theory perspective [14][15][16][17] and using holography [28,13]. Corresponding soliton solutions of the M5-brane equations of motion have been discussed in the literature under the moniker of self-dual strings [1].
The resemblance to Wilson loops is evident in both the field theoretic and the holographic approach. In the former, for N = 1 as is studied in Section 3, we define the surface operator in analogy to the Maldacena-Wilson loops [4] as where B + is the pullback of the chiral 2-form to the surface Σ and Φ i are the scalar fields. Since for N > 1 there is no realisation of the theory in terms of fundamental fields, we cannot give an analogous definition of the surface operator. However, by analogy with Wilson loops [3][4][5], in the large N limit, these operators in the fundamental representation have a nice holographic dual as M2-branes ending on the surface and extending into the AdS 7 × S 4 bulk, as discussed in Section 4. In the absence of a scalar coupling breaking the so(5) R-symmetry, these would be delocalised on the S 4 [29,30]. At leading order, we need only consider minimal 3-volumes [4,13] (similar to the minimal surfaces of interest in the Wilson loop case [3][4][5]), and to find the anomaly, which is a local quantity, it is enough to understand the volume close to the AdS boundary. High-rank (anti-)symmetric representations are dual to configurations involving M5 branes shrinking to the surface on the boundary of AdS 7 and have been considered in [31][32][33][34].
The definition in (1.1) includes BPS operators. Simple examples are the plane or sphere with constant unit n i . Other examples are briefly discussed in Section 3 and will be explored in more detail elsewhere [35]. We call operators with generic Σ and unit length n i "locally BPS", and show that they possess some nice properties, in particular that all power law divergences cancel.
In the next section we recall the structure of surface operator anomalies and introduce the anomaly coefficients. We evaluate these anomaly coefficients for the two known realisations of the N = (2, 0) theory; first as the theory of a single M5-brane (N = 1) [36], for which the equations of motion are known [37], and second, using holography (for the large N limit) from M-theory on the AdS 7 × S 4 background [38] found in [39]. The resulting anomaly coefficients are presented in equations (3.24) and (4.18). After performing the free field and holographic calculations, we address in Section 5 surfaces with singularities. We discuss our results in Section 6 and offer a simple argument for the relation between two of the anomaly coefficients. We collect some technical tools in appendices. Our conventions can be found in Appendix A. Details of the geometry of submanifolds are compiled in Appendix B. Appendix C contains an alternative, more geometric derivation of the field theory results in Section 3.

Surface anomalies
The most natural quantities associated to surface operators in conformal field theories are their anomaly coefficients. To understand their origin, note that, unlike line operators, the expectation values of surface operators typically suffer from ultraviolet divergences, which cannot be removed by the addition of local counterterms. The regularised expectation value satisfies where is a regulator, A Σ is known as the anomaly density, and we suppressed possible power-law divergences. A Σ is scheme independent and indicates an anomalous Weyl symmetry, since for a constant rescaling g → e 2ω g, the expectation value varies as where the subscript • g denotes the background metric. The anomaly is constrained by the Wess-Zumino consistency condition [18,40] to be conformally invariant. In dimensions d ≥ 3, the local geometric conformal invariants for a 2d submanifold, which have been classified in [27], are R Σ : The Ricci scalar of the induced metric h ab on Σ.
H 2 + 4 tr P : H µ is the mean curvature, P ab the pullback of the Schouten tensor (B.2).
tr W : W abcd is the pullback of the Weyl tensor.
Under conformal transformations, the first two change by a total derivative (type A anomalies) and the last is itself conformally invariant (type B).
As we allow for variable couplings to the scalars, parametrised by a unit 5-vector n i , we find an extra potential type B Weyl anomaly associated to it: (∂n) 2 ≡ ∂ a n i ∂ a n i . This is (up to total derivatives) the only quantity of the correct dimension that can be constructed using only n.
The anomaly of a surface operator in any 6d N = (2, 0) theory then takes the form 3) The anomaly coefficients a 1 , a 2 , b and c depend on the theory (that is on N ) and the type of surface operator (which, at least at large N , is specified by the representation of the A N −1 algebra [41,42]), but not on its geometry or n. They are the focus of this paper. Let us mention that there exists another commonly used basis where whereĨI µ ab is the traceless part of the second fundamental form (see (B.8)). These bases are related through the Gauss-Codazzi equation (B.7). The relation between the coefficients is then Some results about these anomaly coefficients are known for surface defects in generic CFTs. The bound b 1 < 0 was derived in [22] by showing that b 1 captures the 2-point function of the displacement operator, which is positive by unitarity. Similarly, it was shown in [43,22] that b 2 is calculated by the one-point function of the stress tensor in the presence of the surface defect (this was also conjectured in [44]). Assuming that the average null energy condition holds in the presence of defects also leads to a bound b 2 > 0 [23].
For the surface operators at hand, these anomaly coefficients were also calculated previously. At large N , the first such result was a calculation of the 1/2-BPS sphere [28], with total anomaly −4N , implying a = b (N ) = 0 [13]. More recently, it was conjectured that N = (2, 0) supersymmetry imposes b = 0 (or b 1 = −b 2 ) for any N [45]. a and b 2 were calculated at any N > 1 (and for any representation) by studying the holographic entanglement entropy in the presence of surface operators [46,34,23,47]. This result is also supported by a recent calculation based on the superconformal index [48], which suggests that it is exact.
The anomaly coefficient c has previously not been discussed, to our knowledge.

Abelian theory with N = 1
In this section we study the anomaly coefficients of the surface operator in the abelian (2, 0) theory. This is the theory of a single M5-brane and the degrees of freedom form the tensor supermultiplet of the osp(8 * |4) symmetry algebra. It consists of three fields [49] (see also [50,51]) • A real closed self-dual 3-form H = dB + .
• Five real scalar fields Φ i .

Surface operators and BPS condition
We define the surface operators V Σ of the abelian theory as in (1.1). To avoid complications arising for null surfaces (which could be interesting, but lie beyond the scope of this work), we restrict to space-like surfaces in flat 6d Minkowski space (with mostly positive signature).
A surface operator is BPS provided that its variation under the supersymmetry transformations (3.1) vanishes Since this is an integral over the insertion of an operator ψ along the surface, this is satisfied only when the integrand vanishes at every point along the surface, leading to the projector equation If we impose that n 2 ≡ n i n i = 1, then Π − is a half rank projector and otherwise it is a full rank matrix. In the case of a planar surface with constant unit n i , this is a single condition, so the surface preserves 16 supercharges, i.e. is 1/2-BPS. 1 In analogy to Wilson loops in 4d theories, it is natural to discuss "locally BPS operators" [5], where the equations (3.3) are satisfied at every point along the surface, but without a global solution. This amounts to the requirement n 2 = 1, and as shown below, leads to the cancellation of all power-like divergences in the evaluation of the surface operator.
One can also look for surfaces, other than planes, that preserve some smaller fraction of the supersymmetry by relating n i (σ) to x µ (σ) and its derivatives. One simple way to realise 1 The BPS condition for a surface operator extended in the time-like direction can be obtained by Wickrotation to this is for surfaces with the geometry R × S, for some curve S ⊂ R 1,4 . Upon dimensional reduction this becomes a Wilson loop in 5d maximally supersymmetric Yang-Mills (or 4d upon further dimension reduction). Then one can choose n i to follow the construction of globally BPS Wilson loops of [52] or [53] to find globally BPS surface operators. Indeed this was realised recently in [54] (see also [55]).
There are further examples of globally BPS surface operators, which do not follow this construction. The simplest is the spherical surface, but there are several other classes of such operators, which will be explored elsewhere [35].

Propagators
Since the abelian theory is non-interacting, the expectation value of V Σ reduces to where h is the determinant of the induced metric on Σ. Evaluating this requires expressions for the propagators of the tensor and scalar fields. While one would preferably derive the propagators from an action, none is readily available. Many actions for the abelian N = (2, 0) theory have been proposed over the years, but they all suffer from some pathologies regarding the self-dual 2-form (see [56,57,36,58,59] for examples of available actions, and [60][61][62] and references therein for recent accounts of the various approaches in the abelian theory). In any case, gauge fixing and inverting the kinetic operator is not straightforward.

Tensor structure
We sidestep these obstacles by determining the propagators in other ways. The scalar propagator in flat 6d is fixed by conformal symmetry to be The proportionality constant depends on the normalisation of the fields. It could be determined from an action, but in its absence it is fixed by supersymmetry below. The more complicated question is the self-dual 2-form propagator. Let us start by considering an unconstrained 2-form field B with a free Maxwell type action were α is a gauge fixing parameter. In Feynman gauge α = 1, this gives the propagator Now we decompose the field into its self-dual and anti-self-dual parts B µν = B + µν + B − µν and try to deduce the propagators for each component.
Since there is no covariant 4-tensor satisfying the self-duality properties of a mixed correlator B + B − , we can decompose The two terms on the right hand side need not be identical, but the difference between them should be parity-odd. 2 The only such term of the right scaling dimension which we can write down is However, terms of this type do not contribute to (3.4), since the integration is symmetric in x and y. Therefore, for the purpose of our calculation we can take B + B + = BB /2. Note that in curved space we can add to the right hand side a term proportional to the Weyl tensor with all the required symmetries.

Normalisation
The normalisation of the tensor field propagator is fixed by the assumption that the surface operator defined in (1.1) corresponds to a single unit of quantised charge. First, for any closed surface Σ, we can rewrite the surface operator (without scalars) in terms of the field strength as where ∂V = Σ. In order for this to be well-defined, any two such V with the same boundary must yield the same result. Equivalently, for every closed and similarly for * H. 2 The two-dimensional analogue is instructive. The propagator of a free boson in complex coordinate z is given by while for a (anti-)chiral boson one finds Indeed the sum reproduces the free boson propagator, but the two differ by a parity-violating imaginary part.
Now consider a flat surface operator in the (x 1 , x 2 ) plane, which we view as a source for the unconstrained B field. The solution to the equations of motion would be given by convoluting the propagator with this source. Using the expression in (3.7), we get Again, because we don't know the self-dual propagator, the field strength we obtain is not self-dual, but the quantisation condition should still be satisfied. Imposing that the charge enclosed in a transverse sphere is quantised leads to The normalisation of the scalar propagator is then fixed by supersymmetry. A simple way to implement that is to compare with the classical BPS solution of the self-dual string [1] which gives 3 We emphasise that this normalisation is obtained by imposing a quantisation condition on the unconstrained B-field, where we treat B + as the self-dual subsector of a general 2form. This follows the discussion in [56], however some caution is warranted. The interplay between the quantisation and self-duality conditions could lead to obstructions, resulting in halving the self-dual source on the right hand side of (3.12). In that case, the overall normalisation of both propagators would increase by a factor of two.
With the flat space propagators we are able to determine the anomaly coefficients a 1 , a 2 , and c. The calculation of b, however, requires a the curved space propagator, where the right-hand side of (3.9) could pick up contributions whose integral does not vanish. Since we do not know how to fix these terms, we cannot determine b.
Note though that we can calculate the contribution of the scalars to the anomaly coefficient b. The propagator of a conformal scalar in a curved background can be expanded in powers of the geodesic distance [17], and the contribution to the anomaly coefficient b is read off as −2/3.
If we give up the requirement of self-duality, we can use the short-distance expansion of an unconstrained 2-form propagator on curved space, which has been computed in [14,16], and again, the Weyl tensor of the background explicitly contributes to the curvature corrections.
Halving that to try to account for self-duality and adding to it the contribution from the scalars, one obtains b = −4/3 [17]. This is in disagreement with the conjecture b = 0 [45] and therefore one may not trust it.

Evaluation of the anomaly
With the propagators at hand, we can compute the expectation value of the surface operator by evaluating the integrals in (3.4). Generically, these integrals are divergent and must be regularised.
In this section we take a rather naive approach of placing a hard UV cutoff on the double integral (3.4), so as to restrict |σ − τ | > (where the distance is measured with the induced metric), the same regularisation that is used in [14]. A different regularisation is employed in [17], where the surface is assumed to be contained within a 5d linear subspace of R 6 and the two copies of the surface are displaced by a distance in the 6th direction. This restriction to R 5 must still yield the correct answer, since even for surfaces in 4d the geometric invariants in the anomaly (2.3) are independent of each other. Still, in Appendix C we redo the calculation removing this assumption by displacing the two copies of the surface along geodesics in the direction of an arbitrary normal vector field. That approach could be important for the calculation of surface operators in four dimensions, where the restriction to a 3d linear subspace does not allow to resolve all the anomaly coefficients.
To find the anomalies we only need the short-distance behaviour of the propagators, so we use normal coordinates η a about a point σ on Σ. The notations and required geometry are presented in Appendix B.
Starting from the scalar contribution to (3.4), the integrand is Using n i n i = 1 and (B.14), (B.12) we have The integral computing the density of the scalar contribution to log V Σ is then (3.17) Using polar coordinates η a = η e a (ϕ), where e is a 2d unit vector, and the identities δ ab δ cd + δ ac δ bd + δ ad δ bc , (3.18) we are left with the radial integral, for which we introduce the cutoff To get the expression in the second line we also used the Gauss-Codazzi equation (B.7). The calculation of the contribution of the 2-form field is very similar. Expanding the tensor structure, we have (3.20) In terms of η a , the differential forms read (see (B.9)) (3.21) Collecting terms and introducing a radial cutoff as above, we find the contribution As discussed above, since we do not know the contribution of the Weyl tensor to the B-field propagator, we cannot determine b (1) . According to the conjecture of [45] however, it should vanish. This relation is the subject of work in progress [19]. Equation (3.23) differs from (2.3) by the absence of the tr P term, which also vanishes in flat space. Since H 2 doesn't vanish in flat space, it determines a 2 unambiguously and in curved space H 2 is necessarily accompanied by 4 tr P , based on the general argument for the form of the anomaly reviewed in Section 2.
Finally, we reiterate that, depending on the form of the quantisation condition, the result for the anomaly coefficients may be multiplied by 2, see the discussion following (3.14). In any case, the abelian theory should have surface operators with an integer multiple of iB + − n i Φ i in (1.1), and for all of them it is still true that a

Generalising the scalar coupling
Note that the preceding calculation is applicable regardless of whether the operator is locally BPS or not, so we may relax the condition n 2 = 1. In that case the result for the anomaly coefficients is a If we replace n i → in i , we recover the expressions for the surface operator studied in [17]. An operator with n 2 = 0 was studied in [14], but assuming a non-self-dual 2-form. The anomaly coefficients computed in [14] are half of the ones we obtain by setting n 2 = 0 in (3.25), due to a difference in the overall normalisation of the propagator. It would be interesting to study this system in the large n 2 limit. This is similar to the "ladder" limit of the cusped Wilson loop in N = 4 SYM in 4d first suggested in [63] which is related to a special scaling limit of that theory, dubbed the "fishnet" model, which also has a 6d version [64].

Holographic description at large N
The holographic calculation of the Weyl anomaly for surface operators was pioneered by Graham and Witten in [13]. Here we present a rewriting of their argument, which we also generalise slightly to include operators extended on the S 4 .

Surface operators
The N = (2, 0) theory is described at large N by 11d supergravity on an asymptotically AdS 7 × S 4 geometry [38] ds 2 = L 2 y 2 dy 2 + g (0) + g (1) y 2 + such that g (0) is the metric of the dual field theory 4 and g S 4 is the metric of S 4 . The background also includes N units of The full form of the metric is determined by the supergravity equations of motion in the presence of fluxes and by requiring the geometry to close smoothly in the interior. While the latter requires nonlocal information, the near-boundary expansion is fixed to the required order by local information about the boundary. Following [65,66], the first term in this expansion was found in [13] as At this order the S 4 is round, so to leading order the solution to (4.2) is simply The holographic description of the surface operators (1.1) is by M2-branes anchored along Σ on the boundary of AdS [4]. UsingΣ for the world-volume of the M2-brane, it has a boundary at y = 0 with ∂Σ = Σ. The expectation value of the surface operators is then given by the minimum of the M2-brane action, reading (in Euclidean signature and with all fermionic terms suppressed) [67] log where T M2 is the tension of the brane, proportional to N . volΣ is the volume form calculated from the induced metric and A 3 is the pullback of the 3-form potential.

Local supersymmetry
Before studying the M2-brane embeddings, let us note that the M2-brane minimizing (4.5) is also locally supersymmetric. The supergravity fields appearing there sit in the supergravity multiplet, which transform as where EM M , Ψ M and A 3 are respectively the vielbein, gravitino and 3-form potential of F 4 (M = 1, . . . , 11 is the frame index). Using these transformations, the variation of (4.5) is We here denote the coordinates on the world-volume byσâ. The projector equation is then The projector is again half-rank, so that the M2-brane locally preserves half of the supersymmetries (16 supercharges). These supercharges can be shown to agree with the field theory BPS condition (3.3) on Σ once we decompose x M into coordinates on the boundary of AdS, x µ , and the S 4 coordinates n i .

Holographic calculation
To find the saddle points of the action (4.5), we parametrise the M2-brane by y, σ a where σ a are coordinates for Σ. We then use the static gauge to describe the embedding by {u a (y, σ), n i (y, σ)}, where u a are the normal directions to the surface Σ at y = 0. In this setup, the boundary conditions are u a (y = 0, σ) = 0 and n i (y = 0, σ) = n i (σ) (where the right hand side has the n i from (1.1)). Because the metric (4.1) diverges at the boundary of AdS, the volume element on the M2-brane diverges as y −3 , which leads to divergences in the action. Finding the shape of the embedding requires knowledge of the full surface and is generally a hard problem. But since we are only interested in the logarithmically divergent part of the action, it is sufficient to solve the equations of motion for small y. We do this perturbatively following [13], mirroring the solution of the background supergravity equations above.
Using (4.3), the lowest order terms in the metric for our coordinates normal and tangent to the surface, are Here h ab = g (0) ab u=0 is the metric on Σ. Note that away from y = 0, this metric depends on u a (for y = 0, generically u a = 0), as in the first line.
To write down the M2-brane action we need the induced metricĥ ab = ∂ a X M ∂ b X N g M N (including also the S 4 directions). We expand the embedding coordinates as u a (y, σ) = O(y 2 ) , n i (y, σ) = n i (σ) + O(y 2 ) . (4.10) It is easy to check that higher order terms are not required. Then the S 4 metric can be replaced with g (0) S 4 = δ ij dn i dn j and the second fundamental form is II a ab = − 1 2 g a b ∂ b g ab . Dropping the explicit O(y ) as well as the subscript | u=0 along with the superscript (0) , since all the quantities are evaluated on the surface, we find (4.11) The determinant of the metric is then detĥ L 6 y 6 1 + ∂ y u a ∂ y u b g a b − 2H a u b g a b + − tr P + 1 4 (∂n) 2 y 2 det h , (4.12) while the pullback of the 3-form does not contribute to the divergences. We thus find the action At order O(y 2 ), we need only solve for u a (y), which has the equation of motion The action evaluated at the classical solution is then where we see that the anomaly indeed takes the form (2.3). The result is where we discarded an irrelevant term proportional to −2 (see the discussion below). This result agrees with the original calculation of [13] and adds to it the coupling to (∂n) 2 . It is also consistent with the explicit calculation of the 1/2-BPS sphere [28], for which the anomaly is −4N . The anomaly coefficients at leading order in N are then As in the case of Wilson loops in N = 4 SYM in 4d, we expect this holographic description to be correct in the locally BPS case when the scalar couplings satisfy n 2 = 1. Following [29,30], the case of n 2 = 0 should be described by the same surface inside AdS 7 , but completely smeared over the S 4 . In this case we find the same result for the geometric anomaly coefficients as above, and, since the corresponding anomaly term vanishes, c (N ) does not apply.

Power-law divergence
Note that in addition to the log divergence in (4.17), (4.16) produces also a power-law divergence While such divergences can be removed by the addition of a local counter-terms, in the field theory result (3.23), they cancelled without extra counter-terms (for the locally BPS operator). A more elegant way of eliminating the power law divergences also in this holographic calculation follows the example of the locally BPS Wilson loops [5]. A careful treatment of the boundary conditions suggests that the natural action is a Legendre transform of (4.5), which differs from the action we used by a total derivative. This modification does not change the equations of motion, but gives a contribution on the boundary, where it precisely cancels the divergence above.
By looking at the M5-brane metric before the decoupling limit, we can identify the coordinate to use in the transform as r i = L 3 n i /2y 2 . Defining its conjugate momentum by differentiating with respect to the boundary value of the coordinate (where y = ) . (4.20) In the last equality we used the value of the classical action (4.16), undoing the integration, so the classical Lagrangian density. The Legendre transformed action is theñ The last term exactly cancels the power law divergence in (4.19).

Surfaces with singularities
An interesting class of surface operators that has received some attention recently is surfaces with conical singularities. For these surfaces, it was found that the regularised expectation value typically diverges as [68,21,69,70] log V Σc ∼ A log 2 + O(log ) .
Let us consider a conical defect (on flat space) of the form We allow here also a "conical singularity" in the scalar couplings, which has s dependence even as r → 0. It is possible to also allow x µ and n i to have higher order terms in r, but since those lead to subleading divergences, they are unimportant. We can try to use the usual formula for the anomaly (2.3) by plugging in the geometric invariants R Σ = Ωδ(r) , where Ω is the deficit angle, κ =γ 2 /|γ| 2 is the curvature of γ. Plugging into (2.3), the Ricci scalar gives a finite contribution, but H 2 and (∂n) 2 diverge as r → 0. Introducing a cutoffˆ on the r integration, this gives This expression is a bit naive, as we should treat all divergences on the same footing and should identifyˆ = . But then we should not use (2.3), rather go back one step and regularise the divergences that gave rise to the original log divergence while also applying it to the r integration. As we show below, this leads to the expression in (5.4) with log logˆ → 1 2 log 2 . In both the free field case and the holographic realisation this factor of 1/2 is a simple consequence of the usual coefficient of the quadratic term in the Taylor expansion, or in other words of an integral of the form log r d log r.
This factor of 1/2 was noticed already in the calculations of [68,21] and justified in [70] by a careful treatment of the holographic calculation, which is repeated below. We think that the comparison of this to the free-field calculation and the universal nature of our result further elucidates this mismatch from the naive expectation. Our calculation is also more generic, for allowing arbitrary conical singularities and incorporating the scalar singularities too.
We should note, as already noticed in [21], that surfaces with "creases", i.e. co-dimension one singularities, do not lead to additional log 2 divergences and the expression (2.3) can be immediately applied to them.

Field theory
Here we do not rely on (3.23), but go further back to where the log arises from an integral of the form (3.19) ρ dη η = − log + finite, (5.5) where η is a radial coordinate around the point x, and ρ is an IR cutoff related to the overall size of the surface, or at least a large smooth patch where we defined our local coordinate. Near the cone the smooth patch is bounded by the distance from x to the apex, which we denote by r. The integral instead gives With this careful treatment of the log, we can go back to (2.3), plug in the expressions from (5.3) and integrate over r and with the same UV cutoff to find (5.7)

Holography
The derivation in holography is similar. We first note that conformal symmetry fixes the form of the solution as y(r, s) = ru(s) (5.8) To get to (4.17), we integrate over y, but the conformal ansatz suggests to impose the range ≤ y ≤ ru max . Plugging the curvatures from (5.3) into equation (4.17) we arrive at which again gives the log 2 divergence with the same 1/2 prefactor, as in the field theory (5.7).

Example: circular cone
As a simple example of a singular surface we compute explicitly the anomaly of a cone. Denoting the deficit angle by φ (see figure 1) and including an internal angle θ for the scalar coupling n i , we parametrise the cone as follows The conformal invariants are explicitly The divergence is then Notice that as long as the anomaly coefficients satisfy the relation a 2 = −c, which we have shown to hold in the abelian and large N case, the anomaly vanishes for configurations θ = ±φ, which correspond generically to 1/8-BPS configurations.
φ r s θ s Figure 1: On the left, the surface wraps a (circular) cone with a deficit angle φ. On the right, the scalar coupling follows a circle at angle θ on S 2 . For a fixed r, we have a curve that simultaneously traces the circles γ(s) and n i (s).

Conclusion
In this paper we calculated the anomaly coefficients of locally supersymmetric surface operators in the N = (2, 0) theory in 6d, refining and generalising the calculations of [13,17]. We first introduced a new anomaly coefficient c (2.3) arising from non-constant dependence on the internal R-symmetry directions. These are explicit scalar couplings in the abelian theory and motion on S 4 in the holographic realisation. We then presented an explicit calculation for the abelian theory and for the large N limit (using holography). The results are in equations (3.24) and (4.18). Although we are not able to compute the anomaly coefficient b at N = 1 because we do not know the general curved space propagator for the self-dual 2-form, we found the others in both cases.
Making all N conjectures based on the asymptotics is a fool's errand, which we carefully tread. This is especially true given that the abelian theory is not the same as the A N −1 theory at N = 1, since the latter is the empty theory. Nevertheless, in both cases we see that a 2 = −c, and we expect this to hold generally. The argument is based on the BPS Wilson loops of [52], where n i is parallel toẋ µ and which have trivial expectation values. If we uplift them to the 6d theory we expect to find surface operators with no anomaly (and no finite part as well). These operators satisfy H 2 = (∂n) 2 and indeed they do not contribute to the anomaly 5 for a 2 = −c. A proof of this relation as well as properties of b, based on defect CFT techniques, will be presented elsewhere [19].
Two more results are the formalism for regularising surface operators presented in Appendix C and the expression for the divergences due to conical singularities over arbitrary curves in Section 5.
All our calculations are for a surface operator in the fundamental representation. It is expected that 1/2-BPS surface operators are classified by representations of the A N −1 algebra 5 In the uplift we find only surfaces with trivial topology, so the anomaly vanishes regardless of a 1 . of the theory. At large N this is proven, since the asymptotically AdS 7 × S 4 solutions of 11d supergravity preserving the symmetry algebra of 1/2-BPS surface operators can be classified in terms of Young diagrams [41,42].
A calculation of anomalies of surface operators in arbitrary representation, based on the bubbling geometries and holographic entanglement entropy was undertaken in [47]. If we assume b = 0, then for a the fundamental representation, their result reads This is supported by an independent calculation using the superconformal index [48]. In the large N limit, our result [13] indeed agrees with theirs. These calculations do not determine the remaining anomaly coefficients in generic representations. But if we believe the b = 0 conjecture of [45] and our argumentation above for c = −a 2 , this fixes the remaining ones. It would be interesting to reproduce these finite N corrections using other methods as well as do direct holographic calculations for higher-dimensional surface operators.
The anomalies studied here are the most basic properties of surface operators, but finding them is only a first step in understanding these observables and the mysterious theory they belong to. Planar/spherical surface operators preserve part of the conformal group (and with the scalar coupling also half the supersymmetries) and their deformations behave like operators in a defect CFT. A natural next step is to study the defect CFT data: spectrum and structure constants.
Another natural question is the classification of globally BPS surface operators (and local operators within the surface operators) beyond the case of the plane/sphere.
We hope to report progress on these questions in the near future.

A.1 d = 11 Clifford algebra
The 11d Clifford algebra is generated by the set of matrices (Γ M ) B A satisfying Here for readability M is used for flat spacetime, unlike (4.6) where it denotes curved spacetime.
The matrices may be chosen such that Γ † 0 = −Γ 0 is antihermitian while the others are hermitian Γ † M = Γ M (M = 0). In addition, there is an orthogonal, real anti-symmetric matrix C AB such that Γ M C = − (Γ M C) T . C naturally defines a real structure by relating Ψ and Ψ † asΨ This is the Majorana condition.

A.2 d = 6 Clifford algebra
An easy way to construct the 6d Clifford algebra is to decompose Γ M = {Γ µ , Γ i } by introducing a chirality matrix Γ * = Γ 0 Γ 1 Γ 2 Γ 3 Γ 4 Γ 5 . The matrices are then (in the chiral basis) where the algebra is Since γ µ andγ i commute, they define independent spinor representations. Explicitly, we decompose A = (α ⊕ α) ⊗α, so that the indices are (γ µ )β α , (γ µ ) β α and (γ i )β α . The chiral and antichiral representations are related through The chirality operator gives 2 additional constraints . . γ ρ] the antisymmetrised product of γ-matrices. 6 The charge conjugation matrix takes the form and is used to lower (or raise) spinor indices. The matrix Ωαβ is the real, antisymmetric symplectic metric of sp(4) and c is unitary: They satisfy A representation of this algebra is given by (A.12)

A.3 Symplectic Majorana condition
In 6d the spinor Ψ decomposes into a chiral and an antichiral 6d spinor as The Majorana condition on Ψ then translates to where in the second equality we use the properties of our representation. The inclusion of the symplectic form Ω in (A.14) is the reason these equations are known as the symplectic Majorana condition. The spinors ε 1 ,ε 2 , and ψ in (3.1) are of this type.

B Geometry of submanifolds
In this appendix we assemble the geometry results used throughout the main text and in Appendix C. Sections B.1 and B.2 contain our conventions for Riemann curvature and the definition of the second fundamental form of an embedded submanifold as well as some standard results relating the two. In Section B.3 the second fundamental form is related to the coefficients of the normal coordinate expansion of the embedding.

B.1 Riemann curvature
We adopt the convention where the Riemann tensor is defined as It is convenient to split it into a conformally invariant Weyl tensor W µνρσ and the Schouten tensor P µν ,

B.2 Extrinsic curvature
We define the second fundamental form to be The second part is the projector to the components orthogonal to the surface (defined by its embedding x µ (σ)), while the first part is the action of the covariant derivative on the (pullback) of x λ (σ

B.3 Embedding in normal coordinates
Using these standard geometry results, we now derive the expressions needed for (3.16) and (3.21). Unlike in Section 3, we state here the result for a generic curved spacetime M . This allows us to perform the calculation in Appendix C on curved space. Let x µ and η a be Riemann normal coordinates on M and Σ about the same point. In terms of these, the embedding Σ → M may be expanded as These coefficients are constrained by the condition that straight lines in normal coordinates correspond to geodesics. In particular, a curve on Σ given by a straight line in η has constant speed and its curvature in M is normal to Σ at every point, which gives the constraints (B.10) Using (B.4) one easily checks that the second order coefficient equals the second fundamental form The geodesic distance between ξ(η) and the origin of the normal frame is found from (B.9) Furthermore, in normal coordinates, the metrics take the form which yields an expansion for the volume factor (B.14)

C Geodesic point-splitting
In this appendix we present an alternative regularisation of (3.4), essentially point splitting, displacing one copy of the surface operator by a distance in an arbitrary normal direction ν. This regularisation is used in [16,17], but there the vector ν is taken to be a constant, and therefore the method is only applicable if the operators are restricted to a codimension-one subspace.
The technology used to define this regularisation scheme applies for generic smooth embedded surfaces in a Riemannian manifold, and we present here a curved space calculation, as opposed to Section 3.3, where for brevity we restricted ourselves to flat space. However, we still have to restrict to conformally flat backgrounds, since otherwise we do not have a short-distance expansion for the propagator and therefore still cannot infer the anomaly coefficient b.
As expected, we recover the result (3.23) exactly, and thus verify scheme-independence.

C.1 Displacement map
We can regularise the integral (3.4) by displacing a copy of the surface a distance along a unit normal vector field ν. Under that map, which we denote by T , the geodesic distance admits an expansion of the form |T (x µ (σ)) − x µ (σ + η)| 2 = 2 + η 2 + We calculate the higher order terms in (C.1) explicitly in (C.7), but first we note that the only terms contributing to the divergent part are f (3) and f (4) . To see that, the integrals computing the expectation value take the form ρ 0 η m+1 dη |T (x µ (σ)) − x µ (σ + η)| 4 , (C. 3) where ρ is an arbitrary but fixed IR cutoff. We can evaluate (C.3) by expanding the integrand in . Writing s ≡ η/ , we obtain By application of Faà di Bruno's formula one checks that the terms in brackets of order n contribute to the divergence only if m + n ≤ 2. We can therefore safely ignore higher orders in . Only a finite number of terms remains to be computed and we find that the only divergent integrals (C.3) are: 1 + 4f 3 ) 2 + 2f The relevant coefficients can be read off of the expansion of the geodesic distance up to combined order of 4 in η and . The second term on the left hand side of (C.1) can be expanded simply using the embedding (B.9). For the first term, we solve the geodesic equation order by order in the displacement to obtain T (x µ ) = x µ + ν µ − 2 2 Γ µ κλ ν κ ν λ + 3 6 −∂ ν Γ µ ρσ + 2Γ µ νλ Γ λ ρσ ν ν ν ρ ν σ + O( 4 ) . (C.6) Combining these expressions, and writing η a = ηe a (ϕ) as in (3.18) and onwards, the only two non-vanishing relevant coefficients read The first contributes to a scheme-dependent divergence −1 , while the second contributes to the anomaly.

C.2 Evaluation of the anomaly
With the displacement map (C.6) in hand, we can evaluate (3.4). The propagators on a conformally flat background can be obtained by considering curved space actions for a conformal scalar and a Maxwell-type 2-form and inverting the kinetic operators order by order, following [17] and [14]. We find: To apply our regularisation, we should replace ξ by (C.1) in the denominator of the propagators before performing the integral over η. A priori, we should also perform the displacement in the numerator, since a term of order O( ) can contribute to the −1 divergence by multiplying (C.5a). However, one easily checks that the only terms of that order are accompanied by nonzero powers of η, and therefore do not contribute to the divergence of (3.4). We therefore drop the in the numerators of the propagators. The expansion of the numerators is then assembled, as before, from (3.16) and (3.21), but in addition, since we are working on curved space, we obtain an additional term at O(η 2 ) explicitly involving tr P from the propagators (C.8). Collecting terms in analogy to Section 3.3, and integrating out the angular coordinate using (3.18), we obtain the scalar contribution 1 2π 2 + H · ν 4π + 1 16π 2R Σ − H 2 + 4 tr P + 4 (∂n) 2 log + finite, (C.10) while the tensor field yields −2R Σ + 3 H 2 + 4 tr P log + finite. (C.11) Combining these terms, we find log V Σ = 1 4π log Σ vol Σ R Σ − H 2 + 4 tr P + (∂n) 2 + finite, (C.12) which agrees exactly with (3.23). Note that the scheme dependence, which is present in the simple pole of both (C.10) and (C.11), cancels in the final result, and the terms H 2 and tr P combine to an anomaly term as in (2.3), as required.