 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 2767 | 回復(fù): 2 | ||
[求助]
如何用amber自帶程序處理軌跡中由于周期性邊界條件，溶劑水分子錯(cuò)亂鏡像成鍵的問題？ 已有1人參與
|
|
有2個(gè)問題求救大家： 第一：如圖，我用amber跑的軌跡，有周期性邊界條件，最后呢，水盒子的水分子就鏡像成鍵了，出現(xiàn)了很多跨越成鍵現(xiàn)象。請(qǐng)問，amber里自帶的程序如何轉(zhuǎn)換或者怎么處理下軌跡，總之我是想把這些水分子的長(zhǎng)程化學(xué)鍵去掉，讓氧原子和周圍正常的2個(gè)氫原子去成鍵。請(qǐng)問，如何處理，在manual里哪塊有說明嗎？或者官網(wǎng)的哪個(gè)練習(xí)有類似操作嗎？比較棘手，多謝各位了！ 我的mdin文件如下： dhfr evb &cntrl ievb = 1, irest = 0, ntx=1, ntb = 1, ntt = 3, gamma_ln=100.0, cut = 10.0, nstlim = 200000, dt = 0.001, ntpr = 200, ntave=1, ntwr = 200, ntwx = 200, ntwe = 200, nscm = 100, temp0 = 300, ntf = 1, ntc = 1, vlimit=20.0, iwrap=1, / &ewald nfft1 = 120, nfft2 = 120, nfft3 = 120, / ~ ~ ~ 第2：我的軌跡mdcrd是200ps，1000個(gè)frames，第一個(gè)frame是正常的，最后一個(gè)就如圖了。然后我用這兩個(gè)frame分別取算單點(diǎn)能，最后一個(gè)frame的能量是第一個(gè)的1000倍，很明顯不對(duì)，它們應(yīng)該差不多的。而且bond energy不對(duì)，我的單點(diǎn)能mdin文件如下： dhfrc1 nvt &cntrl ievb=1, imin=0, irest = 0, ntx=1, nstlim=1, dt = 0.0, ntb=1, ntf = 1, ntc = 1, ntt = 3, gamma_ln=100.0, temp0 = 300, cut=12.0, ntpr = 1, ntave=1, ntwr = 1, ntwx = 1, ntwe = 1, vlimit=20.0, iwrap=1, / &ewald nfft1 = 120, nfft2 = 120, nfft3 = 120, / 輸出的能量文件如下： NSTEP = 0 TIME(PS) = 0.000 TEMP(K) = 0.00 PRESS = 0.0 Etot = ************** EKtot = 0.0000 EPtot = ************** BOND = ************** ANGLE = 1438.2584 DIHED = 1764.9645 1-4 NB = 606.6218 1-4 EEL = 5867.2479 VDWAALS = 33125.4300 EELEC = -268477.1384 EHBOND = 0.0000 RESTRAINT = 0.0000 EVB:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ V_TOT = ************** V_EVB = ************** V_UMB = 0.0000 C_0^2 = 0.9915 0.0085 EVB RC = 0.0000 Ewald error estimate: 0.7715E-01 ------------------------------------------------------------------------------ NSTEP = 1 TIME(PS) = 0.000 TEMP(K) =305151.29 PRESS = 0.0 Etot = ************** EKtot = 51070814.7000 EPtot = ************** BOND = ************** ANGLE = 1438.2584 DIHED = 1764.9645 1-4 NB = 606.6218 1-4 EEL = 5867.2479 VDWAALS = 33125.4300 EELEC = -268477.1384 EHBOND = 0.0000 RESTRAINT = 0.0000 EVB:~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ V_TOT = ************** V_EVB = ************** V_UMB = 0.0000 C_0^2 = 0.9915 0.0085 EVB RC = 0.0000 Ewald error estimate: 0.7715E-01 ------------------------------------------------------------------------------ 我的bond energy老是算不對(duì)，請(qǐng)問這是怎么回事呢？ manual上哪里有提到嗎？ 大家如果了解的請(qǐng)多多幫忙，新學(xué)amber不久 ，不怕遇到問題，只是對(duì)于遇到的問題很沒有頭緒，不知道該怎么解決，多謝了！ 3`O13KF4GZC(@F6F~EWNR%W.jpg |
回復(fù)此樓» 猜你喜歡
274求調(diào)劑
已經(jīng)有7人回復(fù)
-
北科281學(xué)碩材料求調(diào)劑
已經(jīng)有11人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有3人回復(fù)
-
302求調(diào)劑
已經(jīng)有3人回復(fù)
-
312求調(diào)劑
已經(jīng)有5人回復(fù)
-
網(wǎng)絡(luò)空間安全0839招調(diào)劑
已經(jīng)有4人回復(fù)
-
329求調(diào)劑
已經(jīng)有3人回復(fù)
-
【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額
已經(jīng)有8人回復(fù)
-
085602 289分求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- Amber 分子動(dòng)力學(xué) 求助�。。�！
已經(jīng)有10人回復(fù)
匿名
用戶注銷 (小有名氣)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 320.7
- 散金: 1
- 帖子: 86
- 在線: 141.1小時(shí)
- 蟲號(hào): 0
- 注冊(cè): 2011-11-08
- 專業(yè): 理論和計(jì)算化學(xué)
★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
lifeisall(月只藍(lán)代發(fā)): 金幣+3, 鼓勵(lì)交流！ 2014-07-08 16:29:34
感謝參與，應(yīng)助指數(shù) +1
lifeisall(月只藍(lán)代發(fā)): 金幣+3, 鼓勵(lì)交流！ 2014-07-08 16:29:34
 本帖僅樓主可見 |
2樓2014-07-03 23:24:32
3樓2014-07-06 18:16:12
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 lifeisall 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 085602 289分求調(diào)劑
+7
|
WWW西西弗斯 2026-03-24 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研]
材料與化工考研調(diào)劑
+9
|
孅華 2026-03-22 | 9/450 |
|
|
[考研] 359求調(diào)劑 +3 | 王了個(gè)楠 2026-03-25 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +7 | shxchem 2026-03-20 | 9/450 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +12 | 用戶573181 2026-03-20 | 17/850 |
|
|
|
[考研] 321求調(diào)劑 +4 | Ymlll 2026-03-24 | 4/200 |
|
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
|
[考研] 291 求調(diào)劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 5/250 |
|
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 7/350 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
|
[考研]
 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
|
[考研]
材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|