Monte Carlo simulation of the atomic master equation for spontaneous emission

A Monte Carlo simulation of the atomic master equation for spontaneous emission in terms of atomic wave functions is developed. Realizations of the time evolution of atomic wave functions are constructed that correspond to an ensemble of atoms driven by laser light undergoing a sequence of spontaneous emissions. The atomic decay times are drawn according to the photon count distribution of the driven atom. Each quantum jump of the atomic electron projects the atomic wave function to the ground state of the atom. Our theory is based on a stochastic interpretation and generalization of Mollow's pure-state analysis of resonant light scattering, and the Srinivas-Davies theory of continuous measurements in photodetection. An extension of the theory to include mechanical light effects and a generalization to atomic systems with Zeeman substructure are given. We illustrate the method by simulating the solutions of the optical Bloch equations for two-level systems, and laser cooling of a two-level atom in an ion trap where the center-of-mass motion of the atom is described quantum mechanically.