On fitting a gold embedded atom method potential using the force matching method.

We fit a new gold embedded atom method (EAM) potential using an improved force matching methodology which included fitting to high-temperature solid lattice constants and liquid densities. The new potential shows a good overall improvement in agreement to the experimental lattice constants, elastic constants, stacking fault energy, radial distribution function, and fcc/hcp/bcc lattice energy differences over previous potentials by Foiles, Baskes, and Daw (FBD) [Phys. Rev. B 33, 7983 (1986)] Johnson [Phys. Rev. B 37, 3924 (1988)], and the glue model potential by Ercolessi et al. [Philos. Mag. A 50, 213 (1988)]. Surface energy was improved slightly as compared to potentials by FBD and Johnson but as a result vacancy formation energy is slightly inferior as compared to the same potentials. The results obtained here for gold suggest for other metal species that further overall improvements in potentials may still be possible within the EAM framework with an improved fitting methodology. On the other hand, we also explore the limitations of the EAM framework by attempting a brute force fit to all properties exactly which was found to be unsuccessful. The main conflict in such a brute force fit was between the surface energy and the liquid lattice constant where both could not be fitted identically. By intentionally using a very large number of spline sections for the pair potential, electron-density function, and embedding energy function, we eliminated a lack of functional freedom as a possible cause of this conflict and hence can conclude that it must result from a fundamental limitation in the EAM framework.