Hi. Friends
I m study on Al2O3 by simulation in lammps. the Potential function is MEAM.
I want to input the “re” parameter of MEAM, for single element of O and Al in library file. but in library.meam file there is no choice to input this parameter and in AlO.meam I can just input the re(I,J) format for Al-O interaction.
I input the “re” for single element by following format to understand what should I do for this.
re(1)=value , re(2)=value format.
r(1,1)=value , r(2,2)= value format.
without the “re” for single element.
In all of 3 cases the lammps don’t show me any error and all of the answers are different from each other.
please help me how can I input the “re” for single element in MEAM potential.
Hi. Friends
I m study on Al2O3 by simulation in lammps. the Potential function is MEAM.
I want to input the "re" parameter of MEAM, for single element of O and Al in library file. but in library.meam file there is no choice to input this parameter and in AlO.meam I can just input the re(I,J) format for Al-O interaction.
I input the "re" for single element by following format to understand what should I do for this.
re(1)=value , re(2)=value format.
r(1,1)=value , r(2,2)= value format.
without the "re" for single element.
In all of 3 cases the lammps don't show me any error and all of the answers are different from each other.
please help me how can I input the "re" for single element in MEAM potential.
We specify the "re" parameter for single element as "alat" in the "meamf" file.
For the cross (unlike) pairs, e.g. re(Al,O), specify re(1,2) in the
"meam.alsimgcufe".
Dear Bohumir
many thanks for your answer.
1- I am so apologize for write the r(1,1) instead of re(1,1).
I used the re(1,1) and re(2,2) in AlO.meam, and the lammps didn't show any error.
2- in the lammps documentation the "alat" parameter in "meamf" is different from "re". for example
the value of alat for Aluminum is 4.05 and the re is 2.86 (www.osti.gov/energycitations/servlets/purl/224267-allunm/webviewable/Modifiedembeddedatommethodcalculationsofinterfaces.pdf).
3- in my experience, the result of simulation is not depend to alat parameter. because I use different values for Oxygen and the result was quite the same.
I am very confuse about this and please help me.
"re" is the equilibrium nearest neighbor distance for Al
"alat" is the lattice parameter for Al
That is why they differ. Notice that 4.05/sqrt(2.0) ~= 2.86.
Likely, results of your simulations do not depend on "alat" parameter
because "re" specified in the MEAM parameter file ("AlO.meam") will
override the "alat" value from the MEAM library file ("meamf"). That
would follow the statement from http://lammps.sandia.gov/doc/pair_meam.html
"If used, the MEAM parameter file contains settings that override or
complement the library file settings."
For a single element interaction, the "re" parameter (which is the equilibrium atomic spacing) can be set in the library file through the lattice parameter (alat). How this gets computed depends on the lattice type. For example, for an FCC crystal, re = alat/sqrt(2). For a dimer structure, which might be how you're specifying oxygen, re = alat.
Although there should be no need, if you want you can override this value by specifying something in the potential file. In that case, the correct syntax is to provide two indices, e.g. re(1,1)=value. I don't think there's any error checking done if you don't provide both indices, so the behavior may be undefined if you specify re(1)=value.
Regards,
Greg
Greg Wagner
Manager, Thermal/Fluid Science and Engineering Department
Sandia National Laboratories, Livermore, CA
Tel: (925) 294-2180 Fax: (925) 294-3410
Email: [email protected]
