Initial states and apodisation for quantum field simulations in phase-space

TL;DR

This paper reviews phase-space methods for simulating bosonic quantum fields, focusing on initial states, boundary treatments, and apodisation techniques to improve accuracy in quantum dynamical calculations.

Contribution

It introduces complex apodisation methods with phase-shifts to enhance boundary accuracy and discusses efficient sampling of initial quantum states in phase-space simulations.

Findings

Complex power law apodisation reduces boundary errors.

Efficient sampling of multi-mode Gaussian and number states.

Method for analyzing number conservation in apodised simulations.

Abstract

Bosonic quantum fields can be simulated with `quantum software' in phase-space. The positive-P, Wigner and Q-function phase-space methods are reviewed. Initial quantum states and boundaries for infinite domains are considered in detail. The quantum initial conditions treated include both multi-mode Gaussian states and number states, which are sampled using phase-shifted weighting techniques. This is more efficient than direct probabilistic sampling. Algorithms for treating periodic boundaries within infinite domains are developed via apodisation, or absorbers near the boundary. Similar truncation or anti-aliasing issues arise in the frequency domain with Fourier transform methods. Complex apodisation, in which there are phase-shifts as well as absorption, improves accuracy by reducing unwanted errors, and this is illustrated numerically. The complex power law apodisation techniques described here can be applied to numerical calculations or even directly in optics. Phase-space simulations, used for both quantum optical and Bose-Einstein condensate quantum dynamical calculations, require quantum noise terms to be included. A method is described to analyze number conservation in apodised calculations. Some techniques developed here are potentially applicable to all digital simulations, whether using classical or quantum logic.