Next Article in Journal
Monitoring and Recognizing Enterprise Public Opinion from High-Risk Users Based on User Portrait and Random Forest Algorithm
Next Article in Special Issue
On Uniform Stability with Growth Rates of Stochastic Skew-Evolution Semiflows in Banach Spaces
Previous Article in Journal
An Intelligent System for Allocating Times to the Main Activities of Managers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Smooth Stable Manifold for Delay Differential Equations with Arbitrary Growth Rate

1
Department of Mathematics, University of Rijeka, HR-51000 Rijeka, Croatia
2
Department of Mathematics and Statistics, Indian Institute of Technology Kanpur, Kanpur 208016, India
*
Author to whom correspondence should be addressed.
Axioms 2021, 10(2), 105; https://doi.org/10.3390/axioms10020105
Submission received: 19 April 2021 / Revised: 11 May 2021 / Accepted: 20 May 2021 / Published: 25 May 2021
(This article belongs to the Special Issue Nonautonomous and Random Dynamical Systems)

Abstract

:
In this article, we construct a C 1 stable invariant manifold for the delay differential equation x = A x ( t ) + L x t + f ( t , x t ) assuming the ρ -nonuniform exponential dichotomy for the corresponding solution operator. We also assume that the C 1 perturbation, f ( t , x t ) , and its derivative are sufficiently small and satisfy smoothness conditions. To obtain the invariant manifold, we follow the method developed by Lyapunov and Perron. We also show the dependence of invariant manifold on the perturbation f ( t , x t ) .

1. Introduction

The invariant manifold theory started with the work of Hadamard [1] in 1901 when he constructed a manifold in the solution space of a differential equation with the property that if the trajectory of a solution starts in the manifold, it will remain in the manifold for all time t > 0 . This proved to be of great importance for analyzing complex systems as it reduces the relevant dimension significantly. Later, Perron [2] and Lyapunov [3] developed another method to construct the invariant manifolds for autonomous differential equations. In the 1970s, due to the fundamental works of Hirsch et al. [4,5], Sacker and Sell [6], and Pesin [7,8], this theory became an important instrument for various fields like applied mathematics, biology, and engineering.
In the first method for constructing the invariant manifolds, Hadamard [1] used the geometrical properties of differential equations. He constructed the manifold over the linearized stable and unstable subspaces. However, Lyapunov [3] and Perron [2] developed an analytical method to construct the invariant manifolds. They obtained the invariant manifolds, using the variation of constants formula of the differential equations. In his approach, Perron introduced (and assumed) the notion of (uniform) exponential dichotomy for the solution operators and proved the existence of Lipschitz stable invariant manifolds for the small nonlinear perturbation of autonomous differential equations. The smoothness of these invariant manifolds is proved by Pesin [7]. In 1977, Pesin [8] generalized the notion of uniform hyperbolicity to nonuniform hyperbolicity which allows the rate of expansion and contraction to depend on initial time. Later he proved the stable manifold theorem in the finite-dimensional settings for nonhyperbolic trajectories. Pugh and Shub [9] proved a similar result for nonhyperbolic trajectories using the method developed by Hadamard. Ruelle [10] extended the result by Pesin to the Hilbert space settings in 1982.
The exponential dichotomy played an essential role in the development of invariant manifold theory for autonomous differential equations. Barreira and Valls extended the notion of exponential dichotomy for nonautonomous differential equations and called it nonuniform exponential dichotomy. With the assumption of nonuniform exponential dichotomy, they constructed the Lipschitz invariant manifold [11] and smooth invariant manifold [12] for nonautonomous differential equations. They also obtained the essential conditions for the existence of the nonuniform exponential dichotomy. The book [13] contains all the early works of Barreira and Valls.
Barreira and Valls observed that the solution operators corresponding to a class of nonautonomous differential equations show dichotomic behavior and also they have growth or decay rates of e c ρ ( t ) , for some function ρ ( t ) . They named this notion as ρ -nonuniform exponential dichotomy. They showed in article [14] that the class of differential equations for which all the Lyapunov’s exponents are infinite, satisfies ρ -nonuniform exponential dichotomy for ρ ( t ) t . Subsequently they proved the stable manifold theorem for ordinary differential equations assuming ρ -nonuniform exponential dichotomy in [15]. In the article [16], Pan proved the existence of Lipschitz stable invariant manifold for impulsive nonautonomous differential equations with the assumption of ρ -nonuniform exponential dichotomy.
In this article, we consider a differential equation with the infinite delay given by,
x = A x ( t ) + L x t + f ( t , x t ) , x s = ϕ ,
in a Banach space X. We assume that the solution operator associated with the corresponding linear delay differential equation
x = A x ( t ) + L x t , x s = ϕ ,
satisfies ρ -nonuniform exponential dichotomy and the nonlinear perturbation f ( t , x t ) is sufficiently small and smooth. With these assumptions, we prove the existence of a C 1 stable invariant manifold for the delay differential Equation (1) following the approach of Perron and Lyapunov. We also showed the dependence of invariant manifolds on perturbations.
Barreira and Valls used the ρ -nonuniform exponential dichotomy for the nonautonomous differential equations in [17] to construct the Lipschitz stable invariant manifold and in [18] to construct the smooth invariant manifold. Pan considered the impulsive differential equation in [16], to construct the Lipschitz stable invariant manifold assuming ρ -nonuniform exponential dichotomy. In the article [19], we considered the case of delay differential equations with nonuniform exponential dichotomy and constructed a Lipschitz invariant manifold. However, in this article, we are assuming a more general ρ -nonuniform exponential dichotomy for differential equation with infinite delay and we are constructing a C 1 stable invariant manifold. In the later part of the article, we also show that a small change in perturbation gives rise to a small variation in the manifold.
The paper is arranged in the following manner. Our setup and some preliminary results are given in Section 2. In the next section, we provide a few examples of differential equations satisfying the ρ -nonuniform exponential dichotomy. Section 3 contains the proof of the existence of the C 1 stable invariant manifold, and in Section 4, we prove the dependency of the manifold on perturbation. In the end, we present a few more examples satisfying the assumptions of our main theorem.

2. Preliminaries

Let ( X , X ) be a Banach space. For any interval J R : = ( , ) , we denote C ( J , X ) as a space of X-valued continuous function on J. For a function x : ( , a ] X and t a , we define a function x t : R : = ( , 0 ] X by x t ( θ ) : = x ( t + θ ) for θ R . Furthermore, let C γ be a space of continuous functions defined by
C γ : = ψ C ( R , X ) : lim θ ψ ( θ ) X e γ θ = 0 ,
for γ > 0 . We define a norm on the phase space C γ ,
ψ C γ : = sup θ R ψ ( θ ) X e γ θ , ψ C γ .
Finally, we consider a linear delay differential equation in the Banach space X,
x = A x ( t ) + L x t , x s = ϕ ,
for ( s , ϕ ) R + × C γ . The linear operator A : X X generates a strongly continuous compact semigroup { T ( t , s ) } t 0 and L : C γ X is a bounded linear operator. Let the evolution operator corresponding to the above differential equation is denoted by V ( t , s ) and for ϕ C γ , V ( t , s ) is given by
[ V ( t , s ) ϕ ] ( θ ) = x t ( θ , s , ϕ , 0 ) , θ R , t s .
One can easily see that the evolution operator defines a strongly continuous semigroup and for every t s r 0 , it satisfies the semigroup property given by:
V ( t , s ) V ( s , r ) = V ( t , r ) , and V ( t , t ) = I .
Here, I is an identity operator on C γ . For a continuous function p : R X , we consider the perturbed system of delay differential equation:
x = A x ( t ) + L x t + p ( t ) , x s = ϕ .
Let the solution of the above delay differential equation is denoted by ( x t ( s , ϕ , p ) . Now, we give a representation of x t ( s , ϕ , p ) depending on the evolution operator { V ( t , s ) } t s . Let us introduce a function Γ n given by
Γ n ( θ ) = ( n θ + 1 ) I X , 1 n θ 0 , 0 , θ < 1 n ,
where n is any positive integer and I X is the identity operator on X. It is easy to verify that for y X ,
Γ n y C γ and Γ n y C γ max { 1 , e γ n } y X y X .
The next result by Hino and Naito [20] establishes the variation of constants formula for the delay differential Equation (4) in the phase space C γ .
Proposition 1.
Let ( s , ϕ ) R + × C γ be given. Then the segment x t ( s , ϕ , p ) of solution x ( · , s , ϕ , p ) of non-homogeneous functional differential Equation (4) satisfies the following relation in C γ :
x t ( s , ϕ , p ) = V ( t , s ) ϕ + lim n s t V ( t , τ ) Γ n p ( τ ) d τ , t s .
Definition 1.
[ρ-nonuniform exponential dichotomy] Let ρ : [ 0 , ) [ 0 , ) be an increasing function with ρ ( t ) as t . We say that the linear Equation (2) admits a ρ-nonuniform exponential dichotomy if for every t s 0 , there exist projection maps P ( t ) : C γ C γ , constants a < 0 b , ϵ 0 and K > 1 , such that :
(i)
P ( t ) V ( t , s ) = V ( t , s ) P ( s ) ;
(ii)
V Q ( t , s ) : = V ( t , s ) : Q ( s ) C γ Q ( t ) C γ is invertible, where Q ( t ) = I P ( t ) is the complementary projection ;
(iii)
V ( t , s ) P ( s ) K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) , V Q ( t , s ) 1 Q ( t ) K e b ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( t ) .
Now, for each t 0 , we define
E ( t ) : = P ( t ) ( C γ ) and F ( t ) : = Q ( t ) ( C γ ) .
We call E ( t ) and F ( t ) stable and unstable subspace respectively. Let us take examples of some differential equations which satisfy ρ -nonuniform exponential dichotomy.
Remark 1.
Here we are presenting a few examples of differential equations that satisfy the ρ-nonuniform exponential dichotomy. In the first example, we consider a case where differential equation satisfies ρ-nonuniform exponential dichotomy for ρ ( t ) = t .
Example 1.
[21]
x ( t ) = α cos t a + sin t + β t sin t x ( t ) ,
for α > 2 β > 0 .
Let V ( t , s ) be the solution operator of the above problem, then
V ( t , s ) = e α ( t s ) ln ( a + sin t ) + ln ( a + sin s ) + β ( s cos s t cos t ) + β ( sin t sin s ) .
V ( t , s ) = e α ( t s ) ln ( a + sin t ) + ln ( a + sin s ) + β ( s cos s t cos t ) + β ( sin t sin s ) e 2 β + 2 ( ln ( a + 1 ) ) e ( α + β ) ( t s ) + 2 β s = C e ( α + β ) ( t s ) + 2 β s .
Note that ( α + β ) < 0 and C is a constant. Hence, the differential equation satisfies ρ-nonuniform exponential dichotomy with ρ ( t ) = t .
In the next example, we consider a system of two differential equations which satisfies ρ -nonuniform exponential dichotomy with ρ ( t ) = t 2 .
Example 2.
x ( t ) = 2 ω t 1 + e t 2 + 2 ϵ t t 2 sin t 2 1 x ( t ) , y ( t ) = 2 ω t 1 + e t 2 2 ϵ t t 2 cos t 2 1 y ( t ) ,
for some ω > 2 ϵ > 0 .
The evolution operator for linear system above is given by
T ( t , s ) = U ( t , s ) 0 0 V ( t , s ) ,
where
U ( t , s ) = e ω ( t 2 s 2 ) ω ln ( 1 + e t 2 ) ( 1 + e s 2 ) + ϵ ( s 2 cos s 2 t 2 cos t 2 ) + ϵ ( sin t 2 sin s 2 ) ϵ ( t 2 s 2 )
and
V ( t , s ) = e ω ( t 2 s 2 ) + ω ln ( 1 + e t 2 ) ( 1 + e s 2 ) + ϵ ( s 2 sin s 2 t 2 sin t 2 ) + ϵ ( cos t 2 cos s 2 ) ϵ ( t 2 s 2 ) .
Now, for the projection map P ( t ) ( x , y ) = ( x , 0 ) ,
T ( t , s ) P ( t ) = U ( t , s ) e 2 ϵ + ω ln 2 e ω ( t 2 s 2 ) + 2 ϵ s 2 ,
and
T ( t , s ) 1 Q ( t ) = V ( t , s ) 1 e 2 ϵ + ω ln 2 e ω ( t 2 s 2 ) + 2 ϵ t 2 ,
for t s 0 . Hence, the linear system satisfies ρ-nonuniform exponential dichotomy with a = ω < 0 and b = ω > 0 where ρ ( t ) = t 2 .

3. Stable Manifold Theorem

This section is dedicated to the construction of the C 1 stable invariant manifold for the delay differential Equation (1). Let Ω be the space of continuous functions f ( t , ϕ ) : R + × C γ X such that
(i)
f ( t , 0 ) = 0 and f ϕ ( t , 0 ) = 0 for all t 0 ;
(ii)
There exists a positive function b : [ 0 , ) ( 0 , ) such that,
j f ϕ j ( t , ϕ ) j f ϕ j ( t , ψ ) X b ( t ) ϕ ψ C γ ,
for j = 0 , 1 .
Proposition 1 ensures the existence of global solution x t ( · , s , ϕ , f ) of the differential Equation (1) with the above mentioned properties of perturbation f ( t , x t ) . We also assume that the solution operator V ( t , s ) satisfies ρ -nonuniform exponential dichotomy. Therefore, Using the projection maps, we can project the global solution x t = ( u t , v t ) E ( t ) × F ( t ) with initial condition ( u s , v s ) E ( s ) × F ( s ) and it satisfy
u t = V ( t , s ) u s + lim n s t P ( t ) V ( t , τ ) Γ n f ( τ , u τ , v τ ) d τ ,
v t = V ( t , s ) v s + lim n s t Q ( t ) V ( t , τ ) Γ n f ( τ , u τ , v τ ) d τ ,
for t s , here u ( t ) and v ( t ) are called stable solution and unstable solution respectively.
Since we want our manifold to be in the form of a graph of some C 1 function, therefore let us consider a space of C 1 functions χ consisting of Φ ( s , · ) : E ( s ) C γ such that for each s 0 ,
1.
Φ ( s , 0 ) = 0 ; ( Φ / ϕ ) ( s , 0 ) = 0 and Φ ( s , E ( s ) ) F ( s ) .
2.
For every ϕ , ψ E ( s ) and for j = 0 , 1 ,
j Φ ϕ j ( s , ϕ ) j Φ ϕ j ( s , ψ ) C γ ϕ ψ C γ .
Using the result of [22], χ is a Banach space with the norm,
| Φ | : = sup Φ ( s , ϕ ) C γ ϕ C γ : s 0 and ϕ E ( s ) \ { 0 } .
Given Φ χ , we consider the graph
W Φ = ( s , ϕ , Φ ( s , ϕ ) ) : ( s , ϕ ) [ 0 , ) × E ( s ) .
Furthermore, for each κ R + , let Ψ κ be the semiflow generated by the autonomous equation
t = 1 , x = A x ( t ) + L x t + f ( t , x t ) .
Given κ = t s 0 and ( s , u s , v s ) R + × E ( s ) × F ( s ) , we have
Ψ κ ( s , u s , v s ) = ( s + κ , u s + κ , v s + κ ) = ( t , u t , v t ) ,
where u t and v t are solutions on stable and unstable subspaces respectively given by Equations (10) and (11).
Here are our assumptions to obtain the C 1 manifold for the Equation (1).
(H1)
a + ϵ 0 .
(H2)
0 < M = sup s 0 { s b ( τ ) e ϵ ρ ( τ ) d τ } < with 4 K M < 1 .
(H3)
0 < M ¯ = sup s 0 { e ( b a + ϵ ) ρ ( s ) s e ( a b + ϵ ) ρ ( τ ) b ( τ ) d τ } < with 12 K 2 M ¯ < 1 .
In the following Theorem, we give the existence of a C 1 stable invariant manifold for the perturbed Equation (1).
Theorem 1.
Assume that the linear Equation (2) satisfies ρ-nonuniform exponential dichotomy with ( H 1 ) and the perturbation in (1) satisfies (9) with ( H 2 ) , ( H 3 ) , then there exists a C 1 function Φ χ such that the set W Φ , defined by the Equation (14), is forward invariant under the semiflow Ψ κ , in the sense that for each κ 0 , Ψ κ ( W Φ ) W Φ . Furthermore, for every κ = t s 0 ; ϕ , ψ E ( s ) there exists D > 0 such that
Ψ κ ( s , ϕ , Φ ( s , ϕ ) ) Ψ κ ( s , ψ , Φ ( s , ψ ) ) C γ D e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ ψ C γ .
Outline: We want our manifold to be invariant under the semiflow, this means that the trajectory of the solutions should remain in the manifold for all time t 0 , provided it starts in the manifold. Since our choice of manifold is graph of some C 1 function Φ , therefore, in the manifold, the solution x ( t ) must take the form, x ( t ) = ( u ( t ) , Φ ( t , u ( t ) ) ) E ( t ) × F ( t ) with ( u s , v s ) E ( s ) × F ( s ) where u ( t ) , Φ ( t , u ( t ) ) satisfy
u t = V ( t , s ) u s + lim n s t P ( t ) V ( t , τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ ,
Φ ( t , u t ) = V ( t , s ) v s + lim n s t Q ( t ) V ( t , τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ .
Now, we are going to prove some lemmas which will be helpful in giving the existence of C 1 smooth manifold. In our first lemma, we establish the existence of solution in the stable direction given by Equation (17).
Lemma 1.
Given ( s , ϕ , Φ ) [ 0 , ) × E ( s ) × χ , there exists a unique C 1 function u : ( , ) × E ( s ) X with u s = ϕ and u t E ( t ) for each t s . The function u t also satisfies Equation (17) for every t s and for each ϕ E ( s ) the segment u t satisfies:
u t C γ 2 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ C γ .
Proof. 
Let Ω * be the space of all C 1 functions u : ( , ) × E ( s ) X such that the following properties satisfies:
(1)
u s = ϕ , and u t E ( t ) for every t s and for each ϕ E ( s ) with the norm,
u * : = sup t s , ϕ E ( s ) \ { 0 } u t C γ ϕ C γ e σ ( t ) 2 K .
(2)
The function u Ω * satisfies
α ( u t ) : = sup t s , ϕ E ( s ) \ { 0 } u t ϕ e σ ( t ) 2 K ,
β ( u t ) : = sup t s , ϕ 1 , ϕ 2 E ( s ) \ { 0 } u t ϕ ( ϕ 1 ) u t ϕ ( ϕ 2 ) e σ ( t ) ϵ ρ ( s ) ϕ 1 ϕ 2 2 K ,
where σ ( t ) = a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) . It follows from result in [22] that Ω * is a complete metric space with norm in (20). Given ( s , ϕ ) [ 0 , ) × E ( s ) and Φ χ , we define an operator L : Ω * Ω * as
( L u ) t = V ( t , s ) ϕ + lim n s t V ( t , τ ) P ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ .
Note that ( L u ) s = ϕ and ( L u ) t E ( t ) for every t s . Now using Equations (5), (7), (9), (12) and (20), we have
( L u ) t C γ V ( t , s ) ϕ C γ + lim n s t V ( t , τ ) P ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) C γ d τ K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ C γ + s t K e a ( ρ ( t ) ρ ( τ ) ) + ϵ ρ ( τ ) b ( τ ) ( u τ , Φ ( τ , u τ ) ) C γ d τ K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ C γ + 2 K s t e a ( ρ ( t ) ρ ( τ ) ) + ϵ ρ ( τ ) b ( τ ) u τ C γ d τ 2 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ C γ 1 2 + 2 K u * s t b ( τ ) e ϵ ρ ( τ ) d τ .
Therefore,
( L u ) * 1 2 + 2 K M u * < .
Hence, L ( Ω * ) Ω * . Now, let u , v Ω * , then again using Equations (5), (7), (9), (12) and (20),
( L u ) t ( L v ) t C γ s t K e a ( ρ ( t ) ρ ( τ ) ) + ϵ ρ ( τ ) f ( τ , u τ , Φ ( τ , u τ ) ) f ( τ , v τ , Φ ( τ , v τ ) ) X d τ 2 K s t e a ( ρ ( t ) ρ ( τ ) ) + ϵ ρ ( τ ) b ( τ ) u τ v τ C γ d τ 4 K 2 ϕ C γ e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) u v * s t b ( τ ) e ϵ ρ ( τ ) d τ .
Furthermore we can write, L u L v * 2 K M u v * .
Since K M < 1 4 , therefore L is a contraction map in Ω * which ensures a unique fixed point function u Ω * such that L u = u . Using Equation (23), we have,
u * 1 2 + 2 K M u * .
Hence, for every t s ,
u t C γ 2 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ C γ .
Next, we give the alternate form for our manifold map.
Lemma 2.
Given Φ χ , let u t be the stable solution given by Lemma 1 for ( s , ϕ ) [ 0 , ) × E ( s ) . Then, for Φ, the following properties holds,
1.
If ( s , ϕ ) [ 0 , ) × E ( s ) and for every t s
Φ ( t , u t ) = V ( t , s ) Φ ( s , ϕ ) + lim n s t V ( t , τ ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ ,
then for all ( s , ϕ ) [ 0 , ) × E ( s ) ,
Φ ( s , ϕ ) = lim n s V Q 1 ( τ , s ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ .
2.
If Equation (25) holds for all ( s , ϕ ) [ s , ) × E ( s ) , then Equation (24) holds for all t s and u t = u t ( s , ϕ , Φ ) w h i l e ( s , ϕ ) [ s , ) × E ( s ) .
Proof. 
Let us first show the validity of the integral in (25) for each fixed n > 0 . Using ρ -nonuniform exponential dichotomy and Equations (9), (12), (19), we have
s V Q 1 ( τ , s ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ C γ s K e b ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( τ ) f ( τ , u τ , Φ ( τ , u τ ) ) X d τ 2 K s e b ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( τ ) b ( τ ) u τ C γ d τ 4 K 2 ϕ C γ s b ( τ ) e ( b + ϵ + a ) ρ ( τ ) e ( b a + ϵ ) ρ ( s ) d τ 4 K 2 ϕ C γ M ¯ < .
Hence, for each fixed n > 0 , the integral in Equation (25) is valid. Now if Equation (24) holds for every ( s , ϕ ) [ 0 , ) × E ( s )   and   t s , then
Φ ( s , ϕ ) = V Q 1 ( t , s ) Φ ( t , u t ) lim n s t V Q 1 ( τ , s ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ .
Using Equations (7), (12) and (19),
V Q 1 ( t , s ) Φ ( t , u t ) C γ K e b ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( t ) u t C γ 2 K 2 e b ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( t ) e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ C γ 2 K 2 e ( b a + ϵ ) ρ ( s ) e ( a b + ϵ ) ρ ( t ) ϕ C γ .
Note that ( a + ϵ b ) < 0 , and taking t proves our result. Additionally, note that the Equation (25) is obtained by operating the invertible map V Q 1 ( t , s ) on the solution x ( t ) = ( u ( t ) , Φ ( t , u ( t ) ) ) , therefore there is no issue of convergence here.
Now assume Equation (25) holds for all ( s , ϕ ) [ 0 , ) × E ( s ) . Note that ( t , u t ) [ s , ) × E ( t ) , therefore it follows from Equation (25),
V ( t , s ) Φ ( s , ϕ ) = lim n s V ( t , τ ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ = lim n s t V ( t , τ ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ lim n t V 1 ( τ , t ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ = lim n s t V ( t , τ ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ + Φ ( t , u t ) .
Hence we get,
Φ ( t , u t ) = V ( t , s ) Φ ( s , ϕ ) + lim n s t V ( t , τ ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ .
The next result shows the dependency of the stable solution u t on the history function ϕ C γ .
Lemma 3.
Let the stable solutions u t , v t be obtained using Lemma 1 for ( s , ϕ , Φ ) and ( s , ψ , Φ ) [ 0 , ) × E ( s ) × χ respectively. We have the following estimates for every t s :
u t v t C γ 2 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ ψ C γ .
Proof. 
Using Lemma 1 and Equations (7), (9), (19) in (17),
u t v t C γ V ( t , s ) ( ϕ ψ ) C γ + s t V ( t , τ ) P ( τ ) f ( τ , u τ , Φ ( τ , u τ ) ) f ( τ , v τ , Φ ( τ , v τ ) ) X d τ K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ ψ C γ + 2 K s t e a ( ρ ( t ) ρ ( τ ) ) + ϵ ρ ( τ ) b ( τ ) u τ v τ C γ d τ 2 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ ψ C γ 1 2 + 2 K u v * s t e ϵ ρ ( τ ) b ( τ ) d τ .
Therefore we have,
u v * 1 2 + 2 K M u v * , u v * 1 2 ( 1 2 K M ) < 1 .
Hence we get the desired result,
u t v t C γ 2 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ ψ C γ , for   every   t s .
Lemma 4.
Let the stable solutions u t , v t be obtained using Lemma 1 for ( s , ϕ , Φ ) and ( s , ϕ , Ψ ) [ 0 , ) × E ( s ) × χ respectively. We have the following estimates for every t s :
u t v t C γ 2 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) | Φ Ψ | ϕ C γ .
Proof. 
Using Equations (12) and (13), we have
Φ ( τ , u τ ) Ψ ( τ , v τ ) C γ u τ v τ C γ + | Φ Ψ | v τ C γ .
Now, using above estimate and Equations (7), (9) and (19) in (17),
u t v t C γ s t K e a ( ρ ( t ) ρ ( τ ) ) + ϵ ρ ( τ ) b ( τ ) ( u τ , Φ ( τ , u τ ) ) ( v τ , Ψ ( τ , v τ ) ) C γ d τ K s t e a ( ρ ( t ) ρ ( τ ) ) + ϵ ρ ( τ ) b ( τ ) 2 u τ v τ C γ + | Φ Ψ | v τ C γ d τ 2 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ C γ 2 u v * + | Φ Ψ | K s t b ( τ ) e ϵ ρ ( τ ) d τ .
Therefore we have,
u v * 2 u v * + | Φ Ψ | K M , u v * K M 1 2 K M | Φ Ψ | .
Since 2 K M < 1 / 2 , therefore, for every t s , we get our desired result:
u t v t C γ 2 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) | Φ Ψ | ϕ C γ .
The following result proves the existence of a C 1 smooth map Φ χ satisfying Equation (25).
Lemma 5.
Let the assumptions a + ϵ 0 , ( H 2 ) and ( H 3 ) hold. Then, there exists a unique function Φ χ such that Equation (25) holds for each ( s , ϕ ) [ 0 , ) × E ( s ) .
Proof. 
Consider an operator J : χ χ given by,
( J Φ ) ( s , ϕ ) : = lim n s V Q 1 ( τ , s ) Q ( τ ) Γ n f ( τ , u τ , Φ ( τ , u τ ) ) d τ ,
for each ( s , ϕ ) [ 0 , ) × E ( s ) , where u t is the unique function given by Lemma 1 for ( s , ϕ , Φ ) . Using the Equation (19), ( J Φ ) ( s , 0 ) = 0 , s > 0 . For s [ 0 , ) and ϕ , ψ E ( s ) , let u t , v t denote the stable solutions obtained using Lemma 1 for ( s , ϕ , Φ ) and ( s , ψ , Φ ) respectively. Now, using Equations (7), (9), (12) and (26) and Lemma 3 we have,
( J Φ ) ( s , ϕ ) ( J Φ ) ( s , ψ ) C γ s K e b ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( τ ) b ( τ ) ( u τ , Φ ( τ , u τ ) ) ( v τ , Φ ( τ , v τ ) ) C γ d τ 2 K s e b ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( τ ) b ( τ ) u τ v τ C γ d τ 4 K 2 ϕ ψ C γ s e ( a b + ϵ ) ρ ( τ ) b ( τ ) e ( b a + ϵ ) ρ ( s ) d τ 4 K 2 M ¯ ϕ ψ C γ .
Since 4 K 2 M ¯ < 1 ,
( J Φ ) ( s , ϕ ) ( J Φ ) ( s , ψ ) C γ ϕ ψ C γ .
Hence, J ( χ ) χ . For ( s , ϕ ) [ 0 , ) × E ( s ) ; Φ , Ψ χ and u t , v t are the corresponding stable solutions, then,
( J Φ ) ( s , ϕ ) ( J Ψ ) ( s , ϕ ) C γ s K e b ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( τ ) b ( τ ) ( u τ v τ , Φ ( τ , u τ ) Ψ ( τ , v τ ) ) C γ d τ K s e b ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( τ ) b ( τ ) 2 u τ v τ C γ + | Φ Ψ | v τ C γ d τ 4 K 2 | Φ Ψ | ϕ C γ e ( b a + ϵ ) ρ ( s ) s e ( a b + ϵ ) ρ ( τ ) b ( τ ) d τ 4 K 2 M ¯ | Φ Ψ | ϕ C γ .
Since 4 K 2 M ¯ < 1 , hence, J is a contraction map which gives the existence of a unique fixed function Φ χ such that J Φ = Φ .
Now we give the proof of Theorem 1 using the lemmas proved in this section.
Proof. 
To prove Theorem 1 we need to find a function Φ χ satisfying Equations (24) and (25), also, the graph of Φ , i.e., W Φ , should be invariant under the semiflow Ψ κ given by Equation (15). For each ( s , ϕ ) [ 0 , ) × E ( s ) , Lemma 5 gives the existence of Φ χ satisfying Equations (24) and (25). Furthermore, from Lemma 1 for each ( s , ϕ , Φ ) [ 0 , ) × E ( s ) × χ there exists a unique function u t satisfying Equation (17). Note that ( t , u t ) [ 0 , ) × E ( t ) , therefore
Ψ t s ( s , ϕ , Φ ( s , ϕ ) ) = ( t , u t , Φ ( t , u t ) ) W Φ , for   all   t s .
Now for each ( s , ϕ ) , ( s , ψ ) [ 0 , ) × E ( s ) and κ = t s 0 by Lemma 3, Lemma 4 and Equation (15) we have,
Ψ κ ( s , ϕ , Φ ( s , ϕ ) ) Ψ κ ( s , ψ , Φ ( s , ψ ) ) ( t , u t ϕ , Φ ( t , u t ϕ ) ) ( t , u t ψ , Φ ( t , u t ψ ) ) 2 u t ϕ u t ψ C γ 4 K e a ( ρ ( t ) ρ ( s ) ) + ϵ ρ ( s ) ϕ ψ C γ .
This completes the proof of the theorem. □

4. Stable Manifold and Perturbations

In this section we are showing how the manifold varies with the change in perturbation. We defined the space of perturbation Ω in the beginning of Section 3. Now we define a norm on Ω , for f Ω ,
f Ω : = sup f ( t , ϕ ) X ϕ C γ , t 0 , ϕ C γ \ { 0 } .
We also assume one more condition:
(H4)
0 < M ^ = sup s 0 { e ( b a + ϵ ) ρ ( s ) s e ( a b + ϵ ) ρ ( τ ) d τ } with 8 K 2 M ^ < 1 .
Theorem 2.
Assume that the differential Equation (1) admits ρ-nonuniform exponential dichotomy and also the assumptions in ( H 1 ) ( H 4 ) hold, then for each set ( s , ϕ , f ) and ( s , ϕ , g ) the maps Φ and Ψ satisfies
| Φ Ψ | f g Ω .
Proof. 
For every ( s , ϕ , f )   and   ( s , ϕ , g ) , we obtain unique functions Φ and Ψ using Lemma 5. Now using Lemma 1 for ( s , ϕ , f , Φ ) and ( s , ϕ , g , Ψ ) we obtain unique functions u Φ = u and u Ψ = v . From Equation (25) we have,
Φ ( s , ϕ ) Ψ ( s , ϕ ) C γ s V Q 1 ( τ , s ) Q ( τ ) f ( τ , u τ , Φ ( τ , u τ ) ) g ( τ , v τ , Ψ ( τ , v τ ) ) C γ d τ s K e b ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( τ ) f ( τ , u τ , Φ ( τ , u τ ) ) g ( τ , v τ , Ψ ( τ , v τ ) ) X d τ K s e b ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( τ ) b ( τ ) 2 u τ v τ C γ + | Φ Ψ | u τ C γ + f g Ω ( u τ , Φ ( τ , u τ ) ) C γ d τ K s e b ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( τ ) 6 K b ( τ ) e a ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( s ) | Φ Ψ | ϕ C γ + 4 K f g Ω e a ( ρ ( τ ) ρ ( s ) ) + ϵ ρ ( s ) ϕ C γ d τ 6 K 2 ϕ C γ | Φ Ψ | M ¯ + 4 K 2 ϕ C γ f g Ω M ^ .
Therefore we have,
| Φ Ψ | 4 K 2 M ^ 1 6 K 2 M ¯ f g Ω .
Since 4 K 2 M ^ 1 6 K 2 M ¯ < 1 , using (H4), hence we got our result. □
Thus we have shown that a small change in the perturbation gives arise to the small variation in manifold map.

5. Examples

In this example, we consider a delay differential equation admitting the dichotomic behavior with the growth rate of type e ρ ( t ) .
Example 3.
x ( t ) = ω ρ ( t ) cos t a + sin t + ϵ ρ ( t ) ρ ( t ) sin ρ ( t ) x ( t ) + e 3 ϵ ( t 1 ) x 2 ( t 1 ) ,
for ω > 2 ϵ > 0 and a > 2 .
For the associated linear problem, the solution x ( t ) is given by x ( t ) = V ( t , s ) x ( s ) , where,
V ( t , s ) = e ω ( ρ ( t ) ρ ( s ) ) ln ( a + sin t ) + ln ( a + sin s ) + ϵ ( ρ ( s ) cos ρ ( s ) ρ ( t ) cos ρ ( t ) ) + ϵ ( sin ρ ( t ) sin ρ ( s ) ) .
Additionally,
V ( t , s ) = e ω ( ρ ( t ) ρ ( s ) ) ln ( a + sin t ) + ln ( a + sin s ) + ϵ ( ρ ( s ) cos ρ ( s ) ρ ( t ) cos ρ ( t ) ) + ϵ ( sin ρ ( t ) sin ρ ( s ) ) e 2 ϵ + 2 ( ln ( a + 1 ) ) e ( ω + ϵ ) ( ρ ( t ) ρ ( s ) ) + 2 ϵ ρ ( s ) = C e ( ω + ϵ ) ( ρ ( t ) ρ ( s ) ) + 2 ϵ ρ ( s ) .
Note that ( ω + ϵ ) < 0 and C is a constant, therefore the solution operator V ( t , s ) admits dichotomic behavior with growth rate e ρ ( t ) . Additionally, the nonlinear perturbation satisfies Equation (9). Hence, Theorem 1 ensures that the differential Equation (32) admits C 1 stable invariant manifold.
Now, we consider a system of delay differential equations which satisfies all the assumptions of Theorem 1 and admits C 1 stable invariant manifold.
Example 4.
x ( t ) = ω ρ ( t ) 1 + e ρ ( t ) + ϵ ρ ( t ) ρ ( t ) sin ρ ( t ) 1 x ( t ) + sin ( t 1 ) e 5 ϵ ( t 1 ) y 2 ( t 1 ) , y ( t ) = ω ρ ( t ) 1 + e ρ ( t ) ϵ ρ ( t ) ρ ( t ) cos ρ ( t ) 1 y ( t ) + cos ( t 1 ) e 7 ϵ ( t 1 ) x 2 ( t 1 ) ,
for some positive constants ω, ϵ and for ϕ = ( ϕ 1 , ϕ 2 ) ; ϕ 1 , ϕ 2 C γ also the function ρ ( t ) satisfies the conditions in (7). Let the non homogeneous term be,
f ( t , ϕ ) = sin t e 5 ϵ t ϕ 2 2 ( t ) , cos t e 7 ϵ t ϕ 1 2 ( t ) .
The evolution operator for the associated linear system is given by
T ( t , s ) = U ( t , s ) 0 0 V ( t , s ) ,
where
U ( t , s ) = e ω ( ρ ( t ) ρ ( s ) ) ω ln ( 1 + e ρ ( t ) ) ( 1 + e ρ ( s ) ) + ϵ ( ρ ( s ) cos ρ ( s ) ρ ( t ) cos ρ ( t ) ) + ϵ ( sin ρ ( t ) sin ρ ( s ) ) ϵ ( ρ ( t ) ρ ( s ) )
and
V ( t , s ) = e ω ( ρ ( t ) ρ ( s ) ) + ω ln ( 1 + e ρ ( t ) ) ( 1 + e ρ ( s ) ) + ϵ ( ρ ( s ) sin ρ ( s ) ρ ( t ) sin ρ ( t ) ) + ϵ ( cos ρ ( t ) cos ρ ( s ) ) ϵ ( ρ ( t ) ρ ( s ) ) .
Consider the projection map P ( t ) ( x , y ) = ( x , 0 ) for all t 0 . Then,
T ( t , s ) P ( t ) = U ( t , s ) e 2 ϵ + ω ln 2 e ω ( ρ ( t ) ρ ( s ) ) + 2 ϵ ρ ( s ) ,
and
T ( t , s ) 1 Q ( t ) = V ( t , s ) 1 e 2 ϵ + ω ln 2 e ω ( ρ ( t ) ρ ( s ) ) + 2 ϵ ρ ( t ) .
for t s 0 . Hence for a = ω < 0 and b = ω > 0 , the associated linear system satisfies ρ-nonuniform exponential dichotomy.
Furthermore, f ( t , ϕ ) admits the conditions in Equation (9), with f ( t , 0 ) = 0 . Hence, Theorem 1 ensures that the differential Equation (4) admits C 1 stable invariant manifold.

6. Conclusions

In this article, we have constructed a C 1 stable invariant manifold for the differential equation with infinite delay (1). We gave the dependence of manifold on the perturbations. We have also included general examples of differential equations which satisfy ρ -nonuniform exponential dichotomy. Thus, we have shown that our result is applicable to a large class of delay differential equations.

Author Contributions

Both the authors worked together to conceptualize and investigate the project. First author has written the article with the guidance of second author. Both authors have read and agreed to the published version of the manuscript.

Funding

This work is supported in part by the Croatian Science Foundation under the project IP-2019-04-1239 (Lokesh Singh) and by the Indian Institute of Technology Kanpur (D. Bahuguna), MHRD, India.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hadamard, J. Sur l’itération et les solutions asymptotiques des équations différentielles. Bull. Soc. Math. Fr. 1901, 29, 224–228. [Google Scholar]
  2. Perron, O. Die stabilitätsfrage bei differentialgleichungen. Math. Z. 1930, 32, 703–728. [Google Scholar] [CrossRef]
  3. Liapounoff, M.A. Problème Général de la Stabilité du Mouvement; Princeton University Press: Princeton, NJ, USA, 1949. [Google Scholar]
  4. Hirsch, M.W.; Pugh, C.C. Stable manifolds and hyperbolic sets. Global Analysis (Proceedings of the Symposia in Pure Mathematics, Volume XIV, Berkeley, CA, USA, 1968. 1970, pp. 133–163. Available online: https://escholarship.org/uc/item/2rr3r1pz (accessed on 18 April 2021).
  5. Hirsch, M.; Pugh, C.; Shub, M. Invariant Manifolds; Volume 583 of Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 1977. [Google Scholar]
  6. Sacker, R.J.; Sell, G.R. Existence of dichotomies and invariant splittings for linear differential systems I. J. Differ. Equ. 1974, 15, 429–458. [Google Scholar] [CrossRef] [Green Version]
  7. Pesin, J.B. Families of invariant manifolds corresponding to nonzero characteristic exponents. Math. USSR Izv. 1976, 10, 1261. [Google Scholar] [CrossRef]
  8. Pesin, Y.B. Characteristic Lyapunov exponents and smooth ergodic theory. Russ. Math. Surv. 1977, 32, 55–114. [Google Scholar] [CrossRef]
  9. Pugh, C. Shub M.: Ergodic attractors Trans. Amer. Math. Soc. 1989, 312, 1–54. [Google Scholar] [CrossRef]
  10. Ruelle, D. Characteristic exponents and invariant manifolds in Hilbert space. Ann. Math. 1982, 115, 243–290. [Google Scholar] [CrossRef]
  11. Barreira, L.; Valls, C. Stable manifolds for nonautonomous equations without exponential dichotomy. J. Differ. Equ. 2006, 221, 58–90. [Google Scholar] [CrossRef] [Green Version]
  12. Barreira, L.; Valls, C. Smoothness of invariant manifolds for nonautonomous equations. Commun. Math. Phys. 2005, 259, 639–677. [Google Scholar] [CrossRef]
  13. Barreira, L.; Valls, C. Stability of Nonautonomous Differential Equations; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  14. Barreira, L.; Valls, C. Growth rates and nonuniform hyperbolicity. Discret. Contin. Dyn. Syst. A 2008, 22, 509–528. [Google Scholar] [CrossRef]
  15. Barreira, L.; Valls, C. Stable invariant manifolds for parabolic dynamics. J. Funct. Anal. 2009, 257, 1018–1029. [Google Scholar] [CrossRef] [Green Version]
  16. Pan, L. Stable invariant manifolds with impulses and growth rates. Adv. Differ. Equ. 2015, 2015, 221. [Google Scholar] [CrossRef] [Green Version]
  17. Barreira, L.; Valls, C. Parameter dependence of stable manifolds under nonuniform hyperbolicity. J. Math. Anal. Appl. 2009, 358, 419–426. [Google Scholar] [CrossRef] [Green Version]
  18. Barreira, L.; Valls, C. Smooth stable invariant manifolds and arbitrary growth rates. Nonlinear Anal. Theory Methods Appl. 2010, 72, 2444–2456. [Google Scholar] [CrossRef]
  19. Singh, L.; Bahuguna, D. Smooth Invariant Manifolds for Differential Equations with Infinite Delay. J. Dyn. Control Syst. 2021, 27, 107–132. [Google Scholar] [CrossRef]
  20. Hino, Y.; Murakami, S.; Naito, T.; Van Minh, N. A variation-of-constants formula for abstract functional differential equations in the phase space. J. Differ. Equ. 2002, 179, 336–355. [Google Scholar] [CrossRef] [Green Version]
  21. Bahuguna, D.; Singh, L. Stable manifolds for impulsive delay equations and parameter dependence. Electron. J. Differ. Equ. 2019, 1–22. [Google Scholar]
  22. Elbialy, M.S. On sequences of Ck,δb maps which converge in the uniform C0-norm. Proc. Am. Math. Soc. 2000, 128, 3285–3290. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Singh, L.; Bahuguna, D. Smooth Stable Manifold for Delay Differential Equations with Arbitrary Growth Rate. Axioms 2021, 10, 105. https://doi.org/10.3390/axioms10020105

AMA Style

Singh L, Bahuguna D. Smooth Stable Manifold for Delay Differential Equations with Arbitrary Growth Rate. Axioms. 2021; 10(2):105. https://doi.org/10.3390/axioms10020105

Chicago/Turabian Style

Singh, Lokesh, and Dhirendra Bahuguna. 2021. "Smooth Stable Manifold for Delay Differential Equations with Arbitrary Growth Rate" Axioms 10, no. 2: 105. https://doi.org/10.3390/axioms10020105

APA Style

Singh, L., & Bahuguna, D. (2021). Smooth Stable Manifold for Delay Differential Equations with Arbitrary Growth Rate. Axioms, 10(2), 105. https://doi.org/10.3390/axioms10020105

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop