---
title: Carleman Linearization Overview
url: https://www.emergentmind.com/topics/carleman-linearization
type: topic
---

# Carleman Linearization Overview

Carleman linearization, also called Carleman embedding, is a procedure that replaces a nonlinear dynamical system by an infinite-dimensional linear system built from lifted observables of the original state. In the classical formulation, the lift consists of monomials or Kronecker powers of the state; in more recent formulations, the same embedding logic is adapted to Fourier modes, tensor powers on Hilbert spaces, shifted coordinates, and piecewise local charts. The method is exact only at the level of the infinite hierarchy. Practical use therefore centers on finite truncations, their error bounds, and the structural conditions under which the first lifted block approximates the original nonlinear trajectory with controlled accuracy [2207.07755] [2503.01498].

## 1. Classical monomial embedding

For polynomial or analytic ODEs near an equilibrium, Carleman linearization starts from a representation such as
\[
x'(t)=F_1x+F_2x^{[2]}+\cdots+F_kx^{[k]},
\]
or, in scalar form,
\[
\dot x=g(x)=\sum_{m\ge 0} c_m x^m.
\]
The state is lifted to monomials or Kronecker powers,
\[
y_i:=x^{[i]},\qquad i\in\mathbb N,
\]
so that differentiation and Leibniz’ rule yield a linear hierarchy
\[
y_i'(t)=\sum_{j=0}^{k-1} A^i_{i+j} y_{i+j}(t).
\]
Stacking the blocks produces an infinite block upper-triangular linear system. In the scalar Maclaurin setting, the lifted vector is
\[
\mathbf{x}=[x,x^2,x^3,\dots]^\mathsf T,
\]
and the resulting infinite matrix is built directly from the coefficients \(c_m\) [1711.02552] [2503.01498].

The structural reason the method is useful is that degree is preserved by the linear part and raised in a controlled way by the nonlinear terms. In the homogeneous case \(f(0)=0\), the lifted operator is block upper triangular, and the original state appears in the first block. A truncation at order \(N\) keeps only finitely many lifted coordinates and defines the finite-section approximation. In the analytic ODE framework of [2207.07755], the truncated system has dimension
\[
\binom{N+d}{d}-1,
\]
while in quadratic reachability formulations the lifted dimension scales as
\[
n+n^2+\cdots+n^N=\frac{n^{N+1}-n}{n-1}.
\]
This tradeoff—linearity in lifted space versus rapid growth of lifted dimension—is the defining computational feature of the method [2207.07755] [2108.10390].

A common misconception is that Carleman linearization yields a finite exact linear model. The exact equivalence is with the infinite hierarchy; truncation is the approximation step. The resulting linear surrogate can be highly effective, but its validity depends on truncation order, time horizon, and the analytic or spectral assumptions used in the convergence theory [1711.02552] [2207.07755].

## 2. Finite-section approximations and convergence theory

A central problem is to determine when the first block of a finite truncation converges to the true nonlinear trajectory as the truncation order increases. Early explicit results for polynomial ODEs provided local, time-dependent truncation-error bounds. For quadratic systems, one bound in [1711.02552] uses an a priori estimate \(\alpha \ge \sup_{\tau\in[0,t]}\|x(\tau)\|\), while a second, generating-function-based bound depends only on \(\|x_0\|\), \(\|F_1\|\), and \(\|F_2\|\). The paper emphasizes that the second strategy is effectively computable and yields a closed-form convergence radius \(T^*\) [1711.02552].

For analytic time-varying systems, [2207.07755] establishes explicit finite-section bounds under a uniform exponential decay assumption on the Maclaurin coefficients. If
\[
0<\|\mathbf x_0\|_\infty<\frac{R}{e},
\]
then the first block of the finite-section approximation converges exponentially in truncation order on an explicit finite interval \([t_0,t_0+T^*]\). Under an additional Hurwitz-type Jacobian condition at the origin, the same paper proves exponential convergence uniformly for all \(t\ge t_0\). This whole-time result is one of the main distinctions between merely local convergence theory and asymptotic stability-based theory [2207.07755].

Subsequent work weakens the classical dissipativity assumptions. In the quantum-simulation setting of [2405.12714], the system
\[
\frac{d}{dt}\mathbf x = F_1\mathbf x + F_2\,\mathbf x\otimes\mathbf x
\]
is analyzed under spectral non-expansion and a no-resonance condition
\[
\lambda_i \neq \sum_{j=1}^n m_j\lambda_j,\qquad \sum_j m_j\ge 2,
\]
rather than strict dissipation. The resulting estimate
\[
\|\mathbf{x}(T)-\mathbf{y}_1(T)\|_1 \le N\,C\,T\,\mathsf R_r^{\,N-1}
\]
shows convergence controlled by the resonance gap \(\Delta\). This broadens the regime from strongly dissipative systems to oscillatory and dispersive examples such as KdV and FPU chains [2405.12714].

A further refinement appears in the observable-specific analysis of [2606.24473]. Instead of bounding the entire retained tail uniformly, that paper works directly in the original monomial basis and tracks how many degree-raising nonlinear interactions are needed before discarded monomials can influence a chosen observable, such as a state coordinate. The resulting bounds are degree-aware, retain logarithmic-norm information from the original linearization, and yield geometric convergence over certified time horizons. This makes the truncation analysis sensitive to the actual degree path from the observable to the omitted tail, rather than to a worst-case global norm estimate [2606.24473].

## 3. Basis changes, shifting, and globalization

The monomial basis is natural near an equilibrium when the vector field is polynomial-like, but it can be poorly matched to periodic or quasi-periodic structure. Carleman-Fourier linearization addresses this by replacing monomials with Fourier observables. For periodic vector fields,
\[
g(x)=\sum_{n\in\mathbb Z} g_n e^{inx},
\]
the lifted variables are exponential modes rather than powers. In the multi-dimensional periodic setting of [2411.11598], the resulting operator is block upper triangular under an analyticity condition on the Fourier coefficients, and finite truncations approximate \(e^{i\mathbf x(t)}\) rather than \(\mathbf x(t)\) directly. The 2025 quasi-periodic extension introduces an extended state \(\tilde{\mathbf x}\) and block variables built from exponentials \(e^{i\boldsymbol\gamma^\mathsf T\tilde{\mathbf x}}\), recovering a block-upper-triangular hierarchy for multiple fundamental frequencies and proving exponential convergence of the primary block [2411.11598] [2503.01498].

This basis change is not merely cosmetic. The Fourier lift removes the classical dependence on conditions such as \(g(0)=0\) and small initial data that arise in equilibrium-centered monomial theory. In the Kuramoto example of [2503.01498], the first block converges exponentially to \([e^{i\theta_1(t)},e^{-i\theta_1(t)}]^\mathsf T\), and the approximation remains accurate over a broader range of initial phases and frequencies than the classical monomial truncation [2503.01498].

A different strategy is to change coordinates before lifting. The pivot-shifted framework of [2605.20071] replaces \(x\) by \(u=x-s\), then applies a Lyapunov transform \(v=Qu\), and only then performs the Carleman lift. For quadratic ODEs
\[
\partial_t x = F_2 x^{\otimes 2} + F_1 x + F_0,
\]
the shifted coefficients \(F_{1,s}\) and \(F_{0,s}\) can move the dynamics into a regime where the shifted linear part is stable. The paper proves long-time convergence when \(\alpha(F_{1,s})<0\) and derives short-time convergence guarantees even for unstable shifted systems, while removing the conventional lower-bound requirement on the initial condition [2605.20071].

Globalization methods attack a different weakness: the locality of a single expansion chart. The piecewise methods in [2510.15715] reconstruct the Carleman embedding whenever the trajectory reaches the boundary of the current chart, either with fixed chart radius, adaptively changing radius, or with a static precomputed grid. The paper’s central claim is that standard Carleman embedding fails in regions where there are multiple fixed points, whereas piecewise charts make the method usable on systems ranging from multi-equilibrium one-dimensional flows to strange attractors [2510.15715].

At the implementation level, [2605.05901] studies duplicate-aware shift-and-lift assembly. By using symmetry-reduced monomial bases, packed exponent-key indexing, and sparse triplet coalescing, the lifted affine model
\[
\dot z = A_{ZZ} z + b_Z
\]
is assembled without proliferating duplicated monomial contributions. The same paper couples this assembly strategy to a moving-center expansion, so that shift and lift are updated jointly along the trajectory [2605.05901].

## 4. Infinite-dimensional and PDE formulations

Carleman linearization extends beyond finite-dimensional ODEs. A direct PDE generalization is given in [2412.00014] for systems with quadratic nonlinearities. Instead of discrete monomial vectors, the lifted objects are tensor products of fields evaluated at independent copies of the spatial variable,
\[
\mathbf y_i := \mathbf u(\mathbf x_1)\otimes\cdots\otimes \mathbf u(\mathbf x_i).
\]
Differentiation yields a tridiagonal infinite linear PDE hierarchy in which each level couples only to levels \(i-1\), \(i\), and \(i+1\). The paper emphasizes that truncation occurs in the Carleman hierarchy rather than by spatial discretization, and illustrates the construction on Burgers’ equation and a Vlasov-type system [2412.00014].

A more functional-analytic treatment is developed in [2510.00722] for semilinear parabolic PDEs written as
\[
y'(t)+Ay(t)+B(y(t)\otimes y(t))=f(t).
\]
The lifted moments \(y^{(k)}(t)=y(t)^{\otimes k}\) satisfy a linear hierarchy with Kronecker-sum lifts of \(A\), \(B\), and \(f\). The paper proves well-posedness of the truncated hierarchy in a Gelfand triple framework and, crucially, shows that after discretization the total approximation error decomposes into two independent components: the discretization error and the linearization error. This separation motivates structure-exploiting discretizations such as sparse grids and avoids mesh-dependent deterioration of the Carleman error estimate [2510.00722].

Semigroup theory provides an even broader framework. In [2605.03381], nonlinear evolution on a Hilbert space is embedded into a linear semigroup on
\[
\Fock(H):=\bigoplus_{n\ge 1} H^n,
\qquad
E(\phi):=\bigoplus_{n\ge 1}\phi^n.
\]
The Carleman operator is analyzed through dissipativity, Lumer–Phillips, and Trotter–Kato approximation. This replaces norm-ratio convergence conditions by semigroup generation and approximation conditions, and the paper treats unbounded operators through closability and integrated semigroup theory, with hyperviscous Burgers as a detailed example [2605.03381].

The phrase “Carleman-based linearization” also appears in a related but non-identical PDE usage. The method of [2108.07914] repeatedly linearizes an over-determined quasilinear elliptic PDE around the current iterate and solves the linearized subproblem via a Carleman-weighted quasi-reversibility functional. Its relation to classical Carleman linearization is conceptual rather than literal: the convergence mechanism comes from a Carleman estimate and successive linearization, not from monomial lifting [2108.07914].

## 5. Applied uses of the lifted linear model

Carleman linearization is used less as an end in itself than as an enabling representation for analysis, control, approximation, and computation.

| Area | Role of the lift | Representative paper |
|---|---|---|
| Nonlinear control | Bilinear approximation and Volterra kernels | [2101.00495] |
| Reachability | High-dimensional linear set propagation | [2108.10390] |
| System identification | Data-driven estimation of truncated lifted dynamics | [2205.15009] |
| Chemical kinetics | Linearization for stiff reaction ODEs and HHL-compatible updates | [2207.01818] |
| Polynomial recurrences | Infinite transition matrix for nonlinear recursions | [2109.13165] |
| Fluid steady states | Second-order truncation as steady-state surrogate | [2605.23380] |

In nonlinear control, [2101.00495] uses Carleman linearization to derive a bilinear model
\[
\dot z = Az + Nzu + bu,\qquad y=c^Tz,
\]
then computes Volterra kernels up to third order and builds an IMC-Volterra controller. On the van de Vusse reactor benchmark, the third-order controller is reported to achieve the lowest control effort, with about a \(7.02\%\) reduction in ISCI relative to the linear controller [2101.00495].

In reachability analysis, [2108.10390] lifts weakly nonlinear quadratic systems to a truncated linear ODE and propagates sets with support functions. Under weak nonlinearity
\[
R := \frac{\|x_0\|\,\|F_2\|}{|\Re(\lambda_1)|} < 1
\]
and dissipativity \(\Re(\lambda_1)<0\), the truncation error obeys
\[
\|\eta_1(t)\| \le \|x_0\|\,R^N \bigl(1-e^{\Re(\lambda_1)t}\bigr)^N,
\]
which yields a dense-time overapproximation by bloating the projected lifted reach set with an error ball [2108.10390].

In system identification, [2205.15009] uses measured trajectories to estimate the truncated lifted matrix directly by least squares,
\[
\hat A = \big(\Gamma(T)-\Gamma(0)\big)\,I_\Gamma^\dagger,
\]
and selects the truncation order from a prescribed tolerance \(\Delta\). For the Van der Pol example, the framework uses 209 trajectories sampled over \([0,10]\) s and produces a certified finite-horizon trajectory bound of the form
\[
\|x(t)-\hat z|_d(t)\| \le D\mu^N + t\,\overline B\,\overline z\,\overline A
\]
on \([0,\tau^*]\) [2205.15009].

In chemical kinetics, [2207.01818] lifts polynomial reaction ODEs to a linear system \(dX/dt=A_cX\), truncates at order \(n_t\), and discretizes implicitly:
\[
(I-\Delta t\,A_N)X_{N+1}=X_N.
\]
The paper reports that increasing truncation order improves accuracy even with larger time steps, while the lifted dimension grows rapidly; in the H\(_2\)/air case, increasing the truncation order by 1 increases the matrix size by about a factor of 9 [2207.01818].

Carleman linearization also applies to nonlinear recurrences. In [2109.13165], polynomial recurrences are embedded into infinite transition matrices whose powers yield explicit solution formulas once the matrix is shifted to an upper-triangular and diagonalizable form. The method covers uni-variable depth-one recurrences, multi-variable depth-one systems, and arbitrary finite-depth systems after auxiliary-variable reduction [2109.13165].

For dissipative fluid problems, [2605.23380] studies the lowest nontrivial truncation \(K=2\), denoted C2. The paper’s main claim is that C2 can recover the late-time steady state even when it misses transient detail. This is shown analytically for a forced logistic equation and numerically for moderate-Re two-dimensional flows, suggesting that very low-order Carleman truncations may still preserve asymptotic structure [2605.23380].

## 6. Quantum formulations, computational limits, and current outlook

A major contemporary motivation for Carleman linearization is that quantum algorithms are much better developed for linear systems than for nonlinear ones. In chemical kinetics, [2207.01818] uses the HHL motivation directly: nonlinear reaction ODEs are first linearized, then rewritten as linear systems of the form
\[
(I-\Delta t\,A_N)X_{N+1}=X_N.
\]
In the broader quantum-simulation setting of [2405.12714], the truncation level needed for error \(\epsilon\) scales as
\[
N=\mathcal O\!\left(\frac{\log(T/\epsilon)}{\log(1/\mathsf R_r)}\right)
\]
under a non-resonance condition, which is the key route to quantum advantage beyond strictly dissipative dynamics [2405.12714].

Quantum implementations also expose structural constraints of the lift itself. The local Carleman linearization algorithm for quantum lattice Boltzmann dynamics in [2511.13072] uses second-order lifting with variables \(f\) and \(g\), reorganizes the encoding so that the collision remains local in the quantum registers, and reports per-time-step cost
\[
O(\log_2^3(N)+Q^4)
\]
with a correct-branch probability of order \(10^{-2}\). This is a concrete illustration that truncation order, encoding locality, and nonunitary embedding overhead are inseparable in quantum Carleman algorithms [2511.13072].

The semigroup and pivot-shifted lines of work suggest two distinct strategies for enlarging the class of tractable systems. The semigroup approach of [2605.03381] replaces fragile norm-ratio arguments by dissipativity and Trotter–Kato approximation, while the pivot-shifted framework of [2605.20071] stabilizes the coordinates before lifting and proves logarithmic truncation-order dependence on simulation time and precision in the stable shifted regime. The latter is designed specifically for preparing a quantum state proportional to the final solution of a quadratic ODE [2605.03381] [2605.20071].

The limitations remain explicit in the recent literature. Truncation error is unavoidable; higher-order truncations improve accuracy but rapidly increase dimension; local convergence can fail near multiple fixed points; closure and conditioning can dominate performance in high-order or moving-center schemes; and low-order models may preserve steady states while missing transient detail. This suggests that the current frontier is not a single universal Carleman theory, but a family of embeddings whose usefulness depends on how well the lifted basis, coordinate choice, and truncation scheme match the structure of the nonlinear system under study [2510.15715] [2605.05901] [2605.23380].

Source: https://www.emergentmind.com/topics/carleman-linearization