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.

Simulating the structural phase transitions of metal-organic frameworks with control over the volume of nanocrystallites.

Flexible metal-organic frameworks (MOFs) can undergo structural transitions with significant pore volume changes upon guest adsorption or other external triggers while maintaining their porosity. In computational studies of this breathing behavior, molecular dynamics (MD) simulations within periodic boundary conditions (PBCs) are commonly performed. However, to account for the finite size and surface effects affecting the phase transition mechanism, the simulation of non-periodic nanocrystallite (NC) models without the constraint of PBCs is an important alternative. In this study, we present an approach allowing the analysis and control of the volume of finite-size structures during MD simulations by a tetrahedral tessellation of the (deformed) NC's volume. The method allows for defining the current NC's volume during the simulation and manipulating it regarding a particular reference volume to compute free energies for the phase transformation via umbrella sampling. The application on differently sized DMOF-1 and DUT-128 NCs reveals flexible pore closing mechanisms without significant biasing of the transition pathway. The concept provides the theoretical foundation for further research on flexible materials regarding targeted initialization of the structural phase behavior to elucidate the underlying mechanism, which can be used to improve the applications of flexible materials by targeted controlling of the phase transition.

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.

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.

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.

Exploring the Impact of the Linker Length on Heat Transport in Metal-Organic Frameworks.

Metal-organic frameworks (MOFs) are a highly versatile group of porous materials suitable for a broad range of applications, which often crucially depend on the MOFs' heat transport properties. Nevertheless, detailed relationships between the chemical structure of MOFs and their thermal conductivities are still largely missing. To lay the foundations for developing such relationships, we performed non-equilibrium molecular dynamics simulations to analyze heat transport in a selected set of materials. In particular, we focus on the impact of organic linkers, the inorganic nodes and the interfaces between them. To obtain reliable data, great care was taken to generate and thoroughly benchmark system-specific force fields building on ab-initio-based reference data. To systematically separate the different factors arising from the complex structures of MOF, we also studied a series of suitably designed model systems. Notably, besides the expected trend that longer linkers lead to a reduction in thermal conductivity due to an increase in porosity, they also cause an increase in the interface resistance between the different building blocks of the MOFs. This is relevant insofar as the interface resistance dominates the total thermal resistance of the MOF. Employing suitably designed model systems, it can be shown that this dominance of the interface resistance is not the consequence of the specific, potentially weak, chemical interactions between nodes and linkers. Rather, it is inherent to the framework structures of the MOFs. These findings improve our understanding of heat transport in MOFs and will help in tailoring the thermal conductivities of MOFs for specific applications.

Thermicity of the Decomposition of Oxygen Functional Groups on Cellulose-Derived Chars.

The evolution of oxygen functional groups (OFGs) and the associated thermic effects upon heat treatment up to 800 °C were investigated experimentally as well as by theoretical calculations. A synthetic carbon with a carbonaceous structure close to that of natural chars, yet mineral-free, was derived from cellulose and oxidized by HNO3 vapor at different temperatures and for varied durations in order to generate char samples with different concentrations and distributions of OFGs. The functionalized samples were subjected to calorimetric temperature-programmed desorption measurements in correlation with an extensive effluent gas analysis, thereby focusing on the specific heat effects of individual OFG evolution. Interpretation of the experimental results was aided by density functional theory (DFT) calculations which allowed one to infer the thermal stability of different OFGs and the reaction energy associated with their evolution upon heating. Results showed that, with increasing temperature, H2O was released due to the loss of physisorbed water, the decomposition of clusters bound to carboxylic acids, and condensation reactions. The associated heat uptake amounted to about 100 kJ mol-1. Contrarily, the release of CO2, attributed to the decomposition of condensed acids, carboxylic acids, anhydrides, and lactones, resulted in a heat release of about 40 kJ mol-1. The most strongly pronounced thermic effects were detected for the release of CO, comprising highly exothermic effects due to the decomposition of condensed acids and carbonyls/quinones as well as endothermic effects attributed to anhydrides and phenols/ethers. Notably, anhydrides can be formed during the oxidative treatment as well as during heating by condensation of adjacent carboxylic acids. In the latter case, the two-step decomposition is overall highly exothermic, indicating the associated occurrence of pronounced carbon matrix rearrangements. DFT investigations suggest that these rearrangements not only affect the immediate OFG proximity but also involve several carbon sheets.

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.

Molecular Dynamics Simulations of the Breathing Phase Transition of MOF Nanocrystallites II: Explicitly Modeling the Pressure Medium.

One of the most investigated properties of porous crystalline metal-organic frameworks (MOFs) is their potential flexibility to undergo large changes in unit cell size upon guest adsorption or other stimuli, referred to as "breathing". Computationally, such phase transitions are usually investigated using periodic boundary conditions, where the system's volume can be controlled directly. However, we have recently shown that important aspects like the formation of a moving interface between the open and the closed pore form or the free energy barrier of the first-order phase transition and its size effects can best be investigated using non-periodic nanocrystallite (NC) models [Keupp et al. (Adv. Theory Simul., 2019, 2, 1900117)]. In this case, the application of pressure is not straightforward, and a distance constraint was used to mimic a mechanical strain enforcing the reaction coordinate. In contrast to this prior work, a mediating particle bath is used here to exert an isotropic hydrostatic pressure on the MOF nanocrystallites. The approach is inspired by the mercury nanoporosimetry used to compress flexible MOF powders. For such a mediating medium, parameters are presented that require a reasonable additional numerical effort and avoid unwanted diffusion of bath particles into the MOF pores. As a proof-of-concept, NCs of pillared-layer MOFs with different linkers and sizes are studied concerning their response to external pressure exerted by the bath. By this approach, an isotropic pressure on the NC can be applied in analogy to corresponding periodic simulations, without any bias for a specific mechanism. This allows a more realistic investigation of the breathing phase transformation of a MOF NC and further bridges the gap between experiment and simulation.

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.

Frustrated flexibility in metal-organic frameworks.

Stimuli-responsive flexible metal-organic frameworks (MOFs) remain at the forefront of porous materials research due to their enormous potential for various technological applications. Here, we introduce the concept of frustrated flexibility in MOFs, which arises from an incompatibility of intra-framework dispersion forces with the geometrical constraints of the inorganic building units. Controlled by appropriate linker functionalization with dispersion energy donating alkoxy groups, this approach results in a series of MOFs exhibiting a new type of guest- and temperature-responsive structural flexibility characterized by reversible loss and recovery of crystalline order under full retention of framework connectivity and topology. The stimuli-dependent phase change of the frustrated MOFs involves non-correlated deformations of their inorganic building unit, as probed by a combination of global and local structure techniques together with computer simulations. Frustrated flexibility may be a common phenomenon in MOF structures, which are commonly regarded as rigid, and thus may be of crucial importance for the performance of these materials in various applications.

Beyond the Scholl reaction - one-step planarization and edge chlorination of nanographenes by mechanochemistry.

The edge chlorination of the benchmark nanographenes triphenylene and hexa-peri-hexabenzocoronene is conducted mechanochemically. This approach overcomes solubility limitations and eliminates the need for elaborate chlorination conditions. Additionally, the planarization of oligophenylenes and their edge-chlorination can be combined in a one-pot approach requiring as little as 60 minutes.

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.

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.

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.

Evaluating Computational Shortcuts in Supercell-Based Phonon Calculations of Molecular Crystals: The Instructive Case of Naphthalene.

Phonons crucially impact a variety of properties of organic semiconductor materials. For instance, charge- and heat transport depend on low-frequency phonons, while for other properties, such as the free energy, especially high-frequency phonons count. For all these quantities one needs to know the entire phonon band structure, whose simulation becomes exceedingly expensive for more complex systems when using methods like dispersion-corrected density functional theory (DFT). Therefore, in the present contribution we evaluate the performance of more approximate methodologies, including density functional tight binding (DFTB) and a pool of force fields (FF) of varying complexity and sophistication. Beyond merely comparing phonon band structures, we also critically evaluate to what extent derived quantities, like temperature-dependent heat capacities, mean squared thermal displacements, and temperature-dependent free energies are impacted by shortcomings in the description of the phonon bands. As a benchmark system, we choose (deuterated) naphthalene, as the only organic semiconductor material for which to date experimental phonon band structures are available in the literature. Overall, the best performance among the approximate methodologies is observed for a system-specifically parametrized second-generation force field. Interestingly, in the low-frequency regime also force fields with a rather simplistic model for the bonding interactions (like the General Amber Force Field) perform rather well. As far as the tested DFTB parametrization is concerned, we obtain a significant underestimation of the unit-cell volume resulting in a pronounced overestimation of the phonon energies in the low-frequency region. This cannot be mended by relying on the DFT-calculated unit cell, since with this unit cell the DFTB phonon frequencies significantly underestimate the experiments.

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.