Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1509.05534 · IEEE Transactions on Automatic Control

Tangential Interpolatory Projection for Model Reduction of Linear Quantum Stochastic Systems

Onvaree Techakesari, Hendra I. Nurdin

Published Feb 29, 2016• 8 citations
arXivHTMLPDF

Abstract

This paper presents a model reduction method for the class of linear quantum stochastic systems often encountered in quantum optics and related fields. The approach is proposed on the basis of an interpolatory projection ensuring that specific input-output responses of the original and the reduced-order systems are matched at multiple selected points (or frequencies). Importantly, the physical realizability property of the original quantum system imposed by the laws of quantum mechanics is preserved under our tangential interpolatory projection. Error bounds are established for the proposed model reduction method and an avenue to select interpolation points is proposed. A passivity preserving model reduction method is also presented. Examples of both active and passive systems are provided to illustrate the merits of our proposed approach.

DOI
10.1109/tac.2016.2535501

Related papers

Tangential interpolatory projection for model reduction of completely passive linear quantum stochastic systems

Onvaree Techakesari, Hendra I. Nurdin

This paper presents a tangential interpolatory projection model reduction method for the class of completely passive linear quantum stochastic systems, often encountered in quantum optics and related fields. We show that the important physical realizability property of the full-order system is preserved through the interpolatory projection approximation. Sufficient conditions that guarantee asymptotic stability and minimality of the reduced-order approximation are established. We also present an optimization based routine for selection of interpolation points and propose a heuristic approach to select tangent directions, the key elements of the interpolatory projection method for multiple-input multiple-output systems. An example is presented to illustrate the utility of our proposed approach.

IEEE Conference on Decision and Control (CDC)•Dec 1, 2015•0 citations

Structures and Transformations for Model Reduction of Linear Quantum Stochastic Systems

Hendra I. Nurdin

The purpose of this paper is to develop a model reduction theory for linear quantum stochastic systems that are commonly encountered in quantum optics and related fields, modeling devices such as optical cavities and optical parametric amplifiers, as well as quantum networks composed of such devices. It is shown that subsystem truncation preserves the physical realizability property of linear quantum stochastic systems, and that the property of complete passivity is invariant under subsystem truncation. However, generic linear quantum stochastic systems need not have a balanced realization under symplectic transformations. Therefore, alternative notions of balancing, including so-called quasi-balancing, are developed, and necessary and sufficient conditions are derived. A truncation error bound is derived for quasi-balanceable linear quantum stochastic systems and it is shown that all asymptotically stable completely passive linear quantum stochastic systems have a quasi-balanced realization. An example is provided to illustrate model reduction in the context of low-pass optical filtering of coherent light using a network of optical cavities.

arXiv:1308.6062•IEEE Transactions on Automatic Control•May 8, 2014•30 citations

On balanced realization of linear quantum stochastic systems and model reduction by quasi-balanced truncation

Hendra I. Nurdin

This paper derives conditions under which co-diagonalization of the controllability and observability Gramians can be performed in a physically consistent manner for the class of linear quantum stochastic systems that are common in the fields of quantum optics and related areas like quantum nano- and optomechanical systems. Such a class is essentially a quantum analogue of classical (non-quantum) linear stochastic systems. However, unlike their classical counterparts, the laws of quantum mechanics impose nonlinear equality constraints on the system matrices of linear quantum stochastic systems, the so-called physical realizability constraints. Therefore, a meaningful model reduction procedure must yield a physically realizable reduced model satisfying these constraints. Previous work has shown that model reduction by subsystem truncation preserves physical realizability, but did not propose a method for truncating a subsystem. One candidate approach for this truncation is a quantum adaptation of the well-known balanced truncation method for classical linear time-invariant systems. It is shown in this paper that balanced realization for linear quantum stochastic systems is only possible under a strongly restrictive necessary and sufficient condition. However, the paper also derives less restrictive necessary and sufficient conditions for other realizations with simultaneously diagonal controllability and observability Gramians, and introduces the notion of a quasi-balanced realization. An example of the application of the results is provided to demonstrate quasi-balanced truncation in the linear quantum stochastic setting.

Proceedings of the ... American Control Conference/Proceedings of the American Control Conference•Jun 1, 2013•5 citations

$\mathscr {H}_{2}$ Model Reduction for Linear Quantum Systems

Guangpu Wu, Shibei Xue, Guofeng Zhang +1 more

In this paper, an $\mathscr{H}_2$ norm-based model reduction method for linear quantum systems is presented, which can obtain a physically realizable model with a reduced order for closely approximating the original system. The model reduction problem is described as an optimization problem, whose objective is taken as an $\mathscr{H}_2$ norm of the difference between the transfer function of the original system and that of the reduced one. Different from classical model reduction problems, physical realizability conditions for guaranteeing that the reduced-order system is also a quantum system should be taken as nonlinear constraints in the optimization. To solve the optimization problem with such nonlinear constraints, we employ a matrix inequality approach to transform nonlinear inequality constraints into readily solvable linear matrix inequalities (LMIs) and nonlinear equality constraints, so that the optimization problem can be solved by a lifting variables approach. We emphasize that different from existing work, which only introduces a criterion to evaluate the performance after model reduction, we guide our method to obtain an optimal reduced model with respect to the $\mathscr{H}_2$ norm. In addition, the above approach for model reduction is extended to passive linear quantum systems. Finally, examples of active and passive linear quantum systems validate the efficacy of the proposed method.

arXiv:2411.07603•IEEE Transactions on Automatic Control•Apr 29, 2026•0 citations

Model reduction of linear quantum stochastic systems by subsystem truncation

Hendra I. Nurdin

The purpose of the paper is to develop new theoretical results on model reduction of linear quantum stochastic systems by truncation of subsystems of such systems. It is shown that a truncated system is guaranteed to be physically realizable (i.e., it represents the dynamics of a valid open Markov quantum system). In addition, it is also shown that the class of completely passive linear quantum stochastic systems is invariant under such subsystem truncation.

Australian Control Conference•Nov 1, 2012•2 citations

Symplectic H2 Model Reduction for High-Dimensional Linear Quantum Systems

Alfo Borzi, Guofeng Zhang

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.

quant-pharXiv:2605.07152•May 11, 2026
This paper presents a model reduction method for the class of linear quantum stochastic systems often encountered in quantum optics and related fields. The approach is proposed on the basis of an interpolatory projection ensuring that specific input-output responses of the original and the reduced-order systems are matched at multiple selected points (or frequencies). Importantly, the physical realizability property of the original quantum system imposed by the laws of quantum mechanics is preserved under our tangential interpolatory projection. Error bounds are established for the proposed model reduction method and an avenue to select interpolation points is proposed. A passivity preserving model reduction method is also presented. Examples of both active and passive systems are provided to illustrate the merits of our proposed approach.
This paper presents a tangential interpolatory projection model reduction method for the class of completely passive linear quantum stochastic systems, often encountered in quantum optics and related fields. We show that the important physical realizability property of the full-order system is preserved through the interpolatory projection approximation. Sufficient conditions that guarantee asymptotic stability and minimality of the reduced-order approximation are established. We also present an optimization based routine for selection of interpolation points and propose a heuristic approach to select tangent directions, the key elements of the interpolatory projection method for multiple-input multiple-output systems. An example is presented to illustrate the utility of our proposed approach.
The purpose of this paper is to develop a model reduction theory for linear quantum stochastic systems that are commonly encountered in quantum optics and related fields, modeling devices such as optical cavities and optical parametric amplifiers, as well as quantum networks composed of such devices. It is shown that subsystem truncation preserves the physical realizability property of linear quantum stochastic systems, and that the property of complete passivity is invariant under subsystem truncation. However, generic linear quantum stochastic systems need not have a balanced realization under symplectic transformations. Therefore, alternative notions of balancing, including so-called quasi-balancing, are developed, and necessary and sufficient conditions are derived. A truncation error bound is derived for quasi-balanceable linear quantum stochastic systems and it is shown that all asymptotically stable completely passive linear quantum stochastic systems have a quasi-balanced realization. An example is provided to illustrate model reduction in the context of low-pass optical filtering of coherent light using a network of optical cavities.
This paper derives conditions under which co-diagonalization of the controllability and observability Gramians can be performed in a physically consistent manner for the class of linear quantum stochastic systems that are common in the fields of quantum optics and related areas like quantum nano- and optomechanical systems. Such a class is essentially a quantum analogue of classical (non-quantum) linear stochastic systems. However, unlike their classical counterparts, the laws of quantum mechanics impose nonlinear equality constraints on the system matrices of linear quantum stochastic systems, the so-called physical realizability constraints. Therefore, a meaningful model reduction procedure must yield a physically realizable reduced model satisfying these constraints. Previous work has shown that model reduction by subsystem truncation preserves physical realizability, but did not propose a method for truncating a subsystem. One candidate approach for this truncation is a quantum adaptation of the well-known balanced truncation method for classical linear time-invariant systems. It is shown in this paper that balanced realization for linear quantum stochastic systems is only possible under a strongly restrictive necessary and sufficient condition. However, the paper also derives less restrictive necessary and sufficient conditions for other realizations with simultaneously diagonal controllability and observability Gramians, and introduces the notion of a quasi-balanced realization. An example of the application of the results is provided to demonstrate quasi-balanced truncation in the linear quantum stochastic setting.
H2\mathscr {H}_{2}H2​ Model Reduction for Linear Quantum Systems
In this paper, an H2\mathscr{H}_2H2​ norm-based model reduction method for linear quantum systems is presented, which can obtain a physically realizable model with a reduced order for closely approximating the original system. The model reduction problem is described as an optimization problem, whose objective is taken as an H2\mathscr{H}_2H2​ norm of the difference between the transfer function of the original system and that of the reduced one. Different from classical model reduction problems, physical realizability conditions for guaranteeing that the reduced-order system is also a quantum system should be taken as nonlinear constraints in the optimization. To solve the optimization problem with such nonlinear constraints, we employ a matrix inequality approach to transform nonlinear inequality constraints into readily solvable linear matrix inequalities (LMIs) and nonlinear equality constraints, so that the optimization problem can be solved by a lifting variables approach. We emphasize that different from existing work, which only introduces a criterion to evaluate the performance after model reduction, we guide our method to obtain an optimal reduced model with respect to the H2\mathscr{H}_2H2​ norm. In addition, the above approach for model reduction is extended to passive linear quantum systems. Finally, examples of active and passive linear quantum systems validate the efficacy of the proposed method.
The H2\mathcal{H}_2H2​ 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 H2\mathcal{H}_2H2​ model reduction of linear quantum systems can be achieved while strictly preserving the underlying physical structure.