Hello,

I am trying to restrain Na+ to a specific position (0.487, 1.620, 1.620) of
my box of dimensions 3.2418 X 3.2418 X 3.2418 nm.  The box is also full of
253 THF molecules.  I added the following to the bottom of my .itp file:

#ifdef POSRES
#include "posre_Na+.itp"
#endif

I wrote the posre_Na+.itp file myself, it just includes the position
restraints for Na+, which is atom 1 in my .gro file, and my posre_Na+.itp
file deals with atom 1.  The force constants are all 1000.

Under the .mdp file, I also made sure that define = -DPOSRES.

I get no errors when things run, but I do get my Na+ moving around.  After I
do energy minimization, Na+ ends up at the position (0.373, 1.504, 1.226).
I also do not notice anything about position restraints appearing in my .log
file afterwards.

I really appreciate the help.

Thank you,
Jenny
-- 
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

Reply via email to