How should we calculate multi-dimensional potential energy surfaces for an accurate reproduction of partition functions?