Stability of hypersurfaces with constant mean curvature trapped between two parallel hyperplanes

Static equilibrium configurations of continua supported by surface tension are given by constant mean curvature (CMC) surfaces which are critical points of a variational problem to extremize the area while keeping the volume fixed. CMC surfaces are used as mathematical models of a variety of continua, such as tiny liquid drops, stars, and nuclei, to play important roles in both mathematics and physics. Therefore, the geometry of CMC surfaces and their properties such as stability are of special importance in differential geometry and in a variety of physical sciences. In this paper we examine the stability of CMC hypersurfaces in arbitrary dimensions, possibly having boundaries on two parallel hyperplanes, by investigating the second variation of the area. We determine the stability of non-uniform liquid bridges or unduloids for the first time in all dimensions and all parameter (the ratio of the neck radius to bulge radius) regimes. The analysis is assisted by numerical computations.


Introduction
A static equilibrium configuration of a continuum supported by surface tension is known to be well approximated by a constant-mean-curvature (CMC) surface, which extremizes the surface area for given volume and boundary conditions. CMC surfaces are used as mathematical models of a variety of continua, such as liquid drops, stars, and nuclei, to play important roles in both mathematics and physics [1,2]. In the three-dimensional Euclidean space R 3 , one of fundamental problems regarding CMC surfaces is to find stable CMC surfaces which possibly have boundaries on given two parallel planes Π 1 and Π 2 . Here, a CMC surface is said to be stable if the second variation of the area for any volume-preserving variation is non-negative. Since only stable surfaces are stably realized in natural phenomena, it is important to judge the stability for a CMC surface, which is difficult in general. Allowing no self-intersections of surface, it is shown that equilibrium surfaces contained in the region bounded by Π 1 and Π 2 are axially-symmetric CMC surfaces with the straight line perpendicular to Π i as its rotation axis [3] (which will be called the z-axis throughout this paper), and they make contact angles π/2 with Π i (i = 1, 2) [4, Sect. 1.6]. Hence these surfaces are spheres, hemispheres, parts of cylinders and unduloids [5] (Fig. 1). Among them, only spheres, hemispheres, and thick cylinders are stable. Thin cylinders and unduloids are unstable CMC surfaces, i.e., they extremize the area but do not minimize it for a given volume [6,7]. 1 The instability of thin cylinder is known as the Plateau-Rayleigh instability [8,9] in fluid mechanics.
An interesting non-trivial aspect of the above problem is that the stability of unduloids depends on the dimension. There is a higher-dimensional counterpart of the unduloid that is an axially-symmetric CMC hypersurface in R n+2 (n = 1, 2, · · · ) which is periodic and has no self-intersection. We call it also an unduloid. A half period of an unduloid from a neck to the next bulge (such as the left figure in Fig. 1) satisfies the boundary condition. We define its non-uniformness parameter as s := 1−(h min /h max ) ∈ (0, 1), where h min and h max denote the radii of the unduloid at the neck and bulge, respectively. Then, for any negative number H and s ∈ (0, 1), up to rigid motion in R n+2 , there exists exactly one unduloid with mean curvature H with respect to the outward-pointing unit normal and having non-uniformness s. Denote by U = U(H, s) a half period of such unduloid. Then, for the Gromov-Hausdorff distance, lim s→1−0 U(H, s) is a hemisphere with radius 1/|H|, and lim s→0+0 U(H, s) is a cylinder with radius r = n/[(n+ 1)|H|] and height L = ( √ nπ)((n+ 1)|H|) −1 . While the unduloids in higher dimensions were numerically obtained and their geometric quantities were computed [12][13][14], their stability has not been clarified completely so far. Let U be a half period of an unduloid. Then the following results on the stability are known [6,7,10,11].
(i) For any n ≥ 1, if U is sufficiently close to a hemisphere, then U is unstable.
(iv) For n ≥ 8, there exists some U that is stable.
In this paper, which corresponds to an extended version of a letter by the present authors [15], we comprehensively examine the stability of unduloids U in all dimensions and parameter regimes by investigating the second variation of area with the help of numerical computations. The results are summarized as statements (I)-(IV) in Sect. 5. A noteworthy result there is, roughly speaking, as follows: When 7 ≤ n ≤ 9, if U is sufficiently close to either a cylinder or a hemisphere, then U is unstable, and moreover there exists stable U.
Especially, the existence of stable unduloid for n = 7 and instability of a half period of an unduloid close to a hemisphere for 7 ≤ n ≤ 9 are found for the first time in this paper.
The geometric quantities of unduloid such as surface area, bulk volume, and mean curvature are obtained with the help of numerical integration. We will see that the stability is determined by the behaviors of these geometric quantities and stability criteria. There, besides the standard criteria for the stability, we use the bifurcation technique (see Sect. 3.1, A.4) developed in [16] in order to judge the stability, which was not used in the papers mentioned above. It is remarkable that the regions of s where the unduloid is stable (resp. unstable) completely coincide with those where the enclosed volume V (s) is non-increasing (resp. increasing) for any n (see Table 1, Sect. 5).
Before starting analysis, let us mention that the higher-dimensional CMC hypersurfaces attract much attention in the study of general relativity, in particular, black holes. The black-hole counterparts of the cylinder and unduloid are called uniform black strings and non-uniform black strings respectively, and they exhibit various similarities with their counterparts [17,18]. Furthermore, the 'surface' of a black hole (i.e., event horizon) was recently shown to indeed be approximated by a timelike CMC hypersurface in a largedimension limit of general relativity [19]. We will return this point in Sect. 5.
The organization of this paper is as follows. We begin by calculating the variations of surface area and bulk volume for axially symmetric hypersurfaces in Sect. 2.1. Then, an eigenvalue problem associated with the second variation of area is introduced in Sect. 2.2. In Sect. 3, the stability criteria for unduloids is presented in terms of the eigenvalues, mean curvature, and volume. The stability of unduloids is examined in Sect. 4, using the criteria prepared in the previous section. Section 5 is devoted to summary and discussions. The proofs of mathematical propositions and the method to compute geometric quantities of unduloids are presented in Appendices A and B, respectively. Figure 1: Left: An axially-symmetric hypersurface in R n+2 with the z-axis as the axis of rotation. Its profile curve is given by a function h = h(z) that represents the distance from the z-axis. Right: A part of an unduloid.
2 Variation and eigenvalue problem 2.1 Area, volume, and their variations We consider axially symmetric hypersurfaces in the closed domain of R n+2 (n ∈ N) bounded by two hyperplanes Π 1 := {z = z 1 } and Π 2 := {z = z 2 }. The local radius of a hypersurface is represented by height function h(z) (see Fig. 1 It is convenient to consider a one-parameter family of height function h(z, ǫ), where ǫ is a variation parameter. Then, the surface area and bulk volume of the axially symmetric object between z = z 1 and z = z 2 are A = a n z2 z1 Here, a n and v n+1 are the volume of a unit n-sphere and that of a unit (n + 1)-ball, respectively, given by a n = (n + 1)v n+1 = 2π A partial derivative is denoted by a subscript as h z := ∂ z h hereafter. The mean curvature of the hypersurface is For a cylinder, hemisphere, and unduloid, H takes a negative value in the present convention. The calculation of variations is equivalent to obtain the coefficients of the following expansion, The coefficient of expansion is obtained by X ℓ = ∂ ℓ ǫ X| ǫ=0 (ℓ = 0, 1, 2, · · · ). The first variations of area and volume are easily obtained in terms of h 0 (z) and h 1 (z), From Eqs. (6) and (7), one sees that the hypersurface which is the surface of revolution of h 0 is an equilibrium configuration or a critical point if and only if the following conditions hold, The CMC condition for the equilibrium configuration (8) corresponds to the Young-Laplace relation in fluid mechanics [20]. Now, let us focus on the volume-preserving variation (V ǫ ≡ 0) of CMC hypersurface, for which Eqs. (8) and (9) hold. For such a variation, the first derivative of area can be written as = −(n + 1)a n z2 z1 Then, the second derivative of area is Using Eqs. (8), (9) and (12), one obtains the second variation of area in terms of h 0 , h 1 , and H 1 , It is noted that A 2 is independent of h 2 due to the addition of term (n + 1)H 0 V ǫ in Eq. (10). The first variation of mean curvature H 1 in Eq. (13) can be written as by defining the following linear operator Therefore, A 2 is written in a simple form,

Eigenvalue problem associated with second variation of area
An equilibrium is defined to be stable if the second variation is non-negative for all volume-preserving variations. This condition is equivalent to A 2 ≥ 0 for all variations satisfying V 1 = 0. From this viewpoint, let us consider the following eigenvalue problem associated with A 2 .
The general variation of the height function is a linear combination of the eigenfunctions h 1 (z) = ∞ i=1 c i ϕ i (z), c i ∈ R. Then, A 2 and V 1 are written in terms of c i and λ i , where the orthonormality z2 z1 ϕ i ϕ j dz = δ ij is assumed. From Eqs. (19) and (20), one sees that an equilibrium h 0 is stable if λ 1 ≥ 0 since in such a case A 2 > 0 for all non-trivial (i.e., h 1 ≡ 0) volume-preserving variations satisfying V 1 = 0. One the other hand, one sees that an equilibrium h 0 is unstable if λ 2 is negative. Namely, (21) holds since in such a case A 2 < 0 for the volume-preserving variation given by For a uniform cylinder h 0 ≡ r = const., Eq. (18) is If one puts z 1 = 0, z 2 = L (> 0), the eigenvalue of a cylinder λ cyl i is obtained by solving Eq. (23), From Eq. (24), one can see that if λ cyl 2 < 0 holds and such a thin cylinder is unstable from criterion (21) (see also Refs. [21,22] for a dynamical counterpart). More precisely, it is proved that the cylinder with radius r and length L is stable if and only if r ≥ r c holds (cf. [11]). We call the cylinder with critical radius r c a critical cylinder.
The sphere S n+1 and the hemisphere with a boundary in either z = z 1 or z = z 2 are stable because S n+1 is the minimizer of area among all closed hypersurfaces enclosing the same volume.

Stability criteria of unduloids
It is convenient to introduce a quantity parameterizing the family of unduloids. As such a quantity, we adopt the non-uniformness parameter introduced in Sect. 1 where h min and h max denote the radii of an unduloid at the neck and bulge, respectively. One can naturally assign s = 0 and s = 1 to the critical cylinder and the largest hemisphere, that fits the interval, respectively. In the rest of this paper, we denote the half period of unduloid itself, mean curvature, volume, and eigenvalue of such an unduloid by U(s), H(s), V (s), and λ i (s), respectively. 3 For U(s), one can show the negativity (resp. positivity) of λ 1 (resp. λ 3 ). Namely, the following holds, In the rest of this section, we will introduce mathematical theories which play crucial roles in the stability analysis of Sect. 4. In Sect. 3.1, we see how to determine the sign of λ 2 (s) from the behavior of H(s). While λ 2 (s) < 0 immediately implies the instability of U(s) from (21), another criterion is needed to determine the stability of U(s) when λ 2 (s) ≥ 0. Therefore, in Sect. 3.2, we see how the behavior of H(s) and V (s) determines the stability when λ 2 (s) ≥ 0.

Sign of second eigenvalue λ 2
From Eq. (24), one can see that the second eigenvalue of cylinder λ cyl 2 increases and changes sign from negative to positive as radius r increases. From the point where λ 2 (0) = 0, two branches of unduloid 4 emanate (see Fig. 2). For these branches of unduloids bifurcating from the critical cylinder, the sign of λ 2 (s) is determined by the relative value of the mean curvature to that of the critical cylinder. Namely, if the mean curvature of the emanating unduloid H(s) is larger (resp. smaller) than that of the critical cylinder H(0), the second eigenvalue of the unduloid λ 2 (s) is negative (resp. positive), This statement is derived from Theorem A.2 in Appendix A.4. The criterion is visualized in Fig. 2. This criterion can be utilized to determine the sign of λ 2 (s) for U(s) close to the cylinder 0 < s ≪ 1.
In order to know when λ 2 (s) changes sign, the following criteria are quite useful.
Their proofs are presented in Appendix A.5. What criteria (29) and (30) mean is that, under the assumption that H ′′ (s) = 0 and V ′ (s) = 0, λ 2 (s) changes sign when H ′ (s) does. Although the possibility that λ 2 (s) vanishes even when H ′ (s) = 0 is not excluded by (29), it can be proved that See Prop. A.2 in Appendix A.6. Thus, once the sign of λ 2 (s) near s = 0 (the critical cylinder) is determined by (28), the sign of λ 2 in the full range of s is known by investigating H ′′ (s) and V ′ (s) at zeros of H ′ (s).

Criteria when λ 2 ≥ 0
While λ 2 (s) < 0 immediately implies that an unduloid U(s) is unstable from (21), another criteria is necessary to determine the stability of U(s) when λ 2 (s) ≥ 0. From the criteria for the stability given by Lemma A.2 in Appendix A.2, we have the following observations.
When λ 2 (s) = H ′ (s) = 0 and V ′ (s) = 0 hold, U(s) is unstable. Namely, the following holds, When λ 2 (s) > 0, the stability is related to the increase and the decrease of the mean curvature and volume.
We will utilize criteria (32) and (33) to determine the stability of U(s) whose λ 2 (s) is non-negative.

Comment: No iteration is needed
In the next section, we numerically obtain the mean curvature and volume for each U(s). Before starting such an analysis, let us see that obtaining H(s), V (s), and their derivatives numerically is a much simpler task than solving eigenvalue problem (18). H(s) and V (s) can be computed by just obtaining the 'background' solution h 0 (z). The function h 0 (z) is obtained by solving H 0 (z) = const. with boundary conditions h 0z (z 1 ) = h 0z (z 2 ) = 0. At a first glance, this problem seems to be a two-boundary problem requiring an iterative integration. By reducing H 0 (z) = const., which is a second-order ODE (ordinary differential equation), to an equivalent potential problem (a first-order ODE) and introducing an appropriate parameterization, however, no iteration turns out to be needed and the geometric quantities of unduloids, H(s) and V (s), are obtained by just estimating several improper integrals numerically (see Appendix B).
On the other hand, the eigenvalue equation (18) is essentially a two-boundary problem requiring an iteration procedure such as the shooting method [23]. Furthermore, one has to numerically solve the "perturbation equation" (18) for ψ 2 (z) and λ 2 on the numerical background h 0 (z), which is a part of operator L in Eq. (15).
Thus, it is stressed here that the stability criteria presented in Sects. 3.1 and 3.2 are not only easy to use but also enormously reduce the amount of numerical computations required in the analysis. This demonstrates the merit of adopting the geometric variational method throughout in our analysis, rather than ordinary mode-decomposition methods which are the standard for stability analysis in physics.

Stability of unduloids in
What are needs to do in order to examine the stability of all unduloids is to obtain the height function h 0 (z) corresponding to the half period of unduloid by numerically integrating the ODE H 0 (z) = const. with boundary conditions h 0z (z 1 ) = h 0z (z 2 ) = 0, while taking the dimension n ∈ N and non-uniformness s ∈ (0, 1) as free parameters. Then, one can estimate the mean curvature H and volume V as functions of s for each n. 5 Finally, utilizing the stability criteria (21) and (28)-(33), one can determine the stability of every unduloid.
In the rest of this section, we will clarify the stability of unduloids in all dimensions and parameter regimes of s. According to the behaviors of geometric quantities, we classify the dimensions into four classes, A (1 ≤ n ≤ 6), B (n = 7), C (8 ≤ n ≤ 9), and D (n ≥ 10), and examine the stability separately. Qualitative features of diagrams and stability structures are common in each class. The results in a final form are summarized as four statements (I)-(IV) in Sect. 5.
The characteristic area-volume diagrams of the unduloid, cylinder, and hemisphere are shown in Fig. 3. In addition, the numerical plots of H ′ (s) and V ′ (s), the derivatives of mean curvature and volume of unduloid U(s) with respect to s, are presented also in Fig. 3.
The area in the area-volume diagram is normalized in such a way that the area of the hemisphere remains unity in all ranges of the volume. The volume is normalized in such a way that the volume of the largest hemisphere, which fits the interval [z 1 , z 2 ], is unity. In the plots of H ′ (s) and V ′ (s), H ′ (s) and V ′ (s) are normalized by lim s→1−0 |H ′ (s)| = 1 and lim s→1−0 |V ′ (s)| = 1, respectively.  The sign of H ′ (s), λ 2 (s), and V ′ (s) and the stability of unduloid U (s) in R n+2 as functions of the non-uniformness parameter s ∈ (0, 1) in Class A, B, C, and D.
It is characteristic of this class that any unduloid has larger area than the cylinder and hemisphere with the same volume, and the area-volume curve of an unduloid has no cusp. From Fig. 3, one can see that H ′ (s) > 0, ∀s ∈ (0, 1). Therefore, the sign of λ 2 (s) is definite in the range of s by criteria (29) and (31). From the fact that H ′ (s) > 0 in a vicinity of s = 0, H(s) > H(0) holds for the unduloid which just bifurcated from the critical cylinder, which implies λ 2 (s) < 0 near s = 0 with criterion (28) (see also Fig. 2). Thus, λ 2 (s) < 0 holds ∀s ∈ (0, 1), from which (21) implies all unduloids are unstable in this class. This is consistent with known results in the literature, i.e., statement (ii).

Class B: n = 7
It is characteristic of this class that the area-volume curve of an unduloid has two cusps which form a swallowtail shape. We observe that both H ′ (s) and V ′ (s) have two simple zeros, which we denote by s k (k = 0, 1, 2, 3) as From the behavior of H ′ (s), one knows that λ 2 (s) vanishes and changes sign only at s = s 0 and s = s 3 with criteria (29)- (31). From this fact and the behavior of H ′ (s) with criterion (28), one can see that λ 2 (s) < 0 (resp. λ 2 (s) ≥ 0) for s ∈ (0, s 0 ) ∪ (s 3 , 1) (resp. s ∈ [s 0 , s 3 ]). Therefore, U(s) for s ∈ (0, s 0 ) ∪ (s 3 , 1) is unstable with criterion (21). Since λ 2 (s) ≥ 0 for s ∈ [s 0 , s 3 ], we have to see also the behavior of V ′ (s) in order to use criteria (32) and (33). From Fig. 3, V ′ (s) vanishes at neither s = s 0 nor s = s 3 , which with (32) The stability of the unduloid depending on s is summarized in Table 1, and values of s k (k = 0, 1, 2, 3) numerically obtained are presented in Table 2.
As mentioned before, the existence of stable unduloids for n = 7 has not been known. Thus, the stable unduloid for s ∈ [s 1 , s 2 ] is found for the first time in this paper.

Class
In this class the area-volume curve of an unduloid has two cusps as in Class B. The quantity V ′ (s) has two simple zeros, but H ′ (s) has only one simple zero. Taking into account the correspondence with Class B, we denote these zeros as follows.

Class D: n ≥ 10
In this class the area-volume curve of an unduloid has only one cusp. Both H ′ (s) and V ′ (s) have a simple zero. Taking into account the correspondence to the other classes, we denote the zeros as follows.

Summary and discussions
The equilibrium configuration of continuum supported by only surface tension (i.e., ignoring external gravity and self-gravity) is well approximated by a CMC (constant-mean-curvature) surface that is a critical point of the variational problem extremizing the surface area while keeping the volume fixed. We have investigated the stability of CMC hypersurfaces in R n+2 (n ∈ N) that possibly have boundaries on two parallel hyperplanes, by examining if the CMC hypersurfaces not only extremize but also minimize the surface area amongst all nearby surfaces while keeping the volume fixed. In particular, we have focused on the stability of non-uniform liquid bridges, known as the Delaunay unduloids [5], for which stability had been known partially [6,7,10,11] as statements (i)-(iv) presented in Sect. 4.
We have revealed the stability of unduloids for all n ∈ N and for all range of non-uniformness parameter s ∈ (0, 1), defined by Eq. (26). After obtaining mean curvature H and volume V of unduloids as functions of s numerically, the stability of unduloids was determined using their derivatives and stability criteria (21) and (28)-(33) presented in Sect. 3.
Although the behaviors of both H(s) and V (s) have played the central roles in our stability analysis, an interesting point is that the regions of s where the unduloid is stable (resp. unstable) completely coincide Table 2: Values of s k (k = 0, 1, 2, 3) for several n. There exists no s k (k = 0, 1, 2, 3) in Class A (1 ≤ n ≤ 6). The existence of s k (k = 0, 1, 2, 3) is the same in each class. n 7 8 9 10 11 s 0 0.437 n/a n/a n/a n/a s 1 0.   with those where V (s) is non-increasing (resp. increasing) for any n (see Table 1). Therefore, the bottom line of the stability analysis is summarized without mentioning H(s) as follows.
Let s ∈ (0, 1) be the non-uniformness parameter of a half period of an unduloid between two parallel hyperplanes in R n+2 (n ∈ N) defined by Eq. (26). The half period of an unduloid with parameter s and its bulk (n + 2)-volume are denoted by U(s) and V (s), respectively (the distance between the two hyperplanes is fixed). Then, the following (I)-(IV) hold.
The values of s 1 and s 2 are presented in Table 2 with other characteristic values, s 0 and s 3 (see Eqs. (34), (35), and (36) for the definitions).
We have not paid much attention to hemispheres and cylinders since their stability structure is completely understood as mentioned in Sect. 2.2. Nevertheless, let us have a look at them here, from which one can see the inevitability of the region where V ′ (s) ≤ 0 for n ≥ 8. As can be seen in Fig. 3, the area-volume curves of hemisphere and cylinder intersect for n ≤ 7, but not for n ≥ 8. A crucial reason of this is that the ratio of the volume of the critical cylinder V (0) to that of the largest hemisphere (i.e., the hemisphere which fits the interval [z 1 , z 2 ]) V (1), given by increases with n and becomes larger than unity for n ≥ 8 (see Table 3). Namely, for n ≥ 8 the branch of unduloids emanating from the critical cylinder at s = 0 must have a region where the volume decreases to reach the largest hemisphere at s = 1. Although dimension n = 7, at and above which the stable unduloid exists, differs from this critical number of dimension n = 8 by one, their closeness is clearly not a coincidence. In passing, let us point out that the area-volume curve deforms continuously if n changes continuously. As n increases from 1, the 'swallowtail' (two cusps) of the area-volume curve appears at n ≃ 7. As n increases further, the swallowtail becomes large. In other words, s 1 decreases to approach 0 and s 2 increases to approach 1. Indeed, s 1 decreases as n increases to vanish finally at n ≃ 10. As far as we know, s 2 continues to increase but does not vanish for arbitrarily large n ≫ 1, which is consistent with statement (i). If one treats n as a continuous parameter and examines the stability for non-integer n, which seems to bring no technical problem, the behaviors of area-volume curve and stability structure expected and described above would be observed.
In this paper, the stability of unduloids was determined by the behaviors of H(s) and V (s), which were obtained by numerical integration. Therefore, the correctness of the conclusions is based on that of these numerical computations. It is noted that one needs highly accurate computation to show that s 2 < s 3 holds (s 2 and s 3 are defined as the zeros of V ′ (s) and H ′ (s), respectively) for n ≥ 8. For example, s 3 /s 2 −1 ≃ +5.6×10 −6 for n = 9 by our computation, and this quantity seems to decrease further as n increases. Nevertheless, we assumed that s 2 < s 3 continues to hold for arbitrarily large n, otherwise our conclusions on the stability might be different from those presented in the text. Therefore, any analytic method or alternative numerical methods that guarantee accuracy will be helpful to confirm the results in this paper.
Related to the results of this paper, one of the most interesting problems would be to investigate the implications to dynamical problems. While this was partially worked by one of the present authors in [22] using the surface-diffusion equation [24,25], there are still many things to do in this direction.
We remark that the stability of black strings qualitatively exhibits a similar dependence on the dimension. Suppose a D-dimensional vacuum spacetime (D ≥ 5) with one spacelike dimension compactified to a circle S 1 . Then, there exist non-uniform black strings of which horizon topology is S D−3 × S 1 . The stability of such black strings has been examined using the thermodynamic criterion, and argued as follows [18]. If 5 ≤ D ≤ 11, all non-uniform black strings are unstable. If 12 ≤ D ≤ 13, there exists a critical non-uniformness below (resp. above) which the non-uniform black strings are unstable (resp. stable). If D ≥ 14, all black strings are stable. We are not so surprised at the similarity of stability between these black objects and CMC hypersurfaces since it was shown that the event horizon of a black hole is approximated by a CMC hypersurface in the large-dimension limit [19]. Nevertheless, it is still interesting to pursue the similarity from various points of view such as the fluid/gravity correspondence [26,27] and the gauge/gravity correspondence [28].

A Mathematical propositions and their proofs A.1 Stability of a half period and one period of an unduloid
Here, we prove Proposition A.1 If a half period of an unduloid is stable (resp. unstable) between two parallel hyperplanes in R n+2 (n ∈ N), one period of the unduloid is stable (resp. unstable) in R n+1 × S 1 , and vice versa.
Proof. Let X be a half period of an unduloid U with the z-axis as its axis of revolution which is generated by a curve This implies that X is perpendicular to the hyperplanes Π 0 = {z = 0}, Π 2 = {z = z 2 }. Without loss of generality, we may assume that z = 0, z = z 2 corresponds to the bulge and the neck of U, respectively. Denote by Y the one period of U that is generated by Assume that X is unstable. Then, there exists a volume-preserving variation X(ǫ) of X such that holds. By reflection with respect to Π 0 , we get a volume-preserving variation Y (ǫ) of Y which satisfies This implies that Y is also unstable. Assume now that Y is unstable. Then, there exists a volume-preserving variation Y (ǫ) of Y such that holds. LetŶ (ǫ) be the Steiner symmetrization of Y (ǫ) with respect to Π 0 , that is,Ŷ (ǫ) is a hypersurface defined by the conditions (i) and (ii) below. Note that we consider only hypersurfaces close to Y . Set respectively. For each point P ∈ Π 0 , denote by L P the straight line that passes through P and is perpendicular to Π 0 . Define two straight line segments byΓ(ǫ) := L P ∩Ĝ(ǫ), Γ(ǫ) := L P ∩ G(ǫ). Note the following.

A.2 Stability criteria for axially symmetric equilibrium hypersurfaces
define an axially-symmetric equilibrium hypersurface X with the z-axis as its axis of revolution, that is, X is a part of either a cylinder or an unduloid with generating curve (40) and it is perpendicular to the hyperplanes Π i = {z = z i }, (i = 1, 2). Then, one can show the following lemma.
Lemma A.1 X is stable if and only if X is stable for axially-symmetric variations.
(ii) The center ofD(ǫ) lies on the z-axis. Then, it is well-known that (a) V (G(ǫ)) = V (Ĝ(ǫ)) holds, (∀ǫ), (b) A(X(ǫ)) ≥ A(X(ǫ)) holds, (∀ǫ) hold (cf. [29,Note A]). Therefore,X(ǫ) is a volume-preserving axially-symmetric variation of X such that holds. Hence X is unstable for axially-symmetric variations. The opposite direction is trivial. Now we give two criteria for the stability of axially-symmetric equilibrium hypersurfaces. The first criterion (Lemma A.2) will be proved by using the second criterion (Lemma A.3) at the end of this subsection.
Lemma A.2 (First stability criterion) Assume that X(s) is a one-parameter smooth family of axiallysymmetric equilibrium hypersurfaces generated by the curves that is, X(s) are half periods of unduloids, and we assume that s is the parameter defined by (26). Denote by H(s), V (s) the mean curvature, the enclosed (n + 2)-dimensional volume of X(s), respectively. Denote by λ i (s) the i-th eigenvalue of the eigenvalue problem (18) for X(s).
The following observation will be used to prove Lemma A.2.
Lemma A.4 If X(s) defined by (41) are half periods of unduloids and s is the parameter defined by (26), then for any fixed s, h s (z, s) ≡ 0 (43) holds for all z ∈ [z 1 , z 2 ], where h s = ∂h/∂s.

Proof. Set
It is sufficient to prove that for any fixed ρ, h ρ (z, ρ) ≡ 0 holds for all z ∈ [z 1 , z 2 ], Now we may assume that holds. Differentiating both sides of (46) with respect to ρ, we obtain This implies that uh n dz = V ′ (s) (n + 1)a n H ′ (s) , which shows that (II) in Lemma A.2 is equivalent to (II) in Lemma A.3. Next, we prove (III). Let e be an eigenfunction belonging to λ 2 = 0. Then, First, assume that H ′ (s) = 0 holds. Then, since h s ≡ 0 (Lemma A.4), Eq. (49) implies that h s is an eigenfunction belonging to λ 2 = 0. Since each eigenspace is one-dimensional, Hence, Therefore, V ′ (s) = 0 is equivalent to z2 z1 eh n dz = 0. From (III-A) in Lemma A.3, X(s) is unstable. This gives (III-a).
Take the function u that is uniquely defined in (III-B) of Lemma A.3. Using (49), h s is written as Then, V ′ (s) = a n z2 z1 h s h n dz = (n + 1)a n H ′ (s) z2 z1 uh n dz (60) holds, which gives z2 z1 uh n dz = V ′ (s) (n + 1)a n H ′ (s) .
A.3 Negativity (resp. positivity) of the first (resp. third) eigenvalue for unduloids In addition to the eigenvalue problem (18), we consider also the following eigenvalue problem with Dirichlet boundary condition: and denote its i-th eigenvalue by λ 0 i [z 1 , z 2 ] = λ 0 i . Then, it is well-known that λ 0 1 < λ 0 2 < · · · and hold. Also recall that, since L is a Sturm-Liouville operator, each of ϕ i (z) and ψ i (z) has exactly i − 1 zeros in (z 1 , z 2 ). Recall that Eq. (14) holds and consider the parallel translation Then, the mean curvature H(ǫ) is the same as the mean curvature H 0 of U. U(ǫ) can be represented as Hence, we have from (14) that Since U is a half period of an unduloid, h z (z 1 ) = h z (z 2 ) = 0 holds and we may assume that h z > 0 on z 1 < z < z 2 . Hence, h z is an eigenfunction of (62) and the corresponding eigenvalue zero is the first eigenvalue λ 0 1 . Hence, by (63), holds. Next, we assume that holds. Then, eigenfunction ϕ 3 (z) has exactly two zeros ζ 1 , ζ 2 , (ζ 1 < ζ 2 ), in (z 1 , z 2 ). Hence, by the monotonicity of the eigenvalues of the problem (62) with respect to the domain, we have which is a contradiction. Hence, λ 3 must be positive.

A.4 Existence of bifurcation and estimate of eigenvalues in a bifurcation branch
Assume that I ⊂ R is an non-empty open interval, and defines a smooth one-parameter family of axially-symmetric equilibrium hypersurfaces {X(ǫ)} ǫ∈I , that is, each X(ǫ) is a part of either a cylinder or an unduloid and it is perpendicular to the hyperplanes Π i = {z = z i }, (i = 1, 2), and X(ǫ) is of C ∞ in ǫ. Denote by H(ǫ) the mean curvature of X(ǫ). Denote by λ i (X(ǫ)) the i-th eigenvalue of the problem (18) for X(ǫ). Now we define the concept "bifurcation instant".
Definition A.1 For ǫ ∈ I, we say that ǫ is a bifurcation instant for the family {X(ǫ)} ǫ∈I if there exists a sequence {ǫ k } k∈N in I and a sequence {Y k } k∈N such that: (ii) Each Y k is an axially-symmetric equilibrium hypersurface that is defined by and the mean curvature of Y k is equal to H(ǫ k ) for all k.
In other words, ǫ is a bifurcation instant for the family {X(ǫ)} ǫ∈I if X(ǫ) is an accumulation of equilibrium hypersurfaces that are not congruent to any of the hypersurface of the family {X(ǫ)} ǫ∈I .
Remark A.2 The following Theorems A.1, A.2 are proved by modifications of the proofs of Theorems 1.1, 6.4 in [16], respectively.
Theorem A.2 Under the assumptions of Theorem A.1, denote by {Y (s)} s∈Î the bifurcating branch of axiallysymmetric equilibrium hypersurfaces given in Theorem A.1. LetĤ(s) be the mean curvature of Y (s), and µ(s) the eigenvalue for the Jacobi operator L Y (s) which was defined by (15). We may assume that H ′ (0) > 0 holds, by changing the parameter t to −t if necessary. Then, the following statements are true.
A.5 Correspondence of sign change between the eigenvalue and mean-curvature derivative In this section we prove (29) and (30). Assume that {X(ǫ)} ǫ∈I satisfies the same assumptions as that in Section A.4. Denote by H(ǫ) the mean curvature of X(ǫ), and by V (ǫ) the enclosed (n + 2)-dimensional volume by X(ǫ). Denote by λ i (ǫ) the i-th eigenvalue of the problem (18) for X(ǫ).
The criterion (29) is proved as follows. Proof of the criterion (29). Assume that X(s 0 ) is a part of an unduloid. Assume also that H ′ = 0 holds at s 0 . If λ 2 changes sign at s 0 , we can see that s 0 is a bifurcation instant in the same way as in the proof of Proposition 2.14 in [31] which is an application of [32, Theorem 2.1]. However, in our variational problem, there is no bifurcation from any unduloid X(ǫ), which is a contradiction. Therefore, λ 2 does not change sign at s 0 , which proves the criterion (29).
The criterion (30) is given by the following Lemma A.5.
Then, there exists a non-zero real number α such that holds. In particular, λ j changes sign at ǫ.
Proof. The formula (2.8) in [33] is about a functional F : H × R → R, where H is a real Hilbert space. In our case, let H be the space of real-valued C ∞ functions on the interval [z 1 , z 2 ] ⊂ R with inner product For any axially-symmetric (not necessarily equilibrium) hypersurface X h generated by and we define another parameter H as −a := (n + 1)H.
A.6 Equivalence between the vanishing of the eigenvalue and mean-curvature derivative Proof. We prove the following (i), (ii) one by one.
Lemma A. 6 Assume that X(s) satisfies the same assumptions as in Proposition A.2. Then, h z (0, s) = 0 and h zz (0, s) = 0 hold. and h max is a strictly-decreasing function in the family of unduloids with the cylinder as the initial surface. Hence, However, if h zz (0, s) = 0 holds, using h z (0, s) = 0, from (4), we have with the following Lagrangian, Here, H is a constant representing the mean curvature of the hypersurface. Since Lagrangian (94) does not depend on z explicitly, the following quantity is conserved, which is called the Beltrami identity. Substituting (94) into (95), we obtain an equation like the law of conservation of mechanical energy: Introducing a new variable w by we have Denote the zeros of U (w) by w ± (0 < w − < w + ). Then, it is easy to see that where ρ = 1 − s is given by Eqs. (26) and (44). Using U (w ± ) = 0, one can express w ± and K as functions of ρ, Here, the last expressions in Eqs. (101)-(103) are convenient to avoid the round-off errors in the numerical estimation for 0 < ρ ≪ 1 and large n.
From Eq. (98) with the assumption that dw dz ≥ 0, we obtain Integrating the left-hand (resp. right-hand) side this with respect to z (resp. w) from z 1 to z 2 (resp. from w − to w + ), we obtain where L := z + −z − is the half period of the unduloid. This is an integral representation of the mean curvature of an unduloid, which is a function of ρ (or equivalently s) and L. Then, using Eqs. (1), (2), (97), and (104), the integral representations of the area and volume of a half period of an unduloid are obtained as A = a n (−H) n+1 which are also functions of ρ and L.
Without loss of generality, we can fix the interval L (L = 1 for example) throughout the stability analysis. Thus, we have expressed the mean curvature, volume, and area of an unduloid as integrals essentially depending only on the non-uniformness parameter ρ (or s), Eqs. (105), (106), and (107). Therefore, as argued in Sect. 3.3, we do not have to solve H 0 (z) = const. with boundary conditions h 0z (z 1 ) = h 0z (z 2 ) = 0, which needs an iterative integration to satisfy the boundary conditions at the both boundaries.

B.2 Manipulation for accurate numerical integration
What is necessary to obtain A(s), V (s), and H(s), which play essential roles in the stability analysis, is to estimate the integrals numerically in Eqs. (105), (106), and (107) as accurately as possible. Since U (w) vanishes at the both ends of integral range, the following manipulation helps us to estimate the integrals numerically [13]. Those who are not interested in the numerics do not need to read the rest of this section.
Here, the right-hand side defines the polynomial expression of g(w), which is regular at w = w ± . The comparison of coefficients yields the following recursion relation and "boundary conditions" to be satisfied by g p (0 ≤ p ≤ n − 1), (w + + w − )g n−1 − g n−2 = 1, (113) These can be easily solved to give the following expression of general term, Again, the final expression is for the avoidance of round-off error. Finally, we fix the integration range as by changing variable from w to ζ, One can accurately estimate H, V , and A numerically for given s = 1 − ρ ∈ (0, 1) (after fixing L, L = 1 for example) using Eq.

B.3 Area-volume diagrams
We describe here how to draw the area-volume diagrams in Fig. 3.
First, let us normalize the volume V by the volume of the largest hemisphere, which has a radius identical to the length of the interval L := z 2 − z 1 , and normalize the surface area A by the surface area of the hemisphere whose radius is R ∈ (0, L],V Here, V and A are the volume and surface area, respectively, of a hemisphere, cylinder, or half-period of unduloid.
For the hemisphere with radius R ∈ (0, L], the normalized volume and area arê respectively. These give a parametric representation of the area-volume curve of hemisphere in Fig. 3 with R/L ∈ (0, 1] being the parameter. For the cylinder with radius r ∈ (0, +∞), the normalized volume and area are given bŷ A cyl = a n r n L . (127) Equations (125) and (127) give a parametric representation of the area-volume curve of cylinder in Fig. 3 with r/L ∈ (0, +∞) being the parameter.
Denoting the volume and area of a half period of unduloid in the form of Eq. (119) by V (s) and A(s), the normalized volume and area of unduloid are given bŷ SolvingV hem =V und for R to obtain R = L(V und ) 1 n+2 , and substituting this into the right-hand side of Eq. (129), one obtainsÂ Equations (128) and (130) give a parametric representation of the area-volume curve of unduloids in Fig. 3 with s ∈ (0, 1) being the parameter.