# Shortcut to Chemically Accurate Quantum Computing
via Density-based Basis-set Correction.

D. Traore, O. Adjoua, C. Feniou, I.-M. Lygatsika, Y. Maday, E. Posenitskiy,
K. Hammernik, A. Peruzzo, J. Toulouse, E. Giner, et al.

► To cite this version:

D. Traore, O. Adjoua, C. Feniou, I.-M. Lygatsika, Y. Maday, et al. .. Shortcut to Chemically Accurate Quantum Computing via Density-based Basis-set Correction.. Communications Chemistry, 2024, 7, pp.269. 10.1038/s42004-024-01348-3 . hal-04579665

HAL Id: hal-04579665

https://hal.science/hal-04579665v1

Submitted on 18 Nov 2024

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Shortcut to chemically accurate quantum computing via density-based basis-set correction

<figure><img src="image_3.png" /></figure>

Check for updates

Diata Traore<fnref n="1" />,², Olivier Adjoua<fnref n="1" />, César Feniou<fnref n="1" />,², Ioanna-Maria Lygatsika<fnref n="1" />,³,<fnref n="6" />,⁷, Yvon Maday³,⁴, Evgeny Posenitskiy², Kerstin Hammernik Ⓔ⁵, Alberto Peruzzo², Julien Toulouse Ⓔ<fnref n="1" />,⁴, Emmanuel Giner<fnref n="1" /> & Jean-Philip Piquemal Ⓔ<fnref n="1" />,² ☒

Using GPU-accelerated state-vector emulation, we propose to embed a quantum computing ansatz into density-functional theory via density-based basis-set corrections to obtain quantitative quantum-chemistry results on molecules that would otherwise require brute-force quantum calculations using hundreds of logical qubits. Indeed, accessing a quantitative description of chemical systems while minimizing quantum resources is an essential challenge given the limited qubit capabilities of current quantum processors. We provide a shortcut towards chemically accurate quantum computations by approaching the complete-basis-set limit through coupling the density-based basis-set corrections approach, applied to any given variational ansatz, to an on-the-fly crafting of basis sets specifically adapted to a given system and user-defined qubit budget. The resulting approach self-consistently accelerates the basis-set convergence, improving electronic densities, ground-state energies, and first-order properties (e.g. dipole moments), but can also serve as a classical, a posteriori, energy correction to quantum hardware calculations with expected applications in drug design and materials science.

Quantum computing (QC) offers a promising approach to solving electronic-structure problems, with algorithms like quantum phase estimation (QPE) and variational quantum eigensolver (VQE) proving effective in performing ground-state quantum-chemistry wave-function calculations<fnref n="1" />⁻<fnref n="6" />. The electronic Hamiltonian is expressed in second quantization, employing an encoding that maps one spin-orbital to one qubit. This mapping allows one to represent an exponentially large Hilbert space using only a linear number of qubits. However, to achieve accurate and practically valuable predictions for chemical systems in real-world applications, the molecular Hamiltonian should be expressed and solved using extensive basis sets of one-electron orbital functions. The number of qubits required for such calculations quickly exceeds the available capacities on current noisy intermediate scale quantum (NISQ) devices, upcoming early fault-tolerant quantum computing (FTQC) devices, and high-performance classical emulators. Therefore, while quantum chemistry has long been stated as a promising application for quantum computing, the endeavors have so far been limited to small molecular systems and minimal basis sets.

However, in practice, minimal basis sets fail to be predictive for ground-state energies, and chemically useful calculations require at least significantly larger than double-zeta basis sets. Requirements for computing molecular properties are even more drastic as they tend to converge more slowly with the size of the basis set than energies.

In electronic-structure theory, the exact solution is defined by the full-configuration interaction (FCI) method in the (infinite) complete-basis-set (CBS) limit. However, in practice, the solution derived from a finite basis set is used, which inherently suffers from truncation errors. These errors can be substantial for small basis sets, but for large enough basis sets it is possible to reach the target of chemical accuracy on energy differences, i.e., 1 kcal/mol (1.6 mHa). Unfortunately, employing sufficiently large basis sets becomes very rapidly impractical for large systems, particularly those with more than a few dozen atoms. As a result, quantum chemists have developed a variety of traditional (i.e., classical) computational methods to approach chemical accuracy at a reasonable cost⁷⁻¹⁰. Similarly, on the quantum computing side, the diversity of available methods has significantly increased¹¹⁻¹⁴.

<figure><img src="image_4.png" /></figure>

e-mail: jean-philip.piquemal@sorbonne-universite.fr

<footnote n="1">¹Sorbonne Université, LCT, UMR 7616 CNRS, 75005 Paris, France. ²Qubit Pharmaceuticals, Advanced Research Department, 75014 Paris, France. ³Sorbonne Université, LJLL, UMR 7598 CNRS, 75005 Paris, France. ⁴Institut Universitaire de France, 75005 Paris, France. ⁵NVIDIA Corporation, Santa Clara, CA, USA.</footnote>

<footnote n="6">⁶Present address: CEA, DAM, DIF, F-91297 Arpajon, France. ⁷Present address: Université Paris-Saclay, LMCE, 91680 Bruyères-le-Châtel, France.</footnote>
In particular, the hybrid quantum-classical VQE techniques initially designed to solve the general eigenvalue problem³, have been shown to be particularly suited for chemical applications in the present NISQ era⁵,¹⁵,¹⁵,¹⁶. VQE algorithms have evolved over the years in two directions: (i) fixed-length ansätze, often inspired by classical coupled cluster, with approaches such as the unitary coupled cluster (UCC) and its extensions (see refs. 17–19 and references therein); (ii) adaptive methods such as the adaptive derivative-assembled pseudo-trotter ansatz variational quantum eigen-solver (ADAPT-VQE)²⁰ have also started to be particularly popular since they allow tailoring system-specific ansätze with shorter circuits. In recent years, ADAPT-VQE approaches have been systematically improved²¹⁻²⁶ while alternative, resource-saving, adaptive techniques have also been introduced²⁷,²⁸.

For applications to quantum computing, the inspiration coming from classical approaches has not been limited¹⁴ to UCC-like techniques, and one can particularly mention the importance of explicitly correlated approaches²⁹,³⁰. Indeed, several quantum strategies originate from the latter and aim to maintain full accuracy while further minimizing quantum resources for obtaining more compact wave-function representations, which is the key criterion for hardware implementations. We can particularly mention the works on transcorrelated approaches³¹⁻³⁶, which introduce short-range correlation effects such as the electron–electron cusp condition³⁷ through various Hamiltonian modifications. Therefore, explicitly correlated quantum computing ansätze based on the transcorrelated approach has been shown to successfully evaluate ground-state energies while using fewer resources³⁸⁻⁴². Some of them have even been extended to the evaluation of excited-state energies³⁰. However, other strategies exist in classical quantum chemistry and could be particularly suited to quantum computing.

In particular, basis-set correction techniques such as the density-based basis-set correction (DBBSC) method⁴³⁻⁵⁵, relying on density-functional theory (DFT), have proven to be effective for calculating ground- and excited-state energies, and also dipole moments, for a variety of systems, including atoms, small organic molecules, and strongly correlated systems. This approach provides the key benefit of an accelerated convergence to the CBS limit with the basis-set size, which is a valuable asset for quantum computing since minimizing quantum resources is paramount. To date, the DBBSC method has only been applied to still relatively large basis sets, more specifically the family of Dunning basis sets⁵⁶ for which chemical accuracy with respect to the CBS limit can be reached starting from a triple-zeta basis set. We note that the approximations developed within the DBBSC method target the CBS limit within a given wave-function ansatz. This means that they do not address the intrinsic errors of the wave-function method itself, such as the effects of neglected higher excitations in truncated configuration-interaction (CI) or coupled-cluster methods.

In this work, we propose the integration of the DBBSC method with quantum algorithms to expedite reaching the CBS limit and achieve chemical accuracy on complex molecular systems. This strategy, which natively limits the required qubit counts, can be applied to quantum algorithms that tackle the ground-state quantum chemistry problem, such as the QPE or VQE algorithms. While QPE can guarantee ground-state energy with arbitrarily high precision given a carefully chosen initial state, it demands large quantum circuits and will only be viable in the FTQC era. Conversely, VQE lacks convergence guarantees but employs smaller circuits, aligning better with this study's aim of advancing short-term quantum computers toward practical quantum-chemistry applications.

The paper is organized as follows. In Section “Methods”, we present two variants of the application of the DBBSC method to wave-function QC calculations, denoted as Strategy 1 and Strategy 2: (1) a basis-set correction *a posteriori* added to the solution of the quantum algorithm, integrating two contributions, namely a basis-set correlation density-functional correction and a basis-set Hartree–Fock (HF) correction, and (2) a self-consistent scheme integrating the DBBSC method to the quantum algorithm that dynamically modifies the one-electron density used in the basis-set correlation correction. Strategy 1 offers the possibility to correct any wave-

function QC energy calculation through a simple additive correction. Strategy 2 enables one to self-consistently access an improved electronic density, offering both improved energies and first-order molecular properties. Also, we introduce a new type of system-adapted basis sets (SABS) for Gaussian-type orbitals (GTOs), with sizes comparable to minimal basis sets. These new methodologies enable us to perform the first investigation of DBBSC corrections in the minimal basis-set regime.

In the Section “Results and discussion”, we provide relevant tests of our approach for different atomic and molecular systems using several families of basis sets. We carry out numerical simulations using graphics-processing-unit (GPU) accelerated QC sparse emulation on up to 32 qubits, exploring the applicability of this method to converge ground-state energies, dissociation curves, and dipole moments. We consistently observe significant improvements over typical quantum algorithm approaches, reaching an accuracy level that would have otherwise required hundreds of qubits. We thus expect the present approach to become a standard part of quantum-enhanced wave-function calculations, particularly the approach of Strategy 1, which can immediately provide large improvements to existing results with relatively simple efforts.

Finally, the Section “Conclusion” contains our conclusions. Additional details and results are provided in the Supplementary Information (SI).

## Results and discussion

We computed the ground-state energies, dissociation curves, and dipole moments for the N₂, H₂O, LiH, and H₂ molecules using both basis-set correction strategies. For all systems except H₂, the 1s molecular orbitals were frozen, i.e., we use the frozen-core approximation. Correspondingly, we use the frozen-core version of the DBBSC method as defined in ref. 44, in which, in particular, the contribution to the density coming from the core electrons is neglected in the basis-set correlation density functional. The error in each calculation is quantified as the deviation from the CBS limit, which itself is determined using a two-point extrapolation scheme based on the cc-pVQZ and cc-pV5Z basis sets⁵⁷. For the dipole–moment calculations we employed only Strategy 2, using an expectation value over the dipole–moment operator, avoiding a finite-difference approach prone to numerical errors⁵⁸. We remind the reader that these calculations are to be interpreted in the framework of perfect (noiseless) logical qubits, and therefore the discussed errors can be only attributed to the method used. The Appendix presents the quantum circuits that were used for all computations and provides the number of CNOT gates required to perform for the algorithm as well as the number of ADAPT-VQE iterations.

The classical computations, including basis-set corrections and the calculation of reference energies and dipole moments, were carried out using Quantum Package 2.0³⁹. In this software, the FCI energies are approximated by the energy from a CIPSI wave function to which a second-order perturbation theory (PT2) correction is added. Given the demonstrated nearly-FCI quality of this approximation for the systems we studied, we simplify our terminology by referring to this approach as simply FCI rather than CIPSI+PT2.

The SABS are labeled VXZ-Y where X is respectively D, T, Q, 5, or 6 for cc-pVDZ, cc-pVTZ, up to cc-pV6Z basis sets, and Y is the target size. We report in the SI the classically computed energies for the N₂, H₂O, LiH, and H₂ molecules.

## Ground-state energies

Ground-state energies of the H₂, LiH, H₂O, and N₂ molecules close to their respective equilibrium geometries are presented in Fig. 1 and Table 1. Details regarding the ADAPT-VQE iterations and the associated “ADAPT” values are provided in the SI, as well as additional tests on atoms and hydrogen chains.

In Fig. 1, a general trend is observed: the basis-set corrected ground-state energies with the small basis sets, i.e., STO-3G, 6-31G, or pcseg-0⁶⁰, align with values between the cc-pVDZ and cc-pVTZ basis-set levels, while requiring much fewer qubits than these latter basis sets. As we can see from Table 1, both strategies provide the same quantitative improvements. These
<figure><img src="image_5.png" /><figcaption>**Fig. 1** | Ground-state energy errors with respect to the extrapolated CBS limit. Outside the green box: ADAPT-VQE calculations. Inside the green box: FCI calculations from a classical computer. Dot markers correspond to energies without basis-set corrections. Cross markers correspond to energies corrected with Strategy 1</figcaption></figure>

[Eq. (6)]. On the x-axis, the labels *B*+ symbolize basis-set corrected values. The blue box corresponds to a range of 1.6 mHa around the extrapolated CBS limit. The errors for each ADAPT-VQE computation with respect to the FCI energy for a given basis set are reported in the SI.

results are consistent with the conclusion from ref. 49 where the non-self-consistent approximation was found sufficient for calculating energies.

For the best cases, the basis-set corrected ground-state energies have errors in the order of tens of mHa: 40 mHa for H₂O with 24 qubits, 60 mHa for N₂ with 16 qubits, and less than 10 mHa for LiH and H₂ with less than 10 qubits. The results stay consistent for the other systems available in the SI (H₄, H₆, H₈). For H₂ using the cc-pVDZ basis set, the basis-set corrected ground-state energy reaches chemical accuracy relative to the CBS limit while requiring 20 logical qubits. Similar convergence to the CBS limit is observed in the SI for He and Be. Without the basis-set correction, reaching a similar accuracy would have typically required more than a hundred qubits.

In the initial work on the DBBSC method⁴³, it was found that the basis-set corrected ground-state energies of atoms reach chemical accuracy with the cc-pVTZ basis set. Similarly, for the molecules studied here, if one had access to a few hundred logical qubits required for the cc-pVTZ basis set, achieving chemical accuracy would be theoretically feasible for all cases. To support this claim, we report in the SI classically computed basis-set corrected FCI ground-state energies using Dunning basis sets. The error relative to the CBS limit for the LiH molecule is 0.2 mHa with the cc-pVTZ basis set. For H₂O and N₂, we achieve errors of 3 and 1 mHa, respectively, with the cc-pVTZ basis set.

Let us now discuss the results obtained with our SABS for LiH, N₂, and H₂O. They were chosen to match the sizes of the small basis sets previously discussed. Specifically for H₂O, employing the V5Z-10 basis set (24 qubits) achieves results slightly superior to those obtained with the 6-31G basis set.

For H₂O, we observed a reduction in the HF basis-set correction by 40 mHa when moving from the 6-31G to the V5Z-10 basis set, by over 100 mHa when comparing the STO-3G to the V5Z-4 basis set for the LiH molecule, and more than a Hartree for the N₂ molecule moving from the STO-3G to the V5Z-6 and V5Z-11 basis sets. Based on these findings, we anticipate that further exploration of this strategy could lead to a basis-set correction scheme requiring only the correlation basis-set correction term. Of course, the interest of SABS is to be systematically improvable toward the parent Dunning basis set. It is important to note that this SABS systematic improvement is not limited to quantum algorithms such as ADAPT-VQE and is also present for the HF and FCI calculations. Since the SABS approach strongly reduces the computational cost while maintaining accuracy, it should offer further applications of the CIPSI approach in classical quantum chemistry. In practice, a cc-pVTZ-like quality is reached with a V5Z-11 basis set (32 qubits) for N₂, and with a V5Z-11 basis set (30 qubits) for H₂O. A V5Z-10 basis set (28 qubits) achieves a cc-pV5Z-like ground-state energy accuracy for LiH. Overall, the SABS always provides the best ADAPT energies with the self-consistent correlation basis-set correction (before the addition of the HF basis-set correction).

Overall, concerning the basis sets, it is important to point out a key anomaly, i.e., the remarkable performance of the minimal STO-3G basis set. This was expected as Pople already discussed the outperformance of this basis set in the seventies⁶¹. Indeed, Davidson and Feller detailed in their 1986 review⁶² the existence of error compensations and stated that “the smaller the basis the more ab initio calculations assume an empirical flavor”. In practice, among the available minimal basis-set possibilities, STO-3G
<table><caption>Table 1 | Ground-state energies (in Ha) for H₂O, N₂, LiH, H₂, and H₈ calculated by FCI, ADAPT-VQE (denoted as ADAPT), and basis-set corrected ADAPT-VQE according to Strategy 1, denoted as A + PBE + ΔHF, and to Strategy 2, denoted as SC(A + PBE) and SC(A + PBE) + ΔHF (i.e., without and with the HF basis-set correction, respectively)</caption><thead><tr><th></th><th>N<sub>qubits</sub></th><th>FCI</th><th>ADAPT</th><th>A + PBE + ΔHF</th><th>SC(A + PBE)</th><th>SC(A + PBE) + ΔHF</th></tr></thead><tbody><tr><td>H<sub>2</sub>O</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>STO-3G</td><td>12</td><td>-75.01250</td><td>-75.01250</td><td>-76.30232</td><td>-75.197880</td><td>-76.30191</td></tr><tr><td>pcseg-0</td><td>24</td><td>-75.90855</td><td>-75.90843<sup>a</sup></td><td>-76.33681</td><td>-76.03999</td><td>-76.33279</td></tr><tr><td>6-31G</td><td>24</td><td>-76.11995</td><td>-76.11989<sup>b</sup></td><td>-76.28035</td><td>-76.23717</td><td>-76.32025</td></tr><tr><td>V5Z-10</td><td>24</td><td>-76.12626</td><td>-76.12409</td><td>-76.32705</td><td>-76.27418<sup>c</sup></td><td>-76.32505</td></tr><tr><td>V5Z-11</td><td>30</td><td>-76.15902</td><td>-76.15165</td><td>-76.33704</td><td>-76.28622</td><td>-76.33570</td></tr><tr><td>cc-pVDZ</td><td>46</td><td>-76.24165</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pVTZ</td><td>114</td><td>-76.33250</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pVQZ</td><td>228</td><td>-76.35985</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pV5Z</td><td>400</td><td>-76.36877</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>CBS</td><td>-</td><td>-76.37812</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>N<sub>2</sub></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>STO-3G</td><td>16</td><td>-107.65251</td><td>-107.65251</td><td>-109.36630</td><td>-107.86974</td><td>-109.36661</td></tr><tr><td>V5Z-6</td><td>16</td><td>-108.88869</td><td>-108.88869</td><td>-109.34552</td><td>-109.09850</td><td>-109.34608</td></tr><tr><td>V5Z-11</td><td>32</td><td>-108.89413</td><td>-109.11566</td><td>-109.37278</td><td>-109.27385</td><td>-109.37281</td></tr><tr><td>cc-pVDZ</td><td>52</td><td>-109.27698</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pVTZ</td><td>116</td><td>-109.37527</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pVQZ</td><td>216</td><td>-109.40558</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pV5Z</td><td>360</td><td>-109.41505</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>CBS</td><td>-</td><td>-109.42498</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>LiH</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>STO-3G</td><td>10</td><td>-7.88218</td><td>-7.88218</td><td>-8.02160</td><td>-7.89590</td><td>-8.02119</td></tr><tr><td>pcseg-0</td><td>14</td><td>-7.98139</td><td>-7.98139</td><td>-8.0160</td><td>-7.99166</td><td>-8.01561</td></tr><tr><td>6-31G</td><td>20</td><td>-7.99800</td><td>-7.99800</td><td>-8.01668</td><td>-8.00806</td><td>-8.01611</td></tr><tr><td>V5Z-4</td><td>10</td><td>-7.99287</td><td>-7.99287</td><td>-8.01758</td><td>-8.00562</td><td>-8.01690</td></tr><tr><td>V5Z-7</td><td>16</td><td>-7.99793</td><td>-7.99793</td><td>-8.01710</td><td>-8.01109</td><td>-8.01643</td></tr><tr><td>V5Z-10</td><td>28</td><td>-8.01302</td><td>-8.01302</td><td>-8.02575</td><td>-8.02134</td><td>-8.02540</td></tr><tr><td>cc-pVDZ</td><td>26</td><td>-8.01438</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pVTZ</td><td>86</td><td>-8.02234</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pVQZ</td><td>190</td><td>-8.02386</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pV5Z</td><td>290</td><td>-8.02433</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>CBS</td><td>-</td><td>-8.02482</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>H<sub>2</sub></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>STO-3G</td><td>4</td><td>-1.13415</td><td>-1.13415</td><td>-1.17606</td><td>-1.15590</td><td>-1.17594</td></tr><tr><td>6-31G</td><td>8</td><td>-1.15003</td><td>-1.15003</td><td>-1.16911</td><td>-1.16196</td><td>-1.16913</td></tr><tr><td>cc-pVDZ</td><td>20</td><td>-1.16275</td><td>-1.16275</td><td>-1.17239</td><td>-1.16858</td><td>-1.17246</td></tr><tr><td>V5Z-8</td><td>24</td><td>-1.16613</td><td>-1.16613</td><td>-1.17315</td><td>-1.17170</td><td>-1.17320</td></tr><tr><td>cc-pVTZ</td><td>56</td><td>-1.17041</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pVQZ</td><td>120</td><td>-1.17182</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pV5Z</td><td>220</td><td>-1.17223</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>CBS</td><td>-</td><td>-1.17265</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>H<sub>8</sub></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>STO-3G</td><td>16</td><td>-4.24339</td><td>-4.24320</td><td>-4.45483</td><td>-4.32764</td><td>-4.45354</td></tr><tr><td>6-31G</td><td>32</td><td>-4.37032</td><td>-4.35752</td><td>-4.43488</td><td>-4.41275</td><td>-4.43502</td></tr><tr><td>cc-pVDZ</td><td>80</td><td>-4.42756</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pVTZ</td><td>222</td><td>-4.47121</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pVQZ</td><td>474</td><td>-4.47702</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>cc-pV5Z</td><td>864</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>CBS</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

Here, PBE refers to the PBE-based correlation basis-set correction ΔHF refers to the HF basis-set correction, and SC stands for "self-consistent". The frozen-core approximation has been used for H₂O, N₂, and LiH. The CBS limits are estimated by two-point extrapolations from cc-pVQZ and cc-pV5Z calculations.

<footnote n="500">⁵⁰⁰ iterations with a 14,668-determinant CIPSI initial state.</footnote>

<footnote n="6100">⁶¹⁰⁰ iterations with a 10,879-determinant CIPSI initial state.</footnote>

<footnote>The value is -76.27636 when using an 11,016-determinant CIPSI initial state.</footnote>
proved to be the most robust and cost-efficient choice for the basis-set corrections. Also, in the DBBSC schemes, the STO-3G basis set highly benefits from the HF basis-set correction, which reduces as the basis set and the number of qubits increases. When compared to SABS, for example, STO-3G always displays smaller self-consistent basis-set corrected ADAPT values and larger HF basis-set corrections. In practice, Pople also noted that STO-3G is especially good for energies around the equilibrium geometry, and larger basis sets are clearly required to describe accurately the notor-iously more difficult dissociation curves and dipole moments.

Finally, one last aspect of the analysis of such large basis-set simulations is related to the convergence of ADAPT-VQE. Indeed, in VQE-type computations³, there is no formal guarantee of convergence, and ADAPT-VQE belongs to such a heuristic family of methods. If, for systems like H₂, LiH, or N₂ (see SI, Figs. 3–14), convergence is obtained in a few dozen iterations, for more complex systems such as H₂O, the number of required iterations strongly increases when a double-zeta basis set is used. Clearly, Figure 2 in the SI shows that more than a thousand iterations are required to achieve full convergence starting from the HF initial state. This heuristic aspect means that no anticipation of the exact required number of iterations can be made. This may be manageable on a real quantum processor, but the use of classical emulation makes each iteration more costly in terms of time-to-solution than the previous one, limiting the overall convergence capabilities. It is possible to force convergence by replacing the HF starting point with a CIPSI one. Figure 2 in the SI shows that chemical accuracy can be easily reached with a well-converged CIPSI solution. However, ADAPT-VQE struggles to improve the solution, which is expected due to the enormous size of the parameter space. Furthermore, the CIPSI-based initial state already contains a significant amount of correlation. This represents a hard challenge for the ADAPT-VQE procedure, which now needs to pick the next ansatz operator to improve the existing quantum state. In any case, starting from an unconverged CIPSI wave function is a robust solution to strongly reduce the overall computational time. One example is given for H₂O and the V5Z-10 basis set which initially led to a not fully converged result. With a better CIPSI starting point, it is possible to recover a few milli-Hartrees (see footnotes in Table 1). It is also possible to change the operator pool, but the point here is that convergence becomes challenging when tackling complex electronic structures, and such computations would not be possible without GPU-accelerated emulation.

## Dissociation curves

We pursue the dissociation energies reported in Fig. 2. Clearly, the simple ADAPT-VQE/STO-3G level of theory appears quite far from an accurate description of the dissociation (compared to the calculations in the largest basis sets), which is achieved with FCI/triple-zeta (and beyond) levels. The non-corrected ADAPT-VQE values are represented with bold lines, whereas the corrected values are reported with dashed curves. These dissociation curves are extremely challenging as they involve several regimes of correlation going from weak correlation happening typically at short distances to strong correlation effects at large distances. First, we notice that for all the cases, the basis-set corrected values around the equilibrium are always substantially closer to the large cc-pV5Z reference values than the uncorrected ones and improve with the basis-set size and the use of SABS. For larger distances, the basis-set requirements appear even more stringent. In practice, the three dissociation curves manage to converge to high accuracy thanks to the use of SABS: we obtain nearly a triple-zeta quality for N₂ using the basis-set correction with the V5Z-6 basis set (16 qubits). A triple-zeta-like quality can be achieved using the larger V5Z-11 basis-set at the price of more qubits (32 qubits, see Table 2). In the same line, nearly a cc-pV5Z quality for H₂ using the basis-set correction with the V5Z-8 basis set, and nearly a cc-pV5Z quality for LiH using the basis-set correction with a V5Z-7 basis set. For H₂, the basis-set corrected cc-pVDZ curve matches the cc-pV5Z reference up to a distance of 2.5 Å, and a slight discrepancy appears at long distances. This is consistent with classical computations leading to convergence to the CBS limit only when the basis-set correction is applied to the cc-pVTZ basis set (114 qubits). It is possible to fix this issue more

affordably by using a V5Z-8 basis set requiring only 24 qubits. For LiH, it is also possible to fix the small discrepancies observed on the V5Z-7 basis-set dissociation curve using a larger V5Z-10 basis set (28 qubits) which reaches the CBS limit.

Finally, it is important to highlight the good performance of the DBBSC method for the prediction of the dissociation curve of the triple-bonded N₂ molecule. Indeed, such computation is well-documented in the literature and known as extremely difficult as it requires both a multi-reference treatment and various weak correlation effects going from short-range to charge polarization effects⁶³⁻⁶⁶. At the cost of a minimal basis set (i.e., 16 qubits), our N₂ DBBSC computations achieved nearly a cc-pVTZ quality, providing an accuracy that would have required around 100 logical qubits in the context of a brute-force simulation. To the best of our knowledge, these results are the most accurate using a quantum algorithm when compared to the recent results from the literature. Indeed, a group of researchers managed to predict this dissociation curve using a Local Unitary Cluster Jastrow (LUCJ) ansatz coupled to double-zeta basis sets (6-31G and cc-pVDZ)⁶⁷. The computation required the use of massive computational resources, namely hundreds of compute nodes of the Fugaku classical supercomputer coupled to a QPU. Alternatively, in the context of the fermionic quantum emulator, N₂ simulations used ADAPT-VQE coupled to the 6-31G and def2-SVP (56 qubits) basis sets to perform computation via approximate Matrix Product States (MPS)⁶⁸. In the present work, we achieve a better accuracy with far fewer qubits and, again, these results can be systematically improved using larger SABS. Indeed, access to a cc-pVTZ-like basis set accuracy is possible using the next (larger) SABS in terms of size, i.e., V5Z-11 (32 qubits).

## Dipole moments

We also explore the idea of extending the basis-set correction scheme in QC calculations of molecular properties such as the dipole moment. As pointed out by Halkier et al.,⁶⁹, the dipole moment also suffers from a slow basis-set convergence, which can be thought of as an indirect impact of the missing short-range correlation effects in a finite basis set. Therefore, it is relevant to apply the basis-set correction to the dipole moment, as already shown in ref. 49. We report in Table 2 the dipole moments of H₂O and LiH. We calculate the basis-set corrected dipole moment as the expectation value over the dipole-moment operator over the last ADAPT-VQE wave function coming out of the self-consistent basis-set correction method. We also add a *posteriori* the HF basis-set correction to the dipole moment, calculated as the difference between the HF dipole moment in the aug-cc-pV5Z basis set and the HF dipole moment in the considered basis set. From Table 2, we see that the correlation basis-set correction is not sufficient to converge the dipole moments. However, the addition of the HF basis-set correction significantly improves the basis-set convergence. For both molecules, strong improvements are observed with respect to the ADAPT-VQE/STO-3G level. Further improvements are observed as the size of the basis set increases. Finally, we remark the ADAPT-VQE slow convergence in terms of iterations is also the source of small errors in the dipole moments.

## Conclusions

We demonstrated the applicability of the DBBSC method to QC algorithms for quantum chemistry. Using ADAPT-VQE, these approaches were shown to be able to systematically improve minimal basis-set results by predicting ground-state energies that are intermediate between double-zeta and triple-zeta FCI qualities. Overall, the presented self-consistent basis-set corrected ground-state energies are very close to their non-self-consistent counterparts. As the latter a *posteriori* basis-set correction approach can be very easily applied to any wave-function QC calculation performed on real quantum hardware, it provides an affordable improvement strategy for current QC chemistry computations. The self-consistent basis-set correction scheme is still useful since it permits the calculation of properties such as dipole moments thanks to the availability of improved QC densities.

An additional reduction of the required basis-set size is provided by our SABS approach. Besides being fast, as SABS can be generated within
**Fig. 2** | Dissociation curves of H₂, LiH, and N₂ molecules. Blue and orange lines correspond to Dunning basis sets using the FCI method. The notation VXZ stands for the basis set cc-pVXZ. Green and purple plain lines correspond to ADAPT-VQE calculations without basis-set corrections. Green and purple dash-dotted lines correspond to basis-set corrected dissociation curves. Raw data are available in the SI.

<figure><img src="image_6.png" /></figure>

seconds, such a “black-box” pivoted Cholesky strategy for the on-the-fly generation of basis sets has been shown to be competitive and often superior to available standard choices. It offers systematic computational savings for large basis sets, reducing significantly the qubit requirements. Thus the SABS’ usefulness is not limited to QC algorithms as they also offer systematically improvable solutions for performing reference classical

computations. Indeed, since SABS are user-defined truncated versions of the Dunning basis sets that can match a qubit budget, they also offer reduced cost access to improved accuracy for any quantum-chemistry method.

Our methodology allows us to compute chemically meaningful energies and properties on systems that would have required far more than 100 logical qubits. For example, the computation of the H₂ total energy at the
<table><caption>**Table 2 | Dipole moments (in atomic units) of LiH and H₂O calculated by FCI and self-consistently basis-set corrected ADAPT-VQE, denoted as SC(A + PBE) and SC(A + PBE) + ΔHF (without and with and the HF basis-set correction, respectively)**</caption><thead><tr><th>LiH</th><th>N<sub>qubits</sub></th><th>FCI</th><th>SC(A + PBE)</th><th>SC(A + PBE) + ΔHF</th></tr></thead><tbody><tr><td>STO-3G</td><td>10</td><td>-1.81835</td><td>-1.86299</td><td>-2.31321</td></tr><tr><td>pcseg-0</td><td>14</td><td>-2.33313</td><td>-2.37289</td><td>-2.25650</td></tr><tr><td>6-31G</td><td>20</td><td>-2.16646</td><td>-2.20768</td><td>-2.23674</td></tr><tr><td>V5Z-4</td><td>10</td><td>-2.37818</td><td>-2.44145</td><td>-2.22886</td></tr><tr><td>V5Z-7</td><td>16</td><td>-2.23095</td><td>-2.27789</td><td>-2.25618</td></tr><tr><td>V5Z-10</td><td>28</td><td>-2.24997</td><td>-2.27458</td><td>-2.31438</td></tr><tr><td>cc-pVDZ</td><td>26</td><td>-2.25566</td><td>-</td><td>-</td></tr><tr><td>cc-pVTZ</td><td>86</td><td>-2.29998</td><td>-</td><td>-</td></tr><tr><td>cc-pVQZ</td><td>190</td><td>-2.30361</td><td>-</td><td>-</td></tr><tr><td>cc-pV5Z</td><td>290</td><td>-2.30647</td><td>-</td><td>-</td></tr><tr><th>H<sub>2</sub>O</th><th>N<sub>qubits</sub></th><th>FCI</th><th>SC(A + PBE)</th><th>SC(A + PBE) + ΔHF</th></tr><tr><td>STO-3G</td><td>12</td><td>-0.63584</td><td>-0.67084</td><td>-0.77162</td></tr><tr><td>pcseg-0</td><td>24</td><td>-0.95822</td><td>-0.99450</td><td>-0.77065</td></tr><tr><td>6-31G</td><td>24</td><td>-0.99020</td><td>-1.01898</td><td>-0.76342</td></tr><tr><td>V5Z-10</td><td>24</td><td>-0.99305</td><td>-1.01887</td><td>-0.77634</td></tr><tr><td>V5Z-11</td><td>30</td><td>-0.99185</td><td>-1.02170</td><td>-0.77737</td></tr><tr><td>cc-pVDZ</td><td>46</td><td>-0.76073</td><td>-</td><td>-</td></tr><tr><td>cc-pVTZ</td><td>114</td><td>-0.75013</td><td>-</td><td>-</td></tr><tr><td>cc-pVQZ</td><td>228</td><td>-0.74994</td><td>-</td><td>-</td></tr><tr><td>cc-pV5Z</td><td>400</td><td>-0.74241</td><td>-</td><td>-</td></tr></tbody></table>

Here, ΔHF corresponds to the HF basis-set correction to the dipole moment, calculated as the difference between the HF dipole moment in the aug-cc-pV5Z basis set and the HF dipole moment in the considered basis set. More data are available in the SI.

FCI/cc-pV5Z level, which would have required more than 220 logical qubits, can be achieved here with only 24 qubits using our basis-set correction scheme and our SABS technique. Overall, we were able to converge four systems to the FCI/CBS limit, including He, Be, H₂, and LiH. We were also able to provide accurate dissociation curves for H₂, LiH, and N₂. Computations on H₂ and LiH required the use of only a single GPU. Since most quantum-chemistry studies are presently out of reach of quantum computers, this research opens the path to more affordable quantitative quantum-chemistry simulations of small molecules using QC algorithms. In particular, the present a posteriori basis-set corrections can be easily added to any type of STO-3G VQE fermionic computations on real hardware⁵,⁷⁰,⁷¹, allowing one to improve significantly their accuracy at very little computational cost. Since adaptive simulations on real hardware are making some progress²⁷ while the hardware itself improves, basis-set corrected simulations should be progressively possible on future quantum computers providing a route to FCI/CBS quality computations. This strategy is particularly suited for resources demanding computations that converge slowly and require large basis sets associated with large qubit counts. Finally, the DBBSC method is not limited to ground-state computations and can be extended to excited states via a linear-response formalism⁵³. This will be the subject of future QC research.

In this context, besides the state-of-the-art ADAPT-VQE hybrid quantum-classical algorithm, it would be interesting to revisit accurate fixed wave-function ansätze⁷², such as UCCSDT⁷³ and others, to analyze their convergence when coupled to the DBBSC method. In practice, the presented basis-set correction framework is not restricted to a given QC wave-function ansatz and can leverage any future QC algorithmic improvements. Our present DBBSC/SABS methodology still requires the use of qubits either through quantum hardware or through the use of a classical quantum

emulator. The qubit count is, therefore, our main limitation. The computations of this paper used up to 32 logical qubits and represent a proof-of-concept study of what one can presently do with quantum emulation to prepare the advent of fault-tolerant quantum computing⁷⁴, see Appendix for detailed resource estimations. However, state-vector simulations have limitations due to memory. Indeed, if they are theoretically possible up to 40-50 qubits on very large exascale supercomputers, they start to become relatively unpractical for quantum-chemistry simulations when reaching 36 qubits due to the high computational resources and time-to-solution requirements. To explore further the chemical electronic space with quantum algorithms, we are currently upgrading our Hyperion-1 framework to increase our emulated qubit counts by going beyond the state-vector formalism thanks to various elements from the CUDA-Q SDK⁷⁵,⁷⁶ developed by NVIDIA. We are also presently testing various implementations of the DBBSC algorithms on available quantum hardware. To conclude, by reducing the number of qubits required to reach the CBS limit, we expect to tackle predictive real-world quantum chemistry applications with strategies applicable to both NISQ and FTQC algorithms.

## Methods

The overall methodological procedure is illustrated in Fig. 3. The procedure starts with the definition of the system, and a standard basis set or our SABS. After defining the second-quantized Hamiltonian, the quantum state is prepared on a quantum processing unit (QPU) or a GPU-accelerated quantum emulator. The basis-set corrections are calculated on the classical CPU using one of the two DBBSC variants labeled as Strategy 1 and Strategy 2.

### Density-based basis-set correction method

In the infinite-dimensional (antisymmetric) N-electron Hilbert space, $\mathcal{H} = \bigwedge^{N} L^2(\mathbb{R}^3 \times \{\uparrow, \downarrow\}, \mathcal{C})$, we consider an atomic or molecular system with Hamiltonian

$$ \hat{H} = \hat{T} + \hat{W}_{ee} + \hat{V}_{ne}, \quad (1) $$

where $\hat{T}$ is the kinetic-energy operator, $\hat{W}_{ee}$ is the Coulomb electron-electron operator, and $\hat{V}_{ne}$ is the nuclei-electron potential operator. The exact ground-state energy is defined as

$$ E_0 = \min_{\Psi \in W} \langle \Psi, \hat{H}\Psi \rangle, \quad (2) $$

where $W = \{\Psi \in \bigwedge^N H^1(\mathbb{R}^3 \times \{\uparrow, \downarrow\}, \mathcal{C}) | \langle \Psi, \Psi \rangle = 1\}$ is the space of admissible wave functions, $H^1$ is the first-order Sobolev space, and $\langle \cdot, \cdot \rangle$ designates the standard inner product of $\mathcal{H}$. In quantum-chemistry calculations, we normally introduce a one-electron basis set $B \subset H^1(\mathbb{R}^3 \times \{\uparrow, \downarrow\}, \mathcal{C})$ and we work in the finite-dimensional N-electron Hilbert space generated by this basis set, i.e., $\mathcal{H}^B = \bigwedge^N \text{span}(B)$. The FCI ground-state energy is then defined as

$$ E_{\text{FCI}}^{B} = \min_{\Psi \in W^{B}} \langle \Psi, \hat{H}\Psi \rangle, \quad (3) $$

where $W^B = \{\Psi \in \mathcal{H}^B | \langle \Psi, \Psi \rangle = 1\}$ is the space of $B$-representable wave functions. In the CBS limit, the FCI ground-state energy tends to the exact ground-state energy, i.e., $E_{\text{FCI}}^{B} \text{CBS} = E_0$, but the convergence is infamously slow due to short-range electron correlation⁵⁷.

In the DBBSC method¹⁴,⁴⁹, for the given basis set $B$, we introduce the following approximation to the ground-state energy

$$ E_0^B = \min_{\Psi \in W^B} \langle \Psi, \hat{H}\Psi \rangle + \bar{E}^B[n_\Psi], \quad (4) $$

where $\bar{E}^B[n_\Psi]$ is a basis-set correction density functional evaluated at the one-electron density of $\Psi$. This basis-set correction density functional must vanish in the CBS limit so that $E_0^B$ properly converges to the exact ground-
<figure><img src="image_7.png" /><figcaption>**Fig. 3 | Architecture of the hybrid quantum-classical scheme introducing the DBBSC method to wave-function QC calculations.** The QPU workload can be either performed by quantum hardware or replaced by a GPU-accelerated quantum emulator. The input are the definition of the system and the basis set (standard basis sets or our SABS can be used). In Strategy 1, the QPU/GPU computes the expectation</figcaption></figure>

value of the standard Hamiltonian $\hat{H}$ and the non-self-consistent correlation basis-set correction and the HF basis-set correction are shifted to the CPU of the classical computer. In Strategy 2, the model includes a self-consistent basis-set correction potential which is iteratively optimized between QPU/GPU and CPU. The insets depict the effect of the basis-set corrections for the two strategies.

state energy, i.e., $E_0^B \to \text{CBS} = E_0$, but it must be such that it accelerates the convergence to the CBS limit. In ref. 44, based on the Perdew-Burke-Ernzerhof (PBE) correlation density functional$^{77}$, such a basis-set correction density functional was constructed in a semilocal form

$$ \bar{E}^B[n] = \int_{\mathbb{R}^3} \bar{e}^B(n(\mathbf{r}), \nabla n(\mathbf{r})) d\mathbf{r}, \quad (5) $$

where $n(\mathbf{r})$ and $\nabla n(\mathbf{r})$ are the density and density gradient at point $\mathbf{r}$, and the function $\bar{e}^B(n, \nabla n)$, including how the dependence on the basis set $\mathcal{B}$ is included, can be found in ref. 44. We will now discuss how the DBBSC method can be adapted to wave-function QC calculations.

**Strategy 1: a posteriori basis-set correction.** In Strategy 1, we use a non-self-consistent approximation to Eq. (4). The idea is to calculate the FCI energy $E_{\text{FCI}}^B$, or a good approximation to it, on a quantum computer, and then add a posteriori the basis-set correlation correction of Eq. (5) and a HF basis-set correction, both calculated on a classical computer. This leads to the following approximation to the ground-state energy

$$ E_1^B = E_{\text{FCI}}^B + \bar{E}^B[n_{\text{HF}}^B] + \Delta E_{\text{HF}}^B. \quad (6) $$

As in previous works$^{43,44}$, the basis-set correlation correction is evaluated at the HF density in the basis set $\mathcal{B}$, leading to a calculation of $\bar{E}_{\text{HF}}^B[n_{\text{HF}}^B]$ with a marginal computational cost with respect to the cost of calculating $E_{\text{FCI}}^B$. The PBE-based basis-set correlation correction only corrects for basis-set incompleteness errors due to short-range electron correlation. However, when using small basis sets, a significant part of the basis-set error also comes from the fact the HF part of the energy is not converged to the CBS limit. Similarly to Refs. 54,55, we thus add a HF basis-set correction. In the

present work, we choose it simply as the difference between the HF energy in the CBS limit, that we estimate as the value $E_{\text{HF}}^{5Z}$ obtained with the cc-pV5Z basis set$^{79}$, and the HF energy $E_{\text{HF}}^B$ in the basis set $\mathcal{B}$

$$ \Delta E_{\text{HF}}^B = E_{\text{HF}}^{5Z} - E_{\text{HF}}^B. \quad (7) $$

We note that it is also possible to avoid performing a HF calculation with the cc-pV5Z basis set for estimating the CBS limit of the HF energy by using a complementary auxiliary basis set, as in refs. 54,55, but we have not found it necessary for the present work.

Let us emphasize again that the PBE-based basis-set correction mentioned above only takes care of the correlation part and therefore does not correct the basis-set error of the HF energy. Hence, there is no double counting between the two basis-set corrections, as shown in Refs. 54,55.

**Strategy 2: self-consistent basis-set correction.** In Strategy 2, we use the self-consistent version of the DBBSC method in Eq. (4), using the basis-set correlation correction of Eq. (5) in the self-consistent part of the calculation, and we only add a posteriori the fixed HF basis-set correction. This leads to the following approximation to the ground-state energy that can be applied to any variational quantum ansatz

$$ E_2^B = \min_{\Psi \in W^B} \left( \langle \Psi, \hat{H} \Psi \rangle + \bar{E}^B[n_{\Psi}] \right) + \Delta E_{\text{HF}}^B. \quad (8) $$

The minimization in Eq. (8) leads to the following self-consistent Schrödinger equation$^{49}$

$$ \hat{P}^B \hat{H}^B [n_{\Psi^B}] \Psi^B = \varepsilon^B \Psi^B, \quad (9) $$
<table><caption>**Table 3 | Walltime (in minutes) required to grow the wave-function ansatz of size N<sub>adapt</sub> for a given molecular system using Hyperion-1 state-vector emulator on N<sub>gpus</sub> NVIDIA GPUs**</caption><thead><tr><th>Molecule/basis set</th><th>N<sub>adapt</sub></th><th>N<sub>qubits</sub></th><th>N<sub>gpus</sub></th><th>A100 walltime [min]</th><th>H100 walltime [min]</th></tr></thead><tbody><tr><td>H<sub>2</sub>O/6-31G</td><td>500</td><td>24</td><td>1</td><td>644</td><td>503</td></tr><tr><td>H<sub>12</sub>/STO-3G</td><td>500</td><td>24</td><td>1</td><td>174</td><td>134</td></tr><tr><td>H<sub>14</sub>/STO-3G</td><td>300</td><td>28</td><td>8</td><td>283</td><td>184</td></tr><tr><td>H<sub>16</sub>/STO-3G</td><td>100</td><td>32</td><td>128</td><td>450</td><td>147</td></tr></tbody></table>

The results have been obtained on A100 (80 GB) and H100 (80 GB) GPUs using CUDA Toolkit 12.0 and NVIDIA HPC SDK 23.3.

where $\hat{P}^B$ is the projector on the N-electron Hilbert space generated by the basis set $B$, i.e., $\mathcal{H}^B$, and $\hat{\bar{H}}^B[n]$ is the effective Hamiltonian

$$ \hat{\bar{H}}^B[n] = \hat{H} + \hat{\bar{V}}^B[n]. \qquad (10) $$

Here, $\hat{\bar{V}}^B[n]$ is the basis-set correction one-electron potential operator

$$ \hat{\bar{V}}^B[n] = \int_{\mathbb{R}^3} \bar{v}^B[n](\mathbf{r}) \hat{n}(\mathbf{r}) \, d\mathbf{r}, \qquad (11) $$

where $\bar{v}^B[n](\mathbf{r}) = \delta\bar{E}^B[n]/\delta n(\mathbf{r})$ is the derivative of the basis-set correlation correction with respect to the density, and $\hat{n}(\mathbf{r})$ is the density operator.

The idea is now to solve iteratively Eq. (9) on a quantum computer. The potential $\bar{v}^B[n](\mathbf{r})$, and therefore the Hamiltonian $\hat{\bar{H}}^B[n]$, is iteratively updated with the density $n_i^B$ of the wave-function ansatz solution $\Psi_i^B$ of the $i^{th}$ iteration. The convergence criterion is reached when the difference between the last two iterated energy eigenvalues $\mathcal{E}_i^B$ is less than 0.1 mHa. In practice, this is fulfilled after two iterations.

Finally, the self-consistent basis-set corrected energy $E_2^B$ can be directly computed from the last energy eigenvalue $\mathcal{E}^B$ and density $n^B = n_{\psi}^B$ coming out from the quantum solver as follows

$$ E_2^B = \mathcal{E}^B + \bar{E}^B[n^B] - \int_{\mathbb{R}^3} \bar{v}^B[n^B](\mathbf{r}) n^B(\mathbf{r}) \, d\mathbf{r} + \Delta E_{\text{HF}}^B, \qquad (12) $$

where the HF basis-set correction $\Delta E_{\text{HF}}^B$ is identical to the one used in Strategy 1.

In Strategy 2, the density of the QC wave-function ansatz $\Psi$ needs to be calculated. For this, we calculate the one-particle density matrix

$$ n_{pq,\sigma} = \langle \Psi, \hat{a}_{p,\sigma}^{\dagger} \hat{a}_{q,\sigma} \Psi \rangle, \qquad (13) $$

where $\hat{a}_{p,\sigma}^{\dagger}$ and $\hat{a}_{q,\sigma}$ are the creation and annihilation operators, respectively, for the $p$th and $q$th orbitals and with spin $\sigma \in \{\uparrow, \downarrow\}$. This is achieved with a Jordan-Wigner mapping of the operator $\hat{a}_{p,\sigma}^{\dagger} \hat{a}_{q,\sigma}$ and the state $\Psi$ prepared using the parameterized ansatz. The density matrix is then used to classically update the Hamiltonian in Eq. (10) and to calculate dipole moments.

**GPU-accelerated QPU emulation and shift of the basis-set correction to classical resources.** The DBBSC method does not increase the qubit count as the basis-set correction is performed classically using a density functional and a HF correction. In practice, such computations present no advantage to be performed at the QC level since they would consume a large number of qubits without accuracy benefit over its classical counterpart. Therefore a hybrid QC wave-function/classical DBBSC approach is highly preferable. In this work, VQE/DBBSC computations are performed classically on GPU-accelerated computing nodes, the GPUs replacing the quantum processing unit (QPU) through the use of a classical quantum emulator. The basis-set correction contributions can be then shifted to the unused CPUs to maximize the computational efficiency while the computationally challenging VQE

algorithm is fully offloaded on GPUs for optimal time-to-solution performances. Of course, the proposed hybrid schemes are also fully tractable on a real quantum computer, as the basis-set correction contributions can be entirely shifted to classical resources harnessing GPUs/CPUs to preserve the qubit count while the core wave-function ansatz is maintained on the QPU.

## System-adapted basis-set generation

Within the DBBSC method, the natural choice for basis sets is the Dunning correlation-consistent basis-set family that offers a path to a reliable convergence to the CBS limit, as demonstrated in initial DBBSC classical quantum-chemistry studies$^{43–55}$. For QC calculations, such basis sets are unfortunately usually out of reach of quantum hardware/emulators as they are associated with very large qubit counts. In the present section, we address the task of generating atomic-orbital (AO) basis sets under a basis-size budget, with controllable accuracy$^{79}$. We propose a mathematical formulation of this goal in terms of a constraint optimization problem, as follows.

Given a molecule composed of $N_{\text{atm}}$ atoms, with fixed nuclear coordinates $\{\mathbf{r}_a\}_{1 \le a \le N_{\text{atm}}}$, let

$$ \mathcal{B} = \{\chi_{\mu}\}_{1 \le \mu \le N_{\text{bas}}} = \bigcup_{a=1}^{N_{\text{atm}}} \{\chi_{\mu}^{a}(\mathbf{r} - \mathbf{r}_{a})\}_{1 \le \mu \le N_{\text{bas}}^{a}} $$

denote a spatial basis set of atom-centered GTOs. It is assumed that a one-electron density, denoted by $n_0^B$, is available after a converged ground-state energy minimization procedure (e.g., HF or CI). Our main purpose is to a *posteriori* extract subsets of the given “large” basis set $\mathcal{B}$, denoted by $\mathcal{B}_I = \{\chi_{\mu}\}_{\mu \in I} \subseteq \mathcal{B}$ for any index subset $I \subseteq \{1, \dots, N_{\text{bas}}\}$, achieving (i) a target size, and (ii) minimal accuracy loss on the given density $n_0^B$ used as a reference. In the present computing scheme, since the DBBSC method uses the cc-pV5Z basis set for the HF basis-set correction, we choose $\mathcal{B} = \text{cc-pV5Z}$. In practice, this choice is motivated by the fact that such a quintuple-zeta basis set fairly approaches chemical accuracy compared to quasi-exact all-electron numeric atom-centered orbital computations$^{80}$. Hence, given a target basis size $M$ smaller than $N_{\text{bas}}$, we seek the optimal index subset, denoted by $I_M$, that minimizes the best approximation error of the density $n_0^B$ over the set of $\mathcal{B}_I$-representable one-electron densities, for any $I \subseteq \{1, \dots, N_{\text{bas}}\}$ with $|I| = M$, where $|\cdot|$ denotes the cardinal of a set, or, in other words, we solve an optimization problem under constraints:

$$ I_M := \arg \min_{I \subseteq \{1, \dots, N_{\text{bas}}\}} \min_{n^B} \|n_0^B - n^B_I\|, \qquad (14) \\ |I| = M $$

where $\| \cdot \|$ is a given norm for functions over $\mathbb{R}^3$. In the present work, we use the norm induced by the Coulomb inner product $\langle u, v \rangle_c := \int_{\mathbb{R}^3} \int_{\mathbb{R}^3} d\mathbf{r} d\mathbf{r}' u(\mathbf{r})v(\mathbf{r}')|\mathbf{r} - \mathbf{r}'|^{-1}$. Let us emphasize that Eq. (14) requires knowledge of the density $n_0^B$, which is assumed to be pre-computed on a classical computer. Such quantity is available in the framework of the DBBSC scheme since HF computations are systematically performed at the cc-pV5Z level to estimate the CBS limit of the HF energy.
The problem in Eq. (14) can be solved as follows. Our approach is to identify a greedy procedure for discarding elements of the full AO-product set that spans the space containing the reference density $n_0^B$, which admits the expansion

$$n_0^B = \sum_{\mu, \nu=1}^{N_{\text{bas}}} D_{\mu\nu} \chi_{\mu} \chi_{\nu},$$

where $D = CC^T$ is the density matrix and $C$ is the $N_{\text{bas}} \times N_{\text{occ}}$ matrix of the coefficients in the AO basis of the $N_{\text{occ}}$ occupied molecular orbitals. We plan to achieve this by eliminating linear dependencies present in the AO-product set$^{81}$. The pivoted Cholesky decomposition (PCD) of a matrix$^{82}$ is an algebraic tool for eliminating linear dependencies occurring between matrix rows (resp. columns), which may be interpreted as an iterative greedy procedure for discarding elements that do not contribute to the full row (resp. column) space, up to an orthogonal projection error tolerance. PCD has been previously applied to the auxiliary basis-set generation for density fitting$^{81,83,84}$. In the present work, we employ PCD within a new scheme, named system-adapted basis-set (SABS) generation, for solving the problem in Eq. (14).

Let us formulate our scheme in detail. Prior to the AO-product selection and in order to ensure orbital symmetry of the resulting SABS, we pre-process the initial basis $\mathcal{B}$ and first contract all angular components (e.g., all three p-type components $p_\alpha, p_\beta, p_\gamma$) of GTOs. To this end, we consider the partition $\{B_i\}_{i=1}^{N_{\text{orb}}}$ of $\{1, ..., N_{\text{bas}}\}$, each $B_i$ being an index block containing all angular components of a single GTO in $\mathcal{B}$, and the $N_{\text{orb}} \times N_{\text{bas}}$ contraction matrix $P$, defined for any $1 \le i \le N_{\text{orb}}$ as $P_{ij} = 1$ if $j \in B_i$ and zero otherwise. Next, we define the four-index tensor $T$ with entries

$$T_{pqrs} = \sum_{\mu, \nu, \kappa, \lambda=1}^{N_{\text{bas}}} P_{\mu\kappa} P_{\nu\lambda} D_{\mu\nu} \langle \chi_{\mu} \chi_{\nu}, \chi_{\kappa} \chi_{\lambda} \rangle_c D_{\kappa\lambda} P_{\rho\kappa} P_{\rho\lambda},$$

and fold pairwise its first two and its last two dimensions, in order to form the $N_{\text{prod}} \times N_{\text{prod}}$ Gram matrix of weighted AO products, denoted by $G$, with $N_{\text{prod}} = (N_{\text{orb}})^2$. As a last pre-processing step, we discard rows and columns of $G$ corresponding to products made of components not centered on the same atom and denote the resulting submatrix $A$. Now, PCD is applied to $A$, using the machine epsilon as a tolerance threshold, yielding an index set sorted in pivot-descending order. Given a target $M$, we define the selected AO-product index set, denoted by $S_M^{\text{Chol}}$, as the $M$-first pivot indices. Lastly, we recover the underlying AO index set

$$J_M = \bigcup_{(p_1, p_2) \in S_M^{\text{Chol}}} \{p_1, p_2\},$$

and post-process it in order to ensure that the same GTO-parameter set is assigned to all atoms of the same chemical type. For this purpose, the new basis set associated to a chemical type is the union of parameters of those GTOs in $B_{I_M}$ that are centered on any atom of that type. This yields a solution $I_M \supseteq J_M$ to our problem in Eq. (14) and the resulting SABS is $B_{I_M}$.

Note that our generation scheme directly fixes the size $M$ of the selected products, i.e., $|S_M^{\text{Chol}}| = M$. The actual size of the AO basis set, equal to $|I_M|$, is only implicitly controlled during our procedure. In practice, as numerical results show, $|I_M|$ is very close to $M$ for s- and p-type basis sets, while it remains the same order as $M$ for higher angular-momentum orbital types. Overall, the SABS generation approach is extremely fast and offers access to compact basis sets, specifically adapted to a given system and user-defined qubit budget. Examples of SABS generation can be found in the SI.

We note that the adaptation and generation of basis sets to the molecular geometry has been the subject of several publications exploring other strategies. Among them, we can cite some recent works$^{85,86}$, but also some other research in the context of quantum computing focusing on the need to

<table><caption>**Table 4 | CNOT-gate counts for qubit and qubit-excitation-based (QEB) operator pools**</caption><thead><tr><th>H<sub>2</sub></th><th>N<sub>qubits</sub></th><th>N<sub>CNOT</sub><br>(if QEB)</th><th>N<sub>CNOT</sub><br>(if qubit)</th><th>N<sub>op</sub></th><th>N<sub>adapt iter</sub></th></tr></thead><tbody><tr><td>STO-3G</td><td>4</td><td>13</td><td>6</td><td>1</td><td>1</td></tr><tr><td>6-31G</td><td>8</td><td>71</td><td>34</td><td>7</td><td>7</td></tr><tr><td>cc-pVDZ</td><td>20</td><td>233</td><td>110</td><td>21</td><td>21</td></tr><tr><td>V5Z-8</td><td>24</td><td>395</td><td>186</td><td>35</td><td>35</td></tr><tr><th>H<sub>4</sub></th><th>N<sub>qubits</sub></th><th>N<sub>CNOT</sub><br>(if QEB)</th><th>N<sub>CNOT</sub><br>(if qubit)</th><th>N<sub>op</sub></th><th>N<sub>adapt iter</sub></th></tr><tr><td>STO-3G</td><td>8</td><td>207</td><td>98</td><td>19</td><td>19</td></tr><tr><th>H<sub>6</sub></th><th>N<sub>qubits</sub></th><th>N<sub>CNOT</sub><br>(if QEB)</th><th>N<sub>CNOT</sub><br>(if qubit)</th><th>N<sub>op</sub></th><th>N<sub>adapt iter</sub></th></tr><tr><td>STO-3G</td><td>12</td><td>2437</td><td>1134</td><td>199</td><td>199</td></tr><tr><th>H<sub>8</sub></th><th>N<sub>qubits</sub></th><th>N<sub>CNOT</sub><br>(if QEB)</th><th>N<sub>CNOT</sub><br>(if qubit)</th><th>N<sub>op</sub></th><th>N<sub>adapt iter</sub></th></tr><tr><td>STO-3G</td><td>16</td><td>12677</td><td>5870</td><td>999</td><td>999</td></tr><tr><td>6-31G</td><td>32</td><td>8855</td><td>4090</td><td>685</td><td>685</td></tr><tr><th>He</th><th>N<sub>qubits</sub></th><th>N<sub>CNOT</sub><br>(if QEB)</th><th>N<sub>CNOT</sub><br>(if qubit)</th><th>N<sub>op</sub></th><th>N<sub>adapt iter</sub></th></tr><tr><td>pc-seg0</td><td>4</td><td>19</td><td>10</td><td>3</td><td>3</td></tr><tr><td>6-31G</td><td>4</td><td>19</td><td>10</td><td>3</td><td>3</td></tr><tr><td>cc-pVDZ</td><td>10</td><td>58</td><td>28</td><td>6</td><td>6</td></tr><tr><th>Be</th><th>N<sub>qubits</sub></th><th>N<sub>CNOT</sub><br>(if QEB)</th><th>N<sub>CNOT</sub><br>(if qubit)</th><th>N<sub>op</sub></th><th>N<sub>adapt iter</sub></th></tr><tr><td>STO-3G</td><td>8</td><td>39</td><td>18</td><td>3</td><td>3</td></tr><tr><td>pc-seg0</td><td>10</td><td>58</td><td>28</td><td>6</td><td>6</td></tr><tr><td>6-31G</td><td>16</td><td>207</td><td>98</td><td>19</td><td>19</td></tr><tr><td>cc-pVDZ</td><td>26</td><td>292</td><td>136</td><td>26</td><td>26</td></tr><tr><th>LiH</th><th>N<sub>qubits</sub></th><th>N<sub>CNOT</sub><br>(if QEB)</th><th>N<sub>CNOT</sub><br>(if qubit)</th><th>N<sub>op</sub></th><th>N<sub>adapt iter</sub></th></tr><tr><td>STO-3G</td><td>10</td><td>90</td><td>44</td><td>10</td><td>10</td></tr><tr><td>pc-seg0</td><td>14</td><td>258</td><td>124</td><td>26</td><td>26</td></tr><tr><td>6-31G</td><td>20</td><td>511</td><td>242</td><td>47</td><td>47</td></tr><tr><td>VQZ-4</td><td>10</td><td>90</td><td>44</td><td>10</td><td>10</td></tr><tr><td>V5Z-4</td><td>10</td><td>90</td><td>44</td><td>10</td><td>10</td></tr><tr><td>V5Z-7</td><td>16</td><td>291</td><td>138</td><td>27</td><td>27</td></tr><tr><td>V5Z-10</td><td>28</td><td>1083</td><td>506</td><td>91</td><td>91</td></tr><tr><th>H<sub>2</sub>O</th><th>N<sub>qubits</sub></th><th>N<sub>CNOT</sub><br>(if QEB)</th><th>N<sub>CNOT</sub><br>(if qubit)</th><th>N<sub>op</sub></th><th>N<sub>adapt iter</sub></th></tr><tr><td>STO-3G</td><td>12</td><td>729</td><td>342</td><td>63</td><td>63</td></tr><tr><td>6-31G</td><td>24</td><td>12657</td><td>5862</td><td>999</td><td>999</td></tr><tr><td>V5Z-11</td><td>30</td><td>12647</td><td>5858</td><td>999</td><td>999</td></tr><tr><th>N<sub>2</sub></th><th>N<sub>qubits</sub></th><th>N<sub>CNOT</sub><br>(if QEB)</th><th>N<sub>CNOT</sub><br>(if qubit)</th><th>N<sub>op</sub></th><th>N<sub>adapt iter</sub></th></tr><tr><td>STO-3G</td><td>16</td><td>4868</td><td>2256</td><td>386</td><td>386</td></tr><tr><td>V5Z-6</td><td>16</td><td>9714</td><td>4492</td><td>758</td><td>758</td></tr><tr><td>V5Z-11</td><td>32</td><td>11718</td><td>5420</td><td>916</td><td>916</td></tr></tbody></table>

$N_{\text{op}}$ is the number of operators in the final wave-function ansatz, and $N_{\text{adapt iter}}$ is the number of ADAPT-VQE iterations achieved and at which we collect the values in the Table.

The numbers of CNOT gates, $N_{\text{CNOT}}$, are evaluated as $N_{\text{single}} + 13N_{\text{double}}$ for the QEB pool, and as $N_{\text{single}} + 6N_{\text{double}}$ for the qubit pool, where $N_{\text{simple}}$ is the number of single-qubit operators and $N_{\text{double}}$ is the number of two-qubit operators.

limit the qubit requirements through basis-set reoptimization such as refs. 87,88. Alternatively, Kottman and Aspuru$^{89}$ proposed a basis-set-free approach through an adaptive representation using pair-natural orbitals which was tested up to 22 qubits.
## Computational details

In the present study, we perform ADAPT-VQE computations²⁰ using the Qubit-Excitation-Based pool of operators, which is considered a standard²². Additional details about the ADAPT-VQE methodology can be found in the SI. The communication of the Hamiltonian from the CPU to the QPU/ GPU is done by using standard FCIDUMP files to communicate the one- and two-electron integrals to the QPU/GPU software in order to construct the fermion operators. TREXIO files⁹⁰ are also useful to communicate a wider range of relevant information.

ADAPT-VQE computations were performed using the Hyperion-1 GPU-accelerated state-vector sparse emulator⁹¹ up to 32 qubits. Hyperion-1 uses classical computing systems and is grounded on an efficient multi-GPU ensemble of fast custom sparse linear-algebra libraries accelerating Hyperion-1's exact/noiseless simulations. In this paper, computations were performed on NVIDIA DGX A100 nodes (8x 80GB A100 GPUs per node) and NVIDIA DGX H100 nodes (4x & 8x 80 GB H100 GPUs per node). QC calculations being strongly memory-dependent, a single GPU can carry out a 20-qubit ADAPT-VQE simulation depending on the nature (i.e., Hamiltonian sparsity) of the system whereas a single node (8 GPUs) can handle up to 28 qubits. Multi-node computations are required beyond such a qubit count. Further details about Hyperion-1 and its full capabilities will be given in a forthcoming publication. Convergence for all ADAPT-VQE computations were set to $10^{-6}$ Ha. Most computations started from a HF initial state. For selected ones (indicated in the text and Tables), we started the ADAPT-VQE procedure from a rough configuration-interaction perturbatively-selected-iteratively (CIPSI)⁹² initial state (converged to only $10^{-2}$ Ha) to save some computational time within Hyperion-1. This reflects a commonly adopted strategy where a multi-determinant initial state is employed instead of a single HF determinant to increase the ground-state support in the initial state. The *Quantum State Preparation* of such classically-derived CIPSI wave functions has been studied in the context of VQE²⁵ and QPE²⁶. Also, we report in Table 3, the walltime required for several ground-state energy calculations. As one can see, the walltime is not only a function of the number of qubits, sparsity is also important. For example, the Hamiltonian of the $H_{12}$ molecule is sparser than the one of the water molecule resulting in a faster convergence. Besides the increased computing power of H100 GPUs leading to improved time-to-solution, our results also highlight the importance of fast node-to-node interconnects when performing large-scale quantum emulation. Indeed, the benefit of H100 over A100 is striking for the largest 32 qubits simulations on 16 nodes where an improvement of factor 3 was observed on the DGX H100 systems. Such speedup is, therefore, also partially related to higher node-to-node bandwidth observed on DGX H100 versus A100 systems.

## Appendix: CNOT counts

CNOT-gate counts for qubit and qubit-excitation-based (QEB) operator pools are presented in Table 4.

### Data availability

Data generated during the study is available upon request from the authors (E-mail: jean-philip.piquemal@sorbonne-universite.fr).

### Code availability

The code used during the study is available upon request from the authors (E-mail: jean-philip.piquemal@sorbonne-universite.fr).

Received: 3 July 2024; Accepted: 30 October 2024;

Published online: 18 November 2024

## References

1. Kitaev, A. Y. Quantum measurements and the abelian stabilizer problem. *Electron. Colloq. Comput. Complex TR96* https://wwwsemanticscholar.org/paper/Quantum-measurements-and-the-Abelian-Stabilizer-Kitaev/e218de049ab533e7d54f336cadac942effddf139* (1995).

2. Aspuru-Guzik, A., Dutoi, A. D., Love, P. J. & Head-Gordon, M. Simulated quantum computation of molecular energies. *Science* **309**, 1704–1707 (2005).

3. Nielsen, M. A. & Chuang, I. L. *Quantum Computation and Quantum Information* (Cambridge University Press, 2010).

4. Whitfield, J. D., Biamonte, J. & Aspuru-Guzik, A. Simulation of electronic structure hamiltonians using quantum computers. *Mol. Phys.* **109**, 735–750 (2011).

5. Peruzzo, A. et al. A variational eigenvalue solver on a photonic quantum processor. *Nat. Commun.* **5**, 1–7 (2014).

6. Tilly, J. et al. The variational quantum eigensolver: a review of methods and best practices. *Phys. Rep.* **986**, 1–128 (2022).

7. Helgaker, T., Klopper, W. & Tew, D. P. Quantitative quantum chemistry. *Mol. Phys.* **106**, 2107–2143 (2008).

8. Bartlett, R. J. & Musiał, M. Coupled-cluster theory in quantum chemistry. *Rev. Mod. Phys.* **79**, 291–352 (2007).

9. Teale, A. M. et al. DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science. *Phys. Chem. Chem. Phys.* **24**, 28700–28781 (2022).

10. Chan, G. K.-L. & Sharma, S. The density matrix renormalization group in quantum chemistry. *Annu. Rev. Phys. Chem.* **62**, 465–481 (2011).

11. Cao, Y. et al. Quantum chemistry in the age of quantum computing. *Chem. Rev.* **119**, 10856–10915 (2019).

12. Bauer, B., Bravyi, S., Motta, M. & Chan, G. K.-L. Quantum algorithms for quantum chemistry and quantum materials science. *Chem. Rev.* **120**, 12685–12717 (2020).

13. Haidar, M., Rančić, M. J., Ayral, T., Maday, Y. & Piquemal, J.-P. Open source variational quantum eigensolver extension of the quantum learning machine for quantum chemistry. *Wiley Interdiscip. Rev.* **13**, e1664 (2023).

14. Jones, G. M. & Jacobsen, H.-A. Distributed quantum computing for chemical applications. Preprint at arXiv https://doi.org/10.48550/arXiv.2408.05351 (2024).

15. McArdle, S., Endo, S., Aspuru-Guzik, A., Benjamin, S. C. & Yuan, X. Quantum computational chemistry. *Rev. Mod. Phys.* **92**, 015003 (2020).

16. McClean, J. R., Romero, J., Babbush, R. & Aspuru-Guzik, A. The theory of variational hybrid quantum-classical algorithms. *N. J. Phys.* **18**, 023023 (2016).

17. Romero, J. et al. Strategies for quantum computing molecular energies using the unitary coupled cluster ansatz. *Quantum Sci. Technol.* **4**, 014008 (2018).

18. Haidar, M., Rancic, M. J., Maday, Y. & Piquemal, J.-P. Extension of the trotterized unitary coupled cluster to triple excitations. *J. Phys. Chem. A* **127**, 3543–3550 (2023).

19. Haidar, M., Adjoua, O., Baddredine, S., Peruzzo, A. & Piquemal, J.-P. Non-iterative disentangled unitary coupled-cluster based on lie-algebraic structure. Preprint at arXiv https://arxiv.org/abs/2408.14289 (2024).

20. Grimsley, H. R., Economou, S. E., Barnes, E. & Mayhall, N. J. An adaptive variational algorithm for exact molecular simulations on a quantum computer. *Nat. Commun.* **10**, 3007 (2019).

21. Tang, H. L. et al. qubit-adapt-vqe: An adaptive algorithm for constructing hardware-efficient ansätze on a quantum processor. *PRX Quantum* **2**, 020310 (2021).

22. Yordanov, Y. S., Armaos, V., Barnes, C. H. & Arvidsson-Shukur, D. R. Qubit-excitation-based adaptive variational quantum eigensolver. *Commun. Phys.* **4**, 228 (2021).

23. Anastasiou, P. G., Chen, Y., Mayhall, N. J., Barnes, E. & Economou, S. E. TETRIS-ADAPT-VQE: an adaptive algorithm that yields shallower, denser circuit ansätze. *Phys. Rev. Research* **6**, 013254 (2024).

24. Shkolnikov, V. O., Mayhall, N. J., Economou, S. E. & Barnes, E. Avoiding symmetry roadblocks and minimizing the measurement overhead of adaptive variational quantum eigensolvers. *Quantum* **7**, 1040 (2023).
25. Feniou, C. et al. Overlap-ADAPT-VQE: practical quantum chemistry on quantum computers via overlap-guided compact ansätze. *Commun. Phys.* **6** (2023).

26. Feniou, C. et al. Sparse quantum state preparation for strongly correlated systems. *J. Phys. Chem. Lett.* **15**, 3197–3205 (2024).

27. Feniou, C. et al. Greedy gradient-free adaptive variational quantum algorithms on a noisy intermediate scale quantum computer. Preprint at arXiv https://doi.org/10.48550/arXiv.2306.17159 (2023).

28. Burton, H. G. A. Accurate and gate-efficient quantum ansätze for electronic states without adaptive optimization. *Phys. Rev. Res.* **6**, 023300 (2024).

29. Hättig, C., Klopper, W., Köhn, A. & Tew, D. P. Explicitly correlated electrons in molecules. *Chem. Rev.* **112**, 4–74 (2012).

30. Ten-no, S. & Noga, J. Explicitly correlated electronic structure theory from r12/f12 ansätze. *WIREs Comput. Mol. Sci.* **2**, 114–125 (2012).

31. Boys, S. F., Handy, N. C. & Linnett, J. W. A condition to remove the indeterminacy in interelectronic correlation functions. *Proc. R. Soc. Lond. A Math. Phys. Sci.* **309**, 209–220 (1969).

32. Ten-no, S. A feasible transcorrelated method for treating electronic cusps using a frozen Gaussian geminal. *Chem. Phys. Lett.* **330**, 169–174 (2000).

33. Cohen, A. J. et al. Similarity transformation of the electronic Schrödinger equation via Jastrow factorization. *J. Chem. Phys.* **151**, 061101 (2019).

34. Baiardi, A., Lesiuk, M. & Reiher, M. Explicitly correlated electronic structure calculations with transcorrelated matrix product operators. *J. Chem. Theory Comput.* **18**, 4203–4217 (2022).

35. Ammar, A., Scemama, A. & Giner, E. Extension of selected configuration interaction for transcorrelated methods. *J. Chem. Phys.* **157**, 134107 (2022).

36. Liao, K. et al. Density matrix renormalization group for transcorrelated hamiltonians: ground and excited states in molecules. *J. Chem. Theory Comput.* **19**, 1734–1743 (2023).

37. Kato, T. On the eigenfunctions of many-particle systems in quantum mechanics. *Commun. Pure Appl. Math.* **10**, 151–177 (1957).

38. Motta, M. et al. Quantum simulation of electronic structure with a transcorrelated hamiltonian: improved accuracy with a smaller footprint on the quantum computer. *Phys. Chem. Chem. Phys.* **22**, 24270–24281 (2020).

39. Kumar, A. et al. Quantum simulation of molecular electronic states with a transcorrelated hamiltonian: Higher accuracy with fewer qubits. *J. Chem. Theory Comput.* **18**, 5312–5324 (2022).

40. Dobrautz, W. et al. Toward real chemical accuracy on current quantum hardware through the transcorrelated method. *J. Chem. Theory Comput.* **20**, 4146–4160 (2024).

41. Schleich, P., Kottmann, J. S. & Aspuru-Guzik, A. Improving the accuracy of the variational quantum eigensolver for molecular systems by the explicitly-correlated perturbative [2]r12-correction. *Phys. Chem. Chem. Phys.* **24**, 13550–13564 (2022).

42. Sokolov, I. O., Dobrautz, W., Luo, H., Alavi, A. & Tavernelli, I. Orders of magnitude increased accuracy for quantum many-body problems on quantum computers via an exact transcorrelated method. *Phys. Rev. Res.* **5**, 023174 (2023).

43. Giner, E. et al. Curing basis-set convergence of wave-function theory using density-functional theory: a systematically improvable approach. *J. Chem. Phys.* **149**, 194301 (2018).

44. Loos, P.-F., Pradines, B., Scemama, A., Toulouse, J. & Giner, E. A density-based basis-set correction for wave function theory. *J. Phys. Chem. Lett.* **10**, 2931–2937 (2019).

45. Giner, E., Scemama, A., Toulouse, J. & Loos, P.-F. Chemically accurate excitation energies with small basis sets. *J. Chem. Phys.* **151**, 144118 (2019).

46. Loos, P.-F., Pradines, B., Scemama, A., Giner, E. & Toulouse, J. Density-based basis-set incompleteness correction for GW methods. *J. Chem. Theory Comput.* **16**, 1018–1028 (2020).

47. Giner, E., Scemama, A., Loos, P.-F. & Toulouse, J. A basis-set error correction based on density-functional theory for strongly correlated molecular systems. *J. Chem. Phys.* **152**, 174104 (2020).

48. Yao, Y., Giner, E., Li, J., Toulouse, J. & Umrigar, C. J. Almost exact energies for the gaussian-2 set with the semistochastic heat-bath configuration interaction method. *J. Chem. Phys.* **153**, 124117 (2020).

49. Giner, E., Traore, D., Pradines, B. & Toulouse, J. Self-consistent density-based basis-set correction: How much do we lower total energies and improve dipole moments? *J. Chem. Phys.* **155**, 044109 (2021).

50. Yao, Y., Giner, E., Anderson, T. A., Toulouse, J. & Umrigar, C. J. Accurate energies of transition metal atoms, ions, and monoxides using selected configuration interaction and density-based basis-set corrections. *J. Chem. Phys.* **155**, 204104 (2021).

51. Traore, D., Giner, E. & Toulouse, J. Basis-set correction based on density-functional theory: rigorous framework for a one-dimensional model. *J. Chem. Phys.* **156**, 044113 (2022).

52. Traore, D., Toulouse, J. & Giner, E. Basis-set correction for coupled-cluster estimation of dipole moments. *J. Chem. Phys.* **156**, 174101 (2022).

53. Traore, D., Giner, E. & Toulouse, J. Basis-set correction based on density-functional theory: Linear-response formalism for excited-state energies. *J. Chem. Phys.* **158**, 234107 (2023).

54. Mester, D. & Kállay, M. Basis set limit of CCSD(T) energies: explicit correlation versus density-based basis-set correction. *J. Chem. Theory Comput.* **19**, 8210–8222 (2023).

55. Heßelmann, A. et al. A density-fitting implementation of the density-based basis-set correction method. *J. Comput. Chem.* **45**, 1247–1253 (2024).

56. Dunning, T. H. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. *J. Chem. Phys.* **90**, 1007 (1989).

57. Helgaker, T., Klopper, W., Koch, H. & Noga, J. Basis-set convergence of correlated calculations on water. *J. Chem. Phys.* **106**, 9639–9646 (1997).

58. Kassal, I. & Aspuru-Guzik, A. Quantum algorithm for molecular properties and geometry optimization. *J. Chem. Phys.* **131**, 224102 (2009).

59. Garniron, Y. et al. Quantum package 2.0: an open-source determinant-driven suite of programs. *J. Chem. Theory Comput.* **15**, 3591–3609 (2019).

60. Jensen, F. Unifying general and segmented contracted basis sets. Segmented polarization consistent basis sets. *J. Chem. Theory Comput.* **10**, 1074–1085 (2014).

61. Pople, J. A. *Modern Theoretical Chemistry* (Plenum, New York, 1976).

62. Davidson, E. R. & Feller, D. Basis set selection for molecular calculations. *Chem. Rev.* **86**, 681–696 (1986).

63. Tubman, N. et al. Postponing the orthogonality catastrophe: efficient state preparation for electronic structure simulations on quantum devices. *Bull. Am. Phys. Soc.* **64** https://meetings.aps.org/Meeting/MAR19/Session/E42.4 (2019).

64. Almlöf, J., Deleeuw, B. J., Taylor, P. R., Bauschlicher Jr, C. W. & Siegbahn, P. The dissociation energy of N2. *Int. J. Quantum Chem.* **36**, 345–354 (1989).

65. Peterson, K. A. & Dunning Jr, T. H. Intrinsic errors in several ab initio methods: the dissociation energy of N2. *J. Phys. Chem.* **99**, 3898–3901 (1995).

66. Malrieu, J.-P., Guihéry, N., Calzado, C. J. & Angeli, C. Bond electron pair: Its relevance and analysis from the quantum chemistry point of view. *J. Comput. Chem.* **28**, 35–50 (2007).

67. Robledo-Moreno, J. et al. Chemistry beyond exact solutions on a quantum-centric supercomputer. Preprint at arXiv https://doi.org/10.48550/arXiv.2405.05068 (2024).

68. Provazza, J. et al. Fast emulation of fermionic circuits with matrix product states. *Journal of Chemical Theory and Computation* (2024).
69. Halkier, A., Klopper, W., Helgaker, T. & Jørgensen, P. Basis-set convergence of the molecular electric dipole moment. *J. Chem. Phys.* **111**, 4424–4430 (1999).

70. Kandala, A. et al. Hardware-efficient variational quantum eigensolver for small molecules and quantum magnets. *Nature* **549**, 242–246 (2017).

71. Sennane, W., Piquemal, J.-P. & Rančić, M. J. Calculating the ground-state energy of benzene under spatial deformations with noisy quantum computing. *Phys. Rev. A* **107**, 012416 (2023).

72. Haidar, M., Rančić, M. J., Ayral, T., Maday, Y. & Piquemal, J.-P. Open source variational quantum eigensolver extension of the quantum learning machine for quantum chemistry. *WIREs Comput. Mol. Sci.* **13**, e1664 (2023).

73. Haidar, M., Rančić, M. J., Maday, Y. & Piquemal, J.-P. Extension of the trotterized unitary coupled cluster to triple excitations. *J. Phys. Chem. A* **127**, 3543–3550 (2023).

74. Acharya, R. et al. Quantum error correction below the surface code threshold. Preprint at arXiv https://arxiv.org/abs/2408.13687 (2024).

75. Bayraktar, H. et al. cuquantum SDK: A high-performance library for accelerating quantum science. In *2023 IEEE International Conference on Quantum Computing and Engineering (QCE)*, vol. 1, 1050–1061 (IEEE, 2023).

76. Kim, J.-S. et al. CUDA quantum: The platform for integrated quantum-classical computing. In *2023 60th ACM/IEEE Design Automation Conference (DAC)*, 1–4 (IEEE, 2023).

77. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. *Phys. Rev. Lett.* **77**, 3865 (1996).

78. Halkier, A., Helgaker, T., Jørgensen, P., Klopper, W. & Olsen, J. Basis-set convergence of the energy in molecular hartree-fock calculations. *Chem. Phys. Lett.* **302**, 437 (1999).

79. Lygatsika, I.-M. *Méthodes Numériques Pour Les Discrétisations Gaussiennes Des Problèmes En Structure électronique*. Phd thesis, Sorbonne Université, Paris, France (2023).

80. Jensen, S. R. et al. The elephant in the room of density functional theory calculations. *J. Phys. Chem. Lett.* **8**, 1449–1457 (2017).

81. Lehtola, S. Curing basis set overcompleteness with pivoted Cholesky decompositions. *J. Chem. Phys.* **151**, 241102 (2019).

82. Lucas, C. LAPACK-style codes for level 2 and 3 pivoted Cholesky factorizations. https://www.netlib.org/lapack/lawnspdf/lawn161.pdf (2004).

83. Aquilante, F. et al. Cholesky decomposition techniques in electronic structure theory. *Linear-Scaling Techn. Comput. Chem. Phys.* **301–343** (2011).

84. Pedersen, T. B., Lehtola, S., Fdez. Galván, I. & Lindh, R. The versatility of the Cholesky decomposition in electronic structure theory. *Wiley Interdiscip. Rev.* **14**, e1692 (2024).

85. Schutt, O. & VandeVondele, J. Machine learning adaptive basis sets for efficient large scale density functional theory simulation. *J. Chem. Theory Comput.* **14**, 4168–4175 (2018).

86. Mao, Y. et al. Approaching the basis set limit for DFT calculations using an environment-adapted minimal basis with perturbation theory: formulation, proof of concept, and a pilot implementation. *J. Chem. Phys.* **145**, 044109 (2016).

87. Wang, W. & Whitfield, J. D. Basis set generation and optimization in the nisq era with quiqbox.jl. *J. Chem. Theory Comput.* **19**, 8032–8052 (2023).

88. Kwon, H.-Y., Curtin, G. M., Morrow, Z., Kelley, C. & Jakubikova, E. Adaptive basis sets for practical quantum computing. *Int. J. Quantum Chem.* **123**, e27123 (2023).

89. Kottmann, J. S., Schleich, P., Tamayo-Mendoza, T. & Aspuru-Guzik, A. Reducing qubit requirements while maintaining numerical precision for the variational quantum eigensolver: A basis-set-free approach. *J. Phys. Chem. Lett.* **12**, 663–673 (2021).

90. Posenitskiy, E. et al. TREXIO: a file format and library for quantum chemistry. *J. Chem. Phys.* **158**, 174801 (2023).

91. Adjoua, O. & Feniou, C. *Sorbonne Université, CNRS and Qubit Pharmaceuticals* (2023).

92. Huron, B., Malrieu, J. & Rancurel, P. Iterative perturbation calculations of ground and excited state energies from multiconfigurational zeroth-order wavefunctions. *J. Chem. Phys.* **58**, 5745–5759 (1973).

## Acknowledgements

We thank the PEPR EPIQ—Quantum Software (ANR-22-PETQ-0007, J.-P.P.) for funding D. Traore's initial postdoctoral position. Support of the HQI initiative is also acknowledged (J.-P.P. and Y.M.). I.-M. Lygatsika's PhD position has been funded by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (grant No. 810367, project EMC2, J.-P.P. and Y.M.). Computations have been performed at IDRIS (Jean Zay) on GENCI Grants: No. A0150712052 (J.-P.P.) and grant GC010815453 (Grand Challenge H100 Jean Zay, J.-P.P.), at Scaleway (Jero) and on the SuperPOD reference cluster EOS at NVIDIA. The authors thank D. Ruau and C. Hundt (NVIDIA) for continuous support, as well as the NVIDIA Quantum team for technical discussions.

## Author contributions

D.T., O.A., I.M.L., E.P., and K.H. performed the computations. O.A., D.T., C.F., and E.P. wrote the quantum emulation code. D.T., O.A., I.M.L., Y.M., E.G., A.P., J.T., E.G., and J.P.P. analyzed the data. D.T., O.A., C.F., I.M.L., Y.M., E.P., J.T., E.G., and J.P.P. contributed to the methodology. D.T., O.A., I.M.L., A.P., J.T., E.G., and J.P.P. wrote the paper with inputs from all other co-authors. J.P.P. designed and supervised the research.

## Competing interests

J.P.P. is a shareholder and co-founder of Qubit Pharmaceuticals. The remaining authors declare no other competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42004-024-01348-3.

**Correspondence** and requests for materials should be addressed to Jean-Philip Piquemal.

**Peer review information** *Communications Chemistry* thanks Mandeep Kaur Saggi, Anthony Gandon, Weishi Wang, and the other, anonymous, reviewer for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024