| 3 | 1/1 | 返回列表 |
| 查看: 548 | 回復(fù): 2 | |||||
ddwy木蟲 (小有名氣)
|
[求助]
分子動力學(xué)模擬出錯 已有1人參與
|
|
大家好: 我用DL_POLY做分子動力學(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ì)算到第步的時候出現(xiàn)錯誤,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 請各位指點(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? 試試從以上提出的可能性來查找錯誤。另外,怎么沒看見shake tolerance關(guān)鍵字? |
木蟲 (小有名氣)

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +10 | 超級伊昂大王 2026-03-24 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +6 | j的立方 2026-03-29 | 6/300 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +5 | 篤行致遠(yuǎn). 2026-03-28 | 5/250 |
|
|
[考研] 085600,專業(yè)課化工原理,320分求調(diào)劑 +4 | 大饞小子 2026-03-29 | 4/200 |
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|