We propose a fourth-order extension of the geodesic problem arising in the continuity of the study of Riemannian cubics. We consider the variational problem in a Riemannian manifold and derive the Euler-Lagrange equation. For the special situation of Lie groups, we use Euler-Poincaré reduction to obtain the reduced equation on the Lie algebra.
Citation: Margarida Camarinha. A natural 4th-order generalization of the geodesic problem[J]. Electronic Research Archive, 2024, 32(5): 3396-3412. doi: 10.3934/era.2024157
[1] | Yuriĭ G. Nikonorov, Irina A. Zubareva . On the behavior of geodesics of left-invariant sub-Riemannian metrics on the group $ \operatorname{Aff}_{0}(\mathbb{R}) \times \operatorname{Aff}_{0}(\mathbb{R}) $. Electronic Research Archive, 2025, 33(1): 181-209. doi: 10.3934/era.2025010 |
[2] | Guoliang Ju, Can Chen, Rongliang Chen, Jingzhi Li, Kaitai Li, Shaohui Zhang . Numerical simulation for 3D flow in flow channel of aeroengine turbine fan based on dimension splitting method. Electronic Research Archive, 2020, 28(2): 837-851. doi: 10.3934/era.2020043 |
[3] | Vladimir Rovenski . Willmore-type variational problem for foliated hypersurfaces. Electronic Research Archive, 2024, 32(6): 4025-4042. doi: 10.3934/era.2024181 |
[4] | Alessandro Portaluri, Li Wu, Ran Yang . Linear instability of periodic orbits of free period Lagrangian systems. Electronic Research Archive, 2022, 30(8): 2833-2859. doi: 10.3934/era.2022144 |
[5] | Fei Shi . Incompressible limit of Euler equations with damping. Electronic Research Archive, 2022, 30(1): 126-139. doi: 10.3934/era.2022007 |
[6] | Min Li, Xueke Pu, Shu Wang . Quasineutral limit for the compressible two-fluid Euler–Maxwell equations for well-prepared initial data. Electronic Research Archive, 2020, 28(2): 879-895. doi: 10.3934/era.2020046 |
[7] | Xuerong Hu, Yuxiang Han, Junyan Lu, Linxiang Wang . Modeling and experimental investigation of quasi-zero stiffness vibration isolator using shape memory alloy springs. Electronic Research Archive, 2025, 33(2): 768-790. doi: 10.3934/era.2025035 |
[8] | Noelia Bazarra, José R. Fernández, Ramón Quintanilla . Numerical analysis of a problem in micropolar thermoviscoelasticity. Electronic Research Archive, 2022, 30(2): 683-700. doi: 10.3934/era.2022036 |
[9] | Li-Bin Liu, Ying Liang, Jian Zhang, Xiaobing Bao . A robust adaptive grid method for singularly perturbed Burger-Huxley equations. Electronic Research Archive, 2020, 28(4): 1439-1457. doi: 10.3934/era.2020076 |
[10] | Michael Barg, Amanda Mangum . Statistical analysis of numerical solutions to constrained phase separation problems. Electronic Research Archive, 2023, 31(1): 229-250. doi: 10.3934/era.2023012 |
We propose a fourth-order extension of the geodesic problem arising in the continuity of the study of Riemannian cubics. We consider the variational problem in a Riemannian manifold and derive the Euler-Lagrange equation. For the special situation of Lie groups, we use Euler-Poincaré reduction to obtain the reduced equation on the Lie algebra.
In the last few decades, there has been a lot of interest in the study of higher-order variational problems on Riemannian manifolds. One of the most prominent second-order variational problems has the Lagrangian function given by the squared norm of the covariant acceleration. Stationary curves for this problem are the so-called Riemannian cubics or bi-harmonic curves [1,2]. These curves have played an important role in the study of interpolation problems [4,5,6]. Approaches to this second order problem are not restricted to its variational character, it has also been studied as an optimal control problem in which the control system is an affine connection system [7,8,9,10,11,12]. From the point of view of applications, this problem can be seen as an optimal navigation problem with regard to fuel consumption, and may involve several agents, have obstacles to avoid or be conditioned by constraints [13,14,15,16,17].
The theory of Riemannian cubics can be interpreted as a higher-order extension of the theory of geodesics. An extension of this theory to odd-degree Riemannian polynomials has been developed in [18,19,20] (for a more general setting, see also [21] and references therein). Although the concept of odd-degree Riemannian polynomials given in [18] can be seen as a higher-order generalization of the concept of geodesics, it cannot be said to be a generalization of lower-degree polynomials, starting with the cubic ones. Indeed, geodesics are particular cases of higher-degree Riemannian polynomials, but Riemannian cubics cannot in general be seen as higher-degree Riemannian polynomials.
The notion of Riemannian polynomials given in [18] was motivated by the optimization properties of the Euclidean polynomials and, thus, can be seen as a generalization of polynomials from Euclidean spaces to Riemannian manifolds. These polynomials play an important role in the study of interpolation and data fitting problems in curved spaces. Applications range from the control theory of mechanical systems to image processing or statistical analysis of shape data. The polynomials are used to construct interpolation curves (Riemannian splines) that give trajectories of physical systems along some prescribed positions [6,20,22], or they can be combined with least squares fitting to better represent image (or shape) data [23,24]. The degree of the polynomials provides a certain smoothness of the splines and accuracy of the represented data.
In this paper, we propose studying a variational problem of order 4 which gives rise to a definition of higher-degree polynomials that is different from the one presented in [18]. The stationary curves of the proposed problem can be interpreted as the best approximations of Riemannian cubics when their existence is conditioned by boundary or interpolation conditions, and reduce naturally to Riemannian cubics when such curves exist.
The study of this variational problem in the rotation group SO(3) was presented in [25]. In this paper, we consider the problem in a general Riemannian manifold and derive the Euler-Lagrange equation (Section 3). Then, in Section 4, we analyse what happens on Lie groups. Since the Lagrangian function is invariant, we use the higher-order Euler-Poincaré reduction presented in [26] to obtain the corresponding equation on the Lie algebra.
In what follows (M,⟨.,.⟩) is an n-dimensional Riemannian manifold. The Levi-Civita connection is the mapping ∇, which assigns to each two smooth vector fields X and Y in M, a new vector field, ∇XY, satisfying the following properties:
(i) ∇ is IR-bilinear in X and Y,
(ii) ∇fXY=f∇XY,
(iii) ∇X(fY)=(Xf)Y+f∇XY,
(iv) ∇XY is smooth in M,
(v) ∇XY−∇YX=[X,Y] (Symmetry),
(vi) X⟨Y,Z⟩=⟨∇XY,Z⟩+⟨Y,∇XZ⟩ (Compatibility),
where X, Y and Z are smooth vector fields on M, f a smooth real valued function on M and [X,Y] denotes the Lie bracket of the vector fields X and Y defined by
[X,Y]f=X(Yf)−Y(Xf). |
The covariant derivative of a vector field X along a curve x in M is given by (DX/dt)(t)=∇(dx/dt)(t)X(t), where dx/dt is the velocity vector field along x. The jth-order covariant derivative DjX/dtj is the covariant derivative of the vector field Dj−1X/dtj−1, where j≥1. The jth-order covariant derivative Djx/dtj of a curve x is the (j−1)th-order covariant derivative of the velocity vector field along x, for j≥2.
For any vector field X along a 2-parameter curve x:(r,t)→x(r,t), the covariant derivatives DX/∂r and DX/∂t are vector fields along x defined as follows: We fix t (resp., r) and consider X restricted to the curve xt:r↦x(r,t) (resp., xr:t↦x(r,t)), which leads to a vector field along xt (resp., xr). The covariant derivative of the vector field X along the curve xt (resp., xr) gives rises to the vector field DX/∂r (resp., DX/∂t) along x. As examples, we can set X equal to ∂x/∂r or ∂x/∂t and obtain the following higher-order covariant derivatives,
D∂r∂x∂r,D∂r∂x∂t,D∂t∂x∂r,D∂t∂x∂t, |
or even higher-order ones. Note that, from the symmetry of the Levi-Civita connection (condition (5)), we get immediately that
D∂r∂x∂t=D∂t∂x∂r. | (2.1) |
To deal with higher-order partial covariant derivatives, we shall need to define the curvature tensor.
The curvature tensor R is given by the identity
R(X,Y)Z=∇X∇YZ−∇Y∇XZ−∇[X,Y]Z, | (2.2) |
where X, Y, Z are vector fields in M. Therefore, R measures how much covariant derivatives fail to commute. For any vector field X along a 2-parameter curve x:(r,t)→x(r,t), we have
D∂rD∂tX−D∂tD∂rX=R(∂x∂r,∂x∂t)X. | (2.3) |
From Eq (2.3) we see that R is an essential tool to interchange the order of partial covariant derivatives. For this reason, the following symmetries of R play such an important role in this paper.
For all vector fields X, Y, Z and W in M, the curvature tensor R satisfies the skew-symmetry identities
R(X,Y)Z+R(Y,X)Z=0, | (2.4) |
and
⟨R(X,Y)Z,W⟩+⟨R(X,Y)W,Z⟩=0, | (2.5) |
and, as a consequence of (2.4) and (2.5), the symmetry by pairs identity
⟨R(X,Y)Z,W⟩=⟨R(W,Z)Y,X⟩. | (2.6) |
The cyclic permutation property
R(X,Y)Z+R(Y,Z)X+R(Z,X)Y=0, | (2.7) |
is called the first Bianchi identity.
The covariant differential of the curvature tensor R is defined by
(∇XR)(Y,Z)W=∇X(R(Y,Z)W)−R(∇XY,Z)W−R(Y,∇XZ)W−R(Y,Z)∇XW, | (2.8) |
and the symmetry properties of R lead to the following equalities.
(∇XR)(Y,Z)W+(∇XR)(Z,Y)W=0, | (2.9) |
(∇XR)(Y,Z)W+(∇XR)(Z,W)Y+(∇XR)(W,Y)Z=0, | (2.10) |
⟨(∇XR)(Y,Z)W,V⟩=⟨(∇XR)(V,W)Z,Y⟩, | (2.11) |
⟨(∇XR)(Y,Z)W,V⟩+⟨(∇XR)(Y,Z)V,W⟩=0, | (2.12) |
and the second Bianchi identity
(∇XR)(Y,Z)W+(∇YR)(Z,X)W+(∇ZR)(X,Y)W=0, | (2.13) |
for every vector field X, Y, Z, W and V on M.
More details about Riemannian geometry can be found in Milnor [27], Lee [28] and O'Neill [29].
Let x be a C3-piecewise smooth curve on M and Γ(x) the vector field along x defined by
Γ(x)=D4xdt4+R(D2xdt2,dxdt)dxdt. | (3.1) |
It is well-known that the curves x satisfying
Γ(x)=0, | (3.2) |
are the so-called Riemannian cubics. They arose as a natural generalization of geodesics. In fact, the Riemannian cubics are the extremal curves of the bienergy functional
B(x)=12∫T0‖D2xdt2‖2dt, | (3.3) |
and this functional is seen as a second order extension of the energy functional
E(x)=12∫T0‖dxdt‖2dt. | (3.4) |
Geodesics are the curves x satisfying
D2xdt2=0, | (3.5) |
and are critical points of the energy functional (3.3).
The variational problem that consists in finding the curves minimizing the bienergy functional (3.3) is called the Riemannian cubic problem and appears as a natural extension of the geodesic problem. It can be formulated as a optimal control problem in the following way.
minU∫T012⟨U,U⟩dt, | (3.6) |
subject to the control system
dxdt=V,DVdt=U, | (3.7) |
and the boundary conditions
x(0)=x0,x(T)=x1,V(0)=v01,V(T)=v11. | (3.8) |
Studying this problem from the point of view of Hamiltonian systems and optimal control proved to be very interesting and led to fruitful research [7,8,10,12].
Remark 1. Riemannian cubics and geodesics can also be seen as generalizations of Euclidean polynomials (of degree 1 and 3, respectively). This simple interpretation is based on optimization properties of these Euclidean curves. In fact, these properties are not restricted to polynomials of degree 1 and 3, but are valid for arbitrary odd-degree Euclidean polynomials. It is well-know that Euclidean polynomials x of degree 2m−1 satisfy the inequality
12∫T0‖dmxdtm‖2dt≤12∫T0‖dmydtm‖2dt, | (3.9) |
among all Cm functions y satisfying the boundary conditions
y(0)=x0,y(T)=x1, | (3.10) |
and
djydtj(0)=v0j,djydtj(T)=v1j,1≤j≤m−1; | (3.11) |
moreover, the equality occurs iff y coincides with x. The differential equation
d2mxdt2m=0, | (3.12) |
is the Euler-Lagrange equation of the m-energy functional
12∫T0‖dmxdtm‖2dt. |
Seen as a solution of Eq (3.12), each Euclidean polynomial is contained in the class of Euclidean polynomial of higher-order (Property A).
Remark 2. In [18], the Euclidean variational problem of order m is extended to Riemannian manifolds by considering the following generalization of the m-energy functional
12∫T0‖Dmxdtm‖2dt. |
This variational problem gives rise to a notion of odd-degree Riemannian polynomials, for which, in general, property A fails, even when the curvature is constant.
To extend the Riemannian cubic problem we formulate the following fourth-order variational problem.
Consider the functional J defined by
J(x)=12∫T0‖Γ(x)‖2dt. | (3.13) |
The Riemannian heptic problem (P) consists in finding the curves minimizing the functional J among the set Ω of all C3-piecewise smooth curves x satisfying the boundary conditions
x(0)=x0,x(T)=x1, | (3.14) |
and
Djxdtj(0)=v0j,Djxdtj(T)=v1j,1≤j≤3, | (3.15) |
where vij∈TxiM, with i=0,1 and 1≤j≤3.
The variational problem (P) can be formulated as an optimal control problem as follows:
minU∫T012⟨U,U⟩dt, | (3.16) |
subject to the control system
dxdt=V,DVdt=Z,DZdt=W,DWdt+R(Z,V)V=U, | (3.17) |
and the boundary conditions
x(0)=x0,x(T)=x1,V(0)=v01,V(T)=v11,Z(0)=v02,Z(T)=v12,W(0)=v03,W(T)=v13. | (3.18) |
When studying optimal control problems for these higher-order dynamic systems, we encounter even higher-order bundles than in the Riemannian cubic problem.
An admissible variation of the curve x∈Ω is a C3-piecewise smooth one-parameter variation
α:(−ϵ,ϵ)×[0,T]→M(r,t)↦α(r,t)=αr(t), |
such that αr∈Ω, for all r∈(−ϵ,ϵ), ϵ>0. The variational vector field X associated with an admissible variation α of x is given by
X(t)=∂α∂r(0,t),t∈[0,T]. |
Admissible variations of x give rise to the real vector space TxΩ of all the C3 piecewise smooth vector field X along x verifying the boundary conditions
X(0)=0,X(T)=0, | (3.19) |
and
DjXdtj(0)=DjXdtj(T)=0,1≤j≤3. | (3.20) |
The first variation δJ(x) of J at the curve x is the linear form on TxΩ defined by
δJ(x)(X)=ddrJ(αr)|r=0, |
where α is an admissible variation of x with associated variational vector field X.
Necessary conditions for x∈Ω to be a local minimizer of J follow from the Hamilton variational principle through the first variation of J at x,
δJ(x)(X)=0,∀X∈TxΩ, | (3.21) |
which gives rise to the Euler-Lagrange equations for the problem (P).
To obtain the first variation formula of J at x let us consider an admissible variation α of x∈Ω.
First of all, we have to compute (D/∂r)Γ(α). To begin, note that, although we know from (2.1) that
D∂r∂α∂t=D∂t∂α∂r, | (3.22) |
we need to use the curvature tensor R to measure how much the second covariant derivative (D/∂r)(D/∂t)Y is symmetric in r and t, for Y an arbitrary vector field along α. From (2.3), we know that
D∂rD∂tY−D∂tD∂rY=R(∂α∂r,∂α∂t)Y. | (3.23) |
Making use of these tools to interchange the order of partial covariant derivatives, we obtain
D∂rD2α∂t2=D2dt2∂α∂r+R(∂α∂r,∂α∂t)∂α∂t,D∂rD3α∂t3=D3∂t3∂α∂r+D∂t(R(∂α∂r,∂α∂t)∂α∂t)+R(∂α∂r,∂α∂t)D2α∂t2,D∂rD4α∂t4=D4∂t4∂α∂r+D2∂t2(R(∂α∂r,∂α∂t)∂α∂t)+D∂t(R(∂α∂r,∂α∂t)D2α∂t2)+R(∂α∂r,∂α∂t)D3α∂t3,D∂r(R(D2α∂t2,∂α∂t)∂α∂t)=(∇∂α∂rR)(D2α∂t2,∂α∂t)∂α∂t+R(D2∂t2∂α∂r,∂α∂t)∂α∂t+R(R(∂α∂r,∂α∂t)∂α∂t,∂α∂t)∂α∂t+R(D2α∂t2,D∂t∂α∂r)∂α∂t+R(D2α∂t2,∂α∂t)D∂t∂α∂r, |
and, consequently,
D∂rΓ(α)=D4∂t4∂α∂r+D2∂t2(R(∂α∂r,∂α∂t)∂α∂t)+D∂t(R(∂α∂r,∂α∂t)D2α∂t2)+R(∂α∂r,∂α∂t)D3α∂t3,(∇∂α∂rR)(D2α∂t2,∂α∂t)∂α∂t+R(D2∂t2∂α∂r,∂α∂t)∂α∂t+R(R(∂α∂r,∂α∂t)∂α∂t,∂α∂t)∂α∂t+R(D2α∂t2,D∂t∂α∂r)∂α∂t+R(D2α∂t2,∂α∂t)D∂t∂α∂r. |
The previous identity allows to write (d/dr)J(αr) as follows:
ddrJ(αr)=∫T0(⟨D4∂t4∂α∂r+D2∂t2(R(∂α∂r,∂α∂t)∂α∂t)+D∂t(R(∂α∂r,∂α∂t)D2α∂t2)+R(∂α∂r,∂α∂t)D3α∂t3+(∇∂α∂rR)(D2α∂t2,∂α∂t)∂α∂t+R(D2∂t2∂α∂r,∂α∂t)∂α∂t+R(R(∂α∂r,∂α∂t)∂α∂t,∂α∂t)∂α∂t+R(D2α∂t2,D∂t∂α∂r)∂α∂t+R(D2α∂t2,∂α∂t)D∂t∂α∂r,Γ(αr)⟩)dt. |
If we apply the properties of the curvature tensor, we see that
ddrJ(αr)=∫T0(⟨D4dt4∂α∂r+D2∂t2(R(∂α∂r,∂α∂t)∂α∂t)+D∂t(R(∂α∂r,∂α∂t)D2α∂t2),Γ(αr)⟩+⟨∂α∂r,R(Γ(αr),D3α∂t3)∂α∂t⟩+⟨(∇∂α∂rR)(D2α∂t2,∂α∂t)∂α∂t,Γ(αr)⟩+⟨R(Γ(αr),∂α∂t)∂α∂t,D2∂t2∂α∂r⟩+⟨R(Γ(αr),∂α∂t)∂α∂t,R(∂α∂r,∂α∂t)∂α∂t⟩−2⟨R(Γ(αr),∂α∂t)D2α∂t2,D∂t∂α∂r⟩+⟨R(Γ(αr),D2α∂t2)∂α∂t,D∂t∂α∂r⟩)dt. |
Now, integrating by parts, one has
ddrJ(αr)=l∑i=1(⟨D3∂t3∂α∂r,Γ(αr)⟩−⟨D2∂t2∂α∂r,D∂tΓ(αr)⟩+⟨Ddt∂α∂r,D2∂t2Γ(αr)⟩−⟨∂α∂r,D3∂t3Γ(αr)⟩+⟨D∂t(R(∂α∂r,∂α∂t)∂α∂t),Γ(αr)⟩−⟨R(∂α∂r,∂α∂t)∂α∂t,D∂tΓ(αr)⟩ +⟨R(∂α∂r,∂α∂t)D2α∂t2,Γ(αr)⟩+⟨R(Γ(αr),∂α∂t)∂α∂t,D∂t∂α∂r⟩−⟨D∂t(R(Γ(αr),∂α∂t)∂α∂t),∂α∂r⟩−2⟨R(Γ(αr),∂α∂t)D2α∂t2,∂α∂r⟩+⟨R(Γ(αr),D2α∂t2)∂α∂t,∂α∂r⟩)t−it+i−1 |
+∫T0⟨∂α∂r,D4∂t4Γ(αr)+R(D2∂t2Γ(αr),∂α∂t)∂α∂t−R(D∂tΓ(αr),D2α∂t2)∂α∂t+R(Γ(αr),D3α∂t3)∂α∂t+(∇Γ(αr)R)(D2α∂t2,∂α∂t)∂α∂t−(∇∂α∂tR)(Γ(αr),∂α∂t)D2α∂t2+(∇∂α∂tR)(Γ(αr),D2α∂t2)∂α∂t+D2∂t2(R(Γ(αr),∂α∂t)∂α∂t)+R(R(Γ(αr),∂α∂t)∂α∂t,∂α∂t)∂α∂t+2D∂t(R(Γ(αr),∂α∂t)D2α∂t2)−D∂t(R(Γ(αr),D2α∂t2)∂α∂t)⟩dt. |
When we take r=0, most of the terms above vanish due to the conditions (2.1) and the fact that both x and X are C3-piecewise smooth on [0,T]. Hence, the expression of δJ(x)(X) can be simplified as follows:
δJ(x)(X)=l−1∑i=1(−⟨(D3Xdt3+Ddt(R(X,V)V)+2R(DVdt,V)X)t=ti,Γ(x(t+i))−Γ(x(t−i))⟩+⟨(D2Xdt2+R(X,V)V)t=ti,DΓdt(x(t+i)−DΓdt(x(t−i)⟩−⟨DXdt(ti),(D2Γ(x)dt2+R(Γ(x),V)V)t+i−(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i⟩+⟨X(ti),Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t+i−Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i⟩+∫T0⟨X,D2dt2(D2∂t2Γ(x)+R(Γ(x),V)V)+R(D2dt2Γ(x)+R(Γ(x),V)V,V)V+Ddt(R(Γ(x),V)DVdt)+(∇Γ(x)R)(DVdt,V)V−R(DdtΓ(x),DVdt)V+R(Γ(x),V)D2Vdt2+R(DVdt,V)DdtΓ(x)⟩dt. |
Lemma 3. If α is an admissible variation of x∈Ω and X the variational vector field associated with α, then
δJ(x)(X)=l−1∑i=1(−⟨(D3Xdt3+Ddt(R(X,V)V)+2R(DVdt,V)X)t=ti,Γ(x(t+i))−Γ(x(t−i))⟩+⟨(D2Xdt2+R(X,V)V)t=ti,DΓdt(x(t+i)−DΓdt(x(t−i)⟩−⟨DXdt(ti),(D2Γ(x)dt2+R(Γ(x),V)V)t+i−(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i⟩+⟨X(ti),Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t+i−Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i⟩)+∫T0⟨X(t),Π(x(t))⟩dt, | (3.24) |
where
Π(x)=D2dt2(D2∂t2Γ(x)+R(Γ(x),V)V)+R(D2dt2Γ(x)+R(Γ(x),V)V,V)V+Ddt(R(Γ(x),V)DVdt)+(∇Γ(x)R)(DVdt,V)V−R(DdtΓ(x),DVdt)V+R(Γ(x),V)D2Vdt2+R(DVdt,V)DdtΓ(x). | (3.25) |
Proposition 4. A necessary condition for the curve x∈Ω to be a local minimizer for the functional J over Ω is that x is smooth and satisfies
Π(x)=0, | (3.26) |
where Π(x) is given by (3.25).
Proof. Suppose that x is a local minimizer of J over Ω. Then, by applying the Hamilton variational principle to J at x, we have δJ(x)(X)=0, ∀X∈TxΩ (see (3.21)), and then the right-hand side of (3.24) is zero, for arbitrary X∈TxΩ.
Let us first consider X∈TxΩ defined by
X(t)=f(t)Π(x(t)),t∈[0,T], |
where Π(x) is the vector field along x given by (3.25) and f is a smooth real-valued function on [0,T] verifying f(k)(ti)=0, k=0,1,2,3, and f(t)>0, for t∈(ti−1,ti), i=1,…,l−1. It follows from (3.24) and (3.21) that Π(x(t))=0 for t∈[0,T].
To show that x is smooth, let us take the vector field X∈TxΩ so that
X(ti)=Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t+i−Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i,DXdt(ti)=(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i−(D2Γ(x)dt2+R(Γ(x),V)V)t=t+i,(D2Xdt2+R(X,V)V)t=ti=DΓdt(x(t+i))−DΓdt(x(t−i)),(D3Xdt3+Ddt(R(X,V)V)+2R(DVdt,V)X)t=ti=Γ(x(t−i))−Γ(x(t+i)), |
for i=1,…,l−1. Thus, since x satisfies the Eq (3.26), it follows from (3.24) and (3.21) that
0=δJ(x)(X)=l−1∑i=1(‖Γ(x(t+i))−Γ(x(t−i))‖2+‖DΓdt(x(t+i)−DΓdt(x(t−i)‖2+‖(D2Γ(x)dt2+R(Γ(x),V)V)t+i−(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i‖2+‖Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t+i−Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i‖2), |
which implies that
Γ(x(t−i))=Γ(x(t+i)),DΓdt(x(t+i))=DΓdt(x(t−i)),(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i=(D2Γ(x)dt2+R(Γ(x),V)V)t=t+i,Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t+i=Ddt(D2Γ(x)dt2+R(Γ(x),V)V)t=t−i, |
for i=1,…,l−1. Hence, Γ(x) is C3 on [0,T] and, consequently x is C7. Since x is a solution of the differential equation of eighth-order (3.26), it follows that the eighth-order covariant derivative of x can be written in terms of the covariant derivatives of order less than 8, so, in turn, the same argument can be used for covariant derivatives of any order. This leads us to conclude that x is smooth on [0,T].
Definition 5. A curve x in M is said to be a Riemannian heptic if x is a smooth solution of the differential Eq (3.26).
Let G be a connected and compact Lie group endowed with a bi-invariant Riemannian metric ⟨⋅,⋅⟩ and let g be its Lie algebra. We denote by Lx:G→G the left-translation map by an element x of G. The left-invariant vector field associated with an element u of g is uL:x↦TeLxu, where TeLx is the differential of Lx at e and e the identity of G. The Levi-Civita connection ∇ induced by ⟨⋅,⋅⟩ is completely determined by its restriction g∇ to g via left-translations. The bilinear map g∇:g2→g is given by
g∇uz=12[u,z], | (4.1) |
for each u and z∈g, and we have
∇uLzL=(g∇uz)L. | (4.2) |
When we restrict the curvature tensor R of G to the Lie algebra g of G, we obtain the map R:g2→g defined by
R(u,z)w=−14[[u,z],w]. | (4.3) |
The differential of the curvature tensor ∇R is constant and equal to zero.
In addition, the inner product on g corresponding to the Riemannian metric verifies
⟨[u,z],w⟩=⟨u,[z,w]⟩. | (4.4) |
The property (4.4) can be described by means of the adjoint representation ad of g, that is, by saying that the map adz is skew-adjoint, for z∈g.
We denote by ♭ and ♯ the musical isomorphisms between g and its dual g∗ defined by the inner product on g. Thus, we know that
(ad∗uz♭)♯=−aduz,u,z∈g. |
More details about bi-invariant metrics on Lie groups can be found in Milnor [27] and Kobayashi and Nomizu [30].
Let x be a curve in G and V the velocity vector field of x. The lift x1 of the curve x to TG can be reduced to the Lie algebra g by considering the curve v given by V=TeLxv. The reduced curve v is described, with respect to a given basis B={e1,…,en} of g, by v=∑ni=1viei, where v1,…,vn can be seen as pseudo-velocities of x with respect to the basis B. The j-order derivative of the curve v is given by v(j)=∑ni=1v(j)iei, j≥1 (and is independent of the chosen basis).
If x is a Riemannian heptic on G then we call Lie hexic to the corresponding reduced curve v in g. To give a characterization of Lie hexics on a Lie algebra we need to compute the expression of TxLx−1Π(x).
Higher-order covariant derivatives of the curve x can be written in terms of the curve v and its derivatives, as we can see for the first two covariant derivatives of the velocity vector field V below.
DVdt=TeLxv′,D2Vdt2=TeLx(v(2)+12[v,v′]). |
Similar formulas can be obtained for the higher-order covariant derivatives of a vector field Y along x. When we consider the reduced curve y of Y in g, it follows that
D2Ydt2+R(Y,V)V=TeLx(y(2)+12[v′,y]+[v,y′]). |
By taking Y=(D/dt)V, we obtain an expression for Γ(x),
Γ(x)=TeLx(v(3)+[v,v(2)]). |
This gives rise to the reduced equation of Riemannian cubics,
v(3)=[v(2),v], | (4.5) |
and leads to the definition of Lie quadratics proposed by Noakes in [31].
We can also derive
Σ(x)=D2dt2Γ(x)+R(Γ(x),V)VandΔ(x)=D2dt2Σ(x)+R(Σ(x),V)V, |
in terms of the curve v and its derivatives and, after many technical entedious computations, it is possible to obtain the reduction of the Eq (3.26) to the Lie algebra g.
Since the Lagrangian function L of the variational problem (P) is left-invariant, this reduction gives rise to the Euler-Poincaré equation for the reduced Lagrangian function l. Therefore, we chose to obtain the reduced equations by applying a variational principle to the Lagrangian function l as we are showing in the next subsection.
The Lagrangian function L associated with the Riemannian heptic problem (P) on the Lie group G is invariant with respect to the lift of the left translation. Under these conditions, we can define the reduced Lagrangian l:g4→IR given by the restriction of L to the Lie algebra g4. If we apply a variational principle to the Lagrangian function l we obtain the Euler-Poincaré equations for the problem (P).
Let x be a curve in G. The corresponding reduced curve v in the Lie algebra g can be lifted to g4 by considering (v;v(1);v(2);v(3)). If we denote the partial functional derivative of l with respect to z∈g by (δl)/(δz), then we have that the curve x in G is a solution of the Euler-Lagrange equations (3.26) iff the reduced curve v in g is a solution of the Euler-Poincaré equations
(∂∂t−ad∗v)(δlδv−∂∂tδlδ˙v+∂2∂t2δlδv(2)−∂3∂t3δlδv(3))=0. |
The restriction of L to g4 is defined as follows:
l:g4→IR;(v,˙v,v(2),v(3))↦12‖v(3)+[v,v(2)]‖2. |
To obtain the Euler-Poincaré equations, we begin by deriving the partial functional derivatives
δlδv=[v(2),v(3)]♭−ad∗v(2)[v,v(2)]♭,δlδ˙v=0,δlδv(2)=[v(3),v]♭−ad∗v[v,v(2)]♭,δlδv(3)=v(3)♭+[v,v(2)]♭. |
Then, it follows that
δlδv−∂∂tδlδ˙v+∂2∂t2δlδv(2)−∂3∂t3δlδv(3)=−ad∗v(2)[v,v(2)]♭+2[v(5),v]♭+5[v(4),˙v]♭+2[v(3),v(2)]♭−ad∗v(2)[v,v(2)]♭−2ad∗˙v[˙v,v(2)]♭−2ad∗˙v[v,v(3)]♭−2ad∗v[˙v,v(3)]♭−ad∗v[v,v(4)]♭−v(6)♭. | (4.6) |
If we apply the dual of the adjoint representation at v to the expression given in (4.6), we conclude that
ad∗v(δlδv−∂∂tδlδ˙v+∂2∂t2δlδv(2)−∂3∂t3δlδv(3))=2ad∗v[v(5),v]♭+5ad∗v[v(4),˙v]♭+2ad∗v[v(3),v(2)]♭−ad∗vv(6)♭−ad∗vad∗v(2)[v,v(2)]♭−ad∗vad∗v(2)[v,v(2)]♭−2ad∗vad∗˙v[˙v,v(2)]♭−2ad∗vad∗˙v[v,v(3)]♭−2ad∗vad∗v[˙v,v(3)]♭−ad∗vad∗v[v,v(4)]♭. | (4.7) |
On the other hand, deriving (4.6) with respect to time, we obtain
∂∂t(δlδv−∂∂tδlδ˙v+∂2∂t2δlδv(2)−∂3∂t3δlδv(3))=−2ad∗v(3)[v,v(2)]♭−4ad∗v(2)[˙v,v(2)]♭−4ad∗v(2)[v,v(3)]♭+2[v(6),v]♭+7[v(5),˙v]♭+7[v(4),v(2)]♭−6ad∗˙v[˙v,v(3)]♭−3ad∗˙v[v,v(4)]♭−2ad∗v[v(2),v(3)]♭−3ad∗v[˙v,v(4)]♭−ad∗v[v,v(5)]♭−v(7)♭. | (4.8) |
Finally, combining Eqs (4.7) and (4.8) and applying the isomorphism ♯, we get the Euler-Poincaré equation for l as follows:
v(7)=3[v(6),v]+7[v(5),˙v]+7[v(4),v(2)]+4[v(3),[v,v(2)]]+4[v(2),[˙v,v(2)]]+2[v(2),[v,v(3)]]+6[˙v,[˙v,v(3)]]+3[˙v,[v,v(4)]]+[v,[v(5),v]]+5[v,[v(4),˙v]]+2[[[v,v(2)],v(2)],v]+2[[[˙v,v(2)],˙v],v]+2[[[v,v(3)],˙v],v]+2[[[˙v,v(3)],v],v]+[[[v,v(4)],v],v]. | (4.9) |
Therefore, Lie hexics are the smooth curves on the Lie algebra g that are solutions of the Euler-Poincaré equation (4.9). Consequently, solving Eq (4.9) together with equation dx/dt=TeLxv leads to Riemannian heptics on G.
In this paper, we study a fourth-order extension of the geodesic problem that appears as the most natural extension of the Riemannian cubic problem. We derive the Euler-Lagrange equation for this variational problem and obtain the concept of Riemannin heptics as a seventh-degree generalization of Riemannian cubics. These curves may, in fact, be seen as the best approximations to Riemannian cubics when the existence of cubics is conditioned by boundary or interpolation conditions, and reduce naturally to Riemannian cubics when such curves exist. From this point of view, Riemannian heptics are a more conceptual generalization of Euclidean heptic polynomials to Riemannian manifolds than the geometrical polynomials defined in [18,19,20].
When we specialized the problem to Lie groups, we obtained an invariant variational problem, which allowed us to derive the corresponding Euler-Poincaré equation (4.9). Hence, this equation characterizes the Lie hexics and, consequently, the Riemannian heptics on Lie groups.
Furthermore, this study shows the advantages of the Euler-Poincaré reduction. This method allows us to significantly simplify the calculations that we would have to do to reduce the Euler-Lagrange equation to the Lie algebra. In fact, it was possible to avoid deriving all the terms with four Lie brackets that would appear if we applied the reduction directly to the Euler-Lagrange equation (3.26) (see Subsection 4.2).
The authors declare that they have not used artificial intelligence (AI) tools in the creation of this article.
This work was partially supported by the Centre for Mathematics of the University of Coimbra (funded by the Portuguese Government through FCT/MCTES, DOI 10.54499/UIDB/00324/2020).
The authors declare that there are no conflicts of interest.
[1] |
L. Noakes, G. Heinzinger, B. Paden, Cubic splines on curved spaces, IMA J. Math. Control Inf., 6 (1989), 465–473. https://doi.org/10.1093/imamci/6.4.465 doi: 10.1093/imamci/6.4.465
![]() |
[2] |
R. Caddeo, S. Montaldo, C. Oniciuc, P. Piu, The Euler-Lagrange method for biharmonic curves, Mediterr. J. Math., 3 (2006), 449–465. https://doi.org/10.1007/s00009-006-0090-x doi: 10.1007/s00009-006-0090-x
![]() |
[3] |
L. Machado, F. S. Leite, K. Krakowski, Higher-order smoothing splines versus least squares problems on Riemannian manifolds, J. Dyn. Control Syst., 16 (2010), 121–148. https://doi.org/10.1007/s10883-010-9080-1 doi: 10.1007/s10883-010-9080-1
![]() |
[4] |
K. Hüper, F. S. Leite, On the geometry of rolling and interpolation curves on Sn, SOn, and Grassmann manifolds, J. Dyn. Control Syst., 13 (2007), 467–502. https://doi.org/10.1007/s10883-007-9027-3 doi: 10.1007/s10883-007-9027-3
![]() |
[5] |
E. Batzies, K. Hüper, L. Machado, F. S. Leite, Geometric mean and geodesic regression on Grassmannians, Linear Algebra Appl., 466 (2015), 83–101. https://doi.org/10.1016/j.laa.2014.10.003 doi: 10.1016/j.laa.2014.10.003
![]() |
[6] |
P. Crouch, F. S. Leite, The dynamic interpolation problem on Riemannian manifolds, Lie groups and symmetric spaces, J. Dyn. Control Syst., 1 (1995), 177–202. https://doi.org/10.1007/BF02254638 doi: 10.1007/BF02254638
![]() |
[7] |
L. Abrunheiro, M. Camarinha, J. Clemente-Gallardo, Cubic polynomials on Lie groups: Reduction of the Hamiltonian system, J. Phys. A: Math. Theor., 44 (2011), 355203. https://doi.org/10.1088/1751-8113/44/35/355203 doi: 10.1088/1751-8113/44/35/355203
![]() |
[8] |
L. Abrunheiro, M. Camarinha, Optimal control of affine connection control systems from the point of view of Lie algebroids, Int. J. Geom. Methods Mod. Phys., 11 (2014), 1450038. https://doi.org/10.1142/S0219887814500388 doi: 10.1142/S0219887814500388
![]() |
[9] |
C. Altafini, Reduction by group symmetry of second order variational problems on a semidirect product of Lie groups with positive definite Riemannian metric, ESAIM Control Optim. Calc. Var., 10 (2004), 526–548. https://doi.org/10.1051/cocv:2004018 doi: 10.1051/cocv:2004018
![]() |
[10] |
P. Balseiro, T. J. Stuchi, A. Cabrera, J. Koiller, About simple variational splines from the Hamiltonian viewpoint, J. Geom. Mech., 9 (2017), 257–290. https://doi.org/10.3934/jgm.2017011 doi: 10.3934/jgm.2017011
![]() |
[11] | M. Barbero-Liñán, M. C. Muñoz-Lecanda, Optimal control problems for affine connection control systems: Characterization of extremals, in A: AIP Conference Proceedings. Geometry and physics: XVI International Fall Workshop, American Institute of Physics, 1023 (2008), 127–131. https://doi.org/10.1063/1.2958162 |
[12] | P. Crouch, F. S. Leite, M. Camarinha, Hamiltonian structure of generalized cubic polynomials, in Proceedings IFAC Workshop on Lagrangian and Hamiltonian Methods for Nonlinear Control, Princeton University, USA, 33 (2000), 13–18. https://doi.org/10.1016/S1474-6670(17)35541-6 |
[13] | A. Bloch, L. Colombo, R. Gupta, D. M. de Diego, A geometric approach to the optimal control of nonholonomic mechanical systems, in Analysis and Geometry in Control Theory and its Applications (eds. P. Bettiol, P. Cannarsa, G. Colombo, M. Motta and F. Rampazzo), Springer, 11 (2015), 35–64. https://doi.org/10.1007/978-3-319-06917-3_2 |
[14] |
I. I. Hussein, A. M. Bloch, Dynamic coverage optimal control for multiple spacecraft interferometric imaging, J. Dyn. Control Syst., 13 (2007), 69–93. https://doi.org/10.1007/s10883-006-9004-2 doi: 10.1007/s10883-006-9004-2
![]() |
[15] |
A. Bloch, M. Camarinha, L. J. Colombo, Variational point-obstacle avoidance on Riemannian manifolds, Math. Control Signals Syst., 33 (2021), 109–121. https://doi.org/10.1007/s00498-021-00276-0 doi: 10.1007/s00498-021-00276-0
![]() |
[16] |
J. Goodman, L. Colombo, Collision avoidance of multiagent systems on Riemannian manifolds, SIAM J. Control Optim., 60 (2022), 168–188. https://doi.org/10.1137/21M1411056 doi: 10.1137/21M1411056
![]() |
[17] |
E. Stratoglou, A. A. Simoes, L. J. Colombo, Reduction in optimal control with broken symmetry for collision and obstacle avoidance of multi-agent system on Lie groups, Commun. Anal. Mech., 15 (2023), 1–23. http://doi.org/10.3934/cam.2023001 doi: 10.3934/cam.2023001
![]() |
[18] |
M. Camarinha, P. Crouch, F. S. Leite, Splines of class Ck on non-euclidean spaces, IMA J. Math. Control Inf., 12 (1995), 399–410. https://doi.org/10.1093/imamci/12.4.399 doi: 10.1093/imamci/12.4.399
![]() |
[19] | M. Camarinha, F. S. Leite, P. Crouch, Geometrical polynomials of odd degree in Riemannian manifolds, in 1997 European Control Conference (ECC), Brussels, Belgium, (1997), 1872–1878. https://doi.org/10.23919/ECC.1997.7082377 |
[20] |
M. Camarinha, F. S. Leite, P. Crouch, High-order splines on Riemannian manifolds, Proc. Steklov Inst. Math., 321 (2023), 158–178. https://doi.org/10.1134/S0081543823020128 doi: 10.1134/S0081543823020128
![]() |
[21] |
V. Branding, A structure theorem for polyharmonic maps between Riemannian manifolds, J. Differ. Equations, 273 (2021), 14–39. https://doi.org/10.1016/j.jde.2020.11.046 doi: 10.1016/j.jde.2020.11.046
![]() |
[22] | M. Camarinha, Luís Machado, F. S. Leite, k-Splines on SPD manifolds, in International Conference on Geometric Science of Information, Springer Nature Switzerland, 14072 (2023), 624–633. https://doi.org/10.1007/978-3-031-38299-4_64 |
[23] |
J. Hinkle, P. T. Fletcher, S. Joshi, Intrinsic polynomials for regression on Riemannian manifolds, J. Math. Imaging Vision, 50 (2014), 32–52. https://doi.org/10.1007/s10851-013-0489-5 doi: 10.1007/s10851-013-0489-5
![]() |
[24] |
C. Samir, P. A. Absil, A. Srivastava, E. Klassen, A gradient-descent method for curve fitting on Riemannian manifolds, Found. Comput. Math., 12 (2012), 49–73. https://doi.org/10.1007/s10208-011-9091-7 doi: 10.1007/s10208-011-9091-7
![]() |
[25] | M. Camarinha, A 4th-Order Variational Problem on SO(3), in APCA International Conference on Automatic Control and Soft Computing, Caparica, Portugal, 930 (2022), 353–362. https://doi:10.1007/978-3-031-10047-5_31 |
[26] |
F. Gay-Balmaz, D. D. Holm, T. S. Ratiu, Higher-order Lagrange-Poincaré and Hamilton-Poincaré reductions, Bull. Braz. Math. Soc. New Ser., 42 (2011), 579–606. https://doi.org/10.1007/s00574-011-0030-7 doi: 10.1007/s00574-011-0030-7
![]() |
[27] | J. Milnor, Morse Theory, Princeton University, 1969. |
[28] | L. M. Lee, Riemannian Manifolds: An Introduction to Curvature, Springer-Verlag, New York, 1997. |
[29] | B. O'Neill, Semi-Riemannian Geometry: With Applications to Relativity, Academic Press, New York, 1983. |
[30] | S. Kobayashi, K. Nomizu, Foundations of Differential Geometry, Interscience Publishers, New York, 1963. |
[31] |
L. Noakes, T. Popiel, Quadratures and cubics in SO(3) and SO(1, 2), IMA J. Math. Control Inf., 23 (2006), 463–473. https://doi.org/10.1093/imamci/dni069 doi: 10.1093/imamci/dni069
![]() |