> Date: Thu, 19 Mar 2009 13:52:09 +0800 > From: huifang liu <huifangliu1...@gmail.com> > Subject: [gmx-users] how to get the the force plot after pulling > dynamic simulation using GMX-3.3 > To: gmx-users@gromacs.org > Message-ID: > <e523c6be0903182252i114666bam4bbe591a10a30...@mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > I am recently doing a SMD. I think i did it very well. But i don't know how > to get the force plot. Is there some one can give me a suggestion? > Thank you very much. > At the end of the SMD you should have obtained a .pdo file. In that file the time and the coordinates of the reference- (R), pull-group (P) and of the spring (S) are written. If i remember correctly it was: t, Rx, Ry, Rz, Px, Sx, Py, Sy, Pz, Sz. If you have more pullgroups then you have first the x-coordinates of P and S, then y-coordinates ... (look in the manual or the header of the .pdo; or for t=0 you should know the position of all particals and springs and can calculate a little bit to see which entry corresponds to which group). For the force-calculation: The force is the deflection of the spring times the force-constant of the spring -> Fx = k(Sx-Px) ... for the other directions it's analog. Then you can calculate the magnitude of F and plot it against the time (if you are only interested in the magnitude of force). Or you can calculate the force in the pulling direction and plot it against time.
Hope this helps. Thomas > Huifang > > -- > Huifang Liu (Ph.D. Student) > School of Pharmacy > Fudan University > > 138 Yi Xue Yuan Rd. Tel: (86-21)54237419 (O) > Shanghai, China, 200032 Cell phone: +86-13764669357 > E-mail: huifangliu1...@gmail.com Fax: (86-21)54237264 > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > http://www.gromacs.org/pipermail/gmx-users/attachments/20090319/2f249dff/attachment-0001.html > > ------------------------------ > _______________________________________________ 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
