REVIEW 3 major objections 5 minor 39 references
Koopman operator based identification of nonlinear networks
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read The paper claims that a two-step Koopman-operator method recovers both the interaction graph and the local dynamics of heterogeneous nonlinear networks from input-state snapshots, with data cost set by sparsity rather than dictionary size.
desk verdict A solid two-step Koopman network identification method with strong numerics, but the edge-detection guarantee is void for the experiments shown. 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 dual Koopman method: approximate the Koopman operator and its generator on the finite space of samples rather than on a chosen dictionary, use Gaussian radial-basis test functions selected to minimize the vector-field approximation error, and read off the vector-field values at the sample points. Edge detection then uses per-node subspaces $F_{x,k}=\mathrm{span}\{\phi_{kl}\}$ and $F_{u,k}=\mathrm{span}\{\psi_{kl}\}$; the zero-mean conditions (25)-(26) make these subspaces mutually orthogonal, so the projection weight $\Lambda_{ik}$ of node $i$'s vector field onto node $k$'s subspace is nonzero exactly when node $k$ influences node $i$ (Theorem 1). Once neighborhoods are known, the local identification step fits lifted matrices $A_i, E_i, B_i$ and converts them back to the unknown coefficients through the continuous-time formulas (43)-(44) and the extraction rule (47).
What would settle it
Take a two-node network with no true edge and node functions {$x_k, x_k^2$} on $[-1,1]$; since $x_k^2$ has nonzero mean, the orthogonality behind Theorem 1 fails, so a noise-free run that still yields $\Lambda_{ik}$ above threshold for the non-neighbor would show that the zero-mean condition is doing the theoretical work and must be enforced in practice.
Extended reading notes
Core claim
The central claim is that the Boolean interaction graph and the full local dynamics can both be recovered from input-state snapshots without solving a regression over the entire network dictionary. Theorem 1 guarantees that weights $\Lambda_{ik}$ and $\Delta_{ik}$, computed from orthogonal projections of each node's vector-field component onto the spans of node functions, vanish exactly when the corresponding coupling or input term is absent, provided the node functions integrate to zero and the true terms are not orthogonal to those spans. Under that guarantee, a sparse regression on the estimated vector-field samples yields per-node neighbor sets; a second, local Koopman lift then estimates the parameters of $F_i$, $H_{ik}$, and $G_{ik}$. The paper also produces a modular linear lifted model of the whole network, assembled blockwise from the local fits.
Load-bearing premise
The edge-detection guarantee relies on the node functions integrating to zero over the state and input domains so that different nodes' subspaces are orthogonal, and the parameter estimates rely on the chosen local dictionaries being nearly invariant under the Koopman generator; the numerical examples use monomials such as $x_i^2$ and $u_i^2$ that do not integrate to zero, so the no-false-positive guarantee is not operative in those demonstrations.
Editorial extensions
If this is right
- Boolean network reconstruction and local dynamics identification can be performed from the same dataset, so no separate topological probing experiment is needed.
- Under network sparsity, the number of samples needed scales with the number of neighbors actually present rather than with the total dictionary size, making large sparse networks tractable.
- Heterogeneous nodes are handled explicitly: each node gets its own local dictionary, so internal dynamics, coupling functions, and input effects may all differ across nodes.
- A by-product linear lifted model of the whole network, assembled blockwise from local fits, can be used for analysis and control design.
- Sampling time $T_s$ directly controls accuracy because the continuous-time parameter formulas come from a zero-order-hold approximation of the lifted dynamics.
Reading between the lines
- A testable extension is to center and orthogonalize node functions empirically against the sampling distribution before computing projections, which could preserve the no-false-positive property for polynomial dictionaries like the monomials used in the numerical demonstrations.
- Because Theorem 1's converse also requires the true coupling to be non-orthogonal to the node-function subspace, a dictionary that misses a mode of the coupling will misclassify that edge; a diagnostic would be to inspect the residual norm after each projection.
- The modular lifted model opens a route to distributed prediction and control, since each node's block is constructed independently from local neighbors once the graph is known.
- A natural stress test is to run the two-step method on partially observed networks, where hidden nodes would likely be absorbed into the recovered local dynamics and could create spurious edges; the method's guarantees appear to require full-state measurement.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a two-step identification method for continuous-time nonlinear networked dynamics with heterogeneous local dynamics, external inputs, and node-dependent coupling functions. In the first step, the authors extend the dual Koopman method to non-autonomous systems and use it to estimate the vector field, then regress it on 'node functions' to determine each node's neighbor set and input set (Boolean topology reconstruction). In the second step, they build a local finite-dimensional Koopman lifted dynamics for each node from the estimated neighbors and recover the coefficients of local and coupling dictionary functions. The main theoretical result, Theorem 1, states that under zero-mean conditions on the node functions, the projection weights Lambda_ik are nonzero exactly when there is an edge k to i. Numerical experiments on Erdos-Renyi graphs, non-polynomial vector fields, and Hindmarsh-Rose neuronal networks support scalability (up to 1500 nodes) and accuracy, and the paper also constructs a modular lifted representation of the network dynamics.
Significance. If the stated guarantees held for the demonstrated configurations, this would be a valuable contribution to nonlinear network identification: it relaxes the homogeneity assumptions of earlier Koopman-based network identification, provides a data-efficiency argument for sparse networks, and yields a modular lifted model that can be reused for control or analysis. The paper includes a broad set of numerical experiments and comparisons with the dual method of [18]. However, the gap between the assumptions of Theorem 1 and the node functions used in Sections 5.1-5.3, together with the unstated invariance assumption in the local identification step, means that the current manuscript does not yet rigorously support the central claim as demonstrated. The core idea is promising and the necessary fixes are identifiable, so the issues are reparable within the scope of the paper.
major comments (3)
- [Theorem 1, Section 4.2.2 and Sections 5.1-5.2] The edge-detection guarantee is void for the node functions used in the experiments. Equations (25)-(26) require all node functions to have zero mean, and the proof of Theorem 1 uses this to conclude that F_x,k is orthogonal to F_x,l and to F_u,l. However, in all numerical sections the node functions are {x_i, x_i^2} and {u_i, u_i^2} on [-1,1], and x_i^2 and u_i^2 have mean 2/3, not zero. For a two-node example with no edge 2 to 1 and F_1 = x_1^2, the L2 projection onto span{x_2, x_2^2} has coefficient 5/9 on x_2^2, so the analogue of Lambda_12 is nonzero. Thus the no-false-positive statement does not cover the demonstrated results; the reported AUROC values are obtained through thresholding and Lasso rather than through Theorem 1. The authors should either use zero-mean node functions (for example, shifted or Legendre polynomials) in the experiments, or extend Theorem 1 to general node functions with explicit correction or threshold bounds.
- [Section 4.2.2, Eqs. (27)-(33)] The theorem assumes that Lambda_ik is computed from the exact orthogonal projection in (27), but Algorithm 1 computes the coefficients from the Lasso problem (33), which is a different estimator. Even when (25)-(26) hold, a regularized regression does not coincide with the projection, so Theorem 1 does not by itself justify the statement that the estimated Lambda_ik are the projection weights used in the proof. The relation between the Lasso solution and the projection coefficients needs to be quantified, or Theorem 1 needs to be restated for the finite-sample regularized estimator, including the effect of the threshold delta in (34).
- [Section 4.3.2, Eqs. (40)-(44) and (47)] The local parameter recovery assumes that the finite-dimensional dictionary span F_i is invariant, or nearly invariant, under the Koopman operator or its generator, so that the discrete lifted matrices in (40) describe an exact finite-dimensional representation and log(Abar_i) in (43) equals the generator matrix restricted to that span. This invariance assumption is not stated. For general dictionaries such as the monomials, sinusoids, and exponentials used in Section 5.2, the span is not invariant, and then the coefficient estimates in (47) are not guaranteed to converge to the true alpha, beta, and gamma. A bound on the effect of the projection error on (47), or a precise approximation assumption, is needed to make the central identification claim rigorous.
minor comments (5)
- [Section 3.2, Proposition 1] The proof drops O(T_s^2) terms before taking the limit; a uniform remainder bound over the finite set of test-function choices would be needed to justify the interchange of the limit and the argmin.
- [Sections 2.1 and 4.1] The shift operator T^t is called a left shift in Section 2.1 and a right shift in Section 4.1, although both definitions are the same; please make the terminology consistent.
- [Section 4.3.1] The symbol Mbar_i is used both for the estimated input set and for the total number of input dictionary functions; please use distinct notations for these two objects.
- [Throughout Section 5] The acronym 'RSME' appears in figure captions and Table 1; it should be 'RMSE' for consistency with the definition in Section 5.
- [Section 5.2, Eq. (50)] The expressions 'sin(xti)' and 'exti' are typeset in a way that makes it difficult to distinguish the variable x_{t_i} from the exponential e^{x_{t_i}}; please clarify the notation and define all variables in the piecewise dynamics.
Circularity Check
No significant circularity: the topology and local identification results are derived from the dynamics and orthogonality assumptions, not from fitted parameters or self-citation chains.
full rationale
The paper's two-step method is not circular. Theorem 1 derives edge detection from the zero-mean conditions (25)-(26) and the additive structure (21) via orthogonal projection; the conclusion does not presuppose the neighbor sets. The local identification in Section 4.3 estimates parameters by solving a least-squares problem for the lifted dynamics (38)-(47); the estimated parameters are not fed back into the theorem that justifies them. The authors rely on the dual and main methods of [18], but that is a published external method, and the present results are extensions rather than restatements of it. Proposition 1 selects test functions by minimizing a training error, but it is an asymptotic equivalence result, not a prediction claimed against held-out data, and it is not used to produce the headline network identification claims. The numerical experiments use node functions {x_i, x_i^2} and {u_i,u_i^2}, which do not satisfy the zero-mean assumptions of Theorem 1; this is a mismatch between theory and implementation that affects the force of the stated guarantee, but it is a correctness and robustness issue, not a circularity. No central claim reduces by construction to its inputs, so the circularity score is 0.
Assumptions & free parameters
free parameters (4)
- threshold δ =
varies in experiments
- RBF width γ =
0.001, 0.0025 to 0.025 in examples
- Lasso regularization ρ =
not reported
- node function sets {φ_kl, ψ_kl} =
monomials up to degree 2 in examples
assumptions (5)
- domain assumption The true dynamics are exactly representable as linear combinations of known dictionary functions (24).
- domain assumption Node functions satisfy zero-mean conditions (25)-(26).
- domain assumption The finite-dimensional EDMD and dual method yield an accurate approximation of the Koopman operator and generator on the data.
- domain assumption The matrix logarithm (1/Ts)log(A) is well-defined for the approximated operators.
- domain assumption All nodes are measured and inputs are held constant over each sampling interval.
Cite this review
Pith. "Pith review of Koopman operator based identification of nonlinear networks." pith.science (2026). https://pith.science/paper/LDQUXMO7
@misc{pith2026241218492,
author = {Pith},
title = {Pith review of: Koopman operator based identification of nonlinear networks},
year = {2026},
howpublished = {\url{https://pith.science/paper/LDQUXMO7}},
note = {Machine review of arXiv:2412.18492}
}
read the original abstract
In this work, we develop a method to identify continuous-time nonlinear networked dynamics via the Koopman operator framework. The proposed technique consists of two steps: the first step identifies the neighbors of each node, and the second step identifies the local dynamics at each node from a predefined set of dictionary functions. The technique can be used to either identify the Boolean network of interactions (first step) or to solve the complete network identification problem that amounts to estimating the local node dynamics and the nature of the node interactions (first and second steps). Under a sparsity assumption, the data required to identify the complete network dynamics is significantly less than the total number of dictionary functions describing the dynamics. This makes the proposed approach attractive for identifying large dimensional networks with sparse interconnections. The accuracy and performance of the proposed identification technique are demonstrated with several examples.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
-
[18]
Mauroy, A., Goncalvez, J.: Koopman-based lifting techniques for nonlinear systems identification. IEEE Trans. on Automatic Control 65(6), 2550–65 (2020)
work page 2020
-
[1]
Neuroscientist 23(5), 499–516 (October 2017)
Bassett, D.S., Bullmore, E.T.: Small-world brain networks revisited. Neuroscientist 23(5), 499–516 (October 2017)
work page 2017
-
[2]
ACM Computing Surveys (CSUR) 51(2), 1– 39 (2018)
Brugere, I., Gallagher, B., Berger-Wolf, T.Y.: Network structure inference, a survey: Motivations, methods, and applications. ACM Computing Surveys (CSUR) 51(2), 1– 39 (2018)
work page 2018
-
[3]
Automatica 48(8), 1553–1565 (2012) 28 Ramachandran Anantharaman, Alexandre Mauroy Fig
Chiuso, A., Pillonetto, G.: A bayesian approach to sparse dynamic network identifica- tion. Automatica 48(8), 1553–1565 (2012) 28 Ramachandran Anantharaman, Alexandre Mauroy Fig. 9: Predicted vs true values of parameters at node 50 Fig. 10: Local identification of a Watts-Strogatz network with Hindmarsh-Rose dynamics
work page 2012
-
[4]
IEEE Transactions on Automatic Control 61(4), 937–952 (2015)
Dankers, A., Van den Hof, P.M., Bombois, X., Heuberger, P.S.: Identification of dynamic models in complex networks with prediction error methods: Predictor input selection. IEEE Transactions on Automatic Control 61(4), 937–952 (2015)
work page 2015
-
[5]
Automatica 49(4), 1031–1036 (2013)
Franceschelli, M., Gasparri, A., Giua, A., Seatzu, C.: Decentralized estimation of Lapla- cian eigenvalues in multi-agent systems. Automatica 49(4), 1031–1036 (2013)
work page 2013
-
[6]
IEEE Transactions on Automatic Control 53, 1670– 1674 (2008)
Goncalves, J., Warnick, S.: Necessary and sufficient conditions for dynamical structure reconstruction of LTI networks. IEEE Transactions on Automatic Control 53, 1670– 1674 (2008)
work page 2008
-
[7]
Modularized data-driven approximation of the Koopman operator and generator
Guo, Y., Schaller, M., Worthmann, K., Streif, S.: Modularized data-driven approxima- tion of the Koopman operator and generator. arXiv preprint arXiv:2408.12277 (2024)
work page Pith review arXiv 2024
Show all 39 references
-
[8]
Physica A: Statistical Mechanics and Its Applications 392(4), 1038–1049 (2013)
He, T., Lu, X., Wu, X., Lu, J., Zheng, W.X.: Optimization-based structure identification of dynamical networks. Physica A: Statistical Mechanics and Its Applications 392(4), 1038–1049 (2013)
2013
-
[9]
IEEE Transactions on Automatic Control 64 (6) , 2240– 2253 (2018)
Hendrickx, J., Gevers, M., Bazanella, A.: Identifiability of dynamical networks with partial node measurements. IEEE Transactions on Automatic Control 64 (6) , 2240– 2253 (2018)
2018
-
[10]
IEEE Transactions on Aerospace and Electronic Systems 57(4), 2558 – 2564 (2021)
Huang, H., Savkin, A.V., Huang, C.: Decentralized autonomous navigation of a UA V network for road traffic monitoring. IEEE Transactions on Aerospace and Electronic Systems 57(4), 2558 – 2564 (2021)
2021
-
[11]
Izhikevich, E.M.: Which model to use for cortical spiking neurons? IEEE Transactions on Neural Networks 15(5), 1063–1070 (2004)
2004
-
[12]
IET systems biology 3(3), 155–166 (2009)
Julius, A., Zavlanos, M., Boyd, S., Pappas, G.J.: Genetic network identification using convex programming. IET systems biology 3(3), 155–166 (2009)
2009
-
[13]
Koopman, B.: Hamiltonian systems and transformation in hilbert space. Proc. of the National Acad. of Sciences of the USA 17(5), 315–318 (1931)
1931
-
[14]
Automatica 93, 149–60 (2018)
Korda, M., Mezi´ c, I.: Linear predictors for nonlinear dynamical systems: Koopman operator meets model predictive control. Automatica 93, 149–60 (2018)
2018
-
[15]
Journal of The- oretical Biology 228, 195–215 (2004) Koopman operator based identification of nonlinear networks 29
Lipniackia, T., et al: Mathematical model of nf- κb regulatory module. Journal of The- oretical Biology 228, 195–215 (2004) Koopman operator based identification of nonlinear networks 29
2004
-
[16]
IEEE Transactions on Automatic Control 55(8), 1860–1871 (2010)
Materassi, D., Innocenti, G.: Topological identification in networks of dynamical sys- tems. IEEE Transactions on Automatic Control 55(8), 1860–1871 (2010). DOI 10.1109/TAC.2010.2042347
2010
-
[17]
In: 2015 54th IEEE Conference on Decision and Control (CDC), pp
Materassi, D., Salapaka, M.V.: Identification of network components in presence of unobserved nodes. In: 2015 54th IEEE Conference on Decision and Control (CDC), pp. 1563–1568. IEEE (2015)
2015
-
[19]
SIAM Journal on Applied Dynamical Systems 16(1), 479–513 (2017)
Mauroy, A., Hendrickx, J.: Spectral identification of networks using sparse measure- ments. SIAM Journal on Applied Dynamical Systems 16(1), 479–513 (2017)
2017
-
[20]
Springer (2020)
Mauroy, A., Mezi´ c, I., Susuki, Y.: The Koopman Operator in Systems and Control: Concepts, Methodologies, and Applications. Springer (2020)
2020
-
[21]
Nonlinear Theory and Its Applications, IEICE 13(2) (2022)
Mei, Z., Oguchi, T.: Network structure identification via Koopman analysis and sparse identification. Nonlinear Theory and Its Applications, IEICE 13(2) (2022)
2022
-
[22]
Nonlinear Dynamics 41(1-3), 309–325 (2005)
Mezi´ c, I.: Spectral properties of dynamical systems, model reduction and decomposi- tions. Nonlinear Dynamics 41(1-3), 309–325 (2005)
2005
-
[23]
SIAM Journal on Control and Optimization 48 (1) (2009)
Michiels, W., Mor˘ arescu, C.I., Niculescu, S.I.: Consensus problems with distributed de- lays, with application to traffic flow models. SIAM Journal on Control and Optimization 48 (1) (2009)
2009
-
[24]
International Journal of Control 49(1), 219–232 (1989)
Moonen, M., De Moor, B., Vandenberghe, L., Vandewalle, J.: On-and off-line identi- fication of linear state-space models. International Journal of Control 49(1), 219–232 (1989)
1989
-
[25]
Physical Review E 77(2), 026103 (2008)
Napoletani, D., Sauer, T.D.: Reconstructing the topology of sparsely connected dynam- ical networks. Physical Review E 77(2), 026103 (2008)
2008
-
[26]
The European Physical Journal Special Topics 223, 2423–2437 (2014)
Nardelli, P.H.J., Rubido, N., Wang, C., Baptista, M.S., Pomalaza-Raez, C., Cardieri, P., Latva-aho, M.: Models for the modern power grid. The European Physical Journal Special Topics 223, 2423–2437 (2014)
2014
-
[27]
IEEE Global Conference on Signal and Information Processing p
Par´ e, P.E., Chetty, V., Warnick, S.: On the necessity of full-state measurement for state-space network reconstruction. IEEE Global Conference on Signal and Information Processing p. 803–806 (2013)
2013
-
[28]
SIAM Journal on Applied Dynamical Systems 15(1), 142–161 (2016)
Proctor, J.L., Brunton, S.L., Kutz, J.N.: Dynamic mode decomposition with control. SIAM Journal on Applied Dynamical Systems 15(1), 142–161 (2016)
2016
-
[29]
IEEE Transactions on Automatic Control 53(4), 894 – 911 (2008)
Qu, Z., Wang, J., Hull, R.A.: Cooperative control of dynamical systems with application to autonomous vehicles. IEEE Transactions on Automatic Control 53(4), 894 – 911 (2008)
2008
-
[30]
Chaos: An Interdisciplinary Journal of Nonlinear Science 28(7) (2018)
Runge, J.: Causal network reconstruction from time series: From theoretical assumptions to practical estimation. Chaos: An Interdisciplinary Journal of Nonlinear Science 28(7) (2018)
2018
-
[31]
Jour- nal of Fluid Mechanics 656, 5–28 (2010)
Schmid, P.J.: Dynamic mode decomposition of numerical and experimental data. Jour- nal of Fluid Mechanics 656, 5–28 (2010)
2010
-
[32]
Nature 440(7087), 1007–1012 (2006)
Schneidman, E., Berry, M.J., Segev, R., Bialek, W.: Weak pairwise correlations imply strongly correlated network states in a neural population. Nature 440(7087), 1007–1012 (2006)
2006
-
[33]
In: Actuators, vol
Tellez-Castro, D., Garcia-Tenorio, C., Mojica-Nava, E., Sofrony, J., Vande Wouwer, A.: Data-driven predictive control of interconnected systems using the koopman operator. In: Actuators, vol. 11(6), p. 151. MDPI (2022)
2022
-
[34]
Journal of Physics A: Mathematical and Theoretical 47(34), 343001 (2014)
Timme, M., Casadiego, J.: Revealing networks from dynamics: an introduction. Journal of Physics A: Mathematical and Theoretical 47(34), 343001 (2014)
2014
-
[35]
Journal of Theoretical Biology p
Tyson, J.J., Laomettachit, T., Kraikivski, P.: Modeling the dynamic behavior of bio- chemical regulatory networks. Journal of Theoretical Biology p. 514–527 (2019)
2019
-
[36]
arXiv preprint arXiv:2405.07636 (2024)
Vizuete, R., Hendrickx, J.M.: Nonlinear network identifiability with full excitations. arXiv preprint arXiv:2405.07636 (2024)
2024
-
[37]
EPL (Europhysics Letters) 94(4), 48006 (2011)
Wang, W.X., Yang, R., Lai, Y.C., Kovanis, V., Harrison, M.A.F.: Time-series–based prediction of complex oscillator networks via compressive sensing. EPL (Europhysics Letters) 94(4), 48006 (2011)
2011
-
[38]
Automatica 89, 247–258 (2018)
Weerts, H.H.M., den Hof, P.M.J.V., Dankers, A.G.: Identifiability of linear dynamic networks. Automatica 89, 247–258 (2018)
2018
-
[39]
Physical Review Letters 97(18), 188701 (2006)
Yu, D., Righero, M., Kocarev, L.: Estimating topology of networks. Physical Review Letters 97(18), 188701 (2006)
2006
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.