[CP2K-user] CDFT weight function breaks for small ATOMIC_RADII

[Chris Ahart]

Dear all,

I have found that when ionic radii are used for a bulk MgO crystal the 
Becke weight function appears to be calculated incorrectly (see attached 
.png). When using covalent radii however the weight function is calculated 
correctly.

More generally, it appears that the weight function breaks when the 
constrained atoms have an atomic radius smaller than the unconstrained 
atoms. I have confirmed this behaviour for two different structures on two 
different machines using CP2K 7.1

I have attached an example input, output and .png for the system with ionic 
and covalent radii.

Any guidance or insight to resolve this problem would be greatly 
appreciated.

Thank you for your help and time.

Regards,
Chris



Ionic radii (Shannon, R. D. (1976))
Mg 0.72 
O 1.40

Covalent radii (Pyykkö, P., & Atsumi, M. (2009)
Mg 1.39
O 0.63


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200330/a1941b8c/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: covalent.png
Type: image/png
Size: 1237510 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200330/a1941b8c/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ionic.png
Type: image/png
Size: 1939047 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200330/a1941b8c/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: covalent.inp
Type: chemical/x-gamess-input
Size: 3024 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200330/a1941b8c/attachment.inp>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: covalent.out
Type: application/octet-stream
Size: 132446 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200330/a1941b8c/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ionic.inp
Type: chemical/x-gamess-input
Size: 3024 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200330/a1941b8c/attachment-0001.inp>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ionic.out
Type: application/octet-stream
Size: 127584 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200330/a1941b8c/attachment-0001.obj>