Optimization of the positive-P representation for the anharmonic oscillator

We demonstrate how the freedom in the definition of the positive-P representation may be used to reduce stochastic sampling errors, using the single-mode anharmonic oscillator as an example. In the modeling of this system, which is well known for being problematic, we achieve an improvement of more than 20 orders of magnitude in the distribution of the trajectories.