In this talk, I will discuss how to enable such predictions by combining advanced statistical mechanics with data-driven machine learning interatomic potentials. As an example [1], for the omnipresent and technologically essential system of water, a first-principles thermodynamic description not only leads to excellent agreement with experiments, but also reveals the crucial role of nuclear quantum fluctuations in modulating the thermodynamic stabilities of different phases of water. As another example [2], we simulated the high pressure hydrogen system with converged system size and simulation length, and found, contrary to established beliefs, supercritical behaviour of liquid hydrogen above the melting line. Finally, I will talk about transport properties. Ref [3] proposed a method to compute the heat conductivities of liquid just from equilibrium molecular dynamics trajectories.