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.
