计算化学公社

 找回密码 Forget password
 注册 Register
Views: 38412|回复 Reply: 23
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 关于ReaxFF力场文件的疑惑

[复制链接 Copy URL]

256

帖子

0

威望

1124

eV
积分
1380

Level 4 (黑子)

在reaxff反应力场中,涉及到很多的参数,由于要研究的体系没有现成的力场参数可用,所以我想了解力场中各个参数是什么意思,以及这些参数是如何得到的。看了reaxFFl力场的手册后,解决了部分问题,但是还是有些问题不知道,期望能得到前辈们的帮助,谢谢!
(1)在图1中的39个general parameters中,有部分数值后面标记的是not used,既然是没有用到,为什么还要写出来呢?
(2)在图2中,例如像van der Waals dissociation energy(Evdw)、Undercoordination energy(Eunder)等的数值是怎么计算出来的,是通过高斯吗?
(3)在reaxFF力场的开发中,作者会提到用量化计算一堆训练集,然后进行参数拟合就得到了反应力场中的文件,可否有相关的资料推荐,非常感谢!


1.jpg (213.55 KB, 下载次数 Times of downloads: 80)

图1

图1

2.jpg (108.55 KB, 下载次数 Times of downloads: 71)

图2

图2

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

2#
发表于 Post on 2017-5-24 02:15:29 | 只看该作者 Only view this author
1 not used可能是开发者搞参数时候临时搞的一些参数,但没太大实用性,所以发布的时候就给隐藏掉了
2 你可以看看原文介绍的细节。力场的vdW参数一般都是基于实验结果和模拟结果,反复优化出来的,极少是直接通过高精度量化算的(除非那种极小体系,才是直接基于高精度算的相互作用势拟合的)。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

256

帖子

0

威望

1124

eV
积分
1380

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2017-5-24 08:51:18 | 只看该作者 Only view this author
sobereva 发表于 2017-5-24 02:15
1 not used可能是开发者搞参数时候临时搞的一些参数,但没太大实用性,所以发布的时候就给隐藏掉了
2 你可 ...

好的,谢谢社长

99

帖子

0

威望

1092

eV
积分
1191

Level 4 (黑子)

4#
发表于 Post on 2017-5-24 15:10:38 | 只看该作者 Only view this author
可以参考ReaxFF发明者Adri C. T. van Duin的文章,例如这篇(http://pubs.acs.org/doi/abs/10.1021/acs.jpclett.6b02902)做了MoS2的力场拟合。我在用ReaxFF做一些粒子辐照方面的模拟,有空可以多交流

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

256

帖子

0

威望

1124

eV
积分
1380

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2017-5-25 10:25:02 | 只看该作者 Only view this author
qwoop 发表于 2017-5-24 15:10
可以参考ReaxFF发明者Adri C. T. van Duin的文章,例如这篇(http://pubs.acs.org/doi/abs/10.1021/acs.jpc ...

好的,谢谢您,请问您在所做的体系,需要自己拟合参数吗?http://pubs.rsc.org/en/content/articlehtml/1994/ft/ft9949002881 这篇文献中介绍到参数拟合方法,不知您是否看懂?

222

帖子

0

威望

4751

eV
积分
4973

Level 6 (一方通行)

6#
发表于 Post on 2017-5-25 15:34:36 | 只看该作者 Only view this author
可以看看potfit是怎么做的

256

帖子

0

威望

1124

eV
积分
1380

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2017-5-25 15:44:23 | 只看该作者 Only view this author
get-it 发表于 2017-5-25 15:34
可以看看potfit是怎么做的

好的,谢谢您,我马上去看看

230

帖子

0

威望

2653

eV
积分
2883

Level 5 (御坂)

8#
发表于 Post on 2017-6-3 01:39:06 | 只看该作者 Only view this author
get-it 发表于 2017-5-25 15:34
可以看看potfit是怎么做的

好像不支持ReaxFF啊

230

帖子

0

威望

2653

eV
积分
2883

Level 5 (御坂)

9#
发表于 Post on 2017-6-3 02:07:56 | 只看该作者 Only view this author
本帖最后由 一颗赛艇 于 2017-6-3 02:09 编辑

发一下自己对力场参数的研究。时间过去太长,有一些也记得不太清楚了。

灰色的在拟合中应该是不变的,根据lammps中Reax/c模块的源码,图中红色的是被忽略的参数,其他颜色记不太清楚了……

而且,元素参数还被分成了四类,根据元素参数第四行倒数第三个和倒数第一个(是否Inner-wall),第二行第二个(是否Shielding vdWaals)

希望能够抛砖引玉


评分 Rate

参与人数
Participants 4
eV +14 收起 理由
Reason
LLG + 3 赞!
Graphite + 5 精品内容
含光君 + 3 谢谢
sobereva + 3

查看全部评分 View all ratings

220

帖子

0

威望

5710

eV
积分
5930

Level 6 (一方通行)

跳跳猪

10#
发表于 Post on 2017-6-3 08:20:15 | 只看该作者 Only view this author
一颗赛艇 发表于 2017-6-3 02:07
发一下自己对力场参数的研究。时间过去太长,有一些也记得不太清楚了。

灰色的在拟合中应该是不变的,根 ...

您好,请问您知道“switch off some bond parameters during the equilibration simulation with ReaxFF”这个是怎么实现的么?
流年似水,浮生如梦。

230

帖子

0

威望

2653

eV
积分
2883

Level 5 (御坂)

11#
发表于 Post on 2017-6-3 13:46:40 | 只看该作者 Only view this author
Mikasa 发表于 2017-6-3 08:20
您好,请问您知道“switch off some bond parameters during the equilibration simulation with ReaxFF ...

这个就不太清楚了。没有研究过。你上谷歌搜搜看看

220

帖子

0

威望

5710

eV
积分
5930

Level 6 (一方通行)

跳跳猪

12#
发表于 Post on 2017-6-3 13:53:02 | 只看该作者 Only view this author
一颗赛艇 发表于 2017-6-3 13:46
这个就不太清楚了。没有研究过。你上谷歌搜搜看看

好的,谢谢~
还没有Google出来
流年似水,浮生如梦。

222

帖子

0

威望

4751

eV
积分
4973

Level 6 (一方通行)

13#
发表于 Post on 2017-6-3 22:39:35 | 只看该作者 Only view this author
一颗赛艇 发表于 2017-6-3 01:39
好像不支持ReaxFF啊

我看的是参数拟合的方法

256

帖子

0

威望

1124

eV
积分
1380

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2017-6-5 09:44:16 | 只看该作者 Only view this author
一颗赛艇 发表于 2017-6-3 02:07
发一下自己对力场参数的研究。时间过去太长,有一些也记得不太清楚了。

灰色的在拟合中应该是不变的,根 ...

谢谢

210

帖子

0

威望

1691

eV
积分
1901

Level 5 (御坂)

15#
发表于 Post on 2017-10-25 22:45:31 | 只看该作者 Only view this author
楼主,您好!我近来想做聚合物热解,但分子动力学软件只会MS里的Forcite,我看了一些文章,他们都是用Lammps里的ReadFF力场,我想着MS不是能拟合力场么,那么ReadFF力场能否直接使用?
越努力越幸运!

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 03:16 , Processed in 0.175054 second(s), 25 queries , Gzip On.

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