Thanks Justin, vdwradii.dat suggestion worked :)<br><br><br><br><div class="gmail_quote">On Mon, Jun 13, 2011 at 11:05 AM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im"><br>
<br>
shivangi nangia wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hello dear gmx-users,<br>
<br>
I am trying to create a system which has 2,5-dihydobenzoic acid (DHB) as the solvent ( with a positively charged protein and few DHB anions) in a box.<br>
<br>
I have created the .itp and .gro file of DHB using PRODRG.<br>
<br>
</blockquote>
<br></div>
Unrelated advice - the topology is probably useless unless you've corrected the charges. See the paper linked from:<br>
<br>
<a href="http://www.gromacs.org/Downloads/Related_Software/PRODRG#Tips" target="_blank">http://www.gromacs.org/Downloads/Related_Software/PRODRG#Tips</a><div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
I am able to create the system, but when I try to energy minimize the system I get the message:<br>
<br>
Steepest Descents:<br>
Tolerance (Fmax) = 1.00000e+01<br>
Number of steps = 50000<br>
Step= 14, Dmax= 1.2e-06 nm, Epot= 9.04353e+19 Fmax= inf, atom= 2781<br>
Stepsize too small, or no change in energy.<br>
Converged to machine precision,<br>
but not to the requested precision Fmax < 10<br>
<br>
Double precision normally gives you higher accuracy.<br>
<br>
writing lowest energy coordinates.<br>
<br>
Back Off! I just backed up em.gro to ./#em.gro.1#<br>
<br>
Steepest Descents converged to machine precision in 15 steps,<br>
but did not reach the requested Fmax < 10.<br>
Potential Energy = 9.0435262e+19<br>
Maximum force = inf on atom 2781<br>
Norm of force = inf<br>
<br>
<br>
The reason behind this was that DHB molecules were on top of each other, highly dense system ( viewed in VMD)<br>
So, to reduce the density of DHB molecules, I tried using-vdwd option of genbox (changed from default value of 0.105 to 0.5)<br>
But, I got almost same number of DHB molecules using genbox.<br>
<br>
</blockquote>
<br></div>
The -vdwd option has no effect unless you're using atom names that are not found in vdwradii.dat (see the genbox help description). Your DHB molecule should be unaffected by the use of -vdwd.<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
What is the best solution to this problem as I understand, reducing the density of DHB molecules around the protein as a solvent or there is some other problem?<br>
<br>
</blockquote>
<br></div>
You need to generate a configuration in which molecules aren't overlapping. There are several ways to accomplish this. You can manually place molecules within a box at a desired location with editconf -center, otherwise genbox -ci -nmol. If the locations of the molecules are not what you desire with genbox, change the value of -seed and try again.<br>
<br>
-Justin<br>
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | <a href="tel:%28540%29%20231-9080" value="+15402319080" target="_blank">(540) 231-9080</a><br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br><font color="#888888">
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<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">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</font></blockquote></div><br>