I have construct two table files: table_b0.xvg for A-A bond, table_b1.xvg for B-B bond. and generated the bonds part of the topology file like this: [ bonds ] ; ai aj funct table_number k 1 2 8 0 1.000 2 3 8 0 1.000 ...... 55 56 8 1 1.000 56 57 8 2 1.000 ...... as force constant k is a function, so I set k = 1.000 and combine it with f.
Am I correct? Thank you very much. I really appreciate your help, Justin. -- * Jia Haitao State Key Laboratory of Polymer Physics and Chemistry Changchun Institute of Applied Chemistry Chinese Academy of Sciences Changchun 130022, CHINA *
-- 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