Efficient Simulation of Sparse, Non-Local Fermion Models
Reinis Irmejs, J. Ignacio Cirac
TL;DR
The paper tackles the challenge of simulating sparse, non-local fermionic models on all-to-all connected qubit hardware by introducing an auxiliary-fermion encoding that eliminates Jordan–Wigner strings. Stabilizers built from auxiliary Majorana modes render interaction terms into local Pauli operators with constant weight, while initializing the auxiliary sector in a mutual +1 eigenstate preserves the original physics during dynamics. A graph-coloring-based grouping of stabilizers minimizes ancilla requirements to $\nu=\lceil \chi/2 \rceil$ per site, and an ordered-state preparation enables efficient initialization in $\mathcal{O}(\log\nu)$ depth. The results show that long-time dynamics can be implemented with asymptotically optimal circuit depth, matching the performance of fermionic hardware up to constants and offering a practical route for simulating sparse fermionic systems on qubit-based quantum computers. The approach reduces the previous $O(\log N)$ multiplicative overhead to an additive overhead, with total ancilla scaling as $O(dN)$, thereby closing a key gap in the quantum simulation of sparse fermionic models.
Abstract
Efficient simulation of interacting fermionic systems is a key application of near-term quantum computers, but is hindered by the overhead required to encode fermionic operators on qubit hardware. Here, we consider models with $N$ fermionic modes in which each participates in at most a constant number $d$ of interactions and study the circuit depth required to implement the Trotterized time evolution on qubit hardware with all-to-all connectivity. We introduce an encoding that augments each physical fermionic mode with a small number of auxiliary fermions, enabling the removal of Jordan--Wigner strings. Although the preparation of the auxiliary fermion state incurs an initial overhead, this state remains invariant under time evolution. As a result, long-time evolution can be implemented with asymptotically optimal circuit depth, reducing a previously multiplicative $O(\log N)$ overhead to an additive overhead. Our results thus establish that the simulation of sparse fermionic models on qubit hardware matches the performance achievable on ideal fermionic hardware up to constant factors and $O(dN)$ ancillary qubits.
