Previous Thread
Next Thread
Print Thread
#812 01/06/04 08:22 PM
Joined: Oct 2003
Posts: 6
jfm17 Offline OP
Forum Member
OP Offline
Forum Member
Joined: Oct 2003
Posts: 6
I am attempting to reproduce a GROMOS force field that contains VDW terms given in the form (A/r**12) - (B/r**6). One atom type has B=0, giving a repulsive VDW term only. CHARMM has the form e(ij) [ (Rmin/r)**12 - 2(Rmin/r)**6 ], so usually it is not a problem to convert the A and B into epsilon(i,j) and Rmin, but how can I reproduce the term where B=0?

Thanks,
James Matthews

Joined: Sep 2003
Posts: 250
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 250
Interesting question:
One way can be to set the attractive term = 0*(rmin/r)^6, which means that one can use the USERE-hookup in CHARMM for a user-energy defined potential. Have just begun to learn more about this facility, but see usage.doc


Moderated by  alex, lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.5
(Release build 20201027)
Responsive Width:

PHP: 7.3.31-1~deb10u3 Page Time: 0.005s Queries: 18 (0.004s) Memory: 0.7276 MB (Peak: 0.7673 MB) Data Comp: Off Server Time: 2023-06-09 11:17:37 UTC
Valid HTML 5 and Valid CSS