Solving the multi-level Maxwell-Bloch equations using GPGPU computing for the simulation of nonlinear optics in atomic gases

Abstract

We present a numerical implementation of a solver for the Maxwell-Bloch equations to calculate the propagation of a light pulse in a nonlinear medium composed of an atomic gas in one, two and three dimensional systems. This implementation solves the wave equation of light using a finite difference method in the time domain scheme, while the Bloch equations for the atomic population in each point of the simulation domain are integrated using splitting methods. We present numerical simulations of atomic-gas systems and performance benchmarks.