<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman, new york, times, serif;font-size:12pt"><DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">The end structure is the same as the start because gromacs cannot find a lower energy structure than the initial one. This indicates something severely wrong, either with the topology or starting structure. You could get a better idea of what's going on by telling gromacs to output structures every step in the minimisation, and examining these so you can actually observe what is happening to your ATP.<BR>I've occasionally seen similar results from minimisation, but only after doing something REALLY bad, like precisely overlaying two molecules on the same coordinates. Debugging something like this is probably something only you are going to be able to do, as there are just so many potential ways a
suitably, uh, imaginative user can mess things up. <BR></DIV>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">----- Original Message ----<BR>From: MoJie Duan <mjduan@smail.hust.edu.cn><BR>To: gmx-users@gromacs.org<BR>Sent: Friday, August 17, 2007 11:15:40 AM<BR>Subject: Re: [gmx-users] The energy minimization....<BR><BR>Hi, Mark:<BR>I have done the energy minimization and simulation of ATP in vacuum individual ( Maybe you have suggested me to do this yesterday, but actually I did not understand it, and just do this in solution).<BR>There are following problems:<BR>1. in the energy minimization, the potential energy is positive and in 14th step, it's potential energy is "nan". But the ".gro" outfile of "mdrun" is just the same as the original file (i.e. the .gro file before minimization), the return messages of mdrun is (there are not any warning):<BR><BR>------------------------------<BR>Getting Loaded...<BR>Reading file ATP.2_min.tpr, VERSION 3.3.1 (single
precision)<BR>Loaded with Money<BR><BR><BR>Back Off! I just backed up ATP.2_minrun.edr to ./#ATP.2_minrun.edr.1#<BR>Steepest Descents:<BR> Tolerance (Fmax) = 1.00000e+00<BR> Number of steps = & #160; 200<BR>Step= 0, Dmax= 1.0e-02 nm, Epot= 1.06678e+05 Fmax= 3.63368e+06, atom= 26<BR>Step= 14, Dmax= 1.2e-06 nm, Epot= nan Fmax= 3.63313e+06, <BR> ^^^^^^^^^^^^^^^^^^^^^^^<BR>atom= 26<BR>Stepsize too small, or no change in energy.<BR>Converged to machine precision,<BR>but not to the requested precision Fmax < 1<BR><BR>Double precision normally gives you higher accuracy.<BR><BR>writing lowest
energy coordinates.<BR><BR>Back Off! I just backed up ATP.2_minrun.gro to ./#ATP.2_minrun.gro.1#<BR><BR>Steepest Descents converged to machine precision in 15 steps,<BR>but did not reach the requested Fmax < 1.<BR>Potential Energy = 1.0667836e+05<BR>Maximum for ce = 3.6336775e+06 on atom 26<BR>Norm of force = nan<BR>___________________________________________<BR><BR>2. in the full MD simulation, the "warning" coming, the messages is:<BR><BR>--------------------------------------------<BR>Step 0, time 0 (ps) LINCS WARNING<BR>relative constraint deviation after LINCS:<BR>max 0.881911 (between atoms 27 and 28) rms nan<BR>bonds that rotated more than 30 degrees:<BR>^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^<BR> atom 1 atom 2 angle previous, current, constraint length<BR>
27 28 90.0 0.1610 0.3030 0.1610<BR>Wrote pdb files with previous and current coordinates<BR>step 2490, remaining runtime: 0 s & #160; Writing final coordinates.<BR><BR>Back Off! I just backed up ATP.2_runout.gro to ./#ATP.2_runout.gro.1#<BR>step 2500, remaining runtime: 0 s <BR>__________________________________________________<BR><BR>And in the ".gro" file after this step, the coordinates of all atoms are "nan". So it means there are crash in the structure? Is the crash between the atom 27 and 28? How to modify the structure file make it normal?<BR>Thank you very much!<BR><BR>Duan<BR><BR><BR>>MoJie Duan wrote:<BR>> > >So look at your
structures like I said last time! I'm not her e to give<BR>> > >my valuable time giving free advice in order to have it ignored...<BR>> > Thank you very much for your kindness and patience. Maybe sometimes my <BR>> > questions seems to be silly and boring, my knowledge about GROMACS is <BR>> >really lack, sorry.....<BR>> <BR>> >Actually, I really cannot understand what you said yesterday. Did you <BR>> >mean is there any difference between the atom coordinates of ATP before- <BR>> >and after- minimization?<BR><BR><BR>> There would normally be some differences visible. If your topology was<BR>> badly broken, then you would usually see where it was broken.<BR><BR>> >I found there are not any difference between <BR>> >these two structure.<BR>> >(There are also not any obvious collision between atoms of ATP when <BR>> >represent it by Rasmol)<BR><BR>> OK, so that means your
structure is in a flat area of the potential<BR>> surface defined by your topology. If the topology is sound, then you're<BR>> in business.<BR><BR>> My first recommendation was to minimize and/or equilibrate these<BR>> structures on their own, and now I suggest doing them also in solvent.<BR>> This will help you eliminate sources of problems and guide you to <BR>> what the real problem is. Divide and conquer...<BR><BR>> That's fine, then.<BR><BR>> OK, so here's your problem. Work out what's breaking and why. Read the<BR>> error messages and look at the structures. Understand what each of > your .mdp file options does.<BR><BR>Mark<BR><BR></DIV><BR></DIV></div><br>
<hr size=1>Fussy? Opinionated? Impossible to please? Perfect. <a href="http://us.rd.yahoo.com/evt=48516/*http://surveylink.yahoo.com/gmrs/yahoo_panel_invite.asp?a=7 ">Join Yahoo!'s user panel</a> and lay it on us.
</body></html>