MD simulations were performed within the AMS driver program 56 using the Velocity-Verlet algorithm with a time step of 0.25 fs. The geometry of the system was optimized before each MD simulation until the norms of all forces on all atoms were below 0.1 eV Å-1 to avoid simulation artifacts due to high-energy starting configurations.