Classical simulation of circuits with realistic Gottesman-Kitaev-Preskill states

Abstract: Classically simulating circuits with bosonic codes is a challenging task due to the prohibitive cost of simulating quantum systems with many, possibly infinite, energy levels. We propose an algorithm to simulate circuits with encoded Gottesman-Kitaev-Preskill states, specifically for odd-dimensional encoded qudits. Our approach is tailored to be especially effective in the most challenging but practically relevant regime, where the codeword states exhibit high (but finite) squeezing. Our algorithm leverages the Zak-Gross Wigner function introduced by J. Davis et al. [arXiv:2407.18394], which represents infinitely squeezed encoded stabilizer states positively. The runtime of the algorithm scales with the amount of negativity of this Wigner function, enabling fast simulation of certain large-scale circuits with a high degree of squeezing.