Hybrid quantum-classical analog simulation of two-dimensional Fermi-Hubbard models with neutral atoms
Sergi Julià-Farré, Antoine Michel, Christophe Domain, Joseph Mikael, Lafoucriere, Jacques-Charles, Joseph Vovrosh, Ahmed Chahlaoui, Dorian Claveau, Guillaume Villaret, Julius de Hond, Loïc Henriet, Antoine Browaeys, Thomas Ayral, Alexandre Dauphin
We experimentally study the two-dimensional Fermi-Hubbard model using a Rydberg-based quantum processing unit in the analog mode. Our approach avoids encoding directly the original fermions into qubits and instead relies on reformulating the original model onto a system of fermions coupled to spins and then decoupling them in a self-consistent manner. We then introduce the auxiliary spin solver: this hybrid quantum-classical algorithm handles a free-fermion problem, which can be solved efficiently with a few classical resources, and an interacting spin problem, which can be naturally encoded in the analog quantum computer. This algorithm can be used to study both the equilibrium Mott transition as well as non-equilibrium properties of the original Fermi-Hubbard model, highlighting the potential of quantum-classical hybrid approaches to study strongly correlated matter.
We experimentally study the two-dimensional Fermi-Hubbard model using a Rydberg-based quantum processing unit in the analog mode. Our approach avoids encoding directly the original fermions into qubits and instead relies on reformulating the original model onto a system of fermions coupled to spins and then decoupling them in a self-consistent manner. We then introduce the auxiliary spin solver: this hybrid quantum-classical algorithm handles a free-fermion problem, which can be solved efficiently with a few classical resources, and an interacting spin problem, which can be naturally encoded in the analog quantum computer. This algorithm can be used to study both the equilibrium Mott transition as well as non-equilibrium properties of the original Fermi-Hubbard model, highlighting the potential of quantum-classical hybrid approaches to study strongly correlated matter.
We propose a hybrid quantum-classical method to investigate the equilibrium physics and the dynamics of strongly correlated fermionic models with spin-based quantum processors. Our proposal avoids the usual pitfalls of fermion-to-spin mappings thanks to a slave-spin method which allows to approximate the original Hamiltonian into a sum of self-correlated free fermions and spin Hamiltonians. Taking as an example a Rydberg-based analog quantum processor to solve the interacting spin model, we avoid the challenges of variational algorithms or Trotterization methods. We explore the robustness of the method to experimental imperfections by applying it to the half-filled, single-orbital Hubbard model on the square lattice in and out of equilibrium. We show, through realistic numerical simulations of current Rydberg processors, that the method yields quantitatively viable results even in the presence of imperfections: it allows to gain insights into equilibrium Mott physics as well as the dynamics under interaction quenches. This method thus paves the way to the investigation of physical regimes, whether out of equilibrium, doped, or multiorbital, that are difficult to explore with classical processors.
This work is a theoretical study of fermionic models. We focus on problems where highly controllable quantum simulators of these models have an important role, and we utilise both the analogue and the digital paradigm of quantum simulation. In the part on analogue quantum simulation, we focus on the proposed 'spin-asymmetric' Josephson effect where Cooper-paired spins display frequency synchronized Josephson oscillations with spin-dependent amplitudes. We consider different scenarios where the phenomenon could manifest in ultracold atomic Fermi gases. We study a Fermi gas Josephson junction in the recently realized Josephson plasma oscillation regime with an additional spin-dependent potential and show that the asymmetry in the resulting spin-dependent plasma oscillation amplitudes is on the order of a couple of per cent. We also demonstrate numerically that spin-asymmetric Josephson-like currents occur in a one-dimensional spin-dependent optical superlattice, with amplitude asymmetries up to 39%. Finally, we show that at zero temperature the tunable critical current in ferromagnetic Josephson junctions can be explained by the spin-asymmetric Josephson effect. In the part where digital quantum simulation is used, we propose a hybrid quantum-classical approach to studying strongly correlated fermion models. In this approach, a digital quantum simulator works in conjunction with a classical feedback loop to solve the infinite-dimensional Hubbard model directly in the thermodynamic limit. The scheme implements the well-established dynamical mean-field theory (DMFT) method, such that the digital quantum simulator solves the classically hard DMFT impurity problem and self-consistency is taken care of in a classical computer. We first present a few-qubit proof-of-principle setup for equilibrium systems that implements the simplified 'two-site' DMFT. This few-qubit setup is used for a qualitative description of the Mott transition in the half-filled infinite-dimensional Hubbard model. We then describe a scalable setup for simulating non-equilibrium many-body quantum dynamics by proposing the implementation of the non-equilibrium extension of DMFT with the hybrid device.
We consider a hybrid digital-analog quantum computing approach, which allows implementing any quantum algorithm without standard two-qubit gates. This approach is based on the always-on interaction between qubits, which can provide an alternative to such gates. We show how digital-analog approach can be applied to simulate the dynamics of fermionic systems, in particular, the Fermi-Hubbard model, using fermionic SWAP network and refocusing technique. We concentrate on the effects of connectivity topology, the spread of interaction constants as well as on errors of entangling operations. We find that an optimal connectivity topology of qubits for the digital-analog simulation of fermionic systems of arbitrary dimensionality is a chain for spinless fermions and a ladder for spin 1/2 particles. Such a simple connectivity topology makes digital-analog approach attractive for the simulation of quantum materials and molecules.
Quantum simulations of many-body systems are among the most promising applications of quantum computers1. In particular, models based on strongly correlated fermions are central to our understanding of quantum chemistry and materials problems2, and can lead to exotic, topological phases of matter3,4. However, owing to the non-local nature of fermions, such models are challenging to simulate with qubit devices5. Here we realize a digital quantum simulation architecture for two-dimensional fermionic systems based on reconfigurable atom arrays6. We utilize a fermion-to-qubit mapping based on Kitaev’s model on a honeycomb lattice3, in which fermionic statistics are encoded using long-range entangled states7. We prepare these states efficiently using measurement8 and feedforward9, realize subsequent fermionic evolution through Floquet engineering10,11 with tunable entangling gates12 interspersed with atom rearrangement, and improve results with built-in error detection. Leveraging this fermion description of the Kitaev spin model, we efficiently prepare topological states across its complex phase diagram13 and verify the non-Abelian spin-liquid phase3 by evaluating an odd Chern number14,15. We further explore this two-dimensional fermion system by realizing tunable dynamics and directly probing fermion exchange statistics. Finally, we simulate strong interactions and study the dynamics of the Fermi–Hubbard model on a square lattice. These results pave the way for digital quantum simulations of complex fermionic systems for materials science, chemistry16 and high-energy physics17. Digital quantum simulations of Kitaev’s honeycomb model are realized for two-dimensional fermionic systems using a reconfigurable atom-array processor and used to study the Fermi–Hubbard model on a square lattice.
The simulation of quantum many-body systems, relevant for quantum chemistry and condensed matter physics, is one of the most promising applications of near-term quantum computers before fault-tolerance. However, since the vast majority of quantum computing technologies are built around qubits and discrete gate-based operations, the translation of the physical problem into this framework is a crucial step. This translation will often be device specific, and a suboptimal implementation will be punished by the exponential compounding of errors on real devices. The importance of an efficient mapping is already revealed for models of spinful fermions in two or three dimensions, which naturally arise when the relevant physics relates to electrons. Using the most direct and well-known mapping, the Jordan-Wigner transformation, leads to a nonlocal representation of local degrees of freedom, and necessities efficient decompositions of nonlocal unitary gates into a sequence of hardware accessible local gates. In this paper, we provide a step-by-step recipe for simulating the paradigmatic two-dimensional Fermi-Hubbard model on a quantum computer using only local operations. To provide the ingredients for such a recipe, we briefly review the plethora of different approaches that have emerged recently but focus on the Derby-Klassen compact fermion mapping to make our discussion concrete. We provide a detailed recipe for an end-to-end simulation including embedding on a physical device, preparing initial states such as ground states, simulation of unitary time evolution, and measurement of observables and spectral functions. We explicitly compute the resource requirements for simulating a global quantum quench and conclude by discussing the challenges and future directions for simulating strongly correlated fermionic matter on quantum computers.
Quantum simulation of the XXZ model with a two-component Bose- or Fermi-Hubbard model based on a Mott insulator background has been widely used in the investigations of quantum magnetism with ultracold neutral atoms. In most cases, the diagonal spin-spin interaction is always accompanied by a large spin-exchange interaction which hinders the formation of long-range magnetic order at low temperature. Here, we show that the spin-exchange interaction can be strongly reduced in a Mott insulator of paired atoms, while the diagonal spin-spin interaction remains unaffected. Thus, the effective magnetic models become the exact quantum Ising models in a one-dimensional (1D) or two-dimensional (2D) lattice. Meanwhile, we analyzed an experimentally achievable three-component Fermi-Hubbard model of 6Li with two hyperfine levels of atoms paired in the lattice. We find the long-range antiferromagnetism of such a three-component Fermi-Hubbard model can be much stronger than that of a typical two-component Fermi-Hubbard model at low temperature. Our results may be useful for experimental investigations of the quantum phase transition and quantum criticality of the 1D and 2D quantum Ising models with ultracold neutral atoms.