Normal forms and invariant manifolds for nonlinear, non-autonomous PDEs, viewed as ODEs in infinite dimensions
Peter Hochs, A.J. Roberts

TL;DR
This paper develops a method to transform nonlinear, non-autonomous PDEs into a normal form in infinite-dimensional spaces, enabling clearer analysis of invariant manifolds and approximate solutions.
Contribution
It introduces a normal form approach for nonlinear, non-autonomous PDEs viewed as ODEs in Fréchet spaces, facilitating the study of invariant and approximate center manifolds.
Findings
Normal form transformation separates center, stable, and unstable coordinates.
Method applies to a class of nonlinear, non-autonomous PDEs.
Enables analysis of approximate center manifolds with controlled precision.
Abstract
We prove that a general class of nonlinear, non-autonomous ODEs in Fr\'echet spaces are close to ODEs in a specific normal form, where closeness means that solutions of the normal form ODE satisfy the original ODE up to a residual that vanishes up to any desired order. In this normal form, the centre, stable and unstable coordinates of the ODE are clearly separated, which allows us to define invariant manifolds of such equations in a robust way. In particular, our method empowers us to study approximate centre manifolds, given by solutions of ODEs that are central up to a desired, possibly nonzero precision. The main motivation is the case where the Fr\'echet space in question is a suitable function space, and the maps involved in an ODE in this space are defined in terms of derivatives of the functions, so that the infinite-dimensional ODE is a finite-dimensional PDE. We show that our…
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsStability and Controllability of Differential Equations · Numerical methods for differential equations · Advanced Differential Equations and Dynamical Systems
\includeversion
ArXiv \excludeversionjournal
{ArXiv}
{journal}
{ArXiv}
{journal}\address
School of Mathematical Sciences, University of Adelaide \[email protected]
\addressSchool of Mathematical Sciences, University of Adelaide \emailmailto:[email protected] \urladdrhttp://orcid.org/0000-0001-8930-1552
Normal forms and invariant manifolds for nonlinear, non-autonomous PDEs, viewed as ODEs in infinite dimensions
Peter Hochs School of Mathematical Sciences, University of Adelaide, [email protected]
A.J. Roberts School of Mathematical Sciences, University of Adelaide, mailto:[email protected], http://orcid.org/0000-0001-8930-1552
Peter Hochs
A.J. Roberts
Normal forms and invariant manifolds for nonlinear, non-autonomous PDEs, viewed as ODEs in infinite dimensions
Peter Hochs School of Mathematical Sciences, University of Adelaide, [email protected]
A.J. Roberts School of Mathematical Sciences, University of Adelaide, mailto:[email protected], http://orcid.org/0000-0001-8930-1552
Peter Hochs
A.J. Roberts
Abstract
We prove that a general class of nonlinear, non-autonomous odes in Fréchet spaces are close to odes in a specific normal form, where closeness means that solutions of the normal form ode satisfy the original ode up to a residual that vanishes up to any desired order. In this normal form, the centre, stable and unstable coordinates of the ode are clearly separated, which allows us to define invariant manifolds of such equations in a robust way. In particular, our method empowers us to study approximate centre manifolds, given by solutions of odes that are central up to a desired, possibly nonzero precision. The main motivation is the case where the Fréchet space in question is a suitable function space, and the maps involved in an ode in this space are defined in terms of derivatives of the functions, so that the infinite-dimensional ode is a finite-dimensional pde. We show that our methods apply to a relevant class of nonlinear, non-autonomous pdes in this way.
Contents
-
5 Polynomial and differentiable maps on graded Fréchet spaces
-
5.3 Sequences of Banach spaces comparable to sequences of Hilbert spaces
1 Introduction
1.1 Background and motivation
Various invariant manifolds are central to many areas of dynamical systems, including using centre manifolds to construct and justify reduced low-dimensional models of high-dimensional dynamics [29, e.g.]. Many dynamical systems involve pdes in infinite dimensional state spaces of functions, and some applications require infinite dimensional centre manifolds [7, 28, e.g.]. In general we also want to cater for non-autonomous systems, with an aim to subsequently generalise to stochastic/rough dynamics [14, e.g.]. Further, encompassing unstable dynamics with both centre and stable is necessary for application to Saint Venant-like, cylindrical, problems [18, 19, 20, e.g.], and to deriving boundary conditions for approximate pdes [27, e.g.]. Consequently, here we address the general challenge of constructing and justifying various infinite dimensional invariant manifolds for non-autonomous dynamical systems which have stable, unstable and centre modes. A crucial novel feature of the approach is that we further develop a backward theory recently initiated for finite dimensional systems [30]: analogous backward theory has been very useful in other domains [16, e.g.].
Applications of the extant forward theory in such a general setting is often confounded by impractical preconditions. Non-autonomous invariant manifold theories typically require bounded operators, and Lipschitz and/or uniformly bounded nonlinearities, [3, 4, 5, 10, 17, e.g.]. The extant boundedness requirement [17, Hypothesis 2.1(i) and 3.8(i), e.g.] arises from the general necessity of both forward and backward time convolutions with the semigroup (e.g., for systems that linearise to ), convolutions that must be continuous in extant forward theory, but cannot be continuous with unbounded operators. Despite many interesting specific scenarios having rigorous invariant manifolds established via strongly continuous semigroup operators and by mollifying nonlinearity [9, 32, e.g.], extant non-autonomous forward theory fails to rigorously apply in many practical cases.
Our main motivation for studying odes in infinite-dimensional vector spaces is their possible application to analysing invariant manifolds of pdes in finite space-time dimensions. In that setting, the infinite-dimensional vector space in question is a space of functions, and the maps occurring in the ode are differential operators. The linear part of such an ode is a linear partial differential operator, which typically is unbounded in applications. Such an operator can be viewed as a bounded operator between different Banach spaces, with norms adapted to make the operator bounded. (For example, the operator is unbounded on , but becomes bounded if we take its domain to be a first-order Sobolev space.) Centre manifold theory in this setting was developed by several authors [17, 21, 33, e.g.], and applied to pdes.
However, to achieve our goal of developing the desired backward theory, and robustly constructing invariant manifolds via coordinate transformations to approximate normal forms, we need to go beyond this setting. This essentially boils down to the fact that a bounded operator on a single Banach space can be iterated to yield new bounded operators, whereas this is not possible for a bounded operator between different Banach spaces. The necessity of iterating operators in our constructions leads us outside the setting of Banach spaces, to graded Fréchet spaces: intersections of infinite sequences of Banach spaces connected by bounded inclusion maps. These include spaces relevant to the study of pdes, such as spaces of smooth functions.
1.2 Results
The first step in the proposed backward theory is to establish an approximate conjugacy between a given system and a ‘nearby’ system for which we know its invariant manifolds, by its construction. Figure 1 illustrates what this article achieves. Planned future research will then provide novel finite domain and error bounds as illustrated in Figure 1. That is, instead of proving that there exists a reduced dimensional manifold for a specified system, which is then approximately constructed, our main results, Theorems 2.17 and 2.22, establish that there is a system which is both ‘close’ to the specified system, and also has a reduced dimensional manifold which we know exactly. In essence, we invoke an (extended) normal form coordinate transform—related to Hartman–Grobman theory [3, 4, 5, e.g.]—and use it from a new point of view.
An intuitive formulation of our main result on the existence of such a normal form is the following Theorem 1.1. Theorem 2.17 is a precise formulation, and Corollary 2.22 is a special case that applies to nonlinear, non-autonomous pdes. (In that special case, the space is a space of functions, and the maps that occur are partial differential operators.)
Theorem 1.1** (Normal form theorem, intuitive formulation).**
Let be an interval in , and a possibly infinite-dimensional vector space. Consider a nonlinear, non-autonomous ode
[TABLE]
for , where is a linear operator on (independent of ), and and its derivative vanish on . For each , there are both an ode
[TABLE]
such that
- •
if satisfies 1.2a, then setting defines a solution of 1.1 up to a residual term that vanishes to order ;
- •
the map is of order in its second entry ;
- •
*the component of a solution to 1.2a in the stable subspace for decays exponentially quickly to zero as increases in . Its component in the unstable subspace for decays exponentially quickly to zero as decreases in . If the solution starts out in either the centre-stable or the centre-unstable subspace for , then its component in the central subspace for is bounded by a constant for all , or at worst by a specified, small exponential growth rate. *
The three most important things to make more precise in this intuitive formulation is what ‘order ’ means, the related question what topology on is used, and what kinds of maps , , and are.
The last point in Theorem 1.1 means that the centre, stable and unstable manifolds (in this case, linear subspaces) of 1.2a are exactly the centre, stable and unstable spaces of , respectively. (And similarly for the centre-stable and centre-unstable subspaces.) The centre, stable, unstable, centre-stable and centre-unstable subspaces for the dynamics in described by 1.2a and (which becomes an ode in if is invertible for all ) are then obtained from these spaces via an application of the coordinate transform . In this way, we show that any (non-autonomous) system of the form 1.1 is arbitrarily close to a system with robustly defined invariant manifolds.
This definition of these key invariant manifolds is a crucial reformation of the backward theory proposed. Classic definitions of un/stable and centre manifolds require the existence of limits as time goes to [2, 6, 17, 26, e.g.]. This consequently requires solutions of the dynamical system to be well-behaved for all time, which requires constraints that in applications are often not found, or are hard to establish. For example, in stochastic systems very rare events will eventually happen over the infinite time requiring global Lipschitz and boundedness that are oppressive in applications. By modifying definitions we establish results for finite times, which are useful in many applications, and for a wider range of non-autonomous systems.
1.3 Ingredients of the proof
The key ingredients of the proofs of our main results are nested sequences of Banach spaces, whose intersections are graded Fréchet spaces; compact polynomial maps between Banach spaces and graded Fréchet spaces; and compactly differentiable maps between such spaces.
1.3.1 Sequences of Banach spaces
It is important to specify what is meant by a nearby infinite-dimensional mathematical system in Figure 1. Intuitively, we mean by this that solutions of the nearby system 1.2 are solutions of the original system 1.1 up to any desired order in the magnitude of such a solution, as in Theorem 1.1. To be more precise about what that this magnitude is, we need to specify norms or seminorms on spaces containing these solutions. Working with a single Banach space (i.e., a single norm) is too restrictive for applications. This is because in applications to pdes, the maps and in Theorem 1.1 are generally not continuous maps from a single Banach space to itself. This could be remedied by allowing maps between two different Banach spaces, but that would not let us iteratively apply maps involving and , which we do in the proof of Theorem 1.1.
A type of space that is both general enough to apply to various nonlinear pdes, while being close enough to Banach spaces to allow us to define a meaningful notion of a solution of an equation up to a given order, is what are often called graded Fréchet spaces. These are intersections of sequences of Banach spaces, each with a bounded inclusion map into the next. The notion of an operator of a given order on a graded Fréchet space is then defined in terms of the norms on these Banach spaces, see Definition 2.4.
For several convergence questions, it would be useful if the Banach spaces that occur in the definition of a graded Fréchet space are Hilbert spaces. Then we can use orthogonality, for example. However, for applications to nonlinear pdes, it is not enough to use Hilbert spaces. For example, a nonlinear term is a well-defined (and differentiable) map from the Banach space to the Hilbert space . To be able to use Hilbert space techniques in such settings, we use the notion of nested sequences of Banach space that are comparable to nested sequences of Hilbert spaces (see Definition 2.6). This effectively means that a graded Fréchet space that naturally occurs as an intersection of Banach spaces can equivalently be presented as an intersection of Hilbert spaces. Proving such a property in situations relevant to pdes involves the relevant Sobolev embedding theorems.
1.3.2 Compact polynomial maps
We construct coordinate transformations on graded Fréchet spaces to bring odes in such spaces into normal forms that allow us to define invariant manifolds directly and robustly. These transformations are polynomial maps, which we construct by adding (infinitely many) monomial terms with the right properties together. At the level of Banach spaces, a polynomial map can be naturally defined as a finite sum of restrictions to the diagonal of bounded multilinear maps. For example, Taylor polynomials of differentiable maps between normed vector spaces are polynomials of this type.
But not all such polynomial maps (for example, the identity map) can be approximated by sums of monomials. This leads us to define compact polynomial maps, which can be approximated in this way in settings relevant to us. The notion of a compact polynomial map that we use seems natural, but we have not been able to find it elsewhere in the literature. Different notions of compact polynomial maps were developed and studied by Gonzalo, Jaramillo and Pełczyńsky [15, 25].
1.3.3 Compactly differentiable maps
Our construction of the required coordinate transforms involves Taylor polynomials of differentiable maps between Banach spaces, and between graded Fréchet spaces. This construction is possible if those polynomials are compact in the sense just mentioned. That is the case for compactly differentiable maps, which we define for this purpose. We will see in Section 4.4 that, in applications to pdes, the relevant differentiable maps are indeed compactly differentiable. This follows from various Sobolev embedding theorems.
1.4 Outline of this paper
The main results of this paper, on normal forms and invariant manifolds of nonlinear, non-autonomous odes in Fréchet spaces, and of nonlinear, non-autonomous pdes in finite-dimensional spaces, are stated in Section 2. We illustrate our results by applying them to an example pde in Section 9.
In the rest of the paper, we prove our main results. We start by reviewing standard material on differentiable maps and polynomials on normed spaces in Section 3. In Sections 4 and 5, we develop technical tools we need for our proofs. Then in Sections 6 and 7, we use these tools to prove the main Theorems 2.17 and 2.21. We prove some properties of the normal form equation, which allow us to identify its invariant manifolds, in Section 8.
A key ingredient in the proof of a version Taylor’s theorem for compactly differentiable maps between Fréchet spaces, mentioned above, is the fact that a compact operator from a Banach space with the approximation property into another Banach space can be approximated by finite-rank operators in a suitable way. {ArXiv} This is reviewed in Appendix A.
1.5 Notation and conventions
We write for the set of positive integers, and for the set of nonnegative integers. We write for the set of sequences in with finitely many nonzero entries, interpreted as multi-indices. For , or in , we denote the (finite) sum of its elements by .
We denote spaces of bounded linear operators by the letter , and spaces of compact linear operators by the letter .
When we mention a normed vector space , the implicitly given norm is denoted by . Similarly, if is an inner product space, then the inner product is denoted by . Inner products on complex vector spaces are assumed to be linear in their second entries, and antilinear in their first entries. For maps , when we write , we implicitly mean that as in .
If is a normed vector space, and is an open interval, and and , for , are maps, then we say that converges to if converges to in uniformly in in compact subsets of . If and are smooth, then we say that converges to differentiably in if converges to for every , in this sense.
For maps and , the maps and are defined by
[TABLE]
for and .
If is an open subset of and , then the Sobolev space of functions from to with weak derivatives up to order in is denoted by . The norm on this space is
[TABLE]
If , we write .
2 Preliminaries and results
Our main result, Theorem 2.17, asserts that a broad class of nonlinear pdes and odes in infinite-dimensional vector spaces may be effectively approximated by normal form systems via well-chosen, time-dependent, coordinate transformations. In this normal form, the centre, stable and unstable components of the pde and ode are clearly separated, which allows us to define centre manifolds for this class of equations in a robust way (Definition 2.10).
We first state our result on normal forms, and the definition of centre manifolds, for odes in a class of abstract vector spaces (Section 2.5). Our main reason for developing this theory is to apply it to the study of pdes, for which the vector spaces used are spaces of functions, and the relevant maps between them are defined in terms of derivatives of functions. We discuss a relevant class of examples of such function spaces and maps in Section 2.6.
2.1 Nested sequences of Banach spaces
The normal form we obtain in Theorem 2.17 is approximate in the sense that functions satisfying an equation transformed into that form satisfy the original equation up to a residual term. An important point in Theorem 2.17 is that this residual vanishes up to a specified order. To make it precise what this vanishing up to a certain order means, we introduce the type of topological vector spaces we consider in this subsection. More details about these spaces and their properties are given in Section 5. A concrete class of examples of these spaces relevant to the study of pdes is given in Section 2.6.
Definition 2.1**.**
By a nested sequence of Banach spaces, we mean a sequence of Banach spaces such that
- •
for every , , where the inclusion map is bounded, and
- •
the intersection is dense in for every .
We then consider as a Fréchet space111Much of what we write about Fréchet spaces of this form holds for more general projective limits of Banach spaces connected by bounded operators. But we do not need that degree of generality. with the seminorms (now actual norms) that are the restrictions of the norms on the spaces .
A compactly nested sequence of Banach spaces is such a sequence such that for every , there is a such that the inclusion is compact.
A Fréchet space as in this Definition 2.1 is often called a graded Fréchet space.
Let be a nested sequence of Banach spaces.
Definition 2.2**.**
The space of bounded operators on consists of the linear maps such that for every , there is a such that the linear map extends continuously to a map in .
Remark 2.3**.**
In Definition 2.2, if , then the composition
[TABLE]
is a bounded operator on . So we may always take in this context, but this does not need to be assumed a priori. Similar remarks apply in analogous situations, such as Definitions 2.4 and 2.5 below.
Definition 2.4**.**
A map is of order , written as , if for every , there is a such that as in .
If is an open interval, a map is of order , written as , if for every , there is a such that as in , uniformly in in compact subsets of .
Definition 2.5**.**
An times differentiable map from to itself is a map such that for every , there is a such that extends to an times differentiable map from to . If a map is times differentiable for every , then it is infinitely differentiable.
Basic material on differentiable maps between normed vector spaces is reviewed in Section 3.
Definition 2.6**.**
Two nested sequences and of Banach spaces are comparable if for every , there are such that we have bounded inclusions and .
In the setting of this definition, .
2.2 Setup and goal
Let be a compactly nested sequence of Banach spaces, such that is a Hilbert space. Let . Let be an open interval in containing , and let be infinitely differentiable with respect to and . Suppose that , and that for every , there is a such that is differentiable, and
[TABLE]
uniformly in in compact subsets of .
Suppose that is a set of eigenvectors of which is a Hilbert basis of (an orthonormal set that spans a dense subspace of ). We assume that the sequence is comparable to a nested sequence of separable Hilbert spaces in which the vectors are orthogonal. However, we will see in Remark 6.1 that we may equivalently make the seemingly stronger but more concrete assumption that the spaces themselves are separable Hilbert spaces. The assumption that is comparable to a nested sequence of separable Hilbert spaces is easier to check in practice than the condition that every space can be chosen to be a separable Hilbert space itself. For example, Section 2.6 discusses a class of relevant cases where the spaces are not Hilbert spaces for . In this sense, the notion of comparable sequences of Banach spaces is a tool that makes it easier to check the conditions of Theorem 2.17.
We study smooth maps satisfying the non-autonomous dynamical system differential equation
[TABLE]
Since the nonlinearity satisfies 2.1, is an equilibrium of the system 2.2. We provide a novel backward approach to establish invariant manifolds in a finite domain about the equilibrium . For these invariant manifolds to be useful in applications, the time interval will be long enough for transient dynamics to decay to insignificance in the context of the application. The proofs of our main results simplify considerably if the time interval is short, or bounded. But we emphasise that we only aim this theory to support the many applications where the time interval is long enough, or unbounded, so that the theorems are non-trivially useful in the application.
2.3 Dynamics in a normal form
We define invariant manifolds, or sets, for dynamical systems in a particular normal form, and show that this definition captures the essence of such manifolds. In Section 2.5, we show that a very general class of odes of the form 2.2 can be brought into this normal form, modulo residuals that vanish to a desired order.
Spectral gap in an exponential trichotomy
Let be such that , and no eigenvalues of have real parts in the intervals and (depending upon the circumstances, or could be , and/or may be zero). For every , let be the eigenvalue of corresponding to . With respect to the parameters , and , we define the sets of indices of central, stable and unstable eigenvalues and eigenvectors, respectively, as
[TABLE]
For , let be the closure in of the span of the eigenvectors , for . For any map into and , we write for its components in . The sets are respectively called the centre/stable/unstable subspaces. Further, we define the centre-stable subspace , and the centre-unstable subspace .
For and a multi-index , we set
[TABLE]
For , write , for such that if , if and if .
Normal form dynamics
Definition 2.7**.**
A smooth map separates invariant subspaces if the components of in , and are of the forms
[TABLE]
for all and , for smooth maps , where the series converge in , differentiably in .
Consider a polynomial map that separates invariant subspaces, and the ode
[TABLE]
in smooth maps . Because separates invariant subspaces, this ode has very explicit invariant manifolds, by Lemma 2.8 and Proposition 2.9 below.
Lemma 2.8**.**
Suppose that satisfies 2.4, where separates invariant subspaces. Let . If there exists a such that , then for all .
Proposition 2.9**.**
There is a neighbourhood of in , with the following property. Let be a solution of 2.4, for some open interval containing [math], and where sepearates invariant subspaces. Write , with for .
- •
If for all with , then for all such , .
- •
If for all with , then for all such , .
- •
Suppose that or . If for all , then for all , .
Since and are positive, this proposition in particular states that stable solutions decrease to zero exponentially quickly as increases in , while unstable solutions decrease to zero exponentially quickly as decreases in . The numbers and represent bounds on what one takes to be relatively small real parts of eigenvalues of (classically, these numbers are zero), so that the third point in Proposition 2.9 intuitively states that central solutions, at worst, only grow relatively slowly as increases.
Lemmas 2.8 and 2.9 are proved in Section 8. The specific form of the set is also specified there, see 8.7.
2.4 Invariant manifolds
Lemmas 2.8 and 2.9 show that, for every , the set
[TABLE]
is a centre, stable or unstable submanifold of for 2.4, respectively. Furthermore, for and , we obtain centre-stable and centre-unstable manifolds, respectively. (Here we use the cases of the third point in Proposition 2.9 where and , respectively.) This motivates Definition 2.10 of invariant subspaces of dynamical systems of a certain form. To state it precisely, we incorporate existence of solutions of 2.4.
For , we write for the infimum of the set of all such that there is a solution of 2.4, with . Similarly, is the supremum of the set of all such that there is a solution of 2.4, with . If such and exist, we set . (In particular, if such a solution exists for all .) If such an exists but no , we set , and if such a exists but no , we set . If there are no such , we set .
Definition 2.10**.**
Let be a smooth map, and let be a polynomial map that separates invariant subspaces. Consider the dynamical system for smooth maps determined by
[TABLE]
for , for a smooth map satisfying 2.4. Let be as in Proposition 2.9. For every , set
[TABLE]
The set is a centre subset of the dynamical system in ; the set is a stable subset of the system; and the set is an unstable subset of the system. The set is a centre-stable subset, and is a centre-unstable subset of the system. Such spaces are invariant or integral subsets of the dynamical system in .
Remark 2.11**.**
If the map in Definition 2.10 is invertible for all (on a suitable domain), then the dynamical system in in that definition is equivalent to the ode
[TABLE]
Remark 2.12**.**
In general, existence and uniqueness of solutions of 2.4 is not guaranteed, hence the careful definition of . Existence and uniqueness of solutions is an assumption in previous definitions [17, Theorem 2.9, e.g.]; see Hypothesis 2.7 in that reference. There are existence and uniqueness results if satisfies a local Lipschitz condition, but that is not the case in many applications to pdes. Under additional assumptions, Vanderbauwhede & Iooss [33, proof of Theorem 3] showed such a local Lipschitz condition holds.
Remark 2.13**.**
Invariant subsets or submanifolds are not unique in general; here this non-uniqueness is due to various possibilities for , and , and is reflected in the use of the indefinite article in Definition 2.10.
Example 2.14**.**
For one example of the non-uniqueness engendered via , consider the classic example system of and in the role of 2.2 (and let the step function when , and when ). This ode system may be given, for every , as the coordinate transformation, 2.5, and together with the system, 2.4, and (by design, here symbolically identical to the original -system). Lemma 2.8 identifies as the centre subspace of this -system. Definition 2.10 then gives the classic non-uniqueness that, for every , are centre manifolds for the -system.
Remark 2.15**.**
In the setting of Definition 2.10, if is a local diffeomorphism in the Fréchet manifold sense, then the subsets in Definition 2.10 are Fréchet manifolds. This would justify the more specific terminology invariant submanifolds rather than just invariant subsets.
2.5 Main result: an approximate normal form
Our main result, Theorem 2.17, states that for an ode of the form 2.2, there is a dynamical system in the normal form used to define invariant manifolds in Definition 2.10, such that solutions of the normal form system satisfy 2.2 up to a residual term that vanishes to any desired order. In this sense, 2.2 is arbitrarily close to a dynamical system with clearly and robustly defined invariant manifolds.
Definition 2.16**.**
A function grows at most polynomially if there are such that for all , . (This condition holds for all bounded functions if is bounded.)
An infinitely differentiable map has polynomial growth if for every , every , and every , the function
[TABLE]
grows at most polynomially.
We use the term -regular integral for an integral of the form
[TABLE]
where and grows at most polynomially. The larger , the better the convergence properties of -regular integrals.
Theorem 2.17**.**
Let be such that , and . Suppose that has polynomial growth. Then there are three infinitely differentiable maps , such that
- •
* and separates invariant manifolds;*
- •
,
and if a smooth map is given by
[TABLE]
for all , for a smooth map satisfying
[TABLE]
for all , then for all ,
[TABLE]
Finally, there is a construction of the map in which all integrals over that occur are -regular.
We prove this theorem in Sections 3, 4, 5, 6 and 7; see in particular Section 7.3.
Theorem 7.9 shows that the maps and can be chosen to be polynomials of a certain type. The conclusion that the integrals occurring are -regular is more than just convenient: this is clear in the classical case where (see Remark 2.18).
Remark 2.18**.**
In cases where the centre eigenvalue bound equals zero, we can always choose so small that the conditions on and in Theorem 2.17 are satisfied. In these cases, the residual can be made to vanish to arbitrarily large order . Furthermore, the integrals that occur in the construction of (see Definition 7.2) are -regular for some precisely if they converge. Hence -regularity for some is a necessary condition for the construction to make sense.
Choosing the centre eigenvalue bound positive, which imposes a positive lower bound on , restricts the vanishing order of , but also makes the construction of the coordinate transform more robust, in the sense that the integrals over in its construction are -regular. Many researchers choose to phrase problems as singular perturbations [8, 24, 34, e.g.]. In such cases the bounds on the hyperbolic rates as the perturbation parameter . Consequently, choosing (say) then the residual again can be made to vanish to arbitrarily large order for small enough .
However, in applications we generally require an invariant manifold in some chosen domain of interest that resolve phenomena on chosen time scales of interest. Such subjective choices, informed by the governing equations, generally dictate the chosen bound separated by a big enough gap from the bounds so that the centre manifold evolution, constructed to a valid order , provides a useful model over the chosen domain for the desired phenomena.
Remark 2.19**.**
The derivative at of the coordinate transformation is the identity map, and hence invertible. If a suitable generalisation of the inverse function theorem applies to , such as a version of the Nash–Moser theorem, then it follows that is a local diffeomorphism at zero. Then it would be justified to call the invariant subsets of Definition 2.10 invariant submanifolds in this setting (at least in a neighbourhood of zero), see Remark 2.15.
Remark 2.20**.**
In the proof of Theorem 2.17, explicit constructions of the maps and are given. In practice, however, it can be easier to determine these maps in more direct ways. This is illustrated in an example in Section 9. Theorem 2.17 implies that one can always find these maps. We prove this by giving a construction that always leads to an answer, even though more direct constructions may exist in specific situations.
Similarly, the domain in Proposition 2.9, defined in 8.7, is guaranteed to have the properties in Proposition 2.9. In practice, these properties often hold on much larger domains.
2.6 A general class of PDEs in bounded domains
Because Theorem 2.17 applies to abstract Banach spaces , it gives one the flexibility to choose these spaces such that, for specific pde applications,
the residual is of order with respect to norms relevant to the problem, and 2. 2.
the spaces incorporate the relevant boundary conditions.
This subsection explores a class of nonlinear pdes to which Theorem 2.17, and hence Definition 2.10, apply.
Let be the dimension of the domain of the pdes to be considered. Let be a bounded, open subset of , or of a -dimensional manifold, with boundary. Let , and let . For , let be the Sobolev space .
Let be a linear partial differential operator. (Here the subscript denotes compactly supported functions.) Let , with , be the ‘polynomial’ order of the nonlinearities in the pdes. Let
[TABLE]
be linear partial differential operators. For index-vector and , we set
[TABLE]
Let , , and be as in Section 2.5. Fix smooth functions222The real line may be replaced by a smaller open interval. , for , with , such that these functions and all their derivatives grow at most polynomially. Define by , where for each ,
[TABLE]
for and . Suppose that .
We write
[TABLE]
Then . The maps and extend continuously to . Suppose that that the eigenfunctions of this extension of form a Hilbert basis of .
Theorem 2.21**.**
The spaces and the maps and satisfy the hypotheses of Theorem 2.17.
We prove Theorem 2.21 in Section 4.4. Together with Theorem 2.17, it has the following immediate consequence.
Corollary 2.22**.**
Let be such that . Suppose that , , and satisfy and (as in Theorem 2.17). Then there are infinitely differentiable maps
[TABLE]
where is a polynomial map that separates invariant subspaces, such that if and are as in 2.7 and 2.6, then
[TABLE]
for all . Further, for every , there is a such that for all ,
[TABLE]
as in . There is a construction of the map in which all integrals over that occur are -regular.
This corollary shows that any pde of the form 2.2, with and as in this subsection, is equivalent up to a residual of order to a pde with clear invariant manifolds, as in Definition 2.10.
Example 2.23**.**
Suppose that , the circle. This amounts to imposing periodic boundary conditions. Take , and let by any linear partial differential operator with constant coefficients. Its eigenfunctions, for and , are orthogonal in the Sobolev spaces . For a map as in Theorem 2.21, that is, a polynomial expression in derivatives of functions, whose polynomial coefficients increase at most polynomially, Theorem 2.21 implies that the conditions of Theorem 2.17 are satisfied in this case, so Corollary 2.22 applies. This generalises directly to cases where is a higher-dimensional torus; that is, to problems in with periodic boundary conditions. Here we used the case where the domain is a manifold, rather than an open subset of .
Most of the rest of this paper is devoted to proofs of Theorems 2.17 and 2.21, and developing the tools used in these proofs. In Section 8, we prove Lemmas 2.8 and 2.9. In Section 9 we illustrate Corollary 2.22 by working out an example.
3 Derivatives and polynomials
In this section we review standard material on derivatives of maps between normed vector spaces. We also briefly discuss polynomial maps between normed vector spaces. Throughout this section, and are normed vector spaces, possibly infinite-dimensional. Let be an open subset, and let be a map. We fix an element .
3.1 First order derivatives
This subsection and the next contain some standard definitions and facts about derivatives of maps between normed vector spaces. Details and proofs can be found in various textbooks [35, e.g.].
For a map , we use the notation for the statement
[TABLE]
were runs over .
Definition 3.1**.**
The map is differentiable at , if there is an operator such that
[TABLE]
Then is the derivative of at . If is differentiable at every point in , then we say that is differentiable. In that case, the derivative of is the map
[TABLE]
mapping to .
The derivative of a map at a point is unique, if it exists.
Lemma 3.2** (Chain rule).**
Let be a third normed vector space. Let be an open subset containing . If is differentiable at and is differentiable at , then is differentiable at , and
[TABLE]
{ArXiv}
Proof.
For all such that , differentiability of at and of at imply that
[TABLE]
Since and are bounded operators, the second term on the right-hand side equals , while the last term is . ∎
Definition 3.3**.**
The map is a near-identity at if the map 3.1 is continuous in a neighbourhood of , and
[TABLE]
3.2 Higher order derivatives
Fix a positive integer . We write for the space of multilinear maps for which the norm
[TABLE]
is finite. There is a natural isometric isomorphism
[TABLE]
mapping an operator in the left-hand side to the operator given by
[TABLE]
for .
Suppose is differentiable. The map is twice differentiable at if the map 3.1 is differentiable at . Then we write
[TABLE]
Inductively, for , is defined to be times differentiable at if it is times differentiable, and the map
[TABLE]
is differentiable at . We then set
[TABLE]
In this case, we write
[TABLE]
As before, we say that is times differentiable if it is times differentiable at every point in . And infinitely differentiable means times differentiable for every .
Theorem 3.4** (Taylor’s theorem).**
Suppose is times differentiable. Suppose that for all in a closed ball around contained in . Then for every in this ball,
[TABLE]
{ArXiv}
3.3 Example: Burgers’ equation
An example of a map to which we would like to apply the material in this section and the next is the nonlinear term in Burgers’ equation
[TABLE]
Let be a bounded, open interval in . For every , consider the th -Sobolev space , with the inner product
[TABLE]
Consider the map given by
[TABLE]
First of all, for , the Cauchy–Schwartz inequality for (or Hölder’s inequality) implies that
[TABLE]
So indeed maps into .
We claim that is infinitely differentiable. Indeed, for ,
[TABLE]
And by 3.6, so
[TABLE]
The map is bounded, because, analogously to 3.6,
[TABLE]
If , then
[TABLE]
So is the operator in given by
[TABLE]
The term in the definition of the derivative is zero in this case, and that does not depend on . This implies that for every , . So is indeed infinitely differentiable.
3.4 Bounded polynomial maps
An operator in is said to be symmetric if it is invariant under permutations of its arguments. Let be the subspace of symmetric operators in . An example of such a symmetric operator is the th derivative of a map.
Lemma 3.5**.**
If is times differentiable at , then is symmetric.
We denote the permutation group of by . {ArXiv}
Lemma 3.6**.**
The subspace is closed.
Proof.
If , and and are such that set
[TABLE]
Let be such that , for the norm 3.2. Then symmetry of and the triangle inequality imply that
[TABLE]
a contradiction. ∎
By this lemma, is a Banach space if and are.
Let be the symmetrisation operator: for every and ,
[TABLE]
{ArXiv}
(An alternative proof of Lemma 3.6 is to show that is continuous, and to note that is the zero level set of minus the identity.)
{journal} The operator is continuous, and that is the zero level set of minus the identity, and hence closed in . So is a Banach space if and are.
An element defines a map by
[TABLE]
We have , and the map is injective on .
Definition 3.7**.**
A bounded homogeneous polynomial map of degree from to is a map of the form as in 3.8. We write for the space of such maps. It inherits a norm from the space via the linear isomorphism . If , then the degree of is .
A bounded polynomial map from to is a finite sum of bounded homogeneous polynomial maps. The degree of a bounded polynomial map is the degree of its highest-degree homogeneous term.
We write for the space of all bounded polynomial maps from to . This is the algebraic direct sum of the spaces .
{ArXiv}
By Lemma 3.6, is a Banach space if and are.
{journal} Because is a Banach space if and are, so is .
We could define as the space of constant maps into , but we only consider homogeneous polynomials of order at least one.
If is times differentiable at , then we have the map
[TABLE]
Lemma 3.8–3.11 below are basic facts showing that bounded polynomials and their orders and compositions behave as one would expect. {journal} Their proofs are short and straightforward.
Lemma 3.8**.**
Every bounded polynomial map is infinitely differentiable.
{ArXiv}
Proof.
Let , for some . Then for all ,
[TABLE]
Hence is differentiable, and
[TABLE]
where on the right-hand side, the operator is applied to copies of , to give an element of . Hence is a bounded polynomial map in . This proves the claim by induction. ∎
Lemma 3.9**.**
If , then there is a constant such that for all ,
[TABLE]
{ArXiv}
Proof.
Let . By boundedness and multilinearity of , we have for all nonzero ,
[TABLE]
∎
Lemma 3.10**.**
If is a polynomial map from to of order lower than , and
[TABLE]
as in , then .
{ArXiv}
Proof.
Let . Let , and suppose that , as in . Then there is a such that for all with unit norm and small enough,
[TABLE]
If , this implies that . ∎
Lemma 3.11**.**
If and , then .
{ArXiv}
Proof.
For and , define by
[TABLE]
for . Then one checks directly that . This implies the claim about polynomials. ∎
3.5 Standard monomials
Let be the continuous dual of . We denote the pairing between and by . For every , let be given. What follows is most natural if is a Hilbert space and is given by taking inner products with an element of a Hilbert basis, but it applies more generally.
Consider a multi-index . If , and is the largest number for which , then we define the element
[TABLE]
In other words, for all ,
[TABLE]
We write for the corresponding homogeneous polynomial. One could call this the standard -monomial with respect to the set . (If and the elements are the standard coordinates, then the monomial functions in the usual sense are precisely the scalar multiples of the maps .)
For , we write
[TABLE]
This product is finite (since has finitely many nonzero terms) and depends on the set . The following lemma follows from the definition of the derivative.
Lemma 3.12**.**
The derivative of in 3.10 is given by
[TABLE]
for all .
4 Compact derivatives and polynomials
It is a nontrivial question in what sense differentiable maps between normed vector spaces can be approximated by polynomial maps [1, 11, 12, 13, 22, 23, e.g.]. In this section we discuss an approach to this problem that is suitable for our purposes. This discussion includes the further problem of approximating a polynomial by sums of the standard monomials of Section 3.5. The polynomials for which this is possible are the compact polynomials introduced in Section 4.2.
Section 4.3 introduces compactly differentiable maps. We combine these with Taylor’s theorem to express the lowest order parts of such maps in terms of standard monomials. We discuss a class of examples of compactly differentiable maps relevant to the study of pdes.
4.1 Compact multilinear maps
Let and be Banach spaces. Let be the image of the space
[TABLE]
under the isomorphism 3.3. {journal} Using induction on , one can show that is closed in , and hence a Banach space.
{ArXiv}
Lemma 4.1**.**
For every the space is closed in .
Proof.
We use induction on . For the claim is standard. Suppose the claim holds for . Then
[TABLE]
which is a closed subspace of . And that space is closed in since is closed in by the induction hypothesis. ∎
By Lemma 4.1, is a Banach space.
Let and be countable subsets whose spans are dense. (So and are separable.) For any , consider the multilinear map
[TABLE]
A Banach space has the approximation property if every compact operator on the space is a norm-limit of finite-rank operators. This is always true for Hilbert spaces, but we need to consider more general Banach spaces for applications. {journal} The following result is standard in the case where and are Hilbert spaces.
Proposition 4.2**.**
If has the approximation property, then the space is dense in .
Proof.
Since has the approximation property, the space of finite-rank operators (linear operators with finite-dimensional images) is dense in . See for example Proposition 4.12(b) in the book by Ryan [31]. The space is dense in the space of finite-rank operators, so the claim follows. ∎
Lemma 4.3**.**
If has the approximation property, then for every , the span of is dense in .
Proof.
We prove this by induction on . If , then the claim is precisely Proposition A.1{ArXiv} in the appendix. Now suppose that the claim holds for a given . By definition,
[TABLE]
By the induction hypothesis, the set has dense span in . Therefore, Proposition A.1, with replaced by , implies that the set
[TABLE]
has dense span in . This is precisely the claim for . ∎
A Schauder basis of a Banach space is a subset such that for each , there are unique complex numbers such that
[TABLE]
A space with a Schauder basis has the approximation property.
Lemma 4.4**.**
If is a Schauder basis of and is a Schauder basis of , then for every , the set is a Schauder basis of .
Proof.
Lemma 4.3 implies that has dense span. So it remains to show that if are such that
[TABLE]
then for all and . Since is a Schauder basis of , this reduces to the case where . {journal} In that case, one can prove the claim by induction on , using the fact that is a Schauder basis of .
{ArXiv} We prove the claim in that case, by induction on .
If , then the claim follows since is a Schauder basis of . So suppose that the claim holds for a given , and let be such that
[TABLE]
Then for all ,
[TABLE]
Since, is a Schauder basis of , this implies that for every ,
[TABLE]
Because the sum
[TABLE]
converges in , we find that the sum converges to zero in this space. By the induction hypothesis, this implies that for every and .
∎
Remark 4.5**.**
In the induction step in the proof of the special case of Lemma 4.3 where is a Hillbert space, we still need the general version of Proposition A.1, where is a Banach space. This is because is only a Banach space, even if is a Hilbert space.
The subspace of symmetric operators in is closed in , since it is the intersection of the closed subspaces and {ArXiv} (Lemmas 3.6 and 4.1).
Hence is a Banach space with respect to the norm 3.2.
A Schauder basis of a Banach space is unconditional if there is a constant such that for all with , and all ,
[TABLE]
In that case, convergence of implies convergence of , for every .
Lemma 4.6**.**
Suppose that and are Hilbert spaces, and that and are orthogonal sets in and respectively, with dense spans. Let be defined by taking inner products with . Then is an unconditional Schauder basis of .
Proof.
The set is a Schauder basis of by Lemma 4.4. It remains to show that it is unconditional. By rescaling the vectors and , we reduce the proof to the case where and are Hilbert bases. In that case, for all finite subsets and all ,
[TABLE]
∎
Lemma 4.7**.**
Let , and be normed vector spaces, and . Let , and . Define by
[TABLE]
for all . Then .
Proof.
We use induction on . For , because the composition of a compact operator and a bounded operator is compact. Suppose that the claim holds for a given . Let , and . For a fixed , define by
[TABLE]
for . For a fixed , define by
[TABLE]
for . Then for all ,
[TABLE]
So by the induction hypothesis, . In this way, we obtain the map
[TABLE]
mapping to . It remains to show that is a compact operator.
Define by
[TABLE]
for and . (This map takes values in by the induction hypothesis.) Since is compact and is bounded, we find that is a compact operator. ∎
4.2 Compact polynomial maps
Definition 4.8**.**
A compact homogeneous polynomial map of degree from to is a map of the form as in 3.8, for . We write for the space of such maps. This space inherits a norm from the space it is contained in.
If , then the degree of is . A compact polynomial map from to is a finite sum of compact homogeneous polynomial maps. The degree of a compact polynomial map is the degree of its highest-degree homogeneous term. We write for the space of all compact polynomial maps between these spaces.
The isometric isomorphism restricts to an isometric isomorphism . So is a closed subspace of the Banach space , and hence is a Banach space itself.
For every , an operator of the form , with as in 3.9, is an element of . Indeed, is an iteration of rank-one operators, So .
The following proposition is the reason why we are interested in compact polynomial maps.
Proposition 4.9**.**
*Suppose that and are Banach spaces, that has a Schauder basis , and that has a Schauder basis . Then the elements *
[TABLE]
where the multi-index ranges over the elements of with , and ranges over the positive integers, form a Schauder basis of .
Proof.
Consider the space
[TABLE]
Lemma 4.4 implies that the set of where is a Schauder basis of . And is a bounded linear isomorphism with bounded inverse. Since such isomorphisms map Schauder bases to Schauder bases, we find that the set , for non-decreasing as above, is a Schauder basis of .
For with non-decreasing entries, define by
[TABLE]
Then . (Note that , for , and , for , are defined differently; compare 3.9 and 4.1.) Every sequence in occurs in exactly one way as , for alpha as above, so , where and , is a Schauder basis of . Since , the claim follows. ∎
A reformulation of Proposition 4.9 is that for every compact polynomial map , there are unique complex numbers such that
[TABLE]
where the sum converges in the norm on . Conversely, all polynomial maps of this form are compact.
Lemma 4.10**.**
In the setting of Lemma 4.6, the set is an unconditional Schauder basis of .
Proof.
The proof is analogous to the proof of Proposition 4.9, where we now use Lemma 4.6 instead of Lemma 4.4, and we use the fact that bounded linear isomorphisms with bounded inverses map unconditional Schauder bases to unconditional Schauder bases. ∎
Lemma 4.11**.**
If and , then .
Proof.
The proof is similar to the proof of Lemma 3.11, with bounded multilinear maps replaced by compact ones. ∎
Remark 4.12**.**
Other notions of compact polynomial maps were studied by Gonzalo, Jaramillo and Pełczyńsky in [15, 25].
4.3 Compactly differentiable maps
Let and be normed vector spaces, let be an open subset containing a vector , and let be times differentiable at .
Definition 4.13**.**
The map is times compactly differentiable at if
[TABLE]
If is times compactly differentiable at , then by Lemma 3.5,
[TABLE]
Then the map of 3.4 is the compact polynomial map associated to . Together with Theorem 3.4 and Proposition 4.9, this leads to the following conclusion.
Corollary 4.14**.**
Suppose that and are Banach spaces, that has a Schauder basis , and that has a Schauder basis . Suppose is times differentiable, and times compactly differentiable for every . Then there are unique complex numbers such that
[TABLE]
where the part of the sum where converges as a function of in the norm on , for .
(Note that, in 4.3, on the left-hand side is a map from to , whereas on the right-hand side, is an element of .)
Lemma 4.15**.**
A compact polynomial map is infinitely compactly differentiable.
Proof.
We show that the derivative of every homogeneous compact polynomial , for , is a compact polynomial in . This implies the claim by induction on . As in the proof of Lemma 3.8, for all . In other words, , with , where we view as an element of . This shows that . ∎
Lemma 4.16**.**
Let , and be normed vector spaces and let and be differentiable maps. If either or is compactly differentiable, then so is .
Proof.
Lemma 3.2 implies that for all ,
[TABLE]
If is compactly differentiable, then . If is compactly differentiable, then . In either case, we find that . ∎
Lemma 4.17**.**
Let , and be normed vector spaces, let be times differentiable, and suppose that is a subspace of , with compact inclusion map . Then is times compactly differentiable as a map from to .
Proof.
Let . Then
[TABLE]
So the claim follows from Lemma 4.7. ∎
Remark 4.18**.**
It is possible for a map to be times compactly differentiable, but not times. For example, the first derivative of the identity operator on an infinite-dimensional Banach space is the identity map itself, and not a compact operator. But its higher-order derivatives are zero, and hence compact.
4.4 A class of compactly differentiable maps
We end this section by discussing a class of compactly differentiable maps (specifically, compact polynomials) that are relevant to the study of nonlinear pdes. These maps are polynomial expressions in derivatives of functions; see Proposition 4.20 below.
Let be a bounded open subset with boundary. For and , consider the Sobolev space .
Lemma 4.19**.**
Let and . Pointwise multiplication of functions defines a map .
Proof.
Let . By Hölder’s inequality, for all ,
[TABLE]
For , there are combinatorial constants , for , with such that , such that for all , we have the generalised Leibniz rule
[TABLE]
Together with 4.4, this implies that
[TABLE]
∎
Proposition 4.20**.**
Let and , with . Let
[TABLE]
be linear partial differential operators of orders smaller than . Fix complex numbers , for with , and . Define by , where for each ,
[TABLE]
for , and with as in 2.9, defines a compact polynomial map
[TABLE]
Proof.
We first consider the case where , and if , and zero otherwise. By Lemma 4.19, pointwise multiplication defines a map .
By Rellich’s lemma, boundedness of implies that the maps extend to compact operators
[TABLE]
The map
[TABLE]
defined by
[TABLE]
for , is in by Lemma 4.7. Hence is an element of .
Every component of a general map of the form 4.5 is a finite sum of maps of the form as above, applied to the components of . Hence it is in . ∎
Example 4.21**.**
{ArXiv}
Consider the map from Section 3.3, mapping to . We now view as a map from to , for , and .
{journal} Let , and be integers. Let be a bounded open interval. Consider the map from to , mapping to .
Taking , , and , for , in Proposition 4.20, we find that is a compact polynomial in for every . Hence, by Lemma 4.15, is in particular infinitely compactly differentiable. For , the map is only a bounded polynomial in .
Remark 4.22**.**
Proposition 4.20 extends directly to relatively compact open subsets of manifolds. The latter extension is relevant, for example, if one uses periodic boundary conditions, so that one works with with functions on a torus.
5 Polynomial and differentiable maps on graded Fréchet spaces
Apart from polynomial and differentiable maps between normed vector spaces, we also use such maps between graded Fréchet spaces, defined in terms of nested sequences of Banach spaces, as in Definition 2.1. In this section, we discuss some further properties of such spaces, and in particular what it means for two sequences of Banach spaces defining such a space to be comparable.
5.1 Properties of nested sequences of Banach spaces
Let be a nested sequence of Banach spaces, as in Definition 2.1. Their intersection is a graded Fréchet space.
Definition 5.1**.**
The space consists of the multilinear maps (with factors ) such that for every , there is a such that extends continuously to a map in .
Note that .
Definition 5.2**.**
We write for the space of all maps such that for every , there is a such that extends continuously to a polynomial in . The space is defined analogously.
A feature of the spaces and that is useful to us, is that they admit natural compositions. If and , then Lemma 3.11 implies that lies in . Similarly, Lemma 4.11 implies that if and .
Definition 5.3**.**
An times (compactly) differentiable map from to itself is a map such that for every , there is a such that extends to an times (compactly) differentiable map from to .
Lemmas 3.2 and 4.16 imply that the spaces of differentiable and compactly differentiable maps from to itself are closed under composition. Lemma 4.17 implies that if is a compactly nested sequence of Banach spaces, then all times differentiable maps from to itself are times compactly differentiable.
It follows directly from Definition 2.4 that if are of orders and , respectively, then is of order . We also have the following lemma.
Lemma 5.4**.**
Consider two maps , where is differentiable and is of order . Suppose that for every , there is a such that is differentiable, and as in . Then the map , mapping to is of order .
Proof.
Let . Let be such that is differentiable, and its derivative satisfies the estimate in the lemma. Let be such that as in . Then as in . ∎
An example of a situation where the condition on in Lemma 5.4 is satisfied is the following.
Lemma 5.5**.**
Let , for . Then for every , there is a such that as in .
Proof.
Let , and let be such that . As in the proof of Lemma 3.8, . So the claim follows from Lemma 3.9. ∎
Remark 5.6**.**
Everything in this subsection generalises directly to polynomial and (compactly) differentiable maps between two different Fréchet spaces that are given as intersections of nested sequences of Banach spaces. We will not need this generalisation, however.
5.2 Comparable sequences of Banach spaces
In the rest of this section, we discuss some relevant properties and examples of comparable sequences of Banach spaces (Definition 2.6). Particularly relevant to Theorem 2.17 are sequences of Banach spaces comparable to sequences of separable Hilbert spaces, which we discuss in Section 5.3. We will see relevant examples in Section 5.4.
Suppose that and are comparable nested sequences of Banach spaces. Then as sets.
Lemma 5.7**.**
The two spaces and are equal as Frećhet spaces.
Proof.
Let be a sequence in such that for every , . Let , and choose such that we have a bounded inclusion . Then there is a constant such that for every , , which goes to zero as . ∎
The following fact follows directly from the definitions, and the fact that the classes of maps in question are closed under composition with bounded linear maps.
Lemma 5.8**.**
If is a (compact) polynomial map or a (compactly) differentiable map, then it also defines a map of the same type on .
This lemma in particular states that as vector spaces. We will use the fact that this equality includes natural topologies on these spaces (Corollary 5.10 below) to prove Corollary 5.13.
Lemma 5.9**.**
For all , there is an such that for every with , there is a such that we have a bounded inclusion map
Proof.
Let . Choose and such that for every , . Let . Choose and such that for every , . Then for all ,
[TABLE]
∎
For a sequence in , we define in to mean that for every , there is a such that in . (This includes the requirement that for every .)
Corollary 5.10**.**
We have , including topologies.
Proof.
Let be a sequence in converging to zero. Let . Choose as in Lemma 5.9. Choose such that in . Choose as in Lemma 5.9.
For each , write , for . By Lemma 5.9, there is a such that for every ,
[TABLE]
which goes to zero as . Hence in . ∎
5.3 Sequences of Banach spaces comparable to sequences of Hilbert spaces
As before, we suppose that and are comparable nested sequences of Banach spaces. Now we make the additional assumption that the spaces are separable Hilbert spaces. Suppose that is a subset of that is orthogonal in all spaces , with dense span. Then taking inner products with defines bounded functionals, all denoted by , on all spaces , and hence in for large enough.
Corollary 5.11**.**
Suppose is an times differentiable map from to itself, and that is times compactly differentiable for every . Then there are unique complex numbers such that
[TABLE]
where is of order , and the part of the sum where converges as a function of in , for .
Proof.
Let . Choose be such that we have bounded inclusions and , and is times differentiable, and times compactly differentiable for every . Then the same is true for . So Corollary 4.14 implies that 5.1 holds, for unique , where the sum converges in , and hence in by Corollary 5.10, and is of order as a map from to itself, and hence as a map from to itself. ∎
Remark 5.12**.**
A useful feature of Corollary 5.11 is that it is not assumed that the spaces have the approximation property. The point is that the separable Hilbert spaces do have this property.
The following corollary is an important way in which we use comparable sequences of Banach spaces. It is used in the proof of Lemma 7.7.
Corollary 5.13**.**
Let be given such that
[TABLE]
converges in . Then for all subsets , the series
[TABLE]
converges in .
Proof.
By Corollary 5.10, the series 5.2 converges in , and it is enough to show that 5.3 converges in . And that follows from Lemma 4.10. ∎
Remark 5.14**.**
In Corollary 5.13, the two series converge to elements of .
5.4 Example: Sobolev spaces and -spaces
Let be a bounded open subset of or of an -dimensional Riemannian manifold, and suppose that the boundary of is . Set
[TABLE]
Lemma 5.15**.**
The above sequences and of Banach spaces are comparable.
Proof.
Set . By a Sobolev embedding theorem, we have bounded inclusions
[TABLE]
for all and every . Now for all such and ,
[TABLE]
So we have bounded inclusions
[TABLE]
Furthermore, since has finite volume, we have bounded inclusions for every and all .
The above arguments imply that we have bounded inclusions
[TABLE]
∎
Remark 5.16**.**
In this example, the spaces are separable Hilbert spaces.
For another example of comparable sequences of Banach spaces, fix . For , set
[TABLE]
The space is complete in the norm given by the maximum of the sup-norms of the partial derivatives of functions up to order .
Lemma 5.17**.**
The sequences and of Banach spaces are comparable.
Proof.
We have a bounded inclusion for every . So it remains to show that for every , there are such that we have bounded inclusions
[TABLE]
By a Sobolev embedding theorem, we have a bounded inclusion
[TABLE]
for all such that
[TABLE]
For , set and . Then this Sobolev embedding theorem yields the desired inclusions 5.4. ∎
Remark 5.18**.**
If , then we may take , so that the spaces are Hilbert spaces.
Proposition 4.20 and Lemma 5.15 together imply Theorem 2.21. The extension of Proposition 4.20 to coeficients depending on a real (time) parameter in a smooth way, and the extension of Lemma 5.15 to vector-valued functions, are straightforward.
6 A coordinate transform
6.1 A residual
Recall the setting of Section 2.2. In this section and the next, based upon the details of some given dynamical system 2.2 we construct both a coordinate transformation 2.6 and a corresponding ‘normal form’ system 2.7, such that solutions to 2.7, transformed by 2.6, satisfy the given dynamical system 2.2 up to residuals of a specified order . See Theorem 2.17. We do this inductively, by showing how to construct such a transformed system to satisfy 2.2 with residual of order from a version with residual of order .
In Section 7, we construct a more specific choice of the general coordinate transform constructed in this section, in order to establish exact invariant manifolds, and study their properties, for constructed systems arbitrarily close to the given system 2.2.
Remark 6.1**.**
In Section 2.2, we assumed that the sequence is comparable to a nested sequence of separable Hilbert spaces in which the vectors are orthogonal. Lemma 5.8 implies that we may equivalently assume that itself is a nested sequence of separable Hilbert spaces, because all maps from to itself we use transfer to maps from to itself of the same type (e.g. compact polynomial and compactly differentiable maps). However, the formulation where is only comparable to a nested sequence of separable Hilbert spaces makes it clearer that we have the flexibility to consider maps between Banach spaces. This is natural for example in the context of Proposition 4.20.
Let , with . Let be such that and are compact polynomial maps of order at most in the component, and infinitely differentiable in . Suppose, furthermore, that is a near-identity at zero, and that .
Recall that our goal is to relate the dynamics of maps satisfying 2.2 to the dynamics of maps satisfying 2.7 when and are related by the coordinate transform as in 2.6.
For maps , with differentiable, we write for the map from to given by
[TABLE]
for all and . (Note that this is different from .) If is a map from to to itself, then the composition is defined analogously. Also recall the notation for compositions of maps to and from and under Notation and conventions in Section 1.5.
Define the maps by
[TABLE]
The map is the residual of the transformed ode, in the following sense.
Lemma 6.2**.**
For all smooth maps satisfying 2.7, and with determined from by 2.6,
[TABLE]
Proof.
For and as in the lemma, the chain rule (Lemma 3.2) and 2.7 imply that for all ,
[TABLE]
∎
Lemma 6.3**.**
The maps and are infinitely compactly differentiable.
Proof.
Because the Banach spaces are compactly nested, it is enough to show that and are infinitely differentiable. And that is true by the chain rule, because is infinitely differentiable, and so are and , by Lemma 3.8. ∎
To recursively construct 2.7 and 2.6, suppose that . We proceed to show that 2.7 and 2.6 may be refined to make the new residual of . By Lemma 6.3, Corollary 5.11 (where and ), and Lemmas 3.9 and 3.10, there are unique, infinitely differentiable maps
[TABLE]
for all multi-indices with , and a map
[TABLE]
such that for all , and ,
[TABLE]
where the sum converges in , differentiably in , and . The monomial is defined as in 3.10, with respect to the set of functionals , for .
In Section 6.2, we construct maps such that the order term in 6.1 may be replaced by an order term , if 2.7 and 2.6 hold with replaced by .
6.2 Construction of the coordinate transform
For each , let be the eigenvalue of corresponding to . For all with , define by
[TABLE]
(This sum has at most nonzero terms.) For such , let be as in 6.2. Let be smooth maps such that
[TABLE]
Suppose that the sums
[TABLE]
converge in , differentiably in .
Define a new coordinate transform map and corresponding map that replaces in 2.7, by
[TABLE]
The following result is the main step in the construction of the coordinate transform we are looking for.
Proposition 6.4**.**
If and are as in 2.7 and 2.6, with replaced by , then 6.1 holds, with the residual replaced by a residual satisfying
[TABLE]
The maps and are compact polynomials in of order at most , and is a near-identity.
Remark 6.5**.**
The maps and can be found explicitly if we decompose 6.4 with respect to the basis . This will be done in Section 7. One solution to 6.4 is and . However, for our purposes, we need the function to be of a specific form. The main purpose of this work is to find such that the -component of is zero for certain combinations of and , in such a way that an exact separation of stable, centre and unstable modes is maintained. See Proposition 7.1.
6.3 Proof of Proposition 6.4
Define the map by
[TABLE]
Lemma 6.6**.**
We have
[TABLE]
Proof.
The left-hand side of 6.7 equals
[TABLE]
By Lemma 5.5, the derivative satisfies the condition of Lemma 5.4, with . Since , Lemma 5.4 implies that . Similarly, . Now is a polynomial map, and a near-identity. By Lemma 5.5, this implies that satisfies the condition of Lemma 5.4, with . So Lemma 5.4 implies that . ∎
Lemma 6.7**.**
For such that , let be as in 6.2. Then for all and ,
[TABLE]
Proof.
First of all, By Lemma 3.12,
[TABLE]
Now, because the vectors are orthogonal with respect to , we have for every . So
[TABLE]
We find that the left-hand side of 6.8 equals
[TABLE]
with as in 6.3. So the claim follows from 6.4. ∎
Define the maps by
[TABLE]
Lemma 6.8**.**
The residual satisfies
Proof.
By Lemmas 6.6 and 6.7,
[TABLE]
where, for and ,
[TABLE]
By 6.2, the last expression in 6.9 equals
[TABLE]
Let be given, and choose such that is differentiable, and . Using Theorem 3.4, we write
[TABLE]
for and . Lemma 3.9 implies that uniformly in in compact sets. So
[TABLE]
uniformly in in compact sets. And then, as in the proof of Lemma 5.4, the assumption 2.1 and Lemma 5.4 imply that , uniformly in in compact sets. ∎
Proof of Proposition 6.4.
The correction terms and lie in . Hence and are compact polynomials, because and are. By Lemma 3.9, this also implies that is a near-identity because is. The desired property of is Lemma 6.8. ∎
7 Centre, stable and unstable coordinates
There is considerable flexibility in choosing the maps and in Section 6.2. In this section, we discuss how to make specific choices, in terms of the eigenvalues of , so that the normal form 2.7 is useful for detecting invariant manifolds.
7.1 Centre, stable and unstable components of
We use the notation from Section 2.5. In particular, let , , and be spectral gap parameters defined there. Recall the definition of polynomial growth in Definition 2.16.
Proposition 7.1**.**
Suppose that and . Suppose that has polynomial growth. The maps and in Section 6.2 can be chosen such that
- •
if either and , or and , then ;
- •
if , then ; and
- •
if , then .
Proposition 7.1 is proved in Section 7.2, after some preparation done in this subsection.
Definition 7.2**.**
Let , such that . Set and . Let be a continuous function on such that as if and as if . Then we define the function on by
[TABLE]
The integrals occurring in this definition are -regular, in the sense defined in Section 2.5.
Lemma 7.3**.**
In the setting of Definition 7.2,
[TABLE]
Proof.
This is a straightforward computation. ∎
Lemma 7.4**.**
Let be a smooth function, and suppose that and all its derivatives grow at most polynomially. Then, for every as in Definition 7.2, and all its derivatives grow at most polynomially.
Proof.
First of all, because , it is enough to consider the function itself rather than all its derivatives.
Let and be such that for all , . We prove by induction on that there is a constant such that for all , . We consider the case where ; the case where is similar.
If , then for all ,
[TABLE]
Suppose that the claim holds for , and suppose that for a constant . Using integration by parts, we find that
[TABLE]
which implies the claim by the induction hypothesis. ∎
Let , with . Consider the differentiable maps
[TABLE]
where the sum converges in , uniformly and differentiably in in compact sets in .
For with , let be the set of such that either
- •
and either and , or and ;
- •
and ; or
- •
and .
For with , and , write , with as in 6.3.
Lemma 7.5**.**
If and , then for every , .
Proof.
If and , and , then
[TABLE]
If and , and , then
[TABLE]
If , and , then
[TABLE]
And if , and , then
[TABLE]
∎
7.2 Update terms for and
Suppose that has polynomial growth. Then the functions and their derivatives grow at most polynomially, uniformly in and .
For every with and , consider the ode for and
[TABLE]
Define the maps as follows. If , then
[TABLE]
This definition makes sense because of Lemma 7.5 and the growth behaviour of the functions . If , then we set
[TABLE]
Lemma 7.6**.**
With the above definitions, the ode 7.1 is satisfied for all and .
Proof.
If , the statement is immediate from the definitions. If , it follows from Lemma 7.3. ∎
Lemma 7.7**.**
*Suppose that and . Then the sums *
[TABLE]
converge in , differentiably in .
Proof.
The first sum in 7.4 equals
[TABLE]
Since is comparable to a nested sequence of separable Hilbert spaces in which the set is orthogonal, Corollary 5.13 implies that this series converges in .
Write Lemma 7.5 states that . So the second sum in 7.4 equals
[TABLE]
Tonelli’s theorem implies that convergence of the first of these sums is equivalent to convergence of
[TABLE]
The sum inside the brackets converges in , uniformly in , by convergence of 6.2 and Corollary 5.13. Since the functions grow at most polynomially, uniformly in and , the value of that sum grows at most polynomially as well, when viewed as a convergent series in . So the integral over converges in , by completeness of the spaces . By continuity of 7.7 in , the convergence is uniform in on compact subsets of . The derivatives of 7.7 with respect to are linear combinations of 7.7 and
[TABLE]
and therefore converge as well.
By an analogous argument, the second sum in 7.6 converges as well, differentiably in . ∎
Proposition 7.1 follows from Lemmas 7.6 and 7.7.
7.3 Proof of Theorem 2.17
Lemma 7.8**.**
If and have polynomial growth, then so do , , and .
Proof.
If has polynomial growth, then the functions and all their derivatives grow at most polynomially, uniformly in and . Hence, by 7.2 and 7.3, the map has polynomial growth. By Lemma 7.4, the same is true for . So and have polynomial growth.
Polynomial growth is preserved under composition and derivatives in the and directions. Hence the map as in 6.6 has polynomial growth, and therefore so do and . ∎
Combining Lemma 7.8 with Propositions 6.4 and 7.1, we prove the following slightly more explicit version of Theorem 2.17.
Theorem 7.9**.**
*Let be such that , and . Suppose that has polynomial growth. Then there are infinitely differentiable maps *
[TABLE]
where , is a polynomial map that separates invariant subspaces, is a near-identity and and are compact polynomials of orders at most , such that if and are as in 2.7 and 2.6, then 6.1 holds. Finally, there is a construction of the map in which all integrals over that occur are -regular.
Proof.
We use induction on to prove that the claim holds for every , including the auxiliary statement that and have polynomial growth.
If , then we may take and for all and . Then , so and have polynomial growth because does.
The induction step follows from Lemmas 7.8, 6.4 and 7.1. ∎
8 Dynamics of the normal form equation
It remains to prove Lemma 2.8 and Proposition 2.9, which we use to justify Definition 2.10 based on Theorem 2.17. Throughout this section, we suppose that is a smooth map that separates invariant subspaces.
Proof of Lemma 2.8..
First, suppose that . For all and all with and or , we have . So the properties 2.3 of the map imply that . This, in turn, implies that for all maps satisfying 2.7, if for a given then . So for all .
Next, suppose that . If and , then if . So for all . And the components of with are zero for the same reason, while its components with are zero since . Hence . We conclude that . As in the case , this implies the claim for .
The argument for is entirely analogous to the case . ∎
To prove Proposition 2.9, we start with a general comparison estimate for solutions of odes in Hilbert spaces.
Lemma 8.1**.**
Let be a Hilbert space, a subspace, an open interval containing [math], and a map from into the space of linear operators from to . Let be a differentiable map (as a map into ), such that for all ,
[TABLE]
If is such that is negative semidefinite for all and , then for all with ,
[TABLE]
Proof.
First, suppose that . Then for all ,
[TABLE]
So is a nonnegative, non-increasing function on , and the claim for follows.
Next, let be arbitrary. Then
[TABLE]
Applying the claim for , with replaced by and by , now yields the claim for . ∎
For any homogeneous polynomial map between normed vector spaces and , where , define the map by
[TABLE]
Here copies of are inserted into on the right hand side.
For all , the map lies in for some . The operator lies in for some . By replacing the smaller of or by the larger of these two numbers, we henceforth assume . Applying the construction 8.1 to each homogeneous term of and adding the resulting maps, we obtain a map , such that for all ,
[TABLE]
For , we write for composed with orthogonal projection onto .
For , let be defined by if , and otherwise.
Lemma 8.2**.**
Let . Write , where for . Then for all , the components of in , and satisfy
[TABLE]
Proof.
Let and . To prove 8.2, we use the fact that by 2.3b,
[TABLE]
So
[TABLE]
where is the basis of dual to . Hence
[TABLE]
which implies 8.2. The equality 8.3 can be proved anaogously.
To prove 8.4, we note that by 2.3a,
[TABLE]
where
[TABLE]
The right hand side of 8.5 only depends on , and the right hand side of 8.6 is zero if or . So, under that condition, ∎
Proof of Proposition 2.9..
Set
[TABLE]
Because is a sum of polynomials of degrees at least two, we have for all . So contains . It is open by continuity of .
Let be a solution of the constructed system 2.7. As in the proof of Lemma 2.8,
[TABLE]
where we used the first equality in Lemma 8.2 and the fact that preserves . For all , the operator
[TABLE]
is negative semidefinite. Hence the claim about follows from the second part of Lemma 8.1. The claim about can be proved similarly, via a version of Lemma 8.1 for positive-definite operators.
Next, suppose that or . By Lemma 2.8, either for all or for all . Similarly to 8.8, the third equality in Lemma 8.2 implies that
[TABLE]
for all . And for all , the operator
[TABLE]
is negative semidefinite. So by Lemma 8.1,
[TABLE]
for all in . It similarly follows that for all in ,
[TABLE]
∎
9 Example: a non-autonomous version of Burgers’ equation
Let , and consider the non-autononomous, nonlinear pde
[TABLE]
with -periodic boundary conditions in . Then Theorem 2.21 applies, where is the circle.
Using Theorem 2.17, we compute the centre manifold of the normal form system approximating 9.1 up to residuals of order three, in Section 9.1. Via a direct approach, we compute all invariant manifolds for residuals of orders three and four, in Section 9.2. We find that the order three centre manifolds computed in the two ways agree. These computations illustrate Remark 2.20, that the construction from Theorem 2.17 is guaranteed to give a result, while a direct computation may be more efficient in concrete situations.
9.1 Centre manifold via Theorem 2.17
In this setting,
[TABLE]
where a prime denotes the derivative in the -direction. The eigenfunctions of are , for , given by . The eigenvalue corresponding to is (which has multiplicity two when ). Choose and such that , and . Suppose that lies within of an integer of the form , for a nonzero . Then the eigenvalue is central up to precision .
We determine a corresponding centre manifold for a system that approximates 9.1 up to a third-order residual. This involves the coordinate transform . To compute this centre manifold, we only need to apply to elements of . In other words, we only need to compute , for and . (We do not determine the domain here.)
For , the map is the identity map. So
[TABLE]
where
[TABLE]
For , let be defined by if , and otherwise. Then, for with ,
[TABLE]
So
[TABLE]
The map is expressed in terms of the map in
[TABLE]
(The order three term in 6.2 now equals zero.) See 7.2 and 7.3. If , then
[TABLE]
The equality holds precisely if and . Hence
[TABLE]
and if . An analogous argument shows that
[TABLE]
and if . The equality holds precisely if and either or . Hence
[TABLE]
and if .
The relevant numbers as in Section 7.1 equal
[TABLE]
(note that for every ). Because and , the real parts of and are greater than , whereas the real part of is smaller than .
And with as in Section 7.1, we have . Indeed,
[TABLE]
so and . Similarly, . And and , so .
Hence, by 7.2,
[TABLE]
The integral converges since , and is -regular. Similarly,
[TABLE]
And because ,
[TABLE]
We conclude that for all and ,
[TABLE]
(The last term is a scalar multiple of the constant function .) If , this simplifies to
[TABLE]
9.2 Invariant manifolds via direct computations
For order of residual , the map is the identity map, .
Proceeding to order of residual we construct quadratic corrections to the identity to form . In the eigenvector basis the field (all sums in this section are over ), and the pde 9.1 becomes
[TABLE]
Writing
[TABLE]
and solving{journal}333The computer algebra code used for the computations in this section is available on http://www.maths.adelaide.edu.au/anthony.roberts/pBurgers.txt.
for such that satisfies 9.3 up to terms of order three if , we find that is given by
[TABLE]
For and odd , the denominators in are not small. Then this map, combined with the linear , matches the pde 9.1 to third-order errors.
However, for and even , some denominators are small, becoming zero when . Then the divisor being zero becomes and hence has zeros for every pair of integer factors of (including negative pairs). Consequently these terms are excluded from the sum 9.4, and instead lead to nonlinearly modifying the evolution for some via
[TABLE]
Often the centre manifold is of most interest, so in setting all except , gives the quadratic approximate centre manifold to be for all except
[TABLE]
This is the same result as 9.2.
Proceeding to order of residual we may construct cubic corrections to to form . For simplicity, restrict attention to the cases of odd. It is straightforward but tedious to construct that for
[TABLE]
The terms excluded from in the sum 9.5 must cause cubic terms in the evolution. For example, when then444The apparent pattern in these odes becomes more complicated—at for example.
[TABLE]
By construction, in this case of , the coordinate transform 9.5 together with the odes 9.6 creates a dynamical system in which is the same as the pde 9.1 to a residual of order four. In the combined system 9.5 and 9.6, by definition Definition 2.10 three invariant manifolds are: the 1D unstable manifold parametrised by with all other ; the 2D centre manifold parametrised by with all other ; and the stable manifold with .
{ArXiv}
See Appendix B for the computer algebra code used the for the computations in this subsection. It is also available on http://www.maths.adelaide.edu.au/anthony.roberts/pBurgers.txt.
Acknowledgement
Part of this research was supported by the Australian Research Council grant DP150102385.
{ArXiv}
Appendix A Compact and finite-rank operators into Banach spaces
Let and be Banach spaces, and suppose that has the approximation property (this is true for example if is a Hilbert space). Let and be countable subsets with dense spans. (So and are separable.)
In the main text, we use the following, which is standard in the case where and are Hilbert spaces.
Proposition A.1**.**
The space is dense in .
Let be the space of finite-rank linear operators from to ; that is, operators whose images are finite-dimensional.
Lemma A.2**.**
The space is dense in .
Proof.
Let . Since the image of is finite-dimensional, there are and such that
Let . For every , choose and and such that
[TABLE]
Using the triangle and Cauchy–Schwartz inequalities, one finds that for all ,
[TABLE]
∎
Proof of Proposition A.1..
Since has the approximation property, is dense in . See for example Proposition 4.12(b) in the book by Ryan [31]. So the claim follows from Lemma A.2. ∎
{ArXiv}
Appendix B Computer algebra code for Burgers example computation
1Comment periodic Burgers-like, infinite-D coord transform.
2Written in Reduce, see https://reduce-algebra.sourceforge.io/
3AJR, 3 Jun 2019;
4on div; off allfac; on revpri;
5clear sum;
6
7operator xx; depend xx,t;
8let df(xx(~i),t)=>sub(j=i,dxxjdt);
9operator a; write "a(j) for eigenvalues alpha_j";
10operator b; write "b(j,k)=k*(j-k)";
11let b(j,j-~k)=>b(j,k);
12factor o;
13
14operator sum; linear sum;
15let df(sum(~a,~b),~t)=>sum(df(a,t),b);
16
17operator d;
18write "d(j,k)=1/(-a(j)+a(k)+a(j-k))=1/(n^2+2jk-2k^2)";
19write "d(j,k,l)=1/(-a(j)+a(k)+a(l-k)+a(j-l))
20 =1/2/(n^2+jl+kl-k^2-l^2)";
21let { d(j,k)a(k) => 1-d(j,k)(-a(j)+a(j-k))
22 , d(j,k,l)a(k) => 1-d(j,k)(-a(j)+a(l-k)+a(j-l))
23 };
24
25write "Linear approx evolution";
26xj:=xx(j);
27dxxjdt:=a(j)*xx(j);
28
29write "Quadratic solution shape";
30let o^2=>0;% order of error
31xj:=xj+o/2sum(b(j,k)(d(j,k)*t-d(j,k)^2)*xx(k)*xx(j-k),k);
32resj:=-df(xj,t)+a(j)*xj
33 +ot/2sum(b(j,l)*sub(j=l,xj)*sub(j=j-l,xj),l);
34resj:=sub(l=k,resj);
35
36write "However, zero divisors occur in d(j,k)
37=1/(-a(j)+a(k)+a(j-k)) =1/(n^2+2jk-2k^2). But not for odd
38n. For even n>0, find that zero divisors come from all the
39factors of n^2/2 since divisor=0 rewrites as k(k-j)=n^2/2.
40Both the positive and negative factors contribute. For
41example, n=2 gives n^2/2=2 with factors 2,1 and -2,-1 and
421,2 and -1,-2. So zero divisor for k=2,j=1 and k=-2,j=-1
43and k=1,j=-1 and k=-1,j=1.";
44
45write "Cubic manifold shape---pv. zero divisors are avoided
46in the quadratic terms by perhaps requiring n be odd.";
47let o^3=>0;% order of error
48resj:=-df(xj,t)+a(j)*xj
49 +ot/2sum(b(j,l)*sub(j=l,xj)*sub(j=j-l,xj),l)$
50depend k,kl; depend l,kl;
51resj:=(resj where sum(~~asum(~b,~l),~k)=>sum(ab,kl))$
52resj:=(resj where sum(~b,l)=>sum(sub(l=k,b),k))$
53resj:=(resj where sum(~a,kl)=>sum(sub(l=j-l,a),kl)
54 when df(a,xx(l)) neq 0)$
55write "Extract the coefficient of X(j-l)X(l-k)X(k)";
56rhsjkl:=(resj where sum(~a,kl) => coeffn( coeffn( coeffn(
57 a,xx(k),1) ,xx(j-l),1) ,xx(l-k),1));
58write "Iterative soln of ODE for coeff---zero divisors excluded";
59cjkl:=0;
60for it:=1:9 do begin
61 write resc:=d(j,k,l)*(rhsjkl-df(cjkl,t))-cjkl;
62 cjkl:=cjkl+resc;
63 if resc=0 then write "success ",it:=10000+it;
64end;
65cjkl:=cjkl;
66write "Update the xj";
67xj:=xj+sum(cjkl*xx(j-l)*xx(l-k)*xx(k),kl);
68
69write "But some RHS resonante and so must be in the
70evolution instead. Need definite r so set ",n:=1,
71"Explore coefficients to a max wavenumber of ",maxj:=10,
72"The negative j cases are the same with k, l, j-l and l-k
73all also changed sign. It looks like that to get dxx(j)
74correct then we need maxj=j+3, but could be different for
75n>1, guess perhaps maxj=j+3*n?";
76r:=n^2;
77o:=1; % remove order symbol
78factor xx;
79% to get substitution have to change variables!
80rhsjkl:=(rhsjkl where { b(l,k)=>kk*(ll-kk)
81 , b(j,l)=>ll*(jj-ll)
82 , d(l,k)=>1/(r+2llkk-2*kk^2) })$
83% Now sum over all contributions to dXj/dt
84array dxx(maxj);
85for j:=0:maxj do write dxx(j):= (r-j^2)*xx(j)+
86 for k:=-maxj:maxj sum for l:=-maxj:maxj sum
87 if n^2+jl+kl-k^2-l^2=0
88 then sub({jj=j,kk=k,ll=l},rhsjkl)*xx(j-l)*xx(l-k)*xx(k)
89 else 0;
90
91end;
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] R. M. Aron and J. B. Prolla. Polynomial approximation of differentiable functions on Banach spaces. J. Reine Angew. Math. , 313:195–216, 1980.
- 2[2] Bernd Aulbach, Martin Rasmussen, and Stefan Siegmund. Invariant manifolds as pullback attractors of nonautonomous differential equations. Discrete and Continuous Dynamical Systems , 15(2):579–596, 2006.
- 3[3] Bernd Aulbach and Thomas Wanner. Integral manifolds for Caratheodory type differential equations in Banach spaces. In B. Aulbach and F. Colonius, editors, Six Lectures on Dynamical Systems , pages 45–119. World Scientific, Singapore, 1996.
- 4[4] Bernd Aulbach and Thomas Wanner. Invariant foliations for Caratheodory type differential equations in Banach spaces. In V. Lakshmikantham and A. A. Martynyuk, editors, Advances of Stability Theory at the End of XX Century . Gordon & Breach Publishers, 1999. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.45.5229&rep=rep 1&type=pdf .
- 5[5] Bernd Aulbach and Thomas Wanner. The Hartman–Grobman theorem for Caratheodory-type differential equations in Banach spaces. Nonlinear Analysis , 40:91–104, 2000.
- 6[6] Luis Barreira and Claudia Valls. Stability of Nonautonomous Differential Equations , volume 1926 of Lecture Notes in Mathematics . Springer, 2007.
- 7[7] J. E. Bunder and A. J. Roberts. Nonlinear emergent macroscale PD Es, with error bound, for nonlinear microscale systems. Technical report, [ https://arxiv.org/abs/1806.10297 ], June 2018.
- 8[8] V. Bykov and V. Gol’dshtein. Fast and slow invariant manifolds in chemical kinetics. Computers & Mathematics with Applications , 2013.
