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:2405.09225 · npj Quantum Materials

Exploring ground states of Fermi-Hubbard model on honeycomb lattices with counterdiabaticity

Jialiang Tang, Ruoqian Xu, Yongcheng Ding, Xusheng Xu, Yue Ban, Man‐Hong Yung, Axel Pérez-Obiol, Gloria Platero, Xi Chen

Published Nov 7, 2024• 5 citations
arXivHTMLPDF

Abstract

Exploring the ground state properties of many-body quantum systems conventionally involves adiabatic processes, alongside exact diagonalization, in the context of quantum annealing or adiabatic quantum computation. Shortcuts to adiabaticity by counter-diabatic driving serve to accelerate these processes by suppressing energy excitations. Motivated by this, we develop variational quantum algorithms incorporating the auxiliary counter-diabatic interactions, comparing them with digitized adiabatic algorithms. These algorithms are then implemented on gate-based quantum circuits to explore the ground states of the Fermi-Hubbard model on honeycomb lattices, utilizing systems with up to 26 qubits. The comparison reveals that the counter-diabatic inspired ansatz is superior to traditional Hamiltonian variational ansatz. Furthermore, the number and duration of Trotter steps are analyzed to understand and mitigate errors. Given the model’s relevance to materials in condensed matter, our study paves the way for using variational quantum algorithms with counterdiabaticity to explore quantum materials in the noisy intermediate-scale quantum era.

DOI
10.1038/s41535-024-00697-5

Related papers

Programming Optical-Lattice Fermi-Hubbard Quantum Simulators

Cristian Tabares, Christian Kokail, P. Zoller +2 more

Fermionic atoms in optical lattices provide a native implementation of Fermi-Hubbard (FH) models that can be used as analog quantum simulators of many-body fermionic systems. Recent experimental advances include the time-dependent local control of chemical potentials and tunnelings, and thus enable us to operate this platform digitally as a programmable quantum simulator. Here, we explore these opportunities and develop ground-state preparation algorithms for different fermionic models, based on the ability to implement both single-particle and many-body high-fidelity fermionic gates, as provided by the native FH Hamiltonian. In particular, we first design variational precompiled quantum circuits to prepare the ground state of the natively implemented FH model, with significant speedups relative to competing adiabatic protocols. Besides, the versatility of this variational approach enables us to target extended FH models, i.e., including terms that are not natively realized on the platform. As an illustration, we include next-nearest-neighbor tunnelings at finite dopings, relevant in the context of d -wave superconductivity. Furthermore, we discuss how to approximate the imaginary-time evolution using variational fermionic circuits, both as an alternative state-preparation strategy and as a subroutine for the quantum Lanczos algorithm to further improve the energy estimation. We benchmark our protocols for ladder geometries, though they can be readily applied to two-dimensional experimental setups to address regimes beyond the capabilities of current classical methods. These results pave the way for more efficient and comprehensive explorations of relevant many-body phases with existing programmable fermionic quantum simulators.

arXiv:2502.05067•PRX Quantum•Jul 21, 2025•4 citations

Shortcuts to Adiabaticity in the Infinite-Range Ising Model by Mean-Field Counter-Diabatic Driving

Takuya Hatomura

The strategy of shortcuts to adiabaticity enables us to realize adiabatic dynamics in finite time. In the counter-diabatic driving approach, an auxiliary Hamiltonian which is called the counter-diabatic Hamiltonian is appended to an original Hamiltonian to cancel out diabatic transitions. The counter-diabatic Hamiltonian is constructed by using the eigenstates of the original Hamiltonian. Therefore, it is in general difficult to construct the counter-diabatic Hamiltonian for quantum many-body systems. Even if the counter-diabatic Hamiltonian for quantum many-body systems is obtained, it is generally non-local and even diverges at critical points. We construct an approximated counter-diabatic Hamiltonian for the infinite-range Ising model by making use of the mean-field approximation. An advantage of this method is that the mean-field counter-diabatic Hamiltonian is constructed by only local operators. We numerically demonstrate the effectiveness of this method through quantum annealing processes going the vicinity of the critical point. It is also confirmed that the mean-field counter-diabatic Hamiltonian is still well-defined in the limit to the critical point. The present method can take higher order contributions into account and is consistent with the variational approach for local counter-diabatic driving.

arXiv:1705.03168•Journal of the Physical Society of Japan•Aug 23, 2017•24 citations

Adaptive variational preparation of the Fermi-Hubbard eigenstates

Gaurav Gyawali, Michael J. Lawler

Approximating the ground states of strongly interacting electron systems in quantum chemistry and condensed matter physics is expected to be one of the earliest applications of quantum computers. In this paper, we prepare highly accurate ground states of the Fermi-Hubbard model for small grids of up to six sites (12 qubits) by using an interpretable, adaptive variational quantum eigensolver (VQE) called ADAPT-VQE [H. R. Grimsley et al., Nat. Commun. 10, 3007 (2019)]. In contrast with nonadaptive VQE, this algorithm builds a system-specific ansatz by adding an optimal gate built from one-body or two-body fermionic operators at each step. We show this adaptive method outperforms the nonadaptive counterpart in terms of fewer variational parameters, short gate depth, and scaling with the system size. The fidelity and energy of the prepared state appear to improve asymptotically with ansatz depth. We also demonstrate the application of adaptive variational methods by preparing excited states and Green's functions using a proposed ADAPT-SSVQE algorithm. Lower depth, asymptotic convergence, noise tolerance of a variational approach [A. Peruzzo et al., Nat. Commun. 5, 4213 (2014); J. R. McLean et al., New J. Phys. 18, 023023 (2016); K. Sharma et al., New J. Phys. 22, 043006 (2020)], and a highly controllable, system-specific ansatz make the adaptive variational methods particularly well suited for noisy intermediate-scale quantum devices.

arXiv:2109.12126•Physical Review A•Jan 6, 2022•12 citations

Finding the ground state of the Hubbard model by variational methods on a quantum computer with gate errors

Jan-Michael Reiner, Frank K. Wilhelm, Gerd Schön +1 more

Abstract A key goal of digital quantum computing is the simulation of fermionic systems such as molecules or the Hubbard model. Unfortunately, for present and near-future quantum computers the use of quantum error correction schemes is still out of reach. Hence, the finite error rate limits the use of quantum computers to algorithms with a low number of gates. The variational Hamiltonian ansatz (VHA) has been shown to produce the ground state in good approximation in a manageable number of steps. Here we study explicitly the effect of gate errors on its performance. The VHA is inspired by the adiabatic quantum evolution under the influence of a time-dependent Hamiltonian, where the—ideally short—fixed Trotter time steps are replaced by variational parameters. The method profits substantially from quantum variational error suppression, e.g. unitary quasi-static errors are mitigated within the algorithm. We test the performance of the VHA when applied to the Hubbard model in the presence of unitary control errors on quantum computers with realistic gate fidelities.

arXiv:1811.04476•Quantum Science and Technology•May 1, 2019•49 citations

Optimal quantum annealing: A variational shortcut-to-adiabaticity approach

Gianluca Passarelli, Rosario Fazio, Procolo Lucignano

Suppressing unwanted transitions out of the instantaneous ground state is a major challenge in unitary adiabatic quantum computation. A recent approach consists in building counterdiabatic potentials approximated using variational strategies. In this contribution, we extend this variational approach to Lindbladian dynamics, having as a goal the suppression of diabatic transitions between pairs of Jordan blocks in quantum annealing. We show that, surprisingly, unitary counterdiabatic Ans\"atze are successful for dissipative dynamics as well, allowing for easier experimental implementations compared to Lindbladian Ans\"atze involving dissipation. Our approach not only guarantees improvements of open-system adiabaticity but also enhances the success probability of quantum annealing.

arXiv:2109.13043•Physical Review A•Feb 23, 2022•29 citations

Variational counterdiabatic driving of the Hubbard model for ground-state preparation

Qing Xie, Kazuhiro Seki, Seiji Yunoki

Counterdiabatic (CD) protocols enable fast driving of quantum states by invoking an auxiliary adiabatic gauge potential (AGP) that suppresses transitions to excited states throughout the driving process. Usually, the full spectrum of the original unassisted Hamiltonian is a prerequisite for constructing the exact AGP, which implies that CD protocols are extremely difficult for many-body systems. Here, we apply a variational CD protocol recently proposed by Claeys et al. [Phys. Rev. Lett. 123, 090602 (2019)] to a two-component fermionic Hubbard model in one spatial dimension. This protocol engages an approximated AGP expressed as a series of nested commutators. We show that the optimal variational parameters in the approximated AGP satisfy a set of linear equations whose coefficients are given by the squared Frobenius norms of these commutators. We devise an exact algorithm that escapes the formidable iterative matrix-vector multiplications and evaluates the nested commutators and the CD Hamiltonian in analytic representations. We then examine the CD driving of the one-dimensional Hubbard model up to $L=14$ sites with driving order $l\ensuremath{\leqslant}3$. Our results demonstrate the usefulness of the variational CD protocol to the Hubbard model and permit a possible route towards fast ground-state preparation for many-body systems.

arXiv:2206.07597•Physical review. B./Physical review. B•Oct 28, 2022•22 citations
Exploring the ground state properties of many-body quantum systems conventionally involves adiabatic processes, alongside exact diagonalization, in the context of quantum annealing or adiabatic quantum computation. Shortcuts to adiabaticity by counter-diabatic driving serve to accelerate these processes by suppressing energy excitations. Motivated by this, we develop variational quantum algorithms incorporating the auxiliary counter-diabatic interactions, comparing them with digitized adiabatic algorithms. These algorithms are then implemented on gate-based quantum circuits to explore the ground states of the Fermi-Hubbard model on honeycomb lattices, utilizing systems with up to 26 qubits. The comparison reveals that the counter-diabatic inspired ansatz is superior to traditional Hamiltonian variational ansatz. Furthermore, the number and duration of Trotter steps are analyzed to understand and mitigate errors. Given the model’s relevance to materials in condensed matter, our study paves the way for using variational quantum algorithms with counterdiabaticity to explore quantum materials in the noisy intermediate-scale quantum era.
Fermionic atoms in optical lattices provide a native implementation of Fermi-Hubbard (FH) models that can be used as analog quantum simulators of many-body fermionic systems. Recent experimental advances include the time-dependent local control of chemical potentials and tunnelings, and thus enable us to operate this platform digitally as a programmable quantum simulator. Here, we explore these opportunities and develop ground-state preparation algorithms for different fermionic models, based on the ability to implement both single-particle and many-body high-fidelity fermionic gates, as provided by the native FH Hamiltonian. In particular, we first design variational precompiled quantum circuits to prepare the ground state of the natively implemented FH model, with significant speedups relative to competing adiabatic protocols. Besides, the versatility of this variational approach enables us to target extended FH models, i.e., including terms that are not natively realized on the platform. As an illustration, we include next-nearest-neighbor tunnelings at finite dopings, relevant in the context of d -wave superconductivity. Furthermore, we discuss how to approximate the imaginary-time evolution using variational fermionic circuits, both as an alternative state-preparation strategy and as a subroutine for the quantum Lanczos algorithm to further improve the energy estimation. We benchmark our protocols for ladder geometries, though they can be readily applied to two-dimensional experimental setups to address regimes beyond the capabilities of current classical methods. These results pave the way for more efficient and comprehensive explorations of relevant many-body phases with existing programmable fermionic quantum simulators.
The strategy of shortcuts to adiabaticity enables us to realize adiabatic dynamics in finite time. In the counter-diabatic driving approach, an auxiliary Hamiltonian which is called the counter-diabatic Hamiltonian is appended to an original Hamiltonian to cancel out diabatic transitions. The counter-diabatic Hamiltonian is constructed by using the eigenstates of the original Hamiltonian. Therefore, it is in general difficult to construct the counter-diabatic Hamiltonian for quantum many-body systems. Even if the counter-diabatic Hamiltonian for quantum many-body systems is obtained, it is generally non-local and even diverges at critical points. We construct an approximated counter-diabatic Hamiltonian for the infinite-range Ising model by making use of the mean-field approximation. An advantage of this method is that the mean-field counter-diabatic Hamiltonian is constructed by only local operators. We numerically demonstrate the effectiveness of this method through quantum annealing processes going the vicinity of the critical point. It is also confirmed that the mean-field counter-diabatic Hamiltonian is still well-defined in the limit to the critical point. The present method can take higher order contributions into account and is consistent with the variational approach for local counter-diabatic driving.
Approximating the ground states of strongly interacting electron systems in quantum chemistry and condensed matter physics is expected to be one of the earliest applications of quantum computers. In this paper, we prepare highly accurate ground states of the Fermi-Hubbard model for small grids of up to six sites (12 qubits) by using an interpretable, adaptive variational quantum eigensolver (VQE) called ADAPT-VQE [H. R. Grimsley et al., Nat. Commun. 10, 3007 (2019)]. In contrast with nonadaptive VQE, this algorithm builds a system-specific ansatz by adding an optimal gate built from one-body or two-body fermionic operators at each step. We show this adaptive method outperforms the nonadaptive counterpart in terms of fewer variational parameters, short gate depth, and scaling with the system size. The fidelity and energy of the prepared state appear to improve asymptotically with ansatz depth. We also demonstrate the application of adaptive variational methods by preparing excited states and Green's functions using a proposed ADAPT-SSVQE algorithm. Lower depth, asymptotic convergence, noise tolerance of a variational approach [A. Peruzzo et al., Nat. Commun. 5, 4213 (2014); J. R. McLean et al., New J. Phys. 18, 023023 (2016); K. Sharma et al., New J. Phys. 22, 043006 (2020)], and a highly controllable, system-specific ansatz make the adaptive variational methods particularly well suited for noisy intermediate-scale quantum devices.
Abstract A key goal of digital quantum computing is the simulation of fermionic systems such as molecules or the Hubbard model. Unfortunately, for present and near-future quantum computers the use of quantum error correction schemes is still out of reach. Hence, the finite error rate limits the use of quantum computers to algorithms with a low number of gates. The variational Hamiltonian ansatz (VHA) has been shown to produce the ground state in good approximation in a manageable number of steps. Here we study explicitly the effect of gate errors on its performance. The VHA is inspired by the adiabatic quantum evolution under the influence of a time-dependent Hamiltonian, where the—ideally short—fixed Trotter time steps are replaced by variational parameters. The method profits substantially from quantum variational error suppression, e.g. unitary quasi-static errors are mitigated within the algorithm. We test the performance of the VHA when applied to the Hubbard model in the presence of unitary control errors on quantum computers with realistic gate fidelities.
Suppressing unwanted transitions out of the instantaneous ground state is a major challenge in unitary adiabatic quantum computation. A recent approach consists in building counterdiabatic potentials approximated using variational strategies. In this contribution, we extend this variational approach to Lindbladian dynamics, having as a goal the suppression of diabatic transitions between pairs of Jordan blocks in quantum annealing. We show that, surprisingly, unitary counterdiabatic Ans\"atze are successful for dissipative dynamics as well, allowing for easier experimental implementations compared to Lindbladian Ans\"atze involving dissipation. Our approach not only guarantees improvements of open-system adiabaticity but also enhances the success probability of quantum annealing.
Counterdiabatic (CD) protocols enable fast driving of quantum states by invoking an auxiliary adiabatic gauge potential (AGP) that suppresses transitions to excited states throughout the driving process. Usually, the full spectrum of the original unassisted Hamiltonian is a prerequisite for constructing the exact AGP, which implies that CD protocols are extremely difficult for many-body systems. Here, we apply a variational CD protocol recently proposed by Claeys et al. [Phys. Rev. Lett. 123, 090602 (2019)] to a two-component fermionic Hubbard model in one spatial dimension. This protocol engages an approximated AGP expressed as a series of nested commutators. We show that the optimal variational parameters in the approximated AGP satisfy a set of linear equations whose coefficients are given by the squared Frobenius norms of these commutators. We devise an exact algorithm that escapes the formidable iterative matrix-vector multiplications and evaluates the nested commutators and the CD Hamiltonian in analytic representations. We then examine the CD driving of the one-dimensional Hubbard model up to L=14L=14L=14 sites with driving order $l\ensuremath{\leqslant}3$. Our results demonstrate the usefulness of the variational CD protocol to the Hubbard model and permit a possible route towards fast ground-state preparation for many-body systems.