| 查看: 2347 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
abin007银虫 (初入文坛)
|
[交流]
【讨论】分子模拟中力场参数的调节讨论已有3人参与
|
||
|
经验力场在分子模拟中的重要性不用我多说了，但困惑我们的是，当我们的结构在现有力场中找不到参数时，或者是即使找到了力场参数，模拟结果并不准确，这时候调整立场参数成了不可避免的问题。 问题的复杂性在于，这么多的力场参数，从哪儿下手？怎么调整？而且有时候是牵一发而动全身啊。我最近涉及到了这一个问题，就我的经验跟大家交流交流，希望得到大家的响应。 就我的例子而言，如果想要修改二面角参数，首先第一步是取消二面角中1-4原子间的静电力和范德华力，可以这样比方：1-4原子间的非键力像是一个大小可变的盒子的尺寸，如果盒子里面的系统始终小于盒子，那么该盒子大小的变化不会带来多大的系统变化；一旦盒子里面的系统要主动变大时（就像我们调整二面角的参数时），这时该盒子的尺寸就会发挥较明显的限制作用了。因此，如果我们希望调整二面角的力场参数并能看到明显地效果时，应该取消1-4原子间的非键力。不知我这个比喻是否恰当。 接下来，我的办法就比较笨了，逐一调整二面角的力参数，来比较最终结果，有时候会添加或删除一些二面角参数。我的原则是，二面角的个数从少到多，不要刚开始就一股脑儿地把所有的二面角都用上。 以上就是我跟大家分享的浅薄经验，希望能跟大家继续交流这个话题。 |
回复此楼» 猜你喜欢
期末给学生划重点都是什么话术啊
已经有4人回复
-
事业单位还是大学好？
已经有19人回复
-
发现督导在听课，需要跟TA招呼示意吗？有同学看手机、课件有疏漏算大问题吗？
已经有5人回复
-
前同事拿我以前未中的国自然标书申请今年的国自然，现在基金委发来调查函
已经有11人回复
-
杭州国企和浙江高校如何选择？
已经有22人回复
-
小孩7岁，上一年级，不自信，学习较差
已经有36人回复
-
想问一下有没有需要文章但是没时间带学生的研究生导师。
已经有19人回复
-
颜宁:基础研究应顶天立地,立地即从需求端作为问题的出发点,立项背景就是面向真实需求
已经有3人回复
-
催稿后秒拒
已经有6人回复
-
298求调剂
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- AutoDock参数设置问题
已经有6人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【讨论】蛋白质上的非正常氨基酸残基参数如何确定
已经有3人回复
- 【讨论】请教大家几个问题，一直理解不透
已经有4人回复
- 【讨论】分子模拟与工作
已经有7人回复
- 【专家会诊】bay 讨论namd,gmx处理水溶液,蛋白等MD 模拟;讨论水及水通道
已经有50人回复
- 【小卒个人文集】写给初学者:小议形式电荷对量化初学者的误导
已经有27人回复
- 【讨论】Tersoff potential 的参数怎么用Gulp实现
已经有6人回复
- 【求助】MS-sorption中力场编辑的一个基础问题
已经有8人回复
- 【求助】请问有拟合分子力场参数的FORTRAN程序吗？
已经有13人回复
- 【讨论】铜和水之间的LJ势函数的参数是多少？
已经有19人回复
- 【sobereva个人文集】转载一篇原创精品文章：过渡态、反应路径的计算方法及相关问题
已经有42人回复
043114076
木虫 (正式写手)
- 应助: 49 (小学生)
- 金币: 3364.9
- 散金: 52
- 红花: 23
- 帖子: 788
- 在线: 369.5小时
- 虫号: 570592
- 注册: 2008-06-07
- 性别: GG
- 专业: 理论和计算化学
2楼2011-02-14 10:26:04
043114076
木虫 (正式写手)
- 应助: 49 (小学生)
- 金币: 3364.9
- 散金: 52
- 红花: 23
- 帖子: 788
- 在线: 369.5小时
- 虫号: 570592
- 注册: 2008-06-07
- 性别: GG
- 专业: 理论和计算化学
4楼2011-02-15 09:21:12
043114076
木虫 (正式写手)
- 应助: 49 (小学生)
- 金币: 3364.9
- 散金: 52
- 红花: 23
- 帖子: 788
- 在线: 369.5小时
- 虫号: 570592
- 注册: 2008-06-07
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+7): 谢谢 2011-02-16 13:05:09
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+7): 谢谢 2011-02-16 13:05:09
|
In OPLS-AA, such nonbonded intramolecular interactions are described by the same parameters used in intermolecular interactions, only scaled by a factor of 0.5. 这句话来自于文献：Molecular Force Field for Ionic Liquids Composed of Triflate or Bistriflylimide Anions Jose´ N. Canongia Lopes, J. Phys. Chem. B 2004, 108, 16893-16898 我看了OPLS拟合的文献，好像也是有这样的定义。 当然不同的力场，定义1-4相互作用也不同。以下是以“charmm 1-4 interaction”在google中搜到的结果，也许对你有用： LJ and electrostatic interactions between 1-4 neighbors are calculated, but can be modified compared to more distant pairs (different LJ interactions, scaling of electrostatic interactions). The details depend on the specific force field used. E.g., in the older polar hydrogen force field ("param19"  , electrostatic interactions and forces between 1-4 neighbors were scaled by 0.4, whereas in the current all-atom force field(s) ("param22", "param27", electrostatic energies and forces are used as is. This scaling factor is controlled by one of the many options to the ENERgy command, E14Fac. Thus, for calculations with the old, polar-hydrogen "param19" force field, E14Fac should be set to 0.4; for calculations with the current "param22/27" force field it should be set to 1.0. 来源：http://www.charmmtutorial.org/index.php/The_Energy_Function 所以我觉得你应该好好查查charmm的使用说明。 |
7楼2011-02-16 09:06:19
10