Optimization and application of lithium parameters for the reactive force field

To make a practical molecular dynamics (MD) simulation of the large-scale reactive chemical systems of Li-H and Li-C, we have optimized parameters of the reactive force field (ReaxFF) for these systems. The parameters for this force field were obtained from fitting to the results of density functional theory (DFT) calculations on the structures and energy barriers for a number of Li-H and Li-C molecules, including Li2, LiH, Li2H2, H3C-Li, H3C-H2C-Li, H2CdC-LiH, HCtCLi, H6C5-Li, and Li2C2, and to the equations of state and lattice parameters for condensed phases of Li. The accuracy of the developed ReaxFF was also tested by comparison to the dissociation energies of lithium-benzene sandwich compounds and the collision behavior of lithium atoms with a C60 buckyball. 1