计算化学公社

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

[GROMACS] 使用文献中有机物分子参数求助

[复制链接 Copy URL]

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

社长好,我想用一篇文献中的有机物分子的参数,有一些问题请教
1.他的力场形式如下图一。这个形式是CHARMM力场吗,文章说的是用charmm的原子类型。
2.我想在体系中利用Merz的离子参数加入离子,按照文章中这个势能的非键作用项用对Merz参数算出的epsilon与sigma进行变化吗(文献中的静电作用项比平时所见的静电项的分母中少了4πε系数,所以不确定范德华参数改如何设置)
3.CHARMM力场中含有Urey-Bradley项,但是用www.swissparam.ch自动生成的itp中没有Urey-Bradley项,如何加上呢。因为文献中的有Kub S0参数,网站自动生成的itp的[angles]项中
[angles];
ai aj ak fu th0 kth ub0 kub th0 kth ub0 kub   重复两次th0 kth但是没有ub0 kub

4.由于文章使用的是NAMD,因此我要把他的力场转化一下,他在拓扑结构部分只给出了ATOM和BOND但是没有给angles和dihedrals,在力场参数部分又给了这两部分的参数,自己写itp的时候需要这两部分吗
我没有找到所给参数中没有给14作用的[paris],是说明他没有设置吗。
另外他给的NONBONDED参数中不太理解他的意思:
om ignored  eq Rmin/2    ignored  eps,1-4      Rmin/2,1-4分别是什么意思(如图二)

原文件见附件。问题较多,辛苦社长









{1GLGP@J2XB{YIP@ROZR8QO.png (10.6 KB, 下载次数 Times of downloads: 23)

力场形式

力场形式

2.png (24.88 KB, 下载次数 Times of downloads: 24)

图2

图2

力场.docx

44.56 KB, 下载次数 Times of downloads: 2

原文件力场部分

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2021-2-6 22:29:51 | 只看该作者 Only view this author
1 是
2 如果式子里没4π,那应该是等效体现在epsilon里了
3 手动加上
4 我对NAMD、CHARMM派系的东西不熟。NAMD网站上应该有力场参数文件格式的解释
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 02:03 , Processed in 0.348120 second(s), 23 queries , Gzip On.

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