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

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.

[1]  E. Sudarshan,et al.  Completely Positive Dynamical Semigroups of N Level Systems , 1976 .

[2]  A. Edelman,et al.  Random matrix theory , 2005, Acta Numerica.

[3]  William Gropp,et al.  Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis , 2013, HiPC 2013.

[5]  D. Witthaut,et al.  Mean-field dynamics of a two-mode Bose–Einstein condensate subject to noise and dissipation , 2008, 0808.3679.

[6]  L. Pastur,et al.  Eigenvalue Distribution of Large Random Matrices , 2011 .

[7]  R. Xu,et al.  Theory of open quantum systems , 2002 .

[8]  J R Anglin,et al.  Bose-Einstein condensates beyond mean field theory: quantum backreaction as decoherence. , 2001, Physical review letters.

[9]  Desmond J. Higham,et al.  Numerical Methods for Ordinary Differential Equations - Initial Value Problems , 2010, Springer undergraduate mathematics series.

[10]  K. Lendi,et al.  Quantum Dynamical Semigroups and Applications , 1987 .

[11]  J. Lambert.  Numerical Methods for Ordinary Differential Equations , 1991 .

[12]  David J. Luitz,et al.  Hierarchy of Relaxation Timescales in Local Random Liouvillians. , 2020, Physical review letters.

[13]  Peter Schlagheck,et al.  Experimental observation of the Poincaré-Birkhoff scenario in a driven many-body quantum system , 2017 .

[14]  M. Gell-Mann.  Symmetries of baryons and mesons , 1962 .

[15]  G. Kimura.  The Bloch Vector for N-Level Systems , 2003, quant-ph/0301152.

[16]  Thomas Häner,et al.  0.5 Petabyte Simulation of a 45-Qubit Quantum Circuit , 2017, SC17: International Conference for High Performance Computing, Networking, Storage and Analysis.

[17]  Daniel Jaschke,et al.  Open source matrix product states: exact diagonalization and other entanglement-accurate methods revisited in quantum systems , 2018, Journal of Physics A: Mathematical and Theoretical.

[18]  A. Alvermann,et al.  Route to chaos in optomechanics. , 2014, Physical review letters.

[19]  Germany,et al.  Quantum states and phases in driven open quantum systems with cold atoms , 2008, 0803.1482.

[20]  G. Lindblad.  On the generators of quantum dynamical semigroups , 1976 .

[21]  Howard Georgi,et al.  Lie Algebras in Particle Physics , 1982 .

[22]  Steven M. Girvin,et al.  Circuit QED: Superconducting Qubits Coupled to Microwave Photons , 2015 .

[23]  H. Carmichael.  An open systems approach to quantum optics , 1993 .

[24]  H. Neven,et al.  Characterizing quantum supremacy in near-term devices , 2016, Nature Physics.

[25]  Eter,et al.  Faber and Newton Polynomial Integrators for Open-System Density Matrix Propagation , 1998 .

[26]  Dariusz Chruściński,et al.  Universal Spectra of Random Lindblad Operators. , 2018, Physical review letters.

[27]  Antoine Georges,et al.  Interaction-induced impeding of decoherence and anomalous diffusion. , 2012, Physical review letters.

[28]  Niklas Teichmann,et al.  Differences between mean-field dynamics and N-particle quantum dynamics as a signature of entanglement. , 2008, Physical review letters.

[29]  Antonello Scardicchio,et al.  Massively parallel implementation and approaches to simulate quantum dynamics using Krylov subspace techniques , 2017, Comput. Phys. Commun..

[30]  P Hänggi,et al.  Computation of the asymptotic states of modulated open quantum systems with a numerically exact realization of the quantum trajectory method. , 2016, Physical review. E.

[31]  M. Ivanchenko,et al.  Unfolding a quantum master equation into a system of real-valued equations: Computationally effective expansion over the basis of SU(N) generators. , 2018, Physical review. E.

[32]  Pierluigi Amodio,et al.  Parallel solution in time of ODEs: some achievements and perspectives , 2009 .

[33]  Gernot Akemann,et al.  Universal Signature from Integrability to Chaos in Dissipative Open Quantum Systems. , 2019, Physical review letters.

[34]  Franco Nori,et al.  QuTiP 2: A Python framework for the dynamics of open quantum systems , 2012, Comput. Phys. Commun..

[35]  Dariusz Chruscinski,et al.  A Brief History of the GKLS Equation , 2017, Open Syst. Inf. Dyn..