Author Topic: Queerness/Issue observed with RI-basis set  (Read 2595 times)


  • Sr. Member
  • ****
  • Posts: 110
  • Karma: +0/-0
Queerness/Issue observed with RI-basis set
« on: November 20, 2019, 01:26:48 PM »

I encountered the following, at first glance counter-intuitive result for an rimp2 calculation for a deprotonated (anionic) closed-shell molecule.
Since it is anion diffuse functions could be important. So, I set up a mp2 calculation using the aug-cc-pVDZ basis set.
For some reason I changed only the main basis to aug-cc-pVDZ but the auxiliary basis (cbas), which was set to cc-pVDZ. This is case 1.
When I saw the inconsistency, I repeated the calculation with the right, i.e. aug-cc-pVDZ, auxiliary basis set. This is case 2.
When I compared the mp2 energies in case 1 and 2 I found out that the mp2 energy in case 1 (inconsistent cbas basis set)  is slightly
(~0.2 kcal/mol) lower than the energy in case 2 (the main and cbas basis sets are consistent). This was unexpected as I thought that
the larger the cbas basis set, the more closer the result to that obtained without the use of RI. In fact it turned out that the mp2 energy
obtained with the smaller cbas basis set, i.e. case 1, is slightly lower than the mp2 energy obtained without the use of RI.

I wounder whether the above described queerness is known and what the reason of it is (can be)?

Best regards,


  • Full Member
  • ***
  • Posts: 39
  • Karma: +0/-0
Re: Queerness/Issue observed with RI-basis set
« Reply #1 on: September 07, 2023, 07:02:10 PM »

I would think that absolute value of difference between MP2 without RI and RI-MP2 shall converge (?monotonously) with the size of the auxiliary basis set to zero, sure.

But, on the other hand side, is there any theorem such as Ritz variational theorem? MP2 method alone is not variational. Perturbational methods could give energy below FCI and also converge with the size of basis set "from below" to their (atomic) basis set limits.

With auxbasis - could it be, that there is no clear "from below"/"from above" bound and the difference might even be sort of random? Thus perhaps the convergence towards (non-RI-)MP2 even non-monotonous?

Best regards,

P.S.:  I am working with elements which does not have auxbases in the standard library (such as uranium), so I would be interested to know more on this topic too.