Symplectic H2 Model Reduction for High-Dimensional Linear Quantum Systems
Pith reviewed 2026-05-11 00:50 UTC · model grok-4.3
The pith
Symplectic projection reduces high-dimensional quantum models while preserving physical realizability.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Within the symplectic Petrov-Galerkin framework, the Quantum IRKA algorithm generates reduced-order models by extracting a symplectic basis from shifted linear solves and updating interpolation points from mirror images of the reduced poles, with all reduced matrices obtained by structure-preserving projection, thereby ensuring that physical realizability identities hold exactly.
What carries the argument
The symplectic Petrov-Galerkin projection paired with an enriched tangential rational Krylov subspace and symplectic Gram-Schmidt orthogonalization, which enforces the canonical symplectic constraint on the reduced trial space.
If this is right
- Reduced-order models preserve symplecticity and physical realizability to machine precision.
- Accurate H2 approximations are obtained for large-scale linear quantum systems with moderate computational cost.
- Reduction performance varies with dissipation geometry, channel placement, and system heterogeneity.
- Scalable H2 model reduction becomes feasible for linear quantum systems while strictly preserving physical structure.
Where Pith is reading between the lines
- This structure-preserving reduction could support faster simulation and controller design for quantum hardware by providing compact yet faithful models.
- Extensions to systems with time-varying parameters or feedback might follow from the same projection technique if the symplectic structure is maintained.
- Comparative benchmarks against non-symplectic methods on a wider range of quantum networks would clarify the trade-off between error and structure preservation.
Load-bearing premise
That the enriched tangential rational Krylov pool combined with symplectic basis extraction produces reduced models whose approximation error stays competitive while the physical realizability conditions hold exactly regardless of the dissipation geometry or channel arrangement.
What would settle it
A high-dimensional linear quantum system where Q-IRKA's H2 error is substantially larger than that of a standard method that ignores the structure, even though PR is preserved.
Figures
read the original abstract
The $\mathcal{H}_2$ model reduction problem for high-dimensional linear quantum systems is studied under the constraint of physical realizability (PR). This constraint requires preservation of the canonical commutation relations and the quantum input-output structure, and therefore prevents the direct use of standard projection methods. A symplectic Petrov-Galerkin framework is presented, in which reduced-order models automatically satisfy the PR identities by construction. Within this framework, a symplectic variant of the iterative rational Krylov algorithm is developed and referred to as Quantum IRKA (Q-IRKA). At each iteration, an enriched tangential rational Krylov pool is generated from shifted linear solves. A symplectic basis is then extracted by a Gram-Schmidt-type procedure, paired with symplectic conjugates, and normalized so that the reduced trial space satisfies the canonical symplectic constraint. The interpolation points are updated from selected mirror images of the poles of the current reduced-order model, while the reduced-order matrices are obtained exclusively by structure-preserving projection. Numerical experiments on low-channel oscillator-chain systems and on a bosonic Kitaev-chain-inspired benchmark show that Q-IRKA is effective for large-scale linear quantum systems. Symplecticity and PR are preserved to machine precision, and accurate reduced-order models are obtained with moderate computational cost. The results also show that reduction quality depends substantially on dissipation geometry, channel placement, heterogeneity, and reduced order. These findings indicate that scalable $\mathcal{H}_2$ model reduction of linear quantum systems can be achieved while strictly preserving the underlying physical structure.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the H2 model reduction problem for high-dimensional linear quantum systems subject to the physical realizability (PR) constraint, which enforces preservation of canonical commutation relations and the quantum input-output structure. It develops a symplectic Petrov-Galerkin projection framework in which reduced models satisfy the PR identities by construction. Within this framework, the authors introduce Quantum IRKA (Q-IRKA), a symplectic adaptation of the iterative rational Krylov algorithm that generates an enriched tangential rational Krylov pool from shifted linear solves, extracts a symplectic basis via Gram-Schmidt paired with conjugates, and updates interpolation points from mirror images of the reduced-model poles. Numerical experiments on low-channel oscillator chains and a bosonic Kitaev-chain benchmark are reported to show that symplecticity and PR are preserved to machine precision while producing accurate reduced-order models at moderate cost, with reduction quality depending on dissipation geometry, channel placement, heterogeneity, and reduced order.
Significance. If the reported accuracy and competitiveness hold under broader testing, the work would provide a useful structure-preserving reduction technique for large-scale linear quantum systems, enabling efficient simulation and control while automatically satisfying the underlying physical constraints. The by-construction enforcement of symplecticity and PR via the Petrov-Galerkin projection and symplectic Gram-Schmidt procedure is a clear methodological strength, as is the adaptation of IRKA with reported machine-precision adherence in the benchmarks. The explicit dependence of performance on system features such as dissipation geometry is a useful observation that could guide future applications.
major comments (2)
- [§5] §5 (Numerical Experiments): The abstract and results section claim that 'accurate reduced-order models are obtained' and that Q-IRKA yields competitive performance, yet no quantitative H2 error tables, baseline comparisons against classical IRKA or other projection methods, or explicit details on how the H2 error norm was computed and measured are provided. This absence makes it impossible to verify the central claim that the enriched tangential Krylov pool produces H2 errors competitive with non-structure-preserving methods for arbitrary dissipation geometries and channel placements.
- [Abstract and §5] Abstract and §5: Experiments are restricted to low-channel oscillator-chain systems and a single bosonic Kitaev-chain-inspired benchmark. Given the abstract's own statement that reduction quality 'depends substantially on dissipation geometry, channel placement, heterogeneity, and reduced order,' additional numerical tests on multi-channel or heterogeneous systems are required to substantiate the claim of effectiveness for general high-dimensional linear quantum systems; the current limited scope leaves the weakest assumption (H2 competitiveness for arbitrary cases) untested.
minor comments (2)
- [§3] Clarify the precise construction of the 'enriched tangential rational Krylov pool' and the selection criterion for mirror-image poles in the Q-IRKA iteration; a short pseudocode or explicit algorithmic listing would improve reproducibility.
- [§3] The manuscript should include a brief discussion of computational complexity for the shifted linear solves and symplectic Gram-Schmidt step when scaling to very high dimensions.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed feedback on our manuscript. The comments highlight important aspects for improving the clarity and substantiation of our claims regarding Q-IRKA. We address each major comment below and commit to revisions that strengthen the numerical evidence without altering the core methodological contributions.
read point-by-point responses
-
Referee: [§5] §5 (Numerical Experiments): The abstract and results section claim that 'accurate reduced-order models are obtained' and that Q-IRKA yields competitive performance, yet no quantitative H2 error tables, baseline comparisons against classical IRKA or other projection methods, or explicit details on how the H2 error norm was computed and measured are provided. This absence makes it impossible to verify the central claim that the enriched tangential Krylov pool produces H2 errors competitive with non-structure-preserving methods for arbitrary dissipation geometries and channel placements.
Authors: We agree that the current presentation would benefit from explicit quantitative data. In the revised manuscript, we will add tables in §5 reporting the computed H2 error norms for each benchmark, along with a precise description of the computation method (via the error system's controllability Gramian satisfying the Lyapunov equation, evaluated using the standard trace formula adapted to the quantum setting). We will also include comparisons against other structure-preserving methods such as symplectic balanced truncation to demonstrate the benefits of the enriched tangential Krylov subspace. Direct comparisons with classical IRKA are inherently limited because it does not enforce physical realizability, often producing non-physical reduced models; we will explicitly discuss this distinction while showing that Q-IRKA achieves competitive accuracy under the PR constraint. revision: yes
-
Referee: [Abstract and §5] Abstract and §5: Experiments are restricted to low-channel oscillator-chain systems and a single bosonic Kitaev-chain-inspired benchmark. Given the abstract's own statement that reduction quality 'depends substantially on dissipation geometry, channel placement, heterogeneity, and reduced order,' additional numerical tests on multi-channel or heterogeneous systems are required to substantiate the claim of effectiveness for general high-dimensional linear quantum systems; the current limited scope leaves the weakest assumption (H2 competitiveness for arbitrary cases) untested.
Authors: We acknowledge that the numerical experiments are currently focused on specific low-channel and Kitaev-inspired cases, which limits the breadth of evidence for general applicability. In the revised version, we will expand §5 with additional numerical tests on multi-channel oscillator chains and heterogeneous systems with varied dissipation geometries and channel placements. These will include quantitative H2 error results and explicit discussion of performance dependence on the factors noted in the abstract, thereby providing stronger substantiation for the method's effectiveness across a wider range of high-dimensional linear quantum systems. revision: yes
Circularity Check
No significant circularity; Q-IRKA derivation is self-contained via standard symplectic projection and IRKA adaptation.
full rationale
The paper constructs a symplectic Petrov-Galerkin projection that enforces physical realizability (PR) and canonical commutation relations directly from the projection matrices and Gram-Schmidt extraction, without defining the reduced model in terms of its own outputs. The interpolation-point update follows the standard IRKA rule of selecting mirror-image poles from the current reduced-order model, which is an iterative procedure independent of the target H2 error metric. No self-citations are load-bearing for the central claims, no ansatz is smuggled, and no fitted parameters are relabeled as predictions. The numerical effectiveness for H2 error is presented as an empirical outcome on specific benchmarks rather than a consequence derived by construction from the inputs.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The system matrices satisfy the physical realizability identities that encode canonical commutation relations and quantum input-output structure.
- standard math The underlying linear quantum system is stable and the H2 norm is well-defined.
Reference graph
Works this paper leans on
- [1]
-
[2]
N. H. Amini, P. Mason, and I. Ramadan, Feedback stabilization via a quantum projection filter,SIAM Journal on Control and Optimization, 63(1), S128–S147 (2025). DOI:10.1137/23M1591037. 23
-
[3]
A. C. Antoulas,Approximation of Large-Scale Dynamical Systems, SIAM, Philadel- phia, PA (2005)
work page 2005
-
[4]
A. C. Antoulas, C. A. Beattie, and S. Gugercin, Interpolatory model reduction of large-scale dynamical systems, in: J. Mohammadpour and K. M. Grigoriadis (Eds.),Efficient Modeling and Control of Large-Scale Systems, Springer, Boston, MA, pp. 3–58 (2010)
work page 2010
-
[5]
A. C. Antoulas, C. A. Beattie, and S. Gugercin,Interpolatory Methods for Model Reduction, SIAM, Philadelphia, PA (2020)
work page 2020
-
[6]
L. Bao, B. Qi, F. Nori, and D. Dong, Exponential sensitivity revival of noisy non- Hermitian quantum sensing with two-photon drives,Physical Review Research, 6, 023216 (2024). DOI:10.1103/PhysRevResearch.6.023216
-
[7]
V. P. Belavkin, Nondemolition principle of quantum measurement theory,Foun- dations of Physics, 24(5), 685–714 (1994). DOI:0.1007/BF02054669
work page 1994
-
[8]
C. Beattie and S. Gugercin, Structure-preserving model reduction for nonlinear port-Hamiltonian systems, inProceedings of the 50th IEEE Conference on Decision and Control and European Control Conference, Orlando, FL, USA, pp. 6564–6569 (2011). DOI:10.1109/CDC.2011.6161504
-
[9]
P. Benner, S. Gugercin, and K. Willcox, A survey of projection-based model re- duction methods for parametric dynamical systems,SIAM Review, 57(4), 483–531 (2015). DOI:10.1137/130932715
-
[10]
L. Bouten, R. van Handel, and M. R. James, An introduction to quantum fil- tering,SIAM Journal on Control and Optimization, 46(6), 2199–2241 (2007). DOI:10.1137/060651239
-
[11]
Z. Dong, G. Zhang, H.-W. J. Lee, and I. R. Petersen, Linear quantum systems: poles, zeros, invertibility and sensitivity,IEEE Transactions on Automatic Control, to appear, (2026). DOI:10.1109/TAC.2026.3680121
-
[12]
Q. Gao, G. Zhang, and I. R. Petersen, An exponential quantum pro- jection filter for open quantum systems,Automatica, 99, 59–68 (2019). DOI:10.1016/j.automatica.2018.10.014
-
[13]
Q. Gao, G. Zhang, and I. R. Petersen, An improved quantum projection filter, Automatica, 112, 108716 (2020). DOI:10.1016/j.automatica.2019.108716
-
[14]
Q. Gao, D. Dong, I. R. Petersen, and S. X. Ding, Design of a quantum projec- tion filter,IEEE Transactions on Automatic Control, 65(8), 3693–3700 (2020). DOI:10.1109/TAC.2019.2953457. 24
-
[15]
J. E. Gough and G. Zhang, On realization theory of quantum linear systems, Automatica, 59, 139–151 (2015). DOI:10.1016/j.automatica.2015.06.023
-
[16]
T. Grigoletto and F. Ticozzi, Model reduction for quantum systems: discrete-time quantum walks and open Markov dynamics,IEEE Transactions on Information Theory, 71, 8524 –8541 (2025). DOI:10.1109/TIT.2025.3601118
-
[17]
T. Grigoletto, Y. Tao, F. Ticozzi, and L. Viola, Exact model reduction for continuous-time open quantum dynamics,Quantum, 9, 1814 (2025)
work page 2025
-
[18]
T. Grigoletto and F. Ticozzi, Exact model reduction for discrete-time condi- tional quantum dynamics,IEEE Control Systems Letters, 8, 550–555 (2024). DOI:10.1109/LCSYS.2024.3399100
-
[19]
T. Grigoletto, C. Pellegrini, and F. Ticozzi, Quantum model reduc- tion for continuous-time quantum filters,Annales Henri Poincar´ e, (2025). DOI:10.1007/s00023-025-01622-7
-
[20]
S. Gugercin and A. C. Antoulas, A survey of model reduction by balanced trunca- tion and some new results,International Journal of Control, 77(8), 748–766 (2004). DOI:10.1080/00207170410001713448
-
[21]
S. Gugercin, A. C. Antoulas, and C. Beattie,H 2 model reduction for large-scale linear dynamical systems,SIAM Journal on Matrix Analysis and Applications, 30(2), 609–638 (2008). DOI:10.1137/060666123
-
[22]
R. van Handel and H. Mabuchi, Quantum projection filter for a highly nonlinear model in cavity QED,Journal of Optics B: Quantum and Semiclassical Optics, 7, S226–S236 (2005). DOI:10.1088/1464-4266/7/10/005
-
[23]
M. R. James, H. I. Nurdin, and I. R. Petersen,H ∞ control of linear quantum stochastic systems,IEEE Transactions on Automatic Control, 53(8), 1787–1803 (2008). DOI:10.1109/TAC.2008.929378
-
[24]
J. Jonas and B. Bamieh, An iterative tangential interpolation algorithm for model reduction of MIMO systems,arXiv preprint, arXiv:2506.03410 (2025). DOI:10.48550/arXiv.2506.03410
-
[25]
W. Liang and D. Dong, Stabilization of time-varying perturbed quan- tum systems via reduced filters,arXiv preprint, arXiv:2511.07949 (2025). DOI:10.48550/arXiv.2511.07949
-
[26]
A. McDonald, T. Pereg-Barnea, and A.A. Clerk, Phase-dependent chiral trans- port and effective non-Hermitian dynamics in a bosonic Kitaev–Majorana chain, Physical Review X, 8, 041031 (2018). DOI:10.1103/PhysRevX.8.041031. 25
-
[27]
P. Mlinari´ c, C. A. Beattie, Z. Drmaˇ c, and S. Gugercin, IRKA is a Riemannian gradient descent method,IEEE Transactions on Automatic Control, 70(5), 2979– 2991 (2025). DOI:10.1109/TAC.2024.3489416
-
[28]
H. I. Nurdin, Structures and transformations for model reduction of linear quan- tum systems,IEEE Transactions on Automatic Control, 59(9), 2413–2425 (2014). DOI:10.1109/TAC.2014.2322731
-
[29]
L. Peng and K. Mohseni, Symplectic model reduction of Hamiltonian systems,SIAM Journal on Scientific Computing, 38(1), A1–A27 (2016). DOI:10.1137/140978922
-
[30]
I. R. Petersen, Singular perturbation approximations for a class of linear quan- tum systems,IEEE Transactions on Automatic Control, 58(1), 193–198 (2013). DOI:10.1109/TAC.2012.2203030
-
[31]
I. R. Petersen, Quantum linear systems theory,The Open Automation and Control Systems Journal, 8, 67–93 (2016). DOI:10.2174/1874444301608010067
-
[32]
R. V. Polyuga and A. van der Schaft, Structure preserving model reduction of port- Hamiltonian systems by moment matching at infinity,Automatica, 46(4), 665–672 (2010). DOI:10.1016/j.automatica.2010.01.018
-
[33]
In: 2025 IEEE 64th Conference on Decision and Control (CDC)
I. Ramadan, N. H. Amini, and P. Mason, Quantum projection filtering for state estimation and feedback stabilization in open quantum systems, inProceedings of the 2025 IEEE 64th Conference on Decision and Control (CDC), pp. 7629–7634 (2025). DOI:10.1109/CDC57313.2025.11312262
-
[34]
J. J. Slim, C. C. Wanjura, M. Brunelli, J. del Pino, A. Nunnenkamp, and E. Verhagen, Optomechanical realization of the bosonic Kitaev chain,Nature, 627, 767–771 (2024). DOI:10.1038/s41586-024-07174-w
-
[35]
O. Techakesari and H. I. Nurdin, Tangential interpolatory projection for model reduction of linear quantum stochastic systems,IEEE Transactions on Automatic Control, 62(1), 5–19 (2017). DOI:10.1109/TAC.2016.2535501
-
[36]
A. van der Schaft,L 2-Gain and Passivity Techniques in Nonlinear Control, 2nd ed., Springer, London (2000)
work page 2000
-
[37]
G. Wu, S. Xue, G. Zhang, and I. R. Petersen,H 2 model reduction for lin- ear quantum systems,IEEE Transactions on Automatic Control, to appear, arXiv:2411.07603 (2024). DOI:10.1109/TAC.2026.3689010
-
[38]
W.-Y. Yan and J. Lam, An approximate approach toH 2 optimal model re- duction,IEEE Transactions on Automatic Control, 44(7), 1341–1358 (1999). DOI:10.1109/9.774107. 26
-
[39]
G. Zhang and Z. Dong, Linear quantum systems: a tutorial,Annual Reviews in Control, 54, 274–294 (2022). DOI:10.1016/j.arcontrol.2022.04.013
-
[40]
L. Zhang and J. Lam, OnH 2 model reduction of bilinear systems,Automatica, 38(2), 205–216 (2002). DOI:10.1016/S0005-1098(01)00204-7
-
[41]
K. Zhou, J. C. Doyle, and K. Glover,Robust and Optimal Control, Prentice Hall, Englewood Cliffs, NJ (1996). 27 Algorithm 1Symplectic Quantum IRKA (Q-IRKA) Require:Full-order PR model (A, B, C, D) satisfying (5); structure matricesJ n, Jm; reduced dimension 2r; initial shifts{σ (0) i }r i=1 withσ (0) i /∈spec(A); enrichment pa- rameterL per shift ; toleran...
work page 1996
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.