Dielectric response of metal-organic frameworks as a function of confined guest species investigated by molecular dynamics simulations.

When using metal-organic frameworks (MOFs) as electric field-dependent sensor devices, understanding their dielectric response is crucial as the orientation of polar groups is largely affected by confinement. To shed light on this at the molecular level, the response to a static field was computationally investigated for two structurally related MOFs, depending on their loading with guest molecules. The pillared-layer MOFs differ in their pillar moiety, with one bearing a rotatable permanent dipole moment and the other being non-polar. Two guest molecules with and without polarity, namely, methanol and methane, were considered. A comprehensive picture of the response of the guest molecules could be achieved with respect to both the amount and polarity of the confined species. For both MOFs, the dielectric response is very sensitive to the introduction of methanol, showing an anisotropic and non-linear increase in the system's relative permittivity expressed by a strongly increasing polarization response to external electric fields scaling with the number of confined methanol molecules. As expected, the effect of methane in the non-dipolar MOF is negligible, whereas subtle differences can be observed for the dipolar response of the MOF with rotatable dipolar linker groups. Taking advantage of these anisotropic and guest-molecule-specific confinement effects may open pathways for future sensing applications. Finally, methanol-induced global framework dynamics were observed in both MOFs.

Modulating Liquid-Liquid Transitions and Glass Formation in Zeolitic Imidazolate Frameworks by Decoration with Electron-Withdrawing Cyano Groups.

The liquid phase of metal-organic frameworks (MOFs) is key for the preparation of melt-quenched bulk glasses as well as the shaping of these materials for various applications; however, only very few MOFs can be melted and transformed into stable glasses. Here, the solvothermal and mechanochemical preparation of a new series of functionalized derivatives of ZIF-4 (Zn(im)2, where im- = imidazolate and ZIF = zeolitic imidazolate framework) containing the cyano-functionalized imidazolate linkers CNim- (4-cynanoimidazolate) and dCNim- (4,5-dicyanoimidazolate) is reported. The strongly electron-withdrawing nature of the CN groups facilitates low-temperature melting of the materials (below 310 °C for some derivatives) and the formation of microporous ZIF glasses with remarkably low glass-transition temperatures (down to only about 250 °C) and strong resistance against recrystallization. Besides conventional ZIF-4, the CN-functionalized ZIFs are so far the only MOFs to show an exothermic framework collapse to a low-density liquid phase and a subsequent transition to a high-density liquid phase. By systematic adjustment of the fraction of cyano-functionalized linkers in the ZIFs, we derive fundamental insights into the thermodynamics of the unique polyamorphic nature of these glass formers as well as further design rules for the porosity of the ZIF glasses and the viscosity of their corresponding liquids. The results provide new insights into the unusual phenomenon of liquid-liquid transitions as well as a guide for the chemical diversification of meltable MOFs, likely with implications beyond the archetypal ZIF glass formers.

Computing chemical potentials of adsorbed or confined fluids.

The chemical potential of adsorbed or confined fluids provides insight into their unique thermodynamic properties and determines adsorption isotherms. However, it is often difficult to compute this quantity from atomistic simulations using existing statistical mechanical methods. We introduce a computational framework that utilizes static structure factors, thermodynamic integration, and free energy perturbation for calculating the absolute chemical potential of fluids. For demonstration, we apply the method to compute the adsorption isotherms of carbon dioxide in a metal-organic framework and water in carbon nanotubes.

Understanding MOF Flexibility: An Analysis Focused on Pillared Layer MOFs as a Model System.

Flexible porous frameworks are at the forefront of materials research. A unique feature is their ability to open and close their pores in an adaptive manner induced by chemical and physical stimuli. Such enzyme-like selective recognition offers a wide range of functions ranging from gas storage and separation to sensing, actuation, mechanical energy storage and catalysis. However, the factors affecting switchability are poorly understood. In particular, the role of building blocks, as well as secondary factors (crystal size, defects, cooperativity) and the role of host-guest interactions, profit from systematic investigations of an idealized model by advanced analytical techniques and simulations. The review describes an integrated approach targeting the deliberate design of pillared layer metal-organic frameworks as idealized model materials for the analysis of critical factors affecting framework dynamics and summarizes the resulting progress in their understanding and application.

Photochemical Sandmeyer-type Halogenation of Arenediazonium Salts.

Trihalide salts were found to efficiently promote photochemical dediazotizing halogenations of diazonium salts. In contrast to classical Sandmeyer reactions, no metal catalysts are required to achieve high yields and outstanding selectivities for halogenation over competing hydridodediazotization. Convenient protocols are disclosed for synthetically meaningful brominations, iodinations, and chlorinations of diversely functionalized derivatives.

An automatized workflow from molecular dynamic simulation to quantum chemical methods to identify elementary reactions and compute reaction constants.

We present an automatized workflow which, starting from molecular dynamics simulations, identifies reaction events, filters them, and prepares them for accurate quantum chemical calculations using, for example, Density Functional Theory (DFT) or Coupled Cluster methods. The capabilities of the automatized workflow are demonstrated by the example of simulations for the combustion of some polycyclic aromatic hydrocarbons (PAHs). It is shown how key elementary reaction candidates are filtered out of a much larger set of redundant reactions and refined further. The molecular species in question are optimized using DFT and reaction energies, barrier heights, and reaction rates are calculated. The setup is general enough to include at this stage configurational sampling, which can be exploited in the future. Using the introduced machinery, we investigate how the observed reaction types depend on the gas atmosphere used in the molecular dynamics simulation. For the re-optimization on the DFT level, we show how the additional information needed to switch from reactive force-field to electronic structure calculations can be filled in and study how well ReaxFF and DFT agree with each other and shine light on the perspective of using more accurate semi-empirical methods in the MD simulation.

Influence of flexible side-chains on the breathing phase transition of pillared layer MOFs: a force field investigation.

The prototypical pillared layer MOFs, formed by a square lattice of paddle-wheel units and connected by dinitrogen pillars, can undergo a breathing phase transition by a "wine-rack" type motion of the square lattice. We studied this behavior, which is not yet fully understood, using an accurate first principles parameterized force field (MOF-FF) for larger nanocrystallites on the example of Zn2(bdc)2(dabco) [bdc: benzenedicarboxylate, dabco: (1,4-diazabicyclo[2.2.2]octane)], and found clear indications for an interface between a closed and an open pore phase traveling through the system during the phase transformation [J. Keupp and R. Schmid, Adv. Theory Simul., 2019, 2, 1900117]. In conventional simulations in small supercells this mechanism is prevented by periodic boundary conditions (PBCs), enforcing a synchronous transformation of the entire crystal. Here, we extend this investigation to pillared layer MOFs with flexible side-chains, attached to the linker. Such functionalized (fu-)MOFs are experimentally known to have different properties with the side-chains acting as fixed guest molecules. First, in order to extend the parameterization for such flexible groups, a new parameterization strategy for MOF-FF had to be developed, using a multi-structure force based fit method. The resulting parameterization for a library of fu-MOFs is then validated with respect to a set of reference systems and shows very good accuracy. In the second step, a series of fu-MOFs with increasing side-chain length is studied with respect to the influence of the side-chains on the breathing behavior. For small supercells in PBCs a systematic trend of the closed pore volume with the chain length is observed. However, for a nanocrystallite model a distinct interface between a closed and an open pore phase is visible only for the short chain length, whereas for longer chains the interface broadens and a nearly concerted transformation is observed. Only by molecular dynamics simulations using accurate force fields can such complex phenomena can be studied on a molecular level.

Configurational Entropy Driven High-Pressure Behaviour of a Flexible Metal-Organic Framework (MOF).

Flexible metal-organic frameworks (MOFs) show large structural flexibility as a function of temperature or (gas)pressure variation, a fascinating property of high technological and scientific relevance. The targeted design of flexible MOFs demands control over the macroscopic thermodynamics as determined by microscopic chemical interactions and remains an open challenge. Herein we apply high-pressure powder X-ray diffraction and molecular dynamics simulations to gain insight into the microscopic chemical factors that determine the high-pressure macroscopic thermodynamics of two flexible pillared-layer MOFs. For the first time we identify configurational entropy that originates from side-chain modifications of the linker as the key factor determining the thermodynamics in a flexible MOF. The study shows that configurational entropy is an important yet largely overlooked parameter, providing an intriguing perspective of how to chemically access the underlying free energy landscape in MOFs.

Experimental Evidence for the Incorporation of Two Metals at Equivalent Lattice Positions in Mixed-Metal Metal-Organic Frameworks.

Metal-organic frameworks containing multiple metals distributed over crystallographically equivalent framework positions (mixed-metal MOFs) represent an interesting class of materials, since the close vicinity of isolated metal centers often gives rise to synergistic effects. However, appropriate characterization techniques for detailed investigations of these mixed-metal metal-organic framework materials, particularly addressing the distribution of metals within the lattice, are rarely available. The synthesis of mixed-metal FeCuBTC materials in direct syntheses proved to be difficult and only a thorough characterization using various techniques, like powder X-ray diffraction, X-ray absorption spectroscopy and electron paramagnetic resonance spectroscopy, unambiguously evidenced the formation of a mixed-metal FeCuBTC material with HKUST-1 structure, which contained bimetallic Fe-Cu paddlewheels as well as monometallic Cu-Cu and Fe-Fe units under optimized synthesis conditions. The in-depth characterization showed that other synthetic procedures led to impurities, which contained the majority of the applied iron and were impossible or difficult to identify using solely standard characterization techniques. Therefore, this study shows the necessity to characterize mixed-metal MOFs extensively to unambiguously prove the incorporation of both metals at the desired positions. The controlled positioning of metal centers in mixed-metal metal-organic framework materials and the thorough characterization thereof is particularly important to derive structure-property or structure-activity correlations.

Can Small Polyaromatics Describe Their Larger Counterparts for Local Reactions? A Computational Study on the H-Abstraction Reaction by an H-Atom from Polyaromatics.

Hydrogen abstraction is one of the crucial initial key steps in the combustion of polycyclic aromatic hydrocarbons. For an accurate theoretical prediction of heterogeneous combustion processes, larger systems need to be treated as compared to pure gas phase reactions. We address here the question on how transferable activation and reaction energies computed for small molecular models are to larger polyaromatics. The approximate transferability of energy contributions is a key assumption for multiscale modeling approaches. To identify efficient levels of accuracy, we start with accurate coupled-cluster and density functional theory (DFT) calculations for different sizes of polyaromatics. More approximate methods as the reactive force-field ReaxFF and the extended semi-empirical tight binding (xTB) methods are then benchmarked against these data sets in terms of reaction energies and equilibrium geometries. Furthermore, we analyze the role of bond-breaking and relaxation energies, vibrational contributions, and post-Hartree-Fock correlation corrections on the reaction, and for the activation energies, we analyze the validity of the Bell-Evans-Polanyi and Hammond principles. First, we find good transferability for this process and that the predictivity of small models at high theoretical levels is way superior than any approximate method can deliver. Second, ReaxFF can serve as a qualitative exploration method, whereas GFN2-xTB in combination with GFN1-xTB appears as a favorable tool to bridge between DFT and ReaxFF so that we propose a multimethod scheme with employing ReaxFF, GFN1/GFN2-xTB, DFT, and coupled cluster to cope effectively with such a complex reactive system.

TopoFF: MOF structure prediction using specifically optimized blueprints.

For the structure prediction of MOFs and related crystalline framework materials we have proposed the Reversed Topological Approach (RTA), where the default embedding of a topology is used as a blueprint. The optimal rotational insertion of the building blocks (BBs) at the fixed vertex positions of the blueprint is performed by minimizing the target function of the average angle deviation (AAD). Here we extend this idea by pre-optimizing the maximum symmetry embedding of a topology in order to minimize the overall mean AAD for the given set of BBs. By this fast and essentially parameter-free topoFF method, the vertex positions and cell parameters of the blueprint are further optimized in order to fit the structural needs of the BBs, which speeds up the overall search for the most energetically favorable structure. In addition, different topologies can be ranked in a quantitative and intuitive way. The definition and implementation of topoFF is explained and its application for the RTA-based structure prediction of MOFs is demonstrated with a number of instructive examples.

Electrode potential dependent desolvation and resolvation of germanium(100) in contact with aqueous perchlorate electrolytes.

The electrode potential dependence of the hydration layer on an n-Ge(100) surface was studied by a combination of in situ and operando electrochemical attenuated total reflection infrared (ATR-IR) spectroscopy and real space density functional theory (DFT) calculations. Constant-potential DFT calculations were coupled to a modified generalised Poisson-Boltzmann ion distribution model and applied within an ab initio molecular dynamics (AIMD) scheme. As a result, potential-dependent vibrational spectra of surface species and surface water were obtained, both experimentally and by simulations. The experimental spectra show increasing absorbance from the Ge-H stretching modes at negative potentials, which is associated with an increased negative difference absorbance of water-related OH modes. When the termination transition of germanium from OH to H termination occurs, the surface switches from hydrophilic to hydrophobic. This transition is fully reversible. During the switching, the interface water molecules are displaced from the surface forming a "hydrophobic gap". The gap thickness was experimentally estimated by a continuum electrodynamic model to be ≈2 Å. The calculations showed a shift in the centre of mass of the interface water by ≈0.9 Å due to the surface transformation. The resulting IR spectra of the interfacial water in contact with the hydrophobic Ge-H show an increased absorbance of free OH groups, and a decreased absorbance of strongly hydrogen bound water. Consequently, the surface transformation to a Ge-H terminated surface leads to a surface which is weakening the H-bond network of the interfacial water in contact.

Verlet-like algorithms for Car-Parrinello molecular dynamics with unequal electronic occupations.

The ab initio molecular dynamics simulations of metallic, charged, and electrochemical systems require, in principle, the inclusion of unequally occupied electronic states. In this contribution, the general approach to work with fixed but arbitrary occupations within the Car-Parrinello molecular dynamics scheme is revisited, focusing on the procedure which is required to maintain the orthonormality constraints in the commonly used position-Verlet integrator. Expressions to constrain also the orbital velocities, as it is demanded by a velocity-Verlet integrator, are then derived. The generalized unequal-occupation SHAKE algorithm is compared with the standard procedure for damped dynamics (energy optimization) of systems including fully unoccupied electronic states. In turn, the proposed unequal-occupation RATTLE algorithm is validated by the corresponding microcanonical ensemble simulations. It is shown that only with the proper orthogonalization method, a correct ordering of states and energy conserving dynamics can be achieved.

QuickFF: A program for a quick and easy derivation of force fields for metal-organic frameworks from ab initio input.

QuickFF is a software package to derive accurate force fields for isolated and complex molecular systems in a quick and easy manner. Apart from its general applicability, the program has been designed to generate force fields for metal-organic frameworks in an automated fashion. The force field parameters for the covalent interaction are derived from ab initio data. The mathematical expression of the covalent energy is kept simple to ensure robustness and to avoid fitting deficiencies as much as possible. The user needs to produce an equilibrium structure and a Hessian matrix for one or more building units. Afterward, a force field is generated for the system using a three-step method implemented in QuickFF. The first two steps of the methodology are designed to minimize correlations among the force field parameters. In the last step, the parameters are refined by imposing the force field parameters to reproduce the ab initio Hessian matrix in Cartesian coordinate space as accurate as possible. The method is applied on a set of 1000 organic molecules to show the easiness of the software protocol. To illustrate its application to metal-organic frameworks (MOFs), QuickFF is used to determine force fields for MIL-53(Al) and MOF-5. For both materials, accurate force fields were already generated in literature but they requested a lot of manual interventions. QuickFF is a tool that can easily be used by anyone with a basic knowledge of performing ab initio calculations. As a result, accurate force fields are generated with minimal effort. © 2015 Wiley Periodicals, Inc.

Structural complexity in metal-organic frameworks: simultaneous modification of open metal sites and hierarchical porosity by systematic doping with defective linkers.

A series of defect-engineered metal-organic frameworks (DEMOFs) derived from parent microporous MOFs was obtained by systematic doping with defective linkers during synthesis, leading to the simultaneous and controllable modification of coordinatively unsaturated metal sites (CUS) and introduction of functionalized mesopores. These materials were investigated via temperature-dependent adsorption/desorption of CO monitored by FTIR spectroscopy under ultra-high-vacuum conditions. Accurate structural models for the generated point defects at CUS were deduced by matching experimental data with theoretical simulation. The results reveal multivariate diversity of electronic and steric properties at CUS, demonstrating the MOF defect structure modulation at two length scales in a single step to overcome restricted active site specificity and confined coordination space at CUS. Moreover, the DEMOFs exhibit promising modified physical properties, including band gap, magnetism, and porosity, with hierarchical micro/mesopore structures correlated with the nature and the degree of defective linker incorporation into the framework.

Tribenzotriquinacene receptors for C60 †fullerene rotors: towards C3 symmetrical chiral stators for unidirectionally operating nanoratchets.

The synthesis of a stereochemically pure concave tribenzotriquinacene receptor (7) for C60 fullerene, possessing C3 point group symmetry, by threefold condensation of C2 -symmetric 1,2-diketone synthons (5) and a hexaaminotribenzotriquinacene core (6) is described. The chiral diketone was synthesized in a five-step reaction sequence starting from C2h -symmetric 2,6-di-tert-butylanthracene. The highly diastereo-discriminating Diels-Alder reaction of 2,6-di-tert-butylanthracene with fumaric acid di(-)menthyl ester, catalyzed by aluminium chloride, is the relevant stereochemistry introducing step. The structure of the fullerene receptor was verified by (1)H and (13)C NMR spectroscopy, mass spectrometry and single crystal X-ray diffraction. VCD and ECD spectra were recorded, which were corroborated by ab initio DFT calculations, establishing the chiral nature of 7 with about 99.7 % ee, based on the ee (99.9 %) of the chiral synthon (1). The absolute configuration of 7 could thus be established as all-S [(2S,7S,16S,21S,30S,35S)-(7)]. Spectroscopic titration experiments reveal that the host forms 1:1 complexes with either pure fullerene (C60) or fullerene derivatives, such as rotor 1'-(4-nitrophenyl)-3'-(4-N,N-dimethylaminophenyl)-pyrazolino[4',5':1,2][60]fullerene (R). The complex stability constants of the complexes dissolved in CHCl3/CS2 (1:1 vol. %) are K([C60 ⊂7]) = 319(±156) M(-1) and K([R⊂7]) = 110(±50) M(-1). With molecular dynamics simulations using a first-principles parameterized force field the asymmetry of the rotational potential for [R⊂7] was shown, demonstrating the potential suitability of receptor 7 to act as a stator in a unidirectionally operating nanoratchet.

Almost enclosed buckyball joints: synthesis, complex formation, and computational simulations of pentypticene-extended tribenzotriquinacene.

We report the synthesis of a tribenzotriquinacene-based (TBTQ) receptor (3) for C60 fullerene, which is extended by pentiptycene moieties to provide an almost enclosed concave ball bearing. The system serves as a model for a self-assembling molecular rotor with a flexible and adapting stator. Unexpectedly, nuclear magnetic resonance spectroscopic investigations reveal a surprisingly low complex stability constant of K1 =213±37 M(-1) for [C60 ⊂3], seemingly inconsistent with the previously reported TBTQ systems. Molecular dynamics (MD) simulations have been conducted for three different [C60 ⊂TBTQ] complexes to resolve this. Because of the dominating dispersive interactions, the binding energies increase with the contact area between guest and host, however, only for rigid host structures. By means of free-energy calculations with an explicit solvent model it can be shown that the novel flexible TBTQ receptor 3 binds weakly because of hampering entropic contributions.