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 [email protected]
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 [EMAIL PROTECTED]
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

_______________________________________________
gmx-users mailing list    [email protected]
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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to