 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3762 | 回復(fù): 13 | ||||
| 本帖產(chǎn)生 1 個 模擬EPI ，點擊這里進行查看 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||
zhangmao511鐵蟲 (小有名氣)
|
[求助]
Gromacs里關(guān)于截斷半徑的一個問題，謝謝��！ 已有1人參與
|
|||
|
在讀Gromacs手冊的時候，有2段話描述了截斷半徑。對這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的時候才計算這些原子之間的短程相互作用。相互作用的原子對，如果距離在rlist之內(nèi)的就被存儲下來，沒有落到rlist內(nèi)但是落到max(rcoulomb,rvdw)內(nèi)的，會在鄰近搜索期間被計算。 那么rlist和rcoulomb, rvdw具體是什么關(guān)系呢？？不知道能不能解釋的詳細(xì)點，謝謝！ 還有最后一句話，“If rlist = max(rcoulomb,rvdw), no forces are evaluated during neighbor list generation”，不怎么理解，請大家指教，謝謝�。� |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 MDs-Gromacs | gromacs |
» 猜你喜歡
【求調(diào)劑】085601材料工程專碩 | 總分272 |
已經(jīng)有5人回復(fù)
-
275求調(diào)劑
已經(jīng)有13人回復(fù)
-
0856，材料與化工321分求調(diào)劑
已經(jīng)有8人回復(fù)
-
085701環(huán)境工程，267求調(diào)劑
已經(jīng)有12人回復(fù)
-
320分，材料與化工專業(yè)，求調(diào)劑
已經(jīng)有7人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
295材料工程專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
265求調(diào)劑
已經(jīng)有8人回復(fù)
-
322求調(diào)劑
已經(jīng)有4人回復(fù)
-
085602 307分 求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 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ù)
- 截斷半徑選取的標(biāo)準(zhǔn)是什么？
已經(jīng)有6人回復(fù)
awaken2013
銅蟲 (小有名氣)
- 模擬EPI: 2
- 應(yīng)助: 26 (小學(xué)生)
- 金幣: 137.9
- 紅花: 15
- 帖子: 243
- 在線: 149.6小時
- 蟲號: 2854550
- 注冊: 2013-12-07
- 性別: GG
- 專業(yè): 生物大分子結(jié)構(gòu)與功能
|
1 listlong: Cut-off distance for the long-range neighbor list. This parameter is only relevant for a twin- range cut-off setup with switched potentials. 2 Long-range interactions are cut off In this version, GROMACS always uses a cut-off radius for the Lennard-Jones interactions and sometimes for the Coulomb interactions as well 3 Twin-range cut-offs To save computation time, slowly varying forces can be calculated less often than rapidly varying forces. In GROMACS such a multiple time step splitting is possible between short and long range non-bonded interactions. In this integrator the long-range forces are determined every n steps and are then integrated into the velocity in eqn. 3.25 using a time step of tLR = nt Twin range cut-offs with neighbor list cut-off rlist and VdW cut-off rvdw, where rvdw >rlist 4 nstlist Frequency to update the neighbor list (and the long-range forces, when using twin- range cut-offs). |
8樓2015-01-20 13:50:20
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是在短程截斷半斤內(nèi) 所需要計算的非鍵相互作用。（某些力場排除1-4相互作用，或則按照0.5加成計算） rcoulmob是短程庫侖力截斷半徑，在此之內(nèi)的計算靜電力。在此之外的按照PME實空間和虛空間算 rvwd是短程范德華力的截斷半徑，在此之內(nèi)的原子對計算范德華力，在此之外的力為0 或范德華力矯正。 以上為3種不同的力的計算截斷。 因為涉及到鄰近列表搜索NS計算。 rlist=rcoulmob=rvwd 計算快，精度較高 rlist<rc/rw 計算資源高，精度更高 第三種，精度有限。因此第一種用的比較多。 |
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
zhangmao511
鐵蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 33
- 紅花: 2
- 帖子: 136
- 在線: 29.7小時
- 蟲號: 2553312
- 注冊: 2013-07-19
- 性別: GG
- 專業(yè): 蛋白質(zhì)組學(xué)
|
感謝你的講解。  1. 你說rlist是是非共價鍵相互作用的截斷半徑。而非共價鍵就包括了靜電相互作用和范德華力 ，那我能理解rlist包含了rcoulmob和rvwd嗎？那既然是這樣，為什么要把rcoulmob和rvwd這個2個力單獨提出來計算呢？？謝謝 2.我是看到有不少文獻的method里，用的是rlist < rcoulmob = rw，而Gromacs官方的例子用的卻是rlist=rcoulmob=rvwd， 當(dāng)時我還很疑惑呢��！我這里想問一下，為什么當(dāng)rlist < rcoulmob = rw的時候計算資源高，精度更高？？二rlist=rcoulmob=rvwd 計算快，精度較高？？ 感謝！！ |
4樓2015-01-20 09:16:17
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 322求調(diào)劑 +4 | 宋明欣 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +4 | 辰沐5211314 2026-03-26 | 4/200 |
|
|
|
[考研] 085600，材料與化工321分，求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研] 349求調(diào)劑 +5 | 杰斯塔里斯 2026-03-21 | 5/250 |
|
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
|
[考研] 總分293求調(diào)劑 +6 | 加一一九 2026-03-25 | 8/400 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。，。，。，。i 2026-03-23 | 6/300 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|