---
title: Probabilistic Carleman Linearization (PCL)
url: https://www.emergentmind.com/topics/probabilistic-carleman-linearization-pcl
type: topic
---

# Probabilistic Carleman Linearization (PCL)

Searching arXiv for recent papers on probabilistic Carleman linearization and closely related Carleman methods.
Probabilistic Carleman Linearization (PCL) denotes a class of methods that extend Carleman linearization from deterministic polynomial dynamics to uncertainty-aware settings by lifting nonlinear dynamics into linear systems over Kronecker powers and then propagating either moments or high-probability truncation bounds. In the discrete-time stochastic polynomial literature, PCL is the transformation of a stochastic polynomial system into an infinite-dimensional deterministic linear moment system, followed by truncation and probabilistic safety analysis [2201.08648]. In work on nonlinear stochastic differential equations driven by Ornstein–Uhlenbeck noise, the same term refers to Carleman truncation equipped with probabilistic error control rather than purely worst-case bounds [2603.12398]. Across these variants, the central objective is unchanged: replace nonlinear evolution by lifted linear dynamics while retaining quantitatively controlled approximation error.

## 1. Scope, problem classes, and relation to classical Carleman lifting

Classical Carleman linearization starts from polynomial dynamics and introduces lifted coordinates built from tensor or Kronecker powers of the state. For quadratic ordinary differential equations of the form
\[
\frac{d x(t)}{dt} = F_1 x + F_2 x^{\otimes 2},
\]
the lift
\[
\hat y_j := x^{\otimes j}, \qquad j\in \mathbb{N},
\]
produces an infinite-dimensional linear system, which is then truncated at order \(N\) and solved as a finite linear ODE in the lifted variables [2108.10390]. Explicit truncation-error analysis for such polynomial systems was developed in deterministic form by Forets and Pouly, who studied polynomial ODEs with equilibrium at the origin,
\[
x'(t)=f(x(t))=F_1x+F_2x^{[2]}+\cdots+F_kx^{[k]},
\]
and derived explicit time-dependent error bounds for finite Carleman truncations [1711.02552].

Within the probabilistic literature considered here, PCL appears in two closely related senses. The first is moment propagation for discrete-time stochastic polynomial systems with random coefficients, where expectation converts the lifted stochastic recursion into a deterministic linear recursion for raw Kronecker moments [1911.12683]. The second is probabilistic truncation control for nonlinear stochastic differential equations, where the lifted linear system remains random and the truncation error is treated as a random variable controlled by concentration arguments [2603.12398].

This dual usage matters conceptually. In the discrete-time moment-propagation setting, probabilistic structure enters through expectations and safety probabilities. In the stochastic-differential setting, probabilistic structure enters through tail bounds on the truncation residual. A plausible implication is that PCL is better viewed as a methodology family than as a single canonical algorithm.

## 2. Lifted moment dynamics for discrete-time stochastic polynomial systems

The discrete-time PCL framework studies stochastic polynomial dynamics of the form
\[
x(t+1)=\sum_{i=0}^{d} F_i(t)\,x^{[i]}(t),
\]
where \(x(t)\in \mathbb{R}^n\), the \(F_i(t)\) are random matrix-valued processes, \(F_0(t)\) represents additive-noise structure, and \(F_1(t),\dots,F_d(t)\) encode multiplicative random polynomial coefficients [2201.08648]. The method assumes that the initial condition and all coefficient processes are mutually independent, and that the coefficient processes are identically distributed over time [2201.08648].

The lift is constructed from Kronecker powers. For truncation level \(k\), the lifted vector is
\[
y_k(t)= \begin{bmatrix} x^{[0]}(t)^\top & x^{[1]}(t)^\top & \cdots & x^{[k]}(t)^\top \end{bmatrix}^\top.
\]
The combinatorial expansion of
\[
x^{[j]}(t+1)=\big(F_t\,y_d(t)\big)^{[j]}
\]
yields block matrices
\[
A_{j,k}(t)\triangleq \sum_{(i_1,\ldots,i_j)\in H_{j,k}} F_{i_1}(t)\otimes\cdots\otimes F_{i_j}(t),
\]
with the sparsity property
\[
A_{j,k}(t)=0 \quad \text{if } k>jd.
\]
Consequently, the lifted system is an infinite-dimensional linear stochastic system in the monomial coordinates [2201.08648].

Taking expectation and using independence between the current lifted coefficient block and the past lifted state produces an exact deterministic moment recursion,
\[
\mathbb{E}[x^{[j]}(t+1)] = \sum_{k=0}^{jd} E_{j,k}\,\mathbb{E}[x^{[k]}(t)],
\]
where
\[
E_{j,k}\triangleq \mathbb{E}[A_{j,k}(t)].
\]
This is the central structural step of discrete-time PCL: nonlinear stochastic dynamics are recast as an infinite deterministic linear system for moments [2201.08648].

The practical approximation truncates at order \(\ell\), defining
\[
\tilde y(t)= \begin{bmatrix} 1 & \tilde x_1(t)^\top & \tilde x_2(t)^\top & \cdots & \tilde x_\ell(t)^\top \end{bmatrix}^\top,
\]
which evolves according to the finite-dimensional linear time-invariant system
\[
\tilde y(t+1)=E_{\ell,\ell}\,\tilde y(t).
\]
This is the form used for fast offline/online propagation in the safety-analysis papers [2201.08648].

## 3. Truncation, exactness, and computational structure

The defining approximation question in PCL is how much accuracy is lost when the infinite lifted hierarchy is cut at finite order. For the \(j_0\)-th moment, the discrete-time moment error is
\[
e_{j_0}(t)\triangleq \mathbb{E}[x^{[j_0]}(t)]-\tilde x_{j_0}(t).
\]
The 2022 moment-propagation paper derives an exact expansion in which the error consists only of contributions from missing higher-order moments beyond the truncation limit, and then proves the global bound
\[
\|e_{j_0}(t)\|_\infty \le \xi(t)\sum_{j=0}^{j_0 d^t}\|\widetilde E_j\|,
\]
with
\[
\xi(t)=\max_{0\le j\le j_0 d^t}\left\|\mathbb{E}[x_{\mathrm{ini}}^{[j]}]\right\|_\infty.
\]
It also gives tighter partially exact bounds based on selected index subsets \(J\) or coordinate subsets \(K\), which are used for efficient online error estimation [2201.08648].

A distinctive feature of the discrete-time theory is exact finite-horizon recovery for sufficiently large truncation. The paper states that
\[
\text{If } j_0 d^t \le \ell,\ \text{then}\  \tilde x_{j_0}(t)=\mathbb{E}[x^{[j_0]}(t)].
\]
This sharply distinguishes the discrete-time moment framework from continuous-time Carleman truncation, where exact recovery generally does not occur [2201.08648]. The earlier 2019 paper established the same general phenomenon: sufficiently large truncation limits can precisely compute moments for sufficiently small degrees and numbers of time steps [1911.12683].

Because the ambient dimension of unreduced Kronecker powers grows rapidly, the 2022 formulation introduces reduced Kronecker powers. For \(x=[x_1,\dots,x_n]^\top\), the reduced \(m\)-th Kronecker power keeps only distinct monomials indexed by multi-indices of total degree \(m\). The paper illustrates the reduction with
\[
x^{2}=\begin{bmatrix}a^2 & ab & b^2\end{bmatrix}^\top
\]
in place of the unreduced
\[
x^{[2]}=\begin{bmatrix}a^2 & ab & ba & b^2\end{bmatrix}^\top.
\]
This reduces offline matrix construction time, memory use, and online propagation cost [2201.08648].

The same computational logic already appeared in the 2019 predecessor, which organized the method as an iterative propagation of approximate raw Kronecker moments through a truncated deterministic linear system and emphasized that the main cost driver is the rapid growth of lifted dimension with truncation order and state dimension [1911.12683].

## 4. Probabilistic safety analysis and certified regions

In the discrete-time PCL literature, propagated moments are used not only as descriptive statistics but as ingredients for safety certification. The 2022 paper defines an ellipsoid by a positive semidefinite matrix \(P\) and seminorm
\[
\|z\|_P=\sqrt{z^\top Pz},
\]
with region
\[
\{x:\|x-c\|_P\le \alpha\}.
\]
The central quantity is the deviation of the state from the propagated mean proxy \(\tilde x_1(t)\), together with a global mean-approximation error bound of the form
\[
\|\tilde x_1(t)-\mathbb{E}[x(t)]\|_P \le P,t.
\]
Markov’s inequality then yields a probabilistic bound of the form
\[
\Pr\big(\|x(t)-\tilde x_1(t)\|_P\ge \alpha\big)
\le
\frac{\sum_{i,j=1}^n p_{ij}\big(\mathbb{E}[x_i(t)x_j(t)]-\mathbb{E}[x_i(t)]\mathbb{E}[x_j(t)]\big)}
{(\alpha-P,t)^2},
\]
and, when first and second moments are sufficiently accurate, a covariance-based surrogate using \(\tilde x_1(t)\) and \(\tilde x_2(t)\) [2201.08648].

The same paper formulates online safe-ellipsoid construction as an approximate three-step optimization procedure. First, it solves a convex surrogate
\[
\max_Q \det(Q)
\]
subject to a chance-bound constraint expressed through propagated second moments. Second, it scales \(Q\) by a scalar \(s>0\) using a moment-error over-approximation. Third, it enlarges the resulting ellipsoid by the mean-error bound to obtain a certified region that contains the state with probability at least \(1-b\) [2201.08648]. This yields a fully moment-based safety pipeline with an explicit offline/online split.

The earlier 2019 paper used the same moment machinery for safety analysis, though with Euclidean-ball regions centered at the approximate mean \(\widetilde{x}_1(t)\) and tail bounds derived from Chebyshev’s inequality and second moments [1911.12683]. In both formulations, the safety certificate is only as strong as the moment approximation and its truncation-error bound. This suggests that error calibration, rather than lifting alone, is the decisive element in probabilistic use of Carleman methods.

The example systems in this line of work are a stochastic logistic map and a vehicle dynamics model under stochastic disturbance. The logistic example uses
\[
x(t+1)=r(t)x(t)(1-x(t)),
\]
with \(r(t)\sim \mathrm{Unif}[0.4,0.6]\) in the 2022 paper, while the vehicle example uses a polynomialized and discretized kinematic bicycle model with stochastic acceleration and demonstrates that larger truncation limits produce moment trajectories closer to Monte Carlo means [2201.08648].

## 5. Continuous-time stochastic PCL and high-probability truncation error

A different PCL construction appears in the quantum-simulation paper on nonlinear stochastic differential equations driven by Ornstein–Uhlenbeck noise. The system is
\[
\frac{d\mathbf{x}}{dt} = \mathbf{F}_2 \mathbf{x}^{\otimes 2} + \mathbf{F}_1 \mathbf{x} + \mathbf{F}_0(t),
\]
where \(\mathbf{F}_0(t)\) is an \(n\)-dimensional OU process satisfying
\[
d\mathbf{F}_0(t) = -\boldsymbol{\Theta}\mathbf{F}_0(t)\,dt + \boldsymbol{\Sigma}\,d\mathbf{W}(t),
\qquad \mathrm{Re}(\lambda_i(\boldsymbol{\Theta}))>0.
\]
Here the Carleman lift again uses tensor powers,
\[
\mathbf{y}_j = \mathbf{x}^{\otimes j},\qquad j=1,\dots,N,
\]
and truncation yields a finite linear stochastic system in the lifted variables [2603.12398].

The paper characterizes the truncation residual by
\[
\mathbf{R}_N(t) =
\begin{bmatrix}
0\\
\vdots\\
0\\
\mathbf{A}_{N+1}^N \mathbf{x}(t)^{\otimes (N+1)}
\end{bmatrix},
\qquad
\|\mathbf{R}_N(t)\|
\le
N\|\mathbf{F}_2\|\,\|\mathbf{x}(t)\|^{N+1}.
\]
Unlike deterministic Carleman analysis, the state norm is now random because the forcing is random. The paper therefore introduces a Lyapunov-norm framework with \(P\succ 0\), logarithmic stability parameter \(\mu_P(\mathbf{F}_1)\), one-sided quadratic growth parameter \(\beta_P\), and
\[
\kappa_P := 2\mu_P(\mathbf{F}_1)+2\beta_P+\gamma.
\]
Under these conditions it proves the pathwise estimate
\[
\|\mathbf{x}(t)\|_P^2 \le e^{\kappa_P t}\|\mathbf{x}_{\mathrm{in}}\|_P^2
+ \frac{e^{\kappa_P t}-1}{\gamma\kappa_P}
\sup_{0\le s\le t}\|\mathbf{F}_0(s)\|_P^2,
\]
which reduces truncation control to concentration of the OU supremum [2603.12398].

This is the point at which the method becomes explicitly probabilistic. By combining the pathwise state bound with Gaussian concentration for the OU process, the paper derives a high-probability truncation bound with a simplified asymptotic form
\[
\mathbb{P}\big(\|\boldsymbol{\eta}^{(N)}(t)\|_{P_N}^2 \ge \Delta_\eta\big)
\le
\exp\!\big(-c\,\Delta_\eta^{1/(N+1)}\big)
\]
for constants \(c,\Delta_{\eta,0}>0\), and, in the steady-state perturbation setting, a sharper bound
\[
\mathbb{P}\big(\|\boldsymbol{\eta}^{(N)}(t)\|_{P_N}^2 \ge \Delta_\eta\big)
\le
C_t\,\exp\!\big(-c_t[\log \Delta_\eta]^2\big).
\]
The paper interprets this as probabilistic exponential convergence for Carleman linearization, provided the nonlinear stochastic system is stable and reaches a steady state [2603.12398].

In that work, PCL is coupled to stochastic linear combination of Hamiltonian simulations (SLCHS) for quantum simulation. The quantum component is separate from the mathematical definition of PCL itself. The PCL contribution is the nonlinear-to-linear reduction together with probabilistic truncation-error control; the SLCHS contribution is the simulation of the resulting linear stochastic system [2603.12398].

## 6. Deterministic foundations, misconceptions, and limitations

PCL is not identical to the broader Carleman literature. Several prominent recent papers are deterministic counterparts or precursors rather than probabilistic methods. The reachability paper on weakly nonlinear polynomial ODEs embeds the nonlinear system into a high-dimensional linear system via Carleman linearization, propagates sets using support functions, and obtains sound dense-time reachable sets by bloating the projected linear flowpipe with a global truncation-error bound [2108.10390]. That framework already separates uncertainty in the initial condition set, deterministic lifted linear dynamics, and explicit truncation error, and the paper notes that this structure naturally suggests propagation of random initial conditions or probabilistic reachable sets. However, it also states explicitly that there are no stochastic disturbances and that the method as written is not probabilistic [2108.10390].

The same distinction applies to the 2024 quantum-analysis paper, which studies deterministic Carleman linearization beyond the dissipative regime. It proves that under a no-resonance condition, truncation error satisfies
\[
\|\mathbf x(T)-\mathbf y_1(T)\|_1 \le NCT\,\mathsf R_r^{\,N-1},
\]
with a rate controlled by the resonance gap \(\Delta\), and emphasizes that there is no randomized lifting, no Monte Carlo sampling of coefficients, no stochastic truncation rule, and no expectation-based estimator [2405.12714]. This directly addresses a common misconception: not every modern extension of Carleman truncation to broader dynamical regimes is a form of PCL.

Deterministic error theory remains foundational for PCL. Forets and Pouly derived two explicit truncation-error bounds, one requiring an a priori trajectory bound and one based only on \(\|x_0\|\), \(\|F_1\|\), and \(\|F_2\|\), thereby supplying computable worst-case envelopes for finite Carleman truncations [1711.02552]. More recent work has refined this theory by producing degree-aware observable-specific estimates via Dyson–Duhamel expansions [2606.24473], semigroup-based convergence results for unbounded generators and semi-discretized evolution equations [2605.03381], and well-posedness and convergence theorems for parabolic PDEs with a clean separation between linearization error and discretization error [2510.00722]. These works are deterministic, but they identify structural mechanisms—degree propagation, dissipativity, semigroup generation, and discretization independence—that any probabilistic extension must still respect.

The main limitations of PCL as currently represented in the literature are consistent across deterministic and probabilistic variants. The lifted dimension grows rapidly, typically as
\[
n+n^2+\cdots+n^N=\mathcal{O}(n^N)
\]
in unreduced tensor coordinates [2108.10390]. The discrete-time stochastic moment framework assumes polynomial dynamics, temporal independence, and identical distribution of the coefficient processes [2201.08648]. The continuous-time stochastic framework requires stability assumptions and is strongest when the system reaches a steady state [2603.12398]. Safety bounds may become conservative when moment-error estimates are loose [1911.12683]. Exact finite-horizon recovery is a discrete-time phenomenon and does not transfer directly to continuous-time Carleman truncation [2201.08648].

A final misconception concerns the source of probabilistic content. PCL is not, in the papers considered here, a randomized embedding or stochastic truncation heuristic. In the discrete-time literature, the probabilistic aspect lies in moment propagation and safety certification for stochastic systems [2201.08648]. In the OU-driven NSDE literature, it lies in high-probability truncation-error bounds derived from stochastic forcing concentration [2603.12398]. This suggests that the defining feature of PCL is not randomness in the lifting map itself, but the systematic integration of probability into the lifted linear approximation and its error calculus.

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