Fermionic Hamiltonians for quantum simulations: a general reduction scheme

Many-body fermionic quantum calculations performed on analog quantum computers are restricted by the presence of k-local terms, which represent interactions among more than two qubits. These originate from the fermion-to-qubit mapping applied to the electronic Hamiltonians. Current solutions to this problem rely on perturbation theory in an enlarged Hilbert space. The main challenge associated with this technique is that it relies on coupling constants with very different magnitudes. This prevents its implementation in currently available architectures. In order to resolve this issue, we present an optimization scheme that unfolds the k-local terms into a linear combination of 2-local terms, while ensuring the conservation of all relevant physical properties of the original Hamiltonian, with several orders of magnitude smaller variation of the coupling constants.