A simple F12 geminal correction in multi-reference perturbation theory

Efficient Exploitation of Numerical Quadrature with Distance-Dependent Integral Screening in Explicitly Correlated F12 Theory: Linear Scaling Evaluation of the Most Expensive RI-MP2-F12 Term

We present a linear scaling atomic orbital based algorithm for the computation of the most expensive exchange-type RI-MP2-F12 term by employing numerical quadrature in combination with CABS-RI to avoid six-center-three-electron integrals. Furthermore, a robust distance-dependent integral screening scheme, based on integral partition bounds [Thompson, T. H.; Ochsenfeld, C. J. Chem. Phys. 2019, 150, 044101], is used to drastically reduce the number of the required three-center-one-electron integrals substantially. The accuracy of our numerical quadrature/CABS-RI approach and the corresponding integral screening is thoroughly assessed for interaction and isomerization energies across a variety of numerical integration grids. Our method outperforms the standard density fitting/CABS-RI approach with errors below 1 μEh even for small grid sizes and moderate screening thresholds. The choice of the grid size and screening threshold allows us to tailor our ansatz to a desired accuracy and computational efficiency. We showcase the approach's effectiveness for the chemically relevant system valinomycin, employing a triple-ζ F12 basis set combination (C54H90N6O18, 5757 AO basis functions, 10,266 CABS basis functions, 735,783 grid points). In this context, our ansatz achieves higher accuracy combined with a 135× speedup compared to the classical density fitting based variant, requiring notably less computation time than the corresponding RI-MP2 calculation. Additionally, we demonstrate near-linear scaling through calculations on linear alkanes. We achieved an 817-fold acceleration for C80H162 and an extrapolated 28,765-fold acceleration for C200H402, resulting in a substantially reduced computational time for the latter─from 229 days to just 11.5 min. Our ansatz may also be adapted to the remaining MP2-F12 terms, which will be the subject of future work.

Nonunitary projective transcorrelation theory inspired by the F12 ansatz

An alternative nonunitary transcorrelation, inspired by the F12 ansatz, is investigated. In contrast to the Jastrow transcorrelation of Boys-Handy, the effective Hamiltonian of this projective transcorrelation features: 1. a series terminating formally at four-body interactions. 2. no spin-contamination within the non-relativistic framework. 3. simultaneous satisfaction of the singlet and triplet first-order cusp conditions. 4. arbitrary choices of pairs for correlation including frozen-core approximations. We discuss the connection between the projective transcorrelation and F12 theory with applications to small molecules, to show that the cusp conditions play an important role to reduce the uncertainty arising from the nonunitary transformation.

Highly Efficient and Accurate Computation of Multiple Orbital Spaces Spanning Fock Matrix Elements on Central and Graphics Processing Units for Application in F12 Theory

We employ our recently published highly efficient seminumerical exchange (sn-LinK) [Laqua, H.; Thompson, T. H.; Kussmann, J.; Ochsenfeld, C. J. Chem. Theory Comput. 2020, 16, 1456-1468] and integral-direct resolution of the identity Coulomb (RI-J) [Kussmann, J.; Laqua, H.; Ochsenfeld, C. J. Chem. Theory Comput. 2021, 17, 1512-1521] methods to significantly accelerate the computation of the demanding multiple orbital spaces spanning Fock matrix elements present in R12/F12 theory on central and graphics processing units. The errors introduced by RI-J and sn-LinK into the RI-MP2-F12 energy are thoroughly assessed for a variety of basis sets and integration grids. We find that these numerical errors are always below "chemical accuracy" (∼1 mH) even for the coarsest settings and can easily be reduced below 1 μH by employing only moderately large integration grids and RI-J basis sets. Since the number of basis functions of the multiple orbital spaces is notably larger compared with conventional Hartree-Fock theory, the efficiency gains from the superior basis scaling of RI-J and sn-LinK (O(N_bas²) instead of O(N_bas⁴) for both) are even more significant, with maximum speedup factors of 37 000 for RI-J and 4500 for sn-LinK. In total, the multiple orbital spaces spanning Fock matrix evaluation of the largest tested structure using a triple-ζ F12 basis set (5058 AO basis functions, 9267 CABS basis functions) is accelerated over 1575× using CPUs and over 4155× employing GPUs.

Improving the accuracy of the variational quantum eigensolver for molecular systems by the explicitly-correlated perturbative [2]R12correction

We provide an integration of the universal, perturbative explicitly correlated [2]_R12-correction in the context of the Variational Quantum Eigensolver (VQE). This approach is able to increase the accuracy of the underlying reference method significantly while requiring no additional quantum resources. The proposed approach only requires knowledge of the one- and two-particle reduced density matrices (RDMs) of the reference wavefunction; these can be measured after having reached convergence in the VQE. This computation comes at a cost that scales as the sixth power of the number of electrons. We explore the performance of the VQE + [2]_R12 approach using both conventional Gaussian basis sets and our recently proposed directly determined pair-natural orbitals obtained by multiresolution analysis (MRA-PNOs). Both Gaussian orbital and PNOs are investigated as a potential set of complementary basis functions in the computation of [2]_R12. In particular the combination of MRA-PNOs with [2]_R12 has turned out to be very promising - persistently throughout our data, this allowed very accurate simulations at a quantum cost of a minimal basis set. Additionally, we found that the deployment of PNOs as complementary basis can greatly reduce the number of complementary basis functions that enter the computation of the correction at a complexity.

Single-Root Multireference Brillouin-Wigner Perturbative Approach to Excitation Energies

The state-specific Brillouin-Wigner multireference perturbation theory [which employs Jeziorski-Monkhorst parametrization of the wave function] using improved virtual orbitals, denoted as IVO-BWMRPT, is applied to calculate excitation energies (EEs) for methylene, ethylene, trimethylenemethane, and benzyne systems exhibiting various degrees of diradical character. In IVO-BWMRPT, all of the parameters appearing in the wave function ansatz are optimized for a specific electronic state. For these systems, the IVO-BWMRPT method provides EEs that are in close agreement with the benchmark results and experiments, where available, indicating that the method does not introduce imbalance in the target-specific treatment of closed- and open-shell states involved. The good performance of the present methodology is primarily related to structural compactness of the formalism. Overall, present findings are encouraging for both further development of the approach and chemical applications on the energy differences of strongly correlated systems.

Scalable Electron Correlation Methods. 8. Explicitly Correlated Open-Shell Coupled-Cluster with Pair Natural Orbitals PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12

We present explicitly correlated open-shell pair natural orbital local coupled-cluster methods, PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. The methods are extensions of our previously reported PNO-R/UCCSD methods (J. Chem. Theory Comput., 2020, 16, 3135-3151, https://pubs.acs.org/doi/10.1021/acs.jctc.0c00192) with additions of explicit correlation and perturbative triples corrections. The explicit correlation treatment follows the spin-orbital CCSD-F12b theory using Ansatz 3*A, which is found to yield comparable or better basis set convergence than the more rigorous Ansatz 3C in computed ionization potentials and reaction energies using double- to quaduple-ζ basis sets. The perturbative triples correction is adapted from the spin-orbital (T) theory to use triples natural orbitals (TNOs). To address the coupling due to off-diagonal Fock matrix elements, the local triples amplitudes are iteratively solved using small domains of TNOs, and a semicanonical (T0) domain correction with larger domains is applied to reduce the domain errors. The performance of the methods is demonstrated through benchmark calculations on ionization potentials, radical stabilization energies, reaction energies of fragmentations and rearrangements in radical cations, and spin-state energy differences of iron complexes. For a few test sets where canonical calculations are feasible, PNO-RCCSD(T)-F12 results agree with the canonical ones to within 0.4 kcal mol^-1, and this maximum error is reduced to below 0.2 kcal mol^-1 when large local domains are used. For larger systems, results using different thresholds for the local approximations are compared to demonstrate that 1 kcal mol^-1 level of accuracy can be achieved using our default settings. For a couple of difficult cases, it is demonstrated that the errors from individual approximations are only a fraction of 1 kcal mol^-1, and the overall accuracy of the method does not rely on error compensations. In contrast to canonical calculations, the use of spin-orbitals does not lead to a significant increase of computational time and memory usage in the most expensive steps of PNO-R/UCCSD(T)-F12 calculations. The only exception is the iterative solution of the (T) amplitudes, which can be avoided without significant errors by using a perturbative treatment of the off-diagonal coupling, known as (T1) approximation. For most systems, even the semicanonical approximation (T0) leads only to small errors in relative energies. Our program is well parallelized and capable of computing accurate correlation energies for molecules with 100-200 atoms using augmented triple-ζ basis sets in less than a day of elapsed time on a small computer cluster.

Optimal composition of atomic orbital basis sets for recovering static correlation energies

Static correlation energies (Estat) are calculated in a range of basis sets for a chemically diverse collection of atoms and molecules. The reliability of a basis set in capturing Estat is assessed according to the following: mean and maximum absolute deviations from near-exact Estat estimates, monotonic convergence to the complete basis set limit, and ability to capture Estat accurately independent of changes in geometry, molecular size, and electronic configuration. Within the polarization and correlation-consistent basis set series, triple-ζ basis sets are the smallest that can reliably capture Estat. The cc-pVTZ basis set performs particularly well, recovering Estat to chemical accuracy for all atoms and molecules in our data set. A series of customized basis sets are constructed by stripping polarization functions from, and swapping polarization functions among, existing basis sets. Basis sets without polarization functions are incapable of accurately recovering Estat. Basis sets with a near-complete set of s, p, and d functions can approach chemical accuracy in maximum absolute error. However, this may be achieved at lower computational cost by using a well balanced triple-ζ basis set including f functions, along with a smaller number of s, p, and d functions. Recommended basis sets for calculating Estat with increasing accuracy at increasing computational cost are 6-311G(2d,2p), cc-pVTZ, and cc-pVQZ stripped of g functions.

Explicitly correlated second-order Møller-Plesset perturbation theory employing pseudospectral numerical quadratures

We implemented explicitly correlated second-order Møller-Plesset perturbation theory with numerical quadratures using pseudospectral construction of grids. Introduction of pseudospectral approach for the calculation of many-electron integrals gives a possibility to use coarse grids without significant loss of precision in correlation energies, while the number of points in the grid is reduced about nine times. The use of complementary auxiliary basis sets as the sets of dealiasing functions is justified at both theoretical and computational levels. Benchmark calculations for a set of 16 molecules have shown the possibility to keep an error of second-order correlation energies within 1 milihartree (mH) with respect to MP2-F12 method with dense grids. Numerical tests for a set of 13 isogyric reactions are also performed.