Additional AIMD calculations were performed at 300 K with an N V T ensemble with the crystalline structure. This run was to verify the CMD results using the developed force field at 300 K. After an equilibrium run for 10 ps at 300 K, the atomic configurations were sampled for 50 ps.