FROM THE N -CLOCK MODEL TO THE XY MODEL: EMERGENCE OF CONCENTRATION EFFECTS IN THE VARIATIONAL ANALYSIS
MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
Abstract. We investigate the relationship between the N -clock model (also known as planar Potts model or ZN -model) and the XY model (at zero temperature) through a Γ-convergence analysis as both the number of particles and N diverge. By suitably rescaling the energy of the N -clock model, we illustrate how its thermodynamic limit strongly depends on the rate of divergence of N with respect to the number of particles. The N -clock model turns out to be a good approximation of the XY model only for N sufficiently large; in other regimes of N , we show with the aid of cartesian currents that its asymptotic behavior can be described by an energy which may concentrate on geometric objects of various dimensions.
Keywords: Γ-convergence, XY model, N -clock model, cartesian currents, topological singularities.
MSC 2010: 49J45, 49Q15, 26B30, 82B20.
Contents 1. Introduction 1 2. Notation and preliminaries 12 3. Currents 15 4. Proofs in the regime ε θε ε| log ε| 29 5. Proofs in the regime θε ∼ ε| log ε| 66 6. Proofs in the regime ε| log ε| θε 1 68 7. Proofs in the regime θε ε 72 References 79
1. Introduction Classical ferromagnetic spin systems on lattices represent fundamental models to under- stand phase transition phenomena. On the one hand, the study of their properties has motivated the introduction of new mathematical tools which have provided useful insights for a number of problems arising in different fields. On the other hand, many techniques borrowed from probability theory, mathematical analysis, topology, and geometry have con- tributed to a better understanding of the properties of these systems. In this paper we make use of fine concepts in geometric measure theory and in the theory of cartesian currents to understand the relationship between the XY -model and the N -clock model (also known as planar Potts model or ZN -model) within a variational framework. The N -clock model is a two-dimensional nearest neighbors ferromagnetic spin model on the square lattice in which the spin field is constrained to take values in a set of N equi-spaced points of S1 . For N large enough, it is usually considered as an approximation of the XY (planar rotator) model, for which the spin field is allowed to attain all the values of S1 . The asymptotic behavior of the N -clock model for large N has been considered by Fr¨ohlich and Spencer in the seminal paper [30]. There the authors have proved that both the N -clock model (for N large enough) and the XY model present Berezinskii-Kosterlitz-Thouless transitions, i.e., phase transitions mediated by the formation and interaction of topological 1 2 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF singularities. The microscopic picture leading to the emergence of such topological phase transitions (first introduced in [14, 41, 42]) is a result of a nontrivial interplay between entropic and energetic effects that takes place at different length scales. This paper contributes to precisely relating the N -clock model and the XY model at zero temperature. Specifically, we show that the enhancement of symmetry, from the discrete one of the N -clock model to the continuous one of the XY model, comes along with concentration of energy on geometric objects of various dimension. This is achieved by studying several rescaled versions of the energy of the N -clock model as N diverges, through a coarse graining procedure which is made rigorous by Γ-convergence. A crucial step of this analysis is the choice of the topologies which best identify the relevant variables of the coarse grained model and lead to the effective description of the microscopic/mesoscopic geometry of the spin field. In contrast to the XY model, the sole study of the distributional Jacobian of the spin field turns out to provide not enough information on the concentration effects of the energy; we shall see how these effects can be detected by cartesian currents, for the first time introduced in the context of lattice spin models. In what follows we present the model and our main results. We consider a bounded, open set with Lipschitz boundary Ω ⊂ R2 . Given a small parameter ε > 0, we consider the 2 2 square lattice εZ and we define Ωε := Ω ∩ εZ . The classical XY model is defined on spin 1 fields u:Ωε → S by X − ε2u(εi) · u(εj) , (1.1) hi,ji where a · b denotes the scalar product between a, b ∈ R2 and the sum is taken over ordered 2 2 pairs of nearest neighbor hi, ji, i.e., (i, j) ∈ Z ×Z such that |i − j| = 1 and εi, εj ∈ Ωε . The variational analysis of the XY model is part of a larger program devoted to the study of systems of spins with continuous symmetry [4, 7, 8, 13, 23, 26, 27, 49, 25]. 2π Here we consider an additional parameter Nε ∈ or, equivalently, θε := . The N Nε 1 integer Nε represents the number of points in the discretization Sε of the codomain S , while θε is the smallest non-zero angle between two neighboring spin values. More precisely, we set
Sε := {exp(ιkθε): k = 0,...,Nε − 1} , where ι is the imaginary unit. The admissible spin fields we consider here are only those 1 taking values in the discrete set Sε , i.e., we consider the energy defined for every u:Ωε → S by X 2 − ε u(εi) · u(εj) if u:Ωε → Sε , Fε(u) := hi,ji +∞ otherwise.
For Nε = N ∈ N, with N independent of ε, the spin system described by the energy Fε is usually referred to as N -clock model, cf. [30]. The particular case where N = 2 is the so-called Ising system. The analysis of Ising-type systems with short-range interaction has been the object of many recent papers in analysis and statistical mechanics [38, 17, 24, 22, 3]. See also [2, 5, 9, 19] for the long-range case. The minimum of the functional defined in (1.1) is achieved on constant spin fields. For this reason it is customary in this setting to refer the energy to its minimum (cf. [4, Theorem 2]) and to introduce the functional X 1 X XY (u) = − ε2 u(εi) · u(εj) − 1 = ε2|u(εi) − u(εj)|2. ε 2 hi,ji hi,ji
This suggests to study the excess energy Eε(u) := Fε(u) − min Fε , under those scalings 1 κε → 0 for which Eε has a nontrivial variational limit in the sense of Γ-convergence, κε cf. [28, 18]. In order to understand such scalings, it is convenient to recast the energy as ( XYε(u) if u:Ωε → Sε , Eε(u) := (1.2) +∞ otherwise. FROM THE N -CLOCK MODEL TO THE XY MODEL 3
One can expect that the relevant scaling κε is affected by Nε , as it emerges in the 1 two limiting scenarios Nε = 2 and Sε = S (formally corresponding to Nε = +∞). In fact, for the Ising system, i.e., Nε = 2, it has been shown in [3] that, choosing κε = ε, a 1 configuration uε with equibounded energies ε Eε(uε) can be identified, as ε → 0, with a partition of Ω in two regions, in each of which the spin is constantly equal to one of the 1 two values (1, 0),(−1, 0). Moreover, the energy ε Eε(uε) approximates an interfacial-type energy, whose anisotropy reflects the symmetries of the underlying lattice. In contrast, for 1 the XY system, i.e., Sε = S , it has been shown in [4, Example 1] that no interfacial-type 2 energy emerges at any scaling κε ε . To show this, the authors provide an example that 1 we recall here. In this example, let Ω = B1/2(0) be the ball of radius 2 centered at 0, let 1 v1 = exp(ιϕ1), v2 = exp(ιϕ2) ∈ S , and let us define for x = (x1, x2) ( v if x ≤ 0 , u(x) := 1 1 (1.3) v2 if x1 > 0 . 1 1 1 2 Then we can construct uε → u in L (Ω; ) such that XYε(uε) → 0 for all κε ε in S κε the following way: for εi = (εi1, εi2) ∈ Ωε we define v1 if εi1 ≤ 0 , εi1 uε(εi) := exp ι (ϕ1 − ϕ2) 1 − + ϕ2 if 0 < εi1 ≤ ηε , (1.4) ηε v2 if εi1 > ηε , where ηε → 0 is chosen below. Then 2 1 ε ηε ε XYε(uε) ∼ 1 − cos (ϕ1 − ϕ2) ∼ . (1.5) κε ηε κε ηεκε 2 ε2 Thanks to the assumption κε ε , it is always possible to find ηε → 0 such that ηε . κε As a consequence, the angle between two neighboring spins uε(ε(i + e1)) and uε(εi) has ε κε 1 modulus |ϕ1 − ϕ2| . This construction, always possible in the case Sε = , may ηε ε S 1 κε not be feasible when Sε 6= S if the minimal angle θε satisfies θε & ε . Hence one cannot exclude the possibility of finding a nontrivial scaling κε such that interfacial-type energies 1 arise in the limit of Eε(uε). κε To better understand the behavior of Eε in the constrained setting uε :Ωε → Sε , we ε modify the construction in (1.4) by choosing ηε so that |ϕ1 − ϕ2| = θε , i.e., ηε = ηε ε |ϕ1 − ϕ2| , see Figure 1. As a result θε 1 ε εθε Eε(uε) ∼ 1 − cos(θε) |ϕ1 − ϕ2| ∼ |ϕ1 − ϕ2| . κε θεκε κε
This suggests that the nontrivial scaling κε = εθε leads to a finite energy proportional to |ϕ1 − ϕ2|. The construction can be optimized by choosing the angles ϕ1 and ϕ2 in 1 such a way that |ϕ1 − ϕ2| equals the geodesic distance on S between v1 and v2 , namely 1 dS (v1, v2).
v1 = exp(ιϕ1) θε v2 = exp(ιϕ2)
ε
ε ηε = |ϕ2 − ϕ1| θε
Figure 1. Construction which shows that 1 E approximates the geodesic distance between εθε ε the two values v1 and v2 of a pure-jump function. During the transition between v1 and v2 in the strip of size η = |ϕ − ϕ | ε the minimal angle between two adjacent vectors is θ . ε 1 2 θε ε 4 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
1 The fact that dS (v1, v2) is the total variation (in the sense of [11, Formula (2.11)]) of the 1 S -valued pure-jump function u defined in (1.3) suggests that at the scaling κε = εθε the 1 1 1 Γ-limit of Eε might be finite on the class BV (Ω; ) of -valued functions of bounded εθε S S variation. This intuition is confirmed in Proposition 4.1, where we prove that a sequence uε :Ωε → Sε (here and in what follows we identify it with its piecewise constant extension to Ω) satisfying the bound 1 Eε(uε) ≤ C εθε 1 1 is compact in BVloc(Ω; S ). More precisely, there exists a u ∈ BV (Ω; S ) such that, up 1 1 to a subsequence, uε → u strongly in L (Ω; S ). This compactness goes along with the 1 first result of the paper, Theorem 1.1 below, which states that Eε(uε) Γ-converges to a εθε functional that can be expressed on BV (Ω; S1) as an anisotropic total variation reflecting the symmetries of the underlying lattice. As it will be explained below, this result is true only in the regime ε| log ε| θε 1. To state the theorem, we need to introduce some notation, cf. [12]. Given a function u ∈ BV (Ω; S1), its distributional derivative Du can be decomposed as Du = ∇uL2 + (c) + − 1 2 D u + (u − u ) ⊗ νuH Ju , where ∇u denotes the approximate gradient, L is the Lebesgue measure in R2 ,D(c)u is the Cantor part of Du, H1 is the 1-dimensional Hausdorff 1 measure, Ju is the H -countably rectifiable jump set of u oriented by the normal νu , and + − u and u are the traces of u on Ju . By | · |1 we denote the 1-norm on vectors and by | · |2,1 the anisotropic norm on matrices given by the sum of the Euclidean norms of the columns.
Theorem 1.1 (Regime ε| log ε| θε 1). Assume that ε| log ε| θε 1. Then the following results hold:
i) (Compactness) Let uε :Ωε → Sε be such that 1 Eε(uε) ≤ C. εθε Then there exists a subsequence (not relabeled) and a function u ∈ BV (Ω; S1) such 1 2 that uε → u in L (Ω; R ). 1 ii)( Γ-liminf inequality) Assume that uε :Ωε → Sε and u ∈ BV (Ω; S ) satisfy uε → u in L1(Ω; R2). Then Z Z (c) − + 1 1 1 |∇u|2,1 dx + |D u|2,1(Ω) + dS (u , u )|νu|1 dH ≤ lim inf Eε(uε) . ε→0 εθε Ω Ju 1 iii)( Γ-limsup inequality) Let u ∈ BV (Ω; S ). Then there exists a sequence uε :Ωε → Sε 1 2 such that uε → u in L (Ω; R ) and Z Z 1 (c) − + 1 1 lim sup Eε(uε) ≤ |∇u|2,1 dx + |D u|2,1(Ω) + dS (u , u )|νu|1 dH . ε→0 εθε Ω Ju
The optimality of the regime ε| log ε| θε can be explained as follows. If θε ε| log ε| and uε :Ωε → Sε is a sequence with equibounded energies 1 Eε(uε) ≤ C, εθε then, by (1.2), we have
1 εθε 1 θε 2 XYε(uε) = 2 Eε(uε) ∼ ∼ 0 . (1.6) ε | log ε| ε | log ε| εθε ε| log ε| This entails that the formation of vortices is not allowed in this regime. Indeed, it has been 1 proven in [4] that the Γ-limit of the functionals ε2| log ε| XYε agrees with the total variation of a measure µ supported on finitely many points, usually referred to as vortices. (Such result drew inspiration from the literature on the Ginzburg-Landau functional, cf. [16, 46, 39, 40, 47, 1, 48]; a variational equivalence between the two models has been proven in [7].) FROM THE N -CLOCK MODEL TO THE XY MODEL 5
εi
µuε = δεi
2 Figure 2. Example of discrete vorticity measure equal to a Dirac delta on the point εi ∈ εZ . By following a closed path on the square of the lattice with the top-right corner in εi , the spin 1 field covers the whole S . The discrete vorticity measure can only have weights in {−1, 0, 1} .
1 A bound ε2| log ε| XYε(uε) ≤ C yields compactness for the discrete vorticity measure µuε associated to the spin field uε , where µuε counts the winding number of uε at each point of 2 f the lattice εZ (see Figure 2 and cf. (2.7) for the precise definition). More precisely, µuε * µ up to a subsequence in the flat convergence (i.e., in duality with Lipschitz functions with PN compact support, see (2.8)), where µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z. The measure µ represents the vortex-like singularities of the spin field uε as ε goes to zero. Moreover, the Γ-liminf inequality in [4, Theorem 1] yields 1 2π|µ|(Ω) ≤ lim inf XYε(uε) . (1.7) ε→0 ε2| log ε| The inequality above is optimal. This can be proven using the fact that the energy of 1 x−x0 the discretization uε :Ωε → of the prototypical vortex-like function , x0 ∈ Ω, S |x−x0| behaves as 1 XY (u ) ∼ 2π . (1.8) ε2| log ε| ε ε Combining inequality (1.7) with (1.6) we conclude that µ is the zero measure, i.e., f µuε * 0 . (1.9)
x−x0 This agrees with the fact that the the discretization vε :Ωε → Sε of carries an energy |x−x0| 1 E (v ) ∼ 2π , (1.10) ε2| log ε| ε ε cf. (4.78) for the rigorous bound. Notice that, in this case, the discretization vε is meant both in the domain and in the codomain. The emergence of the additional constraint on the absence of vortices in the limit suggests the possibility that the lower bound in Theorem 1.1- (ii) might be improved. This phenomenon does not occur in the regime ε| log ε| θε , where one has no information on the flat limit of µuε . Since smooth maps u:Ω → S1 provide the simplest example of S1 -valued maps with no 1 vortices, the problem of finding a better lower bound for Eε in the regime θε ε| log ε| εθε reminds that of finding the L1 -lower semicontinuous envelope of the functional with linear growth Z 1 1 1,1 1 |∇u| dx , u ∈ C (Ω; S ) ∩ W (Ω; S ) . (1.11) Ω Even though in this paper we are interested in anisotropic energies as the one obtained in Theorem 1.1, we notice here that they share the same features as (1.11). For this reason, not to overburden the reader with additional notation, in what follows we recall the analysis done for (1.11) in [33] (see also [34] for a generalization to manifold-valued maps). There the authors show that the functional Z Z (c) − + 1 1 1 |∇u| dx + |D u|(Ω) + dS (u , u ) dH , u ∈ BV (Ω; S ) ,
Ω Ju 6 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF is, in general, strictly smaller than the L1 -lower semicontinuous envelope of (1.11). For x example, let Ω = B be the ball of radius 1 centered at 0 and u(x) = |x| ; then it is x 1 1 not possible to find a sequence of smooth maps uk such that uk → |x| in L (Ω; S ) and R R x 1 Ω |∇uk| → Ω ∇ |x| . Nonetheless, it is possible to find a sequence of smooth functions x uk → |x| such that Z Z x |∇uk| dx → ∇ |x| dx + 2π , (1.12) Ω Ω where 2π is reminiscent of the energy due to the necessary concentration of ∇uk on a radius of Ω = B . The correct framework to carry out this analysis in its full generality is that of cartesian currents cart(Ω×S1) (cf. [31, 32]), which can be characterized as those currents obtained as weak limits of graphs of smooth maps from Ω to S1 (see Theorem 3.3). Given a sequence 1 1 1,1 1 1 ∗ 2 uk ∈ C (Ω; S ) ∩ W (Ω; S ) and u ∈ BV (Ω; S ) with uk * u in BV (Ω; R ) there exists a 1 T ∈ cart(Ω×S ) such that the graphs Guk of uk converge weakly (in the sense of currents) to T . The currents associated to the graphs of smooth maps are boundaryless. Such a property is inherited by the limit object T , i.e., ∂T = 0. For this reason, in general, the current T is different from the graph Gu of the limit map u, which may have a nontrivial boundary. Nevertheless, T can be represented as
1 T = Gu + L× S , (1.13) J K where S1 is the current defined by integration over S1 and L is an integer multiplicity 1- rectifiableJ K current, which keeps track of the possible concentration of ∇uk on 1-dimensional sets (see also Figure 3). The L1 -lower semicontinuous envelope of (1.11) is then (see, e.g., [34, Theorem 7.4]) Z (c) 1 |∇u| dx + |D u|(Ω) + I(u; Ω) , u ∈ BV (Ω; S ) , (1.14) Ω where I(u; Ω) is a surface energy depending on the jump part of Du and obtained by optimizing in the 1-current L. The latter is reminiscent of the energy due to concentration effects of ∇uk and it satisfies Z − + 1 1 dS (u , u ) dH ≤ I(u; Ω) .
Ju
(See [37] to interpret the energy (1.14) in terms of the optimal lifting problem in BV .) x Note that in the example (1.12) I(u; Ω) 6= 0, even if |x| has no jumps, and, more precisely, I(u; Ω) = 2π , obtained with L supported on a radius of Ω = B . If L = 0, then the inequality above is actually an equality.
1 It can be defined, e.g., by uk(x) = exp(ιϕk(x)), where ϕk(x) is the smooth function defined as follows. Define the function ψek in polar coordinates (ρ, ϑ) by ψek(ρ, ϑ) = ϑ, if ϑ ∈ [0, 2π − 1/k], and ψek(ρ, ϑ) = ∞ (2πk − 1)(2π − ϑ), if ϑ ∈ [2π − 1/k, 2π]. Let ϕek be a suitable regularization of ψek . Let ζk ∈ Cc (B2/k) be a cut-off function such that 0 ≤ ζk ≤ 1, ζk ≡ 1 on B1/k , |∇ζk| ≤ 2k and define ϕk = ϕek(1 − ζk). Then
Z Z Z 1 Z 2π−1/k Z 1 Z 2π |∂ϑψek| |∇uk| dx ∼ dx = dϑ dρ + (2πk − 1) dϑ dρ → 2π + 2π , B\B2/k B\B2/k ρ 2/k 0 2/k 2π−1/k Z Z |∇uk| dx ≤ |∇ϕek|(1 + ζk) + |ϕek||∇ζk| dx → 0 . B2/k\B1/k B2/k\B1/k FROM THE N -CLOCK MODEL TO THE XY MODEL 7
u+
− S1 J K
u− L
u 0
Ω
1 Figure 3. Depiction of a cartesian current of the form T = Gu − L× S . In the picture, Ω is the unit disc centered at the origin. The function u has no jumpsJ andK presents a vortex- like singularity, turning once counterclockwise around the origin. In particular, the graph 1 Gu has a hole, namely, ∂Gu = −δ0× S . The current T features a concentration part 1 J K −L× S . It is supported on a radius of the ball and is characterized by a vertical part (in J K 1 − + gray) that connects clockwise in S the (equal) traces u and u of u on the two sides of the radius. Note that the vertical part is not given by the geodesic connecting u− and u+ . The concentration part is needed to compensate the boundary of the graph Gu , so that ∂T = 0 . 1 1 Indeed, −∂L× S = δ0× S = −∂Gu inside Ω . In conclusion, the current T does not turn J K J K 1 around the origin. In this figure, for H -a.e. x in the support of L , the length `T (x) is 2π .
1 We now discuss the asymptotic behavior of Eε when ε θε ε| log ε|. (The εθε importance of the lower bound ε θε will be highlighted after the statement of Theorem 1.3 1 below.) In this regime, the limit of Eε shares strong similarities with (1.14). To a spin εθε 1 field uε :Ωε → Sε with equibounded energy Eε(uε) ≤ C we associate a current Gu , εθε ε 1 whose main feature is its relation with the discrete vorticity measure ∂Guε = −µuε × S , cf. J K Proposition 3.11. Thanks to the latter property and (1.9), the limit T of the currents Guε satisfies ∂T = 0 and, more precisely, T ∈ cart(Ω×S1). 1 To define the energy obtained as Γ-limit of Eε in this regime, we need to introduce εθε further quantities. It will feature an anisotropic surface term J (u; Ω), playing the same role of I(u; Ω) in (1.14). More precisely, for u ∈ BV (Ω; S1) we consider (see (4.15) for the precise definition) Z 1 1 J (u; Ω) = inf `T (x)|νT (x)|1 dH (x): T ∈ cart(Ω×S ) ,T as in (1.13) for some L .
JT (1.15) Here JT is the 1-dimensional jump-concentration set of T oriented by the normal νT , accounting for both the jump set of u and the support of the concentration part L in the decomposition (1.13). At each point x ∈ JT , the current T has a vertical part, given by a 1 curve in S which connects the traces of u on the two sides of JT ; `T (x) is its length. We are now in a position to state the theorem.
Theorem 1.2 (Regime ε θε ε| log ε| – no vortices). Assume that ε θε ε| log ε|. Then the following results hold:
i) (Compactness) Let uε :Ωε → Sε be such that 1 Eε(uε) ≤ C. εθε f Then µuε * 0. In addition, there exists a subsequence (not relabeled) and a function 1 1 2 u ∈ BV (Ω; S ) such that uε → u in L (Ω; R ). 8 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
1 ii)( Γ-liminf inequality) Assume that uε :Ωε → Sε and u ∈ BV (Ω; S ) are such that 1 2 uε → u in L (Ω; R ). Then Z (c) 1 |∇u|2,1 dx + |D u|2,1(Ω) + J (u; Ω) ≤ lim inf Eε(uε) . ε→0 εθε Ω 1 iii)( Γ-limsup inequality) Let u ∈ BV (Ω; S ). Then there exists a sequence uε :Ωε → Sε 1 2 such that uε → u in L (Ω; R ) and Z 1 (c) lim sup Eε(uε) ≤ |∇u|2,1 dx + |D u|2,1(Ω) + J (u; Ω) . ε→0 εθε Ω 1 In Theorem 1.2 it is proven that the bound Eε(uε) ≤ C forces the absence of vortices εθε in the limit allowing, instead, for BV -type concentration effects. This raises the natural 1 question of finding a bound on Eε(uε) which allows for configurations that also concen- εθε trate energy on finitely many vortices. The energy needed for the formation of a single vortex of type x−x0 can be deduced |x−x0| from (1.10) (cf. (4.78) for the rigorous bound). Indeed, for the discretization vε :Ωε → Sε of x−x0 we have that |x−x0| 1 ε Eε(vε) ∼ 2π| log ε| . (1.16) εθε θε Notice that, since θε ε| log ε|, the right-hand side diverges as ε → 0. Estimate (1.16) suggests that a bound of the type 1 ε Eε(uε) ≤ C| log ε| (1.17) εθε θε is compatible with the formation of finitely many vortices. Indeed, assuming (1.17), we f PN shall prove in Proposition 4.10 that µuε * µ (up to a subsequence), with µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z, and 2π|µ|(Ω) ≤ C . A finer analysis can be done by assuming the more precise bound 1 ε Eε(uε) − 2πM| log ε| ≤ C, (1.18) εθε θε 1 namely, assuming finiteness of the excess energy obtained by removing from Eε(uε) the εθε leading order term corresponding to a vortex configuration µ with |µ|(Ω) = M . Condition (1.18) suggests that a sequence uε may display in the limit simultaneously vortex-type and BV -type concentration effects. More precisely, (1.18) implies an energy f PN bounded as in (1.17) and, in particular, that µuε * µ, with µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z, and |µ|(Ω) ≤ M . In fact, the interesting situation to study is |µ|(Ω) = M . In that case, once the diverging energy 2πM| log ε| ε has been saturated through the formation θε 1 of the vortices µ, a finite energy Eε is still accessible to the system and it might lead εθε to BV -type concentration effects as in Theorem 1.2. As already discussed, to a spin field 1 uε :Ωε → Sε we associate a current Guε with ∂Guε = −µuε × S . This induces a nontrivial J K 1 constraint on the current T given by the limit of the Guε , namely ∂T = −µ× S . This condition couples the vortex-type and BV -type concentration effects displayed byJ theK spin field. 1 The next theorem shows that the limit of Eε features a coupling term J (µ, u; Ω) εθε generalizing J (u; Ω) in (1.15) given by Z 1 J (µ, u; Ω) := inf `T (x)|νT (x)|1 dH (x): T ∈ Adm(µ, u; Ω) ,
JT where Adm(µ, u; Ω), defined in (4.22), is a suitable class of currents T satisfying, in par- ticular, the constraint ∂T = −µ× S1 . Observe that J (µ, u; Ω) reduces to J (u; Ω) for µ = 0. J K
Theorem 1.3 (Regime ε θε ε| log ε| – M vortices). Assume that ε θε ε| log ε|. Then the following results hold: FROM THE N -CLOCK MODEL TO THE XY MODEL 9
i) (Compactness) Let M ∈ N and let uε :Ωε → Sε be such that 1 ε Eε(uε) − 2πM| log ε| ≤ C. εθε θε PN Then there exists a measure µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z, such that (up f to a subsequence) µuε * µ and |µ|(Ω) ≤ M . If, in addition, |µ|(Ω) = M , then 1 there exists a function u ∈ BV (Ω; S ) such that (up to a subsequence) uε → u in L1(Ω; R2). PN ii)( Γ-liminf inequality) Let uε :Ωε → Sε and let µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z f 1 with |µ|(Ω) = M . Assume that µuε * µ. Let u ∈ BV (Ω; S ) be such that uε → u in L1(Ω; R2). Then Z (c) 1 ε |∇u|2,1 dx + |D u|2,1(Ω) + J (µ, u; Ω) ≤ lim inf Eε(uε) − 2πM| log ε| . ε→0 εθε θε Ω PN iii)( Γ-limsup inequality) Let µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z with |µ|(Ω) = M and 1 f let u ∈ BV (Ω; S ). Then there exists a sequence uε :Ωε → Sε such that µuε * µ, 1 2 uε → u in L (Ω; R ), and Z 1 ε (c) lim sup Eε(uε) − 2πM| log ε| ≤ |∇u|2,1 dx + |D u|2,1(Ω) + J (µ, u; Ω) . ε→0 εθε θε Ω
To motivate the condition ε θε , it is convenient to recall the first order analysis of the XY model developed in [8, Section 4] (see also [16, 48, 10] for the Ginzburg-Landau model). f PM Let us suppose that µuε * µ, where µ = h=1 dhδxh , with dh = ±1. By (1.7)–(1.8), the least energy needed to form M vortices is 2πM| log ε| and additional information on the physical system are obtained by carrying out a finer analysis under the bound on the excess energy 1 XY (u ) − 2πM| log ε| ≤ C. ε2 ε ε Specifically, in [8, Theorem 4.2] the following asymptotic expansion of the energy 1 XY (u ) ∼ 2πM| log ε| + (µ) + Mγ (1.19) ε2 ε ε W is rigorously proven in a variational sense. In the formula, W is a Coulomb-type interaction potential referred to as renormalized energy, while γ is the core energy carried by each vortex, cf. (7.2) and (7.4) for the precise definitions. Analogously, the result stated in Theorem 1.3 leads to the formal expansion 1 θ Z E (u ) ∼ 2πM| log ε| + ε |∇u| dx + |D(c)u| (Ω) + J (µ, u; Ω) . (1.20) ε2 ε ε ε 2,1 2,1 Ω
If ε θε , the BV -type energy in (1.20) has a divergent pre-factor, in contrast to the finite term W(µ) + Mγ in (1.19). If, instead, θε ε, we expect the following formal asymptotic behavior to hold true: 1 E (u ) ∼ 2πM| log ε| + (µ) + Mγ . ε2 ε ε W This intuition is confirmed in the next theorem.
Theorem 1.4 (Regime θε ε). Assume that θε ε. Then the following results hold: i) (Compactness) Let M ∈ N and let uε :Ωε → Sε be such that 1 E (u ) − 2πM| log ε| ≤ C. ε2 ε ε PN Then there exists a measure µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z, with |µ|(Ω) ≤ M f such that (up to a subsequence) µuε * µ. Moreover, if |µ|(Ω) = M , then |dh| = 1. 10 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
PM ii)( Γ-liminf inequality) Let µ = h=1 dhδxh , xh ∈ Ω, |dh| = 1, and let uε :Ωε → Sε f be such that µuε * µ. Then 1 W(µ) + Mγ ≤ lim inf Eε(uε) − 2πM| log ε| . ε→0 ε2 PM iii)( Γ-limsup inequality) Let µ = h=1 dhδxh , xh ∈ Ω, |dh| = 1. Then there exists f uε :Ωε → Sε such that µuε * µ and 1 lim sup 2 Eε(uε) − 2πM| log ε| ≤ W(µ) + Mγ . ε→0 ε We remark that, in contrast to Theorem 1.3, only vortices with multiplicity ±1 are allowed in Theorem 1.4. The proof of the Γ-limsup inequality in Theorem 1.4 can be obtained by a careful adap- tation of some of the techniques adopted in the literature about the first order analysis of the Ginzburg-Landau [16, 48, 10] and of the classical XY model [8, 23]. In contrast, the construction of the recovery sequence in the proof of Theorems 1.1, 1.2, 1.3 turns out to be a significantly more delicate task. We highlight here some of its main technical difficulties, referring to Section 4 for a detailed explanation. Given u ∈ BV (Ω; S1), we define its recov- ery sequence following a gradual approximation procedure, which involves a series of steps of increasing complexity. At each of these steps, the map u is modified without essentially changing the energy. The first main issue is to regularize the map u. It is known that, in general, a map belonging to BV (Ω; S1) cannot be approximated in energy by means of S1 -valued smooth functions. Nonetheless, the result in [15] (see also [6]) guarantees the density of S1 -valued maps that are smooth outside finitely many point-singularities (whose number and position are however not provided). Such a regularization suffices to proceed with the construction in Theorem 1.1; in contrast, a more precise approximation is necessary for the purposes of Theorems 1.2, 1.3. Indeed, in that case the limit energy depends on the vorticity measure µ, which needs to be precisely related to singularities emerging in the regularization of u. This is achieved using the Approximation Theorem for cartesian currents, cf. Lemma 4.17. The next main issue is to construct a recovery sequence uε for such a regularization of u. Close to each singularity, uε is defined by discretizing on the ε-lattice a proper translation x of |x| and then by projecting its values on Sε . The energy carried by this discrete spin field close to a singularity is sensitive to the regime of θε : if ε| log ε| θε 1 as in Theorem 1.1, its energy is negligible; if, instead, ε θε ε| log ε| as in Theorems 1.2, 1.3, then the energy diverges according to (1.16). Far from the singularities, the problem reduces to the construction of a recovery sequence for a smooth S1 -valued map. This can be further simplified to the case of a piecewise constant S1 -valued map by introducing a mesoscopic scale into the problem, see Lemma 4.13. For such maps, the construction is a refinement of that in the example presented above and depicted in Figure 1. The most delicate step is to merge the different parts of the recovery sequence close to and far from the singularities. This is achieved in the proof of Proposition 4.22 (Step 2) by a careful interpolation on dyadic layers of mesoscopic squares, whose size is chosen to x be smaller for layers closer to the singularity. At each layer generation, |x| is sampled at a different mesoscopic length-scale. The latter is optimized in order to provide the correct control on the energy in progressing from each layer to the next one. We collect here the results proved in this paper. Theorems 1.1, 1.2, 1.3, and 1.4 cover the regimes of θε from θε finite to θε ε. For each regime we have scaled the Nε -clock energy Eε to highlight concentration effects of the spin field. Accordingly, we have shown that the energy may concentrate on sets of various dimensions. On the one hand, if θε is finite, the physical system is an Ising-type system with finitely many phases, for which it is well known that the energy concentrates on 1-dimensional domain walls separating phases of constant spin, cf. [19, 45]. On the other hand, we prove in Theorem 1.4 that, for θε ε, the Nε -clock model approximates the XY model (in FROM THE N -CLOCK MODEL TO THE XY MODEL 11 the sense of the asymptotic expansion (1.19)). In particular, the spin field develops 0- dimensional singularities, the vortices. For ε θε 1, we observed two intermediate 1 phenomena. When θε 1, every value of S is admissible for the limit spin field, and the limit energy detects both the formation of 1-dimensional domain walls and diffuse transitions. As shown in Theorem 1.1, this completes the description of the limit energy in the case ε| log ε| θε 1. If, instead, ε θε ε| log ε|, the previous phenomenon is coupled in a nontrivial way with the possible formation of vortices. The precise behavior of the energy is described in Theorems 1.2–1.3. All our results are summarized in Table 1. We remark that the coexistence of singularities of two different dimensions has been already included in the framework of XY -type models in [13] and of Ginzburg-Landau-type energies in [35]. While in the aforementioned papers the presence of singularities of various dimensions is enforced by suitably modifying the multi-well potential or by adding a surface term to the usual energy, here they appear as a result of the dependence on ε of the codomain.
Regime Energy bound Limit of µuε Energy behavior
1 θε finite ε Eε ≤ C not relevant interfaces
1 ε| log ε| θε Eε ≤ C not relevant BV εθε
1 f ε θε ε| log ε| Eε ≤ C µu * 0 BV + concentration εθε ε
vortices 1 ε f ε θε ε| log ε| Eε − 2πM| log ε| ≤ C µu * µ + εθε θε ε BV + concentration
1 f θε ε ε2 Eε − 2πM| log ε| ≤ C µuε * µ XY
Table 1. In this table we summarize our results. By “Ising” we mean that the energy concen- trates on 1 -dimensional domain walls that separate the different phases, while “ BV ” denotes a BV -type total variation. The expression “ BV +concentration” indicates the presence in a BV -type energy of a surface term of the form J (u; Ω) or J (µ, u; Ω) which accounts for con- centration effects on 1-dimensional surfaces. By “vortices” we mean that a diverging energy is carried by the system for the creation of vortex-like singularities in the limit. Finally, “ XY ” expresses the fact that the energy is a good approximation (at first order) of the classical XY model.
In this paper we also study the critical regime θε ∼ ε| log ε| in Section 5. In this case, 1 ε 1 0 a bound Eε(uε) − 2M| log ε| ≤ C is equivalent to both 2 Eε(uε) ≤ C and εθε θε ε | log ε| 1 0 0 Eε(uε) ≤ C for some other constant C . The former bound is compatible with the εθε f formation of vortices and, in fact, µuε * µ up to a subsequence; the latter bound implies compactness for uε in BV and for the currents Guε . In this setting, the excess energy 1 ε Eε(uε) − 2M| log ε| cannot diverge to −∞, thus we do not need to require additional εθε θε assumptions on the mass of µ. Indeed, we prove the following result.
Theorem 1.5 (Regime θε ∼ ε| log ε|). Assume that θε = ε| log ε|. Then the following results hold:
i) (Compactness) Let uε :Ωε → Sε be such that 1 E (u ) ≤ C. ε2| log ε| ε ε
PN Then there exists a measure µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z, such that (up to a f 1 subsequence) µuε * µ and there exists a function u ∈ BV (Ω; S ) such that (up to 1 2 a subsequence) uε → u in L (Ω; R ). 12 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
PN ii)( Γ-liminf inequality) Let uε :Ωε → Sε , let µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z, and 1 f 1 2 let u ∈ BV (Ω; S ). Assume that µuε * µ and uε → u in L (Ω; R ). Then Z (c) 1 |∇u|2,1 dx + |D u|2,1(Ω) + J (µ, u; Ω) + 2π|µ|(Ω) ≤ lim inf Eε(uε) . ε→0 ε2| log ε| Ω PN 1 iii)( Γ-limsup inequality) Let µ = h=1 dhδxh , xh ∈ Ω, dh ∈ Z and let u ∈ BV (Ω; S ). f 1 2 Then there exists a sequence uε :Ωε → Sε such that µuε * µ, uε → u in L (Ω; R ), and Z 1 (c) lim sup 2 Eε(uε) ≤ |∇u|2,1 dx + |D u|2,1(Ω) + J (µ, u; Ω) + 2π|µ|(Ω) . ε→0 ε | log ε| Ω
In the critical regime θε ∼ ε one can adapt Step 1 in the proof of Proposition 4.22 to 1 show that the variational limit of ε2| log ε| Eε is given by 2π|µ|(Ω), where µ is the vorticity 1 measure. When the energy is scaled as ε2 Eε one may expect from (1.19)–(1.20) a nontrivial interplay between the BV -type energy and the renormalized and core energies. This appears to deserve a separate analysis and is not treated in this paper.
2. Notation and preliminaries We denote the imaginary unit by ι. We shall often identify R2 with the complex plane C. 2 Given a vector a = (a1, a2) ∈ R , its 1-norm is |a|1 = |a1| + |a2|. We define the (2, 1)-norm 2×2 of a matrix A = (aij) ∈ R as the sum of the Euclidean norms of its columns, i.e., 2 2 1/2 2 2 1/2 |A|2,1 := a11 + a21 + a12 + a22 . 1 1 1 If u, v ∈ S , their geodesic distance on S is denoted by dS (u, v). It is given by the angle 1 in [0, π] between the vectors u and v , i.e., dS (u, v) = arccos(u · v). Observe that 1 1 1 2 |u − v| = sin 2 dS (u, v) , (2.1) and, in particular, we have the equivalence of distances in the sense that π |u − v| ≤ d 1 (u, v) ≤ |u − v| . (2.2) S 2 αε Given two sequences αε and βε , we write αε βε if limε→0 = 0. We will use the βε notation deg(u)(x0) to denote the topological degree of a continuous map u ∈ C(Bρ(x0) \ 1 {x0}; S ), i.e., the topological degree of its restriction u|∂Br (x0) , independent of r < ρ. We denote by Iλ(x) the half-open squares given by 2 Iλ(x) = x + [0, λ) . (2.3) 2.1. BV-functions. In this section we recall basic facts about functions of bounded varia- tion. For more details we refer to the monograph [12]. Let O ⊂ Rd be an open set. A function u ∈ L1(O; Rn) is a function of bounded variation if its distributional derivative Du is given by a finite matrix-valued Radon measure on O. In that case, we write u ∈ BV (O; Rn). n n The space BVloc(O; R ) is defined as usual. The space BV (O; R ) becomes a Banach space when endowed with the norm kukBV (O) = kukL1(O) +|Du|(O), where |Du| denotes the total variation measure of Du. The total variation with respect to the anisotropic norm | · |2,1 is n denoted by |Du|2,1 . When O is a bounded Lipschitz domain, then BV (O; R ) is compactly 1 n ∗ n embedded in L (O; R ). We say that a sequence un converges weakly in BV (O; R ) to 1 n ∗ u if un → u in L (O; R ) and Dun * Du in the sense of measures. We state some fine properties of BV -functions. To this end, we need some definitions. A function u ∈ L1(O; Rn) is said to have an approximate limit at x ∈ O whenever there exists z ∈ Rn such that 1 Z lim |u(y) − z| dy = 0 . ρ→0 d ρ Bρ(x) FROM THE N -CLOCK MODEL TO THE XY MODEL 13
Next we introduce so-called approximate jump points. Given x ∈ O and ν ∈ Sd−1 we set ± Bρ (x, ν) = {y ∈ Bρ(x): ±(y − x) · ν > 0} .
We say that x ∈ O is an approximate jump point of u if there exist a 6= b ∈ Rn and ν ∈ Sd−1 such that 1 Z 1 Z lim |u(y) − a| dy = lim |u(y) − b| dy = 0 . ρ→0 d + ρ→0 d − ρ Bρ (x,ν) ρ Bρ (x,ν) The triplet (a, b, ν) is determined uniquely up to the change to (b, a, −ν). We denote it by + − (u (x), u (x), νu(x)) and let Ju be the set of approximate jump points of u. The triplet + − (u , u , νu) can be chosen as a Borel function on the Borel set Ju . Denoting by ∇u the approximate gradient of u, we can decompose the measure Du as Z Z + − d−1 (c) Du(B) = ∇u dx + (u (x) − u (x)) ⊗ νu(x) dH + D u(B) , B Ju∩B
(c) (j) + − d−1 where D u is the so-called Cantor part and D u = (u − u ) ⊗ νuH Ju is the so-called jump part. We will need the slicing properties of BV -functions. Let d ≥ 2. Given a unit vector ξ ∈ Sd−1 , we denote by Πξ the hyperplane orthogonal to ξ . For every set E ⊂ Rd and z ∈ Πξ , ξ the section of E corresponding to z is the set Ez := {t ∈ R : z + tξ ∈ E}. Accordingly, n ξ ξ n ξ for any function u: E → R , the function uz : Ez → R is defined by uz(t) := u(z + tξ). We recall a characterization of BV functions by slicing [12, Remark 3.104]. Let us fix an open set O ⊂ Rd and u ∈ L1(O; Rn). Then u ∈ BV (O; Rn) if and only if for every d−1 ξ ξ n d−1 ξ ξ ∈ S we have uz ∈ BV (Oz; R ) for H -a.e. z ∈ Π and Z ξ ξ d−1 |Duz|(Oz) dH (z) < ∞ . Πξ Moreover it is possible to reconstruct the distributional gradient Du from the gradients of ξ d−1 ξ ξ the slices Duz through the formula Du ξ = H Π ⊗ Duz , i.e., Z ξ ξ d−1 Du ξ(B) = Duz(Bz ) dH (z) , Πξ for every Borel set B ⊂ Rd . More precisely, the same decomposition holds true for each part of the decomposition of Du, namely Z Z ξ ξ d−1 ∇u ξ dx = ∇uz(Bz ) dH (z) , B Πξ Z (c) (c) ξ ξ d−1 D u ξ(B) = D uz(Bz ) dH (z) , Πξ Z (j) (j) ξ ξ d−1 D u ξ(B) = D uz(Bz ) dH (z) , Πξ
d ξ d−1 ξ ξ ± for every Borel set B ⊂ . Moreover, J ξ = (Ju) for H -a.e. z ∈ Π and (u ) (t) = R uz z z ± ξ ∓ ξ ξ (u )z(t) (= (u )z(t), respectively) for every t ∈ (Ju)z if ξ·νu(z+tξ) > 0 (if ξ·νu(z+tξ) < 0, respectively). If O is a bounded open set, we define the space BV (O; S1) as the space of those functions u ∈ BV (O; R2) such that u(x) ∈ S1 almost everywhere. We remark that a slightly different approach to define the set BV (O; S1) is taken in [11] using S1 as a metric space. The results of this paper (see, for instance, the limit functionals in Theorem 1.1) fit more in that abstract framework, but we decided to take the above definition for simplicity. 14 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
2.2. Results for the classical XY model. We recall here some results about the classi- 1 cal XY model, namely when the spin field uε :Ωε → S is not constrained to take values in a discrete set. Following [7], in order to define the discrete vorticity of the spin variable, it is convenient to introduce the projection Q: R → 2πZ defined by Q(t) := argmin{|t − s| : s ∈ 2πZ} , (2.4) with the convention that, if the argmin is not unique, then we choose the one with minimal modulus. Then for every t ∈ R we define (see Figure 4) Ψ(t) := t − Q(t) ∈ [−π, π] . (2.5)
π
−2π −π 0 π 2π
−π
Figure 4. Graph of the function Ψ for t ∈ (−2π, 2π) . Observe that Ψ is an odd function.
2 1 Let u: εZ → S and let ϕ:Ωε → [0, 2π) be the phase of u defined by the relation u = exp(ιϕ). The discrete vorticity of u is defined for every εi ∈ εZ2 by 1 h du(εi) := Ψ ϕ(εi + εe1) − ϕ(εi) + Ψ ϕ(εi + εe1 + εe2) − ϕ(εi + εe1) 2π (2.6) i + Ψ ϕ(εi + εe2) − ϕ(εi + εe1 + εe2) + Ψ ϕ(εi) − ϕ(εi + εe2) .
As already noted in [7], the discrete vorticity du only takes values in {−1, 0, 1}, i.e., only singular vortices can be present in the discrete setting. We introduce the discrete measure representing all vortices of the discrete spin field defined as X µu := du(εi)δεi+(ε,ε) . (2.7) 2 εi∈εZ
Remark 2.1. In [7, 8] the vorticity measure ˚µu is supported in the centers of the squares completely contained in Ω, i.e., X ˚µu = du(εi)δεi+1/2(ε,ε) . 2 εi∈εZ εi+[0,ε]2⊂Ω In this paper we prefer definition (2.7) since it fits well with our definition of discrete currents in Section 3.5 on the whole set Ω. However, as we will borrow some results from [7, 8], we have to ensure that these definitions are asymptotically equivalent with respect to the flat convergence defined below. Definition 2.2 (Flat convergence). Let O ⊂ R2 be an open set. A sequence of measures µj ∈ Mb(O) converges flat to µ ∈ Mb(O) if Z Z 0,1 ψ dµj → ψ dµ for every ψ ∈ Cc (O) . (2.8) O O f In that case, we denote the convergence by µj * µ. Observe that the flat convergence is weaker than the weak* convergence. The two notions are equivalent when the measures µj have equibounded total variations. The two vorticity measures µu and ˚µu are then close in the following sense. 2 1 1 Lemma 2.3. Assume that uε : εZ → S is a sequence such that ε2 Eε(uε) ≤ C| log ε|. f Then µuε Ω − ˚µuε * 0. FROM THE N -CLOCK MODEL TO THE XY MODEL 15
0,1 Proof. Note that, for any ψ ∈ Cc (Ω) with kψkW 1,∞ ≤ 1, for ε > 0 small enough (depend- ing on the support of ψ ) we have X ε ε 1 |hµ Ω − ˚µ , ψi| ≤ |d (εi)|√ ≤ √ |˚µ |(Ω) ≤ Cε E (u) ≤ Cε| log ε| , uε uε u u ε2 ε 2 2 2 εi∈εZ εi+[0,ε]2⊂Ω where in the last but one inequality we used [7, Remark 3.4]. This proves the claim. We recall the following compactness and lower bound for the XY model.
1 Proposition 2.4. Let uε :Ωε → S and assume that there is a constant C > 0 such that 1 PN ε2| log ε| XYε(uε) ≤ C . Then there exists a measure µ ∈ Mb(Ω) of the form µ = h=1 dhδxh f with dh ∈ Z and xh ∈ Ω, and a subsequence (not relabeled) such that µuε Ω * µ. Moreover 1 2π|µ|(Ω) ≤ lim inf XYε(uε) . ε→0 ε2| log ε| Proof. In [8, Theorem 3.1-(i)] it is proven that (up to a subsequence) the discrete vorticity measures ˚µuε converge flat to a measure of the claimed form satisfying also the lower bound. The claim follows by applying Lemma 2.3.
Remark 2.5. Observe that in the regime θε ε| log ε| the bound 1 Eε(uε) ≤ C εθε f and Proposition 2.4 imply that µuε Ω * 0.
3. Currents For the theory of currents and cartesian currents we refer to the books [31, 32]. We recall here some notation, definitions, and basic facts about currents. We additionally prove some technical lemmata.
3.1. Definitions and basic facts. Given an open set O ⊂ Rd , we denote by Dk(O) the space of k -forms ω : O 7→ ΛkRd that are C∞ with compact support in O.A k -current k T ∈ Dk(O) is an element of the dual of D (O). The duality between a k -current and a k -form ω will be denoted by T (ω). The boundary of a k -current T is the (k−1)-current k−1 ∂T ∈ Dk−1(O) defined by ∂T (ω) := T (dω) for every ω ∈ D (O) (or ∂T := 0 if k = 0). As for distributions, the support of a current T is the smallest relatively closed set K in O 0 such that T (ω) = 0 if ω is supported outside K . Given a smooth map f : O → O0 ⊂ RN such that f is proper2, f #ω ∈ Dk(O) denotes the pull-back of a k -form ω ∈ Dk(O0) through 0 f . The push-forward of a k -current T ∈ Dk(O) is the k -current f#T ∈ Dk(O ) defined by # k f#T (ω) := T (f ω). Given a k -form ω ∈ D (O), we can write it via its components
X α ∞ ω = ωα dx , ωα ∈ Cc (O) , |α|=k where the expression |α| = k denotes all multi-indices α = (α1, . . . , αk) with 1 ≤ αi ≤ d, and dxα = dxα1∧...∧dxαk . The norm of ω(x) is denoted by |ω(x)| and it is the Euclidean norm of the vector with components (ωα(x))|α|=k . The total variation of a k -current T ∈ Dk(O) is defined by |T |(O) := sup{T (ω): ω ∈ Dk(O), |ω(x)| ≤ 1} .
If T ∈ Dk(O) with |T |(Ω) < ∞, then we can define the measure |T | ∈ Mb(O) |T |(ψ) := sup{T (ω): ω ∈ Dk(O), |ω(x)| ≤ ψ(x)}
2that means, f −1(K) is compact in O for all compact sets K ⊂ O0 . 16 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
for every ψ ∈ C0(O), ψ ≥ 0. As a consequence of Riesz’s Representation Theorem (see [31, ~ d ~ 2.2.3, Theorem 1]) there exists a |T |-measurable function T : O 7→ ΛkR with |T (x)| = 1 for |T |-a.e. x ∈ O such that Z T (ω) = hω(x), T~(x)i d|T |(x) (3.1)
O for every ω ∈ Dk(O). We note that if T has finite total variation, then it can be extended to a linear functional acting on all forms with bounded, Borel-measurable coefficients via the dominated convergence theorem. In particular, in this case the push-forward f#T can be defined also for f ∈ C1(O,O0) with bounded derivatives, cf. the discussion in [31, p. 132]). A set M ⊂ O is a countably Hk -rectifiable set if it can be covered, up to an Hk -negligible subset, by countably many k -manifolds of class C1 . As such, it admits at Hk -a.e. x ∈ M a tangent space Tan(M, x) in a measure theoretic sense. A current T ∈ Dk(O) is an integer multiplicity (i.m.) rectifiable current if it is representable as Z T (ω) = hω(x), ξ(x)iθ(x) dHk(x) , for ω ∈ Dk(O) , (3.2)
M where M ⊂ O is a Hk -measurable and countably Hk -rectifiable set, θ : M → Z is locally k d k H M-summable, and ξ : M → ΛkR is a H -measurable map such that ξ(x) spans Tan(M, x) and |ξ(x)| = 1 for Hk -a.e. x ∈ M. We use the short-hand notation T = τ(M, θ, ξ). One can always replace M by the set M ∩ θ−1({0}), so that we may always assume that θ 6= 0. Then the triple (M, θ, ξ) is uniquely determined up to Hk -negligible modifications. Moreover, one can show, according to the Riesz’s representation in (3.1), that T~ = ξ and the total variation3 is given by |T | = |θ|Hk M. If Tj are i.m. rectifiable currents and Tj *T in Dk(O) with supj(|Tj|(V ) + |∂Tj|(V )) < ∞ for every V ⊂⊂ O, then by the Closure Theorem [31, 2.2.4, Theorem 1] T is an i.m. rectifiable current, too. By M we denote the current defined by integration over M. J K 3.2. Currents in product spaces. Let us introduce some notation for currents defined on the product space Rd1 ×Rd2 . We will denote by (x, y) the points in this space. The standard d1 basis of the first space R is {e1, . . . , ed1 }, while {e¯1,..., e¯d2 } is the standard basis of the d2 d1 d2 second space R . Given O1 ⊂ R ,O2 ⊂ R open sets, T1 ∈ Dk1 (O1), T2 ∈ Dk2 (O2) and k1+k2 a (k1 + k2)-form ω ∈ D (O1×O2) of the type X α β ω(x, y) = ωαβ(x, y) dx ∧dy ,
|α|=k1 |β|=k2 the product current T1 × T2 ∈ Dk1+k2 (O1×O2) is defined by X X β α T1×T2(ω) := T1 T2 ωαβ(x, y) dy dx ,
|α|=k1 |β|=k2 α β while T1×T2(φ dx ∧dy ) = 0 if |α| + |β| = k1 + k2 but |α|= 6 k1 , |β|= 6 k2 .
3.3. Graphs. Let O ⊂ Rd be an open set and u:Ω → R2 a Lipschitz map. Then we can consider the d-current associated to the graph of u given by Gu := (id×u)# O ∈ 2 2 J K D2(O×R ), where id×u: O → O×R is the map (id×u)(x) = (x, u(x)). Note that by definition we have Z Gu(ω) = hω(x, u(x)),M(∇u(x))i dx O for all ω ∈ Dd(O × R2), with the d-vector 1 2 1 2 M(∇u) = (e1 + ∂x1 u e¯1 + ∂x1 u e¯2) ∧ ... ∧ (ed + ∂xd u e¯1 + ∂xd u e¯2) . (3.3)
3For i.m. rectifiable currents, the total variation coincides with the so-called mass. Hence, we will not distinguish between these two concepts. FROM THE N -CLOCK MODEL TO THE XY MODEL 17
The above formula can be extended to the class A1(O; R2) defined by 1 2 1 2 1 A (O; R ) := {u ∈ L (O; R ): u approx. diff. a.e. and all minors of ∇u are in L (O)} . 1,2 2 1 2 2 Remark 3.1. We recall that ∂Gu|Ω×R = 0 when u ∈ W (O; R ) ⊂ A (O; R ), see [31, 3.2.1, Proposition 3]. This property however fails for general functions u ∈ A1(O; R2). In Lemma 3.4 we need to interpret the graphs of W 1,1(O; S1) as currents. This can be done because of the following observation.
Lemma 3.2. Let O ⊂ Rd be an open, bounded set. Then W 1,1(O; S1) ⊂ A1(O; R2). Proof. It is well-known that Sobolev functions are approximately differentiable a.e. More- over, all 1-minors of ∇u are in L1(O). We argue that all 2-minors vanish at a.e. point. To this end, denote by P : R2 \{0} → R2 the smooth mapping P (x) = x/|x|. Since for u ∈ W 1,1(O; S1) we have u = P ◦ u almost everywhere, for a.e. x ∈ O the chain rule for approximate differentials yields ∇u(x) = ∇P (u(x))∇u(x). Since ∇P (u(x)) has at most rank 1, also ∇u(x) has at most rank 1 and therefore all 2-minors have to vanish as claimed. Later on we use the orientation of the graph of a smooth function u: O ⊂ R2 → S1 (cf. 2 [31, 2.2.4]). For such maps we have |Gu| = H M, where M = (id×u)(Ω), and p 2 1 + |∇u(x)| G~ u(x, y) = e1 ∧ e2 1 2 + ∂x2 u (x)e1 ∧ e¯1 + ∂x2 u (x)e1 ∧ e¯2 (3.4) 1 2 − ∂x1 u (x)e2 ∧ e¯1 − ∂x1 u (x)e2 ∧ e¯2 , for every (x, y) ∈ M.
3.4. Cartesian currents. Let O ⊂ Rd be a bounded, open set. We recall that the class of cartesian currents in O×R2 is defined by 2 2 2 cart(O×R ) := {T ∈ Dd(O×R ): T is i.m. rectifiable, ∂T |O×R = 0, O π#T = O ,T | dx ≥ 0 , |T | < ∞ , kT k1 < ∞} , J K O 2 where π : O×R → O denotes the projection on the first component, T | dx ≥ 0 means ∞ 2 that T (φ(x, y) dx) ≥ 0 for every φ ∈ Cc (O×R ) with φ ≥ 0, and ∞ 2 kT k1 = sup{T (φ(x, y)|y| dx): φ ∈ Cc (O×R ) , |φ| ≤ 1} . Note that, if for some function u Z Z T (φ(x, y) dx) = φ(x, u(x)) dx then kT k1 = |u| dx . (3.5) O O The class of cartesian currents in O×S1 is 1 2 1 cart(O×S ) := {T ∈ cart(O×R ) : supp(T ) ⊂ O×S } , (cf. [32, 6.2.2] for this definition). We recall the following approximation theorem which explains that cartesian currents in O×S1 are precisely those currents that arise as limits of graphs of S1 -valued smooth maps. The proof, based on a regularization argument on the lifting of T , can be found in [33, Theorem 7].4
Theorem 3.3 (Approximation Theorem). Let T ∈ cart(O×S1). Then there exists a se- ∞ 1 quence of smooth maps uh ∈ C (O; S ) such that 2 Guh *T in Dd(O×R ) and 2 2 |Guh |(O×R ) → |T |(O×R ) .
4Notice that some results in [33] require O to have smooth boundary. This is not the case for this theorem, which is based on a local construction. 18 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
Using the above approximation result, we now prove an extension result for cartesian currents, which we could not find in the literature.
Lemma 3.4 (Extension of cartesian currents). Let O ⊂ Rd be a bounded, open set with Lipschitz boundary and let T ∈ cart(O×S1). Then there exist an open set Oe ⊃⊃ O and a 1 2 2 current T ∈ cart(Oe×S ) such that Te|O×R = T and |Te|(∂O×R ) = 0. ∞ 1 Proof. Applying Theorem 3.3 we find a sequence uk ∈ C (O; S ) such that Guk *T in 2 2 2 2 Dd(O×R ) and |Guk |(O×R ) → |T |(O×R ). In particular, the sequence |Guk |(O×R ) is bounded, which implies that Z sup |∇uk|dx < C . (3.6) k O
Next we extend the functions uk . To this end, note that there exists t > 0 and a bi- Lipschitz map Γ: (∂O×(−t, t)) → Γ(∂O×(−t, t)) such that Γ(x, 0) = x for all x ∈ ∂O, Γ(∂O×(−t, t)) is an open neighborhood of ∂O and
2 Γ(∂O×(−t, 0)) ⊂ O, Γ(∂O × (0, t)) ⊂ R \ O. (3.7) This result is a consequence of [44, Theorem 7.4 & Corollary 7.5]; details can be found for instance in [43, Theorem 2.3]. The extension of uk is then achieved via reflection. More 0 0 0 precisely, for a sufficiently small t > 0 we define it on O with O = O + Bt0 (0) by
( −1 uk(Γ(P (Γ (x)))) if x∈ / O, uek(x) = (3.8) uk(x) otherwise,
1,1 0 1 where P (x, τ) = (x, −τ). Since Γ is bi-Lipschitz, we have that uek ∈ W (O ; S ) and by a change of variables we can bound the L1 -norm of its gradient via Z Z Z Z −1 |∇uek| dx ≤ |∇uk| dx + CΓ |(∇uk) ◦ Γ ◦ P ◦ Γ | dx ≤ CΓ |∇uk| dx , (3.9) O0 O O0\O O where the constant CΓ depends only on the bi-Lipschitz properties of Γ and the dimension. Due to Lemma 3.2 we obtain that u ∈ A1(O0; 2). In particular, the current G ∈ ek R uek 0 2 Dd(O × R ) is well-defined in the sense of Z G (ω) = hω(x, u (x)),M(∇u (x))i dx , uek ek ek O0 with M(∇u ) given by (3.3). We next prove that G ∈ cart(O0 × 1). First note that ek uek S whenever ω ∈ Dd(O0 ×R2) is a form with supp(ω) ⊂⊂ O0 ×R2 \(O0 ×S1), then the definition yields G (ω) = 0. By standard arguments (apply for instance [31, 3.2.1 Proposition 5] to uek a constant sequence) it then suffices to prove that ∂G | 0 2 = 0. We will argue locally. uek O ×R 0 For each x ∈ O we choose a rotation Qx , radii rx > 0, and heights hx > 0 such that the d−1 cylinders Cx := x + Qx (−rx, rx) × (−hx, hx) satisfy
(i) Cx ⊂⊂ O if x ∈ O; 0 0 (ii) Cx ⊂⊂ O \ O if x ∈ O \ O; 0 (iii) Cx ⊂⊂ O if x ∈ ∂O and 0 d 0 d−1 0 Cx ∩ O = Cx ∩ x + Qx{(x , xd) ∈ R : x ∈ (−rx, rx) : −hx < xd < ψ(x )} for some ψ ∈ Lip(Rd−1). ∞ 1 For x ∈ O we have ∂G | 2 = ∂G | 2 = 0 since u ∈ C (C ; ). Next consider uek Cx×R uk Cx×R k x S 0 0 the second case, namely x ∈ O \ O. Since Cx ⊂⊂ O \ O, the properties in (3.7) imply −1 that Γ ◦ P ◦ Γ (Cx) ⊂⊂ O. In particular, by the smoothness of uk on O we have that 1,∞ u ∈ W (C ), so that by Remark 3.1 again ∂G | 2 = 0. Finally, we consider ek x uek Cx×R x ∈ ∂O. Since Cx ∩ Ω is (up to a rigid motion) the subgraph of a Lipschitz function, it is in particular simply connected. By classical lifting theory, we find a sequence of scalar functions FROM THE N -CLOCK MODEL TO THE XY MODEL 19
∞ ϕk ∈ C (Cx ∩ O) such that uk(x) = exp(ιϕk(x)). In particular, using the chain rule we 1,1 see that ϕk ∈ W (Cx ∩ O). Now fix 0 < δ < rx small enough such that Bδx (x) ⊂ Cx and −1 0 (Γ ◦ P ◦ Γ )(Bδx (x) ∩ (O \ O)) ⊂ Cx ∩ O, which can be realized due to (3.7). We then extend the lifting ϕk to a function ϕek ∈ 1,1 W (Bδx (x)) via the same reflection construction as in (3.8), which is well-defined due to the above inclusions. Observe that this definition guarantees that uek(y) = exp(ιϕek(y)) for almost every y ∈ Bδx (x). Expressed in terms of currents this means that
G | 2 = χ G | 2 , uek Bδx (x)×R # ϕek Bδx (x)×R where G ∈ D (B (x)× ) is the current associated to the graph of ϕ and χ: d× → ϕek d δx R ek R R Rd×S1 is the covering map defined by χ(x, ϑ) := (x, cos(ϑ), sin(ϑ)). In particular, by [33, 1 Theorem 2, p. 97 & Proposition 1 (i), p. 100] we have G | 2 ∈ cart(B (x) × ), uek Bδx (x)×R δx S so that by the definition of cartesian currents we have ∂G | 2 = 0. uek Bδx (x)×R 0 0 Thus we have shown that for every x ∈ O there exists a ball Bδx (x) ⊂ O such that ∂G | 2 = 0. Using a partition of unity to localize the support of any form ω ∈ uek Bδx (x)×R d−1 0 2 D (O × ) with respect to the x-variable, we conclude that ∂G | 0 2 = 0 and R uek O ×R therefore G ∈ cart(O0 × 1). As seen in the proof of Lemma 3.2, all 2-minors of Du uek S vanish, so that the bounds (3.6) and (3.9) yield
Z Z 1 |G |(O0× 2) = |M(∇u )| dx ≤ 1 + |∇u |2 2 dx ≤ C. uek R ek ek O0 O0 Hence, up to a subsequence, we can assume that G * T in D (O0× 2), see [31, 2.2.4 uek e d R 0 2 Theorem 2]. From [31, 4.2.2. Theorem 1] it follows that Te ∈ cart(O ×R ). Since uek = uk 2 2 on O, we find that Te|O×R = T . It remains to show that |Te|(∂O × R ) = 0. To this end, note that for 0 < η < η0 < 1, by the bi-Lipschitz continuity of Γ and (3.7) we have that −1 out in (Γ ◦ P ◦ Γ )(Oη ) ⊂ Oη0 , out in where the sets Oη and Oη0 are defined as out 0 in 0 Oη := {x ∈ O \ O : dist(x, ∂O) < η},Oη0 = {x ∈ O : dist(x, ∂O) < η } . Hence, similar to (3.9) we obtain that
Z 1 2 2 2 in 2 |G |((∂O + B (0))× ) ≤ C 1 + |∇u | dx = C |G |(O 0 × ) . (3.10) uek η R Γ k Γ uk η R Oin η0 2 2 0 Since |Guk |(O×R ) → |T |(O×R ) and |T | is a finite measure, there exists a sequence η → 0 in 2 in 2 such that |Guk |(Oη0 ×R ) → |T |(Oη0 ×R ). Applying the lower semicontinuity of the mass with respect to weak convergence of currents in (3.10), we infer that 2 in 2 |Te|((∂O + Bη(0))×R ) ≤ CΓ|T |(Oη0 ×R ) . 0 2 Sending first η → 0 and then η → 0 we conclude that |Te|(∂O×R ) = 0 as claimed. We will also use the structure theorem for cartesian currents in O×S1 that has been proven in [33, Section 3, Theorems 1, 5, 6].5 However, to simplify notation, from now on we focus on dimension two. Recall that Ω ⊂ R2 is a bounded, open set with Lipschitz boundary. To state the theorem, we recall the following decomposition for a current T ∈ cart(Ω×S1). Letting M the countably H2 -rectifiable set where T is concentrated, we denote by M(a) the set of points (x, y) ∈ M at which the tangent plane Tan(M, (x, y)) does not contain vertical vectors (namely, the Jacobian of the projection πΩ restricted to Tan(M, (x, y)) has Ω (jc) (a) 1 dπ#|T | maximal rank), by M := (M\M )∩(JT ×S ), where JT := {x ∈ Ω: dH1 (x) > 0}, and by M(c) := M\ (M(a) ∪ M(jc)). Then we can split the current via T = T (a) + T (c) + T (jc) ,
5As for the Approximation Theorem, no boundary regularity is required for this result. 20 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF where T (a) := T M(a) , T (c) := T M(c) , T (jc) := T M(jc) are mutually singular measures, and we denote by the restriction of the Radon measure T . Hereafter we use 1 2 2 1 the notation xb = x and xb = dx . Theorem 3.5 (Structure Theorem for cart(Ω×S1)). Let T ∈ cart(Ω×S1). Then there 1 exists a unique map uT ∈ BV (Ω; S ) and an (not unique) i.m. rectifiable 1-current LT = ~ (jc) (j) 1 τ(L, k, LT ) ∈ D1(Ω) such that T = T + LT × S and J Z K (a) T (φ(x, y) dx) = T (φ(x, y) dx) = φ(x, uT (x)) dx , (3.11) Ω Z T (a)(φ(x, y) dxl∧dym) = (−1)2−l φ(x, u (x))∂(a)um(x) dx , (3.12) b T xl T Ω Z T (c)(φ(x, y) dxl∧dym) = (−1)2−l φ(x, u (x)) d∂(c)um(x) , (3.13) b eT xl T Ω Z Z T (j)(φ(x, y) dxl∧dym) = (−1)2−l φ(x, y) dym νl (x) dH1(x) (3.14) b uT γ JuT x ∞ 2 1 − for every φ ∈ Cc (Ω×R ), γx being the (oriented) geodesic arc in S that connects uT (x) + to uT (x) and ueT being the precise representative of uT . Remark 3.6. In [33, Theorem 6] the structure of T (j) is formulated in a slightly different way, − − + + using the counter-clockwise arc γϕ−,ϕ+ between (cos(ϕ ), sin(ϕ )) and (cos(ϕ ), sin(ϕ )) and replacing JuT by Jϕ , where ϕ ∈ BV (Ω) is a local lifting of T . More precisely, the 2 2 1 notion of lifting is understood in the sense that T = χ#Gϕ , where χ: R ×R → R ×S is the covering map (x, ϑ) 7→ (x, cos(ϑ), sin(ϑ)) and Gϕ ∈ cart(Ω×R) is the cartesian current given by the boundary of the subgraph of ϕ (hence the push-forward via χ is well-defined as Gϕ has finite mass, see Section 3.1). To explain how to deduce (3.14), we recall the local + 0 construction in [33]: for every x ∈ Jϕ one chooses p (x) ≥ 0 and k (x) ∈ N ∪ {0} such that ϕ+(x) = p+(x) + 2πk0(x), 0 ≤ p+(x) − ϕ−(x) < 2π , where we recall that in the scalar case the traces (and the normal to the jump set) are − + 0 arranged to satisfy ϕ < ϕ on Jϕ . Then, locally, the 1-current LT in [33, Theorem 6] 0 0 0 ~ 0 0 0 is given by LT = τ(L , k (x), LT ), where L ⊂ Jϕ denotes the set of points with k (x) ≥ 1 ~ 0 0 ~ 0 2 1 and LT is the orientation of L defined via LT = νϕe1 −νϕe2 . To obtain the representation via geodesics, we let ( (p+(x), k0(x)) if p+(x) − ϕ−(x) < π , (q+(x), k(x)) = (p+(x) − 2π, k0(x) + 1) if p+(x) − ϕ−(x) > π , The case p+(x) − ϕ−(x) = π , i.e, antipodal points, needs special care. In this case we + ± ± define q (x) and k(x) according to the following rule: let ϕe (x) := ϕ (x) mod 2π ∈ [0, 2π). Then
( + 0 + − + (p (x), k (x)) if Ψ(ϕe (x) − ϕe (x)) = π , (q (x), k(x)) = + 0 + − (p (x) − 2π, k (x) + 1) if Ψ(ϕe (x) − ϕe (x)) = −π , with the function Ψ defined in (2.5). Replacing (p+(x), k0(x)) by (q+(x), k(x)), the modified structure of T (j) can be proven following exactly the lines of [33, p.107-108], noting that by + − the chain rule in BV [12, Theorem 3.96] we have Jϕ = JuT ∪ {x ∈ Jϕ : q (x) = ϕ (x)}. In particular, ~ ~ 2 1 LT = τ(L, k, LT ) , L = {x ∈ Jϕ : k(x) ≥ 1} , LT = νϕe1 − νϕe2 (3.15) still depend on the local lifting ϕ, but in (3.14) the curves γϕ−,ϕ+ are replaced by the − − − + more intrinsic geodesic arcs γx connecting uT (x) = (cos(ϕ (x)), sin(ϕ (x))) to uT (x) = + + (cos(ϕ (x)), sin(ϕ (x))) (these formulas are consistent with the choice νuT (x) = νϕ(x)). FROM THE N -CLOCK MODEL TO THE XY MODEL 21
− + In particular, exchanging uT (x) and uT (x) will change the orientation of the arc (also in 6 (j) the case of antipodal points) and of the normal νuT (x), so that the formula for T is invariant, hence well-defined without the use of local liftings. The structure theorem highlights a peculiar property of cartesian currents in Ω×S1 : the current T cannot be written only in terms of its associated BV function uT because of 1 the presence of the concentration term LT × S . This feature can be seen in the typical 1 J K x example of a current T ∈ cart(Ω×S ) whose associated BV function is uT (x) = |x| a.e. 1 in B1(0). Such a current must be of the form T = G x + LT × S , where LT is an i.m. |x| J K rectifiable 1-current with ∂LT |B1(0) = δ0 . E.g., LT can be concentrated on any curve L which connects the point 0 to the boundary. The current LT is therefore necessary to compensate the possible presence of holes in the graph associated to uT . It is convenient to recast the jump-concentration part of T ∈ cart(Ω×S1) in the following 1 way. Let LT = τ(L, k, L~ T ) as in Theorem 3.5. We introduce for H -a.e. x ∈ JT the normal
νT (x) to the 1-rectifiable set JT = JuT ∪ L as ( ν (x) if x ∈ J , ν (x) = uT uT (3.16) T ~ 2 ~ 1 (−LT (x), LT (x)) if x ∈ L \ JuT , ~ 2 ~ 1 1 where we choose νuT (x) = (−LT (x), LT (x)) if x ∈ L∩JuT . For H -a.e. x ∈ JT we consider T − + the curve γx given by: the (oriented) geodesic arc γx which connects uT (x) to uT (x) if 1 x ∈ JuT \L (in the sense of Remark 3.6 in case of antipodal points); the whole S turning k(x) times if x ∈ L \ JuT , k(x) being the integer multiplicity of LT ; the sum (in the sense 7 1 of currents) of the oriented geodesic arc γx and of S with multiplicity k(x) if x ∈ JuT ∩L. Then Z Z (jc) l m 2−l n mo l 1 T (φ(x, y) dxb∧dy ) = (−1) φ(x, y) dy νT (x) dH (x) . (3.17) T JT γx T m The integration over γx with respect to the form dy in the formula above is intended T 1 with the correct multiplicity of the curve γx defined for H -a.e. x ∈ JT by the integer number ±1 , if x ∈ Ju \L , y ∈ supp(γx) , T k(x) , if x ∈ L \ J , y ∈ 1, m(x, y) := uT S (3.18) k(x) ± 1 , if x ∈ L ∩ J , y ∈ supp(γ ) , uT x T k(x) , if x ∈ L ∩ JuT , y ∈ supp(γx ) \ supp(γx) , where ± = +/− if the geodesic arc γx is oriented counterclockwise/clockwise, respectively. More precisely, Z Z m m m 1 φ(y) dy = (−1) φ(y)yb m(x, y) dH (y) . (3.19) T T γx supp(γx )
Remark 3.7. Note that we constructed m(x, y) based on the orientation (3.16) of νT . As discussed in Remark 3.6, changing the orientation of νuT changes the orientation of the geodesic γx , while a change of the orientation of L~ T switches the sign of k(x). Hence changing the orientation of νT (x) changes m(x, y) into −m(x, y). If we choose locally νT = νϕ as in Remark 3.6, our construction above yields m(x, y) ≥ 0. In the proposition below, we derive an explicit formula for the vector T~ of a cartesian current. It seems that this result is well-known to experts, but since we could not find a precise reference, we include a proof for the reader’s convenience.
6 1 More precisely, assume that u1, u2, ν ∈ S and assume that the geodesic arc from u1 to u2 is coun- + − terclockwise. If (uT (x), uT (x), νuT (x)) = (u2, u1, ν), then γx is oriented counterclockwise. If, instead, + − (uT (x), uT (x), νuT (x)) = (u1, u2, −ν) (equivalent to the first choice, according to the definition of jump point), then γx is oriented clockwise. 7 T 1 In this case, a more elementary way of defining γx is the following: let γx : [0, 1] → S be the geodesic arc, and let ϕx : [0, 1] → R be a continuous function (unique up to translations of an integer multiple of 2π ) T such that γx(t) = exp(ιϕx(t)). Then γx (t) = exp ι(1 − t)ϕx(0) + ιt(ϕx(1) + 2πk(x)) . 22 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
1 Proposition 3.8. Let T ∈ cart(Ω×S ), let uT be the BV function associated to T . Then |T (a)| = H2 M(a) , |T (c)| = H2 M(c) , |T (jc)| = |m|H2 M(jc) , and
p 2 1 + |∇uT (x)| T~(x, y) = e1 ∧ e2 (a) 1 (a) 2 + ∂x2 uT (x)e1 ∧ e¯1 + ∂x2 uT (x)e1 ∧ e¯2 (3.20) (a) 1 (a) 2 − ∂x1 uT (x)e2 ∧ e¯1 − ∂x1 uT (x)e2 ∧ e¯2 , for H2 -a.e. (x, y) ∈ M(a) , d∂(c)u1 d∂(c)u2 T~(x, y) = x2 T (x)e ∧ e¯ + x2 T (x)e ∧ e¯ d|D(c)u | 1 1 d|D(c)u | 1 2 T T (3.21) d∂(c)u1 d∂(c)u2 − x1 T (x)e ∧ e¯ − x1 T (x)e ∧ e¯ , (c) 2 1 (c) 2 2 d|D uT | d|D uT | for H2 -a.e. (x, y) ∈ M(c) , and 2 2 2 1 sign(m(x, y))T~(x, y) = − ν (x)y e1 ∧ e¯1 + ν (x)y e1 ∧ e¯2 T T (3.22) 1 2 1 1 + νT (x)y e2 ∧ e¯1 − νT (x)y e2 ∧ e¯2 , for H2 -a.e. (x, y) ∈ M(jc) , where m(x, y) is the integer defined in (3.18). Proof. Assume Ω simply connected (if not, the following arguments can be repeated locally). Let us consider the covering map χ:Ω×R → Ω×S1 defined by χ(x, ϑ) := (x, cos(ϑ), sin(ϑ)). By [33, Corollary 1, p. 105] there exists a lifting of T , i.e., there is a function ϕ ∈ BV (Ω; R) such that T = χ#Gϕ , where Gϕ ∈ cart(Ω×R) is the cartesian current given by the boundary of the subgraph of ϕ. The fine structure of such currents is well known, compare [29, Theorem 4.5.9], [31, 4.1.5 & 4.2.4]. We recall here that, if we consider the subgraph SGϕ := {(x, y) ∈ Ω×R : y < ϕ(x)}, then SGϕ is a set of finite perimeter; Gϕ is the current Gϕ = ∂ SGϕ . The interior normal to SGϕ is given by J K d(Dϕ, −L2) n(x, ϕ(x)) = (x) , for x ∈ Ω \ Jϕ , d|(Dϕ, −L2)| (3.23) − + n(x, ϑ) = (νϕ(x), 0) , for x ∈ Jϕ , ϑ ∈ [ϕ (x), ϕ (x)] , where νϕ is the normal to the jump set Jϕ . Moreover, the current Gϕ can be represented − as Gϕ = G~ ϕ|Gϕ| where |Gϕ| is concentrated on the reduced boundary ∂ SGu , |Gϕ| = 2 − ~ 3 ~ H ∂ SGu , and Gϕ is the 2-vector in R such that −Gϕ(x, ϑ) ∧ n(x, ϑ) = e1 ∧ e2 ∧ e3 , i.e., 3 2 1 G~ ϕ = −n e1 ∧ e2 + n e1 ∧ e3 − n e2 ∧ e3 . Finally, letting (a) 3 Σ := {(x, ϕe(x)) : x ∈ Ω \ Jϕ, n (x, u(x)) 6= 0} , (c) 3 Σ := {(x, ϕe(x)) : x ∈ Ω \ Jϕ, n (x, u(x)) = 0} , (j) − + 3 Σ := {(x, ϑ): x ∈ Jϕ, ϑ ∈ [ϕ (x), ϕ (x)], n (x, ϑ) = 0} ,
− (a) (c) (j) (a) (a) (c) (c) we have ∂ SGϕ = Σ ∪ Σ ∪ Σ and, denoting Gϕ = Gϕ Σ , Gϕ = Gϕ Σ , (j) (j) Gϕ = Gϕ Σ , by [33, formulas (2) and (16)] we have on the one hand that uT = (cos(ϕ), sin(ϕ)) a.e. and (a) (a) (c) (c) (j) (jc) χ#Gϕ = T , χ#Gϕ = T , χ#Gϕ = T . (3.24) − 4 On the other hand, observe that the Jacobian of dχ: Tan(∂ SGϕ, x) 7→ R equals 1 (in- deed dχ maps any pair of orthonormal vectors of R3 to a pair of orthonormal vectors in R4 ). Hence by the area formula, for σ ∈ {a, c, j} we obtain Z Z (σ) # ~ 2 X ~ 2 χ#Gϕ (ω) = hχ ω, Gϕi dH (x, ϑ) = hω(x, y), dχ(x, ϑ)Gϕ(x, ϑ)i dH (x, y) . (x,ϑ)∈χ−1(x,y) Σ(σ) χ(Σ(σ)) (3.25) FROM THE N -CLOCK MODEL TO THE XY MODEL 23
Next, note that for σ ∈ {a, c} the map χ:Σ(σ) → χ(Σ(σ)) is one-to-one and for any (a) (c) (x, ϕe(x)) ∈ Σ ∪ Σ we have ~ 3 dχ(x, ϕe(x))Gϕ(x, ϕe(x)) = − n (x, ϕe(x))e1 ∧ e2 2 2 − n (x, ϕe(x)) sin(ϕe(x))e1 ∧ e¯1 + n (x, ϕe(x)) cos(ϕe(x))e1 ∧ e¯2 1 1 + n (x, ϕe(x)) sin(ϕe(x))e2 ∧ e¯1 − n (x, ϕe(x)) cos(ϕe(x)e2 ∧ e¯2 . (3.26)
~ 2 Since |n| = 1 we see that | dχ(x, ϕe(x))Gϕ(x, ϕe(x))| = 1, too. Moreover, for H -a.e. (x, y) ∈ (σ) ~ χ(Σ ) the vector dχ(x, ϕe(x))Gϕ(x, ϕe(x)) orients the tangent space at (x, y). Hence (3.24) and the uniqueness of the representation of i.m. rectifiable currents (cf. Section 3.1) implies χ(Σ(σ)) = M(σ) up to a H2 -negligible set, |T (σ)| = H2 M(σ) , and ~ ~ T (χ(x, ϕe(x))) = dχ(x, ue(x))Gϕ(x, ϕe(x)) 2 (σ) for H -almost every (x, y) = χ(x, ϕe(x)) ∈ Σ . By the chain rule in BV [12, Theorem 3.96] we deduce that 2 2 −uT (c) −ueT (c) ∇uT = 1 ⊗ ∇ϕ , D uT = 1 ⊗ D ϕ . uT ueT Combined with the formula for n given by (3.23), the formulas (3.20) and (3.21) then follow from (3.26) by a straightforward calculation. In order to treat the case σ = j , note that due to (3.23) we have for any (x, y) = χ(x, ϑ) ∈ χ(Σ(j))
~ 2 2 2 1 dχ(x, ϑ)Gϕ(x, ϑ) = − νϕ(x)y e1 ∧ e¯1 + νϕ(x)y e1 ∧ e¯2 1 2 1 1 + νϕ(x)y e2 ∧ e¯1 − νϕ(x)y e2 ∧ e¯2 =: ξ(x, y) .
Again |ξ(x, y)| = 1 and ξ(x, y) orients the tangent space at H2 -a.e. (x, y) ∈ χ(Σ(j)). Thus (3.25) and the uniqueness of the representation of i.m. rectifiable currents imply (up to H2 -negligible sets) that M(jc) = χ(Σ(j)), T~ = ξ on M(jc) , and |T (jc)| = N(x, y)H2 M(jc) , with N(x, y) = #{ϑ ∈ [ϕ−(x), ϕ+(x)] : (cos(ϑ), sin(ϑ)) = y} .
To conclude, we have to relate m(x, y) to N(x, y) and νT (x) to νϕ(x). First note that the proof of the structure theorem (sketched in Remark 3.6) yields JuT ∪L = Jϕ and, combined T − + T with the definition of the curves γx (cf. (3.17)), implies that χ[ϕ (x), ϕ (x)] = supp(γx ) for x ∈ Jϕ . Hence
(jc) (j) 2 T M = χ(Σ ) = {(x, y) ∈ Ω×R : x ∈ JuT ∪ L , y ∈ supp(γx )} . (3.27)
Moreover, provided we orient JuT the same way as Jϕ and L according to (3.15), equa- tion (3.16) also yields νT = νϕ and m(x, y) = N(x, y) (a detailed proof of the latter requires to distinguish different cases, which we omit here).8 Inserting this equality in (3.27) con- cludes the proof of (3.22).
Finally, we recall the following result, proven in [33, Section 4].
1 1 Proposition 3.9. If u ∈ BV (Ω; S ), then there exists a T ∈ cart(Ω×S ) such that uT = u a.e. in Ω.
8 As noted in Remark 3.7, the choice νT (x) = νϕ(x) always yields m(x, y) ≥ 0. The factor sign(m(x, y)) in (3.22) makes the formula invariant under the change of νT (x). 24 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
3.5. Currents associated to discrete spin fields. We introduce the piecewise constant interpolations of spin fields. For every set S , we put 2 2 2 PCε(S) := {u: R → S : u(x) = u(εi) if x ∈ εi + [0, ε) for some i ∈ εZ } 1 Given u:Ωε → S , we can always identify it with its piecewise constant interpolation belong- 1 2 ing to PCε(S ), arbitrarily extended to R . Note that the piecewise constant interpolation of u coincides with u on the bottom-left corners of the squares of the lattice εZ2 . 1 2 Given u ∈ PCε(S ), we will associate to u the current Gu ∈ D2(Ω×R ) defined by Z 1 2 Gu(φ(x, y) dx ∧dx ) := φ(x, u(x)) dx , (3.28) Ω Z Z l m 2−l m l 1 Gu(φ(x, y) dxb∧dy ) := (−1) φ(x, y) dy νu(x) dH (x) , (3.29)
Ju γx 1 2 Gu(φ(x, y) dy ∧dy ) := 0 , (3.30)
∞ 2 for every φ ∈ Cc (Ω×R ), where Ju is the jump set of u, νu(x) is the normal to Ju at x, 1 − + and γx ⊂ S is the (oriented) geodesic arc which connects the two traces u (x) and u (x). + − 1 If u (x) and u (x) are opposite vectors, the choice of the geodesic arc γx ⊂ S is done consistently with the choice made in (2.4) for the values Ψ(π) and Ψ(−π) as follows: let ± ± + − ϕ (x) ∈ [0, 2π) be the phase of u (x); if Ψ(ϕ (x) − ϕ (x)) = π , then γx is the arc that − + + − connects u (x) to u (x) counterclockwise; if Ψ(ϕ (x) − ϕ (x)) = −π , then γx is the arc − + that connects u (x) to u (x) clockwise. Note that the choice of the arc γx is independent of the orientation of the normal νu(x). 1 We define for H -a.e. x ∈ Ju the integer number m(x) = ±1, where ± = +/− if the geodesic arc γx is oriented counterclockwise/clockwise, respectively. Then Z Z m m m 1 φ(y) dy = (−1) m(x) φ(y)yb dH (y) . (3.31)
γx supp(γx)
u+ εi u− εi + εe1
Figure 5. The current Gu has vertical parts concentrated on the jump set Ju , where a transition from u− to u+ occurs.
1 2 Proposition 3.10. Let u ∈ PCε(S ) and let Gu ∈ D2(Ω×R ) be the current defined in (3.28)–(3.30). Then Gu is an i.m. rectifiable current and, according to the representation 2 formula (3.1), Gu = G~ u|Gu|, where |Gu| = H M, (a) (j) M = M ∪ M = {(x, u(x)) : x ∈ Ω \ Ju} ∪ {(x, y): x ∈ Ju, y ∈ γx} , and
G~ u(x, y) = e1 ∧ e2 (3.32) for H2 -a.e. (x, y) ∈ M(a) and
2 2 2 1 G~ u(x, y) = sign(m(x)) − ν (x)y e1 ∧ e¯1 + ν (x)y e1 ∧ e¯2 u u (3.33) 1 2 1 1 + νu(x)y e2 ∧ e¯1 − νu(x)y e2 ∧ e¯2 for H2 -a.e. (x, y) ∈ M(j) . FROM THE N -CLOCK MODEL TO THE XY MODEL 25
Proof. First note that the set M is countably H2 -rectifiable. Since u is piecewise constant, for horizontal forms we have Z Z 2 (a) Gu(φ(x, y) dx) = φ(x, u(x)) dx = φ(x, y) dH M (x, y) .
2 Ω Ω×R By (3.31) we deduce that for l, m = 1, 2 Z Z l m 2−l m l 1 Gu(φ(x, y) dxb∧dy ) = (−1) φ(x, y) dy νu(x) dH (x)
Ju γx Z Z 2−l+m m 1 l 1 = (−1) φ(x, y)yb dH (y) νu(x)m(x) dH (x)
Ju supp(γx) Z 2−l+m m l 2 (j) = (−1) φ(x, y)yb νu(x)m(x) dH M (x, y) . 2 Ω×R 2 Then for every ω ∈ D2(Ω×R ) we have Z 2 Gu(ω) = hω, G~ ui dH M
2 Ω×R for G~ u defined as in (3.32)–(3.33) and moreover G~ u(x, y) is associated to the tangent space ~ 2 at (x, y) ∈ M. Since also |Gu(x, y)| = 1 for |Gu|-a.e. (x, y) ∈ Ω×R , we conclude the proof.
1 2 Proposition 3.11. Let u ∈ PCε(S ) and let Gu ∈ D2(Ω×R ) be the current defined in (3.28)–(3.30). Then 1 ∂Gu|Ω× 2 = −µu× S , R J K 2 1 2 where µu is the discrete vorticity measure defined in (2.7) for u|εZ : εZ → S . 1 2 Proof. Let us fix 0 < ρ < min{ε/4, dist(Ωε, ∂Ω)} and η ∈ D (Ω×R ). With a partition of unity we can split η into the sum of 1-forms depending on their supports. We discuss here all the possibilities for the supports. Case 1: supp(η) ⊂ (εi + (0, ε)2)×R2 for some i ∈ Z2 . Since u is constant in (εi + 2 2 (0, ε) )×R , we get automatically ∂Gu(η) = 0 by Remark 3.1. Case 2: Let H be the side of the square εi + [0, ε]2 connecting two vertices p, q ∈ εZ2 2 and let U be the ρ/2-neighborhood of H \ Bρ(p)∪Bρ(q) . Assume that supp(η) ⊂ U×R . We claim that
∂Gu(η) = 0 . (3.34) To prove this, we approximate the pure-jump function u by means of a sequence of Lipschitz ± functions uj . Let u be the traces of u on the two sides of H and let νH be the normal ± ± ± ± to H oriented as νu . We let ϕb ∈ [0, 2π) be the phases of u defined by u = exp(ιϕb ). − − + − + − We set ϕ := ϕb and ϕ := ϕb + Ψ(ϕb − ϕb ) ∈ (−π, 3π), where Ψ is the function given by (2.5). We then define ϕ− , if t ≤ − 1 2 − + − 1 1 1 ϕ(t) := ϕ + ϕ − ϕ (t + 2 ) if − 2 < t < 2 + 1 ϕ if t ≥ 2 , and ϕk(s) := ϕ(ks) for k large enough. Note that the curve t ∈ (−1/2, 1/2) 7→ exp(ιϕ(t)) 1 − + parametrizes the geodesic arc γ± ⊂ S which connects u to u , consistently with the choice done in formula (3.29). Then we put uk(x) := exp ιϕk(νH · (x − p)) for x ∈ U. 26 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
2 ∞ 2 We prove that Guk *Gu in D2(U×R ). Let us fix φ ∈ Cc (U×R ). Since uk → u in measure, we have that Z Z
Guk (φ(x, y) dx) = φ(x, uk(x)) dx → φ(x, u(x)) dx = Gu(φ(x, y) dx) . U U 0 0 Writing x ∈ U as x = x + sνH with x ∈ H , s ∈ R, for l = 1, 2 we further obtain that Z G (φ(x, y) dxl∧dy1) = (−1)2−l φ(x, u (x))∂ u1 (x) dx uk b k xl k U Z 3−l 0 l = (−1) φ(x, uk(x)) sin(ϕk(νH · (x − p)))ϕk(νH · x)νH dx U 1/2k Z Z 3−l 0 0 l 1 0 = (−1) φ x + sνH , exp ιϕk(s) sin ϕk(s) ϕk(s) ds νH dH (x ) H −1/2k 1/2 Z Z 3−l 0 t 0 l 1 0 = (−1) φ x + k νH , exp ιϕ(t) sin ϕ(t) ϕ (t) dt νH dH (x ) H −1/2 Z Z 2−l 0 t 1 l 1 0 = (−1) φ x + k νH , y dy νH dH (x )
H γ± Z Z 2−l 0 1 l 1 0 l 1 → (−1) φ x , y dy νH dH (x ) = Gu(φ(x, y) dxb∧dy ) ,
H γ± 1 − + where γ± ⊂ S is the geodesic arc connecting u to u . With analogous computations one proves G (φ(x, y) dxl∧dy2) → G (φ(x, y) dxl∧dy2). uk b u b Hence, due to Stokes’ Theorem we have that
0 = ∂Guk (η) = Guk (dη) → Gu(dη) = ∂Gu(η) , which proves (3.34). 2 2 Case 3: supp(η) ⊂ Bρ(p)×R , where p = εi + εe1 + εe2 for some i ∈ Z . In this case we will approximate the current Gu with graphs of a sequence of functions uk which are Lipschitz outside the point p. For notation simplicity we let ϕb1, ϕb2, ϕb3, ϕb4 ∈ [0, 2π) be the phases defined by the relations
u(εi + εe1 + εe2) =: u1 = exp(ιϕ1) , u(εi + εe2) =: u2 = exp(ιϕ2) , b b (3.35) u(εi) =: u3 = exp(ιϕb3) , u(εi + εe1) =: u4 = exp(ιϕb4) . We define the auxiliary angles
ϕσ(h+1) := ϕbσ(h) + Ψ(ϕbσ(h+1) − ϕbσ(h)) , (3.36) for h = 1, 2, 3, 4, where σ(h) ∈ {1, 2, 3, 4} is such that σ(h) ≡ h mod 4 (ϕσ(h+1) is the oriented angle in [−π, π] between the two vectors uσ(h) and uσ(h+1) ). We introduce the 2π -periodic function ϕk : R → R ϕ , if − π + h π < ϑ ≤ h π − 1 , bσ(h) 4 2 2 2k ϕ (ϑ) := π 1 π 1 π 1 k ϕbσ(h) + k ϕσ(h+1) − ϕbσ(h) (ϑ − h 2 + 2k ) , if h 2 − 2k < ϑ < h 2 + 2k π 1 π π ϕσ(h+1) , if h 2 + 2k ≤ ϑ < 4 + h 2 . π π π π for ϑ ∈ − 4 + h 2 , 4 + h 2 ), h ∈ Z. The function ϕk might have jumps at the points π π 4 + h 2 , h ∈ Z; note, however, that the amplitude of the jump is given by ϕbσ(h+1) − ϕσ(h+1) = ϕbσ(h+1) − ϕbσ(h) − Ψ(ϕbσ(h+1) − ϕbσ(h)) = Q(ϕbσ(h+1) − ϕbσ(h)) ∈ 2πZ , according to (2.4). FROM THE N -CLOCK MODEL TO THE XY MODEL 27
ϕ1
ϕ4 2π
π π π − 4 0 4 2 ϕ 1 1 k
Figure 6. Example for the definition of ϕk for h = 0 .
1 1 1 We now define a map vk : S → S . Given y ∈ S , let ϑ(y) ∈ [0, 2π) be the angle such that y = exp(ιϑ(y)) and set vk(y) := exp ιϕk(ϑ(y)) . The definition actually does not depend on the choice of the phase ϑ(y), due to the 2π - periodicity of ϕk . Thus we could also choose ϑ(y) ∈ [2πh, 2π(h + 1)) for any h ∈ Z. Note that vk is continuous: indeed the possible jumps of ϕk have amplitude in 2πZ, and thus are not seen by vk . In particular, we can compute the degree of the map vk via the formula Z Z 3 # X deg(v )2π = deg(v ) ω 1 = v ω 1 = ϕ − ϕ k k S k S σ(h+1) bσ(h) 1 1 h=0 S S 3 X = Ψ(ϕbσ(h+1) − ϕbσ(h)) = du(εi)2π , h=0 1 1 where ωS is the volume form on S and du(εi) is the discrete vorticity defined in (2.6). 1 We now define the map uk : Bρ(p) → S by x−p uk(x) := vk |x−p| . Note that, if (r, ϑ) are polar coordinates for the point x − p, then the polar coordinates of uk(x) are (1, ϕk(ϑ)).
u2 u1 u2 u1 1 ∼ k u3 u4 u3 u4
Figure 7. Example of the approximation uk (on the left) of the function u (on the right). The jump set of the function u is expanded and a transition between the jumps of u is 1 constructed using the geodesic arcs in S between the traces. If u has a nontrivial discrete vorticity as in the picture, then the graph Guk of the function uk has a hole in the center, x as it happens for the graph of the map x 7→ |x| . The hole is then preserved in the passage to limit to Gu , see formula (3.37).
By [31, 3.2.2, Example 2] we get that 1 1 1 2 2 ∂Guk |Bρ(p)× = − deg(vk)δp× S = −du(εi)δp× S = −µu× S |Bρ(p)× . (3.37) R J K J K J K R 2 Therefore, to conclude the proof it suffices to show the convergence Guk *Gu in D2(Ω×R ), so that 1 −µu× S (η) = ∂Guk (η) → ∂Gu(η) . ∞ J K2 To do so, let us fix φ ∈ Cc (Bρ(p)×R ). Since uk → u in measure, we have that Z Z
Guk (φ(x, y) dx) = φ(x, uk(x)) dx → φ(x, u(x)) dx = Gu(φ(x, y) dx) .
Bρ(p) Bρ(p) 28 MARCO CICALESE, GIANLUCA ORLANDO, AND MATTHIAS RUF
l m To compute the limit on forms of the type φ(x, y) dxb∧dy , observe that uk is not constant only in the 4 sectors of Bρ(p) given in polar coordinates by
h π 1 π 1 Ak := (r, ϑ): r ∈ (0, ρ), ϑ ∈ h 2 − 2k , h 2 + 2k , h ∈ {0, 1, 2, 3} , thus, for l, m = 1, 2, Z 2−l l m m (−1) G (φ(x, y) dx ∧dy ) = φ(x, u (x))∂ l u (x) dx uk b k x k
Bρ(p) 3 Z X m = φ(x, uk(x))∂xl uk (x) dx .
h=0 h Ak
h The integrals on the sets Ak can be computed in polar coordinates. We show the compu- tations for h = 0 and m = 1, the other cases being analogous. Changing variables in the integral on the interval (−1/2k, 1/2k) we obtain Z 1 φ(x, uk(x))∂x2 uk(x) dx
0 Ak ρ 1/2k Z Z 0 = − φ r exp ιϑ , exp ιϕk(ϑ) sin(ϕk(ϑ))ϕk(ϑ) cos(ϑ) dϑ dr 0 −1/2k ρ 1/2 Z Z t 0 t = − φ r exp ι k , γ41(t) sin(γ41(t))γ41(t) cos k dt dr 0 −1/2 ρ 1/2 Z Z Z Z 0 1 2 1 → − φ (r, 0), γ41(t) sin(γ41(t))γ41(t) dt dr = φ(x, y) dy ν dH (x) ,
0 −1/2 J41 γ41 and
ρ 1/2k Z Z Z 1 0 φ(x, uk(x))∂x1 uk(x) dx = φ r exp ιϑ , exp ιϕk(ϑ) sin(ϕk(ϑ))ϕk(ϑ) sin(ϑ) dϑ dr 0 0 Ak −1/2k ρ 1/2 Z Z t 0 t = φ r exp ι k , γ41(t) sin(γ41(t))γ41(t) sin k dt dr → 0 , 0 −1/2