1. lammpstrj文件
lammpstrj文件应包含周期性、每帧原子数、各原子的id、type_id(原子类型序号)、x、y、z字段。最简单的方式是直接用dump atom指令输出。
ITEM: TIMESTEP
1000
ITEM: NUMBER OF ATOMS
10000
ITEM: BOX BOUNDS pp pp pp
0.0 40.0
0.0 40.0
0.0 40.0
ITEM: ATOMS id type x y z
1 1 0.000 0.000 0.000
2 2 1.000 1.000 1.000
...
id、type_id、x、y、z的左右顺序可以改变。也可以是id、type_id、xs、ys、zs。
由于type_id是数字而不是元素符号,所以在运行ReaxTools时还要加上如-t C,H,O,N定义对应元素符号。
2. xyz (extended)文件
xyz (extended)文件同样应包含周期性、每帧原子数、各原子的id、type(元素符号)、x、y、z字段。
xyz (extended)文件获得的方式包括: GPUMD和CP2K直接输出、OVITO转格式。
10000
Lattice="40.000 0.0 0.0 0.0 40.000 0.0 0.0 0.0 40.000" Properties=species:S:1:pos:R:3
1 C 0.000 0.000 0.000
2 H 1.000 1.000 1.000
...
当xyz文件中不含id信息时,请务必确认各帧原子顺序是一致的。
顺序一致
10000
Lattice="40.000 0.0 0.0 0.0 40.000 0.0 0.0 0.0 40.000" Properties=species:S:1:pos:R:3
C 0.000 0.000 0.000
H 1.000 1.000 1.000
...
顺序不一致,转格式等过程中顺序紊乱
10000
Lattice="40.000 0.0 0.0 0.0 40.000 0.0 0.0 0.0 40.000" Properties=species:S:1:pos:R:3
H 5.000 6.000 7.000
H 1.000 1.000 1.000
...
由于type已经是元素符号,所以不用加上-t选项。
3. xyz文件
xyz文件与xyz (extended)文件的要求相同,只是缺少周期性信息,这种情况下ReaxTools的运行速度更慢,也不会处理跨周期的键。
无周期性信息的xyz文件一般是从ORCA、CP2K进行非周期性分子体系的AIMD时输出的。
当然,由于xyz已经是最简的坐标格式了,所以绝大多数软件的轨迹都可以转为xyz文件,如果稍加修改,添加上Lattice="..."这行,那么就是xyz (extended)格式
10000
Any comment or other information
1 C 0.000 0.000 0.000
2 H 1.000 1.000 1.000
...
同样地,当xyz文件中不含id信息时,请务必确认各帧原子顺序是一致的。
由于type已经是元素符号,所以不用加上-t选项。