<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=utf-8">
<META content="MSHTML 6.00.2800.1589" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=幼圆>
<P class=MsoNormal 
style="MARGIN: 0cm 0cm 0pt; TEXT-ALIGN: left; mso-pagination: widow-orphan" 
align=left><SPAN lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 幼圆; mso-hansi-font-family: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">Tsjerk 
Wassenaar, thank you very much for you suggestion even if most of them have been 
manipulated.</SPAN><SPAN lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt; mso-ascii-font-family: 幼圆; mso-fareast-font-family: 幼圆">&nbsp;</SPAN><SPAN 
lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt"><?xml:namespace 
prefix = o ns = "urn:schemas-microsoft-com:office:office" 
/><o:p></o:p></SPAN></P>
<P class=MsoNormal 
style="MARGIN: 0cm 0cm 0pt; TEXT-ALIGN: left; mso-pagination: widow-orphan" 
align=left><SPAN lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">&nbsp;<o:p></o:p></SPAN></P>
<P class=MsoNormal 
style="MARGIN: 0cm 0cm 0pt; TEXT-ALIGN: left; mso-pagination: widow-orphan" 
align=left><I><SPAN lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">"Write 
out the energies at every step and check which of the energy<BR>contributions 
goes astray. This may give you a good clue to where<BR>you're mistake is. Also, 
run through the topology again, by hand,<BR>checking all bonds, angles, 
dihedrals, exclusion (implicit ones!),<BR>etc."</SPAN></I><SPAN lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt"><o:p></o:p></SPAN></P>
<P class=MsoNormal 
style="MARGIN: 0cm 0cm 0pt; TEXT-ALIGN: left; mso-pagination: widow-orphan" 
align=left><SPAN lang=EN-US 
style="FONT-SIZE: 12pt; COLOR: green; FONT-FAMILY: 幼圆; mso-hansi-font-family: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">when 
the explosion happened, the energy of bond, angle and potential will were 
increase. And the topology has been checked several times. </SPAN><SPAN 
lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt"><o:p></o:p></SPAN></P>
<P class=MsoNormal 
style="MARGIN: 0cm 0cm 0pt; TEXT-ALIGN: left; mso-pagination: widow-orphan" 
align=left><SPAN lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">&nbsp;<o:p></o:p></SPAN></P>
<P class=MsoNormal 
style="MARGIN: 0cm 0cm 0pt; TEXT-ALIGN: left; mso-pagination: widow-orphan" 
align=left><SPAN lang=EN-US 
style="FONT-SIZE: 12pt; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">&nbsp;<o:p></o:p></SPAN></P>
<P class=MsoNormal 
style="MARGIN: 0cm 0cm 0pt; TEXT-ALIGN: left; mso-pagination: widow-orphan" 
align=left><SPAN lang=EN-US 
style="FONT-SIZE: 24pt; COLOR: green; FONT-FAMILY: 幼圆; mso-hansi-font-family: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">Now 
I</SPAN><SPAN lang=EN-US 
style="FONT-SIZE: 24pt; COLOR: green; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt; mso-ascii-font-family: 幼圆; mso-fareast-font-family: 幼圆">&nbsp;</SPAN><SPAN 
lang=EN-US 
style="FONT-SIZE: 24pt; COLOR: green; FONT-FAMILY: 幼圆; mso-hansi-font-family: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">increase 
the cutoff of VDW &amp; Coulomb from 0.9 to 1.5. The system</SPAN><SPAN 
lang=EN-US 
style="FONT-SIZE: 24pt; COLOR: green; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt; mso-ascii-font-family: 幼圆; mso-fareast-font-family: 幼圆">&nbsp;</SPAN><SPAN 
lang=EN-US 
style="FONT-SIZE: 24pt; COLOR: green; FONT-FAMILY: 幼圆; mso-hansi-font-family: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">looks 
ok from last night to this time. If the problem is fixed, I conclude that cutoff 
0.9</SPAN><SPAN lang=EN-US 
style="FONT-SIZE: 24pt; COLOR: green; FONT-FAMILY: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt; mso-ascii-font-family: 幼圆; mso-fareast-font-family: 幼圆">&nbsp;</SPAN><SPAN 
lang=EN-US 
style="FONT-SIZE: 24pt; COLOR: green; FONT-FAMILY: 幼圆; mso-hansi-font-family: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt">is 
not enough for electric system.</SPAN></P>
<P class=MsoNormal 
style="MARGIN: 0cm 0cm 0pt; TEXT-ALIGN: left; mso-pagination: widow-orphan" 
align=left><SPAN lang=EN-US 
style="FONT-SIZE: 24pt; COLOR: green; FONT-FAMILY: 幼圆; mso-hansi-font-family: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt"></SPAN>&nbsp;</P>
<P class=MsoNormal 
style="MARGIN: 0cm 0cm 0pt; TEXT-ALIGN: left; mso-pagination: widow-orphan" 
align=left><SPAN lang=EN-US 
style="FONT-SIZE: 24pt; COLOR: green; FONT-FAMILY: 幼圆; mso-hansi-font-family: 宋体; mso-bidi-font-family: 宋体; mso-font-kerning: 0pt"><FONT 
face=宋体 color=#000000 size=3>Hope this conclusion is helpful for GMX 
users.</FONT></SPAN></P></FONT></DIV>
<DIV><FONT face=幼圆 size=6></FONT>&nbsp;</DIV>
<DIV><FONT face=幼圆></FONT>&nbsp;</DIV>
<DIV><FONT face=幼圆></FONT>&nbsp;</DIV>
<DIV><EM>Message: 2<BR>Date: Wed, 7 Mar 2007 17:49:30 +0100<BR>From: "Tsjerk 
Wassenaar" &lt;</EM><A 
href="mailto:tsjerkw@gmail.com"><EM>tsjerkw@gmail.com</EM></A><EM>&gt;<BR>Subject: 
Re: Re: [gmx-users] vacuum right but melt collapse<BR>To: "kitty ji" &lt;</EM><A 
href="mailto:jiqing@iccas.ac.cn"><EM>jiqing@iccas.ac.cn</EM></A><EM>&gt;, 
"Discussion list for GROMACS<BR>users" &lt;</EM><A 
href="mailto:gmx-users@gromacs.org"><EM>gmx-users@gromacs.org</EM></A><EM>&gt;<BR>Message-ID:<BR>&lt;</EM><A 
href="mailto:8ff898150703070849r59cb8a8dr60b894a881fcc331@mail.gmail.com"><EM>8ff898150703070849r59cb8a8dr60b894a881fcc331@mail.gmail.com</EM></A><EM>&gt;<BR>Content-Type: 
text/plain; charset=GB2312; format=flowed<BR><BR>Hi Ji Qing,<BR><BR>Write out 
the energies at every step and check which of the energy<BR>contributions goes 
astray. This may give you a good clue to where<BR>you're mistake is. Also, run 
through the topology again, by hand,<BR>checking all bonds, angles, dihedrals, 
exclusion (implicit ones!),<BR>etc.<BR><BR>Good luck,<BR><BR>Tsjerk<BR><BR>On 
3/7/07, kitty ji &lt;</EM><A 
href="mailto:jiqing@iccas.ac.cn"><EM>jiqing@iccas.ac.cn</EM></A><EM>&gt; 
wrote:<BR>&gt;<BR>&gt;<BR>&gt; Thanks Mark.<BR>&gt; But when one polymer chain 
was used, non-bonded interaction has already been<BR>&gt; there. Because one 
chain means C44H156O22.<BR>&gt; After all, I'd like to try more 
chains.<BR>&gt;<BR>&gt; Message: 1<BR>&gt; Date: Tue, 06 Mar 2007 17:01:52 
+1100<BR>&gt; From: Mark Abraham &lt;</EM><A 
href="mailto:Mark.Abraham@anu.edu.au"><EM>Mark.Abraham@anu.edu.au</EM></A><EM>&gt;<BR>&gt; 
Subject: Re: [gmx-users] vacuum right but melt collapse<BR>&gt; To: Discussion 
list for GROMACS users &lt;</EM><A 
href="mailto:gmx-users@gromacs.org"><EM>gmx-users@gromacs.org</EM></A><EM>&gt;<BR>&gt; 
Message-ID: &lt;</EM><A 
href="mailto:45ED03D0.6010803@anu.edu.au"><EM>45ED03D0.6010803@anu.edu.au</EM></A><EM>&gt;<BR>&gt; 
Content-Type: text/plain; charset=GB2312<BR>&gt;<BR>&gt; kitty ji wrote:<BR>&gt; 
&gt; Hi GMX user:<BR>&gt; &gt;<BR>&gt; &gt; A model for Poly(vinyl methyl ether) 
was built with OPLS force field (<BR>&gt; &gt; all H ).<BR>&gt; &gt;<BR>&gt; 
&gt;<BR>&gt; &gt;<BR>&gt; &gt; When the model was run alone (a single chain) in 
vacuum, anything goes ok.<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt; 
But with 45 chains box, the system will be collapse with a little bond<BR>&gt; 
&gt; elongates abnormal whatever in NVT or NPT.<BR>&gt; &gt;<BR>&gt; 
&gt;<BR>&gt; &gt;<BR>&gt; &gt; The relax time of pressure and temperature were 
adjust but no helpful.<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt;<BR>&gt; &gt; How 
can I found some possible reason and then fix it ? Any advise will<BR>&gt; &gt; 
be appreciated.<BR>&gt;<BR>&gt; If an isolated molecule doesn't do anything 
strange, and two more more<BR>&gt; do, then you have a problem with 
inter-molecular interactions, i.e. the<BR>&gt; non-bonded ones. Try with just 2 
chains to verify this... and look<BR>&gt; carefully at the atom types and the 
non-bonded parameters for them.<BR>&gt;<BR>&gt; 
Mark<BR>&gt;<BR>&gt;<BR>&gt;<BR>&gt; 
*************************************************<BR>&gt; Ji Qing<BR>&gt; 
Institute of Chemistry, Chinese Academy of Sciences<BR>&gt; Tel: 
0086-10-62562894&nbsp; £¬82618423<BR>&gt; 
*************************************************<BR>&gt; 
_______________________________________________<BR>&gt; gmx-users mailing 
list&nbsp;&nbsp;&nbsp; </EM><A 
href="mailto:gmx-users@gromacs.org"><EM>gmx-users@gromacs.org</EM></A><BR><EM>&gt; 
</EM><A 
href="http://www.gromacs.org/mailman/listinfo/gmx-users"><EM>http://www.gromacs.org/mailman/listinfo/gmx-users</EM></A><BR><EM>&gt; 
Please don't post (un)subscribe requests to the list. Use the<BR>&gt; www 
interface or send it to </EM><A 
href="mailto:gmx-users-request@gromacs.org"><EM>gmx-users-request@gromacs.org</EM></A><EM>.<BR>&gt; 
Can't post? Read<BR>&gt; </EM><A 
href="http://www.gromacs.org/mailing_lists/users.php"><EM>http://www.gromacs.org/mailing_lists/users.php</EM></A><BR><EM>&gt;<BR><BR><BR>-- 
<BR>Tsjerk A. Wassenaar, Ph.D.<BR>Junior UD (post-doc)<BR>Biomolecular NMR, 
Bijvoet Center<BR>Utrecht University<BR>Padualaan 8<BR>3584 CH Utrecht<BR>The 
Netherlands<BR>P: +31-30-2539931<BR>F: 
+31-30-2537623<BR></EM></DIV></BODY></HTML>