The wavefunctions in eq (27) are given as the product of a Gaussian and a Hermite polynomial. One can also sample coordinates and momenta according to a Hermite function whose Gaussian exponent and center are provided by the user. 69 This feature is helpful for accelerating the convergence of certain observables as a function of the number of trajectories 70 and transforming the results from a given PDF to any desired PDF, 69 in both cases, using the importance sampling technique.