|From:||Tony Liu <tonyliu0224@...24...>|
|Date:||Fri, 16 Jun 2017 13:42:48 +0900|
On Wed, Jun 14, 2017 at 1:43 AM, Tony Liu <tonyliu0224@...24...> wrote:
Dear Lammps Community,
Please forgive me to again report a strange phenomenon resulted from fix reax/c/species command.
I built an bicrystal alloy model consisting of total 3010 atoms.
By the atom type setting command: set group alloy type/fraction 1 0.82 12963 type/fraction 2 0.13 23486 type/fraction 3 0.05 23486,
I intend to set 82% Fe atom, 13% Cr atom and 5% Al atom.
And the log.lammps shows the following reasonable results:
2454 settings made for type/fraction
397 settings made for type/fraction
143 settings made for type/fraction
However, the fix reax/c/species command output a different result: Fe2118Cr251Al624.
Although I enabled the reax/c with kokkos package by the following changes, the same problem occurred again.yes, because the problem is *not* with fix reax/c/species but falls into the PEBCAC category, i.e. it resides between the chair and the computer.the "fraction" is always applied to the group of the set command (which is all atoms). so with the second type/fraction set will *also* consider atoms, that have already been changed to type 1. you can confirm, that fix reax/c/species is working correctly, by issuing the following 3 commands and checking the log file:group Fe type 12131 atoms in group Fegroup Cr type 2254 atoms in group Crgroup Al type 3625 atoms in group Alto explain. the first type/fraction command does change 2454 atoms from type 3 to type 1, the second command changes 397 atoms to type 2, however quite a few of them were of type 1 and not type 3. the last command makes things even worse, as you now change 5% of all 3010 atoms to type 3, and - again - it is most likely that you change atoms of type 1 or 2.the correct way to go about this, is the following:set group all type/fraction 1 0.82 12963# now there are about 18% atoms left as type 3# you want a ratio of 13:5 Cr to Al, i.e. a fraction of 13/18 = 0.72222 needs to be converted from Cr to Al# make a group that only contains type 3 atoms.group CrAl type 3set group CrAl type/fraction 2 0.722222222 23486group Fe type 12454 atoms in group Fegroup Cr type 2394atoms in group Crgroup Al type 3162atoms in group Al# that is: 81.5% Fe, 13.1% Cr and 5.4% Alor even better, do it the other way around, i.e. create all atoms as type 1, then convert 18% of those to type 2 and 27.7777% of the type 2 to type 3.with that i get:set group all type/fraction 2 0.18 12963Setting atom values ...508 settings made for type/fractiongroup CrAl type 2508 atoms in group CrAlset group CrAl type/fraction 3 0.2777777777 23486Setting atom values ...149 settings made for type/fractiongroup Fe type 12502 atoms in group Fegroup Cr type 2359 atoms in group Crgroup Al type 3149 atoms in group Al-> 83.1% Fe, 11.9% Cr, 5% Almystery solved. LAMMPS is working fine, but your input script needs an upgrade. ;-)axel.
e.g. pair_style reax/c into pair_style reax/c/kk & fix reax/c/species into fix reax/c/species/kk
My lammps information is as follows:
Lammps version: 31Mar17;
Packages: Kokkos, User-reaxc, User-misc, molecule, qeq, reax, replica, rigid, kspace, manybody
Running command: mpirun -np 4 -k on -sf kk -in in.mini
How do you like such problem?
Any advices will be highly appreciated! I look forward to your positive response,
Thanks in advance,
Xiaolong Liu from Korea Aerospace University.Fax: 82-2-3158-2191Tel.: 82-2-300-005776, Hanggongdaehang Road, Deokyang District, Goyang, Gyeonggi Province, 412-791, R.O. Korea.Korea Aerospace UniversityDepartment of Aerospace and Mechanical EngineeringXiaolong LiuPh.D student
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
lammps-users mailing list