We describe a hybrid molecular dynamics approach for the description of
ultracold neutral plasmas, based on an adiabatic treatment of the electron gas
and a full molecular dynamics simulation of the ions, which allows us to follow
the long-time evolution of the plasma including the effect of the strongly
coupled ion motion. The plasma shows a rather complex relaxation behavior,
connected with temporal as well as spatial oscillations of the ion temperature.
Furthermore, additional laser cooling of the ions during the plasma evolution
drastically modifies the expansion dynamics, so that crystallization of the ion
component can occur in this nonequilibrium system, leading to lattice-like
structures or even long-range order resulting in concentric shells