| 查看: 1951 | 回复: 5 | |||
ruthxu金虫 (正式写手)
|
[求助]
Gromacs 中 PME 的用法
|
|
在看文献时，有一段关于cutoff的设定的内容，如下： For Lennard-Jones interactions, a cutoff at 1.4 nm was applied, while electrostatic interactions were treated with the particle mesh Ewald (PME) method and a realspace cutoff of 0.9 nm. A larger than usual Lennard-Jones cutoff was used in order to increase the accuracy of the van der Waals interactions. Electrostatic interactions were treated with the PME method, which does not introduce artificial ordering like cutoff methods. The real-space cutoff of 0.9 nm used with the PME method is merely a numerical device in order to separate the direct- and reciprocal-space sums. Long-range dispersion corrections for energy and pressure were applied. 通过上述内容，在mdp文件中如何具体的设置呢？ 从这个问题中发现，在使用PME 时对于 rlist， rcoulomb， fourierspacing 的设定很困惑，希望高手可以解答一下！ 有相同问题的，可以一起讨论一下！ |
回复此楼» 收录本帖的淘帖专辑推荐
 gromacs |
» 猜你喜欢
求调剂一志愿武汉理工大学材料工程（085601）
已经有5人回复
-
上海电力大学材料防护与新材料重点实验室招收调剂研究生（材料、化学、电化学，环境）
已经有3人回复
-
急发核心期刊论文
已经有5人回复
-
291 求调剂
已经有5人回复
-
298-一志愿中国农业大学-求调剂
已经有11人回复
-
一志愿南京理工大学085701资源与环境302分求调剂
已经有7人回复
-
316求调剂
已经有7人回复
-
材料与化工考研调剂
已经有4人回复
-
一志愿重庆大学085700资源与环境，总分308求调剂
已经有7人回复
-
一志愿211 初试270分 求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用gromacs做分子动力学模拟，如何把底物分子固定住啊
已经有4人回复
- NAMD软件在代码层面上是如何计算蛋白质势能的？
已经有8人回复
- gromacs中MM-GBSA计算
已经有7人回复
- gromacs中如何使用Buckingham势能参数
已经有13人回复
- 求助：Gromacs的g_energy另外一种用法 的具体步骤
已经有7人回复
- 【求助】用NAMD 做NPT模拟遇到下面的问题
已经有10人回复
fcenf
木虫 (初入文坛)
- 模拟EPI: 1
- 应助: 0 (幼儿园)
- 金币: 3032.2
- 红花: 1
- 帖子: 29
- 在线: 93.7小时
- 虫号: 740553
- 注册: 2009-04-05
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ruthxu(金币+10): 很感谢！ 我的设定也是这样的，可师兄认为 rlist=rcoulom=rvdw=1.4 而fourierspacing=0.9，我觉得不对可又讲不出哪里不对。 2011-07-15 08:17:42
zh1987hs(金币+5): 谢谢 2011-07-15 22:56:37
zh1987hs(金币+5): 谢谢 2011-07-15 22:57:06
ruthxu(金币+10): 很感谢！ 我的设定也是这样的，可师兄认为 rlist=rcoulom=rvdw=1.4 而fourierspacing=0.9，我觉得不对可又讲不出哪里不对。 2011-07-15 08:17:42
zh1987hs(金币+5): 谢谢 2011-07-15 22:56:37
zh1987hs(金币+5): 谢谢 2011-07-15 22:57:06
|
按文献里说的，我感觉mdp文件应该是这样吧。 rlist=0.9 vdwtype=cut-off rvdw=1.4 DispCorr=EnerPres coulombtype=PME rcoulomb=0.9 对于PME方法，应该rlist=rcoulomb。在可选的范围内，rlist越大，计算结果越精确，但过大的rlist值会影响计算效率的。主要看你的体系总体势能函数随距离变化的衰减速度，衰减快可以取小些，衰减慢就要适当取大些。 至于fourierspacing，我也说不太好，一般情况下用默认值。感觉其他条件不变，fourierspacing越小，PME方法的结果应该越精确，但过小就会影响计算的效率。0.1～0.12的设定应该是可以接受的。另为，手册上说同时增大或减小rcoulomb和fourierspacing，静电势的计算精度不变的。 |
2楼2011-07-15 00:42:32
fcenf
木虫 (初入文坛)
- 模拟EPI: 1
- 应助: 0 (幼儿园)
- 金币: 3032.2
- 红花: 1
- 帖子: 29
- 在线: 93.7小时
- 虫号: 740553
- 注册: 2009-04-05
- 性别: GG
- 专业: 理论和计算化学
3楼2011-07-15 00:54:17
ruthxu
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1109.2
- 散金: 62
- 红花: 2
- 帖子: 319
- 在线: 138小时
- 虫号: 1159022
- 注册: 2010-11-30
- 性别: MM
- 专业: 高分子科学
4楼2011-07-15 08:21:19
fcenf
木虫 (初入文坛)
- 模拟EPI: 1
- 应助: 0 (幼儿园)
- 金币: 3032.2
- 红花: 1
- 帖子: 29
- 在线: 93.7小时
- 虫号: 740553
- 注册: 2009-04-05
- 性别: GG
- 专业: 理论和计算化学
5楼2011-07-15 16:36:30
ruthxu
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1109.2
- 散金: 62
- 红花: 2
- 帖子: 319
- 在线: 138小时
- 虫号: 1159022
- 注册: 2010-11-30
- 性别: MM
- 专业: 高分子科学
6楼2011-07-17 09:53:03
