First-Quantized Quantum Simulation of Non-Relativistic QED with Emergent Topologically Protected Coulomb Interactions

TL;DR

This paper introduces a quantum simulation algorithm for non-relativistic QED where Coulomb interactions emerge from Gauss' law constraints, featuring topological protection that reduces errors and potentially offers computational advantages for large systems.

Contribution

The authors develop a first-quantized quantum simulation method for non-relativistic QED with emergent Coulomb interactions and topological error protection, differing from previous explicit Coulomb models.

Findings

Coulomb interactions emerge from Gauss' law constraints.

Topological protection prevents certain electric field errors.

The algorithm scales favorably compared to explicit Coulomb simulations under specific conditions.

Abstract

We provide a simulation algorithm that properly addresses light matter interaction between non-relativistic first-quantized charged particles and quantum electromagnetic fields. Unlike previous work, our Hamiltonian does not include an explicit Coulomb interaction between particles. Rather, the Coulomb interaction emerges from the imposition of Gauss' law as a constraint upon the system in an appropriate non-relativistic limit. Furthermore, a form of topological protection emerges in our formalism, analogous to that of the Toric code Hamiltonian. This mechanism prevents simulation-induced electric field errors that can be contracted to a point from causing any deviations from Coulomb's law in the non-relativistic limit and any error that forms a non-contractable loop is energetically dissallowed in the limit of large volume. We find that, under appropriate continuity assumptions, the number of non-Clifford gates required by our algorithm scales in the thermodynamic limit as $\widetilde{O}(N^{2/3}\eta^{4/3} t \log^5(1/\epsilon))$ for $\eta$ particles, $N$ spatial grid points, simulation time $t$ and error tolerance $\epsilon$. In comparison, the more specific problem of simulating the Coulomb interaction scales as $\widetilde{O}(N^{1/3} \eta^{8/3} t \log^2(1/\epsilon))$. This suggests that if $N \in \tilde{o}(\eta^4)$ that our non-relativistic electrodynamic simulation method could provide a computational advantage for electronic structure problems in the thermodynamic limit under appropriate continuity assumptions as it obviates the need to compute the $O(\eta^2)$ pairwise interactions in the Coulomb Hamiltonian.