LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] Mixing EAM potentilas to create alloy potential in LAMMPS
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[lammps-users] Mixing EAM potentilas to create alloy potential in LAMMPS

From: ankur agrawal <ankur.agrawal1103@...24...>
Date: Tue, 19 Dec 2017 11:01:57 +0530

I am facing a problem in making a alloy EAM potential for Ni-Co by using the individual element EAM potentials from NIST database, and merging them using averaging of phi interactions
I have attached the 2 individual EAM potential files. For combining them, I copy pasted them one below the another. Another thing that I kept in mind was that the Nr phi terms were removed from their place in individual potentials and added at last along with phi 2-1 calculated by averaging in excel.
One of the problems is that in this setfl format, we are allowed to specify Npho, dpho, Nr and dr only once in the beginning. They will be same for all the elements present.
But, the individual files for Ni and Co that are available on NIST have different values of these distances  dr and dpho. In these case how to mix these potentials?
For time being I just used the values of Co in the file, but I am getting errors ( lost atoms) while relaxing the alloy structure. I think this is a bad potential that I have created.
Can somebody please guide me on mixing EAM potentials to create alloy potential?
Yours Sincerely
Ankur Agrawal
FifthYear Dual Degree Student
Department of Chemical Engineering
Indian Institute of Technology, kanpur, India


Attachment: Co_PurjaPun_2012.eam.alloy
Description: Binary data

Attachment: Ni1_Mendelev_2012.eam.fs
Description: Binary data