Modeling the sorption dynamics of NaH using a reactive force field

We have parametrized a reactive force field for NaH, ReaxFFNaH, against a training set of ab initio derived data. To ascertain that ReaxFFNaH is properly parametrized, a comparison between ab initio heats of formation of small representative NaH clusters with ReaxFFNaH was done. The results and trend of ReaxFFNaH are found to be consistent with ab initio values. Further validation includes comparing the equations of state of condensed phases of Na and NaH as calculated from ab initio and ReaxFFNaH. There is a good match between the two results, showing that ReaxFFNaH is correctly parametrized by the ab initio training set. ReaxFFNaH has been used to study the dynamics of hydrogen desorption in NaH particles. We find that ReaxFFNaH properly describes the surface molecular hydrogen charge transfer during the abstraction process. Results on heat of desorption versus cluster size shows that there is a strong dependence on the heat of desorption on the particle size, which implies that nanostructuring enhances desorption process. To gain more insight into the structural transformations of NaH during thermal decomposition, we performed a heating run in a molecular dynamics simulation. These runs exhibit a series of drops in potential energy, associated with cluster fragmentation and desorption of molecular hydrogen. This is consistent with experimental evidence that NaH dissociates at its melting point into smaller fragments.