Calculating Domain of Attraction Boundary of Power Systems Based on the Gentlest Ascent Dynamics
Pith reviewed 2026-05-20 23:17 UTC · model grok-4.3
The pith
The domain of attraction boundary in power systems equals the closure of stable manifolds of index-1 critical elements.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors prove that the DOA boundary is the closure of the union of stable manifolds of index-1 critical elements and establish a stability theory for a perturbed GAD system. They employ GAD for 1-saddle points, adjoint operators for periodic orbits, and stable manifold algorithms to compute the boundary in benchmark power system models.
What carries the argument
Gentlest ascent dynamics for locating index-1 saddle points together with stable manifold construction algorithms and adjoint methods for periodic orbits.
If this is right
- DOA boundaries become computable by locating index-1 critical elements and building their stable manifolds.
- The method handles both saddle-point-only cases and cases containing periodic orbits.
- Transient stability after faults can be assessed by examining whether post-disturbance states lie inside the constructed region.
- Numerical tests on small generator models show the algorithms reproduce the expected boundary geometry.
Where Pith is reading between the lines
- The same manifold-tracing approach could apply to stability questions in other classes of nonlinear oscillators and networks.
- Stability results for the perturbed GAD system may improve understanding of numerical continuation methods in related dynamical problems.
- Scaling the algorithms to higher-dimensional models could support stability screening in grids that include inverter-based resources.
Load-bearing premise
Certain assumptions under which the domain of attraction boundary equals the closure of the union of stable manifolds of index-1 critical elements.
What would settle it
A concrete power system trajectory that starts inside the computed manifold boundary yet diverges from the stable equilibrium, or one that starts outside yet converges.
Figures
read the original abstract
The power system, a fundamental public utility, is increasingly important due to growing global electricity demand. Recent large-scale blackouts (e.g., Iberian Peninsula, UK) have raised concerns about transient stability under impact faults. Transient stability is determined by post-disturbance synchronizing capability of synchronous generators, formulated as identifying the domain of attraction (DOA) boundary of the asymptotically stable equilibrium. Using a benchmark model of synchronous-generator-dominated power systems, this report employs a gentlest ascent dynamics (GAD) method for 1-saddle points, an adjoint operator method for periodic orbits, and stable manifold algorithms to compute the DOA boundary. These algorithms transform DOA boundary determination into constructing unstable critical elements (saddle points and periodic orbits) and their stable manifolds. Theoretically, under certain assumptions we prove that the DOA boundary is the closure of the union of stable manifolds of index-1 critical elements, and establish a stability theory for a perturbed GAD system. Numerical experiments on two-machine and three-machine systems (with only saddle points or with periodic orbits) validate the effectiveness and accuracy. Results show the algorithms accurately capture the geometric structure of the DOA boundary, providing a new numerical tool for transient stability analysis.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes using the gentlest ascent dynamics (GAD) method to locate index-1 saddle points, an adjoint operator approach for periodic orbits, and stable manifold algorithms to compute the domain of attraction (DOA) boundary for the post-fault stable equilibrium in synchronous-generator-dominated power system models. It claims a theoretical result that, under certain assumptions, the DOA boundary equals the closure of the union of stable manifolds of index-1 critical elements, along with a stability theory for a perturbed GAD system, and reports numerical validation on two- and three-machine benchmark systems (with and without periodic orbits).
Significance. If the unspecified assumptions hold and the numerical accuracy is confirmed, the work supplies a new computational framework for transient stability analysis in power systems by reducing DOA boundary determination to the construction of unstable critical elements and their stable manifolds. This could be useful for large-scale blackout prevention, especially if the GAD-based search scales better than existing methods for locating saddles and periodic orbits in high-dimensional swing equations.
major comments (2)
- [Abstract and theoretical development] Abstract and theoretical development: The central claim that the DOA boundary is the closure of the union of stable manifolds of index-1 critical elements (saddles and periodic orbits) is proved only under 'certain assumptions.' These assumptions are not enumerated, so it is impossible to verify whether they hold for the two-machine and three-machine models, particularly the case containing periodic orbits. This is load-bearing because the justification for using GAD to locate the critical elements and then compute their stable manifolds rests entirely on this reduction.
- [Numerical experiments section] Numerical experiments section: The reported validation on two- and three-machine systems asserts that the algorithms 'accurately capture the geometric structure' of the DOA boundary, yet no quantitative error metrics, baseline comparisons against other DOA estimation techniques, or details on post-disturbance trajectory generation are provided. Without these, the empirical support for accuracy and practical utility remains incomplete.
minor comments (1)
- [Abstract] The abstract refers to 'a benchmark model of synchronous-generator-dominated power systems' without giving the explicit differential equations or parameter values; these should be stated early so readers can reproduce the critical-element computations.
Simulated Author's Rebuttal
We thank the referee for the insightful comments on our manuscript. We address the major comments point by point below, indicating where revisions will be made to strengthen the paper.
read point-by-point responses
-
Referee: [Abstract and theoretical development] The central claim that the DOA boundary is the closure of the union of stable manifolds of index-1 critical elements (saddles and periodic orbits) is proved only under 'certain assumptions.' These assumptions are not enumerated, so it is impossible to verify whether they hold for the two-machine and three-machine models, particularly the case containing periodic orbits. This is load-bearing because the justification for using GAD to locate the critical elements and then compute their stable manifolds rests entirely on this reduction.
Authors: We agree that explicitly listing the assumptions would improve clarity and allow readers to verify their applicability. The proof in the manuscript is based on assumptions including the hyperbolicity of the index-1 critical elements (saddles and periodic orbits) and the transversality of their stable manifolds to the DOA boundary. We will revise the manuscript by adding a subsection that enumerates these assumptions clearly and discusses their validity for the benchmark systems, including the three-machine system with periodic orbits. revision: yes
-
Referee: [Numerical experiments section] The reported validation on two- and three-machine systems asserts that the algorithms 'accurately capture the geometric structure' of the DOA boundary, yet no quantitative error metrics, baseline comparisons against other DOA estimation techniques, or details on post-disturbance trajectory generation are provided. Without these, the empirical support for accuracy and practical utility remains incomplete.
Authors: We acknowledge the need for more rigorous quantitative validation. In the revised manuscript, we will include quantitative error metrics, such as the maximum deviation of the computed boundary from reference trajectories obtained via numerical integration of the system dynamics. We will also add comparisons with baseline methods like the energy function approach or closest UEP method where applicable, and provide details on the post-disturbance trajectory generation process used for validation. This will better demonstrate the accuracy and utility of the proposed framework. revision: yes
Circularity Check
No significant circularity; theoretical claim applies standard dynamical systems results to power-system model
full rationale
The paper's strongest claim is a proof, under unspecified assumptions, that the DOA boundary equals the closure of the union of stable manifolds of index-1 critical elements, together with a stability theory for a perturbed GAD system. No equation or derivation in the abstract or described content reduces this result to a fitted quantity, self-citation chain, or renamed input; the GAD and adjoint-operator constructions are presented as independent numerical tools that locate the critical elements whose manifolds are then computed. The work therefore remains self-contained against external dynamical-systems benchmarks, with the 'certain assumptions' serving as a standard caveat rather than a hidden self-definition.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Under certain assumptions the DOA boundary is the closure of the union of stable manifolds of index-1 critical elements
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Theoretically, under certain assumptions we prove that the DOA boundary is the closure of the union of stable manifolds of index-1 critical elements, and establish a stability theory for a perturbed GAD system.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
GAD method … augmented dynamical system … 1-saddle point … asymptotically stable equilibrium point of an augmented system.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Perko.Differential Equations and Dynamical Systems, volume 7 ofTexts in Applied Mathematics
L. Perko.Differential Equations and Dynamical Systems, volume 7 ofTexts in Applied Mathematics. Springer, third edition, 2000. doi:10.1007/978-1-4613-0003-8
-
[2]
Kundur.Power System Stability and Control
P. Kundur.Power System Stability and Control. McGraw-hill, New York, 1994. ISBN 0-07-035958-X
work page 1994
-
[3]
R. Genesio, M. Tartaglia, and A. Vicino. On the estimation of asymptotic stability regions: State of the art and new proposals.IEEE Transactions on Automatic Control, 30(8):747–755, 1985
work page 1985
-
[4]
H. D. Chiang, M. W. Hirsch, and F. F. Wu. Stability regions of nonlinear au- tonomous dynamical systems.IEEE Transactions on Automatic Control, 33(1): 16–27, January 1988. ISSN 1558-2523. doi:10.1109/9.357
-
[5]
H. D. Chiang.Direct Methods for Stability Analysis of Electric Power Systems. Wiley, 2011. ISBN 978-0-470-48440-1. doi:10.1002/9780470872130
-
[6]
B. Stott. Power system dynamic response calculations.Proceedings of the IEEE, 67(2):219–241, 1979
work page 1979
-
[7]
M. K. Bakhshizadeh, S. Ghosh, G. Y. Yang, and L. Kocewika. Transient stability analysis of grid-connected converters in wind turbine systems based on linear lya- punov function and reverse-time trajectory.Journal of Modern Power Systems and Clean Energy, 12(3):782–790, 2024
work page 2024
- [8]
-
[9]
M. Kabalan, P. Singh, and D. Niebur. Nonlinear lyapunov stability analysis of seven models of a dc/ac droop controlled inverter connected to an infinite bus. IEEE Transactions on Smart Grid, 10(1):772–781, 2019
work page 2019
-
[10]
A. A. Fouad and V. Vittal.Power System Transient Stability Analysis Using the Transient Energy Function Method. Prentice Hall, 1991
work page 1991
-
[11]
M. A. Pai.Energy Function Analysis for Power System Stability. Springer, 1989. 35
work page 1989
-
[12]
C. C. Chu and H. D. Chiang. Constructing analytical energy functions for loss- less network-reduction power system models: Framework and new developments. Circuits, Systems and Signal Processing, 18(1):1–16, 1999
work page 1999
-
[13]
C. C. Chu and H. D. Chiang. Constructing analytical energy functions for network- preserving power system models.Circuits, Systems and Signal Processing, 24(4): 363–383, 2005
work page 2005
- [14]
-
[15]
M. A. Pai and S. Varwandkar. On the inclusion of transfer conductances in lya- punov functions for multimachine power systems.IEEE Transactions on Automatic Control, 22(6):983–985, 1977
work page 1977
- [16]
-
[17]
A. H. El-Abiad and K. Nagappan. Transient stability regions of multimachine power systems.IEEE Transactions on Power Apparatus and Systems, PAS-85(2):169–179, 1966
work page 1966
-
[18]
H. D. Chiang and J. S. Thorp. The closest unstable equilibrium point method for power system dynamic security assessment.IEEE Transactions on Circuits and Systems, 36(9):1187–1200, 1989
work page 1989
-
[19]
H. D. Chiang, F. Wu, and P. Varaiya. Foundations of direct methods for power system transient stability analysis.IEEE Transactions on Circuits and Systems, 34 (2):160–173, 1987
work page 1987
-
[20]
H. D. Chiang, F. F. Wu, and P. P. Varaiya. A bcu method for direct analysis of power system transient stability.IEEE Transactions on Power Systems, 9(3): 1194–1208, 1994
work page 1994
-
[21]
N. Kakimoto, Y. Ohnogi, H. Matsuda, and H. Shibuya. Transient stability analysis of large-scale power system by lyapunov’s direct method.IEEE Transactions on Power Apparatus and Systems, PAS-103(1):160–167, 1984
work page 1984
-
[22]
H. D. Chiang, F. F. Wu, and P. P. Varaiya. Foundations of the potential en- ergy boundary surface method for power system transient stability analysis.IEEE Transactions on Circuits and Systems, 35(6):712–728, 1988
work page 1988
-
[23]
P. A. Parrilo.Structured Semidefinite Programs and Semialgebraic Geometry Meth- ods in Robustness and Optimization. PhD thesis, California Institute of Technology, 2000. 36
work page 2000
-
[24]
A. Papachristodoulou and S. Prajna. A tutorial on sum of squares techniques for systems analysis. InProceedings of the 2005 American Control Conference, pages 2686–2700, 2005
work page 2005
-
[25]
M. Anghel, F. Milano, and A. Papachristodoulou. Algorithmic construction of lya- punov functions for power system stability analysis.IEEE Transactions on Circuits and Systems I: Regular Papers, 60(9):2533–2546, September 2013. ISSN 1558-0806. doi:10.1109/tcsi.2013.2246233
-
[26]
S. Izumi, H. Somekawa, X. Xin, and T. Yamasaki. Estimation of regions of attrac- tion of power systems by using sum of squares programming.Electrical Engineering, 100(4):2205–2216, May 2018. ISSN 1432-0487. doi:10.1007/s00202-018-0690-z
-
[27]
Y. Liu, Z. Chen, H. Yao, L. Yi, and Q. H. Wu. Estimating critical clearing time of grid faults using da of state-reduction model of power sys- tems.CSEE Journal of Power and Energy Systems, 10(2):807–820, April 2024. doi:10.17775/CSEEJPES.2022.07170
-
[28]
Chesi.Domain of Attraction: Analysis and Control via SOS Programming
G. Chesi.Domain of Attraction: Analysis and Control via SOS Programming. Springer, 2011
work page 2011
-
[29]
T. Takagi and M. Sugeno. Fuzzy identification of systems and its applications to modeling and control.IEEE Transactions on Systems, Man, and Cybernetics, SMC-15(1):116–132, 1985
work page 1985
-
[30]
T. M. Guerra and L. Vermeiren. Lmi-based relaxed nonquadratic stabilization conditions for nonlinear systems in the takagi-sugeno’s form.Automatica, 40(5): 823–829, 2004
work page 2004
-
[31]
X. L. Li, Z. Wang, Y. X. Liu, Z. G. Wang, L. Zhu, L. Guo, C. Zhang, J. B. Zhu, and C. S. Wang. The largest estimated domain of attraction and its applications for transient stability analysis of pll synchronization in weak-grid-connected vscs. IEEE Transactions on Power Systems, 38(5):4107–4121, 2023
work page 2023
-
[32]
W. Z. Tang, B. Li, S. C. Hou, X. Q. Shao, and H. J. Yu. Transient stability analysis for grid-connected renewable power generation systems based on lmi optimization modelling.Electronics, 13(24):5052, 2024
work page 2024
-
[33]
Y. W. Ding, F. Gao, and M. M. Khan. Transient stability analysis of microgrid con- sidering impact of grid-following converter’s current controller.IEEE Transactions on Power Electronics, 39(8):9100–9105, 2024
work page 2024
-
[34]
Y. C. Chang, N. Roohi, and S. C. Gao. Neural lyapunov control. InProceedings of the 33rd International Conference on Neural Information Processing Systems, pages 3245–3254, 2019. 37
work page 2019
- [35]
-
[36]
T. Q. Zhao, J. H. Wang, X. N. Lu, and Y. H. Du. Neural lyapunov control for power system transient stability: a deep learning-based approach.IEEE Transactions on Power Systems, 37(2):955–966, 2022
work page 2022
-
[37]
Y. F. Liu, L. J. Zhou, S. Chen, J. R. Zhang, and Y. Wang. An improved neural lyapunov method for transient stability assessment of networked microgrids.IEEE Transactions on Smart Grid, 15(2):1410–1422, 2024
work page 2024
-
[38]
S. T. Zhang, Z. X. Zhu, and Y. Li. A critical review of data-driven transient stability assessment of power systems: principles, prospects and challenges.Energies, 14(21): 7238, 2021
work page 2021
-
[39]
stability regions of nonlinear autonomous dynamical systems
M. W. Fisher and I. A. Hiskens. Comments on “stability regions of nonlinear autonomous dynamical systems”.IEEE Transactions on Automatic Control, 66 (12):6194–6196, December 2021. ISSN 2334-3303. doi:10.1109/tac.2021.3061674
-
[40]
M. W. Fisher and I. A. Hiskens. Hausdorff continuity of region of attraction bound- ary under parameter variation with application to disturbance recovery.SIAM Jour- nal on Applied Dynamical Systems, 21(1):327–365, 2022. doi:10.1137/20M1371944
-
[41]
W. N. E and X. Zhou. The gentlest ascent dynamics.Nonlinearity, 24(6):1831, May 2011. doi:10.1088/0951-7715/24/6/008
-
[42]
H. D. Chiang and L. F. C. Alberto.Stability Regions of Nonlinear Dynamical Systems: Theory, Estimation, and Applications. Cambridge, 2015. ISBN 978-1- 107-03540-9. doi:10.1017/CBO9781139548861
-
[43]
S. Azimi, O. Ashtari, and T. M. Schneider. Constructing periodic orbits of high- dimensional chaotic systems by an adjoint-based variational method.Phys. Rev. E, 105:014217, Jan 2022. doi:10.1103/PhysRevE.105.014217
-
[44]
Y. F. Tang, R. L. Zhang, B. B. Zhu, A. Q. Zhu, and S. X. Wu.Geometric Al- gorithms for Dynamical Systems with Applications (in Chinese), volume 215 of Modern Mathematical Basic Series. Science Press, 2026. ISBN 9787030846334
work page 2026
-
[45]
J. M. Lee.Introduction to Smooth Manifolds, volume 218 ofGraduate Texts in Math- ematics. Springer, second edition, 2013. ISBN 978-1-4419-9981-8. doi:10.1007/978- 1-4419-9982-5
-
[46]
S. T. Gu and X. P. Zhou. Simplified gentlest ascent dynamics for saddle points in non-gradient systems.Chaos: An Interdisciplinary Journal of Nonlinear Science, 28(12):123106, December 2018. doi:10.1063/1.5046819. 38
-
[47]
Yoshua Bengio, Aaron Courville, and Pascal Vincent. Representation learning: a review and new perspectives. 35(8):1798–1828. doi:10.1109/TPAMI.2013.50
-
[48]
Gradient-based learning applied to document recognition,
Yann LeCun, L´ eon Bottou, Yoshua Bengio, and Patrick Haffner. Gradient-based learning applied to document recognition. 86(11):2278–2324. doi:10.1109/5.726791
-
[49]
Dahl, Abdel-rahman Mohamed, Navdeep Jaitly, Andrew Senior, Vincent Vanhoucke, Patrick Nguyen, Tara N
Geoffrey Hinton, Li Deng, Dong Yu, George E. Dahl, Abdel-rahman Mohamed, Navdeep Jaitly, Andrew Senior, Vincent Vanhoucke, Patrick Nguyen, Tara N. Sainath, and Brian Kingsbury. Deep neural networks for acoustic modeling in speech recognition. 29(6):82–97. doi:10.1109/MSP.2012.2205597
-
[50]
Dropout: a simple way to prevent neural networks from overfitting
Nitish Srivastava, Geoffrey Hinton, Alex Krizhevsky, Ilya Sutskever, and Ruslan Salakhutdinov. Dropout: a simple way to prevent neural networks from overfitting. 15(1):1929–1958
work page 1929
-
[51]
Alex Krizhevsky, Ilya Sutskever, and Geoffrey E. Hinton. Imagenet classification with deep convolutional neural networks. 25:1090–1098. 39
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.