Hi Dhananjay,<br><br>You have some serious problems with your starting structure. You may want to try to energy minimize your system before solvation. If that doesn't work, inspect your structure for overlapping atoms or crossing bonds.
<br><br>Cheers,<br><br>Tsjerk<br><br><div><span class="gmail_quote">On 14 Dec 2005 08:30:18 -0000, <b class="gmail_sendername">Dhananjay joshi</b> <<a href="mailto:joshidc_in1@rediffmail.com">joshidc_in1@rediffmail.com
</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><p>
<br>
Hello Mark,<br>
<br>
I have again started working on CORE_MTB.pdb<br>
<br>
The first step that I took was <br>
<br>
1. pdb2gmx -f CORE_MTB.pdb -o CORE_MTB.gro -p CORE_MTB.top<br>
<br>
<br>
It has generated fatal error as follows<br>
<br>
Fatal error:<br>
Atom HD1 in residue HISB 61 not found in rtp entry with 12 atoms<br>
while sorting atoms. Maybe different protonation state.<br>
Remove this hydrogen or choose a different protonation state.<br>
Option -ignh will ignore all hydrogens in the input.<br>
<br>
Then I used -ignh option and the steps are as follows..<br>
<br>
<br>
2. editconf -f CORE_MTB.gro -o CORE_MTB.gro -d 0.8<br>
<br>
WARNING: No boxtype specified - distance condition applied in each dimension.<br>
If the molecule rotates the actual distance will be smaller. You might want<br>
to use a cubic box instead, or why not try a dodecahedron today?<br>
<br>
3. genbox -cp CORE_MTB.gro -cs -o CORE_MTB_b4em.gro -p CORE_MTB.top<br>
<br>
4.grompp -f em.mdp -po em_out.mdp -c CORE_MTB_b4em.gro -p CORE_MTB.top -o CORE_MTB_em.tpr<br>
<br>
<br>
5.mdrun -s CORE_MTB_em.tpr -o CORE_MTB_em.trr -c CORE_MTB_b4pr.gro -e em.edr -g em.log<br>
<br>
After this I got some message with output as follows<br>
<br>
Steepest Descents converged to machine precision in 15 steps,<br>
but did not reach the requested Fmax < 2000.<br>
Potential Energy = 2.1868101e+17<br>
Maximum force = inf on atom 28916<br>
Norm of force = nan<br>
<br>
Could you please tell me what is wrong ?<br>
<br>
The mdp file used for grompp amd mdrun is as follows..<br>
<br>
<br>
<br>
; User spoel (236)<br>
; Wed Nov 3 17:12:44 1993<br>
; Input file<br>
;<br>
cpp = /usr/bin/cpp<br>
define = -DFLEX_SPC<br>
constraints = none<br>
integrator = steep<br>
nsteps = 100<br>
;<br>
; Energy minimizing stuff<br>
;<br>
emtol = 2000<br>
emstep = 0.01<br>
<br>
nstcomm = 1<br>
ns_type = grid<br>
rlist = 1<br>
rcoulomb = 1.0<br>
rvdw = 1.0<br>
Tcoupl = no<br>
Pcoupl = no<br>
gen_vel = no<br>
<br>
<br>
<br>
On Tue, 13 Dec 2005 Mark Abraham wrote :<br>
>Dhananjay joshi wrote:<br>
>> Yes Mark,<br>
>><br>
>>You are right.<br>
>><br>
>>using "gmx_dump -s CORE_MTB_em.tpr | less" i found large vector of velocities (v) after the large vector of positions (x) and that they are all non-zero.<br>
>><br>
>>what should I do now ?<br>
><br>
>>>see if there is a large vector of velocities (v) after the large vector of positions (x) and that they are all zero. If they aren't, I suspect you haven't actually done as you describe above.<br>
><br>
>grompp is the only thing that can be generating velocities. Check you've done what you think you've done.<br>
><br>
>Mark<br>
>_______________________________________________<br>
>gmx-users mailing list<br>
><a href="mailto:gmx-users@gromacs.org" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">gmx-users@gromacs.org</a><br>
><a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
>Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">gmx-users-request@gromacs.org
</a>.<br>
</p>
<br><br>
<a href="http://adworks.rediff.com/cgi-bin/AdWorks/sigclick.cgi/www.rediff.com/signature-home.htm/1507191490@Middle5?PARTNER=3" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)"><img border="0" hspace="0" vspace="0">
</a>
<br>_______________________________________________<br>gmx-users mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">
http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don't post (un)subscribe requests to the list. Use the<br>www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org">
gmx-users-request@gromacs.org</a>.<br><br></blockquote></div><br><br clear="all"><br>-- <br><br>Tsjerk A. Wassenaar, M.Sc.<br>Groningen Biomolecular Sciences and Biotechnology Institute (GBB)<br>Dept. of Biophysical Chemistry
<br>University of Groningen<br>Nijenborgh 4<br>9747AG Groningen, The Netherlands<br>+31 50 363 4336<br>