计算化学公社

标题: N2/O2的TraPPE模型参数如何加到itp文件中 [打印本页]

作者
Author: 五月雨 时间: 2022-7-4 22:03
标题: N2/O2的TraPPE模型参数如何加到itp文件中
本帖最后由五月雨于 2022-7-5 10:51 编辑

各位老师好，我打算使用GAGG力场模拟单乙醇胺在气相中的性质，因为主要考察胺的性质，所以就使用了GAFF力场。之前想着N2/O2也通过sobtop在GAFF或者基于Hessian矩阵计算分配力场参数。但是看到sob老师在以前的帖子上说N2/O2这种重要的小分子使用TraPPE模型参数比较好。但是TraPPE的itp文件不会编写。

1. 我只会根据第一张图片里的内容编写itp文件，给的TraPPE的参数实在看不懂。范德华参数是不是ε/kB [K]乘上理想气体摩尔常数8.314472J/（mol·K）？那力常数呢？

2. 还想请问老师，TraPPE的模型能不能和GAFF联用呢？

(, 下载次数 Times of downloads: 27)
TraPPE：
(, 下载次数 Times of downloads: 25)
(, 下载次数 Times of downloads: 17)

作者
Author: sobereva 时间: 2022-7-5 09:18
参考我写的EPM2的CO2的拓扑文件
(, 下载次数 Times of downloads: 93)

如果体系的某些部分不得不用GAFF，那就联用

作者
Author: 五月雨 时间: 2022-7-5 10:06
本帖最后由五月雨于 2022-7-5 19:44 编辑

sobereva 发表于 2022-7-5 09:18
参考我写的EPM2的CO2的拓扑文件

谢谢sob老师的回复。还想请问老师，TraPPE的参数是不是由于是刚性模型，所以没有键和角的力常数？我在编写TraPPE参数的itp文件时，把键和角都加一个约束是不是就可以了？这是我编写的，老师您看有没有问题。
(, 下载次数 Times of downloads: 39)

作者
Author: sobereva 时间: 2022-7-6 07:28

五月雨发表于 2022-7-5 10:06
谢谢sob老师的回复。还想请问老师，TraPPE的参数是不是由于是刚性模型，所以没有键和角的力常数？我在编 ...

是
试一下，跑起来如果形态能维持就行

作者
Author: 高处裹棉被 时间: 2022-11-1 14:05

五月雨发表于 2022-7-5 10:06
谢谢sob老师的回复。还想请问老师，TraPPE的参数是不是由于是刚性模型，所以没有键和角的力常数？我在编 ...

作者大大您好，我看了您的N2拓扑文件，

想请问N2只有两个原子，中间那个是不是应该设置成虚拟原子？不过还想请问虚拟原子该怎么设置？

谢谢您的解答

作者
Author: 五月雨 时间: 2022-11-3 08:46

高处裹棉被发表于 2022-11-1 14:05
作者大大您好，我看了您的N2拓扑文件，

想请问N2只有两个原子，中间那个是不是应该设置成虚拟原子？不 ...

我是参考四点水模型设置的，gromacs官方教程里面有对几种虚原子的介绍。我也不知道我设置的对不对，最后因为gro文件原子不对应问题解决起来很麻烦，我就没再用这个模型

作者
Author: 高处裹棉被 时间: 2022-11-3 10:28

五月雨发表于 2022-11-3 08:46
我是参考四点水模型设置的，gromacs官方教程里面有对几种虚原子的介绍。我也不知道我设置的对不对，最后 ...

非常感谢，想请问您最后对氮气用的模型文件是怎样设置的，

如果可以您能提供一下具体内容吗，我想试一试

作者
Author: 五月雨 时间: 2022-11-6 11:34

高处裹棉被发表于 2022-11-3 10:28
非常感谢，想请问您最后对氮气用的模型文件是怎样设置的，

如果可以您能提供一下具体内容吗，我想试一 ...

你试一试看看可不可行

作者
Author: 喵星大佬 时间: 2023-4-4 14:53

五月雨发表于 2022-11-6 11:34
你试一试看看可不可行

你这里面根本就没有定义虚拟点的类型
我根据TraPPE-Small力场参数，做了一下空气的模拟(78%氮气+22%氧气)

gmx insert-molecules用两个pdb文件，
(, 下载次数 Times of downloads: 68)
(, 下载次数 Times of downloads: 38)

两个分子的itp文件，
(, 下载次数 Times of downloads: 86)
(, 下载次数 Times of downloads: 56)

最终的结构文件和拓扑文件
(, 下载次数 Times of downloads: 28)
(, 下载次数 Times of downloads: 57)

和常温常压模拟的mdp文件在附件里
(, 下载次数 Times of downloads: 47)

作者
Author: 喵星大佬 时间: 2023-4-5 16:19
包含该力场完整参数的gromacs格式力场包(包含结构的pdb文件方便建模)已上传

见这个帖子

http://bbs.keinsci.com/thread-36352-1-1.html

作者
Author: margin 时间: 2024-6-11 09:22

五月雨发表于 2022-7-5 10:06
谢谢sob老师的回复。还想请问老师，TraPPE的参数是不是由于是刚性模型，所以没有键和角的力常数？我在编 ...

请问作者，你联用GAFF力场和TRappe力场的时候，top文件应该怎么写。能稍微介绍一下吗?

作者
Author: 五月雨 时间: 2024-6-30 12:06

margin 发表于 2024-6-11 09:22
请问作者，你联用GAFF力场和TRappe力场的时候，top文件应该怎么写。能稍微介绍一下吗?

N2和O2都是小分子，没有1-4相互作用，所以top我是根据GAFF力场写的
[ defaults ]
; nbfunc comb-rule gen-pairs fudgeLJ fudgeQQ
1 2 yes 0.5 0.8333

作者
Author: sdf 时间: 2024-11-20 20:29

margin 发表于 2024-6-11 09:22
请问作者，你联用GAFF力场和TRappe力场的时候，top文件应该怎么写。能稍微介绍一下吗?

请问作者，我现在和你有同款疑问吗？请问你最后咋解决的，方便提供一下top文件学习一下吗？

作者
Author: raku 时间: 2024-12-5 10:44

喵星大佬发表于 2023-4-4 14:53
你这里面根本就没有定义虚拟点的类型
我根据TraPPE-Small力场参数，做了一下空气的模拟(78%氮气+22%氧气 ...

请教下大佬，在O2_TraPPE和N2_TraPPE中都存在virtual site 名字为MW。发现在您的air.top中 [ atomtypes ]下只有一个MW，我没太理解这点，可以解释解释吗？按照我的理解应该包含O2和N2中的两个virtual site.谢谢啦！

作者
Author: raku 时间: 2024-12-5 10:45

sdf 发表于 2024-11-20 20:29
请问作者，我现在和你有同款疑问吗？请问你最后咋解决的，方便提供一下top文件学习一下吗？

12楼就是呀

作者
Author: 喵星大佬 时间: 2024-12-5 15:18

raku 发表于 2024-12-5 10:44
请教下大佬，在O2_TraPPE和N2_TraPPE中都存在virtual site 名字为MW。发现在您的air.top中 [ atomtypes ] ...

都是没有vdW作用的虚拟点，只有电荷不同，只需要定义一次，电荷在moleculetype里定义就行

作者
Author: raku 时间: 2024-12-5 16:34

喵星大佬发表于 2024-12-5 15:18
都是没有vdW作用的虚拟点，只有电荷不同，只需要定义一次，电荷在moleculetype里定义就行

明白了，谢谢！

还有个问题是在您的pdb中Mw为小写，itp中MW为大写，我试算时候说不对，我修改了就好了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)