Supersymmetric Wilson loops via integral forms

We study supersymmetric Wilson loops from a geometrical perspective. To this end, we propose a new formulation of these operators in terms of an integral form associated to the immersion of the loop into a supermanifold. This approach provides a unifying description of Wilson loops preserving different sets of supercharges, and clarifies the flow between them. Moreover, it allows to exploit the powerful techniques of super- differential calculus for investigating their symmetries. As remarkable examples, we discuss supersymmetry and kappa-symmetry invariance.


Introduction
Wilson loops [1] are among the most important physical observables in gauge theories. They describe the phase developed by a charged particle moving in a gauge background and are physically detectable in quantum theory. In the confining phase of QCD they allow to compute the static quark-antiquark potential. They are also one of the basic notions entering a lattice formulation of strongly coupled gauge theories. In topological theories, like pure Chern-Simons theories in three dimensions, their vacuum expectation values provide information on the topological invariants of the model. In the context of the AdS/CFT correspondence some supersymmetric Wilson loops may have a dual description in terms of minimal string worldsheet solutions. Therefore they play an important role in testing the correspondence, especially when they are amenable to exact evaluation via localization.
In supersymmetric gauge theories a distinguished class of Wilson loops are the socalled BPS Wilson loops. These are operators that are invariant under a fraction of the supersymmetry charges. They are in general formulated in components, and are expressed as the holonomy along ordinary contours of generalized connections which contain matter fields in addition to the ordinary gauge connection. Studying their invariance under supersymmetry and classifying them in terms of their BPS degree is not always an easy task. Therefore a manifestly supersymmetric formulation would be desirable.
The supersymmetric generalization of an ordinary Wilson loop appeared for the first time in [2], for four dimensional gauge theories in N = 1 superspace. Roughly speaking, it corresponds to replacing the ordinary path-ordered exponential as W = Tr P e λ dx µ Aµ −→ Tr P e Λ dz M A M (1.1) where z M = (x a , θ α ,θα) are superspace coordinates running on a supercontour Λ and A M = (A a , A α ,Āα) is the gauge superconnection. Further study of these operators has been done later in the development of supersymmetric field theories [3][4][5][6][7]. More recently, this kind of operators have been investigated within the context of the AdS/CFT correspondence [8][9][10], and integrability and Yangian invariance of the N = 4 SYM theory [11][12][13]. Light-like super-Wilson loops have been studied as dual to super-amplitudes in N = 4 SYM [14][15][16][17], also in a twistor formulation [18][19][20].
In this paper we discuss a geometric formulation of supersymmetric Wilson loops alternative to (1.1), which makes use of integral forms in supermanifolds. 1 Our proposal for a super-Wilson loop in superspace is the following where A (1|0) is the gauge superconnection evaluated on the entire (n|m)-dimensional supermanifold SM and Y (n−1|m) Λ is a superform representing the Poincaré dual (PCO) of the immersion of the supercontour Λ into the supermanifold. Since the integration measure on 1 Integral forms have been already used to develop a geometric formulation of some simple topological theories such as super Chern-Simons theory [21][22][23] towards d = 3 N = 1 supergravity.

JHEP04(2020)161
SM is [dxdθd(dx)d(dθ)] (see [24]) with x and θ the bosonic and fermionic coordinates, respectively, this formulation puts on the same ground even/odd coordinates and even/odd differentials in a unified treatment, so it can potentially help to clarify the role of the fermionic part of the supercontour.
Ordinary (bosonic) loops are naturally built-in. In fact, it is sufficient to set the m fermionic coordinates and the corresponding differentials to zero to recover the bosonic counterpart of (1.2). In that case the bosonic PCO Y (n−1|0) λ describes the immersion of an ordinary path λ into the bosonic slice of the supermanifold.
This formulation can be easily generalised to include a given parametrization of the supercontour coordinates in terms of a variable τ ∈ T ⊂ R. In this case the integration is performed on the extended supermanifold SM × T and all the τ dependence is assigned to the PCO, which lifts to an (n|m)-superform. This becomes very useful especially in nonabelian gauge theories, as the problem of dealing with the path-ordering can be totally encoded in the PCO ordering.
A crucial advantage of our formulation regards the study of symmetries and invariances of the Wilson operator. Since Γ in (1.2) is expressed as the integral of an integral form on the entire supermanifold, invariance under (super)diffeomorphisms is manifest. Moreover, being the integrand factorized into the product of two objects, invariance translates into the important identity which relates the variation in form of Γ (and then of the Wilson loop) to the variation of its supercontour. Therefore, the invariances of the Wilson operator are totally ascribable to the isometries of the PCO, which in turn can be investigated using differential geometry and cohomology. As we discuss in the main text, all the PCOs belong to the same d-cohomological class, i.e. the addition of a d-exact term does not change their defining properties. However different representatives, that is PCOs corresponding to different contours, exhibit in general a different spectrum of isometries. This freedom of choosing a particular representative can be used to algebraically impose a given set of isometries on Y (n−1|m) Λ , leading to a Wilson loop that possesses a given set of symmetries. We exploit this mechanism of d-varying symmetries to investigate the behaviour of a super-Wilson loop under supersymmetry and kappa symmetry. A notable example is the BPS Wilson-Maldacena loop in N = 4 SYM that we prove to be obtainable from the ordinary non-BPS operator by the addition of a suitable d-exact term to the original PCO.
As clearly appears from this discussion, in our formulation the problem of classifying BPS Wilson operators translates into a cohomological problem. In particular, the d-equivalence of all PCOs implies that δY (n−1|m) Λ in (1.3) is always a d-exact term. If we restrict to supersymmetry variations, this gives rise to a Killing spinor equation whose solutions allow to classify the whole spectrum of BPS operators with different degree of supersymmetry. We do not discuss this in general, but recover some known examples in four and ten dimensions. In particular, we find that constraints for kappa-symmetry invariance in ten dimensions correspond to BPS constraints in four dimensions [9,12,13].

JHEP04(2020)161
Our formulation of Wilson operators is ready to be generalized to the case of curved (super)manifolds, so leading to Wilson operators in (super)gravity, which technically is already built-in. It is also easily adaptable to the description of higher dimensional objects like Wilson (hyper)surfaces.
The paper is structured as follows. In section 2 we briefly review the main tools of superdifferential calculus, primarily integral forms and Poincaré duals. Section 3 is focused on the geometrical construction of abelian Wilson loops along the lines described above, both for the bosonic and the supersymmetric cases. Within the present geometric framework, in section 4 we investigate Wilson loop invariance under a reparametrization of the path, superdiffeomorphisms, supersymmetry and kappa symmetry. In particular, we show how Killing spinor equations corresponding to BPS Wilson loops arise in the present formalism. In section 5 the generalization to Wilson loops in non-abelian gauge theories is briefly presented. Finally, section 6 contains a brief discussion about the interesting relation between our geometric construction of Wilson operators and a similar construction in the context of pure spinor string theory. A brief summary of our main results and a discussion on possible follows-up can be found in section 7. Five appendices follow, which provide some technical material to support the main text and the equations therein.

Integral forms and Picture Changing Operators
The geometric formulation of supersymmetric Wilson loops that we propose in sections 3 and 5 heavily relies on differential supercalculus. Therefore, we begin by briefly recalling the main concepts that will be used. We refer to appendices A, B and C for more details.
The basic ingredients of differential supercalculus are differential superforms defined on supermanifolds [24]. In the space of differential superforms there is no notion of top form, that is a form that can be suitably integrated on the supermanifold. This is due to the commuting nature of the fundamental one-forms dθ's corresponding to odd θ-coordinates. As proposed in [24][25][26][27][28], the notion of top form has to be found into a new complex of forms known as integral forms. Here we follow the strategy pioneered by Belopolsky [29], where integral forms are distributional-like forms on which a suitable Cartan calculus can be developed [22,24,30,31].
The strategy that we use for constructing integral forms and the corresponding supermanifold integrals is the following. Given a bosonic p-form ω (p|0) on a supermanifold SM of dimensions (n|m) (n ≥ p), its integration over a p-dimensional submanifold N ⊂ SM can be defined as the integration on the entire supermanifold of the integral [29,32]. Precisely, if we denote ω

JHEP04(2020)161
The second expression is the integral over the whole supermanifold of a (n|m)-dimensional top form to which we can then apply the usual Cartan calculus rules. Operator Y (n−p|m) N is also known as Picture Changing Operator (PCO), being related to a similar concept in string theory (see e.g. [27,33,34]). This is a well-known formula in differential geometry (see for example [35]) which allows to disentangle the geometrical properties of the immersed surface N in the entire supermanifold from the properties of the ω (p|0) integrand. In topological field theories it is a powerful tool used to prove the Duistermaat-Heckman formula [36][37][38] for the localization technique and to implement the computations in that framework using the Thom isomorphism [39].
An interesting advantage of prescription (2.1) is that it converts the integration region from N to the entire supermanifold, so making invariance under superdiffeomorphisms manifest.
The PCO in (2.1) is independent of the coordinates, it only depends on the immersion through its homology class. It has two crucial properties. First, it is closed but not exact Second, by changing the immersion ι to an homologically equivalent surface N , the new Poincaré dual Y (n−p|m) N differs from the original one by d-exact terms. It is important to note that if ω (p|0) is a closed form, then (2.1) is automatically invariant under any change of the embedding (we will always assume there are no boundary contributions).
A notable example of application of this formalism is represented by the action of a rigid supersymmetric model, which can be written as where the (n|0)-form lagrangian L (n|0) (Φ, V, ψ) is built using the rheonomic rules (see [40][41][42]) and turns out to be a function of dynamical superfields Φ and the rigid supervielbeins V a , ψ α defined in eq. (B.3). The PCO Y (0|m) instead contains only geometric data, for instance supervielbeins or coordinates themselves. If dL (n|0) (Φ, V, ψ) = 0 we can change the PCO by exact terms without changing the action. This can be conveniently exploited for choosing for instance a PCO that possesses manifest symmetries.
This example has a natural generalization to supergravity. After the change (A.8), (E a , E α ) are promoted to dynamical fields and the action becomes The closure of the lagrangian and the closure of the PCO imply the conventional supergravity constraints that reduce the spectrum of independent fields to the one of physical fields.

JHEP04(2020)161
3 Geometric construction of a supersymmetric Wilson loop: the abelian case We present a general construction of supersymmetric Wilson loops in terms of integral forms. The main goal is to obtain a general expression suitable for any geometry of the loop and whose invariances are easily analysable. For the time being we restrict to the case of an abelian gauge theory. The generalization to the non-abelian case is discussed in section 5.

Ordinary Wilson loops as integral forms
As a warm-up, we begin by discussing how to write ordinary (i.e. bosonic) Wilson loops in terms of integral forms. Given an abelian gauge theory with gauge connection A (1) defined on a manifold M of arbitrary dimension n, a Wilson loop along a curve λ ⊂ M is given by where A (1) * is the pull-back of the connection one-form A (1) = A a dx a along the curve. The integration of a one-form ensures the parametrization independence of the loop. As usual, by choosing a suitable parametrization, one can compute the integral.
When λ is a closed path the W operator is gauge invariant. This can be made manifest by alternatively expressing the Wilson loop in terms of the curvature two-form F (2) = dA (1) . In fact, using the Stokes theorem we can rewrite Γ as an integral over a two dimensional surface S whose boundary is λ This expression is then manifestly invariant under gauge transformations. We now prove that Γ can be rewritten as the integral of an n-form on the entire manifold M. To this end we introduce the PCO dual to the immersion of the one-dimensional curve λ into the manifold M where {φ i } i=1,...,n−1 is a set of (n − 1) functions whose zero locus defines the curve λ ⊂ M. In the second equality we have used δ(dφ i ) = dφ i , being dφ i anticommuting differential one-forms. As a simple example we consider the unit circle in two dimensions. In this case we have a single function φ(x 0 , x 1 ) = x 2 0 + x 2 1 − 1 whose locus defines the curve. The Poincaré dual to the immersion is then and it is manifestly invariant under the O(2) isometry group of the circle. 3

JHEP04(2020)161
The PCO in (3.3) possesses the following fundamental properties where ι j is the contraction along the vector field ∂ j and acts as ι j δ(dφ i ) = ∂/∂(dφ j )δ(dφ i ), while δφ j is the variation of the constraints. The first identity can be proven by using the chain rule dδ(φ i ) = ( j dφ j ∂ ∂φ j )δ(φ i ) (the differential dφ j is kept on the left hand side of the delta) and the distributional property dφ i δ(dφ i ) = 0. To prove the second identity one needs to list all possible candidates for η (n−2) and then check that there is none. The last identity is more elaborated and makes use of the additional distributional identity (integration by parts) dφ i ι i δ(dφ i ) = −δ(dφ i ) (i is not summed) [29]. In particular, it states that any variation of Y that is as a top form integrated over the entire manifold. The two expressions are clearly equivalent, but their interpretation is rather different. On the left hand side, the connection is computed on a submanifold corresponding to the curve suitably parametrized. On the right hand side instead, the connection is a generically assigned abelian gauge field on M while the geometrical data concerning the path are entirely captured by the PCO. In particular, the latter can be modified as Y (n−1) λ → Y (n−1) λ + dΣ (n−2) , while preserving properties (3.6) and leaving the connection unchanged. This freedom can be exploited to enhance the set of manifest symmetries of Γ; these algebraic properties embody the strength of this method, since it would be much more difficult to ascribe these properties to the curve λ, namely on the homology side.
The Γ integral in (3.7) is manifestly invariant under gauge transformations and deformations of the path within the class of homologically equivalent contours.
Gauge invariance is manifest thanks to the closure property of Poincaré duals (first equation in (3.6)). In fact, under a gauge transformation δA (1) = dα the integral transforms as and the r.h.s. vanishes if ∂M = ∅ or if we impose α to vanish at the intersection λ ∩ ∂M. Invariance of the Wilson loop under a deformation of the path is also easy to study. In fact, from the last identity in (3.6) it turns out that a deformation of the path equations amounts to a shift of Y (n−1) λ by an exact term dη (n−2) . Therefore, integrating by parts, we have JHEP04(2020)161 and the r.h.s. vanishes if the connection has zero curvature on the surface connecting the loop and its deformation, namely if the curve λ has been deformed without encountering singularities. This shows the equivalence between Wilson loops computed on homologically equivalent curves. It is interesting to investigate how to recast in this new framework the identity in (3.2) which states the equivalence between the line integral of the connection A (1) and the surface integral of the field strength F (2) . Given a surface S with ∂S = λ, we call Y (n−2) S the PCO dual to the surface immersed in the space M. Therefore, we can write where we have assumed that d-exact terms integrate to zero. As discussed above, Stokes theorem (or equivalently eq. (3.2)) implies dY is the PCO of the path λ. However, this condition seems to violate the second identity in (3.6).
This apparent contradiction can be sorted out by observing that Y (n−2) S does not have compact support, while Y (n−1) λ is a distribution with compact support. In order to elaborate on this point we assume that locally we can split the manifold as M = M × R + , with the factor R + described by the additional coordinate x . We take λ to be immersed into M only and the surface S to be the union S = λ ∪ {x > 0}. Moreover, we denote by Y where Θ(x ) is the Heaviside theta function, a non-compact support distribution equal to 1 for x > 0, it follows that This is the expected identity which establishes relation (3.2) in the language of integral forms.
Before closing this section, we give a simple formula for the bosonic Wilson loop and the corresponding PCO when the curve is parametrized as τ → x a (τ ), with τ ∈ T ⊆ R.
We enlarge the manifold to M × T with coordinates (x µ , τ ) and we construct the PCO dual to the embedding τ → (x µ (τ ), τ ) as follows It then follows that JHEP04(2020)161 which is the usual expression for a Wilson loop along λ parametrized by τ .
To summarise, we have proposed a new expression for the holonomy of a bosonic Wilson operator as the integral of a top form on the entire manifold (see eq. (3.7)). To our knowledge this is a new formulation, which has never appeared in the literature before. It has the advantage to split the field and the contour dependences, making the investigation of invariances easier. Moreover, it allows for a natural generalization to the supersymmetric case, as we are going to discuss in the next section.

Supersymmetric Wilson loops as integral forms
The supersymmetric version of eq. (3.1) can be defined as [2,43] is the pull-back of the connection superform on a supercurve Λ defined in a supermanifold SM and parametrized by a set of local coordinates z M (τ ) = (x a (τ ), θ α (τ )), a = 1, . . . , n and α = 1, . . . , m. For example, in ten dimensional N = 1 superspace (n = 10, m = 16) the connection superform is given by (C.1), and using definitions (B.3) it can be explicitly written as Generalizing the procedure used in the bosonic case, we construct a super-Poincaré dual which localizes the integrand on the supercurve and allows to rewrite Γ in (3.16) as an integral over the entire supermanifold. Precisely, if the immersion equations of the supercurve Λ in SM are with {φ a } a set of bosonic superfields in SM and {g α } a set of fermionic ones, we introduce a factorized PCO Y The second PCO carries no form degree, but it carries picture number equal to m.
Assigned the PCO, we can rewrite the Wilson loop exponent Γ in (3.16) as The superconnection is generically defined on SM, while the geometrical data featuring the supercurve are captured by the Poincaré dual Y This expression for Γ can be made more explicit if we parametrize the supercurve Λ in terms of smooth functions τ → z M (τ ) on T ⊆ R. For the bosonic part of the PCO we can proceed exactly as done in section 3.1 by including τ as an extra bosonic coordinate and extending the integration to the supermanifold SM × T. A straightforward supersymmetrization of eq. (3.13) leads to For the PCO of the fermionic sector we choose where in the second line we have expanded the Dirac delta functions exploiting the presence of the anticommuting one-form dτ . Here ι β is the contraction along the D β vector field. Using a shorter notation we can then write

JHEP04(2020)161
Focusing on the fermionic part, we can write Γ as where we have used the product m α=1 (θ α − θ α (τ )) to localize the superfield θ-coordinates on the supercurve. Due to the presence of the factor α δ(ψ α ) the only non-vanishing contributions come from terms in the integrand which do not contain any power of ψ α , like for instance A a (x, θ(τ ))dx a from the first term, or terms linear in ψ α where the action of the contraction ι α has the effect to replace ψ α →θ α dτ . Therefore, using the PCO (3.22) to localize also the bosonic coordinates on the supercurve Λ, from eq. (3.25) we easily find In the special case of ten dimensional N = 1 superspace, this expression coincides with (3.16), (3.17) and describes the supersymmetric Wilson operator studied in [9,12]. Similarly, in the four dimensional N = 1 case Γ reduces to the well-known superholonomy and gives rise to the super Wilson loop proposed in [2]. 4 Properties of the fermionic PCO. The fermionic PCO defined in (3.20) satisfies the same properties of the bosonic one, eqs. (3.6). Therefore the total operator Y (n|m) Λ is closed, but not exact and its variations are d-exact.
The last statement is a consequence of a remarkable feature of the fermionic PCO's: given the non-supersymmetric PCO corresponding to immersion functions g α (τ ) = θ α (i.e. θ α (τ ) = 0), then describing an ordinary curve localized at θ α = 0, all the fermionic PCO's turn out to be in the same dcohomological class of Y (0|m) 0 . In order to prove this property we consider a generic Y (m = 1), and using dθδ (dθ) = −δ(dθ), dθδ (dθ) = −2δ (dθ), we can write the following chain of identities so proving the property in the m = 1 case. Since the generalization of the proof to more than one fermionic coordinate is straightforward, we conclude that a generic fermionic PCO is d-equivalent to the non-supersymmetric one, independently of the particular defining function g α (x, θ). It then follows that any pair of PCOs that differ for the choice of the supercontour, i.e. for the choice of the immersion functions, are d-equivalent (clearly, if the two contours are linked by a deformation that does not cross singularities). In particular, this implies that any variation of the PCO induced by a deformation of the path is d-exact, as stated above. The same conclusions remain true when we complete the PCO with its bosonic part Y (n|0) Λ . Although the addition of d-exact terms does not change the cohomological properties of a PCO, it can change its degree of supersymmetry, that is the number of supercharges under which the operator is invariant. We now elaborate on this important point.
Using the geometrical approach, a supersymmetry transformation generated by a spinor acts on the PCO as the Lie derivative (B.6). Exploiting its d-closure property we can write The Y (0|m) 0 operator introduced above breaks supersymmetry completely, δ Y (0|m) 0 = 0. In fact, its defining constraints θ α = 0 are trivially not invariant under supersymmetry transformations, δθ α = α . However, we can perform the shift (we include also the bosonic part) and determine Σ (n−1|m) in such a way that δ Y (n|m) Λ = 0. This condition is equivalent to requiring If this equation is true for arbitrary α , then the shifted PCO is manifestly invariant under all the supersymmetry charges. In general the Killing spinors are a function of τ and supersymmetry is realized locally on the contour. Therefore, by simply adding a d-exact term we can move from a PCO localizing on a non-supersymmetric contour to a PCO localizing on a supersymmetric one. Between these two extreme cases we may have a plethora of intermediate situations where eq. (3.31) holds only for a subset of Killing spinor α components, so defining PCO localizing on partially JHEP04(2020)161 supersymmetric supercontours. An easy way to convince about this fact is to consider for instance the following PCO in ten dimensions Writing V a explicitly as in (B.3), after little algebra one can show that this PCO is d-equivalent to the non-supersymmetric one and is invariant under a supersymmetry transformation generated by the Killing spinor = (1, 0, . . . , 0), that is it preserves only one supercharge. More generally, if in Y (0|16) 0 we replace θ α 1 . . . θ αp with p factors (V a γ a ι) α i we obtain a well-defined fermionic PCO which preserves p supercharges. We note that this procedure can be applied as long as p ≤ 10.
Beyond that limit, we would end up with an exceeding number of V forms that would trivialize the expression. In particular, this construction cannot be used to generate a fully supersymmetric PCO.
In the fully supersymmetric case, we claim that the solution to (3.31) is given by the following expression where we have introduced the Lie derivative along the vector field ∂ τ , the tangent vector along the curve.
In order to support this statement we prove that (3.34) is d-closed and invariant under supersymmetry transformations. To this end, it is convenient to remind the following identities [d, Here we have introduced the super-vector field˜ = (1 − exp(−∂ τ )) α Q α . We note that this super-vector is vanishing in the case of supersymmetry globally defined on the supercontour. From eqs. (3.36) it then follows that it is sufficient to study the closure and the supersymmetry invariance of Y . For sake of clarity, we do the calculation in the simplest case of n = m = 1, being the generalisation lengthy but straightforward. For the d-closure we have

JHEP04(2020)161
whereas for the supersymmetry variation we obtain and the same for δ˜ Y . The results have been obtained by using nilpotence properties like θ 2 = 0 = dx ∧ dx and the usual distributional properties recalled in section 2. Now, inserting back in (3.36) we conclude that Y (n|m) Λ is indeed closed and fully supersymmetric. To close this section it is important to observe that if two PCO's correspond to two different supercontours, and therefore differ by a d-exact term, they give rise in general to two different Wilson operators. In fact, if we start from (3.21) and perform the shift Therefore, by tuning the d-exact term we can flow from one operator to another one. In particular, since different choices of PCO's may correspond to different degrees of supersymmetry preserved by the corresponding supercontours, the d-cohomological equivalence can be used to vary the number of supercharges preserved by the Wilson loop. This will be discussed in detail in section 4.3, whereas in the next subsection we give a first example of this mechanism at work.

The Wilson-Maldacena operator in N = 4 SYM theory
In this section we provide an explicit example of the d-varying supersymmetry mechanism described above by studying the remarkable case of the Wilson-Maldacena loop in four dimensional N = 4 SYM theory [44,45]. We consider the four dimensional N = 4 SYM theory formulated in the (4|16)supermanifold. An ordinary Wilson loop along a curve λ parametrized by τ → x a (τ ), is defined as in eq. (3.25) by taking the non-supersymmetric PCO As already observed, it never preserves any supercharge, no matter is the choice of the contour. Instead, let us consider the d-equivalent PCO

JHEP04(2020)161
Here ι α is the contraction respect to fermionic vector field ∂ α , and N AB is a real vector of the SU(4) R-symmetry group satisfyingN AB = ABCD N CD . Plugging the shifted PCO (3.41) into the general expression for Γ we obtain a shifted holonomy of the form (3.39). If we now replace F (2|0) with its explicit expression (C.18) valid for the N = 4 case, thanks to its non-trivial dependence on the scalar fields, we obtain This expression coincides with the integral of the Wilson-Maldacena generalised connection that includes non-trivial couplings to the six scalars φ [AB] . As is well-known, under a suitable choice of the λ contour and the internal couplings N AB this operator is partially supersymmetric [8]. Therefore, this example proves that d-exact terms can be used to enhance the degree of supersymmetry of a Wilson operator. More generally, if we start from the super-Wilson loop (3.25) corresponding to a generic PCO (3.24) and perform the shift Y , with a similar procedure we find the supersymmetric version of the Wilson-Maldacena operator which has been proposed in [12]. This construction holds for any gauge theory with extended supersymmetry N ≥ 2. In fact, in all these cases the superfield strength F (2|0) contains terms of the form F αIβJ ψ αI ψ βJ , with F αIβJ being proportional to the scalar fields of the gauge multiplet [40][41][42]. Therefore, as in the Wilson-Maldacena example, a careful choice of Σ (n−1|m) leads to an operator which contains non-trivial couplings to the scalar sector.

Variations and symmetries
In this section we study how invariances of a super-Wilson loop can be studied in the language of supermanifolds. As representatives we will consider operators in N = 1 SYM in ten dimensions and N = 4 SYM in four dimensions. We begin by checking invariance under a reparametrization of the path, and then move to the study of invariance under superdiffeomorphisms, supersymmetry and kappa symmetry.

Reparametrisation invariance of the PCO
We start by briefly studying the reparametrisation invariance of the PCO in (3.24). To this end, it is convenient to rewrite it in the following form where we have introduced the volume form

JHEP04(2020)161
Now, under a given reparametrisation τ → σ(τ ), the PCO variation, expressed as usual by a Lie derivative, reads 3) since the object inside the round brackets is odd. This proves the independence of the Γ integral from the contour parametrization.

Variation under superdiffeomorphisms
Given a super-Wilson loop W = e Γ with Γ written as in eq. (3.21), we study its behavior under an infinitesimal superdiffeomorphism generated by a vector field X. This is equivalent to studying how the Γ exponent transforms. Since we have written Γ as a top form integrated on the entire supermanifold and a generic superdiffeomorphism is nothing but a change of coordinates in the supermanifold, we can immediately conclude that by construction Γ, and then W, are manifestly invariant under superdiffeomorphisms. Explicitly, taking into account that for an infinitesimal trasformation the PCO changes by a d-exact If in the second term we integrate by parts and assume that there are no boundary terms, this identity can be equivalently written as for any arbitrary Ω (n|m) form. Identity (4.4) is equivalent to state that in superspace the variation in form of the superconnection induced by the X-tranformation is compensated by the variation of the supercontour Λ encoded in the PCO. In other words, we can write where δ X on the l.h.s. is the X-variation done by keeping the supercontour fixed. 5 When uplifted at the level of the super-Wilson loop, taking into account that a PCO identifies a supercontour uniquely, this implies that (δ X W)(Λ) = −W(δ X Λ). Therefore, the variation of the Wilson operator follows from the X-transformation of the supercontour. In particular, a given X-diffeomorphism is a symmetry for W if (δ X Γ)(Λ) = 0, but from identity (4.6) this is true if and only if δ X Λ = 0. Therefore, the set of W invariances coincides with the set of Λ symmetries. We note that the same reasoning can be applied to bosonic loops defined in ordinary manifolds: (δ X W )(λ) = 0 if and only if δ X λ = 0. JHEP04(2020)161

Supersymmetry invariance
A supersymmetry transformation is a particular superdiffeomorphism generated by the vector X ≡ = α Q α , where Q α are the supersymmetry charges. Therefore, the behavior of a Wilson loop under supersymmetry transformations can be easily infered from the discussion in the previous section. In particular, Γ is manifestly supersymmetric by construction, and from (4.6) we can write (δ Γ)(Λ) = −Γ(δ Λ). This means that its variation is entirely due to the variation of the supercontour. This property has been already discussed in [11,12]. What is interesting to stress here is that in the present formalism, being the Γ's integrand factorized into the product of a contour-independent superfield and a PCO that encloses the whole dependence on the contour, this pattern arises straightforwardly. A Wilson loop preserves a given amount of supersymmetry (it is BPS) when for a particular generator it satisfies (δ W)(Λ) = 0, or equivalently (δ Γ)(Λ) = 0. But, from the previous reasoning this can be traded for the condition Γ(δ Λ) = 0. Therefore, counting the number of supersymmetries preserved by W gets translated into counting the number of supersymmetries preserved by the corresponding supercontour. More precisely, from (4.5) we read up to d-exact terms that we neglect. As discussed in section 3.2, we can exploit the d-equivalence of super-PCO's to vary their degree of supersymmetry. Precisely, given a particular supersymmetry transformation generated by an assigned we can always construct an -preserving PCO from an -breaking operator by performing the shift (3.30), with Σ (n−1|m) satisfying condition (3.31). Therefore, choosing a specific representative within the d-class corresponds to fixing the amount of supersymmetry preserved by the corresponding Wilson loop. Enhancing or de-enhancing supersymmetry can then be done by adding d-exact terms. This result may have important implications in the study of renormalization group flows between Wilson operators preserving different amount of supersymmetry [46][47][48].
Equation (4.7) is the Killing spinor equation selecting the supersymmetry invariances of an assigned Wilson operator. We study it in details, in the ten dimensional case.
First of all, if we express the PCO as in eq. (3.24) and take into account identities (B.8), the ι -contraction on Y (10|16) Λ gives rise to the following two terms Now, according to (4.7), this expression has to be multiplied by F (2|0) . Using the rheonomic parametrization (C.6), it is easy to see that from the first term in (4.8) we obtain a nontrivial contribution both from F ab V a V b and (ψγ a W )V a , whereas from the second term we obtain only one contribution from (ψγ a W )V a , being the V V term trivially zero. Summing all the contributions and factorizing out the volume form (4.2), we finally obtain that the Killing spinor equation reads 2 γ a θΠ b F ab − 2 γ a θ W γ aθ + γ a W Π a Λ = 0 (4.9)

JHEP04(2020)161
where all the quantities are evaluated on the supercontour. When we deal with a supersymmetry preserving PCO, identity (4.5) implies that the following equation has to be automatically satisfied, up to d-terms. There are two possibilities for which this is true. Exploiting the d-closure of the PCO, the first possibility is that ι F (2|0) = dΥ (0|0) on the entire supermanifold, or the even stronger condition ι F (2|0) = 0. These conditions imply a constraint on the gauge field itself and are rarely satisfied. 6 The second possibility is that up to d-terms, which means that ι F (2|0) is vanishing or it is a total derivative on the supercontour only. Using the explicit expression (C.6) for the superfield strength it is easy to check that this condition leads exactly to the Killing spinor equation (4.9). This is a consistency check of the manifest supersymmetry invariance in superspace.
In general, for arbitrary values of the field strengths, equation (4.9) can be solved locally on the contour, leading to a local supersymmetry generated by a Killing spinor (τ ). Remarkably, in the case of a Wilson loop defined on an ordinary bosonic path (θ(τ ) = 0 on the supercontour) it leads to the well-known condition (τ )γ aẋ a (τ ) = 0 (4.12) When reduced to four dimensions, solutions to this equation for constant lead to Zarembolike BPS operators in N = 4 SYM [49]. Instead, in the case of ten dimensional light-like paths, eq. (4.12) has a non-trivial kernel, since it automatically squares to 0. Reduced to four dimensions it defines 1/2-BPS operators in N = 4 SYM if the extra coordinates are identified with the internal couplings to the scalars [8,44]. In this case a systematic classification of solutions to (4.12) has been given in [50], which involves ten dimensional pure spinors. We note that the light-like nature of the contour in ten dimensions is related to kappa-symmetry, as we are going to analyse in the next section.

Kappa symmetry
The superconnection Γ that defines a Wilson loop can be interpreted as the action of a non-dynamical superparticle moving in an electromagnetic field. Since the superparticle in ten dimensions exhibits kappa-symmetry invariance [51], it is sensible to study how the ten dimensional Γ behaves under this symmetry. This has been extensively discussed in [9,12,13]. Here we reformulate the problem in the language of superdifferential forms. In particular, we will confirm the result that kappa-symmetry invariance in ten dimensions is strictly related to BPS properties of the super-Wilson operator in N = 4 SYM theory. 6 We note that this is generically what happens for a supersymmetric invariant action, ι dL (n|0) = 0.
If the action is d-closed (which is possible when auxiliary fields are present), then we have a manifest supersymmetric action. In other cases, the absence of auxiliary fields implies that the action satisfies the weaker condition.

JHEP04(2020)161
A kappa-symmetry transformation is generated by a vector κ ≡ κ α D α , with the kappasymmetry parameter expressed in terms of geometric data as κ α = (γ a ) αβ L a K β (4.13) Here K β is a 0-form carrying a spinorial index and L a is the infinitesimal translation operator. As is well-known, only half of the κ α components are independent. This can be easily understood by proving that the operator (γ a ) αβ L a has a non-trivial kernel, thus allowing to fix half of the fermionic components. An alternative proof, as well as kappa-symmetry transformations of the coordinates, of the basic one-forms and of generic superfields are reviewed in appendix A.

Kappa-symmetry for the super-Wilson loop in 10D
We investigate the action of kappa-symmetry on the Wilson operator W = e Γ , with Γ given in (3.24). Since kappa-symmetry transformations fall into the class of superdiffeomorphisms discussed in section 4.2 the Wilson loop is manifestly invariant under kappa-symmetry by construction. In particular, it has to satisfy identity (4.4) with X = κ, which once again tells us that the Wilson loop variation is entirely due to the variation of its supercontour, i.e. (δ κ W)(Λ) = −W(δ κ Λ). We want to study the WL behavior (δ κ W)(Λ) at fixed Λ and see under which conditions this variation, or equivalently (δ κ Γ)(Λ), vanishes. As just said, and in analogy to what we have done for supersymmetry invariance, this is traded by the following condition Specializing to ten dimensions and using identities (B.14) we have The wedge product with F (2|0) in (C.6) eventually gives with the volume form given in (4.2). Integrating on SM × T we eventually obtain that invariance under kappa-symmetry transformations is ensured by the condition Π a (W γ a κ) Substituting κ with expression (4.13) in momentum representation and localised on the supercontour we end up with δ κ Γ ∝ Π 2 . Therefore, the Wilson loop invariance under kappasymmetry is ensured by the light-like condition, Π 2 (τ ) = 0 at each point of the contour. In the AdS/CFT framework, the worldline kappa-symmetry invariance of the Wilson loop corresponds to the kappa-symmetry invariance of the dual string worldsheet [9]. We note that the fact that we are in ten dimensions has not played any special role in the derivation of this result. Therefore, the same procedure can be applied to super-Wilson loops in 4D without the Wilson-Maldacena terms. Also in that case we find that kappasymmetry is ensured by the light-like condition on the supercovariant momentum [19].

Kappa-symmetry for the Wilson-Maldacena loop in 4D
We now study the kappa-symmetry variation of the four dimensional Wilson-Maldacena connection given in eq. (3.43). As discussed in section 3.3 this connection can be written in terms of an integrable superform associated to PCO (3.41), (3.42), which differs from the PCO localizing the path at θ α = 0 by a d-exact term. Explicitly it is given by where D αA ,D Ȧ α , A = 1, . . . , 4 are the covariant spinorial derivatives in the non-chiral N = 4 superspace.
As discussed above, the kappa-symmetry invariance of the corresponding Wilson loop is ensured when the form is integrated to zero. Being the first term in (4.19) similar to the ten dimensional expression studied in the previous section, its variation can be easily figured out by reducing the previous result (4.17) to four dimensions. We obtain The second term in (4.19) is new and requires a separated analysis. First of all, neglecting d-exact terms, from (4.18) we obtain Now taking the wedge product with the superfield strength given in (C.18), it is easy to realize that only the last two terms there contribute and we are left with We now have to sum the two expressions (4.20) and (4.21), and choose a particular parametrization for the four-dimensional spinors in terms of independent components. The most general expression with the correct index structure is Inserting in the previous equations it is easy to see that mixed Π-N and Π-N contributions cancel, whereas from (4.20) we obtain a term proportional to the four-dimensional Π 2 ≡ Π αα Π αα and from (4.21) an expression proportional to N ABN AB . The total variation δ λ Γ turns out to be proportional to (Π 2 + N ABN AB ). Therefore, invariance under kappa-symmetry requires Π 2 = −N ABN AB . This is the well-known condition that in four dimensional N = 4 SYM theory leads to BPS Wilson loops [8].
Again, this formalism allows for an easy extension to the general case (3.44).

JHEP04(2020)161 5 Generalization to non-abelian gauge groups
The construction of super-Wilson loops in terms of integral forms can be strightforwadly generalized to the case of a non-abelian gauge theory. In fact, it is sufficient to recall that in the non-abelian case the ordinary definition of a gauge invariant Wilson operator reads where λ is a closed path, Tr R is the trace in representation R and the exponential has been generalized to a path ordered exponential. 7 Therefore, in the present set-up it is sufficient to use definition (5.1), but write Γ as in (3.7) for the bosonic operator and (3.21) for the supersymmetric one.
What is interesting to investigate is how in this geometric set-up the invariances of the (super)-Wilson loop discussed in sections 3 and 4 generalize to the case of a non-abelian (super)connection. As prototypical examples, we are going to study gauge invariance of the bosonic Wilson loop and the conditions for supersymmetry invariance of the super-Wilson operator in ten dimensions.

Gauge invariance
For a gauge theory associated to a non-abelian group G, we consider the bosonic W operator expanded as 8 λ is given in (3.13) and localizes the integrands on a closed path λ, while the path-ordered product of PCOs is defined as We note that in (5.2) the path ordering involves only the PCOs, since it is well-defined only for functions living on the contour. Inserting (5.3) in the W expansion and performing the M integrals we are back to the usual path-ordered expansion defined in footnote 7. We consider the gauge variation of (5.2) under (1) , ω] ≡ ∇ω (5.4) 7 We use the convention P λ A (1) * (tn). 8 To simplify the reading we avoid writing explicitly the wedge product symbol. Moreover, we introduce the notation Y

JHEP04(2020)161
where ω is a smooth function on the M manifold with values in the Lie algebra of G. Due to the second term in this transformation, the gauge invariance of (5.2) requires cancellation of terms arising from different orders in the expansion. We are going to check gauge invariance up to cubic order in the connection.
We start discussing the variation of the linear term in (5.2). At this order gauge invariance easily follows from the chain of identities where in the first line we have used ∇Y We now move to the second order term in (5.2). We begin by considering the contribution coming from δA (1) → dω. It is explicitly given by In the last step we have integrated by parts the differentials d 1 and d 2 acting on x 1 and x 2 coordinates, respectively. Now, we can use the following identities (we refer to appendix E for their proof) where the minus sign in the second equation is due to the path ordering. We can then write (5.6) as

JHEP04(2020)161
where in the last line we have integrated in the (x 2 , τ 2 ) variables using the identity δ(τ 1 − τ 2 ) = −δ(τ 2 − τ 1 ) in order to preserve the orientation of the loop. By taking the trace this term eventually vanishes. We now consider the contribution from the variation of the second order term in (5.2) under δA (1) This term is not vanishing itself, but it is expected to compensate the variation of the cubic term in (5.2) under δA (1) (x) → dω(x). In fact, integrating by parts, the variation of the cubic term gives rise to As described in appendix E (see for instance eq. (E.6)), the action of the d j=1,2,3 differential on the path-ordered product of PCO's has the net effect to replace 9 Θ(τ j−1 − τ j ) with −dτ j δ(τ j−1 −τ j ) and correspondingly Y (n) In both cases the delta functions allow to perform the dτ j n a=1 dx a j integrations, so reducing (5.10) to a double integral. Moreover, having the two terms opposite sign, we can easily reconstruct a commutator [ω(x j ), A (1) (x j )] for every j = 1, 2, 3. Exploiting the symmetries of the integrand under the exchange of integration variables the six terms in each path-ordered product give eventually the same contribution, so that we end up with This expression cancels exactly the contribution in (5.9). We have then proved gauge invariance of (5.2), up to cubic order. However, it is an easy task to realize that the same pattern keeps repeating order by order, so ensuring gauge invariance of the complete Wilson operator.
A similar analysis holds also in the case of non-abelian super-Wilson loops. JHEP04(2020)161

Supersymmetry invariance
We now consider the non-abelian generalization of a super-Wilson loop, W = Tr R P e Γ with super-holonomy Γ given in (3.21), and study its variation under supersymmetry transformations. For the abelian case, in section 4.3 we have discussed conditions that can be imposed in order to have either local or global supersymmetry. In this section we briefly show that in the non-abelian case slight differences arise. Given an expansion for W in powers of the supergauge connection similar to (5.2), we will restrict to the study of the linear term. We first recall that the supersymmetry variation of a non-abelian gauge superfield, obtained as usual by the action of a Lie derivative, can be expressed as , ·] is the gauge covariant differential in superspace. The supersymmetry variation of the linear term in the W expansion can then be written as For closed paths the first term vanishes identically. Therefore, the only term that may affect the supersymmetry invariance of W is the second one. Following the discussion in section 4.3 and adapting it to the non-abelian case, this term drops out if we require ι F (2|0) Λ = ∇Υ (5.14) As in the abelian case, neglecting ∇-exact terms this condition reduces to (4.9).

Relating Wilson loops and pure spinor vertex operators
The Wilson loop expectation value W , which describes the motion of a superparticle along a path Λ in a gauge background has a stringy interpretation within the AdS/CFT correspondence [44,45]. In fact, being the particle excited by an open massless vertex operator at the boundary of the string worldsheet, the Wilson loop expectation value equals the string partition function on a worldsheet ending on Λ at the boundary. 10 In particular, in the α → 0 limit the partition function can be computed semiclassically and leads to a prediction for W at strong coupling [8,44,45].
On the other hand, in the pure spinor approach to string theory the integrated vertex operator for the massless spectrum of the open superstring reads [52] where d α is a worldsheet field related to the conjugate momentum to θ α and N ab is the Lorentz generator in the pure spinor space. At quantum level it is invariant under the BRST transformations [52] Qd α = Π a (γ a λ) α , QN ab = d α (γ ab ) αβ λ β (6.2)

JHEP04(2020)161
where the nilpotency conditions Q 2 d α = Q 2 N ab = 0 follow from the requirement for λ α to be a commuting pure spinor, i.e. λγ a λ = 0. Comparing equation (6.1) with the expression for the ten dimensional superholonomy Γ given in (3.26) we see that the first two terms are identical. Therefore, from the perspective of relating Wilson loops to open string worldsheets, we investigate whether it is possible to modify Γ in such a way to obtain an expression formally identical to the string vertex operator. Indeed, we show that this is possible by applying a d-deformation to the PCO (3.24) along the lines described in section 3.2.
To prove this statement we deform the original PCO in (3.24) as where we have introduced the two-vector N ab and the ten dimensional spinor d α . In order to compute dΣ (9|16) we need to specify how the differential acts on these new fields. In analogy with the action of the Q operator in eq. (6.2) we propose We note that, without imposing any pure spinor constraint, these definitions automatically satisfy the Bianchi identities d 2 d α = 0 and d 2 N ab = 0. In particular, for N ab this is guaranteed by the addition of the extra term − 1 2 V a ∧ V b , which is instead absent in (6.2). Given these definitions we can now evaluate how the original superholonomy gets modified. Recalling that we focus only on the new term proportional to F (2|0) . Inserting (6.3) and the superfield strength (C.6), it is explicitly given by (we shortly indicate d α (γ a 10 ) αβ ∂ ∂ψ β ≡ dγ a 10 ι) ..a 10 V a 1 ∧ V a 8 N a 9 a 10 δ 16 (ψ) + dτ a 1 ...a 10 V a 1 ∧ V a 9 dγ a 10 ι δ 16 (ψ) In the first term we simply antisymmetrize the vielbeins to obtain a desidered term proportional to F ab N ab times a factorized volume form V 10 δ 16 (ψ). In the second term we first integrate by parts ι on the ψ spinor and, after a bit of algebra, we produce a contribution proportional to W α d α times a factorized volume. In total, summing the two terms we obtain SM×T (F ab N ab + W α d α ) dτ V 10 δ (16) (ψ) (6.7)

JHEP04(2020)161
We can now project the integrand onto the Wilson path by performing the integrations on the supermanifold coordinates. The obtained contributions, when added to the original Γ as in (6.5), reproduce the pure spinor vertex operator (6.1) written as a supermanifold integral.

Conclusions
We have constructed super-Wilson operators in terms of integral forms describing the immersion of the supercontour in a supermanifold. In such a formulation the corresponding superholonomy is written as an integral over the entire supermanifold and the invariance of the operator under superdiffeomorphisms becomes manifest. As a by-product, we obtain an alternative description also of the ordinary Wilson loops, which can be obtained from the supersymmetric one by setting all the spinorial coordinates to zero. We have reformulated kappa-symmetry in this language and studied the Killing spinor equations associated to supersymmetry invariance.
We have highlighted the role of the d-cohomology in the construction of the Picture Changing Operators (PCO). Different PCOs corresponding to different supercontours are all comohological equivalent. Nevertheless, they may preserve a different amount of supersymmetry and more generally they exhibit a different spectrum of symmetries. In particular, it follows that by adding d-exact terms we can tune the BPS degree of a Wilson operator and we can easily relate two operators which preserve a different fraction of supersymmetries. As a remarkable example, we have shown how the BPS Wilson-Maldacena loop of N = 4 SYM theory can be obtained from the non-BPS one by the addition of a suitable d-exact term to the ordinary non-supersymmetric PCO.
It would be interesting to generalize our formulation at quantum level to compute perturbative corrections to Wilson loops. In particular, it would be nice to understand which is the effect of the d-varying symmetries mechanism, in a frame where d-exact terms could be treated as perturbations. More ambitiously, it could be interesting to understand how to reformulate localization in a geometrical framework and exploit our expression for the Wilson loop to compute its vacuum expectation value exactly.
Finally, as emphasised in the paper, this formalism allows for a straightforward generalization to curved supermanifolds, hence leading to Wilson operators defined in a supergravity framework [56,57]. This geometrical setting might be also applied to Wilson loops in different dimensions, for example to the well known bosonic BPS Wilson loops in three dimensional Chern-Simons-matter theories [58][59][60][61].

A Superdifferential forms
We briefly review definitions and properties regarding superdifferential calculus that have been used in the main body of the paper. For a more extensive introduction to this topic we refer the reader for example to [22,24,29,32,57,62].
We consider a supermanifold SM with n bosonic and m fermionic dimensions. We denote the local coordinates in an open set as (x a , θ α ). A (p|q)-form ω (p|q) has the following structure where the dθ α appearing in the product are independent of those appearing in the deltas (α i = β j for any pair i, j) and From the first identity of the second line we note that δ(dθ) has to be treated as an anticommuting object, unlike the standard δ distribution. This is due to the fact that δ(dθ) is used to compute the oriented volume of the supermanifold. Indeed, δ (dθ) is not a distribution on smooth functions, but rather on "smooth differential forms". This is mathematically called a de Rham current (see [24] for further explanations). The two quantum numbers p and q in eq. (A.1) correspond to the form number and the picture number, respectively, and they range as −∞ < p < +∞ and 0 ≤ q ≤ m. The total form degree is given by p = r + s − i=q i=1 b i since the derivatives act effectively as negative forms and the delta functions do not carry any form degree. The total picture q of ω (p|q) corresponds to the number of delta functions. In particular, we call it superform if q = 0, ω (p|0) (x, θ, dx, dθ) = ω(x, θ)dx a 1 . . . dx ar dθ α 1 . . . dθ αs , p = r + s , Otherwise it is called pseudoform. A top integral form ω (n|m) corresponds to an element of the line bundle known as Berezinian bundle (the transition functions are represented by the superdeterminant of the JHEP04(2020)161 Jacobian) and it can be locally expressed as in eq. (A.4) with p = n. As in conventional geometry, we can define the integral of a top form on the superspace T * SM endowed with a super-measure [dxdθd(dx)d(dθ)] as where the order of the integration variables is kept fixed and the measure is invariant under coordinate transformations. We refer the reader to [24] of the ω(x, θ) superfield. In the present formulation the Stokes theorem for integral forms is also valid. By changing the one-forms dx a , dθ α as a top form ω (n|m) transforms as where Ber(E) is the superdeterminant (i.e. Berezinian) of the supervielbein (E a , E α ).

B Superspace conventions
In this appendix we collect the conventions on supermanifolds that we have used in the main text. To be definite we focus on the N = 1 superspace in ten dimensions described by coordinates z M = (x µ , θ α ), with µ = 0, . . . , 9 and α = 1, . . . , 16. Here θ α are Majorana-Weyl spinors. Introducing ten dimensional 16 × 16 gamma matrices γ a αβ , 12 supercharges and supercovariant derivatives are defined as

JHEP04(2020)161
with Q α = D α + 2θ β γ a αβ ∂ a . They satisfy Flat supervielbeins are defined as and satisfy the Maurer-Cartan equations dV a = ψγ a ψ , dψ α = 0 (B.4) Introducing the ordinary super-differential basis (dx a , dθ α ) through the defining identities ∂ a , dx b = δ b a and ∂ α , dθ β = δ β α , the differential operator d is given by Supersymmetry transformations. In the present framework, a supersymmetry transformation is a superdiffeomorphism whose action on differential forms is represented by a Lie derivative along the vector field = α D α + 2 γ a θ∂ a ≡ α Q α , where in general ι v is the contraction operator defined on a p-form as Using ι θ α = ι x a = 0, it then follows that which are the ordinary supersymmetry transformations of the superspace coordinates. In particular, these defining identities follow ι V a = 2 γ a θ , ι ψ α = α (B.8) As a consistency check of our conventions we find that (V a , ψ α ) are invariant under supersymmetry as a consequence of identities (B.4) and (B.7) L ψ α = ι dψ α + dι ψ α = 0 , L V a = ι dV a + dι V a = 2 γ a ψ + d(2 γ a θ) = 0 (B.9) The supersymmetry variation of a scalar function Φ on the supermanifold is defined as where Q α is the supersymmetry generator introduced in eq. (B.1). Similarly, the supersymmetry variation of a superconnection A (1|0) reads where the first term is a gauge transformation of A (1|0) and in the second term F (2|0) is the superfield strength defined as

JHEP04(2020)161
Kappa-symmetry transformations. A kappa-symmetry transformation is generated by a vector κ ≡ κ α D α which differs from the supersymmetry generator by the simple replacement In particular, it follows that a kappa-symmetry transformation can be formally written as the action of the Lie derivative L κ = {ι κ , d} where When ι κ acts on the flat supervielbeins (V a , ψ α ), using ι a V b = δ b a we find Therefore, kappa-symmetry transformations of the superspace coordinates and the basic one-forms read In the case of rigid symmetry, dκ α = 0. It is interesting to note that the replacement of the Q with the D generators in (B.12) leads to the following dual situation between supersymmetry and kappa-symmetry transformations of the bosonic supervielbein Applied to a generic superfield Φ, a kappa-symmetry transformation reads and can be correctly obtained from a supersymmetry transformation (B.10) by replacing the Q α generator with D α . Parametrizing theκ generator as with K β a 0-form and L a the infinitesimal translation operator, we shift K β as Consequently, the κ α parameter transforms as where we have exploited [L a , L b ] = 0. If K is an harmonic function, then transformation (B.22) is a symmetry of the kappa-symmetry parameter, i.e. K β and K β give rise to the same kappa-symmetry transformation. This degeneracy can be used to halve the number of independent components of the K spinor.

JHEP04(2020)161 C Geometry of supersymmetric gauge fields
We recall some basic facts about the geometrical construction of supersymmetric gauge theories. We discuss the abelian and non-abelian theories separately. We primarily consider the 10D N=1 SYM theory, since other cases can be obtained by dimensional reduction and suitable truncations. The special case of 4D N = 4 abelian theory is reviewed at the end of this section.
The abelian case. The 10D gauge supermultiplet is represented by one vector superfield and one spinor superfield (the gaugino) with degrees of freedom matching on-shell [63][64][65].
No off-shell superspace formulation is known, which includes the correct spectrum of auxiliary fields allowing to construct a superspace action that leads to the correct equations of motion. However, a super-geometric formulation can be developped, which stems from promoting the gauge field to a superfield (1|0)-superform (with V a and ψ α defined in (B.3)) The corresponding field strength, defined as is subject to Bianchi identities supplemented by the conventional gauge invariant constraint from which we obtain A a as a function of the spinorial components A α . As a consequence, the other components turn out to be uniquely expressed in terms of the gaugino (0|0)superform W α F aα = (γ a W ) α , F ab = (Dγ ab W ) (C. 4) and satisfy the additional constraints These constraints automatically imply the equations of motion for all the physical fields. By using suitable gamma matrices identities, one can prove that the previous relations can be recast in the following superform equations The great advantage of using the geometric formulation is that supersymmetry transformations can be expressed as superdiffeomorphisms along the fermionic directions (see appendix B for the geometric definition of supersymmetry transformations). In particular, the gauge superfields transform as

JHEP04(2020)161
These relations give rise to the ordinary supersymmetry transformations up to a gauge transformation of the gauge field A, while the gaugino superfield W α is gauge covariant. We note that these rules remain true also in the case of local transformations.
The non-abelian case. As for the abelian case, an off-shell superspace formulation of gauge superfields with auxiliary fields is not known, but a geometric formulation can be provided [63][64][65]. For a non-abelian gauge group the superfield strength is defined as where the A (1|0) superconnection is still expanded as in (C.1). The superfield strength is subject to the Bianchi identities ∇F (2|0) = 0 (C.9) with the covariant derivative defined as ∇F (2|0) = dF (2|0) + [A (1|0) , F (2|0) ] = 0. This is supplemented by the conventional gauge invariant constraint from which one obtains A a as a function of the spinorial components A α . The other components turn out to be expressed in terms of the gaugino superfield W α as and satisfy the additional constraints Equations (C.11), (C.12) imply the equations of motion, which are then a consequence of the superspace constraints. Supersymmetry transformations are easily expressed as These relations give rise to the ordinary supersymmetry transformations up to a gauge transformation of the gauge field A, while the gaugino superfield W α is gauge covariant.
Dimensional Reduction to 4D. As is well-known, D = 4, N = 4 SYM theory can be obtained by dimensional reduction of the D = 10, N = 1 theory, while preserving the maximal amount of supersymmetry. Here we clarify how to perform the dimensional reduction in the geometric set-up. For simplicity, we restrict to abelian theories.
Starting from the ten dimensional superform (C.1), we first perform the following decompositions Here V αα can be identified with the components of the four-dimensional vielbein, whereas V [AB] is the vielbein along the extra six directions. It satisfies the self-duality constraint . Similarly, ψ = (ψ A,α ,ψα A ) represents the decomposition of the rigid gravitino field. They satisfy the following equations In the same way, we decompose the gaugino superform W α = (W A α ,W Aα ) according to its SL(2, C) × SU(4) representation.
The dimensional reduction is then achieved by removing the dependence of the fields upon the transverse coordinates y [AB] . The four A αα components then describe the gauge connection in four dimensions, φ [AB] are the six real scalars of the N = 4 SYM theory and W A,α give rise to the four gaugini. As a consequence, from the definition of the field strength F (2|0) in (C.6) we obtain As described in [66], in order to complete the dimensional reduction we have to redefine the connection as where Φ [AB] are six chiral superfields containing the φ [AB] scalars. As a consequence, the superfield strength becomes F (2|0) = V αα ∧V ββ F αα,ββ −V αα (ψ A,α W A α +ψ A αWA,α )+( ABCD αβψα Cψα D + αβ ψ A,α ψ B,β )Φ AB (C. 18) and coincides with the expression for the superfield strength of the N = 4 SYM theory obtained directly in four dimensional non-chiral superspace (see for instance [40][41][42]). We note that additional pieces proportional to the flat gravitinos appear, which carry an explicit dependence on the scalar fields φ [AB] . As we explain in the main text, these terms are crucial for the construction of the supersymmetric version of the Wilson-Maldacena loop in four dimensions. The dimensional reduction of non-abelian theories works exactly in the same way, with the obvious covariantization of the equations. where Vol (n) is the volume form on M, given by Vol (n) = A (1) ∧ n−1 a=1 dφ a . This formula can be proved as follows. Given the Poincaré dual Y (n−1) λ in eq. (3.3) that defines the immersion of the λ curve, we choose a (local) basis of vectors {X a } n−1 a=1 of λ ⊥ normalised by ι a dφ b ≡ X a (φ b ) = δ a b . Then, it is easy to prove that where in the last equality we have used the Leibnitz rule and ι Xa 1 dφ a 1 ∧ . . . ∧ dφ a n−1 ι Xa 2 . . . ι Xa n−1 Vol (n) = 0 (D.4) since the form inside the brackets is an (n + 1)-form in an n-dimensional manifold.

E Proof of identity (5.7)
Here we want to compute the following expression where d j ≡ dx a j ∂ ∂x a j + dτ j ∂ ∂τ j , j = 1, . . . , M , denotes the exterior derivative w.r.t. the set of coordinates (x a j , τ j ) acting on the path-ordered wedge product of M PCO's of the form (3.13). This expression is required to prove gauge invariance of the non-abelian Wilson operator, as discussed in section 5.1.
All the PCO's localize on the same contour, but parametrized by different parameters τ i ∈ T. The definition of the corresponding path-ordering is 13 P Y  13 In order to avoid cluttering, in what follows we will neglect the wedge symbol.
We now generalize the calculation to the product of M PCO's. For the sake of clarity, we focus on a single term of (E.2), namely and first consider applying d 1 . Since d 1 acts on a single theta function, we obtain the following chain of identities If we now apply d 2 we have to take into account that this time the differential acts on two different theta functions. Therefore, in this case we obtain We see that the first term is exactly minus the term in (E.4), whereas the second term will coincide with minus one of the two terms which arise when we apply d 3 . This pattern JHEP04(2020)161 repeats itself for any other differential acting on intermediate theta functions. The d j differential will produce two terms the first one being opposite in sign to a term coming from d j−1 and the second one opposite to a term from the application of d j+1 .
The same pattern holds for any other term of (E.2) when we consider the contributions coming from the application of d σ(τ i−1 ) , d σ(τ i ) and d σ(τ i+1 ) on the product Θ (σ(τ i−1 ) − σ(τ i )) Θ (σ(τ i ) − σ(τ i+1 ))). Precisely, the d σ(τ i ) derivative produces two terms which come in pair with opposite signs with one term from d σ(τ i−1 ) and one from d σ(τ i+1 ) . As we discuss in section 5.1, these signs are crucial for reconstructing commutators and ensure cancellation in the gauge variation of the Wilson loop.
It is important to observe that an identical proof works also in the case of a supermanifold SM, i.e. for products of super-Poincaré duals of the form Y (n|0) ∧ Y (0|m) (see eq. (3.20)) localizing super-integrals on supercontours.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.