Development and initial applications of an e-ReaxFF description of Ag nanoclusters

Metal nanocrystals are of considerable scientific interest because of their uses in electronics, catalysis, and spectroscopy, but the mechanisms by which nanocrystals nucleate and grow to achieve selective shapes are poorly understood. Ab initio calculations and experiments have consistently shown that the lowest energy isomers for small silver nanoparticles exhibit two-dimensional (2D) configurations and that a transition into three-dimensional (3D) configurations occurs with the addition of only a few atoms. We parameterized an e-ReaxFF potential for Ag nanoclusters (N ≤ 20 atoms) that accurately reproduces the 2D-3D transition observed between the Ag₅ and Ag₇ clusters. This potential includes a four-body dihedral term that imposes an energetic penalty to 3D structures that is significant for small clusters but is overpowered by the bond energy from out-of-plane Ag-Ag bonds in larger 3D clusters. The potential was fit to data taken from density-functional theory and coupled-cluster calculations and compared to an embedded atom method potential to gauge its quality. We also demonstrate the potential of e-ReaxFF to model redox reactions in silver halides and plasmon motion using molecular dynamics simulations. This is the first case in which e-ReaxFF is used to describe metals. Furthermore, the inclusion of a bond-order dependent dihedral angle in this force field is a unique solution to modeling the 2D-3D transition seen in small metal nanoclusters.