ReaxTools输入文件格式说明

返回工具主页

ReaxTools支持lammpstrj、xyz (extended)、xyz格式的输入文件。

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选项。