Using comparison principles, we examine the asymptotic characteristics of a third-order nonlinear damped neutral differential equation. Our results substantially generalize numerous previously established results as well as drastically improving them. To illustrate the relevance and effectiveness of our results, we use numerical examples.
Citation: Taher S. Hassan, Emad R. Attia, Bassant M. El-Matary. Iterative oscillation criteria of third-order nonlinear damped neutral differential equations[J]. AIMS Mathematics, 2024, 9(8): 23128-23141. doi: 10.3934/math.20241124
[1] | Yuriĭ G. Nikonorov, Irina A. Zubareva . On the behavior of geodesics of left-invariant sub-Riemannian metrics on the group Aff0(R)×Aff0(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 |
Using comparison principles, we examine the asymptotic characteristics of a third-order nonlinear damped neutral differential equation. Our results substantially generalize numerous previously established results as well as drastically improving them. To illustrate the relevance and effectiveness of our results, we use numerical examples.
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] |
R. Agarwal, M. Bohner, T. Li, C. Zhang, Oscillation of third-order nonlinear delay differential equations, Taiwanese J. Math., 17 (2013), 545–558. http://dx.doi.org/10.11650/tjm.17.2013.2095 doi: 10.11650/tjm.17.2013.2095
![]() |
[2] | R. Agarwal, S. Grace, D. O'Regan, Oscillation theory for difference and functional differential equations, Dordrecht: Springer, 2000. http://dx.doi.org/10.1007/978-94-015-9401-1 |
[3] |
R. Agarwal, S. Grace, D. O'Regan, On the oscillation of certain functional differential equations via comparison methods, J. Math. Anal. Appl., 286 (2003), 577–600. http://dx.doi.org/10.1016/S0022-247X(03)00494-3 doi: 10.1016/S0022-247X(03)00494-3
![]() |
[4] |
R. Agarwal, S. Grace, D. O'Regan, The oscillation of certain higher order functional differential equations, Math. Comput. Model., 37 (2003), 705–728. http://dx.doi.org/10.1016/S0895-7177(03)00079-7 doi: 10.1016/S0895-7177(03)00079-7
![]() |
[5] | B. Bacuíková, J. Džurina, Oscillation of third-order functional differential equations, Electron. J. Qual. Theo., 43 (2010), 1–10. |
[6] |
B. Bacuíková, J. Džurina, Oscillation of third-order nonlinear differential equations, Appl. Math. Lett., 24 (2011), 466–470. http://dx.doi.org/10.1016/j.aml.2010.10.043 doi: 10.1016/j.aml.2010.10.043
![]() |
[7] |
B. Bacuíková, E. Elabbasy, S. Saker, J. Džurina, Oscillation criteria for third-order nonlinear differential equations, Math. Slovaca, 58 (2008), 201–220. http://dx.doi.org/10.2478/s12175-008-0068-1 doi: 10.2478/s12175-008-0068-1
![]() |
[8] |
M. Bohner, S. Grace, I. Sağer, E. Tunç, Oscillation of third-order nonlinear damped delay differential equations, Appl. Math. Comput., 278 (2016), 21–32. http://dx.doi.org/10.1016/j.amc.2015.12.036 doi: 10.1016/j.amc.2015.12.036
![]() |
[9] |
G. Chatzarakis, J. Džurina, I. Jadlovska, Oscillatory and asymptotic properties of third-order quasilinear delay differential equations, J. Inequal. Appl., 2019 (2019), 23. http://dx.doi.org/10.1186/s13660-019-1967-0 doi: 10.1186/s13660-019-1967-0
![]() |
[10] |
G. Chatzarakis, S. Grace, I. Jadlovská, Oscillation criteria for third-order delay differential equations, Adv. Differ. Equ., 2017 (2017), 330. http://dx.doi.org/10.1186/s13662-017-1384-y doi: 10.1186/s13662-017-1384-y
![]() |
[11] |
E. Elabbasy, B. Qaraad, T. Abdeljawad, O. Moaaz, Oscillation criteria for a class of third-order damped neutral differential equations, Symmetry, 12 (2020), 1988. http://dx.doi.org/10.3390/sym12121988 doi: 10.3390/sym12121988
![]() |
[12] | E. Elabbasy, T. Hassan, B. Elmatary, Oscillation criteria for third order delay nonlinear differential equations, Electron. J. Qual. Theo., 5 (2012), 1–11. |
[13] |
L. Erbe, T. Hassan, A. Peterson, Oscillation of third order nonlinear functional dynamic equations on time scales, Differ. Equ. Dyn. Syst., 18 (2010), 199–227. http://dx.doi.org/10.1007/s12591-010-0005-y doi: 10.1007/s12591-010-0005-y
![]() |
[14] | L. Erbe, Q. Kong, B. Zhang, Oscillation theory for functional differential equations, New York: Routledge, 1995. http://dx.doi.org/10.1201/9780203744727 |
[15] |
S. Grace, Oscillation criteria for third order nonlinear delay differential equations with damping, Opusc. Math., 35 (2015), 485–497. http://dx.doi.org/10.7494/OpMath.2015.35.4.485 doi: 10.7494/OpMath.2015.35.4.485
![]() |
[16] |
S. Grace, R. Agarwal, R. Pavani, E. Thandapani, On the oscilation certain third order nonlinear functional differential equations, Appl. Math. Comput., 202 (2008), 102–112. http://dx.doi.org/10.1016/j.amc.2008.01.025 doi: 10.1016/j.amc.2008.01.025
![]() |
[17] | K. Gopalsamy, Stability and oscillation in delay differential equations of population dynamics, Dordrecht: Springer, 1992. http://dx.doi.org/10.1007/978-94-015-7920-9 |
[18] |
I. Gyori, F. Hartung, Stability of a single neuron model with delay, J. Comput. Appl. Math., 157 (2003), 73–92. http://dx.doi.org/10.1016/S0377-0427(03)00376-5 doi: 10.1016/S0377-0427(03)00376-5
![]() |
[19] | I. Gyori, G. Ladas, Oscillation theory of delay differential equations with applications, Oxford: Clarendon Press, 1991. http://dx.doi.org/10.1093/oso/9780198535829.001.0001 |
[20] | J. Hale, S. Verduyn Lunel, Introduction to functional differential equations, New York: Springer-Verlag, 1993. http://dx.doi.org/10.1007/978-1-4612-4342-7 |
[21] |
T. Hassan, B. El-Matary, Asymptotic behavior and oscillation of third-order nonlinear neutral differential equations with mixed nonlinearities, Mathematics, 11 (2023), 424. http://dx.doi.org/10.3390/math11020424 doi: 10.3390/math11020424
![]() |
[22] | T. Hassan, B. El-Matary, Oscillation criteria for third order nonlinear neutral differential equation, PLOMS Math., 1 (2021), 00001. |
[23] |
T. Hassan, L. Erbe, A. Peterson, Forced oscillation of second order functional differential equations with mixed nonlinearities, Acta Math. Sci., 31 (2011), 613–626. http://dx.doi.org/10.1016/S0252-9602(11)60261-0 doi: 10.1016/S0252-9602(11)60261-0
![]() |
[24] | T. Hassan, Q. Kong, Interval criteria for forced oscillation of differential equations with p-Laplacian, damping, and mixed nonlinearities, Dynam. Syst. Appl., 20 (2011), 279–294. |
[25] |
I. Jadlovská, G. Chatzarakis, J. Džurina, S. Grace, On sharp oscillation criteria for general third-order delay differential equations, Mathematics, 9 (2021), 1675. http://dx.doi.org/10.3390/math9141675 doi: 10.3390/math9141675
![]() |
[26] | Y. Kitamura, Oscillation of functional differential equations with general deviating arguments, Hiroshima Math. J., 15 (1985), 445–491. |
[27] | R. Koplatadze, T. Chanturiya, Oscillating and monotone solutions of first-order differential equations with deviating argument, Differ. Uravn., 18 (1982), 1463–1465. |
[28] |
B. Karpuz, Ö. Öcalan, New oscillation tests and some refinements for first-order delay dynamic equations, Turk. J. Math., 40 (2016), 850–863. http://dx.doi.org/10.3906/mat-1507-98 doi: 10.3906/mat-1507-98
![]() |
[29] |
G. Ladas, Y. Sficas, I. Stavroulakis, Necessary and sufficient conditions for oscillations of higher order delay differential equations, Trans. Amer. Math. Soc., 285 (1984), 81–90. http://dx.doi.org/10.2307/1999473 doi: 10.2307/1999473
![]() |
[30] | G. Ladas, V. Lakshmikantham, J. Papadakis, Oscillations of higher-order retarded differential equations generated by the retarded argument, In: Delay and functional differential equations and their applications, New York: Academic Press, 1972,219–231. http://dx.doi.org/10.1016/B978-0-12-627250-5.50013-7 |
[31] |
G. Ladas, V. Lakshmikantham, Sharp conditions for oscillations caused by delays, Appl. Anal., 9 (1979), 93–98. http://dx.doi.org/10.1080/00036817908839256 doi: 10.1080/00036817908839256
![]() |
[32] |
T. Li, Y. Rogovchenko, On the asymptotic behavior of solutions to a class of third-order nonlinear neutral differential equations, Appl. Math. Lett., 105 (2020), 106293. http://dx.doi.org/10.1016/j.aml.2020.106293 doi: 10.1016/j.aml.2020.106293
![]() |
[33] |
T. Li, C. Zhang, G. Xing, Oscillation of third-order neutral delay differential equations, Abstr. Appl. Anal., 2012 (2012), 569201. http://dx.doi.org/10.1155/2012/569201 doi: 10.1155/2012/569201
![]() |
[34] |
O. Moaaz, I. Dassios, W. Muhsin, A. Muhib, Oscillation theory for non-linear neutral delay differential equations of third order, Appl. Sci., 10 (2020), 4855. http://dx.doi.org/10.3390/app10144855 doi: 10.3390/app10144855
![]() |
[35] |
O. Moaaz, E. Elabbasy, E. Shaaban, Oscillation criteria for a class of third order damped differential equations, Arab Journal of Mathematical Sciences, 24 (2018), 16–30. http://dx.doi.org/10.1016/j.ajmsc.2017.07.001 doi: 10.1016/j.ajmsc.2017.07.001
![]() |
[36] | S. Padhi, S. Pati, Theory of third-order differential equations, New Delhi: Springer, 2014. http://dx.doi.org/10.1007/978-81-322-1614-8 |
[37] |
Ch. Philos, On the existence of nonoscillatory solutions tending to zero at ∞ to differential equations with positive delays, Arch. Math., 36 (1981), 168–178. http://dx.doi.org/10.1007/BF01223686 doi: 10.1007/BF01223686
![]() |
[38] |
A. Tiryaki, M. Aktas, Oscillation criteria of a certain class of third order nonlinear delay differential equations with damping, J. Math. Anal. Appl., 325 (2007), 54–68. http://dx.doi.org/10.1016/j.jmaa.2006.01.001 doi: 10.1016/j.jmaa.2006.01.001
![]() |
[39] |
M. Wei, C. Jiang, T. Li, Oscillation of third-order neutral differential equations with damping and distributed delay, Adv. Differ. Equ., 2019 (2019), 426. http://dx.doi.org/10.1186/s13662-019-2363-2 doi: 10.1186/s13662-019-2363-2
![]() |
[40] |
L. Yang, Z. Xu, Oscillation of certain third-order quasilinear neutral differential equations, Math. Slovaca, 64 (2014), 85–100. http://dx.doi.org/10.2478/s12175-013-0189-z doi: 10.2478/s12175-013-0189-z
![]() |