<DIV>Dear all,</DIV>
<DIV>&nbsp;</DIV>
<DIV>Following the MD problems I put forward before, I choose another small step and perfrom MD&nbsp;&nbsp;calculation (no thermostat). The problem still met.&nbsp;</DIV>
<DIV>&nbsp;</DIV>
<DIV>&nbsp; 339&nbsp; 0.00120&nbsp;&nbsp;&nbsp; 30.8&nbsp;&nbsp;&nbsp; -294.15250&nbsp;&nbsp;&nbsp; -294.14388&nbsp;&nbsp;&nbsp; -294.14268&nbsp;&nbsp; 0.011&nbsp;&nbsp; 20.55<BR>&nbsp; 340&nbsp; 0.00120&nbsp;&nbsp;&nbsp; 30.4&nbsp;&nbsp;&nbsp; -294.15241&nbsp;&nbsp;&nbsp; -294.14388&nbsp;&nbsp;&nbsp; -294.14268&nbsp;&nbsp; 0.011&nbsp;&nbsp; 22.52<BR>&nbsp; 341&nbsp; 0.00120&nbsp;&nbsp;&nbsp; 30.1&nbsp;&nbsp;&nbsp; -294.15232&nbsp;&nbsp;&nbsp; -294.14388&nbsp;&nbsp;&nbsp; -294.14268&nbsp;&nbsp; 0.011&nbsp;&nbsp; 20.17<BR>&nbsp; 342&nbsp; 0.00121&nbsp;&nbsp;&nbsp; 29.7&nbsp;&nbsp;&nbsp; -294.15222&nbsp;&nbsp;&nbsp; -294.14389&nbsp;&nbsp;&nbsp; -294.14268&nbsp;&nbsp; 0.011&nbsp;&nbsp; 21.72<BR>&nbsp; 343&nbsp; 0.00121&nbsp;&nbsp;&nbsp; 29.4&nbsp;&nbsp;&nbsp; -294.15212&nbsp;&nbsp;&nbsp; -294.14389&nbsp;&nbsp;&nbsp; -294.14268&nbsp;&nbsp; 0.011&nbsp;&nbsp; 21.73<BR>&nbsp; 344&nbsp; 0.00121&nbsp;&nbsp;&nbsp; 29.0&nbsp;&nbsp;&nbsp; -294.15202&nbsp;&nbsp;&nbsp; -294.14389&nbsp;&nbsp;&nbsp;
 -294.14268&nbsp;&nbsp; 0.010&nbsp;&nbsp; 21.10<BR>&nbsp; 345&nbsp; 0.00122&nbsp;&nbsp;&nbsp; 28.6&nbsp;&nbsp;&nbsp; -294.15192&nbsp;&nbsp;&nbsp; -294.14390&nbsp;&nbsp;&nbsp; -294.14268&nbsp;&nbsp; 0.010&nbsp;&nbsp; 21.93<BR>&nbsp; 346&nbsp; 0.00122&nbsp;&nbsp;&nbsp; 28.2&nbsp;&nbsp;&nbsp; -294.15181&nbsp;&nbsp;&nbsp; -294.14390&nbsp;&nbsp;&nbsp; -294.14268&nbsp;&nbsp; 0.010&nbsp;&nbsp; 19.99<BR>&nbsp; 347&nbsp; 0.00123&nbsp;&nbsp;&nbsp; 27.8&nbsp;&nbsp;&nbsp; -294.15171&nbsp;&nbsp;&nbsp; -294.14390&nbsp;&nbsp;&nbsp; -294.14268&nbsp;&nbsp; 0.010&nbsp;&nbsp; 19.55<BR>&nbsp; 348&nbsp; 0.00123&nbsp;&nbsp;&nbsp; 27.4&nbsp;&nbsp;&nbsp; -294.15160&nbsp;&nbsp;&nbsp; -294.14391&nbsp;&nbsp;&nbsp; -294.14268&nbsp;&nbsp; 0.010&nbsp;&nbsp; 20.57<BR>&nbsp; 349&nbsp; 0.00124&nbsp;&nbsp;&nbsp; 27.0&nbsp;&nbsp;&nbsp; -294.15149&nbsp;&nbsp;&nbsp; -294.14391&nbsp;&nbsp;&nbsp; -294.14267&nbsp;&nbsp; 0.010&nbsp;&nbsp; 20.38<BR>&nbsp; 350&nbsp; 0.00125&nbsp;&nbsp;&nbsp; 26.6&nbsp;&nbsp;&nbsp;
 -294.15138&nbsp;&nbsp;&nbsp; -294.14392&nbsp;&nbsp;&nbsp; -294.14267&nbsp;&nbsp; 0.010&nbsp;&nbsp; 21.35</DIV>
<DIV>&nbsp;RESTART INFORMATION WRITTEN ON FILE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ./RESTART.1<BR>&nbsp; 351&nbsp; 0.00125&nbsp;&nbsp;&nbsp; 26.2&nbsp;&nbsp;&nbsp; -294.15126&nbsp;&nbsp;&nbsp; -294.14392&nbsp;&nbsp;&nbsp; -294.14267&nbsp;&nbsp; 0.010&nbsp;&nbsp; 20.84<BR>&nbsp; 352&nbsp; 0.00126&nbsp;&nbsp;&nbsp; 25.8&nbsp;&nbsp;&nbsp; -294.15115&nbsp;&nbsp;&nbsp; -294.14392&nbsp;&nbsp;&nbsp; -294.14266&nbsp;&nbsp; 0.010&nbsp;&nbsp; 20.53<BR>&nbsp; 353&nbsp; 0.00195&nbsp;&nbsp;&nbsp; 25.3&nbsp;&nbsp;&nbsp; -294.15103&nbsp;&nbsp;&nbsp; -294.14393&nbsp;&nbsp;&nbsp; -294.14198&nbsp;&nbsp; 0.010&nbsp;&nbsp; 20.71<BR>&nbsp; 354138.10491&nbsp;&nbsp;&nbsp; 49.1&nbsp;&nbsp;&nbsp; -294.14967&nbsp;&nbsp;&nbsp; -294.13592&nbsp;&nbsp;&nbsp; -156.03100&nbsp;&nbsp; 0.010&nbsp;&nbsp; 23.11<BR></DIV>
<DIV>The fictitious electronic kinetic energy become very stange at the 354th step. The timestep is very small!</DIV>
<DIV>&nbsp;</DIV>
<DIV>Can anybody tell me what is wrong with it?</DIV>
<DIV>&nbsp;</DIV>
<DIV>Many thanks!</DIV>
<DIV>&nbsp;</DIV>
<DIV>Zhou</DIV><p><br><hr size=1><b>Do You Yahoo!?</b><br>
<a href="http://cn.rd.yahoo.com/mail_cn/tag/?http://cn.tech.yahoo.com/zhuanti/laji/index.html" target=blank>国内电邮用户反垃圾调查拉开帷幕</a>