Bridging molecular dynamics and correlated wave-function methods for accurate finite-temperature properties