Abstract
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.
Introduction
Quantum chemistry is one of the first and most successful scientific applications of digital computers (; ; ; ; ; Roothaan, 1951; ; ; ; ; Pople et al., 1965; ; Reeves, 1966; Pulay, 1969). This success has led to a large number of research, open-source (; Stanton et al., 1992; Schmidt et al., 1993; ; ; ; ; ; ; ; ; ), and commercial codes (; te Velde et al., 2001; ; ; Werner et al., 2012; Shao et al., 2015; ; ) (for a larger list of quantum chemistry software, see (Wikipedia, The Free Encyclopedia, 2020)), which are used on a regular basis by tens of thousands of scientists, engineers, and students from a variety of scientific and engineering domains. With Moore’s law as a backdrop (), the cycle of new machines leading to new algorithms stimulated the field for many decades, and as a consequence, a large number of quantum chemistry methods were developed along with a variety of numerical methods to solve them. However, in recent decades, the maturity and success of these codes coupled with the imminent death of Moore’s law (; Rotman, 2020) that made numerical software development much more difficult and less accessible to the average scientist have resulted in the field having priorities other than just new science, such as porting and optimizing these codes to the next generation of computers (, ; Richard et al., 2018; van Dam et al., 2020), standardization of methods (; Wilkins-Diehr and Crawford, 2018), and marketing (; ).
With the advent of quantum computing, there is excitement again, and quantum chemists are beginning to rethink how they carry out quantum chemistry calculations, in particular very accurate and very expensive instances of systems containing strong electron-electron correlations. This is because it is anticipated that quantum computers with 50–100 qubits will be able to surpass classical digital computers for these types of calculations (Preskill, 2018). Quantum computing has thus emerged as an alternative avenue to the continuity of quantum chemistry in the long run (Wasielewski et al., 2020) but poses several challenges that demand careful consideration in order to eventually mature into a viable replacement for classical computers and large, highly parallelizable high-performance computing clusters.
Present quantum devices are plagued by short coherence times and vulnerability to environment interference, i.e., noise. Albeit quantum algorithms have been developed with proved exactness, such as quantum phase estimation, these are not a viable option in the present/near-term time frame. Therefore, it is desirable to limit the operation of quantum processors to a complementary concerted execution with classical counterparts, whereby each of these components is only in charge of those tasks for which it is more suitable. This has materialized into the variational quantum eigensolver (VQE) (Peruzzo et al., 2014) and other hybrid algorithms. Briefly, this class of algorithms strives to find the lowest eigenvalue of a given observable by assuming that the associated quantum state can be accurately represented by a trial wave function and whose parameters are varied according to the Rayleigh—Ritz method (variational principle), with these parameters being updated by the classical computer. The burden on the quantum processor can be further alleviated with strategies such as Trotterization, which in turn introduce other challenges (; ) but can be successfully exploited in the construction of favorable ansatz, as long as a predefined form for the trial wave function is imposed. This is at the heart of the ADAPT-VQE ().
Most high-level methods for strongly correlated systems in use today (e.g., full configuration interaction (CI), coupled cluster (CC) and Green’s function (GF) approaches) are based on second-quantized Hamiltonians, which are written in terms of creation and annihilation operators for fermion orbitals. These methods are amenable to quantum computers because fermionic creation and annihilation operators can be readily mapped to qubits through the use of some established transformation, among which Jordan—Wigner (), Bravyi—Kitaev (), and binary codes (Steudtner and Wehner, 2018) stand out, where the number of qubits scales with the number of orbitals in the second-quantized Hamiltonian. In principle, converting the full many-body electronic Hamiltonian to a second-quantized form is exact and popular CC and GF approximations based on this form are very accurate. However, this conversion has a drawback in that it requires the introduction of a basis set, which, for computational cost reasons, needs to be small. Typically, these basis sets are composed of atomic-like orbitals generated with heuristics based on an atom calculation for each kind of atom in the system. An example of this type of basis set is the popular Dunning correlation consistent basis set (; ) in which the atomic orbitals are optimized at the CISD (configuration interaction method with single and double excitations) level of theory. While the size of this basis set is small compared to other basis sets used in quantum chemistry, such as plane waves, it still needs to contain a large number of atomic orbitals to produce a truly accurate result.
Solving relevant chemistry problems analogously to what is classically done with MCSCF or FCI on near-term quantum computers that contain 10s to 100s of noisy qubits (Reiher et al., 2017), in which only limited numbers of operations can be performed, is a monumental challenge. One way to reduce the cost of these calculations is to develop new procedures for optimizing basis sets. In this manuscript, a new method is presented for generating a plane-wave derived correlation optimized orbital basis sets. These derived basis sets can also be used in other many-body approaches, including CC theory, and can easily be generalized to work with recently developed Filon’s Integration Strategy for two-electron integrals in periodic systems (). This method is different from other plane-wave derived optimized orbital basis sets (Shirley, 1996; Prendergast and Louie, 2009; ); in that, it is based on optimizing small select CI problems rather than fitting one-electron eigenvalue spectra and band structures.
The paper is organized as follows. In Section 2, a brief description of the second-quantized Hamiltonian and the double unitary CC downfolding method that can be used with the pseudopotential plane-wave method is given, followed by comparisons between restricted Hartree–Fock (RHF) calculations using plane-wave and Gaussian basis sets. Using this framework, CI calculations up to 20 virtual orbitals, generated from plane-wave Hartree–Fock and one-electron Hamiltonians, are shown for the H2 molecule. The VQE quantum computing algorithms used in this work are described in Section 3. Section 4 presents a new class of algorithms for generating a virtual space in which the orbitals are generated by minimizing small pairwise CI Hamiltonians, and a complete set of equations for implementing these optimizations is given in Subsections 4.1–4.4. Using this new type of virtual space, CI calculations up to 18 virtual orbitals for the ground state energy curve of the H2 molecule are presented in Section 5 followed by results using quantum computing simulations in Section 6, and lastly, the conclusions are given in Section 7.
Pseudopotential Plane-Wave Many-Body Hamiltonian
The nonrelativistic electronic Schrödinger eigenvalue equation of quantum chemistry can be written aswhere H is the electronic structure Hamiltonian under the Born–Oppenheimer approximation and is the quantum mechanical wave function that is a function of the spatial and spin coordinates of the electrons, . When solving this equation, the Pauli exclusion principle constraint of particle exchange must be enforced, in which the wave function changes sign when the coordinates of two particles, and , are interchanged; i.e.,
For the Born–Oppenheimer Hamiltonian, the interaction between the electrons and nuclei is described by the proper potentials , which for plane-wave solvers can cause trouble with convergence because of the singular behavior at . A standard way to remove this issue in plane-wave calculations is to replace these singular potentials with pseudopotentials. By making this replacement, the Hamiltonian, H, in Eq. 1 can be written aswhere the first term is the kinetic energy operator, the second term contains the local and nonlocal pseudopotentials, and represent the electron-ion interactions, and the last term is the electron-electron repulsion.
Instead of writing the many-electronic Hamiltonian in the traditional Schrödinger form, as in the equations above, it is more common today to write it in an alternative representation, known as the second-quantization form, defined using the creation, , and annihilation, , operators. The second-quantized Hamiltonian is written aswhere represent the one-electron spin-orbital basis. A nice feature about this form of the Hamiltonian is that the antisymmetry of wavefunction requirement as given in Eq. 2 is automatically enforced through the standard fermionic anticommutation relations and .
In this formulation, the choice of the one-electron spin-orbital basis is nebulous and requires some care in its choosing in order to obtain accurate results with this type of Hamiltonian. Typically, in quantum chemistry, one uses the filled and virtual orbitals from a Hartree–Fock calculation. For methods that utilize linear combinations of atomic orbitals (LCAO) as the basis, the size of the basis set and subsequently generated Hartree–Fock orbitals is fairly small. However, for plane-wave solvers and other grid-based solvers, the size of the basis set is very large and the number of the one- and two-electron integrals in Eq. 4 will become prohibitive if all possible Hartree–Fock orbitals are used.
One approach to this problem is to only include virtual orbital up to a certain energy threshold, and another related approach is to use the plane-wave derived optimized orbital basis set, e.g., the Shirley approach. While the number of these orbitals needed to accurately describe eigenvalue spectra over a range of ∼100 eV is significantly smaller than the number of plane waves, it is still significantly larger than the number of orbitals generated by an LCAO method. The reason for this is that the virtual orbitals in a plane-wave Hartree–Fock calculation, because of Coulomb repulsion, are often unbound scattering states that interact very weakly with the filled orbitals. As a result, very little correlation energy is captured from them. In contrast, LCAO basis methods can only describe bound states, and hence, Hartree–Fock calculations on this basis do not generate these types of scattering states.
Many-Body Downfolding Techniques
One technique for reducing the dimensionality of large plane-wave calculations is to construct effective Hamiltonians that capture correlation effects of the full calculation in manageable active spaces. That way, all many-body effects are retained as opposed to simply truncating the orbital space. In ), we introduced downfolding techniques, which utilize the double unitary CC (DUCC) ansatz for exact ground state wave function ,where and are the general type anti-Hermitian operatorsdefined by amplitudes defining action within and outside of the predefined active space, respectively; i.e., the amplitudes defining the operator must carry at least one inactive spin-orbital index whereas all amplitudes defining the operator carry active spin-orbital indices only. In Eq. 5, designates properly chosen reference function (usually chosen as a Hartree–Fock (HF) Slater determinant). The exactness of the expansion 5 has been recently discussed in ), where it was also shown that the standard UCC expansions can provide a basic approximation of the exact and operators, i.e.,where and are single-reference-type internal and external cluster amplitudes (in the sense defined above).
Using ansatz in Eq. 5 we have shown () that the exact energy of the systems can be reproduced by the diagonalization of the effective (or downfolded) Hamiltonian in the corresponding active space:where
In Eq. 11, P and are the projection operators onto the reference function and all active-space excited Slater determinants (with respect to ).
We will discuss the utility of the downfolding techniques in the next section for the ground state calculations of H2. This is just one of the two approaches presented in this paper for reducing the dimensionality of the quantum problem (Figure 1).
FIGURE 1
Results for the 1 Ground State of H2 Using Virtual Space From Hartree–Fock and One-Electron Hamiltonians
The NWChem program package (; Valiev et al., 2010; ; ; ) was used for all calculations in this study, except for the FCI calculations, which used the TINYMRCC suite by Jiří Pittner. The plane-wave calculations used a simple cubic box with L = 26a0 and cutoff energy of 100 Ry. The valence electron interactions with the atomic core are approximated with generalized norm-conserving Hamann () pseudopotentials modified to the separable form suggested by Kleinman and Bylander (). The pseudopotentials used in this study were constructed using the following core radii: H: rcs = 0.8 a.u. and rcp = 0.8 a.u.; Be: rcs = 1.219 a.u., and rcp = 1.219 a.u. (vide infra). The RHF and coupled cluster singles and doubles (CCSD) LCAO calculations used the Dunning cc-pVTZ Gaussian basis set.
As can be seen in Figure 2, the RHF ground state energy curve of the H2 molecule using plane-wave and LCAO Gaussian basis sets gives nearly identical results. However, when we performed plane-wave FCI calculations (not shown) for this system using up to 20 RHF virtual orbitals, the amount of correlation energy calculated was nearly zero (<1.0 e-4 Hartree). This result was not surprising since most of the virtual states were scattering states as shown in Figure 3.
FIGURE 2
FIGURE 3
Instead of using virtual states of the RHF Hamiltonian, virtual states were also generated using the 1-electron part of the RHF Hamiltonian, (i.e., just the kinetic energy and pseudopotential terms). As shown in Figure 3, the Hamiltonian generated virtual orbitals that were bound and looked like the virtual orbitals generated in an LCAO RHF calculation. Using these generated orbitals, we performed plane-wave CI calculations using 19 of these virtual orbitals. As seen in Figure 2 and Table 1, a significant improvement was seen using these orbitals as they were able to capture a nonzero amount of the correlation energy; however, it was still significantly less than that found in LCAO calculations. In addition, results using the quantum phase estimation (QPE) algorithm in the Microsoft QDK package (Svore et al., 2018; ) in which the number of orbitals was reduced to 4 and 6 orbitals using the DUCC method are shown. These results showed that the DUCC QDK QPE method produces total energies that are within a few milli-Hartrees of the 20 orbital FCI result with only 4 or 6 orbitals.
TABLE 1
| R (H-H) | PW FCI | PW QDK | PW QDK | CCSD |
|---|---|---|---|---|
| (Å) | 19 Virt. | DUCC 4 | DUCC 6 | cc-pVTZ |
| 0.423 | −0.99396 | −0.99113 | −0.99052 | −1.01540 |
| 0.529 | −1.10715 | −1.10363 | −1.10440 | −1.12144 |
| 0.741 | −1.15340 | −1.14968 | −1.15010 | −1.17234 |
| 1.058 | −1.11042 | −1.10768 | −1.10785 | −1.13617 |
| 1.588 | −1.01251 | −1.01435 | −1.01417 | −1.05526 |
| 2.117 | −0.94070 | −0.94303 | −0.94318 | −1.01485 |
| 4.233 | −0.83962 | −0.84025 | −0.84156 | −0.99965 |
Total energies as a function of distance from plane-wave FCI calculations for the H2 molecule 19 virtual orbitals.
Variational Quantum Eigensolver Methods
VQE is a method to find the quantum state that minimizes a cost function defined in operator form (Peruzzo et al., 2014; ). This is a hybrid computational approach in which the preparation of the quantum circuit is tuned using feedback from classical evaluations of the cost function. Reduction of a given problem to minimization, such as solving for the ground state energy (lowest energy eigenvalue) of a molecular Hamiltonian, may then rely on the variational principle to affirm that only the true ground state could satisfy the minimum energy (; ).
Formally, we may consider the problem of solving for the ground state energy, , aswhere H represents the second-quantized Hamiltonian of Eq. 4 and labels the electronic configuration. Within the context of quantum computing, the fermionic representations of the Hamiltonian and state are transformed to alternate representations in terms of the spin operators. This transformation recast the molecular Hamiltonian into a representation that is defined in terms of the usual Pauli spin operators. Similarly, the electronic state is represented as a variable unitary operator acting on a fiducial quantum state (|0>). This yields the equivalent representation of Eq. 13 as
The second equality is pertinent to the current context as it makes explicit the fact that 1) , the Hamiltonian in terms of strings of Pauli operators, relates to H through some transformation that maps fermionic creation and annihilation operators to qubits operators, and 2) the trial wave function emerges from the action of the parameterized unitary operator which builds entanglement, usually starting from an unentangled wave function, such as Hartree–Fock. For practical considerations, H is transformed into with the Jordan-Wigner transformation (), but alternatives have been reported in the literature (; Setia and Whitfield, 2018; Steudtner and Wehner, 2018). The classical search for the quantum state that minimizes the energy represents the conventional computing task, while evaluation of the expectation value is performed using the quantum computer. In particular, the quantum state is prepared by executing a quantum circuit, which is expressed formally as a series of unitary operators acting on a well-defined initial state. The action of a specific sequence of unitaries is to prepare a given state that is subject to the measurements necessary to recover the desired expectation value.
In practice, the quantum state that minimizes the energy is unknown, and, therefore, a search over possible unitaries is necessary to find the form that minimizes the energy. This variational approach to circuit synthesis underlies the VQE method and an essential choice is the selection of a quantum circuit ansatz which defines the range of unitaries that may be formed to prepare and evaluate a quantum state. For electronic structure calculations, seemingly randomized unitaries may offer advantages for efficient circuit construction, but they lack much of the intuition available from theoretical chemistry (). Rather, ansatz circuits derived from unitary coupled cluster theory offer a convenient connection to the expected unitary forms of the minimal quantum state (Romero et al., 2017).
VQE has been applied previously to recover the electronic energy from the Hamiltonian presented in Eq. 4. The literature provides several examples of usage of VQE for problems of chemical interest, in terms of both simulation and implementation on actual quantum hardware. Given the current limitations faced by present quantum computers, these instances are usually accompanied by strategies that reduce the effective Hilbert space, thus leading to a decrease in the computational expense, such as the use of active spaces and natural orbitals (Verma et al., 2020), as well as downfolding techniques introduced earlier. Another route is to modify the form of the ansatz; an example of this alternative would be the so-called Trotterization, which can be used in conjunction with Hilbert space-reducing techniques.
Recently, the principle of VQE was extended to use ansatz circuits that are tailored to computational chemistry applications and specifically the unitary coupled cluster singles and doubles (UCCSD) ansatz state. Adaptive ansatz construction is attractive because it obeys the underlying complexity of the electronic structure in question, whereas a predefined form for the trial wave function in Eq. 13 may fall short of the flexibility necessary for intricate problems. The prime example of this class of algorithms is the ADAPT-VQE, which iteratively assembles a circuit according to the expected energy gain signaled by the gradient with respect to the variational parameters.
An important consideration in the performance of both VQE and ADAPT-VQE is the depth of the ansatz circuit and the time required to construct the optimal variational circuit. For electronic structures dominated by weak correlation, ADAPT-VQE tends to be very economical, adding only operators that make a meaningful contribution toward the lowest eigenvalue in the spectrum of the Hamiltonian in Eq. 14. On the other hand, the usual UCCSD, by virtue of being defined ahead of time, may contain operators with little impact on the energy, but the classical optimizer will still need to perform a number of calls to the cost function in order to find their best values. Also, the gates originating from these operators, even if they are deemed unimportant because of a small associated parameter, will nevertheless be present in the circuit, adding to its depth. If high accuracy is sought, then ADAPT-VQE may require an ansatz comprised of a large number of operators, which in turn adds to the depth of the underlying circuit. More operators also mean more variational parameters, leading to an onerous optimization process. A more detailed analysis of this trade-off can be found in .
Algorithm for Defining a Virtual Space With a Small CI Hamiltonian
In this section, we present a downfolding method to define virtual orbitals for expanding the second-quantized Hamiltonian given in
Eq. 4. These new types of orbitals are able to capture significantly more correlation energy than the virtual orbitals coming from Hartree–Fock and one-electron Hamiltonians tested in
Section 2.2. The basis of this method is to define a set of virtual orbitals,
with
, which we call correlation optimized virtual orbitals or COVOs for short, by optimizing a small select CI Hamiltonian with respect to a single virtual orbital and then the next virtual orbitals in sequence, subject to them being orthonormal to the filled and previously computed virtual orbitals. The algorithm to calculate these new types of orbitals can be formulated as follows:
- 1.
Set .
- 2.
Using the ground state one-electron orbital, (or ground state orbitals for many-electron systems), and the virtual orbital to be optimized, , generate a CI matrix.
- 3.
Calculate the select CI expansion coefficients by diagonalizing the CI matrix.
- 4.
Using the CI coefficients associated with the lowest eigenvalue, calculate the gradient with respect to the and then update with a conjugate-gradient or similar method while making sure that is normalized and orthogonal to and for .
- 5.
If the gradient is small, then .
- 6.
If , go to step 2; otherwise, finish.
In the case of the H2 molecule, a small CI wave function for the 2 electron system composed of 2 one-electron orbitals, and , can be written as a linear combination of 6 determinant wave functions, or just 3 determinant wave functions for just singlet (or triplet) states,
Using this small CI ansatz, the energies of the system can be obtained by diagonalizing the following eigenvalue equation:where
Note that the overlap matrix, S, is the identity matrix for orthonormal and . The variation with respect to can be simply obtained using the following formula:
It should be noted that the above formulas can be generalized to work beyond two-electron systems by using corresponding orbitals techniques (; ). The next two Subsections 4.1–4.4 provide formulas that can be used to generate the matrix elements in Eq. 14 and the gradients with respect to in Eq. 15.
We also note that the COVO approach proposed in this work is similar in spirit to the optimized virtual orbital space (OVOS) approach developed over 30 years ago by Adamowicz and Bartlett (; ). The differences in our approach compared to this previous work is that the variational space used by COVOs is significantly bigger because plane-wave basis sets are used instead of LCAO Gaussian basis sets and that a second-order Hylleraas functional (; ; ; ; ) was used to describe the correlation in the OVOS procedure rather than a small CI Hamiltonian. Other differences with the COVOs approach are that the orbitals are optimized one at a time and the cost to generate them is similar to generating regular RHF virtual orbitals (just 4 to 9 times more expensive relative to RHF). Moreover, the resulting electronic gradient is non-Hermitian, which in addition to requiring more involved optimizers can result in extended energy plateaus that occur during the initial stages of the geodesic line searches in a conjugate-gradient or quasi-Newton optimization method.
One-Electron Orbitals for Two-State Hamiltonian
The four one-electron spin orbitals of two-state Hamiltonian are
where the spatial orbitals and spin functions are orthonormalized,
Two-Electron Orbitals for a Two-State Hamiltonian
For the two-state system, there are six two-electron wave functions, two of which are singlet, two are triplet, and two contain a mixture of singlet and triplet character. These wave functions can be written as
Note that and cannot be written as a product of a spatial wave function times a spin function. Moreover, these functions are not eigenfunctions of the spin operators and , and as a result, these determinants contain both singlet and triplet components. However, if we take linear combinations of them, we can get two new wave functions that are separable in spatial and spin functions and at the same time being eigenfunctions of and ,
Matrix Elements From the One-Electron Operators
The operator for H2 molecule iswhere is a function/operator of the coordinate r; i.e.,
Matrix Elements From the Two-Electron Operators
The operator for H2 molecule is
Results for 1 Ground State of the H2 Molecule Using Correlation Optimized Virtual Orbitals (COVOs)
The results for PW FCI calculations of H2 with 1, 4, 8, 12, and 18 COVOs are shown in Figure 4 and Table 2. The average difference error for the 1, 4, 8, and 12 COVOs calculations from the 18 COVOs calculation is 11.8 kcal/mol, 1.4 kcal/mol, 0.9 kcal/mol, and 0.3 kcal/mol, respectively. While the error is significant for 1 virtual, the difference is quite small by 4 virtual orbitals, and the error steadily decreases as the number of virtual orbitals increases. The error seen in the 4 optimized virtual orbitals’ calculations is similar to the 1.6 kcal/mol error seen in the DUCC calculations for the 19 virtual orbitals’ calculations in Section 2.2. Another measure of the error is the extensivity error. The energy for large R should be the same as the energy of twice the energy of an isolated H atom. For the pseudopotential plane-wave method being used, the energy of 2 H atoms is −0.997765 Hartrees (E(1H) = −0.498825 Hartrees). This difference at R = 7 Å is found to be 11.6, 1.2, 1.1, 0.5, and 0.4 kcal/mol for 1, 4, 8, 12, and 18 optimized virtual orbital calculations, respectively.
FIGURE 4
TABLE 2
| R (H-H) | PW FCI | PW FCI | PW VQE | PW FCI | PW FCI | PW FCI |
|---|---|---|---|---|---|---|
| (Å) | 1 COVO | 4 COVOs | 4 COVOs | 8 COVOs | 12 COVOs | 18 COVOs |
| 0.60 | −1.13749 | −1.15729 | −1.15728 | −1.15902 | −1.16028 | −1.16089 |
| 0.70 | −1.15321 | −1.17179 | −1.17178 | −1.17353 | −1.17467 | −1.17525 |
| 0.80 | −1.15128 | −1.16858 | −1.16857 | −1.17033 | −1.17136 | −1.17192 |
| 0.90 | −1.14124 | −1.15726 | −1.15724 | −1.15903 | −1.15995 | −1.16049 |
| 1.00 | −1.12742 | −1.14216 | −1.14213 | −1.14399 | −1.14478 | −1.14533 |
| 1.50 | −1.05311 | −1.06195 | −1.06195 | −1.06473 | −1.06516 | −1.06564 |
| 2.00 | −1.00793 | −1.01225 | −1.01220 | −1.01868 | −1.01916 | −1.01945 |
| 2.50 | −0.98862 | −1.00150 | −1.00150 | −1.00195 | −1.00228 | −1.00301 |
| 3.00 | −0.98137 | −0.99704 | −0.99701 | −0.99737 | −0.99789 | −0.99872 |
| 3.50 | −0.97883 | −0.99573 | −0.99570 | −0.99629 | −0.99698 | −0.99766 |
| 4.00 | −0.97810 | −0.99613 | −0.99611 | −0.99614 | −0.99693 | −0.99736 |
| 4.50 | −0.97817 | −0.99609 | −0.99608 | −0.99609 | −0.99722 | −0.99729 |
| 5.00 | −0.97845 | −0.99604 | −0.99598 | −0.99603 | −0.99716 | −0.99727 |
| 6.00 | −0.97906 | −0.99597 | −0.99596 | −0.99597 | −0.99705 | −0.99719 |
| 7.00 | −0.97928 | −0.99596 | −0.99596 | −0.99596 | −0.99703 | −0.99717 |
| NPE | 11.88 | 6.00 | 6.04 | 1.10 | 0.76 | — |
Total energies as a function of distance for the H2 molecule from plane-wave FCI calculations with 1, 4, 8, 12, and 18 COVOs and ADAPT-VQE simulations with 4 COVOs. Nonparallelity errors (NPE) are evaluated with respect to the calculations with the largest virtual orbital space (PW FCI 18 COVOs) and reported in milli-Hartree.
With only 4 optimized virtual orbitals, the correlation energy at the minimum was found to be −0.035 Hartrees, which is comparable to the −0.039 Hartrees found with CCSD/cc-pVTZ. The correlation energy decreases to −0.037, then −0.038, and finally −0.039 Hartrees as the number of optimized virtual orbitals increases to 8, 12, and 18 orbitals, respectively. These results showed that by 18 COVOs the same amount of correlation energy was recovered as with the cc-pVTZ LCAO basis set calculation. The results also showed that there was systematic convergence toward the benchmark LCAO result as the number of COVOs was increased, and with just 4 COVOs, a significant portion of the correlation energy was recovered.
Quantum Simulations of the 1 Ground State of the H2 Molecule Using COVOs
The previous section provides indisputable evidence for the performance of the proposed virtual orbitals for correlation energy recovery. Besides the possible ramifications in quantum chemistry carried out with classical computers, one immediate application is in the realm of quantum simulations. Because the present quantum hardware has not fully matured, hybrid algorithms that leverage classical resources and restrict the workload delegated to quantum computers, namely, state preparation and measurements of highly entangled states, are essential to meaningful quantum computations. The COVOs meet this requirement by decreasing the dimensionality of the problem, i.e., by enabling simulations with fewer qubits.
In order to probe the performance of COVOs in quantum simulations, we use the Hamiltonian with 4 COVOs and simulate the 1 ground state of H2 in the same bond distances shown in Table 2 and Figure 4. The ansatz circuit for the simulations is generated according to the ADAPT-VQE algorithm as implemented in the XACC (, ) framework for hybrid quantum computing using the tensor network quantum virtual machine () as the noiseless simulator backend (). In the present study, the ADAPT-VQE cycle is repeated until the norm of the gradient vector falls below 1e-2 and we use an operator pool containing all spin-adapted single and double excitation operators (one- and two-body rotations). A detailed account of ADAPT-VQE is exposed elsewhere (). Optimization of the parameterized circuit is conducted with the COBYLA (Powell, 1994) optimizer as implemented in the NLOpt package (). Results for the simulated potential energy curve are plotted in Figure 5.
FIGURE 5
It is evident from Figure 5 that ADAPT-VQE can generate a circuit capable of reproducing the FCI results in the current active space. These simulations deliver a smooth, continuous potential energy curve that tracks the FCI values strikingly well. The deviations from the corresponding FCI energies are all found below 1e-4 Hartree. This means that not only do these simulations deliver results that are well below the conventional chemical accuracy mark, but also more importantly in the current context is that this error is inconsequential compared to the effect of noise in case of deployment on actual quantum hardware.
It is remarkable that the results degrade little throughout the energy scan, which attests to the aptness and flexibility of ADAPT-VQE in determining an ansatz according to the complexity of the underlying electronic structure. The ansatz in the vicinity of the equilibrium bond length 0.5–1.0 Å is comprised solely by pair excitations as would be expected given a restricted HF reference, which means no determinant obtained via one-body rotations can lower the energy below that of HF. As we approach the Coulson-Fischer point (), single excitations start to become part of the ansatz, which signals the inadequacy of a restricted reference wave function and that inclusion of these operators enables the ansatz to remain in the 1 potential energy curve, which means that this flexibility may come at the expense of deeper circuits. Because one-qubit gates tend to be executed in a short timescale and are fairly insensitive to noise, we can use the number of CNOTs present in the circuit as indicative of the complexity in its implementation, which we provide in Figure 6, showing that the ansatzes generated from ADAPT-VQE are much more affordable than those obtained by ordinary UCCSD VQE simulations.
FIGURE 6
Along these lines, once the operator composition of the ansatz is defined, by virtue of introducing more parameters, we are likely to experience a more arduous optimization of the corresponding parameterized gates. This has a compound effect with the circuit depth since more measurements are needed, each of which requires the circuit to be implemented and measurements to take place. Figure 7 gives a profile of the optimization performance along the potential energy scan.
FIGURE 7
It should come as no surprise that the optimization is more difficult in the regime of stronger correlation. This region also demands a more complex ansatz, as the top plot in Figure 7 shows that only in this vicinity (1.5–3.5 Å) we observe ansatzes with more than four operators. Interestingly, the number of objective function calls does not show large deviations for ansatzes with 1–3 operators, regardless of where they are found in the potential energy curve, which is further corroborated by the relatively small error bars in the corresponding columns of the bottom plot. This observation does not hold as more parameters/operators are introduced in the ansatz in order to accommodate a more complex electronic structure. Thus, with four parameters, not only are more calls to the objective function needed, but also there is a more pronounced standard deviation. Ansatzes with five or more operators can only be found in the (1.5–3.5 Å), as we can see that the calls to the objective function coming from them dominate the overall number of optimization cycles. Due to the scarce occurrence of these ansatzes in the current energy scan, the corresponding statistical information that can be derived from these instances is not as reliable. All in all, this plot is valuable in lending additional insight into the resources required to perform these simulations. It is important to mention that, for every new ansatz, the variational parameters are initialized at zero. Alternatively, the parameters corresponding to the previously optimizer ansatzes could be initialized at their optimal values and the new parameter would be introduced in the ansatz, which would accelerate convergence. Moreover, the convergence profile likely displays pronounced dependence on the chosen optimizer, which is not pursued here.
Conclusion
In summary, we have developed a new approach for defining virtual spaces with a pseudopotential plane-wave code for use in many-body methods described by second-quantized Hamiltonians. The method is based on optimizing the virtual orbitals to minimize a small select CI Hamiltonian (i.e., COVOs) that contains configurations containing filled RHF orbitals and the one virtual orbital to be optimized. Subsequent virtual orbitals are optimized in the same way, but with the added constraint of being orthogonal to the previously calculated filled and virtual orbitals. The method was applied to the simple, but nontrivial, H2 molecule. As summarized in Figure 8, these new types of virtual orbitals were significantly better at capturing correlation in plane-wave calculations than from virtual spaces from Hartree–Fock and one-electron Hamiltonian, and moreover, we were able to obtain good agreement with Gaussian cc-pVTZ basis set results with just 4 virtual orbitals for the H2 molecule. Subsequent calculations showed that the correlation energy converged steadily as more virtual orbitals were included in the calculation. With 18 virtual orbitals, the correlation energies were found to be converged to less than 0.5 kcal/mol. The robustness of the proposed basis sets is corroborated by its ready applicability to quantum simulations, which in the case of ADAPT-VQE show remarkable agreement with the classical, exact diagonalization result (FCI) in the same basis set (4 COVOs).
FIGURE 8
Because this study is focused on how one might carry out plane-wave CI calculations on near-term quantum computers in the next few years, we have only shown results for the H2 dimer. However, we are optimistic that these correlation optimized virtual orbitals open up the door to many-body calculations using pseudopotential plane-wave calculations, including coupled cluster, Møller–Plesset, and Green’s function theories as well as other FCI-approaching methods for quantum computers. We hope in future studies to more thoroughly test the effectiveness of the COVOs procedure on larger and more complicated molecules and materials. To lend credence to this assertion, we show the promising results for Be2 dimer with a small number of COVOs in Figure 9. Also as shown in Figure 10, the shapes of the COVOs end up being similar to what is found for the virtual orbitals from LCAO calculations. This suggests that new classes of LCAO basis sets might be able to be generated using a simple rotation of the filled orbitals and COVOs. Future work will focus on using this approach on larger molecular and periodic systems. With the validation granted by our quantum simulations, further studies are called for, including the further reduction of dimension by employing active-space DUCC downfolded Hamiltonians, OVOS, and natural orbitals, as well as work in conjunction with VQE methods.
FIGURE 9
FIGURE 10
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://github.com/ebylaska/PWH2-Data/archive/master.zip.
Author contributions
All authors contributed to the writing and editing of the manuscript. EB and DS developed and implemented in the NWChem the algorithm for defining a virtual space with small CI Hamiltonians. KK wrote a part of the NWChem-DUCC interface and codeveloped the original form of the DUCC formalism. NB wrote a part of the NWChem-DUCC interface, codeveloped the original form of the DUCC formalism, carried out FCI and QPE calculations, and prepared formatted integrals for VQE calculations. DC wrote the algorithm for and performed the VQE quantum simulations.
Funding
This material is based upon work supported by the United States Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division through its “Embedding Quantum Computing into Many-body Frameworks for Strongly Correlated Molecular and Materials Systems” project at Pacific Northwest National Laboratory (PNNL). This work was also supported by the Quantum Science Center (QSC), a National Quantum Information Science Research Center of the United States Department of Energy (DOE). We also would like to thank the DOE BES Chem CCS, DOE BES Geochemistry, and DOE Advanced Scientific Computing Research (ASCR) ECP NWChemEx programs for their support of software development for high-performance computers and computer time needed to carry out the work. PNNL is operated for the United States Department of Energy by the Battelle Memorial Institute under Contract DE-AC06-76RLO-1830. This research was also partially supported, thru their support of software development for high-performance computers, by DOE BES Chem CCS and DOE BES Geosciences programs, as well as the Exascale Computing Project (17-SC-20-SC), a collaborative effort of the United States Department of Energy Office of Science and the National Nuclear Security Administration. Calculations have been performed using computational resources at the Pacific Northwest National Laboratory (PNNL). This research used resources of the National Energy Research Scientific Computing Center (NERSC), a User Facility supported by the Office of Science of the United States DOE under Contract No. DE-AC02-05CH11231, and the Argonne ALCF computing center through their early science program. This research also used resources of the Oak Ridge Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC05-00OR22725. This manuscript has been authored in part by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725 with the United States Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan.
Acknowledgments
We would like to thank the NWChem project team and the people that have helped the progress of the NWChem software over the years. DC would like to thank Alexander McCaskey for discussions and help with the software engineering of the quantum algorithm used in this paper.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AdamowiczL.BartlettR. J. (1987). Optimized virtual orbital space for high‐level correlated calculations. J. Chem. Phys.86, 6314–6324. 10.1063/1.452468
2
AdamowiczL.BartlettR. J.SadlejA. J. (1988). Optimized virtual orbital space for high‐level correlated calculations. II. Electric properties. J. Chem. Phys.88, 5749–5758. 10.1063/1.454721
3
AllenL. C.KaroA. M. (1962). Electronic structure of simple Molecules1. J. Phys. Chem.66, 2329–2331. 10.1021/j100818a011
4
ApràE.BylaskaE. J.De JongW. A.GovindN.KowalskiK.StraatsmaT. P.et al (2020). NWChem: past, present, and future. J. Chem. Phys.152, 184102. 10.1063/5.0004997
5
BaumanN. P.BylaskaE. J.KrishnamoorthyS.LowG. H.WiebeN.GranadeC. E.et al (2019). Downfolding of many-body Hamiltonians using active-space models: extension of the sub-system embedding sub-algebras approach to unitary coupled cluster formalisms. J. Chem. Phys.151, 014107. 10.1063/1.5094643
6
BetteridgeP. W.CarruthersJ. R.CooperR. I.ProutK.WatkinD. J. (2003). CRYSTALSversion 12: software for guided crystal structure analysis. J. Appl. Crystallogr.36, 1487. 10.1107/S0021889803021800
7
BoysS. F.CookG. B.ReevesC. M.ShavittI. (1956). Automatic fundamental calculations of molecular structure. Nature178, 1207–1209. 10.1038/1781207a0
8
BoysS. F. (1950). Electronic wave functions-i. a general method of calculation for the stationary states of any molecular system. Proc. R. Soc. Lond. A Math. Phys. Sci.200, 542–554. 10.1098/rspa.1950.0036
9
BravyiS. B.KitaevA. Y. (2002). Fermionic quantum computation. Ann. Phys.298, 210–226. 10.1006/aphy.2002.6254
10
BriggsE.SullivanD.BernholcJ. (1996). Real-space multigrid-based approach to large-scale electronic structure calculations. Phys. Rev. B Condens. Matter54, 14362. 10.1103/PhysRevB.54.14362
11
BylaskaE. J.RossoK. (2018). Corresponding orbitals derived from periodic bloch states for electron transfer calculations of transition metal oxides. J. Chem. Theor. Comput.14, 4416–4426. 10.1021/acs.jctc.7b01180
12
BylaskaE. J.ApràE.KowalskiK.JacquelinM.De JongW. A.VishnuA.et al (2017a). “Transitioning NWChem to the next generation of manycore machines,” in Exascale scientific applications: scalability and performance portability. Editors StraatsmaT. P.AntypasK. B.Timothy WilliamsT. J. J. (Boca Raton, FL: CRC Press), 165. 10.1201/b21930-8
13
BylaskaE. J.JacquelinM.de JongW. A.HammondJ. R.KlemmM. (2017b). “Performance evaluation of NWChem ab-initio molecular dynamics (AIMD) simulations on the Intel® xeon Phi™ processor,” in High performance computing. Editors KunkelJ. M.YokotaR.TauferM.ShalfJ. (Cham, Switzerland: Springer International Publishing), 404–418. 10.1007/978-3-319-67630-2_30
14
BylaskaE. J. (2017). “Plane-wave DFT methods for chemistry,” in Annuual Reports in Computational Chemistry. Editor DixonD. A. (Amsterdam, Netherlands: Elsevier), 13, 185–228. 10.1016/bs.arcc.2017.06.006
15
BylaskaE. J.WatersK.HermesE. D.ZádorJ.RossoK. M. (2020). A filon-like integration strategy for calculating exact exchange in periodic boundary conditions: a plane-wave dft implementation. Mater. Theory.4, 1–31. 10.1186/s41313-020-00019-9
16
BylaskaE.TsemekhmanK.GovindN.ValievM. (2011). “Large-scale plane-wave-based density functional theory: formalism, parallelization, and applications,” in Computational methods for large systems: electronic structure approaches for biotechnology and nanotechnology. Editor ReimersJ. R. (Hoboken, NJ: Wiley Online Library), 77–116. 10.1002/9780470930779.ch3
17
ChallacombeM. (2000). Linear scaling computation of the fock matrix. v. hierarchical cubature for numerical integration of the exchange-correlation matrix. J. Chem. Phys.113, 10037–10043. 10.1063/1.1316012
18
ChenM.GuoG.HeL. (2011). Electronic structure interpolation via atomic orbitals. J. Phys. Condens. Matter23, 325501. 10.1088/0953-8984/23/32/325501
19
ClarkS. J.SegallM. D.PickardC. J.HasnipP. J.ProbertM. I. J.RefsonK.et al (2005). First principles methods using castep. Z. für Kristallogr. - Cryst. Mater.220, 567–570. 10.1524/zkri.220.5.567.65075
20
CoulsonC. A.FischerI. (1949). Notes on the molecular orbital treatment of the hydrogen molecule. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science.40, 386. 10.1080/14786444908521726
21
CrawfordT. D.ClementiC.HarrisonR.Head-GordonT.JhaS.KrylovA.et al (2017). The molecular sciences software institute. Abstr. Pap. Am. Chem. Soc.Vol. 254. WASHINGTON, DC: AMER CHEMICAL SOC 1155 16TH ST, NW.
22
DeslippeJ.SamsonidzeG.StrubbeD. A.JainM.CohenM. L.LouieS. G. (2012). Berkeleygw: a massively parallel computer package for the calculation of the quasiparticle and optical properties of materials and nanostructures. Comput. Phys. Commun.183, 1269–1289. 10.1016/j.cpc.2011.12.006
23
DubashM. (2005). Moore’s law is dead, says gordon moore. Techworld. com13.
24
DunningT. H. (1989). Gaussian basis sets for use in correlated molecular calculations. i. the atoms boron through neon and hydrogen. J. Chem. Phys.90, 1007–1023. 10.1063/1.456153
25
DunningT. H.HayP. J. (1977). Gaussian basis sets for molecular calculations. Methods of electronic structure theory 3, 1–27. 10.1007/978-1-4757-0887-5_1
26
DupuisM.WattsJ. D.VillarH. O.HurstG. J. B. (1989). The general atomic and molecular electronic structure system hondo: version 7.0. Comput. Phys. Commun.52, 415–425. 10.1016/0010-4655(89)90116-1
27
EvangelistaF. A.ChanG. K.ScuseriaG. E. (2019). Exact parameterization of fermionic wave functions via unitary coupled cluster theory. J. Chem. Phys.151, 244112. 10.1063/1.5133059
28
FrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.et al (2016). Gaussian 16 revision C.01. Wallingford CT: Gaussian Inc.
29
GiannozziP.BaroniS.BoniniN.CalandraM.CarR.CavazzoniC.et al (2009). Quantum espresso: a modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter21, 395502. 10.1088/0953-8984/21/39/395502
30
GoldbeckG. (2017). The scientific software industry: a general overview. Goldbeck Consulting, Technical Report.
31
GonzeX.JolletF.AraujoF. A.AdamsD.AmadonB.ApplencourtT.et al (2016). Recent developments in the abinit software package. Comput. Phys. Commun.205, 106–131. 10.1016/j.cpc.2016.04.003
32
GrimsleyH. R.ClaudinoD.EconomouS. E.BarnesE.MayhallN. J. (2020). Is the trotterized uccsd ansatz chemically well-defined?J. Chem. Theor. Comput.16, 1–6. 10.1021/acs.jctc.9b01083
33
GrimsleyH. R.EconomouS. E.BarnesE.MayhallN. J. (2019). An adaptive variational algorithm for exact molecular simulations on a quantum computer. Nat. Commun.10, 3007. 10.1038/s41467-019-10988-2
34
GygiF. (2008). Architecture of qbox: a scalable first-principles molecular dynamics code. IBM J. Res. Dev.52, 137–144. 10.1147/rd.521.0137
35
HallG. (1951). The molecular orbital theory of chemical valency viii. a method of calculating ionization potentials. Proc. R. Soc. Lond. A Math. Phys. Sci.205, 541–552. 10.1098/rspa.1951.0048
36
HamannD. R. (1989). Generalized norm-conserving pseudopotentials. Phys. Rev. B Condens. Matter40, 2980. 10.1103/physrevb.40.2980
37
HarrisonR. J.BeylkinG.BischoffF. A.CalvinJ. A.FannG. I.Fosso-TandeJ.et al (2016). Madness: a multiresolution, adaptive numerical environment for scientific simulation. SIAM J. Sci. Comput.38, S123–S142. 10.1137/15M1026171
38
HocquetA.WieberF. (2017). “only the initiates will have the secrets revealed”: computational chemists and the openness of scientific software. IEEE Ann. Hist. Comput.39, 40–58. 10.1109/MAHC.2018.1221048
39
HutterJ.IannuzziM.SchiffmannF.VandeVondeleJ. (2014). cp2k: atomistic simulations of condensed matter systems. WIREs Comput. Mol. Sci.4, 15–25. 10.1002/wcms.1159
40
HylleraasE. A. (1929). Neue berechnung der energie des heliums im grundzustande, sowie des tiefsten terms von ortho-helium. Z. Phys.54, 347–366. 10.1007/bf01375457
41
HylleraasE. A. (1964). The schrödinger two-electron atomic problem. Adv. Quant. Chem.1, 1–33. 10.1016/s0065-3276(08)60373-1
42
HylleraasE. A. (1930). Über den Grundterm der Zweielektronenprobleme von H?, He, Li+, Be++ usw. Z. Phys.65, 209–225. 10.1007/bf01397032
43
HylleraasE. A. (1928). Über den Grundzustand des Heliumatoms. Z. Phys.48, 469–494. 10.1007/bf01340013
44
JohnsonS. (2007). The NLopt nonlinear-optimization package. Available at: http://github.com/stevengj/nlopt.10.1364/np.2007.ntub1
45
JordanP.WignerE. (1928). Über das Paulische äquivalenzverbot. Z. Phys.47, 631–651. 10.1007/BF01331938
46
KandalaA.MezzacapoA.TemmeK.TakitaM.BrinkM.ChowJ. M.et al (2017). Hardware-efficient variational quantum eigensolver for small molecules and quantum magnets. Nature549, 242–246. 10.1038/nature23879
47
KendallR. A.ApràE.BernholdtD. E.BylaskaE. J.DupuisM.FannG. I.et al (2000). High performance computational chemistry: an overview of nwchem a distributed parallel application. Comput. Phys. Commun.128, 260–283. 10.1016/S0010-4655(00)00065-5
48
KingH. F.StantonR. E.KimH.WyattR. E.ParrR. G. (1967). Corresponding orbitals and the nonorthogonality problem in molecular quantum mechanics. J. Chem. Phys.47, 1936–1941. 10.1063/1.1712221
49
KleinmanL.BylanderD. M. (1982). Efficacious form for model pseudopotentials. Phys. Rev. Lett.48, 1425. 10.1103/physrevlett.48.1425
50
KogaT. (1992). Hylleraas wave functions revisited. J. Chem. Phys.96, 1276–1279. 10.1063/1.462164
51
KohnW.ShamL. J. (1965). Self-consistent equations including exchange and correlation effects. Phys. Rev.140, A1133. 10.1103/PhysRev.140.A1133
52
KowalskiK.BaumanN. P. (2020). Publisher's Note: “Sub-system quantum dynamics using coupled cluster downfolding techniques” [J. Chem. Phys. 152, 244127 (2020)]. J. Chem. Phys.153, 049902. 10.1063/5.000843610.1063/5.0021137
53
KresseG.FurthmüllerJ. (1996). Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci.6, 15–50. 10.1016/0927-0256(96)00008-0
54
LowG. H.BaumanN. P.GranadeC. E.PengB.WiebeN.BylaskaE. J.et al (2019). Q# and nwchem: tools for scalable quantum chemistry on quantum computers. Preprint: arXiv:1904.01131.
55
McCaskeyA.DumitrescuE.ChenM.LyakhD.HumbleT. (2018a). Validating quantum-classical programming models with tensor network simulations. PloS One13, e0206704–19. 10.1371/journal.pone.0206704
56
McCaskeyA. J.DumitrescuE. F.LiakhD.ChenM.FengW.HumbleT. S. (2018b). A language and hardware independent approach to quantum-classical computing. Software7, 245–254. 10.1016/j.softx.2018.07.007
57
McCaskeyA. J.LyakhD. I.DumitrescuE. F.PowersS. S.HumbleT. S. (2020). XACC: a system-level software infrastructure for heterogeneous quantum-classical computing. Quantum Sci. Technol.5, 024002. 10.1088/2058-9565/ab6bf6
58
McCaskeyA. J.ParksZ. P.JakowskiJ.MooreS. V.MorrisT. D.HumbleT. S.et al (2019). Quantum chemistry as a benchmark for near-term quantum computers. Npj Quantum Inf.5, 1–8. 10.1038/s41534-019-0209-0
59
MooreG. E. (2006). Cramming more components onto integrated circuits, Reprinted from Electronics, volume 38, number 8, April 19, 1965, pp.114 ff. IEEE Solid-State Circuits Soc. Newsl.11, 33–35. 10.1109/N-SSC.2006.4785860
60
MullikenR. S.RiekeC. A.BrownW. G. (1941). Hyperconjugation. J. Am. Chem. Soc.63, 41–56. 10.1021/ja01846a008
61
MullikenR. S.RiekeC. A. (1941). Improved computations on conjugation and hyperconjugation. J. Am. Chem. Soc.63, 1770–1771. 10.1021/ja01851a514
62
NeeseF. (2018). Software update: the orca program system, version 4.0. WIREs Comput Mol Sci.8, e1327. 10.1002/wcms.1327
63
NesbetR. K. (1963). Computer programs for electronic wave-function calculations. Rev. Mod. Phys.35, 552. 10.1103/RevModPhys.35.552
64
NesbetR. K. (1960). Ground state electronic wave function of methane. J. Chem. Phys.32, 1114–1122. 10.1063/1.1730858
65
O’MalleyP. J. J.BabbushR.KivlichanI. D.RomeroJ.McCleanJ. R.BarendsR.et al (2016). Scalable quantum simulation of molecular energies. Phys. Rev. X.6, 031007. 10.1103/PhysRevX.6.031007
66
ParrR. G.MullikenR. S. (1950). LCAOSelf‐consistent field calculation of the π‐electron energy levels of cis‐ andtrans−1,3‐butadiene. J. Chem. Phys.18, 1338–1346. 10.1063/1.1747474
67
PeruzzoA.McCleanJ.ShadboltP.YungM. H.ZhouX. Q.LoveP. J.et al (2014). A variational eigenvalue solver on a photonic quantum processor. Nat. Commun.5, 4213. 10.1038/ncomms5213
68
PopleJ. A.SantryD. P.SegalG. A. (1965). Approximate self‐consistent molecular orbital theory. I. Invariant procedures. J. Chem. Phys.43, S129–S135. 10.1063/1.1701475
69
PowellM. J. D. (1994). A direct search optimization method that models the objective and constraint functions by linear interpolation. Dordrecht, Netherlands: Springer Netherlands, 51–67. 10.1007/978-94-015-8330-5_4
70
PrendergastD.LouieS. G. (2009). Bloch-state-based interpolation: an efficient generalization of the shirley approach to interpolating electronic structure. Phys. Rev. B80, 235126. 10.1103/PhysRevB.80.235126
71
PreskillJ. (2018). Quantum computing in the nisq era and beyond. Quantum.2, 79. 10.22331/q-2018-08-06-79
72
PulayP. (1969). Ab initiocalculation of force constants and equilibrium geometries in polyatomic molecules. Mol. Phys.17, 197–204. 10.1080/00268976900100941
73
ReevesC. M. (1966). An algorithm for generating projective reduction formulas for matrix elements of many-electron wavefunctions. Commun. ACM9, 276–279. 10.1145/365278.365540
74
ReiherM.WiebeN.SvoreK. M.WeckerD.TroyerM. (2017). Elucidating reaction mechanisms on quantum computers. Proc. Natl. Acad. Sci. U.S.A.114, 7555–7560. 10.1073/pnas.1619152114
75
RichardR. M.BertoniC.BoschenJ. S.KeipertK.PritchardB.ValeevE. F.et al (2019). Developing a computational chemistry framework for the exascale era. Comput. Sci. Eng.21, 48–58. 10.1109/MCSE.2018.2884921
76
RomeroJ.BabbushR.McCleanJ. R.HempelC.LoveP. J.Aspuru-GuzikA. (2017). Strategies for quantum computing molecular energies using the unitary coupled cluster ansatz. Quantum Sci. Technol.4, 014008. 10.1088/2058-9565/aad3e4
77
RoothaanC. C. J. (1951). New developments in molecular orbital theory. Rev. Mod. Phys.23, 69. 10.1103/RevModPhys.23.69
78
RotmanD. (2020). We’re not prepared for the end of moore’s law. MIT Technology Review.
79
SchmidtM. W.BaldridgeK. K.BoatzJ. A.ElbertS. T.GordonM. S.JensenJ. H.et al (1993). General atomic and molecular electronic structure system. J. Comput. Chem.14, 1347–1363. 10.1002/jcc.540141112
80
SetiaK.WhitfieldJ. D. (2018). Bravyi-kitaev superfast simulation of electronic structure on a quantum computer. J. Chem. Phys.148, 164104. 10.1063/1.5019371
81
ShaoY.GanZ.EpifanovskyE.GilbertA. T.WormitM.KussmannJ.et al (2015). Advances in molecular quantum chemistry contained in the q-chem 4 program package. Mol. Phys.113, 184–215. 10.1080/00268976.2014.952696
82
ShirleyE. L. (1996). Optimal basis sets for detailed brillouin-zone integrations. Phys. Rev. B Condens. Matter.54, 16464. 10.1103/PhysRevB.54.16464
83
StantonJ. F.GaussJ. r.WattsJ. D.LauderdaleW. J.BartlettR. J. (1992). The aces ii program system. Int. J. Quant. Chem.44, 879–894. 10.1002/qua.560440876
84
SteudtnerM.WehnerS. (2018). Lowering qubit requirements for quantum simulations of fermionic systems. Preprint: arXiv:1712.07067v2.
85
SvoreK.RoettelerM.GellerA.TroyerM.AzariahJ.GranadeC.et al (2018). “Q# enabling scalable quantum computing and development with a high-level dsl,” in Proceedings of the Real World Domain Specific Languages Workshop 2018, Vienna, Austria, February 2018. 1–10. 10.1145/3183895.3183901
86
te VeldeG.BickelhauptF. M.BaerendsE. J.Fonseca GuerraC.van GisbergenS. J. A.SnijdersJ. G.et al (2001). Chemistry with adf. J. Comput. Chem.22, 931–967. 10.1002/jcc.1056
87
ValievM.BylaskaE. J.GovindN.KowalskiK.StraatsmaT. P.Van DamH. J. J.et al (2010). NWChem: a comprehensive and scalable open-source solution for large scale molecular simulations. Comput. Phys. Commun.181, 1477–1489. 10.1016/j.cpc.2010.04.018
88
van DamH.ApraE.BairR.BoschenJ.BylaskaE.De JongW.et al (2020). NWChemEx–computational chemistry for the exascale era. Bull. Am. Phys. Soc.65.
89
VermaP.HuntingtonL.CoonsM.KawashimaY.YamazakiT.ZaribafiyanA. (2020). Scaling up electronic structure calculations on quantum computers: the frozen natural orbital based method of increments. Preprint: arXiv:2002.07901.
90
WasielewskiM. R.ForbesM. D. E.FrankN. L.KowalskiK.ScholesG. D.Yuen-ZhouJ.et al (2020). Exploiting chemistry and molecular systems for quantum information science. Nat. Rev. Chem.4, 490–15. 10.1038/s41570-020-0200-5
91
WernerH.-J.KnowlesP. J.KniziaG.ManbyF. R.SchützM. (2012). Molpro: a general-purpose quantum chemistry program package. WIREs Comput. Mol. Sci.2, 242–253. 10.1002/wcms.82
92
Wikipedia, The Free Encyclopedia (2020). List of quantum chemistry and solid-state physics software. San Francisco: Wikimedia Foundation Inc.
93
Wilkins-DiehrN.CrawfordT. D. (2018). NSF's inaugural software institutes: the science gateways community institute and the molecular sciences software institute. Comput. Sci. Eng.20, 26–38. 10.1109/MCSE.2018.05329813
Summary
Keywords
quantum computing, ADAPT-VQE, many-body calculations, DUCC, pseudopotential plane-wave, correlation optimized virtual orbitals, COVOs, nwchem
Citation
Bylaska EJ, Song D, Bauman NP, Kowalski K, Claudino D and Humble TS (2021) Quantum Solvers for Plane-Wave Hamiltonians: Abridging Virtual Spaces Through the Optimization of Pairwise Correlations. Front. Chem. 9:603019. doi: 10.3389/fchem.2021.603019
Received
04 September 2020
Accepted
11 January 2021
Published
18 March 2021
Volume
9 - 2021
Edited by
Javier Carrasco, CIC energigune, Spain
Reviewed by
Rajat Chaudhuri, Indian Institute of Astrophysics, India
Masato Kobayashi, Hokkaido University, Japan
Updates
Copyright
© 2021 Bylaska, Song, Bauman, Kowalski, Claudino and Humble.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eric J. Bylaska, eric.bylaska@pnnl.gov
This article was submitted to Physical Chemistry and Chemical Physics, a section of the journal Frontiers in Chemistry
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.