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:1812.11626 · Physical review. E

Unfolding a quantum master equation into a system of real-valued equations: Computationally effective expansion over the basis of SU(N) generators

Alexey Liniov, Iosif B. Meyerov, Козинов Евгений Александрович, Valentin Volokitin, Igor Yusipov, Mikhail Ivanchenko, S. Denisov

Published Nov 11, 2019• 9 citations
arXivHTMLPDF

Abstract

Dynamics of an open N-state quantum system is often modeled with a Markovian master equation describing the evolution of the system density operator. By using generators of SU(N) group as a basis, the density operator can be transformed into a real-valued "coherence-vector." A generator of the dissipative evolution, so-called "Lindbladian," can be expanded over the same basis and recast in the form of a real matrix. Together, these expansions result is a nonhomogeneous system of N^{2}-1 real-valued linear ordinary differential equations. Now one can, e.g., implement standard high-performance algorithms to integrate the system of equations forward in time while being sure in exact preservation of the trace (norm) and Hermiticity of the density operator. However, when performed in a straightforward way, the expansion turns to be an operation of the time complexity O(N^{10}). The complexity can be reduced when the number of dissipative operators is independent of N, which is often the case for physically meaningful models. Here we present an algorithm to transform quantum master equation into a system of real-valued differential equations and propagate it forward in time. By using a specific scalable model, we evaluate computational efficiency of the algorithm and demonstrate that it is possible to handle the model system with N=10^{3} states on a single node of a computer cluster.

DOI
10.1103/physreve.100.053305

Related papers

Transforming Lindblad Equations into Systems of Real-Valued Linear Equations: Performance Optimization and Parallelization of an Algorithm

Iosif B. Meyerov, Козинов Евгений Александрович, Alexey Liniov +4 more

With their constantly increasing peak performance and memory capacity, modern supercomputers offer new perspectives on numerical studies of open many-body quantum systems. These systems are often modeled by using Markovian quantum master equations describing the evolution of the system density operators. In this paper, we address master equations of the Lindblad form, which are a popular theoretical tools in quantum optics, cavity quantum electrodynamics, and optomechanics. By using the generalized Gell-Mann matrices as a basis, any Lindblad equation can be transformed into a system of ordinary differential equations with real coefficients. Recently, we presented an implementation of the transformation with the computational complexity, scaling as O(N5logN) for dense Lindbaldians and O(N3logN) for sparse ones. However, infeasible memory costs remains a serious obstacle on the way to large models. Here, we present a parallel cluster-based implementation of the algorithm and demonstrate that it allows us to integrate a sparse Lindbladian model of the dimension N=2000 and a dense random Lindbladian model of the dimension N=200 by using 25 nodes with 64 GB RAM per node.

arXiv:1912.01491•Entropy•Oct 6, 2020•2 citations

Transforming the Lindblad Equation into a System of Linear Equations: Performance Optimization and Parallelization.

Iosif B. Meyerov, Козинов Евгений Александрович, Alexey Liniov +4 more

Rapidly growing performance and memory capacity of modern supercomputers open new perspectives for numerical studies of open quantum systems. These systems are in contact with their environments and their modeling is typically based on Markovian kinetic equations, describing the evolution of the system density operators. Additional to the exponential growth of the system dimension $N$ with the number of the system's parts, the computational complexity scales quadratically with $N$, since we have to deal now with super-operators represented by $N^2 \times N^2$ matrices (instead of standard $N \times N$ matrices of operators). In this paper we consider the so-called Lindblad equation, a popular tool to model dynamics of open systems in quantum optics and superconducting quantum physics. Using the generalized Gell-Mann matrices as a basis, we transform the original Lindblad equation into a system of ordinary differential equations (ODEs) with real coefficients. Earlier, we proposed an implementation of this idea with the complexity of computations scaling as $O(N^5 log(N))$ for dense systems and $O(N^3 log(N))$ for sparse systems. However, infeasible memory costs remained an obstacle on the way to large models. Here we present a new parallel algorithm based on cluster manipulations with large amount of data needed to transform the original Lindblad equation into an ODE system. We demonstrate that the algorithm allows us to integrate sparse systems of dimension $N=2000$ and dense systems of dimension $N=200$ by using up to 25 nodes with 64 GB RAM per node. We also managed to perform large-scale supercomputer sampling to study the spectral properties of an ensemble of random Lindbladians for $N=200$, which is impossible without the distribution of memory costs among cluster nodes.

•Dec 3, 2019•2 citations

Algebraic structures of the Lindblad equation

Leonel Bixano, Guillermo López-Alvarez, Victor Alberto Cruz-Barriguete +4 more

We investigate the algebraic structure underlying the Lindblad equation for finite-dimensional open quantum systems. By introducing a suitable operator representation of the Liouville superoperator, we show that the dynamics can be formulated in terms of a closed algebra of Hermitian operators that is independent of the particular physical model. This formulation reveals that dissipative dynamics requires a substantially richer algebraic structure than purely unitary evolution, thereby providing a clear characterization of the additional complexity introduced by the Lindbladian. The resulting framework naturally leads to parametrizations of the dynamical map and to differential equations governing its evolution. We further derive recursion relations that enable the efficient construction of the algebra for systems of increasing dimension. Because the algebraic basis is universal, while all model-dependent information enters through a single set of coefficients, the proposed approach significantly reduces the computational cost of constructing the Liouville superoperator compared with direct methods. To facilitate the implementation of the method, we provide a Mathematica notebook containing a one-qubit example that can be systematically extended to an arbitrary number of qubits. The proposed framework therefore provides both a general mathematical description of finite-dimensional Lindblad dynamics and a practical foundation for efficient analytical and numerical implementations.

quant-pharXiv:2606.26477•New Journal of Physics•Jun 26, 2026•0 citations

Generic method for integrating Lindblad master equations

Jiayin Gu, Fan Zhang

The time evolution of Markovian open quantum systems is governed by Lindblad master equations, whose solution can be formally written as the Lindbladian exponential acting on the initial density matrix. By expanding this Lindbladian exponential into the Taylor series, we propose a generic method for integrating Lindblad master equations. In this method, the series is truncated, retaining a finite number of terms, and the iterative actions of Lindbladian on the density matrix follow the corresponding master equation. Our method offers significant improvements in numerical efficiencies both in memory cost and computation time, especially for systems with many degrees of freedom. Moreover, our proposed method can be integrated seamlessly with tensor networks. Two illustrative examples, a two-level system exhibiting damped Rabi oscillations and a driven dissipative Heisenberg chain, are used to demonstrate the validity of our method. The superiority of our method is benchmarked with detailed performance tests.

arXiv:2412.13661•Physical review. E•Sep 8, 2025•1 citation

Algebraic aspects of the driven dynamics in the density operator and correlation functions calculation for multi-level open quantum systems

Nikolai N. Bogolubov, Andrey Vladimirovich Soldatov

Exact and approximate master equations were derived by the projection operator method for the reduced statistical operator of a multi-level quantum system with finite number N of quantum eigenstates interacting with arbitrary external classical fields and dissipative environment simultaneously. It was shown that the structure of these equations can be simplified significantly if the free Hamiltonian driven dynamics of an arbitrary quantum multi-level system under the influence of the external driving fields as well as its Markovian and non-Markovian evolution, stipulated by the interaction with the environment, are described in terms of the SU(N) algebra representation. As a consequence, efficient numerical methods can be developed and employed to analyze these master equations for real problems in various fields of theoretical and applied physics. It was also shown that literally the same master equations hold not only for the reduced density operator but also for arbitrary nonequilibrium multi-time correlation functions as well under the only assumption that the system and the environment are uncorrelated at some initial moment of time. A calculational scheme was proposed to account for these lost correlations in a regular perturbative way, thus providing additional computable terms to the correspondent master equations for the correlation functions.

International Journal of Modern Physics B•Oct 13, 2017•0 citations

Markovian evolution from a novel scheme

Fardin Kheirandish, Zahra Iranshahi, Fatemeh Bakhtiari +1 more

The Markovian dynamics of open quantum many-body systems are typically governed by the Lindblad master equation, yet obtaining the reduced density matrix for bosonic and fermionic systems remains a formidable numerical challenge due to the exponential growth of the Hilbert space. Here, we introduce a computationally efficient framework for constructing the reduced density matrix by modeling the environment as a copy of the primary system with a monotonically decaying coupling that enforces unidirectional energy flow. Our method directly yields exact solutions that rigorously satisfy the Lindblad master equation for both bosonic and fermionic cases, bypassing the need for costly Liouvillian diagonalization. We validate our approach through applications to paradigmatic models, demonstrating accurate reproduction of dissipative dynamics across a broad parameter regime. This work provides a straightforward and powerful tool for simulating Markovian open-system evolution, with immediate applicability to quantum transport and control problems in many-body physics.

quant-pharXiv:2607.13809•Jul 16, 2026•0 citations
Dynamics of an open N-state quantum system is often modeled with a Markovian master equation describing the evolution of the system density operator. By using generators of SU(N) group as a basis, the density operator can be transformed into a real-valued "coherence-vector." A generator of the dissipative evolution, so-called "Lindbladian," can be expanded over the same basis and recast in the form of a real matrix. Together, these expansions result is a nonhomogeneous system of N^{2}-1 real-valued linear ordinary differential equations. Now one can, e.g., implement standard high-performance algorithms to integrate the system of equations forward in time while being sure in exact preservation of the trace (norm) and Hermiticity of the density operator. However, when performed in a straightforward way, the expansion turns to be an operation of the time complexity O(N^{10}). The complexity can be reduced when the number of dissipative operators is independent of N, which is often the case for physically meaningful models. Here we present an algorithm to transform quantum master equation into a system of real-valued differential equations and propagate it forward in time. By using a specific scalable model, we evaluate computational efficiency of the algorithm and demonstrate that it is possible to handle the model system with N=10^{3} states on a single node of a computer cluster.
With their constantly increasing peak performance and memory capacity, modern supercomputers offer new perspectives on numerical studies of open many-body quantum systems. These systems are often modeled by using Markovian quantum master equations describing the evolution of the system density operators. In this paper, we address master equations of the Lindblad form, which are a popular theoretical tools in quantum optics, cavity quantum electrodynamics, and optomechanics. By using the generalized Gell-Mann matrices as a basis, any Lindblad equation can be transformed into a system of ordinary differential equations with real coefficients. Recently, we presented an implementation of the transformation with the computational complexity, scaling as O(N5logN) for dense Lindbaldians and O(N3logN) for sparse ones. However, infeasible memory costs remains a serious obstacle on the way to large models. Here, we present a parallel cluster-based implementation of the algorithm and demonstrate that it allows us to integrate a sparse Lindbladian model of the dimension N=2000 and a dense random Lindbladian model of the dimension N=200 by using 25 nodes with 64 GB RAM per node.
Rapidly growing performance and memory capacity of modern supercomputers open new perspectives for numerical studies of open quantum systems. These systems are in contact with their environments and their modeling is typically based on Markovian kinetic equations, describing the evolution of the system density operators. Additional to the exponential growth of the system dimension NNN with the number of the system's parts, the computational complexity scales quadratically with NNN, since we have to deal now with super-operators represented by N2×N2N^2 \times N^2N2×N2 matrices (instead of standard N×NN \times NN×N matrices of operators). In this paper we consider the so-called Lindblad equation, a popular tool to model dynamics of open systems in quantum optics and superconducting quantum physics. Using the generalized Gell-Mann matrices as a basis, we transform the original Lindblad equation into a system of ordinary differential equations (ODEs) with real coefficients. Earlier, we proposed an implementation of this idea with the complexity of computations scaling as O(N5log(N))O(N^5 log(N))O(N5log(N)) for dense systems and O(N3log(N))O(N^3 log(N))O(N3log(N)) for sparse systems. However, infeasible memory costs remained an obstacle on the way to large models. Here we present a new parallel algorithm based on cluster manipulations with large amount of data needed to transform the original Lindblad equation into an ODE system. We demonstrate that the algorithm allows us to integrate sparse systems of dimension N=2000N=2000N=2000 and dense systems of dimension N=200N=200N=200 by using up to 25 nodes with 64 GB RAM per node. We also managed to perform large-scale supercomputer sampling to study the spectral properties of an ensemble of random Lindbladians for N=200N=200N=200, which is impossible without the distribution of memory costs among cluster nodes.
We investigate the algebraic structure underlying the Lindblad equation for finite-dimensional open quantum systems. By introducing a suitable operator representation of the Liouville superoperator, we show that the dynamics can be formulated in terms of a closed algebra of Hermitian operators that is independent of the particular physical model. This formulation reveals that dissipative dynamics requires a substantially richer algebraic structure than purely unitary evolution, thereby providing a clear characterization of the additional complexity introduced by the Lindbladian. The resulting framework naturally leads to parametrizations of the dynamical map and to differential equations governing its evolution. We further derive recursion relations that enable the efficient construction of the algebra for systems of increasing dimension. Because the algebraic basis is universal, while all model-dependent information enters through a single set of coefficients, the proposed approach significantly reduces the computational cost of constructing the Liouville superoperator compared with direct methods. To facilitate the implementation of the method, we provide a Mathematica notebook containing a one-qubit example that can be systematically extended to an arbitrary number of qubits. The proposed framework therefore provides both a general mathematical description of finite-dimensional Lindblad dynamics and a practical foundation for efficient analytical and numerical implementations.
The time evolution of Markovian open quantum systems is governed by Lindblad master equations, whose solution can be formally written as the Lindbladian exponential acting on the initial density matrix. By expanding this Lindbladian exponential into the Taylor series, we propose a generic method for integrating Lindblad master equations. In this method, the series is truncated, retaining a finite number of terms, and the iterative actions of Lindbladian on the density matrix follow the corresponding master equation. Our method offers significant improvements in numerical efficiencies both in memory cost and computation time, especially for systems with many degrees of freedom. Moreover, our proposed method can be integrated seamlessly with tensor networks. Two illustrative examples, a two-level system exhibiting damped Rabi oscillations and a driven dissipative Heisenberg chain, are used to demonstrate the validity of our method. The superiority of our method is benchmarked with detailed performance tests.
Exact and approximate master equations were derived by the projection operator method for the reduced statistical operator of a multi-level quantum system with finite number N of quantum eigenstates interacting with arbitrary external classical fields and dissipative environment simultaneously. It was shown that the structure of these equations can be simplified significantly if the free Hamiltonian driven dynamics of an arbitrary quantum multi-level system under the influence of the external driving fields as well as its Markovian and non-Markovian evolution, stipulated by the interaction with the environment, are described in terms of the SU(N) algebra representation. As a consequence, efficient numerical methods can be developed and employed to analyze these master equations for real problems in various fields of theoretical and applied physics. It was also shown that literally the same master equations hold not only for the reduced density operator but also for arbitrary nonequilibrium multi-time correlation functions as well under the only assumption that the system and the environment are uncorrelated at some initial moment of time. A calculational scheme was proposed to account for these lost correlations in a regular perturbative way, thus providing additional computable terms to the correspondent master equations for the correlation functions.
The Markovian dynamics of open quantum many-body systems are typically governed by the Lindblad master equation, yet obtaining the reduced density matrix for bosonic and fermionic systems remains a formidable numerical challenge due to the exponential growth of the Hilbert space. Here, we introduce a computationally efficient framework for constructing the reduced density matrix by modeling the environment as a copy of the primary system with a monotonically decaying coupling that enforces unidirectional energy flow. Our method directly yields exact solutions that rigorously satisfy the Lindblad master equation for both bosonic and fermionic cases, bypassing the need for costly Liouvillian diagonalization. We validate our approach through applications to paradigmatic models, demonstrating accurate reproduction of dissipative dynamics across a broad parameter regime. This work provides a straightforward and powerful tool for simulating Markovian open-system evolution, with immediate applicability to quantum transport and control problems in many-body physics.