From harmonic maps to the nonlinear supersymmetric sigma model of quantum field theory. At the interface of theoretical physics, Riemannian geometry and nonlinear analysis

Harmonic maps from Riemann surfaces arise from a conformally invariant variational problem. Therefore, on one hand, they are intimately connected with moduli spaces of Riemann surfaces, and on the other hand, because the conformal group is noncompact, constitute a prototype for the formation of singularities, the so-called bubbles, in geometric analysis. In theoretical physics, they arise from the nonlinear $\sigma$-model of quantum field theory. That model possesses a supersymmetric extension, coupling a harmonic map like field with a nonlinear spinor field. In the physical model, that spinor field is anticommuting. In this contribution, we analyze both a mathematical version with a commuting spinor field and the original supersymmetric version. Moreover, this model gives rise to a further field, a gravitino, that can be seen as the supersymmetric partner of a Riemann surface metric. Altogether, this leads to a beautiful combination of concepts from quantum field theory, structures from Riemannian geometry and Riemann surface theory, and methods of nonlinear geometric analysis.


Introduction
At his premature death, Eberhard Zeidler left his magnum opum unfinished, a comprehensive mathematical treatment of quantum field theory. He had planned six volumes, of which three, with more than 3000 pages, [81,82,83], could be realized, and it is planned that one of us (J.T.) can complete the fourth volume on the basis of the notes that he had left behind. His vision consisted into a profound interpenetration of the two fields, leading to deeper insight and new methods in both of them.
Indeed, the action functionals of (quantum) field theory usually contain fundamental geometric structures that are incorporated in their symmetries and invariances. In turn, analyzing their Euler-Lagrange equations-or physically speaking, field equations-from a rigorous mathematical perspective leads to challenging problems in nonlinear partial differential equations and the calculus of variations. And solving those problems can provide us with insights and methods that can find profound novel applications in geometry, as for instance the Yang-Mills and Seiberg-Witten equations have shown. And superstring theory, while controversial among physicists, has triggered important advances in various areas of mathematics.
In this spirit, in this contribution, we discuss some of the basic models of quantum field theory, the nonlinear σ-model. In its most basic version, mathematically, it involves harmonic mappings from a Riemann surface to a sphere or some other model space. Since the target is not Euclidean, but rather a symmetric space, that is, a particular Riemannian manifold, the Euler-Lagrange equations of the model are generically nonlinear. Their mathematical treatment uncovered the phenomenon of bubbling of solutions of conformally invariant variational problems, ultimately leading to the development of the mathematical theory of pseudoholomorphic curves and Gromov-Witten invariants. Also, there are many analogies to the four-dimensional Yang-Mills equations.
From the perspective of physics, the ultimate version of the model is still richer, as it includes supersymmetry. That is, the original bosonic fields are supplemented by fermionic partners. This brings us into the realm of supergeometry and, in particular, the theory of super Riemann surfaces.
Here, on one hand, we confront that challenge and develop the mathematical theory of the supersymmetric nonlinear σ-model. As we shall see, the geometric perspective, that is, the moduli space of super Riemann surfaces, and the physical or variational one, the symmetries and invariances of the action functional, are two sides of the same coin, and either of them deeply reflects the other. More precisely, the Noether currents associated to the invariances of the functional yield the cotangent vectors of the moduli space of super Riemann surfaces.
We have, however, discovered more. There is also a version of the theory that does not need supergeometry. This converts the formal expansion that supergeometry works with into nonlinear elliptic systems of partial differential equations that extend the harmonic map system. These systems on one hand provide us with a test bed for the development of new analytical methods that should eventually find applications in the wider realm of geometric analysis, and on the other hand uncover some novel geometric structure that seems to be worth in exploring further, in order to create new types of invariants of Riemannian manifolds.
We gratefully recall the inspiration that Eberhard Zeidler has given us over the years, through his works, his generous and warm personality, and his scientific vision, and we hope that it is not inappropriate to dedicate the present work to his memory.

The geometric and analytic context: Harmonic mappings between Riemannian manifolds
Harmonic mappings from Riemann surfaces into Riemannian manifolds can be considered as special cases of harmonic mappings between Riemannian manifolds. While, as we shall see and explore, when the domain is a Riemann surface, harmonic mappings enjoy special properties, on which the rest of this article will essentially build, it is nevertheless useful to first understand the wider context. Thus, let M be a Riemannian manifold of dimension m. In the sequel, we shall use basic concepts of Riemannian geometry without further explanation; for example, T M is the tangent bundle of M , and T x M is the tangent space over the point x ∈ M . The reader may consult [34]. We shall usually implicitly assume that all manifolds are compact, without boundary. While boundary value problems are also analytically challenging and geometrically interesting, for this article, we have decided not to go into that topic.
More precisely, the Riemannian metric is the symmetric 2-tensor but for most purposes, we shall just work with its coefficient matrix (1). From (2), the reader can already see that we use the standard Einstein summation convention. We shall also use the following notations The volume form of M , |γ| dx 1 · · · dx m in local coordinates, will often be abbreviated as dvol γ . The Laplace-Beltrami operator of M is We employ here a sign convention different from that in [34], where ∆ M has a minus sign, in order to make it a positive operator.

Harmonic functions
For a function φ : M → R, we have the Dirichlet integral Here, we have used an invariant and a coordinate notation, that is, The critical points of E are the harmonic functions of M , that is, the solutions of the corresponding Euler-Lagrange equation For the relevant results in the calculus of variations, we refer to [38], and for the PDE aspects, to [33]. We need to recall a few facts here, however. The Dirichlet integral (6) is naturally defined on the Sobolev space W 1,2 (M ) of functions that possess square integrable weak derivatives, but need not be any smoother than that. In fact, they need not even be continuous, let alone differentiable. Therefore, for such a Sobolev function, a priori (8) need not make sense. Thus, from minimizing E, one first obtains a so-called weak solution, that is, a φ solving for all test functions η ∈ W 1,2 (M ). One then shows, and this is the prototype of elliptic regularity theory, that any solution of (9) is in fact smooth and therefore solves (8), as one finds by integrating (9) by parts to move the first derivative from η, which is allowed when φ is smooth. We refer to [33] for the systematic theory. Here, smoothness of a weak solution is easy because (8) is a linear elliptic partial differential equation. Of course, this assumes that the metric γ itself is smooth and bounded, and also that |γ| is bounded from below, so that the inverse metric tensor γ αβ is likewise smooth and bounded. More precisely, one needs the uniform ellipticity condition λ ξ 2 ≤ |γ|γ αβ ξ α ξ β ≤ Λ ξ 2 for all ξ = (ξ 1 , . . . ξ m ) and some 0 < λ, Λ < ∞.
When this holds, and the quantities |γ|γ αβ are measurable, then by the De Giorgi-Nash theorem, a weak solution of (9) is Hölder continuous. And if these coefficients possess more regularity, then so does the solution φ. In fact, by Schauder's theorem, if they are of class C k,α for some k ∈ N, 0 < α < 1, then φ ∈ C k+2,α . Proofs of these results can be found in [33]. Anticipating the key discussion of later sections, however, we not only want to treat the function φ as a variable, but also the metric γ. Even though later, for two-dimensional domains, which are those of main interest in this article, we shall also derive equations for variations of the Dirichlet integral w.r.t. the metric γ, these will not be differential equations that will yield any regularity for the metric. Rather, there is some gauge invariance principle. That principle says that, in dimension two, again, we can always choose local coordinates so that the metric becomes the Euclidean metric, that is, constant, and therefore as smooth as one wants. And then φ is also smooth to any order, in fact even analytic. However, when we still move further, to the supersymmetric sigma model, we get further variables, a spinor field that is the partner of the function φ, and a gravitino field that is the partner of the metric γ. And then, subtle regularity issues will arise. But before that, we need to generalize things in another direction.

Harmonic mappings
We shall now move to nonlinear elliptic systems where the regularity for φ is no longer so easy, and in fact, not even true in general. We move into this nonlinear realm by replacing the target R by another Riemannian manifold (N, g) of dimension n, with metric tensor g = g ij dx i ⊗ dx j . Analogously to (3), (4), we define g ij to be the coefficients of the inverse metric and the Christoffel symbols Γ i jk . For a map φ : M → N in the appropriate Sobolev space W 1,2 (M, N ) (see [34] for the precise definition), we define its energy 1 of φ : M → N is where in the last expression, we use local coordinates (x 1 , . . . , x m ) on M and (φ 1 , . . . , φ n ) on N . Analogously to (7), The energy integral thus is defined in terms of the Riemannian metrics of M and N . When N = R, this reduces to (6). By a somewhat tedious, but otherwise straightforward, computation, we can generalize (8).
We usually write this system of equations as and call τ (φ) the tension field of φ.
The Euler-Lagrange equation (11) constitutes a semilinear elliptic system. It is in general nonlinear, because of the geometry of N that enters through the Christoffel symbols Γ i jk (φ(x)) that are contracted with first derivatives of φ, resulting in an expression that is quadratic in those derivatives. Therefore, weak solutions, defined analogously to (9), in general need not be regular. Since the regularity theory becomes much better when the domain is two-dimensional, our exclusive case of interest below, we don't enter here into that general theory, however.
As we have seen, harmonic maps are generalization of harmonic functions. Both the Laplace-Beltrami operator ∆ M and its generalization, the tension field τ , involve the Levi-Civita connection of M , the latter also that of N .
The concept of harmonic maps was first proposed by S. Bochner [3] in 1940. See also [22]. Subsequently, it was developed very extensively, and in fact too extensively to be reviewed here. We refer the reader to [30] for a survey.
In particular, one may apply variational methods to the action functional E. As an alternative, one may use a parabolic method, that is, study the flow where φ now is a map from M × [0, ∞) to N . Of course, one needs to show that a solution exists for all t ∈ [0, ∞) and converges to a solution of (11) as t → ∞. This is not always true, but holds when the target space N has nonpositive curvature. Like the regularity theory for (12), the subject has been too extensively studied to be properly reviewed here.

Harmonic maps from Riemann surfaces and conformal invariance
In this section and the following, we collect some material that is discussed in more detail in [28,29,32,30,31,34,53], where also further references can be found. The energy integral, and therefore harmonic maps, have some special property when the dimension of the domain is two, because of conformal invariance, as we are now going to explain. We denote the domain by Σ and let it be an oriented two-dimensional Riemannian manifold with metric γ αβ . In the sequel, we shall use different versions of the uniformization theorem, and we want to describe them now. First, there is the local uniformization theorem, saying that locally, any metric on a surface can be transformed into the Euclidean one by a suitable coordinate transformation. This was first proved by Gauss for smooth metrics. It was an important step towards the modern regularity theory of solutions of elliptic partial differential equations when Lichtenstein [61] proved such a result for metrics of class C 1,α . Lavrent'ev [59] then showed such a result for continuous metrics, and Morrey [62] obtained it for bounded ones; of course, all metrics have to be uniformly positive definite. See also the presentation in [26]. A global uniformization theorem was first stated by Riemann, and it is known as the Riemann mapping theorem. The final version was proved by Koebe who showed that any non-empty simply connected open subset U of the complex plane C other than C itself is conformally equivalent to the open unit disk. That means that there exists a bijective conformal map from U onto the unit disk. Finally, there is Poincaré's uniformization theorem which says that the universal cover of a compact Riemann surface is conformally equivalent to either the unit sphere (for genus 0), the complex plane (for genus 1), or the hyperbolic plane (for genus p ≥ 2). By the uniformization theorem of Gauss, Σ then obtains the structure of a Riemann surface, that is, we can locally find holomorphic coordinates z = x + i y on Σ such that the metric tensor (1) becomes where λ(z) is a real-valued positive function. Here, dz = dx + i dy, dz = dx − i dy. Conversely, given a compact Riemann surface Σ, by the Poincaré uniformization theorem, we can find a conformal metric as in (13) with constant curvature; we shall be mainly interested in the case of genus ≥ 2, the hyperbolic case, where this curvature is negative, −1. This metric is unique (up to diffeomorphism, see below). As such a surface Σ is of the form H/Γ, where H := {x+i y ∈ C : y > 0} is the upper half plane and Γ is a discrete subgroup of the isometry group PSL(2, R) of H. Here the group PSL(2, R) acts on H by Möbius transformations. The metric on S descends from the hyperbolic metric 1 y 2 dz ⊗ dz on H. As before, let N be a Riemannian manifold with metric tensor (g ij ). The energy (10) of a map φ : Σ → N then becomes The last equation implies that the energy of a map from a Riemann surface Σ into a Riemannian manifold is conformally invariant in the sense that it does not depend on the choice of a metric on Σ, but only on the Riemann surface structure. This only holds when the dimension of the domain is 2, that is, when we look at maps from a surface. The reason is that the volume form behaves like an n 2 -th power of the metric, and therefore, only for n = 2, this can cancel the effect of the inverse metric tensor in the energy integral. Nevertheless, as we shall explore below, it turns out to be useful to keep the domain metric λ 2 in the picture, even though, by the observation we have just made, this introduces some redundancy.
Also, if k : Σ → S is a bijective holomorphic or antiholomorphic map between Riemann surfaces then for any φ : S → N (of class C 1 ) and if φ is harmonic, then so is φ • k.
The harmonic map equation (2.1) is then also independent of the choice of conformal metric on the domain. This means that the map φ : Σ → N of class C 2 is harmonic iff When the target is also a Riemann surface S with metric then the harmonic map equation (15) becomes Thus, holomorphic or antiholomorphic maps between Riemann surfaces are harmonic as they obviously satisfy (16). In the sequel, we shall write ± holomorphic to mean holomorphic or antiholomorphic. The converse does not necessarily hold. Harmonic maps exist in any homotopy class of maps between compact Riemann surfaces (unless the target is a sphere, that is, of genus 0), but ± holomorphic maps h : Σ → S have to satisfy the Riemann-Hurwitz formula where p, q are the genera of Σ, S, m is the degree of h, and v h ≥ 0 is its total ramification index. In particular, as a necessary topological condition for a harmonic map to be ± holomorphic, we have the inequality where χ denoteos the Euler characteristic. In fact, even if (18) holds, in general there will exist no holomorphic map between Σ and S, as this requires that the conformal structure of Σ be a conformal cover of the one of S. Thus, in general, harmonic maps between Riemann surfaces Σ and S cannot be ± holomorphic. Also, as we see from (15), whether a map is harmonic depends on the metric of the target S whereas the property of being ± holomorphic only depends on the conformal structure of S, but not on its metric. Nevertheless, harmonic maps into Riemann surfaces enjoy special properties, as demonstrated by the following result of Schoen-Yau [67] and Sampson [66] (In fact, harmonic diffeomorphism exist even without the curvature restriction, see [25,47,18].) A key ingredient in the proof are the following identities. We put, for a harmonic map φ : Σ → S between Riemann surfaces with curvatures K 1 , K 2 , Then, at points where H resp. L is nonzero, Actually, the Riemann-Hurwitz formula (17) can also be deduced from such identities, noting that for instance when h is holomorphic, L ≡ 0. Wolf [78] has developed a systematic calculus on the basis of these identities. The relationship between harmonic and ± holomorphic maps (or conformal maps, when the target is of higher dimension) is clarified by the following is a holomorphic quadratic differential. Furthermore, The proof of Lemma 3.1 is an easy computation. In the case where the target is also a Riemann surface S, that is, where we have a map φ : Σ → S, we have Then the computation for the proof of Lemma 3.1 becomes by (16).
We obtain, with the above notation λ 2 dz dz and ρ 2 dφ dφ for the metrics on Σ and S, (Of course, when we recall the definition (19) of H and L, the domain metric λ 2 dz dz drops out. Also, in our local coordinates, T (z) dz 2 = ρ 2 φ z φ z dz 2 .) Thus, the quadratic differential T (z) dz 2 is the (2, 0)-part of the pullback of the image metric.
In intrinsic terms, a (holomorphic) quadratic differential is a (holomorphic) Every holomorphic quadratic differential on the two-sphere S 2 vanishes identically, or putting it differently, the canonical bundle of S 2 admits no holomorphic section. Therefore, Lemma 3.1 implies In the next section, we shall see that this corollary reflects the fact that the conformal structure of S 2 admits no deformations.

Harmonic maps and Teichmüller theory
We consider Riemann surfaces of genus p ≥ 2 and the deformations of their complex structure. For Riemann surfaces of genus 0, there is only one complex structure, that of CP 1 . For Riemann surfaces of genus 1, tori, there is a complex one-dimensional family of complex structures. The theory is much simpler than that for genus p ≥ 2.
As already mentioned, such a Riemann surface Σ of genus p ≥ 2 is a quotient H/Γ of the Poincaré upper half plane by a discrete group Γ ∈ PSL(2, R) of isometries with respect to the hyperbolic metric 1 y 2 dz ⊗ dz, with z = x + i y. Since Γ is isomorphic to the fundamental group π 1 (Σ), the Riemann surface is described by a faithful representation ρ of π 1 (Σ) in PSL(2, R). This leads to the approach of Ahlfors and Bers to Teichmüller theory. Representations that only differ by a conjugation with an element of PSL(2, R) yield the same conformal structure. Thus, we consider the space of faithful representations up to conjugacy. A representation can be defined by the images of the generators, that is, by 2p elements of PSL(2, R), and this induces a natural topology on the moduli space. Also, from an easy count, we see that the (real) dimension of the moduli space of representations of π 1 (Σ) in PSL(2, R) modulo conjugations is 6p − 6. This is the dimension of the moduli space M p of Riemann surfaces of genus p. Here, however, we have somewhat more structure than simply the conformal structure of the Riemann surface in question. In fact, the moduli space presently discussed of discrete, faithful representations of π 1 (Σ) in PSL(2, R) modulo conjugations yields the Teichmüller space T p , a simply connected singularity-free infinite cover of M p . Putting it the other way around, M p is obtained as the quotient of T p by the mapping class group Γ p . Γ p is the group of homotopy classes of positively oriented diffeomorphisms of the underlying surface Σ, that is, where Diff + stands for orientation preserving diffeomorphisms and Diff 0 for those that are homotopic to the identity of Σ.
As a moduli of representations of a discrete group, T p also acquires natural structures, like a differentiable and a complex one. It is diffeomorphic, but not biholomorphic to C 3p−3 . The complex structure of Teichmüller space was originally proposed by Teichmüller and investigated by Ahlfors [1] and Bers [2], see [63] for details.
The complex dimension 3p − 3 of T p equals the dimension of the space of holomorphic quadratic differentials on a surface Σ of genus p. In fact, the holomorphic quadratic differentials on Σ are the cotangent vectors to T p at the point represented by Σ.
We shall now explore the geometry of T p more systematically, with the approach of Riemannian geometry. Since a two-dimensional Riemannian manifold defines a conformal structure, that is, a Riemann surface, we can naturally look at all Riemannian metrics on a given compact surface F of genus p and then identify those that induce the same conformal structure. We shall describe here the approach of Tromba and Fischer, see [73]. Thus, we consider the space R p of all smooth Riemannian metrics on F . R p , a space of Riemannian metrics, carries itself a Riemannian metric. Let g = g ij dx i ⊗ dx j (in local real coordinates (x 1 , x 2 )) be some Riemannian metric on F , that is, g is an element of R p . (g ij ) is a positive definite symmetric 2 × 2 tensor. Tangent vectors to R p at g then are given by symmetric 2 × 2 tensors (h ij ), ( ij ). The metric of R p on the tangent space at g is then given by A word of warning: This L 2 -metric on the infinite dimensional space R p is only a weak Riemannian metric, that is, the tangent spaces of R p are not complete w.r.t. this metric. Therefore, the general theory of Riemann-Hilbert manifolds does not apply. Clarke [17] showed that R p becomes a metric space with the distance function induced by this weak metric (23), and so, the situation is not as bad as one might fear.
In view of our discussion of the gravitino below in Section 7, we consider a Riemannian metric as the local difference between the Euclidean coordinate structure and the Riemann surface structure. As the uniformization theorem tells us, locally the difference can be made to vanish by choosing suitable coordinates. Globally, however, the metric on a compact Riemann surface cannot be made Euclidean, unless the surface is a torus, i.e., of genus one. Nevertheless, the metric carries some redundant information about the underlying conformal structure, as we shall now analyze.
Recall that the Riemannian manifolds (F, g 1 ) and (F, g 2 ) are called isometric if there is a diffeomorphism k : F → F , such that k * g 2 = g 1 . In that case we would like to consider the metrics g 1 and g 2 as equivalent. Conversely, any diffeomorphism k : F → F yields the obviously isometric Riemannian manifolds (F, g) and (F, k * g). This is not yet reflected in the definition of R p . We therefore divide out of R p the action of the (orientation preserving) diffeomorphism group D p of F . The diffeomorphism group D p acts isometrically on R p and hence the Riemannian metric (23) descends to the quotient of R p by D p .
However, even metrics that are not isometric might still induce the same conformal structure, and we therefore should also identify all such metrics. Such metrics differ by a some scalar factor. When we multiply a metric g by a positive function λ, the metric λg induces the same conformal structure as g, and conversely.
In order to eliminate the ambiguity of the conformal factor, one seeks a slice in R p transversal to the conformal changes. As explained above, by the Poincaré uniformization theorem, any Riemannian metric on a surface of genus p > 1 is conformally equivalent to a unique hyperbolic metric, that is, one with curvature −1, turning it into a quotient H/Γ. The moduli space M p is then obtained as the space R p,−1 of metrics of curvature −1 divided by the action of the diffeomorphism group D p . The geometric structures on R p then induce corresponding geometric structures on M p , see [73]. We can use the metric (23) to identify those directions on R p that are orthogonal to the conformal slice and the action of the diffeomorphism group. These are the ones that lead to nontrivial deformations of the complex structure. We recall that a tangent vector to R p is a symmetric 2 × 2 tensor (h ij ). Such a tensor is orthogonal to conformal multiplications when trace-free, and orthogonal to the action of D p when divergence-free. Trace-and divergence-free symmetric two-tensors are nothing but holomorphic quadratic differentials on the Riemann surface, and these are precisely the objects produced from harmonic mappings, see Lemma 3.1.
The Riemannian metric on R p then induces a Riemannian metric on the Teichmüller space T p . This theory is described in detail in [73]. And it follows from the preceding considerations that this metric in turn induces a product between holomorphic quadratic differentials on the Riemann surface Σ in question. In complex notation, with the hyperbolic metric on Σ denoted by λ 2 dz dz, let q 1 , q 2 be holomorphic quadratic differentials on Σ. Their product w.r.t. the metric then is given by This is the Weil-Petersson metric g W P . The Weil-Petersson is a Kähler metric w.r.t. the complex structure of the moduli space. This metric can be studied with the help of harmonic mappings. For instance, Tromba [74,75] found that the energy functional for harmonic maps between Riemann surfaces as a function of the conformal structure of the domain yields a strictly convex exhaustion function for the Weil-Petersson metric on T p . For some reviews of the geometry of the Weil-Petersson metric, see [79,80]. Tromba [72] proved that g W P has negative sectional curvature, and its holomorphic sectional curvature even has a negative upper bound k < 0, and this can again be derived with the help of harmonic mappings, see [27,45].
We now explain the harmonic map approach to Teichmüller theory in some more detail. We consider harmonic maps φ : Σ → S between compact Riemann surfaces of the same genus ≥ 2. In contrast to the domain Σ, for which we only need a conformal structure, the target S needs a Riemannian metric, and we can take the hyperbolic metric of constant curvature −1. As explained, we equip the surfaces with an additional structure, a representation of the fundamental group. Equivalently, we fix homotopy classes of diffeomorphisms. We then assume that our harmonic map respects these homotopy classes. In still other words, both Σ and S are considered to be diffeomorphic to some abstract underlying differentiable surface F , and we assume that the harmonic map φ is then homotopic to the identity of F . By Lemma 3.1, the harmonic map φ induces a holomorphic quadratic differential on Σ, that is, a holomorphic section of T * C Σ ⊗ T * C Σ. So, the harmonic map φ yields a cotangent vector of the Teichmüller space T p at the point corresponding to Σ. And this cotangent vector then varies when we vary the conformal structure of either Σ or S. Let us consider the effect of variations of the latter. From Theorem 3.1 and Lemma 3.1 and basic harmonic map theory, we have Lemma 4.1. Let Σ be a compact Riemann surface with local conformal coordinate z = x + i y, S another such surface of the same genus p ≥ 2, equipped with its unique hyperbolic metric, locally written as ρ 2 (h) dh dh. We also fix a homotopy class of (orientation preserving) diffeomorphisms k : Σ → S. Then there exists a unique harmonic diffeomorphism is a holomorphic quadratic differential, and φ is conformal iff T = 0.
We keep Σ fixed. The harmonic map φ then varies with S, and so then does the holomorphic quadratic differential on Σ. We obtain a map into the space of holomorphic quadratic differentials Q(Σ) on Σ. The map q is defined on Teichmüller space T p instead of on the moduli space M p because the harmonic map depends on the choice of a diffeomorphism k : Σ → S. Wolf [78] (see also [28]) then showed The proof makes crucial use of the identities (20), (21). Theorem 4.1 says that with its natural differentiable structure, T p is diffeomorphic to the vector space C 3p−3 . This is Teichmüller's theorem. As explained above, however, it is not biholomorphic to that space with respect to its natural complex structure. Also, q(Σ) is not an isometry w.r.t. the Weil-Petersson metric on T p . Nevertheless, there does exist a relationship between harmonic maps and the Weil-Petersson metric that we are now going to describe. For that purpose, we consider the energy as a function of S and present some formulas obtained in [78,28]. When S = Σ as an element of T p , that is, when the conformal structure of S is such that k is homotopic to a conformal diffeomorphism, and then the harmonic map φ is conformal. Therefore, L = 0, and the integral of the Jacobian, is a topological invariant that does not depend on the conformal structures of Σ and S, but only on their genus. Therefore, comparing with (25), E as a function of S achieves its minimum precisely when φ is conformal, that is, when the conformal structures of Σ and S are the same. At the point S = Σ ∈ T p , the infinitesimal variations of the target surface correspond to the cotangent vectors of T p at Σ, that is, to the holomorphic quadratic differentials on Σ. Since E achieves its minimum here, where E q = d dt t=0 E(φ(Σ, S + tq)) is the derivative of E for a variation of the target in the direction of q, and where we recall the identification of T p with Q(Σ) given by Theorem 4.1. Moreover, Lemma 4.2. S = Σ is the only critical point of E as a function of S ∈ T p . E(φ(Σ, ·)) is a proper function on T p . Its second derivatives at S = Σ for Consequently, E(φ(Σ, ·)) is a proper exhaustion function with a single critical point at which its Hessian is positive definite.
Equation (26) tells us that the second variation of the energy at S = Σ yields the Weil-Petersson product of the holomorphic quadratic differentials q 1 , q 2 . On this basis, one can develop systematic expansions [27].

Variations of conformal structures
We now look at the case where we vary the domain instead of the target. The corresponding calculations extend to the case where we have a general Riemannian manifold N as a target. This is useful when one studies minimal surfaces in Riemannian manifolds, see [27,29]. For simplicity, we assume here that all the harmonic maps into N that we consider here are unique in their homotopy classes, so that we don't have to deal with issues of bifurcation.
There is a new difficulty, however. When we vary the conformal structure of the domain, then z = x + i y no longer remain holomorphic coordinates. Thus, when we wish to express our formulas in the elegant complex notation, we need to account for that. Thus, instead of expressing things in non-holomorphic coordinates, we pull everything back from the varying surfaces to the fixed reference surface that we want to vary. And for our computations, we shall work with a variation of the domain metric instead of the conformal structure. While, as explained, many of the variations of the metric, more precisely those by diffeomorphisms or conformal factors, are redundant, insofar as they don't change the underlying complex structure, this is computationally easier. Thus, let γ = λ 2 dz ⊗ dz be a conformal metric on our domain Σ. And let γ t , t ∈ (−ε, ε) be a smooth variation of it. By what we have said before, we might assume that all the metrics γ t are hyperbolic metrics. Thus, the surfaces Σ t with the conformal structure determined by γ t (via Gauss' uniformization theorem) can be represented as quotients H/Γ t of the upper half plane. While this can help with the computations, here we only display the results, referring to [27], Sections 3.3 and 6.4, for details.
For a mapping u : Σ → N and a metric η on Σ, we let E(u, η) be its energy (14). Again, to repeat an important point, although in the second line of (14), the metric no longer appears explicitly, it is still implicitly contained, because the holomorphic coordinates z = x + i y depend on the metric.
Also, let u(γ t ) : Σ t → N be the harmonic maps that vary the harmonic map u(γ) : Σ → N (it follows from general principles that these maps depend smoothly on t). We let κ t : Σ → Σ t be diffeomorphisms, again depending smoothly on t, with κ 0 = id. Since u = u(γ) is harmonic, we consider u • κ −1 t : Σ t → N as an approximation of the harmonic u(γ t ). Since, as u 0 is harmonic, for any smooth variation u t of u = u 0 , this won't make a difference for the first order computations to follow. Indeed, with (27), we compute with T (z) dz 2 = g ∂φ ∂z , ∂φ ∂z dz 2 being the holomorphic quadratic differential associated with the harmonic map u, see (22), and with being the infinitesimal variation of the underlying structure. Now, transforms as a quadratic differential, and we may write (28) as This is the negative Weil-Petersson product (24) between the quadratic differentials T (z) dz 2 and Q(z) dz 2 . Here, T dz 2 is holomorphic when u is harmonic, but Q needs not be holomorphic. As we have discussed above, however, those variations µ in (29) that do not change the conformal structure lead to quadratic differentials Q dz 2 in (30) that are orthogonal to the holomorphic ones, and therefore, their Weil-Petersson products (31) with a holomorphic T dz 2 vanish. Conceptually, we can also consider the left hand side of (31) as the variation of γ in the direction Q dz 2 and write it as recalling again that for a harmonic u, we have ∂ ∂u E(u, γ) = 0. We can then also proceed to compute second derivatives when u(γ) is again assumed to be harmonic, and in fact, we get the analogue of (26), see [27], Section 6.4.

Dirac-harmonic maps and the gravitino
From the perspective of quantum field theory, the harmonic map and Riemann surface model is only the shadow of a much deeper and richer one, the supersymmetric nonlinear σ-model. The latter, while not directly a model of real physical forces, is one of the most basic and important models of quantum field theory, because of the symmetries it encodes and the conceptual and computational links with other models it provides. And even more fundamentally, it has become the basic action functional of superstring theory. For mathematical introductions, see [20,21,31]. As we have seen, the fundamental symmetry, conformal invariance, only holds when the domain is two-dimensional, that is, a Riemann surface. Another symmetry, supersymmetry, requires an additional field besides the harmonic map, a twisted spinor field along the map. It also requires that this field be anticommuting. This requires the introduction of additional anticommuting, so-called odd variables, in addition to the standard coordinate functions that we have been working with so far. This brings us into the realm of supergeometry, which has a number of features that appear odd from the usual perspective of differential geometry. In particular, for anticommuting variables, one cannot form inequalities, and therefore, the standard tools and estimates of PDE theory no longer apply. In order to prepare the ground, however, we shall first develop a mathematical version of the theory that need only ordinary, commuting variables and fields, introduced in [9,8]. This has since become a very active direction of research in the calculus of variations and geometric analysis, with many difficult and challenging problems that lead to the development of new and powerful mathematical tools that are quite useful also for other problems in geometric analysis. We now want to describe that theory.

Dirac-harmonic maps
(N, g) still is a Riemannian manifold. Let (Σ, γ) be a Riemann surface with a metric γ, as before, but now also equipped with a spin structure. Metrics and connections on the bundles appearing in the sequel will be those induced by γ and its Levi-Civita connection. The Spin(2) principal bundle corresponding to the spin structure is called P Spin (2) . Cl (Σ, −γ) is the corresponding Clifford algebra bundle; it is isomorphic to the quotient of the tensor algebra by the two-sided ideal generated by for X, Y ∈ Γ(T Σ). Its fiber is the Clifford algebra Cl 0,2 . We consider the spinor bundle S = P Spin(2) × Spin(2) Cl 0,2 where Spin(2) ⊂ Cl 0,2 acts via the left-regular representation of Cl 0,2 . We denote the Clifford multiplication of a vector X with a spinor s by X · s. For more details on the spin geometry and the spinor bundle, we refer to [36]. General background material can be found in [34].
We shall study the action functional A DH defined on the space where Γ(S ⊗φ * T N ) is the space of sections of the twisted spinor bundle S ⊗φ * T N . The action functional then is Here, the first term is the same as in (10). In the second term, / D is the twisted spin Dirac operator on S ⊗ φ * T N . In local coordinates y j on N , ψ = ψ j ⊗ φ * ∂ ∂y j and with a local orthonormal frame e α on Σ, it is given by where / ∂ is the ordinary spin Dirac operator on S (see for instance [34]). Furthermore, here and in the following we use the convention that we sum over repeated lower indices of the orthogonal frames e α . Since we choose the Clifford algebra Cl 0,2ˆ ψ, / Dψ S⊗φ * T N dvol =ˆ / Dψ, ψ S⊗φ * T N dvol (34) under integration by parts as well as by symmetry of the product ·, · S⊗φ * T N , that is, the corresponding term in the action functional is symmetric. Had we chosen Cl 2,0 , integration by parts would have introduced a minus sign, and the corresponding term in the action functional would have been zero. From (33), we see that / D depends not only on the spin geometry of the domain Σ, but also on the map φ. This leads to a coupling of the two fields φ and ψ in (32). Therefore, also the Euler-Lagrange equations for (32) are coupled. They are Here, τ (φ) is the tension field (12) of φ, the Dirac operator / D was defined in (33), and R N (φ, ψ) is a term depending on the Riemannian curvature 2 on N : These equations were derived in [9], and the investigation of the regularity of their solutions has been started in [8]. Further aspects and extensions of the model have been studied in [5,6,10,11,12,14,15,16,39,44,42,70,76,85,84] and many other papers. For the regularity theory, the method of Rivière [64] turned out to be very useful.
The existence theory met with the difficulty that, in contrast to the energy (10), the action (32) is not bounded from below. This comes from the term ψ, / Dψ , and ultimately from the fact that the spectrum of the Dirac operator / ∂, in contrast to that of the Laplace operator (5), is not bounded from below. Therefore, variational methods do not apply to show the existence of solutions. As already mentioned at the end of Section 2, an alternative to variational methods consists in heat flow methods. In the present case, however, there is the additional difficulty that the Dirac operator as a first order differential operator does not admit a natural parabolic version, again in contrast to the Laplace operator. Therefore, in [13], a novel elliptic-parabolic system has been introduced. The idea is to convert the second order elliptic system (35) into a parabolic system, but to carry the first order elliptic system (36) along as an elliptic side constraint. This leads to the system with suitable initial conditions, and also suitable boundary conditions when the domain has a boundary, see [16]. The short and long time existence of solutions, their regularity and their asymptotic behavior has been further investigated in [41,40,43,77].

The gravitino
In fact, the full supersymmetric action functional (which we shall discuss in Section 7 below) has a still richer structure and involves additional fields, in particular the gravitino field, a superpartner of the domain metric γ. Recently, in [36], we have constructed and investigated a version of this model that involves only commuting variables, analogous to and generalizing the Dirac-harmonic action function (32). For some further results on this model, see [52,51,37].
We shall now describe that model. 2 We use the sign conventions of [34], i.e., Definition 6.1. A gravitino χ is a section of the bundle S ⊗ T M .
For the moment, we assume χ to be smooth, but smoothness will become an issue below. In [35], a gravitino is introduced as a section of the bundle S ⊗ T * M ; we can, however, use the metric γ to identify T * M with T M .
The bundle S ⊗ T M can be decomposed into irreducible representations of Spin (2). In fact, we get two representations of type 1 2 and two of type 3 2 . We define the operator Q as the projection onto the 3 2 -parts. In local coordinates, it is given by We can now present the action functional, whose form, however, will only be derived and justified and whose content will only become clear in the next Section 7, with the curvature term Since Q is an orthogonal projection, And we observe that in the action functional (37), only the 3 2 -part Qχ of χ enters.
The functional A(φ, ψ; γ, χ) is of geometric interest because of its symmetries. Like the harmonic map energy (10), which it generalizes, it is invariant under generalized conformal transformations of the metric. However, as the spinor bundles depend on the metric we need to be a little careful. It is explained in [4] that for any two metrics γ and γ there is an isometric isometry b : (T Σ, γ) → (T Σ, γ ) that lifts to an isometry β : S → S of the spinor bundles. Using those maps, we obtain for γ = e 2u γ Here we use the scaling behavior of the twisted Dirac operator, / D e 2u γ e −u ψ = e −2u / D γ ψ. However, in our setting of conformal rescaling of the metric, the maps can be realized by b = e −u Id and β = e − 1 2 u Id which leads to the following more explicit invariance formula: A(φ, e −u ψ; e 2u γ, e −2u χ) = A(φ, ψ; γ, χ).
We call this invariance rescaled conformal invariance, since not only the metric but also the spinors are rescaled.
A further invariance comes from the super Weyl transformations, whose role will again become clear in the next Section 7, with Qχ = 0. As we have noted above, only the part Qχ enters into A, and since Q is a projection, we have Q 2 χ = Qχ.
As for the energy functional, we also have diffeomorphism invariance, again to be generalized in the next Section 7, and furthermore, A is Spin(2)-gaugeinvariant.
Although, as we shall see in Section 7, the action functional comes from that of the two-dimensional nonlinear supersymmetric sigma model, we are using here commuting spinors, and therefore (37) is in general no longer supersymmetric. Only in the case χ = 0 and under certain conditions on the curvature of the target, the functional A is supersymmetric also in the case of commuting spinors, see [37]. Thus, we have to pay a price for moving from anticommuting to commuting variables. That price is a loss of symmetry. Nevertheless, as described, our functional still possesses many symmetries, and these are crucial for its analysis and for its geometric content. Furthermore, as described in [58], the functional A is completely determined by the requirement of rescaled conformal and super Weyl invariance, given that the equations of motion are at most of second order.
From a somewhat lengthy computation, in [36], and adopting the notation there, we have obtained Theorem 6.1. The system of Euler-Lagrange equations for the action functional A is These equations already make the growth order transparent with which the various fields enter. The term SR(ψ), which has been defined in (38), is cubic in ψ; S∇R(ψ) involves derivatives of the curvature tensor R N and is quartic in ψ.
As for the harmonic map system, one can define weak solutions. Here, as before, φ needs to be of class W 1,2 , and the spinor ψ needs to be of class W 1,4/3 (this can be seen from the action functional and the Sobolev embedding theorem; the latter guarantees that such a ψ is, in particular, of class L 4 ). In [36], we have shown Theorem 6.2. Assume that the metric γ and the gravitino χ are smooth. Then the critical points of the action functional A are smooth.
Compare the discussion at the end of Section 2.1. Here, however, the regularity theory is much harder. A crucial ingredient is Rivière's method [64].
As before, the symmetries lead to conserved Noether currents, as we have already seen that conformal invariance of the harmonic map energy leads to a holomorphicity of the quadratic differential. Now in our situation here we have more symmetries, hence we may expect more conserved currents. They are actually given by the variation of the action functional with respect to the Riemannian metric and the gravitino respectively. Formally, we have In a local orthonormal frame (e α ), the energy-momentum tensor T is given by T = T αβ e α ⊗ e β with components while the supercurrent J is given by J = J α ⊗ e α with components J α = 2 φ * e β , e β · e α ψ + |ψ| 2 e β · e α · χ β .
The rescaled conformal symmetry prescribes the γ-trace of the energy-momentum tensor trace(T ). However, since the rescaled conformal invariance is not a pure conformal invariance, the trace will not be zero but rather depend on ψ and χ. The super Weyl symmetry prescribes the vanishing of the 1 2 -part of the supercurrent: (Id − Q)J = 0. As in the case of harmonic maps, the diffeomorphism invariance yields a conservation law of divergence-type: Here div χ is defined as the formal adjoint operator of the Lie derivative operator L S⊗T M , just as div γ is formally adjoint to the Lie derivative L on symmetric two tensors on Σ, see [73,37]. This is a law which involves the derivatives of the currents.
As explained above, we cannot expect full supersymmetry in this model. However, in the rare cases, that the action A is supersymmetric, we obtain another conservation law of divergence type. In that case the energy-momentum tensor can be considered as a holomorphic quadratic differential and the metric dual of the supercurrent as a holomorphic section of the 3 2 -part of S * ⊗T * Σ. This issue is analyzed in detail in [37]. The holomorphicity of T and J is somewhat surprising, but obtains a deep geometric explanation in the fully supersymmetric model below.
More generally, the scheme of converting a supersymmetric action principle of quantum field theory into a mathematical version with commuting fields works also in other cases, like the super Liouville action, and leads to a rich mathematical theory, see for instance [48,49,50,54,55,56].
This seems to open a beautiful research direction, with many difficult and challenging analytical problems, concerning the existence and regularity of solutions, and with the potential to profound geometric applications, possibly not unlike the Seiberg-Witten equations (see [68,69] and for instance [34] and [46] for the analytic aspects) or pseudoholomorphic curves giving rise to the Gromov-Witten invariants [24,71].

The supersymmetric action functional
We shall now switch from the Clifford algebra Cl 0,2 to Cl 2,0 . This will introduce a minus sign when we integrate by parts in (34). In order to compensate, for this minus sign, we now choose ψ as an anticommuting field. In physics terminology, ψ then becomes a fermionic field, as it should. And so will the gravitino field χ. In particular, this will allow us to gain an additional symmetry, supersymmetry. Mathematically, this means that we should enter the realm of supergeometry. This means that we can induce symmetries of the fields by transformations of the independent variables. Thus, we can capture supersymmetry as an extension or a partner of diffeomorphism invariance. And thus, the gravitino will become the superpartner of the metric, in the same manner as the spinor field ψ will become a superpartner of the bosonic field φ.

The action functional
In this section, we summarize results of [57,35]. The supersymmetric action functional that we want to investigate is This functional was introduced in early works on string theory, see [19,7], motivated by the search for a supersymmetric extension of the harmonic action functional. In order to realize supersymmetry, physicists always work with anticommuting variables. Indeed, it was argued in [58] that, while the action A can be motivated purely geometrically, the requirement of supersymmetry needs anticommuting spinors in general. Consequently, while A formally looks like A introduced in (37), we do now assume that ψ and χ are anticommuting fields and we shall work with the Cliffordalgebra Cl (Σ, γ) as in the physics literature. Those two changes together assure that the Dirac action is nontrivial, compare Equation (34).
As Cl 2,0 is isomorphic to the algebra of real 2 × 2-matrices, the corresponding spinor bundle S = P Spin(2) ×R 2 is of real rank two. The almost complex structure J Σ of T Σ can be extended to S by setting J S s = −ωs = −e 1 · e 2 · s for spinors s. The compatibility of the almost complex structures J S and J Σ yields J S (X · s) = (J Σ X) · s = −X · J S s for all vectors X and spinors s.

Consequently,
γ (Γ(s, t), X) = X · s, t S defines a complex linear isomorphism of line bundles Γ : S ⊗ S → T Σ. From now 3 on we assume that ψ is a section of S ⊗ φ * T N and χ a section of S ⊗ T Σ.
The bundle S ⊗ T Σ decomposes into irreducible Spin(2)-representations as S ⊕ S ⊗ C S ⊗ C S. We denote again by Q the projection on the 3 2 -part. Despite the change of conventions for the Clifford-algebra the functional A shares the invariances of A. The fact that we work with anticommuting spinor fields yields an additional invariance that gives rise to the so-called supersymmetry of the action. Thus, let us list the invariances of (39).

Diffeomorphisms of Σ:
4. Supersymmetry, that is first order invariance under: where q is a spinor, e α is a γ-orthonormal frame and ∇ S a spin connection with torsion.
Of course, all these symmetries induce Noether currents. In order to understand the action functional (39) and its symmetries in a systematic manner, we shall introduce the formalism of supergeometry, and more precisely, the notion of a super Riemann surface. In particular the local transformations (40), which look somewhat ad-hoc will be geometrically interpreted as particular diffeomorphisms of the underlying supermanifold. We shall then see that the same deep relationship that we have described in Sections 4, 5 between the action functional (10) and the geometry of the moduli space of Riemann surfaces extends to the present setting. In fact, the Noether currents for the symmetries of A will become cotangent vectors to the moduli space of super Riemann surfaces, and we can give a description of the latter in terms of metrics and gravitinos by dividing out symmetries.

Super geometry
We recall some notions from supergeometry. A supermanifold is a locally where Λ n is the real Grassmann algebra of n generators. The elements of Λ n anticommute and are called the odd coordinate directions. A supermanifold of dimension (m|0) is simply an ordinary manifold of dimension m, also called an even supermanifold, whereas one of dimension (0|n) is purely odd. We shall only be interested here in the case (m|n) = (2|2), that is, where we have two even and two odd dimensions. We shall also use a complex notation and write C 1|1 with coordinates z = x + i y (even) and ϑ = ϑ 1 + i ϑ 2 (odd). A function Φ : R 2|2 → R is then given by for even x = (x 1 , x 2 ) and odd η = (η 1 , η 2 ). Here, we already see the principle, that we consider ordinary (smooth) functions of the even variables x and expand in the odd variables η. We observe that because the coordinates of η anticommute, and therefore for instance η 1 η 1 = 0, there can be no further terms in the expansion (41). If we also want to have odd functions f in the expansion, we need to extend the base and consider, for some purely odd parameter domain B, where the f µ are odd and involve odd parameters of B. Colloquially speaking, we simply pull some odd parameters out of the hat whenever we need them. Thus, we leave B unspecified. We mostly consider a complex supermanifold of complex dimension (1|1). Coordinate transformations then are of the form More generally, when we consider a family over base B, we can also use odd functions of z and therefore consider the more general coordinate transformations with odd ξ, ζ. As explained, this is only meaningful if they involve an odd parameter from B.
One possible definition of the Berezin integral is to use the local formulâ and to globalize with the help of a partition of unity. Here f top is the coefficient of η n · · · η 1 , that is, of the highest order term occurring in the expansion of the function f in terms of the η i (recall again that they are anticommuting). Thus, each dη i in the integral simply cancels the coefficient η i in the expansion.

Super Riemann surfaces
A super Riemann surface is a (1|1)-dimensional complex supermanifold with an additional structure, a maximally non-integrable distribution D of the tangent bundle T M of rank (0|1), i.e., an isomorphism This is the description of LeBrun and Rothstein [60]. According to them, we can find local coordinates (z, ϑ) so that D is locally generated by D satisfies the key property which easily follows from the anticommutation rules for the odd variables. Thus, D is a square root of the ordinary derivative. Since the latter can be considered as the generator of spatial translations, that is, of a family of diffeomorphisms, D generates odd translations, or a family of superdiffeomorphisms. A general superdiffeomorphism then combines those two types of translations. Under the coordinate transformation (42), D is transformed to and for this in order to be equal to D , we need because then we get D = g(z) ϑ ∂ ∂z + ∂ ∂ϑ . This means that ϑ transforms as an (odd) section of the spinor bundle S = T |M | 1/2 . This in turn means that M corresponds to (|M |, S) where |M | is a Riemann surface and S is a spinor bundle. There exist 2 2p spin structures on a Riemann surface of genus p, one for each element of H 1 (|M |, Z 2 ), i.e., a choice of sign of the square root along 1-cycles. While this observation completely classifies trivial families of super Riemann surfaces, it is important to study non-trivial familes, in particular those where the base has odd dimensions.
An alternative description of super Riemann surfaces is due to Giddings and Nelson [23]. For them, a super Riemann surface is a (2|2)-dimensional real supermanifold with a reduction of the structure group to with certain integrability conditions. A further reduction is achieved by From that, we obtain a non-degenerate supersymmetric bilinear form on T M , given in U (1)-frames by where Latin indices a, b stand for the even, Greek indices α, β for the odd directions. δ ab is the symmetric Kronecker tensor, αβ the corresponding antisymmetric one. From (44), we get a split exact sequence The differential of the embedding i : |M | → M yields a different splitting of the pull back of (45) along i, and so, we have Here,p induces p s : i * T M → S. We have the identifications and the second equips T |M | with a metric γ, and the first makes S a spinor bundle for that metric, because by (44) The gravitino then is the section of T |M | ⊗ S given by Thus, analogously to our description in Section 4 of a Riemannian metric as encoding the difference between the local Euclidean coordinate structure and the global conformal structure of a Riemann surface, we see here that the gravitino encodes the difference between the local coordinate structure realized by the embedding i and the global super Riemann surface structure. In particular, in the same manner that a metric can be locally made Euclidean by the uniformization theorem, a gravitino can be locally gauged to 0. But, of course, not globally. Of course, since the gravitino here is an odd field, regularity issues as for the uniformization theorem do not arise here. On the other hand, from this perspective, we may ask what the appropriate regularity assumptions on the commuting gravitino field in Theorem 6.2 might be. This depends on whether also in that context, less regular fields can be transformed into more regular ones, by exploiting the redundancy contained in the gravitino field. The point here is that since the Riemannian metric contains redundant information about the underlying conformal structure, we also expect that its superpartner, the gravitino, contains some redundant information. We recall, for instance, that the action functional only involves the projection Qχ onto the spin 3/2 component. The metric, spinor bundle and gravitino entirely prescribe the super Riemann surfaces: Theorem 7.1. A (2|0)-dimensional Riemann surface |M | over a base B with a spinor bundle S, a metric γ and a gravitino χ determines a unique super Riemann surface M with embedding i : |M | → M , such that the above construction gives γ and χ back, up to Weyl (i.e., conformal) transformations of γ and super Weyl transformations of χ.
The proof depends on rather long and tedious computations employing a suitable adapted set of coordinates that were introduced by Wess and Zumino. This corollary is useful for deformation theory, because the right hand side involves no integrability conditions.
In fact, we even expect a global bijection The proof decomposes variations (h, r) of (γ, χ) as h = e u γ + L X γ + susy q γ + η r = χ + L X χ + susy q χ + ρ, where (e u γ, χ ) are infinitesimal (super) Weyl transformations, L X is the Lie derivative in the direction of the vector field X, and (susy q γ, susy q χ) is the supersymmetry transformation of (γ, χ). The remaining terms (η, ρ) are the true deformations of the super Riemann surface.

The symmetries of the action functional
We have already listed the symmetries of the action functional (39) above. With the geometric concepts developed, we can now write the action functional in a manner that makes the both the structure of this functional and the analogies with the functional (10) transparent. In fact, we can write where the integral on the right is a Berezinian (43), and the field Φ contains both φ and ψ as terms in its expansion with respect to the already mentioned Wess-Zumino coordinates Φ = φ + η µ ψ µ + η 1 η 2 F.
The Euler-Lagrange equations for the field F are F = 0, that is, F becomes zero on-shell. As the field F appears in the expansion of Φ, it is needed for the full superdiffeomorphism invariance of A. However, we do not enter this and simply omit the terms involving F . It should be noted that in (47) the action takes values in O B , hence a supercommutative ring possibly with nilpotent elements. This is necessary to allow a mathematical treatment of the anticommuting fields ψ and χ. As a consequence, the adaption of the analytical result for this variational problem requires further development of algebraic techniques.
The Berezin integral on the right-hand of (47) does not only look surprisingly similar to the harmonic action (10), but also its supergeometric properties resemble those of the harmonic action functional. Indeed, the Berezin integral does not depend on the chosen U (1)-reduction given by m, but rather only on the super Riemann surface structure. Every U (1)-reduction of the super Riemann surface will lead to the same value of the integral. This is called superconformal invariance of the integral. The superconformal invariance of the Berezin integral translates into the (super) Weyl invariance of A. Furthermore, the Berezin integral is superdiffeomorphism invariant. Every superdiffeomorphism of M can be decomposed into a diffeomorphism of |M | and a diffeomorphism of M that leaves i : |M | → M invariant. The latter ones induce supersymmetry transformations on the fields φ, ψ, γ and χ. Hence the supersymmetry invariance of A.
The superconformal and superdiffeomorphism invariance of A indicate that A should be a suitable tool to study the quotient of superconformal structures up to superdiffeomorphisms, i.e. the supermoduli space as in (46). This would be in parallel to the harmonic maps approach to ordinary Teichmüller theory.
As in Section 6.2, the Weyl invariance prescribes the trace of T and the super Weyl invariance prescribes J = QJ. Using the Euler-Lagrange equations for φ and ψ, the diffeomorphism invariance of A yields a coupled differential equation for T and J of divergence-type. A second coupled equation for T and J is obtained from supersymmetry.
In the case χ = 0 the equations decouple and give div γ T = 0 and div γ J = 0. In that case, being trace and divergence free, T can be considered as a holomorphic section of (T * M ) ⊗2 , and the metric dual of J as a holomorphic section of S * ⊗ C S * ⊗ C S * . Thus, according to Theorem 7.2, they constitute tangent vectors to the moduli space of super Riemann surfaces.
In conclusion, the action functional A reflects the structure of the moduli space of super Riemann surfaces by its symmetries and invariances in such a manner that its associated Noether currents become cotangent vectors to that space. Thus, we can represent the (infinitesimal) geometry of that moduli space analytically through a physical action functional. In turn, in order to understand the symmetries of that action functional, we can represent them as geometric variations of a superset of independent variables. A further understanding of the critical points is expected to lead to deep insights on the Teichmüller theory of super Riemann surfaces.