REVIEW 4 major objections 5 minor 33 references
Carleman-Fourier Linearization of Complex Dynamical Systems: Convergence and Explicit Error Bounds
T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The paper proves explicit exponential error bounds for a Fourier-exponential version of Carleman linearization of complex nonlinear systems with periodic vector fields, over a computable time horizon.
desk verdict Real extension of Carleman-Fourier linearization with explicit exponential bounds, but the advertised scope overstates what is proven for real periodic fields with Fourier decay radius R≤e. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the Fourier-state lift: new variables $w_\alpha=e^{i\alpha x}$ for nonzero nonnegative multi-indices $\alpha$, grouped in blocks by total degree $k=1,\dots,N$. Under the analyticity condition $g_\alpha(t)=0$ for all negative multi-indices, the derivative of $w_\alpha$ only involves $w_\beta$ with $|\beta|\ge|\alpha|$, so the infinite matrix $B(t)$ is block upper triangular with diagonal blocks $i\operatorname{diag}(\alpha^T g_0(t))$. The finite-section system keeps the first $N$ blocks. The explicit error bound is obtained by writing the error $u_k=v_{k,N}-w_k$ in integral form with the diagonal kernel $K_k(t,s)=\exp(i\int_s^t \alpha^T g_0\,du)$, bounding the coupling blocks by the Schur-norm estimate $\|B_{k,l}\|_S\le D_0 k R^{k-l}$, and applying a discrete Gronwall-type inequality. The augmented state $[x^T,-x^T]^T$ is a second mechanism that manufactures the analyticity condition for general multi-frequency vector fields by pairing every negative index with a positive one.
What would settle it
For the system $\dot{x}=i(1-e^{ix})$ with initial $x_0=iy$, compute the $N$-truncated first block from (5.9) and compare $\max_{t\le T}|v_{1,N}(t)e^{-ix(t)}-1|$ with the right-hand side of (3.11); a discrepancy at fixed $N$, $T$ would refute the bound. Alternatively, adding a small term $\epsilon e^{-ix}$ to the vector field should break the block-triangular structure, so the observed error should stop shrinking exponentially in $N$, confirming that condition (1.9) is the load-bearing premise.
Extended reading notes
Core claim
The central discovery is that replacing monomials $x^{\alpha}$ by Fourier exponentials $e^{i\alpha x}$ in Carleman's lifting scheme converts a complex dynamical system $\dot{x}=g(t,x)$ with a one-sided periodic vector field into an infinite-dimensional linear system with a block upper-triangular matrix $B(t)$. For such systems, the first block $v_{1,N}$ of the $N$-th finite-section approximation satisfies the explicit bound $\max_j |v_{j,N}(t)e^{-ix_j(t)}-1| \le C_0 N^{-3/2} e^{D_0 t N} (e\|\exp(ix_0)\|_\infty/R)^{(e-1)N/(2e-1)}$ for $0\le t\le T^*_{CF}$; taking logarithms and writing $v_{j,N}=e^{i\xi_{j,N}}$ yields an approximation $\xi_{j,N}$ of the original state $x_j$ whose error is at most four times the same bound. When the constant Fourier coefficient $g_0$ has strictly positive imaginary part, the convergence extends to all $t\ge 0$ with the rate $((D_0+\mu_0)\|\exp(ix_0)\|_2/(\mu_0 R))^N$. For vector fields like Kuramoto's that have negative frequencies, the paper shows that lifting the augmented state $[x^T,-x^T]^T$ restores the hypothesis and proves the analogous exponential bound under the condition $R>e$ for real initial states.
Load-bearing premise
The central proofs require the periodic vector field to have only nonnegative Fourier frequencies, condition (1.9), which makes the lifted system triangular; when it fails, the paper's remedy only applies under the additional restriction that the Fourier decay radius exceed $e$ for real initial states.
Editorial extensions
If this is right
- For systems satisfying condition (1.9) and Assumption 1.1, an order-$N$ truncation of the lifted system approximates $e^{ix(t)}$ with the explicit error bound of Theorem 3.1 on $[0,T^*_{CF}]$, so a desired tolerance directly fixes the truncation order.
- Under the positive-imaginary condition (1.14) and the small-initial-state condition (1.15), the same finite-section approximation converges exponentially for all $t\ge 0$ with rate $((D_0+\mu_0)\|\exp(ix_0)\|_2/(\mu_0 R))^N$.
- For vector fields with multiple fundamental frequencies that fail (1.9), the augmented state $[x^T,-x^T]^T$ restores the analyticity condition; with real initial states and $R>e$, exponential convergence holds on $[0,\widetilde{T}^*_{CF}]$ (Corollary 4.2).
- When $1\le\|x_0\|_\infty<e^{-1}\ln R$, the Fourier method's guaranteed time horizon is at least as long and its convergence rate no worse than monomial Carleman linearization, as stated in (1.19).
- The error bound depends on the initial state only through its imaginary parts, so the guaranteed accuracy is insensitive to how large the real parts of $x_0$ are.
Reading between the lines
- The formulas suggest a quantitative notion of a usable basin in complexified phase space: shifting a real initial state upward in the imaginary direction should extend the guaranteed horizon exactly as (1.11) predicts, a prediction one could test numerically on the example $\dot{x}=a(1-e^{ix})$.
- The augmented-state doubling in Section 4 doubles the lifted dimension; for systems with additional symmetry, such as the zero-sum phases of the normalized Kuramoto model, a smaller symmetry-adapted basis might satisfy the analyticity condition with less overhead.
- For quantum simulation of dissipative polynomial dynamics, the explicit $N$-dependence in Theorems 3.1 and 3.3 could be converted into a query or qubit count, though the paper only lists quantum computing as a motivation.
- The comparison with real-valued systems suggests the real case has strictly better guaranteed horizons; extrapolating, a complex system whose vector field is real on the real axis should be approximated more efficiently by first separating real and imaginary parts rather than by the direct complex lift.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces Carleman-Fourier linearization for complex nonlinear dynamical systems with periodic vector fields. The method lifts the system using Fourier monomials e^{iα·x}; under the analyticity condition (1.9) the lifted matrix is block upper-triangular. Theorems 3.1 and 3.3 state explicit exponential error bounds for finite-section approximations on a finite horizon and on the entire time axis, respectively. Section 4 handles multiple fundamental frequencies by augmenting the state vector with its negative, giving Theorem 4.1 and Theorem 4.3. Section 5 presents numerical experiments for the scalar equation ẋ=a(1−be^{ix}) and for the Kuramoto model.
Significance. If the results are correct, the paper gives a useful extension of Carleman linearization to periodic vector fields, with explicit and computable error bounds, finite-section truncation criteria, and a natural treatment of multiple frequencies through state augmentation. The block upper-triangular observation and the global convergence result under positivity are interesting, and the numerical experiments support the qualitative claims. The main weaknesses are the heavy reliance on arguments and lemmas from the authors' earlier work [2], the omission of the proof of Theorem 2.1, a concrete incorrect constant bound in Section 4.2, and a scope limitation for real initial states when the Fourier decay radius R≤e that is not reflected in the abstract.
major comments (4)
- [Section 4.2, Eq. (4.16)] The asserted inequality C1 ≤ R^2/(2π e(e−1)) is false for the C1 defined in Theorem 4.1. Substituting the condition max|ℑ(ω_l x_{0,j})| < ln R − 1 into the displayed formula for C1 yields only C1 ≤ R^2/(√(2π) e(e−1)), and the stronger claimed bound fails, e.g., with R=4 and max|ℑ|=0.1 one obtains C1≈1.22 while R^2/(2π e(e−1))≈0.545. Since Corollary 4.2 and Eq. (4.19) use this simplification, the constants in those statements need to be corrected or the stronger inequality must be proved by a different argument.
- [Section 2, Theorem 2.1] Theorem 2.1 is stated with an explicit error bound but its proof is omitted: the text says the argument in [2] can be followed and the details are omitted. This theorem is listed as a contribution and is used in the comparison (1.19). The paper should either provide the proof in Section 6 or state precisely which theorem of [2] implies the complex-case bound with the same constants.
- [Section 6, Lemmas 6.2 and 6.3] Lemmas 6.2 and 6.3 are quoted from [2] without proof and are load-bearing for the proof of Theorem 3.1. Since the current paper advertises self-contained explicit error bounds, these lemmas should either be proved in the appendix or stated as cited theorems with exact references, so that a reader can verify the constants used in the main derivation.
- [Abstract and Section 4, Corollary 4.2(ii)] The abstract and introduction claim applicability to 'periodic vector fields' without qualification, but for real initial states the proved exponential convergence requires R>e (Corollary 4.2(ii)); for a real periodic vector field whose optimal Fourier decay radius satisfies R≤e, no exponential convergence of the finite-section approximation is established by the theorems in this paper. This is a genuine scope limitation, acknowledged only later in Section 4, and the claims in the abstract and in the comparison (1.19) should be qualified accordingly.
minor comments (5)
- [Theorem 3.1, Eq. (3.11)] The phrase 'exponential convergence' should be qualified: at t=T*_CF the factor e^{D0tN}(e∥exp(ix0)∥∞/R)^{(e−1)N/(2e−1)} equals 1, so the bound degenerates to O(N^{-3/2}) at the endpoint. Exponential-in-N convergence is established for t<T*_CF, as used in Corollary 3.2.
- [Eq. (3.11)] There is a typo in the displayed inequality: 'and and C0' contains a duplicated word.
- [Section 4.2 and Eq. (4.19)] The constant correction in the first major comment should be propagated to Corollary 4.2 and to the bound in Eq. (4.19), which currently inherit the incorrect simplification from Eq. (4.16).
- [Section 6.3] The proof of the comparison (1.19) is very sketchy and relies on numerical constants such as 4.9215 and 0.7076 without a complete derivation; a fully justified argument should be supplied.
- [References] Reference [26] is cited as 'In preparation' and is used for the comparison with real-system results; the paper should cite a published or otherwise publicly verifiable version, or the comparison should be proved directly.
Circularity Check
No significant circularity: the Carleman-Fourier error bounds are explicit functions of the stated assumptions and are not fitted or defined in terms of the predicted quantities.
full rationale
The paper's derivation chain is not circular. The central error bound (3.11) in Theorem 3.1 is an explicit function of the assumed Fourier-decay constants D0 and R, the initial condition x0, the truncation order N, and the time t; no parameter is fitted to the error that the bound is supposed to predict. The analyticity condition (1.9) is a structural assumption that forces the lifted matrix B(t) to be block upper triangular (3.5), and the proof independently controls w1(t)=e^{ix(t)} via Lemma 6.1 and then derives the finite-section error through Lemmas 6.2-6.3. Although Lemmas 6.2, 6.3 and the proof of Theorem 2.1 are cited from the authors' own [2], those are general ODE and combinatorial tools whose stated assumptions do not include the target result of this paper; under Rule 4 they count as independent support rather than circularity. Section 4's augmented-state construction removes (1.9), and the real-initial-state restriction R>e in Corollary 4.2(ii) is a genuine, explicitly acknowledged scope limitation, not a circular move: the bound is still proved from (4.2), (4.6), and (4.7) rather than assumed. The comparison with the unpublished [26] is contextual and not load-bearing for any theorem in this paper. There is no self-definition, no fitted input renamed as a prediction, and no uniqueness claim imported from the authors' prior work to force the choice of method. Thus the analysis is self-contained against the stated assumptions, and no circularity score above zero is warranted.
Assumptions & free parameters
assumptions (5)
- domain assumption Assumption 1.1: sup_t sum_{|alpha|=k} sum_j |g_{j,alpha}(t)| <= D0 R^{-k} for all k>=0
- domain assumption Analyticity condition (1.9): g_alpha(t)=0 for all alpha in Z^d \ Z^d_+ and t>=0
- domain assumption Positivity condition (1.14): min_{1<=j<=d} Im g_{j,0}(t) >= mu0 > 0 for all t>=0
- standard math Lemmas 6.2 and 6.3 (taken from [2, Lemmas 5.3 and 5.4]) giving integral representations and combinatorial bounds for the finite-section error system
- standard math The standard inequality |z mod 2pi| <= 4epsilon for all z in C with |e^{iz}-1| <= epsilon <= 1/2
Cite this review
Pith. "Pith review of Carleman-Fourier Linearization of Complex Dynamical Systems: Convergence and Explicit Error Bounds." pith.science (2026). https://pith.science/paper/FFFKNECY
@misc{pith2026241111598,
author = {Pith},
title = {Pith review of: Carleman-Fourier Linearization of Complex Dynamical Systems: Convergence and Explicit Error Bounds},
year = {2026},
howpublished = {\url{https://pith.science/paper/FFFKNECY}},
note = {Machine review of arXiv:2411.11598}
}
read the original abstract
This paper presents a Carleman-Fourier linearization method for nonlinear dynamical systems with periodic vector fields involving multiple fundamental frequencies. By employing Fourier basis functions, the nonlinear dynamical system is transformed into a linear model on an infinite-dimensional space. The proposed approach yields accurate approximations over extended regions around equilibria and for longer time horizons, compared to traditional Carleman linearization with monomials. Additionally, we develop a finite-section approximation for the resulting infinite-dimensional system and provide explicit error bounds that demonstrate exponential convergence to the original system's solution as the truncation length increases. For specific classes of dynamical systems, exponential convergence is achieved across the entire time horizon. The practical significance of these results lies in guiding the selection of suitable truncation lengths for applications such as model predictive control, safety verification through reachability analysis, and efficient quantum computing algorithms. The theoretical findings are validated through illustrative simulations.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
- [2]
-
[26]
N. Motee and Q. Sun, On exponential convergence of Carleman-Fourier linearization of nonlinear real dynamical systems, In prepartion
- [1]
- [3]
- [4]
- [5]
- [6]
-
[7]
Brockett, The early days of geometric nonlinear control, Automatica, 50(9), 2014, pp
R. Brockett, The early days of geometric nonlinear control, Automatica, 50(9), 2014, pp. 2203–2224
work page 2014
Show all 33 references
-
[8]
J. C. Bronski, T. E. Carty and L. DeVille, Synchronisation conditions in the Kuramoto model and their relationship to seminorms, Nonlinearity, 34(8), 2021, pp. 5399–5433
2021
-
[9]
S. L. Brunton, J. L. Proctor and J. N. Kutz, Discovering governing equations from data by sparse identification of nonlinear dynamical systems, Proceedings of the National Academy of Sciences, 113(15), 2016, pp. 3932–3937
2016
-
[10]
Dietert, Stability and bifurcation for the Kuramoto model, Journal de Math´ ematiques Pures et Appliqu´ ees, 105(4), 2016, pp
H. Dietert, Stability and bifurcation for the Kuramoto model, Journal de Math´ ematiques Pures et Appliqu´ ees, 105(4), 2016, pp. 451–489
2016
-
[11]
Forets and A
M. Forets and A. Pouly, Explicit error bounds for Carleman linearization, 2017, arXiv preprint arXiv:1711.02552
2017 arXiv
-
[12]
Forets and C
M. Forets and C. Schilling, Reachability of weakly nonlinear systems using Carleman linearization, International Conference on Reachability Problems , Springer, 2021, pp. 85–99
2021
-
[13]
Y. Guo, D. Zhang, Z. Li, Q. Wang and D. Yu, Overviews on the applications of the Kuramoto model in modern power system analysis, International Journal of Electrical Power & Energy Systems, 129, 2021, article no. 106804
2021
-
[14]
Hashemian and A
N. Hashemian and A. Armaou, Fast moving horizon estimation of nonlinear processes via Carleman linearization, 2015 American Control Conference (ACC) , IEEE, 2015, pp. 3379–3385
2015
-
[15]
O. A. Heggli, J. Cabral, I. Konvalinka, P. Vuust and M. L. Kringelbach, A Kuramoto model of self-other integration across interpersonal synchronization strategies, PLoS Computational Biology, 15(10), 2019, article no. e1007422, 17pp
2019
-
[16]
P. Ji, T. K. Peron, F. A. Rodrigues and J. Kurths, Low-dimensional behavior of Kuramoto model with inertia in complex networks, Scientific Reports, 4(1), 2014, article no. 4783
2014
-
[17]
Korda and I
M. Korda and I. Mezi´ c, Linear predictors for nonlinear dynamical systems: Koopman operator meets model predictive control, Automatica, 93, 2018, pp. 149–160
2018
-
[18]
Korda and I
M. Korda and I. Mezi´ c, On convergence of extended dynamic mode decomposition to the Koopman operator, Journal of Nonlinear Science , 28, 2018, pp. 687–710
2018
-
[19]
Kowalski and W.H
K. Kowalski and W.H. Steeb, Nonlinear dynamical systems and Carleman linearization, World Scientific, 1991
1991
-
[20]
A. J. Krener, Linearization and bilinearization of control systems, Proceedings of the 1974 Allerton Conference on Circuit and Systems Theory, Urbana III , 1974
1974
-
[21]
A. J. Krener, Bilinear and nonlinear realizations of input-output maps, SIAM Journal on Control, 13(4), 1975, pp. 827–834
1975
-
[22]
Kuramoto, Chemical oscillations, waves, and turbulence, New York, Springer-Verlag, 1984
Y. Kuramoto, Chemical oscillations, waves, and turbulence, New York, Springer-Verlag, 1984
1984
-
[23]
J.P. Liu, H. O. Kolden, H. K. Krovi, N. F. Loureiro, K. Trivisa and A. W. Childs, Efficient quantum algorithm for dissipative nonlinear differential equations, Proceedings of the National Academy of Sciences, 118(35), 2021, article no. e2026805118
2021
-
[24]
Loparo and G
K. Loparo and G. Blankenship, Estimating the domain of attraction of nonlinear feedback systems, IEEE Transactions on Automatic Control , 23(4), 1978, pp. 602–608. CARLEMAN-FOURIER LINEARIZATION 43
1978
-
[25]
Minisini, A
J. Minisini, A. Rauh and E. P. Hofer, Carleman linearization for approximate solutions of nonlinear control problems: Part 1–theory, Proc. of the 14th Intl. Workshop on Dynamics and Control, 2007, pp. 215–222
2007
-
[27]
Pruekprasert, J
S. Pruekprasert, J. Dubut, T. Takisaka, C. Eberhart and A. Cetinkaya, Moment propagation through Carleman linearization with application to probabilistic safety analysis, 2022, arXiv preprint arXiv:2201.08648
2022 arXiv
-
[28]
A. Rauh, J. Minisini and H. Aschemann, Carleman linearization for control and for state and disturbance estimation of nonlinear dynamical processes, IF AC Proceedings Volumes, 42(13), 2009, pp. 455–460
2009
-
[29]
Rotondo, G
D. Rotondo, G. Luta and J. H. U. Aarvag, Towards a Taylor-Carleman bilinearization approach for the design of nonlinear state-feedback controllers, European Journal of Control 68 , 2022, article no. 100670
2022
-
[30]
W. H. Steeb and F. Wilhelm, Non-linear autonomous systems of differential equations and Carleman linearization procedure, Journal of Mathematical Analysis and Applications , 77(2), 1980, pp. 601–611
1980
-
[31]
Surana, A
A. Surana, A. Gnanasekaran and T. Sahai, An efficient quantum algorithm for simulating polynomial dynamical systems, Quantum Information Processing, 23(3), 2024, article no. 105, 22pp
2024
-
[32]
Z. Wang, R. M. Jungers and C. J. Ong, Computation of invariant sets via immersion for discrete-time nonlinear systems, Automatica, 147, 2023, article no. 110686, 9pp
2023
-
[33]
H. C. Wu, J. Wang and X. Li, Quantum algorithms for nonlinear dynamics: revisiting Carleman linearization with no dissipative conditions, 2024, arXiv preprint arXiv:2405.12714 . Department of Mathematics, University of Central Florida, Orlando, Florida 32816 Email address : pa...
2024 arXiv
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.