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

Abstract

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 tool 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. This allows us to use standard high-performance parallel algorithms to integrate the equations and thus to emulate open quantum dynamics in a computationally efficient way. Recently we presented an implementation of the transform with the computational complexity scaling as $O(N^{5}logN)$ for dense Lindbaldians and $O(N^{3}logN)$ for sparse ones. However, infeasible memory costs remain 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.

References in corpus (18)

• Quantum States and Phases in Driven Open Quantum Systems with Cold Atoms
• Characterizing Quantum Supremacy in Near-Term Devices
• Concepts and methods in the theory of open quantum systems
• The Bloch Vector for N-Level Systems
• Condensates beyond mean field theory: quantum backreaction as decoherence
• 0.5 Petabyte Simulation of a 45-Qubit Quantum Circuit
• Complex Spacing Ratios: A Signature of Dissipative Quantum Chaos
• Route to chaos in optomechanics
• Universal Signature from Integrability to Chaos in Dissipative Open Quantum Systems
• A Brief History of the GKLS Equation
• Interaction-induced impeding of decoherence and anomalous diffusion
• Universal spectra of random Lindblad operators
• Differences between mean-field dynamics and N-particle quantum dynamics as a signature of entanglement
• Hierarchy of relaxation timescales in local random Liouvillians
• Mean-field dynamics of a two-mode Bose-Einstein condensate subject to noise and dissipation
• Massively parallel implementation and approaches to simulate quantum dynamics using Krylov subspace techniques
• Open source Matrix Product States: Exact diagonalization and other entanglement-accurate methods revisited in quantum systems
• Unfolding quantum master equation into a system of real-valued equations: computationally effective expansion over the basis of $SU(N)$ generators

Cited by in corpus (1)

• Random generators of Markovian evolution: A quantum-classical transition by superdecoherence