This package contains a solver simulating the time evolution of an open quantum system of coupled two-level qubits, modeled by a Lindblad master equation.
The code uses matrix-product-state (MPS) and matrix-product-operator (MPO) data structures, implemented in C++ (using the ITensor library) for maximizing the performance with multithreaded computations, and wrapped by a Python layer with an easy-to-use interface and rich plotting features. The solver integrates in fixed time steps the Markovian master equation for the density matrix,
The solver supports various initial states and observables, with the Hamiltonian and dissipator terms having time-independent coefficients;
How to cite: A paper introducing the solver and presenting a research performed using this solver is currently available as a preprint at arXiv:2203.05871.
Verifying numerical simulations: The above paper also explains in detail how to construct and verify numerical results with multiqubit simulations. We recommend users of the package to familiarize themselves with different approaches to verify their numerical simulations, such as those explained in this paper. The sources used to perform this research are included in the detailed examples.
Note that the C++ binaries of the solver must be built locally - see the installation guide.
For more information, tutorials, and detailed examples check the documentation:
- Installation guide
- Tutorial: Solving the Lindblad dynamics of a qubit ring (go here for a basic example)
- Background and literature on the solver method
- Supported model
- Python interface (API documentation of the Python input, output, and call interface)
- Simulation examples (advanced Python examples for research using the solver)
- C++ solver interface (only relevant for an advanced, low-level usage)
- C++ code structure (only relevant for contributing to the C++ code)