[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