Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2301.01457 · Journal of Chemical Theory and Computation

Bootstrap Embedding on a Quantum Computer

Yuan Liu, Oinam Romesh Meitei, Zachary E. Chin, Arkopal Dutt, Max Tao, Isaac L. Chuang, Troy Van Voorhis

Published Mar 31, 2023• 20 citations
arXivHTMLPDF

Abstract

We extend molecular bootstrap embedding to make it appropriate for implementation on a quantum computer. This enables solution of the electronic structure problem of a large molecule as an optimization problem for a composite Lagrangian governing fragments of the total system, in such a way that fragment solutions can harness the capabilities of quantum computers. By employing state-of-art quantum subroutines including the quantum SWAP test and quantum amplitude amplification, we show how a quadratic speedup can be obtained over the classical algorithm, in principle. Utilization of quantum computation also allows the algorithm to match─at little additional computational cost─full density matrices at fragment boundaries, instead of being limited to 1-RDMs. Current quantum computers are small, but quantum bootstrap embedding provides a potentially generalizable strategy for harnessing such small machines through quantum fragment matching.

DOI
10.1021/acs.jctc.3c00012

Related papers

Bootstrap Embedding for Molecules

Hong‐Zhou Ye, Nathan D. Ricke, Henry K. Tran +1 more

Fragment embedding is one way to circumvent the high computational scaling of accurate electron correlation methods. The challenge of applying fragment embedding to molecular systems primarily lies in the strong entanglement and correlation that prevent accurate fragmentation across chemical bonds. Recently, Schmidt decomposition has been shown effective for embedding fragments that are strongly coupled to a bath in several model systems. In this work, we extend a recently developed quantum embedding scheme, bootstrap embedding (BE), to molecular systems. The resulting method utilizes the matching conditions naturally arising from using overlapping fragments to optimize the embedding. Numerical simulation suggests that the accuracy of the embedding improves rapidly with fragment size for small molecules, whereas larger fragments that include orbitals from different atoms may be needed for larger molecules. BE scales linearly with system size (apart from an integral transform) and hence can potentially be useful for large-scale calculations.

Journal of Chemical Theory and Computation•Jul 25, 2019•61 citations

Variational Quantum Eigensolver-Based Quantum Bootstrap Embedding for Molecules

Derek Peng

Simulating strongly correlated molecular systems on near-term quantum hardware remains challenging because current hardware offers limited quantum volume and moderate-fidelity qubits. One way to address this challenge is through bootstrap embedding (BE). Bootstrap embedding partitions molecules into smaller fragments that are then embedded in the ``bath'' of other fragments iteratively. Bootstrap embedding is appealing for quantum simulation because fragmentation reduces the qubit requirements of each fragment problem. In this work, we develop a quantum bootstrap embedding (QBE) workflow that uses variational quantum eigensolver (VQE) fragment solvers and examines the algorithmic choices that determine the success of the overall VQE-QBE procedure. To improve efficiency, we introduce FastAdaptVQE, a sparse-matrix-accelerated form of the adaptive variational quantum eigensolver (ADAPT-VQE) that replaces symbolic commutator evaluation with direct statevector calculations, and MatrixFreeAdaptVQE, a matrix-free extension that removes the sparse-matrix memory bottleneck that arises for larger fragments. We also explore modifying the ADAPT-VQE operator-selection step by replacing the purely greedy choice with a lookahead strategy. Benchmarks on H$_4$, F$_2$, and selected geometries of H$_6$ reach chemical accuracy, within 1 kcal/mol of BE results obtained with a full configuration interaction (FCI) solver. For $n$-butane, density-matching residuals converged, but energy errors remained outside chemical accuracy. These results show that combining QBE with VQE can accurately compute molecular energies, while density convergence alone does not ensure energy accuracy. This work lays the foundation for extending these energy calculations to larger molecular systems and quantum materials on near-term quantum hardware.

physics.chem-phCross-listedarXiv:2606.17095•Jun 16, 2026•0 citations

Bootstrap Embedding For Large Molecular Systems

Hong‐Zhou Ye, Henry K. Tran, Troy Van Voorhis

Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [ J. Chem. Theory Comput. 2019, 15, 4497–4506; J. Phys. Chem. Lett. 2019, 10, 6368–6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O ( N 3 ) for the integral transform and O ( N ) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.

Journal of Chemical Theory and Computation•Jun 26, 2020•40 citations

Multiscale Embedding for Quantum Computing

Leah P. Weisburn, Minsik Cho, Moritz Bensberg +3 more

We present a novel multiscale embedding scheme that links conventional QM/MM embedding and bootstrap embedding (BE) to allow simulations of large chemical systems on limited quantum devices. We also propose a mixed-basis BE scheme that facilitates BE calculations on extended systems using classical computers with limited memory resources. Benchmark data suggest the combination of these two strategies as a robust path in attaining the correlation energies of large realistic systems, combining the proven accuracy of BE with chemical and biological systems of interest in a lower computational cost method. Due to the flexible tunability of the resource requirements and systematic fragment construction, future developments in the realization of quantum computers naturally offer improved accuracy for multiscale BE calculations.

arXiv:2409.06813•Journal of Chemical Theory and Computation•Apr 15, 2025•19 citations

Bootstrap embedding: An internally consistent fragment-based method

Matthew Welborn, Takashi Tsuchimochi, Troy Van Voorhis

Strong correlation poses a difficult problem for electronic structure theory, with computational cost scaling quickly with system size. Fragment embedding is an attractive approach to this problem. By dividing a large complicated system into smaller manageable fragments "embedded" in an approximate description of the rest of the system, we can hope to ameliorate the steep cost of correlated calculations. While appealing, these methods often converge slowly with fragment size because of small errors at the boundary between fragment and bath. We describe a new electronic embedding method, dubbed "Bootstrap Embedding," a self-consistent wavefunction-in-wavefunction embedding theory that uses overlapping fragments to improve the description of fragment edges. We apply this method to the one dimensional Hubbard model and a translationally asymmetric variant, and find that it performs very well for energies and populations. We find Bootstrap Embedding converges rapidly with embedded fragment size, overcoming the surface-area-to-volume-ratio error typical of many embedding methods. We anticipate that this method may lead to a low-scaling, high accuracy treatment of electron correlation in large molecular systems.

The Journal of Chemical Physics•Aug 15, 2016•96 citations

Molecular Quantum Computations on a Protein

Akhil Shajan, Danil Kaliakin, Fangchun Liang +7 more

This work presents the implementation of a fragment-based, quantum-centric supercomputing workflow for computing molecular electronic structure using quantum hardware. The workflow is applied to predict the relative energies of two conformers of the 300-atom Trp-cage miniprotein. The methodology employs wave function-based embedding (EWF) as the underlying fragmentation framework, in which all atoms in the system are explicitly included in the CI treatment. CI calculations for individual fragments are performed using either sample-based quantum diagonalization (SQD) for challenging fragments or full configuration interaction (FCI) for trivial fragments. To assess the accuracy of SQD for fragment CI calculations, EWF-(FCI,SQD) results are compared against EWF-MP2 and EWF-CCSD benchmarks. Overall, the results demonstrate that large-scale electronic configuration interaction (CI) simulations of protein systems containing hundreds or even thousands of atoms can be realized through the combined use of quantum and classical computing resources.

quant-pharXiv:2512.17130•Dec 22, 2025•0 citations
We extend molecular bootstrap embedding to make it appropriate for implementation on a quantum computer. This enables solution of the electronic structure problem of a large molecule as an optimization problem for a composite Lagrangian governing fragments of the total system, in such a way that fragment solutions can harness the capabilities of quantum computers. By employing state-of-art quantum subroutines including the quantum SWAP test and quantum amplitude amplification, we show how a quadratic speedup can be obtained over the classical algorithm, in principle. Utilization of quantum computation also allows the algorithm to match─at little additional computational cost─full density matrices at fragment boundaries, instead of being limited to 1-RDMs. Current quantum computers are small, but quantum bootstrap embedding provides a potentially generalizable strategy for harnessing such small machines through quantum fragment matching.
Fragment embedding is one way to circumvent the high computational scaling of accurate electron correlation methods. The challenge of applying fragment embedding to molecular systems primarily lies in the strong entanglement and correlation that prevent accurate fragmentation across chemical bonds. Recently, Schmidt decomposition has been shown effective for embedding fragments that are strongly coupled to a bath in several model systems. In this work, we extend a recently developed quantum embedding scheme, bootstrap embedding (BE), to molecular systems. The resulting method utilizes the matching conditions naturally arising from using overlapping fragments to optimize the embedding. Numerical simulation suggests that the accuracy of the embedding improves rapidly with fragment size for small molecules, whereas larger fragments that include orbitals from different atoms may be needed for larger molecules. BE scales linearly with system size (apart from an integral transform) and hence can potentially be useful for large-scale calculations.
Simulating strongly correlated molecular systems on near-term quantum hardware remains challenging because current hardware offers limited quantum volume and moderate-fidelity qubits. One way to address this challenge is through bootstrap embedding (BE). Bootstrap embedding partitions molecules into smaller fragments that are then embedded in the ``bath'' of other fragments iteratively. Bootstrap embedding is appealing for quantum simulation because fragmentation reduces the qubit requirements of each fragment problem. In this work, we develop a quantum bootstrap embedding (QBE) workflow that uses variational quantum eigensolver (VQE) fragment solvers and examines the algorithmic choices that determine the success of the overall VQE-QBE procedure. To improve efficiency, we introduce FastAdaptVQE, a sparse-matrix-accelerated form of the adaptive variational quantum eigensolver (ADAPT-VQE) that replaces symbolic commutator evaluation with direct statevector calculations, and MatrixFreeAdaptVQE, a matrix-free extension that removes the sparse-matrix memory bottleneck that arises for larger fragments. We also explore modifying the ADAPT-VQE operator-selection step by replacing the purely greedy choice with a lookahead strategy. Benchmarks on H4_44​, F2_22​, and selected geometries of H6_66​ reach chemical accuracy, within 1 kcal/mol of BE results obtained with a full configuration interaction (FCI) solver. For nnn-butane, density-matching residuals converged, but energy errors remained outside chemical accuracy. These results show that combining QBE with VQE can accurately compute molecular energies, while density convergence alone does not ensure energy accuracy. This work lays the foundation for extending these energy calculations to larger molecular systems and quantum materials on near-term quantum hardware.
Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [ J. Chem. Theory Comput. 2019, 15, 4497–4506; J. Phys. Chem. Lett. 2019, 10, 6368–6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O ( N 3 ) for the integral transform and O ( N ) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.
We present a novel multiscale embedding scheme that links conventional QM/MM embedding and bootstrap embedding (BE) to allow simulations of large chemical systems on limited quantum devices. We also propose a mixed-basis BE scheme that facilitates BE calculations on extended systems using classical computers with limited memory resources. Benchmark data suggest the combination of these two strategies as a robust path in attaining the correlation energies of large realistic systems, combining the proven accuracy of BE with chemical and biological systems of interest in a lower computational cost method. Due to the flexible tunability of the resource requirements and systematic fragment construction, future developments in the realization of quantum computers naturally offer improved accuracy for multiscale BE calculations.
Strong correlation poses a difficult problem for electronic structure theory, with computational cost scaling quickly with system size. Fragment embedding is an attractive approach to this problem. By dividing a large complicated system into smaller manageable fragments "embedded" in an approximate description of the rest of the system, we can hope to ameliorate the steep cost of correlated calculations. While appealing, these methods often converge slowly with fragment size because of small errors at the boundary between fragment and bath. We describe a new electronic embedding method, dubbed "Bootstrap Embedding," a self-consistent wavefunction-in-wavefunction embedding theory that uses overlapping fragments to improve the description of fragment edges. We apply this method to the one dimensional Hubbard model and a translationally asymmetric variant, and find that it performs very well for energies and populations. We find Bootstrap Embedding converges rapidly with embedded fragment size, overcoming the surface-area-to-volume-ratio error typical of many embedding methods. We anticipate that this method may lead to a low-scaling, high accuracy treatment of electron correlation in large molecular systems.
This work presents the implementation of a fragment-based, quantum-centric supercomputing workflow for computing molecular electronic structure using quantum hardware. The workflow is applied to predict the relative energies of two conformers of the 300-atom Trp-cage miniprotein. The methodology employs wave function-based embedding (EWF) as the underlying fragmentation framework, in which all atoms in the system are explicitly included in the CI treatment. CI calculations for individual fragments are performed using either sample-based quantum diagonalization (SQD) for challenging fragments or full configuration interaction (FCI) for trivial fragments. To assess the accuracy of SQD for fragment CI calculations, EWF-(FCI,SQD) results are compared against EWF-MP2 and EWF-CCSD benchmarks. Overall, the results demonstrate that large-scale electronic configuration interaction (CI) simulations of protein systems containing hundreds or even thousands of atoms can be realized through the combined use of quantum and classical computing resources.