Center manifold

From Wikipedia, the free encyclopedia
Jump to navigation Jump to search
Center (red) and unstable (green) invariant manifolds of saddle-node equilibrium point of the system .
Randomly selected points of the phase space converge exponentially to center manifold on which dynamics are slow (non exponential). Studying dynamics of the center manifold allows to determine stability of the non-hyperbolic fixed point at the origin.

In mathematics, the center manifold of an equilibrium point of a dynamical system consists of orbits whose behavior around the equilibrium point is not controlled by either the attraction of the stable manifold or the repulsion of the unstable manifold. The first step when studying equilibrium points of dynamical systems is to linearize the system. The eigenvectors corresponding to eigenvalues with negative real part form the stable eigenspace, which gives rise to the stable manifold. Similarly, eigenvalues with positive real part yield the unstable manifold.

This concludes the story if the equilibrium point is hyperbolic (i.e., all eigenvalues of the linearization have nonzero real part). However, if there are eigenvalues whose real part is zero, then these give rise to the center manifold. If the eigenvalues are precisely zero, rather than just real part being zero, then these more specifically give rise to a slow manifold. The behavior on the center (slow) manifold is generally not determined by the linearization and thus is more difficult to study.

Center manifolds play an important role in: bifurcation theory because interesting behavior takes place on the center manifold; and multiscale mathematics because the long time dynamics often are attracted to a relatively simple center manifold.


Let be a dynamical system with equilibrium point . The linearization of the system near the equilibrium point is

The Jacobian matrix defines three main subspaces:

  • the stable subspace, which is spanned by the generalized eigenvectors corresponding to the eigenvalues with ;
  • the unstable subspace, which is spanned by the generalized eigenvectors corresponding to the eigenvalues with ;
  • the center subspace, which is spanned by the generalized eigenvectors corresponding to the eigenvalues with .

Depending upon the application, other subspaces of interest include center-stable, center-unstable, sub-center, slow, and fast subspaces. These subspaces are all invariant subspaces of the linearized equation.

Corresponding to the linearized system, the nonlinear system has invariant manifolds, each consisting of sets of orbits of the nonlinear system.[1]

  • An invariant manifold tangent to the stable subspace and with the same dimension is the stable manifold.
  • The unstable manifold is of the same dimension and tangent to the unstable subspace.
  • A center manifold is of the same dimension and tangent to the center subspace. If, as is common, the eigenvalues of the center subspace are all precisely zero, rather than just real part zero, then a center manifold is often called a slow manifold.

Center manifold theorems[edit]

The center manifold existence theorem states that if the right-hand side function is ( times continuously differentiable), then at every equilibrium point there exists a neighborhood of some finite size in which there is at least one of [2]

  • a unique stable manifold,
  • a unique unstable manifold,
  • and a (not necessarily unique) center manifold.

In example applications, a nonlinear coordinate transform to a normal form can clearly separate these three manifolds.[3] A web service [1] currently undertakes the necessary computer algebra for a range of finite-dimensional systems.

In the case when the unstable manifold does not exist, center manifolds are often relevant to modelling. The center manifold emergence theorem then says that the neighborhood may be chosen so that all solutions of the system staying in the neighborhood tend exponentially quickly to some solution on the center manifold. That is, for some rate .[4] This theorem asserts that for a wide variety of initial conditions the solutions of the full system decay exponentially quickly to a solution on the relatively low dimensional center manifold.

A third theorem, the approximation theorem, asserts that if an approximate expression for such invariant manifolds, say , satisfies the differential equation for the system to residuals as , then the invariant manifold is approximated by to an error of the same order, namely .

However, some applications, such as to shear dispersion, require an infinite-dimensional center manifold.[5] The most general and powerful theory was developed by Aulbach and Wanner.[6][7][8] They addressed non-autonomous dynamical systems in infinite dimensions, with potentially infinite dimensional stable, unstable and center manifolds. Further, they usefully generalised the definition of the manifolds so that the center manifold is associated with eigenvalues such that , the stable manifold with eigenvalues , and unstable manifold with eigenvalues . They proved existence of these manifolds, and the emergence of a center manifold, via nonlinear coordinate transforms. Potzsche and Rasmussen established a corresponding approximation theorem for such infinite dimensional, non-autonomous systems.[9]

Center manifold and the analysis of nonlinear systems[edit]

As the stability of the equilibrium correlates with the "stability" of its manifolds, the existence of a center manifold brings up the question about the dynamics on the center manifold. This is analyzed by the center manifold reduction, which, in combination with some system parameter μ, leads to the concepts of bifurcations.

Correspondingly, two web services currently undertake the necessary computer algebra to construct just the center manifold for a wide range of finite-dimensional systems (provided they are in multinomial form).

  • One web service [2] constructs slow manifolds for systems which are linearly diagonalised, but which may be non-autonomous or stochastic.[10]
  • Another web service [3] constructs center manifolds for systems with general linearisation, but only for autonomous systems.[11]


The Wikipedia entry on slow manifolds gives more examples.

A simple example[edit]

Consider the system

The unstable manifold at the origin is the y axis, and the stable manifold is the trivial set {(0, 0)}. Any orbit not on the stable manifold satisfies an equation of the form for some real constant A. It follows that for any real A, we can create a center manifold by piecing together the curve for x > 0 with the negative x axis (including the origin). Moreover, all center manifolds have this potential non-uniqueness, although often the non-uniqueness only occurs in unphysical complex values of the variables.

Delay differential equations often have Hopf bifurcations[edit]

Another example shows how a center manifold models the Hopf bifurcation that occurs for parameter in the delay differential equation . Strictly, the delay makes this DE infinite-dimensional.

Fortunately, we may approximate such delays by the following trick that keeps the dimensionality finite. Define and approximate the time delayed variable, , by using the intermediaries and .

For parameter near critical, , the delay differential equation is then approximated by the system

Copying and pasting the appropriate entries, the web service [4] finds that in terms of a complex amplitude and its complex conjugate , the center manifold

and the evolution on the center manifold is

This evolution shows the origin is linearly unstable for , but the cubic nonlinearity then stabilises nearby limit cycles as in classic Hopf bifurcation.


  1. ^ Guckenheimer & Holmes (1997), Section 3.2
  2. ^ Guckenheimer & Holmes (1997), Theorem 3.2.1
  3. ^ Murdock, James (2003). Normal forms and unfoldings for local dynamical systems. Springer-Verlag.
  4. ^ Iooss, G.; Adelmeyer, M. (1992). Topics in Bifurcation Theory. p. 7.
  5. ^ Roberts, A. J. (1988). "The application of centre manifold theory to the evolution of systems which vary slowly in space". J. Austral. Math. Soc. B. 29: 480–500. doi:10.1017/S0334270000005968.
  6. ^ Aulbach, B.; Wanner, T. (1996). "Integral manifolds for Caratheodory type differential equations in Banach spaces". In Aulbach, B.; Colonius, F. (eds.). Six Lectures on Dynamical Systems. Singapore: World Scientific. pp. 45–119.
  7. ^ Aulbach, B.; Wanner, T. (1999). "Invariant foliations for Caratheodory type differential equations in Banach spaces". In Lakshmikantham, V.; Martynyuk, A. A. (eds.). Advances of Stability Theory at the End of XX Century. Gordon & Breach.
  8. ^ Aulbach, B.; Wanner, T. (2000). "The Hartman–Grobman theorem for Caratheodory-type differential equations in Banach spaces". Nonlinear Analysis. 40: 91–104. doi:10.1016/S0362-546X(00)85006-3.
  9. ^ Potzsche, C.; Rasmussen, M. (2006). "Taylor approximation of integral manifolds". Journal of Dynamics and Differential Equations. 18: 427–460. Bibcode:2006JDDE...18..427P. doi:10.1007/s10884-006-9011-8.
  10. ^ A.J. Roberts (2008). "Normal form transforms separate slow and fast modes in stochastic dynamical systems". Physica A. 387: 12–38. arXiv:math/0701623. Bibcode:2008PhyA..387...12R. doi:10.1016/j.physa.2007.08.023.
  11. ^ A.J. Roberts (1997). "Low-dimensional modelling of dynamics via computer algebra". Comput. Phys. Commun. 100: 215–230. arXiv:chao-dyn/9604012. Bibcode:1997CoPhC.100..215R. doi:10.1016/S0010-4655(96)00162-2.


External links[edit]