Abstract
In this article, we investigate the influence of anharmonicities and nuclear quantum effects (NQEs) in modelling the structural properties and thermal expansion of the empty MOF-5 metal-organic framework. To introduce NQEs in classical molecular dynamics simulations, two different methodologies are considered, comparing the approximate, but computationally cheap, method of generalised Langevin equation thermostatting to the more advanced, computationally demanding path integral molecular dynamics technique. For both methodologies, similar results were obtained for all the properties under investigation. The structural properties of MOF-5, probed by means of radial distribution functions (RDFs), show some distinct differences with respect to a classical description. Besides a broadening of the RDF peaks under the influence of quantum fluctuations, a different temperature dependence is also observed due to a dominant zero-point energy (ZPE) contribution. For the thermal expansion of MOF-5, by contrast, NQEs appear to be only of secondary importance with respect to an adequate modelling of the anharmonicities of the potential energy surface (PES), as demonstrated by the use of two differently parametrised force fields. Despite the small effect in the temperature dependence of the volume of MOF-5, NQEs do however significantly affect the absolute volume of MOF-5, in which the ZPE resulting from the intertwining of NQEs and anharmonicities plays a crucial role. A sufficiently accurate description of the PES is therefore prerequisite when modelling NQEs.