[CP2K-user] [CP2K:19924] Understanding GAPW gradient implementation: why many terms are not calculated?
zc62
chenzh1993 at gmail.com
Fri Feb 16 21:51:13 UTC 2024
Hi,
Recently I have been trying to derive GAPW nuclear coordinate gradients
(somehow they are not presented in the 1999 and 2000 GAPW papers) and
compare with the actual implementation, but there are a few places I could
not understand.
[image: 截屏2024-02-16 下午3.35.17.png]
Above shows what I have derived for Hartree and Exc (without density matrix
derivative terms). However, the actual implementation does not have any 1c
terms implemented. To be specific, they are:
1. C' derivatives. C' is calculated from the projection, and the overlap
between the projector at A and basis at B certainly will change.
2. For the compensation density (rho0), the integral of g_A^lm gradient
with V_H is calculated in qs_rho0_ggrid.F on the global grid, but the
corresponding local grid integral is not calculated (if I did not miss it).
3. All the local grid integrals involve the gradient of basis functions are
not calculated (at least not in hartree_local_methods.F).
According to my numerical tests, the analytic gradients do agree with
numerical differentiations, but I do feel that I am unable to correctly
understand why these terms should be zero. Can you help me develop a proper
understanding for this?
Thanks!
--
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 on the web visit https://groups.google.com/d/msgid/cp2k/28b358d5-3ccb-4e54-8649-f5ce588c54a2n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240216/7467b08e/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ??2024-02-16 ??3.35.17.png
Type: image/png
Size: 79321 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240216/7467b08e/attachment-0001.png>
More information about the CP2K-user
mailing list