INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS OF NONAUTONOMOUS DIFFERENTIAL EQUATIONS Bernd Aulbach and Martin Rasmussen1 Department of Mathematics University of Augsburg D-86135 Augsburg, Germany
Stefan Siegmund2 Department of Mathematics University of Frankfurt D-60325 Frankfurt am Main, Germany
Abstract. We discuss the relationship between invariant manifolds of nonautonomous differential equations and pullback attractors. This relationship is essential, e.g., for the numerical approximation of these manifolds. In a first step we show that the unstable manifold is the pullback attractor of the differential equation. The main result says that every (hyperbolic or nonhyperbolic) invariant manifold is the pullback attractor of a related system which we construct explicitly using spectral transformations. To illustrate our theorem we present an application to the Lorenz system and approximate numerically the stable as well as the strong stable manifold of the origin.
1. Introduction. In the theory of dynamical systems (e.g., autonomous differen´ tial equations) the stable and unstable manifold theorem – first proved by Poincare [12] and Hadamard [10] – plays an important role in the study of the flow near a hyperbolic rest point or a periodic solution. The idea is to show that the rest point x = 0, say, locally inherits the dynamical structure of its hyperbolic linearization x˙ = Ax, i.e., the spectrum σ(A) = σ s (A) ∪ σ u (A) = {λ1 , . . . , λr } ∪ {λr+1 , . . . , λn } consists of eigenvalues λ1 , . . . , λr with negative real parts and λr+1 , . . . , λn with positive real parts. Let V s and V u denote the sum of all generalized eigenspaces corresponding to the eigenvalues in σ s (A) and σ u (A), respectively. Then V s and V u are the stable and unstable manifold of the linear system, respectively, and the stable manifold theorem provides nonlinear analogues W s = {x0 : x(t, x0 ) → 0 as t → ∞}
and
W u = {x0 : x(t, x0 ) → 0 as t → −∞}
which are tangential to V s and V u , respectively; x(t, x0 ) being the general solution of the nonlinear system. The fact that the linear manifolds V s and V u are direct sums of eigenspaces consisting of solutions with an exponential growth or decay rate given by the real parts of the corresponding eigenvalues also carries over to W s and 2000 Mathematics Subject Classification. 34C30, 34D45, 37B55, 37D10, 65L20. Key words and phrases. Nonautonomous differential equation, Pullback attractor, Invariant manifold, Numerical approximation. 1 Supported by the Deutsche Forschungsgemeinschaft, grant GK 283. 2 Supported by the Deutsche Forschungsgemeinschaft, grant Si801.
1
2
B. AULBACH, M. RASMUSSEN AND S. SIEGMUND
W u . A submanifold of W s which corresponds to some (but not all) eigenvalues with smallest real part is called a strong stable manifold. As for an example, the Lorenz system at the origin (for the usual parameters) has two negative eigenvalues λ1 < λ2 < 0; the unstable manifold therefore has dimension 2 whereas the strong stable manifold (consisting of solutions with decay rate λ1 ) has dimension 1. The analytic computation of invariant manifolds is possible only in rare cases. It is reasonable, therefore, to develop numerical tools for the approximation of this kind of sets. For stable and unstable manifolds of autonomous differential equations this topic is well examined. In Dellnitz and Hohmann [6, 7], a subdivision as well as a continuation algorithm is introduced to approximate the global attractor of the system; and since the global attractor contains all unstable manifolds, these algorithms are suitable for the approximation of unstable manifolds. In this paper, we go one step further and develop an approximation result for more general types of invariant manifolds. In fact, we consider all invariant manifolds of a so-called hierarchy of invariant manifolds which are nonlinear analogues of the eigenspaces of the linearization, including strong stable manifolds and nonhyperbolic manifolds such as center or center-stable manifolds. We make use of the observation that (under suitable global assumptions) the unstable manifold is globally attracting – and therefore amenable to a numerical approximation. Then we use spectral transformations to relate an arbitrary invariant manifold of the above-mentioned hierarchy to the unstable manifold of a transformed system. The price one has to pay for this wide range of applications is that the transformed systems are nonautonomous even if the original system is autonomous. This fact is the major reason for dealing with nonautonomous systems from the outset. This paper is organized as follows. In Section 2 we introduce some notation and define nonautonomous attractors. Section 3 is devoted to nonautonomous invariant manifolds. The relationship between hyperbolic invariant manifolds and attractors is established in Section 4. In Section 5 we describe the nonhyperbolic case leading to the general result in Section 6. To illustrate our theorem, in Section 7 we compute and visualize the stable as well as the strong stable manifold of the origin of the Lorenz system using a nonautonomous extension of the software package GAIO. 2. Preliminaries. As usual we denote by R the set containing all reals and by RN ×N the set of all real N × N matrices. We use the symbol 1 for the unit matrix. With N X k(x1 , . . . , xN )k := |xi | for all (x1 , . . . , xN ) ∈ RN i=1 N
the R is a normed vector space. We write U² (x0 ) = {x ∈ RN : kx − x0 k < ²} for the ²-neighborhood of a point x0 ∈ RN . For arbitrary nonempty sets A, B ⊂ RN and x ∈ RN let d(x, A) := inf{kx − yk : y ∈ A} be the distance of x to A and d(A|B) := sup{d(x, B) : x ∈ A} be the Hausdorff semi-distance of A and B. A function g : R → RN is called γ + -quasibounded if −γt kgk+ : t ≥ τ} < ∞ τ,γ := sup{kg(t)ke
for some τ ∈ R . Accordingly, the function g is called γ − -quasibounded if −γt kgk− : t ≤ τ} < ∞ τ,γ := sup{kg(t)ke
for some τ ∈ R . In this article we are concerned with nonautonomous differential equations x˙ = f (t, x) ,
INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS
3
where f : R × RN → RN is a function which fulfills conditions ensuring global existence and uniqueness of solutions. This equation then gives rise to a general solution λ(t, τ, ξ) which is the solution satisfying the initial condition x(τ ) = ξ with τ and ξ treated as additional parameters. For a nonautonomous linear differential equation x˙ = B(t)x
(1)
with a continuous function B : R → RN ×N the general solution is linear with respect to ξ, i.e., there exists a function Φ : R × R → RN ×N , the so-called evolution operator of (1), such that λ(t, τ, ξ) = Φ(t, τ )ξ
for all t, τ ∈ R and ξ ∈ RN .
We say that a function P : R → RN ×N is an invariant projector of (1) if P (t)2 = P (t) for all t ∈ R , P (t)Φ(t, s) = Φ(t, s)P (s) for all t, s ∈ R . A subset A of R × RN is called a nonautonomous set if for all t ∈ R the so-called t-fibers A(t) := {x ∈ RN : (t, x) ∈ A} are nonempty. We call A closed or compact, if all t-fibers are closed or compact, respectively. Finally, a nonautonomous set A is called invariant if λ(t, τ, A(τ )) = A(t) for all t, τ ∈ R . In the literature pullback attractors are usually defined to be compact (see, e.g., Kloeden & Keller & Schmalfuß [11]), but the (global) nonautonomous manifolds under consideration here are always noncompact. If we want to establish connections between pullback attractors and nonautonomous manifolds we need a more general notion of pullback attractor. This is prepared by the following definition: an invariant nonautonomous set A is said to be compactly generated if there exists a compact set K ⊂ RN , a so-called generator of A, with the following property: for any compact set C ⊂ RN there exists a number T (K, C) > 0 such that for any τ ∈ R we have λ(τ, τ − t, A(τ − t) ∩ K) ⊃ A(τ ) ∩ C
for all t > T (K, C) .
An invariant, closed and compactly generated nonautonomous set A is called a global pullback attractor if for any compact set C ⊂ X we have lim d(λ(τ, t, C)|A(τ )) = 0
t→−∞
for all τ ∈ R .
A global pullback attractor is always unique. We have proved this result in [1], where we have also introduced algorithms to approximate pullback attractors. 3. Theory of nonautonomous manifolds. In the sequel we consider nonautonomous differential equations of the form x˙ 1 = B1 (t)x1 + F1 (t, x1 , x2 ) x˙ 2 = B2 (t)x2 + F2 (t, x1 , x2 ) ,
(2)
where B1 : R → RN ×N , B2 : R → RM ×M , F1 : R × RN × RM → RN and F2 : R × RN × RM → RM are continuous functions satisfying F1 (t, 0, 0) = 0 and F2 (t, 0, 0) = 0 for all t ∈ R. Moreover we assume: (H1) Hypotheses on linear part: The evolution operators Φ1 and Φ2 of the linear equations x˙ 1 = B1 (t)x1 and x˙ 2 = B2 (t)x2 , respectively, satisfy the estimates kΦ1 (t, s)k ≤ Keα(t−s) kΦ2 (t, s)k ≤ Ke
β(t−s)
for all t ≥ s , for all t ≤ s
4
B. AULBACH, M. RASMUSSEN AND S. SIEGMUND
with real constants K ≥ 1 and α < β . (H2) Hypotheses on perturbation: For all (x1 , x2 ), (¯ x1 , x ¯2 ) ∈ RN × RM and t ∈ R we have kF1 (t, x1 , x2 ) − F1 (t, x ¯1 , x ¯2 )k ≤ Lkx1 − x ¯1 k + Lkx2 − x ¯2 k , kF2 (t, x1 , x2 ) − F2 (t, x ¯1 , x ¯2 )k ≤ Lkx1 − x ¯1 k + Lkx2 − x ¯2 k , where the constant L satisfies the estimate 0 ≤ L < We denote the general solution of this system by
β−α 4K
.
λ(t, τ, ξ, η) = (λ1 (t, τ, ξ, η), λ2 (t, τ, ξ, η)) ∈ RN × RM and choose an arbitrary constant δ ∈ (2KL, β−α 2 ]. The theorems in this section are slight modifications of the results obtained in Aulbach & Wanner [3, 4]. First we state the fundamental existence theorem on nonautonomous manifolds which says that system (2) gives rise to two nonautonomous manifolds, the pseudo-stable manifold S0 and the pseudo-unstable manifold R0 . If system (2) is hyperbolic, i.e., α < 0 < β, then S0 and R0 are called stable manifold and unstable manifold, respectively. Theorem 3.1. There exists a uniquely determined continuous mapping s0 : R × RN → RM whose graph S0 := {(τ, ξ, s0 (τ, ξ)) : τ ∈ R, ξ ∈ RN } allows the representation S0 = {(τ, ξ, η) : λ(·, τ, ξ, η) is γ + -quasibounded} for every choice of γ ∈ [α + δ, β − δ]. Moreover, there exists a uniquely determined continuous mapping r0 : R × RM → RN whose graph R0 := {(τ, r0 (τ, η), η) : τ ∈ R, η ∈ RM } allows the representation R0 = {(τ, ξ, η) : λ(·, τ, ξ, η) is γ − -quasibounded} for every choice of γ ∈ [α + δ, β − δ]. The nonautonomous sets S0 and R0 are invariant and their intersection is the trivial solution of (2). We say that S0 and R0 are the nonautonomous manifolds of the trivial solution. Not only the trivial solution but every solution of (2) admits nonautonomous manifolds. This is the central statement of the following theorem. Theorem 3.2. There exists a mapping s : R × RN × R × RN × RM → RM such that for every (τ0 , ξ0 , η0 ) ∈ R × RN × RM the graph Sτ0 ,ξ0 ,η0 := {(τ, ξ, s(τ, ξ, τ0 , ξ0 , η0 )) : τ ∈ R, ξ ∈ RN } allows the representation Sτ0 ,ξ0 ,η0 = {(τ, ξ, η) : λ(·, τ, ξ, η) − λ(·, τ0 , ξ0 , η0 ) is γ + -quasibounded} for every choice of γ ∈ [α + δ, β − δ] . The function s is continuous and has the following properties: (S1) For all τ, τ0 ∈ R, ξ0 , ξ1 , ξ2 ∈ RN and η0 ∈ RM we have ks(τ, ξ1 , τ0 , ξ0 , η0 ) − s(τ, ξ2 , τ0 , ξ0 , η0 )k ≤
K 2 L(δ − KL) kξ1 − ξ2 k . δ(δ − 2KL)
INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS
5
(S2) For all τ, τ0 ∈ R, ξ, ξ0 ∈ RN , η0 ∈ RM and γ ∈ [α + δ, β − δ] we have kλ1 (·, τ, ξ, s(τ, ξ, τ0 , ξ0 , η0 )) − λ1 (·, τ0 , ξ0 , η0 )k+ τ,γ µ ¶ K 3 L2 K+ kξ − λ1 (τ, τ0 , ξ0 , η0 )ke−γτ , δ(δ − 2KL)
≤
kλ2 (·, τ, ξ, s(τ, ξ, τ0 , ξ0 , η0 )) − λ2 (·, τ0 , ξ0 , η0 )k+ τ,γ
≤
2
K L(δ − KL) kξ − λ1 (τ, τ0 , ξ0 , η0 )ke−γτ . δ(δ − 2KL) Moreover, there exists a mapping r : R × RM × R × RN × RM → RN such that for every (τ0 , ξ0 , η0 ) ∈ R × RN × RM the graph Rτ0 ,ξ0 ,η0 := {(τ, r(τ, η, τ0 , ξ0 , η0 ), η) : τ ∈ R, η ∈ RM } allows the representation Rτ0 ,ξ0 ,η0 = {(τ, ξ, η) : λ(·, τ, ξ, η) − λ(·, τ0 , ξ0 , η0 ) is γ − -quasibounded} for every choice of γ ∈ [α + δ, β − δ] . The function r is continuous and has the following properties: (R1) For all τ, τ0 ∈ R, ξ0 ∈ RN and η0 , η1 , η2 ∈ RM we have kr(τ, η1 , τ0 , ξ0 , η0 ) − r(τ, η2 , τ0 , ξ0 , η0 )k ≤
K 2 L(δ − KL) kη1 − η2 k . δ(δ − 2KL)
(R2) For all τ, τ0 ∈ R, ξ0 ∈ RN , η, η0 ∈ RM and γ ∈ [α + δ, β − δ] we have kλ1 (·, τ, r(τ, η, τ0 , ξ0 , η0 ), η) − λ1 (·, τ0 , ξ0 , η0 )k− τ,γ
≤
2
K L(δ − KL) kη − λ2 (τ, τ0 , ξ0 , η0 )ke−γτ , δ(δ − 2KL) kλ2 (·, τ, r(τ, η, τ0 , ξ0 , η0 ), η) − λ2 (·, τ0 , ξ0 , η0 )k+ τ,γ ¶ µ K 3 L2 kη − λ2 (τ, τ0 , ξ0 , η0 )ke−γτ . K+ δ(δ − 2KL)
≤
Remark 3.1. The properties (R1) and (R2) hold also true for the function r0 of Theorem 3.1 instead of r, because of r(·, ·, 0, 0, 0) ≡ r0 (·, ·). This fact is needed later in the proof of Theorem 4.1. The nonautonomous set FH (τ0 , η0 ) := Sτ0 ,r0 (τ0 ,η0 ),η0 is called the horizontal fiber bundle through (τ0 , η0 ). The next theorem says – provided the constant L is small enough – that every point of the extended phase space lies on exactly one horizontal fiber bundle. Theorem 3.3. Suppose that the constant L of system (2) satisfies p β−α (2 + K − 4 + K 2 ) . 0≤L< 2 4K Then there exists a continuous mapping F : R × RN × RM → RM with the following property: every point (τ0 , ξ0 , η0 ) ∈ R × RN × RM lies exactly on one horizontal fiber bundle, namely FH (τ0 , η) with η = F(τ0 , ξ0 , η0 ). Furthermore, for all τ ∈ R, ξ ∈ RN and η ∈ RM the estimate 1 (kξk + kηk) (3) kF(τ, ξ, η)k ≤ 1 − C(L)2
6
B. AULBACH, M. RASMUSSEN AND S. SIEGMUND
holds with the constant C(L) defined by C(L) :=
K 2 L(δ − KL) . δ(δ − 2KL)
4. Hyperbolic manifolds and pullback attractors. We suppose now that system (2) is hyperbolic. For simplicity we refer to (2) also by writing x˙ = f (t, x). In the following theorem we provide sufficient conditions in terms of the spectral gap and the constant L in order to establish connections between the two hyperbolic manifolds and pullback attractors. Theorem 4.1. We suppose that 0 ∈ (α + 2KL, β − 2KL) ,
(4)
and that the constant L satisfies p β−α (2 + K − 0≤L< 4 + K 2) . 4K 2 Then the unstable manifold R0 is the global pullback attractor of (2). Furthermore, we obtain the following connection between the global pullback attractor A of the system under time-reversal x˙ = −f (−t, x)
(5)
and the stable manifold S0 of system (2): A(t) = S0 (−t)
for all t ∈ R .
Proof. First we show that R0 is the global pullback attractor of system (2). Obviously, R0 is an invariant nonautonomous set and graph of the continuous function r0 ; hence a closed nonautonomous set. Let us prove now that R0 is compactly generated with generator U1 (0). For this we choose an arbitrary compact set C ⊂ RN . Since C is bounded, we have an M > 0 with kxk < M for all x ∈ C. Due to (4) there exists a number δ ∈ (2KL, β−α 2 ] with γ := β − δ > 0. Thus we get a T (U1 (0), C) > 0 such that for all t > T (U1 (0), C) we have µ 2 ¶ K L(δ − KL) K 3 L2 +K + M e−γt < 1 . (6) δ(δ − 2KL) δ(δ − 2KL) Choose τ0 ∈ R and (ξ, η) = (r0 (τ0 , η), η) ∈ C ∩ R0 (τ0 ) arbitrarily. Then for all t > T (U1 (0), C) we have kλ(τ0 − t, τ0 , r0 (τ0 , η), η)k = kλ1 (τ0 − t, τ0 , r0 (τ0 , η), η)k + kλ2 (τ0 − t, τ0 , r0 (τ0 , η), η)k ¶ µ 2 (R2) (6) K 3 L2 K L(δ − KL) +K + kηk e−γt < 1 . ≤ δ(δ − 2KL) δ(δ − 2KL) |{z} T (U1 (0), C) the relation {(ξ, η)} = λ(τ0 , τ0 − t, {λ(τ0 − t, τ0 , r0 (τ0 , η), η)}) ⊂ λ(τ0 , τ0 − t, U1 (0)) is fulfilled. Since (ξ, η) ∈ C ∩ R0 (τ0 ) has been chosen arbitrarily we have C ∩ R0 (τ0 ) ⊂ λ(τ0 , τ0 − t, U1 (0)) for all t > T (U1 (0), C) . This means that R0 is compactly generated with generator U1 (0) . Finally we prove that R0 attracts every compact set in the sense of pullback attraction. Therefore
INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS
7
choose τ0 ∈ R and a compact set C ⊂ RN arbitrarily. Due to (4) there exists a number δ ∈ (2KL, β−α 2 ] with γ := α + δ < 0. We see that for all (ξ, η) ∈ C and t ∈ R we have kr0 (t, F(t, ξ, η)) − ξk ≤ kr0 (t, F(t, ξ, η))k + kξk (R1)
≤
(3)
≤
K 2 L(δ − KL) kF(t, ξ, η)k + kξk δ(δ − 2KL) K 2 L(δ − KL) 1 (kξk + kηk) + kξk ≤ M1 δ(δ − 2KL) 1 − C(L)2
with M1 > 0 since C is bounded. Moreover, for all (ξ, η) ∈ C and t ∈ R the relation s(t, r0 (t, F(t, ξ, η)), t, ξ, η) = F(t, ξ, η) is fulfilled since the identity (t, r0 (t, F(t, ξ, η)), F(t, ξ, η)) ∈ St,ξ,η = FH (t, F(t, ξ, η)) holds for all (ξ, η) ∈ C and t ∈ R. Hence, for all (ξ, η) ∈ C and t ≤ τ0 kλ(τ0 , t, r0 (t, F(t, ξ, η)), F(t, ξ, η)) − λ(τ0 , t, ξ, η)k = kλ1 (τ0 , t, r0 (t, F(t, ξ, η)), s(t, r0 (t, F(t, ξ, η)), t, ξ, η)) − λ1 (τ0 , t, ξ, η)k + kλ2 (τ0 , t, r0 (t, F(t, ξ, η)), s(t, r0 (t, F(t, ξ, η)), t, ξ, η)) − λ2 (τ0 , t, ξ, η)k µ ¶ (S2) K 3 L2 ≤ K+ kr0 (t, F(t, ξ, η)) − ξkeγ(τ0 −t) + δ(δ − 2KL) K 2 L(δ − KL) kr0 (t, F(t, ξ, η)) − ξkeγ(τ0 −t) δ(δ − 2KL) µ ¶ K 3 L2 K 2 L(δ − KL) ≤ K+ + M1 eγ(τ0 −t) =: M2 eγ(τ0 −t) δ(δ − 2KL) δ(δ − 2KL) with M2 > 0. We choose an ² > 0. Then there exists a t˜ > 0 with M2 eγ t˜ < ². Thus, for all (ξ, η) ∈ C and t ≥ t˜ d(λ(τ0 , τ0 − t, ξ, η), R0 (τ0 )) ≤ kλ(τ0 , τ0 − t, r0 (τ0 − t, F(τ0 − t, ξ, η)), F(τ0 − t, ξ, η)) − λ(τ0 , τ0 − t, ξ, η)k ≤ M2 eγt < ² . This leads immediately to d(λ(τ0 , τ0 − t, C)|R0 (τ0 )) < ² for all t ≥ t˜ and finishes the proof that R0 is the global pullback attractor of (2). The verification of the second assertion of Theorem 4.1 follows from the first one since system (5) fulfills the hypotheses of Section 3, too, and the t-fiber of the unstable manifold of (5) is S0 (−t).
The following example shows that not every unstable manifold is a pullback attractor.
8
B. AULBACH, M. RASMUSSEN AND S. SIEGMUND
Example 4.1. An easy calculation shows that the function F : R2 → R defined by F (ξ, η) := − 12 ξ 2 + 12 η 2 + 13 ξ 3 is a Hamiltonian of the planar autonomous system x˙ = y y˙ = x − x2 . This implies that the relation F (λ(t, 0, ξ, η)) = F (ξ, η) holds for all (ξ, η) ∈ R2 and t ∈ R. Then the global unstable manifold ª W u := {x ∈ R2 : lim λ(t, 0, x) = 0
Wu
t→−∞
of the hyperbolic equilibrium (0, 0) is not an attractor of the system: first, for all x ∈ W u F (x) =
lim F (λ(t, 0, x))
t→−∞
= F ( lim λ(t, 0, x)) t→−∞
Figure 1. Homoclinic orbit as a part of the unstable manifold
= F (0) = 0 holds. An analysis of the function F shows then that a homoclinic orbit is part of the unstable manifold W u (see Figure 1). Due to phase plane analysis, it is clear that for every point x in the interior of the homoclinic orbit we have ¡ lim inf d λ(t, 0, {x})|W u ) > 0 . t→∞
This means that the set {x} is not attracted by W u . We consider now a numerical application for Theorem 4.1. Example 4.2. The nonautonomous differential equation x˙ = −x + 0.38 cos(t)(sin(x) + sin(y)) y˙ = y + 0.38 sin(t) sin(x) satisfies the assumptions of Section 3 with the constants n = 2, K = 1, α = −1, β = 1 and L = 0.38. R0
R0 S0
Figure 2. Fiber t = −2
S0
Figure 3. Fiber t = 2
INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS
√ 2 Due to L < 0.382 ≈ β−α 4K 2 (2 + K − 4 + K ) the conditions of Theorem 4.1 are fulfilled. Since we are able to approximate pullback attractors (see Aulbach & Rasmussen & Siegmund [1]), Theorem 4.1 enables us to compute the stable and unstable manifold of this system: Figure 2 and Figure 3 visualize the −2- and 2-fibers in the square [−2, 2] × [−2, 2] of the stable and unstable manifold. In Figure 4 we have plotted the nonautonomous manifolds in the extended phase space in the cuboid [−2, 2] × [−2, 2] × [−π, π]. As one can see, the −π- and the πfibers of the manifolds are identical. This is no coincidence, since in Aulbach & Wanner [3, Corollary 4.4] it is proved that the periodicity of the differential equation carries over to the manifolds.
9
Figure 4. Manifolds in the extended phase space
5. Nonhyperbolic manifolds and pullback attractors. In this section we generalize the results of the previous one by omitting the assumption of hyperbolicity. We use spectral transformations to relate this situation to the hyperbolic case. Theorem 5.1. We suppose that there exists a real constant c with 0 ∈ (c + α + 2KL, c + β − 2KL) , and that the constant L satisfies 0≤L
0. We consider a differential equation x˙ = g(t, x) n
(11)
n
with a function g : R × R → R fulfilling conditions ensuring the global existence and uniqueness of solutions. The transformation y = T (t)x yields the differential equation y˙ = T˙ (t)T −1 (t)y + T (t)g(t, T −1 (t)y) .
(12)
Then (11) has a global pullback attractor if and only if (12) has a global pullback attractor. If A is the global pullback attractor of (11) and A˜ the global pullback attractor of (12), then ˜ = T (t)A(t) A(t)
for all t ∈ R .
INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS
13
Proof. It suffices to prove one of the two implications since (11) is the result of the transformation x = T −1 (t)y of system (12). An obvious calculation shows that the general solution of (12) is given by µ ˜(t, τ, ξ) = T (t)µ(t, τ, T −1 (τ )ξ) for all t, τ ∈ R and ξ ∈ Rn , provided that µ is the general solution of (11). Next we prove that the nonautonomous set A˜ defined by ˜ := T (t)A(t) A(t)
for all t ∈ R
is the global pullback attractor of (12) if (11) has the global pullback attractor A. The invariance of A˜ with respect to system (12) is a consequence of the special transformation and the invariance of A with respect to (11). A˜ is also closed since the linear mapping T (t) : Rn → Rn is a homeomorphism for all t ∈ R. Let K be a generator of A. Since K is bounded, there exists an M1 > 0 with kxk < M1 for ˜ Therefore, we all x ∈ K. We show now that the set UΓM1 (0) is a generator of A. n choose a compact set C ⊂ R arbitrarily. Due to the boundedness of C there exists an M2 > 0 with kxk < M2 for all x ∈ C. Since A is compactly generated (with respect to (11)) there exists a number T (K, UΓM2 (0)) > 0 such that for all τ ∈ R and all t > T (K, UΓM2 (0)) we get µ(τ, τ − t, A(τ − t) ∩ K) ⊃ A(τ ) ∩ UΓM2 (0) . Thus, for all τ ∈ R and all t > T (K, UΓM2 (0)) we have ˜ − t) ∩ UΓM (0)) ⊃ µ µ ˜(τ, τ − t, A(τ ˜(τ, τ − t, T (τ − t)(A(τ − t) ∩ K)) 1 ¡ ¢ = T (τ )µ τ, τ − t, A(τ − t) ∩ K ⊃ T (τ )(A(τ ) ∩ UΓM2 (0)) ˜ )∩C. ⊃ T (τ )A(τ ) ∩ C = A(τ This means that A˜ is compactly generated with respect to (12). Completing this ˜ Since C is bounded proof we show that every compact set C ⊂ Rn is attracted by A. there is an M3 > 0 such that kxk < M3 for all x ∈ C. Then, for all τ ∈ R we get ˜ )) ≤ Γd(T −1 (τ )˜ ˜ )) d(˜ µ(τ, τ − t, C)|A(τ µ(τ, τ − t, C)|T −1 (τ )A(τ ¡ −1 = Γd µ(τ, τ − t, T (τ − t)C)|A(τ )) ¡ ≤ Γd µ(τ, τ − t, UΓM3 (0))|A(τ )) . Letting t tend to infinity the assertion follows and the lemma is proved. For simplicity we refer to system (9) by also writing x˙ = f (t, x). The analogue of Theorem 5.1 is then given as follows. Theorem 6.1. We suppose that there exists a real constant c such that 0 ∈ (c + α + 8K 3 L, c + β − 8K 3 L) and that the constant L satisfies p ¢ β − α¡ L< 2 + 2K 2 − 2 1 + K 4 . 5 32K
(13)
We use the spectral transformation y = ec(t−t0 ) x, t0 ∈ R , to transform (9) into y˙ = cy + ec(t−t0 ) f (t, e−c(t−t0 ) y) .
(14)
14
B. AULBACH, M. RASMUSSEN AND S. SIEGMUND
Then the following connection holds between the global pullback attractor A1 of system (14) and the pseudo-unstable manifold R0 of system (9): A1 (t) = ec(t−t0 ) R0 (t)
for all t ∈ R .
In addition, we have the following connection between the global pullback attractor A2 of the time-reversed version of (14), y˙ = −cy − e−c(t+t0 ) f (−t, ec(t+t0 ) y) , and the pseudo-stable manifold S0 of system (9): A2 (t) = e−c(t+t0 ) S0 (−t)
for all t ∈ R .
√ β−α 2 Proof. From (13) it follows that 2KL < 4(2K − 4 + 4K 4 ). Thus, 2 )2 (2 + 2K Theorem 5.1 can be applied to system (10). Therefore, the global pullback attractor A˜ of −1 ˙ x˙ = S(t)S (t)x + ec(t−t0 ) S(t)f (t, ec(t0 −t) S −1 (t)x) (15) satisfies ˜ = ec(t−t0 ) R ˜ 0 (t) for all t ∈ R . A(t) Since system (14) is obtained by the transformation y = S −1 (t)x from (15), Lemma 6.2 implies ˜ = ec(t−t0 ) S −1 (t)R ˜ 0 (t) = ec(t−t0 ) R0 (t) for all t ∈ R . A1 (t) = S −1 (t)A(t) The corresponding assertion for the pullback attractor A2 can be proved similarly. 6.2. Hierarchies of nonautonomous manifolds. In this section we generalize our situation by considering the system x˙ 1 = B1 (t)x1 + F1 (t, x1 , x2 , . . . , xn ) x˙ 2 = B2 (t)x2 + F2 (t, x1 , x2 , . . . , xn ) .. .
(16)
x˙ n = Bn (t)xn + Fn (t, x1 , x2 , . . . , xn ) , Ni ×Ni
where Bi : R → R , Fi : R×RN → RNi (i = 1, . . . , n) are continuous mappings with Fi (t, 0, 0) = 0 for all t ∈ R (N = N1 + · · · + Nn ). Moreover we assume: (H1) Hypotheses on linear part: There exist real constants K ≥ 1 and αi < βi (i = 1, . . . , n − 1) with βi ≤ αi+1 (i = 1, . . . , n − 2) such that the evolution operators Φi of the linear equations x˙ i = Bi (t)xi (i = 1, . . . , n) satisfy the estimates kΦi (t, s)k ≤ Keαi (t−s) kΦi+1 (t, s)k ≤ Ke
βi (t−s)
for all t ≥ s , for all t ≤ s
for all i = 1, . . . , n − 1. (H2) Hypotheses on perturbation: For all x = (x1 , . . . , xn ), x ¯ = (¯ x1 , . . . , x ¯n ) ∈ RN , t ∈ R and i ∈ {1, . . . , n} we have n X kFi (t, x1 , . . . , xn ) − Fi (t, x ¯1 , . . . , x ¯n )k ≤ L kxj − x ¯j k = Lkx − x ¯k , j=1
where the constant L satisfies the estimate ¾ ½ βi − αi : i = 1, . . . , n − 1 . 0 ≤ L < min 4K(n − 1)
INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS
15
We denote the general solution of (16) by λ and choose for i = 1, . . . , n − 1 i constants δi ∈ (2KL(n − 1), βi −α 2 ]. Theorem 6.2. There exist nonautonomous manifolds Wi,j (1 ≤ i ≤ j ≤ n), the socalled hierarchy of nonautonomous manifolds, with the following characterizations: • W1,n = R × RN . • For all i ∈ {1, . . . , n − 1} and γ ∈ [αi + δi , βi − δi ] we have W1,i = {(τ, ξ) : λ(·, τ, ξ) is γ + -quasibounded} , Wi+1,n = {(τ, ξ) : λ(·, τ, ξ) is γ − -quasibounded} . • For all 1 < i ≤ j < n we have Wi,j = W1,j ∩ Wi,n . Hence, for every choice of γ1 ∈ [αj + δj , βj − δj ] and γ2 ∈ [αi−1 + δi−1 , βi−1 − δi−1 ] the relation Wi,j = {(τ, ξ) : λ(·, τ, ξ) is γ1+ - and γ2− -quasibounded} is fulfilled. Proof. See Aulbach & Wanner [3, Theorem 5.1]. The following diagram visualizes the relations between the manifolds Wi,j of the hierarchy: W1,1
⊂
W1,2 ∪ W2,2
⊂ ···
⊂
⊂ ···
⊂
..
W1,n−1 ∪ W2,n−1 ∪ .. .
⊂
Wn−1,n−1
⊂
.
W1,n ∪ W2,n ∪ .. .
⊂
Wn−1,n ∪ Wn,n
We state now the version of Theorem 5.1 which applies to system (16). For simplicity we refer to (16) by also writing x˙ = f (t, x) . Theorem 6.3. We suppose that there exists an i ∈ {1, .., n − 1} and a c > 0 with 0 ∈ (c + αi + 2KL(n − 1), c + βi − 2KL(n − 1)) , and that the constant L satisfies p βi − αi 0≤L< (2 + K − 4 + K 2 ) . 2 4K (n − 1) We use the spectral transformation y = ec(t−t0 ) x, t0 ∈ R, to transform (16) into y˙ = cy + ec(t−t0 ) f (t, e−c(t−t0 ) y) .
(17)
Then the following connection holds between the global pullback attractor A1 of system (17) and the manifold Wi+1,n of system (16): A1 (t) = ec(t−t0 ) Wi+1,n (t)
for all t ∈ R .
In addition, we have the following connection between the global pullback attractor A2 of the time-reversed version of (17), y˙ = −cy − e−c(t+t0 ) f (−t, ee(t+t0 ) y) , and the manifold W1,i of system (16): A2 (t) = e−c(t+t0 ) W1,i (−t)
for all t ∈ R .
16
B. AULBACH, M. RASMUSSEN AND S. SIEGMUND
Proof. This result follows directly from Theorem 5.1 by putting the components 1, . . . , i and i + 1, . . . , n together. Remark 6.1. Provided the constant L of system (16) is small enough, every manifold of the hierarchy can be approximated numerically. For the manifolds W1,i and Wi+1,n , i = 1, . . . , n − 1, this follows from Theorem 6.3. The manifolds Wi,j for 1 < i ≤ j < n can be obtained by intersecting W1,j and Wi,n (see Theorem 6.2). Remark 6.2. In Section 6.1 the coupled linear part of system (9) was “decoupled” by means of the assumed exponential dichotomy. Using the Sacker-Sell or dichotomy spectrum (see [14]) we can apply the results in [15] to transform x˙ = B(t)x + F (t, x) to (16) with decoupled linear part where each component corresponds to a spectral interval of x˙ = B(t)x. 7. Approximation of a strong stable manifold of the Lorenz system. In 1963, the meteorologist Lorenz considered the three dimensional autonomous system x˙ = σ(y − x) y˙ = ρx − y − xz z˙ = −βz + xy . For values and ρ > 1 this system has the r1 := (0, 0, 0), p β > 0p p three equilibria p r2 := ( β(ρ − 1), β(ρ − 1), ρ−1) and r3 := (− β(ρ − 1), − β(ρ − 1), ρ−1) (see Guckenheimer & Holmes [9, Chapter 2.3]). The linearization in the equilibrium r1 is given by −σ σ 0 D = ρ −1 0 . 0 0 −β For the parameter values suggested by Lorenz (σ = 10, ρ = 28 and β = 38 ) the matrix D has the eigenvalues λ1 ≈ −22.8, λ2 ≈ −2.7 and λ3 ≈ 11.8. Thus, the stable manifold of the origin has dimension two. Because of λ1 < λ2 , there also exists a strong stable manifold corresponding to the eigenvalue λ1 . In order to compute the stable and the strong stable manifold in the cuboid [−25, 25]3 we used the subdivision algorithm for the approximation of the pullback attractor of the time-reversed as well as the time-reversed and spectral-transformed system (c = 3). Since the Lorenz system does not fulfill the global assumptions of our theorems it is not (rigorously) guaranteed that the outcome of our computations – shown in Figure 5 – precisely represent the manifolds under consideration. However, using the ideas elaborated in this paper also local versions of our theorems can be derived. It is worth mentioning that – with our approach – the approximation of the strong stable manifold is not possible in a purely autonomous setting. This is due to the fact that the spectral transformation always yields a nonautonomous system. Acknowledgements. The computations have been carried out with an extended version of the computer program GAIO (Global Analysis of Invariant Objects). This software package is being developed since 1995 by Michael Dellnitz and Oliver Junge. The most important commands as well as some implemented algorithms can be found in Dellnitz & Froyland & Junge [8]. The two dimensional graphics in this paper have been created with the computer program Matlab. For the three dimensional pictures we used Grape from the research group ”Nichtlineare Partielle Differentialgleichungen”’ of the University of Bonn.
INVARIANT MANIFOLDS AS PULLBACK ATTRACTORS
17
Figure 5. Stable and strong stable manifold of the origin from two different perspectives REFERENCES [1] B. Aulbach, M. Rasmussen, S. Siegmund, Approximation of Attractors of Nonautonomous Dynamical Systems, to appear in: Discrete and Continuous Dynamical Systems – Series B. [2] B. Aulbach, M. Rasmussen, S. Siegmund, Invariant Manifolds and Pullback Attractors of Nonautonomous Difference Equations, to appear in: Proceedings of ICDEA 2003. eodory Type Differential Equations in [3] B. Aulbach, T. Wanner, Integral Manifolds for Carath´ Banach Spaces, in: B. Aulbach & F. Colonius (eds.), Six Lectures on Dynamical Systems, World Scientific, Singapore (1996). eodory Type Differential Equations [4] B. Aulbach, T. Wanner, Invariant Foliations for Carath´ in Banach Spaces, in: V. Lakshmikantham & A. A. Martynyuk (eds.), Advances of Stability Theory at the End of the 20th Century, Taylor & Francis, London (2003). [5] W. A. Coppel, Dichotomies in Stability Theory (Springer Lecture Notes in Mathematics 629). Springer, Berlin Heidelberg New York (1978). [6] M. Dellnitz, A. Hohmann, The Computation of Unstable Manifolds Using Subdivision and Continuation, in: Progress in Nonlinear Differential Equations and Their Applications 19, Birkh¨ auser (1996), 449–459. [7] M. Dellnitz, A. Hohmann, A Subdivision Algorithm for the Computation of Unstable Manifolds and Global Attractors, Numerische Mathematik 75 (1997), 293–317.
18
B. AULBACH, M. RASMUSSEN AND S. SIEGMUND
[8] M. Dellnitz, G. Froyland, O. Junge, The Algorithms behind GAIO – Set Oriented Numerical Methods for Dynamical Systems, in: B. Fiedler (ed.): Ergodic Theory, Analysis, and Efficient Simulation of Dynamical Systems, Springer (2001). [9] J. Guckenheimer, P. Holmes, Nonlinear Oscillation, Dynamical Systems, and Bifurcations of Vector Fields, Springer, New York (1983). [10] J. Hadamard, Sur l’it´ eration et les solutions asymptotiques des ´ equations diff´ erentielles, Bull. Soc. Math. France 29 (1901), 224–228. [11] P. E. Kloeden, H. Keller, B. Schmalfuß, Towards a Theory of Random Numerical Dynamics, in: H. Crauel & M. Gundlach (eds.), Stochastic Dynamics, Springer, Berlin Heidelberg New York (1999). e, M´ emoire sur les courbes d´ efinie par une ´ equation diff´ erentielle, IV, J. Math. [12] H. Poincar´ Pures Appl. 2 (1886), 151–217. [13] M. Rasmussen, Approximation von Attraktoren und Mannigfaltigkeiten nichtautonomer Systeme, Diploma Thesis, University of Augsburg (2002). [14] S. Siegmund, Dichotomy Spectrum for Non-autonomous Differential Equations, Journal of Dynamics and Differential Equations 14 (2002), 243–258. [15] S. Siegmund, Reducibility of Nonautonomous Linear Differential Equations, Journal of the London Mathematical Society 65 (2002), 397–410. [16] S. Siegmund, Computation of Nonautonomous Invariant Manifolds, in: B. Aulbach & S. Elaydi & G. Ladas (eds.), New Progress in Difference Equations, CRC-Press, Boca Raton (2004). E-mail address:
[email protected] E-mail address:
[email protected] E-mail address:
[email protected]