Previous Thread
Next Thread
Print Thread
#32490 08/25/13 05:12 AM
Joined: Feb 2004
Posts: 68
G
Forum Member
OP Offline
Forum Member
G
Joined: Feb 2004
Posts: 68
I wonder what are recommended cutoff parameters for non-bonded interactions (in particular VDW) when using PME.

In nbonds.doc I read:

For no-cutoff periodic systems:

NBONDS ATOM EWALD PMEWALD KAPPA 0.32 -
FFTX 32 FFTY 32 FFTZ 32 ORDER 6 -
CUTNB 12.0 CTOFNB 10.0 VDW VSHIFT

So, here the recommendation is 9/10/12 for CTONNB/CTOFNB/CUTNB

However, in file c26test/pme_p21.inp I find:

dyna verlet ... -
cutnb 12.0 ctofnb 10. cdie eps 1. -
ctonnb 9.99 vshift cutim 14.0 imgfrq 25 wmin .5 -
ewald pmewald

So, here they have 9.99/10/12 which corresponds to an almost abrupt truncation of VDW.

I would appreciate any advice on this. Thanks!

Last edited by gianluca; 08/25/13 05:15 AM.
gianluca #32491 08/25/13 02:13 PM
Joined: Sep 2003
Posts: 8,605
Likes: 23
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,605
Likes: 23
Do not use the test cases as examples of good usage.

With VSHIFT, the value of CTONNB is irrelevant, as it is not used; only switching functions use this value.

For the VDW terms, we use and recommend

VATOM VFSWITCH CUTNB 14. CTOFNB 12. CTONNB 10.

which is the default in most "all36" parameter sets.


Rick Venable
computational chemist


Moderated by  BRBrooks, 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~deb10u1 Page Time: 0.007s Queries: 18 (0.004s) Memory: 0.7286 MB (Peak: 0.7690 MB) Data Comp: Off Server Time: 2022-09-27 17:32:06 UTC
Valid HTML 5 and Valid CSS