<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=gb2312">
<META content="MSHTML 6.00.2800.1589" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=ΣΧΤ² color=#008000>Polymer exploding</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>Hi:</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>I built Poly(vinyl methyl ether) (PVME)with the 
alkane and ether paramters in OPLS force field. Its repeats number is 22. Tip5p 
was set as solvent.</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>The initial density of PVME melt is 0.8, 1.0 
and 1.05 g/cm3. Its experimental value is 1.05 g/cm3.The initial density of PVME 
solution is 1.0 g/cm3. The concentration is 20% and 60%.</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>The plan runing process is l-bfgs, NVT and 
NPT.</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>However in NVT running, the 3 melt and 60% 
system all exploding. With the help from google, I increase the tau_t from 0.1 
to 0.5. A serial of begin structures were prepared. But it is no helpful. While 
the 20% model is ok. </FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>What is the possible error? How can it fix it? 
</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>Bellow is the md.log file around the exploding 
time point.</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>&nbsp;&nbsp; 
Step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Lambda<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
14140&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
28.28000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>&nbsp;&nbsp; Energies 
(kJ/mol)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Bond&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Angle 
Ryckaert-Bell.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
LJ-14&nbsp;&nbsp;&nbsp;&nbsp; Coulomb-14<BR>&nbsp;&nbsp;&nbsp; 
1.23776e+04&nbsp;&nbsp;&nbsp; 1.97954e+04&nbsp;&nbsp;&nbsp; 
1.57445e+04&nbsp;&nbsp;&nbsp; 6.49124e+03&nbsp;&nbsp; 
-5.38649e+04<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ (SR)&nbsp; Disper. 
corr.&nbsp;&nbsp; Coulomb (SR)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Potential&nbsp;&nbsp;&nbsp; Kinetic En.<BR>&nbsp;&nbsp; -1.32445e+04&nbsp;&nbsp; 
-1.20529e+03&nbsp;&nbsp;&nbsp; 5.64457e+04&nbsp;&nbsp;&nbsp; 
4.25397e+04&nbsp;&nbsp;&nbsp; 3.30808e+04<BR>&nbsp;&nbsp; Total 
Energy&nbsp;&nbsp;&nbsp; Temperature Pressure (bar)<BR>&nbsp;&nbsp;&nbsp; 
7.56204e+04&nbsp;&nbsp;&nbsp; 2.98767e+02&nbsp;&nbsp; -3.82251e+02</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² 
color=#008000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Lambda<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
14150&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
28.30000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>&nbsp;&nbsp; Energies 
(kJ/mol)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Bond&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Angle 
Ryckaert-Bell.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
LJ-14&nbsp;&nbsp;&nbsp;&nbsp; Coulomb-14<BR>&nbsp;&nbsp;&nbsp; 
1.24969e+04&nbsp;&nbsp;&nbsp; 1.97987e+04&nbsp;&nbsp;&nbsp; 
1.57022e+04&nbsp;&nbsp;&nbsp; 6.37298e+03&nbsp;&nbsp; 
-5.38467e+04<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ (SR)&nbsp; Disper. 
corr.&nbsp;&nbsp; Coulomb (SR)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Potential&nbsp;&nbsp;&nbsp; Kinetic En.<BR>&nbsp;&nbsp; -1.22941e+04&nbsp;&nbsp; 
-1.20529e+03&nbsp;&nbsp;&nbsp; 5.64734e+04&nbsp;&nbsp;&nbsp; 
4.34981e+04&nbsp;&nbsp;&nbsp; 4.72135e+04<BR>&nbsp;&nbsp; Total 
Energy&nbsp;&nbsp;&nbsp; Temperature Pressure (bar)<BR>&nbsp;&nbsp;&nbsp; 
9.07116e+04&nbsp;&nbsp;&nbsp; 4.26405e+02&nbsp;&nbsp;&nbsp; 
7.73594e+02</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² 
color=#008000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Lambda<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
14160&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
28.32000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>&nbsp;&nbsp; Energies 
(kJ/mol)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Bond&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Angle 
Ryckaert-Bell.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
LJ-14&nbsp;&nbsp;&nbsp;&nbsp; Coulomb-14<BR>&nbsp;&nbsp;&nbsp; 
4.17034e+16&nbsp;&nbsp;&nbsp; 2.76374e+04&nbsp;&nbsp;&nbsp; 
1.57162e+04&nbsp;&nbsp;&nbsp; 6.48244e+03&nbsp;&nbsp; 
-5.37968e+04<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ (SR)&nbsp; Disper. 
corr.&nbsp;&nbsp; Coulomb (SR)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Potential&nbsp;&nbsp;&nbsp; Kinetic En.<BR>&nbsp;&nbsp; -1.33165e+04&nbsp;&nbsp; 
-1.20529e+03&nbsp;&nbsp;&nbsp; 5.64059e+04&nbsp;&nbsp;&nbsp; 
4.17034e+16&nbsp;&nbsp;&nbsp; 3.87149e+16<BR>&nbsp;&nbsp; Total 
Energy&nbsp;&nbsp;&nbsp; Temperature Pressure (bar)<BR>&nbsp;&nbsp;&nbsp; 
8.04183e+16&nbsp;&nbsp;&nbsp; 3.49651e+14&nbsp;&nbsp; -3.11470e+14</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² 
color=#008000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Lambda<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
14170&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
28.34000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>&nbsp;&nbsp; Energies 
(kJ/mol)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Bond&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Angle 
Ryckaert-Bell.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
LJ-14&nbsp;&nbsp;&nbsp;&nbsp; Coulomb-14<BR>&nbsp;&nbsp;&nbsp; 
4.17033e+16&nbsp;&nbsp;&nbsp; 2.76378e+04&nbsp;&nbsp;&nbsp; 
1.57140e+04&nbsp;&nbsp;&nbsp; 6.48244e+03&nbsp;&nbsp; 
-5.37968e+04<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ (SR)&nbsp; Disper. 
corr.&nbsp;&nbsp; Coulomb (SR)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Potential&nbsp;&nbsp;&nbsp; Kinetic En.<BR>&nbsp;&nbsp; -1.33165e+04&nbsp;&nbsp; 
-1.20529e+03&nbsp;&nbsp;&nbsp; 5.64060e+04&nbsp;&nbsp;&nbsp; 
4.17033e+16&nbsp;&nbsp;&nbsp; 3.87149e+16<BR>&nbsp;&nbsp; Total 
Energy&nbsp;&nbsp;&nbsp; Temperature Pressure (bar)<BR>&nbsp;&nbsp;&nbsp; 
8.04182e+16&nbsp;&nbsp;&nbsp; 3.49651e+14&nbsp;&nbsp; -3.11456e+14</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² 
color=#008000>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Lambda<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
14180&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
28.36000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000>&nbsp;&nbsp; Energies 
(kJ/mol)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Bond&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Angle 
Ryckaert-Bell.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
LJ-14&nbsp;&nbsp;&nbsp;&nbsp; Coulomb-14<BR>&nbsp;&nbsp;&nbsp; 
4.17033e+16&nbsp;&nbsp;&nbsp; 2.76378e+04&nbsp;&nbsp;&nbsp; 
1.57140e+04&nbsp;&nbsp;&nbsp; 6.48244e+03&nbsp;&nbsp; 
-5.37968e+04<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ (SR)&nbsp; Disper. 
corr.&nbsp;&nbsp; Coulomb (SR)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
Potential&nbsp;&nbsp;&nbsp; Kinetic En.<BR>&nbsp;&nbsp; -1.33165e+04&nbsp;&nbsp; 
-1.20529e+03&nbsp;&nbsp;&nbsp; 5.64060e+04&nbsp;&nbsp;&nbsp; 
4.17033e+16&nbsp;&nbsp;&nbsp; 3.87149e+16<BR>&nbsp;&nbsp; Total 
Energy&nbsp;&nbsp;&nbsp; Temperature Pressure (bar)<BR>&nbsp;&nbsp;&nbsp; 
8.04182e+16&nbsp;&nbsp;&nbsp; 3.49651e+14&nbsp;&nbsp; -3.11456e+14</FONT></DIV>
<DIV>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² color=#008000></FONT>&nbsp;</DIV>
<DIV><FONT face=ΣΧΤ² 
color=#008000>*************************************************<BR>Ji 
Qing<BR>Institute of Chemistry, Chinese Academy of Sciences<BR>Tel: 
0086-10-62562894&nbsp; 
£¬82618423<BR>*************************************************</FONT></DIV></BODY></HTML>