<div dir="ltr">Dear all,<div><br></div><div>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.</div><div><br></div><div>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</div><div><br></div><div>I have attached an example input, output and .png for the system with ionic and covalent radii.</div><div><br></div><div><div>Any guidance or insight to resolve this problem would be greatly appreciated.</div><div><br></div><div>Thank you for your help and time.</div><div><br></div><div>Regards,</div><div>Chris</div></div><div><br></div><div><br></div><div><br></div><div>Ionic radii (Shannon, R. D. (1976))</div><div>Mg 0.72 </div><div>O 1.40</div><div><br></div><div>Covalent radii (Pyykkö, P., & Atsumi, M. (2009)</div><div><div>Mg 1.39</div><div>O 0.63</div></div><div><br></div><div><div tabindex="0" class="F0XO1GC-nb-P" style="line-height: normal;"><div style="overflow: auto;"><div style="max-height: 10000px;"><div dir="ltr"><br></div></div></div></div></div></div>