The approach we adopt most closely follows the method of Refs. 63 and 104. Similar to some of the other methods listed above, configurations are statistically sampled from a Gaussian probability distribution defined by the quantum covariance of atomic displacements. Quantum covariance of a given pair of atoms a and b in their respective Cartesian directions i and j :