Density functional theory with plane-wave basis sets is widely employed in computational materials science, including applications to isolated molecular systems. However, the inadequate description of electron correlation remains a fundamental limitation. Accurate correlation treatments based on many-body Hamiltonians require reliable representations of both occupied and virtual orbitals, yet virtual orbitals are often poorly described in conventional computational schemes, resulting in reduced accuracy. In this work, we introduce localized correlation-converged virtual orbitals (LCCVOs) as an efficient basis for constructing accurate many-body Hamiltonians in molecular systems. Using a substantially reduced number of orbitals, the LCCVO framework yields dissociation energies for singlet, doublet, and triplet molecules that are comparable to, and in many cases exceed, those obtained with high-level correlation-consistent basis sets such as cc-pVXZ (X = D, T, Q, 5). These results demonstrate the efficiency, scalability, and robustness of the LCCVO approach for high-accuracy quantum chemical calculations.
Density functional theory with plane-wave basis sets is widely employed in computational materials science, including applications to isolated molecular systems. However, the inadequate description of electron correlation remains a fundamental limitation. Accurate correlation treatments based on many-body Hamiltonians require reliable representations of both occupied and virtual orbitals, yet virtual orbitals are often poorly described in conventional computational schemes, resulting in reduced accuracy. In this work, we introduce localized correlation-converged virtual orbitals (LCCVOs) as an efficient basis for constructing accurate many-body Hamiltonians in molecular systems. Using a substantially reduced number of orbitals, the LCCVO framework yields dissociation energies for singlet, doublet, and triplet molecules that are comparable to, and in many cases exceed, those obtained with high-level correlation-consistent basis sets such as cc-pVXZ (X = D, T, Q, 5). These results demonstrate the efficiency, scalability, and robustness of the LCCVO approach for high-accuracy quantum chemical calculations.
Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron-electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. The transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.
For many-body methods such as MCSCF and CASSCF, in which the number of one-electron orbitals is optimized and independent of the basis set used, there are no problems with using plane-wave basis sets. However, for methods currently used in quantum computing such as select configuration interaction (CI) and coupled cluster (CC) methods, it is necessary to have a virtual space that is able to capture a significant amount of electron-electron correlation in the system. The virtual orbitals in a pseudopotential plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. The use of virtual spaces derived from the one-electron operators has also been tried, and while some correlations are captured, the amount is quite low. To overcome these limitations, we have been developing new classes of algorithms to define virtual spaces by optimizing orbitals from small pairwise CI Hamiltonians, which we term as correlation optimized virtual orbitals with the abbreviation COVOs. With these procedures, we have been able to derive virtual spaces, containing only a few orbitals, which are able to capture a significant amount of correlation. The focus in this manuscript is on using these derived basis sets to target full CI (FCI) quality results for H2 on near-term quantum computers. However, the initial results for this approach were promising. We were able to obtain good agreement with FCI/cc-pVTZ results for this system with just 4 virtual orbitals, using both FCI and quantum simulations. The quality of the results using COVOs suggests that it may be possible to use them in other many-body approaches, including coupled cluster and Møller–Plesset perturbation theories, and open up the door to many-body calculations for pseudopotential plane-wave basis set methods.
We introduce a novel and efficient explicitly correlated implementation of second-order perturbation theory for solids. The required three-electron integrals are computed directly using a plane wave basis set. We parametrize the employed correlation factors using results previously obtained for a finite uniform electron gas simulation cell. We demonstrate for a range of solids that basis set converged correlation energies, equilibrium volumes, and bulk moduli can be obtained efficiently in this theory using a few ten orbitals per atom. To stretch the capabilities of this novel method we compute the Schottky defect formation energy in MgO, studying systems with 54 atoms in the supercell. We verify the accuracy of the calculated formation energies using the more accurate coupled cluster singles and doubles theory. Furthermore, we discuss other potential applications for the derived and implemented expressions such as an occupied orbital only correlation energy functional.
A well-parallelized local singles and doubles coupled-cluster (LCCSD) method using pair natural virtual orbitals (PNOs) is presented. The PNOs are constructed using large domains of projected atomic orbitals (PAOs) and orbital specific virtual orbitals (OSVs). We introduce a hierarchy of close, weak, and distant pairs, based on pair energies evaluated by local Møller–Plesset perturbation theory (LMP2). In contrast to most previous implementations of LCCSD methods, the close and weak pairs are not approximated by LMP2 but treated by higher-order methods. This leads to greatly improved accuracy for large systems, in particular when long-range dispersion interactions are important. Close pair amplitudes are optimized using approximate LCCSD equations, in which slowly decaying terms that mutually cancel at long-range are neglected. For weak pairs, the same approximations are used, but in addition, the nonlinear terms are neglected (coupled electron pair approximation). Distant pairs are treated by spin-component scaled (SCS)-LMP2 using multipole approximations. For efficiency reasons, various projection approximations are also introduced. The impact of these approximations on absolute and relative energies depends on the PNO domain sizes. The errors are found to be negligible, provided that sufficiently large PNO domains are used for close and weak pairs. For the selection of these domains the usual natural orbital occupation number criterion is found to be insufficient, and an additional energy criterion is used. For extended one-dimensional systems, the computational effort of the method scales nearly linearly with the number of correlated electrons, but the linear scaling regime is usually not reached in real-life applications for three-dimensional systems. Nevertheless, due to the parallelization that is efficient up to about 100–200 CPU cores (dependent on the molecular size), accurate calculations for three-dimensional molecules with about 100 atoms and augmented triple-ζ basis sets (e.g., cc-pVTZ-F12) can be carried out in 1–3 h of elapsed time (depending on the molecular structure and the number of CPU cores, excluding the time for Hartree–Fock). The convergence of the results with respect to the thresholds and options that control the domain, pair, and projection approximations is extensively tested. Benchmark examples for several difficult and large cases are presented, which demonstrate that the errors of relative energies (e.g., reaction energies, barrier heights) caused by the pair and projection approximations can be reduced to below 1 kJ mol –1 . The remaining errors are mainly caused by the finite PAO domains. The larger these are made, the more intramolecular or intermolecular basis set superposition errors are introduced, and the canonical results are approached only very slowly. This problem is eliminated in the explicitly correlated variant of our method (PNO-LCCSD-F12), which will be described in a separate paper, along with a larger set of benchmark calculations.
The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.
We propose a general method for constructing system-dependent basis functions for correlated quantum calculations. Our construction combines features from several traditional approaches: plane waves, localized basis functions, and wavelets. In a one-dimensional mimic of Coulomb systems, it requires only 2--3 basis functions per electron to achieve high accuracy, and reproduces the natural orbitals. We illustrate its effectiveness for molecular energy curves and chains of many one-dimensional atoms. We discuss the promise and challenges for realistic quantum chemical calculations.