Abstract \noindent This paper applies the Variational Quantum Linear Solver (VQLS) to the continuous-time Lyapunov equation
ATX+XA=−Q. Using the Kronecker product formulation, the matrix equation is transformed into the linear system
Hz=b, where
H=In⊗AT+AT⊗In and
z=vec(X). VQLS approximates the solution via parameterized quantum circuits optimized in a hybrid quantum-classical loop. Numerical experiments on four control system benchmarks with Kronecker dimensions up to
25×25 demonstrate residual norms as low as
10−9, with logarithmic qubit scaling (2-5 qubits) and shallow circuit depths (6-14 layers) suitable for near-term quantum devices. Compared to the HHL algorithm, VQLS offers reduced circuit complexity and qubit overhead, making it practical for noisy intermediate-scale quantum (NISQ) hardware. This paper demonstrates the potential of variational quantum algorithms for solving Lyapunov equations in control and dynamical systems applications.
The purpose of this paper is to investigate the coherent feedback
H∞ control problem for linear quantum systems. A key contribution is a simplified design methodology that guarantees closed-loop stability and a prescribed level of disturbance attenuation. It is shown that for general linear quantum systems, a physically realizable quantum controller can be obtained by solving at most four Lyapunov equations. In the passive case, a necessary and sufficient condition is provided in terms of two uncoupled pairs of Lyapunov equations. These results represent a significant simplification over the standard approach, which requires solving two coupled algebraic Riccati equations. The effectiveness of the proposed method is demonstrated through two typical quantum optical devices: an empty optical cavity and a degenerate parametric amplifier. These results provide a computationally efficient procedure for the robust and optimal control of quantum optical and optomechanical systems.
This paper investigates a quantum system described by the Schrödinger equation, utilizing the concept of the quantum Lyapunov function. The Lyapunov function is chosen based on the mean value of a virtual mechanical quantity, where different values of P, the mean value of the virtual mechanical quantity in the Lyapunov function, have an impact on the attractive domain of the quantum system. The selected primary optimization algorithms approximating matrix P are the particle swarm optimization (PSO) algorithm and the simulated annealing (SA) algorithm. This study examines the characteristics of the system’s attraction domain under these two distinct algorithms and establishes stability conditions for the nonlinear quantum system. We introduce a method to estimate the size of the attractive domain using the Lyapunov function approach, converting the attractive domain issue into an optimization challenge. Numerical simulations are conducted in various two-dimensional test systems and spin 1/2 particle systems.