Dear all, <div><br /></div><div>I am doing a PBE0 calculation. </div><div>For the primary basis set, I use TZVP-MOLOPT-PBE0-GTH; for the pseudopotential, I use GTH-PBE0; and for the auxiliary basis set, I use FIT3. </div><div><br /></div><div>Using Si as an example, the exact numbers are the following:</div><div><br /></div><div><b>PRIMARY BASIS SET</b></div><div># Si [12s12p4d|3s3p1d]<br /># SOURCE: https://github.com/cp2k/cp2k/raw/ac0226eb549c7ef1ea50d0597d545f29d4c8fc87/data/BASIS_MOLOPT_UZH#L5260<br />Si TZVP-MOLOPT-PBE0-GTH-q4 TZVP-MOLOPT-HYB-GTH-q4<br />1<br />2 0 2 4 3 3 1<br /> 1.213658630823 -0.030815174945 -0.634072471864 -0.542641279521 0.136162907866 0.066271452739 0.848730515959 0.278715935685<br /> 0.428809571621 0.005952863069 0.570670193748 0.150997228423 -0.777588638886 -0.496367537858 0.304840077053 0.788029670430<br /> 0.165906813294 -0.876412895115 0.270751967113 -0.381705514748 -0.613274545168 0.667988337027 -0.181028440213 0.547416077765<br /> 0.042600429814 -0.480536601823 -0.446072867026 0.732830934650 0.026644013324 0.550471564668 0.392374492719 0.040771356636<br /></div><div><br /></div><div><b>PSEUDOPOTENTIAL</b></div><div># Si [10|2s2p]<br /># SOURCE: https://github.com/cp2k/cp2k/raw/ac0226eb549c7ef1ea50d0597d545f29d4c8fc87/data/POTENTIAL_UZH#L6300<br />Si GTH-PBE0-q4 GTH-HYB-q4<br />2 2 0 0<br /> 0.44599315 2 -6.14101714 0.02439420<br /> 2<br /> 0.43467858 2 8.97165588 -2.70628598<br /> 3.49478023<br /> 0.49725736 1 2.43351142<br /></div><div><b>AUXILIARY BASIS SET</b></div><div>Si FIT3<br /> 6<br /> 1 0 0 1 1<br /> 0.09648 1.0<br /> 1 0 0 1 1<br /> 0.24696 1.0<br /> 1 0 0 1 1<br /> 1.27756 1.0<br /> 1 1 1 1 1<br /> 0.08949 1.0<br /> 1 1 1 1 1<br /> 0.29563 1.0<br /> 1 1 1 1 1<br /> 2.20247 1.0<br /></div><div><br /></div><div>I have understood that I should use a primary basis set and a pseudopotential that have been generated considering the same number of valence electrons. In the example above, this is indicated with the -q4 label. </div><div><br /></div><div>However, shouldn't the same rule also apply when choosing the auxiliary basis set?</div><div>The FIT3 is from the file BASIS_ADMM of CP2k version 2023.1, but there is no q label attached to the names of the auxiliary basis sets. </div><div>Given the atom type, can the auxiliary basis sets be used with any pseudopotential?</div><div><br /></div><div>Thanks for your help.</div><div>Best regards,</div><div>Torstein</div><div><br /></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/26c8c114-f50c-4c1f-85e2-4f811c89fc49n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/26c8c114-f50c-4c1f-85e2-4f811c89fc49n%40googlegroups.com</a>.<br />