shivangi nangia wrote:
Hello dear gmx-users,

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.

I have created the .itp and .gro file of DHB using PRODRG.


Unrelated advice - the topology is probably useless unless you've corrected the charges. See the paper linked from:

http://www.gromacs.org/Downloads/Related_Software/PRODRG#Tips

I am able to create the system, but when I try to energy minimize the system I get the message:

Steepest Descents:
   Tolerance (Fmax)   =  1.00000e+01
   Number of steps    =        50000
Step= 14, Dmax= 1.2e-06 nm, Epot= 9.04353e+19 Fmax= inf, atom= 2781
Stepsize too small, or no change in energy.
Converged to machine precision,
but not to the requested precision Fmax < 10

Double precision normally gives you higher accuracy.

writing lowest energy coordinates.

Back Off! I just backed up em.gro to ./#em.gro.1#

Steepest Descents converged to machine precision in 15 steps,
but did not reach the requested Fmax < 10.
Potential Energy  =  9.0435262e+19
Maximum force     =            inf on atom 2781
Norm of force     =            inf


The reason behind this was that DHB molecules were on top of each other, highly dense system ( viewed in VMD) 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)
But, I got almost same number of DHB molecules using genbox.


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.


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?


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.

-Justin

--
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to