Variational Multi-Gaussian Phase-Space Dynamics via Automatic Differentiation
This method simulates interacting open quantum bosonic systems using a multi-dimensional Wigner phasespace representation and a Variational Multi-Gaussian (VMG) ansatz. It successfully reveals critical slowing down with 2D quantum Ising universality class exponents in driven-dissipative systems, offering scalable accuracy comparable to exact diagonalization.