Existence and Regularity of Spheres Minimising the Canham-Helfrich Energy

We prove existence and regularity of minimisers for the Canham-Helfrich energy in the class of weak (possibly branched and bubbled) immersions of the $2$-sphere. This solves (the spherical case) of the minimisation problem proposed by Helfrich in 1973, modelling lipid bilayer membranes. On the way to prove the main results we establish the lower semicontinuity of the Canham-Helfrich energy under weak convergence of (possibly branched and bubbled) weak immersions.


Introduction
The basic structural and functional unit of all known living organisms is the cell. The interior material of a cell, the cytoplasm, is enclosed by biological membranes. Most of the cell membranes of living organisms are made of a lipid bilayer, which is a thin polar membrane consisting of two opposite oriented layers of lipid molecules.
In 1970, in order to explain the biconcave shape of red blood cells, Canham [Can70] proposed a bending energy density dependent on the squared mean curvature.
Three years later Helfrich proposed the following curvature elastic energy per unit area of a closed lipid bilayer [Hel73,Equation (12)] where H is the mean curvature, K is the Gauss curvature, c 0 is the so-called spontaneous curvature, and k c ,k c are the curvature elastic moduli. The values of the parameters can be measured experimentally (e.g. see [EF72], [DH76] for c 0 , and [MH90] for k c ). The constantk c is not important for the purpose of this paper as by the Gauss-Bonnet Theorem, the integrated Gauss curvature is a topological constant. Lipid bilayers are very thin compared to their lateral dimensions, thus are usually modelled as surfaces. Suppose the surface and hence the membrane is represented by a 1.1 Remark. A candidate embedding Φ 0 which achieves the global minimum is called a minimiser. In general it is not unique and, more dramatically, it may not exist: later in the introduction we show that for a suitable choice of parameters the minimum is achieved by a singular immersion and it cannot be achieved by a smooth one. The constraints and the functional H c 0 are invariant under reparametrisation as well as rigid motions in R 3 . Of course, in order to have a non-empty class of competitors, the constraints have to satisfy the Euclidean isoperimetric inequality A 3 0 36πV 2 0 . 1.2 Problem. Let c 0 , A 0 , and V 0 be given constants. Minimise H c 0 ( Φ) in the class of smooth embeddings Φ : S 2 → R 3 subject to the constraints Area Φ = A 0 and Vol Φ = V 0 . (1.3) That is, find an embedding Φ 0 : S 2 → R 3 such that Area Φ 0 = A 0 , Vol Φ 0 = V 0 , and for any other smooth embedding Φ : S 2 → R 3 satisfying the constraints (1.3).
Problem 1.2 is the classical formulation suggested in [Hel73] and [DH76]. According to Bernard, Wheeler and Wheeler [BWW17], many issues for the Canham-Helfrich energy, including Problem 1.2, remain open and form important questions that future research should address. A similar problem in the 2-dimensional case (i.e. closed curves in the Euclidean plane) was formulated and solved by Bellettini, Dal Maso, and Paolini [BDMP93] by a relaxation procedure.
While there was essentially no work on the variational theory of the Willmore energy after Blaschke's seminal work, in 1965 Willmore [Wil65] reintroduced this Lagrangian which is now named after him. He showed that the round sphere is a minimiser of Problem 1.2 in the special case c 0 = 0 without constraints (1.3), see [Wil82]. Simon [Sim86] proved existence of higher genus minimisers for the Willmore energy (see also Kusner [Kus96] and Bauer-Kuwert [BK03]), using the so-called ambient approach, i.e. convergence of surfaces is considered in the measure-theoretic sense. Rivière [Riv08,Riv14] proved the analogous result with the so called parametric approach, i.e. based on PDE theory and functional analysis as opposed to geometric measure theory. In the present paper we shall adapt the parametric approach. The case c 0 = 0 with constraints (1.3) was solved by Schygulla [Sch12] using the ambient approach, and generalised to higher genus surfaces by Keller, the first author, and Rivière [KMR14] using the parametric approach.
From the mathematical point of view, the spontaneous curvature c 0 causes a couple of differences between the Willmore energy and the Canham-Helfrich energy. Most obviously, the Canham-Helfrich energy cannot be bounded below by a strictly positive constant, whereas the Willmore energy is bounded below by 4π, see (2.13). Secondly, while the Willmore functional is invariant under conformal transformations, the Canham-Helfrich energy is not conformally invariant. We will be concerned with yet another property that fails for the Canham-Helfrich energy due to non-negative spontaneous curvature. Namely lower semi-continuity with respect to varifold convergence: while it is well known that the Willmore functional is lower semi-continuous under varifold convergence, the Canham-Helfrich energy in general is not. Indeed, Große-Brauckmann [GB93, Remark (ii) on page 550] constructed a sequence of non-compact infinite genus surfaces Σ 1 , Σ 2 , . . . with constant mean curvature equal to 1 which converges in the varifold sense to a double plane Σ ∞ . Hence, the mean curvature H ∞ of the limit Σ ∞ is zero and for any continuous non-negative, non-zero function η on R 3 of compact support, where H 2 is the 2-dimensional Hausdorff measure. Hence, the general Canham-Helfrich energy is not lower semi-continuous under varifold convergence. However, in order to solve Problem 1.2 by the so-called direct method of calculus of variations, lower semi continuity is required. According to Röger [MR208], it was an open question under which conditions/in which natural weak topology on the space of immersions one obtains lower semi continuity of the Canham-Helfrich energy. This is presumably the reason why Problem 1.2 was only partially solved for non-zero spontaneous curvature c 0 . The axisymmetric case was solved in 2013 by Choksi and Veneroni [CV13], who proved the existence of a minimiser of the Canham-Helfrich energy among a suitable class of axisymmetric (possibly singular) surfaces under fixed surface area and enclosed volume constraints. Five years later, Dalphin [Dal18] showed existence of minimisers in a class of C 1,1 surfaces whose principal curvatures are bounded by a given constant 1/ε. Though, in his setting, it is still unclear how to get compactness and lower semi-continuity as ε tends to zero.
As already alluded to, we tackle Problem 1.2 by the direct method of calculus of variations. The issue is of course to find a suitable class F A 0 ,V 0 of admissible maps (endowed with a suitable topology) having area A 0 and enclosed volume V 0 such that the Canham-Helfrich energy is lower semi-continuous and has (pre-)compact sub-levels. A natural choice is the class of weak (Sobolev) immersions in W 2,2 (S 2 , R 3 ), already employed in the context of the Willmore energy for instance by Rivière [Riv14] or Kuwert and Li [KL12]. We will use the space of bubble trees of weak possibly branched immersions. It will shortly become clear why we have to allow branched points and multiple bubbles.
In the following, we denote by · (resp. ×) the Euclidean scalar (resp. vector) product on R 3 .

Definition.
A map Φ : S 2 → R 3 is called weak (possibly branched) immersion with finite total curvature if Φ ∈ W 1,∞ (S 2 , R 3 ) and the following holds: 1. There exists C > 1 such that, for a.e. p ∈ S 2 , where the norms are taken with respect to the standard metric on S 2 and with respect to the Euclidean metric of R 3 , and where d Φ × d Φ is the tensor given in local coordinates on S 2 by 2. There exist a positive integer N and finitely many points b 1 , . . . , b N ∈ S 2 such that log |d Φ| ∈ L ∞ loc (S 2 \ {b 1 , · · · , b N }); 3. The Gauss map n, defined by in any local chart x of S 2 , satisfies n ∈ W 1,2 (S 2 , R 3 ).
The space of weak (possibly branched) immersions with finite total curvature is denoted by F.
Since by assumption Φ is a Lipschitz map, it induces an L ∞ -metric g given by for elements X, Y of the tangent bundle T S 2 . In the usual way (see for instance [Heb99, 1.2]), the L ∞ -metric g induces a Radon measure µ on S 2 which is mutually absolutely continuous to the 2-dimensional Hausdorff measure on S 2 . Using Müller-Svěrák theory of weak isothermic charts [MŠ95] and Hélein's moving frame technique [Hél02] one can prove the following proposition (see for instance [Riv16]) Then there exists a bilipschitz homeomorphism Ψ of S 2 such that Φ • Ψ is weakly conformal: it satisfies almost everywhere on S 2 where x is a local arbitrary conformal chart on S 2 for the standard metric. Moreover 1.5 Remark. In view of Proposition 1.4, a careful reader could wonder why we do not work with conformal W 2,2 weak, possibly branched, immersions only and why we do not impose for the membership in F, Φ to be conformal from the beginning. The reason why it is technically convenient not to impose conformality from the beginning is to allow general perturbations in the variational problem, which do not have to respect infinitesimally the conformal condition.
The reason why we chose the class F as above is the following theorem of the first author and Rivière [MR14, Theorem 1.5] (see also [CL14] where n k are the Gauss maps, µ k are the corresponding Radon measures, and Then, after passing to a subsequence, there exist a family Ψ k of bilipschitz homeomorphisms of S 2 , a positive integer N , sequences f 1 k , . . . , f N k of positive conformal diffeomorphisms of S 2 , ξ 1 ∞ , . . . , ξ N ∞ ∈ F, non-negative integers N 1 , . . . , N N , and finitely many points on the sphere The theorem already gives (pre-)compactness, a notion of convergence, and lower semi-continuity (actually, continuity) of the third summand in (1.2) of the Canham-Helfrich energy, i.e. of the area functional. Indeed, T := ( f ∞ , ξ 1 ∞ , . . . , ξ N ∞ ) forms a bubble tree, see Definition 3.2. In particular, the limit T is not in the class F anymore. At an informal level, a non expert reader can think of a bubble tree T := ( f , ξ 1 , . . . , ξ N ) as a "pearl necklace" where each "pearl" corresponds to the image of a possibly branched weak immersion ξ i (S 2 ) and f is a Lipschitz map from S 2 to R 3 "parametrising" the whole pearl necklace, in particular f (S 2 ) = N i=1 ξ i (S 2 ). To get a better understanding of why we obtain a bubble tree in the limit, we will look at an example of Problem 1.2. Let Then, the infimum in Problem 1.2 is achieved by the bubble tree T = ( f , Id S 2 , Id S 2 ) of twice the unit sphere. Indeed, H c 0 ( T ) = 0 and H c 0 ( Φ) ≥ 0 for any other smooth immersion Φ of S 2 into R 3 , so T achieves the infimum. A minimising sequence Φ k (S 2 ) of smoothly embedded spheres converging to such a bubble tree can be achieved by glueig (1 + 1/k)S 2 to (1 − 1/k)S 2 via a small catenoidal neck of size 2/k. Notice also that if Φ satisfies H c 0 ( Φ) = 0, then the image Φ[S 2 ] is the unit sphere by a classical theorem of Hopf [Hop83]. Getting a bubble tree in the limit is in accordance with the earlier result on existence of minimisers by Choksi and Veneroni [CV13] in the axisymmetric case: indeed the minimiser in [CV13, Theorem 1] is made by a finite union of axisymmetric surfaces. Moreover, the bubbling phenomenon is also known as budding transition in biology and has been recorded with video microscopy, see Seifert [Sei97] or Seifert, Berndl, and Lipowsky [SBL91].
In Chapter 3 we sharpen Theorem 1.6 in a way that we get lower semi-continuity for the Canham-Helfrich functional. This can be seen as a possible answer to the aforementioned open question raised by Röger [MR208]. In Chapter 4 we compute the Euler-Lagrange equation for the Canham-Helfrich energy in divergence form. Moreover, we prove that all the weak branched conformal immersions of a minimising bubble tree (actually more generally for a critical bubble tree) are smooth away from their branch points. Our proof is based on the regularity theory for Willmore surfaces developed by Rivière [Riv08]. It relies on conservation laws discovered by Rivière [Riv08] in the context of the Willmore energy and adjusted by Bernard [Ber16] for the Canham-Helfrich energy. We get the following final result.
Moreover, for each i ∈ {1, . . . , N } there exist a non-negative integer N i and finitely . . be a minimising sequence of (1.6). There holds where H Φ k is the mean curvature corresponding to Φ k , see (2.3). By the Gauss-Bonnet theorem (see (2.6) for the precise statement in case of weak branched immersions and (2.8) for the estimate below), which means the first inequality of (3.17) is satisfied. Moreover, (2.10) implies the second inequality of (3.17). Hence, we can apply Theorem 3.3, Theorem 4.3, Lemma 4.4, (2.13) and (2.12) to conclude the proof.
1.8 Remark. The arguments in the proof of Theorem 1.7 yield also that the minimum of H c 0 is achieved in the class of bubble trees of possibly branched weak immersions, by a bubble tree of possibly branched immersions which are smooth out of the branch points.
A more general form of the Canham-Helfrich energy is given by for Φ ∈ F where the parameter α 0 is referred to as tensile stress, and ρ 0 as osmotic pressure. We get the following solution of Problem (P2) from the introduction in [BWW17].
Moreover, if the inequality is strict, then there exist Φ 0 ∈ F, a positive integer N 0 , and As α > 0, we have sup k Area Φ k < ∞ and thus, using (1.7), also sup k Therefore, analogously to the proof of Theorem 1.7, we can apply Theorem 3.3 to obtain an integer N and Φ 1 , . . . , Obviously, and since there are no constraints, we simply get N = 1.

Acknowledgements.
A.M. is supported by the EPSRC First Grant EP/R004730/1 "Optimal transport and Geometric Analysis" and by the ERC Starting Grant 802689 "CURVATURE". C.S. is supported by the EPSRC as part of the MASDOC DTC at the University of Warwick, grant No. EP/HO23364/1.

Notation
We adopt the conventions of [Riv16]. To avoid indices and to get clearly arranged equations, we will employ the following suggestive notation. For R 3 valued maps e and f defined on the unit disk D 2 , we write where · denotes the Euclidean inner product and × denotes the usual vector product on R 3 . Similarly, for λ : Moreover, for a vector field The m-dimensional Lebesgue measure is denoted by L m .

Weak (possibly branched) conformal immersions
We adapt the notion of weak immersions which was independently formalized by Rivière [Riv14] and Kuwert and Li [KL12]. Let (Σ, c 0 ) be a smooth closed Riemann surface (in the rest of the paper we will take (Σ, c 0 ) to be the 2-sphere endowed with the standard round metric). Without loss of generality we can assume that (Σ, c 0 ) is endowed with a metric g c 0 of constant curvature and area 4π (see for instance [Jos06]). For the definition of the Sobolev spaces W k,p (Σ, R 3 ) on Σ see for instance Hebey [Heb99]. A map Φ : Σ → R 3 is called a weak branched conformal immersion with finite total curvature if and only if there exists a positive integer N , finitely many points almost everywhere for any conformal chart x of Σ, and its Gauss map n defined by The space of weak branched conformal immersions with finite total curvature is denoted by F Σ or just F in case Σ = S 2 . We define the L ∞ -metric g pointwise for almost every p ∈ Σ by for elements X, Y of the tangent space T p Σ. In the usual way, the L ∞ -metric g induces a Radon measure µ g on Σ. The conformality condition (2.1) implies that g = e 2λ g c 0 for some λ ∈ L ∞ loc (Σ \ {b 1 , . . . , b N }) called conformal factor. Moreover, we define the second fundamental form I pointwise for almost every p ∈ Σ by Then Φ ∈ W 2,2 (Σ) and the conformal factor λ is an element of L 1 (Σ). Moreover, for each singular point b j , j = 1, . . . , N, there exists a strictly positive integer n j ∈ N such that the following holds: for some ω ∈ C 0 ∩ W 1,2 .
• The multiplicity of the immersion Φ at Φ(b j ) is n j . Moreover, if n j = 1, then Φ is a conformal immersion of a neighbourhood of b j .
• The conformal factor λ satisfies the following singular Liouville equation in distributional sense where δ b j is the Dirac delta centred at b j , K Φ is the Gaussian curvature of Φ, and K 0 ∈ R is the (constant) curvature of (Σ, g c 0 ).
By integrating the singular Liouville equation (2.5), we obtain the Gauss-Bonnet Theorem for weak branched immersions: where χ(Σ) is the Euler Characteristic of Σ. Note in particular that, once the topology of Σ is fixed, the number of branch points counted with multiplicity is bounded by the Willmore energy: Moreover, the Willmore energy controls the L 2 norm squared of the second fundamental form: (2.8)

Simon's monotonicity formula and Li-Yau inequality for weak branched immersions
Let Φ ∈ F Σ be any weak branched conformal immersion with finite total curvature and branch points {b 1 , . . . , b N }. In the usual way (by splitting the vector field in its tangential and normal parts and using integration by parts) one shows A simple cut-off argument together with (2.4) shows that the first variation formula (2.9) is true for all X ∈ W 1,2 (Σ, R 3 ). In the following we will gather a couple of facts that are well known for weak unbranched immersions and, due to (2.9), are also valid for weak branched conformal immersions with finite total curvature. Firstly, letting X(p) := Φ(p) − Φ(a 0 ) for p ∈ Σ and some fixed a 0 ∈ Σ, one has div Φ X = 2 and hence, see Simon [Sim86, Lemma 1.1] (2.10) The push forward measure µ : where the weak mean curvature is almost everywhere given by , then Φ is an embedding (compare also with Proposition 2.1).

Canham-Helfrich energy
Given real numbers c 0 ∈ R and α, ρ 0 as well as a weak branched conformal immersion with finite total curvature Φ : Σ → R 3 , we define the Canham-Helfrich energy H c 0 α,ρ ( Φ) in its most general form by (2.14) Note that, in case Φ : Σ → R 3 is a smooth (actually Lipschitz is enough) embedding, by the Divergence Theorem the last integral equals the volume enclosed by Φ(Σ). The parameter α is referred to as tensile stress, ρ as osmotic pressure. Compare this definition for instance with [Ber16, Equation (3.6)] or [BWW17].

Existence of minimisers
In this chapter we will prove compactness of sequences with uniformly bounded Willmore energy and area as well as lower semi-continuity of the Canham-Helfrich energy under this convergence, see Theorem 3.3. The proof of Theorem 3.3 will build on top of [MR14] and the next Lemma 3.1 which establishes the convergence of the constraints and the lower semi-continuity of the Willmore energy away from the branch points (Lemma 3.1 should be compared with [Riv16, Lemma 5.2]).

Lemma (Convergence outside the branch points).
Suppose ξ 1 , ξ 2 , . . . ∈ F S 2 is a sequence of weak branched conformal immersions with finite total curvature of the 2sphere S 2 into R 3 , µ 1 , µ 2 , . . . are the corresponding Radon measures on S 2 , n 1 , n 2 , . . . are the corresponding Gauss maps, there exists ξ ∞ ∈ F S 2 , a positive integer N , and b 1 , . . . , b N ∈ S 2 such that Then, there exists a sequence of positive numbers s 1 , s 2 , . . . converging to zero such that where the balls are taken with respect to the geodesic distance on the standard S 2 , µ ∞ and n ∞ are the Radon measure and the Gauss map corresponding to ξ ∞ , and the H k 's and H ∞ are the mean curvatures corresponding to the ξ k 's and ξ ∞ . Equations (3.4)-(3.6) remain valid for s k replaced by any sequence t k converging to zero and satisfying t k s k , for all k ∈ N. Moreover, for any sequence s 1 , s 2 , . . . of positive numbers converging to zero, there exists a sequence t k s k converging to zero such that Proof. Suppose U is an open subset of S 2 \ {b 1 , . . . , b N }, K is a compact subset of U , and x : U → R 2 is a conformal chart for S 2 . Denote by the conformal factors. Notice that the volume element corresponding to ξ k is given by e 2λ k . In a first step we will show that for any 1 p < ∞, as well as A simple argument by contradiction shows that it is enough to prove the statement after passing to a subsequence of k. Since the ξ k 's and ξ ∞ are conformal and x is a conformal chart, we can write the mean curvature vector as where ∆ is the flat Laplacian with respect to x. By Hypothesis (3.3), we have that as k → ∞ weakly in L 2 (x[K], R 3 ), which implies (3.10). By the Rellich-Kondrachov Compactness Theorem, after passing to a subsequence, there holds for any 1 p < ∞. Therefore, using Hypothesis (3.2) and passing to a further subsequence, it follows It follows , which implies (3.11) by lower semi-continuity of the L 2 -norm under weak convergence. Similarly, from Hypothesis (3.2) and the strong convergence (3.12), we infer (3.8).
Again by the strong convergence (3.12) and Hypothesis (3.2), we can extract a subsequence such that by dominated convergence, Using this and the fact that by the Rellich-Kondrachov Compactness Theorem for any 1 p < ∞, one verifies (3.9). Next, let r k be any sequence of positive numbers converging to zero and abbreviate First, notice that for any Borel function f on S 2 with S 2 |f | dµ ∞ < ∞, there holds which is a consequence of the dominated convergence theorem and the fact that finite sets have µ ∞ measure zero. Let n 0 = 1. For each positive integer j, we use (3.8) to inductively choose n j > n j−1 such that Moreover, define l k = j for all integers k with n j−1 < k n j and define s k = r l k . Then, we have that s k → 0 as k → ∞ as well as which in particular remains valid for s k replaced by any t k s k . Hence, by (3.13) we can deduce (3.4). Using the convergence on compact sets (3.9)-(3.11), Equations (3.5)-(3.7) follow similarly. It only remains to show that Equation (3.5) is still valid after replacing s k by any sequence t k s k converging to zero. Hence, we only have to show that

This follows as by Hölder's inequality
The first factor on the right hand side is bounded by (3.1). To see that the second factor goes to zero as k tends to infinity, we apply (3.14) and the fact that In the following we will define the notion of a bubble tree. The idea is that the different bubbles can be parametrised by decomposing a single 2-sphere. The bubbles can then be attached to each other by a Lipschitz map, see (3.15) and (3.16). There exist open geodesic balls B 1 , . . . , B N ⊂ S 2 such that

Definition
For all i ∈ {1, . . . , N } there exists a positive integer N i and disjoint open geodesic balls B i,1 , . . . , B i,N i ⊂ S 2 whose closures are included in B i such that For all i ∈ {1, . . . , N } there exist distinct points b i,1 , . . . , b i,N i ∈ S 2 and a Lipschitz diffeomorphism Moreover, for all i ∈ {1, . . . , N }, and for all j ∈ {1, . . . , N i } there exists p i,j ∈ R 3 such that The next theorem establishes the weak closure of bubble trees, as well as the convergence of the constraints in the Helfrich problem and the lower semi-continuity of the Willmore energy. The proof builds on top of [MR14].

Theorem (Weak closure and lower semi-continuity of bubble trees). Suppose
is a sequence of bubble trees of weak immersions and Then, there exists a subsequence of T k which we again denote by T k such that N k = N for some positive integer N and there exists a sequence of diffeomorphisms Ψ k of S 2 such that  many points b i,j,1 , . . . b i,j,Q for some branched Lipschitz conformal immersion ξ i,j ∞ ∈ F S 2 . Furthermore, is a bubble tree of weak immersions and Proof. We first consider the special case where N k = 1 for all positive integers k. By and for j = 1, . . . , N i the sets of indices Finally, for any µ Φ k integrable Borel function ϕ on S 2 , we get We notice that by the strong convergence (1.5), for some finite number C > 0. Hence, by Hölder's inequality By (3.17) and (3.19), the right hand side of each line goes to zero as k tends to infinity. That means the last term of Equation (3.20) goes to zero as k tends to infinity when ϕ is replaced by H Φ k as well as when ϕ is replaced by n Φ k · Φ k . Therefore, using (3.5) and (3.6), we can conclude the convergence of the integrated mean curvature and the convergence of the volume from (3.20
for all ω ∈ C ∞ c (Σ \ {b 1 , . . . , b N }, R 3 ). In the following, we will first compute the Canham-Helfrich equation in divergence form, see Lemma 4.1. Then, we will prove that a weak immersion satisfying the Canham-Helfrich equation is smooth away from its branch points, see Theorem 4.3. The proof is based on the regularity theory for weak Willmore immersions developed by Rivière [Riv08,Riv16]. An important step in Riviere's regularity theory is the discovery of hidden conservation laws for weak Willmore immersions. In the framework of Canham-Helfich immersions, the corresponding hidden conservation laws were discovered by Bernard [Ber16].
corresponds to the first variation of the Willmore energy.
Proof. After composing with a conformal chart away from the branch points, we may assume that Φ is a map D 2 → R 3 . Let ω ∈ C ∞ c (D, R 3 ) and define Φ t := Φ + t ω for t ∈ R. The conformal factor λ is given by 2e 2λ = |∇ Φ| 2 and the metric coefficients (g t ) ij by Therefore, using Using that ω has compact support in D 2 , and using the symmetry of the second fundamental form, i.e.

Remark.
Usually in the literature (see for instance [Ber16,Chapter 3.3]) one finds the expression of the first variation for Hdµ g written as (4.10) It is not hard to check the equivalence of (4.10) with (4.6) proved above. The advantage of the expression (4.6) is two fold: first it invokes less regularity of the immersion map Φ, second it is already in divergence form. Both advantages will be useful in establishing the regularity of weak Canham-Helfrich immersions: indeed, (4.10) would correspond to an L 1 term in the Euler-Lagrange equation (which is usually a problematic right hand side for elliptic regularity theory) while (4.6) corresponds to the divergence of an L 2 term (which is a much better right hand side in elliptic regularity).

Theorem (Smoothness of weak Canham-Helfrich immersions)
. Suppose Φ ∈ F Σ is a weak Canham-Helfrich immersion. Then Φ is a C ∞ immersion away from the branch points.
Proof. After composing with a conformal chart of Σ away from the branch points onto the unit disk D 2 , we may assume that Φ is a map D 2 → R 3 without branch points and Φ satisfies the Canham-Helfrich equation (4.2). It is enough to show that Φ ∈ C ∞ (B 1/2 (0)). The proof splits into three parts.
Step 1: Conservation laws. In view of the Canham-Helfrich equation (4.2), we define 1 T ∈ L 2 (D 2 , (R 3 ) 2 ) by letting Then, div T = − W where W is as in (4.3). Hence W ∈ H −1 (D 2 , R 3 ) and there exists a . Therefore, we can find After the breakthrough of Rivière [Riv08], Bernard [Ber16, Chapter 2.2] showed that by invariance of the Willmore functional under conformal transformation and the weak Poincaré Lemma, one can find potentials L, R ∈ W 1,2 (D 2 , R 3 ), and S ∈ W 1,2 (D 2 , R) such that Indeed, from [Ber16, Chapter 3.3] we find that R, S, X, Y , and Φ satisfy the following system of conservation laws (4.12) (4.13) (4.14) (4.15) 1 Note that here, T is not a bubble tree.