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

金蟲 (正式寫手)

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

木蟲 (初入文壇)

木蟲 (初入文壇)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +8 | 何潤采123 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[基金申請(qǐng)]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|