A first-order reduction of the axisymmetric Willmore equation via two first integrals

Abstract

We derive a first-order formulation of the axisymmetric Willmore equation by combining two independent first integrals, yielding the relation

$$\left[\frac{\Psi q^{2}+2q+2C_{1}\rho}{\sqrt{1-\Psi^{2}}}\right]^{2}+(q^{2}-2)^{2}=C_{2},$$

where $C_{1}$ and $C_{2}$ are constants of integration. Here $\rho$ denotes the distance from the axis of revolution, and $q=\rho\Psi^{\prime}-\Psi$ is the signed hyperbolic curvature of the meridian. We have set $\Psi=\sin\psi$, where $\psi$ is the tangent angle of the meridian, and a prime denotes differentiation with respect to $\rho$. The above relation is equivalent to a quartic algebraic equation for $\rho\Psi^{\prime}$ and thus reduces the third-order axisymmetric Willmore equation to a first-order ordinary differential equation on each continuous algebraic branch. This formulation recovers the sphere, axisymmetric minimal surfaces, and the Clifford torus as special cases, while organizing the local solution space according to the pair $(C_{1},C_{2})$. The same hyperbolic formulation also explains how the spontaneous-curvature term destroys the conservation structure underlying the corresponding first integral for axisymmetric Helfrich surfaces.

1 Introduction

The elasticity of membranes and thin shells has long been a central topic in mathematical physics, soft matter physics, and differential geometry. In the nineteenth century, Poisson introduced a functional for the bending energy of a shell [1]:

$$F_{P}=\frac{k_{c}}{2}\int_{M}(2H)^{2}\,\mathrm{d}A,$$ (1)

where $M$ denotes the midsurface of the shell, $H$ is its mean curvature, $\mathrm{d}A$ is the area element on $M$, and $k_{c}$ is the bending rigidity of the shell. Canham later employed the energy functional (1) to investigate the shapes of red blood cells [2]. The corresponding Euler–Lagrange equation is the well-known Willmore equation:

$$\nabla^{2}H+2H(H^{2}-K)=0,$$ (2)

where $K$ is the Gaussian curvature. Willmore conjectured [3] that the Clifford torus, for which the ratio of the major and minor radii is $\sqrt{2}$, gives the global minimum of the functional (1) among embedded tori. This conjecture was proved by Marques and Neves [4].

In the elastic theory of lipid membranes, the above functional is generalized to the Helfrich functional [5]:

$$F_{H}=\frac{k_{c}}{2}\int_{M}\left[(2H+c_{0})^{2}\right]\mathrm{d}A,$$ (3)

where $c_{0}$ is the spontaneous curvature. Including area and volume constraints, Zhong-can and Helfrich obtained the shape equation of lipid vesicles [6, 7]:

$$k_{c}\nabla^{2}(2H)+k_{c}(2H+c_{0})(2H^{2}-2K-c_{0}H)-2\lambda H+p=0$$ (4)

from the variation of the Helfrich functional (3). Here $\lambda$ and $p$ are two Lagrange multipliers. With the Zhong-can–Helfrich equation (4), many configurations of lipid vesicles have been extensively investigated [8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22].

The Zhong-can–Helfrich equation becomes a fourth-order nonlinear partial differential equation (PDE) in the coordinate representation. Under the assumption of axial symmetry, Hu and Ou-Yang [23] transformed the fourth-order PDE into a third-order ordinary differential equation (ODE). Zheng and Liu [24] reduced the third-order ODE to a second-order ODE by finding a first integral. This first integral can also be understood from a Hamiltonian viewpoint [25]. In the special case in which $p$, $\lambda$, $c_{0}$, and the Zheng–Liu first integral vanish, the second-order ODE is integrable and can be transformed into a first-order ODE [27, 28]. The natural question is whether such a first-order reduction remains possible when the Zheng–Liu integral is nonzero.

The Willmore equation (2) is an important special case of Eq. (4) with vanishing $c_{0}$, $\lambda$, and $p$. Langer and Singer [29, 30] showed that the meridian of a Willmore surface of revolution can be interpreted as an elastic curve in the hyperbolic plane. This geometric correspondence was subsequently developed further by Bryant and Griffiths [31] within a variational reduction framework. The central geometric quantity is the geodesic curvature of the associated elastic curve, parametrized by hyperbolic arclength. A fundamental result of this theory is the existence of an energy-type first integral involving the geodesic curvature and its derivative, which we refer to as the Langer–Singer first integral. This correspondence and the associated first integral have been used to classify Willmore surfaces of revolution, derive explicit representations of such surfaces, and study the Willmore flow of tori of revolution [32, 33, 34].

In the present paper, we combine the Zheng–Liu and Langer–Singer first integrals to derive an explicit first-order differential formula for $\Psi(\rho)$. In addition to yielding a direct reduction of order, the two integration constants associated with these first integrals provide a natural basis for the local classification of axisymmetric Willmore surfaces. The remainder of the paper is organized as follows. In section 2, we present the axisymmetric shape equation. In section 3, we revisit the Zheng–Liu and Langer–Singer first integrals and combine them into a first-order ODE. In section 4, we verify the resulting formulation for the sphere and the Clifford torus. In section 5, we make the corresponding local classification of solution branches. In section 6, we examine the implications of the hyperbolic formulation for axisymmetric Helfrich surfaces. The main results are summarized in the final section.

2 Axisymmetric shape equation

Consider an axisymmetric surface generated by rotating a planar profile curve $\gamma$ shown in Fig. 1 around the $z$-axis. The surface can be parametrized as

$$x = \rho\cos\phi$$ (5)

$$y = \rho\sin\phi$$ (6)

$$z = z_{0}+\int^{\rho}_{\rho_{0}}\tan\psi(\rho)\,\mathrm{d}\rho$$ (7)

where $\psi$ is the angle between the tangent to the profile curve and the horizontal direction. Then the mean curvature and the Gaussian curvature can be expressed as

$$H=-(\rho\sin\psi)^{\prime}/2\rho,$$ (8)

$$K=(\sin^{2}\psi)^{\prime}/2\rho,$$ (9)

respectively. Throughout this paper, a prime denotes differentiation with respect to $\rho$. For an axisymmetric scalar function, the Laplace–Beltrami operator is

$$\nabla^{2}=\frac{1}{\rho^{2}}\frac{\partial^{2}}{\partial\phi^{2}}+\frac{\cos\psi}{\rho}\frac{\partial}{\partial\rho}\left(\rho\cos\psi\frac{\partial}{\partial\rho}\right).$$ (10)

Figure 1: Profile of the generated curve.

Substituting Eqs. (8)–(10) into the Zhong-can–Helfrich equation (4), we obtain the shape equation for axisymmetric lipid vesicles [23]:

$$-\frac{\cos\psi}{\rho}\left\{\rho\cos\psi\left[\frac{\left(\rho\sin\psi\right)^{\prime}}{\rho}\right]^{\prime}\right\}^{\prime}-\frac{1}{2}\left[\frac{\left(\rho\sin\psi\right)^{\prime}}{\rho}\right]^{3}$$

$$+\frac{\left(\rho\sin\psi\right)^{\prime}\left(\sin^{2}\psi\right)^{\prime}}{\rho^{2}}-\frac{c_{0}\left(\sin^{2}\psi\right)^{\prime}}{\rho}+\left(\frac{\lambda}{k_{c}}+\frac{c_{0}^{2}}{2}\right)\frac{\left(\rho\sin\psi\right)^{\prime}}{\rho}+\frac{p}{k_{c}}=0.$$ (11)

If $c_{0}=0$, $\lambda=0$, and $p=0$, this equation reduces to the axisymmetric Willmore equation:

$$-\frac{\cos\psi}{\rho}\left\{\rho\cos\psi\left[\frac{\left(\rho\sin\psi\right)^{\prime}}{\rho}\right]^{\prime}\right\}^{\prime}-\frac{1}{2}\left[\frac{\left(\rho\sin\psi\right)^{\prime}}{\rho}\right]^{3}+\frac{\left(\rho\sin\psi\right)^{\prime}\left(\sin^{2}\psi\right)^{\prime}}{\rho^{2}}=0.$$ (12)

3 First integrals

In this section, we will revisit the Zheng–Liu and Langer–Singer first integrals and combine them into a first-order ODE.

3.1 Zheng–Liu first integral

As pointed out by Zheng and Liu, Eq. (11) has a first integral [24], which can be expressed as

$$\frac{\Psi^{3}-\Psi(\rho\Psi^{\prime})^{2}}{2\rho}-\rho(1-\Psi^{2})\left[\frac{(\rho\Psi)^{\prime}}{\rho}\right]^{\prime}-c_{0}\Psi^{2}+\left(\frac{\lambda}{k_{c}}+\frac{c_{0}^{2}}{2}\right)\rho\Psi+\frac{p\rho^{2}}{2k_{c}}=C_{1},$$ (13)

where $C_{1}$ is a constant of integration. Note that in the Zheng–Liu first integral (13), a new notation

$$\Psi\equiv\sin\psi$$ (14)

has been introduced for simplicity. The constant $C_{1}$ is related to the Noether charge associated with axial symmetry [26].

Putting $c_{0}=0$, $\lambda=0$, and $p=0$ in Eq. (13), we obtain the first integral for the axisymmetric Willmore equation:

$$\frac{\Psi^{3}-\Psi(\rho\Psi^{\prime})^{2}}{2\rho}-\rho(1-\Psi^{2})\left[\frac{(\rho\Psi)^{\prime}}{\rho}\right]^{\prime}=C_{1},$$ (15)

A direct rearrangement transforms this equation into an equivalent compact form

$$\frac{(1-\Psi^{2})^{3/2}}{2\rho\Psi^{\prime}}\left[\frac{\Psi^{2}-(\rho\Psi^{\prime})^{2}}{\sqrt{1-\Psi^{2}}}\right]^{\prime}=C_{1}.$$ (16)

If $C_{1}=0$, this equation implies another first integral

$$\frac{\Psi^{2}-(\rho\Psi^{\prime})^{2}}{\sqrt{1-\Psi^{2}}}=I$$ (17)

with a constant $I$. This result was obtained in previous work [27, 28]. This raises the question of whether a second first integral remains available for the axisymmetric Willmore equation (12) when $C_{1}\neq 0$?

3.2 Langer–Singer first integral

The Langer–Singer correspondence [29, 30] between an axisymmetric Willmore surface and a hyperbolic elastic curve provide another first integral. The most relevant point can be translated into the notation used in the present work as follows. The axisymmetric Willmore equation (12) can be rewritten as

$$\rho\sqrt{1-\Psi^{2}}[\rho^{2}\sqrt{1-\Psi^{2}}\,\Psi^{\prime\prime}]^{\prime}+\frac{1}{2}(\rho\Psi^{\prime}-\Psi)^{3}-(\rho\Psi^{\prime}-\Psi)=0.$$ (18)

Using $\rho\Psi^{\prime\prime}=(\rho\Psi^{\prime}-\Psi)^{\prime}$, this equation can be further rewritten as

$$\rho\sqrt{1-\Psi^{2}}[\rho\sqrt{1-\Psi^{2}}\,(\rho\Psi^{\prime}-\Psi)^{\prime}]^{\prime}+\frac{1}{2}(\rho\Psi^{\prime}-\Psi)^{3}-(\rho\Psi^{\prime}-\Psi)=0.$$ (19)

Equivalently, introducing the differential operator

$$D\equiv\rho\sqrt{1-\Psi^{2}}\,\frac{\mathrm{d}}{\mathrm{d}\rho},$$ (20)

and the signed hyperbolic curvature

$$q=\rho\Psi^{\prime}-\Psi,$$ (21)

Eq. (19) is transformed into the form

$$D^{2}q+\frac{1}{2}q^{3}-q=0.$$ (22)

Multiplying the above equation by $Dq$ then we obtain an energy-type first integral:

$$\frac{1}{2}(Dq)^{2}+\frac{1}{8}(q^{2}-2)^{2}=\frac{C_{2}}{4},$$ (23)

which can be further expressed in an equivalent form:

$$4[\rho\sqrt{1-\Psi^{2}}\,(\rho\Psi^{\prime}-\Psi)^{\prime}]^{2}+[(\rho\Psi^{\prime}-\Psi)^{2}-2]^{2}=C_{2},$$ (24)

We refer to Eq. (23) or (24) as the Langer–Singer first integral.

3.3 First-order form of the axisymmetric Willmore equation

The two first integrals stem from distinct geometric origins. Specifically, the Zheng–Liu first integral corresponds to the Noether charge associated with axial symmetry, whereas the Langer–Singer first integral originates from the conformal invariance of the energy functional (1). Hence, these two integrals are generally independent of each other, which naturally allows us to explore the physical consequences of their combination.

Rewriting the Zheng–Liu integral (15) in terms of $q$ and $Dq$ yields

$$2\sqrt{1-\Psi^{2}}\,Dq+\Psi q^{2}+2q+2C_{1}\rho=0.$$ (25)

Eliminating $Dq$ between Eqs. (24) and (25) casts into a sum of perfect squares:

$$\left[\frac{\Psi q^{2}+2q+2C_{1}\rho}{\sqrt{1-\Psi^{2}}}\right]^{2}+(q^{2}-2)^{2}=C_{2}.$$ (26)

Substituting $q=\rho\Psi^{\prime}-\Psi$ into this relation yields the first-order ODE governing axisymmetric Willmore surfaces:

$$\left[\frac{\Psi(\rho\Psi^{\prime}-\Psi)^{2}+2(\rho\Psi^{\prime}-\Psi)+2C_{1}\rho}{\sqrt{1-\Psi^{2}}}\right]^{2}+[(\rho\Psi^{\prime}-\Psi)^{2}-2]^{2}=C_{2}.$$ (27)

Expanding Eq. (27) yields the following quartic polynomial equation:

$$(\rho\Psi^{\prime})^{4}+\alpha_{2}(\rho\Psi^{\prime})^{2}+\alpha_{1}(\rho\Psi^{\prime})+\alpha_{0}=0,$$ (28)

where the variable coefficients are given by

$$\alpha_{2} =4C_{1}\rho\Psi-2\Psi^{2},$$ (29)

$$\alpha_{1} =8C_{1}\rho(1-\Psi^{2}),$$ (30)

$$\alpha_{0} =\Psi^{4}-4C_{1}\rho(\Psi^{3}+\Psi-C_{1}\rho)-(C_{2}-4)(1-\Psi^{2}).$$ (31)

Along any admissible real-valued root branch of this quartic equation, we can formally express the relation as

$$\rho\Psi^{\prime}=G(\rho,\Psi;C_{1},C_{2}).$$ (32)

Subsequent integration of this first-order ODE yields the meridian shape profile $\Psi=\Psi(\rho;C_{1},C_{2})$, where the choice of the root branch must preserve continuity and remain strictly bounded by $|\Psi|\leq 1$. For a nonzero $C_{1}$, the profile height $z$ can be explicitly reconstructed by integrating Eq. (16) once more [25]:

$$z=C_{0}+\frac{\rho\Psi}{\sqrt{1-\Psi^{2}}}-\frac{\Psi^{2}-G^{2}}{2C_{1}\sqrt{1-\Psi^{2}}},$$ (33)

where the arguments $(\rho;C_{1},C_{2})$ of $\Psi$ and $G$ are omitted for brevity.

4 Two simple examples

4.1 Spherical surface

Consider a spherical surface with radius $R$, which can be expressed as $\Psi=\rho/R$. Eqs. (15) and (24) are simultaneously satisfied with $C_{1}=0$ and $C_{2}=4$ when $\Psi=\rho/R$.

Now we check the consistency of Eq. (28). If $C_{1}=0$ and $C_{2}=4$, we have $\alpha_{2}=-2\Psi^{2}$, $\alpha_{1}=0$, and $\alpha_{0}=\Psi^{4}$ from Eqs. (29)–(31). Then Eq. (28) reduces to $(\rho\Psi^{\prime})^{4}-2\Psi^{2}(\rho\Psi^{\prime})^{2}+\Psi^{4}=0$, which implies $\rho\Psi^{\prime}=\pm\Psi$. The solution $\Psi=\rho/R$ corresponds to the positive branch $\rho\Psi^{\prime}=\Psi$. The negative branch corresponds to an axisymmetric minimal surface, $\Psi=A/\rho$, where $A$ is a constant.

4.2 Clifford torus

The Clifford torus is a classical Willmore surface. It may be represented as a torus of revolution whose generating circle has major radius $R$ and minor radius $r$ with ratio

$$\frac{R}{r}=\sqrt{2}.$$ (34)

This torus can be expressed as $\Psi=\rho/r-\sqrt{2}$. In this case, Eqs. (15) and (24) are simultaneously satisfied with $C_{1}=-1/r$ and $C_{2}=0$.

Conversely, since Eq. (27) is a sum of perfect squares, we can derive

$$\Psi(\rho\Psi^{\prime}-\Psi)^{2}+2(\rho\Psi^{\prime}-\Psi)-2\rho/r = 0$$ (35)

$$(\rho\Psi^{\prime}-\Psi)^{2}-2 = 0$$ (36)

when $C_{1}=-1/r$ and $C_{2}=0$. The above equations lead to $\Psi=\rho/r-\sqrt{2}$ (the other solution $\Psi=\rho/r+\sqrt{2}$ is omitted since it exceeds 1).

5 Local classification of axisymmetric Willmore surfaces

Because $C_{1}$ enters Eq. (27) only through the product $C_{1}\rho$, every nonzero $C_{1}$ can be normalized to $C_{1}=\pm 1$ by a rescaling of $\rho$. Thus, up to scale, the local branches fall into the sectors $C_{1}=0,+1,-1$ which correspond to zero-, positive- or negative Noether charge, respectively. The second constant $C_{2}$ is nonnegative and bounds the signed hyperbolic curvature $q$. Indeed, Eq. (26) implies

$$\max\!\left(0,2-\sqrt{C_{2}}\right)\leq q^{2}\leq 2+\sqrt{C_{2}}.$$ (37)

Consequently, $C_{2}=0$ forces $q=\pm\sqrt{2}$; for $0<C_{2}<4$, $q$ cannot vanish; $C_{2}=4$ is the threshold at which $q=0$ becomes possible; and $C_{2}>4$ allows sign changes of $q$.

The zero-charge sector requires separate attention. Setting $C_{1}=0$ in Eq. (27) gives

$$(\rho^{2}\Psi^{\prime 2}-\Psi^{2})^{2}=(C_{2}-4)(1-\Psi^{2}).$$ (38)

For a regular local branch with $|\Psi|<1$, Eq. (38) requires $C_{2}\geq 4$. Hence there are no regular zero-charge branches with $0\leq C_{2}<4$. At $C_{2}=4$, the elementary factors $\rho\Psi^{\prime}=\pm\Psi$ yield the spherical branch $\Psi=\rho/R$ and the axisymmetric minimal-surface branch $\Psi=A/\rho$.

For $C_{1}=-1$ and $C_{2}=0$, the admissible constant-curvature choice is $q=\sqrt{2}$, giving $\Psi=\rho-\sqrt{2}$, i.e. the Clifford torus after restoring the scale. The branch $q=-\sqrt{2}$ gives $\Psi=\rho+\sqrt{2}>1$ and is inadmissible. The $C_{1}=+1$ sector is obtained by the corresponding orientation-reversed branch. Table 1 summarizes all possible classifications of axisymmetric Willmore surfaces.

Table 1: Classification of axisymmetric Willmore surfaces according to Eq. (27)

$C_{1}$	$C_{2}$	Behavior of $q=\rho\Psi^{\prime}-\Psi$	Representative information
$0$	$0\leq C_{2}<4$	No regular branch with $|\Psi|<1$	Excluded by Eq. (38).
	$4$	$q$ may reach zero; $0\leq q^{2}\leq 4$	Includes $\Psi=\rho/R$ and $\Psi=A/\rho$.
	$>4$	$q$ may change sign	More general zero-charge Willmore branches.
$+1$	$0$	$q=-\sqrt{2}$	Orientation-reversed normalized Clifford-torus branch $\Psi=\sqrt{2}-\rho$; the opposite sign is inadmissible.
	$0<C_{2}<4$	$q$ has fixed nonzero sign	Positive-charge branches with fixed-sign hyperbolic curvature.
	$\geq 4$	$q=0$ is possible; sign changes occur for $C_{2}>4$	Critical and sign-changing positive-charge branches.
$-1$	$0$	$q=\sqrt{2}$	Normalized Clifford-torus branch $\Psi=\rho-\sqrt{2}$; the opposite sign is inadmissible.
	$0<C_{2}<4$	$q$ has fixed nonzero sign	Negative-charge branches with fixed-sign hyperbolic curvature.
	$\geq 4$	$q=0$ is possible; sign changes occur for $C_{2}>4$	Critical and sign-changing negative-charge branches.

6 Implications for axisymmetric Helfrich surfaces

We finally examine how the preceding hyperbolic formulation is modified when the spontaneous curvature $c_{0}$ is nonzero in the Helfrich functional (3). Here, an unconstrained stationary configuration of this functional is referred to as the Helfrich surface. Setting $\lambda=0$ and $p=0$ in Eq. (4) yields the governing shape equation for the Helfrich surface:

$$\nabla^{2}(2H)+(2H+c_{0})(2H^{2}-2K-c_{0}H)=0.$$ (39)

In the axisymmetric case, setting $\lambda=p=0$ in Eq. (13), we have the first integral:

$$\frac{\Psi^{3}-\Psi(\rho\Psi^{\prime})^{2}}{2\rho}-\rho(1-\Psi^{2})\left[\frac{(\rho\Psi)^{\prime}}{\rho}\right]^{\prime}-c_{0}\Psi^{2}+\frac{c_{0}^{2}\rho\Psi}{2}=C_{1}.$$ (40)

Under the Langer–Singer correspondence, the meridian of an axisymmetric Willmore surface is interpreted as an elastic curve in the hyperbolic plane [29, 30], from which the first Langer–Singer integral (24) naturally emerges. Nevertheless, this elegant geometric picture is disrupted for axisymmetric Helfrich surfaces, as the $c_{0}$ term in the Helfrich functional spoils the underlying conformal symmetry. This motivates us to apply the Langer–Singer framework to further investigate the Helfrich functional.

The hyperbolic plane $\{(z,\rho)\}$ with $\rho\geq 0$, equipped with metric

$$\mathrm{d}s_{H}^{2}=\frac{\mathrm{d}\rho^{2}+\mathrm{d}z^{2}}{\rho^{2}}=\frac{\mathrm{d}\rho^{2}}{{1-\Psi^{2}}},$$ (41)

has negative Gaussian curvature $-1$. Along the meridian, the hyperbolic arclengh element is

$$\mathrm{d}s_{H}=\frac{\mathrm{d}\rho}{\rho\sqrt{1-\Psi^{2}}}.$$ (42)

Considering the signed hyperbolic geodesic curvature (21), we have $(\rho\Psi)^{\prime}=\rho\Psi^{\prime}+\Psi=q+2\Psi$. Thus the mean curvature (8) can be expressed as

$$H=-\frac{q+2\Psi}{2\rho}.$$ (43)

In addition, we can calculate the area element, which reads

$$\mathrm{d}A=\frac{\rho\mathrm{d}\rho\mathrm{d}\phi}{\sqrt{1-\Psi^{2}}}=2\pi\rho^{2}\,\mathrm{d}s_{H}.$$ (44)

Substituting the above two equations into (3), we arrive in the hyperbolic representation of the axisymmetric Helfrich functional:

$$F_{H}=\pi k_{c}\int_{\gamma}(q+2\Psi-c_{0}\rho)^{2}\,\mathrm{d}s_{H},$$ (45)

where $\gamma$ corresponds to the meridian.

Note that $(q+2\Psi-c_{0}\rho)^{2}-(q-c_{0}\rho)^{2}=4\Psi(q+\Psi-c_{0}\rho)=4\rho\Psi(\Psi^{\prime}-c_{0})$, and $\rho\Psi(\Psi^{\prime}-c_{0})\mathrm{d}s_{H}$ can be expressed in a total differential form as $\rho\Psi(\Psi^{\prime}-c_{0})\,\mathrm{d}s_{H}=-\mathrm{d}(\sqrt{1-\Psi^{2}}+c_{0}z)$. Since $\int_{\gamma}\mathrm{d}(\sqrt{1-\Psi^{2}}+c_{0}z)$ will enter the boundary term, the functional (45) is variationally equivalent to

$$\tilde{F}_{H}=\pi k_{c}\int_{\gamma}(q-c_{0}\rho)^{2}\,\mathrm{d}s_{H},$$ (46)

which implies that the Helfrich meridian is an inhomogeneous hyperbolic elastica with spontaneous curvature $q_{0}=c_{0}\rho$ which is proportional to $\rho$. If $c_{0}=0$, the above functional reduces to the Langer–Singer form. Its variation gives the axisymmetric Willmore equation (22), from which one obtains the Langer–Singer first integral (24). This conservation law ultimately reflects the conformal invariance of the functional (1). The presence of $c_{0}$ breaks the conformal invariance and the functional (46) gives the axisymmetric Helfrich equation:

$$D^{2}q+\frac{1}{2}q^{3}-q+2c_{0}\rho\Psi(q+\Psi)-\frac{c_{0}^{2}\rho^{2}}{2}(q+2\Psi)=0,$$ (47)

which is equivalent to (11) with vanishing $\lambda$ and $p$. The explicit $\rho$- and $\Psi$-dependent terms related to $c_{0}$ prevent Eq. (47) from having the same autonomous energy integral as the Willmore equation (24).

Although the spontaneous-curvature terms prevent the immediate construction of the Langer–Singer integral, the hyperbolic representation still identifies two simple exact branches. First, consider the pointwise zero-density condition for the functional (45) which corresponds to

$$q+2\Psi-c_{0}\rho=0.$$ (48)

Combining this condition with Eq. (21) gives

$$\Psi=\frac{c_{0}\rho}{2}+\frac{A}{\rho},$$ (49)

where $A$ is a constant. This branch has constant mean curvature $H=-c_{0}/2$.

Second, consider the pointwise zero-density condition for the variationally equivalent functional (46) which corresponds to

$$q-c_{0}\rho=0.$$ (50)

Combining this condition with Eq. (21) gives

$$\Psi=c_{0}\rho\ln\left(\frac{\rho}{\rho_{B}}\right),$$ (51)

where $\rho_{B}>0$ is a constant. This solution was previously used to describe a biconcave, disk-like red-blood-cell profile [9].

7 Conclusion

In summary, the combination of the Zheng–Liu and Langer–Singer first integrals yields a first-order differential formulation (27), equivalent to the quartic equation (28), which fully characterizes axisymmetric Willmore surfaces. This framework provides a unified classification of the local solution space parametrized by $(C_{1},C_{2})$ while seamlessly recovering classic exact solutions. When the spontaneous curvature $c_{0}$ is introduced, the Helfrich functional corresponds to an inhomogeneous hyperbolic elastica with preferred curvature $c_{0}\rho$, which breaks the conservation law underlying the Langer–Singer first integral. Nevertheless, this correspondence still reveals the constant-mean-curvature branch and the logarithmic branch associated with biconcave profiles. To further explore the integrability of Eq. (40) for axisymmetric Helfrich surfaces, a promising avenue is to derive additional conservation laws via the Ibragimov–Noether theorem [35].

National Natural Science Foundation of China (Grant Nos. 12475032, 11975050, and 11322543)