The PESs obtained from the MuFi approach is input to the N ×1DHR statistical mechanics treatment. We model the 1DHRs using the python package TAMkin. At least 100 basis functions were used to model the wave functions. A higher number not only gives enough flexibility to the wave functions but also leads to a high number of resulting energy levels.