Quantum simulation of large open systems is a hard task that demands huge computation and memory costs.

The rotational dynamics of non-linear molecules at high-temperature under external fields is such an example. For room temperatures, the initial density matrix populates ~ 10⁴ rotational states, and the whole coupled Hilbert space can reach ~ 10⁶ states. Simulation by neither the direct density matrix nor the full basis set of populated wavefunctions is impossible.

We employ the random phase wave function method to represent the initial state and compute several time dependent and independent observables such as the orientation and the alignment of the molecules. The error of the method was found to scale as N^-1/2, where N is the number of wave function realizations employed. Scaling vs. the temperature was computed for weak and strong fields. As expected, the convergence of the method increase rapidly with the temperature and the field intensity.