Rapid progress in noisy intermediate-scale quantum (NISQ) computing technology has led to the development of novel resource-efficient hybrid quantum-classical algorithms, such as the variational quantum eigensolver (VQE), that can address open challenges in quantum chemistry, physics, and material science. Proof-of-principle quantum chemistry simulations for small molecules have been demonstrated on NISQ devices. While several approaches have been theoretically proposed for correlated materials, NISQ simulations of interacting periodic models on current quantum devices have not yet been demonstrated. Here, we develop a hybrid quantum-classical simulation framework for correlated electron systems based on the Gutzwiller variational embedding approach. We implement this framework on Rigetti quantum processing units (QPUs) and apply it to the periodic Anderson model, which describes a correlated heavy electron band hybridizing with noninteracting conduction electrons. Our simulation results quantitatively reproduce the known ground state quantum phase diagram including metallic, Kondo and Mott insulating phases. This is the first fully self-consistent hybrid quantum-classical simulation of an infinite correlated lattice model executed on QPUs, demonstrating that the Gutzwiller hybrid quantum-classical embedding framework is a powerful approach to simulate correlated materials on NISQ hardware. This benchmark study also puts forth a concrete pathway towards practical quantum advantage on NISQ devices.
Rapid progress in noisy intermediate-scale quantum (NISQ) computing technology has led to the development of novel resource-efficient hybrid quantum-classical algorithms, such as the variational quantum eigensolver (VQE), that can address open challenges in quantum chemistry, physics, and material science. Proof-of-principle quantum chemistry simulations for small molecules have been demonstrated on NISQ devices. While several approaches have been theoretically proposed for correlated materials, NISQ simulations of interacting periodic models on current quantum devices have not yet been demonstrated. Here, we develop a hybrid quantum-classical simulation framework for correlated electron systems based on the Gutzwiller variational embedding approach. We implement this framework on Rigetti quantum processing units (QPUs) and apply it to the periodic Anderson model, which describes a correlated heavy electron band hybridizing with noninteracting conduction electrons. Our simulation results quantitatively reproduce the known ground state quantum phase diagram including metallic, Kondo and Mott insulating phases. This is the first fully self-consistent hybrid quantum-classical simulation of an infinite correlated lattice model executed on QPUs, demonstrating that the Gutzwiller hybrid quantum-classical embedding framework is a powerful approach to simulate correlated materials on NISQ hardware. This benchmark study also puts forth a concrete pathway towards practical quantum advantage on NISQ devices.
Rapid progress in noisy intermediate-scale quantum (NISQ) computing technology has led to the development of novel resource-efficient hybrid quantum-classical algorithms, such as the variational quantum eigensolver (VQE), that can address open challenges in quantum chemistry, physics, and material science.Proof-of-principle quantum chemistry simulations for small molecules have been demonstrated on NISQ devices.While several approaches have been theoretically proposed for correlated materials, NISQ simulations of interacting periodic models on current quantum devices have not yet been demonstrated.Here, we develop a hybrid quantumclassical simulation framework for correlated electron systems based on the Gutzwiller variational embedding approach.We implement this framework on Rigetti quantum processing units (QPUs) and apply it to the periodic Anderson model, which describes a correlated heavy electron band hybridizing with noninteracting conduction electrons.Our simulation results quantitatively reproduce the known ground state quantum phase diagram including metallic, Kondo and Mott insulating phases.This is the first fully self-consistent hybrid quantum-classical simulation of an infinite correlated lattice model executed on QPUs, demonstrating that the Gutzwiller hybrid quantum-classical embedding framework is a powerful approach to simulate correlated materials on NISQ hardware.This benchmark study also puts forth a concrete pathway towards practical quantum advantage on NISQ devices.
The accurate determination of the electronic structure of strongly correlated materials using first principle methods is of paramount importance in condensed matter physics, computational chemistry, and material science. However, due to the exponential scaling of computational resources, incorporating such materials into classical computation frameworks becomes prohibitively expensive. In 2016, Bauer et al. proposed a hybrid quantum-classical approach to correlated materials [B. Bauer et al., Hybrid quantum-classical approach to correlated materials, Phys. Rev. X 6, 031045 (2016).PRXHAE2160-330810.1103/PhysRevX.6.031045] that can efficiently tackle the electronic structure of complex correlated materials. Here, we experimentally demonstrate that approach to tackle the computational challenges associated with strongly correlated materials. By seamlessly integrating quantum computation into classical computers, we address the most computationally demanding aspect of the calculation, namely the computation of the Green's function, using a spin quantum processor. Furthermore, we realize a self-consistent determination of the single impurity Anderson model through a feedback loop between quantum and classical computations. A quantum phase transition in the Hubbard model from the metallic phase to the Mott insulator is observed as the strength of electron correlation increases. As the number of qubits with high control fidelity continues to grow, our experimental findings pave the way for solving even more complex models, such as strongly correlated crystalline materials or intricate molecules.
Scientific simulations demand methods combining scalability with predictive accuracy. Density Functional Theory (DFT) on High-Performance Computing (HPC) enables large-scale electronic-structure simulations but is limited by approximations affecting strongly correlated systems and band-gap predictions. Quantum computing offers a pathway to address this, though current Noisy Intermediate-Scale Quantum (NISQ) hardware remains constrained by qubit resources, noise, and execution cost.
This work presents a hybrid DFT-Quantum Embedding (QDFT) framework integrating classical HPC-based DFT with a quantum electronic-structure solver. Large systems are partitioned to isolate a chemically relevant active space, treated via the Variational Quantum Eigensolver (VQE), while the remaining degrees of freedom are described by DFT. The framework incorporates active-space selection, embedded Hamiltonian construction, symmetry preservation, operator mapping, self-consistent density updating, and modular classical-quantum coupling.
We focus on noiseless quantum simulation to systematically evaluate accuracy, convergence, active-space dependence, computational cost, and HPC scalability without hardware noise. Detailed profiling identifies computational bottlenecks and highlights limitations of CPU-based quantum simulation. A QPU runtime-estimation methodology is additionally developed to assess execution requirements on actual quantum hardware.
Results demonstrate quantum embedding's potential to improve selected electronic-structure properties while retaining classical HPC's scalability. Noisy quantum simulation and QPU execution remain key future directions, providing a pathway toward practical, scalable HPC-quantum hybrid simulations as hardware matures.
Understanding strongly correlated systems is essential for advancing quantum chemistry and materials science, yet conventional methods like Density Functional Theory (DFT) often fail to capture their complex electronic behavior. To address these limitations, we develop a hybrid quantum-classical framework that integrates Multiconfigurational Self-Consistent Field (MCSCF) with the Variational Quantum Eigensolver (VQE). Our initial benchmarks on water dissociation enabled the systematic optimization of key computational parameters, including ansatz selection, active space construction, and error mitigation. Building on this, we extend our approach to investigate the interactions between graphene analogues and water, demonstrating that our framework produces binding energies consistent with high-accuracy quantum methods. Furthermore, we apply this methodology to predict the binding energies of transition metals (Fe, Co, Ni) on both pristine and defective graphene analogues, revealing strong charge transfer effects and pronounced multireference character─phenomena often misrepresented by standard DFT. In contrast to many existing quantum algorithms that are constrained to small molecular systems, our framework achieves chemically accurate predictions for larger, strongly correlated systems such as metal-graphene complexes. This advancement highlights the capacity of hybrid quantum-classical approaches to address complex electronic interactions and demonstrates a practical route toward realizing quantum advantage for real-world materials applications in the Noisy Intermediate-Scale Quantum (NISQ) era.
Abstract Quantum computers promise to solve electronic structure, lattice dynamics, and strongly correlated problems beyond classical reach. However, near-term noisy intermediate scale quantum (NISQ) processors require systematic benchmarking to establish where variational algorithms with error mitigation can achieve chemical accuracy and cor rectly predict physical phenomena. Here we present 28 quantum material simula tions implemented with Qiskit, spanning molecular hydrogen (H 2 ) through the iron molybdenum cofactor (FeMoco), one-dimensional Hubbard chains, Kitaev honeycomb spin liquids, topological insulators (SSH model), BCS superconductors, and disordered Anderson insulators. Using the variational quantum eigensolver (VQE) with UCCSD ansatz, we achieve chemical accuracy (±1.6 mHa) for H 2 across the full dissociation curve and recover exact ground-state energies for H4 within 9 mHa after zero-noise extrapolation (ZNE) — reducing raw VQE errors from 152 mHa to 9 mHa. We accu rately identify the transverse-field Ising quantum critical point (h c = J) with 4 qubits, resolve Majorana Dirac dispersions in the Kitaev spin liquid, and reproduce the Mott transition in the Fermi-Hubbard model at U c = 6t. Resource estimation for FeMoco reveals O(N 2.5 ) scaling, projecting chemical accuracy with ∼200 logical qubits. Across all benchmarks, we quantify the accuracy-circuit depth tradeoff, demonstrating that depth-6 VQE ans¨atze achieve chemical accuracy for small molecules while depth-10 is required for correlated lattices. These results establish a comprehensive benchmark suite for evaluating quantum hardware, error mitigation protocols, and variational al gorithms for materials science on NISQ processors.
Abstract Simulating correlated materials on present-day quantum hardware remains challenging due to limited quantum resources. Quantum embedding methods offer a promising route by reducing computational complexity through the mapping of bulk systems onto effective impurity models, allowing more feasible simulations on pre- and early-fault-tolerant quantum devices. This work develops a quantum-classical embedding framework based on the ghost Gutzwiller approximation to enable quantum-enhanced simulations of ground-state properties and spectral functions of correlated electron systems. Circuit complexity is analyzed using an adaptive variational quantum algorithm on a statevector simulator, applied to the infinite-dimensional Hubbard model with increasing ghost mode numbers from 3 to 5, resulting in circuit depths growing from 16 to 104. Noise effects are examined using a realistic error model, revealing significant impact on the spectral weight of the Hubbard bands. To mitigate these effects, the Iceberg quantum error detection code is employed, achieving up to 40% error reduction in simulations. Finally, the accuracy of the density matrix estimation and the derived spectral function is benchmarked on IBM and Quantinuum quantum hardware, featuring distinct qubit-connectivity and employing multiple levels of error mitigation techniques.