Dear CP2K community,<div>I am trying to calculate the binding energy of the argon dimer using RPA. I used the basis set TZV2P-MOLOPT-PBE-GTH-q8 from the file BASIS_MOLOPT_UZH with auxiliary basis set RI_85_7_6_5_5_0_0_0_8.0e-08 from the file BASIS_RI_MOLOPT:</div>BASIS_SET TZV2P-MOLOPT-PBE-GTH-q8<br /><div>BASIS_SET RI_AUX RI_85_7_6_5_5_0_0_0_8.0e-08</div><div><br /></div><div>With this basis I get -17.9 meV as binding energy. Then I tried to remove RI_AUX in front of the auxiliary basis set:<br />BASIS_SET TZV2P-MOLOPT-PBE-GTH-q8<br /><div>BASIS_SET RI_85_7_6_5_5_0_0_0_8.0e-08</div></div><div><br /></div><div>With these parameters i read in the output file:<br />Orbital Basis Set                  RI_85_7_6_5_5_0_0_0_8.0e-08TZV2P-MOLOPT-</div><div>and</div><div>RI Auxiliary Basis Set             Ar-RI-AUX-RI_85_7_6_5_5_0_0_0_8.0e-08TZV</div><div><br /></div><div>It seems that the program mixes the two basis sets, but I don't know how. The result is -10.6 meV which is closer to the experimental value (which is -12.3 meV) than the result I obtained writing RI_AUX in front of the auxiliary basis set. Could anyone explain me what happens with this basis mixing and tell me if the result based in this parameters is reliable? I have attached to this message the input and output files.</div><div>Thank you in advance.</div><div>Kind regards,</div><div><br /></div><div>Riccardo Pezzetta</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/dff85d15-2f99-485d-9556-3ff9bd58c3adn%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/dff85d15-2f99-485d-9556-3ff9bd58c3adn%40googlegroups.com</a>.<br />