Integrable Flows for Starlike Curves in Centroaffine Space
We construct integrable hierarchies of flows for curves in centroaffine R³ through a natural pre-symplectic structure on the space of closed unparametrized starlike curves. We show that the induced evolution equations for the differential invariants are closely connected with the Boussinesq hierarch...
Збережено в:
Дата: | 2013 |
---|---|
Автори: | , , |
Формат: | Стаття |
Мова: | English |
Опубліковано: |
Інститут математики НАН України
2013
|
Назва видання: | Symmetry, Integrability and Geometry: Methods and Applications |
Онлайн доступ: | http://dspace.nbuv.gov.ua/handle/123456789/149224 |
Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Цитувати: | Integrable Flows for Starlike Curves in Centroaffine Space / A. Calini, T. Ivey, G.M. Beffa // Symmetry, Integrability and Geometry: Methods and Applications. — 2013. — Т. 9. — Бібліогр.: 30 назв. — англ. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-149224 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1492242019-02-20T01:26:45Z Integrable Flows for Starlike Curves in Centroaffine Space Calini, A. Ivey, T. Beffa, G.M. We construct integrable hierarchies of flows for curves in centroaffine R³ through a natural pre-symplectic structure on the space of closed unparametrized starlike curves. We show that the induced evolution equations for the differential invariants are closely connected with the Boussinesq hierarchy, and prove that the restricted hierarchy of flows on curves that project to conics in RP² induces the Kaup-Kuperschmidt hierarchy at the curvature level. 2013 Article Integrable Flows for Starlike Curves in Centroaffine Space / A. Calini, T. Ivey, G.M. Beffa // Symmetry, Integrability and Geometry: Methods and Applications. — 2013. — Т. 9. — Бібліогр.: 30 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 37K10; 53A20; 53C44 DOI: http://dx.doi.org/10.3842/SIGMA.2013.022 http://dspace.nbuv.gov.ua/handle/123456789/149224 en Symmetry, Integrability and Geometry: Methods and Applications Інститут математики НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
English |
description |
We construct integrable hierarchies of flows for curves in centroaffine R³ through a natural pre-symplectic structure on the space of closed unparametrized starlike curves. We show that the induced evolution equations for the differential invariants are closely connected with the Boussinesq hierarchy, and prove that the restricted hierarchy of flows on curves that project to conics in RP² induces the Kaup-Kuperschmidt hierarchy at the curvature level. |
format |
Article |
author |
Calini, A. Ivey, T. Beffa, G.M. |
spellingShingle |
Calini, A. Ivey, T. Beffa, G.M. Integrable Flows for Starlike Curves in Centroaffine Space Symmetry, Integrability and Geometry: Methods and Applications |
author_facet |
Calini, A. Ivey, T. Beffa, G.M. |
author_sort |
Calini, A. |
title |
Integrable Flows for Starlike Curves in Centroaffine Space |
title_short |
Integrable Flows for Starlike Curves in Centroaffine Space |
title_full |
Integrable Flows for Starlike Curves in Centroaffine Space |
title_fullStr |
Integrable Flows for Starlike Curves in Centroaffine Space |
title_full_unstemmed |
Integrable Flows for Starlike Curves in Centroaffine Space |
title_sort |
integrable flows for starlike curves in centroaffine space |
publisher |
Інститут математики НАН України |
publishDate |
2013 |
url |
http://dspace.nbuv.gov.ua/handle/123456789/149224 |
citation_txt |
Integrable Flows for Starlike Curves in Centroaffine Space / A. Calini, T. Ivey, G.M. Beffa // Symmetry, Integrability and Geometry: Methods and Applications. — 2013. — Т. 9. — Бібліогр.: 30 назв. — англ. |
series |
Symmetry, Integrability and Geometry: Methods and Applications |
work_keys_str_mv |
AT calinia integrableflowsforstarlikecurvesincentroaffinespace AT iveyt integrableflowsforstarlikecurvesincentroaffinespace AT beffagm integrableflowsforstarlikecurvesincentroaffinespace |
first_indexed |
2025-07-12T21:07:45Z |
last_indexed |
2025-07-12T21:07:45Z |
_version_ |
1837476842652565504 |
fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 9 (2013), 022, 21 pages
Integrable Flows for Starlike Curves
in Centroaffine Space?
Annalisa CALINI †‡, Thomas IVEY † and Gloria MARÍ BEFFA §
† College of Charleston, Charleston SC, USA
E-mail: calinia@cofc.edu, iveyt@cofc.edu
‡ National Science Foundation, Arlington VA, USA
§ University of Wisconsin, Madison WI, USA
E-mail: maribeff@math.wisc.edu
Received September 07, 2012, in final form February 27, 2013; Published online March 06, 2013
http://dx.doi.org/10.3842/SIGMA.2013.022
Abstract. We construct integrable hierarchies of flows for curves in centroaffine R3 through
a natural pre-symplectic structure on the space of closed unparametrized starlike curves. We
show that the induced evolution equations for the differential invariants are closely connected
with the Boussinesq hierarchy, and prove that the restricted hierarchy of flows on curves
that project to conics in RP2 induces the Kaup–Kuperschmidt hierarchy at the curvature
level.
Key words: integrable curve evolutions; centroaffine geometry; Boussinesq hierarchy; bi-
Hamiltonian systems
2010 Mathematics Subject Classification: 37K10; 53A20; 53C44
In honor of Peter Olver.
1 Introduction
1.1 Integrable evolutions of space curves
Much of the work on integrable curve evolution equations has been guided by the fundamental
role played by the differential invariants of the curve (e.g., curvature and torsion in the Euclidean
setting) in helping identify the curve evolution as an integrable one. Perhaps the most important
example in the case of space curves is that of the Localized Induction Equation (LIE)
γt = γx × γxx, (1.1)
describing the evolution of a curve with position vector γ(x, t) in R3, and Euclidean arclength
parameter x. The complete integrability of equation (1.1) was uncovered by the realization, due
to Hasimoto [9], that the function ψ = κ exp(i
∫
τ dx), of the curvature κ and torsion τ of γ, is
a solution of the cubic focusing nonlinear Schrödinger (NLS) equation
iψt + ψxx + 1
2 |ψ|
2ψ = 0,
one of the two best-known integrable nonlinear wave equations (the other being the KdV equa-
tion).
In this paper, we also use as a guiding principle the observation that many (but not all) inte-
grable curve evolutions have the property of local preservation of arclength, i.e., the associated
?This paper is a contribution to the Special Issue “Symmetries of Differential Equations: Frames, Invariants
and Applications”. The full collection is available at http://www.emis.de/journals/SIGMA/SDE2012.html
mailto:calinia@cofc.edu
mailto:iveyt@cofc.edu
mailto:maribeff@math.wisc.edu
http://dx.doi.org/10.3842/SIGMA.2013.022
http://www.emis.de/journals/SIGMA/SDE2012.html
2 A. Calini, T. Ivey and G. Maŕı Beffa
vector fields satisfy a non-stretching condition. For example, the LIE vector field W = γx×γxx
satisfies the condition δW ‖γx‖ = 0, where δW denotes the variation in the direction of W . Thus
the local arclength parameter x is independent of t, and the compatibility conditions γxt = γtx,
γxxt = γtxx, γxxxt = γtxxx (more commonly written as compatibility conditions of the Frenet
equations and the evolution equations for the Frenet frame) turn out to be equivalent to the
Lax pair of the NLS equation for ψ.
Indeed, many integrable curve evolutions in various geometries have been found by looking
for non-stretching vector fields that produce compatible equations for the moving frame of the
evolving curve; in the case of space curves, the geometries explored include Euclidean [11],
spherical [6], Minkowski [21], affine and centroaffine [4]. (Moreover, integrable curve evolutions
without preservation of arclength have been found in projective [17], conformal [16] and other
parabolic geometries.) The approach in these investigations involves finding suitable choices for
the coefficients of the non-stretching vector fields (relative to a Frenet-type frame) and often
assuming special relations among the differential invariants; thus it can be challenging to identify
integrable hierarchies.
Another approach to investigating the relation between a non-stretching curve evolution and
the integrable PDE system satisfied by the differential invariants is to seek a natural Hamil-
tonian setting for the curve flow. The LIE was shown by Marsden and Weinstein [18] to be
a Hamiltonian flow on a suitable phase space endowed with a symplectic form of hydrodynamic
origin (see also [1, 2]). In a fundamental paper [14] Langer and Perline used this framework to
explore in depth the correspondence between the LIE and NLS equations and, along the way, de-
rived a geometric recursion operator at the curve level that made it easy to obtain the integrable
hierarchies of both curve and curvature flows, as well as meaningful reductions thereof [12, 13].
In this article we study integrable evolution equations for closed curves in centroaffine R3
beginning, as in [14], with a natural pre-symplectic form on an appropriate infinite-dimensional
phase space. The Hamiltonian setting allows us to construct integrable hierarchies of curve flows
and the associated families of integrable evolution equations for the centroaffine differential
invariants (which turn out to be equivalent to the Boussinesq hierarchies). The motivation for
addressing the centroaffine case comes from an interesting article by Pinkall [26], who derived
a Hamiltonian evolution equation on the space of closed nondegenerate curves in the centroaffine
plane. The simple definition of the symplectic form in the planar case (related to the SL(2)-
invariant area form) suggests that an analogous description may be possible in the 3-dimensional
case, where a parallel could be drawn with the more familiar Euclidean case treated by [14].
Before describing the organization of the paper, we briefly discuss Pinkall’s original setting
and some results of ours for the planar case.
1.2 Pinkall’s flow in R2
Centroaffine differential geometry in Rn refers to the study of submanifolds and their proper-
ties that are invariant under the action of SL(n), not including translations1. For example,
a parametrized curve γ : I → Rn (where I is an interval on the real line) is nondegenerate if
det
(
γ(x), γ′(x), . . . , γ(n−1)(x)
)
6= 0
for all x ∈ I, and this property is clearly invariant under the action of SL(n). Thus, for these
curves the integral∫ ∣∣γ, γ′, . . . , γ(n−1)∣∣2/n(n−1) dx (1.2)
1Some authors [23] refer to this geometry as centro-equi-affine due to the choice of the unimodular group
SL(n), while using centro-affine to refer to geometry invariant under the general linear group GL(n).
Integrable Flows for Starlike Curves in Centroaffine Space 3
is SL(n)-invariant, and represents the centroaffine arclength, where for the sake of convenience
we use the notation
|X1, . . . , Xn| := det(X1, . . . , Xn)
for n-tuples of vectors Xi ∈ Rn. (The fractional power in (1.2) is necessary to make the integral
invariant under reparametrization.)
In the case where n = 2, Pinkall [26] defined a geometrically natural flow for nondegen-
erate curves in R2, which he referred to as star-shaped curves, as follows. Suppose that γ
is parametrized by centroaffine arclength s, so that |γ, γ′| = 1 identically. It follows that
γss = −p(s)γ where p(s) is defined as the centroaffine curvature. Along a closed curve γ, one
defines the skew-symmetric form
ω(X,Y ) =
∮
γ
|X,Y | ds, (1.3)
where X and Y are vector fields along γ. This pairing is nondegenerate on the space of vector
fields that locally preserve arclength. Then the symplectic dual with respect to (1.3) of the
functional
∮
γ p(s) ds is the vector field
X = 1
2ps
γ − pγs.
Pinkall’s flow γt = 1
2ps
γ − pγs induces an evolution equation for curvature that coincides
with the KdV equation, up to rescaling. In an earlier paper [3], we showed how to use solutions
of the (scalar) Lax pair for KdV to generate solutions of Pinkall’s flow. In particular, we showed
that varying the spectral parameter in the Lax pair for a fixed KdV potential q corresponds to
constructing a solution to the flow with curvature given by a Galileian KdV symmetry applied
to q. We also derived conditions under which periodic KdV solutions corresponded to smoothly
closed loops (for appropriate values of the spectral parameter) and illustrated this using finite-
gap KdV solutions.
1.3 Organization of the paper
In Section 2 we introduce basic notions concerning the differential geometry of nondegenerate
curves in centroaffine R3, including centroaffine arclength, differential invariants, and non-
stretching curve variations. This section also contains a discussion of the relation between
nondegenerate curves and parametrized maps into RP2. In Section 3 we generalize Pinkall’s
setting to R3 by introducing a pre-symplectic form on the space of closed unparametrized star-
like curves; we also compute Hamiltonian vector fields associated with the total length and
total curvature functionals. Flow by these vector fields induces evolution equations for the
differential invariants; we discuss these equations in Section 4, including their bi-Hamiltonian
formulation, Lax representation, and the connection with the Boussinesq equation. In Section 5
we show that the Poisson operators introduced in Section 4 give rise to the Boussinesq recursion
operator, generating a (double) hierarchy of commuting evolution equations for the differential
invariants. In Theorem 5.4, we relate the Hamiltonian structure for starlike curves and the
Poisson structure for the differential invariants, and obtain a double hierarchy of centroaffine
geometric evolution equations. We conclude Section 5, and the paper, by considering which of
these flows preserve the property that γ corresponds to a conic under the usual projectiviza-
tion map π : R3 → RP2. We show that the sub-hierarchy of conicity-preserving curve evolutions
induces the Kaup–Kuperschmidt hierarchy at the curvature level.
4 A. Calini, T. Ivey and G. Maŕı Beffa
2 Centroaffine curve flows in R3
2.1 Centroaffine invariants
Let γ : I → R3 be nondegenerate. We parameterize γ by centroaffine arclength, so that
|γ, γ′, γ′′| = 1. (2.1)
We assume for the rest of this subsection that x is an arclength parameter.
It follows by differentiating (2.1) with respect to x that
γ′′′ = p0γ + p1γ
′ (2.2)
for some functions p0(x) and p1(x). As explained below, these constitute a complete set of
differential invariants for nondegenerate curves.
Remark 2.1. Huang and Singer [10] refer to nondegenerate curves in centroaffine R3 as starlike.
They define invariants κ and τ which correspond to −p1 and p0 respectively. Labeling p0 as
torsion is appropriate, since nondegenerate curves that lie in a plane in R3 (not containing the
origin) are exactly those for which p0 is identically zero.
Remark 2.2. Some insight into the meaning of the centroaffine curve invariants can be gained
by considering the relationship between γ and the corresponding parametrized curve Υ = π ◦ γ
in RP2, where π : R3 → RP2 is projectivization. The nondegeneracy condition on γ corresponds
to Υ being regular and free of inflection points. Conversely, any such parametrized curve
Υ : R → RP2 has a unique lift to γ : R → R3 which is centroaffine arclength-parametrized; we
refer to γ as the canonical lift of Υ. When written in terms of Υ instead of γ, the invariants p0
and p1 are (up to sign) the well-known Wilczynski invariants [30]. Since these invariants define
a differential equation whose solution determines the curve uniquely up to the action of the
group SL(3), any other differential invariant must be functionally dependent on p0, p1 and their
x-derivatives.
According to Ovsienko and Tabachnikov [25], the cubic differential (p0 − 1
2p
′
1)(dx)3 has the
interesting property that it is invariant under reparametrizations of Υ. Curves in RP2 for which
this differential vanishes identically are conics. For curves for which the coefficient p0 − 1
2p
′
1 is
nowhere vanishing, one can define the projective arclength differential (p0 − 1
2p
′
1)
1/3dx. Those
parametrized curves in RP2 for which p0 − 1
2p
′
1 = C (a nonzero constant) are parametrized
proportional to projective arclength, and we use the same terminology for their canonical lifts
into R3. (Note that, in this case, the projective arclength differential is C1/3 times the cen-
troaffine arclength differential dx.)
Along a nondegenerate curve, an analogue of the Frenet frame is provided by vectors γ, γ′, γ′′.
In fact, if we combine them as columns in an SL(3)-valued matrix W = (γ, γ′, γ′′), then the
analogue of the Frenet equations is
Wx = W
0 0 p0
1 0 p1
0 1 0
.
However, for later use it will be convenient to define a different SL(3)-valued frame F (x) =
(γ, γ′, γ′′ − p1γ) which satisfies the Frenet-type equation
Fx = FK, K =
0 k1 k2
1 0 0
0 1 0
, (2.3)
where k1 = p1, k2 = p0 − p′1. Of course, k1, k2 also constitute a complete set of differential
invariants, and we will come to use these in place of the Wilczynski invariants from Section 4
onwards.
Integrable Flows for Starlike Curves in Centroaffine Space 5
2.2 Non-stretching variations
Suppose that Γ : I × (−ε, ε) → R3 is a smooth mapping such that, for fixed t, Γ(x, t) is
a nondegenerate curve parametrized by x. Without loss of generality, we will assume that
γ(x) = Γ(x, 0) is parametrized by centroaffine arclength. Let X denote the variation of γ in the
t-direction, and expand
X =
∂
∂t
∣∣∣∣
t=0
Γ = aγ + bγ′ + cγ′′.
(We still use primes to denote derivatives with respect to x, although x is not necessarily an
arclength parameter along the curves in the family for t 6= 0.)
To compute the variation of the arclength differential |γ, γ′, γ′′|1/3 dx, we introduce the no-
tation δ for variation in the t-direction along γ. Using the relation (2.2), we compute
δγ′ = X ′ = (a′ + p0c)γ + (a+ b′ + p1c)γ
′ + (b+ c′)γ′′,
δγ′′ = X ′′ = (a′′ + 2p0c
′ + p′0c+ p0b)γ + (2a′ + p0c+ b′′ + 2p1c
′ + p′1c+ p1b)γ
′
+ (a+ 2b′ + p1c+ c′′)γ′′.
Then
δ|γ, γ′, γ′′| = |δγ, γ′, γ′′|+ |γ, δγ′, γ′′|+ |γ, γ′, δγ′′| = 3a+ 3b′ + c′′ + 2p1c.
In particular, the variation X preserves the centroaffine arclength differential if and only if
b′ = −a− 1
3(c′′ + 2p1c), (2.4)
i.e.,
X = aγ −
(∫ (
a+ 1
3c
′′ + 2
3p1c
)
dx
)
γ′ + cγ′′. (2.5)
We refer to vector fields of this form as non-stretching, since not only do such variations preserve
the overall arclength of, say, a closed loop, but also no small portion of the curve is stretched or
compressed.
3 Hamiltonian curve flows
3.1 Symplectic structure on starlike loops
Generalizing Pinkall’s setting [26] for planar star-shaped loops to the three-dimensional case, we
introduce the infinite-dimensional space
M̂ = {γ : S1 → R3 : |γ, γ′, γ′′| = 1},
as a subset of the vector space V = Map(S1,R3) of C∞ maps from S1 to R3. Assume that
γ ∈ M̂ , i.e., γ is a closed starlike curve parametrized by centroaffine arclength; then a vector
field X = aγ+ bγ′+ cγ′′ is in the tangent space TγM̂ if and only if X is of the form (2.5), where
the coefficients a and c are 2π-periodic functions of x and satisfy the “zero mean” condition∮
γ
(
a+ 2
3p1c
)
dx = 0,
ensuring that the coefficient of γ′ in (2.5) is also periodic.
6 A. Calini, T. Ivey and G. Maŕı Beffa
On V , define the skew-symmetric form
ωγ(X,Y ) =
∮
γ
|X, γ′, Y |dx, X, Y ∈ TγV. (3.1)
Note that ω is automatically closed (that is, dω = 0) since the integrand in (3.1) is a volume
form on R3 [1, 2].
Letting X = aγ + bγ′ + cγ′′, Y = ãγ + b̃γ′ + c̃γ′′, we compute
ωγ(X,Y ) =
∮
γ
(ac̃− ãc) dx. (3.2)
Assuming that γ ∈ M̂ and X,Y ∈ TγM̂ , then ωγ(X,Y ) = 0 for all Y if and only if a = −2
3p1c
and c is constant. Thus, the restriction of ω to M̂ is a degenerate closed 2-form (a pre-symplectic
form), with kernel given by the subspace RZ0 +RZ1 of constant-coefficient linear combinations
of the vector fields Z0 = γ′ and Z1 = γ′′− 2
3p1
γ (corresponding to c = 0 and c = 1 respectively).
Note that this degeneracy is the result of restricting the 2-form (3.1) to the space of closed
curves satisfying the arclength constraint. Degeneracy coming from constraints is common
when defining symplectic structures on loop spaces [19] and phase spaces of nonlinear evolution
equations [7].
Remark 3.1. The generators Z0 and Z1 of the kernel of ω will turn out to be the seeds of the
double hierarchy of curve flows discussed in Section 5. A similar situation is encountered for
the LIE hierarchy [14], where the seed γ′ spans the kernel of the natural pre-symplectic form on
loops in Euclidean R3.
One could attempt to remove the degeneracy by constructing a quotient of M̂ with respect
to group actions generated by flowing by Z0 and Z1. Flow by Z0 generates an action of the
additive group R, simply by translation in x. The resulting quotient space M = M̂/R can be
identified with the space of unparametrized starlike loops. Because of translation invariance, ω
descends to a give a well-defined closed 2-form on M , with one-dimensional kernel RZ1. On the
other hand, we do not know of a natural geometric interpretation for the quotient of M by flow
under Z1, on which ω would become non-degenerate.
However, ω can still be used to define a link between vector fields and functionals, and we
will see that Z1 is linked in this way to the arclength functional.
3.2 Examples
Recall that the correspondence between vector fields XH and (differentials of) Hamiltonians
H ∈ C∞(M) on a manifold M with symplectic form ω is defined by the relation
dH[X] = ωγ(X,XH), ∀X ∈ TγM, (3.3)
XH being the Hamiltonian vector field corresponding to H. However, when ω is degenerate
the correspondence is no longer an isomorphism: for those functionals H for which there is
a Hamiltonian vector field, XH is only defined up to addition of elements in the kernel of ω.
We will use (3.3) to compute Hamiltonian vector fields for a few interesting functionals; to do
so, we will initially work in the ambient space V , and use the arclength-preserving condition (2.4)
to rewrite the differential in a form suitable for applying (3.3).
We first consider the arclength functional
L(γ) =
∮
γ
∣∣γ, γ′, γ′′∣∣1/3dx. (3.4)
Integrable Flows for Starlike Curves in Centroaffine Space 7
on the space V . Given an arbitrary vector field X = aγ + bγ′ + cγ′′ (not necessarily arclength
preserving), the variation of the determinant in (3.4) along X is given by
δ|γ, γ′, γ′′| = |X, γ′, γ′′|+ |γ,X ′, γ′′|+ |γ, γ′, X ′′| = (3a+ 3b′ + c′′ + 2p1c)|γ, γ′, γ′′|.
Assume now that γ ∈ M̂ , so that |γ, γ′, γ′′| = 1 and
dL[X] =
∮
γ
(
a+ 2
3p1c
)
dx.
We now seek a vector field XL = ãγ + b̃γ′ + c̃γ′′ ∈ TγM̂ such that dL[X] = ωγ(X,XL) =∮
γ(ac̃ − ãc)dx. Using the non-stretching condition (2.4), we obtain the following Hamiltonian
vector field
XL ≡ Z1 = γ′′ − 2
3p1
γ
(which is unique only up to adding a constant times Z0). In Section 4.4 we will see that the
associated curve flow γt = XL leads to the Boussinesq equation for the curvatures k1, k2.
Next, we introduce the total curvature functional
P (γ) =
∮
γ
p1 dx, γ ∈ M̂. (3.5)
From γ′′′ = p0γ + p1γ
′ and (2.1), it follows that p1 = |γ, γ′′′, γ′′|. Then the variation of p1 along
an arbitrary vector field X = aγ + bγ′ + cγ′′ is given by
δp1 = |X, γ′′′, γ′′|+ |γ,X ′′′, γ′′|+ |γ, γ′′′, X ′′|
= |X, p0γ + p1γ
′, γ′′|+ |γ,X ′′′, γ′′|+ |γ, p0γ + p1γ
′, X ′′|
= 3p1a+ 3c′p0 + 2cp′0 + 3a′′ + b′′′ + 4c′′p1 + 3c′p′1 + cp′′1 + 5b′p1 + bp′1 + 2cp21. (3.6)
Then, up to perfect derivatives,
dP [X] =
∮
γ
δp1 dx =
∮
γ
3p1a+ 3c′p0 + 2cp′0 + 4c′′p1 + 3c′p′1 + cp′′1 + 4b′p1 + 2cp21 dx.
Assuming X is an arclength-preserving vector field, we set b′ = −a− 1
3(c′′+2p1c) and compute∮
γ
δp1 dx =
∮
γ
−ap1 + 3c′p0 + 2cp′0 + 8
3c
′′p1 + 3c′p′1 + cp′′1 − 2
3cp
2
1 dx.
Integrating by parts, we arrive at
dP (X) =
∮
γ
−p1a+
(
−p′0 + 2
3p
′′
1 − 2
3p
2
1
)
cdx. (3.7)
Suppose that XP = ãγ + b̃γ′ + c̃γ′′ is also arclength-preserving. Setting the right-hand side
of (3.7) equal to ωγ(X,XP ) =
∮
γ
(ac̃ − ãc) dx, we get ã = p′0 − 2
3p
′′
1 + 2
3p
2
1 and c̃ = −p1. Using
equation (2.4) we compute b̃′ = −
(
p′0 − 2
3p
′′
1 + 2
3p
2
1
)
− 1
3
(
−p′′1 − 2p21
)
=
(
p′1 − p0
)′
, a perfect
derivative. Thus, a Hamiltonian vector field corresponding to (3.5) is
XP =
(
2
3
(
p21 − p′′1
)
+ p′0
)
γ + (p′1 − p0)γ′ − p1γ′′.
(This is only unique up to adding a linear combination of Z0 and Z1.) Again, we will see that
the associated curve flow γt = XP is also directly related, at the level of the curvatures, to one
of the flows in the Boussinesq hierarchy.
8 A. Calini, T. Ivey and G. Maŕı Beffa
4 Integrable centroaffine curve flows
In this section we will examine the evolution of centroaffine curvatures induced by the curve
flows defined in Section 3.2. We begin by computing the evolution of invariants under more
general curve flows.
4.1 Evolution of invariants
First, we consider how the centroaffine invariants of a starlike curve evolve under a general
non-stretching evolution equation
γt = r0γ + r1γ
′ + r2γ
′′ (4.1)
with r0 = −r′1 − 1
3(r′′2 + 2p1r2). (Thus, we assume from now on that γ(x, t) is parametrized
by arclength x at each time.) Of course, in order for (4.1) to represent a geometric evolution
equation, r1 and r2 should be functions of the invariants p0, p1 and their arclength derivatives.
Proposition 4.1. The evolution equations induced by (4.1) for the Wilzcynski invariants are
(p0)t = −r′′′′1 + p1r
′′
1 + 3p0r
′
1 + p′0r1
− 1
3
(
r′′′′′2 + p1r
′′′
2
)
+
(
(3p0 − 2p′1)r
′
2
)′
+ 2
3
(
p21r
′
2 + (p1p
′
1 − p′′′1 )r2
)
+ p′′0r2, (4.2)
(p1)t = −2r′′′1 + 2p1r
′
1 + p′1r1 − r′′′′2 + p1r
′′
2 + (3p0 − p′1)r′2 + (2p′0 − p′′1)r2. (4.3)
Proof. The second equation (4.3) follows by substituting a = −r′1 − 1
3(r′′2 + 2p1r2), c = r2 in
the last line of (3.6). Similarly, using p0 = |γ′′′, γ′, γ′′|, we obtain (4.2) by computing
(p0)t =
∣∣(γt)′′′, γ′, γ′′∣∣+
∣∣p0γ + p1γ
′, (γt)
′, γ′′
∣∣+ p0
∣∣γ, γ′, (γt)′′∣∣. �
We note that these evolution equations previously appeared in [4].
From now on, we will take k1 = p1 and k2 = p0 − p′1 as fundamental invariants; one reason
for doing this is that the evolution equations for these invariants induced by (4.1) take the form(
k1
k2
)
t
= P
(
r1
r2
)
, (4.4)
where P is the skew-adjoint matrix differential operator
P =
(
−2D3 +Dk1 + k1D −D4 +D2k1 + 2Dk2 + k2D
D4 − k1D2 + 2k2D +Dk2
2
3
(
D5 + k1Dk1 − k1D3 −D3k1
)
+
[
k2, D
2
]) , (4.5)
D stands for the derivative with respect to x and [·, ·] denotes the commutator on pairs of
operators2. This operator P, which arises naturally when using k1, k2 instead of p0, p1, will
play a significant role in the integrable structure of the flows we study.
4.2 Two integrable flows
The vector field XL induces a non-stretching evolution equation
γt = γ′′ − 2
3k1
γ. (4.6)
2Note that expressions like Dk1 and Dk2 denote composition of D with multiplication by k1 and k2, respectively.
The skew-adjointness of P is easy to check, given that D is skew-adjoint.
Integrable Flows for Starlike Curves in Centroaffine Space 9
(This will be the first non-trivial curve evolution in the hierarchy discussed in Section 5, where
the right-hand side is labeled as Z1.) By setting r1 = 0, r2 = 1 in (4.4), we obtain the
corresponding curvature evolution(
k1
k2
)
t
= P
(
0
1
)
=
(
k′′1 + 2k′2
2
3(k1k
′
1 − k′′′1 )− k′′2
)
. (4.7)
This PDE system for curvatures is Hamiltonian, since it can be written in the form(
k1
k2
)
t
= PEk2,
where E denotes the vector-valued Euler operator
Ef =
∑
j≥0
(−D)j
∂f
∂k
(j)
1
,
∑
j≥0
(−D)j
∂f
∂k
(j)
2
T
(4.8)
on scalar functions f of k1, k2 and their higher x-derivatives k
(j)
1 , k
(j)
2 . (One can check that
the Poisson bracket defined using the Hamiltonian operator P on the appropriate function
space – see Section 5.2 below – satisfies the usual requirements of skew-symmetry and the
Jacobi identity.)
Moreover, (4.7) can also be written in Hamiltonian form as(
k1
k2
)
t
= QEρ3, (4.9)
for a different Hamiltonian operator and density
Q =
(
0 D
D 0
)
, ρ3 := 1
3(k′1)
2 + k2k
′
1 + k22 + 1
9k
3
1. (4.10)
(The notation ρ3 is explained below.) Since the curvature evolution can be written in Hamil-
tonian form in two ways (4.7) and (4.9), the integrals
∫
k2 dx and
∫
ρ3 dx are conserved by the
flow (for appropriate boundary conditions).
Remark 4.2. In fact, the curvature evolution here is a bi-Hamiltonian system, because P and Q
are a Hamiltonian pair, i.e., their linear combinations form a pencil of Hamiltonian operators,
and a pencil of compatible Poisson structures. This assertion can be verified mechanically (see,
e.g., Section 7.1 in [22] for details), but it also follows from the fact that, at least in the periodic
case, the Poisson structures are reductions of a well-known compatible pencil of Poisson brackets
on the space of loops in sl(3). (Indeed, when γ is periodic –or more generally has monodromy –
the matrix K in (2.3) provides a lift into this loop space.) The proof of the reduction of these
brackets can be found in [8], where P is linked to the Adler–Gel’fand–Dikii bracket for sl(3)
and Q is associated to its companion. The brackets were later linked to curve evolutions and
differential invariants in [15], where more details are available.
The (negative of the) Hamiltonian vector field XP of Section 3.2 induces the non-stretching
evolution
γt = k1γ
′′ + k2γ
′ + r0γ, (4.11)
where
r0 = −
(
k′2 + 1
3
(
k′′1 + 2k21
))
. (4.12)
10 A. Calini, T. Ivey and G. Maŕı Beffa
(The right-hand side of (4.11) is labeled as Z2 in the hierarchy discussed in Section 5.) We
similarly obtain the curvature evolution equations induced by this flow by setting r1 = k2,
r2 = k1 in (4.4). We remark that the resulting system is also bi-Hamiltonian, since it can be
written as(
k1
k2
)
t
= P Eρ2 = QEρ4 (4.13)
for
ρ2 = k1k2, ρ4 = 1
3(k′′1)2 + k′′1
(
k′2 − k21
)
− k1(k′1)2 + (k′2)
2 − k21k′2 + 1
9k
4
1 + 2k1k
2
2.
Thus,
∫
ρ2 dx and
∫
ρ4 dx are conserved integrals for (4.11). (Because XP corresponds symplec-
tically to the Hamiltonian
∫
k1 dx, it is automatic that this integral is also conserved.)
Remark 4.3. The arclength normalization (2.1) is preserved by the simultaneous rescaling
x 7→ λx, γ 7→ λ−1γ. Under this rescaling, k1 and k2 scale by multiples λ2 and λ3 respectively.
Thus, we may assign scaling weights 2 and 3 respectively to these curvatures, and each x-
derivative taken increases weight by one.
It will turn out (see Section 5.1 below) that the conserved densities for evolution equa-
tions (4.6) and (4.11) are all of homogeneous weight, with one density for each positive weight
not congruent to 1 modulo 3. We will number the densities in order of increasing weight, letting
ρ0 = k1, ρ1 = k2 and so on; thus, the density in (4.9) is denoted by ρ3, since its weight falls
between those of ρ2 and ρ4.
The curve flows (4.6) and (4.11) turn out to share the same conservation laws; for example,∫
k1 dx is conserved by (4.6) because (4.7) implies that
(k1)t = D(k′1 + 2k2).
Similarly, (4.11) conserves
∫
k2 dx because (4.13) implies that
(k2)t = D
(
2
3k
(4)
1 + k′′′2 − 2k1k
′′
1 − (k′1)
2 − 2k1k
′
2 + 4
9k
3
1 + 2k22
)
.
In Section 5 we will show that these flows share an infinite sequence of conservation laws.
4.3 Lax representation
In this subsection we use geometric considerations to derive Lax pairs for curvature evolution
equations induced by (4.6) and (4.11).
In [3] we found that the components of the solution γ(x, t) of Pinkall’s flow satisfied the
scalar Lax pair for the KdV equation. In the same spirit, we seek a system of the form
Ly = 0, yt =My, (4.14)
satisfied by each component of γ, where L andM are differential operators in x with coefficients
involving k1, k2. Using (2.2), we see that every component of γ satisfies the scalar ODE y′′′ =
(k1y)′ + k2y, and so we will let
L := D3 −Dk1 − k2
and seek operators M1 for (4.6) and M2 for (4.11).
In the case of (4.6), the components of γ also satisfy yt = y′′ − 2
3k1y, so we choose
M1 := D2 − 2
3k1.
Integrable Flows for Starlike Curves in Centroaffine Space 11
One can then verify that (4.7) implies that
Lt = [M1,L]. (4.15)
In the case of (4.11), the components of γ satisfy yt = k1y
′′ + k2y
′ + r0y, with r0 as given
by (4.12). So, we might set M2 = k1D
2 + k2D + r0. However, (4.14) would also be satisfied
if we modify M2 by adding NL, where N is an arbitrary differential operator. In fact, the
system (4.13) actually implies that Lt = [M2,L] for
M2 :=
(
k1D
2 + k2D + r0
)
− 3DL.
Writing these systems in Lax form (4.14) enables us to interpolate a spectral parameter into the
linear equations satisfied by the components. Thus, consider solutions of the compatible system
Ly = λy, yt =Mjy, (4.16)
where j = 1 or j = 2. Of course, the components of the evolving curve satisfy (4.16) only when
λ = 0. When λ 6= 0, we can construct solutions of the curve flow using solutions of (4.16):
Proposition 4.4. Let k1, k2 satisfy the evolution equation (4.7) for j = 1 or (4.13) for j = 2.
For fixed λ ∈ R, let y1, y2, y3 be linearly independent solutions of (4.16), with Wronskian W .
Then W is constant in x and t, and γ = W−1/3(y1, y2, y3)
T is arclength-parametrized at each
time t, with centroaffine invariants k1 and k̃2 = k2 + λ. Furthermore, γ satisfies the evolution
equation
γt =
{
γ′′ − 2
3k1
γ, j = 1,
k1γ
′′ + (k̃2 − 4λ)γ′ + r0γ, j = 2.
Proof. If we let y = (y1, y2, y3)
T and form the matrix F = (y, y′, y′′), then F satisfies differential
equations of the form
F−1Fx =
0 0 k2 + k′1 + λ
1 0 k1
0 1 0
, F−1Ft = Nj ,
where both right-hand side matrices have trace zero. For example, when j = 1 one can directly
calculate, by differentiating yt =M1y, that
N1 =
−
2
3k1 k2 + 1
3k
′
1 + λ k′2 + 1
3k
′′
1
0 1
3k1 k2 + 2
3k
′
1 + λ
1 0 1
3k1
.
Thus, the Wronskian W is constant in x and t.
Because γ′′′ = (k1γ)′ + (k2 + λ)γ, the centroaffine invariants of γ are k1 and k̃2. It is
straightforward to compute γt in the j = 1 case, using yt =M1y. In the j = 2 case, we compute
yt =M2 y = k1y
′′ + (k̃2 − λ)y′ + r0y − 3DLy = k1y
′′ + (k̃2 − 4λ)y′ + r0y. �
4.4 Connection with Boussinesq equations
In [4] Chou and Qu note that, under the centroaffine curve flow (4.6), the curvatures k1, k2
satisfy a two-component system of evolution equations that is equivalent to the Boussinesq
equation. This suggests that the other integrable flow (4.11) under discussion may be related
to the Boussinesq hierarchy.
12 A. Calini, T. Ivey and G. Maŕı Beffa
Dickson et al. [5] write the (first) Boussinesq equation as a system
(q0)t + 1
6q
′′′
1 + 2
3q1q
′
1 = 0, (q1)t − 2q′0 = 0. (4.17)
They embed this in a hierarchy of integrable equations, each of which is written in Lax form as
Lt = [Pm, L], L := D3 + q1D + 1
2q
′
1 + q0, (4.18)
where Pm is a differential operator of order m 6≡ 0 mod 3, with coefficients depending on q0, q1
and their x-derivatives. Note that Pm must be chosen so that [Pm, L] has order one. For example,
while P1 = D yields the trivial evolution (q1)t = q′1, (q0)t = q′0, setting P2 = D2 + 2
3q1 gives the
Boussinesq equation (4.17).
Given the resemblance between (4.15) and (4.18), it is tempting to find substitutions to
connect the Boussinesq equation with (4.7). In fact, we can make L and L coincide by setting
k1 = −q1, k2 = 1
2q
′
1 − q0. (4.19)
With this substitution, M1 coincides with P2, so it follows that (4.7) and (4.17) are equivalent.
In [5] it is shown how the coefficients of the operators Pm can be obtained solving a recursive
system of differential equations, and thus these depend on a number of constants of integration.
For example, the expression for P4 is
P4 =
[
f1D
2 + (g1 − 1
2f
′
1)D + (16f
′′
1 − g′1 + 2
3q1f1)
]
+
+
[
f0D
2 + (g0 − 1
2f
′
0)D + (16f
′′
0 − g′0 + 2
3q1f0)
]
L+ k4,0 + k4,1L,
where f0 = 0, g0 = 1, f1 = 1
3q1 + c1, g1 = 1
3q0 + d1, and k4,0, k4,1, c1, d1 are arbitrary constants.
For convenience, we will set all these arbitrary constant to zero, so that
P4 = D4 + 4
3q1D
2 + 4
3(q′1 + q0)D + 5
9q
′′
1 + 2
3q
′
0 + 2
9q
2
1.
Again, if we use the substitutions (4.19), we find that the operator M2 coincides with −3P4.
Thus, (4.13) is equivalent to the second nontrivial flow in the Boussinesq hierarchy, provided we
also rescale time by t→ −3t.
5 Hierarchies
In [22] the Boussinesq hierarchy is discussed as an example of a bi-Hamiltonian system, in
which two sequences of commuting flows (and conservation laws) are generated by applying
recursion operators. Thus, given the equivalences established in Section 4.4, it is not surprising
that the Poisson operators defined in Section 4 can be combined to give a recursion operator
that generates a double hierarchy of commuting evolution equations for k1, k2. In fact, we will
show that our recursion operator is equivalent to the Boussinesq recursion operator as given in
Example 7.28 of [22]. The new information we add is that each of these evolution equations is
induced by a centroaffine geometric evolution equation for curves, which is itself Hamiltonian
relative to the pre-symplectic structure defined in Section 3.1 (see Theorem 5.4 below).
5.1 Recursion operators
We define a sequence of evolution equations for k1, k2
∂
∂tj
(
k1
k2
)
= Fj [k1, k2], (5.1)
Integrable Flows for Starlike Curves in Centroaffine Space 13
via the recursion
Fj+2 = PQ−1Fj , (5.2)
with initial data given by
F0 =
(
k′1
k′2
)
, F1 =
(
k′′1 + 2k′2
2
3(k1k
′
1 − k′′′1 )− k′′2
)
.
(Note that F1 is the right-hand side of (4.7), while for j = 0 (5.1) gives a simple transport
equation for k1, k2, corresponding to flow in the direction of the tangent vector γ′.)
In order to assert that the Fj defined by (5.2) are local functions of k1, k2 and their deriva-
tives – i.e., in calculating each Fj , the operator D−1 is only applied to exact x-derivatives of local
functions – we cite well-known results on the Boussinesq hierarchy. For example, the version of
the first Boussinesq equation used by Olver [22] is
uτ = v′, vτ = 1
3u
′′′ + 8
3uu
′, (5.3)
where τ is the time variable. If one considers linear transformations on the variables, it is
necessary to use some imaginary coefficients to make our version (4.7) of the first Boussinesq
equation for k1, k2 equivalent to (5.3):
x = x, τ = it, k1 = −2u, k2 = u′ − iv. (5.4)
Proposition 5.1. Under the above change of variables, the recursion operator PQ−1 is equiv-
alent to the Boussinesq recursion operator in [22].
Proof. The transformation between k1, k2 and u, v can be written as(
k1
k2
)
= G
(
u
v
)
, G :=
(
−2 0
D −i
)
.
Thus, if ∂/∂t (u, v)T = F [u, v] is an evolution equation for u, v, the right-hand side of the
corresponding evolution for k1, k2 is G ◦ F . Thus, our recursion operator PQ−1 for flows on
the k1, k2 variables corresponds to a recursion operator
G−1PQ−1G (5.5)
on flows in the u, v variables. In fact, when one calculates (5.5) and substitutes for k1, k2 in
terms of u, v using (5.4), the result is exactly −i times the Boussinesq recursion operator given
in [22]. �
Since in [27] (see Section 5.4 in that paper) it is proven that the Boussinesq recursion operator
from [22] always produces local flows when applied to the ‘seed’ evolution equations (i.e., the
tangent flow and first Boussinesq), it follows that the same is true for our recursion operator.
Remark 5.2. Once one checks that the evolution equations (5.1) for j = 0 and j = 1 commute, it
is automatic from the bi-Hamiltonian structure that all evolution equations in the sequence (5.1)
commute in pairs (see, e.g., Theorem 7.24 in [22]).
It is easy to check that the ‘seeds’ F0, F1 for the recursion are related to the initial conserved
densities by
F0 = PEρ0 = QEρ2, F1 = PEρ1 = QEρ3. (5.6)
14 A. Calini, T. Ivey and G. Maŕı Beffa
Table 1.
Non-stretching vector field Conserved density
Z0 = γ′ ρ0 = k1
Z1 = γ′′ − 2
3k1
γ ρ1 = k2
Z2 = k1γ
′′ + k2γ
′ + . . . ρ2 = k1k2
Z3 = (k′1 + 2k2)γ
′′ +
(
1
3k
2
1 − 2
3k
′′
1 − k′2
)
γ′ + . . . ρ3 = 1
3(k′1)
2 + k′1k2 + 1
9k
3
1 + k22
Z4 = (−k′′′1 − 2k′′2 + 2k1k
′
1 + 4k1k2)γ
′′ +
(
2
3k
(4)
1 + k′′′2 ρ4 = 1
3(k′′1)2 + k′′1(k′2 − k21)− k1(k′1)2
−2k1k
′′
1 − (k′1)
2 − 2k1k
′
2 + 4
9k
3
1 + 2k22
)
γ′ + . . . +(k′2)
2 − k21k′2 + 1
9k
4
1 + 2k1k
2
2
(The coefficient of γ in some vector fields is omitted for reasons of space, but can be determined from
the non-stretching condition.)
• The γ′ and γ′′ coefficients of Zj match the components of Eρj .
• Densities satisfy the recursion relation Eρj+2 = Q−1PEρj
• γt = Zj induces curvature evolution
(
k1
k2
)
t
= PEρj = QEρj+2.
• For j ≥ 2, Zj is a Hamiltonian vector field for −
∫
ρj−2 dx.
(The second set of equations was derived in Section 4.2.) While PQ−1 is the recursion operator
for commuting flows, it is evident from (5.6) that Q−1P should be the recursion operator for
conservation law characteristics (i.e., the result of applying the Euler operator E to a density).
In fact, we may define an infinite sequence of conserved densities by
Eρj+2 = Q−1PEρj , j ≥ 0. (5.7)
The first few densities calculated using this recursion appear in Table 1.
We now use these densities to define a sequence of flows for centroaffine curves, and relate
each of them to a curvature evolution equation in the sequence (5.1). Namely, if f is any local
function of k1, k2 and their derivatives, we define
Xf := (Ef)1γ
′ + (Ef)2γ
′′ + r0γ, (5.8)
where the subscripts indicate the components given by (4.8) and r0 is determined by the non-
stretching condition. Then for the sequence of densities defined recursively by (5.7) we define
the vector fields
Zj := Xρj (5.9)
and the corresponding sequence of curve flows
γt = Zj . (5.10)
Proposition 5.3. For each j ≥ 0 the curve flow (5.10) induces the curvature evolution
∂
∂t(k1, k2)
T = Fj.
Proof. From (5.6) and the recursion relations, it follows by induction that
Fj = PEρj , j ≥ 0.
Then the result follows immediately from (4.4). �
Integrable Flows for Starlike Curves in Centroaffine Space 15
5.2 Hamiltonian structure at the curve level
We now consider the question of how the Hamiltonian operator P is related to the Hamilto-
nian structure defined at the curve level in Section 3.1. Recall from Section 3.2 that X is
a Hamiltonian vector field associated to the functional H if
dH[Y ] = ωγ(Y,X)
for any non-stretching vector field Y .
Theorem 5.4. Let H(γ) =
∮
γ
ρdx and assume that
Eρ̂ = Q−1PEρ, (5.11)
i.e., ρ̂ is next after ρ in the sequence of densities generated by the recursion operator Q−1P.
Then −X ρ̂ (as defined by (5.8)) is Hamiltonian for H.
Proof. Based on the definition (3.1) of ω, we need to show that dH[Y ] equals∮
|Y, γ′,−X ρ̂|dx =
∮
|X ρ̂, γ′, Y |dx ∀Y ∈ TγM̂.
If X = aγ + bγ′ + cγ′′ and Y = ãγ + b̃γ′ + c̃γ′′, then from (3.2),∮
γ
|X, γ′, Y |dx =
∮
γ
(ac̃− ãc) dx.
However, using (2.4) to eliminate a and ã, we obtain∮
γ
|X, γ′, Y |dx =
∮
γ
(
−b′c̃+ c b̃′ + 1
3(c c̃′′ − c′′c̃)
)
dx = −
∮
γ
(b′c̃+ c′b̃) dx,
where the last equation follows by integration by parts. Thus,∮
γ
∣∣X ρ̂, γ′, Y
∣∣dx = −
∮
γ
(
b̃
c̃
)
· QEρ̂dx = −
∮
γ
(
b̃
c̃
)
· PEρdx,
using (5.11) in the last step. Then, because P is skew-adjoint,∮
γ
∣∣X ρ̂, γ′, Y
∣∣dx =
∮
γ
Eρ · P
(
b̃
c̃
)
dx.
On the other hand, using the properties of the Euler operator we have
dH[Y ] =
∮
γ
(Eρ)1δY k1 + (Eρ)2δY k2 dx,
where δY denotes the first variation in the direction of Y . Now using (4.4) we have
dH[Y ] =
∮
γ
Eρ ·
(
δY k1
δY k2
)
=
∮
γ
Eρ · P
(
b̃
c̃
)
dx.
This concludes the proof. �
The following corollaries are immediate consequences of the theorem.
16 A. Calini, T. Ivey and G. Maŕı Beffa
Corollary 5.5. Define the Poisson bracket
{H,G} =
∮
Eh · PEg dx,
where G(k1, k2) =
∮
g dx and H(k1, k2) =
∮
hdx, and g, h are functions of periodic k1, k2 and
their derivatives. Then dH[Xg] = {H,G}.
Corollary 5.6. The vector fields Zj defined by (5.9) are Hamiltonian for j ≥ 2.
Corollary 5.7. A closed curve γ is critical for the functional H(γ) =
∮
γ ρj dx with respect to
non-stretching variations if and only if γ is stationary for a constant-coefficient linear combina-
tion Zj+2 + c0Z0 + c1Z1.
Proof. By Theorem 5.4, −Zj+2 is Hamiltonian for H. Thus, a curve γ is H-critical if and only
if ω(Y, Zj+2) = 0 for all Y ∈ TγM̂ . This condition is satisfied if and only if Zj+2 is in the kernel
of ωγ , i.e., along γ it is equal to a constant-coefficient linear combination of vectors Z0 and Z1.
Equivalently, Zj+2 + c0Z0 + c1Z1 = 0 along γ for some constants c0, c1, expressing the property
that γ is stationary for a linear combination of these vector fields. �
5.3 Projective properties
As stated in Remark 2.2, a centroaffine curve γ projects to give a conic in RP2 if and only if the
Wilczynski invariants satisfy p0 − 1
2p
′
1 = 0. (The corresponding condition in terms of k1, k2 is
k2 + 1
2k
′
1 = 0.) In this subsection we will investigate flows in the hierarchy having the property
that, if γ projects to a conic at time zero, then it continues to have a conical projection at
subsequent times. We will show later that the equation of the conic in homogeneous coordinates
is fixed in time. We will also discuss flows that preserve a parametrization that is proportional
to projective arclength; in that case, the corresponding condition in terms of curvatures is that
k2 + 1
2k
′
1 is a nonzero constant along the curve.
These investigations are much easier if, instead of k2, we use an invariant that vanishes
precisely when the condition we are investigating holds. Accordingly, we fix a constant C, and
define an alternative pair of invariants
u = k1, v = k2 + 1
2k
′
1 − C. (5.12)
Thus, the curve is a conic if v = 0 when C = 0, and the curve has a constant-speed parametriza-
tion (relative to projective arclength) if v = 0 when C 6= 0.
We will convert the evolution equations in the hierarchy (at the level of the invariants) to
these variables. Suppose that a curve flow causes the invariants to evolve by(
k1
k2
)
t
= F [k1, k2],
where F is a vector-valued function of k1, k2 and their derivatives. Then the corresponding
evolution equation for the alternative invariants is(
u
v
)
t
= G ◦ F
[
u, v − 1
2u
′ + C
]
, G :=
(
1 0
1
2D 1
)
.
Similarly, if R is the recursion operator generating the hierarchy of evolution equations
for k1, k2, then the recursion operator for the corresponding flows on u, v differs by a gauge
transformation:
R̃ = GRG−1, where G−1 =
(
1 0
−1
2D 1
)
.
Integrable Flows for Starlike Curves in Centroaffine Space 17
(It is understood that, in R on the right-hand side, k1, k2 are substituted for in terms of u, v.)
Specifically, using R = PQ−1 as defined by (4.5), (4.10), we compute
R̃ = R̃0 + F̃0D
−1 (0 1
)
+ F̃1D
−1 (1 0
)
,
where
R̃0 =
(
3(v + C) 2
(
u−D2
)
N 3(v + C)
)
, F̃0 =
(
u′
v′
)
, F̃1 =
(
2v′
2
3uu
′ − 1
6u
′′′
)
,
and N is the scalar differential operator 1
6D
4 − 5
6uD
2 − 5
4u
′D+ 2
3u
2 − 3
4u
′′. One can check that
the vectors F̃0, F̃1 are the time derivatives of u, v, corresponding to the ‘seeds’ Z0 and Z1 for
the hierarchy of curve flows.
By applying R̃ to F̃0, F̃1, one can generate the right-hand sides of the evolution equations in
the hierarchy in terms of u and v. Letting F̃j denote these vectors, we compute (for example) that
F̃2 = R̃F̃0 =
(
−2v′′′ + 4(uv)′ + 3Cu′
1
6u
(5) − uu′′′ − 2u′u′′ + 4
3u
2u′ + (4v + 3C)v′
)
,
F̃3 = R̃F̃1 =
(
1
3u
(5) + 5
3(u2u′ − uu′′′)− 25
6 u
′u′′ + (10v + 6C)v′
1
3v
(5) −
(
5
6v + 1
2C
)
u′′′ + 5
3((u2v)′ − uv′′′ − u′′v′)− 5
2u
′v′′ + 2Cuu′
)
.
Here, when applying D−1 to differential polynomials in u, v, the constant of integration is taken
to be zero.
Notice in particular that if v ≡ 0, then the bottom component of F̃3 − 3CF̃1 vanishes. Thus,
the flow Z3 − 3CZ1 preserves the condition that v is identically zero. In fact, we can calculate
two infinite sequences of evolution equations for u, v that preserve this condition; the right-hand
sides of these are
Gk =
k∑
j=0
(
k
j
)
(−3C)jF̃2(k−j), k even,
k∑
j=0
(
k
j
)
(−3C)jF̃2(k−j)+1, k odd.
(5.13)
While it is routine to verify that any individual curvature evolution equation in these sequences
preserves v ≡ 0, it is easier to observe that the members of these sequences satisfy the recursion
relation Gk+2 = (R̃ − 3C)2Gk. Then the fact that they all preserve v ≡ 0 is a consequence of
the following:
Proposition 5.8. If a curvature evolution (ut, vt)
T = Gk[u, v] in this sequence preserves v ≡ 0,
then so does the evolution (ut, vt)
T = Gk+2[u, v].
Proof. We assume that Gk = (D`1, D`2)
T for local functions `1, `2 of u, v and their derivatives.
(This form for Gk is necessary if we are able to apply operator R̃ to it and produce local
functions.)
Within the ring of polynomials in u, v and their derivatives, let V denote the ideal generated
by v, v′, v′′, etc. By hypothesis, D`2 ∈ V, and the same is true for `2.
We compute
(R̃ − 3C)Gk = `1F̃1 + `2F̃0 +
(
3vD`1 + 2
(
u−D2
)
D`2
ND`1 + 3vD`2
)
. (5.14)
18 A. Calini, T. Ivey and G. Maŕı Beffa
Thus, the bottom component of (R̃ − 3C)2Gk is given by(
N + F̃12D
−1)(`1F̃11 + `2F̃01 + 3vD`1 + 2(u−D2)D`2
)
+
(
3v + F̃02D
−1)(`1F̃12 + `2F̃02 +ND`1 + 3vD`2
)
.
(Here, F̃j1 and F̃j2 denote the top and bottom entries in the vector F̃j .) The second factor in
the top line is the top entry of (R̃ − 3C)Gk. This polynomial lies in V, and the same is true if
we apply N or D−1 to it. On the other hand, because F̃02 = v′ ∈ V, the coefficient in front on
the second line also vanishes when v ≡ 0. �
In the special case when C = 0, we see that the following curve flows (as defined by (5.9))
preserve conicity:
Z0, Z3, Z4, Z7, Z8, Z11, Z12, . . . . (5.15)
However, when C 6= 0, some care needs to be taken in matching the evolution equations for u, v
that preserve v ≡ 0 with the corresponding linear combinations of the curve flows Zj .
In the proof of Proposition 5.8, we used the fact that if an exact derivative D` lies in V, then
by choosing the constant of integration equal to zero, the antiderivative ` also lies in V. However,
if we express D` in terms of k1, k2 instead of u, v, and then take an antiderivative, a particular
constant of integration must be chosen in order to belong in V. Thus, when we compute the kth
evolution equation for u, v by using the recursion operator R̃ (which involves applying D−1),
then convert this to an evolution equation for k1, k2, and finally try to match it with a curve
flow in the hierarchy (5.10), we get a linear combination of Zk with lower-order flows of the same
parity. For example, if we substitute (5.12) into F̃2, and then apply the operator G−1, we get
G−1F̃2
[
k1, k2+ 1
2k
′
1−C
]
=
(
−k′′′′1 −2k′′′2 +2k1k
′′
1 +4(k1k2)
′+2(k′1)
2−Ck′1
2
3k
′′′′′
1 +k′′′′2 −2k1k
′′′
1 −4k′1k
′′
1 +4k2k
′
2−2(k1k
′
2)
′+ 4
3(k1)
2k′1−Ck′2
)
= F2 − CF0.
Thus, F̃2 is induced by the curve flow Z2−CZ0; similarly, F̃3 is induced by Z3−2CZ1, and so on.
Similarly, when we apply the recursion operator (R̃ − 3C)2 to generate higher-order flows
that preserve v ≡ 0, these constants of integration accumulate and change the relatively nice
pattern of the coefficients exhibited by (5.13). Here is what we get when we compute the first
few curve evolutions corresponding to the evolutions Gk:
u, v evolution centroaffine curve flow
G0 Z0
G1 Z3 − 5CZ1
G2 Z4 − 7CZ2 + 14C2Z0
G3 Z7 − 11CZ5 + 44C2Z3 − 220
3 C
3Z1
G4 Z8 − 13CZ6 + 65C2Z4 − 455
3 C
3Z2 + 455
3 C
4Z0
5.4 Conical evolutions and the Kaup–Kuperschmidt hierarchy
In this section we examine special properties of the conicity-preserving flows (5.15) (hence, from
now on we are assuming C = 0). These properties will enable us to connect our hierarchy of
centroaffine curve flows in R3 with the Kaup–Kuperschmidt hierarchy and with curve flows
in centroaffine R2. We begin with the observation that, when restricted to conical curves, the
coefficient of γ′′ vanishes for as many of the vector fields in (5.15) as one cares to check. In other
words, this coefficient belongs to V, the ideal within the ring of differential polynomials in k1, k2
generated by k2+ 1
2k
′
1 and its derivatives. In fact, this is true in general, as shown in the following:
Integrable Flows for Starlike Curves in Centroaffine Space 19
Proposition 5.9. If j ≡ 0 or j ≡ −1 modulo 4, then the bottom component of Eρj belongs in V;
hence, the γ′′ coefficient of Zj vanishes on conical curves.
Proof. The statement can be verified directly for j = 0 and j = 3. For higher values, we use the
recursion relation between the characteristics, which implies that Eρj+4 = (Q−1P)2Eρj . From
Proposition 5.8 we know that the bottom component of F̃j = GPEρj lies in V. By inserting
powers of G and G−1 into the recursion relation, we get
Eρj+4 = Q−1PQ−1G−1GPEρj = Q−1RG−1F̃j .
As in the proof of Proposition 5.8, we can assume that F̃j = (D`1, D`2)
T where `2 ∈ V. Taking
C = 0 in equation (5.14), we see that the top entry of RG−1F̃j = G−1R̃F̃j is
`1F̃11 + `2F̃01 + 3vD`1 + 2
(
u−D2
)
D`2, (5.16)
which clearly is in V. (Note that v = k2 + 1
2k
′
1 here.) Noting the form of Q, we see that the
bottom entry of Eρj+4 is D−1 applied to (5.16), so it must also belong to V. �
Next, we make a connection with curve flows in centroaffine R2 to show that, for the flows
in (5.15), the cone that the curve lies on is preserved by the time evolution.
Proposition 5.10. If γ(x, t) evolves by any of the vector fields in (5.15), and γ(x, 0) lies on
a cone through the origin in R3, then γ(x, t) lies on the same cone at later times.
Proof. Using the action of SL(3) we can, without loss of generality, assume that the equation
of the cone is y1y3 − (y2)
2 = 0. We fix a map V from R2 onto this cone:
V
(
x1
x2
)
= 2−1/3
x21
x1x2
x22
.
Of course, when we projectivize on each end this gives the Veronese embedding of RP1 as
a quadric in RP2. The scale factor of 2−1/3 is chosen so that if X(x) is a parametrized curve in R2
satisfying the centroaffine normalization |X,X ′| = 1, then Γ = V◦X satisfies the normalization
|Γ,Γ′,Γ′′| = 1. Moreover, if p(x) is the curvature of X, then the invariants of Γ are k1 = −4p
and k2 = 2p′. Finally, if X evolves by the non-stretching flow
Xt = rX ′ − 1
2r
′X, (5.17)
then Γ(x, t) = V ◦X(x, t) satisfies Γt = rΓ′ − r′Γ.
By Proposition 5.9 all the flows in (5.15), when restricted to conical curves, are of this form,
for some choice of differential polynomial r in k1. For any initial data γ(x, 0), we can define
a curve X(x, 0) in R2 such that γ(x, 0) = V◦X(x, 0) and make X(x, t) evolve by (5.17). Because
Γ(x, t) = V ◦X(x, t) satisfies the same initial value problem, then γ(x, t) = Γ(x, t) at all times,
and γ(x, t) lies on the cone defined by y1y3 − (y2)
2 = 0 at all times. �
In [4] Chou and Qu discovered a non-stretching flow for curves in centroaffine R3 which
preserves the conicity condition k2 + 1
2k
′
1 = 0 and which causes the curvature k1 to evolve by
the Kaup–Kuperschmidt equation:
ut = u′′′′′ − 5uu′′′ − 25
2 u
′u′′ + 5u2u′
(see Case 3 in Section 3 of their paper, taking λ = 0). In fact, up to a multiplicative factor
of 1/3, Chou and Qu’s flow is the same as the restriction of Z3 to conical curves.
20 A. Calini, T. Ivey and G. Maŕı Beffa
Not only does flow Z3 give a geometric realization of the Kaup–Kuperschmidt equation, but
the entire sequence (5.15) of flows realizes the Kaup–Kuperschmidt hierarchy, when restricted to
conical curves. To see this, note that the square of the recursion operator R̃ relates the evolution
of k1 = u under Zj to its evolution under Zj+4. (Here, we use the notation of Section 5.3 but
with C = 0 and v = 0 because of conicity.) The resulting recursion operator is
−1
3D
6 + 2uD4 + 6u′D3 +
(
49
6 u
′′−3u2
)
D2 +
(
35
6 u
′′′−10uu′′
)
D + 13
6 u
′′′′ − 41
6 uu
′′
− 23
4 (u′)2 + 4
3u
3 + u′D−1 ◦ (13u
2− 1
6u
′′) + 1
3
(
u′′′′′−5uu′′′− 25
2 u
′u′′+5u2u′
)
D−1.
This agrees with the known recursion operator for symmetries of the Kaup–Kuperschmidt hierar-
chy. (See, e.g., Example 2.20 in [29], where the operator differs by changing u to −u and rescaling
time by a factor of 1/3.) Using this, one can check that the curvature flows induced by (5.15)
for conical curves coincide with the commuting flows of the Kaup–Kuperschmidt hierarchy.
Remark 5.11. The curve flow discovered by Chou and Qu is in fact also defined for centroaffine
curves parametrized proportional to projective arclength (i.e., those for which k2 + 1
2k
′
1 is a con-
stant), and nevertheless still induces Kaup–Kuperschmidt evolution for k1. Recently, Musso [20]
has extended this to a hierarchy of flows for arclength-parametrized curves in RP2 which induces
the Kaup–Kuperschmidt hierarchy for curvature evolution. We suspect that these flows coincide
with the restrictions of the flows studied in Section 5.3 (for C 6= 0) to the centroaffine lifts of
such curves in RP2.
Remark 5.12. Schwartz and Tabachnikov [28] showed that certain maps defined on the space
of convex polygons preserve the subset of polygons that are inscribed (or circumscribed) on
a conic: that is, if the vertices of the polygon (or those of its projective dual) lie on a conic, then
the same is true for its image under the map. The building blocks for these maps are elementary
maps Tr that associate to a given polygon another polygon obtained from the intersections of
diagonals joining each vertex to the vertex located r positions to the left or right of it. In fact,
the maps preserving conicity are particular combinations of Tr for certain values of r.
The map corresponding to r = 2 is called the pentagram map and it is known to be a dis-
cretization (in both time and space) of the Boussinesq equation [24]. It is natural to wonder if
the maps in [28] are somehow associated to flows in the Boussinesq hierarchy. We are currently
investigating this.
Acknowledgements
The authors would like to acknowledge the careful and detailed work of the anonymous referees
in helping improve this paper. The authors also gratefully acknowledge support from the Na-
tional Science Foundation: A. Calini through grants DMS-0608587 and DMS-1109017, and as
a current NSF employee; T. Ivey through grant DMS-0608587; and G. Maŕı Beffa through grant
DMS-0804541. T. Ivey also acknowledges support from the College of Charleston Mathematics
Department. G. Maŕı Beffa also acknowledges the support of the Simons Foundation through
their Fellows program.
References
[1] Arnold V.I., Khesin B.A., Topological methods in hydrodynamics, Applied Mathematical Sciences, Vol. 125,
Springer-Verlag, New York, 1998.
[2] Brylinski J.-L., Loop spaces, characteristic classes and geometric quantization, Progress in Mathematics,
Vol. 107, Birkhäuser Boston Inc., Boston, MA, 1993.
[3] Calini A., Ivey T., Maŕı Beffa G., Remarks on KdV-type flows on star-shaped curves, Phys. D 238 (2009),
788–797, arXiv:0808.3593.
http://dx.doi.org/10.1016/j.physd.2009.01.007
http://arxiv.org/abs/0808.3593
Integrable Flows for Starlike Curves in Centroaffine Space 21
[4] Chou K.S., Qu C., Integrable motions of space curves in affine geometry, Chaos Solitons Fractals 14 (2002),
29–44.
[5] Dickson R., Gesztesy F., Unterkofler K., Algebro-geometric solutions of the Boussinesq hierarchy, Rev.
Math. Phys. 11 (1999), 823–879, solv-int/9809004.
[6] Doliwa A., Santini P.M., An elementary geometric characterization of the integrable motions of a curve,
Phys. Lett. A 185 (1994), 373–384.
[7] Dorfman I., Dirac structures and integrability of nonlinear evolution equations, Nonlinear Science: Theory
and Applications, John Wiley & Sons Ltd., Chichester, 1993.
[8] Drinfel’d V.G., Sokolov V.V., Lie algebras and equations of Korteweg–de Vries type, J. Sov. Math. 30
(1985), 1975–2036.
[9] Hasimoto R., A soliton on a vortex filament, J. Fluid Mech. 51 (1972), 477–485.
[10] Huang R., Singer D.A., A new flow on starlike curves in R3, Proc. Amer. Math. Soc. 130 (2002), 2725–2735.
[11] Lamb G.L., Solitons on moving space curves, J. Math. Phys. 18 (1977), 1654–1661.
[12] Langer J., Recursion in curve geometry, New York J. Math. 5 (1999), 25–51.
[13] Langer J., Perline R., Local geometric invariants of integrable evolution equations, J. Math. Phys. 35 (1994),
1732–1737, solv-int/9401001.
[14] Langer J., Perline R., Poisson geometry of the filament equation, J. Nonlinear Sci. 1 (1991), 71–93.
[15] Maŕı Beffa G., Bi-Hamiltonian flows and their realizations as curves in real semisimple homogeneous mani-
folds, Pacific J. Math. 247 (2010), 163–188.
[16] Maŕı Beffa G., Poisson brackets associated to the conformal geometry of curves, Trans. Amer. Math. Soc.
357 (2005), 2799–2827.
[17] Maŕı Beffa G., The theory of differential invariants and KdV Hamiltonian evolutions, Bull. Soc. Math.
France 127 (1999), 363–391.
[18] Marsden J., Weinstein A., Coadjoint orbits, vortices, and Clebsch variables for incompressible fluids, Phys. D
7 (1983), 305–323.
[19] Mokhov O., Symplectic and Poisson geometry on loop spaces of manifolds and nonlinear equations, in
Topics in Topology and Mathematical Physics, Amer. Math. Soc. Transl. Ser. 2, Vol. 170, Amer. Math.
Soc., Providence, RI, 1995, 121–151, hep-th/9503076.
[20] Musso E., Motions of curves in the projective plane inducing the Kaup–Kupershmidt hierarchy, SIGMA 8
(2012), 030, 20 pages, arXiv:1205.5329.
[21] Nakayama K., Motion of curves in hyperboloid in the Minkowski space, J. Phys. Soc. Japan 67 (1998),
3031–3037.
[22] Olver P.J., Applications of Lie groups to differential equations, Graduate Texts in Mathematics, Vol. 107,
2nd ed., Springer-Verlag, New York, 1993.
[23] Olver P.J., Moving frames and differential invariants in centro-affine geometry, Lobachevskii J. Math. 31
(2010), 77–89.
[24] Ovsienko V., Schwartz R., Tabachnikov S., The pentagram map: A discrete integrable system, Comm.
Math. Phys. 299 (2010), 409–446, arXiv:0810.5605.
[25] Ovsienko V., Tabachnikov S., Projective differential geometry old and new. From the Schwarzian deriva-
tive to the cohomology of diffeomorphism groups, Cambridge Tracts in Mathematics, Vol. 165, Cambridge
University Press, Cambridge, 2005.
[26] Pinkall U., Hamiltonian flows on the space of star-shaped curves, Results Math. 27 (1995), 328–332.
[27] Sanders J.A., Wang J.P., Integrable systems and their recursion operators, Nonlinear Anal. 47 (2001),
5213–5240.
[28] Schwartz R.E., Tabachnikov S., Elementary surprises in projective geometry, Math. Intelligencer 32 (2010),
31–34, arXiv:0910.1952.
[29] Wang J.P., A list of 1 + 1 dimensional integrable equations and their properties, J. Nonlinear Math. Phys.
9 (2002), suppl. 1, 213–233.
[30] Wilczynski E.J., Projective differential geometry of curves and ruled surfaces, B.G. Teubner, Leipzig, 1906.
http://dx.doi.org/10.1016/S0960-0779(01)00179-5
http://dx.doi.org/10.1142/S0129055X9900026X
http://dx.doi.org/10.1142/S0129055X9900026X
http://arxiv.org/abs/solv-int/9809004
http://dx.doi.org/10.1016/0375-9601(94)90170-8
http://dx.doi.org/10.1007/BF02105860
http://dx.doi.org/10.1017/S0022112072002307
http://dx.doi.org/10.1090/S0002-9939-02-06631-5
http://dx.doi.org/10.1063/1.523453
http://dx.doi.org/10.1063/1.530567
http://arxiv.org/abs/solv-int/9401001
http://dx.doi.org/10.1007/BF01209148
http://dx.doi.org/10.2140/pjm.2010.247.163
http://dx.doi.org/10.1090/S0002-9947-04-03589-5
http://dx.doi.org/10.1016/0167-2789(83)90134-3
http://arxiv.org/abs/hep-th/9503076
http://dx.doi.org/10.3842/SIGMA.2012.030
http://arxiv.org/abs/1205.5329
http://dx.doi.org/10.1143/JPSJ.67.3031
http://dx.doi.org/10.1007/978-1-4612-4350-2
http://dx.doi.org/10.1134/S1995080210020010
http://dx.doi.org/10.1007/s00220-010-1075-y
http://dx.doi.org/10.1007/s00220-010-1075-y
http://arxiv.org/abs/0810.5605
http://dx.doi.org/10.1016/S0362-546X(01)00630-7
http://dx.doi.org/10.1007/s00283-010-9137-8
http://arxiv.org/abs/0910.1952
http://dx.doi.org/10.2991/jnmp.2002.9.s1.18
1 Introduction
1.1 Integrable evolutions of space curves
1.2 Pinkall's flow in R2
1.3 Organization of the paper
2 Centroaffine curve flows in R3
2.1 Centroaffine invariants
2.2 Non-stretching variations
3 Hamiltonian curve flows
3.1 Symplectic structure on starlike loops
3.2 Examples
4 Integrable centroaffine curve flows
4.1 Evolution of invariants
4.2 Two integrable flows
4.3 Lax representation
4.4 Connection with Boussinesq equations
5 Hierarchies
5.1 Recursion operators
5.2 Hamiltonian structure at the curve level
5.3 Projective properties
5.4 Conical evolutions and the Kaup-Kuperschmidt hierarchy
References
|