计算化学公社

标题: reaxff反应力场模拟碳氢氧化合物燃烧求助 [打印本页]

作者
Author:
zhang0325    时间: 2023-4-19 17:47
标题: reaxff反应力场模拟碳氢氧化合物燃烧求助
本帖最后由 zhang0325 于 2023-4-19 17:51 编辑

老师们好!我用material studio建立了10个聚丙烯分子和200个氧气分子,用reax/c/species生成species文件,发现结果中聚丙烯和氧气几乎没有发生反应,请问是步数太少还是我的in文件有问题?或者是其他地方不对?下面是in文件和species文件截图,感谢!species
(, 下载次数 Times of downloads: 42)
in文件
(, 下载次数 Times of downloads: 45)


作者
Author:
wildon    时间: 2023-4-19 20:30
初始温度为1000k,nvt的起始温度为1500k。这两个温度最好一致啊,然后温度这么低要等很多步才会开始分解的。可以升高温度试试。
作者
Author:
zhang0325    时间: 2023-4-19 21:41
wildon 发表于 2023-4-19 20:30
初始温度为1000k,nvt的起始温度为1500k。这两个温度最好一致啊,然后温度这么低要等很多步才会开始分解的 ...

好的谢谢您!感谢感谢!!
作者
Author:
Kmetsch    时间: 2023-4-19 21:42
你这一共跑了10ps,根本来不及反应,燃烧的活化能很高的
作者
Author:
zhang0325    时间: 2023-4-19 22:54
Kmetsch 发表于 2023-4-19 21:42
你这一共跑了10ps,根本来不及反应,燃烧的活化能很高的

嗷嗷,好的!那请问大约得跑多少步才可能有反应呢?
作者
Author:
Kmetsch    时间: 2023-4-19 23:12
zhang0325 发表于 2023-4-19 22:54
嗷嗷,好的!那请问大约得跑多少步才可能有反应呢?

我不是这个方向的,但我寻思用活化能算一下速率常数应该能够近似估计出需要的时长吧
作者
Author:
sobereva    时间: 2023-4-20 05:00
zhang0325 发表于 2023-4-19 22:54
嗷嗷,好的!那请问大约得跑多少步才可能有反应呢?

reaxFF模拟燃烧的文章数不胜数,随便找几篇看看便知
作者
Author:
zhang0325    时间: 2023-4-20 11:06
Kmetsch 发表于 2023-4-19 23:12
我不是这个方向的,但我寻思用活化能算一下速率常数应该能够近似估计出需要的时长吧

好的感谢!

作者
Author:
zhang0325    时间: 2023-4-20 11:07
sobereva 发表于 2023-4-20 05:00
reaxFF模拟燃烧的文章数不胜数,随便找几篇看看便知

好的感谢!

作者
Author:
LGX    时间: 2023-4-20 17:25
velocity的温度最好与初始NVT的温度一致,模拟燃烧反应的话看你这个应该是时间不够长,而且根据reaxff而言,模拟温度一般要高于试验温度,可提高温度
作者
Author:
zhang0325    时间: 2023-4-20 20:39
LGX 发表于 2023-4-20 17:25
velocity的温度最好与初始NVT的温度一致,模拟燃烧反应的话看你这个应该是时间不够长,而且根据reaxff而言 ...

好的谢谢!

作者
Author:
cyx98    时间: 2023-4-20 20:51
zhang0325 发表于 2023-4-19 22:54
嗷嗷,好的!那请问大约得跑多少步才可能有反应呢?

通常反应时长100-200ps
作者
Author:
zhang0325    时间: 2023-4-23 18:31
cyx98 发表于 2023-4-20 20:51
通常反应时长100-200ps

感谢!!




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3