| 5 | 1/1 | 返回列表 |
| 查看: 3745 | 回復(fù): 13 | ||||
| 本帖產(chǎn)生 1 個(gè) 模擬EPI ,點(diǎn)擊這里進(jìn)行查看 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
zhangmao511鐵蟲 (小有名氣)
|
[求助]
Gromacs里關(guān)于截?cái)喟霃降囊粋(gè)問題,謝謝。 已有1人參與
|
|||
|
在讀Gromacs手冊(cè)的時(shí)候,有2段話描述了截?cái)喟霃。?duì)這2段和rlist,rcoulomb,rvdw比較疑惑,想請(qǐng)教各位大蝦一下。 1. GROMACS employs a pair list that contains those particle pairs for which non-bonded forces must be calculated. The pair list contains atoms i, a displacement vector for atom i, and all particles j that are within rlist of this particular image of atom i. The list is updated every nstlist steps, where nstlist is typically 10. 2. In the neighbor list, all interaction pairs that fall within rlist are stored. Furthermore, the interactions between pairs that do not fall within rlist but do fall within max(rcoulomb,rvdw) are computed during NS. The forces and energy are stored separately and added to short-range forces at every time step between successive NS. If rlist = max(rcoulomb,rvdw), no forces are evaluated during neighbor list generation. 我對(duì)上面部分的理解是,兩個(gè)原子間的距離小于rlist的時(shí)候才計(jì)算這些原子之間的短程相互作用。相互作用的原子對(duì),如果距離在rlist之內(nèi)的就被存儲(chǔ)下來,沒有落到rlist內(nèi)但是落到max(rcoulomb,rvdw)內(nèi)的,會(huì)在鄰近搜索期間被計(jì)算。 那么rlist和rcoulomb, rvdw具體是什么關(guān)系呢??不知道能不能解釋的詳細(xì)點(diǎn),謝謝! 還有最后一句話,“If rlist = max(rcoulomb,rvdw), no forces are evaluated during neighbor list generation”,不怎么理解,請(qǐng)大家指教,謝謝!! |
MDs-Gromacs | gromacs |
鐵蟲 (小有名氣)
|
1. 雙程截?cái),就是將力的范圍分為short range + long range兩部分進(jìn)行計(jì)算吧? 你說當(dāng)rw=rc>rlist時(shí),在rw和rlist間的距離是rlistlong,mdp參數(shù)文件里有一個(gè)參數(shù)也叫rlistlong,你指的是這個(gè)?rw和rlist間的距離難道不就是2者的差(應(yīng)該算long range的部分吧),如果在rw和rlist間力也需要每步計(jì)算,那么nslist(刷新頻率)的作用體現(xiàn)在哪呢……??謝謝 2. 我看有的文獻(xiàn)里也會(huì)有rw=rc<rlist,那這種情況呢??用PME算法的時(shí)候,是不是選用rlist=rw=rc比較好?謝謝 ![]() ![]() |
銅蟲 (小有名氣)
|
rlist是在短程截?cái)喟虢飪?nèi) 所需要計(jì)算的非鍵相互作用。(某些力場(chǎng)排除1-4相互作用,或則按照0.5加成計(jì)算) rcoulmob是短程庫(kù)侖力截?cái)喟霃剑诖酥畠?nèi)的計(jì)算靜電力。在此之外的按照PME實(shí)空間和虛空間算 rvwd是短程范德華力的截?cái)喟霃,在此之?nèi)的原子對(duì)計(jì)算范德華力,在此之外的力為0 或范德華力矯正。 以上為3種不同的力的計(jì)算截?cái)唷?br /> 因?yàn)樯婕暗洁徑斜硭阉鱊S計(jì)算。 rlist=rcoulmob=rvwd 計(jì)算快,精度較高 rlist<rc/rw 計(jì)算資源高,精度更高 第三種,精度有限。因此第一種用的比較多。 |
銅蟲 (小有名氣)
鐵蟲 (小有名氣)
|
感謝你的講解。 ![]() ![]() 1. 你說rlist是是非共價(jià)鍵相互作用的截?cái)喟霃。而非共價(jià)鍵就包括了靜電相互作用和范德華力 ,那我能理解rlist包含了rcoulmob和rvwd嗎?那既然是這樣,為什么要把rcoulmob和rvwd這個(gè)2個(gè)力單獨(dú)提出來計(jì)算呢??謝謝 2.我是看到有不少文獻(xiàn)的method里,用的是rlist < rcoulmob = rw,而Gromacs官方的例子用的卻是rlist=rcoulmob=rvwd, 當(dāng)時(shí)我還很疑惑呢!!我這里想問一下,為什么當(dāng)rlist < rcoulmob = rw的時(shí)候計(jì)算資源高,精度更高??二rlist=rcoulmob=rvwd 計(jì)算快,精度較高?? 感謝。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | Mqqqqqq 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +6 | 申子申申 2026-03-19 | 11/550 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過,有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|
[考研] 265求調(diào)劑 +8 | 梁梁校校 2026-03-17 | 8/400 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +3 | sbdksD 2026-03-19 | 3/150 |
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|