Recently I am trying to apply fix GCMC with EAM potential. As you mentioned, “EAM potentials for metals only include the pair potential portion of the EAM interaction, not the embedding term.”
Why the embedding term does not include currently? For the sake of implement of parallel computing? Or any other reason? Meanwhile, if the embedding term was not included, does the simulation result make sense as you expected?