pith. sign in

arxiv: 2605.22171 · v1 · pith:Y2G2FX6Gnew · submitted 2026-05-21 · 📡 eess.SY · cs.SY

Equilibrium-Free Contraction Stability Analysis for Grid-Forming Converter-Based Microgrids

Shijie Peng , Xiuqiang He , Xi Ru , Hua Geng This is my paper

Pith reviewed 2026-05-22 04:32 UTC · model grok-4.3

classification 📡 eess.SY cs.SY
keywords grid-forming convertersmicrogrid stabilityequilibrium-free analysiscontraction theoryJacobian decompositionprojected state spacepower system dynamicstrajectory guarantees
0
0 comments X

The pith

Symmetry-aware projection and blockwise Jacobian decomposition certify microgrid stability without locating any equilibrium point.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops an equilibrium-free contraction stability method for renewable microgrids dominated by grid-forming converters that experience constant power fluctuations. It removes the rotational symmetry from uniform angle shifts by reformulating the dynamics in a projected state space. A blockwise decomposition of the Jacobian then isolates the coupled active and reactive power flows to produce a computable regional contraction condition. This condition converts directly into forward-invariant certificates that guarantee trajectory behavior for both undisturbed operation and cases with disturbances. Readers care because conventional stability tests require knowing the operating point in advance, which is rarely feasible when injections vary continuously.

Core claim

By formulating the system in a symmetry-aware projected state space that removes the intrinsic rotational mode induced by uniform angle shifts, the paper introduces a blockwise Jacobian decomposition to characterize coupled active and reactive power dynamics. This yields a computable regional contraction condition that is converted into forward-invariant stability certificates. For autonomous operation without disturbances the method provides an equilibrium-free nonlinear stability characterization together with an estimation of the region of attraction. For non-autonomous operation under disturbances it derives explicit bounds for quasi-steady tracking under slowly varying injections andfor

What carries the argument

Blockwise Jacobian decomposition applied inside the symmetry-aware projected state space, which isolates coupled active-reactive dynamics and produces a regional contraction condition convertible to forward-invariant certificates.

If this is right

  • Equilibrium-free nonlinear stability characterization for autonomous microgrid operation.
  • Estimation of the region of attraction for the undisturbed case.
  • Explicit bounds on quasi-steady tracking under slowly varying power injections.
  • Robustness bounds under fast or composite disturbances.
  • Validation on a 9-bus test system.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same projection-plus-decomposition pattern could scale to larger converter-dominated networks by treating subsystems as blocks.
  • Control design that enlarges the certified region of attraction becomes directly testable against the contraction metric.
  • The approach may transfer to other rotationally symmetric systems such as multi-machine swing dynamics.
  • Comparing the size of the certified region against numerically computed basins under step-load changes would quantify conservatism.

Load-bearing premise

The system dynamics admit a symmetry-aware projected state-space formulation that removes the intrinsic rotational mode induced by uniform angle shifts while preserving the coupled active-reactive power behavior.

What would settle it

A simulation on the 9-bus system in which a trajectory leaves the forward-invariant set defined by the computed contraction condition, or remains inside the set when the regional condition is violated, would show whether the contraction condition truly implies the claimed stability certificates.

Figures

Figures reproduced from arXiv: 2605.22171 by Hua Geng, Shijie Peng, Xi Ru, Xiuqiang He.

Figure 1
Figure 1. Figure 1: Schematic of the droop-controlled islanded microgrid. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Semi-contraction and convergence in the R-projected space. The gray dashed balls indicate a distance bound that decays over time. C. Semi-Contraction Preliminaries Contraction theory offers an equilibrium-free view of stabil￾ity by studying convergence between neighboring trajectories rather than toward a fixed equilibrium. Such convergence can be certified through the matrix measure of the system Jacobian… view at source ↗
Figure 3
Figure 3. Figure 3: Impact of heterogeneity on cθ. The droop gains are parameterized as mp(ηh) = ¯mp1 + ηh(m0 p − m¯ p1), where m¯ p := 1 N 1⊤m0 p is the mean droop gain. The baseline case is ηh = 1 (dashed line), and ηh = 0 gives the uniform-droop case. 1) Active Power–Angle Block Characterization: The block Sθθ characterizes the synchronizing action in the active power– angle channel. Although Jθθ(x, t) has zero row sums, i… view at source ↗
Figure 2
Figure 2. Figure 2: within a forward-invariant set, projected distances be [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: Unified framework for contraction stability analysis. [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Reduced 3-bus converter network derived from the modified IEEE [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Time-domain verification for the autonomous case. (a), (b) Projected [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 6
Figure 6. Figure 6: Invariant-set verification for the autonomous case. (a) Invariant ball in [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: Time-domain verification under slowly varying disturbances, with [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Time-domain verification under composite disturbances, with [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
read the original abstract

Renewable-driven microgrids dominated by grid-forming (GFM) converters are subject to persistent power fluctuations, making equilibrium-known stability assessments restrictive. This paper develops an equilibrium-free contraction stability method based on semi-contraction theory. By formulating the system in a symmetry-aware projected state space, the intrinsic rotational mode induced by uniform angle shifts is removed. A blockwise Jacobian decomposition is introduced to characterize the coupled active and reactive power dynamics, yielding a computable regional contraction condition. This condition is then converted into forward-invariant stability certificates that provide trajectory-level performance guarantees. For autonomous operation without disturbances, the method provides an equilibrium-free nonlinear stability characterization together with an estimation of the region of attraction (ROA). For non-autonomous operation under disturbances, it derives explicit bounds for quasi-steady tracking under slowly varying injections and for robustness under fast or composite disturbances. Case studies on a 9-bus system validate the proposed method.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper develops an equilibrium-free contraction stability method for grid-forming converter-based microgrids based on semi-contraction theory. It formulates the system in a symmetry-aware projected state space that removes the intrinsic rotational mode induced by uniform angle shifts, introduces a blockwise Jacobian decomposition to characterize coupled active-reactive power dynamics, and derives a computable regional contraction condition. This condition is converted into forward-invariant stability certificates that provide trajectory-level performance guarantees. For autonomous operation, the method yields an equilibrium-free nonlinear stability characterization and an estimate of the region of attraction; for non-autonomous operation under disturbances, it provides explicit bounds for quasi-steady tracking under slowly varying injections and for robustness under fast or composite disturbances. The approach is illustrated on a 9-bus system.

Significance. If the projection step and subsequent derivations hold, the work supplies a practical tool for stability assessment in renewable-dominated microgrids where equilibria are not known a priori due to fluctuating injections. The trajectory-level guarantees and explicit disturbance bounds distinguish it from equilibrium-dependent methods and could support certification of GFM converter systems. The combination of semi-contraction theory with symmetry-aware projection is a notable technical contribution if the preservation of coupled dynamics is rigorously established.

major comments (2)
  1. [Formulation / symmetry-aware projection] The symmetry-aware projected state-space formulation is presented as removing the rotational mode while exactly preserving the coupled active-reactive power behavior (abstract and formulation section). However, the manuscript must supply an explicit construction of the projection operator, a proof that the projected vector field remains equivalent for stability purposes, and verification that no residual coupling is introduced that would invalidate transfer of the contraction condition back to original coordinates. This step is load-bearing for the central claim.
  2. [Blockwise Jacobian decomposition] § on blockwise Jacobian decomposition: the claim that the decomposition yields a computable regional contraction condition requires explicit equations for the blocks, the resulting contraction metric or rate, and any error bounds or data-exclusion rules used in the 9-bus validation. Without these, it is not possible to confirm that the condition is independent of equilibrium knowledge and that the forward-invariant certificates follow directly.
minor comments (2)
  1. [Abstract] The abstract is concise but could include one or two quantitative outcomes from the 9-bus case studies (e.g., ROA size or tracking error bounds) to give readers an immediate sense of practical performance.
  2. [Notation and figures] Notation for the projected state variables and the decomposition blocks should be introduced with a short table or diagram to improve readability for readers unfamiliar with semi-contraction theory.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough and constructive review. The comments help clarify the presentation of the symmetry-aware projection and blockwise decomposition. We address each major comment below and will incorporate the requested details in the revised manuscript.

read point-by-point responses

  1. Referee: [Formulation / symmetry-aware projection] The symmetry-aware projected state-space formulation is presented as removing the rotational mode while exactly preserving the coupled active-reactive power behavior (abstract and formulation section). However, the manuscript must supply an explicit construction of the projection operator, a proof that the projected vector field remains equivalent for stability purposes, and verification that no residual coupling is introduced that would invalidate transfer of the contraction condition back to original coordinates. This step is load-bearing for the central claim.

    Authors: We agree that greater explicitness will strengthen the manuscript. The projection operator is constructed as the orthogonal projector onto the subspace orthogonal to the all-ones vector in angle coordinates, removing the uniform rotational mode while leaving relative angles and frequencies unchanged. In the revision we will add the explicit matrix form of the projector, a short lemma establishing that the projected vector field is equivalent for contraction analysis (the rotational mode is neutrally stable and decoupled from power-flow dynamics), and a direct verification that the block structure of the Jacobian is preserved with no spurious cross terms. These additions will make the transfer of the contraction condition back to original coordinates fully rigorous. revision: yes

  2. Referee: [Blockwise Jacobian decomposition] § on blockwise Jacobian decomposition: the claim that the decomposition yields a computable regional contraction condition requires explicit equations for the blocks, the resulting contraction metric or rate, and any error bounds or data-exclusion rules used in the 9-bus validation. Without these, it is not possible to confirm that the condition is independent of equilibrium knowledge and that the forward-invariant certificates follow directly.

    Authors: We will supply the missing explicit material. The Jacobian is partitioned into four blocks (active-active, active-reactive, reactive-active, reactive-reactive) whose entries are written in closed form using the network admittance matrix and voltage magnitudes. The regional contraction condition is obtained by requiring that the maximum eigenvalue of a suitably weighted combination of these blocks remains negative inside a forward-invariant set; the metric is a diagonal weighting matrix whose entries are chosen from local voltage bounds. In the revision we will state the four block equations, the resulting contraction rate, and the precise data-exclusion rule (voltage magnitudes outside a 0.9–1.1 p.u. band are excluded from the regional certificate) used for the 9-bus example. This will confirm both equilibrium independence and the direct derivation of the forward-invariant certificates. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation relies on external semi-contraction framework

full rationale

The paper develops an equilibrium-free method based on semi-contraction theory, an external established framework. It introduces a symmetry-aware projected state-space formulation to remove the rotational mode and applies blockwise Jacobian decomposition to obtain a regional contraction condition, which is converted to forward-invariant certificates. These steps are presented as constructive methodological advances without reducing by the paper's own equations to fitted parameters, self-definitions, or self-citation chains. No load-bearing uniqueness theorems or ansatzes from prior author work are invoked in a circular manner. The approach remains self-contained with independent content for trajectory guarantees in microgrid systems.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The method depends on background results from semi-contraction theory and on the modeling choice that the rotational symmetry can be projected out without loss of the essential dynamics.

axioms (1)
  • domain assumption Semi-contraction theory supplies sufficient conditions for regional contraction and forward invariance
    Invoked to convert the Jacobian condition into stability certificates

pith-pipeline@v0.9.0 · 5692 in / 1255 out tokens · 62965 ms · 2026-05-22T04:32:49.741993+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

28 extracted references · 28 canonical work pages

  1. [1]

    Research roadmap on grid-forming inverters,

    Y . Linet al., “Research roadmap on grid-forming inverters,” National Renewable Energy Laboratory (NREL), Technical Report NREL/TP- 5D00-73476, 2020

    work page 2020

  2. [2]

    A review of power system flexibility with high penetration of renewables,

    B. Mohandes, M. S. E. Moursi, N. Hatziargyriou, and S. E. Khatib, “A review of power system flexibility with high penetration of renewables,” IEEE Transactions on Power Systems, vol. 34, no. 4, pp. 3140–3155, 2019

    work page 2019

  3. [3]

    Short term fluctuations of wind and solar power systems,

    M. Anvari, G. Lohmann, M. W ¨achter, P. Milan, E. Lorenz, D. Heine- mann, M. R. R. Tabar, and J. Peinke, “Short term fluctuations of wind and solar power systems,”New Journal of Physics, vol. 18, no. 6, p. 063027, jun 2016

    work page 2016

  4. [4]

    Kundur,Power System Stability and Control

    P. Kundur,Power System Stability and Control. New York, NY , USA: McGraw-Hill, 1994

    work page 1994

  5. [5]

    Chiang,Direct Methods for Stability Analysis of Electric Power Systems: Theoretical Foundation, BCU Methodologies, and Applica- tions, 2nd ed

    H.-D. Chiang,Direct Methods for Stability Analysis of Electric Power Systems: Theoretical Foundation, BCU Methodologies, and Applica- tions, 2nd ed. Hoboken, NJ, USA: Wiley, 2011

    work page 2011

  6. [6]

    Foundations and challenges of low-inertia systems (invited paper),

    F. Milano, F. D ¨orfler, G. Hug, D. J. Hill, and G. Verbi ˇc, “Foundations and challenges of low-inertia systems (invited paper),” in2018 Power Systems Computation Conference (PSCC), 2018, pp. 1–25

    work page 2018

  7. [7]

    Conditions for stability of droop-controlled inverter-based microgrids,

    J. Schiffer, R. Ortega, A. Astolfi, J. Raisch, and T. Sezi, “Conditions for stability of droop-controlled inverter-based microgrids,”Automatica, vol. 50, no. 10, pp. 2457–2469, 10 2014

    work page 2014

  8. [8]

    Synchronization and transient stability in power networks and nonuniform kuramoto oscillators,

    F. D ¨orfler and F. Bullo, “Synchronization and transient stability in power networks and nonuniform kuramoto oscillators,”SIAM Journal on Control and Optimization, vol. 50, no. 3, pp. 1616–1642, 2012

    work page 2012

  9. [9]

    Microgrid stability definitions, analysis, and examples,

    M. Farrokhabadi, C. A. Ca ˜nizares, J. W. Simpson-Porco, E. Nasr, L. Fan, P. A. Mendoza-Arayaet al., “Microgrid stability definitions, analysis, and examples,”IEEE Transactions on Power Systems, vol. 35, no. 1, pp. 13–29, 2020

    work page 2020

  10. [10]

    High-fidelity model order reduction for microgrids stability assess- ment,

    P. V orobev, P.-H. Huang, M. Al Hosani, J. L. Kirtley, and K. Turitsyn, “High-fidelity model order reduction for microgrids stability assess- ment,”IEEE Transactions on Power Systems, vol. 33, no. 1, pp. 874–887, 2018

    work page 2018

  11. [11]

    A partial stability approach to the problem of transient power system stability,

    J. L. Willems, “A partial stability approach to the problem of transient power system stability,”International Journal of Control, vol. 19, no. 1, pp. 1–14, 1974

    work page 1974

  12. [12]

    Power systems as dynamic networks,

    D. Hill and G. Chen, “Power systems as dynamic networks,” in2006 IEEE International Symposium on Circuits and Systems (ISCAS), 2006, pp. 4 pp.–725

    work page 2006

  13. [13]

    Partial stability and control: The state-of-the-art and development prospects,

    V . I. V orotnikov, “Partial stability and control: The state-of-the-art and development prospects,”Automation and Remote Control, vol. 66, pp. 511–561, 2005

    work page 2005

  14. [14]

    N. P. Bhatia and G. P. Szeg ¨o,Stability theory of dynamical systems. Berlin, Heidelberg: Springer-Verlag, 2002

    work page 2002

  15. [15]

    Augmented synchronization of power systems,

    P. Yang, F. Liu, T. Liu, and D. J. Hill, “Augmented synchronization of power systems,”IEEE Transactions on Automatic Control, vol. 69, no. 6, pp. 3673–3688, 2024

    work page 2024

  16. [16]

    Synchronization in complex oscillator networks and smart grids,

    F. D ¨orfler, M. Chertkov, and F. Bullo, “Synchronization in complex oscillator networks and smart grids,”Proceedings of the National Academy of Sciences, vol. 110, no. 6, pp. 2005–2010, 2013

    work page 2005

  17. [17]

    Synchronization and power sharing for droop-controlled inverters in islanded microgrids,

    J. W. Simpson-Porco, F. D ¨orfler, and F. Bullo, “Synchronization and power sharing for droop-controlled inverters in islanded microgrids,” Automatica, vol. 49, no. 9, pp. 2603–2611, 2013

    work page 2013

  18. [18]

    Stability analysis of power systems: A network synchronization perspective,

    L. Zhu and D. J. Hill, “Stability analysis of power systems: A network synchronization perspective,”SIAM Journal on Control and Optimiza- tion, vol. 56, no. 3, pp. 1640–1664, 2018

    work page 2018

  19. [19]

    On contraction analysis for non-linear systems,

    W. LOHMILLER and J.-J. E. SLOTINE, “On contraction analysis for non-linear systems,”Automatica, vol. 34, no. 6, pp. 683–696, 1998

    work page 1998

  20. [20]

    Bullo,Contraction Theory for Dynamical Systems, 1.3 ed

    F. Bullo,Contraction Theory for Dynamical Systems, 1.3 ed. Kindle Direct Publishing, 2026. [Online]. Available: https://fbullo.github.io/ctds

    work page 2026

  21. [21]

    Kron reduction of graphs with applications to electrical networks,

    F. D ¨orfler and F. Bullo, “Kron reduction of graphs with applications to electrical networks,”IEEE Transactions on Circuits and Systems I: Regular Papers, vol. 60, no. 1, pp. 150–163, 2012

    work page 2012

  22. [22]

    Robust distributed finite-time secondary frequency control of islanded ac microgrids with event- triggered mechanism,

    Z. Wang, H. Li, H. He, and Y . Sun, “Robust distributed finite-time secondary frequency control of islanded ac microgrids with event- triggered mechanism,” in2023 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT). IEEE, 2023, pp. 1–5

    work page 2023

  23. [23]

    Transient stability of power systems integrated with inverter-based generation,

    X. He and H. Geng, “Transient stability of power systems integrated with inverter-based generation,”IEEE Transactions on Power Systems, vol. 36, no. 1, pp. 553–556, 2021

    work page 2021

  24. [24]

    Analysis of power system stability under slowly varying perturbations in renewable generation,

    X. Ru, P. Yang, R. Liu, C. Jia, F. Duan, and F. Liu, “Analysis of power system stability under slowly varying perturbations in renewable generation,” in2024 IEEE 7th Student Conference on Electric Machines and Systems (SCEMS). IEEE, 2024, pp. 1–7

    work page 2024

  25. [25]

    Weak and semi- contraction for network systems and diffusively coupled oscillators,

    S. Jafarpour, P. Cisneros-Velarde, and F. Bullo, “Weak and semi- contraction for network systems and diffusively coupled oscillators,” IEEE Transactions on Automatic Control, vol. 67, no. 3, pp. 1285–1300, 2022

    work page 2022

  26. [26]

    R. A. Horn and C. R. Johnson,Matrix Analysis, 2nd ed. Cambridge University Press, 2012

    work page 2012

  27. [27]

    H. K. Khalil,Nonlinear Systems, 3rd ed. Upper Saddle River, NJ: Prentice Hall, 2002

    work page 2002

  28. [28]

    J. R. Munkres,Topology. Prentice Hall, 2000

    work page 2000