Universal perturbative explicitly correlated basis set incompleteness correction

Basis set incompleteness error for an arbitrary approximate electronic wave function is robustly reduced using a second-order perturbative correction into a basis of explicitly correlated, internally contracted geminal functions. The Hylleraas functional for the second-order energy correction is evaluated algebraically involving at most a four-electron reduced density matrix and four-electron integrals. By using the R12 technology in combination with screening approximations such a correction only requires a two-electron reduced density matrix and two-electron integrals. Preliminary investigations of potential energy surfaces of hydrogen fluoride and nitrogen molecules at the multireference configuration interaction singles and doubles indicate that with the perturbative correction only an aug-cc-pVDZ basis is necessary to compute correlation energies of an aug-cc-pVQZ quality, or better. The proposed correction, dubbed [2](R12), can in principle be combined with any single reference and multireference method in use today.