Molecular dynamics simulation of disorder-induced heating in ultracold neutral plasma

Disorder-induced heating (DIH) is one of the main reasons reducing the coupling strength in ultracold plasma. We propose applying an optical lattice as periodic confinement before the ultracold atomic cloud is ionized to eliminate its effect. We demonstrate a numerical simulation for the dynamics of the ultracold plasmas using classical molecular dynamics method with open boundary. DIH is reproduced in the simulation with the random Gaussian initial distribution and is absent in the results with the ordered lattice initial distribution. We further find that the collisional heating from electrons is important for ultracold plasmas with chosen spatial correlations in the optical lattice. Carefully preparing the initial condition (e.g., the ion density, initial electron temperature, and so on), collisional heating for the ions would be significantly reduced and eventually negligible. This allows a much stronger coupling in ultracold plasma to be realized.