Bootstrapping the stationary state of bosonic open quantum systems

Abstract

We propose a method to compute expectation values of observables in the stationary state of a (Markovian) bosonic open quantum system. Using a hierarchy of semi-definite relaxations, we obtain finer and finer upper and lower bounds to any expectation value of interest. The bounds are rigorous, robust to stationary state degeneracies, and numerically improve as the occupation number increases on the examples we considered. This makes it adapted to the simulation of stationary states of bosonic qubits and in particular dissipatively stabilized cat qubits.