On 12/12/12 8:14 AM, Shima Arasteh wrote:
I write the exact commands up to grompp here, not to miss anything;
1. pdb2gmx -f dimer-rotated.pdb -o dimer-processed.pdb -ter -water tip3p
2. grompp -f minim.mdp -c md-popc.gro -p topol_popc.top -o em.tpr
3. trjconv -s em.tpr -f md-popc.gro -o popc-whole.gro -pbc mol -ur compact
4. editconf -f dimer-processed.pdb -o dimer-newbox.pdb -c -box 8.78486
8.51308 12.40746
5. editconf -f dimer-newbox.pdb -o dimer-newbox.gro
6. cat dimer-newbox.gro popc-whole.gro > system.gro
7. genrestr -f dimer-newbox.gro -o strong_posre.itp -fc 100000 100000 100000
8. perl inflategro.pl system.gro 6 POPC 14 system_inflated.gro 5 area.dat
9. grompp -f minim.mdp -c system_inflated.gro -p topol.top -o em.tpr
10. mdrun -deffnm em
The positive potential I have got here, is before the "solvate with water
step". I don't have sol molecules in my .top file, however I see sol molecules in
system.gro . I visualized system_inflated.gro , I don't see any water molecules.
Now, I'm confused if this potential energy of a system with sol molecules or
these are just droplets as you said?
Your results are fine. InflateGRO removes water molecules (as stated in its
documentation). You should only be concerned if the potential energy of the
final system (with water and ions, after all the rounds of shrinking are done)
is positive. The energies of these intermediate steps mean very little.
-Justin
--
========================================
Justin A. Lemkul, Ph.D.
Research Scientist
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