Charge optimized many body (COMB) potentials for Pt and Au

Evaluation of Force Fields for Molecular Dynamics Simulations of Platinum in Bulk and Nanoparticle Forms

Understanding the size- and shape-dependent properties of platinum nanoparticles is critical for enabling the design of nanoparticle-based applications with optimal and potentially tunable functionality. Toward this goal, we evaluated nine different empirical potentials with the purpose of accurately modeling faceted platinum nanoparticles using molecular dynamics simulation. First, the potentials were evaluated by computing bulk and surface properties-surface energy, lattice constant, stiffness constants, and the equation of state-and comparing these to prior experimental measurements and quantum mechanics calculations. Then, the potentials were assessed in terms of the stability of cubic and icosahedral nanoparticles with faces in the {100} and {111} planes, respectively. Although none of the force fields predicts all the evaluated properties with perfect accuracy, one potential-the embedded atom method formalism with a specific parameter set-was identified as best able to model platinum in both bulk and nanoparticle forms.

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid-liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. The use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

Effects of surface charge and cluster size on the electrochemical dissolution of platinum nanoparticles using COMB3 and continuum electrolyte models

We study the site-dependent dissolution of platinum nanoparticles under electrochemical conditions to assess their thermodynamic stability as a function of shape and size using empirical molecular dynamics and electronic-structure models. The third-generation charge optimized many-body potential is employed to determine the validity of uniform spherical representations of the nanoparticles in predicting dissolution potentials (the Kelvin model). To understand the early stages of catalyst dissolution, implicit solvation techniques based on the self-consistent continuum solvation method are applied. It is demonstrated that interfacial charge and polarization can shift the dissolution energies by amounts on the order of 0.74 eV depending on the surface site and nanoparticle shape, leading to the unexpected preferential removal of platinum cations from highly coordinated sites in some cases.

Nanoscale Structure and Dynamics of Water on Pt and Cu Surfaces from MD Simulations

The interaction of liquid water with Pt(111) is investigated with classical molecular dynamics (MD) simulations, where the forces are determined using the third-generation charge optimized many-body (COMB3) interatomic potential. In cases of sub-monolayer water coverage, the parameterized empirical potential predicts experimentally observed and energetically favorable √37 and √39 reconstructed water structures with "575757" di-interstitial defects. At both sub-monolayer and multilayer water coverages, the structure of the first wetting layer of liquid water on Pt(111) exhibits a characteristic distribution where the molecules form two distinct buckled layers as a result of the interplay between water-metal adsorption and water-water hydrogen bonds. The dynamic spreading rate of water nanodroplets on large Pt surfaces (>200 nm²) characterized by molecular kinetic spreading theory is an order of magnitude slower than the molecular kinetic rate of the same droplet on close-packed Cu surfaces due to variation in molecular distributions at the water-metal interface. These nanoscale MD simulation predictions using the COMB3 interatomic potential demonstrate the capability of capturing both many-body interactions between H₂O and Pt or Cu and hydrogen bonding in liquid water.