where k B is the Boltzmann constant. This generation of the initial conditions is very effective for classical ensembles. However, we can also use MBS for crude simulation of NQEs. In simulations at elevated temperatures, pair distribution functions (PDFs) are similar to those in the classical case. The squared vibrational amplitude (Equation ) for the interatomic distance r in the classical harmonic approximation is: