The sampling procedure is simple. First, we choose a reasonable τ . Next, for each coordinate (x, y, z) of all nuclei we generate a displacement and velocity according to the distribution given in Equation . The procedure can also be defined in terms of a multivariate quasi-distribution of all Cartesian coordinates and momentums, similar to that in Equation . In this case, the sampling is done on the basis of the distribution