Table of Contents
Fetching ...

Why are diffuse atomic orbitals needed for accurate electronic wave functions of even neutral molecules?

Samuel R. Powell, Edward F. Valeev

TL;DR

This paper investigates why diffuse atomic orbitals (AOs) are needed for accurate electronic wave functions of neutral molecules, especially when using explicitly correlated F12 methods. It introduces a dual-basis CC-F12 framework that separates occupied and virtual spaces and analyzes how diffuse AOs and F12 terms interact to reduce basis-set incompleteness error (BSIE). The study shows that diffuse AOs remain important for traditional dynamical correlation, but F12 methods dramatically shrink BSIE, making the diffuse AO contribution comparatively larger; dual-basis formulations effectively isolate these effects and offer practical benefits for large systems. The work suggests that dual-basis CC-F12 approaches can enable robust, accurate reduced-scaling calculations in extended systems and solids, aided by system-optimized geminal parameters that separate short-range F12 physics from long-range diffuse correlation.

Abstract

An accurate description of electron correlation energies in molecules requires either basis set extrapolation or the use of explicitly-correlated wave functions that address the deficiencies of standard determinantal expansions at short interelectronic distances. Practical and robust explicitly-correlated F12 methods require the use of standard or specialized atomic orbital (AO) basis sets that include diffuse AOs, even for neutral species. Although modern reduced-scaling formulations of explicitly-correlated many-body methods have become routinely applicable to molecules with hundreds of atoms, application of F12 methods to large molecular systems can be severely hampered due to the onset of ill-conditioning spurred by the presence of diffuse AOs in the F12-appropriate orbital basis sets. Here we re-examine why diffuse AOs are necessary for application of F12 methods. To help such an investigation, we developed a dual-basis formulation of traditional and F12 coupled-cluster singles and doubles (CCSD) methods in which the reference (occupied) and correlating (virtual) orbitals are expanded in separate AO basis sets. Our conclusion is that diffuse AOs are fundamentally important for the traditional (non-F12) description of dynamical correlation; the necessity of diffuse AOs in F12 calculations arises indirectly due to the dramatic reduction of the basis set error by the F12 terms such that the error due to the lack of diffuse AOs becomes comparable to the residual basis set incompleteness. The dual-basis CC methods are suggested as an important candidate formalism for accurate (in particular, F12) reduced-scaling many-body methods in extended systems.

Why are diffuse atomic orbitals needed for accurate electronic wave functions of even neutral molecules?

TL;DR

This paper investigates why diffuse atomic orbitals (AOs) are needed for accurate electronic wave functions of neutral molecules, especially when using explicitly correlated F12 methods. It introduces a dual-basis CC-F12 framework that separates occupied and virtual spaces and analyzes how diffuse AOs and F12 terms interact to reduce basis-set incompleteness error (BSIE). The study shows that diffuse AOs remain important for traditional dynamical correlation, but F12 methods dramatically shrink BSIE, making the diffuse AO contribution comparatively larger; dual-basis formulations effectively isolate these effects and offer practical benefits for large systems. The work suggests that dual-basis CC-F12 approaches can enable robust, accurate reduced-scaling calculations in extended systems and solids, aided by system-optimized geminal parameters that separate short-range F12 physics from long-range diffuse correlation.

Abstract

An accurate description of electron correlation energies in molecules requires either basis set extrapolation or the use of explicitly-correlated wave functions that address the deficiencies of standard determinantal expansions at short interelectronic distances. Practical and robust explicitly-correlated F12 methods require the use of standard or specialized atomic orbital (AO) basis sets that include diffuse AOs, even for neutral species. Although modern reduced-scaling formulations of explicitly-correlated many-body methods have become routinely applicable to molecules with hundreds of atoms, application of F12 methods to large molecular systems can be severely hampered due to the onset of ill-conditioning spurred by the presence of diffuse AOs in the F12-appropriate orbital basis sets. Here we re-examine why diffuse AOs are necessary for application of F12 methods. To help such an investigation, we developed a dual-basis formulation of traditional and F12 coupled-cluster singles and doubles (CCSD) methods in which the reference (occupied) and correlating (virtual) orbitals are expanded in separate AO basis sets. Our conclusion is that diffuse AOs are fundamentally important for the traditional (non-F12) description of dynamical correlation; the necessity of diffuse AOs in F12 calculations arises indirectly due to the dramatic reduction of the basis set error by the F12 terms such that the error due to the lack of diffuse AOs becomes comparable to the residual basis set incompleteness. The dual-basis CC methods are suggested as an important candidate formalism for accurate (in particular, F12) reduced-scaling many-body methods in extended systems.

Paper Structure

This paper contains 13 sections, 29 equations, 4 figures, 3 tables.

Figures (4)

  • Figure 1: The basis set convergence of CCSD and CCSD-F12 energies of the 10-electron benchmark set with and without diffuse AOs.
  • Figure 2: Relative impact of adding diffuse AOs [$R^\text{aug}$ (\ref{['eq:Raug']}), row 1] vs incrementing the OBS cardinal number [$R^{X+1}$ (\ref{['eq:RXp1']}), row 2] on CCSD correlation energy BSIE, without (column 1) and with F12 terms, using standardVRG:peterson:2008:JCP (column 2) and revisedVRG:powell:2025:JCTC (column 3) geminal exponents.
  • Figure 3: Impact of the addition of diffuse AOs on CCSD-F12 correlation energy with three choices of geminal lengthscales (see text) relative to that of the CCSD counterpart.
  • Figure 4: Magnitudes of $\Delta E^\text{DB}$ for CCSD and CCSD-F12 total energies, and the respective differences ($mE_\text{h}$).