In the final scheme, 10,000 configurations from the T = 500 K simulation were optimized at the B3LYP/cc-pVDZ level of theory and vibrational frequencies were computed at the optimized geometry. Configurations were taken from each step in the optimization procedure if the energy was more than 0.25 kcal/mol above the energy of the fully optimized structure.