REVIEW 4 major objections 5 minor 26 references
Physics-informed Learning for Passivity-based Tracking Control
T0 review · 4 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read This paper claims that a Gaussian-process model of a partially unknown port-Hamiltonian system, combined with a modified matching equation and a damping condition, renders the tracking-error equilibrium stable with probability at least…
desk verdict A credible but incremental tracking extension of GP-PHS; the stability theorem is conditional on a strong dissipation inequality that the simulation never verifies. 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 central object is the Gaussian process port-Hamiltonian system (GP-PHS), a Gaussian process whose kernel is built from the system's interconnection and dissipation matrices so that its sample paths are port-Hamiltonian dynamics with an unknown Hamiltonian. It supplies the mean prediction and per-component variances that define the confidence set for the model error. The argument is carried by the modified matching equation for tracking, which makes the designed control assign the desired port-Hamiltonian error dynamics, and by the dissipation inequality (8), which lets the desired Hamiltonian $H_d$ serve as a Lyapunov function despite the perturbation $\eta$.
What would settle it
For the microactuator simulation, evaluate inequality (8) at finitely many states along the closed-loop trajectory, using the trained GP-PHS posterior variances with $\eta_i=\beta_i\operatorname{var}(\dot{x}_i|x,\mathcal{D})$ and the chosen $R_d$. If any sampled state violates the inequality, the conditions of Theorem 1 are not met at that point; a violation on a set of positive measure would refute the claimed probabilistic stability guarantee for that particular design.
Extended reading notes
Core claim
The paper's own terms: data-driven IDA-PBC can be taken from set-point regulation to reference tracking while keeping probabilistic stability and passivity guarantees. With a GP-PHS model, the exact matching equation is replaced by $\hat{G}^{\perp}\mu(\dot{x}|x,\mathcal{D})=\hat{G}^{\perp}\left([J_d-R_d]\nabla_{\bar{x}}H_d+\dot{x}_d\right)$, and the control law cancels the learned dynamics and assigns the desired port-Hamiltonian dynamics. The closed-loop error equation becomes $\dot{\bar{x}}=[J_d-R_d]\nabla_{\bar{x}}H_d+\eta$, with the GP model error $\eta$ bounded by $\beta_i\operatorname{var}(\dot{x}_i|x,\mathcal{D})$ with probability $1-p$. If $[\nabla_{\bar{x}}H_d]^\top\eta\le [\nabla_{\bar{x}}H_d]^\top R_d\nabla_{\bar{x}}H_d$, then $H_d$ is a Lyapunov function and the error equilibrium is stable; weaker forms of the same inequality give boundedness and semi-passivity.
Load-bearing premise
The load-bearing premise is condition (8) (and its relaxed form (15)): the injected damping $R_d$ must dominate the worst-case Gaussian-process model error $\eta$ for every state and every $\eta$ in the confidence set defined by $\beta_i\operatorname{var}(\dot{x}_i|x,\mathcal{D})$. This state-dependent inequality is not derived from the GP model or from $H_d$, is acknowledged as hard to satisfy, and is never verified in the simulation.
Editorial extensions
If this is right
- If Theorem 1 holds, a tracking controller can be designed without knowing the Hamiltonian, provided training data make the GP error small and the damping satisfies inequality (8).
- Stability is probabilistic with an explicit probability $1-p$, so data quantity and quality enter directly into the safety claim.
- Under the relaxed condition (15), the tracking error is ultimately bounded in a neighborhood of the desired trajectory, and the neighborhood size can be estimated by level-set methods.
- Corollary 2 means the closed-loop system is semi-passive with respect to an external input, allowing interconnection with other passive systems.
Reading between the lines
- The paper leaves open whether the modified matching equation (7) can be solved for $H_d$ and $R_d$ in general; the simulation chooses $H_d$ by hand as the posterior mean evaluated at $\bar{x}$, so the general solvability of the PDE is untested.
- A direct test of the load-bearing condition would be to compute $\beta_i\operatorname{var}(\dot{x}_i|x,\mathcal{D})$ along the simulated trajectory and check (8) pointwise; the paper does not report this check.
- The confidence set for $\eta$ is symmetric and componentwise state-independent; if the true model error is biased or heteroscedastic, inequalities (8) and (15) could fail even when the GP marginal variances look small.
- The approach suggests a design loop the paper does not pursue: jointly optimize $H_d$ and $R_d$ to satisfy both (7) and (8), which would make the damping condition easier to meet and the guarantees constructive.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. Summary: The paper proposes a data-driven passivity-based tracking controller for port-Hamiltonian systems with partially unknown dynamics. An unknown Hamiltonian is learned with a Gaussian process port-Hamiltonian system (GP-PHS) model, and a modified matching equation extends IDA-PBC to reference tracking. Theorem 1 states that if the desired Hamiltonian and damping satisfy the dissipation domination inequality (8), the tracking error equilibrium of the closed loop is stable with probability at least 1-p. Lemma 1 relaxes (8) to hold only outside an epsilon-ball, yielding probabilistic boundedness, and Corollary 2 gives semi-passivity under an external input. The approach is demonstrated on an electrostatic microactuator.
Significance. Significance: The paper addresses a genuine gap—extending IDA-PBC to tracking under unstructured uncertainty—and the use of GP-PHS is a principled way to carry uncertainty into the control design. The Lyapunov argument is transparent, and the relaxed result in Lemma 1 is a sensible fallback when the strong inequality (8) cannot be certified. The simulation illustrates the closed-loop behavior. However, the central stable-equilibrium claim is only conditional on (8), which is neither constructively designed nor verified; the demonstrated result is boundedness. With a verification procedure and a correction of the confidence-set scaling, the paper would make a solid contribution.
major comments (4)
- [Section IV, Theorem 1 and Lemma 1] The stability guarantee in Theorem 1 depends entirely on the dissipation domination condition (8), but the paper provides no design rule for choosing H_d and R_d that makes (8) hold, and Section IV never verifies it. The simulation is justified with Lemma 1, yet the relaxed condition (15) is also not evaluated: the paper does not report the epsilon for which (15) holds on the set X=[-2,2]^3, nor the resulting bound on norm(x-x_d). As a consequence, the numerical example supports boundedness but not the stable-equilibrium claim of Theorem 1. Please add a verification procedure for (8) or, for the relaxed claim, compute the region where (15) holds and the corresponding epsilon.
- [Section III-A, Theorem 1 and Lemma 1] The confidence set is stated as |eta_i(x)| <= beta_i var(xdot_i | x,D), but the standard GP error bound cited from [23] is of the form beta_i sigma_i(x), where sigma_i^2 = var(xdot_i | x,D). With the variance in place of the standard deviation, the admissible eta-set is mis-sized and the 'with probability at least 1-p' statement in Theorem 1 and Lemma 1 is not justified by the cited bound. Please replace var(xdot_i | x,D) by the posterior standard deviation, or provide a separate derivation for the variance-scaled bound.
- [Section III-B, proof of Theorem 1] In Eq. (10), the model is decomposed as xdot = mu(xdot | x,D) + Ghat(x)u + eta(x), but mu(xdot | x,D) is the posterior mean of the GP-PHS model (6), which already contains the input-dependent term Ghat(x)u. If mu includes Ghat u, the control input is double-counted in (11), and the matching equation (7) must be re-derived. Please define mu precisely (e.g., as the posterior mean of the Hamiltonian drift with Ghat u separated out) and rewrite (7) and (10)-(12) consistently.
- [Section II-B, Property 1 and Theorem 1] Property 1 requires only that H_d has its minimum at x=x_d, but the proof of Theorem 1 uses H_d as a Lyapunov function to conclude stability. Lyapunov stability at the equilibrium bar-x=0 also requires H_d to be positive definite with respect to bar-x (at least locally) and, for attractivity, to have an isolated minimum; these conditions are not stated or verified. Corollary 1 adds an isolated minimum only for asymptotic stability. Please state explicit positive-definiteness and regularity assumptions on H_d in Property 1.
minor comments (5)
- [Section III-B, Corollary 2] The phrase 'than the conditions in Lemma 1 are satisfied' should read 'then the conditions in Lemma 1 are satisfied.'
- [Section IV] The desired trajectory components x_d,2 and x_d,3 are defined only implicitly as solutions of (23); the paper should give the explicit expressions or the numerical solution procedure so that the simulation is reproducible.
- [Section IV, Eq. (22)] The relationship between the plant damping b and the desired-dynamics parameter bhat used in (22) should be stated; the text says the only change is 1/r_d, but the matrix also contains bhat.
- [Section III-A] The meaning of beta_i is left ambiguous; in [23] the confidence parameter depends on the information gain and the chosen failure probability. Please define beta_i explicitly and state how it is computed for the PHS kernel.
- [Figure 3 caption] The caption 'Top 1 and 2: Closed-loop system with the proposed tracking control law' is grammatically unclear; consider separating the two panels or rewording the caption.
Circularity Check
No significant circularity: the control-law derivation is conditional on explicitly stated matching and dissipation inequalities, and the main guarantees are not obtained by fitting a parameter to the claimed output.
full rationale
The paper's central claims are conditional theorems, not the result of fitting a parameter and then renaming the fit a prediction. In Theorem 1, the control law (9) is constructed from the GP-PHS mean model and the desired dynamics (Jd, Rd, Hd) that satisfy the modified matching equation (7). The stability argument uses Hd as a Lyapunov-like function; the derivative in (13) is non-positive precisely when the stated dissipation inequality (8) holds. That inequality is an explicit sufficient condition on the design variables (Hd, Rd) and the uncertainty set, not an artifact of the learning procedure. The relaxed Lemma 1 and Corollary 2 similarly rest on the explicit condition (15). The simulation section does not verify (8) or (15), and the paper itself acknowledges that (8) 'might be challenging to satisfy'; this is an evidentiary or rigor gap in the numerical demonstration, not circular reasoning. The paper adopts GP-PHS from the authors' prior work [19], [21], and cites [23] for Gaussian-process confidence bounds, but these citations are not used to define the target result or to forbid alternative designs; they supply a modeling framework and an uncertainty bound, which are independent inputs. The chosen desired Hamiltonian is a design choice, and no fitted constant is silently converted into the claimed stability conclusion. Therefore no circular step can be exhibited, and the appropriate finding is no significant circularity.
Assumptions & free parameters
free parameters (3)
- GP hyperparameters =
not reported
- Desired damping r_d =
10 (1/r_d = 10)
- GP bound constants β_i =
not reported
assumptions (4)
- domain assumption The physical system is a port-Hamiltonian system (1) with known structural matrices J, R, G up to unknown parameters.
- domain assumption The unknown Hamiltonian dynamics have bounded RKHS norm with respect to the PHS kernel, ||(J-R)∇H||_{k_phs} < ∞.
- ad hoc to paper The matching equation (7) has a solution Hd, Jd, Rd satisfying Property 1.
- ad hoc to paper The dissipation domination inequality (8) (or its relaxed version (15)) holds for the designed Hd and Rd and the GP confidence set.
Cite this review
Pith. "Pith review of Physics-informed Learning for Passivity-based Tracking Control." pith.science (2026). https://pith.science/paper/QQ2R24F4
@misc{pith2026250501569,
author = {Pith},
title = {Pith review of: Physics-informed Learning for Passivity-based Tracking Control},
year = {2026},
howpublished = {\url{https://pith.science/paper/QQ2R24F4}},
note = {Machine review of arXiv:2505.01569}
}
read the original abstract
Passivity-based control ensures system stability by leveraging dissipative properties and is widely applied in electrical and mechanical systems. Port-Hamiltonian systems (PHS), in particular, are well-suited for interconnection and damping assignment passivity-based control (IDA-PBC) due to their structured, energy-centric modeling approach. However, current IDA-PBC faces two key challenges: (i) it requires precise system knowledge, which is often unavailable due to model uncertainties, and (ii) it is typically limited to set-point control. To address these limitations, we propose a data-driven tracking control approach based on a physics-informed model, namely Gaussian process Port-Hamiltonian systems, along with the modified matching equation. By leveraging the Bayesian nature of the model, we establish probabilistic stability and passivity guarantees. A simulation demonstrates the effectiveness of our approach.
Figures
Reference graph
Works this paper leans on
-
[23]
Information-theoretic regret bounds for Gaussian process optimiza- tion in the bandit setting,
N. Srinivas, A. Krause, S. M. Kakade, and M. W. Seeger, “Information-theoretic regret bounds for Gaussian process optimiza- tion in the bandit setting,” IEEE Transactions on Information Theory , vol. 58, no. 5, pp. 3250–3265, 2012
work page 2012
-
[1]
Passivity-based control of nonlinear systems: A tutorial,
R. Ortega, Z. P. Jiang, and D. J. Hill, “Passivity-based control of nonlinear systems: A tutorial,” in Proceedings of the American Control Conference, vol. 5. IEEE, 1997, pp. 2633–2637
work page 1997
-
[2]
Hatanaka, N
T. Hatanaka, N. Chopra, M. Fujita, and M. W. Spong, Passivity-based control and estimation in networked robotics . Springer, 2015
2015
-
[3]
Passivity-based controllers for the stabilization of dc-to-dc power converters,
H. Sira-Ramirez, R. A. Perez-Moreno, R. Ortega, and M. Garcia- Esteban, “Passivity-based controllers for the stabilization of dc-to-dc power converters,” automatica, vol. 33, no. 4, pp. 499–513, 1997
work page 1997
-
[4]
J. A. Acosta, R. Ortega, A. Astolfi, and A. D. Mahindrakar, “Intercon- nection and damping assignment passivity-based control of mechan- ical systems with underactuation degree one,” IEEE Transactions on Automatic Control, vol. 50, no. 12, pp. 1936–1955, 2005
work page 1936
-
[5]
Van der Schaft, L2-gain and passivity techniques in nonlinear control
A. Van der Schaft, L2-gain and passivity techniques in nonlinear control. Springer, 2000
2000
-
[6]
R. Ortega, A. Van Der Schaft, B. Maschke, and G. Escobar, “Inter- connection and damping assignment passivity-based control of port- controlled Hamiltonian systems,” Automatica, vol. 38, no. 4, pp. 585– 596, 2002
work page 2002
-
[7]
Physical damping in ida- pbc controlled underactuated mechanical systems,
F. G ´omez-Estern and A. J. Van der Schaft, “Physical damping in ida- pbc controlled underactuated mechanical systems,” European Journal of control, vol. 10, no. 5, pp. 451–468, 2004
work page 2004
Show all 26 references
-
[8]
Interconnection and damping assignment passivity-based control: A survey,
R. Ortega and E. Garcia-Canseco, “Interconnection and damping assignment passivity-based control: A survey,” European Journal of control, vol. 10, no. 5, pp. 432–450, 2004
2004
-
[9]
Optimal task-invariant energetic control for a knee-ankle exoskeleton,
J. Lin, N. V . Divekar, G. Lv, and R. D. Gregg, “Optimal task-invariant energetic control for a knee-ankle exoskeleton,” in American Control Conference (ACC), 2021, pp. 5029–5034
2021
-
[10]
Modified energy-balancing-based control for the tracking problem,
Z. Wang and P. Goldsmith, “Modified energy-balancing-based control for the tracking problem,” IET Control Theory & Applications , vol. 2, no. 4, pp. 310–322, 2008
2008
-
[11]
Port- Hamiltonian systems in adaptive and learning control: A survey,
S. P. Nageshrao, G. A. Lopes, D. Jeltsema, and R. Babu ˇska, “Port- Hamiltonian systems in adaptive and learning control: A survey,”IEEE Transactions on Automatic Control , vol. 61, no. 5, pp. 1223–1238, 2015
2015
-
[12]
Simultaneous stabilization of a set of nonlinear port-controlled hamiltonian systems,
Y . Wang, G. Feng, and D. Cheng, “Simultaneous stabilization of a set of nonlinear port-controlled hamiltonian systems,” Automatica, vol. 43, no. 3, pp. 403–415, 2007
2007
-
[13]
Adaptive tracking control of fully actuated port-Hamiltonian mechanical systems,
D. A. Dirksz and J. M. Scherpen, “Adaptive tracking control of fully actuated port-Hamiltonian mechanical systems,” in IEEE International Conference on Control Applications , 2010, pp. 1678–1683
2010
-
[14]
Rein- forcement learning for port-Hamiltonian systems,
O. Sprangers, R. Babu ˇska, S. P. Nageshrao, and G. A. Lopes, “Rein- forcement learning for port-Hamiltonian systems,” IEEE Transactions on Cybernetics , vol. 45, no. 5, pp. 1017–1027, 2014
2014
-
[15]
Iterative learning control of Hamiltonian systems: I/o based optimal control approach,
K. Fujimoto and T. Sugie, “Iterative learning control of Hamiltonian systems: I/o based optimal control approach,” IEEE Transactions on Automatic Control, vol. 48, no. 10, pp. 1756–1761, 2003
2003
-
[16]
On iterative learning control of nonholonomic Hamil- tonian systems,
K. Fujimoto, “On iterative learning control of nonholonomic Hamil- tonian systems,” in Proc. 16th Symposium on Mathematical Theory of Networks and Systems (MTNS2004) . Citeseer, 2004
2004
-
[17]
A robust IDA-PBC approach for handling uncertainties in underactuated mechanical systems,
M. Ryalat and D. S. Laila, “A robust IDA-PBC approach for handling uncertainties in underactuated mechanical systems,” IEEE Transac- tions on Automatic Control , vol. 63, no. 10, pp. 3495–3502, 2018
2018
-
[18]
Hamiltonian-based neural ode networks on the se (3) manifold for dynamics learning and control,
T. Duong and N. Atanasov, “Hamiltonian-based neural ode networks on the se (3) manifold for dynamics learning and control,” arXiv preprint arXiv:2106.12782, 2021
2021 arXiv
-
[19]
Data-driven Bayesian control of port-Hamiltonian sys- tems,
T. Beckers, “Data-driven Bayesian control of port-Hamiltonian sys- tems,” in 62nd IEEE Conference on Decision and Control (CDC) . IEEE, 2023, pp. 8708–8713
2023
-
[20]
Port-Hamiltonian systems theory: An introductory overview,
A. Van Der Schaft and D. Jeltsema, “Port-Hamiltonian systems theory: An introductory overview,” F oundations and Trends in Systems and Control, vol. 1, no. 2-3, pp. 173–378, 2014
2014
-
[21]
Gaussian process port-Hamiltonian systems: Bayesian learning with physics prior,
T. Beckers, J. Seidman, P. Perdikaris, and G. J. Pappas, “Gaussian process port-Hamiltonian systems: Bayesian learning with physics prior,” in IEEE Conference on Decision and Control (CDC) , 2022, pp. 1447–1453
2022
-
[22]
C. E. Rasmussen and C. K. Williams, Gaussian processes for machine learning. MIT press Cambridge, 2006, vol. 1
2006
-
[24]
Stability of dynamical systems,
A. N. Michel, L. Hou, and D. Liu, “Stability of dynamical systems,” Birkha¨user , Boston, MA, 2008
2008
-
[25]
Passivity based design of synchronizing systems,
A. Y . Pogromsky, “Passivity based design of synchronizing systems,” International journal of Bifurcation and Chaos , vol. 8, no. 02, pp. 295–319, 1998
1998
-
[26]
Nonlinear dynamic output feedback stabilization of electrostatically actuated mems,
D. S. Maithripala, J. M. Berg, and W. P. Dayawansa, “Nonlinear dynamic output feedback stabilization of electrostatically actuated mems,” in 42nd IEEE International Conference on Decision and Control, vol. 1. IEEE, 2003, pp. 61–66
2003
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.