A fast, high-order numerical method for the simulation of single-excitation states in quantum optics

Abstract

We consider the numerical solution of a nonlocal partial differential equation which models the process of collective spontaneous emission in a two-level atomic system containing a single photon. We reformulate the problem as an integro-differential equation for the atomic degrees of freedom, and describe an efficient solver for the case of a Gaussian atomic density. The problem of history dependence arising from the integral formulation is addressed using sum-of-exponentials history compression. We demonstrate the solver on two systems of physical interest: in the first, an initially-excited atom decays into a photon by spontaneous emission, and in the second, a photon pulse is used to an excite an atom, which then decays.

References in corpus (1)

• The Quantum Internet

Cited by in corpus (4)

• A fast time domain solver for the equilibrium Dyson equation
• Decomposing imaginary time Feynman diagrams using separable basis functions: Anderson impurity model strong coupling expansion
• Transient dynamics of subradiance and superradiance in open optical ensembles
• Analysis of single-excitation states in quantum optics