A bug in -add-shells about atomsk HOT 4 CLOSED

Dear Sir,

Thank you for your message.

Indeed when writing a LAMMPS data file, Atomsk generates new "types" for the shells. If the atom "types" are already defined but not in increasing order, it leads to this bug. I will try to modify the code to fix this bug.

With that said, I always advise against using LAMMPS data format for intermediate files when constructing an atomic system. This file format only stores atom "types" and not necessarily chemical species, which messes things up. Use another file format (like CFG) for intermediate files, and you will have fewer problems. Use the LAMMPS data format only when writing your final output file.

I'll come back here when I solve the problem.

Best regards

from atomsk.

I just pushed a new version of the code. You may download it from GitHub to test it. I would appreciate if you could confirm that it solves your problem.

Best regards

from atomsk.

Thank you for your prompt reply! The problem has been fixed in the new version.

Regarding using other types as intermediate files, I actually tried to use .gin file as the intermediate (since .cfg does not supply the core-shell model). However, since gin does not specify the type number (only by element name), when I convert .gin file to .lmp file, the order of O and La are always reversed because the first atom in the file is O. For the example I provided earlier,

atomsk input.lmp -add-shells all output.gin
atomsk output.gin lmp

This leads to O and La with the type numbers (1,3) and (2,4), respectively. Although it is not the result I expected, I am not sure whether this should be considered a bug. Just listed here for your reference.

Thank you for your quick reply and prompt fix!

from atomsk.

I apologize for the delayed answer, but I understand this issue is fixed now so I will close it.

Concerning your last remark, it is not a bug. When atom "types" are not defined, Atomsk generates new types in order of appearance of the atoms. If O atoms appear first then they are assigned type 1, and so on. To fix this, just change the order of atoms, or write atom types in a file and use option "-properties" to assign types.

from atomsk.