REVIEW 3 major objections 5 minor 28 references
Quantum computing-based solver for interacting power grids
T0 review · 3 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read This paper proposes a quantum-classical hybrid solver that uses RVVQE to extract complex eigenvalues of a power grid's non-Hermitian admittance matrix, enabling resonance mode analysis with logarithmic memory.
desk verdict The core variational cost function cannot select eigenvectors of a non-Hermitian matrix; the 5-bus validation is unexplained. 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 load-bearing mechanism is the Real Variance-based Variational Quantum Eigensolver (RVVQE), a VQE-style algorithm for non-Hermitian operators. It splits the non-Hermitian admittance matrix into its Hermitian parts H and K, maps each to Pauli strings on N = ⌈log2 n⌉ qubits, and minimises the real cost C_var = Var(H) + Var(K). Zero variance is interpreted as meaning that the trial state is a simultaneous eigenvector of both H and K, so the complex eigenvalue is read out as ⟨H⟩ + i⟨K⟩. This variance-minimization step, together with the logarithmic qubit encoding, is what carries the argument from a non-Hermitian eigenproblem to a near-term quantum measurement.
What would settle it
Take a small non-normal admittance matrix (e.g., a 3×3 complex symmetric matrix) with known eigenvalues and non-commuting H and K. Run RVVQE with the same cost function: if the optimizer cannot drive C_var to zero and the recovered eigenvalues deviate from exact diagonalization beyond numerical tolerance, the central claim fails. A direct check is to compute the residual norm ||(Y − λ I)|ψ⟩|| for the obtained eigenpair on any test instance.
Extended reading notes
Core claim
The central claim is that the RVVQE algorithm reliably extracts the complex eigenvalues of the non-Hermitian admittance matrix that governs resonance mode analysis. For each frequency step, the algorithm writes Y_B,f = H_f + iK_f with Hermitian H_f and K_f, encodes these operators into N-qubit Pauli strings, and prepares a variational state |ψ(θ)⟩. The cost C_var(θ) = ⟨H_f²⟩ − ⟨H_f⟩² + ⟨K_f²⟩ − ⟨K_f⟩² is minimized; when it reaches zero the state is a simultaneous eigenstate of H_f and K_f, and the associated eigenvalue is λ = ⟨H_f⟩ + i⟨K_f⟩. Using a 3-qubit ansatz on a zero-padded 5×5 admittance matrix, the authors report that the critical resonance modal impedances obtained this way are in
Load-bearing premise
The entire method rests on the existence of a quantum state that is simultaneously an eigenvector of both Hermitian components H and K of the non-Hermitian admittance matrix; for a general non-normal matrix this simultaneous eigenvector does not exist, so the variance cost may never reach zero and the extracted 'eigenvalues' could be meaningless.
Editorial extensions
If this is right
- Resonance mode analysis of a 70,000-bus grid would need about 17 qubits, eliminating the classical O(n²) memory bottleneck for admittance matrices.
- The same RVVQE workflow applies directly to larger IEEE test systems (118-bus, 300-bus) and to synthetic continental-scale grids such as those proposed as future benchmarks.
- The reported agreement was obtained without quantum error mitigation; the paper argues that adding ZNE or PEC should further improve accuracy on physical hardware.
- The extracted modal impedances plug directly into stability criteria such as the positive-mode-damping condition, so the quantum output is usable for harmonic mitigation and contingency assessment.
Reading between the lines
- The zero-variance objective implicitly requires H and K to have a common eigenvector; for a general non-normal admittance matrix this is not guaranteed, so the method's validity may be restricted to cases where H and K nearly commute or to matrices that are unitarily diagonalizable. The paper does not address this constraint.
- The empirical validation uses a 5-bus system padded to 8×8, which is a very small non-normal matrix; scaling up may reveal that the variance landscape is rugged and that the classical optimizer struggles to reach global minima, a risk the paper does not quantify.
- The same variance-based technique could be adapted to other non-Hermitian eigenvalue problems arising in open quantum systems (Lindbladians) or linear stability analysis, so the contribution may be a general template rather than a grid-specific solution.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a hybrid quantum-classical framework for Resonance Mode Analysis (RMA) of power transmission grids. The non-Hermitian bus admittance matrix Y_B,f is decomposed into Hermitian parts H_f and K_f, and the Real Variance-based Variational Quantum Eigensolver (RVVQE) is used to minimize C_var = Var(H_f) + Var(K_f) as a proxy for finding eigenstates. The authors claim that this yields the complex eigenvalues of Y_B,f with a logarithmic number of qubits, and they validate the approach on a 5-bus system padded to 8×8, comparing the critical modal impedance |Z_cm| against exact classical diagonalization over 100 Hz to 3000 Hz.
Significance. The application of quantum eigensolvers to power-system resonance analysis is timely, and the authors are right that standard Hermitian VQE is inapplicable to non-Hermitian admittance matrices. The paper deserves credit for benchmarking against exact SciPy diagonalization rather than fitting to the target, and for clearly identifying the non-Hermitian challenge. However, the central algorithmic claim is not supported by the stated mathematics: the cost function in Eq. (7) vanishes only on simultaneous eigenstates of H_f and K_f, which generally do not coincide with the eigenstates of a non-normal Y_B,f. This is not a minor technical gap but a load-bearing flaw affecting Algorithm 1, the convergence criterion, and the interpretation of the reported numerical agreement. The scalability argument also omits the cost of Pauli decomposition and measurement. If the central claim were correct, the paper would be a valuable bridge between quantum algorithms and power-system stability, but in its present form the proposed method does not demonstrably compute the desired eigenvalues.
major comments (3)
- [§II-B, Eq. (7); Algorithm 1, steps 9–15] The deflation/orthogonality penalty described in §III is incompatible with non-Hermitian eigenproblems. Eigenvectors of a non-normal matrix are not mutually orthogonal in the Hilbert-space inner product ⟨·|·⟩. For a complex symmetric Y, the appropriate biorthogonality relation is v_i^T v_j = δ_ij after suitable scaling, not ⟨v_i|v_j⟩=0. Enforcing ⟨ψ_current|ψ_found⟩=0 as a penalty can move the optimizer away from true eigenvectors and generally prevents extraction of the full spectrum. Since the manuscript claims that the sequential deflation technique yields all eigenvalues, this is an independent correctness issue for Algorithm 1's stated output.
- [§III, 'To obtain all the eigenvalues...'] The deflation/orthogonality penalty described in §III is incompatible with non-Hermitian eigenproblems. Eigenvectors of a non-normal matrix are not mutually orthogonal in the Hilbert-space inner product ⟨·|·⟩. For a complex symmetric Y, the appropriate biorthogonality relation is v_i^T v_j = δ_ij after suitable scaling, not ⟨v_i|v_j⟩=0. Enforcing ⟨ψ_current|ψ_found⟩=0 as a penalty can move the optimizer away from true eigenvectors and generally prevents extraction of the full spectrum. Since the manuscript claims that the sequential deflation technique yields all eigenvalues, this is an independent correctness issue for Algorithm 1's stated output.
- [§II-B, Algorithm 1 step 7; §IV] The logarithmic-memory scalability claim is not established. The number of qubits N=⌈log2 n⌉ is logarithmic, but the Hermitian components H_f and K_f must be represented as sums of Pauli strings. For a dense n×n matrix this decomposition already has O(n²)=O(4^N) terms, and for sparse Y_B,f the number of Pauli terms and the associated measurement overhead are not analyzed. A logarithmic qubit count does not imply logarithmic total resources: the classical preprocessing, the number of circuit evaluations, and the memory needed to store the Pauli decomposition may reintroduce the scaling bottleneck the paper claims to bypass.
minor comments (5)
- [Fig. 4] Typo: 'Qauntum Simulation' should be 'Quantum Simulation'.
- [References [23]] The central RVVQE algorithm is delegated to an unpublished, same-author preprint (ref. [23]) with no independent proof or implementation details. Algorithm 1 is a useful summary, but important choices—ansatz design, optimizer, deflation-penalty strength, zero-eigenvalue removal threshold, diagonal-shift amount—are only described verbally, which hampers reproducibility.
- [§III, zero-padding] Padding the 5×5 Ybus to 8×8 introduces artificial zero eigenvalues. The manuscript says an 'automatic validation routine removes zero eigenvalues,' but it does not explain how these are distinguished from physically meaningful near-zero resonance modes without biasing the extracted critical modes.
- [§II-A, Eq. (4)] The identity R_f = L_f^{-1} = L_f^T is not automatically true for every complex symmetric matrix; it requires a particular scaling of the eigenvectors. Please state the normalization or justify the equality, since it underlies the participation-factor formula in Eq. (5).
- [Fig. 4 and §III] The validation shows only the magnitude of the critical modal impedance |Z_cm|. No comparison of the full eigenvalue spectra, eigenvectors, or participation factors is reported, and no convergence history or final C_var values are given. Such data would be necessary to assess whether the quantum simulation actually converged to the intended eigenstates.
Circularity Check
RVVQE's validity is inherited from the authors' own unpublished preprint; the numerical benchmark is external, but the algorithm's core is self-cited.
-
self citation load bearing
[Section II-B (Quantum Embedding and the RVVQE Algorithm); Algorithm 1, lines 9-15; Section III]
"To solve this, we apply a systematic formulation based on a Real Variance-based Variational Quantum Eigensolver (RVVQE) algorithm, which is specifically designed for non-Hermitian operators [23]."
The central premise—that RVVQE extracts complex eigenvalues of non-Hermitian matrices—is not derived in this paper. It is imported from reference [23], an unpublished preprint by the same three authors (D. Pandey, A. K. Das, P. Arumugam). The cost function (Eq. 7) and the convergence criterion Cvar≈0 are the entire selection mechanism, but no proof is offered that zero-variance states are eigenvectors of a general non-normal Y. The 5-bus SciPy comparison checks the output values but does not independently establish the algorithm's theoretical validity; that validity rests on the self-citation.
full rationale
I find no step where a 'prediction' is literally fitted to the target by construction. The resonance modal impedances in Fig. 4 are compared against exact SciPy diagonalization, so the benchmark is external and the reported numbers are not obviously manufactured from the reference values. The principal load-bearing weakness is a self-citation chain: the RVVQE method and its cost function are taken from the authors' own unpublished preprint [23], with no derivation or proof of the simultaneous-eigenvector condition needed for Cvar=0 to correspond to eigenvalues of Y. That is a genuine self-citation that carries the algorithm's validity. The additional mathematical concern—that for a non-normal complex-symmetric Y, the Hermitian parts H and K need not commute, so exact eigenvectors of Y need not be zeros of Cvar—is a correctness risk rather than a circularity: it challenges whether the stated algorithm can produce the plotted agreement, but it does not make the output equal to the input by definition. Accordingly, the circularity score is 3.
Assumptions & free parameters
free parameters (4)
- Deflation penalty coefficient
- Diagonal shift amount for closely spaced eigenvalues
- Zero-eigenvalue removal threshold
- Grid-search initialization step over [0,2π]
assumptions (4)
- ad hoc to paper There exists a state |ψ⟩ with Var(H)+Var(K)=0 for the non-Hermitian Y_B,f.
- domain assumption RVVQE (ref [23]) correctly solves non-Hermitian eigenvalue problems as used in Algorithm 1.
- domain assumption The admittance matrix Y_B,f is complex symmetric with R_f = L_f^{-1} = L_f^T.
- standard math The Hermitian components H and K can be mapped to Pauli strings for the N-qubit system.
Cite this review
Pith. "Pith review of Quantum computing-based solver for interacting power grids." pith.science (2026). https://pith.science/paper/UTGZEENI
@misc{pith2026260729582,
author = {Pith},
title = {Pith review of: Quantum computing-based solver for interacting power grids},
year = {2026},
howpublished = {\url{https://pith.science/paper/UTGZEENI}},
note = {Machine review of arXiv:2607.29582}
}
read the original abstract
The proliferation of power electronics in multi-terminal transmission grids has increasingly led to harmonic distortions and dynamic instabilities. While Resonance Mode Analysis (RMA) provides deep insights into these system resonances, evaluating the critical modes of large-scale grids presents a severe computational bottleneck. Classical iterative techniques must continuously diagonalize massively high-dimensional, non-Hermitian admittance matrices across a wide frequency spectrum, a process that rapidly exhausts classical memory and processing limits. To overcome this scaling barrier, we propose a novel quantum-classical hybrid methodology that natively maps the transmission grid's admittance matrix onto a Quantum Processing Unit (QPU). Because the grid's matrix is non-Hermitian, standard quantum eigensolvers are insufficient; thus, we employ the Real Variance-based Variational Quantum Eigensolver (RVVQE) algorithm to accurately extract the complex eigenvalues that represent the system's modes. Validated against a standard 5-bus transmission system, the quantum-derived critical-resonance modal impedances demonstrate near-perfect alignment with the exact classical frequency responses. Crucially, by encoding the grid's state logarithmically into quantum memory, this methodology bypasses classical RAM limitations. The successful implementation of the RVVQE framework not only bridges the mathematical topologies of dissipative electrical grids and open quantum systems but also provides a profoundly scalable architecture capable of diagnosing resonance instabilities in massive, continental-scale networks that currently exceed classical computational boundaries.
Figures
Reference graph
Works this paper leans on
-
[1]
Renewable energy transmission by hvdc across the continent: system challenges and opportunities,
J. Sun, M. Li, Z. Zhang, T. Xu, J. He, H. Wang, and G. Li, “Renewable energy transmission by hvdc across the continent: system challenges and opportunities,”CSEE Journal of Power and Energy Systems, vol. 3, no. 4, pp. 353–364, 2017. [Online]. Available: https://ieeexplore.ieee.org/document/8233581
arXiv 2017
-
[2]
Modal frequency sensitivity for power system harmonic resonance analysis,
Y . Cui and X. Wang, “Modal frequency sensitivity for power system harmonic resonance analysis,”IEEE Transactions on Power Delivery, vol. 27, no. 2, pp. 1010–1017, 2012. [Online]. Available: https://ieeexplore.ieee.org/document/6164298
arXiv 2012
-
[3]
Computational time efficiency analysis for resonance studies in transmission grids and microgrid clusters,
O. Cartiel, J.-J. Mesas, L. Monjo, and L. Sainz, “Computational time efficiency analysis for resonance studies in transmission grids and microgrid clusters,”Mathematics and Computers in Simulation, vol. 243, pp. 486–498, 2026. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S0378475425005300
2026
-
[4]
Harmonic resonance mode analysis,
W. Xu, Z. Huang, Y . Cui, and H. Wang, “Harmonic resonance mode analysis,”IEEE Transactions on Power Delivery, vol. 20, no. 2, pp. 1182–1190, 2005. [Online]. Available: https://ieeexplore.ieee.org/ document/1413368
arXiv 2005
-
[5]
Y . Zhan, X. Xie, H. Liu, H. Liu, and Y . Li, “Frequency-domain modal analysis of the oscillatory stability of power systems with high-penetration renewables,”IEEE Transactions on Sustainable Energy, vol. 10, no. 3, pp. 1534–1543, 2019. [Online]. Available: https://ieeexplore.ieee.org/document/8643991
arXiv 2019
-
[6]
Application of modal sensitivity for power system harmonic resonance analysis,
Z. Huang, Y . Cui, and W. Xu, “Application of modal sensitivity for power system harmonic resonance analysis,”IEEE Transactions on Power Systems, vol. 22, no. 1, pp. 222–231, 2007. [Online]. Available: https://ieeexplore.ieee.org/document/4077113
arXiv 2007
-
[7]
IEEE 14-bus system,
“IEEE 14-bus system,” Electric Grid Test Case Repository, ac- cessed: Mar. 2026. [Online]. Available: https://electricgrids.engr.tamu. edu/electric-grid-test-cases/ieee-14-bus-system/
2026
-
[8]
IEEE 24-bus system,
“IEEE 24-bus system,” Electric Grid Test Case Repository, ac- cessed: Mar. 2026. [Online]. Available: https://electricgrids.engr.tamu. edu/electric-grid-test-cases/ieee-24-bus-system/
2026
Show all 28 references
-
[9]
IEEE 118-bus system,
“IEEE 118-bus system,” Electric Grid Test Case Repository, ac- cessed: Mar. 2026. [Online]. Available: https://electricgrids.engr.tamu. edu/electric-grid-test-cases/ieee-118-bus-system/
2026
-
[10]
IEEE 300-bus system,
“IEEE 300-bus system,” Electric Grid Test Case Repository, ac- cessed: Mar. 2026. [Online]. Available: https://electricgrids.engr.tamu. edu/electric-grid-test-cases/ieee-300-bus-system/. 6
2026
-
[11]
Stability assessment for multi-infeed grid-connected vscs modeled in the admittance matrix form,
L. Orellana, L. Sainz, E. Prieto-Araujo, and O. Gomis-Bellmunt, “Stability assessment for multi-infeed grid-connected vscs modeled in the admittance matrix form,”IEEE Transactions on Circuits and Systems I: Regular Papers, vol. 68, no. 9, pp. 3758–3771, 2021. [Online]. Availab...
2021
-
[12]
Mechanism analysis of modal resonance in converter-based power systems,
X. Zhou and S. Bu, “Mechanism analysis of modal resonance in converter-based power systems,”IEEE Transactions on Power Systems, vol. 40, no. 4, pp. 3528–3542, 2025. [Online]. Available: https://ieeexplore.ieee.org/document/10833803
2025
-
[13]
Bus participation factor analysis for harmonic instability in power electronics based power systems,
E. Ebrahimzadeh, F. Blaabjerg, X. Wang, and C. L. Bak, “Bus participation factor analysis for harmonic instability in power electronics based power systems,”IEEE Transactions on Power Electronics, vol. 33, no. 12, pp. 10 341–10 351, 2018. [Online]. Available: https://ieeexplor...
2018
-
[14]
Saad,Numerical Methods for Large Eigenvalue Problems
Y . Saad,Numerical Methods for Large Eigenvalue Problems. Philadelphia, PA, USA: SIAM, 2011. [Online]. Available: https: //epubs.siam.org/doi/abs/10.1137/1.9781611970739.bm
2011 doi
-
[15]
G. H. Golub and C. F. Van Loan,Matrix computations (3rd ed.). USA: Johns Hopkins University Press, 1996
1996
-
[16]
A faster resonance mode analysis approach based on a modified shifted- inverse power iteration method,
O. Cartiel, J. Jos ´e Mesas, L. Sainz, and A. Fabregas, “A faster resonance mode analysis approach based on a modified shifted- inverse power iteration method,”IEEE Transactions on Power Delivery, vol. 38, no. 6, pp. 4145–4156, 2023. [Online]. Available: https://ieeexplore.iee...
2023
-
[17]
An efficient implementation of the nonsymmetric lanczos algorithm,
D. Day, “An efficient implementation of the nonsymmetric lanczos algorithm,”SIAM Journal on Matrix Analysis and Applications, vol. 18, no. 3, pp. 566–589, 1997. [Online]. Available: https: //doi.org/10.1137/S0895479895292503
1997 doi
-
[18]
ACTIVSg70k: 70,000 bus synthetic grid on footprint of East- ern United States,
“ACTIVSg70k: 70,000 bus synthetic grid on footprint of East- ern United States,” Electric Grid Test Case Repository, accessed: Mar. 2026. [Online]. Available: https://electricgrids.engr.tamu.edu/ electric-grid-test-cases/activsg70k/
2026
-
[19]
Simulating physics with computers,
R. P. Feynman, “Simulating physics with computers,”International Journal of Theoretical Physics, vol. 21, no. 6, pp. 467–488, 1982. [Online]. Available: https://doi.org/10.1007/BF02650179
1982 doi
-
[20]
Quantum random access memory,
V . Giovannetti, S. Lloyd, and L. Maccone, “Quantum random access memory,”Phys. Rev. Lett., vol. 100, p. 160501, Apr 2008. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevLett.100.160501
2008 doi
-
[21]
A variational eigenvalue solver on a photonic quantum processor,
A. Peruzzo, J. R. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O’brien, “A variational eigenvalue solver on a photonic quantum processor,”Nat. Commun., vol. 5, p. 4213,
-
[22]
Quantum simulation of nuclear hamiltonian with a generalized transformation for gray code encoding,
P. Siwach and P. Arumugam, “Quantum simulation of nuclear hamiltonian with a generalized transformation for gray code encoding,” Phys. Rev. C, vol. 104, p. 034301, Sep 2021. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevC.104.034301
2021 doi
-
[23]
Real variance-based variational quantum eigensolver for non-hermitian matrices,
D. Pandey, A. K. Das, and P. Arumugam, “Real variance-based variational quantum eigensolver for non-hermitian matrices,” 2026. [Online]. Available: https://arxiv.org/abs/2603.28892
2026
-
[24]
Quantum simulations of nuclear resonances with variational methods,
A. Singh, P. Siwach, and P. Arumugam, “Quantum simulations of nuclear resonances with variational methods,”Phys. Rev. C, vol. 112, p. 024323, Aug 2025. [Online]. Available: https://link.aps.org/doi/10. 1103/x8vx-y5mh
2025
-
[25]
A rapid modal analysis method for harmonic resonance using modified power iteration,
Z. Li, H. Hu, Y . Zhou, and Z. He, “A rapid modal analysis method for harmonic resonance using modified power iteration,”IEEE Transactions on Power Delivery, vol. 33, no. 3, pp. 1495–1497, 2018. [Online]. Available: https://ieeexplore.ieee.org/document/7553550
2018
-
[26]
Probabilistic harmonic resonance assessment considering power system uncertainties,
Z. Li, H. Hu, Y . Wang, L. Tang, Z. He, and S. Gao, “Probabilistic harmonic resonance assessment considering power system uncertainties,”IEEE Transactions on Power Delivery, vol. 33, no. 6, pp. 2989–2998, 2018. [Online]. Available: https://ieeexplore.ieee.org/document/8382332
2018
-
[27]
Siso impedance-based stability analysis for system-level small-signal stability assessment of large-scale power electronics-dominated power systems,
D. Yang and Y . Sun, “Siso impedance-based stability analysis for system-level small-signal stability assessment of large-scale power electronics-dominated power systems,”IEEE Transactions on Sustainable Energy, vol. 13, no. 1, pp. 537–550, 2022. [Online]. Available: https://i...
2022
-
[2014]
Available: https://www.nature.com/articles/ncomms5213
[Online]. Available: https://www.nature.com/articles/ncomms5213
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.