Hi Nisha,
you say "And it runs minimization with bad contact errors", but I don't see any indication of that error message here.
Is it perhaps that your EM exits early and then your MD throws an error? In any event, please post complete error information
as output by gromacs.
Chris.
-- original message --
I have a question regarding the simulations that I am running. I tried
creating a pdb file with one methane and water molecules, but the
system keeps showing 'bad contact' error when I run the minimization
step and I don't understand what I am doing wrong.
This is my methane file :
TITLE methane in water
CRYST1 33.010 33.010 33.010 90.00 90.00 90.00 P 1 1
MODEL 1
ATOM 1 CH4 CH4 1 19.520 25.090 16.410 1.00 0.00
So I use this file and run editconf -f methane1.pdb -o box.pdb -bt
cubic. This creates a
box with volume of around 35 nm^3
Then I add water molecules using :
genbox -cp box.pdb -cs spc216 -o methwa.pdb -p topol.top
This adds around 1191 water molecules to the existing methane file. So
my box is now
(33.010 33.010 33.010 with total 1192 molecules). Then I run
minimization on this system
using minim.mdp as follows:
integrator = steep
emtol = 1.0
nsteps = 100
nstenergy = 10
nstxtcout = 10
xtc_grps = CH4 SOL
energygrps = CH4 SOL
nstlist = 5
ns_type = simple
rlist = 1.0
coulombtype = cut-off
rcoulomb = 1.0
rvdw = 1.0
constraints = none
pbc = no
grompp -f minim.mdp -c methwa.pdb -p topol.top -v
mdrun -v
And it runs minimization with bad contact errors and this is the result I get:
Steepest Descents converged to machine precision in 32 steps,
but did not reach the requested Fmax < 100.
Potential Energy = -4.9481789e+04
Maximum force = 3.2827865e+06 on atom 500
Norm of force = 6.9951227e+04
Any suggestions?
Nisha.P
--
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/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/mailing_lists/users.php
