Gaussian dynamics equation in normal product form

In this paper, we discuss the normal product form of the density operator of multimode Gaussian states, and obtain the correlation equation between the kernel matrix R of the Gaussian density operator in the normal product form and its kernel matrix G in the standard quadratic form. Further, we explore the time evolution mechanism of R and obtain the Gaussian dynamical equation under the normal product R=i(RJH-HJR). Our work is devoted to searching for another mechanism for Gaussian dynamics. By exploring the description of the normal ordered density matrix under the coherent state representation, we find that our mechanism is feasible and easy to operate.