 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 5251 | 回復(fù): 11 | |||
longwei1221鐵蟲 (小有名氣)
|
[求助]
forcite做動(dòng)力學(xué)計(jì)算出錯(cuò)了，哪位幫忙看一下是什么原因
|
||
|
我想計(jì)算水的自擴(kuò)散系數(shù)，用幾種方法都出錯(cuò)了。下面是其中一種使用forcite來計(jì)算的結(jié)果，哪位看看什么參數(shù)設(shè)置錯(cuò)誤導(dǎo)致中斷的？ 我分兩種方法建模，一個(gè)是用AC模塊設(shè)置密度為1，另外一種是直接建模用forcite的幾何優(yōu)化來建模再做動(dòng)力學(xué)模擬。同時(shí)在計(jì)算動(dòng)力學(xué)我也嘗試了discovery和forcite，結(jié)果都中斷失敗了，不知道什么原因，哪位幫忙看一下吧，謝謝啦！ 下面分別是forcite和discovery的錯(cuò)誤提示。另外如果誰知道計(jì)算擴(kuò)散系數(shù)方法，麻煩告知一聲，不勝感激，我愿意奉送全部的金幣 Forcite ------- Task : Dynamics Version : 5.0 Build date : Oct 22 2009 Host : meilongwei-PC Threads : Parallel 3 Operating system : Windows Task started : Tue Oct 23 16:48:46 2012 ---- Dynamics parameters ---- Ensemble : NPT Temperature : 293.00 K Control method : Velocity Scale Temp. difference : 10.00 K Pressure : 0.00010 GPa Control method : Andersen Cell time constant : 1.0000000 ps Timestep : 10.00 fs Number of steps : 1000000 Duration : 10000 ps Initial velocities : Random Random number seed : 1350982127 ---- Energy parameters ---- Forcefield : Universal Electrostatic terms: Summation method : Atom based Truncation method : Cubic spline Cutoff distance : 15.5 A Spline width : 1 A Buffer width : 0.5 A van der Waals terms: Summation method : Atom based Truncation method : Cubic spline Cutoff distance : 15.5 A Spline width : 1 A Buffer width : 0.5 A Dynamics calculation - Aborting dynamics. An unphysical integration step (>100 Angstroms) has been detected. Try reducing the time step or relaxing the structure first. discovery錯(cuò)誤提示 DISCOVER Molecular Simulation Program Version: 2009.1 Build: Oct 22 2009 Date: Tue Oct 23 15:49:33 2012 User Name: MEILONGWEI-PC$ Host Name: meilongwei-PC Host Type: Windows Threads: 1 --------------------------------------------------------------- Checked out license feature: MS_discover_MP --------------------------------------------------------------- randomSeed is set to 978573 Line 7:BTCL> autoEcho off --------------------------------------------------------------- Checked out license feature: MS_compass_MP --------------------------------------------------------------- INPUT FILES ___________ File Type Name --------- ---- Forcefield C:\PROGRA~1\Accelrys\MATERI~1.0\etc\Gateway\..\..\share\Discover\res/compass.frc Molecular data Molecule1.mdf Coordinate Molecule1.car Periodic Boundary Conditions ____________________________ Length (A) Angle (degrees) ---------- --------------- a 10.61650 alpha 90.00000 b 10.61650 beta 90.00000 c 10.61650 gamma 90.00000 MOLECULAR TOPOLOGY __________________ Number of molecules: 40 Number of residues: 40 Number of atoms: 120 (asymmetric unit: 120) Number of atom types: 3 Number of bonds: 80 Number of consolidated angles: 40 FORCEFIELD OPTIONS __________________ Filename : compass.frc Definition name : cff91 Version : 2.8 Last modification date : 3/6/2007 # of automatic parameters : 0 NONBOND ENERGY CUTOFFS ______________________ Cutoff (A) Spline Width (A) Buffer Width (A) vdW 9.50 1.00 0.50 Coulomb -- -- -- NONBOND ENERGY EWALD PARAMETERS _______________________________ Accuracy (kcal/mol) Update Width (A) vdW -- -- Coulomb 0.01000000 1.00 Summation method for vdW : Atom based Summation method for Coulomb : Ewald Dielectric : 1.00 NOTE: Pressure and energy tail corrections will be added to all subsequent energy evaluations MOLECULAR DYNAMICS __________________ Ensemble : NPT Temperature : 293.00 K Control Method : Andersen, Collision frequency is every 68 steps Pressure : 0.0001 GPa Stress : Sxx = 0.0000 GPa Syy = 0.0000 GPa Szz = 0.0000 GPa : Syz = 0.0000 GPa Sxz = 0.0000 GPa Sxy = 0.0000 GPa Control Method : Berendsen , tau = 0.1443 psec , compressibility = 0.5000 Timestep : 1.00 fs Duration :1000000.00 fs Integration Method : Velocity Verlet Initial Velocities : Random Velocities from Boltzmann distribution Initial Temp. : 293.00 K Error: Dynamics: Energy difference between successive steps in dynamics is greater than the user defined variable DEVIATION, 5000.000 invoked from within "dynamics \ time = 1000000.00 \ timestep = 1.00 \ initial_temperature = 293.00 \ ensemble = NPT \ temperature_control_method = ande ..." Total time used by DISCOVER: 24 mins 7.05 secs (1447.05 secs) Completion date: Tue Oct 23 16:13:41 2012 Exiting Discover: An ERROR has occurred! |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 分子模擬 |
» 猜你喜歡
288求調(diào)劑 一志愿哈工大 材料與化工
已經(jīng)有31人回復(fù)
-
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑
已經(jīng)有5人回復(fù)
-
085801 總分275 本科新能源 求調(diào)劑
已經(jīng)有17人回復(fù)
-
324求調(diào)劑
已經(jīng)有6人回復(fù)
-
學(xué)碩化學(xué)工程與技術(shù)，一志愿中國海洋大學(xué)320+求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料調(diào)劑
已經(jīng)有5人回復(fù)
-
初試301，代碼085701環(huán)境工程，本碩一致，四六級(jí)已過，有二區(qū)一作，共發(fā)表5篇論文
已經(jīng)有6人回復(fù)
-
化學(xué)308分調(diào)劑
已經(jīng)有15人回復(fù)
-
346求調(diào)劑
已經(jīng)有3人回復(fù)
-
環(huán)境285分，過六級(jí)，求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 用MS中的Forcite計(jì)算冰，結(jié)果發(fā)現(xiàn)溫度波動(dòng)異常的大，感覺像做周期性的震蕩衰減
已經(jīng)有6人回復(fù)
- 【討論】discover或者forcite計(jì)算擴(kuò)散系數(shù)
已經(jīng)有5人回復(fù)
- 【求助】MS forcite動(dòng)力學(xué)問題
已經(jīng)有13人回復(fù)
- 【求助成功】MS模擬MCM-41過程中forcite filed計(jì)算的問題
已經(jīng)有12人回復(fù)
longwei1221
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 116.8
- 散金: 128
- 紅花: 20
- 帖子: 191
- 在線: 81.7小時(shí)
- 蟲號(hào): 1439584
- 注冊: 2011-10-12
- 性別: GG
- 專業(yè): 核技術(shù)及其應(yīng)用
9樓2012-11-03 11:18:11
wangzhikun
金蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 783.5
- 散金: 62
- 紅花: 2
- 帖子: 151
- 在線: 110.9小時(shí)
- 蟲號(hào): 1253117
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
2樓2012-10-23 19:54:35
longwei1221
鐵蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 116.8
- 散金: 128
- 紅花: 20
- 帖子: 191
- 在線: 81.7小時(shí)
- 蟲號(hào): 1439584
- 注冊: 2011-10-12
- 性別: GG
- 專業(yè): 核技術(shù)及其應(yīng)用
|
麻煩幫忙看一下這個(gè)錯(cuò)誤是什么原因?qū)е碌暮脝�？被告知能量太大，不知道該如何設(shè)置可以減小這個(gè)能量 Forcite ------- Task : Dynamics Version : 5.0 Build date : Oct 22 2009 Host : meilongwei-PC Threads : Parallel 3 Operating system : Windows Task started : Wed Oct 24 10:17:10 2012 ---- Dynamics parameters ---- Ensemble : NPT Temperature : 300.00 K Control method : Andersen Collision ratio : 1.0000000 Pressure : 0.00010 GPa Control method : Andersen Cell time constant : 1.0000000 ps Timestep : 1.00 fs Number of steps : 1000000 Duration : 1000 ps Initial velocities : Random Random number seed : 1351045030 ---- Energy parameters ---- Forcefield : Universal Electrostatic terms: Summation method : Atom based Truncation method : Cubic spline Cutoff distance : 18.5 A Spline width : 1 A Buffer width : 0.5 A van der Waals terms: Summation method : Atom based Truncation method : Cubic spline Cutoff distance : 18.5 A Spline width : 1 A Buffer width : 0.5 A ---- Thermostat ---- Andersen collision period: every 222 steps. Dynamics calculation - Deviation in total energy (2.21112e+012) between successive steps was greater than the specified maximum (5000). Task terminated : Wed Oct 24 10:17:17 2012 Total CPU time used by Forcite: 4 seconds (3.82s) Termination status : Error |
3樓2012-10-24 10:20:38
huangshuang1
木蟲 (著名寫手)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 2762.2
- 帖子: 2263
- 在線: 464小時(shí)
- 蟲號(hào): 1886296
- 注冊: 2012-07-10
- 性別: MM
- 專業(yè): 機(jī)械動(dòng)力學(xué)
4樓2012-10-24 10:29:46
wangzhikun
金蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 783.5
- 散金: 62
- 紅花: 2
- 帖子: 151
- 在線: 110.9小時(shí)
- 蟲號(hào): 1253117
- 注冊: 2011-04-02
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
5樓2012-10-24 16:25:00
yangyy11
銀蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 541.1
- 帖子: 102
- 在線: 54.2小時(shí)
- 蟲號(hào): 1942486
- 注冊: 2012-08-17
- 性別: GG
- 專業(yè): 光學(xué)和光電子材料
【答案】應(yīng)助回帖
★ ★
感謝參與，應(yīng)助指數(shù) +1
jiaoyixiong: 金幣+1, 鼓勵(lì)交流 2012-10-26 14:54:37
longwei1221: 金幣+1, 謝謝 2012-11-03 11:18:29
感謝參與，應(yīng)助指數(shù) +1
jiaoyixiong: 金幣+1, 鼓勵(lì)交流 2012-10-26 14:54:37
longwei1221: 金幣+1, 謝謝 2012-11-03 11:18:29
| 這里有個(gè)水分子的擴(kuò)散系數(shù)的模擬教程，應(yīng)該能解決樓主的問題！http://www.gaoyang168.com/bbs/viewthread.php?tid=5096305 |
6樓2012-10-25 17:02:29
sikasun
木蟲 (正式寫手)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 3822.2
- 散金: 696
- 紅花: 1
- 帖子: 800
- 在線: 181.5小時(shí)
- 蟲號(hào): 632362
- 注冊: 2008-10-21
- 性別: MM
- 專業(yè): 環(huán)境污染化學(xué)
加那么多不影響結(jié)果嗎 7樓2012-10-26 13:52:02
憤怒的姿態(tài)
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1168.9
- 散金: 22
- 帖子: 158
- 在線: 78.4小時(shí)
- 蟲號(hào): 1652991
- 注冊: 2012-02-29
- 性別: GG
- 專業(yè): 實(shí)驗(yàn)地球化學(xué)和計(jì)算地球化
【答案】應(yīng)助回帖
★ ★ ★
zh1987hs: 金幣+2, 謝謝 2012-10-30 22:37:58
longwei1221: 金幣+1, 謝謝，我知道錯(cuò)誤在這里，可是不知道如何修改 2012-11-03 11:18:55
zh1987hs: 金幣+2, 謝謝 2012-10-30 22:37:58
longwei1221: 金幣+1, 謝謝，我知道錯(cuò)誤在這里，可是不知道如何修改 2012-11-03 11:18:55
|
同志，你的錯(cuò)誤很明顯啊，那不是明顯的error嗎？Error: Dynamics: Energy difference between successive steps in dynamics is greater than the user defined variable DEVIATION, 5000.000 Dynamics calculation - Deviation in total energy (2.21112e+012) between successive steps was greater than the specified maximum (5000). 樓主自己看下，在對應(yīng)的選項(xiàng)改下就好了！ 希望對你有用！祝好！ |
8樓2012-10-30 21:23:27
憤怒的姿態(tài)
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1168.9
- 散金: 22
- 帖子: 158
- 在線: 78.4小時(shí)
- 蟲號(hào): 1652991
- 注冊: 2012-02-29
- 性別: GG
- 專業(yè): 實(shí)驗(yàn)地球化學(xué)和計(jì)算地球化
10樓2012-11-05 08:50:57
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 longwei1221 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程322分 +8 | 哈哈哈吼吼吼哈 2026-04-01 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 化工調(diào)劑303分，過四級(jí) +5 | 棲梧待風(fēng) 2026-04-02 | 5/250 |
|
|
|
[考研] 337求調(diào)劑 +11 | 《樹》 2026-03-29 | 11/550 |
|
|
|
[考研] 279求調(diào)劑 +6 | 學(xué)而思兮知 2026-04-01 | 6/300 |
|
|
|
[考研] 348環(huán)境工程調(diào)劑 +3 | 吳彥祖24k 2026-04-01 | 3/150 |
|
|
|
[考研] 一志愿北交大材料工程，總分358 +4 | cs0106 2026-04-01 | 4/200 |
|
|
|
[考研] 379求調(diào)劑 +3 | ?苦瓜不苦 2026-04-01 | 3/150 |
|
|
|
[考研] 一志愿西安交大材料學(xué)碩（英一數(shù)二）347，求調(diào)劑到高分子/材料相關(guān)專業(yè) +7 | zju51 2026-03-31 | 9/450 |
|
|
|
[考研] 283求調(diào)劑 +9 | A child 2026-03-28 | 9/450 |
|
|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +9 | 吃的不少 2026-03-30 | 9/450 |
|
|
|
[考研] 生物學(xué)296求調(diào)劑 +10 | 湯圓包 2026-03-29 | 14/700 |
|
|
|
[考研] 0710生物學(xué)求調(diào)劑！ +6 | 敘述文 2026-03-31 | 6/300 |
|
|
|
[考研] 一志愿食品科學(xué)與工程083200求調(diào)劑 +4 | XQTJZ 2026-03-30 | 4/200 |
|
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|