Refined asymptotics of the Teichmüller harmonic map flow into general targets

The Teichmüller harmonic map flow is a gradient flow for the harmonic map energy of maps from a closed surface to a general closed Riemannian target manifold of any dimension, where both the map and the domain metric are allowed to evolve. Given a weak solution of the flow that exists for all time t≥0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t\ge 0$$\end{document}, we find a sequence of times ti→∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t_i\rightarrow \infty $$\end{document} at which the flow at different scales converges to a collection of branched minimal immersions with no loss of energy. We do this by developing a compactness theory, establishing no loss of energy, for sequences of almost-minimal maps. Moreover, we construct an example of a smooth flow for which the image of the limit branched minimal immersions is disconnected. In general, we show that the necks connecting the images of the branched minimal immersions become arbitrarily thin as i→∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\rightarrow \infty $$\end{document}.


Introduction
The harmonic map flow of Eells and Sampson [4] in two dimensions is the gradient flow of the harmonic map energy E(u) of a map u from a smooth closed oriented Riemannian surface (M, g) of genus γ to a smooth compact Riemannian manifold N = (N , G) of any dimension. Here the energy E(u) is defined by If we also allow the domain metric g to vary, and restrict g to have fixed constant curvature 1, 0 or −1 (and fixed area), then we end up with the Teichmüller harmonic map flow introduced in [11]. This is the flow given, for fixed parameter η > 0, by ∂u ∂t = τ g (u); ∂g ∂t = η 2 4 Re(P g ( (u, g))), (1.1) where τ g (u) := tr ∇du is the tension field of u, P g is the L 2 -orthogonal projection from the space of quadratic differentials on (M, g) onto the space of holomorphic quadratic differentials, and is the Hopf differential, written in terms of a local complex coordinate z = x + iy. See [11] for further information. Under the flow (1.1), ∂ t g is the real part of a holomorphic quadratic differential, and is therefore trace-free and divergence-free (see e.g. [21]). This implies that both the area form and the Gauss curvature are constant under the flow [22,Proposition 2.3.9]. For M = S 2 , the only holomorphic quadratic differential is identically zero, so ∂g ∂t ≡ 0 in this case and we are reduced to the classical harmonic map flow. In general, the flow decreases the energy according to The global existence theory for this flow, especially [2,12,15,16], leads naturally to the question of understanding the asymptotics of a given smooth flow as t → ∞, and it is this question that we address in this paper. Assume for the moment that M has genus γ ≥ 2, so g flows within the space M −1 of hyperbolic metrics (i.e. Gauss curvature everywhere −1). In the case that the length (g(t)) of the shortest closed geodesic in the domain (M, g(t)) is bounded from below by some ε > 0 uniformly as t → ∞ (which corresponds to no collar degeneration, even as t → ∞) it was proved in [11] that the maps u(t) subconverge, after reparametrisation, to a branched minimal immersion (or a constant map) with the same action on π 1 as the initial map u 0 . The same result also applies in the case of M being a torus with the length (g(t)) of the shortest closed geodesic bounded away from zero, as follows from the work in [2] when combined with the Poincaré inequality on nondegenerate tori proved in [14] (see also Remark 3.4).
On the other hand, on surfaces of genus γ ≥ 2 and in the case that lim inf t→∞ (g(t)) = 0, an initial description of the asymptotics of the flow was given in [13]. Loosely speaking, it was shown that the surface (M, g(t)) will degenerate into finitely many lower genus surfaces, with the map u(t) subconverging (modulo bubbling) to branched minimal immersions (or constant maps) on each of these components. It is convenient for us here to split the main result in [13] into an initial result extracting a sequence of times t i → ∞ at which the maps (u i , g i ) = (u(t i ), g(t i )) are a sequence of almost-minimal maps, and a separate compactness result for such sequences. For a justification of this terminology, note that using the construction in [11]-i.e. viewing the harmonic map energy E as a functional on the space of maps and constant curvature metrics modulo diffeomorphisms isotopic to the identity, equipped with the natural analogue of the Weil-Petersson metric-the gradient of E at points represented by (u i , g i ) converges to zero, while critical points are branched minimal immersions (see [11]). We are primarily concerned with the case that M has genus at least two, in which case each g i will lie in the space M −1 of hyperbolic metrics. (N , G), and a smooth flow (u, g) solving (1.1) for which lim inf t→∞ (g(t)) = 0, there exists a sequence t i → ∞ such that lim i→∞ (g(t i )) = 0 and (u(t i ), g(t i )) is a sequence of almostminimal maps.

Proposition 1.2 Given an oriented closed surface M, a closed Riemannian manifold
This standard observation is proved, for completeness, in Sect. 3. The sequence found in this proposition can be analysed with the following result, which is effectively what is proved in [13] (modulo minor adjustments to Sect. 3 of that paper). In this paper, we take this analysis of the asymptotics of (u i , g i ) and we refine it in several ways. First, after passing to a further subsequence, we extract all bubbles that can develop. What is well understood is that we can extract bubbles at each of the points in S (where possibly multiple bubbles can develop). In what follows we will call these bubbles {ω k }. Our first task is to isolate a new set of bubbles, called { j } below, that are disappearing into the 2k punctures found in Theorem 1.3, or equivalently (as we describe below), being lost down the one or more collars that degenerate in the domain (M, g i ) as i → ∞.
Having extracted the complete set of bubbles, we show that the chosen subsequence enjoys a no-loss-of-energy property in which the limit lim i→∞ E(u i , g i ) is precisely equal to the sum of the energies (or equivalently areas) of the branched minimal immersions found in Theorem 1. 3 Showing that no loss of energy occurs in intermediate regions around the bubbles developing at points in S is standard, following in particular the work of Ding and Tian [3] we describe in a moment, although one could also use energy decay estimates of the form we prove in this paper. However, showing that no energy is lost near the 2k punctures, away from where the bubbles develop, is different, and a key ingredient is the Poincaré estimate for quadratic differentials discovered in [14], which is applied globally, not locally where the energy is being controlled. In this step we exploit the smallness of P g ( (u, g)) that holds for almostminimal maps. That this is essential is demonstrated by the work of Parker [8] and Zhu [23], which established that energy can be lost along 'degenerating collars' in general sequences of harmonic maps from degenerating domains.
The following is the foundational compactness result when the domain is fixed, cf. [3,7,9,[18][19][20]. Theorem 1.5 Suppose (ϒ, g 0 ) is a fixed surface, possibly noncompact, possibly incomplete, and let u i be a sequence of smooth maps into (N , G) from either (ϒ, g 0 ), or more generally from a sequence of subsets ϒ i ⊂ ϒ that exhaust ϒ. Suppose that E(u i ) ≤ E 0 and that Then there is a subsequence for which the following holds true. There exist a smooth harmonic map u ∞ : (ϒ, g 0 ) → (N , G) (possibly constant) and a finite set of points S ⊂ ϒ such that we have, as i → ∞, N ). At each point in S, a bubble tree develops in the following sense. After picking local isothermal coordinates centred at the given point in S, there exist a finite number of nonconstant harmonic maps ω j : S 2 → (N , G), for j ∈ {1, . . . , J }, J ∈ N (so-called bubbles) which we view as maps from R 2 ∪ {∞} via stereographic projection, and sequences of numbers λ j i ↓ 0 and coordinates a j i → 0 ∈ R 2 , such that Moreover, we do not count bubbles more than once in the sense that
The bubbling has no energy loss in the sense that for each point x 0 ∈ S analysed as above, and each neighbourhood U ⊂⊂ ϒ of x 0 such that U ∩ S = {x 0 } only, we have Moreover, the bubbling enjoys the no-necks property Remark 1. 6 We note that the proof of the first part of Theorem 1.4 (virtually) immediately follows from Theorems 1.3 and 1.5: Away from S we can combine the strong W 1,2 -convergence of the maps with the convergence of the metrics. To analyse the maps U i = u i • f i near points in S we then apply Theorem 1.5 on small geodesic balls B which are of course isometric to one another provided r > 0 is chosen sufficiently small as the metrics g i are all hyperbolic. Finally, the convergence of the metrics allows us to relate the δ-thick part of ( , f * i g i ) to the δ-thick part of ( , h), compare [13,Lemma A.7], as well as the geodesic balls B This completes the proof of the first part of Theorem 1.4.

Remark 1.7
To do more, we must recall more about the structure of sequences of degenerating hyperbolic metrics, and in particular we need the precise description of the metrics g i near to the geodesics σ j i of Theorem 1.3 given by the Collar Lemma A.1 in the appendix. In particular, for δ ∈ (0, arsinh(1)) sufficiently small, the δ-thin part of (M, g i ) is isometric to a finite disjoint union of cylinders C δ, j i := (−X δ ( j i ), X δ ( j i )) × S 1 with the metric from Lemma A.1; each cylinder has a geodesic σ j i at the centre, with length j i → 0 as i → ∞. These initial observations motivate us to analyse in detail almost-harmonic maps from cylinders. Definition 1.8 When we apply Theorem 1.5 in the case that (ϒ, g 0 ) = R × S 1 is the cylinder with its standard flat metric, then we say that the maps u i converge to a bubble branch, and extract bubbles { j } as follows. First we add all the bubbles {ω j } to the list { j }. In the case that u ∞ : R × S 1 → N is nonconstant, we view it (via a conformal map of the domain) as a harmonic map from the twice punctured 2-sphere, remove the two singularities (using the Sacks-Uhlenbeck removable singularity theorem [17]) to give a smooth nonconstant harmonic map from S 2 , and add it to the list { j }. We say that u i converges to a nontrivial bubble branch if the collection { j } is nonempty.
We use the term 'bubble branch' alone to informally refer to the collection of bubbles together with the limit u ∞ .
In the present paper we prove a refinement of the above convergence to a bubble branch. To state this result, we shall use the following notations: For a < b, define C (a, b) := (a, b)× S 1 to be the finite cylinder which will be equipped with the standard flat metric g 0 = ds 2 + dθ 2 unless specified otherwise. For > 0 we write for short C = C (− , ). Furthermore, given sequences a i and b i of real numbers we write a i b i if a i < b i for all i ∈ N and b i − a i → ∞ as i → ∞. Theorem 1.9 Let X i → ∞ and let u i : C X i → N be a sequence of smooth maps with uniformly bounded energy, E(u i ; C X i ) ≤ E 0 < ∞, which are almost harmonic in the sense that (1.7) Then after passing to a subsequence in i, there exist a finite number of sequences s m i (for m ∈ {0, . . . ,m},m ∈ N) with −X i =: s 0 i s 1 i · · · sm i := X i such that the following holds true. for each m ∈ {1, . . . ,m}.

If we suppose in addition that the Hopf-differentials tend to zero
then there is no loss of energy on the connecting cylinders C (s m−1 i , s m i ) in the sense that for each m ∈ {1, . . . ,m} we have (1.10) Definition 1. 10 In the setting of Theorem 1.9, we abbreviate the conclusions of parts 1 and 2 by saying that the maps u i converge to a full bubble branch. In the case that (1.10) also holds (i.e. the conclusion of part 3) we say that the maps u i converge to a full bubble branch with no loss of energy.
Returning to the observations of Remark 1.7, we note that the length of each of the cylinders C δ, j i is converging to infinity, and that any fixed length portion of either end of any of these cylinders will lie within theδ-thick part of (M, g i ) for some smallδ ∈ (0, δ), and thus be captured by the limit u ∞ from Theorem 1.3. Our main no-loss-of-energy result can therefore be stated as the following result about the limiting behaviour on the middle of the collars, which constitutes our main theorem.

Theorem 1.11
In the setting of Theorem 1.3, we fix j ∈ {1, . . . , k} in order to analyse the jth collar surrounding the geodesic σ j i . Now that j is fixed, we drop it as a label for simplicity. Thus we consider the collar Then after passing to a subsequence, the restrictions of the maps u i to the collars C( i ) converge to a full bubble branch with no loss of energy in the sense of Definition 1.10.
We will furthermore analyse almost-minimal maps from degenerating tori in order to obtain the refined asymptotics for the flow (1.1) also in case that the genus of M is one.

Remark 1.12
Recall that any point in the moduli space of the torus can be represented uniquely as a quotient of (C, g eucl ) with respect to a lattice 2π,a+ib generated by 2π and a + ib ∈ C with −π < a ≤ π and |a + ib| ≥ 2π (with strict inequality if a < 0), b > 0. As a consequence, T 2 equipped with any flat unit area metric g is isometric to the quotient of (R × S 1 , 1 2π b g eucl ), with (s, θ) and (s + b, θ + a) identified. The length (g) of the shortest closed geodesic of (T 2 , g) is then given by Given a pair of sequences of maps u i : T 2 → N and of flat unit area metrics g i on T 2 , we let (a i , b i ) be as above so that ( We then lift the maps u i to maps defined on the full cylinder R × S 1 , also called u i and extend the notion of convergence to a full bubble branch as follows. Definition 1.13 Let (u i , g i ) be as above and assume that as We then say that the maps u i converge to a nontrivial full bubble branch if there exists a sequence s 0 i so that the maps u i (s 0 i + ·, ·) converge to a nontrivial bubble branch and so that the restrictions of u i (s 0 converge to a full bubble branch in the sense of Definition 1.10. Then there is a subsequence such that either (i) the maps u i converge to a nontrivial full bubble branch with no loss of energy or (ii) E(u i , g i ) → 0 and the whole torus is mapped near an i-dependent curve in the sense that lim sup i→∞ sup s∈[0,b i ] osc(u i ; {s} × S 1 ) = 0. In particular, we always have that for this subsequence where { j } is the collection of all bubbles.
As usual, this theorem can be applied to the flow (1.1), thanks to Proposition 1.2. For a discussion of existing claims concerning no-loss-of-energy for M = T 2 we refer to Remark 3.4. Both Theorems 1.11 and 1.14 indirectly describe the map u i on 'connecting cylinders' as being close to an i-dependent curve, thanks to (1.8). We are not claiming that this curve has zero length in the limit, as is the case in some similar situations, e.g. for necks in harmonic maps [8] and the harmonic map flow from fixed domains [9]. We are also not claiming that in some limit the curve should satisfy an equation, for example that it might always be a geodesic as would be the case for sequences of harmonic maps from degenerating surfaces, see [1]. The following construction can be used to show that these claims would be false in general. (N , G), a C 2 unit-speed curve α : [−L/2, L/2] → N , and any sequence of degenerating hyperbolic collars C X i , X i → ∞, equipped with their collar metrics g i as in the Collar Lemma A.1, the maps u i :

Proposition 1.15 Given a closed Riemannian manifold
We give the computations in Sect. 4. The proposition can be used to construct a sequence of almost-minimal maps with nontrivial connecting curves. For example, one can take any curve α as above, and any sequence of hyperbolic metrics g i with a separating collar degenerating, and then take the maps u i to be essentially constant on either side of this one degenerating collar where the map is modelled on that constructed in the proposition. A slight variation of the construction would show that the i-dependent connecting curve need not have a reasonable limit as i → ∞ in general, whichever subsequence we take, and indeed that its length can converge to infinity as i → ∞.
Finally, we consider the more specific question of what the connecting curves can look like in the case that we are considering the flow (1.1) and we have applied Proposition 1.2 to get a sequence of almost-minimal maps u(t i ) : (M, g(t i )) → N . One consequence of Theorem 1.9 is that for large λ and very large i, the restriction of u(t i ) to the connecting and in the images of branched minimal immersions that we have already found. (Note that it is not important to take θ = 0 in these limits. Any sequence θ i would give the same limits). Now that we have restricted to the particular case in which our theory is applied to the Teichmüller harmonic map flow, one might hope to rule out or restrict necks from developing. However, these necks do exist, and we do not have to have p m−1  16) i.e. at least one neck connects distinct points.
The simplest way of constructing an example as required in the theorem is to arrange that there can be no nonconstant branched minimal immersions in the limit, while preventing the flow from being homotopic to a constant map. The flow then forces a collar to degenerate in the limit t → ∞, and maps it to a curve in the target as we describe in Theorem 1.11. The precise construction will be given in Sect. 4. A key ingredient is the regularity theory for flows into nonpositively curved targets developed in [15].

Remark 1.17
It would be interesting to prove that in a large class of situations the connecting curves of Theorems 1.11 and 1.14, when applied to the Teichmüller harmonic map flow, will have a limit, and that that limit will necessarily be a geodesic. In the case that M = T 2 , and under the assumption that the total energy converges to zero as t → ∞, Ding et al. [2] proved that the image of the torus indeed converges to a closed geodesic.
The conclusion of the theory outlined above is a much more refined description of how the flow decomposes an arbitrary map into a collection of branched minimal immersions from lower genus surfaces.
As in [15], although we state our flow results for compact target manifolds, the proofs extend to some noncompact situations, for example when N is noncompact but the image of u(0) lies within the sublevel set of a proper convex function on N .
The paper is organised as follows: In the next section we derive bounds on the angular part of the energy of almost harmonic maps on long Euclidean cylinders. The main results about almost-minimal maps are then established in Sect. 3, where we first prove Theorem 1.9, which then allows us to show Theorem 1.11, and as a consequence to complete the proof of Theorem 1.4, and to prove Theorem 1.14. In Sect. 4 we prove the results on the images of the connecting cylinders stated in Proposition 1.15 and Theorem 1.16. In the appendix we include the statements of two well-known results for hyperbolic surfaces, the Collar lemma and the Deligne-Mumford compactness theorem, the statements and notations of which are used throughout the paper.

Angular energy decay along cylinders for almost-harmonic maps
Throughout this section we consider smooth maps u : is a compact Riemannian manifold that we isometrically embed in R N 0 and the cylinder C is equipped with the flat metric (ds 2 + dθ 2 ). The tension τ of u is given by Our goal is to prove a decay result for almost-harmonic maps from cylinders, forcing the angular energy to be very small on the middle of the cylinder C when we apply it in the setting of Theorem 1.11. This will be done by first controlling the angular energy, defined in terms of the angular energy on circles The proof of the following lemma is very similar to [20,Lemma 2.13], which in turn optimised [7]. Energy decay in such situations arose earlier in [9], and such results for harmonic functions are classical. More sophisticated decay results were required in [15]. Furthermore when 1 < λ < , we have the angular energy estimate

2)
and thus We require a standard 'small-energy' estimate, very similar to e.g. [
We now establish a differential inequality for ϑ(s). This is similar to [7, Lemma 2.1], but without requiring a bound on sup |∇u|. It is proved analogously to [20, Lemma 2.13], working on cylinders instead of annuli and considering a general target N .

Lemma 2.4
There exists a constant δ > 0 depending on N such that for u ∈ C ∞ (C (−1, 2), N ) satisfying E(u; C (−1, 2)) < δ and τ 2 L 2 (C (−1,2)) < δ, and for any s ∈ (0, 1), we have the differential inequality Proof of Lemma 2.4 From the proof of [15,Lemma 3.7] we have the expression We can estimate the penultimate term as in [15] using integration by parts and Young's inequality: while the final term of (2.4) requires just Young's inequality: where C is a constant only depending on N , that is revised at each step. Summing gives To apply Corollary 2.3, we can ask that δ < δ 0 , and thus handle the first term on the right-hand side as follows: and thus for δ sufficiently small, depending on N , we can improve (2.5) to It remains to estimate the inner product of u θθ with the tension in (2.4). By Young's inequality, we have by Wirtinger's inequality.
Lemma 2.4 can be applied all along a long cylinder C as arising in Lemma 2.1, and we can analyse the resulting differential inequality as in the next lemma to deduce bounds on ϑ. We then select A = 2E 0 e −S 2 and B = 2E 0 e S 1 to obtain such a solution for whichf ( . The maximum principle impliesf ≥ f and thus the claim.
We now apply the estimate from Lemma 2.5 to establish decay of angular energy.
Proof of Lemma 2.1 First note that we may assume that ≥ 1, otherwise the lemma is vacuous. By definition of ϑ, and the upper bound on the total energy, we have We choose δ > 0 smaller than both the δ of Lemma 2.4 and the δ 0 in Corollary 2.3.

Proofs of the main theorems; convergence to full bubble branches
Our main initial objective in this section is to prove Theorem 1.9, giving convergence of almost-harmonic maps to full bubbles branches. This will then be combined with the Poincaré estimate for quadratic differentials of [14] to give Theorem 1.11. Theorem 1.9 will also lead us to a proof of the analogous result for tori, Theorem 1.14. But we begin with a proof of Proposition 1.2.
Proof of Proposition 1.2 It is a standard idea (see e.g. [13,Section 3]) that by integrating (1.3) from t = 0 to t = ∞, and appealing to the boundedness of the energy, one can see that it is possible to extract a sequence of times t i → ∞ at which (u(t i ), g(t i )) is a sequence of almost-minimal maps. Note here that by a simple computation, we have We have to adjust the times t i to ensure that (g(t i )) → 0 as i → ∞ (by virtue of the hypothesis lim inf t→∞ (g(t)) = 0) and thus must argue that it is impossible for (g(t)) to spend almost all of the time away from zero, but drop quickly and occasionally down near zero. To do this, we pick timest i → ∞ witht i+1 ≥t i + 1 such that (g(t i )) ≤ 1/i for sufficiently large i, and claim that (g(t)) ≤ C/i for t ∈ [t i ,t i +1/i], with C depending only on the genus γ , the coupling constant η and an upper bound E 0 for the energy. If this claim were true then we would be able to pick our sequence of times t i from the set ∪ i [t i ,t i + 1/i], which has infinite measure, in the usual way. When the genus of M is at least 2, the claim follows from [15, Lemma 2.3], which implies in particular the Lipschitz bound whenever < 2 arsinh(1). Similarly, when M is a torus, we shall see in Appendix A.2 that Proof of Theorem 1.9 Let u i : C X i → N be a sequence of smooth almost harmonic maps as considered in Theorem 1.9. The first task is to construct sequences s m i as in the statement of the theorem. We would like to apply (2.3) on the regions C (s m−1 i + λ, s m i − λ) to the maps u i for large i, so we let δ > 0 be as in Lemma 2.1, which will be independent of i, of course.
We proceed to construct auxiliary sequencesŝ m i , where m ∈ {0, . . . ,m + 1} for somem ≥ 0. For each i, consider the overlapping chunks of length 3 of the form (k−1, k+2)×S 1 ⊂ C X i for k ∈ Z, i.e. for integral k such that −X i < k − 1 < k + 2 < X i . These chunks cover C X i except possibly for cylinders of length no more than 1 at the ends.
For each i, we initially choose the numbersŝ m i , for m = 1, 2, . . . , m i , to be the increasing sequence of integers so that (ŝ m i − 1,ŝ m i + 2) × S 1 are precisely the chunks above that have energy at least δ 2 . Note that by the bound on the total energy, there is a uniform bound on the number m i of such chunks, depending only on N and E 0 . Finally, we add inŝ 0 i = −X i andŝ m i +1 i = X i . By passing to a subsequence of the u i we can assume that for each i, we have the same number of sequence elementsŝ m i , i.e. m i =m for each i. Note also that for any region (s − 1, s + 2) × S 1 ⊂ C X i with energy at least δ there is some associated overlapping integer chunk (k − 1, k + 2) × S 1 ⊂ C X i of energy at least δ 2 which is assigned a label in the above construction, except possibly for regions very close to the ends of the cylinder in the sense that s − 1 < −X i + 1 or s + 2 > X i − 1.
From this auxiliary sequence we form s m i . Set s 0 i = −X i =ŝ 0 i , and consider the differencê s 1 i − s 0 i . If this has a subsequence converging to infinity, pass to that subsequence and take s 1 i =ŝ 1 i ; if not, discardŝ 1 i . Proceed iteratively to define s m i (i.e. s 2 i is the nextŝ m i such that the respective differenceŝ m i − s 1 i diverges for some subsequence, after having passed to that subsequence). This process will terminate with the selection of sm i , for somem. Whatever sequence sm i was chosen, redefine it as sm i = X i , which can only change it by an amount that is uniformly bounded in i. This finishes the construction.
By translation we can consider u i on C m i . We denote the shifted maps aŝ For each λ > 1, the estimate (2.3) from Lemma 2.1 applies (as in particular we have no concentration of energy) for sufficiently large i, giving Taking the limit i → ∞, and using the assumption (1.7) we find that Letting λ → ∞ proves that the 'no-loss-of-energy' claim (1.10) holds true provided the maps satisfy the additional assumption (1.9), which completes the proof of Part 3 of the theorem. We remark that this last step is the only part of the proof where (1.9) is used.
and thus we have Taking limits, and using once more (1.7), gives We conclude by observing that by the fundamental theorem of calculus and Cauchy-Schwarz, we can control the oscillation ofû m i on a circle in terms of the angular energy on that circle, by which implies the oscillation bound for u i claimed as Part 2 of the theorem.
The principle of the above proof also allows us to estimate the loss of energy on the connecting cylinders for sequences u i of almost harmonic maps that do not satisfy (1.9), as we shall require in the case M = T 2 .
while for m = 1 =m (respectively m =m = 1, respectively m = 1 =m) the analogue of (3.2) holds true with the integral over We remark that the reason for obtaining a weaker estimate in the case of m = 1 and m =m is that the corresponding cylinder does not connect two bubble branches but rather connects the ends of the cylinder (ifm = 1) or an end of the cylinder to a bubble branch (if m > 1). To prove this lemma we will use Remark 3.2 Recalling that bubble maps have vanishing Hopf differential (they are harmonic, so have holomorphic Hopf differential, which can only be identically zero on S 2 ) we see that in the setting of Theorem 1.5 the strong W 1,2 loc convergence of (1.6) and the fact that the bubbles separate as described in (1.5), tell us that In case of convergence to a bubble branch as in Definition 1.8, the map u ∞ also has vanishing Hopf differential (since it can also be extended to a harmonic map from the sphere) so this further reduces to We then note that Remark 3.2, applied to the shifted maps u m i (s, θ) = u i (s m i + s, θ), m ∈ {1, . . . ,m − 1}, yields that for every λ > 0 lim sup Combined with (3.3) this gives the claim.
The key step needed to derive Theorem 1.11 from Theorem 1.9 is to use the Poincaré inequality for quadratic differentials to get control on the Hopf differential.

Lemma 3.3
In the setting of Theorem 1.11, the Hopf differential decays according to Proof The Poincaré estimate for quadratic differentials [14] states that for any quadratic differential on the domain (M, g), and in particular for the Hopf differential , we have where C depends only on the genus γ ≥ 2 of M and is thus in particular independent of g. By (1.4), as the area of (M, g i ) is fixed, we know that and by direct computation (see e.g. [11, Lemma 3.2]) we know that where E 0 is an upper bound on the energies E(u i , g i ). Therefore, by (1.4) we find that ∂ (u i , g i ) L 1 → 0, and we conclude from (3.4) that Based on Lemma 3.3 and Theorem 1.9 we can now give the Proof of Theorem 1.11 To derive Theorem 1.11 from Theorem 1.9, we want to view the restriction of the maps u i to the collars C( i ) as maps from Euclidean cylinders C X i , X i = X ( i ), which are almost harmonic (with respect to g 0 ). We first remark that E(u i ; C X i ) is bounded uniformly thanks to the conformal invariance of the energy and the assumed uniform bound on E(u i , g i ). We then note that the conformal factors of the metrics ρ 2 (s)(ds 2 + dθ 2 ) of the hyperbolic collars (C( ), ρ 2 g 0 ), ∈ (0, 2 arsinh(1)), described in Lemma A.1 are bounded uniformly by Given that the norm of the tension scales as under a conformal change of the metric g = ρ 2 g 0 , where we continue to abbreviate τ (u) := τ g 0 (u) and equip C with the flat metric unless specified otherwise, we obtain from (1.4) that We furthermore note that the L 1 -norm of quadratic differentials is invariant under conformal changes of metric, compare (A. 3), and that the Hopf-differential depends only on the conformal structure. Lemma 3.3 thus yields Consequently all assumptions of Theorem 1.9, including (1.9), are satisfied and Theorem 1.11 follows.
Proof of Theorem 1.4 Continuing on from Remarks 1.6 and 1.7 it remains to analyse the energy on the degenerating collars C( j i ). After passing to a subsequence, Theorem 1.11 gives convergence to a full bubble branch without loss of energy on each of these collars so that the energy on δ-thin(M, Here we use that λ δ ( Proof of Theorem 1.14 Let (u i , g i ) be as in the theorem. As before we will think of (T 2 , g i ) as a quotient of (R × S 1 , 1 2π b i g 0 ) where (s, θ) ∼ (s + b i , θ + a i ) are identified and where we lift the map u i to the whole cylinder R × S 1 without changing the notation. Note that (g i ) → 0 implies b i → ∞. Since the maps u i satisfy (1.4) we have, by (3.5), The main difference compared with the proof of Theorem 1.11 is that we are not allowed to appeal to the Poincaré inequality since (3.4) is not valid with a uniform constant for tori.
Instead we shall use the fact that all holomorphic quadratic differentials on T 2 are of the form c · dz 2 , z = s + iθ, c ∈ C, so that the projection of the Hopf differential is As the tori (T 2 , g i ) have unit area, we furthermore know that dz 2 L 2 (T 2 ,g i ) = |dz 2 | g i = 4πb i , by (A.2). Thus (1.4) implies that (3.7) Let us first assume that, after passing to a subsequence, we can find a sequence s 0 i ∈ [0, b i ) so that the restrictions of the shifted maps u i (s 0 i + ·, ·) to C b i 2 converge to a nontrivial bubble branch in the sense of Definition 1.8. We can then analyse the restriction of the maps u i (s 0 i + 1 2 b i + ·, ·) to the cylinder C b i 2 with Theorem 1.9 and, after passing to a further subsequence, obtain that the maps u i converge to a nontrivial full bubble branch as described in Definition 1.13.
To compute the loss of energy on the connecting cylinders we apply Lemma 3.1 to conclude that holds true for m = 2, . . . ,m − 1. It also holds for m = 1 and m =m as we shall now explain.
In the general situation of Lemma 3.1 we did not obtain the above estimate for m = 1 and m =m because in that situation the cylinders C (s 0 i , s 1 i ) and C (sm −1 i , sm i ) do not connect two bubble branches but rather connect the boundary of the domain C X i to a bubble branch. In the present situation of M being a torus this issue does not arise and so we obtain (3.8) also for m = 1 and m =m.
To be more precise, if we extend u i periodically to a larger cylinder, say C s 0 i −b i ,s 0 i +2b i , and repeat the above argument then the cylinders C (s 0 i , s 1 i ) and C (sm −1 i , sm i ) obtained above appear as connecting cylinders between two bubble branches, so Lemma 3.1 implies that (3.2) is valid also for m = 1 and m =m.
The total amount of energy lost on the connecting cylinders is thus Assume instead that we cannot find, even for a subsequence, numbers s 0 i so that u i (s 0 i + ·, ·) converges to a nontrivial bubble branch. In this case we apply Theorem 1.9 twice on C 1 2 b i , once for u i itself and once for the shifted mapû i := u i (· + 1 2 b i , ·). The assumption means that the set of points s m i , m = 1, . . . ,m − 1, where nontrivial bubble branches develop, is empty, i.e.m = 1. Part 2 of Theorem 1.9 thus implies that the whole torus is mapped close to an i-dependent curve as described in the theorem. But Theorem 1.5 also implies that, after passing to a subsequence, the maps u i converge to a bubble branch, which must be trivial by assumption, and thus that E(u i ; C λ ) → 0 for every λ > 0. Combined with Lemma 3.1 and (3.7) we thus obtain Remark 3.4 When we apply Theorem 1.14 to the flow (1.1), using Proposition 1.2, we obtain a no-loss-of-energy result from (1.12) for M a torus. An assertion of this form was already made in [2], where a claim is made that the Hopf differential converges to zero in L 1 , based on the knowledge that its 'average' converges to zero. (Note that taking our viewpoint, the fact that the average of the Hopf differential converges to zero would be interpreted as the projection of the Hopf differential onto the space of holomorphic quadratic differentials converging to zero.) However, our take on the genus one case, in the absence of a uniform Poincaré estimate available in the higher genus case, is that we can only deduce that the L 1 norm of the Hopf differential converges to zero as a consequence of no-loss-of-energy, and along the way we need to appeal to the L 1 smallness of the Hopf differential of bubble branches from Remark 3.2.

Construction of a nontrivial neck
The main purpose of this section is to prove Theorem 1.16, but we first record the following elementary computations.
Proof of Proposition 1. 15 To ease notation, we drop all subscripts i for the following computations. We also simplify matters by embedding (N , G) isometrically in some Euclidean space and composing u with that embedding. The energy is conformally invariant, thus we calculate with respect to the flat metric Using (3.5) and the fact that ρ −1 ≤ −1 ≤ C X, which follows from Lemma A.1, we compute Finally, we compute the L 2 -norm of (u, g). With z = s+ iθ , (u, g) = |u s | 2 dz 2 . Recalling that |dz 2 | g = 2ρ −2 (see A.2) we similarly find The remainder of this section is devoted to the proof of Theorem 1.16, constructing a flow that develops a nontrivial neck. We opt for a general approach, although essentially explicit constructions are also possible. To this end, consider any closed oriented surface M of genus at least 2, and take the target N to be S 1 . Choose a smooth initial map u 0 : M → S 1 that maps some closed loop α on M exactly once around S 1 , and take any hyperbolic metric g 0 on M. We claim that the subsequent flow (1.1) develops a nontrivial neck.
The first key point is that since S 1 has nonpositive sectional curvature, the regularity theory from [15, Theorems 1.1, 1.2] applies, so the flow exists for all time.
The second key point is that because the target is S 1 , there do not exist any branched minimal immersions, except if one allows constant maps. In particular, no bubbles can form. If no collar degenerated in this flow, i.e. if there were a uniform positive lower bound for the lengths of all closed geodesics in (M, g(t)), then by the results in [11], the map u 0 would be homotopic to the constant map, which is false by hypothesis.
Therefore there are degenerating collars, and we can analyse them with Theorems 1.3 and 1.11 (using Proposition 1.2). We next demonstrate that a neck forms that is nontrivial in the sense that (1.15) holds. If not, then after passing to a subsequence, the maps from each degenerating collar would become C 0 close to constant maps. By [13, Theorem 1.1] this would imply that u 0 would be C 0 close to a constant map, and thus in particular it would be homotopic to a constant map, which again is false by hypothesis.
We have proved that our flow develops a neck in the sense that (1.15) holds for some degenerating collar, and some m. By Theorem 1.11 the image of the subcollar C (s m−1 i + λ, s m i − λ) will be close to a curve for each i. However, the limiting endpoints (1.13) and (1.14) of the curves may not be distinct, i.e. (1.16) might fail in general.
To make a construction in which (1.16) must hold for some collar and some m, it suffices to adjust our construction so that again all extracted branched minimal immersions must be constant, but so that the union of the images is not just one point. By our theory, the connecting cylinders will thus be mapped close to curves connecting these distinct image points so a nontrivial neck with distinct end points must develop.
To achieve this, we will lift the whole flow to a finite cover M of M. Given such a cover, we need to check that the lifted flow still satisfies (1.1). Locally, the lifting will not affect the tension field of u, so the lifted flow will satisfy the first equation of (1.1). However, care is required with the second equation since when we pass to the cover, new holomorphic quadratic differentials arise in addition to the lifts of the original holomorphic quadratic differentials. The following lemma will establish that this causes no problems.
where P g and P g are the L 2 -orthogonal projections onto the space of holomorphic quadratic differentials on (M, g) and (M, g), respectively.
As a consequence, given a solution (u, g) to (1.1) on M, and a covering map q as in the lemma, the lifted pair (u • q, q * g) will be a solution to (1.1) on M.
Proof If q is an n-fold cover, then the linear map from quadratic differentials on (M, g) to quadratic differentials on (M, g) given by is an L 2 isometry onto its image, because Thus it is clear that q * (P g ) coincides with the L 2 (M, g) projection of q * onto the space q * H, where H is the space of holomorphic quadratic differentials on (M, g). Therefore to establish (4.1), it remains to prove that if is a holomorphic quadratic differential on (M, g), then , q * = 0 for every quadratic differential ∈ H ⊥ . To see this, we consider the adjoint ϒ to the map → q * , which pushes down a quadratic differential on (M, g) to (M, g), adding up the n preimages. In particular, ϒ maps holomorphic quadratic differentials to holomorphic quadratic differentials, and so , q * = ϒ( ), = 0, as required.
Returning to our construction, we may assume that all the branched minimal immersions are mapping to the same limit point p (or we are done already). The aim is now to use a lifting construction, justified by the above, to obtain a lifted flow with the images of the corresponding branched minimal immersions being the two different lifts of p in a double cover of the target.
To this end, fix an arbitrary base point x 0 on M. Now consider the index 2 subgroup H of π 1 (M, x 0 ) consisting of loops whose images under the initial map u 0 go round the target S 1 an even number of times. We can pass to a (double) cover q : M → M of the domain satisfying q * π 1 M, x 0 = H (see e.g.  (u 0 , g 0 ). We can analyse this lifted flow using Theorem 1.11 for a subsequence of the times t i at which we analysed the flow on M.
It suffices to show that the images of the branched minimal immersions we can construct from the lifted flow consist of both of the lifts of p, not just one. These distinct points can then only be connected by nontrivial necks.
To see this, note that from the analysis of the original flow on M with Theorem 1.1 from [13] (i.e. with Proposition 1.2 and Theorem 1.3), we can find some δ > 0 sufficiently small such that for sufficiently large i, the δ-thin part of (M, g(t i )) will consist of a (disjoint) union of (sub)collars that eventually degenerate. For large enough i, the image of the δ-thick part of (M, g(t i )) will be contained in a small neighbourhood of p. For each such large i, we pick a point y i in the δ-thick part, and deform α to pass through y i . We view α then as a path that starts and ends at y i , and by assumption, the composition u 0 • α takes us exactly once round the target S 1 . In particular, as we pass once round the lift of α, we move from one lift of y i to the other, and the flow map moves from being close to one lift of p to being close to the other lift.
In particular, the branched minimal immersions in the lifted picture map to both lifts of p as required.
For δ ∈ (0, arsinh(1)), the δ-thin part of a collar is given by the subcylinder (−X δ ( ), X δ ( )) × S 1 ⊆ C( ), where X δ ( ) = 2π π 2 − arcsin sinh( 2 ) sinh δ (A.1) for δ ≥ /2, respectively zero for smaller values of δ. We make repeated use of the following differential geometric version of the Deligne-Mumford compactness theorem, describing how given a sequence of closed, oriented, hyperbolic surfaces (M, g i ), a subsequence converges to a complete, possibly noncompact, hyperbolic surface ( , h) with cusp ends. A consequence of the proposition is that if lim i→∞ inj g i M = 0 then we know that k ≥ 1 and is noncompact, but otherwise we may assume that k = 0 so = M.

A.2 Evolution of for M a torus
Here we provide a proof of the estimate (3.1) on the evolution of the length of the shortest closed geodesic of (T 2 , g) used in the proof of Proposition 1.2. Given a solution (u, g) of the flow (1.1) on M = T 2 , we let a, b be as in Remark 1.12 so that (T 2 , g(t)) is isometric to the quotient of (R × S 1 ), with (s, θ) ∼ (s + b, θ + a), |a + ib| ≥ 2π, equipped with 1 2π b g eucl . The space of holomorphic quadratic differentials on (T 2 , g) is spanned by dz 2 , z = s +iθ , so we can write We note that |dz 2 | is constant on T 2 , so (1.11) yields and (A.5) results in To estimate the change of along the flow we note that in the coordinates fixed above and for ∂ t g given by (A.4) We note that if |a + ib| > 2π then the only simple closed geodesics in (T 2 , g) of length are the curves {s} × S 1 . This allows us to choose a C 1 family of curves (σ ε ) ε∈(−δ,δ) with σ 0 = {0} × S 1 so that (t + ε) = L g(t+ε) (σ ε ). The resulting formula and (A.6) thus yield (3.1).
If |a + ib| = 2π then the closed geodesics which are represented by the lines s → (s, θ 0 + a b s) on R × S 1 have the same length as the curves {s} × S 1 so we cannot expect that the function t → (t) is differentiable. It is however still Lipschitz, with Lipschitz-constant bounded by 2π 2 |c+id| , as a simple modification of the above argument, carried out separately on the intervals [t, t + ε) and [t − ε, t], shows.