| 查看: 1833 | 回复: 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 |
» 猜你喜欢
AI论文写作工具：是科研加速器还是学术作弊器？
已经有5人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
用gromacs做分子动力学模拟，如何把底物分子固定住啊
已经有4人回复- NAMD软件在代码层面上是如何计算蛋白质势能的？
已经有8人回复
- gromacs中MM-GBSA计算
已经有7人回复
- gromacs中如何使用Buckingham势能参数
已经有13人回复
- 求助：Gromacs的g_energy另外一种用法 的具体步骤
已经有7人回复
- 【求助】用NAMD 做NPT模拟遇到下面的问题
已经有10人回复
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
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
fcenf
木虫 (初入文坛)
- 模拟EPI: 1
- 应助: 0 (幼儿园)
- 金币: 3032.2
- 红花: 1
- 帖子: 29
- 在线: 93.7小时
- 虫号: 740553
- 注册: 2009-04-05
- 性别: GG
- 专业: 理论和计算化学
5楼2011-07-15 16:36:30