| 查看: 1954 | 回复: 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 |
» 猜你喜欢
085600材料与化工调剂
已经有10人回复
-
279分求调剂 一志愿211
已经有18人回复
-
0703化学求调剂
已经有5人回复
-
材料292调剂
已经有3人回复
-
一志愿陕师大生物学071000，298分，求调剂
已经有3人回复
-
306求0703调剂一志愿华中师范
已经有7人回复
-
333求调剂
已经有3人回复
-
335求调剂
已经有4人回复
-
336求调剂
已经有4人回复
-
303求调剂
已经有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
- 专业: 理论和计算化学
5楼2011-07-15 16:36:30
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
