<div>Hi all, </div><div><br /></div><div>I'm using GAPW with pseudos for a system with antimony as I need a high precision on the atomic forces - everything works as intended, but with the default parameters, Sb has no hard part: </div><div><br /></div><div>&KIND Sb<br /> BASIS_SET ORB TZVP-MOLOPT-PBE-GTH-q5<br /> POTENTIAL GTH-GGA-q5<br />&END KIND</div><div><br /></div><div>From the output: </div><div><br /></div><div> QS| GAPW| At least one kind is NOT PAW, i.e. it has only soft AO<br /> QS| GAPW| The NOT PAW atoms are treated fully GPW</div><div><br /></div><div>This might not be a problem (since I get fairly well-converged forces), but the quantities I usually check and which should converge to zero are very much non-zero here: </div><div><br /></div><div> Total Rho_soft + Rho1_hard - Rho1_soft -1440.0000494132<br /><b> Total charge density (r-space): -1290.0000494132</b><br /> Total Rho_soft + Rho0_soft (g-space): -0.0000531123</div><div><br /></div><div>I know it's due to the hard/soft separation, since when I use a smaller pseudo for Sb (such as q23), there is a hard part, and these quantities converge to zero, but I wanted to ask you if this is intended or important. </div><div><br /></div><div>Otherwise, would it make sense to tune HARD_EXP_RADIUS in KIND? </div><div><br /></div><div>Also, is there a reason why there is a q15 pseudo available for Sb (in POTENTIAL_UZH), but no associated basis set? </div><div><br /></div><div>thanks, </div><div>Arthur</div><div><br /></div>
<p></p>
-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/f819c101-df6e-4316-a967-96df20b6b9a9n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/f819c101-df6e-4316-a967-96df20b6b9a9n%40googlegroups.com</a>.<br />