pith. sign in

arxiv: 2605.07152 · v1 · submitted 2026-05-08 · 🪐 quant-ph · cs.NA· cs.SY· eess.SY· math.NA· math.OC

Symplectic H2 Model Reduction for High-Dimensional Linear Quantum Systems

Pith reviewed 2026-05-11 00:50 UTC · model grok-4.3

classification 🪐 quant-ph cs.NAcs.SYeess.SYmath.NAmath.OC
keywords model reductionlinear quantum systemssymplectic methodsH2 normphysical realizabilityKrylov subspace methodsPetrov-Galerkin projectionoscillator chains
0
0 comments X

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.

The paper develops a symplectic Petrov-Galerkin framework that reduces the dimension of linear quantum systems while ensuring the reduced models automatically satisfy the physical realizability constraints of quantum mechanics, including canonical commutation relations and the input-output structure. Within this framework it introduces Quantum IRKA, which builds an enriched tangential rational Krylov subspace, extracts a symplectic basis via Gram-Schmidt, and updates interpolation points from mirror images of the reduced poles. Standard projection methods are ruled out because they break the required quantum identities; the symplectic approach avoids that breakage by construction. Numerical tests on oscillator chains and Kitaev-chain models demonstrate that the resulting models remain accurate at moderate cost and keep symplecticity to machine precision. A reader would care because this supplies a practical route to simulating or controlling large quantum devices with compact yet faithful models.

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

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

  • 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

Figures reproduced from arXiv: 2605.07152 by Alfo Borzi, Guofeng Zhang.

Figure 1
Figure 1. Figure 1: Hankel singular values; homogeneous (left) and heterogeneous (right). [PITH_FULL_IMAGE:figures/full_fig_p016_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Gramian spectra; homogeneous (left) and heterogeneous (right). [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Symplectic diagnostics; homogeneous (left) and heterogeneous (right). [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Shift convergence; homogeneous (left) and heterogeneous (right). [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Hankel singular values for the BKC benchmark; homogeneous (left) and [PITH_FULL_IMAGE:figures/full_fig_p020_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Gramian spectra for the BKC benchmark; homogeneous (left) and heteroge [PITH_FULL_IMAGE:figures/full_fig_p020_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Symplectic and left-inverse defects for the BKC benchmark; homogeneous [PITH_FULL_IMAGE:figures/full_fig_p021_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Reduced full-port PR residuals for the BKC benchmark; homogeneous (left) [PITH_FULL_IMAGE:figures/full_fig_p021_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Shift convergence for the BKC benchmark; homogeneous (left) and hetero [PITH_FULL_IMAGE:figures/full_fig_p022_9.png] view at source ↗
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.

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 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)
  1. [§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.
  2. [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)
  1. [§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.
  2. [§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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 2 axioms · 0 invented entities

The method rests on the existence of a symplectic basis for the state space of linear quantum systems and on the standard properties of the H2 norm for stable systems; no new free parameters or invented entities are introduced.

axioms (2)
  • domain assumption The system matrices satisfy the physical realizability identities that encode canonical commutation relations and quantum input-output structure.
    Invoked throughout the abstract as the constraint that standard projections violate.
  • standard math The underlying linear quantum system is stable and the H2 norm is well-defined.
    Required for the model reduction problem to make sense.

pith-pipeline@v0.9.0 · 5583 in / 1435 out tokens · 47333 ms · 2026-05-11T00:50:44.121249+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    Absil, R

    P.-A. Absil, R. Mahony, and R. Sepulchre,Optimization Algorithms on Matrix Manifolds, Princeton University Press, Princeton, NJ (2008)

  2. [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. [3]

    A. C. Antoulas,Approximation of Large-Scale Dynamical Systems, SIAM, Philadel- phia, PA (2005)

  4. [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)

  5. [5]

    A. C. Antoulas, C. A. Beattie, and S. Gugercin,Interpolatory Methods for Model Reduction, SIAM, Philadelphia, PA (2020)

  6. [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. [7]

    V. P. Belavkin, Nondemolition principle of quantum measurement theory,Foun- dations of Physics, 24(5), 685–714 (1994). DOI:0.1007/BF02054669

  8. [8]

    Beattie and S

    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. [9]

    Benner, S

    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. [10]

    Bouten, R

    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. [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. [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. [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. [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. [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. [16]

    Grigoletto and F

    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. [17]

    Grigoletto, Y

    T. Grigoletto, Y. Tao, F. Ticozzi, and L. Viola, Exact model reduction for continuous-time open quantum dynamics,Quantum, 9, 1814 (2025)

  18. [18]

    Grigoletto and F

    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. [19]

    Grigoletto, C

    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. [20]

    Gugercin and A

    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. [21]

    Gugercin, A

    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. [22]

    van Handel and H

    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. [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. [24]

    Jonas and B

    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. [25]

    Liang and D

    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. [26]

    McDonald, T

    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. [27]

    Mlinari´ c, C

    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. [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. [29]

    Peng and K

    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. [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. [31]

    I. R. Petersen, Quantum linear systems theory,The Open Automation and Control Systems Journal, 8, 67–93 (2016). DOI:10.2174/1874444301608010067

  32. [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. [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. [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. [35]

    Techakesari and H

    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. [36]

    van der Schaft,L 2-Gain and Passivity Techniques in Nonlinear Control, 2nd ed., Springer, London (2000)

    A. van der Schaft,L 2-Gain and Passivity Techniques in Nonlinear Control, 2nd ed., Springer, London (2000)

  37. [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. [38]

    Yan and J

    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. [39]

    Zhang and Z

    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. [40]

    Zhang and J

    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. [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...