| 5 | 1/1 | 返回列表 |
| 查看: 1945 | 回復: 5 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
ruthxu金蟲 (正式寫手)
|
[求助]
Gromacs 中 PME 的用法
|
||
|
在看文獻時,有一段關于cutoff的設定的內(nèi)容,如下: 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. 通過上述內(nèi)容,在mdp文件中如何具體的設置呢? 從這個問題中發(fā)現(xiàn),在使用PME 時對于 rlist, rcoulomb, fourierspacing 的設定很困惑,希望高手可以解答一下! 有相同問題的,可以一起討論一下! |
gromacs |

木蟲 (初入文壇)

木蟲 (初入文壇)
|
按文獻里說的,我感覺mdp文件應該是這樣吧。 rlist=0.9 vdwtype=cut-off rvdw=1.4 DispCorr=EnerPres coulombtype=PME rcoulomb=0.9 對于PME方法,應該rlist=rcoulomb。在可選的范圍內(nèi),rlist越大,計算結(jié)果越精確,但過大的rlist值會影響計算效率的。主要看你的體系總體勢能函數(shù)隨距離變化的衰減速度,衰減快可以取小些,衰減慢就要適當取大些。 至于fourierspacing,我也說不太好,一般情況下用默認值。感覺其他條件不變,fourierspacing越小,PME方法的結(jié)果應該越精確,但過小就會影響計算的效率。0.1~0.12的設定應該是可以接受的。另為,手冊上說同時增大或減小rcoulomb和fourierspacing,靜電勢的計算精度不變的。 |

木蟲 (初入文壇)

金蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 初試 317 +3 | 半拉月丙 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 299求調(diào)劑 +3 | 某某某某位 2026-03-21 | 3/150 |
|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
[考研] 機械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 307求調(diào)劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
[考研] A區(qū)線材料學調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 289求調(diào)劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學666 2026-03-16 | 7/350 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|