I notice, on the attached output, that you haven't neutralized the sistem. A neutral sistem is required for ewald sum, and it is artificially neutralized with uniform plasma (see your output). Maybe this can cause a problem. You problem could also depend on the initial structure. If it is a PDB file or a part of a pdb file, there should be no problems, but if it was drawn by hand, some artifacts may be present that make the energy high and that can't be solved by minimization.
Stefano Pieraccini
Content-type: multipart/alternative;
boundary="Boundary_(ID_oOhDSNnNENFnNM4uD5RcJQ)"
--Boundary_(ID_oOhDSNnNENFnNM4uD5RcJQ)
Content-type: text/plain; charset=us-ascii
Content-transfer-encoding: 7BIT
hello amber,
I am performing Implicit solvent md on small proteins using 4r dielectric, restraing the helices ans sheets. The input file for the initial minimization step is:
&cntrl
imin=1,
ntx=1,ntb=0,
ntr=1,dielc=4,cut=300.0,
maxcyc=1000, ncyc=500,
&end
&ewald
eedmeth=5,
&end
restraint 1
25.0
RES 10 19
RES 25 37
RES 42 55
END
END
My problem is that the total energy at the start of the minimization is positive and after the minimization (1000 steps) is positive. Is there someting which i am doing wrong, because of which total energy is positive. I am attaching the output file.
Thanks.
---------------------------------
Do you Yahoo!?
Take Yahoo! Mail with you! Get it on your mobile phone.
--Boundary_(ID_oOhDSNnNENFnNM4uD5RcJQ)
Content-type: text/html; charset=us-ascii
Content-transfer-encoding: 7BIT
<DIV>
<DIV>hello amber,</DIV>
<DIV>I am performing Implicit solvent md on small proteins using 4r dielectric, restraing the helices ans sheets. The input file for the initial minimization step is:</DIV>
<DIV>&cntrl<BR> imin=1,<BR> ntx=1,ntb=0,<BR> ntr=1,dielc=4,cut=300.0,<BR> maxcyc=1000, ncyc=500,<BR> &end<BR> &ewald<BR> eedmeth=5,<BR> &end<BR>restraint 1<BR> 25.0<BR>RES 10 19<BR>RES 25 37<BR>RES 42 55<BR>END<BR>END</DIV>
<DIV>My problem is that the total energy at the start of the minimization is positive and after the minimization (1000 steps) is positive. Is there someting which i am doing wrong, because of which total energy is positive. I am attaching the output file.</DIV>
<DIV>Thanks.</DIV></DIV><p>
<hr size=1>Do you Yahoo!?<br>
<a href="
http://us.rd.yahoo.com/mail_us/taglines/mobile/*http://mobile.yahoo.com/maildemo">Take Yahoo! Mail with you!</a> Get it on your mobile phone.
--Boundary_(ID_oOhDSNnNENFnNM4uD5RcJQ)--
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Dec 01 2004 - 09:53:00 PST
