Hello, I am relatively new to gromacs and I really would appreciate help with the following question.
I have used a coarse-grained model to represent a DNA fragment with approximately 50 nm of length. I have run a Langevin dynamics simulation of DNA with ions NA+ and CL- to study their distribution of around DNA. During simulation I used the option to remove the center of mass motion (default option) because when there are no external forces the com linear moment must be conserved. Now, I want study the redistribution of ions around DNA when an external electric field is applied. My doubt is the following: when external electric field is applied, should I remove the center of mass motion? I thought that the answer was no, then I tried the following option: comm_mode none but half of DNA goes out of the box during simulation (I see this with ngmx and in out.gro) and nothing enters the another side. I know gromacs use pbc to calculate the interactions between particles, but I wonder if g_rdf program considers that part of DNA is out the box to calculate the radial function distribution of ions around DNA. Thank. Pablo
_______________________________________________ gmx-users mailing list gmx-users@gromacs.org http://www.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