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