M. Waroquier

Efficient Construction of Free Energy Profiles of Breathing Metal-Organic Frameworks Using Advanced Molecular Dynamics Simulations

R. Demuynck, S.M.J. Rogge, L. Vanduyfhuys, J. Wieme, M. Waroquier, V. Van Speybroeck
Journal of Chemical Theory and Computation (JCTC)
13 (12), 5861-5873
2017
A1

Abstract 

In order to reliably predict and understand the breathing behavior of highly flexible metal–organic frameworks from thermodynamic considerations, an accurate estimation of the free energy difference between their different metastable states is a prerequisite. Herein, a variety of free energy estimation methods are thoroughly tested for their ability to construct the free energy profile as a function of the unit cell volume of MIL-53(Al). The methods comprise free energy perturbation, thermodynamic integration, umbrella sampling, metadynamics, and variationally enhanced sampling. A series of molecular dynamics simulations have been performed in the frame of each of the five methods to describe structural transformations in flexible materials with the volume as the collective variable, which offers a unique opportunity to assess their computational efficiency. Subsequently, the most efficient method, umbrella sampling, is used to construct an accurate free energy profile at different temperatures for MIL-53(Al) from first principles at the PBE+D3(BJ) level of theory. This study yields insight into the importance of the different aspects such as entropy contributions and anharmonic contributions on the resulting free energy profile. As such, this thorough study provides unparalleled insight in the thermodynamics of the large structural deformations of flexible materials.

Open Access version available at UGent repository
Gold Open Access

Reliably Modeling the Mechanical Stability of Rigid and Flexible Metal-Organic Frameworks

S.M.J. Rogge, M. Waroquier, V. Van Speybroeck
Accounts of Chemical Research
51 (1), 138-148
2018
A1

Abstract 

Over the past two decades, metal–organic frameworks (MOFs) have matured from interesting academic peculiarities toward a continuously expanding class of hybrid, nanoporous materials tuned for targeted technological applications such as gas storage and heterogeneous catalysis. These oft-times crystalline materials, composed of inorganic moieties interconnected by organic ligands, can be endowed with desired structural and chemical features by judiciously functionalizing or substituting these building blocks. As a result of this reticular synthesis, MOF research is situated at the intriguing intersection between chemistry and physics, and the building block approach could pave the way toward the construction of an almost infinite number of possible crystalline structures, provided that they exhibit stability under the desired operational conditions. However, this enormous potential is largely untapped to date, as MOFs have not yet found a major breakthrough in technological applications. One of the remaining challenges for this scale-up is the densification of MOF powders, which is generally achieved by subjecting the material to a pressurization step. However, application of an external pressure may substantially alter the chemical and physical properties of the material. A reliable theoretical guidance that can presynthetically identify the most stable materials could help overcome this technological challenge.

In this Account, we describe the recent research the progress on computational characterization of the mechanical stability of MOFs. So far, three complementary approaches have been proposed, focusing on different aspects of mechanical stability: (i) the Born stability criteria, (ii) the anisotropy in mechanical moduli such as the Young and shear moduli, and (iii) the pressure-versus-volume equations of state. As these three methods are grounded in distinct computational approaches, it is expected that their accuracy and efficiency will vary. To date, however, it is unclear which set of properties are suited and reliable for a given application, as a comprehensive comparison for a broad variety of MOFs is absent, impeding the widespread use of these theoretical frameworks.

Herein, we fill this gap by critically assessing the performance of the three computational models on a broad set of MOFs that are representative for current applications. These materials encompass the mechanically rigid UiO-66(Zr) and MOF-5(Zn) as well as the flexible MIL-47(V) and MIL-53(Al), which undergo pressure-induced phase transitions. It is observed that the Born stability criteria and pressure-versus-volume equations of state give complementary insight into the macroscopic and microscopic origins of instability, respectively. However, interpretation of the Born stability criteria becomes increasingly difficult when less symmetric materials are considered. Moreover, pressure fluctuations during the simulations hamper their accuracy for flexible materials. In contrast, the pressure-versus-volume equations of state are determined in a thermodynamic ensemble specifically targeted to mitigate the effects of these instantaneous fluctuations, yielding more accurate results. The critical Account presented here paves the way toward a solid computational framework for an extensive presynthetic screening of MOFs to select those that are mechanically stable and can be postsynthetically densified before their use in targeted applications.

Open Access version available at UGent repository
Gold Open Access

Influence of a confined methanol solvent on the reactivity of active sites in UiO-66

C. Caratelli, J. Hajek, S.M.J. Rogge, S. Vandenbrande, E.J. Meijer, M. Waroquier, V. Van Speybroeck
ChemPhysChem
19 (4), 420-429
2018
A1

Abstract 

UiO-66, composed of Zr-oxide bricks and terephthalate linkers, is currently one of the most studied metal-organic frameworks due to its exceptional stability. Defects can be introduced in the structure, creating undercoordinated Zr atoms which are Lewis acid sites. Here, additional Brønsted sites can be generated by coordinated protic species from the solvent. In this contribution, a multilevel modeling approach was applied to unravel the effect of a confined methanol solvent on the active sites in UiO-66. First, active sites were explored with static periodic density functional theory calculations to investigate adsorption of water and methanol. Solvent was then introduced in the pores with grand canonical Monte Carlo simulations, followed by a series of molecular dynamics simulations at operating conditions. A hydrogen-bonded network of methanol molecules is formed, allowing the protons to shuttle between solvent methanol, adsorbed water, and the inorganic brick. Upon deprotonation of an active site, the methanol solvent aids the transfer of protons and stabilizes charged configurations via hydrogen bonding, which could be crucial in stabilizing reactive intermediates. The multilevel modeling approach adopted here sheds light on the important role of a confined solvent on the active sites in the UiO-66 material, introducing dynamic acidity in the system at finite temperatures by which protons may be easily shuttled from various positions at the active sites.

Open Access version available at UGent repository
Gold Open Access

Methane Adsorption in Zr-Based MOFs: Comparison and Critical Evaluation of Force Fields

S. Vandenbrande, T. Verstraelen, J. J. Gutierrez-Sevillano, M. Waroquier, V. Van Speybroeck
Journal of Physical Chemistry C
121 (45), 25309-25322
2017
A1

Abstract 

The search for nanoporous materials that are highly performing for gas storage and separation is one of the contemporary challenges in material design. The computational tools to aid these experimental efforts are widely available and adsorption isotherms are routinely computed for huge sets of (hypothetical) frameworks. Clearly the computational results depend on the interactions between the adsorbed species and the adsorbent, which are commonly described using force fields. In this paper, an extensive comparison and in-depth investigation of several force fields from literature is reported for the case of methane adsorption in the Zr-based Metal-Organic Frameworks UiO-66, UiO-67, DUT-52, NU-1000 and MOF-808. Significant quantitative differences in the computed uptake are observed when comparing different force fields, but most qualitative features are common which suggests some predictive power of the simulations when it comes to these properties. More insight into to the host-guest interactions is obtained by benchmarking the force fields with an extensive number of ab initio computed single molecule interaction energies. This analysis at the molecular level reveals that especially ab initio derived force fields perform well in reproducing the ab initio interaction energies. Finally, the high sensitivity of uptake predictions on the underlying potential energy surface is explored.

Open Access version available at UGent repository
Gold Open Access

Nature of active sites on UiO-66 and beneficial influence of water in the catalysis of Fischer esterification

C. Caratelli, J. Hajek, F. G. Cirujano, M. Waroquier, F. X. Llabres i Xamena, V. Van Speybroeck
Journal of Catalysis
352, 401-414
2017
A1

Abstract 

Zirconium terephthalate UiO-66 type metal organic frameworks (MOFs) are known to be highly active, stable and reusable catalysts for the esterification of carboxylic acids with alcohols. Moreover, when defects are present in the structure of these MOFs, coordinatively unsaturated Zr ions with Lewis acid properties are created, which increase the catalytic activity of the resulting defective solids. In the present work, molecular modeling techniques combined with new experimental data on various defective hydrated and dehydrated materials allow to unravel the nature and role of defective active sites in the Fischer esterification and the role of coordinated water molecules to provide additional Brønsted sites. Periodic models of UiO-66 and UiO-66-NH2 catalysts have been used to unravel the reaction mechanism on hydrated and dehydrated materials. Various adsorption modes of water and methanol are investigated. The proposed mechanisms are in line with experimental observations that amino groups yield a reduction in the reaction barriers, although they have a passive role in modulating the electronic structure of the material. Water has a beneficial role on the reaction cycle by providing extra Brønsted sites and by providing stabilization for various intermediates through hydrogen bonds.

Open Access version available at UGent repository
Gold Open Access

The remarkable amphoteric nature of defective UiO-66 in catalytic reactions

J. Hajek, B. Bueken, M. Waroquier, D. De Vos, V. Van Speybroeck
ChemCatChem
9 (12), 2203-2210
2017
A1

Abstract 

One of the major requirements in solid acids and bases catalyzed reactions is that the reactants, intermediates or activated complexes cooperate with several functions of catalyst support. In this work the remarkable bifunctional behavior of the defective UiO-66(Zr) metal organic framework is shown for acid-base pair catalysis. The active site relies on the presence of undercoordinated zirconium sites, which may be tuned by removing framework linkers and by removal of water from the inorganic bricks using a dehydration treatment. To elucidate the amphoteric nature of defective UiO-66, the Oppenauer oxidation of primary alcohols has been theoretically investigated using density functional theory (DFT) and the periodic approach. The presence of acid and basic centers within molecular distances has been shown crucial for determining the catalytic activity of the material. Hydrated and dehydrated bricks have a distinct influence on modulation of the acidity and basicity of the active sites. In any case both functions need to cooperate in a concerted way to enable the chemical transformation.

Open Access version available at UGent repository
Gold Open Access

The Monomer Electron Density Force Field (MEDFF): A Physically Inspired Model for Non-Covalent Interactions

S. Vandenbrande, M. Waroquier, V. Van Speybroeck, T. Verstraelen
Journal of Chemical Theory and Computation (JCTC)
13 (1), 161–179
2017
A1

Abstract 

We propose a methodology to derive pairwise-additive noncovalent force fields from monomer electron densities without any empirical input. Energy expressions are based on the symmetry-adapted perturbation theory (SAPT) decomposition of interaction energies. This ensures a physically motivated force field featuring an electrostatic, exchange repulsion, dispersion, and induction contribution, which contains two types of parameters. First, each contribution depends on several fixed atomic parameters, resulting from a partitioning of the monomer electron density. Second, each of the last three contributions (exchange-repulsion, dispersion, and induction) contains exactly one linear fitting parameter. These three so-called interaction parameters in the model are initially estimated separately using SAPT reference calculations for the S66x8 database of noncovalent dimers. In a second step, the three interaction parameters are further refined simultaneously to reproduce CCSD(T)/CBS interaction energies for the same database. The limited number of parameters that are fitted to dimer interaction energies (only three) avoids ill-conditioned fits that plague conventional parameter optimizations. For the exchange repulsion and dispersion component, good results are obtained for all dimers in the S66x8 database using one single value for the associated interaction parameters. The values of those parameters can be considered universal and can also be used for dimers not present in the original database used for fitting. For the induction component such an approach is only viable for the dispersion dominated dimers in the S66x8 database. For other dimers (such as hydrogen-bonded complexes), we show that our methodology remains applicable. However, the interaction parameter needs to be determined on a case-specific basis. As an external validation:, the force field predicts interaction energies in good agreement with CCSD(T)/CBS values for dispersion dominated dimers extracted from an HIV-II protease crystal structure with a bound ligand (indinavir). Furthermore, experimental second virial coefficients of small alkanes and alkenes are well reproduced.

Open Access version available at UGent repository
Green Open Access

Thermodynamic Insight in the High-Pressure Behavior of UiO-66: Effect of Linker Defects and Linker Expansion

S.M.J. Rogge, J. Wieme, L. Vanduyfhuys, S. Vandenbrande, G. Maurin, T. Verstraelen, M. Waroquier, V. Van Speybroeck
Chemistry of Materials
28 (16), 5721-5732
2016
A1

Abstract 

In this Article, we present a molecular-level understanding of the experimentally observed loss of crystallinity in UiO-66-type metal–organic frameworks, including the pristine UiO-66 to -68 as well as defect-containing UiO-66 materials, under the influence of external pressure. This goal is achieved by constructing pressure-versus-volume profiles at finite temperatures using a thermodynamic approach relying on ab initio derived force fields. On the atomic level, the phenomenon is reflected in a sudden drop in the number of symmetry operators for the crystallographic unit cell because of the disordered displacement of the organic linkers with respect to the inorganic bricks. For the defect-containing samples, a reduced mechanical stability is observed, however, critically depending on the distribution of these defects throughout the material, hence demonstrating the importance of judiciously characterizing defects in these materials.

This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.
http://pubs.acs.org/doi/abs/10.1021/acs.chemmater.6b01956

Open Access version available at UGent repository
Gold Open Access

Effect of temperature and branching on the nature and stability of alkene cracking intermediates in H-ZSM-5

P. Cnudde, K. De Wispelaere, J. Van der Mynsbrugge, M. Waroquier, V. Van Speybroeck
Journal of Catalysis
345, 53-69
2017
A1

Abstract 

Catalytic cracking of alkenes takes place at elevated temperatures in the order of 773–833 K. In this work, the nature of the reactive intermediates at typical reaction conditions is studied in H-ZSM-5 using a complementary set of modeling tools. Ab initio static and molecular dynamics simulations are performed on different C4single bond C5 alkene cracking intermediates to identify the reactive species in terms of temperature. At 323 K, the prevalent intermediates are linear alkoxides, alkene π-complexes and tertiary carbenium ions. At a typical cracking temperature of 773 K, however, both secondary and tertiary alkoxides are unlikely to exist in the zeolite channels. Instead, more stable carbenium ion intermediates are found. Branched tertiary carbenium ions are very stable, while linear carbenium ions are predicted to be metastable at high temperature. Our findings confirm that carbenium ions, rather than alkoxides, are reactive intermediates in catalytic alkene cracking at 773 K.

Open Access version available at UGent repository

Water coordination and dehydration processes in defective UiO-66 type metal organic frameworks

M. Vandichel, J. Hajek, A. Ghysels, A. De Vos, M. Waroquier, V. Van Speybroeck
CrystEngComm
18 (37), 7056-7069
2016
A1

Abstract 

The UiO-66 metal organic framework is one of the most thermally and chemically stable hybrid materials reported to date. However, it is also accepted that the material contains structurally embedded defects, which may be engineered to enhance properties towards specific applications such as catalysis, sensing, etc. The synthesis conditions determine to a large extent the level of perfection of the material and additionally the catalytic activity may be enhanced by post-synthesis activation at high temperature under vacuum, in which defect coordinating species (H2O, HCl, monocarboxylic modulators, etc.) evaporate. The molecular level characterization of defects is extremely challenging from both theoretical and experimental points of view. Such experimental endeavor was recently proposed via experimental SXRD measurements, also unraveling the coordination of water on the Zr–O–Zr defect sites [Angew. Chem., Int. Ed., 2015, 54, 11162–11167]. The present work provides a theoretical understanding of defect structures in UiO-66(Zr) by means of periodic density functional theory calculations and ab initio molecular dynamics simulations. A range of defect structures are generated with different numbers of missing linkers. For each of the defects, the free energetic and mechanical stability is discussed and the coordination of water and charge balancing hydroxide ions is studied. For catalysis applications, the material is mostly pretreated to remove water by dehydration reactions. For each of the proposed defect structures, mechanistic pathways for dehydration reactions of the Zr-bricks are determined employing nudged elastic band (NEB) calculations. During the dehydroxylation trajectory, loose hydroxyl groups and terephthalate decoordinations are observed. Furthermore, dehydration reactions are lower activated if terephthalate linkers are missing in the immediate environment of the inorganic brick. The creation of defects and the dehydration processes have a large impact on the mechanical properties of the material, which is evidenced by lower bulk moduli and elastic constants for structures with more defects.

DOI 

10.1039/C6CE01027J

Pages

Subscribe to RSS - M. Waroquier