It seems that you can modify the following line in the function
*mk_nbfp() in force.c and achieve what you want.
C6(nbfp,atnr,i,j) = idef->iparams[k].lj.c6; (it is in line 117
of force.c in version 3.3)
and change it to
C6(nbfp,atnr,i,j) = 0;
This will set the attractive interaction between atom i & j to 0.
pb.
On Jun 30, 2006, at 4:14 AM, Alessandro Mattozzi wrote:
Hi All
I simulate PE using AUA. I would like to switch the
attractive term of the vdW force. How can I do it?
Regards
Alessandro Mattozzi
M.Phil., Ph.D. student
Dept. of Fibre and Polymer Technology
Royal Institute of Technology
Stockholm, Sweden
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-request@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
--
Pradip K. Biswas, PhD.
Research Associate, Department of Chemistry;
Cleveland State University, Ohio-44115
Phone: 1-216-875-9723
http://comppsi.csuohio.edu/groups/people/biswas.html