Ab initio Molecular Dynamics for Determination of Structures of Alkali Metal Clusters and Their Temperatures Behavior; An Example of Li9+

It will be shown that an ab initio molecular dynamics procedure based on gradient corrected density functionals for exchange and correlation and using a Gaussian atomic basis (AIMD-GDF) implemented for parallel processing represents a suitable tool for detailed and accurate investigation of structural and dynamical properties of small systems. Gradients of the Born-Oppenheimer ground state energy, obtained by iterative solution of the Kohn-Sham equations, are used to calculate the forces acting on atoms at each instantaneous configuration along trajectories generated by solving classical equations of motion. Dynamics of different isomers of the Li9+ cluster have been investigated as a function of excess energy. It is shown that different isomers, even those similar in energy, can exhibit different structural and dynamical behavior. The analysis of the simulations leads to the conclusion that structures with a central atom, in particular the centered antiprism of Li9+ exhibit concerted mobility of the peripheral atoms at relatively low excess energy. In contrast, compact tetrahedral type structures show much more rigid behavior at low excess energy. However, the former ones need larger excess of internal energy to undergo isomerizations to geometrically different structures than the latter ones. At the time scale of our simulations we found that for the intermediate excess energies it is "easier" to carry the cluster in the basin of the lowest energy isomer than in the reverse direction. It has been found that the liquid-like behavior in small Li clusters becomes apparent at relatively high temperature in spite of large mobility of their atoms.