---
abstract: 'We study the Ginzburg-Landau model of type-I superconductors in the regime of small external magnetic fields. We show that, in an appropriate asymptotic regime, flux patterns are described by a simplified branched transportation functional. We derive the simplified functional from the full Ginzburg-Landau model rigorously via $\Gamma$-convergence. The detailed analysis of the limiting procedure and the study of the limiting functional lead to a precise understanding of the multiple scales contained in the model.'
author:
- 'Sergio Conti[^1]'
- 'Michael Goldman[^2]'
- 'Felix Otto[^3]'
- 'Sylvia Serfaty[^4]'
bibliography:
- 'CGOS.bib'
title: '[A branched transport limit of the Ginzburg-Landau functional]{}'
---

Introduction
============

In 1911, K. Onnes discovered the phenomenon of superconductivity, manifested in the complete loss of resistivity of certain metals and alloys at very low temperature. W. Meissner discovered in 1933 that this was coupled with the expulsion of the magnetic field from the superconductor at the critical temperature. This is now called the Meissner effect. After some preliminary works of the brothers F. and H. London, V. Ginzburg and L. Landau proposed in 1950 a phenomenological model describing the state of a superconductor. In their model (see below), which belongs to Landau’s general theory of second-order phase transitions, the state of the material is represented by the order parameter $u: \Omega\to {\mathbb{C}}$, where $\Omega$ is the material sample. The density of superconducting electrons is then given by $\rho:= |u|^2$. A microscopic theory of superconductivity was first proposed by Bardeen-Cooper-Schrieffer (BCS) in 1957, and the Ginzburg-Landau model was derived from BCS by Gorkov in 1959 (see also [@frankseiringer] for a rigorous derivation).

One of the main achievements of the Ginzburg-Landau theory is the prediction and the understanding of the mixed (or intermediate) state below the critical temperature. This is a state in which, for moderate external magnetic fields, normal and superconducting regions coexist. The behavior of the material in the Ginzburg-Landau theory is characterized by two physical parameters. The first is the coherence length $\xi$ which measures the typical length on which $u$ varies, the second is the penetration length $\lambda$ which gives the typical length on which the magnetic field penetrates the superconducting regions. The Ginzburg-Landau parameter is then defined as $\kappa:= \frac{\lambda}{\xi}$. The Ginzburg-Landau functional is given by $$\label{GLorigin}
 \int_{\Omega} |\nabla_A u|^2 +\frac{\kappa^2}{2} (1-|u|^2)^2 dx +\int_{\R^3} |\nabla\times A- {B_{\mathrm{ext}}}|^2 dx$$ where $A: \R^3\to \R^3$ is the magnetic potential (so that $B:=\nabla\times A$ is the magnetic field), $ \nabla_A u:= \nabla u -i Au$ is the covariant derivative of $u$ and ${B_{\mathrm{ext}}}$ is the external magnetic field. In these units, the penetration length $\lambda$ is normalized to $1$. As first observed by A. Abrikosov this theory predicts two types of superconductors. On the one hand, when $\kappa < 1/\sqrt{2}$, there is a positive surface tension which leads to the formation of normal and superconducting regions corresponding to $\rho\simeq 0$ and $\rho\simeq 1$ respectively, separated by interfaces. These are the so-called type-I superconductors. On the other hand, when $\kappa>1/\sqrt{2}$, this surface tension is negative and one expects to see the magnetic field penetrating the domain through lines of vortices. These are the so-called type-II superconductors. In this paper we are interested in better understanding the former type but we refer the interested reader to [@tinkham; @SanSerf; @Serfrev] for more information about the latter type. In particular, in that regime, there has been an intensive work on understanding the formation of regular patterns of vortices known as Abrikosov lattices.

In type-I superconductors, it is observed experimentally [@ProzHo; @Proz; @Proal] that complex patterns appear at the surface of the sample. It is believed that these patterns are a manifestation of branching patterns inside the sample. Although the observed states are highly history-dependent, it is argued in [@ChokKoOt; @Proal] that the hysteresis is governed by low-energy configurations at vanishing external magnetic field. The scaling law of the ground-state energy was determined in [@ChokConKoOt; @ChokKoOt] for a simplified sharp interface version of the Ginzburg-Landau functional and in [@CoOtSer] for the full energy, these results indicate the presence of a regime with branched patterns at low fields.

This paper aims at a better understanding of these branched patterns by going beyond the scaling law. Starting from the full Ginzburg-Landau functional, we prove that in the regime of vanishing external magnetic field, low energy configurations are made of nearly one-dimensional superconducting threads branching towards the boundary of the sample. In a more mathematical language, we prove $\Gamma-$convergence [@braides; @dalmaso] of the Ginzburg-Landau functional to a kind of branched transportation functional in this regime. We focus on the simplest geometric setting by considering the sample $\Om$ to be a box $Q_{\Lz,T}:=(-\frac{\Lz}{2}, \frac{\Lz}{2})^2\times(-T,T)$ for some $T,\Lz>0$ and consider periodic lateral boundary conditions. The external magnetic field is taken to be perpendicular to the sample, that is ${B_{\mathrm{ext}}}:= {b_{\mathrm{ext}}}e_3$ for some ${b_{\mathrm{ext}}}>0$ and where $e_3$ is the third vector of the canonical basis of $\R^3$. After making an isotropic rescaling, subtracting the bulk part of the energy and dropping lower order terms (see and the discussion after it), minimizing can be seen as equivalent to minimizing

$$\begin{gathered}
\label{Etilde}
  E_T(u,A):=\frac{1}{L^2}
\int_{Q_{L,1}} |\nabla_{TA} u |^2  + \left(B_3 - \a (1-\rho)\right)^2 + |B'|^2 dx
+ \| B_3-\a \b\|_{H^{-1/2}(\{x_3=\pm 1\})}^2,\,\end{gathered}$$

where we have let $B:=(B',B_3):=\nabla \times A$, $$\kappa T:=\sqrt{2}\alpha,  \qquad {b_{\mathrm{ext}}}:=\frac{\beta \kappa}{\sqrt{2}} \qquad \textrm{ and } \qquad L:=\Lz/T.$$ If $u=\rho^{1/2} \exp(i\theta)$, since $|\nabla_{TA} u |^2 =|\nabla \rho^{1/2}|^2 +\rho|\nabla \theta-T A|^2$, in the limit $T\to +\infty$ we obtain, at least formally, that $A$ is a gradient field in the region where $\rho>0$ and therefore the Meissner condition $\rho B=0$ holds. Moreover, in the regime $\alpha\gg 1$, from we see that $B_3 \simeq \a (1-\rho)$ and $\rho$ takes almost only values in $\{0,1\}$. Hence $\Div B=0$ can be rewritten as $$  \partial_3 \chi + \frac1\alpha\Divp \chi B'=0,$$ where $\chi:= (1-\rho)$ and $\Divp$ denotes the divergence with respect to the first two variables. Therefore, from the Benamou-Brenier formulation of optimal transportation [@AGS; @villani] and since from the Meissner condition, $B'\simeq \frac{1}{\chi} B'$, the term $$\int_{Q_{L,1}} |B'|^2 dx\simeq \int_{Q_{L,1}} \frac{1}{\chi}|B'|^2 dx$$ in the energy can be seen as a transportation cost. We thus expect that inside the sample (this is, in $Q_{L,1}$), superconducting domains where $\rho\simeq1$ and $B\simeq 0$ alternate with normal ones where $\rho\simeq 0$ and $B_3\simeq \alpha$. Because of the last term $\| B_3-\a \b\|_{H^{-1/2}(\{x_3=\pm 1\})}^2$ in the energy , one expects $B\simeq \a \b e_3$ outside the sample. This implies that close to the boundary the normal domains have to refine. The interaction between the surface energy, the transportation cost and the penalization of an $H^{-1/2}$ norm leads to the formation of complex patterns (see Figure \[figscales\]).

It has been proven in [@CoOtSer] that in the regime $T\gg 1$, $\a \gg 1$ and $\b \ll 1$, $$\label{scallawintro}
\min E_T(u,A)\sim \min\{ \a^{4/3} \b^{2/3}, \a^{10/7} \b\}\,.$$ The scaling $\min E_T(u,A)\sim \a^{4/3} \b^{2/3}$ (relevant for $\a^{-2/7}\ll \b$) corresponds to uniform branching patterns whereas the scaling $\min E_T(u,A)\sim \a^{10/7} \b$ corresponds to non-uniform branching ones. We focus here for definiteness on the regime $\min E_T(u,A)\sim \a^{4/3} \b^{2/3}$, although we believe that our proof can be extended to the other one. Based on the construction giving the upper bounds in , we expect that in the first regime there are multiple scales appearing (see Figure \[figscales\]): $$\label{scalesepar}\stackrel{\displaystyle\textrm{penetration}}{\textrm{ length}}\,\stackrel{\displaystyle\ll}{\, } \, \stackrel{\displaystyle\textrm{ coherence}}{\textrm{ length}}\, \stackrel{\displaystyle\ll}{\,} 
\, \stackrel{\displaystyle\textrm{ diameter of the}}{\textrm{ threads in the bulk}} \, \,  \stackrel{\displaystyle\ll}{\,} \, \,  \stackrel{\displaystyle\textrm{distance between the}}{\textrm{ threads in the bulk}},$$ which amounts in our parameters to $$T^{-1} \ll \alpha^{-1}\ll \alpha^{-1/3}\beta^{1/3}\ll \alpha^{-1/3}\beta^{-1/6}.$$

\[figscales\]

In order to better describe the minimizers we focus on the extreme region of the phase diagram $T,T \a^{-1}, \beta^{-1}, \a \b^{7/2}\to +\infty$, with $L={\widetilde}L  \alpha^{-1/3}\beta^{-1/6}$ for some fixed ${\widetilde}L>0$. In this regime, we have in particular $\a^{-1}\ll \a^{-1/3}\b^{1/3}$ so that the separation of scales holds. We introduce an anisotropic rescaling (see Section \[Sec:GL\]) which leads to the functional

$$\begin{aligned}
\label{eqdefwe}
\wE(u,A)&:=\frac{1}{\widetilde{L}^2}\Big[\int_{Q_{{\widetilde}L,1}} \a^{-2/3}\b^{-1/3}\lt|\nabla_{\a^{1/3}\b^{-1/3}TA}'u\rt|^2+ 
  \a^{-4/3}\b^{-2/3}\lt|(\nabla_{\a^{1/3}\b^{-1/3}TA} u)_3\rt|^2\\
\nonumber
 +  & \lt.\a^{2/3}\b^{-2/3}\left({B_3} -(1-|{u}|^2)\right)^2
+  \b^{-1}|{B'}|^2 dx +\a^{1/3}\b^{7/6}\|\b^{-1}{B_3}-1\|_{H^{-1/2}(x_3=\pm1)}^2\rt].\end{aligned}$$

Our main result is a $\Gamma-$convergence result of the functional $\wE$ towards a functional defined on measures $\mu$ living on one-dimensional trees. These trees correspond to the normal regions in which $\rho\simeq 0$ and where the magnetic field $B$ penetrates the sample. Roughly speaking, if for a.e. $x_3\in (-1,1)$ the slice of $\mu=\mu_{x_3}\otimes dx_3$ has the form $\mu_{x_3}=\sum_i {\varphi}_i \delta_{X_i(x_3)}$ where the sum is at most countable, then we let (see Section \[Sec:limiting\] for a precise definition) $$\label{eqdefIintro}
I(\mu):=\frac1{\widetilde{L}^2}\int_{-1}^1 K_* \sum_i \sqrt{{\varphi}_i}+ {\varphi}_i |{\dot{X}}_i|^2 dx_3,$$ where $K_*=8\sqrt\pi/3$ and ${\dot{X}}_i$ denotes the derivative (with respect to $x_3$) of $X_i(x_3)$. The $X_i$’s represent the graphs of each branch of the tree (parametrized by height) and the ${\varphi}_i$’s represent the flux carried by the branch. We can now state our main result

\[maintheo\] Let $T_n, \a_n, \b_n^{-1}\to+\infty$ with $T_n \a_n^{-1}, \a_n \b_n^{7/2}\to +\infty$, ${\widetilde}L>0$, then:

1.  \[maintheolb\] For every sequence $(u_n, A_n)$ with $\sup_n \wE(u_n,A_n)<+\infty$, up to subsequence , $\beta_n^{-1}(1-|u_n|^2)$ weakly converges to a measure $\mu$ of the form $\mu=\mu_{x_3}\otimes dx_3$ with $\mu_{x_3}=\sum_i {\varphi}_i \delta_{X_i}$ for a.e. $x_3\in (-1,1)$, $\mu_{x_3} \weaklim dx'$ (where $dx'$ denotes the two dimensional Lebesgue measure on $Q_{{\widetilde}L}$) when $x_3\to \pm 1$ and such that $$\liminf_{n\to +\infty}  \wE(u_n,A_n)\ge I(\mu).$$

2.  \[maintheoub\] If in addition $L_n^{2} \a_n \b_n T_n\in 2\pi \Z$, where $L_n:= {\widetilde}L \a_n^{-1/3}\b_n^{-1/6}$, then for every measure $\mu$ such that $I(\mu)<+\infty$ and $\mu_{x_3} \weaklim dx'$ as $x_3\to \pm1$ , there exists $(u_n,A_n)$ such that $\b_n^{-1}(1-|u_n|^2)\weaklim \mu$ and $$\limsup_{n\to +\infty} \, \wE(u_n,A_n)\le I(\mu).$$

By scaling, it suffices to consider the case ${\widetilde}L=1$. The first assertion follows from Proposition \[gammaliminf\], the second one from Proposition \[gammalimsup-v2\].

Let us stress once again that our result could have been equivalently stated for the full Ginzburg-Landau energy instead of $\wE$ (see Section \[Sec:GL\]).\
Within our periodic setting, the quantization condition $L_n^2 \a_n \b_n T_n\in 2\pi \Z$ for the flux is a consequence of the fact that the phase circulation of the complex-valued function in the original problem is naturally quantized. It is necessary in order to make our construction but we believe that it is also a necessary condition for having sequences of bounded energy (see the discussion in Section \[Sec:GL\] and the construction in Section \[recoveryGL\]). We remark that scaling back to the original variables this condition is the physically natural one $L_0^2{b_{\mathrm{ext}}}\in 2\pi\Z$.\
Before going into the proof of Theorem \[maintheo\] we address the limiting functional $I(\mu)$, which has many similarities with irrigation (or branched transportation) models that have recently attracted a lot of attention (see [@BeCaMo] and more detailed comments in Section \[Sec:reliri\] or the recent papers [@BraWirth; @BraWir16] where the connection is also made to some urban planning models). In Section \[Sec:limiting\], we first prove that the variational problem for this limiting functional is well-posed (Proposition \[existmu\]) and show a scaling law for it (Proposition \[branchingmu\] and Proposition \[branchinglowerbound\]). In Proposition \[def:subsystem\], we define the notion of subsystems which allows us to remove part of the mass carried by the branching measure. This notion is at the basis of Lemma \[noloop\] and Proposition \[finitebranch\] which show that minimizers contain no loops and that far from the boundary, they are made of a finite number of branches. From the no-loop property, we easily deduce Proposition \[reg\] which is a regularity result for minimizers of $I$. The main result of Section \[Sec:limiting\] is Theorem \[theodens\] which proves the density of “regular” measures in the topology given by the energy $I(\mu)$. As in nearly every $\Gamma-$convergence result, such a property is crucial in order to implement the construction for the upper bound \[maintheoub\]. We now comment on the proof of Theorem \[maintheo\]. Let us first point out that if the Meissner condition $\rho B =0$ were to hold, and $A$ could be written as a gradient field in the set $\{\rho>0\}$, then $|\nabla_{TA} u|^2= |\nabla \rho^{1/2}|^2$ and we would have $$\label{introlowerbound}\int_{Q_{L,1}} |\nabla_{TA} u |^2  + \left(B_3 - \a (1-\rho)\right)^2dx=\int_{Q_{L,1}} |\nabla \rho^{1/2} |^2  +  \a^2 \chi_{\rho>0} (1-\rho)^2 +\chi_{\rho=0} (B_3-\a)^2 dx.$$ This is a Modica-Mortola [@ModMort] type of functional with a degenerate double-well potential given by $W(\rho):= \chi_{\rho>0} (1-\rho)^2$. Thanks to Lemma \[lemmameissner\], one can control how far we are from satisfying the Meissner condition. From this, we deduce that almost holds (see Lemma \[lemmafirstlowerbound\]). This implies that the Ginzburg-Landau energy gives a control over the perimeter of the superconducting region $\{\rho>0\}$. In addition, $\b\ll 1$ imposes a small cross-area fraction for $\{\rho>0\}$. Using then isoperimetric effects to get convergence to one-dimensional objects (see Lemma \[lowerbound2d\]), we may use Proposition \[gammaliminf\] to conclude the proof of \[maintheolb\].

In order to prove \[maintheoub\], thanks to the density result in Theorem \[theodens\], it is enough to consider regular measures. Given such a measure $\mu$, we first approximate it with quantized measures (Lemma \[lemmaquantize\]). Far from branching points the construction is easy (see Lemma \[lemmacurve3\]). At a branching point, we need to pass from one disk to two (or vice-versa); this is done passing through rectangles (see Lemma \[Lembranch3\] and Figure \[fig1\]). Close to the boundary we use instead the construction from [@CoOtSer], which explicitly generates a specific branching pattern with the optimal energy scaling; since the height over which this is done is small the prefactor is not relevant here (Proposition \[propboundraylayer\]). The last step is to define a phase and a magnetic potential to get back to the full Ginzburg-Landau functional. This is possible since we made the construction with the Meissner condition and quantized fluxes enforced, see Proposition \[thirdupperbound\].

From and the discussion around , for type-I superconductors, the Ginzburg-Landau functional can be seen as a non-convex, non-local (in $u$) functional favoring oscillations, regularized by a surface term which selects the lengthscales of the microstructures. The appearance of branched structures for this type of problem is shared by many other functionals appearing in material sciences such as shape memory alloys [@KohnMuller92; @KohnMuller94; @Conti00; @KnuepferKohnOtto2013; @BelGol; @Zwicknagl2014; @ChanConti2015; @ContiZwicknagl], uniaxial ferromagnets [@ChoKoOtmicro; @ViehOtt; @KnMu] and blistered thin films [@BCDM00; @JinSternberg; @BCDM02]. Most of the previously cited results on branching patterns (including [@ChokConKoOt; @ChokKoOt; @CoOtSer] for type-I superconductors) focus on scaling laws. Here, as in [@ViehOtt; @CDZ], we go one step further and prove that, after a suitable anisotropic rescaling, configurations of low energy converge to branched patterns. The two main difficulties in our model with respect to the one studied in [@ViehOtt] are the presence of an additional lengthscale (the penetration length) and its sharp limit counterpart, the Meissner condition $\rho B=0$ which gives a nonlinear coupling between $u$ and $B$. Let us point out that for the Kohn-Müller model [@KohnMuller92; @KohnMuller94], a much stronger result is known, namely that minimizers are asymptotically self-similar [@Conti00] (see also [@Viehmanndiss; @AlChokOt] for related results). In [@Gol], the optimal microstructures for a two-dimensional analogue of $I(\mu)$ are exactly computed.

The paper is organized as follows. In section \[sec:not\], we set some notation and recall some notions from optimal transport theory. In Section \[Sec:GL\], we recall the definition of the Ginzburg-Landau functional together with various important quantities such as the superconducting current. We also introduce there the anisotropic rescaling leading to the functional $\widetilde{E}_T$. In Section \[Sec:intermediate\], we introduce for the sake of clarity intermediate functionals corresponding to the different scales of the problem. Let us stress that we will not use them in the rest of the paper but strongly believe that they help understanding the structure of the problem. In Section \[Sec:limiting\], we carefuly define the limiting functional $I(\mu)$ and study its properties. In particular we recover a scaling law for the minimization problem and prove regularity of the minimizers. We then prove the density in energy of ’regular’ measures. This is a crucial result for the main $\Gamma-$ convergence result which is proven in the last two sections. As customary, we first prove the lower bound in Section \[Sec:gammaliminf\] and then make the upper bound construction in Section \[Sec:upperbound\].

 Notation and preliminary results {#sec:not}
=================================

In the paper we will use the following notation. The symbols $\sim$, $\ges$, $\les$ indicate estimates that hold up to a global constant. For instance, $f\les g$ denotes the existence of a constant $C>0$ such that $f\le Cg$, $f\sim g$ means $f\les g$ and $g\les f$. In heuristic arguments we use $a\simeq b$ to indicate that $a$ is close (in a not precisely specified sense) to $b$. We use a prime to indicate the first two components of a vector in $\R^3$, and identify $\R^2$ with $\R^2\times\{0\}{\subseteq}\R^3$. Precisely, for $a\in \R^3$ we write $a'=(a_1,a_2,0)\in\R^2{\subseteq}\R^3$; given two vectors $a,b\in\R^3$ we write $a'\times b'=(a\times b)_3 =
(a'\times b')_3$. We denote by $(e_1,e_2,e_3)$ the canonical basis of $\R^3$. For $L>0$ and $T>0$, $Q_L:=(-\frac{L}{2},\frac{L}{2})^2$ and $Q_{L,T}:=Q_L\times(-T,T)$. For a function $f$ defined on $Q_{L,T}$, we denote $f_{x_3}$ the function $f_{x_3}(x'):=f(x',x_3)$ and we analogously define for $\Omega {\subseteq}Q_{L,T}$, the set $\Omega_{x_3}$. For $x=(x',x_3)$ and $r>0$ we let $\B_r(x)=\B(x,r)$ be the ball of radius $r$ centered at $x$ (in $\R^3$) and $\B'_r(x')=\B'(x',r)$ be the analogue two-dimensional ball centered at $x'$. Unless specified otherwise, all the functions and measures we will consider are periodic in the $x'$ variable, i.e., we identify $Q_L$ with the torus $\R^2/L\Z^2$. In particular, for $x',y'\in Q_L$, $|x'-y'|$ denotes the distance for the metric of the torus, i.e., $|x'-y'|:=\min_{k\in \Z^2} |x'-y'+Lk|$. We denote by $\H^k$ the $k-$dimensional Hausdorff measure. We let $\cP(Q_L)$ be the space of probability measures on $Q_L$ and $\M(Q_{L,T})$ be the space of finite Radon measures on $Q_{L,T}$, and similarly $\M(Q_L)$. Analogously, we define $\M^+(Q_{L,T})$ and $\M^+(Q_L)$ as the spaces of finite Radon measures which are also positive. For a measure $\mu$ and a function $f$, we denote by $f\sharp \mu$ the push-forward of $\mu$ by $f$.

We recall the definition of the (homogeneous) $H^{-1/2}$ norm of a function $f\in L^2(Q_{L})$ with $\int_{Q_{L}} f dx'=0$, $$\label{defHdem}
 \|f\|_{H^{-1/2}}^2:=\inf\lt\{ \int_{Q_{L}\times [0,+\infty)} |B|^2 dx \ : \ \Div B=0, \ B_3(\cdot,0)=f \rt\},$$ which can be alternatively given in term of the two-dimensional Fourier series as (see, e.g., [@ChokKoOt]) $$ \|f\|_{H^{-1/2}}^2=\frac{1}{2\pi} \sum_{ k'\in \lt(\frac{2\pi}{L}\Z\rt)^2\backslash\{0\}} \frac{ |\hat f (k')|^2}{|k'|}.$$ We shall write $\|f\|^2_{H^{-1/2}(Q_L\times\{\pm T\})}$ for $\|f\|^2_{H^{-1/2}(Q_L\times\{T\})}+
 \|f\|^2_{H^{-1/2}(Q_L\times\{-T\})}$.

The $2$-Wasserstein distance between two measures $\mu$ and $\nu\in {\mathcal{M}}^+(Q_L)$ with $\mu(Q_L)=\nu(Q_L)$ is given by $$W_2^2(\mu,\nu):=\min\lt\{ \mu(Q_L)\int_{Q_L\times Q_L} |x-y|^2 \, d\Pi(x,y)  \, : \, \Pi_1=\mu, \ \Pi_2=\nu\rt\},$$ where the minimum is taken over measures on $Q_L\times Q_L$ and $\Pi_1$ and $\Pi_2$ are the first and second marginal of $\Pi$[^5], respectively. For measures $\mu,\nu\in{\mathcal{M}}^+(Q_{L,T})$, the $2$-Wasserstein distance is correspondingly defined. We now introduce some notions from metric analysis, see [@AGS; @villani] for more detail. A curve $\mu:(a,b)\to \cP(Q_L)$, $z\mapsto \mu_z$ belongs to $AC^2(a,b)$ (where $AC$ stands for absolutely continuous) if there exists $m\in L^2(a,b)$ such that $$\label{metricderiv}
W_2(\mu_{z},\mu_{{\widetilde}{z}})\le \int_{z}^{{\widetilde}{z}} m(t) dt \ \qquad \forall a<z\le {\widetilde}{z} <b.$$ For any such curve, [the [*speed*]{}]{} $$|\mu'|(z):=\lim_{{\widetilde}{z}\to z}\frac{W_2(\mu_{{\widetilde}{z}},{\mu}_{z})}{|z-{\widetilde}{z}|}$$ exists for $\H^1-$a.e. $z\in (a,b)$ and $|\mu'|(z)\le m(z)$ for $\H^1$-a.e. $z\in (a,b)$ for every admissible $m$ in . Further, there exists a Borel vector field $B$ such that $$\label{eqexistsB}
 B(\cdot,z)\in L^2(Q_L,\mu_{z}), \qquad \| B(\cdot,z)\|_{L^2(Q_L,\mu_{z})}\le |\mu'|(z) \quad \textrm{for } \H^1\text{-a.e. } \ z\in (a,b)$$ and the continuity equation $$\label{conteqini}
 \partial_3 \mu_{z}+\Divp (B \mu_{z})=0$$ holds in the sense of distributions [@AGS Th. 8.3.1]. Conversely, if a weakly continuous curve $\mu_{z} : (a,b)\to \cP(Q_L)$ satisfies the continuity equation for some Borel vector field $B$ with $\| B(\cdot,z)\|_{L^2( Q_L,\mu_{z})}\in L^2(a,b)$ then $\mu\in AC^2(a,b)$ and $|\mu'|(z)\le \| B(\cdot,z)\|_{L^2(Q_L,\mu_{z})}$ for $\H^1$-a.e. $z\in (a,b)$. In particular, we have $$\label{eqw22b}
W_2^2(\nu,\hat\nu)= \min_{\mu,B} \lt\{ 2T\int_{Q_{L,T}} |B|^2 d\mu_{z} dz \ :\ \mu_{-T}=\nu, \ \mu_T= \hat\nu \textrm{ and \eqref{conteqini} holds}\rt\}\,,$$ where by scaling the right-hand side does not depend on $T$.\
For a (signed) measure $\mu\in \M(Q_L)$, we define the Bounded-Lipschitz norm of $\mu$ as $$\label{BLdefin}
 \|\mu\|_{BL}:=\sup_{\|\psi\|_{Lip}\le 1} \int_{Q_L} \psi d\mu,$$ where for a $Q_L-$periodic and Lipschitz continuous function $\psi$, $\|\psi\|_{Lip}:=\|\psi\|_{\infty}+\|\nabla \psi\|_{\infty}$. By the Kantorovich-Rubinstein Theorem [@villani Th. 1.14], the $1-$Wasserstein and the Bounded-Lipschitz norm are equivalent.

The Ginzburg-Landau functional {#Sec:GL}
==============================

In this section we recall some background material concerning the Ginzburg-Landau functional and introduce the anisotropic rescaling leading to $\wE$.\
For a (non necessarily periodic) function $u:Q_{\Lz,T}\to{\mathbb{C}}$, called the order parameter, and a vector potential $A: Q_{\Lz}\times\R\to \R^3$ (also not necessarily periodic), we define the covariant derivative $$\nabla_A u:= \nabla u-iAu,$$ the magnetic field $$B:=\nabla\times A,$$ and the [superconducting]{} current $$\label{defj}
  j_A:= \frac12\left( -i\bar u (\nabla_A u)+ i u (\overline{\nabla_A  u})\right)=\textrm{Im}(iu \overline{\nabla_A u})\,.$$ Let us first notice that $|\nabla_A u|^2$ and the observable quantities $\rho$, $B$ and $j_A$ are invariant under change of gauge. That is, if we replace $u$ by $u e^{i\varphi}$ and $A$ by $A+\nabla \varphi$ for any function $\varphi$, they remain unchanged. We also point out that if $u$ is written in polar coordinates as $u=\rho^{1/2} e^{i\theta}$, then $$ |\nabla_A u|^2=|\nabla \rho^{1/2}|^2+\rho|\nabla \theta-A|^2.$$ For any admissible pair $(u,A)$, that is such that $\rho$, $B$ and $j_A$ are $Q_{\Lz}$-periodic, we define the Ginzburg-Landau functional as $$  {E_\mathrm{GL}}(u,A):=\int_{Q_{\Lz,T}} |\nabla_A u|^2
+ \frac{\kappa^2}{2} (1-|u|^2)^2 dx
+ \int_{Q_{\Lz}\times \R} |\nabla\times A-{B_{\mathrm{ext}}}|^2 dx\,.$$ We remark that $u$ and $A$ need not be (and, if ${B_{\mathrm{ext}}}\ne 0$, cannot be) periodic. See [@CoOtSer] for more details on the functional spaces we are using. Here ${B_{\mathrm{ext}}}:={b_{\mathrm{ext}}}e_3$ is the external magnetic field and $\kappa\in(0,1/\sqrt2)$ is a [material constant, called the Ginzburg-Landau parameter]{}. From periodicity and $\Div B=0$ it follows that $ \int_{Q_{\Lz}\times\{x_3\}} B_3 \, dx'$ does not depend on $x_3$ and therefore, if the energy is finite, necessarily $$\label{B2bex}
\int_{Q_{\Lz}\times\{x_3\}} B_3 \, dx'=\Lz^2 
{b_{\mathrm{ext}}}\hskip1cm\text{ for all }x_3\in\R\,.$$ We first remove the bulk part from the energy ${E_\mathrm{GL}}$. In order to do so, we introduce the quantity $${\mathcal{D}}_A^3u :=   (\partial_2 u - i A_2 u) - i (\partial_1 u
 - i A_1 u) 
 = (\nabla_A u)_2 - i (\nabla_A u)_1$$ and, more generally, $${\mathcal{D}}_A^ku := (\partial_{k+2} u - i A_{k+2} u) - i (\partial_{k+1} u
 - i A_{k+1} u)  = (\nabla_A u)_{k+2} - i (\nabla_A u)_{k+1} \, ,$$ where components are understood cyclically (i.e., $a_k=a_{k+3}$). The operator ${\mathcal{D}}_A$ (which corresponds to a creation operator for a magnetic Laplacian) was used by Bogomol’nyi in the proof of the self-duality of the Ginzburg-Landau functional at ${\kappa}=\frac{1}{\sqrt{2}} $ (cf. e.g. [@jaffetaubes]). His proof relied on identities similar to the next ones, which will be crucial in enabling us to separate the leading order part of the energy.

Expanding the squares, one sees (for details see [@CoOtSer Lem. 2.1]) that (recall that $\rho:=|u|^2$) $$\label{magic}
  |\nabla'_A u|^2 =  |{\mathcal{D}}_A^3u|^2 + \rho B_3 + \nabla'\times 
   j'_A$$ and, for any $k=1,2,3$, $$\label{magic2}
   |(\nabla_A u)^{k+1}|^2 +    |(\nabla_A u)^{k+2}|^2= |{\mathcal{D}}_A^ku|^2
   + \rho B_k + (\nabla\times j_A)_k  \,.$$ This implies $$|\nabla_A' u |^2 =   (1-{\kappa}\sqrt2) |\nabla_A' u |^2 + {\kappa}\sqrt2 
|{\mathcal{D}}_A^3u|^2 
+{\kappa}\sqrt2\rho B_3 + {\kappa}\sqrt2\nabla'\times j_A'\,.$$ The last term integrates to zero by the periodicity of $j_A$. Therefore, for each fixed $x_3$, using , we have $$\int_{Q_{\Lz}} |\nabla_A' u |^2\,  dx' =    \int_{Q_{\Lz}} (1-{\kappa}\sqrt2) |\nabla_A' u |^2 + {\kappa}\sqrt2 
|{\mathcal{D}}_A^3u|^2  +{\kappa}\sqrt2(\rho-1) B_3 \, dx'+ \Lz^2 {\kappa}\sqrt2 {b_{\mathrm{ext}}}\,.$$ We substitute and obtain, using $\int_{Q_{\Lz}} (B_3-{b_{\mathrm{ext}}})^2 dx'=\int_{Q_{\Lz}} B_3^2-{b_{\mathrm{ext}}}^2 dx'$ [and completing squares]{}, $$\label{eqEGLE}
{E_\mathrm{GL}}(u,A)=2T\Lz^2\left({\kappa}\sqrt2 {b_{\mathrm{ext}}}-{b_{\mathrm{ext}}}^2\right) + E(u,A)
  + {\kappa}\sqrt2
  \int_{Q_{\Lz,T}}|{\mathcal{D}}_A^3u|^2-|\nabla_A'u|^2 dx,$$ where $$\begin{aligned}
  E(u,A)&:=&
\int_{Q_{\Lz,T}} |\nabla_A u |^2  + \left(B_3 - \frac{{\kappa}}{\sqrt2} (1-\rho)\right)^2 dx 
\\&&
+  \int_{Q_{\Lz}\times \R} |B'|^2  dx+
   \int_{Q_{\Lz}\times(\R\setminus(-T,T))} (B_3 -{b_{\mathrm{ext}}})^2 dx\,.\end{aligned}$$ In particular, the bulk energy is $2{\Lz}^2T({\kappa}\sqrt2 {b_{\mathrm{ext}}}-{b_{\mathrm{ext}}}^2) $. Since we are interested in the regime $\kappa\ll 1$ and since $ |{\mathcal{D}}_A^3u|^2\le 2 |\nabla_A'u|^2$, the contribution of the last term in (\[eqEGLE\]) to the energy is (asymptotically) negligible with respect to the first term in $ E$, and therefore it can be ignored in the following.

Applying to $B-{b_{\mathrm{ext}}}e_3  $ and minimizing outside $Q_{\Lz} \times [-T,T]$ if necessary, the last two terms in $E(u,A)$ can be replaced by $$\int_{Q_{{\Lz},T}} |B'|^2 dx+ \| B_3-{b_{\mathrm{ext}}}\|_{H^{-1/2}(\{x_3= T\})}^2+\| B_3-{b_{\mathrm{ext}}}\|_{H^{-1/2}(\{x_3=-T\})}^2\,,$$ so that $ E(u,A)$ becomes $$\label{E}
  E(u,A)=
\int_{Q_{\Lz,T}}  |\nabla_A u |^2 + \left(B_3 - \frac{{\kappa}}{\sqrt2} (1-\rho)\right)^2 +  |B'|^2 dx+\| B_3-{b_{\mathrm{ext}}}\|_{H^{-1/2}(\{x_3=\pm T\})}^2\,.$$

Let us notice that the normal solution $\rho=0$, $B={b_{\mathrm{ext}}}e_3$ (for which we can take $A(x_1,x_2,x_3)={b_{\mathrm{ext}}}x_1 e_2$) is always admissible but has energy equal to $${E_\mathrm{GL}}(u,A)= \Lz^2 T \kappa^2\gg 2\Lz^2T({\kappa}\sqrt2 {b_{\mathrm{ext}}}-{b_{\mathrm{ext}}}^2),$$ in the regime $\kappa\gg {b_{\mathrm{ext}}}$ that we consider here.

The following scaling law is established in [@CoOtSer].

\[theoCoOtSer\] For ${b_{\mathrm{ext}}}<\kappa/8$, $\kappa\le 1/2$, $\kappa
T \ge 1$, $\Lz$ sufficiently large, if the quantization condition $$\label{quantificationini}
 {b_{\mathrm{ext}}}\Lz^2 \in 2\pi \Z,$$ holds then $$\label{minE}
\min E(u,A) \sim \min \left\{{b_{\mathrm{ext}}}^{2/3} \kappa^{2/3} T^{1/3} , {b_{\mathrm{ext}}}\kappa^{3/7} T^{3/7} \right\}\Lz^2.$$

We believe that is also a necessary condition for to hold. Indeed, we expect that if $(u,A)$ is such that $$E(u,A) \les  \left\{{b_{\mathrm{ext}}}^{2/3} \kappa^{2/3} T^{1/3} \Lz^2, {b_{\mathrm{ext}}}\kappa^{3/7} T^{3/7} \Lz^2 \right\},$$ then the normal phase $\rho\simeq 0$ is the minority phase (typically disconnected on every slice) and there exist $x_3\in(-T,T)$ and (periodic) curves $\gamma_1$ and $\gamma_2$ such that $\Gamma_1:=\{(\gamma_1(s), s, x_3) : s\in[ 0, 1]\} {\subseteq}\{\rho\simeq 1\}$ and $\Gamma_2:=\{(s,\gamma_2(s), x_3) : s\in[ 0, 1]\} {\subseteq}\{\rho\simeq 1\}$, with $\gamma_i(1)=\gamma_i(0)+L_0 e_i$. If this holds then using Stokes Theorem on large domains the boundary of which is made of concatenations of the curves $\Gamma_i$, it is possible to prove that must hold. As in [@CoOtSer], we will need to assume in order to build the recovery sequence in Section \[recoveryGL\].\
The first regime in corresponds to uniform branching patterns while the second corresponds to well separated branching trees (see [@ChokConKoOt; @ChokKoOt; @CoOtSer]). We focus here on the first regime, that is ${\kappa}^{5/7}\ll {b_{\mathrm{ext}}}T^{2/7}$, and replace $\kappa$ and ${b_{\mathrm{ext}}}$ by the variables $\alpha$, $\beta$, defined according to $$\kappa T = \sqrt{2}\a \qquad {b_{\mathrm{ext}}}= \frac{\beta \kappa }{\sqrt{2}}=\frac{\alpha\beta}{T}\,,$$ and then rescale

  ------------------------------------------------------- ---------------------------
  $\hat{x}:= x/T$,                                        $  L:=\Lz/T$,
  \[6pt\] $\hat{u}(\hat{x}):= u(x)$,                      $\hat{A}(\hat{x}):= A(x)$
  \[6pt\] $E_T(\hat u,\hat A):=\frac{1}{T L^2} E(u,A)$,   
  ------------------------------------------------------- ---------------------------

so that in particular $ \hat{B}(\hat{x})= \widehat{{\nabla}}\times \hat{A}(\hat{x})= TB(x)$ and ${\nabla}_A u(x)= T^{-1} \widehat{{\nabla}}
\hat{u}( x/T)-i \hat{A}(x/T)\hat u (x/T)$. Changing variables and removing the hats yields $$E_T(u,A)=  \frac{1}{L^2}\lt(\int_{Q_{L,1}}
|{\nabla}_{T{A}} {u}|^2+ \left({B_3} - \a(1-\rho)\right)^2 +  |{B'}|^2 dx+\|{B_3}- \a\b\|_{H^{-1/2}(\{x_3=\pm1\})}^2\rt)\,.$$ as was anticipated in (\[Etilde\]). In these new variables, the scaling law becomes $E_T\sim \min\{\alpha^{4/3}\beta^{2/3},\a^{10/7}\beta \}$ and the uniform branching regime corresponds to $E_T\sim \alpha^{4/3}\beta^{2/3}$ which amounts to $\a^{-2/7}\ll \b\ll 1$, see also (\[scallawintro\]). Constructions (leading to the upper bounds in [@CoOtSer; @ChokKoOt; @ChokConKoOt]), suggest that in this regime, typically, the penetration length of the magnetic field inside the superconducting regions is of the order of $T^{-1}$, the coherence length (or domain walls) is of the order of $\a^{-1}$, the width of the normal domains in the bulk is of the order of $\a^{-1/3} \b^{1/3}$ and their separation of order $\a^{-1/3}\b^{-1/6}$. These various lengthscales motivate the anisotropic rescalings that we will introduce in Section \[Sec:intermediate\].

In closing this section we present the anisotropic rescaling that will lead to the functional defined in (\[eqdefwe\]), postponing to the next section a detailed explanation of its motivation. We set for $x \in Q_{L,1}$,

  --------------------------------------------------------------------------------------------------------------------------------------------------------------------------- --------------------------------------------
  $\begin{pmatrix}{\widetilde}{x}'\\ {\widetilde}{x}_3\end{pmatrix}:= \begin{pmatrix}\a^{1/3} \b^{1/6} x'\\x_3\end{pmatrix}$,                                                 $ {\widetilde}L:= \a^{1/3}\b^{1/6}L$,
  \[8pt\] $\begin{pmatrix}{\widetilde}{A}'\\ {\widetilde}{A}_3\end{pmatrix}({\widetilde}{x}):=\begin{pmatrix}\a^{-2/3}\b^{1/6} A'\\ \a^{-1/3} \b^{1/3}A_3\end{pmatrix}(x)$,   $ {\widetilde}{u}({\widetilde}{x}):=u(x),$
  --------------------------------------------------------------------------------------------------------------------------------------------------------------------------- --------------------------------------------

to get ${\widetilde}B_3({\widetilde}x)=\alpha^{-1}B_3(x)$, ${\widetilde}B'({\widetilde}x)=\alpha^{-2/3}\beta^{1/6}B'(x)$ inside the sample. Outside the sample, i.e. for $|x_3|\ge 1$, we make the isotropic rescaling ${\widetilde}{x}:= \pm e_3+\a^{1/3}\b^{1/6}( x \mp e_3)$ to get ${\widetilde}B({\widetilde}x)=\a^{-1} B(x)$. A straightforward computation leads to $\wE({\widetilde}u, {\widetilde}A)=\a^{-4/3}\b^{-2/3} E_T(u,A)$, where $$\begin{gathered}
  \wE(u,A):=\frac{1}{\widetilde{L}^2}\Big[\int_{Q_{{\widetilde}L,1}} \a^{-2/3}\b^{-1/3}\lt|\nabla_{\a^{1/3}\b^{-1/3}TA}'u\rt|^2+ 
  \a^{-4/3}\b^{-2/3}\lt|(\nabla_{\a^{1/3}\b^{-1/3}TA} u)_3\rt|^2\\
+  \a^{2/3}\b^{-2/3}\left({B_3} -(1-|{u}|^2)\right)^2
+  \b^{-1}|{B'}|^2 dx +\a^{1/3}\b^{7/6}\|\b^{-1}{B_3}-1\|_{H^{-1/2}(x_3=\pm1)}^2\Big],\end{gathered}$$ with ${\nabla}\times A=B$ (and in particular $\Div B=0$). We assume that ${\widetilde}L$ is a fixed quantity of order 1. For simplicity of notation, the detailed analysis is done only for the case ${\widetilde}L=1$.\
Let us point out that in these units, the penetration length is of order $T^{-1}\a^{1/3}\b^{1/6}$, the coherence length of order $\a^{-2/3}\b^{1/6}$, the width of the normal domains in the bulk of order $\b^{1/2}$ and the distance between the threads of order one. That is, the scale separation reads now $$\label{scalesepar2}
 T^{-1}\a^{1/3}\b^{1/6}\ll \a^{-2/3}\b^{1/6}\ll \b^{1/2}\ll1.$$

The intermediate functionals {#Sec:intermediate}
============================

In this section we explain the origin of the rescaling leading from $E_T$ to ${\widetilde}E_T$, and the different functionals which appear at different scales. This material is not needed for the proofs but we think it is important to illustrate the meaning of our results. We carry out the scalings in detail but the relations between the functionals are here discussed only at a heuristic level.

We want to successively send $T\to+\infty$, $\alpha\to +\infty$ and $\beta\to 0$. For this we are going to introduce a hierarchy of models starting from $E_T(u,A)$ and finishing at $I(\mu)$. When sending first $T\to +\infty$ with fixed $\alpha$ and $\beta$, the functional $E_T$ approximates $$ F_{\alpha,\beta}(\rho,B):=\frac{1}{L^2} \lt(\int_{Q_{L,1}} |\nabla \rho^{1/2}|^2 +(B_3-\alpha (1-\rho))^2 +|B'|^2 dx +\|{B_3}- \a\b\|_{H^{-1/2}(x_3=\pm1)}^2\rt),$$ with the constraints $$\label{meissner4}
 \Div B=0 \qquad\textrm{and } \qquad \rho B=0.$$ The main difference between $E_T$ and $F_{\alpha,\beta}$ is that for the latter, since the penetration length (which corresponds to $T^{-1}$) was sent to zero, the Meissner condition is enforced. We now want to send the coherence length (of order $\a^{-1}$) to zero at fixed $\beta$, while keeping superconducting domains of finite size. Since the typical domain diameter is of order $\a^{-1/3}\b^{1/3}$ and their distance is of order $\a^{-1/3}\b^{-1/6}$, we are led to the anisotropic rescaling:

  -------------------------------------------------------------------------------------------------------------------------------- ---------------------------------
  $\begin{pmatrix}\hat{x}'\\ \hat{x}_3\end{pmatrix}:= \begin{pmatrix}\a^{1/3} x'\\x_3\end{pmatrix}$,                               $\hat{L}:= \a^{1/3}L$,
  \[6pt\] $\begin{pmatrix}\hat{B}'\\ \hat{B}_3\end{pmatrix}(\hat{x}):=\begin{pmatrix}\a^{-2/3} B'\\ \a^{-1}B_3\end{pmatrix}(x)$,   $\hat{\rho}(\hat{x}):=\rho(x)$,
  \[6pt\] $\wF=\a^{-4/3} F_{\a,\b}$.                                                                                               
  -------------------------------------------------------------------------------------------------------------------------------- ---------------------------------

In these variables, the coherence length is of order $\a^{-2/3}\ll 1$ (at least horizontally) while the diameter of the normal domains is of order $\beta^{1/3}$ and their separation of order $\beta^{-1/6}$. Dropping the hats (we just keep them on the functional and on $L$ to avoid confusion) we obtain $$\begin{gathered}
 \wF(\rho,B):=\frac{1}{\hL^2} \lt(\int_{Q_{\hL,1}} \a^{-2/3}\left|\begin{pmatrix} \nabla' \rho^{1/2}\\ \a^{-1/3} \partial_3 \rho^{1/2}  
                                                                           \end{pmatrix}\right|^2
 +\a^{2/3}\lt(B_3- (1-\rho)\rt)^2 +|B'|^2 dx\rt.\\
\lt. +\a^{1/3}\|{B_3}- \b\|_{H^{-1/2}(x_3=\pm1)}^2\rt),\end{gathered}$$ with the constraints . The scaling of Theorem \[theoCoOtSer\] indicates that $\wF$ behaves as $ \min\{\beta^{2/3},\a^{2/21}\beta \}$ which is of order $\beta^{2/3}$ if $\a\gg 1$ and $\beta$ is fixed. We remark that, letting $\eta:= \alpha^{-1/3}$ and $\delta:=\eta^2=\alpha^{-2/3}$, one has $$\begin{gathered}
 \wF(\rho,B)= \frac{1}{\hL^2} \left(\int_{Q_{\hL,1}} \delta \lt|\begin{pmatrix}\nabla'\rho^{1/2}\\ \eta \partial_3 \rho^{1/2}\end{pmatrix}\rt|^2+\frac{1}{\delta} \lt(B_3-(1-\rho)\rt)^2 
+|B'|^2 dx \rt.\\
\lt.+ \eta^{-1} \|B_3-\beta\|_{H^{-1/2}(x_3=\pm1)}^2\rt).\end{gathered}$$ In this form, $ \wF(\rho,B)$ is very reminiscent of the functional studied in [@ViehOtt]. Notice however that besides the Meissner condition which makes our functional more rigid, the scaling $\delta=\eta^2$ is borderline for the analysis in [@ViehOtt].

Recalling that $\lt(B_3-(1-\rho)\rt)^2 = \chi_{\rho>0} (1-\rho)^2$, the corresponding term in $\wF$ has the form of a double well-potential, and so in the limit $\alpha\to+\infty$ the functional $\wF$ approximates $$ G_\beta(\chi,B'):=\frac{1}{\hL^2} \lt( \int_{Q_{\hL,1}} \frac{4}{3}|D'\chi| +|B'|^2 dx \rt),$$ with the constraints $\chi\in \{0,1\}$, $\chi(\cdot,x_3) \weaklim \beta dx'$ when $x_3\to \pm1$ and $$ \partial_3 \chi +\Divp B'=0 \qquad \textrm{and } \qquad \chi B'=B'.$$ This is similar to the simplified sharp-interface functional that was studied in [@ChokKoOt; @ChokConKoOt]. In the definition of $G_\beta$, we used the notation $$ \int_{Q_{\hL,1}} |D' u|:= \sup_{\stackrel{\xi \in C^{\infty}(Q_{\hL,1}),}{ |\xi|_\infty \le 1}}  \, \int_{Q_{\hL,1}} u \Divp \xi \, dx,$$ for the horizontal $BV$ norm of a function $u\in L^1(Q_{\hL,1})$. By definition it is lower semicontinuous for the $L^1$ convergence and it is not hard to check that if we let $u_{x_3}(x'):= u(x',x_3)$, then $$\int_{Q_{\hL,1}} |D'u|=\int_{-1}^1 \left(\int_{Q_{\hL}} |D' u_{x_3}|\right) dx_3,$$ where $\int_{Q_{\hL}} |D' u_{x_3}|$ is the usual $BV$ norm of $u_{x_3}$ in $Q_{\hL}$ [@AFP]. From this and the usual co-area formula [@AFP Th. 3.40], we infer that $$\label{coarea}
 \int_{Q_{\hL,1}} |D'u|=\int_{\R} \int_{-1}^1 \H^1( \partial \{u_{x_3}>s\}) dx_3 ds.$$ In , $\partial \{u_{x_3}>s\}$ represents the measure-theoretic boundary of $\{u_{x_3}>s\}$ in $Q_{\hL}$.

We finally want to send the volume fraction of the normal phase to zero and introduce the last rescaling in $\beta$ for which we let

  -------------------------------------------------------------------------------------------------------------------- ----------------------------------------------
  $\begin{pmatrix}{\widetilde}{x}'\\ {\widetilde}{x}_3\end{pmatrix}:= \begin{pmatrix}\b^{1/6} x'\\x_3\end{pmatrix}$,   $ 
                                                                                                                        {\widetilde}L:= \beta^{1/6}\hL$,

  \[8pt\] $\tB'({\widetilde}{x}):=\beta^{1/6} B'(x)$                                                                   $\tchi({\widetilde}{x}):=\beta^{-1}\chi(x)$,

  \[6pt\] $\wG:=\b^{-2/3} G_{\b}$.                                                                                     
  -------------------------------------------------------------------------------------------------------------------- ----------------------------------------------

After this last rescaling, the domain width is of order $\b^{1/2}\ll1$, the separation between domains of order $1$. We obtain (dropping the tildas again) the order-one functional $$ \wG(\chi,B'):= \frac{1}{\widetilde{L}^2}\int_{Q_{{\widetilde}L,1}} \frac{4}{3} \beta^{1/2}|D'\chi| +\chi|B'|^2 dx$$ under the constraints $\chi\in \{0,\beta^{-1}\}$, $\chi(\cdot,x_3) \weaklim  dx'$ when $x_3\to \pm1$, and $$ \partial_3 \chi +\Divp (\chi B')=0 \qquad \textrm{and } \qquad \chi B'=\beta^{-1}B'.$$ This functional converges to $I(\mu)$ as $\beta\to 0$.

Let us point out that since we are actually passing directly from the functional $\wE$ to $I$ in Theorem \[maintheo\], we are covering the whole parameter regime of interest. In particular, our result looks at first sight stronger than passing first from $ E_T$ to $F_{\a,\b}$, then from $\wF$ to $G_\beta$ and finally from $\wG$ to $I$. However, because of the Meissner condition, we do not have a proof of density of smooth objects for $F_{\a,\b}$ and $G_\b$. Because of this, we do not obtain the $\Gamma-$convergence of the intermediate functionals (the upper bound is missing).

The limiting energy {#Sec:limiting}
===================

Before proving the $\Gamma$-limit we study the limiting functional $I$ that was mentioned in (\[eqdefIintro\]) and motivated in the previous section. We give here a self-contained treatment of the functional $I$, which is motivated by the analysis discussed above, and will be crucial in the proofs that follow. However, in this discussion we do not make use of the relation to the Ginzburg-Landau functional.

For $L,T>0$ we denote by ${\mathcal{A}}_{L,T}$ the set of pairs of measures ${\mu}\in{\mathcal{M}}^+(Q_{L,T})$, $m\in {\mathcal{M}}(Q_{L,T};\R^2)$ with $m\ll\mu$, satisfying the continuity equation $$\label{conteq0}
 \pa_3 {\mu}+\Divp m=0 \qquad \textrm{in } Q_{L,T},$$ and such that ${\mu}={\mu}_{x_3}\otimes d{x_3}$ where, for a.e. $x_3\in(-T,T)$, ${\mu}_{x_3}=\sum_i {\varphi}_i \delta_{X_i}$ for some ${\varphi}_i> 0$ and $X_i\in Q_L$. We denote by ${\mathcal{A}}_{L,T}^*:=\{\mu: \exists m, (\mu,m)\in{\mathcal{A}}_{L,T}\}$ the set of admissible $\mu$.

We define the functional $I:{\mathcal{A}}_{L,T}\to[0,+\infty]$ by $$\label{Imum}
 I(\mu,m):=
   \frac{K_*}{L^2} \int_{-T}^T  \sum_{x'\in Q_L} \left(\mu_{x_3}(x')\right)^{1/2}  \, dx_3 + 
   \frac{1}{L^2}\int_{Q_{L,T}}
   \left(\frac{dm}{d{\mu}} \right)^2 d{\mu},$$ where $K_*:= \frac{8 \sqrt{\pi}}{3}$ and (with abuse of notation) $I:{\mathcal{A}}_{L,T}^*\to[0,\infty]$ by $$\label{Imu}
I({\mu}):=\min \{ I(\mu,m)\ : \ m\ll{\mu}, \ \pa_3 {\mu}+ {\Divp} m=0\}.$$

Condition (\[conteq0\]) is understood in a $Q_L$-periodic sense, i.e., for any $\psi\in C^1(\R^3)$ which is $Q_L$-periodic and vanishes outside $\R^2\times (0,T)$ one has $\int_{Q_{L,T}} \partial_3\psi d\mu + \nabla'\psi \cdot dm=0$. If $\mu_{x_3}=\sum_i {\varphi}_i \delta_{X_i}$ then $\sum_{x'\in Q_L} \left(\mu_{x_3}(x')\right)^{1/2}=\sum_i {\varphi}_i^{1/2}$. Because of , $\mu_{x_3}(Q_L)$ does not depend on $x_3$.

Let us point out that the minimum in (\[Imu\]) is attained thanks to (\[eqexistsB\]). Moreover, the minimizer is unique by strict convexity of $m\to\int_{Q_{L,T}} \left(\frac{dm}{d{\mu}} \right)^2 d{\mu}$. As proven in Lemma \[lemmacurves\] below, if $\mu$ is made of a finite union of curves then there is actually only one admissible measure $m$ for . More generally, Since every measure $\mu$ with finite energy is rectifiable (see Corollary \[coroXia\]), we believe that it is actually always the case. For $\mu$ an admissible measure and $z,{\widetilde}{z}\in[-T,T]$, we let $$\label{Iztz}
 I^{(z,{\widetilde}{z})}(\mu):= \frac{K_*}{L^2} \int_{z}^{{\widetilde}{z}}  \sum_{x'\in Q_L} \left(\mu_{x_3}(x')\right)^{1/2}  \, dx_3 + 
   \frac{1}{L^2}\int_{Q_{L}\times[z,{\widetilde}{z}]}
   \left(\frac{dm}{d{\mu}} \right)^2 d{\mu},$$ where $m$ is the optimal measure for ${\mu}$ on $[z,{\widetilde}{z}]$ (which coincides with the restriction to $[z,{\widetilde}{z}]$ of the optimal measure on $[-T,T]$).

From (\[eqw22b\]) one immediately deduces for every measure $\mu$, and every $x_3, {\widetilde}x_3\in[-T,T]$, the following estimate on the Wasserstein distance $$\label{HolderW2}
  W_2^2(\mu_{x_3}, {\mu}_{{\widetilde}x_3})\le L^2 I(\mu) |x_3-{\widetilde}x_3|.$$ In particular for every measure $\mu$ with $I(\mu)<+\infty$, the curve $x_3\mapsto \mu_{x_3}$ is Hölder continuous with exponent $1/2$ in ${\mathcal{M}}^+( Q_L)$ (endowed with the metric $W_2$) and the traces $\mu_{\pm T}$ are well defined.

Existence of minimizers
-----------------------

Given two measures $\bar {\mu}_{\pm}$ in $\M^+(Q_L)$ with $\bar \mu_+(Q_L)=\bar \mu_-(Q_L)$, we are interested in the variational problem $$\label{limitProb}
 \inf\lt\{ I(\mu) \ : \ \mu_{\pm T}= \bar {\mu}_{\pm} \rt\}.$$ We first prove that any pair of measures with equal flux can be connected with finite cost and that there always exists a minimizer. The construction is a branching construction which gives the expected scaling (see [@ChokConKoOt; @CoOtSer]) if the boundary data is such that $\bar {\mu}_+=\bar{\mu}_-$.

\[branchingmu\] For every pair of measures $ \bar {\mu}_{\pm}\in {\mathcal{M}}^+(Q_L)$ with $\bar\mu_+(Q_L)=\bar\mu_-(Q_L)=\Phi$, there is $\mu\in {\mathcal{A}}^*_{L,T}$ such that $\mu_{\pm T}=\bar\mu_\pm$ and $$I(\mu)\les  \frac{T\Phi^{1/2}}{L^2}+\frac{\Phi}{T}.$$ If $\bar\mu_+=\bar\mu_-$, then there is a construction with $$I(\mu)\les  \frac{T \Phi^{1/2}}{L^2} +\frac{  T^{1/3} \Phi^{2/3} }{L^{4/3}},$$ and such that the slice at $x_3=0$ is given by $\mu_0=\Phi N^{-2}\sum_j \delta_{X_j}$, with $X_j$ the $N^2$ points in $[-L/2,L/2)^2\cap ((L/N)\Z)^2$, and $N:=\lfloor1+ \Phi^{1/6}L^{2/3}/T^{2/3}\rfloor$. The measure $\mu$ is supported on countably many segments, which only meet at triple points.

By periodicity we can work on $[0,L)^2$ instead of $[-L/2,L/2)^2$. We first perform the construction for $x_3\ge 0$. The idea is to approximate $\bar\mu_+$ by linear combinations of Dirac masses, which become finer and finer as $x_3$ approaches $T$. Fix $N\in\N$, chosen below. For $n\in\N$, fix $x_{3,n}:=T(1-3^{-n})$, and let $T_n:= x_{3,n}-x_{3,n-1}=\frac{2}{3^{n}}T$ be the distance between two consecutive planes. At level $x_{3,n}$ we partition $Q_L$ into squares of side length $L_n:= \frac{L}{2^nN}$. More precisely, for $i,j=0,...,2^nN-1 $, we let $x'_{ij,n}:=\lt(L_n\, i, L_n \, j\rt)$ be a corner of the square $Q_{ij,n}:= x'_{ij,n}+ [0,L_n)^2$, and we let $\Phi_{ij,n}:=\bar\mu_+(Q_{ij,n})$ be the flux associated to this square.

We define the measures $\mu^\mathrm{br}$ and $m^\mathrm{br}$ (here the suffix $\mathrm{br}$ stands for branching) by $$\mu^\mathrm{br}_{x_3}:= \sum_{ij} \Phi_{ij,n} \delta_{X_{ij,n}(x_3)} \text{ and }
m^\mathrm{br}_{x_3}:= \sum_{ij} \frac{d X_{ij,n}}{dx_3}(x_3) \Phi_{ij,n} \delta_{X_{ij,n}(x_3)} \quad 
\text{for }  x_3\in [x_{3,n-1}, x_{3,n}),$$ where $X_{ij,n}:[x_{3,n-1}, x_{3,n}]\to Q_L$ is a piecewise affine function such that $X_{ij,n}(x_{3,n})=x'_{ij,n}$, $X_{ij,n}(x_{3,n}-\frac12 T_n)=x'_{i_*j,n}$, and $X_{ij,n}(x_{3,n-1})=x'_{i_*j_*,n}$, where $i_*=2\lfloor i/2\rfloor$, $j_*=2\lfloor j/2\rfloor$. Four such curves end in every $i_*$, $j_*$ (which corresponds to the pair $i_*/2$, $j_*/2$ at level $n-1$), but they are pairwise superimposed for $x_3\in [x_{3,n-1},x_{3,n}-\frac12 T_n]$, therefore all junctions are triple points (one curve goes in, two go out).

Using that $\sum_{i j} \Phi_{ij,n}=\Phi$ and $\sum_{ij}\sqrt{\Phi_{ij,n}} \le
(\sum_{ij} \Phi_{ij,n})^{1/2}(\sum_{ij} 1)^{1/2} = \Phi^{1/2} 2^nN$, we get that the energy of $\mu^\mathrm{br}$ is given by $$\begin{aligned}
 I(\mu^\mathrm{br})&=\frac{1}{L^2} \sum_{n=1}^{+\infty} \sum_{ij} \left( K_* T_n   \sqrt{\Phi_{ij,n}}+ \Phi_{ij,n} T_n \frac{2 L_n^2}{ T_n^2}\rt)\\
& \les  L^{-2} TN\Phi^{1/2} \sum_{n=0}^{+\infty} \lt(\frac{2}{3}\rt)^n + \frac{\Phi}{TN^2}\sum_{n=0}^{+\infty} \lt(\frac{3}{4}\rt)^n\,.\end{aligned}$$ If we choose $N=1$, then there is only one point in the central plane, $\mu_0=\Phi \delta_0$. Therefore the top and bottom constructions can be carried out independently, since by assumption the total flux is conserved, and we obtain the first assertion.

If $\bar\mu_+=\bar\mu_-$, we can choose the value of $N$ which makes the energy minimal. Up to constants this is the value given in the statement. Inserting in the estimate above gives the second assertion.

If the boundary densities are maximally spread, in the sense that they are given by the Lebesgue measure, the scaling is optimal, as the following lower bound shows.

\[branchinglowerbound\] For every measure $\mu\in {\mathcal{A}}^*_{L,T}$ such that $\mu_{\pm T}=\Phi L^{-2} dx'$ one has $$\label{lower}
  I(\mu)\ges  \frac{T \Phi^{1/2}}{L^2}+ \frac{ T^{1/3}\Phi^{2/3}}{ L^{4/3}}\, .$$

The bound $I(\mu)\ge L^{-2} T\Phi^{1/2}$ follows at once from the subadditivity of the square root. Hence we only need to prove the other one. We give two proofs of this bound. The first uses only elementary tools while the second is based on an interpolation inequality.

[*First proof:*]{} Let $I:=I(\mu,v\mu)$, where $v:=dm/d\mu$. Fix $\lambda>0$, chosen below. Choose $x_3\in (-T,T)$ such that $\mu_{x_3}=\sum_i\varphi_i \delta_{X_i}$ obeys $$\label{eqlbchoicex3}
  \sum_i \varphi_i^{1/2} \le \frac{L^2 I}{T}\,.$$ For some set ${\mathcal{I}}{\subseteq}\N$ to be chosen below, let $\psi:Q_L\to \R$ be a mollification of the function $\max\{ (\lambda-\dist(x',X_i))_+: i\in {\mathcal{I}}\}$, where as usual the distance is interpreted periodically. By the divergence condition, $$\int_{Q_L} \psi d\mu_{x_3} = \frac{\Phi}{L^{2}} \int_{Q_L} \psi dx'+\int_{-T}^{x_3} \int_{Q_L} \nabla'\psi\cdot  v d\mu$$ (to prove this, pick $\xi_\eps\in C^1_c((-T,x_3))$ which converge pointwise to $1$ and use $\xi_\eps\psi$ as a test function in (\[conteq0\]) and then pass to the limit). Since $|\nabla'\psi|\le 1$, $$\sum_{i\in {\mathcal{I}}} \lambda \varphi_i \le  \frac{\Phi}{L^{2}} \sum_{i\in {\mathcal{I}}}\frac{ \pi}{3} \lambda^3+ \int_{-T}^{x_3} \int_{Q_L}  |v| d\mu
    \le \frac{\Phi}{L^2}\sum_{i\in {\mathcal{I}}} \pi \lambda^3+  L (T\Phi)^{1/2} I^{1/2},$$ where in the second step we used Hölder’s inequality and flux conservation. We choose ${\mathcal{I}}=\{i\in \N: \varphi_i\ge 4 \Phi \lambda^2/L^2\}$. From the definition of ${\mathcal{I}}$, we have $$\frac{\Phi}{L^2}\sum_{i\in {\mathcal{I}}} \pi \lambda^3\le \frac{\pi}{4} \sum_{i\in {\mathcal{I}}} \lambda {\varphi}_i.$$ Therefore, since $\pi<4$, we obtain $$\label{E1}
    \sum_{i\in {\mathcal{I}}} \lambda \varphi_i {\stackrel{<}{\sim}}L(T\Phi I)^{1/2}\,.$$ At the same time, again by the definition of ${\mathcal{I}}$ and (\[eqlbchoicex3\]), $$\label{E2}
  \sum_{i\not\in {\mathcal{I}}} \varphi_i\le \frac{2\lambda\Phi^{1/2}}{L} \sum_{i\not\in {\mathcal{I}}} \varphi_i^{1/2} \le 2\lambda L \Phi^{1/2}\, \frac{I}{T}\,.$$Adding and , we obtain $$\sum_{i\in {\mathcal{I}}} \varphi_i {\stackrel{<}{\sim}}\frac{1}{\lambda} L(T\Phi I)^{1/2}+ \lambda L \Phi^{1/2}\, \frac{I}{T}\, ,$$ hence $$\Phi^{1/2}{\stackrel{<}{\sim}}\frac{1}{\lambda} L(T I)^{1/2}+ \lambda L \, \frac{I}{T}\, ,$$ and optimizing over $\lambda$ by choosing $\lambda = \frac{T^{3/4}}{I^{1/4}}$ yields $I\ges \Phi^{2/3} T^{1/3}L^{-4/3}$.

[*Second Proof:*]{} As before, let $x_3\in (-T,T)$ be such that $\mu_{x_3}=\sum_i\varphi_i \delta_{X_i}$ obeys . By Young’s inequality and , we have $$I\ges  \frac{T}{L^2} \left( \sum_i \varphi_i^{1/2} \right) +\frac{W_2^2(\mu_{x_3}, \Phi L^{-2} dx')}{L^2 T}\ges L^{-2} T^{1/3}\left( \sum_i \varphi_i^{1/2} \right)^{2/3} \left(W_2^2(\mu_{x_3}, \Phi L^{-2} dx')\right)^{1/3}.$$ The desired lower bound would then follow if we can show that for every measure $\mu\in\M^+(Q_L)$ with $\mu=\sum_i {\varphi}_i \delta_{X_i}$ and $\sum_i {\varphi}_i = \Phi$, $$\label{interpolBVW2}
  \left( \sum_i \varphi_i^{1/2} \right)^{2/3} \left(W_2^2(\mu, \Phi L^{-2} dx')\right)^{1/3}\ges \Phi^{2/3} L^{-4/3}.$$ By rescaling it is enough considering $\Phi=L=1$. The optimal transport map is necessarily of the form $\psi(x')=X_i$ if $x'\in E_i$, where $E_i$ is a partition of $Q_1$ with $|E_i|={\varphi}_i$ (the corresponding transport plan is $(Id\times \psi)\sharp dx'$). By definition, it holds $$W_2^2(\mu, dx')\ge \sum_{i} \int_{E_i} |x'-X_i|^2 dx'.$$ But since $|E_i|={\varphi}_i=|\B'(X_i,({\varphi}_i/\pi)^{1/2})|$, $$\sum_{i} \int_{E_i} |x'-X_i|^2 dx'\ge \sum_{i} \int_{\B(X_i,({\varphi}_i/\pi)^{1/2})} |x'|^2 dx'\ges \sum_{i} {\varphi}_i^2 .$$ By Hölder’s inequality, we conclude that $$1=\sum_i {\varphi}_i \le \left(\sum_i {\varphi}_i^{1/2}\right)^{2/3} \left(\sum_i {\varphi}_i^2\right)^{1/3}\les \left( \sum_i \varphi_i^{1/2} \right)^{2/3} \left(W_2^2(\mu, dx')\right)^{1/3},$$ as desired.

The lower bound can also be obtained as a consequence of the scaling law proven in [@CoOtSer] for the Ginzburg-Landau model combined with our lower bound in Section \[Sec:gammaliminf\] (which does not use this lower bound). However, since the proof here is much simpler and contains some of the main ideas behind the proofs of [@ChokConKoOt; @CoOtSer], we decided to include it. Similarly, the interpolation inequality can be obtained by approximation from a similar inequality proven in [@CintiOt] (where it is used in the same spirit as here to re-derive the lower bounds of [@ChokConKoOt]).

We end this section by proving existence of minimizers.

\[existmu\] For every pair of measures $ \bar {\mu}_{\pm }$ with $\bar {\mu}_+(Q_L)=\bar {\mu}_-(Q_L)$, the infimum in is finite and attained.

In this proof we assume $L=T=1$ and $\bar\mu_+(Q_{1})=1$. By Proposition \[branchingmu\] the infimum is finite. Let now $\mu^n$ be a minimizing sequence for $I$. Since $\sup_n  I(\mu^n)<+\infty$, thanks to , the functions $x_3\mapsto \mu^n_{x_3}$ are equi-continuous in $\cP(Q_1)$ (recall that $W_2$ metrizes the weak convergence in $\mathcal{P}(Q_1)$) hence by the Arzelà-Ascoli theorem there exists a subsequence, still denoted $\mu^n$, uniformly converging (in $x_3$) to some measure $\mu$ which also satisfies the given boundary conditions. Moreover, if $m^n$ is an optimal measure in for $\mu^n$, since by the Cauchy-Schwarz inequality we have $$\int_{Q_{1,1}} d|m^n|\le \lt(\int_{Q_{1,1}} \lt(\frac{dm^n}{d\mu^n}\rt)^2 d\mu^n\rt)^\hal \lt(\int_{Q_{1,1}} d\mu^n\rt)^\hal\les 1,$$ there also exists a subsequence $m^n$ converging to some measure $m$ satisfying . By [@AFP Th. 2.34 and Ex. 2.36] we deduce that $m\ll \mu$ and $$\liminf_{n\to +\infty} \int_{Q_{1,1}} \lt(\frac{dm^n}{d\mu^n}\rt)^2 d\mu^n\ge \int_{Q_{1,1}} \lt(\frac{dm}{d\mu}\rt)^2 d\mu.$$ It remains to prove that $\mu_{x_3}=\sum_i {\varphi}_i(x_3) \delta_{X_i(x_3)}$ for a.e. $x_3$ and that $$\label{semicontperi}
 \liminf_{n\to +\infty} \int_{-1}^1 \sum_{x'\in Q_{1}} \left(\mu^n_{x_3}(x')\right)^{1/2} \, dx_3 \ge \int_{-1}^1 \sum_{x'\in Q_{1}} \left(\mu_{x_3}(x')\right)^{1/2} \, dx_3.$$ If $\mu^n_{x_3}=\sum {\varphi}_i^n(x_3) \delta_{X^n_i}(x_3)$, with ${\varphi}_i^n$ ordered in a decreasing order, we let $f_n(x_3):= \sum_i \sqrt{{\varphi}_i^n(x_3)}$ and observe that $\int_{-1}^1 f_n dx_3 \le I(\mu^n)\les 1 $. Hence, by Fatou’s lemma, $$\label{fatouf}
1\ges \liminf_{n\to +\infty} \int_{-1}^1 f_n(x_3) dx_3\ge \int_{-1}^1 \liminf_{n\to +\infty} f_n(x_3) dx_3,$$ from which we infer that $g(x_3):=\liminf_{n\to+\infty} f_n(x_3)$ is finite for a.e. $x_3$. Consider such an $x_3$ and let $\psi(n)$ be a subsequence (which depends on $x_3$) such that $g(x_3)=\lim_{n\to+\infty} f_{\psi(n)}(x_3)$. Up to another subsequence, still denoted $\psi(n)$, we may assume that for every $i\in \N$, ${\varphi}_i^{\psi(n)}(x_3)$ converges to some ${\varphi}_i(x_3)$ and $X_i^{\psi(n)}(x_3)$ converges to some $X_i(x_3)$. By Lemma \[lemsqrt\] (see below), for every $N\in \N$, $$\sum_{i\le N} {\varphi}_i^{\psi(n)}(x_3)\ge 1- \frac{f_{\psi(n)}(x_3)}{\sqrt{N}}.$$ This implies, by tightness, $\mu_{x_3}^{\psi(n)}\weaklim \sum_i {\varphi}_i(x_3) \delta_{X_i(x_3)}$ and $\sum_i ({\varphi}_i(x_3))^{1/2}\le g(x_3)$. But since $\mu_{x_3}^{\psi(n)}\weaklim {\mu}_{x_3}$, we have ${\mu}_{x_3}=\sum_i {\varphi}_i(x_3) \delta_{X_i(x_3)}$. Finally, by the subadditivity of the square root, and the definition of $g$ we obtain .

\[lemsqrt\] If a nonincreasing sequence of positive numbers $\gamma_i$ is such that $$\sum_i \gamma_i = c_0 \qquad \text{ and }\qquad \sum_i
\sqrt{\gamma_i} \le C_0,$$ then for all $N\in \N$ one has $$\sum_{i\le N} \gamma_i \ge c_0- C_0 \sqrt{\frac{c_0}{N}}.$$

Indeed $\sum_{i >N} \gamma_i \le \sqrt{\gamma_N} \sum_{i >N}
\sqrt{\gamma_i} \le C_0 \sqrt{\gamma_N},$ while $c_0 \ge \sum_{i\le
N} \gamma_i \ge N \gamma_N$.

Regularity of minimizers
------------------------

We now want to prove regularity of the minimizing measures $\mu$. In order to prove that we can restrict our attention to measures containing no loops, we first define the notion of subsystem.

\[Existence of subsystems\]\[def:subsystem\] Given a point $x:=(X, x_3) \in Q_{L,T}$ and $\mu\in{\mathcal{A}}^*_{L,T}$ with $I(\mu)<+\infty$, there exists a subsystem ${\widetilde}{{\mu}}$ of ${\mu}$ emanating from $x$, meaning that there exists ${\widetilde}{{\mu}} \in {\mathcal{A}}^*_{L,T}$ such that

1.  \[def:subsystempos\] $ {\widetilde}{{\mu}}\le {\mu}$ in the sense that ${\mu}-{\widetilde}{{\mu}}$ is a positive measure,

2.  \[def:subsystemdelta\] ${\widetilde}{{\mu}}_{x_3}=a\delta_{X}$, where $a={\mu}_{x_3}(X)$,

3.  \[def:subsystemconti\] $$\pa_3 {\widetilde}{{\mu}}+ \div' \lt(\frac{dm}{d{\mu}} {\widetilde}{{\mu}}\rt)=0.$$

In particular, \[def:subsystemdelta\] implies that $({\mu}_{x_3} - {\widetilde}{{\mu}}_{x_3}) \perp  \delta_{X}$ in the sense of the Radon-Nikodym decomposition.

Let us for notational simplicity assume that $x_3=0$, $L=T=1$, $\mu(Q_{1,1})=2$. Let us denote $v= \frac{d m}{d{\mu}}$. According to [@AGS Th. 8.2.1 and (8.2.8)], since $v\in L^2(Q_{1,1},\mu)$, there exists a positive measure $\sigma$ on $ C^0([-1, 1];Q_1)$ (endowed with the sup norm), whose disintegration [@AGS Th. 5.3.1] with respect to $\mu_0$, i.e. $\sigma=\int_{Q_1} \sigma_{x'} d\mu_0(x')$, is made of probability measures $\sigma_{x'}$ concentrated on the set of curves $ \gamma$ solving $$\left\{ \begin{array}{l}
\dot{\gamma}(x_3)= v(\gamma(x_3))\\
\gamma(0)= x',\end{array}\right.$$ and such that for every $x_3\in [-1,1]$, ${\mu}_{x_3}= (e_{x_3})_\# \sigma$, where $e_{x_3}$ denotes the evaluation at $x_3$, in the sense that $$\int_{Q_1}\varphi d\mu_{x_3} = \int_{  C^0([-1, 1];Q_1)} \varphi(\gamma(x_3)) d\sigma (\gamma)
 \hskip1cm \text{ for all }\varphi\in C^0(Q_1)\,.$$ Then, the measure ${\widetilde}{{\mu}}=   {\widetilde}{{\mu}}_{x_3}\otimes dx_3$ with ${\widetilde}{\mu}_{x_3}= (e_{x_3})_\# (a \sigma_{X})$, where $a=\mu_{0}(X)$, satisfies all the required properties.

\[noloop\] Let ${\mu}$ be a minimizer for the Dirichlet problem (\[limitProb\]), $\bar x_3\in (-T,T)$. Let $x_1=(X_1,\bar{x}_3)$, $x_2=(X_2,\bar{x}_3)$ be two points in the plane $\{x:x_3=\bar{x}_3\}$. Let ${\mu}_1$ and ${\mu}_2$ be subsystems of ${\mu}$ emanating from $x_1$ and $x_2$. Let $x_+:=(X_+,z_+)$ be a point with $ z_+  > \bar{x}_3$ and $x_-:=(X_-,z_-)$ a point with $ z_-<\bar{x}_3$, and such that ${\mu}_{1,z_+}, {\mu}_{2,z_+}$ both have Diracs at $X_+$, and ${\mu}_{1,z_-}, {\mu}_{2,z_-}$ both have Diracs at $X_-$ (with nonzero mass). Then $X_1=X_2$.

Let ${\varphi}_1:=\mu_{\bar x_3}(X_1)$ be the mass of ${\mu}_1$ and ${\varphi}_2$ be the mass of $\mu_2$. Let ${\varphi}_{1, +}:={\mu}_{1,\bar x_3}(X_+)$ be the mass of ${\mu}_1$ at $x_+$, ${\varphi}_{2, +}$ the mass of ${\mu}_2$ at $x_+$, ${\varphi}_{1, - } $ the mass of ${\mu}_1$ at $x_-$, ${\varphi}_{2,-}$ the mass of ${\mu}_2$ at $x_-$. Let ${\varphi}:= \min \{{\varphi}_{1, +},{\varphi}_{1, -}, {\varphi}_{2, +}, {\varphi}_{2,-}\}$ which by assumption is positive.

We define ${\mu}_{1, +}$ as the subsystem of ${\mu}_1$ coming from $x_+$, it is thus of mass ${\varphi}_{1, +}$, and at level $\bar{x}_3$ all its mass is at $X_1$ (since it is a subsystem of ${\mu}_1$ for which this is the case). Similarly with ${\mu}_{1, -}$, ${\mu}_{2, +}, {\mu}_{2,-}$. We can now define ${\widetilde}{{\mu}}_1:= \frac{{\varphi}}{{\varphi}_{1, +} } {\mu}_{1, +}$ for $x_3\ge\bar{x}_3$ and ${\widetilde}{{\mu}}_1:=\frac{{\varphi}}{{\varphi}_{1, -} }{\mu}_{1, -}$ for $x_3 <\bar{x}_3$, and the same with ${\widetilde}{{\mu}}_2$. The measures ${\widetilde}{{\mu}}_1$ and ${\widetilde}{{\mu}}_2$ are “systems” of mass ${\varphi}$ that join $X_-$ and $X_+$. By construction, we have $$\partial_{x_3}  ( {\widetilde}{{\mu}}_1- {\widetilde}{{\mu}}_2 ) \restr[z_-,z_+]+\div' \lt( ( {\widetilde}{{\mu}}_1- {\widetilde}{{\mu}}_2 ) \restr[z_-,z_+]\frac{dm}{d\mu}\rt)=0$$ and $$-\frac{\varphi}{\min(\varphi_{2,+}, \varphi_{2,-})}\mu \le  ( {\widetilde}{{\mu}}_1- {\widetilde}{{\mu}}_2 ) \restr[z_-,z_+]\le \frac{\varphi}{\min(\varphi_{1,+}, \varphi_{1,-})}\mu.$$

We now define $\hat{{\mu}}_\eta:=   {\mu}+ \eta ( {\widetilde}{{\mu}}_1- {\widetilde}{{\mu}}_2 ) \restr[z_-,z_+]$, which is admissible for $\eta$ small enough [(and different from $\mu$ unless $X_1=X_2$)]{}, and evaluate $$\begin{gathered}
I(\hat{{\mu}}_\eta)- I({\mu}) =
 K_*\int_{z_-}^{z_+}
 \sum_{x'\in Q_1}\left( \mu_{x_3} (x') +\eta ({\widetilde}{\mu}_1)_{x_3}(x')-\eta ({\widetilde}{\mu}_2)_{x_3}(x')  \right)^\hal
   -  \sum_{x'\in Q_1}\left( \mu_{x_3} (x') \right)^\hal \, dx_3 \\
+\eta  \int_{Q_1 \times [z_-, z_+] }
     \left(\frac{dm}{d{\mu}} \right)^2
     (   d{\widetilde}{{\mu}}_1- \, d{\widetilde}{{\mu}}_2)\,.
     \end{gathered}$$ But the function $\eta \mapsto  \sqrt{a+\eta b}$ is strictly concave for $a>0$ and $b\ne0$, therefore $I(\hat{{\mu}}_\eta)+I(\hat{{\mu}}_{-\eta})<2I({\mu})$ for any $\eta\ne 0$, a contradiction with the minimality of ${\mu}$.

A consequence of this lemma is the following. Consider a minimizing measure $\mu$ of $I$. Let $z_-$ and $z_+$ be any two slices and let $X_-$ be one of the Diracs at slice $z_-$. Let ${\widetilde}{{\mu}}$ be a subsystem emanating from $(X_-,z_-)$. Let $X_+$ be any point in the slice $z_+$ where ${\widetilde}{{\mu}}$ carries mass. Then, there is a unique “path" connecting $X_-$ to $X_+$ (otherwise there would be a loop). Since this is true for any couple of “sources" in two different planes, this means that there are at most a countable number of absolutely continuous curves (absolutely continuous because of the transport term) on which ${\mu}\restr[z_-, z_+]$ is concentrated. So we have a representation of the form $$\label{representationmu}
{\mu}= \sum_i \frac{\varphi_i}{\sqrt{1+|{\dot{X}}_i|^2}} \, \mathcal{H}^1\restr\Gamma_i   \, ,$$ where the sum is countable and $\Gamma_i=\{ (X_i(x_3),x_3) : x_3\in [a_i,b_i]\} $ with $X_i$ absolutely continuous and almost everywhere non overlapping.

Another consequence is that if there are two levels at which ${\mu}$ is a finite sum of Diracs, then it is the case for all the levels in between. So, if there is a slice with an infinite number of points, then either it is also the case for all the slices below or for all the slices above.

For measures which are concentrated on finitely many curves we obtain a simple representation formula for $I(\mu)$.

\[lemmacurves\] Let $\mu=\sum_{i=1}^N \frac{\varphi_i}{\sqrt{1+|{\dot{X}}_i|^2}} \, \mathcal{H}^1\restr \Gamma_i \in{\mathcal{A}}^*_{L,T}$ with $\Gamma_i=\{ (X_i(x_3),x_3) : x_3\in [a_i,b_i]\}$ for some absolutely continuous curves $X_i$, almost everywhere non overlapping. Every ${\varphi}_i$ is then constant on $[a_i,b_i]$ and we have conservation of mass. That is, for $x:=(x',x_3)$, letting $$\begin{aligned}
{\mathcal{I}}^+(x)&:=\{ i\in[1,N] \ : \ x_3=b_i, \ X_i(b_i)=x'\},  \\ {\mathcal{I}}^-(x)&:=\{ i\in[1,N] \ : \ x_3=a_i, \ X_i(a_i)=x'\},  \end{aligned}$$ it holds $$\sum_{i\in{\mathcal{I}}^-(x)} {\varphi}_i=\sum_{i\in{\mathcal{I}}^+(x)} {\varphi}_i.$$ Moreover, $m=\sum_i \frac{\varphi_i}{\sqrt{1+|{\dot{X}}_i|^2}}
{\dot{X}}_i \, \mathcal{H}^1\restr\Gamma_i$ and $$\label{Iparticul}
 I(\mu)=\frac{1}{L^2}\sum_i \int_{a_i}^{b_i} K_* \sqrt{{\varphi}_i}+ {\varphi}_i |{\dot{X}}_i|^2 dx_3.$$

Let $\bar x=(\bar x',\bar x_3)$ with $\bar x_3 \in (-T,T)$ be such that $\mu_{\bar x_3}(\bar x')\neq 0$. Then, by continuity of the $X_i$’s, there exist $\delta>0, \eps>0$ such that every curve $\Gamma_i$ with $\Gamma_i\cap\lt( \B_\eps'(\bar x')\times [\bar x_3-\delta,\bar x_3+\delta]\rt) \neq \emptyset$ satisfies $X_i(\bar x_3)=\bar x',$ and such that $\mu\restr (\B'_{2\eps}(\bar x')\backslash \B_\eps'(\bar x'))\times [\bar x_3-\delta,\bar x_3+\delta]=0$ (and thus also $m\restr (\B'_{2\eps}(\bar x')\backslash \B_\eps'(\bar x'))\times [\bar x_3-\delta,\bar x_3+\delta]=0$ since $m\ll {\mu}$). Consider then $\psi_1\in C^\infty_c(\B'_{2\eps}(\bar x'))$ with $\psi_1=1$ in $\B_\eps(\bar x')$ and $\psi_2\in C^\infty_c( \bar x_3-\delta,\bar x_3+\delta)$ and test with $\psi:=\psi_1 \psi_2$ to obtain $$\begin{aligned}
\int_{\bar x_3-\delta}^{\bar x_3+\delta} \frac{d\psi_2}{dx_3} (x_3)\lt(\sum_{X_i(x_3)\in \B_\eps'(\bar x')} {\varphi}_i\rt) dx_3&=\int_{\bar x_3-\delta}^{\bar x_3+\delta}\int_{\B'_\eps(\bar x')} \frac{d\psi_2}{dx_3}(x_3)\psi_1(x') d\mu\\
&=\int_{Q_{1,1}} \frac{\partial \psi}{\partial x_3} d\mu=-\int_{Q_{1,1}} \nabla'\psi \cdot dm\\
 &=-\int_{\bar x_3-\delta}^{\bar x_3+\delta}\int_{\B_{2\eps}(\bar x')\backslash \B_\eps(x')} \psi_2 \nabla' \psi_1 \cdot dm =0,\end{aligned}$$ from which the first two assertions follow. It can be easily checked that this implies that $\bar m:= \sum_i \frac{\varphi_i}{\sqrt{1+|{\dot{X}}_i|^2}}
{\dot{X}}_i \, \mathcal{H}^1\restr\Gamma_i$ satisfies . Let $m$ be any other measure satisfying and let us prove that $\nu:=\bar m -m=0$. Since $\Divp \nu_{x_3}=0$, we have for every $\psi\in C^{\infty}(Q_1)$ $$\sum_i \nabla' \psi(X_i(x_3))\cdot \nu_i(x_3)=0,$$ where $\nu_{x_3}=\sum_i \nu_i (x_3) \delta_{X_i(x_3)}$, from which the claim follows.

We remark that Corollary \[coroXia\] below will imply that representation holds for every measure $\mu$ with $I(\mu)<+\infty$. The previous results lead to the following.

\[reg\] A minimizer of the Dirichlet problem (\[limitProb\]) with boundary conditions $\bar{\mu}_+= \sum_{i =1}^N {\varphi}_i^+\delta_{X_i^+}$ and $\bar{\mu}_{-} =
\sum_{i=1}^N {\varphi}_i^- \delta_{X_i^-} $ (some ${\varphi}_i$ may be zero) satisfies

1.  $\mu= 
    \sum_{i=1}^M \frac{\varphi_i}{\sqrt{1+|{\dot{X}}_i|^2}} \, \mathcal{H}^1\restr\Gamma_i  $ for some $M\in\N$, where $\Gamma_i=\{ (X_i(x_3),x_3) : x_3\in [a_i,b_i]\}$ are disjoint up to the endpoints, and the $X_i$ are absolutely continuous.

2.  Each $X_i$ is affine.

3.  If $\bar{\mu}_{-}= \bar{\mu}_{+} $ then there exists a symmetric minimizer with respect to the $x_3=0$ plane.

Let $\mu^{ij}$ be the subsystem emanating from $(X_i^-,-T)$ of the subsystem emanating from $(X_j^+,T)$ of $\mu$, so that $\mu=\sum_{ij}\mu^{ij}$ by $\mu^{ij}_{x_3}\le \mu_{x_3}$ and conservation of mass. By Lemma \[noloop\] we have $\mu^{ij}_{x_3}={\varphi}_{ij}(x_3)\delta_{X^{ij}}(x_3)$ for all $x_3$, otherwise there would be loops. By Lemma \[lemmacurves\], ${\varphi}_{ij}(x_3)$ does not depend on $x_3$. By (\[HolderW2\]), if ${\varphi}_{ij}>0$ then $X^{ij}$ is absolutely continuous. After a relabeling, (i) is proven.

Assertion (ii) follows from minimizing $I({\mu})$ as given by with respect to ${\dot{X}}_i$.

Let now $\bar{\mu}_{-}= \bar{\mu}_{+}$. If $I(\mu,(-T,0))\le I(\mu,(0,T))$ we obtain a symmetric minimizer $\hat\mu$ by reflection of $\mu\LL(-T,0)$ across $\{x_3=0\}$, and analogously in the other case. This proves (iii).

We now show that for symmetric minimizers, at arbitrarily small distance from the boundary we have a finite number of Diracs. We already know that at arbitrarily small distance we have a countable number, and then that we have a representation of ${\mu}$ of the form . Let us point out that we will not use this proposition but rather include it for its own interest.

\[finitebranch\] Fix $\bar\mu\in {\mathcal{M}}^+(Q_L)$. Let $\mu$ be a symmetric minimizer of $I$ subject to $\mu_{\pm T}=\bar\mu$. Then for any $\delta>0$ sufficiently small, the number of Diracs in each slice $x_3 \in [-T +
\delta T , T- \delta T]$ is ${\stackrel{<}{\sim}}\delta^{-4}$.

We may assume $L=T=1$, $\mu(Q_{1,1})=2$. By symmetry, we need only to consider the interval $[0,1-\delta]$. If $\mu_{1-\delta}=\sum_{i} \p_i \delta_{X_i}$, it suffices to prove that ${\varphi}_i{\stackrel{>}{\sim}}\delta^{4}$ for every $i$. For the rest of the proof we fix a point $X_i$ and in order to ease notation we write ${\varphi}:={\varphi}_i$ and $X:= X_i$. Let ${\widetilde}{{\mu}}$ be the subsystem emanating from $(X,1-\delta)$. Thanks to the symmetry of ${\mu}$ and to the no-loop condition, $\mu$ and $\mu-{\widetilde}{{\mu}}$ are disjoint for $x_3> 1-\delta$. Indeed, if this was not the case, by symmetry they would meet also for $x_3<- 1+\delta$, and there would be a loop, which is excluded by Lemma \[noloop\]. Therefore $$\label{secondestimI}
  I({\mu})- I({\mu}- {\widetilde}{{\mu}})\ge \int_{1-\delta}^1 \sum_{x'\in Q_1} \left( {\widetilde}{{\mu}}_{x_3}(x')\right)^\hal  \, dx_3\ge \delta \sqrt{\p},$$ where in the second step we used subadditivity of the square root.

Let now ${\widetilde}{\mu}_{1}$ be the trace of ${\widetilde}{{\mu}}$ on $x_3={1}$ and for $z:={\varphi}^{1/4}$, let $\hat\mu$ be the symmetric comparison measure constructed as follows:

-   in $[0, 1-z]$, $\hat\mu_{x_3}:= (1+\frac{\p}{1-\p})(\mu_{x_3}-{\widetilde}{\mu}_{x_3})$ and

-   in $[1-z,1]$, $\hat\mu_{x_3}:=\mu_{x_3}-{\widetilde}{\mu}_{x_3}+\nu_{x_3}$ where $\nu$ is a measure connecting $\frac{\p}{1-\p} (\mu-{\widetilde}{{\mu}})_{1-z}$ to ${\widetilde}{{\mu}}_{1}$ constructed in Proposition \[branchingmu\], so that (recall ) $$ I^{(1-z,1)}(\nu)\les z\sqrt{\p}+\frac{\p}{z}\sim \p^{3/4}.$$

Since ${\mu}$ is a minimizer it follows by subadditivity of the energy that, for some universal (but generic) constant $C$, $$\begin{aligned}
  I(\mu)\le I(\hat\mu)&\le \lt(1+\frac{{\varphi}}{1-\p}\rt) I(\mu-{\widetilde}{{\mu}})+ C\p^{3/4}\\
 &\le I(\mu-{\widetilde}{{\mu}})+C\p^{3/4}.
 \end{aligned}$$ Indeed, $I(\mu-{\widetilde}{{\mu}}) \le I(\mu) {\stackrel{<}{\sim}}1$ while $\varphi \ll 1$ without loss of generality. Recalling , we deduce that $C\p^{3/4}\ge \delta \p^{1/2}$, which yields the result.

\[defpolygonal\] We say that a measure $\mu$ is polygonal if $$\mu= \sum_i \frac{\varphi_i}{\sqrt{1+|{\dot{X}}_i|^2}} \, \mathcal{H}^1\restr\Gamma_i,$$ where the sum is countable, $\Gamma_i$ are segments of the form $\Gamma_i=\{ (X_i(x_3),x_3) : x_3\in [a_i,b_i]\}$ disjoint up to the endpoints, and for any $z\in (0,T)$ only finitely many segments intersect $Q_L\times (-z,z)$. We say it is finite polygonal if the total number of segments is finite.

For any polygonal measure, the representation formula holds.

Let $\bar {\mu}\in {\mathcal{M}}^+(Q_L)$. Then, every symmetric minimizer of $I$ with boundary data $\mu_{\pm T}=\bar {\mu}$ is polygonal.

It suffices to show that for any $z\in (0,T)$ the measure $\mu$ is polygonal in $Q_L\times (-z,z)$. By Proposition \[finitebranch\] the measures $\mu_{z}$ and $\mu_{-z}$ are finite sums of Diracs. Since $\mu$ is a minimizer, it minimizes $I$ restricted to $(-z,z)$ with boundary data $\mu_{\pm z}$. By Proposition \[reg\] we conclude.

Density of regular and quantized measures
-----------------------------------------

In this section we want to prove that when $\bar \mu_{\pm T }=  \Phi L^{-2} dx'$, the set of “regular" measures is dense in energy.

\[defregular\] We denote by $\Mreg(Q_{L,T}){\subseteq}{\mathcal{A}}^*_{L,T}$ the set of regular measures, i.e., of measures $\mu$ such that:

-   The measure $\mu$ is finite polygonal, according to Def. \[defpolygonal\].

-   All branching points are triple points. This means that any $x\in Q_{L,T}$ belongs to the closures of no more than three segments.

For $N\in \N$, we say that $\mu$ is $N$-regular, $\mu\in \MNreg(Q_{L,T}){\subseteq}\Mreg(Q_{L,T})$, if in addition

-   The traces obey $\mu_{T}=\mu_{-T}=\Phi N^{-2}\sum_{j} \delta_{X_j}$, where the $X_j$ are $N^2$ points on a square grid, spaced by $L/N$, and $\Phi\ge0$.

We can now state the main theorem of this section:

\[theodens\] For every measure $\mu\in {\mathcal{A}}^*_{L,T}$ with $I(\mu)<+\infty$ and $ \mu_{\pm T}= \Phi L^{-2} dx'$, there exists a sequence of measures $\mu_N
  \in\MNreg$, with $\mu_N \weaklim \mu$ and such that $\limsup_{N \to +\infty} I(\mu_N)\le I(\mu)$, $\mu_N(Q_{L,T})=\mu(Q_{L,T})$.

The proof will be based on the following intermediate result.

\[lemmamudensityinside\] For every measure $\mu\in {\mathcal{A}}^*_{L,T}$ with $I(\mu)<+\infty$, and such that the traces $\mu_T$ and $\mu_{-T}$ are finite sums of Diracs, there exists a sequence of regular measures $\mu^{(N)}\in \Mreg(Q_{L,T})$ with $\mu^{(N)} \weaklim \mu$, $\mu^{(N)}_{\pm T}=\mu_{\pm T}$, and such that $\limsup_{N \to +\infty} I(\mu^{(N)})\le I(\mu)$.

We shall modify $\mu$ in two steps to make it polygonal: first on finitely many layers, to have finitely many Diracs on each of them, and then in the rest of the volume, using local minimization.

Fix $N\in\N$ and $\delta>0$, both chosen later. We choose levels $z_j\in (Tj/N, T(j+1)/N)$, for $j=-N+1,\dots, N-1$, with the property that $\mu_{z_j}=\sum_{k\in\N} \varphi_{j,k}\delta_{x'_{j,k}}$, with $\sum_{k\in \N} \varphi_{j,k}^{1/2}<+\infty$ for every $j$. We shall iteratively truncate the measure at these levels so that it is supported on finitely many points; for notational simplicity we also define $z_{\pm N}:=\pm T$. The measures $\mu^j$, $j=-N,\dots, N$ will all satisfy $\mu^j\ll\mu$ and $I(\mu^j)\le I(\mu)$.

We start with $\mu^{-N}:=\mu$. In order to construct $\mu^{j+1}$ from $\mu^j$, we first choose $K_j$ such that $\sum_{k\ge K_j} \varphi_{j,k}^{1/2}\le \delta/N$. Then we define $\mu^{j+1}$ as the sum of the subsystems of $\mu^j$ originating from the points $(x'_{j,k},z_j)$ with $k<K_j$. Clearly $I(\mu^{j+1})\le I(\mu^j)\le I(\mu)$. At the same time, since $\sum_{k=K_j}^{+\infty} {\varphi}_{j,k} \le \Phi^{1/2} \sum_{k=K_j}^{+\infty} {\varphi}_{j,k}^{1/2}$, $$\mu^{j+1}_{z_j}(Q_L)
=  \mu^{j}_{z_j}(Q_L)- \sum_{k=K_j}^{+\infty} {\varphi}_{j,k}\ge \mu^{j}_{z_j}(Q_L)-\frac{ \Phi^{1/2} \delta}{N},$$ so that $|\mu^{j+1}-\mu^j|(Q_{L,T})\le 2\Phi^{1/2}\delta T/N$. Therefore $|\mu^N-\mu|(Q_{L,T})\les \delta T \Phi^{1/2}$.

We define $\hat\mu^N$ as the minimizer with boundary data $\mu^N_{z_j}$ and $\mu^N_{z_{j+1}}$ in each stripe $Q_L\times (z_j, z_{j+1})$. Then $I(\hat\mu^N)\le I(\mu^N)\le I(\mu)$, and $\hat\mu^N_{\pm T}=\mu^N_{\pm T}$.

At this point we fix the boundary data. For this, we let ${\widetilde}{\mu}$ be the minimizer on $Q_{L,T}$ with boundary data ${\widetilde}{\mu}_\pm:=\mu_{\pm T}-\mu^N_{\pm T}$. These boundary data are finite sums of Diracs, and their flux is $|{\widetilde}{\mu}_\pm|(Q_L)=|\mu_T-\mu^N_T|(Q_L)\les \delta \Phi^{1/2}$. By Proposition \[reg\] the minimizer is finite polygonal, by Proposition \[branchingmu\] it has energy no larger than a constant times $\delta \Phi^{1/2}T^{-1} + \delta^{1/2} T \Phi^{1/4} L^{-2}$. Finally we set $\mu^{(N)}:=\hat\mu^N+{\widetilde}{\mu}$. Then $\mu^{(N)}_{\pm T}=\mu_{\pm T}$, and $$I(\mu^{(N)})\le I(\hat\mu^N)+I({\widetilde}{\mu}) \le I(\mu) +  C\left(\delta \Phi^{1/2}T^{-1} + \delta^{1/2} T \Phi^{1/4} L^{-2}\rt) .$$

Up to a small perturbation, we may further assume that all junctions are triple. We can now choose for instance $\delta=1/N$. It only remains to show that $\mu^{(N)}\weaklim\mu$ as $N\to+\infty$. Recalling that $|\mu^N-\mu|(Q_{L,T})+|{\widetilde}{\mu}|(Q_{L,T})\les \delta \Phi^{1/2}T$, we only need to show that $\hat\mu^N-\mu^N\weaklim 0$.

For $x_3\in (z_j,z_{j+1})$ we have by definition of $\hat \mu^N$ $$W_2(\mu_{x_3}^N, \hat \mu_{x_3}^N) \le W_2(\mu_{x_3}^N, \mu_{z_j}^N) +W_2(\hat \mu_{x_3}^N, \hat\mu_{z_j}^N),$$ and by $$W_2^2(\mu_{x_3}^N, \mu_{z_j}^N)+W_2^2(\hat \mu_{x_3}^N, \hat \mu_{z_j}^N) \le L^2 (z_{j+1}-z_j) \lt( I(\mu^N)+I(\hat \mu^N)\rt),$$ so that $$\max_{x_3} W_2^2(\mu_{x_3}^N, \hat \mu_{x_3}^N) {\stackrel{<}{\sim}}\frac{L^2T}{N}I(\mu).$$ For $x_3\in (-T,T)$, let $\Pi_{x_3}$ be an optimal transport plan from $\mu^N_{x_3}$ to $\hat{\mu}^N_{x_3}$. Considering then the transport plan $\Pi:=\Pi_{x_3}\otimes dx_3$ between $\mu^N$ and $\hat{\mu}^N$, we get $$W_2^2(\mu^N,\hat{\mu}^N)\le \int_{-T}^T W^2_2(\mu_{x_3}^N,\hat\mu_{x_3}^N) dx_3\les \frac{L^2T^2}{N}I(\mu)
\,,$$ which yields that indeed $\hat\mu^N-\mu^N\weaklim 0$.

Let $\eps\in (0,1/4)$, chosen such that it tends to zero as $N\to\infty$. We define $\hat\mu$ in $Q_{L,(1-2\eps)T}$ as a rescaling by $(1-2\eps)$ in the vertical direction, $\hat\mu_{(1-2\eps)x_3}=\mu_{x_3}$. An easy computation shows that $I(\hat\mu, -(1-2\eps)T, (1-2\eps)T)\le \frac{1}{1-2\eps} I(\mu)$. In particular, $\hat\mu_{(1-2\eps)T}=\Phi L^{-2} dx'$. For $x_3\in ((1-2\eps)T, T)$ we define $\hat\mu$ as the result of Proposition \[branchingmu\]. Then we set ${\widetilde}{\mu}=\hat\mu$ on $(-(1-\eps)T,(1-\eps)T)$ and extend it constant outside, in the sense that ${\widetilde}{\mu}_{x_3}=\hat\mu_{(1-\eps)T}$ for $x_3\in ((1-\eps)T,T)$, and the same on the other side. Since $\hat \mu_{(1-\eps)T}$ is the midplane configuration of the branching measure constructed in Proposition \[branchingmu\], ${\widetilde}{\mu}_{x_3}$ is a finite sum of Diracs for $|x_3|\ge (1-\eps) T$. We obtain $$I({\widetilde}{\mu}) \le \frac{1}{1-2\eps} I(\mu) + C \eps^{1/3} \left(\Phi^{2/3}T^{1/3} L^{-4/3}+\eps^{2/3} \Phi^{1/2} T  L^{-2} \right)\,.$$ By Lemma \[lemmamudensityinside\] applied to the inner domain $(-(1-\eps)T,(1-\eps)T)$ there is a finite polygonal measure $\check\mu$ which is close to ${\widetilde}{\mu}$ and has the same boundary data at $x_3=\pm T(1-\eps)$. The measure given by $\check{\mu}$ inside, and ${\widetilde}{\mu}$ outside, has the required properties.

We now turn to the quantization of the measures.

We say that a regular measure $\mu\in\Mreg(Q_{L,T}){\subseteq}{\mathcal{A}}_{L,T}^*$ is $k$-quantized, for $k>0$, if for all $(x',z)\in Q_{L,T}$ one has $k\mu_z(\{x'\})\in 2\pi\N$.

\[lemmaquantize\] Let $\mu\in \MNreg(Q_{L,T})$ and $\Phi=\mu_T(Q_L)$. For any $k>0$ such that $k\Phi\in 2\pi\N$ there is a $k$-quantized regular measure $\mu^{k}\in \Mreg$ such that $\mu^{k}(Q_{L,T})=\mu(Q_{L,T})$ and $$\left(1-\frac{ C(\mu)}k\right) \mu \le \mu^k \le  \left(1+\frac{ C(\mu)}k\right)\mu\,.$$ This implies in particular $\mu^k\ll\mu$, $\mu^k\to\mu$ strongly, $W_2^2(\mu,\mu_k)\les C(\mu)k^{-1}$ and $I(\mu^k)\to I(\mu)$ as $k\to\infty$.

The measure $\mu$ consists of finitely many segments, each with a flux. To prove the assertion it suffices to round up or down the fluxes to integer multiples of $2\pi/k$ without breaking the divergence condition, and without changing the total flux.

Since $\mu\in\MNreg(Q_{L,T})$, we have $\mu_{T}=\Phi N^{-2}\sum_i \delta_{X_i}$. We select $\varphi_i^k$ as $2\pi \lfloor k \Phi /(2\pi N^2)\rfloor /k$ or $2\pi \lfloor k \Phi /(2\pi N^2)+1\rfloor /k$ , depending on $i$. Precisely, we choose the first value for $i=0$ and then, at each $i$, we choose the lower one if $\sum_{j<i} ({\varphi}_j^k- \Phi N^{-2})>0$, and the upper one otherwise. This concludes the definition of $\mu^k_T$.

The fluxes in the interior of the sample are defined by propagating the rounding. At each point where a bifurcation occurs, if there is more then one outgoing branch we distribute the rounding as discussed for $\mu_T^k$. This increases the maximal error by at most $2\pi/k$, at each branching point. Since $\mu$ is finite polygonal, there is a finite number of branching points, hence the total error is bounded by a constant times $1/k$. Precisely $|{\varphi}_i^k-{\varphi}_i|\le C(\mu)/k$ for any segment $i$. Since ${\varphi}_i$ only takes finitely many values, $|{\varphi}_i^k-{\varphi}_i|\le {\varphi}_i C(\mu)/k$ for any segment $i$, which concludes the proof.

Relation with irrigation problems {#Sec:reliri}
---------------------------------

The functional $I(\mu)$ bears similarities with the so-called irrigation problems which have attracted a lot of interest (see for instance [@Xia; @BeCaMo]). Besides their applications to the modeling of communication networks and other branched patterns (see again [@BeCaMo] and the references therein), they have also been recently used in the study of Sobolev spaces between manifolds [@bethuel]. Let us recall their definition and for this, follow the notation of [@BeCaMo]. For $E(G)$ a set of oriented straight edges and $\p: E(G)\to (0,+\infty)$ we define the irrigation graph $G$ as the vector measure $$G:=\sum_{e\in E(G)} {\varphi}(e) \mathbf{e} \, \H^1\restr e \,$$ where $\mathbf{e}$ is the unit tangent vector to $e$. For $\a\in [0,1]$, we then define the Gilbert energy of $G$ by $$M^{\a}(G):=\sum_{e\in E(G)} \p(e)^\a \H^1(e).$$ Given two atomic probability measures $\mu^+=\sum_{i=1}^k a_i \delta_{X_i}$ and $\mu^-=\sum_{j=1}^l b_j \delta_{Y_j}$, we say that $G$ irrigates $(\mu^+,\mu^-)$ if $\div G= \mu^+-\mu^-$ in the sense of distributions (this implies in particular that $G$ satisfies Kirchoff’s law). If we are now given any two probability measures $(\mu^+,\mu^-)$ and a vector measure $G$, with $\div G=\mu^+-\mu^-$ (sometimes called an irrigation path between $\mu^+$ and $\mu^-$), we define $$M^{\a}(G):=\inf \{ \liminf_{i\to +\infty} M^{\a}(G_i)\},$$ where the infimum is taken among all the sequences of irrigation graphs $G_i$ with $G_i\weaklim G$ in the sense of measures and such that $\div G= \mu^+_i-\mu^-_i$ for some atomic measures $\mu^{\pm}_i$ tending to $\mu^{\pm}$. If no such sequence exists then we set $M^\a(G)=+\infty$. The irrigation problem then consists in minimizing $M^{\a}(G)$ among all the transport paths $G$ between $\mu^+$ and $\mu^-$. For $\a=0$ this is a generalization of the famous Steiner problem while for $\a=1$ it is just the Monge-Kantorovich problem.

Using some powerful rectifiability criterion of B. White, the following theorem was proven by Q. Xia [@Xia].

\[theoXia\] Given $0<\a<1$, any transport path $G$ with $M^{\a}(G)+M^1(G)<+\infty$ is rectifiable in the sense that $$G={\varphi}\tau \, \H^1\restr\Gamma$$ for some density function $\p$ and some $1-$rectifiable set $\Gamma$ having $\tau$ as tangent vector.

For minimal irrigation paths, much more is known about their interior and boundary regularity [@BeCaMo]. For instance, as for our functional $I(\mu)$ (see Proposition \[reg\]), it can also be proven that minimal irrigation paths contain no loops and that for $\a> 1-\frac{1}{n}$ (where $n$ is the dimension of the ambient space i.e. $n=3$ for us), any two probability measures $\mu^{\pm}$ can be irrigated at a finite cost (compare with Proposition \[existmu\]).

Using Theorem \[theoXia\] and Lemma \[lemmamudensityinside\], we can obtain the following rectifiability result.

\[coroXia\] Every measure $\mu$ for which $I(\mu)<+\infty$ is rectifiable.

Using the construction of Lemma \[lemmamudensityinside\], we can find a sequence $\mu^n$ such that $\mu^n\weaklim \mu$, $\limsup_{n\to +\infty} I(\mu^n)\le I(\mu)$ and $\mu^n=\sum_{i=1}^N \frac{{\varphi}_i}{\sqrt{1+|{\dot{X}}_i|^2}} \H^1\restr \Gamma_i$ for some straight edges $\Gamma_i=\{(x_3,X_i(x_3)) \, :\, x_3\in(a_i,b_i)\}$. Letting $\widetilde{\mu}^n:=\sum_{i=1}^N \frac{{\varphi}_i}{\sqrt{1+|{\dot{X}}_i|^2}} \begin{pmatrix} {\dot{X}}_i\\ 1\end{pmatrix} \H^1\restr \Gamma_i$, we have for $\a\ge\frac{3}{4}$, $$\begin{aligned}
M^\a({\widetilde}\mu^n)&=\sum_i \int_{a_i}^{b_i} {\varphi}_i^\alpha \sqrt{1+ |{\dot{X}}_i|^2} dx_3\\
&\les \sum_i \int_{a_i}^{b_i} {\varphi}_i (1+|{\dot{X}}_i|^2) + {\varphi}_i^{2\alpha-1} dx_3\\
&\les \sum_i \int_{a_i}^{b_i} {\varphi}_i (1+|{\dot{X}}_i|^2) + \sqrt{{\varphi}_i} dx_3\les I(\mu^n) +1\end{aligned}$$ so that $\liminf_{n\to +\infty} M^{\a}(\widetilde{\mu}^n)\les I(\mu)+1<+\infty$ and by Theorem \[theoXia\], the claim follows.

In [@OudSant], an approximation of the functional $M^\a$ in the spirit of the Modica-Mortola [@ModMort] approximation of the perimeter was proposed. Even though their proofs and constructions are completely different from ours, this approach bears some similarities with our derivation of the functional $I({\mu})$ from the Ginzburg-Landau functional $E_T(u,A)$.

Lower bound {#Sec:gammaliminf}
===========

In the rest of the paper we consider sequences with $$\label{eqlbassumptcoeff}
T_n\to +\infty, \hskip5mm \alpha_n\to +\infty,\hskip5mm  \beta_n\to 0, \hskip5mm \frac{T_n}{\alpha_n} \to +\infty,\hskip5mm \alpha_n \beta_n^{7/2}\to +\infty\,.$$ No constant appearing in the sequel will depend on the specific choice of the sequence. We observe that (\[eqlbassumptcoeff\]) immediately implies $\alpha_n\beta_n^2\to+\infty$ and $\alpha_n^2\beta_n\to+\infty$. Let us recall that in this proof we set ${\widetilde}L=1$ and that (see (\[eqdefwe\])) $$\begin{gathered}
  \wE(u,A)=\int_{Q_{1,1}} \a^{-2/3}\b^{-1/3}\lt|\nabla_{\a^{1/3}\b^{-1/3}TA}'u\rt|^2+ 
 \a^{-4/3}\b^{-2/3}\lt|(\nabla_{\a^{1/3}\b^{-1/3}TA} u)_3\rt|^2\\
+  \a^{2/3}\b^{-2/3}\left({B_3} -(1-|{u}|^2)\right)^2
+  \b^{-1}|{B'}|^2 dx +\a^{1/3}\b^{7/6}\|\b^{-1}{B_3}-1\|_{H^{-1/2}(x_3=\pm1)}^2.\end{gathered}$$ In this section, we prove the following compactness and lower bound result.

\[gammaliminf\] Fix sequences of positive numbers $\alpha_n$, $\beta_n$, $T_n$ such that (\[eqlbassumptcoeff\]) holds, and let $(u_n, A_n)$ be such that $\sup_{n} \wE(u_n,A_n)<+\infty$. Then up to a subsequence, the following holds :

1.  $\b_n^{-1}(1-\rho_n)\weaklim \mu$ for some measure $\mu$, $\b_n^{-1} B'_n\weaklim m$ for some vector-valued measure $m\ll {\mu}$ satisfying the continuity equation .

2.  For almost every $x_3\in (-1,1)$, there exists some probability measure $\mu_{x_3}$ on $Q_1$ with $\mu=\mu_{x_3}\otimes dx_3$ and such that $\mu_{x_3}\weaklim dx'$ as $x_3\to \pm 1$.

3.  For almost every $x_3\in (-1,1)$, ${\mu}_{x_3}=\sum_{i\in {\mathcal{I}}} \p_i \delta_{X_i}$ with ${\mathcal{I}}$ at most countable and $\p_i>0$.

4.  One has $(\mu,m)\in {\mathcal{A}}_{1,1}$ with $$\liminf_{n\to +\infty} \wE(u_n,A_n)\ge I({\mu},m).$$

Let us first show that the energy gives a quantitative control on the failure of the Meissner condition $\rho B=0$ in a weak sense.

\[lemmameissner\] For every $Q_1$-periodic test function $\psi\in H^1_\mathrm{per}(Q_{1,1})$, if $\|\rho\|_\infty\le 1$ then $$\label{estimmeissner} 
\lt|\int_{Q_{1,1}} \rho B_3 \psi dx\rt|\les 
\frac{\a^{1/3}\b^{2/3}}{T}\wE(u,A)\|\psi\|_{L^\infty}+
\frac{\beta^{1/2}}{T}\wE(u,A)^\hal   \|{\nabla}'
\psi\|_{L^2},$$ and, if additionally $\psi(x',\pm1)=0$, for $k=1,2$ and $\alpha^2\beta\ge 1$, $$\label{estimmeissnerprime} 
\lt|\int_{Q_{1,1}} \rho B_k \psi dx\rt|\les 
\frac{\a^{2/3}\b^{5/6}}{T}\wE(u,A)\|\psi\|_{L^\infty}+
\frac{\alpha^{1/3}\beta^{2/3}}{T}\wE(u,A)^\hal   \|{\nabla}\psi\|_{L^2}.$$ Moreover, if $\xi\in H^1_0(-1,1)$ and $\psi$ is a periodic Lipschitz continuous function on $Q_1$ then $$\label{estimmeissnerter} 
\lt|\int_{Q_{1,1}} \rho B'\cdot \nabla' \psi \, \xi dx\rt|\les 
\frac{\a^{2/3}\b^{5/6}}{T}\wE(u,A)\|\xi \nabla'\psi\|_{L^\infty}+
\frac{\b^{1/2}}{T}\wE(u,A)^\hal   \| \partial_3 \xi \nabla'\psi \|_{L^2}.$$

Let $\lambda:=\a^{1/3}\b^{-1/3}T$. For (\[estimmeissner\]) we use formula with $A$ substituted by $\lambda A$, that is $|{\nabla}'_{\lambda A} u|^2 = |{\mathcal{D}}_{\lambda A}^3 u|^2 +\rho \lambda B_3+
{\nabla}'\times j'_{\lambda A}$. We integrate against a test function $\psi$, $$\begin{aligned}
\lt|\int_{Q_{1,1}} \rho B_3 \psi dx\rt| & = & \frac{\a^{1/3}\b^{2/3}}{T}\lt|\int_{Q_{1,1}}\a^{-2/3}\b^{-1/3} \left(|{\nabla}_{\lambda A}'u|^2 - |{\mathcal{D}}_{\lambda A}^3 u|^2 - {\nabla}'\times j'_{\lambda A}\right) \psi dx\rt|\\
& \les &  \frac{\a^{1/3}\b^{2/3}}{T}\lt(\wE(u,A)\|\psi\|_{L^\infty} + \a^{-1/3}\b^{-1/6}\int_{Q_{1,1}} \a^{-1/3}\b^{-1/6}|j'_{\lambda A}| |{\nabla}' \psi|dx\rt)\\
& \les & \frac{\a^{1/3}\b^{2/3}}{T}\lt(\wE(u,A)\|\psi\|_{L^\infty}+ \a^{-1/3}\b^{-1/6}\wE(u,A)^\hal   \|{\nabla}'
\psi\|_{L^2} \rt),\end{aligned}$$ [where we have used that $|j'_{\lambda A}|\le |\nabla'_{\lambda A} u| $ in view of the definition and the upper bound $\rho\le 1$.]{} We obtain similarly : One first checks from the definition of ${\mathcal{D}}_{\lambda A}u$ that $$\left| |{\mathcal{D}}^1_{\lambda A} u|^2- |(\nabla_{\lambda A} u)_3|^2 - |(\nabla_{\lambda A} u)_2|^2\right|
 \le 2 |(\nabla_{\lambda A} u)_2| \, |(\nabla_{\lambda A} u)_3|.$$ Testing (\[magic2\]) with $\psi$ and integrating by parts the term with $j_{\lambda A}$ as above gives $$\lambda \lt|\int_{Q_{1,1}} \rho B_k \psi dx\rt|\les 
\int_{Q_{1,1}} 2 |\psi|  \, |(\nabla_{\lambda A} u)_2| \, |(\nabla_{\lambda A} u)_3| 
+ |\nabla \psi| \, |\nabla_{\lambda A} u| dx\,.$$ Estimating $ 2|(\nabla_{\lambda A} u)_2| \, |(\nabla_{\lambda A} u)_3| \le  \alpha^{1/3}\beta^{1/6}|(\nabla_{\lambda A} u)_2|^2 +   \alpha^{-1/3}\beta^{-1/6}|(\nabla_{\lambda A} u)_3|^2 $ and $\|\nabla_{\lambda A} u\|_2^2\le \alpha^{4/3}\beta^{2/3} \wE(u,A)$ concludes the proof of (\[estimmeissnerprime\]).\
The proof of is very similar to the proof of . Arguing as above with $\xi\partial_{k}\psi$ playing the role of $\psi$, we get $$\lambda \lt|\int_{Q_{1,1}} \rho B'\cdot \nabla'\psi \xi dx\rt|\les 
\int_{Q_{1,1}} 2 |\xi \nabla' \psi |  \, |(\nabla_{\lambda A} u)_2| \, |(\nabla_{\lambda A} u)_3|  dx
+ \lt|\int_{Q_{1,1}} (\nabla \times j_{\lambda A})\cdot \nabla \psi \xi dx\rt|\,.$$ The first term is estimated exactly as before, while the second one gives after integration by parts of $\nabla \times $ and using $\nabla \times \nabla=0$, $$\lt|\int_{Q_{1,1}} (\nabla \times j_{\lambda A})\cdot \nabla \psi \xi dx\rt|= \lt|\int_{Q_{1,1}} \nabla \xi\cdot  (j_{\lambda A}\times \nabla \psi) dx\rt|\le \|\partial_3 \xi \nabla' \psi\|_{L^2} \lt(\int_{Q_{1,1}} |j'_{\lambda A}|^2\rt)^{1/2},$$ from which we conclude the proof.

We now prove that for admissible pairs $(u,A)$ of bounded energy, the corresponding curves $x_3\mapsto \beta^{-1}B_3(\cdot,x_3)$ satisfy a sort of uniform Hölder continuity. This is the analog of for the limiting energy.

\[HolderE\] For every admissible pair $(u,A)$ with $\|\rho\|_\infty\le 1$, and every $x_3,{\widetilde}{x}_3\in(-1,1)$, letting $E:=\wE(u,A)$ it holds (recall ) $$\begin{gathered}
\label{HolderE2}
   \|\b^{-1}B_3(\cdot,x_3) - \b^{-1}B_3(\cdot,{\widetilde}{x}_3) \|_{BL}\les E^{1/2} |x_3-{\widetilde}{x}_3|^{1/2}
   + \sigma(\alpha,\beta,T)(E^{1/2}+E),
 \end{gathered}$$ where $\sigma(\alpha,\beta,T):= \lt(\frac{\a}{T}\rt)^{1/2} \lt(\a^2\b^{5/2}\rt)^{-1/6}+ (\a^{1/2}\b)^{-1/3} + \lt(\frac{\alpha^{1/3}}{T\b^{5/6}} \rt)^{1/2}$, which goes to zero in the regime . In particular, in that regime, if $E\les 1$, for every $x_3,{\widetilde}{x}_3\in(-1,1)$ with $|x_3-{\widetilde}{x}_3|\ge \sigma^{1/2}(\alpha,\b,T)$, there holds $$\label{equiconB}
  \|\beta^{-1}B_3(\cdot,x_3)-\beta^{-1}B_3(\cdot,{\widetilde}{x}_3)\|_{BL}\les |x_3-{\widetilde}{x}_3|^{1/2}.$$

The proof resembles that of [@CoOtSer Lem. 3.13]. First, we show that for every $Q_1-$periodic and Lipschitz continuous function $\psi$ with $\|\psi\|_{Lip}\le 1$, $$\label{firstestimHolder}
  \lt|\int_{Q_1\times\{x_3\} } \b^{-1}B_3 \psi dx' -\int_{Q_1\times\{{\widetilde}{x}_3\} } \b^{-1}B_3 \psi dx'\rt|\le|x_3-{\widetilde}{x}_3|^{1/2} \b^{-1/2} E^{1/2}.$$ This follows from $\Div B=0$ and integration by parts, which yields $$\begin{aligned}
 \lt|\int_{Q_1\times\{x_3\} } \b^{-1}B_3 \psi dx' -\int_{Q_1\times\{{\widetilde}{x}_3\} } \b^{-1}B_3 \psi dx'\rt|&=\lt|\int_{Q_1\times(x_3,{\widetilde}{x}_3)} \b^{-1}\partial_3 B_3 \psi  dx\rt|\\
 &=\lt|\int_{Q_1\times(x_3,{\widetilde}{x}_3)}\beta^{-1} B' \cdot \nabla' \psi  dx\rt|\\
 &\le |x_3-{\widetilde}{x}_3|^{1/2} \beta^{-1/2} \lt(\int_{Q_1\times(x_3,{\widetilde}{x}_3)} \beta^{-1}|B'|^2 dx\rt)^{1/2}\\
 &\le |x_3-{\widetilde}{x}_3|^{1/2} \beta^{-1/2}E^{1/2}.\end{aligned}$$ For $|x_3-{\widetilde}{x}_3|\le T^{-1} \a^{1/3}\b^{1/6}$, this implies that $$\|\b^{-1}B_3(\cdot,x_3) - \b^{-1}B_3(\cdot,{\widetilde}{x}_3) \|_{BL}\le \lt(\a^{1/3} T^{-1} \b^{-5/6}\rt)^{1/2} E^{1/2}\le \sigma(\a,\b,T) E^{1/2},$$ and is proven. Letting $\hat{\sigma}(\a,\b,T):= \lt(\frac{\a}{T}\rt)^{1/2} \lt(\a^2\b^{5/2}\rt)^{-1/6}+ (\a^{1/2}\b)^{-1/3}$, we are left to prove that for $Q_1-$periodic and Lipschitz continuous $\psi$ with $\|\psi\|_{Lip}\le 1$ and $|x_3-{\widetilde}{x}_3|\ge T^{-1}\a^{1/3}\b^{1/6}$, $$\label{toproveHolder}
  \lt|\int_{Q_{1}\times\{x_3\}} \beta^{-1}B_3 \psi dx' -\int_{Q_{1}\times\{{\widetilde}{x}_3\}} \beta^{-1}B_3 \psi dx' \rt|\les |x_3-{\widetilde}{x}_3|^{1/2}E^{1/2}+   \hat{\sigma}(\alpha,\beta,T)(E^{1/2}+E).$$ Up to translation we may assume ${\widetilde}{x}_3=0$ and $x_3>0$. Let $\delta\le x_3/2$ and define $\xi: \R\to \R$ by $$\xi(z):=\begin{cases}
             \frac{z}{\delta} & \textrm{if } 0<z<\delta\\
             1 & \textrm{if } \delta\le z \le x_3-\delta\\
             \frac{x_3-z}{\delta} & \textrm{if } x_3-\delta\le z\le x_3\\
             0 & \textrm{otherwise}.
            \end{cases}$$ We then have using again $\Div B=0$ and integration by parts $$\label{integpart}
 \int_{Q_{1,1}} \beta^{-1} B_3 \psi \partial_3\xi dx=-\int_{Q_{1,1}} \b^{-1}\rho  B'\cdot \nabla'\psi \, \xi dx -\int_{Q_{1,1}} \b^{-1}(1-\rho) B'\cdot\nabla'\psi\,  \xi dx.$$ The first term on the right-hand side of is estimated by . For the second term, we now estimate $$\label{firsttermright}
 \int_{Q_{1,1}} \b^{-1}(1-\rho)|B'|\xi dx\le \lt(\int_{Q_1\times(0,x_3)} \b^{-1}(1-\rho) dx\rt)^{1/2} \lt(\int_{Q_{1,1}}\b^{-1}(1-\rho)|B'|^2 dx \rt)^{1/2}.$$ We rewrite the first factor as $$\int_{Q_1\times(0,x_3)} \b^{-1}(1-\rho) dx=\int_{Q_1\times(0,x_3)} \b^{-1}B_3 dx  +\int_{Q_1\times(0,x_3)} \b^{-1}(B_3-(1-\rho)) dx,$$ from which we obtain $$\int_{Q_1\times(0,x_3)} \b^{-1}(1-\rho) dx\le \lt|\int_{Q_1\times(0,x_3)} \b^{-1}B_3 dx \rt| +|x_3|^{1/2} \lt(\int_{Q_{1,1}} \b^{-2} (B_3-(1-\rho))^2 dx\rt)^{1/2}.$$ This allows to make use of $$\int_{Q_{1,1}} \b^{-2} (B_3-(1-\rho))^2 dx\le \a^{-2/3}\b^{-4/3} E,$$ and $\int_{Q_1\times\{z\}} \beta^{-1} B_3 dx' =1$, yielding $$\label{firstestimrightbis}
  \int_{Q_1\times(0,x_3)} \b^{-1}(1-\rho) dx\le |x_3|+|x_3|^{1/2}\a^{-1/3}\b^{-2/3} E^{1/2}.$$ The second factor in is directly estimated by $$\int_{Q_{1,1}}\b^{-1}(1-\rho)|B'|^2 dx\le E,$$ so that inserting this and into gives $$\label{firstestimrightter}
\lt| \int_{Q_{1,1}} \b^{-1}(1-\rho) B'\cdot\nabla'\psi \xi dx\rt|\les |x_3|^{1/2} E^{1/2} +|x_3|^{1/4}  \a^{-1/6}\b^{-1/3}E^{3/4}.$$

Letting $f(z):=\int_{Q_1\times\{z\}}\beta^{-1} B_3 \psi dx'$, we thus obtain from , and , $$\label{firstestimf}
\lt|\int_0^{x_3} f \partial_3 \xi dz\rt|\les \frac{\a^{2/3}\b^{-1/6}}{T} E +\frac{\a^{1/3}\b^{-1/3}}{T} E^{1/2} \|\partial_3\xi\|_{L^2}+ |x_3|^{1/2} E^{1/2} +|x_3|^{1/4} E^{3/4} \a^{-1/6}\b^{-1/3}.$$ Since by definition of $\xi$, $\int_0^{x_3} f \partial_3\xi dz=\frac{1}{\delta} \int_0^\delta f dz-\frac{1}{\delta}\int_{x_3-\delta}^{x_3} f dz$, we have $$f(x_3)-f(0)=\int_0^{x_3} f \partial_3 \xi dz  +\frac{1}{\delta}\int_{0}^\delta (f-f(0)) dz+\frac{1}{\delta}\int_{x_3-\delta}^{x_3} (f(x_3)-f) dz,$$ so that $$|f(x_3)-f(0)|\le \lt|\int_0^{x_3} f \partial_3 \xi dz\rt|+\sup_{(0,\delta)}|f-f(0)| +\sup_{(x_3-\delta,x_3)}|f-f(x_3)|.$$ In view of this elementary inequality, the estimates and combine to $$\begin{gathered}
|f(x_3)-f(0)|\les |x_3|^{1/2} E^{1/2} +\frac{\a^{1/3}\b^{-1/3}}{T} E^{1/2} \delta^{-1/2} + \delta^{1/2}\b^{-1/2}E^{1/2} +\frac{\a^{2/3}\b^{-1/6}}{T} E\\ + |x_3|^{1/4} E^{3/4} \a^{-1/6}\b^{-1/3}.  \end{gathered}$$ We now optimize in $\delta $ by choosing $\delta=T^{-1}\a^{1/3}\b^{1/6}$, which combined with $\frac{\a^{2/3}\b^{-1/6}}{T}\ll \a^{-1/6}\b^{-1/3}$, yields in the form of $$|f(x_3)-f(0)|\les |x_3|^{1/2} E^{1/2} + (E^{1/2}+E)\lt( \lt(\frac{\a}{T}\rt)^{1/2} \lt(\a^2\b^{5/2}\rt)^{-1/6}+ \a^{-1/6}\b^{-1/3}\rt).$$

\[remarkequi\] We notice that thanks to the Kantorovich-Rubinstein Theorem [@villani Th. 1.14], if $B_3$ is non-negative then we can substitute the Bounded-Lipschitz norm in by a $1-$Wasserstein distance. In particular, it would imply that if $(\alpha_n,\beta_n,T_n)$ satisfy and if $(u_n,A_n)$ are admissible with $\|\rho_n\|_{\infty}\le 1$, $\wE(u_n,A_n)\les 1$, and $B^n_3\ge 0$, then the corresponding curves $x_3\mapsto \beta_n^{-1} B^n_3(\cdot,x_3)$ would be in some sense equi-continuous in the space of probability measures endowed with the Wasserstein metric.

For $\eps>0$ fixed, we define the following regularization of the singular double well potential $\chi_{\rho>0}(1-\rho)^2$ : $$\label{doublewell}
 W_\eps(\rho):=\eta_\eps(\rho)(1-\rho)^2  \ \text{with } \ \eta_\eps(\rho):=\min\{\rho/\eps,1\} 
  \,,$$ see and Figure \[figweps\]. We next show that the energy controls $W_\eps(\rho)$. Similar ideas have been used in the context of Bose-Einstein condensates [@GolMer].

![The cutoff function $\eta_\eps$ and the two-well potential $W_\eps$ used in Lemma \[lemmafirstlowerbound\].[]{data-label="figweps"}](fig-weps-crop){height="4cm"}

\[lemmafirstlowerbound\] For every $\eps>0$ there exists $C_\eps>0$ such that for every $(u,A)$ with $\|\rho\|_\infty\le 1$ it holds, $$\label{estimlower}
 \int_{Q_{1,1}} \a^{2/3}\b^{-2/3}W_\eps(\rho) dx
 \le \int_{Q_{1,1}} \a^{2/3}\b^{-2/3}\left(B_3 -(1-\rho)\right)^2 dx
 +C_\eps \wE(u,A) \frac{\a}{T}.$$

As above, to lighten notation, we let $E:= \wE(u,A)$. Writing $(1-\rho)= B_3 -(B_3-(1-\rho))$, we obtain by Young’s inequality $$\begin{aligned}
 (1-\rho)^2&= B_3 (1-\rho)-(B_3-(1-\rho))  (1-\rho)\\
&\le  B_3 (1-\rho) +\frac{1}{2}  (1-\rho)^2+\frac{1 }{2} (B_3-(1-\rho))^2.\end{aligned}$$ Multiplying by $2\eta_\eps(\rho)$ and using that $0\le\eta_\eps\le 1$ we obtain for $W_\eps(\rho)$ the estimate $$\label{wer}
  W_\eps(\rho) =\eta_\eps{(\rho)}   (1-\rho)^2\le (B_3-(1-\rho))^2+2\eta_\eps{(\rho)}  (1-\rho)B_3.$$ Let $\psi_\eps(s):=2\frac{\eta_\eps(s)}{s}(1-s)= 2\min \{ \frac{1}{\eps}, \frac{1}{s}\} (1-s) $ then $\psi_\eps$ is bounded by $1/\eps$ and is Lipschitz continuous in $s^{1/2} $ with constant of order $\eps^{-3/2}$ i.e. $\sup_t |(\psi_\eps(t^2))'|\les \eps^{-3/2}$. Since $2\eta_\eps {(\rho)} (1-\rho)B_3= \rho B_3 \psi_\eps(\rho)$, using with $\psi=\psi_\eps{(\rho)} $, we get $$\begin{aligned}
\left| \int_{Q_{1,1}} 2\eta_\eps {(\rho)} (1-\rho)B_3 dx\right|&\les \frac{\a^{1/3}\b^{2/3}}{T}\lt(\eps^{-1} E+  \a^{-1/3}\b^{-1/6}  E^{1/2}  \|{\nabla}' (\psi_\eps(\rho))\|_{L^2} \rt)  \\
&\les \frac{\a^{1/3}\b^{2/3}}{T}\lt(\eps^{-1}E+ \eps^{-3/2} \a^{-1/3}\b^{-1/6}  E^{1/2}  \|{\nabla}'\rho^{1/2}\|_{L^2} \rt) \\
&\les   C_\eps \frac{\a^{1/3}\b^{2/3}}{T}\lt(E+ \int_{Q_{1,1}} \a^{-2/3}\b^{-1/3} |\nabla' \rho^{1/2}|^2 dx\rt) \\
&\les C_\eps \frac{\a^{1/3}\b^{2/3}}{T} E,\end{aligned}$$ where we used that $|\nabla' \rho^{1/2}|\le |\nabla'_{\a^{1/3}\b^{1/3}TA} u|$ and thus $\int_{Q_{1,1}} \a^{-2/3}\b^{-1/3} |\nabla' \rho^{1/2}|^2 dx\le E$. Estimate follows from inserting this estimate into .

To prove the lower bound, we will need the following two dimensional result.

\[lowerbound2d\] Let $\chi_n\in BV(Q_1, \{0,\b_n^{-1}\})$ be such that $\lim_{n\to +\infty}\int_{Q_1}\chi_n dx'=1$ and $$\sup_{n} \int_{Q_1} \b_n^{1/2} |D' \chi_n| <+\infty.$$ Then, up to a subsequence, $\chi_n\weaklim \sum_i {\varphi}_i \delta_{X_i}$ for some at most countable family of ${\varphi}_i> 0$ and $X_i \in Q_1$, and $$\label{liminfperibeta}
\liminf_{n\to +\infty}  \int_{Q_1} \b_n^{1/2} |D' \chi_n|\ge 2\sqrt{\pi}\sum_i \sqrt{{\varphi}_i}.$$

[[*Step 1 (Compactness):*]{}]{} For each $n$ we split the cube $Q_1$ into small cubes of side length $3\b_n^{1/2}$. Let $Q_i^{n}$ be an enumeration of theses cubes such that $$\varphi_i^{n} :=\int_{Q_i^{n}} \chi_n dx'$$ is nonincreasing in $i$. Since $|\{\chi_n=\b_n^{-1}\}|=\b_n \int_{Q_1} \chi_n dx'=\b_n+o({\beta_n})$, we have $|{Q_i^{n}}\cap\{\chi_n=\b_n^{-1}\}|\le\b_n+o(\b_n)\le \frac{1}{2}|{Q_i^{n}}|$ and thus by the relative isoperimetric inequality [@AFP Th. 3.46], we have on each ${Q_i^{n}}$ $$\varphi_i^{n}=\int_{{Q_i^{n}}}  \chi_n  dx'\les\left( \int_{{Q_i^{n}}} \b_n^{1/2} |D'
\chi_n|\right)^2\,.$$ It follows that $$\sum_i \sqrt{\varphi_i^{n}} \les \int_{Q_1} \b_n^{1/2}|D' \chi_n|  \le C,$$ by the energy bound. Arguing as in the proof of Proposition \[existmu\], we deduce from Lemma \[lemsqrt\] that up to extracting a subsequence, $\chi_n \weaklim \sum_{i\in {\mathcal{I}}}{\varphi}_i \delta_{X_i}$ for some ${\varphi}_i> 0$ and $X_i\in Q_1$, pairwise distinct.

[[*Step 2 (Lower bound):*]{}]{} Assume now that the ${\varphi}_i$ are labeled in a decreasing order and fix $N\in \N$. Choose $r\in(0,1/4)$ sufficiently small so that $$\label{bxr}
\B'(X_i,r)\cap \B'(X_j,r)=\emptyset \qquad \qquad \forall i,j \le N \text{ with } i\neq j,$$ and let $\psi\in C^\infty_c(\B'(X_i, r);[0,1])$ be a smooth function such that $\psi=1$ in $\B'(X_i, r/2)$ and $|{\nabla}\psi|\le C/r$. Let $C_\mathrm{iso}=(2\sqrt{\pi})^{-1}$ be the isoperimetric constant in dimension 2, then we may write $$\begin{aligned}
\nonumber
\left(\int_{\B'(X_i, r/2)
} \b_n^{-1} \chi_n dx'\right)^\hal &
 \le &\left( \int_{Q_1} ( \psi \chi_n)^2 dx'\right)^{1/2} 
 \le  C_\mathrm{iso} \int_{Q_1} |D' (\psi \chi_n)|\\
 \nonumber & \le & C_\mathrm{iso} \left( \int_{\B'(X_i, r) }
 |D' \chi_n| + \frac{C}{r} \int_{\B'(X_i,r)\backslash \B'(X_i, r/2)}
  \chi_n dx'\right)\\
   \nonumber & \le & C_\mathrm{iso} \int_{\B'(X_i, r) }
 |D' \chi_n| + \frac{C}{r} 
 \end{aligned}$$ since $\int_{Q_1} \chi_n dx'\to 1$. Multiplying by $\beta_n^{1/2}$ and summing over $N$, we get $$\label{estimisop}
 \sum_{i=1}^N \left(\int_{\B'(X_i, r/2)
}  \chi_n dx'\right)^{1/2} \le C_\mathrm{iso} \int_{Q_1} \beta_n^{1/2}|D' \chi_n| + \beta_n^{1/2}\frac{CN}{r}.$$

Next observe that since $\chi_n\rightharpoonup 
 \sum_{i} {\varphi}_i \delta_{X_i}$, we have for every $i=1,\dots,N$, $$\liminf_{n\to +\infty}\lt(\int_{\B'(X_i, r/2)}  \chi_n dx'\right)^{1/2}=   {\varphi}_i^{1/2}
\,.$$ Therefore, passing to the limit in , we obtain $$\sum_{i=1}^N   {\varphi}^{1/2}_i\le  C_\mathrm{iso} \liminf_{n \to +\infty} \int_{Q_1} \beta_n^{1/2}|D' \chi_n|.$$ Since $N$ was arbitrary this implies .

With this lemma at hand, we can prove the compactness and lower bound result.

We fix for the proof a sequence $(u_n,A_n)$ with $\wE(u_n,A_n)\les 1$. We then let $B_n:=(B'_n,B_3^n):=\nabla \times A_n$ and $\rho_n:=|u_n|^2$.\
[[*Step 1 (Compactness):*]{}]{} Notice first that $$\label{bounddoublewell}
 \int_{Q_{1,1}} \lt( \frac{B^n_3-(1-\rho_n)}
 {\beta_n} 
 \rt)^2 dx\le \a_n^{-2/3} \b_n^{-4/3} \wE(u_n,A_n)
 \to0\,.$$ By [@CoOtSer Lem. 3.7] there is $\hat u_n$ with $\hat\rho_n:=|\hat u_n|^2=\min\{\rho_n,1\}$ such that $\wE(\hat u_n,A_n)\le (1+2\alpha_n/T_n)\wE(u_n,A_n)$ (the error comes from the last two terms in (\[eqEGLE\])). In particular, also $ (B^n_3-(1-\hat\rho_n))/\beta_n\to0$ in $L^2(Q_{1,1})$. Using this and $|B^n_3|\le |B^n_3-(1-\hat\rho_n)|$ on $\{ B^n_3\le 0\}$ we obtain $$\label{convneg}\lt|\int_{\{B^n_3<0\}} \b_n^{-1} B^n_3 dx\rt|\les \lt( \b_n^{-2} \int_{Q_{1,1}} (B^n_3-(1-\hat\rho_n))^2 dx\rt)^\hal   \le \frac{1}{(\alpha_n\beta_n^2)^{2/3} } \wE(u_n,A_n)   \to 0,$$ and since $\int_{Q_{1,1}} \b_n^{-1}B^n_3 dx=2$ the sequence $\b_n^{-1}B^n_3$ is bounded in $L^1$ and, after extracting a subsequence, $\b_n^{-1} B^n_3 \weaklim \mu$ for some measure $\mu$. From we also get $\b_n^{-1}(1-\rho_n)\weaklim \mu$, and the same for $\hat\rho_n$. It also follows from that $$\int_{Q_{1,1}} \b_n^{-1} (1- \hat\rho_n) dx= 2+\left|\beta_n^{-1} \int_{Q_{1,1}} \left(B_3^n-(1-\hat \rho_n) \right) dx\right|\to 2.$$ Moreover, since $\int_{Q_{1,1}} \b_n^{-1}(1-\hat\rho_n) |B_n'|^2 dx\le\int_{Q_{1,1}} \b_n^{-1} |B'_n|^2 dx\le \wE(u_n,A_n)$ it holds $$\int_{Q_{1,1}} \b_n^{-1}(1-\hat\rho_n) | B_n'| dx\le \lt(\int_{Q_{1,1}} \b_n^{-1}(1-\hat\rho_n) |B_n'|^2 dx\rt)^\hal \lt(\int_{Q_{1,1}} \b_n^{-1}(1-\hat\rho_n) dx\rt)^{\hal}\les 1,$$ thus (up to a subsequence) $\b_n^{-1}(1-\hat\rho_n) B_n'\weaklim m$ for some vector-valued measure $m$. By [@AFP Th. 2.34], $$\label{liminftrans}
 \liminf_{n\to +\infty} \int_{Q_{1,1}} \b_n^{-1}|B_n'|^2 dx\ge\liminf_{n\to +\infty} \int_{Q_{1,1}} \b_n^{-1}(1-\hat\rho_n)|B_n'|^2 dx\ge \int_{Q_{1,1}} \lt(\frac{dm}{d{\mu}}\rt)^2 d\mu,$$ and $m\ll \mu$. Moreover, from Lemma \[lemmameissner\] we have that $\b_n^{-1}\hat\rho_n B'_n\weaklim 0$ in a distributional sense and therefore $\b^{-1}_n B_n'$ itself converges to $m$. Letting $n\to +\infty$ in $$\b_n^{-1} \Div B_n=\partial_3 \left[\b_n^{-1} B^n_3\right]+\Divp \lt[\b_n^{-1}B_n'\rt]=0,$$ we obtain $\partial_3 \mu +\Divp m=0$. This proves (i).\
We now prove that $\mu=\mu_{x_3}\otimes dx_3$, that $\beta^{-1}_nB^n_3(\cdot,x_3)\weaklim \mu_{x_3}$ for a.e. $x_3\in(-1,1)$ and that $\mu_{x_3}\weaklim dx'$ as $x_3\to \pm 1$. By we have that, up to a subsequence in $n$, for a.e. $x_3\in (-1,1)$, $$\lim_{n\to +\infty} \int_{Q_1\cap \{B^n_3<0\}} \beta^{-1}_n B^n_3 dx=0.$$ Let $\mathcal{G}{\subseteq}(-1,1)$ be the set of $x_3$ for which this hold. For every $x_3\in \mathcal{G}$, the $L^1$ norm of $\beta^{-1}_nB^n_3(\cdot,x_3)$ is bounded thus if we fix a countable dense set $\mathcal{G}_d{\subseteq}\mathcal{G}$, we can assume up to extraction that for every $x_3\in \mathcal{G}_d$, $\beta^{-1}_n B^n_3(\cdot,x_3)\weaklim \nu_{x_3}$ for some probability measure $\nu_{x_3}$. For $x_3,{\widetilde}{x}_3\in \mathcal{G}_d$, thanks to the weak lower semi-continuity of $\|\cdot \|_{BL}$ and , $$\begin{aligned}
 \|\nu_{x_3}-\nu_{{\widetilde}{x}_3}\|_{BL}&\le \limsup_{n\to +\infty} \|\beta_n^{-1} B_3^n (\cdot,   x_3) - \beta_n^{-1} B_3^n(\cdot, {\widetilde}{x}_3)\|_{BL} \\
 &\les \lim_{n\to +\infty}\lt( |x_3-{\widetilde}{x}_3|^{1/2} +\sigma(\alpha_n,\beta_n,T_n)\rt) \\
 &= |x_3-{\widetilde}{x}_3|^{1/2}. \end{aligned}$$ Therefore, there exists a unique Hölder-continuous extension of $\nu_{x_3}$ to $(-1,1) \ni x_3 \mapsto \nu_{x_3} \in \mathcal{P}(Q_1)$. We claim that $\mu=\nu_{x_3}\otimes dx_3$. For $K\to +\infty$, let $\{z^K_j\}_{j=1}^{K+1} \in \mathcal{G}_d$ be an increasing sequence such that $|z^K_j-z^K_{j+1}|\les K^{-1}$, $|z_1^K+1|\les K^{-1}$ and $|z_{K+1}^K-1|\les K^{-1}$. Notice that for $n$ large enough, we have for every $j$ that $|z^K_j-z^K_{j+1}|\ge \sigma^{1/2}(\a_n,\b_n,T_n)$ (where $\sigma$ is defined in Lemma \[HolderE\]) so that applies. Let $\psi$ be a $Q_1-$periodic and Lipschitz continuous function on $Q_{1,1}$ with $\|\psi\|_{Lip}\le 1$. By the continuity of $x_3 \mapsto  \nu_{x_3}$, we have $$\begin{aligned}
 \int_{Q_{1,1}} \psi (d\mu- d\nu_{x_3}\otimes dx_3)&=\lim_{K\to +\infty}  \int_{Q_{1,1}} \psi d\mu-\sum_{j=1}^K  (z^K_{j+1}-z^K_j) \int_{Q_1} \psi(\cdot ,z_j^K)  d\nu_{z_j^K}\\
 &=\lim_{K\to +\infty} \lim_{n\to +\infty} \sum_{j=1}^K \int_{z_j^K}^{z_{j+1}^K} \int_{Q_1} \psi(x',x_3) \beta_n^{-1} B^n_3(x',x_3)\\
 &\qquad  \qquad \qquad \qquad \qquad  -\psi(x',z_j^K)\beta_n^{-1} B^n_3(x',z_j^K) dx' dx_3.\end{aligned}$$ Using the finite difference version of Leibniz’ rule, $$\begin{gathered}
\psi(x',x_3)  B^n_3(x',x_3)-\psi(x',z_j^K) B^n_3(x',z_j^K)\\
=   B^n_3(x',x_3)(\psi(x',x_3)-\psi(x', z_j^K))+  \psi(x',z_j^K)( B^n_3(x',x_3) - B^n_3(x',z_j^K))\end{gathered}$$ and using that $\|\psi\|_{Lip}\le 1$, we can estimate for fixed $K,j$ and $n$ large enough, $$\begin{aligned}
\lt| \int_{z_j^K}^{z_{j+1}^K} \rt.& \lt.\int_{Q_1} \psi(x',x_3) \beta_n^{-1} B^n_3(x',x_3)-\psi(x',z_j^K)\beta_n^{-1} B^n_3(x',z_j^K) dx' dx_3\rt|\\
&\le  \int_{z_j^K}^{z_{j+1}^K} \int_{Q_1} \beta_n^{-1} |B^n_3| |x_3-z^K_j| dx + \int_{z_j^K}^{z_{j+1}^K} \|\beta^{-1}_n( B^n_3(x',x_3) - B^n_3(x',z_j^K))\|_{BL}\\
&\les K^{-1} \int_{z_j^K}^{z_{j+1}^K} \int_{Q_1} \beta_n^{-1} |B^n_3|dx +K^{-1} K^{-1/2},\end{aligned}$$ where in the last line we have used that $|x_3-z^K_j|\les K^{-1}$ and . Summing this estimate over $j$, we obtain $$\left|\int_{Q_{1,1}} \psi (d\mu- d\nu_{x_3}\otimes dx_3)\right| \les \lim_{K\to +\infty} K^{-1}\lt[\lim_{n\to +\infty} \int_{Q_{1,1}}\beta_n^{-1} |B^n_3| dx + K^{-1/2}\rt]=0.$$ This establishes that $\mu= \nu_{x_3}\otimes dx_3$. Moreover, this proves that for every $x_3\in \mathcal{G}_d$, the whole sequence $\beta_n^{-1} B^n_3(\cdot,x_3)$ weakly converges to $\mu_{x_3}$. Since the set $\mathcal{G}_d$ was arbitrary, this proves the above convergence for all $x_3\in \mathcal{G}$.\
We finally show that the boundary conditions hold. For this we focus on $x_3=1$. For $x_3\in \mathcal{G}$, it holds by the weak lower semi-continuity of $\|\cdot \|_{BL}$, $$\begin{aligned}
 \|1-\mu_{x_3}\|_{BL}&\le \liminf_{n\to +\infty} \|1-\beta_n^{-1}B_3^n(\cdot, x_3)\|_{BL}\\
 & \le \liminf_{n\to +\infty} \|1-\beta_n^{-1}B_3^n(\cdot, 1)\|_{BL}+\limsup_{n\to +\infty} \|\beta_n^{-1}B_3^n(\cdot, 1)-\beta_n^{-1}B_3^n(\cdot, x_3)\|_{BL}.\end{aligned}$$ By , the second right-hand side term is controlled by $|1-x_3|^{1/2}$. For the first right-hand side term we note that because of $\|\psi\|_{H^{1/2}}{\stackrel{<}{\sim}}\|\psi\|_{Lip}$, we have $$\|\beta_n^{-1}B_3^n (\cdot, 1)- 1\|_{BL}^2  {\stackrel{<}{\sim}}\|\beta_n^{-1} B_3^n (\cdot, 1)-1\|_{H^{1/2}}^2 \le \alpha_n^{-1/3}\beta_n^{-7/6} \wE(\hat u_n, A_n).$$ Hence, it is the last assumption in that ensures that this term vanishes in the limit $n\to +\infty$. We thus obtain the desired estimate $$\| 1- \mu_{x_3}\|_{BL}{\stackrel{<}{\sim}}|1-x_3|^{1/2}.$$

[[*Step 2 (Lower bound and structure of $\mu$):*]{}]{} The starting point is an application of the usual Modica-Mortola trick. In this step we only deal with $\hat u_n$ and $\hat\rho_n$, and drop the hats for brevity. By and $|\nabla' \rho^{1/2}|\le |\nabla'_{\lambda A} u|$ we obtain from the Cauchy-Schwarz inequality : $$\label{estimMM}
 \left(1 + C_\eps \frac{\alpha_n}{T_n}\right)\wE(u_n,A_n)\ge \int_{Q_{1,1}} \b_n^{-1/2}2\sqrt{W_\eps(\rho_n)} |\nabla' \rho_n^{1/2}| +\b_n^{-1}|B_n'|^2 dx,$$ for any $\eps>0$. We momentarily fix a small $\delta>0$ and estimate by the co-area formula , $$\begin{aligned}
\int_{Q_{1,1}} 2\sqrt{W_\eps(\rho_n)}|\nabla' \rho_n^{1/2}| dx&\ge \int_{\delta}^{1-\delta} \int_{-1}^1  2\sqrt{W_\eps(s^2)} \H^1(\partial \{\rho_n(\cdot,x_3)>s^2\}) dx_3 ds.\end{aligned}$$ In particular there exists $s_n\in[\delta,1-\delta]$ depending on $n$ such that $$\begin{gathered}
\int_{\delta}^{1-\delta} 2\sqrt{W_\eps(s^2)} \int_{-1}^1  \H^1(\partial \{\rho_n(\cdot,x_3)>s^2\}) dx_3 ds\\
 \ge  \lt(\int_{\delta}^{1-\delta} 2\sqrt{W_\eps(s^2)} ds\rt)\int_{-1}^1  \H^1(\partial \{\rho_n(\cdot,x_3)>s^2_n\}) dx_3.\end{gathered}$$ Letting $\chi_n(x',x_3):=\b_n^{-1}(1-\chi_{\{\rho_n(\cdot,x_3)>s_n^2\}}(x'))$ this reads $$\label{boundBV}
 \int_{Q_{1,1}} \b_n^{-1/2}2\sqrt{W_\eps(\rho_n)}|\nabla' \rho_n^{1/2}| dx\ge   C_{\delta,\eps}\int_{Q_{1,1}} \b_n^{1/2} |D' \chi_n| dx_3,$$ where $C_{\delta,\eps}:=\int_{\delta}^{1-\delta} 2\sqrt{W_\eps(s^2)} ds$.

Let $\gamma_n\to 0$ to be chosen later. For $n$ large enough, if $\rho_n\le \gamma_n$, then $\chi_n=\beta_n^{-1}$ while if $\rho_n\ge 1-\gamma_n$, $\chi_n=0$ so that $$\begin{aligned}
\int_{Q_{1,1}} |\chi_n-\beta_n^{-1}(1-\rho_n)| dx&\le\beta_n^{-1} \int_{\{\rho_n\le \gamma_n\}} \rho_n dx  +\beta_n^{-1} \int_{\{\rho_n\ge 1-\gamma_n\}} (1-\rho_n) dx \\
& \qquad +\beta_n^{-1}|\{\gamma_n<\rho_n<1-\gamma_n\}|
\\ &\le 2\beta_n^{-1}\gamma_n +\beta_n^{-1}|\{\gamma_n<\rho_n<1-\gamma_n\}|.\end{aligned}$$ By definition of $W_\eps$ (recall , $\min_{[\gamma_n,1-\gamma_n]} W_\eps=\min(\frac{\gamma_n}{\eps}, \gamma_n^2)=\gamma_n^2$ so that using that $\int_{Q_{1,1}} \alpha_n^{2/3}\beta_n^{-2/3} W_\eps(\rho_n) dx\les 1$, $$\beta_n^{-1}|\{\gamma_n<\rho_n<1-\gamma_n\}|\le \beta_n^{-1} \gamma_n^{-2} \int_{Q_{1,1}} W_\eps(\rho_n) dx \les \beta^{-1/3}_n\gamma_n^{-2}\alpha_n^{-2/3}.$$ Therefore, if we choose $\gamma_n$ such that $\beta_n\gg \gamma_n\gg \alpha_n^{-1/3}\beta_n^{-1/6}$, which is possible since by hypothesis $\alpha_n \beta_n^{7/2}\to +\infty$, we obtain that $$\lim_{n\to +\infty}\int_{Q_{1,1}} |\chi_n-\beta_n^{-1}(1-\rho_n)| dx=0.$$ Combining this with , we obtain that $$\lim_{n\to +\infty}  \int_{Q_{1,1}}|\chi_n-\b_n^{-1}B_3^n| dx=0.$$ By Fubini, this implies that, after passing to a subsequence in $n$, for a.e. $x_3\in (-1,1)$, if $ \b_n^{-1}B_3^n(\cdot,x_3)\weaklim \mu_{x_3}$ then also $\chi_n(\cdot,x_3)\weaklim  \mu_{x_3}$. Moreover, from $\int_{Q_1} \b_n^{-1} B^n_3(x',x_3) dx'=1$ for a.e. $x_3\in(-1,1)$, we obtain that $\lim_{n\to +\infty}\int_{Q_1} \chi_n(x',x_3) dx'=1$ for a.e. $x_3\in(-1,1)$. We thus can use Lemma \[lowerbound2d\] to prove that $\mu_{x_3}=\sum_i \p_i \delta_{X_i}$ for some $\p_i>0$ and $$\label{lowerboundperim}
 \liminf_{n\to +\infty} \int_{Q_{1,1}} \b_n^{1/2} |D'\chi_n| dx_3 \ge 2\sqrt{\pi} \int_{-1}^1 \sum_i \sqrt{\p_i} dx_3\, ,$$ where we used Fatou lemma. This shows (iii). Putting , , and together we find $$\liminf_{n\to +\infty} \wE(u_n,A_n)\ge \int_{-1}^1 2\sqrt{\pi} C_{\delta,\eps} \sum_i \sqrt{{\varphi}_i} dx_3 +\int_{Q_{1,1}} \lt(\frac{dm}{d\mu}\rt)^2 d\mu,$$ for any $\eps$ and $\delta$. Since $$\lim_{\eps\to0}\lim_{\delta\to0}C_{\delta,\eps}= 2\int_0^1 (1-t^2) dt=\frac{4}{3},$$ and $K_*=\frac{8\sqrt{\pi}}{3}$, this concludes the proof of (iv).

Upper bound {#Sec:upperbound}
===========

In this section we construct a recovery sequence for any sequences $T_n, \alpha_n, \beta_n$ which obey (\[eqlbassumptcoeff\]) and additionally the condition of quantization of the total flux $$\label{defKtotal} 
  L_n^2 T_n\a_n \b_n\in 2\pi \N\,,$$ where $L_n:={\widetilde}L \a^{-1/3}_n \b_n^{-1/6}$. Recalling the form of the gradient term in the functional $\wE$ defined in (\[eqdefwe\]), to discuss quantization of the flux of individual domains it is convenient to introduce $$\label{defK} 
k_n:=\alpha_n^{1/3}\beta_n^{2/3}T_n\,.$$ The global flux quantization (\[defKtotal\]) then reads $k_n {\widetilde}L^ 2\in 2\pi\N$ and, in the ${\widetilde}L=1$ case we are considering here, simplifies to $k_n\in 2\pi \N$. Condition (\[eqlbassumptcoeff\]) implies $k_n\to +\infty$ and in particular $k_n/\beta_n\to+\infty$, so that the quantization condition becomes less and less stringent with increasing $n$. Aim of this section is to prove the following:

\[gammalimsup-v2\] Assume ${\widetilde}L=1$, (\[eqlbassumptcoeff\]) and (\[defKtotal\]). Then, for every $\mu$ with $I(\mu)<+\infty$ and $\mu_1=\mu_{-1}=dx'$ there exist sequences $u_n: Q_{1,1}\to {\mathbb{C}}$ and $A_n: Q_{1,1}\to \R^3$ such that $$\limsup_{n\to +\infty} \,  \wE(u_n,A_n)\le I(\mu).$$ The fields $\rho_n:=|u_n|^2$ and $B_n:= \nabla \times A_n$ are $Q_1$-periodic and it holds $$\b_n^{-1}(1-\rho_n)\weaklim \mu \qquad \textrm{and } \qquad \b_n^{-1} B_n'\weaklim m,$$ where $m$ is the measure such that $I(\mu)=I({\mu},m)$.

The idea of the construction is to use the density result Theorem \[theodens\] to separate the construction in two regions. In the bulk, the measure will be approximated by a finite polygonal measure for which the construction is made in Section \[secconstrbulk\]. In the boundary layer, we plug in the construction of [@CoOtSer], see Section \[secconstrboundary\], which is optimal up to a factor. Since the energy in the boundary layer is small, its suboptimal effect disappears in the limit.

We shall first construct the density $\rho$ and the magnetic field $B$. The appropriate energy is ${\widetilde}F_{\alpha,\beta}:={\widetilde}F_{\alpha,\beta}^{(-T,T)}+{\widetilde}F_{\alpha,\beta}^{\Ext}$, where $$\begin{gathered}
\label{eqdeftildaF}
  {\widetilde}F_{\alpha,\beta}^{(a,b)}(\rho,B):=\int_{Q_1\times(a,b)} \Bigl( \a^{-2/3}\b^{-1/3}\lt|\nabla'\rho^{1/2}\rt|^2+ 
 \a^{-4/3}\b^{-2/3}\lt|\partial_3\rho^{1/2}\rt|^2\\
+  \a^{2/3}\b^{-2/3}\left({B_3} -(1-\rho)\right)^2
+  \b^{-1}|{B'}|^2 \Bigr)dx\,,\end{gathered}$$ and $$\label{eqdeftildaFext} 
  {\widetilde}F_{\alpha,\beta}^{\Ext}(B):=
  \alpha^{1/3}\beta^{7/6}
  \| \beta^{-1}B_3-1\|^2_{H^{-1/2}(Q_1\times \{\pm1\})}\,.$$ Their sum corresponds to the energy $\wE$, up to a reconstruction of $u$ and $A$ that will be discussed in Section \[recoveryGL\]. A pair $(\rho,B)$ is admissible for ${\widetilde}F_{\a,\b}^{(a,b)}$ if $\rho$ and $B$ are $Q_1$-periodic, $\div B=0$ and $\int_{Q_1} B_3(x',x_3)dx'=\beta$ for all $x_3$.

We say that a pair $(\rho,B)$ is $k$-quantized if there is a closed set $\omega{\subseteq}Q_{1,1}$ such that $B=0$ outside $\omega$, $\rho=0$ in $\omega$, and the flux of $ \b^{-1}B_3$ over every connected component of $\omega\cap\{x_3=z\}$ is an integer multiple of $2\pi /k$, for all $z\in(-1,1)$.

Construction in the bulk {#secconstrbulk}
------------------------

This section is concerned with the local construction of flux tubes. For notational simplicity we present this construction in $\R^2\times(a,b)$, without the periodicity assumption; since $\rho=1$ and $B=0$ outside a small region, its periodic extension is immediate (see proof of Proposition \[propubregular\] below). We start from the optimal profile at the boundary of the individual tubes, with the lengths measured in units of the coherence length (recall ). For the purpose of the upcoming constructions, we cut the profile at a lengthscale $R\gg1$ towards the normal region.

\[lemmaVR\] Consider the functional $\displaystyle G(v):=2\sqrt\pi\int_0^{+\infty} |\dot{v}|^2+(1-v^2)^2 dt$. Then, $$\inf\left\{ G(v): v(0)=0, \lim_{t\to+\infty} v(t)=1\right\}=K_*=\frac83\sqrt\pi\,.$$ Furthermore, for all $R\ge3$ there is $v_R\in C^\infty(\R;[0,1])$ with $v_R(t)=0$ for $t\le 0$, $v_R(t)=1$ for $t\ge R$, $|\dot{v}_R|\le 2$, and, setting $K_R:=G(v_R)$, one has $$\lim_{R\to\infty} K_R=\frac83\sqrt\pi,$$ and $\int_0^{+\infty} t \left(|\dot{v}_R|^2+(1-v_R^2)^2\right) dt{\stackrel{<}{\sim}}1$.

The lower bound follows from the usual Modica-Mortola type computation $$\int_0^{+\infty} |\dot{v}|^2+(1-v^2)^2dt\ge 2\int_0^{+\infty} (1-v^2)\dot{v}dt=2\int_0^1 (1-s^2) ds = \frac43\,.$$ To prove the upper bound, we recall that $v(t):= \tanh t$ is the minimizer of $G$ under the constraint $v(0)=0$. A direct computation shows that $G(v)=K_*$. We then define for $R>0$, $$\hat v_R(t):=
 \begin{cases}
 0 & \text{ if } t<1/R,\\
  \displaystyle\frac{\tanh (t-1/R)}{\tanh (R-2/R)}&\text{ if } t\in[1/R,R-1/R],\\
  1 & \text{ if } t>R-1/R,
 \end{cases}$$ and $v_R:=\psi_{1/R}\ast \hat v_R$, with $\psi_{1/R}\in C^\infty_c(-1/R,1/R)$ a mollifier. By construction, this has the desired properties and verifies $G(v_R)\to G(v)$ as $R\to\infty$.

We start with the simple case in which the limiting measure comes from a Lipschitz curve; in practice this will be used only for affine or piecewise affine curves.

\[lemmacurve3\] Let $X: (a,b)\to \R^2$ be a Lipschitz curve, ${\varphi}>0$, $R>0$. We define $$\rho(x):=  v_R^2\left(\frac{|x'-X(x_3)|-\sqrt{\beta{\varphi}/\pi}}{\eta}\right),$$ where $\eta:=\alpha^{-2/3}\beta^{1/6}\ll 1$ is the coherence length (see ) and $v_R$ was introduced in Lemma \[lemmaVR\], and define $B$ through $$B_3(x):=\chi_{\B'(X(x_3),\sqrt{\beta {\varphi}/\pi} )}(x')\,, \hskip1cm
B'(x):=B_3(x) \dot{X}(x_3) \,.$$ Then, $\rho B=0$ almost everywhere, $\div B=0$, and $${\widetilde}F^{(a,b)}_{\alpha,\beta}(\rho,B)\le  \left(1+ \frac{C}{({\varphi}\alpha^{4/3}\beta^{2/3})^{1/2}}
+ \frac{C}{{\varphi}\alpha^{4/3}\beta^{2/3}}\right)
\int_a^b \bigl( K_R\sqrt{{\varphi}}+{\varphi}|\dot{X}|^2\bigr) dx_3.$$ The constant $C$ is universal. Moreover, if $\mu:={\varphi}\delta_{X(x_3)}\otimes dx_3$, $$\label{distestim}
 W_2^2(\b^{-1} B_3, \mu)\le |b-a|\frac{\beta {\varphi}^2}{2\pi}.$$

The condition $\rho B=0$ follows from $v_R(t)=0$ for $t\le 0$. To check the divergence condition, pick $\psi\in C^1_c(\R^2\times(a,b))$ and compute $$\begin{aligned}
 \int_{\R^2\times(a,b)} (B_3\partial_3 \psi  + B'\cdot \nabla'\psi) dx
 &=\int_a^b  \int_{\B'(X(x_3), \sqrt{\beta{\varphi}/\pi})} 
( \partial_3 \psi + \nabla'\psi \cdot \dot{X}(x_3)) dx'dx_3\\
& =\int_a^b  \int_{\B'(0, \sqrt{\beta{\varphi}/\pi})} 
 \frac{d}{dx_3} \psi(X(x_3)+y', x_3) dy' dx_3=0\,.\end{aligned}$$ We now estimate the energy. We start from the interfacial energy at fixed $x_3$ (see ), $$\begin{aligned}
E_1(x_3):=\int_{\R^2} \alpha^{-2/3}\beta^{-1/3} |\nabla'\rho^{1/2}|^2 + \alpha^{2/3}\beta^{-2/3}
 (B_3-(1-\rho))^2  dx'\,.\end{aligned}$$ If $|x'-X(x_3)|<  \sqrt{\beta{\varphi}/\pi}$ then $\rho=0$ and $B_3=1$, whereas if $|x'-X(x_3)|>\eta R +\sqrt{\beta{\varphi}/\pi}$ then $\rho=1$, $B_3=0$. In the intermediate region, we use $|\nabla'\rho^{1/2}|=|\dot{v}_R|/\eta$. Passing to polar coordinates and using $r=|x'-X(x_3)|$ as an integration variable, $$\begin{aligned}
E_1(x_3)=
 \int_{\sqrt{\beta{\varphi}/\pi}}^{\sqrt{\beta{\varphi}/\pi}+\eta R} \Bigl[
 \frac{\alpha^{-2/3}\beta^{-1/3}}{\eta^2} |\dot{v}_R|^2
 + \alpha^{2/3}\beta^{-2/3} (1-v_R^2)^2\Bigr]
 \lt(\frac{r-\sqrt{\beta{\varphi}/\pi}}\eta\rt)2\pi r dr \,.\end{aligned}$$ We change variables according to $r=\sqrt{\beta{\varphi}/\pi}+s\eta$, insert the definition of $\eta$, and by Lemma \[lemmaVR\] obtain $$\begin{aligned}
E_1(x_3)&=
 \int_0^R \beta^{-1/2}
 \left[|\dot{v}_R|^2+ (1-v_R^2)^2\right]2\pi (\sqrt{\beta{\varphi}/\pi}+s\eta) ds \\
& \le  K_R \sqrt{\varphi}+ C  \alpha^{-2/3}\beta^{-1/3}
 \le  \sqrt{\varphi}\left(K_R+\frac{C}{{\varphi}^{1/2}\alpha^{2/3}\beta^{1/3}}\right)\,.\end{aligned}$$ The other contributions to the energy are the cost of transport and the vertical part of the gradient, $$\begin{aligned}
E_2(x_3):=\int_{\R^2} \alpha^{-4/3}\beta^{-2/3} |\partial_3\rho^{1/2}|^2 + \beta^{-1}|B'|^2 dx'\,.\end{aligned}$$ By definition of $B'$, we have for the second term, $$\begin{aligned}
\int_{\R^2}  \beta^{-1}|B'|^2 dx'= |\dot{X}(x_3)|^2 {\varphi}\,.\end{aligned}$$ For the first one we use $|\partial_3\rho^{1/2}|\le |\dot{v}_R||\dot{X}|/\eta$ and change variables as above to obtain $$\begin{aligned}
\int_{\R^2} \alpha^{-4/3}\beta^{-2/3} |\partial_3\rho^{1/2}|^2  dx'
&{\stackrel{<}{\sim}}\frac{ \alpha^{-4/3}\beta^{-2/3}}{\eta} |\dot{X}(x_3)|^2
\int_0^R (  \sqrt{\beta{\varphi}/\pi}+s\eta)|v_R'|^2 ds\\
&{\stackrel{<}{\sim}}\alpha^{-4/3}\beta^{-2/3} |\dot{X}(x_3)|^2 (1 + \sqrt{\beta{\varphi}}/\eta)\\
&=  |\dot{X}(x_3)|^2 {\varphi}\left(
\frac{1}{{\varphi}\alpha^{4/3}\beta^{2/3}}+
\frac{1}{{\varphi}^{1/2}\alpha^{2/3}\beta^{1/3}}
\right)
\,.\end{aligned}$$ To prove we consider the transport map $T(x',x_3)=(X(x_3),x_3)$ which gives $$W_2^2(\b^{-1}B_3,\mu)\le \int_a^b \b^{-1} \int_{\B'(X(x_3),\sqrt{\b {\varphi}/\pi})} |x'-X(x_3)|^2 dx'dx_3=|b-a| \frac{\b {\varphi}^2}{2\pi}.$$

We now turn to the construction around branching points. Since the total length around branching points is small, the construction here does not need to achieve the optimal constant but only the optimal scaling. The idea of the construction is the following. We first transform disks into squares, then split the square into two rectangles and then retransform each rectangle into a disk. The construction is sketched in Figure \[fig1\]. We start with the transformation from a rectangle to a disk.

![Construction around a branching point[]{data-label="fig1"}](figsep2-crop){width="6cm"}

\[lemmacurveendpoints\] Let $a\le b\in \R$, $\gamma>0$, $X\in \R^2$, ${\varphi}>0$, and $R\ge 3$. Let then $Q {\subseteq}\R^2$ be a rectangle with side lengths $w$ and $h$ centered in $X$, such that $wh=\beta{\varphi}$, $w/h+h/w\le \gamma$ and $\a^{-4/3}\b^{-2/3}\le {\varphi}$. Let as before $\eta:=\alpha^{-2/3}\beta^{1/6}$ be the coherence length.

Then there are $\rho\in L^\infty(\R^2\times [a,b];[0,1])$ and $B\in L^2(\R^2\times [a,b];\R^3)$ such that $\div B=0$, $\rho B=0$, with $B=0$ and $\rho=1$ on $(\R^2\setminus \B'(X,r))\times[a,b]$ for some $r\sim \eta R + \sqrt{\beta{\varphi}}$, and $${\widetilde}F^{(a,b)}_{\alpha,\beta}(\rho, B) {\stackrel{<}{\sim}}\sqrt{\varphi} |a-b| R^2+ \frac{r^2}{|a-b|}\varphi R^2,$$ where the implicit constants only depend on $\gamma$. Further, $\rho$ and $B$ satisfy the boundary conditions $$\rho(x',a)=  v_R^2\left(\frac{|x'-X|-\sqrt{\beta{\varphi}/\pi}}{\eta}\right),
\hskip1cm
 \rho(x',b)=  \min\left\{1,\frac{\dist^2(x', Q)}{\eta^2}\right\}\,,$$ and $$B_3(x',a)= \chi_{\B'(X, \sqrt{\beta{\varphi}/\pi})}(x')\,,
\hskip3mm\text{ }\hskip3mm
B_3(x',b)= \chi_{Q}(x')\,.$$

We note that our assumptions on the parameters $w, h, {\varphi},$ just mean $w\sim h\sim \sqrt{\beta {\varphi}}\ge \eta$. That is, the thread diameter is large compared to the coherence length. Note that $r\sim \eta R+\sqrt{\beta {\varphi}}$ behaves like the maximum between the thread diameter $w\sim h$ and the cut-off scale $\eta R$. The proof of Lemma \[lemmacurveendpoints\] is based on an explicit construction for a bilipschitz bijection with unit determinant that transforms a rectangle into a circle, which we first present.

![The construction in Lemma \[lemmasquarecircle\] transforms a rectangle into a circle, keeping $\det \nabla'u=1$. Each curve corresponds to a different value of $x_3$, and plots the solution of $\hat r(r,\theta)=t$, which corresponds to $r(\theta)=t/[\lambda(x_3)\cos(\theta x_3)]$ for $\theta\in (-\pi/4,\pi/4)$.[]{data-label="fig-square"}](fig-square-crop){width="5cm"}

\[lemmasquarecircle\] Assume that $z_-,z_+,h,w>0$, $X\in\R^2$ are given, $z_-< z_+$. Then there is $u:\R^2\times[z_-,z_+]\to\R^2$ such that $$u(x',z_-)=x'\,,\hskip5mm
  u(X+(-\frac12w,\frac12w)\times(-\frac12h,\frac12h),z_+)=\B'(X,\sqrt{hw/\pi})\,,\hskip5mm
  \det \nabla' u =1 \text{ a.e.}\,.$$ The function $x\mapsto (u(x),x_3)$ is bilipschitz, its inverse is of the form $y\mapsto (U(y),y_3)$, and $u(X,x_3)=X$ for all $x_3$. If additionally $h/w+w/h\le \gamma$, then the bounds $|\nabla'u|+|\nabla'U|{\stackrel{<}{\sim}}1$, $|\partial_3 u|+|\partial_3 U|{\stackrel{<}{\sim}}h/|z_+-z_-|$ hold, with constants which only depend on $\gamma$.

By scaling and translation we may assume that $hw=\pi$, $z_+=1$, $z_-=0$, $X=0$. We can further assume $h=w=\sqrt\pi$, as the general case is obtained by taking the composition at each $x_3$ with the linear map ${\mathop{\mathrm {diag}}}(g(x_3), 1/g(x_3))$, where $g(x_3):=\sqrt{h/w}(1-x_3)+x_3$. We work in polar coordinates, and construct functions $\hat r$, $\hat\theta$ of $r$, $\theta$ and $x_3$ such that $$u(r\cos\theta,r\sin\theta,x_3)=(\hat r\cos\hat\theta, \hat r\sin\hat\theta),$$ with $r\ge 0$ and $0\le\theta\le \pi/4$, and then extend by symmetry. We set $$\hat \theta = f(\theta, x_3) \,,\hskip5mm \hat r = r \lambda(x_3) \cos(x_3\theta),$$ where $\lambda$ and $f$ are two functions still to be determined (see Figure \[fig-square\]). The extension of the $\frac{1}{8}-$sectors by reflection is feasible provided that $f(0,x_3)=0$ and $f(\pi/4,x_3)=\pi/4$ for all $x_3$, the boundary data are attained provided that $f(\theta,0)=\theta$, $\lambda(0)=1$, $\lambda(1)=2/\sqrt\pi$. The latter ensures that indeed the straight segment $r \cos \theta=\frac{w}{2}=\frac{\sqrt{\pi}}{2}$ is mapped into the unit circle $\hat{r}=1$. The determinant condition is equivalent to $$1=\frac{\hat r}{r}\partial_r \hat{r}  \partial_\theta f= \lambda^2(x_3)\cos^2(x_3\theta) \partial_\theta f(\theta, x_3),$$ which can be solved (using $f(0,x_3)=0$) to give $$f(\theta,x_3):=\frac{1}{x_3\lambda^2(x_3)}\tan (\theta x_3)\,,$$ smoothly extended to $f(\theta,0)=\theta/\lambda^2(0)$. The condition $f(\pi/4,x_3)=\pi/4$ determines $\lambda$, $$\lambda(x_3):=\left(\frac{\tan (\pi x_3/4)}{\pi x_3/4}\right)^{1/2},$$ which obeys $\lambda(1)=2/\sqrt\pi$ and smoothly extends to $\lambda(0)=1$. Clearly $\lambda\sim 1$ so that $\partial_\theta f\sim 1$. This implies that the change of variables defines a smooth deformation of the $\frac{1}{8}-$sector which smoothly depends on $x_3\in[0,1]$.

After a rotation and a translation, we may assume that $X=0$, $a=0$ and $b>0$; so that $Q=(-w/2,w/2)\times (-h/2,h/2)$. We treat two regions separately.

In the lower region $\R^2\times [0,b/2]$, we interpolate between $\B'_0:=\B'(0, \sqrt{\beta{\varphi}/\pi})$ and $Q$. To do this, let $u$ and $U$ be the functions from Lemma \[lemmasquarecircle\], using it for the rectangle $Q$ and $z_-=0$, $z_+=b/2$. The magnetic field is defined by $$B_3(x):=\chi_{\B'_0}(u(x))\,,\hskip1cm
  B'(x):=\chi_{\B'_0}(u(x))\partial_3 U(u(x),x_3)\hskip3mm \text{ for } x_3\in[0,b/2]\,.$$ The density is defined by $$\rho(x):=v_R^2\left(\frac{|u(x)|-\sqrt{\beta{\varphi}/\pi}}{\eta}\right) \hskip3mm \text{ for } x_3\in[0,b/2]\,.$$ The condition $\rho B=0$ follows immediately, as well as the boundary data at $x_3=0$. Since $B=0$ and $\rho=1$ whenever $|u(x)|\ge \eta R + \sqrt{\beta{\varphi}/\pi}$ and since $|u(x)|\ge |x'|/\|\nabla'U\|_\infty$, we have $B=0$ and $\rho=1$ whenever $|x'| \ge \|\nabla'U\|_\infty  (\eta R + \sqrt{\beta{\varphi}})$. In order to check $\div B=0$, we fix $\psi\in C^1_c(\B'_r\times(0,b/2))$. Performing then a change of variables at each $x_3$ gives, since $\det \nabla'U=1$, $$\begin{aligned}
1
 \int_{\R^2\times(0,b/2)} \left[\partial_3 \psi B_3 + \nabla'\psi\cdot B' \right]dx &=
 \int_{\R^2\times(0,b/2)} (\partial_3\psi B_3 + \nabla' \psi \cdot  B')(U(y),y_3) dy
 \\
 &= \int_{\R^2\times(0,b/2)} \chi_{\B'_0} \left[\partial_3 \psi( U(y),y_3)  + \partial_3 U(y)\cdot \nabla'\psi( U(y),y_3))   \right] dy \\
 &= \int_{\B'_0} \int_{(0,b/2)} \frac{d}{dx_3} (\psi(U(y),y_3)) dy_3 dy'=0\,.\end{aligned}$$ By the properties of $u$ we obtain $B_3(x',b/2)=\chi_Q(x')$.

In the upper region $\R^2\times[b/2,b]$, we keep $B(x',x_3):=e_3 \chi_Q(x')$ concentrated on $Q$ and linearly interpolate the profile between $\rho^{1/2}(x',b/2)$ and the profile $\rho^{1/2}(x',b)$ given in the statement: $$\rho^{1/2}(x):=\frac{2}{b}\lt[
\rho^{1/2}(x',b/2) (b-x_3)+ \rho^{1/2}(x',b)(x_3-b/2)\rt]
 \hskip3mm \text{ for } x_3\in(b/2,b)\,.$$ It is immediate to check that $B_3$ and $\rho$ match continuously at all interfaces, $\rho B=0$ and $\div B=0$. This concludes the construction.

We estimate the energy similarly to the proof of Lemma \[lemmacurve3\]. Lemma \[lemmasquarecircle\] gives $|\nabla'u|{\stackrel{<}{\sim}}1$, $|\partial_3 u|{\stackrel{<}{\sim}}h/b$, $|\partial_3 U|{\stackrel{<}{\sim}}h/b$, with constants depending only on $\gamma$. This yields $|B'|\les h/b$. Furthermore, by Lemma \[lemmaVR\] $|v_R|\le1$ and $|\dot{v}_R|\le 2$, so that $\rho\le 1$, $|\nabla'\rho^{1/2}|{\stackrel{<}{\sim}}1/\eta $ and $|\partial_3\rho^{1/2}|{\stackrel{<}{\sim}}h/(\eta b)$.

We start with the region $x_3\in [0,b/2]$. All integrals in $x'$ can be restricted to the set $$\Omega_{x_3}:=\{x': \sqrt{
 \beta{\varphi}/\pi}\le |u(x',x_3)|\le \eta R+\sqrt{\beta{\varphi}/\pi}\}.$$ Since $\det \nabla u'=1$, we have $$|\Omega_{x_3}|=|\B'(0,\eta R+\sqrt{\beta{\varphi}/\pi})\setminus \B'_0|
= \pi \eta^2R^2+2\sqrt{\beta{\varphi}\pi}\eta R{\stackrel{<}{\sim}}\eta \sqrt{\beta{\varphi}} R^2,$$ (in the last step we used the assumption on ${\varphi}$ and $R\ge 3$). Therefore $$\begin{aligned}
   {\widetilde}F_{\alpha,\beta}^{(0,b/2)}(\rho,B)
   &{\stackrel{<}{\sim}}\int_{0}^{b/2}\int_{\Omega_{x_3}} \left( \frac{\alpha^{-2/3}\beta^{-1/3}}{\eta^2}
   + \frac{h^2\alpha^{-4/3}\beta^{-2/3}}{\eta^2 b^2} + \alpha^{2/3}\beta^{-2/3} + \frac{h^2}{\beta b^2}\right) dx' dx_3\\
   &{\stackrel{<}{\sim}}\int_{0}^{b/2}|\Omega_{x_3}| dx_3 \left( \alpha^{2/3}\beta^{-2/3} +  \frac{h^2}{\beta b^2}\right) \\
   &{\stackrel{<}{\sim}}\left( b + \frac{h^2}b \alpha^{-2/3}\beta^{-1/3}\right) R^2 \sqrt{{\varphi}}\,\\
   &{\stackrel{<}{\sim}}\sqrt{{\varphi}} b R^2 +\frac{r^2}{b} {\varphi}R^2,
   \end{aligned}$$ where in the last line we have used that $h^2\les r^2$ and $\a^{-2/3}\b^{-1/3}\les \sqrt{{\varphi}}$. The region $(b/2,b)$ is simpler, as the $|B'|^2$ term does not appear, the others are the same with the exception of $|\partial_3 \rho^{1/2}|\les 1/b$, which is smaller by the factor $\eta/h\les 1$.

Using this building block we can finally produce the construction that will be used at branching points.

\[Lembranch3\] Let $a<b\in \R$, $\gamma>0$, $\ell>0$, $R\ge 3$, $X_i\in\R^2$ and $\varphi_i>0$ for $i=0,1,2$ with ${\varphi}_0={\varphi}_1+{\varphi}_2$, ${\varphi}_0/{\varphi}_1+{\varphi}_0/{\varphi}_2\le \gamma$ and $\sqrt{\beta {\varphi}_0}\ge \eta$, where as above $\eta:=\alpha^{-2/3}\beta^{1/6}$.

Then, if $|X_0-X_1|,|X_0-X_2|\le \ell/4$, $\sqrt{\beta \varphi_0}\ll |X_1-X_2|$ and $ \eta R\ll \ell$, then there are $\rho\in L^\infty(\R^2\times [a,b];[0,1])$ and $B\in L^2(\R^2\times [a,b];\R^3)$ such that $\div B=0$, $\rho B=0$, $B=0$ and $\rho=1$ on $(\R^2\setminus \B'(X_0,3\ell/4))\times[a,b]$, $${\widetilde}F^{(a,b)}_{\alpha,\beta}(\rho, B) {\stackrel{<}{\sim}}\sqrt{\varphi_0} |a-b| R^2+ \frac{\ell^2}{|a-b|}\varphi_0R^2,$$ where the implicit constants only depend on $\gamma$, and which satisfy the boundary conditions$$\rho(x',a)=v_R^2\left(\frac{|x'-X_0|-\sqrt{\beta{\varphi}_0/\pi}}{\eta}\right)
\hskip1mm\text{  }\hskip1mm
\rho(x',b)=\min_{i=1,2}v_R^2\left(\frac{|x'-X_i|-\sqrt{\beta{\varphi}_i/\pi}}{\eta}\right)\,,$$ and $$B_3(x',a)= \chi_{\B'(X_0, \sqrt{\beta{\varphi}_0/\pi})}(x')\,,
\hskip3mm\text{ }\hskip3mm
B_3(x',b)=\sum_{i=1,2} \chi_{\B'(X_i, \sqrt{\beta{\varphi}_i/\pi})}(x')\,.$$

Note that our assumptions on the parameters ${\varphi}_0, {\varphi}_1, {\varphi}_2$ and $h$ just mean that $\eta\les \sqrt{\beta {\varphi}_i}\ll |X_1-X_2|\les \ell$ and $\eta R\ll \ell$. That is, the thread diameter is at least as large as the coherence length $\eta$ but small compared to the distance between the threads. Likewise, $\ell$ is large compared to the cut-off scale $\eta R$.

After a translation and a rotation we may assume that $a=0$, $b>0$, $X_0=0$, $X_2-X_1=\zeta e_1$, with $\zeta>0$. Let $$w_0:=h:=(\b {\varphi}_0)^{1/2}, \qquad w_1:= \frac{{\varphi}_1}{{\varphi}_0} w_0 \qquad \textrm{ and } \qquad w_2:= \frac{{\varphi}_2}{{\varphi}_0}w_0.$$ Let then $Q^i:=X_i+[-w_i/2,w_i/2]\times[-h/2,h/2]$, for $i=0,1,2$. Notice that since $|X_i|\le \ell/4$ and $w_i\le w_0=h\ll \ell$, $Q^i{\subseteq}\B'(3\ell/8)$.\
We divide the interval $(0,b)$ in three parts (see again Figure \[fig1\]). In $(0,b/3)$, we apply Lemma \[lemmacurveendpoints\] to transform $\B'(\sqrt{\beta{\varphi}_0/\pi})$ into $Q^0$ (in particular we have $\rho=1$ on $\lt(\R^{2}\backslash \B'(\frac{3\ell}{4})\rt)\times (0,b/3)$). In $(b/3,2b/3)$, we connect $Q^0$ to $Q^1$ and $Q^2$ by an explicit construction (see below). Finally, in $(2b/3,b)$ we apply again Lemma \[lemmacurveendpoints\] to transform $Q^1$ and $Q^2$ back into $\B'(X_1, \sqrt{\beta {\varphi}_1/\pi})$ and $\B'(X_2, \sqrt{\beta {\varphi}_2/\pi})$. Notice that in $(2b/3,b)$, if $\rho(x',x_3)\ne 1$ then by Lemma \[lemmacurveendpoints\] and our hypothesis on the parameters, necessarily$$|x'|\le \max_i |X_i|+O(\eta R+\sqrt{\beta{\varphi}_0})\le \frac{\ell}{4} +o(\ell).$$ Thus, $\rho=1$ on $\lt(\R^{2}\backslash \B'(\frac{3\ell}{4})\rt)\times (2b/3,b)$.\

It only remains to discuss the construction in the central region. Let $y_1:=X_1+w_2e_1/2$, $y_2:=X_2-w_1e_1/2$ and, for $i=1,2$, ${\widetilde}Q^i:=Q^i-y_i$, so that up to a null set, $Q^0$ is the disjoint union of ${\widetilde}Q^1$ and ${\widetilde}Q^2$ (see Figure \[figseparsquares\]). Since $y_2-y_1=X_2-X_1-w_0e_1/2$, and by assumption $w_0=(\beta{\varphi}_0)^{1/2} \ll|X_1-X_2|$, we have $(y_2-y_1)\cdot e_1>0$.\
For $i=1,2$ and $x_3\in(b/3,2b/3)$ we set ${\widetilde}{Q}^i(x_3):= {\widetilde}{Q}_i+\frac{x_3-b/3}{b/3} y_i$. Since $(y_2-y_1)\cdot e_1>0$, we have ${\widetilde}{Q}^1(x_3)\cap{\widetilde}{Q}^2(x_3)=\emptyset$ for all $x_3$. Furthermore, since $Q^i{\subseteq}\B'(3\ell/8)$, also ${\widetilde}{Q}^i(x_3){\subseteq}\B'(3\ell/8)$ for $x_3\in(b/3,2b/3)$. We finally let $$\rho^{1/2}(x):=\min \lt\{ 1,\eta^{-1}\dist(x',{\widetilde}{Q}^1(x_3)\cup {\widetilde}{Q}^2(x_3))\rt\}\hskip3mm \text{ for } x_3\in[b/3,2b/3],$$ and correspondingly $$(B',B_3)(x):=
  \sum_{i=1,2}
   \lt( \frac{y_i}{b/3},1\rt)\chi_{{\widetilde}{Q}^i(x_3)}(x')\,.$$ All admissibility conditions are easily checked. In particular, $\rho=1$ if $\dist(x',{\widetilde}{Q}^1\cup{\widetilde}{Q}^2(x_3))\ge  \eta $ which holds if $|x'|\ge 3\ell/4$. The energy estimate is immediate.

Boundary layer {#secconstrboundary}
--------------

\[propboundraylayer\] Let $N\in\N$, $\alpha,\beta,T,t>0$ be given. Let $k:=\a^{1/3}\b^{2/3}T$ and let ${\varphi}_1, \dots, {\varphi}_{N^2}$ be positive numbers such that $$k {\varphi}_i \in 2\pi \N \, \qquad \textrm{and } \qquad\sum_{i=1}^{N^2} {\varphi}_i=1.$$ Assume that in the regime $\beta\le \alpha$, $\alpha \sqrt{2}/T\le 1$, we have $t\gg \alpha^{-1}$, $1/N\ges \alpha^{-2/3}\beta^{-1/3}$ and ${\varphi}_i\sim 1/N^2$ for every $i$. Then there are $\rho$ and $B$, admissible for ${\widetilde}F_{\alpha,\beta}^{(0,t)}$ and $k$-quantized, such that $$\label{eqboundarylint}
{\widetilde}F_{\alpha,\beta}^{(0,t)}(\rho,B)
\les tN + \frac{1}{N^2t},$$ $$\label{eqboundarylext}
\alpha^{1/3}\beta^{7/6}
  \|\beta^{-1}{B_3}- 1\|_{H^{-1/2}(x_3=0)}^2\les \frac{\beta^{1/3}}{\alpha^{1/3}} + \frac{\beta^{1/2}}{N^2t}
  +\alpha^{1/3}\beta^{7/6},$$ and, denoting by $Q^i$ the squares centered in the $N^2$ points of the square grid of spacing $N^{-1}$ with $|Q^i|=\beta {\varphi}_i$, $$\label{eqrandwb3bdly}
  B_3(x',t)=\sum_i  \chi_{Q^i}(x') \text{ and }
 \rho^{1/2}(x',t)=\min\{ 1,\eta^{-1} \dist(x',\cup_i Q^i) \}\,.$$

Note that the assumptions on the parameters $N$ and ${\varphi}_i$ mean that the sidelength $\sqrt{\beta {\varphi}_i}$ of the squares is larger than the coherence length $\a^{-2/3}\b^{1/6}$ (see ) and that both of them are small compared to the distance (equal to $1/N$) between the squares. The assumption on the parameter $t$ means that the thickness $t$ of the boundary layer is large with respect to the coherence length, where we recall that vertical and horizontal lengths have different units.

The key construction is described in [@CoOtSer Lem. 4.7]. However, the notation and the scalings are different in that paper. Indeed, since there was no need to rescale $x_3$ by the thickness and $x'$ by the distance between the threads $T \alpha^{-1/3} \beta^{-1/6}$, in that paper length was measured in terms of the penetration length. We first let $r_i$ be a square of side $1/N$ centered on the square grid of spacing $N^{-1}$, $Q^i{\subseteq}r_i$ be a square with the same center and area given by $ |Q^i|=\beta{\varphi}_i$, and $b_i:= \beta{\varphi}_iN^2$. Denoting with a star the quantities from [@CoOtSer], we set $T_*:=Tt/2$, $\kappa_*:=\alpha \sqrt2/T$, $L_*:=T\alpha^{-1/3}\beta^{-1/6}$, $r_i^*:=L_*r_i$, $\hat r_i^*:=L_*Q^i$, $d_0^*:=L_*/N$, $\rho_0^*:=\beta^{1/2}L_*/N$, $N_*=N$, and $ b_i^*:=(|\hat r_i^*|/|r_i^*|)\kappa_*/\sqrt2  =b_i \kappa_*/\sqrt2$.

[@CoOtSer Lem. 4.7] gives $Q_{L_*}$-periodic fields $\chi_*\in BV_\loc(\R^2\times (0,T_*);\{0,1\})$ and $B_*=(B'_*,B_3^*)\in L^2_\loc(\R^3;\R^3)$ such that $\div B_*=0$, $B_*(1-\chi_*)=0$, and, $$\chi_*(x',T_*)=\sum_i \chi_{\hat r_i^*}(x')  \qquad \textrm{and} \qquad B^*_3(x',T_*)=\frac{\kappa_*}{\sqrt2} \sum_i \chi_{\hat r_i^*}(x') \hskip5mm\text{ for } x'\in Q_{L_*}\,,$$ with energy $$\label{estimstar}
 \frac{1}{L_*^2} \int_{Q_{L_*}\times(0,T_*)} \kappa_* |D\chi_*|+|B_*'|^2+\chi_* \lt(B_3^*-\frac{\kappa_*}{\sqrt2}\rt)^2 dx \les
E_*^{int}:=  \lt( \frac{\kappa_* \rho_0^* T_*}{\sqrt2 } +  \frac{\kappa_*^2(\rho_0^*d_0^*)^2}{2T_*}\rt)\lt(\frac{N}{L_*}\rt)^2$$ and $$\frac{1}{L_*^2}  \|B_3^*-\sum_i b_i^*\chi_{q_j^*}\|^2_{H^{-1/2}(Q_{L_*})}\les
 E_*^{ext}:=\lt(
 \kappa_*(\rho_0^*)^2+
      \frac{\kappa_*^2(\rho_0^*)^3 d_0}{T_*}\rt)\lt(\frac{N}{L_*}\rt)^2\,.$$ We extend $\chi_*$ and $B_*$ to $Q_{L_*}\times(T_*,2T_*)$ by $$\chi_*(x):=\chi_*(x',T_*) \qquad \textrm{and} \qquad B_*=(0,B_3^*(x',T_*)) \hskip5mm\text{ for } x'\in Q_{L_*}\times(T_*,2T_*),$$ so that holds in $Q_{L_*}\times(0,2T_*)$. We set $$\rho_*(x):=\min\{1,\kappa_*^2d(x,\omega_*)^2\},$$ where $\omega_*=\{x: \chi_*(x)=1\}$ (here $d(\cdot,\omega_*)$ is the 3D distance, periodic in the tangential directions). Since $\omega_*$ is invariant in the $x_3-$direction inside $(T_*,2T_*)$ and since $t\gg \a^{-1}$, giving $T_*\gg \kappa_*^{-1}$, for $x_3=2T_*$, $$\min\{1,\kappa_*^2d(x,\omega_*)^2\}=\min\{1,\kappa_*^2\dist(x',\omega_*\cap\{x_3=2T_*\})^2\} \qquad \textrm{for } x=(x',2T_*).$$ Hence, $$\rho_*(x',2T_*)=\min_{i} \lt\{\min(1,\kappa_*^2\dist(x',\hat{r}^i_*)\rt\}.$$ The same computation as in the proof of [@CoOtSer Th. 4.9] leads to $$\begin{aligned}
1
 \frac{1}{L_*^2} \int_{Q_{L_*}\times(0,2T_*)}  |\nabla\rho_*^{1/2}|^2+|B_*'|^2+ (B_3^*-\frac{\kappa_*}{\sqrt2}(1-\rho_*))^2 dx \les
E_*^{int}\,.
 \end{aligned}$$ We scale back in the tangential direction to obtain $$\rho(x',x_3):=\rho^*(L_*x',x_3)\quad \text{ and } \quad 
 (B',B_3)(x',x_3):=\frac{\sqrt2}{\kappa_*}\lt(\frac{B'_*}{L_*},B_3^*\rt)(L_*x',x_3)\,.$$ We have $\div B=0$ and (\[eqrandwb3bdly\]) holds. Changing variables gives $\nabla'\rho^{1/2}(x)=L_*\nabla'\rho_*^{1/2}(L_*x',x_3)$ and $\partial_3\rho^{1/2}(x)=\partial_3 \rho^{1/2}_*(L_*x',x_3)$, so that $$\begin{gathered}
{\widetilde}F^{(0,t)}_{\alpha,\beta}(\rho,B)\le 
\alpha^{-4/3}\beta^{-2/3} \frac{1}{L_*^2} \int_{Q_{L_*}\times(0,2T_*)} |\nabla \rho_*^{1/2}|^2 + |B'_*|^2 +
\lt(B_3^*-\frac{\kappa_*}{\sqrt2}(1-\rho_*)\rt)^2 dx\\
\les\alpha^{-4/3}\beta^{-2/3}  E_*^{int}.\end{gathered}$$ Since $E_*^{int}=\alpha^{4/3}\beta^{2/3} (Nt/2+2/(N^2t))$, this concludes the estimate for ${\widetilde}F^{(0,t)}_{\alpha,\beta}$.

We now estimate the boundary term. By the embedding of $L^\infty$ into $H^{-1/2}$ (recall that since $\sum_i {\varphi}_i =1$, $\int_{Q_1} \sum_i(b_i-\beta) \chi_{r_i} dx'=0$) we have $$\|\sum_i (b_i-\beta)\chi_{r_i}\|^2_{H^{-1/2}(Q_{1})}\les
 \|b_i-\beta\|_\infty^2  \,.$$ Rescaling the boundary estimate for $B_3^*$ leads to $$\|B_3-\sum_i b_i\chi_{r_i}\|^2_{H^{-1/2}(Q_1)}\les
 \frac{1}{L_*^3}  \frac2{\kappa_*^2}\|B_3^*-\sum_i b_i^*\chi_{r_j^*}\|^2_{H^{-1/2}(Q_{L_*})}\les
\frac{\beta^{7/6}}{\alpha^{2/3}} + \frac{\beta^{4/3}}{\alpha^{1/3} N^2t}$$ Adding terms and using that $\sum_i \chi_{r_i}\equiv 1$, we conclude that $$\alpha^{1/3}\beta^{7/6}
  \|\beta^{-1}{B_3}- 1\|_{H^{-1/2}(Q_1)}^2\les \frac{\beta^{1/3}}{\alpha^{1/3}} + \frac{\beta^{1/2}}{N^2t}
  +\frac{\alpha^{1/3}}{\beta^{5/6}} \|b_i-\beta\|_\infty^2  \,.$$

Back to the full GL functional {#recoveryGL}
------------------------------

In this section, we relate the functional ${\widetilde}F_{\alpha,\beta}$, as defined in (\[eqdeftildaF\]–\[eqdeftildaFext\]), with the functional ${\widetilde}E_T$, as defined in (\[eqdefwe\]).

\[thirdupperbound\] Let $k=\alpha^{1/3}\beta^{2/3}T$ and let $\rho,B$ be $k$-quantized admissible functions for ${\widetilde}F_{\alpha,\beta}$ with $\rho B=0$. Let $\omega:=\{\rho=0\}$ and then for $z\in(-1,1)$, $\omega_{z}:=\omega\cap \{x_3=z\}$. Assume that $\omega$ is closed, that $\omega^c$ is connected and that for every $x_3\in(-1,1)$, $\omega_{x_3}^c$ is also connected. Then, there is a pair $(u,A)$, admissible for ${\widetilde}E_T$, such that $\rho= |u|^2$, $B=\nabla\times A$, and $$\label{limsupT}
 {\widetilde}E_T(u,A)={\widetilde}F_{\a,\b}(\rho,B).$$

We shall construct a function $A$ on $\R^2\times (-1,1)$ so that $B=\nabla\times A$ everywhere, and then a multivalued function $\theta$ on the set $\Omega^c:=\Z^2\times\{0\}+\omega^c$ such that $\nabla\theta=\b^{-1}kA$. Here $B$ and $\rho$ are $Q_1$-periodic, but $A$ and $\theta$ not necessarily. Setting $u:=\rho^{1/2} e^{i\theta}$ we then have $\nabla_{\b^{-1}kA}u= e^{i\theta}\nabla\rho^{1/2}$, which directly implies (\[limsupT\]). Notice that thanks to the hypothesis on $\omega$, the set $\Omega^c$ is a connected open set.

We start from the construction of $A$. By [@CoOtSer Lem. 4.8] applied to $B-\beta e_3$ there is a $Q_1$-periodic potential $A_\mathrm{per}$ such that $\nabla \times A_\mathrm{per}=  B- \b e_3$ and $\Div A=0$. We define $ A(x):= A_\mathrm{per}(x)+\b x_1 e_2$ so that $\nabla \times  A=B$. We remark that on any open set where $B=0$ the vector field $A$ is curl-free and divergence-free, therefore harmonic, and in particular smooth. In particular, since $\Omega^c$ is open and since by the Meissner condition, $B=0$ in $\Omega^c$, $A$ is smooth in $\Omega^c$.

We now turn to the existence of $\theta$. For a fixed level $x_3$, let $h+\omega^i_{x_3}$, $h\in\Z^2$, $i=1,..., I(x_3)$ be the connected components of $(\R^2\times\{x_3\}) \cap (\Z^2+\omega_{x_3})$. Denote the flux going through $\omega^i_{x_3}$ by $$\Phi^i(x_3):=\int_{\omega^i_{x_3}} B_3\, dx'.$$ By assumption we have $$\label{condfluxi}
\b^{-1}k \Phi^i(x_3)\in 2\pi \Z \qquad \textrm{ for every $i$. }$$

Fix a smooth curve $\Gamma_0:=\{(\gamma_0(x_3),x_3) \ : \ x_3\in(-1,1)\}{\subseteq}\Omega^c$. For $x_3,y_3\in(-1,1)$, let $\Gamma_0^{x_3,y_3}:=\{(\gamma_0(t),t) \ : \ t\in(x_3,y_3)\}{\subseteq}\Gamma_0$. For $x=(\gamma_0(x_3),x_3)$, let $\theta(x):= \b^{-1}k \int_{\Gamma_0^{0,x_3}}  A\cdot \tau d\H^1$. Now, for a generic $x=(x',x_3)\in \Omega^c$, let $$\theta(x):= \theta(\gamma_0(x_3),x_3)+\b^{-1}k\int_{\Gamma^x}  A\cdot \tau d\H^1,$$ where $\Gamma^x$ is any (horizontal) curve in $\Omega^c\cap(\R^2\times\{x_3\})$ connecting $x$ to $(\gamma_0(x_3),x_3)$. This gives a well defined $\theta:\Omega^c\to\R/2\pi\Z$ since for every closed (horizontal) curve $\Gamma$ in $\Omega^c\cap(\R^2\times\{x_3\})$, $$\label{condfluxcurvebis}
  \b^{-1}k\int_{\Gamma}  A\cdot \tau d\H^1\in 2\pi \Z.$$ Indeed, this follows by Stokes’ Theorem and . Let us show that $\nabla \theta=\b^{-1}k A $ in $\Omega^c$. Let $x=(x',x_3)$ be fixed and let $\Gamma^x$ be a fixed simple smooth curve joining $x$ to $(\gamma_0(x_3),x_3)$ inside $\Omega^c\cap(\R^2\times\{x_3\})$. Since $\Omega^c$ is open, there is a simply connected neighborhood $V$ of $\Gamma^x$ such that $V{\subseteq}\Omega^c$. Let then $y=(y',y_3)\in V$. Upon shrinking $V$, we may assume that $(\gamma_0(y_3),y_3)\in V$. Let $\Gamma^y{\subseteq}V$ be a smooth curve joining $y$ to $ (\gamma_0(y_3),y_3)$ and let $\Gamma^{x,y}{\subseteq}V$ be a smooth curve joining $x$ to $y$. By definition, we have $$\theta(x)=\theta(\gamma_0(x_3),x_3)+ \b^{-1}k \int_{\Gamma^x}A\cdot \tau d\H^1, \quad \theta(y)=\theta(\gamma_0(y_3),y_3)+ \b^{-1}k\int_{\Gamma^y}  A\cdot \tau d\H^1$$ and $$\theta(\gamma_0(y_3),y_3)=\theta(\gamma_0(x_3),x_3)+ \b^{-1}k\int_{\Gamma_0^{x_3,y_3}}  A\cdot \tau d\H^1,$$ so that $$\theta(y)-\theta(x)= \b^{-1}k\int_{\Gamma^y} A\cdot \tau d\H^1+ \b^{-1}k\int_{\Gamma_0^{x_3,y_3}} A\cdot \tau d\H^1-\b^{-1}k\int_{\Gamma^x}  A\cdot \tau d\H^1.$$ However, by Stokes Theorem (and $B=0$ in $V$), $$\int_{\Gamma^y}A\cdot \tau d\H^1+\int_{\Gamma_0^{x_3,y_3}} A\cdot \tau d\H^1-\int_{\Gamma^x}  A\cdot \tau d\H^1=\int_{\Gamma^{x,y}} A\cdot \tau d\H^1,$$ so that $$\theta(y)-\theta(x)=\b^{-1}k\int_{\Gamma^{x,y}}  A\cdot \tau d\H^1,$$ proving that indeed, $\nabla \theta=\b^{-1}k A$ in $\Omega^c$.

Proof of the upper bound
------------------------

We start from a construction for an $N$-regular measure and finite $R$.

\[propubregular\] Let $\mu\in \MNreg(Q_{1,1})$ for some $N\in\N$, $R\ge 3$, and assume (\[eqlbassumptcoeff\]) and (\[defKtotal\]) hold. Then, there exist sequences $u_n: Q_{1,1}\to {\mathbb{C}}$ and $A_n: Q_{1,1}\to \R^3$ such that $$\label{estimenerregular}
\limsup_{n\to +\infty} \,  \wE(u_n,A_n)\le 
\frac{K_R}{K_*}
I(\mu)
+ \frac{C}{N^{1/2}}(1+I(\mu)),$$ with $|u_n|$ and $\nabla\times A_n$ $Q_1-$periodic. Here $K_R$ and $K_*$ are as in Lemma \[lemmaVR\], $C$ is universal. Moreover, $$\label{estimW2regular}
 \limsup_{n\to +\infty} W_2^2(\b^{-1}_n B^n_3,\mu)\les  N^{-3/2}.$$

We first modify slightly $\mu$ in order to be able to use Proposition \[propboundraylayer\] for the boundary layer. Set $\varepsilon_N:=N^{-3/2}$. For $x_3\in [-1+\eps_N,1-\varepsilon_N]$ let $\hat{\mu}_{x_3}:=\mu_{x_3/(1-\varepsilon_N)}$ and for $x_3\in [-1,-1+\eps_N]\cup[1-\varepsilon_N,1]$, let $\hat{\mu}_{x_3}:=\mu_1=\mu_{-1}=N^{-2}\sum_l \delta_{X_l}$, where the $\{X_l\}_{l=1}^{N^2}$ form a regular grid with spacing $1/N$ (recall Def. \[defregular\]). We then have $$I(\hat{\mu})\le \frac{1}{1-\varepsilon_N} I(\mu)+2K_*\varepsilon_N N.$$ Moreover, since $|\hat \mu|(Q_1\times (1-\varepsilon_N,1))= \varepsilon_N$, the same on the other side, and $$W_2^2((1-\varepsilon_N)\mu,\hat{\mu}\LL(Q_1\times (-1+\varepsilon_N,1-\varepsilon_N))\les \eps_N^2,$$ so that $$W_2^2(\mu, \hat{\mu})\les \eps_N+\eps^2_N\les N^{-3/2},$$ it is enough to prove the estimates for $\hat{\mu}$ instead of $\mu$. We start by characterizing the geometry of the construction, which will not depend on $n$. The measure $\hat{\mu}$ is supported on finitely many polygonal curves, parametrized by $X_i:[a_i,b_i]\to Q_1$, which are disjoint up to the endpoints and carry a flux ${\varphi}_i>0$. The endpoints are either on the boundary of $Q_{1,1}$, or they are triple points. Those on the boundary constitute a regular grid. We let $\phimin:=\min_i{\varphi}_i$, $\phimax:=\max_i {\varphi}_i$, and $\gamma:=8\phimax/\phimin$. Since there are finitely many curves, these quantities are finite and positive. We define as before $\eta_n:=\alpha_n^{-2/3}\beta_n^{1/6}$ to be the coherence length; by (\[eqlbassumptcoeff\]) we have $\eta_n\to0$.

Let $y_j:=(Y_j,z_j)\in Q_{1,1}$ denote the internal endpoints of the curves. For $\ell>0$ sufficiently small one has that, for any $j$, the only curves which intersect $\B'_\ell(Y_j)\times(z_j-\ell,z_j+\ell)$ are those with an endpoint in $y_j$. Since $y_j$ is a triple point, there are three such curves $\Gamma_0$, $\Gamma_1$ and $\Gamma_2$, intersecting only at $y_j$. If we let $M$ be the maximal slope of all curves and $t_\ell:=\ell/(8M)$, then no curve intersects $\partial\B'_{\ell/8}(Y_j)\times (z_j-t_\ell,z_j+t_\ell)$ (this means, they all “exit” from the top and bottom faces). Without loss of generality, we can assume that $(X_0,z_j-t_\ell)=\Gamma_0\cap\lt( \B'(Y_j,\ell)\times\{z_j-t_\ell\}\rt)$, $(X_1,z_j+t_\ell)=\Gamma_1\cap \lt(\B'(Y_j,\ell)\times\{z_j+t_\ell\}\rt)$ and $(X_2,z_j+t_\ell)=\Gamma_2\cap \lt(\B'(Y_j,\ell)\times\{z_j+t_\ell\}\rt)$ (see Figure \[figtriplepoint\]). By definition of $t_\ell$, it holds $|X_0-Y_j|\le \frac{\ell}{8}$ and $|X_i-X_0|\le \frac{\ell}{4}$ for $i=1,2$. We set $\omega_j:=\B'_{\ell}(Y_j)\times (z_j-t_\ell,z_j+t_\ell)$ and let $\delta_\ell$ be the minimum distance between any two curves outside $\cup_j\omega_j$.

To explain the strategy, we first carry out a construction which ignores quantization. For sufficiently large $n$ we have $ \eta_nR+\sqrt{\phimax \beta_n/\pi}<\delta_\ell/2$ and can thus use the construction of Lemma \[lemmacurve3\] in a $\delta_\ell/2$-neighborhood of each curve (outside of $\cup_j \omega_j$), extending by $\rho_n=1$ and $B_n=0$ to the complement. In the cylinders $\omega_j$, since the geometry is fixed, for $n$ sufficiently large the conditions $\sqrt{\beta_n \phimin}\ge \eta_n$, $\sqrt{\beta_n\phimax}\ll |X_1-X_2|$ and $ \eta_n R\ll \ell$ are satisfied and we can use Lemma \[Lembranch3\]. We then have that $\rho_n=1$ outside $\B'(X_0, 3\ell/4)\times(z_j-t_\ell,z_j+t_\ell){\subseteq}\omega_j$.

In order to obtain a quantized field, we define $k_n$ as in (\[defK\]), which obeys $ k_n\in 2\pi\N$ and $k_n\to+\infty$. Let $\hat{\mu}^n$ be the $k_n$-quantized approximation of $\hat{\mu}$, as given by Lemma \[lemmaquantize\]. For sufficiently large $n$ we have $\frac12 \hat{\mu}\le \hat{\mu}^n\le 2\hat{\mu}$ and for $|x_3|\ge 1-\varepsilon_N$, $$\hat{\mu}_{x_3}^n=\sum_l{\varphi}_l^n\delta_{X_l},$$ where $|{\varphi}_l^n-1/N^2|\le C(\hat{\mu})/k_n$, ${\varphi}_l^n\le 2/N^2$. The fluxes ${\varphi}_i^n$ obey $\frac12{\varphi}_i\le {\varphi}_i^n\le 2{\varphi}_i$.

We then construct $B_n$ and $\rho_n$ using Lemma \[lemmacurve3\] and Lemma \[Lembranch3\] as discussed above. The geometry is the one determined by $\hat{\mu}$. In particular, the points $y_j$, and the constants $M$, $\phimin$, $\phimax$, $\delta_\ell$ and $t_\ell$ do not depend on $n$, and only $\delta_\ell$ and $t_\ell$ depend on $\ell$. Adding terms gives for sufficiently large $n$ (as a geometry dependent function of $\ell$ and $R$) $$\begin{aligned}
 {\widetilde}F_{\alpha_n,\beta_n}^{(-1,1)}(\rho_n,B_n)\le &
\sum_i   \left(1+ \frac{C}{({\varphi}_i^n \alpha_n^{4/3}\beta_n^{2/3})^{1/2}}
+ \frac{C}{{\varphi}_i^n \alpha_n^{4/3}\beta_n^{2/3}}\right)
\int_{a_i}^{b_i} \bigl( K_R\sqrt{{\varphi}_i^n}+{\varphi}_i^n |\dot{X}_i|^2\bigr) dx_3
\\
&+C(\gamma) \sum_j \sqrt{\phimax} t_\ell R^2+ \frac{\ell^2}{t_\ell}\phimax R^2,\end{aligned}$$ where the first sum runs over all curves and the second over the cylinders and where $C(\gamma)>0$ is a constant depending on $\gamma$. In the limit $n\to+\infty$ we have ${\varphi}_i^n\to{\varphi}_i$, $\alpha_n^2\beta_n\to+\infty$ and therefore, inserting the definition of $t_\ell$, $$\begin{aligned}
\limsup_{n\to+\infty} {\widetilde}F_{\alpha_n,\beta_n}^{(-1,1)}(\rho_n,B_n)\le &
\sum_i   \int_{a_i}^{b_i} \bigl( K_R\sqrt{{\varphi}_i}+{\varphi}_i |\dot{X}_i|^2\bigr) dx_3
+C(\gamma) R^ 2 \ell \sum_j  \left(\frac{\sqrt{\phimax}}{M}  +M\phimax\right)\,.\end{aligned}$$ The sum over $j$ depends only on $\hat{\mu}$. Therefore if $\ell$ is chosen sufficiently small we have $$\limsup_{n\to +\infty} \,  {\widetilde}F_{\alpha,\beta}^{(-1,1)}(\rho_n,B_n)\le \frac{K_R}{K_*} I(\mu)+\frac1{N^{1/2}}\,.$$ Moreover, thanks to and Lemma \[lemmaquantize\], we have $$\begin{aligned}
 W_2^2(\beta^{-1}_n B^n_3,\hat{\mu})&\les W_2^2(\beta^{-1}_n B^n_3,\hat{\mu}^n)+W_2^2(\hat{\mu},\hat{\mu}^n)\\
 &\les C(\hat{\mu})\b_n  + t_\ell \, \sharp\{Y_j\}+ C(\hat{\mu})k_n^{-1}\\
 &\les C(\hat{\mu})\b_n  + \frac{\ell}{M} \, \sharp\{Y_j\}+ C(\hat{\mu})k_n^{-1}.\end{aligned}$$ Again, since $ \sharp\{Y_j\}/M$ only depends on $\hat\mu$, if we choose $\ell$ sufficiently small then $$\label{estimW2inside}
\limsup_{n\to+\infty} W_2^2(\beta^{-1}_n B^n_3,\hat{\mu})\le N^{-3/2}.$$ We finally address the boundary layer, focusing for definiteness on the side $x_3>0$. We first apply once more Lemma \[lemmacurveendpoints\] to each curve for $x_3\in (1-\eps_N, 1-\eps_N/2)$, so that the resulting fields $B_n$ and $\rho_n$ obey for all $x_3\in (1-\eps_N/2,1)$, $$B_n(x)=\sum_l e_3 \chi_{Q^l}(x') \text{ and }
 \rho_n(x)= \min\{1,\eta_n^{-1} \dist(x',\cup_l Q^l)\},$$ where $Q^l$ are squares centered in the $N^2$ points $X_l$ with $|Q^l|= \b{\varphi}_l^n $ , and $|{\varphi}_l^n-1 /N^2|\le C(\hat{\mu})/k_n$. Then we modify $\rho_n$ and $B_n$ in the set $x_3\in (1-\eps_N/2,1)$ using Proposition \[propboundraylayer\]. This results in new fields $\hat \rho_n$, $\hat B_n$ which obey $${\widetilde}F_{\alpha_n,\beta_n}^{(1-\eps_N/2,1)}(\hat\rho_n,\hat B_n)+
  \alpha_n^{1/3}\beta_n^{7/6}
  \|\beta_n^{-1}{(\hat{B}_n)_3}- 1\|_{H^{-1/2}(x_3=1)}^2\les
  \frac{1}{N^{1/2}}+
  \frac{\beta_n^{1/3}}{\alpha_n^{1/3}} + \frac{\beta_n^{1/2}}{N^{1/2}}
  +\frac{\alpha_n^{1/3}}{\beta_n^{5/6}} \frac{C(\hat\mu)N^4}{k_n^2}  \,.$$ By (\[eqlbassumptcoeff\]) and (\[defK\]) all terms up to the first one tend to zero as $n\to+\infty$. Therefore $$\limsup_{n\to +\infty} \,  {\widetilde}F_{\alpha_n,\beta_n}(\hat\rho_n,\hat B_n)\le \frac{K_R}{K_*}
I(\hat{\mu})
+\frac{C}{N^{1/2}}\,.$$ Finally, we pass from $(\hat\rho_n,\hat B_n)$ to $(u_n,A_n)$ via Proposition \[thirdupperbound\] and conclude the proof of . We also obtain since $$W^2_2(\beta^{-1}_n (\hat{B}_n)_3, \hat{\mu})\les W^2_2(\beta^{-1}_n B^n_3, \hat{\mu})+N^{-3/2}.$$

It only remains to combine the different steps.

By density (Theorem \[theodens\]) there is a sequence $\mu^N\in \MNreg(Q_{1,1})$ of $N$-regular measures converging weakly to $\mu$, with $\limsup_{N\to+\infty} I(\mu^N)\le I(\mu)$. Fix $R\ge 3$ and let $(u_n^N,A_n^N)$ be as in Proposition \[propubregular\]. Taking a diagonal subsequence (first with $N$, then with $R$) we obtain $$\limsup_{n\to +\infty} \,  \wE(u_n^{N(n)},A_n^{N(n)})\le I(\mu),$$ and $\b_n^{-1}(B_n^{N(n)})_3\weaklim \mu $. From the compactness statement of Proposition \[gammaliminf\] and uniqueness of $m$ one obtains $\b_n^{-1} (B_n^{N(n)})'\weaklim m$ and $\b_n^{-1}(1-|u_n^{N(n)}|^2)\weaklim \mu $.

Acknowledgment {#acknowledgment .unnumbered}
==============

M. G. thanks E. Esselborn, F. Barret and P. Bella for stimulating discussions about optimal transportation, and the FMJH PGMO fundation for partial support through the project COCA. The work of S.C. was partially supported by the Deutsche Forschungsgemeinschaft through the Sonderforschungsbereich 1060 [*“The mathematics of emergent effects”*]{}.

[^1]: Institut für Angewandte Mathematik, Universität Bonn, Germany, email: sergio.conti@uni-bonn.de

[^2]: LJLL, Université Paris Diderot, CNRS, UMR 7598, France email: goldman@math.univ-paris-diderot.fr, part of his research was funded by a Von Humboldt PostDoc fellowship

[^3]: Max Planck Institute for Mathematics in the Sciences, Leipzig, Germany, email: otto@mis.mpg.de

[^4]: Courant Institute, NYU & Institut Universitaire de France & UPMC, Paris, France, email: serfaty@cims.nyu.edu

[^5]: for $p\ge 1$, we analogously define $W_p^p(\mu,\nu):=\min\lt\{ \mu(Q_L)\int_{Q_L\times Q_L} |x-y|^p \, d\Pi(x,y)  \, : \, \Pi_1=\mu, \ \Pi_2=\nu\rt\}$.
