<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div><span>Dear MArkWow, That was kind of you.</span></div><div><span>I appreciate you.</span></div><div><span>Ok, I will try it.</span></div><div><span>Thanks,</span></div><div><span>Shima<br></span></div><div><br></div> <div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">From:</span></b> Mark Abraham <Mark.Abraham@anu.edu.au><br> <b><span style="font-weight: bold;">To:</span></b> Discussion list for GROMACS users <gmx-users@gromacs.org> <br> <b><span style="font-weight: bold;">Sent:</span></b> Sunday, April 29, 2012 2:11 PM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [gmx-users] Make an index file<br>
</font> </div> <br>
<div id="yiv131850763">
<div>
On 29/04/2012 7:32 PM, Shima Arasteh wrote:
<blockquote type="cite">
<div style="color:rgb(0, 0, 0);background-color:rgb(255, 255,
255);font-family:arial, helvetica, sans-serif;font-size:12pt;">
<div><span>Actually, I am confused somehow. <br>
</span></div>
<div><span>I want to equilibrate the system. It contains popc
and water. To equilibrate it, using .mdp file I use the
nvt.mdp file as below.</span></div>
<div><span><br>
</span></div>
<div><span>title = NVT equilibration for POPC<br>
define = -DPOSRES ; position restrain the protein<br>
; Run parameters<br>
integrator = md ; leap-frog integrator<br>
nsteps = 50000 ; 2 * 50000 = 100 ps<br>
dt = 0.002 ; 2 fs<br>
; Output control<br>
nstxout = 100 ; save coordinates every 0.2 ps<br>
nstvout = 100 ; save velocities every 0.2 ps<br>
nstenergy = 100 ; save energies every 0.2 ps<br>
nstlog = 100 ; update log file every 0.2 ps<br>
; Bond parameters<br>
continuation = no ; first dynamics run<br>
constraint_algorithm = lincs ; holonomic constraints <br>
constraints = all-bonds ; all bonds (even
heavy atom-H bonds) constrained<br>
lincs_iter = 1 ; accuracy of LINCS<br>
lincs_order = 4 ; also related to
accuracy<br>
; Neighborsearching<br>
ns_type = grid ; search neighboring grid cels<br>
nstlist = 5 ; 10 fs<br>
rlist = 1.2 ; short-range neighborlist cutoff
(in nm)<br>
rcoulomb = 1.2 ; short-range electrostatic cutoff
(in nm)<br>
rvdw = 1.2 ; short-range van der Waals cutoff
(in nm)<br>
; Electrostatics<br>
coulombtype = PME ; Particle Mesh Ewald for
long-range electrostatics<br>
pme_order = 4 ; cubic interpolation<br>
fourierspacing = 0.16 ; grid spacing for FFT<br>
; Temperature coupling is on<br>
tcoupl = V-rescale ; modified
Berendsen thermostat<br>
tc-grps = POPC SOL ; two coupling groups - more
accurate<br>
</span></div>
</div>
</blockquote>
<br>
Here you instruct grompp to expect the POPC group...<br>
<br>
<blockquote type="cite">
<div style="color:rgb(0, 0, 0);background-color:rgb(255, 255,
255);font-family:arial, helvetica, sans-serif;font-size:12pt;">
<div><span>tau_t = 0.1 0.1 ; time
constant, in ps<br>
ref_t = 323 323 ; reference
temperature, one for each group, in K<br>
; Pressure coupling is off<br>
pcoupl = no ; no pressure coupling in NVT<br>
; Periodic boundary conditions<br>
pbc = xyz ; 3-D PBC<br>
; Dispersion correction<br>
DispCorr = EnerPres ; account for cut-off vdW scheme<br>
; Velocity generation<br>
gen_vel = yes ; assign velocities from Maxwell
distribution<br>
gen_temp = 323 ; temperature for Maxwell
distribution<br>
gen_seed = -1 ; generate a random seed<br>
<br>
<br>
</span></div>
<div><span><br>
</span></div>
<div><span> I enter this command:</span></div>
<div><span># grompp -f nvt.mdp -c em.gro -p popc.top -o nvt.tpr</span></div>
<div><br>
<span></span></div>
<div><span>Getting this </span><span>fatal error:<br>
Group POPC not found in index file.<br>
Group names must match either [moleculetype] names<br>
or custom index group names,in which case you<br>
must supply an index file to the '-n' option of grompp. <br>
</span></div>
<div><br>
<span></span></div>
<div><span>I don't </span><span>need </span><span>index file,
do I? What is the problem with my .mdp file?</span></div>
</div>
</blockquote>
<br>
... but POPC is not defined implicitly by your .top+.gro
combination, nor explicitly in an index file supplied with grompp
-n. AFAIK if your [moleculetype] is named POPC then it is all OK,
but maybe the fact that your .gro file has the POPC "residues" named
POP is confusing that mechanism. As you can see with your earlier
example with make_ndx, the POP group is defined implicitly by the
residue names in your .gro file, so if you use that output index
file and reference POP in your .mdp file, all will be good. You
might get away without using that index file, but I am not sure
about this. Simpler still is to use Water and non-Water as your
T-coupling groups, but this will need changing if down the road you
need Protein as well (or such).<br>
<br>
Mark<br>
</div>
</div><br>-- <br>gmx-users mailing list <a ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">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 <br>www interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" href="mailto:gmx-users-request@gromacs.org">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><br> </div> </div> </div></body></html>