Ab initio recombination in the evolving ultracold plasmas

Efficiency of recombination is of crucial importance for the existence of ultracold plasmas (UCP), particularly, the ones formed in magneto-optical traps. Unfortunately, an equilibrium thermodynamic treatment of the ionization-recombination processes is inappropriate for the evolving UCP clouds, while a straightforward kinetic simulation encounters a problem of huge difference in spatial and temporal scales for free and bound motion of the electrons. As a result, only the 'virtual' electron-ion pairs are usually reproduced in such modeling, and it is necessary to employ some heuristic criteria to identify them with recombined atoms. The aim of this paper is to present the first successful ab initio simulation of a non-equilibrium recombination in the evolving UCP plasmas. We employ a special algorithm, which is based on using the 'scalable' reference frame, co-moving with an expanding substance. Then, the recombination events are identified by a series of sharp equidistant peaks in the kinetic and/or potential energies, which are caused by captured electrons passing near the pericenters of their orbits; and this is confirmed by a detailed inspection of their trajectories. Thereby, we are able to trace the real, rather than 'virtual', electron-ion pairs; and the total efficiency of their formation was found to be about 20%, which is in agreement with laboratory measurements.