Minimal submanifolds from the abelian Higgs model

Given a Hermitian line bundle $L\to M$ over a closed, oriented Riemannian manifold $M$, we study the asymptotic behavior, as $\epsilon\to 0$, of couples $(u_\epsilon,\nabla_\epsilon)$ critical for the rescalings \begin{align*}&E_\epsilon(u,\nabla)=\int_M\Big(|\nabla u|^2+\epsilon^2|F_\nabla|^2+\frac{1}{4\epsilon^2}(1-|u|^2)^2\Big) \end{align*} of the self-dual Yang-Mills-Higgs energy, where $u$ is a section of $L$ and $\nabla$ is a Hermitian connection on $L$ with curvature $F_{\nabla}$. Under the natural assumption $\limsup_{\epsilon\to 0}E_\epsilon(u_\epsilon,\nabla_\epsilon)<\infty$, we show that the energy measures converge subsequentially to (the weight measure $\mu$ of) a stationary integral $(n-2)$-varifold. Also, we show that the $(n-2)$-currents dual to the curvature forms converge subsequentially to $2\pi\Gamma$, for an integral $(n-2)$-cycle $\Gamma$ with $|\Gamma|\le\mu$. Finally, we provide a variational construction of nontrivial critical points $(u_\epsilon,\nabla_\epsilon)$ on arbitrary line bundles, satisfying a uniform energy bound. As a byproduct, we obtain a PDE proof, in codimension two, of Almgren's existence result of (nontrivial) stationary integral $(n-2)$-varifolds in an arbitrary closed Riemannian manifold.


Introduction
A level set approach for the variational construction of minimal hypersurfaces was born from the work of Modica-Mortola [30], Modica [29], and Sternberg [34]. Starting from a suggestion by De Giorgi [12], they highlighted a deep connection between minimizers u ǫ : M → R of the Allen-Cahn functional and two-sided minimal hypersurfaces in M , showing essentially that the functionals F ǫ Γ-converge to ( 4 3 times) the perimeter functional on Caccioppoli sets. Several years later, Hutchinson and Tonegawa [19] initiated the asymptotic study of critical points v ǫ of F ǫ with bounded energy, without the energy-minimality assumption. They showed, in particular, that their energy measures concentrate along a stationary, integral (n−1)varifold, given by the limit of the level sets v −1 ǫ (0). These developments, together with the deep regularity work by Tonegawa and Wickramasekera on stable solutions [38], opened the doors to a fruitful min-max approach to the construction of minimal hypersurfaces, providing a PDE alternative to the rather involved discretized min-max procedure implemented by Almgren and Pitts ( [5], [31]) in the setting of geometric measure theory. This promising min-max approach based on the Allen-Cahn functionals was recently developed by Guaraco and Gaspar-Guaraco [16,14], and has been used successfully to attack some deep questions concerning the structure of min-max minimal hypersurfaces-most notably in Chodosh and Mantoulidis's work on the multiplicity one conjecture [11].
The initial motivation for this paper is to find, in a similar vein, a natural way to construct minimal varieties of codimension two through PDE methods. Recently, other attempts in this direction have been made by Cheng [10] and the second-named author [33], based on the study of the Ginzburg-Landau functionals on complex-valued maps v : M → C. While the Ginzburg-Landau approach can be employed successfully to produce nontrivial stationary rectifiable (n − 2)-varifolds (building on the analysis of [28], [8], and others), and leads to existence results of independent interest for solutions of the Ginzburg-Landau equations, it is not yet known whether the varifolds produced in this way are integral, nor is it known whether the full energies F ǫ (v ǫ ) of the min-max critical points converge to the mass of the limiting minimal variety in the case b 1 (M ) = 0. While it is possible that these and other technical difficulties may be overcome with sufficient effort-and establishing integrality in particular remains a fascinating open problem-they point to the deeper fact that the Ginzburg-Landau functionals, though intimately related to the (n − 2)-area, do not provide a straightforward regularization of the codimension-two area functional. Indeed, we stress that the Ginzburg-Landau energies should be understood first and foremost as a relaxation of the Dirichlet energy for singular maps to S 1 , and while the limiting singularities of critical points may coincide with minimal varieties, the associated variational problems exhibit substantial qualitative differences at both large and small scales.
In the present paper, we consider instead the self-dual Yang-Mills-Higgs energy while F ∇ ∈ Ω 2 (End(L)) denotes the curvature of ∇.
For the trivial bundle L = C × R 2 on the plane M = R 2 , a detailed study of the functional (1.1) and its critical points can be found in the doctoral work of Taubes [35,36]. In [36], all finite-energy critical points (u, ∇) of (1.1) in the plane are shown to solve the first order system 1 known as the vortex equations-a two-dimensional counterpart of the instanton equations in four-dimensional Yang-Mills theory. In particular, all such solutions (u, ∇) minimize energy among pairs (u, ∇) with fixed vortex number and carry energy exactly E(u, ∇) = 2π|N |. In [35], Taubes shows moreover that there exist solutions of (1.4) with any prescribed zero set which are unique up to gauge equivalence, so that [35] and [36] together give a complete classification of finite-energy critical points of (1.1) in the plane.
In [18], Hong, Jost, and Struwe initiate the study of the rescaled functionals (1.2) in the limit ǫ → 0 for line bundles L → Σ over a closed Riemann surface Σ. The main result of [18] shows that, for solutions (u ǫ , ∇ ǫ ) of the rescaled vortex equations (given by replacing 1 2 (1 − |u| 2 ) with 1 2ǫ 2 (1 − |u ǫ | 2 ) in (1.4)), the curvature * 1 2π F ∇ǫ converges as ǫ → 0 to a finite sum of Dirac masses of total mass | deg(L)|, away from which ∇ ǫ converges to a flat connection ∇ 0 , and u ǫ to a unit section u 0 with ∇ 0 u 0 = 0. While the authors of [18] focus on the vortex equations over Riemann surfaces, they suggest that the asymptotic analysis of the rescaled functionals E ǫ may also yield interesting results in higher dimension, pointing to similarities with the Allen-Cahn functionals for scalar-valued functions.
In the present paper, we develop the asymptotic analysis as ǫ → 0 for critical points of E ǫ associated to line bundles L → M over Riemannian manifolds M n of arbitrary dimension n ≥ 2. The bulk of the paper is devoted to the proof of the following theorem, which describes the limiting behavior as ǫ → 0 of the energy measures µ ǫ := 1 2π e ǫ (u ǫ , ∇ ǫ ) vol g and curvatures F ∇ǫ for critical points (u ǫ , ∇ ǫ ) satisfying a uniform energy bound.
Roughly speaking, Theorem 1.1 says that the energy of the critical points concentrates near the zero sets u −1 ǫ (0) of u ǫ as ǫ → 0, which converge to a (possibly rather singular) minimal submanifold of codimension two. In the case dim(M ) = 3, for instance, it follows from the results above and work of Allard and Almgren [3] that energy concentrates along a stationary geodesic network with integer multiplicities. The convergence of the curvature, moreover, to an integral cycle Poincaré dual to c 1 (L), with mass bounded above by lim ǫ→0 E ǫ (u ǫ , ∇ ǫ ), provides a higher dimensional analog to the limiting behavior described in two dimensions by Hong-Jost-Struwe [18].
At first glance, the obvious advantages of Theorem 1.1 over analogous results for the complex Ginzburg-Landau equations (cf., e.g., [8], [33]) are the integrality of the limit varifold V , and the concentration of the full energy measure to V , independent of the topology of M . Indeed, Theorem 1.1 and the analysis leading to its proof align much more closely with the work of Hutchinson and Tonegawa [19] on the Allen-Cahn equations than they do with related results (e.g. [27], [8]) for the complex Ginzburg-Landau equations. The parallels between the analysis presented here and that of the Allen-Cahn equations in [19] are in fact quite striking in places-a point to which we will draw the reader's attention throughout the paper.
Remark 1.2. We warn the reader, however, that while the qualitative analysis of the Allen-Cahn functionals does not depend on the precise choice of the double-well potential W , the analysis of the abelian Yang-Mills-Higgs functionals (1.1)-(1.2) seems to depend quite strongly on the choice W (u) = 1 4 (1 − |u| 2 ) 2 . Indeed, already in two dimensions, replacing W with a potential W λ (u) := λ 4 (1 − |u| 2 ) 2 for some λ = 1 yields a dramatically different qualitative behavior, breaking the symmetry which leads to the first-order equations (1.4), and introducing interactions between disjoint components of the zero set (see, e.g., [21,). This should serve as one indication that the analysis of the abelian Higgs model is somewhat more delicate than that of related semilinear scalar equations, in spite of the strong parallels.
To get some idea of the role played by gauge invariance, note that unit sections of a Hermitian line bundle are indistinguishable up to change of gauge (when no preferred connection has been selected), and for a given unit section u of L, one can always choose locally a connection with respect to which u appears constant. Thus, while most of the energy of solutions v ǫ to the complex Ginzburg-Landau equations falls on annular regions-relatively far from the zero set-where v ǫ resembles a harmonic S 1 -valued map, the energy e ǫ (u ǫ , ∇ ǫ ) of a critical pair (u ǫ , ∇ ǫ ) for the abelian Yang-Mills-Higgs energy instead concentrates near the zero set u −1 ǫ (0), with |∇ ǫ u ǫ | vanishing rapidly outside this region.
Of course, the results of Theorem 1.1 would be of limited interest if nontrivial critical points (u ǫ , ∇ ǫ ) could be found only in a few special settings. After completing the proof of Theorem 1.1, we therefore establish the following general existence result, showing that nontrivial families satisfying the hypotheses of Theorem 1.1 arise naturally on any line bundle (including, importantly, the trivial bundle) over any oriented Riemannian manifold M n , from variational constructions. Theorem 1.3. For any Hermitian line bundle L → M over an arbitrary closed base manifold M n , there exists a family (u ǫ , ∇ ǫ ) satisfying the hypotheses of Theorem 1.1, with nonempty zero sets u −1 ǫ (0) = ∅. In particular, the energy µ ǫ of these families concentrates (subsequentially) on a nontrivial stationary integral (n − 2)-varifold V as ǫ → 0.
For nontrivial bundles L → M , this follows from a fairly simple argument, showing that the minimizers (u ǫ , ∇ ǫ ) of E ǫ satisfy uniform energy bounds as ǫ → 0. For these energy-minimizing solutions, we expect moreover that the limiting minimal variety µ = θH n−2 Σ, i.e. the weight measure |V | of V , coincides with the weight measure |Γ| of the limiting (n − 2)-cycle Γ = lim ǫ→0 * 1 2π F ∇ǫ , and that Γ minimizes (n − 2)-area in its homology class. While we do not take up this question here, we believe that it would be very interesting to study the convergence of the functionals (1.2) to the (n − 2)area functional in a Γ-convergence framework. Let us mention that an asymptotic study for minimizers of the Ginzburg-Landau functional, on a domain with boundary, was successfully carried out by Lin and Rivière [27], who were able to identify the concentration measure with the weight of an integral current. (See also [1], [22] for related Γ-convergence results in that setting.) Remark 1.4. We remark that a very special class of minimizers for E ǫ are given by solutions (u ǫ , ∇ ǫ ) of the first-order vortex equations in Kähler manifolds (M 2n , ω K ) of higher dimension; these generalize the system (1.4) from the two-dimensional setting by replacing * F ∇ in (1.4) by the inner product F ∇ , ω K with the Kähler form ω K , and requiring additionally that F 0,2 ∇ = 0. As in the two-dimensional setting, solutions of this first-order system minimize the energy E ǫ in appropriate line bundles on Kähler manifolds, and it was shown by Bradlow 2 [9] that the moduli space of solutions corresponds to the space of complex subvarieties in M (of complex codimension one) via the zero locus (u ǫ , ∇ ǫ ) → u −1 ǫ (0). In particular, the zero loci u −1 ǫ (0) in this case are already area-minimizing subvarieties, before passing to the limit ǫ → 0. Note that the analysis of the vortex equations plays a key role in the study of Seiberg-Witten invariants of Kähler surfaces [39], and a similar analysis figures crucially into Taubes's work relating the Seiberg-Witten and Gromov-Witten invariants of symplectic four-manifolds [37]. For a concise introduction to the higher-dimensional vortex equations and connections to Seiberg-Witten theory, we refer the interested reader to the survey [13] by García-Prada. For the trivial bundle L ∼ = C × M , we prove Theorem 1.3 by applying min-max methods to the functionals (1.2), to produce nontrivial families (u ǫ , ∇ ǫ ) satisfying a uniform energy bound as ǫ → 0. While we consider only one min-max construction in the present paper, we remark that many more may be carried out in principle, due to the rich topology of the space where G := Maps(M, S 1 ) is the gauge group. Indeed, on a closed oriented manifold M , one can show that the homotopy groups π i (M) are given by it may be of interest to note that these are isomorphic to the homotopy groups of the space Z n−2 (M ; Z) of integral (n − 2)-cycles in M , as computed by Almgren [4].
As an application of Theorem 1.3, we obtain a new proof of the existence of stationary integral (n − 2)-varifolds in an arbitrary Riemannian manifold-a result first proved by Almgren in 1965 [5] using a powerful, but rather involved geometric measure theory framework. As already mentioned, similar constructions for the Allen-Cahn equations have been carried out successfully by Guaraco [16] and Gaspar-Guaraco [14], yielding new proofs of the existence of minimal hypersurfaces of optimal regularity, and leading to other recent breakthroughs in the min-max theory of minimal hypersurfaces (e.g., [11]).
In [11] and [16] (building on results of [38]), the stability properties of the min-max critical points for the Allen-Cahn functionals play a central role in controlling the regularity and multiplicity of the limit hypersurface. To obtain an improved understanding of min-max families (u ǫ , ∇ ǫ ) and the associated minimal varieties in the abelian Higgs setting, it would likewise be very interesting to refine the conclusions of Theorem 1.1 under the assumption that the families (u ǫ , ∇ ǫ ) satisfy a uniform Morse index bound as ǫ → 0. We hope to take up this line of investigation in future work.
1.1. Organization of the paper. In Section 2 we fix notation and record some basic properties satisfied by critical pairs (u ǫ , ∇ ǫ ) for the energies E ǫ .
In Section 3, we record some useul Bochner identities for the gauge-invariant quantities |u| 2 , |F ∇ | 2 , and |∇u| 2 , and use them to establish an initial rough estimate on , whose role should be compared to that of the discrepancy function in the Allen-Cahn setting. Under suitable assumptions on the curvature of M , the fact that ξ ǫ ≤ 0 follows quickly from the aforementioned Bochner identities and the maximum principle. Without the curvature assumptions, some nontrivial additional work is required to obtain the pointwise upper bound ξ ǫ ≤ C(M, E ǫ (u, ∇)). This estimate is the key ingredient to obtain the sharp (n − 2)-monotonicity of the energy.
In Section 4 we derive the stationarity equation for inner variations, from which an obvious (n − 4)-monotonicity property of the energy follows rather immediately. Using our rough initial bounds on ξ ǫ from Section 3, we deduce an intermediate (n − 3)monotonicity; we use this to reach the pointwise bound ξ ǫ ≤ C(M, E ǫ (u, ∇)), from which we finally infer the sharp (n − 2)-monotonicity.
In Section 5 we show that, similar to the Allen-Cahn setting, the energy density e ǫ (u, ∇) decays exponentially away from the set u −1 (0)-more precisely, away from Section 6, which constitutes the main part of the paper, contains an initial description of the limiting varifold, showing that it is stationary, (n − 2)-rectifiable, and has a lower density bound on the support. Then we establish its integrality with a blow-up analysis, employing the estimates from the preceding sections to reduce the problem to a statement for entire planar solutions, already contained in the work of Jaffe and Taubes [21]. We then use this analysis to show that the level sets u −1 ǫ (0) converge to the support of V in the Hausdorff topology, and conclude the section with a discussion of the asymptotics for the curvature forms 1 2π F ∇ǫ . In Section 7, we show that E ǫ satisfies a variant of the Palais-Smale property on suitable function spaces, allowing us to produce critical points via classical min-max methods. We provide a variational construction to get nontrivial critical points satisfying the assumptions of our main theorem, with energy bounded from above and below, both for nontrivial and trivial line bundles.
Finally, the Appendix addresses the issue of obtaining regularity of critical points, as obtained from Section 7, when they are read in a local or global Coulomb gauge.

Acknowledgements
A hearty thank-you goes to Tristan Rivière for introducing the authors to each other, and for suggesting the line of investigation taken up in the present paper. D.S. also thanks Fernando Codá Marques for his interest in this work, and Francesco Lin for pointing him to the reference [39]. A.P. is partially supported by SNSF grant 172707. During the completion of this work, D.S. was supported in part by NSF grant DMS-1502424.

The Yang-Mills-Higgs equations on U (1) bundles
Let M be a closed, oriented Riemannian manifold, and let L → M n be a complex line bundle over M , endowed with a Hermitian structure ·, · . Denote by W : L → R the nonlinear potential For a Hermitian connection ∇ on L, a section u ∈ Γ(L) and a parameter ǫ > 0, denote by E ǫ (u, ∇) the scaled Yang-Mills-Higgs energy where F ∇ is the curvature of ∇. Throughout, we will identify the curvature F ∇ with a closed real two-form ω via In computing inner products for two-forms, we follow the convention ω(e j , e k ) 2 (2.3) with respect to a local orthonormal basis {e j } n j=1 for T M . It is easy to check that the smooth pair (u, ∇) gives a critical point for the energy E ǫ , with respect to smooth variations, if and only if it satisfies the system Note that, in our convention, the adjoint to d : (∇ e j ω)(e j , e k ).
Since the curvature form ω is closed, taking the exterior derivative of (2.5) gives i.e., where ψ(u)(e j , e k ) := 2 i∇ e j u, ∇ e k u . For future reference, we record the simple bound To confirm (2.7), fix x ∈ M and note that the linear map ∇u(x) : T x M → L x has a kernel of dimension at least n − 2. Take an orthonormal basis {e j } of T x M such that e j ∈ ker ∇u(x) for j > 2. We compute at x that which gives (2.7).

Bochner identities and preliminary estimates
From the equations (2.6) and (2.4), we apply the standard Bochner-Weitzenböck formulas to obtain some identities which will play a central role in our analysis. For the curvature two-form ω, it will be useful to record the Bochner identity where R 2 denotes the Weitzenböck curvature operator for two-forms on the base Riemannian manifold M . For any δ > 0 we have Since |D|ω|| 2 ≤ |Dω| 2 , (3.1) implies Dividing by (|ω| 2 + δ 2 ) 1/2 and letting δ → 0, we obtain , in the distributional sense (and classically on {|ω| > 0}). Note that, by (2.7), the relation (3.2) also gives us the cruder subequation (3.3) ∆|ω| ≥ ǫ −2 |u| 2 |ω| − ǫ −2 |∇u| 2 − |R − 2 ||ω|. For the modulus |u| 2 of the Higgs field u, we record and observe that a simple application of the maximum principle yields the pointwise bound |u| 2 ≤ 1 on M. For the energy density |∇u| 2 of the Higgs field u, we see that Next, we introduce the function and combine (3.3) with (3.4) to see that From a simple application of the maximum principle, we see in particular that if R 2 > 0, then ξ ǫ ≤ 0 everywhere on M , and consequently (cf. [21, Theorem III.8.1]) This balancing of the Yang-Mills and potential terms, which should be compared with Modica's gradient estimate in the asymptotic analysis of the Allen-Cahn equations (cf. [19,Proposition 3.3]), will play a key role in our analysis, allowing us to upgrade the obvious (n − 4)-monotonicity typical of Yang-Mills-Higgs problems to the much stronger (n − 2)-monotonicity d dr (r 2−n Br e ǫ (u ǫ , ∇ ǫ )) ≥ 0. Without the positive curvature assumption, we may still employ the subequation ∆ξ ǫ ≥ |u| 2 ǫ 2 ξ ǫ − C(M )ǫ|F ∇ |, to obtain strong estimates for the positive part ξ + ǫ of ξ ǫ . To begin, denote by G(x, y) the nonnegative Green's function for the Laplacian on M , so that ∆ x G(x, y) = 1 vol(M ) − δ y , and set Taking C ′ to be the constant appearing in (3.7), for the difference ξ ǫ − C ′ h ǫ , we then have Observe that the L 1 norm of ξ ǫ − C ′ h ǫ is bounded by the energy: Thus, applying Moser iteration to the positive part (ξ ǫ − C ′ h ǫ ) + , we deduce that (Where the constant C(M ) may of course change from line to line.) As a simple application of (3.10), we note that by definition (3.8) of h ǫ and the standard estimate (see, e.g., [7,Chapter 4]) If n = 2, this inequality and (3.10) give a pointwise bound In the sequel, we assume n ≥ 3 and aim for a similar pointwise bound. We have Using this in (3.10), we compute at a maximum point for |F ∇ | to see that and, by an application of Young's inequality, it follows that for any δ ∈ (0, 1). Taking δ = ǫ 2/n , we arrive at the crude preliminary estimate where α(ǫ) → 0 as ǫ → 0. Now, consider the function By virtue of the preceding estimate for F ∇ L ∞ , we then see that pointwise. Appealing once again to (3.4) and (3.3), we see that ∆f ≥ |u| 2 ǫ 2 f − Cǫ|F ∇ |, so at a point where f achieves its maximum we have On the other hand, we know that |u| 2 ≥ ǫ C(1+Eǫ(u,∇) 1/2 ) f everywhere, so the preceding computations yield an estimate of the form and we deduce that f ≤ C(M, E ǫ (u, ∇)) everywhere. Putting all this together, we arrive at the following lemma.
In the next section, we will improve the rough preliminary estimate of Lemma 3.1 to a uniform pointwise bound of the form ξ ǫ ≤ C(M, Λ), but this will require some additional effort.

Inner variations and improved monotonicity
In this section, we derive the inner variation equation for solutions of (2.4)-(2.5), and explore the scaling properties of the energy E ǫ (u ǫ , ∇ ǫ ) over balls of small radius. Under the assumption that the curvature operator R 2 appearing in (3.3) is positive-definite (so that (3.6) holds), the analysis simplifies considerably, leading with little effort to the desired monotonicity of the (n − 2)-energy density. Without this curvature assumption, more work is required, first building on the cruder estimates of the preceding section to obtain a uniform pointwise bound for ξ ǫ .
Fixing notation, with respect to a local orthonormal basis {e i } for T M , define the (0, 2)-tensors ∇u * ∇u and ω * ω by Note that tr(∇u * ∇u) = |∇u| 2 and tr(ω * ω) = 2|ω| 2 . Denote by e ǫ (u, ∇) the energy integrand The fact that dω = 0 reads where D is the Levi-Civita connection of M . Using this identity, it is straightforward to check that In particular, defining the stress-energy tensor T ǫ (u, ∇) by for (u, ∇) solving (2.4) and (2.5) it follows that meaning that i (D e i T )(e i , ·) = 0. Integrating (4.4) against a vector field X on some domain Ω ⊆ M , we arrive at the usual inner-variation equation where we identify T ǫ (u, ∇) with a (1, 1)-tensor and denote by ν the outer unit normal to Ω. Taking Ω = B r (p) to be a small geodesic ball of radius r about a point p ∈ M , and taking X = grad( 1 2 d 2 p ), where d p is the distance function to p, (4.5) gives Now, by the Hessian comparison theorem, we know that applying this in the relations above, we see that it follows from the computations above (temporarily throwing out the additional nonnegative boundary terms) that At this point, one easily observes that the right-hand side of (4.7) is bounded below by n−4 r f (p, r), to obtain the monotonicity of the (n − 4)-energy density For general Yang-Mills and Yang-Mills-Higgs problems, this codimension-four energy growth is well known to be sharp (cf., e.g., [32], [40]). For solutions of (2.4) and (2.5) on Hermitian line bundles, however, we show now that this can be improved to (near-) monotonicity of the (n − 2)-density r 2−n f (p, r) on small balls, which constitutes a key technical ingredient in the proof of Theorem 1.1.
To begin, we rearrange (4.7), to see that recalling the notation ξ ǫ := ǫ|F ∇ | − 1 2ǫ (1 − |u| 2 ). Now, by Lemma 3.1, assuming E ǫ (u, ∇) ≤ Λ, we have the pointwise bound Applying this in our preceding computation for ∂f ∂r , we deduce that for some constant C ′′ (M, Λ) and 0 < r < c(M ). Taking ǫ sufficiently small, we arrive next at the following coarse estimate for the (n − 3)-energy density, which we will then use to establish an improved bound for ξ ǫ .
With Lemma 4.1 in hand, we can now improve the bounds of Lemma 3.1 to a uniform pointwise estimate, as follows.
Proof. We can assume n ≥ 3, as we already obtained the claim for n = 2 in Section 3.
Recall from that section the function where G is the nonnegative Green's function on M . As discussed in Section 3, we can deduce from (3.7) a pointwise estimate of the form Thus, to arrive at the desired bound (4.9), it will suffice to establish a pointwise bound of the same form for h ǫ .
To this end, recall again that G(x, y) ≤ C(M )d(x, y) 2−n , so that by definition we have where the last line is a simple application of Young's inequality. Since the integral On the other hand, by Lemma 4.1, we know that r 3−n Br(x) e ǫ (u, ∇) ≤ C(M, Λ) for every r, so we see finally that as desired.
Applying (4.9) in our original computation for f ′ (r), we see now that In fact, bringing in the extra boundary terms that we have been neglecting, and applying Young's inequality to the term r n−2 With this differential inequality in place, a straightforward computation leads us finally to one of our key technical theorems, the monotonicity formula for the (n − 2)-density.
As a simple corollary of the monotonicity result (together with a pointwise bound for |∇u| derived in the following section), we deduce that (u, ∇) must have positive (n − 2)-energy density wherever |u| is bounded away from 1.

Decay away from the zero set
In the preceding section, we obtained the pointwise estimate when ǫ ≤ ǫ m . As a first step toward establishing strong decay of the energy away from the zero set of u, we show in the following proposition that the full energy density e ǫ (u, ∇) is controlled by the potential W (u) ǫ 2 . Proposition 5.1. For (u, ∇) as above, we have the pointwise estimates Proof. To begin, let C 1 = C 1 (M, Λ) be the constant from (5.1), and consider the function Similar to the proof of Lemma 3.1, observe that C 1 |u| 2 ≥ f pointwise, by (5.1), while the computations from Section 3 give so that (max f ) 2 ≤ Cǫ, and consequently f ≤ Cǫ 1/2 everywhere. As a consequence, at any point, we have either f < 0, in which case In either scenario, we obtain a bound of the desired form (5.2).
To bound |∇u| 2 , recall from Section 3 the identity . In view of the estimate (5.1) for |F ∇ | = |ω| and (2.7), we can estimate the term 2 ω, ψ(u) from above by to obtain the existence of C 1 (M, Λ) such that For ∆|∇u|, this then gives Recalling once again the equation (3.4) for ∆ 1 2 |u| 2 , we define and observe that We then have If w has a positive maximum, it follows that at this maximum point; in particular, we deduce then that at this point, and see from (5.6) that here If ǫ ≤ ǫ d (M, Λ) is small enough, it follows that max w ≤ Cǫ; as a consequence, we check that As a simple consequence of the estimates in Proposition 5.1, we obtain the following corollary.
Corollary 5.2. There exist constants 0 < β d (M, Λ) < 1 and C(M, Λ) such that, for (u, ∇) as above, we have Proof. By the formula (3.4) for ∆ 1 2 |u| 2 , we know that Combining this with the estimate (5.3) for |∇u| 2 , we then deduce the existence of a constant C(M, Λ) such that By taking β = β(M, Λ) > 0 sufficiently small, we can arrange that we have an estimate of the form Proof. Fix a point p ∈ M , and let r = r(p) = dist(p, Z β ) as above. We can clearly assume r(p) < 1 2 inj(M ). On the ball B r (p), for some constant a = a d > 0 to be chosen later, consider the function . Now, fix some constant c 2 > 0 to be chosen later, and let Combining the preceding computation with (5.7), we see that, on B r (p), Choosing a = a d (M ) > 0 sufficiently small, we can arrange that 2a 2 + C 1 a ≤ 1, so that the above computation gives On the boundary of the ball ∂B r (p), it follows from definition of r = r(p) that |u| 2 ≥ 1 − β d , and therefore Taking c 2 (M, Λ) := β d e −ar/ǫ , it then follows that f < 0 on ∂B r (p), so we can apply the maximum principle with (5.9) to deduce that Evaluating at p, this gives Combining these estimates with those of Proposition 5.2, we arrive immediately at the following decay estimate for the energy integrand e ǫ (u, ∇).

The energy-concentration varifold
This section is devoted to the proof of the main result of the paper, which we recall now.
Let (u ǫ , ∇ ǫ ) be as in Theorem 6.1, and pass to a subsequence ǫ j → 0 such that the energy measures µ ǫ j converge weakly-* to a limiting measure µ, in duality with C 0 (M ).
With Proposition 6.2 in place, we will invoke a result by Ambrosio and Soner [6] to conclude that the limiting measure µ = lim ǫ→0 µ ǫ coincides with the weight measure of a stationary, rectifiable (n − 2)-varifold. Recall from Section 4 the stress-energy tensors We record first the following lemma; in its statement, we canonically identify (and pair with each other) tensors of rank (2, 0), (1, 1), and (0, 2), using the underlying metric g. Lemma 6.3. As ǫ → 0, the tensors T ǫ converge (subsequentially) as Sym(T M )-valued measures (in duality with C 0 (M, Sym(T M ))) to a limit T satisfying Proof. For each ǫ > 0, note that, by definition of T ǫ , for every continuous vector field Evaluating (2.3) in an orthonormal basis such that X is a multiple of e 1 , we see that As an immediate consequence, we see that the uniform energy bound E ǫ (u ǫ , ∇ ǫ ) ≤ Λ gives a uniform bound on T ǫ (C 0 ) * as ǫ → 0, so we can indeed extract a weak-* subsequential limit T ∈ C 0 (M, Sym(T M )) * , for which (6.7) follows from (6.8).
With this estimate in hand, we then see that Fixing δ and taking the limit as Finally, taking δ → 0, we conclude that M e ǫ (u ǫ , ∇ ǫ ) 1/2 → 0 as ǫ → 0, completing the proof.
Estimate (6.7) says that |T | is absolutely continuous with respect to µ, so by the Radon-Nikodym theorem we can write the limiting Sym(T M )-valued measure T from Lemma 6.3 as for some L ∞ (with respect to µ) section P : M → Sym(T M ). Moreover, it follows from (6.6) and (6.7) that −g ≤ P (x) ≤ g and tr(P (x)) ≥ n − 2 at µ-a.e. x ∈ M , so that T defines in a natural way a generalized (n − 2)-varifold in the sense of Ambrosio and Soner, namely a Radon measure on the bundle We refer the reader to [6,Section 3]. Note that in [6] the authors work in the Euclidean space and require the trace to be equal to n − 2 in (6.10); however, the main result on generalized varifolds, namely [6, Theorem 3.8], still holds in our setting (with the same proof).
Hence, in view of the stationary condition (6.5) and the density bounds of Proposition 6.2, we can apply [6, Theorem 3.8(c)] to conclude that T can be identified with a stationary, rectifiable (n − 2)-varifold with weight measure µ (so, in particular, spt(µ) is (n − 2)-rectifiable), and that P (x) is given µ-a.e. by the orthogonal projection onto the (n − 2)-subspace T x spt(µ) ⊂ T x M . We collect this information in the following statement.

Integrality of the limit varifold and convergence of level sets.
We now show that the varifold V is integer rectifiable. Given x ∈ spt(µ) and s > 0, we define M x,s to be the ball of radius s −1 inj(M ) in the Euclidean space (T x M, g x ) and define ι x,s : M x,s → M by ι x,s (y) := exp x (sy). We endow M x,s with the smooth metric g x,s := s −2 ι * x,s g, which converges locally smoothly to the Euclidean metric g x as s → 0. By rectifiability, for µ-a.e. x the dilated varifolds V x,s := (ι −1 x,s ) * (V B inj(M ) (x)) satisfy V x,s ⇀ v(T x Σ, Θ n−2 (x)) (6.12) as s → 0, in duality with C c (R n ). Fix x ∈ spt(µ) such that (6.12) holds. The integrality of V will follow once we prove that Θ n−2 (µ, x) is an integer.
We can identify (T x M, g x ) with R n by a linear isometry such that T x Σ = {0} × R n−2 . We also call µ x,s the mass measure of V x,s ; equivalently, µ x,s := s 2−n (ι −1 x,s ) * (µ B inj(M ) (x)).
With a diagonal selection, changing our sequence ǫ → 0 accordingly, we can find scales s ǫ → 0 such that we have the convergence of Radon measures where ( u ǫ , ∇ ǫ ) is the pullback of (u sǫǫ , ∇ sǫǫ ) by means of ι x,s , and µ ǫ is the associated energy measure. Note that ( u ǫ , ∇ ǫ ) is stationary for E ǫ in the line bundle ι * x,sǫ L, with respect to the base metric g x,sǫ . We introduce the notation Balls will be denoted by B r (y) or B n r (y), depending on whether they are with respect to g x,sǫ or g R n , respectively. The volume |E| of a set E will be always understood with respect to the Euclidean metric. The next proposition is the analogue of [26,Lemma 2.4] in this setting.

= lim
ǫ→0 (e C M sǫr r 2−n µ ǫ (B r (y)) + C M s ǫ r) → Θ n−2 (µ, x)ω n−2 . (6.13) Pick 3 ≤ i ≤ n and fix R > 0. Choosing y := −2Re i , we can apply (4.12) between the radii R and 3R to obtain that where p i := ι x,sǫ (−2Re i ) and ν R,i := grad d p i . Now (6.13) and the comparability of g x,sǫ with g R n give where ν R,i is the gradient of the distance function d −2Re i , both with respect to the metric g x,sǫ . Since eventually The smooth convergence g x,sǫ → g R n gives ν R,i (y) → Y R,i (y) := y+2Re i |y+2Re i | uniformly on B 2 2 × B n−2

2
. Hence, the bound (6.15) and (6.14) give Now Y R,i → e i = ∂ i as R → ∞, and the statement follows from (6.16) and the uniform bound (6.15).
Before giving the proof, let us see how this implies the integrality of V .

2
, R n−2 ) we can integrate (4.4) against χ( n i=3 Y i ∂ i ) and obtain, in the Euclidean metric, for some sequence λ ǫ → 0, thanks to the smooth convergence g x,sǫ → g R n . Invoking Proposition 6.5 and noting that Y L ∞ ≤ 2 DY L ∞ , we can conclude that the nonnegative function f ǫ (t) := 1 2π R 2 ×{t} χe ǫ satisfies , R n−2 ) (C 0 denoting the closure of C c ), we can find real measures (ν ǫ ) i j such that Since the sets F ǫ of Proposition 6.6 have positive measure, there clearly exists t ∈ F ǫ such that Recalling (6.17), we deduce that Hence, by (6.18), we get dist(Θ n−2 (µ, x), N) = 0, which concludes the proof that V is integral.

Corollary 5.4 gives
where we used Fubini's theorem in the second equality. The statement follows.
The following key technical proposition, used in the proof of Proposition 6.6, relies ultimately on the quantization phenomenon for the energy of entire solutions in the plane, presented in [21,Chapter III]. For the reader's convenience, we give a selfcontained proof, including the relevant arguments from [21]. Proposition 6.7. Given 0 < β, δ < 1 2 and S > 1, there exist K(β, δ, S) > S and 0 < κ(β, δ, S, n) < K −1 such that the following is true. Assume (u, ∇) is smooth and solves (2.4) and (2.5), with |u| ≤ 1 and ǫ = 1, on a line bundle L over a cylinder (Q, g), the energy bounds as well as the decay where p is the degree of u |u| (S·, 0), as a map from the circle to itself. Proof. To begin with, fix a real number K(β, δ, S) > S so big that ∞ K (2πr)Se −S −1 (r−S) < δ. (6.29) Arguing by contradiction, assume there exists a sequence κ j → 0 such that the statement admits a counterexample (u j , ∇ j ) (for κ = κ j ) for a (necessarily trivial) line bundle L j , with respect to a metric g = g j satisfying g − g R n C 2 ≤ κ j .
Fixing a trivialization of L j over Q j , we can write ∇ j = d − iA j for some real one-form A j .
By virtue of the uniform pointwise estimate (6.28) for e 1 (u j , ∇ j ) ≥ |d|u j || 2 , we see that the functions |u j | are locally equi-Lipschitz. In particular, we can apply the Arzelà-Ascoli theorem to extract a subsequence |u j | → ρ ∞ converging in C 0 loc . Since |d|u j || ≤ |∇ j u j |, (6.27) implies that ρ ∞ depends only on the first two variables. Moreover, (6.25) gives The degree p j is uniformly bounded as, for r ≥ S and t ∈ R n−2 , for j sufficiently large, so averaging over S < r < 2S and t ∈ B n−2 1 we get Thus, up to subsequences we can assume p j = p is constant.
We now claim that, up to change of gauge, (u j , ). Let u j = e iθ j u j be the section in the Coulomb gauge on the domain (B n 5S , g j ), with A j (ν) = 0 on the boundary (as described in the Appendix). Note that B n 5S includes the cylinder Q ′ := B 2 4S × B n−2

1
. and observe that, on Q ′′ : , u j has the form u j (re iθ , t) = |u j |e ipθ+iψ j for a unique real function ψ j with 0 ≤ ψ j (2S, 0) < 2π. Hence, u j = |u j |e i(pθ+ψ j −θ j ) on Q ′′ and we can extend ψ j − θ j uniquely to a function σ j : , for all 1 ≤ q < ∞, thanks to the Coulomb gauge specification (per Proposition A.1 in the Appendix.)

Moreover, in the exterior annular region
, we have that u j (re iθ , t) = |u j |e piθ and we can obtain local W 2,q bounds noting that Indeed, since the right-hand side is bounded by e 1 (u j , ∇ j ) 1/2 ≤ S 1/2 and pdθ is a fixed smooth one-form, we immediately obtain uniform L ∞ bounds for A j locally in A j . Next, note that the identity (3.4) applies to give us an estimate in A j , from which it follows that the modulus |u j | satisfies uniform W 2,q bounds for every q ∈ (1, ∞) locally in A j . Taking the imaginary part of (2.4) gives from which it follows that d * A j satisfies uniform L ∞ bounds locally in A j as well; together with the obvious pointwise bound |dA j | ≤ e 1 (u j , ∇ j ) 1/2 ≤ S 1/2 , this in particular yields uniform bounds on the full derivative DA j L q for every q ∈ (1, ∞) on fixed compact subsets of A j . Finally, writing (2.5) as the preceding chain of identities and estimates give a uniform L q bound on the righthand side over any fixed compact subset of A j , for any q ∈ (1, ∞); in particular, this gives us the desired uniform local W 2,q bounds for A j (while we already have the desired W 2,q bounds for u j = |u j |e piθ ).
Remark 6.8. As a consequence, one also finds that if |u| > 0 everywhere on the cylinder Q. Indeed, if |u| > 0 everywhere, then the degree p in the statement of Proposition 6.7 clearly must vanish.
We are now able to address the statement on the convergence of level sets. Proposition 6.9. For any 0 ≤ δ < 1 we have spt(µ) = lim ǫ→0 {|u ǫ | ≤ δ}, in the Hausdorff topology.
6.3. Limiting behavior of the curvature.
Since the two-forms ω ǫ are closed, for any ξ ∈ Ω n−3 (M ) we have so Γ is a cycle. By construction, Γ is Poincaré dual to c 1 (L).
To complete the proof, it remains to show that Γ has integer multiplicity. By means of a diagonal selection of a subsequence, as in the previous subsection, we can deduce integrality at those points p ∈ spt(µ) where µ blows up to Θ n−2 (µ, p)H n−2 T p Σ, using the following lemma. Note that its hypotheses are verified thanks to Corollary 5.4 and the fact that Z β d (u ǫ ) necessarily converges to T p Σ in the local Hausdorff topology, after rescaling (see the proof of Proposition 6.2).
Since µ is (n − 2)-rectifiable, we deduce that the limiting current Γ has integer multiplicity H n−2 -a.e. on its support, as claimed.
Lemma 6.11. On the Euclidean ball B n 4 , let (u ǫ , ∇ ǫ ) be a sequence of sections and connections in a trivial line bundle L → B n 4 (not necessarily satisfying any equation) By assumption, we then have Fixing trivializations of L over B n 2 , we write ∇ ǫ = d − iA ǫ for some one-forms A ǫ , so that ω ǫ = dA ǫ , and the right-hand term in the preceding limit becomes On B n 2 we can choose our trivializations so that d * A ǫ = 0, and A ǫ (ν) = 0 on ∂B n 2 . We then have the L 2 control (6.36) as ǫ → 0, where we have used the fact that dψ(x 1 , x 2 ) = 0 for |(x 1 , x 2 )| ≤ 1 2 , and the assumption that e ǫ (u ǫ , ∇ ǫ ) → 0 in C 0 loc (B n 2 \ P ).

Examples from variational constructions
The goal of this section is to show that, for every closed manifold M and every line bundle L → M endowed with a Hermitian metric, there exist critical couples (u ǫ , ∇ ǫ ) for the Yang-Mills-Higgs functional E ǫ , for ǫ small enough, in such a way that This will be easier when the line bundle is nontrivial, as in this case we can just take (u ǫ , ∇ ǫ ) to be a global minimizer for E ǫ . The upper and lower bounds in (7.1) have the following immediate consequence-proved previously by Almgren [5] using GMT methods.
Proof. We can always equip M with the trivial line bundle L := M × C. As shown in the next subsection, there exists a sequence of critical couples (u ǫ , ∇ ǫ ) satisfying (7.1). The statement now follows from Theorem 6.1.

Min-max families for the trivial line bundle.
In this section we will show how min-max methods may be applied to the functionals E ǫ to produce nontrivial critical points in the trivial bundle L = M × C on an arbitrary closed, oriented manifold M of dimension n ≥ 2. The min-max construction that we consider here is based on two-parameter families parametrized by the unit disk, similar to the constructions employed in [10] and [33] for the Ginzburg-Landau functionalswith several technical adjustments to account for the gauge-invariance and other features particular to the Yang-Mills-Higgs energies. We remark that the families we consider induce a nontrivial class in π 2 (M) for the quotient , ∇ a hermitian connection}/{gauge transformations}, and the analysis that follows can be reformulated in terms of min-max methods applied directly to the Banach manifold M.
Without loss of generality, we assume henceforth that M is connected.
In what follows, X will denote the Banach space of couples (u, A), where u ∈ L p (M, C) and A ∈ Ω 1 (M, R), both of class W 1,2 , with the norm Denote by X := {(u, A) ∈ X : d * A = 0} the subspace consisting of those couples for which the connection form A is co-closed.
Note that for (u, A) ∈ X, the full covariant derivative M |DA| 2 is bounded by C(M ) M (|A| 2 + |dA| 2 ). Note that the continuity of (u, A) → d(e iθ u) = e iθ (du + iudθ), from X to L 2 , follows from the fact that L p · L 2 * ⊆ L 2 , where 2 * = 2n n−2 . Throughout this section, W (u) = f (|u|) will be a smooth radial function given by W (u) = (1−|u| 2 ) 2 4 for |u| ≤ 3/2, and satisfying W (u) > 0 for all |u| > 1. For technical reasons, we also find it convenient to require that W (u) = |u| p for |u| ≥ 2, (G) which evidently gives the additional estimates |u|W ′ (|u|) + |u| 2 W ′′ (|u|) ≤ C|u| p for |u| ≥ 2, for some constant C. For future use, observe also that the potential W (u) then satisfies a simple bound of the form Proposition 7.5. The functional E ǫ is of class C 1 on X. Moreover, a couple (u, A) is critical in X for E ǫ if and only if R((u, A)) is critical in X.
Proof. Given a point (u, A) ∈ X and a pair (v, B) ∈ X with (v, B) X ≤ 1, direct computation gives where we are using the fact that X · X ⊆ L n · L 2 * ⊆ L 2 to see that , and we invoke our assumptions on the structure of W to see that for fixed (u, A) ∈ X. It follows immediately that E ǫ is C 1 on X, with gradient To confirm the second statement, assume without loss of generality that v and B are smooth, and observe that R((u + tv, A + tB)) = (e tiψ u + e iθ+tiψ v, A + tB + tdψ), where ( u, A) := R((u, A)) = (e iθ u, A + dθ) and ψ solves ∆ψ = d * B. This easily gives R((u + tv, A + tB)) = R((u, A)) + t(e iθ v + iψ u, B + dψ) + o(t) in X and, using the gauge invariance E ǫ = E ǫ • R, we deduce that We next show that the functionals E ǫ satisfy a suitable variant of the Palais-Smale condition on X, giving compactness of critical sequences for E ǫ after an appropriate change of gauge. (Cf. [23] for similar results in the Seiberg-Witten setting.) Proposition 7.6. The functional E ǫ satisfies the following form of the Palais-Smale condition: every sequence (u j , A j ) in X with bounded energy and dE ǫ (u j , A j ) → 0 in X * admits a subsequence converging strongly in X to a critical couple (u ∞ , A ∞ ), up to possibly replacing (u j , A j ) with Proof. First, by assumption, we have that ( (u j , 0) ).
The first term is bounded by 2E ǫ (u j , A j ), hence uniformly bounded as j → ∞. Moreover, it is clear from ( Denote by Λ ⊂ H 1 (M ) the lattice in the space of harmonic one-forms given by and let λ j ∈ Λ be a closest integral harmonic one-form to h(A j ) (with respect to the L 2 norm, say, on H 1 (M )). Then λ j = −v * j (dθ) for a suitable harmonic map v j : M → S 1 , and Replacing (u j , A j ) with the change of gauge (v j u j , A j − λ j ) ∈ X, we can then assume that h(A j ) is bounded. By standard Hodge theory we can write for some closed ξ j ∈ W 2,2 satisfying ∆ H ξ j = dA j and d * ξ j W 1,2 ≤ C(M ) dA j L 2 . Thus, given the energy bound E ǫ (u j ) ≤ C, we see that whereby A j is bounded in W 1,2 and, consequently, in L 2 * . As a consequence, we see next that taking into account (7.5), we infer then that du j L 2 is also bounded as j → ∞.
We have therefore shown that (u j , A j ) is uniformly bounded in X as j → ∞, so passing to subsequences we can assume that (u j , A j ) converges pointwise a.e. and weakly (in X) to a limiting couple (u ∞ , A ∞ ). In particular, defining r by where n < q < p is an arbitrary fixed exponent, it follows from the compactness of the embedding W 1,2 ֒→ L r that Moreover, the boundedness of u j in L p and the pointwise convergence to u ∞ give u j → u ∞ strongly in L q . (7.6) By definition of r, this implies in particular that lim j,k→∞ Next, compute and observe that the we then see that Now, by our assumptions (G) on the structure of W (u), it is not difficult to check (see, e.g., [17,Corollary 1]) that the zeroth order term in our computation for D j,k satisfies a lower bound for some constant C. In particular, it follows now from the preceding computations and the L 1 convergence u j → u ∞ that as j, k → ∞. On the other hand, since dE ǫ (u j , A j ) → 0 and (u j − u k , A j − A k ) is bounded in X, we know also that D j,k → 0 as j, k → ∞, and it then follows that (u j , A j ) is Cauchy in X. In particular, (u j , A j ) converges strongly to (u ∞ , A ∞ ), which necessarily satisfies Having confirmed that the energies E ǫ satisfy a Palais-Smale condition, we now argue in roughly the same spirit as [10], [33] to produce nontrivial critical points via min-max methods. To begin, note that the space X splits as C ⊕ Y , where C is identified with the set of constant couples (α, 0) and for all θ ∈ R. Equivalently, under the above identification C ⊂ X, we require F | ∂D = id. We denote by ω ǫ (M ) the "width" of Γ with respect to the energy E ǫ , namely Thanks to Proposition 7.6, we can apply classical min-max theory for C 1 functionals on Banach spaces (see e.g. [15,Theorem 3.2]) to conclude that ω ǫ is achieved as the energy of a smooth critical couple (u ǫ , A ǫ ). In the following proposition, we show that ω ǫ (M ) is positive, so that the corresponding critical couples (u ǫ , A ǫ ) are nontrivial.
Proof. We argue by contradiction, though the proof could be made quantitative. Since we are proving only the positivity ω ǫ (M ) > 0 at this stage-making no reference to the dependence on ǫ-in what follows we take ǫ = 1 for convenience. Assume that we have a family F ∈ Γ with max y∈D E(F (y)) < δ, with δ very small. Writing F (y) = (u, A), this implies that When b 1 (M ) = 0, some additional work is required to deduce that the harmonic part h(A) of A must also be small for all couples (u, A) = F (y) in the family. In particular, we will need to employ the following lemma, showing that h(A) lies close to the integral lattice Λ ⊂ H 1 (M ) when E(u, A) < δ. Lemma 7.9. There exists C(M ) < ∞ such that if (u, A) ∈ X satisfies E(u, A) < δ, then dist(h(A), Λ) ≤ Cδ 1/2 .
Proof. As in [33], it is convenient to define a box-type norm | · | b on the space H 1 (M ) of harmonic one-forms as follows. Fix a collection γ 1 , . . . , γ b 1 (M ) ∈ C ∞ (S 1 , M ) of embedded loops generating H 1 (M ; Q) and, for h ∈ H 1 (M ), set Since H 1 (M ) is finite-dimensional, this is of course equivalent to any other norm on H 1 (M ). Since M is orientable, we may fix a collection of diffeomorphims Φ i : B n−1 Suppose now that (u, A) ∈ X satisfies the energy bound As a consequence of the curvature bound dA L 2 ≤ δ 1/2 and the definition of X, it follows that A − h(A) 2 L 2 ≤ Cδ as well. As in the proof of Proposition 7.6, applying a gauge transformation φ · (u, A) by an appropriate choice of harmonic map φ : M → S 1 , we may assume moreover that which together with the energy bound (7.9) and the definition of X leads us to the estimate (Note that making a harmonic change of gauge preserves not only the energy E(u, A), but also the distance dist(h(A), Λ), so it indeed suffices to establish the desired estimate in this gauge.) Combining these estimates with a simple Fubini argument, we see that there exists a nonempty set S of t ∈ B n−1 1 for which Recalling the pointwise bound (7.2) for W (u), observe next that |d(1 − |u|) 2 | = 2(1 − |u|)|d|u|| ≤ CW (u) + |du − iAu| 2 , so that, along a curve γ t i satisfying (7.11), it follows that (7.14) (1 − |u|) 2 C 0 ≤ C (1 − |u|) 2 W 1,1 ≤ Cδ. Now, choose δ < δ 1 (M ) sufficiently small that (7.14) gives 1 − |u| C 0 ≤ η < 1 2 on γ t i , so that φ = u/|u| defines there an S 1 -valued map φ : γ t i → S 1 , whose degree is given by When (7.11)-(7.13) hold, we observe next that Since |u| ≥ 1 2 on γ t i , it follows that as well. Combining this with (7.12), we then deduce that On the other hand, we already made a gauge transformation so that so for δ chosen sufficiently small that Cδ 1/2 < π, it follows that the degree deg(φ) = 0.
In particular, we can now conclude that giving the desired estimate.
Returning to the proof of Proposition 7.8, suppose again that we have a family D ∋ y → F (y) ∈ X in Γ with max y∈D E(F (y)) < δ.
For δ < δ 1 (M ) sufficiently small, it follows from the lemma that dist b (h(A), Λ) < π for every couple (u, A) = F (y) in the family. In particular, since the assignment (u, A) → h(A) gives a continuous map X → H 1 (M ), and since h(A) = A = 0 for y ∈ ∂D, it follows that 0 is the nearest point in the lattice Λ to h(A) for every y ∈ D, and the estimate therefore becomes h(A) ≤ Cδ 1/2 . In particular, combining this with (7.7), we see now that (7.16) A W 1,2 ≤ Cδ 1/2 for every couple (u, A) = F (y) in the family. Now, for (u, A) = F (y), our structural assumption (G) on W (u) gives u p L p ≤ C + E(u, A) ≤ C + δ, which together with the smallness Combining this with the fact that |du − iuA| 2 ≤ E(u, A) < δ by assumption, we then deduce that M |du| 2 ≤ Cδ as well.
Finally, by (7.2) and the Poincaré inequality, we have As a consequence, we find that M u y is nonzero for all (u y , A y ) = F (y) in the family. But then the averaging map gives a retraction D → ∂D, whose nonexistence is well known. This gives the desired contradiction.
Having shown positivity ω ǫ (M ) > 0 of the min-max energies, we can now deduce the lower bound in (7.1) from the following simple fact.
Remark 7.11. For future reference, we make the obvious observation that the trivial case E ǫ (u, ∇) = 0 can only occur when the bundle L is trivial.
Proof. As discussed in the appendix, it is straightforward to see that critical points are smooth up to change of gauge. We claim that, whenever E ǫ (u, ∇) > 0, u has to vanish at some point x 0 ∈ M . Once we have this, Corollary 4.4 implies that r 2−n E ǫ (u, ∇, B r (x 0 )) has a lower bound independent of ǫ for any ǫ < r < inj(M ), and the statement follows.
Indeed, if the claim fails, then u is nowhere vanishing, so L must be trivial and we can use the section u |u| to identify L isometrically with the trivial line bundle M ×C, equipped with the canonical Hermitian metric. Under this identification, u : M → C takes values into positive real numbers. Writing ∇ = d − iA and observing that ∇u, iu = −|u| 2 A, (2.5) becomes Integrating against A we get M (ǫ 2 |dA| 2 + u 2 |A| 2 ) = 0, so A = 0 and ∇ is the trivial connection. At a minimum point y 0 for u, (3.4) gives , which forces u(y 0 ) ≥ 1 and thus u = 1 everywhere, giving the contradiction E ǫ (u, ∇) = 0.
Finally, we turn to the uniform upper bound. In the next statement, L → M is a Hermitian line bundle with a fixed Hermitian reference connection ∇ 0 . We identify any other Hermitian connection ∇ with the real one-form A such that ∇s = ∇ 0 s − iA ⊗ s for all sections s. Proposition 7.12. Given a smooth section u : M → L, we can find a smooth couple (u ′ , A ′ ) such that for a universal constant C.
Proof. On {u = 0} we let The compatibility of ∇ 0 with the Hermitian metric on L forces ∇ 0 w, w = 0, so that A is a real one-form. Equivalently, viewing w as a map from M to the circle bundle U (L) of L, which is a principal S 1 -bundle with induced connection form ̟ ∈ Ω 1 (U (L), R), Proof of (7.1). The method used in [33, Section 3] gives a continuous map H : D → W 1,2 ∩ C 0 (M, C) such that H(y) ≡ y for y ∈ ∂D and dH(y) L ∞ ≤ Cǫ −1 , {|H(y)|≤ 3 4 } |dH(y)| 2 ≤ C, vol |u| ≤ 1 2 ≤ Cǫ 2 (7.19) for all y ∈ D (the full Dirichlet energy having a worse bound M |dH(y)| 2 ≤ C log ǫ −1 , which is the natural one in the setting of Ginzburg-Landau). By approximation, we can assume that H takes values in C ∞ (M, C), continuously in y, and still satisfies the same uniform bounds (7.19) (possibly increasing C and replacing 3 4 with 1 2 ).
To each section H(y) of the trivial line bundle, Proposition 7.12 assigns in a continuous way an element F (y) ∈ X. From the way F (y) is constructed, it is clear that F ∈ Γ. Finally, applying Proposition 7.12 with (7.19) gives ω ǫ (M ) ≤ max y∈D E ǫ (F (y)) ≤ C.

Minimizers for nontrivial line bundles.
Suppose now that L is a nontrivial line bundle, equipped with a Hermitian metric. Fix a smooth Hermitian connection ∇ 0 and identify any other Hermitian connection ∇ with the real one-form A such that We can define X and X as in the previous subsection. With this notation, observe that the curvature of ∇ is given by Hence, writing F ∇ 0 = −iω 0 , we have Definition 7.13. For a fixed n < p < ∞, we define X to be the Banach space of couples (u, A), where u : M → L is an L p section and A ∈ Ω 1 (M, R), both of class W 1,2 , with the norm (u, A) := u L p + ∇ 0 u L 2 + A L 2 + DA L 2 .
The analogous statements to Remark 7.4 and Propositions 7.5 and 7.6 hold, with identical proofs (replacing du and uA with ∇ 0 u and A ⊗ u, respectively).
Arguing as in the proof of Proposition 7.6, it is easy to see that a minimizing sequence of couples for E ǫ converges-in the appropriate Coulomb gauge-to a global minimizer (u ǫ , A ǫ ). We now show that the energy of these minimizers enjoys uniform upper and lower bounds as ǫ → 0.
Proof of (7.1). The lower bound in (7.1) follows directly from Proposition 7.10. In order to obtain the upper bound, pick a smooth section s : M → L transverse to the zero section (see, e.g., [24, Theorem IV.2.1]) and let N := {s = 0}, which is a smooth embedded (n − 2)-submanifold of M . Proposition 7.12 applied to ǫ −1 s gives a couple (u ′ ǫ , A ′ ǫ ) with By transversality of s, the set {|s| ≤ ǫ 2 } is contained in a C(s)ǫ-neighborhood of N , whose volume is bounded by C(s)ǫ 2 . We infer that Remark 7.14. N can be oriented in such a way that [N ] ∈ H n−2 (M, R) is Poincaré dual to the Euler class e(L) ∈ H 2 (M, R) of the line bundle, which equals the first Chern class c 1 (L). The fact that the energy of our competitors concentrates along N suggests that, given a sequence of global minimizers (u ǫ , A ǫ ), up to subsequences the corresponding energy concentration varifold is induced by an integral mass-minimizing current whose homology class is Poincaré dual to c 1 (L). Theorem 6.10 provides the natural candidate Γ, which also satisfies |Γ| ≤ µ.

Appendix. Interior regularity in the Coulomb gauge
In this short appendix, we describe the essential ingredients needed to establish local regularity in the Coulomb gauge for finite-energy critical points (u, A) of the (ǫ = 1) abelian Higgs energy E(u, A), collecting some estimates which will be of use elsewhere in the paper.
Consider the manifold with boundary (Ω n , g) given by a smooth, contractible domain Ω n ⊂ R n equipped with a C 2 metric g, and let L ∼ = C × Ω be the trivial line bundle over Ω, with the standard Hermitian structure. With respect to the metric g, we then define the Yang-Mills-Higgs energies For the remainder of the section, we will assume that the pair (u, A) is already in the Coulomb gauge on Ω, so that A satisfies (A.6). Note that the last term in the equation (A.3) then vanishes, so that we have (A.7) ∆u = 2 idu, A + |A| 2 u − 1 2 (1 − |u| 2 )u.
Returning finally to the equation (A.2) for A, in the Coulomb gauge, we see that and it therefore follows from the preceding estimates that for some intermediate domain Ω ′ ⊂⊂ Ω ′′ ⊂⊂ Ω. In particular, this gives us upper bounds for ∆A L q (Ω ′′ ) for every q ∈ (1, ∞), and L q regularity theory therefore gives us the desired estimates for A in W 2,q (Ω ′ ).
Finally, we remark that higher regularity of u and A in the Coulomb gauge follows in a standard way-e.g., via Schauder theory-from the W 2,q estimates obtained in the preceding proposition.