Multimode Bogoliubov transformation and Husimi’s Q-function

In this paper, we present numerical schemes for evaluating the matrix elements of Gaussian/non-Gaussian operators in the Fock state basis, which are identified as multivariate Hermite polynomials (MHPs). Using the integral transformation operator to perform the multimode Bogoliubov transformation, Husimi’s Q-functions of Gaussian/non-Gaussian operators are easily derived as the generating functions of MHPs.