  Intense terahertz-frequency pulses induce coherent rotational dynamics and
orientation of polar molecular ensembles. Exact numerical methods for
rotational dynamics are computationally not feasible for the vast majority of
molecular rotors - the asymmetric top molecules at ambient temperatures. We
exemplify the use of Random Phase Wave Functions (RPWF) by calculating the
terahertz-induced rotational dynamics of sulfur dioxide (SO2) at ambient
temperatures and high field strengths and show that the RPWF method gains
efficiency with the increase in temperature and in the THz-field strengths. The
presented method provides wide-ranging computational access to rotational
dynamical responses of molecules at experimental conditions which are far
beyond the reach of exact numerical methods.
