 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3740 | 回復(fù): 13 | ||||
| 本帖產(chǎn)生 1 個 模擬EPI ，點(diǎn)擊這里進(jìn)行查看 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
zhangmao511鐵蟲 (小有名氣)
|
[求助]
Gromacs里關(guān)于截?cái)喟霃降囊粋�問題，謝謝��！ 已有1人參與
|
|||
|
在讀Gromacs手冊的時候，有2段話描述了截?cái)喟霃�。對這2段和rlist,rcoulomb,rvdw比較疑惑，想請教各位大蝦一下。 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. 我對上面部分的理解是，兩個原子間的距離小于rlist的時候才計(jì)算這些原子之間的短程相互作用。相互作用的原子對，如果距離在rlist之內(nèi)的就被存儲下來，沒有落到rlist內(nèi)但是落到max(rcoulomb,rvdw)內(nèi)的，會在鄰近搜索期間被計(jì)算。 那么rlist和rcoulomb, rvdw具體是什么關(guān)系呢？？不知道能不能解釋的詳細(xì)點(diǎn)，謝謝！ 還有最后一句話，“If rlist = max(rcoulomb,rvdw), no forces are evaluated during neighbor list generation”，不怎么理解，請大家指教，謝謝��！ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 MDs-Gromacs | gromacs |
» 猜你喜歡
265求調(diào)劑
已經(jīng)有8人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有21人回復(fù)
-
295材料求調(diào)劑，一志愿武漢理工085601專碩
已經(jīng)有4人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有9人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
已經(jīng)有7人回復(fù)
-
材料學(xué)碩297已過四六級求調(diào)劑推薦
已經(jīng)有6人回復(fù)
-
295復(fù)試調(diào)劑
已經(jīng)有7人回復(fù)
-
本人考085602 化學(xué)工程 專碩
已經(jīng)有20人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- lammps中LJ勢能下單位問題
已經(jīng)有5人回復(fù)
- 分子動力學(xué)模擬時間的意義？
已經(jīng)有15人回復(fù)
- 關(guān)于MS做分子動力學(xué)的問題
已經(jīng)有8人回復(fù)
- MD模擬軟件選擇
已經(jīng)有3人回復(fù)
- gromacs或者VMD中怎么截取特定的面？
已經(jīng)有8人回復(fù)
- 請問周期性邊界條件什么意思？
已經(jīng)有6人回復(fù)
- 【求助】ReaxFF 力場參數(shù)化問題請教
已經(jīng)有9人回復(fù)
- 【轉(zhuǎn)帖】用Amber+Gaussian做小分子力場
已經(jīng)有16人回復(fù)
- 截?cái)喟霃竭x取的標(biāo)準(zhǔn)是什么？
已經(jīng)有6人回復(fù)
zhangmao511
鐵蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 33
- 紅花: 2
- 帖子: 136
- 在線: 29.7小時
- 蟲號: 2553312
- 注冊: 2013-07-19
- 性別: GG
- 專業(yè): 蛋白質(zhì)組學(xué)
|
感謝你的講解。  1. 你說rlist是是非共價(jià)鍵相互作用的截?cái)喟霃�。而非共價(jià)鍵就包括了靜電相互作用和范德華力 ，那我能理解rlist包含了rcoulmob和rvwd嗎？那既然是這樣，為什么要把rcoulmob和rvwd這個2個力單獨(dú)提出來計(jì)算呢？？謝謝 2.我是看到有不少文獻(xiàn)的method里，用的是rlist < rcoulmob = rw，而Gromacs官方的例子用的卻是rlist=rcoulmob=rvwd， 當(dāng)時我還很疑惑呢��！我這里想問一下，為什么當(dāng)rlist < rcoulmob = rw的時候計(jì)算資源高，精度更高？？二rlist=rcoulmob=rvwd 計(jì)算快，精度較高？？ 感謝��！ |
4樓2015-01-20 09:16:17
awaken2013
銅蟲 (小有名氣)
- 模擬EPI: 2
- 應(yīng)助: 26 (小學(xué)生)
- 金幣: 137.9
- 紅花: 15
- 帖子: 243
- 在線: 149.6小時
- 蟲號: 2854550
- 注冊: 2013-12-07
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
abinitio: 金幣+2, 鼓勵交流 2015-01-20 07:10:51
zhangmao511: 金幣+20, ★★★★★最佳答案, 非常清晰的講解，謝謝！ 2015-01-20 22:16:02
abinitio: 模擬EPI+1, 系統(tǒng)而全面的解答，加EPI獎勵 2015-01-21 15:38:25
感謝參與，應(yīng)助指數(shù) +1
abinitio: 金幣+2, 鼓勵交流 2015-01-20 07:10:51
zhangmao511: 金幣+20, ★★★★★最佳答案, 非常清晰的講解，謝謝！ 2015-01-20 22:16:02
abinitio: 模擬EPI+1, 系統(tǒng)而全面的解答，加EPI獎勵 2015-01-21 15:38:25
|
rlist是在短程截?cái)喟虢飪?nèi) 所需要計(jì)算的非鍵相互作用。（某些力場排除1-4相互作用，或則按照0.5加成計(jì)算） rcoulmob是短程庫侖力截?cái)喟霃�，在此之�?nèi)的計(jì)算靜電力。在此之外的按照PME實(shí)空間和虛空間算 rvwd是短程范德華力的截?cái)喟霃�，在此之�?nèi)的原子對計(jì)算范德華力，在此之外的力為0 或范德華力矯正。 以上為3種不同的力的計(jì)算截?cái)唷?br /> 因?yàn)樯婕暗洁徑斜硭阉鱊S計(jì)算。 rlist=rcoulmob=rvwd 計(jì)算快，精度較高 rlist<rc/rw 計(jì)算資源高，精度更高 第三種，精度有限。因此第一種用的比較多。 |
2樓2015-01-20 00:21:05
awaken2013
銅蟲 (小有名氣)
- 模擬EPI: 2
- 應(yīng)助: 26 (小學(xué)生)
- 金幣: 137.9
- 紅花: 15
- 帖子: 243
- 在線: 149.6小時
- 蟲號: 2854550
- 注冊: 2013-12-07
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
3樓2015-01-20 00:46:06
awaken2013
銅蟲 (小有名氣)
- 模擬EPI: 2
- 應(yīng)助: 26 (小學(xué)生)
- 金幣: 137.9
- 紅花: 15
- 帖子: 243
- 在線: 149.6小時
- 蟲號: 2854550
- 注冊: 2013-12-07
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
【答案】應(yīng)助回帖
|
rlist確切的概念是鄰近列表搜索的范圍（跟cutoff的概念有點(diǎn)像 但又不完全一致 ，在rlist內(nèi)的原子對的力必須被每步計(jì)算 而在ns期間是保持的。為了節(jié)省計(jì)算資源 rlist外的力就不必每步計(jì)算。也就是當(dāng)rlist=rw=rc時。只需要計(jì)算臨近列表內(nèi)的力即可。計(jì)算資源節(jié)省 當(dāng)rw=rc》rlist時，也就是gromacs的雙程截?cái)�，在rw和rlist間的距離是rlistlong，這范圍內(nèi)的力也是需要每步計(jì)算，當(dāng)然計(jì)算資源就高了。（相比3個相等時） 當(dāng)然第三種情況占用的資源和精度可想而知。 ps： 非鍵作用 可不止有三種。包括二面角 離面運(yùn)動 鍵伸縮等多種勢能 [ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ] |
5樓2015-01-20 09:46:30
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +7 | Doleres 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[論文投稿]
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 085600材料與化工 +5 | 安全上岸！ 2026-03-16 | 5/250 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|