Porosity and surface area analysis play a prominent role in modern materials science. At the heart of this sits the Brunauer–Emmett–Teller (BET) theory, which has been a remarkably successful contribution to the field of materials science. The BET method was developed in the 1930s for open surfaces but is now the most widely used metric for the estimation of surface areas of micro‐ and mesoporous materials. Despite its widespread use, the calculation of BET surface areas causes a spread in reported areas, resulting in reproducibility problems in both academia and industry. To prove this, for this analysis, 18 already‐measured raw adsorption isotherms were provided to sixty‐one labs, who were asked to calculate the corresponding BET areas. This round‐robin exercise resulted in a wide range of values. Here, the reproducibility of BET area determination from identical isotherms is demonstrated to be a largely ignored issue, raising critical concerns over the reliability of reported BET areas. To solve this major issue, a new computational approach to accurately and systematically determine the BET area of nanoporous materials is developed. The software, called “BET surface identification” (BETSI), expands on the well‐known Rouquerol criteria and makes an unambiguous BET area assignment possible.

Thermal engineering of metal–organic frameworks for adsorption-based applications is very topical in view of their industrial potential, in particular, since heat management and thermal stability have been identified as important obstacles. Hence, a fundamental understanding of the structural and chemical features underpinning their intrinsic thermal properties is highly sought-after. Herein, we investigate the nanoscale behavior of a diverse set of frameworks using molecular simulation techniques and critically compare properties such as thermal conductivity, heat capacity, and thermal expansion with other classes of materials. Furthermore, we propose a hypothetical thermodynamic cycle to estimate the temperature rise associated with adsorption for the most important greenhouse and energy-related gases (CO2 and CH4). This macroscopic response on the heat of adsorption connects the intrinsic thermal properties with the adsorption properties and allows us to evaluate their importance.

Optimal mechanical impact absorbers are reusable and exhibit high specific energy absorption. The forced intrusion of liquid water in hydrophobic nanoporous materials, such as zeolitic imidazolate frameworks (ZIFs), presents an attractive pathway to engineer such systems. However, to harness their full potential, it is crucial to understand the underlying water intrusion and extrusion mechanisms under realistic, high-rate deformation conditions. Herein, we report a critical increase of the energy absorption capacity of confined water-ZIF systems at elevated strain rates. Starting from ZIF-8 as proof-of-concept, we demonstrate that this attractive rate dependence is generally applicable to cage-type ZIFs but disappears for channel-containing zeolites. Molecular simulations reveal that this phenomenon originates from the intrinsic nanosecond timescale needed for critical-sized water clusters to nucleate inside the nanocages, expediting water transport through the framework. Harnessing this fundamental understanding, design rules are formulated to construct effective, tailorable, and reusable impact energy absorbers for challenging new applications.

Metal-organic frameworks show both fundamental interest and great promise for applications in adsorption-based technologies, such as the separation and storage of gases. The flexibility and complexity of the molecular scaffold poses a considerable challenge to atomistic modeling, especially when also considering the presence of guest molecules. We investigate the role played by quantum and anharmonic fluctuations in the archetypical case of MOF-5, comparing the material at various levels of methane loading. Accurate path integral simulations of such effects are made affordable by the introduction of an accelerated simulation scheme and the use of an optimized force 1 arXiv:1901.03770v1 [cond-mat.mtrl-sci] 11 Jan 2019field based on first-principles reference calculations. We find that the level of statistical treatment that is required for predictive modeling depends significantly on the property of interest. The thermal properties of the lattice are generally well described by a quantum harmonic treatment, with the adsorbate behaving in a classical but strongly anharmonic manner. The heat capacity of the loaded framework -which plays an important role in the characterization of the framework and in determining its stability to thermal fluctuations during adsorption/desorption cycles -requires, however, a full quantum and anharmonic treatment, either by path integral methods or by a simple but approximate scheme. We also present molecular-level insight into the nanoscopic interactions contributing to the material's properties and suggest design principles to optimize them.

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.