© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2403.05908 · Quantum

Low-Rank Variational Quantum Algorithm for the Dynamics of Open Quantum Systems

Sara Diana Leal dos Santos, Xinyu Song, Vincenzo Savona

Published Feb 4, 2025• 7 citations

Abstract

The simulation of many-body open quantum systems is key to solving numerous outstanding problems in physics, chemistry, material science, and in the development of quantum technologies. Near-term quantum computers may bring considerable advantage for the efficient simulation of their static and dynamical properties, thanks to hybrid quantum-classical variational algorithms to approximate the dynamics of the density matrix describing the quantum state in terms of an ensemble average. Here, a variational quantum algorithm is developed to simulate the real-time evolution of the density matrix governed by the Lindblad master equation, under the assumption that the quantum state has a bounded entropy along the dynamics, entailing a low-rank representation of its density matrix. The algorithm encodes each pure state of the statistical mixture as a parametrized quantum circuit, and the associated probabilities as additional variational parameters stored classically, thereby requiring a significantly lower number of qubits than algorithms where the full density matrix is encoded in the quantum memory. Two variational ansatze are proposed, and their effectiveness is assessed in the simulation of the dynamics of a 2D dissipative transverse field Ising model. The results underscore the algorithm's efficiency in simulating the dynamics of open quantum systems in the low-rank regime with limited quantum resources on a near-term quantum device.

DOI: 10.22331/q-2025-02-04-1620

Related papers

Variational Quantum Simulation of Lindblad Dynamics via Quantum State Diffusion

Jianming Luo, Kaihan Lin, Xing Bao Gao

Quantum simulation of dynamics in open quantum systems is crucial but poses a significant challenge due to the non-Hermitian nature leading to nonunitary evolution and the limited quantum resources on current quantum computers. Here we introduce a variational hybrid quantum-classical algorithm designed for simulating the time evolution governed by the Lindblad master equation. Our approach involves on a stochastic unveiling of the density matrix, transforming the Lindblad equation into a wave function-based quantum state diffusion (QSD) method with the aim of reducing qubit requirements. We then apply variational quantum simulation (VQS) to efficiently capture the nonunitary evolution in QSD. We demonstrate our QSD-VQS algorithm by investigating the quantum dynamics in a two-level system subjected to an amplitude damping channel and a four-level transverse field Ising model within a dissipative environment including time-independent and periodic Hamiltonian cases. The results reveal its promising utility with upcoming hardware in the near future.

The Journal of Physical Chemistry Letters•Mar 22, 2024•21 citations

Adaptive variational low-rank dynamics for open quantum systems

Luca Gravina, Vincenzo Savona

We introduce a model-independent method for the efficient simulation of low-entropy systems, whose dynamics can be accurately described with a limited number of states. Our method leverages the time-dependent variational principle to efficiently integrate the Lindblad master equation, dynamically identifying and modifying the low-rank basis over which we decompose the system's evolution. By dynamically adapting the dimension of this basis, and thus the rank of the density matrix, our method maintains optimal representation of the system state, offering a substantial computational advantage over existing adaptive low-rank schemes in terms of both computational time and memory requirements. We demonstrate the efficacy of our method through extensive benchmarks on a variety of model systems, with a particular emphasis on multiqubit bosonic codes, a promising candidate for fault-tolerant quantum hardware. Our results highlight the method's versatility and efficiency, making it applicable to a wide range of systems characterized by arbitrary degrees of entanglement and moderate entropy throughout their dynamics. We provide an implementation of the method as a Julia package, making it readily available to use. Published by the American Physical Society 2024

arXiv:2312.13676•Physical Review Research•Apr 19, 2024•11 citations

Variational Quantum Simulation of General Processes

Suguru Endo, Jinzhao Sun, Ying Li +2 more

Variational quantum algorithms have been proposed to solve static and dynamic problems of closed many-body quantum systems. Here we investigate variational quantum simulation of three general types of tasks-generalized time evolution with a non-Hermitian Hamiltonian, linear algebra problems, and open quantum system dynamics. The algorithm for generalized time evolution provides a unified framework for variational quantum simulation. In particular, we show its application in solving linear systems of equations and matrix-vector multiplications by converting these algebraic problems into generalized time evolution. Meanwhile, assuming a tensor product structure of the matrices, we also propose another variational approach for these two tasks by combining variational real and imaginary time evolution. Finally, we introduce variational quantum simulation for open system dynamics. We variationally implement the stochastic Schrödinger equation, which consists of dissipative evolution and stochastic jump processes. We numerically test the algorithm with a 6-qubit 2D transverse field Ising model under dissipation.

arXiv:1812.08778•Physical Review Letters•Jun 29, 2020•237 citations

Variational quantum algorithms for simulation of Lindblad dynamics

Tasneem Watad, Netanel H. Lindner

Abstract We introduce variational hybrid classical-quantum algorithms to simulate the Lindblad master equation and its adjoint for time-evolving Markovian open quantum systems and quantum observables. Our methods are based on a direct representation of density matrices and quantum observables as quantum superstates. We design and optimize low-depth variational quantum circuits that efficiently capture the unitary and non-unitary dynamics of the solutions. We benchmark and test the algorithms on different models and system sizes, showing their potential for utility with near-future hardware.

arXiv:2305.02815•Quantum Science and Technology•Dec 21, 2023•14 citations

Hybrid quantum variational algorithm for simulating open quantum systems with near-term devices

M. Mahdian, H Davoodi Yeganeh

Abstract Hybrid quantum-classical (HQC) algorithms make it possible to use near-term quantum devices supported by classical computational resources by useful control schemes. In this paper, we develop an HQC algorithm using an efficient variational optimization approach to simulate open system dynamics under the Noisy-Intermediate Scale Quantum computer. Using the time-dependent variational principle (TDVP) method and extending it to McLachlan TDVP for density matrix which involves minimization of Frobenius norm of the error, we apply the unitary quantum circuit to obtain the time evolution of the open quantum system in the Lindblad formalism. Finally, we illustrate the use of our methods with detailed examples which are in good agreement with analytical calculations.

arXiv:2008.05321•Journal of Physics A Mathematical and Theoretical•Aug 7, 2020•14 citations

Hybrid algorithm simulating non-equilibrium steady states of an open quantum system

Hongyi Zhou, Mao, Rui, Xiaoming Sun

Non-equilibrium steady states are a focal point of research in the study of open quantum systems. Previous variational algorithms for searching these steady states have suffered from resource-intensive implementations due to vectorization or purification of the system density matrix, requiring large qubit resources and long-range coupling. In this work, we present a novel variational quantum algorithm that efficiently searches for non-equilibrium steady states by simulating the operator-sum form of the Lindblad equation. By introducing the technique of random measurement, we are able to estimate the nonlinear cost function while reducing the required qubit resources by half compared to previous methods. Additionally, we prove the existence of the parameter shift rule in our variational algorithm, enabling efficient updates of circuit parameters using gradient-based classical algorithms. To demonstrate the performance of our algorithm, we conduct simulations for dissipative quantum transverse Ising and Heisenberg models, achieving highly accurate results. Our approach offers a promising solution for effectively addressing non-equilibrium steady state problems while overcoming computational limitations and implementation challenges.

arXiv:2309.06665•Sep 13, 2023•2 citations

The simulation of many-body open quantum systems is key to solving numerous outstanding problems in physics, chemistry, material science, and in the development of quantum technologies. Near-term quantum computers may bring considerable advantage for the efficient simulation of their static and dynamical properties, thanks to hybrid quantum-classical variational algorithms to approximate the dynamics of the density matrix describing the quantum state in terms of an ensemble average. Here, a variational quantum algorithm is developed to simulate the real-time evolution of the density matrix governed by the Lindblad master equation, under the assumption that the quantum state has a bounded entropy along the dynamics, entailing a low-rank representation of its density matrix. The algorithm encodes each pure state of the statistical mixture as a parametrized quantum circuit, and the associated probabilities as additional variational parameters stored classically, thereby requiring a significantly lower number of qubits than algorithms where the full density matrix is encoded in the quantum memory. Two variational ansatze are proposed, and their effectiveness is assessed in the simulation of the dynamics of a 2D dissipative transverse field Ising model. The results underscore the algorithm's efficiency in simulating the dynamics of open quantum systems in the low-rank regime with limited quantum resources on a near-term quantum device.

Quantum simulation of dynamics in open quantum systems is crucial but poses a significant challenge due to the non-Hermitian nature leading to nonunitary evolution and the limited quantum resources on current quantum computers. Here we introduce a variational hybrid quantum-classical algorithm designed for simulating the time evolution governed by the Lindblad master equation. Our approach involves on a stochastic unveiling of the density matrix, transforming the Lindblad equation into a wave function-based quantum state diffusion (QSD) method with the aim of reducing qubit requirements. We then apply variational quantum simulation (VQS) to efficiently capture the nonunitary evolution in QSD. We demonstrate our QSD-VQS algorithm by investigating the quantum dynamics in a two-level system subjected to an amplitude damping channel and a four-level transverse field Ising model within a dissipative environment including time-independent and periodic Hamiltonian cases. The results reveal its promising utility with upcoming hardware in the near future.

We introduce a model-independent method for the efficient simulation of low-entropy systems, whose dynamics can be accurately described with a limited number of states. Our method leverages the time-dependent variational principle to efficiently integrate the Lindblad master equation, dynamically identifying and modifying the low-rank basis over which we decompose the system's evolution. By dynamically adapting the dimension of this basis, and thus the rank of the density matrix, our method maintains optimal representation of the system state, offering a substantial computational advantage over existing adaptive low-rank schemes in terms of both computational time and memory requirements. We demonstrate the efficacy of our method through extensive benchmarks on a variety of model systems, with a particular emphasis on multiqubit bosonic codes, a promising candidate for fault-tolerant quantum hardware. Our results highlight the method's versatility and efficiency, making it applicable to a wide range of systems characterized by arbitrary degrees of entanglement and moderate entropy throughout their dynamics. We provide an implementation of the method as a Julia package, making it readily available to use. Published by the American Physical Society 2024

Variational quantum algorithms have been proposed to solve static and dynamic problems of closed many-body quantum systems. Here we investigate variational quantum simulation of three general types of tasks-generalized time evolution with a non-Hermitian Hamiltonian, linear algebra problems, and open quantum system dynamics. The algorithm for generalized time evolution provides a unified framework for variational quantum simulation. In particular, we show its application in solving linear systems of equations and matrix-vector multiplications by converting these algebraic problems into generalized time evolution. Meanwhile, assuming a tensor product structure of the matrices, we also propose another variational approach for these two tasks by combining variational real and imaginary time evolution. Finally, we introduce variational quantum simulation for open system dynamics. We variationally implement the stochastic Schrödinger equation, which consists of dissipative evolution and stochastic jump processes. We numerically test the algorithm with a 6-qubit 2D transverse field Ising model under dissipation.

Abstract We introduce variational hybrid classical-quantum algorithms to simulate the Lindblad master equation and its adjoint for time-evolving Markovian open quantum systems and quantum observables. Our methods are based on a direct representation of density matrices and quantum observables as quantum superstates. We design and optimize low-depth variational quantum circuits that efficiently capture the unitary and non-unitary dynamics of the solutions. We benchmark and test the algorithms on different models and system sizes, showing their potential for utility with near-future hardware.

Abstract Hybrid quantum-classical (HQC) algorithms make it possible to use near-term quantum devices supported by classical computational resources by useful control schemes. In this paper, we develop an HQC algorithm using an efficient variational optimization approach to simulate open system dynamics under the Noisy-Intermediate Scale Quantum computer. Using the time-dependent variational principle (TDVP) method and extending it to McLachlan TDVP for density matrix which involves minimization of Frobenius norm of the error, we apply the unitary quantum circuit to obtain the time evolution of the open quantum system in the Lindblad formalism. Finally, we illustrate the use of our methods with detailed examples which are in good agreement with analytical calculations.

Non-equilibrium steady states are a focal point of research in the study of open quantum systems. Previous variational algorithms for searching these steady states have suffered from resource-intensive implementations due to vectorization or purification of the system density matrix, requiring large qubit resources and long-range coupling. In this work, we present a novel variational quantum algorithm that efficiently searches for non-equilibrium steady states by simulating the operator-sum form of the Lindblad equation. By introducing the technique of random measurement, we are able to estimate the nonlinear cost function while reducing the required qubit resources by half compared to previous methods. Additionally, we prove the existence of the parameter shift rule in our variational algorithm, enabling efficient updates of circuit parameters using gradient-based classical algorithms. To demonstrate the performance of our algorithm, we conduct simulations for dissipative quantum transverse Ising and Heisenberg models, achieving highly accurate results. Our approach offers a promising solution for effectively addressing non-equilibrium steady state problems while overcoming computational limitations and implementation challenges.