Abstract
A systematic numerical approach to approximate high-dimensional Lindblad equations is described. It is based on a deterministic rank approximation of the density operator, the rank being the only parameter to adjust. From a known initial density operator, this rank approximation gives at each time step an estimate of its largest eigenvalues with their associated eigenvectors. A numerical integration scheme is also proposed. Its numerical efficiency in the case of a rank approximation is demonstrated for oscillation revivals of 50 atoms interacting resonantly with a slightly damped coherent quantized field of 200 photons.
- Received 14 November 2012
DOI:https://doi.org/10.1103/PhysRevA.87.022125
©2013 American Physical Society