Adiabatic elimination for composite open quantum systems: Reduced-model formulation and numerical simulations
François-Marie Le Régent, Pierre Rouchon
Abstract
A numerical method is proposed for simulation of composite open quantum systems. It is based on Lindblad master equations and adiabatic elimination. Each subsystem is assumed to converge exponentially towards a stationary subspace, slightly impacted by some decoherence channels and weakly coupled to the other subsystems. This numerical method is based on a perturbation analysis with an asymptotic expansion. It exploits the formulation of the slow dynamics with reduced dimension. It relies on the invariant operators of the local and nominal dissipative dynamics attached to each subsystem. Second-order expansion can be computed only with local numerical calculations. It avoids computations on the tensor-product Hilbert space attached to the full system. This numerical method is particularly well suited for autonomous quantum error correction schemes. Simulations of such reduced models agree with complete full model simulations for typical gates acting on one and two cat qubits (z, zz, and cnot) when the mean photon number of each cat qubit is less than eight. For larger mean photon numbers and gates with three cat qubits (zzz and ccnot), full model simulations are almost impossible whereas reduced model simulations remain accessible. In particular, one observes numerically the simultaneous capture of the dominant phase-flip error rate and of the very small bit-flip error rate with its exponential suppression versus the mean photon number.