计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: Graphite
打印 Print 上一主题 Last thread 下一主题 Next thread

[Lammps] 利用Lammps ReaxFF研究反应动力学一例

  [复制链接 Copy URL]

37

帖子

0

威望

586

eV
积分
623

Level 4 (黑子)

61#
发表于 Post on 2025-5-22 21:41:19 | 只看该作者 Only view this author
各位老师好,能否请教以下对于参照上述帖子已经获得的lammps.trj轨迹如果rerun调整bond_order cutoffs来重新得到species,在rerun处应该怎么写呢,我根据官网和上述模板在去掉能量最小化和初始速度后对run部分作出了如下修改,一直报错“ERROR: Read_dump must use at least either 'id' or 'type' field (../read_dump.cpp:1175)“,请问正确的方法是什么?
(1)直接将run部分修改为“rerun ${basename}.lammpstrj dump every 1 id type xs ys zs”
(2)将trj文件转化为xyz,将run部分修改为:“rerun ${basename}.xyz dump every 1  x y z” ;
谢谢!

395

帖子

8

威望

3906

eV
积分
4461

Level 6 (一方通行)

石墨

62#
 楼主 Author| 发表于 Post on 2025-5-23 04:21:04 | 只看该作者 Only view this author
本帖最后由 Graphite 于 2025-5-23 04:22 编辑
love_yy 发表于 2025-5-22 21:41
各位老师好,能否请教以下对于参照上述帖子已经获得的lammps.trj轨迹如果rerun调整bond_order cutoffs来重 ...

直接用reax_tools或者reacnetgenerator读轨迹分析得了,说实话基于LAMMPS键级和cutoff,有时候反而结果相当的紊乱

reax_tools用vdw缩放因子调意外地又LOW又好用
自在飞花轻似梦,无边丝雨细如愁。

全自动反应动力学(ReaxFF、AIMD、NEP等)后处理工具网页版:http://cc-portal.xyz/reax_tools

37

帖子

0

威望

586

eV
积分
623

Level 4 (黑子)

63#
发表于 Post on 2025-5-23 12:36:16 | 只看该作者 Only view this author
Graphite 发表于 2025-5-23 04:21
直接用reax_tools或者reacnetgenerator读轨迹分析得了,说实话基于LAMMPS键级和cutoff,有时候反而结果相 ...

谢谢老师回复

1

帖子

0

威望

9

eV
积分
10

Level 1 能力者

64#
发表于 Post on yesterday 13:49 | 只看该作者 Only view this author
老师您好,请问reax_species.py是否有下载渠道呢

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-12 19:04 , Processed in 0.142116 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list