Center for Molecular Modeling - S. Vandenbrande

High-rate nanofluidic energy absorption in porous zeolitic frameworks

sven — Thu, 22 Apr 2021 19:37:50 +0000

Y. Sun, S.M.J. Rogge, A. Lamaire, S. Vandenbrande, J. Wieme, C.R. Siviour, V. Van Speybroeck, J.-C. Tan

Nature Materials

20 (7), 1015–1023

2021

Abstract

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. Here, 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.

DOI

http://dx.doi.org/10.1038/s41563-021-00977-6

Private attachment

s41563-021-00977-6.pdf

Chlorination of a Zeolitic-Imidazolate Framework Tunes Packing and van der Waals Interaction of Carbon Dioxide for Optimized Adsorptive Separation

sven — Sun, 07 Feb 2021 12:54:58 +0000

L.H. Wee, S. Vandenbrande, S.M.J. Rogge, J. Wieme, K. Asselman, E. Jardim, J. Silvestre-Albero, J. Navarro, V. Van Speybroeck, J.A. Martens, C. Kirschhock

JACS (Journal of the American Chemical Society)

143 (13), 4962-4968

2021

Abstract

Molecular separation of carbon dioxide (CO₂) and methane (CH₄) is of growing interest for biogas upgrading, carbon capture and utilization, methane synthesis and for purification of natural gas. Here, we report a new zeolitic-imidazolate framework (ZIF), coined COK-17, with exceptionally high affinity for the adsorption of CO₂ by London dispersion forces, mediated by chlorine substituents of the imidazolate linkers. COK-17 is a new type of flexible zeolitic-imidazolate framework Zn(4,5-dichloroimidazolate)₂ with the SOD framework topology. Below 200 K it displays a metastable closed-pore phase next to its stable open-pore phase. At temperatures above 200 K, COK-17 always adopts its open-pore structure, providing unique adsorption sites for selective CO₂ adsorption and packing through van der Waals interactions with the chlorine groups, lining the walls of the micropores. Localization of the adsorbed CO₂ molecules by Rietveld refinement of X-ray diffraction data and periodic density functional theory calculations revealed the presence and nature of different adsorption sites. In agreement with experimental data, grand canonical Monte Carlo simulations of adsorption isotherms of CO₂ and CH₄ in COK-17 confirmed the role of the chlorine functions of the linkers and demonstrated the superiority of COK-17 compared to other adsorbents such as ZIF-8 and ZIF-71.

Gold Open Access

DOI

http://dx.doi.org/10.1021/jacs.0c08942

Private attachment

jacs.0c08942.pdf

Charting the Metal-Dependent High-Pressure Stability of Bimetallic UiO-66 Materials

sven — Fri, 20 Mar 2020 17:16:34 +0000

S.M.J. Rogge, P.G. Yot, J. Jacobsen, F. Muniz-Miranda, S. Vandenbrande, J. Gosch, V. Ortiz, I. Collings, S. Devautour-Vinot, G. Maurin, N. Stock, V. Van Speybroeck

ACS Materials Letters

2 (4), 438-445

2020

Abstract

In theory, bimetallic UiO-66(Zr:Ce) and UiO-66(Zr:Hf) metal-organic frameworks (MOFs) are extremely versatile and attractive nanoporous materials as they combine the high catalytic activity of UiO-66(Ce) or UiO-66(Hf) with the outstanding stability of UiO-66(Zr). Using in situ high-pressure powder X-ray diffraction, however, we observe that this expected mechanical stability is not achieved when incorporating cerium or hafnium in UiO-66(Zr). This observation is akin to the earlier observed reduced thermal stability of UiO-66(Zr:Ce) compounds. To elucidate the atomic origin of this phenomenon, we chart the loss-of-crystallinity pressures of 22 monometallic and bimetallic UiO-66 materials and systematically isolate their intrinsic mechanical stability from their defect-induced weakening. This complementary experimental/computational approach reveals that the intrinsic mechanical stability of these bimetallic MOFs decreases nonlinearly upon cerium incorporation but remains unaffected by the zirconium:hafnium ratio. Additionally, all experimental samples suffer from defect-induced weakening, a synthesis-controlled effect that is observed to be independent of their intrinsic stability.

Gold Open Access

DOI

http://dx.doi.org/10.1021/acsmaterialslett.0c00042

Private attachment

acsmaterialslett.0c00042.pdf

Light Olefin Diffusion during the MTO Process on H-SAPO-34: a Complex Interplay of Molecular Factors

cnudde — Sun, 15 Mar 2020 16:53:30 +0000

P. Cnudde, R. Demuynck, S. Vandenbrande, M. Waroquier, G. Sastre, V. Van Speybroeck

JACS (Journal of the American Chemical Society)

142 (13), 6007-6017

2020

Abstract

The methanol-to-olefins process over H-SAPO-34 is characterized by its high shape selectivity toward light olefins. The catalyst is a supramolecular system consisting of nanometer-sized inorganic cages, decorated by Brønsted acid sites, in which organic compounds, mostly methylated benzene species, are trapped. These hydrocarbon pool species are essential to catalyze the methanol conversion but may also clog the pores. As such, diffusion of ethene and propene plays an essential role in determining the ultimate product selectivity. Enhanced sampling molecular dynamics simulations based on either force fields or density functional theory are used to determine how molecular factors influence the diffusion of light olefins through the 8-ring windows of H-SAPO-34. Our simulations show that diffusion through the 8-ring in general is a hindered process, corresponding to a hopping event of the diffusing molecule between neighboring cages. The loading of different methanol, alkene, and aromatic species in the cages may substantially slow down or facilitate the diffusion process. The presence of Brønsted acid sites in the 8-ring enhances the diffusion process due to the formation of a favorable π-complex host–guest interaction. Aromatic hydrocarbon pool species severely hinder the diffusion and their spatial distribution in the zeolite crystal may have a significant effect on the product selectivity. Herein, we unveil how molecular factors influence the diffusion of light olefins in a complex environment with confined hydrocarbon pool species, high olefin loadings, and the presence of acid sites by means of enhanced molecular dynamics simulations under operating conditions.

DOI

http://dx.doi.org/10.1021/jacs.9b10249

Private attachment

Cnudde_Diffusion_JACS.pdf

Thermal Engineering of Metal-Organic Frameworks for Adsorption Applications: A Molecular Simulations Perspective

jelle — Thu, 01 Aug 2019 08:17:37 +0000

J. Wieme, S. Vandenbrande, A. Lamaire, V. Kapil, L. Vanduyfhuys, V. Van Speybroeck

ACS Applied Materials & Interfaces

11 (42), 38697-38707

2019

Abstract

Thermal engineering of metal-organic frameworks (MOFs) for adsorption-based applications is very topical in view of their industrial potential, especially 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 material classes. 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.

DOI

http://dx.doi.org/10.1021/acsami.9b12533

Private attachment

Wieme et al., 2019.pdf

Modeling the structural and thermal properties of loaded metal-organic frameworks. An interplay of quantum and anharmonic fluctuations

jelle — Fri, 21 Dec 2018 13:33:25 +0000

V. Kapil, J. Wieme, S. Vandenbrande, A. Lamaire, V. Van Speybroeck, M. Ceriotti

Journal of Chemical Theory and Computation

15 (5), 3237-3249

2019

DOI

http://dx.doi.org/10.1021/acs.jctc.8b01297

Private attachment

Kapil_JChemTheoryComput2019.pdf

Modeling Gas Adsorption in Flexible Metal–Organic Frameworks via Hybrid Monte Carlo / Molecular Dynamics Schemes

sven — Fri, 21 Dec 2018 13:00:52 +0000

S.M.J. Rogge, R. Goeminne, R. Demuynck, J.J. Gutiérrez-Sevillano, S. Vandenbrande, L. Vanduyfhuys, M. Waroquier, T. Verstraelen, V. Van Speybroeck

Advanced Theory and Simulations

2 (4), 1800177

2019

Abstract

Herein, a hybrid Monte Carlo (MC)/molecular dynamics (MD) simulation protocol that properly accounts for the extraordinary structural flexibility of metal–organic frameworks (MOFs) is developed and validated. This is vital to accurately predict gas adsorption isotherms and guest‐induced flexibility of these materials. First, the performance of three recent models to predict adsorption isotherms and flexibility in MOFs is critically investigated. While these methods succeed in providing qualitative insight in the gas adsorption process in MOFs, their accuracy remains limited as the intrinsic flexibility of these materials is very hard to account for. To overcome this challenge, a hybrid MC/MD simulation protocol that is specifically designed to handle the flexibility of the adsorbent, including the shape flexibility, is introduced, thereby unifying the strengths of the previous models. It is demonstrated that the application of this new protocol to the adsorption of neon, argon, xenon, methane, and carbon dioxide in MIL‐53(Al), a prototypical flexible MOF, substantially decreases the inaccuracy of the obtained adsorption isotherms and predicted guest‐induced flexibility. As a result, this method is ideally suited to rationalize the adsorption performance of flexible nanoporous materials at the molecular level, paving the way for the conscious design of MOFs as industrial adsorbents.

Gold Open Access

DOI

http://dx.doi.org/10.1002/adts.201800177

Private attachment

main.pdf

Ab initio evaluation of Henry coefficients using importance sampling

steven — Wed, 31 Oct 2018 08:34:51 +0000

S. Vandenbrande, M. Waroquier, V. Van Speybroeck, T. Verstraelen

Journal of Chemical Theory and Computation

14 (12), 6359–6369

2018

Abstract

We present a new algorithm that allows for an efficient evaluation of the Henry coefficient of a guest molecule inside a porous material, which permits to use ab initio energy calculations. The Widom insertion method, which is currently used to compute these Henry coefficients, typically requires millions of energy evaluations. Our new methodology reduces this number by more than 1 order of magnitude, enabling the use of an ab initio potential energy surface. The methodology we propose is reminiscent of the well-known importance sampling technique which is frequently used in Monte Carlo integrations. First, a conventional Widom insertion simulation is performed using a force field. In the second step, the Widom results are used to select a limited number of configurations and only for these configurations the ab initio evaluation of the energy is required. Finally, by appropriately reweighting the latter energies, an accurate estimation of the ab initio Henry coefficient is possible at a moderate computational cost. We apply our methodology to the adsorption of CO2 in Mg-MOF-74, a prototypical case where interactions of a polar guest molecule with unsaturated metal sites dominate the adsorption mechanism. In this case generic force fields such as UFF or Dreiding are inappropriate and the use of ab initio methods is indispensable. In a second case study, we compute Henry coefficients of methane in UiO-66 using different levels of theory. We pay particular attention to the influence of the dispersion corrections and the role of many-body effects. For the final example, we qualitatively investigate adsorption features for a series of functionalized UiO-66 frameworks. Overall the cases we present show that accurate computations of Henry coefficients is extremely challenging, as different levels of theory provide strongly varying results. At the same time ab initio calculations have added value compared to force fields, as they provide a physically more sound description of the adsorption mechanism and in some cases clearly improve correspondence with experiment.

Open Access version available at UGent repository

Gold Open Access

DOI

http://dx.doi.org/10.1021/acs.jctc.8b00892

i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations

jelle — Sun, 24 Jun 2018 18:17:00 +0000

V. Kapil, M. Rossi, O. Marsalek, R. Petraglia, Y. Litman, T. Spura, B. Cheng, A. Cuzzocrea, R.H. Meißner, D. Wilkins, P. Juda, S.P. Bienvenue, J. Kessler, I. Poltavsky, S. Vandenbrande, J. Wieme, C. Corminboeuf, T. Kühne, D. Manolopoulos, T. Markland, J. Richardson, A. Tkatchenko, G.A. Tribello, V. Van Speybroeck, M. Ceriotti

Computer Physics Communications

236, 214-223

2019

Abstract

Progress in the atomic-scale modelling of matter over the past decade has been tremendous. This progress has been brought about by improvements in methods for evaluating interatomic forces that work by either solving the electronic structure problem explicitly, or by computing accurate approximations of the solution and by the development of techniques that use the Born-Oppenheimer (BO) forces to move the atoms on the BO potential energy surface. As a consequence of these developments it is now possible to identify stable or metastable states, to sample configurations consistent with the appropriate thermodynamic ensemble, and to estimate the kinetics of reactions and phase transitions. All too often, however, progress is slowed down by the bottleneck associated with implementing new optimization algorithms and/or sampling techniques into the many existing electronic-structure and empirical-potential codes. To address this problem, we are thus releasing a new version of the i-PI software. This piece of software is an easily extensible framework for implementing advanced atomistic simulation techniques using interatomic potentials and forces calculated by an external driver code. While the original version of the code was developed with a focus on path integral molecular dynamics techniques, this second release of i-PI not only includes several new advanced path integral methods, but also offers other classes of algorithms. In other words, i-PI is moving towards becoming a universal force engine that is both modular and tightly coupled to the driver codes that evaluate the potential energy surface and its derivatives.

http://arxiv.org/abs/1808.03824

Open Access version available at UGent repository

Green Open Access

DOI

http://dx.doi.org/10.1016/j.cpc.2018.09.020

Private attachment

1808.03824.pdf

Extension of the QuickFF force field protocol for an improved accuracy of structural, vibrational, mechanical and thermal properties of Metal Organic Frameworks

wim — Mon, 27 Nov 2017 07:12:54 +0000

L. Vanduyfhuys, S. Vandenbrande, J. Wieme, M. Waroquier, T. Verstraelen, V. Van Speybroeck

Journal of Computational Chemistry

39 (16), p. 999-1011

2018

Abstract

QuickFF was originally launched in 2015 to derive accurate force fields for isolated and complex molecular systems in a quick and easy way. Apart from the general applicability, the functionality was especially tested for metal-organic frameworks (MOFs), a class of hybrid materials consisting of organic and inorganic building blocks. Herein, we launch a new release of the QuickFF protocol which includes new major features to predict structural, vibrational, mechanical and thermal properties with greater accuracy, without compromising its robustness and transparant workflow. First, the ab initio data necessary for the fitting procedure may now also be derived from periodic models for the molecular system, as opposed to the earlier cluster-based models. This is essential for an accurate description of MOFs with one dimensional metal-oxide chains. Second, cross terms that couple internal coordinates (ICs) and anharmonic contributions for bond and bend terms are implemented. These features are essential for a proper description of vibrational and thermal properties. Third, the fitting scheme was modified to improve robustness and accuracy. The new features are tested on MIL-53(Al), MOF-5, CAU-13 and NOTT-300. As expected, periodic input data is proven to be essential for a correct description of structural, vibrational and thermodynamic properties of MIL-53(Al). Bulk moduli and thermal expansion coefficients of MOF-5 are very accurately reproduced by static and dynamic simulations using the newly derived force fields which include cross terms and anharmonic corrections. For the flexible materials CAU-13 and NOTT-300, the transition pressure is accurately
predicted provided cross terms are taken into account.

Open Access version available at UGent repository

Gold Open Access

DOI

http://dx.doi.org/10.1002/jcc.25173

Private attachment

published_online.pdf