RESUMO
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.
RESUMO
Quantum systems out of equilibrium are presently a subject of active research, both in theoretical and experimental domains. In this work, we consider time-periodically modulated quantum systems that are in contact with a stationary environment. Within the framework of a quantum master equation, the asymptotic states of such systems are described by time-periodic density operators. Resolution of these operators constitutes a nontrivial computational task. Approaches based on spectral and iterative methods are restricted to systems with the dimension of the hosting Hilbert space dimH=Nâ²300, while the direct long-time numerical integration of the master equation becomes increasingly problematic for Nâ³400, especially when the coupling to the environment is weak. To go beyond this limit, we use the quantum trajectory method, which unravels the master equation for the density operator into a set of stochastic processes for wave functions. The asymptotic density matrix is calculated by performing a statistical sampling over the ensemble of quantum trajectories, preceded by a long transient propagation. We follow the ideology of event-driven programming and construct a new algorithmic realization of the method. The algorithm is computationally efficient, allowing for long "leaps" forward in time. It is also numerically exact, in the sense that, being given the list of uniformly distributed (on the unit interval) random numbers, {η_{1},η_{2},...,η_{n}}, one could propagate a quantum trajectory (with η_{i}'s as norm thresholds) in a numerically exact way. By using a scalable N-particle quantum model, we demonstrate that the algorithm allows us to resolve the asymptotic density operator of the model system with N=2000 states on a regular-size computer cluster, thus reaching the scale on which numerical studies of modulated Hamiltonian systems are currently performed.