| 3 | 1/1 | 返回列表 |
| 查看: 547 | 回復(fù): 2 | ||||
ddwy木蟲 (小有名氣)
|
[求助]
分子動(dòng)力學(xué)模擬出錯(cuò) 已有1人參與
|
|
大家好: 我用DL_POLY做分子動(dòng)力學(xué)模擬,CONTROL輸入文件是: temperature 298.15 # pressure 0.001 ensemble nvt hoover 0.5 integrat velocity restart steps 3000000 equilibration 1000000 mult 1 scale 5 stats 20 rdf 10 print rdf traj 1000001 100 0 densvar 20 timestep 0.002 cutoff 12.0 ewald precision 1.0E-6 minimise energy 50 5.0E-03 shake 1.0E-6 mxshak 10000 quaternion tolerance 1.0E-05 job time 99999999999999.0 close time 200000.0 finish 但是計(jì)算到第步的時(shí)候出現(xiàn)錯(cuò)誤,OUTPUT文件為: ---------------------------------------------------------------------------------------------------------------------------------- 174600 -2.8822E+04 2.9815E+02 -3.1606E+04 3.9295E+03 -3.5876E+04 1.7169E+02 1.0398E+02 6.5269E+01 0.0000E+00 349.20000 -2.9405E+04 0.0000E+00 9.0503E+03 -9.0711E+04 3.5875E+04 8.3469E+01 0.0000E+00 6.3803E+04 0.0000E+00 119751.3 4.9579E+04 0.0000E+00 0.0000E+00 0.0000E+00 9.0000E+01 9.0000E+01 9.0000E+01 0.0000E+00 -1.2860E+00 rolling -2.6814E+04 2.9487E+02 -2.9563E+04 4.1468E+03 -3.4341E+04 3.6908E+02 1.8861E+02 7.4187E+01 0.0000E+00 averages -2.6525E+04 0.0000E+00 6.3658E+03 -9.3434E+04 3.4340E+04 1.4041E+03 0.0000E+00 6.4056E+04 0.0000E+00 4.9579E+04 0.0000E+00 0.0000E+00 0.0000E+00 9.0000E+01 9.0000E+01 9.0000E+01 0.0000E+00 -8.0166E-02 ---------------------------------------------------------------------------------------------------------------------------------- *** warning - global constraint number 6 , with particle numbers: 2885 & 2886 , converges to a length of 7.94 Angstroms with a factor 3.1039E+01 ,contributes towards next error !!! *** *** warning - global constraint number 7 , with particle numbers: 2885 & 2887 , converges to a length of 12.88 Angstroms with a factor 8.2412E+01 ,contributes towards next error !!! *** *** warning - global constraint number 8 , with particle numbers: 2886 & 2887 , converges to a length of 18.15 Angstroms with a factor 1.0006E+02 ,contributes towards next error !!! *** DL_POLY_4 terminated due to error 105 error - shake algorithm (constraints_shake) failed to converge "OUTPUT" 64939L, 6075639C 請(qǐng)各位指點(diǎn)一下,謝謝 |

專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
http://www.uark.edu/ua/fengwang/DLPOLY2/node350.html Message 105: error - shake algorithm failed to converge The RD-SHAKE algorithm for bond constraints is iterative. If the maximum number of permitted iterations is exceeded, the program terminates. Possible causes include: a bad starting configuration; too large a time step used; incorrect force field specification; too high a temperature; inconsistent constraints involving shared atoms etc. Action: Corrective action depends on the cause. It is unlikely that simply increasing the iteration number will cure the problem, but you can try: follow the standard user response to increase the control parameter mxshak. But the trouble is much more likely to be cured by careful consideration of the physical system being simulated. For example, is the system stressed in some way? Too far from equilibrium? 試試從以上提出的可能性來查找錯(cuò)誤。另外,怎么沒看見shake tolerance關(guān)鍵字? |
木蟲 (小有名氣)

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +13 | 田洪有 2026-03-26 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
[考研] 309求調(diào)劑 +10 | 誰不是少年 2026-03-29 | 10/500 |
|
|
[考研] 317求調(diào)劑 +10 | 蛋黃咸肉粽 2026-03-26 | 10/500 |
|
|
[考研] 求調(diào)劑 +5 | 青春裁為三截 2026-03-29 | 5/250 |
|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè),初試351分,復(fù)試被刷求調(diào)劑! +5 | 真名有冰 2026-03-29 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研] 329求調(diào)劑 +7 | 星野? 2026-03-26 | 7/350 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|