Self-consistent density-based basis-set correction: How much do we lower total energies and improve dipole moments?

Accelerated basis-set convergence of coupled-cluster excitation energies using the density-based basis-set correction method

We present the first application to real molecular systems of the recently proposed linear-response theory for the density-based basis-set correction method [J. Chem. Phys., 158, 234107 (2023)]. We apply this approach to accelerate the basis-set convergence of excitation energies in the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) method. We use an approximate linear-response framework that neglects the second-order derivative of the basis-set correction density functional and consists in simply adding to the usual Hamiltonian the one-electron potential generated by the first-order derivative of the functional. This additional basis-set correction potential is evaluated at the Hartree-Fock density, leading to a very computationally cheap basis-set correction. We tested this approach over a set of about 30 excitation energies computed for five small molecular systems and found that the excitation energies from the ground state to Rydberg states are the main source of basis-set error. These excitation energies systematically increase when the size of the basis set is increased, suggesting a biased description in favour of the excited state. Despite the simplicity of the present approach, the results obtained with the basis-set-corrected EOM-CCSD method are encouraging as they yield a mean absolute deviation of 0.02 eV for the aug-cc-pVTZ basis set, while it is 0.04 eV using the standard EOM-CCSD method. This might open the path to an alternative to explicitly correlated approaches to accelerate the basis-set convergence of excitation energies.

A density-fitting implementation of the density-based basis-set correction method

This work reports an efficient density-fitting implementation of the density-based basis-set correction (DBBSC) method in the MOLPRO software. This method consists in correcting the energy calculated by a wave-function method with a given basis set by an adapted basis-set correction density functional incorporating the short-range electron correlation effects missing in the basis set, resulting in an accelerated convergence to the complete-basis-set limit. Different basis-set correction density-functional approximations are explored and the complementary-auxiliary-basis-set single-excitation correction is added. The method is tested on a benchmark set of reaction energies at the second-order Møller-Plesset (MP2) level and a comparison with the explicitly correlated MP2-F12 method is provided. The results show that the DBBSC method greatly accelerates the basis convergence of MP2 reaction energies, without reaching the accuracy of the MP2-F12 method but with a lower computational cost.

Basis Set Limit of CCSD(T) Energies: Explicit Correlation Versus Density-Based Basis-Set Correction

A thorough comparison is carried out for explicitly correlated and density-based basis-set correction approaches, which were primarily developed to mitigate the basis-set incompleteness error of wave function methods. An efficient implementation of the density-based scheme is also presented, utilizing the density-fitting approximation. The performance of these approaches is comprehensively tested for the second-order Møller-Plesset (MP2), coupled-cluster singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] methods with respect to the corresponding complete basis set references. It is demonstrated that the density-based correction together with complementary auxiliary basis set (CABS)-corrected Hartree-Fock energies is highly robust and effectively reduces the error of the standard approaches; however, it does not outperform the corresponding explicitly correlated methods. Nevertheless, what still makes the density-corrected CCSD and CCSD(T) methods competitive is that their computational costs are roughly half of those of the corresponding explicitly correlated variants. Additionally, an incremental approach for standard CCSD and CCSD(T) is introduced. In this simple scheme, the total energies are corrected with the CABS correction and explicitly correlated MP2 contributions. As demonstrated, the resulting methods yield surprisingly good results, below 1 kcal/mol for thermochemical properties even with a double-ζ basis, while their computational expenses are practically identical to those of the density-based basis-set correction approaches.

Electronic excited states in deep variational Monte Carlo

Obtaining accurate ground and low-lying excited states of electronic systems is crucial in a multitude of important applications. One ab initio method for solving the Schrödinger equation that scales favorably for large systems is variational quantum Monte Carlo (QMC). The recently introduced deep QMC approach uses ansatzes represented by deep neural networks and generates nearly exact ground-state solutions for molecules containing up to a few dozen electrons, with the potential to scale to much larger systems where other highly accurate methods are not feasible. In this paper, we extend one such ansatz (PauliNet) to compute electronic excited states. We demonstrate our method on various small atoms and molecules and consistently achieve high accuracy for low-lying states. To highlight the method's potential, we compute the first excited state of the much larger benzene molecule, as well as the conical intersection of ethylene, with PauliNet matching results of more expensive high-level methods.

Ground- and Excited-State Dipole Moments and Oscillator Strengths of Full Configuration Interaction Quality

We report ground- and excited-state dipole moments and oscillator strengths (computed in different "gauges" or representations) of full configuration interaction (FCI) quality using the selected configuration interaction method known as Configuration Interaction using a Perturbative Selection made Iteratively (CIPSI). Thanks to a set encompassing 35 ground- and excited-state properties computed in 11 small molecules, the present near-FCI estimates allow us to assess the accuracy of high-order coupled-cluster (CC) calculations including up to quadruple excitations. In particular, we show that incrementing the excitation degree of the CC expansion (from CC with singles and doubles (CCSD) to CC with singles, doubles, and triples (CCSDT) or from CCSDT to CC with singles, doubles, triples, and quadruples (CCSDTQ)) reduces the average error with respect to the near-FCI reference values by approximately 1 order of magnitude.

Basis-set correction for coupled-cluster estimation of dipole moments

The present work proposes an approach to obtain a basis-set correction based on density-functional theory (DFT) for the computation of molecular properties in wave-function theory (WFT). This approach allows one to accelerate the basis-set convergence of any energy derivative of a non-variational WFT method, generalizing previous works on the DFT-based basis-set correction where either only ground-state energies could be computed with non-variational wave functions [Loos et al., J. Phys. Chem. Lett. 10, 2931 (2019)] or properties could be computed as expectation values over variational wave functions [Giner et al., J. Chem. Phys. 155, 044109 (2021)]. This work focuses on the basis-set correction of dipole moments in coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)], which is numerically tested on a set of 14 molecules with dipole moments covering two orders of magnitude. As the basis-set correction relies only on Hartree–Fock densities, its computational cost is marginal with respect to the one of the CCSD(T) calculations. Statistical analysis of the numerical results shows a clear improvement of the basis convergence of the dipole moment with respect to the usual CCSD(T) calculations.

Accurate energies of transition metal atoms, ions, and monoxides using selected configuration interaction and density-based basis-set corrections

The semistochastic heat-bath configuration interaction method is a selected configuration interaction plus perturbation theory method that has provided near-full configuration interaction (FCI) levels of accuracy for many systems with both single- and multi-reference character. However, obtaining accurate energies in the complete basis-set limit is hindered by the slow convergence of the FCI energy with respect to basis size. Here, we show that the recently developed basis-set correction method based on range-separated density functional theory can be used to significantly speed up basis-set convergence in SHCI calculations. In particular, we study two such schemes that differ in the functional used and apply them to transition metal atoms and monoxides to obtain total, ionization, and dissociation energies well converged to the complete-basis-set limit within chemical accuracy.