计算化学公社

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

[GROMACS] 非正常二面角参数应该怎么去拟合?

[复制链接 Copy URL]

22

帖子

0

威望

115

eV
积分
137

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
二面角的参数拟合需要先用QM扫描二面角的能量,然后用MM的能量去拟合。但是对于非正常二面角不知道应该怎么去进行能量扫描,所以想请问老师拟合力场非正常二面角参数的一般步骤大概是怎么样的?

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

2#
发表于 Post on 2021-10-15 05:11:22 | 只看该作者 Only view this author
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,避免有任何歧义,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“非正常二面角参数拟合”改了,以后务必注意。

说清楚什么叫非正常二面角,举个例子,各个图
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

22

帖子

0

威望

115

eV
积分
137

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2021-10-15 10:13:28 | 只看该作者 Only view this author

如图,理论上四个原子应该处于同一个平面,但是在MD的过程中j原子会和另外三个原子形成的平面有角度,这时候需要引入improper dihedral项使得在模拟的过程中这四个原子在同一平面

22

帖子

0

威望

115

eV
积分
137

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2021-10-15 10:18:17 | 只看该作者 Only view this author
[ dihedraltypes ]
; 'improper' dihedrals
;      i        j        k        l  func         phi0         kphi
       O      CP1      NH2       CC     2     0.000000   376.560000
这是charmm力场中对酰胺的非正常二面角参数,其他参数都好给,主要是需要拟合kphi

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

5#
发表于 Post on 2021-10-16 02:38:11 | 只看该作者 Only view this author
VincenzoChen 发表于 2021-10-15 10:13
如图,理论上四个原子应该处于同一个平面,但是在MD的过程中j原子会和另外三个原子形成的平面有角度,这 ...

描述问题时一些关键英文术语不要用有广泛共识的中文,说非正常二面角会带来严重歧义,而说improper内行人都明白

improper二面角主要目的是维持平面性,达到目的就够了,对于力常数不需要太讲究,随便猜一个够大的(但又不能太大,否则造成高频振荡)通常就能满足目的。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

6#
发表于 Post on 2021-10-16 02:40:52 | 只看该作者 Only view this author
如果必须明确地柔性扫描phi角,恐怕只能用Gaussian 16的GIC,靠比如n-i-j和k-i-j角基于几何关系定义phi变量并扫它
看下文相关部分
详谈使用Gaussian做势能面扫描
http://sobereva.com/474http://bbs.keinsci.com/thread-12660-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

22

帖子

0

威望

115

eV
积分
137

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2021-10-18 15:38:54 | 只看该作者 Only view this author
sobereva 发表于 2021-10-16 02:40
如果必须明确地柔性扫描phi角,恐怕只能用Gaussian 16的GIC,靠比如n-i-j和k-i-j角基于几何关系定义phi变量 ...

好的,谢谢sob老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 01:13 , Processed in 0.157164 second(s), 20 queries , Gzip On.

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