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:2506.05031 · EPJ Quantum Technology

Quantum simulation of the Hubbard model on a graphene hexagon: strengths of IQPE and noise constraints

M. Mirzakhani, Kyungsun Moon

Published Nov 18, 2025• 0 citations
arXivHTMLPDF

Abstract

Abstract Quantum computing offers transformative potential for simulating real-world materials, providing a powerful platform to investigate complex quantum systems across quantum chemistry and condensed matter physics. In this work, we leverage this capability to simulate the Hubbard model on a six-site graphene hexagon using Qiskit, employing the Iterative Quantum Phase Estimation (IQPE) and adiabatic evolution algorithms to determine its ground-state properties. Our results show that a single Slater determinant is sufficient to initialize IQPE and accurately recover ground-state energies (GSEs) in small-scale Hubbard systems. In noiseless simulations, IQPE converges within a few iterations to exact GSEs, while adiabatic simulations yield charge and spin densities and correlation functions in excellent agreement with exact diagonalization. However, deploying IQPE and adiabatic evolution on today's noisy quantum hardware remains highly challenging. To investigate these limitations in IQPE, we use the Qiskit Aer simulator with a custom noise model tailored to the characteristics of IBM's real hardware. This model includes realistic depolarizing gate errors, thermal relaxation, and readout noise, allowing us to explore how these factors degrade simulation accuracy. Further, we implement the IQPE algorithm on IBM's ibm_strasbourg and ibm_fez devices for a reduced three-site Hubbard model, enabling direct comparison between simulated and real hardware noise. While ibm_fez runs closely match exact results, discrepancies highlight the gap between modeled and physical noise. This study demonstrates both the IQPE's potential and current limitations for simulating strongly correlated systems under realistic conditions.

DOI
10.1140/epjqt/s40507-025-00439-x

Related papers

Quantum Simulation of the Fermi-Hubbard Model on a Graphene Hexagon

M. Mirzakhani, Kyungsun Moon

Quantum computing offers transformative potential for simulating real-world materials, providing a powerful platform to investigate complex quantum systems across quantum chemistry and condensed matter physics. In this work, we leverage this capability to simulate the Hubbard model on a six-site graphene hexagon using Qiskit, employing the Iterative Quantum Phase Estimation (IQPE) to determine its groundstate energies (GSEs). Noiseless simulations yield accurate GSEs in excellent agreement with exact diagonalization, underscoring the precision and reliability of quantum simulation for strongly correlated electron systems.

IEEE Conference on Pervasive and Intelligent Computing (PICom)•Oct 21, 2025•0 citations

Adiabatic quantum algorithm for artificial graphene

Axel Pérez-Obiol, Adrián Pérez-Salinas, Sergio Sanchez +2 more

We devise a quantum-circuit algorithm to solve the ground state and ground energy of artificial graphene. The algorithm implements a Trotterized adiabatic evolution from a purely tight-binding Hamiltonian to one including kinetic, spin-orbit, and Coulomb terms. The initial state is obtained efficiently using Gaussian-state preparation, while the readout of the ground energy is organized into seventeen sets of measurements, irrespective of the size of the problem. The total depth of the corresponding quantum circuit scales polynomially with the size of the system. A full simulation of the algorithm is performed and ground energies are obtained for lattices with up to four hexagons. Our results are benchmarked with exact diagonalization for systems with one and two hexagons. For larger systems we use the exact state vector and approximate matrix product state simulation techniques. The latter allows us to systematically trade off precision with memory and therefore to tackle larger systems. We analyze adiabatic and Trotterization errors, providing estimates for optimal periods and time discretizations given a finite accuracy. In the case of large systems we also study approximation errors.

arXiv:2204.03013•Physical Review A•Nov 8, 2022•8 citations

Ground-State Energy Estimation on Current Quantum Hardware through the Variational Quantum Eigensolver: A Practical Study

Nacer Eddine Belaloui, Abdellah Tounsi, Abdelmouheymen Rabah Khamadja +4 more

molecule, emphasizing practical implementation and performance on current quantum hardware. Our research presents a comparative study of HEA and UCCSD ansätze on noiseless and noisy simulations and implements VQE on recent IBM quantum computer noise models and a real quantum computer, IBM Fez, providing a fully functional code employing Qiskit 1.2. Our experiments confirm UCCSD's reliability in ideal conditions, while the HEA demonstrates greater robustness to hardware noise, achieving chemical accuracy on state-vector simulation (SVS). The results reveal that achieving ground-state energy within chemical accuracy is feasible without error mitigation during VQE convergence. We demonstrate that current quantum devices effectively optimize circuit parameters despite misestimating simulated energies. The SVS-evaluated energies provide a more accurate representation of the solution quality compared to QPU-estimated energy values, indicating that VQE converges to the correct ground state despite quantum noise. Our study also applies noise mitigation as a postprocessing technique, using zero-noise extrapolation (ZNE) on a real quantum computer. The detailed methodologies presented in this study, including Hamiltonian construction and Fermionic-to-qubit transformations, facilitate flexible adaptation of the VQE approach for various algorithm variants and across different levels of algorithmic implementation.

Journal of Chemical Theory and Computation•Jun 29, 2025•11 citations

Aplicación del algoritmo VQE al cálculo de la estructura de bandas del grafeno

Raúl Medrano-Millán

Graphene, a material of extraordinary interest in condensed matter physics due to its unique electronic properties, presents a fundamental challenge: the precise calculation of its band structure. This Master’s Thesis addresses this problem by proposing the use of one of the most widely used algorithms in quantum computing. To achieve this, a tight-binding model is used that describes the electron hopping between the atoms of the single-layer hexagonal lattice that defines graphene. This model is extended to include not only the interactions of each atom in the lattice with adjacent atoms (first neighbors), but also with the next closest atoms (second neighbors), allowing for the capture of more subtle physical effects. The resulting hamiltonian is then prepared for execution in the quantum simulators of IBM’s Qiskit platform. The results under ideal conditions demonstrate that the algorithm reproduces analytical solutions with high fidelity, capturing essential features such as Dirac points, where valence and conduction bands intersect, causing electrons to behave as massless particles, and electron-hole symmetry breaking, which reflects the difference in the energetic behavior of electrons and their absences (holes). However, when introducing a noise model that emulates current quantum processors (inherently noisy, intermediate-scale devices, the socalled NISQ era), a systematic deviation is observed that underscores the limiting impact of the hardware. The study concludes that VQE is a conceptually viable and powerful tool for materials simulation, but its practical application to achieve quantum advantage will critically depend on the development of advanced error mitigation techniques.

•Jul 23, 2025•0 citations

Classical Benchmarks for Variational Quantum Eigensolver Simulations of the Hubbard Model

Antonios M. Alvertis, Abid Khan, Thomas Iadecola +2 more

Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.

arXiv:2408.00836•Quantum•May 20, 2025•4 citations

Simulating strongly interacting Hubbard chains with the variational Hamiltonian ansatz on a quantum computer

Baptiste Anselme Martin, Pascal Simon, Marko J. Rančić

Hybrid quantum-classical algorithms have been proposed to circumvent noise limitations in quantum computers. Such algorithms delegate only a calculation of the expectation value to the quantum computer. Among them, the variational quantum eigensolver has been implemented to study molecules and condensed matter systems on small size quantum computers. Condensed matter systems described by the Hubbard model exhibit a rich phase diagram alongside exotic states of matter. In this paper we try to answer the question: How much of the underlying physics of a 1D Hubbard chain is described by a problem-inspired variational Hamiltonian ansatz in a broad range of parameter values? We start by probing how much the solution increases fidelity with increasing ansatz complexity. Our findings suggest that even low fidelity solutions capture energy and number of doubly occupied sites well, while spin-spin correlations are not well captured even when the solution is of high fidelity. Our powerful simulation platform allows us to incorporate a realistic noise model and shows a successful implementation of noise-mitigation strategies---postselection and the Richardson extrapolation. Finally, we compare our results with an experimental realization of the algorithm on IBM Quantum's ibmq_quito device.

arXiv:2111.11996•Physical Review Research•Jun 6, 2022•20 citations
Abstract Quantum computing offers transformative potential for simulating real-world materials, providing a powerful platform to investigate complex quantum systems across quantum chemistry and condensed matter physics. In this work, we leverage this capability to simulate the Hubbard model on a six-site graphene hexagon using Qiskit, employing the Iterative Quantum Phase Estimation (IQPE) and adiabatic evolution algorithms to determine its ground-state properties. Our results show that a single Slater determinant is sufficient to initialize IQPE and accurately recover ground-state energies (GSEs) in small-scale Hubbard systems. In noiseless simulations, IQPE converges within a few iterations to exact GSEs, while adiabatic simulations yield charge and spin densities and correlation functions in excellent agreement with exact diagonalization. However, deploying IQPE and adiabatic evolution on today's noisy quantum hardware remains highly challenging. To investigate these limitations in IQPE, we use the Qiskit Aer simulator with a custom noise model tailored to the characteristics of IBM's real hardware. This model includes realistic depolarizing gate errors, thermal relaxation, and readout noise, allowing us to explore how these factors degrade simulation accuracy. Further, we implement the IQPE algorithm on IBM's ibm_strasbourg and ibm_fez devices for a reduced three-site Hubbard model, enabling direct comparison between simulated and real hardware noise. While ibm_fez runs closely match exact results, discrepancies highlight the gap between modeled and physical noise. This study demonstrates both the IQPE's potential and current limitations for simulating strongly correlated systems under realistic conditions.
Quantum computing offers transformative potential for simulating real-world materials, providing a powerful platform to investigate complex quantum systems across quantum chemistry and condensed matter physics. In this work, we leverage this capability to simulate the Hubbard model on a six-site graphene hexagon using Qiskit, employing the Iterative Quantum Phase Estimation (IQPE) to determine its groundstate energies (GSEs). Noiseless simulations yield accurate GSEs in excellent agreement with exact diagonalization, underscoring the precision and reliability of quantum simulation for strongly correlated electron systems.
We devise a quantum-circuit algorithm to solve the ground state and ground energy of artificial graphene. The algorithm implements a Trotterized adiabatic evolution from a purely tight-binding Hamiltonian to one including kinetic, spin-orbit, and Coulomb terms. The initial state is obtained efficiently using Gaussian-state preparation, while the readout of the ground energy is organized into seventeen sets of measurements, irrespective of the size of the problem. The total depth of the corresponding quantum circuit scales polynomially with the size of the system. A full simulation of the algorithm is performed and ground energies are obtained for lattices with up to four hexagons. Our results are benchmarked with exact diagonalization for systems with one and two hexagons. For larger systems we use the exact state vector and approximate matrix product state simulation techniques. The latter allows us to systematically trade off precision with memory and therefore to tackle larger systems. We analyze adiabatic and Trotterization errors, providing estimates for optimal periods and time discretizations given a finite accuracy. In the case of large systems we also study approximation errors.
molecule, emphasizing practical implementation and performance on current quantum hardware. Our research presents a comparative study of HEA and UCCSD ansätze on noiseless and noisy simulations and implements VQE on recent IBM quantum computer noise models and a real quantum computer, IBM Fez, providing a fully functional code employing Qiskit 1.2. Our experiments confirm UCCSD's reliability in ideal conditions, while the HEA demonstrates greater robustness to hardware noise, achieving chemical accuracy on state-vector simulation (SVS). The results reveal that achieving ground-state energy within chemical accuracy is feasible without error mitigation during VQE convergence. We demonstrate that current quantum devices effectively optimize circuit parameters despite misestimating simulated energies. The SVS-evaluated energies provide a more accurate representation of the solution quality compared to QPU-estimated energy values, indicating that VQE converges to the correct ground state despite quantum noise. Our study also applies noise mitigation as a postprocessing technique, using zero-noise extrapolation (ZNE) on a real quantum computer. The detailed methodologies presented in this study, including Hamiltonian construction and Fermionic-to-qubit transformations, facilitate flexible adaptation of the VQE approach for various algorithm variants and across different levels of algorithmic implementation.
Graphene, a material of extraordinary interest in condensed matter physics due to its unique electronic properties, presents a fundamental challenge: the precise calculation of its band structure. This Master’s Thesis addresses this problem by proposing the use of one of the most widely used algorithms in quantum computing. To achieve this, a tight-binding model is used that describes the electron hopping between the atoms of the single-layer hexagonal lattice that defines graphene. This model is extended to include not only the interactions of each atom in the lattice with adjacent atoms (first neighbors), but also with the next closest atoms (second neighbors), allowing for the capture of more subtle physical effects. The resulting hamiltonian is then prepared for execution in the quantum simulators of IBM’s Qiskit platform. The results under ideal conditions demonstrate that the algorithm reproduces analytical solutions with high fidelity, capturing essential features such as Dirac points, where valence and conduction bands intersect, causing electrons to behave as massless particles, and electron-hole symmetry breaking, which reflects the difference in the energetic behavior of electrons and their absences (holes). However, when introducing a noise model that emulates current quantum processors (inherently noisy, intermediate-scale devices, the socalled NISQ era), a systematic deviation is observed that underscores the limiting impact of the hardware. The study concludes that VQE is a conceptually viable and powerful tool for materials simulation, but its practical application to achieve quantum advantage will critically depend on the development of advanced error mitigation techniques.
Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.
Hybrid quantum-classical algorithms have been proposed to circumvent noise limitations in quantum computers. Such algorithms delegate only a calculation of the expectation value to the quantum computer. Among them, the variational quantum eigensolver has been implemented to study molecules and condensed matter systems on small size quantum computers. Condensed matter systems described by the Hubbard model exhibit a rich phase diagram alongside exotic states of matter. In this paper we try to answer the question: How much of the underlying physics of a 1D Hubbard chain is described by a problem-inspired variational Hamiltonian ansatz in a broad range of parameter values? We start by probing how much the solution increases fidelity with increasing ansatz complexity. Our findings suggest that even low fidelity solutions capture energy and number of doubly occupied sites well, while spin-spin correlations are not well captured even when the solution is of high fidelity. Our powerful simulation platform allows us to incorporate a realistic noise model and shows a successful implementation of noise-mitigation strategies---postselection and the Richardson extrapolation. Finally, we compare our results with an experimental realization of the algorithm on IBM Quantum's ibmq_quito device.