Determining ground state energies of quantum systems by hybrid classical/quantum methods has emerged as a promising candidate application for near-term quantum computational resources. Short of large-scale fault-tolerant quantum computers, small-scale devices can be leveraged with current computational techniques to identify important subspaces of relatively large Hamiltonians. Inspired by the work that described the merging of dynamical mean-field theory (DMFT) with a small-scale quantum computational resource as an impurity solver [Bauer et al., arXiv:1510.03859v2], we describe an alternative embedding scheme, density matrix embedding theory (DMET), that naturally fits with the output from the variational quantum eigensolver and other hybrid approaches. This approach is validated using a quantum abstract machine simulator [Smith et al., arXiv:1608.03355] that reproduces the ground state energy of the Hubbard model converged to the infinite limit.
Determining ground state energies of quantum systems by hybrid classical/quantum methods has emerged as a promising candidate application for near-term quantum computational resources. Short of large-scale fault-tolerant quantum computers, small-scale devices can be leveraged with current computational techniques to identify important subspaces of relatively large Hamiltonians. Inspired by the work that described the merging of dynamical mean-field theory (DMFT) with a small-scale quantum computational resource as an impurity solver [Bauer et al., arXiv:1510.03859v2], we describe an alternative embedding scheme, density matrix embedding theory (DMET), that naturally fits with the output from the variational quantum eigensolver and other hybrid approaches. This approach is validated using a quantum abstract machine simulator [Smith et al., arXiv:1608.03355] that reproduces the ground state energy of the Hubbard model converged to the infinite limit.
Calculating the ground state properties of a Hamiltonian can be mapped to the problem of finding the ground state of a smaller Hamiltonian through the use of embedding methods. These embedding techniques have the ability to drastically reduce the problem size, and hence the number of qubits required when running on a quantum computer. However, the embedding process can produce a relatively complicated Hamiltonian, leading to a more complex quantum algorithm. In this paper we carry out a detailed study into how density matrix embedding theory (DMET) could be implemented on a quantum computer to solve the Hubbard model. We consider the variational quantum eigensolver (VQE) as the solver for the embedded Hamiltonian within the DMET algorithm. We derive the exact form of the embedded Hamiltonian and use it to construct efficient ansatz circuits and measurement schemes. We conduct detailed numerical simulations up to 16 qubits, the largest to date, for a range of Hubbard model parameters and find that the combination of DMET and VQE is effective for reproducing ground state properties of the model.
We introduce density matrix embedding theory (DMET), a quantum embedding theory for computing frequency-independent quantities, such as ground-state properties, of infinite systems. Like dynamical mean-field theory, DMET maps the bulk interacting system to a simpler impurity model and is exact in the noninteracting and atomic limits. Unlike dynamical mean-field theory, DMET is formulated in terms of the frequency-independent local density matrix, rather than the local Green's function. In addition, it features a finite, algebraically constructible bath of only one bath site per impurity site, with no bath discretization error. Frequency independence and the minimal bath make DMET a computationally simple and efficient method. We test the theory in the one-dimensional and two-dimensional Hubbard models at and away from half filling, and we find that compared to benchmark data, total energies, correlation functions, and metal-insulator transitions are well reproduced, at a tiny computational cost.
Accurately and efficiently predicting the equilibrium geometries of large molecules remains a central challenge in quantum computational chemistry, even with hybrid quantum–classical algorithms. Two major obstacles hinder progress: the large number of qubits required and the prohibitive cost of conventional nested optimization. In this work, we introduce a co-optimization framework that combines Density Matrix Embedding Theory (DMET) with Variational Quantum Eigensolver (VQE) to address these limitations. This approach substantially reduces the required quantum resources, enabling the treatment of molecular systems significantly larger than previously feasible. We first validate our framework on benchmark systems, such as H 4 and H 2 O 2, before demonstrating its efficacy in determining the equilibrium geometry of glycolic acid (C 2 H 4 O 3 )─a molecule of a size previously considered intractable for quantum geometry optimization. Our results show the method achieves high accuracy while drastically lowering computational cost. This work thus represents a significant step toward practical, scalable quantum simulations, moving beyond the small, proof-of-concept molecules that have historically dominated the field. More broadly, our framework establishes a tangible path toward leveraging quantum advantage for the in silico design of complex catalysts and pharmaceuticals.
Computing ground-state properties of molecules is a promising application for quantum computers operating in concert with classical high-performance computing resources. Quantum embedding methods are a family of algorithms particularly suited to these computational platforms: they combine high-level calculations on active regions of a molecule with low-level calculations on the surrounding environment, thereby avoiding expensive high-level full-molecule calculations and allowing the distribution of computational cost across multiple and heterogeneous computing units. Here, we present the first density matrix embedding theory (DMET) simulations performed in combination with a sample-based quantum diagonalization (SQD) method. We employ the DMET-SQD formalism to compute the ground-state energy of a ring of 18 hydrogen atoms and the relative energies of the chair, half-chair, twist-boat, and boat conformers of cyclohexane. The full-molecule 41- and 89-qubit simulations are decomposed into 27- and 32-qubit active-region simulations, which we carry out on the ibm_cleveland device, obtaining results in agreement with reference classical methods. Our DMET-SQD calculations mark tangible progress in the size of active regions that can be accurately tackled by near-term quantum computers and are an early demonstration of the potential for quantum-centric simulations to accurately treat the electronic structure of large molecules, with the ultimate goal of tackling systems such as peptides and proteins.
Computing ground-state properties of molecules is a promising application for quantum computers operating in concert with classical high-performance computing resources. Quantum embedding methods are a family of algorithms particularly suited to these computational platforms: they combine high-level calculations on active regions of a molecule with low-level calculations on the surrounding environment, thereby avoiding expensive high-level full-molecule calculations and allowing to distribute computational cost across multiple and heterogeneous computing units. Here, we present the first density matrix embedding theory (DMET) simulations performed in combination with the sample-based quantum diagonalization (SQD) method. We employ the DMET-SQD formalism to compute the ground-state energy of a ring of 18 hydrogen atoms, and the relative energies of the chair, half-chair, twist-boat, and boat conformers of cyclohexane. The full-molecule 41- and 89-qubit simulations are decomposed into 27- and 32-qubit active-region simulations, that we carry out on the ibm_cleveland device, obtaining results in agreement with reference classical methods. Our DMET-SQD calculations mark a tangible progress in the size of active regions that can be accurately tackled by near-term quantum computers, and are an early demonstration of the potential for quantum-centric simulations to accurately treat the electronic structure of large molecules, with the ultimate goal of tackling systems such as peptides and proteins.
Advances in quantum computing hold the potential to revolutionize the simulation of complex materials, including those critical for next-generation batteries. However, the complexity of typical battery materials is still prohibitive for exact simulations on current or near-term quantum hardware. One promising approach is the use of quantum embedding techniques, which enable the treatment of correlated electronic states in large systems by partitioning them into manageable subproblems. We explore the application of Density Matrix Embedding Theory (DMET) [1,2] as a quantum embedding method to bridge classical and quantum computational resources [3]. Within DMET, few correlated degrees of freedom can be solved on the quantum computer, while the remaining part is solved with classical methodology, self-consistently connected to the quantum part. This approach paves the way for real-world problems from battery simulation to profit from near-term quantum computers. Literature: [1] Knizia and Chan, Phys. Rev. Lett. 109, 186404 (2012) [2] Wouters et al., J. Chem. Theory Comput. 12, 2706 (2016) [3] Cao et al., npj Comput. Mater. 9, 78 (2023)