Experimental and Theoretical Evidence for the Promotional Effect of Acid Sites on the Diffusion of Alkenes through Small-Pore Zeolites.

The diffusion of saturated and unsaturated hydrocarbons is of fundamental importance for many zeolite-catalyzed processes. Transport of small alkenes in the confined zeolite pores can become hindered, resulting in a significant impact on the ultimate product selectivity and separation. Herein, intracrystalline light olefin/paraffin diffusion through the 8-ring windows of zeolite SAPO-34 is characterized by a complementary set of first-principle molecular dynamics simulations, PFG-NMR experiments, and pulse-response temporal analysis of products measurements, yielding information at different length and time scales.

Structural and Photophysical Properties of Various Polypyridyl Ligands: A Combined Experimental and Computational Study.

Covalent triazine frameworks (CTFs) with polypyridyl ligands are very promising supports to anchor photocatalytic complexes. Herein, we investigate the photophysical properties of a series of ligands which vary by the extent of the aromatic system, the nitrogen content and their topologies to aid in selecting interesting building blocks for CTFs. Interestingly, some linkers have a rotational degree of freedom, allowing both a trans and cis structure, where only the latter allows anchoring.

Unraveling the thermodynamic criteria for size-dependent spontaneous phase separation in soft porous crystals.

Soft porous crystals (SPCs) harbor a great potential as functional nanoporous materials owing to their stimuli-induced and tuneable morphing between different crystalline phases. These large-amplitude phase transitions are often assumed to occur cooperatively throughout the whole material, which thereby retains its perfect crystalline order. Here, we disprove this paradigm through mesoscale first-principles based molecular dynamics simulations, demonstrating that morphological transitions do induce spatial disorder under the form of interfacial defects and give rise to yet unidentified phase coexistence within a given sample.

A Supramolecular View on the Cooperative Role of Brønsted and Lewis Acid Sites in Zeolites for Methanol Conversion.

A systematic molecular level and spectroscopic investigation is presented to show the cooperative role of Brønsted acid and Lewis acid sites in zeolites for the conversion of methanol. Extra-framework alkaline-earth metal containing species and aluminum species decrease the number of Brønsted acid sites, as protonated metal clusters are formed. A combined experimental and theoretical effort shows that postsynthetically modified ZSM-5 zeolites, by incorporation of extra-framework alkaline-earth metals or by demetalation with dealuminating agents, contain both mononuclear [MOH] and double protonated binuclear metal clusters [M(μ-OH)M] (M = Mg, Ca, Sr, Ba, and HOAl).

Dynamic Interplay between Defective UiO-66 and Protic Solvents in Activated Processes.

UiO-66, composed by Zr-oxide inorganic bricks [Zr (μ -O) (μ -OH) ] and organic terephthalate linkers, is one of the most studied metal-organic frameworks (MOFs) due to its exceptional thermal, chemical, and mechanical stability. Thanks to its high connectivity, the material can withstand structural deformations during activation processes such as linker exchange, dehydration, and defect formation. These processes do alter the zirconium coordination number in a dynamic way, creating open metal sites for catalysis and thus are able to tune the catalytic properties.

On the importance of anharmonicities and nuclear quantum effects in modelling the structural properties and thermal expansion of MOF-5.

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.

Ab Initio Evaluation of Henry Coefficients Using Importance Sampling.

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.

Protocol for Identifying Accurate Collective Variables in Enhanced Molecular Dynamics Simulations for the Description of Structural Transformations in Flexible Metal-Organic Frameworks.

Various kinds of flexibility have been observed in metal-organic frameworks, which may originate from the topology of the material or the presence of flexible ligands. The construction of free energy profiles describing the full dynamical behavior along the phase transition path is challenging since it is not trivial to identify collective variables able to identify all metastable states along the reaction path. In this work, a systematic three-step protocol to uniquely identify the dominant order parameters for structural transformations in flexible metal-organic frameworks and subsequently construct accurate free energy profiles is presented.

How Chain Length and Branching Influence the Alkene Cracking Reactivity on H-ZSM-5.

Catalytic alkene cracking on H-ZSM-5 involves a complex reaction network with many possible reaction routes and often elusive intermediates. Herein, advanced molecular dynamics simulations at 773 K, a typical cracking temperature, are performed to clarify the nature of the intermediates and to elucidate dominant cracking pathways at operating conditions. A series of C-C alkene intermediates are investigated to evaluate the influence of chain length and degree of branching on their stability.

The Importance of Cell Shape Sampling To Accurately Predict Flexibility in Metal-Organic Frameworks.

In this work, the influence of cell shape sampling on the predicted stability of the different metastable phases in flexible metal-organic frameworks at finite temperatures is investigated. The influence on the free energy by neglecting cell shape sampling is quantified for the prototypical MIL-53(Al) and the topical DUT-49(Cu). This goal is achieved by constructing free energy profiles in ensembles either in which the phase space associated with the cell shape is sampled explicitly or in which the cell shape is kept fixed.

Extension of the QuickFF force field protocol for an improved accuracy of structural, vibrational, mechanical and thermal properties of metal-organic frameworks.

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 transparent workflow.

Thermodynamic insight into stimuli-responsive behaviour of soft porous crystals.

Knowledge of the thermodynamic potential in terms of the independent variables allows to characterize the macroscopic state of the system. However, in practice, it is difficult to access this potential experimentally due to irreversible transitions that occur between equilibrium states. A showcase example of sudden transitions between (meta)stable equilibrium states is observed for soft porous crystals possessing a network with long-range structural order, which can transform between various states upon external stimuli such as pressure, temperature and guest adsorption.

Influence of a Confined Methanol Solvent on the Reactivity of Active Sites in UiO-66.

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.

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

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.

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

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.

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

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.

The Remarkable Amphoteric Nature of Defective UiO-66 in Catalytic Reactions.

One of the major requirements in solid acid and base 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 coordinatively unsaturated zirconium sites, which may be tuned by removing framework linkers and by removal of water from the inorganic bricks using a dehydration treatment.

The Monomer Electron Density Force Field (MEDFF): A Physically Inspired Model for Noncovalent Interactions.

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.

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

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.

Exploring the Flexibility of MIL-47(V)-Type Materials Using Force Field Molecular Dynamics Simulations.

The flexibility of three MIL-47(V)-type materials (MIL-47, COMOC-2, and COMOC-3) has been explored by constructing the pressure versus volume and free energy versus volume profiles at various temperatures ranging from 100 to 400 K. This is done with first-principles-based force fields using the recently proposed QuickFF parametrization protocol. Specific terms were added for the materials at hand to describe the asymmetry of the one-dimensional vanadium-oxide chain and to account for the flexibility of the organic linkers.

Minimal Basis Iterative Stockholder: Atoms in Molecules for Force-Field Development.

Atomic partial charges appear in the Coulomb term of many force-field models and can be derived from electronic structure calculations with a myriad of atoms-in-molecules (AIM) methods. More advanced models have also been proposed, using the distributed nature of the electron cloud and atomic multipoles. In this work, an electrostatic force field is defined through a concise approximation of the electron density, for which the Coulomb interaction is trivially evaluated.

A Comparison of Barostats for the Mechanical Characterization of Metal-Organic Frameworks.

In this paper, three barostat coupling schemes for pressure control, which are commonly used in molecular dynamics simulations, are critically compared to characterize the rigid MOF-5 and flexible MIL-53(Al) metal-organic frameworks. We investigate the performance of the three barostats, the Berendsen, the Martyna-Tuckerman-Tobias-Klein (MTTK), and the Langevin coupling methods, in reproducing the cell parameters and the pressure versus volume behavior in isothermal-isobaric simulations. A thermodynamic integration method is used to construct the free energy profiles as a function of volume at finite temperature.

PPV Polymerization through the Gilch Route: Diradical Character of Monomers.

Despite various studies on the polymerization of poly(p-phenylene vinylene) (PPV) through different precursor routes, detailed mechanistic knowledge on the individual reaction steps and intermediates is still incomplete. The present study aims to gain more insight into the radical polymerization of PPV through the Gilch route. The initial steps of the polymerization involve the formation of a p-quinodimethane intermediate, which spontaneously self-initiates through a dimerization process leading to the formation of diradical species; chain propagation ensues on both sides of the diradical or chain termination occurs by the formation of side products, such as [2.