---
title: Koopman-Krylov Framework
url: https://www.emergentmind.com/topics/koopman-krylov-framework
type: topic
---

# Koopman-Krylov Framework

The Koopman-Krylov framework denotes a class of constructions in which a Koopman operator, or its infinitesimal generator, is approximated on a Krylov subspace generated by repeated action on a seed observable or on a sequence of snapshots. In the discrete-time setting, this yields finite-section or companion-matrix representations of the Koopman operator; in continuous-time and generator-based settings, it yields Galerkin or Lanczos reductions of the Liouville/Koopman generator. Across the cited literature, the framework appears in numerical Koopman spectral analysis, time-delay and Hankel formulations of DMD, robust estimation for stochastic systems, data-driven control of parabolic PDEs, and generator-based diagnostics of integrability breaking in semiclassical strings [2009.05883] [1808.09557] [2011.00078] [2407.00432] [2602.23421].

## 1. Operator-theoretic formulation

In its most basic discrete-time form, one considers a dynamical system \(T:M\to M\) with Koopman operator \(U\) acting on a function space \(\mathcal F\), fixes a seed observable \(f\in\mathcal F\), and forms the Krylov sequence
\[
f_1=f,\quad f_2=Uf,\quad \dots,\quad f_m=U^{m-1}f.
\]
The associated Krylov subspace is
\[
\mathcal K_m=\mathrm{span}\{f_1,\dots,f_m\},
\]
and the finite-section or Galerkin approximation is
\[
U_m=P_mUP_m:\mathcal K_m\to\mathcal K_m,
\]
where \(P_m\) is the projection onto \(\mathcal K_m\). In the ordered basis \((f_1,\dots,f_m)\), the matrix of \(U_m\) has Frobenius companion form, with a shift structure in the first \(m-1\) columns and a last column determined by the projection of \(U^mf\) onto \(\mathcal K_m\) [2009.05883].

An equivalent state-space formulation arises from a snapshot sequence
\[
x_0,x_1,\dots,x_m\in\mathbb C^n,\qquad x_{k+1}=Ax_k.
\]
Writing
\[
X_m=[x_0\;\;x_1\;\;\dots\;\;x_{m-1}],
\]
one obtains
\[
AX_m=X_mC_m+r_me_m^T,
\]
where \(C_m\in\mathbb C^{m\times m}\) is a Frobenius companion matrix and \(r_m\) is orthogonal to the Krylov subspace \(\mathrm{span}\{x_0,Ax_0,\dots,A^{m-1}x_0\}\). Algebraically,
\[
C_m=(X_m^*X_m)^{-1}(X_m^*AX_m),
\]
so \(C_m\) is the representation of \(A\) on the Krylov subspace in the nonorthonormal basis \(X_m\) [1808.09557].

This identification of Koopman approximation with a Krylov compression is also the basis of time-delay and Hankel constructions. In the stochastic setting, a Hankel matrix built from one scalar observable produces a column space approximating
\[
\mathrm{span}\{f,\mathcal K f,\dots,\mathcal K^m f\},
\]
thereby realizing a Koopman-Krylov approximation from a single observable measured over a single trajectory [2011.00078].

## 2. Spectral representation, modes, and reconstruction weights

When the companion matrix \(C_m\) has simple eigenvalues \(\{\lambda_1,\dots,\lambda_m\}\), its left eigenvectors are the rows of the Vandermonde matrix
\[
V_m=[V_{ij}]_{i,j=1}^m,\qquad V_{ij}=\lambda_i^{\,j-1},
\]
and its right eigenvectors are the columns of \(V_m^{-1}\). Hence
\[
C_m=V_m^{-1}\Lambda_mV_m,\qquad \Lambda_m=\mathrm{diag}(\lambda_1,\dots,\lambda_m).
\]
Defining the Ritz or Koopman modes by
\[
Z_m=X_mV_m^{-1}=[z_1\;\cdots\;z_m],
\]
the snapshots admit the modal reconstruction
\[
x_k=\sum_{j=1}^m z_j\,\lambda_j^{\,k-1},\qquad k=1,\dots,m.
\]
If \(Z_m=W_mD_{\mathfrak a}\) with \(\|w_j\|=1\) and amplitudes \(\mathfrak a_j\ge 0\), then
\[
x_k=\sum_{j=1}^m w_j\,\mathfrak a_j\,\lambda_j^{\,k-1}.
\]
This is the canonical finite-dimensional Koopman modal expansion in the companion-matrix realization [1808.09557].

A reduced reconstruction uses only \(\ell<m\) modes indexed by \(J\). One then seeks coefficients \(\alpha_j\) minimizing the least-squares error across all snapshots,
\[
\min_{\alpha\in\mathbb C^\ell}\sum_{k=1}^m\left\|x_{k-1}-\sum_{j\in J}z_j\,\alpha_j\,\lambda_j^{\,k-1}\right\|_2^2.
\]
If \(Z_J=QR\) and \(g_k=Q^*x_{k-1}\), the problem reduces to a structured least-squares system. The Moore-Penrose solution is \(\alpha=S^\dagger g\), but the introduction of a reflexive \(g\)-inverse \(S^{-}\) yields an explicit formula,
\[
\alpha_{*j}
=
\left(\sum_{k=1}^m|\lambda_j|^{2(k-1)}\right)^{-1}
\sum_{i=1}^m\overline{\lambda}_j^{\,i-1}\,(R^{-1}g_i)_j,
\qquad j\in J.
\]
Equivalently, these coefficients minimize a weighted norm
\[
\|x\|_M=\sqrt{x^*Mx},\qquad M=I\otimes(RR^*)^{-1}.
\]
The resulting approximation is
\[
x_{k-1}\approx \sum_{j\in J}z_j\,\alpha_{*j}\,\lambda_j^{\,k-1}
\equiv \sum_{j\in J}\psi_j\,w_j\,\lambda_j^{\,k-1},
\]
with \(\psi_j=z_j\) as Koopman modes and \(w_j=\alpha_{*j}\) as reconstruction weights [1808.09557].

The same paper makes the relation to Generalized Laplace Analysis precise. In matrix form, for a diagonalizable operator \(A=Z\Lambda Z^{-1}\), the coordinate GLA is
\[
\alpha=Z^{-1}f=\frac1m\sum_{k=1}^m\Lambda^{-(k-1)}Z^{-1}x_{k-1}.
\]
Truncating to \(\ell\) modes gives GLA reconstruction weights
\[
\alpha_{(\mathrm{GLA}),j}
=
\frac1m\sum_{k=1}^m\lambda_j^{-(k-1)}(R^{-1}g_k)_j.
\]
If all \(\lambda_j\) lie on the unit circle, the GLA weights coincide exactly with the reflexive \(g\)-inverse weights. If \(|\lambda_j|<1\), the reflexive \(g\)-inverse weights become a nonuniformly weighted ergodic average. This shows that finite-sample optimal reconstruction and GLA are not separate constructions but two limits of the same modal-weighting problem [1808.09557].

## 3. Numerical conditioning and the Vandermonde-Cauchy remedy

The central numerical obstacle in companion-form Koopman-Krylov methods is the inverse Vandermonde matrix. Direct inversion of \(V_m\) or direct solution of Vandermonde systems suffers \(\mathcal O(\kappa(V_m)\times\epsilon)\) relative error, and \(\kappa(V_m)\) can grow like \(O(m^{m+1})\). In floating-point arithmetic, directly computing \(V_m^{-1}\) often destroys the accuracy of the Koopman modes even when the Ritz values themselves are meaningful [1808.09557].

The stabilization developed in "Data driven Koopman spectral analysis in Vandermonde-Cauchy form via the DFT: numerical method and theoretical insights" [1808.09557] is based on the observation that multiplication by the unitary discrete Fourier transform \(F\) converts the Vandermonde matrix into a generalized Cauchy matrix. With
\[
F_{pq}=\frac1{\sqrt m}e^{-2\pi i(p-1)(q-1)/m},
\qquad
\zeta_j=e^{-2\pi i(j-1)/m},
\]
one has
\[
(V_mF)_{ij}
=
\frac{\lambda_i^m-1}{\sqrt m}\,
\frac{1}{\lambda_i-\zeta_j}\,
\zeta_j,
\]
or, in matrix form,
\[
V_mF=\mathcal D_1\,C\,\mathcal D_2,
\]
where \(\mathcal D_1\) and \(\mathcal D_2\) are diagonal and \(C_{ij}=\frac1{\lambda_i-\zeta_j}\) is a generalized Cauchy matrix. Equivalently,
\[
V_m^{-1}=F\,\mathcal D_2^{-1}\,C^{-1}\,\mathcal D_1^{-1}.
\]

This change of basis does not remove ill-conditioning in the spectral sense, but it relocates the problem to a structured Cauchy factorization for which specially tailored numerical linear algebra is available. In particular, a pivoted LDU factorization
\[
\Pi_1(\mathcal D_1C\mathcal D_2)\Pi_2=L\Delta U
\]
can be computed by entry-wise forward-stable algorithms, with
\[
|\widetilde L_{ij}-L_{ij}|\le\epsilon |L_{ij}|,\qquad
|\widetilde\Delta_{ii}-\Delta_{ii}|\le\epsilon |\Delta_{ii}|,\qquad
|\widetilde U_{ij}-U_{ij}|\le\epsilon |U_{ij}|,
\]
even when \(\mathrm{cond}(V_m)\gg \epsilon^{-1}\). Once these factors are available, \(V_m^{-1}\) is recovered through unitary transforms and triangular solves in \(O(m^2)\), and \(Z_m=X_mV_m^{-1}\) can be obtained to full machine precision across all entries [1808.09557].

The same stability analysis extends to the eigenvalue computation. Rather than using a black-box eigensolver that destroys companion structure, the paper advocates companion-QR or unitary-plus-rank-one QR methods. These produce eigenvalues that are exactly the roots of a nearby companion polynomial
\[
\widetilde c(z)=z^m-\sum_j(c_j+\delta c_j)z^{j-1},
\qquad \|\delta c\|\approx O(\epsilon)\|c\|,
\]
thereby preserving backward stability in the polynomial coefficients and maintaining accuracy even when Ritz values cluster [1808.09557].

## 4. Convergence theory, DMD relations, and failure modes

Theoretical analysis of Koopman-Krylov approximations in finite-section form is developed in "On Numerical Approximations of the Koopman Operator" [2009.05883]. Under the assumptions that the Koopman operator \(U\) has pure point spectrum \(\{\lambda_j\}_{j=1}^\infty\subset\{z:|z|\le 1\}\), that the corresponding normalized eigenfunctions form a Riesz basis, and that the seed observable admits an expansion
\[
f=\sum_{j=1}^\infty c_j\phi_j,
\]
the eigenpairs of the Krylov-Galerkin truncations become pseudoeigenpairs of the full operator. Concretely, if \((\tilde\lambda_m,\tilde\phi_m)\) is an eigenpair of \(U_m\) with \(\|\tilde\phi_m\|=1\), then for every \(\varepsilon>0\) there exists \(M\) such that for all \(m\ge M\),
\[
\|U\tilde\phi_m-\tilde\lambda_m\tilde\phi_m\|<\varepsilon.
\]
The convergence is therefore in the pseudospectral sense rather than as a direct pointwise spectral convergence statement [2009.05883].

The key residual identity is
\[
U\tilde\phi_m-\tilde\lambda_m\tilde\phi_m
=
e_m\bigl(U^mf-P_mU^mf\bigr),
\]
where \(e_m\) is the coefficient of \(f_m=U^{m-1}f\) in \(\tilde\phi_m\). Since \(|e_m|\le 1\), the approximation error is controlled by
\[
r_m=U^mf-P_mU^mf,
\qquad
\|U\tilde\phi_m-\tilde\lambda_m\tilde\phi_m\|\le \|r_m\|.
\]
Under the pure-point assumption, one obtains a rate function
\[
\rho(m)
=
\sum_{|\lambda_j|<1}|c_j||\lambda_j|^m
+
\inf_{g\in\mathcal K_m}
\left\|
\sum_{|\lambda_j|=1}c_j\lambda_j^m\phi_j-g
\right\|
\to 0,
\]
which separates decaying interior-unit-disc contributions from almost-periodic unit-circle contributions [2009.05883].

The same paper places DMD-type methods in the context of finite-section theory. Both Krylov-Galerkin and DMD approximate \(P_mUP_m\); the difference is the basis choice. Krylov-Galerkin uses the time-delay basis
\[
f,\;Uf,\;U^2f,\dots,
\]
whereas standard extended DMD uses a prescribed dictionary. The Krylov choice avoids the "curse of dimensionality" in the sense stated in the paper, because the subspace dimension grows with the number of iterates rather than with the ambient state dimension [2009.05883]. This also clarifies why Hankel-DMD and companion-form DMD are naturally interpreted as Koopman-Krylov constructions.

The same analysis also identifies a sharp limitation. For a mixing map such as \(z\mapsto z^2\) on the unit circle, finite-section eigenvalues do not approach the continuous Koopman spectrum, and the companion matrix collapses to a nilpotent shift. By contrast, for the rotation \(T(z)=e^{i\omega}z\) with \(f(z)=z\), the method recovers the exact eigenvalue \(e^{i\omega}\). This establishes that pure point spectrum is a favorable regime for Koopman-Krylov approximation, whereas continuous-spectrum settings can defeat finite-section convergence [2009.05883].

## 5. Stochastic Koopman-Krylov and robust estimation under noise

For random dynamical systems
\[
x_{n+1}=F(x_n,\omega_n),
\]
with i.i.d. \(\omega_n\), the stochastic Koopman operator acts on deterministic observables by
\[
(\mathcal Kf)(x)=\mathbb E_\omega[f(F(x,\omega))]
=
\int_\Omega f(F(x,\omega))\,dP(\omega).
\]
Its iterates form a semigroup \(\mathcal K^{t+s}=\mathcal K^t\circ\mathcal K^s\) because the maps are i.i.d. [2011.00078].

In this setting, standard DMD can be biased when either the dynamics or the observables are noisy. If \(\tilde f(x,\omega)\) is a noisy observable with mean \(f(x)\), the empirical covariance matrices built by time averaging along one trajectory satisfy
\[
\tilde G_j
=
G_j
+
\int_M
\mathrm{Cov}\!\bigl(\tilde{\mathbf f}_{\theta_j\omega}\circ T_\omega^j,\tilde{\mathbf f}_\omega\bigr)
\,d\mu,
\]
so the usual EDMD estimator \(\tilde G_1\tilde G_0^\dagger\) converges to a biased matrix rather than to the true finite-rank restriction of \(\mathcal K\) [2011.00078].

The robust remedy introduces a second noisy observable \(\tilde{\mathbf g}\) that is uncorrelated with \(\tilde{\mathbf f}\) and with the forward dynamics. With
\[
X=[\hat{\mathbf f}(0)\;\cdots\;\hat{\mathbf f}(N-1)],\quad
Y=[\hat{\mathbf f}(1)\;\cdots\;\hat{\mathbf f}(N)],\quad
Z=[\hat{\mathbf g}(0)\;\cdots\;\hat{\mathbf g}(N-1)],
\]
the cross-covariance matrices are
\[
\tilde G_0=\frac1N XZ^*,\qquad
\tilde G_1=\frac1N YZ^*,
\]
and the robust estimator is
\[
A_{\mathrm{robust}}=\tilde G_1\tilde G_0^\dagger.
\]
Under mild ergodicity and independence conditions, this estimator converges to the true restriction of \(\mathcal K\) on the chosen subspace [2011.00078].

The Koopman-Krylov extension uses time-delayed observables. For a scalar signal \(\{\hat f(0),\dots,\hat f(N)\}\), the Hankel matrix
\[
H_m=
\begin{bmatrix}
\hat f(0)&\hat f(1)&\cdots&\hat f(N-m)\\
\hat f(1)&\hat f(2)&\cdots&\hat f(N-m+1)\\
\vdots&\vdots&\ddots&\vdots\\
\hat f(m)&\hat f(m+1)&\cdots&\hat f(N)
\end{bmatrix}
\]
has columns approximating
\[
[f,\mathcal Kf,\dots,\mathcal K^mf]^T
\]
up to noise, so \(\mathrm{ColSpan}(H_m)\) approximates the Krylov subspace. Using the shifted Hankel matrix \(H_m'\) and a lagged dual observable \(\hat{\mathbf g}(t)=\hat{\mathbf f}(t-m)\), one obtains
\[
A_{\mathrm{robust,Krylov}}
=
\Bigl(\tfrac1{N-m}H_m'Z^*\Bigr)
\Bigl(\tfrac1{N-m}H_mZ^*\Bigr)^\dagger,
\]
which converges to the finite-dimensional companion matrix of \(\mathcal K\) on the Krylov subspace. For numerical stability, the same paper also proposes an SVD reduction
\[
H_m\approx V\Sigma W^*,
\qquad
C=\Sigma^{-1}V^*H_m'W\approx V^*\mathcal KV
\]
as Algorithm 4 [2011.00078].

The worked example of a random rotation on the circle illustrates the distinction. Standard EDMD yields eigenvalues biased toward the real axis, whereas robust DMD with a lagged dual observable recovers the analytic eigenvalues within sampling error. In the Hankel test with \(m=5\), standard Hankel DMD is highly ill-conditioned and biased, while the robust SVD-based algorithm accurately recovers six modes when enough delays are used [2011.00078].

## 6. Data-driven control of parabolic PDEs

A control-theoretic extension appears in "Data-Driven Control of Linear Parabolic Systems using Koopman Eigenstructure Assignment" [2407.00432]. The plant is a boundary-controlled linear parabolic PDE on \(z\in[0,1]\),
\[
\partial_tx(z,t)=\rho\,\partial_z^2x(z,t)+a(z)x(z,t),
\]
with boundary conditions
\[
\partial_zx(0,t)=q_0x(0,t)+u_1(t),\qquad
\partial_zx(1,t)=q_1x(1,t)+u_2(t).
\]
The state space is \(X=L_2(0,1)\), the uncontrolled generator is a self-adjoint Sturm-Liouville operator \(\mathcal A\), and the semigroup is \(\mathcal T(t)=e^{\mathcal At}\). For bounded linear functionals \(g\in X'\), the Koopman operator is
\[
(\mathscr K(t)g)[x]=g[\mathcal T(t)x].
\]
The open-loop Koopman eigenfunctionals are
\[
\varphi_i[x]=\langle x,\phi_i\rangle,
\qquad
\mathscr K(t)\varphi_i[x]=e^{\lambda_it}\varphi_i[x],
\]
where \(\mathcal A\phi_i=\lambda_i\phi_i\) and \(\{\phi_i\}\) is an orthonormal eigenbasis [2407.00432].

The paper extends classical Krylov-DMD to this infinite-dimensional setting using only finitely many point-averaged outputs sampled at times \(t_k=kt_s\),
\[
y_k=\langle x(kt_s),c\rangle\in\mathbb R^M.
\]
From the data matrix
\[
\mathbb D_{n+1}(x)=[y_0\;\;y_1\;\;\dots\;\;y_n],
\]
one seeks coefficients \(f_0,\dots,f_{n-1}\) such that
\[
y_n\approx -\sum_{j=0}^{n-1}f_jy_j,
\]
which yields a companion-DMD matrix \(F\in\mathbb R^{n\times n}\) with least-squares residual \(r_n\). Solving
\[
Fv_i=\mu_iv_i
\]
produces DMD approximations of the Koopman spectrum and modes; the errors are controlled by \(\|r_n\|\) [2407.00432].

These approximations are then used in a Koopman eigenstructure assignment problem. With approximate open-loop data \(\hat\Lambda_n\) and \(\hat B_n\), and desired closed-loop eigenvalues \(\tilde\lambda_1,\dots,\tilde\lambda_n\), the gain is parameterized through vectors \(p_i\in\mathbb C^2\) and
\[
\tilde v_i=(\hat\Lambda_n-\tilde\lambda_iI)^{-1}\hat B_np_i.
\]
If the \(\tilde v_i\) are linearly independent, the unique real feedback gain is
\[
K=[p_1\;\dots\;p_n][\tilde v_1\;\dots\;\tilde v_n]^{-1}.
\]
The resulting controller requires only a finite number of open-loop Koopman eigenvalues and modes of the state, extracted from output data and input samples rather than from an identified PDE model [2407.00432].

The stability result is correspondingly finite-dimensional in construction but infinite-dimensional in consequence. In the ideal case, the assigned finite block evolves under a Hurwitz matrix while the uncontrolled tail remains exponentially decaying. In the data-driven case, Theorem 4 shows that small DMD errors in \(\hat\Lambda_n\) and \(\hat B_n\), together with sufficiently small mode-interpolation error, preserve exponential stability of the full PDE in the \(L_2\) norm [2407.00432]. The numerical example uses
\[
\rho=1,\quad a(z)=7-8(z-0.5)^2,\quad q_0=2,\quad q_1=1,
\]
with \(M=500\) spatial points, \(t_s=0.004\), \(n+1=12\) snapshots, and residual \(\|r_{11}\|\approx 1.6\cdot 10^{-7}\); the closed-loop design shifts the unstable and slow modes to \(\{-7.00,-10.77,-52.73\}\) and yields rapid stabilization in simulation [2407.00432].

## 7. Generator-based Koopman-Krylov space and semiclassical strings

A generator-centric variant is developed in "Integrability breaking in semiclassical strings in Koopman-Krylov space" [2602.23421]. Here the starting point is the Koopman-von Neumann formulation of classical mechanics: observables \(f:X\to\mathbb C\) are elements of
\[
\mathcal H_{\mathrm{obs}}=L^2(X,\mu),
\]
the nonlinear flow \(F^\sigma\) induces the unitary Koopman operator
\[
(U^\sigma f)(x)=f(F^\sigma(x)),
\qquad
U^\sigma=e^{\sigma\mathcal L},
\]
and for a Hamiltonian \(H(q,p)\) the generator is the Liouville operator
\[
\mathcal Lf(q,p)=\{f,H\}.
\]
Writing \(K=-i\mathcal L\), the autocorrelation of a seed observable has the spectral representation
\[
C_f(\sigma)=\langle f,U^\sigma f\rangle
=
\int_{\mathbb R}e^{-i\sigma\omega}\,d\mu_f(\omega),
\]
with \(\mu_f\) the spectral measure of \((K,f)\) [2602.23421].

Since \(\mathcal L\) acts on an infinite-dimensional space, the paper uses generator-extended DMD (gEDMD) to construct a finite-dimensional Galerkin approximation. With a dictionary
\[
\Phi(x)=(\phi_1(x),\dots,\phi_N(x))^T,
\]
empirical Gram and generator matrices are estimated from sampled phase-space points,
\[
G_{ij}\approx \frac1M\sum_{m=1}^M\phi_i(x_m)\phi_j(x_m),
\qquad
A_{ij}\approx \frac1M\sum_{m=1}^M\phi_i(x_m)\,\nabla\phi_j(x_m)\cdot F(x_m),
\]
and the projected generator is
\[
K_N=G^{-1}A.
\]
After spectral decomposition of \(K_N\), one obtains a finite-resolution spectral measure
\[
\mu_f^{(N)}(\omega)=\sum_{j=1}^N w_j\,\delta(\omega-\lambda_j)
\]
for the observable under study [2602.23421].

To analyze spreading under the generator, the framework then constructs a Krylov chain
\[
f_0=f,\quad f_1=\mathcal Lf,\quad f_2=\mathcal L^2f,\dots
\]
and orthonormalizes it by Lanczos recursion to obtain basis vectors \(|n)\) and a tridiagonal matrix
\[
T_m=
\begin{pmatrix}
a_0&b_1&0&\dots\\
b_1&a_1&b_2&\dots\\
0&b_2&a_2&\dots\\
\vdots&\vdots&\vdots&\ddots
\end{pmatrix}.
\]
The time-evolved observable is expanded as
\[
f(\sigma)=e^{\sigma\mathcal L}f_0=\sum_{n=0}^{m-1}\phi_n(\sigma)\,|n),
\qquad
p_n(\sigma)=|\phi_n(\sigma)|^2,
\]
with \(\phi(\sigma)=\exp(\sigma T_m)e_1\) in the Lanczos basis [2602.23421].

From the Krylov probabilities \(p_n(\sigma)\), the paper defines diagnostics including Krylov complexity
\[
C_K(\sigma)=\sum_{n=0}^{m-1}n\,p_n(\sigma),
\]
the inverse participation ratio
\[
\mathrm{IPR}(\sigma)=\sum_n p_n(\sigma)^2,
\qquad
N_{\mathrm{eff}}(\sigma)=\frac1{\mathrm{IPR}(\sigma)},
\]
the Krylov entropy
\[
S_K(\sigma)=-\sum_n p_n(\sigma)\log p_n(\sigma),
\]
and a finite-resolution spectral measure from the Lanczos tridiagonalization,
\[
\mu_f^{(m)}(\omega)=\sum_{j=1}^m |(e_1,u_j)|^2\,\delta(\omega-\lambda_j).
\]
To compare integrable and deformed dynamics, the framework also uses the Wasserstein-1 distance
\[
W_1(\mu_{\mathrm{int}},\mu_{\mathrm{def}})
=
\int_{\mathbb R}|F_{\mathrm{int}}(\omega)-F_{\mathrm{def}}(\omega)|\,d\omega,
\]
and sector-resolved leakage fractions \(f_{\mathrm{prot}}\) and \(f_{\mathrm{rest}}\) derived from participation ratios in protected and complementary sectors [2602.23421].

The applications cover three classes of non-integrable semiclassical string solutions: a two-loop Landau-Lifshitz \(SU(2)\) sector, a Leigh-Strassler \(SU(3)\) Landau-Lifshitz sector, and two near-Penrose-limit ansätze of \({\rm AdS}_5\times T^{1,1}\). Across these examples, the reported pattern is that integrability breaking induces redistribution of spectral weight, observable-dependent delocalization and spreading in Krylov space, and nonzero Wasserstein distance between integrable and deformed spectral measures. The paper also emphasizes that the effect is not an abrupt ballistic growth in Krylov space, but a gradual, resonance-driven redistribution of spectral weight [2602.23421].

A plausible unifying implication of these results is that "Koopman-Krylov" is not a single algorithm but a common reduction principle: approximate the Koopman operator or generator on a subspace generated by iterates of the dynamics, then read spectral, modal, control, or transport information from the resulting finite representation. The specific realization may be a companion matrix, a Hankel/SVD reduction, a robust cross-covariance estimator, a finite-section Galerkin truncation, or a Lanczos tridiagonalization, but the shared mathematical object is the Krylov subspace generated by Koopman action [1808.09557] [2009.05883] [2011.00078] [2407.00432] [2602.23421].

Source: https://www.emergentmind.com/topics/koopman-krylov-framework