I use this syntax to create voids in the structure , for ex:
region void0 sphere 20 20 25 0.5
delete_atoms region void0
However, is there any way to create precipitates in the structure. I am looking to add a group of interstitial atoms in the structure or a region of interstitial atoms ? I don’t know if this is possible or not, I am new to LAMMPS
I use the EAM potential with zbl.
My concern regarding the pair_coeff command. As I understand that atomic number for two atom types should be associated with the pair coefficient command. I use the atomic number of 26 and 24 for both Fe and Cr, respectively. Is this correct to use? is there any restrictions to this command while using the atomic number?
This question is covered in the documentation, under Mixing. You script contains a variety of syntax and conceptual errors. you should instead do something like this (I did not test this, so it may still be incorrect):
hybrid/overlay adds all potential contributions together
pair_style hybrid/overlay eam/cd zbl 1.0 2.05
with hybrid or hybrid/overlay, each argument list must begin with the pair_style name
pair_coeff * * eam/cd …/potentials/FeCr.cdeam Fe Cr
I changed the script and it works fine using hybrid/overlay for mixed pairs.
I have one last question, I am doing zbl to harden/ stiffen the potential, however, the Fe-Fe part is already stiffened in the potential. Do I need to apply the zbl to Fe-Fe part again? or I can just do zbl for Fe-Cr and Cr-Cr parts since Fe-Fe is already stiffened?
in addition, you mentioned that the zbl coefficients are required for all pure pairs.
I can not tell you how to build a correct model. You will have to figure that out yourself. LAMMPS does require that all pure pairs be specified. If you wish the ZBL interaction to be zero for Fe-Fe, the code is not set up to support that. You can trick it by making the Z values small, but not too small, so that for the smallest separations encountered in the simulation, the function phi (see doc page) evaluates to a sufficiently small value (say 1e-6).