EAM potential by Woo Kyun Kim for the Pt-Au based on the model developed by O'Brien et al. (2018) (EAM_Dynamo_OBrienBarrPrice_2018_PtAu__MO_946831081299_000). Original Method of O'Brien et al. ---------------------------------- 1. They used Morse, exponential and Rose functions to get the pair, electron density (edens) and embedding (embed) functions. All the parameters are presented in the original paper. 2. For the pair and edens, an internal cutoff radius (ri) is introduced so that a cutoff function is applied between the internal cutoff and the global cutoff radii to have continuous derivatives up to 2nd order. 3. Line 5 of the original parameter file has the following values: 5000 1.000000000000e-02 5000 1.150000000000e-03 5.750000000000, which are nrho, drho, nr, dr, rcut. However, I found that these drho (1.0e-2) and dr (1.15e-3) are not correct so used corrected values of drho = max_rhobar / (nrho-1), 1.0002e-2, and dr = rcut / (nr-1), 1.15023e-3, to generate the pair, edens, and embed curves for comparison with my revised ones. max_rhobar (50) was found from the paper. Woo Kyun Kim Modifications: --------------------------- 1. The most important change is that, between ri and rcut, a 7th-order polynomial cutoff function was used to get a continuous derivatives up to the third order. 2. Two parameters in the paper do not appear correct. Using the function values obtained from the original parameter files, they were recovered as: (1) am for Pt-Au: 2.2672956035 -> 2.26719547780531455623948652850535261 (2) rho0 for Au: 1.3619497798 -> 0.509993746380817964219435709369485414 3. In order to remove the discontinuity issue, the global cutoff radius was increased from 5.75 A to 6.65 A.