---
title: Isostables in Nonlinear Dynamics
url: https://www.emergentmind.com/topics/isostables
type: topic
---

# Isostables in Nonlinear Dynamics

Isostables are geometric objects in the state space of nonlinear dynamical systems, defined as the level sets of the principal Koopman eigenfunction associated with the slowest decaying mode of attraction. They provide a global foliation of the basin of attraction for both fixed points and periodic orbits, characterizing sets of initial conditions that converge synchronously to the attractor at identical asymptotic rates. The concept of isostables arises in the context of spectral operator theory and has become fundamental in phase-amplitude reduction, model reduction, numerical computation, and the analysis of transient phenomena in high-dimensional and nonlinear systems [1911.11996].

## 1. Definition and Spectral Characterization

Formally, for a $C^1$ flow $\Phi^t$ on the basin $Q$ of a hyperbolic attracting fixed point $x_0$ or periodic orbit $\Gamma$, a Koopman eigenfunction $\psi: Q \to \mathbb{C}$ with exponent $\mu \in \mathbb{C}$ satisfies
\[
\psi \circ \Phi^t(x) = e^{\mu t}\ \psi(x),\ \forall t.
\]
The principal eigenfunction is defined by $\psi(x_0) = 0$ and $D\psi(x_0) \neq 0$ in the fixed point case, or $\psi|_{\Gamma}\equiv 0$ and $D\psi(x_0)\neq 0$ for $x_0\in \Gamma$ in the periodic orbit case. The isostables are the level sets of $|\psi|$, or equivalently, of $\operatorname{Re}\ \psi$ or $\operatorname{Im}\ \psi$ for the real or complex-conjugate principal eigenfunctions. The corresponding Koopman exponent $\mu$ must yield the slowest decaying Floquet multiplier $e^\mu$, which dominates the asymptotic convergence [1911.11996].

Uniqueness and global existence of principal Koopman eigenfunctions (and thus isostable coordinates) are ensured under $k$-nonresonance conditions on the dynamics, as established in [1911.11996]. When $e^\mu$ is simple and nonresonant, any two eigenfunctions differ only by a constant scalar multiple.

## 2. Existence, Uniqueness, and Explicit Construction

The core result is that, given a hyperbolic sink $x_0$ with $D\Phi^1(x_0)$ of spectral radius $\rho<1$, for any complex $\mu$ with $|e^\mu|>\rho^{k+\alpha}$ and appropriate left eigenvector $w$, there exists a unique $\psi \in C^{k,\alpha}(Q, \mathbb{C})$ with $\psi \circ \Phi^t = e^{\mu t}\psi$ and $D\psi(x_0)=w$. Explicitly, for any “approximate” eigenfunction $P$ satisfying $D P(x_0) = w$ and an appropriate approximate eigenrelation, the limit
\[
\psi(x) = \lim_{t\to\infty} e^{-\mu t}P(\Phi^t(x))
\]
converges to the principal eigenfunction in $C^{k,\alpha}$ on compacts [1911.11996].

For systems with an attracting periodic orbit, the construction is via the Poincaré map on a strong-stable leaf, extended globally by
\[
\psi(\Phi^t(x)) = e^{tA}\psi(x),
\]
with $A$ the Floquet exponent matrix. The isostables are thus globally smooth, well-defined, and uniquely determined principal Koopman eigenfunctions.

## 3. Hierarchy of Isostables and Spectral Decomposition

Beyond the slowest mode, "faster" isostable coordinates can be defined recursively. For each stable direction with Koopman exponent $\mu_j$, the Laplace-type limit
\[
\psi_j(x) = \lim_{T\to\infty} e^{-\mu_j T} P(\Phi^T(x)),
\]
for an approximate $P$ vanishing on more strongly stable directions, converges to the unique principal eigenfunction associated to $\mu_j$ whenever $|e^{\mu_j}| > \rho^{k+\alpha}$ and nonresonance conditions hold [1911.11996]. The hierarchy of isostables thus mirrors the stable spectral decomposition of the linearized system.

This spectral hierarchy enables construction of a vector-valued factor map
\[
\Psi: Q \to \mathbb{C}^m,\quad \Psi \circ \Phi^t = e^{tA} \Psi,
\]
where the components of $\Psi$ are commuting Koopman eigenfunctions—in local coordinates, these are precisely the isostable coordinates.

## 4. Global Linearization, Pullback Algebra, and Coordinate Representations

Under the existence and uniqueness theorems, there exists a unique $C^{k,\alpha}$ semiconjugacy $\Psi$ relating the nonlinear flow to its linearization $e^{tA}$ on the basin of attraction, with the algebra $A_\Phi$ generated by the components $\{\psi_j\}$ being independent of the particular choice of conjugacy [1911.11996]. The pullback algebra, constructed from the linearized principal algebra and the nonlinear coordinate transformation, is thus robust and well-defined.

Global linearization via isostable coordinates provides a mechanism for diagonalizing the (nonlinear) dynamics and for constructing global (nonlinear) action-angle or phase-amplitude coordinates, which are essential in the reduction and analysis of transient behavior even far from the attractor.

## 5. Illustrative Examples: Linear and Nonlinear Systems

For the linear system $\dot{x} = -x$, $\dot{y} = -2y$, with $Df(0) = \mathrm{diag}(-1, -2)$, the slow and fast exponents are $\mu_1 = -1$, $\mu_2=-2$. The principal (slowest) eigenfunction is $\psi_1(x,y) = x$, and the fast isostable coordinate is $\psi_2(x,y)=y$, both recovering classical stable subspaces (isostables $x=\text{const}$ and $y=\text{const}$, respectively). For a nonlinear perturbation $\dot{x} = -x,\ \dot{y} = -2y + x^2$, the isostable $\psi_2(x,y)=y-x^2$ is recovered via the semiconjugacy to the linear system, illustrating the persistence of the isostable framework under nonlinear conjugacy [1911.11996].

The general vector-valued isostable coordinates smoothly interpolate between linear subspaces and nonlinear normal forms, providing global coordinates for the dynamics on the entire basin.

## 6. Relation to Other Notions: Isochrons, Koopmanism, and Pullback Structures

Isostables are complementary to isochrons, which are level sets of the phase function for periodic orbits (associated with the zero Floquet exponent). While isochrons foliate the basin by asymptotic phase, isostables organize it by amplitude of decay or approach rate to the attractor. The global action-angle formulation unifies these through the Koopman eigenfunction framework.

The pullback algebra of principal eigenfunctions is shown to be unique under nonresonance hypotheses, eliminating ambiguity in the spectral description and ensuring the robustness of isostable coordinates in applied and computational Koopmanism [1911.11996].

## 7. Significance and Applications

Isostables articulated through the Koopman operator theory provide a natural, global framework for:

- Analysis of transient behavior and reduction of nonlinear systems to phase-amplitude coordinates;
- Computation and control strategies for high-dimensional dynamical systems, including optimal stabilization, synchronization, and suppression of transient deviations;
- Systematic organization of basins of attraction and synchronized convergence for a wide class of dissipative systems.

The existence and uniqueness theorems guarantee rigorous and robust applicability, while explicit constructions via Laplace-type limits and semiconjugacies enable practical computation in both theoretical and engineered systems. The framework generalizes and sharpens classical results such as Sternberg linearization and Floquet theory within a unifying operator-theoretic perspective [1911.11996].

Source: https://www.emergentmind.com/topics/isostables