Papers
Topics
Authors
Recent
Search
2000 character limit reached

Koopman Linearizability

Updated 6 July 2026
  • Koopman Linearizability is the transformation of a nonlinear dynamical system into an exact finite-dimensional linear flow via an appropriate state-space embedding.
  • On compact invariant sets, exact linearization is characterized by torus symmetry, ensuring quasiperiodic or integrable behavior under a smooth or continuous torus action.
  • For attractors, basins, and control-affine systems, additional conditions such as continuous asymptotic phase and Lie bracket commutativity are essential for achieving global linearization.

Searching arXiv for the cited papers to ground the article in current arXiv metadata. Koopman linearizability denotes the existence of a finite-dimensional linear representation of a nonlinear dynamical system through an embedding of its state space into Euclidean space, so that the embedded dynamics evolve exactly under a linear flow. For continuous-time systems with flow Φ\Phi, the defining relation is FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F, where FF is an embedding and BB is a matrix; equivalently, the nonlinear flow is realized as the restriction of a finite-dimensional linear flow on Rn\mathbb R^n (Kvalheim et al., 2023). In the finite-dimensional Koopman viewpoint, the embedding coordinates are observables that evolve linearly under the Koopman semigroup, while in computational settings the same idea appears as a lifted linear system for observables, often only approximately (Shi et al., 2023). Recent work has clarified that exact finite-dimensional Koopman linearization is highly nongeneric: on compact invariant sets it is characterized by torus symmetry, on attractor basins it additionally requires asymptotic phase and transverse linear structure, and for control-affine systems it is constrained by Lie-bracket and distribution-rank conditions (Kvalheim et al., 2023, Deka, 11 Jun 2026).

1. Exact linearizing embeddings and the Koopman formulation

The basic object is a continuous-time dynamical system

x˙=f(x),\dot x=f(x),

with flow Φ:R×MM\Phi:\mathbb R\times M\to M, where Φt(x)\Phi^t(x) denotes the trajectory at time tt. A map F:XRnF:X\to \mathbb R^n defined on a FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F0-invariant set FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F1 is a linearizing embedding if it is an embedding in the usual topological or smooth sense and satisfies

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F2

for some matrix FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F3 (Kvalheim et al., 2023). This intertwining relation is the exact finite-dimensional Koopman representation: the nonlinear dynamics become linear in the embedded coordinates.

In Koopman terms, the state evolution is nonlinear, but observables evolve under a linear operator family. For nonlinear autonomous systems, the Koopman operator family is defined by

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F4

and its generator is

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F5

(Shi et al., 2023). Koopman eigenpairs FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F6 satisfy

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F7

which implies

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F8

(Shi et al., 2023). In this sense, exact Koopman linearizability corresponds to the existence of a finite set of observables whose span is invariant and closes under the dynamics.

The same principle extends to nonlinear infinite-dimensional systems. For an abstract evolution equation

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F9

with semiflow FF0, the Koopman semigroup on bounded continuous functionals is

FF1

and the Lie generator is

FF2

(Mauroy, 2021). This is again a linearization principle on observables rather than on state variables.

2. Compact invariant sets: torus actions as the exact criterion

For compact invariant sets, exact linearizability is characterized sharply by torus symmetry. If FF3 is a compact smooth manifold and FF4 is a smooth flow, then FF5 is linearizable by a smooth embedding if and only if FF6 is a FF7-parameter subgroup of a smooth torus action on FF8 (Kvalheim et al., 2023). Concretely, there must exist a smooth torus action

FF9

and a vector BB0 in the Lie algebra of BB1 such that

BB2

In the compact smooth case, this means that the flow must be generated by a torus action; the dynamics are quasiperiodic or integrable in a strong geometric sense (Kvalheim et al., 2023).

An analogous topological characterization holds for compact finite-dimensional metrizable spaces. If BB3 is compact finite-dimensional metrizable and BB4 is continuous, then BB5 is continuously linearizable if and only if BB6 is a BB7-parameter subgroup of a continuous torus action on BB8 with finitely many orbit types (Kvalheim et al., 2023). The finiteness of orbit types excludes pathological actions in the purely topological setting.

The compact-case proof uses equivariant embedding theorems of Mostow–Palais. If a torus action exists, one can embed equivariantly into BB9, producing a linear flow. Conversely, if the flow is already linearly embedded in Rn\mathbb R^n0, compactness forces the relevant linear dynamics to lie in the purely imaginary spectral part, and the closure of the one-parameter subgroup is a torus (Kvalheim et al., 2023). After a change of basis, the compact linear dynamics reduce to

Rn\mathbb R^n1

on Rn\mathbb R^n2, with Rn\mathbb R^n3; the closure of Rn\mathbb R^n4 is then a torus Rn\mathbb R^n5 (Kvalheim et al., 2023).

A useful recognition criterion is also available. If there exists a smooth map

Rn\mathbb R^n6

and a frequency vector Rn\mathbb R^n7 with rationally independent components such that

Rn\mathbb R^n8

then Rn\mathbb R^n9 is smoothly conjugate to a constant vector field on x˙=f(x),\dot x=f(x),0, hence smoothly linearizable (Kvalheim et al., 2023). The same result is equivalent to having x˙=f(x),\dot x=f(x),1 smooth nowhere-zero Koopman eigenfunctions with rationally independent imaginary parts. This suggests that a finite set of Koopman eigenfunctions can certify exact linearizability, but only in a very restrictive torus-based setting.

3. Attractors and basins: asymptotic phase and transverse linear structure

For noncompact invariant sets that arise as basins of attraction, the conditions are stronger. Let x˙=f(x),\dot x=f(x),2 be the basin of attraction of a globally asymptotically stable compact invariant set x˙=f(x),\dot x=f(x),3, with x˙=f(x),\dot x=f(x),4 finite-dimensional, locally compact, and metrizable. Then x˙=f(x),\dot x=f(x),5 is linearizable by a topological embedding if and only if two conditions hold: x˙=f(x),\dot x=f(x),6 has continuous asymptotic phase, and x˙=f(x),\dot x=f(x),7 is a x˙=f(x),\dot x=f(x),8-parameter subgroup of a continuous torus action on x˙=f(x),\dot x=f(x),9 with finitely many orbit types (Kvalheim et al., 2023). Continuous asymptotic phase means there exists a continuous retraction

Φ:R×MM\Phi:\mathbb R\times M\to M0

Any such topological linearizing embedding is automatically proper (Kvalheim et al., 2023).

The smooth theorem is stricter. If Φ:R×MM\Phi:\mathbb R\times M\to M1 is a smooth manifold and Φ:R×MM\Phi:\mathbb R\times M\to M2 is a globally asymptotically stable compact invariant set, then Φ:R×MM\Phi:\mathbb R\times M\to M3 is linearizable by a smooth embedding if and only if three conditions hold: Φ:R×MM\Phi:\mathbb R\times M\to M4 is a smoothly embedded submanifold with a smooth asymptotic phase map Φ:R×MM\Phi:\mathbb R\times M\to M5; Φ:R×MM\Phi:\mathbb R\times M\to M6 is a Φ:R×MM\Phi:\mathbb R\times M\to M7-parameter subgroup of a smooth torus action on Φ:R×MM\Phi:\mathbb R\times M\to M8; and there exist Φ:R×MM\Phi:\mathbb R\times M\to M9, a matrix Φt(x)\Phi^t(x)0 with all eigenvalues having negative real part, an open neighborhood Φt(x)\Phi^t(x)1, and a smooth map

Φt(x)\Phi^t(x)2

such that

Φt(x)\Phi^t(x)3

whenever Φt(x)\Phi^t(x)4 (Kvalheim et al., 2023). In other words, smooth basin linearization requires smooth phase synchronization to the attractor, torus dynamics on the attractor, and a smooth linearization of the transverse stable directions.

The global embedding for basin problems is assembled from an embedding on the attractor, a transverse coordinate, and a time-to-impact function Φt(x)\Phi^t(x)5. One explicit form is

Φt(x)\Phi^t(x)6

and in the smooth case the intertwining relation becomes

Φt(x)\Phi^t(x)7

(Kvalheim et al., 2023). The requirement that the transverse matrix Φt(x)\Phi^t(x)8 have eigenvalues with negative real part reflects stability off the attractor.

These basin theorems extend the spirit of Hartman–Grobman and Floquet theory beyond their classical local settings. Instead of a local conjugacy near a hyperbolic equilibrium or periodic orbit, one obtains a global embedding of the entire basin into a linear flow, provided the attractor has the required torus structure, asymptotic phase, and transverse linearization (Kvalheim et al., 2023).

4. Topological obstructions, equilibria, and geometric normal forms

Exact smooth linearizability on compact manifolds imposes strong topological restrictions. If Φt(x)\Phi^t(x)9 is smoothly linearizable on a connected compact manifold and has an isolated equilibrium, then tt0 must be even-dimensional and the Hopf index at any isolated equilibrium must be tt1 (Kvalheim et al., 2023). This is obtained by reducing locally to a skew-symmetric linear vector field coming from a torus action.

Several corollaries follow. If tt2 is odd-dimensional and has at least one isolated equilibrium, then it cannot be smoothly linearizable. If a compact smooth manifold supports a smooth linearizable flow with finitely many equilibria, then

tt3

If tt4 is a compact connected surface and there are only isolated equilibria, then smooth linearizability is possible only for

tt5

(Kvalheim et al., 2023). These are checkable obstructions derived from parity, Hopf index, and Euler characteristic.

The topological compact case also admits a geometric normal form. A compact continuous flow is topologically linearizable if and only if it is topologically conjugate to a quasiperiodic pinched torus family (Kvalheim et al., 2023). This description organizes topologically linearizable compact flows as quotients of torus-like families in which some circle factors are collapsed over specified closed subsets.

For basins, asymptotic phase yields another obstruction. A necessary condition for topological linearization is the existence of a continuous asymptotic phase map. The paper gives a planar spiral-to-circle flow where the attractor circle is linearizable, but the basin is not, because asymptotic phase fails (Kvalheim et al., 2023). A plausible implication is that linearizability of the attractor alone is insufficient; the geometry of the stable foliation must also be compatible with the embedding.

5. Computational linearization: spectral approximations and comparison with Carleman methods

Computational work often uses the language of Koopman linearization more broadly, referring to finite-dimensional approximations rather than exact embeddings. In one such formulation, the nonlinear autonomous system

tt6

is lifted to a linear system for observables through a spectral discretization of the Koopman generator using Chebyshev or Gauss–Lobatto collocation and the differentiation matrix (Shi et al., 2023). For a scalar observable tt7, the lifted system is

tt8

with solution

tt9

(Shi et al., 2023).

In one dimension, the finite matrix is

F:XRnF:X\to \mathbb R^n0

In two dimensions,

F:XRnF:X\to \mathbb R^n1

and for general dimension F:XRnF:X\to \mathbb R^n2,

F:XRnF:X\to \mathbb R^n3

(Shi et al., 2023). This is explicitly a linear system for observables rather than a direct state-space linearization.

The method is compared with Carleman linearization. The relationship is precise: the Carleman matrix is a transposed finite section approximation of the Koopman operator generator in the polynomial basis (Shi et al., 2023). Carleman linearization uses polynomial or tensor-product monomials F:XRnF:X\to \mathbb R^n4, truncating the monomial hierarchy at order F:XRnF:X\to \mathbb R^n5, whereas Koopman spectral linearization uses spectral collocation and a Chebyshev–Lagrange basis on Gauss–Lobatto points (Shi et al., 2023). The reported matrix sizes are F:XRnF:X\to \mathbb R^n6 for the Koopman spectral method and F:XRnF:X\to \mathbb R^n7, approximated as F:XRnF:X\to \mathbb R^n8, for Carleman lifting (Shi et al., 2023).

The numerical experiments cover five systems: F:XRnF:X\to \mathbb R^n9, FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F00, the simple pendulum, the Lotka–Volterra model, and the Kraichnan–Orszag model (Shi et al., 2023). Reported observations are that error generally decays exponentially with truncation order FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F01; error versus radius FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F02 has a bowl-shaped or V-shaped curve; Koopman spectral linearization is typically faster than Carleman, especially in FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F03; and for non-polynomial systems Carleman performs much worse because of the polynomial or Taylor truncation issue (Shi et al., 2023). At the same time, the method has important limitations: it is built for nonlinear autonomous dynamical systems, tailored to smooth scalar observables, relies critically on choosing a suitable radius FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F04, and remains expensive in high dimension because of tensor-product growth FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F05 (Shi et al., 2023).

This computational literature therefore concerns finite-dimensional approximation of Koopman linearization, not the exact geometric notion characterized by torus actions and asymptotic phase. A common misconception is to conflate successful lifted linear prediction with exact finite-dimensional Koopman linearizability. The comparison suggests that the two notions should be kept distinct.

6. Infinite-dimensional systems and control-affine extensions

For nonlinear infinite-dimensional systems such as PDEs and integro-differential equations, the Koopman framework operates on bounded continuous functionals rather than on finite-dimensional state coordinates. The finite-dimensional approximation is obtained by compressing the Koopman semigroup or Lie generator onto a subspace

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F06

yielding

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F07

with matrix representations FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F08 and FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F09 (Mauroy, 2021). Given data pairs FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F10, the least-squares approximation is

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F11

which generalizes EDMD to infinite-dimensional systems (Mauroy, 2021). Eigenvalues of FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F12 approximate eigenvalues of FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F13, generator eigenvalues are estimated by

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F14

and Koopman eigenfunctionals are approximated by right eigenvectors of FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F15 (Mauroy, 2021).

The same framework yields a linear identification method for nonlinear PDEs. If

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F16

then with basis functionals FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F17 and FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F18, one estimates

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F19

(Mauroy, 2021). Under the stated linear-independence assumptions, the paper proves that

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F20

(Mauroy, 2021). This identifies the Lie generator as the linear object through which nonlinear infinite-dimensional dynamics can be analyzed and learned.

For control-affine nonlinear systems,

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F21

the question becomes whether there exists a diffeomorphism FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F22 into a lifted space such that

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F23

with the same input FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F24 and no input feedback transformation (Deka, 11 Jun 2026). Recent differential-geometric results show that exact finite-dimensional Koopman linearizability is highly constrained. The central distribution is

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F25

where

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F26

Necessary conditions for Koopman linearization are: constant rank of the span of FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F27 on a neighborhood, and commutativity

FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F28

for all FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F29 on that neighborhood (Deka, 11 Jun 2026). The same conditions are sufficient for a weaker notion of Koopman linearizability on a control-invariant manifold, and with the added condition that the span has dimension FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F30, they become necessary and sufficient for Koopman linearizability to a controllable linear system (Deka, 11 Jun 2026).

These controlled results show that linearizability of the autonomous drift alone is not enough. One example has the same drift as a Koopman-linearizable case but a different input direction, causing the rank condition to fail; the system is then not Koopman linearizable even though the autonomous part can be Koopman-lifted (Deka, 11 Jun 2026). Another example is feedback linearizable but not Koopman linearizable because the bracket condition fails (Deka, 11 Jun 2026). This establishes a sharp distinction between Koopman linearization and classical feedback linearization.

7. Conceptual scope and limits

Across compact flows, attractor basins, infinite-dimensional systems, and control-affine dynamics, a common theme is that exact finite-dimensional Koopman linearizability is a structural property rather than a generic approximation phenomenon. On compact invariant sets, it is equivalent to hidden torus symmetry; on basins, it requires torus dynamics on the attractor together with asymptotic phase and transverse stable linearization; on control-affine systems, it requires constant-rank and commutativity conditions on iterated Lie brackets of the control vector fields with the drift (Kvalheim et al., 2023, Deka, 11 Jun 2026).

The role of topology is central. Compactness forces the relevant spectral part to be purely imaginary in compact linearizable cases, while finite-dimensionality and metrizability are required for the topological characterization (Kvalheim et al., 2023). The role of invariant manifold theory is equally central: for smooth basins, the attractor must be a smoothly embedded invariant manifold with smooth asymptotic phase, and the paper remarks more generally that it must be a normally hyperbolic invariant manifold and even an eventually relatively FΦt=eBtFF\circ \Phi^t=e^{Bt}\circ F31-NHIM (Kvalheim et al., 2023). In the control setting, invariant manifolds reappear as the natural domains on which weaker Koopman linearizations can exist (Deka, 11 Jun 2026).

A plausible implication is that the phrase “Koopman linearizability” has two distinct usages in the literature. One usage refers to exact finite-dimensional embeddings satisfying an intertwining relation, with stringent geometric and topological prerequisites. The other refers to finite-dimensional lifted approximations of the Koopman generator or semigroup, often used for prediction, spectral estimation, or identification (Shi et al., 2023, Mauroy, 2021). The former is an exact equivalence problem; the latter is an approximation methodology. Maintaining this distinction is essential for interpreting both positive computational results and impossibility theorems.

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Koopman Linearizability.