[CP2K-user] [CP2K:21590] Clarification of R_CUTOFF in VDW_POTENTIAL

emerson p l emersonp90 at gmail.com
Thu Jun 26 19:41:04 UTC 2025



Dear cp2k community, I hope you’re all doing well. 

I’m writing to ask for clarification on the choice of the cutoff radius 
R_CUTOFF​ parameter when using periodic boundary conditions in 
VDW_POTENTIAL in AIMD and geometry/cell optimizations:

      &VDW_POTENTIAL
         POTENTIAL_TYPE PAIR_POTENTIAL
         &PAIR_POTENTIAL
            PARAMETER_FILE_NAME dftd3.dat
            TYPE DFTD3
            REFERENCE_FUNCTIONAL PBE
            R_CUTOFF [angstrom] 10
         &END
      &END VDW_POTENTIAL

Suppose I have a cubic simulation box with side length L = 5 angstrom, and 
I (mistakenly) set R_CUTOFF​ = 10 angstrom. Would setting RCUTOFF​ = 10 
angstrom in a 5 angstrom cubic box cause the cutoff sphere to extend beyond 
the primary cell and overlap with periodic images, resulting in unphysical, 
multiply-counted interactions?

My understanding from classical MD is that we must obey the minimum-image 
convention, which requires R_CUTOFF​ ≤ 2 Lmin​​. I don't know if R_CUTOFF 
has this physical meaning of classical MD.

Thank you very much for your time and insights. 

Best,

Emerson


-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/0afee8cd-252c-4c9e-b0c1-e2142551930bn%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20250626/780d824d/attachment.htm>


More information about the CP2K-user mailing list