Symmetric self-shrinkers for the fractional mean curvature flow

TL;DR

This paper demonstrates the existence of symmetric, shrinking solutions to the fractional mean curvature flow with multiple concentric spheres, highlighting their stability properties.

Contribution

It introduces new symmetric solutions for the fractional mean curvature flow and analyzes their stability, expanding understanding of the flow's dynamics.

Findings

Existence of homothetically shrinking solutions with concentric spheres.

All solutions except the ball are dynamically unstable.

Provides insights into the stability of symmetric solutions.

Abstract

We show existence of homothetically shrinking solutions of the fractional mean curvature flow, whose boundary consists in a prescribed numbers of concentric spheres. We prove that all these solutions, except from the ball, are dynamically unstable.

Equations

$$\partial_{t}x\cdot\nu=-H_{s}(x,E_{t}):=-\lim_{\varepsilon\to 0}\int_{\mathbb{R}^{n}\setminus B_{\varepsilon}(x)}\Big{(}\chi_{\mathbb{R}^{n}\setminus E_{t}}(y)-\chi_{E_{t}}(y)\Big{)}\ \frac{1}{|x-y|^{n+s}}\,dy,$$

$$x\cdot \nu =c{H}_s(x,E)\text{for some constant c>0. }$$

$$H_s(x_r,B_r)=\frac{k(n)}{r^s}\text{where }k(n):=H_s(x_1,B_1).$$

$${\int }_{\partial B_{|x|+\delta }}\frac{1}{|y-x{|}^{n+s}}dy=\frac{1}{|\delta {|}^{1+s}}(c+o(1)),$$

$${\int }_{\partial B_{|x|+\delta }}\frac{1}{|y-x{|}^{n+s}}dy={\int }_{\partial B_{|x|+\delta }(|x|e_1)}\frac{1}{|y^{\prime }{|}^{n+s}}d{y}^{\prime }$$

$${\int }_{\partial B_{|x|+\delta }(|x|e_1)\cap \{y^{\prime }\cdot e_1>|x|\}}\frac{1}{|y^{\prime }{|}^{n+s}}d{y}^{\prime }⩽\frac{n{\omega }_n(|x|+\delta {)}^{n-1}}{2|x{|}^{n+s}}⩽C.$$

$${\int }_{\partial B_{|x|+\delta }(-x)\cap \{y^{\prime }\cdot e_1<|x|\}}\frac{1}{|y^{\prime }{|}^{n+s}}d{y}^{\prime }={\int }_{B_{|x|+\delta }^{\prime }}\frac{|x|+\delta }{\sqrt{(|x|+\delta {)}^2-|z{|}^2}}\frac{1}{{\left(f(z{)}^2+|z{|}^2\right)}^{\frac{n+s}{2}}}dz$$

$$g_{\delta }(\rho ):=\frac{R_{\delta }}{\sqrt{R_{\delta }^2-{\rho }^2}}{\rho }^{n-2}{\left[{\left(\frac{\delta }{|\delta |}-\frac{{\rho }^2}{R_{\delta }+\sqrt{R_{\delta }^2-{\rho }^2}}\right)}^2+{\rho }^2\right]}^{-\frac{n+s}{2}}.$$

$${\int }_{R_{\delta }/2}^{R_{\delta }}{g}_{\delta }(\rho )d\rho ⩽\frac{C}{R_{\delta }^{s+1}}=C\frac{|\delta {|}^{s+1}}{(|x|+\delta {)}^{s+1}}.$$

$$g_{\delta }(\rho ){\chi }_{(0,R_{\delta }/2)}(\rho )⩽C{\chi }_{(0,1)}(\rho )+C\frac{1}{{\rho }^{s+2}}{\chi }_{(1,+\infty )}(\rho )\in L^1(0,+\infty ),$$

$$g_{\delta }(\rho )\to \frac{{\rho }^{n-2}}{(1+{\rho }^2{)}^{\frac{n+s}{2}}}\text{as δ→0},$$

$${\int }_{\partial B_{|x|+\delta }(-x)\cap \{y_1^{\prime }<|x|e_1\}}\frac{1}{|y^{\prime }{|}^{n+s}}d{y}^{\prime }=\frac{(n-1){\omega }_{n-1}}{|\delta {|}^{s+1}}\left[{\int }_0^{\frac{R_{\delta }}{2}}{g}_{\delta }(\rho )d\rho +{\int }_{\frac{R_{\delta }}{2}}^{R_{\delta }}{g}_{\delta }(\rho )d\rho \right]=\frac{c+o(1)}{|\delta {|}^{s+1}}$$

$$c:=(n-1){\omega }_{n-1}{\int }_0^{+\infty }\frac{{\rho }^{n-2}}{(1+{\rho }^2{)}^{\frac{n+s}{2}}}d\rho .$$

$$A:=B_R\setminus B_r$$

$$1=c{H}_s(x_1,A)\text{ for all x1 with ∣x1∣=1}\text{ and }r=-c{H}_s(x_r,A)\text{ for all xr with ∣xr∣=r}$$

$$H_s(x_r,A)=-r{H}_s(x_1,A).$$

$$f_s(r)=H_s(x_r,A)+r{H}_s(x_1,A).$$

$$H_s(x_1,A_{1,r^{\prime }})>H_s(x_1,A)⩾H_s(x_1,B_1)=k(n)>0.$$

$$r\in (0,1)\mapsto H_s(x_1,A)\text{ is monotone increasing and positive}.$$

$$H_s(x_r,A)=-H_s(x_r,B_r)+2{\int }_{{\mathbb{R}}^n\setminus B_1}\frac{1}{|x_r-y{|}^{n+s}}dy.$$

$${\int }_{{\mathbb{R}}^n\setminus B_1}\frac{1}{|x_r-y{|}^{n+s}}dy<{\int }_{{\mathbb{R}}^n\setminus B_1}\frac{1}{|x_{r^{\prime }}-y{|}^{n+s}}dy$$

$$r\in (0,1)\mapsto H_s(x_r,A)\text{ is monotone increasing}.$$

$$H_s(x_r,A)={\int }_{{\mathbb{R}}^n\setminus B_1}\frac{1}{|x_r-y{|}^{n+s}}dy+\underset{\epsilon \to 0}{\lim }\left({\int }_{B_{r-\epsilon }}\frac{1}{|x_r-y{|}^{n+s}}dy-{\int }_{B_1\setminus B_{r+\epsilon }}\frac{1}{|x_r-y{|}^{n+s}}dy\right).$$

$$H_{s}(x_{r},A)=H_{s}(x_{r},B_{r})+\lim_{\varepsilon\to 0}\left(2\int_{B_{r-\varepsilon}}\frac{1}{|x_{r}-y|^{n+s}}dy-2\int_{B_{1}\setminus B_{r+\varepsilon}}\frac{1}{|x_{r}-y|^{n+s}}dy\right)\\ =\frac{k(n)}{r^{s}}+\lim_{\varepsilon\to 0}\left(2\int_{\varepsilon}^{r}\int_{\partial B_{r-\delta}}\frac{1}{|x_{r}-y|^{n+s}}dyd\delta-2\int_{\varepsilon}^{1-r}\int_{\partial B_{r+\delta}}\frac{1}{|x_{r}-y|^{n+s}}dyd\delta\right)\\ =\frac{k(n)}{r^{s}}+2\int_{1-r}^{r}\left[\int_{\partial B_{r-\delta}}\frac{1}{|x_{r}-y|^{n+s}}dy-\int_{\partial B_{r+\delta}}\frac{1}{|x_{r}-y|^{n+s}}dy\right]d\delta\\ =\frac{k(n)}{r^{s}}+2\big{(}c+o(1)\big{)}\left(\frac{1}{(1-r)^{s}}-\frac{1}{r^{s}}\right).$$

$$E:=\bigcup _{k=1}^N\left(B_{r_{2k}}\setminus B_{r_{2k-1}}\right)$$

$$\tilde{E}:=B_{{\tilde{r}}_0}\cup \bigcup _{k=1}^N\left(B_{{\tilde{r}}_{2k}}\setminus B_{{\tilde{r}}_{2k-1}}\right)$$

$$f_i(r_1,\dots ,r_{2N-1}):=r_i{H}_s(x_1,E)+(-1{)}^{i-1}{H}_s(x_{r_i},E)=0\forall i\in \{1,\dots ,2N-1\}$$

$$f_i(r_0,\dots ,r_{2N-1}):=r_i{H}_s(x_1,\tilde{E})+(-1{)}^{i-1}{H}_s(x_{r_i},\tilde{E})=0\forall i\in \{0,\dots ,2N-1\}.$$

$$\underset{r_0\to 0}{\lim }{H}_s(x_1,\tilde{E})=H_s(x_1,A_{1,r_1})>0\underset{r_0\to r_1}{\lim }{H}_s(x_1,\tilde{E})=H_s(x_1,B_1)=k(n)>0.$$

$$H_s(x_{r_0},\tilde{E})=-2{\int }_{B_1\setminus B_{r_1}}\frac{1}{|x_{r_0}-y{|}^{n+s}}dy+H_s(x_{r_0},B_{r_0})=-2{\int }_{B_1\setminus B_{r_1}}\frac{1}{|x_{r_0}-y{|}^{n+s}}dy+\frac{k(n)}{r_0^s},$$

Full text

The authors were supported by the Italian INDAM–GNAMPA and by the University of Pisa via grant PRA 2017 Problemi di ottimizzazione e di evoluzione in ambito variazionale.

Annalisa Cesaroni

Department of Statistical Sciences

University of Padova

Via Cesare Battisti 141, 35121 Padova, Italy

Matteo Novaga

Department of Mathematics

University of Pisa

Largo Bruno Pontecorvo 5, 56127 Pisa, Italy

1. Introduction

Let us introduce the geometric evolution which we consider in this paper. Given an initial set $E\subset\mathbb{R}^{n}$, we define its evolution $E_{t}$ according to fractional mean curvature flow as follows: the velocity at a point $x\in\partial E_{t}$ is given by

[TABLE]

where $s\in(0,1)$ is a fixed parameter and $\nu$ is the outer normal at $\partial E_{t}$ in $x$. The fractional mean curvature of a set has been introduced in [MR2675483] as the first variation of the fractional perimeter functional, and it has been proved in [av] that for sufficiently smooth sets $E$ the rescaled fractional mean curvature $(1-s)H_{s}(x,E)$ converges as $s\to 1$ to the classical mean curvature of $E$ at $x$. The evolution law (1.1) can be interpreted as the $L^{2}$-gradient flow of the fractional perimeter.

Existence and uniqueness of viscosity solutions to a level set formulation of (1.1) has been provided in [i, cmp], and qualitative properties of smooth solutions have been studied in [SAEZ]. However, we point out that the short-time existence of smooth solutions has not yet been proved. In [cs] the convergence to the fractional mean curvature flow of a threshold dynamics scheme is proved; this result was adapted to the anisotropic case, even in presence of a driving force in [cnr], where it is also shown that the flow preserves convexity. It has also been observed that the geometric law (1.1) presents some different behavior with respect to the classical mean curvature flow: we refer for instance to the paper [csv1] about the formation of neck-pinch singularities, and to the paper [cdnv] about fattening and non-fattening phenomena.

In this paper we are interested in the homothetically shrinking solutions for the flow (1.1). A homothetic solution to (1.1) is a self-similar solution to (1.1): substituting $E_{t}=\lambda(t)E$ in (1.1), it is easy to see, using scale invariance of the fractional mean curvature, that this is equivalent to $\lambda^{\prime}(t)x\cdot\nu=-\frac{1}{\lambda(t)^{s}}H_{s}(x,E)$ for all $x\in\partial E$. So homothetically shrinking solutions to (1.1) are given by the solutions to (1.1) with initial datum every set $E\subseteq\mathbb{R}^{n}$ of class $C^{1,1}$ which satisfies

[TABLE]

Homothetically shrinking solutions are particularly relevant in the analysis of the classical mean curvature flow, as they are canonical examples of singularities, in the sense that any solution converges to a self-shrinker, if properly rescaled around a singular point. This result follows from an important monotonicity formula established by G. Huisken in [MR1030675] for the mean curvature flow. The analog of such formula in the fractional setting is still an open problem. We recall moreover that, at the moment, the existence theorem for local in time regular solutions of (1.1), even if expected, has not been proved.

It is well-known that the only embedded planar curve which is homothetically shrinking under curvature flow is the circle [MR845704], whereas in higher dimensions there exist other smooth embedded surfaces which are self-shrinkers for the mean curvature flow, starting from the rotationally symmetric torus discovered by Angenent [MR1167827], and then going to more complex configurations as punctured compact surfaces or non-compact asymptotically conical surfaces, see [MR1361726, kl]. However, it is easy to show that the ball is the only self-shrinker which is also radially symmetric.

In the fractional setting the classification of self-shrinkers is still at a very early stage. As far as we know, we provide here the first examples of fractional self-shrinkers which are different from balls and cylinders. More precisely, in Section 2 we show the existence of homothetic solutions to the flow (1.1) which are radially symmetric, and have a prescribed number of boundary spheres (see Theorem 2.3). Moreover, in the case of a single annulus, we show uniqueness of the ratio ${R/r}$ for which the flow starting from the annulus $B_{R}\setminus B_{r}$ self-similarly shrinks to a point. The existence of such radially symmetric self-shrinkers, different from balls, is a new feature compared with the local case, and it is due to the nonlocal nature of the fractional mean curvature.

A natural question arising about self-similar shrinkers is the issue of their dynamic stability. In the case of the classical mean curvature flow, the study of the dynamic stability of self-shrinkers was initiated in [cm], and later developed by other authors. From the convergence results in [MR840401, MR772132] it follows that the balls is dynamically stable under mean curvature flow (see also [MR1637988, np, pp] for a discussion of the stability of the Wulff-Shape as homothetic solution of the anisotropic and crystalline curvature flow). Moreover, in [cm] it is shown that balls and cylinders are the only stable self-shrinkers.

In the fractional case none of such results is currently available, in particular it is not known whether the ball is dynamically stable, and if convex sets shrink to a round point at the singular time. We discuss in this paper the stability issue for the class of solutions that we construct in Theorem 2.3. In particular, in Section 3 we show that the radial self-shrinkers different from the ball are all dynamically unstable (see Theorem 3.1).

Acknowledgements. The authors are members of INDAM-GNAMPA. The second author was partially supported by the University of Pisa Project PRA 2017 "Problemi di ottimizzazione e di evoluzione in ambito variazionale".

2. Existence of symmetric self-shrinkers

We start with a technical result which will be useful in the sequel. We denote by $B_{r}$ the ball of center [math] and radius $r>0$, and we let $B_{r}(x)=x+B_{r}$. Moreover, we recall that, by the scale invariance of fractional mean curvature, for all $x_{r}\in\partial B_{r}$ there holds (see [SAEZ, Lemma 2])

[TABLE]

Lemma 2.1**.**

Let $x\in\mathbb{R}^{n}\setminus\{0\}$ and $\delta\neq 0$. Then, as $\delta\to 0$, the following estimate holds:

[TABLE]

for a constant $c>0$ depending only on $n$ and $s$.

Proof.

Up to a rotation of the reference system, we can assume that $x=-|x|e_{1}$. By the change of coordinates $y^{\prime}=(y-x)$, we get,

[TABLE]

Note that

[TABLE]

Moreover we write $\partial B_{|x|+\delta}(|x|e_{1})\cap\{y^{\prime}\cdot e_{1}<|x|\}=\{(f(z),z)\,|\,z\in\mathbb{R}^{n-1},z\in B^{\prime}_{|x|+\delta}\}$, where $B^{\prime}_{r}\subseteq\mathbb{R}^{n-1}$ denotes the ball of center [math] and radius $r$ in $\mathbb{R}^{n-1}$ and $f(z)=|x|-\sqrt{(|x|+\delta)^{2}-|z|^{2}}$. Therefore, denoting $R_{\delta}:=\frac{|x|+\delta}{|\delta|}$, we get

[TABLE]

Let

[TABLE]

Now we observe that there exists $C=C(n,s)>0$, such that

[TABLE]

Moreover, taking $|\delta|$ sufficiently small such that $R_{\delta}>4$, we get that there exists a dimensional constant $C=C(n,s)>0$ such that

[TABLE]

where $\chi_{(a,b)}$ is the characteristic function of the interval $(a,b)$. Using (2.3) and observing that

[TABLE]

we conclude by (2.2) and by Lebesgue dominated convergence theorem that

[TABLE]

where

[TABLE]

The conclusion follows by this estimate and (2.1). ∎

First of all we look to the simplest example of rotationally symmetric set different from a ball. We show that there exists a unique value of the ratio $\frac{R}{r}$ which depends on the dimension $n$ and on the fractional power $s\in(0,1)$ such that the annulus $B_{R}\setminus B_{r}$ is a self-shrinker.

Proposition 2.2**.**

Let $n\geqslant 1$. Then, for all $R>0$ fixed there exists a unique $r=r(n,s)\in(0,R)$ depending only on $R$, $s\in(0,1)$ and $n$, such that the flow (1.1) with initial datum the annulus

[TABLE]

is a homothetically shrinking solution of the flow.

Proof.

Up to rescaling the set we fix $R=1$. We observe that $A$ is a solution to (1.2) if and only if for some $c>0$,

[TABLE]

and so if and only if

[TABLE]

By rotational invariance, we get that $H_{s}(x_{r},A),H_{s}(x_{1},A)$ do not depend on the points $x_{r},x_{1}$, but only on $0<r<1$. Moreover they are both continuous functions with respect to $r$, due to the continuity of the fractional mean curvature with respect to $C^{2}$-convergence of sets (see [cmp, Section 5.2]). We consider the following function defined for $r\in(0,1)$

[TABLE]

Note that the function $f_{s}$ is continuous on $(0,1)$. To prove the statement it is sufficient to show that there exists a unique $r=r(n,s)$ such that $f_{s}(r(n,s))=0$.

Let $r,r^{\prime}$ such that $0<r<r^{\prime}<1$. By the inclusions $A_{1,r^{\prime}}:=B_{1}\setminus B_{r^{\prime}}\subseteq A\subseteq B_{1}$, we get, by the monotonicity of the fractional mean curvature (see [cmp, Section 5.2]), that

[TABLE]

This implies that

[TABLE]

Moreover, we observe, recalling the definitions, that

[TABLE]

Note that if $r^{\prime}>r$ then $H_{s}(x_{r^{\prime}},B_{r^{\prime}})=\frac{k(n)}{(r^{\prime})^{s}}<\frac{k(n)}{r^{s}}=H_{s}(x_{r},B_{r})$, whereas for $1>r^{\prime}>r$, $|x_{r}-y|\geqslant|\frac{r^{\prime}}{r}x_{r}-y|$ for all $y\in\mathbb{R}^{n}\setminus B_{1}$, and $x_{r}$ with $|x_{r}|=r$. Therefore by symmetry of the kernel we have that

[TABLE]

for all $x_{r},x_{r^{\prime}}$ with $|x_{r}|=r,|x_{r^{\prime}}|=r^{\prime}$. Using these facts we conclude that

[TABLE]

Due to (2.7), (2.8), we notice that the function $f_{s}(r)$ defined in (2.5) is monotone increasing. Now we claim that $\lim_{r\to 0}f_{s}(r)=-\infty$ and that $\lim_{r\to 1}f_{s}(r)=+\infty$. If the claim is true, then the proof is concluded.

First of all we observe that

[TABLE]

This implies that $\lim_{r\to 0}H_{s}(x_{r},A)=-\infty$, and so also $\lim_{r\to 0}f_{s}(r)=-\infty$.

Moreover, recalling Lemma 2.1 we get that

[TABLE]

So, $\lim_{r\to 1}H_{s}(x_{r},A)=+\infty$, which permits to conclude that $\lim_{r\to 1}f_{s}(r)=+\infty$. ∎

We now look for more general symmetric self-shrinkers, given by the union of a finite number of annuli.

Theorem 2.3**.**

For all $N\geqslant 1$ and all $R>0$ there exists an increasing sequence $0<r_{1}<\ldots<r_{2N-1}<r_{2N}=R$, depending only on $n$, $s$ and $N$, such that such that the flow (1.1) with initial datum

[TABLE]

is a homothetically shrinking solution of (1.1).

Similarly, for all $N\geqslant 1$ and $R>0$ there exists an increasing sequence $0<\tilde{r}_{0}<\tilde{r}_{1}<\ldots<\tilde{r}_{2N-1}<r_{2N}=R$, depending only on $n$, $s$ and $N$, such that such that the flow (1.1) with initial datum

[TABLE]

is a homothetically shrinking solution of (1.1).

Proof.

The argument is similar to that in the proof of Proposition 2.2. As before, up to rescaling the sets $E$, $\widetilde{E}$, we can assume $r_{2N}=1$. Then, we want to find radii $r_{i}$ in such a way that, letting $x_{r_{i}}\in\partial B_{r_{i}}$, there hold

[TABLE]

and

[TABLE]

Notice that the functions $f_{i}$ are all continuous in their domain of definition.

We divide the proof into 4 steps. In the first step we deal with the case $N=1$ , and in step 2, 3 and 4 we consider the case $N>1$. For $N>1$ we provide the proof just of (2.9) for the existence of the set $E$, since the analogous assertion (2.10) for $\widetilde{E}$ follows similarly.

Step 1. The case $N=1$ for $E$ has been proved in Propositions 2.2. So, we consider the set $\tilde{E}$.

First of all we fix $r_{1}\in(0,1)$ and we prove that there exists $r_{0}=r_{0}(r_{1})\in(0,r_{1})$ such that $f_{0}(r_{0}(r_{1}),r_{1})=0$ for all $r_{1}\in(0,1)$. Due to the monotonicity properties of the fractional mean curvature, fixed $r_{1}\in(0,1)$ we get

[TABLE]

Moreover, by definition we get that,

[TABLE]

from which we conclude that

[TABLE]

Therefore, we obtain that

[TABLE]

By continuity of the function $f_{0}$, we deduce that for all $r_{1}\in(0,1)$ there exists at least one $r=r(r_{1})\in(0,r_{1})$ such that

[TABLE]

We choose as $r_{0}(r_{1})$ to be the smallest among all possible $r(r_{1})\in(0,r_{1})$ which solve (2.11). Observe that due to this choice the function $r\to f_{1}(r_{0}(r),r)$ is continuous. To conclude it is sufficient to prove that that there exists $r_{1}\in(0,1)$ such that $f_{1}(r_{0}(r_{1}),r_{1})=0$. Indeed, this would imply that $(B_{1}\setminus B_{r_{1}})\cup B_{r_{0}(r_{1})}$ is a solution to (1.2).

Observe that $\lim_{r\to 0}r_{0}(r)=0$, and therefore we get

[TABLE]

We now claim that

[TABLE]

Recalling Lemma 2.1, we observe that as $r\to 1$,

[TABLE]

where the constant $c=c(n,s)>0$ is given by Lemma 2.1. Similarly, we have that

[TABLE]

and

[TABLE]

Therefore as $r\to 1$ by (2.14) and (2.15)

[TABLE]

We claim that

[TABLE]

Note that the claim is equivalent to

[TABLE]

and this implies immediately, recalling (2.17), that $\lim_{r\to 1}f_{1}(r_{0}(r),r)=+\infty$.

To prove (2.18) we recall that $f_{0}(r_{0}(r),r)=0$ and using (2.14) and (2.16) we get

[TABLE]

from which we deduce that

[TABLE]

Recalling that

[TABLE]

from (2.19) we get that

[TABLE]

which gives the claim (2.18).

By continuity of $f_{1}$, from (2.12) and (2.13) it follows that there exists $r_{1}\in(0,1)$ such that $f_{1}(r_{0}(r_{1}),r_{1})=0$, which gives the thesis.

Step 2. We pass now to consider the case $N>1$. We provide a proof of the existence of a sequence of radii $r_{i}$ which solves (2.9). We shall determine $r_{i}$ by induction on $i$.

For $i=1$ we observe that, given a choice of $0<r_{2}<\ldots<r_{2N-1}<1$, we have

[TABLE]

By continuity of the function $f_{1}$ it follows that there exists $\bar{r}_{1}=\bar{r}_{1}(r_{2},\ldots r_{2N-1})\in(0,r_{2})$ such that $f_{1}(\bar{r}_{1},\ldots,r_{2N-1})=0$. As before, in case of multiple solutions we choose the smallest one. Notice that $\bar{r}_{1}$ is continuous as a function of $r_{2},\ldots,r_{2N-1}$. Notice also that, if we fix $r_{3},\ldots,r_{2N-1}$ and let $r_{2}\to r_{3}$, letting $F:=B_{\bar{r}_{1}}\cup A_{r_{3},r_{2}}$ and proceeding as in Step 1, we get

[TABLE]

Since $f_{1}(\bar{r}_{1},\ldots,r_{2N-1})=0$, we also have $H_{s}(x_{\bar{r}_{1}},E)=-\bar{r}_{1}H_{s}(x_{\bar{r}_{1}},E)=O(1)$, whence

[TABLE]

Step 3. Let now $2\leqslant i<2N-1$. By induction assumption, for all $j<i$ there exist continuous functions $\bar{r}_{j}(r_{i},\ldots,r_{2N-1}))$ such that $f_{j}(\bar{r}_{1},\ldots,\bar{r}_{i-1},r_{i},\ldots,r_{2N})=0$. In view of (2.20), we shall also assume that

[TABLE]

which is equivalent to

[TABLE]

Given a choice of $r_{j}$ for $j>i$, we want to find $\bar{r}_{i}$ such that

[TABLE]

and

[TABLE]

We first notice that

[TABLE]

We now consider the limit $r_{i}\to r_{i+1}$. Reasoning as in Step 1, we get

[TABLE]

and therefore, recalling (2.21),

[TABLE]

By continuity of $f_{i}$ it follows that there exists $\bar{r}_{i}$ such that $f_{i}(\bar{r}_{1},\ldots,\bar{r}_{i},r_{i+1},\ldots,r_{2N})=0$. As before, in case of multiple solutions we choose the smallest one.

We now show (2.23). If we fix $r_{i+2},\ldots,r_{2N-1}$ and let $r_{i+1}\to r_{i+2}$, from (2.22) we get $H_{s}(x_{\bar{r}_{i}},E)=O(1)$, which implies

[TABLE]

Multiplying by $|r_{i+1}-\bar{r}_{i}|^{s}$ and recalling (2.21) we then get

[TABLE]

which gives (2.23).

Step 4. Finally, for $i=2N-1$ we still have

[TABLE]

We now consider the limit $r_{2N-1}\to 1$. Recalling (2.23) with $i=2N-2$, as in Step 1 we get

[TABLE]

Therefore, we have

[TABLE]

As before, it follows that there exists $\bar{r}_{2N-1}$ such that $f_{2N-1}(\bar{r}_{1},\ldots,\bar{r}_{2N-1})=0$. ∎

Remark 2.4**.**

An interesting question which is left open by the previous result is the issue of uniqueness for self-shrinkers with a prescribed number of boundary spheres. In the simplest case, that is the annulus, in Proposition 2.2 we prove uniqueness of the ratio $\frac{R}{r}$ for which the annulus $B_{R}\setminus B_{r}$ is a self-similar shrinker.

From Theorem 2.3 we readily obtain the existence of cylindrical self-shrinkers.

Corollary 2.5**.**

Let $k<n$. For all $N\geqslant 1$ and $R>0$ there exists an increasing sequence $0<r_{1}<\ldots<r_{2N-1}<r_{2N}=R$, depending only on $k$, $s$ and $N$, such that such that the flow (1.1) with initial datum

[TABLE]

is a homothetically shrinking solution of (1.1), where $B^{k}_{r}$ denotes the ball of radius $r$ in $\mathbb{R}^{k}$.

Similarly, for all $N\geqslant 1$ and $R>0$ there exists an increasing sequence $0<\tilde{r}_{0}<\tilde{r}_{1}<\ldots<\tilde{r}_{2N-1}<\tilde{r}_{2N}=R$, depending only on $k$, $s$ and $N$, such that such that the flow (1.1) with initial datum

[TABLE]

is a homothetically shrinking solution of (1.1).

Remark 2.6**.**

We observe that the radii $r(n,s)$ in Proposition 2.2, $r_{i}(n,s),\tilde{r}_{i}(n,s)$ in Theorem 2.3 and Corollary 2.5 all satisfy $\lim_{s\to 1}r(n,s)=\lim_{s\to 1}r_{i}(n,s)=\lim_{s\to 1}\tilde{r}_{i}(n,s)=R$.

We give a brief justification of this fact just for the simplest case, that is the case of $r(n,s)$ in Proposition 2.2, the others being completely analogous. We recall that if $E\subset\mathbb{R}^{n}$ is a compact set with $C^{2}$ boundary then $(1-s)H_{s}(x,E)$ converges uniformly as $s\to 1$ to the classical mean curvature $H(x,\partial E)$ (see [av]). Under the same notation as in the proof of Proposition 2.2, we note that for $s=1$ the function $f_{1}(r)$ defined in (2.5) is given by $r-\frac{1}{r}$ (this is also true for the functions $f_{i}$ defined in the proof of Theorem 2.3, that is $f_{i}=r_{i}-\frac{1}{r_{i}}$). So, by uniform convergence of the curvatures, we get that if $(r_{k})_{k}$ is a sequence with $r_{k}\in(0,1)$ and $r_{k}\to 1$, there exists $(s_{k})_{k}$ with $0<s_{k}<1$ such that $f_{s}(t)<0$ for $t\in(0,r_{k}]$ and $s\geqslant s_{k}$. This implies that $r(n,s)>r_{k}$ for all $s\geqslant s_{k}$ and that $s_{k}\to 1$, since $\lim_{r\to 1}f_{s}(r)=+\infty$ for all $s<1$.

3. Stability

We now discuss the dynamic stability of the symmetric self-shrinkers constructed in the previous section. By definition self-shrinkers are stationary solutions to the flow

[TABLE]

If the initial datum is rotationally symmetric as in Theorem 2.3 then (3.1) becomes a system of ODE’s in the radii $r_{i}$, and Theorem 2.3 guarantees the existence of a stationary point for every number of radii. We are interested in the stability of such critical points, with respect to perturbations which are orthogonal to the vector $(r_{1},\ldots r_{2N})$ (or resp. $(r_{0},\ldots r_{2N})$) given by the radii. Indeed this vector corresponds to a rescaling of the initial datum, and therefore gives a direction of instability for the system which is not geometrically significant.

In the symmetric situation, we can rewrite (3.1) as the system of ODE’s

[TABLE]

Theorem 3.1**.**

Fix $N\geqslant 1$, and let $E$ (resp. $\tilde{E}$) be the symmetric shrinker given by Theorem (2.3), corresponding to the stationary point $(\bar{r}_{1},\ldots\bar{r}_{2N})$ (resp. $(\bar{r}_{0},\ldots\bar{r}_{2N})$) for the system (3.2). Then, the Morse index of such point is at least $2$, in particular the corresponding homothetic solution is dynamically unstable.

Proof.

We shall prove the assertion for the shrinker $E$, since the proof for $\tilde{E}$ is analogous.

For the reader convenience, we first present in detail the case $N=1$, corresponding to an annulus $A=B_{\bar{r}_{2}}\setminus B_{\bar{r}_{1}}$. The system (3.2) then becomes

[TABLE]

We define the function $g(r_{1},r_{2})=\left(g_{1}(r_{1},r_{2}),g_{2}(r_{1},r_{2})\right)$ as follows:

[TABLE]

We now compute the Jacobian matrix $Dg$ at the point $(\bar{r}_{1},\bar{r}_{2})$ which is a stationary point for (3.3), that is $g(\bar{r}_{1},\bar{r}_{2})=0$.

We observe the following fact: for $\delta\neq 0$, $\varepsilon>0$, $R>r>|\delta|$, there hold

[TABLE]

So, using these equalities we get that the derivative of $g_{1}$ at $(\bar{r}_{1},\bar{r}_{2})$ are given by

[TABLE]

Analogously, we observe that for $\delta\neq 0$, $R>r>|\delta|$, there holds

[TABLE]

Using this equality, we compute the derivative of $g_{2}$ at $(\bar{r}_{1},\bar{r}_{2})$:

[TABLE]

Note that, using (3.5) and (3.6),

[TABLE]

so that $(\bar{r}_{1},\bar{r}_{2})$ is an eigenvector with eigenvalue $s+1>0$. Moreover, by (3.6), we observe that $\partial_{r_{2}}g_{2}(\bar{r}_{1},\bar{r}_{2})>s+1$. This implies that

[TABLE]

which gives that $Dg(\bar{r}_{1},r_{2})$ has a second eigenvalue bigger than $s+1$, and then in particular positive.

We now consider the general case of a self-shrinker

[TABLE]

We also let $\bar{r}=(\bar{r}_{1},\ldots,\bar{r}_{2N}),\,g(\bar{r})=(g_{1}(\bar{r}),\ldots,g_{2N}(\bar{r}))\in\mathbb{R}^{2}N$, where

[TABLE]

if the index $i$ is even, and

[TABLE]

if $i$ is odd. Notice that, since $\bar{r}$ is a stationary solutions to (3.2), we have $g(\bar{r})=0$.

We compute, for $j\neq i$,

[TABLE]

and

[TABLE]

Notice that

[TABLE]

so that $\bar{r}$ is an eigenvector with eigenvalue $s+1>0$.

Now we claim that

[TABLE]

If the claim is true, then reasoning as in (3.7), we conclude that there exists an eigenvalue of $Dg(\bar{r})$ which is strictly greater than $s+1$ (and then positive), so that the Morse index of $(\bar{r}_{1},\ldots\bar{r}_{2N})$ is at least $2$.

Since

[TABLE]

to get the claim (3.9) it is sufficient to prove that for all $1\leqslant i<j<2N$ there holds

[TABLE]

We shall prove a slightly stronger statement, namely that

[TABLE]

Indeed, we compute

[TABLE]

which shows (3.11), and so proves (3.10).

∎

Remark 3.2**.**

It would be interesting to determine exactly the Morse index of the stationary points $(\bar{r}_{1},\bar{r}_{2},\dots,\bar{r}_{2N})$ (resp. $(\bar{r}_{0},\bar{r}_{1},\dots,\bar{r}_{2N})$) of the flow (3.2). In the simplest case $N=1$, we proved in Theorem 3.1 that the index of $(\bar{r}_{1},\bar{r}_{2})$ is equal to $2$.

It would also be interesting to understand if the ball is dynamically stable for any perturbation, not necessarily radial, as it happens for the standard mean curvature flow [MR772132, cm].