Home Mathematics A new energy method for shortening and straightening complete curves
Article Open Access

A new energy method for shortening and straightening complete curves

  • Tatsuya Miura ORCID logo and Fabian Rupp ORCID logo EMAIL logo
Published/Copyright: November 5, 2025

Abstract

We introduce a novel energy method that reinterprets “curve shortening” as “tangent aligning”. This conceptual shift enables the variational study of infinite-length curves evolving by the curve shortening flow, as well as higher-order flows such as the elastic flow, which involves not only the curve shortening but also the curve straightening effect. For the curve shortening flow, we prove convergence to a straight line under mild assumptions on the ends of the initial curve. For the elastic flow, we establish a global well-posedness theory, and investigate the precise long-time behavior of solutions. In fact, our method applies to a more general class of geometric evolution equations including the surface diffusion flow, Chen’s flow, and the free elastic flow.

1 Introduction

The curve shortening flow is one of the most classical geometric flows (see, e.g., [14, 4]), defined by a one-parameter family of immersed curves γ = γ ( t , x ) : [ 0 , T ) × I R n , for n 2 and an interval I R , satisfying the evolution equation

(CSF) t γ = κ ,

where κ : = s 2 γ denotes the curvature vector and s : = | x γ | 1 x is the arclength derivative along 𝛾. The flow is referred to as “curve shortening” because it arises as the L 2 ( d s ) -gradient flow of the length functional

L [ γ ] : = γ d s ,

where d s : = | x γ | d x . In particular, if the initial curve has finite length, then 𝐿 monotonically decreases along the flow, and this fact is essential in energy methods for asymptotic analysis. However, such methods clearly break down for curves of infinite length.

Identifying gradient flow structures is a fundamental technique in the analysis of parabolic PDEs. Such an underlying variational structure is however not unique – the same equation can be a gradient flow for different energies, possibly with respect to different metrics, as observed, for instance, in the seminal work of Jordan–Kinderlehrer–Otto [32]. In this spirit, this work develops a conceptually new energy method for the curve shortening flow, based on the direction energy; see (2.1) below. Our approach opens the possibility to apply energy methods to study the long-time behavior of curve-shortening-type flows for non-compact complete curves, avoiding the use of maximum principles, and thus provides a reliable and robust tool also in higher codimensions.

We first apply this method to the curve shortening flow. While the well-established theory based on maximum principles, monotonicity properties, and the classification of solitons provides some global existence results even for curves of infinite length, the asymptotic analysis remains challenging. As a direct application of our energy method, we obtain a new result on the convergence of the curve shortening flow to a straight line. Notably, this result applies to initial curves that are neither necessarily planar nor confined to a slab, which are typically difficult to handle using maximum principles.

More importantly, our method is directly applicable to higher-order flows, due to its independence from maximum principles. For higher-order flows, various energy methods have been well developed in the compact case, but the non-compact case remains significantly less explored due to the absence of a canonical finite energy.

Here we apply our direction energy method to a class of fourth-order flows of curve shortening type, including the classical surface diffusion flow [50],

(SDF) t γ = s 2 κ ,

as well as Chen’s flow introduced in [9],

(CF) t γ = s 4 γ = s 2 κ + | κ | 2 κ + 3 κ , s κ s γ ,

to obtain results analogous to those for the curve shortening flow. Here s denotes the normal derivative along 𝛾 defined by s f : = ( s f ) = s f s f , s γ s γ .

Besides the direction energy method, we also develop a new technical tool for controlling higher-order derivatives, extending the interpolation method of Dziuk–Kuwert–Schätzle [20] from the compact to the non-compact setting. This tool is not only necessary for the analysis of curve-shortening-type flows but is particularly well-suited for analyzing curve-straightening-type flows, which decrease quantities involving the bending energy

B [ γ ] : = 1 2 γ | κ | 2 d s .

In this paper, we also address the L 2 ( d s ) -gradient flow of the elastic energy B + λ L for λ 0 , the so-called 𝜆-elastic flow,

(𝜆-EF) t γ = s 2 κ 1 2 | κ | 2 κ + λ κ ,

which was introduced in [58, 20] and has since been extensively studied in the compact case (see also the survey [42]). We establish a foundational global existence theory for this flow in the non-compact case and then investigate the precise long-time behavior of solutions. In the purely straightening case λ = 0 , we show that the flow always converges to a line. On the other hand, when λ > 0 , the interaction between the curve shortening and straightening effects leads to much more complex asymptotic behavior.

Our results also yield new rigidity theorems for solitons as direct corollaries.

This paper is organized as follows. In Section 2, we provide a more detailed explanation of our main results described above, while comparing them with previous works. Section 3 presents preliminaries, particularly basic properties of the direction energy. Section 4 provides the key interpretation of the flows as gradient flows involving the direction energy. In Section 5, we prove the crucial weighted interpolation estimate, and then in Section 6 apply it to deduce global curvature bounds of any order along the flows. Using these estimates, in Section 7, we obtain the main blow-up and convergence results for each flow. Finally, Section 8 discusses rigidity results for solitons. The paper is complemented by an appendix which includes a discussion of local well-posedness in the non-compact case (Appendix A).

2 Main results

2.1 Direction energy and interpolation method

For a curve 𝛾 immersed in R n , we define the direction energy (with respect to the direction e 1 ) as

(2.1) D [ γ ] : = 1 2 γ | s γ ( s ) e 1 | 2 d s ,

where, here and in the sequel, { e j } j = 1 n denotes the canonical basis of R n . Of course, the preferred tangent direction e 1 can be replaced by any other unit vector, after suitable rotation. This energy can be equivalently represented by

D [ γ ] = γ ( 1 s γ ( s ) , e 1 ) d s ,

and hence, if L [ γ ] < , is given by

D [ γ ] = L [ γ ] γ s γ ( s ) , e 1 d s .

In particular, if 𝛾 is closed, then we simply have D = L as the last integral vanishes by integration by parts. In contrast, for non-closed curves, the two functionals may differ.

However, an important observation is that the integrand s γ ( s ) , e 1 is a null Lagrangian. As a result, the gradient flow of the direction energy still produces the curve shortening flow. This perspective offers the new geometrical interpretation

“curve shortening tangent aligning” .

Besides its conceptual interest, this idea is particularly useful in the case of non-compact complete curves, since L = , whereas 𝐷 can remain finite.

Remark 2.1

Remark 2.1 (Admissible ends)

The class of curves with finite direction energy comprises complete curves whose ends suitably converge to two parallel semi-lines with opposite directions. Not only that, this class also allows for slowly diverging ends with or without oscillation (Figure 1), such as the graphs of u ( x ) = x α or u ( x ) = x α sin log x as x , where α ( 0 , 1 2 ) , as well as spatial spiraling-out ends like u ( x ) = ( x α sin log x , x α cos log x ) . In fact, the finiteness of 𝐷 is equivalent to having graphical ends with suitably integrable gradients (see Lemma 3.2).

Figure 1 
                  An immersed planar curve with admissible ends.
Figure 1

An immersed planar curve with admissible ends.

The key idea of the present work is to replace 𝐿 with 𝐷 to produce finite monotone quantities even for non-compact curves. This idea is inspired by previous works by the first author [44] and in collaboration with Wheeler [48] on the analysis of elastica, a stationary problem involving 𝐵 and 𝐿. In those papers, the direction energy is used to variationally handle the so-called borderline elastica (Figure 2), which has infinite length and horizontal ends. Our work here is the first to apply this idea to gradient flows.

Figure 2 
                  The borderline elastica, parametrized as in (7.26).
Figure 2

The borderline elastica, parametrized as in (7.26).

After showing the energy decay, we need to handle higher-order quantities through interpolation estimates. In [20], Dziuk–Kuwert–Schätzle established the key Gagliardo–Nirenberg-type interpolation inequalities to apply energy-type methods to evolutions of closed curves. This has since been applied to a variety of flows in the compact case (even of more than fourth order, e.g., [56, 5, 38]), but does not directly extend to the non-compact case, since the inequalities involve the total length – which is infinite here. In this paper, we develop, for the first time, a variant of this interpolation method which is suitable to treat the non-compact situation. To that end, inspired by the localization procedure in [36], we first derive an interpolation estimate with a power-type weight function (Proposition 5.7). This approach is necessary, since we do not assume any a priori integrability and it is of particular importance to keep track of the relation between the exponents and the order of derivatives. We then establish higher-order estimates under localization, which are finally globalized in both time and space by a time cutoff argument in the spirit of Kuwert–Schätzle’s work [35] (Section 6).

In what follows, we discuss concrete geometric flows. Note that we only consider classical solutions in this work. However, using parabolic smoothing, our results easily transfer to appropriate weak solutions, e.g., in suitable Sobolev classes.

2.2 Curve shortening flow

The celebrated works of Gage–Hamilton [26] and Grayson [28] proved that the curve shortening flow of closed planar embedded curves must shrink to a round point in finite time.

The case of non-compact complete curves is more delicate. Following Polden [57] and Huisken [31], Chou–Zhu [13] established global-in-time existence for a wide class of complete embedded planar initial curves. Polden [57] and Huisken [31, Theorem 2.5] also proved convergence to self-similar solutions, assuming that the ends of the initial curve are asymptotic to semi-lines. The limit curve depends on the angle formed by the asymptotic semi-lines; in particular, if the angle is 𝜋, then the solution locally smoothly converges to a straight line parallel to the asymptotic semi-lines. Their proof relies on the fact that the solution remains inside an initial slab (see also Remark 7.7). Recently, Choi–Choi–Daskalopoulos [11] obtained a refined convergence to the grim reaper under convexity. See also [51, 52, 65, 66] for the planar stability of lines.

The above works strongly rely on several types of maximum principles, thus are based on planarity, since, in higher codimensions, some important maximum principles are not available; for example, embeddedness may be lost [2]. This makes the analysis much more complicated, and indeed, the authors are not aware of any relevant convergence results for infinite-length curves if n 3 .

Our main result here extends Polden and Huisken’s result on convergence to the straight line, not only to more general ends (possibly not contained in any slab) but also to general codimensions by a purely energy-based approach. We first formulate our result in the form of the following general dichotomy theorem. Hereafter, we use the notation C ̇ for the space of smooth functions with bounded derivatives, as defined in Section 3.1; this regularity assumption on the initial datum is required in our proof of local well-posedness, see Appendix A for a detailed discussion.

Theorem 2.2

Let γ : [ 0 , T ) × R R n be a curve shortening flow (CSF) with initial datum γ 0 C ̇ ( R ; R n ) such that inf R | x γ 0 | > 0 and maximal existence time T ( 0 , ] . Suppose that D [ γ 0 ] < . Then the direction energy D [ γ ( t , ) ] continuously decreases in t 0 . In addition, the following dichotomy holds.

  1. If T < , then

    lim inf t T ( T t ) 1 / 2 γ ( t , ) | κ | 2 d s > 0 .

  2. If T = , then

    lim t γ ( t , ) | κ | 2 d s = 0 ,

    and moreover,

    (2.2) lim t sup R | s γ ( t , ) e 1 | = 0 , lim t sup R | s m κ ( t , ) | = 0

    for all m N 0 . In particular, after reparametrization by arclength and translation, the solution locally smoothly converges to a horizontal line.

Remark 2.3

Here the translation is with respect to any base points; see (7.14) for the precise statement of convergence. The convergence of derivatives is global in space. In particular, this directly ensures that every global-in-time solution is eventually graphical; this is in stark contrast to the elastic flow discussed below. Even for the curve shortening flow, the curve itself may not converge uniformly (see Example 7.6). Also, translation is essentially needed, because the solution curve can keep oscillating or even escape to infinity (see Examples 7.5 and 7.6).

Remark 2.4

The L blow-up rate lim inf t T ( T t ) 1 / 2 κ ( t , ) > 0 is well known at least in the compact case and easily extended to the non-compact case, since the proof is based on pointwise estimates. Integral-type estimates are also known in the compact case; in particular, Angenent [6] obtained the sharp L p blow-up rate

lim inf t T ( T t ) ( p 1 ) / 2 γ ( t , ) | κ | p d s > 0

for all p ( 1 , ) in the case of closed planar curves. Our result seems the first integral-type blow-up estimate in the non-compact case.

Now we recall previous work providing sufficient conditions for global existence. There are two typical cases: planar embedded curves, and (possibly non-planar) graphical curves. We say that an immersed curve γ : R R n is graphical if inf R s γ , e 1 > 0 . This together with inf R | x γ | > 0 implies that 𝛾 may be represented by the graph of a function u : R R n 1 over the e 1 -axis.

Corollary 2.5

Let γ 0 C ̇ ( R ; R n ) with inf R | x γ 0 | > 0 and D [ γ 0 ] < . Suppose that γ 0 is planar and embedded (resp. γ 0 is graphical). Then the unique curve shortening flow (CSF) starting from γ 0 remains planar and embedded (resp. graphical), exists globally in time, and thus converges as in Theorem 2.2ii.

Proof

The planar embedded case directly follows by Huisken’s distance comparison argument for [31, Theorem 2.5] since D [ γ 0 ] < implies that γ 0 has appropriate graphical ends (Lemma 3.2) so that the extrinsic/intrinsic distance ratio is bounded away from zero, also near infinity. The graphical case follows since the arguments of Altschuler–Grayson [3, Theorem 1.13, Theorem 2.6 (1), (2)] for periodic graphical curves in R 3 directly extend to non-periodic graphical curves in R n (see also [30, 69]), where we need a maximum principle on the whole line as in [60, Proposition 52.4]. ∎

We can also give a new, simple blow-up criterion in the planar case.

Corollary 2.6

Let γ 0 C ̇ ( R ; R 2 ) with inf R | x γ 0 | > 0 and D [ γ 0 ] < . Suppose that γ 0 has nonzero rotation number. Then the unique curve shortening flow (CSF) starting from γ 0 blows up in finite time as in Theorem 2.2i.

Proof

Since the rotation number is preserved (see Lemma B.2), the flow keeps having a leftward tangent somewhere. If the flow exists globally, this contradicts the uniform convergence of the tangent vector in Theorem 2.2ii. ∎

2.3 Fourth-order flows of curve shortening type

As mentioned in the introduction, our argument also works even for a class of higher-order flows including the surface diffusion flow (SDF) and Chen’s flow (CF). These flows also decrease the length in the compact case, and some energy methods are already developed in previous work; see, e.g., [24, 20, 12, 67, 47] for the surface diffusion flow and [16] for Chen’s flow.

For these flows, we obtain quite analogous results to (CSF); see Theorem 7.8 for dichotomy and Corollary 7.9 for blow-up. To the authors’ knowledge, the latter is the first rigorous blow-up result for the non-compact surface diffusion as well as Chen’s flow. Although no results corresponding to Corollary 2.5 are known due to the lack of maximum principles, global existence holds at least for suitable small initial data (cf. [34]).

In addition, the above results directly imply new classification results for solitons to (CSF), (SDF), and (CF) in a unified manner. For (CSF), comprehensive classification results for solitons have been obtained in R 2 (see [29]) as well as in R n (see [1]), while for the higher-order flows, the analysis is much more complicated and such classification results are not available. However, some rigidity (i.e., partial classification) results are already known for (SDF) in the planar case [23, 61]. Here we obtain a rigidity theorem (Corollary 8.1), which is not only new in R 2 but also the first rigidity result in higher codimensions R n with n 3 .

2.4 Fourth-order flows of curve straightening type

Now we turn to flows involving the curve straightening effect, namely (𝜆-EF) with λ 0 . Up to rescaling, without loss of generality we may only consider the two cases λ = 0 , 1 .

We first address the case λ = 0 , called the free elastic flow,

(FEF) t γ = s 2 κ 1 2 | κ | 2 κ ,

which is also equivalent to the Willmore flow with initial surfaces Σ 0 invariant under translation, i.e., Σ 0 = γ 0 ( R ) × R R n + 1 . In the compact case, thanks to the decay property of the bending energy, the flow always exists globally in time only assuming finiteness of the bending energy [20], and also some asymptotic properties have recently been obtained [68, 49]. Here we establish the non-compact counterpart for the first time. The limit curve is again a straight line, but the direction is not determined due to the absence of the tangent aligning effect.

Theorem 2.7

Let γ 0 C ̇ ( R ; R n ) be a properly immersed curve with inf R | x γ 0 | > 0 . Suppose that B [ γ 0 ] < . Then there exists a unique, smooth, properly immersed, global-in-time solution γ : [ 0 , ) × R R n to the free elastic flow (FEF) with initial datum γ 0 . Along the flow, the energy B [ γ ( t , ) ] continuously decreases in t 0 .

In addition,

lim t sup R | s m κ ( t , ) | = 0

holds for all m N 0 . In particular, after reparametrization by arclength, translation, and rotation, the solution locally smoothly converges to a straight line.

Finally, we discuss the case λ = 1 , which we simply call the elastic flow,

(EF) t γ = s 2 κ 1 2 | κ | 2 κ + κ .

There are a number of known results about global existence and convergence in the compact case; see, e.g., [58, 20, 18, 17, 27, 39, 43, 54, 53, 19, 40].

To the authors’ knowledge, the only known result for the infinite-length elastic flow is a pioneering work by Novaga–Okabe [53]. Roughly speaking, they show that if the initial curve is planar and if the ends are asymptotic to the e 1 -axis and represented by graphs of class W m , 2 for all m 0 , then there exists a global-in-time elastic flow, and the solution locally smoothly sub-converges to either a line or a borderline elastica (Figure 2). Here and in the sequel, sub-convergence means that, for any time sequence t j , there exists a subsequence t j along which the flow converges (in a suitable sense). Their result does not include any uniqueness because their construction is based on approximating an infinite-length solution by finite-length solutions using the Arzelà–Ascoli theorem.

In this paper, we apply our energy method to improve the result of Novaga–Okabe in several aspects: we prove uniqueness, extend the analysis to general codimension, cover more general ends (possibly unbounded), and ensure the horizontality of the convergence limits. To this end, we apply our direction energy method and regard the elastic flow as the gradient flow of the adapted energy

E [ γ ] : = B [ γ ] + D [ γ ] = γ ( 1 2 | κ | 2 + 1 2 | s γ e 1 | 2 ) d s .

The finiteness of 𝐸 is equivalent to having graphical ends 𝑢 with u W 1 , 2 (see Lemma 3.2), so the diverging examples as in Remark 2.1 are still admissible.

Our main result can be summarized as follows.

Theorem 2.8

Theorem 2.8 (Theorems 7.107.12)

Let γ 0 C ̇ ( R ; R n ) with inf R | x γ 0 | > 0 . Suppose that E [ γ 0 ] < . Then there exists a unique, smooth, properly immersed, global-in-time solution γ : [ 0 , ) × R R n to the elastic flow (EF) with initial datum γ 0 . Along the flow, the energy E [ γ ( t , ) ] continuously decreases in t 0 .

In addition, after any reparametrization by arclength and some translation, the solution sub-converges to an elastica γ locally smoothly. The limit curve γ is either a straight line or a (planar) borderline elastica, and in each case, the curve γ has asymptotically horizontal ends: s γ ( s ) e 1 as s ± .

Remark 2.9

In stark contrast to the previous flows, global-in-space convergence does not hold, in general. Indeed, we can construct a peculiar example with loops “escaping” to infinity (see Example 7.17 and Figure 3).

Figure 3 
                  Two loops moving apart along the elastic flow.
Figure 3

Two loops moving apart along the elastic flow.

Furthermore, we find several conditions for initial data such that we can identify the type of the limit elastica, which in particular implies full (locally smooth) convergence without taking a time subsequence. Among other results (see Corollary 7.15), we obtain the optimal energy threshold below which the same convergence holds as in the case of the curve shortening flow.

Corollary 2.10

Let γ 0 C ̇ ( R ; R n ) with inf R | x γ 0 | > 0 . Suppose that

E [ γ 0 ] < 8 .

Then the elastic flow (EF) starting from γ 0 satisfies

lim t sup R | s γ ( t , ) e 1 | = 0 , lim t sup R | s m κ ( t , ) | = 0

for all m N 0 . In particular, after reparametrization by arclength and translation, the solution locally smoothly converges to a horizontal line.

Remark 2.11

The energy threshold is optimal. Indeed, a borderline elastica γ 0 with E [ γ 0 ] = 8 is a stationary solution.

We finally mention that, although the curve shortening flow preserves many positivity properties such as planar embeddedness or graphicality (cf. Corollary 2.5), a wide class of higher-order flows does not due to the lack of maximum principles [10, 25]. Recently, optimal energy thresholds for preserving embeddedness of elastic flows of closed curves have been obtained in [46]. Our new energy method is also useful in this direction; we can obtain optimal thresholds in terms of 𝐸 for preserving planar embeddedness (resp. graphicality) in the non-compact case. This will be addressed in our future work.

3 Preliminaries

3.1 Notation

Let 𝐍 denote the set of positive integers and N 0 : = N { 0 } . For k N 0 and α ( 0 , 1 ] , we denote by C k , α ( R ; R n ) the usual space of 𝑘-times differentiable functions u : R R n with finite norm

u C k , α ( R ; R n ) : = i = 0 k x i u + [ x i u ] α ,

where [ ] α denotes the 𝛼-Hölder seminorm if α < 1 and the best Lipschitz constant for α = 1 . Since we work with non-compact curves, we consider, for k N and α ( 0 , 1 ] , the function spaces

C ̇ k , α ( R ; R n ) : = { u : R R n x u C k 1 , α ( R ; R n ) } .

In similar fashion, we define C ̇ k ( R ; R n ) for k 1 . We also set

C ̇ ( R ; R n ) : = k = 1 C ̇ k ( R ; R n ) .

If the codomain is clear from the context, we also just write C k , α ( R ) or C ̇ k , α ( R ) .

3.2 Basic properties of the direction energy

The bending energy 𝐵 is invariant with respect to isometries of R n , and also any reparametrization. On the other hand, the direction energy 𝐷 has less invariances; translation, rotation around the e 1 -axis, reflection in directions orthogonal to e 1 , and orientation-preserving reparametrization.

Lemma 3.1

Let γ : R R n be a smooth immersion with inf R | x γ | > 0 and D [ γ ] < . Then

lim x ± γ ( x ) , e 1 = ± ,

and hence 𝛾 is proper. If in addition γ C ̇ 1 , 1 ( R ; R n ) , then

lim x ± s γ ( x ) = e 1 .

Proof

The first assertion follows by taking the limits x ± in the identity

0 x | x γ ( x ) | d x γ ( x ) , e 1 = 0 x ( | x γ ( x ) | x γ ( x ) , e 1 ) d x γ ( 0 ) , e 1

since the right-hand side is bounded thanks to the assumption D [ γ ] < .

The second assertion follows since the integrability of | x γ | x γ , e 1 and its Lipschitz continuity imply that | x γ | x γ , e 1 0 as x ± . ∎

Lemma 3.2

Let γ C ̇ 1 , 1 ( R ; R n ) be a smooth immersion with inf R | x γ | > 0 .

  1. We have D [ γ ] < if and only if outside a slab S R : = { p R n R p , e 1 R } the curve 𝛾 is represented by a graph curve

    { ( x , u ( x ) ) R n x R }

    of u C ̇ 1 , 1 ( R ; R n 1 ) with u L 2 ( R ; R n 1 ) .

  2. We have B [ γ ] + D [ γ ] < if and only if outside a slab S R the curve 𝛾 is represented by a graph curve of u C ̇ 1 , 1 ( R ; R n 1 ) with u W 1 , 2 ( R ; R n 1 ) .

Proof

We first prove the “only if” part of (i). Suppose D [ γ ] < . By Lemma 3.1, we deduce that, outside a slab S R , the curve 𝛾 is a graph curve with bounded gradient; more precisely, there are R > 0 and u C ̇ 1 , 1 ( R ; R n 1 ) such that

γ ( R ) S R = { ( x , u ( x ) ) R n x R } S R .

Since D [ γ ] < , and since the energy inside the slab S R does not affect the finiteness, the graph of 𝑢 also has finite direction energy D ̃ [ u ] < , where

D ̃ [ u ] : = R ( 1 1 1 + | u | 2 ) 1 + | u | 2 d x = R | u | 2 1 + 1 + | u | 2 d x c u L 2 2 ,

with c > 0 depending only on u . This completes the “only if” part of (i). The “if” part is easier since D ̃ [ u ] 1 2 u L 2 2 < , implying that D [ γ ] < . This completes the proof of (i).

Now we turn to (ii). Suppose B [ γ ] + D [ γ ] < . Let u C ̇ 1 , 1 be as above; note that u ′′ L . Then, in addition to u L 2 obtained above, we have the finiteness of the bending energy B ̃ [ u ] < , where

B ̃ [ u ] : = R | u ′′ | 2 2 ( 1 + | u | 2 ) 3 1 + | u | 2 d x c u ′′ L 2 2 .

This implies the “only if” part of (ii). The “if” part is again easier since B ̃ [ u ] 1 2 u ′′ L 2 2 < . The proof is now complete. ∎

4 Gradient flow structures

In this section, we obtain key energy-decay properties for the flows under consideration, through a cutoff argument.

From this section onward, we will primarily use γ : ( t , x ) γ ( t , x ) to denote a flow of curves. In this case, we sometimes abbreviate γ ( t , ) as γ ( t ) . Furthermore, we write f d s to denote an integral of the following form:

f d s : = γ ( t , ) f ( t , ) d s = R f ( t , x ) | x γ ( t , x ) | d x .

Finally, we define the norm

f L p ( d s ) : = ( | f | p d s ) 1 / p ,

which we may also express as f L p ( d s ( t ) ) to emphasize the time dependence.

Throughout this section, we use 𝜉 and 𝑉 to denote the tangential and normal velocity of a flow 𝛾, respectively; namely,

ξ : = t γ , s γ , V : = t γ = t γ ξ s γ .

4.1 Evolution of localized energy functionals

We first compute general time-derivative formulae for the direction and bending energy.

Lemma 4.1

Let T ( 0 , ) and γ : [ 0 , T ] × R R n be a smooth family of immersed curves. Let η C c ( R ) and define D [ γ | η ] : = ( 1 s γ , e 1 ) η d s . Then

d d t D [ γ | η ] + κ , V η d s = ξ ( 1 s γ , e 1 ) s η d s + V , e 1 s η d s .

Proof

Standard geometric evolution formulae (cf. [20, Lemma 2.1]) yield

(4.1) t ( d s ) = ( s ξ κ , V ) d s , t ( s γ ) = s V + ξ κ .

We thus compute

t [ ( 1 s γ , e 1 ) d s ] = ( s V + ξ κ , e 1 + ( 1 s γ , e 1 ) ( s ξ κ , V ) ) d s .

Now we differentiate D [ γ | η ] . By the regularity assumption, the integral is finite, and we may differentiate under the integral. Hence we have

d d t D [ γ | η ] = ( s V + ξ κ , e 1 + ( 1 s γ , e 1 ) ( s ξ κ , V ) ) η d s .

After rearranging,

d d t D [ γ | η ] + κ , V η d s = s ξ ( 1 s γ , e 1 ) η d s s V , e 1 η d s + ( ξ κ , e 1 + s γ , e 1 κ , V ) η d s .

Integrating by parts in the first and second term of the right-hand side, and in particular using s ( e 1 ) = s γ , e 1 κ for the second term, we obtain the desired formula. ∎

Lemma 4.2

Let T ( 0 , ) and let γ : [ 0 , T ] × R R n be a smooth family of immersed curves. Let η C c ( R ) and define

B [ γ | η ] : = 1 2 | κ | 2 η d s .

Then

d d t B [ γ | η ] + s 2 κ 1 2 | κ | 2 κ , V η d s = 1 2 | κ | 2 ξ s η d s + 2 s κ , V s η d s + κ , V s 2 η d s .

Proof

In addition to (4.1), recalling from [20, Lemma 2.1] that

(4.2) t κ = s 2 V + κ , V κ + ξ s κ ,

we argue along the lines of the proof of Lemma 4.1 to obtain

d d t B [ γ | η ] = ( κ , s 2 V + 1 2 | κ | 2 κ , V + ξ κ , s κ + 1 2 | κ | 2 s ξ ) η d s .

Integrating by parts in the first term yields

κ , s 2 V η d s = s κ , s V η d s κ , s V s η d s = s 2 κ , V η d s + 2 s κ , V s η d s + κ , V s 2 η d s ,

so that, after some rearranging,

d d t B [ γ | η ] + s 2 κ 1 2 | κ | 2 κ , V η d s = s ( 1 2 | κ | 2 ξ ) η d s + 2 s κ , V s η d s + κ , V s 2 η d s .

Integration by parts for the first term yields the assertion. ∎

4.2 Energy identity

Now we turn to the proof of the energy identity. We first address the flows having the structure of an L 2 -gradient flow, namely (CSF) and (𝜆-EF). Given σ , λ 0 , we define the functional

F [ γ ] = F σ , λ [ γ ] : = σ B [ γ ] + λ D [ γ ] = γ ( σ 2 | κ | 2 + λ ( 1 s γ , e 1 ) ) d s .

Proposition 4.3

Let σ , λ 0 , T ( 0 , ) , and let γ : [ 0 , T ] × R R n be a smooth solution to

t γ = σ ( s 2 κ + 1 2 | κ | 2 κ ) + λ κ

such that

(4.3) inf [ 0 , T ] × R | x γ | > 0 , k = 1 M x k γ L ( [ 0 , T ] × R ) + t γ L ( [ 0 , T ] × R ) < ,

where M = 4 if σ > 0 , while M = 2 if σ = 0 . If F [ γ ( 0 ) ] < , then

F [ γ ( T ) ] F [ γ ( 0 ) ] = 0 T | t γ | 2 d s d t .

Proof

Suppose σ > 0 . Let η C c ( R ) and define D [ γ | η ] and B [ γ | η ] as above. Combining Lemmas 4.1 and 4.2, we find

(4.4) d d t F [ γ | η ] + | V | 2 η d s = σ 1 2 | κ | 2 ξ s η d s + 2 σ s κ , V s η d s + σ κ , V s 2 η d s λ ξ ( 1 s γ , e 1 ) s η d s + λ V , e 1 s η d s .

Take a cutoff function φ 0 C c ( R ) with

χ [ 1 / 2 , 1 / 2 ] φ 0 χ [ 1 , 1 ] ,

let R > 1 , define φ ( x ) : = φ 0 ( R 1 x ) , and set η = φ 4 . In the following, C ( 0 , ) is a constant depending on λ , σ , φ 0 , and the bounds on 𝛾 in (4.3), that may change from line to line. Then, by (4.3),

(4.5) | s η | C R 1 φ 3 , | s 2 η | C R 1 φ 2 ,

and also

(4.6) R R d s C R .

Using (4.3), we can bound the integrands on the right-hand side of (4.4) by

C ( | s η | + | s 2 η | ) d s .

Since this is bounded by 𝐶 thanks to (4.5) and (4.6), we conclude

d d t F [ γ | φ 4 ] + | V | 2 φ 4 d s C .

Integrating over [ 0 , t ] [ 0 , T ] and sending R , we conclude from the monotone convergence theorem that F [ γ ] = lim R F [ γ | φ 4 ] < for all t [ 0 , T ] with

F [ γ ( t ) ] + 0 t | V | 2 d s d τ F [ γ ( 0 ) ] + C T .

In particular, it follows that

(4.7) t F [ γ ( t ) ] L ( 0 , T ) and t V L 2 ( d s ( t ) ) L 2 ( 0 , T ) .

We have η 1 a.e. in ( 0 , T ) × R as R and, by (4.5) and (4.6),

(4.8) s η L 2 ( d s ( t ) ) 0 , s 2 η L 2 ( d s ( t ) ) 0 in L 2 ( 0 , T ) .

Therefore, integrating (4.4) in time and sending R imply the assertion; indeed, by (4.3), the absolute value of the right-hand side is bounded by

C ξ s η 0 T ( σ 2 | κ | 2 + λ ( 1 s γ , e 1 ) ) d s d t + C ( s κ + 1 ) 0 T | V | | s η | d s d t + C κ 0 T | V | | s 2 η | d s d t C s η 0 T F [ γ ] d t + C 0 T | V | ( | s η | + | s 2 η | ) d s d t ,

which converges to zero thanks to (4.5), (4.6), (4.7), and (4.8).

Finally, if σ = 0 , the exact same argument applies after removing all terms involving 𝜎, which results in only terms involving derivatives up to order 2. Thus the statement follows. ∎

Our method also works for gradient flows which are not of L 2 -type. We now address fourth-order flows of curve shortening type including (SDF) and (CF).

Proposition 4.4

Let T ( 0 , ) and let γ : [ 0 , T ] × R R n be a smooth solution to

t γ = s 2 κ + μ | κ | 2 κ

for some constant μ 0 such that

inf [ 0 , T ] × R | x γ | > 0 , k = 1 4 x k γ L ( [ 0 , T ] × R ) + t γ L ( [ 0 , T ] × R ) < .

If D [ γ ( 0 ) ] < , then

D [ γ ( T ) ] D [ γ ( 0 ) ] = 0 T ( | s κ | 2 + μ | κ | 4 ) d s d t .

Proof

We use Lemma 4.1 with the same cutoff function η = φ 4 as in the proof of Proposition 4.3. Integration by parts for terms involving s 2 κ yields

(4.9) d d t D [ γ | η ] + ( | s κ | 2 + μ | κ | 4 ) η d s = ξ ( 1 s γ , e 1 ) s η d s + μ | κ | 2 κ , e 1 s η d s + s κ , e 1 s 2 η d s s κ , κ s γ , e 1 s η d s s κ , κ s η d s ,

where we again used s ( e 1 ) = s γ , e 1 κ . As before, bounding the right-hand side by 𝐶, integrating in time over [ 0 , t ] [ 0 , T ] , and sending R , we obtain

D [ γ ( t ) ] + 0 t ( | s κ | 2 + μ | κ | 4 ) d s D [ γ ( 0 ) ] + C T ,

which implies

t D [ γ ( t ) ] L ( 0 , T ) , t s κ L 2 ( d s ( t ) ) 2 + μ κ L 4 ( d s ( t ) ) 4 L 1 ( 0 , T ) .

Integrating (4.9) over [ 0 , T ] , the right-hand side may be bounded by

ξ s η 0 T D [ γ ] d t + 0 T ( μ | κ | 3 | s η | + | s κ | | s 2 η | + 2 | κ | | s κ | | s η | ) d s d t .

Using Hölder’s inequality and the fact that, by (4.5), for any p > 1 , we have

s η L p ( d s ( t ) ) + s 2 η L p ( d s ( t ) ) 0 in L p ( 0 , T )

as R yield the claim. ∎

5 Localized interpolation estimates

The goal of this section is to prove the fundamental Gagliardo–Nirenberg-type estimate, Proposition 5.7 below. Throughout this section, we fix a smooth immersion γ : R R n and a cutoff function

(5.1) ζ C c ( R ) with ζ 0 and | s ζ | Λ ,

where the arclength derivative is along 𝛾, i.e., s ζ = | x γ | 1 x ζ , and Λ > 0 .

We begin with establishing L 2 -type interpolation estimates with the weight 𝜁. Our approach is inspired by the interpolation estimates in [18, 20] and the spatial cutoff in [36].

Lemma 5.1

Let N with 2 . There exists C = C ( Λ , ) > 0 such that, for any smooth vector field 𝑋 normal along 𝛾 and any ε ( 0 , 1 ] , we have

| s X | 2 ζ d s ε | s 2 X | 2 ζ + 2 d s + C ε [ ζ > 0 ] | X | 2 ζ 2 d s .

Proof

By integration by parts and Young’s inequality,

| s X | 2 ζ d s = s 2 X , X ζ d s s X , X ζ 1 s ζ d s ε 2 | s 2 X | 2 ζ + 2 d s + C ε [ ζ > 0 ] | X | 2 ζ 2 d s + 1 2 | s X | 2 ζ d s + C ( Λ , ) [ ζ > 0 ] | X | 2 ζ 2 d s .

The claim follows after absorbing and using ε ( 0 , 1 ] . ∎

Lemma 5.2

Let k N with k 2 , and let i N with i < k . There exists

C = C ( Λ , i , k ) > 0

such that, for any smooth vector field 𝑋 normal along 𝛾 and any ε ( 0 , 1 ] , we have

| s i X | 2 ζ 2 i d s ε | s k X | 2 ζ 2 k d s + C ε i i k [ ζ > 0 ] | X | 2 d s .

Proof

We proceed by induction on 𝑘. First, for k = 2 and i = 1 , this is Lemma 5.1 with = 2 .

Suppose the statement is true for k 2 . Let 0 < i < k + 1 . In the following,

C = C ( Λ , i , k ) > 1

is a constant that is allowed to change from line to line.

We first consider the case i < k . Then, by the induction hypothesis, we obtain

| s k 1 X | 2 ζ 2 k 2 d s δ | s k X | 2 ζ 2 k d s + C δ 1 k [ ζ > 0 ] | X | 2 d s

for all δ ( 0 , 1 ] . On the other hand, Lemma 5.1 (with = 2 k ) yields

(5.2) | s k X | 2 ζ 2 k d s η 2 | s k + 1 X | 2 ζ 2 k + 2 d s + C η | s k 1 X | 2 ζ 2 k 2 d s η 2 | s k + 1 X | 2 ζ 2 k + 2 d s + C δ η | s k X | 2 ζ 2 k d s + C δ 1 k η [ ζ > 0 ] | X | 2 d s

for any η ( 0 , 1 ] . Choosing δ = η 2 C < 1 in the above equation, we may absorb the second term. Now, again, by the induction hypothesis, for any ε ( 0 , 1 ] , we have

| s i X | 2 ζ 2 i d s ε | s k X | 2 ζ 2 k d s + C ε i i k [ ζ > 0 ] | X | 2 d s , ε η | s k + 1 X | 2 ζ 2 k + 2 d s + C ε δ k [ ζ > 0 ] | X | 2 d s + C ε i i k [ ζ > 0 ] | X | 2 d s ,

plugging in (5.2) and using the choice of 𝛿. We now take

η = ε 1 k i ( 0 , 1 ] and ε ̃ = ε k + 1 i k i ( 0 , 1 ] ,

and observe that ε η = ε ̃ , that ε δ k = C ε η k = C ε i i k , and that ε i i k = ε ̃ i i ( k + 1 ) . This yields the statement (for 𝑘 replaced by k + 1 ) in the case 0 < i < k for 𝜀 replaced by ε ̃ .

It remains to consider the case i = k . In this case, for any ε , δ ( 0 , 1 ] , by Lemma 5.1 (with = 2 k ) and the induction hypothesis, we have

| s k X | 2 ζ 2 k d s ε 2 | s k + 1 X | 2 ζ 2 k + 2 d s + C ε | s k 1 X | 2 ζ 2 k 2 d s ε 2 | s k + 1 X | 2 ζ 2 k + 2 d s + C ε ( δ | s k X | 2 ζ 2 k d s + C δ 1 k [ ζ > 0 ] | X | 2 d s ) .

Taking δ = ε 2 C < 1 here and absorbing, we thus obtain

| s k X | 2 ζ 2 k d s ε | s k + 1 X | 2 ζ 2 k + 2 d s + C ε k [ ζ > 0 ] | X | 2 d s .

Since ε k k ( k + 1 ) = ε k , the case i = k is proven and the statement follows. ∎

We now establish a multiplicative version. For 𝑋 normal along 𝛾, we define the seminorms (using the weight 𝜁)

| X | 0 , 2 := X L 2 ( [ ζ > 0 ] , d s ) , | X | k , 2 := ζ k s k X L 2 ( d s ) + | X | 0 , 2 ( k 1 ) .

Lemma 5.3

Let k N and i N 0 with i k . There exists C = C ( Λ , i , k ) > 0 such that, for any smooth vector field 𝑋 normal along 𝛾, we have

ζ i s i X L 2 ( d s ) C | X | k , 2 i k | X | 0 , 2 k i k .

Proof

Without loss of generality, we may assume k 2 and 0 < i < k ; otherwise, there is nothing to show. Lemma 5.2 implies that, for all ε ( 0 , 1 ] , we have

ζ i s i X L 2 ( d s ) C ( ε | X | k , 2 + ε i i k | X | 0 , 2 ) .

If | X | 0 , 2 = 0 , the statement is trivial. Otherwise, choosing

ε = ( | X | k , 2 | X | 0 , 2 ) i k k ( 0 , 1 ] ,

the statement follows. ∎

This will enable us to estimate nonlinearities arising in higher-order energy evolutions. In order to apply L 2 -type estimates for L p -type nonlinearities, we now prove the following key Gagliardo–Nirenberg-type interpolation inequality. The crucial observation is that this estimate distributes the weight in a way that is compatible with the seminorms | X | k , 2 , | X | 0 , 2 defined above.

Lemma 5.4

Let r 0 and p [ 2 , ] . There exists C = C ( Λ , p , r ) > 0 such that, for any smooth normal vector field 𝑋 along 𝛾, we have

(5.3) ( | X | p ζ p ( r + θ ) d s ) 1 p C ( | s X | 2 ζ 2 r + 2 d s ) θ 2 ( | X | 2 ζ 2 r d s ) 1 θ 2 + C ( | X | 2 ζ 2 r d s ) 1 2

for p [ 2 , ) , where θ = p 2 2 p [ 0 , 1 2 ) , and, for p = ,

ζ r + 1 2 X C ( | s X | 2 ζ 2 r + 2 d s ) 1 4 ( | X | 2 ζ 2 r d s ) 1 4 + C ( | X | 2 ζ 2 r d s ) 1 2 .

Remark 5.5

Defining, for p [ 2 , ] , the seminorm

| X | 0 , p : = ζ 1 2 1 p X L p ( [ ζ > 0 ] , d s ) ( = ζ θ X L p ( [ ζ > 0 ] , d s ) ) ,

with the interpretation | X | 0 , : = ζ 1 2 X , Lemma 5.4 resembles the classical Gagliardo–Nirenberg inequality, since (5.3) yields

| ζ r X | 0 , p C | ζ r X | 1 , 2 θ | ζ r X | 0 , 2 1 θ .

By taking 𝛾 in Lemma 5.4 to be an arclength parametrization of a straight line in R 2 , we also obtain a non-geometric version for functions on the real line. Although we will not use it here, we include the statement for the reader’s convenience.

Corollary 5.6

Let ξ C c ( R ) with ξ 0 and | ξ | Λ . Let r 0 and p [ 2 , ] . There exists C = C ( Λ , p , r ) > 0 such that, for all u W loc 1 , 2 ( R ) , we have

( | u | p ξ p ( r + θ ) d x ) 1 p C ( | u | 2 ξ 2 r + 2 d x ) θ 2 ( | u | 2 ξ 2 r d x ) 1 θ 2 + C ( | u | 2 ξ 2 r d x ) 1 2

for p [ 2 , ) , where θ = p 2 2 p , and for p = ,

ξ r + 1 2 u C ( | u | 2 ξ 2 r + 2 d x ) 1 4 ( | u | 2 ξ 2 r d x ) 1 4 + C ( | u | 2 ξ 2 r d x ) 1 2 .

Proof of Lemma 5.4

If p = 2 , there is nothing to prove. Let p > 2 . Since 𝑋 is normal, we have | s | X | | | s X | a.e. Using that 𝜁 has compact support, we have

| X | 2 ζ 2 r + 1 ( 2 | X | | s X | ζ 2 r + 1 + ( 2 r + 1 ) | X | 2 ζ 2 r | s ζ | ) d s 2 ( | X | 2 ζ 2 r d s ) 1 2 ( | s X | 2 ζ 2 r + 2 d s ) 1 2 + ( 2 r + 1 ) Λ | X | 2 ζ 2 r d s .

This yields the statement for p = . For p ( 2 , ) , we further estimate

| X | p ζ p ( r + 1 2 ) 1 d s | X | p 2 ζ ( r + 1 2 ) ( p 2 ) | X | 2 ζ 2 r d s = | X | 2 ζ 2 ( r + 1 2 ) p 2 2 | X | 2 ζ 2 r d s C ( | X | 2 ζ 2 r d s ) p + 2 4 ( | s X | 2 ζ 2 r + 2 d s ) p 2 4 + C ( | X | 2 ζ 2 r d s ) p 2 .

The statement follows. ∎

Following the notation in [18] (cf. [20]) for a , b , c N 0 , we denote by P b a , c any linear combination (with universal coefficients) of terms of the form s i 1 κ s i b κ that may be scalar or normal vector fields along 𝛾. Here i j c for all j = 1 , , b and j = 1 b i j = a .

Proposition 5.7

Let k N , a , b , c N 0 , b 2 , and suppose c k , a + b 2 1 < 2 k . Then, for any ε > 0 , there exists C = C ( ε , Λ , a , b , k ) such that

| P b a , c | ζ a + b 2 1 d s ε | s k κ | 2 ζ 2 k d s + C ( [ ζ > 0 ] | κ | 2 d s ) b δ 2 δ + C ( [ ζ > 0 ] | κ | 2 d s ) b 2

with δ = 1 k ( a + b 2 1 ) < 2 .

Proof

If suffices to consider the case P b a , c = s i 1 κ s i b κ . Again, we denote by 𝐶 a constant that is allowed to vary from line to line, depending only on ε , Λ , a , b , k .

We first prove

(5.4) | P b a , c | ζ a + b 2 1 d s C | κ | k , 2 δ | κ | 0 , 2 b δ + C | κ | k , 2 a k | κ | 0 , 2 b a k .

We may assume | κ | 0 , 2 > 0 ; otherwise, the statement is trivial.

For proving (5.4), we first assume c k 1 . Using Hölder’s inequality and

a + b 2 1 = j = 1 b ( i j + 1 2 1 b ) ,

we obtain

(5.5) | s i 1 κ s i b κ | ζ a + b 2 1 d s j = 1 b ζ i j + 1 2 1 b s i j κ L b ( d s ) .

We now estimate the factors in (5.5) by using the Gagliardo–Nirenberg-type inequality. Indeed, Lemma 5.4 with r = i j yields

( | s i j κ | b ζ b i j + b 2 1 d s ) 1 b C ( | s i j + 1 κ | 2 ζ 2 i j + 2 ) θ 2 ( | s i j κ | 2 ζ 2 i j d s ) 1 θ 2 + C ( | s i j κ | 2 ζ 2 i j d s ) 1 2 ,

where θ = b 2 2 b . Applying Lemma 5.3 and combining exponents, we find

( | s i j κ | b ζ b i j + b 2 1 d s ) 1 b C | κ | k , 2 i j + θ k | κ | 0 , 2 k i j θ k + C | κ | k , 2 i j k | κ | 0 , 2 k i j k = C | κ | k , 2 i j k | κ | 0 , 2 k i j k ( | κ | k , 2 θ k | κ | 0 , 2 θ k + 1 ) ,

using that | κ | 0 , 2 > 0 . Plugging this into (5.5), we obtain

| s i 1 κ s i b κ | ζ a + b d s C ( j = 1 b | κ | k , 2 i j k | κ | 0 , 2 k i j k ) ( | κ | k , 2 θ k | κ | 0 , 2 θ k + 1 ) b C | κ | k , 2 a k | κ | 0 , 2 b k a k ( | κ | k , 2 b θ k | κ | 0 , 2 b θ k + 1 ) = C | κ | k , 2 a + b θ k | κ | 0 , 2 b k a b θ k + C | κ | k , 2 a k | κ | 0 , 2 b k a k .

Noting that 1 k ( a + b θ ) = 1 k ( a + b 2 1 ) = δ , we obtain (5.5) and thus (5.4) when c k 1 .

On the other hand, suppose c = k . Then P b a , c contains a term of the form s k κ , so we have k a a + b 2 1 < 2 k (as b 2 ), and we may write P b a , c = P b 1 a k , k 1 s k κ . The Cauchy–Schwarz inequality yields

| P b a , c ζ a + b 2 1 | d s ( | s k κ | 2 ζ 2 k d s ) 1 2 ( | P b ̃ a ̃ , k 1 | ζ a ̃ + b ̃ 2 1 d s ) 1 2 ,

where a ̃ : = 2 ( a k ) , b ̃ : = 2 ( b 1 ) . Since a ̃ + b ̃ 2 1 = 2 ( a + b 2 1 k ) < 2 k and as we have already proven (5.4) in the case c k 1 , we may apply (5.4) to the second factor, yielding

| P b a , c ζ a + b 2 1 | d s C | κ | k , 2 ( | κ | k , 2 δ ̃ | κ | 0 , 2 b ̃ δ ̃ + C | κ | k , 2 a ̃ k | κ | 0 , 2 b ̃ a ̃ k ) 1 2 .

As δ = 1 + δ ̃ 2 , 1 2 ( b ̃ δ ̃ ) = b δ , 1 + a ̃ 2 k = a k , and 1 2 ( b ̃ a ̃ k ) = b a k , we have thus proven (5.4) also in the case c = k .

Lastly, we estimate the first term of the right-hand side of (5.4) using Young’s inequality and δ [ 0 , 2 ) , yielding

| κ | k , 2 δ | κ | 0 , 2 b δ C ( | s k κ | 2 ζ 2 k d s ) δ 2 ( [ ζ > 0 ] | κ | 2 d s ) b δ 2 + C ( [ ζ > 0 ] | κ | 2 d s ) b 2 ε | s k κ | 2 ζ 2 k d s + C ( [ ζ > 0 ] | κ | 2 d s ) b δ 2 δ + C ( [ ζ > 0 ] | κ | 2 d s ) b 2 .

Using the same argument with 𝛿 replaced by a k , we additionally estimate the second term in (5.4) and get the term

( [ ζ > 0 ] | κ | 2 d s ) b a k 2 a k { ( [ ζ > 0 ] | κ | 2 d s ) b 2 if [ ζ > 0 ] | κ | 2 d s 1 , ( [ ζ > 0 ] | κ | 2 d s ) b δ 2 δ if [ ζ > 0 ] | κ | 2 d s 1 ,

where we used that x b x 2 x is monotonically increasing and the estimate 0 a k δ in the last step. Thus (5.4) implies the statement. ∎

6 Global curvature control

In this section, we prove global curvature estimates for a flow with uniformly bounded bending energy 𝐵. Having established the key interpolation estimate, Proposition 5.7, we now consider a rather general evolution law and present a unified approach to treat all the flows we consider in this work simultaneously by a combination of the strategies in [20, 35].

6.1 Evolution of localized curvature integrals

We consider a general class of geometric flows in this section, including (CSF), (SDF), (CF), (𝜆-EF), i.e., we assume that

γ : [ 0 , T ) × R R n

is a smooth family of proper immersions, satisfying

(6.1) t γ = σ s 2 κ + λ κ + μ | κ | 2 κ + ϑ κ , s κ s γ .

Here λ , μ , σ , ϑ R are fixed parameters. In the following key lemma, in abuse of notation, we write, for α , β R ,

( α + β ) ( P b a , c + P b a , c ) : = α P b a , c + β P b a , c + α P b a , c + β P b a , c .

Lemma 6.1

Let ζ ̃ C c ( R n ) and set ζ ( t , ) : = ζ ̃ γ ( t , ) . For all m N 0 and r N with r 2 , we have

d d t 1 2 | s m κ | 2 ζ 2 m + r d s + σ | s m + 2 κ | 2 ζ 2 m + r d s + λ | s m + 1 κ | 2 ζ 2 m + r d s = ( ( σ + μ + ϑ ) ( P 4 2 m + 2 , m + 2 + P 6 2 m , m ) + λ P 4 2 m , m ) ζ 2 m + r d s + ( 2 m + r ) σ P 3 2 m + 2 , m + 2 + λ P 3 2 m , m , D ζ ̃ γ ζ 2 m + r 1 d s + ( 2 m + r ) μ P 5 2 m , m + ϑ P 4 2 m + 1 , m + 1 s γ , D ζ ̃ γ ζ 2 m + r 1 d s 2 σ ( 2 m + r ) s m + 2 κ , s m + 1 κ ζ 2 m + r 1 s ζ d s σ ( 2 m + r ) s m + 2 κ , s m κ ζ 2 m + r 2 ( ( 2 m + r 1 ) ( s ζ ) 2 + ζ s 2 ζ ) d s λ ( 2 m + r ) s m + 1 κ , s m κ ζ 2 m + r 1 s ζ d s .

Proof

We first show that the derivatives of the curvature satisfy

(6.2) t s m κ + σ s m + 4 κ λ s m + 2 κ = ( σ + μ + ϑ ) ( P 3 m + 2 , m + 2 + P 5 m , m ) + λ P 3 m , m .

Indeed, for m = 0 , this follows by (4.2), where

V : = t γ = σ s 2 κ + λ κ + μ P 3 0 , 0 and ξ : = t γ , s γ = ϑ P 2 1 , 1 .

For m 1 , recall from [20, (2.8)] that we have

( t s s t ) s m κ = κ V s m + 1 κ + s ξ s m + 1 κ + κ s m κ s V ,

so (6.2) follows by induction on m N 0 , using s ξ = ϑ P 2 2 , 2 .

We now compute the time derivative of the localized integral. Note that, since γ ( t ) is proper by assumption, ζ ( t , ) has compact support for all t [ 0 , T ) . By (4.1), we have

(6.3) d d t 1 2 | s m κ | 2 ζ 2 m + r d s = s m κ , t s m κ ζ 2 m + r d s + 2 m + r 2 | s m κ | 2 ζ 2 m + r 1 D ζ ̃ γ , t γ d s + 1 2 | s m κ | 2 ζ 2 m + r ( s ξ κ , V ) d s .

By (6.2), the leading order term arising from the first term satisfies

(6.4) s m + 4 κ , s m κ ζ 2 m + r d s = s m + 3 κ , s m + 1 κ ζ 2 m + r d s ( 2 m + r ) s m + 3 κ , s m κ ζ 2 m + r 1 s ζ d s = | s m + 2 κ | 2 ζ 2 m + r d s + 2 ( 2 m + r ) s m + 2 κ , s m + 1 κ ζ 2 m + r 1 s ζ d s + ( 2 m + r ) s m + 2 κ , s m κ ζ 2 m + r 2 × ( ( 2 m + r 1 ) ( s ζ ) 2 + ζ s 2 ζ ) d s ,

using integration by parts twice. Similarly, by integration by parts,

(6.5) s m + 2 κ , s m κ ζ 2 m + r d s = | s m + 1 κ | 2 ζ 2 m + r d s ( 2 m + r ) s m + 1 κ , s m κ ζ 2 m + r 1 s ζ d s .

Moreover, for the second term, we have

(6.6) | s m κ | 2 ζ 2 m + r 1 D ζ ̃ γ , t γ d s = σ P 3 2 m + 2 , max { m , 2 } + λ P 3 2 m , m , D ζ ̃ γ ζ 2 m + r 1 d s + μ P 5 2 m , m + ϑ P 4 2 m + 1 , max { m , 1 } s γ , D ζ ̃ γ ζ 2 m + r 1 d s .

Using that s ξ = P 2 2 , 2 , the third term is

(6.7) | s m κ | 2 ζ 2 m + r ( s ξ κ , V ) d s = ( σ P 4 2 m + 2 , max { m , 2 } + λ P 4 2 m , m + μ P 6 2 m , m + ϑ P 4 2 m + 2 , max { m , 2 } ) ζ 2 m + r d s .

Inserting (6.2) into (6.3), using (6.4), (6.5), (6.6), (6.7), and noting that the terms of the form P b a , max { c , c } can also be regarded as P b a , c + c , we conclude the desired identity. ∎

6.2 Curvature control for fourth-order flows

Now fix ζ ̃ C c ( R n ) with 0 ζ ̃ 1 , | D ζ ̃ | Λ , and | D 2 ζ ̃ | Λ 2 . If γ C ( R ; R n ) is a proper immersion, then

ζ : = ζ ̃ γ C c ( R )

satisfies

(6.8) ζ C c ( R ) , 0 ζ 1 , | s ζ | Λ , and | s 2 ζ | Λ 2 + | κ | Λ .

In particular, 𝜁 satisfies the assumptions in (5.1), so that Proposition 5.7 is applicable. Then, for the flow in (6.1), we obtain the following key energy estimate.

Lemma 6.2

Let σ > 0 . Suppose that B [ γ ( t ) ] M holds at some time t [ 0 , T ) . Then, at the time 𝑡,

d d t | s m κ | 2 ζ 2 m + 4 d s + | s m κ | 2 ζ 2 m + 4 d s + σ | s m + 2 κ | 2 ζ 2 m + 4 d s C

holds for all m N 0 , where C = C ( λ , μ , σ , ϑ , Λ , m , M ) .

Proof

We estimate the terms on the right-hand side of the identity in Lemma 6.1 for r = 4 . Using (6.8) for the derivatives of 𝜁, there is C = C ( λ , μ , σ , ϑ , Λ , m ) such that

d d t 1 2 | s m κ | 2 ζ 2 m + 4 d s + σ | s m + 2 κ | 2 ζ 2 m + 4 d s C | P b a , c | ζ q d s ,

where the sum runs over the finitely many terms with c m + 2 and a + b 2 1 < 2 ( m + 2 ) as well as q a + b 2 1 . Noting also 0 ζ 1 , which yields 0 ζ q ζ a + b 2 1 , we may thus apply Proposition 5.7 with k = m + 2 to conclude that

d d t 1 2 | s m κ | 2 ζ 2 m + 4 d s + 2 σ 3 | s m + 2 κ | 2 ζ 2 m + 4 d s C ,

where C = C ( λ , μ , σ , ϑ , Λ , m , M ) . Since Proposition 5.7 also applies to

| s m κ | 2 ζ 2 m + 4 d s = P 2 2 m , m ζ 2 m + 4 d s ,

the statement follows. ∎

A smooth initial datum with finite bending energy 𝐵 only needs to satisfy κ L 2 ( d s ) , so following the argument of [35, Theorem 3.5], we obtain instantaneous higher Sobolev integrability for positive time.

Lemma 6.3

Let σ > 0 and let T T < . Suppose B [ γ ( t ) ] M for all t [ 0 , T ) . Then, for all m N and t ( 0 , T ) ,

| s m κ | 2 ζ 2 m + 4 d s C ( λ , μ , σ , ϑ , Λ , m , M , T ) t m 2 .

Proof

Fix any t 0 ( 0 , T ) . Define Lipschitz cutoff functions in time ξ j : R [ 0 , 1 ] via

(6.9) ξ j ( t ) : = { 0 for t ( j 1 ) t 0 m , m t 0 ( t ( j 1 ) t 0 m ) for ( j 1 ) t 0 m t j t 0 m , 1 for t j t 0 m ,

where 0 j m . We also define ξ 1 ( t ) : = 0 . We note that

(6.10) 0 d d t ξ j m t 0 ξ j 1 for all 0 j m and a.e. t 0 .

We further define

E j ( t ) : = | s 2 j κ | 2 ζ 4 j + 4 d s , e j ( t ) : = ξ j ( t ) E j ( t ) .

Now Lemma 6.2, (6.10), and 0 ζ 1 imply that, for 0 j m and a.e. t [ 0 , T ) ,

(6.11) d d t e j + σ ξ j E j + 1 m t 0 ξ j 1 E j + C ( λ , μ , σ , ϑ , Λ , j , M ) .

We now show that, for all 0 j m and t [ 0 , t 0 ] , we have

(6.12) e j ( t ) + σ 0 t ξ j E j + 1 d τ C ( λ , μ , σ , ϑ , Λ , j , m , M , T ) t 0 j .

We prove (6.12) by induction on 𝑗. For j = 0 , (6.11) yields

e 0 ( t ) + σ 0 t ξ 0 E 1 d τ C ( λ , μ , σ , ϑ , Λ , M ) t + e 0 ( 0 ) C ( λ , μ , σ , ϑ , Λ , M , T ) .

Let j 1 . Then, integrating (6.11) and using ξ j ( 0 ) = 0 (for j 1 ), we obtain

e j ( t ) + σ 0 t ξ j E j + 1 d τ 0 t m t 0 ξ j 1 E j d τ + C ( λ , μ , σ , ϑ , Λ , j , M ) t 0 j + 1 t 0 j m t 0 C ( λ , μ , σ , ϑ , Λ , j , m , M , T ) t 0 j 1 + C ( λ , μ , σ , ϑ , Λ , j , M , T ) t 0 j ,

where we also used the induction hypothesis in the last step. This proves (6.12). Evaluating (6.12) at t = t 0 for j = m and using ξ m ( t 0 ) = 1 , we conclude

(6.13) | s 2 m κ | 2 ζ 4 m + 4 d s | t = t 0 C ( λ , μ , σ , ϑ , Λ , m , M , T ) t 0 m .

To get the same for 2 m replaced by 2 m + 1 , we argue as in Lemma 5.1. By integration by parts, Hölder, and Young’s inequality,

| s 2 m + 1 κ | 2 ζ 4 m + 6 d s ( | s 2 m + 2 κ | 2 ζ 4 m + 8 d s ) 1 2 ( | s 2 m κ | 2 ζ 4 m + 4 d s ) 1 2 + 1 2 | s 2 m + 1 κ | ζ 4 m + 6 d s + C ( Λ , m ) | s 2 m κ | ζ 4 m + 4 d s .

Absorbing, evaluating at t = t 0 , and using (6.13), we find

(6.14) | s 2 m + 1 κ | 2 ζ 4 m + 6 d s | t = t 0 C ( λ , μ , σ , ϑ , Λ , m , T ) t 0 m + 1 2 t 0 m 2 + C ( λ , μ , σ , ϑ , Λ , m , T ) t 0 m C ( λ , μ , σ , ϑ , Λ , m , T ) t 0 2 m + 1 2 .

Renaming t 0 into 𝑡, the statement now follows from (6.13) and (6.14). ∎

Lemma 6.4

Let σ > 0 . Suppose that B [ γ ( t ) ] M for all t [ 0 , T ) . Then, for all m N 0 and t ( 0 , T ) ,

(6.15) s m κ L 2 ( d s ) C ( λ , μ , σ , ϑ , m , M ) ( 1 + t m 4 ) ,
(6.16) s m κ C ( λ , μ , σ , ϑ , m , M ) ( 1 + t 2 m + 1 8 ) .

Proof

We first prove (6.15). For any R > 0 , there exists ζ ̃ C c ( R n ) with

χ B R ( 0 ) ζ ̃ χ B 2 R ( 0 ) ,

and such that | D ζ ̃ | c R 1 and | D 2 ζ ̃ | c R 2 for some universal c > 0 , in particular independent of 𝑛. Hence, taking R 1 sufficiently large (depending on the universal c > 0 ), we may thus assume Λ = 1 for ζ = ζ ̃ γ in (6.8) in the sequel.

If T 1 or t 1 T , then the statement follows directly from Lemma 6.3 choosing T = 1 , sending R , and using the monotone convergence theorem.

Suppose T > 1 and t 1 . Lemma 6.3 with R yields

(6.17) | s m κ | 2 d s | t = 1 C ( λ , μ , σ , ϑ , m , M ) .

For 𝑅 sufficiently large, Lemma 6.2 and Gronwall’s inequality imply

| s m κ | 2 ζ 2 m + 4 d s C ( λ , μ , σ , ϑ , m , M ) + | s m κ | 2 ζ 2 m + 4 d s | t = 1 e 1 t C ( λ , μ , σ , ϑ , m , M ) ( 1 + t 2 m 4 ) ,

using (6.17) and t 1 . Sending R and using a + b a + b for a , b 0 give (6.15).

For (6.16), by Lemma 5.4 (with Λ = 1 ), we have

ζ m + 1 2 s m κ C ( m ) ( | s m + 1 κ | 2 ζ 2 m + 2 d s ) 1 4 ( | s m κ | 2 ζ 2 m d s ) 1 4 + C ( m ) ( | s m κ | 2 ζ 2 m d s ) 1 2 C ( λ , μ , σ , ϑ , m , M ) [ ( 1 + t m + 1 4 ) 1 2 ( 1 + t m 4 ) 1 2 + 1 + t m 4 ] ,

where we used (6.15) in the last step. Sending R yields (6.16). ∎

We next establish bounds on the derivatives s m κ (without taking the normal projection).

Lemma 6.5

Under the assumptions of Lemma 6.4, for m N 0 and t ( 0 , T ) ,

(6.18) s m κ C ( λ , μ , σ , ϑ , m , M ) ( 1 + t 2 m + 1 8 ) ,
(6.19) s m κ L 2 ( d s ) C ( λ , μ , σ , ϑ , m , M ) ( 1 + t m 4 ) .

Proof

We recall from [20, Lemma 2.6] that

(6.20) s m κ s m κ = i = 1 [ m 2 ] Q 2 i + 1 m 2 i ( κ ) + i = 1 [ m + 1 2 ] Q 2 i m + 1 2 i ( κ ) s γ .

Here [ ] is the floor function, and Q b a ( κ ) is a linear combination of terms of the form

s i 1 κ s i b κ , with j = 1 b i j = a .

Estimate (6.18) thus follows by induction from (6.16) and (6.20). For (6.19), one may also proceed by induction using (6.15), (6.18), and (6.20), as well as | κ | 2 d s M . ∎

6.3 Curvature control for the curve shortening flow

We now consider the case σ = μ = ϑ = 0 and λ > 0 in (6.1), yielding essentially the curve shortening flow.

Lemma 6.6

Let σ = μ = ϑ = 0 and λ > 0 . Suppose that B [ γ ( t ) ] M holds at some time t [ 0 , T ) . Then, at the time 𝑡,

d d t | s m κ | 2 ζ 2 m + 2 d s + | s m κ | 2 ζ 2 m + 2 d s + λ | s m + 1 κ | 2 ζ 2 m + 2 d s C

holds for all m N 0 , where C = C ( λ , Λ , m , M ) .

Proof

For this choice of parameters, Lemma 6.1 with r = 2 reads

d d t 1 2 | s m κ | 2 ζ 2 m + 2 d s + λ | s m + 1 κ | 2 ζ 2 m + 2 d s = λ P 4 2 m , m ζ 2 m + 2 d s + ( 2 m + 2 ) λ P 3 2 m , m , D ζ ̃ γ ζ 2 m + 1 d s λ ( 2 m + 2 ) s m κ , s m + 1 κ ζ 2 m + 1 s ζ d s ,

where C = C ( λ , Λ , m ) . As in the proof of Lemma 6.2, the right-hand side consists of terms of the form | P b a , c | ζ q d s , where the parameters satisfy a + b 2 1 < 2 ( m + 1 ) , q a + b 2 1 . Hence Proposition 5.7 with k = m + 1 applies and, after absorbing, we have

d d t 1 2 | s m κ | 2 ζ 2 m + 2 d s + 2 λ 3 | s m + 1 κ | 2 ζ 2 m + 2 d s C .

Applying Proposition 5.7 (also with k = m + 1 ) to

| s m κ | 2 ζ 2 m + 2 d s = P 2 2 m , m ζ 2 m + 2 d s

yields the statement. ∎

We have the following analogue of Lemma 6.3.

Lemma 6.7

Let σ = μ = ϑ = 0 and λ > 0 , and let 0 < T T < . Suppose that B [ γ ( t ) ] M for all t [ 0 , T ) . Then, for all m N and t ( 0 , T ) ,

| s m κ | 2 ζ 2 m + 2 d s C ( λ , Λ , m , M , T ) t m .

Proof

We proceed exactly as in the proof of Lemma 6.3, redefining E j by

E j ( t ) : = | s j κ | 2 ζ 2 j + 2 d s , e j ( t ) : = ξ j ( t ) E j ( t ) ,

where ξ j is as in (6.9). Now, we may use Lemma 6.6 to conclude that (6.11) holds with σ > 0 replaced by λ > 0 . We may thus prove (6.12) by induction and evaluate at t = t 0 and for j = m to conclude the statement. ∎

Arguing exactly as in Lemmas 6.4 and 6.5, we obtain the following.

Lemma 6.8

Let σ = μ = ϑ = 0 and λ > 0 . Suppose B [ γ ( t ) ] M for all t [ 0 , T ) . Then, for all m N 0 and t ( 0 , T ) ,

s m κ L 2 ( d s ) + s m κ L 2 ( d s ) C ( λ , Λ , m , M ) ( 1 + t m 2 ) , s m κ + s m κ C ( λ , Λ , m , M ) ( 1 + t 2 m + 1 4 ) .

7 Long-time behavior

In this section, we discuss the long-time behavior, proving the main theorems in Section 2. Throughout this section, we will repeatedly use, without explicitly mentioning, the fact that, under the finiteness of the direction energy, we always have properness (Lemma 3.1) and hence we can use the results from Section 6.

7.1 Blow-up rates

Consider a general flow (6.1) with σ > 0 . Given γ 0 C ̇ ( R ) with inf R | x γ 0 | > 0 , we define the maximal existence time T ( 0 , ] by the supremum of T ( 0 , ) such that there is a unique smooth solution γ : [ 0 , T ] × R R n to (6.1) starting from γ 0 such that γ γ 0 C ( [ 0 , T ] × R ) and inf [ 0 , T ] × R | x γ | > 0 . The maximal existence time is well-defined thanks to Theorem A.3. We call a solution with maximal existence time 𝑇maximal solution.

Of course, not all flows as in (6.1) will develop a singularity. However, our curvature estimates enable us to characterize the behavior of the bending energy whenever this happens in finite time.

Theorem 7.1

Suppose σ > 0 . Let γ 0 C ̇ ( R ) be a properly immersed curve with inf R | x γ 0 | > 0 and B [ γ 0 ] < . Suppose that the maximal solution γ : [ 0 , T ) × R R n to (6.1) with initial datum γ 0 exists for a finite time T < . Then we have

lim inf t T ( T t ) 1 / 4 B [ γ ( t ) ] > 0 .

To prove this, we first show the continuity of the bending energy along the flow.

Lemma 7.2

Let γ : [ 0 , T ) × R R n be a properly immersed maximal solution to (6.1) with σ > 0 such that B [ γ 0 ] < . Then the map t B [ γ ( t ) ] is finite-valued, continuous on [ 0 , T ) , and locally Lipschitz continuous on ( 0 , T ) .

Proof

Let T [ 0 , T ) . Let R > 1 and let η = φ 4 with φ C c ( R ) as in the proof of Proposition 4.3. For all t [ 0 , T ] , Lemma 4.2 yields

(7.1) d d t B [ γ | η ] + σ | s 2 κ | 2 η d s = λ s 2 κ , κ η d s + μ | κ | 2 s 2 κ , κ η d s + 1 2 | κ | 2 κ , V η d s 1 2 | κ | 2 ξ s η d s + 2 s κ , V s η d s + κ , V s 2 η d s .

Estimating

(7.2) | λ s 2 κ , κ η d s | σ 2 | s 2 κ | 2 η d s + λ 2 2 σ | κ | 2 η d s

and rearranging, we deduce

d d t B [ γ | η ] + σ 2 | s 2 κ | 2 η d s 1 2 | κ | 2 ( λ 2 σ 1 + 2 μ s 2 κ , κ + κ , V ) η d s 1 2 | κ | 2 ξ s η d s + 2 s κ , V s η d s + κ , V s 2 η d s .

Using that the derivatives of 𝛾 are uniformly bounded by C = C ( γ , T ) (independently of 𝑅) on [ 0 , T ] × R , and that, by (4.5) and (4.6), we have ( | s η | + | s 2 η | ) d s C , we estimate

(7.3) d d t B [ γ | η ] + σ 2 | s 2 κ | 2 η d s C ( γ , λ , μ , σ , T ) B [ γ | η ] + C ( γ , T ) .

Dropping the term involving σ > 0 , applying Gronwall’s lemma to the smooth function

t B [ γ ( t ) | η ] ,

and sending R yield

sup [ 0 , T ] B [ γ ] < .

We now integrate (7.1) in time over [ t 1 , t 2 ] [ 0 , T ] , rearrange so that the left-hand side is only B [ γ ( t 2 ) | η ] B [ γ ( t 1 ) | η ] , take absolute values of both sides, use (7.2), and send R to get

(7.4) | B [ γ ( t 2 ) ] B [ γ ( t 1 ) ] | t 1 t 2 ( 3 σ 2 | s 2 κ | 2 d s + C B [ γ ( t ) ] + C ) d t 3 σ 2 t 1 t 2 | s 2 κ | 2 d s d t + C | t 2 t 1 | ,

where C = C ( γ , λ , μ , σ , T ) . The L 2 -bounds in Lemma 6.4 imply that the first term on the right is controlled if [ t 1 , t 2 ] [ t 0 , T ] ( 0 , T ] , in which case we conclude the Lipschitz estimate

| B [ γ ( t 2 ) ] B [ γ ( t 1 ) ] | C ( γ , λ , μ , σ , t 0 , T ) | t 2 t 1 | .

On the other hand, integrating (7.3) on [ 0 , T ] , we have

B [ γ ( T ) ] B [ γ ( 0 ) ] + σ 2 0 T | s 2 κ | 2 d s d t C T ,

which implies that 0 T | s 2 κ | 2 d s d t < . Hence, taking t 1 = 0 and t 2 = t 0 in (7.4), we conclude continuity at t = 0 . ∎

Proof of Theorem 7.1

We first note that (although the finiteness of the bending energy does not imply properness), since here we assume the properness of γ 0 , the maximal solution 𝛾 is also proper on [ 0 , T ) by Theorem A.3.

By Lemma 7.2, we have B [ γ ( t ) ] < for all t [ 0 , T ) . We first prove that

lim sup t T B [ γ ( t ) ] = .

Indeed, if this was not the case, then there exists M ( 0 , ) such that B [ γ ( t ) ] M for all t [ 0 , T ) . By Lemma 6.5, for all m N 0 , we have

sup t [ T / 2 , T ) s m κ C ( λ , μ , σ , ϑ , m , M , T ) .

From (6.1), we conclude that, for t [ T / 2 , T ) , we have

s m t γ C ( λ , μ , σ , ϑ , m , M , T ) ,

and further deduce

sup t [ T / 2 , T ) γ ( t ) γ ( T / 2 ) L ( R ) T 2 sup t [ T / 2 , T ) t γ L ( R ) C ( λ , μ , σ , ϑ , M , T ) .

In addition, by the differential equation t | x γ | = κ , V | x γ | and the boundedness of κ , V , on [ T / 2 , T ) , we also have the uniform estimate

(7.5) C 1 | x γ | C , where  C = C ( x γ 0 , λ , μ , σ , ϑ , M , T ) > 0 .

Differentiation and induction as in the proof of [20, Theorem 3.1] yield that, for t [ T / 2 , T ) and all m N 0 , we have

x m | x γ | C ( x γ 0 , λ , μ , σ , ϑ , m , M , T ) .

Combining the above estimates, we finally deduce that, for all t [ T / 2 , T ) and m N 0 ,

x m ( γ γ ( T / 2 ) ) C ( x γ 0 , λ , μ , σ , ϑ , m , M , T ) .

For any T / 2 t 1 t 2 < T , we have

sup x R | x m γ ( t 1 , x ) x m γ ( t 2 , x ) | sup x R max t [ t 1 , t 2 ] | t x m γ ( t , x ) | | t 1 t 2 | sup t [ 0 , T ) x m t γ | t 1 t 2 | C ( x γ 0 , λ , μ , σ , ϑ , m , M , T ) | t 1 t 2 | .

This together with (7.5) implies that γ ( t ) γ ( T / 2 ) is Cauchy in C m ( R ) as t T and thus there exists a smooth immersion γ ( T ) C ̇ ( R ) with inf R | x γ ( T ) | > 0 such that

γ ( t ) γ ( T / 2 ) γ ( T ) γ ( T / 2 )

as t T in C m ( R ) for all m N 0 . Now we may apply Theorem A.3 with initial datum γ ( T ) and glue the resulting solution together with γ : [ 0 , T ] × R R n to extend 𝛾 beyond 𝑇, contradicting the maximality of 𝑇.

Hence there exist t j T with B [ γ ( t j ) ] . By Lemma 6.1 with m = 0 and r = 4 , we have

d d t 1 2 | κ | 2 ζ 4 d s + σ | s 2 κ | 2 ζ 4 d s | P b a , c | ζ q d s ,

where the sum is finite. The indices satisfy c 2 and the triple ( a , b , q ) runs over the set of ( 2 , 4 , 4 ) , ( 0 , 6 , 4 ) , ( 0 , 4 , 4 ) , ( 2 , 3 , 3 ) , ( 0 , 3 , 3 ) , ( 0 , 5 , 3 ) , ( 1 , 4 , 3 ) , ( 3 , 2 , 3 ) , ( 2 , 2 , 2 ) , ( 2 , 3 , 3 ) , ( 1 , 2 , 3 ) . Applying Proposition 5.7 (with k = 2 ), we may estimate

(7.6) d d t 1 2 | κ | 2 ζ 4 d s + σ 2 | s 2 κ | 2 ζ 4 d s C ( B [ γ ] b δ 2 δ + B [ γ ] b 2 ) ,

where δ = 1 2 ( a + b 2 1 ) . Using the inequality x p C ( p , p 0 , p 1 ) ( x p 0 + x p 1 ) for all x 0 and 0 < p 0 p p 1 < , we may estimate the right-hand side of (7.6) by finding the smallest and largest exponent, respectively. We thus obtain

d d t 1 2 | κ | 2 ζ 4 d s + σ 2 | s 2 κ | 2 ζ 4 d s C ( B [ γ ] + B [ γ ] 5 ) .

Integrating and taking η 1 , we obtain, for all 0 t 1 t 2 < T ,

(7.7) B [ γ ( t 2 ) ] B [ γ ( t 1 ) ] + C t 1 t 2 ( B [ γ ( τ ) ] + B [ γ ( τ ) ] 5 ) d τ .

The statement then follows from Lemma B.1i. ∎

In the case of the curve shortening flow ( σ = μ = ϑ = 0 and λ = 1 ), the situation is similar but algebraically less involved. The maximal existence time is defined analogously as in the case σ > 0 . The following result is a part of Theorem 2.2.

Lemma 7.3

Let γ 0 C ̇ ( R ) with inf R | x γ 0 | > 0 and D [ γ 0 ] < . Suppose that the maximal solution γ : [ 0 , T ) × R R n to (CSF) with initial datum γ 0 exists for a finite time T < . Then

lim inf t T ( T t ) 1 / 2 B [ γ ( t ) ] > 0 .

As in the fourth-order case, we first show the continuity of B [ γ ( t ) ] . Here the initial time t = 0 has to be excluded since the initial bending energy may be infinite.

Lemma 7.4

Let γ : [ 0 , T ) × R R n be a maximal solution to (CSF) with D [ γ 0 ] < . Then the map ( 0 , T ) t B [ γ ( t ) ] is finite-valued and locally Lipschitz continuous.

Proof

Using Lemma 4.2 with V = κ and ξ = 0 , and integrating by parts, we obtain

(7.8) d d t B [ γ | η ] + | s κ | 2 η d s = 1 2 | κ | 4 η d s s κ , κ s η d s .

Let T ( 0 , T ) . Let R > 1 and η = φ 4 with φ C c ( R ) as in the proof of Proposition 4.3.

Dropping the second term on the left-hand side of (7.8) and using that the derivatives of 𝛾 are bounded by C = C ( γ , T ) on [ 0 , T ] × R , and that, by (4.5) and (4.6), we have | s η | d s C , we deduce that, for all t [ 0 , T ] ,

(7.9) d d t B [ γ | η ] C B [ γ | η ] + C .

Integrating over [ t 1 , t 2 ] [ 0 , T ] and sending R yield, for all 0 t 1 t 2 T ,

(7.10) B [ γ ( t 2 ) ] B [ γ ( t 1 ) ] + C t 1 t 2 ( B [ γ ( t ) ] + 1 ) d t .

At this moment, B [ γ ( t 2 ) ] and B [ γ ( t 1 ) ] are possibly infinite, but by Proposition 4.3, the last integral term is always finite. By this integrability, we have B [ γ ( t 1 ) ] < for a.e. t 1 [ 0 , T ] , and for any such t 1 , we deduce from (7.10) that B [ γ ( t 2 ) ] < for all t 2 [ t 1 , T ] . This implies that B [ γ ( t ) ] < holds for all t ( 0 , T ] .

Now going back to (7.9), using Gronwall’s lemma for the smooth function t B [ γ ( t ) | η ] and sending R yield, for any t 0 ( 0 , T ] ,

sup t [ t 0 , T ] B [ γ ( t ) ] C ( γ , T , t 0 ) .

This allows us to use the L 2 - and L -bounds in Lemma 6.8 on the time interval [ t 0 , T ] . Thanks to this, integrating (7.8) over any interval [ t 1 , t 2 ] [ t 0 , T ] , sending R , and using (4.5), we obtain

(7.11) B [ γ ( t 1 ) ] B [ γ ( t 2 ) ] + t 1 t 2 | s κ | 2 d s d t = 1 2 t 1 t 2 | κ | 4 d s d t ,

and in particular the desired Lipschitz continuity

| B [ γ ( t 1 ) ] B [ γ ( t 2 ) ] | C ( γ , t 0 , T ) | t 1 t 2 | .

Proof of Lemma 7.3

We first establish an integral inequality for the bending energy (independent of the finiteness of 𝑇). With R > 1 and φ C c ( R ) as in the proof of Proposition 4.3, Proposition 5.7 implies

1 2 | κ | 4 φ 2 d s | s κ | 2 φ 4 d s + C 0 ( ( R R | κ | 2 d s ) 3 + ( R R | κ | 2 d s ) 2 )

for some universal C 0 > 0 . Sending R and inserting the resulting estimate into (7.11), we conclude that, for a.e. 0 < t 1 t 2 < T , we have

(7.12) B [ γ ( t 2 , ) ] B [ γ ( t 1 , ) ] + C 0 t 1 t 2 ( B [ γ ( t ) ] 3 + B [ γ ( t ) ] 2 ) d t .

Now, suppose that T < . If lim sup t T B [ γ ( t ) ] < , then using Lemma 6.8 and proceeding exactly as in Theorem 7.1, we can extend the flow past 𝑇, contradicting maximality. Hence lim sup t T B [ γ ( t ) ] = , and the assertion follows from (7.12) and Lemma B.1i. ∎

Hereafter, we examine precise convergence properties for the flows discussed in Section 2. Since in general the parametrization speed | x γ | of the solution is hard to control as t , we work with the arclength reparametrization as usual. Note carefully that, in our non-compact case, the choice of the base point of the reparametrization is important, because depending on the choice the limit may change.

For a smooth family of immersions

γ : [ 0 , T ) × R R n with L [ γ ( t ) ] = for all t [ 0 , T ) ,

we define the canonical reparametrization γ ̃ : [ 0 , T ) × R R n of 𝛾 by

γ ̃ ( t , s ) : = γ ( t , Φ t ( s ) ) ,

where Φ t : R R denotes the inverse map of x 0 x | x γ ( t , x ) | d x , meaning that the base point of the reparametrization is x = 0 . By construction, Φ t gives an orientation-preserving unit-speed reparametrization of 𝛾. In particular, both the direction energy and the bending energy are invariant under this reparametrization.

7.2 Convergence for the curve shortening flow

We first discuss convergence for the curve shortening flow, completing the proof of Theorem 2.2.

Proof of Theorem 2.2

We have already shown the blow-up rate (Lemma 7.3) and the decreasing property of the direction energy (Proposition 4.3). Hereafter, we show the remaining convergence properties, so suppose that T = .

We first prove that B [ γ ( t ) ] 0 . By Proposition 4.3, we have B [ γ ] L 1 ( 0 , ) . Hence

(7.13) lim sup j inf t [ j , j + 1 ] B [ γ ( t ) ] lim j j j + 1 B [ γ ( t ) ] d t = 0 ,

and thus there exists t j [ j , j + 1 ] such that B [ γ ( t j ) ] 0 as j . Moreover, by Lemma 7.4, B [ γ ] is locally Lipschitz continuous on ( 0 , ) . The desired convergence now follows from (7.12) and Lemma B.1ii.

For smooth convergence, Lemma 7.4 and the above convergence B [ γ ( t ) ] 0 yield sup t [ 1 , ) B [ γ ( t ) ] < , and hence we can use the bounds in Lemma 6.8 as t . Fix any family of base points { s t } t 0 R and take any time sequence t j . By Lemma 6.8 and a standard compactness argument – such as the Arzelà–Ascoli theorem applied in { | s | R } for all m 0 and a diagonal argument – we deduce that there is a time subsequence t j such that the sequence of translated arclength reparametrizations

γ ̃ ( t j , + s t j ) γ ̃ ( t j , s t j )

locally smoothly converges to an arclength parametrized curve

γ C ̇ ( R ; R n ) with γ ( 0 ) = 0 .

By Fatou’s lemma and B [ γ ( t j ) ] 0 , we have B [ γ ] = 0 , i.e., γ is a straight line. Also by Fatou’s lemma and Proposition 4.3, we deduce that

D [ γ ] lim inf j D [ γ ( t j ) ] D [ γ 0 ] < ,

so the line is uniquely given by γ ( s ) = s e 1 . The uniqueness of the limit implies that

(7.14) lim t ( γ ̃ ( t , + s t ) γ ̃ ( t , s t ) ) = γ in C loc m ( R ) for all m N 0 .

Moreover, since { s t } is arbitrary, we have shown (2.2) for any m N 0 . ∎

In order to complement our result, we mention two nontrivial examples. Both are given in the framework of planar graphs, and thus we use the fact that a planar graph remains so under the curve shortening flow.

Example 7.5

Example 7.5 (Oscillation)

There exists a graphical planar initial curve

γ 0 ( x ) = ( x , u 0 ( x ) ) with u 0 C ( R )

and D [ γ 0 ] < from which the curve shortening flow 𝛾 keeps oscillating. More precisely, if x t R denotes the unique point such that γ ( t , x t ) , e 1 = 0 , then

lim inf t γ ( t , x t ) , e 2 = 1 , lim sup t γ ( t , x t ) , e 2 = 1 .

This follows by Nara–Taniguchi’s result [52, Theorem 1.4] on a certain equivalence of the asymptotic behaviors of solutions to the heat equation and to the graphical curve shortening equation, combined with oscillating solutions to the heat equation by Collet–Eckmann [15, Lemma 8.6]. The function u 0 can be taken to have first derivative in L p ( R ) for all p > 1 , thus having finite direction energy by Lemma 3.2.

In this case, it is also easy to show that γ ( t ) cannot uniformly converge to a line, even after translation.

Example 7.6

Example 7.6 (Escaping to infinity)

There exists an unbounded graphical planar initial curve γ 0 with D [ γ 0 ] < from which the curve shortening flow 𝛾 remains unbounded at each time, and escapes as time goes to infinity, in the sense that

lim x ± γ ( t , x ) , e 2 = for any t 0 , lim t inf x R γ ( t , x ) , e 2 = .

In fact, we may take any embedded curve γ 0 C ̇ ( R ; R 2 ) with graphical ends given by | x | α as x ± , where α ( 0 , 1 2 ) , so that D [ γ 0 ] < by Lemma 3.2. Then we can argue by comparing the solution with suitable barriers.

For the first assertion, given any t 0 > 0 and any M > 0 , we site a right-going grim reaper in a slab { p R 2 M 1 p , e 2 M } , sufficiently far from the e 2 -axis (depending on γ 0 , t 0 , M ) so that, for all t 0 , the grim reaper lies below γ 0 . Then the comparison between those two curves implies that lim inf x γ ( t 0 , x ) , e 2 M . By arbitrariness of 𝑀 and t 0 , we deduce the first assertion as x (the case x is parallel).

For the second assertion, let h > 0 be an arbitrary height. Then we can prepare a smooth graphical curve γ lying below the initial curve γ 0 and agreeing with the horizontal line : = { p , e 2 = h } outside a compact set. Classical heat kernel estimates imply that the solution to the heat equation from γ uniformly converges to the line , which together with [52, Theorem 1.4] implies the same convergence of the curve shortening flow from γ . By the comparison principle, we have lim inf t inf x R γ ( t , x ) , e 2 h . The arbitrariness of h implies the assertion.

Remark 7.7

Remark 7.7 (Convergence proof using maximum principles)

Here, for comparison purposes, we sketch how maximum principles ensure convergence of the planar curve shortening flow to a line. Consider an initial planar embedded curve with ends C 1 -asymptotic to horizontal semi-lines with opposite directions. Then, by the comparison principle, the solution remains contained in a horizontal slab. Also, outside a vertical slab, the initial curve is graphical. Applying Angenent’s intersection number principle [7] to the solution and vertical lines implies that the graphical part remains graphical, so the non-graphical part stays in a compact set. Then an extension of Grayson’s argument [28] (cf. [57]) implies that, at some positive time, the solution becomes entirely graphical. Note also that the C 1 -asymptotic property of the ends is preserved along the flow, by simple comparisons with grim reapers and interior gradient estimates [22] (see also [13, Lemma 0.1]). Hence, by Ecker–Huisken’s convergence result [21, Theorem 5.1], the solution locally smoothly converges to a self-similar solution (in a rescaled sense). Therefore, this limit must be a horizontal line, which is the only graphical self-similar solution contained in the slab. In particular, identifying the limit in a suitable sense using maximum principles seems substantially hard, especially without the slab property or planarity.

7.3 Convergence for the surface diffusion and Chen’s flow

We now turn to fourth-order flows of curve shortening type including (SDF) and (CF), and obtain an analogous result to Theorem 2.2.

Theorem 7.8

Let μ 0 , ϑ R , and let γ : [ 0 , T ) × R R n be a maximal solution to

(7.15) t γ = s 2 κ + μ | κ | 2 κ + ϑ κ , s κ s γ ,

with initial datum γ 0 C ̇ ( R ; R n ) such that inf R | x γ 0 | > 0 . Suppose that D [ γ 0 ] < . Then the energy D [ γ ( t , ) ] continuously decreases in t 0 . In addition, the following dichotomy holds.

  1. If T < , then

    lim inf t T ( T t ) 1 / 4 γ ( t , ) | κ | 2 d s > 0 .

  2. If T = , then

    lim t γ ( t , ) | κ | 2 d s = 0 ,

    and moreover,

    lim t sup R | s γ ( t , ) e 1 | = 0 , lim t sup R | s m κ ( t , ) | = 0

    for all m N 0 . In particular, after reparametrization by arclength and translation, the solution locally smoothly converges to a horizontal line.

Proof

The decay property for 𝐷 is already shown in Proposition 4.4. Case i is a direct consequence of Theorem 7.1. We prove case ii, so suppose that T = . Once we establish the convergence

(7.16) lim t B [ γ ( t ) ] = 0 ,

then the remaining argument is completely parallel to the proof of Theorem 2.2ii, where we use Lemma 7.2 instead of Lemma 7.4, as well as Lemma 6.5 instead of Lemma 6.8.

We prove (7.16). By Proposition 4.4, we have D [ γ ( t ) ] D [ γ 0 ] for all t [ 0 , ) and | s κ | 2 d s L 1 ( 0 , ) . Let R > 1 and let η C c ( R ) be as in the proof of Proposition 4.3. Using s κ , s γ = | κ | 2 and integration by parts, we find

(7.17) s γ e 1 , s κ η d s = s γ e 1 , s κ + | κ | 2 s γ η d s = | κ | 2 s γ , e 1 η d s s γ e 1 , κ s η d s .

Moreover, since 𝜂 has compact support, we have

(7.18) | κ | 2 | s γ e 1 | 2 η d s | s γ e 1 | 2 d s | κ | 2 η 2 D [ γ 0 ] ( 2 | κ | | s κ | η d s + | κ | 2 s η d s ) 4 D [ γ 0 ] 2 | s κ | 2 η d s + | κ | 2 η d s + 2 D [ γ 0 ] | κ | 2 s η d s ,

Combining (7.17) and (7.18), and using Hölder, we find

(7.19) 2 | κ | 2 η d s = | κ | 2 ( | s γ e 1 | 2 + 2 s γ , e 1 ) η d s 4 D [ γ 0 ] 2 | s κ | 2 η d s + | κ | 2 η d s + 2 D [ γ 0 ] | κ | 2 s η d s + 2 ( 2 D [ γ 0 ] ) 1 / 2 ( | s κ | 2 η d s ) 1 / 2 + 2 ( 2 D [ γ 0 ] ) 1 / 2 ( | κ | 2 | s η | 2 d s ) 1 / 2 .

For almost all t 0 > 0 , we have κ ( t 0 ) L ( R ) and s κ ( t 0 ) L 2 ( d s ( t 0 ) ) . Using (4.5) and (4.6), after absorbing, it thus follows from sending R , η 1 in (7.19) that B [ γ ( t 0 ) ] < . Lemma 7.2 implies that B [ γ ( t ) ] < for all t t 0 , and hence for all t > 0 . Therefore, for t > 0 , sending R in (7.19) yields

(7.20) | κ | 2 4 D [ γ 0 ] 2 | s κ | 2 + 2 ( 2 D [ γ 0 ] ) 1 / 2 ( | s κ | 2 d s ) 1 / 2 .

As | s κ | 2 d s L 1 ( 0 , ) , we may argue as in (7.13) to find a sequence t j [ j , j + 1 ] with s κ L 2 ( d s ( t j ) ) 0 . It follows from (7.20) that also B [ γ ( t j ) ] 0 . By (7.7) (which is valid even if T = ) and Lemma 7.2, we may use Lemma B.1ii to get B [ γ ( t ) ] 0 as t . ∎

Analogously to Corollary 2.6, we obtain the following result by the same proof.

Corollary 7.9

Let γ 0 C ̇ ( R ; R 2 ) with inf R | x γ | > 0 and D [ γ 0 ] < . Suppose that γ 0 has nonzero rotation number. Then the unique solution to (7.15) with μ 0 starting from γ 0 blows up in finite time as in Theorem 7.8i.

7.4 Global existence and convergence for the 𝜆-elastic flow

Now we consider the 𝜆-elastic flow (𝜆-EF), which may be regarded as a gradient flow of E λ [ γ ] : = B [ γ ] + λ D [ γ ] . We first establish a general global existence theorem.

Theorem 7.10

Let λ 0 and γ 0 C ̇ ( R ; R n ) with inf R | x γ 0 | > 0 and E λ [ γ 0 ] < . If λ = 0 , suppose in addition that γ 0 is proper. Then there exists a unique, properly immersed, global-in-time solution γ : [ 0 , ) × R R n to (𝜆-EF) with initial datum γ 0 such that

γ γ 0 C ( [ 0 , T ] × R ) and inf [ 0 , T ] × R | x γ | > 0 for any T ( 0 , ) .

Proof

By Theorem 7.1, it suffices to show that the unique maximal solution

γ : [ 0 , T ) × R R n

from γ 0 satisfies sup t [ 0 , T ) B [ γ ( t ) ] < . This follows since, by Proposition 4.3, we have

B [ γ ( t ) ] E λ [ γ ( t ) ] E λ [ γ 0 ] < for all t [ 0 , T ) .

We now prove that any global solution to the 𝜆-elastic flow converges to a stationary solution, thus being an elastica. Recall that a curve is called an elastica if it solves the equation s 2 κ + 1 2 | κ | 2 κ λ κ = 0 for some λ R . We also call it a 𝜆-elastica to specify the value of 𝜆. A 𝜆-elastica corresponds to a stationary solution to the 𝜆-elastic flow. The λ = 0 case is also called a free elastica.

Theorem 7.11

Let λ 0 and let 𝛾 be the global-in-time solution to (𝜆-EF) obtained in Theorem 7.10. Then 𝛾 satisfies the following properties.

  1. For any sequence of times t j and of base points { s j } j R , there exist subsequences { t j } { t j } and { s j } { s j } such that the sequence of translated reparametrizations γ ̃ ( t j , + s j ) γ ̃ ( t j , s j ) locally smoothly converges as j to a 𝜆-elastica γ C ̇ ( R ; R n ) with γ ( 0 ) = 0 .

  2. The above limit satisfies the lower semi-continuity lim t E λ [ γ ̃ ( t ) ] E λ [ γ ] .

Proof

By Proposition 4.3, we have

(7.21) E λ [ γ ( t ) ] E λ [ γ 0 ] = 0 t | s 2 κ + 1 2 | κ | 2 κ λ κ | 2 d s d t .

In particular, sup t [ 0 , ) B [ γ ( t ) ] E λ [ γ 0 ] , and thus Lemmas 6.4 and 6.5 yield

(7.22) s m κ + s m κ C ( λ , m , E λ [ γ 0 ] ) ( 1 + t 2 m + 1 8 ) ,
(7.23) s m κ L 2 ( d s ) + s m κ L 2 ( d s ) C ( λ , m , E λ [ γ 0 ] ) ( 1 + t m 4 ) .

We now prove property i. Fix any sequences t j and { s j } R . By (7.22) and a standard compactness argument (as in the proof of Theorem 2.2ii), up to taking a further subsequence (without relabeling), the sequence of curves γ ̃ ( t j , + s j ) γ ̃ ( t j , s j ) locally smoothly converges to an arclength parametrized curve γ C ̇ ( R ; R n ) with γ ( 0 ) = 0 . To establish that γ is indeed a 𝜆-elastica, it suffices to show that

lim t V L 2 ( d s ) ( t ) = 0 , where V = s 2 κ + 1 2 | κ | 2 κ λ κ .

Indeed, by Fatou’s lemma, this would then imply V L 2 ( d s ) = 0 for γ , which confirms that γ is a 𝜆-elastica. In order to show that V L 2 ( d s ) ( t ) 0 , it suffices to prove that t V L 2 ( d s ) 2 ( t ) is globally Lipschitz on [ 1 , ) , since V L 2 ( d s ) L 2 ( 0 , ) by property (7.21). To that end, let t 1 . Using the definition of 𝑉 and also (6.2) (with m = 0 , 2 ), properties (7.22) and (7.23) above imply, for t 1 ,

(7.24) V + V L 2 ( d s ) + t V L 2 ( d s ) C ( λ , E λ [ γ 0 ] ) .

Taking any cutoff function η C c ( R ) , we find that

(7.25) d d t | V | 2 η d s = 2 V , t V η d s | V | 2 κ , V η d s .

Integrating (7.25), taking absolute values, sending η 1 while using (7.24) yield

| V L 2 ( d s ) 2 ( t ) V L 2 ( d s ) 2 ( t ) | C | t t |

for all t , t 1 , which is the desired global Lipschitz continuity.

Property ii then easily follows from (7.21), the invariance of the energy E λ under orientation-preserving reparametrizations, and Fatou’s lemma. ∎

We now examine the precise convergence behavior of the elastic flow with λ = 1 , and the free elastic flow ( λ = 0 ). We first address the free elastic flow, proving Theorem 2.7.

Proof of Theorem 2.7

Thanks to Theorems 7.10 and 7.11, we only need to discuss the asymptotic behavior. By Theorem 7.11ii, the free elastica γ in Theorem 7.11i satisfies B [ γ ] < . The well-known classification of elasticae (see e.g. [37, 45]) implies that any free elastica is either a straight line with B = 0 , or a periodic elastica with B = . Hence γ must be a line. Therefore, by Theorem 7.11i, | s m κ ( t j , s j ) | always converges to zero for any t j , { s j } R , and m N 0 . As in the proof of Theorem 2.2, the arbitrariness of t j , s j implies the assertion. ∎

From now on, we discuss the case λ > 0 . As mentioned in the introduction, up to rescaling, we may assume that λ = 1 , and only consider E = E 1 = B + D .

In this case, the stationary solution may not be a line even in the finite energy regime – it may be a so-called borderline elastica. We recall basic facts on the (planar) borderline elastica. A prototypical arclength parametrization is given by

(7.26) γ b ( s ) : = ( s 2 tanh s 2 sech s ) ;

cf. Figure 2. The tangential angle of γ b , satisfying s γ b = ( cos θ b , sin θ b ) , is given by

θ b ( s ) = 4 arctan ( e s ) ,

which increases from 0 to 2 π as 𝑠 varies from to ∞. The signed curvature is then given by k b ( s ) = s θ b ( s ) = 2 sech s = 4 e s ( 1 + e 2 s ) 1 . This satisfies k s s + 1 2 k 3 k = 0 . In particular, we can explicitly compute

(7.27) E [ γ b ] = R ( 1 2 k b 2 + 1 cos θ b ) d s = R ( 1 2 k b 2 + 8 tan 2 θ b 4 ( 1 + tan 2 θ b 4 ) 2 ) d s = R 16 e 2 s ( 1 + e 2 s ) 2 d s = [ 8 1 + e 2 s ] = 8 .

Theorem 7.12

Let λ = 1 and let 𝛾 be the global-in-time solution to (EF) obtained in Theorem 7.10. Then the sequence of curves and the limit curve in Theorem 7.11i further satisfy the following properties.

  1. The limit curve γ is given by either the straight line γ ( s ) = s e 1 , or the (planar) borderline elastica of the form

    γ ( s ) = ( s 2 tanh ( s + s 0 ) ) e 1 ( 2 sech ( s + s 0 ) ) ω + v 0

    for some s 0 R and ω span { e 2 , , e 2 } with | ω | = 1 , and

    v 0 : = ( 2 tanh s 0 ) e 1 ( 2 sech s 0 ) ω

    so that γ ( 0 ) = 0 .

  2. If s γ ̃ ( t j , s j ) e 1 as j , then γ is a straight line. Otherwise, γ is a borderline elastica.

  3. If γ is a straight line, then E [ γ ] = 0 , while if γ is a borderline elastica, then | κ ( s ) | = 2 sech ( s + s 0 ) and E [ γ ] = 8 .

Proof

Property i is shown by using the known classification of elasticae [37, 45] as follows. The only infinite-length elasticae with bounded bending energy are a straight line and a borderline elastica, because otherwise the squared curvature | κ | 2 is a nonzero periodic function. We also note that, since γ satisfies the elastica equation s 2 κ + 1 2 | κ | 2 κ κ = 0 , in the borderline case, the size of the loop (or equivalently, the curvature maximum) is uniquely determined; in particular, the curvature must be of the form | κ | = 2 sech ( s + s 0 ) , and thus the image γ is congruent to γ b . Finally, since E [ γ ] < follows by Theorem 7.11ii, both for γ being a line and a borderline elastica, the asymptotic tangent direction must be e 1 . This with γ ( 0 ) = 0 implies the desired form of γ .

Property ii follows from smooth convergence and the fact that the tangent direction of our borderline elastica is never rightward, i.e., s γ e 1 everywhere.

Finally, property iii follows by direct computation; cf. (7.27). ∎

Remark 7.13

In the elastic flow case λ > 0 , the choice of the base points { s j } in Theorem 7.11i is particularly important because the limit does depend on this choice. Indeed, since we always have s γ ̃ ( t , s ) e 1 as s ± (see Lemma 3.1), for any solution 𝛾, we can choose s t R large enough depending on t 0 so that γ ̃ ( t , + s t ) γ ̃ ( t , s t ) locally smoothly converges to the e 1 -axis as t . This even occurs when 𝛾 is a stationary solution given by a borderline elastica.

Remark 7.14

On the other hand, it is not clear whether the limit also depends on the choice of the time sequence t j even if we choose a constant sequence of base points. It is an interesting problem if the limit in Theorem 7.11i is always unique, i.e., independent of t j , provided that s j = 0 for all j N .

Next we give several conditions on initial data for full convergence as t such that we can identify the type of the limit elastica. We first complete the proof of global-in-space convergence results in Corollary 2.10.

Proof of Corollary 2.10

Theorem 7.11ii and Theorem 7.12i, iii imply that the only limit γ in Theorem 7.11i with E [ γ ] < 8 is the line γ ( s ) = s e 1 . Therefore, the arbitrariness of t j , s j yields the assertion. ∎

We also exhibit some conditions on initial data such that, for specific choices of base points, we obtain full convergence results, locally in space.

Corollary 7.15

Let λ = 1 and let 𝛾 be the global-in-time solution to (EF) obtained in Theorem 7.10. Let γ ̃ be the canonical reparametrization of 𝛾.

  1. If γ 0 has point symmetry such that γ 0 ( x ) = γ 0 ( x ) for x R , then the curve γ ̃ ( t , ) locally smoothly converges to a straight line as t .

  2. Suppose that γ 0 has reflection symmetry such that γ 0 ( x ) = γ 0 ( x ) 2 γ 0 ( x ) , e 1 e 1 for x R . If s γ 0 ( 0 ) = e 1 , then the curve γ ̃ ( t , ) γ ̃ ( t , 0 ) locally smoothly converges to a straight line as t . If otherwise s γ 0 ( 0 ) = e 1 , then for some { R t } t 0 O ( n ) , the curve R t [ γ ̃ ( t , ) γ ̃ ( t , 0 ) ] locally smoothly converges to a borderline elastica γ with s γ ( 0 ) = e 1 as t .

  3. Suppose that n = 2 (planar). If N [ γ 0 ] > 0 (resp. N [ γ 0 ] < 0 ), then there are sequences of base points { s t i } t 0 R , where i = 1 , , | N [ γ 0 ] | , such that, for each 𝑖, the curve γ ̃ ( t , + s t i ) γ ̃ ( t , s t i ) locally smoothly converges as t to a borderline elastica γ with s γ ( 0 ) = e 1 and N [ γ ] = 1 (resp. N [ γ ] = 1 ). If in addition | N [ γ 0 ] | 2 , then the base points can be chosen so that

    lim t ( s t i + 1 s t i ) = for each i = 1 , , | N [ γ 0 ] | 1 .

Proof

We first prove case i. By unique solvability, the point symmetry of the initial curve is preserved under the flow. Hence, in particular, γ ̃ ( t , 0 ) is fixed at the origin, and thus, by Theorem 7.11i (with s j 0 ), the curve γ ̃ ( t , ) sub-converges (without translation) to an elastica γ with point symmetry. No borderline elastica in Theorem 7.12i has such symmetry, so the limit is uniquely determined as the straight line γ ( s ) = s e 1 . Thanks to the uniqueness of the limit, the convergence is now upgraded to the full limit t .

We prove case ii. By symmetry, s γ 0 ( 0 ) { e 1 , e 1 } . As above, the reflection symmetry is preserved under the flow. This with smoothness of the flow implies that the initial condition on s γ 0 ( 0 ) is preserved under the flow, and hence s γ 0 ( 0 ) = s γ ( 0 ) { e 1 , e 1 } holds for the limit γ in Theorem 7.11i with s j 0 . If s γ ( 0 ) = e 1 , then by Theorem 7.12ii, the limit is again the unique straight line as above and thus the proof is complete. If s γ ( 0 ) = e 1 , then by Theorem 7.12ii, the limit is a borderline elastica. Since s γ is then injective, the point γ ( 0 ) must be the tip of the loop, so that the image of the limit curve is unique up to taking a suitable sequence of orthogonal matrices { R t } O ( n ) ; more precisely, in view of the parametrization in Theorem 7.12i, if n = 2 , then each matrix R t is either the identity or reflection across the e 1 -axis, while if n 3 , then R t is rotation around the e 1 -axis. Then again the uniqueness of the limit curve implies full convergence.

We finally prove case iii. We only discuss the case N : = N [ γ 0 ] > 0 ; the other is parallel. Since the rotation number N [ γ ( t , ) ] is preserved along the flow (Lemma B.2), there exists a smooth tangential angle θ : [ 0 , ) × R R , i.e., s γ ̃ ( t , s ) = ( cos θ ( t , s ) , sin θ ( t , s ) ) , such that

lim s θ ( t , s ) = 0 , lim s θ ( t , s ) = 2 π N .

Then, for each t 0 , there are 𝑁 points s t 1 < < s t N such that θ ( t , s t i ) = 2 π ( i 1 2 ) and s θ ( t , s t i ) 0 for i = 1 , , N . Since s γ ̃ ( t , s t i ) = e 1 and the signed curvature is nonnegative there, for each 𝑖, the curve γ ̃ ( t , + s t i ) γ ̃ ( t , s t i ) locally converges to the desired unique borderline elastica, namely γ = γ b γ b ( 0 ) (resp. γ is the vertical reflection of γ b γ b ( 0 ) ). If N 2 , then in view of the compatibility between those local convergences, the distance between any adjacent base points needs to diverge as t . ∎

Remark 7.16

Novaga–Okabe [53, Theorem 3.3] also claimed convergence to a borderline elastica assuming nonzero rotation number, but their argument seems not ruling out the possibility that the limit curve is a non-rightward straight line. Our argument not only improves their result, but also clearly rules out non-rightward lines using the finiteness of the direction energy.

The above results already reveal interesting dynamical aspects of the non-compact elastic flow. In particular, we have the following peculiar example.

Example 7.17

Example 7.17 (Escaping loops)

Consider the planar elastic flow 𝛾 starting from a planar initial curve γ 0 : R R 2 with E [ γ 0 ] < and reflection symmetry in the sense that

γ 0 ( x ) = γ 0 ( x ) 2 γ 0 ( x ) , e 1 e 1 for x R , and s γ 0 ( 0 ) = e 1 .

Then Corollary 7.15ii implies that, around the origin, the curve γ ̃ ( t , ) locally converges to a line (after translation).

In addition, note that, in this case, N [ γ 0 ] 2 Z . We further assume that N [ γ 0 ] = 2 . A typical situation is that two loops are positioned symmetrically, as in Figure 3. Then Corollary 7.15iii, combined with symmetry, implies the existence of a sequence { s t } t 0 ( 0 , ) such that

γ ̃ ( t , ± s t ) γ ̃ ( t , ± s t )

locally converges to a borderline elastica γ with N [ γ ] = 1 .

In view of the compatibility of the above local convergences, we need to have s t as t . This roughly describes the dynamics in which the two loops escape to infinity.

We also observe that the above solution directly shows that we cannot have global-in-space convergence. Indeed, now we have s γ ̃ ( t , ) e 1 = 2 since | s γ ̃ ( t , s t ) e 1 | = 2 , although s γ ̃ ( t , ) e 1 locally uniformly converges to zero. A similar argument works for any order of derivative.

Motivated by the above observations, we finally pose a conjecture on the general behavior of planar elastic flows. For a general planar curve γ : R R 2 with

inf R | x γ | > 0 and E [ γ ] < ,

we can obtain the general lower bound of the energy 𝐸 in terms of the absolute rotation number (in the spirit of the phase transition aspect [44]),

(7.28) E [ γ ] = R ( 1 2 s θ 2 + 1 cos θ ) d s | R 2 s θ 1 cos θ d s | = 2 | 0 2 π N [ γ ] 1 cos θ d θ | = 8 | N [ γ ] | .

Therefore, any planar elastic flow with E [ γ 0 ] < satisfies

(7.29) lim t E [ γ ( t ) ] 8 | N [ γ 0 ] | .

In fact, we conjecture that the energy of a general solution will be asymptotically “quantized” by the number of loops captured in Corollary 7.15iii in the sense that equality is attained in (7.29).

Conjecture 7.18

If n = 2 , then the elastic flow in Theorem 2.8 satisfies

lim t E [ γ ( t ) ] = 8 | N [ γ 0 ] | .

In view of (7.27), equality in (7.28) is attained for the borderline elastica γ b (with N = 1 ); hence the limit could be interpreted as N [ γ 0 ] copies of γ b .

8 Rigidity theorems for solitons

In this final section, we observe that simple applications of our convergence as well as energy decay results yield new rigidity theorems for solitons with finite direction energy.

Let γ : [ 0 , T ) × R R n be a solution to a flow with maximal existence time T ( 0 , ] and initial datum γ 0 . We call 𝛾 an expanding soliton (resp. shrinking soliton) if there is an increasing (resp. decreasing) function α : [ 0 , T ) ( 0 , ) with α ( 0 ) = 1 and lim t T α ( t ) = (resp. lim t T α ( t ) = 0 ) such that γ ( t ) is isometric to α ( t ) γ 0 for all t [ 0 , T ) . We also call 𝛾 a non-scaling soliton if T = and γ ( t ) is isometric to γ 0 for all t 0 ; typical examples are translators and rotators. In all cases, we call γ 0 the profile curve of the soliton 𝛾.

We first address the flows of curve shortening type.

Corollary 8.1

Let γ 0 C ̇ ( R ; R n ) be the profile curve of either a non-scaling or expanding soliton to one of (CSF), (SDF), (CF). Suppose that D [ γ 0 ] < . Then γ 0 is a straight line.

Proof

The expanding case is easily dealt with thanks to the energy decay results in Theorems 2.2 and 7.8, since an expanding dilation increases the direction energy. In the non-scaling case, the convergence results there imply that γ 0 must be a line. ∎

For (SDF), several non-stationary solitons are known in compact and non-compact cases, e.g., shrinking Bernoulli’s figure-eight [23] and a non-compact complete translator [55]. Also, in [34], an abstract existence theorem is obtained for expanding solitons of entire graphs with conical ends. See also [8, 33] for solitons with free boundary. For (CF), the only known non-stationary solitons would be the circles and the figure-eight [16]. It is not known whether there are shrinking solitons with finite direction energy for these flows.

As for the free elastic flow, we obtain a similar type of new rigidity result.

Corollary 8.2

Let γ 0 C ̇ ( R ; R n ) be the profile curve of either a non-scaling or shrinking soliton to (FEF). Suppose that B [ γ 0 ] < . Then γ 0 is a straight line.

Proof

The shrinking case can be treated by the energy decay in Theorem 2.7, while in the non-scaling case, the convergence result implies the assertion. ∎

Concerning expanding solitons to (FEF), the circles and the figure-eight are again explicit examples [49], and also the result in [34] implies existence of expanding entire graphs. It is unknown if the latter examples have finite bending energy.

Finally, we also obtain a new rigidity result for (EF), where we can rule out both expanding and shrinking solitons.

Corollary 8.3

Let γ 0 C ̇ ( R ; R n ) be the profile curve of either a non-scaling, expanding, or shrinking soliton to (EF). Suppose that E [ γ 0 ] < . Then γ 0 is stationary, i.e., either a straight line or a borderline elastica.

Proof

The expanding case (with D [ γ 0 ] > 0 ) is ruled out since if γ : [ 0 , ) × R R n is an expanding soliton, then the limit property α ( t ) implies D [ γ ( t ) ] , which contradicts E [ γ ( t ) ] E [ γ 0 ] < . A similar argument rules out the shrinking case (with B [ γ 0 ] > 0 ), where α ( t ) 0 implies B [ γ ( t ) ] . In the non-scaling case, again the convergence results imply the assertion. ∎

For (CF), (FEF), and (EF), non-stationary non-scaling solitons (with infinite energy) are not known.

Award Identifier / Grant number: JP21H00990

Award Identifier / Grant number: JP23H00085

Award Identifier / Grant number: JP24K00532

Funding source: Austrian Science Fund

Award Identifier / Grant number: 10.55776/ESP557

Funding statement: The first author is supported by JSPS KAKENHI Grant Numbers JP21H00990, JP23H00085, JP24K00532. The second author is funded in whole, or in part, by the Austrian Science Fund (FWF), grant number 10.55776/ESP557. Part of this work was done when the second author was visiting Kyoto University supported by a Mobility Fellowship of the Strategic Partnership Program between the University of Vienna and Kyoto University.

A Local well-posedness

In this section, we discuss local well-posedness for a class of arbitrary order geometric flows of complete curves. The analogous results in the case of compact curves are standard. The usual parabolic Hölder space approach is however doomed to fail in the non-compact situation. For the convenience of the reader, we provide some details here.

A similar result for arbitrary order flows on compact manifolds has been proven by Mantegazza–Martinazzi [41]. For m N 0 , T > 0 , and a suitable initial datum γ 0 : R R n , we consider the general initial value problem

(A.1) { t γ = ( 1 ) m s 2 m κ + F ( s γ , , s 2 m + 1 γ ) in ( 0 , T ) × R , γ ( 0 , ) = γ 0 , in R .

Here F : ( R n ) 2 m + 1 R n is smooth. We normalize (A.1) by assuming the coefficient of the leading order term to be ( 1 ) m . Of course, all the results in this section remain valid if the leading order term takes the form ( 1 ) m σ for some fixed σ > 0 by transforming time via t σ t . Examples include the flows discussed in this paper, i.e., (CSF), (SDF), (CF), and (𝜆-EF).

A.1 Transformation into a parabolic problem

Equation (A.1) is not parabolic due to its geometric invariance. In the spirit of the “De Turck trick”, this issue can be overcome by adding a suitable tangential velocity to the evolution law. We first make the following observation which is easily proven by induction on m N 0 .

Lemma A.1

For each m N 0 , there exists a polynomial 𝒫 such that

( 1 ) m s 2 m κ = [ ( 1 ) m x 2 m + 2 γ | x γ | 2 m + 2 + P ( | x γ | 1 , x γ , , x 2 m + 1 γ ) ] .

We now consider the operator of order 2 m + 2 given by

A [ γ ] : = ( 1 ) m x 2 m + 2 γ | x γ | 2 m + 2 + P ( | x γ | 1 , x γ , , x 2 m + 1 γ ) ,

so that ( 1 ) m s 2 m κ = A [ γ ] by Lemma A.1. Writing F [ γ ] : = F ( s γ , , s 2 m + 1 γ ) , we study the so-called analytic problem

(A.2) { t γ = A [ γ ] + F [ γ ] in ( 0 , T ) × R , γ ( 0 , ) = γ 0 in R .

We will see that sufficiently regular solutions to (A.2) can be reparametrized to give solutions of (A.1); see the proof of Theorem A.3 below. The leading order term on the right-hand side of (A.2) is now of the form

( 1 ) m | x γ | ( 2 m + 2 ) x 2 m + 2 ,

making (A.2) a quasilinear parabolic system of order 2 m + 2 .

We now formulate our short-time well-posedness result using the notation of parabolic Hölder spaces C l + α l , l + α ( [ 0 , T ] × R ) ; see [62]. We also denote by C ( [ 0 , T ] × R ; R n ) the space of u : [ 0 , T ] × R R n such that all derivatives t α x β u are bounded in [ 0 , T ] × R for any α , β N 0 .

Theorem A.2

Let α ( 0 , 1 ) , m N 0 , and l : = 2 m + 2 . Let γ 0 C ̇ l , α ( R ) such that inf R | x γ 0 | > 0 . Then there exist T > 0 and a unique γ : [ 0 , T ] × R R n with

(A.3) inf [ 0 , T ] × R | x γ | > 0

and

(A.4) γ γ 0 C l + α l , l + α ( [ 0 , T ] × R )

which solves (A.2). Moreover, there exist η > 0 and C > 0 such that, for any γ ̃ 0 with

γ 0 γ ̃ 0 C l , α ( R ) η ,

the corresponding solution γ ̃ of (A.2) exists on [ 0 , T ] and satisfies

(A.5) γ γ ̃ C l + α l , l + α ( [ 0 , T ] × R ) C γ 0 γ ̃ 0 C l , α ( R ) .

Further, if even γ 0 C ̇ ( R ) , then γ γ 0 C ( [ 0 , T ] × R ) .

This result can be established with classical tools from the theory of quasilinear parabolic systems. Since it is not the main focus of this work, we will use Theorem A.2 without proof, but instead, we give a rough sketch of the main idea. It can be proven by applying the contraction principle and parabolic maximal regularity in Hölder spaces, in particular [62, Theorem 4.10]. A subtlety is of course that the curves we consider are unbounded, thus naturally not in the parabolic Hölder spaces C l + α l , l + α ( [ 0 , T ] × R ) . Instead, as in [13], we consider the affine space given by (A.4) and use that the quantities involved in equation (A.1) do not depend on zeroth-order quantities of 𝛾, but only on its derivatives which are bounded in time and space in the regularity class (A.4). To obtain the continuous dependence (A.5), one may proceed as in [59, Theorem 5.1].

A.2 Local solutions to the geometric problem

By construction, Theorem A.2 gives a flow with a particular tangential velocity. We will now construct a solution to (A.1). It is worth pointing out that a subtlety of our construction based on an ODE argument is that it does not preserve a finite degree of smoothness; see [63, 64]. This is why we focus on the smooth category in the sequel.

Theorem A.3

Let γ 0 C ̇ ( R ) with inf R | x γ 0 | > 0 . Then there exist T > 0 and a solution γ : [ 0 , T ] × R R n of (A.1) with

(A.6) inf [ 0 , T ] × R | x γ | > 0

and

(A.7) γ γ 0 C ( [ 0 , T ] × R ) .

The solution is given by a time-dependent reparametrization of the solution constructed in Theorem A.2, and it is unique in the class of solutions satisfying (A.6) and (A.7). Moreover, if γ 0 is proper, then γ ( t , ) is proper for all t [ 0 , T ] .

Proof

Let 𝛾 be the solution of (A.2) constructed in Theorem A.2. Then, setting

ξ ̂ : = | x γ | 2 A [ γ ] , x γ ,

we have t γ = A [ γ ] + F [ γ ] + ξ ̂ x γ . From (A.3) and since γ γ 0 C ( [ 0 , T ] × R ) , we deduce ξ ̂ C ( [ 0 , T ] × R ) . For all x R , we now consider the ODE problem

(A.8) { t Φ t ( x ) = ξ ̂ ( t , Φ t ( x ) ) , t > 0 , Φ 0 ( x ) = x .

Local existence and uniqueness for fixed x R follows from classical ODE theory. Since ξ ̂ is uniformly bounded, the solution exists on [ 0 , T ] for all x R . Moreover, since

ξ ̂ C ( [ 0 , T ] × R ) ,

we have ( t , x ) Φ t ( x ) x C ( [ 0 , T ] × R ) . Differentiating (A.8) with respect to 𝑥 yields

(A.9) x Φ t ( x ) = exp ( 0 t x ξ ̂ ( t , Φ t ( x ) ) d t ) ,

so that Φ t : R R is an orientation-preserving diffeomorphism for all t [ 0 , T ] . In abuse of notation, in the following, we denote by X Φ t the composition of a time-dependent function 𝑋 with Φ t in the spatial variable, i.e., ( X Φ t ) ( t , x ) = X ( t , Φ t ( x ) ) . We set γ ̃ : = γ Φ t . Then γ ̃ ( 0 ) = γ 0 , γ ̃ γ 0 C ( [ 0 , T ] × R ; R n ) , and

t γ ̃ = t γ Φ t + ( x γ Φ t ) t Φ t = ( A [ γ ] + F [ γ ] ) Φ t + ξ ̂ Φ t ( x γ Φ t ) + ( x γ Φ t ) t Φ t = A [ γ ̃ ] + F [ γ ̃ ]

using (A.8) and the transformation of the geometric objects. We conclude that γ ̃ is a solution to (A.1) with (A.7). Also, (A.3) and (A.9) yield (A.6).

For the uniqueness part, let 𝛾 be any solution to (A.1) satisfying (A.6) and (A.7). By the definition of 𝒜, if Φ : R R satisfies inf R x Φ > 0 , then

A [ γ Φ ] = ( 1 ) m x 2 m + 2 ( γ Φ ) | x ( γ Φ ) | 2 m + 2 + P ( | x ( γ Φ ) | 1 , x ( γ Φ ) , , x 2 m + 1 ( γ Φ ) ) = ( 1 ) m ( x γ Φ ) x 2 m + 2 Φ | ( x γ Φ ) x Φ | 2 m + 2 + H [ γ , Φ ] .

As before, the composition with Φ is only in the spatial variable. By Faà di Bruno’s formula,

H [ γ , Φ ] = H ( | x γ Φ | 1 , | x Φ | 1 , x γ Φ , , x 2 m + 2 γ Φ , x Φ , , x 2 m + 1 Φ )

is a polynomial in its arguments. We now consider the PDE

(A.10) { t Φ = ( 1 ) m x 2 m + 2 Φ | x γ Φ | 2 m + 2 ( x Φ ) 2 m + 2 + H [ γ , Φ ] , x γ Φ | x γ Φ | 2 in ( 0 , T ) × R , Φ ( 0 , x ) = x for x R .

The structure of (A.10) is similar to (A.2), and thus the same techniques as for proving Theorem A.2 can be used to show that there exist T ̂ T and a unique solution Φ of (A.10) on [ 0 , T ̂ ] with

( t , x ) Φ ( t , x ) x C ( [ 0 , T ̂ ] × R )

satisfying inf [ 0 , T ̂ ] × R x Φ > 0 . Writing Φ t ( x ) : = Φ ( t , x ) , by construction,

t Φ t = A [ γ Φ t ] , x γ Φ t | x γ Φ t | 2 ,

so that, defining γ ̃ : = γ Φ t , we compute

t γ ̃ = t γ Φ t + ( x γ Φ t ) t Φ t = ( A [ γ ] γ + F [ γ ] ) Φ t + ( x γ Φ t ) t Φ t = A [ γ ̃ ] γ ̃ + F [ γ ̃ ] + A [ γ ̃ ] γ ̃ = A [ γ ̃ ] + F [ γ ̃ ] ,

where A [ γ ̃ ] γ ̃ is the tangential projection of A [ γ ̃ ] along γ ̃ . Thus γ ̃ solves (A.2) with

γ ̃ γ 0 C ( [ 0 , T ̂ ] × R )

and therefore coincides (on [ 0 , T ̂ ] ) with the unique solution from Theorem A.2 with initial datum γ 0 .

Hence, if γ 1 , γ 2 are two solutions to (A.1) on [ 0 , T ] with initial datum γ 0 satisfying (A.6) and (A.7), there exists 0 < T ̂ T and a smooth family of reparametrizations Φ t such that Φ 0 ( x ) = x and γ 1 = γ 2 Φ t . Since both satisfy (A.1), this implies

A [ γ 1 ] γ 1 + F [ γ 1 ] = t γ 1 = t γ 2 Φ t + ( x γ 2 Φ t ) t Φ t = ( A [ γ 2 ] γ 2 + F [ γ 2 ] ) Φ t + ( x γ 2 Φ t ) t Φ t = A [ γ 1 ] γ 1 + F [ γ 1 ] + ( x γ 2 Φ t ) t Φ t ,

where again we used the transformation of the geometric objects in the last step. It follows t Φ t = 0 , and thus γ 1 = γ 2 on [ 0 , T ̂ ] . Hence

T 0 : = sup { τ [ 0 , T ] γ 1 ( t ) = γ 2 ( t ) for t [ 0 , τ ] } > 0 .

If T 0 < T , one may repeat the same argument for the flow with initial datum γ 1 ( T 0 ) = γ 2 ( T 0 ) , contradicting maximality of T 0 . We conclude γ 1 = γ 2 on [ 0 , T ] .

Lastly, the properness of γ ( t , ) follows by observing that (A.7) and properness of γ 0 imply that | γ ( t , x ) | | γ 0 ( x ) | t t γ L ( [ 0 , T ] × R ) as | x | . ∎

Remark A.4

The finite degree of smoothness in Theorem A.2 is natural in view of the order of the PDE (A.2) and is enough for short-time existence of (A.1) if one allows for a tangential velocity in the flow. The class C ̇ in Theorem A.3 appears to prevent a loss of regularity in the ODE argument. However, it is not entirely clear if restricting to this class is necessary to have local well-posedness for (A.1).

Moreover, it is worth pointing out that, even if an arclength parametrized initial curve has derivatives of all orders, it may not be possible to have a short-time existence result. An instructive example for (CSF) is a locally convex curve with infinitely many loops such that the enclosed area converges to zero at infinity, which would however have unbounded curvature. This is in stark contrast to the compact case, where pointwise smoothness automatically yields C ̇ -regularity (or even C ).

B Technical lemmas

B.1 A general blow-up and convergence lemma

Lemma B.1

Let T ( 0 , ] and let f : [ 0 , T ) [ 0 , ) be a continuous function. Suppose that there are C > 0 and 1 p < q < such that, for all 0 t 1 t 2 < T ,

(B.1) f ( t 2 ) f ( t 1 ) C t 1 t 2 ( f ( t ) p + f ( t ) q ) d t .

  1. If T < and lim sup t T f ( t ) = , then

    (B.2) lim inf t T ( T t ) 1 q 1 f ( t ) > 0 .

  2. If T = and there exist t j with lim sup j ( t j + 1 t j ) < and f ( t j ) 0 , then

    lim t f ( t ) = 0 .

Proof

Using the inequality x p C ( p , q ) ( x + x q ) , x 0 , we may without loss of generality assume that (B.1) is satisfied with p = 1 .

To prove i, fix t 0 [ 0 , T ) and define F ( t ) : = f ( t 0 ) + C t 0 t ( f ( τ ) + f ( τ ) q ) d τ . Then 𝐹 is locally absolutely continuous on [ 0 , T ) and f ( t ) F ( t ) by (B.1). Hence

F ( t ) = C ( f ( t ) + f q ( t ) ) C ( F ( t ) + F ( t ) q ) .

For 𝑡 sufficiently close to 𝑇, we have F ( t ) > 0 . Indeed, otherwise, there exist t j T with F ( t j ) = 0 , implying f 0 on [ t 0 , t j ] and contradicting lim sup t T f ( t ) = . Since

1 F ( t ) + F ( t ) q = 1 F ( t ) F q 2 ( t ) 1 + F q 1 ( t ) ,

choosing t j T with f ( t j ) , we find that

C ( T t 0 ) lim j t 0 t j F ( t ) F ( t ) + F ( t ) q d t = lim j t 0 t j d d t ( log F ( t ) 1 q 1 log ( 1 + F ( t ) q 1 ) ) d t = 1 q 1 log ( 1 + F ( t 0 ) q 1 F ( t 0 ) q 1 ) = 1 q 1 log ( 1 + f ( t 0 ) q 1 f ( t 0 ) q 1 ) .

Taking t 0 T , we conclude that lim t T f ( t ) = . For 𝑡 close to 𝑇, we thus have f ( t ) 1 and (B.1) implies

f ( t 2 ) f ( t 1 ) 2 C t 1 t 2 f ( t ) q d t .

As above, we conclude that G ( t ) : = f ( t 0 ) + 2 C t 0 t f ( τ ) q d τ satisfies f ( t ) G ( t ) and

G ( t ) 2 C G ( t ) q .

Hence, taking t 0 close to 𝑇 and integrating from t 0 to 𝑇, we have

2 C ( T t 0 ) t 0 T d d t ( 1 1 q G ( t ) 1 q ) d t = 1 q 1 G ( t 0 ) 1 q = 1 q 1 f ( t 0 ) 1 q .

Renaming t 0 into 𝑡 and rearranging, (B.2) follows.

We prove part ii. By assumption, there are j 0 N and r > 0 such that t j + 1 t j r < for all j j 0 and such that

(B.3) sup j j 0 ( f ( t j ) + 1 j ) e 2 C r .

Fix any j j 0 . Let U ( t ) : = f ( t + t j ) , V ( t ) : = ( f ( t j ) + 1 j ) e 2 C t , and W : = U V . Then W ( 0 ) = 1 j < 0 and, since 𝑓 is continuous, the number

t : = sup { t [ 0 , r ] W ( τ ) 0 for all τ [ 0 , t ] }

is strictly positive, i.e., t ( 0 , r ] . For 0 t t r , we have 0 V ( t ) 1 by (B.3). It follows that

V ( t ) = V ( 0 ) + 2 C 0 t V ( τ ) d τ V ( 0 ) + C 0 t ( V ( τ ) + V q ( τ ) ) d τ .

On the other hand, by (B.1), we have

U ( t ) U ( 0 ) + C 0 t ( U ( τ ) + U q ( τ ) ) d τ .

Combining the above estimates gives W ( t ) < 0 for all t [ 0 , t ] , so by continuity, t = r . Hence we conclude f ( t ) ( f ( t j ) + 1 j ) e 2 C ( t t j ) for t [ t j , t j + r ] , which implies

sup t [ t j , t j + 1 ] f ( t ) sup t [ t j , t j + r ] f ( t ) e 2 C r ( f ( t j ) + 1 j ) .

This yields that f ( t ) 0 as t . ∎

B.2 Rotation number

We define the rotation number for planar complete curves and prove its preservation along smooth flows.

For a planar curve γ C ̇ 2 ( R ; R 2 ) such that

inf R | x γ | > 0 and D [ γ ] < ,

since s γ ( x ) e 1 as x ± by Lemma 3.1, we can define the rotation number by

N [ γ ] : = lim R 1 2 π R R k d s = lim R ( θ ( R ) θ ( R ) ) ,

where k : R R denotes the signed curvature and θ : R R the tangential angle of 𝛾. Note that, by the boundary condition at infinity, as in the closed curve case, N [ γ ] is always an integer.

Lemma B.2

Let γ : [ 0 , T ] × R R 2 be a family of smooth curves with

inf [ 0 , T ] × R | x γ | > 0

such that γ ( t ) C ̇ 2 ( R ) , lim t t γ ( t ) γ ( t ) C 1 ( R ) = 0 , and lim x ± s γ ( t , x ) = e 1 hold for each t [ 0 , T ] . Then, for all t [ 0 , T ] , we have N [ γ ( t , ) ] = N [ γ ( 0 , ) ] .

Proof

Since the range of 𝑁 is discrete, it is sufficient to show that, at each t 0 [ 0 , T ] , the map t N [ γ ( t , ) ] is continuous. Since s γ ( t 0 , x ) e 1 , there is R 0 > 0 such that, for all | x | R 0 , we have s γ ( t 0 , x ) , e 1 1 2 . By C 1 -continuity in time and the assumption on | x γ | , there is ε 0 > 0 such that if | t t 0 | ε 0 and | x | R 0 , then s γ ( t , x ) , e 1 1 3 . This implies that the curve γ ( t , ) is graphical outside [ R 0 , R 0 ] . Now we define γ ̂ ( t , ) by cutting off the second component of 𝛾, namely γ ̂ : = ( γ 1 , ζ γ 2 ) with ζ C c ( R ) such that χ [ R 0 , R 0 ] ζ χ [ 2 R 0 , 2 R 0 ] . In view of its topological nature, the rotation number is preserved by this procedure whenever | t t 0 | ε 0 , because it only replaces graphical parts with other graphical curves having the same tangent directions at x = ± R 0 and x ± . Therefore, if | t t 0 | ε 0 , then the curves γ ( t , ) and γ ̂ ( t , ) , and hence even γ ̂ ( t , ) | [ 2 R 0 , 2 R 0 ] , have the same rotation number 𝑁. Thus it follows that the problem is now reduced to showing that N [ γ ̂ ( t , ) | [ 2 R 0 , 2 R 0 ] ] is continuous at t = t 0 , but this is a standard matter since the tangents are always rightward at the endpoints, s γ ̂ ( t , ± 2 R 0 ) = e 1 , for t [ 0 , T ] ( t 0 ε 0 , t 0 + ε 0 ) . ∎

Acknowledgements

The authors would like to thank Ben Andrews for discussions about the curve shortening flow.

References

[1] D. J. Altschuler, S. J. Altschuler, S. B. Angenent and L. F. Wu, The zoo of solitons for curve shortening in R n , Nonlinearity 26 (2013), no. 5, 1189–1226. 10.1088/0951-7715/26/5/1189Search in Google Scholar

[2] S. J. Altschuler, Singularities of the curve shrinking flow for space curves, J. Differential Geom. 34 (1991), no. 2, 491–514. 10.4310/jdg/1214447218Search in Google Scholar

[3] S. J. Altschuler and M. A. Grayson, Shortening space curves and flow through singularities, J. Differential Geom. 35 (1992), no. 2, 283–298. 10.4310/jdg/1214448076Search in Google Scholar

[4] B. Andrews, B. Chow, C. Guenther and M. Langford, Extrinsic geometric flows, Grad. Stud. Math. 206, American Mathematical Society, Providence 2020. 10.1090/gsm/206Search in Google Scholar

[5] B. Andrews, J. McCoy, G. Wheeler and V.-M. Wheeler, Closed ideal planar curves, Geom. Topol. 24 (2020), no. 2, 1019–1049. 10.2140/gt.2020.24.1019Search in Google Scholar

[6] S. Angenent, Parabolic equations for curves on surfaces. I. Curves with 𝑝-integrable curvature, Ann. of Math. (2) 132 (1990), no. 3, 451–483. 10.2307/1971426Search in Google Scholar

[7] S. Angenent, Parabolic equations for curves on surfaces. II. Intersections, blow-up and generalized solutions, Ann. of Math. (2) 133 (1991), no. 1, 171–215. 10.2307/2944327Search in Google Scholar

[8] T. Asai and Y. Giga, On self-similar solutions to the surface diffusion flow equations with contact angle boundary conditions, Interfaces Free Bound. 16 (2014), no. 4, 539–573. 10.4171/ifb/329Search in Google Scholar

[9] Y. Bernard, G. Wheeler and V.-M. Wheeler, Concentration-compactness and finite-time singularities for Chen’s flow, J. Math. Sci. Univ. Tokyo 26 (2019), no. 1, 55–139. Search in Google Scholar

[10] S. Blatt, Loss of convexity and embeddedness for geometric evolution equations of higher order, J. Evol. Equ. 10 (2010), no. 1, 21–27. 10.1007/s00028-009-0038-2Search in Google Scholar

[11] B. Choi, K. Choi and P. Daskalopoulos, Convergence of curve shortening flow to translating soliton, Amer. J. Math. 143 (2021), no. 4, 1043–1077. 10.1353/ajm.2021.0027Search in Google Scholar

[12] K.-S. Chou, A blow-up criterion for the curve shortening flow by surface diffusion, Hokkaido Math. J. 32 (2003), no. 1, 1–19. 10.14492/hokmj/1350652421Search in Google Scholar

[13] K.-S. Chou and X.-P. Zhu, Shortening complete plane curves, J. Differential Geom. 50 (1998), no. 3, 471–504. 10.4310/jdg/1214424967Search in Google Scholar

[14] K.-S. Chou and X.-P. Zhu, The curve shortening problem, Chapman & Hall/CRC, Boca Raton 2001. 10.1201/9781420035704Search in Google Scholar

[15] P. Collet and J.-P. Eckmann, Space-time behaviour in problems of hydrodynamic type: A case study, Nonlinearity 5 (1992), no. 6, 1265–1302. 10.1088/0951-7715/5/6/004Search in Google Scholar

[16] M. K. Cooper, G. Wheeler and V.-M. Wheeler, Theory and numerics for Chen’s flow of curves, J. Differential Equations 362 (2023), 1–51. 10.1016/j.jde.2023.02.064Search in Google Scholar

[17] A. Dall’Acqua, C.-C. Lin and P. Pozzi, Elastic flow of networks: Long-time existence result, Geom. Flows 4 (2019), no. 1, 83–136. 10.1515/geofl-2019-0005Search in Google Scholar

[18] A. Dall’Acqua and P. Pozzi, A Willmore–Helfrich L 2 -flow of curves with natural boundary conditions, Comm. Anal. Geom. 22 (2014), no. 4, 617–669. 10.4310/CAG.2014.v22.n4.a2Search in Google Scholar

[19] A. Diana, Elastic flow of curves with partial free boundary, NoDEA Nonlinear Differential Equations Appl. 31 (2024), no. 5, Paper No. 96. 10.1007/s00030-024-00984-xSearch in Google Scholar

[20] G. Dziuk, E. Kuwert and R. Schätzle, Evolution of elastic curves in R n : Existence and computation, SIAM J. Math. Anal. 33 (2002), no. 5, 1228–1245. 10.1137/S0036141001383709Search in Google Scholar

[21] K. Ecker and G. Huisken, Mean curvature evolution of entire graphs, Ann. of Math. (2) 130 (1989), no. 3, 453–471. 10.2307/1971452Search in Google Scholar

[22] K. Ecker and G. Huisken, Interior estimates for hypersurfaces moving by mean curvature, Invent. Math. 105 (1991), no. 3, 547–569. 10.1007/BF01232278Search in Google Scholar

[23] M. Edwards, A. Gerhardt-Bourke, J. McCoy, G. Wheeler and V.-M. Wheeler, The shrinking figure eight and other solitons for the curve diffusion flow, J. Elasticity 119 (2015), no. 1–2, 191–211. 10.1007/s10659-014-9502-5Search in Google Scholar

[24] C. M. Elliott and H. Garcke, Existence results for diffusive surface motion laws, Adv. Math. Sci. Appl. 7 (1997), no. 1, 467–490. Search in Google Scholar

[25] C. M. Elliott and S. Maier-Paape, Losing a graph with surface diffusion, Hokkaido Math. J. 30 (2001), no. 2, 297–305. 10.14492/hokmj/1350911955Search in Google Scholar

[26] M. Gage and R. S. Hamilton, The heat equation shrinking convex plane curves, J. Differential Geom. 23 (1986), no. 1, 69–96. 10.4310/jdg/1214439902Search in Google Scholar

[27] H. Garcke, J. Menzel and A. Pluda, Long time existence of solutions to an elastic flow of networks, Comm. Partial Differential Equations 45 (2020), no. 10, 1253–1305. 10.1080/03605302.2020.1771364Search in Google Scholar

[28] M. A. Grayson, The heat equation shrinks embedded plane curves to round points, J. Differential Geom. 26 (1987), no. 2, 285–314. 10.4310/jdg/1214441371Search in Google Scholar

[29] H. P. Halldorsson, Self-similar solutions to the curve shortening flow, Trans. Amer. Math. Soc. 364 (2012), no. 10, 5285–5309. 10.1090/S0002-9947-2012-05632-7Search in Google Scholar

[30] J. Hättenschweiler, Curve shortening flow in higher dimension, Ph.D. Thesis, ETH Zürich, 2015. Search in Google Scholar

[31] G. Huisken, A distance comparison principle for evolving curves, Asian J. Math. 2 (1998), no. 1, 127–133. 10.4310/AJM.1998.v2.n1.a2Search in Google Scholar

[32] R. Jordan, D. Kinderlehrer and F. Otto, The variational formulation of the Fokker–Planck equation, SIAM J. Math. Anal. 29 (1998), no. 1, 1–17. 10.1137/S0036141096303359Search in Google Scholar

[33] T. Kagaya and Y. Kohsaka, Existence of non-convex traveling waves for surface diffusion of curves with constant contact angles, Arch. Ration. Mech. Anal. 235 (2020), no. 1, 471–516. 10.1007/s00205-019-01426-0Search in Google Scholar

[34] H. Koch and T. Lamm, Geometric flows with rough initial data, Asian J. Math. 16 (2012), no. 2, 209–235. 10.4310/AJM.2012.v16.n2.a3Search in Google Scholar

[35] E. Kuwert and R. Schätzle, The Willmore flow with small initial energy, J. Differential Geom. 57 (2001), no. 3, 409–441. 10.4310/jdg/1090348128Search in Google Scholar

[36] E. Kuwert and R. Schätzle, Gradient flow for the Willmore functional, Comm. Anal. Geom. 10 (2002), no. 2, 307–339. 10.4310/CAG.2002.v10.n2.a4Search in Google Scholar

[37] J. Langer and D. A. Singer, The total squared curvature of closed curves, J. Differential Geom. 20 (1984), no. 1, 1–22. 10.4310/jdg/1214438990Search in Google Scholar

[38] L. Langer, Dynamics of elastic wires: Preserving area without nonlocality, Calc. Var. Partial Differential Equations 64 (2025), no. 3, Paper No. 106. 10.1007/s00526-025-02948-0Search in Google Scholar

[39] C.-C. Lin, L 2 -flow of elastic curves with clamped boundary conditions, J. Differential Equations 252 (2012), no. 12, 6414–6428. 10.1016/j.jde.2012.03.010Search in Google Scholar

[40] C.-C. Lin, H. Schwetlick and D. T. Tran, An elastic flow for nonlinear spline interpolations in R n , Trans. Amer. Math. Soc. 375 (2022), no. 7, 4893–4942. 10.1090/tran/8639Search in Google Scholar

[41] C. Mantegazza and L. Martinazzi, A note on quasilinear parabolic equations on manifolds, Ann. Sc. Norm. Super. Pisa Cl. Sci. (5) 11 (2012), no. 4, 857–874. 10.2422/2036-2145.201010_001Search in Google Scholar

[42] C. Mantegazza, A. Pluda and M. Pozzetta, A survey of the elastic flow of curves and networks, Milan J. Math. 89 (2021), no. 1, 59–121. 10.1007/s00032-021-00327-wSearch in Google Scholar

[43] C. Mantegazza and M. Pozzetta, The Łojasiewicz–Simon inequality for the elastic flow, Calc. Var. Partial Differential Equations 60 (2021), no. 1, Paper No. 56. 10.1007/s00526-020-01916-0Search in Google Scholar

[44] T. Miura, Elastic curves and phase transitions, Math. Ann. 376 (2020), no. 3–4, 1629–1674. 10.1007/s00208-019-01821-8Search in Google Scholar

[45] T. Miura, Elastic curves and self-intersections, preprint (2024), https://arxiv.org/abs/2408.03020; to appear in MATRIX Ann. Search in Google Scholar

[46] T. Miura, M. Müller and F. Rupp, Optimal thresholds for preserving embeddedness of elastic flows, Amer. J. Math. 147 (2025), no. 1, 33–80. 10.1353/ajm.2025.a950273Search in Google Scholar

[47] T. Miura and S. Okabe, On the isoperimetric inequality and surface diffusion flow for multiply winding curves, Arch. Ration. Mech. Anal. 239 (2021), no. 2, 1111–1129. 10.1007/s00205-020-01591-7Search in Google Scholar

[48] T. Miura and G. Wheeler, Uniqueness and minimality of Euler’s elastica with monotone curvature, J. Eur. Math. Soc. (JEMS) (2025), 10.4171/JEMS/1708. 10.4171/JEMS/1708Search in Google Scholar

[49] T. Miura and G. Wheeler, The free elastic flow for closed planar curves, J. Funct. Anal. 289 (2025), no. 7, Article ID 111030. 10.1016/j.jfa.2025.111030Search in Google Scholar

[50] W. W. Mullins, Theory of thermal grooving, J. Appl. Phys. 28 (1957), 333–339. 10.1063/1.1722742Search in Google Scholar

[51] M. Nara and M. Taniguchi, Stability of a traveling wave in curvature flows for spatially non-decaying initial perturbations, Discrete Contin. Dyn. Syst. 14 (2006), no. 1, 203–220. 10.3934/dcds.2006.14.203Search in Google Scholar

[52] M. Nara and M. Taniguchi, The condition on the stability of stationary lines in a curvature flow in the whole plane, J. Differential Equations 237 (2007), no. 1, 61–76. 10.1016/j.jde.2007.02.012Search in Google Scholar

[53] M. Novaga and S. Okabe, Curve shortening-straightening flow for non-closed planar curves with infinite length, J. Differential Equations 256 (2014), no. 3, 1093–1132. 10.1016/j.jde.2013.10.009Search in Google Scholar

[54] M. Novaga and S. Okabe, Convergence to equilibrium of gradient flows defined on planar curves, J. reine angew. Math. 733 (2017), 87–119. 10.1515/crelle-2015-0001Search in Google Scholar

[55] W. J. Ogden and M. Warren, Grim raindrop: A translating solution to curve diffusion flow, preprint (2024), https://arxiv.org/abs/2412.15934. Search in Google Scholar

[56] S. Parkins and G. Wheeler, The polyharmonic heat flow of closed plane curves, J. Math. Anal. Appl. 439 (2016), no. 2, 608–633. 10.1016/j.jmaa.2016.02.033Search in Google Scholar

[57] A. Polden, Evolving curves, Honours Thesis, Australian National University, 1991. Search in Google Scholar

[58] A. Polden, Curves and Surfaces of least total curvature and fourth-order flows, Ph.d. thesis, Universität Tübingen, 1996. Search in Google Scholar

[59] J. Prüss and G. Simonett, Moving interfaces and quasilinear parabolic evolution equations, Monogr. Math. 105, Birkhäuser, Cham 2016. 10.1007/978-3-319-27698-4Search in Google Scholar

[60] P. Quittner and P. Souplet, Superlinear parabolic problems, Birkhäuser Adv. Texts Basler Lehrbücher, Birkhäuser/Springer, Cham 2019. 10.1007/978-3-030-18222-9Search in Google Scholar

[61] P. Rybka and G. Wheeler, A classification of solitons for the surface diffusion flow of entire graphs, Phys. D 477 (2025), Article ID 134702. 10.1016/j.physd.2025.134702Search in Google Scholar

[62] V. A. Solonnikov, On boundary value problems for linear parabolic systems of differential equations of general form, Trudy Mat. Inst. Steklov. 83 (1965), 3–163. Search in Google Scholar

[63] A. Spener, Short time existence for the elastic flow of clamped curves, Math. Nachr. 290 (2017), no. 13, 2052–2077. Search in Google Scholar

[64] A. Spener, Short time existence for the elastic flow of clamped curves, Math. Nachr. 291 (2018), no. 13, 2115–2116; erratum to Math. Nachr. 290 (2017), no. 13, 2052–2077. Search in Google Scholar

[65] X. Wang and W. Wo, On the stability of stationary line and grim reaper in planar curvature flow, Bull. Aust. Math. Soc. 83 (2011), no. 2, 177–188. 10.1017/S0004972710001942Search in Google Scholar

[66] X.-L. Wang and W.-F. Wo, On the asymptotic stability of stationary lines in the curve shortening problem, Pure Appl. Math. Q. 9 (2013), no. 3, 493–506. 10.4310/PAMQ.2013.v9.n3.a6Search in Google Scholar

[67] G. Wheeler, On the curve diffusion flow of closed plane curves, Ann. Mat. Pura Appl. (4) 192 (2013), no. 5, 931–950. 10.1007/s10231-012-0253-2Search in Google Scholar

[68] G. Wheeler and V.-M. Wheeler, Curve diffusion and straightening flows on parallel lines, Comm. Anal. Geom. 32 (2024), no. 7, 1979–2034. 10.4310/CAG.241204023527Search in Google Scholar

[69] Y. Y. Yang and X. X. Jiao, Curve shortening flow in arbitrary dimensional Euclidian space, Acta Math. Sin. (Engl. Ser.) 21 (2005), no. 4, 715–722. 10.1007/s10114-004-0426-zSearch in Google Scholar

Received: 2025-06-03
Revised: 2025-09-24
Published Online: 2025-11-05
Published in Print: 2026-01-01

© 2025 Walter de Gruyter GmbH, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 12.1.2026 from https://www.degruyterbrill.com/document/doi/10.1515/crelle-2025-0077/html
Scroll to top button