计算化学公社

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

[GROMACS] 关于ClayFF力场中参数取值的疑问

[复制链接 Copy URL]

14

帖子

0

威望

107

eV
积分
121

Level 2 能力者

本帖最后由 KnightWE 于 2019-7-26 00:30 编辑

ClayFF力场常用于黏土及水分子/阴阳离子的模拟,最早于2004年发布于JPCB(https://doi.org/10.1021/jp0363287),引用千次有余。但网上关于该力场的讨论似乎较少,特别是其在GROMACS中的应用。楼主初学MD,试图基于论文中给出的参数编写GROMACS适用的力场文件,但在此过程中发现一些参数取值略奇怪,:

ClayFF力场中,利用简谐振动模型描述羟基OH键的伸缩振动和弯曲振动(见图中右半公式),但将文中Force Constant转为GROMACS中所用单位时,我发现其值与其他力场相差甚大(仅有一半左右),这令我让我感到非常费解。

随后我查到了另一篇应用ClayFF进行计算的文献(doi:10.1016/j.gca.2011.01.026),发现列于SI的力场参数也与原文不符,取值为原文两倍左右(接近常见FF所用参数)。

能否请坛里前辈帮忙看看,正好2倍的差值是我对力场参数的理解在哪里出问题了么?作为一篇被广泛引用的经典论文,小生觉得问题出在自己身上的概率更大些。

望不吝赐教,感谢!




4

帖子

0

威望

65

eV
积分
69

Level 2 能力者

19#
发表于 Post on 2025-3-29 10:31:40 | 只看该作者 Only view this author
阿di达思 发表于 2019-7-27 15:54
您好,请问您现在是否成功使用过clayff在gromacs进行计算?能否共享一下力场文件。

大佬您好,请问您现在是否成功使用过clayff在gromacs进行计算?请问可以共享一下力场文件吗。

1481

帖子

0

威望

4554

eV
积分
6035

Level 6 (一方通行)

18#
发表于 Post on 2023-10-16 18:45:46 | 只看该作者 Only view this author
又菜又爱玩

14

帖子

0

威望

107

eV
积分
121

Level 2 能力者

17#
 楼主 Author| 发表于 Post on 2023-7-18 17:08:40 | 只看该作者 Only view this author
chieko 发表于 2023-1-10 23:21
Advances in Clayff Molecular Simulation of Layered and Nanoporous Materials and Their Aqueous Interf ...

多谢!

33

帖子

0

威望

726

eV
积分
759

Level 4 (黑子)

16#
发表于 Post on 2023-1-10 23:21:49 | 只看该作者 Only view this author
Advances in Clayff Molecular Simulation of Layered and Nanoporous Materials and Their Aqueous Interfaces, J. Phys. Chem. C 2021, 125, 32, 17573–17589
https://pubs.acs.org/doi/10.1021/acs.jpcc.1c04600
关于CLAYFF力场的综述,SI里面有CLAYFF和gmx单位转化的详细过程和gromacs力场文件。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
tjuptz + 5 谢谢

查看全部评分 View all ratings

30

帖子

0

威望

928

eV
积分
958

Level 4 (黑子)

15#
发表于 Post on 2019-9-18 17:01:06 | 只看该作者 Only view this author
KnightWE 发表于 2019-7-30 02:02
原作者的力场文件里species没写全(我看了他的毕业论文,因为他的体系不需要所有参数),你按照你需要的 ...

您好,请问可以分享一下原作者的毕业论文吗?

516

帖子

1

威望

4703

eV
积分
5239

Level 6 (一方通行)

14#
发表于 Post on 2019-7-30 07:09:53 | 只看该作者 Only view this author
KnightWE 发表于 2019-7-30 02:02
原作者的力场文件里species没写全(我看了他的毕业论文,因为他的体系不需要所有参数),你按照你需要的 ...

我说原子类型咋那么少,看着不对劲,我也看看原文去

14

帖子

0

威望

107

eV
积分
121

Level 2 能力者

13#
 楼主 Author| 发表于 Post on 2019-7-30 02:02:25 | 只看该作者 Only view this author
tjuptz 发表于 2019-7-29 19:00
谢谢,期待您的分享~

原作者的力场文件里species没写全(我看了他的毕业论文,因为他的体系不需要所有参数),你按照你需要的体系补充下就可以了。继续交流哈

评分 Rate

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

查看全部评分 View all ratings

516

帖子

1

威望

4703

eV
积分
5239

Level 6 (一方通行)

12#
发表于 Post on 2019-7-29 19:00:20 | 只看该作者 Only view this author
KnightWE 发表于 2019-7-29 18:47
作者似乎是根据他自己的体系做了一定的修正(有个0.89倍的系数),总体来说没问题,我给作者发了邮件还没 ...

谢谢,期待您的分享~

14

帖子

0

威望

107

eV
积分
121

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2019-7-29 18:47:48 | 只看该作者 Only view this author
tjuptz 发表于 2019-7-29 16:34
您好,我在https://github.com/thomasunderwood/ClayFF/issues/1 上看到了您的提问,请问该力场文件跟原 ...

作者似乎是根据他自己的体系做了一定的修正(有个0.89倍的系数),总体来说没问题,我给作者发了邮件还没回。

我也是基于他之前的工作在写力场文件,这应该是我找到的唯一公开的比较有价值的参考了

516

帖子

1

威望

4703

eV
积分
5239

Level 6 (一方通行)

10#
发表于 Post on 2019-7-29 16:34:27 | 只看该作者 Only view this author
KnightWE 发表于 2019-7-28 22:13
抱歉,还在尝试中,还存在一些问题。等弄完之后我会上传到github

您好,我在https://github.com/thomasunderwood/ClayFF/issues/1 上看到了您的提问,请问该力场文件跟原文相比是有问题吗?

14

帖子

0

威望

107

eV
积分
121

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2019-7-28 22:13:17 | 只看该作者 Only view this author
阿di达思 发表于 2019-7-27 15:54
您好,请问您现在是否成功使用过clayff在gromacs进行计算?能否共享一下力场文件。

抱歉,还在尝试中,还存在一些问题。等弄完之后我会上传到github

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
阿di达思 + 5 精品内容

查看全部评分 View all ratings

14

帖子

0

威望

107

eV
积分
121

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2019-7-28 22:12:29 | 只看该作者 Only view this author
ruanyang 发表于 2019-7-26 21:35
可以试试这个力场 https://bionanostructures.com/interface-md/

嗯嗯,等我把ClayFF搞懂之后我再去看下。感谢分享

516

帖子

1

威望

4703

eV
积分
5239

Level 6 (一方通行)

7#
发表于 Post on 2019-7-28 10:33:40 | 只看该作者 Only view this author
同求

30

帖子

0

威望

4529

eV
积分
4559

Level 6 (一方通行)

6#
发表于 Post on 2019-7-27 15:54:24 | 只看该作者 Only view this author
您好,请问您现在是否成功使用过clayff在gromacs进行计算?能否共享一下力场文件。

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

GMT+8, 2025-8-18 10:35 , Processed in 1.126544 second(s), 25 queries , Gzip On.

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