计算化学公社

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

[综合交流] ReaxFF反应力场参数优化时如何保证合理的反应温度?

[复制链接 Copy URL]

6

帖子

0

威望

93

eV
积分
99

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
使用MCFF进行ReaxFF反应力场优化,训练集中包括Geo和Eng。对多个随机初始力场样本进行优化后得到了多个误差相近的最优力场,每个力场参数中的sigma键解离能也差别不大,但分别应用这些最优力场进行MD模拟发现,不同力场中相同分子的分解温度差异巨大。
请问在ReaxFF力场中主要是什么参数影响分解温度?在力场优化中如何才能保证得到合理的反应温度?
谢谢!

128

帖子

0

威望

2556

eV
积分
2684

Level 5 (御坂)

2#
发表于 Post on 2020-12-16 13:09:02 | 只看该作者 Only view this author
MCFF是AMS上带的吗?

561

帖子

0

威望

3410

eV
积分
3971

Level 5 (御坂)

3#
发表于 Post on 2020-12-16 14:04:32 | 只看该作者 Only view this author
力场是在拟合势能面。你现在的情况说明这个拟合任务可以收敛到不同的解上,并且相互之间不能统一。

增加训练集可破。策略可以参照https://github.com/deepmodeling/dpgen

6

帖子

0

威望

93

eV
积分
99

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2020-12-16 14:53:56 | 只看该作者 Only view this author
fineren 发表于 2020-12-16 13:09
MCFF是AMS上带的吗?

是的

6

帖子

0

威望

93

eV
积分
99

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2020-12-16 15:33:45 | 只看该作者 Only view this author
k64_cc 发表于 2020-12-16 14:04
力场是在拟合势能面。你现在的情况说明这个拟合任务可以收敛到不同的解上,并且相互之间不能统一。

增加 ...

非常感谢!您说的意思我大致理解。但是我对dpgen不熟悉,看了一下说明还是不太明白。您可以稍详细解释一下的增加训练集的具体思路吗?
我的想法是,比如优化CO2和H2O的力场,需要在训练集中加入CO2和H2O各自的参数,此外还需要加入CO2和H2O同时存在时候的相对能量参数,用来把两种分子的势能面限制到相似的水平。请问这样理解对吗?

561

帖子

0

威望

3410

eV
积分
3971

Level 5 (御坂)

6#
发表于 Post on 2020-12-16 18:01:59 | 只看该作者 Only view this author
本帖最后由 k64_cc 于 2020-12-16 18:03 编辑
liuheng268 发表于 2020-12-16 15:33
非常感谢!您说的意思我大致理解。但是我对dpgen不熟悉,看了一下说明还是不太明白。您可以稍详细解释一 ...

他那个思路其实挺简单的。你先并行训练几个模型,用每个模型再跑一些构象,然后算新构象们在不同模型下的能量,选出最不一致的一批新构象加入训练集。
你都想算分子间相互作用了,训练集里肯定得放包含了它的数据啊……

128

帖子

0

威望

2556

eV
积分
2684

Level 5 (御坂)

7#
发表于 Post on 2020-12-18 08:42:43 | 只看该作者 Only view this author
k64_cc 发表于 2020-12-16 14:04
力场是在拟合势能面。你现在的情况说明这个拟合任务可以收敛到不同的解上,并且相互之间不能统一。

增加 ...

是你们课题组写的吗?有没有中文教程或者培训啊?我要做一类物质的反应动力学,但是没有相应的力场参数

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 22:54 , Processed in 0.695548 second(s), 20 queries , Gzip On.

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